TopNEXt: automatic DDA exclusion framework for multi-sample mass spectrometry experiments.

MOTIVATION: Liquid Chromatography Tandem Mass Spectrometry experiments aim to produce high-quality fragmentation spectra, which can be used to annotate metabolites. However, current Data-Dependent Acquisition approaches may fail to collect spectra of sufficient quality and quantity for experimental outcomes, and extend poorly across multiple samples by failing to share information across samples or by requiring manual expert input. RESULTS: We present TopNEXt, a real-time scan prioritization framework that improves data acquisition in multi-sample Liquid Chromatography Tandem Mass Spectrometry metabolomics experiments. TopNEXt extends traditional Data-Dependent Acquisition exclusion methods across multiple samples by using a Region of Interest and intensity-based scoring system. Through both simulated and lab experiments, we show that methods incorporating these novel concepts acquire fragmentation spectra for an additional 10% of our set of target peaks and with an additional 20% of acquisition intensity. By increasing the quality and quantity of fragmentation spectra, TopNEXt can help improve metabolite identification with a potential impact across a variety of experimental contexts. AVAILABILITY AND IMPLEMENTATION: TopNEXt is implemented as part of the ViMMS framework and the latest version can be found at https://github.com/glasgowcompbio/vimms. A stable version used to produce our results can be found at 10.5281/zenodo.7468914.

On-line targeted metabolomics for real-time monitoring of relevant compounds in fermentation processes.

Fermentation monitoring is a powerful tool for bioprocess development and optimization. On-line metabolomics is a technology that is starting to gain attention as a bioprocess monitoring tool, allowing the direct measurement of many compounds in the fermentation broth at a very high time resolution. In this work, targeted on-line metabolomics was used to monitor 40 metabolites of interest during three Escherichia coli succinate production fermentation experiments every 5 min with a triple quadrupole mass spectrometer. This allowed capturing high-time resolution biological data that can provide critical information for process optimization. For nine of these metabolites, simple univariate regression models were used to model compound concentration from their on-line mass spectrometry peak area. These on-line metabolomics univariate models performed comparably to vibrational spectroscopy multivariate partial least squares regressions models reported in the literature, which typically are much more complex and time consuming to build. In conclusion, this work shows how on-line metabolomics can be used to directly monitor many bioprocess compounds of interest and obtain rich biological and bioprocess data.

On-line untargeted metabolomics monitoring of an Escherichia coli succinate fermentation process.

The real-time monitoring of metabolites (RTMet) is instrumental for the industrial production of biobased fermentation products. This study shows the first application of untargeted on-line metabolomics for the monitoring of undiluted fermentation broth samples taken automatically from a 5 L bioreactor every 5 min via flow injection mass spectrometry. The travel time from the bioreactor to the mass spectrometer was 30 s. Using mass spectrometry allows, on the one hand, the direct monitoring of targeted key process compounds of interest and, on the other hand, provides information on hundreds of additional untargeted compounds without requiring previous calibration data. In this study, this technology was applied in an Escherichia coli succinate fermentation process and 886 different m/z signals were monitored, including key process compounds (glucose, succinate, and pyruvate), potential biomarkers of biomass formation such as (R)-2,3-dihydroxy-isovalerate and (R)-2,3-dihydroxy-3-methylpentanoate and compounds from the pentose phosphate pathway and nucleotide metabolism, among others. The main advantage of the RTMet technology is that it allows the monitoring of hundreds of signals without the requirement of developing partial least squares regression models, making it a perfect tool for bioprocess monitoring and for testing many different strains and process conditions for bioprocess development.

Ranking microbial metabolomic and genomic links in the NPLinker framework using complementary scoring functions.

Specialised metabolites from microbial sources are well-known for their wide range of biomedical applications, particularly as antibiotics. When mining paired genomic and metabolomic data sets for novel specialised metabolites, establishing links between Biosynthetic Gene Clusters (BGCs) and metabolites represents a promising way of finding such novel chemistry. However, due to the lack of detailed biosynthetic knowledge for the majority of predicted BGCs, and the large number of possible combinations, this is not a simple task. This problem is becoming ever more pressing with the increased availability of paired omics data sets. Current tools are not effective at identifying valid links automatically, and manual verification is a considerable bottleneck in natural product research. We demonstrate that using multiple link-scoring functions together makes it easier to prioritise true links relative to others. Based on standardising a commonly used score, we introduce a new, more effective score, and introduce a novel score using an Input-Output Kernel Regression approach. Finally, we present NPLinker, a software framework to link genomic and metabolomic data. Results are verified using publicly available data sets that include validated links.

GraphOmics: an interactive platform to explore and integrate multi-omics data.

BACKGROUND: An increasing number of studies now produce multiple omics measurements that require using sophisticated computational methods for analysis. While each omics data can be examined separately, jointly integrating multiple omics data allows for deeper understanding and insights to be gained from the study. In particular, data integration can be performed horizontally, where biological entities from multiple omics measurements are mapped to common reactions and pathways. However, data integration remains a challenge due to the complexity of the data and the difficulty in interpreting analysis results. RESULTS: Here we present GraphOmics, a user-friendly platform to explore and integrate multiple omics datasets and support hypothesis generation. Users can upload transcriptomics, proteomics and metabolomics data to GraphOmics. Relevant entities are connected based on their biochemical relationships, and mapped to reactions and pathways from Reactome. From the Data Browser in GraphOmics, mapped entities and pathways can be ranked, sorted and filtered according to their statistical significance (p values) and fold changes. Context-sensitive panels provide information on the currently selected entities, while interactive heatmaps and clustering functionalities are also available. As a case study, we demonstrated how GraphOmics was used to interactively explore multi-omics data and support hypothesis generation using two complex datasets from existing Zebrafish regeneration and Covid-19 human studies. CONCLUSIONS: GraphOmics is fully open-sourced and freely accessible from https://graphomics.glasgowcompbio.org/ . It can be used to integrate multiple omics data horizontally by mapping entities across omics to reactions and pathways. Our demonstration showed that by using interactive explorations from GraphOmics, interesting insights and biological hypotheses could be rapidly revealed.

Rapid Development of Improved Data-Dependent Acquisition Strategies.

Tandem mass spectrometry (LC-MS/MS) is widely used to identify unknown ions in untargeted metabolomics. Data-dependent acquisition (DDA) chooses which ions to fragment based upon intensities observed in MS1 survey scans and typically only fragments a small subset of the ions present. Despite this inefficiency, relatively little work has addressed the development of new DDA methods, partly due to the high overhead associated with running the many extracts necessary to optimize approaches in busy MS facilities. In this work, we first provide theoretical results that show how much improvement is possible over current DDA strategies. We then describe an in silico framework for fast and cost-efficient development of new DDA strategies using a previously developed virtual metabolomics mass spectrometer (ViMMS). Additional functionality is added to ViMMS to allow methods to be used both in simulation and on real samples via an Instrument Application Programming Interface (IAPI). We demonstrate this framework through the development and optimization of two new DDA methods that introduce new advanced ion prioritization strategies. Upon application of these developed methods to two complex metabolite mixtures, our results show that they are able to fragment more unique ions than standard DDA strategies.

Changing environments and genetic variation: natural variation in inbreeding does not compromise short-term physiological responses.

Selfing plant lineages are surprisingly widespread and successful in a broad range of environments, despite showing reduced genetic diversity, which is predicted to reduce their long-term evolutionary potential. However, appropriate short-term plastic responses to new environmental conditions might not require high levels of standing genetic variation. In this study, we tested whether mating system variation among populations, and associated changes in genetic variability, affected short-term responses to environmental challenges. We compared relative fitness and metabolome profiles of naturally outbreeding (genetically diverse) and inbreeding (genetically depauperate) populations of a perennial plant, Arabidopsis lyrata, under constant growth chamber conditions and an outdoor common garden environment outside its native range. We found no effect of inbreeding on survival, flowering phenology or short-term physiological responses. Specifically, naturally occurring inbreeding had no significant effects on the plasticity of metabolome profiles, using either multivariate approaches or analysis of variation in individual metabolites, with inbreeding populations showing similar physiological responses to outbreeding populations over time in both growing environments. We conclude that low genetic diversity in naturally inbred populations may not always compromise fitness or short-term physiological capacity to respond to environmental change, which could help to explain the global success of selfing mating strategies.

MetaboCraft: building a Minecraft plugin for metabolomics.

Motivation: The rapid advances in metabolomics pose a significant challenge in presentation and interpretation of results. Development of new, engaging visual aids is crucial to advancing our understanding of new findings. Results: We have developed MetaboCraft, a Minecraft plugin which creates immersive visualizations of metabolic networks and pathways in a 3D environment and allows the results of user experiments to be viewed in this context, presenting a novel approach to exploring the metabolome. Availability and implementation: https://github.com/argymeg/MetaboCraft/; https://hub.docker.com/r/ronandaly/metabocraft/. Supplementary information: Supplementary data are available at Bioinformatics online.

ShinyKGode: an interactive application for ODE parameter inference using gradient matching.

Motivation: Mathematical modelling based on ordinary differential equations (ODEs) is widely used to describe the dynamics of biological systems, particularly in systems and pathway biology. Often the kinetic parameters of these ODE systems are unknown and have to be inferred from the data. Approximate parameter inference methods based on gradient matching (which do not require performing computationally expensive numerical integration of the ODEs) have been getting popular in recent years, but many implementations are difficult to run without expert knowledge. Here, we introduce ShinyKGode, an interactive web application to perform fast parameter inference on ODEs using gradient matching. Results: ShinyKGode can be used to infer ODE parameters on simulated and observed data using gradient matching. Users can easily load their own models in Systems Biology Markup Language format, and a set of pre-defined ODE benchmark models are provided in the application. Inferred parameters are visualized alongside diagnostic plots to assess convergence. Availability and implementation: The R package for ShinyKGode can be installed through the Comprehensive R Archive Network (CRAN). Installation instructions, as well as tutorial videos and source code are available at https://joewandy.github.io/shinyKGode. Supplementary information: Supplementary data are available at Bioinformatics online.

Ms2lda.org: web-based topic modelling for substructure discovery in mass spectrometry.

MOTIVATION: We recently published MS2LDA, a method for the decomposition of sets of molecular fragment data derived from large metabolomics experiments. To make the method more widely available to the community, here we present ms2lda.org, a web application that allows users to upload their data, run MS2LDA analyses and explore the results through interactive visualizations. RESULTS: Ms2lda.org takes tandem mass spectrometry data in many standard formats and allows the user to infer the sets of fragment and neutral loss features that co-occur together (Mass2Motifs). As an alternative workflow, the user can also decompose a data set onto predefined Mass2Motifs. This is accomplished through the web interface or programmatically from our web service. AVAILABILITY AND IMPLEMENTATION: The website can be found at http://ms2lda.org, while the source code is available at https://github.com/sdrogers/ms2ldaviz under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PiMP my metabolome: an integrated, web-based tool for LC-MS metabolomics data.

SUMMARY: The Polyomics integrated Metabolomics Pipeline (PiMP) fulfils an unmet need in metabolomics data analysis. PiMP offers automated and user-friendly analysis from mass spectrometry data acquisition to biological interpretation. Our key innovations are the Summary Page, which provides a simple overview of the experiment in the format of a scientific paper, containing the key findings of the experiment along with associated metadata; and the Metabolite Page, which provides a list of each metabolite accompanied by 'evidence cards', which provide a variety of criteria behind metabolite annotation including peak shapes, intensities in different sample groups and database information. AVAILABILITY AND IMPLEMENTATION: PiMP is available at http://polyomics.mvls.gla.ac.uk, and access is freely available on request. 50 GB of space is allocated for data storage, with unrestricted number of samples and analyses per user. Source code is available at https://github.com/RonanDaly/pimp and licensed under the GPL. CONTACT: karl.burgess@glasgow.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enhanced Quality Factor Label-free Biosensing with Micro-Cantilevers Integrated into Microfluidic Systems.

Microelectromechanical systems (MEMS) have enabled the development of a new generation of sensor platforms. Acoustic sensor operation in liquid, the native environment of biomolecules, causes, however, significant degradation of sensing performance due to viscous drag and relies on the availability of capture molecules to bind analytes of interest to the sensor surface. Here, we describe a strategy to interface MEMS sensors with microfluidic platforms through an aerosol spray. Our sensing platform comprises a microfluidic spray nozzle and a microcantilever array operated in dynamic mode within a closed loop oscillator. A solution containing the analyte is sprayed uniformly through picoliter droplets onto the microcantilever surface; the micrometer-scale drops evaporate rapidly and leave the solutes behind, adding to the mass of the cantilever. This sensing scheme results in a 50-fold increase in the quality factor compared to operation in liquid, yet allows the analytes to be introduced into the sensing system from a solution phase. It achieves a 370 femtogram limit of detection, and we demonstrate quantitative label-free analysis of inorganic salts and model proteins. These results demonstrate that the standard resolution limits of cantilever sensing in dynamic mode can be overcome with the integration of spray microfluidics with MEMS.

Incorporating peak grouping information for alignment of multiple liquid chromatography-mass spectrometry datasets.

MOTIVATION: The combination of liquid chromatography and mass spectrometry (LC/MS) has been widely used for large-scale comparative studies in systems biology, including proteomics, glycomics and metabolomics. In almost all experimental design, it is necessary to compare chromatograms across biological or technical replicates and across sample groups. Central to this is the peak alignment step, which is one of the most important but challenging preprocessing steps. Existing alignment tools do not take into account the structural dependencies between related peaks that coelute and are derived from the same metabolite or peptide. We propose a direct matching peak alignment method for LC/MS data that incorporates related peaks information (within each LC/MS run) and investigate its effect on alignment performance (across runs). The groupings of related peaks necessary for our method can be obtained from any peak clustering method and are built into a pair-wise peak similarity score function. The similarity score matrix produced is used by an approximation algorithm for the weighted matching problem to produce the actual alignment result. RESULTS: We demonstrate that related peak information can improve alignment performance. The performance is evaluated on a set of benchmark datasets, where our method performs competitively compared to other popular alignment tools. AVAILABILITY: The proposed alignment method has been implemented as a stand-alone application in Python, available for download at http://github.com/joewandy/peak-grouping-alignment.

MetAssign: probabilistic annotation of metabolites from LC-MS data using a Bayesian clustering approach.

MOTIVATION: The use of liquid chromatography coupled to mass spectrometry has enabled the high-throughput profiling of the metabolite composition of biological samples. However, the large amount of data obtained can be difficult to analyse and often requires computational processing to understand which metabolites are present in a sample. This article looks at the dual problem of annotating peaks in a sample with a metabolite, together with putatively annotating whether a metabolite is present in the sample. The starting point of the approach is a Bayesian clustering of peaks into groups, each corresponding to putative adducts and isotopes of a single metabolite. RESULTS: The Bayesian modelling introduced here combines information from the mass-to-charge ratio, retention time and intensity of each peak, together with a model of the inter-peak dependency structure, to increase the accuracy of peak annotation. The results inherently contain a quantitative estimate of confidence in the peak annotations and allow an accurate trade-off between precision and recall. Extensive validation experiments using authentic chemical standards show that this system is able to produce more accurate putative identifications than other state-of-the-art systems, while at the same time giving a probabilistic measure of confidence in the annotations. AVAILABILITY AND IMPLEMENTATION: The software has been implemented as part of the mzMatch metabolomics analysis pipeline, which is available for download at http://mzmatch.sourceforge.net/.

Cross-Linking the Fibers of Supramolecular Gels Formed from a Tripodal Terpyridine Derived Ligand with d-Block Metal Ions.

The tripodal terpyridine ligand, L, forms 1D helical supramolecular polymers/gels in H2O-CH3OH solution mediated through hydrogen bonding and π-π interactions. These gels further cross-link into 3D supramolecular metallogels with a range of metal ions (M) such as Fe(II), Ni(II), Cu(II), Zn(II), and Ru(III); the cross-linking resulting in the formation of colored or colorless gels. The fibrous morphology of these gels was confirmed using scanning electron microscopy (SEM); while the self-assembly processes between L and M were investigated by absorbance and emission spectroscopy from which their binding constants were determined by using a nonlinear regression analysis.

Quantitative monitoring and modelling of retrodialysis drug delivery in a brain phantom.

A vast number of drug molecules are unable to cross the blood-brain barrier, which results in a loss of therapeutic opportunities when these molecules are administered by intravenous infusion. To circumvent the blood-brain barrier, local drug delivery devices have been developed over the past few decades such as reverse microdialysis. Reverse microdialysis (or retrodialysis) offers many advantages, such as a lack of net volume influx to the intracranial cavity and the ability to sample the tumour's micro-environment. However, the translation of this technique to efficient drug delivery has not been systematically studied. In this work, we present an experimental platform to evaluate the performance of microdialysis devices in reverse mode in a brain tissue phantom. The mass of model drug delivered is measured by computing absorbance fields from optical images. Concentration maps are reconstructed using a modern and open-source implementation of the inverse Abel transform. To illustrate our method, we assess the capability of a commercial probe in delivering methylene blue to a gel phantom. We find that the delivery rate can be described by classical microdialysis theory, except at low dialysate flow rates where it is impacted by gravity, and high flow rates where significant convection to the gel occurs. We also show that the flow rate has an important impact not only on the overall size of the drug plume, but also on its shape. The numerical tools developed for this study have been made freely available to ensure that the method presented can be used to rapidly and inexpensively optimise probe design and protocol parameters before proceeding to more in-depth studies.

Streamlined identification of strain engineering targets for bioprocess improvement using metabolic pathway enrichment analysis.

Metabolomics is a powerful tool for the identification of genetic targets for bioprocess optimisation. However, in most cases, only the biosynthetic pathway directed to product formation is analysed, limiting the identification of these targets. Some studies have used untargeted metabolomics, allowing a more unbiased approach, but data interpretation using multivariate analysis is usually not straightforward and requires time and effort. Here we show, for the first time, the application of metabolic pathway enrichment analysis using untargeted and targeted metabolomics data to identify genetic targets for bioprocess improvement in a more streamlined way. The analysis of an Escherichia coli succinate production bioprocess with this methodology revealed three significantly modulated pathways during the product formation phase: the pentose phosphate pathway, pantothenate and CoA biosynthesis and ascorbate and aldarate metabolism. From these, the two former pathways are consistent with previous efforts to improve succinate production in Escherichia coli. Furthermore, to the best of our knowledge, ascorbate and aldarate metabolism is a newly identified target that has so far never been explored for improving succinate production in this microorganism. This methodology therefore represents a powerful tool for the streamlined identification of strain engineering targets that can accelerate bioprocess optimisation.

Simulated-to-real benchmarking of acquisition methods in untargeted metabolomics.

Data-Dependent and Data-Independent Acquisition modes (DDA and DIA, respectively) are both widely used to acquire MS2 spectra in untargeted liquid chromatography tandem mass spectrometry (LC-MS/MS) metabolomics analyses. Despite their wide use, little work has been attempted to systematically compare their MS/MS spectral annotation performance in untargeted settings due to the lack of ground truth and the costs involved in running a large number of acquisitions. Here, we present a systematic in silico comparison of these two acquisition methods in untargeted metabolomics by extending our Virtual Metabolomics Mass Spectrometer (ViMMS) framework with a DIA module. Our results show that the performance of these methods varies with the average number of co-eluting ions as the most important factor. At low numbers, DIA outperforms DDA, but at higher numbers, DDA has an advantage as DIA can no longer deal with the large amount of overlapping ion chromatograms. Results from simulation were further validated on an actual mass spectrometer, demonstrating that using ViMMS we can draw conclusions from simulation that translate well into the real world. The versatility of the Virtual Metabolomics Mass Spectrometer (ViMMS) framework in simulating different parameters of both Data-Dependent and Data-Independent Acquisition (DDA and DIA) modes is a key advantage of this work. Researchers can easily explore and compare the performance of different acquisition methods within the ViMMS framework, without the need for expensive and time-consuming experiments with real experimental data. By identifying the strengths and limitations of each acquisition method, researchers can optimize their choice and obtain more accurate and robust results. Furthermore, the ability to simulate and validate results using the ViMMS framework can save significant time and resources, as it eliminates the need for numerous experiments. This work not only provides valuable insights into the performance of DDA and DIA, but it also opens the door for further advancements in LC-MS/MS data acquisition methods.

Outcomes of Uterine Rupture in the Setting of the Unscarred Compared With the Scarred Uterus.

Uterine rupture is a rare obstetric complication that is associated with maternal and neonatal morbidity and mortality. The aim of this study was to examine uterine rupture and its outcomes in the setting of the unscarred compared with the scarred uterus. A retrospective observational cohort study was performed examining all cases of uterine rupture in three tertiary care hospitals in Dublin, Ireland, over a 20-year period. The primary outcome was perinatal mortality rate with uterine rupture, which was 11.02% (95% CI 6.5-17.3). There was no significant difference in perinatal mortality between cases of scarred and unscarred uterine rupture. Unscarred uterine rupture was associated with higher maternal morbidity , defined as major obstetric hemorrhage or hysterectomy.