A simple multi-scale Gaussian smoothing-based strategy for automatic chromatographic peak extraction

Generating realistic data through modeling and parametric probability for the numerical evaluation of data processing algorithms in two-dimensional chromatography

BACKGROUND

Comprehensive two-dimensional chromatography generates complex data sets, and numerous baseline correction and noise removal algorithms have been proposed in the past decade to address this challenge. However, evaluating their performance objectively is currently not possible due to a lack of objective data.

RESULT

To tackle this issue, we introduce a versatile platform that models and reconstructs single-trace two-dimensional chromatography data, preserving peak parameters. This approach balances real experimental data with synthetic data for precise comparisons. We achieve this by employing a Skewed Lorentz-Normal model to represent each peak and creating probability distributions for relevant parameter sampling. The model's performance has been showcased through its application to two-dimensional gas chromatography data where it has created a data set with 458 peaks with an RMSE of 0.0048 or lower and minimal residuals compared to the original data. Additionally, the same process has been shown in liquid chromatography data.

SIGNIFICANCE

Data analysis is an integral component of any analytical method. The development of new data processing strategies is of paramount importance to tackle the complex signals generated by state-of-the-art separation technology. Through the use of probability distributions, quantitative assessment of algorithm performance of new algorithms is now possible. Therefore, creating new opportunities for faster, more accurate, and simpler data analysis development.

AntDAS-GCMS: A New Comprehensive Data Analysis Platform for GC-MS-Based Untargeted Metabolomics with the Advantage of Addressing the Time Shift Problem

Over the years, a number of state-of-the-art data analysis tools have been developed to provide a comprehensive analysis of data collected from gas chromatography-mass spectrometry (GC-MS). Unfortunately, the time shift problem remains unsolved in these tools. Here, we developed a novel comprehensive data analysis strategy for GC-MS-based untargeted metabolomics (AntDAS-GCMS) to perform total ion chromatogram peak detection, peak resolution, time shift correction, component registration, statistical analysis, and compound identification. Time shift correction was specifically optimized in this work. The information on mass spectra and elution profiles of compounds was used to search for inherent landmarks within analyzed samples to resolve the time shift problem across samples efficiently and accurately. The performance of our AntDAS-GCMS was comprehensively investigated by using four complex GC-MS data sets with various types of time shift problems. Meanwhile, AntDAS-GCMS was compared with advanced GC-MS data analysis tools and classic time shift correction methods. Results indicated that AntDAS-GCMS could achieve the best performance compared to the other methods.

Automated method development in high-pressure liquid chromatography

Method development in liquid chromatography is a crucial step in the optimization of analytical separations for various applications. However, it is often a challenging endeavour due to its time-consuming, resource intensive and costly nature, which is further hampered by its complexity requiring highly skilled and experienced scientists. This review presents an examination of the methods that are required for a completely automated method development procedure in liquid chromatography, aimed at taking the human out of the decision loop. Some of the presented approaches have recently witnessed an important increase in interest as they offer the promise to facilitate, streamline and speed up the method development process. The review first discusses the mathematical description of the separation problem by means of multi-criteria optimization functions. Two different strategies to resolve this optimization are then presented; an experimental and a model-based approach. Additionally, methods for automated peak detection and peak tracking are reviewed, which, upon integration in an instrument, allow for a completely closed-loop method development process. For each of these approaches, various currently applied methods are presented, recent trends and approaches discussed, short-comings pointed out, and future prospects highlighted.

Automated tuning of denoising algorithms for noise removal in chromatograms

Different algorithms, such as the Savitzky-Golay filter and Whittaker smoother, have been proposed to improve the quality of experimental chromatograms. These approaches avoid excessive noise from hampering data analysis and as such allow an accurate detection and quantification of analytes. These algorithms require fine-tuning of their hyperparameters to regulate their roughness and flexibility. Traditionally, this fine-tuning is done manually until a signal is obtained that removes the noise while conserving valuable peak information. More objective and automated approaches are available, but these are usually method specific and/or require previous knowledge. In this work, the L-and V-curve, k-fold cross-validation, autocorrelation function and residual variance estimation approach are introduced as alternative automated and generally applicable parameter tuning methods. These methods do not require any previous information and are compatible with a multitude of denoising methods. Additionally, for k-fold cross-validation, autocorrelation function and residual variance estimation, a novel implementation based on median estimators is proposed to handle the specific shape of chromatograms, typically composed of alternating flat baselines and sharp peaks. These tuning methods are investigated in combination with four denoising methods; the Savitsky-Golay filter, Whittaker smoother, sparsity assisted signal smoother and baseline estimation and denoising using sparsity approach. It is demonstrated that the median estimators approach significantly improves the denoising and information conservation performance of relevant smoother-tuner combinations up to a factor 4 for simulated datasets and even up to a factor 10 for an experimental chromatogram. Moreover, the parameter tuning methods relying on residual variance estimation, k-fold cross-validation and autocorrelation function lead to similar small root-mean squared errors on the different simulated datasets and experimental chromatograms. The sparsity assisted signal smoother and baseline estimation and denoising using sparsity approach, which both rely on the use of sparsity, systematically outperform the two other methods and are hence most appropriate for chromatograms.

Using UHPLC-HRMS-based comprehensive strategy to efficiently and accurately screen and identify illegal additives in health-care foods

Accurately and high-thoroughly screening illegal additives in health-care foods continues to be a challenging task in routine analysis for the ultrahigh performance liquid chromatography-high resolution mass spectrometry based techniques. In this work, we proposed a new strategy to identify additives in complex food matrices, which consists of both experimental design and advanced chemometric data analysis. At first, reliable features in the analyzed samples were screened based on a simple but efficient sample weighting design, and those related to illegal additives were screened with robust statistical analysis. After the MS¹ in-source fragment ion identification, both MS¹ and MS/MS spectra were constructed for each underlying compound, based on which illegal additives can be precisely identified. The performance of the developed strategy was demonstrated by using mixture and synthetic sample datasets, indicating an improvement of data analysis efficiency up to 70.3 %. Finally, the developed strategy was applied for the screening of unknown additives in 21 batches of commercially available health-care foods. Results indicated that at least 80 % of false-positive results can be reduced and 4 additives were screened and confirmed.

A study of flavor variations during the flaxseed roasting procedure by developed real-time SPME GC-MS coupled with chemometrics

Volatile compound variations during the roasting procedure play an essential role in the flaxseed-related product. In this work, we proposed a new strategy to high-throughput characterize the dynamic variations of flavors in flaxseed. Volatile compounds released at various roasting times were comprehensively investigated by a newly developed real-time solid-phase microextraction coupled with gas chromatography-mass spectrometry (GC-MS). Raw data files were analyzed by our advanced GC-MS data analysis software AntDAS-GCMS. Chemometric methods such as principal component analysis and partial least squares-discrimination analysis have realized the differences of samples with various roasting times. Finally, a total of 51 compounds from 11 aromas were accurately identified and confirmed with standards, and their variations as a function of roasting time were studied. In conclusion, we provided a new solution for the online monitoring of volatile compounds during the industrial roasting process.

Deep convolutional autoencoder for the simultaneous removal of baseline noise and baseline drift in chromatograms

Enhancement of chromatograms, such as the reduction of baseline noise and baseline drift, is often essential to accurately detect and quantify analytes in a mixture. Current methods have been well studied and adopted for decades and have assisted researchers in obtaining reliable results. However, these methods rely on relatively simple statistics of the data (chromatograms) which in some cases result in significant information loss and inaccuracies. In this study, a deep one-dimensional convolutional autoencoder was developed that simultaneously removes baseline noise and baseline drift with minimal information loss, for a large number and great variety of chromatograms. To enable the autoencoder to denoise a chromatogram to be almost, or completely, noise-free, it was trained on data obtained from an implemented chromatogram simulator that generated 190.000 representative simulated chromatograms. The trained autoencoder was then tested and compared to some of the most widely used and well-established denoising methods on testing datasets of tens of thousands of simulated chromatograms; and then further tested and verified on real chromatograms. The results show that the developed autoencoder can successfully remove baseline noise and baseline drift simultaneously with minimal information loss; outperforming methods like Savitzky-Golay smoothing, Gaussian smoothing and wavelet smoothing for baseline noise reduction (root mean squared error of 1.094 mAU compared to 2.074 mAU, 2.394 mAU and 2.199 mAU) and Savitkzy-Golay smoothing combined with asymmetric least-squares or polynomial fitting for baseline noise and baseline drift reduction (root mean absolute error of 1.171 mAU compared to 3.397 mAU and 4.923 mAU). Evidence is presented that autoencoders can be utilized to enhance and correct chromatograms and consequently improve and alleviate downstream data analysis, with the drawback of needing a carefully implemented simulator, that generates realistic chromatograms, to train the autoencoder.

Rapid detection of five pesticide residues using complexes of gold nanoparticle and porphyrin combined with ultraviolet visible spectrum

BACKGROUD

Pesticides are widely used to control insect infestation and weeds in agriculture. However, concerns about the pesticide residues in agricultural products have been raised in recent years because of public interest in health and food quality and safety. Thus, rapid, convenient, and accurate analytical methods for the detection and quantification of pesticides are urgently required.

RESULTS

A nanohybrid system composed of gold nanoparticles (AuNPs) and tetrakis(N-methyl-4-pyridiniumyl) porphyrin (TMPyP) was used as an optical probe for the detection and quantification of five pesticides (Paraquat, Dipterex, Dursban, methyl thiophanate and Cartap). The method is based on the aggregation effect of pesticides on the carboxyl group modified by AuNPs. Subsequently, with the help of particle swarm optimization-optimized sample weighted least squares-support vector machine (PSO-OSWLS-SVM), all the pesticides could be successfully quantified. In addition, partial least squares discriminant analysis (PLS-DA) was applied and the five pesticides were satisfactorily recognized based on data array obtained from the ultraviolet visible (UV-visible) spectra of AuNP-TMPyP complex. Furthermore, the quantitative and qualitative analysis of the five pesticides could be also achieved in the complex real samples, in which all the relative standard deviations (RSDs) were less than 0.3‰ and all the linear absolute correlation coefficients were more than 0.9990. Furthermore, recognition rate of the training set and the prediction set based on multiplicative scatter correction (MSC), or second-order derivative (2nd derivative) UV-visible spectra in PLS-DA model could reach 100%.

CONCLUSION

This method was successfully applied for the rapid and accurate determination of multicomponent pesticide residues in real food samples. © 2020 Society of Chemical Industry.

Food Phenotyping: Recording and Processing of Non-Targeted Liquid Chromatography Mass Spectrometry Data for Verifying Food Authenticity

Experiments based on metabolomics represent powerful approaches to the experimental verification of the integrity of food. In particular, high-resolution non-targeted analyses, which are carried out by means of liquid chromatography-mass spectrometry systems (LC-MS), offer a variety of options. However, an enormous amount of data is recorded, which must be processed in a correspondingly complex manner. The evaluation of LC-MS based non-targeted data is not entirely trivial and a wide variety of strategies have been developed that can be used in this regard. In this paper, an overview of the mandatory steps regarding data acquisition is given first, followed by a presentation of the required preprocessing steps for data evaluation. Then some multivariate analysis methods are discussed, which have proven to be particularly suitable in this context in recent years. The publication closes with information on the identification of marker compounds.

Differentiating Westlake Longjing tea from the first- and second-grade producing regions using ultra high performance liquid chromatography with quadrupole time-of-flight mass spectrometry-based untargeted metabolomics in combination with chemometrics

There are numerous articles published for geographical discrimination of tea. However, few research works focused on the authentication and traceability of Westlake Longjing green tea from the first- and second-grade producing regions because the tea trees are planted in a limited growing zone with identical cultivate condition. In this work, a comprehensive analytical strategy was proposed by ultrahigh performance liquid chromatography-quadrupole time-of-flight mass spectrometry-based untargeted metabolomics coupled with chemometrics. The automatic untargeted data analysis strategy was introduced to screen metabolites that expressed significantly among different regions. Chromatographic features of metabolites can be automatically and efficiently extracted and registered. Meanwhile, those that were valuable for geographical origin discrimination were screened based on statistical analysis and contents in samples. Metabolite identification was performed based on high-resolution mass values and tandem mass spectra of screened peaks. Twenty metabolites were identified, based on which the two-way encoding partial least squares discrimination analysis was built for geographical origin prediction. Monte Caro simulation results indicated that prediction accuracy was up to 99%. Our strategy can be applicable for practical applications in the quality control of Westlake Longjing green tea.

Recent applications of chemometrics in one- and two-dimensional chromatography

The proliferation of increasingly more sophisticated analytical separation systems, often incorporating increasingly more powerful detection techniques, such as high-resolution mass spectrometry, causes an urgent need for highly efficient data-analysis and optimization strategies. This is especially true for comprehensive two-dimensional chromatography applied to the separation of very complex samples. In this contribution, the requirement for chemometric tools is explained and the latest developments in approaches for (pre-)processing and analyzing data arising from one- and two-dimensional chromatography systems are reviewed. The final part of this review focuses on the application of chemometrics for method development and optimization.

An ultra-robust fingerprinting method for quality assessment of traditional Chinese medicine using multiple reaction monitoring mass spectrometry

Chromatographic fingerprinting has been perceived as an essential tool for assessing quality and chemical equivalence of traditional Chinese medicine. However, this pattern-oriented approach still has some weak points in terms of chemical coverage and robustness. In this work, we proposed a multiple reaction monitoring (MRM)-based fingerprinting method in which approximately 100 constituents were simultaneously detected for quality assessment. The derivative MRM approach was employed to rapidly design MRM transitions independent of chemical standards, based on which the large-scale fingerprinting method was efficiently established. This approach was exemplified on QiShenYiQi Pill (QSYQ), a traditional Chinese medicine-derived drug product, and its robustness was systematically evaluated by four indices: clustering analysis by principal component analysis, similarity analysis by the congruence coefficient, the number of separated peaks, and the peak area proportion of separated peaks. Compared with conventional ultraviolet-based fingerprints, the MRM fingerprints provided not only better discriminatory capacity for the tested normal/abnormal QSYQ samples, but also higher robustness under different chromatographic conditions (i.e., flow rate, apparent pH, column temperature, and column). The result also showed for such large-scale fingerprints including a large number of peaks, the angle cosine measure after min-max normalization was more suitable for setting a decision criterion than the unnormalized algorithm. This proof-of-concept application gives evidence that combining MRM technique with proper similarity analysis metrices can provide a highly sensitive, robust and comprehensive analytical approach for quality assessment of traditional Chinese medicine.

•

MRM fingerprints are proposed for quality assessment of traditional medicine.

•

MRM fingerprints show favorable robustness, coverage and discriminatory capacity.

•

Similarity analysis methods for such large-scale fingerprints are proposed.

Multi-scale Gaussian/Haar wavelet strategies coupled with sub-window factor analysis for an accurate alignment in nontargeted metabolic profiling to enhance herbal origin discrimination capability

Metabolic dataset can provide an overview of different herbal origin, which is conducted by some statistical procedures. Such results often deviate to a certain degree, due to peaks shifts in chromatographic signals. In order to solve this problem, an improved algorithm of combining sub-window factor analysis with the mass spectrum information is proposed. The algorithm uses a peak detection approach derived either from multi-scale Gaussian function or Haar wavelet to locate the peaks with different application scope; the candidate drift points at each peak are estimated by Fast Fourier transform cross correlation; Specifically, the best drift points at each candidate peaks are confirmed by sub-window factor analysis and mass spectrum information in nontargeted metabolic profiling. Finally, the peak regions were aligned against a reference chromatogram, and the non-peak regions were used linear interpolation. The chromatographic signals of 30 Bupleurum samples were aligned as an illustration of this algorithm, and they could be well distinguished using some statistical procedures. The result demonstrates that the presented method is stronger than other mass-spectra based algorithms, when facing the alignment of some co-eluted peaks.

Automatic data analysis workflow for ultra-high performance liquid chromatography-high resolution mass spectrometry-based metabolomics

Data analysis for ultra-performance liquid chromatography high-resolution mass spectrometry-based metabolomics is a challenging task. The present work provides an automatic data analysis workflow (AntDAS2) by developing three novel algorithms, as follows: (i) a density-based ion clustering algorithm is designed for extracted-ion chromatogram extraction from high-resolution mass spectrometry; (ii) a new maximal value-based peak detection method is proposed with the aid of automatic baseline correction and instrumental noise estimation; and (iii) the strategy that clusters high-resolution m/z peaks to simultaneously align multiple components by a modified dynamic programing is designed to efficiently correct time-shift problem across samples. Standard compounds and complex datasets are used to study the performance of AntDAS2. AntDAS2 is better than several state-of-the-art methods, namely, XCMS Online, Mzmine2, and MS-DIAL, to identify underlying components and improve pattern recognition capability. Meanwhile, AntDAS2 is more efficient than XCMS Online and Mzmine2. A MATLAB GUI of AntDAS2 is designed for convenient analysis and is available at the following webpage: http://software.tobaccodb.org/software/antdas2.

Chemometric Analysis of the Volatile Compounds Generated by Aspergillus carbonarius Strains Isolated from Grapes and Dried Vine Fruits

Ochratoxin A (OTA) contamination in grape production is an important problem worldwide. Microbial volatile organic compounds (MVOCs) have been demonstrated as useful tools to identify different toxigenic strains. In this study, Aspergillus carbonarius strains were classified into two groups, moderate toxigenic strains (MT) and high toxigenic strains (HT), according to OTA-forming ability. The MVOCs were analyzed by GC-MS and the data processing was based on untargeted profiling using XCMS Online software. Orthogonal projection to latent structures discriminant analysis (OPLS-DA) was performed using extract ion chromatogram GC-MS datasets. For contrast, quantitative analysis was also performed. Results demonstrated that the performance of the OPLS-DA model of untargeted profiling was better than the quantitative method. Potential markers were successfully discovered by variable importance on projection (VIP) and t-test. (E)-2-octen-1-ol, octanal, 1-octen-3-one, styrene, limonene, methyl-2-phenylacetate and 3 unknown compounds were selected as potential markers for the MT group. Cuparene, (Z)-thujopsene, methyl octanoate and 1 unknown compound were identified as potential markers for the HT groups. Finally, the selected markers were used to construct a supported vector machine classification (SVM-C) model to check classification ability. The models showed good performance with the accuracy of cross-validation and test prediction of 87.93% and 92.00%, respectively.

Joint Bounding of Peaks Across Samples Improves Differential Analysis in Mass Spectrometry-Based Metabolomics

As mass spectrometry-based metabolomics becomes more widely used in biomedical research, it is important to revisit existing data analysis paradigms. Existing data preprocessing efforts have largely focused on methods which start by extracting features separately from each sample, followed by a subsequent attempt to group features across samples to facilitate comparisons. We show that this preprocessing approach leads to unnecessary variability in peak quantifications that adversely impacts downstream analysis. We present a new method, bakedpi, for the preprocessing of both centroid and profile mode metabolomics data that relies on an intensity-weighted bivariate kernel density estimation on a pooling of all samples to detect peaks. This new method reduces this unnecessary quantification variability and increases power in downstream differential analysis.

Automatic time-shift alignment method for chromatographic data analysis

Time shift among samples remains a significant challenge in data analysis, such as quality control of natural plant extracts and metabolic profiling analysis, because this phenomenon may lead to invalid conclusions. In this work, we propose a new time shift alignment method, namely, automatic time-shift alignment (ATSA), for complicated chromatographic data analysis. This technique comprised the following alignment stages: (1) automatic baseline correction and peak detection stage for providing useful chromatographic information; (2) preliminary alignment stage through adaptive segment partition to correct alignment for the entire chromatogram; and (3) precise alignment stage based on test chromatographic peak information to accurately align time shift. In ATSA, the chromatographic peak information of both reference and test samples can be completely employed for time-shift alignment to determine segment boundaries and avoid loss of information. ATSA was used to analyze a complicated chromatographic dataset. The obtained correlation coefficients among samples and data analysis efficiency indicated that the influences of time shift can be considerably reduced by ATSA; thus accurate conclusion could be obtained.