Jamin C. Hoggard

Review of chemometric analysis techniques for comprehensive two dimensional separations data

Comprehensive two-dimensional (2D) separations, such as comprehensive 2D gas chromatography (GC×GC), liquid chromatography (LC×LC), and related instrumental techniques, provide very large and complex data sets. It is often up to the software to assist the analyst in transforming these complex data sets into useful information, and that is precisely where the field of chemometric data analysis plays a pivotal role. Chemometric tools for comprehensive 2D separations are continually being developed and applied as researchers make significant advances in novel state-of-the-art algorithms and software, and as the commercial sector continues to provide user friendly chemometric software. In this review, we build upon previous reviews of this topic, by focusing primarily on advances that have been reported in the past five years. Most of the reports focus on instrumental platforms using GC×GC with either flame ionization detection (FID) or time-of-flight mass spectrometry (TOFMS) detection, or LC×LC with diode array absorbance detection (DAD). The review covers the following general topics: data preprocessing techniques, target analyte techniques, comprehensive nontarget analysis techniques, and software for chemometrics in multidimensional separations.

Comprehensive analysis of yeast metabolite GC×GC–TOFMS data: combining discovery-mode and deconvolution chemometric software

The first extensive study of yeast metabolite GC x GC-TOFMS data from cells grown under fermenting, R, and respiring, DR, conditions is reported. In this study, recently developed chemometric software for use with three-dimensional instrumentation data was implemented, using a statistically-based Fisher ratio method. The Fisher ratio method is fully automated and will rapidly reduce the data to pinpoint two-dimensional chromatographic peaks differentiating sample types while utilizing all the mass channels. The effect of lowering the Fisher ratio threshold on peak identification was studied. At the lowest threshold (just above the noise level), 73 metabolite peaks were identified, nearly three-fold greater than the number of previously reported metabolite peaks identified (26). In addition to the 73 identified metabolites, 81 unknown metabolites were also located. A Parallel Factor Analysis graphical user interface (PARAFAC GUI) was applied to selected mass channels to obtain a concentration ratio, for each metabolite under the two growth conditions. Of the 73 known metabolites identified by the Fisher ratio method, 54 were statistically changing to the 95% confidence limit between the DR and R conditions according to the rigorous Students t-test. PARAFAC determined the concentration ratio and provided a fully-deconvoluted (i.e. mathematically resolved) mass spectrum for each of the metabolites. The combination of the Fisher ratio method with the PARAFAC GUI provides high-throughput software for discovery-based metabolomics research, and is novel for GC x GC-TOFMS data due to the use of the entire data set in the analysis (640 MB x 70 runs, double precision floating point).

A principal component analysis based method to discover chemical differences in comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry (GC × GC-TOFMS) separations of metabolites in plant samples

Comprehensive two-dimensional gas chromatography combined with time-of-flight mass spectrometry (GCxGC-TOFMS) provides high resolution separations of complex samples with a mass spectrum at every point in the separation space. The large volumes of multidimensional data obtained by GCxGC-TOFMS analysis are analyzed using a principal component analysis (PCA) method described herein to quickly and objectively discover differences between complex samples. In this work, we submitted 54 chromatograms to PCA to automatically compare the metabolite profiles of three different species of plants, namely basil (Ocimum basilicum), peppermint (Mentha piperita), and sweet herb stevia (Stevia rebaudiana), where there were 18 chromatograms for each type of plant. The 54 scores of the m/z 73 data set clustered in three groups according to the three types of plants. Principal component 1 (PC 1) separated the stevia cluster from the basil and peppermint clusters, capturing 61.84% of the total variance. Principal component 2 (PC 2) separated the basil cluster from the peppermint cluster, capturing 16.78% of the total variance. The PCA method revealed that relative abundances of amino acids, carboxylic acids, and carbohydrates were responsible for differentiating the three plants. A brief list of the 16 most significant metabolites is reported. After PCA, the 54 scores of the m/z 217 data set clustered in three groups according to the three types of plants, as well, yielding highly loaded variables corresponding with chemical differences between plants that were complementary to the m/z 73 information. The PCA data mining method is applicable to all of the monitored selective mass channels, utilizing all of the collected data, to discover unknown differences in complex sample profiles.

Development and Application of a Comprehensive Two-Dimensional Gas Chromatography with Time-of-Flight Mass Spectrometry Method for the Analysis of L-β-methylamino-alanine in Human Tissue

L-Beta-methylamino-alanine (BMAA) has been proposed as a worldwide contributor to neurodegenerative diseases, including Parkinson dementia complex (PDC) of Guam and Alzheimers disease (AD). Recent conflicting reports of the presence of this amino acid in human brain from patients affected by these diseases have made it necessary to develop methods that provide unambiguous detection in complex samples. Comprehensive two-dimensional gas chromatography coupled with time-of-flight-mass-spectrometry analysis (GCxGC-TOFMS) followed by a targeted Parallel Factor Analysis (PARAFAC) deconvolution method has been used recently in metabolomic investigations to separate, identify, and quantify components of complex biological specimens. We have extended and applied this methodology to the toxicological problem of detecting BMAA in extracts of brain tissue. Our results show that BMAA can be isolated from closely eluting compounds and detected in trace amounts in extracts of brain tissue spiked with low levels of this analyte, ranging from 2.5ppb to 50ppb, with a limit of detection (LOD) of 0.7ppb. This new method was sufficiently sensitive to detect BMAA in cerebral extracts of mice fed BMAA. This optimized approach was then applied to analyze tissue from humans; however, no BMAA was detected in the brain extracts from controls or patients with PDC or AD. Our results demonstrate the application of multidimensional chromatography-mass spectrometry methods and computational deconvolution analysis to the problem of detecting trace amounts of a potential toxin in brain extracts from mice and humans.

Impurity Profiling of a Chemical Weapon Precursor for Possible Forensic Signatures by Comprehensive Two-Dimensional Gas Chromatography/Mass Spectrometry and Chemometrics

In this report we present the feasibility of using analytical and chemometric methodologies to reveal and exploit the chemical impurity profiles from commercial dimethyl methylphosphonate (DMMP) samples to illustrate the type of forensic information that may be obtained from chemical-attack evidence. Using DMMP as a model compound of a toxicant that may be used in a chemical attack, we used comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry (GC x GC/TOF-MS) to detect and identify trace organic impurities in six samples of commercially acquired DMMP. The GC x GC/TOF-MS data was analyzed to produce impurity profiles for all six DMMP samples using 29 analyte impurities. The use of PARAFAC for the mathematical resolution of overlapped GC x GC peaks ensured clean spectra for the identification of many of the detected analytes by spectral library matching. The use of statistical pairwise comparison revealed that there were trace impurities that were quantitatively similar and different among five of the six DMMP samples. Two of the DMMP samples were revealed to have identical impurity profiles by this approach. The use of nonnegative matrix factorization indicated that there were five distinct DMMP sample types as illustrated by the clustering of the multiple DMMP analyses into five distinct clusters in the scores plots. The two indistinguishable DMMP samples were confirmed by their chemical supplier to be from the same bulk source. Sample information from the other chemical suppliers supported the idea that the other four DMMP samples were likely from different bulk sources. These results demonstrate that the matching of synthesized products from the same source is possible using impurity profiling. In addition, the identified impurities common to all six DMMP samples provide strong evidence that basic route information can be obtained from impurity profiles. Finally, impurities that may be unique to the sole bulk manufacturer of DMMP were found in some of the DMMP samples.

Handling within run retention time shifts in two-dimensional chromatography data using shift correction and modeling

The use of PARAFAC for modeling GC x GC-TOFMS peaks is well documented. This success is due to the trilinear structure of these data under ideal, or sufficiently close to ideal, chromatographic conditions. However, using temperature programming to cope with the general elution problem, deviations from trilinearity within a run are more likely to be seen for the following three cases: (1) compounds (i.e., analytes) severely broadened on the first column hence defined by many modulation periods, (2) analytes with a very high retention factor on the second column and likely wrapped around in that dimension, or (3) with fast temperature program rates. This deviation from trilinearity is seen as retention time-shifted peak profiles in subsequent modulation periods (first column fractions). In this report, a relaxed yet powerful version of PARAFAC, known as PARAFAC2 has been applied to handle this shift within the model step by allowing generation of individual peak profiles in subsequent first column fractions. An alternative approach was also studied, utilizing a standard retention time shift correction to restore the data trilinearity structure followed by PARAFAC. These two approaches are compared when identifying and quantifying a known analyte over a large concentration series where a certain shift is simulated in the successive first column fractions. Finally, the methods are applied to real chromatographic data showing severely shifted peak profiles. The pros and cons of the presented approaches are discussed in relation to the model parameters, the signal-to-noise ratio and the degree of shift.

Automated resolution of nontarget analyte signals in GC x GC-TOFMS data using parallel factor analysis.

We previously reported a method for the automated (objective) selection of a PARAFAC model having an appropriate number of factors for mathematical resolution of signal from a target analyte in GC x GC-TOFMS data (i.e., for an analysis in which the identity of the analyte is known a priori). While the previous target method has been successfully applied in several studies, the target method requires that the identity of the analyte be known. Also, multiple applications of the target method are required in cases where several analytes of interest are present in a single subsection of the chromatogram. Thus, having to know the analyte identity a priori restricts the applicability of the automated implementation of PARAFAC. The method presented in this report generalizes the previous method to allow analysis of one or more nontarget analyte signals in a subsection of a GC x GC-TOFMS chromatogram (i.e., for analyses when identities of analyte and interferences are not known a priori), thereby addressing and overcoming the limitations of the target method. Herein, we put the nontarget analyte PARAFAC method into theoretical context and illustrate the mechanics of the method using simulated data. We use real experimental GC x GC-TOFMS data to demonstrate the broad applicability of the method, with various analysis situations selected to illustrate challenging chemical analysis scenarios.

Time-Dependent Profiling of Metabolites from Snf1 Mutant and Wild Type Yeast Cells

The effect of sampling time in the context of growth conditions on a dynamic metabolic system was investigated in order to assess to what extent a single sampling time may be sufficient for general application, as well as to determine if useful kinetic information could be obtained. A wild type yeast strain (W) was compared to a snf1Delta mutant yeast strain (S) grown in high-glucose medium (R) and in low-glucose medium containing ethanol (DR). Under these growth conditions, different metabolic pathways for utilizing the different carbon sources are expected to be active. Thus, changes in metabolite levels relating to the carbon source in the growth medium were anticipated. Furthermore, the Snf1 protein kinase complex is required to adapt cellular metabolism from fermentative R conditions to oxidative DR conditions. So, differences in intracellular metabolite levels between the W and S yeast strains were also anticipated. Cell extracts were collected at four time points (0.5, 2, 4, 6 h) after shifting half of the cells from R to DR conditions, resulting in 16 sample classes (WR, WDR, SR, SDR) x (0.5, 2, 4, 6 h). The experimental design provided time course data, so temporal dependencies could be monitored in addition to carbon source and strain dependencies. Comprehensive two-dimensional (2D) gas chromatography coupled to time-of-flight mass spectrometry (GC x GC-TOFMS) was used with discovery-based data mining algorithms ( Anal. Chem. 2006, 78, 5068-5075 (ref 1); J. Chromatogr., A 2008, 1186, 401-411 (ref 2)) to locate regions within the 2D chromatograms (i.e., metabolites) that provided chemical selectivity between the 16 sample classes. These regions were mathematically resolved using parallel factor analysis to positively identify the metabolites and to acquire quantitative results. With these tools, 51 unique metabolites were identified and quantified. Various time course patterns emerged from these data, and principal component analysis (PCA) was utilized as a comparison tool to determine the sources of variance between these 51 metabolites. The effect of sampling time was investigated with separate PCA analyses using various subsets of the data. PCA utilizing all of the time course data, averaged time course data, and each individual time point data set independently were performed to discern the differences. For the yeast strains examined in the current study, data collection at either 4 or 6 h provided information comparable to averaged time course data, albeit with a few metabolites missing using a single sampling time point.

Achieving high peak capacity production for gas chromatography and comprehensive two-dimensional gas chromatography by minimizing off-column peak broadening

By taking into consideration band broadening theory and using those results to select experimental conditions, and also by reducing the injection pulse width, peak capacity production (i.e., peak capacity per separation time) is substantially improved for one dimensional (1D-GC) and comprehensive two dimensional (GC×GC) gas chromatography. A theoretical framework for determining the optimal linear gas velocity (the linear gas velocity producing the minimum H), from experimental parameters provides an in-depth understanding of the potential for GC separations in the absence of extra-column band broadening. The extra-column band broadening is referred to herein as off-column band broadening since it is additional band broadening not due to the on-column separation processes. The theory provides the basis to experimentally evaluate and improve temperature programmed 1D-GC separations, but in order to do so with a commercial 1D-GC instrument platform, off-column band broadening from injection and detection needed to be significantly reduced. Specifically for injection, a resistively heated transfer line is coupled to a high-speed diaphragm valve to provide a suitable injection pulse width (referred to herein as modified injection). Additionally, flame ionization detection (FID) was modified to provide a data collection rate of 5kHz. The use of long, relatively narrow open tubular capillary columns and a 40°C/min programming rate were explored for 1D-GC, specifically a 40m, 180μm i.d. capillary column operated at or above the optimal average linear gas velocity. Injection using standard auto-injection with a 1:400 split resulted in an average peak width of ∼1.5s, hence a peak capacity production of 40peaks/min. In contrast, use of modified injection produced ∼500ms peak widths for 1D-GC, i.e., a peak capacity production of 120peaks/min (a 3-fold improvement over standard auto-injection). Implementation of modified injection resulted in retention time, peak width, peak height, and peak area average RSD%s of 0.006, 0.8, 3.4, and 4.0%, respectively. Modified injection onto the first column of a GC×GC coupled with another high-speed valve injection onto the second column produced an instrument with high peak capacity production (500-800peaks/min), ∼5-fold to 8-fold higher than typically reported for GC×GC.

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 (2D) gas chromatography coupled with time-of-flight mass spectrometry (GC × GC-TOFMS) is a versatile instrumental platform capable of collecting highly informative, yet highly complex, chemical data for a variety of samples. Fisher-ratio (F-ratio) analysis applied to the supervised comparison of sample classes algorithmically reduces complex GC × GC-TOFMS data sets to find class distinguishing chemical features. F-ratio analysis, using a tile-based algorithm, significantly reduces the adverse effects of chromatographic misalignment and spurious covariance of the detected signal, enhancing the discovery of true positives while simultaneously reducing the likelihood of detecting false positives. Herein, we report a study using tile-based F-ratio analysis whereby four non-native analytes were spiked into diesel fuel at several concentrations ranging from 0 to 100 ppm. Spike level comparisons were performed in two regimes: comparing the spiked samples to the nonspiked fuel matrix and to each other at relative concentration factors of two. Redundant hits were algorithmically removed by refocusing the tiled results onto the original high resolution pixel level data. To objectively limit the tile-based F-ratio results to only features which are statistically likely to be true positives, we developed a combinatorial technique using null class comparisons, called null distribution analysis, by which we determined a statistically defensible F-ratio cutoff for the analysis of the hit list. After applying null distribution analysis, spiked analytes were reliably discovered at ∼1 to ∼10 ppm (∼5 to ∼50 pg using a 200:1 split), depending upon the degree of mass spectral selectivity and 2D chromatographic resolution, with minimal occurrence of false positives. To place the relevance of this work among other methods in this field, results are compared to those for pixel and peak table-based approaches.