Ontologies for increasing the FAIRness of plant research data.

The importance of improving the FAIRness (findability, accessibility, interoperability, reusability) of research data is undeniable, especially in the face of large, complex datasets currently being produced by omics technologies. Facilitating the integration of a dataset with other types of data increases the likelihood of reuse, and the potential of answering novel research questions. Ontologies are a useful tool for semantically tagging datasets as adding relevant metadata increases the understanding of how data was produced and increases its interoperability. Ontologies provide concepts for a particular domain as well as the relationships between concepts. By tagging data with ontology terms, data becomes both human- and machine- interpretable, allowing for increased reuse and interoperability. However, the task of identifying ontologies relevant to a particular research domain or technology is challenging, especially within the diverse realm of fundamental plant research. In this review, we outline the ontologies most relevant to the fundamental plant sciences and how they can be used to annotate data related to plant-specific experiments within metadata frameworks, such as Investigation-Study-Assay (ISA). We also outline repositories and platforms most useful for identifying applicable ontologies or finding ontology terms.

Leaf metabolic traits reveal hidden dimensions of plant form and function.

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.

Impact of in vitro phytohormone treatments on the metabolome of the leafy liverwort Radula complanata (L.) Dumort.

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.

Untargeted Metabolomics for Integrative Taxonomy: Metabolomics, DNA Marker-Based Sequencing, and Phenotype Bioimaging.

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.

Host Tree and Geography Induce Metabolic Shifts in the Epiphytic Liverwort Radula complanata.

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.

Proteomics Standards Initiative at Twenty Years: Current Activities and Future Work.

The Human Proteome Organization (HUPO) Proteomics Standards Initiative (PSI) has been successfully developing guidelines, data formats, and controlled vocabularies (CVs) for the proteomics community and other fields supported by mass spectrometry since its inception 20 years ago. Here we describe the general operation of the PSI, including its leadership, working groups, yearly workshops, and the document process by which proposals are thoroughly and publicly reviewed in order to be ratified as PSI standards. We briefly describe the current state of the many existing PSI standards, some of which remain the same as when originally developed, some of which have undergone subsequent revisions, and some of which have become obsolete. Then the set of proposals currently being developed are described, with an open call to the community for participation in the forging of the next generation of standards. Finally, we describe some synergies and collaborations with other organizations and look to the future in how the PSI will continue to promote the open sharing of data and thus accelerate the progress of the field of proteomics.

The community ecology perspective of omics data.

The measurement of uncharacterized pools of biological molecules through techniques such as metabarcoding, metagenomics, metatranscriptomics, metabolomics, and metaproteomics produces large, multivariate datasets. Analyses of these datasets have successfully been borrowed from community ecology to characterize the molecular diversity of samples (ɑ-diversity) and to assess how these profiles change in response to experimental treatments or across gradients (ß-diversity). However, sample preparation and data collection methods generate biases and noise which confound molecular diversity estimates and require special attention. Here, we examine how technical biases and noise that are introduced into multivariate molecular data affect the estimation of the components of diversity (i.e., total number of different molecular species, or entities; total number of molecules; and the abundance distribution of molecular entities). We then explore under which conditions these biases affect the measurement of ɑ- and ß-diversity and highlight how novel methods commonly used in community ecology can be adopted to improve the interpretation and integration of multivariate molecular data. Video Abstract.

Minimum Information Standards in Chemistry: A Call for Better Research Data Management Practices.

Research data management (RDM) is needed to assist experimental advances and data collection in the chemical sciences. Many funders require RDM because experiments are often paid for by taxpayers and the resulting data should be deposited sustainably for posterity. However, paper notebooks are still common in laboratories and research data is often stored in proprietary and/or dead-end file formats without experimental context. Data must mature beyond a mere supplement to a research paper. Electronic lab notebooks (ELN) and laboratory information management systems (LIMS) allow researchers to manage data better and they simplify research and publication. Thus, an agreement is needed on minimum information standards for data handling to support structured approaches to data reporting. As digitalization becomes part of curricular teaching, future generations of digital native chemists will embrace RDM and ELN as an organic part of their research.

Networks and Graphs Discovery in Metabolomics Data Analysis and Interpretation.

Both targeted and untargeted mass spectrometry-based metabolomics approaches are used to understand the metabolic processes taking place in various organisms, from prokaryotes, plants, fungi to animals and humans. Untargeted approaches allow to detect as many metabolites as possible at once, identify unexpected metabolic changes, and characterize novel metabolites in biological samples. However, the identification of metabolites and the biological interpretation of such large and complex datasets remain challenging. One approach to address these challenges is considering that metabolites are connected through informative relationships. Such relationships can be formalized as networks, where the nodes correspond to the metabolites or features (when there is no or only partial identification), and edges connect nodes if the corresponding metabolites are related. Several networks can be built from a single dataset (or a list of metabolites), where each network represents different relationships, such as statistical (correlated metabolites), biochemical (known or putative substrates and products of reactions), or chemical (structural similarities, ontological relations). Once these networks are built, they can subsequently be mined using algorithms from network (or graph) theory to gain insights into metabolism. For instance, we can connect metabolites based on prior knowledge on enzymatic reactions, then provide suggestions for potential metabolite identifications, or detect clusters of co-regulated metabolites. In this review, we first aim at settling a nomenclature and formalism to avoid confusion when referring to different networks used in the field of metabolomics. Then, we present the state of the art of network-based methods for mass spectrometry-based metabolomics data analysis, as well as future developments expected in this area. We cover the use of networks applications using biochemical reactions, mass spectrometry features, chemical structural similarities, and correlations between metabolites. We also describe the application of knowledge networks such as metabolic reaction networks. Finally, we discuss the possibility of combining different networks to analyze and interpret them simultaneously.

A Modular and Expandable Ecosystem for Metabolomics Data Annotation in R.

Liquid chromatography-mass spectrometry (LC-MS)-based untargeted metabolomics experiments have become increasingly popular because of the wide range of metabolites that can be analyzed and the possibility to measure novel compounds. LC-MS instrumentation and analysis conditions can differ substantially among laboratories and experiments, thus resulting in non-standardized datasets demanding customized annotation workflows. We present an ecosystem of R packages, centered around the MetaboCoreUtils, MetaboAnnotation and CompoundDb packages that together provide a modular infrastructure for the annotation of untargeted metabolomics data. Initial annotation can be performed based on MS1 properties such as m/z and retention times, followed by an MS2-based annotation in which experimental fragment spectra are compared against a reference library. Such reference databases can be created and managed with the CompoundDb package. The ecosystem supports data from a variety of formats, including, but not limited to, MSP, MGF, mzML, mzXML, netCDF as well as MassBank text files and SQL databases. Through its highly customizable functionality, the presented infrastructure allows to build reproducible annotation workflows tailored for and adapted to most untargeted LC-MS-based datasets. All core functionality, which supports base R data types, is exported, also facilitating its re-use in other R packages. Finally, all packages are thoroughly unit-tested and documented and are available on GitHub and through Bioconductor.

Metabolic drift in the aging nervous system is reflected in human cerebrospinal fluid.

Chronic diseases affecting the central nervous system (CNS) like Alzheimer's or Parkinson's disease typically develop with advanced chronological age. Yet, aging at the metabolic level has been explored only sporadically in humans using biofluids in close proximity to the CNS such as the cerebrospinal fluid (CSF). We have used an untargeted liquid chromatography high-resolution mass spectrometry (LC-HRMS) based metabolomics approach to measure the levels of metabolites in the CSF of non-neurological control subjects in the age of 20 up to 74. Using a random forest-based feature selection strategy, we extracted 69 features that were strongly related to age (page < 0.001, rage = 0.762, R2Boruta age = 0.764). Combining an in-house library of known substances with in silico chemical classification and functional semantic annotation we successfully assigned putative annotations to 59 out of the 69 CSF metabolites. We found alterations in metabolites related to the Cytochrome P450 system, perturbations in the tryptophan and kynurenine pathways, metabolites associated with cellular energy (NAD+, ADP), mitochondrial and ribosomal metabolisms, neurological dysfunction, and an increase of adverse microbial metabolites. Taken together our results point at a key role for metabolites found in CSF related to the Cytochrome P450 system as most often associated with metabolic aging.

Mass spectrometry-based metabolomics: a guide for annotation, quantification and best reporting practices.

Mass spectrometry-based metabolomics approaches can enable detection and quantification of many thousands of metabolite features simultaneously. However, compound identification and reliable quantification are greatly complicated owing to the chemical complexity and dynamic range of the metabolome. Simultaneous quantification of many metabolites within complex mixtures can additionally be complicated by ion suppression, fragmentation and the presence of isomers. Here we present guidelines covering sample preparation, replication and randomization, quantification, recovery and recombination, ion suppression and peak misidentification, as a means to enable high-quality reporting of liquid chromatography- and gas chromatography-mass spectrometry-based metabolomics-derived data.

Modulation of Phosphate Deficiency-Induced Metabolic Changes by Iron Availability in Arabidopsis thaliana.

Concurrent suboptimal supply of several nutrients requires the coordination of nutrient-specific transcriptional, phenotypic, and metabolic changes in plants in order to optimize growth and development in most agricultural and natural ecosystems. Phosphate (Pi) and iron (Fe) deficiency induce overlapping but mostly opposing transcriptional and root growth responses in Arabidopsis thaliana. On the metabolite level, Pi deficiency negatively modulates Fe deficiency-induced coumarin accumulation, which is controlled by Fe as well as Pi deficiency response regulators. Here, we report the impact of Fe availability on seedling growth under Pi limiting conditions and on Pi deficiency-induced accumulation of amino acids and organic acids, which play important roles in Pi use efficiency. Fe deficiency in Pi replete conditions hardly changed growth and metabolite profiles in roots and shoots of Arabidopsis thaliana, but partially rescued growth under conditions of Pi starvation and severely modulated Pi deficiency-induced metabolic adjustments. Analysis of T-DNA insertion lines revealed the concerted coordination of metabolic profiles by regulators of Fe (FIT, bHLH104, BRUTUS, PYE) as well as of Pi (SPX1, PHR1, PHL1, bHLH32) starvation responses. The results show the interdependency of Pi and Fe availability and the interplay between Pi and Fe starvation signaling on the generation of plant metabolite profiles.

Untargeted In Silico Compound Classification-A Novel Metabolomics Method to Assess the Chemodiversity in Bryophytes.

In plant ecology, biochemical analyses of bryophytes and vascular plants are often conducted on dried herbarium specimen as species typically grow in distant and inaccessible locations. Here, we present an automated in silico compound classification framework to annotate metabolites using an untargeted data independent acquisition (DIA)-LC/MS-QToF-sequential windowed acquisition of all theoretical fragment ion mass spectra (SWATH) ecometabolomics analytical method. We perform a comparative investigation of the chemical diversity at the global level and the composition of metabolite families in ten different species of bryophytes using fresh samples collected on-site and dried specimen stored in a herbarium for half a year. Shannon and Pielou's diversity indices, hierarchical clustering analysis (HCA), sparse partial least squares discriminant analysis (sPLS-DA), distance-based redundancy analysis (dbRDA), ANOVA with post-hoc Tukey honestly significant difference (HSD) test, and the Fisher's exact test were used to determine differences in the richness and composition of metabolite families, with regard to herbarium conditions, ecological characteristics, and species. We functionally annotated metabolite families to biochemical processes related to the structural integrity of membranes and cell walls (proto-lignin, glycerophospholipids, carbohydrates), chemical defense (polyphenols, steroids), reactive oxygen species (ROS) protection (alkaloids, amino acids, flavonoids), nutrition (nitrogen- and phosphate-containing glycerophospholipids), and photosynthesis. Changes in the composition of metabolite families also explained variance related to ecological functioning like physiological adaptations of bryophytes to dry environments (proteins, peptides, flavonoids, terpenes), light availability (flavonoids, terpenes, carbohydrates), temperature (flavonoids), and biotic interactions (steroids, terpenes). The results from this study allow to construct chemical traits that can be attributed to biogeochemistry, habitat conditions, environmental changes and biotic interactions. Our classification framework accelerates the complex annotation process in metabolomics and can be used to simplify biochemical patterns. We show that compound classification is a powerful tool that allows to explore relationships in both molecular biology by "zooming in" and in ecology by "zooming out". The insights revealed by our framework allow to construct new research hypotheses and to enable detailed follow-up studies.

Empowering large chemical knowledge bases for exposomics: PubChemLite meets MetFrag.

Compound (or chemical) databases are an invaluable resource for many scientific disciplines. Exposomics researchers need to find and identify relevant chemicals that cover the entirety of potential (chemical and other) exposures over entire lifetimes. This daunting task, with over 100 million chemicals in the largest chemical databases, coupled with broadly acknowledged knowledge gaps in these resources, leaves researchers faced with too much-yet not enough-information at the same time to perform comprehensive exposomics research. Furthermore, the improvements in analytical technologies and computational mass spectrometry workflows coupled with the rapid growth in databases and increasing demand for high throughput "big data" services from the research community present significant challenges for both data hosts and workflow developers. This article explores how to reduce candidate search spaces in non-target small molecule identification workflows, while increasing content usability in the context of environmental and exposomics analyses, so as to profit from the increasing size and information content of large compound databases, while increasing efficiency at the same time. In this article, these methods are explored using PubChem, the NORMAN Network Suspect List Exchange and the in silico fragmentation approach MetFrag. A subset of the PubChem database relevant for exposomics, PubChemLite, is presented as a database resource that can be (and has been) integrated into current workflows for high resolution mass spectrometry. Benchmarking datasets from earlier publications are used to show how experimental knowledge and existing datasets can be used to detect and fill gaps in compound databases to progressively improve large resources such as PubChem, and topic-specific subsets such as PubChemLite. PubChemLite is a living collection, updating as annotation content in PubChem is updated, and exported to allow direct integration into existing workflows such as MetFrag. The source code and files necessary to recreate or adjust this are jointly hosted between the research parties (see data availability statement). This effort shows that enhancing the FAIRness (Findability, Accessibility, Interoperability and Reusability) of open resources can mutually enhance several resources for whole community benefit. The authors explicitly welcome additional community input on ideas for future developments.

LC-MS based plant metabolic profiles of thirteen grassland species grown in diverse neighbourhoods.

In plants, secondary metabolite profiles provide a unique opportunity to explore seasonal variation and responses to the environment. These include both abiotic and biotic factors. In field experiments, such stress factors occur in combination. This variation alters the plant metabolic profiles in yet uninvestigated ways. This data set contains trait and mass spectrometry data of thirteen grassland species collected at four time points in the growing season in 2017. We collected above-ground vegetative material of seven grass and six herb species that were grown in plant communities with different levels of diversity in the Jena Experiment. For each sample, we recorded visible traits and acquired shoot metabolic profiles on a UPLC-ESI-Qq-TOF-MS. We performed the raw data pre-processing in Galaxy-W4M and prepared the data for statistical analysis in R by applying missing data imputation, batch correction, and validity checks on the features. This comprehensive data set provides the opportunity to investigate environmental dynamics across diverse neighbourhoods that are reflected in the metabolomic profile.

ELIXIR and Toxicology: a community in development.

Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities.

Reshaping of the Arabidopsis thaliana Proteome Landscape and Co-regulation of Proteins in Development and Immunity.

Proteome remodeling is a fundamental adaptive response, and proteins in complexes and functionally related proteins are often co-expressed. Using a deep sampling strategy we define core proteomes of Arabidopsis thaliana tissues with around 10 000 proteins per tissue, and absolutely quantify (copy numbers per cell) nearly 16 000 proteins throughout the plant lifecycle. A proteome-wide survey of global post-translational modification revealed amino acid exchanges pointing to potential conservation of translational infidelity in eukaryotes. Correlation analysis of protein abundance uncovered potentially new tissue- and age-specific roles of entire signaling modules regulating transcription in photosynthesis, seed development, and senescence and abscission. Among others, the data suggest a potential function of RD26 and other NAC transcription factors in seed development related to desiccation tolerance as well as a possible function of cysteine-rich receptor-like kinases (CRKs) as ROS sensors in senescence. All of the components of ribosome biogenesis factor (RBF) complexes were found to be co-expressed in a tissue- and age-specific manner, indicating functional promiscuity in the assembly of these less-studied protein complexes in Arabidopsis.Furthermore, we characterized detailed proteome remodeling in basal immunity by treating Arabidopsis seeldings with flg22. Through simultaneously monitoring phytohormone and transcript changes upon flg22 treatment, we obtained strong evidence of suppression of jasmonate (JA) and JA-isoleucine (JA-Ile) levels by deconjugation and hydroxylation by IAA-ALA RESISTANT3 (IAR3) and JASMONATE-INDUCED OXYGENASE 2 (JOX2), respectively, under the control of JASMONATE INSENSITIVE 1 (MYC2), suggesting an unrecognized role of a new JA regulatory switch in pattern-triggered immunity. Taken together, the datasets generated in this study present extensive coverage of the Arabidopsis proteome in various biological scenarios, providing a rich resource available to the whole plant science community.

Feature-based molecular networking in the GNPS analysis environment.

Molecular networking has become a key method to visualize and annotate the chemical space in non-targeted mass spectrometry data. We present feature-based molecular networking (FBMN) as an analysis method in the Global Natural Products Social Molecular Networking (GNPS) infrastructure that builds on chromatographic feature detection and alignment tools. FBMN enables quantitative analysis and resolution of isomers, including from ion mobility spectrometry.

Chemical Diversity and Classification of Secondary Metabolites in Nine Bryophyte Species.

The central aim in ecometabolomics and chemical ecology is to pinpoint chemical features that explain molecular functioning. The greatest challenge is the identification of compounds due to the lack of constitutive reference spectra, the large number of completely unknown compounds, and bioinformatic methods to analyze the big data. In this study we present an interdisciplinary methodological framework that extends ultra-performance liquid chromatography coupled to electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC/ESI-QTOF-MS) with data-dependent acquisition (DDA-MS) and the automated in silico classification of fragment peaks into compound classes. We synthesize findings from a prior study that explored the influence of seasonal variations on the chemodiversity of secondary metabolites in nine bryophyte species. Here we reuse and extend the representative dataset with DDA-MS data. Hierarchical clustering, heatmaps, dbRDA, and ANOVA with post-hoc Tukey HSD were used to determine relationships of the study factors species, seasons, and ecological characteristics. The tested bryophytes showed species-specific metabolic responses to seasonal variations (50% vs. 5% of explained variation). Marchantia polymorpha, Plagiomnium undulatum, and Polytrichum strictum were biochemically most diverse and unique. Flavonoids and sesquiterpenoids were upregulated in all bryophytes in the growing seasons. We identified ecological functioning of compound classes indicating light protection (flavonoids), biotic and pathogen interactions (sesquiterpenoids, flavonoids), low temperature and desiccation tolerance (glycosides, sesquiterpenoids, anthocyanins, lactones), and moss growth supporting anatomic structures (few methoxyphenols and cinnamic acids as part of proto-lignin constituents). The reusable bioinformatic framework of this study can differentiate species based on automated compound classification. Our study allows detailed insights into the ecological roles of biochemical constituents of bryophytes with regard to seasonal variations. We demonstrate that compound classification can be improved with adding constitutive reference spectra to existing spectral libraries. We also show that generalization on compound classes improves our understanding of molecular ecological functioning and can be used to generate new research hypotheses.