Role of female-predominant MYB39-bHLH13 complex in sexually dimorphic accumulation of taxol in Taxus media

Abstract

Taxus trees are major natural sources for the extraction of taxol, an anti-cancer agent used worldwide. Taxus media is a dioecious woody tree with high taxol yield. However, the sexually dimorphic accumulation of taxoids in T. media is largely unknown. Our study revealed high accumulation of taxoids in female T. media trees using a UPLC–MS/MS method. Thereafter, many differential metabolites and genes between female and male T. media trees were identified using metabolomic and transcriptomic analyses, respectively. Most of the taxol-related genes were predominantly expressed in female trees. A female-specific R2R3-MYB transcription factor gene, TmMYB39, was identified. Furthermore, bimolecular fluorescence complementation and yeast two-hybrid assays suggested the potential interaction between TmMYB39 and TmbHLH13. Several taxol biosynthesis-related promoter sequences were isolated and used for the screening of MYB recognition elements. The electrophoretic mobility shift assay indicated that TmMYB39 could bind to the promoters of the GGPPS, T10OH, T13OH, and TBT genes. Interaction between TmMYB39 and TmbHLH13 transactivated the expression of the GGPPS and T10OH genes. TmMYB39 might function in the transcriptional regulation of taxol biosynthesis through an MYB-bHLH module. Our results give a potential explanation for the sexually dimorphic biosynthesis of taxol in T. media.

References

[1] Jennewein S, Croteau R . Taxol: biosynthesis, molecular genetics, and biotechnological applications. Appl Microbiol Biotechnol. 2001; 57: 13-9.

[2] Sanchez-Munoz R, Bonfill M, Cusido RM et al. Advances in the regulation of in vitro paclitaxel production: methylation of a Y-patch promoter region alters BAPT gene expression in Taxus cell cultures . Plant Cell Physiol. 2018; 59: 2255-67.

[3] Howat S, Park B, Oh IS et al. Paclitaxel: biosynthesis, production and future prospects. New Biotechnol. 2014; 31: 242-5.

[4] Kuang X, Sun S, Wei J et al. Iso-Seq analysis of the Taxus cuspidata transcriptome reveals the complexity of taxol biosynthesis. BMC Plant Biol. 2019; 19: 210.

[5] Croteau R, Ketchum RE, Long RM et al. Taxol biosynthesis and molecular genetics. Phytochem Rev. 2006; 5: 75-97.

[6] Ansbacher T, Freud Y, Major DT . Slow-starter enzymes: role of active-site architecture in the catalytic control of the biosynthesis of taxadiene by taxadiene synthase. Biochemistry. 2018; 57: 3773-9.

[7] Kaspera R, Croteau R . Cytochrome P450 oxygenases of taxol biosynthesis. Phytochem Rev. 2006; 5: 433-44.

[8] Long RM, Lagisetti C, Coates RM et al. Specificity of the N-benzoyl transferase responsible for the last step of taxol biosynthesis. Arch Biochem Biophys. 2008; 477: 384-9.

[9] Zhang M, Chen Y, Nie L et al. Transcriptome-wide identification and screening of WRKY factors involved in the regulation of taxol biosynthesis in Taxus chinensis . Sci Rep. 2018; 8: 5197.

[10] Li S, Zhang P, Zhang M et al. Functional analysis of a WRKY transcription factor involved in transcriptional activation of the DBAT gene in Taxus chinensis . Plant Biol. 2013; 15: 19-26.

[11] Zhang M, Jin X, Chen Y et al. TcMYC2a, a basic helix-loop-helix transcription factor, transduces JA-signals and regulates taxol biosynthesis in Taxus chinensis . Front Plant Sci. 2018; 9: 863.

[12] Zhang M, Li S, Nie L et al. Two jasmonate-responsive factors, TcERF12 and TcERF15, respectively act as repressor and activator of tasy gene of taxol biosynthesis in Taxus chinensis . Plant Mol Biol. 2015; 89: 463-73.

[13] Cao Y, Li K, Li Y et al. MYB transcription factors as regulators of secondary metabolism in plants. Biology (Basel). 2020; 9: 61.

[14] Yu C, Luo X, Zhang C et al. Tissue-specific study across the stem of Taxus media identifies a phloem-specific TmMYB3 involved in the transcriptional regulation of paclitaxel biosynthesis. Plant J. 2020; 103: 95-110.

[15] Song H, Lei Y, Zhang S . Differences in resistance to nitrogen and phosphorus deficiencies explain male-biased populations of poplar in nutrient-deficient habitats. J Proteome. 2018; 178: 123-7.

[16] Zhang S, Tang D, Korpelainen H et al. Metabolic and physiological analyses reveal that Populus cathayana males adopt an energy-saving strategy to cope with phosphorus deficiency. Tree Physiol. 2019; 39: 1630-45.

[17] Zhang S, Zhou R, Zhao H et al. iTRAQ-based quantitative proteomic analysis gives insight into sexually different metabolic processes of poplars under nitrogen and phosphorus deficiencies. Proteomics. 2016; 16: 614-28.

[18] Wu X, Liu J, Meng Q et al. Differences in carbon and nitrogen metabolism between male and female Populus cathayana in response to deficient nitrogen. Tree Physiol. 2021; 41: 119-33.

[19] Liao Q, Du R, Gou J et al. The genomic architecture of the sex-determining region and sex-related metabolic variation in Ginkgo biloba . Plant J. 2020; 104: 1399-409.

[20] Mukherjee S, Ghosh B, Jha TB et al. Variation in content of taxol and related taxanes in eastern Himalayan populations of Taxus wallichiana . Planta Med. 2002; 68: 757-9.

[21] Nadeem M, Rikhari HC, Kumar A et al. Taxol content in the bark of Himalayan yew in relation to tree age and sex. Phytochemistry. 2002; 60: 627-31.

[22] An JP, Wang X-F, Zhang X-W et al. An apple MYB transcription factor regulates cold tolerance and anthocyanin accumulation and undergoes MIEL1-mediated degradation. Plant Biotechnol J. 2020; 18: 337-53.

[23] Li L, Zhang Y, Luo J et al. Sex-specific responses of Populus yunnanensis exposed to elevated CO2 and salinity. Physiol Plant. 2013; 147: 477-88.

[24] Korgiopoulou C, Bresta P, Nikolopoulos D et al. Sex-specific structural and functional leaf traits and sun-shade acclimation in the dioecious tree Pistacia vera (Anacardiaceae) . Funct Plant Biol. 2019; 46: 649-59.

[25] Cossard GG, Toups MA, Pannell JR . Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb. Ann Bot. 2019; 123: 1119-31.

[26] Yu C, Zhang C, Xu X et al. Omic analysis of the endangered Taxaceae species Pseudotaxus chienii revealed the differences in taxol biosynthesis pathway between Pseudotaxus and Taxus yunnanensis trees . BMC Plant Biol. 2021; 21: 104.

[27] Zhou T, Luo X, Yu C et al. Transcriptome analyses provide insights into the expression pattern and sequence similarity of several taxol biosynthesis-related genes in three Taxus species. BMC Plant Biol. 2019; 19: 33-3.

[28] Sun YQ, Zhao W, Xu C-Q et al. Genetic variation related to high elevation adaptation revealed by common garden experiments in Pinus yunnanensis . Front Genet. 2019; 10: 1405.

[29] Mutanda I, Li J, Xu F et al. Recent advances in metabolic engineering, protein engineering, and transcriptome-guided insights toward synthetic production of taxol. Front Bioeng Biotechnol. 2021; 9: 632269.

[30] Sabzehzari M, Zeinali M, Naghavi MR . Alternative sources and metabolic engineering of taxol: advances and future perspectives. Biotechnol Adv. 2020; 43: 107569.

[31] Baebler S, Camloh M, Kovac M et al. Jasmonic acid stimulates taxane production in cell suspension culture of yew (Taxus × media) . Planta Med. 2002; 68: 475-6.

[32] Walker K, Long R, Croteau R . The final acylation step in taxol biosynthesis: cloning of the taxoid C13-side-chain N-benzoyltransferase from Taxus . Proc Natl Acad Sci USA. 2002; 99: 9166-71.

[33] Dubos C, Stracke R, Grotewold E et al. MYB transcription factors in Arabidopsis . Trends Plant Sci. 2010; 15: 573-81.

[34] Kirik V, Kolle K, Wohlfarth T et al. Ectopic expression of a novel MYB gene modifies the architecture of the Arabidopsis inflorescence . Plant J. 1998; 13: 729-42.

[35] Mu RL, Cao Y-R, Liu Y-F et al. An R2R3-type transcription factor gene AtMYB59 regulates root growth and cell cycle progression in Arabidopsis . Cell Res. 2009; 19: 1291-304.

[36] He C, Teixeira da Silva JA, Wang H et al. Mining MYB transcription factors from the genomes of orchids (Phalaenopsis and Dendrobium) and characterization of an orchid R2R3-MYB gene involved in water-soluble polysaccharide biosynthesis. Sci Rep. 2019; 9: 13818.

[37] Jeena GS, Kumar S, Shukla RK . Characterization of MYB35 regulated methyl jasmonate and wound responsive geraniol 10-hydroxylase-1 gene from Bacopa monnieri . Planta. 2021; 253: 89.

[38] Zhou W, Shi M, Deng C et al. The methyl jasmonate-responsive transcription factor SmMYB1 promotes phenolic acid biosynthesis in Salvia miltiorrhiza . Hortic Res. 2021; 8: 10.

[39] Li Y, Chen X, Wang J et al. Two responses to MeJA induction of R2R3-MYB transcription factors regulate flavonoid accumulation in Glycyrrhiza uralensis Fisch. PLoS One. 2020; 15: e0236565.

[40] Ramsay NA, Glover BJ . MYB-bHLH-WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005; 10: 63-70.

[41] Wang XC, Wu J, Guan M-L et al. Arabidopsis MYB4 plays dual roles in flavonoid biosynthesis . Plant J. 2020; 101: 637-52.

[42] Meng Y, Wang Z, Wang Y et al. The MYB activator WHITE PETAL1 associates with MtTT8 and MtWD40-1 to regulate carotenoid-derived flower pigmentation in Medicago truncatula . Plant Cell. 2019; 31: 2751-67.

[43] Lloyd A, Brockman A, Aguirre L et al. Advances in the MYB-bHLH-WD repeat (MBW) pigment regulatory model: addition of a WRKY factor and co-option of an anthocyanin MYB for betalain regulation. Plant Cell Physiol. 2017; 58: 1431-41.

[44] Yu C, Zhan X, Zhang C et al. Comparative metabolomic analyses revealed the differential accumulation of taxoids, flavonoids and hormones among six Taxaceae trees. Sci Hortic. 2021; 285: 110196.

[45] Hao J, Guo H, Shi X et al. Comparative proteomic analyses of two Taxus species (Taxus × media and Taxus mairei) reveals variations in the metabolisms associated with paclitaxel and other metabolites . Plant Cell Physiol. 2017; 58: 1878-90.

[46] Eisenreich W, Menhard B, Hylands PJ et al. Studies on the biosynthesis of taxol: the taxane carbon skeleton is not of mevalonoid origin. Proc Natl Acad Sci USA. 1996; 93: 6431-6.

[47] Koepp AE, Hezari M, Zajicek J et al. Cyclization of geranylgeranyl diphosphate to taxa-4(5),11(12)-diene is the committed step of taxol biosynthesis in Pacific yew. J Biol Chem. 1995; 270: 8686-90.

[48] Vandermoten S, Haubruge E, Cusson M . New insights into short-chain prenyltransferases: structural features, evolutionary history and potential for selective inhibition. Cell Mol Life Sci. 2009; 66: 3685-95.

[49] Yu C, Guo H, Zhang Y et al. Identification of potential genes that contributed to the variation in the taxoid contents between two Taxus species (Taxus media and Taxus mairei) . Tree Physiol. 2017; 37: 1659-71.

[50] Grabherr MG, Haas BJ, Yassour M et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29: 644-52.

[51] Zhou T, Luo X, Zhang C et al. Comparative metabolomic analysis reveals the variations in taxoids and flavonoids among three Taxus species. BMC Plant Biol. 2019; 19: 529.

[52] Yu C, Luo X, Zhan X et al. Comparative metabolomics reveals the metabolic variations between two endangered Taxus species (T. fuana and T. yunnanensis) in the Himalayas . BMC Plant Biol. 2018; 18: 197.