# Gene expression profiling reveals that PXR activation inhibits Hepatic PPAR alpha activity and decreases FGF21 secretion in male C57Bl6/J mice

Sharon Barretto, Frédéric Lasserre, Anne Fougerat, Lorraine Smith, Tiffany Fougeray, Céline Lukowicz, Arnaud Polizzi, Sarra Smati, Marion Régnier, Claire Naylies, et al.

► To cite this version:

Sharon Barretto, Frédéric Lasserre, Anne Fougerat, Lorraine Smith, Tiffany Fougeray, et al. Gene expression profiling reveals that PXR activation inhibits Hepatic PPAR alpha activity and decreases FGF21 secretion in male C57Bl6/J mice. International Journal of Molecular Sciences, 2019, 20 (15), 10.3390/ijms20153767. [DOI: 10.3390/ijms20153767](https://doi.org/10.3390/ijms20153767)

HAL Id: hal-02617929

https://hal.inrae.fr/hal-02617929v1

Submitted on 25 May 2020

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 CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Article

# Gene Expression Profiling Reveals that PXR Activation Inhibits Hepatic PPARα Activity and Decreases FGF21 Secretion in Male C57Bl6/J Mice

Sharon Ann Barretto † ID, Frédéric Lasserre †, Anne Fougerat, Lorraine Smith, Tiffany Fougeray, Céline Lukowicz, Arnaud Polizzi, Sarra Smati, Marion Régnier, Claire Naylies, Colette Bétoulières, Yannick Lippi ID, Hervé Guillou ID, Nicolas Loiseau ID, Laurence Gamet-Payrastre, Laila Mselli-Lakhal and Sandrine Ellero-Simatos * ID

Institut National de la Recherche Agronomique (INRA), UMR1331 Toxalim, F31-027 Toulouse CEDEX 3, France

* Correspondence: sandrine.ellero-simatos@inra.fr

† These authors contributed equally to this work.

Received: 21 June 2019; Accepted: 31 July 2019; Published: 1 August 2019

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

check for
updates

**Abstract:** The pregnane X receptor (PXR) is the main nuclear receptor regulating the expression of xenobiotic-metabolizing enzymes and is highly expressed in the liver and intestine. Recent studies have highlighted its additional role in lipid homeostasis, however, the mechanisms of these regulations are not fully elucidated. We investigated the transcriptomic signature of PXR activation in the liver of adult wild-type vs. *Pxr*⁻/⁻ C57Bl6/J male mice treated with the rodent specific ligand pregnenolone 16α-carbonitrile (PCN). PXR activation increased liver triglyceride accumulation and significantly regulated the expression of 1215 genes, mostly xenobiotic-metabolizing enzymes. Among the down-regulated genes, we identified a strong peroxisome proliferator-activated receptor α (PPARα) signature. Comparison of this signature with a list of fasting-induced PPARα target genes confirmed that PXR activation decreased the expression of more than 25 PPARα target genes, among which was the hepatokine fibroblast growth factor 21 (*Fgf21*). PXR activation abolished plasmatic levels of FGF21. We provide a comprehensive signature of PXR activation in the liver and identify new PXR target genes that might be involved in the steatogenic effect of PXR. Moreover, we show that PXR activation down-regulates hepatic PPARα activity and FGF21 circulation, which could participate in the pleiotropic role of PXR in energy homeostasis.

**Keywords:** nuclear receptors; hepatokines; transcriptomics

## 1. Introduction

Pregnane X receptor (PXR, systematic name NR1I2) is a member of the nuclear receptor superfamily and is highly expressed in the liver and intestine of mammals [1]. PXR was characterized as a xenosensor that regulates the expression of xenobiotic-metabolizing enzymes and transporters, thereby facilitating the elimination of xenobiotics and endogenous toxic chemicals such as bile acids [2]. Upon ligand-binding, PXR translocates to the nucleus, heterodimerizes with retinoid X receptor (RXR, NR2B1) and binds to PXR direct repeat 4 (DR-4) response elements (PXRE) that are usually located upstream of target genes. Because of an unusually large and flexible binding pocket, PXR can be activated by a variety of structurally diverse chemicals, including pharmaceutical drugs, dietary supplements, herbal medicines, environmental pollutants, and endogenous molecules [3]. In line with the role of PXR as a master regulator of xenobiotic metabolism, its first described target gene was cytochrome P450 (CYP) 3A4 in humans [4], which represents 10% of all clinically relevant drug-metabolizing CYPs in the human liver and up to 75%-85% in the intestine [5] and is responsible for the metabolization of 60% of marketed drugs [6].
Besides its original function as part of the detoxification machinery, recent studies have also unveiled functions for PXR in intermediary metabolism. There is an increasing amount of clinical evidence showing that PXR agonists cause hyperglycemia in humans [7] and pre-clinical work suggesting that PXR regulates hepatic glucose metabolism, however, there is still no solid understanding of the consequences, or of the mechanisms involved. Activated PXR has been shown to repress expression of the gluconeogenic genes glucose-6-phosphatase (G6Pase) and phosphoenolpyruvate carboxykinase (PCK1) [8], and of genes involved in glucose uptake such as GLUT2 and of glucokinase (GCK) [9]. Although there is limited data on the relationship between PXR and fatty liver in humans in vivo, many studies have demonstrated that PXR activation also causes hepatic lipid accumulation in human cell models, and in vitro and in vivo mouse models [7,10]. This pro-steatotic effect is thought to result from both the activation of lipogenesis and inhibition of β-oxidation [7]. However, the mechanisms by which PXR activation induces perturbations of lipid metabolism are not fully elucidated. Recently, it was shown that the activation of intestinal PXR signaling induced dyslipidemia and intestinal cholesterol accumulation [11], while activation of hepatic PXR signaling was sufficient to promote hypercholesterolemia and hepatic lipid accumulation [12].

Here, we aimed to gain insights into the potential metabolic dysregulations induced upon PXR activation and performed a transcriptomic comparison of the hepatic gene profiles of wild type (WT) vs. *Pxr*⁻/⁻ male mice treated with the rodent specific PXR ligand pregnenolone 16α-carbonitrile (PCN). As expected, we observed that PCN treatment-induced hepatic steatosis. We unraveled several previously unknown PXR target genes involved in liver lipid accumulation and discovered a very robust peroxisome proliferator-activated receptor α (PPARα) signature amongst the PXR down-regulated target genes. The PXR-induced decrease in PPARα activity included the regulation of the hepatokine FGF21, a liver-derived hormone with major endocrine roles [13]. This cross-talk between PXR and PPARα in the regulation of FGF21 may contribute to endocrine disruption by xenobiotics acting as ligands for PXR.

## 2. Results

### 2.1. *Effect of PXR Activation on Physiological Parameters and Liver Lipids*

We investigated the effect of PXR activation by its pharmacological ligand PCN in WT and *Pxr*⁻/⁻ male mice. PCN treatment did not affect body weight but increased relative liver mass in a PXR-dependent way (Figure 1a). In the liver, PXR activation significantly increased cholesterol esters and triglyceride levels but did not significantly impact free cholesterol (Figure 1a). In the plasma, PXR activation increased alanine transaminase (ALT) and decreased total cholesterol levels but did not impact free fatty acids, triglycerides (Figure 1a), HDL, LDL, or glucose levels (Figure S1).

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Effect of pregnenolone 16α-carbonitrile (PCN) treatment on liver parameters (a) and plasma biochemistry (b). Data are shown as mean ± SEM of *n* = 5–6 per group. \* *p* ≤ 0.05, \*\* *p* ≤ 0.01, \*\*\* *p* ≤ 0.005 for PCN effect using 2-way ANOVA and Tukey's post-tests. \$ *p* ≤ 0.05 for genotype effect. ALT: Alanine amino-transferase; FFA: Free fatty acids.</figcaption></figure>
## 2.2. Effects of PXR Activation on the Hepatic Transcriptome

Using microarrays, we obtained global transcriptional profiles. Principal component analysis (PCA) first illustrated that PCN treatment significantly impacted the hepatic transcriptome (Figure 2a). The discrimination of WT PCN vs. WT Cont seems stronger than that of the *Pxr*⁻/⁻ PCN vs. *Pxr*⁻/⁻ Cont, confirming, as expected, a significant PXR-dependent transcriptional effect of PCN. We next used linear models and considered genes to be significantly regulated with a fold-change >1.5 and a false discovery rate (FDR) <0.05. Heatmap clustering confirmed the PCA results (Figure S2). It indeed revealed five gene clusters with the largest cluster (1602 probes) comprised of genes up-regulated by PCN in WT mice only (cluster 5). Another cluster (cluster 2) showed 407 probes down-regulated upon PCN treatment in WT mice only. Cluster 4 contained 498 probes that showed genes differentially regulated in WT vs. *Pxr*⁻/⁻ mice, independently of PCN. Finally, cluster 3 (605 probes) illustrated a PCN effect in both WT and *Pxr*⁻/⁻ mice.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Impact of pregnane X receptor (PXR) activation on the hepatic transcriptome. (a) Principal component analysis (PCA) score plots of the whole transcriptomic dataset. (b) qPCR confirmation on selected genes. \* *p* ≤ 0.05, \*\* *p* ≤ 0.01, \*\*\* *p* ≤ 0.005 for PCN effect using 2-way ANOVA and Tukey's post-tests. (c,f) Venn diagram representing the number of genes affected by PCN treatment. (d,g) Gene enrichment analyses of the PXR-target genes. (e,h) The 30 genes with the highest fold-changes upon PCN treatment.</figcaption></figure>

We next sought to decipher the biological functions affected by PXR activation. PCN treatment significantly up-regulated the expression of 1258 genes in WT animals, and of 333 genes in *Pxr*⁻/⁻ mice (Figure 2c). Using the 1029 "prototypical" PXR target genes (those that were up-regulated only in WT animals), we conducted a pathway enrichment analysis, which revealed seven functional clusters significantly enriched (Figure 2d and Table S1) with genes involved in cell cycle, cell division and mitosis, glutathione metabolism, cytochromes P450, lipid metabolism, chemotaxis, and positive regulation of inflammatory response. Figure 2e confirms these results by illustrating the fold-changes of the top 30 most up-regulated genes. These results first confirmed the well-described influence
of PXR activation on hepatic xenobiotic-metabolizing enzymes, mainly those from the Cyp3 family. Table S2 provides a full description of the impact of PCN treatment on all xenobiotic-metabolizing enzymes. Induction of two of the most well-described PXR targets, *Cyp2c55* and *Cyp3a11* were further confirmed using RT-qPCR (Figure 2b). Interestingly, the “lipid metabolism” pathway was also highly significantly enriched upon PXR activation and, among the 30 genes with the highest fold-change, the patatin-like phospholipase domain containing 3 (*Pnpla3*), the thyroid hormone-responsive spot 14 (*Thrsp* or *Spot14*), and the growth/differentiation factor 15 (*Gdf15*) belonged to this pathway. Induction of *Gdf15* was also confirmed by RT-qPCR (Figure 2b). Finally, the regulation of genes involved in de novo lipogenesis was also confirmed by qPCR and showed a significant increase of the SREBP-1 lipogenic pathway in *Pxr*⁻/⁻ mice compared to WT mice (Figure S3).

We next investigated the effect of PCN on gene down-regulation. PCN treatment significantly decreased the expression of 186 genes in a PXR-dependent manner (Figure 2f). GO analyses revealed that these genes were involved in lipid metabolic process, biological rhythms, transforming growth factor-β (TGFβ) signaling pathway, glucose metabolism, and cytochromes P450 (Figure 2g). The 30 genes with the highest fold-changes are illustrated in Figure 2h. Interestingly, among these 30 genes, five (namely *Fgf21*, *Cyp4a10*, *Cyp4a31*, *Acot1*, and *Plin4*) are well-described target genes of PPARα, a key hepatic transcriptional regulator involved in lipid homeostasis.

## 2.3. *Comparison of PXR and PPARα-Dependent Transcriptome*

This prompted us to investigate the intersection between PXR and PPARα activation to test the hypothesis that PXR activation influenced PPARα activity. We took advantage of our previously published microarray dataset [14], in which C57Bl6/J male mice carrying an hepatocyte-specific deletion of *Ppara* (*Ppara*<sup>hep-/-</sup>) were fasted for 24 h to induce PPARα activity and compared to their wild-type littermates (*Ppara*<sup>hep+/+</sup>). We have indeed previously shown that, during fasting, PPARα senses increased levels of free fatty acids released from adipocytes, and in response, controls the expression of hundreds of genes involved in fatty acid uptake, transport, and catabolism in hepatocytes [14,15]. Figure 3a indeed illustrates that a large number of the fasting-induced hepatic genes are PPARα sensitive, with 538 genes significantly up-regulated in a PPARα-dependent manner. Also, 461 genes were significantly down-regulated in a PPARα-dependent manner upon fasting (Figure 3d). We compared these genes with those regulated upon PXR activation. We found 27 genes that were both up-regulated upon PPARα activation and down-regulated upon PXR activation (Figure 3b). These genes are illustrated in Figure 3c and include, among others, *Ppara* itself, *Cyp4a14*, *Cyp4a10*, *Cyp4a31*, and *Fgf21*. There were also 46 genes that were regulated in the opposite direction, i.e., that were down-regulated upon PPARα activation and up-regulated upon PXR activation (Figure 3e,f).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Comparison between PXR and peroxisome proliferator-activated receptor α (PPARα) target genes. (a,d) Venn diagrams representing the number of genes up-(a) or down-(d) regulated upon fasting in Pparα^hep+/+ vs. Pparα^hep-/- mice. (b,e) Venn diagrams representing the number of genes regulated upon PPARα (red) or PXR (blue) activation. (c,f) Fold-changes for the genes that are shared in the previous Venn diagrams.</figcaption></figure>

## 2.4. *Regulation of FGF21*

Using RT-qPCR analyses, we confirmed that PXR activation down-regulated *Pparα* and its target genes expression (Figure 4a), among which was *Fgf21*. FGF21 is a recently described hepatokine with systemic metabolic effects [16]. We measured plasmatic FGF21 and confirmed that circulating FGF21 was decreased upon PCN treatment, since its levels were not detectable anymore in WT-treated mice (Figure 4b). Surprisingly, PXR deletion also influenced FGF21 level since *Pxr*⁻/⁻ mice also showed no detectable levels of circulating FGF21. These differences were not due to different fasting states since glycemia was not significantly different between the four groups (Figure S1).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Impact of PXR activation on hepatic PPARα activity. Gene expression in the liver (a) derived from the microarray and from complementary qPCR experiments. (b) Plasma levels of FGF21. Data are mean ± SEM of *n* = 5–6 per group. \* *p* ≤ 0.05, \*\* *p* ≤ 0.01, \*\*\* *p* ≤ 0.005 for PCN effect, \# *p* ≤ 0.05, \## *p* ≤ 0.01 for genotype effect using 2-way ANOVA and Tukey's post-tests. UDL: Under the detection limit.</figcaption></figure>

**3. Discussion**

The liver is one of the major organs involved in energy production. Hepatic lipid metabolism plays a crucial role during fasting and/or prolonged exercise. Upon lowering of blood glucose, the liver increases glucose production by augmenting gluconeogenesis and glycogenolysis to maintain blood glucose levels, increases fatty acid oxidation and ketogenesis to provide extra-hepatic tissues with ketone bodies, and decreases lipogenesis to attenuate triglyceride storage. These processes are under tight transcriptional control and, in response to hormones such as glucagon and glucocorticoids, many transcription factors cooperate to regulate various genes involved in metabolic pathways aimed at restoring homeostasis [17]. Among those, hepatic PPARα has been well described as crucial for this adaptation. However, recent data have highlighted that other nuclear receptors, such as the aryl H receptor (AhR), the constitutive androstane receptor (CAR), and PXR, which were historically described as xenobiotic sensors, can also interact with the hormone-responsive transcription factors to regulate the liver metabolic processes [18].

Here, we investigated the transcriptomic effects of a pharmacological activation of PXR. The expression of PXR was not described as highly circadian, however, its activity, as measured by the expression of its prototypical target gene *Cyp3a11*, has been shown to be influenced by the time of the day, and is highest as zeitgeber time (ZT)6 [19]. Therefore, we decided to investigate the effects of PXR activation at ZT6, a time at which mice were in a physiological semi-fasted state.

Several studies have already investigated the hepatic signature of PXR activation in vivo [20–22] or in vitro [23]. However, most of these studies focused on the effect of PXR activation on xenobiotic-metabolizing enzymes. Here, we confirm that the regulation of xenobiotic metabolism is one of PXR's most potent functions in the hepatocytes (Figure 2; Table S3). However, our gene enrichment analyses
also revealed that lipid metabolism was among the top-dysregulated pathways upon PXR activation, considering both the up-regulated, as well as the down-regulated genes.

First, PXR activation induced a very significant decrease in plasma cholesterol levels and a significant increase in liver triglycerides and cholesterol esters (Figure 1). The pro-steatotic effects of acute PXR activation have been shown in many studies. However, its role in the regulation of cholesterol homeostasis is more controversial. The anti-HIV drug Efavirenz has been recently shown to induce steatosis and hypercholesterolemia, an effect that was absent in a model of hepatic deletion of PXR [12]. These perturbations were mediated through increased fatty acid transport and cholesterol synthesis, via the PXR-dependent regulation of *Cd36* and *Sqle*. In our data, we confirmed that PXR activation significantly affected *Cd36* and other transporters involved in cholesterol transport, but did not observe any regulation of genes involved in cholesterol biosynthesis, such as *Cyp7a1*, *Sqle*, and *Hmgcr* (Figure S4). This resulted in decreased circulating cholesterol.

Among the up-regulated genes in the liver, we observed that PXR activation increased the expression of several genes that correlate with lipogenesis, such as the patatin-like phospholipase domain containing 3 (*Pnpla3*) and the thyroid hormone-responsive spot 14 (*Thrsp* or *Spot14*). *Spot14* was first identified as a thyroid-responsive gene and is known to transduce hormone- and nutrient-related signals to genes involved in lipogenesis [24]. Regulation of *SPOT14* by PXR was previously described in human hepatocytes [25] and led to increased fatty acid synthase (FASN) expression and triglyceride accumulation. The PNPLA3 protein has lipase activity towards triglycerides in hepatocytes and a loss-of-function polymorphism of this gene has been shown to be strongly associated with nonalcoholic fatty liver disease [26]. However, to our knowledge, the regulation of *Pnpla3* expression by PXR has not been previously described. Among the lipid-metabolic-related genes, we also observed that the expression of the growth/differentiation factor 15 (*Gdf15*), also known as *MIC-1*, was increased by a factor of four upon PCN treatment, in a PXR-dependent way. GDF15 is a distant member of the transforming growth factor-β (TGF-β) superfamily that is considered a crucial hormone in regulating lipid and carbohydrate metabolism. In animal models, overexpression of GDF15 leads to a lean phenotype and improvements of metabolic parameters by increasing the expression of key thermogenic and lipolytic genes in brown and white adipose tissue [27]. Hepatic and circulating GDF15 levels were also increased in animals with blunted β-oxidation (*Cpt2*<sup>hep/-</sup> mice) to maintain systemic energy homeostasis upon fasting [28]. Whether the observed increase in *Gdf15* mRNA upon PCN treatment results from direct regulation of *Gdf15* by PXR or represents a secondary adaptation to decreased β-oxidation remains to be determined. In both cases, regulation of GDF15 levels upon PXR activation might be of physiological relevance since GDF15 has been implicated in a wide variety of biological functions including control of food intake and body weight [29].

Among the genes that were down-regulated upon PXR activation, we observed a very consistent PPARα-like signature, with the decreased expression of many *Cyp4* genes, which are highly sensitive PPARα target genes [14,15]. These results coincide with previous findings in which PCN decreased the hepatic expression of *Ppara*, *Cyp4a10*, and *Cyp4a14* [21]. Neonatal exposure to a single dose of PCN also persistently down-regulated *Cyp4a* expression and decreased PPARα binding to the *Cyp4a* gene loci in adult mice [20]. By comparing the list of genes down-regulated upon PXR activation to a list of genes up-regulated upon PPARα activation, we here extend these previous findings and demonstrate that the inhibition of PPARα activity by PXR affects more than the expression of *Cyp4* genes. For example, the PXR-PPARα interaction probably inhibited the expression of the acetyl-Coenzyme A acyltransferase 1B (*Acaa1b*), of the acyl-coA thioesterase 3 (*Acot3*), of *Krt23* and *Rab30*, of the rate-limiting enzyme in ketogenesis 3-hydroxy-3-methylglutaryl-CoenzymeA synthase 2 (*Hmgcs2*) and of the hepatokine *Fgf21*, all of which are well-described PPARα targets [14]. Using a similar approach in human primary hepatocytes treated with the hPXR ligand rifampicine and the hPPARα ligand WY14643, Kandel et al. had previously shown that more than 14 genes were responsive to both WY14643 (up-regulated) and to rifampicine (down-regulated), among which ACAA2, CYP4A11, and HMGCS2 [23], therefore suggesting the human relevance of our results.
FGF21 is predominantly produced in the liver [30] and exerts pleiotropic effects on the body to maintain overall metabolic homeostasis. FGF21 metabolic benefits range from reducing body weight to alleviating hyperglycemia, insulin resistance, and improvement of lipid profiles [16]. In animal models of obesity, as well as in obese patients, FGF21 has been shown to induce body weight loss and to increase insulin sensitivity and lipid homeostasis [30]. The effects of FGF21 on fertility, growth, and longevity are also well documented [31,32]. Finally, FGF21 seems to be involved in food preferences. For example, FGF21 production in response to carbohydrate intake significantly decreases sugar preferences [33].

Although PXR is mainly expressed in the liver and in the intestine, and not in adipose tissue [34], deletion of *Pxr* appears to influence insulin sensitivity in white adipose tissue and in the muscle [35], serum leptin, and adiponectin levels [36] and PXR activation regulates gene expression in both white and brown adipose tissues [37]. This suggests systemic effects of *Pxr* deletion and activation for which mechanisms have not been described yet. White and brown adipose tissues are among the most described target tissues of FGF21 [16]. Whether FGF21 could be an effector of the systemic effects of PXR remains an open question. Here, we demonstrate that both PXR-activation and PXR deletion decrease the hepatic *Fgf21* mRNA levels and completely abolished the circulating FGF21 levels. This apparent contradictory effect was not limited to the regulation of FGF21 but was also observed in other PPARα target genes (Figure 4). Therefore, it seems that both PXR activation and silencing result in the inhibition of PPARα activity, probably through distinct mechanisms that would need additional investigations. However, it is worth noticing that the same apparent contradictory effect was observed for the regulation of de novo lipogenesis. In human HepG2 cells, PXR activation by rifampicin promoted steatosis via induction of SREBP-1 pathway (mainly SREBP-1α), whereas PXR silencing enhanced AKR1B10 expression, which subsequently stabilized the acetyl-CoA carboxylase, thereby promoting de novo lipogenesis [10]. However, these mechanisms are probably species-specific as, in our data, we did not observe this increase in AKR1B10 expression, whereas the SREBP-1 pathway was increased by PXR ablation and not by PCN treatment (Figure S3). Overall, this demonstrated that complex species-specific mechanisms occur in the regulation of lipogenic pathways by PXR activation and ablation, and our results suggest that this might also be true for the regulation of β-oxidation and PPARα activity.

Perspectives and limitations of our study include the use of male mice only, while PXR activation has been shown to impact both xenobiotic-metabolizing enzymes and glucose and lipid metabolism in a sexually-dimorphic way [38,39]. Therefore, it would be interesting to decipher whether the signature of PXR activation described in our study is also valid in female mice. Second, our study focused on short-term changes. An important remaining question is to determine the effect of multiple weak PXR agonists such as those present in our environment on the observed regulations, especially on FGF21 secretion. Indeed, PXR's main target gene *CYP3A4* is known to be involved in the metabolism of more than 60% of the currently marketed drugs [6] and several hundreds of environmental, occupational, and natural products are demonstrated PXR agonists in both mice and humans [3]. Therefore, regulation of hepatic lipid accumulation by acute or chronic PXR activation might be an important mechanism of xenobiotic-induced steatosis. Finally, the fact that we did not generate the PXR and PPARα dependent transcriptomes in a parallel fashion might have underestimated the number of genes affected by the cross-talk between the two receptors. Therefore, it would be interesting to investigate the effect of PXR activation upon prolonged fasting such as the one used to trigger PPARα. It could also be interesting to decipher whether the pro-steatotic effect of PCN depends on PPARα by treating PPARα knock-out mice with PCN.

Altogether, our results present an additional resource of transcriptome analyses that confirm and extend previous findings on the genes involved in the pro-steatotic effects of PXR. As previously observed in various models [7], we confirm that the observed pro-steatotic effect of PXR activation probably results from both induction of lipogenesis and repression of β-oxidation, and further highlight that this repression is certainly mediated, at least in part, through inhibition of PPARα. We also provide
new hypotheses regarding the yet poorly explored pleiotropic effects of PXR that could result from
the regulation of recently discovered hepatokines, such as GDF15 and/or FGF21. More studies are
needed to confirm the physiological relevance of these regulations. Our findings might have clinical
and public health relevance given the wide range of drugs and environmental xenobiotics that have
been described as PXR ligands and potential endocrine disruptors.

**4. Materials and Methods**

4.1. *Animals*

In vivo studies were performed in a conventional laboratory animal room following the
European Union guidelines for laboratory animal use and care. The current project was approved
by an independent ethics committee (CEEA-86 Toxcométhique) under the authorization number
2018062810452910. The animals were treated humanely with due consideration to the alleviation of
distress and discomfort. All mice were housed at 21-23 °C on a 12 h light (ZT0-ZT12) 12 h dark
(ZT12-ZT24) cycle and allowed free access to the diet (Teklad Global 18% Protein Rodent Diet) and tap
water. ZT stands for Zeitgeber time; ZT0 is defined as the time when the lights are turned on. Twelve
six-week-old wild-type (WT) C57BL/6J male mice were purchased from Charles River and 12 Pxr⁻/⁻
animals (backcrossed on the C57Bl/6J background) were engineered in Pr. Meyer's laboratory [40] and
were bred for 10 y in our animal facility. Mice were acclimatized for two weeks, then randomly allocated
to the different experimental groups: Wild-type control (WT CONT, n = 6), wild-type PCN-treated
(WT PCN, n = 6), Pxr⁻/⁻ control (Pxr⁻/⁻ CONT, n = 6), Pxr⁻/⁻ PCN-treated (Pxr⁻/⁻ PCN, n = 6).
PCN-treated mice received a daily intraperitoneal injection of PCN (100 mg/kg) in corn oil for 4 days while control
mice received corn oil only. Mice were killed at ZT6, 6 h after the last PCN injection.

4.2. *Blood and Tissue Samples*

Bodyweight was monitored at the beginning and at the end of the experimental period. Prior to
sacrifice, the submandibular vein was lanced, and blood was collected into lithium heparin-coated
tubes (BD Microtainer, Franklin Lake, NJ, USA). Plasma was prepared by centrifugation (1500 g, 10 min,
4 °C) and stored at –80 °C. At sacrifice, the liver was removed and snap-frozen in liquid nitrogen and
stored at –80 °C until used for RNA extraction.

4.3. *Gene Expression*

Total RNA was extracted with TRIzol reagent (Invitrogen, Carlsbad, CA, USA). Gene expression
profiles were obtained at the GeT-TRiX facility (GénoToul, Génopole Toulouse Midi-Pyrénées, France)
using Sureprint G3 Mouse GE v2 microarrays (8 × 60 K; design 074,809; Agilent Technologies, Santa
Clara, CA, USA) following the manufacturer's instructions Microarray data and experimental details
are available in NCBI's Gene Expression Omnibus [41] and are accessible through GEO Series accession
numbers GSE123804. For real-time quantitative polymerase chain reaction (qPCR), 2 µg RNA samples
were reverse-transcribed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems,
Foster City, CA, USA). Table S3 presents the SYBR Green assay primers. Amplifications were performed
using an ABI Prism 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). qPCR data
were normalized to TATA-box-binding protein mRNA levels and analyzed with LinRegPCR.v2015.3.

4.4. *Plasma Analysis*

Alanine transaminase (ALT), total cholesterol, triglycerides and free fatty acids (FFA) were
determined using a Pentra 400 biochemical analyzer (Anexplo facility, Toulouse, France). Plasma
FGF21 was assayed using the rat/mouse FGF21 ELISA kit (EMD Millipore, Billerica, MA, USA)
following the manufacturer's instructions.
## 4.5. Liver Neutral Lipid Analysis

Tissue samples were homogenized in methanol/5 mM EGTA (2:1, *v/v*); then, lipids (corresponding to an equivalent of 2 mg tissue) were extracted according to the Bligh and Dyer method [42] with chloroform/methanol/water (2.5:2.5:2.1, *v/v/v*), in the presence of the following internal standards: glyceryl trinonadecanoate, stigmasterol, and cholesteryl heptadecanoate (Sigma, Saint-Louis, MO, USA). Triglycerides, free cholesterol, and cholesterol esters were analyzed with gas-liquid chromatography on a Focus Thermo Electron system equipped with a Zebron-1 Phenomenex fused-silica capillary column (5 m, 0.25 mm i.d., 0.25 mm film thickness). The oven temperature was programmed to increase from 200 to 350 °C at 5 °C/min, and the carrier gas was hydrogen (0.5 bar). Injector and detector temperatures were 315 °C and 345 °C respectively.

## 4.6. Statistical Analysis

Microarray data were processed using R (http://www.r-project.org, accessed at 22 September 2017) and Bioconductor packages (www.bioconductor.org, accessed at 22 September 2017, v 3.0). Raw data (median signal intensity) were filtered, log2 transformed, corrected for batch effects (microarray washing bath), and normalized using CrossNorm method [43]. Normalized data were first analyzed using Matlab (v2014.8). The principal component analysis was performed using an in-house function. The linear model was fitted using the limma lmFit function [44]. Pair-wise comparisons between biological conditions were applied using specific contrasts. A correction for multiple testing was applied using the Benjamini-Hochberg procedure for false discovery rate (FDR). Probes with FDR ≤ 0.05 and |fold-change| > 1.5 were considered to be differentially expressed between conditions. Gene-annotation enrichment analysis and functional annotation clustering were evaluated using DAVID [45]. For non-microarray data, differential effects were analyzed by analysis of variance followed by Tukey's post-hoc tests. A *p*-value < 0.05 was considered significant.

**Supplementary Materials:** Supplementary materials can be found at http://www.mdpi.com/1422-0067/20/15/3767/s1.

**Author Contributions:** Conceptualization, H.G. and S.E.-S.; methodology, S.A.B., F.L., L.S., C.N., C.B., S.S., C.L., and T.F.; software, Y.L.; formal analysis, S.A.B., Y.L., and S.E.-S.; investigation, S.A.B., F.L., M.R., and A.P.; writing—original draft preparation, S.A.B. and S.E.-S.; writing—review and editing, A.F., C.L., H.G., and S.E.-S.; supervision, H.G., L.G.-P., L.M.-L., N.L., and S.E.-S.; project administration, H.G., L.G.-P., L.M.-L., N.L., and S.E.-S.; funding acquisition, H.G., L.G.-P., L.M.-L., N.L., and S.E.-S.

**Funding:** S.A.B. is supported by a Ph.D. grant from Région Occitanie and INRA AlimH department. This work was supported by grants from Agence Nationale de la Recherche (ANR), Fond Européen de Développement Régional (FEDER), and Région Occitanie. SES is supported by a Joint Programming Initiative (JPI) grant Fatmal.

**Acknowledgments:** We thank all members of the EZOP staff for their careful help with this project. We thank the staff from the Genotoul: Anexplo, Get-TriX, and Metatoul-Lipidomic facilities. *Pxr*⁻/⁻ mice are a generous gift from Pr Steven Kliewer (University of Texas Southwestern Medical School, Dallas, TX, USA). Colony founders were kindly provided by Pr Urs A Meyer (Biozentrum, University of Basel, Basel, Switzerland).

**Conflicts of Interest:** The authors declare no conflict of interest.
## Abbreviations

<table><tr><td>ACAA</td><td>acetyl-Coenzyme A acyltransferase</td></tr><tr><td>Acot</td><td>acyl-coA thioesterase</td></tr><tr><td>AhR</td><td>aryl H receptor</td></tr><tr><td>CAR</td><td>constitutive androstane receptor</td></tr><tr><td>Cpt</td><td>carnitine palmitoyltransferase</td></tr><tr><td>CYP</td><td>cytochrome P450</td></tr><tr><td>DR4</td><td>direct repeat 4</td></tr><tr><td>FASN</td><td>fatty acid synthase</td></tr><tr><td>FDR</td><td>false discovery rate</td></tr><tr><td>FGF21</td><td>fibroblast growth factor 21</td></tr><tr><td>G6Pase</td><td>glucose-6-phosphatase</td></tr><tr><td>GDF15</td><td>growth/differentiation factor 15</td></tr><tr><td>GCK</td><td>glucokinase</td></tr><tr><td>GO</td><td>gene ontology</td></tr><tr><td>GLUT2</td><td>glucose transporter 2</td></tr><tr><td>HMGCS</td><td>3-hydroxy-3-methylglutaryl-CoenzymeA synthase</td></tr><tr><td>hPPARα</td><td>Human PPARα</td></tr><tr><td>hPXR</td><td>human PXR</td></tr><tr><td>Krt23</td><td>keratin 23</td></tr><tr><td>PXR</td><td>pregnane X receptor</td></tr><tr><td>PXRE</td><td>pregnane X receptor response elements</td></tr><tr><td>PCA</td><td>principal component analysis</td></tr><tr><td>PCN</td><td>pregnenolone 16α-carbonitrile</td></tr><tr><td>PEPCK</td><td>phosphoenopyruvate carboxykinase</td></tr><tr><td>PNPLA3</td><td>patatin-like phospholipase domain containing 3</td></tr><tr><td>PPARα</td><td>peroxisome proliferator-activated receptor α</td></tr><tr><td>Rab30</td><td>ras-related protein rab-30</td></tr><tr><td>RT-qPCR</td><td>quantitative reverse transcription PCR</td></tr><tr><td>RXR</td><td>retinoid X receptor</td></tr><tr><td>SPOT14</td><td>thyroid hormone-responsive spot 14</td></tr><tr><td>TGFβ</td><td>transforming growth factor-β</td></tr><tr><td>Thrsp</td><td>thyroid hormone-responsive spot 14</td></tr><tr><td>WT</td><td>wild-type</td></tr><tr><td>ZT</td><td>Zeitgeber time</td></tr></table>

## References

1. Bookout, A.L.; Jeong, Y.; Downes, M.; Yu, R.T.; Evans, R.M.; Mangelsdorf, D.J. Anatomical profiling of nuclear receptor expression reveals a hierarchical transcriptional network. *Cell* **2006**, *126*, 789–799. [CrossRef] [PubMed]

2. Kliewer, S.A.; Moore, J.T.; Wade, L.; Staudinger, J.L.; Watson, M.A.; Jones, S.A.; McKee, D.D.; Oliver, B.B.; Willson, T.M.; Zetterström, R.H.; et al. An orphan nuclear receptor activated by pregnanes defines a novel steroid signaling pathway. *Cell* **1998**, *92*, 73–82. [CrossRef]

3. Hernandez, J.P.; Mota, L.C.; Baldwin, W.S. Activation of CAR and PXR by Dietary, Environmental and Occupational Chemicals Alters Drug Metabolism, Intermediary Metabolism, and Cell Proliferation. *Curr. Pharm. Pers. Med.* **2009**, *7*, 81–105. [CrossRef] [PubMed]

4. Guengerich, F.P. Cytochrome P-450 3A4: Regulation and role in drug metabolism. *Annu. Rev. Pharm. Toxicol.* **1999**, *39*, 1–17. [CrossRef] [PubMed]

5. Drozdzik, M.; Busch, D.; Lapczuk, J.; Müller, J.; Ostrowski, M.; Kurzawski, M.; Oswald, S. Protein Abundance of Clinically Relevant Drug-Metabolizing Enzymes in the Human Liver and Intestine: A Comparative Analysis in Paired Tissue Specimens. *Clin. Pharmacol. Ther.* **2018**, *104*, 515–524. [CrossRef] [PubMed]
6. Yu, J.; Petrie, I.D.; Levy, R.H.; Ragueneau-Majlessi, I. Mechanisms and Clinical Significance of Pharmacokinetic-based Drug-drug Interactions with Drugs Approved by the U.S. Food and Drug Administration in 2017. *Drug Metab. Dispos.* **2019**, *47*, 135-144. [CrossRef]

7. Hakkola, J.; Rysä, J.; Hukkanen, J. Regulation of hepatic energy metabolism by the nuclear receptor PXR. *Biochim. Biophys. Acta* **2016**, *1859*, 1072-1082. [CrossRef]

8. Kodama, S.; Koike, C.; Negishi, M.; Yamamoto, Y. Nuclear receptors CAR and PXR cross talk with FOXO1 to regulate genes that encode drug-metabolizing and gluconeogenic enzymes. *Mol. Cell. Biol.* **2004**, *24*, 7931-7940. [CrossRef]

9. Rysä, J.; Buler, M.; Savolainen, M.J.; Ruskoaho, H.; Hakkola, J.; Hukkanen, J. Pregnane X receptor agonists impair postprandial glucose tolerance. *Clin. Pharmacol. Ther.* **2013**, *93*, 556-563. [CrossRef]

10. Bitter, A.; Rümmele, P.; Klein, K.; Kandel, B.A.; Rieger, J.K.; Nüssler, A.K.; Zanger, U.M.; Trauner, M.; Schwab, M.; Burk, O. Pregnane X receptor activation and silencing promote steatosis of human hepatic cells by distinct lipogenic mechanisms. *Arch. Toxicol.* **2015**, *89*, 2089-2103. [CrossRef]

11. Meng, Z.; Gwag, T.; Sui, Y.; Park, S.-H.; Zhou, X.; Zhou, C. The atypical antipsychotic quetiapine induces hyperlipidemia by activating intestinal PXR signaling. *JCI Insight* **2019**, *4*, 1. [CrossRef]

12. Gwag, T.; Meng, Z.; Sui, Y.; Helsley, R.N.; Park, S.-H.; Wang, S.; Greenberg, R.N.; Zhou, C. Non-nucleoside reverse transcriptase inhibitor efavirenz activates PXR to induce hypercholesterolemia and hepatic steatosis. *J. Hepatol.* **2019**, *70*, 930-940. [CrossRef]

13. Kliewer, S.A.; Mangelsdorf, D.J. A Dozen Years of Discovery: Insights into the Physiology and Pharmacology of FGF21. *Cell Metab.* **2019**, *29*, 246-253. [CrossRef] [PubMed]

14. Régnier, M.; Polizzi, A.; Lippi, Y.; Fouché, E.; Michel, G.; Lukowicz, C.; Smati, S.; Marrot, A.; Lasserre, F.; Naylies, C.; et al. Insights into the role of hepatocyte PPARα activity in response to fasting. *Mol. Cell. Endocrinol.* **2017**, *471*, 75-88. [CrossRef]

15. Montagner, A.; Polizzi, A.; Fouché, E.; Ducheix, S.; Lippi, Y.; Lasserre, F.; Barquissau, V.; Régnier, M.; Lukowicz, C.; Benhamed, F.; et al. Liver PPARα is crucial for whole-body fatty acid homeostasis and is protective against NAFLD. *Gut* **2016**, *65*, 1202-1214. [CrossRef] [PubMed]

16. BonDurant, L.D.; Potthoff, M.J. Fibroblast Growth Factor 21: A Versatile Regulator of Metabolic Homeostasis. *Annu. Rev. Nutr.* **2018**, *38*, 173-196. [CrossRef] [PubMed]

17. Goldstein, I.; Hager, G.L. Transcriptional and Chromatin Regulation during Fasting-The Genomic Era. *Trends Endocrinol. Metab.* **2015**, *26*, 699-710. [CrossRef]

18. Konno, Y.; Negishi, M.; Kodama, S. The roles of nuclear receptors CAR and PXR in hepatic energy metabolism. *Drug Metab. Pharmacokinet.* **2008**, *23*, 8-13. [CrossRef] [PubMed]

19. Montagner, A.; Korecka, A.; Polizzi, A.; Lippi, Y.; Blum, Y.; Canlet, C.; Tremblay-Franco, M.; Gautier-Stein, A.; Burcelin, R.; Yen, Y.-C.; et al. Hepatic circadian clock oscillators and nuclear receptors integrate microbiome-derived signals. *Sci. Rep.* **2016**, *6*, 20127. [CrossRef]

20. Li, C.Y.; Cheng, S.L.; Bammler, T.K.; Cui, J.Y. Editor's Highlight: Neonatal Activation of the Xenobiotic-Sensors PXR and CAR Results in Acute and Persistent Down-regulation of PPARα-Signaling in Mouse Liver. *Toxicol. Sci.* **2016**, *153*, 282-302. [CrossRef]

21. Cui, J.Y.; Klaassen, C.D. RNA-Seq reveals common and unique PXR- and CAR-target gene signatures in the mouse liver transcriptome. *Biochim. Biophys. Acta* **2016**, *1859*, 1198-1217. [CrossRef]

22. Nagahori, H.; Nakamura, K.; Sumida, K.; Ito, S.; Ohtsuki, S. Combining Genomics To Identify the Pathways of Post-Transcriptional Nongenotoxic Signaling and Energy Homeostasis in Livers of Rats Treated with the Pregnane X Receptor Agonist, Pregnenolone Carbonitrile. *J. Proteome Res.* **2017**, *16*, 3634-3645. [CrossRef]

23. Kandel, B.A.; Thomas, M.; Winter, S.; Damm, G.; Seehofer, D.; Burk, O.; Schwab, M.; Zanger, U.M. Genomewide comparison of the inducible transcriptomes of nuclear receptors CAR, PXR and PPARα in primary human hepatocytes. *Biochim. Biophys. Acta* **2016**, *1859*, 1218-1227. [CrossRef]

24. LaFave, L.T.; Augustin, L.B.; Mariash, C.N. S14: Insights from knockout mice. *Endocrinology* **2006**, *147*, 4044-4047. [CrossRef]

25. Moreau, A.; Téruel, C.; Beylot, M.; Albalea, V.; Tamasi, V.; Umbdenstock, T.; Parmentier, Y.; Sa-Cunha, A.; Suc, B.; Fabre, J.-M.; et al. A novel pregnane X receptor and S14-mediated lipogenic pathway in human hepatocyte. *Hepatology* **2009**, *49*, 2068-2079. [CrossRef]

26. Dai, G.; Liu, P.; Li, X.; Zhou, X.; He, S. Association between PNPLA3 rs738409 polymorphism and nonalcoholic fatty liver disease (NAFLD) susceptibility and severity: A meta-analysis. *Medicine* **2019**, *98*, e14324. [CrossRef]
27. Chrysovergis, K.; Wang, X.; Kosak, J.; Lee, S.-H.; Kim, J.S.; Foley, J.F.; Travlos, G.; Singh, S.; Baek, S.J.; Eling, T.E. NAG-1/GDF-15 prevents obesity by increasing thermogenesis, lipolysis and oxidative metabolism. *Int. J. Obes. (Lond)* **2014**, *38*, 1555–1564. [CrossRef]

28. Lee, J.; Choi, J.; Scafidi, S.; Wolfgang, M.J. Hepatic Fatty Acid Oxidation Restrains Systemic Catabolism during Starvation. *Cell Rep.* **2016**, *16*, 201–212. [CrossRef]

29. Mullican, S.E.; Lin-Schmidt, X.; Chin, C.-N.; Chavez, J.A.; Furman, J.L.; Armstrong, A.A.; Beck, S.C.; South, V.J.; Dinh, T.Q.; Cash-Mason, T.D.; et al. GFRAL is the receptor for GDF15 and the ligand promotes weight loss in mice and nonhuman primates. *Nat. Med.* **2017**, *23*, 1150–1157. [CrossRef]

30. Markan, K.R.; Naber, M.C.; Ameka, M.K.; Anderegg, M.D.; Mangelsdorf, D.J.; Kliewer, S.A.; Mohammadi, M.; Potthoff, M.J. Circulating FGF21 is liver derived and enhances glucose uptake during refeeding and overfeeding. *Diabetes* **2014**, *63*, 4057–4063. [CrossRef]

31. Zhang, Y.; Xie, Y.; Berglund, E.D.; Coate, K.C.; He, T.T.; Katafuchi, T.; Xiao, G.; Potthoff, M.J.; Wei, W.; Wan, Y.; et al. The starvation hormone, fibroblast growth factor-21, extends lifespan in mice. *Elife* **2012**, *1*, e00065. [CrossRef]

32. Owen, B.M.; Bookout, A.L.; Ding, X.; Lin, V.Y.; Atkin, S.D.; Gautron, L.; Kliewer, S.A.; Mangelsdorf, D.J. FGF21 contributes to neuroendocrine control of female reproduction. *Nat. Med.* **2013**, *19*, 1153–1156. [CrossRef]

33. Iroz, A.; Montagner, A.; Benhamed, F.; Levavasseur, F.; Polizzi, A.; Anthony, E.; Régnier, M.; Fouché, E.; Lukowicz, C.; Cauzac, M.; et al. A Specific ChREBP and PPARα Cross-Talk Is Required for the Glucose-Mediated FGF21 Response. *Cell Rep.* **2017**, *21*, 403–416. [CrossRef]

34. Ellero-Simatos, S.; Chakhtoura, G.; Barreau, C.; Langouët, S.; Benelli, C.; Penicaud, L.; Beaune, P.; de Waziers, I. Xenobiotic-metabolizing cytochromes p450 in human white adipose tissue: Expression and induction. *Drug Metab. Dispos.* **2010**, *38*, 679–686. [CrossRef]

35. He, J.; Gao, J.; Xu, M.; Ren, S.; Stefanovic-Racic, M.; O'Doherty, R.M.; Xie, W. PXR ablation alleviates diet-induced and genetic obesity and insulin resistance in mice. *Diabetes* **2013**, *62*, 1876–1887. [CrossRef]

36. Spruiell, K.; Jones, D.Z.; Cullen, J.M.; Awumey, E.M.; Gonzalez, F.J.; Gyamfi, M.A. Role of human pregnane X receptor in high fat diet-induced obesity in pre-menopausal female mice. *Biochem. Pharm.* **2014**, *89*, 399–412. [CrossRef]

37. Ma, Y.; Liu, D. Activation of pregnane X receptor by pregnenolone 16 α-carbonitrile prevents high-fat diet-induced obesity in AKR/J mice. *PLoS ONE* **2012**, *7*, e38734. [CrossRef]

38. Lu, Y.-F.; Jin, T.; Xu, Y.; Zhang, D.; Wu, Q.; Zhang, Y.-K.J.; Liu, J. Sex differences in the circadian variation of cytochrome p450 genes and corresponding nuclear receptors in mouse liver. *Chronobiol. Int.* **2013**, *30*, 1135–1143. [CrossRef]

39. Spruiell, K.; Gyamfi, A.A.; Yeyeodu, S.T.; Richardson, R.M.; Gonzalez, F.J.; Gyamfi, M.A. Pregnane X Receptor-Humanized Mice Recapitulate Gender Differences in Ethanol Metabolism but Not Hepatotoxicity. *J. Pharm. Exp.* **2015**, *354*, 459–470. [CrossRef]

40. Staudinger, J.L.; Goodwin, B.; Jones, S.A.; Hawkins-Brown, D.; MacKenzie, K.I.; LaTour, A.; Liu, Y.; Klaassen, C.D.; Brown, K.K.; Reinhard, J.; et al. The nuclear receptor PXR is a lithocholic acid sensor that protects against liver toxicity. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 3369–3374. [CrossRef]

41. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. *Nucleic Acids Res.* **2002**, *30*, 207–210. [CrossRef]

42. Bligh, E.G.; Dyer, W.J. A RAPID METHOD OF TOTAL LIPID EXTRACTION AND PURIFICATION. *Can. J. Biochem. Physiol.* **2011**, *37*, 911–917. [CrossRef]

43. Cheng, L.; Lo, L.-Y.; Tang, N.L.S.; Wang, D.; Leung, K.-S. CrossNorm: A novel normalization strategy for microarray data in cancers. *Sci. Rep.* **2016**, *6*, 18898. [CrossRef]

44. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.* **2004**, *3*, [CrossRef]

45. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. *Nucleic Acids Res.* **2009**, *37*, 1–13. [CrossRef]

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