Abstract: ABSTRACT Changes in seasonal photoperiod provides an important environmental signal that affects the timing of winter dormancy in perennial, deciduous, temperate tree species, such as hybrid aspen (Populus tremula x Populus tremuloides). In this species, growth cessation, cold acclimation and dormancy are induced in the autumn by the detection of day-length shortening that occurs at a given critical day length. Important components in the detection of such day-length changes are photoreceptors and the circadian clock, and many plant responses at both the gene regulation and metabolite levels are expected to be diurnal. To directly examine this expectation and study components in these events, here we report transcriptomic and metabolomic responses to a change in photoperiod from long to short days in hybrid aspen. We found about 16% of genes represented on the arrays to be diurnally regulated, as assessed by our pre-defined criteria. Furthermore, several of these genes were involved in circadian-associated processes, including photosynthesis and primary and secondary metabolism. Metabolites affected by the change in photoperiod were mostly involved in carbon metabolism. Taken together, we have thus established a molecular catalog of events that precede a response to winter.
Abstract: The investigation presented here describes a protocol designed to perform high-throughput metabolic profiling analysis on human blood plasma by ultra-performance liquid chromatography/mass spectrometry (UPLC/MS). To address whether a previous extraction protocol for gas chromatography (GC)/MS-based metabolic profiling of plasma could be used for UPLC/MS-based analysis, the original protocol was compared with similar methods for extraction of low-molecular-weight compounds from plasma via protein precipitation. Differences between extraction methods could be observed, but the previously published extraction method was considered the best. UPLC columns with three different stationary phases (C8, C18, and phenyl) were used in identical experimental runs consisting of a total of 60 injections of extracted male and female plasma samples. The C8 column was determined to be the best for metabolic profiling analysis on plasma. The acquired UPLC/MS data of extracted male and female plasma samples was subjected to principal component analysis (PCA) and orthogonal projections to latent structures discriminant analysis (OPLS-DA). Furthermore, a strategy for compound identification was applied here, demonstrating the strength of high-mass-accuracy time-of-flight (TOF)/MS analysis in metabolic profiling.
Abstract: A strategy for processing of metabolomic GC/MS data is presented. By considering the relationship between quantity and quality of detected profiles, representative data suitable for multiple sample comparisons and metabolite identification was generated. Design of experiments (DOE) and multivariate analysis was used to relate the changes in settings of the hierarchical multivariate curve resolution (H-MCR) method to quantitative and qualitative characteristics of the output data. These characteristics included number of resolved profiles, chromatographic quality in terms of reproducibility between analytical replicates, and spectral quality defined by purity and number of spectra containing structural information. The strategy was exemplified in two datasets: one containing 119 common metabolites, 18 of which were varied according to a DOE protocol; and one consisting of rat urine samples from control rats and rats exposed to a liver toxin. It was shown that the performance of the data processing could be optimized to produce metabolite data of high quality that allowed reliable sample comparisons and metabolite identification. This is a general approach applicable to any type of data processing where the important processing parameters are known and relevant output data characteristics can be defined. The results imply that this type of data quality optimization should be carried out as an integral step of data processing to ensure high quality data for further modeling and biological evaluation. Within metabolomics, this degree of optimization will be of high importance to generate models and extract biomarkers or biomarker patterns of biological or clinical relevance.
Abstract: We have developed a comprehensive method combining analytical techniques of one-dimensional (1D) and two-dimensional (GC x GC) gas chromatography-time-of-flight (TOF)-mass spectrometry. This method was applied to the metabolic phenotyping of natural variants in rice for the 68 world rice core collection (WRC) and two other varieties. Ten metabolites were selected as metabolite representatives, and the selected ion current of each metabolite peak obtained from both techniques were statistically compared. Our method of combining 1D- and GC x GC-TOF/MS is useful for the metabolic phenotyping of natural variants in rice for further studies in breeding programs.
Abstract: A novel hypothesis-free multivariate screening methodology for the study of human exercise metabolism in blood serum is presented. Serum gas chromatography/time-of-flight mass spectrometry (GC/TOFMS) data was processed using hierarchical multivariate curve resolution (H-MCR), and orthogonal partial least-squares discriminant analysis (OPLS-DA) was used to model the systematic variation related to the acute effect of strenuous exercise. Potential metabolic biomarkers were identified using data base comparisons. Extensive validation was carried out including predictive H-MCR, 7-fold full cross-validation, and predictions for the OPLS-DA model, variable permutation for highlighting interesting metabolites, and pairwise t tests for examining the significance of metabolites. The concentration changes of potential biomarkers were verified in the raw GC/TOFMS data. In total, 420 potential metabolites were resolved in the serum samples. On the basis of the relative concentrations of the 420 resolved metabolites, a valid multivariate model for the difference between pre- and post-exercise subjects was obtained. A total of 34 metabolites were highlighted as potential biomarkers, all statistically significant (p < 8.1E-05). As an example, two potential markers were identified as glycerol and asparagine. The concentration changes for these two metabolites were also verified in the raw GC/TOFMS data.The strategy was shown to facilitate interpretation and validation of metabolic interactions in human serum as well as revealing the identity of potential markers for known or novel mechanisms of human exercise physiology. The multivariate way of addressing metabolism studies can help to increase the understanding of the integrative biology behind, as well as unravel new mechanistic explanations in relation to, exercise physiology.
Abstract: BACKGROUND: Metabolites are not only the catalytic products of enzymatic reactions but also the active regulators or the ultimate phenotype of metabolic homeostasis in highly complex cellular processes. The modes of regulation at the metabolome level can be revealed by metabolic networks. We investigated the metabolic network between wild-type and 2 mutant (methionine-over accumulation 1 [mto1] and transparent testa4 [tt4]) plants regarding the alteration of metabolite accumulation in Arabidopsis thaliana. RESULTS: In the GC-TOF/MS analysis, we acquired quantitative information regarding over 170 metabolites, which has been analyzed by a novel score (ZMC, z-score of metabolite correlation) describing a characteristic metabolite in terms of correlation. Although the 2 mutants revealed no apparent morphological abnormalities, the overall correlation values in mto1 were much lower than those of the wild-type and tt4 plants, indicating the loss of overall network stability due to the uncontrolled accumulation of methionine. In the tt4 mutant, a new correlation between malate and sinapate was observed although the levels of malate, sinapate, and sinapoylmalate remain unchanged, suggesting an adaptive reconfiguration of the network. Gene-expression correlations presumably responsible for these metabolic networks were determined using the metabolite correlations as clues. CONCLUSION: Two Arabidopsis mutants, mto1 and tt4, exhibited the following changes in entire metabolome networks: the overall loss of metabolic stability (mto1) or the generation of a metabolic network of a backup pathway for the lost physiological functions (tt4). The expansion of metabolite correlation to gene-expression correlation provides detailed insights into the systemic understanding of the plant cellular process regarding metabolome and transcriptome.
Abstract: A method for predictive metabolite profiling based on resolution of GC-MS data followed by multivariate data analysis is presented and applied to three different biofluid data sets (rat urine, aspen leaf extracts, and human blood plasma). Hierarchical multivariate curve resolution (H-MCR) was used to simultaneously resolve the GC-MS data into pure profiles, describing the relative metabolite concentrations between samples, for multivariate analysis. Here, we present an extension of the H-MCR method allowing treatment of independent samples according to processing parameters estimated from a set of training samples. Predictions or inclusion of the new samples, based on their metabolite profiles, into an existing model could then be carried out, which is a requirement for a working application within, e.g., clinical diagnosis. Apart from allowing treatment and prediction of independent samples the proposed method also reduces the time for the curve resolution process since only a subset of representative samples have to be processed while the remaining samples can be treated according to the obtained processing parameters. The time required for resolving the 30 training samples in the rat urine example was approximately 13 h, while the treatment of the 30 test samples according to the training parameters required only approximately 30 s per sample (approximately 15 min in total). In addition, the presented results show that the suggested approach works for describing metabolic changes in different biofluids, indicating that this is a general approach for high-throughput predictive metabolite profiling, which could have important applications in areas such as plant functional genomics, drug toxicity, treatment efficacy and early disease diagnosis.
Abstract: A multivariate strategy for studying the metabolic response over time in urinary GC/MS data is presented and exemplified by a study of drug-induced liver toxicity in the rat. The strategy includes the generation of representative data through hierarchical multivariate curve resolution (H-MCR), highlighting the importance of obtaining resolved metabolite profiles for quantification and identification of exogenous (drug related) and endogenous compounds (potential biomarkers) and for allowing reliable comparisons of multiple samples through multivariate projections. Batch modelling was used to monitor and characterize the normal (control) metabolic variation over time as well as to map the dynamic response of the drug treated animals in relation to the control. In this way treatment related metabolic responses over time could be detected and classified as being drug related or being potential biomarkers. In summary the proposed strategy uses the relatively high sensitivity and reproducibility of GC/MS in combination with efficient multivariate curve resolution and data analysis to discover individual markers of drug metabolism and drug toxicity. The presented results imply that the strategy can be of great value in drug toxicity studies for classifying metabolic markers in relation to their dynamic responses as well as for biomarker identification.
Abstract: A strategy for robust and reliable mechanistic statistical modelling of metabolic responses in relation to drug induced toxicity is presented. The suggested approach addresses two cases commonly occurring within metabonomic toxicology studies, namely; 1) A pre-defined hypothesis about the biological mechanism exists and 2) No such hypothesis exists. GC/MS data from a liver toxicity study consisting of rat urine from control rats and rats exposed to a proprietary AstraZeneca compound were resolved by means of hierarchical multivariate curve resolution (H-MCR) generating 287 resolved chromatographic profiles with corresponding mass spectra. Filtering according to significance in relation to drug exposure rendered in 210 compound profiles, which were subjected to further statistical analysis following correction to account for the control variation over time. These dose related metabolite traces were then used as new observations in the subsequent analyses. For case 1, a multivariate approach, named Target Batch Analysis, based on OPLS regression was applied to correlate all metabolite traces to one or more key metabolites involved in the pre-defined hypothesis. For case 2, principal component analysis (PCA) was combined with hierarchical cluster analysis (HCA) to create a robust and interpretable framework for unbiased mechanistic screening. Both the Target Batch Analysis and the unbiased approach were cross-verified using the other method to ensure that the results did match in terms of detected metabolite traces. This was also the case, implying that this is a working concept for clustering of metabolites in relation to their toxicity induced dynamic profiles regardless if there is a pre-existing hypothesis or not. For each of the methods the detected metabolites were subjected to identification by means of data base comparison as well as verification in the raw data. The proposed strategy should be seen as a general approach for facilitating mechanistic modelling and interpretations in metabolomic studies.
Abstract: LC/MS is an analytical technique that, due to its high sensitivity, has become increasingly popular for the generation of metabolic signatures in biological samples and for the building of metabolic data bases. However, to be able to create robust and interpretable (transparent) multivariate models for the comparison of many samples, the data must fulfil certain specific criteria: (i) that each sample is characterized by the same number of variables, (ii) that each of these variables is represented across all observations, and (iii) that a variable in one sample has the same biological meaning or represents the same metabolite in all other samples. In addition, the obtained models must have the ability to make predictions of, e.g. related and independent samples characterized accordingly to the model samples. This method involves the construction of a representative data set, including automatic peak detection, alignment, setting of retention time windows, summing in the chromatographic dimension and data compression by means of alternating regression, where the relevant metabolic variation is retained for further modelling using multivariate analysis. This approach has the advantage of allowing the comparison of large numbers of samples based on their LC/MS metabolic profiles, but also of creating a means for the interpretation of the investigated biological system. This includes finding relevant systematic patterns among samples, identifying influential variables, verifying the findings in the raw data, and finally using the models for predictions. The presented strategy was here applied to a population study using urine samples from two cohorts, Shanxi (People's Republic of China) and Honolulu (USA). The results showed that the evaluation of the extracted information data using partial least square discriminant analysis (PLS-DA) provided a robust, predictive and transparent model for the metabolic differences between the two populations. The presented findings suggest that this is a general approach for data handling, analysis, and evaluation of large metabolic LC/MS data sets.
Abstract: In metabolomics, the objective is to identify differences in metabolite profiles between samples. A widely used tool in metabolomics investigations is gas chromatography-mass spectrometry (GC/MS). More than 400 compounds can be detected in a single analysis, if overlapping GC/MS peaks are deconvoluted. However, the deconvolution process is time-consuming and difficult to automate, and additional processing is needed in order to compare samples. Therefore, there is a need to improve and automate the data processing strategy for data generated in GC/MS-based metabolomics; if not, the processing step will be a major bottleneck for high-throughput analyses. Here we describe a new semiautomated strategy using a hierarchical multivariate curve resolution approach that processes all samples simultaneously. The presented strategy generates (after appropriate treatment, e.g., multivariate analysis) tables of all the detected metabolites that differ in relative concentrations between samples. The processing of 70 samples took similar time to that of the GC/TOFMS analyses of the samples. The strategy has been validated using two different sets of samples: a complex mixture of standard compounds and Arabidopsis samples.
Abstract: Analysis of the entire set of low molecular weight compounds (LMC), the metabolome, could provide deeper insights into mechanisms of disease and novel markers for diagnosis. In the investigation, we developed an extraction and derivatization protocol, using experimental design theory (design of experiment), for analyzing the human blood plasma metabolome by GC/MS. The protocol was optimized by evaluating the data for more than 500 resolved peaks using multivariate statistical tools including principal component analysis and partial least-squares projections to latent structures (PLS). The performance of five organic solvents (methanol, ethanol, acetonitrile, acetone, chloroform), singly and in combination, was investigated to optimize the LMC extraction. PLS analysis demonstrated that methanol extraction was particularly efficient and highly reproducible. The extraction and derivatization conditions were also optimized. Quantitative data for 32 endogenous compounds showed good precision and linearity. In addition, the determined amounts of eight selected compounds agreed well with analyses by independent methods in accredited laboratories, and most of the compounds could be detected at absolute levels of approximately 0.1 pmol injected, corresponding to plasma concentrations between 0.1 and 1 microM. The results suggest that the method could be usefully integrated into metabolomic studies for various purposes, e.g., for identifying biological markers related to diseases.
Abstract: Here we present an approach for on-line control and monitoring of pulpwood chip properties based on near infrared (NIR) spectroscopy and multivariate data analysis. In addition, this paper suggests how to deal with large multivariate data sets in order to extract information which can be used as a basis for changes in raw material or process conditions in the drive towards more optimal intermediate or end product properties within the pulp and paper industry. The pulpwood chips used as raw material in a pulp and paper making process were characterized at- and on-line using NIR spectroscopic measurements. Collected NIR spectra were used in multivariate calibration models for prediction of the moisture content as well as the between- and within-species variation in the studied raw material. Statistical experimental design was used to form a calibration data set including most of the variation occurring in a real on-line situation. NIR spectra for all designed samples were measured at-line and the estimated calibration models were used for carrying out predictions on-line. Predictions of the moisture content (% dry weight) as well as the percentage contents of pine and sawmill chips in the raw material were carried out using partial least squares projections to latent structures (PLS) methodology. NIR spectra were collected subsequently on-line once every minute, and, to reduce the problem with noise in the time series predictions, the measured signals were filtered using a moving average of 100 predicted values. This provided smoother predictions more suitable for process monitoring and control. To validate the quality of the predictions, wood chips from the studied process were sampled and analysed in the laboratory before being subjected to predictions in the on-line model. Comparison of the filtered on-line predictions with the results obtained from the laboratory measurements indicated that moisture and pine chip contents could be well predicted by the on-line model, while predictions of sawmill chip content showed less promising results.
Abstract: The usual aim in metabolomic studies is to quantify the entire metabolome of each of a series of biological samples. To do this for complex biological matrices, e.g., plant tissues, efficient and reproducible extraction protocols must be developed. However, derivatization protocols must also be developed if GC/MS (one of the mostly widely used analytical methods for metabolomics) is involved. The aim of this study was to investigate how different chemical and physical factors (extraction solvent, derivatization reagents, and temperature) affect the extraction and derivatization of the metabolome from leaves of the plant Arabidopsis thaliana. Using design of experiment procedures, variation was systematically introduced, and the effects of this variation were analyzed using regression models. The results show that this approach allows a reliable protocol for metabolomic analysis of Arabidopsis to be determined with a relatively limited number of experiments. Following two different investigations an extraction and derivatization protocol was chosen. Further, the reproducibility of the analysis of 66 endogenous compounds was investigated, and it was shown that both hydrophilic and lipophilic compounds were detected with high reproducibility.
Abstract: In metabolomics, the purpose is to identify and quantify all the metabolites in a biological system. Combined gas chromatography and mass spectrometry (GC/MS) is one of the most commonly used techniques in metabolomics together with 1H NMR, and it has been shown that more than 300 compounds can be distinguished with GC/MS after deconvolution of overlapping peaks. To avoid having to deconvolute all analyzed samples prior to multivariate analysis of the data, we have developed a strategy for rapid comparison of nonprocessed MS data files. The method includes baseline correction, alignment, time window determinations, alternating regression, PLS-DA, and identification of retention time windows in the chromatograms that explain the differences between the samples. Use of alternating regression also gives interpretable loadings, which retain the information provided by m/z values that vary between the samples in each retention time window. The method has been applied to plant extracts derived from leaves of different developmental stages and plants subjected to small changes in day length. The data show that the new method can detect differences between the samples and that it gives results comparable to those obtained when deconvolution is applied prior to the multivariate analysis. We suggest that this method can be used for rapid comparison of large sets of GC/MS data, thereby applying time-consuming deconvolution only to parts of the chromatograms that contribute to explain the differences between the samples.
Abstract: One trend in the "omics" sciences is the generation of increasing amounts of data, describing complex biological samples. To cope with this and facilitate progress towards reliable diagnostic tools, it is crucial to develop methods for extracting representative and predictive information. In global metabolite analysis (metabolomics and metabonomics) NMR, GC/MS and LC/MS are the main platforms for data generation. Multivariate projection methods (e.g. PCA, PLS and O-PLS) have been recognized as efficient tools for data analysis within subjects such as biology and chemistry due to their ability to provide interpretable models based on many, correlated variables. In global metabolite analysis, these methods have been successfully applied in areas such as toxicology, disease diagnosis and plant functional genomics.
This thesis describes the development of processing methods for the unbiased extraction of representative and predictive information from metabolic GC/MS and LC/MS data characterizing biofluids, e.g. plant extracts, urine and blood plasma. In order to allow the multivariate projections to detect and highlight differences between samples, one requirement of the processing methods is that they must extract a common set of descriptors from all samples and still retain the metabolically relevant information in the data. In Papers I and II this was done by applying a hierarchical multivariate compression approach to both GC/MS and LC/MS data. In the study described in Paper III a hierarchical multivariate curve resolution strategy (H-MCR) was developed for simultaneously resolving multiple GC/MS samples into pure profiles. In Paper IV the H-MCR method was applied to a drug toxicity study in rats, where the method?s potential for biomarker detection and identification was exemplified. Finally, the H-MCR method was extended, as described in Paper V, allowing independent samples to be processed and predicted using a model based on an existing set of representative samples. The fact that these processing methods proved to be valid for predicting the properties of new independent samples indicates that it is now possible for global metabolite analysis to be extended beyond isolated studies. In addition, the results facilitate high through-put analysis, because predicting the nature of samples is rapid compared to the actual processing. In summary this research highlights the possibilities for using global metabolite analysis in diagnosis.
