Lipid transport function is the main target of oral oleoylethanolamide to reduce adiposity in high-fat-fed mice.

Abstract We evaluated the biological basis of reduced fat gain by oleoylethanolamide (OEA) in high-fat-fed mice and sought to determine how degradation of OEA affected its efficiency by comparing its effects to those of KDS-5104, a nonhydrolyzable lipid OEA analog. Mice were given OEA or KDS-5104 by the oral route (100 mg/kg body weight). Sixty-eight variables per mouse, describing six biological processes (lipid transport, lipogenesis, energy intake, energy expenditure, endocannabinoid signaling, and glucose metabolism), spanning gene expression of biochemical and physiological parameters were examined to determine the primary target whereby OEA reduces fat gain. Although KDS-5104 but not OEA was resistant to fatty acid amide hydrolase hydrolysis, OEA was degraded by an unidentified hydrolysis system in the liver. Nevertheless, both compounds equally decreased body fat pads after 5 weeks (20%; P < 0.05). The six biological functions constructed from the 68 initial variables predicted up to 58% of adipose fat variations. Lipid transport appeared central to the explanation for body fat deposition (16%; P < 0.0001), in which decreased expression of the FAT/CD36 gene was the component most related to adipose depots. Lipid transport appears to be a determinant player in the OEA fat-lowering response, with adipose tissue FAT/CD36 expression being the most relevant bioindicator of OEA action.

washed three times using warm saline solution to eliminate secondary products. KDS-5104 purity was checked using a 1200 series Agilent Rapid Resolution Liquid Chromatography (RRLC) (Waldbronn, Germany) coupled through an ESI interface to a Micro Quadrupole Time of Light (microQTOF TM , Bruker Daltonik, Bremen, Germany) (99% purity).

Animals, diet, and experimental design
All experiments were conducted according to the French Regulations for Animal Experimentation (Art 19. Oct 1987, Ministry of Agriculture) and in conformity with the Public Health Service Policy after approval by our institutional Animal Care and Use Committee. Male C57Bl6j mice purchased at 8 weeks of age from Janvier Elevage (Le Genest-St.-Isle, France) were individually housed and fed ad libitum on a high-fat diet (50% of energy intake as lipids) for 2 weeks. For each kilogram, the composition of the high-fat diet was: 235 g of casein, 201 g of sucrose, 3.5 g of L -cysteine, 85 g of starch, 116 g of maltodextrin, 60 g of cellulose, 12 g of vitamins mix, 51.5 g of mineral mix, 18.4 g of canola oil, and 236 g of lard [Unité de Préparation des Aliments Expérimentaux (UPAE) INRA Jouy en Josas, France]. Subsequently, 30 mice were blocked by body weight (BW) and randomly allocated to three groups (n = 10 per group) to be fed the same high-fat diet for the subsequent 5 weeks' experimental period. Two groups of mice were given OEA (100 mg/kg BW) or KDS-5104 (100 mg/kg BW) by addition to the diet, while the third group received the basal diet alone. Daily food intake was monitored, and mice were weighed three times per week. A second group of 24 mice was used to determine fat absorption. These mice were fed the same high-fat diet but with either 0.0, 0.1, or 1.0 g of OEA per kilogram of BW by addition to the diet for 5 weeks. During the fourth week of feeding, the mice were kept for 5 days in metabolic cages to collect feces for evaluation of fat absorption (see supplementary material).

Sampling
At the end of the experimental period, mice were euthanized by exsanguination via cardiac puncture after anesthesia with isofl urane (Abbot France, Rungis, France). Plasma was obtained by centrifugation (1,000 g for 10 min at 4°C). Mesenteric, epididymal, inguinal, and peritoneal adipose depots, as well as liver, stomach, small intestine mucosa, pancreas, and gastrocnemius muscle were removed and frozen in liquid nitrogen. Plasma, feces, and organ samples were kept at Ϫ 80°C until analysis.

Biological analysis
All biological parameters and methods related to biochemical, physiological (including indirect calorimetry), gene expression, in vitro and ex vivo assays, and statistical validations are fully detailed in the supplementary material.

Statistical analysis
Univariate analysis. Results are presented as means ± standard error of the means (SEM). Treatment effects on physiological parameters and gene expression data were tested by one-way ANOVA. A mixed-effects model for repeated measures with a fi rst-order autoregressive variance covariance matrix was used, with each daily food intake, lipid oxidation (LOX), carbohydrate oxidation (GOX), and EE as dependent variables and group (OEA or KDS-5104) as independent variable. The effect of the interaction between Group and Time variables was tested with and without total activity adjustment (except for food intake). Statistical signifi cance was set at a P level of <0.05. All analyses were performed with Statview and SAS version 9.1 software (SAS institute, Cary, NC). content in hepatocytes, as well as through a decrease in serum cholesterol and triglyceride levels ( 9 ). In addition, increased lipid ␤ oxidation in muscle and hepatic cells has been described ( 11 ), suggesting that OEA induces fat utilization as an energy substrate. Nonetheless, a direct evaluation of OEA on whole-body energy expenditure (EE) has remained unexplored.
Similar to results shown with intraperitoneal administration, previous rodent studies have shown that orally administered ( 8 ) or diet-supplemented OEA (e.g., see Ref. 14) is biologically effective, although the extent of its breakdown and, thus, its bioeffi ciency is not known. The only clinical study performed did not fi nd an anorectic effect of OEA, but the dietary supplement used was phosphatidylethanolamine given as a parent OEA metabolic precursor ( 12 ).
From the above work, data regarding the various possible OEA modes of action to reduce body fat accretion appear only partial and fragmented. Furthermore, most of the reports make the assumption that the primary OEA response can be explained by the modulated expression of one or a few gene products such as PPAR-␣ , FAT/CD36, UCP2 ( 8 ), whereas the body fat modulation by OEA seems highly multifactorial ( 13 ). Examining together a large number of representative molecular indicators of biological functions assumed to be targeted by the treatment appears more likely to determine the complex mechanisms involved ( 14 ). However, with increasing complexity, it is diffi cult to interpret biological information if the numerous biomodulators are examined independently. The picture can be simplifi ed and more easily interpreted through the use of a biological modules model ( 15 ). This perspective can be achieved by grouping and modeling the individual molecular players that belong to the same biological string into one composite bioindicator that integrates and summarizes the information ( 16 ).
By comparing the effects of long-term oral OEA supplementation with those of its nonhydrolyzable analog KDS-5104 and using an integrated modeling approach embedding molecular, biochemical, and physiological assessments, our goal was also to evaluate the infl uence of enzymatic breakdown of OEA on its ability to reduce body fat gain. Our primary focus was to identify the possible molecular targets of OEA that are relevant to body fat mass, as well as the underlying metabolic functions most likely affected in this process. and R-( Ϫ )-2-aminopropanol (Sigma-Aldrich) under stoichiometric conditions. Acyl chloride was preferred to the triglyceride form for oleic acid in order to induce a total reaction and to avoid heating conditions. Then, the synthesized molecule was (Q2 after cross-validation; P values after cross-validation, ANOVA) are shown in supplementary Table II. Hierarchical clustering of signifi cant variables was performed using Permutmatrix software (http://www.lirmm.fr/~caraux/ PermutMatrix/PermutMatrixHome.htm) with Pearson correlation as distance and the Ward method as clustering condition. In addition, a correlation network displaying the strongest associations among the biological modules computed from the hierarchical PLS was built using Cytoscape open source freeware ( 18 ). The number of edges of the network, which represent the Pearson correlation coeffi cient between the biological module variables, was defi ned from the simplest network to display the entire biological modules. Thus, of the 21 possible relationships theoretically possible, only the 7 strongest are shown (correlation coeffi cient among the variables ranged from 0.703 to 0.85), validated by the false discovery rate for multiple testing.
Oral OEA and KDS-5104 both induced a signifi cant decrease in daily food intake throughout the experimental period ( P = 0.001). Consequently, fi nal cumulative food intake was signifi cantly decreased after OEA ( Ϫ 5.4%, P = 0.0033) and KDS-5104 ( Ϫ 6.3%, P = 0.0004) administration ( Table 2 ). Indeed, correction of adipose fat weight gain for cumulative food intake by a generalized linear modeling abrogated the statistical differences among the treatment groups.
Based on respiratory quotient, OEA and KDS-5104 administration induced a signifi cant increase in LOX relative to GOX , resulting in a 13% increase in LOX ratio over time (LOX/GOX + LOX; 2-way ANOVA; P <0.0001 and P = 0.0001, respectively) ( Table 3 ). Analysis of EE in different treatment groups, calculated from the volume of oxygen consumed (VO 2 ) and the volume of carbon dioxide produced (VCO 2 ), revealed a signifi cant increase in EE over the 24-h period (+6%, P < 0.0001) and during the postprandial period (+7%, P < 0.0001) for the OEA group but not for the KDS-5104 group, compared with that of the control group ( Table 3 ). Total 24-h spontaneous activity was increased in groups treated with OEA (+216%, P = 0.0008) and KDS-5104 (+226%, P = 0.0002) compared with controls ( Table 3 ). The relative lipid oxidative consumption levels and the EE were strongly correlated with the total spontaneous activity in every group ( P < 0.0001 for LOX and EE, data not shown). Consequently, our statistical model was adjusted to account for the total activity level. Results obtained with the adjusted model revealed that OEA and KDS-5104 similarly, signifi cantly, and inde-Multivariate analysis. We used a multivariate statistical approach to obtain an integrated view of the biological impact of the treatments. The fl owchart of multivariate statistical analysis is shown in Fig. 1 . The objective was to fi nd among selected biological processes ( Table 1 ) those processes affected mostly by both OEA and KDS administration and also associated with the reduction of body fat gain (see a detailed description below). This included linking of gene expression with the biochemical and physiological data that belong to the same string of biological functions across various systems and specifi cally grouped as lipid transport, EE, energy intake, endocannabinoid signaling, lipogenesis, and glucose metabolism ( Table 1 ). Such a statistical approach ( 17 ) allows a semiquantitation of the activity of each biological function under our nutritional challenge.
Due to the large number of biological variables compared with the low number of mice, partial least squares discriminant analysis (PLS-DA) and partial least-square regressions (PLS) were used to identify factors that signifi cantly contributed to the treatment effect. PLS-DA was used to select among the original 68 variables those that were treatment-responsive according to their variable importance index (VIP) (51 of 68 variables) ( Fig. 1 ). Concurrently, PLS was used to identify the variables that were related to adipose fat gain (36 of 68 variables). The variables in common (i.e., both treatment-sensitive and adipose fat-related variables determined independently) were retained (31 of 68 variables). These were reassigned to their respective biological modules. For each biological module, signifi cance of pathway enrichment was tested by comparing the PLS-DA model consisting of all the variables of a biological module to that calculated from the variables retained after the PLS-DA and PLS fi ltering ( Fig. 1 ). Signifi cant differences between the nonfi ltered and fi ltered models were performed using both Cooman's plot critical distance and ANOVA for the hierarchical PLS scores from nonfi ltered and fi ltered models (see supplementary material).
Calculations were performed using Simcap+ version 12.0 software (Umetrics, Umea, Sweden) using unit-of-variance scaling of the data. In each instance, model validation was performed using an R2 value (total variance explained) over 0.5, and a Q2 value (variance predicted after cross validation) at least equal to 2/3 of R2. Additional permutation tests (at least 100 permutations) allowed evaluation of model robustness. Variable selection was performed based on the signifi cant (95% jackknife confi dence interval) VIP with values over 0.8. All model validation indexes The combined contribution of the treatment-sensitive and adipose fat-related biological processes ( Table 1 ) explained 58% of the fat mass changes (PLS model, adipose tissue variance explained R2Y = 0.58; adipose tissue variance predicted after cross validation Q2Y = 0.52; decrease of R2Y and Q2Y after 100 permutations to 0.008 and Ϫ 0.14 at the intercept, respectively; P = 0.00021 after cross-validation ANOVA). Although the contribution of each function cannot be statistically compared to adipose fat mass, the R2 and VIP values (multivariate index of the pathway contribution to adipose fat mass) allowed us to rank their relationship to body fat mass variations ( Fig. 5 ). This contribution was, in decreasing order: regulation of lipid transport, then both energy intake and lipogenesis, and similarly endocannabinoid signaling, EE, and glucose metabolism.
We also performed a pair-wise correlation analysis to show the relationships among the biological events associated with body fat. The pathway correlation analysis ( Fig.  6 ) indicates that the lipid transport function was the most closely associated to adipose fat content. The comparison pendently increased LOX ( P = 0.0071) but did not affect GOX or EE compared with control.
We compared the extent of OEA biological degradation with that of its reportedly nonhydrolyzable analog KDS-5104 by using both in vitro and ex vivo assays. Using recombinant FAAH, we confi rmed that KDS-5104 is indeed stable under FAAH hydrolysis conditions, as its inhibition potential with FAAH activity remains unchanged after 2 h of preincubation in the presence of the enzyme. By comparison, for OEA, being a substrate of FAAH, 2 h of preincubation resulted in decreased inhibition of the enzyme ( Fig. 3A, B ). However, we found that the KDS-5104 inhibition potential was reduced after 2 h of preincubation with mouse tissue, suggesting it is hydrolyzed (data not shown). To confi rm this fi nding, we incubated KDS-5104 in the presence of mouse brain or liver homogenates for 2 h and found that 25% and 70% of the KDS-5104 was hydrolyzed by the tissue preparations, respectively ( Fig. 2D ). A similar experiment with OEA resulted in almost complete degradation of OEA ( Fig. 3C ).
Of note, when we looked at the intestinal FAAH activity of the treated mice, we noted that KDS-5104 supplementation reduced the FAAH activity by 57% compared with that of the vehicle-treated mice, whereas in OEA-fed mice, FAAH activity was decreased by 20% ( P = 0.0003 and 0.0422, respectively; see supplementary Table I).
A complementary separate experiment specifi cally examined fecal fat excretion. The level of fecal lipids was signifi cantly increased by OEA administration compared with that of control (a similar effect was observed for both the 0.1 and the 1 g/kg BW doses, P = 0.0042 and 0.0442, respectively) ( Fig. 4 ). An in vitro assay of pancreatic lipase activity showed a signifi cant decrease of the lipase hydrolytic activity after the addition of two different concentrations of OEA compared with that of control ( Ϫ 21.6% for 300 µmol of OEA/liter and Ϫ 25.3% for 600 µmol of OEA/ liter) ( P = 0.02, not shown). No signifi cant difference was observed between the two OEA concentrations.

DISCUSSION
We investigated which molecular and physiological processes were involved with the reduced body fat gain from oral OEA administration in a mouse model of diet-induced obesity. We also sought to determine the extent to which OEA could be broken down and how this would affect its bioeffi ciency by comparing its effects to those of OEA's more stable analog KDS-5104. KDS-5104 has a structure similar to OEA, and like OEA, effi ciently activates PPAR-␣ and induces a reduction in food intake ( 19 ).
Here, both OEA and KDS-5104 were similarly effective in signifi cantly reducing body fat gain in diet-induced obesity of C57Bl/6J male mice (16% to 20%). The biological processes involved in this outcome were also similar, although they occurred with some specifi city between OEA and KDS-5104. One may have predicted a somewhat higher effi ciency in the decrease of adipose fat for KDS-5104, based on its in vivo resistance to hydrolysis compared with that of OEA ( 19 ). However, although we found that KDS-5104 is resistant to FAAH degradation, we showed here that it is extensively broken down ex vivo (up to 70%) by using mouse liver homogenates. This implies that a system of the individual scores of each function indicated that both OEA and KDS-5104 modifi ed the regulation of the biological functions in the same direction, but the strength of the effect differed: the effect was equal for lipid transport and lipogenesis; there was a tendency to be higher (but not statistically different) in OEA fed-mice for energy intake and EE; and the effect was greater in KDS-5104 fedmice for endocannabinoid signaling and glucose metabolism. Finally, the individual PLS VIP values ( Fig. 6 ) showed which individual components of the biological functions studied were more likely affected by both OEA and KDS-5104 in relationship with adipose tissue expansion: downregulation of the epididymal adipose tissue FAT/CD36 gene expression in the lipid transport function, respectively; down-and upregulation of both liver SCD1 and FAS gene expression in the lipogenesis pathway; reduced cumulative food intake in the energy intake process; upregulation of the intestine CPT-I gene, along with increased blood ␤ -hydroxybutyrate and liver free glycerol; upregulation of the liver G6P gene expression in the glucose metabolism function; and upregulation of intestine FAAH enzyme activity and the GPR119 gene expression in the endocannabinoid-signaling module. Results are means ± SEM (N = 10 mice per group). Postprandial period lasting from feeding completion to 600 min afterward (from 396 to 1,008 min after the introduction of the mouse in the calorimetric system).
a For total activity, the P value for time effect was 0.0192, and the P value for treatment effect was 0.0039; time × treatment is nonsignifi cant. Total activity corresponds to the sum of rearing and spontaneous activity, normalized to the control mice mean value. response to treatment could be calculated ( 17 ). We determined that all of the biological modules evaluated were affected by the treatments and could explain up to 58% (PLS regression) of body fat variation. The part not captured by our model could occur for other reasons, as we measured only part of the biological components of the functions investigated. Instead of attempting to measure all functions, we chose to focus on dedicated processes that we wanted to assess as targets of either OEA or KDS-5104 ( Table 1 ). Other biologically relevant events could also be involved in body fat lowering, such as fat excretion. This was evaluated in a complementary study that other than FAAH is involved in catabolism of OEA. The important observation from this study is that the amount of OEA fed (100 mg/kg BW) is high enough to exert a biological effect despite a high degradation rate compared with that of KDS-5104. This is substantiated by our published results using two oral doses of OEA (10 mg and 100 mg/kg BW) and showing that both doses exerted the same effect on fat storage, despite a 10-fold intake difference ( 13 ). A dose extrapolation to the human, taking into account the metabolic rate difference with the mouse ( 20 ), would correspond to roughly 1 g/70 kg adult/day. If we also consider that OEA is effi ciently catabolized in vivo, this would suggest that it is active at very low levels in cells and that its level in cells is tightly controlled, thus emphasizing its biological importance.
Our study is the fi rst to investigate the impact of OEA on EE and its various components by indirect calorimetry. It revealed that OEA and KDS-5104 shifted the energy balance in favor of lipid oxidation. In our diagram of the EE pathway ( Fig. 6 ), the LOX component clustered with the intestine CPTI gene, a gene committed in fatty acid oxidation but also with metabolic indicators of a greater lipid substrate utilization (blood ␤ -hydroxybutyrate occurring from LOX [ 21 ]). Thus, the channeling of fat toward the oxidation pathway contributes to lower body fat expansion upon high-fat feeding. Lipid oxidation also clustered with spontaneous activity that is dramatically increased by OEA, thus OEA can contribute to nonexercise activity thermogenesis, a factor that also contributes to obesity resistance ( 22 ).
To obtain a simple picture of the biological effects associated with body fat gain, the various biological indicators were grouped as biological modules from which a  sion was upregulated by the treatments, as observed by others with OEA ( 2,13 ). It would be interesting to study whether this acts to signal the brain to inhibit overfeeding and trigger satiety. Finally, in the lipid transport cluster, we identifi ed the fact that the adipose tissue FAT/CD36 gene expression was the indicator most strongly associated with reduced body fat gain in response to OEA (and KDS-5104). This makes it a highly putative target for oral OEA action in the control of body fat.
Among the other most relevant targets found here, food intake control by OEA was also identifi ed as important to the control of body fatness. In fact, after correction for food intake, the differences in adipose fat were no longer signifi cant among treatment groups, thus emphasizing the important role of this factor in the antiobesity action of OEA and confi rming other similar fi ndings ( 6,13,23,24 ). In addition, downregulation of the stearoyl-CoA desaturase 1 (SCD1) gene activity in the liver was observed in both the OEA and KDS-5104 groups. SCD1 is a key enzyme that coordinates fat storage and fat oxidation ( 25 ). Decreased SCD1 expression has been demonstrated to be strongly associated with obesity resistance and increased lipid ␤ -oxidation ( 25,26 ) and also lowers food intake upon stearoylethanolamine administration ( 27 ). A similar fi nding was found in Zucker rats but only when 5 mg/kg injected OEA was combined with 1 mg/kg SR141716A (rimonabant) ( 28 ). We also found that downregulation of the SCD1 gene was one of the strongest factors (along with FAT/CD36 in adipose tissue) affecting cumulative food intake in a PLS regression model (not shown). The relationship found among the lipogenesis functions, the regulation of energy intake, and adipose tissue variations, as shown in Fig. 6, is also consistent with the above-described fi ndings and makes SCD1 a possible important target of OEA antiobesity action. In contrast to downregulation of SCD1, FAS expression, which is another de novo lipogenic confi rmed that intestinal fat absorption is decreased by OEA feeding. Although our in vitro results indicated a direct inhibitory effect of OEA on pancreatic lipase, it is not known if this mechanism is important under the conditions prevailing in vivo in the intestinal lumen. In our study, measured adiposity refl ected the cumulative exposure to either OEA or KDS-5104 over 5 weeks, whereas our snapshot of biological evaluations can be subject to time variations. Such nonlinearity can weaken our predictions. Nevertheless, we identifi ed important processes and some of their individual components targeted by either or both OEA and KDS-5104, which, in combination, are associated with decreased fat mass. For instance, from our model (R2 values in Fig. 5 ), we identifi ed the fact that the lipid transport function can predict most of the body fat variation (R2 = 0.60, e.g., 16% of adipose mass variation, P < 0.0001); followed by both energy intake and lipogenesis (R2 = 0.4, 11% of adipose mass variations, P = 0.0006; and R2 = 0.38, 10% of adipose mass variation, P = 0.0008, respectively); then endocannabinoid-signaling regulation (R2 = 0.30, 8% of adipose mass variation, P = 0.0039); EE (R2 = 0.26, 7.0% of adipose mass variation, P = 0.0075); and glucose metabolism and (R2 = 0.23, 6% of adipose mass variation, P = 0.0078). Additionally, the pathways analyses ( Fig. 6 ) focusing on pair-wise relationships indicated that the lipid transport process is the cornerstone function most strongly associating the other functions to adipose fat. This emphasizes its functional importance in the inhibition of body fat gain by OEA intake.
Although fecal fat cannot be modeled by our equation, as it was determined in a complementary study, its excretion also appears as an indirect expression of lipid transport in the intestine. This again underlines its importance in the control of body fatness by oral OEA. Paradoxically, in contradiction to greater fecal fat excretion, the upperintestine fatty acid transporter FAT/CD36 gene expres- Last, it is interesting that the reported nonhydrolyzable OEA analog KDS-5104 did not exert more pronounced effects on fat storage in our study. The differences in responses between KDS-5104 and OEA were mainly related to glucose metabolism and endocannabinoid signaling for which the effect of KDS-5104 was greater. The impact on adipose weight gain was relatively weak, as both glucose metabolism and endocannabinoid signaling were not the primary factors explaining body fatness in our model. In contrast, OEA and KDS-5104 exerted the same effects on lipid transport, energy intake, and lipogenesis, the three main modules related to body fat deposition. The similar effects of both KDS-5104 and OEA on these three main modules are likely to explain the identical outcome on adipose fat gain. enzyme strongly linked to body fat, was upregulated by both OEA and KDS-5104 in our study. Decreased SCD1 activity produces lower oleic acid, which is the building block of stored triglycerides ( 29 ), thus, decreased SCD1 activity would direct de novo synthesized fatty acids from FAS to oxidation instead of storage, provided lipid oxidation is effi cient ( 30 ). This is also substantiated by a decreased concentration of postprandial plasma triglyceride in OEA-(and KDS-5104-) treated mice, accompanied by higher LOX, plasma ␤ -hydroxybutyrate, and CPT-1 gene transcription in the intestine, and lower SCD1 gene transcription. Finally, our present fi ndings confi rmed that endocannabinoid signaling, mainly through both intestinal and adipose tissue FAAH, appears to be a signifi cant player in modulating OEA effects ( 13 ). Fig. 6. Correlation network among processes with activity modifi ed by both OEA and KDS-5104 treatments is shown. Each node depicts individual biological activity. Histograms represent means ± SEM of semiquantitative activity of biological functions according to treatments, based on the hierarchical PLS computation performed from the individual components, with letters a, b, and c the ANOVA signifi cance threshold ( P < 0.05). Node size represents the importance of the process to predict adipose tissue mass (based on VIP values of the PLS analysis). The edges represent the absolute Pearson coeffi cient values among the nodes. Edge thickness represents the strength of the correlation among connected nodes. The correlation network was built on the maximum stringency calculated from the minimum structure displaying all nodes. The signifi cance of the Pearson correlation coeffi cient was checked by false discovery rate with a q value of <0.05. This threshold corresponded to Pearson correlation values of >0.701. By convention, pathway score values of the control mice were always positive, thus emphasizing the contrasting effect of treatments. Individual pathway components are displayed as a heat map sorted according to hierarchical clustering (red box, relative increase; green box, relative decrease [see supplementary material for signifi cant levels across treatments). The asterisk indicates which individual components in the biological module weighs more in the biological score process and is more likely to predict adipose tissue fat pads mass (greatest VIP index of the PLS analysis). Adipose tissue mass is shown as means ± SEM of unit of variance. AT, adipose tissue; G, gene; B, biochemical; P, physiological; ep, epididymal; per, peritoneal; intest., upper intestine; muscle, gastrocnemius muscle; ES, endocannabinoid signaling; GM, glucose metabolism; LT, lipid transport; EI, energy intake; LG, lipogenesis.
In conclusion, we confi rmed that OEA feeding decreased body fat gain upon high-fat feeding in the mouse. We found that at least six different pathways (lipid transport, energy intake, regulation of EE, endocannabinoid signaling, lipogenesis, and glucose metabolism) can be involved in this effect and explained up to 58% of the body fat variation. Other biologically relevant effects such as fecal fat excretion and other as-yet unidentifi ed mechanisms also contributed to reduce body fat gain upon OEA feeding. Our study also pointed out that the physiological level of OEA degradation does not compromise its bioefficiency at the oral dose used. Seemingly, a degradation system other than the fatty acyl anandamide hydrolyze is able to process both OEA and its nonhydrolyzable analog KDS-5104.
Finally, our results highlight the multifactorial aspect of body fat changes that occur following OEA nutritional intervention and that take place at several levels of biological systems. Although we are aware that we captured only part of these changes, our multilevel integrated biological analysis brings new mechanical hypotheses to explain the effect of OEA on body fat control in the diet-induced obesity mouse model and allow ranking their effects in importance with regard to adipose tissue weight gain. With this perspective, we conclude that lipid transport appears as a determinant player of the OEA biological response, in which adipose tissue FAT/CD36 expression was found to be the most relevant bioindicator/effector of OEA action.