| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Stroke. 2009;40:1245.)
© 2009 American Heart Association, Inc.
Original Contributions |
From the Laboratory for Genotyping Development (T.M., M.K., K.A., B.L.), Center for Genomic Medicine, The Institute of Physical and Chemical Research (RIKEN), Yokohama, Kanagawa, Japan; the Departments of Clinical Sciences (T.M., M.K., T.K., S.I., M.I.) and Environmental Medicine (M.K., K.Y., T.N., J.H., Y.D., Y.K.), Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan; and the Laboratory for Molecular Medicine (Y.N.), Human Genome Center, The Institute of Medical Science, University of Tokyo, Tokyo, Japan.
Correspondence to Michiaki Kubo, MD, PhD, Laboratory for Genotyping Development, Center for Genomic Medicine, RIKEN, 1-7-22, Suehiro-cho, Tsurumi, Yokohama, Kanagawa, 230-0045, Japan. E-mail mkubo{at}src.riken.jp
| Abstract |
|---|
|
|
|---|
Methods— We investigated the association between variations of PDE4D and ischemic stroke by 3 methods: single-marker, haplotype, and tag-single nucleotide polymorphism (SNP) analyses. In the single-marker analysis, we evaluated the association using 2 large case–control samples (1112 cases and 1112 control subjects in a sample obtained from Kyushu, Japan, and 1711 cases and 1786 control subjects in BioBank Japan) and a prospective cohort with 14 years of follow-up. These samples were analyzed both separately and pooled. Haplotype and tag-SNP analyses were performed using the 2 case–control samples together.
Results— In single-marker association tests, we found no significant association in the same direction among the 6 SNP reported in the initial study and ischemic stroke subtypes. Haplotype analysis revealed no significant association between the region around the 5'-end of the gene and combined atherothrombotic and cardioembolic infarction. Rs7730070, a SNP located around the 3'-end of PDE4D, showed the lowest nominal probability value by tag-SNP analysis but was not significant after adjustment for multiple testing (adjusted probability value =0.36).
Conclusions— These results suggest that variations in PDE4D are not associated with ischemic stroke risk in the Japanese population.
Key Words: cerebral infarct genetics PDE4D
| Introduction |
|---|
|
|
|---|
To our knowledge, 15 replication studies have been published on the association between SNPs in PDE4D and ischemic stroke.4–18 However, the results are still controversial. Of the 4 studies that examined associations between the 2 markers (AC008818–1 and SNP45) and combined ATI and CEI, none replicated the original findings.4–7 Among the 14 studies that examined at least SNP45,4–17 only one found a nominal association with combined ATI and CEI.6 Other groups reported significant associations between different phenotypes and different SNPs in the 1.5-Mb region of the gene.6–13,19 There are thought to be several reasons for these inconsistencies among the results. The sample sizes in most studies were too small and had insufficient power to detect associations.20 Sampling biases of cases and controls may have distorted true associations. Several positive findings in different SNPs might reflect associations among hidden causative variants linked to the SNPs or to the G0 haplotype. The association between variants in PDE4D and ischemic stroke risk might differ among ethnic groups.
According to the recent published criteria, replication studies should examine the same SNP or a SNP in perfect or very high linkage disequilibrium with the prior SNP on the same or a very similar phenotype. They also should show similar magnitude of effect and significance in the same direction.21 Therefore, we performed single-marker association tests between the 6 SNPs and the same subtypes of ischemic stroke as in the initial study and used a sufficient sample size. We also performed haplotype analyses in blocks B and C using tag-SNPs selected from the same regions. To examine the possibility of hidden causative SNPs, we additionally genotyped 190 tag-SNPs that covered a 2.2-Mb region, including PDE4D, and performed association analyses.
| Materials and Methods |
|---|
|
|
|---|
Another case–control sample was selected from the BioBank Japan project.24 This project was started in 2003 to collect a total of 300 000 cases who have at least one of 47 diseases by a collaborative network of 66 hospitals located throughout Japan. The registration of cases was based on diagnoses made by physicians at the affiliated hospitals. From June 2003 to March 2006, 7974 cases with ischemic stroke were registered. We selected 1711 cases diagnosed with ischemic stroke subtypes by brain imaging, the same as with the Kyushu sample. The subtypes included 1143 with lacunar infarction, 355 with ATI, and 213 with CEI. Control subjects were randomly selected from the subjects who were registered with BioBank Japan for other diseases.
For the prospective cohort study, we used a cohort population of the Hisayama study established in 1988.25 In this cohort, 2634 Hisayama residents aged
40 years and who had no history of stroke or coronary heart disease were enrolled in 1988 and continuously followed up for 14 years until the occurrence of cardiovascular disease or death. Among them, 1656 subjects participated in the examination between 2002 and 2003 and were used in the present study. During the 14-year follow-up, 67 events of first-ever ischemic stroke were observed.
Written informed consent was obtained from all study subjects. The study was approved by the ethics committees of the Graduate School of Medical Sciences at Kyushu University and the Institute of Physical and Chemical Research.
Clinical characteristics of 2 case–control samples are shown in Supplemental Table I, available online at http://stroke.ahajournals.org. In both samples, hypertension was defined as systolic blood pressure of
140 mm Hg or diastolic blood pressure of
90 mm Hg or current treatment with hypertensive medication.
|
SNP Selection and Genotyping
For the association study, we selected 6 SNPs that were significantly associated with ischemic stroke in the initial study: SNP41 (rs12153798), SNP45 (rs12188950), SNP56 (rs702553), SNP83 (rs966221), SNP87 (rs2910829), and SNP89 (rs1396476). In the haplotype analysis, we selected 16 additional tag-SNPs from the regions of blocks B and C defined in the initial study. For tag-SNP analysis, we selected 190 tag-SNPs from the 2.2-Mb region, including PDE4D. Tag-SNPs were selected from the Hapmap JPT data by the pairwise tagging method with the following criteria: r2>0.8, minor allele frequency >5%, and call rate >75%.
Genomic DNA was extracted from peripheral blood leukocytes by a standard method. We genotyped SNPs using the multiplex polymerase chain reaction-based Invader assay26 (Third Wave Technologies) or TaqMan assays (Applied Biosystems) in a blind fashion to the clinical information of study samples. All genotypes were called by visual inspection, and we determined genotype success as <10 undetermined samples in a 384-well plate. When we failed to genotype more than one 384-well plate in a total of 16 plates, we excluded the SNP from further analyses. To validate the genotyping data, we genotyped 10 SNPs in 48 subjects using direct sequencing, and the concordance rate was 99.6%.
Statistical Analysis
We examined the association both by each population and by meta-analysis. We assessed case–control association analysis and Hardy-Weinberg equilibrium by
2 test or Fisher exact test, as appropriate. In the association analysis, we mainly used an additive model and also referenced dominant and recessive models. For an easy understanding of the risk direction, we calculated the OR and 95% CI of each SNP according to the risk allele in the initial study. In a meta-analysis of the single-marker association test, pooled estimates of the ORs for 2 case–control studies and one prospective study were obtained using a fixed-effect model. Heterogeneities across the population were estimated formally using Cochrans Q test and the I2 statistic. Haplotype analysis was performed using Haploview version 4.0 (Broad Institute). For the adjustment for multiple testing, we performed a random permutation test with 10 000 replications. linkage disequilibrium was calculated as D', and haplotype blocks were defined by Gabriels criteria.27
| Results |
|---|
|
|
|---|
|
|
|
Stratified Analysis by Hypertension Status
Some replication studies showed significant associations between the SNPs in PDE4D and ischemic stroke in subjects without hypertension.11,17 Thus, we evaluated the association between the 6 SNPs and the subtypes of ischemic stroke among the combined samples stratified by hypertension status (Table 2). However, none of the SNPs were associated with ATI or the combined ATI and CEI phenotypes even in the subjects without hypertension.
|
Haplotype Analysis
Because SNP45 and SNP41, which are key SNPs for haplotype construction in block B, were monomorphic in our population, we constructed haplotypes using SNP56 and 16 additional tag SNPs selected from the regions of blocks B and C (Table 3
). In block B, none of the haplotypes were significantly associated with the combined ATI and CEI phenotypes. In block C, the most common haplotype, G-C-C-A-G, showed the lowest probability value, but the association was not significant after adjustment for multiple testing (adjusted P=0.33). There was no significant haplotype in the combined region of blocks B and C (data not shown).
|
|
Tag-SNP Analysis
To determine the possibility of a hidden causative SNP, we attempted to examine the associations between tag-SNPs in PDE4D and ischemic stroke. We selected 190 additional tag-SNPs from the 2.2-Mb region that included PDE4D and genotyped in combined samples of 2823 cases and 2898 control subjects. Because 14 SNPs did not pass our criteria, we finally analyzed 198 SNPs (the 6 reported in the initial study and 192 tag-SNPs). The genomic structure, case–control results, and linkage disequilibrium map of the 2.2-Mb region are shown in the Figure. Although the initial study showed a strong association around the region of blocks A to C, none of the SNPs in this region showed any association. The rs7730070 SNP, located around the 3'-end of PDE4D, showed the lowest probability value (OR, 1.21; 95% CI, 1.06 to 1.37; P=0.0037). However, this SNP was not linked to the 5'-end of the gene that was the causative region in the initial study (Figure, C). Moreover, this association was not significant after adjustment for multiple testing (adjusted P=0.36).
|
| Discussion |
|---|
|
|
|---|
Among the replication studies that examined variations of PDE4D and ischemic stroke, the most probable reason for the inconsistent findings is that the small sample sizes missed true associations of modest effect. Assuming our sample size, the allele frequencies of the SNPs in our control subjects, and the relative risks of the SNPs in the initial study, the power to detect associations at a significance level of 0.05 would be greater than 99% for SNP83 and SNP56, 98.3% for SNP87, and 69.7% for SNP89 in the case–control samples. In contrast, the statistical power of the prospective cohort was <30% for the 6 SNPs. However, a meta-analysis of these samples should increase the statistical power to detect the association. Therefore, if a true association exists, our study could detect the association between SNPs or haplotypes in PDE4D and ischemic stroke with high probability. A recent meta-analysis of 5216 cases and 6615 control subjects also showed that allele 0 of AC008818 and haplotype G0 carriers were associated with increased risk of ischemic stroke, but these associations become nonsignificant after exclusion of the initial study.28 These results indicate that the effect size of PDE4D variants on ischemic stroke, if it exists, may be small.
Because the initial study could not determine a causative SNP or haplotype in PDE4D, many replication studies have reported positive associations between different SNPs in PDE4D and various ischemic stroke subtypes.19 This indicates the possibility that hidden causative SNPs for ischemic stroke might exist in PDE4D. We analyzed a total of 198 tag-SNPs that covered the 2.2-Mb region, including PDE4D, but none of the SNPs were significant after adjustment for multiple testing. Because we selected tag-SNPs according to strict criteria, this analysis was able to capture the most common SNPs in PDE4D. Therefore, the previous positive findings of different SNPs may be attributable to chance.
One possible reason for the lack of association between PDE4D and ischemic stroke in our study was the difference in the ethnic background. Indeed, SNP45 and SNP41, which showed the most significant association with the combined ATI and CEI phenotypes in an Icelandic population, were monomorphic and all of the Japanese populations studied were homozygotes of the risk alleles in both SNPs. If SNP45 or SNP41 or absolutely linked variations are causative, we cannot estimate the effects of these variations on ischemic stroke, because all causative variations are homozygotes of risk alleles in both cases and control subjects.
Several limitations of this study should be discussed. First, we did not genotype the microsatellite marker, AC008818-1, in this study. However, we genotyped 16 tag-SNPs selected from the regions of blocks B and C according to strict criteria. Therefore, we believe that the effect of AC008818-1 could be sufficiently covered by haplotype analysis using tag-SNPs. Second, we could use only 1656 of 2634 subjects in the prospective cohort. Subjects who developed ischemic stroke would have a higher mortality rate than subjects who did not, and this may have resulted in the lower participation rate in this study. There is a possibility that the results of the prospective cohort might have been distorted by a survivorship bias. Third, the criteria used for classifying ischemic stroke were different between the initial study and ours. For classification of ischemic stroke, the initial study used the Trial of Org 10172 in Acute Stroke Treatment research criteria29 and we used NINDS-III.23 However, these 2 classifications are similar to each other, and we diagnosed the subtypes of ischemic stroke by adequate laboratory examinations. We believe that there is no large difference in the phenotype definition.
In conclusion, although we performed a replication study between the variations of PDE4D and ischemic stroke risk using 2 independent large case–control samples and a population-based prospective cohort, we failed to replicate the associations. We suggest that variations of PDE4D do not confer risk for ischemic stroke in the Japanese population.
| Acknowledgments |
|---|
Source of Funding
This work was supported in part by Ministry of Education, Culture, Sports, Science and Technology.
Disclosures
None.
Received June 1, 2008; revision received September 8, 2008; accepted October 22, 2008.
| References |
|---|
|
|
|---|
2. Kiely DK, Wolf PA, Cupples LA, Beiser AS, Myers RH. Familial aggregation of stroke. The Framingham Study. Stroke. 1993; 24: 1366–1371.
3. Gretarsdottir S, Thorleifsson G, Reynisdottir ST, Manolescu A, Jonsdottir S, Jonsdottir T, Gudmundsdottir T, Bjarnadottir SM, Einarsson OB, Gudjonsdottir HM, Hawkins M, Gudmundsson G, Gudmundsdottir H, Andrason H, Gudmundsdottir AS, Sigurdardottir M, Chou TT, Nahmias J, Goss S, Sveinbj
rnsdottir S, Valdimarsson EM, Jakobsson F, Agnarsson U, Gudnason V, Thorgeirsson G, Fingerle J, Gurney M, Gudbjartsson D, Frigge ML, Kong A, Stefansson K, Gulcher JR. The gene encoding phosphodiesterase 4D confers risk of ischemic stroke. Nat Genet. 2003; 35: 131–138.[CrossRef][Medline]
[Order article via Infotrieve]
4. L
hmussaar E, Gschwendtner A, Mueller JC, Org T, Wichmann E, Hamann G, Meitinger T, Dichgans M. ALOX5AP gene and the PDE4D gene in a central European population of stroke patients. Stroke. 2005; 36: 731–736.
5. Kostulas K, Gretarsdottir S, Kostulas V, Manolescu A, Helgadottir A, Thorleifsson G, Gudmundsson LJ, Thorsteinsdottir U, Gulcher JR, Stefansson K, Hillert J. PDE4D and ALOX5AP genetic variants and risk for ischemic cerebrovascular disease in Sweden. J Neurol Sci. 2007; 263: 113–117.[CrossRef][Medline] [Order article via Infotrieve]
6. Meschia JF, Brott TG, Brown RD Jr, Crook R, Worrall BB, Kissela B, Brown WM, Rich SS, Case LD, Evans EW, Hague S, Singleton A, Hardy J. Phosphodiesterase 4D and 5–lipoxygenase activating protein in ischemic stroke. Ann Neurol. 2005; 58: 351–361.[CrossRef][Medline] [Order article via Infotrieve]
7. Bevan S, Porteous L, Sitzer M, Markus HS. Phosphodiesterase 4D gene, ischemic stroke, and asymptomatic carotid atherosclerosis. Stroke. 2005; 36: 949–953.
8. Nilsson-Ardnor S, Wiklund PG, Lindgren P, Nilsson AK, Janunger T, Escher SA, Hallbeck B, Stegmayr B, Asplund K, Holmberg D. Linkage of ischemic stroke to the PDE4D region on 5q in a Swedish population. Stroke. 2005; 36: 1666–1671.
9. Nakayama T, Asai S, Sato N, Soma M. Genotype and haplotype association study of the STRK1 region on 5q12 among Japanese: a case–control study. Stroke. 2006; 37: 69–76.
10. Song Q, Cole JW, O'Connell JR, Stine OC, Gallagher M, Giles WH, Mitchell BD, Wozniak MA, Stern BJ, Sorkin JD, McArdle PF, Naj AC, Xu Q, Gibbons GH, Kittner SJ. Phosphodiesterase 4D polymorphisms and the risk of cerebral infarction in a biracial population: the Stroke Prevention in Young Women Study. Hum Mol Genet. 2006; 15: 2468–2478.
11. Brophy VH, Ro SK, Rhees BK, Lui LY, Lee JM, Umblas N, Bentley LG, Li J, Cheng S, Browner WS, Erlich HA. Association of phosphodiesterase 4D polymorphisms with ischemic stroke in a US population stratified by hypertension status. Stroke. 2006; 37: 1385–1390.
12. van Rijn MJE, Slooter AJC, Schut AFC, Isaacs A, Aulchenko YS, Snijders PJLM, Kappelle LJ, van Swieten JC, Oostra BA, van Duijn CM. Familial aggregation, the PDE4D gene, and ischemic stroke in a genetically isolated population. Neurology. 2005; 65: 1203–1209.
13. Staton JM, Sayer MS, Hankey GJ, Attia J, Thakkinstian A, Yi Q, Cole VJ, Baker R, Eikelboom JW. Association between phosphodiesterase 4D gene and ischaemic stroke. J Neurol Neurosurg Psychiatry. 2006; 77: 1067–1069.
14. Woo D, Kaushal R, Kissela B, Sekar P, Wolujewicz M, Pal P, Alwell K, Haverbusch M, Ewing I, Miller R, Kleindorfer D, Flaherty M, Chakraborty R, Deka R, Broderick J. Association of phosphodiesterase 4D with ischemic stroke: a population-based case–control study. Stroke. 2006; 37: 371–376.
15. Fidani L, Clarimon J, Goulas A, Hatzitolios AI, Evans W, Tsirogianni E, Hardy J, Kotsis A. Association of phosphodiesterase 4D gene G0 haplotype and ischaemic stroke in a Greek population. Eur J Neurol. 2007; 14: 745–749.[CrossRef][Medline] [Order article via Infotrieve]
16. Kuhlenbäumer G, Berger K, Huge A, Lange E, Kessler C, John U, Funke H, Nabavi DG, Stögbauer F, Ringelstein EB, Stoll M. Evaluation of single nucleotide polymorphisms in the phosphodiesterase 4D gene (PDE4D) and their association with ischaemic stroke in a large German cohort. J Neurol Neurosurg Psychiatry. 2006; 77: 521–524.
17. Zee RYL, Brophy VH, Cheng S, Hegener HH, Erlich HA, Ridker PM. Polymorphisms of the phosphodiesterase 4D, cAMP-specific (PDE4D) gene and risk of ischemic stroke: a prospective, nested case–control evaluation. Stroke. 2006; 37: 2012–2017.
18. Saleheen D, Bukhari S, Haider SR, Nazir A, Khanum S, Shafqat S, Anis MK, Frossard P. Association of phosphodiesterase 4D gene with ischemic stroke in a Pakistani population. Stroke. 2005; 36: 2275–2277.
19. Rosand J, Bayley N, Rost N, de Bakker PI. Many hypotheses but no replication for the association between PDE4D and stroke. Nat Genet. 2006; 38: 1091–1092.[CrossRef][Medline] [Order article via Infotrieve]
20. Gulcher JR, Kong A, Gretarsdottir S, Thorleifsson G, Stefansson K. Reply to Many hypotheses but no replication for the association between PDE4D and stroke. Nat Genet. 2006; 38: 1092–1093.[CrossRef]
21. NCI-NHGRI Working Group on Replication in Association Studies. Replicating genotype–phenotype associations. Nature. 2007; 447: 655–670.[CrossRef][Medline] [Order article via Infotrieve]
22. Kubo M, Hata J, Ninomiya T, Matsuda K, Yonemoto K, Nakano T, Matsushita T, Yamazaki K, Ohnishi Y, Saito S, Kitazono T, Ibayashi S, Sueishi K, Iida M, Nakamura Y, Kiyohara Y. A nonsynonymous SNP in PRKCH (protein kinase C
) increases the risk of cerebral infarction. Nat Genet. 2007; 39: 212–217.[CrossRef][Medline]
[Order article via Infotrieve]
23. Special report from the National Institute of Neurological Disorders and Stroke. Classification of cerebrovascular diseases III. Stroke. 1990; 21: 637–676.
24. Nakamura Y. The BioBank Japan Project. Clin Adv Hematol Oncol. 2007; 5: 696–697.[Medline] [Order article via Infotrieve]
25. Kubo M, Kiyohara Y, Kato I, Tanizaki Y, Arima H, Tanaka K, Nakamura H, Okubo K, Iida M. Trends in the incidence, mortality, and survival rate of cardiovascular disease in a Japanese community. The Hisayama Study. Stroke. 2003; 34: 2349–2354.
26. Ohnishi Y, Tanaka T, Ozaki K, Yamada R, Suzuki H, Nakamura Y. A high-throughput SNP typing system for genome-wide association studies. J Hum Genet. 2001; 46: 471–478.[CrossRef][Medline] [Order article via Infotrieve]
27. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D. The structure of haplotype blocks in the human genome. Science. 2002; 296: 2225–2229.
28. Bevan S, Dichgans M, Gschwendtner A, Kuhlenbäumer G, Ringelstein EB, Markus HS. Variation in the PDE4D gene and ischemic stroke risk. A systematic review and meta-analysis on 5200 cases and 6600 controls. Stroke. 2008; 39: 1966–1971.
29. Adams HP Jr, Bendixen BH, Kappelle LJ, Biller J, Love BB, Gordon DL, Marsh EE III. Classification of subtype of acute ischemic stroke definitions for use in a multicenter clinical trial. TOAST. Trial of Org 10172 in Acute Stroke Treatment. Stroke. 1993; 24: 35–41.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Stroke Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2009 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |