Next Article in Journal
Neuroprotective Strategies during Cardiac Surgery with Cardiopulmonary Bypass
Next Article in Special Issue
The Effect of High Pressure Techniques on the Stability of Anthocyanins in Fruit and Vegetables
Previous Article in Journal
Immunomodulatory Function of the Tumor Suppressor p53 in Host Immune Response and the Tumor Microenvironment
Previous Article in Special Issue
p-Hydroxyphenyl-pyranoanthocyanins: An Experimental and Theoretical Investigation of Their Acid—Base Properties and Molecular Interactions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 17
Issue 11
10.3390/ijms17111881
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies

by
Leifeng Xu
1,
Panpan Yang
1,2,
Suxia Yuan
1,
Yayan Feng
1,
Hua Xu
1,
Yuwei Cao
1 and
Jun Ming
1,*
1
Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
2
College of Landscape Architecture, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2016, 17(11), 1881; https://doi.org/10.3390/ijms17111881
Submission received: 2 September 2016 / Revised: 1 November 2016 / Accepted: 3 November 2016 / Published: 20 November 2016
(This article belongs to the Special Issue Anthocyanins)
Browse Figures
Versions Notes

Abstract

:
Lily tepals have a short lifespan. Once the tepals senesce, the ornamental value of the flower is lost. Some cultivars have attractive purple ovaries and fruits which greatly enhance the ornamental value of Asiatic hybrid lilies. However, little is known about the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries. To investigate the transcriptional network that governs purple ovary coloration in Asiatic hybrid lilies, we obtained transcriptome data from green ovaries (S1) and purple ovaries (S2) of Asiatic “Tiny Padhye”. Comparative transcriptome analysis revealed 4228 differentially expressed genes. Differential expression analysis revealed that ten unigenes including four CHS genes, one CHI gene, one F3H gene, one F3′H gene, one DFR gene, one UFGT gene, and one 3RT gene were significantly up-regulated in purple ovaries. One MYB gene, LhMYB12-Lat, was identified as a key transcription factor determining the distribution of anthocyanins in Asiatic hybrid lily ovaries. Further qPCR results showed unigenes related to anthocyanin biosynthesis were highly expressed in purple ovaries of three purple-ovaried Asiatic hybrid lilies at stages 2 and 3, while they showed an extremely low level of expression in ovaries of three green-ovaried Asiatic hybrid lilies during all developmental stages. In addition, shading treatment significantly decreased pigment accumulation by suppressing the expression of several unigenes related to anthocyanin biosynthesis in ovaries of Asiatic “Tiny Padhye”. Lastly, a total of 15,048 Simple Sequence Repeats (SSRs) were identified in 13,710 sequences, and primer pairs for SSRs were designed. The results could further our understanding of the molecular mechanisms of anthocyanin biosynthesis in Asiatic hybrid lily ovaries.

1. Introduction

Lily (Lilium spp.) is one of the most important ornamental plant species worldwide, owing to its large flowers, which have diverse tepal colors, unique shapes, and attractive floral scents. However, lily tepals have a short lifespan, which ends with wilting or abscission. Lily inflorescence longevity is approximately 2 weeks, and individual flower longevity is limited to approximately 1 week. Once the tepals senesce, the ornamental value of the flower is lost. However, some cultivars have purple ovaries, which develop into attractive purple fruits. These fruits maintain their color for approximately 2 months, which greatly enhances the ornamental value of these hybrid lilies. This ornamental trait can be developed for container-grown and outdoor-planted lily cultivars, and we believe that it has great potential to increase the ornamental value of hybrid lilies.
Anthocyanins, a class of flavonoid compounds, are the main pigments in flowers and fruits where they act as visual cues for pollinators and seed distributors [1]. Large variations in the intensity, hue, and patterning of anthocyanins are observed among angiosperm flowers [2]. In floricultural crops, these variations are one of the most important factors determining the market value of flowers. The biosynthetic pathway of anthocyanins has been studied extensively in different plant species and is one of the best-known pathways of specialized metabolism [3,4]. The genes in the anthocyanin biosynthetic include chalcone synthase (CHS), chalcone isomerase (CHI), flavonoid 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3'H), flavonoid 3′5′ hydroxylase (F3′5′H), dihydroflavonol 4-reductase (DFR), anthyocyanidin synthase (ANS), and UDP-glucose:flavonoid 3-O-glucosyltransferase (UFGT) [3,5]. The expression of these biosynthetic genes is mainly controlled at the transcriptional level and is regulated by a complex consisting of R2R3-MYB, basic helix–loop–helix (bHLH), and WD-repeat transcription factors [5,6]. Among these transcription factors, the R2R3-MYB transcription factors play a major role in determining the spatio-temporal expression of anthocyanin biosynthetic genes.
There have been some reports on the molecular mechanisms of flower colors and color patterning in Asiatic hybrid lilies in recent years. Three CHS genes (LhCHSA, LhCHSB, and LhCHSC) and one DFR gene (LhDFR) were isolated from Asiatic “Montreux” [7]. The full-length cDNA sequences of LhPAL, LhF3H, LhF3′H, and LhANS were isolated from the tepals of Asiatic lily “Montreux” [8], and two R2R3-MYB transcription factors, LhMYB6 and LhMYB12, were also isolated from Asiatic “Montreux” [9]. LhMYB12 was shown to directly bind to the promoters of LhCHSA and LhDFR to activate their expression [8]. In another study, LhMYB12-Lat, which contributes to the splatter pigmentation of the Tango Series cultivars of Asiatic hybrid lilies, was isolated from anthocyanin-accumulating tepals of Asiatic “Latvia” [2]. In addition to the large variations in flower colors and color patterning, Asiatic hybrids also show variations in ovary colors. Most Asiatic hybrids have green ovaries (e.g., “Yellow Pixels”), but a few (e.g., “Tiny Padhye”) have purple ones. However, little is known about the molecular mechanisms of purple ovary coloration in Lilium spp.
In the present study, the transcriptomes of stage 1 (green) and stage 2 (purple) ovaries of the purple-ovaried cultivar “Tiny Padhye” were sequenced using an Illumina Hiseq 4000 sequencer. After analyzing the data, we identified some key candidate genes related to anthocyanin biosynthesis in Lilium ssp. These transcriptome sequences may provide a valuable genomic resource to further our understanding of the molecular mechanisms of purple ovary coloration in Asiatic hybrid lilies.

2. Results

2.1. Measurement of Anthocyanin Content

A single anthocyanin pigment (cyanidin 3-O-β-rutinoside) was detected in the pigmented ovaries of the cultivar “Tiny Padhye”. During the two stages (S1 and S2) of ovary development of Asiatic “Tiny Padhye”, the ovaries underwent a rapid change in colors from green to purple (Figure 1b). This change was due to the increase of anthocyanin content from 0.05 to 4.75 μmol·g−1 fresh weight (FW) (Figure 2a). Microscopic examination of fresh hand-cut sections of ovaries of Asiatic “Tiny Padhye” showed that anthocyanins accumulated in both epidermal and sub-epidermal layers in purple ovaries (S2), but not in green ones (S1) (Figure 2b).

2.2. Sequencing and Sequence Assembly

The transcriptomes of green ovaries (S1) and purple ovaries (S2) of Asiatic “Tiny Padhye” were obtained by Illumina technology. Approximately 157 million and 145 million raw reads were generated from green ovaries and purple ovaries, respectively. After removing low-quality sequences, adapters, and ambiguous reads, we obtained approximately 153 million (S1) and 141 million (S2) high-quality clean reads (Table S1). More than 94% of bases had a q-value > 20 for both samples (Table S1). All raw reads were archived at NCBI Sequence Reads Archive (SRA) under accession number SRP081210. In total, approximately 294 million clean reads were used to assemble the transcriptome data using the Trinity method.

2.3. Identification of Differentially Expressed Genes

The genes showing differences in expression between the two samples were identified and filtered. A total of 4228 unigenes were identified as being significantly differentially expressed, of which 2270 showed higher expression in purple ovaries, and 1958 showed higher expression in green ovaries (Table S2 and Figure 3).

2.4. Functional Analysis of Differentially Expressed Genes

To further illustrate the main biological functions of the transcriptome, we conducted Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. A GO enrichment analysis provides a description of gene products in terms of their associated Biological Process (BP), Cellular Component (CC), and Molecular Function (MF) [10]. The obtained differentially expressed genes (DEGs) were categorized into 60 major functional groups, and metabolic process (GO:0008152), protein complex (GO:0043234), and catalytic activity (GO:0003824) were the most highly represented GO terms in BP, CC, and MF, respectively (Figure 4).
In the KEGG pathway analysis, the unigenes were mapped to 127 pathways, and 57 unigenes were involved in the anthocyanin biosynthesis pathway (Table 1). The DEGs up-regulated in purple ovaries were mapped to 102 KEGG pathways. The DEGs up-regulated in green ovaries were mapped to 94 KEGG pathways. In total, 19 of the DEGs up-regulated in purple ovaries were involved in the anthocyanin biosynthesis pathway (Table 1 and Table 2).

2.5. Expression Patterns of Genes Involved in Anthocyanin Biosynthesis Pathway

The anthocyanin biosynthesis pathway produced the purple anthocyanin pigment that gave the purple ovaries their color. Therefore, unigenes involved in the anthocyanin biosynthesis pathway were selected and studied in detail. We identified 57 unigenes encoding seven putative enzymes involved in anthocyanin biosynthesis from the transcriptomes. Among the 57 unigenes, 19 were DEGs. After further analysis of the expression patterns of these DEGs, we found that LhCHS1 (c124520_g1), LhCHS2 (c114453_g1), LhCHS3 (c122851_g2), LhCHS4 (c124520_g2), LhCHI (c123012_g2), LhF3H (c121630_g1), LhF3'H (c113932_g1), LhDFR (c124316_g1), LhUFGT1 (c113956_g1), and Lh3RT (c117809_g1) showed significantly higher transcript levels in purple ovaries than in green ones (Table 2 and Figure 5). The expression of most of the structural genes in the anthocyanin biosynthetic pathway is coordinately regulated by a ternary complex comprising three groups of transcription factors: R2R3-MYB, bHLH, and WDR transcription factors [5]. Thus, we analyzed the expression levels of regulatory factors that potentially control anthocyanin biosynthesis. We identified 36 DEGs encoding candidate transcription factors from the DEG data, including eight unigenes encoding R2R3-MYBs, 19 encoding bHLHs, and nine encoding WD40s (Table 1). Among these transcription factor unigenes, 19 were highly expressed in purple ovaries and 17 were highly expressed in green ones (Table 1).

2.6. Phylogenetic Analysis of R2R3-MYB Transcription Factors

To study the R2R3-MYB transcription factors that regulate anthocyanin biosynthesis in ovaries, we identified eight R2R3-MYB transcription factors among the DEGs. In the phylogenetic analysis, c51332_g1, c72124_g1, and c116256_g1 were in the same cluster as plant anthocyanin-promoting R2R3-MYBs, including AmROSEA1, AmROSEA2, and AtMYB75/90/113/114 (Figure 6). The FPKM value of LhMYB (c51332_g1) was significantly higher in purple ovaries than in green ones (Table 2). In the dendrogram, c36531_g1, c121270_g3 and c96386_g1 were grouped together, c116626_g1 was in another clade with AtMYB55, and c104994_g2 was grouped with AtCPC (Figure 6).

2.7. qPCR Validation of RNA-Seq Data

To validate the reliability of the comparative transcriptional data and to further analyze the differences in the expression profiles between purple-ovaried and green-ovaried Asiatic hybrid lilies during the three different development stages, 11 unigenes related to anthocyanin biosynthesis (LhCHS1, LhCHS2, LhCHS3, LhCHS4, LhCHI, LhF3H, LhF3'H, LhDFR, LhUFGT1, LhMYB, and Lh3RT) were selected for qPCR analyses. In the three purple-ovaried cultivars, all the chosen unigenes were highly expressed at stages 2 and 3, while their transcripts were almost undetectable or showed a low level of expression at stage 1 (Figure 7). These qPCR results were consistent with those obtained from the DGE expression profiling (Table 1 and Figure 7). In the three green-ovaried cultivars, transcripts of these unigenes showed an extremely low level of expression during all developmental stages (Figure 7). These qPCR results were consistent with the accumulation of anthocyanins in these cultivars (Figure 2a and Figure 7).

2.8. Effect of Shading Treatment on Pigment Accumulation and Gene Expression

To investigate whether shade treatment affects pigment accumulation and gene transcription, 2.50 cm Asiatic lily “Tiny Padhye” flower buds were covered with silver paper. Relative anthocyanin levels were very low in the shaded ovaries compared with those in the control group, indicating that the shade treatment significantly decreased the relative anthocyanin content. Furthermore, qPCR showed that expression of four anthocyanin structural genes, LhF3'H, LhDFR, LhUFGT1, and Lh3RT, was high in the control ovaries, but was extremely low or undetectable in the shaded ovaries (Figure 8). LhCHS1, LhCHI, LhF3H, and LhMYB expression levels were lower in the shaded ovaries than in the control ovaries, whereas LhCHS2, LhCHS3, and LhCHS4 expression in ovaries was not affected by the shade treatment (Figure 8).

2.9. Simple Sequence Repeat (SSR) Marker Identification and Characterization

In total, 15,048 SSRs were identified in 13,710 sequences (Table 3). Of these sequences, 1201 sequences contained more than one Simple Sequence Repeat (SSR) and 444 SSRs were present in the compound formation. Among the identified SSRs, the largest category of the SSRs was mononucleotide (8163, 54.2%), followed by dinucleotide (4446, 29.5%), trinucleotide (2309, 15.3%), tetranucleotide (96, 0.6%), hexanucleotide (19, 0.1%), and pentanucleotide (15, 0.1%) repeats. Using Primer 3.0 software, primer pairs for each SSR were designed (Table S4).

3. Discussion

3.1. Co-Expression of Genes Involved in Anthocyanin Biosynthesis Determines Anthocyanin Accumulation in Purple Ovaries

In this study, 57 unigenes involved in the anthocyanin biosynthetic pathway were selected and studied. To the best of our knowledge, this is the first time that all of the structural genes in the anthocyanin biosynthesis pathway have been identified in Lilium ssp. These data will help us to elucidate the regulatory mechanisms of anthocyanin biosynthesis in lily ovaries and tepals. There were several instances where more than one unigene was annotated as the same enzyme (Table 2), which might be because these unigenes represented different members of the same gene family or different fragments of a single transcript [11]. Among the 57 unigenes, we identified that four CHS genes, one CHI gene, one F3H gene, one F3'H gene, one DFR gene, one UFGT gene, and one 3RT gene were significantly up-regulated in purple ovaries, compared with green ovaries. We speculate that the high co-expression of these enzymes involved in anthocyanin biosynthesis is responsible for anthocyanin accumulation in the purple ovaries. Further studies are required to determine whether changes in the transcript levels of these candidate genes are related to anthocyanin metabolism.

3.2. LhMYB12-Lat Is Responsible for Purple Ovaries

In recent years, several studies have characterized regulatory genes involved in red and purple fruit coloration in eudicots and monocots. In apple (Malus × domestica), two MYB genes, MdMYB1 and MdMYBA, are responsible for anthocyanin biosynthesis in red-skinned apple cultivars [12,13]. Similarly, in grapevine (Vitis vinifera), VvMYBA1 and VvMYBA2 are responsible for anthocyanin biosynthesis in red-skinned cultivars [14]. GmMYB10, MrMYB1, PcMYB1, and LcMYB1 are involved in the anthocyanin biosynthetic pathway in mangosteen (Garcinia mangostana), Chinese bayberry (Myrica rubra), pear (Pyrus pyrifolia), and litchi (Litchi chinensis), respectively [15,16,17,18]. In maize (Zea mays) seeds, anthocyanin biosynthesis is regulated by C1 [19]. In this study, the phylogenetic analysis showed that c51332_g1, c72124_g1, and c116256_g1 were grouped into the same cluster with plant anthocyanin-promoting R2R3-MYBs, and the FPKM value of c51332-g1 was significantly higher in purple ovaries than in green ones (Table 2). In addition, the sequence of c51332_g1 was the same as that of LhMYB12-Lat, which contributes to the splatter pigmentation of the Tango Series cultivars of Asiatic hybrid lilies [2]. Therefore, we speculate that LhMYB12-Lat is the main anthocyanin-promoting protein that determines the distribution of anthocyanin pigments in Asiatic hybrid lily ovaries. These results indicate that, in both monocots and eudicots, R2R3-MYB genes have crucial roles in determining the distribution of anthocyanin pigments in fruits. Interestingly, in the Tango Series of Asiatic hybrid lily cultivars, some cultivars have green ovaries even though their tepals have splattered pigmentation due to the expression of LhMYB12-Lat. We suggest three possible explanations for this phenomenon. First, some negative regulators of the anthocyanin biosynthesis pathway may exist in ovaries of green-ovaried cultivars (e.g., a R2R3-MYB MYB repressor, a competitive R3-MYB repressor, or a competitive bHLH repressor). Alternatively, the expression of LhMYB12-Lat may be influenced by other factors such as microRNAs and/or NAC/SPL transcription factors. Last, methylation of the promoter of LhMYB12-Lat may affect its expression.

3.3. Shading Treatment Decreased Pigment Accumulation

Light is one of the important environmental factors that affect the anthocyanin biosynthesis [20]. Numerous studies have shown that light exposure can enhance the anthocyanin accumulation in fruit skin [21,22]. In addition to light intensity, light quality also influences the concentration of anthocyanins in fruits of many plant species. Anthocyanin accumulation is enhanced by ultraviolet-B light in apple [23,24,25] and red Chinese sand pear (P. pyrifolia) [26]. Blue light induces anthocyanin accumulation and expression of associated genes in apple fruits [27], strawberry (Fragaria × ananassa) fruits [28], and Chinese bayberry fruits [29]. In this study, compared with the control group, the shading treatment significantly decreased pigment accumulation by suppressing the expression of several anthocyanin structural genes and a related regulatory gene in ovaries of Asiatic “Tiny Padhye”. However, whether light quality and the wavelength of light affect anthocyanin biosynthesis in Asiatic lily ovaries remains unclear. This issue should be investigated further in future studies.

3.4. SSR Marker Identification and Characterization

Microsatellites or Simple Sequence Repeats (SSRs) are important for gene mapping and molecular breeding in plant research. Previously, approximate 5000 SSR markers had been reported in Lilium [30,31,32,33]. Considering the large and complicated lily genome, these makers are not sufficient for effective genetic mapping and marker-assisted selection breeding. Here, a total of 15,048 SSRs were identified in 13,710 sequences. The majority of SSRs were mononucleotide followed by dinucleotide and trinucleotide repeats. Using Primer 3.0, primer pairs for each SSR were designed. These SSRs can be used for genetic map development and molecular breeding within Lilium.

4. Materials and methods

4.1. Plant Materials

Three purple-ovaried Asiatic lily cultivars (“Tiny Padhye”, “Black Eye”, and “Dixie Jazz”) and three green-ovaried Asiatic lily cultivars (“Dot.com”, “Pink Pixels”, and “Yellow Pixels”) (Figure S1) were used in this study. Plants were grown in a greenhouse at the Chinese Academy of Agricultural Sciences (Beijing, China). Ovary sampling was performed at three different developmental stages: stage 1 (S1, bud length of about 2.50 cm); stage 2 (S2, 3 days post-anthesis); stage 3 (S3, 9 days post-anthesis) (Figure 1). Peels were sampled from ten ovaries at each stage and pooled to improve sample homogeneity and decrease sampling bias. Three biological replicates were used for each sample. All samples were immediately frozen in liquid nitrogen and stored at −80 °C.

4.2. Measurement of Anthocyanin Content

Anthocyanins were extracted with a solvent mixture of trifluoroacetic acid, methanol, methane acid, and water (1:70:2:27, v:v:v:v). The extract solution of each sample was then analyzed using HPLC according to the methods of Cunja et al. [34].

4.3. cDNA Library Construction, Sequencing, and Transcriptome Assembly

Total RNA was extracted from ovary wall using an RNAprep pure Plant Kit (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer’s protocol. A total of six RNA-seq libraries (three biological replicates for each sample) were generated according to the manufacturer’s protocol of the ScriptSeq mRNA-Seq Library Preparation Kit (Epicentre Biotechnologies, Madison, WI, USA) and were sequenced to generate 150 paired-end raw reads on an Illumina Hiseq 4000 platform. After clean data were generated according to [35], transcriptome assembly was performed using Trinity software as previously described for de novo transcriptome assembly without a reference genome [36].

4.4. Functional Annotation and Classification

The assembled unigenes were searched against public databases, including the NCBI Non-Redundant Protein Sequences (NR) database, NCBI Nucleotide Sequences (NT) database, Protein Family (PFAM) database, Eukaryotic Ortholog Groups (KOG) database, Swiss-Prot protein database, Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and Gene Ontology (GO) database. Gene Ontology (GO) enrichment analysis was performed on the differentially expressed genes (DEGs) using the GOseq R package [37].

4.5. Differentially Expressed Genes Analysis

The differentially expressed genes (DEGs) between green ovaries and purple ovaries were identified with the R package DESeq [38]. |log2 (foldchange)| > 1 and q-value < 0.01 were set as the thresholds for significantly different gene expression between the two samples. The heat maps showing expression profiles were generated based on the log10 (FPKM + 1) values using HemI (Heat map illustrator) with the default values [39].

4.6. Phylogenetic Analysis

A phylogenetic tree of selected R2R3-MYB transcription factors was constructed using MEGA 5 software by the neighbor-joining method with bootstrap analysis of 1000 replicates [40].

4.7. Quantitative Real-Time PCR Analyses

A total of 11 unigenes related to anthocyanin biosynthesis were chosen for quantitative real-time PCR (qPCR) analyses. qPCR analyses were performed using SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) (Takara, Dalian, China) and a Bio-Rad iQ5 Gradient Real Time PCR system with the following reaction conditions: one cycle of 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. The lily Actin gene was used as an internal control to normalize the expression data. Relative gene expression levels were calculated using the 2−ΔΔCt Method [41]. All primers used are listed in Table S3. Experiments were performed with three biological replicates and three technical replicates.

4.8. Effect of Shading Treatment on Pigment Accumulation and Gene Expression

To investigate the effect of shading treatment on pigment accumulation and gene transcription, 2.50 cm Asiatic lily “Tiny Padhye” flower buds were covered with silver paper. When ovary development reached stage 2, the paper was removed from the buds, and pigments and RNA were immediately extracted from the ovary wall. Plants exposed to light were used as the control group. Pigment measurement, RNA isolation, and qPCR were performed as described above.

4.9. Identification SSRs in ESTs and Primer Design

MISA (MIcroSAtellite; http://pgrc.ipk-gatersleben.de/misa) was used for microsatellite mining, and the minimum number of repeats used for selecting the SSRs was ten repeats for mononucleotide motifs and five repeats for dinucleotide, trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide motifs. Primer pairs for SSRs were designed using Primer 3.0 software (Premier Biosoft International, Palo Alto, CA, USA), with the following criteria: primer length ranging from 17 to 23 bp, melting temperature of 40–60 °C, GC content of 40%–60%, and PCR product size from 150 to 300 bp, and other parameter values as default.

5. Conclusions

In the present study, the transcriptomes of green ovaries (S1) and purple ovaries (S2) of the purple-ovaried cultivar “Tiny Padhye” were surveyed using the RNA-seq technology. Using this approach, ten putative anthocyanin structural genes and one R2R3-MYB TF were selected as the candidate genes most likely involved in determining the distribution of anthocyanins in Asiatic hybrid ovaries, and the expression level of these genes in green ovaries (S1) was much higher than purple ovaries (S2). In addition, shading treatment significantly decreased pigment accumulation by suppressing the expression of several unigenes related to anthocyanin biosynthesis in the ovaries of Asiatic “Tiny Padhye”. The transcriptome data from this study will facilitate our understanding of the molecular mechanisms of anthocyanin biosynthesis in the ovaries of Asiatic hybrid lilies.

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/1422-0067/17/11/1881/s1.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31672196) and the Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences. This research was conducted at the Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, China.

Author Contributions

Jun Ming designed the research. Leifeng Xu conducted the experiments. Leifeng Xu, Panpan Yang, Suxia Yuan, Yayan Feng, Hua Xu, and Yuwei Cao analyzed the data. Leifeng Xu wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davies, K.M.; Albert, N.W.; Schwinn, K.E. From landing lights to mimicry: The molecular regulation of flower colouration and mechanisms for pigmentation patterning. Funct. Plant Biol. 2012, 39, 619–638. [Google Scholar] [CrossRef]
  2. Yamagishi, M.; Toda, S.; Tasaki, K. The novel allele of the LhMYB12 gene is involved in splatter-type spot formation on the flower tepals of Asiatic hybrid lilies (Lilium spp.). New Phytol. 2014, 201, 1009–1020. [Google Scholar] [CrossRef] [PubMed]
  3. Winkel-Shirley, B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001, 126, 485–493. [Google Scholar] [CrossRef] [PubMed]
  4. Grotewold, E. The genetics and biochemistry of floral pigments. Annu. Rev. Plant Biol. 2006, 57, 761–780. [Google Scholar] [CrossRef] [PubMed]
  5. Koes, R.; Verweij, W.; Quattrocchio, F. Flavonoids: A colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005, 10, 236–242. [Google Scholar] [CrossRef] [PubMed]
  6. Ferreyra, M.L.F.; Rius, S.P.; Casati, P. Flavonoids: Biosynthesis, biological functions, and biotechnological applications. Front. Plant Sci. 2012. [Google Scholar] [CrossRef]
  7. Nakatsuka, A.; Izumi, Y.; Yamagishi, M. Spatial and temporal expression of chalcone synthase and dihydroflavonol 4-reductase genes in the Asiatic hybrid lily. Plant Sci. 2003, 165, 759–767. [Google Scholar] [CrossRef]
  8. Lai, Y.S.; Shimoyamada, Y.; Nakayama, M.; Yamagishi, M. Pigment accumulation and transcription of LhMYB12 and anthocyanin biosynthesis genes during flower development in the Asiatic hybrid lily (Lilium spp.). Plant Sci. 2012, 193, 136–147. [Google Scholar] [CrossRef] [PubMed]
  9. Yamagishi, M.; Shimoyamada, Y.; Nakatsuka, T.; Masuda, K. Two R2R3-MYB Genes, Homologs of petunia AN2, regulate anthocyanin biosyntheses in flower tepals, tepal spots and leaves of Asiatic hybrid lily. Plant Cell Physiol. 2010, 51, 463–474. [Google Scholar] [CrossRef] [PubMed]
  10. Berardini, T.Z.; Mundodi, S.; Reiser, L.; Huala, E.; Garcia-Hernandez, M.; Zhang, P.F.; Mueller, L.A.; Yoon, J.; Doyle, A.; Lander, G.; et al. Functional annotation of the Arabidopsis genome using controlled vocabularies. Plant Physiol. 2004, 135, 745–755. [Google Scholar] [CrossRef] [PubMed]
  11. Zhao, D.Q.; Jiang, Y.; Ning, C.L.; Meng, J.S.; Lin, S.S.; Ding, W.; Tao, J. Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall.). BMC Genom. 2014. [Google Scholar] [CrossRef] [PubMed]
  12. Takos, A.M.; Jaffe, F.W.; Jacob, S.R.; Bogs, J.; Robinson, S.P.; Walker, A.R. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 2006, 142, 1216–1232. [Google Scholar] [CrossRef] [PubMed]
  13. Ban, Y.; Honda, C.; Hatsuyama, Y.; Igarashi, M.; Bessho, H.; Moriguchi, T. Isolation and functional analysis of a MYB transcription factor gene that is a key regulator for the development of red coloration in apple skin. Plant Cell Physiol. 2007, 48, 958–970. [Google Scholar] [CrossRef] [PubMed]
  14. Walker, A.R.; Lee, E.; Bogs, J.; McDavid, D.A.J.; Thomas, M.R.; Robinson, S.P. White grapes arose through the mutation of two similar and adjacent regulatory genes. Plant J. 2007, 49, 772–785. [Google Scholar] [CrossRef] [PubMed]
  15. Palapol, Y.; Ketsa, S.; Kui, L.W.; Ferguson, I.B.; Allan, A.C. A MYB transcription factor regulates anthocyanin biosynthesis in mangosteen (Garcinia mangostana L.) fruit during ripening. Planta 2009, 229, 1323–1334. [Google Scholar] [CrossRef] [PubMed]
  16. Niu, S.S.; Xu, C.J.; Zhang, W.S.; Zhang, B.; Li, X.; Lin-Wang, K.; Ferguson, I.B.; Allan, A.C.; Chen, K.S. Coordinated regulation of anthocyanin biosynthesis in Chinese bayberry (Myrica rubra) fruit by a R2R3 MYB transcription factor. Planta 2010, 231, 887–899. [Google Scholar] [CrossRef] [PubMed]
  17. Feng, S.Q.; Wang, Y.L.; Yang, S.; Xu, Y.T.; Chen, X.S. Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10. Planta 2010, 232, 245–255. [Google Scholar] [CrossRef] [PubMed]
  18. Lai, B.; Li, X.J.; Hu, B.; Qin, Y.H.; Huang, X.M.; Wang, H.C.; Hu, G.B. LcMYB1 Is a Key Determinant of Differential Anthocyanin Accumulation among Genotypes, Tissues, Developmental Phases and ABA and Light Stimuli in Litchi chinensis. PLoS ONE 2014. [Google Scholar] [CrossRef] [PubMed]
  19. Cone, K.C.; Burr, F.A.; Burr, B. Molecular analysis of the maize anthocyanin regulatory locus C1. Proc. Natl. Acad. Sci. USA 1986, 83, 9631–9635. [Google Scholar] [CrossRef] [PubMed]
  20. Zhao, D.Q.; Tao, J. Recent advances on the development and regulation of flower color in ornamental plants. Front. Plant Sci. 2015. [Google Scholar] [CrossRef] [PubMed]
  21. Azuma, A.; Yakushiji, H.; Koshita, Y.; Kobayashi, S. Flavonoid biosynthesis-related genes in grape skin are differentially regulated by temperature and light conditions. Planta 2012, 236, 1067–1080. [Google Scholar] [CrossRef] [PubMed]
  22. Shinomiya, R.; Fujishima, H.; Muramoto, K.; Shiraishi, M. Impact of temperature and sunlight on the skin coloration of the “Kyoho” table grape. Sci. Hortic. 2015, 193, 77–83. [Google Scholar] [CrossRef]
  23. Ubi, B.E.; Honda, C.; Bessho, H.; Kondo, S.; Wada, M.; Kobayashi, S.; Moriguchi, T. Expression analysis of anthocyanin biosynthetic genes in apple skin: Effect of UV-B and temperature. Plant Sci. 2006, 170, 571–578. [Google Scholar] [CrossRef]
  24. Li, Y.Y.; Mao, K.; Zhao, C.; Zhao, X.Y.; Zhang, H.L.; Shu, H.R.; Hao, Y.J. MdCOP1 ubiquitin E3 ligases interact with MdMYB1 to regulate light-induced anthocyanin biosynthesis and red fruit coloration in apple. Plant Physiol. 2012, 160, 1011–1022. [Google Scholar] [CrossRef] [PubMed]
  25. Peng, T.; Saito, T.; Honda, C.; Ban, Y.; Kondo, S.; Liu, J.H.; Hatsuyama, Y.; Moriguchi, T. Screening of UV-B-induced genes from apple peels by SSH: Possible involvement of MdCOP1-mediated signaling cascade genes in anthocyanin accumulation. Physiol. Plant. 2013, 148, 432–444. [Google Scholar] [CrossRef] [PubMed]
  26. Zhang, D.; Qian, M.J.; Yu, B.; Teng, Y.W. Effect of fruit maturity on UV-B-induced post-harvest anthocyanin accumulation in red Chinese sand pear. Acta Physiol. Plant. 2013, 35, 2857–2866. [Google Scholar] [CrossRef]
  27. Li, Y.Y.; Mao, K.; Zhao, C.; Zhao, X.Y.; Zhang, R.F.; Zhang, H.L.; Shu, H.R.; Hao, Y.J. Molecular cloning and functional analysis of a blue light receptor gene MdCRY2 from apple (Malus domestica). Plant Cell Rep. 2013, 32, 555–566. [Google Scholar] [CrossRef] [PubMed]
  28. Kadomura-Ishikawa, Y.; Miyawaki, K.; Noji, S.; Takahashi, A. Phototropin 2 is involved in blue light-induced anthocyanin accumulation in Fragaria x ananassa fruits. J. Plant Res. 2013, 126, 847–857. [Google Scholar] [CrossRef] [PubMed]
  29. Shi, L.Y.; Cao, S.F.; Chen, W.; Yang, Z.F. Blue light induced anthocyanin accumulation and expression of associated genes in Chinese bayberry fruit. Sci. Hortic. 2014, 179, 98–102. [Google Scholar] [CrossRef]
  30. Lee, S.I.; Park, K.C.; Song, Y.S.; Son, J.H.; Kwon, S.J.; Na, J.K.; Kim, J.H.; Kim, N.S. Development of expressed sequence tag derived-simple sequence repeats in the genus Lilium. Genes Genom. 2011, 33, 727–733. [Google Scholar] [CrossRef]
  31. Shahin, A.; van Kaauwen, M.; Esselink, D.; Bargsten, J.W.; van Tuyl, J.M.; Visser, R.G.F.; Arens, P. Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa. BMC Genom. 2012. [Google Scholar] [CrossRef] [PubMed]
  32. Yuan, S.X.; Ge, L.; Liu, C.; Ming, J. The development of EST-SSR markers in Lilium regale and their cross-amplification in related species. Euphytica 2013, 189, 393–419. [Google Scholar] [CrossRef]
  33. Du, F.; Jiang, J.; Jia, H.M.; Zhao, X.Y.; Wang, W.H.; Gao, Q.K.; Mao, W.H.; Wu, Y.; Zhang, L.; Grierson, D.; et al. Selection of generally applicable SSR markers for evaluation of genetic diversity and identity in Lilium. Biochem. Syst. Ecol. 2015, 61, 278–285. [Google Scholar] [CrossRef]
  34. Cunja, V.; Mikulic-Petkovsek, M.; Stampar, F.; Schmitzer, V. Compound identification of selected rose species and cultivars: An insight to petal and leaf phenolic profiles. J. Am. Soc. Hortic. Sci. 2014, 139, 157–166. [Google Scholar]
  35. Liu, S.S.; Chen, J.; Li, S.C.; Zeng, X.; Meng, Z.X.; Guo, S.X. Comparative Transcriptome analysis of genes involved in GA-GID1-DELLA regulatory module in symbiotic and asymbiotic seed germination of Anoectochilus roxburghii (Wall.) Lindl. (Orchidaceae). Int. J. Mol. Sci. 2015, 16, 30190–30203. [Google Scholar] [CrossRef] [PubMed]
  36. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  37. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010. [Google Scholar] [CrossRef] [PubMed]
  38. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010. [Google Scholar] [CrossRef] [PubMed]
  39. Deng, W.K.; Wang, Y.B.; Liu, Z.X.; Cheng, H.; Xue, Y. HemI: A toolkit for illustrating heatmaps. PLoS ONE 2014. [Google Scholar] [CrossRef] [PubMed]
  40. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  41. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCt method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Ijms 17 01881 g001 550
Figure 1. Phenotype of ovaries of Asiatic “Yellow Pixels” (a) and “Tiny Padhye” (b) at three ovary developmental stages.
Figure 1. Phenotype of ovaries of Asiatic “Yellow Pixels” (a) and “Tiny Padhye” (b) at three ovary developmental stages.
Ijms 17 01881 g001
Ijms 17 01881 g002 550
Figure 2. Relative anthocyanin content and cross sections of ovaries. (a) The relative anthocyanin content of lily ovaries at different developmental stages. Error bars show the standard error (SE) of the mean; (b) Cross sections of ovaries of Asiatic “Tiny Padhye” at developmental stage 1 and stage 2. Scale bars: 100 μm.
Figure 2. Relative anthocyanin content and cross sections of ovaries. (a) The relative anthocyanin content of lily ovaries at different developmental stages. Error bars show the standard error (SE) of the mean; (b) Cross sections of ovaries of Asiatic “Tiny Padhye” at developmental stage 1 and stage 2. Scale bars: 100 μm.
Ijms 17 01881 g002
Ijms 17 01881 g003 550
Figure 3. The number of up- and down-regulated unigenes in a comparison of S2 to S1.
Figure 3. The number of up- and down-regulated unigenes in a comparison of S2 to S1.
Ijms 17 01881 g003
Ijms 17 01881 g004 550
Figure 4. The distribution of enriched Gene Ontology (GO) term categories of differentially expressed genes (DEGs).
Figure 4. The distribution of enriched Gene Ontology (GO) term categories of differentially expressed genes (DEGs).
Ijms 17 01881 g004
Ijms 17 01881 g005 550
Figure 5. Heatmap of DEGs related to anthocyanin biosynthesis pathway. The expression pattern of each unigene was shown using different grids from blue to red, with the values representing the log10 (FPKM + 1) values. FPKM standed for Fragment Per Kilobase of transcript per Million mapped reads.
Figure 5. Heatmap of DEGs related to anthocyanin biosynthesis pathway. The expression pattern of each unigene was shown using different grids from blue to red, with the values representing the log10 (FPKM + 1) values. FPKM standed for Fragment Per Kilobase of transcript per Million mapped reads.
Ijms 17 01881 g005
Ijms 17 01881 g006 550
Figure 6. Phylogenetic relationships of selected R2R3-MYB transcription factors in plants. The GenBank accession numbers are as follows: AmMYB308 (P81393), AmMYB340 (P81396), AmMYB306 (P81392), AmMYB305 (P81391), AmROSEA1 (ABB83826), AmROSEA2 (ABB83827), AmMYB330 (P81395), and all arabidopsis R2R3-MYB proteins were retrieved from TAIR (Available online: https://www.arabidopsis.org). Plant anthocyanin-promoting R2R3-MYBs were shown in the red circle.
Figure 6. Phylogenetic relationships of selected R2R3-MYB transcription factors in plants. The GenBank accession numbers are as follows: AmMYB308 (P81393), AmMYB340 (P81396), AmMYB306 (P81392), AmMYB305 (P81391), AmROSEA1 (ABB83826), AmROSEA2 (ABB83827), AmMYB330 (P81395), and all arabidopsis R2R3-MYB proteins were retrieved from TAIR (Available online: https://www.arabidopsis.org). Plant anthocyanin-promoting R2R3-MYBs were shown in the red circle.
Ijms 17 01881 g006
Ijms 17 01881 g007 550
Figure 7. Expression profiles of LhMYB and candidate anthocyanin structural genes in ovaries of Asiatic “Tiny Padhye”, “Black Eye”, “Dixie Jazz”, “Dot.com”, “Pink Pixels”, and “Yellow Pixels”. Error bars show the SE of the mean.
Figure 7. Expression profiles of LhMYB and candidate anthocyanin structural genes in ovaries of Asiatic “Tiny Padhye”, “Black Eye”, “Dixie Jazz”, “Dot.com”, “Pink Pixels”, and “Yellow Pixels”. Error bars show the SE of the mean.
Ijms 17 01881 g007
Ijms 17 01881 g008 550
Figure 8. Effect of shading treatment on gene expression in ovaries of Asiatic “Tiny Padhye”. Error bars show the SE of the mean.
Figure 8. Effect of shading treatment on gene expression in ovaries of Asiatic “Tiny Padhye”. Error bars show the SE of the mean.
Ijms 17 01881 g008
Table
Table 1. Candidate anthocyanin structural genes and transcription factor genes in purple ovaries.
Table 1. Candidate anthocyanin structural genes and transcription factor genes in purple ovaries.
GeneProteinNo. All aNo. Up bNo. Down c
CHSChalcone synthase38110
CHIChalcone isomerase110
F3HFlavanone 3-hydroxylase210
F3′HFlavonoid 3′-hydroxylase410
DFRDihydroflavonol 4-reductase310
ANSAnthocyanidin synthase300
UFGTAnthocyanidin 3-O-glucosyltransfersae530
3RTanthocyanidin-3-glucoside rhamnosyltransferase110
R2R3-MYBR2R3-MYB transcription factors2880
bHLHBasic helix–loop–helix transcription factors19811
WDWD-repeat proteins1936
a The total number of uni-transcripts analyzed; b The number of uni-transcripts with expression significantly up-regulated in purple ovaries compared with green ovaries; c The number of uni-transcripts with expression significantly down-regulated in purple ovaries compared with green ovaries.
Table
Table 2. Putative anthocyanin structural genes and transcription factor R2R3-MYBs identified among differentially expressed genes. Unigenes, showing significantly higher transcript levels in purple ovaries than in green ones, were highlighted in bold.
Table 2. Putative anthocyanin structural genes and transcription factor R2R3-MYBs identified among differentially expressed genes. Unigenes, showing significantly higher transcript levels in purple ovaries than in green ones, were highlighted in bold.
AnnotationUnigene IDS1 FPKMS2 FPKM
CHSc124520_g11.071255.60
c114453_g11.911495.41
c122851_g22.161996.29
c124520_g22.152926.54
c201597_g11.17106.51
c115882_g30.96614.86
c122851_g10.00342.93
c115882_g20.95375.32
c81738_g10.0012.03
c114453_g20.063.35
c124865_g10.071.22
CHIc123012_g293.54743.77
F3Hc121630_g125.791134.48
F3'Hc113932_g125.641099.95
DFRc124316_g18.435147.89
UFGTc113956_g15.673286.62
c114752_g18.03171.41
c92724_g213.5749.55
3RTc117809_g11.20584.93
R2R3-MYBc51332_g12.81264.36
c72124_g10.001.75
c116256_g12.2610.81
c36531_g12.477.92
c121270_g342.98256.59
c96386_g165.52182.05
c116626_g12.215.60
c104994_g20.0016.16
Table
Table 3. Statistics of the Simple Sequence Repeats (SSRs) identified in Lilium transcriptome.
Table 3. Statistics of the Simple Sequence Repeats (SSRs) identified in Lilium transcriptome.
ParameterValueParameterValue
SSR IdentificationDistribution to Different Repeat Type Classes
Total number of sequences examined295,787Mononucleotide8163 (54.2%)
Total size of examined sequences (bp)160,816,025Dinucleotide4446 (29.5%)
Total number of identified SSRs15,048Trinucleotide2309 (15.3%)
Number of SSR containing sequences13,710Tetranucleotide96 (0.6%)
Number of sequences containing more than one SSR1201Pentanucleotide15 (0.1%)
Number of SSRs present in compound formation444Hexanucleotide19 (0.1%)

Share and Cite

MDPI and ACS Style

Xu, L.; Yang, P.; Yuan, S.; Feng, Y.; Xu, H.; Cao, Y.; Ming, J. Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies. Int. J. Mol. Sci. 2016, 17, 1881. https://doi.org/10.3390/ijms17111881

AMA Style

Xu L, Yang P, Yuan S, Feng Y, Xu H, Cao Y, Ming J. Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies. International Journal of Molecular Sciences. 2016; 17(11):1881. https://doi.org/10.3390/ijms17111881

Chicago/Turabian Style

Xu, Leifeng, Panpan Yang, Suxia Yuan, Yayan Feng, Hua Xu, Yuwei Cao, and Jun Ming. 2016. "Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies" International Journal of Molecular Sciences 17, no. 11: 1881. https://doi.org/10.3390/ijms17111881

APA Style

Xu, L., Yang, P., Yuan, S., Feng, Y., Xu, H., Cao, Y., & Ming, J. (2016). Transcriptome Analysis Identifies Key Candidate Genes Mediating Purple Ovary Coloration in Asiatic Hybrid Lilies. International Journal of Molecular Sciences, 17(11), 1881. https://doi.org/10.3390/ijms17111881

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop