Integrated Metabolome and Transcriptome Analysis Uncovers the Role of Anthocyanin Metabolism in Michelia maudiae

Integrated transcriptome and metabolome analysis provides insights into anthocyanin biosynthesis in Cichorium intybus L

BACKGROUND

Chicory is a unique and nutritious vegetable crop. However, the molecular mechanisms underlying anthocyanin biosynthesis in chicory remain poorly understood. We combined transcriptomics and metabolomics analyses to explore the molecular basis of anthocyanin biosynthesis in red-budded (Z1) and yellow-budded (Z7) chicory.

RESULTS

Integrated transcriptomics and metabolomics analyses were performed to investigate the molecular basis of anthocyanin biosynthesis in chicory. A total of 26 key structural genes, including F3'H, DFR, CHS, and ANS, were identified and enriched in pathways such as flavonoid and anthocyanin biosynthesis. Additionally, 29 transcription factors were identified, including 11 MYB, five bHLH, and two WD40 transcription factors, with seven MYB genes upregulated and four genes downregulated, indicating their roles in regulating anthocyanin biosynthesis. Notably, the MYB transcription factor, CI35997, which is homologous to RLL2A in lettuce, was predicted to positively regulate anthocyanin biosynthesis. Other transcription factors, such as AP2/ERF, bZIP, NAC, and Trihelix, have also been implicated. Metabolomics analysis revealed that cyanidin derivatives were the main contributors to the red coloration of chicory buds, with cyanidin-3-O-(6-O-malonyl)-glucoside being the most abundant. Furthermore, a competitive relationship between lignin and anthocyanin biosynthesis was observed, wherein the downregulation of lignin-related genes enhanced anthocyanin accumulation.

CONCLUSIONS

This study identified key structural genes and transcription factors that offer molecular-level insights into anthocyanin biosynthesis in chicory. These findings provide valuable guidance for genetic improvement of chicory and other crops with high anthocyanin content.

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

Coloration differences in three Camellia reticulata Lindl. cultivars: 'Tongzimian', 'Shizitou' and 'Damanao'

Camellia reticulata Lindl., also known as Yunnan Camellia, is an important ornamental plant in China, especially for its large and stunning flowers. A comprehensive understanding of their coloration mechanisms can aid breeders in developing new cultivars and improving their ornamental value; however, it is still unclear in Yunnan Camellia, especially in mixed-color flowers. In this study, we conducted metabolic and transcriptomic comparison analyses to investigate the coloration differences in three Yunnan Camellia cultivars: C. reticulata 'Shizitou' (SZT), C. reticulata 'Damanao' (MN), and C. reticulata 'Tongzimian' (TZM). Our results revealed that the initial flowering stage may play a critical role in the color change of MN. Metabolome analysis demonstrated that cyanidin was the primary anthocyanin in SZT and MN's red region, while its content was low in TZM and MN's white region. According to the transcriptome analysis, the anthocyanins biosynthesis pathway was reconstructed in Yunnan Camellia, and the low expression of CHS was detected in TZM and MN's white region, while ANR maintained a high expression level, which may lead to the low content of cyanidin in them. Transcription factors MYBs, bHLH, and bZIP may play a key role in regulating anthocyanin-structural genes. The co-expression analysis showed that the meristem tissue may play a crucial role in the formation of the mixed white-red color in MN. Our study enriched the genetic basis of flower coloration differences in Yunnan Camellia which will be a valuable genomic resource to understanding the biology of coloration formation and for breeding the Camellia cultivars.

Phenolic Acids and Flavonoids Play Important Roles in Flower Bud Differentiation in Mikania micrantha: Transcriptomics and Metabolomics

Mikania micrantha is a highly invasive vine, and its ability to sexually reproduce is a major obstacle to its eradication. The long-distance dissemination of M. micrantha depends on the distribution of seeds; therefore, inhibiting M. micrantha flowering and seed production is an effective control strategy. The number of blooms of M. micrantha differs at different altitudes (200, 900, and 1300 m). In this study, we used a combination of metabolomics and transcriptomics methods to study the patterns of metabolite accumulation in the flower buds of M. micrantha. Using LC-MS/MS, 658 metabolites were found in the flower buds of M. micrantha at three different altitudes (200, 900, and 1300 m). Flavonoids and phenolic acids were found to be the main differential metabolites, and their concentrations were lower at 900 m than at 200 m and 1300 m, with the concentrations of benzoic acid, ferulic acid, and caffeic acid being the lowest. The biosynthesis pathways for flavonoids and phenolic compounds were significantly enriched for differentially expressed genes (DEGs), according to the results of transcriptome analysis. The production of flavonoid and phenolic acids was strongly linked with the expressions of phenylalanine ammonia-lyase (PAL), caffeoyl-CoA O-methyltransferase (COMT), and 4-coumarate-CoA ligase (4CL), according to the results of the combined transcriptome and metabolome analysis. These genes' roles in the regulation of distinct phenolic acids and flavonoids during M. micrantha bud differentiation are still unknown. This study adds to our understanding of how phenolic acids and flavonoids are regulated in M. micrantha flower buds at various altitudes and identifies regulatory networks that may be involved in this phenomenon, offering a new approach for the prevention and management of M. micrantha.

Identification of key genes responsible for green and white colored spathes in Anthurium andraeanum (Hort.)

Modern anthuriums, Anthurium andraeanum (Hort.) are among the most popular flowering plants and widely used for interior decoration. Their popularity is largely attributed to the exotic spathes with different colors. Previous studies have reported color development in red spathe cultivars, but limited information is available on key genes regulating white and green colored spathes. This study analyzed anthocyanin, chlorophyll, and carotenoid contents as well as transcript differences in spathes of eight cultivars that differed in spathe colors ranging from red to white and green. Results showed that increased expression of a transcription factor AaMYB2 was associated with elevated levels of anthocyanin in spathes, but decreased expression of AaMYB2 and increased expression of AaLAR (leucoanthocyanidin reductase) and AaANR (anthocyanidin reductase) were accompanied with the accumulation of colorless proanthocyanidin, thus the white spathe. As to the green colored spathe, chlorophyll content in the green spathe cultivar was substantially higher than the other cultivars. Correspondingly, transcripts of chlorophyll biosynthesis-related genes AaHemB (porphobilinogen synthase) and AaPor (protochlorophyllide oxidoreductase) were highly upregulated but almost undetectable in white and red spathes. The increased expression of AaHemB and AaPor was correlated with the expression of transcription factor AaMYB124. Subsequently, qRT-PCR analysis confirmed their expression levels in nine additional cultivars with red, white, and green spathes. A working model for the formation of white and green spathes was proposed. White colored spathes are likely due to the decreased expression of AaMYB2 which results in increased expression of AaLAR and AaANR, and the green spathes are attributed to AaMYB124 enhanced expression of AaHemB and AaPor. Further research is warranted to test this working model.

Transcriptomic and metabolomic analyses unveil the growth advantage mechanism conferred by heterosis of Michelia 'Zhongshanhanxiao'

Michelia compressa (Maxim.) Sarg. is one of the important timber trees in Taiwan province, P. R. China. Michelia 'Zhongshanhanxiao' is a group of variants found among the progeny of M. compressa that exhibit higher growth rates compared with normal individuals, with a significantly increased stem diameter and height, as well as enlarged leaves and flowers. However, the molecular mechanisms fostering the growth advantage and morphological variations are unknown and deserve further study. Through analysing the transcriptome, metabolome and physiological processes of leaves, we identified remarkable differences in gene expression and metabolic profiles between Michelia 'Zhongshanhanxiao' and both the maternal M. compressa and its normal progeny. These differences were widely associated with a plant-pathogen interaction, phenylpropanoid biosynthesis, cyanoamino acid metabolism, carbon fixation in photosynthetic organisms and plant hormone signal transduction. Additionally, physiological measurements showed that Michelia 'Zhongshanhanxiao' possesses stronger photosynthetic capacity and higher plant hormone content. These results suggest that the heterosis of Michelia 'Zhongshanhanxiao' is regulated by candidates related to cell division, resistance to pathogens and the accumulation of organic compounds. The findings of this study provide crucial information on the molecular mechanisms underlying the growth advantages conferred by heterosis in trees.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

Transcriptomic and Metabolomic Analyses Providing Insights into the Coloring Mechanism of Docynia delavayi

The metabolome and transcriptome profiles of three different variations of mature Docynia delavayi fruit were synthesized to reveal their fruit color formation mechanism. A total of 787 secondary metabolites containing 149 flavonoid metabolites, most of which were flavonoids and flavonols, were identified in the three variations using ultra performance liquid chromatography- tandem mass spectrometry (UPLC-MS/MS), and we found that the secondary metabolites cyanidin-3-O-galactoside and cyanidin-3-O-glucoside were the major coloring substances in D. delavayi. This was associated with the significant upregulation of the structural genes F3H and F3′H in the anthocyanin synthesis pathway and the control genes WRKY, MYB, bZIP, bHLH, and NAC in RP. F3′H expression may play a significant role in the selection of components for anthocyanin synthesis. Our results contribute to breeding and nutritional research in D. delavayi and provide insight into metabolite studies of the anthocyanin biosynthetic pathway.

Comparative Metabolome and Transcriptome Analysis of Anthocyanin Biosynthesis in White and Pink Petals of Cotton (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is one of the important fiber crops. Cotton flowers usually appear white (or cream-colored) without colored spots at the petal base, and turn pink on the next day after flowering. In this study, using a mutant showing pink petals with crimson spots at their base, we conducted comparative metabolome and transcriptome analyses to investigate the molecular mechanism of coloration in cotton flowers. Metabolic profiling showed that cyanidin-3-O-glucoside and glycosidic derivatives of pelargonidins and peonidins are the main pigments responsible for the coloration of the pink petals of the mutant. A total of 2443 genes differentially expressed (DEGs) between the white and pink petals were identified by RNA-sequencing. Many DEGs are structural genes and regulatory genes of the anthocyanin biosynthesis pathway. Among them, MYB21, UGT88F3, GSTF12, and VPS32.3 showed significant association with the accumulation of cyanidin-3-O-glucoside in the pink petals. Taken together, our study preliminarily revealed the metabolites responsible for the pink petals and the key genes regulating the biosynthesis and accumulation of anthocyanins in the pink petals. The results provide new insights into the biochemical and molecular mechanism underlying anthocyanin biosynthesis in upland cotton.

Comparative transcriptome analysis of Lupinus polyphyllus Lindl. provides a rich molecular resource for research on coloration mechanism

Lupinus polyphyllus is rich in color, making it a well-known horticultural ornamental plant. However, little is known about the genes related to anthocyanin and carotenoid biosynthesis in L. polyphyllus. In this study, transcriptome sequencing was performed on eight different colors of L. polyphyllus. A total of 1.13 billion clean reads were obtained and assembled into 89,124 unigenes, which were then aligned with six databases, resulting in the identification of 54,823 annotated unigenes. Among these unigenes, 76 and 101 were involved in the biosynthetic pathway of carotenoids and anthocyanins, respectively. In addition, 505 transcription factors were revealed, which belonged to the MYB, R2R3-MYB, NAC, bHLH, and WD40 families. A total of 6,700 differentially expressed genes (DEGs) were obtained by comparative transcriptome analysis. Among them, 17 candidate unigenes (four carotenoid genes, seven anthocyanin genes, and six TFs) were specifically up-regulated for one or more colors of L. polyphyllus. Eight representative candidate unigenes were analyzed by qRT-PCR. The findings enrich the transcriptome database of lupine, and provide a rich molecular resource for research on the coloration mechanism of L. polyphyllus.

Proteotranscriptomics - A facilitator in omics research

Applications in omics research, such as comparative transcriptomics and proteomics, require the knowledge of the species-specific gene sequence and benefit from a comprehensive high-quality annotation of the coding genes to achieve high coverage. While protein-coding genes can in simple cases be detected by scanning the genome for open reading frames, in more complex genomes exonic sequences are separated by introns. Despite advances in sequencing technologies that allow for ever-growing numbers of genomes, the quality of many of the provided genome assemblies do not reach reference quality. These non-contiguous assemblies with gaps and the necessity to predict splice sites limit accurate gene annotation from solely genomic data. In contrast, the transcriptome only contains transcribed gene regions, is devoid of introns and thus provides the optimal basis for the identification of open reading frames. The additional integration of proteomics data to validate predicted protein-coding genes further enriches for accurate gene models. This review outlines the principles of the proteotranscriptomics approach, discusses common challenges and suggests methods for improvement.

Combined Analysis of the Transcriptome and Metabolome Revealed the Mechanism of Petal Coloration in Bauhinia variegata

Bauhinia variegata petals are colorful, rich in anthocyanins, and have ornamental, nutritional, and medicinal value. However, the regulatory mechanism of anthocyanin accumulation in B. variegata remains unclear. In this study, a combined analysis of the metabolome and transcriptome was performed in red and white B. variegata cultivars in the early, middle, and blooming stages. A total of 46 different anthocyanins were identified, of which 27 showed marked differences in accumulation between the two cultivars, and contribute to their different petal colors. Malvidin 3-O-galactoside, peonidin 3-O-galactoside, cyanidin 3-O-glucoside, cyanidin 3-O-galactoside, and malvidin 3-O-glucoside were much more abundant in the second stage of flowering. In the blooming stage, except for the anthocyanins mentioned, delphinidin 3-O-galactoside and petunidin 3-O-galactoside were the most abundant anthocyanins in the red flowers, indicating that malvidin, peonidin, cyanidin, delphinidin, and petunidin were all responsible for the red color of petals in B. variegata. RNA sequencing identified 2,431 differentially expressed genes (DEGs), of which 26 were involved in the anthocyanin synthesis pathway. Correlations between the anthocyanin biosynthesis-related DEGs and anthocyanin contents were explored, and the DEGs involved in anthocyanin accumulation in B. variegata petals were identified. Eighteen of these DEGs encoded key catalytic enzymes, such as anthocyanidin reductase (ANR) and flavonoid-3'5'-hydroxylase (F3'5'H), and 17 of them encoded transcription factors (TFs) belonging to 14 families (including MYB, NAC, SPL, ERF, and CHR28). These results improve our understanding of the roles of anthocyanins, catalytic enzymes, and TFs in B. variegata petal-color expression.

Gene-Metabolite Network Analysis Revealed Tissue-Specific Accumulation of Therapeutic Metabolites in Mallotus japonicus

Mallotus japonicus is a valuable traditional medicinal plant in East Asia for applications as a gastrointestinal drug. However, the molecular components involved in the biosynthesis of bioactive metabolites have not yet been explored, primarily due to a lack of omics resources. In this study, we established metabolome and transcriptome resources for M. japonicus to capture the diverse metabolite constituents and active transcripts involved in its biosynthesis and regulation. A combination of untargeted metabolite profiling with data-dependent metabolite fragmentation and metabolite annotation through manual curation and feature-based molecular networking established an overall metabospace of M. japonicus represented by 2129 metabolite features. M. japonicus de novo transcriptome assembly showed 96.9% transcriptome completeness, representing 226,250 active transcripts across seven tissues. We identified specialized metabolites biosynthesis in a tissue-specific manner, with a strong correlation between transcripts expression and metabolite accumulations in M. japonicus. The correlation- and network-based integration of metabolome and transcriptome datasets identified candidate genes involved in the biosynthesis of key specialized metabolites of M. japonicus. We further used phylogenetic analysis to identify 13 C-glycosyltransferases and 11 methyltransferases coding candidate genes involved in the biosynthesis of medicinally important bergenin. This study provides comprehensive, high-quality multi-omics resources to further investigate biological properties of specialized metabolites biosynthesis in M. japonicus.

Unraveling the Regulatory Mechanism of Color Diversity in Camellia japonica Petals by Integrative Transcriptome and Metabolome Analysis

Camellia japonica petals are colorful, rich in anthocyanins, and possess important ornamental, edible, and medicinal value. However, the regulatory mechanism of anthocyanin accumulation in C. japonica is still unclear. In this study, an integrative analysis of the metabolome and transcriptome was conducted in five C. japonica cultivars with different petal colors. Overall, a total of 187 flavonoids were identified (including 25 anthocyanins), and 11 anthocyanins were markedly differentially accumulated among these petals, contributing to the different petal colors in C. japonica. Moreover, cyanidin-3-O-(6″-O-malonyl) glucoside was confirmed as the main contributor to the red petal phenotype, while cyanidin-3-O-rutinoside, peonidin-3-O-glucoside, cyanidin-3-O-glucoside, and pelargonidin-3-O-glucoside were responsible for the deep coloration of the C. japonica petals. Furthermore, a total of 12,531 differentially expressed genes (DEGs) and overlapping DEGs (634 DEGs) were identified by RNA sequencing, and the correlation between the expression level of the DEGs and the anthocyanin content was explored. The candidate genes regulating anthocyanin accumulation in the C. japonica petals were identified and included 37 structural genes (especially CjANS and Cj4CL), 18 keys differentially expressed transcription factors (such as GATA, MYB, bHLH, WRKY, and NAC), and 16 other regulators (mainly including transporter proteins, zinc-finger proteins, and others). Our results provide new insights for elucidating the function of anthocyanins in C. japonica petal color expression.

Integrated metabolome and transcriptome analysis of the anthocyanin biosynthetic pathway in relation to color mutation in miniature roses

BACKGROUND

Roses are famous ornamental plants worldwide. Floral coloration is one of the most prominent traits in roses and is mainly regulated through the anthocyanin biosynthetic pathway. In this study, we investigated the key genes and metabolites of the anthocyanin biosynthetic pathway involved in color mutation in miniature roses. A comparative metabolome and transcriptome analysis was carried out on the Neptune King rose and its color mutant, Queen rose, at the blooming stage. Neptune King rose has light pink colored petals while Queen rose has deep pink colored petals.

RESULT

A total of 190 flavonoid-related metabolites and 38,551 unique genes were identified. The contents of 45 flavonoid-related metabolites, and the expression of 15 genes participating in the flavonoid pathway, varied significantly between the two cultivars. Seven anthocyanins (cyanidin 3-O-glucosyl-malonylglucoside, cyanidin O-syringic acid, cyanidin 3-O-rutinoside, cyanidin 3-O-galactoside, cyanidin 3-O-glucoside, peonidin 3-O-glucoside chloride, and pelargonidin 3-O-glucoside) were found to be the major metabolites, with higher abundance in the Queen rose. Thirteen anthocyanin biosynthetic related genes showed an upregulation trend in the mutant flower, which may favor the higher levels of anthocyanins in the mutant. Besides, eight TRANSPARENT TESTA 12 genes were found upregulated in Queen rose, probably contributing to a high vacuolar sequestration of anthocyanins. Thirty transcription factors, including two MYB and one bHLH, were differentially expressed between the two cultivars.

CONCLUSIONS

This study provides important insights into major genes and metabolites of the anthocyanin biosynthetic pathway modulating flower coloration in miniature rose. The results will be conducive for manipulating the anthocyanin pathways in order to engineer novel miniature rose cultivars with specific colors.

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.

Integrated Analysis of the Transcriptome and Metabolome Reveals Genes Involved in Terpenoid and Flavonoid Biosynthesis in the Loblolly Pine (Pinus taeda L.)

Loblolly pine (Pinus taeda L.) is an important tree for afforestation with substantial economic and ecological value. Many metabolites with pharmacological activities are present in the tissues of P. taeda. However, the biosynthesis regulatory mechanisms of these metabolites are poorly understood. In the present study, transcriptome and metabolome analyses were performed on five tissues of P. taeda. A total of 40.4 million clean reads were obtained and assembled into 108,663 unigenes. These were compared with five databases, revealing 39,576 annotated unigenes. A total of 13,491 differentially expressed genes (DEGs) were observed in 10 comparison groups. Of these, 487 unigenes exhibited significantly different expressions in specific tissues of P. taeda. The DEGs were explored using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes metabolic pathway analysis. We identified 343 and 173 candidate unigenes related to the biosynthesis of terpenoids and flavonoids, respectively. These included 62 R2R3-MYB, 30 MYB, 15 WRKY, seven bHLH, seven ERF, six ZIP, five AP2, and one WD40 genes that acted as regulators in flavonoid and/or terpenoid biosynthesis. Additionally, metabolomics analysis detected 528 metabolites, among which 168 were flavonoids. A total of 493 differentially accumulated metabolites (DAMs) were obtained in 10 comparison groups. The 3,7-Di-O-methyl quercetin was differentially accumulated in all the comparison groups. The combined transcriptome and metabolome analyses revealed 219 DEGs that were significantly correlated with 45 DAMs. Our study provides valuable genomic and metabolome information for understanding P. taeda at the molecular level, providing a foundation for the further development of P. taeda-related pharmaceutical industry.

Transcriptome and metabolome analyses reveal the regulation of peel coloration in green, red Chinese prickly ash (Zanthoxylum L.)

Peel colour is an important external economic characteristic of Chinese prickly ash cultivars (Zanthoxylum bungeanum Maxim.). To gain insight into their coloration mechanisms, we performed an integrated analysis of green and red peels using combined metabolomic and transcriptomic analyses. Pelargonin-O-hexoside-O-rhamnoside-O-hexoside, pelargonidin 3,5-diglucoside, peonidin O-hexoside, cyanidin O-syringic acid and peonidin 3-O-glucoside were found to be the key anthocyanins. Transcriptome data indicated that the anthocyanidin synthase genes and UDP-glucose flavonoid 3-O-glucosytransferase genes were significantly increased to promote the redness of the peels. In addition, we discussed the role of R₂R₃-MYB transcription factors in coloration, of which the c80935 and c226097 genes may be the key regulatory factors for anthocyanin biosynthesis. Generally, this is the first study to identify and reveal the main anthocyanins in Chinese prickly ash peels during different developmental periods. The results of this research lay the foundation for understanding the regulation of coloration in Chinese prickly ash peels.

The complete chloroplast genome sequence of Magnolia maudiae

The first complete chloroplast genome (cpDNA) sequence of Magnolia maudiae was determined from Illumina HiSeq pair-end sequencing data in this study. The cpDNA is 160,205 bp in length, contains a large single-copy region (LSC) of 88,249 bp and a small single-copy region (SSC) of 18,806 bp, which were separated by a pair of inverted repeats (IR) regions of 26,575 bp. The genome contains 132 genes, including 87 protein-coding genes, 8 ribosomal RNA genes, and 37 transfer RNA genes. Further phylogenomic analysis showed that M. maudiae was close to Magnolia odora and Magnolia laevifolia in Magnolia genus.

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.