Melania Melis1†
Gianfranca Carta1†
Stefano Pintus2Paolo Pintus2Carla A. Piras2
Elisabetta Murru1
Claudia Manca1
Vincenzo Di Marzo3
Sebastiano Banni1*
Iole Tomassini Barbarossa1*- 1Department of Biomedical Sciences, Section of Physiology, University of Cagliari, Monserrato, Italy
- 2Department of Internal Medicine, Center for Metabolic Diseases, Azienda Ospedaliera G. Brotzu, Cagliari, Italy
- 3Endocannabinoid Research Group, Institute of Biomolecular Chemistry, Consiglio Nazionale delle Ricerche, Pozzuoli, Italy
The multifunctional CD36 scavenger receptor facilitates fatty acid (FA) uptake and oxidation and it has been involved in the pathophysiology related to dysfunctional FA metabolism. The common variant in the CD36 gene, rs1761667 (A/G), whose allele A is characterized by a reduced protein expression, has been associated with taste sensitivity to and preference for fat. We therefore aimed at evaluating whether the CD36 polymorphism may influence fatty acid metabolism and endocannabinoid biosynthesis in normal weight (NW) and obese (OB) subjects. Red blood cell (RBC) fatty acid composition, and plasma endocannabinoid levels were determined. In NW subjects with AA genotype was found a marked reduction of RBC saturated fatty acids and palmitic/linoleic ratio (PA/LA), considered as de novo lipogenesis (DNL) biomarkers. Remarkably, to the reduction of DNL biomarkers corresponded an increase of omega-6 index, an indirect marker of the impact on fatty acid metabolism of dietary omega-6 fatty acids, endocannabinoid levels and a higher waist/hip ratio. The presence of the G allele was instead associated with increased endocannabinoid plasma levels and a trend for increased waist/hip ratio in obese subjects, even though exhibited decreased BMI with respect to those with AA genotype. These data indicate that the CD36 polymorphism, rs1761667, leads to a distinct metabolic pattern in NW and in OB subjects. Therefore, their determination may be crucial in developing personalized therapeutic strategies for ameliorating dyslipidemia and other metabolic disorders.
Introduction
The multifunctional CD36 scavenger receptor is an 88kDa membrane glycosilated protein of 471 amino acid which has been originally identified in platelet membrane (Ikeda et al., 1989). CD36 is expressed on different cell types, such as adipocytes, myocytes, monocytes, macrophages, hepatocytes, vascular endothelial cells and intestinal enterocytes (Silverstein and Febbraio, 2009; Su and Abumrad, 2009). Topology of this receptor predicts an extracellular domain with one hydrophobic sequence, which may loop back into the membrane, and two shorts cytoplasmic (Su and Abumrad, 2009). In the last decades, several papers have been focused on the role of the scavenger receptor CD36 in facilitating the uptake and oxidation of fatty acid (FA) in humans and rodents and its involvement in the pathophysiological mechanisms associated with dysfunctional FA metabolism (Ma et al., 2004; Love-Gregory et al., 2008; Schwenk et al., 2008; Zhou et al., 2008). In addition to facilitate FA uptake, the membrane protein CD36 acts as a signaling molecule and a receptor for a wide range of ligands (Febbraio and Silverstein, 2007), binds to resident lipoproteins facilitating the cholesteryl ester uptake and enables the uptake of oxidized low/high-density lipoproteins and cholesterol (Calvo et al., 1998; Nassir et al., 2007; Thorne et al., 2007). As a consequence of its many ligands and functions, CD36 could impact on a wide range of dysmetabolic conditions associated with obesity, such as dyslipidaemia, insulin resistance, diabetes, inflammation, atherosclerosis and cancer (Hirano et al., 2003; Pascual et al., 2017).
CD36 plays a decisive role in the orosensory perception of dietary lipids in mammals since long-chain fatty acids appear to be mainly accountable for dietary fat taste perception in the oral cavity (Fukuwatari et al., 2003; Martin et al., 2011; Pepino et al., 2012; Melis et al., 2015; Besnard et al., 2016). CD36 is expressed also on taste bud cells of the circumvallate papillae where it has been demonstrated to mediate the perception of FA and the initiation of the cephalic phase (Laugerette et al., 2005; Sclafani et al., 2007; Martin et al., 2011). CD36 expression in circumvallate papillae has been found to be significantly decreased in high-fat diet-induced obese rats (Zhang et al., 2011). This suggested that decreased expression of CD36 diminishes sensitivity to fat taste, leading to increase fatty food intake as a compensatory reaction. On the other hand, absolute deficiency seems to decrease fat intake, probably as a post-ingestive effect of late absorption (Sclafani et al., 2007).
The presence of CD36 has been found in gustatory papillae in humans (Simons et al., 2011; Ozdener et al., 2014) as the main long-chain fatty acid receptor in taste bud cells, contributing to the orosensory perception of dietary lipid and fat preference (Keller et al., 2012; Pepino et al., 2012). Greater preference for, and increased eating of fatty foods have been documented in obese subjects (Stewart et al., 2010), which may reflect a reduction of oral and gastrointestinal FA sensitivity in obesity (Stewart et al., 2011).
The polymorphisms of CD36 gene, rs1761667, rs1527483, and rs3840546, have been associated with taste perception to and preference for fat (rs1761667, rs1527483), and obesity (rs3840546) (Keller et al., 2012; Pepino et al., 2012; Daoudi et al., 2015; Melis et al., 2015; Karmous et al., 2017). Several reports highlighted the effect of the polymorphism rs1761667 on CD36 protein expression levels (Ghosh et al., 2011; Love-Gregory et al., 2011), that could justify changes in oro-sensory perception of fats (Ma et al., 2004; Madden et al., 2008; Baranowski et al., 2011).
Recently, we have found that taste sensitivity for the bitter taste of 6-n-propylthiouracil (PROP), regarded as a general oral marker for chemosensory perception, food preference and BMI (Tepper, 2008; Padiglia et al., 2010; Tepper et al., 2011, 2014; Melis and Tomassini Barbarossa, 2017; Melis et al., 2017), affects lipid metabolism in normal weight (NW) (Tomassini Barbarossa et al., 2013) and obese (OB) (Carta et al., 2017) subjects. Peculiarly, in subjects with obesity, we showed an increased BMI and decreased levels of plasmatic endocannabinoids in subjects with extremely elevated PROP sensitivity, also named super-testers, vs. subjects with reduced PROP sensitivity (non-tasters). However, OB super-tasters showed lower circulating non-esterified fatty acid (NEFA) and retinol levels, which may indicate a more favorable lipid metabolism and body fat distribution than those of OB non-tasters. Since circulating levels of endocannabinoids are strongly associated to visceral obesity, less efficient mitochondrial function and ectopic fat deposition (Cristino et al., 2014), these results suggested that PROP taste sensitivity may predispose to metabolic changes and thereby body mass composition. Therefore, due to the relevant influence of the CD36 polymorphism rs1761667 on macronutrient preference and cellular lipid metabolism in humans, we aimed at evaluating whether it may impact fatty acid metabolism and endocannabinoid biosynthesis in NW and OB subjects.
Materials and Methods
Participants
126 Caucasian volunteers (46 M and 80 F) were recruited through public advertisements in the area of Cagliari, Italy. Subjects were divided into two groups based on their BMI: Normal weight subjects (NW) had a BMI ranging from 18 to 25 kg/m2 (N = 64; 21 M, 43 F) and subjects with obesity (OB) had a BMI ranging from 30 to 50 kg/m2 (N = 62; 25 M, 37 F).
Participants with major diseases (e.g., diabetes and kidney disease), pregnancy or lactation, food allergies and those who used medications that could affect lipid metabolism (e.g., steroids, antihistamines, and certain antidepressants) were excluded. All subjects were informed concerning the procedure and the purpose of the study. All participants provided a signed informed consent form. The study was performed according to the guidelines of the Declaration of Helsinki of 1975 (revised in 1983), and the procedures involving human participants were approved by the ethical committee of the Brotzu Institution. The present study has been registered at ClinicalTrials.gov (Identifer: NCT02729584).
Study Design
All subjects were requested to refrain from eating and drinking (except water) for at least 8 h prior to the sample withdrawal. They had to be in the test room at 8.00 a.m. In each subject, weight (kg) and height (m) were measured and the body mass index (BMI) (kg/m2) calculated, and waist and hip circumferences were determined. In each subject, a sample of blood (4-mL) was collected, promptly centrifuged and stored at −80°C until the analyses were completed as described below. Afterwards, a 2-mL sample of saliva was collected and transferred into an acid-washed polypropylene test tube. The saliva samples were stored at −80°C until molecular analyses were completed, as described below.
Molecular Analyses
DNA was extracted from each saliva sample by using the QIAamp ® DNA Mini Kit (QIAGEN Hilden, Germany) according to the manufacturer's instructions. Its concentration was assessed by measurements at an optical density of 260 nm with an Agilent Cary 60 Uv-Vis Spectrophotometer (Agilent, Palo Alto, CA).
Subjects were genotyped for the single nucleotide polymorphism (SNP), rs1761667 (G/A), of CD36, located at the −31118 promoter region of exon 1A.
A polymerase chain reaction (PCR) was used to amplify the CD36 region, including the polymorphism. The forward 5′-CAAAATCACAATCTATTCAAGACCA-3′ and reverse 5′-TTTTGGGAGAAATTCTGAAGAG-3′primers were synthesized by Invitrogen (50 nmol scale, desalted) (Europrim, Invitrogen Cambridge, UK). To genotype the CD36 rs1761667 polymorphism molecular analyses were carried out by PCR followed by analysis with restriction enzyme of the fragments obtained according to Banerjee et al. (2010), as briefly described below. A190-bp fragment was amplified with the primers and digested with the HhaI restriction enzyme (Thermo Scientific Inc, Waltham, USA). This restriction enzyme recognizes GCG∧C site and cut if incubated at 37°C overnight. Digestion products were separated by electrophoresis on a 2% agarose gel and the DNA bands were visualized by ethidium bromide staining and ultraviolet light to score the deletion. PCR 50 bp Low Ladder DNA was used as a molecular mass marker (Gene Ruler™ -Thermo Scientific).
Samples Preparation, Lipid, and Endocannabinoid Extraction
Blood samples were taken from the antecubital veins of volunteers, EDTA (10 IU/mL) was added to blood samples and promptly centrifuged at 2,000 g for 10 min, and the resulting plasma and erythrocytes were stored at −80°C until analysis.
Plasma or erythrocytes (red blood cells, RBC) samples were homogenized and extracted with a trichloromethane-methanol solution (2:1, v/v) (Folch et al., 1957) containing 2 μg of vitamin E and internal deuterated standards, anandamide (AEA), 2-arachidonoyl-monoacylglycerol (2-AG), palmitoylethanolamide (PEA) and oleoylethanolamide (OEA), for quantification by isotope dilution ([2H]8 AEA, [2H]5 2-AG, [2H]4 PEA, and [2H]4 OEA as deuterated isotopes; Cayman Chemical, MI, USA; Piscitelli et al., 2011).
Analysis of the Endocannabinoids and Their Congeners in Plasma
Aliquots from the lipid-containing organic phase, from plasma, were used for AEA, 2-AG, PEA, and OEA quantification by liquid chromatography–atmospheric pressure chemical ionization–mass spectrometry [1100 HPLC system (Agilent Technologies, Santa Clara, CA, USA) equipped with MS Detector 6110 single quadrupole] by selected ion monitoring at M+1 values for the four compounds and their deuterated homologs, as described previously (Piscitelli et al., 2011).
Analysis of Erythrocytes Fatty Acids
An aliquot of the lipid fraction, from erythrocyte, was mildly saponified as described by Banni et al. (1996) in order to obtain free fatty acids (FFA) for HPLC analysis. Separation of FFA was performed with Agilent 1100 HPLC system equipped with a diode array detector (Agilent, Palo Alto, CA), a C-18 Inertsil 5 ODS-2 Chrompack column (Chrompack International BV, The Netherlands), 5 μm particle size, 150 4.6 mm, was used with a mobile phase of CH3CN/H2O/CHCOOH (70/30/0.12, v/v/v) at a flow rate of 1.5 mL/min as described in Melis et al. (2001). Unsaturated FFA were detected at 200 nm; spectra (195–315 nm) of the eluate were obtained every 1.28 s and were electronically stored, HP Chemstation software (Agilent, Palo Alto, CA). These spectra were taken to validate the identification of the HPLC peaks (Angioni et al., 2002).
Because saturated fatty acid (SAFA) are transparent to UV, after derivatization, they were analyzed as fatty acid methyl esters by a gas chromatograph with flame ionization detector (FID) (Agilent, Model 6890, Palo Alto, CA) as described previously in Batetta et al. (2009). Different fatty acid metabolic parameters were calculated. The ratio of the concentration of PA/LA and PA/ARA were taken as an indirect measurement of de novo lipogenesis (Jacobs et al., 2015). Omega-6 index, taken as an index of the impact of dietary intake of omega-6 on erythrocyte concentrations of omega-6 highly polyunsaturated fatty acids (HPUFAs), was calculated as the ratio between the sum of the concentrations of omega-6 HPUFAs (those with 20 or more carbon atoms and three or more double bonds) and the sum of total fatty acids: (20:4n-6+20:3n-6+22:4n-6+22:5n-6)/(sum of total fatty acids), similarly to how the omega-3 index is calculated (Harris and Von Schacky, 2004).
Statistical Analyses
Genotype distribution and allele frequencies of the CD36 SNP between NW and OB was tested by the Fisher method. The data are expressed as the mean ± SEM. A two-way analysis of covariance (ANCOVA) was used to analyze the differences for all parameters associated to genotypes of the CD36 SNP in NW and OB, where age was used as a covariate. Post hoc comparisons were performed with the Fisher least significant difference (LSD) test, unless the assumption of homogeneity of variance was violated, in which case the Duncan's test was used. Statistical analyses were conducted using the software STATISTICA for WINDOWS (version 7; StatSoft Inc, Tulsa, OK, USA). P ≤ 0.05 were considered significant.
Results
Molecular analysis at the rs1761667 (A/G) polymorphism of CD36 identified in the group of NW seven subjects being homozygous AA, 39 heterozygous AG, and 18 homozygous GG, while in OB, 15 subjects were homozygous AA, 31 heterozygous AG, and 16 homozygous GG. No differences were found between NW and OB according to the genotype distribution and allele frequency of this polymorphism (χ2 < 2.7282; p > 0.2308). However, it is interesting to note that the number of subjects who had the AA genotype was more than two-fold in OB (24%) than in NW (10%) (χ2 = 3.384; p = 0.05).
The mean values ± SEM of BMI and the waist/hip ratio, in NW and OB subjects, according to genotypes of rs1761667 locus are shown in Figure 1. Post hoc comparisons subsequent to two-way ANCOVA showed that OB subjects homozygous for the G-allele exhibited a BMI lower than OB subjects being homozygous AA (p = 0.019, Fisher LSD test) (Figure 1A). An opposite trend, although not significant, was found in NW (Figure 1A). A pairwise comparison also showed that the waist/hip ratio of homozygous AA and NW subjects was significantly higher than that of NW subjects that were heterozygous or homozygous for the G-allele (p < 0.040, Duncan test), and not significantly different from that determined in OB subjects with the same genotype, who exhibited values not different from those of other genotype groups of OB subjects (Figure 1B). On the other hand, heterozygous or homozygous for the G-allele OB subjects still showed a higher waist/hip ratio than NW subjects with the same genotype (p < 0.000085, Duncan test) (Figure 1B).
Figure 1. BMI and waist/hip ratio according to the rs1761667 (A/G) polymorphism of CD36 in normal weight and obese subjects. Mean values ± SEM of body mass index (BMI) (A) and the waist/hip ratio (B) of homozygous AA at the polymorphism of CD36 (N = 7), heterozygous (N = 39) and homozygous GG (N = 18) subjects with normal weight (NW) and homozygous AA (N = 15), heterozygous (N = 31) and homozygous GG (N = 16) subjects with obesity (OB). The different letters indicate significant differences (BMI: p < 0.0419; waist/hip ratio: p < 0.040; Duncan test subsequent two-way ANCOVA). *a significant difference between corresponding NW or OB subjects (waist/hip ratio: p < 0.000085; Duncan test subsequent two-way ANCOVA).
The plasma levels of the endocannabinoids 2-AG and AEA were measured in NW and OB according to genotypes of rs1761667 locus and are presented in Figure 2. An ANCOVA revealed a significant two-way interaction of genotype group × NW/OB status for the 2-AG value [F(2, 119) = 4.1645; p = 0.0179] (Figure 2A). Post hoc comparisons revealed that 2-AG plasma levels of homozygous AA and NW subjects were significantly higher than those determined in NW subjects that were heterozygous or homozygous for the G-allele (p < 0.050, Fisher LSD test), and not significantly different from those of AA and OB subjects, who had values lower than those of homozygous GG OB subjects (p = 0.034, Fisher LSD test; Figure 2A). In addition, OB subjects heterozygous or homozygous for the G-allele showed higher 2-AG values than corresponding NW subjects (p < 0.000085, Fisher LSD test; Figure 2A). Post hoc comparisons also showed that the AEA plasma levels in homozygous AA and OB subjects were significantly lower than those determined in OB subjects who were heterozygous or homozygous for the G-allele (p < 0.050, Fisher LSD test; Figure 2B). An opposite trend, although not significant, was found in NW subjects (Figure 2B). OB subjects heterozygous or homozygous for the G-allele, showed higher AEA values than corresponding NW subjects (p < 0.00089, Fisher LSD test; Figure 2B), while the values of homozygous AA OB did not vary with respect to those of homozygous AA NW subjects (Figure 2B). No significant changes of endocannabinoid congeners OEA and PEA were found between NW and OB and rs1761667 genotype (data not shown).
Figure 2. 2-AG and AEA according to the rs1761667 (A/G) polymorphism of CD36 in normal weight and obese subjects. Mean values ± SEM of endocannabinoids (AEA and 2-AG) plasma levels measured in homozygous AA at the polymorphism of CD36 (N = 7), heterozygous (N = 39) and homozygous GG (N = 18) subjects with normal weight (NW) and homozygous AA (N = 15), heterozygous (N = 31) and homozygous GG (N = 16) subjects with obesity (OB). (A): 2-arachidonoylglycerol (2-AG) values. (B): anandamide (AEA) values. The different letters indicate significant differences (2-AG: p < 0.05; AEA: p < 0.05; Fisher LSD test subsequent two-way ANCOVA). *a significant difference between corresponding NW or OB subjects (2-AG: p < 0.0016; AEA: p < 0.00089; Fisher LSD test subsequent two-way ANCOVA).
The erythrocyte SAFA levels (mean values ± SEM), PA/LA ratio, palmitic/arachidonic acid (PA/ARA) ratio, and the omega-6 index measured in NW and OB according to genotypes of the rs1761667 locus are shown in Figure 3. A significant two-way interaction of genotype group × NW/OB status on SAFA values was found [F(2, 119) = 3.7075, p = 0.02752; two-way ANCOVA]. Post hoc comparisons showed that the SAFA values (Figure 3A), the PA/LA ratio (Figure 3B) and the PA/ARA ratio (Figure 3C) of NW homozygous for the A-allele were significantly lower than those of NW subjects that were heterozygous or homozygous GG (p < 0.029 and p < 0.031, respectively; Duncan test). Accordingly, the omega-6 index (Figure 3D) of homozygous for the A-allele NW subjects was significantly higher than that of heterozygous or homozygous GG NW subjects (p = 0.045 and p = 0.0069, respectively; Duncan test). No significant differences in SAFA values, PA/LA ratios, PA/ARA ratios or the omega-6 index related to rs1761667 locus were found in OB subjects (p > 0.05). The SAFA values, the PA/LA ratio and the PA/ARA ratio of homozygous for the A-allele NW subjects were significantly lower, while the omega-6 index was significantly higher, than those of the corresponding OB subjects (SAFA values: p = 0.000026; PA/LA: p = 0.000029; PA/ARA: p = 0.00037; Omega-6 index: p = 0.0020; Duncan test). The PA/LA ratio of OB subjects that were heterozygous or homozygous GG OB was higher than that of the corresponding NW subjects (p < 0.010; Duncan test). Concentrations of other fatty acids, as well related metabolic parameters such as omega-3 index, were not found significantly different between NW and OB and rs1761667 genotype (data not shown).
Figure 3. SAFA, PA/LA, PA/ARA, and omega-6 index according to the rs1761667 (A/G) polymorphism of CD36 in normal weight and obese subjects. Mean values ± SEM of levels of erythrocyte saturated fatty acid (SAFA), the palmitic/linoleic acids ratio (PA/LA), the palmitic/arachidonic acids ratio (PA/ARA), omega-6 index [(20:4n-6+20:3n-6+22:4n-6+22:5n-6)/(sum of total fatty acids)] measured in homozygous AA at the polymorphism of CD36 (N = 7), heterozygous (N = 39) and homozygous GG (N = 18) subjects with normal weight (NW) and homozygous AA (N = 15), heterozygous (N = 31) and homozygous GG (N = 16) subjects with obesity (OB). (A): erythrocyte SAFA values. (B): erythrocyte palmitic/linoleic acids ratio (PA/LA). (C): erythrocyte palmitic/arachidonic acids ratio (PA/ARA). (D): omega-6 index values. The different letters indicate significant differences (SAFA: p < 0.033; PA/LA: p < 0.031; PA/ARA: p < 0.0070; omega-6 index: p = 0.045 Duncan test). *a significant difference between corresponding NW or OB subjects (SAFA: p < 0.000026; PA/LA: p < 0.010; PA/ARA: p < 0.00037; omega-6 index: p = 0.0020 Duncan test).
Discussion
Our results confirm our previous reports (Tomassini Barbarossa et al., 2013; Carta et al., 2017) that genetic factors influencing dietary preferences affect fatty acid metabolism and circulating levels of endocannabinoids, differently in NW and OB subjects. Interestingly, in NW with the AA genotype of the rs1761667 polymorphism in CD36, characterized by a lower expression of CD36 protein with respect to subjects carrying the G allele, as previously shown (Ghosh et al., 2011; Love-Gregory et al., 2011), we found a marked reduction of RBC SAFA levels and PA/LA and PA/ARA ratios which have been used as markers of DNL (Jacobs et al., 2015). DNL is a metabolic pathway characterized by PA formation from glucose (Vernieri et al., 2016). Thereby biomarkers of DNL has also been taken as markers of carbohydrate intake (Jacobs et al., 2015). In addition, detection of these biomarkers in fatty acid profile of RBC has the advantage to represent long term metabolic changes because it is more steady than that of the plasma (Harris and Thomas, 2010). Therefore, we may hypothesize that NW individuals with AA genotype have a lower intake of carbohydrates compared to individuals NW carrying the G allele in the same locus.
Remarkably, to the reduction of DNL biomarkers corresponded an increase of the omega-6 index, which has been proposed as a marker of fat intake (Hudgins et al., 1996; Raatz et al., 2001; Chong et al., 2008; Jacobs et al., 2015). Indeed, individuals with the AA genotype show a preference for dietary fat (Keller et al., 2012). Thus, higher 2-AG levels in NW subjects with AA genotype may indicate that 2-AG biosynthesis is influenced by a high fat intake, particularly omega-6 fatty acids. In fact, 2-AG is ultimately derived from the omega-6 polyunsaturated fatty acid arachidonic acid and its biosynthetic precursor LA, and it has been shown in different experimental models (Matias et al., 2008; Batetta et al., 2009; Di Marzo et al., 2010; Piscitelli et al., 2011; Alvheim et al., 2012; Demizieux et al., 2016) and in humans (Banni et al., 2011), that changes in the precursor availability significantly affect endocannabinoids tissue levels. AEA circulating levels showed a similar trend, which did not reach statistical significance suggesting that other variables may influence the circulating levels of this endocannabinoid. Importantly, NW subjects with the AA genotype were also characterized by a higher waist/hip ratio, in agreement with data in the literature showing that in humans and experimental models, elevated endocannabinoid levels are associated to visceral fat deposition (Côté et al., 2007; Di Marzo et al., 2009). In obese subjects, however, AA genotype, was not associated to changes in fatty acid profile biomarkers of intake of carbohydrates or fat, in particular the omega-6 index. A possible explanation could be a deranged glucose disposal in obese subjects, irrespective of the genotype, which may enhance DNL, as revealed by the specific biomarker PA/LA. On the other hand, it is less intuitive that in obese subjects carrying the G allele again exhibit increased endocannabinoids plasma levels, with respect to those of OB subjects with the AA genotype. CD36 GG genotype is characterized by a highly expressed CD36 receptor which has a significant role in the regulation of fatty acid entry into the cell in physiological conditions. However, it has been shown that in conditions of an excess fatty acid entry, as it occurs in obese subjects with higher circulating levels of NEFA, CD36 activity is impaired, inducing ectopic fat deposition and decreasing mitochondrial activity (Pepino et al., 2014; Samovski et al., 2015) leading to a reduction of lean mass and increasing fat mass particularly in the visceral adipose tissue. This hypothesis is supported by the decrease of BMI, possibly due to a reduction in lean mass, which corresponded to a trend for increased waist/hip ratio in GG and GA OB subjects when compared to those carrying the AA genotype, which is also more frequent in OB subject than in NW. Thus, while in NW subjects with the AA genotype of CD36 the increased 2-AG plasma levels might be due at least in part to increased omega-6 intake and reduced DNL, in OB subjects with the GG or GA phenotype the excess visceral fat might become a source of endocannabinoids.
We may therefore conclude that a reduced function due to a low expression of CD 36 protein, either for genetic reason, as in the AA genotype, or due to an excess exposure to fatty acids as it occurs in metabolically dysfunctional OB subjects, may lead to impaired fat disposal at the cellular and tissue level and increased endocannabinoid signaling at cannabinoid CB1 receptors. The latter effect may predispose NW subjects carrying the AA allele of CD36 to become obese, and set into motion a feedforward vicious circle leading to more visceral fat accumulation in OB subjects carrying the GG allele. At any rate, our data indicate that in the obese state, dysmetabolic conditions influence CD36 activity differently according to its genotype. It has been shown that CD36 activity may also influence the metastatic properties of tumor cells (Pascual et al., 2017), implicated in the incorporation of oxidized LDL in macrophages during atherosclerotic plaque formation (Park et al., 2009) and has been associated to Alzheimer disease (Šerý et al., 2017). Conversely, in physiological conditions a functional CD36 regulates, in a concerted action with AMPK, fatty acid entry and mitochondrial beta-oxidation in the muscle and heart (Pepino et al., 2014; Samovski et al., 2015), even though excessive fatty acid exposure impairs AMPK activity (Samovski et al., 2015) and may thereby lead to metabolic derangement. This implies that CD36 plays a key role in cellular fatty acid metabolism in physiological and pathological conditions in different tissues. Indeed, CD36 knockout mice are resistant to obesity and dysmetabolic disorders (Clugston et al., 2014).
In conclusion, the A/G allele of the rs1761667 polymorphism of CD36 was found associated to a distinct metabolic pattern in NW and obese subjects. Therefore, the identification of the latter may be crucial in developing personalized therapeutic strategies for ameliorating dyslipidemia and other metabolic disorders.
Author Contributions
Designed the work: SB, IT; data acquisition: GC, MM, CM, EM, CP, SP, and PP; analysis and interpretation of data: MM, IT, SB, and GC; drafting the work: SB, IT; revising critically the work: SB, IT, MM, EM, CP, GC, and VD; final approval of the version to be published: SB, IT, MM, GC, CM, EM, CP, and VD; agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work were appropriately investigated and resolved: SB, IT, MM, GC, CM, EM, CP, and VD.
Funding
This work was supported by grant from the University of Cagliari (Fondo Integrativo per la Ricerca, FIR 2016).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank the volunteers, without whose contribution this study would not have been possible.
Abbreviations
2-AG, 2-arachidonoyl-monoacylglycerol; AEA, Anandamide; ARA, arachidonic acid; DNL, de novo lipogenesis; FA, fatty acid; FFA, free fatty acid; LA, linoleic acid; NEFA, non esterified fatty acid; NW, normal weight subjects; OB obese subjects; OEA, oleoylethanolamide; PA, palmitic acid; PEA, palmitoylethanolamide; PROP, 6-n-propylthiouracil; RBC, red blood cell; SAFA, saturated fatty acid.
References
Alvheim, A. R., Malde, M. K., Osei-Hyiaman, D., Lin, Y. H., Pawlosky, R. J., Madsen, L., et al. (2012). Dietary linoleic acid elevates endogenous 2-AG and anandamide and induces obesity. Obesity 20, 1984–1994. doi: 10.1038/oby.2012.38
Angioni, E., Lercker, G., Frega, N. G., Carta, G., Melis, M. P., Murru, E., et al. (2002). UV spectral properties of lipids as a tool for their identification. Eur. J. Lipid Sci. Techenol. 104.59–64. doi: 10.1002/1438-9312(200201)104:1<59::AID-EJLT59>3.0.CO;2-I
Banerjee, M., Gautam, S., Saxena, M., Bid, H. K., and Agrawal, C. G. (2010). Association of CD36 gene variants rs1761667 (G > A) and rs1527483 (C > T) with Type 2 diabetes in North Indian population. Int. J. Diab. Mellitus 2, 179–183. doi: 10.1016/j.ijdm.2010.08.002
Banni, S., Carta, G., Contini, M. S., Angioni, E., Deiana, M., Dessi, M. A., et al. (1996). Characterization of conjugated diene fatty acids in milk, dairy products, and lamb tissues. J. Nutr. Biochem. 7, 150–155. doi: 10.1016/0955-2863(95)00193-X
Banni, S., Carta, G., Murru, E., Cordeddu, L., Giordano, E., Sirigu, A. R., et al. (2011). Krill oil significantly decreases 2-arachidonoylglycerol plasma levels in obese subjects. Nutr. Metab. 8:7. doi: 10.1186/1743-7075-8-7
Baranowski, T., Baranowski, J. C., Watson, K. B., Jago, R., Islam, N., Beltran, A., et al. (2011). 6-n-propylthiouracil taster status not related to reported cruciferous vegetable intake among ethnically diverse children. Nutr. Res. 31, 594–600. doi: 10.1016/j.nutres.2011.07.004
Batetta, B., Griinari, M., Carta, G., Murru, E., Ligresti, A., Cordeddu, L., et al. (2009). Endocannabinoids may mediate the ability of (n-3) fatty acids to reduce ectopic fat and inflammatory mediators in obese Zucker rats. J. Nutr. 139, 1495–1501. doi: 10.3945/jn.109.104844
Besnard, P., Passilly-Degrace, P., and Khan, N. A. (2016). Taste of Fat: a Sixth Taste Modality? Physiol. Rev. 96, 151–176. doi: 10.1152/physrev.00002.2015
Calvo, D., Gomez-Coronado, D., Suarez, Y., Lasuncion, M. A., and Vega, M. A. (1998). Human CD36 is a high affinity receptor for the native lipoproteins HDL, LDL, and VLDL. J. Lipid Res. 39, 777–788.
Carta, G., Melis, M., Pintus, S., Pintus, P., Piras, C. A., Muredda, L., et al. (2017). Participants with normal weight or with obesity show different relationships of 6-n-propylthiouracil (PROP) taster status with bmi and plasma endocannabinoids. Sci. Rep. 7:1361. doi: 10.1038/s41598-017-01562-1
Chong, M. F., Hodson, L., Bickerton, A. S., Roberts, R., Neville, M., Karpe, F., et al. (2008). Parallel activation of de novo lipogenesis and stearoyl-CoA desaturase activity after 3 d of high-carbohydrate feeding. Am. J. Clin. Nutr. 87, 817–823.
Clugston, R. D., Yuen, J. J., Hu, Y., Abumrad, N. A., Berk, P. D., Goldberg, I. J., et al. (2014). CD36-deficient mice are resistant to alcohol- and high-carbohydrate-induced hepatic steatosis. J. Lipid Res. 55, 239–246. doi: 10.1194/jlr.M041863
Côté, M., Matias, I., Lemieux, I., Petrosino, S., Almeras, N., Despres, J. P., et al. (2007). Circulating endocannabinoid levels, abdominal adiposity and related cardiometabolic risk factors in obese men. Int. J. Obes. 31, 692–699. doi: 10.1038/sj.ijo.0803539
Cristino, L., Becker, T., and Di Marzo, V. (2014). Endocannabinoids and energy homeostasis: an update. Biofactors 40, 389–397. doi: 10.1002/biof.1168
Daoudi, H., Plesnik, J., Sayed, A., Sery, O., Rouabah, A., Rouabah, L., et al. (2015). Oral fat sensing and cd36 gene polymorphism in algerian lean and obese teenagers. Nutrients 7, 9096–9104. doi: 10.3390/nu7115455
Demizieux, L., Piscitelli, F., Troy-Fioramonti, S., Iannotti, F. A., Borrino, S., Gresti, J., et al. (2016). Early low-fat diet enriched with linolenic acid reduces liver endocannabinoid tone and improves late glycemic control after a high-fat diet challenge in mice. Diabetes 65, 1824–1837. doi: 10.2337/db15-1279
Di Marzo, V., Côté, M., Matias, I., Lemieux, I., Arsenault, B. J., Cartier, A., et al. (2009). Changes in plasma endocannabinoid levels in viscerally obese men following a 1 year lifestyle modification programme and waist circumference reduction: associations with changes in metabolic risk factors. Diabetologia 52, 213–217. doi: 10.1007/s00125-008-1178-6
Di Marzo, V., Griinari, M., Carta, G., Murru, E., Ligresti, A., Cordeddu, L., et al. (2010). Dietary krill oil increases docosahexaenoic acid and reduces 2-arachidonoylglycerol but not N-acylethanolamine levels in the brain of obese Zucker rats. Int. Dairy J. 20, 231–235. doi: 10.1016/j.idairyj.2009.11.015
Febbraio, M., and Silverstein, R. L. (2007). CD36: implications in cardiovascular disease. Int. J. Biochem. Cell Biol. 39, 2012–2030. doi: 10.1016/j.biocel.2007.03.012
Folch, J., Lees, M., and Sloane Stanley, G. H. (1957). A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem. 226, 497–509.
Fukuwatari, T., Shibata, K., Iguchi, K., Saeki, T., Iwata, A., Tani, K., et al. (2003). Role of gustation in the recognition of oleate and triolein in anosmic rats. Physiol. Behav. 78, 579–583. doi: 10.1016/S0031-9384(03)00037-4
Ghosh, A., Murugesan, G., Chen, K., Zhang, L., Wang, Q., Febbraio, M., et al. (2011). Platelet CD36 surface expression levels affect functional responses to oxidized LDL and are associated with inheritance of specific genetic polymorphisms. Blood 117, 6355–6366. doi: 10.1182/blood-2011-02-338582
Harris, W. S., and Thomas, R. M. (2010). Biological variability of blood omega-3 biomarkers. Clin. Biochem. 43, 338–340. doi: 10.1016/j.clinbiochem.2009.08.016
Harris, W. S., and Von Schacky, C. (2004). The Omega-3 Index: a new risk factor for death from coronary heart disease? Prev. Med. 39, 212–220. doi: 10.1016/j.ypmed.2004.02.030
Hirano, K., Kuwasako, T., Nakagawa-Toyama, Y., Janabi, M., Yamashita, S., and Matsuzawa, Y. (2003). Pathophysiology of human genetic CD36 deficiency. Trends Cardiovasc. Med. 13, 136–141. doi: 10.1016/S1050-1738(03)00026-4
Hudgins, L. C., Hellerstein, M., Seidman, C., Neese, R., Diakun, J., and Hirsch, J. (1996). Human fatty acid synthesis is stimulated by a eucaloric low fat, high carbohydrate diet. J. Clin. Invest. 97, 2081–2091. doi: 10.1172/JCI118645
Ikeda, H., Mitani, T., Ohnuma, M., Haga, H., Ohtzuka, S., Kato, T., et al. (1989). A new platelet-specific antigen, Naka, involved in the refractoriness of HLA-matched platelet transfusion. Vox Sang. 57, 213–217. doi: 10.1111/j.1423-0410.1989.tb00826.x
Jacobs, S., Jager, S., Jansen, E., Peter, A., Stefan, N., Boeing, H., et al. (2015). Associations of Erythrocyte Fatty Acids in the De Novo Lipogenesis Pathway with Proxies of Liver Fat Accumulation in the EPIC-Potsdam Study. PLoS ONE 10:e0127368. doi: 10.1371/journal.pone.0127368
Karmous, I., Plesnik, J., Khan, A. S., Sery, O., Abid, A., Mankai, A., et al. (2017). Orosensory detection of bitter in fat-taster healthy and obese participants: genetic polymorphism of CD36 and TAS2R38. Clin. Nutr. doi: 10.1016/j.clnu.2017.06.004. [Epub ahead of print].
Keller, K. L., Liang, L. C., Sakimura, J., May, D., van Belle, C., Breen, C., et al. (2012). Common variants in the CD36 gene are associated with oral fat perception, fat preferences, and obesity in African Americans. Obesity 20, 1066–1073. doi: 10.1038/oby.2011.374
Laugerette, F., Passilly-Degrace, P., Patris, B., Niot, I., Febbraio, M., Montmayeur, J. P., et al. (2005). CD36 involvement in orosensory detection of dietary lipids, spontaneous fat preference, and digestive secretions. J. Clin. Invest. 115, 3177–3184. doi: 10.1172/JCI25299
Love-Gregory, L., Sherva, R., Schappe, T., Qi, J. S., McCrea, J., Klein, S., et al. (2011). Common CD36 SNPs reduce protein expression and may contribute to a protective atherogenic profile. Hum. Mol. Genet. 20, 193–201. doi: 10.1093/hmg/ddq449
Love-Gregory, L., Sherva, R., Sun, L., Wasson, J., Schappe, T., Doria, A., et al. (2008). Variants in the CD36 gene associate with the metabolic syndrome and high-density lipoprotein cholesterol. Hum. Mol. Genet. 17, 1695–1704. doi: 10.1093/hmg/ddn060
Ma, X., Bacci, S., Mlynarski, W., Gottardo, L., Soccio, T., Menzaghi, C., et al. (2004). A common haplotype at the CD36 locus is associated with high free fatty acid levels and increased cardiovascular risk in Caucasians. Hum. Mol. Genet. 13, 2197–2205. doi: 10.1093/hmg/ddh233
Madden, J., Carrero, J. J., Brunner, A., Dastur, N., Shearman, C. P., Calder, P. C., et al. (2008). Polymorphisms in the CD36 gene modulate the ability of fish oil supplements to lower fasting plasma triacyl glycerol and raise HDL cholesterol concentrations in healthy middle-aged men. Prostaglandins Leukot. Essent. Fatty Acids 78, 327–335. doi: 10.1016/j.plefa.2008.04.003
Martin, C., Chevrot, M., Poirier, H., Passilly-Degrace, P., Niot, I., and Besnard, P. (2011). CD36 as a lipid sensor. Physiol. Behav. 105, 36–42. doi: 10.1016/j.physbeh.2011.02.029
Matias, I., Carta, G., Murru, E., Petrosino, S., Banni, S., and Di Marzo, V. (2008). Effect of polyunsaturated fatty acids on endocannabinoid and N-acyl-ethanolamine levels in mouse adipocytes. Biochim. Biophys. Acta. 1781, 52–60. doi: 10.1016/j.bbalip.2007.11.001
Melis, M. P., Angioni, E., Carta, G., Murru, E., Scanu, P., Spada, S., et al. (2001). Characterization of conjugated linoleic acid and its metabolites by RP-HPLC with diode array detector. Eur. J. Lipid Sci. Technol. 103, 617–621. doi: 10.1002/1438-9312(200109)103:9<617::AID-EJLT6170>3.0.CO;2-C
Melis, M., Sollai, G., Muroni, P., Crnjar, R., and Barbarossa, I. T. (2015). Associations between orosensory perception of oleic acid, the common single nucleotide polymorphisms (rs1761667 and rs1527483) in the CD36 gene, and 6-n-propylthiouracil (PROP) tasting. Nutrients 7, 2068–2084. doi: 10.3390/nu7032068
Melis, M., and Tomassini Barbarossa, I. (2017). Taste perception of sweet, sour, salty, bitter, and umami and changes due to l-Arginine supplementation, as a function of genetic ability to taste 6-n-Propylthiouracil. Nutrients 9:541. doi: 10.3390/nu9060541
Melis, M., Yousaf, N. Y., Mattes, M. Z., Cabras, T., Messana, I., Crnjar, R., et al. (2017). Sensory perception of and salivary protein response to astringency as a function of the 6-n-propylthioural (PROP) bitter-taste phenotype. Physiol. Behav. 173, 163–173. doi: 10.1016/j.physbeh.2017.01.031
Nassir, F., Wilson, B., Han, X., Gross, R. W., and Abumrad, N. A. (2007). CD36 is important for fatty acid and cholesterol uptake by the proximal but not distal intestine. J. Biol. Chem. 282, 19493–19501. doi: 10.1074/jbc.M703330200
Ozdener, M. H., Subramaniam, S., Sundaresan, S., Sery, O., Hashimoto, T., Asakawa, Y., et al. (2014). CD36- and GPR120-mediated Ca(2)(+) signaling in human taste bud cells mediates differential responses to fatty acids and is altered in obese mice. Gastroenterology 146, 995–1005. doi: 10.1053/j.gastro.2014.01.006
Padiglia, A., Zonza, A., Atzori, E., Chillotti, C., Cal,ò, C., Tepper, B. J., et al. (2010). Sensitivity to 6-n-propylthiouracil is associated with gustin (carbonic anhydrase VI) gene polymorphism, salivary zinc, and body mass index in humans. Am. J. Clin. Nutr. 92, 539–545. doi: 10.3945/ajcn.2010.29418
Park, Y. M., Febbraio, M., and Silverstein, R. L. (2009). CD36 modulates migration of mouse and human macrophages in response to oxidized LDL and may contribute to macrophage trapping in the arterial intima. J. Clin. Invest. 119, 136–145. doi: 10.1172/JCI35535
Pascual, G., Avgustinova, A., Mejetta, S., Martin, M., Castellanos, A., Attolini, C. S., et al. (2017). Targeting metastasis-initiating cells through the fatty acid receptor CD36. Nature 541, 41–45. doi: 10.1038/nature20791
Pepino, M. Y., Kuda, O., Samovski, D., and Abumrad, N. A. (2014). Structure-function of CD36 and importance of fatty acid signal transduction in fat metabolism. Annu. Rev. Nutr. 34, 281–303. doi: 10.1146/annurev-nutr-071812-161220
Pepino, M. Y., Love-Gregory, L., Klein, S., and Abumrad, N. A. (2012). The fatty acid translocase gene CD36 and lingual lipase influence oral sensitivity to fat in obese subjects. J. Lipid Res. 53, 561–566. doi: 10.1194/jlr.M021873
Piscitelli, F., Carta, G., Bisogno, T., Murru, E., Cordeddu, L., Berge, K., et al. (2011). Effect of dietary krill oil supplementation on the endocannabinoidome of metabolically relevant tissues from high-fat-fed mice. Nutr. Metab. 8:51. doi: 10.1186/1743-7075-8-51
Raatz, S. K., Bibus, D., Thomas, W., and Kris-Etherton, P. (2001). Total fat intake modifies plasma fatty acid composition in humans. J. Nutr. 131, 231–234.
Samovski, D., Sun, J., Pietka, T., Gross, R. W., Eckel, R. H., Su, X., et al. (2015). Regulation of AMPK activation by CD36 links fatty acid uptake to beta-oxidation. Diabetes 64, 353–359. doi: 10.2337/db14-0582
Schwenk, R. W., Luiken, J. J., Bonen, A., and Glatz, J. F. (2008). Regulation of sarcolemmal glucose and fatty acid transporters in cardiac disease. Cardiovasc. Res. 79, 249–258. doi: 10.1093/cvr/cvn116
Sclafani, A., Ackroff, K., and Abumrad, N. A. (2007). CD36 gene deletion reduces fat preference and intake but not post-oral fat conditioning in mice. Am. J. Physiol. Regul. Integr. Comp. Physiol. 293, R1823–1832. doi: 10.1152/ajpregu.00211.2007
Šerý, O., Janoutova, J., Ewerlingova, L., Halova, A., Lochman, J., Janout, V., et al. (2017). CD36 gene polymorphism is associated with Alzheimer's disease. Biochimie 135, 46–53. doi: 10.1016/j.biochi.2017.01.009
Silverstein, R. L., and Febbraio, M. (2009). CD36, a scavenger receptor involved in immunity, metabolism, angiogenesis, and behavior. Sci. Signal. 2:re3. doi: 10.1126/scisignal.272re3
Simons, P. J., Kummer, J. A., Luiken, J. J., and Boon, L. (2011). Apical CD36 immunolocalization in human and porcine taste buds from circumvallate and foliate papillae. Acta Histochem. 113, 839–843. doi: 10.1016/j.acthis.2010.08.006
Stewart, J. E., Feinle-Bisset, C., Golding, M., Delahunty, C., Clifton, P. M., and Keast, R. S. (2010). Oral sensitivity to fatty acids, food consumption and BMI in human subjects. Br. J. Nutr. 104, 145–152. doi: 10.1017/S0007114510000267
Stewart, J. E., Seimon, R. V., Otto, B., Keast, R. S., Clifton, P. M., and Feinle-Bisset, C. (2011). Marked differences in gustatory and gastrointestinal sensitivity to oleic acid between lean and obese men. Am. J. Clin. Nutr. 93, 703–711. doi: 10.3945/ajcn.110.007583
Su, X., and Abumrad, N. A. (2009). Cellular fatty acid uptake: a pathway under construction. Trends Endocrinol. Metab. 20, 72–77. doi: 10.1016/j.tem.2008.11.001
Tepper, B. J. (2008). Nutritional implications of genetic taste variation: the role of PROP sensitivity and other taste phenotypes. Annu. Rev. Nutr. 28, 367–388. doi: 10.1146/annurev.nutr.28.061807.155458
Tepper, B. J., Banni, S., Melis, M., Crnjar, R., and Tomassini Barbarossa, I. (2014). Genetic Sensitivity to the Bitter Taste of 6-n-Propylthiouracil (PROP) and Its Association with Physiological Mechanisms Controlling Body Mass Index (BMI). Nutrients 6, 3363–3381. doi: 10.3390/nu6093363
Tepper, B. J., Neilland, M., Ullrich, N. V., Koelliker, Y., and Belzer, L. M. (2011). Greater energy intake from a buffet meal in lean, young women is associated with the 6-n-propylthiouracil (PROP) non-taster phenotype. Appetite 56, 104–110. doi: 10.1016/j.appet.2010.11.144
Thorne, R. F., Mhaidat, N. M., Ralston, K. J., and Burns, G. F. (2007). CD36 is a receptor for oxidized high density lipoprotein: implications for the development of atherosclerosis. FEBS Lett. 581, 1227–1232. doi: 10.1016/j.febslet.2007.02.043
Tomassini Barbarossa, I., Carta, G., Murru, E., Melis, M., Zonza, A., Vacca, C., et al. (2013). Taste sensitivity to 6-n-propylthiouracil is associated with endocannabinoid plasma levels in normal-weight individuals. Nutrition 29, 531–536. doi: 10.1016/j.nut.2012.09.018
Vernieri, C., Casola, S., Foiani, M., Pietrantonio, F., de Braud, F., and Longo, V. (2016). Targeting cancer metabolism: dietary and pharmacologic interventions. Cancer Discov. 6, 1315–1333. doi: 10.1158/2159-8290.CD-16-0615
Zhang, X. J., Zhou, L. H., Ban, X., Liu, D. X., Jiang, W., and Liu, X. M. (2011). Decreased expression of CD36 in circumvallate taste buds of high-fat diet induced obese rats. Acta Histochem. 113, 663–667. doi: 10.1016/j.acthis.2010.09.007
Keywords: CD36 gene, obesity, fatty acids, metabolism, endocannabinoids
Citation: Melis M, Carta G, Pintus S, Pintus P, Piras CA, Murru E, Manca C, Di Marzo V, Banni S and Tomassini Barbarossa I (2017) Polymorphism rs1761667 in the CD36 Gene Is Associated to Changes in Fatty Acid Metabolism and Circulating Endocannabinoid Levels Distinctively in Normal Weight and Obese Subjects. Front. Physiol. 8:1006. doi: 10.3389/fphys.2017.01006
Received: 27 September 2017; Accepted: 21 November 2017;
Published: 06 December 2017.
Edited by:
Naim Akhtar Khan, Université de Bourgogne, FranceCopyright © 2017 Melis, Carta, Pintus, Pintus, Piras, Murru, Manca, Di Marzo, Banni and Tomassini Barbarossa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sebastiano Banni, banni@unica.it
Iole Tomassini Barbarossa, tomassin@unica.it
†These authors have contributed equally to this work.