Persistent organic pollutants modify gut microbiota-host metabolic homeostasis in mice through aryl hydrocarbon receptor activation.
Obesity risk factors include changes in diet and lifestyle and, in some rare instances, can be explained by genetic predisposition (Jaaskelainen et al. 2013); however, these factors alone seem unlikely to explain the growing worldwide obesity epidemic. There is increasing evidence that chronic exposure to environmental chemicals through the diet, especially persistent organic pollutants (POPs), may promote the development of obesity and type 2 diabetes in humans (Lee et al. 2014). Of particular interest is the role of the aryl hydrocarbon receptor (AHR), which is bound and activated by a variety of POPs including coplanar polychlorinated biphenyls (PCBs) and halogenated aromatic hydrocarbons, in mediating the toxic effects of these compounds at high doses (Alaluusua and Lukinmaa 2006; Ren et al. 2011) and promoting obesity at low doses (Arsenescu et al. 2008; Howell and Mangum 2011). More recent evidence suggests that gut microbiota, which can be modulated by the AHR (Jin et al. 2014), may play a pivotal role in POP-induced obesity (Myre and Imbeault 2014).
The AHR is a xenobiotic sensor that was recently shown to be a pivotal regulator of gut and immune homeostasis (Monteleone et al. 2012). AHR agonists, including both endogenous (e.g., several tryptophan catabolites), dietary (e.g., indole-3-carbinol metabolite indolo[3,2b] carbazole), and xenobiotic [e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)] ligands, can dramatically impact host metabolism and immunity (Angrish et al. 2013; Bjeldanes et al. 1991; Claus et al. 2008; Zelante et al. 2013). For example, microbiota-derived tryptophan metabolites, including indole-3-aldehyde, were shown to modulate the AHR-IL-22 axis, thus impacting mucosal immune homeostasis in the gut (Zelante et al. 2013). Furthermore, alteration of gut microbiota through dietary manipulation and the presence of environmental contaminants may disturb physiological homeostasis, leading to various diseases such as obesity (Greenblum et al. 2012), diabetes (Mathis and Benoist 2012), and inflammatory bowel disease (Nagalingam and Lynch 2012). Recent studies suggest that interactions between gut microbiota and environmental toxicants may contribute in part to the development of obesity and diabetes (Snedeker and Hay 2012). Therefore, bacteria in the gut may be affected by dietary exposure to TCDD and other TCDD-like compounds, thus altering not only bacterial populations but also their inherent metabolic activity (Maurice et al. 2013). In the present study, we used TCDF instead of TCDD, given its considerably shorter half-life in rodents (DeVito and Birnbaum 1995), and we used doses higher than those of polychlorinated dibenzo-p-dioxins and dibenzofurans typically found in the diet of adults [0.3-3.0 pg toxic equivalents (TEQs)/kg body weight] (Nakatani et al. 2011; Schecter et al. 1994). Further, TCDF was provided through the diet over 5 days to represent the major route of exposure in humans (Jones and Bennett 1989).
In the present study we used a combination of 16S rRNA gene sequencing, metabolite profiling [ultra-performance liquid chromatography coupled with triplequadrupole mass spectrometry (UPLC-TQMS)], and [sup.1]H nuclear magnetic resonance (NMR)-based metabolomics to investigate the alteration of gut microbiota and the host metabolome in mice treated with TCDF through the diet. In addition, we investigated possible correlations between gut microbiota composition and signaling pathways, such as AHR and farnesoid X receptor (FXR), in response to TCDF exposure with the goal of identifying a specific host-microbiota signaling axis. Our findings provide new evidence that exposure to persistent environmental contaminants mediated through activation of the AHR strongly impacts the gut microbiota and host metabolism.
Materials and Methods
Animals and diets. Animal experimental procedures were performed according to the National Institutes of Health (NIH) guidelines and were reviewed and approved by the Pennsylvania State University Institutional Animal Care and Use Committee. Mice were treated humanely and with regard for the alleviation of suffering. A total of twelve 8-week-old male C57BL/6J mice ([Ahr.sup.+/+]) were purchased from the Jackson Laboratory and housed together in groups of two to four mice per cage with a 12-hr light/dark cycle and constant temperature (23 [+ or -] 1[degrees]C) and humidity (45-65%); mice were allowed free access to food and water. Bedding material consisted of Alpha Dri (Shepherd Specialty Papers). After acclimatization for 1 week to normal flora in our facility, mice were randomly divided into two groups of six and trained to eat transgenic dough pills (Bio-Serve) that contained no TCDF (Cambridge Isotope Laboratories Inc.) for 4 days. Dough pills containing TCDF (acetone solution as vehicle) or acetone alone were prepared using tablet molds, and one pill uniformly contained 0.6 [micro]g TCDF, thus providing a final dose of 24 [micro]g/kg. In the TCDF studies, mice were fed the dough pills containing TCDF continuously for 5 days in the morning (one pill per mouse per day). Urine and feces were collected every other day (i.e., days 2, 4, and 6) over the TCDF treatment period. Blood, liver, intestine, cecum, and cecal content samples were collected immediately following C[O.sub.2] asphyxiation on day 7. All samples were stored at -80[degrees]C until analysis. Ten male [Ahr.sup.-/-] mice (C57BL/6J congenic mice) were obtained from C. Bradfield (University of Wisconsin) and inbred to maintain the line. All mice were killed by C[O.sub.2] asphyxiation at 1300 hours, after they were transported to the laboratory.
Quantitative real-time polymerase chain reaction (qPCR). RNA was extracted from frozen liver and intestine (~ 50 mg) using TRIzol reagent (Invitrogen). cDNA was synthesized from 1 [micro]g of total RNA using Superscript II reverse transcriptase (Invitrogen), and the products were diluted to 1:10 before use in subsequent reactions. Gene-specific primers were used in each reaction, and all results were normalized to ribosomal protein L13 or [beta]-actin mRNA (for primer sequences, see Supplemental Material Table S1). qPCR assays were carried out using SYBR Green qPCR Master Mix on an ABI Prism 7900HT Fast Real-Time PCR sequence detection system (Applied Biosystems). The reactions were analyzed according to the [DELTA][DELTA]CT method. qPCR conditions were 40 cycles of 95[degrees]C for 20 sec, 95[degrees]C for 0.01 sec, 60[degrees]C for 20 sec, 95[degrees]C for 15 sec, 60[degrees]C for 15 sec, and 95[degrees]C for 15 sec.
16S rRNA gene sequencing of the micro biota. The bacteria in the cecal contents were extracted using the PowerSoil DNA Isolation Kit (MO BIO Laboratories Inc.). The PCR products (~ 1,000 bp) were purified using Agencourt AMPure technology (Beckman Coulter) as described for the Short Fragment Removal Procedure (Roche Diagnostics 2013). After purification, the products were quantified by both Qubit (Lifetech) and qPCR, using the KAPA Biosystems Library Quantification Kit (Kapa-Biosystems). Products were pooled based on molar concentrations, subjected to electrophoresis on a 1% agarose gel, and extracted. After purification with a QIAquick PCR Purification kit (QIAGEN), the quality and quantity of the preparations were assessed using a DNA 7500LabChip on the Agilent 2100 Bioanalyzer (Agilent Technologies) and Qubit quantification. Sequencing was performed using a quarter Picotiter Plate on a 454 Life Sciences Genome Sequencer FLX (Roche Diagnostics). 16S rRNA gene sequencing analysis was performed using the mothur platform, as previously described (Schloss et al. 2009).
Bile acid profiling by TQMS. Targeted analysis of bile acids was performed using an Acquity UPLC system coupled to a Waters Xevo TQS MS (Waters) with a C18 BEH column (2.1 x 100 mm, 1.7 [micro]m; Waters). External bile acid standards, including cholic acid (CA), lithocholic acid (LCA), chenodeoxycholic acid (CDCA), deoxycholic acid (DCA), ursodeoxycholic acid (UDCA), [alpha]-muricholic acid ([alpha]MCA), [beta]-muricholic acid ([beta]MCA), [omega]-muricholic acid ([omega]MCA), glycocholic acid (GCA), glycochenodeoxycholic acid (GCDCA), taurocholic acid (TCA), taurodeoxycholic acid (TDCA), tauro chenodeoxycholic acid (TCDCA), taurolithocholic acid (TLCA), tauro-[alpha]-muricholic acid (T[alpha]MCA), and tauro-[beta]-muricholic acid (T[beta]MCA), were purchased from Sigma-Aldrich. The deuterated internal standards [cholic acid-2,2,4,4-D4 (CA-d4), lithocholic acid-2,2,4,4-D4 (LCA-d4), ursodeoxycholic acid-2,2,4,4-D4 (UDCA-d4), chenodeoxycholic acid-2,2,4,4 (CDCA-d4), and deoxycholic acid-2,2,4,4-D4 (DCA-d4)] were also obtained from Sigma-Aldrich. Tissues and feces were weighed (25 mg) and internal standards (5 [micro]L, 1 [micro]g/mL) were added before extraction with 500 [micro]L ethanol by homogenization, followed by shaking at 60[degrees]C for 30 min. After cooling, the samples were incubated at 100[degrees]C for 3 min and centrifuged at 1,600 x g for 10 min at 4[degrees]C, and supernatants were transferred to a new tube. The pellets were extracted twice with 500 [micro]L ethanol. After centrifugation at 11,200 x g for 1 min at 4[degrees]C, the combined supernatants were dried down and later dissolved with 500 [micro]L methanol. Following vortexing and centrifugation, 300 [micro]L of the supernatant was transferred to an autosampler vial. Analytes were detected by selected ion monitoring (SIM) or multiple reaction monitoring (MRM) and normalized by internal standard. Analytes of CA, LCA, UDCA, CDCA, and DCA were quantified with their respective deuterated internal standards. Others were quantified with CA-d4. The results were calculated according to individual standard curves established as follows: [area.sub.analyte]/ [area.sub.internal std]. Retention times of analytes and their SIM/MRM are presented in Supplemental Material, Table S2.
Fecal IgA and serum LAL (Limulus amebocyte lysate) measurements. We obtained IgA mouse enzyme-linked immunosorbent assay (ELISA) kits from Abcam, and Pierce LAL chromogenic endotoxin quantitation kits and lipocalin-2 (LCN2) ELISA kits from Thermo Scientific. Fecal IgA, LCN2, and serum LAL were measured according to manufacturer instructions.
Scanning electron microscopy (SEM). SEM images were used for observation of gut microbiota morphology, particularly to visualize bacteria associated with the intestinal epithelium. The ileum was opened and fixed with 2% glutaraldehyde overnight at 4[degrees]C. Samples were washed with 0.1 M sodium cacodylate buffer, pH 7.4, three times and then postfixed with 1% osmium tetroxide in 0.1 M sodium cacodylate buffer. The samples were washed and dehydrated through a graded ethanol series, then critical-point dried, mounted, and coated with 10 nm gold palladium alloy. Images were acquired on a JSM 5400 SEM (JEOL).
TCDF reporter assay. TCDF treatment in [Ahr.sup.-/-] mice was confirmed using the Hepa1.1 mouse hepatoma reporter cell line. After killing the [Ahr.sup.-/-] mice, 50 mg of liver was extracted using the Folch method (Folch et al. 1957), lyophylized, and reconstituted in 20 [micro]L of DMSO; 0.1 [micro]L of the DMSO-reconstituted extracted was used for the reporter assay. The Hepa 1.1 mouse hepatoma cell line containing the stably integrated pGudluc 1.1 luciferase reporter construct was obtained from M. Denison (University of California, Davis). The Hepa 1.1 mouse hepatoma cell line was maintained in [alpha]-modified essential media (Sigma-Aldrich) supplemented with 8% fetal bovine serum (Hyclone Laboratories) and 100 IU/mL penicillin/100 [micro]g/mL streptomycin (Sigma-Aldrich). Cells were cultured at 37[degrees]C in a humidified atmosphere composed of 95% air and 5% C[O.sub.2]. The reporter cells (Hepa 1.1) were seeded in 12-well plates and cultured to approximately 90% confluence. The AHR responsiveness of extracts was examined after treatment for 4 hr, followed by lysing with 400 [micro]L lysis buffer.
Histopathology and clinical biochemistry. Formalin-fixed liver tissue was embedded in paraffin, sectioned (3-4 [micro]m), stained with hematoxylin and eosin (H&E), and examined under a microscope by a qualified liver pathologist (N.T.). Serum alanine transaminase (ALT) and alkaline phosphatase (ALP) were measured using the VetScan VS2 and the Mammalian Liver Profile rotor (Abaxis) according to the manufacturer's instructions.
NMR-based metabolomics experiment. Sample preparation. Liver and intestinal tissues (~ 50 mg) were extracted three times with 600 [micro]L of precooled methanol-water mixture (2/1, vol/vol) using the PreCellys Tissue Homogenizer (Bertin Technologies). Fecal and cecal content samples (~ 50 mg) were subjected to three consecutive freeze--thaws and directly extracted by precooled phosphate buffer with homogenization using the Precellys Tissue Homogenizer. After extraction, 550 [micro]L of each extract was centrifuged and then transferred to a 5-mm NMR tube (for more detailed methods, see Supplemental Material, "NMR-based metabolomics experiment: Sample preparation").
[sup.1]H NMR spectroscopy. [sup.1]H NMR spectra of all the biological samples were recorded at 298 K on a Bruker Avance III 600 MHz spectrometer (operating at 600.08 MHz for [sup.1]H) equipped with a Bruker inverse cryogenic probe (Bruker Biospin). Typical one-dimensional NMR spectra were acquired for each of the samples employing the first increment of NOESY (Nuclear Overhauser effect spectroscopy) pulse sequence (NOESYPR1D). To facilitate NMR signal assignments, a range of two-dimensional (2D) NMR spectra were acquired and processed for selected samples (for more detailed methods, see Supplemental Material, "NMR-based metabolomics experiment: [sup.1]H NMR Spectroscopy").
NMR data processing and multivariate data analysis. [sup.1]H NMR spectra were corrected manually for phase and baseline distortions, and spectral region [delta] 0.5-9.5 was integrated into regions with equal width of 0.004 ppm (2.4 Hz) using AMIX software (V3.8; Bruker-Biospin). Multivariate data analyses, including principal component analysis and orthogonal projection to latent structures with discriminant analysis (OPLS-DA), were performed on the NMR data using SIMCA-P+ software (version 13.0; Umetrics). To facilitate interpretation of the results, we performed color-coded correlation coefficient plots with back-transformation of the loadings generated from the OPLS-DA using an in-house-developed script for MATLAB (Mathworks) (see Supplemental Material, "NMR-based metabol omics experiment: Spectral data processing and multivariate data analysis").
Data analysis. All experimental values are presented as mean [+ or -] SD. Graphical illustrations and statistical analysis were performed using GraphPad Prism (version 6.0; GraphPad). p-Values < 0.05 were considered as significant.
Effects of dietary TCDF on liver enzymes and morphology. The effect of sustained AHR activation on the gut was assessed using a potent AHR ligand, TCDF. [Ahr.sup.+/+] and [Ahr.sup.-/-] mice treated with dietary TCDF (24 [micro]g/kg body weight) exhibited mild and no histopathological changes in the liver, respectively (Figure 1A,B). Mildly elevated serum ALT (19.3-37.5 U/L) and ALP (36.5-52.2 U/L) levels (Figure 1C), but no significant differences in body weight (Figure 1D), were observed in the [Ahr.sup.+/+] mice after TCDF treatment, thus indicating minimal hepatic toxicity. Further, no evidence of steatosis, bile duct injury, or cellular degeneration was found, as is often observed with high-dose TCDD treatment (Ozeki et al. 2011). [Ahr.sup.-/-] mice exhibited no significant changes in the liver enzymes ALT and ALP (Figure 1C). TCDF exposure in the [Ahr.sup.-/-] mice was confirmed using dichloromethane extracts of liver and assessing the induction of AHR signaling in a Hepa 1.1 stable reporter cell line (see Supplemental Material, Figure S1). Transcriptional targets of AHR, including Cyp1al, Cyp1a2, and Cyp2e1, were significantly induced in the liver and across the small intestine (Figure 1E) and were AHR dependent (Figure 1F). These observations suggest that the subsequent metabolic changes described below were specific to AHR activation and not due to overt liver toxicity.
Effects of dietary TCDF on the morphology, population, and composition of gut microbiota. Given recent evidence that the AHR is important for regulating gut homeostasis (Monteleone et al. 2013), and that in gut microbiota changed after oral dietary exposure to POPs (Snedeker and Hay 2012), we examined this relationship. Weighted UniFrac principal coordinate analysis (for assessing changes in abundance) of 16S rRNA sequencing results indicated that dietary TCDF induced a remarkable change in the overall population of gut microbiota (Figure 2A). Firmicutes and Bacteroidetes exhibited significant changes, with reduction of Firmicutes/Bacteroidetes ratio after dietary TCDF exposure (Figure 2B; see also Supplemental Material, Figure S2). Furthermore, cecal contents from TCDF-treated mice were enriched with Butyrivibrio spp. but depleted in Oscillobacter spp., compared with those from vehicle-treated mice (Figure 2C). Further, we observed an increase in the class Flavobacteria but a decrease in the class Clostridia. SEM images showed that dense segmented filamentous bacteria (SFB) formed a network of segmented filaments in the ileum of vehicle-treated mice, whereas they were dramatically depleted in the ileum of TCDF-treated mice (Figure 2D). qPCR analyses revealed that dietary TCDF significantly reduced SFB levels in the ileum, with no significant changes observed between TCDF-treated [Ahr.sup.-/-] mice and vehicle-treated [Ahr.sup.-/-] mice (Figure 2E).
Inflammatory signaling and immune responses. TCDF led to a significant increase in mRNAs that encode factors in the ileum involved in inflammatory signaling, such as Il-1[beta], Il-10, Tnf-[alpha], Saa1, and Saa3 (Figure 3A); these changes were AHR dependent (Figure 3B). Lcn-2 mRNA, encoding lipocalin-2, a sensitive biomarker for intestinal inflammation (Chassaing et al. 2012), was dramatically increased in the ileum (Figure 3C) and in the feces (Figure 3D) of mice exposed to TCDF. No significant change in Lcn-2 mRNA (Figure 3C) was observed in the ileum of [Ahr.sup.-/-] mice after TCDF treatment. Dietary TCDF resulted in a significant elevation of serum lipopoly-saccharide (LPS) (Figure 3F) and fecal IgA (Figure 3G), and also resulted in decreased Myosin Vb and Ptprh in the ileum (Figure 3E), both of which are closely associated with gut permeability and the host immune system. Taken together, these results suggest that dietary TCDF triggers robust intestinal inflammation and inflammatory signaling in mice in an AHR-dependent manner.
Bile acid metabolism. It is well known that gut microbiota has profound effects on bile acid metabolism (Nicholson et al. 2012). To determine whether the composition of bile acids was affected by dietary TCDF exposure, we used UPLC-TQMS to examine bile acids and their levels in the liver, intestine, cecum, and feces, according to published methods (Jiang et al. 2015). Compared with vehicle-treated mice, TCDF-treated [Ahr.sup.+/+] mice had significantly higher levels of DCA in the small intestine and feces (Figure 4A and 4C), and they also had significantly higher levels of three conjugated bile acids (GCA, TCDCA, and T[beta]MCA) (Figure 4B,D) that were affected by gut microbial metabolism, as previously reported (Sayin et al. 2013). The size of the bile acid pool was significantly increased in TCDF-treated mice (Figure 4E). In addition, levels of LCA (intestine and cecum), CDCA (liver), and TLCA (liver and feces) in TCDF-treated mice were higher than those in vehicle-treated mice (see Supplemental Material, Figure S3A-E).
Effects of dietary TCDF on gene expression profiles and FXR signaling. FXR serves as a critical modulator of enterohepatic circulation and thus plays a key role in the regulation of bile acid synthesis and homeostasis (Claudel et al. 2005). TCDF treatment was associated with significant down-regulation of Fxr mRNA and its target gene small heterodimer partner (Shp) mRNA in both the ileum and liver of mice (Figure 4F,G). Intestinal fibroblast growth factor 15 (Fgf15) mRNA, which encodes a growth factor that is released by small intestine epithelial cells, was also markedly decreased by TCDF treatment (Figure 4F). In addition, no significant differences in the expression of Fxr, Fgf15, and Shp in the ileum (Figure 4F) or liver (Figure 4G) were observed between TCDF-treated and vehicle-treated [Ahr.sup.-/-] mice. Taken together, these observations indicate that TCDF affects FXR signaling pathways through an AHR-dependent mechanism.
Effects of dietary TCDF on gene expression profiles: genes associated with bile acid synthesis, conjugation, and reabsorption. Next, we addressed the question of whether altered bile acid profiles in TCDF-treated [Ahr.sup.+/+] mice were associated with changes in gene expression of enzymes involved in bile acid synthesis and/or transport. TCDF treatment resulted in a 10-fold elevation in the expression of Cyp7a1 mRNA in the liver (Figure 4H), and CYP7A1 protein levels were also increased (see Supplemental Material, Figure S4). Expression of Cyp7b1 mRNA was also increased by 1.8-fold in the liver (Figure 4H). In contrast, no significant change in Cyp7a1 mRNA expression was observed in TCDF-treated [Ahr.sup.-/-] mice (see Supplemental Material, Figure S3E). The expression of mRNAs encoded by other genes, such as Cyp8b1, Cyp27a1, and Akr1dl, all involved in bile acid synthesis, were not significantly affected by TCDF exposure (Figure 4H).
We observed a significant decrease in mRNAs encoding the taurine biosynthetic enzymes cysteine dioxygenase (CDO) and cysteine sulfinate decarboxylase (CSD), the rate-limiting enzyme for the conversion of cysteine to taurine (Bitoun and Tappaz 2000a, 2000b), in the livers of TCDF-treated [Ahr.sup.+/+] mice (Figure 4I). These observations were confirmed by the TCDF-induced down-regulation of taurine levels in gut tissues revealed by [sup.1]H NMR metabolomics (see Supplemental Material, Figure S6A-C and Table S3). However, expression of taurine transporter (Taut) mRNA, the main regulator of intracellular taurine levels (Bitoun and Tappaz 2000b), was up-regulated in the livers of TCDF-treated mice (Figure 4I), which may indicate an attempt to compensate for the reduced hepatic taurine content. We also noted a significant elevation in the expression of bile acid-coenzyme A synthetase (Bacs) mRNA--encoding an important enzyme involved in the process of bile acid conjugation with taurine--in the livers of TCDF-treated mice (Figure 4I). These results are in agreement with the observation of elevated taurine-conjugated bile acids.
For enterohepatic circulation of bile acids, both the hepatocyte and the enterocyte must efficiently transport bile acids. Thus, transporters play key roles in the transport and reabsorption of bile acids. In this study, the expression of genes involved in bile acid transport and homeostasis in liver and ileum were determined (Figure 4J-K). The multidrug-resistance protein (Mrp3) gene was significantly up-regulated in the livers of TCDF-treated [Ahr.sup.+/+] mice (Figure 4J). Previous studies reported that TCDD exposure resulted in up-regulation of Mrp3 in mouse liver (Maher et al. 2005) and induced expression of inflammatory mediators as well as necrosis- and/or apoptosis-related genes (Hayashi et al. 2005). Expression of Mrp2 and Ibat mRNAs (encoding apical bile acid transporters) were down-regulated, whereas Mrp3 mRNA (encoding one of the basolateral bile acid transports) was up-regulated in the ileum of TCDF-treated mice (Figure 4K).
Effects of dietary TCDF on the host metabolome. [sup.1]H NMR-based metabolomics coupled with multivariate statistical analysis was employed to evaluate the metabolic changes induced by TCDF exposure. Typical [sup.1]H NMR spectra of liver, fecal, and cecal content extracts obtained from vehicle- and TCDF-treated [Ahr.sup.+/+] mice are presented in Supplemental Material, Figure S5A-F. Metabolite assignments were carried out as described previously (Tian et al. 2012; Wu et al. 2010) and confirmed with a series of 2D NMR experiments (see Supplemental Material, Table S4). To obtain the metabolic variations associated with different biological sample from vehicle- and TCDF-treated mice, we performed pair-wise OPLS-DA of data obtained from feces (Figure 5A), cecal contents (Figure 5B), liver (Figure 5C), and gut tissues, including duodenum, jejunum, ileum, and cecum (see Supplemental Material, Figure S6A-D). The model quality indicators ([R.sup.2]X and [Q.sup.2]; see Supplemental Material, Table S3) clearly showed that the extracts obtained from the biological matrices were distinctive in terms of their metabolite profiles. However, the [Q.sup.2] values (0.007, -0.12, and 0.29 for the models of feces, cecal contents, and liver samples, respectively) of the OPLS-DA models of data from [Ahr.sup.-/-] mice revealed no significant differences between TCDF and vehicle treatment (see Supplemental Material, Figure S6E-G and Table S5). These observations were further supported by results from the model evaluation with analysis of variance of the cross-validated residuals (CV-ANOVA) (p < 0.05) and permutation test for the OPLS-DA models (see Supplemental Material, Figure S7 and Table S5). The metabolites with statistically significant contributions between TCDF-and vehicle-treated groups are identified in the corresponding color-coded coefficient plots (Figure 5A-C; see also Supplemental Material, Figure S6A-D), and the correlation coefficient values were also tabulated (see Supplemental Material, Table S3). Compared with the vehicle-treated mice, dietary TCDF significantly elevated the levels of the short chain fatty acids (SCFAs) propionate and n-butyrate, but significantly decreased the levels of oligosaccharides and glucose in feces and cecal contents (Figure 5A,B; see also Supplemental Material, Table S3). Reliable assignments of n-butyrate and propionate were further confirmed by 2D [sup.1]H-[sup.1]H total correlation spectroscopy (TOCSY) NMR spectroscopy (see Supplemental Material, Figure S8). Similar changes in SCFAs and oligosaccharides were also obtained from calculating their relative concentration from the NMR peaks integration to internal standard (sodium 3-trimethylsilyl [2,2,3,3-d4] propionate; TSP-d4) in feces and cecal contents of the TCDF-treated wild-type mice against those from the respective controls (see Supplemental Material, Figure S9A,B). However, no significant differences in the levels of SCFAs and oligos accharides were observed in feces and cecal content between TCDF-treated and vehicle-treated [Ahr.sup.-/-] mice (see Supplemental Material, Figure S9A,B), supporting the dependence on AHR status. Significant increases in the expression of the G-protein-coupled receptors (Gpr41 and Gpr43) were found in the colons of TCDF-treated mice (see Supplemental Material, Figure S10A).
Dietary TCDF significantly elevated the levels of lipid, unsaturated fatty acid, and glycogen in the livers and significantly decreased the levels of glucose, lactate, fumarate, a range of amino acids, nucleotide metabolites (e.g., inosine and hypoxanthine), nicotinurate, and cell membrane-related metabolites [choline, phosphatidylcholine (PC), and glycerylphosphorylcholine (GPC)] in comparison with the vehicle-treated mice (Figure 5C). However, no significant change in lipid or glucose metabolism was observed in the liver of [Ahr.sup.-/-] mice after TCDF treatment (see Supplemental Material, Figure S6G). Down-regulation of Gck, G6pase, Glut2, and Pepck expression was found in the livers of TCDF-treated [Ahr.sup.+/+] mice (see Supplemental Material, Figure S10B).
In gut tissues, dietary TCDF exposure also resulted in a significant elevation of amino acids, including branched chain amino acids (leucine, isoleucine, and valine), lysine, glutamine, tyrosine, and phenylalanine, but reduction of taurine, choline, PC, GPC, and some nucleic acids such as hypoxanthine and uridine (see Supplemental Material, Figure S6A-D and Table S3).
Historically, POPs have been studied for their ability to cause adverse effects, including liver toxicity (Lee et al. 2014). We anticipate, however, that because most exposure to POPs occurs through the diet, the host gastrointestinal tract and gut microbiota are likely to be exposed to and modulated by POPs. In the present study, TCDF treatment increased gut inflammation, modulated the gut microbiota population, increased bacterial fermentation, and had a profound impact on host metabolism in an AHR-dependent manner.
Based on allometric scaling calculations, the TCDF dose (24 [micro]g/kg body weight) used in this study is equivalent to approximately 3,000 ng/kg body weight in humans. Several epidemiological studies have reported that TCDD blood levels were, on average, 1,434 ng/kg (range 301-3,683 ng/kg) among workers involved in the remediation of the 2,4,5-trichlorophenol reactor accident in Seveso, Italy (World Health Organization 1998). Accounting for the reduced toxicity of TCDF compared with TCDD (the toxic equivalency factor for TCDF is 0.1), the dose used in the present study is considerably higher than the estimated daily oral intake of polychlorinated dibenzo-p-dioxins and dibenzofurans for adults (0.3-3.0 pg TEQ/kg body weight) (Nakatani et al. 2011; Schecter et al. 1994). Furthermore, we have not accounted for total body burden, and as mentioned above, the absorbed TCDF dose is considerably higher than what is expected from normal dietary intakes. However, our study reinforces the idea that gut microbiota are important targets and/or mediators of the toxicologic response and helps to establish possible end points for studies using lower doses of TCDF.
One of the most prominent findings in our investigation was the significant alteration of gut microbiota by TCDF and changes in co-metabolites of the host and gut microbiota, including bile acids and SCFAs. Within the intestinal microbiota, at the phylum level, TCDF treatment significantly shifted the microbial community from Firmicutes toward Bacteroidetes, indicating that the total population of the gut microbiota was modulated by TCDF. Similar observations have been observed with oral exposure to PCBs in mice, which interestingly, could be reversed with voluntary exercise (Choi et al. 2013). The predominant members, Firmicutes and Bacteroidetes, provide key metabolic functions to their host (Hooper et al. 2002), as well as important developmental (Stappenbeck et al. 2002) and immunologic properties (Hooper et al. 2003; Mazmanian et al. 2005). Further, TCDF treatment enriched Butyrivibrio spp., common butyrate-producing gut microbes, that are able to ferment a wide range of sugars and cellodextrins (Russell 1985). The increased levels of Butyrivibrio spp. was confirmed by the activation of fermentation by TCDF treatment, which induced significant depletion of glucose and oligo-saccharides coupled with elevation of SCFAs such as butyrate and propionate in feces and cecal contents. Interestingly, certain members of the class Flavobacteria, which includes Flavobacterium spp., were reported to possess dehalogenase activity, and although it has not been demonstrated in isolates derived from the mammalian gut, they may contribute to the metabolism of TCDF and other halogenated compounds (Xun et al. 1992). A clear role for Oscillibacter spp. has not been established, but has been reported to be associated with the development of fatty liver, an effect observed with TCDD exposure (Angrish et al. 2013; Matsubara et al. 2012). Members of the Clostridia class and Lactobacillus spp. were reported to be associated with strong enzymatic activity of bile salt hydrolase (Li et al. 2013), and the elevation of DCA, [alpha]MCA, and T[beta]MCA suggests that TCDF remodeling of gut microbiota contributes to changes in bile acid pools in the gut and liver (Jin et al. 2014).
After TCDF treatment of mice, SFB were drastically depleted, thus providing further evidence that TCDF impacted the gut microbiota. One of the most outstanding features of SFB as commensal bacteria is their specific attachment to host epithelial cells in the terminal ileum (Caselli et al. 2010). Previous studies demonstrated that SFB were present only in mice that had Th17 cellinducing microbiota, which contain relatively abundant [CD4.sup.+], IFN[[gamma].sup.+], and [CD25.sup.+] [Foxp3.sup.+] cells (Ivanov et al. 2009). SFB lead to an increase in the relative proportions of Th17 cells and Th17 cell cytokines, which provide protection against mucosal infection through modulation of the immune response of the host. Furthermore, other high-affinity ligands of AHR, TCDD and 6-formylindolo [3,2-b] carbazole (FICZ), have been shown to correlate with an increased frequency of Th17 cells (Hautekiet et al. 2008; Marshall and Kerkvliet 2010). Thus, TCDF exposure may be associated with heightened immune response that is dependent on AHR status. Additional studies focused on the direct impact of TCDF and other POPs on the gut microbiota are indeed warranted and should shed some light on how these compounds might influence the gut microbiota population directly.
Supporting the overall assertion that dietary TCDF exposure leads to a broad inflammatory response in mice is the significant elevation of fecal IgA, fecal LCN-2, serum LPS, and mRNA expression of intestinal innate immune factors including Il-1[beta], Il-10, Tnf-[alpha], Saa1, and Saa3 reported here. One possible explanation for these results is that the AHR plays a direct role in the regulation of inflammatory gene expression. There is support for this notion from studies that found functional dioxin response elements in the Il-6 promoter that are capable of synergistically increasing transcription in the presence of an inflammatory signal and an AHR ligand (DiNatale et al. 2010). The binding of the AHR to the Il-6 promoter results in jettison of HDAC1 (histone deacetylase 1) from the promoter, leading to increased NF-[kappa]B (nuclear factor-[kappa]B) occupation and acetylation. Whether the same type of combinatorial regulation occurs on other inflammatory genes remains to be explored. IgA, a major class of immunoglobulin, is produced in mucosal tissues and was reported to reduce intestinal proinflammatory signaling and bacterial epitope expression when binding to the commensal Bacteroides thetaiotaomicron (Peterson et al. 2007). This suggests that IgA inhibits innate immune responses and regulates the composition and the function of gut microbiota (Macpherson et al. 2012). LPS is found in the outer membrane of gram-negative bacteria, contributes greatly to the structural integrity of the bacteria, acts as an endotoxin, and elicits strong immune responses in animals. LPS can mediate an IL-10-dependent inhibition of CD4 T-cell expansion and function by up-regulating levels of PD-1 (programmed cell death protein 1) on monocytes, which leads to IL-10 production by monocytes after binding of PD-1 by PD-L (Said et al. 2010). Previous studies have shown that LPS in the portal vein or translocation of bacteria promotes intestinal fibrogenesis via TLR4 (toll-like receptor 4), predominantly from gram-negative bacteria (Chan et al. 1997; Lin et al. 1995). In the present study, an elevated level of serum LPS would suggest an increase in gut permeability or altered immune surveillance, leading to higher levels of inflammatory signaling in TCDF-treated mice compared with vehicle-treated mice. In addition, levels of LCN-2, a sensitive and dynamic non invasive biomarker for intestinal inflammation, are dramatically increased after TCDF exposure (Chassaing et al. 2012). It is important to note that the observed innate immune responses in mice induced by TCDF required AHR expression. Whether the heightened host inflammatory signaling observed through potent AHR activation leads to changes in the gut microbiota, which in turn causes additional alterations in gut homeostasis, remains to be investigated.
Another prominent finding in the present study was the profound changes of the host metabolome induced by TCDF. In particular, TCDF induced significant reduction of glucose and oligosaccharides with elevation of SCFAs such as n-butyrate and propionate in feces and cecal content extracts, indicating activation of bacterial fermentation by TCDF. SCFAs are clearly one of the most important gut microbial products and affect a range of host processes including energy utilization, host-microbe signaling, gut composition, and motility (Nicholson et al. 2012). Schwiertz et al. (2010) reported that obesity is linked to the composition of human microbiota and the production of SCFAs, providing an additional source of energy for the body. In addition to being energy sources, SCFAs control colonic gene expression by inhibiting the enzyme HDAC and metabolic regulation by signaling through G-protein-coupled receptors, such as GPR41 and GPR43 (Tremaroli and Backhed 2012). Significant elevation of colonic Gpr41 and Gpr43 mRNAs further indicated activation of bacterial fermentation by TCDF.
Significantly elevated levels of lipid and unsaturated fatty acid in the liver of TCDF-treated mice indicate hepatic lipogenesis, which appears to be a common cellular response to a number of liver-toxic compounds such as allyl formate (Yap et al. 2006), acetaminophen (Coen et al. 2003), and aflatoxin-B1 (Zhang et al. 2011). Further, in the present study, depletion of hepatic glucose and lactate coupled with elevation of glycogen in TCDF-treated mice (Figure 5) suggests that one of the consequences of TCDF exposure is inhibition of gluconeogenesis and glycogenolysis. Similar observations were reported in the case of suppression of hepatic gluconeogenesis by TCDD through the AHR target gene TiPARP [TCDD-inducible poly(ADP-ribose) polymerase; Diani-Moore et al. 2010]. We also consistently observed markedly decreased expression and activity of Pepck, Gck, G6pase, and Glut2, all enzymes that control gluconeogenic flux in the liver of TCDF-treated mice. Interestingly TCDF exposure influenced the FXR signaling pathway, but the connection with AHR remains unclear. Further studies using tissue-specific Ahr-knockout mice should yield important insights into how AHR status in different tissues (e.g., liver or intestinal epithelium) impacts the gut microbiota.
Dietary TCDF induced alterations in the gut microbiota accompanied by modulation of bile acid metabolism and host immune responses. Depletion of Clostridia by TCDF was associated with accumulation of DCA, which is involved in liver cancer (Yoshimoto et al. 2013). The increases in TCDCA and T[beta]MCA by TCDF were associated with the inhibition of the FXR signaling pathway. Furthermore, through activation of AHR signaling, dietary TCDF altered many host metabolic pathways involved in hepatic lipogenesis, gluconeogenesis and glycogenolysis, bacterial fermentation, and amino acid and nucleic acid metabolism. These findings provide new evidence that exposure to persistent environmental contaminants strongly impacts the gut microbiota-host metabolic axis in mice.
Alaluusua S, Lukinmaa PL. 2006. Developmental dental toxicity of dioxin and related compounds--a review. Int J Dent 56:323-331.
Angrish MM, Dominici CY, Zacharewski TR. 2013. TCDD-elicited effects on liver, serum, and adipose lipid composition in C57BL/6 mice. Toxicol Sci 131:108-115.
Arsenescu V, Arsenescu RI, King V, Swanson H, Cassis LA. 2008. Polychlorinated biphenyl-77 induces adipocyte differentiation and proinflammatory adipokines and promotes obesity and atherosclerosis. Environ Health Perspect 116:761-768; doi:10.1289/ehp.10554.
Bitoun M, Tappaz M. 2000a. Gene expression of taurine transporter and taurine biosynthetic enzymes in brain of rats with acute or chronic hyperosmotic plasma. A comparative study with gene expression of myo-inositol transporter, betaine transporter and sorbitol biosynthetic enzyme. Brain Res Mol Brain Res 77:10-18.
Bitoun M, Tappaz M. 2000b. Gene expression of taurine transporter and taurine biosynthetic enzymes in hyperosmotic states: a comparative study with the expression of the genes involved in the accumulation of other osmolytes. Adv Exp Med Biol 483:239-248.
Bjeldanes LF, Kim JY, Grose KR, Bartholomew JC, Bradfield CA. 1991. Aromatic hydrocarbon responsiveness-receptor agonists generated from indole-3-carbinol in vitro and in vivo: comparisons with 2,3,7,8-tetrachlorodibenzo-p-dioxin. Proc Natl Acad Sci USA 88:9543-9547.
Caselli M, Holton J, Boldrini P, Vaira D, Calo G. 2010. Morphology of segmented filamentous bacteria and their patterns of contact with the follicle-associated epithelium of the mouse terminal ileum. Gut Microbes 1:367-372.
Chan CC, Hwang SJ, Lee FY, Wang SS, Chang FY, Li CP, et al. 1997. Prognostic value of plasma endotoxin levels in patients with cirrhosis. Scand J Gastroentero 32:942-946.
Chassaing B, Srinivasan G, Delgado MA, Young AN, Gewirtz AT, Vijay-Kumar M. 2012. Fecal lipocalin 2, a sensitive and broadly dynamic non-invasive biomarker for intestinal inflammation. PloS One 7:e44328; doi:10.1371/journal.pone.0044328.
Choi JJ, Eum SY, Rampersaud E, Daunert S, Abreu MT, Toborek M. 2013. Exercise attenuates PCB-induced changes in the mouse gut microbiome. Environ Health Perspect 121:725-730; doi:10.1289/ ehp.1306534.
Claudel T, Staels B, Kuipers F. 2005. The Farnesoid X receptor: a molecular link between bile acid and lipid and glucose metabolism. Arterioscler Thromb Vasc Biol 25:2020-2030.
Claus SP, Tsang TM, Wang Y, Cloarec O, Skordi E, Martin FP, et al. 2008. Systemic multicompartmental effects of the gut microbiome on mouse metabolic phenotypes. Mol Syst Biol 4:219; doi:10.1038/ msb.2008.56.
Coen M, Lenz EM, Nicholson JK, Wilson ID, Pognan F, Lindon JC. 2003. An integrated metabonomic investigation of acetaminophen toxicity in the mouse using NMR spectroscopy. Chem Res Toxicol 16:295-303.
DeVito MJ, Birnbaum LS. 1995. The importance of pharmacokinetics in determining the relative potency of 2,3,7,8-tetrachlorodibenzo-p-dioxin and 2,3,7,8-tetrachlorodibenzofuran. Fundam Appl Toxicol 24:145-148.
Diani-Moore S, Ram P, Li X, Mondal P, Youn DY, Sauve AA, et al. 2010. Identification of the aryl hydrocarbon receptor target gene TiPARP as a mediator of suppression of hepatic gluconeogenesis by 2,3,7,8-tetrachlorodibenzo-p-dioxin and of nicotinamide as a corrective agent for this effect. J Biol Chem 285:38801-38810.
DiNatale BC, Schroeder JC, Francey LJ, Kusnadi A, Perdew GH. 2010. Mechanistic insights into the events that lead to synergistic induction of interleukin 6 transcription upon activation of the aryl hydrocarbon receptor and inflammatory signaling. J Biol Chem 285:24388-24397.
Folch BJ, Lees M, Sloane Stanley GH. 1957. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem 226:497-509.
Greenblum S, Turnbaugh PJ, Borenstein E. 2012. Metagenomic systems biology of the human gut microbiome reveals topological shifts associated with obesity and inflammatory bowel disease. Proc Natl Acad Sci USA 109:594-599.
Hautekiet V, Geert V, Marc V, Rony G. 2008. Development of a sanitary risk index for Salmonella seroprevalence in Belgian pig farms. Prev Vet Med 86:75-92.
Hayashi H, Takada T, Suzuki H, Onuki R, Hofmann AF, Sugiyama Y. 2005. Transport by vesicles of glycine-and taurine-conjugated bile salts and taurolithocholate 3-sulfate: a comparison of human BSEP with rat Bsep. Biochim Biophys Acta 1738:54-62.
Hooper LV, Midtvedt T, Gordon JI. 2002. How hostmicrobial interactions shape the nutrient environment of the mammalian intestine. Annu Rev Nutr 22:283-307.
Hooper LV, Stappenbeck TS, Hong CV, Gordon JI. 2003. Angiogenins: a new class of microbicidal proteins involved in innate immunity. Nat Immunol 4:269-273.
Howell G III, Mangum L. 2011. Exposure to bioaccumulative organochlorine compounds alters adipogenesis, fatty acid uptake, and adipokine production in NIH3T3-L1 cells. Toxicol In Vitro 25:394-402.
Ivanov II, Atarashi K, Manel N, Brodie EL, Shima T, Karaoz U, et al. 2009. Induction of intestinal Th17 cells by segmented filamentous bacteria. Cell 139:485-498.
Jaaskelainen T, Paananen J, Lindstrom J, Eriksson JG, Tuomilehto J, Uusitupa M, et al. 2013. Genetic predisposition to obesity and lifestyle factors--the combined analyses of twenty-six known BMI- and fourteen known waist:hip ratio (WHR)-associated variants in the Finnish Diabetes Prevention Study. Br J Nutr 110:1856-1865.
Jiang C, Xie C, Li F, Zhang L, Nichols RG, Krausz KW, et al. 2015. Intestinal farnesoid X receptor signaling promotes nonalcoholic fatty liver disease. J Clin Invest 125:386-402.
Jin UH, Lee SO, Sridharan G, Lee K, Davidson LA, Jayaraman A, et al. 2014. Microbiome-derived tryptophan metabolites and their aryl hydrocarbon receptor-dependent agonist and antagonist activities. Mol Pharmacol 85:777-788.
Jones KC, Bennett BG. 1989. Human exposure to environmental polychlorinated dibenzo-p-dioxins and dibenzofurans: an exposure commitment assessment for 2,3,7,8-TCDD. Sci Total Environ 78:99-116.
Lee DH, Porta M, Jacobs DR Jr, Vandenberg LN. 2014. Chlorinated persistent organic pollutants, obesity, and type 2 diabetes. Endocr Rev 35:557-601.
Li F, Jiang C, Krausz KW, Li Y, Albert I, Hao H, et al. 2013. Microbiome remodelling leads to inhibition of intestinal farnesoid X receptor signalling and decreased obesity. Nat Commun 4:2384; doi:10.1038/ncomms3384.
Lin RS, Lee FY, Lee SD, Tsai YT, Lin HC, Lu RH, et al. 1995. Endotoxemia in patients with chronic liver diseases: relationship to severity of liver diseases, presence of esophageal varices, and hyperdynamic circulation. J Hepatol 22:165-172.
Macpherson AJ, Geuking MB, McCoy KD. 2012. Homeland security: IgA immunity at the frontiers of the body. Trends Immunol 33:160-167.
Maher JM, Cheng X, Slitt, AL, Dieter MZ, Klaassen CD. 2005. Induction of the multidrug resistance-associated protein family of transporters by chemical activators of receptor-mediated pathways in mouse liver. Drug Metab Dispos 33:956-962.
Marshall NB, Kerkvliet NI. 2010. Dioxin and immune regulation: emerging role of aryl hydrocarbon receptor in the generation of regulatory T cells. Ann NY Acad Sci 1183:25-37.
Mathis D, Benoist C. 2012. The influence of the microbiota on type-1 diabetes: on the threshold of a leap forward in our understanding. Immunol Rev 245:239-249.
Matsubara T, Tanaka N, Krausz KW, Manna SK, Kang DW, Anderson ER, et al. 2012. Metabolomics identifies an inflammatory cascade involved in dioxin- and diet-induced steatohepatitis. Cell Metab 16:634-644.
Maurice CF, Haiser HJ, Turnbaugh PJ. 2013. Xenobiotics shape the physiology and gene expression of the active human gut microbiome. Cell 152:39-50.
Mazmanian SK, Liu CH, Tzianabos AO, Kasper DL. 2005. An immunomodulatory molecule of symbiotic bacteria directs maturation of the host immune system. Cell 122:107-118.
Monteleone I, MacDonald TT, Pallone F, Monteleone G. 2012. The aryl hydrocarbon receptor in inflammatory bowel disease: linking the environment to disease pathogenesis. Curr Opin Gastroenterol 28:310-313.
Monteleone I, Pallone F, Monteleone G. 2013. Aryl hydrocarbon receptor and colitis. Semin Immunopathol 35:671-675.
Myre M, Imbeault P. 2014. Persistent organic pollutants meet adipose tissue hypoxia: does cross-talk contribute to inflammation during obesity? Obes Rev 15:19-28.
Nagalingam NA, Lynch SV. 2012. Role of the microbiota in inflammatory bowel diseases. Inflamm Bowel Dis 18:968-984.
Nakatani T, Yamamoto A, Ogaki S. 2011. A survey of dietary intake of plychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans, and dioxin-like coplanar polychlorinated biphenyls from food during 2000-2002 in Osaka City, Japan. Arch Environ Contam Toxicol 60:543-555.
Nicholson JK, Holmes E, Kinross J, Burcelin R, Gibson G, Jia W, et al. 2012. Host-gut microbiota metabolic interactions. Science 336:1262-1267.
Ozeki J, Uno S, Ogura M, Choi M, Maeda T, Sakurai K, et al. 2011. Aryl hydrocarbon receptor ligand 2,3,7,8-tetrachlorodibenzo-p-dioxin enhances liver damage in bile duct-ligated mice. Toxicology 280:10-17.
Peterson DA, McNulty NP, Guruge JL, Gordon JI. 2007. IgA response to symbiotic bacteria as a mediator of gut homeostasis. Cell Host Microbe 2:328-339.
Ren A, Qiu X, Jin L, Ma J, Li Z, Zhang L, et al. 2011. Association of selected persistent organic pollutants in the placenta with the risk of neural tube defects. Proc Natl Acad Sci USA 108:12770-12775.
Roche Diagnostics. 2013. Short Fragment Removal for the Amplicon Library Preparation Procedure. TCB No. 2011-007. Available: http://454.com/downloads/my454/ documentation/technical-buMetins/TCB-11007_ ShortFragmentRemovalProcedureAmplicon LibraryPrep_May2013.pdf [accessed 12 June 2015].
Russell JB. 1985. Fermentation of cellodextrins by cellulolytic and noncellulolytic rumen bacteria. Appl Environ Microbiol 49:572-576.
Said EA, Dupuy FP, Trautmann L, Zhang Y, Shi Y, El-Far M, et al. 2010. Programmed death-1-induced interleukin-10 production by monocytes impairs CD4+ T cell activation during HIV infection. Nat Med 16:452-459.
Sayin SI, Wahlstrom A, Felin J, Jantti S, Marschall HU, Bamberg K, et al. 2013. Gut microbiota regulates bile acid metabolism by reducing the levels of tauro-beta-muricholic acid, a naturally occurring FXR antagonist. Cell Metab 17:225-235.
Schecter A, Startin J, Wright C, Kelly M, Papke O, Lis A, et al. 1994. Congener-specific levels of dioxins and dibenzofurans in U.S. food and estimated daily dioxin toxic equivalent intake. Environ Health Perspect 102:962-966.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537-7541.
Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, et al. 2010. Microbiota and SCFA in lean and overweight healthy subjects. Obesity (Silver Spring) 18:190-195.
Snedeker SM, Hay AG. 2012. Do interactions between gut ecology and environmental chemicals contribute to obesity and diabetes? Environ Health Perspect 120:332-339; doi:10.1289/ehp.1104204.
Stappenbeck TS, Hooper LV, Gordon JI. 2002. Developmental regulation of intestinal angiogenesis by indigenous microbes via Paneth cells. Proc Natl Acad Sci USA 99:15451-15455.
Tian Y, Zhang L, Wang Y, Tang H. 2012. Age-related topographical metabolic signatures for the rat gastroi ntestinal contents. J Proteome Res 11:1397-1411.
Tremaroli V, Backhed F. 2012. Functional interactions between the gut microbiota and host metabolism. Nature 489:242-249.
World Health Organization. 1998. Executive Summary: Assessment of the Health Risk of Dioxins: Re-Evaluation of the Tolerable Daily Intake (TDI). WHO Consultation, May 25-29 1998, Geneva, Switzerland. Geneva:WHO European Centre for Environment and Health, International Programme on Chemical Safety. Available: http://www.who.int/ ipcs/publications/en/exe-sum-final.pdf [accessed 28 May 2015].
Wu JF, An YP, Yao JW, Wang YL, Tang HR. 2010. An optimised sample preparation method for NMR-based faecal metabonomic analysis. Analyst 135:1023-1030.
Xun L, Topp E, Orser CS. 1992. Glutathione is the reducing agent for the reductive dehalogenation of tetrachloro-p-hydroquinone by extracts from a Flavobacterium sp. Biochem Biophys Res Commun 182:361-366.
Yap IK, Clayton TA, Tang H, Everett JR, Hanton G, Provost JP, et al. 2006. An integrated metabonomic approach to describe temporal metabolic disregulation induced in the rat by the model hepatotoxin allyl formate. J Proteome Res 5:2675-2684.
Yoshimoto S, Loo TM, Atarashi K, Kanda H, Sato S, Oyadomari S, et al. 2013. Obesity-induced gut microbial metabolite promotes liver cancer through senescence secretome. Nature 499:97-101.
Zelante T, Iannitti RG, Cunha C, De Luca A, Giovannini G, Pieraccini G, et al. 2013. Tryptophan catabolites from microbiota engage aryl hydrocarbon receptor and balance mucosal reactivity via interleukin-22. Immunity 39:372-385.
Zhang L, Ye Y, An Y, Tian Y, Wang Y, Tang H. 2011. Systems responses of rats to aflatoxin B1 exposure revealed with metabonomic changes in multiple biological matrices. J Proteome Res 10:614-623.
Limin Zhang, (1,2) Robert G. Nichols, (1) Jared Correll, (1) Iain A. Murray, (1) Naoki Tanaka, (3) Philip B. Smith, (1) Troy D. Hubbard, (1) Aswathy Sebastian, (4) Istvan Albert, (4) Emmanuel Hatzakis, (5) Frank J. Gonzalez, (3) Gary H. Perdew, (1) and Andrew D. Patterson (1)
(1) Center for Molecular Toxicology and Carcinogenesis, Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, Pennsylvania, USA; (2) Chinese Academy of Sciences (CAS) Key Laboratory of Magnetic Resonance in Biological Systems, State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Centre for Magnetic Resonance, Wuhan Institute of Physics and Mathematics, CAS, Wuhan, China; laboratory of Metabolism, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, USA; (4) Bioinformatics consulting Center, and (5) Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania, USA
Address correspondence to A.D. Patterson, Center for Molecular Toxicology and Carcinogenesis, Department of Veterinary and Biomedical Sciences, The Pennsylvania State University, University Park, PA 16801 USA. Telephone: (814) 867-4565. E-mail: firstname.lastname@example.org
Supplemental Material is available online (http:// dx.doi.org/10.1289/ehp.1409055).
This work was supported by the National Institute of Environmental Health Sciences, National Institutes of Health [NIH; grants ES004869 (G.H.P.), ES019964 (G.H.P.), and ES022186 (A.D.P.)] and by the Intramural Research Program of the Center for Cancer Research, National Cancer Institute, NIH. L.Z. is a Penn State Institutes of Energy and the Environment Research Fellow.
The authors declare they have no actual or potential competing financial interests.
Received: 7 August 2014; Accepted: 9 March 2015; Advance Publication: 13 March 2015; Final Publication: 1 July 2015.
Caption: Figure 1. Xenobiotic responses in mice after dietary exposure to vehicle or TCDF (24 [micro]g/kg). Light microscopic examination of H&E-stained liver sections from [Ahr.sup.+/+] (A) and [Ahr.sup.-/-] (B) mice; arrows indicate inflammatory foci (bars = 200 [micro]m). (C) Serum concentrations of ALT (left) and ALP (right) from [Ahr.sup.+/+] and [Ahr.sup.-/-] mice. (D) Percent change in body weight of [Ahr.sup.+/+] mice recorded every other day during the adaptation and treatment periods. (E) Cyp1a1, Cyp1a2, Cyp2e1, and Cyp2a1 expression in the liver (left) and Cyp1a1 mRNA expression in the intestine (right) of [Ahr.sup.+/+] mice. (F) Cyp1a1 and Cyp1a2mRNA expression in the liver of [Ahr.sup.-/-] mice. Data are presented as mean [+ or -] SD; n = 5 or 6/group. NS, not significant.
* p < 0.05, and ** p < 0.01 by two-tailed Student's t-test or Mann-Whitney test.
Caption: Figure 2. Effects of dietary TCDF on the morphology, population, and composition of gut microbiota of mice; animals were treated for 5 days with vehicle or TCDF (24 [micro]g/kg) and sampled on day 7. (4) Weighted UniFrac principal coordinate analysis of the total population of the gut microbiome of cecal content from [Ahr.sup.+/+] and [Ahr.sup.-/-] mice. (B) 16S rRNA gene sequencing analysis of the cecal content of [Ahr.sup.+/+] mice at the phylum level. (C) 16S rRNA gene sequencing analysis of cecal contents of [Ahr.sup.+/+] mice at the class and genus levels. (D) Scanning electron microscopy images of ileum (bars = 50 [micro]m) and (E) qPCR analysis of ileum SFB from [Ahr.sup.+/+] and [Ahr.sup.-/-] mice. Images in (D) represent replicates from two mice in each group. Data are presented as mean [+ or -] SD; n = 5 or 6/group. NS, not significant by two-tailed Student's f-test or Mann-Whitney test.
Caption: Figure 3. AHR-dependent inflammation in mice after dietary exposure to vehicle or TCDF (24 [micro]g/kg). (A) qPCR analysis of inflammatory cytokine (IL-1[beta], TNF-[alpha], IL-10, Saa1, and Saa3) mRNA expression in the ileum of [Ahr.sup.+/+] mice. (B) IL-1[beta] and Tnf-[alpha] expression in the ileum of [Ahr.sup.-/-] mice. (C) qPCR analysis of Lcn-2 mRNA expression in the ileum of [Ah.sup.r+/+] and [Ahr.sup.-/-] mice. (D) Quantification of fecal LCN2 in [Ahr.sup.+/+] mice by ELISA. (E) qPCR analysis of Myosin Vb and Ptprh mRNA in the ileum of [Ahr.sup.+/+] mice after TCDF treatment. (F) Quantification of serum LPS in [Ahr.sup.+/+] mice. (G) Quantification of IgA in [Ahr.sup.+/+] mice by ELISA. Data are presented as mean [+ or -] SD; n = 5/group. NS, not significant.
* p < 0.05, ** p < 0.01, and ***p < 0.001, by two-tailed Student's t-test.
Caption: Figure 4. Bile acid metabolism in animals treated with vehicle or TCDF (24 ([micro]g/kg). Abbreviations: CA, cholic acid; CDCA, chenodeoxycholic acid; DCA, deoxycholic acid; G, glycine-conjugated; LCA, lithocholic acid; MCA, muricholic acid; NS, not significant; T, taurine-conjugated species; UDCA, ursodeoxycholic acid. Quantification of specific bile acids levels in intestinal tissue (A,B) and feces (C,D) throughout the enterohepatic circulation of [Ahr.sup.+/+] mice. (E) Quantification of total bile acids in liver, intestine, cecum, and feces of [Ahr.sup.+/+] mice; the bile acid profile in the small intestine shows the data from jejunum segment. qPCR analysis of Fgf15, Fxr, and Shp mRNAs in the ileum (F) and Fxrand Shp mRNA expression in the liver (G) of [Ahr.sup.+/+] and [Ahr.sup.-/-] mice. (H) qPCR analysis of Cyp7a1, Cyp8b1, Cyp27a1, Akr1d1, and Cyp7b1 mRNAs in the liver of [Ahr.sup.+/+] mice. (I-K) mRNA encoding bile acid transporters involved in taurine biosynthesis and bile acid conjugation in the ileum (I), and mRNA encoding bile acid transporters in the distal liver (J) and ileum (K) of vehicle- and TCDF-treated [Ahr.sup.+/+] mice. Data are presented as mean [+ or -] SD; n = 6/group. See also Supplemental Material, Tables S1 and S2.
*p < 0.05, and **p < 0.01, by two-tailed Student's f-test.
Caption: Figure 5. Host metabolism and the bacterial fermentation process in [Ahr.sup.+/+] mice after dietary exposure to vehicle or TCDF (24 [micro]g/kg). OPLS-DA scores (left) and coefficient-coded loadings plots for the models (right) from NMR spectra of aqueous fecal (A), cecal content (6), and liver extracts (C), discriminating between the vehicle (black circles) and TCDF-treated mice (red squares). These models are cross-validated with CV-ANOVA: p = 1.64 x [10.sup.-3], p = 0.033, and p = 0.0018 for feces, cecal content, and liver, respectively. Metabolite assignment is shown in Supplemental Material, Figure S5 and Table S4, and correlation coefficient values for the significantly changed metabolites are shown in Supplemental Material, Table S3.
Please note: Illustration(s) are not available due to copyright restrictions.
|Printer friendly Cite/link Email Feedback|
|Author:||Zhang, Limin; Nichols, Robert G.; Correll, Jared; Murray, Iain A.; Tanaka, Naoki; Smith, Philip B.;|
|Publication:||Environmental Health Perspectives|
|Date:||Jul 1, 2015|
|Previous Article:||A case-only study of vulnerability to heat wave-related mortality in New York City (2000-2011).|
|Next Article:||Understanding the radioactive ingrowth and decay of naturally occurring radioactive materials in the environment: an analysis of produced fluids from...|