Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comprehensive Evaluation of the Association of APOE Genetic Variation with Plasma Lipoprotein Traits in U.S. Whites and African Blacks

  • Zaheda H. Radwan,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Xingbin Wang,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Fahad Waqar,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Dilek Pirim,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • Vipavee Niemsiri,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • John E. Hokanson,

    Affiliation Department of Epidemiology, Colorado School of Public Health, University of Colorado Denver, Aurora, Colorado, United States of America

  • Richard F. Hamman,

    Affiliation Department of Epidemiology, Colorado School of Public Health, University of Colorado Denver, Aurora, Colorado, United States of America

  • Clareann H. Bunker,

    Affiliation Department of Epidemiology, Graduate School of Public Health, University Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • M. Michael Barmada,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • F. Yesim Demirci,

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

  • M. Ilyas Kamboh

    kamboh@pitt.edu

    Affiliation Department of Human Genetics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America

Abstract

Although common APOE genetic variation has a major influence on plasma LDL-cholesterol, its role in affecting HDL-cholesterol and triglycerides is not well established. Recent genome-wide association studies suggest that APOE also affects plasma variation in HDL-cholesterol and triglycerides. It is thus important to resequence the APOE gene to identify both common and uncommon variants that affect plasma lipid profile. Here, we have sequenced the APOE gene in 190 subjects with extreme HDL-cholesterol levels selected from two well-defined epidemiological samples of U.S. non-Hispanic Whites (NHWs) and African Blacks followed by genotyping of identified variants in the entire datasets (623 NHWs, 788 African Blacks) and association analyses with major lipid traits. We identified a total of 40 sequence variants, of which 10 are novel. A total of 32 variants, including common tagSNPs (≥5% frequency) and all uncommon variants (<5% frequency) were successfully genotyped and considered for genotype-phenotype associations. Other than the established associations of APOE*2 and APOE*4 with LDL-cholesterol, we have identified additional independent associations with LDL-cholesterol. We have also identified multiple associations of uncommon and common APOE variants with HDL-cholesterol and triglycerides. Our comprehensive sequencing and genotype-phenotype analyses indicate that APOE genetic variation impacts HDL-cholesterol and triglycerides in addition to affecting LDL-cholesterol.

Introduction

Coronary heart disease (CHD), a multifactorial disease modulated by multiple genetic and environmental factors, continues to be a leading cause of morbidity and mortality worldwide [1]. Dyslipidemia with high low-density lipoprotein cholesterol (LDL-C) and low high-density lipoprotein cholesterol (HDL-C) is associated with high risk of CHD [1]. Genes involved in lipid metabolism are considered to be candidate genes for CHD risk, and their genetic variation could contribute, in part, to the inter-individual variation in plasma lipoprotein-lipid levels.

Apolipoprotein E (ApoE, protein; APOE, gene) is a major constituent of very low-density lipoproteins (VLDL) and high-density lipoproteins (HDL) [2], [3] and plays a crucial role in lipid metabolism through enhancing hepatic uptake of triglyceride-rich lipoproteins (TGRL) and participating in reverse cholesterol transport mechanism (RCT) [4]. Besides its significant contribution in lipid metabolism, ApoE is involved in multiple functions in the human body, including nerve growth and regeneration [5][8], cognitive function [9], [10], immunoregulation and influencing susceptibility to infectious diseases [11][13]. The APOE gene is located on chromosome 19q13.32 as part of the APOE-C1-C4-C2 gene cluster, and is composed of 4 exons and 3 introns that span 3.6 kb [14] and encodes for 299 amino acids [3].

APOE is one of the most extensively studied candidate genes and the influence of its genetic variation on plasma lipid levels and CHD risk has been well investigated [15][16]. The epsilon polymorphism of APOE is defined by the rs7412 and rs429358 SNPs which leads to the generation of ApoE2, ApoE3 and Apo E4 isoforms and are coded by three codominant alleles (designated as E*2 E*3 and E*4). The three isoforms differ by an amino acid substitution at position 112 or position 158 in the 299-amino-acid peptide chain. Although the major effect of APOE genetic variation has been reported to be on LDL-C levels, recent genome-wide association studies (GWAS) on lipid traits also identified statistically significant associations of APOE common variants with HDL-C and triglyceride (TG) levels [17][18]. Thus, deep resequencing of the APOE gene in selected individuals with high/low lipid levels is warranted in order to characterize both rare and common variants that might affect plasma lipid profile.

In this study, we resequenced the entire APOE gene region (total 5.5 kb), including all four exons (1,180 bp), three introns (2,432 bp), and ∼1 kb of each of the flanking regions in selected individuals with extreme HDL-C levels (falling within the upper and lower 10th percentiles) from two ethnically-distinct populations (95 US non-Hispanic Whites (NHWs) and 95 African Blacks). Following the sequencing-based discovery step, we genotyped all identified common tagSNPs (r2≥0.9) with minor allele frequency (MAF) ≥5%, and relevant uncommon and rare variants with MAF<5% in the entire sample sets (623 NHWs and 788 African Blacks) to evaluate their associations with lipid traits. The association of APOE genetic variation was examined with three lipid traits (LDL-C, HDL-C and TG) and apolipoprotein B (ApoB) using single-site association analysis for variants with MAF≥1%, gene–based and haplotype-based association analyses for all variants, and SKAT-O (sequencing Kernel association optimal test) for uncommon and rare variants (MAF<5%).

Materials and Methods

Study Samples

The study was conducted on two epidemiologically well-characterized population samples comprising 623 US non-Hispanic Whites (NHWs) and 788 African Blacks. NHW samples were collected as part of the San Luis valley Diabetes Study that was designed as geographical case-control study of non-insulin dependent diabetes mellitus and cardiovascular disease in Alamosa and Conejos counties of South Colorado [19). All NHWs used in this study were non-diabetic controls and the basic characteristics of this study are described elsewhere [19][20]. African Blacks were recruited from Benin City, Nigeria as part of a study on CHD risk factors in Blacks and the study details have been described in Bunker et al. [21][22]. While LDL-C, HDL-C and TG were measured in all subjects, ApoB was measured only in a subset of NHW individuals [23][24]. The demographic and lipid characteristics of these study samples can be found in our previous publications [24][26]. The study was approved by the University of Pittsburgh and University of Colorado Denver Institutional Review Boards and all study participants provided written informed consent.

DNA Extraction

The genomic DNA used for sequencing and genotyping was extracted from blood clots in Blacks and from buffy coats in NHWs using standard procedures.

DNA Sequencing

Ninety-five individuals with high HDL-C levels falling within the upper 10th percentile (47 NHWs, and 48 African Blacks) and 95 individuals with low HDL-C levels falling in the lower 10th percentile (48 NHWs, and 47 African Blacks) were selected for Sanger sequencing. The characteristics of the selected samples in both ethnic groups are summarized in Table S1 in S1 File.

A total of ∼5.5 kb of the APOE gene region, including all 4 exons and 3 introns, 1,034 bp in 5′ flanking region, and 845 bp in 3′ flanking region were PCR-amplified using M13 tagged forward and reverse primers. Publicly available information at SeattleSNPs database (http://pga.mbt.washington.edu/) was used to order M13 tagged primers, which generated nine overlapping PCR amplicons. PCR reaction and cycling conditions are available upon request. The PCR-amplified samples were sent to a commercial lab (Beckman Coulter Genomics, Danvers, MA) for automated fluorescence-based cycle sequencing and capillary electrophoresis on ABI 3730x1DNA Analyzers. Variant Reporter version 1.0 (Applied Biosystems, Foster City, CA) and Sequencher version 4.8 (Gene Codes Corporation, Ann Arbor, MI) were used for sequencing analysis and variant detection.

DNA Genotyping

Common tagSNPs (MAF≥5%) were determined by Tagger analysis of the sequencing data in each ethnic group using Haploview software and an r2 cut-off of 0.9. All common tagSNPs and uncommon/rare variants (MAF<5%) identified in each ethnic group by our sequencing, as well as the suspicious variants with low sequencing quality and/or low coverage that warrant validation and the previously reported common variants not detected in our sequencing, were selected for follow-up genotyping.

TaqMan (Applied Biosystems) or iPLEX Gold (Sequenom, San Diego, CA) genotyping methods were used for genotyping following manufacturer's protocols and recommendations. Whole genome amplified DNAs dried in 384-well plates were used for genotyping. Endpoint fluorescence reading of custom or pre-made TaqMan assays was done using the ABI Prism 7900HT Sequence Detection System. The iPLEX Gold genotyping was performed in the Core laboratories of the University of Pittsburgh. Sequences of primers and probes used for genotyping are available upon request. All the samples used in sequencing were also included in genotyping as a quality control measure. The comparison of sequencing and genotyping calls was conducted to check the concordance as well as to increase the call rate in both sequencing and genotyping sets.

Statistical Analysis

Analyses for NHWs and African Blacks were performed separately. For sequencing subsets, the Haploview software (www.broadinstitute.org/haploview) was used to analyze allele frequencies, their distributions in the two extreme HDL-C groups, their concordance with Hardy-Weinberg Equilibrium (HWE), and their linkage disequilibrium (LD) patterns.

SNPs with extensive missing data (>20%) and/or deviating highly from HWE (P<0.01) were excluded from association analyses. A total of 15 variants in NHWs and 23 variants in Blacks remained for downstream analysis. The associations between SNPs and lipid traits were analyzed using additive linear regression model. We took the best power Box-Cox transformation such that the transformed lipid traits achieved normality. Stepwise regression in both directions was performed to identify significant covariates for each lipid trait. The covariates included were gender, age, BMI and smoking in NHWs and gender, age, BMI, waist measurement, smoking, exercise (minutes walking or bicycling to work each day), and staff level (junior or senior) in Blacks. Detailed information on those covariates and their effects can be found elsewhere [24]. Since the epsilon APOE E2/E3/E4 polymorphism has an established effect on cholesterol levels, we also adjusted the effects of novel associations for the epsilon APOE polymorphism. Single-site, haplotype-based and rare variants analyses were implemented in R and the versatile gene-based associations (VEGAS) [27] were also performed. For single-site analysis, we applied Benjamini-Hochberg procedure [28] to control for false discovery rate (FDR) and considered an FDR (q-value) of <0.20 as statistically significant.

For haplotype association analysis, the generalized linear model (GLM) was used [29]. Including too many haplotypes can make above model inefficient and impractical. To reduce the number of haplotypes considered in association analysis, we used the sliding window, 4 SNPs per window, and assessed evidence for association within each window. Specifically, a global p-value for testing overall effects of the haplotypes with frequency greater than 0.01 was used to assess the associations between the traits and haplotypes in each window. Sliding-window haplotype analysis was performed using the haplo.glm function in the Haplo.Stats R package (version 1.5.0).

We analyzed the cumulative effects of uncommon/rare variants by using the SKAT-O method [30], which has been proposed to be the optimal test for rare variant analysis and outperformed the SKAT and burden tests in several ways. The analysis was performed by using three different minor allele frequency bin thresholds (≤1%, ≤2% and <5%). The SKAT method was implemented using the “SKAT” R package.

Results

APOE Sequencing Results

Sequencing of ∼5.5 kb genomic region of APOE (including all 4 exons, 3 introns, 1,034 bp in the 5′ and 845 bp in the 3′ flanking regions), in 190 selected individuals (95 NHWs and 95 African Blacks) with extreme HDL-C levels revealed a total of 40 variants in both population groups, including 30 known and 10 novel variants (as compared to NCBI dbSNP human Build 141) (Table S2 in S1 File). All novel variants identified in this study have been submitted to dbSNP database: (http://www.ncbi.nlm.nih.gov/SNP/snp_viewTable.cgi?handle=KAMBOH).

The codon position used for specifying the coding variants corresponds to the premature protein that also includes the first 18 amino acids of signal peptide. The distribution of the 40 variants is as follows: 10 in 5′ flanking region, 7 in exons (including 2 in 3′ UTR), 16 in introns (including 1 in splice site), and 7 in the 3′ flanking region. Four of the 5 coding variants (80%) were non-synonymous. Ten of the 40 variants were present in both groups, while 9 variants were unique to NHWs and 21 variants were specific to African Blacks. Four of the ten shared-variants showed statistically significant allele frequency differences between the two ethnic groups (see Table S2 in S1 File for variants at positions 560, 624, 832, and 1163).

Distribution of APOE sequence variants in two extreme HDL-C groups

Comparison of sequencing variants distribution between the two extreme HDL-C groups in NHWs and African Blacks is presented in Table S3 in S1 File and Table S4 in S1 File, respectively.

Among the 8 rare/uncommon variants (overall MAF<5%) in NHWs, 6 were unique to the high HDL-C group, 1 was unique to the low HDL-C group, and 1 was present in both lipid groups. In parallel with observing more unique rare variants in the high HDL-C group, 21% (10/47 subjects) of this group had at least one unique rare variant as compared to 2% (1/48 subjects) of the other lipid group (Fisher exact test p-value = 0.0037). Furthermore, the two rare coding variants observed in this study (Ala23Ala; Val254Glu) were present only in the high HDL-C group.

Among the 21 rare/uncommon variants (overall MAF<5%) observed in African Blacks, 6 were unique to the high HDL-C group, 5 were unique to the low HDL-C group, and the remaining 10 were equally distributed among the two extreme HDL-C groups. Unlike NHWs, the distribution of the unique rare variants was similar in the two extreme lipid groups among African Blacks. Fifteen percent (7/48 subjects) of the high HDL-C group had at least one unique rare variant as compared to 6% (3/47 subjects) of the low HDL-C group (Fisher exact test p-value = 0.316).

Single-site association analysis of the SNPs in the entire NHW and Black samples

Following the identification of genetic variation in the sequencing step, common tagSNPs covering the entire APOE gene and rare variants were genotyped in the total sample of NHWs (n = 623) and African Blacks (n = 788) for genotype-phenotype association analyses. Initially, 20 variants in NHWs (9 tagSNPs, 8 rare variants, 2 suspicious SNPs, and 1 database SNP) and 32 variants in African Blacks (9 tagSNPs, 21 rare variants, 1 suspicious SNP, and 1 database SNP) were selected for genotyping. In NHWs, 2 of the 20 variants (APOE2294; MAF = 0.005, and APOE4951/rs1081105;MAF = 0.042) failed in both Sequenom and TaqMan designs or runs, 2 suspicious variants (APOE4489, and APOE4490) were confirmed as not being genuine and one variant (APOE624/rs769446) with low call rate was excluded from the association analyses. So, a total of 15 variants (14 sequencing variants and 1 database SNP APOE3106/rs769452) were successfully genotyped in the entire NHW sample. In African Blacks, 6 of 32 variants (APOE471/rs439382;MAF = 0.132, APOE494;MAF = 0.005, APOE526;MAF = 0.005, APOE2576;MAF = 0.005, APOE4951/rs1081105; MAF = 0.042, and APOE5229/rs80125357; MAF = 0.059) failed in both Sequenom and TaqMan designs or runs, the database SNP (APOE1586/rs74625294) and the suspicious variant (APOE91) were excluded because they turned out to be non-polymorphic in our population and an additional variant (APOE1591/rs147236548) was excluded from the statistical analyses because it was out of HWE. Thus, a total of 23 variants were successfully genotyped in the entire African Black sample.

The LD plot of the genotyped variants with MAF>1% in NHWs is shown in Fig. 1, the association results for all genotyped variants with the three lipid traits (LDL-C, TG, and HDL-C) and ApoB are presented in Table 1 and the adjusted mean distributions of all the evaluated lipid traits among the genotype groups are given in Table S5 in S1 File. As expected, the two known and well-established SNPs as part of the APOE epsilon polymorphism, E*4 (rs429358) and E*2 (rs7412) were significantly associated with plasma levels of LDL-C (β = 8.10; p = 0.0103, and β = −21.84; p = 1.84E-07, respectively) and ApoB (β = 2.14; p = 0.0005, and β = −5.60; p = 9.65E-13, respectively). Four additional LDL-C associations were observed independent of E*2/E*4: APOE832/rs405509 in 5′ flanking (β = −5.17; p = 0.0345; FDR = 0.139), APOE1163/rs440446 in intron 1 (β = 6.11; p = 0.018; FDR = 0.139), APOE2440/rs769450 in intron2 (β = 5.52; p = 0.0275; FDR = 0.139), and APOE4310/rs199768005 (Val254Glu) in exon4 (β = −35.36; p = 0.043; FDR = 0.139). These same four SNPs were also associated with TG (p = 0.0019 and FDR = 0.01, p = 0.0012 and FDR = 0.01, p = 0.002 and FDR = 0.01, and p = 0.028 and FDR = 0.074, respectively). An additional SNP, APOE4528/rs374329439 in 3′UTR, was also associated with TG (p = 0.022; FDR = 0.071).

thumbnail
Figure 1. LD Plot of the genotyped variants with MAF>1% in NHWs.

The values in the cells are the pairwise degree of LD indicated by r2×100. r2 = 0 is shown as white, 0<r2<1 is shown in gray and r2 = 1 is shown in black.

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

thumbnail
Table 1. Single-site association analysis with lipid traits in NHWs.

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

In African Blacks, 23 variants with high call rate and in compliance with HWE were included in the association analyses and their single-site association results are shown in Table 2. The LD plot of the genotyped variants with MAF>1% is shown in Fig. 2 and the adjusted mean of the evaluated lipid traits among the genotype groups are presented in Table S6 in S1 File. As expected, the E*4 (rs429358) and E*2 (rs7412) SNPs were associated with LDL-C (β = 0.46; p = 0.0317 and β = −2.05; p = 5.35E-07, respectively). Four additional variants also showed association with LDL-C independent of E*2/E*4: APOE2269/rs61357706 in intron 2 (β = −2.23; p = 0.0034; FDR = 0.02), APOE2544/rs115299243 in intron 2 (β = −2.54; p = 0.0008; FDR = 0.008), APOE4036/rs769455(Arg163Cys) in exon 4 (β = −2.41; p = 0.0004; FDR = 0.008) and a novel association in 3′UTR, APOE4569 (β = 8.35; p = 0.024; FDR = 0.124). Two of these variants were also associated with TG APOE4036/rs769455 (Arg163Cys) (p = 0.0343; FDR = 0.199) and APOE2544/rs115299243 (p = 0.0378; FDR = 0.199). Two additional variants were also found to be associated with TG: APOE73/rs1081101 (p = 0.0115; FDR = 0.145), and APOE1279/rs877973 (p = 0.014; FDR = 0.145). One novel rare variant (APOE618) located in 5′ flanking region and observed in one individual was associated with extremely low HDL-C (13.5 mg/dl vs. 47.8 mg/dl; (β = −12.18; p =  0.001; FDR = 0.020); Table S6 in S1 File).

thumbnail
Figure 2. LD Plot of the genotyped variants with MAF>1% in African Blacks.

The values in the cells are the pairwise degree of LD indicated by r2×100. r2 = 0 is shown as white, 0<r2<1 is shown in gray and r2 = 1 is shown in black.

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

thumbnail
Table 2. Single-site association analysis with lipid traits in African Blacks.

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

Gene-based association analysis

Gene-based tests including all APOE common and rare variants simultaneously within each ethnic group were performed (Table 3 and Table 4). Gene-based association analysis showed significant associations (p<0.05) with TG, LDL-C and ApoB in NHWs and with LDL-C in African Blacks.

thumbnail
Table 3. Gene-based association analysis with lipid traits in NHWs.

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

thumbnail
Table 4. Gene-based association analysis with lipid traits in African Blacks.

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

Haplotype-based association analysis

The adjacent variants were evaluated as a group of four variants instead of relying on the effect of a single variant. The p-values for 4-SNP haplotype windows for each evaluated lipid trait are given in Fig. 3, and Fig. 4 for NHWs and Blacks, respectively. For the haplotype-based association results please see the Tables S7-S12 in S1 File for NHWs and Tables S13-S16 in S1 File for Blacks.

thumbnail
Figure 3. Haplotype analysis with lipid traits in NHWs.

Haplotype windows for LDL-C (a), for ApoB (b), for HDL-C (c), and for TG (d). X-axis has the genotyped markers names and the Y-axis has the –log (global p-value), horizontal lines represent the 4-SNP windows, red-line represents the p-value threshold (p = 0.05) and everything below the threshold is considered non-significant and vice versa.

https://doi.org/10.1371/journal.pone.0114618.g003

thumbnail
Figure 4. Haplotype analysis with lipid traits in African Blacks.

Haplotype windows for LDL-C (a), for ApoB (b), for HDL-C (c), and for TG (d). X-axis has the genotyped markers name and the Y-axis has the –log (global p-value), horizontal lines represent the 4-SNP windows, red-line represents the p-value threshold (p = 0.05) and everything below the threshold is considered non-significant and vice versa.

https://doi.org/10.1371/journal.pone.0114618.g004

In NHWs, the strongest haplotype associations were observed with ApoB followed by LDL-C. The region covered by five consecutive haplotype windows, including windows 7, 8, 9, 10 and 11 that harbor the variation in exon 4, showed the most significant global p-value with LDL-C (p-value ranges between 1.12E-07 and 0.0339), most likely due to the effect of E*2 (rs7412) and E*4 (rs429358) SNPs present in these windows. Additional four windows (1, 3, 4 and 5) showed nominally significant global p-value with LDL-C confirming the independent effect of APOE2440/rs769450 (p = 0.038). Similarly, the consecutive windows 7, 8, 9, and 10 that harbor variation in exon 4 showed significant haplotype global p-values with ApoB (p = 0.0027, 4.37E-14, 8.32E-13, and 5.47E-12) more likely due to the significant contribution of E*2 (rs7412) and E*4 (rs429358) on ApoB variation. Additionally the first five windows (windows 1, 2, 3, 4, and 5) showed significant global p-values (p = 1.57E-05, 0.0004, 8.05E-07, 0.0.0265, and 0.0176) more likely due to the effects of APOE832/rs405509 and APOE1998/rs769449 variants on ApoB Moreover, four windows (windows 1, 2, 3, and 5) showed independent evidence of association with TG (p = 0.0043, 0.0196, 0.0194, and 0.0344) likely to be mediated by the following three variants; APOE832/rs405509 (p = 0.003), APOE1163/rs440446 (p = 0.0018), and APOE2440/rs769450 (p = 0.0082) confirming their single-site effects. Only the last window, window 12 showed significant haplotype association with HDL-C (p = 0.0301), more likely due to the effect APOE4737/rs117656888.

In Blacks, the strongest haplotype associations were observed with LDL-C. Similar to NHWs, the last four windows (17, 18, 19, and 20), which include common polymorphisms in exon 4 showed the most significant p-values with LDL-C in African Blacks (p-values range between 8.86E-09 and 8.2E-06). Additional twelve windows showed significant effect on LDL-C (p-values range between 0.0022 and 0.036) including windows 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, and 16 and confirming the single-site effects of multiple variants (APOE560rs449647, APOE624/rs769446, APOE832/rs405509, APOE2269/rs61357706, and APOE2544/rs115299243) on LDL-C Unlike NHWs, only two windows (18, and 19) showed significant global p-value (0.035, and 0.038) with ApoB more likely due to the significant effect of E*2. Only the first window showed significant global p-value with TG (p = 0.035), most likely due to the effect of APOE73/rs1081101 as seen in the single-site analysis (p = 0.0093). Findings from haplotype-based association analyses confirm the single-site association results.

Uncommon/Rare variants association analysis

Uncommon/rare variants association analysis was performed to examine the cumulative effect of uncommon/rare variants (MAF<5%) on lipid traits (HDL-C, LDL-C, and TG) using SKAT-O test. We found significant association with HDL-C in NHWs after including all 7 uncommon/rare variants in the analysis (p = 0.0061), and APOE1575/rs769448 with MAF = 0.021 contributed largely to this significance (Table 5) as it also showed the most significant association in single-site analysis (p = 0.0197). In Blacks, rare variants analysis (Table 6) showed significant association with LDL-C (p = 0.00018) and the significant association was driven by three variants with MAF between 0.017 and 0.020 (APOE2269/rs61357706, APOE2544/rs115299243, and APOE4036/rs769455), all of which showed significant association in single-site analysis (p range = 0.0009–0.0064).

thumbnail
Table 5. Rare/uncommon (MAF<5%) variants analysis with lipid traits in NHWs.

https://doi.org/10.1371/journal.pone.0114618.t005

thumbnail
Table 6. Rare/uncommon (MAF<5%) variants analysis with lipid traits in African Blacks.

https://doi.org/10.1371/journal.pone.0114618.t006

Functional annotation of the sequence variation

We used open-access database RegulomeDB (http://regulome.stanford.edu) to predict the potential implication of the identified genetic variation on the gene expression regulation. The RegulomeDB score of 1-5 is based on its strength of association with the gene regulation process; the lowest score represents the highest significant impact on regulation process (based on these features; expression quantitative trait loci (eQTL), transcription binding site or DNase hypersensitivity) while the highest score represents the least significant implication in regulation process. The RegulomeDB score for each variant is given in Tables 1 and 2. According to the RegulomeDB score, three variants in the 5′ flanking region (ApE173, ApE624/rs769446, and ApE832/rs405509), one intronic variant (ApE1231), and three variants in the 3′flanking region (ApE5223, ApE5231, and ApE5361/rs1081106) seem to affect gene expression as they have small scores (RegulomeDB score = 1–3). However, only two of these variants (APOE 624/rs769446 and APOE832/rs405509) had significant effect on LDL-C or ApoB or TG, and these two variants, APOE5223 and APOE5361/rs1081106, showed borderline effects on LDL-C and HDL-C, respectively. Although the remaining variants with strong regulatory effects (APOE173, APOE1231, and APOE5231) were not associated with lipid variation, they may yet have other biological consequences.

Discussion

The role of common APOE genetic variation in affecting interindividual variation in plasma cholesterol, especially LDL-C, in the general population is well established. Less clear, however, is if APOE genetic variation has also an impact on other major lipid traits, like plasma HDL-C and TG. Recent lipid GWAS indicate that in addition to LDL-C, APOE common variants are also associated with HDL-C and TG levels [17][18]. Since common variants explain only ∼25–30% of the genetic variance of each major lipid trait [18], it has been hypothesized that uncommon low-frequency and rare variants in candidate genes may explain part of the missing heritability, as it has already been shown for some lipid genes [30][36]. Thus, deep resequencing of the APOE gene is warranted to identify both uncommon and common variants that might affect plasma lipid profile. The objective of this study was to evaluate the ‘common disease common variants’ (CDCV) and ‘common disease rare variants’ (CDRV) hypotheses by sequencing the entire APOE gene in selected individuals (n = 190) with extreme HDL-C levels from two ethnic groups in the variant discovery stage and then genotyping common tagSNPs and relevant uncommon/rare variants in the full datasets (NHWs = 623, and Blacks = 788) to evaluate their association with lipid traits. To our knowledge, this is the first population-based association study designed to evaluate the effect of the full spectrum of APOE genetic variation on major plasma lipid traits and ApoB levels. Previously, sequencing of the APOE gene has been reported in two different studies [38][39] and by the 1000 Genome project in order to characterize its genetic variation in unselected individuals without regards to lipid levels. Furthermore, most of the previous studies have only evaluated the influence of APOE coding and promoter variants on lipid traits [30][42].

By sequencing ∼5.5 kb of the APOE gene region, including all four exons, three introns, and ∼1 kb in each flanking region in selected individuals with extreme HDL-C levels in both population groups, we identified a total of 40 variants, including 10 novel variants not previously reported. As expected African Blacks tend to have more population-specific variants (21/31 = 68%) as compared to NHWs (9/19 = 47%). In NHWs, the proportion of common and uncommon variants was similar (56% vs. 44%), while in African Blacks more uncommon variants were observed than common ones (70% vs. 30%) (Table S2 in S1 File). We observed more subjects carrying group-specific uncommon variants in the high HDL-C group than in the low HDL-C group in NHWs (21% vs. 2%; p = 0.0037) and in African Blacks (15% vs. 6%; p = 0.316), although the difference in Blacks was not statistically significant. Likewise, the cumulative uncommon/rare variant analysis using SKAT-O also showed significant association with HDL-C in NHWs (p = 0.0061; Table 5).

The established association of the E*2 (rs7412) and E*4 (rs426538) SNPs with LDL-C and ApoB [15], [43][44] was confirmed in our study in which E*2 was associated with lowering effect on LDL-C (p = 1.84E-07 in NHWs, p = 5.35E-07 in Blacks), and ApoB (p = 9.65E-13 in NHWs, p = 0.0356 in Blacks), while E*4 was associated with elevating effect on LDL-C in both population groups (p = 0.0103 in NHWs; p = 0.0317 in Blacks) and elevating ApoB in NHWs (p = 0.0005). Although E*4 did not achieve the nominal significance with ApoB in Blacks, it showed similar trend of association. Moreover, we have identified 8 additional variants (4 in NHWs and 4 in Blacks) that were associated with LDL-C independent of the E*2 and E*4 SNPs. The 4 LDL-significant variants in NHWs include APOE832/rs405509,APOE1163/rs4405509, APOE4310/rs199768005(Val254Glu) and APOE2440/rs769450. While the first 3 variants are associated with lowering effect on LDL-C, the last variant was associated with elevating effect (see Table 1). Among the 4 LDL-associated variants in Blacks, 3 (APOE2269/rs61357706, APOE2544/rs115299243, and APOE4036/rs769455 (Arg163Cys) were associated with low LDL-C and this association is more likely mediated by APOE4036/rs769455 (Arg163Cys) that has previously been associated with type III hyperlipoproteinemia [50][51]. The fourth variant, APOE4569 (exon4), was associated with high LDL-C (see Table 2). While the 4 LDL-significant variants observed in Blacks were not detected in NHWs, 3 of the 4 LDL-significant variants in NHWs were observed in Blacks (APOE832/rs405509, APOE1163/rs4405509, and APOE2440/rs769450) and they also showed suggestive associations with LDL-C.

Of the above-mentioned 8 significant variants independent of the E*2 and E*4 SNPs only 3 (APOE832/rs405509, APOE1163/rs440446, and APOE4036/rs769455 (Arg163Cys)) have been examined previously in relation to lipid traits. APOE832/rs405509 located in the putative promoter region has previously been shown to be associated with LDL-related traits (LDL-C, TC, and ApoB) [16][17], [45], APOE gene expression [46], myocardial infarction risk [47], and premature CHD [48]. Our findings confirm the potentially important role of this variant in LDL metabolism by observing significant associations with LDL-C. APOE1163/rs440446 was earlier reported to be associated with CHD risk [49] and our current finding with its association with LDL-C validates this link given the relation between LDL-C levels and CHD risk. The non-synonymous variant APOE4036/rs769455 (Arg163Cys) has previously been reported to be associated with type III hyperlipoproteinemia [50][51] and is probably the main contributor to the significance signal of the two other closely linked variants (APOE2269/rs61357706 and APOE2544/rs115299243) with LDL-C.

In addition to the known contribution of APOE to LDL-C, we have found multiple associations of common and uncommon variants with TG and HDL-C. One NHW-specific uncommon variant (APOE1575/rs769448) was associated with elevating effect on HDL-C (p = 0.0223) and one rare Black-specific variant (APOE618) was associated with extremely low HDL-C, (p = 0.001), implying the significant contribution of APOE uncommon/rare variants on plasma HDL-C variation. To our knowledge, these are novel associations and need to be confirmed in independent studies. Based on their locations (intron 1, and 5′ flanking region, respectively), and RegulomeDB scores [4], they may be moderately involved in gene expression regulation. Nine variants showed significant association with TG, including five in NHWs: APOE832/rs405509 (p = 0.002), APOE1163/rs440446 (p = 0.0012), APOE2440/rs769450 (p = 0.0022), APOE4310/rs199768005 (p = 0.028), and APOE4528/rs374329439 (p = 0.0218) and four in Blacks: APOE73/rs1081101 (p = 0.0115), APOE1279/rs877973 (p = 0.014) APOE2544/rs115299243 (p = 0.038), and APOE4036/rs769455 (p = 0.0343). Four of these variants are uncommon, including two present only in NHWs (APOE4310/rs199768005/Val254Glu, and APOE4528/rs374329439) and two present only in African Blacks (APOE2544/rs115299243, and APOE4036/rs769455/Arg163Cys). Two of these population-specific variants involving non-synonymous changes (Arg163Cys, and Val254Glu) have previously been reported to be associated with type III hyperlipoproteinemia either in E*2-independent (rs769455/Arg163Cys) [50][51] or E*2-dependent (rs199768005/Val254Glu) [41] fashion. In our population-based samples while Arg163Cys was associated with higher TG levels, Val254Glu was associated with lower TG levels. The latter observation may not be surprising given that this variant was associated with hypertriglyceridemia only among E*2 carriers [41] and all our 5 subjects with this mutation in our study were non-E*2 carriers. This also implies that Val254Glu variant may be protective in the absence of E*2. In accordance with our observations, APOE832/rs405509 [17] has been found previously to be associated with VLDL as an indicator of TG variation and APOE1163/rs440446 [49] has previously been found to be associated with TG variation. To our knowledge, the remaining five TG associations observed in this study have not been reported previously and await confirmation in future studies.

In summary, this is the first comprehensive study that has evaluated the association of APOE common and rare variation with plasma lipid traits in two ethnic groups. In addition to the known association of common APOE variation with LDL-C, we have found that uncommon APOE variants also affect LDL-C levels. Our data also indicate the contribution of APOE genetic variation in affecting HDL-C and TG levels in the general population. Strengths of our study include the use of two extreme lipid groups for resequencing from two ethnic groups and then genotyping of the entire sample sets for genotype-phenotype association analyses. Limitations of our study include the use of relatively small sample sizes for resequencing. Many of our significant findings with uncommon/rare variants should be considered provisional until replicated in independent and large data sets.

Supporting Information

S1 File.

Table S1. Demographic characteristics of the resequencing samples; Table S2. APOE sequencing variants identified in 95 NHWs and 95 African Blacks (n = 95); Table S3. Distribution of the sequence variants in the extreme HDL-C groups in NHWs (n = 95); Table S4. Distribution of the sequence variants in the extreme HDL-C groups in African Blacks (n = 95); Table S5. Single-site association analysis in NHWs (n = 623); Table S6. Single-site association analysis in African Blacks (n = 788); Table S7. 4-SNPs window haplotype-based association results for LDL-C, TG and HDL-C in NHWs (n = 623); Table S8. 4-SNPs window haplotype-based association results for ApoB in NHWs (n = 623); Table S9. Haplotype-based association summary of significant windows with LDL-C in NHWs (n = 623); Table S10. Haplotype-based association summary of significant windows with ApoB in NHWs (n = 623); Table S11. Haplotype-based association summary of significant windows with TG in NHWs (n = 623); Table S12. Haplotype-based association summary of significant windows with HDL-C in NHWs (n = 623); Table S13. 4-SNPs window haplotype-based association results for lipid traits in African Blacks (n = 788); Table S14. Haplotype-based association summary of significant windows with LDL-C in African Blacks (n = 788); Table S15. Haplotype-based association summary of significant windows with ApoB in African Blacks (n = 788); Table S16. Haplotype-based association summary of significant windows with TG in African Blacks (n = 788).

https://doi.org/10.1371/journal.pone.0114618.s001

(DOCX)

Author Contributions

Conceived and designed the experiments: MIK FYD. Performed the experiments: ZHR FYD FW. Analyzed the data: ZHR XBW MIK FYD. Contributed reagents/materials/analysis tools: CHB JEH RFH MMB FYD MIK. Wrote the paper: ZHR MIK FYD. Provided critical revisions: XBW DP VN CHB JEH RFH MMB. Interpreted the results: ZHR XBW DP VN CHB JEH RFH MMB MIK FYD.

References

  1. 1. Ordovas JM, Shen J (2008) Gene-environment interactions and susceptibility to metabolic syndrome and other chronic diseases. J Periodontol 79:1508–1513.
  2. 2. Mahley RW, Innerarity TL, Rall SC Jr, Weisgraber KH (1984) Plasma lipoproteins: apolipoprotein structure and function. J Lipid Res 25:1277–1294.
  3. 3. Mahley RW (1988) Apolipoprotein E: Cholesterol transport protein with expanding role in cell biology. Science 240:622–630.
  4. 4. Miettinen TA, Gylling H, Vanhanen H, Ollus A (1992) Cholesterol absorption, elimination, and synthesis related to LDL kinetics during fat intake in men with different apoprotein E phenotypes. Arterioscler Thromb 12:1044–1052.
  5. 5. Pitas RE, Boyles JK, Lee SH, Foss D, Mahley RW (1987) Astrocytes synthesize apolipoprotein E and metabolize apolipoprotein E-containing lipoproteins. Biochim Biophys Acta 917:148–161.
  6. 6. Boyles JK, Zoellner CD, Anderson LJ, Kosik LM, Pitas RE, et al. (1989) A role for apolipoprotein E, apolipoprotein A-I, and low-density lipoprotein receptors in cholesterol transport during regeneration and remyelination of the rat sciatic nerve. J Clin Invest 83:1015–1031.
  7. 7. Westlye LT, Reinvang I, Rootwelt H, Espeseth T (2012) Effects of APOE on brain white matter microstructure in healthy adults. Neurology 79:1961–1969.
  8. 8. Dumanis SB, DiBattista AM, Miessau M, Moussa CE, Rebeck GW (2013) APOE genotype affects the pre-synaptic compartment of glutamatergic nerve terminals. J Neurochem 124:4–14.
  9. 9. De Blasi S, Montesanto A, Martino C, Dato S, De Rango F, et al. (2009) APOE polymorphism affects episodic memory among non demented elderly subjects. Experimental Gerontology 44:224–227.
  10. 10. Savitz J, Solms M, Ramesar R (2006) Apolipoprotein E variants and cognition in healthy individuals: a critical opinion. Brain Research Reviews 51:125–135.
  11. 11. Mahley RW, Rall SC Jr (2000) Apolipoprotein E: far more than a lipid transport protein. Annu Rev Genomics Hum Genet 1:507–537.
  12. 12. Kuhlmann I, Minihane AM, Huebbe P, Nebel A, Rimbach G (2010) Apolipoprotein E genotype and hepatitis C, HIV and herpes simplex disease risk: a literature review. Lipids Health Dis
  13. 13. Nikodemova M, Finn L, Mignot E, Salzieder N, Peppard PE (2013) Association of Sleep Disordered Breathing and Cognitive Deficit in APOE ε4 Carriers. Sleep 36:873–880.
  14. 14. Das HK, McPherson J, Bruns GA, Karathanasis SK, Breslow JL (1985) Isolation, characterization and mapping to chromosome 19 of the human apolipoprotein E gene. J Biol Chem 260:6240–6247.
  15. 15. Sanna S, Li B, Mulas A, Sidore C, Kang HM, et al. (2011) Fine mapping of five loci associated with low-density lipoprotein cholesterol detects variants that double the explained heritability. PLoS Genet 7:e1002198.
  16. 16. Ken-Dror G, Talmud PJ, Humphries SE, Drenos F (2010) APOE/C1/C4/C2 gene cluster genotypes, haplotypes and lipid levels in prospective coronary heart disease risk among UK healthy men. Mol Med 16:389–399.
  17. 17. Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, et al. (2013) Discovery and refinement of loci associated with lipid levels. Nat Genet 11:1274–1283.
  18. 18. 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.
  19. 19. Hamman RF, Marshall JA, Baxter J, Kahn LB, Mayer EJ, et al. (1989) Methods and prevalence of non-insulin dependent diabetes mellitus in a biethnic Colorado population. The San Luis Valley Diabetes Study. Am J Epidemiol 129:295–311.
  20. 20. Rewers M, Shetterly SM, Hoag S, Baxter J, Marshall J, et al. (1993) Is the risk of coronary heart disease lower in Hispanics than in non-Hispanic whites? The San Luis Valley Diabetes Study. Ethn Dis 3:44–54.
  21. 21. Bunker CH, Ukoli FA, Matthews KA, Kriska AM, Huston SL, et al. (1995) Weight threshold and blood pressure in a lean black population. Hypertension 2:616–623.
  22. 22. Bunker CH, Ukoli FA, Okoro FI, Olomu AB, Kriska AM, et al. (1996) Correlates of serum lipids in a lean black population. Atherosclerosis 123:215–225.
  23. 23. Kamboh MI, Rewers M, Aston CE, Hamman RF (1997) Plasma apolipoprotein A-I, apolipoprotein B, and lipoprotein(a) concentrations in normoglycemic Hispanics and non-Hispanic whites from the San Luis Valley, Colorado. Am J Epidemiol 146:1011–1018.
  24. 24. Harris MR, Bunker CH, Hamman RF, Sanghera DK, Aston CE, et al. (1998) Racial differences in the distribution of a low-density lipoprotein receptor-related protein (LRP) polymorphism and its association with serum lipoprotein, lipid and apolipoprotein levels. Atherosclerosis 137:187–195.
  25. 25. Demirci FY, Dressen AS, Hamman RF, Bunker CH, Kammerer CM, et al. (2010) Association of a common G6PC2 variant with fasting plasma glucose levels in non-diabetic individuals. Ann Nutr Metab 56:59–64.
  26. 26. Bryant EK, Dressen AS, Bunker CH, Hokanson JE, Hamman RF, et al. (2013) A multiethnic replication study of plasma lipoprotein levels-associated SNPs identified in recent GWAS. PLoS One 8:e63469
  27. 27. Liu JZ, McRae AF, Nyholt DR, Medland SE, Wray NR, et al. (2010) A versatile gene-based test for genome-wide association studies. Am J Hum Genet 87:139–145.
  28. 28. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 57:289–300.
  29. 29. Lake SL, Lyon H, Tantisira K, Silverman EK, Weiss ST, et al. (2003) Estimation and tests of haplotype-environment interaction when linkage phase is ambiguous. Hum Hered 55:56–65.
  30. 30. Lee S, Wu MC, Lin X (2012) Optimal tests for rare variant effects in sequencing association studies. Biostatistics 13(4):762–775.
  31. 31. Johansen CT, Wang J, Lanktree MB, Cao H, McIntyre AD, et al. (2010) Excess of rare variants in genes identified by genome-wide association study of hypertriglyceridemia. Nat Genet 42:684–687.
  32. 32. Johansen CT, Hegele RA (2012) The complex genetic basis of plasma triglycerides. Curr Atheroscler Rep 14:227–234.
  33. 33. Cohen JC, Kiss RS, Pertsemlidis A, Marcel YL, McPherson R, et al. (2004) Multiple rare alleles contribute to low plasma levels of HDL cholesterol. Science. 305:869–872.
  34. 34. Cohen J, Pertsemlidis A, Kotowski IK, Graham R, Garcia CK, et al. (2005) Low LDL cholesterol in individuals of African descent resulting from frequent nonsense mutations in PCSK9. Nat Genet 37:161–165.
  35. 35. Kotowski IK, Pertsemlidis A, Luke A, Cooper RS, Vega GL, et al. (2006) A spectrum of PCSK9 alleles contributes to plasma levels of low-density lipoprotein cholesterol. Am J Hum Genet 78:410–422.
  36. 36. Romeo S, Pennacchio LA, Fu Y, Boerwinkle E, Tybjaerg-Hansen A, et al. (2007) Population-based resequencing of ANGPTL4 uncovers variations that reduce triglycerides and increase HDL. Nat Genet 39:513–516.
  37. 37. Pollin TI, Damcott CM, Shen H, Ott SH, Shelton J, et al. (2008) A null mutation in human APOC3 confers a favorable plasma lipid profile and apparent cardioprotection. Science 322:1702–1705.
  38. 38. Nickerson DA, Taylor SL, Fullerton SM, Weiss KM, Clark AG, et al. (2000) Sequence diversity and large-scale typing of SNPs in the human apolipoprotein E gene. Genome Res 10:1532–1545.
  39. 39. Fullerton SM, Clark AG, Weiss KM, Nickerson DA, Taylor SL, et al. (2000) Apolipoprotein E variation at the sequence haplotype level: implications for the origin and maintenance of a major human polymorphism. Am J Hum Genet 67:881–900.
  40. 40. Ozturk Z, Enkhmaa B, Shachter NS, Berglund L, Anuurad E (2010) Integrated role of two apoliprotein E polymorphisms on apolipoprotein B levels and coronary artery disease in a biethnic population. Metab Syndr Relat Disord 8:531–538.
  41. 41. van den Maagdenberg AM, Weng W, de Bruijn IH, de Knijff P, Funke H, et al. (1993) Characterization of five new mutants in the carboxyl-terminal domain of human apolipoprotein E: no cosegregation with severe hyperlipidemia. Am J Hum Genet 52:937–946.
  42. 42. Akanji AO, Suresh CG, Fatania HR, Al-Radwan R, Zubaid M (2007) Associations of apolipoprotein E polymorphism with low-density lipoprotein size and subfraction profiles in Arab patients with coronary heart disease. Metabolism 56:484–490.
  43. 43. Medina-Urrutia AX, Cardoso-Saldaña GC, Zamora-González J, Liria YK, Posadas-Romero C (2004) Apolipoprotein E polymorphism is related to plasma lipids and apolipoproteins in Mexican adolescents. Hum Biol 76:605–614.
  44. 44. Lucatelli JF, Barros AC, Silva VK, Machado Fda S, Constantin PC, et al. (2011) Genetic influences on Alzheimer's disease: Evidence of interactions between the genes APOE, APOC1 and ACE in a sample population from the South of Brazil. Neurochem Res 36:1533–1539.
  45. 45. Klos KL, Sing CF, Boerwinkle E, Hamon SC, Rea TJ, et al. (2006) Consistent effects of genes involved in reverse cholesterol transport on plasma lipid and apolipoprotein levels in CARDIA participants. Arterioscler Thromb Vasc Biol 26:1828–1836.
  46. 46. Artiga MJ, Bullido MJ, Sastre I, Recuero M, García MA, et al. (1998) Allelic polymorphisms in the transcriptional regulatory region of apolipoprotein E gene. FEBS Lett 421:105–108.
  47. 47. Lambert JC, Brousseau T, Defosse V, Evans A, Arveiler D, et al. (2000) Independent association of an APOE gene promoter polymorphism with increased risk of myocardial infarction and decreased APOE plasma concentrations-the ECTIM Study. Hum Mol Genet 9:57–61.
  48. 48. Viitanen L, Pihlajamäki J, Miettinen R, Kärkkäinen P, Vauhkonen I, et al. (2001) Apolipoprotein E gene promoter (-219G/T) polymorphism is associated with premature coronary heart disease. J Mol Med (Berl) 79:732–737.
  49. 49. Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, et al. (2008) Gender differences in genetic risk profiles for cardiovascular disease. PLoS One 3:e3615.
  50. 50. Havel RJ, Kotite L, Kane JP, Tun P, Bersot T (1983) Atypical familial dysbetalipoproteinemia associated with apolipoprotein phenotype E3/3. J Clin Invest 72:379–387.
  51. 51. Rall SC Jr, Newhouse YM, Clarke HR, Weisgraber KH, McCarthy BJ, et al. (1989) Type III hyperlipoproteinemia associated with apolipoprotein E phenotype E3/3. Structure and genetics of an apolipoprotein E3 variant. J Clin Invest 83:1095–1101.