CNR1 Genotype Influences HDL-Cholesterol Response to Change in Dietary Fat Intake

Abstract

Background

Success in further reducing the burden of cardiovascular disease (CVD) is threatened by the increasing prevalence of obesity-related atherogenic dyslipidemia. HDL-cholesterol (HDL-C) level is inversely correlated with CVD risk; each 1 mg/dl decrease in HDL-C is associated with a 6% reduction in risk. We previously showed that a common CNR1 haplotype, H3 (frequency 20%), is protective against the reduction in HDL-C that typically accompanies weight gain. In the present study, we extend that observation by reporting the effect of CNR1 haplotype on HDL-C response to modification of dietary fat intake in weight maintenance and weight loss.

Methods

Six haplotype tagging SNPs that cover the CNR1 gene locus were genotyped in 590 adults of varying body mass index (cohort 1 is 411 males with BMI 18.5–30.0 kg/m2; cohort 2 is 71 females with BMI18.5–30.0 kg/m2; and cohort 3 is 108 females with BMI 30–39.9 kg/m2). Dietary intakes were modified so that fat intake in the “high fat” condition was 15–20% greater than in the “low fat” condition, and lipid profiles were compared between carriers versus noncarriers for each of the five commonly observed CNR1 haplotypes (H1–H5).

Results

In normal to overweight subjects on eucaloric diets, the H3 haplotype was significantly associated with short-term high fat diet induced changes in HDL-C level in females (carriers 5.9 mg/dl>noncarriers, p = 0.007). The H3 haplotype was also significantly associated with HDL-C level after 16 weeks on high fat calorie restricted diet in obese females (carriers 6.8 mg/dl>noncarriers, p = 0.009).

Conclusion

Variability within the CNR1 gene locus contributes to gender-related differences in the HDL-cholesterol response to change in dietary fat intake. Functional characterization of this relationship in vitro may offer insights that potentially yield therapeutic guidance targeting dietary macronutrient composition, a direction much needed in the current epidemic of obesity.

About the Authors:

Heidi J. Silver

Contributed equally to this work with: Heidi J. Silver, Kevin D. Niswender

* E-mail: [email protected]

Affiliation: Division of Gastroenterology, Hepatology and Nutrition, Department of Medicine, Vanderbilt University School of Medicine, Nashville, Tennessee, United States of America

Kevin D. Niswender

Contributed equally to this work with: Heidi J. Silver, Kevin D. Niswender

Affiliations Tennessee Valley Healthcare System, Nashville, Tennessee, United States of America, Division of Diabetes, Endocrinology and Metabolism, Department of Medicine, Vanderbilt University School of Medicine, Nashville Tennessee, United States of America

Charles D. Keil

Affiliation: Division of Gastroenterology, Hepatology and Nutrition, Department of Medicine, Vanderbilt University School of Medicine, Nashville, Tennessee, United States of America

Lan Jiang

Affiliation: Center for Human Genetics Research, Department of Molecular Physiology and Biophysics, Vanderbilt University School of Medicine, Nashville, Tennessee, United States of America

Qiping Feng

Affiliation: Division of Clinical Pharmacology, Department of Medicine, Vanderbilt University School of Medicine, Nashville, Tennessee, United States of America

Sally Chiu

Affiliation: Children's Hospital Oakland Research Institute, Oakland, California, United States of America

Ronald M. Krauss

Affiliation: Children's Hospital Oakland Research Institute, Oakland, California, United States of America

Russell A. Wilke

Affiliation: Division of Clinical Pharmacology, Department of Medicine, Vanderbilt University School of Medicine, Nashville, Tennessee, United States of America

Introduction

Obesity is associated with an atherogenic lipoprotein profile characterized by increased very low-density lipoprotein triglycerides (VLDL-TG), small dense low density lipoprotein cholesterol (LDL-C), and decreased circulating high density lipoprotein cholesterol (HDL-C) [1]. The anti-inflammatory properties of HDL-C and its role in reverse cholesterol transport make it a significant modulator of risk for the development of cardiovascular disease (CVD). As event rate is inversely correlated with HDL-C concentration, each 1 mg/dl decrease in HDL-C is associated with a 6% reduction in CVD risk [2]. Decreased HDL-C is especially important in the current epidemic of obesity as higher body mass index (BMI) is associated with lower HDL-C levels in all racial/ethnic groups [3]. Notably, the prevalence of subthreshold (“at risk”) HDL-C levels has improved three-fold more in men than women (↓5.6% versus ↓1.8%) over the past 30 years [4]. In fact, the strongest clinically-relevant demographic predictors of plasma HDL-C concentration are BMI and gender [5], [6].

While modest weight loss of 5–10% improves plasma lipoprotein profile [7] and modifying dietary macronutrient intake can increase HDL-C level [8], scientific and clinical debate continues regarding the importance of the relative amount of fat in a calorically restricted diet [9], [10]. Low-fat (LF) diets can increase circulating VLDL-TG levels, reduce HDL-C, and thus, adversely affect atherogenic dyslipidemia [11], [12]. Although meta-analyses show greater improvements in HDL-C on high fat (HF) diets [13], [14], some studies yield conflicting results [15], [16]. Overall, the lipoprotein response to dietary macronutrient composition is a highly variable and complex trait [11], influenced by BMI [17], gender [18], and individual genetic characteristics [19], [20].

Preclinical and clinical studies targeting the endocannabinoid system (ECS) reveal potent effects on body weight, body fat distribution, and lipoprotein metabolism. The ECS consists primarily of two phospholipid-derived signaling molecules, N-arachidonoylethanolamide (anandamide) and 2-arachidonoylglycerol (2-AG), and two G protein-coupled receptors, CB1 and CB2 [21]. Most of the metabolic effects of ECS signaling appear to be mediated through CB1 receptors in central and peripheral tissues. In randomized clinical trials, CB1 receptor blockade resulted in weight loss, with increased HDL-C levels that exceeded the magnitude of change expected from weight loss alone [22], [23]. Collectively, these observations suggest that the ECS regulates lipid homeostasis via peripheral mechanisms that are at least partially independent of weight change.

Consistent with this hypothesis, we previously reported that sequence variation in CB1 (gene name CNR1) is associated with variability in circulating HDL-C levels independent of BMI [24], [25]. Importantly, a common CNR1 haplotype, H3, was further protective against the marked reduction in HDL-C that occurs in severe obesity (BMI≥40 kg/m2) [25]. However, it remains unknown whether the CNR1 H3 haplotype influences the effect of dietary fat intake on lipoprotein profiles. We tested this question in adults across a range of BMI who were prescribed various percentages of dietary fat intakes. Herein, we provide evidence that the response of HDL-C to increased dietary fat intake during weight maintenance is strongly influenced by the CNR1 H3 haplotype in women. We also report that this same H3 haplotype may protect women against low circulating HDL-C during voluntary weight loss.

Methods

Three distinct cohorts comprising 590 adults who participated in clinical trials designed to modify the percentage of dietary fat consumed were studied. Cohort #1 (normal to overweight males) and cohort #2 (normal to overweight females) were convenience samples from prior trials using eucaloric diets designed for weight maintenance. Because we had previously found an association between CNR1 and HDL-C in a population sample that was predominantly female [26], [27], we enrolled more males (cohort #1) than females (cohort #2) in the present study (at a ratio of 4∶1) to assure that we would not be missing a significant relationship in males. As the strength of the genotype-phenotype associations in cohorts #1 & #2 were consistent with our prior work, cohort #3 consisted of a sample of obese females on a calorically restricted diet, which not only broadened our range of BMI in females to better match current population demographics, but also allowed investigation of the effects of dietary fat intake within the context of weight loss.

For all subjects, the amount of fat prescribed in the “high fat” (HF) conditions were 15% higher than in the “low fat” (LF) conditions. Dietary intakes were assessed by Registered Dietitians (RDs) from 4-day food diaries and 24-hour diet recalls using multi-pass methodology [28]. Energy and macronutrient intake analysis was performed using Nutrition Data System for Research (NDS, Nutrition Coordinating Center, University of Minnesota, Minneapolis). Written informed consent was obtained and study protocols were approved by the Institutional Review Boards of Children's Hospital and Research Center of Oakland and the E.O. Lawrence Berkeley National Laboratory at the University of California for cohorts #1 & #2 and the Vanderbilt University Medical Center for cohort #3.

Cohort #1: Males on Eucaloric Diet

Males age 21–60 years with BMI 18.9–30 kg/m2 who had participated in one of four trials [26], [29], [30], [31], [32] were included if they had no active chronic disease, did not take lipid altering medications, and were nonsmokers. Males were randomly assigned in a crossover design to eucaloric HF diet (prescribed as 35% fat [∼18% saturated, 11% monounsaturated, 6% polyunsaturated], 50% carbohydrate, 15% protein) and LF diet (prescribed as 20% fat [∼5% saturated, 11% monounsaturated, 4% polyunsaturated], 65% carbohydrate [half simple and half complex], 15% protein) for 6 weeks each. Cholesterol content of the prescribed diets was 125 mg/1000 kcal and fiber content was 5 g/1000 kcal. Subjects maintained habitual physical activity and weighed themselves daily. RDs adjusted menu energy content when necessary to maintain weight. Diet compliance was assessed from 4-day weighed food intake records and daily recording of menu items consumed. Height, weight and venous blood samples were collected at baseline and end of study in an overnight fasted state. Plasma total cholesterol and triglyceride levels were determined by enzymatic procedures on a Gilford Impact 400E or Express 550 Plus analyzer (Ciba Corning, Oberlin, OH). HDL-C was measured after dextran sulfate precipitation of plasma and LDL-C was calculated from the Friedewald formula [33].

Cohort #2: Females on Eucaloric Diet

Premenopausal females age 21–50 years with BMI 18.9–30 kg/m2 were included if they had no active chronic disease, did not take lipid altering medications, were nonsmokers, and were not pregnant, breastfeeding or taking oral contraceptives [27]. Although not randomly crossed over, females were prescribed similar eucaloric HF and LF diets as the males in cohort #1 (above) for 8 week periods. The same procedures as described for cohort #1 were followed for compliance and outcome measurements.

Cohort #3: Females on Calorie Restricted Diet

Premenopausal females age 21–50 years with BMI 30–39.9 kg/m2 were included if they had no active chronic disease, did not take medications for lipid metabolism, were nonsmokers, and were not pregnant, breastfeeding or taking oral contraceptives. After habituation on the lower fat diet condition (35% fat), subjects were prescribed a 15% calorically restricted HF diet (50% fat [1/3 saturated, 1/3 monounsaturated, 1/3 polyunsaturated], 30% carbohydrate [half simple and half complex], 20% protein) for 16 weeks. Cholesterol content of the HF diet was 150 mg/1000 kcal and fiber content was 7 g/1000 kcal. Standard assays were performed for triglyceride, total cholesterol, LDL-C and HDL-C by selective enzymatic hydrolysis at the Vanderbilt Department of Pathology Clinical Laboratory. Serum glucose was quantified by colorimetric timed endpoint method and serum insulin by chemiluminescent immunoassay. The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated as [fasting glucose (mM) x fasting insulin (mU/L)]/22.5 [34].

Genotyping

Genomic DNA was extracted from whole blood according to standard procedures. Based on our prior work [24], [25], we tagged the entire CNR1 gene plus 5 kilobases (kb) upstream and 5 kb downstream using Haploview (http://www.broad.mit.edu/mpg/haploview/) for linkage disequilibrium analysis [35] and selection of tagSNPs. Gene locus, haplotype structure and allele frequencies of the tagSNPs covering the variability within our locus of interest (Chromosome 6, physical position 88, 901, 306–88, 916, 775) are shown in Figure 1. These six tagSNPs include two promoter SNPs (rs806370, rs806369), one coding SNP (rs1049353), two SNPs in the 3′ untranslated region (3′-UTR: rs12720071, rs806368), and one SNP distal to the polyA tail (rs806366). The tagSNPs were genotyped for cohorts #1 & #2 using multiplex polymerase chain-reaction (PCR), single base primer extension (SBE), and generation of mass spectra (Sequenom, San Diego, CA). For cohort #3, genotyping was conducted using commercial Taqman Allelic discrimination assays (Applied Biosystems, Inc., Foster City, CA) [36]. Each assay contained two Taqman probes uniquely labeled to bind separate (major versus minor) alleles.

[Figure omitted. See PDF.]

Figure 1. The CNR1 gene: chromosomal location, tagSNPs, haplotype structure and alleles.

*The CNR1 gene locus covers ∼15 kb on Chromosome 6 (chromosome 6 is shown in the reverse orientation, with arm q on the left). Six haplotype tagging SNPs were genotyped, based upon our prior work [25], and reference sequence numbers for the tag SNPs are presented. These SNPs are inherited in combination (haplotypes).

https://doi.org/10.1371/journal.pone.0036166.g001

Statistical Analysis

Statistical analyses were conducted using SPSS version 17 (SPSS Inc., Chicago, IL), JMP version 7.0 (SAS Institute Inc., Cary, NC) and PLINK genetic analyses toolset version 1.07 (http://pngu.mgh.harvard.edu/~purcell/plink) [37]. Data are expressed as means and standard deviations unless otherwise noted. Baseline lipid variables among cohorts were assessed using the Kruskal-Wallis test and relationships between demographics and lipids were examined using Spearman Rho correlations. Minor allele frequency and Hardy-Weinberg (H-W) equilibrium statistics were calculated for each SNP. As we have previously demonstrated that haplotypes provide greater power than single SNPs for identifying associations related to CNR1 [25], we performed haplotype analyses based on a dominant model (carrier versus noncarrier). Five CNR1 haplotypes (Figure 1) were identified at frequency over 5%. From 5′ to 3′, they are CCAATT (23%), CTGATC (30%), CCGATC (20%), TCGACT (13%), CCGGCT (7%). These five are the most common CNR1 haplotypes found in the Centre d'Etude du Polymorphisme Humain (CEPH) cohort within the International HapMap Consortium and in our prior cohorts [24]. Although the tagging SNPs were, by definition, in tight linkage dysequilibrium, we corrected for multiple hypothesis testing of the five haplotypes [38] and considered results statistically significant only with p value<0.01.

For cohorts #1 and #2, the five CNR1 haplotypes were estimated using PHASE software. As expected, the same five haplotypes were observed in cohort #3. However, the greater racial diversity of this cohort introduced more structural variability at the CNR1 gene locus, and thus, we identified two novel CNR1 haplotypes (CCGGTT, frequency 0.7%; and CCGATT, frequency 2%) in the African American (AA) females (data not presented). Importantly, these two new haplotypes were not significantly associated with HDL-C levels. For cohort #3, haplotypes were estimated by applying standard E-M algorithms located in the –hap option of PLINK.

All genotype-phenotype association analyses were performed using the five CNR1 haplotypes. In cohorts #1 and #2, linear regression models were created to examine relationships between haplotypes and lipid variables. Data were adjusted for age and baseline BMI. To be consistent with this approach, for cohort #3 a linear regression model was applied to test each haplotype (versus all other haplotypes) using the –hap-linear option in PLINK. For cohort #3, which was obese by definition, baseline BMI and HOMA-IR were significantly associated with baseline serum HDL-C, and thus, included as covariates in the model using the –covar option. This approach allowed us to replicate our initial findings from cohorts #1 &#2 in a second independent cohort (cohort #3) to reduce the likelihood of a spurious association between CNR1 and HDL-C levels.

Results

Subject Characteristics

The combined sample included 411 males and 179 females with a mean age of 39.9±4.0 years (Table 1). As defined by eligibility criteria, subjects in cohorts #1 & #2 were normal to overweight (mean BMI 24.8±1.7 kg/m2), while cohort #3 subjects were obese (mean BMI 34.7±2.7 kg/m2). Ninety percent of subjects in cohorts #1 & #2 self-identified as being of European ancestry (EA). In cohort #3, 75 (69%) females self-identified as EA and 33 (31%) as AA. Baseline TG, total cholesterol, LDL-C and HDL-C were similar across cohorts. There were no significant associations between minor allele frequency and diet adherence (data not presented).

[Figure omitted. See PDF.]

Table 1. Baseline characteristics of study subjects (n = 590).

https://doi.org/10.1371/journal.pone.0036166.t001

Cohort #1

As expected with eucaloric diets, no significant changes occurred in body weight in cohort #1 subjects. In these normal to overweight males, HF diet significantly decreased TG (−30.2±63.6 mg/dl, p<0.001), increased total cholesterol (+12.0±21.9 mg/dl, p<0.001), increased LDL-C (+11.5±20.6 mg/dl, p<0.001) and increased HDL-C (+6.4±5.7 mg/dl, p<0.001). Consistent with previous report in a predominantly female cohort [24], the reduction in TG was significantly greater in male noncarriers of the CNR1 H4 haplotype compared to male carriers (at least one copy), p = 0.005. Although prior work also indicated that the CNR1 H3 haplotype was associated with HDL-C in both genders, in the present study H3 was not significantly associated with HDL-C response to high fat diet in normal to overweight males (Table 2); male carriers of the H3 haplotype had only a 1.3 mg/dl higher HDL-C level (from 40.5±9.2 to 47.8±11.4 mg/dl) than male noncarriers (from 40.1±9.4 to 46.0±11.5 mg/dl).

[Figure omitted. See PDF.]

Table 2. Effect of haplotype on triglycerides, LDL-cholesterol and HDL-cholesterol in 411 normal to overweight males (Cohort #1) on eucaloric diets.

https://doi.org/10.1371/journal.pone.0036166.t002

Cohort #2

Like the males in cohort #1, females in cohort #2 experienced no significant changes in body weight. In this cohort, HF diet significantly decreased TG (−17.5±34.1 mg/dl, p<0.001), increased total cholesterol (+12.2±19.8 mg/dl, p<0.001), increased LDL-C (+4.8±17.0 mg/dl, p = 0.02) and increased HDL-C (+10.9±8.7 mg/dl, p<0.001). In contrast to the males in cohort #1, the CNR1 H3 haplotype was strongly associated (p = 0.007) with high fat diet induced changes in serum HDL-C levels in these normal to overweight females (Table 3). The increase in HDL-C level after HF diet in normal to overweight female carriers of the H3 haplotype (from 49±8 to 62.8±11.2 mg/dl) was 5.9 mg/dl greater than the increase in noncarriers (from 49.8±10.9 to 57.7±10.9 mg/dl) (Figure 2a).

[Figure omitted. See PDF.]

Figure 2. Effect of H3 haplotype on serum HDL-cholesterol levels in women.

Values shown are mean+standard error. In normal to overweight females, changing from low fat (LF) to high fat (HF) diet increased HDL-C from 49.0±8.0 mg/dl to 62.8±11.2 mg/dl in H3 haplotype carriers and from 49.8±10.9 mg/dl to 57.7±10.9 mg/dl in H3 non-carriers. In obese females, calorically restricted (−15% kcal) HF diet increased HDL-C from 50.11±11.60 mg/dl to 52.29±10.38 mg/dl in H3 haplotype carriers and decreased HDL-C from 47.99±12.90 mg/dl to 45.49±8.95 mg/dl in H3 non-carriers.

https://doi.org/10.1371/journal.pone.0036166.g002

[Figure omitted. See PDF.]

Table 3. Effect of haplotype on triglycerides, LDL-cholesterol and HDL-cholesterol in 71 normal to overweight females (cohort #2) on eucaloric diet.

https://doi.org/10.1371/journal.pone.0036166.t003

Cohort #3

The obese females experienced a significant weight loss of 6.2±3.5 kg over the 16-week intervention (p<0.001) due to the 15% caloric restriction prescribed in the HF diet condition. Overall, cohort #3 experienced significantly decreased BMI (−2.5±1.5 units, p<0.001), total cholesterol (−9.6±19.2 mg/dl, p<0.001), LDL-C (−8.0±17.6 mg/dl, p<0.001), plasma insulin (−2.03±4.79 mg/dl, p = 0.001), and had non-significant changes in TG (−5.58±33.6 mg/dl, p = 0.178), HDL-C (0.30±7.90 mg/dl, p = 0.758), and plasma glucose (−1.16±10.3 mg/dl, p = 0.360). As with normal to overweight females, the CNR1 H3 haplotype was strongly associated with serum HDL-C levels in obese females, p = 0.009 (Table 4). After 16 weeks of high fat diet, obese female carriers of the H3 haplotype had significantly higher HDL-C levels (52.3±2.3 mg/dl) than noncarriers (45.5±1.4 mg/dl) (Figure 2b). In re-testing the EA obese females alone, the H3 haplotype remained significantly associated (p = 0.005) with HDL-C level, and thus, the relationship was independent of race. The relationship between H3 haplotype and HDL-C was also independent of circulating TG level (p = 0.06 for the association between H3 haplotype and TG/HDL-C ratio).

[Figure omitted. See PDF.]

Table 4. Effect of haplotype on triglycerides, LDL-cholesterol and HDL-cholesterol in 108 obese females (cohort #3) on calorie restricted diet.a

https://doi.org/10.1371/journal.pone.0036166.t004

Discussion

The present findings demonstrate that variability in the CNR1 gene influences HDL-C response to changes in dietary fat intake. In normal to overweight adults on eucaloric diets, subjects expressing the CNR1 H3 haplotype had a more robust increase in HDL-C level during high fat intake and the magnitude of the effect of the H3 haplotype on HDL-C was four-fold greater in females compared to males. We observed a similar response to high fat intake in females with class I and II obesity who were prescribed caloric restriction. After 16 weeks of high fat diet, obese female carriers of the H3 haplotype had HDL-C levels that were ∼7 mg/dl greater than noncarriers (52.3±2.3 mg/dl in carriers versus 45.5±1.4 mg/dl in noncarriers, p = 0.009). The current findings are consistent with prior evidence that the CNR1 H3 haplotype was strongly associated with HDL-C levels in individuals with class III obesity [25].

It is important to recognize that the frequency of the H3 haplotype does not differ by gender and is common (20% frequency) in the general population. Given that women disproportionately seek dietary approaches for weight loss [39], along with the increasing recognition that women are at significant risk for cardiovascular disease - particularly in the current obesogenic environment, it is imperative to understand factors that determine favorable lipoprotein responses to dietary weight loss interventions. Prior investigation indicates that women may be more sensitive to diet induced changes in HDL-C during weight loss than men [40]. For example, a recent comparison of popular weight loss diets showed women had markedly increased HDL-C on 50% dietary fat intake (Atkins diet: ∼5 mg/dl increase in HDL-C), moderately increased HDL-C on 30% fat intake (Zone and Learn diets: ∼2.5 mg/dl increase in HDL-C) and no change in HDL-C on 10% fat intake (Ornish diet) [41]. The gender distinction observed in HDL-C response to popular weight loss diets may in part be attributable to underlying variability in the CNR1 gene.

The specific physiological mechanism linking the CNR1 gene product to lipid homeostasis is not completely understood. The CB1 receptor modulates food intake through an effect on central neuronal circuits [42]. Central CB1 receptor signaling is also known to directly influence peripheral metabolic processes involved in VLDL-TG secretion and clearance [43], [44]. Although changes in the abundance of VLDL could alter circulating HDL-C levels, by shifting the flux of free fatty acids between larger and smaller lipoproteins [45], our data do not support such a conclusion. In the present study, it appears that CNR1 modulated HDL-C response to high fat diet independent of BMI and fasting TG level. Peripherally, the CB1 receptor is expressed in the liver [46], pancreas, gastrointestinal tract [47] and adipose tissue [48]. HDL particle composition, and therefore HDL clearance rate, is modified within these tissues through the action of lipolytic enzymes [49]. Studies conducted using conditional (hepatic) CB1 knockout mice have shown that the CB1 receptor directly regulates hepatic lipid metabolism [50]. Lipoprotein lipase (LPL) and hepatic lipase (HL) convert larger HDL particles into smaller HDL remnants, resulting in HDL particles that are cleared more rapidly from the circulation. Endothelial lipase (EL) is also a strong predictor of circulating HDL-C concentration [51], [52], and EL expression is modulated by adipocytokines [53]. Because CB1 antagonists modulate the synthesis and secretion of adiponectin [54] and pro-inflammatory cytokines like TNFα [55], [56], genetic variability in CNR1 may be affecting HDL-C level through a mechanism involving its intravascular modeling. While our study was limited to assessment of HDL-C concentration, future investigation is needed to address more mechanistic hypotheses regarding the distribution of particle subfraction and changes in inflammatory biomarkers within the HDL proteome. In addition, large datasets are needed to condition the genetic effect of CNR1 on LPL, HL and EL, as well as known functional variants in the ATP-binding cassette transporter A1 (ABCA1), cholesteryl ester transfer protein (CETP), and scavenger receptor class B type 1 (SRB1). While all are known to influence HDL-C level [5], [51], [57], combinatorial models assessing interactions [58] will require targeted genotyping within cohorts of enormous size.

As these studies proceed, it is important to recognize that the causative allele linking CNR1 to HDL-C level remains undetermined. It is likely that the H3 haplotype examined in the present study is only in partial linkage (ie, D′<1.0) with the actual functional allele, a group of common and rare variants with the capacity to alter CNR1 gene function [59], [60]. To fully characterize these variants, future investigation must target deep re-sequencing of the CNR1 gene in individuals expressing H3, from cohorts of racially diverse populations enriched for the traits of interest. In conclusion, we have demonstrated that the common H3 haplotype at the CNR1 gene locus potently influences HDL-C response to change in dietary fat intake in a gender-specific manner. Functional characterization of this relationship may further our understanding of the biology underlying lipid homeostasis, and potentially yield therapeutic guidance to improve public health by optimizing dietary macronutrient modifications, a direction much needed in the current epidemic of obesity.

Acknowledgments

The authors thank Taren B Spence, R.D. for study implementation and Will Bush, Ph.D. for assistance with manuscript preparation. In addition, we thank the following for donating study food supplies: Atkins Nutritionals, Inc; Oregon Hazelnut Commission; National Peanut Board; Setton Farms; National Sunflower Association; Almond Board of California.

Author Contributions

Conceived and designed the experiments: HJS KDN RAW. Performed the experiments: HJS CDK QF SC. Analyzed the data: HJS CDK QF SC LJ KDN RMK RAW. Contributed reagents/materials/analysis tools: RAW. Wrote the paper: HJS CDK QF SC KDN RMK RAW. Administrative and technical support: CDK QF SC LJ.

Citation: Silver HJ, Niswender KD, Keil CD, Jiang L, Feng Q, Chiu S, et al. (2012) CNR1 Genotype Influences HDL-Cholesterol Response to Change in Dietary Fat Intake. PLoS ONE7(5): e36166. https://doi.org/10.1371/journal.pone.0036166

References

1. Ford ES, Giles WH, Dietz WH (2002) Prevalence of the metabolic syndrome among US adults: findings from the third National Health and Nutrition Examination Survey. Jama 287: 356–359.ES FordWH GilesWH Dietz2002Prevalence of the metabolic syndrome among US adults: findings from the third National Health and Nutrition Examination Survey.Jama287356359

2. Ashen MD, Blumenthal RS (2005) Clinical practice. Low HDL cholesterol levels. N Engl J Med 353: 1252–1260.MD AshenRS Blumenthal2005Clinical practice. Low HDL cholesterol levels.N Engl J Med35312521260

3. Brown CD, Higgins M, Donato KA, Rohde FC, Garrison R, et al. (2000) Body mass index and the prevalence of hypertension and dyslipidemia. Obes Res 8: 605–619.CD BrownM. HigginsKA DonatoFC RohdeR. Garrison2000Body mass index and the prevalence of hypertension and dyslipidemia.Obes Res8605619

4. Kim JK, Alley D, Seeman T, Karlamangla A, Crimmins E (2006) Recent changes in cardiovascular risk factors among women and men. J Womens Health (Larchmt) 15: 734–746.JK KimD. AlleyT. SeemanA. KarlamanglaE. Crimmins2006Recent changes in cardiovascular risk factors among women and men.J Womens Health (Larchmt)15734746

5. Wilke RA (2011) High-density lipoprotein (HDL) cholesterol: leveraging practice-based biobank cohorts to characterize clinical and genetic predictors of treatment outcome. Pharmacogenomics J 11: 162–173.RA Wilke2011High-density lipoprotein (HDL) cholesterol: leveraging practice-based biobank cohorts to characterize clinical and genetic predictors of treatment outcome.Pharmacogenomics J11162173

6. Wilke RA, Berg RL, Linneman JG, Peissig P, Starren J, et al. (2010) Quantification of the clinical modifiers impacting high-density lipoprotein cholesterol in the community: Personalized Medicine Research Project. Prev Cardiol 13: 63–68.RA WilkeRL BergJG LinnemanP. PeissigJ. Starren2010Quantification of the clinical modifiers impacting high-density lipoprotein cholesterol in the community: Personalized Medicine Research Project.Prev Cardiol136368

7. Kanders B, Peterson FJ, Lavin P (1994) Long term health effects associated with significant weight loss: A study of the dose response effect. In: Kanders GLBaBS, editor. Obesity Pathophysiology, Psychology and Treatment. New York: Chapman & Hall. pp. 167–181.B. KandersFJ PetersonP. Lavin1994Long term health effects associated with significant weight loss: A study of the dose response effect.GLBaBS KandersObesity Pathophysiology, Psychology and TreatmentNew YorkChapman & Hall167181

8. Stein O, Stein Y (1999) Atheroprotective mechanisms of HDL. Atherosclerosis 144: 285–301.O. SteinY. Stein1999Atheroprotective mechanisms of HDL.Atherosclerosis144285301

9. Bray GA (2011) Is dietary fat important? Am J Clin Nutr 93: 481–482.GA Bray2011Is dietary fat important?Am J Clin Nutr93481482

10. Willett WC (1998) Dietary fat and obesity: an unconvincing relation. Am J Clin Nutr 68: 1149–1150.WC Willett1998Dietary fat and obesity: an unconvincing relation.Am J Clin Nutr6811491150

11. Siri-Tarino PW, Sun Q, Hu FB, Krauss RM (2010) Meta-analysis of prospective cohort studies evaluating the association of saturated fat with cardiovascular disease. Am J Clin Nutr 91: 535–546.PW Siri-TarinoQ. SunFB HuRM Krauss2010Meta-analysis of prospective cohort studies evaluating the association of saturated fat with cardiovascular disease.Am J Clin Nutr91535546

12. Sieri S, Krogh V, Berrino F, Evangelista A, Agnoli C, et al. (2010) Dietary glycemic load and index and risk of coronary heart disease in a large italian cohort: the EPICOR study. Arch Intern Med 170: 640–647.S. SieriV. KroghF. BerrinoA. EvangelistaC. Agnoli2010Dietary glycemic load and index and risk of coronary heart disease in a large italian cohort: the EPICOR study.Arch Intern Med170640647

13. Sacks FM, Katan M (2002) Randomized clinical trials on the effects of dietary fat and carbohydrate on plasma lipoproteins and cardiovascular disease. Am J Med 113: Suppl 9B13S–24S.FM SacksM. Katan2002Randomized clinical trials on the effects of dietary fat and carbohydrate on plasma lipoproteins and cardiovascular disease.Am J Med113Suppl 9B13S24S

14. Nordmann AJ, Nordmann A, Briel M, Keller U, Yancy WS Jr, et al. (2006) Effects of low-carbohydrate vs low-fat diets on weight loss and cardiovascular risk factors: a meta-analysis of randomized controlled trials. Arch Intern Med 166: 285–293.AJ NordmannA. NordmannM. BrielU. KellerWS Yancy Jr2006Effects of low-carbohydrate vs low-fat diets on weight loss and cardiovascular risk factors: a meta-analysis of randomized controlled trials.Arch Intern Med166285293

15. Mensink RP, Zock PL, Kester AD, Katan MB (2003) Effects of dietary fatty acids and carbohydrates on the ratio of serum total to HDL cholesterol and on serum lipids and apolipoproteins: a meta-analysis of 60 controlled trials. Am J Clin Nutr 77: 1146–1155.RP MensinkPL ZockAD KesterMB Katan2003Effects of dietary fatty acids and carbohydrates on the ratio of serum total to HDL cholesterol and on serum lipids and apolipoproteins: a meta-analysis of 60 controlled trials.Am J Clin Nutr7711461155

16. Volek JS, Sharman MJ, Gomez AL, DiPasquale C, Roti M, et al. (2004) Comparison of a very low-carbohydrate and low-fat diet on fasting lipids, LDL subclasses, insulin resistance, and postprandial lipemic responses in overweight women. J Am Coll Nutr 23: 177–184.JS VolekMJ SharmanAL GomezC. DiPasqualeM. Roti2004Comparison of a very low-carbohydrate and low-fat diet on fasting lipids, LDL subclasses, insulin resistance, and postprandial lipemic responses in overweight women.J Am Coll Nutr23177184

17. Jansen S, Lopez-Miranda J, Salas J, Castro P, Paniagua JA, et al. (1998) Plasma lipid response to hypolipidemic diets in young healthy non-obese men varies with body mass index. J Nutr 128: 1144–1149.S. JansenJ. Lopez-MirandaJ. SalasP. CastroJA Paniagua1998Plasma lipid response to hypolipidemic diets in young healthy non-obese men varies with body mass index.J Nutr12811441149

18. Weggemans RM, Zock PL, Urgert R, Katan MB (1999) Differences between men and women in the response of serum cholesterol to dietary changes. Eur J Clin Invest 29: 827–834.RM WeggemansPL ZockR. UrgertMB Katan1999Differences between men and women in the response of serum cholesterol to dietary changes.Eur J Clin Invest29827834

19. Denke MA, Adams-Huet B, Nguyen AT (2000) Individual cholesterol variation in response to a margarine- or butter-based diet: A study in families. JAMA 284: 2740–2747.MA DenkeB. Adams-HuetAT Nguyen2000Individual cholesterol variation in response to a margarine- or butter-based diet: A study in families.JAMA28427402747

20. Ordovas JM (2006) Nutrigenetics, plasma lipids, and cardiovascular risk. J Am Diet Assoc 106: 1074–1081.JM Ordovas2006Nutrigenetics, plasma lipids, and cardiovascular risk.J Am Diet Assoc10610741081quiz 1083. quiz 1083.

21. Piomelli D (2003) The molecular logic of endocannabinoid signalling. Nat Rev Neurosci 4: 873–884.D. Piomelli2003The molecular logic of endocannabinoid signalling.Nat Rev Neurosci4873884

22. Despres JP, Golay A, Sjostrom L (2005) Effects of rimonabant on metabolic risk factors in overweight patients with dyslipidemia. N Engl J Med 353: 2121–2134.JP DespresA. GolayL. Sjostrom2005Effects of rimonabant on metabolic risk factors in overweight patients with dyslipidemia.N Engl J Med35321212134

23. Ruilope LM, Despres JP, Scheen A, Pi-Sunyer X, Mancia G, et al. (2008) Effect of rimonabant on blood pressure in overweight/obese patients with/without co-morbidities: analysis of pooled RIO study results. J Hypertens 26: 357–367.LM RuilopeJP DespresA. ScheenX. Pi-SunyerG. Mancia2008Effect of rimonabant on blood pressure in overweight/obese patients with/without co-morbidities: analysis of pooled RIO study results.J Hypertens26357367

24. Baye TM, Zhang Y, Smith E, Hillard CJ, Gunnell J, et al. (2008) Genetic variation in cannabinoid receptor 1 (CNR1) is associated with derangements in lipid homeostasis, independent of body mass index. Pharmacogenomics 9: 1647–1656.TM BayeY. ZhangE. SmithCJ HillardJ. Gunnell2008Genetic variation in cannabinoid receptor 1 (CNR1) is associated with derangements in lipid homeostasis, independent of body mass index.Pharmacogenomics916471656

25. Feng Q, Jiang L, Berg RL, Antonik M, MacKinney E, et al. (2010) A common CNR1 (cannabinoid receptor 1) haplotype attenuates the decrease in HDL cholesterol that typically accompanies weight gain. PLoS One 5: e15779.Q. FengL. JiangRL BergM. AntonikE. MacKinney2010A common CNR1 (cannabinoid receptor 1) haplotype attenuates the decrease in HDL cholesterol that typically accompanies weight gain.PLoS One5e15779

26. Krauss RM, Dreon DM (1995) Low-density-lipoprotein subclasses and response to a low-fat diet in healthy men. Am J Clin Nutr 62: 478S–487S.RM KraussDM Dreon1995Low-density-lipoprotein subclasses and response to a low-fat diet in healthy men.Am J Clin Nutr62478S487S

27. Dreon DM, Fernstrom HA, Williams PT, Krauss RM (1997) LDL subclass patterns and lipoprotein response to a low-fat, high-carbohydrate diet in women. Arterioscler Thromb Vasc Biol 17: 707–714.DM DreonHA FernstromPT WilliamsRM Krauss1997LDL subclass patterns and lipoprotein response to a low-fat, high-carbohydrate diet in women.Arterioscler Thromb Vasc Biol17707714

28. Moshfegh AJ, Borrud L, Perloff B, LaComb R (1999) Improved Method for the 24-hour Dietary Recall for Use in National Surveys. FASEB Journal 13: A603.AJ MoshfeghL. BorrudB. PerloffR. LaComb1999Improved Method for the 24-hour Dietary Recall for Use in National Surveys.FASEB Journal13A603

29. Dreon DM, Fernstrom HA, Campos H, Blanche P, Williams PT, et al. (1998) Change in dietary saturated fat intake is correlated with change in mass of large low-density-lipoprotein particles in men. Am J Clin Nutr 67: 828–836.DM DreonHA FernstromH. CamposP. BlanchePT Williams1998Change in dietary saturated fat intake is correlated with change in mass of large low-density-lipoprotein particles in men.Am J Clin Nutr67828836

30. Dreon DM, Krauss RM (1997) Diet-gene interactions in human lipoprotein metabolism. J Am Coll Nutr 16: 313–324.DM DreonRM Krauss1997Diet-gene interactions in human lipoprotein metabolism.J Am Coll Nutr16313324

31. Krauss RM (1997) Genetic, metabolic, and dietary influences on the atherogenic lipoprotein phenotype. World Rev Nutr Diet 80: 22–43.RM Krauss1997Genetic, metabolic, and dietary influences on the atherogenic lipoprotein phenotype.World Rev Nutr Diet802243

32. Pennacchio LA, Olivier M, Hubacek JA, Cohen JC, Cox DR, et al. (2001) An apolipoprotein influencing triglycerides in humans and mice revealed by comparative sequencing. Science 294: 169–173.LA PennacchioM. OlivierJA HubacekJC CohenDR Cox2001An apolipoprotein influencing triglycerides in humans and mice revealed by comparative sequencing.Science294169173

33. Friedewald WT, Levy RI, Fredrickson DS (1972) Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 18: 499–502.WT FriedewaldRI LevyDS Fredrickson1972Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge.Clin Chem18499502

34. Wallace TM, Levy JC, Matthews DR (2004) Use and abuse of HOMA modeling. Diabetes Care 27: 1487–1495.TM WallaceJC LevyDR Matthews2004Use and abuse of HOMA modeling.Diabetes Care2714871495

35. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229.SB GabrielSF SchaffnerH. NguyenJM MooreJ. Roy2002The structure of haplotype blocks in the human genome.Science29622252229

36. Holland PM, Abramson RD, Watson R, Gelfand DH (1991) Detection of specific polymerase chain reaction product by utilizing the 5′----3′ exonuclease activity of Thermus aquaticus DNA polymerase. Proc Natl Acad Sci U S A 88: 7276–7280.PM HollandRD AbramsonR. WatsonDH Gelfand1991Detection of specific polymerase chain reaction product by utilizing the 5′----3′ exonuclease activity of Thermus aquaticus DNA polymerase.Proc Natl Acad Sci U S A8872767280

37. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

38. Russo P, Strazzullo P, Cappuccio FP, Tregouet DA, Lauria F, et al. (2007) Genetic variations at the endocannabinoid type 1 receptor gene (CNR1) are associated with obesity phenotypes in men. J Clin Endocrinol Metab 92: 2382–2386.P. RussoP. StrazzulloFP CappuccioDA TregouetF. Lauria2007Genetic variations at the endocannabinoid type 1 receptor gene (CNR1) are associated with obesity phenotypes in men.J Clin Endocrinol Metab9223822386

39. Serdula MK, Mokdad AH, Williamson DF, Galuska DA, Mendlein JM, et al. (1999) Prevalence of attempting weight loss and strategies for controlling weight. Jama 282: 1353–1358.MK SerdulaAH MokdadDF WilliamsonDA GaluskaJM Mendlein1999Prevalence of attempting weight loss and strategies for controlling weight.Jama28213531358

40. Sharman MJ, Gomez AL, Kraemer WJ, Volek JS (2004) Very low-carbohydrate and low-fat diets affect fasting lipids and postprandial lipemia differently in overweight men. J Nutr 134: 880–885.MJ SharmanAL GomezWJ KraemerJS Volek2004Very low-carbohydrate and low-fat diets affect fasting lipids and postprandial lipemia differently in overweight men.J Nutr134880885

41. Gardner CD, Kiazand A, Alhassan S, Kim S, Stafford RS, et al. (2007) Comparison of the Atkins, Zone, Ornish, and LEARN diets for change in weight and related risk factors among overweight premenopausal women: the A TO Z Weight Loss Study: a randomized trial. Jama 297: 969–977.CD GardnerA. KiazandS. AlhassanS. KimRS Stafford2007Comparison of the Atkins, Zone, Ornish, and LEARN diets for change in weight and related risk factors among overweight premenopausal women: the A TO Z Weight Loss Study: a randomized trial.Jama297969977

42. Di Marzo V, Matias I (2005) Endocannabinoid control of food intake and energy balance. Nat Neurosci 8: 585–589.V. Di MarzoI. Matias2005Endocannabinoid control of food intake and energy balance.Nat Neurosci8585589

43. Ruby MA, Nomura DK, Hudak CS, Mangravite LM, Chiu S, et al. (2008) Overactive endocannabinoid signaling impairs apolipoprotein E-mediated clearance of triglyceride-rich lipoproteins. Proc Natl Acad Sci U S A 105: 14561–14566.MA RubyDK NomuraCS HudakLM MangraviteS. Chiu2008Overactive endocannabinoid signaling impairs apolipoprotein E-mediated clearance of triglyceride-rich lipoproteins.Proc Natl Acad Sci U S A1051456114566

44. Stafford JM, Yu F, Printz R, Hasty AH, Swift LL, et al. (2008) Central nervous system neuropeptide Y signaling modulates VLDL triglyceride secretion. Diabetes 57: 1482–1490.JM StaffordF. YuR. PrintzAH HastyLL Swift2008Central nervous system neuropeptide Y signaling modulates VLDL triglyceride secretion.Diabetes5714821490

45. Ginsberg HN (1998) Lipoprotein physiology. Endocrinol Metab Clin North Am 27: 503–519.HN Ginsberg1998Lipoprotein physiology.Endocrinol Metab Clin North Am27503519

46. Osei-Hyiaman D, DePetrillo M, Pacher P, Liu J, Radaeva S, et al. (2005) Endocannabinoid activation at hepatic CB1 receptors stimulates fatty acid synthesis and contributes to diet-induced obesity. J Clin Invest 115: 1298–1305.D. Osei-HyiamanM. DePetrilloP. PacherJ. LiuS. Radaeva2005Endocannabinoid activation at hepatic CB1 receptors stimulates fatty acid synthesis and contributes to diet-induced obesity.J Clin Invest11512981305

47. Tyler K, Hillard CJ, Greenwood-Van Meerveld B (2000) Inhibition of small intestinal secretion by cannabinoids is CB1 receptor-mediated in rats. Eur J Pharmacol 409: 207–211.K. TylerCJ HillardB. Greenwood-Van Meerveld2000Inhibition of small intestinal secretion by cannabinoids is CB1 receptor-mediated in rats.Eur J Pharmacol409207211

48. Bensaid M, Gary-Bobo M, Esclangon A, Maffrand JP, Le Fur G, et al. (2003) The cannabinoid CB1 receptor antagonist SR141716 increases Acrp30 mRNA expression in adipose tissue of obese fa/fa rats and in cultured adipocyte cells. Mol Pharmacol 63: 908–914.M. BensaidM. Gary-BoboA. EsclangonJP MaffrandG. Le Fur2003The cannabinoid CB1 receptor antagonist SR141716 increases Acrp30 mRNA expression in adipose tissue of obese fa/fa rats and in cultured adipocyte cells.Mol Pharmacol63908914

49. Lewis GF, Rader DJ (2005) New insights into the regulation of HDL metabolism and reverse cholesterol transport. Circ Res 96: 1221–1232.GF LewisDJ Rader2005New insights into the regulation of HDL metabolism and reverse cholesterol transport.Circ Res9612211232

50. Osei-Hyiaman D, Liu J, Zhou L, Godlewski G, Harvey-White J, et al. (2008) Hepatic CB1 receptor is required for development of diet-induced steatosis, dyslipidemia, and insulin and leptin resistance in mice. J Clin Invest 118: 3160–3169.D. Osei-HyiamanJ. LiuL. ZhouG. GodlewskiJ. Harvey-White2008Hepatic CB1 receptor is required for development of diet-induced steatosis, dyslipidemia, and insulin and leptin resistance in mice.J Clin Invest11831603169

51. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, et al. (2010) Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466: 707–713.TM TeslovichK. MusunuruAV SmithAC EdmondsonIM Stylianou2010Biological, clinical and population relevance of 95 loci for blood lipids.Nature466707713

52. Yasuda T, Ishida T, Rader DJ (2010) Update on the role of endothelial lipase in high-density lipoprotein metabolism, reverse cholesterol transport, and atherosclerosis. Circ J 74: 2263–2270.T. YasudaT. IshidaDJ Rader2010Update on the role of endothelial lipase in high-density lipoprotein metabolism, reverse cholesterol transport, and atherosclerosis.Circ J7422632270

53. Badellino KO, Wolfe ML, Reilly MP, Rader DJ (2008) Endothelial lipase is increased in vivo by inflammation in humans. Circulation 117: 678–685.KO BadellinoML WolfeMP ReillyDJ Rader2008Endothelial lipase is increased in vivo by inflammation in humans.Circulation117678685

54. Gary-Bobo M, Elachouri G, Scatton B, Le Fur G, Oury-Donat F, et al. (2006) The cannabinoid CB1 receptor antagonist rimonabant (SR141716) inhibits cell proliferation and increases markers of adipocyte maturation in cultured mouse 3T3 F442A preadipocytes. Mol Pharmacol 69: 471–478.M. Gary-BoboG. ElachouriB. ScattonG. Le FurF. Oury-Donat2006The cannabinoid CB1 receptor antagonist rimonabant (SR141716) inhibits cell proliferation and increases markers of adipocyte maturation in cultured mouse 3T3 F442A preadipocytes.Mol Pharmacol69471478

55. Gary-Bobo M, Elachouri G, Gallas JF, Janiak P, Marini P, et al. (2007) Rimonabant reduces obesity-associated hepatic steatosis and features of metabolic syndrome in obese Zucker fa/fa rats. Hepatology 46: 122–129.M. Gary-BoboG. ElachouriJF GallasP. JaniakP. Marini2007Rimonabant reduces obesity-associated hepatic steatosis and features of metabolic syndrome in obese Zucker fa/fa rats.Hepatology46122129

56. Jourdan T, Djaouti L, Demizieux L, Gresti J, Verges B, et al. (2010) CB1 antagonism exerts specific molecular effects on visceral and subcutaneous fat and reverses liver steatosis in diet-induced obese mice. Diabetes 59: 926–934.T. JourdanL. DjaoutiL. DemizieuxJ. GrestiB. Verges2010CB1 antagonism exerts specific molecular effects on visceral and subcutaneous fat and reverses liver steatosis in diet-induced obese mice.Diabetes59926934

57. Knoblauch H, Bauerfeind A, Toliat MR, Becker C, Luganskaja T, et al. (2004) Haplotypes and SNPs in 13 lipid-relevant genes explain most of the genetic variance in high-density lipoprotein and low-density lipoprotein cholesterol. Hum Mol Genet 13: 993–1004.H. KnoblauchA. BauerfeindMR ToliatC. BeckerT. Luganskaja2004Haplotypes and SNPs in 13 lipid-relevant genes explain most of the genetic variance in high-density lipoprotein and low-density lipoprotein cholesterol.Hum Mol Genet139931004

58. Turner SD, Berg RL, Linneman JG, Peissig PL, Crawford DC, et al. (2011) Knowledge-driven multi-locus analysis reveals gene-gene interactions influencing HDL cholesterol level in two independent EMR-linked biobanks. PLoS One 6: e19586.SD TurnerRL BergJG LinnemanPL PeissigDC Crawford2011Knowledge-driven multi-locus analysis reveals gene-gene interactions influencing HDL cholesterol level in two independent EMR-linked biobanks.PLoS One6e19586

59. Goldstein DB (2009) Common genetic variation and human traits. N Engl J Med 360: 1696–1698.DB Goldstein2009Common genetic variation and human traits.N Engl J Med36016961698

60. Manolio TA (2010) Genomewide association studies and assessment of the risk of disease. N Engl J Med 363: 166–176.TA Manolio2010Genomewide association studies and assessment of the risk of disease.N Engl J Med363166176