MetabR: an R script for linear model analysis of quantitative metabolomic data

Background Metabolomics is an emerging high-throughput approach to systems biology, but data analysis tools are lacking compared to other systems level disciplines such as transcriptomics and proteomics. Metabolomic data analysis requires a normalization step to remove systematic effects of confounding variables on metabolite measurements. Current tools may not correctly normalize every metabolite when the relationships between each metabolite quantity and fixed-effect confounding variables are different, or for the effects of random-effect confounding variables. Linear mixed models, an established methodology in the microarray literature, offer a standardized and flexible approach for removing the effects of fixed- and random-effect confounding variables from metabolomic data. Findings Here we present a simple menu-driven program, “MetabR”, designed to aid researchers with no programming background in statistical analysis of metabolomic data. Written in the open-source statistical programming language R, MetabR implements linear mixed models to normalize metabolomic data and analysis of variance (ANOVA) to test treatment differences. MetabR exports normalized data, checks statistical model assumptions, identifies differentially abundant metabolites, and produces output files to help with data interpretation. Example data are provided to illustrate normalization for common confounding variables and to demonstrate the utility of the MetabR program. Conclusions We developed MetabR as a simple and user-friendly tool for implementing linear mixed model-based normalization and statistical analysis of targeted metabolomic data, which helps to fill a lack of available data analysis tools in this field. The program, user guide, example data, and any future news or updates related to the program may be found at http://metabr.r-forge.r-project.org

MetabR: an R script for linear model analysis of quantitative metabolomic data

Background Metabolomics is an emerging high-throughput approach to systems biology, but data analysis tools are lacking compared to other systems level disciplines such as transcriptomics and proteomics. Metabolomic data analysis requires a normalization step to remove systematic effects of confounding variables on metabolite measurements. Current tools may not correctly normalize every metabolite when the relationships between each metabolite quantity and fixed-effect confounding variables are different, or for the effects of random-effect confounding variables. Linear mixed models, an established methodology in the microarray literature, offer a standardized and flexible approach for removing the effects of fixed- and random-effect confounding variables from metabolomic data. Findings Here we present a simple menu-driven program, “MetabR”, designed to aid researchers with no programming background in statistical analysis of metabolomic data. Written in the open-source statistical programming language R, MetabR implements linear mixed models to normalize metabolomic data and analysis of variance (ANOVA) to test treatment differences. MetabR exports normalized data, checks statistical model assumptions, identifies differentially abundant metabolites, and produces output files to help with data interpretation. Example data are provided to illustrate normalization for common confounding variables and to demonstrate the utility of the MetabR program. Conclusions We developed MetabR as a simple and user-friendly tool for implementing linear mixed model-based normalization and statistical analysis of targeted metabolomic data, which helps to fill a lack of available data analysis tools in this field. The program, user guide, example data, and any future news or updates related to the program may be found at http://metabr.r-forge.r-project.org/ webcite

Metabolomics reveals distinct neurochemical profiles associated with stress resilience

Acute social defeat represents a naturalistic form of conditioned fear and is an excellent model in which to investigate the biological basis of stress resilience. While there is growing interest in identifying biomarkers of stress resilience, until recently, it has not been feasible to associate levels of large numbers of neurochemicals and metabolites to stress-related phenotypes. The objective of the present study was to use an untargeted metabolomics approach to identify known and unknown neurochemicals in select brain regions that distinguish susceptible and resistant individuals in two rodent models of acute social defeat. In the first experiment, male mice were first phenotyped as resistant or susceptible. Then, mice were subjected to acute social defeat, and tissues were immediately collected from the ventromedial prefrontal cortex (vmPFC), basolateral/central amygdala (BLA/CeA), nucleus accumbens (NAc), and dorsal hippocampus (dHPC). Ultra-high performance liquid chromatography coupled with high resolution mass spectrometry (UPLC-HRMS) was used for the detection of water-soluble neurochemicals. In the second experiment, male Syrian hamsters were paired in daily agonistic encounters for 2 weeks, during which they formed stable dominant-subordinate relationships. Then, 24 h after the last dominance encounter, animals were exposed to acute social defeat stress. Immediately after social defeat, tissue was collected from the vmPFC, BLA/CeA, NAc, and dHPC for analysis using UPLC-HRMS. Although no single biomarker characterized stress-related phenotypes in both species, commonalities were found. For instance, in both model systems, animals resistant to social defeat stress also show increased concentration of molecules to protect against oxidative stress in the NAc and vmPFC. Additionally, in both mice and hamsters, unidentified spectral features were preliminarily annotated as potential targets for future experiments. Overall, these findings suggest that a metabolomics approach can identify functional groups of neurochemicals that may serve as novel targets for the diagnosis, treatment, or prevention of stress-related mental illness

Bone Metabolite Profile Differs between Normal and Femur Head Necrosis (FHN/BCO)-Affected Broilers: Implications for Dysregulated Metabolic Cascades in FHN Pathophysiology

Femur head necrosis (FHN), also known as bacterial chondronecrosis with osteomyelitis (BCO), has remained an animal welfare and production concern for modern broilers regardless of efforts to select against it in primary breeder flocks. Characterized by the bacterial infection of weak bone, FHN has been found in birds without clinical lameness and remains only detectable via necropsy. This presents an opportunity to utilize untargeted metabolomics to elucidate potential non-invasive biomarkers and key causative pathways involved in FHN pathology. The current study used ultra-performance liquid chromatography coupled with high-resolution mass spectrometry (UPLC–HRMS) and identified a total of 152 metabolites. Mean intensity differences at p \u3c 0.05 were found in 44 metabolites, with 3 significantly down-regulated and 41 up-regulated in FHN-affected bone. Multivariate analysis and a partial least squares discriminant analysis (PLS-DA) scores plot showed the distinct clustering of metabolite profiles from FHN-affected vs. normal bone. Biologically related molecular networks were predicted using an ingenuity pathway analysis (IPA) knowledge base. Using a fold-change cut off of −1.5 and 1.5, top canonical pathways, networks, diseases, molecular functions, and upstream regulators were generated using the 44 differentially abundant metabolites. The results showed the metabolites NAD+, NADP+, and NADH to be downregulated, while 5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR) and histamine were significantly increased in FHN. Ascorbate recycling and purine nucleotides degradation were the top canonical pathways, indicating the potential dysregulation of redox homeostasis and osteogenesis. Lipid metabolism and cellular growth and proliferation were some of the top molecular functions predicted based on the metabolite profile in FHN-affected bone. Network analysis showed significant overlap across metabolites and predicted upstream and downstream complexes, including AMP-activated protein kinase (AMPK), insulin, collagen type IV, mitochondrial complex, c-Jun N-terminal kinase (Jnk), extracellular signal-regulated kinase (ERK), and 3β-hydroxysteroid dehydrogenase (3β HSD). The qPCR analysis of relevant factors showed a significant decrease in AMPKα2 mRNA expression in FHN-affected bone, supporting the predicted downregulation found in the IPA network analysis. Taken as a whole, these results demonstrate a shift in energy production, bone homeostasis, and bone cell differentiation that is distinct in FHN-affected bone, with implications for how metabolites drive the pathology of FHN

Evolution of the Toxins Muscarine and Psilocybin in a Family of Mushroom-Forming Fungi

Mushroom-forming fungi produce a wide array of toxic alkaloids. However, evolutionary analyses aimed at exploring the evolution of muscarine, a toxin that stimulates the parasympathetic nervous system, and psilocybin, a hallucinogen, have never been performed. The known taxonomic distribution of muscarine within the Inocybaceae is limited, based only on assays of species from temperate regions of the northern hemisphere. Here, we present a review of muscarine and psilocybin assays performed on species of Inocybaceae during the last fifty years. To supplement these results, we used liquid chromatography–tandem mass spectrometry (LC–MS/MS) to determine whether muscarine was present in 30 new samples of Inocybaceae, the majority of which have not been previously assayed or that originated from either the tropics or temperate regions of the southern hemisphere. Our main objective is to test the hypothesis that the presence of muscarine is a shared ancestral feature of the Inocybaceae. In addition, we also test whether species of Inocyabceae that produce psilocybin are monophyletic. Our findings suggest otherwise. Muscarine has evolved independently on several occasions, together with several losses. We also detect at least two independent transitions of muscarine-free lineages to psilocybin-producing states. Although not ancestral for the family as a whole, muscarine is a shared derived trait for an inclusive clade containing three of the seven major lineages of Inocybaceae (the Inocybe, Nothocybe, and Pseudosperma clades), the common ancestor of which may have evolved ca. 60 million years ago. Thus, muscarine represents a conserved trait followed by several recent losses. Transitions to psilocybin from muscarine-producing ancestors occurred more recently between 10–20 million years ago after muscarine loss in two separate lineages. Statistical analyses firmly reject a single origin of muscarine-producing taxa. DOI: 10.1371/journal.pone.006464

Evolution of the Toxins Muscarine and Psilocybin in a Family of Mushroom-Forming Fungi

Mushroom-forming fungi produce a wide array of toxic alkaloids. However, evolutionary analyses aimed at exploring the evolution of muscarine, a toxin that stimulates the parasympathetic nervous system, and psilocybin, a hallucinogen, have never been performed. The known taxonomic distribution of muscarine within the Inocybaceae is limited, based only on assays of species from temperate regions of the northern hemisphere. Here, we present a review of muscarine and psilocybin assays performed on species of Inocybaceae during the last fifty years. To supplement these results, we used liquid chromatography–tandem mass spectrometry (LC–MS/MS) to determine whether muscarine was present in 30 new samples of Inocybaceae, the majority of which have not been previously assayed or that originated from either the tropics or temperate regions of the southern hemisphere. Our main objective is to test the hypothesis that the presence of muscarine is a shared ancestral feature of the Inocybaceae. In addition, we also test whether species of Inocyabceae that produce psilocybin are monophyletic. Our findings suggest otherwise. Muscarine has evolved independently on several occasions, together with several losses. We also detect at least two independent transitions of muscarine-free lineages to psilocybin-producing states. Although not ancestral for the family as a whole, muscarine is a shared derived trait for an inclusive clade containing three of the seven major lineages of Inocybaceae (the Inocybe, Nothocybe, and Pseudosperma clades), the common ancestor of which may have evolved ca. 60 million years ago. Thus, muscarine represents a conserved trait followed by several recent losses. Transitions to psilocybin from muscarine-producing ancestors occurred more recently between 10–20 million years ago after muscarine loss in two separate lineages. Statistical analyses firmly reject a single origin of muscarine-producing taxa. DOI: 10.1371/journal.pone.006464

Functional Characteristics of the Gut Microbiome in C57BL/6 Mice Differentially Susceptible to Plasmodium yoelii

C57BL/6 mice are widely used for in vivo studies of immune function and metabolism in mammals. In a previous study, it was observed that when C57BL/6 mice purchased from different vendors were infected with Plasmodium yoelii, a causative agent of murine malaria, they exhibited both differential immune responses and significantly different parasite burdens: these patterns were reproducible when gut contents were transplanted into gnotobiotic mice. To gain insight into the mechanism of resistance, we removed whole ceca from mice purchased from two vendors, Taconic Biosciences (low parasitemia) and Charles River Laboratories (high parasitemia), to determine the combined host and microflora metabolome and metatranscriptome. With the exception of two Charles River samples, we observed 90% similarity in overall bacterial gene expression within vendors and 80% similarity between vendors. In total 33 bacterial genes were differentially expressed in Charles River mice (p-value \u3c 0.05) relative to the mice purchased from Taconic. Included among these, fliC, ureABC, and six members of the nuo gene family were overrepresented in microbiomes susceptible to more severe malaria. Moreover, 38 mouse genes were differentially expressed in these purported genetically identical mice. Differentially expressed genes included basigin, a cell surface receptor required for P. falciparum invasion of red blood cells. Differences in metabolite pools were detected, though their relevance to malaria infection, microbial community activity, or host response is not yet understood. Our data have provided new targets that may connect gut microbial activity to malaria resistance and susceptibility phenotypes in the C57BL/6 model organism

Urea Is Both a Carbon and Nitrogen Source for Microcystis aeruginosa: Tracking 13C Incorporation at Bloom pH Conditions

The use of urea as a nitrogenous fertilizer has increased over the past two decades, with urea itself being readily detected at high concentrations in many lakes. Urea has been linked to cyanobacterial blooms as it is a readily assimilated nitrogen (N) - source for cyanobacteria that possess the enzyme urease. We tested the hypothesis that urea may also act as a carbon (C) source to supplemental growth requirements during the alkaline conditions created by dense cyanobacterial blooms, when concentrations of dissolved CO2 are vanishingly low. High rates of photosynthesis markedly reduce dissolved CO2 concentrations and drive up pH. This was observed in Lake Erie during the largest bloom on record (2015) over long periods (months) and short periods (days) of time, suggesting blooms experience periods of CO2-limitation on a seasonal and daily basis. We used 13C-urea to demonstrate that axenic cultures of the model toxic cyanobacterium, Microcystis aeruginosaNIES843, assimilated C at varying environmentally relevant pH conditions directly into a spectrum of metabolic pools during urea hydrolysis. Primarily, 13C from urea was assimilated into central C metabolism and amino acid biosynthesis pathways, including those important for the production of the hepatotoxin, microcystin, and incorporation into these pathways was at a higher percentage during growth at higher pH. This corresponded to increased growth rates on urea as the sole N source with increasing pH. We propose this ability to incorporate C from urea represents yet another competitive advantage for this cyanobacterium during dense algal blooms

Pancreatic β-Cell Death in Response to Pro-Inflammatory Cytokines Is Distinct from Genuine Apoptosis

A reduction in functional β-cell mass leads to both major forms of diabetes; pro-inflammatory cytokines, such as interleukin-1beta (IL-1β) and gamma-interferon (γ-IFN), activate signaling pathways that direct pancreatic β-cell death and dysfunction. However, the molecular mechanism of β-cell death in this context is not well understood. In this report, we tested the hypothesis that individual cellular death pathways display characteristic phenotypes that allow them to be distinguished by the precise biochemical and metabolic responses that occur during stimulus-specific initiation. Using 832/13 and INS-1E rat insulinoma cells and isolated rat islets, we provide evidence that apoptosis is unlikely to be the primary pathway underlying β-cell death in response to IL-1β+γ-IFN. This conclusion was reached via the experimental results of several different interdisciplinary strategies, which included: 1) tandem mass spectrometry to delineate the metabolic differences between IL-1β+γ-IFN exposure versus apoptotic induction by camptothecin and 2) pharmacological and molecular interference with either NF-κB activity or apoptosome formation. These approaches provided clear distinctions in cell death pathways initiated by pro-inflammatory cytokines and bona fide inducers of apoptosis. Collectively, the results reported herein demonstrate that pancreatic β-cells undergo apoptosis in response to camptothecin or staurosporine, but not pro-inflammatory cytokines. DOI: 10.1371/journal.pone.002248

Root-Associated Streptomyces Isolates Harboring melC Genes Demonstrate Enhanced Plant Colonization

Streptomycetaceae assemble into the internal, root endophytic compartment of a wide variety of plants grown in soils worldwide, suggesting their ability to survive during root microbiome assembly. A previous study found that among four nonpathogenic, root-isolated Streptomyces strains (303, 299, CL18, and 136), only 303 and 299 colonized endophytic root tissue of the majority of Arabidopsis thaliana roots when inoculated with 34 other bacterial isolates. Here we demonstrate that 303 and 299 also colonize significantly more in singly inoculated A. thaliana seedlings. The genomes of melanin-producing 303 and 299 each contain two copies of the gene encoding tyrosinase (melC2 and melD2), an enzyme necessary for melanin biosynthesis in Streptomyces. These genes were not found in the genomes of 136 or CL18. Tyrosinase activity was detected in 303 and 299 whole cell and supernatant protein extracts, suggesting functional intracellular and extracellular enzymes.. Because tyrosinase oxidizes phenolic compounds and Streptomyces colonization of A. thaliana appears to be influenced by the phenolic compound salicylic acid (SA), we measured direct sensitivity of Streptomyces isolates to the phenolic compounds catechol, ferulic acid (FA), and SA in vitro. While both 303 and 299 showed higher numbers of surviving colonies than CL18 and 136 in the presence of catechol, only 303 demonstrated a higher number of surviving colonies when isolates were challenges with FA and SA. Finally, when seedlings were singly inoculated with a collection of related plant-associated Streptomyces isolates, colonization was significantly higher in isolates possessing two tyrosinase gene copies than isolates with zero or one gene copy. Overall, we describe a connection between microbial tyrosinase activity and increased seedling colonization of nonpathogenic Streptomyces isolates in A. thaliana. We propose tyrosinase activity in Streptomyces partially protects against harmful plant-produced phenolic compounds as they transition into an endophytic lifestyle