- Open Access
The beta-3 adrenergic agonist (CL-316,243) restores the expression of down-regulated fatty acid oxidation genes in type 2 diabetic mice
Nutrition & Metabolismvolume 12, Article number: 8 (2015)
The hallmark of Type 2 diabetes (T2D) is hyperglycemia, although there are multiple other metabolic abnormalities that occur with T2D, including insulin resistance and dyslipidemia. To advance T2D prevention and develop targeted therapies for its treatment, a greater understanding of the alterations in metabolic tissues associated with T2D is necessary. The aim of this study was to use microarray analysis of gene expression in metabolic tissues from a mouse model of pre-diabetes and T2D to further understand the metabolic abnormalities that may contribute to T2D. We also aimed to uncover the novel genes and pathways regulated by the insulin sensitizing agent (CL-316,243) to identify key pathways and target genes in metabolic tissues that can reverse the diabetic phenotype.
Male MKR mice on an FVB/n background and age matched wild-type (WT) FVB/n mice were used in all experiments. Skeletal muscle, liver and fat were isolated from prediabetic (3 week old) and diabetic (8 week old) MKR mice. Male MKR mice were treated with CL-316,243. Skeletal muscle, liver and fat were isolated after the treatment period. RNA was isolated from the metabolic tissues and subjected to microarray and KEGG database analysis.
Significant decreases in the expression of mitochondrial and peroxisomal fatty acid oxidation genes were found in the skeletal muscle and adipose tissue of adult MKR mice, and the liver of pre-diabetic MKR mice, compared to WT controls. After treatment with CL-316,243, the circulating glucose and insulin concentrations in the MKR mice improved, an increase in the expression of peroxisomal fatty acid oxidation genes was observed in addition to a decrease in the expression of retinaldehyde dehydrogenases. These genes were not previously known to be regulated by CL-316,243 treatment.
This study uncovers novel genes that may contribute to pharmacological reversal of insulin resistance and T2D and may be targets for treatment. In addition, it explains the lower free fatty acid levels in MKR mice after treatment with CL-316,243 and furthermore, it provides biomarker genes such as ACAA1 and HSD17b4 which could be further probed in a future study.
The global prevalence of diabetes is rising . As a result of population ageing, increasing rates of obesity and a sedentary lifestyle, there is an increasing incidence of Type 2 diabetes (T2D), which comprises 90% of diabetes cases worldwide [1,2]. There is an interplay between genetic susceptibility and environmental influences in this epidemic , and studies such as the Diabetes Prevention Program (DPP) have demonstrated that lifestyle intervention, particularly weight loss, leads to a significant reduction in the risk of diabetes [4,5]. It is generally believed that insulin resistance in metabolic tissues (liver, fat, and skeletal muscle) is a key factor in the development of T2D . What causes the development of insulin resistance in metabolic tissues is unclear. However, reducing insulin resistance by using medications, such as thiazolidinediones, decreases the incidence of diabetes [5,7]. Beta-3 agonists improve insulin resistance in animals, and short-term human studies, with evidence of browning of white adipose tissue in rodent studies [8,9]; but the mechanisms through which this occurs are incompletely understood.
Inter-tissue cross talk between skeletal muscle, adipose tissue, the liver, brain, and pancreatic beta cells may contribute to insulin resistance. Increasing circulating levels of free fatty acids (FFAs) and triglycerides (TG) from adipose tissue lipolysis are frequently observed in pre-diabetes and T2D and contribute to insulin resistance. Abnormalities in fatty acid metabolism contribute to the accumulation of lipids in tissues such as skeletal muscle and liver and to the development or worsening of insulin resistance . Other molecules, such as myokines and adipokines may contribute to inter-tissue cross talk. Skeletal muscle, for example releases a number of myokines that can act in an autocrine, paracrine, or endocrine manner to regulate metabolic processes [11,12]. Adipose tissue releases multiple adipokines, including leptin, adiponectin and resistin that act through their receptors on many other organs to regulate metabolism.
In this study we have used an animal model of insulin resistance and T2D to understand the gene expression changes that occur in pre-diabetes, T2D, and pharmacological treatment of T2D. The LeRoith lab developed a transgenic mouse model of T2D that overexpresses a dominant-negative IGF-1R (KR-IGF-1R) specifically in the skeletal muscle under the muscle creatine kinase (MCK) promoter . The male MKR mice exhibit hyperinsulinemia by 2–3 weeks of age, and hyperglycemia by 5–6 weeks of age with decreased whole body glucose uptake, failure to suppress hepatic gluconeogenesis in response to insulin and pancreatic beta cell dysfunction [13,14]. In addition, the male MKR mice have elevated circulating free fatty acids, triglycerides, and hepatic and muscle triglyceride content, compared to wild-type mice [13,15]. We chose the MKR mouse for these studies as this mouse model has a well-characterized progression from “pre-diabetes” to overt diabetes with hyperglycemia, without the need for high fat diet feeding. Therefore, this allows us to study the effects of insulin resistance in different tissues without the confounding effects of different diets with different fat contents. Therefore, the male MKR mouse is an excellent mouse model for studying the metabolic derangements associated with the pre-diabetic insulin resistant state and the diabetic hyperglycemic condition.
Previous gene expression studies on T2D have focused on single metabolic tissues [16-22], have studied target genes or proteins [13,23-27], or have been performed in vitro in cell culture studies [28-31]. None of these studies have explored the global gene expression changes in multiple tissues, or those caused by the administration of the insulin sensitizing beta-3 adrenergic agonist (CL–316, 243). In this study, we aimed to uncover novel changes in metabolic tissues (skeletal muscle, liver and adipose tissue) between pre-diabetic and diabetic MKR mice, using microarrays and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database analyses. The KEGG database is a comprehensive database constructed from well-known molecular interaction networks and is extensively used to study biological pathways . The enrichment of KEGG pathways was used to encode all significantly differentially expressed genes in this study. The study of extracted KEGG pathways related to T2D, indicate that they may help building effective computational tools in the study of T2D. In addition to examining differences in the metabolic tissues in the pre-diabetic and diabetic models, we also aimed to uncover novel genes and pathways that were altered by the pharmacological treatment with a CL-316,243 . Using the same methods, we uncovered novel genes and pathways that could be targets for future therapeutics for insulin resistance and T2D.
All animal studies were approved by the Mount Sinai School of Medicine Institutional Animal Care and Use Committee. Mice were housed in The Mount Sinai School of Medicine Center for Comparative Medicine and Surgery, Association for Assessment and Accreditation of Laboratory Animal Care International (AALAC) and Office of Laboratory Animal Welfare (OLAW) accredited facility, where animal care and maintenance were provided. Mice were kept on a 12 hour light/dark cycle, had free access to diet (Picolab Rodent Diet 20, 5053) and fresh water. All MKR and WT mice used in these studies were male, on the FVB/N background and were 3–12 weeks of age. The generation and characterization of the MKR mice have been previously described . All tissues were collected from mice in the fed state. Three animals per group were used in each microarray experiment. Liver, gonadal fat and skeletal muscle (quadriceps) were collected and flash frozen in liquid nitrogen for subsequent RNA extraction.
Treatment with CL-316,243
We have previously demonstrated that the beta-3 agonist CL-316,243 improves insulin sensitivity and lowers glucose levels in the MKR mice [26,33]. Nine to ten week old male WT and MKR mice were injected intraperitoneally with CL-316,243 (1 mg/kg BW/day) or with an equivalent volume of vehicle (sterile water and phosphate buffered saline) for three weeks, n = 5-8 per group. Body weight was measured before treatment and twice weekly during treatment. Body composition analysis was performed using the EchoMRI 3-in-1 NMR system (Echo Medical Systems, Houston, TX, USA) before treatment, and 13 days after the start of treatment. Fed blood glucose measurements were performed on tail vein blood during tumor studies using a Bayer Contour Glucometer (Bayer Healthcare, Mishawaka, IN, USA), prior to commencing treatment and weekly thereafter. Plasma fed insulin levels were measured at the end of the studies using the Sensitive Rat Insulin RIA kit (Millipore, St. Charles, MO, USA). Serum triglycerides were measured using the Point Scientific Liquid Reagent (Pointe Scientific, Canton, MI, USA). Means for each group and standard error of the means were calculated. Differences between groups were analyzed by a One-Way ANOVA Test with a Tukey-Kramer Post-Hoc Test, assuming a significance level of 5% using SigmaStat 3.5 (Systat Software Inc, San Jose, CA, USA).
RNA isolation and microarray analysis procedures
Total RNA was isolated from the liver, fat, and skeletal muscle tissues of three mice from each group using the RNeasy Microarray Tissue Mini kit (Qiagen, Valencia, CA, USA), according to the manufacturer’s instructions. The concentration and quality of RNA was determined using the NanoDrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA), the Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA), and RNA Integrity Number (RIN, Agilent, Santa Clara, CA, USA). All samples used for reverse transcription and microarray analysis had 260/280 ratios greater than 1.8 on the Nanodrop Spectrophotometer and RIN values greater than 8.0. 100 ng of RNA from each sample was reverse transcribed to cDNA and amplified using NUGEN Applause WT–Amp ST system (NuGEN Technologies, San Carlos, CA, USA), according to the manufacturer’s protocol. 2.5 μg of cDNA was fragmented and 3′-biotinylated using Encore Biotin Module (NuGEN Technologies, San Carlos, CA, USA). The resultant sample mix with hybridization reagents were loaded into the GeneChip Mouse Gene 1.0 ST array and incubated for 16–20 hours in hybridization oven rotating at 60 rpm at 45°C (Affymetrix, Santa Clara, CA, USA). Arrays were processed using GeneChip Fluidics Station 450 (Affymetrix, Santa Clara, CA, USA). Chips were scanned using the GeneChip Scanner 3000 (Affymetrix, Santa Clara, CA, USA), operated by Gene Chip Operating Software, version 1. 4 and generated .CEL, .CHP and RPT files. Poly-A controls (dap, lys, phe, thr) and hybridization controls (bioB, bioC, bioD and cre) were used as spike controls for cDNA synthesis and hybridization, respectively, using methods described in the manufacturer’s instructions (Affymetrix, Santa Clara, CA, USA).
The microarray raw data was analyzed using Partek software, version 6.3 (Partek. St. Louis, MO, USA). Raw data was subjected to Robust Multichip Average (RMA) quantile normalization to remove biases introduced by technical and experimental effects. All expression data were log base 2 -transformed to get near normal distribution for accurate statistical inference. Quality control by visualizing the data using Principal Component Analysis cluster plot ensured that no outliers were included for the analysis. Next, two-way ANOVA analysis was performed to obtain a set of differentially expressed genes. A filter of p value < 0.05 and fold-change > 1.5 times was applied to identify the significantly differentially expressed genes.
The significantly differentially expressed genes list was exported to Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City, CA) for finding biological inference. Statistically significant genes from Partek analysis were overlaid on the Ingenuity Pathway Analysis (IPA) global molecular network after applying a filter on species (mouse) and tissue type.
KEGG pathways analysis
Kyoto Encyclopedia of Genes and Genomes (KEGG) database pathways and their corresponding genes were downloaded from KEGG website (http://www.genome.jp/kegg/) on 11 June 2014 [34,35]. Differentially expressed genes for each study type, i.e. MKR vs WT (fat, skeletal muscle, and liver tissues) and CL 316,243 treated MKR vs vehicle treated MKR fat tissues and genes from entire mouse genome were mapped to each KEGG pathway. This provided a list of count of genes in each of the above four datasets and the entire mouse genome annotated to each KEGG pathway. Calculation of enrichment P values is based on the following mathematical approach: in a dataset of N genes from entire genome, where M genes are annotated to a particular KEGG pathway, and n genes are differentially expressed, then the probability of having k genes to be differentially expressed and also included in the above KEGG pathway is given by a hypergeometric distribution as described the equation below 
Based on this approach, calculations of enrichment P values and programming tasks were performed using MATLAB version 2010a [Natick, Massachusetts: The MathWorks Inc., 2010]. Enrichment P values were calculated using MATLAB’s hygecdf and hygepdf functions. The enrichment P value cutoff for KEGG analysis was 0.0001 based on Bonferroni correction applied for testing of multiple pathways.
KEGG pathways analysis of the microarray results
The results of the microarray analysis are summarized in Venn diagrams shown in Figure 1a and 1b. KEGG pathways analysis was done to determine overrepresented pathways (containing high number of differentially expressed genes). The results of KEGG pathway analysis of lipid metabolism is shown in Table 1. The following pathways were found to be significantly overrepresented in fat tissues of MKR vs WT mice: fatty acid degradation, glycerolipid metabolism, and glycerophospholipid metabolism. When fat tissues of CL-316,243 treated MKR mice were compared with non-treated MKR mice, the glycerolipid metabolism, glycerophospholipid metabolism, fat digestion and absorption, and ether lipid metabolism were found to be significantly overrepresented. Circulating ether lipid levels have previously found to be significantly higher in obese and diabetic subjects . As shown in Figure 1c when fat tissues of treated MKR were compared with non-treated MKR difference in dysregulation of ether lipid metabolism was observed. Most of the genes in fat tissues related to ether lipid metabolism were found to be down-regulated in MKR compared with WT mice. Following treatment this trend was reversed. A complete list of KEGG enrichment analysis results is provided in Additional file 1: Table S1. In addition, pathways related to carbohydrate metabolism such as TCA cycle, pentose phosphate pathway, fructose and mannose metabolism, pyruvate metabolism, and propanoate metabolism were found to be overrepresented only in fat tissues of MKR vs WT mice and not in the fat tissues of CL-316,243 treated MKR compared with non-treated MKR. Moreover, carbohydrate metabolism pathway amino sugar and nucleotide sugar metabolism was found to be overrepresented only after treatment with CL-316,243 in fat tissues of MKR mice.
Fatty acid oxidation genes are dysregulated in the adipose tissues of the diabetic MKR mice compared with WT mice
Although the diabetic MKR mice were generated by overexpression of a dominant-negative human IGF-1R in skeletal muscle, it has previously been found that these mice also develop insulin resistance in liver and adipose tissue as well as a decrease in adipose tissue mass. Pathway analysis of the gene expression in the adipose tissue of MKR and WT mice by microarrays revealed dysregulated expression of the TCA cycle, branched chain amino acid (valine) degradation, tryptophan degradation, N-glycan biosynthesis, pyrimidine metabolism, ether lipid metabolism, and fatty acid oxidation. As the MKR mice have been previously found to have changes in fatty acid metabolism in the skeletal muscle, with increased adipose tissue and hepatic insulin resistance, we further investigated the specific alterations in the expression of genes in the fatty acid oxidation pathway in the adipose tissue of the MKR mice, compared to WT mice. We found down regulation of number of genes involved in fatty acid oxidation in adipose tissue (Figure 2, Table 2), with a few notable exceptions, including the aldehyde dehydrogenases ALDH3A2 and ALDH7A1, the acyl CoA synthetases ACSL4 and ACSBG2 and the peroxisomal fatty acid oxidase ACOX3, that were upregulated in the adipose tissue of the MKR mice.
Insulin resistance leads to downregulation of fatty acid oxidation genes in the skeletal muscle of MKR mice
We next examined the expression of fatty acid oxidation genes in the skeletal muscle of the MKR and WT mice by microarray analysis. Significant differences were found in the fatty acid oxidation pathway between MKR and WT mice. The MKR mice had a significant downregulation of many of the genes involved in fatty acid beta-oxidation (Figure 3, Table 3). Notable exceptions to the general downregulation of fatty acid oxidation genes in the skeletal muscle of the MKR mice include carnitine palmitoyl transferase 1a (CPT1A), the peroxisomal acyl CoA oxidase ACOX3, the alcohol dehydrogenase ADH6A, dicarbonyl/L-xylulose reductase DCXR, and the aldehyde/retinaldehyde dehydrogenase ALDH1A2, which were upregulated approximately 2 fold (Table 3). ACOX3 was upregulated in both the adipose tissue and the skeletal muscle of the MKR mice. It is one of the three acyl CoA oxidases that perform the first step of fatty acid oxidation in mouse peroxisomes, specifically the oxidation of branched-chain fatty acids [38-45]. ACOX1, which was downregulated in the skeletal muscle and adipose tissue of MKR mice, metabolizes very long-chain fatty acids and long-chain dicarboxylic acids (DCAs). The primary genetic defect in the MKR mice is in the skeletal muscle, and the muscle of the MKR mice have been previously shown to have greater accumulation of fatty acid intermediates compared to WT mice. Therefore, the results of our microarray analysis on the skeletal muscle were consistent with our previously published data showing a decrease in fatty acid oxidation; although the microarray data demonstrated previously unidentified novel changes in gene expression in this pathway in the skeletal muscle of the MKR mice.
No changes in hepatic fatty acid oxidation gene expression were found in the liver of adult diabetic MKR mice
The liver in the MKR mice is known to have increased triglyceride deposits compared to WT mice , and demonstrate hepatic insulin resistance with failure to suppress gluconeogenesis . Therefore, we hypothesized that we would also find a downregulation of fatty acid oxidation genes in the liver of the MKR mice. However, analysis of hepatic gene expression revealed no significant changes in fatty acid oxidation genes in the adult MKR mouse, compared to WT mice (Figure 4, Table 4). Decreased expression of only two genes in this pathway, the acyl CoA synthetase ACSBG1 and serine dehydratase (SDS) was found, while notably, an increase in acetyl CoA acyltransferase 1 (ACAA1) and the acyl CoA synthetase, ACSL6 was found (both of which had decreased expression in the skeletal muscle and adipose tissue of MKR mice). Therefore, despite the insulin resistance in the liver of the MKR mice and known hepatic TG accumulation, no significant decrease in the genes of hepatic FA oxidation found.
Hepatic fatty acid oxidation genes were significantly downregulated in the pre-diabetic MKR mouse
We then examined the gene expression of the liver from 3-week-old pre-diabetic MKR mice to age-matched WT mice. The 3-week old MKR mice have insulin resistance and hyperinsulinemia, but are not hyperglycemic, and are therefore a model of pre-diabetes. Our previous studies have demonstrated that the MKR mice at 3 weeks of age have significant increases in hepatic triglyceride content when compared to WT mice . Using the same methods, our microarray analysis revealed that the FA oxidation pathway was the most significantly dysregulated pathway in the liver of the pre-diabetic MKR mice with significant and substantial decreases in the expression of a large number of FA oxidation genes, compared to the age-matched control mice (Table 5).
Treatment of adult MKR mice with the beta 3-adrenergic receptor agonist CL-316,243 led to improvement in metabolic parameters and altered expression of fatty acid oxidation genes
WT and MKR mice were treated for two weeks with CL-316,243. Consistent with our previous studies, a decrease in random fed blood glucose, fed plasma insulin, serum triglyceride concentration and body fat was observed in the MKR mice after two weeks of treatment (Table 6). Microarray analysis of adipose tissue from MKR mice treated with CL-316,243 revealed a significant change in the expression of a number of the genes in the FA oxidation pathway that were differentially regulated in the MKR compared to WT mice (Figure 5, Table 7). Genes that were downregulated in the MKR adipose tissue, but then were upregulated after CL-316,243 treatment included acyl CoA acyltransferase (ACAA1), the acyl CoA synthetase ACSL6, the alcohol dehydrogenase ADHFE1, and the hydroxysteroid dehydrogenase HSD17B4. Genes that were upregulated in the adipose tissue of MKR mice compared to WT mice and were downregulated after CL-316,343 treatment included the peroxisomal fatty acid oxidation enzyme ACOX3, and the aldehyde dehydrogenases ALDH3A2, ALDH7A1. Chronic CL-316,243 treatment led to a significant downregulation of retinaldehyde dehydrogenases, including ALDH1A1 in the adipose tissue of MKR mice, the deficiency of which has previously been associated with browning of adipose tissue, a phenomenon observed with CL-316,243 treatment . As shown in Table 7, a number of genes in the FA oxidation pathway were altered by CL-316,243 treatment but were not differentially regulated between the WT and MKR mice, and many genes were further downregulated by CL-316,243 treatment. These genes may be downregulated after chronic administration of CL-316,243, but may be increased in the acute setting, or their altered expression may be related to the browning of white adipose tissue observed after CL-316,243 treatment.
The key finding in this study is the significant down-regulation of FA oxidation genes in the skeletal muscle and adipose tissues of MKR mouse model of Type 2 diabetes, and in the pre-diabetic insulin resistant liver of these mice. The finding shows dysregulation of FA oxidation in the tested tissues of the insulin resistant, pre-diabetic and diabetic state and altered expression in genes that regulate peroxisomal beta-oxidation of fatty acids. Treatment of the mice with a beta-3 adrenergic receptor agonist improved the diabetic state and led to the differential regulation of genes involved in peroxisomal fatty acid oxidation and genes that have been previously reported to regulate browning of white adipose tissue in genetic knockout mice, but have not been previously shown to be regulated by a pharmacological agent such as the one used in this study. This robust study on changes in gene expressions after treatment with beta 3-adrenergic receptor agonist, highlights novel findings for better understanding of pathophysiology of T2D in MKR mice, and identifies potential treatment targets.
Computational efforts to understand diseases, such as T2D, with enormous complexity affecting more than one organ, have accelerated the understanding of the perturbed metabolism in such pathophysiological conditions. We have previously utilized computational approaches to develop a multi-confidence level multi-tissue model to study T2D . To understand the mechanism of CL 316,243 compound on diabetic mice, a thorough computational analysis was performed using KEGG pathways analysis based on hypergeometric test and Ingenuity Pathway Analysis (IPA – www.ingenuity.com) based on Fisher’s exact test. The results revealed that treatment with CL-316,243 led to dysregulation of ether lipid metabolism as well as pathways such as fat digestion and absorption, glycerolipid metabolism, and glycerophospholipid metabolism. Additionally, the TCA cycle that was highly overrepresented in MKR fat tissues was back to the same level as the healthy animals following CL-316,243 treatment. In short, the overrepresentation of TCA cycle during T2D was found to be reversed after treatment with the drug.
Defects in fatty acid oxidation may cause the accrual of fatty acyl CoAs and diacylglycerol in the skeletal muscle of obese subjects, resulting in inhibition of insulin-mediated glucose uptake and insulin resistance [48,49]. In the setting of hyperinsulinemia in insulin-treated rodents, hepatic expression of sterol regulatory binding protein 1c (SREBP-1c) increased, along with increased expression of acetyl CoA carboxylase (ACC) and other lipogenic enzymes, suggesting that insulin stimulates lipid synthesis in the liver . Individuals with T2D have increased adipose tissue lipolysis that leads to increased circulating free fatty acids . In a study of 72 human subjects who underwent insulin clamp studies, skeletal muscle and adipose tissue biopsies, impaired fatty acid and branched chain amino acid metabolism were proportional to the degree of insulin resistance . While initial reports stated that white adipose tissue is not an important tissue for fatty acid oxidation , later studies have shown that although the rate of fatty acid oxidation is low in white adipose tissue, normal fatty acid oxidation in adipocytes is important to regulate circulating free fatty acid concentrations and to prevent the development of fatty liver and beta cell dysfunction .
Previous genome wide association studies have identified subsets of gene candidates for T2D . Of the genes identified in that study, some were involved in the metabolism of fatty acids and amino acids, and were also found in our present study to be differentially regulated in the skeletal muscle (ACAA2, ECHS1), and adipose tissue (ECHS1) of the adult MKR mice, and in the liver of pre-diabetic MKR mice (ACAA2, ECHS1). ACAA2 encodes the protein acetyl CoA acyltransferase 2 that catalyzes the final step of mitochondrial fatty acid oxidation along with acetyl coA thiolase (ACAT1), which is also downregulated in the skeletal muscle of adult MKR mice, and in the liver of pre-diabetic mice. ECHS1 encodes enoyl CoA hydratase short chain 1 and catalyzes the second step of mitochondrial FA oxidation.
Beta-oxidation of fatty acids occurs in the peroxisomes and the mitochondria. The mitochondria oxidize the majority of long chain fatty acids (LCFAs) from diet and fat stores, while the peroxisomes oxidize specific carboxylic acids, including very long-chain fatty acids (VLCFAs), branched-chain fatty acids (BCFAs), bile acids and fatty dicarboxylic acids. Mouse peroxisomes have three acyl CoA oxidases that perform the first step of beta-oxidation. ACOX1 the first step oxidizes the straight chain substrates such as very long-chain fatty acids and long-chain dicarboxylic acids (DCAs), ACOX2 is active with the bile acid intermediates and the branched-chain fatty acids and ACOX3 accepts the branched-chain fatty acids [38-45]. The second and third steps of peroxisomal beta-oxidation involve Hsd17b4 (also known as D-peroxisomal bifunctional enzyme) and EHHADH (also known as L-peroxisomal bifunctional enzyme). The peroxisomal beta-oxidation genes that are differentially regulated in the tissues of the MKR mice are the fatty acid transporter SLC27A1 (also known as FATP1 or ACSVL4) in adipose tissue, and SLC27A4 (FATP4, ACSVL5) in adipose tissue, skeletal muscle and pre-diabetic liver. In human studies, links between both SLC27A1 and SLC27A4, and obesity and T2D have been proposed [56,57]. SLC27A1 overexpression in primary human myocytes and HEK293 cells has been reported to trigger increased incorporation of FAs into triacylglycerol, and an increase in 1,2-diacylglycerol acyltransferase activity [58,59]. Reduced triacylglycerol deposition has been observed in 3 T3-L1 adipocytes with shRNA- induced knockdown of SLC27A1 . Insulin has been shown to have varying effects on SLC27A1 expression in adipocytes, with an early repression, followed by a stimulation during differentiation and a subsequent repression in mature adipocytes . Insulin may also stimulate the translocation of SLC27A1 to the cell surface in adipocytes, however SLC27A4 is not recruited to the plasma membrane by insulin [62,63]. Of all the FA oxidation genes found to have lower expression in the pre-diabetic MKR mouse liver, EHHADH expression decreased more than seven fold compared with the WT mice. Adult MKR adipose tissue, skeletal muscle, and the liver of prediabetic MKR mice also showed significant decrease in the expression of HSD17b4 which appears to be able to handle all of the peroxisomal beta oxidation substrates . Recent studies have shown that L-PBE is protective against the over-accumulation of DCAs in mice fed MCFAs in the form of coconut oil . L-PBE knockout mice fed with diet high in MCFAs developed DCA accumulation in the liver, inflammation, hepatic fibrosis and death . HSD17b4 deficient patients and Hsd17b4 knockout mice accumulate VLCFAs, BCFAs and bile acid intermediates [65,66]. Our microarray analysis thus demonstrates significant differences in peroxisomal FA oxidation genes between the MKR and WT mice in different tissues that may indicate that abnormalities in peroxisomal beta-oxidation, contribute to the accumulation of DCAs and play a role in the development of insulin resistance.
The relative lack of difference in fatty acid oxidation gene expression in the liver of adult MKR compared with WT mice was an unexpected finding, as the MKR mice have increased hepatic TG content compared to WT mice . In addition, the FA oxidation pathway was the most significantly down-regulated pathway in the pre-diabetic MKR mice. Possible explanations for the normalization of FA oxidation in the adult MKR mouse are that the described abnormalities in peroxisomal FA oxidation may lead to accumulation of DCAs. DCAs induce PPARα FA oxidation genes , and therefore the accumulation of DCAs in the pre-diabetic liver which may lead to an up-regulation (or normalization) of FA oxidation genes in the adult mice. Alternatively, if FA oxidation is impaired, then gluconeogenesis will not occur and fasting hypoglycemia will result. Therefore, the lack of down-regulated fatty acid oxidation in the adult MKR mice may provide the fuel for gluconeogenesis and the maintenance of hyperglycemia . This hypothesis is supported by the up-regulation of fatty acid transporter genes ACSBG2, ACSL4 as well as acetyl CoA acyltransferase 1 (ACAA1) and the retinaldehyde dehydrogenase 2 (ALDH1A2), which may be involved in hepatic gluconeogenesis .
Treatment with CL-316,243 showed enrichment of many biochemical pathways and an increase in the expression of the peroxisomal FA oxidation enzymes ACAA1 and HSD17b4 in adipose tissue of MKR mice. It is interesting that, chronic CL-316,243 treatment led to a significant down-regulation of several retinaldehyde dehydrogenases, including ALDH1A1 in the white adipose tissue of MKR mice. Aldehyde dehydrogenases convert retinaldehyde to retinoic acid. ALDH1A1 deficiency induces a brown adipose tissue-like transcriptional program in white adipose tissue, with uncoupling of respiration and adaptive thermogenesis . In obese mice, ALDH1A1 knockdown led to a reduction in weight gain and improvement in glucose homeostasis . No previous studies have demonstrated that beta 3 adrenergic receptor agonists decrease the expression of aldehyde dehydrogenases and the function of the Aldh family of enzymes is still being elucidated [69,70]. These changes may contribute to the systemic metabolic improvements that occur with CL-316,243 treatment, and may be targets for future T2D treatments.
Overall, this study uses microarray and computational analysis to demonstrate that FA oxidation genes are significantly altered in the metabolic tissues of diabetic and pre-diabetic MKR mice. In addition, it uncovers biomarker genes such as ACAA1 and HSD17b4 which could be further probed in future studies of pre-diabetes and T2D. Furthermore, it identifies novel changes that occur upon treatment with CL-316,243, which not only explain the lower free fatty acid levels in MKR mice after treatment but also may provide targets for future therapies for T2D.
Whiting DR, Guariguata L, Weil C, Shaw J. IDF diabetes atlas: global estimates of the prevalence of diabetes for 2011 and 2030. Diabetes Res Clin Pract. 2011;94(3):311–21.
Wild S, Roglic G, Green A, Sicree R, King H. Global prevalence of diabetes: estimates for the year 2000 and projections for 2030. Diabetes Care. 2004;(275):1047-53.
Tong Y, Lin Y, Zhang Y, Yang J, Zhang Y, Liu H, et al. Association between TCF7L2 gene polymorphisms and susceptibility to type 2 diabetes mellitus: a large Human Genome Epidemiology (HuGE) review and meta-analysis. BMC Med Genet. 2009;10:15.
Hamman RF, Wing RR, Edelstein SL, Lachin JM, Bray GA, Delahanty L, et al. Effect of weight loss with lifestyle intervention on risk of diabetes. Diabetes Care. 2006;29(9):2102–7.
Diabetes Prevention Program Research G, Knowler WC, Fowler SE, Hamman RF, Christophi CA, Hoffman HJ, et al. 10-year follow-up of diabetes incidence and weight loss in the Diabetes Prevention Program Outcomes Study. Lancet. 2009;374(9702):1677–86.
LeRoith D. Beta-cell dysfunction and insulin resistance in type 2 diabetes: role of metabolic and genetic abnormalities. Am J Med. 2002;113(Suppl 6A):3S–11.
Dream Trial Investigators, Gerstein HC, Yusuf S, Bosch J, Pogue J, Sheridan P, et al. Effect of rosiglitazone on the frequency of diabetes in patients with impaired glucose tolerance or impaired fasting glucose: a randomised controlled trial. Lancet. 2006;368(9541):1096–105.
Weyer C, Tataranni PA, Snitker S, Danforth Jr E, Ravussin E. Increase in insulin action and fat oxidation after treatment with CL 316,243, a highly selective beta3-adrenoceptor agonist in humans. Diabetes. 1998;47(10):1555–61.
Ghorbani M, Himms-Hagen J. Appearance of brown adipocytes in white adipose tissue during CL 316,243-induced reversal of obesity and diabetes in Zucker fa/fa rats. Int J Obes Relat Metab Disord. 1997;21(6):465–75.
McGarry JD. Banting lecture 2001: dysregulation of fatty acid metabolism in the etiology of type 2 diabetes. Diabetes. 2002;51(1):7–18.
Norheim F, Raastad T, Thiede B, Rustan AC, Drevon CA, Haugen F. Proteomic identification of secreted proteins from human skeletal muscle cells and expression in response to strength training. Am J Physiol Endocrinol Metab. 2011;301(5):E1013–21.
Henningsen J, Rigbolt KT, Blagoev B, Pedersen BK, Kratchmarova I. Dynamics of the skeletal muscle secretome during myoblast differentiation. Mol Cell Proteomics. 2010;9(11):2482–96.
Fernandez AM, Kim JK, Yakar S, Dupont J, Hernandez-Sanchez C, Castle AL, et al. Functional inactivation of the IGF-I and insulin receptors in skeletal muscle causes type 2 diabetes. Genes Dev. 2001;15(15):1926–34.
Asghar Z, Yau D, Chan F, Leroith D, Chan CB, Wheeler MB. Insulin resistance causes increased beta-cell mass but defective glucose-stimulated insulin secretion in a murine model of type 2 diabetes. Diabetologia. 2006;49(1):90–9.
Vaitheesvaran B, LeRoith D, Kurland IJ. MKR mice have increased dynamic glucose disposal despite metabolic inflexibility, and hepatic and peripheral insulin insensitivity. Diabetologia. 2010;53(10):2224–32.
Coort SL, van Iersel MP, van Erk M, Kooistra T, Kleemann R, Evelo CT. Bioinformatics for the NuGO proof of principle study: analysis of gene expression in muscle of ApoE3*Leiden mice on a high-fat diet using PathVisio. Genes Nutr. 2008;3(3–4):185–91.
Seda O, Sedova L, Oliyarnyk O, Kazdova L, Krenova D, Corbeil G, et al. Pharmacogenomics of metabolic effects of rosiglitazone. Pharmacogenomics. 2008;9(2):141–55.
Skov V, Glintborg D, Knudsen S, Jensen T, Kruse TA, Tan Q, et al. Reduced expression of nuclear-encoded genes involved in mitochondrial oxidative metabolism in skeletal muscle of insulin-resistant women with polycystic ovary syndrome. Diabetes. 2007;56(9):2349–55.
Luo TH, Zhao Y, Li G, Zhang HL, Li WY, Ldu M. Identification of genes that are differentially expressed in omental fat of normal weight subjects, obese subjects and obese diabetic patients. Zhongguo Ying Yong Sheng Li Xue Za Zhi. 2007;23(2):229–34.
Hulver MW, Berggren JR, Carper MJ, Miyazaki M, Ntambi JM, Hoffman EP, et al. Elevated stearoyl-CoA desaturase-1 expression in skeletal muscle contributes to abnormal fatty acid partitioning in obese humans. Cell Metab. 2005;2(4):251–61.
Fu WJ, Haynes TE, Kohli R, Hu J, Shi W, Spencer TE, et al. Dietary L-arginine supplementation reduces fat mass in Zucker diabetic fatty rats. J Nutr. 2005;135(4):714–21.
Corominola H, Conner LJ, Beavers LS, Gadski RA, Johnson D, Caro JF, et al. Identification of novel genes differentially expressed in omental fat of obese subjects and obese type 2 diabetic patients. Diabetes. 2001;50(12):2822–30.
Kim H, Haluzik M, Asghar Z, Yau D, Joseph JW, Fernandez AM, et al. Peroxisome proliferator-activated receptor-alpha agonist treatment in a transgenic model of type 2 diabetes reverses the lipotoxic state and improves glucose homeostasis. Diabetes. 2003;52(7):1770–8.
Heron-Milhavet L, Haluzik M, Yakar S, Gavrilova O, Pack S, Jou WC, et al. Muscle-specific overexpression of CD36 reverses the insulin resistance and diabetes of MKR mice. Endocrinology. 2004;145(10):4667–76.
Kim H, Haluzik M, Gavrilova O, Yakar S, Portas J, Sun H, et al. Thiazolidinediones improve insulin sensitivity in adipose tissue and reduce the hyperlipidaemia without affecting the hyperglycaemia in a transgenic model of type 2 diabetes. Diabetologia. 2004;47(12):2215–25.
Kim H, Pennisi PA, Gavrilova O, Pack S, Jou W, Setser-Portas J, et al. Effect of adipocyte beta3-adrenergic receptor activation on the type 2 diabetic MKR mice. Am J Physiol Endocrinol Metab. 2006;290(6):E1227–36.
Zhao H, Yakar S, Gavrilova O, Sun H, Zhang Y, Kim H, et al. Phloridzin improves hyperglycemia but not hepatic insulin resistance in a transgenic mouse model of type 2 diabetes. Diabetes. 2004;53(11):2901–9.
Malmgren S, Spegel P, Danielsson AP, Nagorny CL, Andersson LE, Nitert MD, et al. Coordinate changes in histone modifications, mRNA levels, and metabolite profiles in clonal INS-1 832/13 beta-cells accompany functional adaptations to lipotoxicity. J Biol Chem. 2013;288(17):11973–87.
Sparks LM, Moro C, Ukropcova B, Bajpeyi S, Civitarese AE, Hulver MW, et al. Remodeling lipid metabolism and improving insulin responsiveness in human primary myotubes. PLoS One. 2011;6(7):e21068.
Bikopoulos G, da Silva PA, Lee SC, Lakey JR, Der SD, Chan CB, et al. Ex vivo transcriptional profiling of human pancreatic islets following chronic exposure to monounsaturated fatty acids. J Endocrinol. 2008;196(3):455–64.
Swagell CD, Henly DC, Morris CP. Expression analysis of a human hepatic cell line in response to palmitate. Biochem Biophys Res Commun. 2005;328(2):432–41.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.
Fierz Y, Novosyadlyy R, Vijayakumar A, Yakar S, LeRoith D. Insulin-sensitizing therapy attenuates type 2 diabetes-mediated mammary tumor progression. Diabetes. 2010;59(3):686–93.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42(Database issue):D199–205.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Wang K, Hu F, Xu K, Cheng H, Jiang M, Feng R, et al. CASCADE_SCAN: mining signal transduction network from high-throughput data based on steepest descent method. BMC Bioinformatics. 2011;12:164.
Donovan EL, Pettine SM, Hickey MS, Hamilton KL, Miller BF. Lipidomic analysis of human plasma reveals ether-linked lipids that are elevated in morbidly obese humans compared to lean. Diabetol Metab Syndr. 2013;5(1):24.
Poirier Y, Antonenkov VD, Glumoff T, Hiltunen JK. Peroxisomal beta-oxidation–a metabolic pathway with multiple functions. Biochim Biophys Acta. 2006;1763(12):1413–26.
Van Veldhoven PP. Biochemistry and genetics of inherited disorders of peroxisomal fatty acid metabolism. J Lipid Res. 2010;51(10):2863–95.
Wanders RJ, Waterham HR. Biochemistry of mammalian peroxisomes revisited. Annu Rev Biochem. 2006;75:295–332.
Van Veldhoven PP, Vanhove G, Assselberghs S, Eyssen HJ, Mannaerts GP. Substrate specificities of rat liver peroxisomal acyl-CoA oxidases: palmitoyl-CoA oxidase (inducible acyl-CoA oxidase), pristanoyl-CoA oxidase (non-inducible acyl-CoA oxidase), and trihydroxycoprostanoyl-CoA oxidase. J Biol Chem. 1992;267(28):20065–74.
Antonenkov VD, Van Veldhoven PP, Waelkens E, Mannaerts GP. Substrate specificities of 3-oxoacyl-CoA thiolase A and sterol carrier protein 2/3-oxoacyl-CoA thiolase purified from normal rat liver peroxisomes. Sterol carrier protein 2/3-oxoacyl-CoA thiolase is involved in the metabolism of 2-methyl-branched fatty acids and bile acid intermediates. J Biol Chem. 1997;272(41):26023–31.
Antonenkov VD, Van Veldhoven PP, Waelkens E, Mannaerts GP. Comparison of the stability and substrate specificity of purified peroxisomal 3-oxoacyl-CoA thiolases A and B from rat liver. Biochim Biophys Acta. 1999;1437(2):136–41.
Qi C, Zhu Y, Pan J, Usuda N, Maeda N, Yeldandi AV, et al. Absence of spontaneous peroxisome proliferation in enoyl-CoA Hydratase/L-3-hydroxyacyl-CoA dehydrogenase-deficient mouse liver. Further support for the role of fatty acyl CoA oxidase in PPARalpha ligand metabolism. J Biol Chem. 1999;274(22):15775–80.
Dieuaide-Noubhani M, Novikov D, Baumgart E, Vanhooren JC, Fransen M, Goethals M, et al. Further characterization of the peroxisomal 3-hydroxyacyl-CoA dehydrogenases from rat liver. Relationship between the different dehydrogenases and evidence that fatty acids and the C27 bile acids di- and tri-hydroxycoprostanic acids are metabolized by separate multifunctional proteins. Eur J Biochem. 1996;240(3):660–6.
Kiefer FW, Vernochet C, O’Brien P, Spoerl S, Brown JD, Nallamshetty S, et al. Retinaldehyde dehydrogenase 1 regulates a thermogenic program in white adipose tissue. Nat Med. 2012;18(6):918–25.
Kumar A, Harrelson T, Lewis NE, Gallagher EJ, LeRoith D, Shiloach J, et al. Multi-tissue computational modeling analyzes pathophysiology of type 2 diabetes in MKR mice. PLoS One. 2014;9(7):e102319.
Lowell BB, Shulman GI. Mitochondrial dysfunction and type 2 diabetes. Science. 2005;307(5708):384–7.
Simoneau JA, Veerkamp JH, Turcotte LP, Kelley DE. Markers of capacity to utilize fatty acids in human skeletal muscle: relation to insulin resistance and obesity and effects of weight loss. FASEB J. 1999;13(14):2051–60.
Shimomura I, Bashmakov Y, Ikemoto S, Horton JD, Brown MS, Goldstein JL. Insulin selectively increases SREBP-1c mRNA in the livers of rats with streptozotocin-induced diabetes. Proc Natl Acad Sci U S A. 1999;96(24):13656–61.
Coppack SW, Jensen MD, Miles JM. In vivo regulation of lipolysis in humans. J Lipid Res. 1994;35(2):177–93.
Sears DD, Hsiao G, Hsiao A, Yu JG, Courtney CH, Ofrecio JM, et al. Mechanisms of human insulin resistance and thiazolidinedione-mediated insulin sensitization. Proc Natl Acad Sci U S A. 2009;106(44):18745–50.
Coppack SW, Fisher RM, Gibbons GF, Humphreys SM, McDonough MJ, Potts JL, et al. Postprandial substrate deposition in human forearm and adipose tissues in vivo. Clin Sci. 1990;79(4):339–48.
Maassen JA, Romijn JA, Heine RJ. Fatty acid-induced mitochondrial uncoupling in adipocytes as a key protective factor against insulin resistance and beta cell dysfunction: a new concept in the pathogenesis of obesity-associated type 2 diabetes mellitus. Diabetologia. 2007;50(10):2036–41.
Tiffin N, Adie E, Turner F, Brunner HG, van Driel MA, Oti M, et al. Computational disease gene identification: a concert of methods prioritizes type 2 diabetes and obesity candidate genes. Nucleic Acids Res. 2006;34(10):3067–81.
Gertow K, Pietilainen KH, Yki-Jarvinen H, Kaprio J, Rissanen A, Eriksson P, et al. Expression of fatty-acid-handling proteins in human adipose tissue in relation to obesity and insulin resistance. Diabetologia. 2004;47(6):1118–25.
Meirhaeghe A, Martin G, Nemoto M, Deeb S, Cottel D, Auwerx J, et al. Intronic polymorphism in the fatty acid transport protein 1 gene is associated with increased plasma triglyceride levels in a French population. Arterioscler Thromb Vasc Biol. 2000;20(5):1330–4.
Garcia-Martinez C, Marotta M, Moore-Carrasco R, Guitart M, Camps M, Busquets S, et al. Impact on fatty acid metabolism and differential localization of FATP1 and FAT/CD36 proteins delivered in cultured human muscle cells. Am J Physiol Cell Physiol. 2005;288(6):C1264–72.
Hatch GM, Smith AJ, Xu FY, Hall AM, Bernlohr DA. FATP1 channels exogenous FA into 1,2,3-triacyl-sn-glycerol and down-regulates sphingomyelin and cholesterol metabolism in growing 293 cells. J Lipid Res. 2002;43(9):1380–9.
Lobo S, Wiczer BM, Smith AJ, Hall AM, Bernlohr DA. Fatty acid metabolism in adipocytes: functional analysis of fatty acid transport proteins 1 and 4. J Lipid Res. 2007;48(3):609–20.
Man MZ, Hui TY, Schaffer JE, Lodish HF, Bernlohr DA. Regulation of the murine adipocyte fatty acid transporter gene by insulin. Mol Endocrinol. 1996;10(8):1021–8.
Stahl A, Evans JG, Pattel S, Hirsch D, Lodish HF. Insulin causes fatty acid transport protein translocation and enhanced fatty acid uptake in adipocytes. Dev Cell. 2002;2(4):477–88.
Wu Q, Ortegon AM, Tsang B, Doege H, Feingold KR, Stahl A. FATP1 is an insulin-sensitive fatty acid transporter involved in diet-induced obesity. Mol Cell Biol. 2006;26(9):3455–67.
Ding J, Loizides-Mangold U, Rando G, Zoete V, Michielin O, Reddy JK, et al. The Peroxisomal Enzyme L-PBE Is Required to Prevent the Dietary Toxicity of Medium-Chain Fatty Acids. Cell Reports. 2013;5(1):248–58.
Ferdinandusse S, Denis S, Mooyer PA, Dekker C, Duran M, Soorani-Lunsing RJ, et al. Clinical and biochemical spectrum of D-bifunctional protein deficiency. Ann Neurol. 2006;59(1):92–104.
Baes M, Huyghe S, Carmeliet P, Declercq PE, Collen D, Mannaerts GP, et al. Inactivation of the peroxisomal multifunctional protein-2 in mice impedes the degradation of not only 2-methyl-branched fatty acids and bile acid intermediates but also of very long chain fatty acids. J Biol Chem. 2000;275(21):16329–36.
Sun Z, Lazar MA. Dissociating fatty liver and diabetes. Trends Endocrinol Metab. 2013;24(1):4–12.
Kiefer FW, Orasanu G, Nallamshetty S, Brown JD, Wang H, Luger P, et al. Retinaldehyde dehydrogenase 1 coordinates hepatic gluconeogenesis and lipid metabolism. Endocrinology. 2012;153(7):3089–99.
Reichert B, Yasmeen R, Jeyakumar SM, Yang F, Thomou T, Alder H, et al. Concerted action of aldehyde dehydrogenases influences depot-specific fat formation. Mol Endocrinol. 2011;25(5):799–809.
Demozay D, Rocchi S, Mas JC, Grillo S, Pirola L, Chavey C, et al. Fatty aldehyde dehydrogenase: potential role in oxidative stress protection and regulation of its gene expression by insulin. J Biol Chem. 2004;279(8):6261–70.
Funding for the study was provided by the Intramural program of the National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health. This work was also funded in part by the National Cancer Institute grant 1K08CA190770-01 to EJG. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We would like to thank Dr. Derek LeRoith, MD PhD, who generously allowed us to use the MKR mice for these experiments and for his endless advice and support while conducting and analyzing the experiments described in this manuscript, and during the compilation of this manuscript.
The authors declare that they have no competing interests.
All authors were involved in the design of experiments, manuscript preparation, reviewing and editing. In addition, AK and JS performed the RNA isolation, microarray and pathway analysis. EJG performed the animal experiments. All authors read and approved the final manuscript.
Complete list of KEGG enrichment analysis results.