Obesity resistant mechanisms in the Lean polygenic mouse model as indicated by liver transcriptome and expression of selected genes in skeletal muscle
© Simončič et al; licensee BioMed Central Ltd. 2011
Received: 22 August 2010
Accepted: 3 February 2011
Published: 3 February 2011
Divergently selected Lean and Fat mouse lines represent unique models for a polygenic form of resistance and susceptibility to obesity development. Previous research on these lines focused mainly on obesity-susceptible factors in the Fat line. This study aimed to examine the molecular basis of obesity-resistant mechanisms in the Lean line by analyzing various fat depots and organs, the liver transcriptome of selected metabolic pathways, plasma and lipid homeostasis and expression of selected skeletal muscle genes.
Expression profiling using our custom Steroltalk v2 microarray demonstrated that Lean mice exhibit a higher hepatic expression of cholesterol biosynthesis genes compared to the Fat line, although this was not reflected in elevation of total plasma or liver cholesterol. However, FPLC analysis showed that protective HDL cholesterol was elevated in Lean mice. A significant difference between the strains was also found in bile acid metabolism. Lean mice had a higher expression of Cyp8b1, a regulatory enzyme of bile acid synthesis, and the Abcb11 bile acid transporter gene responsible for export of acids to the bile. Additionally, a higher content of blood circulating bile acids was observed in Lean mice. Elevated HDL and upregulation of some bile acids synthesis and transport genes suggests enhanced reverse cholesterol transport in the Lean line - the flux of cholesterol out of the body is higher which is compensated by upregulation of endogenous cholesterol biosynthesis. Increased skeletal muscle Il6 and Dio2 mRNA levels as well as increased activity of muscle succinic acid dehydrogenase (SDH) in the Lean mice demonstrates for the first time that changes in muscle energy metabolism play important role in the Lean line phenotype determination and corroborate our previous findings of increased physical activity and thermogenesis in this line. Finally, differential expression of Abcb11 and Dio2 identifies novel strong positional candidate genes as they map within the quantitative trait loci (QTL) regions detected previously in crosses between the Lean and Fat mice.
We identified novel candidate molecular targets and metabolic changes which can at least in part explain resistance to obesity development in the Lean line. The major difference between the Lean and Fat mice was in increased liver cholesterol biosynthesis gene mRNA expression, bile acid metabolism and changes in selected muscle genes' expression in the Lean line. The liver Abcb11 and muscle Dio2 were identified as novel positional candidate genes to explain part of the phenotypic difference between the Lean and Fat lines.
The 20th century is marked by an explosive increase in obesity and metabolic syndrome commencing in populations of the Western world in the 1960's and is becoming a significant burden also in developing nations . Metabolic syndrome encompasses a cluster of deleterious metabolic events characterized by obesity, hyperglycemia, insulin resistance, hypercholesterolemia, hypertriglyceridemia and hypertension. Environmental determinants influencing the development of the metabolic syndrome include low-cost highly palatable foods and the increasing mechanization of transport and manual tasks following the industrial revolution. While the intake of calories gradually declined over the past decades, physical inactivity and its interaction with genetic predisposition to obesity seem to be an increasingly prevalent determinant of metabolic syndrome development [2, 3].
Previous research in the field of metabolic syndrome and obesity primarily studied monogenic knockout rodent models . Although murine models with single-gene mutations are useful to study the effects of a few genes with major effect on human obesity, monogenic defects account for only a minority of patients with morbid obesity. In contrast, epidemic proportions of the metabolic syndrome and obesity are ascribed to the polygenic effects, where multiple genes interact with many environmental factors over time [5, 6]. Thus, single gene animal models can not explain a large portion of genetic variation that exists in the polygenic form making polygenic obesity mouse models [7, 8], such as the lines used herein, a more suitable genetic resource.
Accordingly, we employed Lean and Fat mouse lines developed from the same base population by long-term (over 60 generations) divergent selection for low or high body fat % . These Fat and Lean mouse lines may represent a unique model of the prevalent polygenic form of human obesity and leanness. Because previous research focused mainly on investigating obesity-susceptibility in the Fat line (e.g., [10–12]), we here aimed to elucidate the molecular basis of obesity-resistance in the Lean line.
We performed a thorough phenotypic characterization of the Lean and Fat lines. In addition, a gene-expression analysis was conducted on the liver tissue using the Steroltalk v2 microarray that was previously used to study hepatic lipid homeostasis [13–15]. This custom microarray is designed for detailed studies of genes from the cholesterol homeostasis and drug metabolism networks, including bile acid metabolism, and several genes from the fatty acid and glucose pathways as well as related transporters for the aforementioned processes.
The results of our study combined with previous findings help to identify, at least in part, a molecular basis for the contrasting phenotype between the Lean and Fat mice and substantiate a hypothesis that differences in metabolic interactions between skeletal muscle and liver tissues play an important part in the phenotype divergence between these lines. The Lean mouse line exhibits characteristics that uncover some of the factors involved in genetic predisposition for obesity-resistance that may also be relevant for obesity and metabolic syndrome basic and applied research in humans.
Animals and maintenance
Inbred lean (Lean) and fat (Fat) mouse lines used here were developed by long-term (over 60 generations) divergent selection for low or high body fat % [9, 16]. Our previous genetic studies demonstrated that these lines differ in large part due to cumulative action of multiple quantitative trait loci (QTLs) of small to medium sized effects, confirming that these lines may present a valuable model of the polygenic form of leanness/obesity [17–19].
Animals were weaned at 3 weeks of age and one pair of males (usually brothers) was housed per cage. Cages were randomly distributed on the rack in order to minimize localized environmental effects and mice were maintained at 21.5°C (± 0.5) with 12 h light/dark cycle. Mice were fed maintenance chow diet (Altromin 1324, Lage, Germany), containing: crude oil, 4%; crude protein, 19%; crude fibre, 6%; ash, 7%; moisture, 13.5%; metabolizable energy, 11.9 MJ/kg where 24 %, 11 % and 65 % of calories from protein, fat and carbohydrates, respectively. Body weights were recorded at 3, 6, 10 and 14 weeks of age. All procedures on animals were performed according to local ethical and regulatory guidelines, which are all in compliance with the EU regulations regarding research on experimental animals.
Phenotype measurements and analyses
At 14 weeks of age (± 3 days) a total of 43 Lean and 42 Fat males were sacrificed after four hour fast by decapitation. Samples of blood, liver, entire left epididymal, left abdominal (perirenal and retroperitoneal), left femoral fat, mesenterial fat, heart and striated muscle from regio cruris (m.gastrocnemius) and thigh muscle (m. quadriceps) were weighed, promptly frozen in liquid nitrogen and stored at -80°C.
Blood was collected in Li-heparin tubes (Sarstedt, Germany) and centrifuged within 3 h of collection at 3000 rpm/min, for 10 min at +4°C. Plasma samples of 10 Lean and 10 Fat mice were individually analyzed for total plasma cholesterol, HDL-cholesterol, plasma triglyceride and glucose concentrations. Plasma total cholesterol and triglyceride concentrations were determined by enzymatic assays adapted to microtiter plates using commercially available reagents (BioMerieux, Lyon, France). Plasma HDL-cholesterol levels were measured after precipitation of apolipoprotein (Apo) B-containing lipoproteins with phosphotungstic acid/Mg (Roche Diagnostics GmbH, Mannheim, Germany). Non-HDL cholesterol was obtained by subtraction of HDL-cholesterol values from total plasma cholesterol.
Plasma lipoproteins were separated by gel filtration chromatography using a Superose 6 HR 10/30 column (Pharmacia, Sweden). The gel was allowed to equilibrate with 10 mM phosphate-buffered saline (PBS) containing 0.01% (wt/vol) EDTA and 0.01% (wt/vol) sodium azide. Plasma was eluted with the buffer at room temperature at a flow rate of 0.2 ml per minute. The effluents were collected in 0.22 ml fractions. Cholesterol and triglyceride concentrations were determined in the eluted fractions as described above.
Cholesterol distribution among lipoproteins was obtained by separation of the major lipoprotein classes (VLDL, IDL+LDL, and HDL) by fast protein liquid chromatography (FPLC) carried out on 200 μl of plasma sample for Fat and Lean line, according to procedure previously detailed . Cholesterol concentrations were determined in the eluted fractions.
Analysis of hepatic lipids in individual samples of 10 Lean and 10 Fat mice was carried out using frozen liver tissue (50 mg) which was homogenized in SET buffer (1 mL ; sucrose 250 mM, EDTA 2 mM and Tris 10 mM), followed by two freeze-thaw cycles and three passages through a 27-gauge syringe needle and a final freeze-thaw cycle to ensure complete cell lysis. Protein content was determined with the BCA method and triglyceride and cholesterol was measured as described above.
Differences between the lines in all measured phenotypic parameters were tested by an unpaired Student t-test on untransformed data as the data exhibited normal distribution.
SDH activity of striated muscles
Oxidative capacity of the muscle fibres (mid-portion of m. gastrocnemius) was determined by the mitochondrial succinate dehydrogenase (SDH) assay as described previously . Briefly, 10 μm thick cryostat sections were incubated for 45 minutes at 37°C in a moisture chamber with the medium composed of equal volumes of 0.2 M phosphate buffer saline (pH 7.6), 0.2 M disodium succinate and aqueous solution of nitro-BT (1 mg/ml). After incubation the slides were washed in tap water, dehydrated and mounted with Synthetic Mountant (Shandon, USA). SDH has an important role in metabolism of tricarbonic acids by the Krebs cycles. In the histochemical reaction the enzyme takes over the hydrogen from disodium succinate and transfers it to the tetrasolium salt (NBT), which turns visible in the form of blue formazan granules.
To analyze the SDH activity in muscle fibres a Nikon Microphot FXA microscope (Nikon Instruments Europe BV, Badhoevendorp, The Netherlands) and the Lucia-G analyzing software (Laboratory Imaging, Prague, Czech Republic) were used. In total 12 of each Lean and Fat line mice were analyzed. 150 muscle fibres of each mouse were measured by densitometry. Higher % densitometry indicated higher transparency resulting from lower SDH-enzyme activity.
10 liver samples were individually homogenized for 10 s using an Ultra-Turrax T8 homogenizer (IKA Labortechnik, Germany) and RNA extracted in »batches« (a pair of Fat and Lean line samples simultaneously) using TRIzol reagent (TRIzol® reagent, Invitrogen, Life technologies). RNA quality was evaluated by ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, USA) and Agilent 2100 bioanalyzer. According to the given RINs (RNA integrity number) a pairwise comparison between Fat and Lean line RNA samples was made in order to couple samples with least quality divergence. RIN values of all samples ranged from 7.5-8.5 suggesting low variability and high overall quality. Maximal tolerated RNA quality divergence for a pair of Lean and Fat-line samples co-hybridized on the same microarray was ± 0.5 RIN.
The Steroltalk v2 cDNA microarray used here was a next generation of the previously described Sterolgene v0 array  with additional cDNA probes detecting a total of 278 mouse genes from selected metabolic pathways relevant to obesity research. The design and preparation of these additional probes was the same as described . The Steroltalk v2 cDNA microarray was focused on targeted functional groups such as enzymes from cholesterol biosynthesis (18/21 probes; coverage 85.7%), nuclear receptor families (45/49 probes; coverage 93.75%; http://www.nursa.org/) and members of the cytochrome P450 (38/103 probes, coverage 36.9%) . Additionally, key enzymes of glucose and fatty acid metabolism; transporters of xenobiotics, cholesterol and bile acids; and proteins involved in cholesterol plasma transport, circadian rhythm, and selected cell signalling pathways involved in the regulation of cholesterol homeostasis were included. A full list of genes contained on the custom Steroltalk v2 array along with associated functional groups is provided in Additional file 1.
Individual total RNA samples from 10 Lean and 10 Fat-line animals were labelled using Cy3 and Cy5 dye. To 20 μg of Lean and Fat line total RNA sample, spike of control RNAs was added: 250 pg of Firefly Luciferase mRNA (Promega, Madison, WI, USA) and 0.5 μL of either test or reference spike mix from Lucidea Universal Scorecard kit (Amersham Biosciences, GE Healthcare UK limited, Little Chalfont, UK). mRNA was reverse transcribed to amino-allyl cDNA using 2.5 μg of Oligo dT (Invitrogen, Carlsbad, CA, USA), 400 U of SuperScript™ III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) and 1 μL of 10 mM amino-allyl dUTP (Sigma, St Louis, MI, USA) according to the manufacturer's protocol. Reactions were stopped after 2 hours by addition of 10 μL of 0.5 M EDTA and of 1 M NaOH and were further incubated at 65°C for 15 min. 10 μL of 1 M HCl was added and cDNA was purified using MinElute PCR Purification Kit (Qiagen GmBH, Hilden, D) according to the manufacturer's protocol with exception of using a phosphate buffer (5 mM KPO4 pH 8.5 in 80% ethanol) for the washing step and MilliQ water for elution. Purified amino-allyl cDNA was dried and resuspended in 4.5 μL of 0.2 M Na2CO3 (pH 9.0) and 4.5 μL of cyanin-3 or cyanin-5 dye in DMSO (Amersham Biosciences, GE Healthcare UK limited, Little Chalfont, UK). The labelling reaction was incubated at room temperature for two hours, precipitated by addition of 35 μL of 0.1 M Na acetate (pH 5.2) and purified reaction using a MinElute PCR Purification Kit (Qiagen GmBH, Hilden, D) were added. Labelled cDNA was eluted in water and 1 μL was used for evaluation on ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, Delaware, USA). 10 Lean and Fat line cDNA samples, respectively, were co-hybridized in pairs and dye swaps were performed. Using LifterSlip cover glasses (Erie Scientific Company, Portsmouth, NH, USA) samples were hybridized for 16 h at 65°C using buffers 3×SSC and 0.2% SDS in humidified hybridization chambers (HybChambers, GeneMachines, San Carlos, CA, USA). The slides were washed and scanned as described .
Images were analyzed using Array-Pro Analyzer 4.5 (Media Cybernetics, Bethesda, MD, USA). Three arrays did not pass quality control parameters (mainly too high signal to noise ratios) and hence 17 arrays were included in downstream statistical analyses. Data were normalized using OWNormalize, a widget for explorative normalization of focused microarrays, implemented within Orange software . Data were first filtered to exclude spots of low quality and then normalized using LOWESS fit to "spike in" control RNAs (Firefly luciferase and Lucidea Universal Scorecard) according to their average intensity. The raw and normalized gene expression data of 17 arrays together with experimental information are deposited in Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/, accession number GSE24967) in compliance with MIAME standards . Differential expression was assessed using Student's t-test at probability of type I error α=0.05. Expression of genes was considered on individual basis and therefore no multiple testing correction has been applied.
Quantitative reverse transcription polymerase chain reaction (qRT-PCR)
Primer sequences used in quantitative qRT-PCR analyses
3-hydroxy-3-methylglutaryl-Coenzyme A reductase
Low density lipoprotein receptor
Insulin induced gene 1
Peroxisome proliferator activated receptor α
Deiodinase iodothyronine type II
18 S ribosomal RNA
Differential gene expression in skeletal muscle between the Lean and Fat mice
Deiodinase iodothyronine type II
Lean mice exhibit a liver transcriptome profile suggesting improved cholesterol homeostasis, bile acid, glucose and lipoprotein metabolism
Differential expression between the Lean and Fat mice in liver genes
NAD(P) dependent steroid dehydrogenase-like
Sterol-C5-desaturase (fungal ERG3, delta-5-desaturase) homolog (S. cerevisae)
Mevalonate (diphospho) decarboxylase
Farnesyl diphosphate synthetase
Farnesyl diphosphate farnesyl transferase 1
Acyl-CoA synthetase short-chain family member 2
Bile acid and glucose metabolism
Pyruvate dehydrogenase (lipoamide) beta
Solute carrier organic anion transporter family, member 1a4
Solute carrier organic anion transporter family, member 1a1
Solute carrier organic anion transporter family, member 1b2
Solute carrier family 10 (sodium/bile acid cotransporter family), member 1
ATP-binding cassette, sub-family G (WHITE), member 5
ATP-binding cassette, sub-family B (MDR/TAP), member 11
Cytochrome P450, family 8, subfamily b, polypeptide 1
Cytochrome P450, family 2, subfamily a, polypeptide 5
Cytochrome P450, family 2, subfamily b, polypeptide 10
Cytochrome P450, family 2, subfamily b, polypeptide 9
Cytochrome P450, family 2, subfamily b, polypeptide 9
Cytochrome P450, family 3, subfamily a, polypeptide 41A
Cytochrome P450, family 2, subfamily c, polypeptide 69
Cytochrome P450, family 1, subfamily a, polypeptide 2
Complement component 9
Serum amyloid P-component
Complement component 4B (Childo blood group)
C-reactive protein, pentraxin-related
Nuclear receptor subfamily 5, group A, member 2 (LRH-1)
Acetyl-Coenzyme A acetyltransferase 2
Lysosomal acid lipase A
Low density lipoprotein receptor
Insulin induced gene 1
Nitric oxide synthase 1, neuronal
Adiponectin receptor 2
Transporter genes involved in the uptake of plasma bile acids into hepatocytes such as Slco1a1, Slco1a4, Slco1b2, Slc10a1, and the cholesterol efflux transporter Abcg5 were down-regulated in Lean mice. However, Abcb11, which controls bile acid excretion, was up-regulated. Cyp8b1, a gene encoding a bile acid synthesis enzyme, was also expressed at a higher level in Lean mice (Table 3).
Gluconeogenesis genes G6Pc ( glucose-6-phosphatase, catalytic) and Pdhb (pyruvate dehydrogenase beta) were markedly down-regulated in Lean mice supporting previously identified predisposition of Lean line to normoglycemia and resistance to diabetes development (Table 3).
The Ldlr, Lip1 and Adipor2 receptors were up-regulated in livers of Lean mice suggesting that livers of Lean mice have an increased ability to capture atherogenic LDL-cholesterol as well as to respond to the white adipose tissue derived hormone adiponectin. In contrast, the transcription factor Ppara, involved in triglyceride homeostasis, was down-regulated in Lean mice (Table 3).
Many of the cytochromes P450 functioning in drug and steroid hormone metabolism were down-regulated (Cyp2a5, Cyp2b10, Cyp2b9 and Cyp3a41) although two (Cyp2c69 and Cyp1a2) were up-regulated (Table 3). However, it is important to note that the cDNA probes used on the microarrays do not enable efficient differentiation between the aforementioned CYP members, which share a relatively high sequence homology. For example, probes for Cyp2b10 and Cyp2b9 cross-react with other Cyp2b family members, and the same holds true for the Cyp2c, Cyp3a and Cyp2a probes. Therefore, to determine exactly which gene from the CYP2,3 families is differentially expressed in our lines, future experiment based on allele-specific hybridization probes or gene specific RT-PCR assays should be conducted.
Lean mice display an improved profile of adiposity, plasma lipoproteins, plasma bile acids and liver triglycerides parameters
Weights of fat depots and heart normalized to body weight in Lean and Fat mice
No. of individuals
BODY WEIGHT at 14 weeks
28.17 ± 1.15
41.76 ± 3.22*
L = 23; F = 24
ratio (mg/g body wt)
1.86 ± 069
9.97 ± 1.08a
L = 43; F = 42
ratio (mg/g body wt)
4.43 ± 1.2
17.67 ± 2.32 a
L = 43; F = 42
ratio (mg/g body wt)
3.34 ± 0.85
23.24 ± 3.12 a
L = 43; F = 42
ratio (mg/g body wt)
6.43 ± 1.26
22.33 ± 3.87 a
L = 18; F = 15
ratio (mg/g body wt)
6.13 ± 0.87
4.58 ± 0.69 a
L = 43; F = 42
Lean mice show increased skeletal muscle oxidative capacity
The main objective of the present study was to identify candidate genes and molecular mechanisms responsible for obesity-resistance in the Lean line. Additionally, we aimed to investigate metabolic parameters of skeletal muscle and liver tissue and how they relate to blood lipoproteins, triglycerides, bile acids and liver transcriptome. Phenotypic characterization and search of differentially expressed genes was also carried out to identify novel positional candidate genes that map to the regions of quantitative trait loci determined in previous crosses between the Fat and Lean mice.
Non-atherogenic lipoprotein profile in Lean mice
Plasma and liver cholesterol concentrations are, under balanced diet conditions, determined by the affluence from hepatic de-novo cholesterol biosynthesis. Additionally, excretion rate of cholesterol and bile acids also contribute to inter-individual variation of total plasma and liver cholesterol.
In the current experiment a large number of liver cholesterol biosynthesis genes were upregulated in the Lean line, but this was not reflected in an elevation of total liver and plasma cholesterol. This inconsistency could be explained in two ways: 1) either upregulation of cholesterol biosynthesis genes was not sufficient to increase the metabolic flow in cholesterol biosynthesis pathway that would result in an increased level of de novo hepatic cholesterol; 2) or that novel and excess hepatic cholesterol was excreted from the liver or metabolized to bile acids. Since we have not measured the metabolic rate of hepatic cholesterol biosynthesis or determined protein levels, we can not confirm yet that the increase in mRNA expression resulted in increased metabolic rate of cholesterol biosynthesis. Regulation of the cholesterol synthesis genes is under control of the SREBP signalling pathway  and upregulation observed in the Lean line is probably a result of activation of the SREBP transcription factor cascade. This is supported by our finding that Insig-1 and Ldlr, both being SREBP target genes, were up-regulated in the Lean mice. The hypothesis that cholesterol in the Lean mice would be excreted from the liver or metabolized to bile acids at an enhanced rate has limited supportive results. Higher efflux of cholesterol from the liver to bile is not likely since one of the genes involved in this process, Abcg5, is downregulated in the Lean mice. Increase in bile acid synthesis is plausible, since Cyp8b1 was up-regulated. However; Cyp7a1, a rate-limiting enzyme of bile acid synthesis, was not changed (data from RT-PCR not shown). Nevertheless, Cyp7a1 gene exhibits large interindividual differences in the level of expression and this hinders statistical significance. Also, we found that the nuclear receptor subfamily 5, group A, member 2 (Nr5a2, also known as LRH-1), a regulator of bile acid synthesis genes, is upregulated in Lean mice. In order to fully explain why up-regulation of cholesterol pathway genes in the Lean line did not lead to increased levels of total liver and plasma cholesterol, future detailed studies should be focused on protein and metabolic level of cholesterol management in Lean and Fat lines.
Although we did not find statistically significant differences in total liver and plasma cholesterol between the lines, the Lean mice livers did contain significantly less triglycerides per mg of liver tissue than the Fat mice - Lean mice had only about 50% of Fat mice triglyceride content (Figure 4). Abcb11, a gene overexpressed in the Lean line liver, is known to enhance the rate of triglyceride excretion from the liver [26, 27]. Also, in a transgenic Abcb11 overexpression model , increased liver Abcb11 expression showed ameliorating effects on hepatic steatosis by increasing excretion of biliary phospholipids and cholesterol. It is possible therefore that upregulation of Abcb11 in the Lean mice might account for their diminished content of liver triglycerides. Abcb11 is also a strong positional candidate because it maps within a QTL Fob1 on chromosome 2, detected in our previous cross between Lean and Fat mice . Also supportive is a study in humans detecting 30 mutations in the ABCB11 gene that were associated with different clinical phenotypes and triglyceride measurements . Therefore, Abcb11 could be one of the potential causal allelic variants between the Lean and Fat mice, which can aid in explaining a part of the phenotype divergence between the strains. Elevated HDL and upregulation of some bile acids synthesis and transport genes suggests enhanced reverse cholesterol transport in the Lean line. As the cholesterol excretion rate is higher, Lean line compensates this by upregulating of endogenous cholesterol biosynthesis as demonstrated by the array experiment.
We established that Lean mice have significantly elevated plasma HDL. Interestingly, results from a previous study on humans imply that more exercise leads to increased total HDL and the average size of HDL particles . FPLC analysis of plasma revealed similar observations - Lean mice have higher HDL plasma concentrations and also tend to have HDL particles of increased size. Interpretation for improved HDL profile in Lean mice can be further substantiated by genetic studies on HDL metabolism, where several QTLs were demonstrated to determine plasma HDL concentration, some also mapping to QTL regions detected in Fat and Lean mice [17, 31]. It is therefore possible that some of the obesity QTLs identified in our previous genetic studies have pleiotropic effects affecting both body fat % and HDL levels in our polygenic animal model [17–19]. The explanation for the observed higher concentration of HDL and lower concentration of LDL-cholesterol in Lean mice may therefore lie in allelic differences in some of the HDL QTLs between the Fat and Lean mice.
Plasma bile acids and changes in skeletal muscle metabolism
The bile acid pool in plasma is maintained by the enterohepatic recirculation. The uptake of bile acids that return to the liver after intestinal absorption is mainly mediated by a family of SLCO transporters (Slco1a1, Slco1a4 , Slco10a1)  whose expression can be modulated by plasma bile acids and cytokines. In our Lean mice down-regulation of genes encoding the aforementioned hepatic bile acid transporters along with significantly elevated concentrations of plasma bile acids was detected. Other studies have also demonstrated that bile acids can down-regulate expression of Slco genes and up-regulate Abcb11, which is in agreement with our results obtained in the Lean line. Additional support comes from a study identifying significant dose-dependent mRNA suppression of Slco1a1 (Oatp1) and Slco1a4 (Oatp2) in mice treated with cytokine IL6 . Detailed in-vitro and in-vivo studies defined IL6 as by far the most pronounced transcriptional regulator of, for example, Slco1a1, Slco1a4 and Abcb11[33, 34]. IL6, also termed as »exercise factor«, is highly expressed in contracting (exercising) skeletal muscle and is the most highly secreted systemic muscle cytokine during exercise . The expression level of Il6 mRNA serves as a predictor of plasma IL6 protein . In line with these reports, our analysis demonstrates high levels of expression of skeletal muscle Il6 mRNA in Lean mice. This increased Il6 expression in the Lean line could be a consequence of increased physical activity of Lean mice reported in our previous experiments . Thus, relatively high IL6 plasma levels in Lean mice might suppress the expression of the Slco family of genes subsequently leading to the increase in plasma bile acids. However, future studies are needed to confirm that the muscle-derived cytokine IL6 in the Lean mice is directly responsible for down-regulation of Slco genes and increase in plasma bile acids.
Because Lean mice contain elevated levels of plasma bile acids we sought to interpret other potential metabolic actions involving bile acids. Elevated levels of plasma bile acids were found to lead to endocrine actions with consequent impact on whole body energy homeostasis . This effect is elicited through the increased mRNA expression of type II iodothyronine deiodinase (Dio2) which was exemplified in the mouse brown adipose tissue and in human skeletal muscle myoblasts . Dio2 expression was elevated in our Lean mice and also maps close to the Fob2 QTL, previously identified as one of four loci significantly contributing to the divergent body fat % in our polygenic mouse model . In humans, polymorphisms at the Dio2 locus determines susceptibility to hypertension and diabetes [39, 40]. Given the expression difference in the present study and its location to previously mapped QTL in the same mice, it is plausible that the Dio2 allelic variant also exists in this animal model, and may consequently determine differential expression and downstream consequences on phenotype in the Lean line.
The overall results of our current study are in line with our previous energy budget results which suggested that Lean mice expend more energy on physical activity  and our follow-up study demonstrating higher running wheel activity as well as higher exercise-independent posture allocation to more energy demanding positions (Non-Exercise Activity Thermogenesis; NEAT--standing, fidgeting, etc.) in Lean mice . Increased oxidative capacity of muscle determined in the current study and resulting increased levels of Il6 expression may therefore indeed stem from higher physical and NEAT activities leading to downstream favourable effects on the liver transcriptome, cholesterol, bile acid, glucose and lipoprotein homeostasis.
However, the identified differences in liver transcriptome and other parameters and their associations with muscle metabolism and physical activity can explain only a part of a large phenotypic divergence between the Lean and Fat line. In a polygenic model it is expected that many alleles at different loci will contribute cumulatively to obesity resistance in various components of energy intake and energy expenditure. Since we showed in this study that Lean mice have similar food intake to the Fat line, we can exclude hyperphagia or genetic differences in energy intake as a major component for the phenotypic divergence between the lines. However, we have already shown earlier that thermoregulation  may also be increased in the Lean line and potential differences in resting metabolic rate can not be excluded. Therefore, more comprehensive experiments monitoring behavioural traits in combination with metabolic changes to measure differences in fuel oxidation will be needed to provide a better basis for future quantitative evaluations of energy expenditure components and fuel utilisation. Our current study reveals several candidate genes as well as important target tissues on which to base new targeted hypotheses to be tested in the future.
Polygenic mouse models of leanness/obesity used in our study comprise complex obesity-resistant (Lean line) and obesity-susceptible (Fat line) gene networks. This study was designed to help identify candidate genes and reveal some metabolic mechanisms responsible for the phenotypic difference between these lines. Results of liver transcriptome analysis along with markers of oxidative capacity (SDH enzyme activity, Il6 expression) in skeletal muscle suggests of multiple tissue metabolic interactions in the obesity resistant Lean line. Candidate genes, particularly those mapping to previously detected QTL regions (Abcb11 and Dio2) between Lean and Fat mice, and new phenotypic parameters obtained in this study will give us leads to identify novel genes underlying this complex network. The results of the present study combined with previous findings help to elucidate, at least in part, a molecular basis for the contrasting phenotype between the Lean and Fat mice and substantiate a hypothesis concerning metabolic interactions between the physical activity, skeletal muscle and liver tissue. We demonstrate that the Lean mouse line is a potentially good model for identifying novel genes and mechanisms in obesity resistance.
We thank the Slovenian Research Agency for co-financing this work (core funding program P4-0220, project Syntol and Young Investigator grant). We are grateful to the Centre for functional genomics and biochips, Ljubljana, Slovenia http://cfgbc.mf.uni-lj.si/ for collaboration in the part of microarray analysis. We would also like to thank Zala Prevoršek, Nataša Debeljak, Nataša Toplak, Gorazd Tompa, Andrej Razpet, Ana Zanjkovič and Irena Oven for technical support and Emmanuelle Vallez for expert technical assistance and Jim McWhir for proofreading the manuscript. Finally, we are grateful to anonymous reviewers for meticulous review and constructive criticism that helped us to improve this manuscript.
- Batsis JA, Nieto-Martinez RE, Lopez-Jimenez F: Metabolic syndrome: From global epidemiology to individualized medicine. Clin Pharmacol Ther. 2007, 82 (5): 509-524. 10.1038/sj.clpt.6100355.PubMedView ArticleGoogle Scholar
- Weinsier RL, Hunter GR, Heini AF, Goran MI, Sell SM: The etiology of obesity: Relative contribution of metabolic factors, diet, and physical activity. Am J Med. 1998, 105 (2): 145-150. 10.1016/S0002-9343(98)00190-9.PubMedView ArticleGoogle Scholar
- Neel JV: The "thrifty genotype" in 1998. Nutr Rev. 1999, 57 (5): S2-S9. 10.1111/j.1753-4887.1999.tb01782.x.PubMedView ArticleGoogle Scholar
- Butler AA, Cone RD: The melanocortin receptors: Lessons from knockout models. Neuropeptides. 2002, 36 (2-3): 77-84. 10.1054/npep.2002.0890.PubMedView ArticleGoogle Scholar
- Romao I, Roth J: Genetic and environmental interactions in obesity and type 2 diabetes. J Am Diet Assoc. 2008, 108 (4): S24-S28. 10.1016/j.jada.2008.01.022.PubMedView ArticleGoogle Scholar
- Ichihara S, Yamada Y: Genetic factors for human obesity. Cell Mol Life Sci. 2008, 65 (7-8): 1086-1098. 10.1007/s00018-007-7453-8.PubMedView ArticleGoogle Scholar
- Brockmann GA, Bevova MR: Using mouse models to dissect the genetics of obesity. Trends Genet. 2002, 18 (7): 367-376. 10.1016/S0168-9525(02)02703-8.PubMedView ArticleGoogle Scholar
- Speakman J, Hambly C, Mitchell S, Krol E: The contribution of animal models to the study of obesity. Laboratory animals. 2008, 42 (4): 413-432. 10.1258/la.2007.006067.PubMedView ArticleGoogle Scholar
- Bunger L, Hill WG: Inbred lines of mice derived from long-term divergent selection on fat content and body weight. Mamm Genome. 1999, 10 (6): 645-648. 10.1007/s003359901063.PubMedView ArticleGoogle Scholar
- Morton NM, Densmore V, Wamil M, Ramage L, Nichol K, Bunger L, Seckl JR, Kenyon CJ: A polygenic model of the metabolic syndrome with reduced circulating and intra-adipose glucocorticoid action. Diabetes. 2005, 54 (12): 3371-3378. 10.2337/diabetes.54.12.3371.PubMedView ArticleGoogle Scholar
- Bunger L, Forsting J, McDonald KL, Horvat S, Duncan J, Hochscheid S, Baile CA, Hill WG, Speakman JR: Long-term divergent selection on body fatness in mice indicates a regulation system that is independent of leptin production and reception. Faseb J. 2002, 16 (13): 85-Google Scholar
- Simoncic M, Rezen T, Juvan P, Fievet C, Staels B, Rozman D, Horvat S: Transcriptome analysis revealed association of some P450 genes with obesity in a polygenic obese mouse model. Acta Chim Slov. 2008, 55 (1): 101-110.Google Scholar
- Tacer KF, Kuzman D, Seliskar M, Pompon D, Rozman D: TNF-alpha interferes with lipid homeostasis and activates acute and proatherogenic processes. Physiol Genomics. 2007, 31 (2): 216-227. 10.1152/physiolgenomics.00264.2006.View ArticleGoogle Scholar
- Rezen T, Juvan P, Fon Tacer K, Kuzman D, Roth A, Pompon D, Aggerbeck LP, Meyer UA, Rozman D: The Sterolgene v0 cDNA microarray: a systemic approach to studies of cholesterol homeostasis and drug metabolism. BMC Genomics. 2008, 9: 76-10.1186/1471-2164-9-76.PubMedPubMed CentralView ArticleGoogle Scholar
- Conklin LS O-HM: Nutritional considerations in pediatric inflammatory bowel disease. Expert Review of Gastroenterology & Hepatology. 2010, 4 (3): 305-317.View ArticleGoogle Scholar
- Sharp GL, Hill WG, Robertson A: Effects of Selection on Growth, Body-Composition and Food-Intake in Mice .1. Responses in Selected Traits. Genet Res. 1984, 43 (1): 75-92. 10.1017/S0016672300025738.PubMedView ArticleGoogle Scholar
- Horvat S, Bunger L, Falconer VM, Mackay P, Law A, Bulfield G, Keightley PD: Mapping of obesity QTLs in a cross between mouse lines divergently selected on fat content. Mamm Genome. 2000, 11 (1): 2-7. 10.1007/s003350010002.PubMedView ArticleGoogle Scholar
- Stylianou IM, Christians JK, Keightley PD, Bunger L, Clinton M, Bulfield G, Horvat S: Genetic complexity of an obesity QTL (Fob3) revealed by detailed genetic mapping. Mamm Genome. 2004, 15 (6): 472-481. 10.1007/s00335-004-3039-z.PubMedView ArticleGoogle Scholar
- Prevorsek Z, Gorjanc G, Paigen B, Horvat S: Congenic and bioinformatics analyses resolved a major-effect Fob3b QTL on mouse Chr 15 into two closely linked loci. Mamm Genome. 2010, 21 (3-4): 172-185. 10.1007/s00335-010-9252-z.PubMedView ArticleGoogle Scholar
- Nachlas MM TK, De Souza E, Cheng CS, Seligman AM: Cytochemical demonstration of succinic dehydrogenase by the use of a new p-nitrophenyl substituted ditetrazole. J Histochem Cytochem. 1957, 5: 420-436.PubMedView ArticleGoogle Scholar
- Curk T, Demsar J, Xu QK, Leban G, Petrovic U, Bratko I, Shaulsky G, Zupan B: Microarray data mining with visual programming. Bioinformatics. 2005, 21 (3): 396-398. 10.1093/bioinformatics/bth474.PubMedView ArticleGoogle Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, et al: Minimum information about a microarray experiment (MIAME) - toward standards for microarray data. Nature Genet. 2001, 29 (4): 365-371. 10.1038/ng1201-365.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- Bunger L, Forsting J, McDonald KL, Horvat S, Duncan J, Hochscheid S, Baile CA, Hill WG, Speakman JR: Long-term divergent selection on fatness in mice indicates a regulation system independent of leptin production and reception. Faseb J. 2003, 17 (1): 85-87.PubMedGoogle Scholar
- Engelking LJ, Liang GS, Hammer RE, Takaishi K, Kuriyama H, Evers BM, Li WP, Horton JD, Goldstein JL, Brown MS: Schoenheimer effect explained - feedback regulation of cholesterol synthesis in mice mediated by Insig proteins. J Clin Invest. 2005, 115 (9): 2489-2498. 10.1172/JCI25614.PubMedPubMed CentralView ArticleGoogle Scholar
- Stieger B, Meier Y, Meier PJ: The bile salt export pump. Pflugers Arch. 2007, 453 (5): 611-620. 10.1007/s00424-006-0152-8.PubMedView ArticleGoogle Scholar
- Geier A, Wagner M, Dietrich CG, Trauner M: Principles of hepatic organic anion transporter regulation during cholestasis, inflammation and liver regeneration. Biochim Biophys Acta-Mol Cell Res. 2007, 1773 (3): 283-308. 10.1016/j.bbamcr.2006.04.014.View ArticleGoogle Scholar
- Figge A, Lammert F, Paigen B, Henkel A, Matern S, Korstanje R, Shneider BL, Chen F, Stoltenberg E, Spatz K, et al: Hepatic overexpression of murine Abcb11 increases hepatobiliary lipid secretion and reduces hepatic steatosis. J Biol Chem. 2004, 279 (4): 2790-2799. 10.1074/jbc.M307363200.PubMedView ArticleGoogle Scholar
- Pauli-Magnus C, Stieger B, Meier Y, Kullak-Ublick GA, Meier PJ: Enterohepatic transport of bile salts and genetics of cholestasis. J Hepatol. 2005, 43 (2): 342-357. 10.1016/j.jhep.2005.03.017.PubMedView ArticleGoogle Scholar
- Kraus WE, Houmard JA, Duscha BD, Knetzger KJ, Wharton MB, McCartney JS, Bales CW, Henes S, Samsa GP, Otvos JD, et al: Effects of the amount and intensity of exercise on plasma lipoproteins. N Engl J Med. 2002, 347 (19): 1483-1492. 10.1056/NEJMoa020194.PubMedView ArticleGoogle Scholar
- Wang XS, Paigen B: Genetics of variation in HDL cholesterol in humans and mice. CircRes. 2005, 96 (1): 27-42.Google Scholar
- Fickert P, Zollner G, Fuchsbichler A, Stumptner C, Pojer C, Zenz R, Lammert F, Stieger B, Meier PJ, Zatloukal K, et al: Effects of ursodeoxycholic and cholic acid feeding on hepatocellular transporter expression in mouse liver. Gastroenterology. 2001, 121 (1): 170-183. 10.1053/gast.2001.25542.PubMedView ArticleGoogle Scholar
- Hartmann G, Cheung AKY, Piquette-Miller M: Inflammatory cytokines, but not bile acids, regulate expression of murine hepatic anion transporters in endotoxemia. J Pharmacol Exp Ther. 2002, 303 (1): 273-281. 10.1124/jpet.102.039404.PubMedView ArticleGoogle Scholar
- Siewert E, Dietrich CG, Lammert F, Heinrich PC, Matern S, Gartung C, Geier A: Interleukin-6 regulates hepatic transporters during acute-phase response. Biochem Biophys Res Commun. 2004, 322 (1): 232-238. 10.1016/j.bbrc.2004.07.102.PubMedView ArticleGoogle Scholar
- Pedersen BK, Akerstrom TCA: Role of myokines in exercise and metabolism. J Appl Physiol. 2007, 103 (3): 1093-1098. 10.1152/japplphysiol.00080.2007.PubMedView ArticleGoogle Scholar
- Febbraio MA, Pedersen BK: Muscle-derived interleukin-6: mechanisms for activation and possible biological roles. Faseb J. 2002, 16 (11): 1335-1347. 10.1096/fj.01-0876rev.PubMedView ArticleGoogle Scholar
- Simoncic M, Horvat S, Stevenson PL, Bunger L, Holmes MC, Kenyon CJ, Speakman JR, Morton NM: Divergent physical activity and novel alternative responses to high fat feeding in polygenic fat and lean mice. Behav Genet. 2008, 38 (3): 292-300. 10.1007/s10519-008-9199-y.PubMedView ArticleGoogle Scholar
- Watanabe M, Houten SM, Mataki C, Christoffolete MA, Kim BW, Sato H, Messaddeq N, Harney JW, Ezaki O, Kodama T, et al: Bile acids induce energy expenditure by promoting intracellular thyroid hormone activation. Nature. 2006, 439 (7075): 484-489. 10.1038/nature04330.PubMedView ArticleGoogle Scholar
- Mentuccia D, Proietti-Pannunzi L, Tanner K, Bacci V, Pollin TI, Poehlman ET, Shuldiner AR, Celi FS: Association between a novel variant of the human type 2 deiodinase gene Thr92Ala and insulin resistance - Evidence of interaction with the Trp64Arg variant of the beta-3-adrenergic receptor. Diabetes. 2002, 51 (3): 880-883. 10.2337/diabetes.51.3.880.PubMedView ArticleGoogle Scholar
- Gumieniak O, Perlstein TS, Williams JS, Hopkins PN, Brown NJ, Raby BA, Williams GH: Ala92 type 2 deiodinase allele increases risk for the development of hypertension. Hypertension. 2007, 49 (3): 461-466. 10.1161/01.HYP.0000256295.72185.fd.PubMedView ArticleGoogle Scholar
- Yang J, Goldstein JL, Hammer RE, Moon YA, Brown MS, Horton JD: Decreased lipid synthesis in livers of mice with disrupted Site-1 protease gene. Proc Natl Acad Sci USA. 2001, 98 (24): 13607-13612. 10.1073/pnas.201524598.PubMedPubMed CentralView ArticleGoogle Scholar
- Yabe D, Komuro R, Liang GS, Goldstein JL, Brown MS: Liver-specific mRNA for Insig-2 down-regulated by insulin: Implications for fatty acid synthesis. Proc Natl Acad Sci USA. 2003, 100 (6): 3155-3160. 10.1073/pnas.0130116100.PubMedPubMed CentralView ArticleGoogle Scholar
- Patel DD, Knight BL, Wiggins D, Humphreys SM, Gibbons GF: Disturbances in the normal regulation of SREBP-sensitive genes in PPAR alpha-deficient mice. J Lipid Res. 2001, 42 (3): 328-337.PubMedGoogle Scholar
- Gutala RV, Reddy PH: The use of real-time PCR analysis in a gene expression study of Alzheimer's disease post-mortem brains. J Neurosci Methods. 2004, 132 (1): 101-107. 10.1016/j.jneumeth.2003.09.005.PubMedView ArticleGoogle Scholar