Solving chromatographic challenges in comprehensive two-dimensional gas chromatography–time-of-flight mass spectrometry using multivariate curve resolution–alternating least squares

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.

Random swapping, an effective and efficient way to boost the intensities of cross peaks in a 2D asynchronous spectrum

Analysis of mixture via chromatographic-spectroscopic and analogous experiments is a common task in analytical chemistry. A 2D/nD asynchronous spectrum is effective in retrieving spectra of pure substances even if different components cannot be separated. However, noise in the 2D/nD asynchronous spectrum becomes a bottleneck in the analysis. Finding a suitable sequence of the 1D spectra used in constructing the 2D/nD asynchronous spectrum is helpful to improve the signal-to-noise level. A 2D/nD asynchronous spectrum is often produced via a large number of 1D spectra. The resultant colossal number of the possible sequences makes stochastic search the only possible way to find a suitable sequence. Random changing (RC) and random swapping (RS) are two ways to obtain a new sequence. We found that the possibility of finding a better sequence via an RS is significantly higher than that via an RC in the advanced stage of stochastic searching. This is the reason why the performance of RS is superior to that of RC in two model systems where 2D asynchronous spectra are used. We applied the RS approach on the analysis of water/isopropanol mixtures, and satisfactory sequences are acquired with affordable computational cost. Thus, the RS approach brings about an opportunity increase the signal-to-noise level of a 2D asynchronous spectrum in the analysis of the bilinear data from complex mixed samples.

Antiretroviral Drugs in African Surface Waters: Prevalence, Analysis, and Potential Remediation

The sources, ecotoxicological impact, and potential remediation strategies of antiretroviral drugs (ARVDs) as emerging contaminants in surface waters are reviewed based on recent literature. The occurrence of ARVDs in water bodies raises concern because many communities in Africa depend on rivers for water resources. Southern Africa is a potential hotspot regarding ARVD contamination due to relatively high therapeutic application and detection thereof in water bodies. Efavirenz and nevirapine are the most persistent in effluents and are prevalent in surface water based on environmental concentrations. Whereas the highest concentration of efavirenz reported in Kenya was 12.4 µg L^-1 , concentrations as high as 119 and 140 µg L^-1 have been reported in Zambia and South Africa, respectively. Concentrations of ARVDs ranging from 670 to 34 000 ng L^-1 (influents) and 540 to 34 000 ng L^-1 (effluents) were determined in wastewater treatment plants in South Africa, compared with Europe, where reported concentrations range from less than limit of detection (LOD) to 32 ng L^-1 (influents) and less than LOD to 22 ng L^-1 (effluents). The present African-based review suggests the need for comprehensive toxicological and risk assessment of these emerging pollutants in Africa, with the intent of averting environmental hazards and the development of sustainable remediation strategies. Environ Toxicol Chem 2022;41:247-262. © 2021 SETAC.

A novel systematic absence of cross peaks-based 2D-COS approach for bilinear data

A novel approach to use two-dimensional correlation spectroscopy (2D-COS) to analyze bilinear data is proposed. A phenomenon called Systematic Absence of Cross Peaks (SACPs) is observed in a 2D asynchronous spectrum. Two theorems relevant to SACPs have been derived. The SACP-based 2D-COS method has been successfully applied on analyzing bilinear data from mixed samples (including one model system and two real systems). Implicit isolated peaks can be identified and assigned to different components based on characteristic pattern of SACPs even if the time-related profiles of different components are severely overlapped. Based on the results of SACPs, spectra of pure components can be retrieved. Identification of SACPs can still be achieved in the presence of artifacts. Thus, neither noise nor baseline drift can produce significant influence on the results obtained from the approach described in this paper. We have used several well-established chemometric methods, including N-Findr, VCA, and MCR with various initial settings, on two systems that can be successfully solved using the 2D-COS method. The chemometric methods mentioned above cannot provide correct spectra of pure components because of severe problem of rotational ambiguity derived from severe overlapping of the time-related profiles. Only when the information from SACPs in 2D-COS is used as additional constraints in MCR calculation, correct spectra can be obtained. That is to say, the SACP-based 2D-COS method provides intrinsic information which is crucial in the analysis of chromatographic-spectroscopic and analogous data even if the time-related profiles of different components overlap severely.

A preliminary study on constructing a high-dimensional asynchronous spectrum to analyze bilinear data

A novel approach to constructing high-dimensional asynchronous spectra (nD-Asyn) is proposed. Three theorems relevant to 1D slices of nD-Asyn are revealed. nD-Asyn is used to analyze bilinear data from mixtures containing multiple components obtained via hyphenated techniques. The spectral contribution of different components can be removed in a stepwise manner by increasing the dimensions of asynchronous spectra. As a result, the spectra of different components can be faithfully recovered even if the time-related profiles of different components severely overlap. Moreover, correct results can still be obtained via the nD-Asyn even if a considerable level of noise and baseline drift are present. The nD-Asyn approach is compared with MCR-ALS using different constraints in analyzing the data for a simulated and also for a real system. The nD-Asyn produced correct spectrum of every component. Only when complete constraints obtained from nD-Asyn method is utilized in the MCR-ALS calculation, correct spectra of all the components can be obtained. Thus, nD-Asyn can be used alone or in conjunction with MCR-ALS to analyze bilinear data containing contributions of multiple components.

Linking the morphological and metabolomic response of Lactuca sativa L exposed to emerging contaminants using GC × GC-MS and chemometric tools

The occurrence of contaminants of emerging concern (CECs) in irrigation waters (up to low μg L-1) and irrigated crops (ng g-1 in dry weight) has been reported, but the linkage between plant morphological changes and plant metabolomic response has not yet been addressed. In this study, a non-targeted metabolomic analysis was performed on lettuce (Lactuca sativa L) exposed to 11 CECs (pharmaceuticals, personal care products, anticorrosive agents and surfactants) by irrigation. The plants were watered with different CEC concentrations (0-50 µg L-1) for 34 days under controlled conditions and then harvested, extracted, derivatised and analysed by comprehensive two-dimensional gas chromatography coupled to a time-of-flight mass spectrometer (GC × GC-TOFMS). The resulting raw data were analysed using multivariate curve resolution (MCR) and partial least squares (PLS) methods. The metabolic response indicates that exposure to CECs at environmentally relevant concentrations (0.05 µg L-1) can cause significant metabolic alterations in plants (carbohydrate metabolism, the citric acid cycle, pentose phosphate pathway and glutathione pathway) linked to changes in morphological parameters (leaf height, stem width) and chlorophyll content.

Chemometric analysis of comprehensive two dimensional gas chromatography-mass spectrometry metabolomics data

Comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry (GC×GC-TOFMS) is a powerful tool for separation and identification of analytes in complex natural samples. In this paper, different chemometric methods were compared for fast non-targeted GC×GC-TOFMS metabolomic profiling of the crustaceous species Daphnia magna and a general chemometric strategy and workflow is proposed. The strategy proposed in this work combined the compression of GC×GC-TOFMS data matrices in the retention time direction using wavelets and the appropriate column-wise data matrix augmentation arrangement, and its modeling by Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS). Using the proposed strategy, eighty different D. magna metabolites were resolved and identified. After calculation of the peak capacities of different columns and peak area changes of these metabolites, the best instrumental configuration and column combination for the GC×GC-TOFMS metabolomic study of D. magna are proposed and discussed. The procedure described in this work can be applied as a general method for untargeted GC×GC-TOFMS data processing and metabolomic profiling.

A road map for multi-way calibration models

A road map is proposed for the selection of a multi-way calibration model according to the data properties.

Effect of different drying methods on the myosin structure, amino acid composition, protein digestibility and volatile profile of squid fillets

The impacts of freeze drying (FD), hot-air drying (AD), and heat pump drying (HPD) on myosin structure, amino acid composition, protein digestibility and volatile compounds of squid (Todarodes pacificus) fillets were evaluated. Freeze-dried squids showed similar amino acid composition to that of raw squids, but differed from that of AD and HPD samples. The percentage of in vitro digestibility followed the order of FD (76.81%)>HPD (70.51%)>raw (67.99%)>AD (61.47%) samples. AD caused more damage to squid myosin structure than HPD, while FD effectively retained the myosin integrity. Drying decreased total number of volatile compounds, but increased the content of total volatile compounds based on GC × GC-TOFMS results. HPD and AD samples had the highest and lowest total numbers and contents of volatiles, respectively. In general, FD provided squids with the best quality, followed by HPD. Considering the production cost and product quality, HPD demonstrated the potential for industrial application.

Chemometrics-assisted effect-directed analysis of crude and refined oil using comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry

An effect-directed analysis (EDA) of fresh and artificially weathered (evaporated, photooxidized) samples of North Sea crude oil and residual heavy fuel oil is presented. Aliphatic, aromatic, and polar oil fractions were tested for the presence of aryl hydrocarbon receptor (AhR) agonist and androgen receptor (AR) antagonist, demonstrating for the first time the AR antagonist effects in the aromatic and, to a lesser extent, polar fractions. An extension of the typical EDA strategy to include an N-way partial least-squares (N-PLS) model capable of relating the comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry (GC × GC-TOFMS) data set to the bioassay data obtained from normal-phase LC fractions is proposed. The predicted AhR binding effects in the fresh and artificially weathered aromatic oil fractions facilitated the identification of alkyl-substituted three- and four-ring aromatic systems in the active fractions through the weighting of their contributions to the observed effects. A N-PLS chemometric model is demonstrated as a potentially useful strategy for future EDA studies that can streamline the compound identification process and provide additional reduction of samples' complexity. The AhR binding effects of the suspected compounds predicted by N-PLS and identified by GC × GC-TOFMS were confirmed using quantitative structure-activity relationship (QSAR) estimates.