Pattern Recognition and Pathway Analysis with Genetic Algorithms in Mass Spectrometry Based Metabolomics

A unique volatile signature distinguishes malaria infection from other conditions that cause similar symptoms

Recent findings suggest that changes in human odors caused by malaria infection have significant potential as diagnostic biomarkers. However, uncertainty remains regarding the specificity of such biomarkers, particularly in populations where many different pathological conditions may elicit similar symptoms. We explored the ability of volatile biomarkers to predict malaria infection status in Kenyan schoolchildren exhibiting a range of malaria-like symptoms. Using genetic algorithm models to explore data from skin volatile collections, we were able to identify malaria infection with 100% accuracy among children with fever and 75% accuracy among children with other symptoms. While we observed characteristic changes in volatile patterns driven by symptomatology, our models also identified malaria-specific biomarkers with robust predictive capability even in the presence of other pathogens that elicit similar symptoms.

Automatic peak assignment and visualisation of copolymer mass spectrometry data using the 'genetic algorithm'

Copolymer analysis is vitally important as the materials have a wide variety of applications due to their tunable properties. Processing mass spectrometry data for copolymer samples can be very complex due to the increase in the number of species when the polymer chains are formed by two or more monomeric units. In this paper, we describe the use of the genetic algorithm for automated peak assignment of copolymers synthesised by a variety of polymerisation methods. We find that in using this method we are able to easily assign copolymer spectra in a few minutes and visualise them into heat maps. These heat maps allow us to look qualitatively at the distribution of the chains, by showing how they alter with different polymerisation techniques, and by changing the initial copolymer composition. This methodology is simple to use and requires little user input, which makes it well suited for use by less expert users. The data outputted by the automatic assignment may also allow for more complex data processing in the future.

Metabolomics analysis of metabolic effects of nicotinamide phosphoribosyltransferase (NAMPT) inhibition on human cancer cells

Nicotinamide phosphoribosyltransferase (NAMPT) plays an important role in cellular bioenergetics. It is responsible for converting nicotinamide to nicotinamide adenine dinucleotide, an essential molecule in cellular metabolism. NAMPT has been extensively studied over the past decade due to its role as a key regulator of nicotinamide adenine dinucleotide–consuming enzymes. NAMPT is also known as a potential target for therapeutic intervention due to its involvement in disease. In the current study, we used a global mass spectrometry–based metabolomic approach to investigate the effects of FK866, a small molecule inhibitor of NAMPT currently in clinical trials, on metabolic perturbations in human cancer cells. We treated A2780 (ovarian cancer) and HCT-116 (colorectal cancer) cell lines with FK866 in the presence and absence of nicotinic acid. Significant changes were observed in the amino acids metabolism and the purine and pyrimidine metabolism. We also observed metabolic alterations in glycolysis, the citric acid cycle (TCA), and the pentose phosphate pathway. To expand the range of the detected polar metabolites and improve data confidence, we applied a global metabolomics profiling platform by using both non-targeted and targeted hydrophilic (HILIC)-LC-MS and GC-MS analysis. We used Ingenuity Knowledge Base to facilitate the projection of metabolomics data onto metabolic pathways. Several metabolic pathways showed differential responses to FK866 based on several matches to the list of annotated metabolites. This study suggests that global metabolomics can be a useful tool in pharmacological studies of the mechanism of action of drugs at a cellular level.

Nucleotides released from palmitate-challenged muscle cells through pannexin-3 attract monocytes

Obesity-associated low-grade inflammation in metabolically relevant tissues contributes to insulin resistance. We recently reported monocyte/macrophage infiltration in mouse and human skeletal muscles. However, the molecular triggers of this infiltration are unknown, and the role of muscle cells in this context is poorly understood. Animal studies are not amenable to the specific investigation of this vectorial cellular communication. Using cell cultures, we investigated the crosstalk between myotubes and monocytes exposed to physiological levels of saturated and unsaturated fatty acids. Media from L6 myotubes treated with palmitate-but not palmitoleate-induced THP1 monocyte migration across transwells. Palmitate activated the Toll-like receptor 4 (TLR4)/nuclear factor-κB (NF-κB) pathway in myotubes and elevated cytokine expression, but the monocyte chemoattracting agent was not a polypeptide. Instead, nucleotide degradation eliminated the chemoattracting properties of the myotube-conditioned media. Moreover, palmitate-induced expression and activity of pannexin-3 channels in myotubes were mediated by TLR4-NF-κB, and TLR4-NF-κB inhibition or pannexin-3 knockdown prevented monocyte chemoattraction. In mice, the expression of pannexin channels increased in adipose tissue and skeletal muscle in response to high-fat feeding. These findings identify pannexins as new targets of saturated fatty acid-induced inflammation in myotubes, and point to nucleotides as possible mediators of immune cell chemoattraction toward muscle in the context of obesity.

A comprehensive workflow of mass spectrometry-based untargeted metabolomics in cancer metabolic biomarker discovery using human plasma and urine

Current available biomarkers lack sensitivity and/or specificity for early detection of cancer. To address this challenge, a robust and complete workflow for metabolic profiling and data mining is described in details. Three independent and complementary analytical techniques for metabolic profiling are applied: hydrophilic interaction liquid chromatography (HILIC-LC), reversed-phase liquid chromatography (RP-LC), and gas chromatography (GC). All three techniques are coupled to a mass spectrometer (MS) in the full scan acquisition mode, and both unsupervised and supervised methods are used for data mining. The univariate and multivariate feature selection are used to determine subsets of potentially discriminative predictors. These predictors are further identified by obtaining accurate masses and isotopic ratios using selected ion monitoring (SIM) and data-dependent MS/MS and/or accurate mass MSn ion tree scans utilizing high resolution MS. A list combining all of the identified potential biomarkers generated from different platforms and algorithms is used for pathway analysis. Such a workflow combining comprehensive metabolic profiling and advanced data mining techniques may provide a powerful approach for metabolic pathway analysis and biomarker discovery in cancer research. Two case studies with previous published data are adapted and included in the context to elucidate the application of the workflow.

Proteomics and metabolomics analysis of a trait anxiety mouse model reveals divergent mitochondrial pathways

BACKGROUND

Although anxiety disorders are the most prevalent psychiatric disorders, no molecular biomarkers exist for their premorbid diagnosis, accurate patient subcategorization, or treatment efficacy prediction. To unravel the neurobiological underpinnings and identify candidate biomarkers and affected pathways for anxiety disorders, we interrogated the mouse model of high anxiety-related behavior (HAB), normal anxiety-related behavior (NAB), and low anxiety-related behavior (LAB) employing a quantitative proteomics and metabolomics discovery approach.

METHODS

We compared the cingulate cortex synaptosome proteomes of HAB and LAB mice by in vivo (15)N metabolic labeling and mass spectrometry and quantified the cingulate cortex metabolomes of HAB/NAB/LAB mice. The combined data sets were used to identify divergent protein and metabolite networks by in silico pathway analysis. Selected differentially expressed proteins and affected pathways were validated with immunochemical and enzymatic assays.

RESULTS

Altered levels of up to 300 proteins and metabolites were found between HAB and LAB mice. Our data reveal alterations in energy metabolism, mitochondrial import and transport, oxidative stress, and neurotransmission, implicating a previously nonhighlighted role of mitochondria in modulating anxiety-related behavior.

CONCLUSIONS

Our results offer insights toward a molecular network of anxiety pathophysiology with a focus on mitochondrial contribution and provide the basis for pinpointing affected pathways in anxiety-related behavior.

Advances in structure elucidation of small molecules using mass spectrometry

The structural elucidation of small molecules using mass spectrometry plays an important role in modern life sciences and bioanalytical approaches. This review covers different soft and hard ionization techniques and figures of merit for modern mass spectrometers, such as mass resolving power, mass accuracy, isotopic abundance accuracy, accurate mass multiple-stage MS(n) capability, as well as hybrid mass spectrometric and orthogonal chromatographic approaches. The latter part discusses mass spectral data handling strategies, which includes background and noise subtraction, adduct formation and detection, charge state determination, accurate mass measurements, elemental composition determinations, and complex data-dependent setups with ion maps and ion trees. The importance of mass spectral library search algorithms for tandem mass spectra and multiple-stage MS(n) mass spectra as well as mass spectral tree libraries that combine multiple-stage mass spectra are outlined. The successive chapter discusses mass spectral fragmentation pathways, biotransformation reactions and drug metabolism studies, the mass spectral simulation and generation of in silico mass spectra, expert systems for mass spectral interpretation, and the use of computational chemistry to explain gas-phase phenomena. A single chapter discusses data handling for hyphenated approaches including mass spectral deconvolution for clean mass spectra, cheminformatics approaches and structure retention relationships, and retention index predictions for gas and liquid chromatography. The last section reviews the current state of electronic data sharing of mass spectra and discusses the importance of software development for the advancement of structure elucidation of small molecules. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12566-010-0015-9) contains supplementary material, which is available to authorized users.

Deducing corticotropin-releasing hormone receptor type 1 signaling networks from gene expression data by usage of genetic algorithms and graphical Gaussian models

BACKGROUND

Dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis is a hallmark of complex and multifactorial psychiatric diseases such as anxiety and mood disorders. About 50-60% of patients with major depression show HPA axis dysfunction, i.e. hyperactivity and impaired negative feedback regulation. The neuropeptide corticotropin-releasing hormone (CRH) and its receptor type 1 (CRHR1) are key regulators of this neuroendocrine stress axis. Therefore, we analyzed CRH/CRHR1-dependent gene expression data obtained from the pituitary corticotrope cell line AtT-20, a well-established in vitro model for CRHR1-mediated signal transduction. To extract significantly regulated genes from a genome-wide microarray data set and to deduce underlying CRHR1-dependent signaling networks, we combined supervised and unsupervised algorithms.

RESULTS

We present an efficient variable selection strategy by consecutively applying univariate as well as multivariate methods followed by graphical models. First, feature preselection was used to exclude genes not differentially regulated over time from the dataset. For multivariate variable selection a maximum likelihood (MLHD) discriminant function within GALGO, an R package based on a genetic algorithm (GA), was chosen. The topmost genes representing major nodes in the expression network were ranked to find highly separating candidate genes. By using groups of five genes (chromosome size) in the discriminant function and repeating the genetic algorithm separately four times we found eleven genes occurring at least in three of the top ranked result lists of the four repetitions. In addition, we compared the results of GA/MLHD with the alternative optimization algorithms greedy selection and simulated annealing as well as with the state-of-the-art method random forest. In every case we obtained a clear overlap of the selected genes independently confirming the results of MLHD in combination with a genetic algorithm. With two unsupervised algorithms, principal component analysis and graphical Gaussian models, putative interactions of the candidate genes were determined and reconstructed by literature mining. Differential regulation of six candidate genes was validated by qRT-PCR.

CONCLUSIONS

The combination of supervised and unsupervised algorithms in this study allowed extracting a small subset of meaningful candidate genes from the genome-wide expression data set. Thereby, variable selection using different optimization algorithms based on linear classifiers as well as the nonlinear random forest method resulted in congruent candidate genes. The calculated interacting network connecting these new target genes was bioinformatically mapped to known CRHR1-dependent signaling pathways. Additionally, the differential expression of the identified target genes was confirmed experimentally.