- Open Access
Genome-wide association study for variants that modulate relationships between cerebrospinal fluid amyloid-beta 42, tau, and p-tau levels
Alzheimer's Research & Therapyvolume 10, Article number: 86 (2018)
A relationship quantitative trait locus exists when the correlation between multiple traits varies by genotype for that locus. Relationship quantitative trait loci (rQTL) are often involved in gene-by-gene (G×G) interactions or gene-by-environmental interactions, making them a powerful tool for detecting G×G.
We performed genome-wide association studies to identify rQTL between tau and Aβ42 and ptau and Aβ42 with over 3000 individuals using age, gender, series, APOE ε2, APOE ε4, and two principal components for population structure as covariates. Each significant rQTL was separately screened for interactions with other loci for each trait in the rQTL model. Parametric bootstrapping was used to assess significance.
We found four significant tau/Aβ42 rQTL from three unique locations and six ptau/Aβ42 rQTL from five unique locations. G×G screens with these rQTL produced four significant G×G interactions (one Aβ42, two ptau, and one tau) with four rQTL where each second locus was from a unique location. On follow-up, rs1036819 and rs74025622 were associated with Alzheimer’s disease (AD) case/control status; rs15205 and rs79099429 were associated with rate of decline.
The two most significant rQTL (rs8027714 and rs1036819) for ptau/Aβ42 are on different chromosomes and both are strong hits for pelvic organ prolapse. While diseases of the nervous system can cause pelvic organ prolapse, it is unlikely related to the ptau/Aβ42 relationship but may suggest that these two loci share a pathway. In addition to a ptau/Aβ42 rQTL and association with AD case/control status, rs1036819 is a strong rQTL for case/control status/Aβ42 and for tau/Aβ42. It resides in the ZFAT gene, which is related to autoimmune thyroid disease. For tau, rs9817620 interacts with the tau/Aβ42 rQTL rs74025622. It is in the CHL1 gene, which is a neural cell adhesion molecule and may be involved in signal transduction pathways. CHL1 is related to BACE1, which is a β-secretase enzyme that initiates production of the β-amyloid peptide involved in AD and is a primary drug target. Overall, there are numerous loci that affect the relationship between these important AD endophenotypes and some are due to interactions with other loci. Some affect the risk of AD and/or rate of progression.
Aβ42, tau, and ptau are important biomarkers for Alzheimer’s disease (AD). Cerebrospinal fluid (CSF) levels of amyloid-beta (Aβ) and tau change before clinical symptoms of AD are observed. Aβ42 is decreased in CSF as the disease progresses, and tau levels increase. This combination of changes appears to be specific to AD . CSF levels of Aβ42 and tau are important indicators of AD pathology. Here, we seek to identify genetic loci that control these levels and modify the relationships between them. A better understanding of the genes and pathways that regulate the relationships between these biomarkers may provide important insights into AD pathology, even at the earliest stages of the disease. Many large genome-wide association studies (GWAS) have been successful in finding loci that affect AD biomarkers such as Aβ42, tau, and ptau. In a series of papers, Hohman et al. [2,3,4] performed a variety of analyses that screened for loci that modify the relationship between various AD-related risk factors. More formally, loci that affect the relationship between two or more traits can be referred to as relationship loci or relationship quantitative trait loci (rQTL) . Differential gene-by-gene (G×G) and/or gene-by-environment (G×E) interactions are often responsible for rQTL . Screening for rQTL is a way to identify important loci that would typically be “invisible” to normal GWAS because they often do not show marginal effects on either trait . It is also a powerful avenue to identify G×G interactions without paying the statistical penalty of performing all possible two-locus tests. rQTL also provide a window into pleiotropy, pleiotropic variability, and selection on traits in complex interconnecting systems .
Here, we present the results of genome-wide screens for rQTL using CSF levels of Aβ42, tau, and ptau subsequent screens for loci that interact with them (G×G). We then present analyses designed to connect the significant rQTL and G×G loci from the primary analyses with AD itself, first by testing the hypothesis that if a locus modulates the relationship between two AD biomarkers (i.e., an rQTL) it may modulate the risk relationship between AD and either biomarker (i.e., an AD/biomarker rQTL), and finally we follow-up with associations of these markers with AD risk and rate of progression (rate of cognitive decline).
These analyses are designed to identify context-dependent effects that are often not “seen” in traditional marginal effect screens (i.e., a typical GWAS). We demonstrate that they exist and that they may point us to other pathways. While interactions can be hard to detect and interpret, they are expected in complex biological systems and they are important as they can create unexplained heterogeneity and can identify subgroups that respond differently to treatment (G×E) or in the context of other genes (G×G).
The data and analyses can be grouped into primary analyses and follow-up. The primary analyses (rQTL and G×G screens) use a large CSF AD biomarker dataset (n = 3146) derived from nine separate studies. The follow up analyses use different datasets to follow-up on the loci found in the primary rQTL and G×G screens. There are three follow-up analyses. The first uses the original CSF biomarker data to determine if any of the significant biomarker rQTL act as rQTL between disease (AD risk) and either biomarker. In the second follow-up, the loci were tested for association with AD risk using a large AD case/control dataset (n = 28,730) from the Alzheimer’s Disease Genetic Consortium (ADGC) and secondarily with results from the International Genomics of Alzheimer’s Project (IGAP) consortium (n = 54,162). The last follow-up is based on the association with rate of progression (cognitive decline) which was performed using data (n = 1499) from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) and the Charles F. and Joanne Knight Alzheimer’s Disease Research Center (Knight ADRC) from Washington University in St. Louis, USA. Direct replication of these results is not possible at this time as there are no other CSF AD biomarker datasets sufficient for replication; we hope to remedy this in the future. Our best approach akin to replication is to establish a connection between these loci and relevant AD endpoints such as risk and rate of progression.
Data and methods for the rQTL and G×G analyses
For the rQTL and G×G analyses, the dataset for CSF biomarker analysis consisted of 3146 participants from nine different studies and is described in detail in our recent publication . Included were 805 individuals (29.34% cases) enrolled in studies at the Knight ADRC, 787 individuals (more than 71% cases) from the ADNI (390 from ADNI1 and 397 from ADNI2), 184 individuals (5.43% cases) from Predictors of Cognitive Decline Among Normal Individuals (BIOCARD), 105 individuals (no AD status) from Saarland University in Homburg/Saar, Germany, 433 individuals (22.17% cases) from the Mayo Clinic, 293 individuals (all cases) from Skåne University Hospital, Sweden, 164 (62.8% cases) from studies at Perelman School of Medicine at the University of Pennsylvania, and 375 (33.33% cases) from studies at the University of Washington. Table 1 shows the demographic data for each study. Clinical assessments, CSF collection, and proteins were measured by each site. Prior to combining data for analyses, CSF levels of tau, p-tau, and Aβ42 were log10-transformed to approximate a normal distribution and the mean from each dataset was standardized to zero to account for the different platforms used by the different studies to measure protein levels. There were no significant differences in the transformed and standardized values for the different studies. Study, age, sex, and the first two principal components were identified as confounding factors by stepwise regression analyses for each protein and corrected for in applicable analyses ; other papers have been published that have combined data in similar ways [10, 11].
rQTL and G×G genetic data: for the rQTL and G×G analyses
Five million single nucleotide polymorphisms (SNPs) (5,088,365) were included for analyses from the 12 million available imputed SNPs after excluding those with a minor allele frequency (MAF) less than 0.01. Data was imputed using the 1000 Genomes Phase 3 data (October 2014) as the reference panel . For the rQTL screens we used a standard 5 × 10−8 genome-wide significance threshold and for the G×G screens we corrected for 5 × 10−8 divided by the number of independent and significant rQTL from each screen (tau/Aβ42 and ptau/Aβ42; see Table 3). Most analyses were performed using the Julia programming language (v0.3.12) [12, 13] and some plots were made in R .
As a strategy to detect G×G interactions, we performed genome-wide single-locus screens for rQTL  with the CSF biomarker data. All models used the same covariates: age, sex, series, APOε2, APOε4, and the first two principle components for ancestry. Significant rQTL were used as independent a-priori hypotheses to identify G×G interactions. Two rQTL models were screened:
where the SNP genotypes are treated as factors and only genotypes with a count of 20 or more were included, leading to models with sufficient genotype counts for tests; either a three-genotype model with 2 degrees of freedom (DF) or a two-genotype model with 1 DF. A putative rQTL exists when the Aβ42*SNP interaction term is significant when compared with a null model excluding that term. Tests of interaction terms sometimes exhibit type I error inflation due to underestimates of the covariance matrix . Work performed by Bůžková et al.  and Voorman et al.  showed that sandwich estimators and the parametric bootstrap can be used to obtain valid p values, with the parametric bootstrap as the gold standard. For tests that initially reached the significance threshold, empirical p values were obtained via 200 million parametric bootstraps . A standard 5 × 10−8 genome-wide significance threshold was used for each rQTL screen.
Context-dependent interactions such as G×E and G×G create single locus patterns such as rQTL. Therefore, each significant rQTL create independent a-priori hypotheses for G×E or G×G interactions. Previous papers on two-stage designs for G×G and G×E studies have shown that each a-priori locus (i.e., loci found in single locus screens) can be treated separately for multiple testing in genome screens for interaction [17,18,19,20,21]. For each significant rQTL, we performed separate screens for loci that interact (G×G) with each rQTL and their specific traits. We corrected for 5 × 10−8 divided by the number of independent and significant rQTL from each screen (tau/Aβ42 and ptau/Aβ42).
Traditionally, in human genetics tests of G×G focus only on the additive-by-additive (AA) contrast while interactions between two biallelic loci can have up to four “orthogonal” contrasts that are traditionally parameterized as AA, additive-by-dominance (AD), dominance-by-additive (DA), and dominance-by-dominance (DD). These other types of interactions exist and would be missed with only an AA test. Unfortunately, in natural populations, allele frequencies are not 50/50 which results in unbalanced and sometimes missing cells (two-locus genotypes). The unbalanced cell counts makes the “orthogonal” contrasts no longer orthogonal and sometimes inestimable when there are missing cells [22,23,24]. We used a simple interaction test that avoids the need for explicit parameterization of the interaction contrasts and only accounts for interactions with sufficient data for estimation. Our full model incorporates covariates and treats the two-locus genotypes as factors similar to a one-way analysis of covariance (ANCOVA). The reduced model includes the covariates and treats the genotypes from each individual locus as factors . The full model accounts for all single and two-locus interaction effects without explicitly parameterizing them while the reduced model explicitly accounts for just the single locus effects. With sufficient counts for all two-locus genotypes the test has 4 DF. All two-locus genotypes with counts less than five were excluded which ensures that estimation of the interaction effect is appropriate while avoiding the issues of inestimable contrasts or trying to determine explicit contrasts for each particular case. Below is an example of the full model:
where the two-locus genotypes are treated as factors, and below is an example of the reduced model where the genotypes of each SNP are treated as factors:
We attempted to connect the significant rQTL and G×G loci identified in the primary biomarker analyses to AD etiology in three different ways. First, we determined if each rQTL between biomarkers acts as an rQTL between AD risk and either of the two biomarkers. This hypothesis extends from the observation that if a locus affects the relationship between two risk factors, it may also modify the relationship between those risk factors and disease. Therefore, each rQTL was fitted for an analogous rQTL model using logistic regression between Alzheimer’s case/control status and that risk factor using the same CSF biomarker data. Second, we used a large AD case/control dataset to assess whether any of the rQTL or G×G loci are directly associated with AD risk. Finally, we use a different dataset to assess if any of these loci are associated with disease progression via rate of cognitive decline.
Data (ADGC and IGAP) and methods for AD risk analyses
We used two datasets for connecting the loci to AD risk through AD case/control data. We used the ADGC as our primary analyses but also report the results from IGAP. While the IGAP study is larger (it includes ADGC), we feel that the ADGC data alone is more homogenous and is a closer match in terms of ethnic composition to our discovery dataset. The ADGC is an NIH-funded collection of GWAS data created for the goal of identifying genetic contributions to late-onset AD. Participants included in this study were from 30 merged datasets  combined by Boehme et al.  and included 25,666 unrelated individuals carrying either an AD (n = 12,532) or cognitively normal control clinical diagnosis (n = 13,134). Participants were recruited and seen between 1984 and 2012. The genetic data are based on merging the different imputed ADGC datasets. Before merging, SNPs within each dataset with low imputation quality (info < 0.5) were filtered out. Dosage information was used to convert to best guess PLINK allele call format files with an uncertainty cutoff of 0.1. Duplicate samples were identified and removed. Other quality control issues such as physical location discrepancies and strand flipping were identified and appropriately dealt with. A kinship coefficient of 0.0442 was used to screen out third-degree relatives or closer. The resulting number of SNPs in the unrelated dataset was ~ 8.6 million. We only used the SNPs corresponding to our significant rQTL and G×G loci, as simple logistic regression analysis was performed to assess the association between our SNPs and case control status using age, sex, cohort, and two principle components as covariates. Full details on the datasets and the merging process are available at .
The IGAP consortium organized AD case/control data from many different studies across the world. For stage 1 of their analysis they used meta-analysis methods for the association between ~ 7 million SNPs and AD case/control data (n = 54,162; 17,008 cases, 37,154 control) from four studies (ADGC, CHARGE, EADI, and GERAD) . We used the results from the stage 1 analyses to determine whether any of the significant rQTL or G×G loci are associated with AD risk. In the preparation of the genetic data, Lambert et al.  excluded SNPs with call rates < 95%, imputed using samples of European ancestry in the 1000 Genome Project, and excluded SNPs with MAF < 1%. In each case/control dataset, the association of late-onset Alzheimer’s disease with genotype dosage was analyzed by a logistic regression model including covariates for age, sex, and principal components to adjust for possible population stratification. For the three CHARGE cohorts with incident Alzheimer’s disease data, Cox proportional hazards models were used. For the meta-analysis they undertook a fixed-effects inverse variance-weighted meta-analysis with the standard errors of the beta coefficient scaled by the square roots of study-specific genomic inflation factors estimated before combining the summary statistics across datasets .
About half of the samples from the CSF AD biomarker data can be found in both the ADGC and IGAP datasets and only represent ~ 5.5% and ~ 3% of the ADGC and IGAP samples. It is also important to note that the analyses leading to these follow-up tests for association with case/control status in ADGC and IGAP are based on rQTL and G×G analyses using Aβ42, tau, and ptau and not any explicit single-locus test for association with case/control status in the CSF samples.
Data and methods for rate of progression analyses
A total of 1499 individuals with longitudinal cognitive data and genetic data were used to assess the association of the significant rQTL and G×G loci with rate of progression or the rate of cognitive decline. More extensive descriptions of the data and methods are described elsewhere . Participants from this study were enrolled in two different longitudinal studies: the Knight ADRC at Washington University (n = 778) and the ADNI (n = 712). Participants were evaluated in accordance with the Clinical Dementia Rating (CDR) scale, where 0 indicates cognitive normality, 0.5 is very mild dementia, 1 is mild dementia, 2 is moderate dementia, and 3 is severe dementia [30, 31]. The scores in each of the six areas are summed to yield a sum of box (CRD-SB) scores ranging from 0 (no impairment) to 18 (maximal impairment). Only participants that had an AD diagnosis and CDR > 0 at their last visit were included in our analyses. Individuals with dementia caused by neurological diseases other than AD were excluded. For those samples with data available, individuals with Aβ42 values equal or greater than 192 pg/mL (ADNI)  or 500 pg/mL (WU)  were also excluded. Noninformative longitudinally measured CDR-SB was removed for each participant (noninformative longitudinal data is defined as data where the CDR-SB is either 0 or 18 and remains constant over a period of time). Only individuals with at least two visits and 1.5 years of follow-up were included.
Participants were genotyped with the Illumina 610 or Omniexpress chip (Illumina, San Diego, CA, USA). IMPUT2 v2.3.2 software and the 1000 genome (phase3 NCBI build 37) data were used as reference to impute up to 6 million SNPs. To avoid the possibility of spurious associations of population structure and to confirm ethnicity of each sample, the two principal components scores were used as covariates in the analysis. Only individuals that clustered with the European-American cluster were included in the study.
A linear mixed-model analysis was carried out using R statistical software  and the nlme package . A linear mixed-model repeated measure framework was used to account for correlation between repeated measures in the same individual. Change in CDR-SB per year was treated as the independent variable including the following covariants: baseline CDR, baseline age, gender, time (follow-up), level of education, the interaction between baseline CDR and time, and, to avoid the possibility of spurious association due to population substructure, the two first principal components scores were included as covariates, and a random effect for time and individual was included in the model with an AR(1) covariance structure.
From our rQTL analyses we identified four genome-wide significant tau/Aβ42 rQTL from three unique locations and six ptau/Aβ42 rQTL from five unique locations (Table 2). Using the significant rQTL (one from each unique location) as a priori hypotheses, subsequent G×G screens identified loci from four different regions involved in significant two-locus interactions with four of the rQTL (one Aβ42, one tau, and two ptau) (Table 3, Fig. 1) (see “Data and methods for the rQTL and G×G Analyses” in Methods). The respective threshold for G×G screens with the tau/Aβ42 rQTL was p < 1.67 × 10−8 and p < 1 × 10−8 for G×G screens with the ptau/Aβ42 rQTL.
In follow-up, two ptau/Aβ42 rQTL were significantly associated with AD case/control status in the ADGC dataset (rs1036819, p = 1.9 × 10−5; rs74025622, p = 0.044); however, they were not significant in IGAP (rs1036819, p = 0.091; rs74025622, p = 0.427) (see the first few paragraphs of “Follow-up” in Methods). In addition, rs1036819 is also a very strong rQTL for both AD case/control status/Aβ42 (p = 6.2 × 10−8) and tau/Aβ42 (p = 6.17 × 10−5), while rs74025622 is involved in G×G interactions with three other loci (Table 3) (see the first paragraph of “Follow up” in Methods). While not significant in ADGC, rs1558634 (a ptau G×G locus; Table 3) is significantly associated with AD risk in IGAP (p = 0.0058). Two loci (rs15205, p = 0.019; rs79099429, p = 0.034) are significant for rate of decline (see “Data and methods for Rate of Progression Analyses” in Methods). rs15205 is a ptau/Aβ42 rQTL (Table 2) involved in a ptau G×G interaction and rs79099429 is a ptau/Aβ42 rQTL (Table 2). Plots similar to Fig. 1 for each interaction in Table 3 can be found in Additional file 1 (Figures S1–S3) along with brief comments. Further annotations for the variants in Tables 2 and 3 can be found in Additional file 2 (Tables S1 and S2).
Discussion and conclusions
From our analyses, we have found evidence for rQTL that modulate the relationship between AD biomarkers. Furthermore, we have found that some of these rQTL interact with other loci (G×G), and some of these interactions directly contribute to the rQTL pattern. Ultimately, some of these loci affect the relationship between these biomarkers and AD and also directly affect the risk for AD and rate of decline.
The most intriguing result is the rs1036819 ptau/Aβ42 rQTL on chromosome 8 (Table 2), which is the second most significant rQTL found in our study. In follow-up a-priori tests, there is strong evidence that rs1046819 is also an rQTL for Aβ42 and AD case/control status and an rQTL for tau/Aβ42. This means that rs1036819 (or something associated with it nearby) modifies the relationship between Aβ42 and ptau, tau, and AD risk. In the separate and much larger AD risk dataset, rs1036819 is genotypically significantly associated with AD case/control status (p = 1.95 × 10−5) and there is a mild (but not significant) suggestion that it may impact the rate of decline (p = 0.0685). However, rs1036819 is not associated with case/control status in IGAP. It may be associated in ADGC because the CSF samples represent a higher concentration of samples in ADGC (~ 5.5%) than IGAP (~ 3%) or because ADGC represents a more homogeneous sample from populations more closely matching the CSF AD biomarker samples. The ptau/Aβ42 rQTL pattern shows a moderate negative covariate-corrected correlation between ptau and Aβ42 in the common AA and AC genotypes (−0.169 and −0.19), but a dramatically stronger negative correlation in the rare CC homozygote (−0.414). This appears to mirror the within genotype relationships between AD case/control status and Aβ42 where the rare homozygote has a much stronger negative relationship. The rare homozygote is also what drives the direct association with AD case/control status in the much larger independent dataset.
It is not entirely clear why rs1036819 or the region around it is important; however, it is a very strong GWAS hit for pelvic organ prolapse along with the most significant rQTL (rs8027714) in the study . We are not suggesting that pelvic organ prolapse is related to AD or to the ptau/Aβ42 relationship; instead, this suggests that rs1036819 and rs8027714 may share some biological and potentially disease-related pathway. rs1036819 lies in a noncoding RNA (ncRNA) in an intron of the ZFAT gene and has a RegDB score of four, binds to the POL2 and POL24H8 proteins, and alters a Foxf2 regulatory motif. ZFAT is related to autoimmune thyroid disease . This is of interest since various papers have asserted a relationship between the immune system and AD [37,38,39]. In addition, rs8027714 also interacts (G×G) with another locus for Aβ42 (Table 3, Fig. 1).
Another interesting finding involves the tau/Aβ42 rQTL, rs74025622. It interacts with three different loci for tau (Table 3, Additional file 1: Figure S2). It is upstream (5’) of a large GWAS hit for sphingolipids in the CNTNAP4 gene, which has a role in neurotransmission of the dopaminergic and GABAergic systems, and mutations may be related to certain psychiatric illnesses. One particular locus it interacts with, rs9817620, is in the CHL1 gene. The G×G is genome-wide significant for tau, has a nominally significant interaction for Aβ42 (p = 0.0013), and the pattern of interaction differs for the two traits in a way that creates the rQTL pattern of tau/Aβ42 correlation changes across genotypes (see Fig. 1a). CHL1 is a neural cell adhesion molecule and may be involved in signal transduction pathways . CHL1 is related to BACE1 as a substrate , which is a β-secretase enzyme  that initiates production of the β-amyloid peptide involved in AD and is a primary drug target . BACE1-deficient mice showed phenotypic similarities to mice with CHL1 loss of function .
Only two of the 10 rQTL and one of the second G×G loci show any nominal direct association (i.e., a marginal/mean effect) with either of the traits in the rQTL or G×G models. This follows a trend where rQTL usually do not (but can) show a direct marginal (mean) effect with any of the traits in the rQTL making them invisible to typical association studies . Theoretically, rQTL are often due to interactions with other loci or environments . Specifically, differential epistasis occurs when the pattern of interaction between two loci differs across traits leading to trait1/trait2 correlation heterogeneity across the genotypes of one or both of the interacting loci. While interactions were not found for all rQTL, a number were found and, of those, many had a pattern of interaction with the two traits that directly contribute to the rQTL pattern (see Fig. 1 and Additional file 1: Figures S1–S3).
The greatest limitation of this study is the lack of an appropriate CSF biomarker study for replication of the rQTL and G×G. The purpose of this work is to generate hypotheses for future, large CSF biomarker studies to test. Consequently, we are not making grand assertions about our findings and our group is actively working to produce and assemble another CSF dataset that can be used for validation in the future. In fact, we are in a similar situation to the original authors of the recent paper  that performed the first analysis of the rate of progression in the ADNI and Knight ADRC data. They did not have a validation set and stated that their findings were meant to create hypotheses for future studies. Here we are using a small subset of their data and methods to connect our findings to make meaningful connections to disease progression. Subsequently, we have shown that some of these loci directly relate to AD through associations with AD risk and with rate of progression.
While context-dependent effects are hard to detect and to interpret, they are expected in biological systems and can help explain some of the heterogeneity among individuals. One hope is that with interactions we can identify subgroups that act differently from the “norm” and that these “exceptions” may help us to understand the more general responses or lead to different biological pathways.
- Aβ :
Alzheimer’s Disease Genetic Consortium
Alzheimer’s Disease Network Initiative
Alzheimer’s Disease Research Center
Clinical Dementia Rating
Degrees of freedom
Genome-wide association studies
International Genomics of Alzheimer’s Project
Minor allele frequency
Relationship quantitative trait loci
single nucleotide polymorphism
Jack CR, Knopman DS, Jagust WJ, Petersen RC, Weiner MW, Aisen PS, et al. Tracking pathophysiological processes in Alzheimer’s disease: an updated hypothetical model of dynamic biomarkers. Lancet Neurol. 2013;12:207–16.
Hohman TJ, Koran MEI, Thornton-Wells TA, Alzheimer’s Disease Neuroimaging Initiative. Genetic modification of the relationship between phosphorylated tau and neurodegeneration. Alzheimers Dement. 2014;10:637–645.e1.
Hohman TJ, Koran MEI, Thornton-Wells TA, Alzheimer’s Neuroimaging Initiative. Genetic variation modifies risk for neurodegeneration based on biomarker status. Front Aging Neurosci. 2014;6:183.
Hohman TJ, Chibnik L, Bush WS, Jefferson AL, De Jaeger PL, Thornton-Wells TA, et al. GSK3β interactions with amyloid genes: an autopsy verification and extension. Neurotox Res. 2015;28:232–8.
Maxwell TJ, Ballantyne CM, Cheverud JM, Guild CS, Ndumele CE, Boerwinkle E. APOE modulates the correlation between triglycerides, cholesterol, and CHD through pleiotropy, and gene-by-gene interactions. Genetics. 2013;195:1397–405.
Pavlicev M, Kenney-Hunt JP, Norgard EA, Roseman CC, Wolf JB, Cheverud JM. Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight. Evolution. 2008;62:199–213.
Pavlicev M, Norgard EA, Fawcett GL, Cheverud JM. Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map. J Exp Zool B Mol Dev Evol. 2011;316:371–85.
Pavlicev M, Cheverud JM, Wagner GP. Evolution of adaptive phenotypic variation patterns by direct selection for evolvability. Proc Biol Sci. 2011;278:1903–12.
Deming Y, Li Z, Kapoor M, Harari O, Del-Aguila JL, Black K, et al. Genome-wide association study identifies four novel loci associated with Alzheimer’s endophenotypes and disease modifiers. Acta Neuropathol. 2017;133:839–56.
Cruchaga C, Kauwe JSK, Harari O, Jin SC, Cai Y, Karch CM, et al. GWAS of cerebrospinal fluid tau levels identifies risk variants for Alzheimer’s disease. Neuron. 2013;78:256–68.
Cruchaga C, Kauwe JSK, Nowotny P, Bales K, Pickering EH, Mayo K, et al. Cerebrospinal fluid APOE levels: an endophenotype for genetic studies for Alzheimer’s disease. Hum Mol Genet. 2012;21:4558–71.
Bezanson J, Edelman A, Karpinski S, Shah V. Julia: a fresh approach to numerical computing. SIAM Rev. 2017;59:65–98.
Bezanson J, Karpinski S, Shah VB, Edelman A. Julia: A fast dynamic language for technical computing. arXiv:1209.5145 [cs] [Internet]. 2012 [cited 2017 Oct 2]; Available from: http://arxiv.org/abs/1209.5145
R Core Development Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017. Available from: http://www.R-project.org/
Voorman A, Lumley T, McKnight B, Rice K. Behavior of QQ-plots and genomic control in studies of gene-environment interaction. PLoS One. 2011;6:e19416.
Bůžková P, Lumley T, Rice K. Permutation and parametric bootstrap tests for gene-gene and gene-environment interactions. Ann Hum Genet. 2011;75:36–45.
Kooperberg C, Leblanc M. Increasing the power of identifying gene x gene interactions in genome-wide association studies. Genet Epidemiol. 2008;32:255–63.
Wei W-H, Knott S, Haley CS, de Koning D-J. Controlling false positives in the mapping of epistatic QTL. Heredity (Edinb). 2010;104:401–9.
Dai JY, Kooperberg C, Leblanc M, Prentice RL. Two-stage testing procedures with independent filtering for genome-wide gene-environment interaction. Biometrika. 2012;99:929–44.
Hsu L, Jiao S, Dai JY, Hutter C, Peters U, Kooperberg C. Powerful cocktail methods for detecting genome-wide gene-environment interaction. Genet Epidemiol. 2012;36:183–94.
Gauderman WJ, Zhang P, Morrison JL, Lewinger JP. Finding novel genes by testing G × E interactions in a genome-wide association study. Genet Epidemiol. 2013;37:603–13.
Cheverud JM, Routman EJ. Epistasis and its contribution to genetic variance components. Genetics. 1995;139:1455–61.
Routman EJ, Cheverud JM. Gene effects on a quantitative trait: two-locus epistatic effects measured at microsatellite markers and at estimated QTL. Evolution. 1997;51:1654–62.
Hamon SC, Stengard JH, Clark AG, Salomaa V, Boerwinkle E, Sing CF. Evidence for non-additive influence of single nucleotide polymorphisms within the apolipoprotein E gene. Ann Hum Genet. 2004;68:521–35.
Hemani G, Shakhbazov K, Westra H-J, Esko T, Henders AK, McRae AF, et al. Detection and replication of epistasis influencing transcription in humans. Nature. 2014;508:249–53.
Steele NZR, Carr JS, Bonham LW, Geier EG, Damotte V, Miller ZA, et al. Fine-mapping of the human leukocyte antigen locus as a risk factor for Alzheimer disease: a case-control study. PLoS Med. 2017;14:e1002272.
Boehme KL, Mukherjee S, Crane PK, Kauwe JS. ADGC 1000 Genomes combined data workflow. 2014 found at: http://kauwelab.byu.edu/Portals/22/adgc_combined_1000G_09192014.pdf.
Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45:1452–8.
Del-Aguila JL, Fernández MV, Schindler S, Ibanez L, Deming Y, Ma S, et al. Assessment of the genetic architecture of Alzheimer’s disease risk in rate of memory decline. J Alzheimers Dis. 2018;62:745–56.
Morris JC, McKeel DW, Fulling K, Torack RM, Berg L. Validation of clinical diagnostic criteria for Alzheimer’s disease. Ann Neurol. 1988;24:17–22.
Morris JC. The clinical dementia rating (CDR): current version and scoring rules. Neurology. 1993;43:2412–4.
Shaw LM, Vanderstichele H, Knapik-Czajka M, Clark CM, Aisen PS, Petersen RC, et al. Cerebrospinal fluid biomarker signature in Alzheimer’s disease neuroimaging initiative subjects. Ann Neurol. 2009;65:403–13.
Morris JC, Roe CM, Xiong C, Fagan AM, Goate AM, Holtzman DM, et al. APOE predicts amyloid-beta but not tau Alzheimer pathology in cognitively normal aging. Ann Neurol. 2010;67:122–31.
Pinheiro, J, Bates, D, Debroy, S, R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1–131. 2017; Available from: https://cran.r-project.org/web/packages/nlme/index.html.
Allen-Brady K, Cannon-Albright L, Farnham JM, Teerlink C, Vierhout ME, van Kempen LCL, et al. Identification of six loci associated with pelvic organ prolapse using genome-wide association analysis. Obstet Gynecol. 2011;118:1345–53.
Nakabayashi K, Shirasawa S. Recent advances in the association studies of autoimmune thyroid disease and the functional characterization of AITD-related transcription factor ZFAT. Nihon Rinsho Meneki Gakkai Kaishi. 2010;33:66–72.
Gagliano SA, Pouget JG, Hardy J, Knight J, Barnes MR, Ryten M, et al. Genomics implicates adaptive and innate immunity in Alzheimer’s and Parkinson’s diseases. Ann Clin Transl Neurol. 2016;3:924–33.
Hu Y-S, Xin J, Hu Y, Zhang L, Wang J. Analyzing the genes related to Alzheimer’s disease via a network and pathway-based approach. Alzheimers Res Ther. 2017;9:29.
International Genomics of Alzheimer’s Disease Consortium (IGAP). Convergent genetic and expression data implicate immunity in Alzheimer’s disease. Alzheimers Dement. 2015;11:658–71.
Zhou L, Barão S, Laga M, Bockstael K, Borgers M, Gijsen H, et al. The neural cell adhesion molecules L1 and CHL1 are cleaved by BACE1 protease in vivo. J Biol Chem. 2012;287:25927–40.
Kizuka Y, Kitazume S, Fujinawa R, Saito T, Iwata N, Saido TC, et al. An aberrant sugar modification of BACE1 blocks its lysosomal targeting in Alzheimer’s disease. EMBO Mol Med. 2015;7:175–89.
Barão S, Gärtner A, Leyva-Díaz E, Demyanenko G, Munck S, Vanhoutvin T, et al. Antagonistic effects of BACE1 and APH1B-γ-secretase control axonal guidance by regulating growth cone collapse. Cell Rep. 2015;12:1367–76.
Yan R. Stepping closer to treating Alzheimer’s disease patients with BACE1 inhibitor drugs. Transl Neurodegener. 2016;5:13.
Hitt B, Riordan SM, Kukreja L, Eimer WA, Rajapaksha TW, Vassar R. β-Site amyloid precursor protein (APP)-cleaving enzyme 1 (BACE1)-deficient mice exhibit a close homolog of L1 (CHL1) loss-of-function phenotype involving axon guidance defects. J Biol Chem. 2012;287:38408–25.
The authors gratefully acknowledge the study subjects for their participation and the many researchers whose efforts contributed to the assembly and characterization of these datasets.
The primary analyses, methodological development, and writing were supported by NIA/NIH grant 1R01AG042611 to JSKK, and NHLBI/NIH grant HL105502 to TJM. NIH grants R01AG044546, P01AG003991, RF1AG053303, R01AG035083, and R01NS085419 supported the data.
Availability of data and materials
Data related to this study can be accessed via NIAGADS (https://www.niagads.org), ADNI (http://adni.loni.usc.edu), and dbGAP (https://www.ncbi.nlm.nih.gov/gap).
Ethics approval and consent to participate
The work reported here was conducted under approval of the Institutional Review Boards at Brigham Young University, Utah State University, and the Washington University School of Medicine.
Consent for publication
This manuscript does not contain any individual person’s data in any form.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.