- Research article
- Open Access
Genetic variants in LPL, OASL and TOMM40/APOE-C1-C2-C4 genes are associated with multiple cardiovascular-related traits
BMC Medical Genetics volume 12, Article number: 123 (2011)
Genome-wide association studies (GWAS) have become a major strategy for genetic dissection of human complex diseases. Analysing multiple phenotypes jointly may improve both our ability to detect genetic variants with multiple effects and our understanding of their common features. Allelic associations for multiple biochemical traits (serum alanine aminotransferase, aspartate aminotransferase, butrylycholinesterase (BCHE), C-reactive protein (CRP), ferritin, gamma glutamyltransferase (GGT), glucose, high-density lipoprotein cholesterol (HDL), insulin, low-density lipoprotein cholesterol (LDL), triglycerides and uric acid), and body-mass index, were examined.
We aimed to identify common genetic variants affecting more than one of these traits using genome-wide association analysis in 2548 adolescents and 9145 adults from 4986 Australian twin families. Multivariate and univariate associations were performed.
Multivariate analyses identified eight loci, and univariate association analyses confirmed two loci influencing more than one trait at p < 5 × 10-8. These are located on chromosome 8 (LPL gene affecting HDL and triglycerides) and chromosome 19 (TOMM40/APOE-C1-C2-C4 gene cluster affecting LDL and CRP). A locus on chromosome 12 (OASL gene) showed effects on GGT, LDL and CRP. The loci on chromosomes 12 and 19 unexpectedly affected LDL cholesterol and CRP in opposite directions.
We identified three possible loci that may affect multiple traits and validated 17 previously-reported loci. Our study demonstrated the usefulness of examining multiple phenotypes jointly and highlights an anomalous effect on CRP, which is increasingly recognised as a marker of cardiovascular risk as well as of inflammation.
Genome-wide association studies (GWAS) have become a major strategy for genetic dissection of human complex diseases. There is substantial overlap, both phenotypically and in allelic associations, between biomarkers and/or risk factors and between related diseases, and it is becoming important to understand the ways in which polymorphisms affect multiple phenotypes. Many phenotypes may be available from a single study population but current GWAS approaches usually examine them separately within a univariate framework. This strategy ignores potential genetic correlation between different traits.
From the perspective of maximising power for a given size of dataset, it has been shown that joint analyses of correlated traits in linkage analysis have substantially improved power in localizing genes [1–4]. Similarly, multivariate approaches in association studies can theoretically improve the ability to detect genetic variants whose effects are too small to be detected in univariate tests . Multivariate association tests have been proposed for unrelated samples  and for family data . Most of these tend to be inefficient and/or computationally intensive, especially at the genome-wide level. The approach proposed by Ferreira and Purcell has been shown to be powerful when traits have moderate to high correlation and efficient when applied to samples of unrelated individuals .
Genetically complex (multifactorial) diseases such as cardiovascular disease and type 2 diabetes often have common risk factors. A number of biochemical markers are known to be associated with obesity, pre-diabetic states, or risk of cardiovascular disease. Lipid traits such as triglycerides, and the low-density lipoprotein (LDL) and high-density lipoprotein (HDL) components of cholesterol, are well-known risk factors for cardiovascular disease. Other biochemical markers such as C-reactive protein (CRP) , the enzymes used as liver function tests (gammaglutamyl transferase, GGT [9–11], alanine aminotransferase, ALT; and aspartate aminotransferase, AST), butyrylycholinesterase (BCHE) [12, 13]), serum ferritin  and uric acid [15, 16] have also been shown to be associated with the risk of cardiovascular disease, hypertension, obesity, insulin resistance or metabolic syndrome. These biochemical markers are correlated so we may gain power, insight or both from a multivariate approach. For example, serum GGT is significantly correlated with total or LDL cholesterol, HDL (inversely) and particularly with triglycerides [17, 18]. Also, GGT is significantly correlated with other liver enzymes AST and ALT [17, 19]. Serum triglyceride is correlated with the liver enzymes  and uric acid and also associated with cardiovascular risk.
The importance of genetic variation has been shown previously through univariate analyses of serum lipids , uric acid [21–23], GGT , ALT  and AST [17, 24], BCHE , ferritin  and for CRP [27, 28]. Nevertheless, little is known about common genetic influences on these variables and joint analysis may reveal whether the same gene influences multiple traits.
The aim of our study is to identify genes and regions associated with multiple biochemical traits related to cardiovascular risk, type 2 diabetes or metabolic syndrome. We used a recently described multivariate association test  to perform genome-wide association analysis. This approach was used initially to screen for multivariate trait-SNP association using a subset of unrelated individuals. To confirm findings from the multivariate test, univariate association tests were conducted making use of the full dataset by including all family members.
Biochemical traits were measured in serum samples from twins and their families, and genome-wide SNP markers were genotyped. The study participants comprise:
(1) Adolescent twins and their non-twin siblings living in south-east Queensland (Australia) who had participated in the Brisbane Longitudinal Twin study [29–32]. Full details are described in Middelberg et al. . A total of 2548 participants (1317 females and 1231 males; mean age of 14.8 years) were genotyped.
(2) Adult twins consisting of twins and their family members who participated in studies of: (i) alcohol and nicotine dependence and metabolic risk for alcoholic liver disease ; (ii) anxiety and depression ; and (iii) endometriosis . A total of 9145 individuals (5703 females and 3442 males; mean age of 46.2 years) were genotyped.
Combining these studies, 20,230 individuals had biochemical measurements and 11,683 (from 4986 families) had both genotype and phenotype data. Out of the 11,683, there are only 1483 (from 1015 families) who had data for all the 13 traits. Where multiple measurements of the same trait in an individual were available, the average of the values was used.
For each of these studies, participants (and, for subjects aged < 18 years, their parents) gave informed consent to the questionnaire, interview, and blood collection, and all studies were approved by the QIMR Human Research Ethics Committee.
Serum was separated from the blood samples and stored at -70°C until analyzed. Serum cholesterol, HDL cholesterol, triglycerides, BCHE, glucose, uric acid, ferritin, CRP, AST, ALT and GGT were measured using Roche methods on a Roche 917 or Modular P analyzer (Roche Diagnostics, Basel, Switzerland). LDL cholesterol was calculated using the Friedewald equation. Insulin was measured on an Abbott Architect. BMI was calculated from measured or self-reported weight and height for the adults and from measured weight and height for the adolescents.
DNA was extracted from blood samples using standard methods and genotyped with Illumina 610K, 317K or 370K chips at CIDR or deCODE Genetics. Data cleaning for SNP genotypes included checking the expected relationships between individual family members and resolving Mendelian errors . Imputed genotypes for non-typed HapMap SNPs were generated using MACH1.0 (http://www.sph.umich.edu/csg/abecasis/mach/index.html) [38–40] program. Any imputed SNP which had r 2 ≥ 0.3 was included in the genotype data.
Distributions of all biochemical variables were examined. Serum AST, ALT and GGT, CRP, triglycerides and BMI were log-transformed. For each trait, individuals who were more than five standard deviations from the mean of that trait were excluded. Results for glucose and insulin in adults were adjusted for fasting time based on the reported time of last meal and time of blood collection. Prior to genetic association analysis, the variables were also adjusted for the effects of age, squared age (age2), sex, sex × age and sex × age2. Standardized residuals were obtained and used in the association analysis. All data pre-processing and descriptive analyses were done using STATA version 7.0  and SPSS version 17.0.2 (Mar 11, 2009). Multivariate association analysis was performed using the PLINK (v1.07) implementation of the multivariate test described by Ferreira and Purcell . This test is computationally too intensive when applied to family data to be efficient for genome-wide analysis. Therefore the analysis was performed in two stages. First, we selected one individual per family (using the person with data for the greatest number of phenotypes) from each of the 4986 families and applied the multivariate test as a screening tool. Next, for each locus with a multivariate p-value of less than 5 × 10-8, we identified the traits that showed evidence for association with that locus (that is, with a canonical correlation weight > |0.2|) and confirmed that specific trait-SNP association with a univariate association test using all relatives for each family. The univariate association test was performed using "fastassoc" in MERLIN 1.1.2  which takes the average of two results in MZ twin pairs.
Means and standard deviations of all the traits for males and females in adolescent and adult genotyped cohorts are listed in Additional file 1, Table S1. Generally, the means of the biochemical traits are lower in the adolescents than the adults, as expected. Phenotypic correlations between each pair of age-corrected traits separately for males and females in the combined sample are shown in Additional file 1, Table S2. The strongest correlations (r > 0.5) observed in males were between glucose and insulin (0.53), between AST and ALT (0.66) and between GGT and ALT (0.57). In the females, the strongest correlations observed were between glucose and insulin (0.59), between BCHE and glucose (-0.59), between BCHE and CRP (-0.53) and between AST and ALT (0.63). Given that most of the other pair-wise correlations (Additional file 1, Table S2) are low to moderate (r < 0.3), the multivariate approach is expected to provide comparable or slightly improved power to detect pleiotropic loci when compared to univariate analysis followed by correction for the number of traits tested .
Genome-wide association analyses
The multivariate analysis identified a total of 766 SNPs in 11 independent (r2 < 0.1) loci associated with biochemical traits with a p-value of less than 5 × 10-8 (Table 1 and Figure 1). Of these, there are eight loci potentially associated with more than one trait (Table 1). Three loci (on chromosomes 8, 12 and 19) showed strong or close to genome-wide significant evidence of associations with more than one trait in the all-subject univariate analyses.
The most strongly associated SNP at the chromosome 8p21.3 locus was rs17091905 (multivariate p = 2.8 × 10-13). HDL, CRP, triglycerides and BCHE had trait loadings of greater than |0.2|. To confirm the multivariate result, we individually tested each of these traits using a univariate test in the full sample of 11,683 individuals. The univariate tests confirmed the association with HDL (p = 5.7 × 10-12) and triglycerides (p = 5.1 × 10-15) but not at genome-wide significance for CRP (p = 0.008) and non-significant (p = 0.069) for BCHE. This variant is in strong or partial LD with previously-reported variants for HDL or triglycerides [43–45] (Additional file 1, Table S3).
The second variant rs3213545 (multivariate p-value = 3.9 × 10-14) which is located on chromosome 12q24.2 (OASL) was confirmed to be significantly associated with GGT (p = 3.6 × 10-15)  and also showed moderately strong significance for LDL (p = 2.9 × 10-5) and CRP (p = 8.8 × 10-5) (Table 1).
The third variant was rs2075650 (multivariate p-value = 5.7 × 10-10) located on chromosome 19q13.32 (TOMM40/APOE-C1-C2-C4 gene cluster) where HDL, LDL, CRP and triglycerides had trait loadings of greater than |0.2|. Significant univariate associations were observed for LDL (p = 1.6 × 10-14) and CRP (p = 4.2 × 10-8) and close to genome-wide significant univariate associations were seen for HDL (p = 8.1 × 10-8) and triglycerides (p = 9.6 × 10-7) (Table 1 and Figure 2). This SNP has previously been reported to be associated with LDL , LDL buoyancy  and CRP , and there is an association between LDL (or TG or HDL) and rs4420638 which is in partial LD (r2 = 0.4) with this SNP (Additional file 1, Table S3).
To determine whether there are any further unreported genes/regions to be detected by multivariate analysis, a lower p-value threshold of multivariate p < 9 × 10-5 was used. No new loci were found but a further six previously reported loci were replicated as listed in Additional file 1, Table S4.
The Q-Q plot from multivariate analysis was also examined closely to determine whether there are any excess association signals detected by multivariate analysis which have not already been detected by univariate analysis. SNPs that were found in significant regions (genes) in univariate analyses were removed (Figure 3) from the plot. The Q-Q plot with excluded SNPs showed that there is no excess of significant p-values hence indicating there are additional loci that have not already been detected by univariate analysis.
Examination of the directions of the allelic effects on the different phenotypes showed unexpected results. At LPL on chromosome 8, the minor allele A at rs17091905 increased HDL-C and decreased triglycerides, but the direction of the nominally significant effect on CRP was to increase it. At the chromosome 12 locus the minor allele A at rs3213545 tended to increase LDL-C but it significantly decreased GGT and tended to decrease CRP. Similarly at the chromosome 19 locus, the effect of the minor allele (G for rs2075650) was to increase LDL-C and triglycerides and to decrease HDL-C, consistent with an increase in cardiovascular risk, but to decrease CRP, again suggesting opposite allelic effects on the markers of different aspects of cardiovascular risk.
We have applied a multivariate approach to identify variants associated with more than one trait, initially using 4986 unrelated individuals across 13 biochemical traits. Univariate testing of the significant or near-significant loci, on the full sample of 11,683 individuals, was then used to confirm these findings. We are interested firstly in the usefulness of multivariate analysis as a substitute for the more laborious and potentially less powerful approach of conducting multiple univariate analyses and comparing the results, and secondly in the details of the loci which are found to have effects on multiple variables in our data.
Testing one individual per family identified three known loci that were significantly or near-significantly associated with more than one trait, and replicated 11 loci in previously published genes that that passed a genome-wide threshold of 5 × 10-8 for single variables (Table 1). When a lower genome-wide threshold (p < 9 × 10-5) was used, a further six published loci were also identified (Additional file 1, Table S4). The three loci in previous publications using univariate association analysis (highlighted in Table 1) had evidence of significant or close to significant associations with more than one trait in our data, hence indicating benefits of detecting pleiotropic loci in multivariate analysis.
We have identified polymorphisms showing strong evidence of allelic associations with HDL and triglycerides on chromosome 8 (LPL gene MIM 609708); with GGT and possibly LDL and CRP on chromosome 12 (OASL gene MIM 603281); and with HDL and LDL and possibly CRP and triglycerides on chromosome 19 (TOMM40 (MIM 608061) /APOE (MIM 107741)-C1 (MIM 107710)-C2 (MIM 608083)-C4 (MIM 600745) gene cluster). Each gene has been previously recognised in genome-wide association studies concentrating on a few of these variables . The function of these genes is reasonably well-established. LPL plays a key role in lipid metabolism and is responsible for hydrolysis of triglyceride molecules present in circulating lipoprotein. APOE and APOC genes also play a key role in lipid metabolism and cholesterol transport by helping to stabilise and solubilize lipoproteins as they circulate in the blood [50, 51]. Both LPL and APOE polymorphisms have been found to be significantly associated with increases in LDL and decreases in HDL . The functional connection between the OASL gene (2',5'-oligoadenylate synthetase-like, also known as "thyroid hormone receptor interactor") and these phenotypes is unclear. However, nearby genes in linkage disequilibrium with the lead SNP in OASL include HNF1A and c12orf43. HNF1A is expressed in liver, kidney and endocrine pancreas and regulates a number of genes involved in lipoprotein metabolism including apolipoproteins, cholesterol synthesis enzymes and bile acid transporters . HNF1A also has allelic associations with type 2 diabetes , CRP [55–57] and coronary heart disease . The findings for the lipids, in particular, were similar to those obtained in previous genetic association studies on general population. However the relationships between inflammation (as presumptively measured by CRP) and the traits associated with obesity and cardiovascular risk are of particular interest. CRP was significantly, though not always strongly, correlated with each of the other traits at the phenotypic level and it also showed up in the multivariate association findings.
The multivariate approach helps us to understand the connections between variables. For example, for rs2075650 on chromosome 19, the multivariate approach suggested LDL, CRP, HDL and triglycerides might be associated with this particular SNP. Although the effects on LDL, HDL and triglycerides are consistent with what was expected (that is, the LDL effect is inversely associated with the HDL effect, positively associated with the triglycerides effect, and the HDL effect is inversely associated with that on triglycerides), the effects on CRP are contrary to expectation. The effect direction is opposite to those for LDL and triglycerides, and the same as that for HDL. This suggests that the alleles or haplotype which have risk-increasing effects on lipids have a potentially protective effect on CRP and (so far as the effect on CRP is reflecting the degree of inflammation) on the inflammatory process. The effect estimates of LDL for rs2075650 obtained in our study were similar to obtained by Aulchenko et al. . The effect of rs2075650 [G] on LDL was estimated as 0.160 ± 0.018 by Aulchenko et al. and 0.153 ± 0.020 in our analysis. The effect estimates of CRP for this SNP was not available from previous study for comparison. In addition, it was interesting to observe that rs4429638 which is in partial LD (r2 = 0.4) with rs2075650 has allelic effects in the opposite direction on LDL [52, 59–61] and CRP . Similarly on chromosome 12, rs3213545 affects LDL and CRP (and GGT) but not HDL or triglycerides. Again, the allelic effects on LDL and CRP are in opposite directions. This shows the usefulness of the multivariate approach to help understand the connections between several trait-SNP associations, which can then be modelled and evaluated in more detail.
Our study differs from previous investigations as it examines a large number of correlated biochemical traits, initially using unrelated individuals and following up the findings in other members of the families. It confirms some published associations and identifies new ones. As our cohort consisted of adolescents and adults, results from adults, adolescents and combined (adults and adolescents) cohorts were examined and compared. Because of the larger number of adults studied, results from adults were similar to the combined data. Results from the adolescents were not notably different from the combined data.
One main limitation of our approach is that only a subset of the data (from unrelated subjects) can be used for the initial multivariate analysis. Although it would add power, it is too computationally intensive to use all the available data (that is, taking account of the family structure) in genome-wide multivariate analysis. Although a subset of data was used, the method applied in our study was very efficient and easy to perform. A more specific limitation in our data is that the glucose and insulin measurements were not made on fasting blood samples. In adults, we made adjustments for the time since the last meal but in adolescents we had to rely on the fact they were seen at the same time of day and blood was taken around three hours after the expected time of breakfast.
Another set of limitations is related to the use of biomarkers of risk or, for CRP, of systemic inflammation. It seems that some loci may affect HDL-C or triglycerides without affecting cardiovascular risk , and it is possible that some loci might affect serum CRP without affecting inflammation. Nevertheless the divergence between allelic effects on risk factors deserves further examination.
Our study demonstrated that it is useful to examine multiple phenotypes jointly in order to better understand the connections between them and to make the distinction between common and unique genetic effects. Our efficient approach (a combination of multivariate and univariate analysis) was able to identify three possible loci that might affect multiple traits, and validated 17 loci that have previously been reported. It highlighted anomalous effects on CRP, which is increasingly recognised as a marker of cardiovascular risk as well as of inflammation. Confirmation and extension of our findings will require studies which measure multiple phenotypes in each genotyped subject, and will benefit from combination of data from multiple studies to achieve sufficient power.
Amos CI, Elston RC, Bonney GE, Keats BJ, Berenson GS: A multivariate method for detecting genetic linkage, with application to a pedigree with an adverse lipoprotein phenotype. Am J Hum Genet. 1990, 47 (2): 247-254.
Almasy L, Dyer TD, Blangero J: Bivariate quantitative trait linkage analysis: pleiotropy versus co-incident linkages. Genet Epidemiol. 1997, 14 (6): 953-958. 10.1002/(SICI)1098-2272(1997)14:6<953::AID-GEPI65>3.0.CO;2-K.
Williams JT, Van Eerdewegh P, Almasy L, Blangero J: Joint multipoint linkage analysis of multivariate qualitative and quantitative traits. I. Likelihood formulation and simulation results. Am J Hum Genet. 1999, 65 (4): 1134-1147. 10.1086/302570.
Jiang C, Zeng ZB: Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics. 1995, 140 (3): 1111-1127.
Liu J, Pei Y, Papasian CJ, Deng HW: Bivariate association analyses for the mixture of continuous and binary traits with the use of extended generalized estimating equations. Genet Epidemiol. 2009, 33 (3): 217-227. 10.1002/gepi.20372.
Lange C, Silverman EK, Xu X, Weiss ST, Laird NM: A multivariate family-based association test using generalized estimating equations: FBAT-GEE. Biostatistics. 2003, 4 (2): 195-206. 10.1093/biostatistics/4.2.195.
Ferreira MA, Purcell SM: A multivariate test of association. Bioinformatics. 2009, 25 (1): 132-133. 10.1093/bioinformatics/btn563.
Ridker PM, Hennekens CH, Buring JE, Rifai N: C-reactive protein and other markers of inflammation in the prediction of cardiovascular disease in women. N Engl J Med. 2000, 342 (12): 836-843. 10.1056/NEJM200003233421202.
Nakanishi N, Nishina K, Li W, Sato M, Suzuki K, Tatara K: Serum gamma-glutamyltransferase and development of impaired fasting glucose or type 2 diabetes in middle-aged Japanese men. J Intern Med. 2003, 254 (3): 287-295. 10.1046/j.1365-2796.2003.01198.x.
Cutrin JC, Zingaro B, Camandola S, Boveris A, Pompella A, Poli G: Contribution of gamma glutamyl transpeptidase to oxidative damage of ischemic rat kidney. Kidney Int. 2000, 57 (2): 526-533. 10.1046/j.1523-1755.2000.00871.x.
Lee DH, Silventoinen K, Hu G, Jacobs DR, Jousilahti P, Sundvall J, Tuomilehto J: Serum gamma-glutamyltransferase predicts non-fatal myocardial infarction and fatal coronary heart disease among 28,838 middle-aged men and women. Eur Heart J. 2006, 27 (18): 2170-2176. 10.1093/eurheartj/ehl086.
Calderon-Margalit R, Adler B, Abramson JH, Gofin J, Kark JD: Butyrylcholinesterase activity, cardiovascular risk factors, and mortality in middle-aged and elderly men and women in Jerusalem. Clin Chem. 2006, 52 (5): 845-852. 10.1373/clinchem.2005.059857.
Iwasaki T, Yoneda M, Nakajima A, Terauchi Y: Serum butyrylcholinesterase is strongly associated with adiposity, the serum lipid profile and insulin resistance. Intern Med. 2007, 46 (19): 1633-1639. 10.2169/internalmedicine.46.0049.
Milman N, Kirchhoff M: Relationship between serum ferritin and risk factors for ischaemic heart disease in 2235 Danes aged 30-60 years. J Intern Med. 1999, 245 (5): 423-433. 10.1046/j.1365-2796.1999.00465.x.
Tomita M, Mizuno S, Yamanaka H, Hosoda Y, Sakuma K, Matuoka Y, Odaka M, Yamaguchi M, Yosida H, Morisawa H, et al: Does hyperuricemia affect mortality? A prospective cohort study of Japanese male workers. J Epidemiol. 2000, 10 (6): 403-409.
Liese AD, Hense HW, Lowel H, Doring A, Tietze M, Keil U: Association of serum uric acid with all-cause and cardiovascular disease mortality and incident myocardial infarction in the MONICA Augsburg cohort. World Health Organization Monitoring Trends and Determinants in Cardiovascular Diseases. Epidemiology. 1999, 10 (4): 391-397. 10.1097/00001648-199907000-00009.
Whitfield JB, Zhu G, Nestler JE, Heath AC, Martin NG: Genetic covariation between serum gamma-glutamyltransferase activity and cardiovascular risk factors. Clin Chem. 2002, 48 (9): 1426-1431.
Whitfield JB: Gamma glutamyl transferase. Crit Rev Clin Lab Sci. 2001, 38 (4): 263-355. 10.1080/20014091084227.
Arndt V, Brenner H, Rothenbacher D, Zschenderlein B, Fraisse E, Fliedner TM: Elevated liver enzyme activity in construction workers: prevalence and impact on early retirement and all-cause mortality. Int Arch Occup Environ Health. 1998, 71 (6): 405-412. 10.1007/s004200050299.
Beekman M, Heijmans BT, Martin NG, Pedersen NL, Whitfield JB, DeFaire U, van Baal GC, Snieder H, Vogler GP, Slagboom PE, et al: Heritabilities of apolipoprotein and lipid levels in three countries. Twin Res. 2002, 5 (2): 87-97. 10.1375/1369052022956.
Friedlander Y, Kark JD, Stein Y: Family resemblance for serum uric acid in a Jerusalem sample of families. Hum Genet. 1988, 79 (1): 58-63. 10.1007/BF00291711.
Whitfield JB, Martin NG: Plasma lipids in twins. Environmental and genetic influences. Atherosclerosis. 1983, 48 (3): 265-277. 10.1016/0021-9150(83)90044-8.
Kalousdian S, Fabsitz R, Havlik R, Christian J, Rosenman R: Heritability of clinical chemistries in an older twin cohort: the NHLBI Twin Study. Genet Epidemiol. 1987, 4 (1): 1-11. 10.1002/gepi.1370040102.
Bathum L, Petersen HC, Rosholm JU, Hyltoft Petersen P, Vaupel J, Christensen K: Evidence for a substantial genetic influence on biochemical liver function tests: results from a population-based Danish twin study. Clin Chem. 2001, 47 (1): 81-87.
Valle A, O'Connor DT, Taylor P, Zhu G, Montgomery GW, Slagboom PE, Martin NG, Whitfield JB: Butyrylcholinesterase: association with the metabolic syndrome and identification of 2 gene loci affecting activity. Clin Chem. 2006, 52 (6): 1014-1020. 10.1373/clinchem.2005.065052.
Njajou OT, Alizadeh BZ, Aulchenko Y, Zillikens MC, Pols HA, Oostra BA, Swinkels DW, van Duijn CM: Heritability of serum iron, ferritin and transferrin saturation in a genetically isolated population, the Erasmus Rucphen Family (ERF) Study. Hum Hered. 2006, 61 (4): 222-228. 10.1159/000094777.
Su S, Snieder H, Miller AH, Ritchie J, Bremner JD, Goldberg J, Dai J, Jones L, Murrah NV, Zhao J, et al: Genetic and environmental influences on systemic markers of inflammation in middle-aged male twins. Atherosclerosis. 2008, 200 (1): 213-220. 10.1016/j.atherosclerosis.2007.12.009.
Lange LA, Burdon K, Langefeld CD, Liu Y, Beck SR, Rich SS, Freedman BI, Brosnihan KB, Herrington DM, Wagenknecht LE, et al: Heritability and expression of C-reactive protein in type 2 diabetes in the Diabetes Heart Study. Ann Hum Genet. 2006, 70 (Pt 6): 717-725.
Zhu G, Duffy DL, Eldridge A, Grace M, Mayne C, O'Gorman L, Aitken JF, Neale MC, Hayward NK, Green AC, et al: A major quantitative-trait locus for mole density is linked to the familial melanoma gene CDKN2A: a maximum-likelihood combined linkage and association analysis in twins and their sibs. Am J Hum Genet. 1999, 65 (2): 483-492. 10.1086/302494.
Wright M, De Geus E, Ando J, Luciano M, Posthuma D, Ono Y, Hansell N, Van Baal C, Hiraishi K, Hasegawa T, et al: Genetics of cognition: outline of a collaborative twin study. Twin Res. 2001, 4 (1): 48-56. 10.1375/1369052012146.
McGregor B, Pfitzner J, Zhu G, Grace M, Eldridge A, Pearson J, Mayne C, Aitken JF, Green AC, Martin NG: Genetic and environmental contributions to size, color, shape, and other characteristics of melanocytic naevi in a sample of adolescent twins. Genet Epidemiol. 1999, 16 (1): 40-53. 10.1002/(SICI)1098-2272(1999)16:1<40::AID-GEPI4>3.0.CO;2-1.
Wright MJ, Martin NG: Brisbane Adolescent Twin Study: Outline of study methods and research projects. Australian Journal of Psychology. 2004, 56 (2): 65-78. 10.1080/00049530410001734865.
Middelberg RP, Martin NG, Whitfield JB: A longitudinal genetic study of plasma lipids in adolescent twins. Twin Res Hum Genet. 2007, 10 (1): 127-135. 10.1375/twin.10.1.127.
Heath AC, Whitfield JB, Martin NG, Pergadia ML, Goate AM, Lind PA, McEvoy BP, Schrage AJ, Grant JD, Chou YL, et al: A Quantitative-Trait Genome-Wide Association Study of Alcoholism Risk in the Community: Findings and Implications. Biol Psychiatry. 2011, Epub April 27
Sullivan PF, de Geus EJ, Willemsen G, James MR, Smit JH, Zandbelt T, Arolt V, Baune BT, Blackwood D, Cichon S, et al: Genome-wide association for major depressive disorder: a possible role for the presynaptic protein piccolo. Mol Psychiatry. 2009, 14 (4): 359-375. 10.1038/mp.2008.125.
Painter JN, Anderson CA, Nyholt DR, Macgregor S, Lin J, Lee SH, Lambert A, Zhao ZZ, Roseman F, Guo Q, et al: Genome-wide association study identifies a locus at 7p15.2 associated with endometriosis. Nat Genet. 2011, 43 (1): 51-54. 10.1038/ng.731.
Medland SE, Nyholt DR, Painter JN, McEvoy BP, McRae AF, Zhu G, Gordon SD, Ferreira MA, Wright MJ, Henders AK, et al: Common variants in the trichohyalin gene are associated with straight hair in Europeans. Am J Hum Genet. 2009, 85 (5): 750-755. 10.1016/j.ajhg.2009.10.009.
Ellinghaus D, Schreiber S, Franke A, Nothnagel M: Current software for genotype imputation. Hum Genomics. 2009, 3 (4): 371-380.
Marchini J, Howie B: Genotype imputation for genome-wide association studies. Nat Rev Genet. 2010, 11 (7): 499-511. 10.1038/nrg2796.
Nothnagel M, Ellinghaus D, Schreiber S, Krawczak M, Franke A: A comprehensive evaluation of SNP genotype imputation. Hum Genet. 2009, 125 (2): 163-171. 10.1007/s00439-008-0606-5.
StataCorp: Stata Statistical Software: Release 7.0. 1997, College Station, Tex: StataCorp
Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin--rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2002, 30 (1): 97-101. 10.1038/ng786.
Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, Pramstaller PP, Penninx BW, Janssens AC, Wilson JF, Spector T, et al: Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet. 2009, 41 (1): 47-55. 10.1038/ng.269.
Talmud PJ, Drenos F, Shah S, Shah T, Palmen J, Verzilli C, Gaunt TR, Pallas J, Lovering R, Li K, et al: Gene-centric association signals for lipids and apolipoproteins identified via the HumanCVD BeadChip. Am J Hum Genet. 2009, 85 (5): 628-642. 10.1016/j.ajhg.2009.10.014.
Ma L, Yang J, Runesha HB, Tanaka T, Ferrucci L, Bandinelli S, Da Y: Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham heart study data. BMC Med Genet. 2010, 11: 55-
Yuan X, Waterworth D, Perry JR, Lim N, Song K, Chambers JC, Zhang W, Vollenweider P, Stirnadel H, Johnson T, et al: Population-based genome-wide association studies reveal six loci influencing plasma levels of liver enzymes. Am J Hum Genet. 2008, 83 (4): 520-528. 10.1016/j.ajhg.2008.09.012.
Ronald J, Rajagopalan R, Ranchalis JE, Marshall JK, Hatsukami TS, Heagerty PJ, Jarvik GP: Analysis of recently identified dyslipidemia alleles reveals two loci that contribute to risk for carotid artery disease. Lipids Health Dis. 2009, 8: 52-10.1186/1476-511X-8-52.
Reiner AP, Barber MJ, Guan Y, Ridker PM, Lange LA, Chasman DI, Walston JD, Cooper GM, Jenny NS, Rieder MJ, et al: Polymorphisms of the HNF1A gene encoding hepatocyte nuclear factor-1 alpha are associated with C-reactive protein. Am J Hum Genet. 2008, 82 (5): 1193-1201. 10.1016/j.ajhg.2008.03.017.
Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ, et al: Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010, 466 (7307): 707-713. 10.1038/nature09270.
Mahley RW: Apolipoprotein E: cholesterol transport protein with expanding role in cell biology. Science. 1988, 240 (4852): 622-630. 10.1126/science.3283935.
Eichner JE, Dunn ST, Perveen G, Thompson DM, Stewart KE, Stroehla BC: Apolipoprotein E polymorphism and cardiovascular disease: a HuGE review. Am J Epidemiol. 2002, 155 (6): 487-495. 10.1093/aje/155.6.487.
Kathiresan S, Melander O, Guiducci C, Surti A, Burtt NP, Rieder MJ, Cooper GM, Roos C, Voight BF, Havulinna AS, et al: Six new loci associated with blood low-density lipoprotein cholesterol, high-density lipoprotein cholesterol or triglycerides in humans. Nat Genet. 2008, 40 (2): 189-197. 10.1038/ng.75.
Odom DT, Zizlsperger N, Gordon DB, Bell GW, Rinaldi NJ, Murray HL, Volkert TL, Schreiber J, Rolfe PA, Gifford DK, et al: Control of pancreas and liver gene expression by HNF transcription factors. Science. 2004, 303 (5662): 1378-1381. 10.1126/science.1089769.
Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G, et al: Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet. 2010, 42 (7): 579-589. 10.1038/ng.609.
Okada Y, Takahashi A, Ohmiya H, Kumasaka N, Kamatani Y, Hosono N, Tsunoda T, Matsuda K, Tanaka T, Kubo M, et al: Genome-wide association study for C-reactive protein levels identified pleiotropic associations in the IL6 locus. Hum Mol Genet. 2011, 20 (6): 1224-1231. 10.1093/hmg/ddq551.
Elliott P, Chambers JC, Zhang W, Clarke R, Hopewell JC, Peden JF, Erdmann J, Braund P, Engert JC, Bennett D, et al: Genetic Loci associated with C-reactive protein levels and risk of coronary heart disease. Jama. 2009, 302 (1): 37-48. 10.1001/jama.2009.954.
Dehghan A, Dupuis J, Barbalic M, Bis JC, Eiriksdottir G, Lu C, Pellikka N, Wallaschofski H, Kettunen J, Henneman P, et al: Meta-Analysis of Genome-Wide Association Studies in > 80 000 Subjects Identifies Multiple Loci for C-Reactive Protein Levels. Circulation. 2011, 731-738.
Erdmann J, Grosshennig A, Braund PS, Konig IR, Hengstenberg C, Hall AS, Linsel-Nitschke P, Kathiresan S, Wright B, Tregouet DA, et al: New susceptibility locus for coronary artery disease on chromosome 3q22.3. Nat Genet. 2009, 41 (3): 280-282. 10.1038/ng.307.
Sandhu MS, Waterworth DM, Debenham SL, Wheeler E, Papadakis K, Zhao JH, Song K, Yuan X, Johnson T, Ashford S, et al: LDL-cholesterol concentrations: a genome-wide association study. Lancet. 2008, 371 (9611): 483-491. 10.1016/S0140-6736(08)60208-1.
Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, Clarke R, Heath SC, Timpson NJ, Najjar SS, Stringham HM, et al: Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet. 2008, 40 (2): 161-169. 10.1038/ng.76.
Waterworth DM, Ricketts SL, Song K, Chen L, Zhao JH, Ripatti S, Aulchenko YS, Zhang W, Yuan X, Lim N, et al: Genetic Variants Influencing Circulating Lipid Levels and Risk of Coronary Artery Disease. Arterioscler Thromb Vasc Biol. 2010, 30 (11): 2264-2276. 10.1161/ATVBAHA.109.201020.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/12/123/prepub
We are grateful to the twins and their families for their generous participation in these studies. We would like to thank Ann Eldridge and Marlene Grace for the collection of data; Leanne Wallace for sample handling and preparation for the serum biochemistry and genotyping; Lisa Bowdler, Steven Crooks, Melinda Richter and Sara Smith for sample processing and preparation for the GWAS genotyping; Harry Beeby for IT support; and the following people in the QIMR GWAS group: Scott Gordon, Dale Nyholt, Stuart Macgregor and Naomi Wray for generating and checking of integrated genome-wide typing data and Sarah Medland for imputing the data. Financial support was provided from the National Health and Medical Research Council (241944, 389875), NIH (AA07535, AA13320, AA13321, AA14041, AA11998, AA17688, DA012854, DA019951) and Wellcome Trust for funding (WT084766/Z/08/Z). GWM is supported by an NHMRC Fellowship (619667) and RPSM is supported by NHMRC Training Fellowship (611512).
The authors declare that they have no competing interests.
RPSM analysed data, interpreted and prepared the manuscript. MARF and JBW contributed to data interpretation and manuscript preparation. AKH and GWM coordinated genotyping studies. ACH, PAFM and NGM coordinated the phenotype studies. All authors read and approved this manuscript.
Electronic supplementary material
Additional file 1: Additional tables and references. Additional Tables S1-S4 and references cited in these tables. Table S1 gives a descriptive statistics for males and females in the two genotyped cohorts. Table S2 gives the phenotypic correlations between each pair of age-corrected traits separately for males and females in combined sample. Table S3 compares the multivariate and univariates results from our study with published results. Table S4 gives a summary of borderline significant associations from multivariate analysis. (PDF 170 KB)
About this article
Cite this article
Middelberg, R.P., Ferreira, M.A., Henders, A.K. et al. Genetic variants in LPL, OASL and TOMM40/APOE-C1-C2-C4 genes are associated with multiple cardiovascular-related traits. BMC Med Genet 12, 123 (2011). https://doi.org/10.1186/1471-2350-12-123
- Multivariate Approach
- Allelic Effect
- Allelic Association
- Biochemical Trait
- Multiple Phenotype