Whole transcriptome profiling of the vernalization process in Lilium longiflorum (cultivar White Heaven) bulbs

Functional analysis of a lily SHORT VEGETATIVE PHASE ortholog in flowering transition and floral development

Lilium is a commercially important genus of bulbous flowers, investigating the flowering molecular mechanisms is important for flowering regulation of lily. MADS-box SHORT VEGETATIVE PHASE (SVP) orthologs are involved in the flowering transition and floral organ differentiation in many plants. In this study, we identified an SVP ortholog from L. × formolongi (LfSVP), which was closely related to Arabidopsis SVP according to phylogenetic analysis. Tissue-specific expression patterns indicated that LfSVP expression levels peaked in the leaves and showed low expression levels in flowering tepals. Stage-dependent expression patterns of LfSVP showed high transcription level in the flowering induction stage under different photoperiods and exhibited transcription peak in the floral budding development stage under long days. Overexpressed LfSVP led to delayed flowering and floral organ defects in Arabidopsis independent of photoperiod. Tobacco rattle virus -induced gene silencing of LfSVP caused a strongly earlier flowering time and floral organ defects of L. × formolongi. Moreover, LfSVP can interact with L. × formolongi APETALA1 (AP1) in both yeast and tobacco cells, and the two may interact to regulate floral organ differentiation. In conclusion, LfSVP is a flowering repressor and may be involved in the regulation of floral organ differentiation. This study will be helpful for the molecular breeding of short-life-period and rich floral patterns lily varieties.

Designing of future ornamental crops: a biotechnological driven perspective

With a basis in human appreciation of beauty and aesthetic values, the new era of ornamental crops is based on implementing innovative technologies and transforming symbols into tangible assets. Recent advances in plant biotechnology have attracted considerable scientific and industrial interest, particularly in terms of modifying desired plant traits and developing future ornamental crops. By utilizing omics approaches, genomic data, genetic engineering, and gene editing tools, scientists have successively explored the underlying molecular mechanism and potential gene(s) behind trait regulation such as floral induction, plant architecture, stress resistance, plasticity, adaptation, and phytoremediation in ornamental crop species. These signs of progress lay a theoretical and practical foundation for designing and enhancing the efficiency of ornamental plants for a wide range of applications. In this review, we briefly summarized the existing literature and advances in biotechnological approaches for the improvement of vital traits in ornamental plants. The future ornamental plants, such as light-emitting plants, biotic/abiotic stress detectors, and pollution abatement, and the introduction of new ornamental varieties via domestication of wild species are also discussed.

Transcriptome analysis of critical genes related to flowering in Mikania micrantha at different altitudes provides insights for a potential control

BACKGROUND

Mikania micrantha is a vine with strong invasion ability, and its strong sexual reproduction ability is not only the main factor of harm, but also a serious obstacle to control. M. micrantha spreads mainly through seed production. Therefore, inhibiting the flowering and seed production of M. micrantha is an effective strategy to prevent from continuing to spread.

RESULT

The flowering number of M. micrantha is different at different altitudes. A total of 67.01 Gb of clean data were obtained from nine cDNA libraries, and more than 83.47% of the clean reads were mapped to the reference genome. In total, 5878 and 7686 significantly differentially expressed genes (DEGs) were found in E2 vs. E9 and E13 vs. E9, respectively. Based on the background annotation and gene expression, some candidate genes related to the flowering pathway were initially screened, and their expression levels in the three different altitudes in flower bud differentiation showed the same trend. That is, at an altitude of 1300 m, the flower integration gene and flower meristem gene were downregulated (such as SOC1 and AP1), and the flowering inhibition gene was upregulated (such as FRI and SVP). Additionally, the results showed that there were many DEGs involved in the hormone signal transduction pathway in the flower bud differentiation of M. micrantha at different altitudes.

CONCLUSIONS

Our results provide abundant sequence resources for clarifying the underlying mechanisms of flower bud differentiation and mining the key factors inhibiting the flowering and seed production of M. micrantha to provide technical support for the discovery of an efficient control method.

Expression Patterns of Key Genes in the Photoperiod and Vernalization Flowering Pathways in Lilium longiflorum with Different Bulb Sizes

Lilium longiflorum is a wild Lilium, and its flowering transition requires a long period of cold exposure to meet the demand of vernalization. The responses of different sized bulbs to cold exposure and photoperiod are different, and the floral transition pathways of small and large bulbs are different. In this study, small and large bulbs were placed in cold storage for different weeks and then cultured at a constant ambient temperature of 25 °C under long day (LD) and short day (SD) conditions. Then, the flowering characteristics and expression patterns of key genes related to the vernalization and photoperiod pathways in different groups were calculated and analyzed. The results showed that the floral transition of Lilium longiflorum was influenced by both vernalization and photoperiod, that vernalization and LD conditions can significantly improve the flowering rate of Lilium longiflorum, and that the time from planting to visible flowering buds’ appearance was decreased. The flowering time and rate of large bulbs were greatly influenced by cold exposure, and the vernalization pathway acted more actively at the floral transition stage. The floral transition of small bulbs was affected more by the photoperiod pathway. Moreover, it was speculated that cold exposure may promote greater sensitivity of the small bulbs to LD conditions. In addition, the expression of LlVRN1, LlFKF1, LlGI, LlCO5, LlCO7, LlCO16, LlFT1, LlFT3 and LlSOC1 was high during the process of floral transition, and LlCO13, LlCO14 and LlCO15 were highly expressed in the vegetative stage. The expression of LlCO13 and LlCO14 was different under different lighting conditions, and the flowering induction function of LlCO9 and LlFT3 was related to vernalization. Moreover, LlFKF1, LlGI, LlCO5, LlCO16, LlSOC1 and LlFT2 were involved in the entire growth process of plants, while LlCO6, LlCO16 and LlFT1 are involved in the differentiation and formation of small bulblets of plants after the inflorescence stage, and this process is also closely related to LD conditions. This study has great significance for understanding the molecular mechanisms of the vernalization and photoperiod flowering pathways of Lilium longiflorum.

Flowering Mechanisms and Environmental Stimuli for Flower Transition: Bases for Production Scheduling in Greenhouse Floriculture

The scheduling of plant production is a critical aspect in modern floriculture since nowadays, sales are not oriented toward the recurring holidays as in the past, but always more toward impulse buying, implying a more diverse and constant demand on the market. This requires continuous production, often regulated by precise commercial agreements between growers and buyers, and between buyers and dealers, particularly in large-scale retail trade. In this scenario, diverse techniques to modulate the duration of the growing cycle, by hastening or slowing down plant growth and development, have been developed to match plant flowering to the market demand. Among the numerous approaches, the manipulation of climatic parameters in the growth environment is one of the most common in greenhouse floriculture. In this review, we summarize the physiological and biochemical bases underlying the main mechanisms of flowering, depending on the plant reaction to endogenous signals or environmental stimuli. In addition, the strategies based on the control of temperature (before or after planting) and light environment (as light intensity and spectrum, and the photoperiod) in the scheduling of flower and ornamental crop production are briefly described.

When vegetation indicates reproduction: The affinity between leaf morphology and flowering commitment in the lily meristem

At the reproductive stage, lily plants bear two morphological types of mature leaves, one at the lower and one at the upper part of the stem. At the vegetative stage, all the leaves are similar to each other and to the reproductive plant's lower leaves. This heterophylly has not yet been explored. In this study, we show that it is not a result of the plant's age but rather an outcome of floral induction. The induction appears as an on-going process, during which the meristem still produces leaves but progressively becomes committed to reproduction. This intermediate period lasts until the ultimate switch to flower primordia occurs. The leaves produced during floral induction, termed here as "inductive," appear at the upper part of the stem. Besides their typical higher stomata density, these leaves have a poly-layered palisade mesophyll, whose cells exhibit a unique morphology and contain more chlorophyll than leaves of vegetative plants. These leaves display higher carbon assimilation, soluble sugar production, and chloroplast-lipid accumulation. Accordingly, genes associated with stomata, chloroplast, and photosynthesis are upregulated in these leaves. Our results were obtained when floral induction was achieved either by vernalization or photoperiod signals, ruling out a mere environmental effect. We suggest that lily plants prepare themselves for the high-energy-demanding bloom by producing leaves with enhanced photosynthetic capacity, leading to an increase in soluble sugars. These novel findings introduce an adjacent affinity between photosynthesis and flowering and provide a nondestructive tool for identifying the plant's developmental stage-vegetative or reproductive.

SMRT and Illumina RNA Sequencing and Characterization of a Key NAC Gene LoNAC29 during the Flower Senescence in Lilium oriental 'Siberia'

Lily (Lilium spp.) is an important cut flower around the world. Flower senescence in lilies is characterized by the wilting and abscission of tepals, which results in a decrease in flower quality and huge economic loss. However, the mechanism underlying flower senescence in lilies is largely unknown. In this study, single-molecule, real-time (SMRT) and Illumina sequencing were carried out in L. oriental ‘Siberia’. Sequencing yielded 73,218 non-redundant transcripts, with an N50 of 3792 bp. These data were further integrated with three published transcriptomes through cogent analysis, which yielded 62,960 transcripts, with an increase in N50 of 3935 bp. Analysis of differentially expressed genes showed that 319 transcription factors were highly upregulated during flower senescence. The expression of twelve NAC genes and eleven senescence-associated genes (SAGs) showed that LoNAC29 and LoSAG39 were highly expressed in senescent flowers. Transient overexpression of LoNAC29 and LoSAG39 in tepals of lily notably accelerated flower senescence, and the promoter activity of LoSAG39 was strongly induced by LoNAC29. This work supported new evidence for the molecular mechanism of flower senescence and provided better sequence data for further study in lilies.

Genome-Wide Transcriptomic Identification and Functional Insight of Lily WRKY Genes Responding to Botrytis Fungal Disease

Genome-wide transcriptome analysis using RNA-Seq of Lilium longiflorum revealed valuable genes responding to biotic stresses. WRKY transcription factors are regulatory proteins playing essential roles in defense processes under environmental stresses, causing considerable losses in flower quality and production. Thirty-eight WRKY genes were identified from the transcriptomic profile from lily genotypes, exhibiting leaf blight caused by Botrytis elliptica. Lily WRKYs have a highly conserved motif, WRKYGQK, with a common variant, WRKYGKK. Phylogeny of LlWRKYs with homologous genes from other representative plant species classified them into three groups- I, II, and III consisting of seven, 22, and nine genes, respectively. Base on functional annotation, 22 LlWRKY genes were associated with biotic stress, nine with abiotic stress, and seven with others. Sixteen unique LlWRKY were studied to investigate responses to stress conditions using gene expression under biotic and abiotic stress treatments. Five genes-LlWRKY3, LlWRKY4, LlWRKY5, LlWRKY10, and LlWRKY12-were substantially upregulated, proving to be biotic stress-responsive genes in vivo and in vitro conditions. Moreover, the expression patterns of LlWRKY genes varied in response to drought, heat, cold, and different developmental stages or tissues. Overall, our study provides structural and molecular insights into LlWRKY genes for use in the genetic engineering in Lilium against Botrytis disease.

Physiological and transcriptomic responses of Lanzhou Lily (Lilium davidii, var. unicolor) to cold stress

Low temperature induces changes in plants at physiological and molecular levels, thus affecting growth and development. The Lanzhou lily (Lilium davidii, var. unicolor) is an important medicinal plant with high economic value. However, the molecular mechanisms underlying its photosynthetic and antioxidation responses to low temperature still remain poorly understood. This study subjected the Lanzhou lily to the two temperatures of 20°C (control) and 4°C (low temperature) for 24 h. Physiological parameters related to membrane integrity, photosynthesis, antioxidant system, and differentially expressed genes were investigated. Compared with control, low temperature increased the relative electrical conductivity by 43.2%, while it decreased net photosynthesis rate, ratio of variable to maximal fluorescence, and catalase activity by 47.3%, 10.1%, and 11.1%, respectively. In addition, low temperature significantly increased the content of soluble protein, soluble sugar, and proline, as well as the activity of superoxide dismutase and peroxidase. Comparative transcriptome profiling showed that a total of 238,109 differentially expressed genes were detected. Among these, 3,566 were significantly upregulated while 2,982 were significantly downregulated in response to low temperature. Gene Ontology enrichment analysis indicated that in response to low temperature, the mostly significantly enriched differentially expressed genes were mainly involved in phosphorylation, membrane and protein kinase activity, as well as photosynthesis, light harvesting, light reaction, and alpha,alpha-trehalose-phosphate synthase activity. Kyoto Encyclopedia of Genes and Genomes enrichment analysis also indicated that the most significantly enriched pathways involved ribosome biogenesis in eukaryotes, phenylalanine metabolism, circadian rhythm, porphyrin and chlorophyll metabolism, photosynthesis of antenna proteins, photosynthesis, and carbon fixation in photosynthetic organisms. Moreover, the expression patterns of 10 randomly selected differentially expressed genes confirmed the RNA-Seq results. These results expand the understanding of the physiological and molecular mechanisms underlying the response of the Lanzhou lily to low temperature stress.

Expression of LhFT1, the Flowering Inducer of Asiatic Hybrid Lily, in the Bulb Scales

Asiatic hybrid lily leaves emerge from their bulbs in spring, after cold exposure in winter, and the plant then blooms in early summer. We identified four FLOWERING LOCUS T (FT)-like genes, LhFT1, LhFT4, LhFT6, and LhFT8, from an Asiatic hybrid lily. Floral bud differentiation initiated within bulbs before the emergence of leaves. LhFT genes were mainly expressed in bulb scales, and hardly in leaves, in which the FT-like genes of many plants are expressed in response to environmental signals. LhFT1 was expressed in bulb scales after vernalization and was correlated to flower bud initiation in two cultivars with different flowering behaviors. LhFT8 was upregulated in bulb scales after cold exposure and three alternative splicing variants with a nonsense codon were simultaneously expressed. LhFT6 was upregulated in bulb scales after flower initiation, whereas LhFT4 was expressed constantly in all organs. LhFT1 overexpression complemented the late-flowering phenotype of Arabidopsis ft-10, whereas that of LhFT8 did so partly. LhFT4 and LhFT6 overexpression could not complement. Yeast two-hybrid and in vitro analyses showed that the LhFT1 protein interacted with the LhFD protein. LhFT6 and LhFT8 proteins also interacted with LhFD, as observed in AlphaScreen assay. Based on these results, we revealed that LhFT1 acts as a floral activator during floral bud initiation in Asiatic hybrid lilies. However, the biological functions of LhFT4, LhFT6, and LhFT8 remain unclear.

Virus-Induced Gene Silencing in Lilies Using Cucumber Mosaic Virus Vectors

Virus-induced gene silencing (VIGS) systems are effective for rapid analysis of gene functions in plants that require a long period of growth such as Lilium. We successfully developed a VIGS system using the cucumber mosaic virus (HL strain, CMV-HL) vector to induce RNA silencing of the L. leichtlinii phytoene desaturase gene (LlPDS), where at 30 days postinoculation (dpi), photo-bleaching was observed in the upper leaves of L. leichtlinii, and at 57 dpi, white regions appeared on flower tepals that accumulate orange carotenoids. This vector spreads in bulbs, and it could induce silencing on emerged shoots in the following year. The CMV-HL vector can be easily constructed by insertion of a 30-60 nt fragment into the cloning site of the RNA3 genome. In this chapter, we describe how to use the CMV-HL vector system in the context of Lilium plants.

De novo transcriptome sequencing and gene expression profiling with/without B-chromosome plants of Lilium amabile

Supernumerary B chromosomes were found in Lilium amabile (2n = 2x = 24), an endemic Korean lily that grows in the wild throughout the Korean Peninsula. The extra B chromosomes do not affect the host-plant morphology; therefore, whole transcriptome analysis was performed in 0B and 1B plants to identify differentially expressed genes. A total of 154,810 transcripts were obtained from over 10 Gbp data by de novo assembly. By mapping the raw reads to the de novo transcripts, we identified 7,852 differentially expressed genes (log₂FC > |10|), in which 4,059 and 3,794 were up-and down-regulated, respectively, in 1B plants compared to 0B plants. Functional enrichment analysis revealed that various differentially expressed genes were involved in cellular processes including the cell cycle, chromosome breakage and repair, and microtubule formation; all of which may be related to the occurrence and maintenance of B chromosomes. Our data provide insight into transcriptomic changes and evolution of plant B chromosomes and deliver an informative database for future study of B chromosome transcriptomes in the Korean lily.

Transcriptome analysis reveals the different compatibility between LAAA × AA and LAAA × LL in Lilium

To unveil the mechanism of the compatibility of odd-allotetraploid lily (LAAA) as female with diploid male lily, the differences of expressed unigenes in the ovaries and leaves between LAAA × AA and LAAA × LL were investigated using transcriptome analysis. The results showed the fruits of LAAA × AA well developed, while those of LAAA × LL aborted. The number of differentially expressed genes was less in the ovaries of LAAA × AA than those of LAAA × LL, but it showed opposite trend in those of leaves. The unigenes related with auxins, cytokinins, gibberellins, antioxidants, expansins, chlorophylls, carbohydrates, transport proteins were usually up-expressed in the ovaries and leaves of LAAA × AA but not in LAAA × LL; while those of abscisic acid, ethylene, jasmonic acid, and salicylic acid were increased in the ovaries or leaves of LAAA × LL but not in LAAA × AA. The up-expressed unigenes in the ovaries and leaves of LAAA × AA played positive roles in its fruit development because the products of the genes, like phytohormones and antioxidants, had functions protecting leaves from senescence or scavenging ROS, and thus LAAA was compatible with AA, while those of LAAA × LL played negative roles and caused its fruits aborted, and hence LAAA was incompatible with LL.

Comparative Transcriptome Analysis of Temperature-Induced Green Discoloration in Garlic

Green discoloration is one of the most important problems that cause low quality of product in the processing of garlic, which can be induced by low-temperature stress. But the mechanism of low temperature-induced green discoloration is poorly understood. In the present study, the control garlic and three low temperature-treated garlic samples (stored at 4°C with 10, 15, and 40 days, respectively) were used for genome-wide transcriptome profiling analysis. A total of 49280 garlic unigenes with an average length of 1337 bp were de novo assembled, 20231 of which were achieved for functional annotation. When being suffered from 10, 15, and 40 days of low-temperature treatment, an increased degree of discoloration was observed, and a total of 4757, 4401, and 2034 unigenes showed a differential expression, respectively. Finally, 5923 differentially expressed genes (DEGs) were found to respond to the low-temperature stress, of which 3921 were identified in at least two treatments. Among these stress-responsive unigenes, there were large numbers of enzyme-encoding genes, which significantly enriched the pathway “proteasome,” many genes of which are potentially involved in the garlic discoloration, such as 7 alliinase-encoding genes, 5 γ-glutamyltranspeptidase-encoding genes, and 1 δ-aminolevulinic acid dehydratase-encoding gene. These stress-responsive enzyme-encoding genes are possibly responsible for the low-temperature-induced garlic discoloration. The identification of large numbers of DEGs provides a basis for further elucidating the mechanism of low-temperature-induced green discoloration in garlic.

The proof is in the bulb: glycerol influences key stages of lily development

A bulb is a whole plant condensed into an underground organ. A geophyte's bulb comprises both food reserves and important developmental history that may affect its whole growth. In Easter lily (Lilium longiflorum), bulb size is associated with the plant's flowering pathway - vernalization or photoperiod - and also affects sprouting, flower quality and abortion rate. The aim of this study was to investigate the reasons for the major physiological differences between large and small bulbs. Lily bulbs start their development from secondary meristems along the stem, with large bulbs being heavier and bear more scales than small ones. Peeling the outer scales of a large bulb converts its physiological responses into those of a small bulb, implying that the physiological discrepancies in plants developing from large or small bulbs are mediated by factors inherent to the bulb. We therefore performed broad analyses of the metabolite composition in the scales of bulbs subjected to temperature regimes affecting further plant development. We found a striking association between the level of glycerol, a primary metabolite mostly synthesized in the outer scales, and a delay in sprouting and flowering time, and reduction in abortion rate. Exogenous glycerol application to the bulbs before planting corroborated these results. Moreover, transcriptome analyses showed that flowering-promoting gene expression was downregulated in the bulb after glycerol treatment, while potential flowering inhibitor as well as a dormancy-related gene expressions were upregulated. Based on these studies, we postulate that glycerol is a major factor influencing both vegetative and reproductive development in lily.

Floral Induction and Flower Development of Orchids

Orchids comprise one of the largest, most highly evolved angiosperm families, and form an extremely peculiar group of plants. Various orchids are available through traditional breeding and micro-propagation since they are valuable as potted plants and/or cut flowers in horticultural markets. The flowering of orchids is generally influenced by environmental signals such as temperature and endogenous developmental programs controlled by genetic factors as is usual in many flowering plant species. The process of floral transition is connected to the flower developmental programs that include floral meristem maintenance and floral organ specification. Thanks to advances in molecular and genetic technologies, the understanding of the molecular mechanisms underlying orchid floral transition and flower developmental processes have been widened, especially in several commercially important orchids such as Phalaenopsis, Dendrobium and Oncidium. In this review, we consolidate recent progress in research on the floral transition and flower development of orchids emphasizing representative genes and genetic networks, and also introduce a few successful cases of manipulation of orchid flowering/flower development through the application of molecular breeding or biotechnology tools.

Functional and Evolutionary Characterization of the CONSTANS-like Family in Lilium�formolongi

Lilium�formolongi is a facultative long-day (LD) plant. Aiming to dissect the molecular regulation of the photoperiodic pathway, largely unknown in Lilium, we explored the CONSTANS/FLOWERING LOCUS T (CO/FT) module, a major regulatory factor in the external coincidence model of the photoperiodic flowering pathway in lily. We identified eight CONSTANS-LIKE (COL) family members in L.�formolongi, which could be divided into three types, according to their zinc-finger (B-box) protein domains. Type I included only LfCOL5, containing two B-box motifs. Type II contained six LfCOLs members that had only one B-box motif. Type III contained only LfCOL9 that showed a normal B-box and a second divergent B-box motif. Phylogenic analyses revealed that LfCOL5 was the closest to Arabidopsis CO. LfCOL5, LfCOL6 and LfCOL9 were up-regulated at the flowering induction stage under LDs, coinciding with an increase in LfFT1 expression. LfCOL5, LfCOL6 and LfCOL9 also showed obvious diurnal expression pattern for 3 d under LDs. However, under short-day (SD) conditions, the expression patterns of LfCOL5, LfCOL6 and LfCOL9 were variable and complex, with regard to the developmental stages and circadian rhythm. LfCOL5, LfCOL6 and LfCOL9 complemented the late flowering phenotype of the co mutant in Arabidopsis. Taken together, the results suggest that LfCOL5, LfCOL6 and LfCOL9 are involved in triggering flowering induction under LDs. LfCOL6 and LfCOL9 belong to types different from functional COL homologs in other plant species, illustrating the variation in phylogeny, evolution and gene function among LfCOL family members.

The metabolic (under)groundwork of the lily bulb toward sprouting

Large bulbs of Lilium longiflorum have an obligatory cold requirement to flower. Bulb cooling is widely used to induce and accelerate flowering. However, in-depth investigations of the effect of bulb cooling on major landmarks of plant development are lacking. It has been demonstrated that low temperature induces carbohydrate degradation, yet integrative studies on metabolic changes occurring in the bulb are not available. We detected that cold exposure mainly hastened bulb sprouting, rather than floral transition or blooming. Metabolite profiling of cooled and non-cooled bulbs was carried out, revealing cold-induced accumulation of soluble sugars, lipids and specific amino acids, and a significant reduction in tricarboxylic acid (TCA)-cycle elements. We observed that metabolic pathways located in the cytosol - including glycolysis, lipid synthesis and part of the gamma-Aminobutyric acid (GABA) shunt - were enhanced by cold exposure, while mitochondrial metabolism - namely the TCA cycle - was reduced by cold. We suggest a physiological model accounting for this metabolic discrepancy.

Overexpression of Lilium formosanumMADS-box (LFMADS) Causing Floral Defects While Promoting Flowering in Arabidopsis thaliana, Whereas Only Affecting Floral Transition Time in Nicotiana tabacum

The Formosa lily (Lilium formosanum) is one of the most common horticultural species in Taiwan. To explore gene regulation involved in this species, we used transcriptome analysis to generate PH-FB (mixed floral buds) and PH-LF (mature leaves) datasets. Combination of the PH-FB and PH-LF constructed a de novo assembly of the ALL dataset, including 18,041 contigs and 23,807 unigenes by Nr, GO, COG, and KEGG databases. The differential gene expression (DGE) analysis revealed 9937 genes were upregulated while 10,383 genes were downregulated in the developing floral buds compared to mature leaves. Seven putative genes (LFMADS1 to 7) encoding floral organ identity proteins were selected for further analysis. LFMADS1-6 genes were specifically expressed in the floral organ, while LFMADS7 in the floral buds and mature leaves. Phylogenetic analysis revealed that LFMADS1-3 is classified into B-class, LFMADS4 into C-class, LFMADS5 into D-class, and LFMADS6-7 into E-class, respectively. LFMADS-GFP fusion proteins appeared to localize in the nucleus, supporting their roles as transcription factors (TFs). Overexpression of the LFMADS2, LFMADS4, and LFMADS6 genes in Arabidopsis resulted in early flowering and floral defect, however, only early flowering in transgenic tobacco was observed. Highly expressed floral integrator genes, including AtFT, AtLFY, and AtFUL in transgenic Arabidopsis and NtFUL and NtSOC1 in transgenic tobacco, resulted in early flowering phenotype through qRT-PCR analysis. Yeast two-hybrid analysis suggested that LFMADSs may form higher order complexes with the B-, C-, D, and/or E-class proteins to determine the floral organ identity. Furthermore, E-class LFMADS proteins may function as a glue to mediate and strengthen the protein-protein interactions. Therefore, our de novo datasets would provide information for investigating other differentially expressed candidate transcripts. In addition, functional conservation of LFMADSs appears to be vital in floral transition and floral organ identity.

Trees, fungi and bacteria: tripartite metatranscriptomics of a root microbiome responding to soil contamination

BACKGROUND

One method for rejuvenating land polluted with anthropogenic contaminants is through phytoremediation, the reclamation of land through the cultivation of specific crops. The capacity for phytoremediation crops, such as Salix spp., to tolerate and even flourish in contaminated soils relies on a highly complex and predominantly cryptic interacting community of microbial life.

METHODS

Here, Illumina HiSeq 2500 sequencing and de novo transcriptome assembly were used to observe gene expression in washed Salix purpurea cv. 'Fish Creek' roots from trees pot grown in petroleum hydrocarbon-contaminated or non-contaminated soil. All 189,849 assembled contigs were annotated without a priori assumption as to sequence origin and differential expression was assessed.

RESULTS

The 839 contigs differentially expressed (DE) and annotated from S. purpurea revealed substantial increases in transcripts encoding abiotic stress response equipment, such as glutathione S-transferases, in roots of contaminated trees as well as the hallmarks of fungal interaction, such as SWEET2 (Sugars Will Eventually Be Exported Transporter). A total of 8252 DE transcripts were fungal in origin, with contamination conditions resulting in a community shift from Ascomycota to Basidiomycota genera. In response to contamination, 1745 Basidiomycota transcripts increased in abundance (the majority uniquely expressed in contaminated soil) including major monosaccharide transporter MST1, primary cell wall and lamella CAZy enzymes, and an ectomycorrhiza-upregulated exo-β-1,3-glucanase (GH5). Additionally, 639 DE polycistronic transcripts from an uncharacterised Enterobacteriaceae species were uniformly in higher abundance in contamination conditions and comprised a wide spectrum of genes cryptic under laboratory conditions but considered putatively involved in eukaryotic interaction, biofilm formation and dioxygenase hydrocarbon degradation.

CONCLUSIONS

Fungal gene expression, representing the majority of contigs assembled, suggests out-competition of white rot Ascomycota genera (dominated by Pyronema), a sometimes ectomycorrhizal (ECM) Ascomycota (Tuber) and ECM Basidiomycota (Hebeloma) by a poorly characterised putative ECM Basidiomycota due to contamination. Root and fungal expression involved transcripts encoding carbohydrate/amino acid (C/N) dialogue whereas bacterial gene expression included the apparatus necessary for biofilm interaction and direct reduction of contamination stress, a potential bacterial currency for a role in tripartite mutualism. Unmistakable within the metatranscriptome is the degree to which the landscape of rhizospheric biology, particularly the important but predominantly uncharacterised fungal genetics, is yet to be discovered.

Breeding of lilies and tulips-Interspecific hybridization and genetic background

Lilies and tulips (Liliaceae family) are economically very important ornamental bulbous plants. Here, we summarize major breeding goals, the role of an integrated method of cut-style pollination and fertilization followed by embryo rescue and mitotic and meiotic polyploidization involved in new assortment development. Both crops have been subjected to extensive interspecific hybridization followed by selection. Additionally, spontaneous polyploidization has played a role in their evolution. In lilies, there is a tendency to replace diploids with polyploid cultivars, whereas in tulip a majority of the cultivars that exist today are still diploid except for triploid Darwin hybrid tulips. The introduction of molecular cytogenetic techniques such as genomic in situ hybridization (GISH) permitted the detailed studies of genome composition in lily and tulip interspecific hybrids and to follow the chromosome inheritance in interspecific crosses. In addition, this review presents the latest information on phylogenetic relationship in lily and tulip and recent developments in molecular mapping using different DNA molecular techniques.

Identification of differentially expressed genes in flower, leaf and bulb scale of Lilium oriental hybrid 'Sorbonne' and putative control network for scent genes

BACKGROUND

Lily is an economically important plant, with leaves and bulbs consisting of overlapping scales, large ornamental flowers and a very large genome. Although it is recognized that flowers and bulb scales are modified leaves, very little is known about the genetic control and biochemical differentiation underlying lily organogenesis and development. Here we examined the differentially expressed genes in flower, leaf and scale of lily, using RNA-sequencing, and identified organ-specific genes, including transcription factors, genes involved in photosynthesis in leaves, carbohydrate metabolism in bulb scales and scent and color production in flowers.

RESULTS

Over 11Gb data were obtained and 2685, 2296, and 1709 differentially expressed genes were identified in the three organs, with 581, 662 and 977 unique DEGs in F-vs-S, L-vs-S and L-vs-F comparisons. By functional enrichment analysis, genes likely to be involved in biosynthetic pathways leading to floral scent production, such as 1-deoxy-D-xylulose-5-phosphate synthase (DXS), 3-ketoacyl-CoA thiolase (KAT), hydroperoxide lyase (HPL), geranylgeranyl pyrophosphate (GGPP) 4-hydroxy-3-methylbut-2-en-1-yl diphosphate (HDS) and terpene synthase (TPS), and floral color genes, such as dihydroflavonol 4-reductase (DFR), chalcone synthase (CHS), chalcone isomerase (CHI), flavonol synthase (FLS) were identified. Distinct groups of genes that participate in starch and sucrose metabolism, such as sucrose synthase (SS), invertase (INV), sucrose phosphate synthase (SPS), starch synthase (SSS), starch branching enzyme (SBE), ADP-glucose pyrophosphorylase (AGP) andβ-amylase (BAM) and photosynthesis genes (Psa, Psb, Pet and ATP) were also identified. The expression of six floral fragrance-related DGEs showed agreement between qRT-PCR results and RPKM values, confirming the value of the data obtained by RNA-seq. We obtained the open reading frame of the terpene synthase gene from Lilium 'Sorbonne', designated LsTPS, which had 99.55% homology to transcript CL4520.Contig5_All. In addition, 54, 48 and 50 differently expressed transcription factor were identified by pairwise comparisons between the three organs and a regulatory network for monoterpene biosynthesis was constructed.

CONCLUSIONS

Analysis of differentially expressed genes in flower, leaf and bulb scale of lily, using second generation sequencing technology, yielded detailed information on lily metabolic differentiation in three organs. Analysis of the expression of flower scent biosynthesis genes has provided a model for the regulation of the pathway and identified a candidate gene encoding an enzyme catalyzing the final step in scent production. These digital gene expression profiles provide a valuable and informative database for the further identification and analysis of structural genes and transcription factors in different lily organs and elucidation of their function.

Analysis of global gene expression profiles during the flowering initiation process of Lilium × formolongi

The onset of flowering is critical for the reproductive development of plants. Lilium × formolongi is a lily hybrid that flowers within a year after sowing. We successfully identified four important stages during vegetative growth and flowering initiation of L. × formolongi under long day conditions. The plant tissues from the four stages were used in a genome-wide transcriptional analysis to investigate stage-specific changes of gene expression in L. × formolongi. In total, the sequence reads of the four RNA-sequencing libraries were assembled into 52,824 unigenes, of which 37,031 (70.10%) were differentially expressed. The global expression dynamics of the differentially expressed genes were predominant in flowering induction phase I and the floral differentiation stage, but down-regulated in flowering induction phase II. Various transcription factor families relevant to flowering were elucidated, and the members of the MADS-box, SBP and CO-like transcription factor families were the most represented. There were 85 differentially expressed genes relevant to flowering. CONSTANS-LIKE, FLOWERING LOCUS T, TREHALOSE-6-PHOSPHATE SYNTHASE and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE homologs were discovered and may play significant roles in the flowering induction and transition process of L. × formolongi. A putative gene regulatory network, including photoperiod, age-dependent and trehalose-6-phosphate flowering pathways, was constructed. This is the first expression dataset obtained from a transcriptome analysis of photoperiod-mediated flowering pathway in lily, and it is valuable for the exploration of the molecular mechanisms of flowering initiation and the short vegetative stage of L. × formolongi.

Biotechnological Advancements for Improving Floral Attributes in Ornamental Plants

Developing new ornamental cultivars with improved floral attributes is a major goal in floriculture. Biotechnological approach together with classical breeding methods has been used to modify floral color, appearance as well as for increasing disease resistance. Transgenic strategies possess immense potential to produce novel flower phenotypes that are not found in nature. Adoption of Genetic engineering has supported the idea of floral trait modification. Ornamental plant attributes like floral color, fragrance, disease resistance, and vase life can be improved by means of genetic manipulation. Therefore, we witness transgenic plant varieties of high aesthetic and commercial value. This review focuses on biotechnological advancements in manipulating key floral traits that contribute in development of diverse ornamental plant lines. Data clearly reveals that regulation of biosynthetic pathways related to characteristics like pigment production, flower morphology and fragrance is both possible and predictable. In spite of their great significance, small number of genetically engineered varieties of ornamental plants has been field tested. Today, novel flower colors production is regarded as chief commercial benefit obtained from transgenic plants. But certain other floral traits are much more important and have high commercial potential. Other than achievements such as novel architecture, modified flower color, etc., very few reports are available regarding successful transformation of other valuable horticultural characteristics. Our review also summarized biotechnological efforts related to enhancement of fragrance and induction of early flowering along with changes in floral anatomy and morphology.

Molecular Regulation of Temperature-Dependent Floral Induction in Tulipa gesneriana

The vegetative-to-reproductive phase change in tulip (Tulipa gesneriana) is promoted by increasing temperatures during spring. The warm winters of recent years interfere with this process and are calling for new adapted cultivars. A better understanding of the underlying molecular mechanisms would be of help, but unlike the model plant Arabidopsis (Arabidopsis thaliana), very little is known about floral induction in tulip. To shed light on the gene regulatory network controlling flowering in tulip, RNA sequencing was performed on meristem-enriched tissue collected under two contrasting temperature conditions, low and high. The start of reproductive development correlated with rounding of the shoot apical meristem and induction of TGSQA expression, a tulip gene with a high similarity to Arabidopsis APETALA1 Gene Ontology enrichment analysis of differentially expressed genes showed the overrepresentation of genes potentially involved in floral induction, bulb maturation, and dormancy establishment. Expression analysis revealed that TERMINAL FLOWER1 (TgTFL1) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1-like1 (TgSOC1-like1) might be repressors, whereas TgSOC1-like2 likely is an activator, of flowering. Subsequently, the flowering time-associated expression of eight potential flowering time genes was confirmed in three tulip cultivars grown in the field. Additionally, heterologous functional analyses in Arabidopsis resulted in flowering time phenotypes in line with TgTFL1 being a floral repressor and TgSOC1-like2 being a floral activator in tulip. Taken together, we have shown that long before morphological changes occur in the shoot apical meristem, the expression of floral repressors in tulip is suppressed by increased ambient temperatures, leading either directly or indirectly to the activation of potential flowering activators shortly before the commencement of the phase change.

Virus-induced gene silencing (VIGS) in Lilium leichtlinii using the Cucumber mosaic virus vector

Lilies (Lilium) are among the most important floriculture crops, and to accelerate research regarding lily genetics, the development of reverse-genetics tools is necessary. However, Agrobacterium-mediated transformation in Lilium is time-consuming, since the plants require several years to progress from acclimation to flowering. Thus, virus-induced gene silencing (VIGS) is an attractive method for assaying gene function. In the present study, we modified a lily-derived strain of Cucumber mosaic virus (CMV-HL) as a VIGS vector and evaluated its effectiveness for inducing gene silencing in Lilium leichtlinii by introducing L. leichtlinii phytoene desaturase (LlPDS) gene fragments into an intercistronic region between the 3a and 3b genes of the CMV-HL RNA3 genome. At 30 days after inoculation (dpi) with LlPDS-containing CMV-HL, photo-bleaching was observed in the upper leaves of L. leichtlinii, and at 57 dpi, we observed that the natural orange color in flower tepals had faded. Reduced LlPDS expression and the detection of small interfering LlPDS RNA indicated that the color changes were the result of LlPDS gene silencing. In addition, the leaves also exhibited a mild photo-bleaching phenotype in the following year. Therefore, our results indicate that CMV-HL spreads systemically in the leaves and flowers of Lilium during the first year of infection, as well as in new shoots during the following year, and that the vector system can be successfully applied to induce short-term endogenous gene silencing in lilies.

Elucidating and mining the Tulipa and Lilium transcriptomes

Genome sequencing remains a challenge for species with large and complex genomes containing extensive repetitive sequences, of which the bulbous and monocotyledonous plants tulip and lily are examples. In such a case, sequencing of only the active part of the genome, represented by the transcriptome, is a good alternative to obtain information about gene content. In this study we aimed to generate a high quality transcriptome of tulip and lily and to make this data available as an open-access resource via a user-friendly web-based interface. The Illumina HiSeq 2000 platform was applied and the transcribed RNA was sequenced from a collection of different lily and tulip tissues, respectively. In order to obtain good transcriptome coverage and to facilitate effective data mining, assembly was done using different filtering parameters for clearing out contamination and noise of the RNAseq datasets. This analysis revealed limitations of commonly applied methods and parameter settings used in de novo transcriptome assembly. The final created transcriptomes are publicly available via a user friendly Transcriptome browser ( http://www.bioinformatics.nl/bulbs/db/species/index ). The usefulness of this resource has been exemplified by a search for all potential transcription factors in lily and tulip, with special focus on the TCP transcription factor family. This analysis and other quality parameters point out the quality of the transcriptomes, which can serve as a basis for further genomics studies in lily, tulip, and bulbous plants in general.

RNA-seq-based evaluation of bicolor tepal pigmentation in Asiatic hybrid lilies (Lilium spp.)

BACKGROUND

Color patterns in angiosperm flowers are produced by spatially and temporally restricted deposition of pigments. Identifying the mechanisms responsible for restricted pigment deposition is a topic of broad interest. Some dicots species develop bicolor petals, which are often caused by the post-transcriptional gene silencing (PTGS) of chalcone synthase (CHS) genes. An Asiatic hybrid lily (Lilium spp.) cultivar Lollypop develops bicolor tepals with pigmented tips and white bases. Here, we analyzed the global transcription of pigmented and non-pigmented tepal parts from Lollypop, to determine the main transcriptomic differences.

RESULTS

De novo assembly of RNA-seq data yielded 49,239 contigs (39,426 unigenes), which included a variety of novel transcripts, such as those involved in flavonoid-glycosylation and sequestration and in regulation of anthocyanin biosynthesis. Additionally, 1258 of the unigenes exhibited significantly differential expression between the tepal parts (false discovery rates <0.05). The pigmented tepal parts accumulated more anthocyanins, and unigenes annotated as anthocyanin biosynthesis genes (e.g., CHS, dihydroflavonol 4-reductase, and anthocyanidin synthase) were expressed 7-30-fold higher than those in non-pigmented parts. These results indicate that the transcriptional regulation of biosynthesis genes is more likely involved in the development of bicolor lily tepals rather than the PTGS of CHS genes. In addition, the expression level of a unigene homologous to LhMYB12, which often regulates full-tepal anthocyanin pigmentation in lilies, was >2-fold higher in the pigmented parts. Thus, LhMYB12 should be involved in the transcriptional regulation of the biosynthesis genes in bicolor tepals. Other factors that potentially suppress or enhance the expression of anthocyanin biosynthesis genes, including a WD40 gene, were identified, and their involvement in bicolor development is discussed.

CONCLUSIONS

Our results indicate that the bicolor trait of Lollypop tepals is caused by the transcriptional regulation of anthocyanin biosynthesis genes and that the transcription profile of LhMYB12 provides a clue for elucidating the mechanisms of the trait. The tepal transcriptome constructed in this study will accelerate investigations of the genetic controls of anthocyanin color patterns, including the bicolor patterns, of Lilium spp.

Current status and biotechnological advances in genetic engineering of ornamental plants

Cut flower markets are developing in many countries as the international demand for cut flowers is rapidly growing. Developing new varieties with modified characteristics is an important aim in floriculture. Production of transgenic ornamental plants can shorten the time required in the conventional breeding of a cultivar. Biotechnology tools in combination with conventional breeding methods have been used by cut flower breeders to change flower color, plant architecture, post-harvest traits, and disease resistance. In this review, we describe advances in genetic engineering that have led to the development of new cut flower varieties.

De novo assembly and analysis of the transcriptome of Ocimum americanum var. pilosum under cold stress

Background

Ocimum americanum var. pilosum is a chilling-sensitive, widely distributed plant that is consumed as a vegetable in central and southern China. To increase our understanding of cold stress responses in this species, we performed de novo transcriptome assembly for O. americanum var. pilosum and compared the transcriptomes of plants grown under normal and low temperatures.

Results

A total of 115,022,842 high quality, clean reads were obtained from four libraries (two replicates of control samples and two replicates of chilling-treated samples) and were used to perform de novo transcriptome assembly. After isoforms were considered, 42,816 unigenes were generated, 30,748 of which were similar to known proteins as determined by a BLASTx search (E-value < =1.0E-05) against NCBI non-redundant, Swiss-Prot, Gene Ontology, KEGG, and Cluster of COG databases. Comparative analysis of transcriptomes revealed that 5179 unigenes were differentially expressed (with at least 2-fold changes, FDR < 0.01) in chilling-treated samples, and that 2344 and 2835 unigenes were up- and down-regulated by chilling stress, respectively. Expression of the 10 most up-regulated and the five most down-regulated unigenes was validated by qRT-PCR. To increase our understanding of these differentially expressed unigenes, we performed Gene ontology and KEGG pathway enrichment analyses. The CBF-mediated transcriptional cascade, a well-known cold tolerance pathway, was reconstructed using our de novo assembled transcriptome.

Conclusion

Our study has generated a genome-wide transcript profile of O. americanum var. pilosum and a de novo assembled transcriptome, which can be used to characterize genes related to diverse biological processes. This is the first study to assess the cold-responsive transcriptome in an Ocimum species. Our results suggest that cold temperature significantly affects genes related to protein translation and cellular metabolism in this chilling sensitive species. Although most of the CBF pathway genes have orthologs in O. americanum var. pilosum, none of the identified cold responsive (COR) gene orthologs was induced by cold, which is consistent with the lack of cold tolerance in this plant.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2507-7) contains supplementary material, which is available to authorized users.

Identification of Flowering-Related Genes Responsible for Differences in Bolting Time between Two Radish Inbred Lines

Late bolting after cold exposure is an economically important characteristic of radish (Raphanus sativus L.), an important Brassicaceae root vegetable crop. However, little information is available regarding the genes and pathways that govern flowering time in this species. We performed high-throughput RNA sequencing analysis to elucidate the molecular mechanisms that determine the differences in flowering times between two radish lines, NH-JS1 (late bolting) and NH-JS2 (early bolting). In total, 71,188 unigenes were identified by reference-guided assembly, of which 309, 788, and 980 genes were differentially expressed between the two inbred lines after 0, 15, and 35 days of vernalization, respectively. Among these genes, 218 homologs of Arabidopsis flowering-time (Ft) genes were identified in the radish, and 49 of these genes were differentially expressed between the two radish lines in the presence or absence of vernalization treatment. Most of the Ft genes up-regulated in NH-JS1 vs. NH-JS2 were repressors of flowering, such as RsFLC, consistent with the late-bolting phenotype of NH-JS1. Although, the functions of genes down-regulated in NH-JS1 were less consistent with late-bolting characteristics than the up-regulated Ft genes, several Ft enhancer genes, including RsSOC1, a key floral integrator, showed an appropriate expression to the late-bolting phenotype. In addition, the patterns of gene expression related to the vernalization pathway closely corresponded with the different bolting times of the two inbred lines. These results suggest that the vernalization pathway is conserved between radish and Arabidopsis.