Diatoms (Bacillariophyceae) are aquatic photosynthetic microalgae with an ecological role as primary producers in the aquatic food web. They account substantially for global carbon, nitrogen, and silicon cycling. Elucidating the chemical space of diatoms is crucial to understanding their physiology and ecology. To expand the known chemical space of a cosmopolitan marine diatom, Skeletonema marinoi, we performed High-Resolution Liquid Chromatography-Tandem Mass Spectrometry (LC-MS2) for untargeted metabolomics data acquisition. The spectral data from LC-MS2 was used as input for the Metabolome Annotation Workflow (MAW) to obtain putative annotations for all measured features. A suspect list of metabolites previously identified in the Skeletonema spp. was generated to verify the results. These known metabolites were then added to the putative candidate list from LC-MS2 data to represent an expanded catalog of 1970 metabolites estimated to be produced by S. marinoi. The most prevalent chemical superclasses, based on the ChemONT ontology in this expanded dataset, were organic acids and derivatives, organoheterocyclic compounds, lipids and lipid-like molecules, and organic oxygen compounds. The metabolic profile from this study can aid the bioprospecting of marine microalgae for medicine, biofuel production, agriculture, and environmental conservation. The proposed analysis can be applicable for assessing the chemical space of other microalgae, which can also provide molecular insights into the interaction between marine organisms and their role in the functioning of ecosystems.

Mapping the chemical space of compounds to chemical structures remains a challenge in metabolomics. Despite the advancements in untargeted liquid chromatography-mass spectrometry (LC–MS) to achieve a high-throughput profile of metabolites from complex biological resources, only a small fraction of these metabolites can be annotated with confidence. Many novel computational methods and tools have been developed to enable chemical structure annotation to known and unknown compounds such as in silico generated spectra and molecular networking. Here, we present an automated and reproducible Metabolome Annotation Workflow (MAW) for untargeted metabolomics data to further facilitate and automate the complex annotation by combining tandem mass spectrometry (MS2) input data pre-processing, spectral and compound database matching with computational classification, and in silico annotation. MAW takes the LC-MS2 spectra as input and generates a list of putative candidates from spectral and compound databases. The databases are integrated via the R package Spectra and the metabolite annotation tool SIRIUS as part of the R segment of the workflow (MAW-R). The final candidate selection is performed using the cheminformatics tool RDKit in the Python segment (MAW-Py). Furthermore, each feature is assigned a chemical structure and can be imported to a chemical structure similarity network. MAW is following the FAIR (Findable, Accessible, Interoperable, Reusable) principles and has been made available as the docker images, maw-r and maw-py. The source code and documentation are available on GitHub (https://github.com/zmahnoor14/MAW). The performance of MAW is evaluated on two case studies. MAW can improve candidate ranking by integrating spectral databases with annotation tools like SIRIUS which contributes to an efficient candidate selection procedure. The results from MAW are also reproducible and traceable, compliant with the FAIR guidelines. Taken together, MAW could greatly facilitate automated metabolite characterization in diverse fields such as clinical metabolomics and natural product discovery.

The metabolome is the biochemical basis of plant form and function, but we know little about its macroecological variation across the plant kingdom. Here, we used the plant functional trait concept to interpret leaf metabolome variation among 457 tropical and 339 temperate plant species. Distilling metabolite chemistry into five metabolic functional traits reveals that plants vary on two major axes of leaf metabolic specialization—a leaf chemical defense spectrum and an expression of leaf longevity. Axes are similar for tropical and temperate species, with many trait combinations being viable. However, metabolic traits vary orthogonally to life-history strategies described by widely used functional traits. The metabolome thus expands the functional trait concept by providing additional axes of metabolic specialization for examining plant form and function.

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.

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.