Skip main navigation

Genetic Heritability of Ischemic Stroke and the Contribution of Previously Reported Candidate Gene and Genomewide Associations

Originally publishedhttps://doi.org/10.1161/STROKEAHA.112.665760Stroke. 2012;43:3161–3167

Abstract

Background and Purpose—

The contribution of genetics to stroke risk, and whether this differs for different stroke subtypes, remainsuncertain. Genomewide complex trait analysis allows heritability to be assessed from genomewide association study (GWAS) data. Previous candidate gene studies have identified many associations with stoke but whether these are important requires replication in large independent data sets. GWAS data sets provide a powerful resource to perform replication studies.

Methods—

We applied genomewide complex trait analysis to a GWAS data set of 3752 ischemic strokes and 5972 controls and determined heritability for all ischemic stroke and the most common subtypes: large-vessel disease, small-vessel disease, and cardioembolic stroke. By systematic review we identified previous candidate gene and GWAS associations with stroke and previous GWAS associations with related cardiovascular phenotypes (myocardial infarction, atrial fibrillation, and carotid intima-media thickness). Fifty associations were identified.

Results—

For all ischemic stroke, heritability was 37.9%. Heritability varied markedly by stroke subtype being 40.3% for large-vessel disease and 32.6% for cardioembolic but lower for small-vessel disease (16.1%). No previously reported candidate gene was significant after rigorous correction for multiple testing. In contrast, 3 loci from related cardiovascular GWAS studies were significant: PHACTR1 in large-vessel disease (P=2.63e−6), PITX2 in cardioembolic stroke (P=4.78e−8), and ZFHX3 in cardioembolic stroke (P=5.50e−7).

Conclusions—

There is substantial heritability for ischemic stroke, but this varies for different stroke subtypes. Previous candidate gene associations contribute little to this heritability, but GWAS studies in related cardiovascular phenotypes are identifying robust associations. The heritability data, and data from GWAS, suggest detecting additional associations will depend on careful stroke subtyping.

Introduction

Stroke is a common cause of death and the major cause of adult chronic disability. Conventional risk factors such as hypertension fail to explain all stroke risk. The most robust data for a genetic component to disease, defined as the proportion of the variation in phenotype explained by genotype, so-called genetic heritability, comes from twin studies. Such studies suggest stroke is more common in monozygotic than dizygotic twins, but the number of stroke cases in these -prospective studies is small and therefore the estimates of heritability have wide CIs.1,2 Much more data are available from family history studies, which show family history of stroke is a risk factor for stroke, supporting a role for genetic risk factors.3 However, such an association could also be accounted for by shared early life environmental risk factors.

A further complicating factor is that stroke is a clinical syndrome resulting from a number of different disease processes. Approximately 80% of stroke is ischemic, whereas 20% is due to primary hemorrhage. Ischemic stroke itself has a number of subtypes with the most common being large-vessel atherosclerotic stroke (LVD), small-vessel disease (SVD), and cardioembolism (CE). It has been suggested that genetic predisposition may differ for these subtypes,4 and of note, most monogenic forms of stroke predispose to individual stroke subtypes.4 However, whether this is also the case for the majority of common ischemic stroke has been debated. This has direct relevance in deciding whether to look for genetic variants in larger, but perhaps less well-phenotyped samples or smaller but better phenotyped cohorts.

Although twin studies can obtain estimates of heritability, and differentiate environmental and genetic risk, -performing such studies in stroke is challenging. In twin cohorts the number of stroke cases is small. Furthermore, in studies performed to date, subtyping is not available, and therefore no estimates of heritability for specific stroke subtypes are -available.1,2 Recently, a method has been developed for determining the genetic component of disease heritability in complex diseases such as stroke using genotype data derived from genomewide association studies (GWAS). Genomewide complex trait analysis (GCTA)5 allows estimation of heritability by assessing the proportion of variation in case–control status explained by genotyped single nucleotide polymorphisms (SNPs). This tool estimates the variance explained by all the SNPs entered against a phenotypic trait rather than individual SNPs. An advantage of this technique is that it can be applied to large case–control data well-phenotyped data sets allowing estimates of heritability to be obtained for specific stroke subtypes.

To date 2 main approaches have been used in studies trying to identify the underlying molecular genetic basis of multifactorial ischemic stoke: the candidate gene and GWAS methods. Until recently, almost all studies have been candidate gene association studies in which the frequency of a variant in a known gene is compared between stroke cases and controls.6 A large number of such studies have been performed, but these have produced inconclusive results. The reasons for this have been explored in detail and include small sample sizes, failure to replicate, and in some cases lack of stroke subtyping.7,8 Meta-analyses of candidate gene associations have suggested a number of these associations represent real findings.9 However, these conclusions may be confounded by publication bias,10 and whether these associations really do contribute to stroke risk requires confirmation in large prospective studies with statistical correction for multiple comparisons.

More recently the GWAS approach has transformed the genetics of complex diseases.11 In GWAS, a large number of variants, which can exceed a million, and are randomly distributed throughout the genome, are genotyped. This approach does not depend on a prior hypothesis and therefore novel associations can be detected. Combined with large sample sizes, primary replication of findings, and rigorous statistical approaches, the GWAS approach has identified novel associations in a large number of complex diseases. Using this approach, novel associations have been identified with cardiovascular phenotypes related to stroke, namely coronary artery disease and atrial fibrillation, a few of which have been tested for and replicated in LVD12 and CE13,14 stroke, respectively. In addition, the GWAS approach has identified other novel associations with related cardiovascular phenotypes, which have not been tested for replication in stroke.

In this study we used a large case–control data set of >3500 cases and 5700 controls from the recent Wellcome Trust Case Control Consortium 2 (WTCCC2) ischemic stroke GWAS15 to estimate the heritability of ischemic stroke and its subtypes using the GCTA method. We then used the same data set to determine how many of the previously reported candidate gene association with ischemic stroke replicated in this data set. Lastly we determined whether other recently reported GWAS associations with phenotypes related to stroke such as coronary artery disease and atrial fibrillation were also risk factors for all ischemic stroke or for specific stroke subtypes.

Our aim was to investigate whether there were heritability differences between ischemic stroke subtypes and whether previous candidate gene studies were informative in understanding the genetic architecture of ischemic stroke, particularly when considering ischemic stroke subtypes individually.

Methods

Patient Population

We used data from the WTCCC2 ischemic stroke GWAS.15 A total of 2343 cases of imaging confirmed ischemic stroke from the United Kingdom and 1197 cases from Germany were compared with 5397 UK and 824 German controls. All cases and control subjects were white of European ancestry. All cases were extensively phenotyped with brain imaging in 100%, imaging of extracerebral vessels in 100%, electrocardiography in 100%, and echocardiography in 43%. Stroke subtyping was performed using the Trial of ORG 10172 in Acute Stroke Treatment system.16 Cases with other specified causes, including monogenic causes of stroke, were excluded. Mean (SD) age was 69.9 (13.4) and 69.0 (14.1) years in the UK and German data sets, respectively. The number of cases with the major subtypes LVD, CE, and SVD were 494, 451, and 467 in the UK cohort and 352, 334, and 110 in the German cohorts, respectively. Full details of case and control populations, genotyping, quality control, and analysis have been published15 and are summarized in the online-only Data Supplement methods. Cases were genotyped using an Illumina 660W array and control subjects using an Illumina 1 million array, which covered those SNPs included on the 660W array.

Heritability Estimates Using GCTA

WTCCC2 raw genotypes were called using GENCALL and data then went through quality control as previously described.15 Further stringent quality control filters were then applied to this data set as follows: SNPs with a minor allele frequency <1% were excluded as were cases and controls showing deviation from Hardy-Weinberg equilibrium of P=0.0001. Quality control filters of 1% individual missingness, 0.01% genotyping missingness, and estimated relatedness (using GCTA) > 0.125 were also applied. The centers were combined using PLINK, leaving a common subset of 307 110 SNPs passing quality control in both data sets, which were subsequently used for GCTA analysis. Principal components were calculated using EIGENSTRAT on an linkage disequilibrium-pruned subset of the data set. The first 20 principal components were then used as covariates in the heritability analysis to reduce genomic inflation (λ) as a consequence of population stratification. Although genomic inflation was minimal when 2 principal components were included (λ=1.06), failure to account for population stratification adequately can lead to overinflated heritability estimates using the GCTA tool. Therefore, 20 principal components were included to rule out any overestimation of heritability due to population structure.

GCTA analysis was performed on the phenotypes of all ischemic stroke, LVD, SVD, and CE to obtain heritability estimates. We also permuted random assignment of case and control status 100 times for each analysis to obtain the null distribution between cases and controls as a comparison.

Identification of Previous Associations With Stroke and Related Cardiovascular Phenotypes

Published Association With Ischemic Stroke

For selection of candidate gene associations with ischemic stroke we undertook a systematic review of the published literature using PubMed on August 16, 2011, using the search terms “Stroke” and “Genetics.” The search was limited to adult humans and studies reported in the English language. Abstracts were reviewed by 2 independent reviewers who identified associations meeting inclusion criteria. Inclusion criteria were any of the following: (1) a positive meta-analysis of >1000 cases for a single gene; (2) a reported genomewide significant association (threshold of 1e−8); (3) a reported association with subsequent replication, either in the same or multiple papers for the same gene (identified by name) or genetic variant (identified by rs identifier or amino acid change), in which both studies were >400 cases or where the initial study was >1000 cases and the replication was of any size. Mutations associated with monogenic stroke such as NOTCH3 causing cerebral autosomal-dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) were excluded. Studies were included regardless of ethnicity and stroke subtype. Variants initially in GWAS studies of related cardiovascular phenotypes but which were subsequently found to be associated with ischemic stroke (eg, PITX2, ZFHX3, 9p21 locus) were considered as GWAS hits from related cardiovascular phenotypes.

Published Associations With Related Cardiovascular Phenotypes

We included associations detected with atrial fibrillation, coronary artery disease, and myocardial infarction. In addition, we included studies of the intermediate stroke phenotype carotid intima-media thickness. We identified previously reported loci in these phenotypes identified by GWAS at a statistical threshold of 1e−8 as reported on a database of GWAS associations (at www.genome.gov/gwastudies) in August 2011. The search terms identified 3222 publications. From these we identified 32 genes from previous candidate gene association studies in ischemic stroke, which met our inclusion criteria, and a further 22 loci from related cardiovascular phenotype GWAS studies. Four loci were common to both cohorts, resulting in 50 loci in total. These are shown in Table 1. The full list of genes examined, specific SNPs, and window boundaries are shown in online-only Data Supplement Table I. All loci were examined for association in all ischemic stroke, LVD, SVD, and CE stroke.

Table 1. List of the 50 Genes and/or Loci Meeting the Inclusion Criteria.

Published Candidate Gene and GWAS Associations With Ischemic StrokePublished GWAS Associations With Related Cardiovascular Phenotypes
PhenotypeGeneChromosomeSNP IDPhenotypeGeneChromosomeSNP ID
Stroke…candidate geneACE17Atrial fibrillationKCNN31rs13376333
ALOX5AP13rs17216473PITX24rs2200733
Angiopoietin-18rs2507800ZFHX316rs7193343
APOA (LPA)6
APOE19Myocardial infarction5rs11748327
CRP1rs2794521BTN2A16rs6929846
CYP4AII1CDKN2A/B9rs10757278
CYP4F219rs2108622CELSR1/PSRC1/SORT11rs646776
CYPIIB28rs1799998CXCL1210rs1746048
DDAH11SLC5A3/MRPS6/KCNE221rs9982601
eNOS (NOS3)7rs1799983LDLR19rs1122608
Factor V Leiden1MIA31rs17465637
Fibrinogen4PCSK91rs11206510
GP1BA17PHACTR16rs12526453
GPIIIa17WDR122rs6725887
IL-67
LTC4S25rs730012Coronary artery diseaseCXCL1210rs501120
MTHFR1rs2274976HNFIA12rs2259816
NPY7rs16147KIAA146210rs3739998
PAI-1 (Serpine)7MRAS3rs9818870
Paroxonase-17rs662MTHFDIL6rs6922269
PDE4D5rs12188950SLC22A3/LPAL2/LPA6rs2048327
Prothrombin11
SGK16rs1057293Carotid intima-media thicknessGRM87rs17691394
TNF-Alpha10RYR315rs2229116
VKORC116rs9923231
Stroke…GWAS-relatedCDKN2A/B9rs4977574
NINJ212rs11833579
PCSK91rs11206510
PITX24rs1906591
PRKCH14rs2230500
ZFHX316rs7193343

rs ID denotes specific SNP for the respective locus if reported in the literature. References supporting each gene or rs ID's inclusion are listed in online-only Data Supplement Table I and in the supplementary reference bibliography.

GWAS indicates genomewide association study; SNP, single nucleotide polymorphism.

Replication of Identified Loci in the WTCCC2 Data Set

All analyses were conducted on meta-analyzed UK and German data. UK and German results were combined in a meta-analysis using METAL under a fixed effects inverse variance model weighted on the inverse of the square of the SE to produce a final probability value per SNP tested. Separate analyses were performed for all ischemic stroke and the subtypes LVD, CE, and SVD.

Where a specific rs# identifier or unambiguous position was identified in the literature, we established a 100-kb window central to that location and examined all directly genotyped variants within that window for significance. Where >1 variant had been identified as significant, we took the most significantly associated variant as the central point for establishment of the window. Additionally, where a gene name has been associated with the variant reported as significant, we established a 100-kb window around the gene boundaries as listed within HAPMAP release 28 and examined all directly genotyped variants within these boundaries. We report the single most significant variant from our GWAS data from either window in this study.

The only exception is for the PITX2 locus associated with atrial fibrillation. For this locus the most significantly associated variants lie 200 kb distal to the gene PITX2 and as such a window of 200 kb was used to capture both the gene and associated loci together with the intervening DNA sequence. Details of all regions and loci windows are contained in the online-only Data Supplement.

Correction for Multiple Testing

We initially assessed significance using a Bonferroni correction for the total number of genes investigated (n=50, P=0.001). This does not take into account the number of variants examined at each locus and therefore we also applied the Nyholt method17 with the modification proposed by Li.18 This method takes account of the linkage disequilibrium structure within regions to determine the number of tests that need to be corrected for. Applying this method to each locus independently, we summed the Nyholt values to establish a global threshold for significance. In total we examined 2872 variants, which using the Nyholt correction method was calculated to be the equivalent of 1496 independent tests, giving a threshold of 3.34e−5.

Analysis of QQ Plots

To gain an estimate of how much of the genetic contribution to stroke was accounted for by both the stroke candidate gene loci and the related cardiovascular GWAS associations, we plotted QQ plots for all variants from within 100-kb gene and SNP windows around significant candidate gene and GWAS SNP loci. These were plotted using R (http://cran.ma.imperial.ac.uk/) and the package ggplot2 and the qqman script available from http://dL.dropbox.com/u/66281/0_Permanent/qqman.r.

Results

Heritability Estimates for Ischemic Stroke and Its Subtypes

The heritability estimates are shown in Table 2. For all ischemic stroke it was estimated at 37.9% (Table 2). It was at a similar level for LVD (40.3%) and CE (32.6%) but lower for SVD (16.1%). In contrast, estimates under the null hypothesis from permuting 100 rounds of randomly assigning case–control status were 1.46% (SD 2.27%) for all ischemic stroke, 2.92% (SD 3.87%) for LVD, 2.18% (SD 3.12%) for CE, and 2.34% (SD 3.46%) for SVD.

Table 2. Estimates of Genetic Heritability for All Ischemic Stroke and Ischemic Stroke Subtypes

PhenotypeCasesControlsHeritability σ2g2p (SE)P Value*
All ischemic stroke333950260.379 (0.052)2.4×10−14
Large-vessel disease77950260.403 (0.076)8.2×10−8
Small-vessel disease55350260.161 (0.077)0.019
Cardioembolic stroke74450260.326 (0.074)2.8×10−6

*P value for likelihood ratio test against null hypothesis that σ2g=0.

Estimated proportion of variation in case–control status explained by all single nucleotide polymorphisms.

Contribution of Previously Reported Associations With Stroke and Cardiovascular Phenotypes to This Heritability

Candidate Gene Associations

Of the 32 previously reported candidate genes associated with ischemic stroke, 4 were associated using a Bonferroni correction for the number of genes examined of P=0.001: ALOX5AP with CE, APOA (LPA) with SVD, fibrinogen with all ischemic stroke, and paroxonase-1 with SVD. Applying the more conservative Nyholt correction, no association met significance. See Table 3 for full results.

Table 3. Genes Associated With Stroke in Candidate Gene Studies and Their Significance in the WTCCC2 Ischemic Stroke GWAS Data

PhenotypeGeneChromosomeAll Stroke (n=3752)LVD (n=877)SVD (n=595)CE (n=827)
Stroke candidate geneACE170.023960.010650.018850.2538
ALOX5AP130.0045560.018990.0029750.000858
Angiopoietin-180.046350.015360.0049560.02718
APOA (LPA)60.017710.0031370.0003110.004928
APOE190.028570.014730.015420.1465
CRP10.0055850.0063090.08690.2723
CYP4AII10.11370.11910.13380.3547
CYP4F2190.010270.016720.0020270.1129
CYPIIB280.2570.021910.032840.3295
DDAH110.10440.072060.035110.006028
eNOS (NOS3)70.0046350.13330.022610.01493
Factor V10.047340.015190.044870.0105
Fibrinogen40.0004360.0060970.0099020.05472
GP1BA170.10430.019720.16140.1278
GPIIIa0.023120.020280.0056880.02229
IL-670.0065750.010460.059850.03084
LTC4S50.093250.17640.19040.004574
MTHFR10.036130.002540.018090.02189
NPY70.0058210.19620.18050.004226
PAI-170.062760.029130.12990.02369
Paroxonase-170.0013040.0021180.000710.02808
PCSK910.015380.014860.11940.0363
PDE4D50.013810.0043550.003710.00317
Prothrombin110.22390.040750.066920.2787
SGK160.091840.013840.096370.02889
TNF-alpha100.037660.0060730.43040.1923
VKORC1160.51740.28110.080380.09158
Stroke…GWAS-relatedNINJ2120.018490.014210.12410.008102
PRKCH140.032330.015090.019920.001438
Coronary artery disease GWASCDKN2A/B90.029840.0025970.0013720.02849
Atrial fibrillation GWASPITX240.0021280.0048170.010064.78E-08*
ZFHX3160.0037580.023860.056495.50e-07*

SNP denotes the reported variant associated within the candidate gene if identified from the literature search. Reference(s) denoting the publication that lead to the inclusion of the variant in this study are included in online-only Data Supplement Table I and in the supplementary bibliography. Reported P values are for the most significant SNP from either the gene window or SNP window as described in the methods.

WTCCC2 indicates Wellcome Trust Case Control Consortium 2; GWAS, genomewide association study; LVD, large-vessel disease; SVD, small-vessel disease; CE, cardioembolic; SNP, single nucleotide polymorphism.

*Significance at a threshold of 3.34e−5 as calculated by the modified Nyholt method.

Significance at a nominal correction of 0.001 for the no. of genes examined in this study.

GWAS Associations of Related Cardiovascular Phenotypes

Of the 18 novel genes identified from GWAS of related cardiovascular phenotypes, 3 showed association with specific stroke subtypes at the modified Nyholt threshold: PHACTR1 in LVD (P=2.63e−6), PITX2 in CE stroke (P=4.78e−8), and ZFHX3 in CE stroke (P=5.50e−7). The population-attributable risk for the lead SNPs in HDAC9, PITX2, and ZFHX3 was 3.6%, 5.1%, and 3.3%, respectively. The previously reported association with the chromosome 9p21 locus and LVD was significant using the Bonferroni correction based on the number of genes but not significant at the Nyholt-corrected level (P=0.001372). No loci were significant in all ischemic stroke or SVD.

QQ Plot Analysis

Figure A shows the variants identified within 100 kb of candidate gene loci; Figure B shows the variants identified within 100 kb of GWAS significant associations in related cardiovascular phenotypes. The lack of deviation from the expected distribution for the ischemic stroke candidate gene loci (Figure A) suggests there is no strong effect on all ischemic stroke risk from these loci. In contrast for related cardiovascular phenotypes, the deviation from the expected distribution as shown in Figure B suggests evidence of shared genetic factors between these phenotypes and ischemic stroke.

Figure.

Figure. QQ plot of observed probability values from the Wellcome Trust Case Control Consortium 2 ischemic stroke genomewide association study data against expected probability value. A, Candidate gene loci from previously reported studies against all ischemic stroke; (B) genomewide association study significant hits from the related phenotypes of myocardial infarction, coronary artery disease, and carotid intima-media thickness against all ischemia stroke.

Discussion

Using a novel method of assessing the genetic risk of stroke, our findings suggest a moderate heritability for stroke but that the contribution varies by stroke subtype, being higher for LVD and CE stroke than for SVD. This finding is in line with previous family history studies on ischemic stroke subtypes, showing a greater risk associated with LVD and CE stroke than SVD stroke. This may in part be due to the heterogeneity of SVD as a distinct phenotypic subtype. The GCTA approach offers a novel way of estimating heritability, however, by considering all genetic information from a GWAS experiment rather than estimated heritability from familial studies. As highlighted in the introduction, previous data from twin and family history studies investigating the heritability of ischemic stroke have significant limitations, but these data from GCTA suggest genetic factors are important in stroke risk. When using this approach, we controlled for population stratification by including the first 20 principal components as covariates in our analysis, because population stratification can lead to falsely high heritability estimates.

Many studies have suggested candidate gene associations with ischemic stroke, but these are potentially subject to publication bias. In this large data set, with appropriate correction for multiple testing, we found little evidence for these candidate gene associations being important in stroke risk. Our systematic review of the literature identified 32 associations with ischemic stroke from candidate gene studies. Using a typical “candidate gene correction” based on a number of genes investigated, 4 associations were significant (ALOX5AP with CE, APOA [LPA] with SVD, fibrinogen with all ischemic stroke, and paroxonase-1 with SVD). These merit examination in even larger data sets. However, when -applying more stringent correction taking into account the number of SNPs examined, all were nonsignificant. Consistent with this, our QQ plots suggested there was no strong effect on ischemic stroke risk from the aggregated candidate gene SNPs.

These results are in contrast to those trying to replicate loci identified by GWAS studies of related cardiovascular phenotypes. Of these, 3 showed association at the more stringent Nyholt threshold; factor PHACTR1 in LVD and PITX2 and ZFHX3 in CE stroke. The 2 associations with CE stroke were first identified in atrial fibrillation14,19 and were subsequently associated with ischemic stroke,13,14 and a very limited portion of the WTCCC2 cohort used here were also used in replication of 3 SNPs in this finding. The PHACTR1 association was first reported in myocardial infarction20 and therefore might be expected to be associated with large artery stroke, although this has not previously been confirmed in stroke.

We also examined a number of other previously reported GWAS associations with related phenotypes and intermediate phenotypes for stroke. Although none was significantly associated in our population, the QQ plot for the accumulation of these variants deviated markedly from the null in all ischemic stroke, suggesting that these variants in combination do make a significant contribution to ischemic stroke risk and may be individually significant in larger sample sizes.

Three previous GWAS studies have found novel associations with ischemic stroke. We were unable to replicate the association of NINJ2 with all ischemic stroke,21 consistent with a similar failure to replicate in a larger case–control study.22 We were also unable to replicate an association with PRKCH, which has been associated with SVD in a Japanese population.23 The individual SNP implicated (rs2230500) is of very low frequency in white populations but we were unable to identify any association in the region surrounding the implicated SNP. The third association with a region on 7p21.1 encompassing the HDAC9 gene was reported for the first time in the WTCCC2 population15 and therefore we could not attempt independent replication in this data set.

Although we carried out this replication in a large cohort of >3000 cases, which is much larger than most previous stroke genetic studies, we still had only moderate power to detect associations, particularly with stroke subtypes. With 3548 cases and 5972 control subjects, we are able to detect an OR of ≥1.24. This is within the range of candidate gene associations reported in large-scale meta-analyses of up to 18 000 cases and 58 000 controls.9 However, the power is lower when considering subtypes where we are powered to detect risks in the region of ≥1.43. It is possible that a number of these variants may predispose to specific stroke subtypes when examined in larger populations. Another potential limitation of the study is that whether an association is detected depends on the method of multiple correction used, and the best approach is uncertain. For this reason we presented data both with a conventional approach controlling for the number of genes studied and a more stringent approach correcting for the number of variants studied.

It should be noted that the GCTA tool uses only GWAS data in its heritability estimates. This is not all genetic variation in the genome, but rather covers all directly genotyped SNPs and those in linkage disequilibrium with the genotyped SNPs. This will miss heritability due to rare variants not typed in the GWAS array. In addition, GCTA represents narrow-sense heritability, because it does not capture gene–environment interactions or epistatic (gene–gene interactions), denoted as broad-sense heritability. Any estimate of heritability has limitations, however, with pedigree estimates potentially biased by shared environmental effects. As a consequence, estimates based on genetic information and estimates based on familial information will differ because they are using different measures to arrive at estimates of heritability. Neither should be considered more accurate than the other.

We could find no previously published reports of heritability estimated by GTCA in ischemic stroke. However, a recent study has used GTCA to estimate narrow-sense heritability for cardiovascular risk factors using data from the Atherosclerosis Risk in Communities (ARIC). It found heritability estimates to be 34% for body mass index, 28% for waist–hip ratio, 33% for fasting glucose, 23% for fasting insulin, 47% for fasting triglyceride, 48% for fasting high-density lipoprotein, and 30% for systolic blood pressure.24

In conclusion, our results suggest most associations from previous candidate gene studies are likely to be of little or low very risk in stroke pathogenesis. In contrast, associations derived from GWAS studies of related cardiovascular phenotypes appear to be more replicable. Associations identified to date appear to be with specific stroke subtypes. Our data suggest that there is substantial heritability for ischemic stroke but that this varies markedly for different stroke subtypes. This suggests that GWAS studies are likely to identify further associations with ischemic stroke but that their success will depend on careful subtyping and adequate power within each individual stroke subtype.

Acknowledgments

We acknowledge use of the British 1958 Birth Cohort DNA collection, funded by Medical Research Council grant G0000934 and Wellcome Trust grant 068545/Z/02, and of the UK National Blood Service controls funded by the Wellcome Trust. A full list of WTCCC2 authors is provided in the supplementary material.

Footnotes

The online-only Data Supplement is available with this article at http://stroke.ahajournals.org/lookup/suppl/doi:10.1161/STROKEAHA.112.665760/-/DC1.

Correspondence to Steve Bevan, PhD, Stroke and Dementia Research Group, St Georges, University of London, Cranmer Terrace, Tooting, London, SW17 0RE, UK. E-mail

References

  • 1. Brass LM, Isaacsohn JL, Merikangas KR, Robinette CD. A study of twins and stroke.Stroke. 1992; 23:221–223.LinkGoogle Scholar
  • 2. Bak S, Gaist D, Sindrup SH, Skytthe A, Christensen K. Genetic liability in stroke: a long term follow up study of Danish twins.Stroke. 2002; 33:769–774.LinkGoogle Scholar
  • 3. Jerrard-Dunne P, Cloud G, Hassan A, Markus HS. Evaluating the genetic component of ischemic stroke subtypes. A family history study.Stroke. 2003; 34:1364–1369.LinkGoogle Scholar
  • 4. Markus HS. Unravelling the genetics of ischaemic stroke.PLoS Med.2010; 7:e1000225.CrossrefMedlineGoogle Scholar
  • 5. Yang J, Lee S, Goddard M, Visscher P. GCTA: a tool for genome wide complex trait analysis.Am J Hum Genet. 2011; 88:76–82.CrossrefMedlineGoogle Scholar
  • 6. Hirschhorn J, Lohmueller K, Byrne E, Hirschhorn K. A comprehensive review of genetic association studies.Genet Med. 2002; 4:45–61.CrossrefMedlineGoogle Scholar
  • 7. Ionnidis J, Ntzani E, Trikalinos T, Contopoulos-Ioannidis D. Replication validity of genetic association studies.Nat Genet. 2001; 29:306–309.CrossrefMedlineGoogle Scholar
  • 8. Dichgans M, Markus HS. Genetic association studies in stroke: methodological issues and proposed standard criteria.Stroke. 2005; 36:2027–2031.LinkGoogle Scholar
  • 9. Casas JP, Hingorani AD, Bautista LS, Sharma P. Meta-analysis of genetic studies in ischaemic stroke: thirty two genes involving approximately 18 000 cases and 58 000 controls.Arch Neurol. 2004; 61:1652–1661.CrossrefMedlineGoogle Scholar
  • 10. Lin B, Clyne M, Walsh M, Gomez O, Yu W, Gwinn M, et al. Tracking the epidemiology of human genes in the literature: the HuGE published literature database.Am J Epidemiol. 2006; 164:1–4.CrossrefMedlineGoogle Scholar
  • 11. McCarthy M, Abecasis G, Cardon L, Goldstein DB, Little J, Ioannidis JP, et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges.Nat Rev Genet. 2008; 9:356–369.CrossrefMedlineGoogle Scholar
  • 12. Gschwendtner A, Bevan S, Cole J, Plourde A, Matarin M, Ross-Adams H, et al. Sequence variants on chromosome 9p21.3 confer risk for atherosclerotic stroke.Ann Neurol. 2009; 65:531–539.CrossrefMedlineGoogle Scholar
  • 13. Gretarsdottir S, Thorleifsson G, Manolescu A, Styrkarsdottir U, Helgadottir A, Gschwendtner A, , et al. Risk variants for atrial fibrillation on chromosome 4q25 associate with ischaemic stroke.Ann Neurol. 2008; 64:402–409.CrossrefMedlineGoogle Scholar
  • 14. Gudbjartsson D, Holm H, Gretarsdottir S, Thorleifsson G, Walters GB, Thorgeirsson G, et al. A sequence variant in ZFHX3 on 16q22 associates with arterial fibrillation and ischaemic stroke.Nat Genet. 2009; 41:876–878.CrossrefMedlineGoogle Scholar
  • 15. Bellenguez C, Bevan S, Gschwendtnew A, Spencer CC, Burgess AI, Pirinen M, et al. Genetic analysis identifies a new susceptibility locus in HDAC9 for large vessel ischemic stroke, and supports genetic heterogeneity across stroke subtypes.Nat Genet. 2012; 44:328–333.CrossrefMedlineGoogle Scholar
  • 16. Adams HP, Bendixen BH, Kappelle LJ, Biller J, Love BB, Gordon DL, et al. Classification of subtype of acute ischemic stroke. Definitions for use in a multicentre clinical trial.TOAST. Trial of Org 10172 in acute stroke treatment.Stroke. 1993; 24:35–41.CrossrefMedlineGoogle Scholar
  • 17. Nyholt DR. A simple correction for multiple testing for single nucleotide polymorphisms in linkage disequilibrium with each other.Am J Hum Genet. 2004; 74:765–769.CrossrefMedlineGoogle Scholar
  • 18. Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix.Heredity. 2005; 95:221–227.CrossrefMedlineGoogle Scholar
  • 19. Gudbjartsson D, Arnar D, Helgadottir A, Gretarsdottir S, Holm H, Sigurdsson A, et al. Variants conferring risk of atrial fibrillation on chromosome 4q25.Nature. 2007; 448:353–357.CrossrefMedlineGoogle Scholar
  • 20. Myocardial Infarction Genetics Consortium,Kathiresan S, Voight BF, Purcell S, Musunuru K, Ardissino D, et al. Genome-wide association of early onset myocardial infarction with single nucleotide polymorphisms and copy number variants.Nat Genet.2009; 41:334–341.CrossrefMedlineGoogle Scholar
  • 21. Ikram A, Seshadri S, Bis J, Fornage M, DeStefano AL, Aulchenko YS, et al. Genomewide association studies of stroke.N Eng J Med.2009; 360:1718–1728.CrossrefMedlineGoogle Scholar
  • 22. International Stroke Genetic Consortium and Wellcome Trust Case Control Consortium 2.Failure to validate association between 12p13 variants and ischaemic stroke.N Eng J Med. 2010; 362:1547–1550.CrossrefMedlineGoogle Scholar
  • 23. Kubo M, Hata J, Ninomiya T, Matsuda K, Yonemoto K, Nakano T, et al. A nonsynonymous SNP in PRKCH (protein kinase C eta) increases the risk of cerebral infarction.Nat Genet. 2007; 39:212–217.CrossrefMedlineGoogle Scholar
  • 24. Vattikuti S, Guo J, Chow CC. Heritability and genetic correlations explained by common SNPs for metabolic syndrome traits.PLoS Genet.2012; 8:e1002637.CrossrefMedlineGoogle Scholar

eLetters(0)

eLetters should relate to an article recently published in the journal and are not a forum for providing unpublished data. Comments are reviewed for appropriate use of tone and language. Comments are not peer-reviewed. Acceptable comments are posted to the journal website only. Comments are not published in an issue and are not indexed in PubMed. Comments should be no longer than 500 words and will only be posted online. References are limited to 10. Authors of the article cited in the comment will be invited to reply, as appropriate.

Comments and feedback on AHA/ASA Scientific Statements and Guidelines should be directed to the AHA/ASA Manuscript Oversight Committee via its Correspondence page.