- Software article
- Open Access
- Open Peer Review
Computational genetics analysis of grey matter density in Alzheimer’s disease
BioData Miningvolume 7, Article number: 17 (2014)
Alzheimer’s disease is the most common form of progressive dementia and there is currently no known cure. The cause of onset is not fully understood but genetic factors are expected to play a significant role. We present here a bioinformatics approach to the genetic analysis of grey matter density as an endophenotype for late onset Alzheimer’s disease. Our approach combines machine learning analysis of gene-gene interactions with large-scale functional genomics data for assessing biological relationships.
We found a statistically significant synergistic interaction among two SNPs located in the intergenic region of an olfactory gene cluster. This model did not replicate in an independent dataset. However, genes in this region have high-confidence biological relationships and are consistent with previous findings implicating sensory processes in Alzheimer’s disease.
Previous genetic studies of Alzheimer’s disease have revealed only a small portion of the overall variability due to DNA sequence differences. Some of this missing heritability is likely due to complex gene-gene and gene-environment interactions. We have introduced here a novel bioinformatics analysis pipeline that embraces the complexity of the genetic architecture of Alzheimer’s disease while at the same time harnessing the power of functional genomics. These findings represent novel hypotheses about the genetic basis of this complex disease and provide open-access methods that others can use in their own studies.
Alzheimer’s disease (AD) is a progressive brain disorder and the most common form of dementia. Genetic studies have revealed a number of polymorphisms associated with risk of Alzheimer’s disease. Many of these are summarized in the online AlzGene database (http://www.alzgene.org/). However, there are many additional genetic risk factors that have not been discovered using standard association methods with Alzheimer’s disease as a discrete endpoint. One approach is to use neuroimaging methods to measure brain structure and function as endophenotypes for Alzheimer’s disease. The working hypothesis is that measures of brain structure will make it easier to identify some of the undiscovered genetic risk factors for Alzheimer’s disease. The goal of the present study was to reanalyze genome-wide association study (GWAS) data from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) using grey matter density as an endophenotype. More specifically, we present a bioinformatics approach that considers the joint effects of all polymorphisms and their aggregation in biologically-defined pathways. Moving beyond the standard one-polymorphism-at-a-time analysis paradigm will allow the formulation of new hypotheses about the genetic architecture of late-onset Alzheimer’s disease.
The data used in this study comes from the Alzheimer’s Disease Neuroimaging Initiative (ADNI), which began on October 1, 2004. The stated goal of this multisite study is to define the rate of progress of mild cognitive impairment and Alzheimer’s disease in order to create better treatments for these conditions. The study carried out functional magnetic resonance imaging (fMRI) every six to twelve months for 818 patients. A total of 733 with genetic data across three categories were studied here: 204 who are neurotypical, 354 with mild cognitive impairment, and 175 with Alzheimer’s disease. A total of 530,992 single-nucleotide polymorphisms (SNPs) were measured across the human genome and passes quality control as part of a previous genome-wide association study (GWAS). The combination of brain imaging and GWAS data makes it possible to carry out voxel-wise genome-wide association studies (vGWAS) creating a many-to-many mapping problem. In addition to the individual voxels, there are many different phenotypes that can be extracted from the brain images. Here, we analyzed grey matter density to identify new candidate genes for Alzheimer’s disease. The details of the genotypic and phenotypic data has been previously described.
A bioinformatics pipeline
The motivation for this analysis approach is to identify gene-gene interactions in Alzheimer’s disease that are not predicted by univariate effects. We combine powerful machine learning methods for detecting synergistic interactions with functional genomics data to reduce the likelihood of identifying false-positive results. Figure 1 provides an overview of our bioinformatics analysis pipeline. All methods are implemented in freely available software packages making this analysis accessible to anyone with basic bioinformatics skills. The goal of Phase I was to carry out a joint analysis of all pairs of SNPs within each gene (i.e. gene-level analysis) to allow identification of both additive effects and non-additive genetic interaction effects. The goal of Phase II of the analysis was to use a bioinformatics approach with functional genomics data (i.e. pathway-level analysis) to further address the possibility of false-positives followed by a final QMDR analysis to assess gene-gene interactions.
Phase I: SNP-SNP interaction analysis
We first mapped each of the SNPs from the GWAS to individual genes using a window of 500 kb upstream and downstream from the gene sequence. This window size was selected to capture as many regulatory SNPs as possible without assigning any one SNP to too many different nearby genes as has been previously used for these kinds of studies. Next, we carried out an exhaustive joint SNP analysis within each gene region using the Quantitative Multifactor Dimensionality Reduction (QMDR) method. MDR is a nonparametric and genetic model-free machine learning approach for detecting genetic associations that exhibit additive or non-additive effects. MDR uses a constructive induction approach to map genotypes combinations from two or more SNPs to a new single variable that makes interactions easier to detect. The QMDR extension allows modeling of quantitative traits such as grey matter density by collapsing multilocus genotypic means into those above and those below the global mean. The means of the two new groups are then compared using a two-sample t-test. In this stage of the analysis we ran QMDR only on all pairs of SNPs within a gene region. The most informative pair of SNPs for each gene was selected and the statistical significance of their corresponding QMDR models determined using a 1000-fold permutation test. The p-value of this pair of SNPs was used to assign a p-value to its corresponding gene. We used a significance level of 0.001 to select genes for the next step of the analysis. This significance level was selected to minimize type II error (false-negatives) while providing moderate control of type I error (false-positives) due to multiple testing. At this stage of the analysis we were more concerned about type II errors than type I errors. We further address the possibility of false-positives in the Phase II analysis using functional genomics data and bioinformatics analysis.
Phase II: functional genomics analysis
The list of genes selected in Phase I were used as input for the integrative multi-species prediction (IMP) webserver (http://imp.princeton.edu) that infers gene relationships using a Bayesian analysis of functional genomics data including thousands of publically available gene expression datasets. The ultimate goal of this analysis is reduce the likelihood of including false-positive genes from the statistical analysis in phase I by focusing on those genes with strong biological evidence for biological interaction. For this analysis we used a very high confidence level of 0.9 for inferring that any two genes are functionally connected. We also allowed IMP to add up to 20 additional genes in the network that were connected to our list of genes with a high confidence. This is a standard option in the software. The output of IMP is a functional gene-gene interaction network. In addition, IMP performs a gene set enrichment analysis on the genes in the network to identify those pathways with more genes than expected by chance. We used a statistical significance level of 0.05. A final list of genes appearing in the IMP network and the significantly enriched pathways were selected for final analysis with QMDR. Here, we used QMDR to model all pairwise, three-way and four-way gene-gene interactions among the SNPs in this gene list that were identified in Phase I. Finally, we assessed the nature of the gene-gene interactions (i.e. independent, redundant or synergistic) from the QMDR results by performing entropy-based analyses[10, 11] using the visualization of statistical epistasis networks (ViSEN) software package.
Results and discussion
Using QMDR, phase I of the analysis revealed a total of 20 genes with 34 unique SNPs that passed the statistical significance threshold for SNP-SNP interactions (p < 0.001). The functional relationships of these 20 genes were inferred in Phase II using a bioinformatics approach that considers the correlation of gene pairs across thousands of gene expression datasets in addition to other information such as protein-protein interactions. In addition, a gene set enrichment analysis was performed to determine whether genes appearing in the functional network occurred more frequently than expected by chance in particular biological processes as defined by Gene Ontology. This latter analysis revealed enrichment for three visual perception pathways (p < 0.05) representing three genes from the gene network (CACNA1C, FKBP4, and TRPC4) as well as two pathways for DNA repair and replication (p < 0.05) representing two genes (MCM5, MCM7). An additional five genes from the olfactory pathway were present in the gene network (OR8K1, OR8K3, OR8K5, OR5R1, and OP8U1). Thus, the Phase II bioinformatics and functional genomics analysis reduced the list of 20 genes identified in Phase I to just six. Four additional genes were added in the gene network analysis based on their functional relationships for a total of 10 genes. Only 10 SNPs were present in these genes due to some SNPs mapping to more than one nearby gene. In fact, the same two SNPs represented all five olfactory genes because they are all clustered together within the 500 kb windows that were used.
An exhaustive QMDR analysis of the 10 SNPs revealed an overall best model consisting of three SNPs (rs661090, rs12222334, and rs1570612). This model was significant based on a 1000-fold permutation test (p < 0.001). The first two SNPs are located in intergenic regions within the olfactory gene cluster while the third is located in an intron of the TRPC4 gene. It is important to note that we tried and failed to replicate this finding in an independent ADNI cohort of similar size. Lack of replication could be due to different data collection methods between the two studies or other factors such as differences in allele frequencies that are known to significantly impact the replication of gene-gene interactions.
The ViSEN analysis of these SNPs revealed that that the two olfactory SNPs had a very strong synergistic interaction in the absence of strong independent effects while the TRPC4 SNP appeared to have an effect on grey matter density that was independent of the other two (Additional file1: Figure S1). Further, there was no evidence of correlation or linkage disequilibrium in the ViSEN analysis. Interestingly, the olfactory genes all had high-confidence functional relationships with at least one other olfactory gene suggesting that the observed synergistic effect on grey matter density might have a functional genomics basis. This is an important supporting biological piece of evidence. In addition, it is known that sensory processing, especially the sense of smell, is among the first aspects to disappear at the onset of Alzheimer’s Disease[15–20]. Interestingly, copy number variants in the olfactory gene region have been previously associated with age at onset of Alzheimer’s disease. Our study is consistent with the idea that olfactory genes might play a role in the genetic architecture of Alzheimer’s disease thus making it a stronger hypothesis that needs to be further tested.
We presented here a bioinformatics pipeline for identifying gene-gene interactions in Alzheimer’s disease. As with any pipeline, a number of analysis decisions had to be made. For example, we selected a significance level of 0.001 in the phase I analysis and a confidence limit of 0.9 in the phase II analysis. These were selected to place more emphasis on using biological interactions to reduce false-positives due to multiple testing. Others might prefer to put emphasis on more stringent statistical criteria thus relying on statistical hypothesis testing for revealing true patterns. For example, a significance cutoff of 0.00001 in phase I would have eliminated the SNPs that were identified in our final best model. These are decisions that each user of the method will need to make based on their own experience and their own concerns about false-positives and power. In addition, it is important to qualify the p-value of less than 0.001 from the permutation testing for the final QMDR analysis since it is not entirely independent of the QMDR analyses performed in the first phase. The reader may want to take this into consideration when interpreting the final significance. Finally, although we modeled two-way and three-way interactions in this study it is possible that the genetic architecture of Alzheimer’s disease is even more complex with higher-order gene-gene and gene-environment interactions. As such, it is possible that our study is overly simplistic and that more advanced methods might be necessary.
The MDR software package is freely available from the authors. More information can be found at http://www.epistasis.org. The IMP software is freely available at imp.preinceton.edu.
Data used in preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at: http://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.
Bertram L, McQueen MB, Mullin K, Blacker D, Tanzi RE: Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet. 2007, 39: 17-23. 10.1038/ng1934.
Weiner MW, Veitch DP, Aisen PS, Beckett LA, Cairns NJ, Green RC, Harvey D, Jack CR, Jagust W, Liu E, Morris JC, Petersen RC, Saykin AJ, Schmidt ME, Shaw L, Siuciak JA, Soares H, Toga AW, Trojanowski JQ, Alzheimer’s Disease Neuroimaging Initiative: The Alzheimer’s Disease Neuroimaging Initiative: a review of papers published since its inception. Alzheimers Dement J Alzheimers Assoc. 2012, 8 (1 Suppl): S1-68.
Shen L, Kim S, Risacher SL, Nho K, Swaminathan S, West JD, Foroud T, Pankratz N, Moore JH, Sloan CD, Huentelman MJ, Craig DW, Dechairo BM, Potkin SG, Jack CR, Weiner MW, Saykin AJ: Whole genome association study of brain-wide imaging phenotypes for identifying quantitative trait loci in MCI and AD: A study of the ADNI cohort. NeuroImage. 2010, 53: 1051-1063. 10.1016/j.neuroimage.2010.01.042.
Stein JL, Hua X, Lee S, Ho AJ, Leow AD, Toga AW, Saykin AJ, Shen L, Foroud T, Pankratz N, Huentelman MJ, Craig DW, Gerber JD, Allen AN, Corneveaux JJ, Dechairo BM, Potkin SG, Weiner MW, Thompson P, Alzheimer’s Disease Neuroimaging Initiative: Voxelwise genome-wide association study (vGWAS). NeuroImage. 2010, 53: 1160-1174. 10.1016/j.neuroimage.2010.02.032.
Kim NC, Andrews PC, Asselbergs FW, Frost HR, Williams SM, Harris BT, Read C, Askland KD, Moore JH: Gene ontology analysis of pairwise genetic associations in two genome-wide studies of sporadic ALS. BioData Min. 2012, 5: 9-10.1186/1756-0381-5-9.
Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH: Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001, 69: 138-147. 10.1086/321276.
Moore JH, Gilbert JC, Tsai C-T, Chiang F-T, Holden T, Barney N, White BC: A flexible computational framework for detecting, characterizing, and interpreting statistical patterns of epistasis in genetic studies of human disease susceptibility. J Theor Biol. 2006, 241: 252-261. 10.1016/j.jtbi.2005.11.036.
Gui J, Moore JH, Williams SM, Andrews P, Hillege HL, van der Harst P, Navis G, Van Gilst WH, Asselbergs FW, Gilbert-Diamond D: A Simple and Computationally Efficient Approach to Multifactor Dimensionality Reduction Analysis of Gene-Gene Interactions for Quantitative Traits. PLoS One. 2013, 8: e66545-10.1371/journal.pone.0066545.
Wong AK, Park CY, Greene CS, Bongo LA, Guan Y, Troyanskaya OG: IMP: a multi-species functional genomics portal for integration, visualization and prediction of protein functions and networks. Nucleic Acids Res. 2012, 40 (Web Server issue): W484-W490.
Hu T, Chen Y, Kiralis JW, Collins RL, Wejse C, Sirugo G, Williams SM, Moore JH: An information-gain approach to detecting three-way epistatic interactions in genetic association studies. J Am Med Inform Assoc JAMIA. 2013, 20: 630-636. 10.1136/amiajnl-2012-001525.
Hu T, Sinnott-Armstrong NA, Kiralis JW, Andrew AS, Karagas MR, Moore JH: Characterizing genetic interactions in human disease association studies using statistical epistasis networks. BMC Bioinformatics. 2011, 12: 364-10.1186/1471-2105-12-364.
Hu T, Chen Y, Kiralis JW, Moore JH: ViSEN: methodology and software for visualization of statistical epistasis networks. Genet Epidemiol. 2013, 37: 283-285. 10.1002/gepi.21718.
Trojanowski JQ, Vandeerstichele H, Korecka M, Clark CM, Aisen PS, Petersen RC, Blennow K, Soares H, Simon A, Lewczuk P, Dean R, Siemers E, Potter WZ, Weiner MW, Jack CR, Jagust W, Toga AW, Lee VM-Y, Shaw LM, Alzheimer’s Disease Neuroimaging Initiative: Update on the biomarker core of the Alzheimer’s Disease Neuroimaging Initiative subjects. Alzheimers Dement J Alzheimers Assoc. 2010, 6: 230-238. 10.1016/j.jalz.2010.03.008.
Greene CS, Penrod NM, Williams SM, Moore JH: Failure to replicate a genetic association may provide important clues about genetic architecture. PLoS One. 2009, 4: e5639-10.1371/journal.pone.0005639.
Peters JM, Hummel T, Kratzsch T, Lötsch J, Skarke C, Frölich L: Olfactory function in mild cognitive impairment and Alzheimer’s disease: an investigation using psychophysical and electrophysiological techniques. Am J Psychiatry. 2003, 160: 1995-2002. 10.1176/appi.ajp.160.11.1995.
Morgan CD, Murphy C: Olfactory event-related potentials in Alzheimer’s disease. J Int Neuropsychol Soc JINS. 2002, 8: 753-763. 10.1017/S1355617702860039.
Wilson RS, Arnold SE, Schneider JA, Boyle PA, Buchman AS, Bennett DA: Olfactory impairment in presymptomatic Alzheimer’s disease. Ann N Y Acad Sci. 2009, 1170: 730-735. 10.1111/j.1749-6632.2009.04013.x.
Murphy C, Solomon ES, Haase L, Wang M, Morgan CD: Olfaction in aging and Alzheimer’s disease: event-related potentials to a cross-modal odor-recognition memory task discriminate ApoE epsilon4+ and ApoE epsilon 4- individuals. Ann N Y Acad Sci. 2009, 1170: 647-657. 10.1111/j.1749-6632.2009.04486.x.
Ansoleaga B, Garcia-Esparcia P, Llorens F, Moreno J, Aso E, Ferrer I: Dysregulation of brain olfactory and taste receptors in AD, PSP and CJD, and AD-related model. Neuroscience. 2013, 248C: 369-382.
Makizako M, Makizako H, Doi T, Uemura K, Tsutsumimoto K, Miyaguchi H, Shimada H: Olfactory identification and cognitive performance in community-dwelling older adults with mild cognitive impairment. Chem Senses. 2014, 39: 39-46. 10.1093/chemse/bjt052.
Shaw CA, Li Y, Wiszniewska J, Chasse S, Zaidi SNY, Jin W, Dawson B, Wilhelmsen K, Lupski JR, Belmont JW, Doody RS, Szigeti K: Olfactory copy number association with age at onset of Alzheimer disease. Neurology. 2011, 76: 1302-1309. 10.1212/WNL.0b013e3182166df5.
This work was funded by NIH R01 grants LM009012, LM010098 and LM011360. The computational analysis was made possible by high-performance computing infrastructure funded by NIH P20 grants GM103534 and P20 GM103506. ALZ was funded by the Women in Science Program (WISP) of Dartmouth College.
Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; BioClinica, Inc.; Biogen Idec Inc.; Bristol-Myers Squibb Company; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; GE Healthcare; Innogenetics, N.V.; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Medpace, Inc.; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Synarc Inc.; and Takeda Pharmaceutical Company. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (http://www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Disease Cooperative Study at the University of California, San Diego. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.
The authors declare they have no competing interests.
ALZ and JHM developed the method, carried out the analyses, interpreted the results and drafted the manuscript. JMF, TH and PCA assisted with method development and with the programming necessary to implement the method. CSG assisted with the implementation of the IMP method and its interpretation. LS and AJS provided the data, interpreted the results and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: A SNP-SNP interaction network derived from the ViSEN analysis. Each node or vertex in the network is a SNP with a main effects proportional to the size of the circle. Lines connecting two SNPs are proportional to the size of the synergistic interaction effects after removing the one-way effects. Triangles connecting three SNPs are proportional in size to the degree of pure three-way synergistic interaction after removing the two-way and one-way effects. Note that SNPs rs661090 and rs12222334 from our final best model have a pairwise synergistic interactions with only indirect interactions with the third SNP in the model, rs1570612. (PNG 58 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.