Tile-Based Fisher Ratio Analysis of Comprehensive Two-Dimensional Gas Chromatography Time-of-Flight Mass Spectrometry (GC × GC–TOFMS) Data Using a Null Distribution Approach

Comprehensive two-dimensional gas chromatography as a tool for targeted and non-targeted analysis of contaminants of emerging concern in wastewater

Wastewater is a major reservoir for chemical contaminants, both anthropogenic and biogenic. Recent chemical and toxicological analysis reveals the abundance and impact of these compounds, often termed contaminants of emerging concern (CECs). Concurrently, incomplete removal of these compounds in wastewater treatment plants sets a precedent for detailed characterisation and monitoring of such substances. Although liquid chromatography (LC) is frequently used for analysis of CECs in wastewater, gas chromatography (GC) maintains its significance for non-polar to mid-polar analytes. GC offers advantages such as increased separation efficiency, fewer matrix effects, and greater availability and reliability of reference mass spectra compared to LC. Comprehensive two-dimensional gas chromatography (GC × GC) delivers unmatched peak capacity and separational capabilities, critical in the resolution of diverse compound groups present within wastewater. When coupled with high resolution mass spectrometry, it provides a powerful identification tool with spectral databases and both 1st and 2nd dimensional retention indices, and has allowed for the separation, reliable annotation and characterisation of diverse CECs within wastewater in recent years. Herein, on the basis of recent studies from the last fifteen years, we outline cutting-edge methodologies and strategies for wastewater analysis using GC × GC. This includes sample preparation, derivatization of polar analytes, instrumental setup, and data analysis, ultimately providing the reader a framework for future non-targeted analysis of wastewater and other complex environmental matrices.

Comprehensive characterization of European house dust contaminants: Concentrations and profiles, geographical variability, and implications for chemical regulation and health risk

This study investigated the concentration profiles and geographical variability of contaminants in house dust across Europe. A collaborative trial (CT) was organized by the NORMAN network using pooled dust and advanced chromatographic and mass spectrometric techniques combined with suspect screening and non-target screening (NTS). Over 1200 anthropogenic compounds were tentatively identified. Additionally, seventy-five individual samples were subjected to target analysis and NTS. The median concentrations of most contaminants varied <3-fold across Europe, and the contaminant profile of European dust was similar to that of North American dust, which was investigated in a previous CT. This similarity may be attributed to the use of similar consumer articles and building materials throughout the developed world. Multivariate data analysis revealed geographical trends in contaminant distribution, with north-south gradients across Europe. Geographical trends were more frequently found for compounds with rapid release (pharmaceuticals, personal care products, fragrances, pesticides, biocides) and smoke-related compounds. The concentrations of chlorinated paraffins, polycyclic aromatic hydrocarbons (PAHs), perfluorinated alkyl substances and stimulants generally increased from north to south, whereas the biocides levels decreased from north to south. Despite widespread presence of in-use contaminants in dusts, some of the highest risks come from compounds that have been restricted for decades or more. These include di(2-ethylhexyl) phthalate (DEHP), polychlorinated biphenyl (PCB) 118 and polybrominated diphenyl ethers 47, 99, and 153. DEHP remains the most abundant contaminant in European house dust, while the other compounds are classified as persistent organic pollutants (POPs). Moreover, there is a striking lack of reliable toxicity data, particularly for emerging compounds. For instance, although acceptable daily intakes (ADIs) were examined for 202 compounds, only 46 had consensus-based ADI values. The results highlight the need for proactive measures to prevent hazardous chemicals from entering the market and for careful selection of substitute chemicals, when such are needed, to avoid regrettable substitutions.

Simultaneous discovery of compounds dominated by either molding kinetics or geographical region of origin for moisture damaged cacao beans using orthogonally applied tile-based fisher ratio analysis of GC×GC-TOFMS data

Herein, two "orthogonal" characteristics of moisture damaged cacao beans (temporally dependent molding kinetics versus the time-independent geographical region of origin) are simultaneously analyzed in a comprehensive two-dimensional (2D) gas chromatography time-of-flight mass spectrometry (GC×GC-TOFMS) dataset using tile-based Fisher ratio (F-ratio) analysis. Cacao beans from six geographical regions were analyzed once a day for six days following the initiation of moisture damage to trigger the molding process. Thus, there are two "extremes" to the experimental sample class design: six time points for the molding kinetics versus the six geographical regions of origin, resulting in a 6 × 6 element signal array referred to as a composite chemical fingerprint (CCF) for each analyte. Usually, this study would involve initial generation of two separate hit lists using F-ratio analysis, one hit list from inputting the data with the six time point classes, then another hit list from inputting the dataset from the perspective of geographic region of origin. However, analysis of two separate hit lists with the intent to distill them down to one hit list is extremely time-consuming and fraught with shortcomings due to the challenges associated with attempting to match analytes across two hit lists. To address this challenge, tile-based F-ratio analysis is "orthogonally applied" to each analyte CCF to simultaneously determine two F-ratios at the chromatographic 2D location (F-ratiokinetic and F-ratioregion) for each hit, by ranking a single hit list using the higher of the two F-ratios resulting in the discovery of 591 analytes. Further, using a pseudo-null distribution approach, at the 99.9% threshold over 400 analytes were deemed suitable for PCA classification. Using a more stringent 99.999% threshold, over 100 analytes were explored more deeply using PARAFAC to provide a purified mass spectrum.

Characterization of semi-volatile compounds in 56 Italian ciders using GC×GC-TOF-MS and multivariate analysis

Fifty-six samples of differently produced commercial Italian ciders were analysed for semi-volatile organic compounds (SVOCs) profiling, using comprehensive two-dimensional gas chromatography coupled to mass spectrometry (GC×GC-TOF-MS) technique for the very first time. To properly support the compositional investigation of this emerging beverage, a chemometric approach through Principal Component Analysis (PCA) was employed. This revealed a sample distribution in agreement with results of the sensory tasting panel performed on such ciders, highlighting an excellent correlation between variables and perceived odorants. In particular, the positions of peculiar and anomalous objects in the Principal Components (PCs) space are explicitly evaluated, exploring the associated loadings (i.e., the importance of the identified chemical compounds), paying attention to their biochemical origin along the cider-making process and their impact on the sample olfactory analysis. Besides this, the t-distributed Stochastic Neighbor Embedding (t-SNE) method was shown to be an efficient tool for gathering pear ciders from the other samples (apple ciders), better than PCA. This study stands for the first survey on Italian commercial craft cider, and its results are aimed to be a milestone for its characterization and to start and promote cider culture in this country.

Volatile organic compounds produced during postmortem processes can be linked via chromatographic profiles to individual postmortem bacterial species

Decomposition odor is produced during postmortem mammalian tissue breakdown by bacteria, insects, and intrinsic chemical processes. Past research has not thoroughly investigated which volatile organic compounds (VOCs) can be linked directly to individual bacterial species on decomposing remains. The purpose of this study was to profile the VOCs produced over time by individual species of bacteria using comprehensive two-dimensional gas chromatography (GC×GC) to expand our foundational knowledge of what each bacterial species contributes to decomposition odor. Five different species of bacteria (Bacillus subtilis, Ignatzschineria indica, Ignatzschineria ureiclastica, Curtobacterium luteum, and Vagococcus lutrae) were cultured on standard nutrient agar individually and monitored daily using solid phase microextraction arrow (SPME Arrow) and GC×GC in combination with quadrupole mass spectrometry (qMS) and flame ionization detection (FID). The GC×GC-qMS/FID approach was used to generate rich VOC profiles that represented the bacterial species' metabolic VOC production longitudinally. The data obtained from the chromatographic output was used to compare with a prior study using one-dimensional GC-qMS, and also between each of the five species to investigate the extent of overlap between species. No single VOC could be found in all five bacterial species investigated, and there was little overlap in the profile between species. To further visualize these differences, chromatographic peak data was investigated using two different ordination strategies, principal component analysis (PCA) and principal coordinate analysis (PCoA). The two ordination strategies were compared with each other using a Procrustes analysis. This was performed to understand differences in ordination strategies between the separation science community and chemical ecological community. Overall, ordination strategies were found to produce similar results, as evidenced by the correlation of PCA and PCoA in the Procrustes analysis. All analysis strategies yielded distinct VOC profiles for each species. Further study of additional species will support understanding of the holistic view of decomposition odor from a chemical ecology perspective, and further support our understanding of the production of decomposition odor that culminates from such a complex environment.

Region of interest selection in heterogeneous digital image: Wine age prediction by comprehensive two-dimensional gas chromatography

This study integrates genetic algorithm (GA) with partial least squares regression (PLSR) and various variable selection methods to identify impactful regions of interest (ROI) in heterogeneous 2D chromatogram images for predicting wine age. As wine quality and aroma evolve over time, transitioning from youthful fruitiness to mature, complex flavors, which leads to alterations in the composition of essential aroma-contributing compounds. Chromatograms are segmented into subimages, and the GA-PLSR algorithm optimizes combinations based on grayscale, red-green-blue (RGB), and hue-saturation-value (HSV) histograms. The selected subimage histograms are further refined through interval selection, highlighting the compounds with the most significant influence on wine aging. Experimental validation involving 38 wine samples demonstrates the effectiveness of this approach. Cross-validation reduces the PLS model error from 2.8 to 2.4 years within a 10 × 10 subset, and during prediction, the error decreases from 2.5 to 2.3 years. The study presents a novel approach utilizing the selection of ROI for efficient processing of 2D chromatograms focusing on predicting wine age.

Global Urinary Volatolomics with (GC×)GC-TOF-MS

Urinary volatolomics offers a noninvasive approach for disease detection and monitoring. Herein we present an improved methodology for global volatolomic profiling. Wide coverage was achieved by utilizing a multiphase sorbent for volatile organic compound (VOC) extraction. A single, midpolar column gas chromatography (GC) assay yielded substantially higher numbers of monitored VOCs compared to our previously reported single-sorbent method. Multidimensional GC (GC×GC) enhanced further biomarker discovery while data analysis was simplified by using a tile-based approach. At the same time, the required urine volume was reduced 5-fold from 2 to 0.4 mL. The applicability of the methodology was demonstrated in a pancreatic ductal adenocarcinoma cohort where previous findings were confirmed while a series of additional VOCs with diagnostic potential were discovered.

Addressing new chemicals of emerging concern (CECs) in an indoor office

Indoor pollutants change over time and place. Exposure to hazardous organics is associated with adverse health effects. This work sampled gaseous organics by Tenax TA tubes in two indoor rooms, i.e., an office set as samples, and the room of chassis dynamometer (RCD) set as backgrounds. Compounds are analyzed by a thermal desorption comprehensive two-dimensional gas chromatography-quadrupole mass spectrometer (TD-GC × GC-qMS). Four new chemicals of emerging concern (CECs) are screened in 469 organics quantified. We proposed a three-step pipeline for CECs screening utilizing GC × GC including 1) non-target scanning of organics with convincing molecular structures and quantification results, 2) statistical analysis between samples and backgrounds to extract useful information, and 3) pixel-based property estimation to evaluate the contamination potential of addressed chemicals. New CECs spotted in this work are all intermediate volatility organic compounds (IVOCs), containing mintketone, isolongifolene, β-funebrene, and (5α)-androstane. Mintketone and sesquiterpenes may be derived from the use of volatile chemical products (VCPs), while (5α)-androstane is probably human-emitted. The occurrence and contamination potential of the addressed new CECs are reported for the first time. Non-target scanning and the measurement of IVOCs are of vital importance to get a full glimpse of indoor organics.

Enabling cuboid-based fisher ratio analysis using total-transfer comprehensive three-dimensional gas chromatography with time-of-flight mass spectrometry

Comprehensive three-dimensional (3D) gas chromatography with time-of-flight mass spectrometry (GC3-TOFMS) is a promising instrumental platform for the separation of volatiles and semi-volatiles due to its increased peak capacity and selectivity relative to comprehensive two-dimensional gas chromatography with TOFMS (GC×GC-TOFMS). Given the recent advances in GC3-TOFMS instrumentation, new data analysis methods are now required to analyze its complex data structure efficiently and effectively. This report highlights the development of a cuboid-based Fisher ratio (F-ratio) analysis for supervised, non-targeted studies. This approach builds upon the previously reported tile-based F-ratio software for GC×GC-TOFMS data. Cuboid-based F-ratio analysis is enabled by constructing 3D cuboids within the GC3-TOFMS chromatogram and calculating F-ratios for every cuboid on a per-mass channel basis. This methodology is evaluated using a GC3-TOFMS data set of jet fuel spiked with both non-native and native components. The neat and spiked jet fuels were collected on a total-transfer (100 % duty cycle) GC3-TOFMS instrument, employing thermal modulation between the first (1D) and second dimension (2D) columns and dynamic pressure gradient modulation between the 2D and third dimension (3D) columns. In total, cuboid-based F-ratio analysis discovered 32 spiked analytes in the top 50 hits at concentration ratios as low as 1.1. In contrast, tile-based F-ratio analysis of the corresponding GC×GC-TOFMS data only discovered 28 of the spiked analytes total, with only 25 of them in the top 50 hits. Along with discovering more analytes, cuboid-based F-ratio analysis of GC3-TOFMS data resulted in fewer false positives. The increased discoverability is due to the added peak capacity and selectivity provided by the 3D column with GC3-TOFMS resulting in improved chromatographic resolution.

A non-target geographical origin screening of botrytized wines through comprehensive two-dimensional gas chromatography coupled with high-resolution mass spectrometry

One of the most effective methods for gaining insight into the composition of trace-level volatile organic characteristics of wine products is through the use of a comprehensive two-dimensional gas chromatography-high resolution mass spectrometry (GC × GC-HRMS) technique. The vast amount of data generated by this method, however, can often be overwhelming requiring exhaustive and time-consuming analysis to identify significant statistical characteristics. The use of advanced chemometric software can achieve the same or even higher efficiency. This study aimed to identify differences based on geographical locations by analyzing the volatile organic compounds in the composition of botrytized wines from Slovakia, Hungary, France, and Austria. The volatile organic compounds were extracted by solid-phase microextraction and analyzed using GC × GC-HRMS. The data obtained from the analysis underwent Fisher-ratio (F-ratio) tile-based analysis to identify statistically significant differences. Principal component analysis demonstrated a significant distinction between wine samples based on geographical location, using only 10 statistically significant features with the highest F-ratio. In the samples, the following compounds were analyzed: methyl-octadecanoate, 2-cyanophenyl-β-phenylpropionate, α-ionone, n-octanoic acid, 1,2-dihydro-1,1,6-trimethyl-naphthalene, methyl-hexadecanoate, ethyl-pentadecanoate, ethyl-decanoate, and γ-nonalactone. These, all play an important role in cluster pattern observed on principal component analysis results. Additionally, hierarchical cluster analysis confirmed this.

Tile-Based Random Forest Analysis for Analyte Discovery in Balanced and Unbalanced GC × GC-TOFMS Data Sets

In this study, we introduce a new nontargeted tile-based supervised analysis method that combines the four-grid tiling scheme previously established for the Fisher ratio (F-ratio) analysis (FRA) with the estimation of tile hit importance using the machine learning (ML) algorithm Random Forest (RF). This approach is termed tile-based RF analysis. As opposed to the standard tile-based F-ratio analysis, the RF approach can be extended to the analysis of unbalanced data sets, i.e., different numbers of samples per class. Tile-based RF computes out-of-bag (oob) tile hit importance estimates for every summed chromatographic signal within each tile on a per-mass channel basis (m/z). These estimates are then used to rank tile hits in a descending order of importance. In the present investigation, the RF approach was applied for a two-class comparison of stool samples collected from omnivore (O) subjects and stored using two different storage conditions: liquid (Liq) and lyophilized (Lyo). Two final hit lists were generated using balanced (8 vs Eight comparison) and unbalanced (8 vs Nine comparison) data sets and compared to the hit list generated by the standard F-ratio analysis. Similar class-distinguishing analytes (p < 0.01) were discovered by both methods. However, while the FRA discovered a more comprehensive hit list (65 hits), the RF approach strictly discovered hits (31 hits for the balanced data set comparison and 29 hits for the unbalanced data set comparison) with concentration ratios, [OLiq]/[OLyo], greater than 2 (or less than 0.5). This difference is attributed to the more stringent feature selection process used by the RF algorithm. Moreover, our findings suggest that the RF approach is a promising method for identifying class-distinguishing analytes in settings characterized by both high between-class variance and high within-class variance, making it an advantageous method in the study of complex biological matrices.

Augmented visualization by computer vision and chromatographic fingerprinting on comprehensive two-dimensional gas chromatographic patterns: Unraveling diagnostic signatures in food volatilome

Computer Vision is an approach of Artificial Intelligence (AI) that conceptually enables "computers and systems to derive useful information from digital images" giving access to higher-level information and "take actions or make recommendations based on that information". Comprehensive two-dimensional chromatography gives access to highly detailed, accurate, yet unstructured information on the sample's chemical composition, and makes it possible to exploit the AI concepts at the data processing level (e.g., by Computer Vision) to rationalize raw data explorations. The goal is the understanding of the biological phenomena interrelated to a specific/diagnostic chemical signature. This study introduces a novel workflow for Computer Vision based on pattern recognition algorithms (i.e., combined untargeted and targeted UT fingerprinting) which includes the generation of composite Class Images for representative samples' classes, their effective re-alignment and registration against a comprehensive feature template followed by Augmented Visualization by comparative visual analysis. As an illustrative application, a sample set originated from a Research Project on artisanal butter (from raw sweet cream to ripened butter) is explored, capturing the evolution of volatile components along the production chain and the impact of different microbial cultures on the finished product volatilome. The workflow has significant advantages compared to the classical one-step pairwise comparison process given the ability to realign and pairwise compare both targeted and untargeted chromatographic features belonging to Class Images resembling chemical patterns from many different samples with intrinsic biological variability.

PARAFAC2×N: Coupled decomposition of multi-modal data with drift in N modes

Analysis of GC×GC-TOFMS data for large numbers of poorly-resolved peaks, and for large numbers of samples remains an enduring problem that hinders the widespread application of the technique. For multiple samples, GC×GC-TOFMS data for specific chromatographic regions manifests as a 4^th order tensor of I mass spectral acquisitions, J mass channels, K modulations, and L samples. Chromatographic drift is common along both the first-dimension (modulations), and along the second-dimension (mass spectral acquisitions), while drift along the mass channel is for all practical purposes nonexistent. A number of solutions to handling GC×GC-TOFMS data have been proposed: these involve reshaping the data to make it amenable to either 2^nd order decomposition techniques based on Multivariate Curve Resolution (MCR), or 3^rd order decomposition techniques such as Parallel Factor Analysis 2 (PARAFAC2). PARAFAC2 has been utilised to model chromatographic drift along one mode, which has enabled its use for robust decomposition of multiple GC-MS experiments. Although extensible, it is not straightforward to implement a PARAFAC2 model that accounts for drift along multiple modes. In this submission, we demonstrate a new approach and a general theory for modelling data with drift along multiple modes, for applications in multidimensional chromatography with multivariate detection. The proposed model captures over 99.9% of variance for a synthetic data set, presenting an extreme example of peak drift and co-elution across two modes of separation.

Investigating analyte breakthrough under non-linear isotherm conditions during solid phase extraction facilitated by non-targeted analysis with comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry

Solid phase extraction (SPE) sample preparation for the analysis of complex organic mixtures is often applied assuming all analytes of interest will preconcentrate on the stationary phase. This assumption ignores the reality that extraction is a dynamic interactive process and a diverse range of affinities for the stationary phase will result in equally diverse breakthrough volumes due to competitive sorption processes. To study this dynamic interactive process, and further to take advantage of it, we extracted a JP-8 jet fuel spiked with 40 ppm of a polar compound mix with silica and alumina SPE cartridges and analyzed sequential extracted fractions of the fuel to both assess the shifting chemical landscape present in the extraction and the impact of both SPE stationary phases on this process. Tile-based 1v1 comparative analysis (a recently reported extension of tile-based Fisher ratio analysis) was used to discover the (polar) compounds whose concentrations change between extracted fractions, discovering 21 compounds extracted with silica and 27 compounds extracted with alumina with at least a 2-fold change in concentration from the neat sample relative to the first 1 mL pass fraction sample. These compounds were quantified in each fraction to construct concentration ratio profiles, defined as the concentration ratio for a given SPE fraction per analyte compound relative to the analyte concentration in the neat fuel, for which the extraction behavior for each analyte could be assessed. These analyte compounds were found to breakthrough at different rates, with some analytes remaining on the column indefinitely (until extracted with a subsequent polar solvent) and other analytes eluting before the extraction is complete. Furthermore, in a comparison of the effect of selected stationary phase, alumina was found to retain oxygen-containing phenolic compounds to a greater extent than silica. Principal component analysis (PCA) was used to analyze the concentration ratio profiles of the various trace analytes in the JP8 fuel (phenols, indoles, etc.) in the context of their stationary phase affinity (silica or alumina) and competitive sorption behavior.

Enhancing partial least squares modeling of comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry data by tile-based variance ranking

Chemometric methods like partial least squares (PLS) regression are valuable for correlating sample-based differences hidden in comprehensive two-dimensional gas chromatography (GC × GC) data to independently measured physicochemical properties. Herein, this work establishes the first implementation of tile-based variance ranking as a selective data reduction methodology to improve PLS modeling performance of 58 diverse aerospace fuels. Tile-based variance ranking discovered a total of 521 analytes with a square of the relative standard deviation (RSD²) in signal between 0.07 to 22.84. The goodness-of-fit for the models were determined by their normalized root-mean-square error of cross-validation (NRMSECV) and normalized root-mean-square error of prediction (NRMSEP). PLS models developed for viscosity, hydrogen content, and heat of combustion using all 521 features discovered by tile-based variance ranking had a respective NRMSECV (NRMSEP) equal to 10.5 % (10.2 %), 8.3 % (7.6 %), and 13.1 % (13.5 %). In contrast, use of a single-grid binning scheme, a common data reduction strategy for PLS analysis, resulted in less accurate models for viscosity (NRMSECV = 14.2 %; NRMSEP = 14.3 %), hydrogen content (NRMSECV = 12.1 %; NRMSEP = 11.0 %), and heat of combustion (NRMSECV = 14.4 %; NRMSEP = 13.6 %). Further, the features discovered by tile-based variance ranking can be optimized for each PLS model with RReliefF analysis, a machine learning algorithm. RReliefF feature optimization selected 48, 125, and 172 analytes out of the original 521 discovered by tile-based variance ranking to model viscosity, hydrogen content, and heat of combustion, respectively. The RReliefF optimized features developed highly accurate property-composition models for viscosity (NRMSECV = 7.9 %; NRMSEP = 5.8 %), hydrogen content (NRMSECV = 7.0 %; NRMSEP = 4.9 %), heat of combustion (NRMSECV = 7.9 %; NRMSEP = 8.4 %). This work also demonstrates that processing the chromatograms with a tile-based approach allows the analyst to directly identify the analytes of importance in a PLS model. Coupling tile-based feature selection with PLS analysis allows for deeper understanding in any property-composition study.

Evaluation of different getter substrates using two-dimensional gas chromatography with time of flight mass spectrometry

While understanding hydrogen uptake by organic based getters such as 1,4-bis(phenylethynyl)benzene (DEB) combined with a palladium(0)bis(dibenzylideneacetone) (Pd(dba)₂) catalyst is essential, another crucial element to understand is the decomposition of the DEB, Pd(dba)₂, and/or substrate material. The breakdown of these materials may create unwanted volatiles, which may interact with and lead to deterioration of sensitive materials. Moreover, it is critical to understand if different substrates cause the getter and/or catalyst to degrade in different manners. Utilizing comprehensive two-dimensional gas chromatography (GC×GC) with time-of-flight mass spectrometry (TOFMS), the presence of volatiles located in the headspace of various DEB/Pd(dba)₂ getter substrates is examined. These samples include a getter infused silicone foam, a hydrogenated getter infused silicone foam, an activated carbon getter pellet, and a hydrogenated activated carbon getter pellet. Application of Fisher ratio (F-ratio) analyses lead to the identification of several compounds that are generated or consumed through the hydrogenation process. These include benzene derivatives such as bibenzyl, benzaldehyde, and vinyl benzoate in the activated carbon pellets and 1,5-diphenyl-3-pentanone, toluene, styrene, and 1-1'(2-pentene 1,5-diyl)bis benzene in the silicone foams, and alkane/alkene derivatives such undecane, 4-tridecene, and decane in the activated carbon pellets and 2,6-dimethyl undecane in the silicone foams. Further comparison of the different hydrogenated getter substrates (e.g. activated carbon pellet and silicone foam) indicates that the different substrates alter the decomposition products created from the degradation of the DEB and Pd(dba)₂.

Optimization of Parameters for ROI Data Compression for Nontargeted Analyses Using LC-HRMS

Nontargeted analyses of low-concentration analytes in the information-rich data collected by liquid chromatography with high-resolution mass spectrometry detection can be challenging to accomplish in an efficient and comprehensive manner. The aim of this study is to demonstrate a workflow involving targeted parameter optimization for entire chromatograms using region of interest (ROI) data compression uncoupled from a subsequent tile-based Fisher ratio (F-ratio) analysis, a supervised discovery-based method, for the discovery of low-concentration analytes. Soil samples spiked with 18 pesticides at nominal concentrations ranging from 0.1 to 50 ppb for a total of six sample classes served as challenging samples to demonstrate the overall workflow. Optimization of two parameters proved to be the most critical for ROI data compression: the signal threshold parameter and the admissible mass deviation parameter. The parameter optimization method workflow we introduce is based upon spiking known analytes into a representative sample and determining the number of detectable spikes and the Δppm for various combinations of the signal threshold and admissible mass deviation, where Δppm is the absolute value of the difference between the theoretical m/z and the ROI m/z. Once optimal parameters are determined providing the lowest average Δppm and the greatest number of detectable analytes, the optimized parameters can be utilized for the intended analysis. Herein, tile-based F-ratio analysis was performed on the ROI compressed data of all spiked soil samples first by applying ROI parameters recommended in the literature, referred to herein as the initial ROI parameters, and finally by the combination of the two optimized parameters. Using the initial ROI parameters, three pesticides were discovered, whereas all 18 spiked pesticides were discovered by optimizing both ROI parameters.

Radiolytic degradation of dodecane substituted with common energetic functional groups†

Explosives exist in and are expected to withstand a variety of harsh environments up to and including ionizing radiation, though little is known about the chemical consequences of exposing explosives to an ionizing radiation field. This study focused on the radiation-induced chemical changes to a variety of common energetic functional groups by utilizing a consistent molecular backbone. Dodecane was substituted with azide, nitro, nitrate ester, and nitramine functional groups and γ-irradiated with ⁶⁰Co in order to study how the functional group degraded along with what the relative stability to ionizing radiation was. Chemical changes were assessed using a combination of analysis techniques including: nuclear magnetic resonance (NMR) spectroscopy, gas chromatography of both the condensed and gas phases, Raman spectroscopy, and Fourier transform infrared (FTIR) spectroscopy. Results revealed that much of the damage to the molecules was on the energetic functional group and often concentrated on the trigger linkage, also known as the weakest bond in the molecule. The general trend from most to least susceptible to radiolytic damage was found to be D–ONO₂ → D–N₃ → D–NHNO₂ → D–NO₂. These results also appear to be in line with the relative stability of these functional groups to things such as photolysis, thermolysis, and explosive insults.

The relative radiolytic stability of dodecane functionalized with common energetic functional groups was explored with gamma irradiation and probed by various analytical techniques.

Fisher ratio feature selection by manual peak area calculations on comprehensive two-dimensional gas chromatography data using standard mixtures with variable composition, storage, and interferences

Comprehensive two-dimensional gas chromatography (GC×GC) is becoming increasingly more common for non-targeted characterization of complex volatile mixtures. The information gained with higher peak capacity and sensitivity provides additional sample composition information when one-dimensional GC is not adequate. GC×GC generates complex multivariate data sets when using non-targeted analysis to discover analytes. Fisher ratio (FR) analysis is applied to discern class markers, limiting complex GC×GC profiles to the most discriminating compounds between classes. While many approaches for feature selection using FR analysis exist, FR can be calculated relatively easily directly on peak areas after any native software has performed peak detection. This study evaluated the success rates of manual FR calculation and comparison to a critical F-value for samples analyzed by GC×GC with defined concentration differences. Long-term storage of samples and other spiked interferences were also investigated to examine their impact on analyzing mixtures using this FR feature selection strategy. Success rates were generally high with mostly 90-100% success rates and some instances of percentages between 80 and 90%. There were rare cases of false positives present and a low occurrence of false negatives. When errors were made in the selection of a compound, it was typically due to chromatographic artifacts present in chromatograms and not from the FR approach itself. This work provides foundational experimental data on the use of manual FR calculations for feature selection from GC×GC data.

Application of clusterization algorithms for analysis of semivolatile pollutants in Arkhangelsk snow

The best way to understand the environmental status of a certain region involves thorough non-target analysis, which will result in a list of pollutants under concern. Arkhangelsk (64° 32' N 40° 32' E, pop. ~ 344,000) is the largest city in the world to the north of the 60^th parallel. Several industrial enterprises and the "cold finger" effect represent the major sources of air contamination in the city. Analysis of snow with comprehensive two-dimensional gas chromatography-high-resolution mass spectrometry allows detecting and quantifying the most hazardous volatile and semivolatile anthropogenic pollutants and estimating long-term air pollution. Target analysis, suspect screening, and non-target analysis of snow samples collected from ten sites within the city revealed the presence of several hundreds of organic compounds including 18 species from the US EPA list of priority pollutants. Fortunately, the levels of these compounds appeared to be much lower than the safe levels established in Russia. Phenol and dioctylphthalate could be considered as the pollutants of concern because their levels were about 20% of the safe thresholds. ChromaTOF® Tile, MetaboAnalyst software platform, and open-source software protocols were applied to process the obtained data. The obtained clusterization results of the samples were generally similar for various tools; however, each of them had certain peculiarities. Bis(2-ethylhexyl) hexanedioate, benzyl alcohol, phthalates, aniline, dinitrotoluenes, and fluoranthene showed the strongest influence on the clusterization of the studied samples. Possible sources of the major pollutants were proposed: car traffic and pulp and paper mills.

Using solid-phase extraction to facilitate a focused tile-based Fisher ratio analysis of comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry data: comparative analysis of aerospace fuel composition

Tile-based Fisher ratio (F-ratio) analysis of comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry (GC × GC-TOFMS) data is a powerful, supervised discovery methodology for pinpointing sample class-distinguishing analytes between two or more sample classes. Herein, we extend this analytical methodology to focus upon specific chemical groups in kerosene-based aerospace fuel using solid-phase extraction (SPE). Treating samples with SPE removes specific compounds depending on the SPE stationary phase (i.e., silica), creating an altered "pass" sample, identical to the original "neat" sample except for the extracted compounds. Application of F-ratio analysis to the neat samples against the pass samples provides global discovery with a numerically sorted hit list of all analytes affected by the SPE procedure. Sections of GC × GC-TOFMS data from the top analyte hits are reconstructed to form a "stitch" chromatogram to visualize the sample class-distinguishing compounds, revealing excellent agreement with the extract chromatogram. Additionally, utilizing the four-grid tiling scheme developed for tile-based F-ratio analysis, we demonstrate a tile-based pairwise analysis method, referred to as 1v1 analysis, to discover analytes that differ in concentration between two fuel chromatograms. Application of 1v1 analysis is highly efficient since replicates do not necessarily need to be run on the GC × GC-TOFMS instrument, which is beneficial for sample-limited applications. The 1v1 analyses discovered most of the same features as F-ratio analysis, ranging from 69 to 81% of the features discovered by F-ratio analysis while requiring one-sixth the data. Lastly, the overall methodology is applied to three candidate rocket fuels to better understand the compound class-distinguishing differences. The separate hit lists produced for high-concentration bulk hydrocarbon differences and low-concentration level polar compound differences provided valuable insight into these candidate rocket fuels.

Untargeted region of interest selection for gas chromatography - mass spectrometry data using a pseudo F-ratio moving window

There are many challenges associated with analysing gas chromatography - mass spectrometry (GC-MS) data. Many of these challenges stem from the fact that electron ionization (EI) can make it difficult to recover molecular information due to the high degree of fragmentation with concomitant loss of molecular ion signal. With GC-MS data there are often many common fragment ions shared among closely-eluting peaks, necessitating sophisticated methods for analysis. Some of these methods are fully automated, but make some assumptions about the data which can introduce artifacts during the analysis. Chemometric methods such as Multivariate Curve Resolution (MCR), or Parallel Factor Analysis (PARAFAC/PARAFAC2) are particularly attractive, since they are flexible and make relatively few assumptions about the data - ideally resulting in fewer artifacts. These methods do require expert user intervention to determine the most relevant regions of interest and an appropriate number of components, k, for each region. Automated region of interest selection is needed to permit automated batch processing of chromatographic data with advanced signal deconvolution. Here, we propose a new method for automated, untargeted region of interest selection that accounts for the multivariate information present in GC-MS data to select regions of interest based on the ratio of the squared first, and second singular values from the Singular Value Decomposition (SVD) of a window that moves across the chromatogram. Assuming that the first singular value accounts largely for signal, and that the second singular value accounts largely for noise, it is possible to interpret the relationship between these two values as a probabilistic distribution of Fisher Ratios. The sensitivity of the algorithm was tested by investigating the concentration at which the algorithm can no longer pick out chromatographic regions known to contain signal. The algorithm achieved detection of features in a GC-MS chromatogram at concentrations below 10 pg on-column. The resultant probabilities can be interpreted as regions that contain features of interest.

Simulating comprehensive two-dimensional gas chromatography mass spectrometry data with realistic run-to-run shifting to evaluate the robustness of tile-based Fisher ratio analysis

Untargeted analysis of comprehensive two-dimensional (2D) gas chromatography time-of-flight mass spectrometry (GC×GC-TOFMS) data has the potential to be hindered by run-to-run retention time shifting. To address this challenge, tile-based Fisher ratio (F-ratio) analysis (FRA) has been developed, which utilizes a supervised, untargeted approach involving a chromatographic segmentation routine termed "tiling" combined with the ANOVA F-ratio statistic to discover class-distinguishing analytes while minimizing false positives arising from shifting. The tiling algorithm is designed to account for retention shifting in both separation dimensions. Although applications of FRA have been reported, there remains a need to thoroughly evaluate the robustness of FRA for different levels of run-to-run retention shifting in order to broaden the scope of its application. To this end, a novel method of simulating GC×GC-TOFMS chromatograms with realistic run-to-run shifting is presented by random generation of low-frequency "shift functions". The dimensionless retention-time precision, <δ_r>, which is four times the standard deviation in retention time normalized to the peak width-at-base is used as a key modeling variable along with the 2D chromatographic saturation, α_e,2D, and within-class relative standard deviation in peak area, RSD_wc. We demonstrate that all three of these variables operate together to impact true positive discovery. To quantify the "success" of true positive discovery, GC×GC-TOFMS datasets for various combinations of <δ_r>, α_e,2D, and RSD_wc were simulated and then analyzed by FRA using a wide range of relative tile areas (RTA), which is a dimensionless measure of tile size. Since each hit in the FRA hit list was known a priori as either a true or false positive based on the simulation inputs, receiver operating characteristic (ROC) curves were readily constructed. Then, the area under the ROC curve (AUROC) was used as a metric for discovery "success" for various combinations of the modeling variables. Based on the results of this study, recommendations for tile size selection and experimental design are provided, and further supported by comparison to previous tile-based FRA applications. For instance, values for <δ_r>, α_e,2D, and RSD_wc obtained from a GC×GC-TOFMS dataset of yeast metabolites suggested an optimum RTA of 6.25, corresponding closely to the RTA of 4.00 employed in the study, implying the simulation results obtained here can be generalized to real datasets.

New Perspectives on Comparative Analysis for Comprehensive Two-Dimensional Gas Chromatography

Because of the growing number of analysis scenarios involving complex samples, comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry (GC×GC–TOF-MS) is now a prominent technique for characterization. However, the limitations on time, expenses, and sample quantities, as well as the need for specialized expertise in comparative analysis, can prevent the discovery of analytes that distinguish multiple samples. This article provides an overview of the development and current status of comparative analysis for GC×GC–TOF-MS data and how key limitations can be overcome with a novel tile-based pairwise analysis method.

Profiling Olefins in Gasoline by Bromination Using GC×GC-TOFMS Followed by Discovery-Based Comparative Analysis

An analytical workflow for the analysis of olefins in gasoline that combines selective bromination and comprehensive two-dimensional (2D) gas chromatography time-of-flight mass spectrometry (GC×GC-TOFMS) with discovery-based analysis is reported. First, a standard mix containing n-alkanes, 1-alkenes, and aromatic species was brominated and quantified using % reacted as a metric for each compound class, defined as the difference in the total peak area between the brominated and original samples normalized to the original sample. The average % reacted (1 s.d.) values were -1.45% (2.8%) for the alkanes, 99.5% (0.4%) for the alkenes, and 6.7% (11.6%) for the aromatics, demonstrating excellent selectivity toward the alkenes with only minor aromatic bromination. The bromination chemistry was then applied to gasoline, followed by GC×GC-TOFMS analysis of the original and brominated gasoline. This GC×GC-TOFMS data set was then submitted to the supervised discovery tool tile-based F-ratio analysis (FRA), which reduced the large data set to only the chromatographic regions which distinguish between the original and brominated gasoline samples. FRA discovered 314 hits, 56 of which were determined statistically significant using combinatorial null distribution analysis (CNDA), a permutation-based significance test. Since the brominated olefins elute in an uncrowded region of the 2D chromatogram and have no signal in the original sample, their discoverability was greatly increased relative to the original olefins. By combining the information gained from brominated olefin standards and the structured patterns of the GC×GC separations, the top hits were identified as the dibromoalkane products of mono-olefins, with five C5 mono-olefins identified on a species level. The analytical workflow has broad implications for using selective reaction chemistries to facilitate supervised discovery by targeting desired compound classes.

Chemical Evaluation and Performance Characterization of Pentaerythritol Tetranitrate (PETN) under Melt Conditions

Pentaerythritol tetranitrate (PETN) is a nitrate ester explosive commonly used in commercial detonators. Although its degradation properties have been studied extensively, very little information has been collected on its thermal stability in the molten state due to the fact that its melting point is only ∼20 °C below its onset of decomposition. Furthermore, studies that have been performed on PETN thermal degradation often do not fully characterize or quantify the decomposition products. In this study, we heat PETN to melt temperatures and identify thermal decomposition products, morphology changes, and mass loss by ultrahigh-pressure liquid chromatography coupled to quadrupole time of flight mass spectrometry, scanning electron microscopy, nuclear magnetic resonance spectroscopy, and differential scanning calorimetry. For the first time, we quantify several decomposition products using independently prepared standards and establish the resulting melting point depression after the first melt. We also estimate the amount of decomposition relative to sublimation that we measure through gas evolution and evaluate the performance behavior of the molten material in commercial detonator configurations.

Tile-Based Pairwise Analysis of GC × GC-TOFMS Data to Facilitate Analyte Discovery and Mass Spectrum Purification

A new tile-based pairwise analysis workflow, termed 1v1 analysis, is presented to discover and identify analytes that differentiate two chromatograms collected using comprehensive two-dimensional (2D) gas chromatography coupled with time-of-flight mass spectrometry (GC × GC-TOFMS). Tile-based 1v1 analysis easily discovered all 18 non-native analytes spiked in diesel fuel within the top 30 hits, outperforming standard pairwise chromatographic analyses. However, eight spiked analytes could not be identified with multivariate curve resolution-alternating least-squares (MCR-ALS) nor parallel factor analysis (PARAFAC) due to background contamination. Analyte identification was achieved with class comparison enabled-mass spectrum purification (CCE-MSP), which obtains a pure analyte spectrum by normalizing the spectra to an interferent mass channel (m/z) identified from 1v1 analysis and subtracting the two spectra. This report also details the development of CCE-MSP assisted MCR-ALS, which removes the identified interferent m/z from the data prior to decomposition. In total, 17 out of 18 spiked analytes had a match value (MV) > 800 with both versions of CCE-MSP. For example, MCR-ALS and PARAFAC were unable to decompose the pure spectrum of methyl decanoate (MVs < 200) due to its low 2D chromatographic resolution (∼0.34) and high interferent-to-analyte signal ratio (∼30:1). By leveraging information gained from 1v1 analysis, CCE-MSP and CCE-MSP assisted MCR-ALS obtained a pure spectrum with an average MV of 908 and 964, respectively. Furthermore, tile-based 1v1 analysis was applied to track moisture damage in cacao beans, where 86 analytes with at least a 2-fold concentration change were discovered between the unmolded and molded samples. This 1v1 analysis workflow is beneficial for studies where multiple replicates are either unavailable or undesirable to save analysis time.

Class comparison enabled mass spectrum purification for comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry

Tile-based Fisher ratio (F-ratio) analysis is emerging as a versatile data analysis tool for supervised discovery-based experimentation using comprehensive two-dimensional (2D) gas chromatography coupled with time-of-flight mass spectrometry (GC × GC-TOFMS). None the less, analyte identification can often be marred by poor 2D resolution and low analyte abundance relative to overlapping compounds. Linear algebra-based chemometric methods, in particular multivariate curve resolution alternating least squares (MCR-ALS), parallel factor analysis (PARAFAC) and PARAFAC2, are often applied in an effort to address this situation. However, these chemometric methods can fail to produce an accurate spectrum when the analyte is at low 2D resolution and/or in low relative abundance. To address this challenge, we introduce class comparison enabled mass spectrum purification (CCE-MSP), a method that utilizes the underlying requirement for signal consistency of the background interference compounds between the two classes in the F-ratio analysis to purify the mass spectrum of the analyte hits. CCE-MSP is validated using a dataset obtained for a neat JP-8 jet fuel spiked with 14 sulfur containing compounds at two levels (15 ppm and 30 ppm), using the p-value and lack-of-fit (LOF) for each analyte hit as consistency metrics. A purified mass spectrum was produced for each spiked analyte hit and their mass spectrum match value (MV) was compared to the MV obtained by MCR-ALS, PARAFAC, and PARAFAC2. The resulting MV for CCE-MSP were found to be as good or better than these chemometric methods, eg., for 2-butyl-5-ethylthiophene with an analyte-to-interference relative signal abundance of 1:87 and a 2D resolution of 0.2, CCE-MSP produced a MV of 831, compared to 476 for MCR-ALS, 403 for PARAFAC, and 336 for PARAFAC2. CCE-MSP is also extended to obtain the purified spectrum for more than one analyte, eg., two analyte hits in overlapping hit locations. The spectra produced by CCE-MSP can also be utilized as estimates to facilitate quantitative signal decomposition using MCR-ALS.

Untargeted profiling and differentiation of geographical variants of wine samples using headspace solid-phase microextraction flow-modulated comprehensive two-dimensional gas chromatography with the support of tile-based Fisher ratio analysis

The volatile fraction of food, also called the food volatilome, is increasingly used to develop new fingerprinting approaches. The characterization of the food volatilome is important to achieve desired flavor profiles in food production processes, or to differentiate different products, with winemaking being one popular area of interest. In the present research, headspace solid-phase microextraction (HS SPME) coupled to flow-modulated comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry (FM GC×GC-TOFMS) was used to characterize geographical-based differences in the volatilome of five white "Grillo" wines (of Sicilian origin), comprising the five sample classes. All wines were produced with the same vinification method in 2019. To minimize the influence of minor bottle-to-bottle differences, three bottles of the same wine were randomly selected, and three samples were collected per bottle, resulting in nine sample replicates per wine. Particular emphasis was devoted to the operational conditions of a novel low duty cycle flow modulator. A fast FM GC×GC-TOFMS method with a modulation period of 700 ms and a re-injection period of 80 ms was developed. Following, the instrumental software was exploited to identify class-distinguishing analytes in the dataset via tile-based Fisher ratio analysis (i.e., ChromaTOF Tile). A tile size of 10 modulations (7 s) on the first dimension and 45 spectra (300 ms) on the second dimension was used to encompass average peak widths and to account for minor retention time shifting. Off-line software was used to apply an ANOVA test. A p-value of 0.01 was applied in order to select the most important class-distinguishing analytes, which were input to principal component analysis (PCA). The PCA scores plot showed distinct clustering of the wines according to geographical origin, although the loadings revealed that only a few analytes were necessary to differentiate the wines. However, a comprehensive flavor profile assessment underscored the importance of all the information output by the ChromaTOF Tile software.

Development of variance rank initiated-unsupervised sample indexing for gas chromatography-mass spectrometry analysis

Traditional non-targeted chemometric workflows for gas chromatography-mass spectrometry (GC-MS) data rely on using supervised methods, which requires a priori knowledge of sample class membership. Herein, we propose a simple, unsupervised chemometric workflow known as variance rank initiated-unsupervised sample indexing (VRI-USI). VRI-USI discovers analyte peaks exhibiting high relative variance across all samples, followed by k-means clustering on the individual peaks. Based upon how the samples cluster for a given peak, a sample index assignment is provided. Using a probabilistic argument, if the same sample index assignment appears for several discovered peaks, then this outcome strongly suggests that the samples are properly classified by that particular sample index assignment. Thus, relevant chemical differences between the samples have been discovered in an unsupervised fashion. The VRI-USI workflow is demonstrated on three, increasingly difficult datasets: simulations, yeast metabolomics, and human cancer metabolomics. For simulated GC-MS datasets, VRI-USI discovered 85-90% of analytes modeled to vary between sample classes. Nineteen out of 53 peaks in the peak table developed for the yeast metabolome dataset had the same sample index assignments, indicating that those indices are most likely due to class-distinguishing chemical differences. A t-test revealed that 22 out of 53 peaks were statistically significant (p < 0.05) when using those sample index assignments. Likewise, for the human cancer metabolomics study, VRI-USI discovered 25 analytes that were statistically different (p < 0.05) using the sample index assignments determined to highlight meaningful sample-based differences. For all datasets, the sample index assignments that were deduced from VRI-USI were the correct class-based difference when using prior knowledge. VRI-USI holds promise as an exploratory data analysis workflow for studies in which analysts do not readily have a priori class information or want to uncover the underlying nature of their dataset.

Unique ion filter-A data reduction tool for chemometric analysis of raw comprehensive two-dimensional gas chromatography-mass spectrometry data

Comprehensive gas chromatography with time of flight mass spectrometry is a powerful tool in the analysis of complex samples. Chemometric analysis of raw chromatographic data is more useful in one- and two-dimensional separations relative to peak tables. The data volume from such experiments generally necessitates the use of data reduction tools. Such tools often sacrifice some of the multivariate information in the mass to charge ratio dimension. The unique ion filter reduces the over-redundancy in two-dimensional gas chromatography-mass spectrometry data by limiting the data to a few unique/pseudo-unique ions, sub-peaks/slices in the first dimension, and spectra in the second dimension. We explore the performance of this algorithm through careful inspection of two-dimensional gas chromatography-mass spectrometry data before and after application of the filter. A reduction (99%) in the number of variables in a two-dimensional gas chromatography-mass spectrometry chromatogram passed on to subsequent analysis was observed. Feature selection times for model optimization reduced from 229 (±13) to 6.8 (±0.5) min when the filter was applied. An estimate of two unique/pseudo-unique ions, one sub-peak in the first dimension and five spectra in the second dimension were considered to provide a true representation of each chromatogram and provided enough information to achieve 100% model prediction accuracy.

Investigation of the limit of discovery using tile-based Fisher ratio analysis with comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry

Comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry (GC×GC-TOFMS) is followed by tile-based Fisher ratio (F-ratio) analysis to investigate the "limit of discovery" for low concentration levels of sulfur-containing compounds in JP8 jet fuel. A mixture of 14 sulfur-containing compounds was spiked at 30 ppm, 15 ppm, 3 ppm and 1.5 ppm into the neat fuel prior to GC×GC-TOFMS analysis with a "reversed" column format (aka polar first dimension (¹D) and non-polar second dimension (²D) column). Prior standard implementation of tile-based F-ratio analysis utilized an average F-ratio requiring a minimum of 3 mass channels (m/z) with the highest F-ratios. Herein, we explore the notion that use of the top F-ratio m/z for hitlist ranking is superior to the standard implementation for analytes near their limit-of-quantitation (LOQ), defined as an analyte concentration that produces a signal equal to ten times the standard deviation of the baseline noise (10σ_n). Hitlist ranking comparisons revealed that using only the top F-ratio m/z resulted in impressive improvements in discoverability for the low concentration comparisons. Specifically, for the 3 ppm versus neat hitlist, 1,4-oxathiane (LOQ = 2.5 ppm) improved from hit 114 via standard F-ratio analysis, to hit 25. For the 1.5 ppm versus neat hitlist, 2-propylthiophene (LOQ = 0.64 ppm) improved from hit 59 to 17, benzo[b]thiophene (LOQ = 1.1 ppm) from hit 98 to 28, and 2,5-dimethylthiophene (LOQ = 1.3 ppm) from hit 262 to 39. Additional hitlist ranking comparisons revealed the importance of proper tile size selection, as analyte discoverability deteriorated upon using either an inappropriately too small or too large of a tile.

Automating and Extending Comprehensive Two-Dimensional Gas Chromatography Data Processing by Interfacing Open-Source and Commercial Software

Comprehensive two-dimensional gas chromatography (GC×GC) is a powerful analytical tool for both nontargeted and targeted analyses. However, there is a need for more integrated workflows for processing and managing the resultant high-complexity datasets. End-to-end workflows for processing GC×GC data are challenging and often require multiple tools or software to process a single dataset. We describe a new approach, which uses an existing underutilized interface within commercial software to integrate free and open-source/external scripts and tools, tailoring the workflow to the needs of the individual researcher within a single software environment. To demonstrate the concept, the interface was successfully used to complete a first-pass alignment on a large-scale GC×GC metabolomics dataset. The analysis was performed by interfacing bespoke and published external algorithms within a commercial software environment to automatically correct the variation in retention times captured by a routine reference standard. Variation in ¹t_R and ²t_R was reduced on average from 8 and 16% CV prealignment to less than 1 and 2% post alignment, respectively. The interface enables automation and creation of new functions and increases the interconnectivity between chemometric tools, providing a window for integrating data-processing software with larger informatics-based data management platforms.

Statistical inference of mass channel purity from Fisher ratio analysis using comprehensive two-dimensional gas chromatography with time of flight mass spectrometry data

Tile-based Fisher ratio (F-ratio) analysis has recently been developed and validated for discovery-based studies of highly complex data collected using comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry (GC×GC-TOFMS). In previous studies, interpretation and utilization of F-ratio hit lists has relied upon manual decomposition and quantification performed by chemometric methods such as parallel factor analysis (PARAFAC), or via manual translation of the F-ratio hit list information to peak table quantitative information provided by the instrument software (ChromaTOF). Both of these quantification approaches are bottlenecks in the overall workflow. In order to address this issue, a more automatable approach to provide accurate relative quantification for F-ratio analyses was investigated, based upon the mass spectral selectivity provided via the F-ratio spectral output. Diesel fuel spiked with 15 analytes at four concentration levels (80, 40, 20, and 10 ppm) produced three sets of two class comparisons that were submitted to tile-based F-ratio analysis to obtain three hit lists, with an F-ratio spectrum for each hit. A novel algorithm which calculates the signal ratio (S-ratio) between two classes (eg., 80 ppm versus 40 ppm) was applied to all mass channels (m/z) in the F-ratio spectrum for each hit. A lack of fit (LOF) metric was utilized as a measure of peak purity and combined with F-ratio and p-values to study the relationship of each of these metrics with m/z purity. Application of a LOF threshold coupled with a p-value threshold yielded a subset of the most pure m/z for each of the 15 spiked analytes, evident by the low deviations (< 5%) in S-ratio relative to the true concentration ratio. A key outcome of this study was to demonstrate the isolation of pure m/z without the need for higher level signal decomposition algorithms.

Sedentariness and Urinary Metabolite Profile in Type 2 Diabetic Patients, a Cross-Sectional Study

Recent findings indicate a significant association between sedentary (SED)-time and type 2 diabetes mellitus (T2DM). The aim of this study was to investigate whether different levels of SED-time could impact on biochemical and physiological processes occurring in sedentary and physically inactive T2DM patients. In particular, patients from the “Italian Diabetes and Exercise Study (IDES)_2 trial belonging to the first and fourth quartile of SED-time were compared. Urine samples were analyzed by comprehensive two-dimensional gas chromatography (GC × GC) with parallel detection by mass spectrometry and flame ionization detection (GC × 2GC-MS/FID). This platform enables accurate profiling and fingerprinting of urinary metabolites while maximizing the overall information capacity, quantitation reliability, and response linearity. Moreover, using advanced pattern recognition, the fingerprinting process was extended to untargeted and targeted features, revealing diagnostic urinary fingerprints between groups. Quantitative metabolomics was then applied to analytes of relevance for robust comparisons. Increased levels of glycine, L-valine, L-threonine, L-phenylalanine, L-leucine, L-alanine, succinic acid, 2-ketoglutaric acid, xylitol, and ribitol were revealed in samples from less sedentary women. In conclusion, SED-time is associated with changes in urine metabolome signatures. These preliminary results suggest that reducing SED-time could be a strategy to improve the health status of a large proportion of diabetic patients.

Non-targeted discovery-based analysis for gas chromatography with mass spectrometry: A comparison of peak table, tile, and pixel-based Fisher ratio analysis

The ability to discover minute differences between samples or sample classes for gas chromatography coupled to mass spectrometry (GC-MS) can be a challenging endeavor, especially when those differences are not a priori. Fisher ratio (F-ratio) analysis is an apt technique to probe the differences between GC-MS chromatograms. F-ratio analysis is a supervised, non-targeted, discovery-based method that compares two different samples (or sample classes) to reduce the GC-MS dataset into a hit list composed of class distinguishing compounds. Three different F-ratio techniques, peak table, tile, and pixel-based were used to "discover" nine non-native analytes that were spiked into gasoline at four different nominal concentrations of 250, 85, 25, 5 parts-per-million (ppm). For the tile and pixel-based F-ratio calculations, a novel methodology is introduced to improve the sensitivity of the F-ratio calculations while reducing false positives. Furthermore, we use a combinatorial technique using null class comparisons, termed null distribution analysis, to determine a statistical F-ratio cutoff for analysis of the hit lists. The pixel-based algorithm was the most sensitive method and was able to "discover" all nine spiked analytes at a nominal concentration of 250 ppm albeit with one false positive interspersed towards the bottom of the hit list. The pixel-based software was also able to "discover" more of the spiked analytes at the lower concentrations with seven of the spiked analytes "discovered" at 85 ppm, four of the spiked analytes "discovered" at 25 ppm, and one analyte "discovered" at 5 ppm.

Automating data analysis for two-dimensional gas chromatography/time-of-flight mass spectrometry non-targeted analysis of comparative samples

Non-targeted analysis of environmental samples, using comprehensive two-dimensional gas chromatography coupled with time-of-flight mass spectrometry (GC × GC/ToF-MS), poses significant data analysis challenges due to the large number of possible analytes. Non-targeted data analysis of complex mixtures is prone to human bias and is laborious, particularly for comparative environmental samples such as contaminated soil pre- and post-bioremediation. To address this research bottleneck, we developed OCTpy, a Python™ script that acts as a data reduction filter to automate GC × GC/ToF-MS data analysis from LECO® ChromaTOF® software and facilitates selection of analytes of interest based on peak area comparison between comparative samples. We used data from polycyclic aromatic hydrocarbon (PAH) contaminated soil, pre- and post-bioremediation, to assess the effectiveness of OCTpy in facilitating the selection of analytes that have formed or degraded following treatment. Using datasets from the soil extracts pre- and post-bioremediation, OCTpy selected, on average, 18% of the initial suggested analytes generated by the LECO® ChromaTOF® software Statistical Compare feature. Based on this list, 63-100% of the candidate analytes identified by a highly trained individual were also selected by OCTpy. This process was accomplished in several minutes per sample, whereas manual data analysis took several hours per sample. OCTpy automates the analysis of complex mixtures of comparative samples, reduces the potential for human error during heavy data handling and decreases data analysis time by at least tenfold.