Comparative analysis of transcription factor gene families from Papaver somniferum: identification of regulatory factors involved in benzylisoquinoline alkaloid biosynthesis

Promising approaches for simultaneous enhancement of medicinally significant benzylisoquinoline alkaloids in opium poppy

Benzylisoquinoline alkaloids (BIAs) produced in opium poppy have been evidenced to heal patients suffering from various diseases. They, therefore, hold an integral position in the herbal drug industry. Despite the adoption of several approaches for the large-scale production of BIAs, opium poppy remains the only platform in this purpose. The only disadvantage associated with producing BIAs in the plant is their small quantity. Thus, recruiting strategies that boost their levels is deemed necessary. All the methods which have been employed so far are just able to enhance a maximum of two BIAs. Thus, if these methods are utilized, a sizable amount of time and budget must be spent on the synthesis of all BIAs. Hence, the exploitation of strategies which increase the content of all BIAs at the same time is more commercially effective and time-saving, avoiding the laborious step of resolving the biosynthetic pathway of each compound. Exposure to biotic and abiotic elicitors, development of a synthetic auto-tetraploid, overexpression of a WRKY transcription factor, formation of an artificial metabolon, and suppression of a gene in the shikimate pathway and miRNA are strategies that turn opium poppy into a versatile bioreactor for the concurrent and massive production of BIAs. The last three strategies have never been applied for BIA biosynthetic pathways.

Transcriptome and Metabolome Analysis of Isoquinoline Alkaloid Biosynthesis of Coptis chinensis in Different Years

Coptis chinensis is a perennial herb of the Ranunculaceae family. The isoquinoline alkaloid is the main active component of C. chinensis, mainly exists in its rhizomes and has high clinical application potential. The in vitro synthesis of isoquinoline alkaloids is difficult because their structures are complex; hence, plants are still the main source of them. In this study, two-year and four-year rhizomes of C. chinensis were selected to investigate the effect of growth years on the accumulation of isoquinoline alkaloids. Two-year and four-year C. chinensis were selected for metabolomics detection and transcriptomic analysis. A total of 413 alkaloids were detected by metabolomics analysis, of which 92 were isoquinoline alkaloids. (S)-reticuline was a significantly different accumulated metabolite of the isoquinoline alkaloids biosynthetic pathway in C. chinensis between the two groups. The results of transcriptome analysis showed that a total of 464 differential genes were identified, 36 of which were associated with the isoquinoline alkaloid biosynthesis pathway of C. chinensis. Among them, 18 genes were correlated with the content of important isoquinoline alkaloids. Overall, this study provided a comprehensive metabolomic and transcriptomic analysis of the rapid growth stage of C. chinensis rhizome from the perspective of growth years. It brought new insights into the biosynthetic pathway of isoquinoline alkaloids and provided information for utilizing biotechnology to improve their contents in C. chinensis.

Dissection of transcriptome and metabolome insights into the isoquinoline alkaloid biosynthesis during stem development in Phellodendron amurense (Rupr.)

Phellodendron amurense (Rupr.) is a well-known medicinal plant with high medicinal value, and its various tissues are enriched in various active pharmaceutical ingredients. Isoquinoline alkaloids are the primary medicinal component of P. amurense and have multiple effects, such as anti-inflammation, antihypertension, and antitumor effects. However, the potential regulatory mechanism of isoquinoline alkaloid biosynthesis during stem development in P. amurense is still poorly understood. In the present study, a total of eight plant hormones for each stem development stage were detected; of those, auxin, gibberellins and brassinosteroids were significantly highly increased in perennial stems and played key roles during stem development in P. amurense. We also investigated the content and change pattern of secondary metabolites and comprehensively identified some key structural genes involved in the isoquinoline alkaloid biosynthesis pathway by combining the transcriptome and metabolomics. A total of 39,978 DEGs were identified in the present study, and six of those had candidate structural genes (NCS, GOT2, TYNA, CODM, TYR, TAT and PSOMT1) that were specifically related to isoquinoline alkaloid biosynthesis in P. amurense. Corydalmine, cyclanoline, dehydroyanhunine, (S)-canadine and corybulbine were the most significantly upregulated metabolites among the different comparative groups. Three differentially expressed metabolites, dopamine, (S)-corytuberine and (S)-canadine, were enriched in the isoquinoline alkaloid biosynthesis pathway. Furthermore, bHLH and WRKY transcription factors play key roles in the isoquinoline alkaloid biosynthesis pathway in P. amurense. The results not only provide comprehensive genetic information for understanding the molecular mechanisms of isoquinoline alkaloid biosynthesis but also lay a foundation for the combinatory usage of the medicinal active ingredient of P. amurense.

Interplay of transcription factors orchestrating the biosynthesis of plant alkaloids

Plants produce a range of secondary metabolites primarily as defence molecules. A plant has to invest considerable energy to synthesise alkaloids, and sometimes they are even toxic to themselves. Hence, the biosynthesis of alkaloids is a spatiotemporally regulated process under quantitative feedback regulation which is accomplished by the signal reception, transcriptional/translational regulation, transport, storage and accumulation. The transcription factors (TFs) initiate the biosynthesis of alkaloids after appropriate cues. The present study recapitulates last decade understanding of the role of TFs in alkaloid biosynthesis. The present review discusses TF families, viz. AP2/ERF, bHLH, WRKY, MYB involved in the biosynthesis of various types of alkaloids. It also highlights the role of the jasmonic acid cascade and post-translational modifications of TF proteins. A thorough understanding of TFs will help us to decide a strategy to facilitate successful pathway manipulation and in vitro production.

Comparative transcriptome analysis of tubers, stems, and flowers of Gastrodia elata Blume reveals potential genes involved in the biosynthesis of phenolics

Orchidaceae, well known for its fascinating flowers, is one of the largest and most diverse families of flowering plants. There are many kinds of plants in this family; these are distributed practically globally and have high ornamental and medicinal values. Gastrodia elata Blume, a traditional Chinese medicinal herb, is a rootless and leafless achlorophyllous orchid. Phenolic compounds are considered to be the major bioactive constituents in G. elata, with antioxidant, antiangiogenic, neuroprotective, antidepressant, anxiolytic, and sedative activities. In this study, we determined the contents of six main phenolic components in tubers, stems and flowers from G. elata. Meanwhile, the transcriptomes of the tuber, stem and flower tissues of G. elata were obtained using the BGISEQ-500 platform. A total of 58.29 Gb of data and 113,067 unigenes were obtained, of which 74,820 unigenes were functionally annotated against seven public databases. Differentially expressed genes between tuber, stem and flower tissues were identified. A total of 76 DEGs encoding eight key enzymes were identified as candidate genes involved in the biosynthesis of phenolics in G. elata. For further validation, the expression levels of unigenes were measured using quantitative real-time PCR. Our results greatly enrich the transcriptomic data of G. elata and provide valuable information for the identification of candidate genes involved in the biosynthesis of secondary metabolites.

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.

Full-length transcriptome analysis and identification of genes involved in asarinin and aristolochic acid biosynthesis in medicinal plant Asarum sieboldii

Asarum sieboldii, a well-known traditional Chinese medicinal herb, is used for curing inflammation and ache. It contains both the bioactive ingredient asarinin and the toxic compound aristolochic acid. To address further breeding demand, genes involved in the biosynthetic pathways of asarinin and aristolochic acid should be explored. Therefore, the full-length transcriptome of A. sieboldii was sequenced using PacBio Iso-Seq to determine the candidate transcripts that encode the biosynthetic enzymes of asarinin and aristolochic acid. In this study, 63 023 full-length transcripts were generated with an average length of 1371 bp from roots, stems, and leaves, of which 49 593 transcripts (78.69%) were annotated against public databases. Furthermore, 555 alternative splicing (AS) events, 10 869 long noncoding RNAs (lncRNAs) as well as their 11 291 target genes, and 17 909 simple sequence repeats (SSRs) were identified. The data also revealed 97 candidate transcripts related to asarinin metabolism, of which six novel genes that encoded enzymes involved in asarinin biosynthesis were initially reported. In addition, 56 transcripts related to aristolochic acid biosynthesis were also identified, including CYP81B. In summary, these transcriptome data provide a useful resource to study gene function and genetic engineering in A. sieboldii.

Interpol review of controlled substances 2016-2019

This review paper covers the forensic-relevant literature in controlled substances from 2016 to 2019 as a part of the 19th Interpol International Forensic Science Managers Symposium. The review papers are also available at the Interpol website at: https://www.interpol.int/content/download/14458/file/Interpol%20Review%20Papers%202019.pdf.

Mining lycodine-type alkaloid biosynthetic genes and genetic markers from transcriptome of Lycopodiastrum casuarinoides

Objective

Lycopodiastrum casuarinoides, a fern of the Lycopodiaceae family, is a traditional Chinese medicine, which has similar efficacy to that of Huperzia serrata in treating rheumatoid arthritis (RA). However, they are different in the contents and compositions of lycopodium alkaloids. In this study, the biosynthesis related genes of lycopodium alkaloids and genetic markers are discovered in L. casuarinoides transcriptome.

Methods

The plant of L. casuarinoides was collected and was subjected to the RNA isolation, cDNA library construction, high throughput RNA sequencing and bioinformatics analysis.

Results

Totally 124, 524 high-quality unigenes were assembled from RNA sequencing reads, with an average sequence length of 601 bp. Among the L. casuarinoides transcripts, 61,304 shared the significant similarity (E-value < 10⁻⁵) with existing protein sequences in the public databases. From 124,524 unigenes, 47,538 open reading frames (ORFs) were predicted. Based on the bioinformatics analysis, all possible enzyme genes involved in the lycodine-type alkaloids biosynthetic pathway of L. casuarinoides were identified, including lysine decarboxylase (LDC), primary amine oxidase (PAO), malonyl-CoA decarboxylase, etc. Sixty-four putative cytochrome p450 (CYP) and 827 putative transcription factors were selected from the transcriptome unigenes as the candidates of lycodine-type alkaloids biosynthesis modifiers. Furthermore, 13,352 simple sequence repeats (SSRs) were identified from 124,524 unigenes, of which dinucleotide motifs AG/CT were the most abundant (50.1%). Meanwhile, we confirmed the amplification effectiveness of 25 PCR primer pairs for randomly selected SSRs.

Conclusion

We obtained the comprehensive transcriptomic information from the high throughput RNA sequencing and bioinformatics analysis, which provided a valuable resource of transcript sequences of L. casuarinoides in public databases.

3'O-Methyltransferase, Ps3'OMT, from opium poppy: involvement in papaverine biosynthesis

Using, in silico, in vitro and in planta functional assays, we demonstrate that Ps3'OMT, an 3'-O methyl transferase is linked to papaverine biosynthesis in opium poppy. Papaverine, one of the benzylisoquinoline alkaloids (BIA) synthesized in the medicinally important plant, Papaver somniferum, is known for the potent pharmacological properties. Papaverine biosynthesis has remained debatable as two different pathways, NH (involving N-desmethylated intermediates) and the NCH₃ (involving N-methylated intermediates), have been proposed. In addition, there are several intermediate steps in both the proposed pathways that are not very well characterized in terms of specific enzymes. In this study, we report the identification and functional characterization of 3'O-methyltransferase (Ps3'OMT) which might participate in the 3'O-methylation of the intermediates in the papaverine biosynthesis. Comparison of transcript and metabolite profiles of high and low papaverine producing cultivar revealed the occurrence of a 3'O-methyltransferase, Ps3'OMT, which was abundant in aerial organs and shared 72% identity with the GfLOMT7 predicted to have 3'OMT activity. In silico studies based on homology modeling, docking and MD simulations predicted (S)-norlaudanine as the potential substrate forming a stable complex with Ps3'OMT. Suppression of Ps3'OMT through virus-induced gene silencing resulted in a remarkable decrease in the level of papaverine in comparison to control plants. The characterization of the functionally unique Ps3'OMT involved in BIA metabolism suggests an involvement of the NH pathway leading to papaverine biosynthesis.

Identification of candidate genes involved in isoquinoline alkaloids biosynthesis in Dactylicapnos scandens by transcriptome analysis

Dactylicapnos scandens (D. Don) Hutch (Papaveraceae) is a well-known traditional Chinese herb used for treatment of hypertension, inflammation, bleeding and pain for centuries. Although the major bioactive components in this herb are considered as isoquinoline alkaloids (IQAs), little is known about molecular basis of their biosynthesis. Here, we carried out transcriptomic analysis of roots, leaves and stems of D. scandens, and obtained a total of 96,741 unigenes. Based on gene expression and phylogenetic relationship, we proposed the biosynthetic pathways of isocorydine, corydine, glaucine and sinomenine, and identified 67 unigenes encoding enzymes potentially involved in biosynthesis of IQAs in D. scandens. High performance liquid chromatography analysis demonstrated that while isocorydine is the most abundant IQA in D. scandens, the last O-methylation biosynthesis step remains unclear. Further enzyme activity assay, for the first time, characterized a gene encoding O- methyltransferase (DsOMT), which catalyzes O-methylation at C7 of (S)-corytuberine to form isocorydine. We also identified candidate transcription factor genes belonging to WRKY and bHLH families that may be involved in the regulation of IQAs biosynthesis. Taken together, we first provided valuable genetic information for D. scandens, shedding light on candidate genes involved in IQA biosynthesis, which will be critical for further gene functional characterization.

Genome-Wide Analysis of the Musa WRKY Gene Family: Evolution and Differential Expression during Development and Stress

The WRKY gene family plays an important role in the development and stress responses in plants. As information is not available on the WRKY gene family in Musa species, genome-wide analysis has been carried out in this study using available genomic information from two species, Musa acuminata and Musa balbisiana. Analysis identified 147 and 132 members of the WRKY gene family in M. acuminata and M. balbisiana, respectively. Evolutionary analysis suggests that the WRKY gene family expanded much before the speciation in both the species. Most of the orthologs retained in two species were from the γ duplication event which occurred prior to α and β genome-wide duplication (GWD) events. Analysis also suggests that subtle changes in nucleotide sequences during the course of evolution have led to the development of new motifs which might be involved in neo-functionalization of different WRKY members in two species. Expression and cis-regulatory motif analysis suggest possible involvement of Group II and Group III WRKY members during various stresses and growth/development including fruit ripening process respectively.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HDZ IV) Gene Family from Musa accuminata

The homeodomain zipper family (HD-ZIP) of transcription factors is present only in plants and plays important role in the regulation of plant-specific processes. The subfamily IV of HDZ transcription factors (HD-ZIP IV) has primarily been implicated in the regulation of epidermal structure development. Though this gene family is present in all lineages of land plants, members of this gene family have not been identified in banana, which is one of the major staple fruit crops. In the present work, we identified 21 HDZIV encoding genes in banana by the computational analysis of banana genome resource. Our analysis suggested that these genes putatively encode proteins having all the characteristic domains of HDZIV transcription factors. The phylogenetic analysis of the banana HDZIV family genes further confirmed that after separation from a common ancestor, the banana, and poales lineages might have followed distinct evolutionary paths. Further, we conclude that segmental duplication played a major role in the evolution of banana HDZIV encoding genes. All the identified banana HDZIV genes expresses in different banana tissue, however at varying levels. The transcript levels of some of the banana HDZIV genes were also detected in banana fruit pulp, suggesting their putative role in fruit attributes. A large number of genes of this family showed modulated expression under drought and salinity stress. Taken together, the present work lays a foundation for elucidation of functional aspects of the banana HDZIV encoding genes and for their possible use in the banana improvement programs.

Comparative transcriptome analysis of different chemotypes elucidates withanolide biosynthesis pathway from medicinal plant Withania somnifera

Withania somnifera is one of the most valuable medicinal plants synthesizing secondary metabolites known as withanolides. Despite pharmaceutical importance, limited information is available about the biosynthesis of withanolides. Chemo-profiling of leaf and root tissues of Withania suggest differences in the content and/or nature of withanolides in different chemotypes. To identify genes involved in chemotype and/or tissue-specific withanolide biosynthesis, we established transcriptomes of leaf and root tissues of distinct chemotypes. Genes encoding enzymes for intermediate steps of terpenoid backbone biosynthesis with their alternatively spliced forms and paralogous have been identified. Analysis suggests differential expression of large number genes among leaf and root tissues of different chemotypes. Study also identified differentially expressing transcripts encoding cytochrome P450s, glycosyltransferases, methyltransferases and transcription factors which might be involved in chemodiversity in Withania. Virus induced gene silencing of the sterol ∆7-reductase (WsDWF5) involved in the synthesis of 24-methylene cholesterol, withanolide backbone, suggests role of this enzyme in biosynthesis of withanolides. Information generated, in this study, provides a rich resource for functional analysis of withanolide-specific genes to elucidate chemotype- as well as tissue-specific withanolide biosynthesis. This genomic resource will also help in development of new tools for functional genomics and breeding in Withania.