Abstract
Amyloid plaques are one of the hallmarks of Alzheimer’s disease (AD). The main constituent of amyloid plaques is amyloid-β peptides, but a complex interplay of other infiltrating proteins also co-localizes. We hypothesized that proteomic analysis could reveal differences between amyloid plaques and adjacent control tissue in the transgenic mouse model of AD (APPPS1-21) and in similar regions from non-transgenic littermates. Our microproteomic strategy included isolation of regions of interest by laser capture microdissection and analysis by liquid chromatography mass spectrometry-based label-free relative quantification. We consistently identified 183, 224, and 307 proteins from amyloid plaques, adjacent control and non-tg samples, respectively. Pathway analysis revealed 27 proteins that were significantly regulated when comparing amyloid plaques and corresponding adjacent control regions. We further elucidated that co-localized proteins were subjected to post-translational modifications and are the first to report 193 and 117 unique modifications associated to amyloid plaques and adjacent control extracts, respectively. The three most common modifications detected in proteins from the amyloid plaques were oxidation, deamidation, and pyroglutamylation. Together, our data provide novel information about the biological processes occurring within and around amyloid plaques in the APPPS1-21 mouse model of AD.
INTRODUCTION
Alzheimer’s disease (AD) is a neurodegenerative disorder and the most common cause of dementia [1, 2]. The underlying mechanisms of the disease is not fully understood although two pathological hallmarks are deposition of extracellular plaques, mainly composed of amyloid-β (Aβ) peptides [3, 4], and intracellular neurofibrillary tangles, mainly composed of hyperphosphorylated tau [5, 6]. Aβ peptides originate from sequential cleavage of the membrane-spanning amyloid-β protein precursor (AβPP) by different secretases. The widely accepted ‘amyloid cascade hypothesis’ postulates that deposition of Aβ into senile plaques (Aβ plaques) initiates a pathological cascade leading to AD [7–11]. The Aβ deposition is thought to begin prior to tau-associated pathology and clinical symptoms of AD underlining the pathological relevance of Aβ proteoforms in AD [12]. The molecular composition of Aβ plaques encompasses a variety of different amyloid-binding and infiltrating proteins [13]. Identification of these co-localizing proteins are reasoned to contain valuable information about processes occurring within and around Aβ plaques. It has been shown that aggregated Aβ peptides contain several post-translational modifications (PTMs) such as oxidation- and pyroglu-modifications [14–16]. However, it remains to be elucidated whether plaque-infiltrating and associating proteins also are subjected to the same degree of PTMs as aggregated Aβ peptides.
Mass spectrometry (MS)-based technologies have in recent years proven to be a valuable tool to achieve in-depth characterization of the proteomic profile of biological specimens and characterization of PTMs [17, 18]. A commonly used method for proteomic discovery is a bottom-up approach in combination with data-dependent acquisition (DDA) [19]. The method is label-free and enables large coverage of the proteome that presents several advantages for elucidating complex biological systems, including identification of proteins containing PTMs. Although detection and validation of PTMs in MS presents several challenges such as low stoichiometry, site specific assignment and poor stability, recent improvements of algorithms enables a strong PTM identification and quantitative analysis by integrating the database search with results from de novo based analysis [20]. Furthermore, a recent study has shown in-depth PTM identification from a DDA approach which supports the applicability to our study [21].
By using laser capture microdissection (LCM), targets of interest within tissue sections can be isolated and subjected to MS analysis providing localized information. A recent study has illustrated the capability of using LCM to capture different Aβ plaque profiles when staining with tetra- and heptameric formyl-thiophene acetic acid (q/h-FTAA), thioflavin-S, or Congo red [22]. Another example of this approach was applied to Aβ plaques from rapidly progressive AD and sporadic AD patients [23].
Transgenic (tg) mice, like the APPPS1-21 strain [24], are useful mechanistic models that are based on familial AD (FAD) mutations which result in an increased Aβ load and early-onset inherited AD. The APPPS1-21 model carries mutations in the APP gene (KM670/671NL, Swedish) and in the PSEN1 gene (L166P), leading to excessive production of Aβ [24]. Despite differences in the senile plaque composition when comparing mouse and human brain as described by others [25] and us [26] previously, the main response mechanism to Aβ accumulation which is of interest for our study here are considered very similar. For example, Aβ plaque deposition and subsequent activation of microglia and astrogliosis begins at approximately six weeks of age in the mouse cerebral cortex [24], demonstrating glial engagement that similarly is observed in AD. Furthermore, dendritic spine loss is observed approximately four weeks after plaque deposition [27], indicating that the tg mouse model also resembles part of the neurodegenerative aspect of AD. Based on these comparable features, it is reasoned that APPPS1-21 mice are relevant to study Aβ-related pathologies and model how abnormal Aβ production affects the brain.
We used a discovery-based proteomic strategy in combination with LCM to profile the proteome of Aβ plaques and closely surrounding regions extracted from the cerebral cortex of APPPS1-21 mice. We discovered several differentially regulated proteins in Aβ plaques compared to adjacent control tissue by using two independent search algorithms. Additionally, we are the first to elucidate the unique PTM profile of the Aβ plaque associated proteins compared to the surrounding tissue.
MATERIALS AND METHODS
Materials
Acetonitrile (ACN; ≥99%), tris(2-carboxyethyl)phosphine (TCEP), 2-chloroacetamide (CAA; ≥98%), and triethylamonium bicarbonate (TEAB) were purchased from Sigma Aldrich, Denmark. Formic acid (FA; ≥98%) was purchased from Fluka, Denmark, thioflavin-t from Merck Eurolab, Denmark and trifluoroacetic acid (TFA; ≥99%) from Fisher Scientific, Denmark. ProteaseMax™ was purchased from Promega, Denmark.
Animals
Six 12-month-old heterozygous female tg APPPS1-21 and corresponding non-transgenic (non-tg) littermates were used for the analysis (Charles River, Germany). The mice were group-housed (five/cage). Water and food were supplied ad libitum, clean cages were provided twice a week with enrichment and mice were kept on a 12 h/12 h light/dark schedule. Room humidity was 55% ±5%, and room temperature (RT) was 21±2°C. The animal experiments were performed in accordance with the European Communities Council Directive no. 86/609, the directives of the Danish National Committee on Animal Research Ethics, and Danish legislation on experimental animals (license no. 2014-15-0201-00339). Mice were anesthetized with avertin (i.p. administration of 2% solution; 1 ml/100 g) and perfused with KPBS/heparin. Brains were fresh frozen on dry ice and stored at -80°C until further use.
Tissue sectioning
The mouse brains were embedded with Tissue-Tek OCT (Sakura, Denmark) and left for temperature adaption at -16°C for 30 min. The samples were cut into 25-μm thick sagittal sections using a cryostat microtome (Leica CM3050 S, Denmark) and mounted on polyethylene naphthalate (PEN) membrane glass slides (Thermo Fisher, Denmark). Three sections from each hemisphere were applied to a PEN membrane glass slide and fixed with 70% EtOH and 99% EtOH for one minute in each solution. Next, they were dried at RT and stored at -80°C until further use. Adjacent sections from both non-tg littermates and APPPS1-21 mice were similarly attached to SuperFrost Plus (VWR, Denmark) glass slides for thioflavin-T and immunohistochemistry analysis. Here, sections used for thioflavin-T staining were similarly fixed in EtOH while sections used for immunohistochemistry were not.
Thioflavin-T staining and immunohistochemistry
Brain sections from both APPPS1-21 and non-tg littermates were stained in 1% thioflavin-T solution for 10 min. The thioflavin-t solution was removed by rinsing the slides in ddH2O three times for 3 min each. The tissue was air dried at RT and mounted in pertex (Histolab, Denmark). Adjacent sections from the APPPS1-21 mice were immunolabeled using the Vectastain Elite ABC HPR kit (Vector Laboratories Inc., Denmark). In brief, the fresh frozen tissue sections were fixed in 4% PFA for 10 min and pretreated with 1% hydrogen peroxide (H2O2) solution to block endogenous peroxidase activity and 90% formic acid solution for antigen retrieval. Washes in KPBS was performed in between. Sections were permeabilized in blocking buffer (1% bovine serum albumin (BSA), 0.3% triton X-100, 5% normal swine serum in KPBS) followed by over-night (O/N) staining with the anti-Aβ antibody (mouse, 6E10; 1 : 500; BioLegend) at 4°C. Tissue was counterstained with biotinylated anti-mouse biotinylated antibodies (DAKO E0354; 1 : 500), followed by incubation in the avidin-biotin complex (ABC) solution (Vector Laboratories Inc., Denmark). Sections were washed in washing buffer (0.25% BSA, 0.1% triton X-100 in KPBS) in between. The immunoreactivity was visualized with the chromogen, diaminobenzidine (DAB, Sigma, Denmark). Imaging was performed on a Zeiss Axio Scope.A1 equipped with a Zeiss Axiocam 506 color camera.
LCM of mouse brain tissue sections
The PEN membrane glass slides containing tissue sections were defrosted by RT adjustment and a dip wash in ddH2O prior to staining with a 1% thioflavin-T solution, as stated in the previous section (except mounting in pertex). Three types of regions (Aβ plaques and adjacent control regions from the APPPS1-21 mice, and similar control areas from the non-tg littermates) were microdissected and collected onto CapSure™ Macro LCM caps (Thermo Fisher, Denmark) using the Arcturus Veritas LCM system. A 20X magnification was utilized and 50 regions were extracted from each section with an average area of 5791.2 μm2 per extracted plaque, 6449.4 μm2 per extracted adjacent control, and 5560.7 μm2 per extracted non-tg littermate regions. The microdissection focused on core plaques, and not diffuse plaques, because the latter are mainly thioflavin-T negative. The caps containing the microdissected material were incubated in 50 μl digestion buffer (0.01% ProteaseMax™ in 50 mM TEAB buffer, pH 8.5) O/N at 37°C. 10 μl of fresh digestion buffer was manually resuspended five times on the cap surface to ensure release of potential remaining material (total volume = 100 μl) and collected in 1.5 ml Eppendorf LoBind microcentrifuge tubes (Eppendorf). An in-solution digest protocol was utilized for each tube. In brief, the tubes were sonicated for 2 min, heat treated at 95°C for 5 min and allowed to cool to RT. The samples were reduced and alkylated in 1 : 50 TCEP (v:v; 0.5M) and 1 : 10 CAA (v:v; 0.5M) and incubated in darkness for 30 min at 37°C. One μl trypsin (0.1 μg/μl) was added and incubated O/N at 37°C. TFA (final concentration = 1%) was added and incubated for 30 min at 37°C to stop the reaction of ProteaseMax™ and trypsin. The tubes were inserted into a SpeedVac (miVac, GeneVac), dried for 6 h and stored at -80°C until analysis. The dried material was reconstituted in 12 μl loading buffer (2% ACN, 0.1% TFA, 0.1% FA) prior to MS analysis.
Data acquisition by LC-MS/MS
A pooled sample was created for each subregional group by adding two μl from each reconstituted sample into one 0.5 ml tube (LoBind, Eppendorf). These samples were utilized to create two pooled higher intensity samples per subregion in the data analysis. Nine μl sample material was added to 96-well plate with a random order. Samples were analyzed on quadrupole Orbitrap MS (Q-Exactive HF-X, Thermo Scientific) connected to an RSLC UHPLC nanopump system (Dionex). Runs of each sample were loaded onto the C18 reversed phase column (Dionex; Acclaim PepMap100 C18, 5 μm precolumn and 75 cm Acclaim Pepmap RSLC, 75 μm ID main column, Thermo Scientific) and eluted with a linear gradient of 96% solvent A (1% formic acid) and 4% solvent B (ACN) [28], increasing solvent B to 35% on a 30-min ramp gradient. The MS was operated in DDA, selecting the 10 precursor-ions with the highest intensity for higher energy collisional dissociation (HCD) fragmentation. Full MS settings were as follows: resolutions = 60 000; AGC target = 3e6 (max 60 ms); scan range 375 to 1500 m/z. The following dd-MS2 settings were used: resolution = 30 000; AGC target = 1e5 (max 120 ms); loop count = 10; isolation window = 1.5 m/z; scan range 200 – 2000 m/z. The maximum IT was set to 120 ms for all samples except the two pooled higher intensity samples which were acquired at 250 ms and 500 ms, respectively (Supplementary Figure 1). The following dd settings were used: underfill ratio = 30.0%; charge exclusion unassigned, 1, 7, 8,>8; peptide match = preferred; exclude isotopes = on. Intensity threshold was set to 2.5e5 for all samples except the two pooled higher intensity samples which were acquired at 1.2e5 or 6.0e4. Dynamic exclusion was set to 10.0 s for all samples except the two pooled higher intensity samples which were either excluded or set 5.0 s, respectively. Data are available via ProteomeXchange with identifier PXD014274.
Data processing
Datafiles were processed in MaxQuant (v1.6.5.0) [29]. The FASTA file, containing the UniProt mus musculus database (UP_000000589), was downloaded 4 September 2017 and contained 16.911 entries. The human amyloid beta A4 protein/amyloid precursor protein (APP) and presenilin 1 (PSEN1) sequences were added due to the co-expression of these proteins in the APPPS1-21 mouse model while the mouse amyloid beta A4 protein/amyloid precursor protein was removed. The Swedish mutation (KM670/671NL) was not included due to its location prior to the Aβ1-42 peptide sequence. Bioinformatic and statistical analysis was performed in Perseus (v1.6.5.0) [30]. Potential contaminants and reverse sequences were removed followed by filtering the dataset by minimum 50% valid values and presence of > 1 unique peptide. Hierarchical clustering was applied together with Pearson correlation analysis to identify potential outliers. A sample was defined as an outlier when two criteria were met; 1) the average Pearson correlation coefficients (<0.8, compared to the corresponding samples in the group (Aβ plaque versus Aβ plaque, adjacent control versus adjacent control, and non-tg versus non-tg), and 2) a median absolute deviation (MAD) analysis (if the value was five times away from the median it was defined as an outlier). The normal distributed data (histograms not shown) was analyzed with two-sampled test (two-sided student’s t-test; minimum valid value filter = 70%; Permutation-based FDR = 0.05). The same MaxQuant settings, except the match-between-run function, were used to elucidate the consequence of excluding this feature. STRING (v11.0) enrichment analysis was used for pathway analysis [31]. The five Biological Process (GO)-associations with lowest FDR were selected for both up- and downregulated proteins. Variation of the second component in the PCA plot was analyzed with student t-test and Fisher exact test using the Gene Ontology Biological Process (GOBP) class.
A parallel data analysis was performed in BSI PEAKS Studio X (v10; Bioinformatics Solution; Canada). All data samples were imported, analyzed with PEAKS Denovo, PEAKS DB search, PEAKS PTM and PEAKS SPIDER. The following settings were used: parent mass error tolerance 10.00 ppm; fragment mass error tolerance = 0.02 Da; fragment mass search type monoisotopic; trypsin with [D, P]; max missed cleavages 3; semi specific digestion mode; fixed modification: Carbamidomethylation (57.02); Variable modification: citrullination (0.98), oxidation (M; 15.99), phosphorylation (STY; 79.97); max variable PTM per peptide = 3. The database that was used for the MaxQuant previously was similarly used for the PEAKS analysis. The following settings were applied for visualization: FDR = 0.1%; PTM ion intensity≥1%; Proteins = 10logP≥20; unique peptides≥1; De novo ALC≥50%; -10logP≤18.4. Similar settings were applied to PEAKS PTM analysis although the variable modifications now included 312 different modifications.
The label-free quantification (LFQ) analysis of Aβ plaques and adjacent control was performed by using the PEAKS DB and PEAKS PTM searches as references (using 0.1% FDR; significance≥20; fold change≥1; significance method = PEAKSQ; At least 1 unique peptide). LFQ modification outputs were processed using R as reported in [21] to identify unique modifications. In brief, LFQ modification outputs (with minimum a-score of 100) were exported and log2 transformed followed by calculation of the LFQ ratio between the modified and unmodified peptides. A non-zero quantitation value was required in more than half of the samples for inclusion.
RESULTS
Proteomic fingerprint of Aβ plaques in APPPS1-21 model mice
We investigated the protein composition of Aβ plaques from 12-month-old APPPS1-21 mice by a combination of tissue microdissection and label-free LC-MS/MS quantification. APPPS1-21 mice overexpress human APP and PS1 with FAD mutations leading to overexpression of human Aβ species and rapid progression of Aβ plaques in the mouse brain [24]. Tissue was dissected from the cerebral cortex region containing a high Aβ plaques load in the 12-month old APPPS1-21 mice (Fig. 1A). Aβ plaques were identified by thioflavin-T staining exclusively in brains from APPPS1-21 and not in brains from non-tg littermates (Supplementary Figure 2). A 20X objective lens was used for microdissection to ensure a true-plaque morphology and avoid inclusion of potential contaminant caused by the fluorescent background. Additionally, we validated the thioflavin-T staining by orthogonal detection of Aβ plaques by immunohistochemistry with the anti-Aβ 6E10 antibody. To this end, we stained adjacent brain sections with either thioflavin-T or 6E10 and observed that the same plaques in the adjacent sections were detected by both thioflavin-T and 6E10 (Supplementary Figure 2).

A) Experimental workflow. Brain tissue sections from right and left hemisphere were collected from six transgenic (tg) APPPS1-21 and six non-tg mice (total = 12 samples per genotype). Laser microdissection was used to extract regions of interest (50 per brain slice; three slices per hemisphere); amyloid plaques (Aβ plaque, green), adjacent control (blue), or non-tg (red), respectively. Extracted proteins were digested, analyzed by liquid chromatography-tandem mass spectrometry, and processed by bioinformatic analysis. B) Principal component analysis plot of LFQ intensities associated with samples (after Pearson correlation correction and median absolute deviation analysis). Samples from both right and left hemispheres are presented (n = 6 mice; samples from both hemispheres). Components 1 and 2 are presented. C) Venn-diagram of the number of proteins observed in different subgroups; Aβ plaque extracted material (green), adjacent control (blue), corresponding regions from non-tg littermate mice (red). D) Volcano plot of differentiating log2 transformed proteins when comparing Aβ plaques and adjacent control regions. The location of amyloid beta A4 protein (App), apolipoprotein E (Apoe), and clusterin (Clu) are visualized.
We microdissected two different regions from the APPPS1-21 mice either containing Aβ plaques or without Aβ plaques of an adjacent control region (Fig. 1A). From non-tg littermates, we microdissected corresponding control regions. Female mice were used in this study because they have been reported to exhibit a more robust Aβ deposition than male APPPS1-21 mice [32].
For this study, we implemented a detergent enhanced solubilization strategy using low volume acid-labile surfactant, ProteaseMax™, to increase solubility without use of detergents that interferes with the MS analysis. Due to the limited sample material in LCM plaque extracts, we applied a match-between-run strategy matching ID to a pooled higher intensity sample to increase the protein identification rate [29]. This approach included the generation of mixed samples from either of the three subgroups. By adjusting the MS settings for these samples, we could affect the number of full- and MS2 scans and thereby establish spectral references that almost doubled the protein identification (Supplementary Figure 1A-C). The strategy was supported by sequential sample analysis with the same chromatographic method. Further, sample concentrations were in the nanogram range, meaning that high loading effect was not an issue as it was significantly lower than the column capacity. The applied software uses normalization in the LFQ algorithm for label-free quantification ensuring a reliable comparison between protein levels [33].
Investigation of the three different microdissected regions by principal component analysis (PCA) implied a difference between the three respective subgroups: non-tg, adjacent control, and Aβ plaques from APPPS1-21. We included samples from both right and left hemisphere from each individual animal and used 6 animals per genotype (Fig. 1B; n = 6 mice; samples from both hemispheres; total = 12 data points). One Aβ plaque sample associated to the left hemisphere of one mouse was lost during the MS analysis. Additionally, the MS analysis of another APPPS1-21 mouse (adjacent control and plaque extract from same hemisphere) did not meet the test criteria (Person correlation minimum value and median absolute deviation) whereby it was defined as an outlier and excluded from the further analysis. Therefore, we have 10 datapoints from the Aβ plaque group (green circles in the PCA plot) and 11 datapoints from the adjacent control group (blue squares in the PCA plot). A partial overlap was observed on the first component between non-tg control and adjacent control regions from APPPS1-21, while no overlap was observed between Aβ plaques from APPPS1-21 and non-tg samples. This separation between non-tg and APPPS1-21 plaque samples on the first component was anticipated. Additionally, we observed a separation on the second component in top and bottom location, which was not dependent on the genotype by affecting samples from both non-tg and APPPS1-21 mice equally. Interestingly, both samples from the left and right hemispheres from each individual animal (except for one non-tg sample) were either in the top or in the bottom location for the second component, indicating that this separation was associated to the individual animals (Supplementary Figure 3). We identified significant protein differences in pathways associated to “response to stimulus” and “metabolic processes” by Gene Ontology Biological Process (GOBP) class between the top and bottom samples in all subgroups (data not shown). Based on these findings, we suggest that the observed separation on the second component was dependent on the individual state (e.g., hormonal and stress related factors) of the mice prior to the MS experiment and consider this separation as normal variation detected within mouse cohorts. Our further pathway analysis focused on tissue-associated differences between the Aβ plaque and adjacent control regions from APPPS1-21.
A total of 555 proteins were detected in extracted Aβ plaque regions from APPPS1-21 mice, following a database search against the mouse reference proteome. Of these, 183 proteins were consistently observed across Aβ plaque samples. In the adjacent control and non-tg control regions, 224 and 307 proteins were identified consistently, respectively.
A total of 427 proteins were shared within all groups (Fig. 1 C, middle of circles). Fourteen proteins, Pafah1b2, Arfgef3, Camk2d, Spock2, Ncam2, Slc1a4, Pgrmc1, Sept2, Apoa1bp, Uqcr10, Lamtor1, Clstn1, Hapln1, and Twf2, were unique to Aβ plaques. Four proteins, Atp2b4, App, Apoe, and Vtn, were unique to the Aβ plaque and adjacent control regions (Fig. 1C; overlapping area between the green and blue circle). Eleven proteins were unique to the adjacent control regions: Cmas, Rpn2, Dlg3, Baiap2, Mtnd5, P4hb, Prdx3, Inpp1, Psmc5, Rplp2, and Prkca.
The inclusion of the adjacent control regions allowed a subsequent identification of relative proteomic differences. By comparing adjacent control and Aβ plaques, we identified App, Apoe, and Clu as the top three regulated proteins illustrated by the volcano plot (Fig. 1D). Two-sided t-test analysis identified 27 proteins that were significantly different when comparing the Aβ plaque proteome data to the adjacent control samples from APPPS1-21 (Tables 1 and 2). The differentiating proteins are listed in (Tables 1 and 2). The regulated proteins (both up and down) contained between 2 and 83 unique peptides (Tables 1 and 2).
Label-free quantification (LFQ) analysis of Aβ plaque extracts compared to adjacent control regions (Upregulated proteins identified by MaxQuant and Perseus analysis; two-sided student’s t-test; minimum valid value filter = 70%; Permutation-based FDR = 0.05; p-value cut-off≤0.05). (*)=Proteins identified as significantly regulated using PEAKS Quantification analysis
Label-free quantification (LFQ) analysis of Aβ plaque extracts compared to adjacent control regions (Downregulated proteins identified by MaxQuant and Perseus analysis; two-sided student’s t-test; minimum valid value filter = 70%; Permutation-based FDR = 0.05; p-value cut-off≤0.05). (*)=Proteins identified as significantly regulated using PEAKS Quantification analysis
To support our protein identification analysis, we performed a parallel analysis with an de novo and tandem MS based identification algorithm, including all sample extracts. The approach identified 2,032 proteins when analyzing the three subgroups together. Furthermore, the additional unbiased PTM module identified 1,397 proteins that contained PTMs (Ascore≥20; FDR = 0.1%; unique peptides≥1; Supplementary Table 1). Moreover, the analysis enabled a ratio determination of truncated Aβ peptides between Aβ plaques and adjacent control regions (Supplementary Figure 4). We found that truncated Aβ species, such as Aβ29 - 40 oxidized (Ox), Aβ29 - 42, Aβ29 - 42 Ox, Aβ29 - 43, and Aβ29 - 43 Ox, were abundant within Aβ plaques and low in adjacent control extracts, indicating that true-plaque content was microdissected in the Aβ plaque group.
The second LFQ analysis identified 95 significantly regulated proteins when comparing the Aβ plaque and adjacent control regions (Supplementary Table 2). Ten proteins were among the previously 27 significantly regulated proteins identified by the MaxQuant and Perseus analysis and thus supported their presence (Indicated by (*) in Tables 1 and 2). The proteins identified by both search algorithms were Eno1, App, Apoe, Clu, Rab2a, Rab7a, Matr3, Map1b, Map2, and Sptan1. The App, Clu, and Apoe proteins were among the most regulated proteins in both LFQ analysis. App, Apoe, and Clu have previously been identified as susceptible genes for AD in different genome-wide association study (GWAS) [34] and identification of these three top upregulated proteins illustrated the validity of our approach in APPPS1-21 mice.
Regulated plaque-proteins are associated to different biological processes
We investigated a possible relationship between the up- and downregulated proteins by using literature based biological process (GO) enrichment analysis of known protein interactions [31]. In general, we observed that up- and downregulated proteins were associated to different biological processes (Table 3).
Enrichment analysis association of significant differentiating proteins identified by label-free quantification (LFQ) in MaxQuant and Perseus. FDR = false discovery rate
The upregulated proteins were associated to the following biological processes: glycolytic process (n = 5), purine ribonucleoside monophosphate biosynthetic process (n = 5), regulation of amyloid fibril formation (n = 3), endocytosis (n = 7), and vesicle-mediated transport (n = 9). Several of these proteins (e.g., Rab1a, Rab2a, Rab7a, App, Apoe, and Vtn) are hypothesized to be involved in endocytosis and/or phagocytosis and thus indicated degradational processes [35–37].
The downregulated proteins were associated to: negative regulation of protein depolymerization (n = 3), dendrite development (n = 3), negative regulation of supramolecular fiber organization (n = 3), negative regulation of cytoskeleton organization (n = 3), and negative regulation of microtubule depolymerization (n = 2). Both proteins, Map1b and Map2, were represented in each of the associations. Map1b and Map2 are constituents of microtubules, suggesting that part of the downregulated proteins was associated with the cytoskeleton of neurons.
Altogether our analysis therefore indicated a functional difference between the up- and downregulated Aβ plaque proteins where the upregulated proteins were associated to endocytosis and glycolysis and downregulated proteins associated to neuronal cytoskeleton.
Proteins co-localizing with Aβ plaques are subjected to unique modifications
PTMs can alter protein structures and thus have an impact on the biological function. This has previously been observed to occur on aggregated Aβ peptides in AD; however, it remains unclear whether co-localizing proteins are subjected to similar degree of modifications. We are the first to report detection of unique PTMs on proteins co-localizing with Aβ plaque extracts isolated from APPPS1-21 mice. Investigation of the modification profile of detected proteins in the Aβ plaques and adjacent control extracts showed 426 shared PTMs (Fig. 2A). Further, there was no clear discriminative difference when depicting the shared PTM data in a PCA plot (Fig. 2B). Interestingly, a total of 193 and 117 modifications were uniquely identified on proteins co-localizing with Aβ plaques and the adjacent control regions, respectively (Ascore > 100 and minimum representation in half of the subgroup; Fig. 2A, C). The identified unique modifications of co-localizing proteins in Aβ plaque extracts were: 2-OH-ethyl thio-Ser, 3-sulfobenzoic succinimidyl ester, acetylation (N-term; Protein N-term), amidation, amidination of lysines or N-terminal amines with methyl acetimidate, aminoethylcysteine, ammonia-loss (N), biotinylation, carbamylation, carboxylation (DKW), carboxymethyl (KW X@N-term), citrullination, crotonaldehyde, deamidation (NQ), deoxy, demethylation (KR), ethanolation (KR), formylation, ISD (z + 2)-series, methionine replacement by azido homoalanine, methylation (KR), methylpyrroline, oxidation (HW), oxidation (M), proline oxidation to pyroglutamic acid, proline oxidation to pyrrolidinone, pyro-glu from Q, replacement of proton with ammonium ion, sulphone, tryptophan oxidation to kynurenin, tryptophan oxidation to oxolactone (Fig. 2 C; Table 4). The three most common unique modifications detected in proteins from the Aβ plaques extracts were oxidation (n = 56), deamidation (n = 26), and pyro-glu (n = 25).

A) Venn-diagram of detected modifications in amyloid-β (Aβ) plaque (green) and adjacent control (blue) extracted regions. B) Principal component analysis plot of the shared modification profile (LFQ intensities) on proteins associated with Aβ plaques (green) and adjacent control (blue) regions. Components 1 and 2 are presented. C) Heatmap representation of unique modifications from Aβ plaque (n = 193; column 1) and adjacent control material (n = 117; column 2). The detected modifications are listed to the left of the heatmap. The number of unique modifications is listed in each cell and illustrated with a rainbow color scale (purple = 0; red = 60). Non-biological modifications, such as ISD (z + 2)-series, were also included to depict all modification identified by the analysis.
Identified modifications on proteins co-localizing with Aβ plaques extracts
The identified modifications in adjacent control regions were: acetylation (N-term), acetylation (Protein N-term), alpha-amino adipic acid, amidation, carboxylation (DKW), citrullination, deamidation (NQ), glycidamide adduct, hexose (NSY), iminobiotinylation, ISD (z + 2)-series, maleimide, membrane protein extraction, methionine replacement by homopropargylglycine, methylation (KR), oxidation (M), phosphorylation (STY), pyro-glu from Q, sodium adduct, and ubiquitin (Fig. 2 C). The three most common modifications in the adjacent control group were oxidation (n = 49), deamidation (n = 22), and citrullination (n = 9).
Our findings indicated that oxidation and deamidation were the most common unique modification of proteins co-localizing with both Aβ plaque and adjacent control regions. The presence of the pyro-glu modification was identified to be much higher in proteins co-localizing with Aβ plaque compared to adjacent control regions.
DISCUSSION
Aβ plaques have previously been observed to encompass hundreds of different infiltrating proteins [23]. We identified 27 up- and downregulated proteins that associated with Aβ plaques in APPPS1-21 mice. To account for the current pathological state of the APPPS1-21 mice, we compared Aβ plaques to the adjacent control areas instead of comparing to brain areas in non-tg littermate mice reflecting normal conditions. Using this approach, we elucidated local proteomic differences and enabled identification of unique modifications associated to Aβ plaques and the surrounding closely related brain tissue of APPPS1-21 mice. On the technical side, we combined a few innovative methods to enhance the data quality. We used a detergent enhanced solubilization strategy with limited sample transfer, minimizing the potential risk of losing sample material during the workflow while enhancing protein solubility and trypsin digestion efficiency [17]. We further implemented a match-between-run strategy to increase the identification rate.
Our MS analysis identified 183 proteins that were consistently observed across Aβ plaque samples in the cerebral cortex of tg APPPS1-21 mice. These included proteins such as Clu, Apoe, Gapdh, and Pgam1. The respective human genes have been identified as susceptible genes for AD in different GWAS [34, 39]. The APOE gene has been described as the most susceptibility gene for SAD. The gene exist in three types of alleles (ɛ2, ɛ3, and ɛ4) were the ɛ4 allele has been reported to be a strong risk factor for both FAD and SAD [40]. Despite reports of reduced GAPDH activity in AD brain, studies have observed interaction with fibrillar Aβ1 - 42 [41, 42]. CLU and PGAM1 have been reported to be elevated in specific brain regions in AD patients [43, 44].
We detected 10 unique peptides for Clu, 8 for Apoe, 16 for Gapdh, and 10 for Pgam1 in comparison with App itself which was detected with 10 unique peptides, indicating a relevance and presence of these proteins within the Aβ plaques. Based on our findings, a relationship may be indicated in which an overexpression of Aβ in the APPPS1-21 mice can cause co-localization of these proteins to Aβ plaques. However, the proteins that we have detected in our study are translated from genes that are reasoned not to include the single nucleotide polymorphisms (SNPs) that have been identified in GWAS studies. The fact that these proteins co-localize with the plaques in the APPPS1-21 model may therefore indicate a physiological mechanism that occurs in response to the Aβ plaque deposition. For example, Apoe has been found to co-localize with Aβ plaques in tg APP mice and hypothesized to influence the degree of Aβ fibrillization [45] whereas Clu has been reported to interact directly with Aβ peptides [46]. The presence of SNPs may change the amino acid sequence of these proteins, impair the biological function and thereby impair the engagement of Aβ plaques. Notably, several differences between human and mouse Apoe has been reported [47]. However, despite the differences, we observed high degree of co-localization of Apoe within the Aβ plaque extracts compared to adjacent control regions in our study. Previous studies have generated mice expressing the human APOE, and crossed them with APP tg mice [49]. Such a double tg mouse strain would be a valuable tool to test the impact of human Apoe on the plaque pathology in the context of our study.
The majority of identified proteins (96.7%; 177/183; list included 2,422 proteins) has previously been detected in plaque extracts from AD patients when comparing to a complete list of identified proteins [23, 50]. In addition, the consistently observed proteins (69.9%; 128/183; list included 277 proteins) were in line with an earlier study [23] supporting the validity of our findings (data not shown). The previous comparison was based on findings in SAD patients whereas the mouse model, that was used in our study, was based on aggressive FAD mutations [24, 51]. The underlying cause for Aβ plaque generation in SAD is still insufficiently understood, and knowledge gained from mouse models with FAD mutations is a necessary approach. Comparing the identified proteins in the APPPS1-21 (FAD) mouse model here to plaque-associated proteins in SAD described previously, reveals high degree of alignment, despite the potential pathological- and species-related differences. The proteomic similarity on a Aβ plaque level in the APPPS1-21 mouse model compared to SAD therefore underlines the relevance of at least this tg mouse model in a human Aβ plaque context.
Our approach identified significantly upregulated proteins associated to two general processes; metabolism and endocytosis. In contrast, the downregulated proteins were associated to neuronal cytoskeleton. Glial cells have frequently been observed to engage and surround plaques in post mortem tissue analysis of AD patients and tg mouse models of AD [52–54]. Astrocytes are hypothesized to release mediators that affect the inflammatory response, for example by production of proteases that induce the Aβ degradation and Apoe that lipidate Aβ which increase the clearance rate [55]. Similarly, activated microglia can migrate to Aβ plaques and initiate phagocytosis and/or release proteases [55]. We reason that the relationship between the up- and downregulated proteins identified in our study could be reflected in the plaque engagement by glial cells and degeneration of neuronal processes (neuritic plaques). The consistent identification of astrocytic markers such as GFAP, Glul, and Slc1a2 in the plaque extracts further supports this reasoning. Engaging glial cells are considered to exhibit increased metabolism that is associated with the degradation process. For example, Pgk1 and Eno1 have previously been reported as enzymes involved in glycolysis and gluconeogenesis but also linked to AD which supports their identification [56, 57].
A recent study used a similar approach to compare the proteomic similarity of Aβ plaques in AD and APPPS1-21 mice [58]. The study identified 4,092 ortholog pairs when comparing AD and APPPS1-21 brain samples indicating high degree of translatability between mouse and human brain proteome. The study also compared differential pathways (proteomic patterns) between plaque and adjacent control regions in both APPPS1-21 mice and AD [58]. They identified pathways that were either similar, such as proteins associated with cholesterol efflux and humoral immune responses, or different; such as proteins associated with AMPA receptor activity and glutathione metabolic processes. Among the most regulated pathways in APPPS1-21 mice were proteins associated with the TYRO protein tyrosine kinase-binding protein (TYROBP) network and microglia pathogen phagocytosis [58]. TYROBP (also known as DAP12) is a sensing and phagocytosis-related, microglia-enriched, phosphotyrosine phosphoprotein that forms functional complexes with TREM2 and the complement receptor CR3, among others [59, 60]. We observed a high number of proteins associated with glial sensing, indicating similarity to the findings of Xiong and colleagues, who also reported that proteins associated with complement activation and inflammatory response were highly regulated in AD. We observed cellular mediators of inflammation (e.g., peroxiredoxin-2 and 14-3-3 protein) among the consistent plaque-proteins in the APPPS1-21 mice, also supporting comparability with the findings from Xiong and colleagues. However, we did not observe proteins linked to the complement activation cascade associated with plaque extracts in our study. In contrast to Xiong and colleagues, we focused the microdissection to the cerebral cortex wherefore regional differences could explain the discrepancy on identifying proteins associated with the complement system.
Studies have indicated a notable involvement of different PTMs in the pathology of AD. These include for example, pyroglutamylation (pyro-glu), phosphorylation, tyrosine nitration, and oxidation of Aβ peptides [14]. However, it remains to be elucidated whether co-localizing proteins are subjected to similar degree of PTMs as Aβ plaques. We are the first to describe unique modification associated to proteins co-localizing with Aβ plaques. We identified 426 common modifications on proteins co-localizing with Aβ plaques and adjacent control extracts. Here, no difference could discriminate the two subgroups in the PCA plot of the LFQ intensities associated to the modifications.
A clear difference could be observed in the number of unique modifications in the two subgroups. The number of unique modifications on proteins co-localizing with Aβ plaque extracts were higher compared to adjacent control regions (193 versus 117), indicating that proteins co-localized with aggregated Aβ are more susceptible for modifications.
On a technical note, to limit false-positive identification of unique PTMs by random stochastic selection of ions during MS/MS, we applied the following validity criteria: 1) The modification was present in more than 50% of samples within a group, and 2) had a minimum a-score of 100 determined by PEAKS X Studio. The software uses a feature-based peptide identification method based on MS1 detection with a femtomole to sub femtomole detection range on the Q-Exactive HF-X instrument that support a true-positive identification.
We observed that oxidation and deamidation were the most abundant modifications on proteins within the Aβ plaque and adjacent control extracts in the APPPS1-21 mice. Oxidative stress and generation of free radicals have been observed in AD brains. These have been hypothesized to originate from Aβ peptides [61]. Our findings suggest that oxidative stress is also occurring in the brain of APPPS1-21 mice, reflected by the increase levels of oxidized proteins. Increased levels of deamidation have similarly been linked to dementia [62]. Our detection of deamidated proteins in APPPS1-21 mice is therefore also in line with a previous study using human brain tissue. On the other hand, deamidation and oxidation of peptides have been reported to be generated during sample preparation and should therefore be interpreted with caution [63, 64]. We therefore handled samples from all subgroups in a parallel approach, suggesting a similar degree of experimentally-generated oxidation and deamidation across the subgroups.
Citrullination was observed as a unique modification in proteins co-localizing with both Aβ plaques and adjacent control regions. This protein modification was identified as the third most abundant modification after oxidation and deamidation in the adjacent control group. Accumulation of citrullinated proteins have been observed in AD patients [65]. Citrullination has also been associated with autoimmune diseases such as rheumatoid arthritis and multiple sclerosis [66, 67], and detection of citrullination in the APPPS1-21 brains could indicate widespread inflammatory response in this mouse model.
The number of unique pyro-glu modification of proteins co-localizing with Aβ plaques was found to be much higher compared to the number in adjacent control regions (25 versus 6) in APPPS1-21 mice. We observed the unique pyro-glu modification on the following plaque co-localizing proteins: Acly, Atp2b3, Atp2b4, Clu, Dync1h1, Dlat, Dpysl2, Etfa, Itpr1, Krt6a, Krt6b, Atp2b1, Rab1A, Rab5c, Sept11, Atp1a1, Atp1a2, Atp1a3, Sptan1, Thy1, Nsf, and Atp6v1c1. Of these, Rab1A, Clu, and Sptan1 were among the significantly regulated proteins identified in our LFQ analysis in MaxQuant and Perseus. The pyro-glu modification of Aβ peptides has been observed to localize within the core of Aβ plaques and hypothesized to serve as seeds for new plaque deposition [68]. The pyro-glu modification is generated by cleavage of the N-terminus followed by a conversion catalyzed by the enzyme glutaminyl cyclase. The levels of glutaminyl cyclase have been found to correlate with pyro-glu modified Aβ peptides in AD patients [69]. Based on our findings from the APPPS1-21 mice, we suggest that glutaminyl cyclase activity is also increased in Aβ plaque regions in the tg mice, leading to increased pyro-glu modifications on co-localizing proteins.
Several immunotherapies targeting Aβ proteoforms have failed in the clinic, for many reasons, one of which could be a poor understanding of the molecular nature of aberrant neurotoxic Aβ species, processes, and PTMs driving pathogenesis. We identified unique modifications on App, Mapk3, Itpr1, Gnaq, and Cyc1 from Aβ plaque regions in the APPPS1-21 mice. These five proteins are reported in Kyoto Encyclopedia of Genes and Genomes (KEGG) AD pathway, highlighting their relevance in the human disease. Furthermore, proteins associated with endocytosis, such as CD47, Dnm1, and Rab1 were also identified to express unique modifications in Aβ plaque regions in the APPPS1-21 mice. As PTMs can affect the biological function of proteins, these modifications could have a potential impact and importance of the AD pathology. Therefore, findings from this study could increase our knowledge of the complex plaque-brain interaction and how PTMs might be involved in the process. Our study may support efforts to design and develop molecules that could target PTMs in AD.
In conclusion, this study has optimized a protocol for extracting and investigating the Aβ plaque proteome composition by LC-MS/MS. We extracted Aβ plaques from the APPPS1-21 mouse model and identified several proteins that support the current biological view of glial cell engagement in response to Aβ plaque deposition. The data further implicated that the level of proteins associated to degradation were increased within Aβ plaque extracts compared to the surrounding tissue in APPPS1-21 mice. Furthermore, we identified unique PTM differences between Aβ plaque and adjacent control extracts in APPPS1-21 mice which also could have relevant impact in AD pathogenesis.
ACKNOWLEDGMENTS
This work was partly funded by the Innovation Fund Denmark (IFD) under File No. 5189-00044B and the Lundbeck Foundation grant (R247-2017-239). The Danish National Mass Spectrometry Platform for Functional Proteomics (PRO-MS; grant no. 5072-00007B); the Obelske Family Foundation and the Svend Andersen Foundation are acknowledged for grants to the analytical platform enabling parts of this study.
Authors’ disclosures available online (https://www.j-alz.com/manuscript-disclosures/19-0652r1).
