Volume 3 Supplement 7
Classification tree for detection of single-nucleotide polymorphism (SNP)-by-SNP interactions related to heart disease: Framingham Heart Study
© Yao et al; licensee BioMed Central Ltd. 2009
Published: 15 December 2009
The aim of this study was to detect the effect of interactions between single-nucleotide polymorphisms (SNPs) on incidence of heart diseases. For this purpose, 2912 subjects with 350,160 SNPs from the Framingham Heart Study (FHS) were analyzed. PLINK was used to control quality and to select the 10,000 most significant SNPs. A classification tree algorithm, Generalized, Unbiased, Interaction Detection and Estimation (GUIDE), was employed to build a classification tree to detect SNP-by-SNP interactions for the selected 10 k SNPs. The classes generated by GUIDE were reexamined by a generalized estimating equations (GEE) model with the empirical variance after accounting for potential familial correlation. Overall, 17 classes were generated based on the splitting criteria in GUIDE. The prevalence of coronary heart disease (CHD) in class 16 (determined by SNPs rs1894035, rs7955732, rs2212596, and rs1417507) was the lowest (0.23%). Compared to class 16, all other classes except for class 288 (prevalence of 1.2%) had a significantly greater risk when analyzed using GEE model. This suggests the interactions of SNPs on these node paths are significant.
Coronary heart disease (CHD) is a common and complex disease that is likely to involve many different genes interacting with each other and with the environment. Many studies published so far have only considered single-nucleotide polymorphisms (SNPs) in a single gene, with little consideration given to the interactions between genes. The Genetic Analysis Workshop 16 (GAW16) Framingham Heart Study (FHS) dataset provided through the database of Genotype and Phenotype (dbGaP) includes data for 550,000 SNPs and provides us with a unique opportunity to investigate this issue .
Genome-wide association studies (GWAS) systematically investigate SNPs in the entire human genome. This process allows identification of SNPs that may be associated with the disease of interest. Although GWAS are often criticized for not being hypothesis driven and can be described as data mining, they have identified unexpected and unpredictable genetic links that have advanced scientific knowledge substantially on heart disease as well as other diseases [2–6].
Challenges of statistical analysis of GWAS data have been addressed . The most difficult problems associated with GWAS analysis are: 1) how to handle extremely large data sets, often times with >10 gigabytes and 2) how to deal with a large p, small n problem due to the immense number of SNPs accompanied by a relatively small sample size. Different methods have been proposed to reduce the dimension of the data. One method is to use machine-learning approaches to select SNPs that could best explain a phenotype. A classification-tree algorithm called GUIDE , which stands for Generalized, Unbiased, Interaction Detection and Estimation, was employed in the present study. It is specifically designed to eliminate variable selection bias, a problem that can undermine the reliability of inferences from a tree structure. The algorithm of GUIDE is unbiased and is sensitive to local interactions during split selection. The FHS data from the GAW16 (Problem 2) was utilized to examine genes associated with CHD.
Data set and initial data quality checking
The FHS is a family-based study that enrolled three generations. CHD is defined as any of the following: recognized myocardial infarction diagnosed through an EKG or enzymes, coronary insufficiency, or death attributed to CHD. The third generation was excluded from the present analysis because most of them were too young to develop CHD, creating the potential for misclassification of the outcome. For the remaining 2941 study subjects, if he/she was ever diagnosed having CHD during the entire study period, this subject was classified as having CHD (case). Otherwise, the subject was classified as CHD free (control).
Dense genotyping for each study subject was performed using approximately 550,000 SNPs across 22 autosomal chromosomes (GeneChip Human Mapping 500 k Array Set and the 50 k Human Gene Focused Panel). Affymetrix conducted all genotyping for the FHS, using the 250 k Sty, 250 k Nsp, and the supplemental 50 k platforms. Quality control checks for the SNPs were performed in PLINK software  (PLINK v1.03, http://pngu.mgh.harvard.edu/purcell/plink). SNPs with >5% missing genotypes (n = 31,975) and with minor allele frequency <5% (n = 111,290) were excluded. Another 20,646 SNPs failed Hardy-Weinberg equilibrium test (p-value < 0.001). Subjects with >5% of SNPs missing (n = 29) were excluded. The remaining 2912 subjects (228 cases and 2684 controls) with 350,160 SNPs were included in subsequent analyses.
The tree algorithm GUIDE, version 7.0, was used for building a classification tree . GUIDE develops a tree by three steps: 1) a chi-square test selects the most significant split variable to split a node; 2) the split set is selected to minimize a node impurity measure (the impurity measure in GUIDE includes entropy and Gini index); 3) Steps 1 and 2 are recursively repeated until too few observations are in each node. After building a complete tree, three methods including cross-validation pruning (default), test-sample pruning, and no pruning are used to decide how much of the tree to retain. The criteria for pruning is to minimize unbiased estimate of misclassification cost. GUIDE allows fast computational speed, natural extension to data sets with categorical variables, and direct detection of local two-variable interactions. It has four useful properties: i) negligible selection bias; ii) sensitivity to curvature and local pairwise interactions between regressor variables; iii) inclusion of categorical predictor variables; and iv) choice of three roles for each ordered predictor variable: split selection only, regression modeling only, or both.
GUIDE can process a large number of SNPs in one run. However, it is still not feasible to run the entire data set with 10 GB and 350,160 SNPs at one time due to computation limitations (i.e., GUIDE stopped running, potentially due to a read buffer that is too small). To overcome this problem, the top 10,000 SNPs associated with CHD using a chi-square test implemented in PLINK were analyzed in GUIDE.
Evaluation of classes identified by GUIDE using GEE model
Current tree algorithms cannot handle dependent data such as that in the FHS, where family members are dependent. To accommodate this limitation, study subjects were treated as independent using GUIDE. The classes of SNPs identified by GUIDE were re-evaluated using a generalized estimating equations (GEE) model with the empirical variance to account for potential familial correlation.
Descriptive statistics of selected traits by generation: range, percentage distribution or mean, and standard deviation at baseline
Generation 1 (n= 356)
Generation 2 (n= 2556)
Overall (n= 2912)
t = 2.31b
Body Mass Index
t = -4.95c
t = 4.86c
t = 1.23
t = 0.68
t = -0.03
χ2 = 25.97c
χ2 = 2.24
Classification tree build by GUIDE
Functions of 16 SNPs that determined the splits in the classification tree
TBC1D14 TBC1 domain family member 14
KCTD20 potassium channel tetramerization domain containing 20
MAGI2 membrane associated guanylate kinase, WW and PDZ domain containing 2
ZNF195 zinc finger protein 195
KRT86 keratin 86
TRHDE thyrotropin-releasing hormone degrading enzyme
ERG v-ets erythroblastosis virus E26 oncogene homolog (avian)
Risk of CHD in each of 17 classes of classification tree
Prevalence of CHD in each of 17 classes identified by GUIDEa
No. CHD (%)
Prevalence of CHD in people older than 45 years was estimated to be 1.6%-16.8%, depending on age and sex according to NHANES 1999-2002 data . In the FHS, about 8% of the first and second generation participants had CHD. The prevalence of CHD in Class 16 (determined by SNPs rs1894035, rs7955732, rs2212596, and rs1417507) was the lowest (0.23%), while the prevalence of CHD in Class 19 (determined by SNPs rs1894035, rs7955732, rs2212596, and rs41009) was the highest (27%). The results suggest that individuals with certain combinations of genotypes are more resistant to developing CHD, while others are more susceptible. The SNPs involved in this risk gradient have not been reported in the CHD literature; they are novel candidates for future CHD research.
One important limitation of this study is that, like other machine-learning approaches, the computational capacity of GUIDE is limited by its inability to handle dependent data found in family-based studies, although the findings from GUIDE were replicated in a GEE model that accounts for familial correlation. Nevertheless, our analysis demonstrated that the unbiased selection tree algorithm, GUIDE, could be useful in reducing the number of SNPs in GWAS. It can be used to study gene-gene interactions associated with complex diseases, such as CHD. In addition, GUIDE can use any categorical or quantitative variables to do the classification (environmental or genetic).
List of abbreviations used
Coronary heart disease
Framingham Heart Study
Genetic Analysis Workshop 16
Genome-wide association studies
Generalized estimating equations
Generalized, Unbiased, Interaction Detection and Estimation
The Genetic Analysis Workshops are supported by NIH grant R01 GM031575 from the National Institute of General Medical Sciences.
We are grateful to Dr. Wen-Yin Loh for his help on GUIDE, other instructors (Dr. Karl Broman, Dr. Bret Payseur, and Ms. Kristin Meyers), and classmates in the "Analytic Methods in Genetic Epidemiology" course at the University of Wisconsin-Madison.
This article has been published as part of BMC Proceedings Volume 3 Supplement 7, 2009: Genetic Analysis Workshop 16. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/3?issue=S7.
- Kaiser J: Genomic: databases - NIH goes after whole genome in search of disease genes. Science. 2006, 311: 933-10.1126/science.311.5763.933a.View ArticlePubMedGoogle Scholar
- Bentham J, Bhattacharya S: Genetic mechanisms controlling cardiovascular development. Ann NY Acad Sci. 2008, 1123: 10-19. 10.1196/annals.1420.003.View ArticlePubMedGoogle Scholar
- Engert JC, Lemire M, Faith J, Brisson D, Fujiwara TM, Roslin NM, Brewer G, Montpetit A, Darmond-Zwaig C, Renaud Y, Doré C, Bailey SD, Verner A, Tremblay G, St-Pierre J, Bétard C, Platko J, Rioux JD, Morgan K, Hudson TJ, Gaudet D: Identification of a chromosome 8p locus for early-onset coronary heart disease in a French Canadian population. Eur J Hum Genet. 2008, 16: 105-114. 10.1038/sj.ejhg.5201920.View ArticlePubMedGoogle Scholar
- Larson MG, Atwood LD, Benjamin EJ, Cupples LA, D'Agostino RB, Fox CS, Govindaraju DR, Guo CY, Heard-Costa NL, Hwang SJ, Murabito JM, Newton-Cheh C, O'Donnell CJ, Seshadri S, Vasan RS, Wang TJ, Wolf PA, Levy D: Framingham Heart Study 100 K project: genome-wide associations for cardiovascular disease outcomes. BMC Med Genet. 2007, 8 (suppl 1): S5-10.1186/1471-2350-8-S1-S5.PubMed CentralView ArticlePubMedGoogle Scholar
- Larson MG, Atwood LD, Benjamin EJ, Cupples LA, D'Agostino RB, Fox CS, Govindaraju DR, Guo CY, Heard-Costa NL, Hwang SJ, Murabito JM, Newton-Cheh C, O'Donnell CJ, Seshadri S, Vasan RS, Wang TJ, Wolf PA, Levy D: A common allele on chromosome 9 associated with coronary heart disease. Science. 2007, 316: 1488-1491. 10.1126/science.1142447.View ArticleGoogle Scholar
- McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, Boerwinkle E, Hobbs HH, Cohen JC: Chromosome 9p21.3 coronary heart disease locus genotype and prospective risk of CHD in healthy middle-aged men. Clin Chem. 2008, 54: 467-474. 10.1373/clinchem.2007.095489.View ArticleGoogle Scholar
- Ziegler A, Konig IR, Thompson JR: Biostatistical aspects of genome-wide association studies. Biom J. 2008, 50: 8-28. 10.1002/bimj.200710398.View ArticlePubMedGoogle Scholar
- Loh W-Y: Regreesion trees with unbiased variable selection and interaction detection. Stat Sinica. 2002, 12: 361-386.Google Scholar
- Purcell SNB, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- American Heart Association: Prevalence of Coronary Heart Disease by Age and Sex based on National Health and Nutrition Examination Survey (NHANES 1999-2002), National Center for Health Statistics and NHLBI. [http://www.americanheart.org]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.