Large-scale human metabolomics studies: a strategy for data (pre-) processing and validation

A large metabolomics study was performed on 600 plasma samples taken at four time points before and after a single intake of a high fat test meal by obese and lean subjects. All samples were analyzed by a liquid chromatography-mass spectrometry (LC-MS) lipidomic method for metabolic profiling. A pragmatic approach combining several well-established statistical methods was developed for processing this large data set in order to detect small differences in metabolic profiles in combination with a large biological variation. Such metabolomics studies require a careful analytical and statistical protocol. The strategy included data preprocessing, data analysis, and validation of statistical models. After several data preprocessing steps, partial least-squares discriminant analysis (PLS-DA) was used for finding biomarkers. To validate the found biomarkers statistically, the PLS-DA models were validated by means of a permutation test, biomarker models, and noninformative models. Univariate plots of potential biomarkers were used to obtain insight in up- or downregulation. The strategy proposed proved to be applicable for dealing with large-scale human metabolomics studies.     

Optimizing the use of quality control samples for signal drift correction in large-scale urine metabolic profiling studies

The evident importance of metabolic profiling for biomarker discovery and hypothesis generation has led to interest in incorporating this technique into large-scale studies, e.g., clinical and molecular phenotyping studies. Nevertheless, these lengthy studies mandate the use of analytical methods with proven reproducibility. An integrated experimental plan for LC-MS profiling of urine, involving sample sequence design and postacquisition correction routines, has been developed. This plan is based on the optimization of the frequency of analyzing identical quality control (QC) specimen injections and using the QC intensities of each metabolite feature to construct a correction trace for all the samples. The QC-based methods were tested against other current correction practices, such as total intensity normalization. The evaluation was based on the reproducibility obtained from technical replicates of 46 samples and showed the feature-based signal correction (FBSC) methods to be superior to other methods, resulting in ~1000 and 600 metabolite features with coefficient of variation (CV) < 15% within and between two blocks, respectively. Additionally, the required frequency of QC sample injection was investigated and the best signal correction results were achieved with at least one QC injection every 2 h of urine sample injections (n = 10). Higher rates of QC injections (1 QC/h) resulted in slightly better correction but at the expense of longer total analysis time.     

Solid-phase extraction and quantitative measurement of omega and omega-1 metabolites of JWH-018 and JWH-073 in human urine

The aminoalkylindole agonists JWH-018 and JWH-073 are contained in "K2/SPICE" products sold as "legal marijuana". Previous human metabolic studies have identified (ω)-hydroxyl and (ω)-carboxyl metabolites as biomarkers that are indicative of product use. However, other primary metabolites exhibiting similar chromatographic properties and mass spectra are also excreted in human urine. Analytical standards were used in this study to identify new primary metabolites as (ω-1)-hydroxyl derivatives of JWH-018 and JWH-073. The liquid chromatography tandem mass spectrometry (LC-MS/MS) procedure, coupled with an automated solid-phase extraction procedure incorporating deuterium-labeled internal standards, provides rapid resolution of the (ω)- and (ω-1) metabolites with adequate sensitivity, precision, and accuracy for trace analysis in human urine. Results from four urine specimens collected after individuals reportedly self-administered either JWH-018 or a mixture of JWH-018 and JWH-073 showed the following: (1) all tested metabolites were excreted in high concentrations, (2) (ω)- and (ω-1)-hydroxyl metabolites were exclusively excreted as glucuronic acid conjugates, and (3) ～5%-80% of the (ω)-carboxyl metabolites was excreted as glucuronic acid conjugates. This is the first report to identify and quantify (ω-1)-hydroxyl metabolites of JWH-018 and JWH-073 and the first to incorporate automated extraction procedures using deuterium-labeled internal standards. Full clinical validation awaits further testing.     

Controlled protein precipitation in combination with chip-based nanospray infusion mass spectrometry. An approach for metabolomics profiling of plasma

Liquid chromatography-mass spectrometry (LC-MS) is a common method for profiling biological samples in metabolomics. However, LC-MS data of metabolomic studies are often affected by high noise levels, retention time shifts, and high variability in signal intensities. With a new chip-based nanoelectrospray source it becomes possible to directly infuse complex biological samples such as plasma without any chromatographic separation beforehand. In combination with highly diluted samples and long data acquisition times, the parallel analysis of hundreds of compounds is now possible. In a proof-of-concept study, 10 human plasma samples from females and males were analyzed with the intention to separate the two groups by their different metabolomes. The reproducibility was so high that statistical analysis of the data could be performed without prior normalization. Two groups of female and male samples were separated by a supervised machine learning algorithm, principal component analysis, and hierarchical clustering. Peaks contributing to the group separation were characterized by accurate mass measurement and MS-MS fragmentation and by spiking experiments. The feasibility of direct sample infusion using the new chip-based nanoelectrospray device opens a new dimension for the rapid parallel analysis of complex biological mixtures.     

Enhancing the coverage of the urinary metabolome by sheathless capillary electrophoresis-mass spectrometry

Sheathless capillary electrophoresis-mass spectrometry (CE-MS), using a porous tip sprayer, is proposed for the first time for highly sensitive metabolic profiling of human urine. A representative metabolite mixture and human urine were used for evaluation of the sheathless CE-MS platform. For test compounds, relative standard deviations (RSDs) for migration times and peak areas were below 2% and 12%, respectively, and an injection volume of only ～8 nL resulted in detection limits between 11 and 120 nM. Approximately 900 molecular features were detected in human urine by sheathless CE-MS whereas about 300 molecular features were found with classical sheath-liquid CE-MS. This difference can probably be attributed to an improved ionization efficiency and increased sensitivity at low flow-rate conditions. The integration of transient-isotachophoresis (t-ITP) as an in-capillary preconcentration procedure in sheathless CE-MS further resulted in subnanomolar limits of detection for compounds of the metabolite mixture, and more than 1300 molecular features were observed in urine. Compared to the classical CE-MS approaches, the integration of t-ITP combined with the use of a sheathless interface provides up to 2 orders of magnitude sensitivity improvement. Hence, sheathless CE-MS can be used for in-depth metabolic profiling of biological samples, and we anticipate that this approach will yield unique information in the field of metabolomics.     

Impact of analytical bias in metabonomic studies of human blood serum and plasma

Concurrent with the explosion in the number of publications reporting biomarker discovery by profiling technologies, such as proteomics and pattern recognition, has been the increase in evidence highlighting the susceptibility of these approaches to analytical and experimental bias. The work presented here addresses these timely issues by delivering a detailed characterization of the effect of common sources of bias in clinical studies on serum and plasma profiles generated by a key technology in metabonomics, NMR spectroscopy. Specifically, differences in composition when blood samples were collected onto and in the absence of ice, over a series of serum-clot contact times, the stability of NMR-prepared samples over time and the effect on the metabolic profile of freeze-thawing were examined. While differences between individuals were far greater than variation from any other experimental factor, each of the conditions examined did cause slight alterations to the NMR profile that could produce a systematic bias. Variation due to clotting time caused changes in energy metabolites, which were delayed by ice with no other spectral effects. Room-temperature stability and hence NMR spectral repeatability were high (<1% intrasample variation). Higher molecular weight species such as lipoproteins were more susceptible to the variations present in the examined factors. These observations have implications for profiling study design, and hence, our results form a new and valuable resource for those attempting clinical metabolic profiling, for regulatory agencies involved in the licensing of clinical tests and in the generation of international reporting standards for metabonomics.     

In situ surface sampling of biological objects and preconcentration of their volatiles for chromatographic analysis

This report describes a rolling stir bar sampling procedure for volatile organic compounds (VOCs) present on various biological surfaces. In combination with thermal desorption/gas chromatography/mass spectrometry, this analytical technique was initially tested for quantitative profiling of human skin VOCs. It is also applicable to additional hydrophobic surfaces such as agricultural products, plant materials, and bird feathers. Use of embedded internal standards provides highly reproducible and quantitative results for a wide variety of sampled trace components. The samples of collected human skin VOCs and standards were found stable under cool storage conditions for at least 14 days, making this approach suitable for field biological and agricultural studies. Additionally, this methodology appears to have potential for forensic and toxicological investigations, as suggested through the analyses of VOC profiles of the human thumb prints recovered from a nonbiological smooth surface.     

Intra- and interlaboratory reproducibility of ultra performance liquid chromatography-time-of-flight mass spectrometry for urinary metabolic profiling

Liquid chromatography coupled to mass spectrometry (LC-MS) is a major platform in metabolic profiling but has not yet been comprehensively assessed as to its repeatability and reproducibility across multiple spectrometers and laboratories. Here we report results of a large interlaboratory reproducibility study of ultra performance (UP) LC-MS of human urine. A total of 14 stable isotope labeled standard compounds were spiked into a pooled human urine sample, which was subject to a 2- to 16-fold dilution series and run by UPLC coupled to time-of-flight MS at three different laboratories all using the same platform. In each lab, identical samples were run in two phases, separated by at least 1 week, to assess between-day reproducibility. Overall, platform reproducibility was good with median mass accuracies below 12 ppm, median retention time drifts of less than 0.73 s and coefficients of variation of intensity of less than 18% across laboratories and ionization modes. We found that the intensity response was highly linear within each run, with a median R(2) of 0.95 and 0.93 in positive and negative ionization modes. Between-day reproducibility was also high with a mean R(2) of 0.93 for a linear relationship between the intensities of ions recorded in the two phases across the laboratories and modes. Most importantly, between-lab reproducibility was excellent with median R(2) values of 0.96 and 0.98 for positive and negative ionization modes, respectively, across all pairs of laboratories. Interestingly, the three laboratories observed different amounts of adduct formation, but this did not appear to be related to reproducibility observed in each laboratory. These studies show that UPLC-MS is fit for the purpose of targeted urinary metabolite analysis but that care must be taken to optimize laboratory systems for quantitative detection due to variable adduct formation over many compound classes.     

Hydrophilic interaction chromatography for mass spectrometric metabonomic studies of urine

High-performance liquid chromatography (LC) coupled to mass spectrometry (MS) is increasingly being used for urinary metabonomic studies. Most studies utilize reversed-phase separation techniques, which are not suited to retaining highly polar analytes. Metabonomic studies should encompass a representative "fingerprint" that contains the largest amount of information possible. In this work, we have analyzed human urine samples with LC-MS, comparing traditional reversed-phase separation with hydrophilic interaction chromatography (HILIC), using both positive and negative electrospray ionization modes. The resulting data were analyzed using principal components analysis and partial least-squares-discriminant analysis. Discriminant models were developed for the response variables gender, diurnal variation, and age and were evaluated using external test sets to classify their predictive ability. The developed models using both positive and negative ionization mode data for reversed-phase and HILIC separations were very comparable, indicating that HILIC is a suitable method for increasing the fingerprint coverage for LC-MS metabonomic studies.     

Metabolic profiling of glucuronides in human urine by LC-MS/MS and partial least-squares discriminant analysis for classification and prediction of gender

Mass spectrometry (MS) is increasingly being used for metabolic profiling, but detection modes such as constant neutral loss or multiple reaction monitoring have not often been reported. These modes allow focusing on structurally related compounds, which could be advantageous for situations in which the trait under investigation is associated with a particular class of metabolites. In this study, we analyzed endogenous glucuronides excreted in human urine by monitoring characteristic transitions of putative steroid glucuronides by LC-MS/MS for discrimination of females from males. Two methods for data extraction were used: (i) a manual procedure based on visual inspection of the chromatograms and selection of 23 peaks and (ii) a software-supported method (MarkerView) set to extract 100 peaks. Data from 10 female and 10 male students were analyzed by principal component analysis (PCA) and partial least-squares discriminant analysis (PLS-DA) using software SIMCA. With PCA, only the manual peak selection resulted in clustering males and females. With PLS-DA, the manual method provided full separation on the basis of one single discriminant; the software-supported approach required a two-component model for complete separation. Loading plots were analyzed for their ability to reveal peaks with high discriminating power, that is, potential biomarkers. The PLS-DA models were validated with urine samples collected from five new females and five new males. Gender was correctly assigned for all. Our results indicate that inclusion of biological criteria for variable selection coupled to class-specific MS analysis and data extraction by appropriate software may constitute a valuable addition to the methods available for metabolomics.     

Quantitative measurement of JWH-018 and JWH-073 metabolites excreted in human urine

"K2/SPICE" products are commonly laced with aminoalkylindole synthetic cannabinoids (i.e., JWH-018 and JWH-073) and are touted as "legal" marijuana substitutes. Here we validate a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for measuring urinary concentrations of JWH-018, JWH-073, and several potential metabolites of each. The analytical procedure has high capacity for sample throughput and does not require solid phase or liquid extraction. Evaluation of human urine specimens collected after the subjects reportedly administered JWH-018 or a mixture of JWH-018 and JWH-073 provides preliminary evidence of clinical utility. Two subjects that consumed JWH-018 primarily excreted glucuronidated conjugates of 5-(3-(1-naphthoyl)-1H-indol-1-yl)-pentanoic acid (>30 ng/mL) and (1-(5-hydroxypentyl)-1H-indol-3-yl)(naphthalene-1-yl)-methanone (>50 ng/mL). Interestingly, oxidized metabolites of both JWH-018 and JWH-073 were detected in these specimens, suggesting either metabolic demethylation of JWH-018 to JWH-073 or a nonreported, previous JWH-073 exposure. Metabolic profiles generated from a subject who consumed a mixture of JWH-018 and JWH-073 were similar to profiles generated from subjects who presumably consumed JWH-018 exclusively. Oxidized metabolites of JWH-018 and JWH-073 were of the same pattern, but JWH-018 metabolites were excreted at lower concentrations. These results begin clinically validating the LC-MS/MS assay for detecting and quantifying aminoalkylindole metabolites. Full validation awaits further testing.     

Validated method for the determination of the ethanol consumption markers ethyl glucuronide, ethyl phosphate, and ethyl sulfate in human urine by reversed-phase/weak anion exchange liquid chromatography-tandem mass spectrometry

A liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for the determination of the ethanol consumption markers ethyl glucuronide (EtG), ethyl phosphate (EtP), and ethyl sulfate (EtS) in human urine was developed. A reversed-phase/weak anion exchange type stationary phase demonstrated particular suitability for the analysis of these highly polar acids. Sample preparation was minimized to centrifugation and dilution of urine prior to injection. The method was validated in the range of 5-750 microg.L(-1) with 1:20 and 1:1000 diluted urine, which corresponds to actual concentration ranges of 0.1-15 and 5-750 mg.L(-1) in undiluted samples. Method validation was carried out using six different lots of human urine. Over the entire calibration range intraday and interday precision (each n = 5, three concentration levels per dilution factor) adopted values between 0.6 and 4.7% and 0.8 and 12.1%, relative standard deviation, respectively. Corresponding accuracy values ranged between 94.2 and 113.5% and 86.6 and 110.9%, respectively. Matrix effects (absolute/relative) were found to be present in minor extent (approximately -30% to +15% MS/MS signal alterations) and were well corrected by the employed isotopically labeled internal standards. The validated assay was applied to urine samples of a drinking study as well as postmortem specimens. It was possible to assess the principal potential of EtP as ethanol consumption marker. Elevated concentration levels were found in real samples; however, EtP seems to be less sensitive compared to the previously known ethanol phase II conjugates EtG and EtS.     

Investigations of the effects of gender, diurnal variation, and age in human urinary metabolomic profiles

Metabolomics may have the capacity to revolutionize disease diagnosis through the identification of scores of metabolites that vary during environmental, pathogenic, or toxicological insult. NMR spectroscopy has become one of the main tools for measuring these changes since an NMR spectrum can accurately identify metabolites and their concentrations. The predominant approach in analyzing NMR data has been through the technique of spectral binning. However, identification of spectral areas in an NMR spectrum is insufficient for diagnostic evaluation, since it is unknown whether areas of interest are strictly caused by metabolic changes or are simply artifacts. In this paper, we explore differences in gender, diurnal variation, and age in a human population. We use the example of gender differences to compare traditional spectral binning techniques (NMR spectral areas) to novel targeted profiling techniques (metabolites and their concentrations). We show that targeted profiling produces robust models, generates accurate metabolite concentration data, and provides data that can be used to help understand metabolic differences in a healthy population. Metabolites relating to mitochondrial energy metabolism were found to differentiate gender and age. Dietary components and some metabolites related to circadian rhythms were found to differentiate time of day urine collection. The mechanisms by which these differences arise will be key to the discovery of new diagnostic tests and new understandings of the mechanism of disease.     

Two-dimensional method for time aligning liquid chromatography-mass spectrometry data

We describe a new time alignment method that takes advantage of both dimensions of LC-MS data to resolve ambiguities in peak matching while remaining computationally efficient. This approach, Warp2D, combines peak extraction with a two-dimensional correlation function to provide a reliable alignment scoring function that is insensitive to spurious peaks and background noise. One-dimensional alignment methods are often based on the total-ion-current elution profile of the spectrum and are unable to distinguish peaks of different masses. Our approach uses one-dimensional alignment in time, but with a scoring function derived from the overlap of peaks in two dimensions, thereby combining the specificity of two-dimensional methods with the computational performance of one-dimensional methods. The peaks are approximated as two-dimensional Gaussians of varying width. This approximation allows peak overlap (the measure of alignment quality) to be calculated analytically, without computationally intensive numerical integration in two dimensions. To demonstrate the general applicability of Warp2D, we chose a variety of complex samples that have substantial biological and analytical variability, including human serum and urine. We show that Warp2D works well with these diverse sample sets and with minimal tuning of parameters, based on the reduced standard deviation of peak elution times after warping. The combination of high computational speed, robustness with complex samples, and lack of need for detailed tuning makes this alignment method well suited to high-throughput LC-MS studies.     

Comprehensive profiling of histone modifications using a label-free approach and its applications in determining structure-function relationships

A two-pronged approach using specialized peptide detection and clustering tools was developed to profile changes in histone post-translational modifications (PTMs). The extent and nature of modification was inferred by comparing the mass profiles of intact core histones from nano LC-MS experiments. Histones displaying changes in their intact mass profiles were fractionated, derivatized with propionic anhydride, and digested with trypsin prior to nano LC-MS analyses. Our methodology was validated by comparing the abundance of histone PTMs in wild type and mutant strains of Saccharomyces cerevisiae lacking the histone acetyltransferase Rtt109 and a nucleosome assembly factor known as Asf1. Both Rtt109 and Asf1 were previously found to be essential for acetylation of histone H3 lysine 56 (H3K56ac), a modification that plays an important role in the response to genotoxic agents that interfere with DNA replication. The generation of ion abundance distribution plots enabled a rapid and comprehensive profiling of histone peptides. Our analytical methodology and data mining approach showed that most common histone PTMs were unaffected in mutant yeast cells lacking Rtt109 and Asf1. However, a subpopulation representing 17% of all H3 histones in wild type cells showed an acetylated K56 residue that was significantly reduced in both mutant strains. Our generic strategy for histone PTM profiling can be applied to assess global changes in histone PTMs across sample sets and to establish structure-function relationships.     

Identifying, evaluating, and controlling bioanalytical risks resulting from nonuniform matrix ion suppression/enhancement and nonlinear liquid chromatography-mass spectrometry assay response

Matrix ion suppression/enhancement is a well-observed and discussed phenomenon in electrospray ionization mass spectrometry. Nonuniform matrix ion suppression/enhancement across different types of samples in an analytical run is widely believed to be well compensated for by using a stable isotope-labeled internal standard (SIL-IS) in bioanalysis using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). Therefore, the risk of nonuniform matrix ion suppression/enhancement is usually deemed low when an SIL-IS is used. Here, we have identified, evaluated, and proposed solutions to control bioanalytical risks from nonuniform matrix ion suppression/enhancement even with an SIL-IS through a case study using omeprazole. Two lots of human blank urine were tested, and ion enhancement of about 500% for omeprazole was observed in one lot but not in the other. When a quadratic regression model had to be used, the assay failed the industry acceptance criteria due to unacceptable positive bias for the middle and high quality control (QC) samples. The failure was attributed to different extents of matrix ion enhancement between the standards (STDs) and QCs, which resulted in the misaligned results from the regression model. It was concluded that, for the same amount of drug, nonuniform ion enhancement for different types of samples (STD or QC) resulted in different ion intensities, therefore leading to different response behaviors (linear or nonlinear) at the mass spectrometer detector. A simplified mathematical model was used to evaluate the risk when unmatched response models occurred for different types of samples. A diagnostic factor Q (Q = X(ULOQ)(-A/B)) was proposed to monitor the risks, where X(ULOQ) is the upper limit of quantitation of the assay, A is the quadratic slope of the curve, and B is the linear slope of the curve. The potential maximum errors were estimated on the basis of the mathematical model for different scenarios, and Q values were given to control the risks under these conditions for bioanalysis using LC-MS/MS.     

High-throughput bioanalytical LC/MS/MS determination of benzodiazepines in human urine: 1000 samples per 12 hours.

The analytical capabilities of liquid chromatography tandem mass spectrometry for sensitive and highly selective determination of target compounds in complex biological samples makes it well suited for high-throughput analysis. We report the fast separation of six benzodiazepines isolated from human urine via selected reaction monitoring liquid chromatography/mass spectrometry using short dwell times to accommodate fast-eluting chromatographic peaks. The analytes were extracted from human urine samples along with their deuterium-labeled internal standards by a simple liquid-liquid extraction in 96-well plates. Using four autosamplers coupled to one chromatographic column and one tandem mass spectrometer operated in the turbo ion spray mode with positive ion detection, 1152 samples (12 96-well plates) were analyzed in less than 12 h. Through an electronic switching box designed and constructed in-house, the autosamplers were synchronized with the mass spectrometer so that injections were made as soon as the mass spectrometer was ready to collect data. Each run required 30 s to complete with another 7-8 s for the data system to load the next data file to be collected. Chromatographic integrity and ion current response remained relatively constant for the duration of the analyses. The results show acceptable precision and accuracy and demonstrate the feasibility of using fast separations with tandem mass spectrometry for high-throughout analysis of biological samples containing multiple analytes.     

Use of a microplate scintillation counter as a radioactivity detector for miniaturized separation techniques in drug metabolism

In miniaturized separation techniques, such as capillary electrophoresis (CE) or capillary liquid chromatography (LC), conventional on-line radioactivity detection of labeled compounds is restricted, because of insufficient sensitivity. It will be shown that a microplate scintillation counter for 96-well plates (TopCount) can be used as a sensitive and easy-to-handle radioactivity detector for capillary LC and CE. The attractive combination of capillary LC, eluent fractionation, and subsequent off-line counting is described. The new method is applied for rapid and sensitive separation and detection of 3H-labeled parent drug and its metabolites at levels between 25 and 700 cpm in rat urine. The advantages of capillary LC coupled to the TopCount, and combined with LC-MS data, can be of benefit in many analytical areas, including the characterization of metabolites at low concentration within complex biological fluids. With the same setup, the fractionation with subsequent off-line counting is equally applicable to CE. This is demonstrated with electrophoretically separated 14C-labeled impurities, nicely resolved from a negatively charged main compound, at low levels.