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

Metabolome and Transcriptome Analyses Reveal Flower Color Differentiation Mechanisms in Various Sophora japonica L. Petal Types

Sophora japonica L. is an important landscaping and ornamental tree species throughout southern and northern parts of China. The most common color of S. japonica petals is yellow and white. In this study, S. japonica flower color mutants with yellow and white flag petals and light purple-red wing and keel petals were used for transcriptomics and metabolomics analyses. To investigate the underlying mechanisms of flower color variation in S. japonica 'AM' mutant, 36 anthocyanin metabolites were screened in the anthocyanin-targeting metabolome. The results demonstrated that cyanidins such as cyanidin-3-O-glucoside and cyanidin-3-O-rutinoside in the 'AM' mutant were the key metabolites responsible for the red color of the wing and keel petals. Transcriptome sequencing and differentially expressed gene (DEG) analysis identified the key structural genes and transcription factors related to anthocyanin biosynthesis. Among these, F3'5'H, ANS, UFGT79B1, bHLH, and WRKY expression was significantly correlated with the cyanidin-type anthocyanins (key regulatory factors affecting anthocyanin biosynthesis) in the flag, wing, and keel petals in S. japonica at various flower development stages.

Comprehensive analysis of metabolome and transcriptome reveals the mechanism of color formation in different leave of Loropetalum Chinense var. Rubrum

BACKGROUND

Loropetalum chinense var. rubrum (L. chinense var. rubrum) is a precious, coloured-leaf native ornamental plant in the Hunan Province. We found an L. chinense var. rubrum tree with three different leaf colours: GL (green leaf), ML (mosaic leaf), and PL (purple leaf). The mechanism of leaf coloration in this plant is still unclear. Therefore, this study aimed to identify the metabolites and genes involved in determining the colour composition of L. chinense var. rubrum leaves, using phenotypic/anatomic observations, pigment content detection, and comparative metabolomics and transcriptomics.

RESULTS

We observed that the mesophyll cells in PL were purple, while those in GL were green and those in ML were a mix of purple-green. The contents of chlorophyll a, b, carotenoids, and total chlorophyll in PL and ML were significantly lower than those in GL. While the anthocyanin content in PL and ML was significantly higher than that in GL. The metabolomics results showed the differences in the content of cyanidin 3-O-glucoside, delphinidin 3-O-glucoside, cyanidin 3,5-O-diglucoside, pelargonidin, and petunidin 3,5-diglucoside in ML, GL, and PL were significant. Considering that the change trend of anthocyanin content change was consistent with the leaf colour difference, we speculated that these compounds might influence the colour of L. chinense var. rubrum leaves. Using transcriptomics, we finally identified nine differentially expressed structural genes (one ANR (ANR1217); four CYP75As (CYP75A1815, CYP75A2846, CYP75A2909, and CYP75A1716); four UFGTs (UFGT1876, UFGT1649, UFGT1839, and UFGT3273) and nine transcription factors (two MYBs (MYB1057 and MYB1211), one MADS-box (MADS1235), two AP2-likes (AP2-like1779 and AP2-like2234), one bZIP (bZIP3720), two WD40s (WD2173 and WD1867) and one bHLH (bHLH1631) that might be related to flavonoid biosynthesis and then impacted the appearance of colour in L. chinense var. rubrum leaves.

CONCLUSION

This study revealed potential molecular mechanisms associated with leaf coloration in L. chinense var. rubrum by analyzing differential metabolites and genes related to the anthocyanin biosynthesis pathway. It also provided a reference for research on leaf colour variation in other ornamental plants.

Transcriptome profiles of yellowish-white and fuchsia colored flowers in the Rheum palmatum complex reveal genes related to color polymorphism

Flower color variation is ubiquitous in many plant species, and several studies have been conducted to elucidate the underlying molecular mechanism. There are two flower color variants (yellowish-white and fuchsia) in the Rheum palmatum complex, however, few studies have investigated this phenomenon. Here, we used transcriptome sequencing of the two color variants to shed light on the molecular and biochemical basis for these color morphs. Comparison of the two transcriptomes identified 9641 differentially expressed unigenes (DEGs), including 6477 up-regulated and 3163 down-regulated genes. Functional analyses indicated that several DEGs were related to the anthocyanin biosynthesis pathway, and the expression profiles of these DEGs were coincident with the qRT-PCR validation results, indicating that expression levels of structural genes have a profound effect on the color variation in the R. palmatum complex. Our results suggested that the interaction of transcription factors (MYB, bHLH and WRKY) also regulated the anthocyanin biosynthesis in the R. palmatum complex. Estimation of selection pressures using the dN/dS ratio showed that 1106 pairs of orthologous genes have undergone positive selection. Of these positively selected genes, 21 were involved in the anthocyanin biosynthetic pathway, indicating that they may encode the proteins for structural alteration and affect flower color in the R. palmatum complex.

Integrative analysis of transcriptome and metabolome reveals flavonoid biosynthesis regulation in Rhododendron pulchrum petals

BACKGROUND

Color is the major ornamental feature of the Rhododendron genus, and it is related to the contents of flavonoid in petals. However, the regulatory mechanism of flavonoid biosynthesis in Rhododendron pulchrum remains unknown. The transcriptome and metabolome analysis of Rhododendron pulchrum with white, pink and purple color in this study aimed to reveal the mechanism of flavonoid biosynthesis and to provide insight for improving the petal color.

RESULTS

Flavonoids and flavonols are the major components of flavonoid metabolites in R.pulchrum, such as laricitrin, apigenin, tricin, luteolin, isoorientin, isoscutellarein, diosmetin and their glycosides derivatives. With transcriptome and metabolome analysis, we found CHS, FLS, F3'H, F3'5'H, DFR, ANS, GT, FNS, IFR and FAOMT genes showed significantly differential expression in cultivar 'Zihe'. FNS and IFR were discovered to be associated with coloration in R.pulchrum for the first time. The FNS gene existed in the form of FNSI. The IFR gene and its related metabolites of medicarpin derivatives were highly expressed in purple petal. In cultivar 'Fenhe', up-regulation of F3'H and F3'5'H and down-regulation of 4CL, DFR, ANS, and GT were associated with pink coloration. With the transcription factor analysis, a subfamily of DREBs was found to be specifically enriched in pink petals. This suggested that the DREB family play an important role in pink coloration. In cultivars 'Baihe', flavonoid biosynthesis was inhibited by low expression of CHS, while pigment accumulation was inhibited by low expression of F3'5'H, DFR, and GT, which led to a white coloration.

CONCLUSIONS

By analyzing the transcriptome and metabolome of R.pulchrum, principal differential expression genes and metabolites of flavonoid biosynthesis pathway were identified. Many novel metabolites, genes, and transcription factors associated with coloration have been discovered. To reveal the mechanism of the coloration of different petals, a model of the flavonoid biosynthesis pathway of R.pulchrum was constructed. These results provide in depth information regarding the coloration of the petals and the flavonoid metabolism of R.pulcherum. The study of transcriptome and metabolome profiling gains insight for further genetic improvement in Rhododendron.

Integrated Transcriptome and Metabolome Analyses Reveal the Anthocyanin Biosynthesis Pathway in AmRosea1 Overexpression 84K Poplar

Populus alba × Populus glandulosa (84K poplar) is model material with excellent genetic engineering resource and ornamental value. In our study, AmRosea1 (Antirrhinum majus) was overexpressed in 84K poplar, and the transgenic 84K (AM) poplar with high content of anthocyanin exhibited red pigmentation leaves. The transcriptome analysis between wild type (WT) and AM showed that 170 differentially expressed genes (DEGs) (86 up-regulated and 84 down-regulated) were found, and some DEGs were involved in flavone and flavonol biosynthesis, flavonoid biosynthesis and anthocyanin biosynthesis. The metabolome analysis showed that 13 anthocyanins-related differentially accumulated metabolites (DAMs) were detected in AM. The correlation analysis between DEGs and DAMs were performed, and the results revealed that 18 DEGs, including 11 MYB genes, two BZ1 genes, one FG2 gene, one ANS gene, and three IF7MAT genes, were negatively or positively correlated with 13 DAMs. The phylogenetic analysis demonstrated that there was high homology between AmRosea1 and PagMYB113, and MYB113 co-expressed with BZ1, ANS and DFR directly. Our results elucidated the molecular mechanism of plant color change mediated by anthocyanin biosynthesis pathway, which laid the foundation for the development and utilization of colorful woody plant.

A single-base insertion in BoDFR1 results in loss of anthocyanins in green-leaved ornamental kale

A CRISPR/Cas9-based knockout assay verified that BoDFR1 drives anthocyanin accumulation in ornamental kale and that BoDFR2, an ortholog of BoDFR1, is redundant. Anthocyanins are widely distributed in nature and give plants their brilliant colors. Leaf color is an important trait for ornamental kale. In this study, we measured anthocyanin contents and performed transcriptome deep sequencing (RNA-seq) of leaves from pink and green ornamental kale. We observed substantial differences in the expression levels of the two DIHYDROFLAVONOL 4-REDUCTASE-encoding genes BoDFR1 (Bo9g058630) and its ortholog BoDFR2 (Bo2g116380) between green-leaved and pink-leaved kale by RNA-seq and RT-qPCR. We cloned and sequenced BoDFR1 and BoDFR2 from both types of kale. We identified a 1-bp insertion in BoDFR1 and a 2-bp insertion in BoDFR2 in green-leaved kale compared to the sequences obtained from pink-leaved kale, both mapping to the second exon of their corresponding gene and leading to premature termination of translation. To confirm the genetic basis of the absence of anthocyanins in green kale, we used CRISPR/Cas9 genome editing to separately knock out BoDFR1 or BoDFR2 in the pink-leaved ornamental kale inbred line P23. We detected very low accumulation of anthocyanins in the resulting mutants Bodfr1-1 and Bodfr1-2, while Bodfr2-1 and Bodfr2-2 had anthocyanin levels comparable to those of the wild-type. We conclude that the insertion in BoDFR1, rather than that in BoDFR2, underlies the lack of anthocyanins in green-leaved ornamental kale. This work provides insight into the function of DFR and will contribute to germplasm improvement of ornamental plants.

Multi-Omics and miRNA Interaction Joint Analysis Highlight New Insights Into Anthocyanin Biosynthesis in Peanuts (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is one of the most important economic and oil crops in the world. At present, peanut varieties with rich anthocyanin in testa are rare in the market, but the selection and breeding of varieties with the related traits has always attracted the attention of breeders. In this study, two peanut varieties with the pink and purple testa, G110 (G) and Z18-40 (Z) were used to conduct interaction joint analysis of multi-omics and miRNA-target gene. The anthocyanin content of Z18-40 was 7.49-8.62-folds higher than G110 on 30 DAF (days after flowering) and 45 DAF via Ultraviolet-visible Spectrophotometer (UV-5800, Shanghai, China). And then, a total of 14 candidate genes related with the anthocyanin biosynthesis were identified for correlation in different comparison groups (R ² ≥ 0.80), among of a novel gene Ah21440 related with hydroxycinnamoyl transferase (HCT) biosynthesis was identified. In addition, Cyanidin 3-O-glucoside (Kuromanin, pmb0550) was the only common differentially accumulated metabolite (DAM) identified using multi-omics joint analysis in G1_vs_G2, Z1_vs_Z2, G1_vs_Z1, and G2_vs_Z2, respectively. Correlation analysis of miRNA-target genes and DEGs in the transcriptome shows that, AhmiR2950, AhmiR398, AhmiR50, and AhmiR51 regulated to HCT and chalcone biosynthesis related candidate genes (Ah21440, AhCHS, AhCHI). Lastly, all of 14 candidate genes and 4 differentially expressed miRNAs were validated using quantitative real-time PCR (qRT-PCR), which trends were consistent with that of the former transcriptome data. The results provide important reference for in-depth research on the anthocyanin metabolism mechanism in peanut testa.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Functional Annotation of a Full-Length Transcriptome and Identification of Genes Associated with Flower Development in Rhododendronsimsii (Ericaceae)

Rhododendronsimsii is one of the top ten famous flowers in China. Due to its historical value and high aesthetic, it is widely popular among Chinese people. Various colors are important breeding objectives in Rhododendron L. The understanding of the molecular mechanism of flower color formation can provide a theoretical basis for the improvement of flower color in Rhododendron L. To generate the R.simsii transcriptome, PacBio sequencing technology has been used. A total of 833,137 full-length non-chimeric reads were obtained and 726,846 high-quality full-length transcripts were found. Moreover, 40,556 total open reading frames were obtained; of which 36,018 were complete. In gene annotation analyses, 39,411, 18,565, 16,102 and 17,450 transcriptions were allocated to GO, Nr, KEGG and COG databases, correspondingly. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with Protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT) and Coding Non Coding Index (CNCI) databases and observed 6170, 2265, 4084 and 1240 lncRNAs, respectively. Based on the results, most genes were enriched in the flavonoid biosynthetic pathway. The eight key genes on the anthocyanin biosynthetic pathway were further selected and analyzed by qRT-PCR. The F3'H and ANS showed an upward trend in the developmental stages of R. simsii. The highest expression of F3'5'H and FLS in the petal color formation of R. simsii was observed. This research provided a huge number of full-length transcripts, which will help to proceed genetic analyses of R.simsii. native, which is a semi-deciduous shrub.

Changes at a Critical Branchpoint in the Anthocyanin Biosynthetic Pathway Underlie the Blue to Orange Flower Color Transition in Lysimachia arvensis

Anthocyanins are the primary pigments contributing to the variety of flower colors among angiosperms and are considered essential for survival and reproduction. Anthocyanins are members of the flavonoids, a broader class of secondary metabolites, of which there are numerous structural genes and regulators thereof. In western European populations of Lysimachia arvensis, there are blue- and orange-petaled individuals. The proportion of blue-flowered plants increases with temperature and daylength yet decreases with precipitation. Here, we performed a transcriptome analysis to characterize the coding sequences of a large group of flavonoid biosynthetic genes, examine their expression and compare our results to flavonoid biochemical analysis for blue and orange petals. Among a set of 140 structural and regulatory genes broadly representing the flavonoid biosynthetic pathway, we found 39 genes with significant differential expression including some that have previously been reported to be involved in similar flower color transitions. In particular, F3'5'H and DFR, two genes at a critical branchpoint in the ABP for determining flower color, showed differential expression. The expression results were complemented by careful examination of the SNPs that differentiate the two color types for these two critical genes. The decreased expression of F3'5'H in orange petals and differential expression of two distinct copies of DFR, which also exhibit amino acid changes in the color-determining substrate specificity region, strongly correlate with the blue to orange transition. Our biochemical analysis was consistent with the transcriptome data indicating that the shift from blue to orange petals is caused by a change from primarily malvidin to largely pelargonidin forms of anthocyanins. Overall, we have identified several flavonoid biosynthetic pathway loci likely involved in the shift in flower color in L. arvensis and even more loci that may represent the complex network of genetic and physiological consequences of this flower color polymorphism.

Transcriptomic and chemical analyses to identify candidate genes involved in color variation of sainfoin flowers

BACKGROUND

Sainfoin (Onobrychis viciifolia Scop) is not only a high-quality legume forage, but also a nectar-producing plant. Therefore, the flower color of sainfoin is an important agronomic trait, but the factors affecting its flower phenotype are still unclear. To gain insights into the regulatory networks associated with metabolic pathways of coloration compounds (flavonoids or anthocyanins) and identify the key genes, we conducted a comprehensive analysis of the phenotype, metabolome and transcriptome of WF and AF of sainfoin.

RESULTS

Delphinidin, petunidin and malvidin derivatives were the main anthocyanin compounds in the AF of sainfoin. These substances were not detected in the WF of sainfoin. The transcriptomes of WF and AF in sainfoin at the S1 and S3 stages were obtained using the Illumina HiSeq4000 platform. Overall, 10,166 (4273 upregulated and 5893 downregulated) and 15,334 (8174 upregulated and 7160 downregulated) DEGs were identified in flowers at S1 and S3 stages, respectively (WF-VS-AF). KEGG pathway annotations showed that 6396 unigenes were annotated to 120 pathways and contained 866 DEGs at S1 stages, and 6396 unigenes were annotated to 131 pathways and included 1546 DEGs at the S3 stage. Nine DEGs belonging to the "flavonoid biosynthesis"and "phenylpropanoid biosynthesis" pathways involved in flower color formation were identified and verified by RT-qPCR analyses. Among these DEGs, 4CL3, FLS, ANS, CHS, DFR and CHI2 exhibited downregulated expression, and F3H exhibited upregulated expression in the WF compared to the AF, resulting in a decrease in anthocyanin synthesis and the formation of WF in sainfoin.

CONCLUSIONS

This study is the first to use transcriptome technology to study the mechanism of white flower formation in sainfoin. Our transcriptome data will be a great enrichment of the genetic information for sainfoin. In addition, the data presented herein will provide valuable molecular information for genetic breeding and provide insight into the future study of flower color polymorphisms in sainfoin.

Genome-wide transcriptome variation landscape in Ruta chalepensis organs revealed potential genes responsible for rutin biosynthesis

Ruta chalepensis L., most commonly known as 'fringed rue,' is an excellent and valuable bioactive plant that produces a range of complex flavonoids, of which rutin is the major compound present in this plant of great pharmaceutical and medicinal significance. The present study is a pioneering attempt to examine the changes in the transcriptomic landscape of leaf, stem, and root tissues and correlate this with rutin quantity in each tissue in order to identify the candidate genes responsible for rutin biosynthesis and to increase genomic resources in fringed rue. Comparative transcriptome sequencing of leaves, stems and roots were performed using the NovaSeq 6000 platform. The de novo transcriptome assembly generated 254,685 transcripts representing 154,018 genes with GC content of 42.60 % and N50 of 2280 bp. Searching assembled transcripts against UniRef90 and SwissProt databases annotated 79.7 % of them as protein coding. The leaf tissues had the highest rutin content followed by stems and roots. Several differentially expressed genes and transcripts relating to rutin biosynthesis were identified in leaves comparing with roots or stems comparing with roots. All the genes known to be involved in rutin biosynthesis showed up-regulation in leaves as compared with roots. These results were confirmed by gene ontology (GO) and pathway enrichment analyses. Up-regulated genes in leaves as compared with roots enriched GO terms with relation to rutin biosynthesis e.g. action of flavonol synthase, biosynthetic mechanism of malonyl-CoA, and action of monooxygenase. Phylogenetic analysis of the rhamnosyltransferase (RT) gene showed that it was highly homologues with RT sequence from Citrus species and all were located in the same clade. This transcriptomic dataset will serve as an important public resource for future genomics and transcriptomic studies in R. chalepensis and will act as a benchmark for the identification and genetic modification of genes involved in the biosynthesis of secondary metabolites.

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

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.