Simple automatic strategy for background drift correction in chromatographic data analysis

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.

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.

AntDAS-DDA: A New Platform for Data-Dependent Acquisition Mode-Based Untargeted Metabolomic Profiling Analysis with Advantage of Recognizing Insource Fragment Ions to Improve Compound Identification

Data-dependent acquisition (DDA) mode in ultra-high-performance liquid chromatography-high-resolution mass spectrometry (UHPLC-HRMS) can provide massive amounts of MS¹ and MS/MS information of compounds in untargeted metabolomics and can thus facilitate compound identification greatly. In this work, we developed a new platform called AntDAS-DDA for the automatic processing of UHPLC-HRMS data sets acquired under the DDA mode. Several algorithms, including extracted ion chromatogram extraction, feature extraction, MS/MS spectrum construction, fragment ion identification, and MS¹ spectrum construction, were developed within the platform. The performance of AntDAS-DDA was investigated comprehensively with a mixture of standard and complex plant data sets. Results suggested that features in complex sample matrices can be extracted effectively, and the constructed MS¹ and MS/MS spectra can benefit in compound identification greatly. The efficiency of compound identification can be improved by about 20%. AntDAS-DDA can take full advantage of MS/MS information in multiple sample analyses and provide more MS/MS spectra than single sample analysis. A comparison with advanced data analysis tools indicated that AntDAS-DDA may be used as an alternative for routine UHPLC-HRMS-based untargeted metabolomics. AntDAS-DDA is freely available at http://www.pmdb.org.cn/antdasdda.

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.

A general-purpose signal processing algorithm for biological profiles using only first-order derivative information

BACKGROUND

Automatic signal-feature extraction algorithms are crucial for profile processing in bioinformatics. Both baseline drift and noise seriously affect the position and peak area of signals. An efficient algorithm named the derivative passing accumulation (DPA) method for simultaneous baseline correction and signal extraction is presented in this article. It is an efficient method using only the first-order derivatives which are obtained through taking the simple differences.

RESULTS

We developed a new signal feature extracting procedure. The vector representing the discrete first-order derivative was divided into negative and positive parts and then accumulated to build a signal descriptor. The signals and background fluctuations are easily separated according to this descriptor via thresholding. In addition, the signal peaks are simultaneously located by checking the corresponding intervals in the descriptor. Therefore, the eternal issues of parsing the 1-dimensional output of detectors in biological instruments are solved together. Thereby, the baseline is corrected, and the signal peaks are extracted.

CONCLUSIONS

We have introduced a new method for signal peak picking, where baseline computation and peak identification are performed jointly. The testing results of both authentic and artificially synthesized data illustrate that the new method is powerful, and it could be a better choice for practical processing.

Metabolomics for a Millenniums-Old Crop: Tea Plant ( Camellia sinensis)

Tea cultivation and utilization dates back to antiquity. Today it is the most widely consumed beverage on earth due to its pleasant taste and several beneficial health properties attributed to specific metabolites. Metabolomics has a tremendous potential to correlate tea metabolites with taste and health properties in humans. Our review on the current application of metabolomics in the science of tea suggests that metabolomics is a promising frontier in the evaluation of tea quality, identification of functional genes responsible for key metabolites, investigation of their metabolic regulation, and pathway analysis in the tea plant. Furthermore, the challenges, possible solutions, and the prospects of metabolomics in tea science are reviewed.

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.

Salient space detection algorithm for signal extraction from contaminated and distorted spectrum

The proposed algorithm minimizes the influence of baseline distortion and exhibits good anti-noise ability and high real-time performance.

Evaluation of intensity drift correction strategies using MetaboDrift, a normalization tool for multi-batch metabolomics data

In recent years, mass spectrometry-based metabolomics has increasingly been applied to large-scale epidemiological studies of human subjects. However, the successful use of metabolomics in this context is subject to the challenge of detecting biologically significant effects despite substantial intensity drift that often occurs when data are acquired over a long period or in multiple batches. Numerous computational strategies and software tools have been developed to aid in correcting for intensity drift in metabolomics data, but most of these techniques are implemented using command-line driven software and custom scripts which are not accessible to all end users of metabolomics data. Further, it has not yet become routine practice to assess the quantitative accuracy of drift correction against techniques which enable true absolute quantitation such as isotope dilution mass spectrometry. We developed an Excel-based tool, MetaboDrift, to visually evaluate and correct for intensity drift in a multi-batch liquid chromatography - mass spectrometry (LC-MS) metabolomics dataset. The tool enables drift correction based on either quality control (QC) samples analyzed throughout the batches or using QC-sample independent methods. We applied MetaboDrift to an original set of clinical metabolomics data from a mixed-meal tolerance test (MMTT). The performance of the method was evaluated for multiple classes of metabolites by comparison with normalization using isotope-labeled internal standards. QC sample-based intensity drift correction significantly improved correlation with IS-normalized data, and resulted in detection of additional metabolites with significant physiological response to the MMTT. The relative merits of different QC-sample curve fitting strategies are discussed in the context of batch size and drift pattern complexity. Our drift correction tool offers a practical, simplified approach to drift correction and batch combination in large metabolomics studies.

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.