Abstract
Polybrominated diphenyl ethers (PBDEs) are persistent environmental toxicants associated with increased risk for metabolic syndrome. Intermediary metabolism is influenced by the intestinal microbiome. To test the hypothesis that PBDEs reduce host-beneficial intermediary metabolites in an intestinal microbiome–dependent manner, 9-week old male conventional (CV) and germ-free (GF) C57BL/6 mice were orally gavaged once daily with vehicle, BDE-47, or BDE-99 (100 μmol/kg) for 4 days. Intestinal microbiome (16S rDNA sequencing), liver transcriptome (RNA-Seq), and intermediary metabolites in serum, liver, as well as small and large intestinal contents (SIC and LIC; LC-MS) were examined. Changes in intermediary metabolite abundances in serum, liver, and SIC, were observed under basal conditions (CV vs. GF mice) and by PBDE exposure. PBDEs altered the largest number of metabolites in the LIC; most were regulated by PBDEs in GF conditions. Importantly, intestinal microbiome was necessary for PBDE-mediated decreases in branched-chain and aromatic amino acid metabolites, including 3-indolepropionic acid, a tryptophan metabolite recently shown to be protective against inflammation and diabetes. Gene-metabolite networks revealed a positive association between the hepatic glycan synthesis gene α-1,6-mannosyltransferase (Alg12) mRNA and mannose, which are important for protein glycosylation. Glycome changes have been observed in patients with metabolic syndrome. In LIC of CV mice, 23 bacterial taxa were regulated by PBDEs. Correlations of certain taxa with distinct serum metabolites further highlight a modulatory role of the microbiome in mediating PBDE effects. In summary, PBDEs impact intermediary metabolism in an intestinal microbiome–dependent manner, suggesting that dysbiosis may contribute to PBDE-mediated toxicities that include metabolic syndrome.
Introduction
Polybrominated diphenyl ethers (PBDEs) have been added to a myriad of consumer products as flame retardants (2001). Although there are 209 predicted PBDE congeners stemming from the degree of bromination and stereochemistry, three commercial formulations were manufactured and used in the United States (2017). One mixture, known as commercial pentaBDE and sold under several trade names, was primarily applied to furniture cushion foam and carpet pads. OctaBDE and decaBDE mixtures, also sold under different trade names, were primarily applied to plastic components of computers, televisions, and other electronics (2017).
Despite the fact that commercial formulations are no longer produced in, or imported into, the United States, PBDEs are ubiquitous in human serum and breast milk samples. Recent data from the Methods Advancement for Milk Analysis (MAMA) Study showed that the median serum concentration of BDE-47 was 18.6 ng/g lipid and 3.9 ng/g lipid for BDE-99 (Marchitti et al., 2017). Median breast milk concentrations for BDE-47 were observed in the same study to be 31.5 ng/g lipid, and 6 ng/g lipid for BDE-99. One route of exposure is from household and workplace dust that is contaminated from PBDE-containing products ((Hites, 2004, US-EPA, 2012, 2017)). In addition, PBDEs are environmentally persistent and bioaccumulative, facilitating human exposures through diet, especially dairy, fish, and meat ((Schecter et al., 2004, 2010, 2017)). Two commercial PBDE formulations (e.g., DE-71 and Bromkal 70-5DE) that were analyzed for individual congeners consisted of 36%–42% BDE-47 (tetrabrominated), as well as 42%–48% BDE-99 (pentabrominated); thus, BDE-47 and BDE-99 account for approx. 80%–90% of total PBDEs in these mixtures (Sanders et al., 2005; LaA Guardia et al., 2006). BDE-47 and BDE-99 are also two of the most prevalent PBDE congeners found in human samples (Sjödin et al., 2001; Schecter et al., 2003; Imm et al., 2009; Marchitti et al., 2017).
Furthermore, recent epidemiologic and animal studies have associations between PBDEs and signs of diabetes and metabolic syndrome. In two independent case-control studies, serum concentrations of BDE-47 in participants with diabetes were significantly associated with increased risk of diabetes (Zhang et al., 2016). In a study among obese individuals, adipose tissue concentrations of both BDE-47 and BDE-99 were significantly elevated in insulin-resistant participants compared with insulin-sensitive participants (Helaleh et al., 2018). BDE-47 and BDE-99 concentrations were also associated with adipose tissue concentrations of insulin and interleukin-6. BDE-47 was associated with alkaline phosphatase, which is a biomarker for dysfunction in hepatobiliary secretion, and has also been shown to interfere with rodent glucose homeostasis (McIntyre et al., 2015; Zhang et al., 2016). Rats exposed to DE-71, one of the commercial PBDE mixtures made up primarily of BDE-47 and BDE-99, as well as other minor constituents exhibited changes in liver gene expression patterns that matched patterns associated with metabolic syndrome (Sanders et al., 2005; Dunnick et al., 2012). However, research into mechanisms for PBDE-mediated changes in intermediary metabolism is ongoing.
The dynamic and diverse populations of microbes that populate the lower gastrointestinal tract, collectively known as the intestinal or gut microbiome, are known to interact with host intermediary and xenobiotic metabolism (Nicholson et al., 2012; Li et al., 2017). We have previously shown that the intestinal microbiome modifies PBDE-mediated effects on bile acids as well as the expression of genes involved in xenobiotic biotransformation (Li et al., 2017, 2018). However, little is known regarding PBDE-mediated intermediary metabolic changes in rodents, and a potential modulatory role for PBDE in the intestinal microbiome. Thus, in this study we integrated transcriptomics, targeted metabolomics, and microbial 16S rDNA profiling to evaluate the effects of PBDEs on intermediary metabolism, and we tested our hypothesis that PBDEs downregulate beneficial intermediary metabolites in an intestinal microbiome–dependent manner.
Materials and Methods
Chemicals.
BDE-47 (cat. no. N-105220-10mg, CAS no. 5346-43-1) and BDE-99 (cat. no. FF-BDE-099N-80mg, CAS no. 60348-60-9) were purchased from Chem Service, Inc. (West Chester, PA), and Accustandard, Inc. (New Haven, CT), respectively. Phosphate-buffered saline (pH 7.4, cat. no. 10010023) was purchased from Thermo Fisher Scientific (Waltham, MA). Deionized (DI) water was provided in-house by a Synergy Ultrapure Water System from EMD Millipore (Billerica, MA). The unlabeled metabolite standards were purchased from Sigma-Aldrich (St. Louis, MO) and Fisher Scientific (Pittsburgh, PA) and labeled standards were purchased from Cambridge Isotopes (Tewksbury, MA) (Zhu et al., 2014). Ammonium acetate [cat. no. A11450, CAS no. 631-61-8; liquid chromatography–mass spectrometry (LC-MS) grade] and acetic acid (cat. no. A11350, CAS no. 64-19-7; LC-MS grade) came from Fisher Scientific. Acetonitrile (ACN; cat. no. 1.00029, CAS no. 75-05-8) and methanol (cat. no. 1.06035, CAS no. 67-56-1) both LC-MS grade, ammonium bicarbonate (cat. no. 09830, CAS no. 1066-33-7; purity >99.99%), and all other chemicals and reagents, unless indicated otherwise, were purchased from Sigma Aldrich.
Animals.
As previously reported (Li et al., 2017, 2018) and briefly summarized here, 8-week-old male C57BL/6J CV mice were purchased from the Jackson Laboratory (Bar Harbor, ME). Mice were acclimated at the University of Washington for 1 week prior to experiments. The initial germ-free (GF) mice, which are on a C57BL/6J background, were purchased from the National Gnotobiotic Rodent Resource Center (University of North Carolina, Chapel Hill, NC). Breeding and colony management of GF mice was performed by the University of Washington Gnotobiotic Animal Core. All mice were housed according to the Association for Assessment and Accreditation of Laboratory Animal Care International guidelines. At 9 weeks of age, the mice (n = 5 per group) were exposed to vehicle (corn oil, 10 ml/kg), BDE-47 [48.5 mg/kg (100 μmol/kg)], or BDE-99 [56.5 mg/kg (or 100 μmol/kg)] via oral gavage once daily for 4 consecutive days. All mice were euthanized using CO2 narcosis followed by cardiac puncture 24 hours after the final PBDE dose. During the cardiac puncture, blood was collected from the heart with a 25-gauge needle and put into a serum separator tube. Serum was collected and stored at −80°C after centrifugation for 10 minutes at 2000g at 4°C. Liver tissue was frozen in liquid nitrogen and stored at −80°C for further analysis. Small (SIC) and large (LIC) intestinal contents were collected as a pellet first by flushing the lumen with 10 mM dithiothreitol in PBS (Sigma Aldrich) and then the luminal contents were centrifuged at 20,000g for 30 minutes at 4°C. All animal studies were approved by the Institutional Animal Care and Use Committee at the University of Washington.
RNA Isolation and Sequencing.
RNA isolation and sequencing methods were previously reported (Li et al., 2017). Briefly, total RNA was isolated from frozen liver tissue (n = 2–4/group) using RNA-Bee (Tel-Test Inc., Friendswood, TX). Construction of cDNA libraries and sequencing was performed as previously reported (Li et al., 2017). Libraries were constructed with an Illumina TruSeq Stranded mRNA kit (Illumina, San Diego, CA) and integrity was quantified using an Agilent 2100 Bioanalyzer. Libraries were sequenced using 50-bp paired-end reads on an Illumina HiSEq. 2000.
16S rDNA Sequencing.
The sequencing of the LIC microbial DNA from CV mice (n = 3/group) was previously reported (Li et al., 2018). Briefly, DNA from the V4 region of the 16S rRNA gene was amplified and 250-bp paired-end sequencing was performed on a HiSEq. 2500 sequencer.
Serum Preparation for Aqueous Metabolite Measurements.
Frozen serum samples were first thawed overnight under 4°C, and 50 μl of each serum sample was placed in a 2-ml Eppendorf vial. The initial step for protein precipitation and metabolite extraction was performed by adding 250 μl of MeOH. The mixture was then vortexed for 10 seconds and stored at −20°C for 30 minutes, followed by centrifugation at 14,000 rpm for 10 minutes at 4°C. The supernatants (200 μl) were collected into a new Eppendorf vial and dried under vacuum using an Eppendorf Vacufuge (Eppendorf, Hauppauge, NY). The dried samples were reconstituted in 600 μl of 40% solvent A/60% solvent B (see details in the LC-MS method section).
Tissue Preparation for Aqueous Metabolite Quantifications.
The sample preparation is the same for liver, as well as small and large intestinal contents. Briefly, each tissue sample (approx. 10 mg) was homogenized in 200 μl water in an Eppendorf tube using a Bullet Blender homogenizer (Next Advance, Averill Park, NY). Then 800 μl of MeOH was added, and after vortexing for 10 second, the samples were stored on dry ice for 30 minutes. The samples were then sonicated in an ice bath for 10 minutes. The samples were centrifuged at 14,000 rpm for 10 minutes (4°C), and 800 μl of supernatant was transferred to a new Eppendorf tube. The samples were then dried under vacuum using an Eppendorf Vacufuge (Eppendorf). Prior to MS analysis, the residue obtained was reconstituted in 600 μl of 40% solvent A/60% solvent B (see details in the LC-MS methods section next).
Measurements of Aqueous Metabolites Using LC-MS.
A robust targeted liquid chromatography–tandem mass spectrometry (LC-MS/MS) method has been developed (Zhu et al., 2014) and used in a growing number of studies (Barton et al., 2015; Carroll et al., 2015; Gu et al., 2015, 2016; Reyes et al., 2015; Sood et al., 2015; Sperber et al., 2015; Zhu et al., 2015; Deng et al., 2016). Briefly, the LC-MS/MS experiments were performed on a Waters Acquity I-Class UPLC TQS-micro MS (Waters, Milford, MA) system. Each prepared serum or tissue sample was injected twice with 2 and 5 μl for analysis using positive and negative ionization mode, respectively. Chromatographic separation was performed on a Waters Xbridge BEH Amide column (2.5 μm, 2.1 × 150 mm) at 40°C. The flow rate was 0.3 ml/min. For positive mode, the mobile phase was composed of solvent A (5 mM ammonium acetate in H2O with 0.1% formic acid) and solvent B (ACN with 0.1% acetic acid). For negative mode, solvent A was 10 mM ammonium bicarbonate in H2O, and solvent B was ACN. The LC gradient conditions were the same for both positive and negative ionization modes. After an initial 1.5-minute isocratic elution of 10% A, the percentage of solvent A was increased to 65% at t = 9 minutes. The percentage of A was then kept the same (65%) for 5 minutes (t = 14 minute), and the percentage of A was decreased back to 10% at t = 15 minutes to prepare for the next injection. The total experimental time for each injection was 30 minutes. The metabolite identities were confirmed by spiking mixtures of standard compounds. The extracted peaks were integrated using the TargetLynx software (Waters). A list of all metabolites detected across biocompartments can be found in Supplemental Table 1. Supplemental Tables 2–5 show mean abundance values and S.E. for all metabolites by exposure group that were detectable in each biocompartment.
Data Analysis
RNA-Seq.
Preliminary analyses of the FASTQ files used in this study, which were deposited by our group in NCBI’s Gene Expression Omnibus (GEO) (GEO accession number GSE101650) during publication of a previous manuscript, have been previously reported (Li et al., 2017). Abundance was expressed as fragments per kilobase of transcript per million mapped reads. Differential expression analysis was performed using Cuffdiff (Trapnell et al., 2010). The differentially expressed genes were defined as false discovery rate (FDR) adjusted P-value < 0.05, in at least one of the PBDE-exposed groups compared with the vehicle-exposed group of the same enterotype, which in this study refers to CV as opposed to GF mice, as well as genes that were differentially expressed by the lack of intestinal microbiome under basal conditions. Assignment of genes into intermediary metabolism categories was performed using Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolism pathways in R 3.4.3 (Kanehisa and Goto, 2000; Kanehisa et al., 2016; Kanehisa et al., 2017; R Core Team, 2017). Mean FPKM S.E. and metabolism categories for intermediary metabolism genes for each treatment group can be found in Supplemental Table 6. Fold changes and associated FDR adjusted P values for comparisons between PBDE-treated mice and enterotype-matched vehicle controls and between CV and GF mice under basal conditions can be found in Supplemental Table 7. Heatmaps were also built in R 3.4.3 (base, gplots package) (Warnes et al., 2013; R Core Team, 2017). Because many gene products have multiple functions, many of the differentially expressed genes classified as intermediary metabolism genes here were reported as drug processing genes previously, mainly cytochrome P450s by Li et al. (2017).
Aqueous Metabolites.
Relative abundance values for each metabolite in each biocompartment were compared using two-way analysis of variance (ANOVA) followed by Tukey’s post-hoc test using R 3.4.3 (base, lsmeans, and multcomp packages) (Hothorn et al., 2008; Lenth, 2016; R Core Team, 2017). Heatmaps were also built in R (base, gplots package) (Warnes et al., 2013; R Core Team, 2017). Gene-metabolite network analyses were performed using MetaboAnalyst (Chong et al., 2018) (http://www.metaboanalyst.ca) using all differentially regulated intermediary metabolism genes with a fold-change >2 (implemented to emphasize most highly changed genes) and significantly altered metabolites in the liver. Pearson’s correlations for gene and metabolite network members were calculated using R 3.4.3 (base, corrplot, and Hmisc packages) (Harrell and Dupont, 2016; Wei and Simko, 2016, R Core Team, 2017). Significant gene-metabolite pair correlation plots and bar plots were created using R3 .4.3 (base, ggplot2 package) (Wickham, 2009; R Core Team, 2017).
16S rDNA Sequencing.
The 16S rDNA sequencing data were analyzed using QIIME v1.9.1 as previously reported (Li et al., 2018). Relative abundance was expressed as percentage of operational taxonomical units (% OTUs) and was examined for PBDE treatment effects using a one-way ANOVA followed by Tukey’s HSD post-hoc test in R 3.4.3 (P < 0.05) (R Core Team, 2017). In this study, these data were used for investigating potential associations between individual microbial taxa and serum aqueous metabolites. Spearman’s correlations for differentially regulated microbial OTUs and serum metabolites by PBDEs were performed using R (base, corrplot package) (Wei and Simko, 2016; R Core Team, 2017). Significant gene-metabolite pair correlation plots and bar plots were created using R (base, ggplot2 package) (Wickham, 2009).
Glutathione.
Levels of total glutathione in liver were quantified using a fluorescence-based microplate assay, as previously described (White et al., 2003; McConnachie et al., 2007; Weldy et al., 2011; Scoville et al., 2019).
Quantitative Reverse Transcription–Polymerase Chain Reaction.
Expression of α-1,6-mannosyltransferase (Alg12), glutamate-cysteine ligase catalytic subunit (GCLC), and NAD(P)H dehydrogenase (Nqo1) were measured using quantitative reverse transcription–polymerase chain reaction (RT-qPCR). Total RNA was reverse transcribed to cDNA using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA). The cDNAs were then amplified using Power SYBR Green PCR Master Mix (Applied Biosystems) using a BioRad CFX384 Real-TimePCR Detection System (Bio-Rad, Hercules, CA). Primers were made by Integrated DNA Technologies (Coralville, IA) and primer sequences can be found in Supplemental Table 8.
Results
Effects of PBDEs on Aqueous Metabolites in Serum of CV and GF Mice.
To determine the impact of PBDE exposure on the dynamic systems of intermediary metabolism, aqueous metabolites from several biocompartments were quantified using LC-MS. Of the total 169 aqueous metabolites, 139 were detectable in serum (Supplemental Tables 1 and 2). As shown in Fig. 1, the relative abundances of 13 metabolites in serum were decreased and two increased by lack of intestinal microbiome, PBDE exposure, or both, compared with CV control mice.
Two-way hierarchical clustering dendrogram of aqueous metabolites measured using LC-MS that were differentially regulated across treatment groups in serum of GF and CV mice following exposures to BDE-47 and BDE-99. Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized mean metabolite relative abundances (z-scores) and blue lower. Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE treated groups within enterotype (adjusted P value <0.05, Tukey’s HSD post-hoc test). Caret signs (^) represent statistically significant baseline differences between CV and GF mice. Percent change from CVCO are shown for treatment groups with significant changes.
Six metabolites were decreased in the serum of GF mice under basal conditions, including the tryptophan metabolites 3-indoxylsulfate, 3-indolepropionic acid (3-IPA), the phenylalanine metabolites phenylpyruvic acid and 3-phenyllactic acid, the branched short-chain fatty acid isobutyric acid, and the fatty acid amide N-acetylethanolamine. Acetylglycine and the saturated long-chain fatty acid palmitic acid were increased in serum under basal conditions in GF mice.
PBDEs differentially regulated 12 serum metabolites. Interestingly, the relative abundances of three tryptophan metabolites, namely the microbial metabolites 3-IPA and indole-3-acetic acid, as well as the host metabolite kynurenine, were all decreased by both PBDE congeners in CV but not GF mice. PBDEs also decreased N-acetylethanolamine and phenyllactic acid in CV but not GF mice. BDE-47 decreased the relative abundances of the amino acid valine and its nonproteogenic isomer norvaline, whereas BDE-99 increased phenylpyruvic acid, isobutyric acid, and the carbohydrate metabolite sorbitol also in CV but not GF mice. In contrast, the lack of gut microbiota predisposed GF mice to PBDE-mediated decreases in the relative abundance of metabolites derived from the branched-chain amino acid (BCAA) valine. BDE-47 increased α-ketoisovaleric, whereas both BDE-47 and BDE-99 mediated increased α-hydroxyisovaleric acid.
Effects of PBDEs on Aqueous Metabolites in Liver of CV and GF Mice.
In liver, 156 aqueous metabolites were detectable (Supplemental Table 3). As shown in Fig. 2, under basal conditions, 11 out of the 16 metabolite changes occurred between CV and GF mice. The relative abundances of the sialic acid N-acetylneuraminic acid, the purine metabolite urate, and the amino acid glycine were all increased in GF mice. The relative abundances of isobutyric acid, 3-phenyllactic acid, the amino acid derivative carnitine, the amino acid derivative ketoleucine, the carbohydrate metabolite glucuronic acid, and diet-derived metabolites (benzoic acid, phthalic acid, 4-methoxyphenylacetic acid) were all decreased in GF mice.
Two-way hierarchical clustering dendrogram of aqueous metabolites measured using LC-MS that were differentially regulated across treatment groups in liver of GF and CV mice following exposures to BDE-47 and BDE-99. Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized treatment group mean metabolite relative abundances (z-scores) and blue lower. Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE treated groups within enterotype (adjusted P value <0.05, Tukey’s HSD post-hoc test). Caret signs (^) represent statistically significant baseline differences between CV and GF mice. Percent change from CVCO are shown for treatment groups with significant changes.
In CV mice, BDE-47 had minimal effect on the hepatic aqueous metabolite profiles. However, BDE-99 decreased relative abundances of carnitine, the carbohydrate mannose, the antioxidant glutathione, the homocysteine metabolite cystathionine, and the vitamin nicotinamide, all in CV but not GF mice. The lack of gut microbiota predisposed GF mice to BDE-47-mediated decreases in the glutamine derivative acetyl-l-glutamine.
Effects of PBDEs on Aqueous Metabolites in SIC of CV and GF Mice.
A total of 147 aqueous metabolites were detectable in the SIC (Supplemental Table 4). As shown in Fig. 3, there were 15 metabolites that were decreased in SIC of GF mice under basal conditions, including several amino acids (phenylalanine, tryptophan, arginine, and cysteine, citrulline, asparagine), amino acid derivatives (carnitine, acetyl-l-glutamine, creatine, acetylornithine, and 4-imidazoleacetic acid), the carbohydrate metabolite sorbitol, as well as other organic compounds (folinic acid, and 4-aminophenol). Conversely, the relative abundance of dextrose (also known as d-glucose) was increased in the SIC of GF mice under basal conditions.
Two-way hierarchical clustering dendrogram of aqueous metabolites measured using LC-MS that were differentially regulated across treatment groups in SIC of GF and CV mice following exposures to BDE-47 and BDE-99. Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized treatment group mean metabolite relative abundances (z-scores) and blue lower. Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE-treated groups within enterotype (adjusted P value <0.05, Tukey’s HSD post-hoc test). Caret signs (^) represent statistically significant baseline differences between CV and GF mice. Percent change from CVCO are shown for treatment groups with significant changes.
PBDE exposure differentially regulated six metabolites in the SIC of CV mice but not GF mice. Asparagine, creatine, sorbitol, and cysteine were decreased by both PBDE congeners, folinic acid was decreased by BDE-47, and acetyl-l-glutamine was decreased by BDE-99. In contrast, the lack of gut microbiota predisposed GF mice to BDE-47-mediated increases in the relative abundances of amino acid derivatives (3-phenyllactic acid, acetyl-l-tyrosine), the redox active coenzyme NAD, and the glycerol derivative glyceric acid; BDE-99 mediated an increase in carnitine; and BDE-99 mediated a decrease in 4-hydroxybenzoic acid.
Effects of PBDEs on Aqueous Metabolites in LIC of CV and GF Mice.
A total of 147 aqueous metabolites were detectable in the SIC (Supplemental Table 5). As shown in Fig. 4, PBDEs differentially regulated 26 metabolites in the LIC—the highest among the four biocompartments examined. Interestingly, the relative abundances for the majority of metabolites (24 out of 26) in the LIC were changed by PBDEs only in GF conditions but not in CV conditions. The two exceptions are the tryptophan metabolite indole-3-pyruvic acid, which was reduced in the LIC of CV mice by BDE-47 (the basal relative abundance in GF mice was also lower), and folinic acid, which was increased by BDE-99 in CV but not GF mice.
Two-way hierarchical clustering dendrogram of aqueous metabolites measured using LC-MS that were differentially regulated across treatment groups in LIC of GF and CV mice following exposures to BDE-47 and BDE-99. Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized treatment group mean metabolite relative abundances (z-scores) and blue lower. Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE treated groups within enterotype (adjusted P value <0.05, Tukey’s HSD post-hoc test). Caret signs (^) represent statistically significant baseline differences between CV and GF mice. Percent change from CVCO are shown for treatment groups with significant changes.
The lack of gut microbiome predisposed GF mice to PBDE-mediated increases in 17 LIC metabolites. Compared with vehicle-exposed GF mice, both BDE-47 and BDE-99 increased the relative abundances of the amino acids leucine and methionine, amino acid derivatives (4-hydroxyproline, 5-aminoleulinic acid, and acetylcarnitine), the choline metabolism intermediate and methionine precursor betaine, and the coenzyme A precursor pantothenic acid. BDE-47 increased relative abundances of the nucleic acid 2-deoxycytidine and amino cysteine in LIC of GF mice, whereas BDE-99 increased several amino acids (norleucine, isoleucine, l-alloisoleucine, tyrosine, phenylalanine, and tryptophan), the tryptophan metabolite indole-3-lactic acid, and the carbohydrate metabolite glycerol-3-phosphate.
The lack of gut microbiota also predisposed GF mice to PBDE-mediated decreases in seven metabolites. Both BDE-47 and BDE-99 decreased the relative abundances of carbohydrate metabolites sorbitol and raffinose in the LIC of GF mice compared with control GF mice. In addition, BDE-47 decreased relative abundances of the nucleic acid adenine and carbohydrate metabolite pyruvate, whereas BDE-99 decreased the NAD precursor nicotinic acid, the histamine metabolite 4-imidazoleacetic acid, and the fatty acid amide N-acetylethanolamine in the LIC of GF mice.
Differentially abundant metabolites unique and common to different biocompartments are summarized in Supplemental Fig. 4. Treatment group-specific and common differentially abundant metabolites within each biocompartment are summarized in Supplemental Fig. 5.
Liver Intermediary-Metabolism–Gene Expression and Gene-Metabolite Networks.
As shown in Figs. 5 and 6, 133 intermediary metabolism-related genes belonging to one or more of the KEGG pathways (carbohydrate, glycan, amino acid, vitamin and cofactor, nucleotide, and energy metabolism) were differentially expressed owing to either PBDE exposure, lack of gut microbiota, or both in liver. Specifically, carbohydrate metabolism and glycan biosynthesis and metabolism were represented by 38 and 13 differentially expressed genes, respectively (Fig. 5A). Amino acid metabolism was represented by 70 differentially expressed genes (Fig. 6A) and vitamin and cofactor metabolism contained 44 genes (Fig. 6C). Nucleotide metabolism and energy metabolism contained 11 and 10 differentially expressed genes, respectively (Supplemental Fig. 1). Metabolism category labels for genes in Fig. 5A, Fig. 6, A and C, and Supplemental Fig. 1 can be found in Supplemental Table 6.
Two-way hierarchical clustering dendrogram of carbohydrate and glycan metabolism genes that were differentially regulated across treatment groups in liver of GF and CV mice following exposures to BDE-47 and BDE-99 (A). Gene expression was quantified using RNA-Seq as previously reported (Li et al., 2017) and can be accessed through NCBI Gene Expression Omnibus (GSE101650). Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized mean metabolite relative abundances and blue lower. Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE-treated groups (FDR adjusted P value <0.05). Caret signs (^) represent statistically significant baseline differences between CV and GF mice. Network of differentially regulated carbohydrate (GCK) and glycan (ALG12) metabolism genes (circles) and the significantly altered carbohydrate metabolite mannose (square), which was quantified using LC-MS, in liver (B). Mouse orthologs for human genes are listed in mixed case. Pearson correlation of log2 abundances of mannose and Alg2 across treatment group means (C). Means and S.E. of Alg12 and mannose for each treatment group (D). Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE-treated groups (FDR adjusted P value <0.05, Cuffdiff). Caret signs (^) represent statistically significant baseline differences between CV and GF mice.
Two-way hierarchical clustering dendrogram of amino acid–metabolism genes that were differentially regulated across treatment groups in liver of GF and CV mice following exposures to BDE-47 and BDE-99 (A). Gene expression was quantified using RNA-Seq as previously reported (Li et al., 2017) and can be accessed through NCBI Gene Expression Omnibus (GEO; GSE101650). Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized mean metabolite relative abundances and blue lower. Network of differentially regulated amino acid metabolism genes [dopa decarboxylase (DDC), NQO1, GNMT, CYP1A2, and CYP3A5] (circles) and significantly altered aqueous metabolites (glycine and glutathione) (squares), which were measured using LC-MS, in liver (B). Mouse orthologs for human genes are listed in mixed case. Two-way hierarchical clustering dendrogram of vitamin and cofactor metabolism genes that were significantly altered across treatment groups in liver of GF and CV mice following exposures to BDE-47 and BDE-99 (C). Gene expression was measured using RNA-Seq as previously reported (Li et al., 2017) and can be accessed through NCBI GEO (GSE101650). Euclidian distance and complete linkage were used to generate the dendrogram. Red indicates higher standardized mean metabolite relative abundances and blue lower. Network of differentially regulated vitamin metabolism genes (CYP2D6 and CYP3A4) (circles) and the significantly altered aqueous metabolite nicotinamide (square), which was measured using LC-MS, in liver (square) (D). Asterisks (*) represent statistically significant differences between corn oil-treated and PBDE-treated groups (FDR adjusted P value <0.05). Caret signs (^) represent statistically significant baseline differences between CV and GF mice.
To integrate the liver intermediary metabolism transcriptomic and metabolomic data and identify potential gene-metabolite regulatory networks impacted by PBDEs, network analyses were performed using MetaboAnalyst of all differentially regulated intermediary metabolism-related genes and aqueous metabolites in liver. As MetaboAnalyst was built for human genes, Ensembl Biomart was used to retrieve human homologs.
The Network Explorer module revealed three networks using the list of differentially regulated liver genes and metabolites. One of the networks generated by MetaboAnalyst centered on the carbohydrate metabolite mannose with connections (representing interactions from literature found by Network Explorer) to the carbohydrate metabolism gene glucokinase (GCK), the glycan synthesis gene ALG12, and the carbohydrate mannose (Fig. 5B). In addition, a significant correlation was observed in our data between ALG12 expression and the relative abundance of mannose (Figs. 5C), and BDE-99 decreased both in CV but not GF mice (Fig. 5D). Expression of Alg12 measured using RT-qPCR matched RNA-seq expression patterns across treatment groups (Supplemental Fig. 2).
The second gene-metabolite network that was generated by MetaboAnalyst centered around glutathione (Fig. 6B), which is synthesized from three amino acids (glycine, glutamate, and cysteine). Glutathione was connected to the amino acid metabolism genes dopa decarboxylase (DDC), cytochrome P450 1A2 (CYP1A2 representing the mouse Cyp1a2), CYP3A5 (representing several mouse Cyp3a isomers), and GCLC (representing the mouse Gclc), which is a glutathione-synthesis gene utilizing glutamate and cysteine as substrates (Fig. 6B). Expression of Cyp1a2 and Cyp3a11, which was the most highly expressed Cyp3a isoform, have been validated using RT-qPCR as previously reported (Li et al., 2017). NAD(P)H dehydrogenase, quinone 1 (NQO1), which is classified as a vitamin and cofactor metabolism gene as well as a prototypical target gene of the antioxidative stress sensor nuclear factor erythroid 2–like 2 (Nrf2), is also connected to glutathione (Fig. 6, B and C) (Ross and Siegel, 2018). Levels of liver glutathione were also measured independently using a fluorescence-based microplate assay. Interestingly, the pattern across exposure groups was similar to that found using the LC-MS–based aqueous metabolite panel, suggesting that although the transcriptomic changes in glutathione synthesis/metabolism-related genes occur following PBDE exposure and the absence of gut microbiome, the liver serves as an efficient repertoire to buffer the oxidative stress by maintaining the total glutathione levels (Supplemental Fig. 3, A and B). To validate the RNA-Seq results, RT-qPCR was performed for the mRNA expression of Gclc and Nqo1, and the expression patterns were consistent among the two assays across exposure groups (Supplemental Fig. 3, C–F).
The third network generated was centered on the vitamin nicotinamide, which was connected to vitamin- and cofactor-processing genes CYP2D6 and CYP3A5, represented by several mouse Cyp2d and Cyp3a isoforms, respectively (Fig. 6D).
Correlations between Serum Metabolites and Gut Microbiota.
To investigate the intersection of PBDE regulation of serum metabolites and the gut microbiome at the individual taxa level, the relative abundances of 23 gut microbial taxa that were changed by PBDEs were tested for correlation with PBDE-regulated serum metabolites. Serum metabolites were prioritized because the changes in this biocompartment could have systemic effects. Of the 23 taxa, 19 were significantly correlated with at least one serum aqueous metabolite (Fig. 7).
Pearson correlation matrix of differentially regulated gut-microbiome operational taxonomical units (identified and quantified using 16S rRNA sequencing) and serum metabolites (quantified using LC-MS). Red indicates stronger positive correlations and blue stronger negative.
Regarding amino acid metabolism, the microbial tryptophan metabolite 3-IPA was positively correlated with the largest number of taxa (8), all of which belong to the phylum Firmicutes, and apart from one taxon from the Cloistridia class, all belong to the class Bacilli (Fig. 7). Four of the Bacilli belong to the Bacillales order and three are Lactobacillales. Within the Bacillales order, two taxa were in the Staphylococcaceaes family, one is a Baciliaceae, and the other is unresolved at the family level and beyond. Both Staphylococcaceaes are Staphylococcus at the genus level, but neither is resolved at the species level. The taxon from the Baciliaceae family is Bacillus at the genus level and also unresolved at the species level. Within the Lactobacillale order, two are in the Lactobacillaceae family, and the other taxon is unresolved at the family level and beyond. Both Lactobacillaceae are Lactobacillus at the genus level. One is the reuteri species and the other is unresolved at the species level. The taxon in the Clostridia class is Clostridiales (order), Lachnospiraceae (family), and unresolved at the genus and species levels.
Another microbial tryptophan metabolite, namely indole-3-acetic acid, is positively correlated with two microbial taxa and inversely correlated with one (Fig. 7). The first positively correlated taxon is Firmicutes (phylum), Clostridia (class), Clostridiales (order), Lachnospiraceae (family), Butyrivibrio (genus), and unresolved at the species level. The second positively correlated taxon is Firmicutes (phylum), Erysipelotrichi (class), Erysipelotrichales (order), Erysipelotrichaceae (family), and unresolved at the genus and species level. The inversely correlated taxon has the same taxonomy as the Erysipelotrichaceae (family), belongs to the Allobaculum genus, and is unresolved at the species level.
The phenylalanine metabolite phenyllactic acid is positively correlated with seven gut microbial taxa (Fig. 7). Two of the seven taxa are the same Staphyllococcus genus as described above. An additional two are the same Lactobacillus genus taxa, one identified as reuteri, and the other unresolved at the species level. The fifth taxon associated with phenyllactic acid is Bacillales at the order level and is also previously described above. One of the remaining taxa associated with phenyllactic acid is Actinobacteria (phylum), Coriobacteriia (class), Coriobacteriales (order), Coriobacteriaceae (family), Adlercreutzia (genus), and unresolved at the species level. The other taxon is Firmicutes (phylum), Clostridia (class), Clostridiales (order), Mogibacteriaceae (family), and unresolved at the genus and species levels.
BCAAs valine and norvaline are positively correlated with the same taxon from the Butyrivibrio genus as indole-3-acetic acid (Fig. 7). They are both also inversely correlated with the same taxon from the Allobaculum genus as indole-3-acetic acid (Fig. 7).
The fatty acid amide N-acetylethanolamine is positively correlated with four microbial taxa (Fig. 7). The full taxonomy of three have already been described: Adlercreutzia (genus), Lactobacillus (genus, unresolved at species level), and Mogibacteriaceae (family). The fourth taxa associated with N-acetylethanolamine has the same taxonomy as the Lactobacillus (genus), and also unresolved at the species level, but is distinct from the Lactobacillus (genus) associated with 3-IPA and 3-phenyllactic acid.
Three metabolites, namely the branched short-chain fatty acid isobutyric acid, the phenylalanine metabolite and phenyllactic acid precursor phenylpyruvic acid, and the carbohydrate sorbitol, are all positively correlated with a microbial taxon that is Firmicutes (phylum), Clostridia (class), Clostridiales (order), Lachnospiraceae (family), Dorea (genus), and unresolved at the species level (Fig. 7). These same three metabolites are also inversely associated with Akkermansia (genus) mucinophila (species), which belongs to Verrucomicrobia phylum, Verrucomicrobiae class, Verrucomicrobiales (order), and Verrucomicrobiaceae family.
Discussion
Although the commercial pentaBDE (composed mainly of BDE-47 and BDE-99) and octaBDE formulations of PBDEs were voluntarily taken off the market in 2004 and decaBDE by 2014 (2017), PBDEs are still readily detectable in human samples. Although concentrations of pentaBDE constituents have previously been observed to decline since their withdrawal, the trend may have flattened out and concentrations may now be on the rise (Ma et al., 2013; Zota et al., 2013; Hurley et al., 2017; Parry et al., 2018; Cowell et al., 2019). This is probably the result of the persistent and bioaccumulative nature of PBDEs, the recycling of PBDE-containing products worldwide, and the potential for more highly brominated PBDEs to degrade in the environment into lower brominated forms (Bezares-Cruz et al., 2004; Zeng et al., 2008).
As originally reported in Li et al. (2017), male mice were dosed once daily for four consecutive days with 48.5 mg/kg BDE-47 and 56.5 mg/kg BDE-99 (100 μmol/kg). Although these doses are high compared with estimated human PBDE exposures, they were selected to allow for comparison of acute xenobiotic metabolism gene expression changes with previously published studies (Pacyniak et al., 2007; Sueyoshi et al., 2014; Li et al., 2017). In this study, metabolomic data from archived samples collected by Li et al. was integrated with previously reported liver gene expression (Li et al., 2017) and intestinal microbiome data (Li et al., 2018) to explore acute PBDE effects on intermediary metabolism.
Limitations of this study include sample size and the evaluation of only male mice only, and that it was mainly associative, all of which stem from the use of archived samples and previously reported datasets. As reported in Li et al. (2017), male mice were used exclusively owing to changes in hormone levels during the estrous cycle. Such hormonal variation is also relevant to this study as it could independently impact glucose metabolism and other intermediary metabolic pathways (Mauvais-Jarvis et al., 2013). Sample size for the three data sets used in this study was also a limitation. However, even with relatively small sample sizes, we were able to identify statistically significant PBDE and gut microbiome–associated changes in intermediary metabolites, expression of associated genes, and correlations among metabolites and several previously reported human health and disease relevant microbes (Li et al., 2018). Importantly, future mechanistic studies will be needed to fully validate the observed findings of this study.
Tryptophan Pathway and Involvement of PBDEs in Metabolic Syndrome.
Notably, serum 3-IPA, which is a tryptophan microbial metabolite known to be protective against diabetes and inflammation, was decreased by both PBDE congeners in CV mice but not in GF mice (Fig. 1). This highlights the potential involvement of the tryptophan metabolism pathway in PBDE-mediated metabolic disorders, because 3-IPA is produced by intestinal bacteria and accumulates in the blood to modulate systemic effects (Dodd et al., 2017) such as metabolic syndrome. Specifically, 3-IPA is associated in epidemiology studies with reduced risk of type 2 diabetes and may preserve pancreatic β cell function (de Mello et al., 2017). These observations in humans were further supported in rats because 3-IPA improved glucose metabolism, as evidenced by reduced fasting plasma insulin and the homeostatic model assessment index of insulin resistance (Abildgaard et al., 2018). In humans, BDE-47 was associated with diabetes in two independent community-based epidemiologic studies (Zhang et al., 2016), and total PBDEs were linked to increased risk of diabetes during pregnancy (Eslami et al., 2016; Smarr et al., 2016). PBDE-exposed rats developed hyperglycemia and aberrant expression of genes involved in glucose homeostasis (Zhang et al., 2016; Krumm et al., 2018), corresponding to disturbed glucose and insulin metabolism in adipose tissue, disturbed lipid metabolism in liver, and increased ketone concentrations in blood (Cowens et al., 2015). PBDEs disturbed glucose signaling in hepatocytes (Søfteland et al., 2011). In vitro studies showed that PBDEs disrupted insulin-producing pancreatic β cells (Karandrea et al., 2017). Therefore, PBDE-mediated decrease in 3-IPA, by disrupting glucose and insulin signaling, may be among the toxification mechanisms for the pathogenesis of diabetes.
Another potential involvement of 3-IPA in metabolic syndrome is inflammation—an important contributor to diabetes (Donath and Shoelson, 2011; Goldfine et al., 2011)—as 3-IPA is known to suppress inflammatory genes by antagonizing aryl hydrocarbon receptor functions (Yisireyili et al., 2017). In addition, 3-IPA is a pregnane X receptor (PXR) ligand in intestine, and activation of PXR reduces intestinal inflammation (Cheng et al., 2012; Venkatesh et al., 2014). Conversely, PBDEs have been shown to affect the immune system during development (Zhou et al., 2006; Martin et al., 2007; Liu et al., 2012; Bondy et al., 2013; Koike et al., 2014; Lv et al., 2015), as well as promote inflammation and oxidative stress (Zhang et al., 2013), which are key contributors for the pathogenesis of diabetes (Farah et al., 2008; Vikram et al., 2014; Domingueti et al., 2016; Hojs et al., 2016). Therefore, a PBDE-mediated decrease in IPA may also contribute to the pathogenesis of diabetes by promoting inflammation.
Branched-Chain Amino Acid Metabolites, PBDEs, and Adipose Tissue Metabolism.
Isobutyric acid is a branched short-chain fatty acid derived from the BCAA valine and can be produced by gut bacteria or by host (Liebich and Först, 1984; Zarling and Ruchim, 1987; Cole, 2015) . Shown to potentially benefit adipocyte metabolism, isobutyric acid was decreased in this study in GF mice under basal conditions and by BDE-99 in CV mice (Heimann et al., 2016) . Although there is controversy as to a potentially beneficial or detrimental role for BCAAs in obesity and type 2 diabetes, preliminary evidence suggests a potential protective role for their metabolites(Lynch and Adams, 2014; Heimann et al., 2016) ).
In obesity, adipocytes exhibit inter-related changes, including increased basal rates of lipolysis and insulin resistance, that are thought to contribute to systemic insulin resistance (Girousse et al., 2013; Arner and Langin, 2014; Morigny et al., 2016). Treatment of rat adipocytes with isobutyric acid reduced lipolysis and improved basal and insulin-induced glucose uptake (Heimann et al., 2016) (). Nonbranched short-chain fatty acids have also been observed to have similar effects on rat adipocytes (Heimann et al., 2014). In contrast, PBDEs have been associated with increased rat adipocyte lipolysis and decreased glucose oxidation (marker of uptake and metabolism) (Hoppe and Carey, 2007). Furthermore, PBDEs have been associated with insulin resistance among obese humans (Helaleh et al., 2018).
The decreases of the potentially beneficial valine metabolite isobutyric acid and, among GF mice, its upstream metabolic intermediates α-ketoisovaleric and α-hydroxyisovaleric acids, suggests that PBDEs may promote metabolic syndrome in part by interference with valine/BCAA metabolism.
PBDEs, Protein Glycosylation, and Metabolic Syndrome.
Liver intermediary metabolism–gene expression and aqueous metabolite data were integrated using the Network Explorer module in MetaboAnalyst (Chong et al., 2018). One of the resulting networks, consisting of the genes ALG12, GCK, and the carbohydrate mannose, stood out owing to a correlation found in the study data between Alg12 and mannose. As part of the N-glycosylation pathway, ALG12 is responsible for the addition of the eighth mannose residue to a growing oligosaccharide that ultimately possesses nine mannoses and three glucoses, all of which are added in the endoplasmic reticulum lumen (Burda and Aebi, 1999; Thiel et al., 2002; Aebi, 2013). As an N-glycan, this molecule is subsequently attached to the nitrogen of the side chain of asparagine in a protein by oligosaccharyltransferase (Burda and Aebi, 1999; Aebi, 2013). N-glycans on proteins participate in the luminal endoplasmic reticulum–associated degradation pathway that detects and degrades misfolded proteins (Quan et al., 2008).
Epidemiologic studies have found that patients diagnosed with metabolic syndrome have an altered N-glycan profile, with some components, including those with seven and eight mannose residues, being positively associated with signs of metabolic syndrome and some inversely (Lu et al., 2011; McLachlan et al., 2016). In this study, BDE-99 reduced Alg12 expression and the relative abundance of mannose in CV but not GF mice, suggesting that a PBDE microbial metabolite could also contribute to these changes. To our knowledge, the potential for PBDEs and/or PBDE host/microbial metabolites to interfere with protein glycosylation has not been explored. However, future glycomics studies will be necessary to evaluate this potential novel mechanism of PBDE toxicity.
Microbiota Correlations with Serum Metabolites.
The present study showed that PBDE-mediated decreases in Lactobacillus reuteri correlate with serum decreases in the tryptophan metabolite 3-IPA. Even though L. reuteri is not known to produce 3-IPA, it does produce indole-3-lactic acid, which is a 3-IPA precursor (Dodd et al., 2017). PBDEs also reduced serum indole-3-acetic acid, which is another tryptophan metabolite produced by L. reuteri (Cervantes-Barragan et al., 2017). Interestingly, in mouse models into which L. reuteri has been introduced in combination with tryptophan supplementation, indole-3-acetic and indole-3-lactic acid were shown to promote the development of intestinal T cells that improve immune tolerance to antigens in the diet (Cervantes-Barragan et al., 2017). However, given the small sample size, microbe-metabolite associations should be viewed as preliminary until follow up studies with larger numbers are completed.
Summary.
The present study demonstrates that oral PBDE exposure profoundly altered many important intermediary metabolites involved in amino acid, carbohydrate, and lipid metabolism, in a biocompartment-specific manner. Although liver metabolites were influenced more by enterotype than PBDEs, metabolites from large intestinal contents were influenced more by PBDEs than by enterotype and such influence occurred almost exclusively in GF conditions. Metabolites from serum and SIC were regulated by both enterotype and PBDEs. Two interesting regulatory patterns were discovered: 1) The presence of gut microbiota is necessary for PBDE-mediated decreases in metabolites, such as those from the tryptophan pathway in serum and those from the metabolism of amino acids and carbohydrates in the SIC; 2) lack of gut microbiota predisposed mice to PBDE-mediated changes in metabolites, including the majority of those from amino acid and carbohydrate metabolism in the LIC. Correlation analyses among targeted metabolomics, hepatic RNA-Seq, and intestinal 16S rDNA sequencing further established potential mechanisms for PBDE-mediated regulation of intermediary metabolites. Together our observations highlight the potential involvement of dysregulation of intermediary metabolism as part of the mechanisms underlying the multiorgan toxicities of PBDEs and unveils the important role of gut microbiota as a key internal modifier for PBDE-mediated host response. Follow-up studies will include mechanistic investigation using male and female wild-type and nuclear receptor knockout mice to assess the impact of key metabolites identified in this study on PBDE toxicity, dysbiosis, and potential interactions with nuclear receptors. Inflammation, glucose homeostasis, and the effects of high-fat diet will also be explored.
Acknowledgments
The authors thank GNAC technical support for the germ-free mouse experiment, Dr. Debbie Nickerson’s laboratory at UW Genome Sciences for sequencer access and support for the A-Seq experiments, Brian High from UW DEOHS IT team for providing server access and advice to improve our A-Seq bioinformatics pipeline, the Northwest Metabolomics Research Center for targeted metabolomics, Dianne Botta and Collin White from Dr. Terrance Kavanagh’s laboratory for their guidance measuring glutathione, and the members of Dr. Cui’s laboratory for help in tissue collection and manuscript revision.
Authorship Contributions
Participated in research design: Li, Cui.
Conducted experiments: Li, Dempsey, Wang, Gu, Cui.
Contributed new reagents or analytic tools: Raftery.
Performed data analysis: Scoville, Li, Wang, Gu, Cui.
Wrote or contributed to the writing of the manuscript: Scoville, Li, Dempsey, Wang, Raftery, Mani, Gu, Cui.
Footnotes
- Received January 23, 2019.
- Accepted May 7, 2019.
↵* Primary laboratories of origin.
This work was supported by the National Institutes of Health [Grants R01 GM111381, R01 ES030197, R01 ES025708], the University of Washington Center for Exposures, Diseases, Genomics, and Environment [P30 ES0007033], and the Sheldon Murphy Endowment.
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- ACN
- acetonitrile
- Alg12
- α-1,6-mannosyltransferase
- BCAA
- branched-chain amino acid
- BSCFA
- branched short-chain fatty acid
- CV
- conventional
- FDR
- false discovery rate
- Gck
- glucokinase
- Gclc
- glutamate-cysteine ligase catalytic subunit
- GEO
- Gene Expression Omnibus
- GF
- germ-free
- Gnmt
- glycine N-methyltransferase
- 3-IPA
- 3-indoleproprionic acid
- LIC
- large intestinal contents
- Nqo1
- NAD(P)H dehydrogenase, quinone 1
- PBDE
- polybrominated diphenyl ether
- RT-qPCR
- quantitative reverse transcription–polymerase chain reaction
- SIC
- small intestinal contents
- Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics
