Association study of genetic variants at newly identified lipid gene TRIB1 with coronary heart disease in Chinese Han population
- Long Wang†1,
- Jinjin Jing†2,
- Qianxi Fu1,
- Xiaojun Tang1,
- Li Su2,
- Shishi Wu3,
- Ge Li1, 4Email author and
- Li Zhou1, 4Email author
© Wang et al. 2015
Received: 18 December 2014
Accepted: 8 May 2015
Published: 19 May 2015
Recent genome-wide association studies (GWAS) have identified the variants near TRIB1 gene affecting blood lipid levels. However, the association between the reported variants and risk of coronary heart disease (CHD) was not confirmed.
We conducted two independent case–control studies. The first study consisted of 300 CHD patients and 300 controls and the second study had 1,332 CHD patients and 2,811 controls. The genotypes of two variants rs3201475 and rs17321515 in TRIB1 were determined by TaqMan assay. The dual-luciferase reporter assay was performed for evaluating the function of the SNP rs3201475.
The statistical analysis indicated that single nucleotide polymorphism (SNP) rs17321515 was replicated to be associated with triglyceride (TG) level, which was also significantly associated with CHD risk when using the stratified analysis after adjusting for conventional risk factors. Compared with GG genotype, AA carriers of SNP rs17321515 had higher risk in males (odds ratio (OR) = 1.28, 95 %CI = 1.01–1.61; P = 0.03) and smokers (OR = 1.41, 95%CI = 1.09–1.88; P = 0.01). We did not find significantly association between genotypes of rs3201475 and CHD risk. In addition, no significant difference was found in the luciferase activity assay of SNP rs3201475.
Our findings indicated that SNP rs17321515 is significantly associated with plasma TG level and the increasing risk of CHD among males and smokers in Chinese, whereas there is no positive association between SNP rs3201475 and CHD risk. Smoking could modify the effects of TRIB1 on CHD risk.
KeywordsGenetic variant Lipid Coronary heart disease
Coronary heart disease (CHD) is the leading cause of death worldwide . World Health Organization (WHO) estimated that more than 700,000 people died from CHD each year. Recent genome-wide association studies (GWAS) have identified SNPs at several loci such as chromosomes 1p13, 3q22, 6p24, 9p21 and 15q22 that are associated with risk of CHD –. Meanwhile, a few loci were also shown to be strongly associated with plasma lipid levels, reinforcing the close mechanistic association between the variability in lipid levels and CHD risk.
TRIB1 gene was firstly identified as a lipid-associated gene by Kathiresan through GWAS in the European population . Then some other GWAS confirmed the association between TRIB1 and lipids levels. Kathiresan et al. identified that the variants in TRIB1 gene were associated with plasma triglyceride (TG) level . Aulchenko et al. found that single nucleotide polymorphism (SNP) rs17321515 in TRIB1 gene was associated with plasma total cholesterol (TC) level . Recently, we conducted a GWAS of lipid levels in Chinese Han population and the data from SNP rs17321515 in TRIB1 confirmed the association previously reported in Europeans .
Blood lipid levels have been consistently associated with risk of CHD. However, the relationship of TRIB1 and CHD risk was indefinite. Meanwhile, the real functional SNP in TRIB1 was not identified in previous studies. To further establish the associations between the genetic variant in TRIB1, lipid levels and risk of CHD and to improve understanding of the mechanisms underlying susceptibility to CHD, we conducted two independent CHD case–control studies in Chinese Han populations to test two SNPs rs17321515 and rs3201475 in TRIB1 gene. SNP rs3201475, in the 5′ untranslated region (UTR) of TRIB1, was predicted to be a transcription factor binding site (TFBS) according to the transcription factor prediction tools: SNP Function Prediction and the JASPAR database. The prediction results showed that the SNP rs3201475 might have regulatory function in TRIB1. So we also examined whether the SNP rs3201475 affected gene expression by performing a reporter gene luciferase activity assay in two types of cell lines.
Materials and methods
We performed two independent case–control analyses. The first population group consisted of 300 patients with CHD and 300 age- and sex-frequency matched healthy control subjects. The second population group consisted of 1,332 CHD cases and 2,811 controls. Patients of the former were consecutively recruited from three hospitals in Chongqing city between March 2012 and May 2013, and the latter were from hospitals in Wuhan and Chongqing cities between May 2008 and October 2013. The inclusion criteria for CHD were either: 1. the presence of a stenosis ≥ 50 % in at least one of the major segments of coronary arteries (the right coronary artery, left circumflex, or left anterior descending arteries) on coronary angiography; 2. based on World Health Organization criteria in terms of elevations of cardiac enzymes, electrocardiographic changes and clinical symptoms ; 3. a documented history of coronary artery bypass graft or percutaneous coronary intervention. Patients with congenital heart disease, cardiomyopathy and valvular disease were excluded. After cases were diagnosed with CHD, they were interviewed in person by a trained interviewer within 3 days. The control subjects, residing in the same communities as the cases, were determined to be free of CHD and peripheral atherosclerotic arterial disease by medical history, clinical examinations, and electrocardiography. Subjects with severe liver and/or kidney disease were excluded. Medical history, family history of CHD among first degree relatives, medication use, home environment, and lifestyle factors were obtained through questionnaire interview.
Subjects were classified as smokers and nonsmokers. Those who had smoked less than 100 cigarettes in the lifetime were defined as nonsmokers; otherwise, they were defined as smokers. The smokers were categorized into groups of cigarettes per day: < 20 and ≥20. BMI was calculated as weight in kilograms divided by the square of height in meters. All subjects of the two populations gave written consent after receiving a full explanation of the study. The Ethics Committee of Chongqing Medical University approved the present study, and written informed consent was obtained from all subjects.
Determination of lipid levels
The plasma TC, TG, LDL and HDL levels were measured by the ARCHITECT Ci8200 automatic analyzer (ABBOTT Laboratories. Abbott Park, Illinois, U.S.A.) using the Abbott Diagnostics reagents according to the manufacturer’s instructions.
Fasting venous blood was collected in 5-ml EDTA tubes, and genomic DNA was isolated with a Puregene kit (Gentra Systems, Inc., Minneapolis, MN, USA). Genotyping was performed with TaqMan SNP allelic discrimination by means of an ABI 7900HT (Applied Biosystems, Foster City, CA, USA), in 384-well format. The TaqMan Assay kit was purchased from Applied Biosystems (Foster City, CA, USA). It included the forward target-specific polymerase chain reaction (PCR) primer, the reverse primer, and the TaqMan MGB probes labeled with two special dyes: FAM and VIC. PCR reactions were carried out in reaction volume of 5 μl containing 5 ng DNA, 2.5 μl 2 × Taqman universal PCR Master MixNo AmpErase UNG (Applied Biosystems, Foster City, CA, USA), 0.125 μl 40 × Assay MIX. PCR conditions included 95 °C for 10 min followed by 40 cycles of 15 s at 92 °C and 1 min at 60 °C. Two blank controls (DNA hydration) and two replicate quality control samples were included in each 384-well format, and two replicate samples were genotyped with 100 % concordance. Automatic allele calling, with the default settings (the quality value of auto caller ≥95.0), was carried out by ABI 7900HT data collection and analysis software version 2.2.1 (SDS 2.2.1).
Reporter plasmid construct
A 1,929 bp-sized fragment (from −1316 to +613) which includes the promoter and 5′ UTR of the TRIB1 gene was polymerase chain reaction (PCR) amplified using either +164 T (rs3201475) homozygous or +164 C homozygous genomic DNA as a template and the following primers: forward primer: 5′– CGGGGTACC GGTGCCCAGGGACTCCAAAC −3′; the italic bold characters represent the KpnI site, reverse primer: 5′– CCCAAGCTT TCATGGCAGAGCGCACCGGA −3′; the italic bold characters denote the HindIII site. PCR conditions were as follows: Initial denaturation for 3 min at 95 °C, 30 cycles of denaturing at 98 °C for 40 s, annealing at 66 °C for 15 s and extension at 72 °C for 5 min. PCR products were purified using High pure PCR product purification Kit1 (Roche, Switzerland). The purified PCR product was sub-cloned into the KpnI–HindIII site of the pGL3-Basic luciferase reporter vector (Promega, Madison, WI, USA) to generate the pGL3−+164 plasmid. Site-directed mutagenesis at the +164 C/T site of the pGL3−+164 plasmid was carried out using a Muta-direct sitedirected mutagenesis kit (Saibaisheng, Beijing, China) following the manufacturer’s protocol. This process generated new construct with the +164 C/T site changed to +164 T/C site. These mutations were confirmed by DNA sequencing.
Cell culture and luciferase assays
HepG2 and HeLa cells were maintained in Dulbecco’s modified Eagle’s medium (Gibco BRL, USA) supplemented with 10 % fetal bovine serum, penicillin (100 U/ml), and streptomycin (100 U/ml). Each constructed pGL3−+164 plasmid was transfected into HepG2 and HeLa cells with the pRL-SV40 plasmid (Promega, USA) by using Lipofectamine 2000 (Invitrogen, USA). Firefly luciferase and renilla luciferase activities were sequentially measured by a luminometer 48 h after transfection, utilizing a Dual-Luciferase reporter assay system (Promega, USA). Results were expressed as relative light units of firefly luciferase activity over relative light units of renilla luciferase activity. All experiments were performed in triplicate and repeated three times.
Each continuous trait (TC, TG, LDL and HDL) was tested for normality and TG values were log-transformed. Continuous variables were reported as mean ± standard deviation (SD). Normal distribution of data was analyzed using the Kolmogorov-Smirnov normality test. Data with a normal distribution were compared by Student’s t-test, and those with unequal variance or without a normal distribution were analyzed by a Mann–Whitney rank sum test. Categorical values were compared by the chi-square test, which was also used to test for deviation of genotype distribution from Hardy-Weinberg equilibrium. The association between SNP and CHD risk was estimated by computing odds ratios (ORs) and 95 % confidence intervals (CIs) from the multivariate logistic regression analyses. The effects of genotypes on plasma lipid levels were assessed by multiple linear regression models. The Haploview program was used to analyze pairwise linkage disequilibrium (LD) with data extracted from International HapMap Project . The probability level accepted for significance was P < 0.05. All data analyses were carried out with the statistical analysis software package SPSS 16.0 (SPSS Inc., Chicago, IL, USA).
General characteristics of the subjects
General characteristics of CHD patients and controls
Cases (N = 300)
Controls (N = 300)
Cases (N = 1332)
Controls (N = 2811)
Sex, m/f, (%)
Age, mean ± SD
68.8 ± 10.3
67.6 ± 15.0
60.8 ± 10.0
61.2 ± 9.1
Body mass index(kg/m2)
23.8 ± 3.1
22.7 ± 4.6
24.5 ± 3.6
24.3 ± 3.5
Smoking, no/yes, (%)
6.9 ± 2.8
7.2 ± 4.2
6.6 ± 3.5
5.3 ± 2.0
4.3 ± 1.3
4.4 ± 1.2
4.7 ± 1.1
4.4 ± 1.0
1.7 ± 1.5
1.5 ± 1.0
1.7 ± 1.2
1.6 ± 1.3
0.9 ± 0.3
1.1 ± 0.3
1.0 ± 0.4
1.2 ± 0.3
2.5 ± 1.0
2.6 ± 0.9
2.7 ± 0.9
2.6 ± 0.8
Relation of the two SNPs in TRIB1 and CHD risk
Association of two SNPs in TRIB1 with CHD cases and controls in the Chinese Han populations
Crude OR (95%CI)
Adjusted ORa (95%CI)
rs3201475 C > T
rs17321515 G > A
Stratification analysis for association between SNP rs17321515 genotypes and risk of CHD
OR (95 % CI)a
Association between SNP rs17321515 in TRIB1 and TG level
Reporter gene luciferase activity assay
In this study, we conducted two independent case–control studies to investigate whether a new genetic susceptibility locus TRIB1 was associated with lipid levels and involved in the development of CHD. However, only one SNP rs17321515 near TRIB1 had modestly association with CHD susceptibility in Chinese Han population. Another SNP rs3201475 which was in the 5′ UTR region of TRIB1 had no significant association in case–control study and luciferase assays.
We identified the SNP rs17321515 in TRIB1 that influence plasma TG level in Chinese, and also provided evidence that this SNP was significantly associated with CHD risk in a large independent case–control study in Chinese. Our results showed that genotype AA of the SNP rs17321515 was associated with increased risk of CHD among males and smokers. It is possible that genetic factors would exert a greater influence in males, and smoking may exacerbate the influence of the genetic factors. However, there was no significant relationship between the SNP rs3201475 and CHD susceptibility in Chinese Han population. In addition, no significant difference between C and T allele transcriptional activity was observed when using HepG2 and HeLa cell lines.
Significant associations of the SNP rs17321515 at 8q24 near TRIB1 with lipid levels have been replicated in several studies. The minor G allele at this SNP was associated with lower TG, lower LDL and higher HDL levels. Kathiresan et al. identified SNP rs17321515 was strongly associated with TG and was also associated with LDL and HDL levels . Wang et al. found that TRIB1 had significant high risk for severe hypertriglyceridemia . Garcia-Rios et al. reported that rs17321515 was associated with familial hypercholesterolemia in Spanish population . Recently, our GWAS also showed that SNP rs17321515 in TRIB1 was significantly associated with TG level in Chinese Han population . However, few studies have examined the associations between the SNPs in TRIB1 and CHD risk. Willer et al. reported that SNP rs17321515 in TRIB1 which associated with increased TG level was also associated with increased risk of cardiovascular disease . Recently, a Danish study replicated the association between another SNP rs2954029 in TRIB1 and risk of ischemic heart disease in White descent in more than 71,000 individuals . Also, in an Asian Malay study, Tai et al. found rs17321515 to be associated with an increased risk of cardiovascular disease , which is consistent with our results.
The frequency of the risk allele A of rs17321515 is different between the Chinese Han population and European population. The data from the HapMap showed that the rs17321515-A allele differs from 0.305 in Europeans to 0.178 in Han Chinese. The SNPs which were strongly associated lipids levels and CHD risk lie considerably downstream of TRIB1 gene, and the most significantly SNP rs17321515 is about 44 kb downstream of it.
Rs3201475, a polymorphism in 5′ UTR of the TRIB1 gene, predicted to be a TFBS with the minor allele frequency (MAF) to be 0.122 in Chinese Han population according to the NCBI database. However, there was no literature had studied on it. In our study, no significant association was observed between the SNP rs3201475 and CHD susceptibility in Chinese Han population and no significant difference between C and T allele transcriptional activity levels in both HepG2 and HeLa cell lines, suggesting the variant is not CHD causal variant. Considering the possible causes, the lipid-associated SNPs in TRIB1 identified by GWAS were located 30 kb downstream from TRIB1 and this gene has not yet been identified as an expression trait locus (eQTL). It suggested that the genomic effects of this locus may not be limited to alterations in expression of the Trib1 protein . Furthermore, the binding of transcription factors to a specific DNA sequence is dependent on various intracellular and extracellular stimuli within a cell. Therefore, although some transcription factors are predicted to bind to these sequences, this does not take into account in vivo conditions such as chromatin structure, and are hence not ideal representations of what occurs in the cell .
TRIB1 encodes a G-protein-coupled receptor-induced protein, a member of the family of proteins that act as secondary messengers in the mitogen-activated protein kinases (MAPK)-related signaling cascade ,  and may regulate lipid metabolism through this pathway . It has also been suggested that TRIB1 expression is regulated by inflammatory stimulation . TRIB1 controls chemotaxis and proliferation of smooth muscle cells in the arterial intima . However, the exact mechanism of TRIB1 in the development of atherosclerosis is still not known. It is interesting to note that smoking acts as modifiers of the association between TRIB1 and CHD risk. Garcia-Rios et al. observed that the variants in TRIB1 was only associated with familial hypercholesterolemia in smokers . It has been demonstrated that nicotine induces MAPK dependent vascular smooth muscle cell migration resulting in the development of atherosclerosis and that cigarette smoking enhances MAPK activation , . These findings suggest that the smoking-induced associations between TRIB1 and CHD risk may involve MAPK.
Another study in Han and Mulao populations in China found that the lipid levels were associated with rs17321515 in males. Serum lipids were also associated with age, gender, BMI, blood pressure, blood glucose and alcohol consumption . All of them were risk factors of the CHD and we adjusted these traditional factors in our study.
There is limitation of our study that must be acknowledged. The subjects used to test the association between the genotypes and CHD risk did not exclude the individuals who were receiving lipid-lowering medications. About 17.4 % CHD cases and 0.2 % controls in our study reported taking cholesterol-lowering medications, so we could not explore the association between the genotypes and lipid levels in all subjects. Fortunately, the very low frequency in control subjects (0.2 %) did not significantly affect the results of the association between the genotype and lipid levels. We identified the effect of genetics on lipid levels in the control subjects.
In summary, SNP rs17321515 near TRIB1 was associated with plasma TG level, and this SNP was significantly associated with increased risk of CHD in males and smokers in Chinese. Smoking could modify the effects of TRIB1 on CHD risk. SNP rs3201475 in TRIB1 showed no significant association with CHD risk. Further studies are needed to validate these findings and investigate potential mechanisms underlying the links between the variation of TRIB1 gene and CHD risk.
SNP Function Prediction, http://snpinfo.niehs.nih.gov/snpinfo/snpfunc.htm; JASPAR database, http://jaspar.genereg.net/cgi-bin/jaspar_db.pl; International HapMap Project, http://hapmap.ncbi.nlm.nih.gov/.
The authors particularly thank all volunteers for assisting in collecting the data and samples and participating in the present study in Second Affiliated Hospital of Chongqing medical university, the Third People’s Hospital of Chongqing, the First People’s Hospital of Jiulongpo District in Chongqing city. This work was supported by the grants from the National Natural Science Foundation of China (grant NSFC-81202274).
- Lopez AD, Mathers CD, Ezzati M, Jamison DT, Murray CJ: Global and regional burden of disease and risk factors, 2001: systematic analysis of population health data. Lancet. 2006, 367: 1747-57.View ArticlePubMedGoogle Scholar
- Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008, 40: 189-97.PubMed CentralView ArticlePubMedGoogle Scholar
- Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, Pramstaller PP: Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet. 2009, 41: 47-55.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Ban MR, Zou GY, Cao H, Lin T, Kennedy BA: Polygenic determinants of severe hypertriglyceridemia. Hum Mol Genet. 2008, 17: 2894-9.View ArticlePubMedGoogle Scholar
- Zhou L, He M, Mo Z, Wu C, Yang H, Yu D: A genome wide association study identifies common variants associated with lipid levels in the Chinese population. PLoS One. 2013, 8: e82420-PubMed CentralView ArticlePubMedGoogle Scholar
- Nomenclature and criteria for diagnosis of ischemic heart disease. Report of the Joint International Society and Federation of Cardiology/World Health Organization task force on standardization of clinical nomenclature. Circulation. 1979; 59:607–609.View ArticleGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-5.View ArticlePubMedGoogle Scholar
- Garcia-Rios A, Perez-Martinez P, Mata P, Fuentes F, Lopez-Miranda J, Alonso R: Polymorphism at the TRIB1 gene modulates plasma lipid levels: insight from the Spanish familial hypercholesterolemia cohort study. Nutr Metab Cardiovasc Dis. 2011, 21: 957-63.View ArticlePubMedGoogle Scholar
- Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008, 40: 161-9.View ArticlePubMedGoogle Scholar
- Varbo A, Benn M, Tybjaerg-Hansen A, Grande P, Nordestgaard BG: TRIB1 and GCKR polymorphisms, lipid levels, and risk of ischemic heart disease in the general population. Arterioscler Thromb Vasc Biol. 2011, 31: 451-7.View ArticlePubMedGoogle Scholar
- Tai ES, Sim XL, Ong TH, Wong TY, Saw SM, Aung T: Polymorphisms at newly identified lipid-associated loci are associated with blood lipids and cardiovascular disease in an Asian Malay population. J Lipid Res. 2009, 50: 514-20.View ArticlePubMedGoogle Scholar
- Khamis A, Palmen J, Lench N, Taylor A, Badmus E, Leigh S: Functional analysis of four LDLR 5′UTR and promoter variants in patients with familial hypercholesterolaemia. Eur J Hum Genet. 2015, 23: 790-5.View ArticlePubMedGoogle Scholar
- Kiss-Toth E, Bagstaff SM, Sung HY, Jozsa V, Dempsey C, Caunt JC: Human tribbles, a protein family controlling mitogen-activated protein kinase cascades. J Biol Chem. 2004, 279: 42703-8.View ArticlePubMedGoogle Scholar
- Hegedus Z, Czibula A, Kiss-Toth E: Tribbles: a family of kinase-like proteins with potent signalling regulatory function. Cell Signal. 2007, 19: 238-50.View ArticlePubMedGoogle Scholar
- Burkhardt R, Toh SA, Lagor WR, Birkeland A, Levin M, Li X: Trib1 is a lipid- and myocardial infarction-associated gene that regulates hepatic lipogenesis and VLDL production in mice. J Clin Invest. 2010, 120: 4410-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Sung HY, Francis SE, Crossman DC, Kiss-Toth E: Regulation of expression and signalling modulator function of mammalian tribbles is cell-type specific. Immunol Lett. 2006, 104: 171-7.View ArticlePubMedGoogle Scholar
- Sung HY, Guan H, Czibula A, King AR, Eder K, Heath E: Human tribbles-1 controls proliferation and chemotaxis of smooth muscle cells via MAPK signaling pathways. J Biol Chem. 2007, 282: 18379-87.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Luozzo G, Pradhan S, Dhadwal AK, Chen A, Ueno H, Sumpio BE: Nicotine induces mitogen-activated protein kinase dependent vascular smooth muscle cell migration. Atherosclerosis. 2005, 178: 271-7.View ArticlePubMedGoogle Scholar
- Koch A, Giembycz M, Stirling RG, Lim S, Adcock I, Wassermann K: Effect of smoking on MAP kinase-induced modulation of IL-8 in human alveolar macrophages. Eur Respir J. 2004, 23: 805-12.View ArticlePubMedGoogle Scholar
- Aung LH, Yin RX, Wu DF, Li Q, Yan TT, Wang YM: Association of the TRIB1 tribbles homolog 1 gene rs17321515 A>G polymorphism and serum lipid levels in the Mulao and Han populations. Lipids Health Dis. 2011, 10: 230-PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.