Transcriptional Profiles of Secondary Metabolite Biosynthesis Genes and Cytochromes in the Leaves of Four Papaver Species

Poppies are well-known plants in the family Papaveraceae that are rich in alkaloids. This family contains 61 species, and in this study we sequenced the transcriptomes of four species’ (Papaver rhoeas, Papaver nudicaule, Papaver fauriei, and Papaver somniferum) leaves. These transcripts were systematically assessed for the expression of secondary metabolite biosynthesis (SMB) genes and cytochromes, and their expression profiles were assessed for use in bioinformatics analyses. This study contributed 265 Gb (13 libraries with three biological replicates) of leaf transcriptome data from three Papaver plant developmental stages. Sequenced transcripts were assembled into 815 Mb of contigs, including 226 Mb of full-length transcripts. The transcripts for 53 KEGG pathways, 55 cytochrome superfamilies, and benzylisoquinoline alkaloid biosynthesis (BIA) were identified and compared to four other alkaloid-rich genomes. Additionally, 22 different alkaloids and their relative expression profiles in three developmental stages of Papaver species were assessed by targeted metabolomics using LC-QTOF-MS/MS. Collectively, the results are given in co-occurrence heat-maps to help researchers obtain an overview of the transcripts and their differential expression in the Papaver development life cycle, particularly in leaves. Moreover, this dataset will be a valuable resource to derive hypotheses to mitigate an array of Papaver developmental and secondary metabolite biosynthesis issues in the future.

[1]  Colin N. Dewey,et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis , 2013, Nature Protocols.

[2]  Seok-Geun Lee,et al.  Transcriptome Profiling of Two Ornamental and Medicinal Papaver Herbs , 2018, International journal of molecular sciences.

[3]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[4]  C. Veeresham Natural products derived from plants as a source of drugs , 2012, Journal of advanced pharmaceutical technology & research.

[5]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[6]  Yuxian Zhu,et al.  Cotton GhPOX1 encoding plant class III peroxidase may be responsible for the high level of reactive oxygen species production that is related to cotton fiber elongation. , 2009, Journal of genetics and genomics = Yi chuan xue bao.

[7]  Heather M. Wallace,et al.  Update of the Scientific Opinion on opium alkaloids in poppy seeds , 2018, EFSA journal. European Food Safety Authority.

[8]  Renate A Weizbauer,et al.  Role of aromatic aldehyde synthase in wounding/herbivory response and flower scent production in different Arabidopsis ecotypes. , 2011, The Plant journal : for cell and molecular biology.

[9]  D. Nelson,et al.  Diversity and Evolution of Plant P450 and P450-Reductases , 1995, Drug metabolism and drug interactions.

[10]  Yanran Li,et al.  Engineering biosynthesis of the anticancer alkaloid noscapine in yeast , 2016, Nature Communications.

[11]  B. Møller,et al.  CYP703 Is an Ancient Cytochrome P450 in Land Plants Catalyzing in-Chain Hydroxylation of Lauric Acid to Provide Building Blocks for Sporopollenin Synthesis in Pollen[W] , 2007, The Plant Cell Online.

[12]  Turgay Unver,et al.  Regulation of the alkaloid biosynthesis by miRNA in opium poppy. , 2015, Plant biotechnology journal.

[13]  S. Subramaniyam,et al.  Transcript expression profiling for adventitious roots of Panax ginseng Meyer. , 2014, Gene.

[14]  Dustin J Cram,et al.  Integration of deep transcript and targeted metabolite profiles for eight cultivars of opium poppy , 2012, Plant Molecular Biology.

[15]  C. Sensen,et al.  Transcriptome analysis of 20 taxonomically related benzylisoquinoline alkaloid-producing plants , 2015, BMC Plant Biology.

[16]  Michael Ott,et al.  De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with Trinity , 2013 .

[17]  Gerard Talavera,et al.  Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. , 2007, Systematic biology.

[18]  C. Stoeckert,et al.  OrthoMCL: identification of ortholog groups for eukaryotic genomes. , 2003, Genome research.

[19]  C. Fu,et al.  Transcriptome Assembly and Systematic Identification of Novel Cytochrome P450s in Taxus chinensis , 2017, Front. Plant Sci..

[20]  Z. Lorigooini,et al.  Medicinal plants: Past history and future perspective , 2017 .

[21]  D. Barceloux Heroin and the Opium Poppy Plant (Papaver somniferum L.) , 2012 .

[22]  Yasuyuki Yamada,et al.  Mining of the Uncharacterized Cytochrome P450 Genes Involved in Alkaloid Biosynthesis in California Poppy Using a Draft Genome Sequence , 2017, Plant & cell physiology.

[23]  R. Zhong,et al.  Essential role of caffeoyl coenzyme A O-methyltransferase in lignin biosynthesis in woody poplar plants. , 2000, Plant physiology.

[24]  Sònia Garcia,et al.  Recent updates and developments to plant genome size databases , 2013, Nucleic Acids Res..

[25]  D. Ro,et al.  The Genome of Medicinal Plant Macleaya cordata Provides New Insights into Benzylisoquinoline Alkaloids Metabolism. , 2017, Molecular plant.

[26]  Yi Li,et al.  Morphinan biosynthesis in opium poppy requires a P450-oxidoreductase fusion protein , 2015, Science.

[27]  D. Nelson,et al.  A P450-centric view of plant evolution. , 2011, The Plant journal : for cell and molecular biology.

[28]  C. Smolke,et al.  Complete biosynthesis of opioids in yeast , 2015, Science.

[29]  Ute Roessner,et al.  A Historical Overview of Natural Products in Drug Discovery , 2012, Metabolites.

[30]  Shanfa Lu,et al.  Biosynthesis and Regulation of Phenylpropanoids in Plants , 2017 .

[31]  J. Pleiss,et al.  Identification of universal selectivity‐determining positions in cytochrome P450 monooxygenases by systematic sequence‐based literature mining , 2015, Proteins.

[32]  Ramya Mathiyalagan,et al.  Ginseng nanoparticles: a budding tool for cancer treatment. , 2017, Nanomedicine.

[33]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[34]  Kazutaka Katoh,et al.  Application of the MAFFT sequence alignment program to large data—reexamination of the usefulness of chained guide trees , 2016, Bioinform..

[35]  Jeremy S. Morris,et al.  Family portraits: the enzymes behind benzylisoquinoline alkaloid diversity , 2017, Phytochemistry Reviews.

[36]  Saravanaraj N. Ayyampalayam,et al.  Phylotranscriptomic analysis of the origin and early diversification of land plants , 2014, Proceedings of the National Academy of Sciences.

[37]  R. Al-Hasani,et al.  Pain and Poppies: The Good, the Bad, and the Ugly of Opioid Analgesics , 2015, The Journal of Neuroscience.

[38]  Stefan Götz,et al.  Blast2GO: A Comprehensive Suite for Functional Analysis in Plant Genomics , 2007, International journal of plant genomics.

[39]  Jinmi Yoon,et al.  Roles of lignin biosynthesis and regulatory genes in plant development , 2015, Journal of integrative plant biology.

[40]  D. Schriemer,et al.  Integration of deep transcriptome and proteome analyses reveals the components of alkaloid metabolism in opium poppy cell cultures , 2010, BMC Plant Biology.

[41]  Sean R. Eddy,et al.  Accelerated Profile HMM Searches , 2011, PLoS Comput. Biol..

[42]  Minoru Kanehisa,et al.  KEGG Bioinformatics Resource for Plant Genomics and Metabolomics. , 2016, Methods in molecular biology.

[43]  Sumya Pathak,et al.  Comparative Transcriptome Analysis Using High Papaverine Mutant of Papaver somniferum Reveals Pathway and Uncharacterized Steps of Papaverine Biosynthesis , 2013, PloS one.

[44]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[45]  Rob Patro,et al.  Salmon provides fast and bias-aware quantification of transcript expression , 2017, Nature Methods.

[46]  I. Sharma,et al.  Versatile roles of brassinosteroid in plants in the context of its homoeostasis, signaling and crosstalks , 2015, Front. Plant Sci..

[47]  Sungyoung Lee,et al.  Classification and characterization of putative cytochrome P450 genes from Panax ginseng C. A. Meyer , 2011, Biochemistry (Moscow).

[48]  J. Bernáth Poppy: The Genus Papaver , 1999 .

[49]  Yi Li,et al.  The opium poppy genome and morphinan production , 2018, Science.

[50]  Y. Li-Beisson,et al.  Fatty Acid and Lipid Transport in Plant Cells. , 2016, Trends in plant science.

[51]  U. C. Lavania,et al.  Evolutionary DNA variation in Papaver , 1991 .

[52]  D. Nelson Cytochrome P450 diversity in the tree of life. , 2018, Biochimica et biophysica acta. Proteins and proteomics.

[53]  Robert M. Waterhouse,et al.  BUSCO Applications from Quality Assessments to Gene Prediction and Phylogenomics , 2017, bioRxiv.

[54]  F. Beisson,et al.  Cytochrome P450 metabolizing fatty acids in plants: characterization and physiological roles , 2011, The FEBS journal.

[55]  J. Gershenzon,et al.  Herbivore-induced poplar cytochrome P450 enzymes of the CYP71 family convert aldoximes to nitriles which repel a generalist caterpillar. , 2014, The Plant journal : for cell and molecular biology.

[56]  J. Thevelein,et al.  Bioengineering of plant (tri)terpenoids: from metabolic engineering of plants to synthetic biology in vivo and in vitro. , 2013, The New phytologist.