Pathway Analysis Using Genome-Wide Association Study Data for Coronary Restenosis – A Potential Role for the PARVB Gene

Abstract

Background

Coronary restenosis after percutaneous coronary intervention (PCI) still remains a significant limitation of the procedure. The causative mechanisms of restenosis have not yet been fully identified. The goal of the current study was to perform gene-set analysis of biological pathways related to inflammation, proliferation, vascular function and transcriptional regulation on coronary restenosis to identify novel genes and pathways related to this condition.

Methods

The GENetic DEterminants of Restenosis (GENDER) databank contains genotypic data of 556,099SNPs of 295 cases with restenosis and 571 matched controls. Fifty-four pathways, related to known restenosis-related processes, were selected. Gene-set analysis was performed using PLINK, GRASS and ALIGATOR software. Pathways with a p<0.01 were fine-mapped and significantly associated SNPs were analyzed in an independent replication cohort.

Results

Six pathways (cell-extracellular matrix (ECM) interactions pathway, IL2 signaling pathway, IL6 signaling pathway, platelet derived growth factor pathway, vitamin D receptor pathway and the mitochondria pathway) were significantly associated in one or two of the software packages. Two SNPs in the cell-ECM interactions pathway were replicated in an independent restenosis cohort. No replication was obtained for the other pathways.

Conclusion

With these results we demonstrate a potential role of the cell-ECM interactions pathway in the development of coronary restenosis. These findings contribute to the increasing knowledge of the genetic etiology of restenosis formation and could serve as a hypothesis-generating effort for further functional studies.

About the Authors:

Jeffrey J. W. Verschuren

Affiliation: Department of Cardiology, Leiden University Medical Center, Leiden, The Netherlands

Stella Trompet

Affiliations Department of Cardiology, Leiden University Medical Center, Leiden, The Netherlands, Department of Gerontology and Geriatrics, Leiden University Medical Center, Leiden, The Netherlands

M. Lourdes Sampietro

Affiliations Department Human Genetics, Leiden University Medical Center, Leiden, The Netherlands, Interuniversity Cardiology Institute of the Netherlands (ICIN), Utrecht, The Netherlands

Bastiaan T. Heijmans

Affiliations Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands, Netherlands Consortium for Healthy Ageing, Leiden University Medical Center, Leiden, The Netherlands

Werner Koch

Affiliation: Deutsches Herzzentrum München, Technische Universität München, Munich, Germany

Adnan Kastrati

Affiliation: Deutsches Herzzentrum München, Technische Universität München, Munich, Germany

Jeanine J. Houwing-Duistermaat

Affiliation: Department of Medical Statistics and Bioinformatics, Leiden University Medical Center, Leiden, The Netherlands

P. Eline Slagboom

Affiliations Molecular Epidemiology, Leiden University Medical Center, Leiden, The Netherlands, Netherlands Consortium for Healthy Ageing, Leiden University Medical Center, Leiden, The Netherlands

Paul H. A. Quax

Affiliation: Department of Vascular Surgery, Leiden University Medical Center, Leiden, The Netherlands

J. Wouter Jukema

* E-mail: [email protected]

Affiliations Department of Cardiology, Leiden University Medical Center, Leiden, The Netherlands, Netherlands Consortium for Healthy Ageing, Leiden University Medical Center, Leiden, The Netherlands, Durrer Center for Cardiogenetic Research, Amsterdam, The Netherlands, Interuniversity Cardiology Institute of the Netherlands (ICIN), Utrecht, The Netherlands

Introduction

Although interventional cardiology has advanced tremendously in recent years, restenosis after percutaneous coronary intervention (PCI) still remains a significant limitation of the procedure. [1] Restenosis is a complex disease for which the causative mechanisms have not yet been fully identified. Genetic susceptibility is known to play a role in the development of this complication. [1] Many studies have focused to identify these genetic markers associated with restenosis. Candidate gene approaches investigating restenosis have resulted in identification of genetic variation in genes involved in the inflammatory response [2], [3] and genes related to vascular remodeling [4], [5] and proliferation [6]. In 2012, our group analyzed the joint effect of the genetic variation of 36 previously described candidate genes of restenosis, demonstrating that the joint effect of all the single nucleotide polymorphisms (SNPs) in these genes together was indeed significantly associated with restenosis development. [7] This association was mainly driven by SNPs in 6 individual genes, together spanning a wide range of different functions.

The first genome-wide association study (GWAS) on restenosis only identified one intergenic region on chromosome 12. [8] Although this finding, after further functional confirmation, potentially helps our understanding of the molecular background of restenosis, it also complicates things since none of the previously described markers were identified at a GWAS significant level. The main explanation for this is likely the multifactorial nature of the disease and thus the small individual effect size of the SNPs. It has been described that instead of examining individual SNPs, analyzing the joint effect of multiple SNPs is a better method to find genetic variation in complex diseases.[9]–[11] Incorporating biological knowledge of the interplay of several genes within functional pathways may improve the power to detect associated genes not previously identified for their involvement. Thus far, the reported genetic variants only explain a small part of the expected genetic component of restenosis development. [7] This so-called “missing heritability” is not uncommon in complex diseases, and pathway approaches are expected to fill at least part of this gap. [12], [13] During the last years, several studies have been reported that successfully applied this approach, for instance, identifying biological pathways related to human longevity [14], pancreatic cancer [15], crohn’s disease [16] and psychiatric disorders. [17], [18].

The goal of the current study was to examine the joint effect of genetic variation in biological pathways related to inflammation, proliferation, vascular function and transcriptional regulation on coronary restenosis in the GENetic DEterminants of Restenosis (GENDER) study population to identify novel genes and pathways related to this condition.

Methods

GENDER Study Population

The main results of the GENDER study [19] and the design and results of the genome-wide association study (GWAS) [8], [20], have been described previously. In brief, GENDER included 3,104 consecutive unrelated symptomatic patients treated successfully by PCI for angina. During a follow-up period of 9 months, the endpoint clinical restenosis, defined as renewed symptoms requiring target vessel revascularization (TVR) either by repeated PCI or CABG, by death from cardiac causes or myocardial infarction not attributable to another coronary event than the target vessel, was recorded. During follow-up, 346 patients developed clinical restenosis. Blood samples were collected at the index procedure for DNA isolation. The study protocol conforms to the Declaration of Helsinki and was approved by the ethics committees of each participating institution (the ethics committee of the Leiden University Medical Center, the ethics committee of the Academic Medical Center Amsterdam, the ethics committee of the University Medical Center Groningen and the ethics committee of the Maastricht University Medical Center). Written informed consent was obtained from each participant before the PCI procedure.

The GWAS was performed in a subpopulation of GENDER of 325 cases of restenosis and 630 controls matched by gender, age, and several other known risk factors for restenosis like total occlusion, diabetes, current smoking and residual stenosis. Genotyping was performed using the Illumina Human 610-Quad Beadchips following the manufacturer’s instructions. After genotyping, samples and genetic markers were subjected to a stringent quality control protocol. [20] The population substructure was analyzed by means of multidimensional scaling. The vast majority of the individuals fell in the same cluster along with the European population but 67 individuals were outside this cluster and were considered as genetic outliers and therefore excluded from further analyses. The final dataset consisted of 866 individuals (295 cases, 571 controls) and 556,099 SNPs that passed all quality control criteria, together covering 89% of the common genetic variation in the European population. [20].

Pathway Selection

Pathways related to known restenosis-related processes (inflammation, vascular function, proliferation and transcription) were selected from the MSigDB database [21], which included curated gene sets from KEGG [22], BioCarta [23] and Reactome [24], [25]. A total of 54 pathways were selected (Table S1).

Replication Cohort

The replication cohort consisted of a German restenosis cohort, i.e. patients presenting with ischemic symptoms or evidence of myocardial ischemia in the presence of ≥50% de novo stenosis located in native coronary vessels. All patients were treated with PCI and stent implantation at Deutsches Herzzentrum München or 1. Medizinische Klinik rechts der Isar der Technischen Universität München. The study protocol conformed to the Declaration of Helsinki and was approved by the institutional ethics committee responsible for both participating centers, Deutsches Herzzentrum München and 1. Medizinische Klinik, Klinikum rechts der Isar, Technische Universität München. All subjects gave their written, informed consent for participation in the study. The main characteristics of this cohort and the protocols of stent placement and post-stenting therapy have been described previously [26], [27]. Briefly, this population included 1537 patients treated with implantation of bare-metal stents (BMS) and 1920 patients treated with implantation of drug-eluting stents (DES). Target lesion revascularization (TLR) within 1 year after PCI was considered as the primary endpoint. Re-hospitalization for repeat angiography was scheduled between 6 and 8 months or earlier if non-invasive evaluation or clinical presentation suggested the presence of ischemia. For the current study we only analyzed the patients receiving BMS, since this was the only type of stent implanted during the GENDER study.

The samples during the replication stage were genotyped by means of iPLEX assays [28]. The SNPs were selected for replication after step-down fine-mapping of the significant associated pathways. Finally only SNPs within these pathways showing a p-value <0.10 and linkage disequilibrium (LD) <0.5 in the GENDER study were selected. Two assays were designed using MassArray design software (Sequenom, San Diego, CA). Genotyping was performed using iPLEX assays with use of the MassARRAY methodology (Sequenom), with alleles discriminated by mass spectrometry, following manufacturer’s instructions. One SNP, rs5764560, did not fit the assay. All the genotyped SNPs had a call rate >95%. One SNP, rs7292425, was out of Hardy Weinberg equilibrium (p-value<0.001), and was excluded for further analysis. Duplicate samples (∼2.5%) showed identical genotypes. Finally, 39 SNPs passed all quality criteria and were further analyzed.

Statistical Analysis

Analyses were performed using PLINK [29], GRASS [30], and ALIGATOR [31] software. We will briefly describe all three methods. During the set-based test of PLINK the joint effect of all genetic variation, fulfilling the test constraints, within the set of genes of pathway of interest is evaluated. First a single SNP analysis of all SNPs within the pathway set is performed. Subsequently, a mean SNP statistic is calculated from the single SNP statistics of a maximum amount of independent SNPs with a p-value <0.2. SNPs are considered independent when the LD expressed in R2 is <0.5. Of the SNPs that are in LD with R2>0.5, the SNP with the lowest p-value in the single SNP analysis is selected. This analysis is repeated 10,000 times in simulated datasets with permutation of the phenotype. An empirical p-value for the SNP set is computed by calculating the number of times the test statistic of the simulated SNP sets exceeds that of the original SNP set.

GRASS calculates “eigenSNPs” for each gene in the pathway set by summarizing the variation of a gene using principal component analysis. Subsequently, one or more of these “eigenSNPs” per gene are selected using regularized logistic regression to calculate a test statistic for each pathway set. This analysis is repeated 10,000 times in simulated datasets with permutation of the phenotype. The p-value per pathway SNP set is calculated by comparing the test statistic of the original pathway SNP set with that of the combined simulated pathway SNP sets.

ALIGATOR (Association LIst Go AnnoTatOR) analyses gene sets for genes enriched with significant SNPs. Enrichment is defined as a gene set containing a larger number of significant genes than expected by chance. Replicate gene lists of the same length as the original are generated by randomly sampling SNPs (thus correcting for variable gene size). The lists are used to obtain p-values for enrichment for each gene set (by comparing the number of significant genes observed on the actual gene list to that observed on each replicate list), to correct these for testing multiple non-independent categories, and to test whether the number of significantly enriched categories is higher than expected. ALIGATOR uses data from all the SNPs tested in a gene and corrects for the variable numbers of SNPs per gene. Each gene is counted once regardless of how many significant SNPs it contains, thus eliminating the influence of LD between SNPs within genes. We used p-value cutoff <0.005 for SNPs, 5000 replicate gene lists and 1000 permutations as parameters to run ALIGATOR. Pathways were included when 2 or more individual genes contained significantly associated SNPs. [31], [32].

Considering the exploratory nature of the analysis and the considerate overlap between the pathways, associations of pathways with restenosis were considered worthwhile exploring with P<0.01. Since it has been suggested that PLINK, GRASS, and ALIGATOR provide complementary information [14], [33], we proceeded with the secondary analysis when a pathway was associated with restenosis with p<0.01 in at least one of the analyses. The pathways meeting this criterion were explored in more detail by fine-mapping of firstly the genes within those pathways and secondly the individual SNPs. For these secondary analyses, as well as during the replications stage p<0.05 was considered significant. When applying a strict Bonferroni correction, correcting for the 54 tested pathways and three different tests, the threshold for statistical significance was set at p<0.0003 ( = 0.05/(54*3).

Results

The 295 patients with restenosis were well matched with the 571 patients without restenosis, for the known risk factors for restenosis, like age, diabetes, current smoking, stenting of the culprit lesion and previous restenosis (Table 1). Hypertension and multivessel disease were more common in the cases compared to the controls.

[Figure omitted. See PDF.]

Table 1. Baseline characteristics of the GENDER study population.

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

A total of 54 pathways were analyzed in the current study. Half of these pathways (n = 27) were related to inflammation, 13 pathways were related to proliferation, six pathways were related to vascular function and finally, eight pathways were related to transcriptional regulation. The pathways associated with restenosis with a p-value below 0.01 in either PLINK, GRASS or ALIGATOR are shown in Table 2. An overview of the results of all the pathways can be found in Table S1. When applying Bonferroni correction none of the pathways were significantly associated.

[Figure omitted. See PDF.]

Table 2. Pathways associated with restenosis with PLINK or GRASS.

https://doi.org/10.1371/journal.pone.0070676.t002

In the set-based analysis using PLINK, the IL6 signaling pathway described by BioCarta was associated most significantly with clinical restenosis with a p-value of 0.0008. Using the GRASS software, this pathway had a p-value of 0.08 and the ALIGATOR analysis resulted in p = 0.01. The strongest associated pathway using GRASS was the BioCarta mitochondria pathway with a p-value of 0.006. However, this pathway was not associated with restenosis using PLINK (p = 0.80) or ALIGATOR. The insulin signaling pathway, described by BioCarta, was the strongest associated pathway using ALIGATOR (p = 0.001). Using PLINK this pathway had a p-value of 0.02, whereas GRASS resulted in p = 0.14. Other pathways that were associated in one of the analyses, considering the p<0.01 cut-off, include the inflammatory IL2 signaling pathway (PLINK p = 0.006, GRASS p = 0.012 and ALIGATOR p = 0.08) and the IL4 signaling pathway (PLINK p = 0.02, GRASS p = 0.04 and ALIGATOR p = 0.008, the proliferation platelet derived growth factor (PDGF) pathway (PLINK p = 0.008, GRASS p = 0.03 and ALIGATOR p = 0.02) and a transcriptional regulatory pathway, the vitamin D receptor (VDR) pathway (PLINK p = 0.008, GRASS p = 0.07 and ALIGATOR p = 0.003). Finally, the cell-extracellular (ECM) matrix interactions pathway, as described by Reactome, was demonstrated to be associated in PLINK (p = 0.009), but not in GRASS (p = 0.17) or ALIGATOR (p = 0.14).

Subsequently, fine mapping analysis was performed of six pathways with a p-value <0.01 in one of the three analyses. The individual genes of each pathway were analyzed using set-based analysis of PLINK (Table S2). Five of the 22 genes of the IL6 signaling pathway were significantly associated, as shown in Table 3. These genes were the RAF1 gene with a p-value of 0.001, the JAK2 gene with p = 0.006, the STAT3 gene with p = 0.006, the IL6 gene with p = 0.04, and the GRB2 gene with a p-value of 0.04 (Table 3). The IL2 signaling pathway and the PDGF signaling pathway showed considerable overlap with the IL6 pathway. Thirteen of the 22 genes (59%) of the IL2 pathway and 14 of the 32 genes (44%) of the PDGF pathway are also involved in the IL6 pathway. In addition to RAF1, JAK2 and GRB2, mentioned above, STAT5A, STAT5B and PIK3R1 were also significantly associated, p = 0.04, p = 0.01 and p = 0.04 respectively. The association of the VDR pathway was driven by three genes, NCOR1 (p = 0.003), VDR (p = 0.01) and BAZ1B (p = 0.02). The association of the cell-ECM interactions pathway was caused by two genes. The PARVB gene was associated with restenosis with a p-value of 0.001 and the LIMS2 gene with a p-value of 0.01. Finally, DIABLO was the only associated gene of the mitochondria pathway with a p-value of 0.03.

[Figure omitted. See PDF.]

Table 3. Fine mapping of the associated pathways and replication of the top SNP.

https://doi.org/10.1371/journal.pone.0070676.t003

In Table 3 the results of the single SNP analysis of the significantly associated genes is reported showing which SNPs are responsible for the significant association of these genes. Of the in total 54 individual SNPs, 34 were associated with restenosis with a p-value below 0.05. The strongest associated SNP was located in the PARVB gene of the cell-ECM interactions pathway (rs139107, p<0.001).

Baseline characteristics of the replication cohort can be found in Table S3. In the replication stage, two SNPs in the PARVB gene, which was part of the cell-ECM interactions pathway, were the only SNPs found to be significantly associated with restenosis (Table 3). The minor allele of rs5764455 was significantly more prevalent in the control group than in the cases, 42% and 36% respectively, p = 0.007. In GENDER the minor allele frequency (MAF) of this SNP was 44% and 39% in controls and cases respectively, p = 0.03. For rs5764066 the MAF was 22% in controls and 26% in cases, p = 0.05, compared to 18% versus 22% in GENDER. In addition, two SNPs in the JAK1 gene, which was part of the IL6 pathway, the IL2 pathway and the PDGF pathway, showed a borderline significant association with target lesion revascularization in the replication cohort. No wet-lab experiments were performed to elucidate the function of the two PARVB SNPs, rs5764455 and rs5764066, but online databases were explored for available data. Both SNPs were located in introns and the UCSC genome browsers database (http://genome.ucsc.edu) did not indicate that they were located in a DNAse I hypersensitive site or in other important regulatory regions. Moreover, neither of two SNPs were reported in two publically available eQTL databases (VarySysDB [34] and the eQTL database of the Pritchard lab [35]).

Discussion

We explored the joint effect of genetic variation in biological pathways on coronary restenosis development. When using a strict Bonferroni correction, none of the pathways reached the threshold of significance. However, we show a potential role for the cell-ECM interactions pathway as described by Reactome. This pathway, related to vascular remodeling, was associated with coronary restenosis in the GENDER study, in particular the PARVB gene and several SNPs in this gene were replicated in an independent study cohort. The SNPs of the five other pathways, associated in the GENDER study, did not replicate.

ECM remodeling and vascular smooth muscle cell migration and proliferation are thought to be the key processes in restenosis development. [1] Cell-ECM interactions, as described by the cell-ECM interactions pathway, play a critical role in regulating a variety of cellular processes in multicellular organisms including motility, shape change, survival, proliferation and differentiation. The PARVB (beta-parvin, also known as affixin, Gene ID: 29780) in particular, encodes a member of the parvin family of actin-binding proteins, which play a role in cytoskeleton organization and cell adhesion. Unlike other family members like alpha-parvin and gamma-parvin, beta-parvin has been shown to be expressed the highest in heart and skeletal muscle tissue. [36] PARVB is known to interact with integrin-linked kinase (ILK) and this complex is involved in the early stage of cell-substrate interaction through integrins. [37] Moreover, it has been suggested that it is a crucial modulator of cell survival [38] and that it has an important role in the angiogenic process [39], making it a likely candidate for restenosis development.

Since the introduction of GWAS many new insights on a wide array of disease have been reported. [40] However, extending these findings to biological mechanisms that explain disease development is a major ongoing challenge. The most likely explanation for this problem is that GWAS only focuses on single SNP associations. The reported loci often are located in intergenic regions or gene deserts and often are of unknown function. Furthermore, these loci have exhibited only modest effect sizes. Since biological pathways are the natural causative mechanisms resulting in disease and it is conceivable that genes functioning within these pathways interact within biological networks, it is only logical that the next step in the post-GWAS era is the exploration of these pathways. The rationale of set-based analysis is to overcome the marginally weak effect of single SNPs by analyzing a set of SNPs, since these SNPs could jointly have strong genetic effects. Examination of genetic variation of complete biological pathways, like in this set-based analysis, is therefore more likely to detect an associated mechanism by combining multiple SNPs with a possible marginal individual effect. [41], [42] In this study, the SNP with the strongest association with restenosis in the GENDER population was located in the PARVB gene of the cell-ECM interactions pathway (rs139107, p<0.001). As a reference, none of the SNPs currently examined, were located in the top 100 strongest associated SNPs of the complete GWAS analysis.(8).

For the current study we performed the analyses with PLINK, GRASS and ALIGATOR software. To date, many different statistical programs have been described to perform this type of analysis with. All methods have advantages and disadvantages. [9] In the case of PLINK and GRASS, both methods use the raw genotypes as input data. An advantage of this method is that they can take into account the LD structure within genes during their computations, which is not the case with methods performing their analyses using p-values as input data. Furthermore they both are self-contained tests that deal with the complexity of the SNP sets by permutation of the case-control status. However, differences between PLINK and GRASS include that PLINK analyzes SNP sets as a whole, and does not take the SNP distribution over the genes within the pathway into account, whereas GRASS first calculates “eigenSNPs” for each gene in the pathway set by summarizing the variation of a gene using principal component analysis. Compared to PLINK and GRASS, ALIGATOR is not a self-contained test but a competitive test, using SNP p-values instead of raw-genotypes, aiming to test whether genes in one pathway are enriched with a greater number of associated SNPs. Several comparative analyses have been performed, not resulting in a consensus of the best method, but rather suggesting to use two or more algorithms resulting in complementary findings. [14], [33].

For the current study we analyzed the data using a p-value threshold for including individual SNPs in the PLINK analyses of p<0.2 for their association with restenosis in the primary single SNP analysis. A standard setting in PLINK is often p<0.05, on relatively arbitrary grounds. Since we aimed at analyzing the combined effect of SNPs with small individual effect sizes we decided to increase this threshold to p<0.2. By increasing this threshold, more SNPs will be selected for further analyses ensuring that also the SNPs that were not associated with restenosis on their own, but that potentially did contribute to the overall effect of a pathway were taken into account. Before the start of this study we decided to apply a threshold of p<0.01 for the association of a pathway with restenosis, to identify pathways of interest. Although this threshold could be considered arbitrary, we decided that is was appropriate considering the considerate overlap between the pathways and the exploratory nature of this study. When using a strict Bonferroni correction for the 54 analyzed pathways and 3 applied tests (0.05/(54*3) = p<0.0003), none of the pathways would have qualified for further investigation. This could indicate that either the population is too small to detect the subtle genetic influence on this complex phenotype or that no solid genetic association is present. However, since to our knowledge, other (larger) restenosis study populations with GWAS data do not exist, we are limited in exploring this hypothesis regarding the sample size issue. A more strict multiple testing correction than currently used increases the risk of false negative results. Future functional studies are needed to further elucidate the biological background of restenosis development. The first step of the functional analyses would likely be to see whether the SNPs affect expression or function of the gene they are located in. If this is not the case, they could influence on other (nearby) genes, which could for instance be explored by mRNA expression analyses or allelic expression imbalance. [43] Furthermore, sequencing of the genes of interest could result in identifying functional variants in LD with the SNPs identified in the current study.

Limitations

Proper replication of these analyses requires other datasets with GWAS data. To our knowledge, the GWAS performed in the GENDER study population is the first and unfortunately is currently the only GWAS on coronary restenosis. The desired replication cohorts are therefore unavailable. Since we did want to analyze our findings in an independent study cohort we decided to genotype only the top SNPs of the associated pathways instead. Considering the minimal replication of the analyzed SNPs, the robustness of the reported associations is limited. However, the presented data do provide evidence that ECM interactions and inflammation contribute to coronary restenosis development. And the results of the current study can serve as a hypothesis-generating effort for further functional studies.

All patients in the GENDER study and in the replication cohort were either treated with balloon-angioplasty (33% in GENDER) or BMS implantation. It is unknown whether these findings are generalizable to patients treated with DES or endothelial-progenitor-cell covered stents.

Conclusion

In conclusion, although these results do not comply to the required strict multiple testing correction, we do demonstrate a potential role of the cell-ECM interactions pathway in the development of coronary restenosis. After performing additional studies, these findings could contribute to the increasing knowledge of the genetic etiology of restenosis formation. Unraveling the mechanisms leading to restenosis development remains challenging. There is a need for larger study populations specifically designed for genetic association analyses. In the select population that will still develop coronary restenosis despite advances in therapy, it will remain an imported cause of morbidity and therefore further understanding of the condition will be essential if we want to achieve personalized treatment and eradicate this costly and clinically very relevant limitation of PCI.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Pathway analysis of pathways potentially related to restenosis.

https://doi.org/10.1371/journal.pone.0070676.s001

(XLS)

Table S2.

SNP set analysis of all the separate genes of the top pathways.

https://doi.org/10.1371/journal.pone.0070676.s002

(XLS)

Table S3.

Baseline characteristics of German BMS cohort (N = 1537).

https://doi.org/10.1371/journal.pone.0070676.s003

(XLS)

Acknowledgments

We want to thank Q. Helmer for her technical assistance in performing the statistical analyses.

Author Contributions

Conceived and designed the experiments: JJWV ST BTH JJH-D PES PHAQ JWJ. Performed the experiments: JJWV ST MLS. Analyzed the data: JJWV ST JJH-D. Contributed reagents/materials/analysis tools: MLS WK AK JJH-D JWJ. Wrote the paper: JJWV ST MLS BTH WK AK JJH-D PES PHAQ JWJ.

Citation: Verschuren JJW, Trompet S, Sampietro ML, Heijmans BT, Koch W, Kastrati A, et al. (2013) Pathway Analysis Using Genome-Wide Association Study Data for Coronary Restenosis – A Potential Role for the PARVB Gene. PLoS ONE 8(8): e70676. https://doi.org/10.1371/journal.pone.0070676

References

1. Jukema JW, Verschuren JJ, Ahmed TA, Quax PH (2012) Restenosis after PCI. Part 1: pathophysiology and risk factors. Nat Rev Cardiol 9: 53–62.

2. Monraats PS, Pires NM, Agema WR, Zwinderman AH, Schepers A, et al. (2005) Genetic inflammatory factors predict restenosis after percutaneous coronary interventions. Circulation 112: 2417–2425.

3. Monraats PS, Kurreeman FA, Pons D, Sewgobind VD, de Vries FR, et al. (2007) Interleukin 10: a new risk marker for the development of restenosis after percutaneous coronary intervention. Genes Immun 8: 44–50.

4. Oguri M, Kato K, Hibino T, Yokoi K, Segawa T, et al. (2007) Genetic risk for restenosis after coronary stenting. Atherosclerosis 194: e172–e178.

5. Pons D, Monraats PS, Zwinderman AH, de Maat MP, Doevendans PA, et al. (2009) Metabolic background determines the importance of NOS3 polymorphisms in restenosis after percutaneous coronary intervention: A study in patients with and without the metabolic syndrome. Dis Markers 26: 75–83.

6. van Tiel CM, Bonta PI, Rittersma SZ, Beijk MA, Bradley EJ, et al. (2009) p27kip1–838C>A single nucleotide polymorphism is associated with restenosis risk after coronary stenting and modulates p27kip1 promoter activity. Circulation 120: 669–676.

7. Verschuren JJ, Trompet S, Postmus I, Sampietro ML, Heijmans BT, et al. (2012) Systematic testing of literature reported genetic variation associated with coronary restenosis: results of the GENDER Study. PLoS ONE 7: e42401.

8. Sampietro ML, Trompet S, Verschuren JJ, Talens RP, Deelen J, et al. (2011) A genome-wide association study identifies a region at chromosome 12 as a potential susceptibility locus for restenosis after percutaneous coronary intervention. Hum Mol Genet 20: 4748–4757.

9. Wang K, Li M, Hakonarson H (2010) Analysing biological pathways in genome-wide association studies. Nat Rev Genet 11: 843–854.

10. Sun YV (2012) Integration of biological networks and pathways with genetic association studies. Hum Genet 131: 1677–1686.

11. McClellan J, King MC (2010) Genetic heterogeneity in human disease. Cell 141: 210–217.

12. Zaitlen N, Kraft P (2012) Heritability in the genome-wide association era. Hum Genet 131: 1655–1664.

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

14. Deelen J, Uh HW, Monajemi R, van HD, Thijssen PE et al.. (2011) Gene set analysis of GWAS data for human longevity highlights the relevance of the insulin/IGF-1 signaling and telomere maintenance pathways. Age (Dordr ) In press. 10.1007/s11357-011-9340-3 [doi].

15. Wei P, Tang H, Li D (2012) Insights into pancreatic cancer etiology from pathway analysis of genome-wide association study data. PLoS ONE 7: e46887.

16. Ballard D, Abraham C, Cho J, Zhao H (2010) Pathway analysis comparison using Crohn’s disease genome wide association studies. BMC Med Genomics 3: 25.

17. Jia P, Kao CF, Kuo PH, Zhao Z (2011) A comprehensive network and pathway analysis of candidate genes in major depressive disorder. BMC Syst Biol 5 Suppl 3S12.

18. Jia P, Wang L, Fanous AH, Chen X, Kendler KS, et al. (2012) A bias-reducing pathway enrichment analysis of genome-wide association data confirmed association of the MHC region with schizophrenia. J Med Genet 49: 96–103.

19. Agema WR, Monraats PS, Zwinderman AH, de Winter RJ, Tio RA, et al. (2004) Current PTCA practice and clinical outcomes in The Netherlands: the real world in the pre-drug-eluting stent era. Eur Heart J 25: 1163–1170.

20. Sampietro ML, Pons D, de KP, Slagboom PE, Zwinderman A, et al. (2009) A genome wide association analysis in the GENDER study. Neth Heart J 17: 262–264.

21. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550.

22. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M (2012) KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res 40: D109–D114.

23. BioCarta website. Available: http://www.biocarta.com/. Accessed 1-4-22012.

24. D’Eustachio P (2011) Reactome knowledgebase of human biological pathways and processes. Methods Mol Biol 694: 49–61.

25. Croft D, O’Kelly G, Wu G, Haw R, Gillespie M, et al. (2011) Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res 39: D691–D697.

26. Koch W, Kastrati A, Mehilli J, Bottiger C, von BN, et al. (2000) Insertion/deletion polymorphism of the angiotensin I-converting enzyme gene is not associated with restenosis after coronary stent placement. Circulation 102: 197–202.

27. Hoppmann P, Erl A, Turk S, Tiroch K, Mehilli J, et al. (2009) No association of chromosome 9p21.3 variation with clinical and angiographic outcomes after placement of drug-eluting stents. JACC Cardiovasc Interv 2: 1149–1155.

28. Gabriel S, Ziaugra L, Tabbaa D (2009) SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Protoc Hum Genet Chapter 2: Unit.

29. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.

30. Chen LS, Hutter CM, Potter JD, Liu Y, Prentice RL, et al. (2010) Insights into colon cancer etiology via a regularized approach to gene set analysis of GWAS data. Am J Hum Genet 86: 860–871.

31. Holmans P, Green EK, Pahwa JS, Ferreira MA, Purcell SM, et al. (2009) Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. Am J Hum Genet 85: 13–24.

32. Kendler KS, Kalsi G, Holmans PA, Sanders AR, Aggen SH, et al. (2011) Genomewide association analysis of symptoms of alcohol dependence in the molecular genetics of schizophrenia (MGS2) control sample. Alcohol Clin Exp Res 35: 963–975.

33. Gui H, Li M, Sham PC, Cherny SS (2011) Comparisons of seven algorithms for pathway analysis using the WTCCC Crohn’s Disease dataset. BMC Res Notes 4: 386.

34. Shimada MK, Matsumoto R, Hayakawa Y, Sanbonmatsu R, Gough C, et al. (2009) VarySysDB: a human genetic polymorphism database based on all H-InvDB transcripts. Nucleic Acids Res 37: D810–D815.

35. eQTL database of the Pritchard lab. Available: http://eqtl.uchicago.edu. Accessed 29-12-22012.

36. Yamaji S, Suzuki A, Sugiyama Y, Koide Y, Yoshida M, et al. (2001) A novel integrin-linked kinase-binding protein, affixin, is involved in the early stage of cell-substrate interaction. J Cell Biol 153: 1251–1264.

37. Rosenberger G, Jantke I, Gal A, Kutsche K (2003) Interaction of alphaPIX (ARHGEF6) with beta-parvin (PARVB) suggests an involvement of alphaPIX in integrin-mediated signaling. Hum Mol Genet 12: 155–167.

38. Kimura M, Murakami T, Kizaka-Kondoh S, Itoh M, Yamamoto K, et al. (2010) Functional molecular imaging of ILK-mediated Akt/PKB signaling cascades and the associated role of beta-parvin. J Cell Sci 123: 747–755.

39. Javerzat S, Franco M, Herbert J, Platonova N, Peille AL, et al. (2009) Correlating global gene regulation to angiogenesis in the developing chick extra-embryonic vascular system. PLoS ONE 4: e7856.

40. 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.

41. Fridley BL, Biernacka JM (2011) Gene set analysis of SNP data: benefits, challenges, and future directions. Eur J Hum Genet 19: 837–843.

42. Torkamani A, Topol EJ, Schork NJ (2008) Pathway analysis of seven common diseases assessed by genome-wide association. Genomics 92: 265–272.

43. Pastinen T, Sladek R, Gurd S, Sammak A, Ge B, et al. (2004) A survey of genetic and epigenetic variation affecting human gene expression. Physiol Genomics 16: 184–193.