| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Atherosclerosis and Lipoproteins |
From Human Genetics Center (K.L.E.K., E.B., J.E.H.), University of Texas Health Science Center, Houston, Tex; Department of Human Genetics (C.F.S., T.J.R.), University of Michigan, Ann Arbor, Mich; Institute of Molecular Medicine (E.B., M.F.), University of Texas Health Science Center, Houston, Tex; Rockefeller University (S.C.H.), New York, NY; Department of Molecular Biology and Genetics (A.C.), Cornell University, Ithaca, NY.
Correspondence to Dr Kathy L.E. Klos, University of Texas Health Science Center at Houston, School of Public Health, Human Genetics Center, P.O. Box 20186, Houston, TX 77225.E-mail kathy.klos{at}uth.tmc.edu
| Abstract |
|---|
|
|
|---|
Methods and Results An average of 5 single nucleotide polymorphisms (SNPs) were genotyped within each of 45 genomic regions (54 genes) in blacks (1131 females and 812 males) and whites (1102 females and 954 males) from the Coronary Artery Risk Development in Young Adults (CARDIA) study. SNPs and gene-based 3-SNP haplotypes were evaluated for their ability to predict variation in plasma apolipoproteins (apo) A-I and apoB, total cholesterol (TC), high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglycerides (TG). We identified 14 SNPs in 6 candidate gene regions that explained statistically significant variation in the same trait in both genders of at least one race and with evidence of consistent genotype mean trend across gender within race. Haplotype analyses identified 9 candidate gene regions that explained statistically significant variation in one or both races.
Conclusion Four gene regions, ABCA1, APOA1/C3/A4/A5, APOE/C1/C4/C2, and CETP, explained plasma lipoprotein variation most consistently across strata. Other gene regions that influence plasma lipid and apolipoprotein levels within race include CYP7A1, LPL, PPARA, SOAT1, and SREBF2.
To identify common gene variations with a consistent influence we evaluated 262 SNPs in 54 genes for association with 6 plasma lipid and apolipoprotein traits. We identified 14 SNPs and 9 gene regions with genotype effects.
Key Words: lipids genetics of cardiovascular disease lipid and lipoprotein metabolism
| Introduction |
|---|
|
|
|---|
Reverse cholesterol transport (RCT) is the process by which cholesterol is effluxed from peripheral tissues onto acceptor particles, primarily HDL, in the plasma for uptake by the liver. RCT involves numerous lipid transfer proteins, enzymes, apolipoproteins, and membrane-bound receptors.11 The genes encoding these proteins, as well as genes encoding proteins that regulate their transcription, are candidates for influencing variation in plasma levels of apoA-I, apoB, HDL-C, LDL-C, TC, and TG. Based on a model RCT pathway,12,13 we identified a set of 54 genes involved in RCT for evaluating the impact of genetic variation on variation in plasma lipid and lipoprotein levels.
Genetic association studies have been plagued by inconsistent results. Two sources of these inconsistencies are false-positives (ie, type I error) and context-dependency of genotype effects.1416 In this study, we sought to identify gene regions having nongender-specific genotype effects on plasma lipid and apolipoprotein levels within population-based samples of black and white young adults (aged 18 to 30).
| Methods |
|---|
|
|
|---|
24 and >24 years) and education (high school or less versus higher than high school) groups. Study participants were given 6 sequential examinations from the time of the study initiation (1985 to 1986). Results shown here pertain to data collected at the first examination. All participants gave written informed consent, including consent for genetic studies. The CARDIA study was approved by the Institutional Review Boards of the 4 participating field centers, and this ancillary study was approved by additional Institutional Review Boards. Venous blood was drawn after a 12-hour fast. TC and TG were measured enzymatically.19 HDL-C was determined by precipitation with dextran sulfate/magnesium chloride.20 LDL-C was calculated using the Friedewald equation, after excluding individuals whose plasma triglycerides exceeded 400 mg/dL.21 Body mass index was calculated as weight (kg) divided by height squared (m2). apoA-I and apoB were determined using radioimmunoassay.22
Single Nucleotide Polymorphisms Selection and Genotyping
Single nucleotide polymorphisms (SNPs) were identified from NCBI dbSNP and Celera Discovery System databases in 54 genes belonging to a pathway model of RCT.13 An average of 5 SNPs per gene were selected for genotyping by identifying SNPs spaced evenly across the gene with relative frequencies of the rare allele reported in the dbSNP database (http://www.ncbi.nlm.nih.gov/SNP/) to be
5% for whites or blacks. Of 272 SNPs selected by this method and genotyped as below, 262 (96.4%) were retained after general quality control measures were applied (missing data <10%, and >98.5% concordance in blind duplicate analysis). Average spacing between SNPs within a gene was 9.06 Kb (range, 41 bp to 90.64 Kb).
For the present study, genotypes were obtained on 3993 individuals (1132 black females, 807 black males, 1101 white females, and 953 white males). Genotyping was performed using polymerase chain reaction amplification of genomic DNA, a short extension reaction across the polymorphic site, and mass spectrometry to detect allele-specific mass differences. Allele detection and genotype calling was performed using a MassARRAY System from Sequenom (San Diego, Calif). The sequences of the polymerase chain reaction and extension primers are available from the authors on request.
Statistical Analysis
Allele frequencies were obtained by direct counting. Hardy-Weinberg equilibrium was evaluated using a chi-square goodness of fit test. Individuals were removed from the analyses if they had not fasted 8 hours or more before the examination. TG levels were ln transformed to reduce skewness. Individuals with extreme values (±3 SD) were also removed. All lipid and apolipoprotein values were adjusted before analysis by fitting a race-, field center- and gender-specific linear regression model containing age, age2, age3, and body mass index. The residuals from the regression model were added back to the race- and gender-specific grand mean to produce an adjusted phenotype value for each individual. For the haplotype analyses, values were further adjusted for gender effects by adding the residuals from a model containing only gender back to the race-specific grand mean.
Differences in quantitative trait means among genotypes were evaluated by a 1-way ANOVA.23 Homogeneity of variances among genotypes was evaluated using Levene test.24 Where variances were unequal, Welch modified F statistic was used in assessing significance.25 To partially control for type 1 error in the single-SNP tests, a probability value from the 1-way ANOVA
0.01 in both genders of a race, or P
0.001 in one gender with P
0.05 in the other, was required for further consideration of genotype mean trends. Homogeneity of genotype mean trends was evaluated by ANOVA using a model that included genotype, gender, and genotype-by-gender interaction terms. Homogeneity of genotype mean trends across genders was rejected if P
0.05 for the interaction term. Consistency has become the sin quo non for evaluating significance of genetic association studies in light of continuing difficulties with multiple comparisons. In this study, we defined consistency as statistically significant and directionally similar effects in both genders within a race.
Haplotype frequencies within race were inferred from unphased genotype data using the expectation maximization algorithm.26,27 Haplotype trend regression, as implemented in the Helixtree (Golden Helix, Bozeman, Mont) software package, was used to evaluate association between plasma lipid and apolipoprotein levels and the haplotype probabilities of individuals.28 In this method, a matrix, D, of the inferred conditional probabilities of haplotypes given genotypes for each individual in the sample is created. A regression equation, E(Y)=Dß, is then used to relate inferred haplotype probabilities to the dependent variable. In this equation, Y is the matrix of individual response values, and D is the matrix of inferred conditional probabilities of haplotypes given the genotype. The conditional probability of a pair of haplotypes (for example h2,h3) for the ith individual with genotype Gi is
|
|
| Results |
|---|
|
|
|---|
0.05) plasma levels of TG, whereas females had higher levels of apoA-I and HDL-C. Black females had higher plasma TC than males. White males had higher plasma apoB and LDL-C than females.
|
Characteristics of the genotyped SNPs are summarized in Table 2. Of the 262 SNPs, 16% were exonic (61% nonsynonymous). Genotype frequencies for 2 SNPs in blacks and 9 SNPs in whites were out of Hardy-Weinberg equilibrium (P
0.001). In all except 2 cases, this was because of fewer than expected heterozygotes. In no case did results for SNPs not in Hardy-Weinberg equilibrium fit our criteria for consistency of association, nor were they involved in defining a 3-SNP haplotype having a statistically significant phenotypic effect.
|
Seventy-four of 262 SNPs were associated (P
0.05) with quantitative traits in black females, 68 in black males, 61 in white females, and 75 in white males (7% of 6280 valid statistical tests). The same trait was associated in both genders (P
0.01 in both genders, or P
0.001 in one gender with P
0.05 in the other) with variation in 11 SNPs in blacks and 8 SNPs in whites. These SNPs were evaluated for evidence of homogeneity of mean trends across gender strata using a nominal probability value of
0.05 as the threshold of significance. In blacks, homogeneity of apoB, LDL-C, and TC mean trends across gender was rejected (P
0.05) for apolipoprotein E/C1/C2/C4 gene cluster (APOE/C1/C2/C4) SNP rs741780 and homogeneity of LDL-C mean trends was rejected for APOE/C1/C2/C4 SNP rs405509.
Tables 3 and 4
summarize the pattern of significant associations between the 6 lipid traits and the 14 SNPs with evidence of a consistent effect across genders. Genotype means for the 14 SNPs in Tables 3 and 4
are presented in online supplemental Table I (see http://atvb.ahajournals.org). For
1 SNPs in 5 gene regions, genotype mean trends were consistent across genders within race for at least one HDL-related plasma trait (Table 3). For
1 SNPs in 2 gene regions, genotype mean trends were consistent across genders within race for at least 1 LDL-related plasma trait (Table 4).
|
|
Results of the gene-based haplotype sliding window analyses are plotted in Figure 1 for the HDL-related traits, plasma apoA-I, HDL-C, and lnTG, and in Figure 2 for the LDL-related traits plasma apoB, LDL-C, and TC. These results are plotted as the log10 false discovery rate-adjusted probability values from the haplotype trend regression analyses. 3-SNP haplotypes are arranged on the x-axis by the location from the p-terminal of chromosome 1 to the q-terminal of chromosome 21. Haplotypes were not estimated across SNPs >150 000 bp apart. For ease of comparison, all gene regions are plotted on the same axis. The threshold for statistical significance, an false discovery rate-adjusted P
0.05, is indicated by a horizontal red line.
|
|
Nine gene regions were statistically significantly associated with at least 1 plasma lipid or apolipoprotein trait in the haplotype trend regression analyses. In both races, HDL-C was associated with ATP-binding cassette transporter A1 (ABCA1), apolipoprotein A-I/apolipoprotein C-III/apolipoprotein A-IV/apolipoprotein A-V (APOA1/C3/A4/A5) cluster, and cholesteryl ester transfer protein (CETP) gene-based haplotypes. Also in both races, plasma apoB, LDL-C, and TC were associated with APOE/C1/C4/C2 cluster haplotypes. Unique to blacks were plasma apoA-I associations with cytochrome p450, family7, subfamily A, polypeptide 1 (CYP7A1), plasma HDL-C with APOE/C1/C4/C2, lipoprotein lipase (LPL), and sterol o-acyltransferase 1 (SOAT1), and plasma LDL-C with peroxisome proliferator-activated receptor alpha (PPARA) haplotypes. Unique to whites were plasma lnTG associations with APOA1/C3/A4/A5 and a plasma apoB association with sterol regulatory element binding transcription factor 2 (SREBF2) haplotypes.
| Discussion |
|---|
|
|
|---|
ABCA1 transports cellular cholesterol across the cell membrane to plasma acceptor particles such as apoA-I.30,31 In CARDIA blacks, HDL-C was consistently associated with ABCA1 SNP rs2515602, located at one peak of significance in the haplotype analyses. This region overlaps with the haplotype analysis peak location in whites. Together, significant haplotypes span approximately the 9th through the 18th exons of ABCA1, a region with high average LD in both races (average D'=0.90 in blacks and 0.99 in whites). SNPs rs2066718 (V771 mol/L) and rs2065412 were in the best haplotypes of both races. An independent effect of V771 mol/L on HDL-C has been previously reported,32 but not consistently.33 In CARDIA, genotype at V771 mol/L does not fully explain the HDL-C mean differences among individuals containing 2, 1, or 0 copies of common haplotypes (not shown).
apoA-I, apoC-III, apoA-IV, and apoA-V are components of plasma HDL and triglyceride-rich lipoproteins. In lipid metabolism, apoA-I serves as a cholesterol acceptor particle, apoC-III is a lipoprotein lipase inhibitor, and apoA-IV and apoA-V have been associated with variation in plasma TG, HDL-C, and TC levels.34 We previously reported the effect of APOA5 gene variation on plasma lnTG levels in CARDIA,35 and evaluation of additional SNPs in the region did not alter our conclusions. These analyses indicate that variation in the APOA1/C3/A4/A5 gene region influences HDL-C levels in CARDIA, consistent with previously reported associations of dyslipidemia36 and HDL-C.3740 The 3 haplotypes significant in the sliding window analysis spanned the APOA1 and APOA4 genes in whites, and the APOC3 gene in blacks. SNP coverage was sparse in this region (an average of just 4 SNPs per gene) and additional genotyping will be required to differentiate among genes in the region for effect on HDL-C.
Results from APOE/C1/C4/C2 gene region analyses in this study were in accord with the well-reported combined effect of APOE SNPs rs7412 and rs429358 (the
2/3/4 protein isoforms) on LDL-related plasma traits4144 and with APOE/C1/C4/C2 gene region effects on HDL-C.44.
Cholesteryl ester transfer protein is an important regulator of the movement of cholesteryl esters and triglycerides among lipoprotein particles.45 In these results, most prominent was the region of CETP containing the rs708272 (Taq1B) polymorphism, previously shown to be in strong LD with at least 1 promoter region polymorphism that may influence CETP concentration.46,47 In CARDIA, haplotype analyses were consistent with this finding (HDL-C means differed by copy number of common haplotypes containing the A allele of rs708272, P<0.05), but were also consistent with the presence of additional functional polymorphisms (HDL-C means also differed by copy number of some common haplotypes with the G allele at that locus P<0.05). Meta analysis has confirmed more than one CETP region polymorphism with an effect on HDL-C level consistent across samples.45
The protein encoded by CYP7A1 is the rate-limiting enzyme in bile acid synthesis.48 In blacks, consistent CYP7A1 associations with LDL-C and TC, characterized by increased levels in carriers of the rare alleles, were in accord with the previous reports of 5' region SNP effects.4951 Plasma apoA-I was associated with a 5' region CYP7A1 haplotype (rs1023649, rs8192871, rs11786580) in blacks (P=0.0001). Tests of equality of apoA-I means in individuals classed as having 0, 1, or 2 copies were significant (P<0.05) for haplotypes distinguished by an A allele at rs1023649. Although not consistently significant by the criteria applied in this study, association with apoA-I level was nominally significant (P<0.05) for rs8192871 in European-American males. Associations with plasma apoA-I and with HDL-C have been previously shown for CYP7A1 gene region polymorphisms in Pacific Islanders52 and in Asians,53 respectively. These findings support evidence of a CYP7A1 influence on HDL-related traits in addition to the influence on LDL-C and TC.
Lipoprotein lipase hydrolyzes triglycerides in plasma chylomicron and very-low-density lipoprotein particles.54 LPL gene variation has been reported to influence both plasma HDL-C and TG levels.5556 In CARDIA blacks, HDL-C was associated with 5 consecutive 3-SNP combinations spanning the first intron through the ninth codon (involving 6 of 8 measured SNPs). The nominal probability values from the haplotype trend regression analyses of these 5 haplotypes were uniformly lower (nominal P=0.0005 to P=0.0008) than probability values for single-SNP tests in the gender-stratified analyses (nominal P=0.0043 to P=0.8355 in females, and P=0.0086 to P=0.9043 in males). Although not consistently significant, associations between HDL-C and LPL variation were also observed in whites (hCV9642874 P=0.034 in females and P=0.008 in males; rs13266204 P=0.033 in females and P=0.015 in males). SNPs nominally significant in whites are in the first intron, whereas SNPs consistently significant in blacks are nearly 25 kb distal within the eighth intron and the 3' untranslated region. Sequence analysis of the LPL gene in blacks, whites, and Finns indicates a recombination hotspot between these 2 regions.57,58 Multiple LPL variable sites may function to influence HDL-C variation.
The peroxisome proliferator-activated receptor PPARA regulates fatty acid catabolism in the liver and other tissues.59 LDL-C was associated with a PPARA haplotype in blacks, but this "haplotype" consisted of a single 3' SNP, rs1053332, for which LDL-C means increased with copy number of the rare allele (P=0.0013). Coding region polymorphisms in PPARA have been associated with plasma levels of apoB, LDL-C, and TC in the Framingham Offspring Study,60 but those polymorphisms were not measured in this study. LD (measured as D') between the 3 coding region SNPs measured and rs1053332 averaged 0.52 in blacks. Further investigation will be required to assess whether rs1053332 associations were attributable to LD with previously studied polymorphic sites.
SOAT1 functions intracellularly to convert cholesterol into cholesteryl esters.61,62 HDL-C was associated in blacks with a haplotype of SOAT1 consisting of rs2783391, rs2247071, and rs2493121, a promoter region and 2 intronic polymorphisms, respectively. A "C" allele at both rs2783391 and rs2247071 distinguished haplotypes with significantly different copy number class means from other common haplotypes. Although not consistently significant, HDL-C was associated with rs2862616 in black males (P=0.0294), whereas apoA-I level was associated with rs2783391, rs2247071, and rs2255375 in white males (P=0.0008, P=0.0281, and P=0.0022, respectively). Although LD across the region was not high (average D'=0.68 in blacks and 0.78 in whites), it remains most likely that these associations are attributable to LD with
1 unmeasured polymorphic sites that function to alter gene expression or structure. Outside of this study, little work has been done to characterize the relationship between SOAT1 polymorphism and plasma lipid metabolism in humans.
The SREBF2 transcription factor regulates genes involved in cholesterol homeostasis.63 In this study, an association between plasma apoB level and a single SREBF2 haplotype (rs1009544, rs1883205, and rs2284083), entirely in the first intron, was unique to whites. Associations with TC and LDL-C have been previously reported for SREBF2 region polymorphisms,64,65 and these finding are further support of a role for this gene in LDL metabolism.
Although the sparse SNP coverage of many of the genes in this study precludes the ability to rule out any gene as having an effect on plasma lipid metabolism, these results suggest the presence of influential variations in at least 9 gene regions. Despite the complexity of genetic and environmental control of lipid metabolism, some genes are likely to influence these CAD risk factors in a manner that is consistent across genders. A more comprehensive evaluation of potential functional variations in and around these genes is warranted.
| Acknowledgments |
|---|
This work was supported by NIH grants HL072905 HL072810, HL072904, and GM065509.
Disclosures
None.
| Footnotes |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Klos, L. Shimmin, C. Ballantyne, E. Boerwinkle, A. Clark, J. Coresh, C. Hanis, K. Liu, S. Sayre, and J. Hixson APOE/C1/C4/C2 hepatic control region polymorphism influences plasma apoE and LDL cholesterol levels Hum. Mol. Genet., July 1, 2008; 17(13): 2039 - 2046. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Thompson, E. Di Angelantonio, N. Sarwar, S. Erqou, D. Saleheen, R. P. F. Dullaart, B. Keavney, Z. Ye, and J. Danesh Association of Cholesteryl Ester Transfer Protein Genotypes With CETP Mass and Activity, Lipid Levels, and Coronary Risk JAMA, June 18, 2008; 299(23): 2777 - 2788. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. R. Graham, C. Kyogoku, S. Sigurdsson, I. A. Vlasova, L. R. L. Davies, E. C. Baechler, R. M. Plenge, T. Koeuth, W. A. Ortmann, G. Hom, et al. Three functional variants of IFN regulatory factor 5 (IRF5) define risk and protective haplotypes for human lupus PNAS, April 17, 2007; 104(16): 6758 - 6763. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
ATVB Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2006 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |