Reverse metabolomics for the discovery of chemical structures from humans.

Determining the structure and phenotypic context of molecules detected in untargeted metabolomics experiments remains challenging. Here we present reverse metabolomics as a discovery strategy, whereby tandem mass spectrometry spectra acquired from newly synthesized compounds are searched for in public metabolomics datasets to uncover phenotypic associations. To demonstrate the concept, we broadly synthesized and explored multiple classes of metabolites in humans, including N-acyl amides, fatty acid esters of hydroxy fatty acids, bile acid esters and conjugated bile acids. Using repository-scale analysis1,2, we discovered that some conjugated bile acids are associated with inflammatory bowel disease (IBD). Validation using four distinct human IBD cohorts showed that cholic acids conjugated to Glu, Ile/Leu, Phe, Thr, Trp or Tyr are increased in Crohn's disease. Several of these compounds and related structures affected pathways associated with IBD, such as interferon-γ production in CD4+ T cells3 and agonism of the pregnane X receptor4. Culture of bacteria belonging to the Bifidobacterium, Clostridium and Enterococcus genera produced these bile amidates. Because searching repositories with tandem mass spectrometry spectra has only recently become possible, this reverse metabolomics approach can now be used as a general strategy to discover other molecules from human and animal ecosystems.

Elucidating the role of lipid interactions in stabilizing the membrane protein KcsA.

The significant effects of lipid binding on the functionality of potassium channel KcsA have been validated by brilliant studies. However, the specific interactions between lipids and KcsA, such as binding parameters for each binding event, have not been fully elucidated. In this study, we employed native mass spectrometry to investigate the binding of lipids to KcsA and their effects on the channel. The tetrameric structure of KcsA remains intact even in the absence of lipid binding. However, the subunit architecture of the E71A mutant, which is constantly open at low pH, relies on tightly associated copurified lipids. Furthermore, we observed that lipids exhibit weak binding to KcsA at high pH when the channel is at a closed/inactivation state in the absence of permeant cation K+. This feeble interaction potentially facilitates the association of K+ ions, leading to the transition of the channel to a resting closed/open state. Interestingly, both anionic and zwitterionic lipids strongly bind to KcsA at low pH when the channel is in an open/inactivation state. We also investigated the binding patterns of KcsA with natural lipids derived from E. coli and Streptomyces lividans. Interestingly, lipids from E. coli exhibited much stronger binding affinity compared to the lipids from S. lividans. Among the natural lipids from S. lividans, free fatty acids and triacylglycerols demonstrated the tightest binding to KcsA, whereas no detectable binding events were observed with natural phosphatidic acid lipids. These findings suggest that the lipid association pattern in S. lividans, the natural host for KcsA, warrants further investigation. In conclusion, our study sheds light on the role of lipids in stabilizing KcsA and highlights the importance of specific lipid-protein interactions in modulating its conformational states.

The Lipidomics Reporting Checklist A framework for transparency of lipidomic experiments and repurposing resource data.

The rapid increase in lipidomic studies has led to a collaborative effort within the community to establish standards and criteria for producing, documenting, and disseminating data. Creating a dynamic easy-to-use checklist that condenses key information about lipidomic experiments into common terminology will enhance the field's consistency, comparability, and repeatability. Here, we describe the structure and rationale of the established Lipidomics Minimal Reporting Checklist to increase transparency in lipidomics research.

Deciphering ApoE Genotype-Driven Proteomic and Lipidomic Alterations in Alzheimer's Disease Across Distinct Brain Regions.

Alzheimer's disease (AD) is a neurodegenerative disease with a complex etiology influenced by confounding factors such as genetic polymorphisms, age, sex, and race. Traditionally, AD research has not prioritized these influences, resulting in dramatically skewed cohorts such as three times the number of Apolipoprotein E (APOE) ε4-allele carriers in AD relative to healthy cohorts. Thus, the resulting molecular changes in AD have previously been complicated by the influence of apolipoprotein E disparities. To explore how apolipoprotein E polymorphism influences AD progression, 62 post-mortem patients consisting of 33 AD and 29 controls (Ctrl) were studied to balance the number of ε4-allele carriers and facilitate a molecular comparison of the apolipoprotein E genotype. Lipid and protein perturbations were assessed across AD diagnosed brains compared to Ctrl brains, ε4 allele carriers (APOE4+ for those carrying 1 or 2 ε4s and APOE4- for non-ε4 carriers), and differences in ε3ε3 and ε3ε4 Ctrl brains across two brain regions (frontal cortex (FCX) and cerebellum (CBM)). The region-specific influences of apolipoprotein E on AD mechanisms showcased mitochondrial dysfunction and cell proteostasis at the core of AD pathophysiology in the post-mortem brains, indicating these two processes may be influenced by genotypic differences and brain morphology.

Assessing Per- and Polyfluoroalkyl Substances in Fish Fillet Using Non-Targeted Analyses.

Per- and polyfluoroalkyl substances (PFAS) are a class of thousands of man-made chemicals that are persistent and highly stable in the environment. Fish consumption has been identified as a key route of PFAS exposure for humans. However, routine fish monitoring targets only a handful of PFAS, and non-targeted analyses have largely only evaluated fish from heavily PFAS-impacted waters. Here, we evaluated PFAS in fish fillets from recreational and drinking water sources in central North Carolina to assess whether PFAS are present in these fillets that would not be detected by conventional targeted methods. We used liquid chromatography, ion mobility spectrometry, and mass spectrometry (LC-IMS-MS) to collect full scan feature data, performed suspect screening using an in-house library of 100 PFAS for high confidence feature identification, searched for additional PFAS features using non-targeted data analyses, and quantified perfluorooctanesulfonic acid (PFOS) in the fillet samples. A total of 36 PFAS were detected in the fish fillets, including 19 that would not be detected using common targeted methods, with a minimum of 6 and a maximum of 22 in individual fish. Median fillet PFOS levels were concerningly high at 11.6 to 42.3 ppb, and no significant correlation between PFOS levels and number of PFAS per fish was observed. Future PFAS monitoring in this region should target more of these 36 PFAS, and other regions not considered heavily PFAS contaminated should consider incorporating non-targeted analyses into ongoing fish monitoring studies.

From big data to big insights: statistical and bioinformatic approaches for exploring the lipidome.

The goal of lipidomic studies is to provide a broad characterization of cellular lipids present and changing in a sample of interest. Recent lipidomic research has significantly contributed to revealing the multifaceted roles that lipids play in fundamental cellular processes, including signaling, energy storage, and structural support. Furthermore, these findings have shed light on how lipids dynamically respond to various perturbations. Continued advancement in analytical techniques has also led to improved abilities to detect and identify novel lipid species, resulting in increasingly large datasets. Statistical analysis of these datasets can be challenging not only because of their vast size, but also because of the highly correlated data structure that exists due to many lipids belonging to the same metabolic or regulatory pathways. Interpretation of these lipidomic datasets is also hindered by a lack of current biological knowledge for the individual lipids. These limitations can therefore make lipidomic data analysis a daunting task. To address these difficulties and shed light on opportunities and also weaknesses in current tools, we have assembled this review. Here, we illustrate common statistical approaches for finding patterns in lipidomic datasets, including univariate hypothesis testing, unsupervised clustering, supervised classification modeling, and deep learning approaches. We then describe various bioinformatic tools often used to biologically contextualize results of interest. Overall, this review provides a framework for guiding lipidomic data analysis to promote a greater assessment of lipidomic results, while understanding potential advantages and weaknesses along the way.

Recent additions and access to a multidimensional lipidomic database containing liquid chromatography, ion mobility spectrometry, and tandem mass spectrometry information.

The importance of lipids in biology continues to grow with their recent linkages to more diseases and conditions, microbiome fluctuations, and environmental exposures. These associations have motivated researchers to evaluate lipidomic changes in numerous matrices and studies. Lipidomic analyses, however, present numerous challenges as lipid species have broad chemistries that require different extraction methods and instrumental analyses to evaluate and separate their many isomers and isobars. Increasing knowledge about different lipid characteristics is therefore crucial for improving their separation and identification. Here, we present a multidimensional database for lipids analyzed on a platform combining reversed-phase liquid chromatography, drift tube ion mobility spectrometry, collision-induced dissociation, and mass spectrometry (RPLC-DTIMS-CID-MS). This platform and the different separation characteristics it provides enables more confident lipid annotations when compared to traditional tandem mass spectrometry platforms, especially when analyzing highly isomeric molecules such as lipids. This database expands on our previous publication containing only human plasma and bronchoalveolar lavage fluid lipids and provides experimental RPLC retention times, IMS collision cross section (CCS) values, and m/z information for 877 unique lipids from additional biofluids and tissues. Specifically, the database contains 1504 precursor [M + H]+, [M + NH4]+, [M + Na]+, [M-H]-, [M-2H]2-, [M + HCOO]-, and [M + CH3COO]- ion species and their associated CID fragments which are commonly targeted in clinical and environmental studies, in addition to being present in the chloroform layer of Folch extractions. Furthermore, this multidimensional RPLC-DTIMS-CID-MS database spans 5 lipid categories (fatty acids, sterols, sphingolipids, glycerolipids, and glycerophospholipids) and 24 lipid classes. We have also created a webpage (tarheels.live/bakerlab/databases/) to enhance the accessibility of this resource which will be populated regularly with new lipids as we identify additional species and integrate novel standards.

Aggregated Molecular Phenotype Scores: Enhancing Assessment and Visualization of Mass Spectrometry Imaging Data for Tissue-Based Diagnostics.

Mass spectrometry imaging (MSI) has gained increasing popularity for tissue-based diagnostics due to its ability to identify and visualize molecular characteristics unique to different phenotypes within heterogeneous samples. Data from MSI experiments are often assessed and visualized using various supervised and unsupervised statistical approaches. However, these approaches tend to fall short in identifying and concisely visualizing subtle, phenotype-relevant molecular changes. To address these shortcomings, we developed aggregated molecular phenotype (AMP) scores. AMP scores are generated using an ensemble machine learning approach to first select features differentiating phenotypes, weight the features using logistic regression, and combine the weights and feature abundances. AMP scores are then scaled between 0 and 1, with lower values generally corresponding to class 1 phenotypes (typically control) and higher scores relating to class 2 phenotypes. AMP scores, therefore, allow the evaluation of multiple features simultaneously and showcase the degree to which these features correlate with various phenotypes. Due to the ensembled approach, AMP scores are able to overcome limitations associated with individual models, leading to high diagnostic accuracy and interpretability. Here, AMP score performance was evaluated using metabolomic data collected from desorption electrospray ionization MSI. Initial comparisons of cancerous human tissues to their normal or benign counterparts illustrated that AMP scores distinguished phenotypes with high accuracy, sensitivity, and specificity. Furthermore, when combined with spatial coordinates, AMP scores allow visualization of tissue sections in one map with distinguished phenotypic borders, highlighting their diagnostic utility.

Using Multidimensional Separations to Distinguish Isomeric Amino Acid-Bile Acid Conjugates and Assess Their Presence and Perturbations in Model Systems.

Bile acids play key roles in nutrient uptake, inflammation, signaling, and microbiome composition. While previous bile acid analyses have primarily focused on profiling 5 canonical primary and secondary bile acids and their glycine and taurine amino acid-bile acid (AA-BA) conjugates, recent studies suggest that many other microbial conjugated bile acids (or MCBAs) exist. MCBAs are produced by the gut microbiota and serve as biomarkers, providing information about early disease onset and gut health. Here we analyzed 8 core bile acids synthetically conjugated with 22 proteinogenic and nonproteogenic amino acids totaling 176 MCBAs. Since many of the conjugates were isomeric and only 42 different m/z values resulted from the 176 MCBAs, a platform coupling liquid chromatography, ion mobility spectrometry, and mass spectrometry (LC-IMS-MS) was used for their separation. Their molecular characteristics were then used to create an in-house extended bile acid library for a combined total of 182 unique compounds. Additionally, ∼250 rare bile acid extracts were also assessed to provide additional resources for bile acid profiling and identification. This library was then applied to healthy mice dosed with antibiotics and humans having fecal microbiota transplantation (FMT) to assess the MCBA presence and changes in the gut before and after each perturbation.

Impacts of Gestational FireMaster 550 Exposure on the Neonatal Cortex Are Sex Specific and Largely Attributable to the Organophosphate Esters.

INTRODUCTION: Flame retardants (FRs) are common bodily and environmental pollutants, creating concern about their potential toxicity. We and others have found that the commercial mixture FireMaster® 550 (FM 550) or its individual brominated (BFR) and organophosphate ester (OPFR) components are potential developmental neurotoxicants. Using Wistar rats, we previously reported that developmental exposure to FM 550 or its component classes produced sex- and compound-specific effects on adult socioemotional behaviors. The underlying mechanisms driving the behavioral phenotypes are unknown. METHODS: To further mechanistic understanding, here we conducted transcriptomics in parallel with a novel lipidomics approach using cortical tissues from newborn siblings of the rats in the published behavioral study. Inclusion of lipid composition is significant because it is rarely examined in developmental neurotoxicity studies. Pups were gestationally exposed via oral dosing to the dam to FM 550 or the BFR or OPFR components at environmentally relevant doses. RESULTS: The neonatal cortex was highly sexually dimorphic in lipid and transcriptome composition, and males were more significantly impacted by FR exposure. Multiple adverse modes of action for the BFRs and OPFRs on neurodevelopment were identified, with the OPFRs being more disruptive than the BFRs via multiple mechanisms including dysregulation of mitochondrial function and disruption of cholinergic and glutamatergic systems. Disrupted mitochondrial function by environmental factors has been linked to a higher risk of autism spectrum disorders and neurodegenerative disorders. Impacted lipid classes included ceramides, sphingomyelins, and triacylglycerides. Robust ceramide upregulation in the OPFR females could suggest a heightened risk of brain metabolic disease. CONCLUSIONS: This study reveals multiple mechanisms by which the components of a common FR mixture are developmentally neurotoxic and that the OPFRs may be the compounds of greatest concern.

Assessing How Chemical Exposures Affect Human Health.

Measuring chemical exposure is extremely challenging due to the range and number of anthropogenic molecules encountered in our daily lives, as well as their complex transformations throughout the body. To broadly characterize how chemical exposures influence human health, a combination of genomic, transcriptomic, proteomic, endogenous metabolomic, and xenobiotic measurements must be performed. However, while genomic, transcriptomic, and proteomic analyses have rapidly progressed over the last two decades, advancements in instrumentation and computations for nontargeted xenobiotic and endogenous small molecule measurements are still greatly needed.

A Preprocessing Tool for Enhanced Ion Mobility-Mass Spectrometry-Based Omics Workflows.

The ability to improve the data quality of ion mobility-mass spectrometry (IM-MS) measurements is of great importance for enabling modular and efficient computational workflows and gaining better qualitative and quantitative insights from complex biological and environmental samples. We developed the PNNL PreProcessor, a standalone and user-friendly software housing various algorithmic implementations to generate new MS-files with enhanced signal quality and in the same instrument format. Different experimental approaches are supported for IM-MS based on Drift-Tube (DT) and Structures for Lossless Ion Manipulations (SLIM), including liquid chromatography (LC) and infusion analyses. The algorithms extend the dynamic range of the detection system, while reducing file sizes for faster and memory-efficient downstream processing. Specifically, multidimensional smoothing improves peak shapes of poorly defined low-abundance signals, and saturation repair reconstructs the intensity profile of high-abundance peaks from various analyte types. Other functionalities are data compression and interpolation, IM demultiplexing, noise filtering by low intensity threshold and spike removal, and exporting of acquisition metadata. Several advantages of the tool are illustrated, including an increase of 19.4% in lipid annotations and a two-times faster processing of LC-DT IM-MS data-independent acquisition spectra from a complex lipid extract of a standard human plasma sample. The software is freely available at https://omics.pnl.gov/software/pnnl-preprocessor.

Development and Application of Multidimensional Lipid Libraries to Investigate Lipidomic Dysregulation Related to Smoke Inhalation Injury Severity.

The implication of lipid dysregulation in diseases, toxic exposure outcomes, and inflammation has brought great interest to lipidomic studies. However, lipids have proven to be analytically challenging due to their highly isomeric nature and vast concentration ranges in biological matrices. Therefore, multidimensional techniques such as those integrating liquid chromatography, ion mobility spectrometry, collision-induced dissociation, and mass spectrometry (LC-IMS-CID-MS) have been implemented to separate lipid isomers as well as provide structural information and increased identification confidence. These data sets are however extremely large and complex, resulting in challenges for data processing and annotation. Here, we have overcome these challenges by developing sample-specific multidimensional lipid libraries using the freely available software Skyline. Specifically, the human plasma library developed for this work contains over 500 unique lipids and is combined with adapted Skyline functions such as indexed retention time (iRT) for retention time prediction and IMS drift time filtering for enhanced selectivity. For comparison with other studies, this database was used to annotate LC-IMS-CID-MS data from a NIST SRM 1950 extract. The same workflow was then utilized to assess plasma and bronchoalveolar lavage fluid (BALF) samples from patients with varying degrees of smoke inhalation injury to identify lipid-based patient prognostic and diagnostic markers.

Cupric Ions Selectively Modulate TRAAK-Phosphatidylserine Interactions.

TRAAK and TREK2 are two-pore domain K+ (K2P) channels and are modulated by diverse factors including temperature, membrane stretching, and lipids, such as phosphatidic acid. In addition, copper and zinc, both of which are essential for life, are known to regulate TREK2 and a number of other ion channels. However, the role of ions in the association of lipids with integral membrane proteins is poorly understood. Here, we discover cupric ions selectively modulate the binding of phosphatidylserine (PS) to TRAAK but not TREK2. Other divalent cations (Ca2+, Mg2+, and Zn2+) bind both channels but have no impact on binding PS and other lipids. Additionally, TRAAK binds more avidly to Cu2+ and Zn2+ than TREK2. In the presence of Cu2+, TRAAK binds similarly to PS with different acyl chains, indicating a crucial role of the serine headgroup in coordinating Cu2+. High-resolution native mass spectrometry (MS) enables the determination of equilibrium binding constants for distinct Cu2+-bound stoichiometries and uncovered the highest coupling factor corresponds to a 1:1 PS-to-Cu2+ ratio. Interestingly, the next three highest coupling factors had a ∼1.5:1 PS-to-Cu2+ ratio. Our findings bring forth the role of cupric ions as an essential cofactor in selective TRAAK-PS interactions.

Combining Isotopologue Workflows and Simultaneous Multidimensional Separations to Detect, Identify, and Validate Metabolites in Untargeted Analyses.

While the combination of liquid chromatography and tandem mass spectrometry (LC-MS/MS) is commonly used for feature annotation in untargeted omics experiments, ensuring these prioritized features originate from endogenous metabolism remains challenging. Isotopologue workflows, such as isotopic ratio outlier analysis (IROA), mass isotopomer ratio analysis of U-13C labeled extracts (MIRACLE), and credentialing incorporate isotopic labels directly into metabolic precursors, guaranteeing that all features of interest are unequivocal byproducts of cellular metabolism. Furthermore, comprehensive separation and annotation of small molecules continue to challenge the metabolomics field, particularly for isomeric systems. In this paper, we evaluate the analytical utility of incorporating ion mobility spectrometry (IMS) as an additional separation mechanism into standard LC-MS/MS isotopologue workflows. Since isotopically labeled molecules codrift in the IMS dimension with their 12C versions, LC-IMS-CID-MS provides four dimensions (LC, IMS, MS, and MS/MS) to directly investigate the metabolic activity of prioritized untargeted features. Here, we demonstrate this additional selectivity by showcasing how a preliminary data set of 30 endogeneous metabolites are putatively annotated from isotopically labeled Escherichia coli cultures when analyzed by LC-IMS-CID-MS. Metabolite annotations were based on several molecular descriptors, including accurate mass measurement, carbon number, annotated fragmentation spectra, and collision cross section (CCS), collectively illustrating the importance of incorporating IMS into isotopologue workflows. Overall, our results highlight the enhanced separation space and increased annotation confidence afforded by IMS for metabolic characterization and provide a unique perspective for future developments in isotopically labeled MS experiments.

High-Resolution Demultiplexing (HRdm) Ion Mobility Spectrometry-Mass Spectrometry for Aspartic and Isoaspartic Acid Determination and Screening.

Isomeric peptide analyses are an analytical challenge of great importance to therapeutic monoclonal antibody and other biotherapeutic product development workflows. Aspartic acid (Asp, D) to isoaspartic acid (isoAsp, isoD) isomerization is a critical quality attribute (CQA) that requires careful control, monitoring, and quantitation during the drug discovery and production processes. While the formation of isoAsp has been implicated in a variety of disease states such as autoimmune diseases and several types of cancer, it is also understood that the formation of isoAsp results in a structural change impacting efficacy, potency, and immunogenic properties, all of which are undesirable. Currently, lengthy ultrahigh-performance liquid chromatography (UPLC) separations are coupled with MS for CQA analyses; however, these measurements often take over an hour and drastically limit analysis throughput. In this manuscript, drift tube ion mobility spectrometry-mass spectrometry (DTIMS-MS) and both a standard and high-resolution demultiplexing approach were utilized to study eight isomeric Asp and isoAsp peptide pairs. While the limited resolving power associated with the standard DTIMS analysis only separated three of the eight pairs, the application of HRdm distinguished seven of the eight and was only unable to separate DL and isoDL. The rapid high-throughput HRdm DTIMS-MS method was also interfaced with both flow injection and an automated solid phase extraction system to present the first application of HRdm for isoAsp and Asp assessment and demonstrate screening capabilities for isomeric peptides in complex samples, resulting in a workflow highly suitable for biopharmaceutical research needs.

CCS Predictor 2.0: An Open-Source Jupyter Notebook Tool for Filtering Out False Positives in Metabolomics.

Metabolite annotation continues to be the widely accepted bottleneck in nontargeted metabolomics workflows. Annotation of metabolites typically relies on a combination of high-resolution mass spectrometry (MS) with parent and tandem measurements, isotope cluster evaluations, and Kendrick mass defect (KMD) analysis. Chromatographic retention time matching with standards is often used at the later stages of the process, which can also be followed by metabolite isolation and structure confirmation utilizing nuclear magnetic resonance (NMR) spectroscopy. The measurement of gas-phase collision cross-section (CCS) values by ion mobility (IM) spectrometry also adds an important dimension to this workflow by generating an additional molecular parameter that can be used for filtering unlikely structures. The millisecond timescale of IM spectrometry allows the rapid measurement of CCS values and allows easy pairing with existing MS workflows. Here, we report on a highly accurate machine learning algorithm (CCSP 2.0) in an open-source Jupyter Notebook format to predict CCS values based on linear support vector regression models. This tool allows customization of the training set to the needs of the user, enabling the production of models for new adducts or previously unexplored molecular classes. CCSP produces predictions with accuracy equal to or greater than existing machine learning approaches such as CCSbase, DeepCCS, and AllCCS, while being better aligned with FAIR (Findable, Accessible, Interoperable, and Reusable) data principles. Another unique aspect of CCSP 2.0 is its inclusion of a large library of 1613 molecular descriptors via the Mordred Python package, further encoding the fine aspects of isomeric molecular structures. CCS prediction accuracy was tested using CCS values in the McLean CCS Compendium with median relative errors of 1.25, 1.73, and 1.87% for the 170 [M - H]-, 155 [M + H]+, and 138 [M + Na]+ adducts tested. For superclass-matched data sets, CCS predictions via CCSP allowed filtering of 36.1% of incorrect structures while retaining a total of 100% of the correct annotations using a ΔCCS threshold of 2.8% and a mass error of 10 ppm.

Stability and Dissociation of Adeno-Associated Viral Capsids by Variable Temperature-Charge Detection-Mass Spectrometry.

Adeno-associated viral (AAV) vectors have emerged as gene therapy and vaccine delivery systems. Differential scanning fluorimetry or differential scanning calorimetry is commonly used to measure the thermal stability of AAVs, but these global methods are unable to distinguish the stabilities of different AAV subpopulations in the same sample. To address this challenge, we combined charge detection-mass spectrometry (CD-MS) with a variable temperature (VT) electrospray source that controls the temperature of the solution prior to electrospray. Using VT-CD-MS, we measured the thermal stabilities of empty and filled capsids. We found that filled AAVs ejected their cargo first and formed intermediate empty capsids before completely dissociating. Finally, we observed that pH stress caused a major decrease in thermal stability. This new approach better characterizes the thermal dissociation of AAVs, providing the simultaneous measurement of the stabilities and dissociation pathways of different subpopulations.

Utilizing Pine Needles to Temporally and Spatially Profile Per- and Polyfluoroalkyl Substances (PFAS).

As concerns over exposure to per- and polyfluoroalkyl substances (PFAS) are continually increasing, novel methods to monitor their presence and modifications are greatly needed, as some have known toxic and bioaccumulative characteristics while most have unknown effects. This task however is not simple, as the Environmental Protection Agency (EPA) CompTox PFAS list contains more than 9000 substances as of September 2020 with additional substances added continually. Nontargeted analyses are therefore crucial to investigating the presence of this immense list of possible PFAS. Here, we utilized archived and field-sampled pine needles as widely available passive samplers and a novel nontargeted, multidimensional analytical method coupling liquid chromatography, ion mobility spectrometry, and mass spectrometry (LC-IMS-MS) to evaluate the temporal and spatial presence of numerous PFAS. Over 70 PFAS were detected in the pine needles from this study, including both traditionally monitored legacy perfluoroalkyl acids (PFAAs) and their emerging replacements such as chlorinated derivatives, ultrashort chain PFAAs, perfluoroalkyl ether acids including hexafluoropropylene oxide dimer acid (HFPO-DA, "GenX") and Nafion byproduct 2, and a cyclic perfluorooctanesulfonic acid (PFOS) analog. Results from this study provide critical insight related to PFAS transport, contamination, and reduction efforts over the past six decades.

Uncovering PFAS and Other Xenobiotics in the Dark Metabolome Using Ion Mobility Spectrometry, Mass Defect Analysis, and Machine Learning.

The identification of xenobiotics in nontargeted metabolomic analyses is a vital step in understanding human exposure. Xenobiotic metabolism, transformation, excretion, and coexistence with other endogenous molecules, however, greatly complicate the interpretation of features detected in nontargeted studies. While mass spectrometry (MS)-based platforms are commonly used in metabolomic measurements, deconvoluting endogenous metabolites from xenobiotics is also often challenged by the lack of xenobiotic parent and metabolite standards as well as the numerous isomers possible for each small molecule m/z feature. Here, we evaluate a xenobiotic structural annotation workflow using ion mobility spectrometry coupled with MS (IMS-MS), mass defect filtering, and machine learning to uncover potential xenobiotic classes and species in large metabolomic feature lists. Xenobiotic classes examined included those of known high toxicities, including per- and polyfluoroalkyl substances (PFAS), polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), and pesticides. Specifically, when the workflow was applied to identify PFAS in the NIST SRM 1957 and 909c human serum samples, it greatly reduced the hundreds of detected liquid chromatography (LC)-IMS-MS features by utilizing both mass defect filtering and m/z versus IMS collision cross sections relationships. These potential PFAS features were then compared to the EPA CompTox entries, and while some matched within specific m/z tolerances, there were still many unknowns illustrating the importance of nontargeted studies for detecting new molecules with known chemical characteristics. Additionally, this workflow can also be utilized to evaluate other xenobiotics and enable more confident annotations from nontargeted studies.