The medicinal properties of cannabis and the its legal status in several countries and jurisdictions has spurred the massive growth of the cannabis economy around the globe. The value of cannabis stems from its euphoric activity offered by the unique phytocannabinoid tetrahydrocannabinol (THC). However, this is rapidly expanding beyond THC owing to other non-psychoactive phytocannabinoids with new bioactivities that will contribute to their development into clinically useful drugs. The discovery of the biosynthesis of major phytocannabinoids has allowed the exploration of their heterologous production by synthetic biology, which may lead to the industrial production of rare phytocannabinoids or novel synthetic cannabinoid pharmaceuticals that are not easily offered by cannabis plants. This review summarizes the biosynthesis of major phytocannabinoids in detail, the most recent development of their metabolic engineering in various systems, and the engineering approaches and strategies used to increase the yield.

Integrative taxonomy is a fundamental part of biodiversity and combines traditional morphology with additional methods such as DNA sequencing or biochemistry. Here, we aim to establish untargeted metabolomics for use in chemotaxonomy. We used three thallose liverwort species Riccia glauca, R. sorocarpa, and R. warnstorfii (order Marchantiales, Ricciaceae) with Lunularia cruciata (order Marchantiales, Lunulariacea) as an outgroup. Liquid chromatography high-resolution mass-spectrometry (UPLC/ESI-QTOF-MS) with data-dependent acquisition (DDA-MS) were integrated with DNA marker-based sequencing of the trnL-trnF region and high-resolution bioimaging. Our untargeted chemotaxonomy methodology enables us to distinguish taxa based on chemophenetic markers at different levels of complexity: (1) molecules, (2) compound classes, (3) compound superclasses, and (4) molecular descriptors. For the investigated Riccia species, we identified 71 chemophenetic markers at the molecular level, a characteristic composition in 21 compound classes, and 21 molecular descriptors largely indicating electron state, presence of chemical motifs, and hydrogen bonds. Our untargeted approach revealed many chemophenetic markers at different complexity levels that can provide more mechanistic insight into phylogenetic delimitation of species within a clade than genetic-based methods coupled with traditional morphology-based information. However, analytical and bioinformatics analysis methods still need to be better integrated to link the chemophenetic information at multiple scales.

Bibenzyls are a specialized metabolite class found throughout the plant kingdom. One of the most prolific producers of bibenzyls are liverworts, specifically plants of the Radula genera. These plants possess an incredible diversity of bibenzyls, prenylated bibenzyls, and a few (bis)bibenzyls, several of which have medicinal properties, including perrottetinene, an analog of tetrahydrocannabinol from cannabis. To provide insight into the bibenzyls’ biosynthesis in planta, exogenous phytohormones were applied to in vitro grown Radula complanata and bibenzyl metabolite production was monitored with targeted and untargeted metabolomics. The targeted metabolomic analysis of six prenylated bibenzyls revealed that production of these metabolites was largely reduced when plants were treated with abscisic acid (AA), salicylic acid (SA), 1-naphthaleneacetic acid (NAA), or 6-benzylaminopurine (BAP). The reduction of these metabolites in the BAP and NAA treatment suggests that prenylated bibenzyl production is negatively correlated with vegetative plant growth. The reduction of bibenzyls at low AA and SA concentrations and mild increase at higher AA and SA concentrations suggest that their production is regulated by these stress hormones. In addition, six other bibenzyl metabolites were tentatively identified from the untargeted analysis. These results provide insight into the influence of phytohormones on the bioactive bibenzyl content of R. complanata.

Introduction: Liverworts are a group of non-vascular plants that possess unique metabolism not found in other plants. Many liverwort metabolites have interesting structural and biochemical characteristics, however the fluctuations of these metabolites in response to stressors is largely unknown. Objectives: To investigate the metabolic stress-response of the leafy liverwort Radula complanata. Methods: Five phytohormones were applied exogenously to in vitro cultured R. complanata and an untargeted metabolomic analysis was conducted. Compound classification and identification was performed with CANOPUS and SIRIUS while statistical analyses including PCA, ANOVA, and variable selection using BORUTA were conducted to identify metabolic shifts. Results: It was found that R. complanata was predominantly composed of carboxylic acids and derivatives, followed by benzene and substituted derivatives, fatty acyls, organooxygen compounds, prenol lipids, and flavonoids. The PCA revealed that samples grouped based on the type of hormone applied, and the variable selection using BORUTA (Random Forest) revealed 71 identified and/or classified features that fluctuated with phytohormone application. The stress-response treatments largely reduced the production of the selected primary metabolites while the growth treatments resulted in increased production of these compounds. 4-(3-Methyl-2-butenyl)-5-phenethylbenzene-1,3-diol was identified as a biomarker for the growth treatments while GDP-hexose was identified as a biomarker for the stress-response treatments. Conclusion: Exogenous phytohormone application caused clear metabolic shifts in Radula complanata that deviate from the responses of vascular plants. Further identification of the selected metabolite features can reveal metabolic biomarkers unique to liverworts and provide more insight into liverwort stress responses.

Bryophytes are prolific producers of unique, specialized metabolites that are not found in other plants. As many of these unique natural products are potentially interesting, for example, pharmacological use, variations in the production regarding ecological or environmental conditions have not often been investigated. Here, we investigate metabolic shifts in the epiphytic Radula complanata L. (Dumort) with regard to different environmental conditions and the type of phorophyte (host tree). Plant material was harvested from three different locations in Sweden, Germany, and Canada and subjected to untargeted liquid chromatography high-resolution mass-spectrometry (UPLC/ESI-QTOF-MS) and data-dependent acquisition (DDA-MS). Using multivariate statistics, variable selection methods, in silico compound identification, and compound classification, a large amount of variation (39%) in the metabolite profiles was attributed to the type of host tree and 25% to differences in environmental conditions. We identified 55 compounds to vary significantly depending on the host tree (36 on the family level) and 23 compounds to characterize R. complanata in different environments. Taken together, we found metabolic shifts mainly in primary metabolites that were associated with the drought response to different humidity levels. The metabolic shifts were highly specific to the host tree, including mostly specialized metabolites suggesting high levels of ecological interaction. As R. complanata is a widely distributed generalist species, we found it to flexibly adapt its metabolome according to different conditions. We found metabolic composition to also mirror the constitution of the habitat, which makes it interesting for conservation measures.