- Open Access
Lack of association between four SNPs in the SLC22A3-LPAL2-LPA gene cluster and coronary artery disease in a Chinese Han population: a case control study
Lipids in Health and Diseasevolume 11, Article number: 128 (2012)
Lipoprotein (a) (Lp [a]) is known being correlated with coronary artery disease (CAD). The SLC22A3-LPAL2-LPA gene cluster, relating with modulating the level of plasma Lp (a), has recently been reported to be associated with CAD in Caucasians. The purpose of this study was to verify whether this finding can be expanded to the Chinese Han population.
Methods and Results
Using a Chinese Han sample, which consisted of 1012 well-characterized CAD patients and 889 healthy controls, we tested the associations of four SNPs (rs2048327, rs3127599, rs7767084 and rs10755578) in the SLC22A3-LPAL2-LPA gene cluster, and their inferred haplotypes with the risk of CAD. Allelic, genotypic and haplotype association analyses all showed that the gene cluster was not associated with CAD in this Chinese Han sample.
We for the first time explored the association of the four SNPs in the SLC22A3-LPAL2-LPA gene cluster with CAD in a large Chinese Han sample. Nevertheless, this study did not reveal any significant evidence of this gene cluster to increase the risk of CAD in this population.
Coronary artery disease (CAD) is now widespread and becoming a heavy burden for both developed and developing countries. A long list of susceptibility loci for CAD has been identified in previous genetic studies. Most recently, the gene cluster aligned by solute carrier family 22 member 3 (SLC22A3), lipoprotein(a)-like 2 (LPAL2), and lipoprotein(a) (LPA), known as SLC22A3-LPAL2-LPA on chromosome 6q26-27, attracted much attention, possibly due to its capability to regulate the plasma level of lipoprotein(a) (Lp(a)) .
Lp(a), a low-density lipoprotein (LDL)-like particle synthesized in the liver , has been well known as an independent risk factor for CAD [3–6]. Its protein component-apolipoprotein(a)[apo(a)], accounts for 91% variation of the plasma Lp(a) concentration . LPA, one member of the SLC22A3-LPAL2-LPA gene cluster, encodes apo(a) , and is associated with Lp(a) levels, explaining up to 36% of Lp(a) variance in European-descent [5, 8, 9]. In addition, the hapoltype formed by four SNPs (rs2048327, rs3127599, rs7767084 and rs10755578) in this region is also related to the plasma level of Lp(a) .
Researchers tried to explore the association between the SLC22A3-LPAL2-LPA gene cluster and CAD. Tregouet et al. identified this region as a risk cluster for CAD in the genome-wide haplotype study (GWHS) in six white populations . Koch et al. demonstrated that the gene cluster was a strong susceptive locus for MI in the European . However, the research of Qi et al. did not confirm the association between haplotypes in the SLC22A3-LPAL2-LPA region and nonfatal acute MI risk in Hispanics .
These contrary data indicate that more large-scale and independent studies should be performed to confirm the association between this cluster and CAD and verify whether this finding can be expanded to other populations. Up to date, there is no report on the correlation between SLC22A3-LPAL2-LPA and CAD in Chinese population, which takes up one fifth of the human population.
We performed a statistical power analysis using the PS program to verify whether the recruited samples could provide adequate power in identifying the association between the SNPs and CAD. Under the assumption of odd ratio being 1.236, and the risk allele frequency being 0.341, as previously reported , our sample size with 1012 well-characterized CAD cases and 889 healthy controls can provide a statistical power of 88.0% and 71.2% at the nominal type I error rate of 0.05 and 0.01 respectively. The power analysis indicated that our sample size is sufficient for identifying the modest-effect-size SNP.
Characteristics of participants
This case–control study included 1012 CAD patients and 889 healthy controls. Characteristics of the participants are summarized in Table 1. Compared with controls, the CAD patients were older, more likely to have a higher BMI, and to be current smokers.
Characteristics of four SNPs
All the call rates for the 4 SNPs were above 99%. All the four SNPs tested were polymorphic, with minor allele frequency ranging from 0.124 to 0.412, and in agreement with Hardy-Weinberg equilibrium. The linkage disequilibrium between rs2048327 and rs10755578 was relatively low (D’= 0.43 in cases, and D’= 0.48 in controls). D’ values between other SNPs ranged from 0.94 to 1. The information about the four SNPs was shown in Table 2 and Figure 1.
Allelic association and genotypic association
First, the associations between polymorphisms of rs2048327, rs3127599, rs7767084, rs10755578 and CAD were evaluated. OR and 95% CI for each SNP were calculated with the ancestral allele as the reference, and the statistical significance was defined by the permutation test mentioned in the method section. There was no significant association between any allele of the four SNPs and CAD (P-obs = 0.23-0.50, OR: 0.92-1.10). Permutation P values calculated using 100 000 Monte Carlo simulations were of no significance, either (as shown in Table 3).
We further examined the association between individual SNP and CAD risk under additive, dominant or recessive model, respectively. In the univariate logistic regression, none of the four SNPs showed significant association with CAD risk (P = 0.09-0.84), and further adjustment for the conventional risk factors such as age, sex, smoking and BMI did not change the results (P = 0.12-0.77). The results were list in Table 4.
Haplotype association analysis
The frequencies of haplotypes from different blocks were estimated and all the haplotypes with a frequency above 1% were included in the following analysis. In the haplotype association analysis, one haplotype was treated as a single variant, and all the other haplotypes were collapsed into the alternative allele to test its association with CAD. Age, sex, smoking and BMI were adjusted. The results were shown in Table 5.
Only two of all the tested haplotypes, TCC (which was made up by rs2048327, rs3127599 and rs10755578), and TC (which was made up by rs2048327 and rs10755578), reached the single-point significance level (P-obs = 0.04 and 0.03). However, the two findings could not survive over multiple test corrections. The P-values for the two haplotypes after the permutation test using 100 000 Monte Carlo simulations were 0.33 and 0.29, respectively. Moreover, the overall tests for all haplotypes formed by the three and the two SNPs did not achieve significance (P-obs = 0.17, with df = 4 and P-obs = 0.16, with df = 3), either. Haplotypes TC and TG formed by rs2048327 and rs10755578 were reported to be the risk haplotypes in a Japanese sample . However, haplotype TG showed no significant association with CAD (P-obs = 0.41, P-emp = 0.98) in this Chinese Han sample, as well as haplotype TC mentioned above (as shown in Table 5).
In previous studies of European populations, the most common haplotype formed by the four SNPs was TCTC, and haplotypes CCTC, CTTG and TTTC were found to be associated with increased risk of CAD [1, 10]. In this study, the most common haplotype was TCTG and the haplotype CTTG was not found. Haplotypes CCTC and TTTC showed no significant association with CAD (P-obs = 0.11, P-emp = 0.64 and P-obs = 0.30, P-emp = 0.95, separately). With the most common haplotype TCTG as the reference, OR and 95% CI for haplotypes CCTC and TTTC were 0.77 (0.45-1.32) and 1.15 (0.86-1.53). Furthermore, the corresponding OR and 95% CI for haplotypes CCTC and TTTC were 0.47 (0.17-1.30) and 0.694 (0.28-1.74), when the haplotype TCTC was used as the reference, which could not reach statistical significance, either. Data were not shown here.
In 2009, Tregouet et al. identified the SLC22A3-LPAL2-LPA gene cluster as a risk cluster and haplotypes CTTG and CCTC formed by rs2048327, rs3127599, rs7767084 and rs10755578 as risk haplotypes for CAD in six White populations . From then on, several GWHS have focused on this hot spot. In a study consisted of 3657 patients with MI and 1211 control individuals, Koch et al. observed significant association between haplotypes formed by the same four SNPs in the SLC22A3-LPAL2-LPA region and MI (P = 0.0005), and found 3 risk haplotypes (CTTG, CCTC, and TTTC) . Later, Sawabe M etal analyzed rs2048327 (C/T) and rs10755578 (C/G) in 1,150 Japanese autopsy cases, and ascertained that haplotypes TC and TG worked as risk factors for both coronary sclerosis and CAD . In addition, Shaw et al. found that genetic variants at the SLC22A3-LPAL2-LPA locus were associated with decreased early-outgrowth colony-forming units, thereby increased the risk of MI , which may support the findings in population studies mentioned above. However, Qi et al. did not confirm the association of haplotypes at the SLC22A3-LPAL2-LPA locus with nonfatal MI risk in Hispanics .
For gene association studies, repeating previous findings across different populations is essential for exploring the full scape of their pathogenic nature. To date, there is no study focusing on the association between CAD and the SLC22A3-LPAL2-LPA gene cluster in Chinese people. Our study for the first time attempted to explore such association in Chinese Hans. We evaluated the association between four SNPs in this gene cluster and CAD by examming all kinds of associations (allelic, genotypic and haplotype). Nevertheless, we did not identify any significant evidence to link this gene cluster to CAD risk in this Chinese Han sample. The genotypic and allelic association between individual SNP and CAD drawn from our data were consistent with results from previous GWHS . Whereas, there are differences exist between our study and previous studies. The most common haplotype we found was TCTG instead of TCTC, which was reported in European populations [1, 10, 11]. Moreover, we did not confirm the association of haplotypes CTTG, CCTC, and TTTC with CAD reported in European populations [1, 10]. In addition, we did not find any association between CAD and two haplotypes TC and TG composed by rs2048327 and rs10755578, which was inconsistent with the results from a Japanese study .
There are many reasons for heterogeneity in genetic association studies. Ethnic differences in genetic structure may produce different LD, thereby differences in the significance of the association test, which also exist in other genetic association studies . Besides, differences in environmental, dietary or behavioral factors may also partially explain the heterogeneity in the genetic associations across ethnicities [15, 16]. Furthermore, different disease definitions under different criteria may also be partly responsible for the variation between studies.
To limit the potential influence of factors mentioned above, we carefully designed and implemented this study. First, we used an adequate sample with enough statistical power, to detect the genetic association, therefore, the discrepancies between our study and others in different populations may be more likely due to the ethnic differences in genetic structure. Second, we performed multivariate logistic model to adjust several possible covariates, such as age, gender, smoking status and BMI. Moreover, we identified the case subjects in a strict accordance with a generally accepted definition of CAD and excluded patients taking niacin which could decrease the plasma level of Lp(a) and/or patients with diabetes since diabetes status was reported to attenuate the relation between Lp(a) and cardiovascular risk .
Despite our study was well organized, several limitations still exist in this exploratory study. First, as a complex disease, many factors may contribute to CAD, such as environmental and polygenic backgrounds, dietary and behavioral factors, hence, the genetic parameter estimates (odd ratios, risk allelic or genotype frequencies) may be biased. In addition, uncontrolled confounding factors may lead to spurious associations. Although many important confounding factors were controlled or adjusted in our analysis, some potential confounders, such as lipid level, were unavailable for a large number of subjects and thus not controlled.
Our study for the first time explored the association between CAD and the four SNPs in the SLC22A3-LPAL2-LPA gene cluster in Chinese Hans. We found no allelic, genotypic and haplotype association between the four SNPs in the SLC22A3-LPAL2-LPA gene cluster and CAD.
The method we recruited participants was described prevously . Briefly, all subjects were of the ethnic Han origin and are not related to each other. Cases and controls are gender-frequency matched. Patients with CAD were recruited from the Guangzhou Military Region General Hospital and the healthy controls were randomly selected from several communities in Guangzhou. CAD was diagnosed if any of the following criteria was met: (1) confirmed MI; (2) ≥50% stenosis in at least one coronary vessel at angiography; (3) a validated history of percutaneous transluminal coronary angioplasty or coronary artery bypass graft surgery and pectoris; (4) primarily diagnosed by symptoms and later confirmed by at least one non-invasive provocation test, e.g. scintigraphy or treadmill test . Patients with diabetes and/or those taking niacin were excluded. Full hospital records were reviewed to confirm the diagnosis. At the enrolment, anthropometric measures and drug uses for both cases and controls were collected by well-trained interns and physicians. Body mass index (BMI) was calculated using the formula: BMI=weight in kilograms/ (height in meters×height in meters). The study was approved by the Ethics Committees of the participating hospitals and institutions. All participants have signed the written informed consent forms. The investigation conformed to the principles outlined in the Declaration of Helsinki.
Genomic DNA was isolated from peripheral blood lymphocytes using the TIANamp Blood DNA Kit (Tiangen Inc; Beijing; China). The four SNPs rs204832, rs3127599, rs7767084 and rs10755578 were genotyped using Illumina Golden Gate Genotyping Bead Chips (Illumina Inc; San Diego; USA) (http://www.Illumina.Com/), which uses illumiCodes, unique 23-bp single stranded DNA oligos, to correctly identify each DNA as well as the loci being interrogated . First, prepared DNA samples were amplified using universal PCR primers labeled with Cy3 and Cy5 fluorescent dyes, and the resulting fluorescently labeled PCR product was then hybridized to a Universal BeadChip, which contained randomly assembled universal beads, each displaying an illumiCode corresponding to specific loci.
The quality for SNP genotyping was assured by independently replicating the genotyping and allelic calls of 30 randomly selected samples. The results from quality control were in perfect agreement with the initial genotyping results. In addition, all the DNA samples for cases and controls were run in the same batches.
Continuous covariates were expressed as meaxn ± SD, and the differences between cases and controls were analyzed by independent t-test or Mann–Whitney U test. Categorical variables were summarized as frequency (percentage) and analyzed by Chi-square test. Genotype coding method of Lewis was used to build additive, dominant and recessive genetic models . Odds ratio (OR) and Wald 95% confidence interval (CI) were calculated using the homozygote of ancestral allele as reference. The ancestral alleles were defined by the dbSNP database (http://www.ncbi.nlm.nih.gov/snp). Genotypic association between each SNP and CAD was analyzed by using univariate (with only SNP included) and multivariate BMI logistic regression under additive, dominant and recessive genetic models. In the multivariate analysis, age, gender, smoking and BMI were also modeled and adjusted. All statistical analyses mentioned above were performed with the SPSS 13.0 (SPSS Inc; Chicago; USA).
Prior to the experiment, the case–control design and the sample size were evaluated by the PS program, a software for power analysis . The Hardy-Weinberg equilibrium test and linkage disequilibrium (LD) between the SNPs were analyzed by employing Haploview 4.2 (http://www.broad.mit.edu/mpg/haploview). The extent of disequilibrium was expressed by D’. Haplotype frequency estimation, allelic and haplotype association analyses were implemented and then followed by permutation analysis with 100 000 Monte Carlo simulations by using PLINK software (http://pngu.mgh.harvard.edu/~purcell/plink). In the allelic association analysis, OR and 95%CI were calculated using the ancestral allele as reference. The associations between all the haplotypes drawn from every two, every three or all four SNPs and CAD were analyzed, except for those with frequency below 1%.
Coronary artery disease
Solute carrier family 22 member 3
Genome-wide haplotype study
Body mass index
Trégouët DA, König IR, Erdmann J, Munteanu A, Braund PS, Hall AS, Grosshennig A, Linsel-Nitschke P, Perret C, DeSuremain M, Meitinger T, Wright BJ, Preuss M, Balmforth AJ, Ball SG, Meisinger C, Germain C, Evans A, Arveiler D, Luc G, Ruidavets JB, Morrison C, van der Harst P, Schreiber S, Neureuther K, Schäfer A, Bugert P, El Mokhtari NE, Schrezenmeir J, Stark K: Genome-wide haplotype association study identifies the SLC223A-LPL2-LPA gene cluster as a risk locus for coronary artery disease. Nature Genet. 2009, 41: 283-285. 10.1038/ng.314
Kraft HG, Menzel HJ, Hoppichler F, Vogel W, Utermann G: Changes of genetic apolipoprotein phenotypes caused by liver transplantation. Implications for apolipoprotein synthesis. J Clin Invest. 1989, 83: 137-142. 10.1172/JCI113849
Erqou S, Kaptoge S, Perry PL, Di Angelantonio E, Thompson A, White IR, Marcovina SM, Collins R, Thompson SG, Danesh J, : Lipoprotein(a) concentration and the risk of coronary heart disease, stroke, and nonvascular mortality. JAMA. 2009, 302: 412-423.
Nordestgaard BG, Chapman MJ, Ray K, Borén J, Andreotti F, Watts GF, Ginsberg H, Amarenco P, Catapano A, Descamps OS, Fisher E, Kovanen PT, Kuivenhoven JA, Lesnik P, Masana L, Reiner Z, Taskinen MR, Tokgözoglu L, Tybjærg-Hansen A: Lipoprotein(a) as a cardiovascular risk factor: Current status. Eur Heart J. 2010, 31: 2844-2853. 10.1093/eurheartj/ehq386
Clarke R, Peden JF, Hopewell JC, Kyriakou T, Goel A, Heath SC, Parish S, Barlera S, Franzosi MG, Rust S, Bennett D, Silveira A, Malarstig A, Green FR, Lathrop M, Gigante B, Leander K, de Faire U, Seedorf U, Hamsten A, Collins R, Watkins H, Farrall M: Genetic variants associated with lp(a) lipoprotein level and coronary disease. N Engl J Med. 2009, 361: 2518-2528. 10.1056/NEJMoa0902604
Dahlen GH, Stenlund H: Lp(a) lipoprotein is a major risk factor for cardiovascular disease: Pathogenic mechanisms and clinical significance. Clin Genet. 1997, 52: 272-280.
Boerwinkle E, Leffert CC, Lin J, Lackner C, Chiesa G, Hobbs HH: Apolipoprotein(a) gene accounts for greater than 90% of the variation in plasma lipoprotein(a) concentrations. J Clin Invest. 1992, 90: 52-60. 10.1172/JCI115855
Ober C, Nord AS, Thompson EE, Pan L, Tan Z, Cusanovich D, Sun Y, Nicolae R, Edelstein C, Schneider DH, Billstrand C, Pfaffinger D, Phillips N, Anderson RL, Philips B, Rajagopalan R, Hatsukami TS, Rieder MJ, Heagerty PJ, Nickerson DA, Abney M, Marcovina S, Jarvik GP, Scanu AM, Nicolae DL: Genome-wide association study of plasma lipoprotein(a) levels identifies multiple genes on chromosome 6q. J Lipid Res. 2009, 50: 798-806.
Lanktree MB, Anand SS, Yusuf S, Hegele RA: Comprehensive analysis of genomic variation in the lpa locus and its relationship to plasma lipoprotein(a) in south asians, chinese, and european caucasians. Circ Cardiovasc Genet. 2010, 3: 39-46. 10.1161/CIRCGENETICS.109.907642
Koch W, Hoppmann P, Erl A, Schrempf M, Schatke A, Tuerk S, Mueller JC, Schoemig A, Kastrati A: Risk polymorphisms and haplotypes identified in genome-wide association studies and their relation with myocardial infarction. Eur Heart J. 2010, 31: 801-
Qi L, Ma J, Qi Q, Hartiala J, Allayee H, Campos H: Genetic risk score and risk of myocardial infarction in hispanics. Circulation. 2011, 123: 374-380. 10.1161/CIRCULATIONAHA.110.976613
Sawabe M, Tanaka N, Arai T, Mieno MN, Ikeda S, Muramatsu M: Genetic polymorphisms of SLC22A3-LPL2-LPA gene cluster has direct and indirect effects on coronary atherosclerosis and coronary heart disease. HUGO J. 2009, 4: 107-
Shaw SY, Cheng S, Cupples LA, Larson MG, McCabe EL, Ngwa JS, Wang YA, Martin RP, Klein RJ, Hashmi B, Ajijola OA, Lau E, O’Donnell CJ, Vasan RS, Cohen KS, Wang TJ: Genetic and clinical correlates of early-outgrowth colony-forming units. Circ Cardiovasc Genet. 2011, 4: 296-304. 10.1161/CIRCGENETICS.110.958470
Lanktree MB, Anand SS, Yusuf S, Hegele RA: Replication of genetic associations withplasma lipoprotein traits in a multiethnic sample. J Lipid Res. 2009, 50: 1487-1496. 10.1194/jlr.P900008-JLR200
Qi L, Hu FB, Hu G: Genes, environment, and interactions in prevention of type 2 diabetes: A focus on physical activity and lifestyle changes. Curr Mol Med. 2008, 8: 519-532. 10.2174/156652408785747915
Qi L, Cho YA: Gene-environment interaction and obesity. Nutr Rev. 2008, 66: 684-694. 10.1111/j.1753-4887.2008.00128.x
Qi Q, Workalemahu T, Zhang C, Hu FB, Qi L: Genetic variants, plasma lipoprotein(a) levels, and risk of cardiovascular morbidity and mortality among two prospective cohorts of type 2 diabetes. Eur Heart J. 2011, Epub ahead of print..
Lv X, Zhang Y, Rao S, Qiu J, Wang M, Luo X, Zuo X, Su D, Feng X, Yang Y, Ouyang P, Chen Y, Li X, Xiao Y, Ling W: Joint effects of genetic variants in multiple Loci on the risk of coronary artery disease in chinese han subjects. Circ J. 2012, 76 (8): 1987-1992. 10.1253/circj.CJ-12-0156
Samani NJ, Deloukas P, Erdmann J, Hengstenberg C, Kuulasmaa K, McGinnis R, Schunkert H, Soranzo N, Thompson J, Tiret L, Ziegler A: Large scale association analysis of novel genetic loci for coronary artery disease. Arterioscler Thromb Vasc Biol. 2009, 29: 774-780.
Pask R, Rance HE, Barratt BJ, Nutland S, Smyth DJ, Sebastian M, Twells RC, Smith A, Lam AC, Smink LJ, Walker NM, Todd JA: Investigating the utility of combining phi29 whole genome amplification and highly multiplexed single nucleotide polymorphism beadarray genotyping. BMC Biotechnol. 2004, 4: 15- 10.1186/1472-6750-4-15
Lewis CM: Genetic association studies: Design, analysis and interpretation. Brief Bioinform. 2002, 3: 146-153. 10.1093/bib/3.2.146
Dupont WD, Plummer WD: Power and sample size calculations for studies involving linear regression. Control Clinical Trials. 1998, 19: 589-601. 10.1016/S0197-2456(98)00037-3.
We give our sincere thanks to all individuals who participated in the present study and to the medical personnel in Department of Cardiovascular of the Guangzhou Military Region General Hospital for their assistance. We are grateful to Yibing Chen and Jianrong He for revising the paper. This research was supported by grants from the National Natural Science Foundation of China (No. 81130052, 81010017 and 31071166 ), Doctor's Launching Project, Guangdong Province Natural Science Fund (Grant No.10451008901006170), and Science and Technology Planning Project of Guangdong Province (Grant No. 2009A030301004).
The author’s declare that they have no competing interests'.
JQ, SR, WL and MX conceived of the project and participated in its design, revised the paper, and they gave final approval of the version to be published. YZ, XL, MW, DS and XL collected clinical data and extracted DNA. XL and JQ carried out statistical analysis and wrote the manuscript. FL and XF took charge of genotyping, and helped to draft the manuscript. XZ, DF, HG, DL,Y Y and PO helped to statistical analysis and revised the paper. XL, YZ and SR contributed equally to this work. All authors read and approved the final manuscript.
Xiaofei Lv, Yuan Zhang, Shaoqi Rao contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.