A Review of Chemometrics Applied to Comprehensive Two-dimensional Separations from 2008–2010

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.

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.

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.

Chemometric methods in data processing of mass spectrometry-based metabolomics: A review

This review focuses on recent and potential advances in chemometric methods in relation to data processing in metabolomics, especially for data generated from mass spectrometric techniques. Metabolomics is gradually being regarded a valuable and promising biotechnology rather than an ambitious advancement. Herein, we outline significant developments in metabolomics, especially in the combination with modern chemical analysis techniques, and dedicated statistical, and chemometric data analytical strategies. Advanced skills in the preprocessing of raw data, identification of metabolites, variable selection, and modeling are illustrated. We believe that insights from these developments will help narrow the gap between the original dataset and current biological knowledge. We also discuss the limitations and perspectives of extracting information from high-throughput datasets.

WITHDRAWN: Recent advances in chemometric methods for plant metabolomics: A review

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Oxygenated weathering products of Deepwater Horizon oil come from surprising precursors

Following the release of crude oil from the Macondo well in 2010, a wide range of weathering processes acted on the spilled oil. A recent study revealed that samples from this spill were oxidized into oxygenated hydrocarbons (OxHC) comprising more than 50% of the extracted hydrocarbons. The precursors of these compounds were not identified despite using a wide range of analytical tools, including gas chromatography (GC). To search for these precursors, over 40 samples were analyzed by comprehensive two-dimensional gas chromatography (GC×GC), one of the largest studies of its kind to date. Partial least squares regression was employed to elucidate the GC×GC peaks that could be the precursors of OxHC in our samples. We found that the formation of OxHC correlated with the disappearance of saturated hydrocarbons, including alkylcyclopentanes, alkyl cyclohexanes, alkylated bicyclic saturated compounds, tricyclic terpanpoids, and alkylbenzenes. These results indicate a previously under-reported chemodynamic process in oil spill weathering.