Floral pigmentation pattern in Oriental hybrid lily (Lilium spp.) cultivar 'Dizzy' is caused by transcriptional regulation of anthocyanin biosynthesis genes

Integrative Metabolomic and Transcriptomic Analyses Reveal the Mechanism of Petal Blotch Formation in Rosa persica

Petal blotch is a specific flower color pattern commonly found in angiosperm families. In particular, Rosa persica is characterized by dark red blotches at the base of yellow petals. Modern rose cultivars with blotches inherited the blotch trait from R. persica. Therefore, understanding the mechanism for blotch formation is crucial for breeding rose cultivars with various color patterns. In this study, the metabolites and genes responsible for the blotch formation in R. persica were identified for the first time through metabolomic and transcriptomic analyses using LC-MS/MS and RNA-seq. A total of 157 flavonoids were identified, with 7 anthocyanins as the major flavonoids, namely, cyanidin 3-O-(6″-O-malonyl) glucoside 5-O-glucoside, cyanidin-3-O-glucoside, cyanidin 3-O-galactoside, cyanidin O-rutinoside-O-malonylglucoside, pelargonidin 3-O-glucoside, pelargonidin 3,5-O-diglucoside, and peonidin O-rutinoside-O-malonylglucoside, contributing to pigmentation and color darkening in the blotch parts of R. persica, whereas carotenoids predominantly influenced the color formation of non-blotch parts. Zeaxanthin and antheraxanthin mainly contributed to the yellow color formation of petals at the semi-open and full bloom stages. The expression levels of two 4-coumarate: CoA ligase genes (Rbe014123 and Rbe028518), the dihydroflavonol 4-reductase gene (Rbe013916), the anthocyanidin synthase gene (Rbe016466), and UDP-flavonoid glucosyltransferase gene (Rbe026328) indicated that they might be the key structural genes affecting the formation and color of petal blotch. Correlation analysis combined with weighted gene co-expression network analysis (WGCNA) further characterized 10 transcription factors (TFs). These TFs might participate in the regulation of anthocyanin accumulation in the blotch parts of petals by modulating one or more structural genes. Our results elucidate the compounds and molecular mechanisms underlying petal blotch formation in R. persica and provide valuable candidate genes for the future genetic improvement of rose cultivars with novel flower color patterns.

LhANS-rr1, LhDFR, and LhMYB114 Regulate Anthocyanin Biosynthesis in Flower Buds of Lilium ‘Siberia’

The bulb formation of Lilium is affected by many physiological and biochemical phenomena, including flower bud differentiation, starch and sucrose accumulation, photoperiod, carbon fixation, plant hormone transduction, etc. The transcriptome analysis of flower buds of Lilium hybrid ‘Siberia’ at different maturity stages showed that floral bud formation is associated with the accumulation of anthocyanins. The results of HPLC-MS showed that cyanidin is the major anthocyanin found in Lilium ‘Siberia’. Transcriptome KEGG enrichment analysis and qRT-PCR validation showed that two genes related to flavonoid biosynthesis (LhANS-rr1 and LhDFR) were significantly up-regulated. The functional analysis of differential genes revealed that LhMYB114 was directly related to anthocyanin accumulation among 19 MYB transcription factors. Furthermore, the qRT-PCR results suggested that their expression patterns were very similar at different developmental stages of the lily bulbs. Virus-induced gene silencing (VIGS) revealed that down-regulation of LhANS-rr1, LhDFR, and LhMYB114 could directly lead to a decrease in anthocyanin accumulation, turning the purple phenotype into a white color. Moreover, this is the first report to reveal that LhMYB114 can regulate anthocyanin accumulation at the mature stage of lily bulbs. The accumulation of anthocyanins is an important sign of lily maturity. Therefore, these findings have laid a solid theoretical foundation for further discussion on lily bulb development in the future.

Wild and weedy Hesperis matronalis hosts turnip mosaic virus across heterogeneous landscapes in upstate New York

Turnip mosaic virus (TuMV) is a widespread and economically important pathogen in agricultural crops and has the widest known host range in the virus family Potyviridae. While management of the virus and its aphid vectors in agricultural fields decreases virus incidence, many alternative wild hosts for TuMV may serve as source populations for crop infection through spillover. Over thirty years ago, research demonstrated that the introduced brassica, Dame's Rocket (Hesperis matronalis) hosts several viruses, including TuMV. Here, we use both enzyme-linked immunosorbent assays (ELISA) and next generation sequencing to document the frequent infection by TuMV of Dame's Rocket, which is common and widespread in disturbed areas around crop fields in upstate New York. Deep sequencing of multiple tissue types of symptomatic hosts indicate that the infection is systemic and causes diagnostic, visible symptoms. In a common garden experiment using host populations from across upstate New York, we found evidence for genetic tolerance to TuMV infection in H. matronalis. Field surveys show that TuMV prevalence varies across populations, but is generally higher in agricultural areas. Examining disease dynamics in this and other common alternative hosts will enhance our understanding of TuMV epidemiology and, more broadly, virus distribution in wild plants.

Integrated mRNA and small RNA sequencing reveals a regulatory network associated with flower color in oriental hybrid lily

Anthocyanins are one of the main components of pigments, that are responsible for a wide range of colors in plants. To clarify the regulatory mechanism of anthocyanin biosynthesis in oriental hybrid lily, UPLC/MS analysis was performed to identify the pigments in two cultivars (white and pink). Four major anthocyanins were identified in pink cultivar, and no anthocyanins were detected in white cultivar. Transcriptome and small RNA sequencing (sRNAseq) analyses were performed using tepal tissues at two floral developmental stages from the two cultivars. In total, 55,698 transcripts were assembled, among which 233 were annotated as putative anthocyanin-related transcripts. Differential expression analysis and qRT-PCR results confirmed that most of the anthocyanin-related structural genes had higher expression levels in pink cultivar than in white cultivar. Conversely, LhANR showed a significantly high expression level in white cultivar. Annotated transcription factors (TFs), including MYB activators, MYB repressors and bHLHs, that putatively inhibit or enhance the expression of anthocyanin-related genes were identified. LhMYBA1, an anthocyanin activator, was isolated, and its heterologous expression resulted in a remarkably high level of anthocyanin accumulation. Additionally, 73 differentially expressed microRNAs (miRNAs), including 23 known miRNAs, were detected through sRNAseq. The miRNA target prediction showed that several anthocyanin-related genes might be targeted by miRNAs. Expression profile analysis revealed that these miRNAs showed higher expression levels at later floral developmental stages in white cultivar than in pink cultivar. The results indicated that anthocyanin deficiency in white cultivar might be influenced by multiple levels of suppressive mechanisms, including mRNAs and sRNAs.

Spatial patterning of scent in petunia corolla is discriminated by bees and involves the ABCG1 transporter

Floral guides are patterned cues that direct the pollinator to the plant reproductive organs. The spatial distribution of showy visual and olfactory traits allows efficient plant-pollinator interactions. Data on the mechanisms underlying floral volatile patterns or their interactions with pollinators are lacking. Here we characterize the spatial emission patterns of volatiles from the corolla of the model plant Petunia × hybrida and reveal the ability of honeybees to distinguish these patterns. Along the adaxial epidermis, in correlation with cell density, the petal base adjacent to reproductive organs emitted significantly higher levels of volatiles than the distal petal rim. Volatile emission could also be differentiated between the two epidermal surfaces: emission from the adaxial side was significantly higher than that from the abaxial side. Similar emission patterns were also observed in other petunias, Dianthus caryophyllus (carnation) and Argyranthemum frutescens (Marguerite daisy). Analyses of transcripts involved in volatile production/emission revealed lower levels of the plasma-membrane transporter ABCG1 in the abaxial versus adaxial epidermis. Transient overexpression of ABCG1 enhanced emission from the abaxial epidermis to the level of the adaxial epidermis, suggesting its involvement in spatial emission patterns in the epidermal layers. Proboscis extension response experiments showed that differences in emission levels along the adaxial epidermis, that is, petal base versus rim, detected by GC-MS are also discernible by honeybees.

Pathway-based analysis of anthocyanin diversity in diploid potato

Anthocyanin biosynthesis is one of the most studied pathways in plants due to the important ecological role played by these compounds and the potential health benefits of anthocyanin consumption. Given the interest in identifying new genetic factors underlying anthocyanin content we studied a diverse collection of diploid potatoes by combining a genome-wide association study and pathway-based analyses. By using an expanded SNP dataset, we identified candidate genes that had not been associated with anthocyanin variation in potatoes, namely a Myb transcription factor, a Leucoanthocyanidin dioxygenase gene and a vacuolar membrane protein. Importantly, a genomic region in chromosome 10 harbored the SNPs with strongest associations with anthocyanin content in GWAS. Some of these SNPs were associated with multiple anthocyanin compounds and therefore could underline the existence of pleiotropic genes or anthocyanin biosynthetic clusters. We identified multiple anthocyanin homologs in this genomic region, including four transcription factors and five enzymes that could be governing anthocyanin variation. For instance, a SNP linked to the phenylalanine ammonia-lyase gene, encoding the first enzyme in the phenylpropanoid biosynthetic pathway, was associated with all of the five anthocyanins measured. Finally, we combined a pathway analysis and GWAS of other agronomic traits to identify pathways related to anthocyanin biosynthesis in potatoes. We found that methionine metabolism and the production of sugars and hydroxycinnamic acids are genetically correlated to anthocyanin biosynthesis. The results contribute to the understanding of anthocyanins regulation in potatoes and can be used in future breeding programs focused on nutraceutical food.

R2R3-MYB genes control petal pigmentation patterning in Clarkia gracilis ssp. sonomensis (Onagraceae)

Petal pigmentation patterning is widespread in flowering plants. The genetics of these pattern elements has been of great interest for understanding the evolution of phenotypic diversification. Here, we investigate the genetic changes responsible for the evolution of an unpigmented petal element on a colored background. We used transcriptome analysis, gene expression assays, cosegregation in F₂ plants and functional tests to identify the gene(s) involved in petal coloration in Clarkia gracilis ssp. sonomensis. We identified an R2R3-MYB transcription factor (CgsMYB12) responsible for anthocyanin pigmentation of the basal region ('cup') in the petal of C. gracilis ssp. sonomensis. A functional mutation in CgsMYB12 creates a white cup on a pink petal background. Additionally, we found that two R2R3-MYB genes (CgsMYB6 and CgsMYB11) are also involved in petal background pigmentation. Each of these three R2R3-MYB genes exhibits a different spatiotemporal expression pattern. The functionality of these R2R3-MYB genes was confirmed through stable transformation of Arabidopsis. Distinct spatial patterns of R2R3-MYB expression have created the possibility that pigmentation in different sections of the petal can evolve independently. This finding suggests that recent gene duplication has been central to the evolution of petal pigmentation patterning in C. gracilis ssp. sonomensis.

Molecular cloning, characterization and expression analysis of LoTPS2 and LoTPS4 involved in floral scent formation in oriental hybrid Lilium variety 'Siberia'

Lilies are a commercially significant cut flower worldwide due not only to their elegant shape but also to their appealing scent. Among Lilium varieties, Lilium 'Siberia' is a cultivar that is prominent and highly favored by consumers due to its snowy white color and strong floral scent. Here, two terpene synthase genes (LoTPS2 and LoTPS4) that are responsible for floral scent production in Lilium 'Siberia' were cloned and functionally characterized. Recombinant LoTPS2 specifically catalyzed the formation of (E, E)-α-farnesene from FPP. Recombinant LoTPS4 is a multiproduct enzyme that produces D-limonene and β-myrcene as major volatile compounds and β-phellandrene, (+)-4-carene and 3-carene as minor products from GPP. Furthermore, LoTPS4 generates trans-α-bergamotene as a major product and di-epi-α-cedrene, α-cubebene and (E)-β-farnesene as minor compounds from FPP. Subcellular localization analysis using GFP fusion constructs revealed that LoTPS2 was localized in the cytosol, whereas LoTPS4 was localized in plastids. Real-time PCR analysis showed that LoTPS2 was highly expressed in the petals and sepals of the flower, while LoTPS4 was highly expressed in the filament of the flower. Moreover, mechanical wounding of flowers revealed that LoTPS2 showed a strong response to wounding via a rapid increase in its mRNA transcript level. Our results will assist scientists in exploring the molecular mechanisms of terpene biosynthesis in this species and will provide new insight into the biotechnological modification of the floral bouquet in Lilium.

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.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts

BACKGROUND

Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species.

RESULTS

The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration.

CONCLUSIONS

In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Repression of anthocyanin biosynthesis by R3-MYB transcription factors in lily (Lilium spp.)

Lily R3-MYB transcription factors are involved in negative regulation to limit anthocyanin accumulation in lily flowers and leaves and create notable color patterns on ectopically expressed petunia flowers. In eudicots, both positive and negative regulators act to precisely regulate the level of anthocyanin accumulation. The R3-MYB transcription factor is among the main factors repressing anthocyanin biosynthesis. Although, in monocots, the positive regulators have been well characterized, the negative regulators have not been examined. Two R3-MYBs, LhR3MYB1 and LhR3MYB2, which were identified in lily transcriptomes, were characterized in this study to understand the regulatory mechanisms of anthocyanin biosynthesis. LhR3MYB1 and LhR3MYB2 had a C2 suppressor motif downstream of a single MYB repeat; the similar amino acid motif appears only in AtMYBL2 among the eudicot R3-MYB proteins. Stable and transient overexpression of LhR3MYB1 and LhR3MYB2 in tobacco plants showed suppression of anthocyanin biosynthesis by both; however, suppression by LhR3MYB2 was stronger than that by LhR3MYB1. In the lily plant, the LhR3MYB2 transcript was detected in leaves with light stimulus-induced anthocyanin accumulation and in pink tepals. Although LhR3MYB1 was expressed in some, but not all tepals, its expression was not linked to anthocyanin accumulation. In addition, LhR3MYB1 expression levels in the leaves remained unchanged by the light stimulus, and LhR3MYB1 transcripts predominantly accumulated in the ovaries, which did not accumulate anthocyanins. Thus, although LhR3MYB1 and LhR3MYB2 have an ability to repress anthocyanin accumulation, LhR3MYB2 is more strongly involved in the negative regulation to limit the accumulation than that by LhR3MYB1. In addition, the overexpression of LhR3MYB2 generated notable color patterns in petunia flowers; thus, the usefulness of the LhR3MYB genes for creating unique color patterns by genetic engineering is discussed.

Proteomics Reveal the Profiles of Color Change in Brunfelsia acuminata Flowers

Brunfelsia acuminata is a popular ornamental plant with different colors resulted from the rapid change of color after blooming. The petals at day one (purple), day three (white and purple) and day five (white) were used to analyze the reason of flower color change by a comparative proteomics approach, gas chromatography coupled to a time-of-flight mass analyzer (GC-TOF-MS) and quantitative real-time PCR (qRT-PCR). The results showed that the 52 identified proteins were classified into eight functional groups, 6% of which were related to the anthocyanin metabolic pathway. The expression levels of all anthocyanin proteins from the first day to fifth day were remarkably down-regulated, which was consistent with the changing patterns of the key genes (CHS, CHI and F3′5′H) in petals. Simultaneously, the main floral volatile components including Linalool and 2-Hexenal (E) were identified, and the contents of 2-Hexenal at day five increased dramatically. Moreover, the content of flavonoids and total phenolic increased at day five. The majority of the proteins associated with stress defense and senescence proteins were up-regulated and the activities of peroxidase (POD), superoxide dismutase (SOD) and catalase (CAT) in the petals at day five were significantly higher than others. It was concluded that the competition in the precursors of metabolic pathways occurs and causes the flow of metabolite to the pathways of floral scent and lignin derived from the shikimate pathway or degrade into others. Therefore, the anthocyanin content significantly decreased, and the petal color changed from deep purple to white.