Identification of WRKY transcription factors involved in regulating the biosynthesis of the anti-cancer drug camptothecin in Ophiorrhiza pumila

Abstract

Camptothecin is a chemotherapeutic drug widely used to treat various cancers. Ophiorrhiza pumila is an ideal plant model for the study of camptothecin production, with various advantages for studying camptothecin biosynthesis and regulation. The DNA-binding WRKY transcription factors have a key regulatory role in secondary metabolite biosynthesis in plants. However, little is currently known about their involvement in camptothecin biosynthesis in O. pumila. We identified 46 OpWRKY genes unevenly distributed on the 11 chromosomes of O. pumila. Phylogenetic and multiple sequence alignment analyses divided the OpWRKY proteins into three subfamilies. Based on spatial expression and co-expression, we targeted the candidate gene OpWRKY6. Overexpression of OpWRKY6 significantly reduced the accumulation of camptothecin compared with the control. Conversely, camptothecin accumulation increased in OpWRKY6 knockout lines. Further biochemical assays showed that OpWRKY6 negatively regulates camptothecin biosynthesis from both the iridoid and shikimate pathways by directly downregulating the gene expression of OpGES, Op10HGO, Op7DLH, and OpTDC. Our data provide direct evidence for the involvement of WRKYs in the regulation of camptothecin biosynthesis and offer valuable information for enriching the production of camptothecin in plant systems.

References

[1] Rai A, Hirakawa H, Nakabayashi R et al. Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis. Nat Commun. 2021; 12: 405.

[2] Martino E, Volpe SD, Terribile E et al. The long story of camptothecin: from traditional medicine to drugs. Bioorg Med Chem Lett. 2017; 27: 701-7.

[3] Sirikantaramas S, Yamazaki M, Saito K . Mutations in topoisomerase I as a self-resistance mechanism coevolved with the production of the anticancer alkaloid camptothecin in plants. Proc Natl Acad Sci U S A. 2008; 105: 6782-6.

[4] Sirikantaramas S, Meeprasert A, Rungrotmongkol T et al. Structural insight of DNA topoisomerases I from camptothecin-producing plants revealed by molecular dynamics simulations. Phytochemistry. 2015; 113: 50-6.

[5] Ruan Q, Patel G, Wang J et al. Current advances of endophytes as a platform for production of anti-cancer drug camptothecin. Food Chem Toxicol. 2021; 151: 112113.

[6] Yang MQ, Wang Q, Liu Y et al. Divergent camptothecin biosynthetic pathway in Ophiorrhiza pumila . BMC Biol. 2021; 19: 122.

[7] Hao X, Xie C, Ruan Q et al. The transcription factor OpWRKY2 positively regulates the biosynthesis of the anticancer drug camptothecin in Ophiorrhiza pumila . Hortic Res. 2021; 8: 7.

[8] Kai GY, Wu C, Gen L et al. Biosynthesis and biotechnological production of anti-cancer drug camptothecin. Phytochem Rev. 2015; 14: 525-39.

[9] Sadre R, Magallanes-Lundback M, Pradhan S et al. Metabolite diversity in alkaloid biosynthesis: a multilane (Diastereomer) highway for camptothecin synthesis in Camptotheca acuminata . Plant Cell. 2016; 28: 1926-44.

[10] Asano T, Kobayashi K, Kashihara E et al. Suppression of camptothecin biosynthetic genes results in metabolic modification of secondary products in hairy roots of Ophiorrhiza pumila . Phytochemistry. 2013; 91: 128-39.

[11] Yamazaki Y, Kitajima M, Arita M et al. Biosynthesis of camptothecin. In silico and in vivo tracer study from [1-13C] glucose. Plant Physiol. 2004; 134: 161-70.

[12] Miettinen K, Dong L, Navrot N et al. The seco-iridoid pathway from Catharanthus roseus . Nat Commun. 2014; 5: 3606.

[13] Shi M, Gong H, Cui L et al. Targeted metabolic engineering of committed steps improves anti-cancer drug camptothecin production in Ophiorrhiza pumila hairy roots. Ind Crop Prod. 2020; 148: 112277.

[14] Yamazaki M, Mochida K, Asano T et al. Coupling deep transcriptome analysis with untargeted metabolic profiling in Ophiorrhiza pumila to further the understanding of the biosynthesis of the anti-cancer alkaloid camptothecin and anthraquinones. Plant Cell Physiol. 2013; 54: 686-96.

[15] Cui L, Ni X, Teng X et al. Co-overexpression of geraniol-10-hydroxylase and strictosidine synthase improves anti-cancer drug camptothecin accumulation in Ophiorrhiza pumila . Sci Rep. 2015; 5: 8227.

[16] Suttipanta N, Pattanaik S, Kulshrestha M et al. The transcription factor CrWRKY1 positively regulates the terpenoid indole alkaloid biosynthesis in Catharanthus roseus . Plant Physiol. 2011; 157: 2081-93.

[17] Chen F, Hu Y, Vannozzi A et al. The WRKY transcription factor family in model plants and crops. Crit Rev Plant Sci. 2017; 36: 311-35.

[18] Eulgem T, Rushton PJ, Robatzek S et al. The WRKY superfamily of plant transcription factors. Trends Plant Sci. 2000; 5: 199-206.

[19] Rushton PJ, Somssich IE, Ringler P et al. WRKY transcription factors. Trends Plant Sci. 2010; 15: 247-58.

[20] Xie Z, Zhang ZL, Zou X et al. Annotations and functional analyses of the Rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells. Plant Physiol. 2005; 137: 176-89.

[21] Wu KL, Guo ZJ, Wang HH et al. The WRKY family of transcription factors in Rice and Arabidopsis and their origins. DNA Res. 2005; 12: 9-26.

[22] Eulgem T, Somssich IE . Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007; 10: 366-71.

[23] Yamasaki K, Kigawa T, Inoue M et al. Solution structure of an Arabidopsis WRKY DNA binding domain. Plant Cell. 2005; 17: 944-56.

[24] Chi YJ, Yang Y, Zhou Y et al. Protein-protein interactions in the regulation of WRKY transcription factors. Mol Plant. 2013; 6: 287-300.

[25] Schluttenhofer C, Yuan L . Regulation of specialized metabolism by WRKY transcription factors. Plant Physiol. 2015; 167: 295-306.

[26] Cheng MN, Huang ZJ, Hua QZ et al. The WRKY transcription factor HpWRKY44 regulates CytP450-like1 expression in red pitaya fruit (Hylocereus polyrhizus) . Hortic Res. 2017; 4: 17039.

[27] Kato N, Dubouzet E, Kokabu Y et al. Identification of a WRKY protein as a transcriptional regulator of benzylisoquinoline alkaloid biosynthesis in Coptis japonica . Plant Cell Physiol. 2007; 48: 8-18.

[28] Wang C, Wu C, Wang Y et al. Transcription factor OpWRKY3 is involved in the development and biosynthesis of camptothecin and its precursors in Ophiorrhiza pumila hairy roots. Int J Mol Sci. 2019; 20: 3996.

[29] Xu M, Wu C, Zhao L et al. WRKY transcription factor OpWRKY1 acts as a negative regulator of camptothecin biosynthesis in Ophiorrhiza pumila hairy roots. Plant Cell Tissue Organ Cult. 2020; 142: 69-78.

[30] Ross CA, Liu Y, Shen QJ . The WRKY gene family in rice (Oryza sativa) . J Integr Plant Biol. 2007; 49: 827-42.

[31] Schmutz J, Cannon SB, Schlueter J et al. Genome sequence of the palaeopolyploid soybean. Nature. 2010; 463: 178-83.

[32] He HS, Dong Q, Shao Y et al. Genome-wide survey and characterization of the WRKY gene family in Populus trichocarpa . Plant Cell Rep. 2012; 31: 1199-217.

[33] Wang M, Vannozzi A, Wang G et al. Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family . Hortic Res. 2014; 1: 14016.

[34] Huang X, Li K, Xu X et al. Genome-wide analysis of WRKY transcription factors in white pear (Pyrus bretschneideri) reveals evolution and patterns under drought stress . BMC Genomics. 2015; 16: 1104.

[35] Zhang J, Yang E, He Q et al. Genome-wide analysis of the WRKY gene family in drumstick (Moringa oleifera lam.) . Peer J. 2019; 7: e7063.

[36] Yan H, Li M, Xiong Y et al. Genome-wide characterization, expression profile analysis of WRKY family genes in Santalum album and functional identification of their role in abiotic stress . Int J Mol Sci. 2019; 20: 5676.

[37] Huang SX, Gao Y, Liu J et al. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum . Mol Gen Genomics. 2012; 287: 495-513.

[38] Li D, Liu P, Yu J et al. Genome-wide analysis of WRKY gene family in the sesame genome and identification of the WRKY genes involved in responses to abiotic stresses. BMC Plant Biol. 2017; 17: 152.

[39] Chen CH, Chen X, Han J et al. Genome-wide analysis of the WRKY gene family in the cucumber genome and transcriptome-wide identification of WRKY transcription factors that respond to biotic and abiotic stresses. BMC Plant Biol. 2020; 20: 443.

[40] Zou Z, Yang L, Wang D et al. Gene structures, evolution and transcriptional profiling of the wrky gene family in castor bean (Ricinus communis L.) . PLoS One. 2016; 11: e0148243.

[41] Dong XS, Yang Y, Zhang Z et al. Genome-wide identification of WRKY genes and their response to cold stress in Coffea canephora . Forests. 2019; 10: 355.

[42] Mangelsen E, Kilian J, Berendzen KW et al. Phylogenetic and comparative gene expression analysis of barley (Hordeum vulgare) WRKY transcription factor family reveals putatively retained functions between monocots and dicots . BMC Genomics. 2008; 9: 194.

[43] Xie T, Chen C, Li C et al. Genome-wide investigation of WRKY gene family in pineapple: evolution and expression profiles during development and stress. BMC Genomics. 2018; 19: 490.

[44] Ding ZJ, Yan JY, Li GX et al. WRKY41 controls Arabidopsis seed dormancy via direct regulation of ABI3 transcript levels not downstream of ABA. Plant J. 2014; 79: 810-23.

[45] Chen J, Nolan TM, Ye H et al. Arabidopsis WRKY46, WRKY54, and WRKY70 transcription factors are involved in brassinosteroid-regulated plant growth and drought responses. Plant Cell. 2017; 29: 1425-39.

[46] Fu X, Peng B, Hassani D et al. AaWRKY9 contributes to light- and jasmonate-mediated to regulate the biosynthesis of artemisinin in Artemisia annua. New Phytol. 2021; 231: 1858-74.

[47] Spyropoulou EA, Haring MA et al. RNA sequencing on Solanum lycopersicum trichomes identifies transcription factors that activate terpene synthase promoters. BMC Genomics. 2014; 15: 402.

[48] Wilkins MR, Gasteiger E, Bairoch A et al. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999; 112: 531-52.

[49] Cheng Y, Yao ZP, Ruan MY et al. In silico identification and characterization of the WRKY gene superfamily in pepper (Capsicum annuum L.) . Genet Mol Res. 2016; 15: gmr.15038675.

[50] Kumar S, Stecher G, Tamura K . MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 33: 1870-4.

[51] Hu B, Jin J, Guo AY et al. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015; 31: 1296-7.

[52] Bailey TL, Boden M, Buske FA et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009; 37: W202-8.

[53] Lescot M, Dehais P, Thijs G et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002; 30: 325-7.

[54] Chow CN, Lee TY, Huang YC et al. PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants. Nucleic Acids Res. 2019; 47: D1155-63.

[55] Wang C, Hao X, Wang Y et al. Genome-wide identification and comparative analysis of the teosinte branched 1/cycloidea/proliferating cell factors 1/2 transcription factors related to anti-cancer drug camptothecin biosynthesis in Ophiorrhiza pumila . Front Plant Sci. 2021; 12: 746648.

[56] Shannon P, Markiel A, Ozier O et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13: 2498-504.

[57] Deng CP, Hao X, Shi M et al. Tanshinone production could be increased by the expression of SmWRKY2 in salvia miltiorrhiza hairy roots. Plant Sci. 2019; 284: 1-8.