CVD is the leading cause of death in North America (Novak, Reference Novak1998; Heart and Stroke Foundation of Canada, 2003), emphasizing the need to understand its aetiology and to improve prevention and treatment strategies. Elevated plasma cholesterol concentrations are a well-known risk factor (Wilson et al. Reference Wilson, Garrison, Castelli, Feinleib, McNamara and Kannel1980; Grundy, Reference Grundy1995). The importance of dietary cholesterol in determining plasma cholesterol concentrations has been well documented (Dayton et al. Reference Dayton, Pearce, Goldman, Harnish, Plotkin, Shickman, Winfield, Zager and Dixon1968; Turpeinen, Reference Turpeinen1979; Clarke et al. Reference Clarke, Frost, Collins, Appleby and Peto1997; Hopkins, Reference Hopkins1992). There exists, however, a large interindividual variability in the plasma cholesterol response to changes in dietary cholesterol intake (Dreon & Krauss, Reference Dreon and Krauss1997; Grundy & Denke, Reference Grundy and Denke1990). Jacobs et al. (Reference Jacobs, Anderson, Hannan, Keys and Blackburn1983) have previously demonstrated that 9 % of subjects could be characterized as ‘hyperresponders’ when referring to the effect of a change in diet on plasma cholesterol levels, whereas 9 % were characterized as ‘hyporesponders’. Scientific evidence suggests that genetic factors might explain some of the interindividual variability in the plasma cholesterol response following modifications in cholesterol intake (Ordovas & Schaefer, Reference Ordovas and Schaefer1999).
The liver X receptor α (LXRα) is a transcription factor expressed predominantly in the liver but also in the kidney, intestine, macrophages, adipose tissue, spleen and adrenal glands (Kohro et al. Reference Kohro, Nakajima and Wada2000). LXRα regulates the expression of target genes by binding DNA sequence elements, termed LXR response elements (Willy et al. Reference Willy, Umesono, Ong, Evans, Heyman and Mangelsdorf1995). Several genes encoding proteins involved in cholesterol metabolism, for example SREBP-1, ABCA1, ABCG5 and ABCG8, are regulated by LXRα (Peet et al. Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998). In addition, oxysterols, which are important in steroid hormone biosynthesis, bile acid synthesis and the conversion of lanosterol to cholesterol, are also potent activators of LXRα, suggesting that LXRα may be an important sensor of cholesterol metabolites (Janowski et al. Reference Janowski, Willy, Devi, Falck and Mangelsdorf1996; Lehmann et al. Reference Lehmann, Kliewer and Moore1997) and might be involved in the interindividual variability observed in the plasma cholesterol response to changes in dietary cholesterol intake.
The first aim of the present study was first to identify genetic variants in the LXRα gene and evaluate their effects on the plasma lipoprotein/lipid profile. Second, we examined whether polymorphisms in the LXRα gene could modulate the association between cholesterol intake and plasma lipoprotein/lipid profile.
Subjects and methods
Subjects
The study sample originates from the Saguenay-Lac-St-Jean region, located in the north-eastern part of the province of Quebec. A total of 732 subjects (614 men, 118 women) were recruited through the Chicoutimi Hospital Lipid Clinic. Subjects were excluded if they were diagnosed with type 2 diabetes, type III dysbetalipoproteinaemia, familial hypercholesterolaemia or familial combined hyperlipidaemia. Type 2 diabetes was diagnosed according to the WHO criteria (Alberti & Zimmet, Reference Alberti and Zimmet1998). Written informed consent was obtained from all participating subjects, and the Medical Ethics Committee of Laval University and Chicoutimi Hospital Lipid Clinic approved the protocol.
Metabolic and anthropometric variables
Blood samples from subjects free of medication were obtained in the morning after a 12 h overnight fast. Blood was taken from an antecubital vein into vacutainer tubes containing EDTA. Blood samples were centrifuged within 1 h, and the plasma was frozen ( − 80°C) until analysis. Plasma cholesterol and triacylglycerol were measured using enzymatic assays (Burstein & Samaille, Reference Burstein and Samaille1960; McNamara & Schaefer, Reference McNamara and Schaefer1987). LDL-cholesterol was calculated using the Friedewald formula (Friedewald et al. Reference Friedewald, Levy and Frederickson1972). The HDL-cholesterol fraction was obtained after precipitation of LDL particles in the infranatant with heparin and MnCl2 (Burstein & Samaille, Reference Burstein and Samaille1960; Havel et al. Reference Havel, Eder and Bragdon1955). Apo B concentrations were measured in plasma by the rocket immunoelectrophoretic method of Laurell (Reference Laurell1966), as previously described (Moorjani et al. Reference Moorjani, Dupont, Labrie, Lupien, Brun, Gagné, Giguère and Bélanger1987). Serum standards were prepared and calibrated against reference sera obtained from the Center for Disease Control (Atlanta, GA, USA). Body weight, height and BMI were recorded.
DNA analysis
Genomic DNA was extracted using the Qiagen extraction kit (San Francisco, USA). The exons, the exon-intron splicing boundaries as well as the 5′ and 3′ regions of the human LXRα gene were sequenced to screen for DNA variants in thirty-five men exhibiting high plasma total-cholesterol and LDL-cholesterol concentrations (>5 and >3·5 mmol/l, respectively). Table 1 shows the primers used to amplify the different parts of the gene. Primers were designed using sequences available on GenBank (Accession numbers: AC024045) and Primer 3.0 software available on the Whitehead Institute/MIT Center for Genome Research server (Rozen & Skaletsky, Reference Rozen, Skaletsky, Krawetz and Misener2000). PCR conditions were as follows: reaction volume 50 μl, 1 U RedTaq DNA polymerase (Sigma; St Louis, MO, USA), 5 μl 10 × PCR buffer recommended by the manufacturer, 1·5 mm-MgCl2, 0·2 mm-dNTP, 8·4 μl of each primer at a final concentration of 7·5 μm, 0·1 μg genomic DNA. Betaine was added for sequencing promoter fragments 4, 3 and 1, and exon 4. The annealing temperature for each fragment is shown in Table 1.
Sequencing reactions were performed using BigDyeTerminator v3.0 cycle sequencing (ABI Prism; Applied Biosystems, Foster City, CA, USA) and the products were analysed on ABI 3100 automated DNA sequencer (PE Applied Biosystems). The gel files were processed using the ABI Prism3100 data collection software Applied Biosystems version 1.1 and ABI Prism DNA sequencing analysis software (PE Applied Biosystems) and then assembled and analysed using STADEN preGap4 and Gap4 software (http://staden.sourceforge.net/staden_home.html). Newly identified single nucleotide polymorphisms were genotyped using the same primers and same methods as for sequencing.
Nutritional assessment
A quantitative food-frequency questionnaire was used to evaluate dietary cholesterol intake and was available for 337 subjects. A trained dietitian administered this forty-eight-item questionnaire to participants. The dietitian asked the subjects to recall average use over the previous year. The frequency of food consumption was based on the number of times items were consumed per day, per week or per month. Nutrition Data System for Research version 4.02 (developed by the Nutrition Coordinating Center, University of Minnesota, Minneapolis, MN, USA; Food and Nutrient Database 30, published in November 1999) was used to calculate nutrients. The food-frequency questionnaire has been validated, in twenty-five men and women within the same population, using a 3 d food record (V Provencher et al. unpublished results). Furthermore, dietary cholesterol intake from food-frequency questionnaires correlated significantly with cholesterol intake as assessed by 3 d food records (r = 0·75, P = 0·0001).
Statistical analyses
Variables not normally distributed (triacylglycerol, fasting insulin, cholesterol and total-cholesterol:HDL-cholesterol ratio) were log10 transformed. To establish possible linkage disequilibrium between two polymorphisms, we used the software programs EH and 2LD (Linkage Utility Programs, Rockefeller University, USA: Mathematics & Statistics on the WWW). The Bayesian statistical method implemented in the program PHASE (version 2.1) was used to estimate haplotype frequencies in this study sample (Stephens et al. Reference Stephens, Smith and Donnelly2001). The association between polymorphisms in the LXRα gene and the plasma lipoprotein/lipid profile was evaluated by analyses of covariance with confounding variables (age and sex or age, sex and BMI) included in the model. Partial Pearson correlation coefficients were used to quantify the age-adjusted interrelationships between the plasma lipoprotein/lipid profile and dietary cholesterol intake in the entire cohort as well as within each genotype group. Gene–diet interactions were evaluated by analyses of covariance using general linear model procedures. Included in the model were the LXRα polymorphism, cholesterol intake, the interaction term and confounding variables (age, gender, energy intake). The source of variation in metabolic parameters was computed using the type III sum of squares. This sum of squares applies to unbalanced study designs and quantifies the effects of an independent variable after adjusting for all other variables included in the model. All the analyses were performed with SAS version 8.2 (SAS Institute, Cary, NC, USA). A value of P ≤ 0·05 was used to identify a statistically significant result.
Results
Molecular screening of the LXRα gene revealed the presence of sixteen polymorphisms (Fig. 1(A)), located principally in the 5′ and intronic regions. One single nucleotide polymorphism was located in exon 5, but it did not alter the predicted serine for this gene at position 99 (c.297C>T or p.Ser99Ser). The relative allele frequency of these variants is shown in Table 2. Database analysis of sequences containing these polymorphisms revealed potential binding sites for the transcription factors Sp-1, HES-1, USF and GATA-1 on the common -115G, -840C and -1830T alleles (Fig. 1(B); using Transcription Element Search Software; http://www.cbil.upenn.edu.tess). The likely functional polymorphisms -115G>A, -840C>A and -1830T>C were therefore selected for further analysis. The genotype distribution of each genetic variant was in Hardy–Weinberg equilibrium. We also observed that these polymorphisms were in perfect linkage disequilibrium (D′ = 1·0, P < 0·0001 for each pair of single nucleotide polymorphism). Haplotypes were further inferred using PHASE software. Owing to complete linkage disequilibrium between each pair of single nucleotide polymorphisms, two haplotypes were constructed. Therefore, haplotype analyses were not pursued any further.
Single nucleotide polymorphism s without ‘rs’ numbers were newly identified.
In order to explore the association of these polymorphisms in the LXRα gene with the plasma lipoprotein/lipid profile, their characteristics were compared between genotype groups (Table 3). Owing to their low number, homozygotes for the rare alleles (-115A, -840A, -1830C; n 19), were combined with heterozygotes (n 168) in the analyses. Analyses were performed combining results from men and women as there was no sex × genotype interaction for each variable tested (data not shown). We found that plasma total cholesterol concentrations were higher in carriers of the -115A allele than in -115G/G homozygotes after adjusting for age and sex (P = 0·05; Table 3. The results remained unchanged after adjusting for age, sex and BMI. When age, sex and BMI were included in the model, we observed a significant increase in the triacylglycerol level in carriers of the -115A allele compared with -115G/G homozygotes. Similar associations were found for the LXRα-840C>A and -1830T>C polymorphisms. We also observed similar trends when homozygotes for the rare allele and heterozygotes were analysed separately (data not shown).
* Mean values were significantly different from homozygotes after adjusting for age and sex (P < 0·05).
† Mean values were significantly different from homozygotes after adjusting for age, sex and BMI.
As cholesterol metabolites are potent ligands of LXRα, we examined whether polymorphisms in the LXRα gene could modulate the association between dietary cholesterol intake and plasma lipoprotein/lipid profile. Age-adjusted Pearson correlation coefficients were calculated. In the whole group, cholesterol intake was negatively correlated with plasma HDL-cholesterol concentration (r = − 0·16, P = 0·005) but was not correlated with the total cholesterol:HDL-cholesterol ratio or plasma triacylglycerol, apo B, total cholesterol or LDL-cholesterol concentrations (range r = 0·00 to r = 0·10, P = NS). When these analyses were performed for each genotype group, dietary cholesterol intake was correlated with plasma total cholesterol and LDL-cholesterol concentrations only in carriers of the rare LXRα polymorphisms -115G>A, -840C>A and -1830T>C (range r = 0·28 to r = 0·31, P < 0·05). No correlation with dietary cholesterol was observed in LXRα -115G/G, -840C/C or -1830T/T homozygotes (range r = − 0·10 to r = − 0·08, P = NS).
The interaction between dietary cholesterol and each polymorphism was therefore tested with plasma total cholesterol and LDL-cholesterol concentrations. In a model including the polymorphism, cholesterol intake, the interaction term (polymorphism × cholesterol (mg)), age, sex and energy intake, the LXRα-115G>A polymorphism and the interaction explained 1·82 % and 2·90 % of the variance in plasma total cholesterol concentration (Table 4). Furthermore, 2·12 % and 2·77 % of the variance in plasma LDL-cholesterol concentration was explained by the polymorphism and the interaction between the LXRα-115G>A polymorphism and dietary cholesterol, respectively. Similar results were obtained for the -840C>A and the -1830T>C polymorphism (Table 4). Including saturated fat intake, apo E genotype or BMI in the model did not modify the results.
* Age, sex and energy intake are included in the model.
In order to illustrate these interaction effects, cholesterol intake was separated into two subgroups according to the median value (290·8 mg; Fig. 2). We observed that carriers of the -115A allele were characterized by higher plasma total cholesterol and LDL-cholesterol concentrations when consuming a diet rich in cholesterol (>290·8 mg). However, in -115G/G homozygotes, plasma total cholesterol and LDL-cholesterol concentrations were similar irrespective of the amount of cholesterol provided from the diet. Similar results were obtained for the effects of the LXRα-840C>A and LXRα-1830T>C polymorphisms on plasma total cholesterol and LDL-cholesterol concentrations (Fig. 2).
Discussion
In the present study, sixteen genetic variants of the LXRα gene were identified. Among them, three polymorphic sites contained potential binding sites for transcription factors on the common allele. The LXRα-115G>A, -840C>A and -1830T>C polymorphisms were associated with moderate elevation of plasma cholesterol and triacylglycerol levels. Moreover, significant interactions between dietary cholesterol intake and the -115G>A, -840C>A and -1830T>C polymorphisms on plasma total cholesterol and LDL-cholesterol concentrations were observed, suggesting that the association between dietary cholesterol and the plasma lipoprotein/lipid profile was influenced by LXRα gene variants.
To our knowledge, the present study is the first to identify molecular variants in the gene encoding LXRα. Of the newly identified single nucleotide polymorphisms, three were located in regions where transcription factors recognized these sequences to regulate the transcription of LXRα. These sites seem to be abolished when the common allele is substituted for the rare allele at positions -115, -840 and -1830 on the LXRα sequence, suggesting a potential functional role for these genetic variants. In addition, all three polymorphisms are located in well-conserved regions in the mouse LXRα gene. Laffitte et al. (Reference Laffitte, Joseph, Walczak, Pei, Wilpitz, Collins and Tontonoz2001) have also shown that the promoter region of the human LXRα gene is similar to the mouse LXRα promoter region, suggesting that regions containing these variants are conserved between species and might thus be important in regulating the transcription of LXRα. Functional analyses must, however, be performed to test the functional significance of these three polymorphisms.
Several studies have demonstrated that LXRα is a key regulator of hepatic cholesterol metabolism. Hepatic cholesterol concentrations were increased 15–20-fold in LXRα− / − mice fed a high-cholesterol diet compared with wild-type mice on a similar diet (Peet et al. Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998). In addition, mice lacking the LXRα gene did not express the cholesterol 7α-hydroxylase gene, the rate-limiting enzyme in bile acid synthesis normally induced by LXRα (Peet et al. Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998), suggesting that LXRα acts as a cholesterol sensor to activate the cholesterol catabolism pathway. In the present study, the plasma lipoprotein/lipid profile was moderately modulated by the presence of molecular variants of the LXRα gene. LXRα also regulates the transcription of genes involved in fatty acid metabolism such as SREBP-1c and FAS, and has been demonstrated to play a role in triacylglycerol metabolism (Peet et al. Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998; Jakel et al. Reference Jakel, Nowak, Moitrot, Dehondt, Hum, Pennacchio, Fruchart-Najib and Fruchart2004). According to data obtained in the present study, variants within the LXRα gene promote the accumulation of plasma triacylglycerol, a finding that contrasts with previous results showing that LXRα activation induces lipogenesis. A possible interaction of LXRα with other genes or with dietary fat intake might explain this apparent contradiction as fatty acids induce the expression of LXRα (Tobin et al. Reference Tobin, Steineger, Alberti, Spydevold, Auwerx, Gustafsson and Nebb2000). Further studies are needed to confirm these results and to understand the effect of LXRα gene variants on plasma triacylglycerol levels.
Interest in studying the interaction between LXRα polymorphisms and dietary cholesterol intake is motivated by the high binding affinity of oxysterols for LXRα (Janowski et al. Reference Janowski, Willy, Devi, Falck and Mangelsdorf1996; Lehmann et al. Reference Lehmann, Kliewer and Moore1997). Some of the rate-limiting steps of important pathways in cholesterol metabolism are potent activators of LXRα (Janowski et al. Reference Janowski, Willy, Devi, Falck and Mangelsdorf1996). As shown by Peet et al. (Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998), mice lacking the LXRα gene (LXRα− / − ) were identical to wild-type mice when fed a standard chow diet. However, the hepatic and plasma cholesterol of LXRα− / − mice fed a diet rich in cholesterol increased dramatically (Peet et al. Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998) owing to decreased transcription of the gene encoding the rate-limiting enzyme in bile synthesis, cholesterol 7α-hydroxylase.
According to results obtained in the present study and by Peet et al. (Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998), we hypothesized that the impact of a change in LXRα transcriptional activity (owing to genetic variants) would be relatively small in a ‘normal-cholesterol’ environment. In the presence of excess of cholesterol from the diet, however, LXRα is not able to handle this surplus, leading to an increase in plasma total cholesterol and LDL-cholesterol concentrations. We could not rely on the absence of induction of cholesterol 7α-hydroxylase, as was seen by Peet et al. (Reference Peet, Turley, Ma, Janowski, Lobaccaro, Hammer and Mangelsdorf1998), to explain our results since cholesterol 7α-hydroxylase is not regulated by LXRα in human subjects (Menke et al. Reference Menke, Macnaul and Hayes2002; Chiang et al. Reference Chiang, Kimmel and Stroup2001). Considering the large number of genes involved in cholesterol metabolism that are regulated by LXRα, genetic variants affecting the transcriptional regulation of LXRα could still have an impact on the plasma lipid profile in the presence of a diet rich in cholesterol.
The impact of plasma total cholesterol and LDL-cholesterol concentrations on the risk of CVD is well established (Wilson et al. Reference Wilson, Garrison, Castelli, Feinleib, McNamara and Kannel1980; Grundy, Reference Grundy1995). Several strategies have been proposed to manage cholesterol concentration. Among non-pharmacological strategies, a reduction in dietary cholesterol intake appears to be associated with a decreased risk of CVD (Turpeinen, Reference Turpeinen1979; Dayton et al. Reference Dayton, Pearce, Goldman, Harnish, Plotkin, Shickman, Winfield, Zager and Dixon1968). A meta-analysis of well-controlled studies confirmed that a diet rich in cholesterol is associated with elevated plasma LDL-cholesterol concentrations (Clarke et al. Reference Clarke, Frost, Collins, Appleby and Peto1997; Hopkins, Reference Hopkins1992). However, the plasma cholesterol response to a change in cholesterol intake is highly variable and may be the result of polymorphisms in genes involved in cholesterol metabolism, such as LXRα. An identification of these genetic factors could help us to understand why studies examining the relationship between cholesterol intake and plasma cholesterol concentration are not always consistent. In the present study, we observed that three genetic variants in the LXRα gene modulated the relationship between cholesterol intake and plasma lipoprotein/lipid profile.
Other polymorphisms of genes involved in cholesterol metabolism have been shown to modulate the plasma lipid profile response to cholesterol intake. Among these, it is now well demonstrated that subjects with the apo E4 allele show a greater plasma cholesterol responsiveness to changes in cholesterol intake (Sarkkinen et al. Reference Sarkkinen, Korhonen, Erkkila, Ebeling and Uusitupa1998). In the present study, an interaction between dietary cholesterol and polymorphisms in the gene encoding LXRα was observed independently of the apo E genotype. In addition to cholesterol intake, physical activity could interact with genes involved in lipid metabolism to modulate total cholesterol and LDL-cholesterol concentrations (Boer et al. Reference Boer, Kuivenhoven, Feskens, Schouten, Havekes, Seidell, Kastelein and Kromhout1999; Taimela et al. Reference Taimela, Lehtimaki, Porkka, Rasanen and Viikari1996). Physical activity was not assessed in this study, but the level of physical activity as reported by Bruce and Katzmarzyk (Reference Bruce and Katzmarzyk2002) was relatively low. Indeed, 77 % of women and 74 % of men were physically inactive in Canada (Bruce & Katzmarzyk, Reference Bruce and Katzmarzyk2002).
Nutritional intervention studies will be needed to confirm the role of LXRα as a modulator of the plasma cholesterol response to cholesterol intake. A randomized and cross-over intervention study allows factors such as food intake to be controlled for. In the present study, cholesterol intake was self-reported and might lead to selective underreporting/overreporting. However, the food-frequency questionnaire has been widely used in epidemiological studies and represents a rapid and inexpensive tool for evaluating nutrient intake in large population studies (Willett et al. Reference Willett, Reynolds, Cottrell-Hoehner, Sampson and Browne1987). The presence of other environmental factors such as cigarette smoking might have influenced the observed effects. Statistical adjustment for this covariate did not attenuate the statistical significance of the gene–diet interaction effect (data not shown). Population admixture has been suggested to influence genetic association studies. In the present study, however, subjects were all of French-Canadian descent and population admixture is unlikely to account for the observed effect. Finally, we cannot rule out the possibility that an epistasis effect, not yet identified, might have influenced the effect of the gene–diet interaction.
In conclusion, the results of the present study suggest for the first time an association between human LXRα gene variations and plasma cholesterol concentrations in the presence of a diet rich in cholesterol. The identification of genetic variants of the LXRα gene provides helpful tools for further examining the role of LXRα in terms of the plasma lipoprotein/lipid profile. Although these results need confirmation in other studies, they underline the importance of gene–diet interaction effects in the aetiology of CVD.
Acknowledgements
The authors would like to express their gratitude to the subjects for their excellent collaboration. We would like to thank the staff of the CHUL Lipid Research Center and the Lipid Clinic, as well as of the Department of Biochemistry and the Cardiology Service of the Chicoutimi Hospital, for their dedicated support and assistance. This study was supported by a grant from the Canadian Institutes of Health Research (MOP-44 074) and the Heart and Stroke Foundation of Canada. J. R. received a doctoral studentship from the Canadian Institutes of Health Research. M.C. V. and S. L. are research scholars from the Fonds de la recherche en santé du Québec. D. G. is the holder of the Canada Research Chair in preventive genetics and community genomics (www.chaires.gc.ca).