We aimed to improve on the single-atlas ventricular segmentation method of (Carmichael, O.T., Thompson, P.M., Dutton, R.A., Lu, A., Lee, S.E., Lee, J.Y., Kuller, L.H., Lopez, O.L., Aizenstein, H.J., Meltzer, C.C., Liu, Y., Toga, A.W., Becker, J.T., 2006. Mapping ventricular changes related to dementia and mild cognitive impairment in a large community-based cohort. IEEE ISBI. 315-318) by using multi-atlas segmentation, which has been shown to lead to more accurate segmentations (Chou, Y., Leporé, N., de Zubicaray, G., Carmichael, O., Becker, J., Toga, A., Thompson, P., 2008. Automated ventricular mapping with multi-atlas fluid image alignment reveals genetic effects in Alzheimer's disease, NeuroImage 40(2): 615-630); with this method, we calculated minimal numbers of subjects needed to detect correlations between clinical scores and ventricular maps. We also assessed correlations between emerging CSF biomarkers of Alzheimer's disease pathology and localizable deficits in the brain, in 80 AD, 80 mild cognitive impairment (MCI), and 80 healthy controls from the Alzheimer's Disease Neuroimaging Initiative. Six expertly segmented images and their embedded parametric mesh surfaces were fluidly registered to each brain; segmentations were averaged within subjects to reduce errors. Surface-based statistical maps revealed powerful correlations between surface morphology and 4 variables: (1) diagnosis, (2) depression severity, (3) cognitive function at baseline, and (4) future cognitive decline over the following year. Cognitive function was assessed using the mini-mental state exam (MMSE), global and sum-of-boxes clinical dementia rating (CDR) scores, at baseline and 1-year follow-up. Lower CSF Abeta(1-42) protein levels, a biomarker of AD pathology assessed in 138 of the 240 subjects, were correlated with lateral ventricular expansion. Using false discovery rate (FDR) methods, 40 and 120 subjects, respectively, were needed to discriminate AD and MCI from normal groups. 120 subjects were required to detect correlations between ventricular enlargement and MMSE, global CDR, sum-of-boxes CDR and clinical depression scores. Ventricular expansion maps correlate with pathological and cognitive measures in AD, and may be useful in future imaging-based clinical trials.

Rates of brain atrophy derived from serial magnetic resonance (MR) studies may be used to assess therapies for Alzheimer's disease (AD). These measures may be confounded by changes in scanner voxel sizes. For this reason, the Alzheimer's Disease Neuroimaging Initiative (ADNI) included the imaging of a geometric phantom with every scan. This study compares voxel scaling correction using a phantom with correction using a 9 degrees of freedom (9DOF) registration algorithm. We took 129 pairs of baseline and 1-year repeat scans, and calculated the volume scaling correction, previously measured using the phantom. We used the registration algorithm to quantify any residual scaling errors, and found the algorithm to be unbiased, with no significant (p=0.97) difference between control (n=79) and AD subjects (n=50), but with a mean (SD) absolute volume change of 0.20 (0.20) % due to linear scalings. 9DOF registration was shown to be comparable to geometric phantom correction in terms of the effect on atrophy measurement and unbiased with respect to disease status. These results suggest that the additional expense and logistic effort of scanning a phantom with every patient scan can be avoided by registration-based scaling correction. Furthermore, based upon the atrophy rates in the AD subjects in this study, sample size requirements would be approximately 10-12% lower with (either) correction for voxel scaling than if no correction was used.

Despite the central role of amyloid deposition in the development of Alzheimer's disease (AD), the pathogenesis of AD still remains elusive at the molecular level. Increasing evidence suggests that compromised mitochondrial function contributes to the aging process and thus may increase the risk of AD. Dysfunctional mitochondria contribute to reactive oxygen species (ROS) which can lead to extensive macromolecule oxidative damage and the progression of amyloid pathology. Oxidative stress and amyloid toxicity leave neurons chemically vulnerable. Because the brain relies on aerobic metabolism, it is apparent that mitochondria are critical for the cerebral function. Mitochondrial DNA sequence changes could shift cell dynamics and facilitate neuronal vulnerability. Therefore we postulated that mitochondrial DNA sequence polymorphisms may increase the risk of AD. We evaluated the role of mitochondrial haplogroups derived from 138 mitochondrial polymorphisms in 358 Caucasian Alzheimer's Disease Neuroimaging Initiative (ADNI) subjects. Our results indicate that the mitochondrial haplogroup UK may confer genetic susceptibility to AD independently of the apolipoprotein E4 (APOE4) allele.

The ability to pool data from multiple MRI scanners is becoming increasingly important with the influx in multi-site research studies. Fast spin echo (FSE) dual spin echo sequences are often chosen for such studies based principally on their short acquisition time and the clinically useful contrasts they provide for assessing gross pathology. The practicality of measuring FSE-T2 relaxation properties has rarely been assessed. Here, FSE-T2 relaxation properties are examined across the three main scanner vendors (General Electric (GE), Philips, and Siemens). The American College of Radiology (ACR) phantom was scanned on four 1.5T platforms (two GE, one Philips, and one Siemens) to determine if the dual echo pulse sequence is susceptible to vendor-based variance. In addition, data from 85 subjects spanning the spectrum of normal aging, mild cognitive impairment (MCI), and Alzheimer's disease (AD) was obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) to affirm the presence of any phantom based between vendor variance and determine the relationship between this variance and disease. FSE-T2 relaxation properties, including peak FSE-T2 and histogram width, were calculated for each phantom and human subject. Direct correspondence was found between the phantom and human subject data. Peak FSE-T2 of Siemens scanners was consistently at least 20ms prolonged compared to GE and Philips. Siemens scanners showed broader FSE-T2 histograms than the other scanners. Greater variance was observed across GE scanners than either Philips or Siemens. FSE-T2 differences were much greater with scanner vendor than between diagnostic groups, as no significant changes in peak FSE-T2 or histogram width between normal aged, MCI, and AD subject groups were observed. These results indicate that whole brain histogram measures are not sensitive enough to detect FSE-T2 changes between normal aging, MCI, and AD and that FSE-T2 is highly variable across scanner vendors.

Imaging traits provide a powerful and biologically relevant substrate to examine the influence of genetics on the brain. Interest in genome-wide, brain-wide search for influential genetic variants is growing, but has mainly focused on univariate, SNP-based association tests. Moving to gene-based multivariate statistics, we can test the combined effect of multiple genetic variants in a single test statistic. Multivariate models can reduce the number of statistical tests in gene-wide or genome-wide scans and may discover gene effects undetectable with SNP-based methods. Here we present a gene-based method for associating the joint effect of single nucleotide polymorphisms (SNPs) in 18,044 genes across 31,662 voxels of the whole brain in 731 elderly subjects (mean age: 75.56±6.82SD years; 430 males) from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Structural MRI scans were analyzed using tensor-based morphometry (TBM) to compute 3D maps of regional brain volume differences compared to an average template image based on healthy elderly subjects. Using the voxel-level volume difference values as the phenotype, we selected the most significantly associated gene (out of 18,044) at each voxel across the brain. No genes identified were significant after correction for multiple comparisons, but several known candidates were re-identified, as were other genes highly relevant to brain function. GAB2, which has been previously associated with late-onset AD, was identified as the top gene in this study, suggesting the validity of the approach. This multivariate, gene-based voxelwise association study offers a novel framework to detect genetic influences on the brain.

Weight changes are common in aging and Alzheimer's disease (AD) and postmortem findings suggest a relation between lower body mass index (BMI) and increased AD brain pathology. In the current multicenter study, we tested whether lower BMI is associated with higher core AD brain pathology as assessed by cerebrospinal fluid (CSF)-based biological markers of AD in 751 living subjects: 308 patients with AD, 296 subjects with amnestic mild cognitive impairment (MCI), and 147 elderly healthy controls (HC). Based upon a priori cutoff values on CSF concentration of total tau and beta-amyloid (Aβ(1-42)), subjects were binarized into a group with abnormal CSF biomarker signature (CSF+) and those without (CSF-). Results showed that BMI was significantly lower in the CSF+ when compared with the CSF- group (F = 27.7, df = 746, p < 0.001). There was no interaction between CSF signature and diagnosis or apolipoprotein E (ApoE) genotype. In conclusion, lower BMI is indicative of AD pathology as assessed with CSF-based biomarkers in demented and nondemented elderly subjects.

Recent GWAS studies focused on uncovering novel genetic loci related to AD have revealed associations with variants near CLU, CR1, PICALM and BIN1. In this study, we conducted a genome-wide association study in an independent set of 1034 cases and 1186 controls using the Illumina genotyping platforms. By coupling our data with available GWAS datasets from the ADNI and GenADA, we replicated the original associations in both PICALM (rs3851179) and CR1 (rs3818361). The PICALM variant seems to be non-significant after we adjusted for APOE e4 status. We further tested our top markers in 751 independent cases and 751 matched controls. Besides the markers close to the APOE locus, a marker (rs12989701) upstream of BIN1 locus was replicated and the combined analysis reached genome-wide significance level (p = 5E-08). We combined our data with the published Harold et al. study and meta-analysis with all available 6521 cases and 10360 controls at the BIN1 locus revealed two significant variants (rs12989701, p = 1.32E-10 and rs744373, p = 3.16E-10) in limited linkage disequilibrium (r²  =  0.05) with each other. The independent contribution of both SNPs was supported by haplotype conditional analysis. We also conducted multivariate analysis in canonical pathways and identified a consistent signal in the downstream pathways targeted by Gleevec (P = 0.004 in Pfizer; P = 0.028 in ADNI and P = 0.04 in GenADA). We further tested variants in CLU, PICALM, BIN1 and CR1 for association with disease progression in 597 AD patients where longitudinal cognitive measures are sufficient. Both the PICALM and CLU variants showed nominal significant association with cognitive decline as measured by change in Clinical Dementia Rating-sum of boxes (CDR-SB) score from the baseline but did not pass multiple-test correction. Future experiments will help us better understand potential roles of these genetic loci in AD pathology.

This article introduces a hypometabolic convergence index (HCI) for the assessment of Alzheimer's disease (AD); compares it to other biological, cognitive and clinical measures; and demonstrates its promise to predict clinical decline in mild cognitive impairment (MCI) patients using data from the AD Neuroimaging Initiative (ADNI). The HCI is intended to reflect in a single measurement the extent to which the pattern and magnitude of cerebral hypometabolism in an individual's fluorodeoxyglucose positron emission tomography (FDG-PET) image correspond to that in probable AD patients, and is generated using a fully automated voxel-based image-analysis algorithm. HCIs, magnetic resonance imaging (MRI) hippocampal volume measurements, cerebrospinal fluid (CSF) assays, memory test scores, and clinical ratings were compared in 47 probable AD patients, 21 MCI patients who converted to probable AD within the next 18months, 76 MCI patients who did not, and 47 normal controls (NCs) in terms of their ability to characterize clinical disease severity and predict conversion rates from MCI to probable AD. HCIs were significantly different in the probable AD, MCI converter, MCI stable and NC groups (p=9e-17) and correlated with clinical disease severity. Using retrospectively characterized threshold criteria, MCI patients with either higher HCIs or smaller hippocampal volumes had the highest hazard ratios (HRs) for 18-month progression to probable AD (7.38 and 6.34, respectively), and those with both had an even higher HR (36.72). In conclusion, the HCI, alone or in combination with certain other biomarker measurements, has the potential to help characterize AD and predict subsequent rates of clinical decline. More generally, our conversion index strategy could be applied to a range of imaging modalities and voxel-based image-analysis algorithms.

Fluorodeoxyglucose positron emission tomography (FDG-PET) studies report characteristic patterns of cerebral hypometabolism in probable Alzheimer's disease (pAD) and amnestic mild cognitive impairment (aMCI). This study aims to characterize the consistency of regional hypometabolism in pAD and aMCI patients enrolled in the AD neuroimaging initiative (ADNI) using statistical parametric mapping (SPM) and bootstrap resampling, and to compare bootstrap-based reliability index to the commonly used type-I error approach with or without correction for multiple comparisons. Batched SPM5 was run for each of 1000 bootstrap iterations to compare FDG-PET images from 74 pAD and 142 aMCI patients, respectively, to 82 normal controls. Maps of the hypometabolic voxels detected for at least a specific percentage of times over the 1000 runs were examined and compared to an overlap of the hypometabolic maps obtained from 3 randomly partitioned independent sub-datasets. The results from the bootstrap derived reliability of regional hypometabolism in the overall data set were similar to that observed in each of the three non-overlapping sub-sets using family-wise error. Strong but non-linear association was found between the bootstrap-based reliability index and the type-I error. For threshold p=0.0005, pAD was associated with extensive hypometabolic voxels in the posterior cingulate/precuneus and parietotemporal regions with reliability between 90% and 100%. Bootstrap analysis provides an alternative to the parametric family-wise error approach used to examine consistency of hypometabolic brain voxels in pAD and aMCI patients. These results provide a foundation for the use of bootstrap analysis characterize statistical ROIs or search regions in both cross-sectional and longitudinal FDG-PET studies. This approach offers promise in the early detection and tracking of AD, the evaluation of AD-modifying treatments, and other biologically or clinical important measurements using brain images and voxel-based data analysis techniques.

As the most common type of dementia, Alzheimer's disease (AD) is a neurodegenerative disorder initially manifested by impaired memory performances. While the diagnosis information indicates a dichotomous status of a patient, memory scores have the potential to capture the continuous nature of the disease progression and may provide more insights into the underlying mechanism. In this work, we performed a targeted genetic study of memory scores on an AD cohort to identify the associations between a set of genes highly expressed in the hippocampal region and seven cognitive scores related to episodic memory. Both main effects and interaction effects of the targeted genetic markers on these correlated memory scores were examined. In addition to well-known AD genetic markers APOE and TOMM40, our analysis identified a new risk gene NAV2 through the gene-level main effect analysis. NAV2 was found to be significantly and consistently associated with all seven episodic memory scores. Genetic interaction analysis also yielded a few promising hits warranting further investigation, especially for the RAVLT list B Score.

Little is known about how genetic variation contributes to neuroanatomical variability, and whether particular genomic regions comprising genes or evolutionarily conserved elements are enriched for effects that influence brain morphology. Here, we examine brain imaging and single-nucleotide polymorphisms (SNPs) data from ∼2,700 individuals. We show that a substantial proportion of variation in cortical surface area is explained by additive effects of SNPs dispersed throughout the genome, with a larger heritable effect for visual and auditory sensory and insular cortices (h(2)∼0.45). Genome-wide SNPs collectively account for, on average, about half of twin heritability across cortical regions (N=466 twins). We find enriched genetic effects in or near genes. We also observe that SNPs in evolutionarily more conserved regions contributed significantly to the heritability of cortical surface area, particularly, for medial and temporal cortical regions. SNPs in less conserved regions contributed more to occipital and dorsolateral prefrontal cortices.

We introduce a mass-univariate framework for the analysis of whole-brain structural trajectories using longitudinal Voxel-Based Morphometry data and Bayesian inference. Our approach to developmental and aging longitudinal studies characterizes heterogeneous structural growth/decline between and within groups. In particular, we propose a probabilistic generative model that parameterizes individual and ensemble average changes in brain structure using linear mixed-effects models of age and subject-specific covariates. Model inversion uses Expectation Maximization (EM), while voxelwise (empirical) priors on the size of individual differences are estimated from the data. Bayesian inference on individual and group trajectories is realized using Posterior Probability Maps (PPM). In addition to parameter inference, the framework affords comparisons of models with varying combinations of model order for fixed and random effects using model evidence. We validate the model in simulations and real MRI data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) project. We further demonstrate how subject specific characteristics contribute to individual differences in longitudinal volume changes in healthy subjects, Mild Cognitive Impairment (MCI), and Alzheimer's Disease (AD).

More and more large-scale imaging genetic studies are being widely conducted to collect a rich set of imaging, genetic, and clinical data to detect putative genes for complexly inherited neuropsychiatric and neurodegenerative disorders. Several major big-data challenges arise from testing genome-wide (NC>12 million known variants) associations with signals at millions of locations (NV~10(6)) in the brain from thousands of subjects (n~10(3)). The aim of this paper is to develop a Fast Voxelwise Genome Wide Association analysiS (FVGWAS) framework to efficiently carry out whole-genome analyses of whole-brain data. FVGWAS consists of three components including a heteroscedastic linear model, a global sure independence screening (GSIS) procedure, and a detection procedure based on wild bootstrap methods. Specifically, for standard linear association, the computational complexity is O (nNVNC) for voxelwise genome wide association analysis (VGWAS) method compared with O ((NC+NV)n(2)) for FVGWAS. Simulation studies show that FVGWAS is an efficient method of searching sparse signals in an extremely large search space, while controlling for the family-wise error rate. Finally, we have successfully applied FVGWAS to a large-scale imaging genetic data analysis of ADNI data with 708 subjects, 193,275voxels in RAVENS maps, and 501,584 SNPs, and the total processing time was 203,645s for a single CPU. Our FVGWAS may be a valuable statistical toolbox for large-scale imaging genetic analysis as the field is rapidly advancing with ultra-high-resolution imaging and whole-genome sequencing.

Diffusion imaging can assess the white matter connections within the brain, revealing how neural pathways break down in Alzheimer's disease (AD). We analyzed 3-Tesla whole-brain diffusion-weighted images from 202 participants scanned by the Alzheimer's Disease Neuroimaging Initiative-50 healthy controls, 110 with mild cognitive impairment (MCI) and 42 AD patients. From whole-brain tractography, we reconstructed structural brain connectivity networks to map connections between cortical regions. We tested whether AD disrupts the "rich club" - a network property where high-degree network nodes are more interconnected than expected by chance. We calculated the rich club properties at a range of degree thresholds, as well as other network topology measures including global degree, clustering coefficient, path length, and efficiency. Network disruptions predominated in the low-degree regions of the connectome in patients, relative to controls. The other metrics also showed alterations, suggesting a distinctive pattern of disruption in AD, less pronounced in MCI, targeting global brain connectivity, and focusing on more remotely connected nodes rather than the central core of the network. AD involves severely reduced structural connectivity; our step-wise rich club coefficients analyze points to disruptions predominantly in the peripheral network components; other modalities of data are needed to know if this indicates impaired communication among non rich club regions. The highly connected core was relatively preserved, offering new evidence on the neural basis of progressive risk for cognitive decline.

We investigated fluctuations in brain volume throughout the day using statistical modeling of magnetic resonance imaging (MRI) from large populations. We applied fully automated image analysis software to measure the brain parenchymal fraction (BPF), defined as the ratio of the brain parenchymal volume and intracranial volume, thus accounting for variations in head size. The MRI data came from serial scans of multiple sclerosis (MS) patients in clinical trials (n=755, 3269 scans) and from subjects participating in the Alzheimer's Disease Neuroimaging Initiative (ADNI, n=834, 6114 scans). The percent change in BPF was modeled with a linear mixed effect (LME) model, and the model was applied separately to the MS and ADNI datasets. The LME model for the MS datasets included random subject effects (intercept and slope over time) and fixed effects for the time-of-day, time from the baseline scan, and trial, which accounted for trial-related effects (for example, different inclusion criteria and imaging protocol). The model for ADNI additionally included the demographics (baseline age, sex, subject type [normal, mild cognitive impairment, or Alzheimer's disease], and interaction between subject type and time from baseline). There was a statistically significant effect of time-of-day on the BPF change in MS clinical trial datasets (-0.180 per day, that is, 0.180% of intracranial volume, p=0.019) as well as the ADNI dataset (-0.438 per day, that is, 0.438% of intracranial volume, p<0.0001), showing that the brain volume is greater in the morning. Linearly correcting the BPF values with the time-of-day reduced the required sample size to detect a 25% treatment effect (80% power and 0.05 significance level) on change in brain volume from 2 time-points over a period of 1year by 2.6%. Our results have significant implications for future brain volumetric studies, suggesting that there is a potential acquisition time bias that should be randomized or statistically controlled to account for the day-to-day brain volume fluctuations.

In this paper we present a method to segment four brainstem structures (midbrain, pons, medulla oblongata and superior cerebellar peduncle) from 3D brain MRI scans. The segmentation method relies on a probabilistic atlas of the brainstem and its neighboring brain structures. To build the atlas, we combined a dataset of 39 scans with already existing manual delineations of the whole brainstem and a dataset of 10 scans in which the brainstem structures were manually labeled with a protocol that was specifically designed for this study. The resulting atlas can be used in a Bayesian framework to segment the brainstem structures in novel scans. Thanks to the generative nature of the scheme, the segmentation method is robust to changes in MRI contrast or acquisition hardware. Using cross validation, we show that the algorithm can segment the structures in previously unseen T1 and FLAIR scans with great accuracy (mean error under 1mm) and robustness (no failures in 383 scans including 168 AD cases). We also indirectly evaluate the algorithm with a experiment in which we study the atrophy of the brainstem in aging. The results show that, when used simultaneously, the volumes of the midbrain, pons and medulla are significantly more predictive of age than the volume of the entire brainstem, estimated as their sum. The results also demonstrate that the method can detect atrophy patterns in the brainstem structures that have been previously described in the literature. Finally, we demonstrate that the proposed algorithm is able to detect differential effects of AD on the brainstem structures. The method will be implemented as part of the popular neuroimaging package FreeSurfer.

We investigated the neural basis of hallucinations Alzheimer's disease (AD) by applying voxel-based morphometry (VBM) to anatomical and functional data from the AD Neuroimaging Initiative.

Quantitative assessment of medial temporal lobe atrophy has been proposed as a biomarker for Alzheimer's disease (AD) diagnostic and prognostic in mild cognitive impairment (MCI) due to AD. We present the first results of our high-dimensional morphometry technique, tracking tissue composition, and atrophy changes on T1-weighted magnetic resonance imaging at various time points. We selected 187 control subjects, 17 control subjects having progressed to MCI and/or AD, 178 subjects with stable MCI, 165 subjects with MCI having progressed to AD, and 147 AD subjects from the Alzheimer's Disease Neuroimaging Initiative study. Results show statistically significant differences between almost every diagnostic and time point comparison pairs (0-12, 12-24, and 24-36 months), including controls having progressed to either MCI or AD and trajectory dynamics that demonstrate the algorithm's ability at tracking specific pathology-related neurodegeneration.

Several studies support the relation between leptin and Alzheimer's disease (AD). We show that leptin levels in CSF are unchanged as subjects progress to AD. However, in AD hippocampus, leptin signalling was decreased and leptin localization was shifted, being more abundant in reactive astrocytes and less in neurons. Similar translocation of leptin was found in brains from Tg2576 and apoE4 mice. Moreover, an enhancement of leptin receptors was found in hippocampus of young Tg2576 mice and in primary astrocytes and neurons treated with Aβ₁₋₄₂. In contrast, old Tg2576 mice showed decreased leptin receptors levels. Similar findings to those seen in Tg2576 mice were found in apoE4, but not in apoE3 mice. These results suggest that leptin levels are intact, but leptin signalling is impaired in AD. Thus, Aβ accumulation and apoE4 genotype result in a transient enhancement of leptin signalling that might lead to a leptin resistance state over time.

The primary constituents of plaques (Aβ42/Aβ40) and neurofibrillary tangles (tau and phosphorylated forms of tau [ptau]) are the current leading diagnostic and prognostic cerebrospinal fluid (CSF) biomarkers for AD. In this study, we performed deep sequencing of APP, PSEN1, PSEN2, GRN, APOE and MAPT genes in individuals with extreme CSF Aβ42, tau, or ptau levels. One known pathogenic mutation (PSEN1 p.A426P), four high-risk variants for AD (APOE p.L46P, MAPT p.A152T, PSEN2 p.R62H and p.R71W) and nine novel variants were identified. Surprisingly, a coding variant in PSEN1, p.E318G (rs17125721-G) exhibited a significant association with high CSF tau (p = 9.2 × 10(-4)) and ptau (p = 1.8 × 10(-3)) levels. The association of the p.E318G variant with Aβ deposition was observed in APOE-ε4 allele carriers. Furthermore, we found that in a large case-control series (n = 5,161) individuals who are APOE-ε4 carriers and carry the p.E318G variant are at a risk of developing AD (OR = 10.7, 95% CI = 4.7-24.6) that is similar to APOE-ε4 homozygous (OR = 9.9, 95% CI = 7.2.9-13.6), and double the risk for APOE-ε4 carriers that do not carry p.E318G (OR = 3.9, 95% CI = 3.4-4.4). The p.E318G variant is present in 5.3% (n = 30) of the families from a large clinical series of LOAD families (n = 565) and exhibited a higher frequency in familial LOAD (MAF = 2.5%) than in sporadic LOAD (MAF = 1.6%) (p = 0.02). Additionally, we found that in the presence of at least one APOE-ε4 allele, p.E318G is associated with more Aβ plaques and faster cognitive decline. We demonstrate that the effect of PSEN1, p.E318G on AD susceptibility is largely dependent on an interaction with APOE-ε4 and mediated by an increased burden of Aβ deposition.