Leukotrienes are potent inflammatory and lipid mediators that participate in atherosclerosis. We analyzed the association of Leukotriene gene (ALOX5, ALOX5AP, LTA4H, and LTC4S) polymorphisms and plasma Leukotriene B4 (LTB4) levels with coronary artery disease (CAD) in a representative cohort of Asian Indians. In all, 136 functional single nucleotide polymorphisms (SNPs) were selected using in silico tools. Forty-five polymorphic SNPs were ranked for predicted functional effect using FastSNP. Finally, 14 functional SNPs along with 10 SNPs identified from the literature were genotyped in 340 CAD patients and 340 controls. Plasma LTB4 levels were measured in 150 cases and 150 controls. None of the 24 SNPs showed significant association with CAD. Plasma LTB4 levels were higher in cases than in controls (76.42 ± 4.46 pg/mL versus 60.89 ± 2.61 pg/mL) (P=0.003), with greater risk being associated with the top quartile as compared to the bottom quartile after adjusting for potential confounders (OR 8.94, 95% CI 2.56–31.95; P=0.001). Four SNPs in the LTA4H gene showed significant association with LTB4 levels (P<0.05) of which rs1978331 (P=0.035) remained significant after correction for multiple testing. LTB4 showed strong correlation with lipids (r=0.24–34) only in cases. Our pilot study suggests that the association between Leukotrienes gene polymorphisms and CAD risk may be modulated through plasma LTB4 levels.
1. Introduction
Both genetics and environmental factors combined with unhealthy lifestyle play an important role in the development of coronary artery disease (CAD) [1]. Hypercholesterolemia and inflammation are two key players in CAD [2]. Leukotrienes act as mediators in the lipid and inflammatory pathway and are synthesized from Arachidonic acid through a cascade process in immune-competent cells. During this process, 5-lipoxygenase (ALOX5) in association with 5-LO activating protein (ALOX5AP), metabolizes Arachidonic acid to Leukotriene A4 (LTA4). LTA4 is unstable and is further converted to Leukotriene B4 (LTB4) and Cysteinyl Leukotrienes (LTC4, LTD4, and LTE4) by LTA4 Hydrolase (LTA4H) and LTC4 Synthase (LTC4S) enzymes, respectively [3]. These Leukotrienes then modulate immune response by triggering the synthesis and release of cytokines and interacting with their cognate receptors [4].
Leukotrienes are known to play a significant role in CAD development, based on biochemical, genetic, and pharmacological studies [5–7]. Immunohistochemical studies demonstrate abundant expression of Leukotrienes on the arterial walls [8, 9], which correlate with symptoms of plaque instability [10]. Population studies have also explored the association of Leukotriene gene polymorphisms with cardiovascular disease. The deCODE Genetics group identified two distinct haplotypes, Hap A and Hap B in the ALOX5AP gene through linkage study, that were associated with elevated risk of myocardial infarction (MI) and stroke in two independent cohorts. Further, functional studies provided evidence that these variants induce excess production of LTB4 [11, 12]. The same group also reported modest increase in the risk of MI in carriers of haplotype, Hap K in the LTA4H gene [13]. In addition, there are reports of ALOX5 promoter variants containing multiple Sp1 binding sites that are associated with higher risk of MI [14, 15] and promoter variant rs730012 of LTC4S gene which is associated with greater risk of early carotid atherosclerosis [16].
Despite these significant findings, many other studies do not concur with these reports [17–19]. Such inconsistent reports on different ethnic populations tend to create ambiguity on the direct role of Leukotriene genetic variants in the pathogenesis of atherosclerosis. A noteworthy aspect is that most studies have been carried out on Caucasians while, to our knowledge, there are no studies relating to Leukotrienes and CAD on Asian Indians. Given the high rate of morbidity and mortality compounded by premature disease onset in Asian Indians, there is an urgent need to identify novel genes and other predisposing factors that contribute to higher disease prevalence in this population. In addition, the most frequently studied haplotype, Hap A studied by Helgadottir et al., included SNPs located outside the region of the gene while Hap B was constructed using SNPs without known functional effects. This triggered our interest to screen for functional variants in these genes. The aim of this pilot case-control study was therefore to evaluate the genetic contribution of Leukotriene variants with CAD and to test the association of these variants with plasma LTB4 levels in a representative cohort of Asian Indians.
2. Materials and Methods2.1. Study Subjects
A total of 680 participants comprising of 340 angiographically confirmed CAD patients (cases) and 340 age and gender matched controls were selected from the ongoing Indian Atherosclerosis Research Study (IARS). An overview of the IARS study design has been previously described [20]. Briefly, it is an ongoing epidemiological study with an objective to investigate the genetic factors and biomarkers against a backdrop of conventional cardiovascular risk factors in Asian Indians living in India. Recruitment of cases and controls was based on predefined inclusion/exclusion criteria. Inclusion criteria for cases were as follows: angiographically confirmed presence of CAD with >70% stenosis in any major epicardial artery, diagnosis or >50% stenosis in two or more arteries, patients undergoing percutaneous transluminal coronary angioplasty or bypass graft surgery, and age at onset ≤60 years for men and ≤65 years for women. Control subjects were clinically asymptomatic for CAD, enrolled from the same geographical area as that of the proband, did not have a family history of cardiovascular disease and showed normal electrocardiogram readings. None of the participants had any concomitant infection. Detailed demographics, anthropometrics, vital parameters, medical history, medication, and pedigree information were recorded for each participant through personal interview. Prevalence of type 2 diabetes, hypertension, and cardiovascular disease was ascertained based on self-report of physician’s diagnosis and/or use of prescription medications along with perusal of their medical records. Body mass index was calculated as a ratio of weight in kg to height in meter2. All participants provided signed, voluntary informed consent and the study was approved by the Ethics Committee of TRI and follows the guidelines of the Indian Council of Medical Research (ICMR) on bioethics [21].
2.2. Clinical Laboratory Measurements
Plasma levels of total cholesterol (TC), triglyceride (TG), and HDL cholesterol were measured using reagents from Siemens Dimension Flex reagent cartridge (Siemens Healthcare Diagnostics Ltd, UK) and standards from Randox laboratories (Crumlin, UK) using Siemens dimension Xpand plus instrument (Siemens Dade Behring, Liederbach, Germany). Liquicheck lipid control was used as an assayed quality control (Bio-Rad) that was run along with every experiment. Low-density lipoprotein-cholesterol (LDL-c) was calculated using the Friedwald’s formula [22].
2.3. SNP Selection Process
Figure 1 provides a schematic overview of the SNP selection process. We selected a total of 136 functional variants belonging to four genes, ALOX5, ALOX5AP, LTA4H, and LTC4S from the Leukotriene pathway using Genomatix:SNP analysis (http://www.genomatix.de) and SNPinfo in silico tools, which pick up SNPs that are more likely to affect the biological function [23]. SNPs located in the coding region, transcription-factor-binding sites (TFBS), splice sites, microRNA-binding sites, and untranslated regions (UTRs) were selected using SNPinfo, while Genomatix was used to select SNPs affecting the TFBS and protein sequence. SNPs generated by both tools were combined to compile the final list of 136 functional SNPs (Table S1) (see Supplementary Materials available online at http://dx.doi.org/10.1155/2013/985743). These SNPs were further short listed based on their polymorphic status and risk scores using FastSNP [24], which resulted in the final selection of 14 functional SNPs with low-medium risk scores for genotyping and located in the 5′ or 3′ region (Table 1). The FastSNP prioritization result is based on predicted functional effects and their estimated risk.
Final list of functional SNPs (N = 14) selected for genotyping based on the predicted functional effect and estimated risk score.
Gene
rs ID
Position (bp)
Type of SNP
Fast SNP rank
Literature Tagging SNP
ALOX5
rs2291427
45256230
intronic
(1-2)
rs12762303
45189177
5′ upstream
(1–3)
rs1864414
45190567
Intronic
(1-2)
rs41314980
45261505
3′ UTR
(0-0)
ALOX5AP
rs4769055
30207830
intronic
(1-2)
rs10507391 (Hap A/B)
rs17222919
30206329
5′ upstream
(0-0)
rs4360791
30216020
intronic
(0-0)
rs9551963 (Hap A)
rs3803277
30216308
intronic
(0-0)
rs9551963 (Hap A)
rs1132340
30236462
3′ UTR
(0-0)
rs4769874 (Hap A)
LTA4H
rs2660899
94954551
5′upstream
(1–3)
rs17525495
94953508
5′ upstream
(0-0)
rs2660898
94950128
intronic
(1-2)
Hap K
rs1978331
94933332
intronic
(1-2)
Hap K, rs2540482, rs2660845
LTC4S
rs730012
179153244
5′ upstream
(0-0)
1-2: very low to low; 1–3: very low to medium; 0: unknown.
Flow chart showing the SNP selection process.
With regard to the literature-based SNPs, we selected Hap A and Hap B blocks of ALOX5AP and Hap K block of LTA4H [12, 13], and the haplotypes were defined by the following SNPs, Hap A: rs17222814 (G), rs10507391 (T), rs4769874 (G), rs9551963 (A); Hap B: rs17216473 (A), rs10507391 (A), rs9315050 (A), rs17222842 (G); Hap K: rs61937881 (C), rs2660880 (G), rs6538697 (T), rs1978331 (A), rs17677715 (T), rs2247570 (T), rs2660898 (T), rs2540482 (C), rs2660845 (G), and rs2540475 (G). To minimize the number of SNPs required to be genotyped, tagging SNPs were selected based on the linkage disequilibrium (LD) distribution pattern between the functional and literature based SNPs. Finally, 10 additional SNPs from literature were selected from ALOX5AP (Hap A-rs17222814, Hap B-rs17216473, rs9315050, and rs17222842) and LTA4H (rs2540475, rs2247570, rs17677715, rs6538697, rs2660880, and rs61937881) genes, respectively.
2.4. SNP Genotyping
Genomic DNA was extracted using a modified salting out procedure [25] and quantitated on ND1000 spectrophotometer (NanoDrop Technologies Inc., Wilmington, DE, USA). Primers for the 14 SNPs were designed using RExPrimer, an integrated graphical web-based method (http://www4a.biotec.or.th/rexprimer), and the oligonucleotides were purchased from Sigma-Aldrich (St. Louis, USA). Forty-six primer pairs were designed to annotate for 136 SNPs (Table S2) and their polymorphic status was established on pooled DNA samples [26], based on targeted sequencing. Genomic DNA from 8 randomly selected individuals was pooled in equal volumes at 20 ng/μL final concentration to construct a single DNA pool. In total, 6 DNA pools for cases and controls were constructed separately from 48 cases and 48 controls, respectively. Regions of interest were amplified on a 9700 GeneAmp PCR instrument and the products were sequenced bi-directionally using Big Dye Terminator v3·1 sequencing kit from ABI as per manufacturer’s recommendations and analyzed with 3130XL automated genetic analyzer (Applied Biosystems, Foster City, USA) using the SeqScape v 2.5 software (Applied Biosystems, Foster City, USA).
After establishing the polymorphic status of SNPs followed by ranking, 14 functional SNPs were genotyped in 680 subjects using allele specific melting temperature Tm shift method [27]. The two allele specific forward primers and one common reverse primer were designed using Primer Express V2.0 (Applied Biosystems, Foster City, CA, USA). The primer sequence is shown in Table S3. Genotyping of 10 literature-based SNPs for 680 subjects was performed by TaqMan probe-based assay. Tm-shift genotyping and allelic discrimination analysis was performed on ABI 7900 HT Sequence Detection System using SDSv2.1 software (Applied Biosystems, Foster City, CA, USA). Ten percent of samples were analyzed in duplicates to confirm genotype accuracy and further verified by sequencing. We obtained 99.5% accuracy in the genotype call rates.
2.5. LTB4 Immunoassay
LTB4 immunoassay was performed in EDTA plasma samples of 150 cases and 150 matched controls using ELISA kit (Cayman Chemical Company, Ann Arbor, MI, USA) according to the manufacturer’s instructions. Absorbance was read at 412 nm with Bio-Tek power waveX microplate reader (Bio-Tek instrument, Vermont, USA). The intraassay and interassay variability were less than 10%, with 100% specificity, and 13 pg/mL assay detection limits.
2.6. Statistical Analysis
All variants were tested for Hardy-Weinberg equilibrium (HWE) in the control group. SNPs deviating from HWE (P<0.05) were excluded from further analysis. HWE, odds ratio (ORs) at 95% confidence intervals (CIs), haplotype frequency, and corresponding P values for each SNP were estimated separately for cases and controls using SNPstat tool (http://bioinfo.iconcologia.net/index.php?module=Snpstats). Student’s t-test, ANOVA, and univariate analysis were used to compare the means for continuous variable while χ² test was used to estimate the association of categorical variables, namely, genotypes with CAD. Normality distribution of plasma LTB4 levels was assessed using Q-Q plot and outliers were removed from further analysis. The final data was log transformed before performing logistic regression to estimate the association between LTB4 quartiles and CAD. Age, gender, hypertension, diabetes, TC, and HDL-c were used as covariates and appropriately adjusted for during the analysis. An interactive linkage disequilibrium (LD) block between SNPs was created using Haploview version 4.2 (http://www.broadinstitute.org/haploview/haploview). Routine analysis was performed using SPSS v17.0 (SPSS Inc., IL, USA). Data were represented as mean ± standard error of mean, unless otherwise indicated. A nominal P value of ≤ 0.05 less was regarded as statistically significant.
3. Results3.1. Clinical Characteristics of Study Subjects
The clinical profile of the study cohort is summarized in Table 2. Out of 340 cases, 225 (66.2%) had positive family history of CVD with at least one affected first degree relative, while 115 subjects (33.8%) did not have a family history of CVD. Average age was 46.93±6.68 years for cases and 47.39±6.66 years for controls, with age and gender being comparable across the two groups. Nearly 50.3% and 40.8% of cases were hypertensive and diabetic, respectively, which was significantly higher than the corresponding prevalence among controls (14.4% and 15.9%, resp.) (P<0.001). There were 179 (54.2%) myocardial infarction (MI) cases and 147 (44.5%) non-MI cases. Frequency distribution of subjects based on the number of disease vessels was as follows: 0 or 1 vessel-59 (18.4%), 2 vessels-83 (25.95%), and 3 or more vessels-179 (55.8%). TC and LDL-cholesterol levels were lower in CAD patients than in controls, which might be attributed to the higher usage of lipid lowering drugs in this group.
Clinical characteristics of the study population.
Controls (n = 340)
Cases (n = 340)
P value
Age (in years)
47.39 ± 6.66
46.93 ± 6.68
0.368
Male n (%)/female n (%)
250 (73.5)/90 (26.5)
246 (73)/91 (27)
BMI (kg/m2)
25.74 ± 4.74
25.97 ± 4.48
0.517
Total cholesterol (mg/dL)
172.74 ± 39.71
139.83 ± 39.24
<0.0001
HDL-cholesterol (mg/dL)
41.44 ± 9.78
34.8 ± 9.95
<0.0001
LDL-cholesterol (mg/dL)
99.88 ± 31.05
72.52 ± 32.45
<0.0001
Triglycerides (mg/dL)
162.11 ± 112.47
166.51 ± 83.07
<0.0001
Hypertension n (%)
49 (14.4)
169 (50.3)
<0.0001
Diabetes n (%)
54 (15.9)
137 (40.8)
<0.0001
SBP (mmHg)
124.2 ± 16.66
119.9 ± 15.28
<0.0001
DBP (mmHg)
81.63 ± 10.75
78.15 ± 9.04
<0.0001
Age at onset of CAD
…
45.77 ± 7.03
…
Positive family history
…
225 (66.2)
…
Negative family history
…
115 (33.8)
…
Stable angina
…
147 (44.5)
…
Unstable angina
…
179 (54.2)
…
Number of disease vessels n (%)
1
59 (18.4)
…
2
83 (25.9)
…
3
179 (55.8)
…
Continuous variables are expressed as mean ± standard deviation. BMI: body mass index; HDL-c: high density lipoprotein cholesterol; LDL-c: low density lipoprotein cholesterol; SBP: systolic blood pressure; DBP: diastolic blood pressure.
3.2. Assessment of Polymorphic Status in the Leukotriene Genes
We evaluated the polymorphic status of 136 functional SNPs by sequencing pooled DNA samples. While 91 markers were found to be invariant, the remaining 45 polymorphic SNPs were ranked to obtain a list of 14 functional SNPs with low to medium risk. 10 SNPs from literature were included and a total of 24 SNPs were genotyped in 680 subjects.
3.3. Association of Leukotriene SNPs with CAD
Except rs2660880, all other SNPs were in HWE in the controls (P>0.05). There was no significant difference in allele or genotype frequencies for any of the SNPs between cases and controls (Table 3). Also, there was no association between the SNPs and disease severity, based on the number of disease vessels. A SNP in the ALOX5 gene, rs12762303, showed significant association with CAD only in males, with the minor allele, “C” showing a protective effect (OR 0.31, 95% CI=0.12–0.81; P=0.007).
Distribution of genotypes of the short-listed SNP panel in cases and controls and estimated association with risk of CAD.
SNPs
Genotype
Frequency (number of samples)
OR (95% CI)
P value
Control (n%)
Cases (n%)
ALOX5
TT
210 (61.8%)
211 (62.4%)
rs12762303
TC
110 (32.4%)
113 (33.4%)
1.01 (0.73–1.40)
0.58
CC
20 (5.9%)
14 (4.1%)
0.70 (0.34–1.42)
GG
118 (34.7%)
131 (38.8%)
rs2291427
GA
169 (49.7%)
155 (45.9%)
0.83 (0.60–1.16)
0.53
AA
53 (15.6%)
52 (15.4%)
0.89 (0.56–1.41)
AA
212 (62.4%)
217 (64.2%)
rs1864414
AC
112 (32.9%)
109 (32.2%)
0.94 (0.68–1.30)
0.72
CC
16 (4.7%)
12 (3.5%)
0.73 (0.34–1.59)
CC
321 (94.4%)
318 (94.4%)
rs41314980
CA
19 (5.6%)
19 (5.6%)
1.01 (0.52–1.94)
0.98
AA
…
…
ALOX5AP
AA
257 (75.6%)
258 (76.8%)
rs9315050
AG
78 (22.9%)
73 (21.7%)
0.93 (0.65–1.34)
0.93
GG
5 (1.5%)
5 (1.5%)
1.00 (0.28–3.48)
GG
212 (62.4%)
208 (61.7%)
rs17216473
GA
115 (33.8%)
107 (31.8%)
0.95 (0.68–1.31)
0.27
AA
13 (3.8%)
22 (6.5%)
1.72 (0.85–3.51)
GG
280 (82.3%)
282 (83.7%)
rs17222814
GA
60 (17.6%)
52 (15.4%)
0.86 (0.57–1.29)
0.09
AA
0 (0%)
3 (0.9%)
…
GG
295 (86.8%)
307 (91.1%)
rs17222842
GA
43 (12.7%)
30 (8.9%)
0.67 (0.41–1.10)
0.07
AA
2 (0.6%)
0 (0%)
…
TT
207 (60.9%)
204 (60.4%)
rs17222919
TG
116 (34.1%)
116 (34.3%)
1.02 (0.74–1.41)
0.98
GG
17 (5%)
18 (5.3%)
1.08 (0.54–2.15)
CC
95 (27.9%)
85 (25.3%)
rs4769055
CA
163 (47.9%)
157 (46.7%)
1.10 (0.77–1.59)
0.48
AA
95 (27.9%)
85 (25.3%)
1.30 (0.85–1.97)
AA
109 (32.1%)
110 (32.5%)
rs4360791
AG
182 (53.5%)
168 (49.7%)
0.91 (0.65–1.28)
0.48
GG
49 (14.4%)
59 (17.5%)
1.19 (0.75–1.89)
AA
109 (32.1%)
117 (34.6%)
rs3803277
AC
181 (53.2%)
162 (47.9%)
0.83 (0.60–1.17)
0.38
CC
50 (14.7%)
58 (17.2%)
1.08 (0.68–1.71)
AA
261 (76.8%)
261 (77.5%)
rs1132340
AG
75 (22.1%)
72 (21.4%)
0.96 (0.67–1.39)
0.98
GG
4 (1.2%)
4 (1.2%)
1.00 (0.25–4.06)
LTA4H
TT
280 (82.6%)
281 (84.1%)
rs17677715
TC
55 (16.2%)
51 (15.3%)
0.92 (0.61–1.40)
0.67
CC
4 (1.2%)
2 (0.6%)
0.50 (0.09–2.74)
AA
169 (49.7%)
175 (52.1%)
rs2540475
GA
140 (41.2%)
135 (40.2%)
0.93 (0.68–1.28)
0.74
GG
31 (9.1%)
26 (7.7%)
0.81 (0.46–1.42)
AA
233 (68.5%)
244 (72.4%)
rs2247570
AG
97 (28.5%)
87 (25.8%)
0.86 (0.61–1.20)
0.41
GG
10 (2.9%)
6 (1.8%)
0.57 (0.21–1.60)
TT
205 (60.3%)
204 (60.5%)
rs6538697
TC
115 (33.8%)
117 (34.7%)
1.02 (0.74–1.41)
0.8
CC
20 (5.9%)
16 (4.8%)
0.80 (0.41–1.60)
GG
263 (77.3%)
263 (78%)
rs61937881
GA
71 (20.9%)
71 (21.1%)
1.00 (0.69–1.45)
0.60
AA
6 (1.8%)
3 (0.9%)
0.50 (0.12–2.02)
TT
140 (41.2%)
140 (41.5%)
rs2660898
TG
155 (45.6%)
160 (47.3%)
1.03 (0.75–1.42)
0.65
GG
45 (13.2%)
37 (11%)
0.82 (0.50–1.35)
GG
189 (55.6%)
183 (54.1%)
rs2660899
GT
123 (36.2%)
126 (37.3%)
1.06 (0.77–1.47)
0.92
TT
28 (8.2%)
29 (8.6%)
1.08 (0.62–1.88)
CC
204 (60.2%)
202 (59.8%)
rs17525495
CT
116 (34.1%)
119 (35.2%)
1.04 (0.75–1.43)
0.79
TT
20 (5.9%)
16 (4.8%)
0.81 (0.41–1.60)
TT
96 (28.2%)
101 (30%)
rs1978331
TC
172 (50.6%)
170 (50.5%)
0.94 (0.66–1.33)
0.79
CC
72 (21.2%)
65 (19.4%)
0.86 (0.55–1.33)
LTC4S
AA
234 (68.8%)
221 (65.4%)
rs730012
AC
96 (28.2%)
104 (30.8%)
1.15 (0.82–1.60)
0.65
CC
10 (2.9%)
12 (3.6%)
1.27 (0.54–3.00)
3.4. Haplotype Analysis
Figures S1(a)–S1(c) show the LD distribution among the Leukotriene SNPs of ALOX5, ALOX5AP, and LTA4H genes. There was strong pairwise LD between rs12762303 and rs1864414 SNPs in the ALOX5 gene (r2=0.99). rs4360791 and rs3803277 as well as rs9315050 and rs1132340 in ALOX5AP gene formed two distinct haplotype blocks. rs17222814 appeared to be the lead SNP in the ALOX5AP gene as it was in strong LD with the other SNPs in this gene (r2=0.90–0.96). There were two distinct haplotype blocks in the LTA4H gene, which included block 1 (5 SNPs) and block 2 (3 SNPs). Modest pair-wise LD was observed between rs2540475 and other SNPs in this gene (r2=0.41–0.73).
Distribution of haplotype frequencies is shown in Table S4. The TGAC haplotype in ALOX5 gene defined by 4 SNPs, rs12762303 rs2291427, rs1864414, rs41314980 and AGGGTCGCA haplotype in ALOX5AP gene defined by 9 SNPs, rs9315050, rs17216473, rs17222814, rs17222842, rs17222919, rs4769055 rs4360791, rs3803277, and rs1132340 and the TGATGTGCT haplotype in LTA4H gene constituted by 9 SNPs, rs17677715, rs2540475, rs2247570, rs6538697, rs61937881, rs2660898, rs2660899, rs17525495, and rs1978331 showed a frequency of 0.42, 0.26 and 0.32 in cases, respectively and were the best represented haplotypes in our cohort. However, none of the haplotypes showed significant association with CAD.
We also estimated the haplotype frequency of Hap A and Hap B of the ALOX5AP gene and compared them with the frequencies reported in the literature. Hap A comprised of rs172228149 (G), rs4769055(C), rs1132340(A), rs4360791(A) and Hap B comprised of rs17216473(A), rs4769055(A), and rs9315050(A), and rs17222842(G). Neither Hap A nor Hap B showed significant association with CAD. The distribution frequency is shown in Table S5.
3.5. Association of Plasma LTB4 Levels with CAD
Plasma LTB4 levels were measured in 150 cases and 150 matched controls. After removing the outliers, data on 141 cases and 144 controls was analyzed. LTB4 levels were 20% higher in cases as compared to controls (76.42±4.46 pg/mL versus 60.89±2.61 pg/mL; P=0.003) (Figure 2). Four SNPs rs6538697, rs2660898, rs17525495, and rs1978331 in the LTA4H gene showed significant association with LTB4 levels (Figure 3(a)), with the minor allele being associated with lower LTB4 levels (Figure 3(b)). However, only one SNP, rs1978331 (P=0.035), retained significance following Bonferroni adjustment for multiple comparisons.
Distribution of mean LTB4 levels between affected and unaffected subjects. Data expressed as mean±standard deviation (SD).
(a) Association results of all 24 SNPs with LTB4 are displayed. Genomic position (x-axis) and negative log10 of P values (y-axis) are noted. SNPs meeting statistical significance have been annotated. (b) Bar diagram showing distribution of LTB4 in significant SNPs (rs6538697, rs2660898, rs17525495, and rs1978331) with different genotype. Bar diagram shows the distribution of LTB4 levels across the four LTA4H SNPs-rs6538697 (T>C), rs2660898 (T>G), rs17525495 (C>T), rs1978331 (T>C).
Plasma LTB4 levels were distributed into quartiles, Q1 (<37.59 pg/mL), Q2 (37.59–54.68 pg/mL), Q3 (54.69–85.59 pg/mL), and Q4 (>85.59 pg/mL). Frequency of cases was higher in Q4 LTB4 levels (n=43, 30.5%) as compared to the controls (n=28, 19.4%). Subjects in the top LTB4 quartile showed significant higher risk of CAD as compared to the bottom quartile (OR = 8.935, 95% CI=2.559–3.195; P=0.001), after adjusting for potential confounders. Plasma LTB4 levels showed significant correlation with TC and LDL-cholesterol (r>0.30, P<0.0001), after adjusting for age, gender, and statins, particularly in cases but not in controls (Table S6).
4. Discussion
Leukotrienes are potent proinflammatory mediators belonging to the 5-lipoxygenase pathway, involved in Arachidonic acid metabolism. ALOX5, ALOX5AP, LTA4H, and LTC4S genes and their metabolic products have been investigated with regard to their role in the pathogenesis of cardiovascular disease [28]. In this pilot study, we examined the association of a select panel of Leukotriene SNPs encoded by these genes. In addition, we measured plasma LTB4 levels and tested its association with Leukotriene variants and CAD. Leukotriene SNPs did not show direct association with CAD but significant association was noted between plasma LTB4 levels and CAD as well as with four select SNPs in the LTA4H gene.
The strength of our study lies in the use of a robust SNP selection process. Most of the genetic association studies carried out on Leukotriene gene variants have focused on markers without functional properties [11–13]. We have adopted a systematic process of selecting functional variants using in silico approach followed by targeted sequencing to establish their polymorphic status. In fact, reports suggest that it might be rewarding to study rare rather than common variants in order to understand the etiology of complex diseases [29]. Such rare variants were selected using SNPinfo and Genomatix tools along with the common variants found in the literature. To our knowledge, this is the first study to use a comprehensive in silico approach for SNP prioritization on the basis of functional relevance in the Leukotriene pathway. Of the 136 SNPs investigated, only 45 SNPs (30%) were polymorphic, which could be attributed to a low minor allele frequency (MAF) of the rare functional variants in the Indian population. In fact, the average MAF was considerably higher for SNPs in the noncoding regions (20%–25%) than in the coding regions (12%–16%) [30], making it difficult to pick them up in our rather small-sized cohort. A recent review of GWAS suggests that more than 70% of genetic variants identified to date for CAD are located in the noncoding regions [31].
Family history is considered to be an independent risk factor for CAD [32]. Hence, in order to negate this bias, we included 225 CAD subjects and 115 subjects with and without family history, respectively. None of the SNPs in either group showed significant association with the disease. The genetic contribution of Hap A and Hap B and their individual SNPs to the cardiovascular phenotypes in Whites, Blacks, and East Asian populations have yielded both positive [11, 12], and negative results [17–19, 33]. In the Icelandic population, Hap A showed association with MI, with a case/control frequency of 0.158/0.095 [12], but no association in the Scottish (0.184/0.142) [11] and Caucasian (0.15/0.15) [34] populations. The case-control frequency obtained in our study (0.13/0.15) was comparable to the German study (14.5/15.2) [18]. On the other hand, prevalence of Hap B in case controls was 0.07/0.07 in Iceland whites [12], 0.07/0.06 in Scottish whites [11], 0.09/0.08 in Caucasians [35], and 0.10/0.08 in African Americans [35]. The present study showed a frequency of 0.18/0.15. None of the studies could demonstrate a significant association between Hap B and CAD. On the other hand, positive association was seen only in the Icelandic cohort [12], and this could be attributed to population-specific distribution in the allele/haplotype frequencies. Other factors such as phenotypic heterogeneity, age/gender bias, and, importantly, ethnic differences can also contribute to the genotype-phenotype association. Further, studies have dealt with varied phenotypes such as CAD, Myocardial Infarction (MI), and stroke, which makes comparisons and interpretations difficult.
The role of LTB4 levels in CAD pathogenesis may be mediated through the pro-inflammatory molecules such as IL-6, TNF-α, and MCP-1 [36]. Previous studies have demonstrated increased LTB4 production in stimulated peripheral mononuclear cells in MI patients [12]. We observed higher LTB4 levels in cases than in controls and putative association between plasma LTB4 levels and 4 SNPs in the LTA4H gene. One of the SNPs, rs17525495, is in fact located in the promoter region and may modulate LTB4 production. The other 3 SNPs are located in the intronic region but exhibit strong LD with the promoter SNP which could explain their association with LTB4 levels. rs1978331 showed persistent association even after Bonferroni correction suggesting a robust association of this variant with plasma LTB4 levels. Thus, it appears that the LTA4H genetic variants regulate LTB4 production under the catalytic influence of the LTA4H enzyme. Evidence suggests that the rare T allele of rs17525495 increases LTA4H gene transcription [37], while the rare alleles of rs6538697, rs2660898, and rs1978331 are associated with lower LTB4 levels [13], which is comparable to the present study findings. Some studies have shown that LTB4 biosynthesis is regulated by cholesterol [38], and inhibition of Leukotriene may reduce plasma TG levels [39]. Our data also showed high correlation of LTB4 levels with TC, TG, and LDL-cholesterol, after adjusting for age, gender, and statins. Interestingly, this correlation was recorded in cases but not in controls reflecting the presence of an atherothrombotic state in CAD.
The short comings of the study are small sample size that might not be able to capture the true magnitude of the modest effect sizes of the genetic variants on CAD. Nonetheless, the interesting findings in this pilot study do provide a basis for planning larger studies on the Leukotriene pathway. Additionally, there may be sampling bias since the controls were selected based on their self-report of absence of any clinical symptoms and normal electrocardiogram readings which does not preclude the presence of a silent underlying disease. It has been demonstrated that 80% statistical power can be obtained under the assumption of 5% prevalence, complete LD, 5% α level given a sample size of 248 cases, and 248 controls for any particular SNP [40]. Under these circumstances, given a sample size of 340 cases and 340 controls and SNPs showing >5% MAF, our study should have >80% power to detect associations. It is however common knowledge that the power to detect small effect sizes of the alleles is duly enhanced in a larger cohort.
In conclusion, to our knowledge this is the first report of a case-control study designed to assess the association between functional polymorphisms in four key Leukotriene genes and CAD in a representative cohort of Asian Indians. Although there was no significant association between Leukotriene gene variants and disease, our preliminary finding showed higher LTB4 levels in CAD individuals and strong correlation between LTB4 levels and lipids in cases alone, reiterating the occurrence of an atherothrombotic state in these subjects. The association between LTA4H SNPs and plasma LTB4 levels suggests that these functional variants might modulate atherosclerosis by regulating LTB4 production. Future studies aimed at understanding the interactions between Leukotrienes and other key inflammatory biomarkers such as interleukins (IL 4, IL5, IL 6), TNF-α, TGF-β, and IFN-γ can throw better light on the inflammation-driven atherosclerotic disease process.
Acknowledgments
This work was supported by the Thrombosis Research Institute, London, the Department of Biotechnology, Ministry of Science and Technology, Government of India (Grant no. BT/01/CDE/08/07), Tata Social Welfare Trust, India (Grant no. TSWT/IG/SNB/JP/Sdm), the Weston foundation, and UK and Bey Foundation, Switzerland. The authors thank all the investigators, staff, administrative teams and participants of the IARS from Narayana Hrudayalaya, Bangalore, for their valuable contributions. The sponsors did not participate in the design, conduct, sample collection analysis and interpretation of the data or in the preparation, review, or approval of the paper. The authors have no conflicts of interests to disclose.
YusufS.HawkenS.OunpuuS.Effect of potentially modifiable risk factors associated with myocardial infarction in 52 countries (the INTERHEART study): case-control study200436494509379522-s2.0-994424275610.1016/S0140-6736(04)17513-2SteinbergD.Atherogenesis in perspective: hypercholesterolemia and inflammation as partners in crime2002811121112172-s2.0-003685354310.1038/nm1102-1211Peters-GoldenM.HendersonW. R.Jr.Leukotrienes200735718184118542-s2.0-3574896172010.1056/NEJMra071371SadikC. D.LusterA. D.Lipid-cytokine-chemokine cascades orchestrate leukocyte recruitment in inflammation20129122072152-s2.0-8485644462610.1189/jlb.0811402BäckM.HanssonG. K.Leukotriene receptors in atherosclerosis20063874935022-s2.0-3375111328610.1080/07853890600982737MehrabianM.AllayeeH.5-Lipoxygenase and atherosclerosis20031454474572-s2.0-014209315110.1097/00041433-200310000-00005RiccioniG.ZanasiA.VitulanoN.ManciniB.D'OrazioN.Leukotrienes in atherosclerosis: new target insights and future therapy perspectives2009200962-s2.0-7674913276910.1155/2009/737282737282Di GennaroA.WagsáterD.MáyránpááM. I.GabrielsenA.SwedenborgJ.HamstenA.SamuelssonB.ErikssonP.HaeggströmJ. Z.Increased expression of leukotriene C4 synthase and predominant formation of cysteinyl-leukotrienes in human abdominal aortic aneurysm20101074921093210972-s2.0-7865046557210.1073/pnas.1015166107SpanbroekR.GräbnerR.LötzerK.HildnerM.UrbachA.RühlingK.MoosM. P. W.KaiserB.CohnertT. U.WahlersT.ZieskeA.PlenzG.RobenekH.SalbachP.KühnH.RådmarkO.SamuelssonB.HabenichtA. J. R.Expanding expression of the 5-lipoxygenase pathway within the arterial wall during human atherogenesis20031003123812432-s2.0-003741797310.1073/pnas.242716099QiuH.GabrielsenA.AgardhH. E.WanM.WetterholmA.WongC.HedinU.SwedenborgJ.HanssonG. K.SamuelssonB.Paulsson-BerneG.HaeggströmJ. Z.Expression of 5-lipoxygenase and leukotriene A4 hydrolase in human atherosclerotic lesions correlates with symptoms of plaque instability200610321816181662-s2.0-3374449755110.1073/pnas.0602414103HelgadottirA.GretarsdottirS.St. ClairD.ManolescuA.CheungJ.ThorleifssonG.PasdarA.GrantS. F. A.WhalleyL. J.HakonarsonH.ThorsteinsdottirU.KongA.GulcherJ.StefanssonK.MacLeodM. J.Association between the gene encoding 5-lipoxygenase-activating protein and stroke replicated in a Scottish population20057635055092-s2.0-1384428914210.1086/428066HelgadottirA.ManolescuA.ThorleifssonG.GretarsdottirS.JonsdottirH.ThorsteinsdottirU.SamaniN. J.GudmundssonG.GrantS. F. A.ThorgeirssonG.SveinbjornsdottirS.ValdimarssonE. M.MatthiassonS. E.JohannssonH.GudmundsdottirO.GurneyM. E.SainzJ.ThorhallsdottirM.AndresdottirM.FriggeM. L.TopolE. J.KongA.GudnasonV.HakonarsonH.GulcherJ. R.StefanssonK.The gene encoding 5-lipoxygenase activating protein confers risk of myocardial infarction and stroke20043632332392-s2.0-1074422079410.1038/ng1311HelgadottirA.ManolescuA.HelgasonA.ThorleifssonG.ThorsteinsdottirU.GudbjartssonD. F.GretarsdottirS.MagnussonK. P.GudmundssonG.HicksA.JonssonT.GrantS. F. A.SainzJ.O'BrienS. J.SveinbjornsdottirS.ValdimarssonE. M.MatthiassonS. E.LeveyA. I.AbramsonJ. L.ReillyM. P.VaccarinoV.WolfeM. L.GudnasonV.QuyyumiA. A.TopolE. J.RaderD. J.ThorgeirssonG.GulcherJ. R.HakonarsonH.KongA.StefanssonK.A variant of the gene encoding leukotriene A4 hydrolase confers ethnicity-specific risk of myocardial infarction200638168742-s2.0-2944444474810.1038/ng1692AllayeeH.BaylinA.HartialaJ.WijesuriyaH.MehrabianM.LusisA. J.CamposH.Nutrigenetic association of the 5-Lipoxygenase gene with myocardial infarction20088849349402-s2.0-53849103285DwyerJ. H.AllayeeH.DwyerK. M.FanJ.WuH.MarR.LusisA. J.MehrabianM.Arachidonate 5-lipoxygenase promoter genotype, dietary arachidonic acid, and atherosclerosis2004350129372-s2.0-034741827710.1056/NEJMoa025079IovannisciD. M.LammerE. J.SteinerL.ChengS.MahoneyL. T.DavisP. H.LauerR. M.BurnsT. L.Association between a leukotriene C4 synthase gene promoter polymorphism and coronary artery calcium in young women: the Muscatine study20072723943992-s2.0-3384654510810.1161/01.ATV.0000252680.72734.10KochW.HoppmannP.MuellerJ. C.SchömigA.KastratiA.No association of polymorphisms in the gene encoding 5-lipoxygenase-activating protein and myocardial infarction in a large central European population2007931872-s2.0-3424715438410.1097/01.gim.0000263700.53258.b0LõhmussaarE.GschwendtnerA.MuellerJ. C.OrgT.WichmannE.HamannG.MeitingerT.DichgansM.ALOX5AP gene and the PDE4D gene in a central European population of stroke patients20053647317362-s2.0-1684436883410.1161/01.STR.0000157587.59821.87ZeeR. Y. L.ChengS.HegenerH. H.ErlichH. A.RidkerP. M.Genetic variants of arachidonate 5-lipoxygenase-activating protein, and risk of incident myocardial infarction and ischemic stroke: a nested case-control approach2006378200720112-s2.0-3374718398310.1161/01.STR.0000229905.25080.01ShankerJ.MaitraA.RaoV. S.MundkurL.DhanalakshmiB.HebbagodiS.KakkarV. V.Rationale, design and preliminary findings of the Indian atherosclerosis research study20106242862952-s2.0-79952275439KumarN. K.Bioethics activities in India2006121S56S652-s2.0-33749494395FriedewaldW. T.LevyR. I.FredricksonD. S.Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge19721864995022-s2.0-0015348189XuZ.TaylorJ. A.SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies2009372W600W6052-s2.0-6784910118210.1093/nar/gkp290YuanH.ChiouJ.TsengW.LiuC.LiuC.LinY.WangH.YaoA.ChenY.HsuC.FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization200634W635W6412-s2.0-3374783218310.1093/nar/gkl236MillerS. A.DykesD. D.PoleskyH. F.A simple salting out procedure for extracting DNA from human nucleated cells1988163, article 12152-s2.0-002428402810.1093/nar/16.3.1215ShamP.BaderJ. S.CraigI.O'DonovanM.OwenM.DNA pooling: a tool for large-scale association studies20023118628712-s2.0-003684488110.1038/nrg930WangJ.ChuangK.AhluwaliaM.PatelS.UmblasM.MirelD.HiguchiR.GermerS.High-throughput SNP genotyping by single-tube PCR with Tm-shift primers20053968858922-s2.0-3364481459110.2144/000112028HaribabuB.Leukotrienes: novel targets for vascular disease20044281287BodmerW.BonillaC.Common and rare variants in multifactorial susceptibility to common diseases20084066957012-s2.0-4434913270810.1038/ng.f.136GoddardK. A. B.HopkinsP. J.HallJ. M.WitteJ. S.Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms five populations20006612162342-s2.0-003392613010.1086/302727RobertsR.StewartA. F. R.9p21 and the genetic revolution for coronary artery disease20125811041122-s2.0-8485606877310.1373/clinchem.2011.172759SheaS.OttmanR.GabrieliC.Family history as an independent risk factor for coronary artery disease1984447938012-s2.0-0021219582WangG.WangY.SunH.CaoW.ZhangJ.XiaoH.ZhangJ.Variants of the arachidonate 5-lipoxygenase-activating protein (ALOX5AP) gene and risk of ischemic stroke in Han Chinese of eastern China20112553193272-s2.0-8005486987410.1016/S1674-8301(11)60043-2HartialaJ.LiD.ContiD. V.VikmanS.PatelY.Wilson TangW. H.BrennanM.NewmanJ. W.StephensenC. B.ArmstrongP.HazenS. L.AllayeeH.Genetic contribution of the leukotriene pathway to coronary artery disease201112966176272-s2.0-7995975024510.1007/s00439-011-0963-3CrosslinD. R.ShahS. H.NelsonS. C.HaynesC. S.ConnellyJ. J.GadsonS.Goldschmidt-ClermontP. J.VanceJ. M.RoseJ.GrangerC. B.SeoD.GregoryS. G.KrausW. E.HauserE. R.Genetic effects in the leukotriene biosynthesis pathway and association with atherosclerosis200912522172292-s2.0-6054908520810.1007/s00439-008-0619-0Stanke-LabesqueF.PepinJ. L.de JouvencelT.ArnaudC.BaguetJ. P.Leukotriene B4 pathway activation and atherosclerosis in obstructive sleep apnea20125319441951TobinD. M.RocaF. J.OhS. F.McFarlandR.VickeryT. W.RayJ. P.KoD. C.ZouY.BangN. D.ChauT. T. H.VaryJ. C.HawnT. R.DunstanS. J.FarrarJ. J.ThwaitesG. E.KingM.SerhanC. N.RamakrishnanL.Host genotype-specific therapies can optimize the inflammatory response to mycobacterial infections201214834344462-s2.0-8485669317110.1016/j.cell.2011.12.023ZagryagskayaA. N.AleksandrovD. A.PushkarevaM. A.GalkinaS. I.GrishinaZ. V.Sud'InaG. F.Biosynthesis of leukotriene B4 in human polymorphonuclear leukocytes: regulation by cholesterol and other lipids2008543473522-s2.0-6524916655110.1080/15476910802482888ChoiJ.JeonH. J.ParkJ.SonnS. K.LeeM.LeeM.YouH. J.KimG.KimJ.LeeM. H.KwonO.NamK.KimH.JeongT.LeeW. S.OhG. T.Anti-atherogenic effect of BHB-TZD having inhibitory activities on cyclooxygenase and 5-lipoxygenase in hyperlipidemic mice201021211461522-s2.0-7795623231110.1016/j.atherosclerosis.2010.05.003HongE. P.ParkJ. W.Sample size and statistical power calculation in genetic association studies201210117122