Transcriptome Sequencing Reveals Regulatory Mechanisms of Taxol Synthesis in Taxus wallichiana var. Mairei

Transcriptome and comparative chloroplast genome analysis of Taxus yunnanensis individuals with high and low paclitaxel yield

Paclitaxel is a potent anti-cancer drug that is mainly produced through semi-synthesis, which still requires plant materials as precursors. The content of paclitaxel and 10-deacetyl baccatin III (10-DAB) in Taxus yunnanensis has been found to differ from that of other Taxus species, but there is little research on the mechanism underlying the variation in paclitaxel content in T. yunnanensis of different provenances. In this experiment, the contents of taxoids and precursors in twigs between a high paclitaxel-yielding individual (TG) and a low paclitaxel-yielding individual (TD) of T. yunnanensis were compared, and comparative analyses of transcriptomes as well as chloroplast genomes were performed. High-performance liquid chromatography (HPLC) detection showed that 10-DAB and baccatin III contents in TG were 18 and 47 times those in TD, respectively. Transcriptomic analysis results indicated that genes encoding key enzymes in the paclitaxel biosynthesis pathway, such as taxane 10-β-hydroxylase (T10βH), 10-deacetylbaccatin III 10-O-acetyltransferase (DBAT), and debenzoyl paclitaxel N-benzoyl transferase (DBTNBT), exhibited higher expression levels in TG. Additionally, qRT-PCR showed that the relative expression level of T10βH and DBAT in TG were 29 and 13 times those in TD, respectively. In addition, six putative transcription factors were identified that may be involved in paclitaxel biosynthesis from transcriptome data. Comparative analysis of plastid genomes showed that the TD chloroplast contained a duplicate of rps12, leading to a longer plastid genome length in TD relative to TG. Fifteen mutation hotspot regions were identified between the two plastid genomes that can serve as candidate DNA barcodes for identifying high-paclitaxel-yield individuals. This experiment provides insight into the difference in paclitaxel accumulation among different provenances of T. yunnanensis individuals.

Mass spectrometry imaging and single-cell transcriptional profiling reveal the tissue-specific regulation of bioactive ingredient biosynthesis in Taxus leaves

Taxus leaves provide the raw industrial materials for taxol, a natural antineoplastic drug widely used in the treatment of various cancers. However, the precise distribution, biosynthesis, and transcriptional regulation of taxoids and other active components in Taxus leaves remain unknown. Matrix-assisted laser desorption/ionization-mass spectrometry imaging analysis was used to visualize various secondary metabolites in leaf sections of Taxus mairei, confirming the tissue-specific accumulation of different active metabolites. Single-cell sequencing was used to produce expression profiles of 8846 cells, with a median of 2352 genes per cell. Based on a series of cluster-specific markers, cells were grouped into 15 clusters, suggesting a high degree of cell heterogeneity in T. mairei leaves. Our data were used to create the first Taxus leaf metabolic single-cell atlas and to reveal spatial and temporal expression patterns of several secondary metabolic pathways. According to the cell-type annotation, most taxol biosynthesis genes are expressed mainly in leaf mesophyll cells; phenolic acid and flavonoid biosynthesis genes are highly expressed in leaf epidermal cells (including the stomatal complex and guard cells); and terpenoid and steroid biosynthesis genes are expressed specifically in leaf mesophyll cells. A number of novel and cell-specific transcription factors involved in secondary metabolite biosynthesis were identified, including MYB17, WRKY12, WRKY31, ERF13, GT_2, and bHLH46. Our research establishes the transcriptional landscape of major cell types in T. mairei leaves at a single-cell resolution and provides valuable resources for studying the basic principles of cell-type-specific regulation of secondary metabolism.

Full-Length Transcriptome Sequencing Combined with RNA-Seq to Analyze Genes Related to Terpenoid Biosynthesis in Cinnamomum burmannii

Cinnamomum burmannii is a cinnamomum plant rich in natural D-borneol. Natural D-borneol is a bicycle monoterpenoid compound widely used in the food, pharmaceutical, and cosmetic industries. Therefore, analyzing the biosynthesis mechanism of natural D-borneol in C. burmannii at the molecular level is helpful for directional breeding in the future and further development and utilization of C. burmannii and its related gene resources. In our study, 76 genes related to terpene metabolism were analyzed through third-generation sequencing and second-generation sequencing. Of these genes, 57 were associated with the synthesis of the terpenoid skeleton, and 19 belonged to terpenoid synthase, including four monoterpenoid synthases, seven sesquiterpenoid synthases, and eight diterpenoid synthases. Two genes in diterpenoid synthase were differentially expressed in high D-borneol and low D-borneol plants. It was speculated that these two genes might be related to D-borneol synthesis. How these two genes participate in the synthesis of D-borneol needs further study.

TcMYB29a, an ABA-Responsive R2R3-MYB Transcriptional Factor, Upregulates Taxol Biosynthesis in Taxus chinensis

Paclitaxel (Taxol), a highly modified diterpene agent mainly obtained from Taxus species, is the most widely used anticancer drug. Abscisic acid (ABA) is a well-known stress hormone that plays important roles in the secondary metabolism of plants, and it can also induce the accumulation of taxol in Taxus cell suspension cultures. However, the mechanism behind the regulation of taxol biosynthesis by ABA remains largely unknown. In previous research, a R2R3 MYB transcription factor (TF) TcMYB29a was observed to show a significant correlation with taxol biosynthesis, indicative of its potential role in the taxol biosynthesis. In this study, the TcMYB29a encoded by its gene was further characterized. An expression pattern analysis revealed that TcMYB29a was highly expressed in the needles and roots. Overexpression of TcMYB29a in Taxus chinensis cell suspension cultures led to an increased accumulation of taxol, and upregulated expression of taxol-biosynthesis-related genes, including the taxadiene synthase (TS) gene, the taxane 5α-hydroxylase (T5OH) gene, and the 3'-N-debenzoyl-2'-deoxytaxol-N-benzoyltransferase (DBTNBT) gene as compared to the controls. Chromatin immunoprecipitation (ChIP) assays, yeast one-hybrid (Y1H) assays, electrophoretic mobility shift assays (EMSAs), and dual-luciferase reporter assays verified that TcMYB29a could bind and activate the promoter of TcT5OH. Promoter sequence analysis of TcMYB29a revealed that its promoter containing an AERB site from -313 to -319 was a crucial ABA-responsive element. Subsequently, the ABA treatment assay showed that TcMYB29a was strongly upregulated at 6 h after ABA pretreatment. Furthermore, TcMYB29a was strongly suppressed at 3 h after the methyl jasmonate (MeJA) treatment and was depressed to the platform at 12 h. Taken together, these results reveal that TcMYB29a is an activator that improves the accumulation of taxol in Taxus chinensis cells through an ABA-medicated signaling pathway which is different from JA-medicated signaling pathways for the accumulation of taxol. These findings provide new insights into the potential regulatory roles of MYBs on the expression of taxol biosynthetic genes in Taxus.

Optimization of Taxol Extraction Process Using Response Surface Methodology and Investigation of Temporal and Spatial Distribution of Taxol in Taxus mairei

Taxus mairei is an important source for industrial extraction of taxol in China. However, the standard and steps of extraction are currently not uniform, which seriously affects the taxol yield. In the present study, the influence of four factors (methanol concentration, solid-liquid ratio, ultrasonic extraction temperature, and ultrasonic extraction time) on the taxol yield was successively explored in T. mairei. A response surface methodology (RSM) was used to optimize the extraction process based on the single-factor experiments above. The optimal conditions were as follows: methanol concentration was 90%, solid-liquid ratio was 1:15 (g/mL), ultrasonic extraction temperature was 40 °C and ultrasonic extraction time was 60 min. Moreover, the twigs and needles from T. mairei with different tree ages were treated by the optimum extraction process, which further revealed temporal and spatial distribution of taxol in the reproducible tissues. Interestingly, the taxol content was relatively higher in needles of T. ‘Jinxishan’ (a cultivar from T. mairei with yellow aril, FY), but was less in FY twigs. The accumulation of taxol in twigs and leaves of females (with red aril, FR) was significantly higher than that of males (M); however, the content showed a decreasing trend with the increasing tree ages. Therefore, it is suitable to increase the proportion of female trees especially the FY leaves as raw materials for the industrial production of taxol from T. mairei, and the tree ages should be better controlled at 3–7 years.

Identification and in-silico characterization of taxadien-5α-ol-O-acetyltransferase (TDAT) gene in Corylus avellana L

Paclitaxel^® (PC) is one of the most effective and profitable anti-cancer drugs. The most promising sources of this compound are natural materials such as tissue cultures of Taxus species and, more recently, hazelnut (Corylus avellana L.). A large part of the PC biosynthetic pathway in the yew tree and a few steps in the hazelnut have been identified. Since understanding the biosynthetic pathway of plant-based medicinal metabolites is an effective step toward their development and engineering, this paper aimed to identify taxadiene-5α-ol-O-acetyltransferase (TDAT) in hazelnut. TDAT is one of the key genes involved in the third step of the PC biosynthetic pathway. In this study, the TDAT gene was isolated using the nested-PCR method and then characterized. The cotyledon-derived cell mass induced with 150 μM of methyl jasmonate (MeJA) was utilized to isolate RNA and synthesize the first-strand cDNA. The full-length cDNA of TDAT is 1423 bp long and contains a 1302 bp ORF encoding 433 amino acids. The phylogenetic analysis of this gene revealed high homology with its ortholog genes in Quercus suber and Juglans regia. Bioinformatics analyses were used to predict the secondary and tertiary structures of the protein. Due to the lack of signal peptide, protein structure prediction suggested that this protein may operate at the cytoplasm. The homologous superfamily of the T5AT protein, encoded by TDAT, has two domains. The highest and lowest hydrophobicity of amino acids were found in proline 142 and lysine 56, respectively. T5AT protein fragment had 24 hydrophobic regions. The tertiary structure of this protein was designed using Modeler software (V.9.20), and its structure was verified based on the results of the Verify3D (89.46%) and ERRAT (90.3061) programs. The T5AT enzyme belongs to the superfamily of the transferase, and the amino acids histidine 164, cysteine 165, leucine 166, histidine 167, and Aspartic acid 168 resided at its active site. More characteristics of TDAT, which would aid PC engineering programs and maximize its production in hazelnut, were discussed.

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

Recent Research Progress in Taxol Biosynthetic Pathway and Acylation Reactions Mediated by Taxus Acyltransferases

Taxol is one of the most effective anticancer drugs in the world that is widely used in the treatments of breast, lung and ovarian cancer. The elucidation of the taxol biosynthetic pathway is the key to solve the problem of taxol supply. So far, the taxol biosynthetic pathway has been reported to require an estimated 20 steps of enzymatic reactions, and sixteen enzymes involved in the taxol pathway have been well characterized, including a novel taxane-10β-hydroxylase (T10βOH) and a newly putative β-phenylalanyl-CoA ligase (PCL). Moreover, the source and formation of the taxane core and the details of the downstream synthetic pathway have been basically depicted, while the modification of the core taxane skeleton has not been fully reported, mainly concerning the developments from diol intermediates to 2-debenzoyltaxane. The acylation reaction mediated by specialized Taxus BAHD family acyltransferases (ACTs) is recognized as one of the most important steps in the modification of core taxane skeleton that contribute to the increase of taxol yield. Recently, the influence of acylation on the functional and structural diversity of taxanes has also been continuously revealed. This review summarizes the latest research advances of the taxol biosynthetic pathway and systematically discusses the acylation reactions supported by Taxus ACTs. The underlying mechanism could improve the understanding of taxol biosynthesis, and provide a theoretical basis for the mass production of taxol.

Investigation of bioactivities of Taxus chinensis, Taxus cuspidata, and Taxus × media by gas chromatography-mass spectrometry

Taxus species have attracted much attention for their potency in cancer treatment. However, investigating the bioactivities of Taxus species is a complex task, due to their diversity, slow growth, and endangered state. The most important Taxus species in China are Taxus chinensis (T. chinensis), Taxus cuspidata (T. cuspidata), and Taxus × media (T. media), which mainly grow in the northeastern region. This article probes deep into the differences among the leaves of T. chinensis, T. cuspidata, and T. media, with the aid of gas chromatography-mass spectrometry (GC-MS). Through GC-MS, 162 compounds were detected in the samples and found to contain 35 bioactive metabolites. On this basis, 20 metabolites with significant bioactivities (antibiotic, antioxidant, anticancer, and antiaging effects) were identified via unsupervised learning of principal component analysis and supervised learning of partial least squares-discriminant analysis. The results show that T. media has the most prominent antibiotic, antioxidant, and anticancer effects, while T. cuspidata has the most diverse and abundant metabolites that slow down aging.

Alternative sources and metabolic engineering of Taxol: Advances and future perspectives

Paclitaxel is one of the strong plant-derived anti-cancer drugs that was first isolated from the Pacific yew. Despite many paclitaxel's clinical successes, the limited accessibility of paclitaxel for clinical trials is recognized as the most important challenge. Thus, researchers are continuously trying to find the innovative ways to meet the community's need for this medicine. In the first step, the alternative sources for Taxol supply were recognized, such as Taxus genus, other plant genera, and endophytic fungi. In the next step, the biosynthetic pathways of Taxol or related metabolites were manipulated in the original organisms, or introduced to heterologous systems and then were manipulated in them. Here, a range of metabolic manipulating approaches have been successfully developed to redirect the metabolic flux toward Taxol, including promoter engineering, enzyme engineering, overexpressing the bottleneck enzymes, over- or down-regulation of transcription factors, activation of the cryptic genes, removing/minimizing the flux for competing pathways, tunable regulation of the metabolic pathway, and increasing the supplies of precursors. In this review, we discuss research progress on the alternative Taxol sources and its metabolic manipulating, and we suggest recent challenges and future perspectives.

The R2R3-MYB transcription factor family in Taxus chinensis: identification, characterization, expression profiling and posttranscriptional regulation analysis

The MYB transcription factor family is one of the largest gene families playing regulatory roles in plant growth and development. The MYB family has been studied in a variety of plant species but has not been reported in Taxus chinensis. Here we identified 72 putative R2R3-MYB genes in T. chinensis using a comprehensive analysis. Sequence features, conversed domains and motifs were characterized. The phylogenetic analysis showed TcMYBs and AtMYBs were clustered into 36 subgroups, of which 24 subgroups included members from T. chinensis and Arabidopsis thaliana, while 12 subgroups were specific to one species. This suggests the conservation and specificity in structure and function of plant R2R3-MYBs. The expression of TcMYBs in various tissues and different ages of xylem were investigated. Additionally, miRNA-mediated posttranscriptional regulation analysis revealed that TcMYBs were the targets of miR858, miR159 and miR828, suggesting the posttranscriptional regulation of MYBs is highly conserved in plants. The results provide a basis for further study the role of TcMYBs in the regulation of secondary metabolites of T. chinensis.

Metabolic Variations of Flavonoids in Leaves of T. media and T. mairei Obtained by UPLC-ESI-MS/MS

The needles of Taxus species contain a large number of bioactive compounds, such as flavonoids. In the present study, the total flavonoid content in leaves of Taxus media and Taxus mairei was 19.953 and 14.464 mg/g, respectively. A total of 197 flavonoid metabolites (70 flavones, 42 flavonols, 26 flavone C-glycosides, 20 flavanones, 15 anthocyanins, 13 isoflavones, 6 flavonolignans, and 5 proanthocyanidins) were identified for the first time by a widely targeted Ultra Performance Liquid Chromatography-Electrospray Ionization-Tandem Mass Spectrometry (UPLC-ESI-MS/MS) method within the two Taxus species, containing 160 common metabolites, with 37 unique metabolites merely determined in T. mairei or T. media. Moreover, 42 differential flavonoid metabolites were screened in the two Taxus species, which showed specific metabolic patterns in isoflavonoid biosynthesis, anthocyanin biosynthesis, and flavone and flavonol biosynthesis pathways. Compared to T. mairei, a more activated phenylpropanoid pathway was found in T. media, which could be responsible for the higher content of total flavonoids in T. media. Our results provide new insights into the diversity of flavonoid metabolites between T. mairei and T. media, and provide a theoretical basis for the sufficient utilization of Taxus species and the development of novel drugs.