Genomic mosaicism with increased amyloid precursor protein (APP) gene copy number in single neurons from sporadic Alzheimer's disease brains

  1. Diane M Bushman
  2. Gwendolyn E Kaeser
  3. Benjamin Siddoway
  4. Jurgen W Westra
  5. Richard R Rivera
  6. Stevens K Rehen
  7. Yun C Yung
  8. Jerold Chun  Is a corresponding author
  1. Dorris Neuroscience Center, The Scripps Research Institute, United States
  2. University of California, San Diego, United States

Abstract

Previous reports have shown that individual neurons of the brain can display somatic genomic mosaicism of unknown function. In this study, we report altered genomic mosaicism in single, sporadic Alzheimer's disease (AD) neurons characterized by increases in DNA content and amyloid precursor protein (APP) gene copy number. AD cortical nuclei displayed large variability with average DNA content increases of ∼8% over non-diseased controls that were unrelated to trisomy 21. Two independent single-cell copy number analyses identified amplifications at the APP locus. The use of single-cell qPCR identified up to 12 copies of APP in sampled neurons. Peptide nucleic acid (PNA) probes targeting APP, combined with super-resolution microscopy detected primarily single fluorescent signals of variable intensity that paralleled single-cell qPCR analyses. These data identify somatic genomic changes in single neurons, affecting known and unknown loci, which are increased in sporadic AD, and further indicate functionality for genomic mosaicism in the CNS.

https://doi.org/10.7554/eLife.05116.001

eLife digest

The instructions for living cells are contained in certain stretches of DNA, called genes, and these instructions have been largely considered to be invariant, such that every cell in the body has the same DNA. However, research has revealed that many neurons in the human brain can contain different amounts of DNA compared to other cells. When cells with varied DNA are present in the same person, it is referred to as mosaicism. The effects of this mosaicism are unknown, although by altering the instructions in brain cells, it is suspected to influence both the normal and diseased brain.

The brains of patients with Alzheimer's disease often contain deposits of proteins called amyloids. The precursor of the protein that makes up most of these deposits is produced from a gene called the amyloid precursor protein gene, or APP. Having an extra copy of the APP gene can cause rare ‘familial’ Alzheimer's disease, wherein the APP duplication can be passed on genetically and is present in all the cells of a patient's body. By contrast, ‘sporadic’ Alzheimer's disease, which constitutes around 95% of cases, does not show any difference in the number of APP genes found in tissue samples, including whole brain. The early studies that discovered this were conducted before an appreciation of brain mosaicism, and thus single neurons were not investigated. This raises the possibility that the number of APP genes may be mosaically increased, which would not be detected by examining non-brain or bulk brain tissue.

Bushman, Kaeser et al. used five different types of experiments to examine the DNA content of single neurons and investigate whether mosaicism could explain the discrepancy in the results of the previous studies. The neurons from people with Alzheimer's disease contained more DNA—on average, hundreds of millions of DNA base pairs more—and more copies of the APP gene, with some neurons containing up to 12 copies.

Bushman, Kaeser et al.'s findings present evidence of a way that mosaicism can affect how the brain works by altering the number of gene copies, and how this impacts the most common form of Alzheimer's disease. Many questions arise from the work, including when does mosaicism arise, and what promotes its formation? How does this relate to age? What parts of the genome are changed, what genes are affected, and how do these changes alter neuronal function?

Furthermore, Bushman, Kaeser et al.'s work suggests that mosaicism may also play a role in other brain diseases, and could also provide new insights into the normal, complex functions of the brain. In the future, this knowledge could help to identify new treatments for brain diseases; for example, by identifying new molecular targets for therapy hidden in the extra DNA or by understanding how to alter mosaicism.

https://doi.org/10.7554/eLife.05116.002

Introduction

The genome has been classically viewed as being constant from cell to cell in the same individual, with genomic differences passed on through the germline. However, within neurons of the brain, numerous studies have reported somatic variability producing complex genomic mosaicism but having unknown function. Identified forms of somatically arising genomic mosaicism include aneuploidy (reviewed in Bushman and Chun, 2013), LINE elements (Muotri et al., 2005; Baillie et al., 2011; Evrony et al., 2012), copy number variations (CNVs) (Gole et al., 2013; McConnell et al., 2013; Cai et al., 2014), and DNA content variation (DCV) (Westra et al., 2010; Fischer et al., 2012).

AD is the most common form of dementia and is characterized by the presence of amyloid plaques, synaptic loss, and cell death (Alzheimer's Association, 2013), notably affecting the prefrontal cortex. The major component of these plaques is β-amyloid (Aβ), a protein encoded by APP (Goldgaber et al., 1987; St George-Hyslop et al., 1987; Tanzi et al., 1987). Familial AD accounts for less than 5% of all cases and has been genetically linked to mutations in APP and two presenilins (PSEN), PSEN1 and PSEN2, the catalytic components of γ-secretase and the units responsible for cleavage of APP (Price and Sisodia, 1998; Bertram et al., 2010). In addition, APP gene dosage is strongly associated with AD pathogenesis based on multiple lines of evidence. First, Down syndrome (DS), with three copies of APP, produces neuropathology virtually identical to AD (Glenner and Wong, 1984; Delabar et al., 1987) and APP locus duplications are sufficient to cause familial AD (Rovelet-Lecrux et al., 2006; Sleegers et al., 2006; McNaughton et al., 2012). Moreover, AD-protective effects have been reported in DS with APP deletion via partial trisomy (Prasher et al., 1998), as well as in familial AD with an APP partial loss-of-function mutation (Jonsson et al., 2012). However, seminal studies in the 1980s failed to detect evidence of APP amplification in sporadic Alzheimer's disease peripheral blood and whole brain (Podlisny et al., 1987; St George-Hyslop et al., 1987; Tanzi et al., 1987) despite strong linkage in familial AD, thus linkage between sporadic AD and APP remains unclear.

The existence of region specific genomic mosaicism in the normal brain (Westra et al., 2010) raised the possibility that DCV, defined as variations in the total DNA amount present in a single cell or population, might play a functional role in sporadic brain diseases by altering pathogenic loci in individual cells. The validated pathogenicity of APP in familial AD suggested that mosaic alterations in APP copy number within single neurons may play a role in producing sporadic AD. Through the use of five independent experimental approaches, we report increased somatic genomic variation within individual sporadic AD neurons involving mosaic increases in both DNA content and APP copy number.

Results

Methodologies to identify genomic mosaicism in AD brain nuclei utilized multiple independent approaches (Figure 1). First, neuronal nuclei were isolated from paired prefrontal cortex and cerebellum of postmortem human brains as previously described (Figure 1A) (Westra et al., 2008, 2010). Nuclei were analyzed for: DNA content changes (Figure 1B), APP CNVs using qPCR in small populations (Figure 1C) and single-nuclei (Figure 1D), trisomy 21 using standard fluorescence in situ hybridization (FISH) (Figure 1E), and amplified APP loci using PNA-FISH (Figure 1F).

Methodologies used in assessing genomically mosaic AD.

(A) Neuronal nuclei were isolated from the prefrontal cortex and cerebellum of postmortem human brain (see ‘Materials and methods’ for samples used) as described (Westra et al., 2008). (B) Nuclear DNA was stained with propidium iodide (PI) and DNA content was quantified using flow cytometric analysis. (C) APP copy number variations were analyzed in small populations of nuclei (∼75 genomes) using custom primers for exon 14 of APP. (D) Single-cell qPCR assessed APP copy number variations in individual neuronal nuclei via TaqMan probes and a modified Biomark integrated fluidic chip system (Fluidigm Corporation, South San Francisco). (E) FISH paints against the whole q arm of chromosome 21 and a point probe against a region on the q arm of 21 (21q22.13-q22.2) were used to double-label and call aneusomies in AD samples. (F) Peptide nucleic acid (PNA) FISH was combined with super resolution microscopy for threshold detection of APP copy number above ∼2 occurring at a single locus.

https://doi.org/10.7554/eLife.05116.003

AD cortices show increased DNA content variation (DCV)

Previous work identified DNA content changes with regional variability in the non-diseased brain, where prefrontal cortical nuclei—particularly neuronal nuclei—displayed increased DCV compared to nuclei of the cerebellum and non-brain controls (Westra et al., 2010). The current study utilizes the same techniques for DNA content analysis by flow cytometry using propidium iodide (PI). PI staining is the predominant methodology for quantitatively differentiating nuclei or cells with variable DNA content and is routinely utilized as a gold standard in multiple fields including genomic comparisons across species in botany (Dolezel et al., 2007), studies of the cell cycle (Krishan, 1975), and DNA degradation produced by apoptosis (Riccardi and Nicoletti, 2006). Prior analyses ruled out the effects of DNA dyes, nuclear size, mitochondrial contamination, and autofluorescent lipofuscin on DNA content and validated genomic increases in DNA content using quantitation of CENP-B PNA probes against human centromere repeats (Westra et al., 2010) which have been further substantiated by identification of copy number gains in single human neurons (Gole et al., 2013; McConnell et al., 2013).

To further validate DCV in human neurons, whole genome amplification (WGA) was used on cortical neuronal nuclei sorted into populations of high or low DNA content based upon PI intensity. Nuclei with high or low DNA content were sorted into 12 replicates of 1000, 500, or 100 and were then subject to DNA content assessment by WGA to assess starting amounts of DNA template in each sample (Figure 2A). Nuclei were denatured and amplified by multiple displacement amplification (MDA) during which DNA synthesis was continually measured by SYBR Green fluorescence (Figure 2B). In every case, nuclei with high PI intensity also showed increased DNA synthesis over those with low PI intensity. These results independently confirm, as expected, that PI staining intensity faithfully reports DNA content.

Figure 2 with 2 supplements see all
AD cortical nuclei show increased DNA content variation (DCV) by flow cytometry.

(A) Histogram displaying gating parameters used in sorting ‘high’ and ‘low’ DNA content populations for validation of DNA content. (B) Validation of DNA content analyses using semi-quantitative MDA whole-genome amplification (WGA) on ‘high’ and ‘low’ DNA content populations of 1000, 500, and 100 nuclei. (C and D) Representative DNA content histograms for lymphocytes (LYM), AD cerebellum (CBL), and AD prefrontal cortex (CTX). Each colored histogram represents a separate sample in each set; CTX and CBL samples are from paired brains. Chicken erythrocyte nuclei (CEN) were used as internal calibration controls. (E) Representative orthogonal view of DNA content vs forward scatter width (FSC-W). For each brain sample, the area to the right of the vertical line indicates a DNA content increase of the population of nuclei. AD-6 CTX is a representative right-hand peak shift and AD-7 is a representative right-hand shoulder (see Figure 3A for more examples). (F) DNA content changes for all human LYM, ND, and AD brain samples examined (AD CTX N = 32, AD CBL N = 16, LYM N = 15 [20 meta analysis], ND CTX = 21 [36 meta analysis], ND CBL = 11 [12 meta analysis]). Red bars denote average for each group relative to lymphocytes. Averages are as follows (including metadata from Westra et al. (2010)): AD CTX 8.219%; AD CBL −0.1104%; LYM −0.2915%; ND CTX 2.239%; ND CBL −3.358%. (G) DNA content changes of the current study (AVOVA p < 0.0001). (H) DNA content changes of the current study combined with metadata from Westra et al. (2010) (ANOVA p < 0.0001). (I) Comparison of mean coefficient of variation (CV statistic from FlowJo of the population, included metadata from Westra et al., 2010) demonstrates that there is an average increase in the variation of AD samples (ANOVA p < 0.0001). *p = 0.05, **p = 0.01, ***p = 0.001, ****p < 0.0001, See Figure 2—source data 1 for exact p values. See Figure 2—figure supplement 1 for age, PMI and Braak score correlations. See Figure 2—figure supplement 2 for control of nuclear size analysis.

https://doi.org/10.7554/eLife.05116.004
Figure 2—source data 1

DNA Index (DI) and percent change values and statistics.

https://doi.org/10.7554/eLife.05116.005

AD neuropathology strongly affects the prefrontal cortex. We therefore first interrogated the DNA content of pathologically confirmed prefrontal cortices (N = 32), using previously described methodologies for DCV analyses in neuronal nuclei (Figure 1B) (Westra et al., 2010). In control experiments, AD cerebellar nuclei showed DNA content profiles similar to lymphocytes, characterized by histograms with sharp peaks and narrow bases (Figure 2C,E). By comparison, AD cortices displayed high variability characterized by right hand shoulders (Figure 2D,E [AD-7]), large right hand peak shifts (Figure 2D,E [AD-6]), and wide bases. Of the AD cortex samples examined, greater than 90% displayed a net DNA content increase, averaging approximately 8% gain over human lymphocyte controls (Figure 2F–H). Notably, the AD cortex displayed increases beyond those observed in non-diseased cortices, with an average gain of approximately 6% over age and sex-matched samples (Figure 2F–H). AD cortices also displayed an increased coefficient of variation over non-diseased cortical nuclei and all cerebellar nuclei (Figure 2I) as well as a consistent skewed distribution compared to AD cerebellum (Figure 2—figure supplement 1A). In addition, we examined 14 paired sets of AD cortex and cerebellum from the same individual (Figure 3A) and 12 paired sets from non-diseased individuals (Figure 3B). Each AD cortex showed unique cortical histograms and increases in total cortical DNA compared to the cerebellum. In AD and non-diseased brain samples, DNA content changes did not correlate with age, Braak score, or postmortem index (Figure 2—figure supplement 1B–H), and DNA content was independent of nuclear size (Figure 2—figure supplement 2) (Westra et al., 2010). Importantly, the increase observed was significantly less than a 4N tetraploid genome (∼12,800 Mb), yet significantly more than what would be expected from a hypersomy of even the largest chromosome (chr1: ∼250 Mb or 3.9%). These results indicate that neuroanatomically selective increases in DNA content represent a distinct, reproducible, and prevalent characteristic of the AD brain.

Pairwise DNA content analyses in AD cortical nuclei vs AD cerebellum.

(A) Pairwise analysis of overlaid DNA content histograms (CTX = solid red, CBL = black dashed lines) in the same AD individual (each graph represents a unique AD individual). (B) Pairwise analysis of overlaid DNA content histograms (CTX = solid blue, CBL = black dashed lines) in the same ND individual.

https://doi.org/10.7554/eLife.05116.008

Neurons produce DNA content increases

In the non-diseased cortex, neurons were identified as the predominant cell type contributing to increased DCV (Westra et al., 2010). To determine whether neurons were also responsible for increased DCV in sporadic AD, nuclei were immunolabeled for the neuronal nuclear antigen, NeuN, and flow cytometry was used to analyze DCV in NeuN-positive and NeuN-negative nuclei (Figure 4A–C). NeuN-positive AD cortical neurons showed right hand DNA content shifts (Figure 4D) and increased DCV, displaying ∼9% gain over that observed in NeuN-negative nuclei (Figure 4E). A comparison of NeuN-positive nuclei in paired samples from the same brain also showed significantly increased DNA content in cortical neurons over cerebellar neurons (Figure 4F). These distinctions do not rule out AD-specific effects on DCV from non-neuronal cells, but do implicate neurons as the major cellular locus for increased DCV.

DNA content increases observed in AD cortical nuclei are attributable to neurons.

(AC) The gating procedure used for NeuN-positive flow cytometry analysis. (A) DNA content peak for identified nuclei. (B) A sample of unlabeled neuronal nuclei that display no NeuN-positive signal. (C) Selection for NeuN-positive nuclei for downstream analysis. (D) DNA content histograms of four AD samples displaying NeuN-positive nuclei (solid purple) vs NeuN-negative nuclei (black dashed line). NeuN-positive populations display distinct cortical histograms with prominent right-shifted peaks (arrows). (E) Comparison of DNA index (DI) increases from NeuN-positive nuclei (solid purple) vs NeuN-negative nuclei (black dashes) from AD CTX samples. NeuN-positive nuclei (DI = 1.10) showed an average gain of 9% over NeuN-negative nuclei (DI = 1.01), **p = 0.0011. (F) Comparison of DNA content in NeuN-positive nuclei from AD CTX (DI = 1.10) (red) vs AD CBL (DI = 0.94) (pink) from the same individual; CTX nuclei displayed a 15.6% gain over CBL nuclei, *p = 0.0335. Statistics are paired two-tailed t-test. Bars indicate ± SEM.

https://doi.org/10.7554/eLife.05116.009

APP copy number variably increases in small cell cohorts

The DNA content increases in sporadic AD neurons raised the question of what genomic loci may be specifically increased in AD. The most highly validated gene associated with AD is APP, wherein constitutive gains have previously been linked to AD through familial AD and DS (Prasher et al., 1998; Rovelet-Lecrux et al., 2006; Jonsson et al., 2012; McNaughton et al., 2012). However, prior large-scale analyses in peripheral, non-neural tissues from sporadic AD cases have not reported APP copy number changes (Bertram and Tanzi, 2009; Harold et al., 2009). The normal existence of genomic mosaicism might produce ‘hot-spots’ of cells with increased APP copy number, and we therefore interrogated small cohorts of brain nuclei from paired cortex and cerebellum. Three of six AD samples displayed significant increases in APP copy number (3.9–6.3 copies compared to reference genes) over paired cerebellar samples (Figure 5A (AD-1, AD-3, AD-6)). Similar APP amplification was not detected in non-diseased brains (Figure 5B), and DS controls along with reference gene assessments demonstrated expected results (Figure 5—figure supplement 1A,B). Notably, AD samples displayed substantial variability in APP copies (Figure 5C), revealing inherent limitations in precisely quantifying copy numbers in small, genomically mosaic cell populations and highlighting the necessity for single-cell analyses.

Figure 5 with 1 supplement see all
Mosaic amplification of the APP locus in small cohorts of AD cortical neurons unrelated to trisomy 21.

(A) Comparison of relative copy number of APP in CBL and CTX fractions from six AD brains. APP locus-specific amplification was determined relative to reference gene SEMA4A; paired CBL nuclei were used as a calibrator sample for each brain, normalized to 2.00 for a diploid cell. Differences in ΔΔCt ± SEM of APP in the cortex vs cerebellum were assessed in each individual using an unpaired, two-tailed t-test (****p = 0.0001, *p = 0.0165, *p = 0.0489) (B) Comparison of relative copy number of APP in CBL and CTX fractions from 4 non-diseased brains. (C) Average relative copy number in non-diseased vs AD brains. Control genes and DS individuals were also examined (Figure 5—figure supplement 1). (DJ) FISH strategy of chromosome 21 counting through simultaneous labeling using chr 21 q arm ‘whole’ chromosome paint (WCP, green) and chr 21 regional FISH probe for 21q22.13-q22.2 (red) (see Figure 5—source data 1 for raw counts). (D and E) The ability to detect aneuploidy was validated using interphase nuclei from a human trisomy 21 brain, where three regional spots (red, encompassing the APP gene) were seen, despite WCP spatial variation (see also Rehen et al. (2005)). (FI) Chromosome 21 aneusomy was examined in prefrontal cortical nuclei. Examples of chr 21 (F) monosomy, (G) disomy, (H) trisomy, and (I) tetrasomy (please note tetrasomy is not an example of aneuploidy). (J) Quantification of individual FISH signals showed no significant differences in monosomy, disomy, trisomy, or tetrasomy. 5 control brains and 9 AD brains were used. At least 450 nuclei were quantified per brain sample. Scale bar = 10 um. 4974 total nuclei examined.

https://doi.org/10.7554/eLife.05116.010
Figure 5—source data 1

Raw dual point-paint probe FISH counts.

https://doi.org/10.7554/eLife.05116.011

APP copy number increases in AD are not due to trisomy 21

APP gene dosage in familial AD and DS has driven hypotheses connecting AD pathogenesis with increased incidence of trisomy 21 (Heston and Mastri, 1977; Potter, 1991; Geller and Potter, 1999). We therefore examined AD cortical nuclei using a highly liberal protocol for calling aneusomies whereby borderline FISH profiles suggestive of aneusomy were always included in quantitative assessments, in an effort to detect possible differences between AD and control brains. Importantly, all analyses were conducted blind to the identity of samples, an approach made possible by interrogating purified nuclei rather than cells or tissues that themselves show identifying increases of plaques and tangles in AD. Nuclei were double-labeled with a commercially available chromosome 21 q-arm ‘whole chromosome paint’ (WCP) and a regional ‘point’ probe for 21q22.13-q22.2 (220 Kb) (Figure 5D–I). The ability of this technique to detect aneuploidy was validated using interphase nuclei from a human trisomy 21 brain (Svendsen et al., 1998) revealing three nuclear signals (Figure 5D,E) (see also (Rehen et al., 2005)). Three separate, blinded observers counted each sample. Despite using liberal counting criteria, no statistically significant changes in chromosome 21 aneuploidy rates, including trisomies, between AD (N = 4974 nuclei, N = 9 brains) and non-diseased brains (N = 2576 nuclei, N = 5 brains) were observed (Figure 5J). Comparably high levels of aneuploidy have been reported for chromosome 21 displaying no difference between AD and non-diseased hippocampal cells (Thomas and Fenech, 2008). Thus, results from dual probe FISH analyses do not support increased trisomy 21 in AD but are consistent with alternative CNV mechanisms that could produce APP copy number values exceeding 3, as was observed in 75-genome qPCR analyses (Figure 5).

APP copy number is increased and enriched in single AD neurons

Mosaic single-cell increases in APP copy number could explain variations observed in small cohort qPCR data. To assess APP copy number in single neurons, an optimized microfluidic protocol to assess genomic copy number via qPCR using a Biomark HD 48.48 Dynamic Array integrated fluidic circuit (IFC) (Fluidigm, South San Francisco) was adapted from gene expression protocols. Nuclei were isolated and interrogated at two reference genes (see Materials and methods) and two APP exons that flanked the majority of the coding sequence at the 5′ and 3′ ends of APP (exons 3 and 14). Samples were run in triplicate and assays in sextuplicate, which generated high numbers of replicate data points allowing for improved copy number resolution to an extent not possible by conventional qPCR (Weaver et al., 2010; Whale et al., 2012). A total of 154 neuronal nuclei were individually examined from three AD and three non-diseased brains, within which nuclei from the cerebellum and prefrontal cortex were separately analyzed. In single AD cortical nuclei, significant increases in APP copy number, ranging up to 12 copies in a single nucleus (Figure 6A,B), were observed with a high concordance rate between exons (Figure 6—figure supplement 1). AD cortical nuclei displayed increased average copy numbers (∼3.8–4 copies) over control samples (∼1.7–2.2 copies) (Figure 6A,B), with increased frequencies of high copy number nuclei (six or more copies) (Figures 6C and 5E, red lines) primarily occurring in prefrontal cortex samples (Figure 6C–F). The AD cortex showed an approximately fourfold increase over the non-diseased cortex of nuclei with greater than 2 APP copies (55–66% vs 12–15%) (Figure 6G,I). In addition, nuclei with fewer than 2 APP copies were also identified. In neuronal nuclei with greater than two copies (Figure 6G,I, gold bars), AD cortical nuclei showed statistically significant increases (∼5 copies for both exons) over the other samples (∼3 copies for both exons) (Figure 6H,J). These data identify mosaic, neuroanatomically enriched and disease-associated increases in APP copy number in single, sporadic AD neurons.

Figure 6 with 1 supplement see all
Mosaic APP locus amplification in single neurons from AD brains.

(A) Single nuclei relative copy numbers for exon 3 of APP from non-diseased (ND) CBL, ND CTX, AD CBL, and AD CTX; each black diamond represents one neuron. For each group, the mean is displayed in red and bars represent 95% confidence intervals. AD CTX showed a mean APP copy number of 3.80; this is significantly higher than AD CBL (2.23), ND CTX (1.60), and ND CBL (2.28). *p = 0.0147, **p = 0.0015, **p < 0.0012, ANOVA p < 0.0001 (see Figure 6—source data 1 for raw numbers and statistics). (B) Single nuclei relative copy numbers for exon 14 of APP, similar to (A). The two exons showed a high concordance (Figure 6—figure supplement 1) where the AD CTX showed a mean APP copy number of 3.40 while the AD CBL (2.34), ND CTX (1.44), and ND CBL (1.92) remained closer to 2 copies. *p = 0.0163, **p = 0.0016, ****p < 0.0001, ANOVA p < 0.0001. (CF) Distribution of copy number calls for exon 3 (C and D) and exon 14 (E and F) binned by relative copy number. The AD CTX for both exons displayed unique distributions, with more nuclei falling into the high copy number bins. (G and I) Distribution of nuclei with copy numbers less than, equal to, and greater than two copies. (H and J) Average copy number increases in nuclei binned with greater than two copies (gold columns in G) (AD CTX: Exon 3 = 5.01, Exon 14 = 4.96, *p = 0.0361). All statistics represent an ANOVA with a Tukey's multiple comparison test. Bars indicate ± SEM.

https://doi.org/10.7554/eLife.05116.013
Figure 6—source data 1

Single Cell qPCR Data and Statistics.

https://doi.org/10.7554/eLife.05116.014

APP amplified loci can be visualized by PNA-FISH

The identification of high APP copy number nuclei by single-cell qPCR suggested the possibility that amplified loci could produce detectable FISH signals. In chromosome 21 non-quantitative regional point probe FISH experiments, we occasionally observed variable puncta sizes (Figure 5H,I) (Rehen et al., 2005) that were initially dismissed as technical hybridization variability but which might potentially indicate sub-chromosomal variation of loci containing target sequences. However, conventional point-probe FISH analysis, while useful for assessments of aneusomy or rearrangements, cannot quantitatively evaluate copy number changes occurring as contiguous copies in close proximity to one another. Moreover, other conventional techniques for identifying copy number in single cells, such as single cell sequencing, rely heavily on genome amplification, which may introduce bias or mask CNVs. We therefore developed a detection assay based upon PNA chemistry. PNA probes hybridize with single base discrimination and have been used to quantify short repeat sequences such as those found on telomeres (Buchardt et al., 1993; Lansdorp et al., 1996). This raised the possibility that amplified copies of APP could be identified using multiple PNA probes simultaneously hybridized to unique gene sequences, provided that amplification occurred in a physically constrained locus rather than being dispersed throughout the genome.

Nine PNA probes of 12–18 residues, each conjugated to a single Alexa-488 fluor, were designed against multiple sites of the same APP exons examined by single-cell qPCR, and validated both in silico and by blotting for specificity and linear quantitative behaviors (Figure 7A, Figure 7—figure supplement 1A–C). Relative probe sequence locations were also designed to avoid fluorochrome proximity quenching. Single PNA probes did not produce a detectable signal in any samples using standard fluorescence microscopy used for aneusomy FISH (Rehen et al., 2001, 2005), contrasting with clearly detectable telomeric signals, which led to evaluation of more sensitive microscopic techniques including confocal, deconvolution (Westra et al., 2010), and ultimately super resolution structured illumination microscopy (N-SIM, Nikon) (Gustafsson, 2000), which all failed to detect possible signals in AD nuclei. However, hybridization of increasing numbers of distinct PNA probes combined with N-SIM visualization resulted in the empirical determination of a threshold that identified one, and rarely two, punctate signals that were much more frequent in AD nuclei (Figure 7B), and with a frequency consistent with single-cell qPCR data (Figure 6). AD neuronal nuclei showed single purely green puncta (Figure 7B, #1 arrow, Video 1–6) and rarer two puncta signals (Figure 7B, #3 arrow, Video 6), all of which could be differentiated from lipofuscin puncta that fluoresced in all channels (Figure 7B, #2 arrow, Video 5). Internal positive control signals identified telomeres (Figure 7, red fluorescent puncta). Use of 3D-SIM enabled acquisition and analysis of super resolution projections of individual green puncta (Figure 7B) revealing a range of morphologies and intensities (Figure 7B,D, Figure 7—figure supplement 1D,E) consistent with varied APP copy number and possibly distinct, intrachromosomal genomic organization. This threshold detection using multiple PNA probes targeting APP visualized green puncta in 56% of AD neuronal nuclei compared to 22% in non-diseased and 14% in DS (Figure 7C). There was no evidence of single-copy detection in any sample, contrasting with standard FISH data (Figure 5D–I), despite always showing telomeric signals. Compared to single-cell qPCR data, which reports both endogenous alleles in addition to amplification, copy numbers at a single locus below ∼2 were not visualized with PNA-FISH, allowing focused analyses on positive profiles resulting in fluorescence intensity plots ( Figure 7D and Figure 7—figure supplement 1D) that were highly reminiscent of APP copy number plots produced by single-cell qPCR (Figure 6).

Figure 7 with 1 supplement see all
Visualization of APP copy number increases in neuronal nuclei from AD brain samples.

(A) Peptide nucleic acid probes (PNA) were developed against nine separate sites on APP (4 sites within exon 3 and 5 sites within exon 14). Each PNA probe consists of a peptide backbone conjugated to a single fluorophore, with separately conjugated nucleotides, substantially increasing specificity (Lansdorp et al., 1996). Single copies of APP are not detectable because of fluorophore detection limits. Detection of increased copy number by PNA probes can be visualized as copies of APP increase (Figure 7—figure supplement 1B,C). Positive internal controls using PNA probes directed against telomere sequences were simultaneously hybridized. (B) Visualization of copy number increases in neuronal nuclei. Green puncta (arrow 1, insets) indicate visualized APP increases. Telomere labeling (red puncta) was present in all nuclei, demonstrating probe accessibility and template fidelity. Lipofuscin (arrow 2, orange puncta) was detected in nuclei, visualized by extensive fluorescence signal in all channels, but was eliminated from quantifications. Limited nuclei displayed two green puncta (arrow 3). V1-6 Refers to the supplemental videos where 3-D projections can be visualized. (C) Graphic representation of non-diseased (blue) and DS (grey) brains displayed limited numbers of threshold-detected increases in APP (Figure 7—source data 1). AD (red) brains displayed significant and consistent threshold-detected increases in APP. (D) Individual threshold-detected APP increases were quantified and plotted on a relative intensity scale (blue diamonds: non-diseased, red diamonds: AD). Dotted line represents the threshold below which APP copy number was undetectable, only limited puncta were identified in non-diseased nuclei. Bars indicate ± SEM, *p < 0.05.

https://doi.org/10.7554/eLife.05116.016
Figure 7—source data 1

Data and statistics for PNA-FISH counts.

https://doi.org/10.7554/eLife.05116.017
Video 1
PNA-Fish analysis of APP in nuclei from non-diseased cortical neurons.

Video of 3-D projection from Figure 7B, V1. Green puncta indicating APP increases were infrequently visualized in non-diseased brains. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.019
Video 2
PNA-Fish analysis of APP in nuclei from non-diseased cortical neurons with lipofuscin.

Video of 3-D projection from Figure 7B, V2. Green puncta indicating APP increases were infrequently visualized in non-diseased brains. Lipofuscin (orange puncta), visualized by extensive fluorescence signal in all channels, was detected in some nuclei, but was excluded from analysis. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.020
Video 3
PNA-Fish analysis of APP in nuclei from non-diseased cortical neurons.

Video of 3-D projection from Figure 7B, V3. Green puncta indicating APP increases were infrequently visualized in non-diseased brains. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.021
Video 4
PNA-Fish analysis of APP in nuclei from AD cortical neurons.

Video of 3-D projection from Figure 7B, V4. Green puncta indicate visualized APP increases. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.022
Video 5
PNA-Fish analysis of APP in nuclei from AD cortical neurons with lipofuscin.

Video of 3-D projection from Figure 7B, V5. Green puncta indicate visualized APP increases. Lipofuscin (orange puncta), visualized by extensive fluorescence signal in all channels, was detected in some nuclei, but was excluded from analysis. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.023
Video 6
PNA-Fish analysis of APP in nuclei from AD cortical neurons.

Video of 3-D projection from Figure 7B, V6. Green puncta indicate visualized APP increases. Limited nuclei in AD displayed two green puncta. Red puncta indicate telomere labeling with separate telomere-specific PNA probes and were visualized as a positive control for PNA hybridization.

https://doi.org/10.7554/eLife.05116.024

Discussion

Our results support a mechanism for the development of sporadic AD whereby the validated familial pathogenic APP gene is mosaically amplified in neurons. Mosaicism was documented based on multiple independent criteria: increased DCV in 90% of examined AD brains, varied and increased APP copy number in small cohorts, increased APP CNVs in single neurons identified by single-cell qPCR, and APP copy number increases in single nuclei visualized by PNA-FISH. These data are consistent with our inability to substantiate increases in trisomy 21 in cells of the sporadic AD brain. The concomitant increase in DCV suggests effects on other loci that are also altered in AD, the most obvious of which would be copy gains, but which could also involve losses based upon the presence of increased DCV range observed in AD; such loci might contribute to the progeric presentation of sporadic AD and possibly familial forms of AD that still require decades for disease to manifest. Our working hypothesis ties the currently accepted pathogenicity of APP gene dosage, established by familial AD and DS, to sporadic AD through a mechanism of somatic, mosaically increased APP copy number in some neurons. Our data do not exclude the possibility that the changes observed are a downstream effect of causative factors in the disease. However, if mosaic genomic changes are downstream of disease onset, APP copy number changes would likely play a significant role in disease progression.

Advances in single-cell genomic sequencing had initially suggested its use for identifying APP CNVs in single AD neurons. However, existent technology is limited to published sequence resolution between ∼2-5 Mb and ∼0.025X genome coverage (Evrony et al., 2012; Gole et al., 2013; McConnell et al., 2013; Cai et al., 2014). Considering the <0.3 Mb size of the APP locus, single cell sequencing would be incapable of identifying all but the rarest APP CNVs observed here. In addition, whole genome single cell sequencing is limited by throughput, demonstrated by the low number of neurons reported in this field, and is complicated further by variable results (Cai et al. N = 82 (QC neurons), Evrony et al. N = 6, Gole et al. N = 6, McConnell et al. N = 110). Notably, DNA losses predominated in both Cai et al. (2014) and McConnell et al. (2013) however, Gole et al. (2013) observed that two thirds of somatic CNVs were gains, consistent with increases in DNA content reported previously, and all three studies support a range of DCV amongst neurons, although none of these prior studies assessed AD neurons. While future advances will improve sequence resolution, throughput, and amplification fidelity in single cells, the distinct single cell strategy employed here allowed copy number assessment of a single targeted gene, APP, an approach which may be generalizable to other loci and diseases. Compared to the prior single-neuron reports, 320 single, neuronal nuclei were assessed here for APP CNVs wherein single cell qPCR data markedly resembled PNA-FISH data. Notably, these single cell techniques also possess limitations. For example, single cell qPCR requires normalization to single copy reference genes and control populations, which presents unique difficulties in assessing mosaic neuronal populations. For example, changes in copy number could reflect changes in both the reference and target gene, which might result in artifactual increases (or decreases) in a target gene. Similarly, assessment of copy number from an amplified template could also produce artifactual changes, and thus results require independent techniques for verification, especially since the original template is consumed in the single cell reaction and cannot be further assessed. The use of PNA-FISH does not require any amplification or normalization, and yet PNA-FISH and single-cell qPCR produced a highly similar distribution of APP copy number amplifications. The employed techniques were not capable of assessing intact neurons or their histological organization, which will require technological modifications. Similarly, assessments of genomic and expression data in a single neuron also await technological advancements. Finally, these techniques are currently capable of interrogating only parts of the entire APP locus; therefore, at this time the precise structure of APP CNVs and the mechanisms leading to their existence are unknown.

Our data have bearing on at least 3 AD hypotheses. First, the prevailing amyloid hypothesis in AD posits that Aβ deposition drives AD (Hardy and Selkoe, 2002). Increased incidences of amyloid senile plaques are observed in all forms of AD and appear to be directly linked to APP dosage in DS and familial AD. The mosaic increases in neuronal APP copy number reported here provide an explanation for the universal presence of Aβ senile plaques in sporadic forms—indeed, all forms—of AD despite an absence of constitutive copy number gain. The presence of rare neurons with APP amplifications in non-diseased brains also provides an explanation for the observed presence of senile plaques in otherwise normal, aged brains (Gibson, 1983; Cras et al., 1991; Mackenzie et al., 1996), consistent with the ages of non-diseased brains examined here, which exceeded 75 years. Along with the occurrence and augmentation of DCV, APP amplification in AD is consistent with dysregulation of normally occurring processes in the etiology of AD pathogenesis. A second hypothesis proposes increased APP dosage through trisomy 21 based largely upon the neuropathology of DS (Heston and Mastri, 1977; Potter, 1991; Geller and Potter, 1999). However, dual chromosome point-paint FISH for chromosome 21 using liberal calling criteria on ∼5000 AD brain cells (N = 14 brains) showed no relationship between total mosaic aneusomies and AD. Moreover, the lack of observed increases of trisomy 21 in sporadic AD are consistent with prior studies in both peripheral non-brain and intact brain tissues, which reported APP levels approximating 2N (Podlisny et al., 1987; St George-Hyslop et al., 1987; Tanzi et al., 1987; Bertram et al., 2010), and is further consistent with mosaic CNVs observed here. The purposeful use of liberal counting criteria resulted in generally higher percentages of aneusomies than reported previously (Rehen et al., 2005; Iourov et al., 2009), however the increases were observed in both AD and non-diseased samples without linkage to AD. Our results do not eliminate roles for mosaic aneuploidy in AD, but do not support trisomy 21 as a specific mechanism for increasing APP copy number in sporadic AD. A third hypothesis is that abnormal cell cycle reentry in post-mitotic neurons contributes to AD (Yang et al., 2001; Herrup and Arendt, 2002; Kruman II et al., 2004; Copani et al., 2006; Mosch et al., 2007; Herrup, 2012). Our DCV analyses detected maximum gains of ∼21% and average gains of ∼8.2% over lymphocytes, representing a fraction of the 100% increase expected in a 4N cell. The subgenomic increases observed here are also consistent with the reported absence of adult neurogenesis in the normal cerebral cortex (Rakic, 2002; Bhardwaj et al., 2006), and could also be relevant to reports of nucleotide incorporation in adult cortical neurons (Gould et al., 1999). Concepts involving changes in DNA synthesis in AD without cell-cycle progression or neurogenesis could be rectified through subgenomic DNA synthesis in post-mitotic neurons that might involve DNA synthetic proteins reported in AD (Copani et al., 2002, 2006).

The observed genomic alterations in AD could occur through both developmental as well as aging processes as reported for at least one form of mosaicism, aneuploidies (Bushman and Chun, 2013) that have been linked to caspase-mediated cell death (Peterson et al., 2012) and could have relevance to developmental APP functions (Nikolaev et al., 2009). Perhaps tellingly, these developmental cell death processes also involve DNA fragmentation and double strand breaks (Blaschke et al., 1996; Staley et al., 1997; Blaschke et al., 1998) that have recently been reported in neurons exposed to physiological and AD stimuli (Suberbielle et al., 2013). It is possible that somatic mosaicism may arise through dysfunction of normally operative DNA repair mechanisms that produce DCV in the normal brain, with augmentation in the diseased brain. We speculate that increased DCV and APP CNVs observed in AD may in part reflect neurons that have delayed or averted cell death, wherein somatic genomic changes could provide a survival advantage at a cost of altered neurophysiological functions.

The DNA content changes observed here showed not only disease relationships but also varied with brain region and cell type. It is important to note that the interrogated neurons are likely terminal in their fate: they will not divide further, contrasting with lymphocytes that, perhaps surprisingly, showed the least amount of DCV despite being a proliferating population. Unlike neurons, proliferating lymphocytes may be similar to stem and progenitor populations that can maintain germline genomes to promote population expansion. By contrast, neurons are by definition post-mitotic, and therefore not subject to this possible requirement. We think it likely that many dividing cells that deviate significantly from the germline genome—and thus would be incapable of further division—are eliminated by cell death and would be removed from the population interrogated for DNA content. Support for this view includes reports on aneuploidies in neuroprogenitor cells, whereby the most extremely aneuploid forms appear to be eliminated by caspase-mediated cell death (Blaschke et al., 1996, 1998; Peterson et al., 2011; Peterson et al., 2012). DNA content changes may therefore be contextual, enabling diverse functions in different cell types. The large differences observed even between brain regions and amongst neuronal nuclei would suggest the existence of distinct neuronal functions produced by DCV, and possibly specific processes or stimuli that can genomically alter neurons.

Diverse stimuli previously implicated as AD risk factors including age and trauma (Sparks et al., 1990; Plassman et al., 2000), might share common endpoints via influences on genomic mosaicism. The consequences of somatic APP amplification in AD may support functionality of other somatically altered genomic loci observed in single neurons that could contribute to the progeric presentation of AD as well as aspects of the disease itself. We further hypothesize that other sporadic or idiopathic brain diseases could arise through altered genomic mosaicism that includes somatic variations at both known and unknown pathogenic loci.

Materials and methods

Human tissue samples

Request a detailed protocol

All human tissue protocols were approved by the Scripps Office for the Protection of Research Subjects (SOPRS) at The Scripps Research Institute (TSRI) and conform to National Institutes of Health guidelines. Fresh-frozen brain tissue was provided by the NICHD Brain and Tissue Bank for Developmental Disorders at the University of Maryland, the University of California Alzheimer's Disease Research Center (UCI-ADRC), and the Institute for Memory Impairments and Neurological Disorders, the Johns Hopkins School of Medicine Alzheimer's Disease Research Center, and Dr Edward Koo at the University of California, San Diego. Lymphocytes from human peripheral blood were obtained from healthy donors at TSRI's Normal Blood Donor Services. Detailed information about the samples used can be found in the table below (Table 1).

Table 1

Human samples used in each experiment

https://doi.org/10.7554/eLife.05116.025
SexAgePaper codeSampleExperimentsPost mortem intervalBraak score
Alzheimer's disease prefrontal cortex, N = 32F811521D23VI
F831562D6U
F74AD-131866DTUU
F79AD-51868DSTUU
F82AD-101875DSTUU
F831893D9U
F871899D5VI
F62AD-91912DSTUU
F80AD-81913DTUU
F54AD-41916DTBUU
F72AD-21921DTBUU
F772400D3.7V
F80AD-112500DP2.3VI
F10150341D19V
F9161788D11V
F9862405D11V
F8962439D11V
F7762509D18VI
M88AD-1102DSBP3IV
M90268D91V
M83736D27U
M821211D18U
M79AD-151252D9U
M921748D5.5V
M85AD-121861DTUU
M85AD-141870DTUU
M82AD-32401DSBP3VI
M84AD-62499DS3.4VI
M63AD-74199DP3VI
M8013173D22IV
M9430022D20V
M9160987D22.5V
Mean82.1
Alzheimer's disease cerebellum, N = 15F74AD-131866DUU
F79AD-51868DSUU
F82AD-101875DSUU
F80AD-81913DUU
F62AD-91912DSUU
F54AD-41916DBUU
F72AD-21921DBUU
M88AD-1102DSB3IV
M79AD-151252D9U
M701625D1U
M85AD-121861DUU
M85AD-141870DUU
M82AD-32401DSB3VI
M84AD-62499DS3.4VI
M63AD-74199D3VI
Mean75.9
Non-diseased prefrontal cortex, N = 40F74ND-21901DSBP2.3II
F74ND-8299D2.8II
F84ND-3703S5.8III
F53ND-111379D15III
F73713*TUU
F9560831D9II
F511568P22U
F171230P16U
F87ND-101502D5II
F8060728D13II
M79ND-9827*DTUU
M96ND-11102DSB3.4II
M83ND-52501DB1.7II
M95ND-41301DS3.5I
M871471*TUU
F711571*TUU
M531344*DTUU
F93318D2.3VI
F92955D20.5III
F564238D12
M704534D28
F9111488D16II
F7913188D12.5
F9013204D9.5II
F10360329D5III
F8560428D8.5III
F9960524D15II
F9562043D20.5
M71389M15
F83719M17III
M69946M12
M872039M6.3III
F864546M22
M9160772M16II
F8061218M5.5
M8761334M8II
M88PDC2MU
M80PDC5MU
M75PDC8MU
Mean79.54
Non-diseased cebellum, N = 15F74ND-21901DSB2.3II
F74ND-8299D2.8II
F84ND-3703S5.8III
F771569D8III
F83719DUU
F531379D15III
F711571DTUU
F87ND-101502D5II
M531344DTUU
M96ND-11102DSB3.4II
M83ND-52501DB1.7II
M95ND-41301DS3.5I
M79ND-9827DUU
M871471DUU
F864546M22U
Mean78.8
DS/AD, N = 3F51DS-1M1864B19U
F47DS-2M3233S24U
F44DS-31258S13U
Mean47.3
LYM, N = 21F405162DN/AN/A
F403963DN/AN/A
F634984DN/AN/A
F604519DN/AN/A
M55Lym 1DN/AN/A
M354651DN/AN/A
29DN/AN/A
187DN/AN/A
4781DN/AN/A
4801DN/AN/A
4903DN/AN/A
M285259DN/AN/A
M5683M
F521344M
F564603M
M544609M
M58Lym 2M
F51Lym 3M
M52Lym 4M
Mean50.0M
  1. Samples in bold are paired CBL and CTX.

  2. *

    Denotes mid frontal gyrus (MFG), D = DNA content analyses, S= Small population qPCR, T = FISH Analysis, B = single cell qPCR on Biomark HD, P=PNA FISH, M = Westra et al. DNA content metadata.

Flow cytometry (FCM) and fluorescence-activated cell sorting (FACS)

Request a detailed protocol

Human brain nuclei were isolated and prepared for FCM and FACS as previously described (Westra et al., 2008, 2010). Isolated nuclei were fixed with 2% paraformaldehyde (or 70% ethanol for microfluidic qPCR), labeled with mouse anti-NeuN antibody (1:100) (Millipore, Germany) and Alexa Fluor 488 goat anti-mouse IgG secondary (1:250) (Life Technologies, Carlsbad, CA), and counterstained with propidium iodide (50 μg/ml) (Sigma, St. Louis, MO) in solution containing 50 μg/ml RNase A (Qiagen, Valencia, CA) and chicken erythrocyte nuclei (CEN) (Biosure, Grass Valley, CA). Electronically gated diploid neuronal nuclei, determined by PI fluorescence and immunolabeling, were analyzed and sorted either in bulk for standard qPCR and PNA-FISH or singly in 96-well plates for microfluidics-based qPCR. FCM and FACS were performed at the TSRI Flow Cytometry Core using a Becton Dickinson (BD Biosciences, San Jose, CA) LSRII and FACS-Aria II, respectively. Post hoc DNA content analyses were performed using FlowJo software (TreeStar Inc., Ashland, OR).

DNA content assessment by whole genome amplification

Request a detailed protocol

For validation of DNA content analyses, nuclei were sorted by FACS in populations of 1000, 500, and 100 into 96-well plates according to relative DNA content (Figure 2A). Nuclei were denatured by multiple freeze–thaw cycles and potassium hydroxide, the solution was neutralized, and phi29 DNA polymerase (Illumina, San Diego, CA), dNTPs, SYBR green, and random hexamers (IDT, Coralville, IA) were added (Gole et al., 2013). The MDA reaction was performed at 30°C and SYBR Green intensity was recorded every 2 min on a QuantStudio RT-PCR 12K Flex (Life Technologies, Carlsbad, CA). 12 replicates were run and averaged for each group.

Quantitative PCR primers

Request a detailed protocol

All qPCR was performed according to MIQE guidelines (Bustin et al., 2009). For standard qPCR, primers against the APP gene exon 14 and reference genes SEMA4A, CCL18, and PCDH11X were designed in house and optimized to an annealing temperature of 59°C using Primer3 software (University of Massachusetts Medical School) and synthesized by Valuegene (San Diego, CA) (Table 2). For microfluidics-based qPCR, TaqMan assays against APP exons 3 and 14 (178 kb apart) (designed from NCBI Reference Sequence: NG_007376.1), SEMA4A and PCDH11X were synthesized by Applied Biosystems (Life Technologies, Carlsbad, CA) or Integrated DNA Technologies (APP exon 3) (San Diego, CA) and optimized to an annealing temperature of 60°C. Primers were assessed in silico to determine specificity of the primer set for a single genomic region and for the presence of SNPs in the targeted genomic region which could decrease primer binding and amplification efficiency. The specificity of all qPCR assays was assessed by gel electrophoresis to confirm a single PCR product of the expected length and sequenced to confirm amplification of the expected product.

Table 2

Primers used for small population and single nuclei qPCR

https://doi.org/10.7554/eLife.05116.026
GeneProteinLocusAssay typePrimer sequenceProbeProduct lengthEfficiency
APP, Exon 14Amyloid precursor protein21q21.3SYBR GreenF-TGCACGTGAAAGCAGTTGAAG, R-AAAGATGGCATGAGAGCATCGN/A2140.973
SEMA4ASemaphorin 4A1q22SYBR GreenF- ATGCCCAGGGTCAGATACTAT, R-TTCTCCGAGATCCTCTGTTTCN/A1770.997
CCL18Chemokine (C–C motif) ligand 1817q11.2SYBR GreenF-TTCCTGACTCTCAAGGAAAGG, R-CTGGCACTTACATGACACCTGN/A2091.006
PCDH11XProtocadherin 11 X-linkedXq21.3SYBR GreenF-TCTTTTGGTCAGTGTTGTGCG, R-CAACAAGTCGCCTATCAGGACN/A1880.993
APP, Exon 14Amyloid precursor protein21q21.3TaqManCGGTCAAAGATGGCATGAGAGCATC*, Assay Hs01255859_cnFAM-MGB911.040
APP, Exon 3Amyloid precursor protein21q21.3TaqManF-GCACTTCTGGTCCCAAGCAT, R-CCAGTTCTGGATGGTCACTGROX-IB1400.992
SEMA4ASemaphorin 4A1q22TaqManGTTCAAGGGTATGTGAGGTGAGATG*, Assay Hs00329046_cn_VICVIC-MGB901.016
  1. *

    Denotes probe sequence provided by Life Technologies.

Standard curves were used to determine the amplification efficiency of all primer sets used for qPCR. For standard qPCR primers, curves were created by serially diluting purified pGEM-T Easy plasmid DNA (Promega, Madison, WI) containing a single copy of the gene of interest (Pfaffl, 2001); serial dilutions of genomic DNA were used for TaqMan primers (D'Haene et al., 2010). DNA concentrations were converted to gene copy number by calculating the weight (in g/mol) of the DNA used for generating the standard curve. A linear regression of the curve comparing the log of the gene copy number vs the crossing threshold (Ct) of the primer set was determined from primer efficiency (E) = 10−1/slope − 1 (Pfaffl, 2001). Only standard curves with R2 values of greater than 0.99 were used.

Standard qPCR on small neuronal populations

Request a detailed protocol

Genomic DNA from bulk-sorted nuclei (described above) was isolated using the DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA) and quantified using Quant-iT PicoGreen (Life Technologies, Carlsbad, CA); genomic DNA was stored at −20°C before use. Standard qPCR reactions using SYBR Green (Promega, Madison, WI) fluorescence detection (ex: 494 nm; em: 529 nm) were performed in triplicate using 0.5 ng of sample gDNA per reaction. Reactions were run on a Rotor-Gene RG-3000 72-well thermocycler (Qiagen, Valencia, CA) using GoTaq qPCR master mix (Promega, Madison, WI) and the following parameters: denaturation (95°C for 5 min), amplification (95°C for 25 s, 59°C for 30 s and 72°C for 30 s), and quantification through 40 cycles; and a melting curve determination (55–99°C, 30 s on the first step, 5 s for each subsequent step). The crossing threshold (Ct) was determined for each primer set within the linear region of the amplification curve. Down Syndrome nuclei were used as a control.

Microfluidics-based qPCR on single neuronal nuclei

Request a detailed protocol

Single neuronal nuclei from FACS were sorted directly into a 96-well plate containing QuickExtract DNA Extraction Solution (Epicentre, Illumina, San Diego, CA) according to the manufacturer’s instructions. Multiple independent sorts were completed for each group and each individual. Prior to analysis on the Fluidigm Biomark HD (South San Francisco, CA), single neuron genomic DNA was pre-amplified as per Fluidigm protocols (Fan and Quake, 2007; Dube et al., 2008; Qin et al., 2008; Jones et al., 2011; White et al., 2011; Whale et al., 2012) on a Veriti thermocycler (Life Technologies, Carlsbad, CA) using locus-specific Taq-man primer sets (primer Table above) (20× initial concentration, 18 µM primer, 5 µM probe; combined and diluted to 0.2×) (95°C denaturation for 5 min; 18 amplification cycles of a 95°C denaturation for 15 s, followed by a 60°C annealing and extension step for 4 min; and a final extension step at 72°C for 7 min). Locus-specific pre-amplification was confirmed on a Roche LightCycler (Roche Applied Science, Indianapolis, IN) using one targeted primer set prior to analysis on the Biomark 48.48 Dynamic Array integrated fluidic circuit (IFC) (Fluidigm, South San Francisco, CA). Samples were diluted 1:5 and loaded into the 48.48 Dynamic Array IFC according to the manufacturer's protocol. DNA was loaded in triplicate and assays in sextuplicate for a total of 18 replicates per assay per nucleus. Samples were run across multiple arrays for quality control between runs (Table 3) and multiple individuals were run on each array. The thermocycling program was performed on the Biomark: 95°C for 10 min, then 55 cycles of 95°C denaturation, and 60°C annealing and extension. Fluorescent probes used for these assays were 5′-FAM or 5′-VIC with a 3′-minor groove binding (MGB) non-fluorescent quencher, or 5′-ROX with a 3′ Iowa Black quencher. Ct values were determined using Fluidigm's Real-Time PCR Analysis Software. Only nuclei with 10 or more replicates per assay were used for analysis.

Table 3

Quality control between 48.48 Dynamic Array runs

https://doi.org/10.7554/eLife.05116.027
Cell 1Cell 2Cell 3Cell 4Cell 5Cell 6Cell 7Cell 8
APP 14Run 121.3718.7718.0218.3417.1119.3320.1218.02
Run 220.6018.7718.1418.4117.2519.5020.6018.26
APP3Run 122.7920.3518.5319.11
Run 221.9719.8418.6519.28
SEMA4ARun 125.5424.2922.4424.2322.8125.7424.4624.67
Run 225.2624.5323.9224.3923.0325.4224.1124.30

Quantitative PCR analysis

Request a detailed protocol

For both qPCR strategies, relative copy number (RCN) for a diploid sample was calculated as RCN = 2 × (1 + E)−ΔΔCt where E is primer efficiency (E = 10−1/slope − 1) (Weaver et al., 2010; Pfaffl, 2001; D'Haene et al., 2010; Livak and Schmittgen, 2001). Paired cerebellar samples were used as a calibrator to determine ΔΔCt values. For single cell RCN determinations, cerebellar ΔCt values for the paired cerebellum were averaged. RCNs were modeled for each copy number 1–6, assuming a system standard deviation of 0.25 and a 95% CI equal to the standard error of the mean multiplied by the critical t value for a two-tailed t-distribution (degrees of freedom = 68) with p = 0.05 (Weaver et al., 2010). 95% confidence interval (CI) was also determined for the RCN of each assay for each single nucleus giving RCN = 2 × (1 + E)−ΔΔCt ± CI, and the upper and lower bounds were used to call copy numbers for APP exons 3 and 14 such that a CI that overlapped with the modeled RCN range was considered as belonging in that copy number bin (See Table 4 for modeled upper and lower bounds). Copy number bins were 1–6 and >6, as beyond 6 copies the CIs for the given degrees of freedom begin to overlap, allowing for assessments of significant increases in copy number but limiting the ability to distinguish between nuclei with, for example, 8 vs 9 copies. There was high concordance between the APP exons 3 and 14 demonstrating minimal amplification bias between the primer sets.

Table 4

Confidence Intervals for (CI) calling Copy Number (CN)

https://doi.org/10.7554/eLife.05116.028
CNCIRCN, value
1Lower0.92156
Upper1.08512
2Lower1.84312
Upper2.17023
3Lower2.76468
Upper3.25535
4Lower3.68624
Upper4.34047
5Lower4.60780
Upper5.42559
6Lower5.52936
Upper6.51070

Statistical analysis

Request a detailed protocol

For DNA content FCM, DNA indices (DIs) were determined by taking the ratio of the mean from the diploid (2N) peak from the brain sample to the mean of the lymphocyte control peak, both normalized to the mean of the CEN standard (Darzynkiewicz and Huang, 2004; Darzynkiewicz et al., 2004; Westra et al., 2010). The percent change was calculated from DI values assuming a 2N diploid would have a DI of 1. p values for comparison of mean percent change, skew and coefficients of variation comparisons were determined by one-way ANOVA and Tukey's multiple comparison tests. Linear regression analysis was used to determine age-percent change correlation. For comparison of average percent change from NeuN-positive vs negative cortical nuclei, or NeuN-positive cortical and cerebellar nuclei, p values were determined by unpaired, two-tailed t-test. For standard qPCR analyses, differences in ΔΔCt ± SEM of APP in the cortex vs cerebellum were assessed in each individual using an unpaired, two-tailed t-test. For single cell qPCR, p values were determined by one-way ANOVA and Tukey's multiple comparison tests.

Fluorescent in situ hybridization

Request a detailed protocol

Isolated nuclei were stained with DAPI and hybridized using dual color FISH as described previously (Rehen et al., 2001; Kaushal et al., 2003). FISH paints against the whole q arm of chromosome 21 and a point probe against a region on the q arm of 21 (21q22.13-q22.2) (Vysis. Downer's Grove, IL) were used. The mounted slides were examined on a Zeiss Axioskop microscope and Axiocam CCD camera (Carl Zeiss, Thornwood, NY). Approximately, 500 nuclei were blindly counted for each brain, by two independent observers, on 14 samples (5 non-diseased, 9 diseased). Total AD cortical nuclei were examined using a highly liberal protocol for calling aneusomies whereby borderline FISH profiles suggestive of aneusomy were always included in quantitative assessments. All analyses were conducted blind to the identity of samples by interrogating purified nuclei. The ability of this technique to detect aneuploidy was validated using interphase nuclei from a human trisomy 21 brain (Svendsen et al., 1998) revealing three nuclear signals (Figure 5D,E).

PNA probe design and dot blot confirmation

Request a detailed protocol

Peptide nucleic acid probes were custom designed in coordination with PNA Bio, Inc. Nucleic acid sequences were identified and analyzed in silico to ensure binding to only one specific genomic region on APP. Nine unique probes (4 on exon 3 and 5 on exon 14) were designed and conjugated to a single fluorophore of Alexa-488 (Table 5). Specificity of probes was confirmed using dot blot DNA detection paired with immunoblot via antibodies against the Alexa-488 fluorophore (AF-488) (Figure 5—figure supplement 1A). Increasing DNA concentration (1.8, 3.6 and 5.4 μg) of plasmids containing one copy of each exon (all nine PNA binding sites) was used to verify linear increases in AF-488 signal (Figure 7—figure supplement 1B). For the APP copy number curve, DNA concentration remained consistent but plasmids containing 0, 3, 6, and 9 copies of all APP PNA binding sites were used. DNA was denatured in 0.1 M NaCl at 50°C and dot blotted onto a positively charged nylon membrane in distilled deionized water (DDW) and washed in 6× SSC. DNA was cross linked to the membrane, and PNA probes were hybridized to the membrane using conditions consistent with probe hybridization to nuclei on slides. Briefly, probes were prepared in 20 mM Tris, 60% formamide, and 0.1 μg/ml salmon sperm. Probes and membrane were heated 85°C for 5 min, and probes were added to the membrane and incubated at 85°C for 10 min, followed by room temperature overnight. Membrane was then washed twice at 60°C in 2× SSC + 0.1% Tween-20, and then in successive washes of 2× SSC, 1× SSC, and DDW. Probes were visualized on the membrane using a Typhoon fluorescence scanner from General Electric (Fairfield, CT). For quantification, color data were removed, image was inverted, brightness/contrast was adjusted evenly across the entire image, and average pixel intensity was acquired for a region of interest with a standardized area across sample comparisons.

Table 5

Peptide nucleic acid (PNA) probe sequences

https://doi.org/10.7554/eLife.05116.029
GeneProteinLocusAssaySequence
APP Exon 3Amyloid precursor protein21q21.3PNA FISHA488-GATGGGTCTTGCACTG, A488-CCCCGCTTGCACCAGTT, A488-GGTTGGCTTCTACCACA, A488-CAGTTCAGGGTAGAC
APP Exon 14Amyloid precursor protein21q21.3PNA FISHA488-CTCCATTCACGG, A488-GTGGTTTTCGTTTCGGT, A488-ACTGATCCTTGGTTCAC, A488-ACTGATCCTTGGTTCAC, A488-ACGTCATCTGAATAGTT
TelomereN/ATelomeresPNA FISHTelC-Cy3 (F1002, PNA BIO)

Peptide nucleic acid fluorescent in situ hybridization

Request a detailed protocol

Probes were hybridized according to the manufacturer's instructions (Lansdorp et al., 1996). Neuronal nuclei from non-diseased and AD brains were sorted for NeuN positivity and dried onto slides. Slides were washed in PBS, fixed with 4% PFA, and treated with RNase (Qiagen, Valencia, CA) for 20 min at 37°C. Slides were then digested with 200 µg/ml proteinase K (Roche Applied Science, Indianapolis, IN) for 5 min at 37°C. Slides were dehydrated in ethanol series and denatured at 85°C for 5 min. PNA probes in hybridization buffer (20 mM Tris pH 7.4, 60% formamide, 0.1 µg/ml salmon sperm DNA) were then added for 10 min at 85°C. Slides were then removed and placed at room temperature for 2 hr. Slides were then washed twice in 2× SSC + 0.1% tween at 60°C for 30 min each and mounted using progold mounting media. Z-stacks were acquired using a Nikon N-SIM (structured illumination microscopy) super resolution microscope with an Andor iXon3 back-illuminated high sensitivity EMCCD camera with single photon detection capability. Projections were rendered using 3D-SIM Elements software from Nikon. Counts were averaged and analyzed by unpaired, two-tailed t-test

References

    1. Blaschke AJ
    2. Staley K
    3. Chun J
    (1996)
    Widespread programmed cell death in proliferative and postmitotic regions of the fetal cerebral cortex
    Development 122:1165–1174.
    1. Copani A
    2. Sortino MA
    3. Caricasole A
    4. Chiechio S
    5. Chisari M
    6. Battaglia G
    7. Giuffrida-Stella AM
    8. Vancheri C
    9. Nicoletti F
    (2002)
    Erratic expression of DNA polymerases by beta-amyloid causes neuronal death
    FASEB Journal 16:2006–2008.
    1. Darzynkiewicz Z
    2. Huang X
    (2004)
    Analysis of cellular DNA content by flow cytometry
    Current Protocols in Immunology, Chapter 5:Unit 5 7, 10.1002/0471142735.im0507s60.
    1. Darzynkiewicz Z
    2. Juan G
    3. Srour EF
    (2004)
    Differential staining of DNA and RNA
    Current Protocols in Cytometry, Chapter 7:Unit 7 3, 10.1002/0471142956.cy0703s30.
  1. Book
    1. Dube S
    2. Mir A
    3. Jones RC
    4. Ramakrishnan R
    5. Sun G
    (2008)
    Computation of maximal resolution of copy number variation on a nanofluidic device using digital PCR
    South San Francisco, CA: Fluidigm corporation.
    1. Herrup K
    2. Arendt T
    (2002)
    Re-expression of cell cycle proteins induces neuronal cell death during Alzheimer's disease
    Journal of Alzheimer's Disease 4:243–247.
    1. Kaushal D
    2. Contos JJ
    3. Treuner K
    4. Yang AH
    5. Kingsbury MA
    6. Rehen SK
    7. McConnell MJ
    8. Okabe M
    9. Barlow C
    10. Chun J
    (2003)
    Alteration of gene expression by chromosome loss in the postnatal mouse brain
    The Journal of Neuroscience 23:5599–5606.
    1. Potter H
    (1991)
    Review and hypothesis: Alzheimer disease and Down syndrome–chromosome 21 nondisjunction may underlie both disorders
    American Journal of Human Genetics 48:1192–1200.
    1. Yang Y
    2. Geldmacher DS
    3. Herrup K
    (2001)
    DNA replication precedes neuronal cell death in Alzheimer's disease
    Journal of Neuroscience 21:2661–2668.

Article and author information

Author details

  1. Diane M Bushman

    1. Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    2. Biomedical Sciences Graduate Program, University of California, San Diego, La Jolla, United States
    Contribution
    DMB, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Gwendolyn E Kaeser
    Competing interests
    The authors declare that no competing interests exist.
  2. Gwendolyn E Kaeser

    1. Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    2. Biomedical Sciences Graduate Program, University of California, San Diego, La Jolla, United States
    Contribution
    GEK, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Contributed equally with
    Diane M Bushman
    Competing interests
    The authors declare that no competing interests exist.
  3. Benjamin Siddoway

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Contribution
    BS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  4. Jurgen W Westra

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Present address
    Genedata Inc., Lexington, United States
    Contribution
    JWW, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  5. Richard R Rivera

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Contribution
    RRR, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  6. Stevens K Rehen

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Present address
    National Laboratory for Embryonic Stem Cell Research, Federal University of Rio de Janeiro, Institute of Biomedical Sciences, Rio de Janeiro, Brazil
    Contribution
    SKR, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  7. Yun C Yung

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Contribution
    YCY, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  8. Jerold Chun

    Department of Molecular and Cellular Neuroscience, Dorris Neuroscience Center, The Scripps Research Institute, La Jolla, United States
    Contribution
    JC, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    jchun@scripps.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institutes of Health (MH076145)

  • Jerold Chun

Shaffer Family Foundation

  • Benjamin Siddoway
  • Richard R Rivera
  • Jerold Chun

National Institutes of Health (T32 DA07315)

  • Diane M Bushman
  • Gwendolyn E Kaeser

National Institutes of Health (T32 AG 000216)

  • Gwendolyn E Kaeser

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Drs Don Cleveland, Edward Koo, and Kun Zhang for helpful comments and suggestions; Drs Olga Pletnikova and Juan Troncoso at the Johns Hopkins School of Medicine Alzheimer's Disease Research Center for providing human brains; the Flow Cytometry Core at TSRI for FACS expertise; Hope Mirendil for help with qPCR and manuscript comments; Mr Mathew McCreight for FISH quantitation; Uma Dandakar and Hemani Wijesuriya of the UCLA GenoSeq Core; Kathy Spencer, Michael Thieleking, and Ariana Lorenzana of The Nikon Center of Excellence at TSRI; Caroline Dando and Susan Robelli of Fluidigm for advice on Fluidgm usage; Andrew Richards for advise on WGA experiments; and Danielle Jones for manuscript editing. This work was partially funded by the NIH (MH076145 to JC), (T32 DA07315 to DMB and GEK), (T32 AG 000216 to GEK), The Shaffer Family Foundation and private sources (to BS, RRR and JC). The UCI-ADRC is funded by NIH/NIA Grant P50 AG16573.

Copyright

© 2015, Bushman et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 7,831
    views
  • 986
    downloads
  • 99
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Diane M Bushman
  2. Gwendolyn E Kaeser
  3. Benjamin Siddoway
  4. Jurgen W Westra
  5. Richard R Rivera
  6. Stevens K Rehen
  7. Yun C Yung
  8. Jerold Chun
(2015)
Genomic mosaicism with increased amyloid precursor protein (APP) gene copy number in single neurons from sporadic Alzheimer's disease brains
eLife 4:e05116.
https://doi.org/10.7554/eLife.05116

Share this article

https://doi.org/10.7554/eLife.05116

Further reading

    1. Computational and Systems Biology
    2. Neuroscience
    Anna Cattani, Don B Arnold ... Nancy Kopell
    Research Article

    The basolateral amygdala (BLA) is a key site where fear learning takes place through synaptic plasticity. Rodent research shows prominent low theta (~3–6 Hz), high theta (~6–12 Hz), and gamma (>30 Hz) rhythms in the BLA local field potential recordings. However, it is not understood what role these rhythms play in supporting the plasticity. Here, we create a biophysically detailed model of the BLA circuit to show that several classes of interneurons (PV, SOM, and VIP) in the BLA can be critically involved in producing the rhythms; these rhythms promote the formation of a dedicated fear circuit shaped through spike-timing-dependent plasticity. Each class of interneurons is necessary for the plasticity. We find that the low theta rhythm is a biomarker of successful fear conditioning. The model makes use of interneurons commonly found in the cortex and, hence, may apply to a wide variety of associative learning situations.

    1. Neuroscience
    Bharath Krishnan, Noah Cowan
    Insight

    Mice can generate a cognitive map of an environment based on self-motion signals when there is a fixed association between their starting point and the location of their goal.