OUP user menu

Eight genetic loci associated with variation in lipoprotein-associated phospholipase A2 mass and activity and coronary heart disease: meta-analysis of genome-wide association studies from five community-based studies

Harald Grallert, Josée Dupuis, Joshua C. Bis, Abbas Dehghan, Maja Barbalic, Jens Baumert, Chen Lu, Nicholas L. Smith, André G. Uitterlinden, Robert Roberts, Natalie Khuseyinova, Renate B. Schnabel, Kenneth M. Rice, Fernando Rivadeneira, Ron C. Hoogeveen, João Daniel Fontes, Christa Meisinger, John F. Keaney, Rozenn Lemaitre, Yurii S. Aulchenko, Ramachandran S. Vasan, Stephen Ellis, Stanley L. Hazen, Cornelia M. van Duijn, Jeanenne J. Nelson, Winfried März, Heribert Schunkert, Ruth M. McPherson, Heide A. Stirnadel-Farrant, Bruce M. Psaty, Christian Gieger, David Siscovick, Albert Hofman, Thomas Illig, Mary Cushman, Jennifer F. Yamamoto, Jerome I. Rotter, Martin G. Larson, Alexandre F.R. Stewart, Eric Boerwinkle, Jacqueline C.M. Witteman, Russell P. Tracy, Wolfgang Koenig, Emelia J. Benjamin, Christie M. Ballantyne
DOI: http://dx.doi.org/10.1093/eurheartj/ehr372 238-251 First published online: 14 October 2011

Abstract

Aims Lipoprotein-associated phospholipase A2 (Lp-PLA2) generates proinflammatory and proatherogenic compounds in the arterial vascular wall and is a potential therapeutic target in coronary heart disease (CHD). We searched for genetic loci related to Lp-PLA2 mass or activity by a genome-wide association study as part of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium.

Methods and results In meta-analyses of findings from five population-based studies, comprising 13 664 subjects, variants at two loci (PLA2G7, CETP) were associated with Lp-PLA2 mass. The strongest signal was at rs1805017 in PLA2G7 [P = 2.4 × 10−23, log Lp-PLA2 difference per allele (beta): 0.043]. Variants at six loci were associated with Lp-PLA2 activity (PLA2G7, APOC1, CELSR2, LDL, ZNF259, SCARB1), among which the strongest signals were at rs4420638, near the APOE–APOC1–APOC4–APOC2 cluster [P = 4.9 × 10−30; log Lp-PLA2 difference per allele (beta): −0.054]. There were no significant gene–environment interactions between these eight polymorphisms associated with Lp-PLA2 mass or activity and age, sex, body mass index, or smoking status. Four of the polymorphisms (in APOC1, CELSR2, SCARB1, ZNF259), but not PLA2G7, were significantly associated with CHD in a second study.

Conclusion Levels of Lp-PLA2 mass and activity were associated with PLA2G7, the gene coding for this protein. Lipoprotein-associated phospholipase A2 activity was also strongly associated with genetic variants related to low-density lipoprotein cholesterol levels.

  • Genome-wide association
  • Inflammation
  • Lipoprotein-associated phospholipase A2

Introduction

Lipoprotein-associated phospholipase A2 (Lp-PLA2) is a 45.4-kDa calcium-independent member of the phospholipase A2 family, which is secreted by leucocytes1 and has been detected in rabbit and human atherosclerotic lesions.2 In the bloodstream, two-thirds of Lp-PLA2 circulates primarily bound to LDL; the remaining third is distributed between HDL and VLDL. Circulating Lp-PLA2 can be measured by different assays ascertaining mass or activity of Lp-PLA2. However, there is only moderate correlation between mass-based and activity-based measurements of Lp-PLA2 (r = 0.51),3 and the independent role of these two measures in cardiovascular disease (CVD) is unclear.

Several lines of evidence suggest that Lp-PLA2 is associated with the development of atherosclerotic disease. Lipoprotein-associated phospholipase A2 generates pro-inflammatory and pro-atherogenic compounds in the arterial vascular wall. A large meta-analysis with almost 80 000 participants in 32 prospective studies showed that high levels of Lp-PLA2 mass and activity were associated with the risk for coronary heart disease (CHD), stroke, and cardiovascular mortality.3 Therefore, Lp-PLA2 might represent an emerging biomarker for improved cardiovascular risk assessment in clinical practice and a potential therapeutic target for primary or secondary prevention of CVD.

Although familial factors explain about one-half and one-quarter of the variance in Lp-PLA2 activity and mass, respectively,4,5 few genetic determinants of Lp-PLA2 have been identified so far.6 Therefore, to better understand genetic control of Lp-PLA2, we investigated genetic loci related to Lp-PLA2 mass or activity by conducting genome-wide association (GWA) analyses in five population-based studies, as part of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium.7

Methods

For more detailed information on the Methods, see Supplementary material online.

Study data for meta-analysis

We used data from five community-based cohorts in the USA and Europe: Atherosclerosis Risk in Communities (ARIC) Study, Cardiovascular Health Study (CHS), Framingham Heart Study (FHS), Rotterdam Study (RS), and Monitoring of Trends and Determinants in Cardiovascular Disease (MONICA)/Cooperative Health Research in the Region of Augsburg Study (KORA). Baseline clinical and demographic characteristics were assessed at the time of cohort entry for CHS, KORA, and RS and at the time of biomarker collection for ARIC and FHS. All participants gave written informed consent, including consent for genetic studies. All studies received approval from local ethical oversight committees.

Associations between significant single nucleotide polymorphisms (SNPs) and CHD/coronary artery disease (CAD) were assessed in the Coronary Artery Disease Genome-Wide Replication and Meta-Analysis (CARDIoGRAM) consortium data. CARDIoGRAM includes over 22 000 cases with CAD, MI, or both and over 60 000 controls from all published and several unpublished genome-wide association studies (GWAS) in individuals of European descent.8

Lipoprotein-associated phospholipase A2 mass concentrations were measured in ARIC, CHS, FHS, and KORA using a commercially available sandwich enzyme immunoassay (PLAC® test; diaDexus, Inc., San Francisco, CA, USA). Lipoprotein-associated phospholipase A2 activity was measured in CHS, FHS, and RS on microtiter plates by colourimetric or radioactive substrate methods (diaDexus CAM Kit, diaDexus, Inc., San Francisco, CA, USA or Perkin Elmer Life Sciences, Inc., Waltham, MA, USA).

Genotyping was performed in each cohort using high-density SNP marker platforms. Genotypes were imputed to ∼2.5 million HapMap SNPs (platform and genetic analysis details provided in Supplementary material online).

Statistical analyses

Associations between genotype and Lp-PLA2 mass and activity were analysed separately within each cohort, using linear regression of natural log-transformed phenotype on number (or imputed dose) of reference alleles, i.e. an additive model (for more details, see Supplementary material online). All analyses were adjusted for age, sex, and if applicable, recruitment site. To assess whether classical cardiovascular risk factors confounded the associations, analyses were additionally adjusted for diabetes, lipid-lowering medication, antihypertensive treatment, aspirin intake ≥3 per week, current smoking, hormone replacement therapy (HRT; coded men, women without HRT, women with HRT), body mass index (BMI), systolic blood pressure, diastolic blood pressure, triglycerides, waist circumference, HDL-C, LDL-C, and prevalent CVD. For the most significant SNPs, we also analysed interactions between genetic variant and age, sex, BMI, or smoking status.

To combine results across cohorts, we performed an inverse variance–weighted meta-analysis using the software METAL,9 which was specifically developed for meta-analysing GWAS results. A meta-analysis is a statistical method that combines analyses from different independent studies to give an overall effect measure. In the case of meta-analysing different GWAS, this approach increases the power to detect significant genetic variants compared with analysing each GWAS alone.10 Cohort-specific standard errors were adjusted using genomic control.11 In GWAS, we chose P = 5 × 10−8 as the threshold for significance12 (see Supplementary material online). To determine whether the multiple SNPs associated in the respective region are due to linkage disequilibrium (LD) with the top SNP or if there are multiple independent signals, we performed a meta-analysis based on models adjusted additionally for the SNPs with the smallest P-values.

Results

Sample size, demographics, and laboratory characteristics of each cohort are presented in Table 1. Median values for Lp-PLA2 mass (241.8 ng/mL in FHS to 413.0 ng/mL in ARIC) and activity (38.2 nmol/min/mL in CHS to 153.9 nmol/min/mL in FHS) varied substantially across cohorts, but much of these differences may be explained by the different assays used.

View this table:
Table 1

Characteristics of study participants (n = 13 664)

Study sample
ARICCHSFHSRSKORA
Number7983217690915381202
Age (years)58.6 ± 5.572.3 ± 5.349.3 ± 13.869.1 ± 9.053.5 ± 8.9
Male (%)6238.546.939.552.6
Current smoking (%)24.211.315.723.517.7
BMI (kg/m2)28.0 ± 4.826.3 ± 4.527.4 ± 5.426.2 ± 3.627.4 ± 4.1
Waist (cm)99.8 ± 13.493.1 ± 12.996.0 ± 15.090.2 ± 11.091.0 ± 12.7
Systolic blood pressure (mmHg)125 ± 19135 ± 21121 ± 17138 ± 22134 ± 19
Diastolic blood pressure (mmHg)73 ± 1070 ± 11 75 ± 1073 ± 1182 ± 11
Hypertension treatment (%)34.230.219.331.117.5
Lipid-lowering medication (%)9.44.413.3 2.2 4.6
Aspirin ≥3 per week (%)28.118.5 5.1
Total cholesterol (mg/dL)214 ± 42213 ± 39194 ± 37256 ± 50238 ± 44
HDL-C (mg/dL)44 ± 1555 ± 16 54 ± 1752 ± 1554 ± 16
LDL-C (mg/dL)138 ± 36130 ± 35115 ± 32147 ± 41
Total/HDL-C ratio5.3 (2.2)4.1 ± 1.23.9 ± 1.45.3 ± 1.64.8 ± 1.8
Triglycerides (mg/dL)159 ± 130125 ± 90187 ± 150
Prevalent diabetes (%)18.611.77.5 9.9 4.7
Prevalent CVDa (%) 0 06.2 7.4 1.9
Hormone-replacement therapy (in women) (%)29.6 8.616.421.328.6
Lp-PLA2 mass (ng/mL)413.0 (328.5/513.5)329.9 (261.4/407.4)241.8 (208.0/299.3)269.0 (214.0/330.5)
Lp-PLA2 activity (nmol/mL/min)38.2 (30.9/47.0)153.9 (129.2/181.2)44.5 (36.3/51.3)
  • Mean ± SD for continuous, per cent for categorical variables, median (25th/75th percentile) for Lp-PLA2 concentrations.

  • aHistory of myocardial infarction, angina, coronary revascularization, stroke, or transient ischemic attack.

Genome-wide association of lipoprotein-associated phospholipase A2 mass and activity

The meta-analysis included 2 661 766 SNPs from one or more studies. Genomic control (λgc) parameters were small (all λgc≤ 1.058), suggesting negligible inflation in the type-I error rate.

Figure 1A and B illustrates the primary findings from the meta-analysis; details are in Table 2. In total, 49 SNPs were significantly associated with Lp-PLA2 mass and 59 with Lp-PLA2 activity. For Lp-PLA2 mass, these signals clustered at two loci on chromosomes 6 and 16. For activity, the SNPs clustered at loci on chromosomes 1, 6, 11, 12, and (two clusters) 19 (Supplementary material online, Tables S1 and S2).

View this table:
Table 2

Association of the top single nucleotide polymorphisms in eight loci with log-transformed lipoprotein-associated phospholipase A2, adjusted for age, sex, and site (if necessary)

SNPMassActivity
ARICCHSFHSKORAMeta-analysisCHSFHSRSMeta-analysis
rs4420638 (G/A)Sample size79832176899120212 11632176899153811 654
Chromosome: 19P-value4.3E–041.9E–017.9E–015.1E–028.6E–025.4E–061.6E–193.9E–104.9E–30
Location: 50114786Beta−0.075−0.0230.002−0.034−0.0090.0710.0480.0990.054
Gene: APOC1Standard error0.0210.0170.0060.0170.0050.0160.0050.0160.005
Effect allele: Ar20.01400.00050.00090.00320.00530.01620.0240
Effect allele frequency0.820.830.840.820.830.840.80
I20.8030.816
rs1805017 (C/T)Sample size79832176896120212 11332176896153811 651
Chromosome: 6P-value2.6E–029.5E–085.0E–141.4E–032.4E–231.8E–021.8E–047.1E–022.4E–06
Location: 46792181Beta0.0430.0510.0400.0520.043−0.021−0.017−0.019−0.018
Gene: PLA2G7Standard error0.0190.0100.0050.0160.0040.0090.0040.0100.004
Effect allele: Tr20.00560.00850.00970.00850.00160.00620.0014
Effect allele frequency0.250.270.260.230.270.260.26
I200
rs7528419 (G/A)Sample sizeNA32176907NA10 12432176907153811 662
Chromosome: 1P-valueNA2.6E–034.9E–03NA7.1E–051.2E–059.4E–122.6E–031.3E–17
Location: 109618715BetaNA0.0290.017NA0.0200.0410.0340.0330.035
Gene: CELSR2Standard errorNA0.0100.006NA0.0050.0090.0050.0110.004
Effect allele: Ar2NA0.00260.0022NA0.00540.01070.0051
Effect allele frequencyNA0.780.79NA0.780.790.76
I20.0960
rs6511720 (G/T)Sample size79832176907120212 12432176907153811 662
Chromosome: 19P-value2.4E–011.6E–033.1E–022.8E–015.5E–059.9E–103.1E–035.2E–032.6E–11
Location: 11063306Beta−0.043−0.039−0.025−0.035−0.0320.0710.0290.0400.045
Gene: LDLRStandard error0.0360.0120.0120.0330.0080.0120.0100.0140.007
Effect allele: Tr20.00160.00290.00160.00100.01010.00450.0043
Effect allele frequency0.100.130.100.080.130.100.12
I200.739
rs964184 (G/C)Sample size79832176907120212 12432176907153811 662
Chromosome: 11P-value2.7E–011.6E–021.3E–012.8E–015.3E–031.4E–036.3E–084.1E–028.4E-11
Location: 116154127Beta−0.026−0.034−0.010−0.020−0.0160.0440.0300.0290.032
Gene: ZNF259Standard error0.0230.0140.0070.0190.0060.0140.0060.0140.005
Effect allele: Cr20.00140.00170.00130.00100.00300.00880.0020
Effect allele frequency0.850.870.860.870.870.860.87
I200
rs7756935 (C/A)Sample size79832176907120212 12432176907153811 662
Chromosome: 6P-value2.8E–019.2E–021.5E–013.3E–017.2E–011.5E–034.9E–062.1E–041.3E–10
Location: 46782984Beta0.0210.017−0.0090.016−0.0020.0300.0230.0430.027
Gene: PLA2G7Standard error0.0200.0100.0060.0160.0050.0100.0050.0120.004
Effect allele: Ar20.00130.00090.00120.00080.00300.00720.0081
Effect allele frequency0.780.790.810.790.790.810.8
I20.5680.255
rs10846744 (G/C)Sample sizeNANA690412028106NA690415388442
Chromosome: 12P-valueNANA3.3E–057.1E–014.9E–05NA2.0E–089.9E–026.1E–09
Location: 123878378BetaNANA0.0260.0070.025NA0.0300.0230.029
Gene: SCARB1Standard errorNANA0.0060.0200.006NA0.0050.0140.005
Effect allele: Cr2NANA0.00410.0001NA0.00860.0011
Effect allele frequencyNANA0.160.16NA0.160.14
I200.066
rs247616 (C/T)Sample size79832176907120212 12432176907153811 662
Chromosome: 16P-value3.3E–013.1E–016.0E–076.1E–042.5E–086.4E–018.4E–012.1E–014.2E–01
Location: 55547091Beta0.0220.0080.0260.0590.023−0.004−0.001−0.013−0.003
Gene: CETPStandard error0.0220.0080.0050.0170.0040.0080.0040.0100.004
Effect allele: Tr20.00110.00030.00430.00980.00010.00340.0004
Effect allele frequency0.330.340.320.340.340.320.33
I20.6180
  • Significant SNPs in the meta-analysis with respective data from each cohort are marked bold.

  • r2: increase in r2 after adding the SNP to a linear regression model with log Lp-PLA2 as outcome and age and sex (+site, if necessary) in the model.

  • I2: measure of heterogeneity.

Figure 1

Association of log-transformed lipoprotein-associated phospholipase A2 mass (A) and lipoprotein-associated phospholipase A2 activity (B) concentrations and 2 661 766 single nucleotide polymorphisms displayed per chromosome and region. The dashed line indicates the significance threshold of P < 5 × 10−8.

Lipoprotein-associated phospholipase A2 mass meta-analysis

Of the 49 SNPs significantly associated with Lp-PLA2 mass, 47 were within the Lp-PLA2 gene (PLA2G7), at chromosome 6p21.2-p12. The strongest association [P = 2.4 × 10−23; log Lp-PLA2 difference per allele (beta): 0.043] was for rs1805017 (Figure 2A), which is located within the coding region of PLA2G7 in exon 4 and leads to an arginine-to-histidine substitution at position 92 (Arg92His). Supplementary material online, Table S3 shows unadjusted median values of Lp-PLA2 mass by rs1805017 genotype for each cohort.

Figure 2

Regional plots for top single nucleotide polymorphisms rs4420638 (A), rs1805017 (B), rs964184 (C), rs7756935 (D), rs10846744 (E), and rs247616 (F). LD: red: r2 (with primary SNP) >0.8; orange: r2 between 0.5 and 0.8; yellow: r2 between 0.2 and 0.5; white: r2 < 0.2.

Based on pairwise LD (based on HapMap CEU) coefficients, 27 of the 47 signals with PLA2G7 lie within the same highly preserved haplotype block (Lewontin's D′ > 0.99). The remaining 21 associated SNPs in this region were correlated with rs1805017 (r² < 0.83).

Lp-PLA2 mass was also associated with SNPs on chromosome 16 within the cholesteryl ester transfer protein (CETP) gene (CETP; Figure 2F): rs247616 (P = 2.5 × 10−8; beta: 0.023).

There were no significant interactions between age, sex, BMI, or smoking and the two most significant SNPs in each region (rs1805017 and rs247616) in relation to Lp-PLA2 mass (Supplementary material online, Table S4). The association with Lp-PLA2 remained significant for these two SNPs after multivariableadjustment with similar effect sizes and P-values (Supplementary material online, Table S5).

Lipoprotein-associated phospholipase A2 activity meta-analysis

The smallest P-value for association with Lp-PLA2 activity (P = 4.9 × 10−30; beta: −0.054) was for rs4420638 (Figure 2B), located on chromosome 19 near the APOE–APOC1–APOC4–APOC2 cluster and ∼11 kb proximal to the APOE-ɛ4 allele. The association was similar in all three studies, with beta coefficients of −0.099 to −0.048 (Table 2). Other SNPs with significant associations were rs7528419 (CELSR2), rs6511720 (LDLR), rs964184 (ZNF259), and rs10846744 (SCARB1) (Table 2). Lipoprotein-associated phospholipase A2 activity was also associated with PLA2G7. The strongest association was for rs7756935 (P = 1.3 × 10−10; beta: −0.027), which was in perfect LD with another functional SNP within the PLA2G7 gene, rs1051931 in exon 11 (Ala379Val).

There were no significant gene–environment interactions of age, sex, BMI, and smoking status with the top SNPs in six loci in relation to log-transformed Lp-PLA2 activity (Supplementary material online, Table S4).

Association of top single nucleotide polymorphisms with prevalent coronary heart disease/coronary artery disease

For all but two (rs247616, CETP; rs6511720, LDLR) of the top SNPs, the association to CHD/CAD could be assessed in the CARDIoGRAM project. Four of the top SNPs related to Lp-PLA2 activity [rs964184 (ZNF259), rs4420638 (APOC1), rs7528419 (CELSR2), rs10846744 (SCARB1)] were significantly associated with prevalent CHD/CAD (Figure 3). The strongest association was found for rs964184 with a P = 8.0 × 10−10, with an odds ratio of 1.13 per G allele (95% confidence interval 1.09–1.18). The two SNPs in the PLA2G7 gene most significantly associated with Lp-PLA2 mass or activity were not significantly associated with prevalent CHD/CAD in CARDIoGRAM. The estimated changes in Lp-PLA2 mass or activity per risk allele revealed in our meta-analysis were combined with the CHD risk estimations related to Lp-PLA2 mass or activity drawn from a published meta-analysis3 to estimate the expected CHD risk changes per allele. For SD values (necessary to calculate the expected risks), we chose the respective data from FHS, as of the five cohorts FHS had the largest number of participants. This approach revealed CHD risk increases between 0.8 and 2.1% per allele and were found to be mostly smaller than that found in the CARDIoGRAM data. The CARDIoGRAM study had ≥70% power to detect CHD risk changes of 3.5% or more (Supplementary material online, Table S6).

Figure 3

Association of the top single nucleotide polymorphisms in six loci with prevalent CHD/CAD in the CARDIoGRAM consortium.

Discussion

Meta-analysis of GWAS from five population-based studies comprising 13 664 subjects suggests that genes influencing Lp-PLA2 mass and activity differ. While enzyme concentration was regulated mainly by the gene coding for Lp-PLA2 (PLA2G7), genetic variants involved in lipid metabolism were strongly associated with Lp-PLA2 activity.

Lipoprotein-associated phospholipase A2 mass

Lipoprotein-associated phospholipase A2 mass was dependent on variants within PLA2G7. Previous studies identified numerous SNPs within PLA2G7: variants in exon 9 (Val279Phe; rs45619133), exon 11 (Val379Ala; rs1051931), exon 7 (Ile198Thr, Iso195Thr; rs1805018), and exon 4 (Arg92His; rs1805017).13 A prior GWA of Framingham data accessed through dbGaP also identified multiple SNPs clustered in or near PLA2G7 that were associated with Lp-PLA2 mass.6

Lipoprotein-associated phospholipase A2 mass was also influenced by rs247616 (CETP).14 The T allele of rs247616 was associated with higher HDL-C concentration14 and higher mean Lp-PLA2 mass. This finding is in contrast to most epidemiological studies, including ARIC15and Framingham,6 which showed an inverse relation between HDL-C level and Lp-PLA2 mass.3 The mechanism(s) by which reduced CETP activity leading to increased HDL-C would also increase Lp-PLA2 mass is unclear.

Lipoprotein-associated phospholipase A2 activity

All the SNPs associated with Lp-PLA2 activity in our meta-analysis also determine plasma lipoprotein concentrations except for those in PLA2G7. The strongest association for Lp-PLA2 activity was for rs4420638 (APOC1), previously associated with LDL-C concentrations.14,16 rs964184 (ZNF259), near the APOA5–APOA4–APOC3–APOA1 cluster, was associated with increased Lp-PLA2 activity and previously with increased triglycerides;14 elevated triglycerides are associated with increased small dense LDLs, which contain increased Lp-PLA2 mass.17

rs7528419 (CELSR2) located at 1q13.3, which also contains the genes PSRC1 and SORT1, has strong associations with LDL-C levels and CHD.14,16,18 Based on HapMap data, rs7528419 was highly correlated with rs599839 (PSRC1) and rs646776 (CELSR2); minor alleles of rs599839 and rs646776 have been associated with lower LDL-C14,16,18 and lower CHD risk.19 The minor allele of rs7528419 also was associated with lower Lp-PLA2 activity, consistent with the reduction in LDL-C in other studies.14,16,18 rs7528419, or a SNP in LD, might have functional effects due to its location in 3′ UTR of the CELSR2 gene (encoding for cadherin, EGF LAG seven-pass G-type receptor 2), which includes binding sites for transcriptional factors (e.g. Oct-1).

rs6511720 (LDLR) was also associated with lower Lp-PLA2 activity and lower LDL-C.18 It was highly correlated with rs2228671 (LDLR) (r2 = 0.734; D′ = 0.899), which was related to lower LDL-C and reduced CAD risk.20

rs10846744, an intronic SNP within the SCARB1 gene, which encodes for scavenger receptor type B class 1 (SRB1), a major receptor for HDL, was associated with Lp-PLA2 activity. Although this SNP was not associated with lipids in this study, other SNPs in SCARB1 have been shown to be associated with levels of HDL-C.21,22 SCARB1 may also modulate Lp-PLA2 activity through its role as a scavenger receptor for oxidized LDL, which has been shown to increase Lp-PLA2 secretion by human macrophages.23 SCARB1 is expressed in macrophages and has been shown to bind avidly to oxidized LDL.24

Finally, rs7756935 (PLA2G7), associated with increased Lp-PLA2 activity, was in complete LD with rs1051931 (PLA2G7; Ala379Val).24 The V379 allele resulted in two-fold lower Lp-PLA2 activity in in vitro studies25 but increased plasma Lp-PLA2 activity in epidemiological studies.26,27

Association of top single nucleotide polymorphisms with prevalent coronary heart disease/coronary artery disease

Of the top polymorphisms associated with Lp-PLA2 activity, four (in ZNF259, APOC1, CELSR2, SCARB1) were significantly associated with CHD or CAD, as also reported in previous studies.14,16,18,20,28 However, the SNPs in the PLA2G7 gene that were most significantly associated with Lp-PLA2 mass (rs1805017) or activity (rs7756935) were not associated with prevalent CHD or CAD. These findings are consistent with Casas et al.,27 who found a significant association between Lp-PLA2 activity and CHD but did not find any significant association with PLA2G7 genotype (12 tagged SNPs including rs1805017). We had >97% power to detect a 5% increase in risk with these two SNPs, but only 35–40% power to detect a 2.5% increase in risk (Supplementary material online, Table S6); therefore, if there is an increased risk, it is most likely in the range of a few percent or less. However, considering the modest effect of these common alleles of the PLA2G7 locus on Lp-PLA2 mass and activity (Supplementary material online, Table S5) and the modest strength of the association of Lp-PLA2 mass and activity with CHD,3 we may have had insufficient power to detect a genetic association that would still be compatible with a causal role for Lp-PLA2 in CHD, despite our large sample size. In contrast to our findings and those of Casas et al., Sutton et al.29 reported significant associations between SNPs in PLA2G7, including rs1805017, with CAD in case–control and family data sets. This suggests that the nature of the population, or in fact having familial data, may influence the strength and likelihood of observing a significant association.

The Val279Phe substitution is located within the catalytical domain of the encoded enzyme and therefore leads to reduction (heterozygotes) or complete loss (homozygotes) of enzyme activity. In Japan, this mutation occurs in 30% of the population (27% heterozygous, 4% homozygous and therefore completely lacking plasma Lp-PLA2 activity and mass).30 Although Val279Phe-mediated loss of activity was an independent CHD risk factor in Japanese men with hypercholesterolaemia, myocardial infarction, stroke, or non-familial dilated and hypertrophic cardiomyopathies and protective from CAD in South Korean men,31 no association with CHD was observed in a Chinese population.32

Clinical implications

In the present meta-analysis, we identified eight SNPs associated with Lp-PLA2 mass or activity which might contribute to the regulation of plasma levels and give insight as to whether Lp-PLA2 is a risk factor or risk marker for CHD. Most of these genetic loci were in regions that have already been associated with levels of lipoproteins. Of the six top SNPs that could be tested in the CARDIoGRAM consortium, four located in ZNF259, APOC1, CELSR2, and SCARB1 were significantly associated with CHD, whereas the two SNPs in PLA2G7, the gene that encodes for Lp-PLA2, did not have any significant association with CHD. This supports the hypothesis that the genetic regions which are known to affect lipoprotein levels also influence both the level of Lp-PLA2 and the development of CHD and is consistent with the hypothesis that Lp-PLA2 level may be a marker of atherogenic lipoproteins; levels of Lp-PLA2 have been shown to be increased in subpopulations of lipoproteins such as electronegative LDL. Although the lack of association does not provide any evidence in support of the hypothesis that Lp-PLA2 is a risk factor for CHD, our data cannot rule out the possibility that Lp-PLA2 is a risk factor. First, despite the large number of CHD cases in CARDIoGRAM, we lacked sufficient power to rule out a small effect. The SNPs identified in PLA2G7 had modest effects on the levels of Lp-PLA2 mass or activity and would be expected to have even smaller effects on risk for development of CHD. Furthermore, we measured total Lp-PLA2 mass and activity in plasma, not the mass or activity associated with LDL or HDL. It is possible that genes which alter lipoprotein metabolism may also alter the binding of Lp-PLA2 to subspecies of lipoproteins,17 and Lp-PLA2 may only be atherogenic when it is associated with atherogenic lipoproteins such as LDL, not with HDL. Our observation that a SNP in the CETP locus, which is associated with lower CETP activity and higher levels of Lp-PLA2 mass, may have implications for clinical programmes that are developing drugs that inhibit CETP.33

Levels of Lp-PLA2 mass and activity were strongly associated with CHD in a large meta-analysis, supporting the view that circulating Lp-PLA2 indeed represents a biomarker of CHD. Experimental data in the pig model indicate that Lp-PLA2 has important proinflammatory effects in the vessel wall.34 Further, clinical trial data using virtual histology (an intravascular ultrasound–derived modality) as a secondary endpoint suggest that elevated level or increased activity of Lp-PLA2 is related to the progression of atherosclerotic disease.35 With regard to a potential causal role in the pathophysiology of atherosclerosis and its complications, only large randomized clinical trials assessing the effect of Lp-PLA2 inhibition on cardiovascular endpoints (STABILITY36 and SOLID–TIMI 5237) can provide a definitive answer.

Strengths and limitations

The routine ascertainment of Lp-PLA2 covariates and GWA data in five community-based cohorts based on >13 000 European ancestry subjects are strengths of our study. However, caution should be taken when generalizing these findings to populations with non-European ancestry. In addition, even the most significant SNPs may be in LD with as yet unknown causal variants, and the functional basis of the relation of the identified SNPs with variations in Lp-PLA2 concentrations requires study. Also, it is possible that missense genetic changes may cause a difference in assay immunoreactivity, warranting caution in interpreting genetic effects on levels of Lp-PLA2. Lastly, because many of the cohorts were community based, CHD was largely symptomatic CAD; we acknowledge that asymptomatic CAD may have been misclassified.

Conclusions

We extended the previously reported Framingham findings6 by conducting a meta-analysis with four additional cohorts. Whereas levels of Lp-PLA2 mass were primarily associated with the gene coding for the protein itself, Lp-PLA2 activity was associated with genetic variants related to LDL-C levels, which were also associated with CHD/CAD. Larger association studies and clinical trials will be needed to determine whether Lp-PLA2 is a causal risk factor for CHD.

Funding

The Atherosclerosis Risk in Communities (ARIC) Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, N01-HC-55022, R01HL087641, R01HL59367, and R01HL086694; National Human Genome Research Institute contract U01HG004402; and National Institutes of Health contract HSN268200625226C. Infrastructure was partly supported by grant number UL1RR025005, a component of the National Institutes of Health and NIH Roadmap for Medical Research. The Cardiovascular Health Study research reported in this article was supported by contract numbers N01-HC-85079 through N01-HC-85086, N01-HC-35129, N01 HC-15103, N01 HC-55222, N01-HC-75150, and N01-HC-45133 and grant numbers U01 HL080295 and R01 HL087652 from the National Heart, Lung, and Blood Institute, with additional contribution from the National Institute of Neurological Disorders and Stroke. A full list of principal CHS investigators and institutions can be found at http://www.chs-nhlbi.org/pi.htm. DNA handling and genotyping was supported in part by the Cedars-Sinai Board of Governors’ Chair in Medical Genetics (J.L.R.), the National Center for Research Resources grant M01RR00425 to the Cedars-Sinai General Clinical Research Center Genotyping core, and the National Institute of Diabetes and Digestive and Kidney Diseases grant DK063491 to the Southern California Diabetes Endocrinology Research Center. Lp-PLA2 measurements were funded by a grant from GlaxoSmithKline. The Framingham Heart Study was supported by the National Heart, Lung, and Blood Institute's Framingham Heart Study (contract number N01-HC-25195) and its contract with Affymetrix, Inc., for genotyping services (contract number N02-HL-6-4278). A portion of this research utilized the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center. The analyses reflect intellectual input and resource development from the Framingham Heart Study investigators participating in the SNP Health Association Resource (SHARe) project; grants HL64753, HL076784, and AG028321 (E.J.B.); and Deutsche Forschungsgemeinschaft (German Research Foundation) Research Fellowship SCHN 1149/1-1 (R.B.S.). Lp-PLA2 activity measurements were provided by GlaxoSmithKline and mass measurements by diaDexus at no cost to the FHS. The Rotterdam Study is supported by the Erasmus Medical Center and Erasmus University Rotterdam; the Netherlands Organization for Scientific Research; the Netherlands Organization for Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly; The Netherlands Heart Foundation; the Ministry of Education, Culture and Science; the Ministry of Health Welfare and Sports; the European Commission; and the Municipality of Rotterdam. Support for genotyping was provided by the Netherlands Organization for Scientific Research (NWO) (175.010.2005.011, 911.03.012) and Research Institute for Diseases in the Elderly (RIDE2) (014-93-015). This study was supported by the Netherlands Genomics Initiative (NGI)/NWO project number 050-060-810. Abbas Dehghan is supported by a grant from NWO (Vici, 918-76-619). The MONICA/KORA Augsburg studies were financed by the Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany, and supported by grants from the German Federal Ministry of Education and Research (BMBF). Part of this work was financed by the German National Genome Research Network (NGFNPlus, project number 01GS0834) and through additional funds from the University of Ulm. Furthermore, the research was supported within the Munich Center of Health Sciences (MC Health) as part of LMU innovative.

Conflict of interest: Lp-PLA2 activity measurements were provided by GlaxoSmithKline and mass measurements by diaDexus at no cost to the FHS. R.R. has been consultant to Cumberland Pharmaceuticals and Celera. R.C.H. has received research/grant support from diaDexus for Lp-PLA2 work in ARIC. J.J.N. is an employee of GlaxoSmithKline. H.S. has received research grants from the EU, project Cardiogenics; NGFN, project Atherogenomics; and CADnet BMBF. H.A.S.-F. is an employee of GlaxoSmithKline. B.M.P. serves on a data and safety monitoring board for a device clinical trial funded by Zoll Lifecor. M.C. has been consultant for GlaxoSmithKline in relation to Lp-PLA2 and had grant funding from them for same. R.P.T.'s  laboratory has received support from GlaxoSmithKline in the form of free assay kits for Lp-PLA2 measurements. W.K. has received honoraria for lectures from diaDexus and GlaxoSmithKline, and is a member of the Steering Committee of the GlaxoSmithKline-sponsored STABILITY and SOLID trials. C.M.B. has received research/grant support from GlaxoSmithKline and diaDexus, and consulting honoraria from GlaxoSmithKline. All other authors have no conflict of interest to declare.

Acknowledgments

The authors thank the staff and participants of the individual studies for their important contributions, and acknowledge the essential role of the Cohorts for Heart and Aging Research in Genome Epidemiology (CHARGE) Consortium in the development and support of this manuscript. CHARGE members include NHLBI's ARIC Study, Cardiovascular Health Study, and Framingham Heart Study; NIA's Iceland Age, Gene/Environment Susceptibility (AGES) Study; and the Netherland's Rotterdam Study.

Footnotes

  • †,‡ These authors contributed equally.

References

View Abstract