Skip to main content
  • Research article
  • Open access
  • Published:

Functional and genetic analysis in type 2 diabetes of Liver X receptor alleles – a cohort study

Abstract

Background

Liver X receptor alpha (LXRA) and beta (LXRB) regulate glucose and lipid homeostasis in model systems but their importance in human physiology is poorly understood. This project aimed to determine whether common genetic variations in LXRA and LXRB associate with type 2 diabetes (T2D) and quantitative measures of glucose homeostasis, and, if so, reveal the underlying mechanisms.

Methods

Eight common single nucleotide polymorphisms in LXRA and LXRB were analyzed for association with T2D in one French cohort (N = 988 cases and 941 controls), and for association with quantitative measures reflecting glucose homeostasis in two non-diabetic population-based samples comprising N = 697 and N = 1344 adults. Investigated quantitative phenotypes included fasting plasma glucose, serum insulin, and HOMAIR as measure of overall insulin resistance. An oral glucose tolerance test was performed in N = 1344 of adults. The two alleles of the proximal LXRB promoter, differing only at the SNP rs17373080, were cloned into reporter vectors and transiently transfected, whereupon allele-specific luciferase activity was measured. rs17373080 overlapped, according to in silico analysis, with a binding site for Nuclear factor 1 (NF1). Promoter alleles were tested for interaction with NF1 using direct DNA binding and transactivation assays.

Results

Genotypes at two LXRB promoter SNPs, rs35463555 and rs17373080, associated nominally with T2D (P values 0.047 and 0.026). No LXRA or LXRB SNP associated with quantitative measures reflecting glucose homeostasis. The rs17373080 C allele displayed higher basal transcription activity (P value < 0.05). The DNA-mobility shift assay indicated that oligonucleotides corresponding to either rs17373080 allele bound NF1 transcription factors in whole cell extracts to the same extent. Different NF1 family members showed different capacity to transactivate the LXRB gene promoter, but there was no difference between promoter alleles in NF1 induced transactivation activity.

Conclusion

Variations in the LXRB gene promoter may be part of the aetiology of T2D. However, the association between LXRB rs35463555 and rs17373080, and T2D are preliminary and needs to be investigated in additional larger cohorts. Common genetic variation in LXRA is unlikely to affect the risk of developing T2D or quantitative phenotypes related to glucose homeostasis.

Peer Review reports

Background

Type 2 diabetes (T2D) is a common cause of morbidity and mortality and is caused by a western sedentary life style with high calory diet in combination with susceptibility genes. Without life style changes, available pharmaceutical treatments of T2D are only partially efficient. Recently, large scale allelic association analyses have begun to identify the genetic factors underlying susceptibility to T2D [1]. This new knowledge may permit unravelling of pathways leading to development of T2D and thereby aid in the development of more efficient prevention and therapy of this disease.

Liver X receptor alpha (LXRA, NR1H3) and beta (LXRB, NR1H2) are established sensors of lipid and cholesterol homeostasis [2]. Recently, a large body of literature has indicated a role of LXRs also in glucose metabolism and homeostasis. Several studies have demonstrated that LXR agonists reduce plasma glucose concentrations and increase insulin sensitivity in different models of diabetes and insulin resistance [3, 4]. LXRs seem to improve glucose metabolism at different levels. In experimental models, LXR activation inhibits hepatic gluconeogenesis and glucose output [5, 6]. Other results support that LXRs regulate peripheral glucose uptake. Activation of LXRs promotes glucose uptake and oxidation in muscle, and expression of the insulin responsive glucose transporter Glucose transporter type 4 (GLUT4) in adipocytes [5, 7, 8]. LXRs are also implicated in regulation of insulin secretion. Activation of LXRs increases glucose dependent insulin secretion in vitro from pancreatic beta-cell line cultures [9]. This effect may be mediated through LXRB since Lxrb -/- mice have lower basal insulin levels and are glucose intolerant due to impaired glucose induced insulin secretion [10]. Together, these observations point towards a physiological role of LXRs in glucose homeostasis.

We have previously reported that two LXRB SNPs associate with obesity in a Swedish cohort [rs35463555 (previously known as LB44732G>A), rs2695121, and marginally, rs17373080 (P value 0.06)] [11]. Furthermore, one LXRA SNP (rs11039155) was reported to associate with the metabolic syndrome in two French cohorts [12]. LXRA is encoded on human chromosome 11p11.2 and LXRB on chromosome 19q13.3.

Here we analyze if eight common SNPs in LXRA and LXRB associate with T2D in one French case-control cohort, and to phenotypes related to glucose homeostasis in two non-diabetic population-based cohorts comprising adult French and Swedish subjects. We also investigate the impact of one LXRB promoter SNP on gene transcription in vitro.

Methods

Subjects

The present study includes one large French cohort that was used for analysis of association between LXR gene SNPs and T2D. In addition, a subset of the non-diabetics in the French cohort, that is D.E.S.I.R., as well as the non-diabetics in the Swedish Stockholm Diabetes Prevention Program (SDPP) cohort, were analyzed for association between LXR SNPs and quantitative phenotypes related to glucose homeostasis. The cohorts used for quantitative trait analyses were population based. The study was approved by the ethical committee of the Karolinska University Hospital, the French National Commission for Informatics and Liberty and Hotel Dieu Ethics Committee. Informed written consent was received from all participants.

French cohort

The French T2D group included 988 unrelated T2D patients (N = 662 males and 326 females, age at examination 61.0 ± 10.2 years, BMI 26.3 ± 2.5 kg/m2) recruited at the Endocrinology-Diabetology Department of the Corbeil-Essonnes Hospital [13]. The type 2 diabetes status was defined according to the World Health Organization 1999 criteria [14]. The T2D patients included in the study were selected on a BMI < 30 kg/m2 and all were unrelated. 346 T2D subjects were treated with insulin, 624 were treated with oral anti-diabetic drugs, and 18 subjects had other anti-diabetic treatment. The 941 control subjects were selected on fasting plasma glucose < 6.1 mmol/l and BMI < 27 kg/m2 (N = 380 males and 561 females, age at examination 62.2.0 ± 7.5 years, BMI 23.2 ± 2.0 kg/m2) and all were unrelated. Of the 941 controls, 697 were selected from the D.E.S.I.R. general French population [15] and 244 were selected from the French families recruited at the CNRS unit in Lille. All subjects are of French Caucasian origin. We have compared the genotypic distributions within the control group depending on the original selections. No significant difference was observed. For 1021 subjects (case and controls) of this study we could estimate the proportion of European ancestry from a previous genome-wide genotyping [16]. We performed logistic regression taking into account also the gender of the individuals, and observed at least 90% (in 1005 subjects) or 99% (in 835 subjects) of European ancestry (results below not included in the text). Thus, we can potentially exclude a possible stratification in the French study sample. Moreover, these results are supported by several other association studies [1624]. Information about family history of T2D was available for the French family cohort and for most of the subjects that were included as controls there was no T2D in the ancestry. Family history of T2D was not available from the D.E.S.I.R. participants. The size of the French cohort provided 70% (90%) power to detect an allele with a frequency of 20% among controls and odds ratio (OR) 1.25 (1.35) to develop T2D assuming a threshold P value of 0.05 and dominant mode of inheritance (QUANTO, http://hydra.usc.edu/gxe/) [25].

Swedish cohort

For the Swedish SDPP study, a short questionnaire was sent to all men born between 1938 and 1957 living in four municipalities in the outskirts of Stockholm and all women born between 1942 and 1961 living in the same municipalities and one additional municipality, asking about country of birth and presence of diabetes in subjects and relatives. Persons were excluded as a result of diabetes and foreign origin. In a second step, subjects with family history of diabetes, together with subjects randomly selected from those without family history of diabetes, matched to the first group by age and municipality were invited to a health examination. Subjects with DNA available were included in this study.

SDPP included 1024 normoglycemic controls with normal glucose tolerance (N = 778 males and 246 females, age at examination 46.7 ± 4.9 years, BMI 25.5 ± 4.4 kg/m2) and 320 subjects with impaired glucose tolerance (IGT) (N = 174 males and 146 females, age at examination 48.2 ± 4.5 years, BMI 29.2 ± 5.3 kg/m2) [2629].

The size of the population-based non-diabetic Swedish SDPP and French D.E.S.I.R. cohorts together provided 53% (98%) power to detect an allele with a frequency of 15% and marginal βG 0.2 (0.4) for HOMAIR as measure of insulin resistance, assuming a threshold P value of 0.05 and dominant mode of inheritance (QUANTO, http://hydra.usc.edu/gxe/) [25].

Clinical evaluation

In the morning a fasting venous blood sample was obtained for determination of fasting levels of insulin and glucose by the routine chemistry laboratory of the hospitals, and for extraction of DNA using standard protocols. An oral glucose tolerance test was performed in the SDPP subjects. Diabetes was defined according to the World Health Organization 1999 criteria: fasting plasma glucose ≥ 7.0 mmol/l and/or 2-h plasma glucose ≥ 11.1 mmol/l, or ongoing treatment with oral antidiabetic agents and/or insulin [14]. IGT was defined as fasting plasma glucose < 7.0 mmol/l and 2-h plasma glucose 7.8–11.0 mmol/l. Insulin resistance index HOMAIR (homeostasis model assessment) was calculated as fasting serum insulin (μU/ml) × fasting plasma glucose (mmol/l)/22.5 [30]. The diagnosis of T2D as opposed to type 1 diabetes was based on clinical records and judgment by the investigating physicians.

SNP selection and genotyping

We genotyped altogether five SNPs in the LXRA gene (Figure 1A). Three of the genotyped LXRA SNPs cover the common variation in this gene in a French cohort according to Legry et al [12]. According to Legry et al a, LXRA comprises three common haplotypes, which are tagged by SNPs genotyped in the present study. We were unable to design a multiplex assay including the SNP rs11039155, which in French cohorts has been associated with the metabolic syndrome. However, we genotyped instead rs2279238 that is in complete linkage disequilibrium (LD) with rs11039155 [12]. Since the SNPs genotyped by Legry did not fully cover the two untranslated 5' exons of LXRA we added two more SNPs in this region from our previous study of LXRA in obesity, that is rs4752822 and rs61896015 (the latter previously labelled LA9462C>A) [11]. For LXRB, we genotyped three SNPs that covered 57% of alleles as r2 ≥ 0.80 according to our own sequencing and genotyping of the gene in a Swedish cohort (Figure 1B) [11]. The reason for the relatively low coverage was that we were unable to design a high throughput genotyping assay for two variations; a repeat and a SNP located in this repeat. The three SNPs cover the three common LXRB gene haplotypes according to our own previous genotyping of this gene in a Swedish cohort [11].

Figure 1
figure 1

Schematic figure of LXRA and LXRB genes with position of genotyped SNPs according to NCBI as indicated.

Both samples were genotyped using matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry [31]. Primers are provided on request. Due to initial low call rate, rs2248949 was re-genotyped by Taqman in a subset of the Swedish cohort (C_1993248_20, Applied Biosystems, Foster City, CA, USA). The final genotype call rate was ≥ 93% for all SNPs. The accuracy was 99.99% according to duplicate analysis of, on average, 2% of the total genotypes. Hardy-Weinberg equilibrium (HWE) calculations were performed in Haploview to ensure that each marker was within population equilibrium [32].

Cloning of the LXRB promoter into reporter vectors and in vitromutagenesis

The -244/+1163 sequence of the human LXRB gene was cloned into the pGL3 basic luciferase reporter vector (Promega) using the Kpn I and Mlu I sites applying primers 5'-ATCAGGTACCGGCCGCAGGCTCAGAGAAGCGCATGAATGAGCTAA-3' and 5'-ATCACTCGAGGGTGGGGTCACGGAGCAGCCTGTAGAATACAGGGGATTGAGAG-3' with the restriction enzyme sites underlined. The G allele (bold) of rs17373080 was generated using the following primer (5'-TAAAGCCAGAAAGCGCGGGGCTGGAGGTTT-3' using the QuikChange® XL Site-Directed Mutagenesis Kit (Stratagene). DNA sequencing confirmed the identity of all clones.

Transient transfections

The mouse pancreatic beta-cell line MIN6 was maintained in Dulbecco's modified Eagle's medium (DMEM, 4.5 g/L glucose, 41965-039), the human hepatoma cell line HuH7 was grown in DMEM with 0.11 G/L NaPyr (41966-029), and African green monkey kidney cells Cos7 in DMEM with 1 g/L glucose (31885-023). All media were from GIBCO-BRL. The cell lines were supplemented with fetal bovine serum (Cos7, HuH7: 10%, MIN6: 15%), 2 mM L-glutamine (only MIN6), and penicillin/streptomycin at a final concentration of 100 U/ml and 100 μg/ml, respectively. Additionally, 50 μM β-mercaptoethanol was added to the MIN6 cells. Cells were grown under 5% CO2 at 37°C. Approximately 4 × 104 cells per well were plated in 24 well plates and transiently transfected using Lipofectamine 2000 (Invitrogen, Carlsbad, CA) according to the manufacturer's protocol. Each well received 50 ng of reporter vector and 125 ng of expression vector, respectively. Vehicle vector was added to ensure equal amounts of DNA in each transfection. Cells were transfected for 24 h and thereafter lysed in 25 mM TAE, 1 mM EDTA, 10% glycerol, 1% Triton X-100, and 2 mM DTT. Luciferase activities were measured using a Luciferase Assay Kit (BioThema, Umeå, Sweden) in a luminometer (Luminoscan Ascent, Thermo electron Corporation, Waltham, MA). The NF1 transcription factors were cloned into the pCH-expression vectors and were a kind gift of Dr. Gronostajski, Cleveland Clinic Foundation Research Institute.

Whole cell extracts

MIN6 cells were grown to confluence in 100 mm dishes, washed with PBS and incubated in TEN buffer (40 mM Tris-HCl, 1 mM EDTA, 150 mM NaCl) for 4 min. Cells were then removed with a cell scraper and pelleted by centrifugation at 3,500 rpm for 2 min at 4°C. Cell pellets were freeze dried on dry ice and resuspended in 50 μl ice-cold buffer C (10 mM Hepes-KOH pH 7.9, 0.4 M NaCl, 0.1 mM EDTA, 5% glycerol, 1 mM DTT, 0.5 mM PMSF). After another round of freeze drying, cell debris were removed by centrifugation for 5 min at 13,000 rpm at 4°C. The supernatant corresponds to MIN6 whole cell extracts.

Electro mobility shift assay (EMSA)

The oligo including a consensus site for nuclear factor 1 (NF1) is 5'-GATCTTATTTTGGATTGAAGCCAATATGAG-3'. The rs17373080 C allele oligo is 5'-GCTAAAGCCACAAAGCGCGGGG-3' with the C allele in bold. The rs17373080 G allele oligo is 5'-GCTAAAGCCAGAAAGCGCGGGG-3' with G allele in bold. 5.0 μg of the respective forward and reverse oligos were annealed in 20 mM Tris-HCl pH 7.8; 2 mM MgCl2; 50 mM NaCl by heating to 95°C for 5 min and slow cooling by 1.5°C/min for 47 cycles. Oligonucleotide probes were labeled by mixing 0.2 μg annealed oligo with 250 μM non-radioactive dATP, dGTP, dTTP respectively, 1× Klenow buffer, 20 μCi 32P-labeled dCTP (Amersham Pharmacia) and 1 Unit Klenow polymerase. Samples were incubated for 20 min at room temperature. Reactions were stopped by adding 0.5 M EDTA. Probes were purified using G-25 Nick Columns (Amersham Pharmacia) and the efficiency of labeling determined using the 1214 Rackbeta liquid scintillation counter (LKB Wallac). For binding reactions, 5 μl of MIN6 whole cell extracts were incubated with 4 × 104 cpm of radiolabeled oligonucleotide in binding buffer (4 μg BSA, 2 μg poly dI-dC, 12 mM HEPES (pH 7.9), 12% glycerol, 0.12 mM EDTA, 0.9 mM MgCl2, 0.6 mM dithiothreitol, 0.6 mM phenylmethylsulfonyl fluoride). Binding reactions were incubated for 20 min at RT, electrophoresed at 150 V for 3 h in 6% polyacrylamide gels at 4°C, dried and finally analyzed by autoradiography. In competition assays, the unlabeled competitor was added prior to the nuclear extracts.

Statistical analyses

The Finetti software and Pearson's Chi2 (d.f. = 1) were used to compare allele and genotype frequencies between cases and controls and to calculate odds ratios (OR) [33]. The genetic models analyzed by the Finetti software were defined as homozygous/recessive: (Case_11*Control_22)/(Case_22*Control_11) and allele positivity/dominant: ((Case_12+Case_11)*Control_22)/(Case_22*(Control_12+Control_11)). To take into account the impact of gender on susceptibility to T2D, we used logistic regression with genotype and sex as independent variables. Haploview was used to calculate LD between SNPs (D' and r2), estimate haplotypes, and to test for association between haplotypes and T2D [32]. The program Thesias was used to test two-by-two sliding window SNPs haplotypes [34]. This program determines both the haplotype frequency and its effect on the phenotype through an estimation-maximisation algorithm. In this analysis, we also considered gender as a covariate.

ANCOVA with SNP genotypes and gender as main effects and age, and when appropriate BMI, as simple regressor was used to analyze differences in HOMAIR, fasting plasma glucose, serum insulin, and BMI between LXR genotypes in the non-diabetic population-based cohorts, SDPP and normoglycemic non-obese samples from the D.E.S.I.R. respectively. The analyses were repeated in joint analyses of both cohorts including cohorts as an additional main effect. 2-h plasma glucose was available and analyzed as above in the Swedish cohort only. Quantitative phenotypes were normalized by natural logarithm transformation before analysis to become normally distributed. In the transient transfections functional studies the student's t-test was applied.

Results

Analysis of association between LXRgene SNPs and T2D

We first analyzed LXRs SNPs for association with T2D in the French case-control cohort. Altogether we genotyped five LXRA SNPs and three LXRB SNPs. All SNPs had a call rate of ≥ 93% and were in HWE in the French cohort (Table 1). We analyzed LXR SNPs for association with T2D under different genetic models. No SNP allele associated to T2D (Table 1). However, genotype at the LXRB promoter SNPs rs35463555 and rs17373080 associated nominally to T2D assuming a dominant T2D-predisposing common allele, P value 0.047 and 0.026 respectively (Table 2). This was due to an increased risk of T2D among heterozygous and protection against T2D among subjects homozygous for the rare allele (Table 1). None of analyzed SNPs showed significant association with T2D when taking into account the gender effect. SNPs within LXRA and LXRB were in strong LD and together built one haploblock for the respective genes (Figure 2). No LXRA and no LXRB haplotype associated with T2D (Table 3). Furthermore, no two-by-two sliding window SNPs haplotypes associated with T2D (results not shown).

Table 1 Analysis of association between LXR alleles and T2D in French cases vs control
Table 2 Association of LXRB promoter SNP genotypes with T2D in the French cohort
Table 3 Analysis of LXR haplotypes in T2D in the French cohort
Figure 2
figure 2

Linkage disequilibrium between LXRA and LXRB SNPs as indicated. Shown are r2 values. D' was 1.00 between all SNPs within each gene.

Analysis of association between LXRSNPs and quantitative phenotypes reflecting glucose homeostasis

We next analyzed LXR gene SNPs for impact on quantitative phenotypes related to glucose homeostasis in non-diabetic subjects. For this purpose, we analyzed the normoglycemic non-obese population-based samples selected from the French D.E.S.I.R. cohort and the Swedish SDPP cohort. Phenotypes of these cohorts are shown in Table 4. In the Swedish cohort an oral glucose tolerance test was performed and the analyzed cohort included 320 subjects with IGT. Genotyping result for the D.E.S.I.R. samples is included in the controls in Table 1. Genotype results for SDPP are shown in Table 5. All SNPs had call rates ≥ 96%. All SNPs were in HWE. Allele frequencies for LXRB SNPs were almost identical in the Swedish and French cohorts whereas allele frequencies for LXRA SNPs differed 4–10% between the populations. LXR SNPs had no impact on fasting plasma glucose or serum insulin, insulin resistance measured as HOMAIR or BMI in either the D.E.S.I.R. or the SDPP cohort, nor in the joint analysis of both cohorts (results not shown). Nor did LXR SNPs associate to the response to an OGTT in SDPP (results not shown).

Table 4 Phenotypic distribution of cohorts
Table 5 LXR genotypes in the Swedish population based SDPP cohort

Functional analysis of the LXRBrs17373080 SNP

The C allele of rs17373080 associates with higher transcription activity in vitro

We have previously determined the transcription start site of the human LXRB gene promoter [35]. The LXRB SNP rs17373080 is located in the proximal promoter close to the transcription start site and could potentially have an impact on the transcription potential of the LXRB gene promoter. Hence we performed functional studies investigating this SNP. The -244/+1163 gene region including the rs17373080 with either the C or the G allele was PCR cloned in front of the luciferase gene reporter (Figure 3A). Otherwise the promoters were identical. The activities of the two different alleles were analyzed in transient transfection assays in three different cell lines. There was a significantly higher basal activity for the C allele in all three cell lines investigated (Figure 3B) indicating that the C allele is able to maintain a higher basal transcription of the LXRB gene.

Figure 3
figure 3

A. Schematic view of the reporter gene construct driven by the endogenous -244/+1163 LXRB gene region. The black dot indicates rs17373080. B. Basal activities of the -244/+1163 rs17373080 G allele and C allele reporter constructs. Using in vitro mutagenesis with the G allele construct as a template a new C allele construct was generated. The results presented are mean differences of more than 20 independent transfection experiments. The activities of the G allele in each cell line analyzed were set to 100% ± SEM. * P < 0.05 by Student's t test.

A functional NF1 binding site overlaps rs17373080

A theoretical search for transcription factor binding sites using the transcription element search system (TESS) suggested a binding site for the NF1 transcription family members overlapping the rs17373080 position (Figure 4A) [36]. The different alleles were tested for interaction with NF1 using direct DNA binding and transactivation assays. The DNA-mobility shift assay indicated equal binding of components in whole cell extracts to oligonucleotides corresponding to the C and G alleles, respectively (Figure 4B, lanes 2 and 8) and the DNA-protein interaction was competed away by increasing amounts of an unlabelled NF1 consensus binding site (lanes 3–5 and 9–11). NF1A1, NF1B and NF1X produced by in vitro transcription/translation gave protein-DNA complexes of the same size as that produced by whole cell extracts (results not shown). Different NF1 family members (NF1A1, NF1B, NF1C and NF1X) showed different transactivation capacity on the LXRB gene promoter (Figure 5). However, no difference in induced transactivation activity of the NF1s was observed between the G (Figure 5A) and C (Figure 5B) promoter alleles.

Figure 4
figure 4

A. A putative NF1 binding site in the promoter region of the LXRB gene as suggested by TESS (underlined) [36]. The rs17373080 SNP is indicated in bold. B. EMSA with MIN6 whole cell extracts (WCE) and labelled oligonucleotides with the rs17373080 C and G variants, respectively. Lane 1: Free oligo C, 2: Oligo C + MIN6 WCE, 3–5: Oligo C + MIN6 WCE + increasing amounts of unlabeled NF1 consensus oligonucleotide, 6: Empty 7: Free oligo G, 8: Oligo G + MIN6 WCE, 9–11: Oligo G + MIN6 WCE + increasing amounts of unlabeled NF1 consensus oligonucleotide. C. EMSA with rs17373080 SNP oligo C and MIN6 whole cell extract or in vitro translated (IVT) NF1. Lane 1: Free oligo C, 2: Oligo C + MIN6 whole cell extract, 3: Oligo C + MIN6 whole cell extract + 250× unlabeled NF1 consensus oligonucleotide, 4: Oligo C + Empty IVT expression vector (vehicle) and 5–7: Oligo C + IVT NF1A1, NF1B and NF1X, respectively. The arrows indicate the positions of the NF1 interactions.

Figure 5
figure 5

Transient co-transfections with the -244/+1163 C (A) or G (B) reporter constructs and members of the NF1 family of transcription factors in three different cell lines as indicated. Each transient transfection was made in triplicates and the results presented are mean differences of three independent experiments. C or G alleles with co-transfection of empty pCH-vector were set as controls. Both C and G allele controls are set to 100% and the effect of NF1 expression compared to either the C-Control (A) or G-Control (B). * indicated significant difference (P < 0.05) ± SEM.

Discussion

We have analyzed common genetic variation in LXRA and LXRB for association with T2D in one French cohort and for association with quantitative measures of glucose homeostasis in two non-diabetic population-based samples from France and Sweden. Genotype at the LXRB promoter SNPs rs17373080 and rs35463555 nominally associates to T2D. rs17373080, which regulates transcription in vitro, may directly affect the risk of developing T2D whereas the association between T2D and rs35463555, which is located in a repeat, could be due to strong LD with 17373080. No LXRA SNPs associate to the investigated phenotypes.

We have previously reported that LXRB SNPs associate with obesity in a Swedish cohort [rs35463555 (previously known as LB44732G>A), rs2695121, and marginally, rs17373080 (P value 0.06)] [11]. The genotype at rs35463555 and rs17373080 predisposing to obesity in the previous study was found to also associate to T2D in this study. The association between rs35463555, rs17373080 and T2D reported here would become non-significant following Bonferroni-correction for analysis of several LXR SNPs and phenotypes. However, since our results are confirmatory in nature, that is the rs35463555 and, marginally, rs17373080 SNPs have previously been associated with another metabolic phenotype, namely obesity, we think it is too stringent to perform Bonferroni correction. Similarly, since analysed phenotypes are dependent, Bonferroni correction may be too stringent. On the other hand, it is important to note that observed association P values for rs35463555 and rs17373080 (0.047 and 0.026) are weak and would become non-significant following adjustment for just two independent analyses. We did not observe any association between LXRB SNPs and BMI in this study. This may be due to phenotypic differences between cohorts. The population-based cohorts D.E.S.I.R. and SDPP that we analyzed here contain markedly less obese people than the previous Swedish study. Together, previous and present results support that LXRB promoter genotype regulates susceptibility to adiposity as well as T2D. However, replications of these associations in additional large cohorts are necessary before establishing an association between LXRB genotype, obesity and T2D. Published genome wide association studies of T2D have not reported association with LXR genes [1624, 3739]. However, it is uncertain to what extent these studies cover the rather small LXRA and LXRB genes since the full datasets are unavailable. In the French cohort, gender distribution was the opposite in T2D cases versus controls. The main concern with different gender distribution in cases and controls is reduced power to identify T2D genes due to gender impact on T2D susceptibility and potential gene-gender interaction.

An LXRB promoter SNP could contribute to development of T2D by altering mRNA levels of LXRB in one or more of the organs regulating glucose homeostasis; LXRB is expressed in pancreas, liver, muscle, as well as adipose tissue [40]. It has recently been reported that LXRB is the main LXR paralogue in pancreatic beta-cells, that activated LXR induces insulin secretion from beta-cells and that mice lacking the Lxrb gene have less glucose stimulated insulin secretion [9]. The T2D-associated LXRB promoter SNP rs35463555 is located in a large repeat region of ~2500 base pairs and is therefore difficult to clone and study functionally. An in silico search using TESS revealed no binding sites for human transcription regulators covering rs35463555 [36]. Therefore, we did not perform any functional studies of this SNP. As for the proximal LXRB promoter SNP rs17373080, according to our results the GG genotype protected against T2D whereas the G allele showed reduced reporter gene activity, which suggests that the G allele is associated with lower LXRB mRNA levels. Thus, these results are contrary to what one would expect based on the results in the Lxrb knock-out mice. However, rs17373080 may regulate overall glucose homeostasis through other mechanisms and in other organs in humans.

We observed that the rs17373080 SNP overlapped an NF1 site. The NF1 site was functional in recruiting NF1. However, we did not observe any differences between the G and C alleles in DNA-binding and transactivation of a reporter gene by NF1 family members. It cannot be excluded that rs17373080 overlaps additional transcription factor binding sites and that the binding and/or activity of these sites is affected by the various rs17373080 alleles. It might also be that the differences between the C and G alleles with regard to transcription factor-induced transcription activity are too subtle to be detected with existing methods. However, our results indicate that the transcriptional activity of NF1 could be tissue or cell specific as similar results were observed in Cos7 and MIN6 cells, but not in HuH7 cells. Speculatively, these results could suggest that the relevance of the SNP rs17373080 in developing T2D could be tissue specific and further studies should be focusing on the role of LXRB in these cell types where NF1 was shown to affect the transcription rate of the LXRB gene.

There was no association between SNPs covering the common variation in LXRA and susceptibility to T2D or levels of plasma glucose, serum insulin, and insulin resistance measured as HOMAIR. According to power calculations our cohorts were large enough to detect also a modest impact of LXRA on these phenotypes and it is therefore unlikely that this gene directly regulates glucose homeostasis. We did not have access to blood lipids in the SDPP cohort. This prohibited us from determining which subjects fulfilled the criteria for the metabolic syndrome and from analysing LXR genotypes for influence on this phenotype. We hypothesize, based on previously reported association of LXRA alleles with the metabolic syndrome and HDL levels, together with the lack of association with T2D and measures of glucose homeostasis reported here, that the impact of LXRA on susceptibility to the metabolic syndrome is mediated via an impact on lipid turnover rather than glucose homeostasis [12].

Conclusion

The common LXRB promoter SNPs rs35463555 and rs17373080 may regulate susceptibility to T2D, whereas common genetic variation in LXRA is unlikely to affect the risk of developing T2D or quantitative phenotypes related to glucose homeostasis. The association between LXRB rs35463555 and rs17373080 and T2D needs to be investigated in additional large cohorts.

References

  1. Frayling TM, McCarthy MI: Genetic studies of diabetes following the advent of the genome-wide association study: where do we go from here?. Diabetologia. 2007, 50 (11): 2229-2233. 10.1007/s00125-007-0825-7.

    Article  CAS  PubMed  Google Scholar 

  2. Zelcer N, Tontonoz P: Liver X receptors as integrators of metabolic and inflammatory signaling. J Clin Invest. 2006, 116 (3): 607-614. 10.1172/JCI27883.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Cao G, Liang Y, Broderick CL, Oldham BA, Beyer TP, Schmidt RJ, Zhang Y, Stayrook KR, Suen C, Otto KA, et al: Antidiabetic action of a liver x receptor agonist mediated by inhibition of hepatic gluconeogenesis. J Biol Chem. 2003, 278 (2): 1131-1136. 10.1074/jbc.M210208200.

    Article  CAS  PubMed  Google Scholar 

  4. Grefhorst A, van Dijk TH, Hammer A, Sluijs van der FH, Havinga R, Havekes LM, Romijn JA, Groot PH, Reijngoud DJ, Kuipers F: Differential effects of pharmacological liver X receptor activation on hepatic and peripheral insulin sensitivity in lean and ob/ob mice. Am J Physiol Endocrinol Metab. 2005, 289 (5): E829-838. 10.1152/ajpendo.00165.2005.

    Article  CAS  PubMed  Google Scholar 

  5. Laffitte BA, Chao LC, Li J, Walczak R, Hummasti S, Joseph SB, Castrillo A, Wilpitz DC, Mangelsdorf DJ, Collins JL, et al: Activation of liver X receptor improves glucose tolerance through coordinate regulation of glucose metabolism in liver and adipose tissue. Proc Natl Acad Sci USA. 2003, 100 (9): 5419-5424. 10.1073/pnas.0830671100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Cha JY, Repa JJ: The liver X receptor (LXR) and hepatic lipogenesis. The carbohydrate-response element-binding protein is a target gene of LXR. J Biol Chem. 2007, 282 (1): 743-751. 10.1074/jbc.M605023200.

    Article  CAS  PubMed  Google Scholar 

  7. Kase ET, Wensaas AJ, Aas V, Hojlund K, Levin K, Thoresen GH, Beck-Nielsen H, Rustan AC, Gaster M: Skeletal muscle lipid accumulation in type 2 diabetes may involve the liver X receptor pathway. Diabetes. 2005, 54 (4): 1108-1115. 10.2337/diabetes.54.4.1108.

    Article  CAS  PubMed  Google Scholar 

  8. Dalen KT, Ulven SM, Bamberg K, Gustafsson JA, Nebb HI: Expression of the insulin-responsive glucose transporter GLUT4 in adipocytes is dependent on liver X receptor alpha. J Biol Chem. 2003, 278 (48): 48283-48291. 10.1074/jbc.M302287200.

    Article  CAS  PubMed  Google Scholar 

  9. Efanov AM, Sewing S, Bokvist K, Gromada J: Liver X receptor activation stimulates insulin secretion via modulation of glucose and lipid metabolism in pancreatic beta-cells. Diabetes. 2004, 53 (Suppl 3): S75-78. 10.2337/diabetes.53.suppl_3.S75.

    Article  CAS  PubMed  Google Scholar 

  10. Gerin I, Dolinsky VW, Shackman JG, Kennedy RT, Chiang SH, Burant CF, Steffensen KR, Gustafsson JA, MacDougald OA: LXRbeta is required for adipocyte growth, glucose homeostasis, and beta cell function. J Biol Chem. 2005, 280 (24): 23024-23031. 10.1074/jbc.M412564200.

    Article  CAS  PubMed  Google Scholar 

  11. Dahlman I, Nilsson M, Jiao H, Hoffstedt J, Lindgren CM, Humphreys K, Kere J, Gustafsson JA, Arner P, Dahlman-Wright K: Liver X receptor gene polymorphisms and adipose tissue expression levels in obesity. Pharmacogenet Genomics. 2006, 16 (12): 881-889. 10.1097/01.fpc.0000236334.49422.48.

    Article  CAS  PubMed  Google Scholar 

  12. Legry V, Cottel D, Ferrieres J, Chinetti G, Deroide T, Staels B, Amouyel P, Meirhaeghe A: Association between liver X receptor alpha gene polymorphisms and risk of metabolic syndrome in French populations. Int J Obes (Lond). 2008, 32 (3): 421-428. 10.1038/sj.ijo.0803705.

    Article  CAS  Google Scholar 

  13. Cheyssac C, Lecoeur C, Dechaume A, Bibi A, Charpentier G, Balkau B, Marre M, Froguel P, Gibson F, Vaxillaire M: Analysis of common PTPN1 gene variants in type 2 diabetes, obesity and associated phenotypes in the French population. BMC Med Genet. 2006, 7: 44-10.1186/1471-2350-7-44.

    Article  PubMed  PubMed Central  Google Scholar 

  14. World Health Organization Report: Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1 Diagnosis and classification of diabetes mellitus. Report of WHO consultation. 1999, [http://whqlibdoc.who.int/hq/1999/]

    Google Scholar 

  15. Vaxillaire M, Veslot J, Dina C, Proenca C, Cauchi S, Charpentier G, Tichet J, Fumeron F, Marre M, Meyre D, et al: Impact of common type 2 diabetes risk polymorphisms in the DESIR prospective study. Diabetes. 2008, 57 (1): 244-254. 10.2337/db07-0615.

    Article  CAS  PubMed  Google Scholar 

  16. Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, et al: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007, 445 (7130): 881-885. 10.1038/nature05616.

    Article  CAS  PubMed  Google Scholar 

  17. Florez JC, Manning AK, Dupuis J, McAteer J, Irenze K, Gianniny L, Mirel DB, Fox CS, Cupples LA, Meigs JB: A 100K genome-wide association scan for diabetes and related traits in the Framingham Heart Study: replication and integration with other genome-wide datasets. Diabetes. 2007, 56 (12): 3063-3074. 10.2337/db07-0451.

    Article  CAS  PubMed  Google Scholar 

  18. Hanson RL, Bogardus C, Duggan D, Kobes S, Knowlton M, Infante AM, Marovich L, Benitez D, Baier LJ, Knowler WC: A search for variants associated with young-onset type 2 diabetes in American Indians in a 100K genotyping array. Diabetes. 2007, 56 (12): 3045-3052. 10.2337/db07-0462.

    Article  CAS  PubMed  Google Scholar 

  19. Hayes MG, Pluzhnikov A, Miyake K, Sun Y, Ng MC, Roe CA, Below JE, Nicolae RI, Konkashbaev A, Bell GI, et al: Identification of type 2 diabetes genes in Mexican Americans through genome-wide association studies. Diabetes. 2007, 56 (12): 3033-3044. 10.2337/db07-0482.

    Article  CAS  PubMed  Google Scholar 

  20. Rampersaud E, Damcott CM, Fu M, Shen H, McArdle P, Shi X, Shelton J, Yin J, Chang YP, Ott SH, et al: Identification of novel candidate genes for type 2 diabetes from a genome-wide association scan in the Old Order Amish: evidence for replication from diabetes-related quantitative traits and from independent populations. Diabetes. 2007, 56 (12): 3053-3062. 10.2337/db07-0457.

    Article  CAS  PubMed  Google Scholar 

  21. Salonen JT, Uimari P, Aalto JM, Pirskanen M, Kaikkonen J, Todorova B, Hypponen J, Korhonen VP, Asikainen J, Devine C, et al: Type 2 diabetes whole-genome association study in four populations: the DiaGen consortium. Am J Hum Genet. 2007, 81 (2): 338-345. 10.1086/520599.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, et al: A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007, 316 (5829): 1341-1345. 10.1126/science.1142382.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Steinthorsdottir V, Thorleifsson G, Reynisdottir I, Benediktsson R, Jonsdottir T, Walters GB, Styrkarsdottir U, Gretarsdottir S, Emilsson V, Ghosh S, et al: A variant in CDKAL1 influences insulin response and risk of type 2 diabetes. Nat Genet. 2007, 39 (6): 770-775. 10.1038/ng2043.

    Article  CAS  PubMed  Google Scholar 

  24. Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, Timpson NJ, Perry JR, Rayner NW, Freathy RM, et al: Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science. 2007, 316 (5829): 1336-1341. 10.1126/science.1142364.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Gauderman WJ, Morrison JM: QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies. 2006, [http://hydra.usc.edu/gxe/]

    Google Scholar 

  26. Carlsson S, Persson PG, Alvarsson M, Efendic S, Norman A, Svanstrom L, Ostenson CG, Grill V: Weight history, glucose intolerance, and insulin levels in middle-aged Swedish men. Am J Epidemiol. 1998, 148 (6): 539-545.

    Article  CAS  PubMed  Google Scholar 

  27. Andersson CM, Bjaras GE, Ostenson CG: A stage model for assessing a community-based diabetes prevention program in Sweden. Health Promot Int. 2002, 17 (4): 317-327. 10.1093/heapro/17.4.317.

    Article  PubMed  Google Scholar 

  28. Agardh EE, Ahlbom A, Andersson T, Efendic S, Grill V, Hallqvist J, Norman A, Ostenson CG: Work stress and low sense of coherence is associated with type 2 diabetes in middle-aged Swedish women. Diabetes Care. 2003, 26 (3): 719-724. 10.2337/diacare.26.3.719.

    Article  PubMed  Google Scholar 

  29. Hilding A, Eriksson AK, Agardh EE, Grill V, Ahlbom A, Efendic S, Ostenson CG: The impact of family history of diabetes and lifestyle factors on abnormal glucose regulation in middle-aged Swedish men and women. Diabetologia. 2006, 49 (11): 2589-2598. 10.1007/s00125-006-0402-5.

    Article  CAS  PubMed  Google Scholar 

  30. Bonora E, Targher G, Alberiche M, Bonadonna RC, Saggiani F, Zenere MB, Monauni T, Muggeo M: Homeostasis model assessment closely mirrors the glucose clamp technique in the assessment of insulin sensitivity: studies in subjects with various degrees of glucose tolerance and insulin sensitivity. Diabetes Care. 2000, 23 (1): 57-63. 10.2337/diacare.23.1.57.

    Article  CAS  PubMed  Google Scholar 

  31. Jurinke C, Boom van den D, Cantor CR, Koster H: Automated genotyping using the DNA MassArray technology. Methods Mol Biol. 2002, 187: 179-192.

    CAS  PubMed  Google Scholar 

  32. Haploview. [http://www.broad.mit.edu/mpg/haploview/]

  33. software F. [http://ihg.gsf.de/cgi-bin/hw/hwa1.pl]

  34. Tregouet DA, Escolano S, Tiret L, Mallet A, Golmard JL: A new algorithm for haplotype-based association analysis: the Stochastic-EM algorithm. Ann Hum Genet. 2004, 68 (Pt 2): 165-177. 10.1046/j.1529-8817.2003.00085.x.

    Article  CAS  PubMed  Google Scholar 

  35. Nilsson M, Dahlman-Wright K, Karelmo C, Ebeling M, Gustafsson JA, Steffensen KR: Elk1 and SRF transcription factors convey basal transcription and mediate glucose response via their binding sites in the human LXRB gene promoter. Nucleic Acids Res. 2007, 35 (14): 4858-4868. 10.1093/nar/gkm492.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. (TESS) Ttess. [http://www.cbil.upenn.edu/tess/]

  37. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678. 10.1038/nature05911.

  38. Zeggini E, Scott LJ, Saxena R, Voight BF, Marchini JL, Hu T, de Bakker PI, Abecasis GR, Almgren P, Andersen G, et al: Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes. Nat Genet. 2008, 40 (5): 638-645. 10.1038/ng.120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, et al: Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007, 316 (5829): 1331-1336. 10.1126/science.1142358.

    Article  CAS  PubMed  Google Scholar 

  40. Steffensen KR, Gustafsson JA: Putative metabolic effects of the liver X receptor (LXR). Diabetes. 2004, 53 (Suppl 1): S36-42. 10.2337/diabetes.53.2007.S36.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This work was supported by The Swedish Diabetes Association, the Royal Physiographic Society in Lund, Stiftelsen Lars Hiertas Minne (both to MN), and the Swedish Research Council. We wish to thank Dr Beverley Balkau and the D.E.S.I.R. Study group for providing a proportion of the French normoglycemic control samples and phenotype data; and Dr Guillaume Charpentier from the Department of Endocrinology-Diabetology, Centre Hospitalier Sud-Francilien, Corbeil-Essonnes, France for his continuous support and collaboration in the recruitment of the diabetic patients.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Knut R Steffensen.

Additional information

Competing interests

Jan-Åke Gustafsson is consultant, shareholder and research grant receiver of KaroBio AB. The other authors declare that they have no competing interests.

Authors' contributions

ID analyzed the data and wrote the manuscript. MN designed the study and performed the in vitro studies. HFG, MV and CL participated in the design of the study, data analysis, and in the writing of the manuscript. CCP, SE, CGO, KB, GC and PF ascertained the patient cohorts. JÅG helped to draft the manuscript. KDW and KRS designed the study and helped to draft the manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Dahlman, I., Nilsson, M., Gu, H.F. et al. Functional and genetic analysis in type 2 diabetes of Liver X receptor alleles – a cohort study. BMC Med Genet 10, 27 (2009). https://doi.org/10.1186/1471-2350-10-27

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-10-27

Keywords