Skip to main content

Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham Heart Study data



Cholesterol concentrations in blood are related to cardiovascular diseases. Recent genome-wide association studies (GWAS) of cholesterol levels identified a number of single-locus effects on total cholesterol (TC) and high-density lipoprotein cholesterol (HDL-C) levels. Here, we report single-locus and epistasis SNP effects on TC and HDL-C using the Framingham Heart Study (FHS) data.


Single-locus effects and pairwise epistasis effects of 432,096 SNP markers were tested for their significance on log-transformed TC and HDL-C levels. Twenty nine additive SNP effects reached single-locus genome-wide significance (p < 7.2 × 10-8) and no dominance effect reached genome-wide significance. Two new gene regions were detected, the RAB3GAP1-R3HDM1-LCT-MCM6 region of chr02 for TC identified by six new SNPs, and the OSBPL8-ZDHHC17 region (chr12) for HDL-C identified by one new SNP. The remaining 22 single-locus SNP effects confirmed previously reported genes or gene regions. For TC, three SNPs identified two gene regions that were tightly linked with previously reported genes associated with TC, including rs599839 that was 10 bases downstream PSRC1 and 3.498 kb downstream CELSR2, rs4970834 in CELSR2, and rs4245791 in ABCG8 that slightly overlapped with ABCG5. For HDL-C, LPL was confirmed by 12 SNPs 8-45 kb downstream, CETP by two SNPs 0.5-11 kb upstream, and the LIPG-ACAA2 region by five SNPs inside this region. Two epistasis effects on TC and thirteen epistasis effects on HDL-C reached the significance of "suggestive linkage". The most significant epistasis effect (p = 5.72 × 10-13) was close to reaching "significant linkage" and was a dominance × dominance effect of HDL-C between LMBRD1 (chr06) and the LRIG3 region (chr12), and this pair of gene regions had six other D × D effects with "suggestive linkage".


Genome-wide association analysis of the FHS data detected two new gene regions with genome-wide significance, detected epistatic SNP effects on TC and HDL-C with the significance of suggestive linkage in seven pairs of gene regions, and confirmed some previously reported gene regions associated with TC and HDL-C.

Peer Review reports


Total cholesterol (TC) is related to coronary diseases and high-density lipoprotein (HDL-C) cholesterol is anti-atherogenic. Genome-wide association studies (GWAS) and human genetic studies have identified a number of genes and gene regions affecting cholesterol phenotypes including TC and HDL-C [111]. A meta-analysis of HDL-C levels that include the FHS data has previously been published [2]. An early report on FHS [12] analyzed TC and HDL-C but used 100 k SNPs and a sample size that was much smaller than the current FHS sample size. Epistasis analysis of TC and HDL-C was unavailable. Here, we apply a quantitative genetics approach to detect additive or dominance single-locus effects and epistasis effects on log-transformed TC and HDL-C using 432,096 SNP markers and over 6000 individuals in FHS. The epistasis effects we tested included additive × additive (A × A), additive × dominance (A × D) or dominance × additive (D × A), and dominance × dominance (D × D) effects, with genetic interpretations of allele × allele, allele × genotype or genotype × allele, and genotype × genotype interactions. The single-locus analysis was intended to detect new targets or confirm existing targets using a method of analysis different from those used in previous reports based on an extended Kempthorne model that allows Hardy-Weinberger disequilibrium and linkage disequilibrium [13] for GWAS analysis of the FHS data while the epistasis analysis of TC and HDL-C was the first such attempt using the FHS data and the 500 k SNP panel.


The single-locus tests detected nine SNPs with additive (or allelic) effects on TC and twenty SNPs with additive effects on HDL-C that reached genome-wide significance (Tables 1-2). No dominance effect reached genome-wide significance. Among the twenty nine SNP effects, twenty were new effects that were not reported in previous studies and nine were previously reported to be associated with various cholesterol phenotypes [112]. Seven SNPs identified two new gene regions while the remaining twenty two SNPs confirmed previously reported gene regions. Two epistasis effects on TC and thirteen epistasis effects on HDL-C representing seven pairs of gene regions reached the significance of "suggestive linkage".

Table 1 Single-locus SNP effects for TC with genome control (GC) adjusted P < 7.2 × 10-8.
Table 2 Single-locus SNP effects for HDL-C with genome control (GC) adjusted P < 7.2 × 10-8.

Single-locus effects

For TC, nine SNPs with additive (or allelic) effects reached genome-wide significance with p < 7.2 × 10-8 (Table 1). Six SNPs inside or near four genes identified a new chr02 region containing RAB3GAP1, R3HDM1, LCT and MCM6 to be associated with TC (Figure 1A). Of the six SNPs in the RAB3GAP1-R3HDM1-LCT-MCM6 region, five SNPs were inside genes and one SNP was 4.2 kb upstream MCM6. The most significant SNP in this region was rs2322660 in intron 12 of LCT (Table 1). The RAB3GAP1-R3HDM1-LCT-MCM6 region contained two other genes (ZRANB3 and UBXD2) that did not have significant SNPs. Eleven other SNPs spanning a 1.23 Mb region (Figure 1A) that includes RAB3GAP1-R3HDM1-LCT-MCM6 had p-values between 1.27 × 10-5 and 7.13 × 10-7, including one SNP upstream ACMSD, one SNP in ACMSD, two SNPs in YSK4, one SNP in R3HDM1, one SNP in UBXD2 (also named UBXN4 according to NCBI [14]), two SNPs in LCT, and three SNPs downstream DARS (data not shown). These less significant results in the same neighborhood should add to the significance of the RAB3GAP1-R3HDM1-LCT-MCM6 region to TC. Three SNPs identified two genes that were tightly linked with previously reported genes associated with TC [1]. These three SNPs were rs599839 that was 10 bases downstream PSRC1 (chr01) and 3.498 kb downstream CELSR2, rs4970834 in intron 28 of CELSR2, and rs4245791 in intron 3 of ABCG8 (chr02) that slightly overlapped with ABCG5, where CELSR2 and ABCG5 regions were reported to be associated with TC in a recent GWAS report [1]. PSRC1 and ABCG8 also were reported to affect low-density lipoprotein cholesterol (LDL-C) [24, 6, 7, 9, 10]. The SNP (rs599839) that was 10 bases downstream PSRC1 had the most significant single-locus effect on TC (p = 3.7 × 10-16), while the SNP inside CELSR2 (rs4970834) had the second most significant single-locus effect on TC (p = 1.29 × 10-10).

Figure 1

Gene regions associated with total cholesterol (TC). A) A 1.23 Mb region containing RAB3GAP1-R3HDM1-LCT-MCM6 with multiple SNP effects on TC. B) One Mb region containing OSBPL8-ZDHHC17 associated with HDL-C. C) One Mb region containing LMBRD1 that had multiple SNPs interacting with an SNP near LRIG3 for TC. D) One Mb region containing LRIG3 which was near an SNP interacting with multiple SNPs in LMBRD1 for TC.

For HDL-C, twenty SNPs with additive effects reached genome-wide significance (Table 2). SNP rs17259942 identified a new gene region associated with HDL-C, the OSBPL8-ZDHHC17 region (q21.2, Figure 1B), with rs17259942 being 117 kb downstream OSBPL8 and 85 kb upstream ZDHHC17 [15]. According to NCBI [14], the OSBPL8-ZDHHC17 region contained three pseudo-genes (RPL7AP59, RPL21P98 and RPL7P43) and rs17259942 was 18 kb downstream RPL21P98 and 43 kb upstream RPL7P43. The other nineteen SNPs confirmed previously reported gene regions, including twelve SNPs 8-45 kb downstream LPL, two SNPs 0.5-11 kb upstream CETP, and five SNPs in the LIPG-ACAA2 region (39.812 kb downstream LIG and 65.963 kb upstream ACAA2) [14, 7, 11, 12]. LPL, CETP and LIPG were reported to be associated with HDL-C in four recent GWAS reports [14] while ACAA2 was reported in [3]. The SNP nearest to CETP (rs1800775) was the most significant effect (p = 8.61 × 10-34) in this study.

QQ plot for single SNP tests on TC and HDL-C showed that p-values of significant results all deviated from the expected p-values under the null hypothesis (Figure 2).

Figure 2

QQ plots for single-SNP whole genome association tests of total cholesterol (TC) and high-density lipoprotein cholesterol (HDL-C). A) TC. B) HDL-C.

Epistasis effects

Two epistasis effects on TC and thirteen epistasis effects on HDL-C reached the significance of suggestive linkage defined in [16] (Table 3). The two epistasis effects on TC involved two different pairs of gene regions while the thirteen epistasis effects on HDL-C involved five different pairs of gene regions, so that the fifteen epistasis effects identified seven pairs of gene regions. Eight SNPs in introns 1, 5, 7, 9, and 14 of LMBRD1 (chr06) (Figure 1C) interacted with a chr12 SNP about 53 kb from LRIG3 (q14.1, Figure 1D) and all these eight pairs had D × D effects on HDL-C. One of the eight epistasis effects involving intron 14 of LMBRD1 was the most significant epistasis effect that was close to reaching "significant linkage" defined in [16] or genome-wide significance with 5% Bonferroni corrected type-I error.

Table 3 Epistasis effects for TC and HDL-C with the significance of suggestive linkage.

Among the seven different pairs of gene regions with epistasis effects, four pairs had A × A effects, one pair had A × D effect, and two pairs had D × D effects (Table 3). For the A × A effect on TC involving chr04 and chr10, the A-T gamete had the highest TC value while the A-G gamete had the lowest TC value (Table 4). This showed that the G and T alleles of rs705169 on chr10 had significantly different effects when combined with the A allele of rs4437278 on chr04, noting that rs705169 did not have significant single-locus effect. The same phenomenon was also observed for the other three A × A effects in Table 4. For the A × D effect, the A-GG allele-genotype combination had the highest TC value while the G-GG allele-genotype combination had the lowest TC value. The two D × D effects were on HDL-C. For the D × D effect of rs4706271 × rs6581219 representing the eight pairs of D × D effects involving the same gene regions, GT-GG had the highest HDL-C value while GG-GG had the lowest HDL-C value. For the remaining D × D effect of rs12596869 × rs6506699 representing the two D × D effects of the same gene regions, CC-AG had the highest HDL-C value while CC-AA had the lowest HDL-C value (Table 4).

Table 4 Frequency and effect of gamete, allele-genotype or genotype-genotype combination in each epistasis effect with statistical significance of suggestive linkage.


The single-locus results in this study had strong confirmations with existing studies. For TC, we confirmed CELSR2 and ABCG5 reported in [1, 8]. These confirmed TC results should be considered as strong confirmation because our study had no overlapping samples with studies of [1, 8]. We detected seven effects on TC in the RAB3GAP1-R3HDM1-LCT-MCM6 region with the SNP in LCT being the most significant. Six of these seven effects had p-values for LDL-C in the range of 0.007-0.056 from a meta-analysis ([2], Table 1). This could be an indication about the significance on TC from a meta-analysis because LDL-C is calculated from TC [17]. A study in FINRISK cohorts with 14,140 individuals reported that LCT was associated with both TC and LDL-C with p-values in the range of 0.0005-0.005 [18]. In Silico replication using 1231 Italian subjects from the InCHIANTI cohort [19] generally lacked confirmation for the TC results in Table 1. The first three markers had p-values in the range of 0.005-0.07 while the other effects had p-values greater than 0.14 from the InCHIANTI cohort. The biological function of LCT for digesting lactose could be a reason for agreements and disagreements in replicating LCT effects on cholesterol. LCT affects lactose digestion and long-term consumption of lactose in rats was found to affect aortic cholesterol levels [20]. Therefore, dietary lactose levels that have not been considered by human GWAS could have affected the LCT results of different studies. MCM6 contains two of the regulatory regions for LCT [21] so that the significant effects in or near MCM6 (Table 1) could be due to MCM6's regulatory role to LCT. HDL-C had twenty significant SNP effects, but only one SNP identified a new gene region (OSBPL8-ZDHHC17) while all the other SNPs confirmed previously reported gene regions, although only seven of the twenty significant SNPs for HDL-C were reported previously (Table 2). OSBPL8 encodes a group of intracellular lipid receptors and suppresses ABCA1 [22], and ABCA1 was found to affect HDL-C level [23, 24]. For HDL-C, the InCHIANTI cohort did not confirm the effects in the LIPG-ACAA region (p > 0.55) but confirmed the other effects. In light of different samples and different methods of data analysis between our study and those in previous reports, the confirmations of gene results we observed for TC and HDL-C should be considered strong confirmations. This study used log-transformed TC and HDL-C values while recent GWAS on TC [1] and HDL-C [14] used the original observations of TC and HDL-C that deviated from normal distribution. However, single-locus effects from our study and previous studies [14] had remarkable mutual confirmation, indicating that single-locus analysis was somewhat robust to data distribution and possibly to methods of analysis.

Epistasis effects on TC and HDL-C were not reported in other GWAS so that a comparison between our epistasis results and those from others was unavailable. We detected eight SNP pairs indicating the interaction between gene LMBRD1 and gene LRIG3 with the significance of suggestive linkage. Both LMBRD1 and LRIG3 encode membrane proteins. LMBRD1 gene is involved in the transportation and metabolism of vitamin B12 which is important for metabolism of branched chain amino acids and odd chain fatty acids [25]. Replication using the InCHIANTI cohort did not confirm the epistasis results (p > 0.15).

The statistical power of epistasis testing is less than that for testing a single-locus effect, particularly for epistasis effects involving dominance such as A × D and D × D effects, with D × D effect being the most difficult to detect. The reason for this difficulty was due to the fact that higher-order effects explain less phenotypic variation even if the effect sizes were the same as lower-order effects [13]. The reduced power for epistasis testing could have contributed to the fact that the epistasis effects we detected only reached 'suggestive linkage' although the sample size was over 6000. The data analysis of this study showed that pairwise analysis was sensitive to outliers. This was due to the fact that artificially significant epistasis effects could occur when rare combinations of loci had extreme genotypic values by chance. This may happen when outliers exist due to the large number of pairwise effects arising from the large number of pairwise combinations. For example, over 466 billion pairwise effects (93,353,260,560 pairs × 5 effects per pair = 466,766,302,800 pairwise effects) were tested per trait in this study. A small fraction of random association between rare frequencies and outliers in opposite directions among a large number of pairs could yield a long list of artificially significant epistasis results. Therefore, dealing with outliers such as removing outliers and using data transformation is important in pairwise analysis. Pairwise analysis is computationally intensive but timely analysis is possible using parallel computing. Using 784 processor cores on the SGI Altix XE 1300 Linux cluster system with 2.66 GHz Intel Clovertown processor at the Minnesota Supercomputer Institute, the completion of pairwise epistasis analysis required about 15 hours per trait.


Genome-wide association analysis of the FHS data detected new single-locus and epistasis effects on TC and HDL-C and confirmed some previously reported effects. Additive effects were the primary single-locus effects of TC and HDL-C while epistasis effects involved allele × allele, allele × genotype (or genotype × allele), and genotype × genotype interactions.


Phenotype and SNP data

The FHS GWAS data (version 2) had 6575 individuals with SNP genotypes of the 500 k SNP panel from dbGAP. Of the 6575 individuals, 6431 had observations on TC and 6078 individuals had observations on HDL-C. A total of 496,858 SNP markers had known chromosome locations and 432,096 of these SNP markers with minor allele frequencies greater than or equal to 0.01 were analyzed.

Statistical Analysis

Original TC and HDL-C observations deviated from normality and had outliers (Figure 3A, D). The Box-Cox transformation analysis [26] implemented by the R statistical package [27] showed that the log-transformation was approximately the best transformation to achieve normality for those two traits (Figure 3B-C, E-F). One TC outlier, the highest TC value, was removed from the data analysis while no HDL-C outlier was removed. Log-transformed TC values were adjusted for blood sugar, body mass index, smoking status, and sex that had significant effects on log(TC). Age, age-squared, cholesterol treatment, and alcohol consumption were also tested for significant effects on log(TC) but were not included in the phenotypic model because they were insignificant. Log(HDL-C) was adjusted for age, age-squared, cholesterol treatment, blood sugar, body mass index, smoking status, number of cigars smoked, alcohol consumption and sex. Age was insignificant for HDL-C but was included because age-squared was nearly significant (p < 0.0543). Single-locus and epistasis effects for both traits were tested using the extended Kempthorne model that allows Hardy-Weinberg disequilibrium and linkage disequilibrium [13]. For each SNP, three effects were tested, genotypic, additive (A) and dominance (D) effects. For each SNP pair, five effects were tested, two-locus genotypic effect, A × A, A × D, D × A, and D × D epistasis effects. The EPISNPmpi parallel computing program [28] with a modification to implement a generalized least squares (GLS) analysis to account for sib correlations [29] was used to implement the statistical tests of single-locus and pairwise epistasis effects. For single-locus tests, p = 7.2 × 10-8 was used as the threshold p-value to declare genome-wide significance [30]. To assess genome-wide significance of pairwise epistasis results, we used 5% type-I error with the Bonferroni correction as the genome-wide significance. The 500 k SNP data was estimated to have 276,666 independent SNPs [31]. Each pairwise test was considered to have four independent tests although five effects were tested, because the two-locus marker genotypic effect was confounded with one of the four epistasis effects in reporting significant results. Therefore, the genome-wide 5% type-I errors with the Bonferroni correction was calculated as p = 0.05 [4(276,666)(276,665)/2]-1 = 3.266 × 10-13. This 5% significance level is equivalent to "significant linkage" defined in [16]. Since the Bonferroni correction is generally considered too severe, we also reported epistasis effects reaching "suggestive linkage" with statistical evidence that would be expected to occur one time at random in a genome-wide analysis [16]. In addition to the GLS method to account for sib correlation, the genomic control (GC) method [32] was used to account for potential sub-population structures in the three generation cohort of the FHS data set. For single-locus tests, all p-values were used to estimate inflation parameters for TC and HDL-C, yielding inflation parameter estimates of 1.14 and 1.11 respectively, and test statistics from the GLS tests were then adjusted by the estimates of inflation parameters and p-values were recalculated using the GC adjusted test statistics, which resulted in fewer significant effects. For the pairwise epistasis testing, we randomly selected 50,000 p-values and test statistics from over 466 billion pairwise tests for computational efficiency. Then we estimated the inflation parameters using two samples of 50,000 data points each for TC and HDL-C, yielding inflation parameter estimates of 1.01 and 1.05 respectively. All p-values were then adjusted using the inflation parameters and such adjustments also resulted in fewer significant epistasis results. Frequency of each subclass in an epistasis effect was calculated and each subclass was required to have a minimal number of five observations. After GC adjustment, QQ plots were made to show deviations of the observed p-values from the expected p-values under the null hypothesis for significant test results for single-locus tests only. QQ plot for epistasis effects were not made because the number of p-values for epistasis tests was too large. Gene locations of significant SNPs were identified according to ENSEMBL [15] and NCBI [14] based on Build 37.0 of the human genome.

Figure 3

Distributions of total cholesterol (TC) and high-density lipoprotein cholesterol (HDL-C) in original scales and in log-transformed scales. A) Distribution of TC in original scale deviated from normality and had an outlier to the far right. B) The Box-Cox maximum likelihood analysis showed that log-transformation (λ ≈ 0) was the best transformation to achieve normality for TC. C) Log-transformed TC values achieved normality. One outlier to the far right was removed from the data analysis. D) Distribution of HDL-C in original scale deviated from normality and had some outliers to the right. E) The Box-Cox maximum likelihood analysis showed that log-transformation (λ ≈ 0) was the best transformation to achieve normality for HDL-C. F) Log-transformed HDL-C values achieved normality without serious outliers.



genome-wide association study. SNP:single nucleotide polymorphism. TC:total cholesterol. HDL-C:high-density lipoprotein cholesterol. LDL-C:low-density lipoprotein cholesterol. LD:linkage disequilibrium. HWD:Hardy-Weinberg disequilibrium. A × A:additive × additive epistasis effect. A × D:additive × dominance epistasis effect. D × A:dominance × additive epistasis effect. D × D:dominance × dominance epistasis effect. GLS:generalized least squares.


  1. 1.

    Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, Pramstaller PP, Penninx BWJH, Janssens ACJW, Wilson JF, Spector T, Martin NG, Pedersen NL, Kyvik KO, Kaprio J, Hofman A, Freimer NB, Jarvelin MR, Gyllensten U, Campbell H, Rudan I, Johansson A, Marroni F, Hayward C, Vitart V, Jonasson I, Pattaro C, Wright A, Hastie N, Pichler I, Hicks AA, Falchi M, Willemsen G, Hottenga JJ, De Geus EJC, Montgomery GW, Whitfield J, Magnusson P, Saharinen J, Perola M, Silander K, Isaacs A, Sijbrands EJG, Uitterlinden AG, Witteman JCM, Oostra BA, Elliott P, Ruokonen A, Sabatti C, Gieger C, Meitinger T, Kronenberg F, Döring A, Wichmann HE, Smit JH, McCarthy MI, Duijn CM, Leena : Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet. 2008, 41: 47-55. 10.1038/ng.269.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, Schadt EE, Kaplan L, Bennett D, Li Y, Tanaka T, Voight BF, Bonnycastle LL, Jackson AU, Crawford G, Surti A, Guiducci C, Burtt NP, Parish S, Clarke R, Zelenika D, Kubalanza KA, Morken MA, Scott LJ, Stringham HM, Galan P, Swift AJ, Kuusisto J, Bergman RN, Sundvall J, Laakso M, Ferrucci L, Scheet P, Sanna S, Uda M, Yang Q, Lunetta KL, Dupuis J, De Bakker PIW, O'Donnell CJ, Chambers JC, Kooner JS, Hercberg S, Meneton P, Lakatta EG, Scuteri A, Schlessinger D, Tuomilehto J, Collins FS, Groop L, Altshuler D, Collins R, Lathrop GM, Melander O, Salomaa V, Peltonen L, Orho-Melander M, Ordovas JM, Boehnke M, Abecasis GR, Mohlke KL, Cupples LA: Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet. 2008, 41: 56-65. 10.1038/ng.291.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, Wahlstrand B, Hedner T, Corella D, Tai ES, Ordovas JM, Berglund G, Vartiainen E, Jousilahti P, Hedblad B, Taskinen M-R, Newton-Cheh C, Salomaa V, Peltonen L, Groop L, Altshuler DM, Orho-Melander M: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008, 40: 189-197. 10.1038/ng.75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM, Strait J, Duren WL, Maschio A, Busonero F, Mulas A, Albai G, Swift AJ, Morken MA, Narisu N, Bennett D, Parish S, Shen H, Galan P, Meneton P, Hercberg S, Zelenika D, Chen WM, Li Y, Scott LJ, Scheet PA, Sundvall J, Watanabe RM, Nagaraja R, Ebrahim S, Lawlor DA, Ben-Shlomo Y, Davey-Smith G, Shuldiner AR, Collins R, Bergman RN, Uda M, Tuomilehto J, Cao A, Collins FS, Lakatta E, Lathrop GM, Boehnke M, Schlessinger D, Mohlke KL, Abecasis GR: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008, 40: 161-169. 10.1038/ng.76.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Karvanen J, Silander K, Kee F, Tiret L, Salomaa V, Kuulasmaa K, Wiklund PG, Virtamo J, Saarela O, Perret C, Perola M, Peltonen L, Cambien F, Erdmann J, Samani NJ, Schunkert H, Evans A: The impact of newly identified loci on coronary heart disease, stroke and total mortality in the MORGAM prospective cohorts. Genet Epidemiol. 2009, 33: 237-246. 10.1002/gepi.20374.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Sandhu MS, Waterworth DM, Debenham SL, Wheeler E, Papadakis K, Zhao JH, Song K, Yuan X, Johnson T, Ashford S, Inouye M, Luben R, Sims M, Hadley D, McArdle W, Barter P, Kesäniemi YA, Mahley RW, McPherson R, Grundy SM, Wellcome Trust Case Control Consortium, Bingham SA, Khaw KT, Loos RJ, Waeber G, Barroso I, Strachan DP, Deloukas P, Vollenweider P, Wareham NJ, Mooser V: LDL-cholesterol concentrations: a genome-wide association study. Lancet. 2008, 371: 483-491. 10.1016/S0140-6736(08)60208-1.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Wallace C, Newhouse SJ, Braund P, Zhang F, Tobin M, Falchi M, Ahmadi K, Dobson RJ, Marcano AC, Hajat C, Burton P, Deloukas P, Brown M, Connell JM, Dominiczak A, Lathrop GM, Webster J, Farrall M, Spector T, Samani NJ, Caulfield MJ, Munroe PB: Genome-wide association study identifies genes for biomarkers of cardiovascular disease: Serum urate and dyslipidemia. Am J Hum Genet. 2008, 82: 139-149. 10.1016/j.ajhg.2007.11.001.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Samani NJ, Braund PS, Erdmann J, Gotz A, Tomaszewski M, Linsel-Nitschke P, Hajat C, Mangino M, Hengstenberg C, Stark K, Ziegler A, Caulfield M, Burton PR, Schunkert H, Tobin MD: The novel genetic variant predisposing to coronary artery disease in the region of the PSRC1 and CELSR2 genes on chromosome 1 associates with serum cholesterol. J Mol Med. 2008, 86: 1233-1241. 10.1007/s00109-008-0387-2.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Nakayama K, Bayasgalan T, Yamanaka K, Kumada M, Gotoh T, Utsumi N, Yanagisawa Y, Okayama M, Kajii E, Ishibashi S, Iwamoto S, The Jichi Community Genetics Team (JCOG): Large scale replication analysis of loci associated with lipid concentrations in a Japanese population. J Med Genet. 2009, 46: 370-374. 10.1136/jmg.2008.064063.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Schadt EE, Molony C, Chudin E, Hao K, Yang X, Lum PY, Kasarskis A, Zhang B, Wang S, Suver C, Zhu J, Millstein J, Sieberts S, Lamb J, Guha Thakurta D, Derry J, Storey JD, Avila-Campillo I, Kruger MJ, Johnson JM, Rohl CA, van Nas A, Mehrabian M, Drake TA, Lusis AJ, Smith RC, Guengerich FP, Strom SC, Schuetz E, Rushmore TH, Ulrich R: Mapping the genetic architecture of gene expression in human liver. PLoS Biol. 2008, 6: e107-10.1371/journal.pbio.0060107.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Kathiresan S, Manning AK, Demissie S, D'Agostino RB, Surti A, Guiducci C, Gianniny L, Burtt NP, Melander O, Orho-Melander M, Arnett DK, Peloso GM, Ordovas JM, Cupples LA: A genome-wide association study for blood lipid phenotypes in the Framingham Heart Study. BMC Med Genet. 2007, 19:8 (Suppl 1): S17-10.1186/1471-2350-8-S1-S17.

    Article  Google Scholar 

  12. 12.

    Murray A, Cluett C, Bandinelli S, Corsi AM, Ferrucci L, Guralnik J, Singleton A, Frayling T, Melzer D: Common lipid-altering gene variants are associated with therapeutic intervention thresholds of lipid levels in older people. Eur Heart J. 2009, 30: 1711-1719. 10.1093/eurheartj/ehp161.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Mao Y, London NR, Ma L, Dvorkin D, Da Y: Detection of SNP epistasis effects of quantitative traits using an extended Kempthorne model. Physiol Genomics. 2007, 28: 46-52.

    Article  Google Scholar 

  14. 14.

    National Center for Biotechnology Information: Last accessed 2/13/2010, []

  15. 15.

    ENSEMBL Genome Browser: Last accessed 2/13/2010, []

  16. 16.

    Lander ES, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995, 11: 241-247. 10.1038/ng1195-241.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Friedewald WT, Levy RI, Fredrickson DS: Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972, 18: 499-502.

    CAS  PubMed  Google Scholar 

  18. 18.

    Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, Auro K, Karvanen J, Kulathinal S, Niemelä M, Ellonen P, Vartiainen E, Jousilahti P, Saarela J, Kuulasmaa K, Evans A, Perola M, Salomaa V, Peltonen L: Gender Differences in Genetic Risk Profiles for Cardiovascular Disease. PLoS ONE. 2008, 3: e3615-10.1371/journal.pone.0003615.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Tanaka T, Shen J, Abecasis GR, Kisialiou A, Ordovas JM, Guralnik JM, Andrew Singleton A, Bandinelli S, Cherubini A, Arnett D, Tsai MY, Ferrucci L: Genome-Wide Association Study of Plasma Polyunsaturated Fatty Acids in the InCHIANTI Study. PLoS Genet. 2009, 5: e1000338-10.1371/journal.pgen.1000338.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Wostmann BS, Bruckner-Kardoss E: The effect of long-term feeding of 10% dietary lactose on serum, liver and aortic cholesterol of the rat and the gerbil. J Nutr. 1980, 110: 82-89.

    CAS  PubMed  Google Scholar 

  21. 21.

    Enattah NS, Sahi T, Savilahti E, Terwilliger JD, Leena Peltonen L, Järvelä I: Identification of a variant associated with adult-type hypolactasia. Nat Genet. 2002, 30: 233-237. 10.1038/ng826.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Yan D, Mäyränpää MI, Wong J, Perttilä J, Lehto M, Jauhiainen M, Kovanen PT, Ehnholm C, Brown AJ, Olkkonen VM: OSBP-related protein 8 (ORP8) suppresses ABCA1 expression and cholesterol efflux from macrophages. J Biol Chem. 2007, 283: 332-340. 10.1074/jbc.M705313200.

    Article  PubMed  Google Scholar 

  23. 23.

    Wang J, Burnett JR, Near S, Young K, Zinman B, Hanley AJ, Connelly PW, Harris SB, Hegele RA: Common and rare ABCA1 variants affecting plasma HDL cholesterol. Arterioscler Thromb Vasc Biol. 2000, 20: 1983-1989.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Cohen JC, Kiss RS, Pertsemlidis A, Marcel YL, McPherson R, Hobbs HH: Multiple rare alleles contribute to low plasma levels of HDL cholesterol. Science. 2004, 305: 869-872. 10.1126/science.1099870.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Rutsch F, Gailus S, Miousse IR, Suormala T, Sagné C, Toliat MR, Nurnberg G, Wittkampf T, Buers I, Sharifi A, Stucki M, Becker C, Baumgartner M, Robenek H, Marquardt T, Hohne W, Gasnier B, Rosenblatt DS, Fowler B, Nurnberg P: Identification of a putative lysosomal cobalamin exporter altered in the cblF defect of vitamin B(12) metabolism. Nat Genet. 2009, 41: 234-239. 10.1038/ng.294.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Box GEP, Cox DR: An analysis of transformations (with discussion). J Roy Stat Soc B. 1964, 26: 211-252.

    Google Scholar 

  27. 27.

    R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, ISBN 3-900051-07-0, 2008. Last accessed 6/22/09, []

  28. 28.

    Ma L, Runesha HB, Dvorkin D, Garbe JR, Da Y: Parallel and serial computing tools for testing single-locus and epistatic SNP effects of quantitative traits in genome-wide association studies. BMC Bioinformatics. 9: 315-10.1186/1471-2105-9-315.

  29. 29.

    Ma L, Amos CI, Yang Da Y: Accounting for correlations among individuals for testing SNP single-locus and epistasis effects in Genome-wide association analysis (Abstract). Plant & Animal Genomes XVI Conference, January 12-16, San Diego, CA. 2008, Last accessed 6/22/09, []

    Google Scholar 

  30. 30.

    Dudbridge F, Gusnanto A: Estimation of significance thresholds for genomewide association scans. Genet Epidemiol. 2008, 32: 227-234. 10.1002/gepi.20297.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Moskvina V, Schmidt K: On multiple-testing correction in genome-wide association studies. Genet Epidemiol. 2008, 32: 567-573. 10.1002/gepi.20331.

    Article  PubMed  Google Scholar 

  32. 32.

    Devlin B, Roeder K: Genomic control for association studies. Biometrics. 1999, 55: 997-1004. 10.1111/j.0006-341X.1999.00997.x.

    CAS  Article  PubMed  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


This research was supported in part by project MN-16-043 of the Agricultural Experiment Station at the University of Minnesota and by a grant from the Minnesota Supercomputer Institute. The Framingham Heart Study is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with Boston University. This manuscript was not prepared in collaboration with investigators of the Framingham Heart Study and does not necessarily reflect the opinions or views of the Framingham Heart Study, Boston University, or the NHLBI. The InCHIANTI study baseline (1998-2000) was supported as a "targeted project" (ICS110.1/RF97.71) by the Italian Ministry of Health and in part by the U.S. National Institute on Aging (Contracts: 263 MD 9164 and 263 MD 821336).

Author information



Corresponding author

Correspondence to Yang Da.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

LM processed the FHS data, implemented the data analysis and prepared the manuscript with YD. JY contributed to the data analysis and manuscript preparation. HBR directed the parallel computing implementation of the GWAS analysis. TT, LF and SB conducted the In Silico replication using the InCHIANTI cohort. YD coordinated this research and prepared the manuscript with LM. All authors read and approved this manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Ma, L., Yang, J., Runesha, H.B. et al. Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham Heart Study data. BMC Med Genet 11, 55 (2010).

Download citation


  • Total Cholesterol
  • Framingham Heart Study
  • Generalize Little Square
  • High Total Cholesterol
  • Suggestive Linkage