Secondary metabolite profiling of Pseudomonas aeruginosa isolates reveals rare genomic traits.

Pseudomonas aeruginosa is a ubiquitous Gram-negative opportunistic pathogen with remarkable phylogenetic and phenotypic variabilities. In this work, we applied classical molecular networking analysis to secondary metabolite profiling data from seven Pseudomonas aeruginosa strains, including five clinical isolates from the lung secretions of people with cystic fibrosis (CF). We provide three vignettes illustrating how secondary metabolite profiling aids in the identification of rare genomics traits in P. aeruginosa. First, we describe the identification of a previously unreported class of acyl putrescines produced by isolate mFLRO1. Secondary analysis of publicly available metabolomics data revealed that acyl putrescines are produced by <5% of P. aeruginosa strains. Second, we show that isolate SH3A does not produce di-rhamnolipids. Whole-genome sequencing and comparative genomics revealed that SH3A cannot produce di-rhamnolipids because its genome belongs to clade 5 of the P. aeruginosa phylogenetic tree. Previous phylogenetic analysis of thousands of P. aeruginosa strains concluded that <1% of publicly available genome sequences contribute to this clade. Last, we show that isolate SH1B does not produce the phenazine pyocyanin or rhamnolipids because it has a one-base insertion frameshift mutation (678insC) in the gene rhlR, which disrupts rhl-driven quorum sensing. Secondary analysis of the tens of thousands of publicly available genomes in the National Center for Biotechnology Information (NCBI) and the Pseudomonas Genome Database revealed that this mutation was present in only four P. aeruginosa genomes. Taken together, this study highlights that secondary metabolite profiling combined with genomic analysis can identify rare genetic traits of P. aeruginosa isolates.IMPORTANCESecondary metabolite profiling of five Pseudomonas aeruginosa isolates from cystic fibrosis sputum captured three traits present in <1%-5% of publicly available data, pointing to how our current library of P. aeruginosa strains may not represent the diversity within this species or the genetic variance that occurs in the CF lung.

ModiFinder: Tandem Mass Spectral Alignment Enables Structural Modification Site Localization.

Untargeted tandem mass spectrometry (MS/MS) has become a high-throughput method to measure small molecules in complex samples. One key goal is the transformation of these MS/MS spectra into chemical structures. Computational techniques such as MS/MS library search have enabled the reidentification of known compounds. Analog library search and molecular networking extend this identification to unknown compounds. While there have been advancements in metrics for the similarity of MS/MS spectra of structurally similar compounds, there is still a lack of automated methods to provide site specific information about structural modifications. Here we introduce ModiFinder which leverages the alignment of peaks in MS/MS spectra between structurally related known and unknown small molecules. Specifically, ModiFinder focuses on shifted MS/MS fragment peaks in the MS/MS alignment. These shifted peaks putatively represent substructures of the known molecule that contain the site of the modification. ModiFinder synthesizes this information together and scores the likelihood for each atom in the known molecule to be the modification site. We demonstrate in this manuscript how ModiFinder can effectively localize modifications which extends the capabilities of MS/MS analog searching and molecular networking to accelerate the discovery of novel compounds.

Network Topology Evaluation and Transitive Alignments for Molecular Networking.

Untargeted tandem mass spectrometry (MS/MS) is an essential technique in modern analytical chemistry, providing a comprehensive snapshot of chemical entities in complex samples and identifying unknowns through their fragmentation patterns. This high-throughput approach generates large data sets that can be challenging to interpret. Molecular Networks (MNs) have been developed as a computational tool to aid in the organization and visualization of complex chemical space in untargeted mass spectrometry data, thereby supporting comprehensive data analysis and interpretation. MNs group related compounds with potentially similar structures from MS/MS data by calculating all pairwise MS/MS similarities and filtering these connections to produce a MN. Such networks are instrumental in metabolomics for identifying novel metabolites, elucidating metabolic pathways, and even discovering biomarkers for disease. While MS/MS similarity metrics have been explored in the literature, the influence of network topology approaches on MN construction remains unexplored. This manuscript introduces metrics for evaluating MN construction, benchmarks state-of-the-art approaches, and proposes the Transitive Alignments approach to improve MN construction. The Transitive Alignment technique leverages the MN topology to realign MS/MS spectra of related compounds that differ by multiple structural modifications. Combining this Transitive Alignments approach with pseudoclique finding, a method for identifying highly connected groups of nodes in a network, resulted in more complete and higher-quality molecular families. Finally, we also introduce a targeted network construction technique called induced transitive alignments where we demonstrate effectiveness on a real world natural product discovery application. We release this transitive alignment technique as a high-throughput workflow that can be used by the wider research community.