Ion-Pairing Chromatography and Amine Derivatization Provide Complementary Approaches for the Targeted LC-MS Analysis of the Polar Metabolome.

Liquid chromatography coupled to mass spectrometry is a key metabolomics/metabonomics technology. Reversed-phase liquid chromatography (RPLC) is very widely used as a separation step, but typically has poor retention of highly polar metabolites. Here, we evaluated the combination of two alternative methods for improving retention of polar metabolites based on 6-aminoquinoloyl-N-hydroxysuccinidimyl carbamate derivatization for amine groups, and ion-pairing chromatography (IPC) using tributylamine as an ion-pairing agent to retain acids. We compared both of these methods to RPLC and also to each other, for targeted analysis using a triple-quadrupole mass spectrometer, applied to a library of ca. 500 polar metabolites. IPC and derivatization were complementary in terms of their coverage: combined, they improved the proportion of metabolites with good retention to 91%, compared to just 39% for RPLC alone. The combined method was assessed by analyzing a set of liver extracts from aged male and female mice that had been treated with the polyphenol compound ampelopsin. Not only were a number of significantly changed metabolites detected, but also it could be shown that there was a clear interaction between ampelopsin treatment and sex, in that the direction of metabolite change was opposite for males and females.

Pathway analysis in metabolomics: Recommendations for the use of over-representation analysis.

Over-representation analysis (ORA) is one of the commonest pathway analysis approaches used for the functional interpretation of metabolomics datasets. Despite the widespread use of ORA in metabolomics, the community lacks guidelines detailing its best-practice use. Many factors have a pronounced impact on the results, but to date their effects have received little systematic attention. Using five publicly available datasets, we demonstrated that changes in parameters such as the background set, differential metabolite selection methods, and pathway database used can result in profoundly different ORA results. The use of a non-assay-specific background set, for example, resulted in large numbers of false-positive pathways. Pathway database choice, evaluated using three of the most popular metabolic pathway databases (KEGG, Reactome, and BioCyc), led to vastly different results in both the number and function of significantly enriched pathways. Factors that are specific to metabolomics data, such as the reliability of compound identification and the chemical bias of different analytical platforms also impacted ORA results. Simulated metabolite misidentification rates as low as 4% resulted in both gain of false-positive pathways and loss of truly significant pathways across all datasets. Our results have several practical implications for ORA users, as well as those using alternative pathway analysis methods. We offer a set of recommendations for the use of ORA in metabolomics, alongside a set of minimal reporting guidelines, as a first step towards the standardisation of pathway analysis in metabolomics.

Calm on the surface, dynamic on the inside. Molecular homeostasis of Anabaena sp. PCC 7120 nitrogen metabolism.

Nitrogen sources are all converted into ammonium/ia as a first step of assimilation. It is reasonable to expect that molecular components involved in the transport of ammonium/ia across biological membranes connect with the regulation of both nitrogen and central metabolism. We applied both genetic (i.e., Δamt mutation) and environmental treatments to a target biological system, the cyanobacterium Anabaena sp PCC 7120. The aim was to both perturb nitrogen metabolism and induce multiple inner nitrogen states, respectively, followed by targeted quantification of key proteins, metabolites and enzyme activities. The absence of AMT transporters triggered a substantial whole-system response, affecting enzyme activities and quantity of proteins and metabolites, spanning nitrogen and carbon metabolisms. Moreover, the Δamt strain displayed a molecular fingerprint indicating nitrogen deficiency even under nitrogen replete conditions. Contrasting with such dynamic adaptations was the striking near-complete lack of an externally measurable altered phenotype. We conclude that this species evolved a highly robust and adaptable molecular network to maintain homeostasis, resulting in substantial internal but minimal external perturbations. This analysis provides evidence for a potential role of AMT transporters in the regulatory/signalling network of nitrogen metabolism and the existence of a novel fourth regulatory mechanism controlling glutamine synthetase activity.

Remodelling of microRNAs in colorectal cancer by hypoxia alters metabolism profiles and 5-fluorouracil resistance.

Solid tumours have oxygen gradients and areas of near and almost total anoxia. Hypoxia reduces sensitivity to 5-fluorouracil (5-FU)-chemotherapy for colorectal cancer (CRC). MicroRNAs (miRNAs) are hypoxia sensors and were altered consistently in six CRC cell lines (colon cancer: DLD-1, HCT116 and HT29; rectal cancer: HT55, SW837 and VACO4S) maintained in hypoxia (1 and 0.2% oxygen) compared with normoxia (20.9%). CRC cell lines also showed altered amino acid metabolism in hypoxia and hypoxia-responsive miRNAs were predicted to target genes in four metabolism pathways: beta-alanine; valine, leucine, iso-leucine; aminoacyl-tRNA; and alanine, aspartate, glutamate. MiR-210 was increased in hypoxic areas of CRC tissues and hypoxia-responsive miR-21 and miR-30d, but not miR-210, were significantly increased in 5-FU resistant CRCs. Treatment with miR-21 and miR-30d antagonists sensitized hypoxic CRC cells to 5-FU. Our data highlight the complexity and tumour heterogeneity caused by hypoxia. MiR-210 as a hypoxic biomarker, and the targeting of miR-21 and miR-30d and/or the amino acid metabolism pathways may offer translational opportunities.

Identification of hydroxyapatite spherules provides new insight into subretinal pigment epithelial deposit formation in the aging eye.

Accumulation of protein- and lipid-containing deposits external to the retinal pigment epithelium (RPE) is common in the aging eye, and has long been viewed as the hallmark of age-related macular degeneration (AMD). The cause for the accumulation and retention of molecules in the sub-RPE space, however, remains an enigma. Here, we present fluorescence microscopy and X-ray diffraction evidence for the formation of small (0.5-20 µm in diameter), hollow, hydroxyapatite (HAP) spherules in Bruch's membrane in human eyes. These spherules are distinct in form, placement, and staining from the well-known calcification of the elastin layer of the aging Bruch's membrane. Secondary ion mass spectrometry (SIMS) imaging confirmed the presence of calcium phosphate in the spherules and identified cholesterol enrichment in their core. Using HAP-selective fluorescent dyes, we show that all types of sub-RPE deposits in the macula, as well as in the periphery, contain numerous HAP spherules. Immunohistochemical labeling for proteins characteristic of sub-RPE deposits, such as complement factor H, vitronectin, and amyloid beta, revealed that HAP spherules were coated with these proteins. HAP spherules were also found outside the sub-RPE deposits, ready to bind proteins at the RPE/choroid interface. Based on these results, we propose a novel mechanism for the growth, and possibly even the formation, of sub-RPE deposits, namely, that the deposit growth and formation begin with the deposition of insoluble HAP shells around naturally occurring, cholesterol-containing extracellular lipid droplets at the RPE/choroid interface; proteins and lipids then attach to these shells, initiating or supporting the growth of sub-RPE deposits.

GlnK Facilitates the Dynamic Regulation of Bacterial Nitrogen Assimilation.

Ammonium assimilation in Escherichia coli is regulated by two paralogous proteins (GlnB and GlnK), which orchestrate interactions with regulators of gene expression, transport proteins, and metabolic pathways. Yet how they conjointly modulate the activity of glutamine synthetase, the key enzyme for nitrogen assimilation, is poorly understood. We combine experiments and theory to study the dynamic roles of GlnB and GlnK during nitrogen starvation and upshift. We measure time-resolved in vivo concentrations of metabolites, total and posttranslationally modified proteins, and develop a concise biochemical model of GlnB and GlnK that incorporates competition for active and allosteric sites, as well as functional sequestration of GlnK. The model predicts the responses of glutamine synthetase, GlnB, and GlnK under time-varying external ammonium level in the wild-type and two genetic knock-outs. Our results show that GlnK is tightly regulated under nitrogen-rich conditions, yet it is expressed during ammonium run-out and starvation. This suggests a role for GlnK as a buffer of nitrogen shock after starvation, and provides a further functional link between nitrogen and carbon metabolisms.

Rapid screening of cellular stress responses in recombinant Pichia pastoris strains using metabolite profiling.

Heterologous protein production in the yeast Pichia pastoris can be limited by biological responses to high expression levels; the unfolded protein response (UPR) is a key determinant of the success of protein production in this organism. Here, we used untargeted NMR metabolic profiling (metabolomics) of a number of different recombinant strains, carried out in a miniaturized format suitable for screening-level experiments. We identified a number of metabolites (from both cell extracts and supernatants) which correlated well with UPR-relevant gene transcripts, and so could be potential biomarkers for future high-throughput screening of large numbers of P. pastoris clones.

Identification of the Elusive Pyruvate Reductase of Chlamydomonas reinhardtii Chloroplasts.

Under anoxic conditions the green alga Chlamydomonas reinhardtii activates various fermentation pathways leading to the creation of formate, acetate, ethanol and small amounts of other metabolites including d-lactate and hydrogen. Progress has been made in identifying the enzymes involved in these pathways and their subcellular locations; however, the identity of the enzyme involved in reducing pyruvate to d-lactate has remained unclear. Based on sequence comparisons, enzyme activity measurements, X-ray crystallography, biochemical fractionation and analysis of knock-down mutants, we conclude that pyruvate reduction in the chloroplast is catalyzed by a tetrameric NAD(+)-dependent d-lactate dehydrogenase encoded by Cre07.g324550. Its expression during aerobic growth supports a possible function as a 'lactate valve' for the export of lactate to the mitochondrion for oxidation by cytochrome-dependent d-lactate dehydrogenases and by glycolate dehydrogenase. We also present a revised spatial model of fermentation based on our immunochemical detection of the likely pyruvate decarboxylase, PDC3, in the cytoplasm.

Influence of the Crc regulator on the hierarchical use of carbon sources from a complete medium in Pseudomonas.

The Crc protein, together with the Hfq protein, participates in catabolite repression in pseudomonads, helping to coordinate metabolism. Little is known about how Crc affects the hierarchy of metabolite assimilation from complex mixtures. Using proton Nuclear Magnetic Resonance (NMR) spectroscopy, we carried out comprehensive metabolite profiling of culture supernatants (metabolic footprinting) over the course of growth of both Pseudomonas putida and P. aeruginosa, and compared the wild-type strains with deletion mutants for crc. A complex metabolite consumption hierarchy was observed, which was broadly similar between the two species, although with some important differences, for example in sugar utilization. The order of metabolite utilization changed upon inactivation of the crc gene, but even in the Crc-null strains some compounds were completely consumed before late metabolites were taken up. This suggests the presence of additional regulatory elements that determine the time and order of consumption of compounds. Unexpectedly, the loss of Crc led both species to excrete acetate and pyruvate as a result of unbalanced growth during exponential phase, compounds that were later consumed in stationary phase. This loss of carbon during growth helps to explain the contribution of the Crc/Hfq regulatory system to evolutionary fitness of pseudomonads.

Statistical Correlations between NMR Spectroscopy and Direct Infusion FT-ICR Mass Spectrometry Aid Annotation of Unknowns in Metabolomics.

NMR spectroscopy and mass spectrometry are the two major analytical platforms for metabolomics, and both generate substantial data with hundreds to thousands of observed peaks for a single sample. Many of these are unknown, and peak assignment is generally complex and time-consuming. Statistical correlations between data types have proven useful in expediting this process, for example, in prioritizing candidate assignments. However, this approach has not been formally assessed for the comparison of direct-infusion mass spectrometry (DIMS) and NMR data. Here, we present a systematic analysis of a sample set (tissue extracts), and the utility of a simple correlation threshold to aid metabolite identification. The correlations were surprisingly successful in linking structurally related signals, with 15 of 26 NMR-detectable metabolites having their highest correlation to a cognate MS ion. However, we found that the distribution of the correlations was highly dependent on the nature of the MS ion, such as the adduct type. This approach should help to alleviate this important bottleneck where both 1D NMR and DIMS data sets have been collected.

Metabolic Youth in Middle Age: Predicting Aging in Caenorhabditis elegans Using Metabolomics.

Many mutations and allelic variants are known that influence the rate at which animals age, but when in life do such variants diverge from normal patterns of aging? Is this divergence visible in their physiologies? To investigate these questions, we have used (1)H NMR spectroscopy to study how the metabolome of the nematode Caenorhabditis elegans changes as it grows older. We identify a series of metabolic changes that, collectively, predict the age of wild-type worms. We then show that long-lived mutant daf-2(m41) worms are metabolically youthful compared to wild-type worms, but that this relative youth only appears in middle age. Finally, we show that metabolic age predicts the timing and magnitude of differences in age-specific mortality between these strains. Thus, the future mortality of these two genotypes can be predicted long before most of the worms die.

Species interactions alter evolutionary responses to a novel environment.

Studies of evolutionary responses to novel environments typically consider single species or perhaps pairs of interacting species. However, all organisms co-occur with many other species, resulting in evolutionary dynamics that might not match those predicted using single species approaches. Recent theories predict that species interactions in diverse systems can influence how component species evolve in response to environmental change. In turn, evolution might have consequences for ecosystem functioning. We used experimental communities of five bacterial species to show that species interactions have a major impact on adaptation to a novel environment in the laboratory. Species in communities diverged in their use of resources compared with the same species in monocultures and evolved to use waste products generated by other species. This generally led to a trade-off between adaptation to the abiotic and biotic components of the environment, such that species evolving in communities had lower growth rates when assayed in the absence of other species. Based on growth assays and on nuclear magnetic resonance (NMR) spectroscopy of resource use, all species evolved more in communities than they did in monocultures. The evolutionary changes had significant repercussions for the functioning of these experimental ecosystems: communities reassembled from isolates that had evolved in polyculture were more productive than those reassembled from isolates that had evolved in monoculture. Our results show that the way in which species adapt to new environments depends critically on the biotic environment of co-occurring species. Moreover, predicting how functioning of complex ecosystems will respond to an environmental change requires knowing how species interactions will evolve.

Metabolite profiling to characterize disease-related bacteria: gluconate excretion by Pseudomonas aeruginosa mutants and clinical isolates from cystic fibrosis patients.

Metabolic footprinting of supernatants has been proposed as a tool for assigning gene function. We used NMR spectroscopy to measure the exometabolome of 86 single-gene transposon insertion mutant strains (mutants from central carbon metabolism and regulatory mutants) of the opportunistic pathogen Pseudomonas aeruginosa, grown on a medium designed to represent the nutritional content of cystic fibrosis sputum. Functionally related genes had similar metabolic profiles. E.g. for two-component system mutants, the cognate response regulator and sensor kinase genes clustered tightly together. Some strains had metabolic phenotypes (metabotypes) that could be related to the known gene function. E.g. pyruvate dehydrogenase mutants accumulated large amounts of pyruvate in the medium. In other cases, the metabolic phenotypes were not easily interpretable. The rpoN mutant, which lacks the alternative σ factor RpoN (σ(54)), accumulated high levels of gluconate in the medium. In addition, endometabolome profiling of intracellular metabolites identified a number of systemic metabolic changes. We linked this to indirect regulation of the catabolite repression protein Crc via the non-coding RNA crcZ and found that a crcZ (but not crc) mutant also shared the high-gluconate phenotype. We profiled an additional set of relevant metabolic enzymes and transporters, including Crc targets, and showed that the Crc-regulated edd mutant (gluconate-6-phosphate dehydratase) had similar gluconate levels as the rpoN mutant. Finally, a set of clinical isolates showed patient- and random amplification of polymorphic DNA (RAPD) type-specific differences in gluconate production, which were associated significantly with resistance across four antibiotics (tobramycin, ciprofloxacin, aztreonam, and imipenem), indicating that this has potential clinical relevance.

Identifying biochemical phenotypic differences between cryptic species.

Molecular genetic methods can distinguish divergent evolutionary lineages in what previously appeared to be single species, but it is not always clear what functional differences exist between such cryptic species. We used a metabolomic approach to profile biochemical phenotype (metabotype) differences between two putative cryptic species of the earthworm Lumbricus rubellus. There were no straightforward metabolite biomarkers of lineage, i.e. no metabolites that were always at higher concentration in one lineage. Multivariate methods, however, identified a small number of metabolites that together helped distinguish the lineages, including uncommon metabolites such as Nε-trimethyllysine, which is not usually found at high concentrations. This approach could be useful for characterizing functional trait differences, especially as it is applicable to essentially any species group, irrespective of its genome sequencing status.

Combining spectral ordering with peak fitting for one-dimensional NMR quantitative metabolomics.

One-dimensional (1)H NMR spectra are widely used for metabolic profiling. Such data sets often contain hundreds or thousands of spectra, which typically have variation in their sample chemistry, which leads to chemical shift variation "positional noise". This is a severe problem for metabolite quantification and data analysis, as peak integrals do not necessarily correspond across all spectra in a set. Various alignment algorithms have been developed to address this problem, but different studies have taken different approaches to evaluating the performance of NMR alignment routines and can be subjective or rely on an arbitrary cutoff. Furthermore, most alignment approaches completely fail to deal with peaks that overlap. We adopt the simple and robust method of ordering spectra with respect to an internally varying peak and use this to compare different alignment algorithms. Furthermore, we use the information from this procedure to help improve a Bayesian approach to automated peak deconvolution by restricting the prior probability distribution of the peak position in a model-free manner and compare the performance to manual peak deconvolution and to binning. This combination of spectral ordering and compound deconvolution improved the quality of the data for quantitative metabolomics.

Biochemical diversity of betaines in earthworms.

The ability to accumulate osmoprotectant compounds, such as betaines, is an important evolutionary feature in many organisms. This is particularly the case for organisms that live in variable environments, which may have fluctuations in moisture and salinity levels. There is, surprisingly, very little known about betaines in soil invertebrates in general, and there is almost no information about earthworms - a group that are important 'ecosystem engineers' and key indicators of soil health. Here, we describe a fast and reliable (1)H-(13)C heteronuclear single quantum coherence (HSQC) 2D NMR approach for the metabolic profiling of a series of betaines and related metabolites in tissue extracts, and list (1)H and (13)C chemical shifts for the trimethylammonium signal for 23 such compounds. The analysis of ten different species from three different families (Lumbricidae, Megascolecidae and Glossoscolecidae) showed an unexpected diversity of betaines present in earthworms. In total ten betaines were identified, including hydroxyproline-betaine, proline-betaine, taurine-betaine, GABA-betaine and histidine-betaine, and a further eleven as-yet unassigned putative betaine metabolites detected. The findings clearly indicate a hitherto-unappreciated important role for betaine metabolism in earthworms.

Direct assessment of metabolite utilization by Pseudomonas aeruginosa during growth on artificial sputum medium.

We grew Pseudomonas aeruginosa in LB and artificial sputum medium (ASM) (filtered and unfiltered) and quantified metabolite utilization and excretion by nuclear magnetic resonance (NMR) spectroscopy (metabolic footprinting or extracellular metabolomics). Utilization rates were similar between media, but there were differences in excretion-e.g., acetate was produced only in unfiltered ASM.

Epigenetic and metabolic reprogramming of fibroblasts in Crohn's disease strictures reveals histone deacetylases as therapeutic targets.

BACKGROUND & AIMS: No effective therapeutic intervention exists for intestinal fibrosis in Crohn's disease [CD]. We characterised fibroblast subtypes, epigenetic and metabolic changes, and signalling pathways in CD fibrosis to inform future therapeutic strategies. METHODS: We undertook immunohistochemistry, metabolic, signalling pathway and Epigenetic [Transposase-Accessible Chromatin using sequencing] analyses associated with collagen production in CCD-18Co intestinal fibroblasts and primary fibroblasts isolated from stricturing [SCD] and non-stricturing [NSCD] CD small intestine. SCD/ NSCD fibroblasts were cultured with TGFß and valproic acid [VPA]. RESULTS: Stricturing CD was characterised by distinct histone deacetylase [HDAC] expression profiles, particularly HDAC1, HDAC2, and HDAC7. As a proxy for HDAC activity, reduced numbers of H3K27ac+ cells were found in SCD compared to NSCD sections. Primary fibroblasts had increased extracellular lactate [increased glycolytic activity] and intracellular hydroxyproline [increased collagen production] in SCD compared to NSCD cultures. The metabolic effect of TGFß-stimulation was reversed by the HDAC inhibitor VPA. SCD fibroblasts appear "metabolically primed" and responded more strongly to both TGFß and VPA. Treatment with VPA revealed TGFß-dependent and independent Collagen-I production in CCD-18Co cells and primary fibroblasts. VPA altered the epigenetic landscape with reduced chromatin accessibility at the COL1A1 and COL1A2 promoters. CONCLUSIONS: Increased HDAC expression profiles, H3K27ac hypoacetylation, a significant glycolytic phenotype, and metabolic priming, characterise SCD-derived as compared to NSCD fibroblasts. Our results reveal a novel epigenetic component to Collagen-I regulation and TGFß-mediated CD fibrosis. HDAC inhibitor therapy may 'reset' the epigenetic changes associated with fibrosis.

Free glucosylglycerate is a novel marker of nitrogen stress in Mycobacterium smegmatis.

Nitrogen is an essential element for bacterial growth, and as such, bacteria have evolved several pathways to assimilate nitrogen and adapt to situations of nitrogen limitation. However, the adaptation of mycobacteria to nitrogen stress and the regulation of the stress response pathways is unknown. Identification of key metabolites produced by mycobacteria during nitrogen stress could therefore provide important insights into mycobacterial survival strategies. Here we used NMR-based metabolomics to monitor and quantify intracellular and extracellular metabolite levels (metabolic footprinting) in Mycobacterium smegmatis grown under nitrogen-limiting and nitrogen-rich conditions. There were several metabolic differences between the two conditions: following nitrogen run-out, there was an increase in intracellular α-ketoglutarate and a decrease in intracellular glutamine and glutamate levels. In addition, a sugar-derived compound accumulated in nitrogen-starved cells that was subsequently assigned as glucosylglycerate (GGA). Free GGA production was responsive to nitrogen stress in M. smegmatis but not to oxidative or osmotic stress; lack of a functional GGA synthesis pathway slightly reduced growth and decreased ammonium uptake rates under nitrogen-limiting conditions. Hence, GGA could contribute to the fitness of mycobacteria under nitrogen limitation.