# Identification of the novel role of butyrate as AhR ligand in human intestinal epithelial cells

Ludovica Marinelli, Camille Martin-Gallausiaux, Jean-Marie Bourhis, Fabienne Beguet-Crespel, Herve Blottiere, Nicolas Lapaque

► To cite this version:

Ludovica Marinelli, Camille Martin-Gallausiaux, Jean-Marie Bourhis, Fabienne Beguet-Crespel, Herve Blottiere, et al. Identification of the novel role of butyrate as AhR ligand in human intestinal epithelial cells. Scientific Reports, 2019, 9 (1), pp.643. 10.1038/s41598-018-37019-2. hal-01996381

HAL Id: hal-01996381

https://hal.sorbonne-universite.fr/hal-01996381

Submitted on 28 Jan 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Identification of the novel role of butyrate as AhR ligand in human intestinal epithelial cells

Received: 26 June 2018

Accepted: 28 November 2018

Published online: 24 January 2019

Ludovica Marinelli<sup>1,2</sup>, Camille Martin-Gallausiaux<sup>1,2</sup>, Jean-Marie Bourhis<sup>1,3</sup>, Fabienne Béguet-Crespel<sup>1</sup>, Hervé M. Blottière<sup>1,4</sup> & Nicolas Lapaque<sup>1</sup>

The ligand activated transcription factor, aryl hydrocarbon receptor (AhR) emerged as a critical regulator of immune and metabolic processes in the gastrointestinal tract. In the gut, a main source of AhR ligands derives from commensal bacteria. However, many of the reported microbiota-derived ligands have been restricted to indolyl metabolites. Here, by screening commensal bacteria supernatants on an AhR reporter system expressed in human intestinal epithelial cell line (IEC), we found that the short chain fatty acid (SCFA) butyrate induced AhR activity and the transcription of AhR-dependent genes in IECs. We showed that AhR ligand antagonists reduced the effects of butyrate on IEC suggesting that butyrate could act as a ligand of AhR, which was supported by the nuclear translocation of AhR induced by butyrate and *in silico* structural modelling. In conclusion, our findings suggest that (i) butyrate activates AhR pathway and AhR-dependent genes in human intestinal epithelial cell-lines (ii) butyrate is a potential ligand for AhR which is an original mechanism of gene regulation by SCFA.

The mammalian gastrointestinal (GI) tract is colonised by a complex microbial community, referred as gut microbiota. It is well established that host-commensal bacteria crosstalk provides numerous functions for the overall host wellbeing, through the production of microbial metabolites. The host-microbiota interaction is particularly substantial for mucosal barrier functions as well as the development and maintenance of the mucosal immune system<fnref n="1" />. Metabolites derived from commensal bacteria described to strongly impact mucosal homeostasis include the short-chain fatty acids (SCFA), originated from the fermentation of dietary fibres, and indoles, from the degradation of dietary tryptophan. The targeted host receptors of these bacterial products include the cell-surface G-protein-coupled receptors GPR41, GPR43, and GPR109A and nuclear receptors, such as the aryl hydrocarbon receptor (AhR), the pregnane X receptor (PXR), and the farnesoid X receptor (FXR). The receptor-metabolite interactions induce signalling pathways that modulate host gene expression and collectively impact on host metabolism and immune responses<fnref n="1" />.

Recently AhR, a ligand activated transcription factor, has gained considerable attention as a crucial modulator of mucosal immune and metabolic processes, especially in the context of diet and microbiota crosstalk with the host²⁻⁴. AhR is a member of the basic helix-loop-helix PER-ARNT-SIM (bHLH-PAS) family, initially identified as a hepatic intracellular protein that bounds with high affinity the environmental halogenated contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, extensive studies revealed that AhR interacts with a wide range of structurally diverse molecules that originate from the diet, the environment, the microbiota or are endogenously produced by the host itself⁵,⁶. Many of the microbiota-derived AhR ligands result from tryptophan catabolism including indole, indole-3-acetic acid, indole-3-aldehyde⁷,⁸. Upon ligand-binding, cytoplasmic AhR translocates in the nucleus, dimerizes with AhR nuclear translocator (ARNT) and initiates the transcription of target genes with promoters containing a xenobiotic-response element (XRE) sequence, such as cytochrome P450 family 1A1 (CYP1A1) and the repressor AHRR⁹,¹⁰.

AhR activation has been largely reported to be implicated in colonic stem cells proliferation, epithelial barrier functions and the regulation of a wide range of immune cell populations including intraepithelial RORγt⁺ innate lymphoid cell-3 (ILC3), γδT cells, T helper (Th)17/22 cells, intraepithelial, FOXP3⁺ regulatory T cells and antigen presenting cells²⁻⁴,¹¹,¹². The AhR-dependent IL22 expression by ILC3 regulates the release of antimicrobial

<footnote n="1">¹Micalis Institute, INRA, AgroParisTech, Université Paris-Saclay, 78350, Jouy-en-Josas, France. ²Sorbonne Universités, UPMC Univ Paris 06, IFD, 4 place Jussieu, 75252, Paris, cedex 05, France. ³University Grenoble Alpes, CNRS, CEA, IBS, F-38000, Grenoble, France. ⁴MetaGenoPolis, INRA, Université Paris-Saclay, 78350, Jouy en Josas, France. Correspondence and requests for materials should be addressed to N.L. (email: nicolas.lapaque@inra.fr)</footnote>
peptides such as RegIIIγ and the expression of tight junctions molecules in IECs thus reinforcing the barrier and defence functions¹³,¹⁴. Interestingly, AhR signalling and known AhR ligands are low in inflammatory bowel diseases (IBD) patients, highlighting the clinical relevance of the AhR pathway in these pathologies¹⁵,¹⁶. AhR activation by ligand administration and AhR knock-down in mouse suppressed and enhanced DSS-induced colitis, respectively suggesting a beneficial effect of AhR activation in IBD¹⁷,¹⁸. An increasing number of bacterial metabolites have been shown to drive the AhR activation, with described protective effects against intestinal inflammation and pathogens colonisations, suggesting a possible role of this signalling pathway in the intestinal homeostasis⁷. This hypothesis has been further supported by a recent study showing that one IBD-associated single nucleotide polymorphism (SNP) within the *CARD9* gene, affects microbiota composition thus altering the production of bacterial AhR ligands and consequently intestinal inflammation¹⁵.

The relevant role of AhR in the maintenance of mucosal homeostasis has been largely focused on immune cells from the *lamina propria*, and the impact of AhR activation in intestinal epithelial cells (IEC) has only been starting to be unravelled. Emerging evidences highlight that AhR activation impacts also IECs responses contributing to intestinal homeostasis. AhR activation in IECs protects mice from colitis by enhancing IECs barrier functions via the increase of the IL10 receptor expression and the enhancement of tight junctions integrity through the regulation of *Notch1*¹⁹,²⁰. Importantly, the excess of AhR ligands in the gut is detrimental to the host. In this context, a recent study in mice model nicely demonstrates that IECs play the role of gatekeeper *via* the expression of the AhR-regulated gene *Cyp1a1* encoding cytochrome P450 1A1 involved in AhR ligands clearance²¹,²². These studies highlighted the importance of AhR ligands in modulating host gut immune homeostasis and prompted us to identify new microbiota-derived activators of the AhR pathway in IECs. We thus tested the bacterial supernatants of over 100 bacterial species of the human microbiota on an AhR reporter system in human intestinal cell lines and found that butyrate-producing bacteria activate the AhR-dependent pathway. A recent work previously highlighted the ability of butyrate to enhance the expression of AhR-dependent genes through its histone deacetylase inhibitor (HDACi) properties²³. We confirmed the butyrate-activating role at the transcriptional level on AhR-dependent genes in Caco-2 and HT-29 cell lines. Interestingly, other HDACi that could enhance AhR-dependent gene expression were unable to mimic the butyrate dependent-activation of the AhR reporter system, suggesting the existence of a second mechanism. We showed that (i) AhR ligand antagonists impaired the butyrate-induced activation of AhR reporter system, (ii) butyrate enhanced AhR nuclear translocation and (iii) *in silico* modelling of butyrate interaction with human AhR, highlighting for the first time that butyrate could act as a ligand of AhR.

## Results

**Metabolites derived from commensal bacteria enhanced AhR activity.** In the gut, Aryl hydrocarbon receptor (AhR) ligands derive from diverse origins that include the intestinal microbiota as one of the main sources. To decipher which bacteria from human gut activate the AhR pathway, we performed a screening of commensal bacteria on a human intestinal epithelial cell line (HT-29-AhR) stably expressing an AhR-dependant reporter system. The AhR reporter system contains three copies of the DNA-binding domain motif recognized by AhR (xenobiotic-responsive element; XRE) driving the transcription of the luciferase reporter gene. Since AhR activators produced by bacteria are secreted in the intestinal lumen, we performed the screening of bacterial culture supernatants⁸. The screening included 132 bacterial strains, belonging to the major phyla of the human intestinal microbiota (Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria and Verrucomicrobia) (Supplementary Fig. S1, Supplementary Table S1) and grown under appropriate bacterial growth conditions (Supplementary Table S2). When possible, we cultured the bacteria in different media, to avoid biased results due to the rich composition of culture media, at least for almost all the bacteria activating the AhR reporter system. In our experimental set-up, AhR activation was detected in HT-29 cells challenged with some supernatants derived from Proteobacteria, Firmicutes, Fusobacteria and few Actinobacteria (Fig. 1).

**Butyrate activates AhR pathway in intestinal epithelial cell lines.** Despite the huge literature on AhR ligands, only few are commensal-derived molecules. Amongst them, indole and other tryptophan derivatives are produced by a wide range of bacteria including *Lactobacillus* (Firmicutes) and *Proteobacteria*¹⁴,²⁴. In our screening we confirmed the AhR activation induced by some lactobacilli as well as by members of the Proteobacteria, probably due to the production of indole derivatives. Interestingly, among the bacteria not reported or not predicted to be indole producers, we identified some members of genera *Ruminoclostridium* and *Roseburia*, belonging to the Firmicutes phylum. These genera, together with other AhR-activating bacteria predicted to produce indole (*Clostridium* and *Lachnoclostridium*) share the ability to degrade diet-derived fibres, leading to the production of short-chain fatty acids (SCFAs)²⁵. We thus hypothesized that SCFAs concentration in the supernatants of these bacteria could explain the activation of the AhR pathway as previously reported²³. We therefore quantified the concentration of SCFAs in some bacterial supernatants (Supplementary Table S1 and Fig. 2A). A principal component analysis (PCA) of the complete dataset (SCFA concentrations and AhR activity) revealed that AhR activity and butyrate concentration were major contributors of total variability (Fig. 2B). Further Spearman correlation analysis showed a strong positive correlation between AhR activation and the production of butyrate quantified in commensal bacteria supernatants (rho = 0.4966, Fig. 2C).

To reinforce experimentally the observed correlations, we tested pure SCFAs at different physiological concentrations found in the intestine (ranging from 0.125 mM to 8 mM) on HT-29-AhR reporter cells (Fig. 2D,F)²⁶,²⁷. Among the tested SCFAs, we confirmed that butyrate was able to induce AhR activation in a dose-dependent manner, in HT-29-AhR cell line as well as in another intestinal reporter cell line, Caco2-AhR (Fig. 2D,E). Propionate, another abundant SCFA produced by bacteria, was also able to activate AhR in a dose-dependent manner in HT-29-AhR cells at a concentration starting from 1 mM (Fig. 2F). Interestingly, we also observed that valerate and the branched chain fatty acids iso-valerate, both described to influence epithelial physiology, were
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Screening of in-house strain collection of commensal bacteria on HT-29-AhR cell line. HT-29-AhR reporter cells were incubated with bacterial supernatants or relative non-inoculated bacterial media for 24 h (10% vol/vol). AhR activation was measured by luciferase activity and expressed as fold increase (±SD from triplicates) toward its control (non-inoculated bacterial media), sorted by Phyla.</figcaption></figure>

also able to activate AhR pathway in HT-29-AhR cell line at similar concentrations²⁸,²⁹ Iso-butyrate activated AhR only at 8 mM a concentration rarely reached in the intestine while acetate, the most abundant SCFA produced by commensal bacteria showed activation at high concentration (20 mM, Fig. 2F and Supplementary Fig. S2). Apart for acetate, the SCFA concentrations inducing AhR activity were consistent with the final SCFAs concentrations on bacterial supernatants measured in the screen thus reinforcing our hypothesis (Supplementary Table S1).

Moreover, we showed by qRT-PCR that AhR-regulated genes, *CYP1A1*, *AHR* and *AHRR*, were up-regulated by butyrate both in HT-29 and Caco2 cells confirming that this SCFA activated AhR pathway (Fig. 3A,B). We confirmed the butyrate-driven up-regulation of AHR at the protein level by Western-Blot analysis in HT-29 (Fig. 3C).

**AhR activation by butyrate is independent of the SCFA receptors GPR41, GPR43, GPR109a and the SCFAs transporter MCT-1.** Butyrate, like other SCFAs, activates eukaryotic cells through two main mechanisms: activation of specific G-protein coupled receptors (GPCR: GPR41, GPR43, GPR109a) and inhibition of histone deacetylases (HDAC)³⁰⁻³⁴. Considering that the three G-protein coupled receptors are expressed in HT-29 and Caco-2 cells (Martin-Gallausiaux *et al. in press*), we firstly tested if the butyrate-induced AhR activation could be mediated by GPRs signalling by stimulating HT-29-AhR and Caco2-AhR cells with known GPRs agonists, targeting GPR41, GPR43 and GPR109a. For each GPR, two agonists were tested (GPR41: 4-CMTB and Tiglic acid; GPR43: AR420626 and 1-MCPC; GPR109a: Niacine and MK1903). Interestingly, none of the tested agonists induced AhR activation in luciferase reporter system suggesting that these GPRs are not involved in the butyrate-induced activation in HT-29-AhR and Caco-2-AhR cell lines (Fig. 4A,B). GPR41, GPR109a are both Goi coupled receptors whereas GPR43 is a Goi and Goq coupled receptor. To further confirm our observation, we used an inhibitor of the Goi pathway inhibitor: the pertussis toxin (Ptx) in HT-29-AhR (Fig. 4C). No impact on the butyrate-induced AhR activation was detected in cells when the Goi subunit was inhibited, further confirming that these GPRs were not involved in the observed AhR activation.

Butyrate uptake by the monocarboxylate transporter MCT-1 has been described as crucial for the GPR-independent regulation of a wide range of genes by butyrate³⁵. To assess if MCT-1 participated in the activation of AhR by butyrate, a well-described inhibitor of MCT-1 transporter, pCMB, was tested on HT-29-AhR cells (Fig. 4D). The inhibition of MCT-1 transporter did not affect the ability of butyrate to activate AhR reporter system, suggesting that MCT-1 was not involved in this process.

**Activation of the AhR signalling pathway by butyrate is independent of its role as inhibitor of HDAC.** SCFAs, *via* their ability to inhibit lysine and histone deacetylases (HDAC), are potent modulators of histones and transcription factors acetylation that are well-documented regulatory mechanisms of gene transcription³⁰⁻³². A recent study showed that SCFAs and other HDAC inhibitors (HDACi) enhanced the expression of *AHR*, *AHRR* and *CYP1A1* *via* the increased level of histone acetylation of AhR-dependent genes²⁵. These results prompted us to investigate the direct role of HDACi in the activation of AhR signalling pathway. To assess if butyrate impacted AhR pathway through its HDACi property, we tested three HDACi targeting a wide range of HDAC, trichostatin A (TSA), Vorinostat (SAHA) belonging to the hydroxamic acids family, structurally and
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Figure 2.** Butyrate activates AhR pathway. (A) HPLC quantification of butyrate produced by tested bacterial supernatants, sorted by Phyla (N.D.: Not Determined). (B) Principal Component Analysis (PCA) of quantified SCFAs and AhR activation for the entire data set. The axes PC1 and PC2 accounted for 34,83% and 20.75% of total data variation, respectively. (C) Spearman correlation analysis showing a positive relationship (rho = 0.4966) between AhR activity and butyrate concentration produced in commensal supernatants classified by rank values. HT29-AhR (D) or Caco2-AhR (E) reporter cells were incubated with a range of concentration of butyrate (0.125 mM to 8 mM). (F) HT-29-AhR reporter cells were incubated with SCFAs at concentrations rising from 0.125 mM to 8 mM. Data are expressed as luciferase fold (±SEM) of at least three independent experiments, normalised on untreated cells. ns: P > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Butyrate activates the expression of AhR-regulated genes. The expression of AhR-regulated genes: AHR, *CYP1A1* and *AHRR* on HT-29 cells (A) and Caco-2 cells (B) treated with butyrate 2 mM or TCDD 10 nM during 6 h was determined by qRT-PCR. *AHR*, *CYP1A1* and *AHRR* relative expression to control is determined by the 2^−ΔΔCt^ method using *β-actin* for normalisation. (C) Up-regulation of AhR protein level by butyrate. HT-29 cells were incubated with butyrate 2 mM or TCDD 10 nM for 24 h. Total cytoplasmic extracts were blotted (Western Blot) for AhR (left panel) GAPDH was used as loading control. Relative quantifications from three independent experiments are expressed as fold-change to un-stimulated control of AhR protein normalised on GAPDH level (right panel). Full-length blots are presented in Supplementary Fig. S4. Data are expressed as means ± SEM of at least three distinct experiments, performed in triplicate. ns: *P* > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.</figcaption></figure>

metabolically unrelated to SCFAs and sodium valproate (VAP) belonging to the fatty acid family³⁶. Interestingly, TSA, SAHA and VAP did not reproduce the butyrate-induced AhR activation in both HT-29-AhR and Caco2-AhR reporter system, suggesting that an additional mechanism leading to butyrate activation of the AhR pathway might exist (Fig. 5A,B). However, we could reproduce the activation of *CYP1A1* by TSA as observed by Jin *et al*. (Fig. 5C)²³. Altogether, these results suggested that HDACi property of SCFAs was not involved in the butyrate-dependent activation of the AhR signalling pathway as monitored with reporter systems although HDAC inhibition might participate in the regulation of AhR-induced genes.

**Butyrate acts as an AhR ligand.** We showed that the two main mechanisms described for butyrate to activate transcription factors, HDACi and the activation pathway via GPRs, are not directly involved in the butyrate-induced activation of AhR signalling pathway monitored by reporter cell assay. Two groups described an original mechanism of gene regulation where butyrate acts as a ligand of the transcription factor PPARγ³⁷,³⁸. As AhR is a ligand binding transcription factors, we hypothesized that butyrate could act as an AhR ligand similarly to what was described for PPARγ. To demonstrate this hypothesis, we incubated butyrate-stimulated HT-29-AhR and Caco-2-AhR reporter cell-lines with two well-characterized antagonists of the AhR-ligand binding (CH-223191 and GNF-351) and one HSP90 inhibitor, described to avoid the downstream AhR-XRE binding, by blocking the release of HSP90 chaperon from the AhR complex (epigallocatechine-3-gallate, EGCG)³⁹⁻⁴². The three tested inhibitors significantly decreased the butyrate-induced activation of AhR signalling pathway in both HT-29 and Caco-2 reporter cell lines (Fig. 6A,B). In addition, we confirmed in HT-29 cells that CH-223191 and GNF-351 were able to block *CYP1A1* expression induced by butyrate (Fig. 6C). The pivotal role of AhR in *CYP1A1* upregulation by butyrate was assayed using siRNA. As shown in Fig. 6D, we observed a diminished *CYP1A1* activation by butyrate and TCDD in cells treated with AhR siRNA compared to control siRNA. These results suggested that the activation of AhR signalling pathway by butyrate was due to a ligand-dependent AhR-XRE interaction, indicating a possible role of butyrate as a direct AhR ligand.
<figure><img src="image_11.png" /><figcaption>**Figure 4.** Mechanism involved in AhR activation by butyrate is independent of GPR and MCT-1. HT-29 (A) and Caco-2 (B) AhR reporter cells were stimulated for 24 h with GPR agonists. GPR41 (orange): 1-MCPC (1 mM), AR-420626 (1 µM); GPR43 (blue): Tiglic acid (1 mM), MK1903 (1 µM); GPR109a (green): Niacine (1 mM), 4-CMTB (1 µM). (C) HT-29-AhR reporter cells were incubated with $G\alpha_i$-subunit inhibitor, Pertussis toxin (Ptx, 0.2 µg/mL) prior a 24 h butyrate stimulation (2 mM). Vehicle (Glycerol) was used as control. (D) HT29-AhR reporter cells were incubated during 24 h with a MCT-1 inhibitor, *p*-chloromercuribenzoic acid (pCMB, 10 uM), prior stimulation with butyrate (2 mM). AhR activation was measured by luciferase activity and expressed as fold increase means (±SEM) of at least three independent experiments, normalised on untreated cells. ns: P > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.</figcaption></figure>

It is well established that ligand binding triggers the accessibility of the nuclear localisation signal on the AhR N-terminus that consequently initiates the AhR nuclear translocation event$^{43}$. Thus, additional evidence for the human AhR agonist potential of butyrate was obtained by performing nuclear translocation assay. Sub-cellular localisation of AhR in HT-29 cells incubated with butyrate (1, 3 and 6 h) and TCDD (1 h) was assessed by immunoblotting assay. We showed an accumulation of AhR proteins in the nucleus upon treatment with butyrate starting at 3 h and with a peak at 6 h, consistent with the action of an AhR agonist such as TCDD (Fig. 7).

Altogether, our results suggested that butyrate activation of the AhR reporter system relied on AhR translocation to the nucleus and consequent DNA-binding and that butyrate might act as a ligand of AhR.

**In silico modelling of butyrate interaction with AhR.** Numerous AhR modulators, binding in the large central pocket of the AhR PAS-B domain, have been extensively studied using a combination of homology modelling and docking simulations (for review see$^{44}$). To support our results suggesting that butyrate acts as an AhR ligand, we performed a molecular docking simulation of the butyrate/AhR binding process. We first generated a homology model of human AhR ligand binding domain (AhR PAS-B) using HIF-2α PAS-B domain as template similarly to Bisson et al.$^{45,46}$. The AhR PAS-B homology domain was obtained with RaptorX
<figure><img src="image_12.png" /><figcaption>**Figure 5.** Impact of HDAC inhibitors on AhR reporter system and *CYP1A1*. HT29 (A) and Caco-2 (B) -AhR reporter cells were stimulated with HDAC inhibitors for 24 h. Trichostatin A (TSA, 0.1 µM and 1 µM), Vorinostat (SAHA 5 µM), Valproic acid (VAP 5 mM), butyrate (But, 2 mM) and TCDD (10 nM). (C) Relative expression to control of *CYP1A1* on HT-29 cells treated with butyrate 2 mM, TCDD 10 nM or TSA 0.1 µM during 6 h. *CYP1A1* expression induced by drugs is expressed as relative to control expression and is determined by the 2^(-ΔΔCt)^ method using *β-ACTIN* as control gene. AhR activation was measured by luciferase activity and expressed as fold increase means (±SEM) of at least three independent experiments, normalised on un-treated cells. ns: P > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.</figcaption></figure>

modelling software and the binding cavity was defined using HOLLOW and Carver softwares (Fig. 8A). Then, the model was minimized to avoid clashes between side chains. From the best scoring solutions cluster using HADDOCK, we observed different orientations of the butyrate inside the AhR PAS-B cavity. Based on four different docking experiments, two main orientations were found for butyrate in the pocket of human AhR PAS-B domain. Interactions between butyrate and AhR are dominated by polar contacts with the side chains of Q383 and S365, for the first orientation (Fig. 8B) and with side chains of Q383 and H291, for the second one (Fig. 8C). Interestingly, by comparing our results with published docking analyses of other ligands on human AhR PAS-B, the first proposed orientation for butyrate docking shared the same interacting side chains (Q383 and S365) with the docking of FICZ, proposed by Bisson *et al.*⁴⁵. Overall these molecular docking analyses are coherent with our findings of butyrate as a direct modulator of AhR by ligand binding.

Discussion

AhR is a ligand-activated transcription factor that is crucial for intestinal homeostasis by repressing inflammation and by maintaining the epithelial barrier in the gastrointestinal tract³,⁴⁷. The gastrointestinal tract is a rich source of AhR ligands, which have been shown to induce AhR dependent responses and to protect the gut upon infection or induced colitis. Despite the huge literature on the structurally pleiotropic nature of AhR ligands,
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 6.** AhR antagonists inhibit AhR activation by butyrate. HT-29-AhR reporter cells (A) and Caco-2-AhR reporter cells (B) were incubated with AhR ligand antagonists (CH-223191 1 µM and GNF-351 1 µM) and an HSP90 inhibitor ((−)-Epigallocatechin-3-gallate, EGCG, 20 µM) prior stimulation with butyrate (But, 2 mM) and TCDD (10 nM) for 24 h. (C) Relative expression of *CYP1A1* on HT-29 cells treated during 6 h with 2 mM butyrate or TCDD 10 nM in presence or absence of AhR antagonists. *CYP1A1* expression is normalised to control expression (unstimulated cells) and is determined by the 2^−ΔΔCt^ method using *β-actin* as control gene. (D) Relative expression of *CYP1A1* on HT-29 cells transfected with control siRNA (white bars) or AhR siRNA (grey bars) and incubated during 6 h with 2 mM butyrate or TCDD 10 nM. *CYP1A1* expression is normalised to stimulated control expression (cells transfected with control siRNA and stimulated with either butyrate or TCDD), expressed in percentage and is determined by the 2^−ΔΔCt^ method using *β-actin* as control gene. Data are means ± SEM of at least three distinct experiments, performed in triplicate. ns: P > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 7.** Butyrate stimulates AhR protein translocation in the nucleus. HT-29 cells were incubated during 1h, 3h or 6h in presence of 2 mM butyrate or 1h with TCDD. Nuclear and cytoplasmic extracts were blotted (Western Blot) for AhR protein expression. GAPDH and Lamp were used as control for cytoplasmic extract preparation; Lamin A/C as control for nuclear extract preparation (left panel). Relative quantification of AhR protein level in nuclear extracts from least six independent experiments is expressed as fold-change (±SEM) to un-stimulated cells (control) after normalisation on Lamin A/C level. Full-length blots are presented in Supplementary Fig. S5. ns: P > 0.05, \*P ≤ 0.05, \*\*P ≤ 0.01, \*\*\*P ≤ 0.001, \*\*\*\*P < 0.0001.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 8.** Modelling of butyrate binding to the ligand-binding pocket of human AhR. (A) Model of PAS-B domain of human AhR represented in cartoon coloured in wheat, the central cavity is represented as a red surface. Model of PAS-B domain of human AhR using 4ZP4 Crystal Structure of the Heterodimeric HIF-2α:ARNT complex obtained using RaptorX. The cavity has been defined by CARVER and HOLLOW. (B) Q383/S365 orientation. The residues Q383 and S365 are displayed as sticks and coloured by atom type with carbon in wheat. Butyrate is displayed as sticks and coloured by atom type with carbon in magenta. Hydrogen bonds are represented by white dashed lines. (C) Q383.H291 orientation. The residues Q383 and H291 are displayed as sticks and coloured by atom type with carbon in wheat. Butyrate is displayed as sticks and coloured by atom type with carbon in green. Hydrogen bonds are represented by white dashed lines. The figures were generated by PyMol.</figcaption></figure>
only few are commensal-derived molecules. Amongst them, indoles and other tryptophan derivatives are pro-
duced by a variety of bacteria including some Lactobacilli (Firmicutes) and Proteobacteria¹⁴,²⁴. Interestingly,
non-commensal bacteria have been reported to produce non-indole AhR ligands such as the phenazine deriv-
ative from *Mycobacteria* and the 1,4-dihydroxy-2-naphthoic acid from the probiotic *Propionibacterium freud-
enreichii*¹³,⁴⁸. Given the abundance and extensive metabolic capacity of gut microbiome, it is thus likely that
metabolites apart from tryptophan derivatives are present in high concentration to stimulate AhR pathway.

By screening commensal bacteria supernatants, we identified butyrate as a potent activator of AhR pathway and AhR-dependent genes. Butyrate is a short chain fatty acid (SCFA) that derived from bacterial fermentation of dietary fibers. Butyrate, via its role as inhibitor of lysine/histone deacetylase (K/HDACi), influences expression of a large variety of host genes in the colon, including some encoding for immune proteins⁴⁹. Previous studies described butyrate as a regulator of AhR-dependent genes thought its role as histone deacetylase inhibitors (HDACi)²³,⁵⁰. However, these studies assessed the butyrate effect through induction of AhR transcriptional targets and not directly on the AhR receptor. Our results are in accordance with these studies as we showed that HDACi and butyrate upregulate *CYP1A1* expression. However, we demonstrated that butyrate does not impact AhR-dependent gene expression solely by its HDACi properties. Indeed by using an AhR reporter system in HT-29 and Caco2 cells, we showed that butyrate, and not HDAC inhibitors, activates AhR signalling pathway.

To decipher the mechanism of butyrate-dependent AhR activation, we investigated the implication of butyrate specific G-protein coupled receptors (GPR41, GPR109a and GPR43) and main transporter, monocarboxylate transporter 1 (MCT1)⁵¹⁻⁵³. By using agonists and inhibitors, we showed that AhR activation by butyrate was not mediated by those GPRs or MCT1. Many studies evidenced that the predominant biological activities of the AHR are through ligand binding. Despite the AhR association with xenobiotic compounds, structurally diverse metabolites from the diet, bacteria or produced by the host have been reported as capable of binding to human AhR⁶,⁵⁴. Our experimental results suggest that butyrate acts as an AhR ligand leading to the activation of the AhR pathway. Moreover, structural modelling of the binding of butyrate to human AhR PAS-B supports this hypothesis as we observed two main orientations of butyrate in the ligand-binding cavity of AhR. Interestingly, one of these orientations showed polar contact of butyrate with the same side chains of PAS-B reported for the binding of FICZ to human AhR supporting our hypothesis that butyrate is a ligand of AhR⁴⁵. Confirmation of these results by direct amino acids point mutations in the ligand-binding cavity of human AHR similarly to previous publication on mice Ahr would have support our hypothesis⁵⁵. This is an original mechanism of butyrate-dependent modulation of host gene expression, which has only been reported for another ligand-dependent transcription factor, PPARγ. Indeed, two groups reported that PPARγ-dependent genes regulation by butyrate is independent of its HDACi properties and is mediated through a direct binding of the SFCA to this transcription factor³⁷,³⁸. The strong impact of butyrate on mucosal immune homeostasis has been largely documented in mice and in vitro models. In human, studies have shown lower concentration of SCFA and butyrate-producing bacteria (i.e. *Faecalibaterium and Clostridium*) in the gut of inflammatory bowel diseases (IBD) patients suggesting a relevant role of butyrate in intestinal health⁵⁶⁻⁵⁸. Butyrate and other SCFAs elicit most of their biological activities by binding to GPRs and by acting as HDAC inhibitors⁴⁹. Our results demonstrate an additional mechanism where butyrate binds to AhR thus activating its signalling pathway and in synergy with its HDACi property leads to the expression of AhR-dependent genes. These results might be generalized to other SCFAs as most of them activated the AhR signalling pathway (Fig. 2 and Supplementary Fig. S2) and in silico modelling of propionate interaction with AhR showed an orientation in the pocket of human AhR PAS-B domain similar to butyrate involving H291 and Q383 (Supplementary Fig. S3). However additional experiments are needed to validate this point. Interestingly, SCFAs such as butyrate exhibit overlapping activities with AhR ligands on intestinal homeostasis and in IBD⁵⁹. AhR ligands regulate epithelial IL-10 receptor α subunit (IL-10RA) expression that damped colitis by promoting epithelial wound healing⁶⁰. Similarly, butyrate promotes epithelial barrier formation through IL-10RA induction on IECs⁶¹. Hence, it is possible that AhR mediates the butyrate-induced IL-10RA up-regulation. Considering the high quantities of SCFAs in the gut, it is likely that butyrate possibly in synergy with other bacterial-derived AhR ligands have a role in the physiological functions of AhR. Nevertheless, the validation of our results in human primary IECs or colonic organoids is a future challenge to further decipher the mechanisms of butyrate impact on AhR signalling pathway, as well as other ligand-dependent transcription factors.

In conclusion, we show that butyrate stimulates AhR-dependent genes through a direct AhR activation and
probably in complement to its HDACi property in human intestinal cell. Our results suggest that butyrate acts as a
ligand of AhR which is, to our knowledge, an original mechanism only been reported for another ligand-binding
transcription factor PPARγ³⁸.

**Materials and Methods**

**Cell Culture of human colonic cell lines.** The human epithelial cell lines HT-29 and Caco-2 were obtained from the American Type Culture Collection (ATCC, Rockville, MD). HT-29 were grown in RPMI 1640 GlutaMAX™ and Caco-2 in DMEM GlutaMAX™ medium supplemented with 10% and 20% of heat-inactivated fetal bovine serum (FBS, Lonza), respectively. Both media were supplemented with 50 U/mL penicillin, 50 U/mL streptomycin and 10%, 100 mM Hepes, 10 mM nonessential amino acids. HT-29 and Caco-2 were grown at 37 °C in a humidified 5% and 10% CO2 atmosphere, respectively. All culture media and supplements were supplied by Gibco (ThermoFisher). Mycoplasma contamination was regularly tested using MycoAlert (Lonza) and PlasmoTest (Invivogen).

**Production of Stable AhR-luciferase Reporter Cell-Lines.** pGL4.43[luc2P/XRE/Hygro] (Promega) was used to establish HT-29-AhR and Caco-2-AhR reporter cell-lines by electroporation using the Nucleofector® device (Lonza) according to the manufacturer's recommendations. Stable AhR reporter cell lines were selected
using Hygromycin (600 µg/mL for HT-29 and 200 µg/mL for the Caco-2 cell line, InvivoGen) and validated using TCDD at 10 nM final concentration.

**Culture of commensal bacteria, preparation of supernatants and SCFA concentration assessment.** 132 human intestinal commensal bacterial strains (106 different species) from the in-house INRA-Micalis collection or from DSMZ were grown. Anaerobic culture conditions were done accordingly to the Hungate method⁶². Screened strains, corresponding growth media, optical densities (OD600), SCFA concentrations are listed in Appendix Table S1 and composition of home-made growth media is listed in Supplementary Table S2. Bacterial cultures were cultured to reach the maximum OD. Bacterial supernatants were harvest after centrifugation at 5,000 × g for 10 min and filtered on a 0.22 µm PES filters and stored at -80 °C. Quality controls were performed using Gram staining method, aerobic growth test and fresh observation on microscope. Non-inoculated bacteria culture medium served as control. Concentrations of SCFAs produced by cultured bacteria were measured by HPLC and gas chromatography as described by⁶³,⁶⁴.

**Luciferase Reporter and Cell Viability Assays.** For the bacterial screening, HT-29-AhR cells were seeded at 3 × 10⁴ cells per well in white 96-well plates (Corning). After 24 h from seeding, cells were stimulated during 24 hours with 10 µL of bacterial supernatant or non-inoculated media in a total culture-volume of 100 µL per well (10% vol/vol). The screening was performed in triplicates and for almost all the samples; experiments were performed at least with two biological replicates. Additionally, when possible, some strains were grown in different bacterial media. For testing the effect of reagents on AhR activity on HT-29 and Caco-2 cells, 24 h after seeding the culture media was replaced with a non-FBS-supplemented RPMI or DMEM. The cells were then stimulated with 10 µL of reagents diluted in non-FBS conditions in a total culture-volume of 100 µL per well (10% vol/vol). Follow-up experiments were performed in triplicates and repeated at least three times. Luciferase activity was quantified as relative luminescence units (RLU) using a microplate reader (infinite® 200 plate reader, TECAN) and the Neolite™ (PerkinElmer) Luciferase Assay System according to the manufacturer's instructions. The AhR activation was normalised on non-inoculated bacterial media or untreated/vehicle-treated cells for bacterial supernatants and tested reagents, respectively. The results were expressed as luciferase fold change. Cell viability was assessed by MTS measurement using the CellTiter 96 Aqueous One solution (Promega) according to the manufacturer's recommendations.

**Reagents.** All agonists, antagonists and drugs tested were dissolved in a proper vehicle (DMSO, glycerol, water, PBS or ethanol) following the manufacturer's recommendations. The final concentration used for vehicles had not detectable effect on metabolic activity of the cells. Sodium salts of tested SCFAs were from Sigma and used in a range of concentration from 0.125 to 8 mM (20 mM for acetate). AhR agonist: 2,3,7,8-Tetrachlorodibenzodioxin (TCDD 10 nM, Sigma). GPRs agonists: GPR41: 4-chloro-α-(1-methylethyl)-N-2-thiazolylbenzeneacetamide (4-CMTB 1 µM, Tocris) and Tiglic acid (1 mM, Sigma); GPR43: N-(2,5-Dichlorophenyl)-4-(furan-2-yl)-2-methyl-5-oxo-1,4,5,6,7,8-hexahydro-quinoline-3-carboxamide (AR420626 1 µM, Cayman) and 1-methylcyclopropane carboxylate (MCPC 1 mM, Sigma); GPR109a: Niacine (1 mM, Sigma) and (4aR, 5aR)-4,4a,5,5a-Tetrahydro-1H-cyclopropa[4,5]cyclopenta[1,2]pyrazole-3-carboxylic acid (MK1903 1 µM, Tocris). Pertussis toxin (Ptx at 0.2 µg/mL, Sigma) was used as $G\alpha_q$-subunit inhibitor. MCT1 inhibitor used was p-Chloromercuribenzoate acid (pCMB 100 µM, Sigma). HDAC inhibitors: Trichostatin A (TSA 0.1 and 1 µM, Sigma), vorinostat (SAHA 5 µM, Sigma) and valproic acid (VPA 5 mM, Sigma). AhR antagonists: CH-22319 (1 µM, Millipore/Calbiochem), GNF-351 (1 µM, Millipore/Calbiochem), (–)Epigallocatechin gallate (20 µM, EGCG, Sigma).

**Real-Time PCR.** Cell lines were seeded in 12-well culture plates at densities of 0.5 × 10⁶ cells per well. The cells were seeded in FBS-supplemented media then, after 24 h, the media was replaced with a non-FBS-supplemented and cells incubated during 24 h before stimulation. After stimulation time of 6 h, total RNA was extracted using RNeasy mini-Kit (Qiagen) according to manufacturer's recommendations. cDNA was synthesized from 2 µg of RNA using the High-Capacity cDNA Archive Kit (Applied Biosystems). qPCRs were carried out using a StepOne (Applied Biosystems) thermal cycler in a reaction volume of 20 µL with Taqman gene expression assay probes: **AHR**: Hs00169233_m1; **AHRR**: Hs01005075_m1; **CYP1A1**: Hs01054796_g1; **β-ACTIN**: Hs99999903_m1. **AHR**, **CYP1A1** and **AHRR** expression relative to control expression was determined by the 2^(-ΔΔCt) method using **β-actin** as control gene. Data are means ± SEM of at least three distinct experiments, performed in triplicate.

**siRNA assays.** 0.5 x 10⁶ HT-29 cells were seeded in 6-well culture plates on day 1 and 25 nM of siRNA were transfected with DharmaFect I on day 2 and 3, following the manufacturer's instructions (Dharmacon). After 48 h of culture, 0.5 x 10⁶ HT-29 cells were seeded in 6-well culture plates (day 4). Incubation with butyrate (2 mM) and TCDD (10 mM) was at day 5 for 6 h prior qRT-PCR assay. siRNA SMARTpool ON-TARGETplus AhR siRNA (L-004990-00-0005) and Non-targeting Pool (D-001810-10-05) were from Dharmacon.

**Cytoplasmic and nuclear protein extraction.** HT-29 cells were seeded at densities of 0.5 × 10⁶ cells per well in 12-well-plates. 24 h after seeding the media was replaced with a non-FBS-supplemented RPMI and cells incubated during 24 h prior to stimulation. When nuclear extracts were not needed, cells were washed twice and lysed in buffer (1% NP40, 150 mM NaCl, 50 mM Tris-HCL pH8, 5 mM EDTA, 1x Complete Protease Inhibitor Cocktail (Roche)). Nucleus were eliminated by centrifugation for 10 minutes 4 °C at 17500 g. For compartments separation, nuclear and cytoplasmic extracts were prepared with NE-PER Nuclear and Cytoplasmic Extraction Reagent Kit (ThermoFisher) according to the manufacture instructions. CER I and NER buffers were supplemented with protease inhibitor cocktail (cOmplete™ ULTRA Tablets, Mini, EASYpack Protease Inhibitor Cocktail, Sigma) prior to use.
**Western Blot analysis.** Protein extracts were run in 10% SDS-PAGE gels and transferred onto PVDF membranes by liquid transfer (Transfer buffer: 192 mM Glycine, 25 mM TrisBase, 20% methanol) at 200 mA during 90 minutes. Membranes were blocked overnight in TBST + 4% BSA (Sigma). Primary antibodies were incubated overnight at 4 °C: anti-AhR (1:500, mouse mAb, clone RTP1, ThermoFisher), anti-Lamin A/C (1:2000, mouse mAb, Cell Signaling), anti-GAPDH (1:4000, mouse mAb, Santa Cruz), anti-Lamp1 (1:2000, mouse mAb, H4A3 from the Developmental Studies Hybridoma Bank (DSHB), H4A3 was deposited to the DSHB by August, J.T./ Hildreth J.E.K. (DSHB hybridoma product H4A3). Secondary mouse horseradish peroxidase-coupled antibody (DAKO) was successively incubated at room temperature for 2 h before detection with the Clarity Western ECL Substrate using the Chemidoc MP System (Bio-Rad). Quantifications were performed using the image Lab software (Bio-Rad). AhR nuclear protein levels were normalised to Lamin A/C protein levels. Lamp1 and GAPDH were used as purification controls for the cytoplasmic proteins.

**Modelling of butyrate binding to AhR.** For modelling the structure of the complex between AhR and butyrate or propionate, HADDOCK software was used<sup>65,66</sup>. HADDOCK is a highly successful modelling approach that makes use of structural knowledge when available to drive the docking procedure. In this case the crystal structure of the heterodimeric HIF-2α:ARNT complex (PDB code 4ZP4)<sup>67</sup> was used as template in order to dock butyrate to human AhR similarly to Bisson *et al.*, using RaptorX<sup>45,46,68,69</sup>. The cavity of AhR PAS-B was defined by CARVER and HOLLOW<sup>70,71</sup>. The figures were generated with PyMOL Molecular Graphic System, version 1.8 Schrödinger, LLC.

**Statistical analysis.** Presented results are representative of at least three independent experiments. The PCA analysis and Spearman correlation were performed using R and RStudio software. Graphics were produced using Prism GraphPad software. The data distribution was tested using D'Agostino-Person omnibus normality test. Normally distributed data was checked using two-sided t test, otherwise, non parametric Wilcoxon signed-rank or Mann-Whitney tests were performed according to the data set. In all tests, ns: $P > 0.05$, *$P \le 0.05$, **$P \le 0.01$, ***$P \le 0.001$, ****$P < 0.0001$.

**References**

1. Postler, T. S. & Ghosh, S. Understanding the Holobiont: How Microbial Metabolites Affect Human Health and Shape the Immune System. *Cell metabolism* **26**, 110–130, https://doi.org/10.1016/j.cmet.2017.05.008 (2017).

2. Gao, J. et al. Impact of the Gut Microbiota on Intestinal Immunity Mediated by Tryptophan Metabolism. *Front Cell Infect Microbiol* **8**, 13, https://doi.org/10.3389/fcimb.2018.00013 (2018).

3. Quintana, F. J. & Sherr, D. H. Aryl hydrocarbon receptor control of adaptive immunity. *Pharmacol Rev* **65**, 1148–1161, https://doi.org/10.1124/pr.113.007823 (2013).

4. Veldhoen, M. & Ferreira, C. Influence of nutrient-derived metabolites on lymphocyte immunity. *Nature medicine* **21**, 709–718, https://doi.org/10.1038/nm.3894 (2015).

5. Barouki, R., Aggerbeck, M., Aggerbeck, L. & Coumoul, X. The aryl hydrocarbon receptor system. *Drug Metabol Drug Interact* **27**, 3–8, https://doi.org/10.1515/dmdi-2011-0035 (2012).

6. Hubbard, T. D. et al. Adaptation of the human aryl hydrocarbon receptor to sense microbiota-derived indoles. *Scientific reports* **5**, 12689, https://doi.org/10.1038/srep12689 (2015).

7. Hubbard, T. D., Murray, I. A. & Perdew, G. H. Indole and Tryptophan Metabolism: Endogenous and Dietary Routes to Ah Receptor Activation. *Drug Metab Dispos* **43**, 1522–1535, https://doi.org/10.1124/dmd.115.064246 (2015).

8. Jin, U. H. et al. Microbiome-derived tryptophan metabolites and their aryl hydrocarbon receptor-dependent agonist and antagonist activities. *Mol Pharmacol* **85**, 777–788, https://doi.org/10.1124/mol.113.091165 (2014).

9. Fujii-Kuriyama, Y., Ema, M., Mimura, J. & Sogawa, K. Ah receptor: a novel ligand-activated transcription factor. *Exp Clin Immunogenet* **11**, 65–74 (1994).

10. Hankinson, O. The aryl hydrocarbon receptor complex. *Annu Rev Pharmacol Toxicol* **35**, 307–340, https://doi.org/10.1146/annurev.pa.35.040195.001515 (1995).

11. Qiu, J. et al. Group 3 innate lymphoid cells inhibit T-cell-mediated intestinal inflammation through aryl hydrocarbon receptor signaling and regulation of microflora. *Immunity* **39**, 386–399, https://doi.org/10.1016/j.immuni.2013.08.002 (2013).

12. Qiu, J. & Zhou, L. Aryl hydrocarbon receptor promotes RORgammat(+) group 3 ILCs and controls intestinal immunity and inflammation. *Seminars in immunopathology* **35**, 657–670, https://doi.org/10.1007/s00281-013-0393-5 (2013).

13. Fukumoto, S. et al. Identification of a probiotic bacteria-derived activator of the aryl hydrocarbon receptor that inhibits colitis. *Immunol Cell Biol* **92**, 460–465, https://doi.org/10.1038/icb.2014.2 (2014).

14. Zelante, T. et al. Tryptophan catabolites from microbiota engage aryl hydrocarbon receptor and balance mucosal reactivity via interleukin-22. *Immunity* **39**, 372–385, https://doi.org/10.1016/j.immuni.2013.08.003 (2013).

15. Lamas, B. et al. CARD9 impacts colitis by altering gut microbiota metabolism of tryptophan into aryl hydrocarbon receptor ligands. *Nature medicine* **22**, 598–605, https://doi.org/10.1038/nm.4102 (2016).

16. Monteleone, I. et al. Aryl hydrocarbon receptor-induced signals up-regulate IL-22 production and inhibit inflammation in the gastrointestinal tract. *Gastroenterology* **141**, 237–248, 248 e231, https://doi.org/10.1053/j.gastro.2011.04.007 (2011).

17. Furumatsu, K. et al. A role of the aryl hydrocarbon receptor in attenuation of colitis. *Dig Dis Sci* **56**, 2532–2544, https://doi.org/10.1007/s10620-011-1643-9 (2011).

18. Goettel, J. A. et al. AHR Activation Is Protective against Colitis Driven by T Cells in Humanized Mice. *Cell reports* **17**, 1318–1329, https://doi.org/10.1016/j.celrep.2016.09.082 (2016).

19. Liu, Z. et al. Aryl hydrocarbon receptor activation maintained the intestinal epithelial barrier function through Notch1 dependent signaling pathway. *Int J Mol Med* **41**, 1560–1572, https://doi.org/10.3892/ijmm.2017.3341 (2018).

20. Yu, M. et al. Aryl Hydrocarbon Receptor Activation Modulates Intestinal Epithelial Barrier Function by Maintaining Tight Junction Integrity. *Int J Biol Sci* **14**, 69–77, https://doi.org/10.7150/ijbs.22259 (2018).

21. Schiering, C., Vonk, A., Das, S., Stockinger, B. & Wincent, E. Cytochrome P4501-inhibiting chemicals amplify aryl hydrocarbon receptor activation and IL-22 production in T helper 17 cells. *Biochem Pharmacol* **151**, 47–58, https://doi.org/10.1016/j.bcp.2018.02.031 (2018).

22. Schiering, C. et al. Feedback control of AHR signalling regulates intestinal immunity. *Nature* **542**, 242–245, https://doi.org/10.1038/nature21080 (2017).

23. Jin, U. H. et al. Short Chain Fatty Acids Enhance Aryl Hydrocarbon (Ah) Responsiveness in Mouse Colonocytes and Caco-2 Human ColonCancer Cells. *Scientific reports* **7**, 10163, https://doi.org/10.1038/s41598-017-10824-x (2017).
24. Sonowal, R. et al. Indoles from commensal bacteria extend healthspan. *Proceedings of the National Academy of Sciences of the United States of America* **114**, E7506-E7515, https://doi.org/10.1073/pnas.1706464114 (2017).

25. Vital, M., Howe, A. C. & Tiedje, J. M. Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data. *mBio* **5**, e00889, https://doi.org/10.1128/mBio.00889-14 (2014).

26. Cummings, J. H. Short chain fatty acids in the human colon. *Gut* **22**, 763-779 (1981).

27. Cummings, J. H., Pomare, E. W., Branch, W. J., Naylor, C. P. & Macfarlane, G. T. Short chain fatty acids in human large intestine, portal, hepatic and venous blood. *Gut* **28**, 1221-1227 (1987).

28. Blachier, F., Mariotti, F., Huneau, J. F. & Tome, D. Effects of amino acid-derived luminal metabolites on the colonic epithelium and physiopathological consequences. *Amino Acids* **33**, 547-562, https://doi.org/10.1007/s00726-006-0477-9 (2007).

29. Schaible, U. E. & Kaufmann, S. H. A nutritive view on the host-pathogen interplay. *Trends in microbiology* **13**, 373-380, https://doi.org/10.1016/j.tim.2005.06.009 (2005).

30. Schilderink, R., Verseijden, C. & de Jonge, W. J. Dietary inhibitors of histone deacetylases in intestinal immunity and homeostasis. *Frontiers in immunology* **4**, 226, https://doi.org/10.3389/fimmu.2013.00226 (2013).

31. Candido, E. P., Reeves, R. & Davie, J. R. Sodium butyrate inhibits histone deacetylation in cultured cells. *Cell* **14**, 105-113 (1978).

32. Sealy, L. & Chalkley, R. The effect of sodium butyrate on histone modification. *Cell* **14**, 115-121 (1978).

33. Hinnebusch, B. F., Meng, S., Wu, J. T., Archer, S. Y. & Hodin, R. A. The effects of short-chain fatty acids on human colon cancer cell phenotype are associated with histone hyperacetylation. *J Nutr* **132**, 1012-1017 (2002).

34. Tolhurst, G. et al. Short-chain fatty acids stimulate glucagon-like peptide-1 secretion via the G-protein-coupled receptor FFAR2. *Diabetes* **61**, 364-371, https://doi.org/10.2337/db11-1019 (2012).

35. Cuff, M., Dyer, J., Jones, M. & Shirazi-Beechey, S. The human colonic monocarboxylate transporter Isoform 1: its potential importance to colonic tissue homeostasis. *Gastroenterology* **128**, 676-686 (2005).

36. Gallinari, P., Di Marco, S., Jones, P., Pallaoro, M. & Steinkuhler, C. HDACs, histone deacetylation and gene transcription: from molecular biology to cancer therapeutics. *Cell Res* **17**, 195-211, https://doi.org/10.1038/sj.cr.7310149 (2007).

37. Kinoshita, M., Suzuki, Y. & Saito, Y. Butyrate reduces colonic paracellular permeability by enhancing PPARgamma activation. *Biochem Biophys Res Commun* **293**, 827-831, https://doi.org/10.1016/S0006-291X(02)00294-2 (2002).

38. Alex, S. et al. Short-chain fatty acids stimulate angiopoietin-like 4 synthesis in human colon adenocarcinoma cells by activating peroxisome proliferator-activated receptor gamma. *Mol Cell Biol* **33**, 1303-1316, https://doi.org/10.1128/MCB.00858-12 (2013).

39. Yin, Z., Henry, E. C. & Gasiewicz, T. A. (-)-Epigallocatechin-3-gallate is a novel Hsp90 inhibitor. *Biochemistry* **48**, 336-345, https://doi.org/10.1021/bi801637q (2009).

40. Perkins, A. et al. A Structural Switch between Agonist and Antagonist Bound Conformations for a Ligand-Optimized Model of the Human Aryl Hydrocarbon Receptor Ligand Binding Domain. *Biology (Basel)* **3**, 645-669, https://doi.org/10.3390/biology3040645 (2014).

41. Smith, K. J. et al. Identification of a high-affinity ligand that exhibits complete aryl hydrocarbon receptor antagonism. *J Pharmacol Exp Ther* **338**, 318-327, https://doi.org/10.1124/jpet.110.178392 (2011).

42. Palermo, C. M., Westlake, C. A. & Gasiewicz, T. A. Epigallocatechin gallate inhibits aryl hydrocarbon receptor gene transcription through an indirect mechanism involving binding to a 90 kDa heat shock protein. *Biochemistry* **44**, 5041-5052, https://doi.org/10.1021/bi047433p (2005).

43. Henry, E. C. & Gasiewicz, T. A. Agonist but not antagonist ligands induce conformational change in the mouse aryl hydrocarbon receptor as detected by partial proteolysis. *Mol Pharmacol* **63**, 392-400 (2003).

44. Bonati, L., Corrada, D., Tagliabue, S. G. & Motta, S. Molecular modeling of the AhR structure and interactions can shed light on ligand-dependent activation and transformation mechanisms. *Curr Opin Toxicol* **2**, 42-49, https://doi.org/10.1016/j.cotox.2017.01.011 (2017).

45. Bisson, W. H. et al. Modeling of the aryl hydrocarbon receptor (AhR) ligand binding domain and its utility in virtual ligand screening to predict new AhR ligands. *J Med Chem* **52**, 5635-5641, https://doi.org/10.1021/jm900199u (2009).

46. Erbel, P. J. et al. Solution structure of the alpha-subunit of human chorionic gonadotropin. *Eur J Biochem* **260**, 490-498 (1999).

47. Stockinger, B., Di Meglio, P., Gialitakis, M. & Duarte, J. H. The aryl hydrocarbon receptor: multitasking in the immune system. *Annu Rev Immunol* **32**, 403-432, https://doi.org/10.1146/annurev-immunol-032713-120245 (2014).

48. Moura-Alves, P. et al. AhR sensing of bacterial pigments regulates antibacterial defence. *Nature* **512**, 387-392, https://doi.org/10.1038/nature13684 (2014).

49. Daly, K. & Shirazi-Beechey, S. P. Microarray analysis of butyrate regulated genes in colonic epithelial cells. *DNA Cell Biol* **25**, 49-62, https://doi.org/10.1089/dna.2006.25.49 (2006).

50. Garrison, P. M., Rogers, J. M., Brackney, W. R. & Denison, M. S. Effects of histone deacetylase inhibitors on the Ah receptor gene promoter. *Arch Biochem Biophys* **374**, 161-171, https://doi.org/10.1006/abbi.1999.1620 (2000).

51. Gill, R. K. et al. Expression and membrane localization of MCT isoforms along the length of the human intestine. *Am J Physiol Cell Physiol* **289**, C846-852, https://doi.org/10.1152/ajpcell.00112.2005 (2005).

52. Hadjiagioupi, C., Schmidt, L., Dudeja, P. K., Layden, T. J. & Ramaswamy, K. Mechanism(s) of butyrate transport in Caco-2 cells: role of monocarboxylate transporter 1. *American journal of physiology. Gastrointestinal and liver physiology* **279**, G775-G780 (2000).

53. Thibault, R. et al. Butyrate utilization by the colonic mucosa in inflammatory bowel diseases: a transport deficiency. *Inflammatory bowel diseases* **16**, 684-695, https://doi.org/10.1002/ibd.21108 (2010).

54. Flaveny, C. A., Murray, I. A., Chiaro, C. R. & Perdew, G. H. Ligand selectivity and gene regulation by the human aryl hydrocarbon receptor in transgenic mice. *Mol Pharmacol* **75**, 1412-1420, https://doi.org/10.1124/mol.109.054825 (2009).

55. Pandini, A., Denison, M. S., Song, Y., Soshilov, A. A. & Bonati, L. Structural and functional characterization of the aryl hydrocarbon receptor ligand binding domain by homology modeling and mutational analysis. *Biochemistry* **46**, 696-708, https://doi.org/10.1021/bi061460t (2007).

56. Sokol, H. et al. Low counts of Faecalibacterium prausnitzii in colitis microbiota. *Inflammatory bowel diseases* **15**, 1183-1189 (2009).

57. Machiels, K. et al. A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. *Gut* **63**, 1275-1283, https://doi.org/10.1136/gutjnl-2013-304833 (2014).

58. Kumari, R., Ahuja, V. & Paul, J. Fluctuations in butyrate-producing bacteria in ulcerative colitis patients of North India. *World J Gastroenterol* **19**, 3404-3414, https://doi.org/10.3748/wj.v19.i22.3404 (2013).

59. Cushing, K., Alvarado, D. M. & Ciorba, M. A. Butyrate and Mucosal Inflammation: New Scientific Evidence Supports Clinical Observation. *Clin Transl Gastroenterol* **6**, e108, https://doi.org/10.1038/ctg.2015.34 (2015).

60. Alexeev, E. E. et al. Microbiota-Derived Indole Metabolites Promote Human and Murine Intestinal Homeostasis through Regulation of Interleukin-10 Receptor. *Am J Pathol*. https://doi.org/10.1016/j.aipath.2018.01.011 (2018).

61. Zheng, L. et al. Microbial-Derived Butyrate Promotes Epithelial Barrier Function through IL-10 Receptor-Dependent Repression of Claudin-2. *Journal of immunology* **199**, 2976-2984, https://doi.org/10.4049/jimmunol.1700105 (2017).

62. Hungate, R. E. The anaerobic mesophilic cellulolytic bacteria. *Bacteriol Rev* **14**, 1-49 (1950).

63. Lakhdari, O. et al. Identification of NF-kappaB modulation capabilities within human intestinal commensal bacteria. *Journal of biomedicine & biotechnology* **2011**, 282356, https://doi.org/10.1155/2011/282356 (2011).

64. Bourriaud, C. et al. Lactate is mainly fermented to butyrate by human intestinal microflora but inter-individual variation is evident. *J Appl Microbiol* **99**, 201-212, https://doi.org/10.1111/j.1365-2672.2005.02605.x (2005).
65. van Zundert, G. C. P. et al. The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. *J Mol Biol* **428**, 720–725, https://doi.org/10.1016/j.jmb.2015.09.014 (2016).

66. Wassenaar, T. et al. WeNMR: Structural Biology on the Grid. *Journal of Grid Computing* **10**, 743–767 (2012).

67. Wu, D., Potluri, N., Lu, J., Kim, Y. & Rastinejad, F. Structural integration in hypoxia-inducible factors. *Nature* **524**, 303–308, https://doi.org/10.1038/nature14883 (2015).

68. Kallberg, M. et al. Template-based protein structure modeling using the RaptorX web server. *Nat Protoc* **7**, 1511–1522, https://doi.org/10.1038/nprot.2012.085 (2012).

69. Ma, J., Peng, J., Wang, S. & Xu, J. A conditional neural fields model for protein threading. *Bioinformatics* **28**, i59–66, https://doi.org/10.1093/bioinformatics/bts213 (2012).

70. Chovancova, E. et al. CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures. *PLoS Comput Biol* **8**, e1002708, https://doi.org/10.1371/journal.pcbi.1002708 (2012).

71. Ho, B. K. & Gruswitz, F. HOLLOW: generating accurate representations of channel and interior surfaces in molecular structures. *BMC Struct Biol* **8**, 49, https://doi.org/10.1186/1472-6807-8-49 (2008).

## Acknowledgements

The authors are grateful to Dr Marion Espeli (INSERM UMR-S996) and to the member of the team Functionality of the Intestinal Ecosystem (Micalis-INRA Jouy-en-Josas) for helpful discussions. The authors are grateful to Agnès David and Catherine Michel for the quantification of SCFA in bacterial cultures. This work was supported by the Institut National de la Recherche Agronomique (INRA), by grants funded by EU-FP7 METACARDIS (HEALTH-F4-2012-305312).

## Author Contributions

Conceived and designed the experiments: L.M., N.L., H.M.B.; performed most of the experiments: L.M.; performed some experiments: C.M.G., F.B.C., J.M.B., N.L.; analysed the data: L.M., C.M.G., N.L.; wrote the paper: L.M., N.L.; edited and revised the manuscript: C.M.G., J.M.B., H.M.B.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-37019-2.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_19.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019