An automated workflow that generates atom mappings for large-scale metabolic models and its application to Arabidopsis thaliana.

Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large-scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large-scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana. We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)-containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N-MFA studies of A. thaliana. The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph-theoretic structural analyses with large-scale metabolic networks.

Metabolic engineering of the anthocyanin biosynthetic pathway in Artemisia annua and relation to the expression of the artemisinin biosynthetic pathway.

MAIN CONCLUSION: Four types of cells were engineered from Artemisia annua to produce approximately 17 anthocyanins, four of which were elucidated structurally. All of them expressed the artemisinin pathway. Artemisia annua is the only medicinal crop to produce artemisinin for the treatment of malignant malaria. Unfortunately, hundreds of thousands of people still lose their life every year due to the lack of sufficient artemisinin. Artemisinin is considered to result from the spontaneous autoxidation of dihydroartemisinic acid in the presence of reactive oxygen species (ROS) in an oxidative condition of glandular trichomes (GTs); however, whether increasing antioxidative compounds can inhibit artemisinin biosynthesis in plant cells is unknown. Anthocyanins are potent antioxidants that can remove ROS in plant cells. To date, no anthocyanins have been structurally elucidated from A. annua. In this study, we had two goals: (1) to engineer anthocyanins in A. annua cells and (2) to understand the artemisinin biosynthesis in anthocyanin-producing cells. Arabidopsis Production of Anthocyanin Pigment 1 was used to engineer four types of transgenic anthocyanin-producing A. annua (TAPA1-4) cells. Three wild-type cell types were developed as controls. TAPA1 cells produced the highest contents of total anthocyanins. LC-MS analysis detected 17 anthocyanin or anthocyanidin compounds. Crystallization, LC/MS/MS, and NMR analyses identified cyanidin, pelargonidin, one cyanin, and one pelargonin. An integrative analysis characterized that four types of TAPA cells expressed the artemisinin pathway and TAPA1 cells produced the highest artemisinin and artemisinic acid. The contents of arteannuin B were similar in seven cell types. These data showed that the engineering of anthocyanins does not eliminate the biosynthesis of artemisinin in cells. These data allow us to propose a new hypothesis that enzymes catalyze the formation of artemisinin from dihydroartemisinic acid in non-GT cells. These findings show a new platform to increase artemisinin production via non-GT cells of A. annua.

Artemisinin Biosynthesis in Non-glandular Trichome Cells of Artemisia annua.

Artemisinin-based combination therapy (ACT) forms the first line of malaria treatment. However, the yield fluctuation of artemisinin has remained an unsolved problem in meeting the global demand for ACT. This problem is mainly caused by the glandular trichome (GT)-specific biosynthesis of artemisinin in all currently used Artemisia annua cultivars. Here, we report that non-GT cells of self-pollinated inbred A. annua plants can express the artemisinin biosynthetic pathway. Gene expression analysis demonstrated the transcription of six known pathway genes in GT-free leaves and calli of inbred A. annua plants. LC-qTOF-MS/MS analysis showed that these two types of GT-free materials produce artemisinin, artemisinic acid, and arteannuin B. Detailed IR-MALDESI image profiling revealed that these three metabolites and dihydroartemisinin are localized in non-GT cells of leaves of inbred A. annua plants. Moreover, we employed all the above approaches to examine artemisinin biosynthesis in the reported A. annua glandless (gl) mutant. The resulting data demonstrated that leaves of regenerated gl plantlets biosynthesize artemisinin. Collectively, these findings not only add new knowledge leading to a revision of the current dogma of artemisinin biosynthesis in A. annua but also may expedite innovation of novel metabolic engineering approaches for high and stable production of artemisinin in the future.

Overexpression of Artemisia annua Cinnamyl Alcohol Dehydrogenase Increases Lignin and Coumarin and Reduces Artemisinin and Other Sesquiterpenes.

Artemisia annua is the only medicinal crop that produces artemisinin for malarial treatment. Herein, we describe the cloning of a cinnamyl alcohol dehydrogenase (AaCAD) from an inbred self-pollinating (SP) A. annua cultivar and its effects on lignin and artemisinin production. A recombinant AaCAD was purified via heterogeneous expression. Enzyme assays showed that the recombinant AaCAD converted p-coumaryl, coniferyl, and sinapyl aldehydes to their corresponding alcohols, which are key intermediates involved in the biosynthesis of lignin. Km, Vmax, and Vmax/Km values were calculated for all three substrates. To characterize its function in planta, AaCAD was overexpressed in SP plants. Quantification using acetyl bromide (AcBr) showed significantly higher lignin contents in transgenics compared with wild-type (WT) plants. Moreover, GC-MS-based profiling revealed a significant increase in coumarin contents in transgenic plants. By contrast, HPLC-MS analysis showed significantly reduced artemisinin contents in transgenics compared with WT plants. Furthermore, GC-MS analysis revealed a decrease in the contents of arteannuin B and six other sesquiterpenes in transgenic plants. Confocal microscopy analysis showed the cytosolic localization of AaCAD. These data demonstrate that AaCAD plays a dual pathway function in the cytosol, in which it positively enhances lignin formation but negatively controls artemisinin formation. Based on these data, crosstalk between these two pathways mediated by AaCAD catalysis is discussed to understand the metabolic control of artemisinin biosynthesis in plants for high production.

Evaluation of Digital Image Recognition Methods for Mass Spectrometry Imaging Data Analysis.

Analyzing mass spectrometry imaging data can be laborious and time consuming, and as the size and complexity of datasets grow, so does the need for robust automated processing methods. We here present a method for comprehensive, semi-targeted discovery of molecular distributions of interest from mass spectrometry imaging data, using widely available image similarity scoring algorithms to rank images by spatial correlation. A fast and powerful batch search method using a MATLAB implementation of structural similarity (SSIM) index scoring with a pre-selected reference distribution is demonstrated for two sample imaging datasets, a plant metabolite study using Artemisia annua leaf, and a drug distribution study using maraviroc-dosed macaque tissue. Graphical Abstract ᅟ.