Functional characterization and correlation analysis of phenylalanine ammonia-lyase (PAL) in coumarin biosynthesis from Peucedanum praeruptorum Dunn

Genome-Wide Identification and Expression Profile Analysis of the Phenylalanine Ammonia-Lyase Gene Family in Hevea brasiliensis

The majority of the world's natural rubber comes from the rubber tree (Hevea brasiliensis). As a key enzyme for synthesizing phenylpropanoid compounds, phenylalanine ammonia-lyase (PAL) has a critical role in plant satisfactory growth and environmental adaptation. To clarify the characteristics of rubber tree PAL family genes, a genome-wide characterization of rubber tree PALs was conducted in this study. Eight PAL genes (HbPAL1-HbPAL8), which spread over chromosomes 3, 7, 8, 10, 12, 13, 14, 16, and 18, were found to be present in the genome of H. brasiliensis. Phylogenetic analysis classified HbPALs into groups I and II, and the group I HbPALs (HbPAL1-HbPAL6) displayed similar conserved motif compositions and gene architectures. Tissue expression patterns of HbPALs quantified by quantitative real-time PCR (qPCR) proved that distinct HbPALs exhibited varying tissue expression patterns. The HbPAL promoters contained a plethora of cis-acting elements that responded to hormones and stress, and the qPCR analysis demonstrated that abiotic stressors like cold, drought, salt, and H2O2-induced oxidative stress, as well as hormones like salicylic acid, abscisic acid, ethylene, and methyl jasmonate, controlled the expression of HbPALs. The majority of HbPALs were also regulated by powdery mildew, anthracnose, and Corynespora leaf fall disease infection. In addition, HbPAL1, HbPAL4, and HbPAL7 were significantly up-regulated in the bark of tapping panel dryness rubber trees relative to that of healthy trees. Our results provide a thorough comprehension of the characteristics of HbPAL genes and set the groundwork for further investigation of the biological functions of HbPALs in rubber trees.

Selection and Verification of Standardized Reference Genes of Angelica dahurica under Various Abiotic Stresses by Real-Time Quantitative PCR

In traditional Chinese medicine, Angelica dahurica is a valuable herb with numerous therapeutic applications for a range of ailments. There have not yet been any articles on the methodical assessment and choice of the best reference genes for A. dahurica gene expression studies. Real-time quantitative PCR (RT-qPCR) is widely employed as the predominant method for investigating gene expression. In order to ensure the precise determination of target gene expression outcomes in RT-qPCR analysis, it is imperative to employ stable reference genes. In this study, a total of 11 candidate reference genes including SAND family protein (SAND), polypyrimidine tract-binding protein (PTBP), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), TIP41-like protein (TIP41), cyclophilin 2 (CYP2), elongation factor 1 α (EF1α), ubiquitin-protein ligase 9 (UBC9), tubulin β-6 (TUB6), thioredoxin-like protein YLS8 (YLS8), and tubulin-α (TUBA) were selected from the transcriptome of A. dahurica. Subsequently, three statistical algorithms (geNorm, NormFinder, and BestKeeper) were employed to assess the stability of their expression patterns across seven distinct stimulus treatments. The outcomes obtained from these analyses were subsequently amalgamated into a comprehensive ranking using RefFinder. Additionally, one target gene, phenylalanine ammonia-lyase (PAL), was used to confirm the effectiveness of the selected reference genes. According to the findings of this study, the two most stable reference genes for normalizing the expression of genes in A. dahurica are TIP41 and UBC9. Overall, our research has determined the appropriate reference genes for RT-qPCR in A. dahurica and provides a crucial foundation for gene screening and identifying genes associated with the biosynthesis of active ingredients in A. dahurica.

A chromosome-scale genome of Peucedanum praeruptorum provide insights into Apioideae evolution and medicinal ingredient biosynthesis

Peucedanum praeruptorum Dunn, a traditional Chinese medicine rich in coumarin, belongs to the Apiaceae family. A high-quality assembled genome of P. praeruptorum is lacking, which has posed obstacles to functional identification and molecular evolution studies of genes associated with coumarin production. Here, a chromosome-scale reference genome of P. praeruptorum, an important medicinal and aromatic plant, was first sequenced and assembled using Oxford Nanopore Technologies and Hi-C sequencing. The final assembled genome size was 1.83 Gb, with a contig N50 of 11.12 Mb. The entire BUSCO evaluation and second-generation read comparability rates were 96.0 % and 99.31 %, respectively. Furthermore, 99.91 % of the genome was anchored to 11 pseudochromosomes. The comparative genomic study revealed the presence of 18,593 orthogroups, which included 476 species-specific orthogroups and 1211 expanded gene families. Two whole-genome duplication (WGD) events and one whole-genome triplication (WGT) event occurred in P. praeruptorum. In addition to the γ-WGT shared by core eudicots or most eudicots, the first WGD was shared by Apiales, while the most recent WGD was unique to Apiaceae. Our study demonstrated that WGD events that occurred in Apioideae highlighted the important role of tandem duplication in the biosynthesis of coumarins and terpenes in P. praeruptorum. Additionally, the expansion of the cytochrome P450 monooxygenase, O-methyltransferase, ATP-binding cassette (ABC) transporter, and terpene synthase families may be associated with the abundance of coumarins and terpenoids. Moreover, we identified >170 UDP-glucosyltransferase members that may be involved in the glycosylation post-modification of coumarins. Significant gene expansion was observed in the ABCG, ABCB, and ABCC subgroups of the ABC transporter family, potentially facilitating the transmembrane transport of coumarins after bolting. The P. praeruptorum genome provides valuable insights into the machinery of coumarin biosynthesis and enhances our understanding of Apiaceae evolution.

Functional and kinetics of two efficient phenylalanine ammonia lyase from Pyrus bretschneideri

BACKGROUND

The enzyme phenylalanine ammonia lyase (PAL) controls the transition from primary to secondary metabolism by converting L-phenylalanine (L-Phe) to cinnamic acid. However, the function of PAL in pear plants (Pyrus bretschneideri) has not yet been fully elucidated.

RESULTS

We identified three PAL genes (PbPAL1, PbPAL2 and PbPAL3) from the pear genome by exploring pear genome databases. The evolutionary tree revealed that three PbPALs were classified into one group. We expressed PbPAL1 and PbPAL2 recombinant proteins, and the purified PbPAL1 and PbPAL2 proteins showed strict substrate specificity for L-Phe, no activity toward L-Tyr in vitro, and modest changes in kinetics and enzyme characteristics. Furthermore, overexpression of PbAL1 and PbPAL1-RNAi, respectively, and resulted in significant changes in stone cell and lignin contents in pear fruits. The results of yeast one-hybrid (Y1H) assays that PbWLIM1 could bind to the conserved PAL box in the PbPAL promoter and regulate the transcription level of PbPAL2.

CONCLUSIONS

Our findings not only showed PbPAL's potential role in lignin biosynthesis but also laid the foundation for future studies on the regulation of lignin synthesis and stone cell development in pear fruit utilizing molecular biology approaches.

Genome sequencing-based transcriptomic analysis reveals novel genes in Peucedanum praeruptorum

BACKGROUND

Peucedanum praeruptorum Dunn, a traditional Chinese herbal medicine, contains coumarin and volatile oil components that have clinical application value. However, early bolting often occurs in the medicinal materials of Apiaceae plants. The rhizomes of the medicinal parts are gradually lignified after bolting, resulting in a sharp decrease in the content of coumarins. At present, the link between coumarin biosynthesis and early bolting in P. praeruptorum has not been elucidated.

RESULTS

Combining the genome sequencing and the previous transcriptome sequencing results, we reanalyzed the differential transcripts of P. praeruptorum before and after bolting. A total of 62,088 new transcripts were identified, of which 31,500 were unknown transcripts. Functional classification and annotation showed that many genes were involved in the regulation of transcription, defense response, and carbohydrate metabolic processes. The main domains are the pentatricopeptide repeat, protein kinase, RNA recognition motif, leucine-rich repeat, and ankyrin repeat domains, indicating their pivotal roles in protein modification and signal transduction. Gene structure analysis showed that skipped exon (SE) was the most dominant alternative splicing, followed by the alternative 3' splice site (A3SS) and the alternative 5' splice site (A5SS). Functional enrichment of differentially expressed genes showed that these differentially expressed genes mainly include transmembrane transporters, channel proteins, DNA-binding proteins, polysaccharide-binding proteins, etc. In addition, genes involved in peroxisome, hexose phosphate pathway, phosphatidylinositol signaling system, and inositol phosphate metabolism pathway were greatly enriched. A protein-protein interaction network analysis discoverd 1,457 pairs of proteins that interact with each other. The expression levels of six UbiA genes, three UGT genes, and four OMT genes were higher during the bolting stage. This observation suggests their potential involvement in the catalytic processes of prenylation, glycosylation, and methylation of coumarins, respectively. A total of 100 peroxidase (PRX) genes were identified being involved in lignin polymerization, but only nine PRX genes were highly expressed at the bolting stage. It is worth noting that 73 autophagy-related genes (ATGs) were first identified from the KEGG pathway-enriched genes. Some ATGs, such as BHQH00009837, BHQH00013830, and novel8944, had higher expression levels after bolting.

CONCLUSIONS

Comparative transcriptome analysis and large-scale genome screening provide guidance and new opinions for the identification of bolting-related genes in P. praeruptorum.

Full-length transcriptome combined with RNA sequence analysis of Fraxinus chinensis

BACKGROUND

The dry root or stem bark of Fraxinus chinensis is a famous herb Qin Pi which is known for its anti-inflammatory, analgesic, anti-tumor, liver protective and diuretic pharmacological effects, the fundamental chemical components are coumarin, phenylethanol glycosides and flavonoids. However, it is difficult to clarify the secondary metabolite synthesis pathway and key genes involved in the pathway because of lack genome information of Fraxinus chinensis.

OBJECTIVE

To generate a complete transcriptome of Fraxinus chinensis and to clarify the differentially expressed genes (DEGs) in leaves and stem barks.

METHODS

In this study, full-length transcriptome analysis and RNA-Seq were combined to characterize Fraxinus chinensis transcriptome.

RESULTS

A total of 69,145 transcripts were acquired and regarded as reference transcriptome, 67,441 transcripts (97.47%) were annotated to NCBI non-redundant protein (Nr), SwissProt, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and eukaryotic orthologous groups (KOG) databases. A total of 18,917 isoforms were annotated to KEGG database and classified to 138 biological pathways. In total, 10,822 simple sequence repeat (SSRs) and 11,319 resistance (R) gene were classified to 18 types, and 3947 transcription factors (TFs) were identified in full-length transcriptome analysis. Additionally, 15,095 DEGs were detected by RNA-seq in leaves and barks, including 4696 significantly up-regulated and 10,399 significantly down-regulated genes. And 254 transcripts were annotated into phenylpropane metabolism pathway containing 86 DEGs and ten of these enzyme genes were verified by qRT-PCR.

CONCLUSION

It laid the foundation for further exploration of the biosynthetic pathway of phenylpropanoids and related key enzyme genes.

A phenylalanine ammonia lyase from Fritillaria unibracteata promotes drought tolerance by regulating lignin biosynthesis and SA signaling pathway

Drought is one of the key threatening environmental factors for plant and agriculture. Phenylalanine ammonia lyase (PAL) is a key enzyme involved in plant defense against abiotic stress, however, the role of PAL in drought tolerance remains elusive. Here, a PAL member (FuPAL1) containing noncanonical Ala-Ser-Gly triad was isolated from Fritillaria unibracteata, one important alpine pharmaceutical plant. FuPAL1, mainly distributed in cytosol, was more conserved than FuCOMT and FuCHI at both nucleotide and amino acid levels. FuPAL1 was overexpressed in Escherichia coli and the purified recombinant FuPAL1 protein showed catalytic preference on L-Phe than L-Tyr. Homology modeling and site-mutation of FuPAL1 exhibited FuPAL1 took part in the ammonization process by forming MIO-like group, and Phe141, Ser208, Ileu218 and Glu490 played key roles in substrate binding and (or) catalysis. HPLC analysis showed that lignin and salicylic acid levels increased but total flavonoid levels decreased in FuPAL1 transgenic Arabidopsis compared to wild-type plants. Moreover, FuPAL1 transgenic Arabidopsis significantly enhanced its drought tolerance, which suggested that FuPAL1 mediated tolerance to drought by inducing the biosynthesis and accumulation of salicylic acid and lignin. Taken together, our results confirmed that the FuPAL1 played an important role in drought tolerance, and FuPAL1 might be a valuable target for genetic improvement of drought resistance in future.

Comparative Transcriptomics Unveil the Crucial Genes Involved in Coumarin Biosynthesis in Peucedanum praeruptorum Dunn

Peucedanum praeruptorum Dunn is a commonly used traditional Chinese medicine that is abundant in furano- and dihydropyrano coumarins. When P. praeruptorum reaches the bolting stage, the roots gradually lignified, and the content of coumarins declines rapidly. Non-bolting has always been a decisive factor for harvesting the P. praeruptorum materials. To evaluate the amount of coumarin components in unbolted and bolted P. praeruptorum, the variations of praeruptorin A, praeruptorin B, praeruptorin E, peucedanocoumarin I, and peucedanocoumarin II were determined. Additionally, 336,505 transcripts were obtained from the comparative transcriptome data. Among them, a total of 1,573 differentially expressed genes were screened out. To identify the critical genes involved in coumarin biosynthesis, comparative transcriptomics coupled with co-expression associated analysis was conducted. Finally, coumarin biosynthesis-related eighteen candidate genes were selected for the validation of qPCR. Additionally, a phylogenetic tree and the expression profile of ATP-binding cassette (ABC) transporters were constructed. To clarify the main genes in the regulation of coumarin biosynthesis, the interaction network of the co-expression genes from thirteen modules was constructed. Current results exhibited the significant increment of praeruptorin A, praeruptorin B and praeruptorin E in the bolted P. praeruptorum. Although, peucedanocoumarin I and peucedanocoumarin II were slightly increased. Besides the content of coumarins, the essential genes involved in the coumarin biosynthesis also exhibited an overall downward trend after bolting. Three peroxidases (PRXs) involved in the production of lignin monomers had been demonstrated to be downregulated. PAL, C4H, HCT, COMT, CCoAOMT, and some ABC transporters were dramatically downregulated at the bolting stage. These results indicated that the downregulation of coumarin biosynthetic genes in the bolted P. praeruptorum ultimately reduced the formation of coumarins. However, the mechanism through which bolting indirectly affects the formation of coumarin still needs extra functional verification.

Anthocyanin regulatory networks in Solanum tuberosum L. leaves elucidated via integrated metabolomics, transcriptomics, and StAN1 overexpression

BACKGROUND

Anthocyanins, which account for color variation and remove reactive oxygen species, are widely synthesized in plant tissues and organs. Using targeted metabolomics and nanopore full-length transcriptomics, including differential gene expression analysis, we aimed to reveal potato leaf anthocyanin biosynthetic pathways in different colored potato varieties.

RESULTS

Metabolomics analysis revealed 17 anthocyanins. Their levels varied significantly between the different colored varieties, explaining the leaf color differences. The leaves of the Purple Rose2 (PurpleR2) variety contained more petunidin 3-O-glucoside and malvidin 3-O-glucoside than the leaves of other varieties, whereas leaves of Red Rose3 (RedR3) contained more pelargonidin 3-O-glucoside than the leaves of other varieties. In total, 114 genes with significantly different expression were identified in the leaves of the three potato varieties. These included structural anthocyanin synthesis-regulating genes such as F3H, CHS, CHI, DFR, and anthocyanidin synthase and transcription factors belonging to multiple families such as C3H, MYB, ERF, NAC, bHLH, and WRKY. We selected an MYB family transcription factor to construct overexpression tobacco plants; overexpression of this factor promoted anthocyanin accumulation, turning the leaves purple and increasing their malvidin 3-o-glucoside and petunidin 3-o-glucoside content.

CONCLUSIONS

This study elucidates the effects of anthocyanin-related metabolites on potato leaves and identifies anthocyanin metabolic network candidate genes.

Dynamic metabolic and transcriptomic profiling reveals the biosynthetic characteristics of hydroxycinnamic acid amides (HCAAs) in sunflower pollen

Sunflower pollen is a natural nutritious food with a long history and multiple functions, however, the main chemical components apart from flavonoids and their biosynthesis processes have not been thoroughly investigated. In this study, seven hydroxycinnamic acid amides (HCAAs) (1-7) abundant in sunflower pollen were isolated and identified as one type of the pollen's main chemicals. For a comprehensive understanding of HCAA biosynthesis in Helianthus annuus flowers, RNA-seq, metabolomics, and key genes related to biosynthesis in the sunflower were studied. A large number of compounds at different sunflower growth stages (the 7th, 14th, 21st, and 28th days) and high expression levels of related genes in the transcriptome were detected. A molecular network was constructed to clarify the synthetic pathway of HCAAs, which revealed high transcriptional levels of spermidine hydroxycinnamoyl transferase genes (HaSHT2795 and HaSHT2436) in 14-21-days-old flowers. HaSHT2795 enzymes catalyze tri-coumaroylspermidine formation, and virus-induced gene silencing to inhibit HaSHT2795 and HaSHT2436 could significantly reduce the synthesis of hydroxycinnamic acid amides in sunflower pollen. HCAAs were inferred to be related to the formation of pollen walls and the health effects of pollen. Analyzing HCAA biosynthesis and accumulation in H. annuus pollen will be helpful to understand the functions of HCAAs in the development of pollen and its nutritional value.

Consequences of LED Lights on Root Morphological Traits and Compounds Accumulation in Sarcandra glabra Seedlings

This study evaluated the effects of different light spectra (white light; WL, blue light; BL and red light; RL) on the root morphological traits and metabolites accumulation and biosynthesis in Sarcandra glabra. We performed transcriptomic and metabolomic profiling by RNA-seq and ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS), respectively. When morphological features were compared to WL, BL substantially increased under-ground fresh weight, root length, root surface area, and root volume, while RL inhibited these indices. A total of 433 metabolites were identified, of which 40, 18, and 68 compounds differentially accumulated in roots under WL (WG) vs. roots under BL (BG), WG vs. roots under RL (RG), and RG vs. BG, respectively. In addition, the contents of sinapyl alcohol, sinapic acid, fraxetin, and 6-methylcoumarin decreased significantly in BG and RG. In contrast, chlorogenic acid, rosmarinyl glucoside, quercitrin and quercetin were increased considerably in BG. Furthermore, the contents of eight terpenoids compounds significantly reduced in BG. Following transcriptomic profiling, several key genes related to biosynthesis of phenylpropanoid-derived and terpenoids metabolites were differentially expressed, such as caffeic acid 3-O-methyltransferase) (COMT), hydroxycinnamoyl-CoA shikimate hydroxycinnamoyl transferase (HCT), O-methyltransferase (OMT), and 1-deoxy-D-xylulose-5-phosphate synthetase (DXS). In summary, our findings showed that BL was suitable for growth and accumulation of bioactive metabolites in root tissue of S. glabra. Exposure to a higher ratio of BL might have the potential to improve the production and quality of S. glabra seedlings, but this needs to be confirmed further.

De novo Transcriptome Sequencing Coupled With Co-expression Analysis Reveal the Transcriptional Regulation of Key Genes Involved in the Formation of Active Ingredients in Peucedanum praeruptorum Dunn Under Bolting Period

Peucedanum praeruptorum Dunn is a perennial and one-off flowering plant of the Peucedanum genus in Umbelliferae. The cultivated P. praeruptorum Dunn usually grows nutritionally in the first year and then moves into the reproductive growth in the second year. The lignification of the roots caused by bolting leads to the quality decline of crude materials. Since most of the previous studies have dealt with coumarin biosynthesis and identification of functional genes in P. praeruptorum, the scientific connotation of the inability that the bolted P. praeruptorum cannot be used medically is still unclear. Here, we employed a transcriptome sequencing combined with coexpression analysis to unearth the regulation mechanism of key genes related to coumarin synthesis in pre- and postbolting period, and to explore the mechanisms underlying the effects of bolting on the formation and transport of coumarins between the annual and biennial plants. Six cDNA libraries were constructed, and the transcripts were sequenced and assembled by Illumina Hiseq platform. A total of 336,505 unigenes were obtained from 824,129 non-redundant spliced transcripts. Unigenes (114,488) were annotated to the NCBI nr database, 119,017 and 10,475 unigenes were aligned to Gene Ontology (GO) functional groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. Differential expression analysis screened out a series of upregulated and downregulated genes related to the phenylpropanoid pathway. The heatmap clustering showed that the similar expression patterns were both observed in groups C vs. D and groups C vs. F. The WGCNA-based coexpression was performed to elucidate the module and trait relationship to unearth important genes related to the bolting process. Seven pivotal modules on the KEGG functional annotations suggested these genes were mainly enriched in the process of plant–pathogen interaction, plant hormone signal transduction, MAPK signaling pathway, α-linolenic acid metabolism, circadian rhythm, and phenylpropanoid pathway. Further analysis provided clues that the key genes of the phenylpropanoid pathway, the ABC transporters, the apoptosis-related and circadian rhythm regulatory genes may play pivotal roles in regulating bolting signaling, biosynthesis, and transportation of coumarins.

Transcriptomic analysis reveals the parallel transcriptional regulation of UV-B-induced artemisinin and flavonoid accumulation in Artemisia annua L

UV-B radiation is a pivotal photomorphogenic signal and positively regulates plant growth and metabolite biosynthesis. In order to elucidate the transcriptional regulation mechanism underlying UV-B-induced artemisinin and flavonoid biosynthesis in Artemisia annua, the transcriptional responses of A. annua L. leaves to UV-B radiation were analyzed using the Illumina transcriptome sequencing. A total of 10705 differentially expressed genes (DEGs) including 533 transcription factors (TFs), were identified. Based on the expression trends of the differentially expressed TFs as well as artemisinin and flavonoid biosynthesis genes, we speculated that TFs belonging to 6 clusters were most likely to be involved in the regulation of artemisinin and/or flavonoid biosynthesis. The regulatory relationship between TFs and artemisinin/flavonoid biosynthetic genes was further studied. Dual-LUC assays results showed that AaMYB6 is a positive regulator of AaLDOX which belongs to flavonoid biosynthesis pathway. In addition, we identified an R2R3 MYB TF, AaMYB4 which potentially mediated both artemisinin and flavonoid biosynthesis pathways by activating the expression of AaADS and AaDBR2 in artemisinin biosynthesis pathway and AaUFGT in flavonoid biosynthesis pathway. Overall, our findings would provide an insight into the elucidation of the parallel transcriptional regulation of artemisinin and flavonoid biosynthesis in A. annua L. under UV-B radiation.

De novo transcriptome assembly of Angelica dahurica and characterization of coumarin biosynthesis pathway genes

Angelica dahurica (Hoffm.) Benth. & Hook.f. ex Franch. & Sav (A. dahurica) is a famous Chinese herb known for the production of coumarins, important secondary metabolites with wide-ranging pharmacological activities. In particular, the methoxylated coumarins like those produced by A. dahurica are known for their anti-inflammatory, anti-cancer, and anti-oxidant pharmacological effects. However, the molecular mechanism of coumarin biosynthesis in A. dahurica has not been studied. Such investigation could help scientists harness the biosynthesis potential of methoxylated coumarins. Here we present, three transcriptomes corresponding to leaf, root, and stem tissues of A. dahurica. A total of 114,310 unigenes with an average length of 1118 bp were de novo assembled, and 81,404 (71.21%) of those unigenes were annotated. Then, 181 unigenes encoding the seven key enzymes involved were identified, for which COMT (Caffeic acid 3-O-methyltransferase) was spatially used in a phylogenetic analysis, and some of these key enzyme genes were verified by qRT-PCR. Differentially expressed genes and root-specific-expressed genes were identified, by comparing genes' profile activity between roots and other tissues. Furthermore, multiple genes encoding key enzymes or transcription factors related to coumarin biosynthesis were identified and analyzed. This study is the first to report comprehensive gene information of A. dahurica at the transcriptional level, and to distinguish candidate genes related to its biosynthesis of coumarin, thus laying a foundation for this pathway's further exploration in A. dahurica.

De novo transcriptome analysis of Cnidium monnieri (L.) Cuss and detection of genes related to coumarin biosynthesis

Cnidium monnieri (L.) Cuss (C. monnieri) is one of the most widely used traditional herbal medicines, exhibiting a wide range of pharmacological functions for treating asynodia, trichomonas vaginitis, and osphyalgia. Its important medicinal value comes from its abundance of coumarins. To identify genes involved in coumarin biosynthesis and accumulation, we analyzed transcriptome data from flower, leaf, root and stem tissues of C. monnieri. A total of 173,938 unigenes with a mean length of 1,272 bp, GC content of 38.79%, and N50 length of 2,121 bp were assembled using the Trinity program. Of these, 119,177 unigenes were annotated in public databases. We identified differentially expressed genes (DEGs) based on expression profile analysis. These DEGs exhibited higher expression levels in flower tissue than in leaf, stem or root tissues. We identified and analyzed numerous genes encoding enzymes involved in coumarin biosynthesis, and verified genes encoding key enzymes using quantitative real-time PCR. Our transcriptome data will make great contributions to research on C. monnieri and provide clues for identifying candidate genes involved in coumarin metabolic pathways.

Functional characterization of cinnamate 4-hydroxylase from Helianthus annuus Linn using a fusion protein method

Sunflower (Helianthus annuus L.) is an important oil crop, the secondary metabolites of it include many compounds such as flavonoids and lignin. However, the research on the biosynthesis of phenolic compounds in sunflowers is still scarce. Cinnamate 4-hydroxylase (C4H) belongs to the cytochrome P450-dependent monooxygenase family and is involved in the synthesis of many phenolic compounds, but C4H in sunflowers has not yet been cloned and functionally characterized. In this study, we screened three C4H genes from the sunflower transcriptome and genomic databases, named HaC4H1, HaC4H2, and, HaC4H3, respectively. In heterologous expression experiments, we had improved a method from previous studies by the addition of restriction sites to make it easier to express multiple C4H functions and suitable for in vitro activity verification. HaC4Hs without the N-terminal membrane anchor region was fused with a redox partner of Arabidopsis thaliana cytochrome P450 enzyme (CYP450) by the method and functionally expressed in E. coli and the results showed that these three enzymes catalyzed the formation of p-coumaric acid. To further investigate whether our fusion protein approach is applicable to other C4Hs, we used this method to explore the functions of C4H from Peucedanum praeruptorum and Angelica decursiva, and they can also convert trans-cinnamic acid to p-coumaric acid. The gene expression profile showed that all three HaC4H genes showed the highest transcription levels in the roots and might be up-regulated by MeJA. In summary, these results reveal the function of HaC4Hs in sunflower and provide a simpler way to explore C4H and even other cytochrome P450 enzymes in prokaryotic expression systems.

Two CYP71AJ enzymes function as psoralen synthase and angelicin synthase in the biosynthesis of furanocoumarins in Peucedanum praeruptorum Dunn

Psoralen synthase and angelicin synthase responsible for the formation of psoralen and angelicin in Peucedanum praeruptorum Dunn were identified and functionally characterized, respectively. Furanocoumarins were reported to possess several activities such as anticancer, anti-inflammatory and neuroprotective, and function as phytotoxin and allelochemical in plants. Furanocoumarins are the main bioactive ingredient in P. praeruptorum which is a commonly used traditional Chinese medicine. Phenylalanine ammonia lyase (PAL), 4-coumarate: CoA ligase (4CL), p-coumaroyl CoA 2'-hyfroxylase (C2'H) were cloned previously to elucidate the biosynthetic mechanism of coumarin lactone ring. However, the genes involved in complex coumarins in P. praeruptorum have not been explored. Herein, putative psoralen synthase CYP71AJ49 and angelicin synthase CYP71AJ51 were cloned from P. praeruptorum. In vivo and in vitro yeast assays were conducted to confirm their activities. Furthermore, the results of High Performance Liquid Chromatography-Electrospray Ionization Mass Spectrometry (HPLC-ESI-MS) verified that CYP71AJ49 catalyzed the conversion of marmesin to psoralen, and CYP71AJ51 catalyzed columbianetin to angelicin. Subsequently, the expression profile showed that CYP71AJ49 and CYP71AJ51 were easily affected by environmental conditions, especially UV and temperature. The genes tissue-specific expression and compounds tissue-specific distribution pattern indicated the existence of substance transport in P. praeruptorum. Phylogenetic analysis was conducted with 27 CYP71AJs, CYP71AJ49 and CYP71AJ51 were classified in I-4 and I-2, respectively. These results provide further insight to understand the biosynthetic mechanism of complex coumarins.