The potential of untargeted metabolomics to answer important questions across the life sciences is hindered because of a paucity of computational tools that enable extraction of key biochemically relevant information. Available tools focus on using mass spectrometry fragmentation spectra to identify molecules whose behavior suggests they are relevant to the system under study. Unfortunately, fragmentation spectra cannot identify molecules in isolation but require authentic standards or databases of known fragmented molecules. Fragmentation spectra are, however, replete with information pertaining to the biochemical processes present, much of which is currently neglected. Here, we present an analytical workflow that exploits all fragmentation data from a given experiment to extract biochemically relevant features in an unsupervised manner. We demonstrate that an algorithm originally used for text mining, latent Dirichlet allocation, can be adapted to handle metabolomics datasets. Our approach extracts biochemically relevant molecular substructures ("Mass2Motifs") from spectra as sets of co-occurring molecular fragments and neutral losses. The analysis allows us to isolate molecular substructures, whose presence allows molecules to be grouped based on shared substructures regardless of classical spectral similarity. These substructures, in turn, support putative de novo structural annotation of molecules. Combining this spectral connectivity to orthogonal correlations (e.g., common abundance changes under system perturbation) significantly enhances our ability to provide mechanistic explanations for biological behavior.metabolomics | mass spectrometry | fragmentation | bioinformatics | topic modeling M ass spectrometry (MS)-based metabolomics aims to capture the entire small-molecule composition of biological systems. Analysis of MS metabolomics data are challenging as many molecules cannot be identified from their mass (e.g., isobaric molecules, and isomers) (1-3). Separation by liquid chromatography before MS (LC-MS) can add discriminatory information but does not solve the problem as isomers can exhibit similar chromatographic behavior, and chromatographic retention time is currently unpredictable.Fragmentation spectra have been used to partially overcome this problem (4-6). Most tools compare individual fragmentation spectra to reference spectra (5, 7) stored in public databases, for example, MassBank (8) or Human Metabolome Database (9), and are thus constrained by the limited number of reference spectra (10-12). Poor identification coverage can result in poor biochemical insight. We propose a method that analyzes all acquired fragmentation spectra to expose underlying biochemistry without relying on metabolite identification, inspired by machine-learning techniques developed initially for text processing (13).The paucity of techniques that share information across fragmentation spectra can be explained by the complexity of fragmentation data (14). One example, "Molecular Networking," clusters MS1 peaks by th...