Identifying genetic interactions associated with late-onset Alzheimer’s disease
© Floudas et al.; licensee BioMed Central. 2014
Received: 30 May 2014
Accepted: 6 December 2014
Published: 19 December 2014
Identifying genetic interactions in data obtained from genome-wide association studies (GWASs) can help in understanding the genetic basis of complex diseases. The large number of single nucleotide polymorphisms (SNPs) in GWASs however makes the identification of genetic interactions computationally challenging. We developed the Bayesian Combinatorial Method (BCM) that can identify pairs of SNPs that in combination have high statistical association with disease.
We applied BCM to two late-onset Alzheimer’s disease (LOAD) GWAS datasets to identify SNPs that interact with known Alzheimer associated SNPs. We also compared BCM with logistic regression that is implemented in PLINK. Gene Ontology analysis of genes from the top 200 dataset SNPs for both GWAS datasets showed overrepresentation of LOAD-related terms. Four genes were common to both datasets: APOE and APOC1, which have well established associations with LOAD, and CAMK1D and FBXL13, not previously linked to LOAD but having evidence of involvement in LOAD. Supporting evidence was also found for additional genes from the top 30 dataset SNPs.
BCM performed well in identifying several SNPs having evidence of involvement in the pathogenesis of LOAD that would not have been identified by univariate analysis due to small main effect. These results provide support for applying BCM to identify potential genetic variants such as SNPs from high dimensional GWAS datasets.
KeywordsGenome-wide association study Epistasis Alzheimer’s disease Bayesian networks
Elucidating the genetic basis of common diseases will lead to an understanding of the biological mechanisms that underlie such diseases and can help in risk assessment, diagnosis, prognosis and development of new therapies. During the past several decades genetic linkage studies have been effective in mapping genetic loci responsible for many Mendelian diseases that are caused by a single genetic variant . More recently, genetic studies have indicated that most common diseases are likely to be polygenic where multiple genetic variants acting singly and in combination underlie the expression of disease .
The commonest type of genetic variation is the single nucleotide polymorphism (SNP) that results when a single nucleotide is replaced by another in the genome sequence. The development of high-throughput genotyping technologies has led to a flurry of genome-wide association studies (GWASs) with the aim of discovering SNPs that are associated with common diseases. GWASs have been moderately successful in identifying SNPs associated with common diseases and traits. However, in most cases the identified SNPs have small effect sizes, and the proportion of heritability explained is quite modest. One view is that SNPs may interact in subtle ways that lead to substantially greater effects than the effect due to any single SNP. Another view is that common diseases may be due to rare and usually deleterious SNPs that cause disease in individual patients and that in different individuals or subpopulations the disease is caused by different deleterious SNPs.
This paper addresses the challenge of identifying interacting SNPs that may have small effects and describes a Bayesian combinatorial method (BCM) for identifying such interactions that are associated with disease. This method has been shown empirically to perform well on low dimensional synthetic data . However, to our knowledge BCM has not been applied to real-world datasets with a large number of SNPs. In this paper we apply BCM to two late-onset Alzheimer’s disease GWAS datasets to identify SNPs that interact with known Alzheimer associated SNPs.
As background, we provide brief summaries about GWASs, genetic interactions, and Alzheimer’s disease in the following sections.
Genome-wide association studies
The development of high-throughput genotyping technologies that assay hundreds of thousands of SNPs or more, along with the identification of SNPs in the human genome by the International HapMap Project led to the emergence of GWASs. GWASs are typically case–control studies aimed at discovering SNPs – either as disease causing variants or as markers of disease – that are associated with a common disease or trait. The success of GWASs is based in large part on the common disease - common variant hypothesis. This hypothesis posits that common diseases in most individuals are caused by relatively common genetic variants that have low penetrance and hence have small to moderate influence in causing disease. An alternative hypothesis is the common disease - rare variant hypothesis, which posits that many rare variants underlie common diseases and each variant causes disease in relatively few individuals with high penetrance. Both these hypotheses likely contribute to common diseases with genetic variants may range from rare to the common SNPs.
GWAS data is typically analyzed for univariate associations between SNPs and the disease of interest; the statistical tests used include the Pearson’s chi-square test, the Fisher’s exact test, the Cochran-Armitage trend test, and odds ratios . SNPs identified as significant by univariate analyses may be further examined for interactions among them using methods such as logistic regression.
Genetic interactions, also known as epistasis, can be defined biologically as well as statistically. Biologically, epistasis refers to gene-gene interaction when the action of one gene is modified by one or several other genes. Statistically, epistasis refers to interaction between variants at multiple loci in which the total effect of the combination of variants at the different loci may differ considerably from a linear combination of the effects of individual loci. The detection of statistical epistasis has the potential to indicate genetic loci that have a biological interaction .
Statistical methods for identifying genetic interactions can be broadly divided into exhaustive and non-exhaustive methods. Exhaustive methods examine all possible SNP-subsets and examples include Multifactor Dimensionality Reduction (MDR) , and the BCM  that we describe in the next section. Examples of non-exhaustive methods include BOolean Operation-based Screening (BOOST), SNPHarvester, and SNPRuler. We briefly describe these methods below.
The software package PLINK that is used widely for the analysis of GWAS datasets implements logistic regression for the detection of SNP-SNP interactions in either all or specific sets of SNPs in a dataset .
MDR exhaustively evaluates all 1-, 2-, 3-, ..n-SNP subsets where n is specified by the user. It combines the variables in a SNP subset to construct a single binary variable and uses classification accuracy of the binary variable to evaluate a SNP-subset. Since MDR does not scale up beyond a few hundred SNPs, for high dimensional data a multivariate filtering algorithm called ReliefF is applied to reduce the number of SNPs to a few hundred ,,,.
BOOST uses a two-step procedure . In the screening step, it uses an approximate likelihood ratio statistic that is computationally efficient and computes it for all pairs of SNPs. Only those SNPs that pass a threshold in the first step are examined for significant interaction effect using the classical likelihood ratio test that is computationally more expensive.
SNPHarvester is a stochastic search algorithm that uses a two-step procedure to identify epistatic interactions . In the first step it identifies 40–50 significant SNP groups using a stochastic search strategy, and in the second step, it fits a penalized logistic regression model to each group.
SNPRuler searches in the space of SNP rules and uses a branch-and-bound strategy to prune the huge number of possible rules in GWAS data . An example of a rule is X 1 = 0 ^ X 2 = 2 ⇒ Z = 1 (X 1 and X 2 are SNPs, the three genotypes that a SNP can take are coded as 0, 1 and 2 and Z is a binary outcome variable). The quality of a rule is evaluated with the chi-square statistic.
Alzheimer’s disease (AD) is the commonest neurodegenerative disease associated with aging and the commonest cause of dementia . AD affects about 3% of all people between ages 65 and 74, about 19% of those between 75 and 84, and about 47% of those over 85. AD is characterized by adult onset of progressive dementia that typically begins with subtle memory failure and progresses to a slew of cognitive deficits like confusion, language disturbance and poor judgment .
AD is typically divided into early-onset Alzheimer’s disease (EOAD) in which the onset of disease is before 60 years of age and late-onset Alzheimer’s disease (LOAD) in which the onset is at or after 60 years of age. EOAD is rare and exhibits an autosomal dominant mode of inheritance. The genetic basis of EOAD is well established, and mutations in one of three genes (amyloid precursor protein gene - APP, presenilin 1 gene - PSEN1, or presenilin 2 gene - PSEN2) account for most cases of EOAD .
LOAD is widespread and is estimated to strike almost half of all people over the age of 85. LOAD is believed to be a disease with both genetic and environmental influences, and elucidating the role of genetic factors in the pathogenesis and development of LOAD has been a major focus of research for more than a decade. One genetic risk factor for LOAD that has been consistently replicated is the apolipoprotein E (APOE) locus determined by the combined genotypes at the loci rs429358 (APOE*4) and rs7412 (APOE*2) . In the past few years, GWASs have identified several additional genetic loci associated with LOAD -.
Bayesian combinatorial method
BCM uses a Bayesian network (BN) to model a set of SNPs and interactions among them and their association with disease, and the model is evaluated with a Bayesian score. It then exhaustively searches a space of all possible models to identify high scoring models.
Bayesian network model and score
In addition to computational efficiency and statistical power BCM has several additional advantages. The BCM score can be adapted to combine knowledge with data which has the potential to enhance the analysis of high dimensional GWAS datasets. Biological knowledge or results from analyses of earlier studies can be encoded in P(M) as a prior distribution over the models. The BN model used by BCM can be used for non-genetic variables such as environmental effects of disease. Thus, BCM can be used to identify combinations of genetic and environmental effects on disease. Finally, BCM can handle unbalanced datasets and its applicability is not restricted to datasets with approximately equal cases and controls.
Materials & methods
This section describes the GWAS datasets, the experimental methods, and previously identified LOAD SNPs.
We used two different LOAD GWAS datasets in our experiments. The first dataset was part of the University of Pittsburgh Alzheimer’s Disease Research Center (ADRC) that is described elsewhere . This dataset consists of 2,245 individuals, of which 1290 had LOAD and 955 did not. For each individual, the genotype data consists of 682,685 SNPs on autosomal chromosomes.
The second dataset was collected by the Translational Genomics Research Institute (TGen) . This dataset consists of 1,411 individuals, of which. 861 had LOAD and 550 did not. For each individual, the imputed genotype data consists of 234,665 SNPs on autosomal chromosomes. For each individual, the genotype data consists of 502,627 SNPs; the original investigators analyzed 312,316 SNPs after applying quality controls. We used those 312,316 SNPs, plus two additional APOE SNPs from the same study namely, rs429358 and rs7412.
The encoding of the genotypes followed the raw file format of the whole-genome analysis package PLINK, which contains a header line and then one line per individual. Each line has v + 6 fields (where v is the number of variants); the first six fields contain identification codes, the sex and the phenotype of the individual. The remaining fields contain variant genotypes, coded as a single allele dosage number (0, 1, 2 of minor allele) . The University of Pittsburgh IRB approved the use of the datasets for the study.
BCM searches exhaustively over all possible SNP-BN models in a dataset. For a GWAS dataset with half a million SNPs, the number of SNP-BN models is 2 n = 9.95 × 10150514 and the number of SNP-BN models with just 2 SNPs is = 1.25 × 1011. Thus, the search space is very large and it is computationally infeasible to evaluate every model in the space .
We addressed this challenge by applying BCM to a restricted space of SNP-BN models that consisted of a subset of all possible 2-SNP models. We considered only those 2-SNP models where one of the SNPs in a model is a member of a set of SNPs previously known to be associated with LOAD and the second SNP is any SNP (excluding the first SNP) in the dataset of interest. Since the number of known LOAD associated SNPs is much smaller than the number of SNPs in a dataset, it was computationally tractable to search this space of SNP-BN models. The selection of the previously identified LOAD SNPs that we used is described in the next section.
We applied BCM to each of the two GWAS datasets separately and analyzed in detail the top scoring 200 SNP-BN models. We chose to examine the top 200 models because the model score decreased substantially after the 200 models. From each SNP-BN model, we extracted the SNP that was not in the set of previously identified LOAD SNPs. We mapped these SNPs to genes and considered only intragenic SNPs for further analyses. We performed the SNP to gene mapping with BioQ, a web-service which uses dbSNP build 135 and Genome Assembly GRCh37.p5 . We performed enrichment analysis of the annotations of the associated genes in the Gene Ontology (GO) with the web-based tool GeneCoDis. For a set of genes GeneCoDis retrieves the associated GO terms, and identifies and ranks those GO terms that are significantly enriched in the set of genes ,. Enriched functional descriptors facilitate the interpretation of the gene set. The hierarchical nature of the GO annotations however means that the set of enriched GO terms may contain terms closely related in a parent–child relationship . Such redundant terms confound the interpretation. Therefore, we further examined the GO terms associated with the intragenic SNPs using the REViGo webserver. The REViGo software evaluates the semantic similarity between the enriched terms, identifies the most informative common ancestors and the related redundant GO terms and groups the latter under their ancestors . The resulting set facilitates simultaneous examination of the enriched GO terms at two levels: a detailed one, at the lowest level overrepresented term and a more abstract one at the highest level common ancestor of overrepresented terms. The detailed level can reveal specific genes of interest whereas the abstract level serves a compact overview of the processes, functions and cellular compartments associated with the genes in the set.
In addition to the ontology analysis of the top scoring 200 SNP-BN models, we performed additional analyses of the top scoring 30 SNP-BN models. We analyzed the genes associated with the intragenic SNPs for differential expression in AD, through the ArrayExpress web server  and biological function analysis. Differential gene expression in relation to AD aims to integrate experimental evidence from transcriptomic analysis with those of genomic analysis. Up-regulation or down-regulation in AD of a gene in our results indicates increased biological plausibility for the reported genetic interaction. Finally, elements from the functional description of a gene (expression site, function related to the nervous system or pathways of LOAD, previous literature) were considered as supporting the biological relevance of an identified interaction.
We also compared BCM with logistic regression (LR) that is implemented in PLINK, since LR is typically used in genetic epidemiology and association studies. In PLINK, we used the --epistasis option that provides an LR test for interaction that assumes an allelic model ,. We applied LR to the ADRC and the TGEN datasets to identify pairwise interactions between the dataset SNPs and the set of previously identified LOAD SNPs that was used with BCM.
Previously identified LOAD SNPs
Previously identified LOAD SNPs
Odds ratio (95% CI)
Results and discussion
This section describes the results that were obtained from applying BCM to the ADRC LOAD dataset and from applying BCM to the ADRC and the TGen GWAS datasets.
Top scoring SNP-BN models
In the ADRC dataset, the known SNP in each of the top scoring 200 SNP-BN models is rs429358 (APOE*4). The dataset SNPs from the top scoring 200 models included 92 intragenic SNPs mapping to 77 genes.
In the TGen dataset, the known SNP in SNP-BN models ranked 1 and 10 to 200 is rs429358 and the dataset SNPs included 82 intragenic SNPs that mapped to 69 genes. In the 8 SNP-BN models ranked 2–9 the known SNP was located on genes GAB2, MS4A6A, MS4A4E, CR1, PICALM, SORL1, TF and the dataset SNP for all 8 models is rs7412. SNPs rs429358 and rs7412 are located on the APOE gene and their combined genotypes determine the APOE allelic status which is known to be the strongest genetic variant that is predictive of LOAD.
A comparison of the top scoring 200 SNP-BN models for the ADRC and TGEN datasets shows that they have in common two dataset SNPs, rs7412 (APOE gene) and rs4420638 (APOC1 gene). Additionally, two of the intragenic dataset SNPs from each dataset map to genes CAMK1D (rs11257738 in ADRC and rs17151584 in TGen) and FBXL13 (rs7779121 in ADRC and rs17475512 in TGen).
GO term analysis
The most informative common ancestors of the overrepresented GO terms obtained from GeneCoDis for the ADRC dataset are given in Additional file 1: Table S3 and for the TGen dataset are given in Additional file 1: Table S4 in the Supplemental Tables. In both sets nervous system-related terms are enriched (e.g., regulation of dendrite development, nervous system development, regulation of axon extension, short term memory), as well as terms related to cholesterol and lipid metabolism (e.g., lipid metabolic process, chylomicron), beta amyloid (beta amyloid binding) cell membranes (e.g., integral to membrane, plasma membrane, postsynaptic, clathrin-coated endocytic vesicle), calmodulin and intracellular calcium homeostasis (e.g., calmodulin binding, cytosolic calcium ion transport) and the immune system (immunoglobulin binding). Overrepresentation of these terms shows that the identified genes from both datasets include genes that are members of biochemical pathways involved in LOAD pathophysiology ,.
Functional description and expression of genes associated with the top 30 dataset SNPs in the ADRC dataset
Expression in AD
1-SNP model rank
p value of pair from PLINK
Appears to modulate the interaction of APOE with beta-migrating VLDL. Binds free fatty acids.
Translocase of outer mitochondrial membrane 40 homolog
Channel-forming subunit of the translocase of the mitochondrial outer membrane (TOM) complex, essential for protein import into mitochondria.
APOE is essential for the normal catabolism of triglyceride-rich lipoprotein constituents. Known risk factor for LOAD.
Specifically expressed in the brain, highly expressed in the cortex. Organizes the subcellular localization of a variety of proteins.
Transmembrane protein 217
Expressed in the brain
Mothers against DPP homolog 6
Negative regulation of BMP and TGF-beta/activin-signaling. BMP-6 is increased in AD brains and leads to impaired neurogenesis . Reduced TGF-beta signaling is involved in neurodegeneration and promotes AD like changes in mice .
Neuronal PAS domain protein 3.
Transcription factor. May regulate genes involved in neurogenesis. Associated with schizophrenia and mental retardation
Phosphatidic acid phosphatase type 2 domain containing 1A
Neuropeptide FF-amide peptide precursor
Modulation of morphine-induced antinociception.
Solute carrier family 25
Known also as ornithine decarboxylase (ODC). Mitochondrial oxoadipate carrier, part of polyamine synthesis pathway.
Member RAS oncogene family
Intracellular protein transportation. Regulated by miRNA155, which also regulates PICALM (a known AD association).
unc-5 homolog D (C. elegans)
Netrin receptor: netrins are secreted proteins that direct axon extension and cell migration during neural development. APP also binds Netrin-1 and in transgenic mice this suppresses amyloid beta peptide production .
Chromodomain helicase DNA binding protein 9, PPARA -interacting complex 320 kDa protein
Transcriptional co-activator for PPARA. The APOE gene promoter has a binding site for PPAR alpha. Low CHD9 activity could reduce APOE levels. Increase in APOE transcription has been shown to clear amyloid beta in AD mouse models .
Contactin 4, Brain-derived immunoglobulin superfamily protein 2
Mainly expressed in brain. Neuronal membrane protein that may play a role in the formation of axon connections in the developing nervous system. Associated with Spinocerebellar Ataxia, Amyotrophic Lateral Sclerosis, 3p deletion syndrome.
Functional description and expression of genes associated with the top 30 dataset SNPs in the TGen dataset
Gene symbol (SNP)
Differential expression in AD
1-SNP model rank
p value of pair from PLINK
APOE is essential for the normal catabolism of triglyceride-rich lipoprotein constituents. Known risk factor for LOAD.
Appears to modulate the interaction of APOE with beta-migrating VLDL. Binds free fatty acids.
Chromosome 10 open reading frame 11
A brain-expressed gene. Haploinsufficiency of C10orf11 contributes to the cognitive defects in 10q22 syndrome .
von Willebrand factor C domain-containing protein 2 (Brorin, Brain-specific chordin-like protein)
Encodes a secreted bone morphogenic protein (BMP) antagonist. The encoded protein is possibly involved in neural function and development and may have a role in cell adhesion. BMP-6 is increased in AD brains and leads to impaired neurogenesis .
Pleckstrin and Sec7 domain containing 3
Guanine nucleotide exchange factor for ARF6 that contributes to the regulation of dendritic branching .
Glucoside xylosyltransferase 2
Elongates the O-linked glucose attached to EGF-like repeats in the extracellular domain of Notch proteins , which are substrates of γ-secretase, the enzyme involved in amyloid beta production .
Underexpressed in a murine AD model 
Gamma-aminobutyric acid (GABA) B receptor, 2
Target for autophagy regulation in neurodegenerative diseases .
Ectonucleotide pyrophosphatase/phosphodiesterase 2
Hydrolyzes lysophospholipids to produce lysophosphatidic acid (LPA) in extracellular fluids. Predominantly expressed in brain, placenta, ovary, and small intestine. Secreted by most body fluids including serum and cerebrospinal fluid .
Glucagon-like peptide 1 receptor
Member of the glucagon receptor family (also includes glucagon, GLP-2, secretin, GHRH and GIP receptors).In the brain located in hypothalamus and brainstem. Protective against amyloid beta accumulation in rats .
MOCO sulphurase C-terminal domain containing 1
A mitochondrial oxidoreductase, cofactor: molybdenum, is expressed in the brain. MOSC1 is a target for miR-129-5p, like GABBR2, and miR-155, like PICALM.
Transmembrane 4 L six family member 20
Tetraspannin superfamily member. Tetraspanins are often thought to act as scaffolding proteins, anchoring multiple proteins to one area of the cell membrane. Other tetraspanin superfamily members have been implicated in Notch signaling and g-secretase activity modulation .
Comparison of BCM with logistic regression applied to the ADRC and TGEN datasets
In the ADRC dataset, comparing the 200 top-ranking SNP-BN models of BCM with the top 200 SNP pairs of LR we have 5 models in common. In the case of the TGEN dataset there are 9 models in common between the 200 top-ranking SNP-BN models of BCM and the top 200 SNP pairs of LR. A comparison of the top 200 SNP pairs obtained with LR from the ADRC and TGEN datasets reveals no SNP pair in common.
Examining all pairs of SNPs in a GWAS dataset for identifying interacting SNP pairs is usually not computationally tractable due to the large number of SNP pairs. We addressed this challenge by examining only a subset of all pairs of SNPs where one member of the pair is drawn from a small set of previously known disease-associated SNPs and by using a Bayesian score to evaluate that statistical association of a SNP pair with the disease. We applied this strategy to two LOAD GWAS datasets and our results show that it can identify interacting SNPs of plausible biological significance. Moreover, this strategy finds SNPs that would be overlooked in a univariate analysis because they exhibit small main effects; however, they are detected when paired with another SNP due to interaction effects.
Evaluation of the identified SNPs did not include experimental validation in the laboratory. Instead, we gathered evidence from the literature that supports that the SNPs we identified may play a role in the biological mechanisms underlying LOAD. To do so, we selected the intragenic SNPs and identified the relevant genes. We used ontology enrichment analysis of the list of genes to provide an overview of the functions of the identified genes. We then used descriptions in gene databases to examine each gene for appropriate context (tissue expression, function, product protein interactions). Finally, we examined the literature for studies reporting association of the gene, of the gene transcript (expression studies) or the encoded protein with LOAD.
In both LOAD GWAS datasets that we examined, the previously known disease-associated SNP that was identified is either rs429358 (APOE*4) or rs7412 (APOE*2); these SNPs reside in the APOE gene which is known to be the strongest genetic determinant for LOAD. GO term enrichment analysis of the dataset SNPs identified terms that are relevant to biochemical pathways implicated in the pathogenesis of LOAD such as lipid metabolic process, calmodulin binding, nervous system development and multiple membrane-related terms.
Gene expression analysis of the dataset SNPs showed that for each dataset studied a majority of the genes corresponding to the top 30 dataset SNPs are differentially expressed in LOAD. Functional annotations and literature evidence that are presented with the expression data in the relevant tables further support the role of these genes in the pathogenesis of LOAD.
Among the genes corresponding to the top 200 dataset SNPs, besides the APOE gene, three other genes are common to both datasets: APOC1, CAMK1D and FBXL13. Evidence supporting the interaction of APOC1 with APOE is presented in the analysis for the top 30 dataset SNPs. The second gene in common is the CAMK1D (calcium/calmodulin-dependent protein kinase ID) gene that encodes a member of the Ca2+/calmodulin-dependent protein kinase 1 subfamily of serine/threonine kinases family of calmodulin kinases , which modulate neuronal development and plasticity . The gene is overexpressed in AD and is expressed in the brain especially during hippocampal formation with high expression in the pyramidal cell layers ,. CAMK1D interacts physically with CALM1 (calmodulin), which has been associated with AD risk . The CAMK1D protein may regulate calcium-mediated granulocyte function and activates MAPK3 (Mitogen-activated protein kinase 3). It phosphorylates -in vitro- the transcription factor CREM (cAMP responsive element binding) isoform beta and probably CREB1 . The CREB pathway is involved in memory formation and CREB phosphorylation has been proposed as a signalling pathway involved in the pathogenesis of AD , while CREB pathway down-regulation may have a role in exacerbations of AD . Another member of the serine/threonine kinases family of calmodulin kinases, the neuronal CaM kinase II phosphorylates tau protein on ser262, an important step in the formation of neurofibrillary tangles in AD . The third gene in common between the two datasets is the FBXL13 (F-box and leucine-rich repeat protein 13) gene that encodes a protein belonging to the F-box protein family. Members of this family have a characteristic approximately 40-amino acid F-box motif and take part in SCF (SKP1-CUL1-F-box protein) complexes that act as protein-ubiquitin ligases . The ubiquitin-proteasome system is involved in protein turnover and degradation and is perturbed in AD . An SCF complex of another F box protein (FBXW7) is probably involved in the degradation of PSEN1 protein .
In addition to the genes corresponding to the top 30 top scoring SNP-BN models in the ADRC dataset, we found other genes in lower scoring SNP-BN models with plausible associations with LOAD. In the 80th scoring model (dataset SNP rs7793977), gene PION [pigeon homolog (Drosophila)], also known as GSAP (gamma-secretase-activating protein), is known to increase amyloid beta production . In the 196th scoring model, (dataset SNP rs6534145), the PDE5A (phosphodiesterase 5A, cGMP-specific) gene could be implicated to LOAD pathogenesis via two different mechanisms. PDE5A is a substrate of CASP3 (caspase 3) , which in turn has been shown to be involved in the early synaptic dysfunction in a mouse model of AD . It has also been shown that inhibition of PDE5A results in a decrease in the transcription of Wnt/β-catenin . A reduction in Wnt signalling has been implicated in the amyloid beta-dependent neurodegeneration in LOAD .
While BCM has been applied to low dimensional synthetic data with good results , in this paper we have applied it to GWAS datasets. BCM has several advantages. It is computationally more efficient than the widely used MDR . Since BCM uses the Bayesian paradigm, the BCM score represents a coherent way to combine knowledge with data. Biological knowledge or results from analyses of earlier studies can be encoded as a prior distribution over the models that can then be used in Equation 1. Use of informative priors is becoming common in the analysis of microarray expression studies, and a similar strategy can be employed for genomic data.
In this paper we applied our method to two GWAS datasets, as replication is considered to be a necessary step in the validation of GWAS findings . However, many markers identified through GWAS have failed to replicate. Proposed explanations for this include environmental interactions, genetic heterogeneity, inadequate statistical power and population differences, i.e. inter-population differences in the minor allele frequencies (MAF). It has been shown that racial differences directly influence the odds ratios of validated markers, an explanation being the differing MAF among populations . In addition, results of simulation studies have shown that even small differences in allelic frequencies can affect the detection of main as well as of interaction effect, and may lead to non-replication or even reversal of the direction of the association . Therefore, failure to replicate a finding across two different populations should not necessarily be viewed as proof of lack of true association of the initial finding.
A limitation of our study is the use of GWAS datasets related to a single disease, although it is an important disease. In future research, we plan to apply and investigate the utility of BCM on GWAS datasets related to additional diseases. Another limitation is the use just 10 previously known LOAD-associated SNPs. In future work, we plan to explore the use of a larger set of known LOAD associated SNPs that will include SNPs with weaker evidence of being associated with LOAD. In addition, we plan to study the effect of excluding the APOE SNPs rs429358 and rs7412 which are present in every SNP pair we examined for biological plausibility. Another limitation is that we did not use informative prior probabilities for encoding prior knowledge from the literature and previous GWASs. BCM can be extended easily to allow the incorporation of informative priors and inclusion of informative priors in the analysis is an interesting area for study.
We applied BCM to two LOAD GWAS datasets to identify pairs of SNPs that in combination have high statistical association with development of LOAD. To reduce the large search space of all possible parts of SNPs in a GWAS dataset we restricted BCM to evaluate those SNP pairs where one of the SNP was drawn from a set of 10 previously known LOAD associated SNPs. Our results identified several SNPs that have biological evidence of being involved in the pathogenesis of LOAD that would not have been identified by univariate analysis alone due to small main effect but were identified in conjunction with another SNP. These results provide support for applying BCM to identify potential genetic variants such as SNPs from high dimensional GWAs datasets.
We thank Mr. Kevin Bui for his help in data preparation and in implementing the BCM algorithm in software.
CSF was supported in part by an award from the Gerondelis Foundation, SV was supported in part by NLM grant HHSN276201000030C, and M. Ilyas Kamboh was supported by National Institutes of Health grants AG030653, AG005133 and AG041718.
- Hardy J, Singleton A: Genomewide association studies and human disease. N Engl J Med. 2009, 360: 1759-1768. 10.1056/NEJMra0808700.View ArticlePubMedPubMed CentralGoogle Scholar
- Thornton-Wells TA, Moore JH, Haines JL: Genetics, statistics and human disease: analytical retooling for complexity. Trends Genet. 2004, 20: 640-647. 10.1016/j.tig.2004.09.007.View ArticlePubMedGoogle Scholar
- Visweswaran S, Wong A-KI, Barmada MM: A Bayesian method for identifying genetic interactions. AMIA Annu Symp Proc. 2009, 2009: 673-677.PubMedPubMed CentralGoogle Scholar
- Balding DJ: A tutorial on statistical methods for population association studies. Nat Rev Genet. 2006, 7: 781-791. 10.1038/nrg1916.View ArticlePubMedGoogle Scholar
- Cordell HJ: Detecting gene-gene interactions that underlie human diseases. Nat Rev Genet. 2009, 10: 392-404. 10.1038/nrg2579.View ArticlePubMedPubMed CentralGoogle Scholar
- Hahn LW, Ritchie MD, Moore JH: Multifactor dimensionality reduction software for detecting gene-gene and gene-environment interactions. Bioinformatics. 2003, 19: 376-382. 10.1093/bioinformatics/btf869.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Moore J, White B: Tuning ReliefF for genome-wide genetic analysis.Evol Comput Mach Learn Data Min Bioinformatics 2007, 4447:166–175.,Google Scholar
- Wan X, Yang C, Yang Q, Xue H, Fan X, Tang NLS, Yu W: BOOST: A fast approach to detecting gene-gene interactions in genome-wide case–control studies. Am J Hum Genet. 2010, 87: 325-340. 10.1016/j.ajhg.2010.07.021.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang C, He Z, Wan X, Yang Q, Xue H, Yu W: SNPHarvester: a filtering-based approach for detecting epistatic interactions in genome-wide association studies. Bioinformatics. 2009, 25: 504-511. 10.1093/bioinformatics/btn652.View ArticlePubMedGoogle Scholar
- Wan X, Yang C, Yang Q, Xue H, Tang NLS, Yu W: Predictive rule inference for epistatic interaction detection in genome-wide association studies. Bioinformatics. 2010, 26: 30-37. 10.1093/bioinformatics/btp622.View ArticlePubMedGoogle Scholar
- Goedert M, Spillantini MG: A century of Alzheimer’s disease. Science. 2006, 314: 777-781. 10.1126/science.1132814.View ArticlePubMedGoogle Scholar
- Bertram L, Lill CM, Tanzi RE: The genetics of Alzheimer disease: back to the future. Neuron. 2010, 68: 270-281. 10.1016/j.neuron.2010.10.013.View ArticlePubMedGoogle Scholar
- Avramopoulos D: Genetics of Alzheimer’s disease: recent advances. Genome Med. 2009, 1: 34-10.1186/gm34.View ArticlePubMedPubMed CentralGoogle Scholar
- Holtzman DM, Morris JC, Goate AM: Alzheimer’s disease: the challenge of the second century. Sci Transl Med. 2011, 3: 77sr1-PubMedPubMed CentralGoogle Scholar
- Wijsman EM, Pankratz ND, Choi Y, Rothstein JH, Faber KM, Cheng R, Lee JH, Bird TD, Bennett DA, Diaz-Arrastia R, Goate AM, Farlow M, Ghetti B, Sweet RA, Foroud TM, Mayeux R: Genome-wide association of familial late-onset Alzheimer’s disease replicates BIN1 and CLU and nominates CUGBP2 in interaction with APOE. PLoS Genet. 2011, 7: e1001308-10.1371/journal.pgen.1001308.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu X, Pickering E, Liu YC, Hall S, Fournier H, Katz E, Dechairo B, John S, Van Eerdewegh P, Soares H: Meta-analysis for genome-wide association study identifies multiple variants at the BIN1 locus associated with late-onset Alzheimer’s disease. PLoS One. 2011, 6: e16616-10.1371/journal.pone.0016616.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamboh MI, Demirci FY, Wang X, Minster RL, Carrasquillo MM, Pankratz VS, Younkin SG, Saykin AJ, Jun G, Baldwin C, Logue MW, Buros J, Farrer L, Pericak-Vance MA, Haines JL, Sweet RA, Ganguli M, Feingold E, DeKosky ST, Lopez OL, Barmada MM: Genome-wide association study of Alzheimer’s disease. Transl Psychiatry. 2012, 2: e117-10.1038/tp.2012.45.View ArticlePubMedPubMed CentralGoogle Scholar
- Hollingworth P, Harold D, Sims R, Gerrish A, Lambert J-C, Carrasquillo MM, Abraham R, Hamshere ML, Pahwa JS, Moskvina V, Dowzell K, Jones N, Stretton A, Thomas C, Richards A, Ivanov D, Widdowson C, Chapman J, Lovestone S, Powell J, Proitsi P, Lupton MK, Brayne C, Rubinsztein DC, Gill M, Lawlor B, Lynch A, Brown KS, Passmore PA, Craig D: Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer’s disease. Nat Genet. 2011, 43: 429-435. 10.1038/ng.803.View ArticlePubMedPubMed CentralGoogle Scholar
- Reiman EM, Webster JA, Myers AJ, Hardy J, Dunckley T, Zismann VL, Joshipura KD, Pearson JV, Hu-Lince D, Huentelman MJ, Craig DW, Coon KD, Liang WS, Herbert RH, Beach T, Rohrer KC, Zhao AS, Leung D, Bryden L, Marlowe L, Kaleem M, Mastroeni D, Grover A, Heward CB, Ravid R, Rogers J, Hutton ML, Melquist S, Petersen RC, Alexander GE: GAB2 alleles modify Alzheimer’s risk in APOE epsilon4 carriers. Neuron. 2007, 54: 713-720. 10.1016/j.neuron.2007.05.022.View ArticlePubMedPubMed CentralGoogle Scholar
- Heckerman D, Geiger D, Chickering DM: Learning Bayesian Networks: The Combination of Knowledge and Statistical Data.Mach Learn 1995, 20:197–243.,Google Scholar
- Visweswaran S, Wong A-KI: Bayesian combinatorial partitioning for detecting interactions among genetic variants. Summit Transl Bioinformatics. 2009, 2009: 133-Google Scholar
- Jiang X, Neapolitan RE, Barmada MM, Visweswaran S: Learning genetic epistasis using Bayesian network scoring criteria. BMC Bioinformatics. 2011, 12: 89-10.1186/1471-2105-12-89.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang X, Barmada MM, Visweswaran S: Identifying genetic interactions in genome-wide data using Bayesian networks. Genet Epidemiol. 2010, 34: 575-581. 10.1002/gepi.20514.View ArticlePubMedPubMed CentralGoogle Scholar
- Ritchie MD: Using biological knowledge to uncover the mystery in the search for epistasis in genome-wide association studies. Ann Hum Genet. 2011, 75: 172-182. 10.1111/j.1469-1809.2010.00630.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Saccone SF, Quan J, Jones PL: BioQ: tracing experimental origins in public genomic databases using a novel data provenance model. Bioinformatics. 2012, 28: 1189-1191. 10.1093/bioinformatics/bts117.View ArticlePubMedPubMed CentralGoogle Scholar
- Carmona-Saez P, Chagoyen M, Tirado F, Carazo JM, Pascual-Montano A: GENECODIS: a web-based tool for finding significant concurrent annotations in gene lists. Genome Biol. 2007, 8: R3-10.1186/gb-2007-8-1-r3.View ArticlePubMedPubMed CentralGoogle Scholar
- Nogales-Cadenas R, Carmona-Saez P, Vazquez M, Vicente C, Yang X, Tirado F, Carazo JM, Pascual-Montano A: GeneCodis: interpreting gene lists through enrichment analysis and integration of diverse biological information. Nucleic Acids Res. 2009, 37 (Web Server issue): W317-W322. 10.1093/nar/gkp416.View ArticlePubMedPubMed CentralGoogle Scholar
- Khatri P, Drăghici S: Ontological analysis of gene expression data: current tools, limitations, and open problems. Bioinformatics. 2005, 21: 3587-3595. 10.1093/bioinformatics/bti565.View ArticlePubMedPubMed CentralGoogle Scholar
- Supek F, Bošnjak M, Škunca N, Šmuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011, 6: e21800-10.1371/journal.pone.0021800.View ArticlePubMedPubMed CentralGoogle Scholar
- Parkinson H, Sarkans U, Kolesnikov N, Abeygunawardena N, Burdett T, Dylag M, Emam I, Farne A, Hastings E, Holloway E, Kurbatova N, Lukk M, Malone J, Mani R, Pilicheva E, Rustici G, Sharma A, Williams E, Adamusiak T, Brandizi M, Sklyar N, Brazma A: ArrayExpress update–an archive of microarray and high-throughput sequencing-based functional genomics experiments. Nucleic Acids Res. 2011, 39 (Database issue): D1002-D1004. 10.1093/nar/gkq1040.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ioannidis JPA, Boffetta P, Little J, O’Brien TR, Uitterlinden AG, Vineis P, Balding DJ, Chokkalingam A, Dolan SM, Flanders WD, Higgins JPT, McCarthy MI, McDermott DH, Page GP, Rebbeck TR, Seminara D, Khoury MJ: Assessment of cumulative evidence on genetic associations: interim guidelines. Int J Epidemiol. 2008, 37: 120-132. 10.1093/ije/dym159.View ArticlePubMedGoogle Scholar
- Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O’Donnell CJ, de Bakker PIW: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008, 24: 2938-2939. 10.1093/bioinformatics/btn564.View ArticlePubMedPubMed CentralGoogle Scholar
- Morgan K: The three new pathways leading to Alzheimer’s disease. Neuropathol Appl Neurobiol. 2011, 37: 353-357. 10.1111/j.1365-2990.2011.01181.x.View ArticlePubMedGoogle Scholar
- Lukk M, Kapushesky M, Nikkilä J, Parkinson H, Goncalves A, Huber W, Ukkonen E, Brazma A: A global map of human gene expression. Nat Biotechnol. 2010, 28: 322-324. 10.1038/nbt0410-322.View ArticlePubMedPubMed CentralGoogle Scholar
- Crews L, Adame A, Patrick C, Delaney A, Pham E, Rockenstein E, Hansen L, Masliah E: Increased BMP6 levels in the brains of Alzheimer’s disease patients and APP transgenic mice are accompanied by impaired neurogenesis. J Neurosci. 2010, 30: 12252-12262. 10.1523/JNEUROSCI.1305-10.2010.View ArticlePubMedPubMed CentralGoogle Scholar
- Tesseur I, Zou K, Esposito L, Bard F, Berber E, Can JV, Lin AH, Crews L, Tremblay P, Mathews P, Mucke L, Masliah E, Wyss-Coray T: Deficiency in neuronal TGF-beta signaling promotes neurodegeneration and Alzheimer’s pathology. J Clin Invest. 2006, 116: 3060-3069. 10.1172/JCI27341.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu F, Arias-Vásquez A, Sleegers K, Aulchenko YS, Kayser M, Sanchez-Juan P, Feng B-J, Bertoli-Avella AM, van Swieten J, Axenovich TI, Heutink P, van Broeckhoven C, Oostra BA, van Duijn CM: A genomewide screen for late-onset Alzheimer disease in a genetically isolated Dutch population. Am J Hum Genet. 2007, 81: 17-31. 10.1086/518720.View ArticlePubMedPubMed CentralGoogle Scholar
- Blacker D, Bertram L, Saunders AJ, Moscarillo TJ, Albert MS, Wiener H, Perry RT, Collins JS, Harrell LE, Go RCP, Mahoney A, Beaty T, Fallin MD, Avramopoulos D, Chase GA, Folstein MF, McInnis MG, Bassett SS, Doheny KJ, Pugh EW, Tanzi RE: Results of a high-resolution genome screen of 437 Alzheimer’s disease families. Hum Mol Genet. 2003, 12: 23-32. 10.1093/hmg/ddg007.View ArticlePubMedGoogle Scholar
- Nilsson T, Bogdanovic N, Volkman I, Winblad B, Folkesson R, Benedikz E: Altered subcellular localization of ornithine decarboxylase in Alzheimer’s disease brain. Biochem Biophys Res Commun. 2006, 344: 640-646. 10.1016/j.bbrc.2006.03.191.View ArticlePubMedGoogle Scholar
- Bernstein HG, Müller M: The cellular localization of the L-ornithine decarboxylase/polyamine system in normal and diseased central nervous systems. Prog Neurobiol. 1999, 57: 485-505. 10.1016/S0301-0082(98)00065-3.View ArticlePubMedGoogle Scholar
- Lourenço FC, Galvan V, Fombonne J, Corset V, Llambi F, Müller U, Bredesen DE, Mehlen P: Netrin-1 interacts with amyloid precursor protein and regulates amyloid-beta production. Cell Death Differ. 2009, 16: 655-663. 10.1038/cdd.2008.191.View ArticlePubMedPubMed CentralGoogle Scholar
- Cramer PE, Cirrito JR, Wesson DW, Lee CYD, Karlo JC, Zinn AE, Casali BT, Restivo JL, Goebel WD, James MJ, Brunden KR, Wilson DA, Landreth GE: ApoE-directed therapeutics rapidly clear β-amyloid and reverse deficits in AD mouse models. Science. 2012, 335: 1503-1506. 10.1126/science.1217697.View ArticlePubMedPubMed CentralGoogle Scholar
- Tzschach A, Bisgaard A-M, Kirchhoff M, Graul-Neumann LM, Neitzel H, Page S, Ahmed A, Müller I, Erdogan F, Ropers H-H, Kalscheuer VM, Ullmann R: Chromosome aberrations involving 10q22: report of three overlapping interstitial deletions and a balanced translocation disrupting C10orf11. Eur J Hum Genet. 2010, 18: 291-295. 10.1038/ejhg.2009.163.View ArticlePubMedGoogle Scholar
- Webster JA, Gibbs JR, Clarke J, Ray M, Zhang W, Holmans P, Rohrer K, Zhao A, Marlowe L, Kaleem M, McCorquodale DS, Cuello C, Leung D, Bryden L, Nath P, Zismann VL, Joshipura K, Huentelman MJ, Hu-Lince D, Coon KD, Craig DW, Pearson JV, Heward CB, Reiman EM, Stephan D, Hardy J, Myers AJ: Genetic control of human brain transcript expression in Alzheimer disease. Am J Hum Genet. 2009, 84: 445-458. 10.1016/j.ajhg.2009.03.011.View ArticlePubMedPubMed CentralGoogle Scholar
- Consortium TU: Reorganizing the protein space at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2012, 40 (Database issue): D71-D75. 10.1093/nar/gkr981.View ArticleGoogle Scholar
- Frykman S, Teranishi Y, Hur J-Y, Sandebring A, Goto Yamamoto N, Ancarcrona M, Nishimura T, Winblad B, Bogdanovic N, Schedin-Weiss S, Kihara T, Tjernberg LO: Identification of two novel synaptic γ-secretase associated proteins that affect amyloid β-peptide levels without altering Notch processing. Neurochem Int. 2012, 61: 108-118. 10.1016/j.neuint.2012.03.016.View ArticlePubMedGoogle Scholar
- D’Amelio M, Cavallucci V, Middei S, Marchetti C, Pacioni S, Ferri A, Diamantini A, De Zio D, Carrara P, Battistini L, Moreno S, Bacci A, Ammassari-Teule M, Marie H, Cecconi F: Caspase-3 triggers early synaptic dysfunction in a mouse model of Alzheimer’s disease. Nat Neurosci. 2011, 14: 69-76. 10.1038/nn.2709.View ArticlePubMedGoogle Scholar
- Lipinski MM, Zheng B, Lu T, Yan Z, Py BF, Ng A, Xavier RJ, Li C, Yankner BA, Scherzer CR, Yuan J: Genome-wide analysis reveals mechanisms modulating autophagy in normal brain aging and in Alzheimer’s disease. Proc Natl Acad Sci U S A. 2010, 107: 14164-14169. 10.1073/pnas.1009485107.View ArticlePubMedPubMed CentralGoogle Scholar
- Perry T, Greig NH: Enhancing central nervous system endogenous GLP-1 receptor pathways for intervention in Alzheimer’s disease. Curr Alzheimer Res. 2005, 2: 377-385. 10.2174/1567205054367892.View ArticlePubMedGoogle Scholar
- Dunn CD, Sulis ML, Ferrando AA, Greenwald I: A conserved tetraspanin subfamily promotes Notch signaling in Caenorhabditis elegans and in human cells. Proc Natl Acad Sci U S A. 2010, 107: 5907-5912. 10.1073/pnas.1001647107.View ArticlePubMedPubMed CentralGoogle Scholar
- Maglott D, Ostell J, Pruitt KD, Tatusova T: Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res. 2011, 39 (Database issue): D52-D57. 10.1093/nar/gkq1237.View ArticlePubMedGoogle Scholar
- Wayman GA, Lee Y-S, Tokumitsu H, Silva AJ, Soderling TR: Calmodulin-kinases: modulators of neuronal development and plasticity. Neuron. 2008, 59: 914-931. 10.1016/j.neuron.2008.08.021.View ArticlePubMedPubMed CentralGoogle Scholar
- Pugazhenthi S, Wang M, Pham S, Sze C-I, Eckman CB: Downregulation of CREB expression in Alzheimer’s brain and in Aβ-treated rat hippocampal neurons. Mol Neurodegener. 2011, 6: 60-10.1186/1750-1326-6-60.View ArticlePubMedPubMed CentralGoogle Scholar
- Lambert J-C, Grenier-Boley B, Chouraki V, Heath S, Zelenika D, Fievet N, Hannequin D, Pasquier F, Hanon O, Brice A, Epelbaum J, Berr C, Dartigues J-F, Tzourio C, Campion D, Lathrop M, Amouyel P: Implication of the immune system in Alzheimer’s disease: evidence from genome-wide pathway analysis. J Alzheimers Dis. 2010, 20: 1107-1118.PubMedGoogle Scholar
- Müller M, Cárdenas C, Mei L, Cheung K-H, Foskett JK: Constitutive cAMP response element binding protein (CREB) activation by Alzheimer’s disease presenilin-driven inositol trisphosphate receptor (InsP3R) Ca2+ signaling. Proc Natl Acad Sci U S A. 2011, 108: 13293-13298. 10.1073/pnas.1109297108.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamauchi T: Neuronal Ca2+/calmodulin-dependent protein kinase II–discovery, progress in a quarter of a century, and perspective: implication for learning and memory. Biol Pharm Bull. 2005, 28: 1342-1354. 10.1248/bpb.28.1342.View ArticlePubMedGoogle Scholar
- Riederer BM, Leuba G, Vernay A, Riederer IM: The role of the ubiquitin proteasome system in Alzheimer’s disease. Exp Biol Med (Maywood). 2011, 236: 268-276. 10.1258/ebm.2010.010327.View ArticleGoogle Scholar
- Frame M, Wan KF, Tate R, Vandenabeele P, Pyne NJ: The gamma subunit of the rod photoreceptor cGMP phosphodiesterase can modulate the proteolysis of two cGMP binding cGMP-specific phosphodiesterases (PDE6 and PDE5) by caspase-3. Cell Signal. 2001, 13: 735-741. 10.1016/S0898-6568(01)00193-0.View ArticlePubMedGoogle Scholar
- Tinsley HN, Gary BD, Keeton AB, Lu W, Li Y, Piazza GA: Inhibition of PDE5 by sulindac sulfide selectively induces apoptosis and attenuates oncogenic Wnt/β-catenin-mediated transcription in human breast tumor cells. Cancer Prev Res (Phila). 2011, 4: 1275-1284. 10.1158/1940-6207.CAPR-11-0095.View ArticleGoogle Scholar
- Inestrosa NC, Toledo EM: The role of Wnt signaling in neuronal dysfunction in Alzheimer’s Disease. Mol Neurodegener. 2008, 3: 9-10.1186/1750-1326-3-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Chanock SJ, Manolio T, Boehnke M, Boerwinkle E, Hunter DJ, Thomas G, Hirschhorn JN, Abecasis G, Altshuler D, Bailey-Wilson JE, Brooks LD, Cardon LR, Daly M, Donnelly P, Fraumeni JF, Freimer NB, Gerhard DS, Gunter C, Guttmacher AE, Guyer MS, Harris EL, Hoh J, Hoover R, Kong CA, Merikangas KR, Morton CC, Palmer LJ, Phimister EG, Rice JP, NCI-NHGRI Working Group on Replication in Association Studies: Replicating genotype-phenotype associations. Nature. 2007, 447: 655-660. 10.1038/447655a.View ArticlePubMedGoogle Scholar
- Ioannidis JPA: Non-replication and inconsistency in the genome-wide association setting. Hum Hered. 2007, 64: 203-213. 10.1159/000103512.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.