Biochemical and Comparative Transcriptomic Analyses Identify Candidate Genes Related to Variegation Formation in Paeonia rockii

Integrative Omics Analysis Reveals Mechanisms of Anthocyanin Biosynthesis in Djulis Spikes

Djulis (Chenopodium formosanum Koidz.), a member of the Amaranthaceae family plant, is noted for its vibrant appearance and significant ornamental value. However, the mechanisms underlying color variation in its spikes remain unexplored. This research initially detected the anthocyanin content at different developmental stages of the spike and subsequently utilized an integrative approach, combining targeted metabolomics, transcriptomics, and untargeted metabolomics analyses, to elucidate the mechanisms of anthocyanin biosynthesis in the spikes of djulis. The results of the combined multi-omics analysis showed that the metabolites associated with anthocyanin synthesis were mainly enriched in the flavonoid biosynthesis pathway (ko00941) and the anthocyanin biosynthesis pathway (ko00942). With the maturation of djulis spikes, a total of 28 differentially expressed genes and 17 differentially expressed metabolites were screened during the transition of spike color from green (G) to red (R) or orange (O). Twenty differentially expressed genes were selected for qRT-PCR validation, and the results are consistent with transcriptome sequencing. The upregulation of seven genes, including chalcone synthase (CfCHS3_1, CfCHS3_2, CfCHS3_3), flavanone 3-hydroxylase (CfF3H_3), flavonoid 3'5'-hydroxylase (CfCYP75A6_1), dihydroflavonol reductase (CfDFRA), and glucosyltransferase (Cf3GGT), promotes the formation and accumulation of delphinidin 3-sambubioside and peonidin 3-galactoside. The research results also showed that anthocyanins and betalains can coexist in the spike of djulis, and the reason for the change in spike color during development may be the result of the combined action of the two pigments. A possible regulatory pathway for anthocyanin biosynthesis during the spike maturation was constructed based on the analysis results. The results provide a reference and theoretical basis for further studying the molecular mechanism of anthocyanin regulation of color changes in Amaranthaceae plants.

Functional identification of anthocyanin glucosyltransferase genes: a Ps3GT catalyzes pelargonidin to pelargonidin 3-O-glucoside painting the vivid red flower color of Paeonia

Glycosylation from an anthocyanidin 3-O-glucosyltransferase Ps3GT (PsUGT78A27) facilitates the accumulation of pelargonidin 3-O-glucoside, which defines the vivid red flower color and occurs only in specific peony tree cultivars. Although tree peony cultivars of Chinese and Japanese both originated from China, vivid red color is only found in flowers of Japanese cultivars but not of Chinese cultivar groups. In this study, a Japanese tree peony cultivar 'Taiyoh' with vivid red petals and a Chinese tree peony cultivar 'Hu Hong' with reddish pink petals were chosen as the experimental materials. Flavonoids profiling indicated that pelargonidin 3-O-glucoside (Pg3G) detected only in Japanese cultivar contributed to vivid red color of tree peony petals, while pelargonidin 3,5-di-O-glucoside (Pg3G5G) found in both of Japanese and Chinese cultivars was responsible for pink flower color. Through the integration of full-length transcriptome sequencing and in vitro enzymatic activity analysis, two anthocyanin glucosyltransferase genes PsUGT78A27 and PsUGT75L45 were isolated from the petals of tree peony, and their encoding products exhibited enzymatic activities of pelargonidin 3-O-glucosyltransferase and anthocyanin 5-O-glucosyltransferase, respectively. Further quantitative real-time PCR revealed that PsUGT78A27 displayed high expression in petals of both cultivars and PsUGT75L45 was expressed at high levels in cultivar 'Hu Hong' only. Using a gene gun technique, the GFP fusion proteins of PsUGT78A27 and PsUGT75L45 were visualized to be cytoplasmic and nuclear localization in the epidermal cells of tree peony petals, and the glucosylation function of PsUGT78A27 and PsUGT75L45 to alter petal color of tree peony and herbaceous peony had been directly validated in vivo. These results demonstrated that PsUGT78A27 and PsUGT75L45 are key players for the presence or absence of vivid red flower color in tree peony cultivars. Our findings further elucidated the chemical and molecular mechanism of petal pigmentation of Paeonia and could help breed the Paeonia cultivars possessing novel flower colors.

PhUGT78A22, a novel glycosyltransferase in Paeonia 'He Xie', can catalyze the transfer of glucose to glucosylated anthocyanins during petal blotch formation

BACKGROUND

Flower color patterns play an important role in the evolution and subsequent diversification of flowers by attracting animal pollinators. This interaction can drive the diversity observed in angiosperms today in many plant families such as Liliaceae, Paeoniaceae, and Orchidaceae, and increased their ornamental values. However, the molecular mechanism underlying the differential distribution of anthocyanins within petals remains unclear in Paeonia.

RESULTS

In this study, we used an intersectional hybrid between the section Moutan and Paeonia, hereafter named Paeonia 'He Xie', which has purple flowers with dark purple blotches. After Ultra-high performance liquid chromatography-diode array detector (UPLC-DAD) analysis of blotched and non-blotched parts of petals, we found the anthocyanin content in the blotched part was always higher than that in the non-blotched part. Four kinds of anthocyanins, namely cyanidin-3-O-glucoside (Cy3G), cyanidin-3,5-O-glucoside (Cy3G5G), peonidin-3-O-glucoside (Pn3G), and peonidin-3,5-O-glucoside (Pn3G5G) were detected in the blotched parts, while only Cy3G5G and Pn3G5G were detected in the non-blotched parts. This suggests that glucosyltransferases may play a vital role in the four kinds of glucosylated anthocyanins in the blotched parts. Moreover, 2433 differentially expressed genes (DEGs) were obtained from transcriptome analysis of blotched and non-blotched parts, and a key UDP-glycosyltransferase named PhUGT78A22 was identified, which could use Cy3G and Pn3G as substrates to produce Cy3G5G and Pn3G5G, respectively, in vitro. Furthermore, silencing of PhUGT78A22 reduced the content of anthocyanidin 3,5-O-diglucoside in P. 'He Xie'.

CONCLUSIONS

A UDP-glycosyltransferase, PhUGT78A22, was identified in P. 'He Xie', and the molecular mechanism underlying differential distribution of anthocyanins within petals was elucidated. This study provides new insights on the biosynthesis of different kinds of anthocyanins within colorful petals, and helps to explain petal blotch formation, which will facilitate the cultivar breeding with respect to increasing ornamental value. Additionally, it provides a reference for understanding the molecular mechanisms responsible for precise regulation of anthocyanin biosynthesis and distribution patterns.

PrMYB5 activates anthocyanin biosynthetic PrDFR to promote the distinct pigmentation pattern in the petal of Paeonia rockii

Paeonia rockii is well-known for its distinctive large dark-purple spot at the white petal base and has been considered to be the main genetic source of spotted tree peony cultivars. In this study, the petal base and petal background of Paeonia ostii (pure white petals without any spot), P. rockii, and other three tree peony cultivars were sampled at four blooming stages from the small bell-like bud stage to the initial blooming stage. There is a distinct difference between the pigmentation processes of spots and petal backgrounds; the spot pigmentation was about 10 days earlier than the petal background. Moreover, the cyanin and peonidin type anthocyanin accumulation at the petal base mainly contributed to the petal spot formation. Then, we identified a C1 subgroup R2R3-MYB transcription factor, PrMYB5, predominantly transcribing at the petal base. This is extremely consistent with PrDFR and PrANS expression, the contents of anthocyanins, and spot formation. Furthermore, PrMYB5 could bind to and activate the promoter of PrDFR in yeast one-hybrid and dual-luciferase assays, which was further verified in overexpression of PrMYB5 in tobacco and PrMYB5-silenced petals of P. rockii by comparing the color change, anthocyanin contents, and gene expression. In summary, these results shed light on the mechanism of petal spot formation in P. rockii and speed up the molecular breeding process of tree peony cultivars with novel spot pigmentation patterns.

Comparative transcriptomic analysis of genes involved in stem lignin biosynthesis in woody and herbaceous Paeonia species

Paeonia is recognized globally due to its ornamental value. However, the mechanisms behind the formation of distinct levels of lignification in Paeonia stems remain largely unknown. In this study, we selected three representative Paeonia species, namely P. ostii (shrub), P. lactiflora (herb), and P. × 'Hexie' (semi-shrub), to evaluate and contrast their respective anatomical structure, phytochemical composition and transcriptomic profile. Our results showed that the degree of lignin deposition on the cell wall, along with the total amount of lignin and its monomers (especially G-lignin) were higher in P. ostii stems compared to the other two species at almost all development stages except 80 days after flowering. Furthermore, we estimated a total number of unigenes of 60,238 in P. ostii, 43,563 in P. × 'Hexie', and 40,212 in P. lactiflora from stem transcriptome. We then built a co-expression network of 25 transcription factors and 21 enzyme genes involved in lignin biosynthesis and identified nine key candidate genes. The expression patterns of these genes were positively correlated with the transcription levels of PAL, C4H, 4CL2, CCR, and COMT, as well as lignin content. Moreover, the highest relative expression levels of CCR, 4CL2, and C4H were found in P. ostii. This study provides an explanation for the observed differences in lignification between woody and herbaceous Paeonia stems, and constitutes a novel reference for molecular studies of stem-specific lignification process and lignin biosynthesis that can impact the ornamental industry.

Characterization of Metabolites and Transcripts Involved in Flower Pigmentation in Primula vulgaris

Primula vulgaris exhibits a wide range of flower colors and is a valuable ornamental plant. The combination of flavonols/anthocyanins and carotenoids provides various colorations ranging from yellow to violet-blue. However, the complex metabolic networks and molecular mechanisms underlying the different flower colors of P. vulgaris remain unclear. Based on comprehensive analysis of morphological anatomy, metabolites, and gene expression in different-colored flowers of P. vulgaris, the mechanisms relating color-determining compounds to gene expression profiles were revealed. In the case of P. vulgaris flower color, hirsutin, rosinin, petunidin-, and cyanidin-type anthocyanins and the copigment herbacetin contributed to the blue coloration, whereas peonidin-, cyandin-, and delphinidin-type anthocyanins showed high accumulation levels in pink flowers. The color formation of blue and pink were mainly via the regulation of F3'5'H (c53168), AOMT (c47583, c44905), and 3GT (c50034). Yellow coloration was mainly due to gossypetin and carotenoid, which were regulated by F3H (c43100), F3 1 (c53714), 3GT (c53907) as well as many carotenoid biosynthetic pathway-related genes. Co-expression network and transient expression analysis suggested a potential direct link between flavonoid and carotenoid biosynthetic pathways through MYB transcription factor regulation. This work reveals that transcription changes influence physiological characteristics, and biochemistry characteristics, and subsequently results in flower coloration in P. vulgaris.

PsbHLH1, a novel transcription factor involved in regulating anthocyanin biosynthesis in tree peony (Paeonia suffruticosa)

Flower color is one of the most important features of ornamental plants. Anthocyanin composition and concentration are usually closely related to flower color formation. The biosynthesis of anthocyanin is regulated by a series of structural genes and regulatory genes. The basic helix-loop-helix proteins (bHLHs) are considered as one of the key transcription factors known as the regulators of anthocyanin biosynthesis. However, the bHLH transcription factor family of tree peony (Paeonia suffruticosa) has not been systematically studied in previous studies, especially for the regulation of petal pigmentation. The aim of this study was to identify bHLH genes and unravel their underlying molecular mechanism involved in the regulation of anthocyanin biosynthesis in tree peony. Based on transcriptome profiling analysis, we identified three bHLHs candidate anthocyanin regulators, PsbHLH1, PsbHLH2, and PsbHLH3. PsbHLH1-3 were phylogenetically clustered in the IIIf bHLH subgroup, which is involved in anthocyanin biosynthesis in other plant species. In addition, three bHLH proteins were localized in the nucleus and displayed transcriptional activation activity in a yeast hybrid system. Through a series of functional experiments, we further demonstrated that PsbHLH1 could transcriptionally activate the expression of PsDFR and PsANS via directly binding to their promoters. These results laid a solid foundation to better understand the regulatory mechanisms of anthocyanin biosynthesis in P. suffruticosa and to benefit molecular breeding of tree peony cultivars with novel color.

Genetic Diversity of Paeonia rockii (Flare Tree Peony) Germplasm Accessions Revealed by Phenotypic Traits, EST-SSR Markers and Chloroplast DNA Sequences

Research Highlights: This study, based on the first collection of cultivated Paeonia rockii (flare tree peony, FTP) germplasm across the main distribution area by our breeding desires, comprehensively evaluates these accessions by using phenotypic traits, expressed sequence tag (EST)-simple sequence repeat (SSR) markers and chloroplast DNA sequences (cpDNA). The results show that these accessions collected selectively by us can represent the genetic background information of FTP as a germplasm of tree crops. Background and Objectives: FTP has high cultural, ornamental and medicinal value traditionally, as well as recently presenting a significance as an emerging edible oil with high α-linolenic acid contents in the seeds. The objectives of this study are to reveal the characteristics of the genetic diversity of FTP, as well as to provide scientific suggestions for the utilization of tree peony breeding and the conservation of germplasm resource. Materials and Methods: Based on the phenotypic traits, EST-SSR markers and chloroplast DNA sequence variation, we studied the diversity of a newly established population of 282 FTP accessions that were collected and propagated by ourselves in our breeding project in recent years. Results: (1) There was an abundant variation in phenotype of the accessions, and the phenotypic variation was evenly distributed within the population, without significant hierarchical structure, (2) the EST-SSR data showed that these 282 accessions had relatively high genetic diversity, in which a total of 185 alleles were detected in 34 pairs of primers. The 282 accessions were divided into three distinct groups, and (3) the chloroplast DNA sequences (cpDNA) data indicated that these accessions had a higher genetic diversity than the population level and a lower genetic diversity than the species level of wild P. rockii, and the existing spatial genetic structure of these accessions can be divided into two branches. Conclusions: From the results of the three analyses, we found that these accessions can fully reflect the genetic background information of FTP germplasm resources, so their protection and utilization will be of great significance for genetic improvement of woody peonies.

Advances in molecular biology of Paeonia L

Molecular biology can serve as a tool to solve the limitations of traditional breeding and cultivation techniques related to flower patterns, the improvement of flower color, and the regulation of flowering and stress resistance. These characteristics of molecular biology ensured its significant role in improving the efficiency of breeding and germplasm amelioration of Paeonia. This review describes the advances in molecular biology of Paeonia, including: (1) the application of molecular markers; (2) genomics, transcriptomics, proteomics, metabolomics, and microRNA studies; (3) studies of functional genes; and (4) molecular biology techniques. This review also points out select limitations in current molecular biology, analyzes the direction of Paeonia molecular biology research, and provides advice for future research objectives.

Chalcone synthase is ubiquitinated and degraded via interactions with a RING-H2 protein in petals of Paeonia 'He Xie'

Flavonoids are secondary metabolites widely distributed among angiosperms, where they play diverse roles in plant growth, development, and evolution. The regulation of flavonoid biosynthesis in plants has been extensively studied at the transcriptional level, but post-transcriptional, translational, and post-translational control of flavonoid biosynthesis remain poorly understood. In this study, we analysed post-translational regulation of flavonoid biosynthesis in the ornamental plant Paeonia, using proteome and ubiquitylome profiling, in conjunction with transcriptome data. Three enzymes involved in flavonoid biosynthesis were identified as being putative targets of ubiquitin-mediated degradation. Among these, chalcone synthase (PhCHS) was shown to have the greatest number of ubiquitination sites. We examined PhCHS abundance in petals using PhCHS-specific antibody and found that its accumulation decreased at later developmental stages, resulting from 26S proteasome-mediated degradation. We further identified a ring domain-containing protein (PhRING-H2) that physically interacts with PhCHS and demonstrated that PhRING-H2 is required for PhCHS ubiquitination. Taken together, our results suggest that PhRING-H2-mediates PhCHS ubiquitination and degradation is an important mechanism of post-translational regulation of flavonoid biosynthesis in Paeonia, providing a theoretical basis for the manipulation of flavonoid biosynthesis in plants.

Comparative transcriptome analyses reveal genes related to pigmentation in the petals of red and white Primula vulgaris cultivars

Primula vulgaris is an important ornamental plant species with various flower color. To explore the molecular mechanism of its color formation, comparative transcriptome analyses of the petals in red and white cultivars was performed. A total of 4451 differentially expressed genes were identified and annotated into 128 metabolic pathways. Candidate genes FLS, F3'H, DFR, ANS and AOMT in the anthocyanin pathway were expressed significantly higher in the red cultivar than the white and may be responsible for the red coloration. In the red petals, a putative transcription factors bHLH (c52273.graph_c0) was up-regulated about 14-fold, while a R2R3-MYB unigene (c36140.graph_c0) was identified as a repressor involved in anthocyanin regulation and was significantly down-regulated. In addition, the anatomy analyses and pigments composition in the red and white petals were also analyzed. The papillae on the adaxial epidermis of the red petals of P. vulgaris display a triangle-shapes, in contrast with a spherical shape for the white petals. Although flavonoids were detected in both cultivars, anthocyanins could only be identified in the red cultivar. Gossypetin and peonidin/rosinin were the most abundant pigments in red petals. This study shed light on the genetic and biochemistry mechanisms underlying the flower coloration in Primula.

A Novel R2R3-MYB Transcription Factor Contributes to Petal Blotch Formation by Regulating Organ-Specific Expression of PsCHS in Tree Peony (Paeonia suffruticosa)

Flower color patterns play critical roles in plant-pollinator interactions and represent one of the most common adaptations during angiosperm evolution. However, the molecular mechanisms underlying flower color pattern formation are less understood in non-model organisms. The aim of this study was to identify genes involved in the formation of petal blotches in tree peony (Paeonia suffruticosa) through transcriptome profiling and functional experiments. We identified an R2R3-MYB gene, PsMYB12, representing a distinct R2R3-MYB subgroup, with a spatiotemporal expression pattern tightly associated with petal blotch development. We further demonstrated that PsMYB12 interacts with a basic helix-loop-helix (bHLH) and a WD40 protein in a regulatory complex that directly activates PsCHS expression, which is also specific to the petal blotches. Together, these findings advance our understanding of the molecular mechanisms of pigment pattern formation beyond model plants. They also benefit molecular breeding of tree peony cultivars with novel color patterns and promote germplasm innovation.

Identification of Two Novel R2R3-MYB Transcription factors, PsMYB114L and PsMYB12L, Related to Anthocyanin Biosynthesis in Paeonia suffruticosa

Flower color is a charming phenotype with very important ornamental and commercial values. Anthocyanins play a critical role in determining flower color pattern formation, and their biosynthesis is typically regulated by R2R3-MYB transcription factors (TFs). Paeonia suffruticosa is a famous ornamental plant with colorful flowers. However, little is known about the R2R3-MYB TFs that regulate anthocyanin accumulation in P. suffruticosa. In the present study, two R2R3-MYB TFs, namely, PsMYB114L and PsMYB12L, were isolated from the petals of P. suffruticosa 'Shima Nishiki' and functionally characterized. Sequence analysis suggested that PsMYB114L contained a bHLH-interaction motif, whereas PsMYB12L contained two flavonol-specific motifs (SG7 and SG7-2). Subsequently, the in vivo function of PsMYB114L and PsMYB12L was investigated by their heterologous expression in Arabidopsis thaliana and apple calli. In transgenic Arabidopsis plants, overexpression of PsMYB114L and of PsMYB12L caused a significantly higher accumulation of anthocyanins, resulting in purple-red leaves. Transgenic apple calli overexpressing PsMYB114L and PsMYB12L also significantly enhanced the anthocyanins content and resulted in a change in the callus color to red. Meanwhile, gene expression analysis in A. thaliana and apple calli suggested that the expression levels of the flavonol synthase (MdFLS) and anthocyanidin reductase (MdANR) genes were significantly downregulated and the dihydroflavonol 4-reductase (AtDFR) and anthocyanin synthase (AtANS) genes were significantly upregulated in transgenic lines of PsMYB114L. Moreover, the expression level of the FLS gene (MdFLS) was significantly downregulated and the DFR (AtDFR/MdDFR) and ANS (AtANS/MdANS) genes were all significantly upregulated in transgenic lines plants of PsMYB12L. These results indicate that PsMYB114L and PsMYB12L both enhance anthocyanin accumulation by specifically regulating the expression of some anthocyanin biosynthesis-related genes in different plant species. Together, these results provide a valuable resource with which to further study the regulatory mechanism of anthocyanin biosynthesis in P. suffruticosa and for the breeding of tree peony cultivars with novel and charming flower colors.

Transcriptomic Analysis Reveals Transcription Factors Related to Leaf Anthocyanin Biosynthesis in Paeonia qiui

Paeonia qiui is a wild species of tree peony. P. qiui has good ornamental value owing to its leaf color change in spring. So far, the molecular mechanism of leaf color change in P. qiui is unclear. This study analyzes the anthocyanin level and transcriptome of two different color stages in P. qiui leaf. The purplish-red leaf stage is rich in anthocyanin, while the green leaf has very low levels of anthocyanin. Transcriptome analysis reveals that 6678 differentially-expressed genes (DEGs) are up-regulated, and 14,667 are down-regulated in the purplish-red leaf. Among these DEGs, 40 MYB (v-myb avian myeloblastosis viral oncogene homolog) genes, 40 bHLH (MYC-like basic helix–loop–helix) genes, and 15 WD40 (WD-repeat protein) genes were found. Based on phylogenetic and alignment analysis with the deduced amino acid sequences with known transcription factors, Unigene0024459 (MYB1) is likely the R2R3-MYB that promotes anthocyanin biosynthesis; Unigene0050761 (MYB2) is likely the R2R3-MYB that represses anthocyanin biosynthesis; Unigene0005081 (bHLH1) and Unigene0006146 (WD40-1) are likely the bHLH and WD40 that participate in regulating anthocyanin biosynthesis. Additionally, quantitative RT-PCR results confirmed the transcriptome analyses for key genes.