Analysis of Multiple Association Studies Provides Evidence of an Expression QTL Hub in Gene-Gene Interaction Network Affecting HDL Cholesterol Levels

Abstract

Epistasis has been suggested to underlie part of the missing heritability in genome-wide association studies. In this study, we first report an analysis of gene-gene interactions affecting HDL cholesterol (HDL-C) levels in a candidate gene study of 2,091 individuals with mixed dyslipidemia from a clinical trial. Two additional studies, the Atherosclerosis Risk in Communities study (ARIC; n = 9,713) and the Multi-Ethnic Study of Atherosclerosis (MESA; n = 2,685), were considered for replication. We identified a gene-gene interaction between rs1532085 and rs12980554 (P = 7.1×10−7) in their effect on HDL-C levels, which is significant after Bonferroni correction (Pc = 0.017) for the number of SNP pairs tested. The interaction successfully replicated in the ARIC study (P = 7.0×10−4; Pc = 0.02). Rs1532085, an expression QTL (eQTL) of LIPC, is one of the two SNPs involved in another, well-replicated gene-gene interaction underlying HDL-C levels. To further investigate the role of this eQTL SNP in gene-gene interactions affecting HDL-C, we tested in the ARIC study for interaction between this SNP and any other SNP genome-wide. We found the eQTL to be involved in a few suggestive interactions, one of which significantly replicated in MESA. Importantly, these gene-gene interactions, involving only rs1532085, explain an additional 1.4% variation of HDL-C, on top of the 0.65% explained by rs1532085 alone. LIPC plays a key role in the lipid metabolism pathway and it, and rs1532085 in particular, has been associated with HDL-C and other lipid levels. Collectively, we discovered several novel gene-gene interactions, all involving an eQTL of LIPC, thus suggesting a hub role of LIPC in the gene-gene interaction network that regulates HDL-C levels, which in turn raises the hypothesis that LIPC's contribution is largely via interactions with other lipid metabolism related genes.

About the Authors:

Li Ma

Current address: Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland, United States of America

Affiliation: Department of Biological Statistics and Computational Biology, Cornell University, Ithaca, New York, United States of America

Christie Ballantyne

Affiliation: Section of Cardiovascular Research, Department of Medicine, Baylor College of Medicine, Houston, Texas, United States of America

Ariel Brautbar

* E-mail: [email protected] (AK); [email protected] (AB)

Affiliations Section of Cardiovascular Research, Department of Medicine, Baylor College of Medicine, Houston, Texas, United States of America, Department of Medical Genetics, Marshfield Clinic, Marshfield, Wisconsin, United States of America

Alon Keinan

* E-mail: [email protected] (AK); [email protected] (AB)

Affiliation: Department of Biological Statistics and Computational Biology, Cornell University, Ithaca, New York, United States of America

Introduction

Genome-wide association studies (GWAS) have enabled the identification of thousands of single-nucleotide polymorphisms (SNPs) that are associated with complex human diseases and traits [1]. However, most of these associated SNPs exhibit small effect sizes and collectively explain only a minor fraction of heritable variation [2]. Epistasis, or gene-gene interactions, plays an important role underlying the genetic basis of complex diseases and traits [3], [4], and has been suggested to underlie some of the “missing heritability” in human GWAS [2], [5]. High-density lipoprotein cholesterol (HDL-C) is an important risk factor for cardiovascular disease and is highly heritable [6]. In the last few years, there has been a surge of GWAS of lipid and lipoprotein traits, which have identified hundreds of loci associated with the levels of HDL-C [1], [7]–[11]. However, the combination of all significantly associated loci in large-scale meta-analyses still explains only 9.9∼12.1% of the total variance in HDL-C levels [10], [11], while 54∼70% is heritable [6]. Several recent studies have re-analyzed GWAS datasets and reported a few gene-gene interactions associated with human HDL-C levels [9], [12]–[15].

While epistasis or gene-gene interactions have long been studied in plants and animals [4], [16], they have proven difficult to identify in human GWAS mainly due to the limited statistical power caused by small effect size, moderate sample size, and the multiple-testing correction burden inherent in the large number of interaction tests required [5], [12], [17]–[20]. Prior biological insight has been incorporated into analysis of epistasis, aiming to increase power by focusing only on a selected subset of genes and SNPs that are likely to be enriched for putative interactions. Such biological insight can include, among others, previous GWAS results, protein-protein interactions, and pathway information [12], [14], [20]–[22]. By testing interactions between genetic loci that are marginally associated with lipid traits, we recently detected and replicated an interaction between HMGCR and a locus upstream of LIPC in their effect on HDL-C levels in several GWAS from different ethnicities [12].

In this study, we follow a similar approach by testing interactions between a focused set of SNPs in a Randomized Candidate Gene (RCG) study, in which the genes have been targeted as related to lipid and lipoprotein levels [13], [23]. We identified and replicated one significant gene-gene interaction affecting HDL-C levels, which involves the exact same SNP as the one previously reported to be involved in a different gene-gene interaction underlying HDL-C levels [12]. Interestingly, this SNP (rs1532085) is an expression QTL (eQTL) of LIPC, a gene that plays a key role in the related pathway. To further investigate this SNP's role in gene-gene interactions affecting HDL-C, we tested for gene-gene interactions between it and all other, genome-wide SNPs in large-scale data sets from the Atherosclerosis Risk in Communities (ARIC) study [24] and the Multi-Ethnic Study of Atherosclerosis (MESA; [25]. Our results combined point to several suggestive gene-gene interactions all involving rs1532085 as one of the two interacting SNPs, two of which significantly replicated in independent datasets. We discuss the possibility of LIPC as a hub in the gene-gene interaction network underlying HDL-C levels and the biological context of such an interaction networks.

Materials and Methods

Study Descriptions

All work done in this paper was approved by Cornell institutional review boards (IRB) committee. The three studies we analyzed were also approved by their local IRB committees.

Randomized Candidate Gene (RCG) Study

We used data from a Randomized Candidate Gene (RCG) study, which utilized a randomized clinical trial (www.clinicaltrials.gov; NCT00300482, NCT00300456, NCT00300469, NCT00300430) to assess the responses to medical treatments in individuals with mixed dyslipidemia [23], [26]. In brief, men and women with triglyceride (TG) ≥150 mg/dl, HDL-C <40 mg/dl for men or <50 mg/dl for women, and low-density lipoprotein cholesterol (LDL-C) ≥130 mg/dl were considered for the study. Baseline levels of HDL-C, TG, Apolipoprotein A1 (APOA1), Apolipoprotein B (APOB), Apolipoprotein C-3 (APOC3) and other necessary covariates were measured at the beginning of the study. A targeted set of SNPs were genotyped in or near candidate genes related to TG, HDL-C, and APOC3 pathways [27]. SNPs identified in published GWAS to be associated with HDL-C and TG [7], [28] were further genotyped. In total, 350 SNPs were genotyped in 2,091 European American (EA) individuals. We excluded samples with call rate <90%, SNPs with a call rate <90%, as well as SNPs deviating from Hardy–Weinberg equilibrium (HWE) at P value <0.001. After quality control (QC), 304 SNPs were kept for interaction analysis. We performed a principal component analysis using all 304 SNPs and identified no apparent population stratification [23], [26], which was also supported by observing no inflation in the QQ-plot of the P values for the interaction tests (Figure 1A).

[Figure omitted. See PDF.]

Figure 1. QQ plot and Box plot of the interactions identified in the RCG study.

A) QQ plot for interaction tests between all pairs of 304 candidate SNPs in RCG. B) Boxplot of HDL-C levels for the nine genotype combinations of rs1532085 and rs12980554, the pair of SNPs that exhibit the most significant gene-gene interaction in this study. The interaction is dominated by the combination of AA genotype in both SNPs.

https://doi.org/10.1371/journal.pone.0092469.g001

Atherosclerosis Risk in Communities (ARIC) Study

The ARIC Study is a multi-center prospective study of atherosclerotic disease [24], which recruited European Americans and African Americans from four communities: Forsyth County, North Carolina; Jackson, Mississippi; suburban areas of Minneapolis, Minnesota; and Washington County, Maryland. A total of 15,792 individuals participated in the baseline examination in 1987–1989, with three triennial follow-up examinations. We analyzed Affymetrix 6.0 SNP array genotyping of 9,713 EA from the study, and for whom data for both HDL-C levels and covariates were available. These data was first used as a replication dataset. Following our initial discovery, these data was also used for epistasis hub analysis of SNP rs1532085. For the latter analysis, after standard QC described previously (Ma et al. 2012a), we first applied linkage disequilibrium (LD) based pruning on the genome-wide autosomal SNPs using PLINK (with settings, —indep-pairwise 50 5 0.5) [29], which resulted in 136,881 SNPs that are in approximate linkage equilibrium (pairwise r2>0.5). We then tested for pairwise gene-gene interaction in effect on HDL-C between rs1532085 and each of these 136,881 SNPs using a linear model approach described in the following.

Multi-Ethnic Study of Atherosclerosis (MESA)

MESA is a prospective cohort study recruiting individuals aged 45–84 years from six US locations (Baltimore, MD; Chicago, IL; Forsyth County, NC; Los Angeles County, CA; northern Manhattan, NY; and St. Paul, MN) [25]. MESA aimed to determine the characteristics of subclinical cardiovascular disease and its progression. Participants were enrolled between July 2000 and August 2002 and self-reported their race group as European American, African American, Hispanic, or Chinese American. We obtained Affymetrix 6.0 SNP array genotyping of the EA samples from the database of Genotypes and Phenotypes (dbGaP; MESA SHARe, downloaded in May 2011; [30]. We considered a sample of 2,685 European Americans, for which both genotypes and HDL-C measurements were available, as a replication dataset.

Testing of Statistical Interactions

We tested for interaction between pairs of SNPs using a previously-described linear model approach [12], [16]. Briefly, an F-test with four degrees of freedom was conducted to compare two models, one with marginal additive and dominance effects of the two SNPs and the other with these as well as and interactions between the two SNPs, while accounting for covariates. Hence, this test is one for pure interaction effect, on top and beyond any marginal effect of either or both SNPs. In order to increase power and false positive rate, we applied additional filtering by only considering SNP pairs for which the least frequent genotype combination, out of the 9 possible combinations, is carried by at least 10 individuals.

For each pair of SNPs, the two linear models were fitted as follows.where β0 denotes a vector of intercept and covariate effects as described above, ai and di denote the additive and dominance effects of SNP i, and iaa, iad, ida, and idd are the four interaction effects between the two SNPs, with a denoting additive effect and d dominance effect. Y is a quantitative trait of interest, Z0 consists of the covariates, including age, sex and body mass index, xi and zi are two indicator variables of the genotype of SNP i, and is a random error. R square value of each model was calculated and the difference of the two measures the effect of a pure gene-gene interaction, on top of any marginal effects of each SNP alone.

Results

An interaction between rs1532085 and rs12980554 affects HDL-C levels

We tested pairwise SNP×SNP interactions on HDL-C levels between 304 candidate SNPs in a sample of 2,091 European American (EA) individuals with mixed dyslipidemia from the RCG study. One interaction, between rs1532085 and rs12980554, was significantly associated with HDL-C levels (Figure 1A; P = 7.1×10−7; Bonferroni corrected Pc = 0.017). Importantly, rs1532085 is the exact same SNP involved in a recently-reported gene-gene interaction, though there is was in interaction with SNPs in the HMGCR gene [12]. While the average individual HDL-C level in the RCG study is 38.5±7.1 mg/dl, individuals that are homozygous for the minor allele of both rs1532085 and rs12980554 exhibit the highest HDL-C level of 43±5.9 mg/dl (Figure 1B).

After discovery of the gene-gene between rs1532085 and rs12980554, we further tested whether it underlies other lipid levels measured in the same study. We found it to be nominally associated with the levels of APOA1 (P = 2.2×10−6) and TG (P = 0.01), and almost associated with APOC3 levels (P = 0.08). Using a locus-based replication approach [12], we also successfully replicated this interaction on HDL-C levels in an independent European American sample from the ARIC study (P = 7.0×10−4; Pc = 0.02). Replication failed in the MESA study, which has a much smaller sample size. The replicated interaction in the ARIC study is between the same SNP upstream of LIPC, rs1532085, and another SNP in the second locus, rs2241589, which is proximate (22.5 kb) and in moderate LD (r2 = 0.1) with rs12980554. This replicated interaction is also nominally affecting the level of LDL-C (P = 0.04) and almost that of total cholesterol (P = 0.07) in ARIC. Both rs2241589 and rs12980554 are located in the KANK3–ANGPTL4 locus, which has been associated with HDL-C levels in European American populations [31], [32]. Interestingly, ANGPTL4 encodes a protein that inhibits hepatic lipase's triglyceride hydrolyzing activity [33].

Epistasis hub effect of rs1532085 on HDL-C

In a previous study [12], we reported an independent and well-replicated interaction on HDL-C levels of the same SNP (rs1532085) involved in the novel gene-gene interaction reported above. Rs1532085 has been shown to cis-regulate the expression of LIPC [10], [34], which encodes hepatic lipase that plays a key role in the metabolism of lipids, lipoproteins, and HDL in particular [35]. LIPC and rs1532085 have been associated with HDL-C levels [10]. Combining all these information, we hypothesize that LIPC serves as a potential hub of the epistatic interaction network regulating HDL-C levels, namely carrying out its role largely in interaction with other gene products.

To test this hypothesis, we tested interactions between rs1532085 and each of the other SNPs in the ARIC study in their effect on HDL-C levels. To reduce the number of dependent tests, we pruned the genome-wide autosomal SNPs based on LD [29], resulting in a subset of 136,881 independent SNPs for the interaction testing. We detected four suggestive independent interactions at a false discovery rate (FDR) of 0.25 (Table 1; Figure 2). While rs1532085 alone explains 0.65% of the total variation in HDL-C, each of these four interactions explains 0.33∼0.42% of the variation in HDL-C and combined explain an additional 1.4% variation in HDL-C for a total of 2.05% of variation in HDL-C levels. While none of the four interactions is significant by itself following a conservative Bonferroni correction, only one false positive is expected out of the four at the applied FDR of 0.25. One of the interactions, involving HMGCR (Table 1), has been previously reported and replicated [12]. We note that the other three interactions were not tested in that previous study since it has focused on candidate SNPs only. Hence, we attempted replicating these three interactions in an independent EA sample from MESA, and significantly replicated one of them, between rs1532085 and rs1867732 (Table 1; P = 0.03). The latter SNP is located in an intron of C5orf64, encoding an uncharacterized protein and as far as we know has not been associated with any lipid trait. Furthermore, while the interaction with another SNP, rs2901656, intronic in C1orf105, does not replicate, we note that C1orf105 has been shown to be co-expressed with C5orf64 [36], as well as with LDLR and APOB that are important players in the pathway of metabolism of lipids and lipoproteins [37]–[40].

[Figure omitted. See PDF.]

Figure 2. QQ plot for interaction tests between rs1532085 and genome-wide LD-pruned SNPs in the ARIC study.

The plot and analysis in main text point to a deviation of the P values for four independent SNPs, which are detailed in Table 1.

https://doi.org/10.1371/journal.pone.0092469.g002

[Figure omitted. See PDF.]

Table 1. Top SNPs (FDR = 0.25) interacting with rs1532085 in affecting HDL-C levels in ARIC.

https://doi.org/10.1371/journal.pone.0092469.t001

Discussion

In this study we have shown that rs1532085 takes part in three replicated gene-gene interactions that affect the levels of HDL-C and other lipids, with two additional interactions involving it discovered but not replicated. Beyond the statistical significance of each interaction, the fact that the same SNPs, rs1532085, has independently been discovered as part of the only significant gene-gene interaction in each of two datasets, interacting with a different SNP in each, underscores its important role in regulating HDL-C. Rs1532085 is located in the regulatory region of LIPC and is an eQTL of LIPC [34], i.e. it cis-regulates LIPC expression. LIPC is known to play an important role in the lipid and lipoprotein metabolism pathway [35]. Furthermore, LIPC, and rs1532085 in particular, has been associated with the level of HDL-C [8], [10], [41], total cholesterol [8], [10] and triglyceride [10], as well as with metabolic syndrome [42]. However, it explains only 0.65% of the variation in HDL-C levels, while together with its gene-gene interactions reported herein, 2.05% of that variation is explained.

One of the interactions of rs1532085, which we discovered in the RCG study and replicated in the ARIC study, is with rs2241589. The latter is intronic in KANK3, which is next to ANGPTL4. This interaction resulted in increased levels of HDL-C and APOA1, which points to increased inhibition of hepatic lipase by ANGPTL4 when the two minor variants are present simultaneously. This interaction may also affect reverse cholesterol transport (RCT), a process in which cholesterol is removed from macrophages in the arterial wall, transferred to the HDL particle, and back to the liver for reprocessing or elimination in bile. It is one of the endogenous mechanisms to reduce atherosclerosis via the HDL particle, which plays a role in RCT, and inhibition of its function may have an effect on atherosclerosis reduction [43]. LIPC, ANGPTL4, as well as HMGCR that is involved in our previously reported interaction with rs1532085, have all been shown to physically interact or co-express with LDLR, a low density lipoprotein receptor [44].

Hubs, namely nodes in a network with a high number of connections, have been shown to be omnipresent in various types of networks, including social networks [45], protein-protein interaction networks [46], [47], and genetic interaction networks [48]. Hub genes in a genetic network tend to be important in evolution [49], and are suggested to play key roles in the biological regulation and functions [48]. Epistasis hubs have been suggested to exist in humans [50], and we propose that LIPC regulation constitutes such a hub in the gene-gene network affecting HDL-C levels. It is perhaps no major surprise that if a hub exists in the gene-gene interaction network underlying the levels of HDL-C, it will be related to LIPC and its regulation. LIPC encodes for the enzyme hepatic lipase which has an important role in HDL-C metabolism and RCT. Hepatic lipase hydrolyzes the HDL's phospholipids and triglycerides making it susceptible for clearance by the kidney and resulting in a net effect of HDL-C reduction.

In conclusion, while tests for detecting gene-gene interactions are not as powerful as the routine detection of marginal associations, judicious strategies of targeting SNPs that have been previously associated with lipid levels and candidate genes more generally, allowed us to identify and replicate several gene-gene interactions underlying HDL-C and other lipid levels. All these interactions involve a regulatory variant that cis-regulate LIPC expression, which we propose is a potential hub in the functional gene-gene interaction network underlying HDL-C metabolism and perhaps the metabolism of other lipids and lipoproteins.

Author Contributions

Conceived and designed the experiments: LM AB AK. Analyzed the data: LM. Contributed reagents/materials/analysis tools: CB AB. Wrote the paper: LM AB AK.

Citation: Ma L, Ballantyne C, Brautbar A, Keinan A (2014) Analysis of Multiple Association Studies Provides Evidence of an Expression QTL Hub in Gene-Gene Interaction Network Affecting HDL Cholesterol Levels. PLoS ONE 9(3): e92469. https://doi.org/10.1371/journal.pone.0092469

References

1. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, et al. (2009) Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A 106: 9362–9367.

2. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, et al. (2009) Finding the missing heritability of complex diseases. Nature 461: 747–753.

3. Moore JH, Williams SM (2005) Traversing the conceptual divide between biological and statistical epistasis: systems biology and a more modern synthesis. Bioessays 27: 637–646.

4. Phillips PC (2008) Epistasis—the essential role of gene interactions in the structure and evolution of genetic systems. Nature Reviews Genetics 9: 855–867.

5. Zuk O, Hechter E, Sunyaev SR, Lander ES (2012) The mystery of missing heritability: Genetic interactions create phantom heritability. Proceedings of the National Academy of Sciences 109: 1193–1198.

6. Weiss LA, Pan L, Abney M, Ober C (2006) The sex-specific genetic architecture of quantitative traits in humans. Nat Genet 38: 218–222.

7. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, et al. (2008) Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nature genetics 40: 161–169.

8. Aulchenko YS, Ripatti S, Lindqvist I, Boomsma D, Heid IM, et al. (2009) Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts. Nat Genet 41: 47–55.

9. Ma L, Yang J, Runesha HB, Tanaka T, Ferrucci L, et al. (2010) Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham Heart Study data. BMC Medical Genetics 11: 55.

10. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, et al. (2010) Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466: 707–713.

11. Asselbergs FW, Guo YR, van Iperen EPA, Sivapalaratnam S, Tragante V, et al. (2012) Large-Scale Gene-Centric Meta-analysis across 32 Studies Identifies Multiple Lipid Loci. American Journal of Human Genetics 91: 823–838.

12. Ma L, Brautbar A, Boerwinkle E, Sing CF, Clark AG, et al.. (2012) Knowledge-Driven Analysis Identifies a Gene-Gene Interaction Affecting High-Density Lipoprotein Cholesterol Levels in Multi-Ethnic Populations. PLoS genetics 8..

13. Ma L, Ballantyne CM, Belmont JW, Keinan A, Brautbar A (2012) Interaction between SNPs in the RXRA and near ANGPTL3 gene region inhibit apolipoprotein B reduction following statin-fenofibric acid therapy in individuals with mixed dyslipidemia. Journal of Lipid Research.

14. Turner SD, Berg RL, Linneman JG, Peissig PL, Crawford DC, et al. (2011) Knowledge-driven multi-locus analysis reveals gene-gene interactions influencing HDL cholesterol level in two independent EMR-linked biobanks. PLoS One 6: e19586.

15. Maxwell TJ, Ballantyne CM, Cheverud JM, Guild CS, Ndumele CE, et al.. (2013) APOE Modulates the Correlation Between Triglycerides, Cholesterol, and CHD Through Pleiotropy and Gene-by-Gene Interactions. Genetics.

16. Cordell HJ (2009) Detecting gene-gene interactions that underlie human diseases. Nature Reviews Genetics 10: 392–404.

17. Moore JH, Williams SM (2009) Epistasis and its implications for personal genetics. American Journal of Human Genetics 85: 309–320.

18. Moore JH, Asselbergs FW, Williams SM (2010) Bioinformatics challenges for genome-wide association studies. Bioinformatics 26: 445–455.

19. Wei W-H, Hemani G, Gyenesei A, Vitart V, Navarro P, et al.. (2012) Genome-wide analysis of epistasis in body mass index using multiple human populations. European Journal of Human Genetics.

20. Emily M, Mailund T, Hein J, Schauser L, Schierup MH (2009) Using biological networks to search for interacting loci in genome-wide association studies. European Journal of Human Genetics 17: 1231–1240.

21. Sun Y (2012) Integration of biological networks and pathways with genetic association studies. Human genetics: 1–10.

22. Ritchie MD (2011) Using Biological Knowledge to Uncover the Mystery in the Search for Epistasis in Genome-Wide Association Studies. Annals of Human Genetics 75: 172–182.

23. Brautbar A, Virani SS, Belmont J, Nambi V, Jones PH, et al. (2012) LPL gene variants affect apoC-III response to combination therapy of statins and fenofibric acid in a randomized clinical trial of individuals with mixed dyslipidemia. Journal of Lipid Research 53: 556–560.

24. Williams OD (1989) The Atherosclerosis Risk in Communities (ARIC) Study - Design and Objectives. American Journal of Epidemiology 129: 687–702.

25. Bild DE, Bluemke DA, Burke GL, Detrano R, Roux AVD, et al. (2002) Multi-ethnic study of atherosclerosis: Objectives and design. American Journal of Epidemiology 156: 871–881.

26. Brautbar A, Covarrubias D, Belmont J, Lara-Garduno F, Virani SS, et al. (2011) Variants in the APOA5 gene region and the response to combination therapy with statins and fenofibric acid in a randomized clinical trial of individuals with mixed dyslipidemia. Atherosclerosis 219: 737–742.

27. van Raalte DH, Li M, Pritchard PH, Wasan KM (2004) Peroxisome Proliferator-Activated Receptor (PPAR)-α: A Pharmacological Target with a Promising Future. Pharmaceutical Research 21: 1531–1538.

28. Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, et al. (2009) Common variants at 30 loci contribute to polygenic dyslipidemia. Nature genetics 41: 56–65.

29. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. The American Journal of Human Genetics 81: 559–575.

30. Mailman MD, Feolo M, Jin Y, Kimura M, Tryka K, et al. (2007) The NCBI dbGaP database of genotypes and phenotypes. Nature genetics 39: 1181–1186.

31. Romeo S, Pennacchio LA, Fu Y, Boerwinkle E, Tybjaerg-Hansen A, et al. (2007) Population-based resequencing of ANGPTL4 uncovers variations that reduce triglycerides and increase HDL. Nature genetics 39: 513–516.

32. Kathiresan S, Willer CJ, Peloso GM, Demissie S, Musunuru K, et al. (2009) Common variants at 30 loci contribute to polygenic dyslipidemia. Nat Genet 41: 56–65.

33. Köster A, Chao YB, Mosior M, Ford A, Gonzalez-DeWhitt PA, et al. (2005) Transgenic angiopoietin-like (angptl) 4 overexpression and targeted disruption of angptl4 and angptl3: regulation of triglyceride metabolism. Endocrinology 146: 4943–4950.

34. Folkersen L, van't Hooft F, Chernogubova E, Agardh HE, Hansson GK, et al. (2010) Association of genetic risk variants with expression of proximal genes identifies novel susceptibility genes for cardiovascular disease. Circ Cardiovasc Genet 3: 365–373.

35. Santamarina-Fojo S, Haudenschild C, Amar M (1998) The role of hepatic lipase in lipoprotein metabolism and atherosclerosis. Current Opinion in Lipidology 9: 211–219.

36. Toedter G, Li K, Marano C, Ma KY, Sague S, et al. (2011) Gene Expression Profiling and Response Signatures Associated With Differential Responses to Infliximab Treatment in Ulcerative Colitis. American Journal of Gastroenterology 106: 1272–1280.

37. Jones J, Otu H, Spentzos D, Kolia S, Inan M, et al. (2005) Gene signatures of progression and metastasis in renal cell cancer. Clinical Cancer Research 11: 5730–5739.

38. Raue U, Trappe TA, Estrem ST, Qian HR, Helvering LM, et al. (2012) Transcriptome signature of resistance exercise adaptations: mixed muscle and fiber type specific profiles in young and old adults. Journal of Applied Physiology 112: 1625–1636.

39. Berchtold NC, Cribbs DH, Coleman PD, Rogers J, Head E, et al. (2008) Gene expression changes in the course of normal brain aging are sexually dimorphic. Proceedings of the National Academy of Sciences of the United States of America 105: 15605–15610.

40. Noble CL, Abbas AR, Cornelius J, Lees CW, Ho GT, et al. (2008) Regional variation in gene expression in the healthy colon is dysregulated in ulcerative colitis. Gut 57: 1398–1405.

41. Sabatti C (2009) Service SK, Hartikainen AL, Pouta A, Ripatti S, et al (2009) Genome-wide association analysis of metabolic traits in a birth cohort from a founder population. Nature genetics 41: 35–46.

42. Kristiansson K, Perola M, Tikkanen E, Kettunen J, Surakka I, et al. (2012) Genome-wide screen for metabolic syndrome susceptibility Loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits. Circ Cardiovasc Genet 5: 242–249.

43. Von Eckardstein A, Nofer JR, Assmann G (2001) High density lipoproteins and arteriosclerosis role of cholesterol efflux and reverse cholesterol transport. Arteriosclerosis, Thrombosis, and Vascular Biology 21: 13–27.

44. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, et al. (2010) The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Research 38: W214–W220.

45. Barabasi AL, Jeong H, Neda Z, Ravasz E, Schubert A, et al. (2002) Evolution of the social network of scientific collaborations. Physica a-Statistical Mechanics and Its Applications 311: 590–614.

46. Rual JF, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, et al. (2005) Towards a proteome-scale map of the human protein–protein interaction network. Nature 437: 1173–1178.

47. Jeong H, Mason SP, Barabasi AL, Oltvai ZN (2001) Lethality and centrality in protein networks. Nature 411: 41–42.

48. Costanzo M, Baryshnikova A, Bellay J, Kim Y, Spear ED, et al. (2010) The Genetic Landscape of a Cell. Science 327: 425–431.

49. Lu CQ, Zhang Z, Leach L, Kearsey MJ, Luo ZW (2007) Impacts of yeast metabolic network structure on enzyme evolution. Genome Biology 8.

50. Hu T, Sinnott-Armstrong NA, Kiralis JW, Andrew AS, Karagas MR, et al.. (2011) Characterizing genetic interactions in human disease association studies using statistical epistasis networks. BMC bioinformatics 12.