Next Article in Journal
The Role of EjSPL3, EjSPL4, EjSPL5, and EjSPL9 in Regulating Flowering in Loquat (Eriobotrya japonica Lindl.)
Next Article in Special Issue
Molecular Identification of Endophytic Fungi and Their Pathogenicity Evaluation Against Dendrobium nobile and Dendrobium officinale
Previous Article in Journal
Carbonic Anhydrase IX—Mouse versus Human
Previous Article in Special Issue
Comparative Analysis and Expression Patterns of the PLP_deC Genes in Dendrobium officinale
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 1
10.3390/ijms21010247
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

Comparative Transcriptomics Provides Insight into Floral Color Polymorphism in a Pleione limprichtii Orchid Population

by
Yiyi Zhang
1,
Tinghong Zhou
2,
Zhongwu Dai
1,
Xiaoyu Dai
1,
Wei Li
1,
Mengxia Cao
1,
Chengru Li
1,
Wen-Chieh Tsai
1,3,
Xiaoqian Wu
1,
Junwen Zhai
1,
Zhongjian Liu
1,* and
Shasha Wu
1,*
1
Key Laboratory of National Forestry and Grassland Administration for Orchid Conservation and Utilization, College of Landscape Architecture, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
Huanglong National Scenic Reserve, Songpan 623300, China
3
Institute of Tropical Plant Sciences and Microbiology, National Cheng Kung University, Tainan City 701, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(1), 247; https://doi.org/10.3390/ijms21010247
Submission received: 4 December 2019 / Revised: 25 December 2019 / Accepted: 26 December 2019 / Published: 30 December 2019
(This article belongs to the Special Issue Orchid Biochemistry)
Browse Figures
Review Reports Versions Notes

Abstract

:
Floral color polymorphism can provide great insight into species evolution from a genetic and ecological standpoint. Color variations between species are often mediated by pollinators and are fixed characteristics, indicating their relevance to adaptive evolution, especially between plants within a single population or between similar species. The orchid genus Pleione has a wide variety of flower colors, from violet, rose-purple, pink, to white, but their color formation and its evolutionary mechanism are unclear. Here, we selected the P. limprichtii population in Huanglong, Sichuan Province, China, which displayed three color variations: Rose-purple, pink, and white, providing ideal material for exploring color variations with regard to species evolution. We investigated the distribution pattern of the different color morphs. The ratio of rose-purple:pink:white-flowered individuals was close to 6:3:1. We inferred that the distribution pattern may serve as a reproductive strategy to maintain the population size. Metabolome analysis was used to reveal that cyanindin derivatives and delphidin are the main color pigments involved. RNA sequencing was used to characterize anthocyanin biosynthetic pathway-related genes and reveal different color formation pathways and transcription factors in order to identify differentially-expressed genes and explore their relationship with color formation. In addition, qRT-PCR was used to validate the expression patterns of some of the genes. The results show that PlFLS serves as a crucial gene that contributes to white color formation and that PlANS and PlUFGT are related to the accumulation of anthocyanin which is responsible for color intensity, especially in pigmented flowers. Phylogenetic and co-expression analyses also identified a R2R3-MYB gene PlMYB10, which is predicted to combine with PlbHLH20 or PlbHLH26 along with PlWD40-1 to form an MBW protein complex (MYB, bHLH, and WDR) that regulates PlFLS expression and may serve as a repressor of anthocyanin accumulation-controlled color variations. Our results not only explain the molecular mechanism of color variation in P. limprichtii, but also contribute to the exploration of a flower color evolutionary model in Pleione, as well as other flowering plants.

1. Introduction

Flower color is one of the most attractive characteristics of plants in nature. With such massive variation, flower color is regarded as an evolutionarily labile trait and has been shown to contribute to plant evolution [1,2,3]. In particular, flower color adaptive mutations mediated through pollinators are directly relevant to phenotypic evolution [4]. Color variation is considered a fixed difference between species and promotes the formation of population polymorphism [5,6]. Among angiosperms, sister species always display differences in flower hue and intensity, This pattern of flower color polymorphism is used as a model trait in the study of ecology, evolution, and gene regulation [7]. Color changes are related to flower pigment content. To date, the molecular mechanism of flower color transition has been investigated in several species, owing to the main floral pigments having been well characterized in many plants [8,9,10,11,12,13] providing sufficient information for studying floral color formation in non-model species and the opportunity to explore the relationship between phenotypic evolution and color variations.
Although flower color is influenced by many factors, flavonoids, especially anthocyanins which are produced by the anthocyanin biosynthesis pathway (ABP), are the primary components that contribute to floral pigments and they are produced by highly conservative structural and regulatory components [14]. The ABP involves multi-metabolic processes which mainly consist of seven core structural genes: CHS, CHI, F3H, F3′H, F3′5′H, DFR, and ANS, and several branch-enzyme genes [15]. Due to the instability of anthocyanidins, they exist mainly as anthocyanins, which are formed by anthocyanidins and various glycosides [16]. They play an irreplaceable role in the color development of plants and are primarily derived from three main anthocyanidins: Pelargonidin (brick red to scarlet), cyaniding (red to magenta), and delphinidin (purple to violet) [16]. Studies have shown that blocking the ABP can lead directly to variations in pigment production and affect flower color [17]. In addition to the structural genes in the ABP, transcription factors also contribute to flower color transition by regulating the spatial and temporal expression of the structural genes [18,19]. The ABP is regulated by three complex, interacting transcription factors: R2R3-MYB, basic helix–loop–helix (bHLH), and WD40-repeat (WDR) [20]. These transcription factors activate or suppress the transcription and expression of target genes, thereby regulating anthocyanin synthesis [21]. Generally, the structural and regulatory genes involved in the ABP have provided a number of targets to reveal the diversity of mutations that could block the ABP [22]. For flower polymorphism within populations, locating the blockage could elucidate the cause of flower color transition at the biochemical and molecular scales [8,10,11]. In addition, understanding their specific ABP is of benefit for predicting evolutionary influences from a genetic perspective.
The genus Pleione (Orchidaceae) comprises nearly 30 species of terrestrial, lithophytic, and epiphytic plants with high ornamental value [23]. There are 27 species in China, while Yunnan is the world biodiversity distribution center of this genus [24,25]. Pleione possesses different flower colors ranging between white, pink, lavender, magenta, light purple, and yellow [26]. In particular, populations of Huanglong, Sichuan Province, there remain a color polymorphism population, consisting of pink flowers of different intensities along with white mutant individuals, which can be considered an ideal situation to study the polymorphism formation mechanisms of Pleione, as well as benefit to explore potential correlation between color pattern and the species evolution.
The focus of our study was to understand the molecular mechanism of color polymorphism, including how the white individuals formed and the main reason caused pink flowers intensities, as well as summarize the general rules of the color distribution pattern. We aimed to investigate distribution of color monomorphic in the Huanglong P. limprichtii population and examine the transcriptome and biochemistry of their color polymorphic petals. RNA sequencing (RNA-seq) and ultra-performance liquid chromatography (UPLC) were used to identify the variation of related genes and the differences in flavonoid intermediates in the ABP that cause color transition, respectively.

2. Results

2.1. Color Differences and Quantity Distribution Pattern

After applying quantitative statistics to the individuals with the three distinct flowers at the full bloom stage randomly distributed in three rock populations, we concluded that the rose-purple individuals account for nearly 60%, the pink ones occupied about 30%, and the white ones not more than 10% (Table 1). The ratio of the number of the three differently flowered plants (rose-purple:pink:white) is roughly 6:3:1. Thus rose-purple is the main flower color in the Huanglong population, while white is the rarest color.
According to principle component analysis, there were significant differences in the a*(redness and greenness) values of the three colors (Figure 1. A comparison of L* (lightness) and C* (chroma) values among the three color groups indicated that C* values can also be used as an indicator to distinguish these three colors, since the rose-purple group has the highest C* value, the white group the lowest, and the pink group an intermediate value (Table 2).

2.2. Major Anthocyanin Compounds in P. limprichtii

UPLC analysis revealed that four anthocyanins and derivatives were detected in both the pigmented and white flower petals: Cyanidin 3-O-glucosyl-malonylglucoside, cyanidin 3-O-malonylhexoside, cyanidin, and delphinidin (Table 3). The cyanidin accumulation and delphinidin accumulation branches of the ABP in P. limprichtii indicate that anthocyanins and their derivatives are the main flower color pigments.

2.3. RNA-Seq and Annotation of Unigenes

To understand the molecular basis of flower polymorphism in P. limprichtii, three distinctly colored flower petals and lips were used for RNA-seq. A total of 80,525 unigenes with a mean length of 856 bp were obtained by de novo assembly. We assessed the quality of the unigenes in the transcriptome library. The length of N50 (sequence length of the shortest contig at 50%) was 1470 bp, the Q20 and Q30 percentages were 98% and 95%, the GC content was 40%, and the unigenes were generally distributed between 200 bp and 3000 bp in length. A total of 69,004,304 base pairs were aligned. These data show that the throughputs and sequencing quality were high enough to ensure further analysis.
According to the BLASTx results, a total of 33,724 unigenes were annotated, accounting for 41.88% of all the unigenes. Among these, 33,459 could be annotated using the Nr database (Non-redundant protein database, 41.55%), 21,177 using the Swiss-prot database (Swiss-Protein protein database, 26.30%), 11,067 using the KEGG database (Kyoto encyclopedia of genes and genomes, 24.68%), and 19,871 using the KOG database (Eukaryotic orthologous groups, 13.74%). In addition, there were 8396, 151, 42, and 23 unigenes annotated only using the Nr, Swiss-protein, KEGG, and KOG databases, respectively (Figure 2).
Statistical analysis of the E-value (The probability due to chance) characteristics of the Nr annotations revealed that 32.05% of the mapped sequences showed strong homology (E-value < 1 × 10-3), while 29.48% in the Swiss-protein database, 44.77% in the KEGG database, and 30.07% in the KOG database showed strong homology (Figure 3).
Based on the Nr annotations, 6495 unigenes were classified into 53 functional categories, belonging to three functional terms: molecular function, cellular component, and biological process. The largest percentages of unigenes identified within each of the three functional terms were metabolic process (3602 unigenes), cell and cell part (2237 and 2233 unigenes), and binding and catalytic activity (2481 and 3409 unigenes), corresponding to biological process, cellular component, and molecular function, respectively (Figure 4).
To exhaustively explore the potential functions of the annotated unigenes, 11,067 unigenes were mapped onto 131 KEGG pathways, including metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of antibiotics, and many other important metabolic pathways. Metabolic pathways and biosynthesis of secondary metabolites were enriched in the most unigenes, with 2450 and 1294 unigenes annotating to these two pathways, respectively, while the anthocyanin pathway was enriched in only one unigene. Among these 131 pathways, there are two color formation-related pathways, the flavonoid biosynthesis pathway (ko00941) to which 36 unigenes were mapped, and the ABP (ko00942) to which one unigene was mapped (the pathways are listed in Supplementary Table S3).

2.4. The Expression Patterns of ABP Genes

Genes involved in the flavonoid biosynthesis and ABPs that are related to color formation were isolated, including structural and regulatory genes, in order to identify those involved in the regulation of flower polymorphism in P. limprichtii. There were 21 ABP-related structural unigenes, which encoded 10 enzymes, expressed both in pigmented flowers and white flowers, with a total of 14 R2R3-MYB, 32 bHLH, and 23 WD40 regulatory unigenes. With regard to structural unigenes, six of the 10 belong to multigene families, with the exception of F3H, F3′5′H, ANR, and FLS, which were single copy genes (Table 4). All of these unigenes were used to analyze the expression pattern of the flower color polymorphism of P. limprichtii.
Since petal anthocyanins are detectable in different color petals, we inferred that floral color differences were caused by different expression patterns of ABP-related genes. The results show that the genes encoding flavonol synthese (FLS, PlFLS), anthocyanin synthese (ANS, PlANS1, PlANS2), and UDP-glucose anthocyanidin 3-O-glucosyltransferase (UFGT, PlUFGT1, PlUFGT2) were expressed at a higher level in pigmented flowers (rose-purple and pink) than in white flowers, and were more up-regulated in the rose-purple flowers than in the pink flowers (Figure 5a and Figure 6). These genes were correlated with flower color intensity and color phenotypes. Besides these directly affected genes, other genes that have different expression patterns between pigmented flowers and white flowers also influence color formation, such as the flavanone 3′-hydroxylase gene (F3′H; Pl F3′H3) and dihydroflavonol 4- reductase (DFR, PlDFR3). These were both up-regulated in rose-purple flowers and contribute to red color formation. According to the ABP-related gene expression patterns and metabolites detected in the three distinct flower groups, we drew a putative ABP of P. limprichtii.
The expression patterns of anthocyanin regulatory genes, including R2R3-MYB, bHLH, and WD40 were also investigated (Figure 5b–d). Phylogenetic analysis (Figure 7) shown that PlMYB13 (unigene0062421) and PlMYB4 (unigene0039181) were clustered with AtMYB75, AtMYB90, and AtMYB113, which belong to subgroup 6 of A. thaliana [27], and have been demonstrated to activate anthocyanin accumulation, while PlMYB10 (Unigene0058559) was homologous to AtMYB11, AtMYB12, and AtMYB111, which belong to subgroup 7 in A. thaliana, and have been suggested to control flavonol biosynthesis [28]. The expression pattern of PlMYB13 was consistent with anthocyanin accumulation, while PlMYB10 exhibited an inverse relationship between its expression and flower color intensity.

2.5. The Relationship between Structure Genes and TFs

Through the analysis of the expression patterns of ABP-related genes and the phylogenetic tree of R2R3-MYB, we obtained several candidate genes that correlate with floral color intensity. We constructed a co-expression network (Figure 8) to identify the interactions between ABP-related genes and MBW complex proteins (MYB, bHLH, and WDR). The results show that ten structure unigenes and five MBW complex transcriptional unigenes composed of two R2R3-MYB, two bHLH and one WD40 unigene exhibited interaction relationships. The expression pattern of PlMYB10 coincided with those of PlbHLH20 (Unigene0062784) and PlbHLH26 (Unigene00660), while PlbHLH20 and PlbHLH26 also correlated with PlWD40-1 (Unigene0002153). PlMYB10 also showed a linear relationship with PlFLS. We therefore inferred that PlMYB10 may interact with PlbHLH20 or PlbHLH26, and PlWD40-1 to form a MBW transcriptional complex, and regulate the expression pattern of PlFLS, finally affecting flower polymorphism in P. limprichtii.

3. Discussion

We investigated the distribution of color polymorphic individuals in three rock populations of P. limprichtii within Huanglong District and combined chemical detection and transcriptomic analysis to isolate the main pigment compounds and candidate genes that determine flower color intensity. We present a putative biosynthesis pathway and discuss the regulatory mechanisms of color formation.
To explore color variation formation factors, we used the CIELAB evaluation system to distinguish rose-purple, pink, and white flower color, and then counted the number of individuals of each of these phenotypes in each rock population. In all three populations the color distribution pattern was nearly consistent, with a color ratio of rose-purple 6:pink 3:white 1. This distribution pattern is very rare in natural color polymorphic populations. Studies have shown that intraspecific flower color variation is often attributed to genetic drift, pollination-mediated selection, environmental conditions, or herbivory [13,29,30,31]. In our experiments, the populations grow across a small range, with the whole population occurring on similar rocks and exposed to the same climatic conditions. Environmental elements such as temperature, drought stress, and ultraviolet radiation were therefore not considered to be crucial promoters of color variation. Pollinator-mediated selection plays an important role in color variation, especially for deceptive pollination species in which competition for pollinators in sympatry promotes flower color divergence [32], and shifts in pollinators also contribute to the macro-evolution of flowers color [22]. Thus, we inferred that the flower color polymorphism within these populations might have been induced by pollinators. Color polymorphism may be a consequence of pollination competition or specific adaptations to pollinators, and pollinator behavior exerts strong selection stress on color variations. Some research also showed adaptive selection for pigmented flowers because colored flowers are less likely to be disrupted by herbivories than colorless ones [3]. According to our field observations, we found that white-flowered individuals were more susceptible to damage than individuals with pigmented flowers, and that white flower petals and cores were severely foraged when blooming. We thus inferred that the dominant pigmented color was beneficial to avoid herbivory, and reduce damage by herbivores to the population. When individuals were damaged they also suffered reduced attractiveness to pollinators, which is not conductive to the stability and development of the population. This explains the distribution pattern of the number of rose-purple flowers in the population. There is another view that flower color may not be the primary goal of natural selection, nor the initial choice of pollinators. Indeed, the biosynthetic precursors of pigments not only display color variations, but also serve other physiological functions [14]. Studies have shown that secondary metabolites associated with plant defense functions share the same biosynthetic pathway, the flavonoid synthesis pathway, correlating pigment with defense ability [33]. Therefore, colorful-flowered individuals were more resistant to some adversities. In summary, the phenomenon that rose-purple flowers were frequent and white ones rare within the population may be mediated by the pollinators and herbivores, and also related to survival adaptability of the P. limprichtii. The 6:3:1 distribution pattern of color polymorphism might be a reproductive strategy for the population to maintain the maximum population density, but further evidences should be investigated.
In theory, the transition from pigmented to white flowers could involve any mutations that block one or more steps in the anthocyanin pathway. This includes loss-of-function mutations in any pathway enzyme-coding genes, as well as the cis-regulatory genes that influence any of the pathway enzymes [34,35]. In our study, expression analysis identified several obvious differentially expressed genes in the petal which were down-regulated in white samples compared to pigmented samples, but metabolite detection found that Cy- and Del-derivatives existed in both white petals and pigmented petals, indicating that the color variations, especially the white petals, do not lacking any steps of the anthocyanin pathway. The cis-regulation of transcription factors is a crucial element to promote color divergence. This result is similar with the white color formation in Primula vulgaris which is caused by different genes expression pattern rather than loss- of-function mutations leading to the lack of anthocyanin [36].
For ABP-related genes, we isolated 21 transcripts which encode ten enzymes. Seven of the ten were flavonoid synthase genes, including PlCHS, PlCHI, PlF3H, PlF3′H, PlF3′5′H, PlDFR, and PlANS; one was a proanthocyanidin synthase gene, PlANR; one was a flavonol synthase gene, PlFLS; and one was an anthocyanin synthase gene, PlUFGT. Most of them are multi-gene families, only PlF3H, PlF3′5′H, PlFLS, and PlANR are single copy. To clearly illustrate the catalyzation steps of the ABP, we regarded Cy-related and Del-related biosynthesis processes as independent branches [8]. Thus we did not have to measure the content of each anthocyanin compound. Each branch of anthocyanin synthase was considered to make an equivalent contribution to pigmentation. Considering that both the Cy-related and Del-related branches may share the majority of enzymes, here, we selected ABP-related genes on the Cy-related branches to analyze their expression patterns.
The expression analysis of early step structural genes revealed a high level of PlCHS and PlCHI expression in white petals and a low level in pigmented petals, suggesting that white petals can produce a large amount of naringenin but cannot eventually flux this into anthocyanin synthase. Meanwhile, the expression levels of PlF3′H, PlF3′5′H, PlDFR, PlANS, and PlUFGT in pigmented flower petals was high, compared with white ones, and we inferred that downstream structural genes make a large contribution to coloration. Analysis also showed that PlFLS was significantly more upregulated in white petals than in pigmented petals. FLS encoding enzymes lead substrate into the flavones and flavonols pathway [37]. It has been suggested that the competition between the anthocyanin synthesis pathway and the flavone and flavonols pathways mainly results in substrate competition between FLS with DFR, while the FLS enzyme strengthens the metabolic flux toward the flavonols and limits anthocyanin accumulation [38]. This situation has also been reported in other species, such as in Paeonia ostii, a higher expression of PoFLS4 in the nearly white flowers promotes dihydroflavonols transition into flavonols [11]. In onions, enhanced AcFLS could maximize flavonol production in the sheath [39]. Finally, in Muscari armeniacun, the conversion of substrate between FLS and DFR facilitates the elimination of blue pigmentation [8]. Thus, we confirmed PlFLS as one of the candidate genes for white color formation in P. limprichtii. The up-regulation of PlDFR in pigmented flowers is closely accompanied by a decrease of PlFLS; hence, more dihydroflavonols flow into anthocyanin production in pigmented flowers. From our analysis, the expression patterns of PlANS and PlUFGT are correlated with color intensity, and they showed their highest expression levels in the rose-purple flowers and their lowest in the white flowers. The ANS (Anthocyanidin synthase) encoding enzyme catalyzes the conversion of colorless leucocyanidin into colored anthocyanin [40], and anthocyanin is further glycosylated by different UFGT (UDP flavonoid glucosyl transferase) encoding enzymes that convert the anthocyanidins to different anthocyanin derivatives, exhibiting the final color [41]. We therefore speculate that PlANS and PlUFGT are two crucial genes that determine color intensity in P. limprichtii.
It has been revealed that the MBW protein complex, a combination of R2R3-MYB and bHLH transcription factors, along with WD40 proteins, play an important role in regulating the transcription of structural genes [41,42,43,44,45]. The activities of R2R3-MYB factors have distinct roles in determining the action of the complexes, either to promote or inhibit the transcription of anthocyanin biosynthesis genes [46,47]. By combining R2R3-MYB phylogenetic and co-expression network analyses, we isolated PlMYB10, which was homologous with the AtMYB11, AtMYB12, and AtMYB111 belonging to S7 in Arabidopsis and that have been demonstrated to contribute to the regulation of genes that account for anthocyanin accumulation in all tissues [27,28]. The expression pattern was high expression in nearly white flower petals, gradually reducing in flowers with increasing color intensity. It also corresponded with the color polymorphism phenotypes. In other species, overexpression of AmMYB330, a negative regulator of the flavonoid biosynthesis, has been proven to inhibit phenylpropanoid metabolism in transgenic tobacco (Nicotiana tabacum) plants [48]. The co-expression network showed that PlbHLH20 and PlbHLH26, along with PlWD40-1 have a strong relationship to PlMYB10. It is likely that this potential MBW complex PlMYB10/PlbHLH20/PlWD40-1 or PlMYB10/PlbHLH26/PlWD40-1 may serve as a repressor responsible for the variation in color intensity in P. limprichtii. Co-expression network revealed PlFLS is the most likely target gene interacting with PlMYB10. Previous studies have verified that R2R3-MYB can regulate the expression pattern of FLS through the overexpression of PsMYB114L (from Paeonia suffruticosa) in Arabidopsis [49], which is consistent with our results. Thus, we tentatively speculate about the ABP of P. limprichtii (Figure 9). Further studies examining sequencing differences in these candidate genes and PlMYB10 are necessary to assess our speculations. Metabolic substance quantification is also necessary to confirm the leading pigments in P. limprichtii.

4. Materials and Methods

4.1. Plant Materials

Pleione limprichtii, belongs to the Orchidaceae family, and is distributed in the Huanglong Nature Reverse in Sichuan Province, China, with three distinct color variations including rose-purple, pink, and white (Figure 10a–c). The study site was located at 32.68°N, 104.04°E, and 1851.61 m above sea level (Figure 10d). The plants are generally found living on rocks along streams. The sexual reproduction of Pleione is commonly dependent on deceptive pollination and each individual usually grows only one flower. The majority of individuals in the population have rose-purple flowers, while white flowers are relatively less common.

4.2. Quantitative Statistics and Flower Colorimeter Analysis

Three rocks with P. limprichtii populations which contained individuals with all three color variations were selected, the number of each color individuals were counted. In addition, in order to evaluate the flower color objectively, a hand spectrophotometer (CS-280, Hangzhou Color Spectrum Technology Co., Ltd., Hangzhou, China) was used to measure the CIE L*a*b* color components with five technical repetitions using five different parts of the selected petals, in order to evaluate the flower color objectively. The formula C* = a ∗ 2   +   b ∗ 2 ) , the L* (lightness) and C* (chroma) through the petals were measured using the lightness coefficient ‘a*’, which indicates greenness to redness as the value increases from negative to positive, and ‘b*’, which represents blueness to yellowness. Principle component analysis was performed using the a* and b* values as the principle components [32].

4.3. Sampling and RNA Extraction and cDNA Synthesis

Petals from rose-purple, pink, and white flowers at the full bloom stage were sampled, immediately frozen in liquid nitrogen, and preserved at −80 °C before pigment analysis and RNA extraction. RNA was isolated using a Quick RNA isolation Kit (Huayueyang Co., Ltd., Beijing, China), according to the manufacturer’s instructions. The concentration and purity of the RNA was measured with a NanoDrop 2000 (Thermo Fisher Scientific Co., Ltd., Waltham, MA, USA) and Agilent 2100 (Agilent Technologies Co., Ltd., Palo Alto, CA, USA) to verify RNA integrity. A total of nine samples, including three biological replicates, with high concentrations of RNA for each of the three color morphs were selected, then the strand cDNA synthesis was performed using a Revert Aid First Strand cDNA Synthesis Kit (Thermo Fisher, Foster City, CA, USA), according to the manufacturer’s instructions, and was stored at −80 °C for RT-qPCR assays.

4.4. Measurement of Flower Anthocyanin

Sample preparation and extraction methods were as follow, the freeze dried tissues were crushed using a mixer mill (MM 400, Verder Shanghai Instruments and Equipment Co., Ltd., Shanghai, China) with a zirconia bead for 1.5 min at 30 Hz. A 100 mg sample of the powder was weighed and extracted overnight at 4 °C with 1.0 mL 70% aqueous methanol. Following centrifugation at 10, 000× g for 10 min, the extracts were absorbed and filtered. Then the sample extracts were analyzed using a UPLC system, the analytical conditions were as follow, HPLC: column, Waters ACQUITY UPLC HSS T3 C18 (1.8 µm, 2.1 mm × 100 mm); solvent system, water (0.04% acetic acid): acetonitrile (0.04% acetic acid); gradient program, 95:5 v/v at 0 min, 5:95 v/v at 11.0 min, 5:95 v/v at 12.0 min, 95:5 v/v at 12.1 min, 95:5 v/v at 15.0 min; flow rate,. 0.40 mL/min; temperature, 40 °C; injection volume: 2 μL [32,50]. The mass spectrometry data was analyzed using the software Analyst 1.6.1 (AB Sciex Pte. Ltd., Concord, Ontario, Canada) and based on the Metware Database (MWDB, Metware Biotechnology Co., Ltd., Wuhan, China), a local self-built database, and a public metabolite information database, to identify anthocyanin compounds. The experiments were repeated three times. 1.8 µm

4.5. Library Preparation and Sequencing

The construction of the libraries and RNA-seq were performed by the Omicshare Biotechnology Corporation (Guangzhou, China). The mRNA was enriched with Oligo (dT) beads, then broken into short fragments and reverse transcribed into cDNA as templates. First and second-strand cDNA were then synthesized. The cDNA fragments were purified with a QiaQuick PCR extraction kit (Qiagen, Valencia, CA, USA) and end repaired. Poly (A) was added and ligated to Illumina sequencing adapters. They were then sequenced using the Illumina HiSeqTM 4000 (Illumina Inc., Centre Drive, San Diego, CA, USA).

4.6. De Novo Transcriptome Assembly Annotation

Transcriptome de novo assembly was performed with clean data, filtered from the raw data by removing adaptors and unknown nucleotides (>10%), and those with low quality reads. The data were assembled using the Trinity platform [51] with the parameters ‘K-mer = 31′ and ‘K-mer cover = 6′. First, short reads of a certain length were combined with overlap to form longer contigs. Then, based on their paired-end information, clean reads were mapped back to the corresponding contigs. Thus the sequences of the transcripts were finished, and defined as unigenes. All assembled unigenes were then annotated using BLASTx (E-value ≤ 1 × 10-5) against protein databases, including the National Center for Biotechnology Information non-redundant (Nr, ftp://ftp.ncbi.nih.gov/blast/db/), Swiss-Protein (https://www.uniprot.org/), Kyoto Encyclopedia of Genes and Genomes (KEGG, https://www.genome.jp/kegg/), and Gene Ontology (GO, http://geneontology.org/) databases. while Nr and Swiss-Protein annotate gene function, KEGG is used to understand biological systems and GO divides genes into different categories. BLASTx was used to search for the unigenes against the public databases with the following order of priority: Nr, Swiss-protein, KEGG, and COG. When a unigene could not be aligned to any of these protein databases, the protein code sequence and sequence direction was confirmed using the ESTscan program.

4.7. Expression Profile and RT-qPCR

To compare color-related unigenes expression divergence, they were first normalized to RPKM (reads per kb per million reads). After that, different sample ratios of RPKM values were calculated. The FDR (false discovery rate) value was used to identify the threshold of the p-value in multiple tests in order to compute the significance of the differences among unigenes. Here, only FDR significance scores < 0.05, and log2 ratios > 1 were regarded as differentially expressed genes and used in subsequent analysis. In order to validate the expression pattern of the RNA-Seq results, ten important different expression genes (the primers used are listed in Table S1) were selected and measured using RT-qPCR on a Quant Studio 5 Real-Time PCR System (Thermo Fisher, Foster City, CA, USA) using the PowerUpTM SYBRTM Green Master Mix (Thermo Fisher, Foster City, CA, USA), according to the manufacturer’s instructions. PlUBC34 and PlUBC35 (ubiquitin-conjugating enzyme; primers are also listed in Table S1) actin genes from Pleione were used as the internal control for the normalization of gene expression. Each sample (including three biological repetitions) was quantified in triplicate.

4.8. Genes Related to the ABP and Phylogenetic Analyses

The ABP-related structure genes and transcription factors including CHS, CHI, F3H, F3′5′H, DFR, FLS, ANS, FNS, UFGT, MYB, bHLH, and WD40 were used as queries to retrieve the corresponding unigenes from the P. limprichtii libraries. Meanwhile, every reference gene from model plants obtained from KEGG were aligned with the libraries using BLASTx to search for more related unigenes.
For the MYB genes, to determine which of the genes belonged to which R2R3-MYB subfamily, the R2R3-MYB genes from Arabidopsis thaliana were used to conduct phylogenetic analyses (Amino acid sequence obtained from Gene Bank, accession number is listed in Supplementary Table S2), and MEGA7.0 software (Institute of Molecular Evolutionary Genetics, PA, USA) [52] was used to perform sequence aligning and construct a circular phylogenetic tree according to neighbor joining method with 1000 interactions.

4.9. Structure Genes and Tfs Co-Expression Network

The String online database (https://string-db.org/) was used to search for the interaction relationships between structure genes and transcription factors. Then we used Cytoscape software (National Institute of General Medical Sciences, Bethesda, MD, USA) to construct a co-expression network to identify their main impact factors [53].

5. Conclusions

The ratio of 6:3:1 distribution patterns of P. limprichtii within the population in Huanglong District seems to reveal that pollinators, herbivores, and survival adaptability could promote the development of such a reproductive strategy in the population to maintain the maximum population density. Function and expression patterns indicated that PlFLS probably is a crucial gene in the formation of white and pigmented flowers. PlANS and PlUFGT were found to determine the color intensity in pigmented flowers. In addition, a putative MBW complex, PlMYB10/PlbHLH20/PlWD40-1 or PlMYB10/PlbHLH26/PlWD40-1 may serve as a repressor of regulated PlFLS expression that is responsible for variation in the color intensity of P. limprichtii. Our results provide valuable molecular information on floral color variations in Pleione, and also provide inspiration to further explore the relationship between color polymorphism and species evolution and to study its contribution to color evolution.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/1/247/s1.

Author Contributions

This study was conceived by S.W. and Z.L. Y.Z. performed the transcriptome and gene expression analysis and wrote the manuscript. T.Z. and Z.D. helped with the sampling and field work. X.D., W.L., and M.C. helped to perform the experiment. C.L. and X.W. revised the manuscript. W.-C.T. and J.Z. provided writing guidance. All authors contributed to revising the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fujian Natural Science Foundation Project (grant no. 2019J01410), the Fujian Agriculture and Forestry University 2015 Outstanding Youth Fund Project (grant no. xjq201620), and the Fujian Agriculture and Forestry University Science and Technology Innovation Special Fund Project (grant no. KFA17331A).

Acknowledgments

The authors thank Bai-Jun Li, Jie-Yu Wang, Yan-Qiong Chen, Guo-Qiang Zhang, Si-jin Zeng, Zhi-yao Ren, and Xueyan Yuan who helped with this research. We also thank Qi-Xuan Peng and other staff of Huanglong National Scenic Reserve for assisting sampling.

Conflicts of Interest

The authors declare no conflict of interest. Huanglong National Scenic Reserve help to sample, and OmicShare company help to analysis raw data.

Abbreviations

UPLCUltra performance liquid chromatography
ABPAnthocyanin biosynthesis pathway
CHSChalcone synthase
CHIChalcone isomerase
F3HFlavanone 3-hydroxylase
F3′HFlavonoid 3’-hydroxylase
F3′5′HFlavonoid 3’5′-hydroxyla
FLSFlavonol synthase
DFRDihydroflavonol 4-reduct
ANSAnthocyanidin synthase
ANRAnthocyanidin reductase
UFGTUDP flavonoid glucosyl transferase
RHSRoyal horticulture society
RHSCC The royal horticulture society color chart
NRNon-redundant protein database
Swiss-Prot Swiss-Protein protein database
KEGG Kyoto encyclopedia of genes and genomes
KOGEukaryotic orthologous groups
GOGene ontology
RPKMreads per kb per million reads
FDRfalse discovery rate
qRT-PCR Real-time reverse transcription-PCR

References

  1. Cooley, A.M.; Carvallo, G.; Willis, J.H. Is floral diversification associated with pollinator divergence? flower shape, flower colour and pollinator preference in Chilean Mimulus. Ann. Bot. 2008, 101, 641–650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. 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]
  3. Vaidya, P.; Ansley, M.; Emily, M.; Michaela, K.; Lauren, C.; Lee, C.R.; Bingham, R.A.; Anderson, J.T. Ecological causes and consequences of flower color polymorphism in a self-pollinating plant (Boechera Stricta). New Phytol. 2018, 218, 380–392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Streinzer, M.; Roth, N.; Paulus, H.F.; Spaethe, J. Color preference and spatial distribution of glaphyrid beetles suggest a key role in the maintenance of the color polymorphism in the peacock anemone (Anemone pavonina, Ranunculaceae) in Northern Greece. J. Comp. Physiol. 2019, 205, 735–743. [Google Scholar] [CrossRef] [PubMed]
  5. Kellenberger, R.T.; Byers, K.J.R.P.; De Brito Francisco, R.M.; Staedler, Y.M.; LaFountain, A.M.; Schönenberger, J.; Schiestl, F.P.; Schlüter, P.M. Emergence of a floral colour polymorphism by pollinator-mediated overdominance. Nat. Commun. 2019, 10, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sobel, J.M.; Matthew, A.S. Flower color as a model system for studies of plant evo-devo. Front. Plant Sci. 2013, 4, 321. [Google Scholar] [CrossRef] [PubMed]
  7. Rausher, M.D. Evolutionary transitions in floral color. Int. J. Plant Sci. 2008, 169, 7–21. [Google Scholar] [CrossRef] [Green Version]
  8. Lou, Q.; Liu, Y.L.; Qi, Y.Y.; Jiao, S.Z.; Tian, F.F.; Jiang, L.; Wang, Y.J. Transcriptome sequencing and metabolite analysis reveals the role of delphinidin metabolism in flower color in grape hyacinth. J. Exp. Bot. 2014, 65, 3157–3164. [Google Scholar] [CrossRef] [Green Version]
  9. Butler, T.; Dick, C.; Carlson, M.L.; Whittall, J.B. Transcriptome analysis of a petal anthocyanin polymorphism in the Arctic Mustard, Parrya Nudicaulis. PLoS ONE 2014, 9, e101338. [Google Scholar] [CrossRef] [Green Version]
  10. Casimiro-Soriguer, I.; Eduardo, N.; Buide, M.L.; Valle, J.C.; Whittall, J.B. Transcriptome and biochemical analysis of a flower color polymorphism in Silene Littorea (Caryophyllaceae). Front. Plant Sci. 2016, 7, 204. [Google Scholar] [CrossRef] [Green Version]
  11. Gao, L.X.; Yang, H.X.; Liu, H.F.; Yang, J.; Hu, Y.H. Extensive transcriptome changes underlying the flower color intensity variation in Paeonia Ostii. Front. Plant Sci. 2016, 6, 1205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, H.; Lucie, C.C.; Bessière, J.M.; Guillaume, C.; Bertrand, S.; Eric, I. Flower color polymorphism in Iris Lutescens (Iridaceae): Biochemical analyses in light of plant–insect interactions. Phytochemistry 2013, 94, 123–134. [Google Scholar] [CrossRef] [PubMed]
  13. Wu, C.A.; Matthew, A.S.; Laura, I.N.; Kaitlyn, A.C. The genetic basis of a rare flower color polymorphism in Mimulus Lewisii provides insight into the repeatability of evolution. PLoS ONE 2013, 8, e81173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Grotewold, E. The genetic and biochemistry of floral pigments. Annu. Rev. Plant Biol. 2006, 57, 761–780. [Google Scholar] [CrossRef] [PubMed]
  15. Holton, T.A.; Edwina, C.C. Genetics and biochemistry of anthocyanin biosynthesis. Plant Cell 1995, 7, 1071–1083. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, D.; Tao, J. Recent advances on the development and regulation of flower color in ornamental plants. Front. Plant Sci. 2015, 6, 261. [Google Scholar] [CrossRef] [Green Version]
  17. Tanaka, Y.; Akemi, O. Seeing is believing: Engineering anthocyanin and carotenoid biosynthetic pathways. Curr. Opin. Biotechnol. 2008, 19, 190–197. [Google Scholar] [CrossRef]
  18. Xie, D.Y.; Shashi, B.S.; Elane, W.; Wang, Z.Y.; Richard, A.D. Metabolic engineering of proanthocyanidins through co-expression of anthocyanidin reductase and the PAP1 MYB transcription factor. Plant J. 2006, 45, 895–907. [Google Scholar] [CrossRef] [Green Version]
  19. Twyford, A.D.; Aaron, M.C.; Pratibha, C.; Ramesh, R.; Jannice, F. Loss of color pigmentation is maintained at high frequency in a monkey flower population. Am. Nat. 2018, 191, 135–145. [Google Scholar] [CrossRef]
  20. Hichri, I.; Barrieu, F.; Bogs, J.; Kappel, C.; Delrot, S.; Lauvergeat, V. Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway. J. Exp. Bot. 2011, 62, 2465–2483. [Google Scholar] [CrossRef] [Green Version]
  21. Albert, N.W.; Kevin, M.D.; David, H.L.; Zhang, H.B.; Mirco, M.; Cyril, B.; Boase, M.R.; Hanh, N.; Jameson, P.E.; Schwinn, K.E. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in Eudicots. Plant Cell 2014, 26, 962–980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wessinger, C.A.; Rausher, M.D. Lessons from flower color evolution on targets of selection. J. Exp. Bot. 2012, 63, 695–709. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chen, X.; Cribb, P.J.; Gale, S.W. “Pleione”. In the Flora of China 25 (Orchidaceae); Wu, Z.Y., Raven, P.H., Hong, D.Y., Eds.; Beijing Science Press and Missouri Botanical Garden Press: Beijing, China, 2009; pp. 325–333. [Google Scholar]
  24. Cribb, P.; Butterfield, I. The genus Pleione, 2nd ed.; Royal Botanic Gardens/Kew Press: London, UK, 1999; p. 27. [Google Scholar]
  25. Govaerts, R.; Campacci, M.A.; Baptista, D.H.; Cribb, J.; George, A.; Kreutz, K. Data from: World Checklist of Orchidaceae; The Board of Trustees of the Royal Botanic Gardens/Kew: London, UK, 2016; Available online: http://apps.kew.org/wcsp (accessed on 20 August 2019).
  26. Gravendeel, B.; Eurlings, M.C.M.; Berg, C.V.D.; Cribb, P.J. Phylogeny of Pleione (Orchidaceae) and parentage analysis of its wild hybrids based on plastid and nuclear ribosomal its sequences and morphological data. Syst. Bot. 2004, 29, 50–63. [Google Scholar] [CrossRef]
  27. Stracke, R.; Werber, M.; Weisshaar, B. The R2R3-MYB gene family in Arabidopsis Thaliana. Curr. Opin. Biotechnol. 2001, 4, 447–456. [Google Scholar] [CrossRef]
  28. Dubos, C.; Ralf, S.; Grotewold, E.; Weisshaar, B.; Martin, C.; Loïc, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef] [PubMed]
  29. Arista, M.; Talavera, M.; Berjano, R.; Ortiz, P.L. Abiotic factors may explain the geographical distribution of flower colour morphs and the maintenance of colour polymorphism in the scarlet pimpernel. J. Ecol. 2013, 101, 1613–1622. [Google Scholar] [CrossRef]
  30. Imbert, E.; Wang, H.; Conchou, L.; Vincent, H.; Talavera, M.; Schatz, B. Positive effect of the yellow morph on female reproductive success in the flower colour polymorphic Iris lutescens (Iridaceae), a deceptive species. J. Evol. Biol. 2014, 27, 1965–1974. [Google Scholar] [CrossRef]
  31. Chalker-Scott, L. Environmental significance of anthocyanins in plant stress responses. Photochem. Photobiol. 1999, 70, 1–9. [Google Scholar] [CrossRef]
  32. Muchhala, N.; Sönke, J.; Smith, S.D. Competition for hummingbird pollination shapes flower color variation in Andean solanaceae: Competition for pollination shapes flower color variation. Evolution 2014, 68, 2275–2286. [Google Scholar] [CrossRef]
  33. Parachnowitsch, A.L.; Christina, M.C. Predispersal seed herbivores, not pollinators, exert selection on floral traits via female fitness. Ecology 2008, 89, 1802–1810. [Google Scholar] [CrossRef]
  34. Jin, X.H.; Huang, H.; Wang, L.; Sun, Y.; Dai, S.L. Transcriptomics and metabolite analysis reveals the molecular mechanism of anthocyanin biosynthesis branch pathway in different senecio cruentus cultivars. Front. Plant Sci. 2016, 7, 1307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Le Maitre, N.C.; Pirie, M.D.; Bellstedt, D.U. Floral color, anthocyanin synthesis gene expression and control in cape erica species. Front. Plant Sci. 2019, 10, 1565. [Google Scholar] [CrossRef] [PubMed]
  36. Li, L.; Zhai, Y.; Luo, X.; Zhang, Y.; Shi, Q.Q. Comparative transcriptome analyses reveal genes related to pigmentation in the petals of red and white Primula vulgaris cultivars. Physiol. Mol. Biol. Plants 2019, 25, 1029–1041. [Google Scholar] [CrossRef] [PubMed]
  37. Yuan, Y.W.; Alexandra, B.R.; Janelle, M.S.; Lauren, E.S.; Harvey, D.B. Competition between anthocyanin and flavonol biosynthesis produces spatial pattern variation of floral pigments between Mimulus species. Proc. Natl. Acad. Sci. USA 2016, 113, 2448–2453. [Google Scholar] [CrossRef] [Green Version]
  38. Davies, K.M.; Schwinn, K.E.; Deroles, S.C.; Manson, D.G.; Lewis, D.H.; Bloor, S.J.; Bradley, J.M. Enhancing anthocyanin production by altering competition for substrate between flavonol synthase and dihydroflavonol 4-Reductase. Euphytica 2003, 131, 259–268. [Google Scholar] [CrossRef]
  39. Sangkyu, P.; Kim, D.H.; Lee, J.Y.; Ha, S.H.; Lim, S.H. Comparative analysis of two flavonol synthases from different-colored onions provides insight into flavonoid biosynthesis. J. Agric. Food Chem. 2017, 65, 5287–5298. [Google Scholar]
  40. Shi, S.G.; Li, S.J.; Kang, Y.X.; Liu, J.J. Molecular characterization and expression analyses of an anthocyanin synthase gene from magnolia sprengeri Pamp. Int. J. Appl. Biotechnol. Biochem. 2015, 175, 477–488. [Google Scholar] [CrossRef]
  41. Wu, X.X.; Gong, Q.H.; Ni, X.P.; Zhou, Y.; Gao, Z.H. UFGT: The key enzyme associated with the petals variegation in Japanese Apricot. Front. Plant Sci. 2017, 8, 108. [Google Scholar] [CrossRef] [Green Version]
  42. Feller, A.; Machemer, K.; Edward, L.B.; Grotewold, E. Evolutionary and comparative analysis of MYB and BHLH plant transcription factors. Plant J. 2011, 66, 94–116. [Google Scholar] [CrossRef]
  43. Heim, M.A. The basic Helix-Loop-Helix transcription factor family in plants: A genome-wide study of protein structure and functional diversity. Mol. Biol. Evol. 2003, 20, 735–747. [Google Scholar] [CrossRef] [Green Version]
  44. Baudry, A.; Heim, M.A.; Dubreucq, B.; Caboche, M.; Weisshaar, B.; Lepiniec, L. TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis Thaliana. Plant J. 2004, 39, 366–380. [Google Scholar] [CrossRef] [PubMed]
  45. Zimmermann, I.; Heim, M.; Weisshaar, B.; Uhrig, J. Comprehensive identification of Arabidopsis thaliana MYB transcription factors interacting with R/B-like bHLH proteins: Systematic analysis of MYB/BHLH-Interaction. Plant J. 2004, 40, 22–34. [Google Scholar] [CrossRef] [PubMed]
  46. Akagi, T.; Ikegami, A.; Tsujimoto, T.; Kobayashi, S.; Sato, A.; Kono, A.; Yonemori, K. DkMyb4 is a myb transcription factor involved in proanthocyanidin biosynthesis in persimmon fruit. Plant Physiol. 2009, 151, 2028–2045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Matsui, K.; Yoshimi, U.; Masaru, O.T. AtMYBL2, a protein with a single MYB domain, acts as a negative regulator of anthocyanin biosynthesis in Arabidops. Plant J. 2008, 55, 954–967. [Google Scholar] [CrossRef] [PubMed]
  48. Lang, X.A.; Li, N.; Li, L.F.; Zhang, S.Z. Integrated metabolome and transcriptome analysis uncovers the role of anthocyanin metabolism in Michelia maudiae. Int. J. Genom. 2019, 2019, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Gu, Z.; Zhu, J.; Hao, Q.; Yuan, Y.; Duan, Y.; Men, S.; Wang, Q.; Hou, Q.; Liu, Z.; Shu, Q.; et al. A novel R2R3-MYB transcription factor contributes to petal blotch formation by regulating organ-Sspecific expression of PsCHS in tree peony (Paeonia Suffruticosa). Plant Cell Physiol. 2018, 60, 599–611. [Google Scholar] [CrossRef] [PubMed]
  50. Chen, W.; Gong, L.; Guo, Z.L.; Wang, W.S.; Zhang, H.Y.; Liu, X.Q.; Yu, S.B.; Xiong, L.Z.; Luo, J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: Application in the study of rice metabolomics. Mol. Plant 2013, 6, 1769–1780. [Google Scholar] [CrossRef] [Green Version]
  51. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  52. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  53. Shannon, P. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
Ijms 21 00247 g001 550
Figure 1. Principal component analysis of petal color of P. limprichtii. Distribution based on bivariate values of a* (redness and greenness) and b* (yellowness and blueness). a* means difference in red and green, b* means difference in yellow and blue.
Figure 1. Principal component analysis of petal color of P. limprichtii. Distribution based on bivariate values of a* (redness and greenness) and b* (yellowness and blueness). a* means difference in red and green, b* means difference in yellow and blue.
Ijms 21 00247 g001
Ijms 21 00247 g002 550
Figure 2. Venn diagram of the number of unigenes annotated by BLASTx with four protein databases.
Figure 2. Venn diagram of the number of unigenes annotated by BLASTx with four protein databases.
Ijms 21 00247 g002
Ijms 21 00247 g003 550
Figure 3. E-value distribution of top BLASTx hits against four protein databases for each unigene. (a) distribution E-values in non-redundant protein database; (b) distribution E-values in Swiss-protein database; (c) distribution E-values in the Kyoto Encyclopedia of Genes and Genomes database; (d) distribution E-values in the Eukaryotic Orthologous Group database.
Figure 3. E-value distribution of top BLASTx hits against four protein databases for each unigene. (a) distribution E-values in non-redundant protein database; (b) distribution E-values in Swiss-protein database; (c) distribution E-values in the Kyoto Encyclopedia of Genes and Genomes database; (d) distribution E-values in the Eukaryotic Orthologous Group database.
Ijms 21 00247 g003
Ijms 21 00247 g004 550
Figure 4. Function classification of the gene ontology of all unigenes based on Nr annotation.
Figure 4. Function classification of the gene ontology of all unigenes based on Nr annotation.
Ijms 21 00247 g004
Ijms 21 00247 g005 550
Figure 5. Expression pattern analysis base on RNA sequencing of anthocyanin biosynthesis pathway-related genes and transcription factors in P. limprichtii. (a) Anthocyanin biosynthesis pathway-related unigenes; (b) bHLH unigenes; (c) R2R3-MYB unigenes; (d) WD40 unigenes. W-pe (petal of white flower), P-petal (petal of pink flower), R-petal (petal of rose -purple flower).
Figure 5. Expression pattern analysis base on RNA sequencing of anthocyanin biosynthesis pathway-related genes and transcription factors in P. limprichtii. (a) Anthocyanin biosynthesis pathway-related unigenes; (b) bHLH unigenes; (c) R2R3-MYB unigenes; (d) WD40 unigenes. W-pe (petal of white flower), P-petal (petal of pink flower), R-petal (petal of rose -purple flower).
Ijms 21 00247 g005
Ijms 21 00247 g006 550
Figure 6. Real time quantitative reverse transcription-PCR of several genes in P. limprichtii. Each value is shown as average ± standard deviation from three biological replicate sampling.
Figure 6. Real time quantitative reverse transcription-PCR of several genes in P. limprichtii. Each value is shown as average ± standard deviation from three biological replicate sampling.
Ijms 21 00247 g006
Ijms 21 00247 g007 550
Figure 7. Phylogenetic analysis of R2R3-MYB DNA binding domains for P. limprichtii and Arabidopsis thaliana. (a) Circular phylogenetic tree; (b) amplification of S4, S6, and S7 branches. The R2R3 domains of the 14 MYBs identified in P. limprichtii petal transcriptome were aligned and analyzed using neighbor-joining phylogenetic methods.
Figure 7. Phylogenetic analysis of R2R3-MYB DNA binding domains for P. limprichtii and Arabidopsis thaliana. (a) Circular phylogenetic tree; (b) amplification of S4, S6, and S7 branches. The R2R3 domains of the 14 MYBs identified in P. limprichtii petal transcriptome were aligned and analyzed using neighbor-joining phylogenetic methods.
Ijms 21 00247 g007
Ijms 21 00247 g008 550
Figure 8. A co-expression network of anthocyanin biosynthesis pathway-related genes and transcription factors involved in pigmentation. Rose-purple circle represents structural genes; dark-gray circle represents R2R3-MYB; light-gray circle represents bHLH; pink circle represents WD40.
Figure 8. A co-expression network of anthocyanin biosynthesis pathway-related genes and transcription factors involved in pigmentation. Rose-purple circle represents structural genes; dark-gray circle represents R2R3-MYB; light-gray circle represents bHLH; pink circle represents WD40.
Ijms 21 00247 g008
Ijms 21 00247 g009 550
Figure 9. Tentative pathways for P. limprichtii color variations. The colored bar is the value of log2 (RPKM + 1), represented using the depth of color, with purple representing the up-regulated expression genes and pink representing the down-regulated expression genes. RPKM means the reads per kb per million reads mapped. CHS, chalcone synthase; CHI, chalcone-flavanone isomerase; F3H, flavanone-3-hydroxylase; F3’H, flavonoid 3’-hydroxylase; F3’5’H, flavonoid 3’5’-hydroxyla; DFR, dihydroflavonols 4-reductase; ANS, anthocyanidin synthase; UFGT, UDP flavonoid glucosyl transferase; FLS, flavonol synthase; The three gene complex consist of a MYB, bHLH and WD40 in most angiosperm.
Figure 9. Tentative pathways for P. limprichtii color variations. The colored bar is the value of log2 (RPKM + 1), represented using the depth of color, with purple representing the up-regulated expression genes and pink representing the down-regulated expression genes. RPKM means the reads per kb per million reads mapped. CHS, chalcone synthase; CHI, chalcone-flavanone isomerase; F3H, flavanone-3-hydroxylase; F3’H, flavonoid 3’-hydroxylase; F3’5’H, flavonoid 3’5’-hydroxyla; DFR, dihydroflavonols 4-reductase; ANS, anthocyanidin synthase; UFGT, UDP flavonoid glucosyl transferase; FLS, flavonol synthase; The three gene complex consist of a MYB, bHLH and WD40 in most angiosperm.
Ijms 21 00247 g009
Ijms 21 00247 g010 550
Figure 10. Flower polymorphism of Pleione limprichtii in Huanglong population. (a) rose-purple flower; (b) pink flower; (c) white flower; (d) one of the polymorphic populations on a rock.
Figure 10. Flower polymorphism of Pleione limprichtii in Huanglong population. (a) rose-purple flower; (b) pink flower; (c) white flower; (d) one of the polymorphic populations on a rock.
Ijms 21 00247 g010
Table
Table 1. Number of P. limprichtii individuals of each flower color polymorph in three selected populations at the opening stage.
Table 1. Number of P. limprichtii individuals of each flower color polymorph in three selected populations at the opening stage.
PopulationRose-PurplePinkWhite
population1409 (58%)227 (32%)70 (10%)
population231 (65%)15 (31%)2 (4%)
population392 (60%)54 (35%)8 (5%)
The percentage indicates the proportion of individuals with a given color account for the total number of individuals in the population.
Table
Table 2. Color parameters of three distinct colors of P. limprichtii petal.
Table 2. Color parameters of three distinct colors of P. limprichtii petal.
Flower ColorRHSL*C*b*/a*
purple redPurple-group N78B67.92 ± 10.2749.64 ± 9.41−1.68 ± 0.61
76.66 ± 4.5152.88 ± 4.18−2.35 ± 0.23
67.35 ± 10.3457.13 ± 6.24−1.68 ± 0.61
81.31 ± 8.7251.89 ± 7.42−1.91 ± 0.44
79.68 ± 15.5546.28 ± 11.76−1.85 ± 0.79
PinkPurple-group 75A61.08 ± 5.2926.71 ± 1.82−1.20 ± 0.36
66.05 ± 6.6335.45 ± 9.60−1.77 ± 0.71
66.35 ± 2.7331.72 ± 3.14−1.57 ± 0.29
58.80 ± 5.2029.66 ± 6.05−2.20 ± 1.07
62.80 ± 5.7825.43 ± 4.99−2.56 ± 0.89
WhitePurple-group 76D93.08 ± 5.5713.31 ± 9.67−4.04 ± 2.06
88.12 ± 8.1318.53 ± 7.31−13.60 ± 5.70
95.85 ± 4.4822.41 ± 6.77−7.97 ± 3.88
87.52 ± 11.4525.27 ± 10.21−6.11 ± 2.61
87.07 ± 7.9721.36 ± 6.12−9.28 ± 4.73
RHS, Royal horticulture society color chart evaluation index; L*, lightness; a*, b*, chromatic components; C*, chroma (brightness).
Table
Table 3. Anthocyanin content of petal in P. limprichtii.
Table 3. Anthocyanin content of petal in P. limprichtii.
IndexIon ModeMolecular Weight Substance
pme3609positive287.24Cyanidin
pme0442positive303.24Delphinidin
pmb0541positive697.1Cyanidin 3-O-glucosyl-malonylglucoside
pmb0542positive535.1Cyanidin 3-O-malonylhexoside
Table
Table 4. Tentative anthocyanin biosynthesis pathway-related genes in P. limprichtii.
Table 4. Tentative anthocyanin biosynthesis pathway-related genes in P. limprichtii.
PathwayGeneEncoding EnzymeNumber
Flavonoid biosynthesisCHSChalcone synthase5
CHIChalcone isomerase2
F3HFlavanone 3-dioxygenase1
F3′HFlavonoid 3’-hydroxylase3
F3′5′HFlavonoid 3’,5’-hydroxylase1
DFRDihydroflavonol 4-reductase3
ANSAnthocyanidin synthase2
ANRAnthocyanidin reductase1
FLSFlavonol synthase1
Anthocyanin biosynthesisUFGTAnthocyanidin 3-O-glucosyltransferase2

Share and Cite

MDPI and ACS Style

Zhang, Y.; Zhou, T.; Dai, Z.; Dai, X.; Li, W.; Cao, M.; Li, C.; Tsai, W.-C.; Wu, X.; Zhai, J.; et al. Comparative Transcriptomics Provides Insight into Floral Color Polymorphism in a Pleione limprichtii Orchid Population. Int. J. Mol. Sci. 2020, 21, 247. https://doi.org/10.3390/ijms21010247

AMA Style

Zhang Y, Zhou T, Dai Z, Dai X, Li W, Cao M, Li C, Tsai W-C, Wu X, Zhai J, et al. Comparative Transcriptomics Provides Insight into Floral Color Polymorphism in a Pleione limprichtii Orchid Population. International Journal of Molecular Sciences. 2020; 21(1):247. https://doi.org/10.3390/ijms21010247

Chicago/Turabian Style

Zhang, Yiyi, Tinghong Zhou, Zhongwu Dai, Xiaoyu Dai, Wei Li, Mengxia Cao, Chengru Li, Wen-Chieh Tsai, Xiaoqian Wu, Junwen Zhai, and et al. 2020. "Comparative Transcriptomics Provides Insight into Floral Color Polymorphism in a Pleione limprichtii Orchid Population" International Journal of Molecular Sciences 21, no. 1: 247. https://doi.org/10.3390/ijms21010247

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