Detecting Regulatory Mechanisms in Endocrine Time Series Measurements

The regulatory mechanisms underlying pulsatile secretion are complex, especially as it is partly controlled by other hormones and the combined action of multiple agents. Regulatory relations between hormones are not directly observable but may be deduced from time series measurements of plasma hormone concentrations. Variation in plasma hormone levels are the resultant of secretion and clearance from the circulation. A strategy is proposed to extract inhibition, activation, thresholds and circadian synchronicity from concentration data, using particular association methods. Time delayed associations between hormone concentrations and/or extracted secretion pulse profiles reveal the information on regulatory mechanisms. The above mentioned regulatory mechanisms are illustrated with simulated data. Additionally, data from a lean cohort of healthy control subjects is used to illustrate activation (ACTH and cortisol) and circadian synchronicity (ACTH and TSH) in real data. The simulation and the real data both consist of 145 equidistant samples per individual, matching a 24-hr time span with 10 minute intervals. The results of the simulation and the real data are in concordance

Multivariate paired data analysis: multilevel PLSDA versus OPLSDA

Metabolomics data obtained from (human) nutritional intervention studies can have a rather complex structure that depends on the underlying experimental design. In this paper we discuss the complex structure in data caused by a cross-over designed experiment. In such a design, each subject in the study population acts as his or her own control and makes the data paired. For a single univariate response a paired t-test or repeated measures ANOVA can be used to test the differences between the paired observations. The same principle holds for multivariate data. In the current paper we compare a method that exploits the paired data structure in cross-over multivariate data (multilevel PLSDA) with a method that is often used by default but that ignores the paired structure (OPLSDA). The results from both methods have been evaluated in a small simulated example as well as in a genuine data set from a cross-over designed nutritional metabolomics study. It is shown that exploiting the paired data structure underlying the cross-over design considerably improves the power and the interpretability of the multivariate solution. Furthermore, the multilevel approach provides complementary information about (I) the diversity and abundance of the treatment effects within the different (subsets of) subjects across the study population, and (II) the intrinsic differences between these study subjects

Between Metabolite Relationships: an essential aspect of metabolic change

Not only the levels of individual metabolites, but also the relations between the levels of different metabolites may indicate (experimentally induced) changes in a biological system. Component analysis methods in current ‘standard’ use for metabolomics, such as Principal Component Analysis (PCA), do not focus on changes in these relations. We therefore propose the concept of ‘Between Metabolite Relationships’ (BMRs): common changes in the covariance (or correlation) between all metabolites in an organism. Such structural changes may indicate metabolic change brought about by experimental manipulation but which are lost with standard data analysis methods. These BMRs can be analysed by the INdividual Differences SCALing (INDSCAL) method. First the BMR quantification is described and subsequently the INDSCAL method. Finally, two studies illustrate the power and the applicability of BMRs in metabolomics. The first study is about the induced plant response of cabbage to herbivory, of which BMRs are a considerable part. In the second study—a human nutritional intervention study of green tea extract—standard data analysis tools did not reveal any metabolic change, although the BMRs were considerably affected. The presented results show that BMRs can be easily implemented in a wide variety of metabolomic studies. They provide a new source of information to describe biological systems in a way that fits flawlessly into the next generation of systems biology questions, dealing with personalized responses

Individual differences in metabolomics: individualised responses and between-metabolite relationships

Many metabolomics studies aim to find ‘biomarkers’: sets of molecules that are consistently elevated or decreased upon experimental manipulation. Biological effects, however, often manifest themselves along a continuum of individual differences between the biological replicates in the experiment. Such differences are overlooked or even diminished by methods in standard use for metabolomics, although they may contain a wealth of information on the experiment. Properly understanding individual differences is crucial for generating knowledge in fields like personalised medicine, evolution and ecology. We propose to use simultaneous component analysis with individual differences constraints (SCA-IND), a data analysis method from psychology that focuses on these differences. This method constructs axes along the natural biochemical differences between biological replicates, comparable to principal components. The model may shed light on changes in the individual differences between experimental groups, but also on whether these differences correspond to, e.g., responders and non-responders or to distinct chemotypes. Moreover, SCA-IND reveals the individuals that respond most to a manipulation and are best suited for further experimentation. The method is illustrated by the analysis of individual differences in the metabolic response of cabbage plants to herbivory. The model reveals individual differences in the response to shoot herbivory, where two ‘response chemotypes’ may be identified. In the response to root herbivory the model shows that individual plants differ strongly in response dynamics. Thereby SCA-IND provides a hitherto unavailable view on the chemical diversity of the induced plant response, that greatly increases understanding of the system

Quantitative proteomics analysis of an ethanol- and a lactate-producing mutant strain of Synechocystis sp. PCC6803

BACKGROUND: This study aimed at exploring the molecular physiological consequences of a major redirection of carbon flow in so-called cyanobacterial cell factories: quantitative whole-cell proteomics analyses were carried out on two (14)N-labelled Synechocystis mutant strains, relative to their (15)N-labelled wild-type counterpart. Each mutant strain overproduced one specific commodity product, i.e. ethanol or lactic acid, to such an extent that the majority of the incoming CO2 in the organism was directly converted into the product. RESULTS: In total, 267 proteins have been identified with a significantly up- or down-regulated expression level. In the ethanol-producing mutant, which had the highest relative direct flux of carbon-to-product (>65%), significant up-regulation of several components involved in the initial stages of CO2 fixation for cellular metabolism was detected. Also a general decrease in abundance of the protein synthesizing machinery of the cells and a specific induction of an oxidative stress response were observed in this mutant. In the lactic acid overproducing mutant, that expresses part of the heterologous l-lactate dehydrogenase from a self-replicating plasmid, specific activation of two CRISPR associated proteins, encoded on the endogenous pSYSA plasmid, was observed. RT-qPCR was used to measure, of nine of the genes identified in the proteomics studies, also the adjustment of the corresponding mRNA level. CONCLUSION: The most striking adjustments detected in the proteome of the engineered cells were dependent on the specific product formed, with, e.g. more stress caused by lactic acid- than by ethanol production. Up-regulation of the total capacity for CO2 fixation in the ethanol-producing strain was due to hierarchical- rather than metabolic regulation. Furthermore, plasmid-based expression of heterologous gene(s) may induce genetic instability. For selected, limited, number of genes a striking correlation between the respective mRNA- and the corresponding protein expression level was observed, suggesting that for the expression of these genes regulation takes place primarily at the level of gene transcription

UvA-DARE (Digital Academic Repository) A classification model for the Leiden proteomics competition A Classification Model for the Leiden Proteomics Competition A Classification Model for the Leiden Proteomics Competition

Abstract A strategy is presented to build a discrimination model in proteomics studies. The model is built using cross-validation. This cross-validation step can simply be combined with a variable selection method, called rank products. The strategy is especially suitable for the low-samplesto-variables-ratio (undersampling) case, as is often encountered in proteomics and metabolomics studies. As a classification method, Principal Component Discriminant Analysis is used; however, the methodology can be used with any classifier. A data set containing serum samples from breast cancer patients and healthy controls is analysed. Double cross-validation shows that the sensitivity of the model is 82% and the specificity 86%. Potential putative biomarkers are identified using the variable selection method. In each cross-validation loop a classification model is built. The final classification uses a majority voting scheme from the ensemble classifier