Integrated proteomics and network analysis identifies protein hubs and network alterations in Alzheimer’s disease
Acta Neuropathologica Communications volume 6, Article number: 19 (2018)
Although the genetic causes for several rare, familial forms of Alzheimer’s disease (AD) have been identified, the etiology of the sporadic form of AD remains unclear. Here, we report a systems-level study of disease-associated proteome changes in human frontal cortex of sporadic AD patients using an integrated approach that combines mass spectrometry-based quantitative proteomics, differential expression analysis, and co-expression network analysis. Our analyses of 16 human brain tissues from AD patients and age-matched controls showed organization of the cortical proteome into a network of 24 biologically meaningful modules of co-expressed proteins. Of these, 5 modules are positively correlated to AD phenotypes with hub proteins that are up-regulated in AD, and 6 modules are negatively correlated to AD phenotypes with hub proteins that are down-regulated in AD. Our study generated a molecular blueprint of altered protein networks in AD brain and uncovered the dysregulation of multiple pathways and processes in AD brain, including altered proteostasis, RNA homeostasis, immune response, neuroinflammation, synaptic transmission, vesicular transport, cell signaling, cellular metabolism, lipid homeostasis, mitochondrial dynamics and function, cytoskeleton organization, and myelin-axon interactions. Our findings provide new insights into AD pathogenesis and suggest novel candidates for future diagnostic and therapeutic development.
Alzheimer’s disease (AD) is the most common neurodegenerative disorder and the leading cause of dementia in the elderly [53, 70]. Neuropathologically, AD is characterized by the presence of amyloid plaques and neurofibrillary tangles in the brain. The vast majority (95%) of AD cases are sporadic, and the remaining 5% are familial AD . The causative genetic defects for several familial forms of AD have been identified, however, the etiology of sporadic AD remains unknown. The lack of effective means to prevent or treat AD and the failure of recent clinical trials [23, 36, 74] emphasize the need for better understanding AD pathogenic mechanisms to find novel targets for AD therapeutic intervention.
Human postmortem AD brain tissues provide a unique and valuable resource for discovery research to identify specific molecular abnormalities and disease processes associated with sporadic AD. We and others have previously used two-dimensional gel electrophoresis (2-DE)-based proteomics to study differential protein expression in AD versus control brains [17,18,19, 41]. Although these studies have found some proteins with altered expression in AD [17,18,19, 41], a limitation of 2-DE proteomics is its relatively low resolution, which limits the number of proteins that can be identified using this approach [6, 30, 79]. Recent advances in high-resolution, high-mass-accuracy mass spectrometry-based proteomics technologies provide powerful, new tools for in-depth profiling and quantitative analysis of protein expression in complex biological samples such as human brain tissues [21, 69].
With the advanced proteomics technologies enabling simultaneous, quantitative measurement of expression profiles for thousands of proteins, how to analyze such large proteomic data sets at the systems level becomes a major challenge. Weighted gene co-expression network analysis (WGCNA) is a systems biology approach originally developed for analysis of high-throughput transcriptomic data to provide an unbiased systems-level organization of the transcriptome into a network of biologically meaningful modules of co-expressed genes [45, 62, 92]. The use of WGCNA in studying transcriptome changes in a number of human diseases has led to the identification of disease-associated network modules and hub genes, which are the most highly connected genes that are key determinants of module function and represent important molecular targets for understanding and treating diseases [12, 27, 33, 46, 54, 82, 88]. Recent studies have begun to show that WGCNA can also be used in analyzing large proteomic data sets to gain systems-level insights into disease-associated proteome changes [37, 71, 80, 93].
In the present study, we performed large-scale, unbiased proteomic analyses of human AD and control frontal cortex tissues to determine disease-associated brain proteome changes by using a liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based, label-free quantitative proteomic approach. In addition to differential expression analysis to identify brain proteins with significantly altered abundance in AD, we performed WGCNA-based systems-level analysis of our entire proteomic data set and identified a network of disease-associated protein modules and intra-modular hub proteins in AD brain. Our study reveals dysregulation of multiple pathways and processes in AD brain and provides novel insights into the pathogenic mechanisms of sporadic AD.
Materials and methods
Human brain tissues
Postmortem frontal cortex tissues from neuropathologically confirmed AD cases and age-matched control subjects were obtained from Emory Center for Neurodegenerative Disease Brain Bank. Amyloid plaque pathology was assessed using the Consortium to Establish a Registry for Alzheimer’s Disease (CERAD) protocol for neuritic plaque scoring , and neurofibrillary tangle pathology was assessed using the Braak staging system . All AD cases meet the criteria of high level of AD neuropathological change based on the ABC scores according to the National Institute on Aging-Alzheimer’s Association guidelines for the neuropathological assessment of Alzheimer’s disease . ApoE genotypes were determined as previously described . Control subjects had no known history of neurological disease and showed no significant neurodegenerative changes at autopsy. Clinical and neuropathological data of all cases, including age, gender, disease status, age at onset, amyloid plaque pathology, neurofibrillary tangle pathology, ApoE genotype, and postmortem interval, are provided in Additional file 1: Table S1. Power analysis showed that the sample size used in this study (the total number of subjects = 16; n = 8 in each AD or control group) has > 80% power at a two-sided Type I error rate of 5% to detect effect size of > 1.6.
Brain tissue homogenization and protein extraction
Approximately 25 mg of human frontal cortex tissue from each AD or control case was homogenized as described  in 150 μl of lysis buffer containing 4% SDS, 100 mM DTT, and 100 mM Tris–HCl, pH 7.6, followed by incubation at 95 °C for 5 min. After cooling to room temperature, the homogenate was centrifuged at 16,000 x g for 5 min to obtain supernatant containing extracted proteins. Because the presence of SDS efficiently inactivates protease activity , no protease inhibitors were included during the brain tissue homogenization and protein extraction process. Protein concentrations of brain protein extracts were measured by UV spectrometry at 280 nm with NanoDrop spectrophotometer (ThermoFisher) using an extinction coefficient of 1.1 for 0.1% (g/L) solution .
Filter-aided sample preparation (FASP)
Human brain protein extracts were processed by using the FASP protocol as described . Briefly, 30 μl of each protein extract was mixed with 200 μl of 8 M urea in 100 mM Tris-HCl, pH 8.5 (UA solution), and the mixture was transferred into a Microcon 30-kDa centrifugal filter unit (MRCF0R030, Merck) and centrifuged at 14,000 x g for 15 min. Cysteine residues were alkylated by adding 100 μl of UA solution containing 50 mM iodoacetamide to the filter unit and incubation in darkness for 30 min at room temperature. After centrifugation at 14,000 x g for 10 min, 100 μl of UA solution was added to the filter unit and centrifuged again. This UA washing step was repeated twice, and the filter unit was then washed with 100 μl of 50 mM NH4HCO3 two times. Next, protein digestion was carried out by adding 40 μl of 50 mM NH4HCO3 solution containing sequencing-grade trypsin (enzyme to protein ratio 1:100) in the filter unit and incubation at 37 °C for 12 h. Digested peptides were eluted by adding 100 μl of 50 mM NH4HCO3 and collected by centrifugation at 14,000 x g for 10 min as a filtrate, and this step was repeated five times. The collected peptides were further purified by using a self-packed C18 ZipTip micro-column. The final concentration of peptides was determined by UV-spectrometry using an extinction coefficient of 1.1 for 0.1% (g/L) solution at 280 nm . All peptides were dried under vacuum at room temperature.
Liquid chromatography-tandem mass spectrometry
LC-MS/MS proteomic analyses were performed using the LTQ-Orbitrap Elite mass spectrometer (ThermoFisher) equipped with an EASY-Spray source and a nano-LC UltiMate 3000 high-performance liquid chromatography system (ThermoFisher). Human brain-derived peptides (2 μg) from each sample were separated by online reversed phase (RP)-HPLC fractionation on an EASY-Spray PepMap C18 column (length, 50 cm; particle size, 2 μm; pore size, 100 Å; ThermoFisher), using a 240-min gradient from 2% to 50% solvent B at a flow rate of 300 nL/min (mobile phase A, 1.95% acetonitrile, 97.95% H2O, 0.1% formic acid; mobile phase B, 79.95% acetonitrile, 19.95% H2O, 0.1% formic acid). A full-scan survey MS experiment (m/z range from 375 to 1600; automatic gain control target, 1,000,000 ions; resolution at 400 m/z, 60,000; maximum ion accumulation time, 50 ms) was performed using the Orbitrap mass analyzer. The ten most intense ions were selected and fragmented in the LTQ mass spectrometer (automatic gain control target value, 10,000) via collision-induced dissociation (CID) with maximum ion accumulation time of 100 ms. Raw data were analyzed by using Proteome Discoverer 1.4 (ThermoFisher) to search against the human Uniprot TrEMBL database (2016_02 Release, 20,198 reviewed entries). The modifications were set as follows: static modification of carbamidomethyl (Cys, + 57.0214 Da); dynamic modification of deamination (Asn, + 0.9840 Da), oxidation (Met, + 15.9949 Da), and acetylation (Lys, + 42.0106 Da). Trypsin was selected as the proteolytic enzyme, and up to two missed cleavages were allowed. The mass tolerance was set to 20 ppm for the precursor ions and 0.5 Da for the fragment ions. The false discovery rate (FDR) for peptide and protein identification was set to 1%.
Label-free protein quantification
Label-free protein quantitative analysis was performed by using Proteome Discoverer 1.4 to quantify precursor ion peak area (i.e., area under the curve), which is linearly proportional to protein abundance . A limitation of “shotgun” label-free quantitative proteomics is that protein identification or abundance data can be missing in some samples . Therefore, we restricted quantitative analysis to the proteins with complete data in all 16 brain samples, excluding proteins with missing data in any sample. In each sample, relative protein abundance for each protein was determined by normalizing the peak area of the protein to the total peak area of all proteins in the sample as described . To account for technical variability present in filter-aided sample preparation and LC-MS/MS analyses, each protein extract was spiked with bovine alpha-2-HS-glycoprotein (fetuin) at 0.1% (μg/μg total protein) as an internal control. The relative protein abundances determined by normalizing each protein peak area to that of the spike-in fetuin protein were similar to the abundances determined by normalizing to the total protein peak area, confirming the validity of the ‘Total Protein Approach’-based protein quantification analysis . The technical variation of the FASP sample processing and LC-MS/MS quantification system, estimated from the relative abundances of the spike-in fetuin protein after normalization of its peak area to the total protein peak area in eight control brain samples, had a coefficient of variation (CV) of 6%.
Differential expression analysis
Differentially expressed proteins in AD versus control were identified by using unpaired two-tailed Student’s t test with the thresholds of ±1.3-fold change over the control (i.e., AD/control ratio > 1.3 or < 0.77) and a P value < 0.05. The q values were calculated by using the q value R package  to correct for multiple comparisons and estimate the false discovery rates . Significantly altered proteins in AD with corresponding P values and q values are provided in Additional file 2: Table S2.
Hierarchical clustering analysis
Unsupervised hierarchical clustering of individual clinical cases and the identified differentially expressed proteins was performed based on their relative protein abundances in each samples by using Heatmapper online tool with an average linkage clustering and Kendall’s tau distance measurement method . Protein expression heat map with dendrograms showing clustering results were generated and visualized by the Heatmapper.
Protein co-expression network analysis
Protein co-expression network analysis was performed with the R package WGCNA as described  using the entire proteomic data set of all identified proteins with no missing values. Briefly, a correlation matrix for all pair-wise correlations of proteins across all samples was generated and then transformed into a matrix of connection strengths, i.e., a weighted adjacency matrix, as described [45, 92] with a soft threshold power β = 16. The connection strengths were then used to calculate topological overlap (TO), a robust, pairwise measure which indicates two proteins’ similarity based on their co-expression relationships with all other proteins in the network . Proteins were hierarchically clustered using 1 − TO as the distance measure to generate a cluster dendrogram, and modules of proteins with similar co-expression relationships were identified by using a dynamic tree-cutting algorithm  with the following parameters: minimal module size = 23, deepSplit = 4, and merge cut height = 0.07. For each module, a module eigenprotein was defined as the first principal component of the module which is a weighted summary of protein expression in the module and explains the maximal possible variability for all proteins within the module . Module membership (kME) was determined by calculating Pearson correlation between each protein and each module eigenprotein and the corresponding P-values . Proteins were (re)assigned to the module for which they had the highest module membership with a reassignment threshold of P < 0.05. Module-trait relationships were determined by using the WGCNA package  to calculate the biweight midcorrelations between each module eigenprotein and a clinical or neuropathological trait and the corresponding P-values. Module networks were graphically depicted by using the igraph package in R .
Gene ontology enrichment analysis and functional annotation of modules and proteins
Gene ontology (GO) enrichment analysis of the generated datasets of differentially expressed proteins and WGCNA module proteins was performed using MetaCore bioinformatics software (Version 6.29, build 68,613; https://portal.genego.com/). The total list of all proteins identified in human frontal cortex samples was used as the background. The hypergeometric test after the Benjamini-Hochberg false discovery rate (FDR) correction was used to assess statistical significance. Enriched GO terms with FDR-corrected P < 0.05 were considered statistically significant. In addition to the use of functional annotation tools, we also searched PubMed manually to gain insights into the functions of the identified differentially expressed proteins and WGCNA module proteins.
Western blot analysis
Human frontal cortex tissues from individual AD or control cases were homogenized in SDS lysis buffer, and protein extracts were subjected to SDS-PAGE. The proteins were then transferred onto PVDF membranes (EMD Millipore) and probed with anti-Smac antibody (1:1000 dilution; Cell Signaling Technology), anti-STK39 antibody (1:1000 dilution; Cell Signaling Technology), or anti-β-actin antibody (1:5000 dilution; EMD Millipore) followed by horseradish-peroxidase-conjugated secondary antibodies (Jackson ImmunoResearch Laboratories) and visualization using enhanced chemiluminescence as described previously . The expression levels of each protein were quantified by measuring protein intensities on immunoblot images using the Image J software (National Institutes of Health) and normalized to the corresponding level of β-actin in each sample. The normalized protein abundances across AD and control cases were compared by using unpaired two-tailed Student’s t test, and P < 0.05 was considered statistically significant.
Analysis of AD-associated proteome changes by quantitative proteomics
To investigate brain proteome alterations associated with sporadic AD, we analyzed brain samples from eight clinically and neuropathologically characterized AD patients and eight age-matched control subjects (Additional file 1: Table S1). Proteins were extracted from the dorsolateral prefrontal cortex tissues of these individuals by using the detergent sodium dodecyl sulfate (SDS) because it is the most effective reagent for solubilizing tissues and cells to achieve complete extraction of proteins . We used a recently developed, filter-aided sample preparation (FASP) method [85, 87] for detergent removal and protein digestion to obtain high-purity peptides from the brain samples. Subsequent LC-MS/MS proteomic analysis using the high-resolution high-mass-accuracy LTQ-Orbitrap Elite mass spectrometer identified a total of 39,819 distinct peptides, corresponding to 6679 unique proteins. Due to stochastic nature of “shotgun” label-free quantitative proteomics, protein identification or abundance data are sometimes missing in certain samples . The proteins with missing data in any sample were excluded in our analysis, resulting in the final quantification of 1968 proteins with complete data across all 16 brain samples from AD and control cases (Additional file 2: Table S2).
Differential expression analysis identifies proteins with altered abundance in AD
We performed differential expression analysis of quantitative proteomics data using the thresholds of ±1.3-fold change in AD over the control (P < 0.05) and identified 487 differentially expressed proteins (262 up-regulated proteins and 225 down-regulated proteins) in AD at FDR < 0.11 (Fig. 1a and Additional file 2: Table S2). Unsupervised hierarchical clustering analysis based on the protein abundances in the 16 individual brain samples showed that the identified differentially expressed proteins can serve as a proteomic signature for distinguishing AD versus control cases (Fig. 1b). The heat map illustrated an overall reproducibility as well as individual heterogeneity of protein expression profiles among different subjects within the AD or control group (Fig. 1b).
The list of the identified dysregulated proteins in AD (Additional file 2: Table S2) includes a number of proteins that have been previously shown by our group and others to be differentially expressed in AD brain, such as DJ-1, APOE, clusterin (CLU), and UCH-L1 [1, 17, 19, 55]. In addition, our proteomic analysis also identified 322 novel proteins that have not been previously reported as differentially expressed in AD, such as serine/threonine protein kinase 39 (STK39) and DIABLO/Smac (Additional file 3: Table S3). To validate our proteomic results, we performed Western blot analysis of STK39 and Smac expression in AD and control brains (Fig. 1c-g). We found that, in accordance with the proteomic data (Additional file 2: Table S2), STK39 protein level was significantly decreased in AD versus control (Fig. 1c, d). STK39 is an important kinase that has been associated with hypertension, Parkinson’s disease, and autism [50, 67, 84]. Our results indicate, for the first time, a link between STK39 and AD. In addition, our Western blot analysis also validated Smac, a key regulator of apoptosis , as an up-regulated protein in AD brain (Fig. 1e-f). Together, these results provide support for the robustness of our label-free quantitative proteomic analysis.
Next, we performed gene ontology (GO) enrichment analysis of the identified differentially expressed proteins to gain insights into the cellular functions and biological processes that are affected in AD brain (Fig. 2; Additional file 4: Table S4). We found that down-regulated proteins in AD were significantly enriched with GO categories linking to ion transport, mitochondrial function, synaptic transmission, myelin sheath, cell-cell adhesion, cytoskeleton organization, and endocytosis, whereas up-regulated proteins in AD were overrepresented with GO terms associated with metabolic process, immune response, cell-cell adhesion, exocytosis, vesicle-mediated transport, response to oxidative stress, translation, and regulation of apoptotic signaling (Fig. 2; Additional file 4: Table S4).
Co-expression network analysis uncovers AD-associated protein network alterations
To gain systems-level insights into the brain proteome changes in AD, we performed protein co-expression network analysis by using WGCNA, a data-driven network approach which uses pairwise correlation relationships of proteins and their topological overlap to organize the proteome into a network of biologically meaningful modules of co-expressed proteins [45, 90, 92]. We applied WGCNA to our entire proteomic data set of all proteins with no missing values (n = 1968 proteins) and constructed a protein co-expression network from protein expression profiles across all AD and control samples. Our WGCNA analysis identified 24 network modules of strongly co-expressed proteins (Fig. 3a; Additional file 5: Table S5). These modules, color coded according to the convention of WGCNA [45, 92], were labeled M1 to M24 based on the module size, ranging from the largest (M1: 223 proteins) to the smallest (M24: 30 proteins) (Fig. 3b). We found that several modules were significantly enriched for brain-specific GO categories, including mitochondria and synaptic transmission (M4), neuron part (M6), nervous system development (M7), myelin sheath and axonal organization (M12), and action potential (M24), whereas other modules were associated with GO categories linked to discrete cellular structures and functions, such as proteostasis and RNA homeostasis (M1), metabolism and lipid homeostasis (M2), cell morphogenesis (M3), mitochondria and cell adhesion (M5), hormone activity (M8), membrane assembly (M9), ion and protein transport (M10), signaling and cytoskeleton regulation (M11), hydrolase activity (M13), ribosome (M14), immune response (M15), inflammatory response (M16), and extracellular region (M17) (Fig. 3b; Additional file 6: Table S6).
To identify disease-relevant modules associated with AD phenotypic traits, we assessed the module-trait relationships by determining the biweight midcorrelations between each module eigenprotein (the module representative which summarizes protein expression profiles in the module ) and various disease-related traits or sample variables (Fig. 4). We identified 11 modules that were significantly correlated with AD status, amyloid plaque pathology (frontal cortex neuritic plaque frequency), and/or neurofibrillary tangle pathology (Braak stage), including 5 positive correlated modules (M1, M2, M15, M16, and M19) and 6 negatively correlated modules (M4, M5, M10, M11, M12, and M13). None of the modules showed significant correlation with age, gender, ApoE genotype, or postmortem interval (Fig. 4), confirming that the identified AD-correlated modules are not due to any of the potential confounding factors. Our analysis showed that most of the positively correlated modules (M1, M2, M15, and M16) had significantly increased module expression levels in AD (Fig. 5b), whereas most of the negatively correlated modules (M4, M5, M10, M11, and M13) had significantly decreased module expression levels in AD (Fig. 5c).
We then assessed the inter-modular relationships by performing eigenprotein network analysis as described [32, 44] to construct a higher-order meta-network based on pairwise correlation relationships of module eigenproteins. The module eigenprotein meta-network revealed the inter-modular connectivity of 24 co-expression modules in brain proteome, showing a hierarchical organization of highly interconnected modules into meta-modules, i.e., groups of highly correlated module eigenproteins (Fig. 5a). Interestingly, the eigenproteins of all modules positively correlated with AD phenotypes (M1, M2, M15, M16, and M19) were clustered in a single meta-module (Fig. 4 and Fig. 5a), suggesting close relationships among the pathways and processes associated with these positively correlated modules. In addition, we identified a meta-module containing eigenproteins from 5 out of the 6 modules negatively correlated with AD phenotypes (M4, M5, M10, M11, and M13), indicating that the corresponding pathways and processes for these negatively correlated modules may also be related.
AD-associated network modules and hub proteins reveal multiple dysregulated pathways in AD brain
Highly connected hub nodes are central to a network’s architecture and function [2, 7], and intramodular hub proteins in disease-related WGCNA modules have emerged as key targets for biomarker and therapeutic development [12, 27, 33, 46, 54, 82, 88]. Intramodular hub proteins can be identified by using module membership (kME), a measure of intramodular connectivity [32, 46]. The top 10 highly connected hub proteins for each of the identified AD-related modules are shown in the center of network plots (Figs. 6 and 7). Unsupervised hierarchical clustering analysis based on the hub protein expression profiles showed that the identified top hub proteins serve as a molecular signature to differentiate AD and control cases (Fig. 8c). We found that the top hub proteins of the modules with positive correlation to AD phenotypes were often up-regulated in AD (Fig. 8a,c), whereas the top hub proteins of the negative correlated modules were often down-regulated in AD (Fig. 8b,c), consistent with the proposed role of hub proteins as key drivers of protein co-expression modules [32, 33]. We assessed the molecular and functional characteristics of each AD-associated module based on its top hub proteins and gene ontology enrichment analysis of module proteins to gain insights into the biological roles of AD-related modules (Additional file 6: Table S6).
Our analyses revealed that M1, the largest module positively correlated with AD phenotypes (Fig. 4), was significantly enriched with GO categories and hub proteins linked to pathways that control protein homeostasis, or “proteostasis” (Fig. 6 and Additional file 6: Table S6), including 11 protein translation machinery components (EIF2S2, EIF3A, EIF4A2, EIF4B, RPLP1, RPL3, RPL10, RPS6, RPS7, RPS14, and RPS17) with 40S ribosome subunit RPS7 as a top hub protein; 19 molecular chaperones and cochaperones (AHSA1, CDC37, BAG5, CANX, DNAJA2, DNAJA3, DNAJC13, FKBP4, ERO1A, GNB4, GANAB, PDIA3, PFDN5, PFDN6, TBCD, CCT4/TCP1-delta, CCT5/TCP1-epsilon, CCT6A/TCP1-zeta-1, and CLU) with Hsp70 cochaperone DNAJA2 as a top hub protein; and 11 proteasome complex components (PSMA2, PSMC1, PSMC2, PSMC4, PSMC6, PSMD1, PSMD12, PSMD13, PSMD14, RAD23A, and RAD23B) with 26S proteasome regulatory subunit PSMD1 as a top hub protein (Fig. 6 and Additional file 5: Table S5). The overrepresentation of the proteostasis machinery components in this AD-related module supports the involvement of dysregulated proteostasis in AD pathophysiology [43, 78, 89].
In addition to proteostasis, RNA homeostasis-related proteins and pathways were also enriched in the M1 module, as demonstrated by the presence of 12 ribonucleoproteins involved in RNA processing (HNRNPC, HNRNPK, HNRNPL, ALYREF, GCN1L1, SSB, NPM1, LUC7L3, TROVE2, EFTUD2, RUVBL1, and SNRPE) with heterogeneous nuclear ribonucleoprotein K (HNRNPK) as a top hub protein (Fig. 6 and Additional file 5: Table S5). Our finding of HNRNPK, a major RNA-binding protein which functions in regulation of transcription, RNA splicing, mRNA stability, and translation , as an up-regulated M1 hub protein in AD (Fig. 8) reveals a previously unrecognized role of HNRNPK in AD pathophysiology. Corroborating with our results, another related M1 module member, HNRNPC, has been reported to be increased in AD and promote APP translation [10, 66]. Additionally, we identified pro-apoptotic factors HTRA2 and AIFM1 as top hub proteins up-regulated in AD (Fig. 6, Additional file 2: Table S2, and Additional file 5: Table S5), indicating enhanced apoptotic signaling is another key feature of this module.
The relevance of the M1 module to AD is further strengthened by its association with APOE and CLU (Fig. 6), two well-established, genetic risk factors for sporadic AD . Our analyses showed that both APOE and CLU proteins were up-regulated in AD (Additional file 2: Table S2) and had high intramodular connectivity values (Additional file 5: Table S5), supporting their role as important determinants of M1 module functions. In addition, we found the fat mass and obesity-associated protein FTO, an AD risk factor which genetically interacts with APOE [38, 68], was the most highly connected hub protein of the M1 module (Fig. 6 and Additional file 5: Table S5). FTO, a demethylase which regulates 6-methyladenosine modifications of mRNAs, has also been linked to increased risk for obesity and type 2 diabetes . Another M1 hub protein, SORBS1 (Fig. 6), which functions in insulin signaling, has also been associated with obesity and type 2 diabetes . The finding of obesity and diabetes-associated FTO and SORBS1 as top hub proteins in AD-related M1 module is consistent with increasing evidence indicating the presence of shared pathways in the pathogenesis of AD, obesity, and diabetes .
M2, a 152-member module with positive correlation to AD phenotypes (Fig. 4), was highly enriched with GO categories, enzymes, and hub proteins linked to metabolic processes and pathways (Fig. 6 and Additional file 6: Table S6). The most prominent feature of this module is the presence of over 40 proteins that function in the carboxylic acid metabolism with serine racemase (SRR) and enolase 1 (ENO1) as top hub proteins (Fig. 6 and Additional file 5: Table S5). SRR, an enzyme for catalyzing the conversion of L-serine to D-serine (an essential co-agonist of the NMDA receptor) , was up-regulated by more than two folds in AD (Additional file 2: Table S2), which may lead to over-activated NMDA receptors, thereby contributing to AD pathophysiology. The M2 module was also highly enriched with proteins involved in the unsaturated fatty acid metabolic process (ACAA1, ACOX1, EPHX2, HSD17B4, LTA4H, PTGDS, PTGR1, PTGR2, GSTM2, GSTM3, GSTP1, and MIF), highlighting a link between dysregulated unsaturated fatty acid metabolism and AD pathophysiology. Furthermore, the M2 module was also significantly enriched with regulators of lipid metabolism (AGK, ACAA2, ALDH3A2, ANXA1, ANXA2, ANXA4, ANXA5, ASAH1, APPL2, DBI, ESYT1, GM2A, HADHA, INPP1, PAFAH1B3, ERLIN2, SLC44A2, PCYT2, PLCD3, and PRDX6) with annexin A5 (ANXA5) as a top hub protein (Fig. 6 and Additional file 5: Table S5). These findings provide new insights into the molecular basis of dysregulated lipid homeostasis in AD brain [26, 60].
The identified top M2 hub proteins also include all three members of the ezrin-radixin-moesin (ERM) family, ezrin (EZR), radixin (RDX), and moesin (MSN), which were up-regulated in AD (Fig. 6 and Additional file 2: Table S2), suggesting a role of ERM proteins in AD. The ERM proteins are FERM (4.1 protein, ezrin, radixin, moesin) domain-containing proteins that function as plasma membrane–cytoskeleton linkers to regulate membrane dynamics, cell adhesion, migration, signal transduction, and immune response . Interestingly, another FERM domain-containing protein, FERMT2, was also identified as an up-regulated M2 hub protein with high intramodular connectivity (Fig. 6, Additional file 2: Table S2, and Additional file 5: Table S5). Our finding, together with the reports of FERMT2 as a genetic risk factor for AD  and a modulator of APP metabolism and tau neurotoxicity [16, 72], supports the involvement of FERMT2 in AD pathogenesis.
M15, a 57-member module positively correlated with AD phenotypes (Fig. 4), was significantly enriched with GO terms and proteins linked to immune response (ALCAM, ALAD, GAPDH, CYB5R3, DDX3X, CAPN1, PPIA, PYGB, EIF2AK2, CAB39, TTR, PDAP1, HIST1H2BK, QARS, VAPA, and PNP) with GAPDH, PPIA, CYB5R3, and PYGB as top hub proteins (Fig. 6 and Additional file 6: Table S6). Our finding of PPIA, CYB5R3, and PYGB, which are associated with neutrophil activation in immune response , as up-regulated M15 hub proteins (Fig. 6 and Additional file 2: Table S2) supports a role of neutrophil-dependent immune response in AD pathophysiology . The enrichment of several aminoacyl-tRNA synthetases for protein translation (SARS, QARS, and NARS) in this module (Fig. 6) is in agreement with AD-associated, protein translation alteration identified in the M1 module.
M16, a 55-member, AD-positively correlated module (Fig. 4), was significantly enriched with GO categories and proteins linked to inflammatory response (C3, FHL1, A2M, CD44, FN1, HP, SERPINA1, and SERPINA3) with complement C3 and alpha-2-macroglobulin (A2M) as top hub proteins (Fig. 6 and Additional file 6: Table S6). Our finding of C3 and A2M, two key components of inflammatory response [42, 63], as up-regulated hub proteins in AD brain (Fig. 6 and Additional file 2: Table S2) supports their potential as candidate AD biomarkers  and the link between neuroinflammation and AD pathogenesis . Another key molecular feature of M16 module is significant overrepresentation of extracellular matrix proteins (COL6A1, COL6A2, COL6A3, COL18A1, LAMA5, FLNA, and FLNB) as hub proteins (Fig. 6), providing evidence for the involvement of extracellular matrix dysfunction in AD .
M19, a 48-member module with positive correlation to AD phenotypes (Fig. 4), was highly enriched with GO terms and proteins linked to small GTPase-mediated trafficking and signaling (RAB1A, RAB1B, RAB3C, RAB4A, RAB4B, RAB6B, RAB8B, RAB10, RAB12, RAB14, RAB23, RAB35, ARF4, and ARF5) with Rab GTPases RAB1A, RAB1B, RAB4A, and RAB4B as top hub proteins (Fig. 6 and Additional file 6: Table S6). The enriched Rab and ARF GTPases function as key regulators of the following trafficking pathways: ER-to-Golgi transport (RAB1A and RAB1B), synaptic vesicle exocytosis and neurotransmitter release (RAB3C), endosome-to-plasma membrane recycling (RAB4A, RAB4B, RAB23, and RAB35), intra-Golgi traffic and endosome-to-Golgi transport (RAB6B), trans-Golgi network (TGN)-to-plasma membrane transport (RAB8B, RAB10, RAB12, RAB14, ARF4, and ARF5), Golgi-to-ER retrograde transport (ARF4 and ARF5), and autophagosome formation (RAB1A, RAB1B, RAB12, and RAB23) [24, 34, 75]. In addition, this module also contained endocytic trafficking regulators, NECAP1 and SORT1 (Fig. 6). The enrichment of the vesicular trafficking machinery components in the AD-correlated M19 module highlights the dysregulation of multiple trafficking pathways in AD brain.
M4, a 119-member module with negative correlation to AD phenotypes (Fig. 4), was highly enriched with GO categories and proteins linked to mitochondrial processes (Fig. 7, Additional file 5: Table S5, and Additional file 6: Table S6), including the mitochondrial contact site and cristae organizing system (MICOS) components (IMMT/MIC60 and CHCHD3/MIC19), mitochondrial import machinery components (TOMM70 and TIMM9), mitochondrial carrier system components (SLC25A3 and SLC25A12), mitochondrial inner membrane fusion GTPase OPA1, electron transport chain subunits (MT-ND5, NDUFA12, NDUFS5, NDUFS7, SDHA, COX4I1, COX5B, COX6C), and mitochondrial ATP synthase subunits (ATP5A1, ATP5B, and ATP5J). The identification of IMMT and SLC25A3 as down-regulated hub proteins in AD brain (Fig. 7 and Additional file 2: Table S2) reveals a previously unrecognized role of the MICOS system and mitochondrial carrier system in AD pathophysiology. The M4 module was also enriched with GO categories and proteins associated with synaptic structure and function (Fig. 7, Additional file 5: Table S5, and Additional file 6: Table S6), including key presynaptic proteins involved in synaptic vesicle trafficking and neurotransmitter release (SNAP25, STX1B, SYP, SV2B, VAMP1, RPH3A, and RAB3GAP2), glutamate receptor subunits (GRIA2/GluR2 and GRIA4/GluR4), GABAA receptor β1 subunit (GABRB1), and postsynaptic density proteins (SYNGAP1, DLG1/SAP97, DLG3/SAP102, and SRGAP3). The finding of these synaptic proteins in the AD-down-regulated M4 module is consistent with a loss of synaptic function in AD .
M5, a 106-member module negatively correlated with AD phenotypes (Fig. 4), was highly enriched with GO terms and proteins linked to oxidative phosphorylation (MT-ND3, NDUFA2, NDUFA7, NDUFAB1, NDUFB11, NDUFS1, NDUFS8, CYC1, UQCRC1, COX5A, COX6B1, and ATP5D) with MT-ND3 and UQCRC1 as top hub proteins, synaptic cell adhesion (CADM2/SynCAM2, CADM3/SynCAM3, NCAM1, LSAMP, NTM, OPCML) with CADM2/SynCAM2 and LSAMP as top hub proteins, synaptic vesicle exocytosis (STX1A, CPLX2, and SYT12), and signal transduction (PPP1CB, MARK2, PAK1, PRKCE, PRKCG, SLK, STK39, RHOA, RHOB, RHOC, RHOG, and CDC42) with PPP1CB, RHOA, and STK39 as top hub proteins (Fig. 7, Additional file 5: Table S5, and Additional file 6: Table S6). These findings further support the involvement of impaired mitochondrial and synaptic functions and dysregulated signaling in AD pathophysiology.
M10, a 62-member module with the most significant negative correlation to AD phenotypes (Fig. 4), was highly enriched with ion-transporting ATPases, such as Na+/K+ ATPase subunits (ATP1A2, ATP1B1, and ATP1B2) for establishing the electrochemical gradients of Na and K ions across the plasma membrane and the H+-transporting, vacuolar ATPase subunit ATP6V1G2 for lysosomal acidification (Fig. 7 and Additional file 6: Table S6), supporting a loss of brain cell ion homeostasis in AD pathogenesis [20, 81]. The M10 module was also significantly enriched with GO terms and proteins linked to transmembrane transport and vesicle-mediated transport (Fig. 7 and Additional file 6: Table S6), including mitochondrial protein import (MTX2), endocytosis (DNM2, DNM3), endosome-to-lysosome trafficking and synaptic vesicle biogenesis (AP3D1), ER-to-Golgi transport (BCAP31), intra-Golgi trafficking (USO1, NAPA, NAPG), TGN-to-plasma membrane transport (AP1M1), endosome-to-plasma membrane recycling (SNX27), and autophagy (GABARAPL2, RTN3, and UBQLN1). The enrichment of these transport machinery components in the AD-down-regulated M10 module indicates impairment of multiple transport pathways in AD brain.
M11, a 62-member module negatively correlated with AD phenotypes (Fig. 4), has heterotrimeric G-protein subunits (GNAI1, GNAI2, GNAI3, and GNAO1) and Src family of tyrosine kinases (FYN and YES1) as top hub proteins (Fig. 7), highlighting the involvement of altered intracellular signaling in AD pathophysiology. In addition, the M11 module was significantly enriched with GO terms and proteins linked to regulation of actin cytoskeleton (CRMP1, CRMP4/DPYSL3, ABI1, FSCN1, and WASF1) and microtubule cytoskeleton (TUBAL3, TBCB, NDRG1, NDRG2, SHTN1, and SNCG) (Fig. 7 and Additional file 6: Table S6), consistent with impaired actin and microtubule dynamics in AD brain [5, 25]. The association of the retromer complex components (SNX1 and SNX2) with the M11 module supports a link between retromer dysfunction and AD pathogenesis .
M12, a 61-member module with negative correlation to neurofibrillary tangle pathology but not amyloid plaque pathology (Fig. 4), is characterized by highly significant enrichment of GO terms and proteins linked to myelin sheath (CNP, MAG, MBP, OMG, PLP1, MOG, PMP2, CLDN11, and ERMN) and the organization of paranodal and juxtaparanodal regions of axon at the node of Ranvier (MAG, ERMN, CNTNAP1, CNTN2, and KCNA2) with CNP, MAG, OMG, and PLP1 as top hub proteins (Fig. 7 and Additional file 6: Table S6). These results, together with our finding of OMG and PLP1 as down-regulated hub proteins in AD (Fig. 7 and Additional file 2: Table S2), support the involvement of myelin degeneration, impaired myelin-axon interactions, and node of Ranvier dysfunction in AD pathogenesis . The M12 module was also significantly enriched with neurofilament proteins (NEFL, NEFM, and INA) and microtubule-binding proteins involved in the control of microtubule polymerization or stabilization (CRYAB, MAPRE1, DST, CRMP2/DPYSL2, CLASP2, and MAP1B) and axonal transport (DCTN1 and DCTN4), indicating an association of impaired neurofilament and microtubule functions with Tau aggregation in AD. Our finding of BIN1, the second most prevalent genetic risk factor for sporadic AD , as a member of the M12 module with negative correlation to neurofibrillary tangle pathology (Fig. 7) is consistent with recent evidence indicating that BIN1 negatively regulates the propagation of Tau pathology .
M13 is a 61-member module down-regulated in AD (Fig. 5) with negative correlation to AD phenotypes (Fig. 4). More than one-third of proteins in this module are associated with hydrolase activity, represented by top hub proteins such as deubiquitinating enzyme OTUB1, small GTPase RAC1, adenosylhomocysteinase-like proteins AHCYL1 and AHCYL2, protein tyrosine phosphatase PTPRZ1, and platelet-activating factor acetylhydrolase subunit PAFAH1B1 (Fig. 7 and Additional file 5: Table S5). In addition, the M13 module was also significantly enriched with GO terms and proteins linked to organization of the actin cytoskeleton (ARPC1A, CFL1, CFL2, CTTN, RAC1, PAFAH1B1, NF1, SPTAN1, and SPTBN2) and microtubule cytoskeleton (MAPRE3, KLC2, NIT2, TBCA, TUBA8, MAP7D1) (Fig. 7 and Additional file 6: Table S6), supporting the involvement of impaired actin and microtubule dynamics in AD pathophysiology [5, 25].
This study shows that integration of quantitative proteomics, differential expression analysis, and co-expression network analysis provides a useful approach for gaining systems-level insights into AD pathogenesis. A critical step in quantitative proteomic analysis is sample preparation, which is a key determinant of the quality of generated proteomic data set. Previous proteomic studies of AD brains used detergent-free, protein extraction with a chaotropic reagent such as urea, which is unable to completely solubilize brain tissue and extract all proteins [3, 56, 83, 87]. To overcome this limitation, we used the strong detergent SDS for complete solubilization and extraction of proteins followed by the filter-aided sample preparation procedure [85, 87] to obtain high-purity peptides for LC-MS/MS-based quantitative proteomic analysis. Our proteomic results support that the SDS-based, filter-aided sample preparation method is highly effective for achieving high proteome coverage and reliable measures of protein expression levels in human AD and control brain tissues.
Differential expression analysis, which compares expression levels for individual proteins between AD and control groups, is a commonly used method in proteomic studies to identify AD-associated protein changes [3, 56, 59]. Using this method, we have identified 487 differentially expressed proteins with significantly altered protein levels (> 1.3-fold change; P < 0.05) in AD versus control, including 262 up-regulated proteins and 225 down-regulated proteins involved in multiple biological processes. The identification of a wide spectrum of protein alterations is consistent with the multifactorial and complex etiology of AD. Our identified differentially expressed proteins include 322 novel proteins that are not previously known to be altered in AD (Additional file 3: Table S3), providing new insights into protein changes in AD brain. Due to the small sample size, this study is expected to have false positives as well as well as false negatives. The identified differentially expressed proteins in AD have a false discovery rate of < 11% based on the estimation by q values (Additional file 2: Table S2). Therefore, our findings will need to be confirmed in larger samples. The independent validation of altered expression in AD of two identified novel proteins, STK39 and DIABLO/Smac, by Western blot analysis highlights the robustness of our label-free quantitative proteomic analysis.
In contrast to differential expression analysis which determines expression changes of single proteins independently, co-expression network analysis relates proteins to each other using pairwise correlation relationships between protein expression profiles to illuminate higher-order molecular organization and define modules of co-expressed proteins that are functionally related and/or coordinately regulated [45, 90, 92]. Using this network analysis, we have identified 11 disease-associated, protein co-expression modules that are significantly correlated with AD phenotypes, including 5 positively correlated modules (M1, M2, M15, M16, and M19) and 6 negatively correlated modules (M4, M5, M10, M11, M12, and M13). The identified, AD-associated modules reveal a number of previously unrecognized co-expression relationships among proteins involved in distinct biological processes and provide a novel view of cellular mechanisms. For example, the M1 module shows that proteins controlling various processes of proteostasis (e.g., protein translation, protein folding, and proteasome-mediated degradation) and RNA homeostasis (e.g., RNA processing, transcription initiation, mRNA modification and stability) are highly connected at a co-expression level (Fig. 6, Additional file 5: Table S5, and Additional file 6: Table S6), indicating coordinate control or interactions among these different processes. The M12 module reveals a strong co-expression relationship linking myelin proteins, neurofilament proteins, and axonal proteins involved in microtubule-based transport (Fig. 7 and Additional file 6: Table S6), highlighting the glia-neuron interactions and coupling between myelin and axonal processes. The M19 module uncovers a previously unknown co-expression relationship connecting Rab GTPases, ARF proteins, and other key regulators of various intracellular membrane trafficking processes (Fig. 6 and Additional file 6: Table S6), suggesting co-regulation of multiple trafficking processes and their involvement in AD pathophysiology. Further studies of the identified protein co-expression relationships and their regulation will advance our knowledge of the cellular mechanisms governing coordinate control and concerted actions of various biological processes in health and Alzheimer’s disease.
Our proteomics-driven network analysis has generated a molecular blueprint of dysregulated protein networks in AD brain and has uncovered many new proteins and pathways in processes implicated in AD, including altered proteostasis, RNA homeostasis, immune response, neuroinflammation, synaptic transmission, vesicular transport, cell signaling, cellular metabolism, lipid homeostasis, mitochondrial dynamics and function, cytoskeleton organization, and myelin-axon interactions. The identified hub proteins of AD-associated protein network modules are particularly useful for biomarker and therapeutic development, as hub proteins are often key drivers of disease-related co-expression modules or key determinants of module function [12, 27, 33, 46, 54, 82, 88]. Our finding that the identified top hub proteins can serve as a molecular signature for differentiating AD and control cases (Fig. 8c) supports their potential as novel AD biomarkers. Furthermore, the hub proteins of AD-related modules uncovered in this study provide attractive drug targets for developing novel therapeutics to shift disease-specific changes of protein networks and cellular functions back to their normal range.
In summary, our integrated proteomics and network analysis provides a systems-level view of proteome changes in AD brain and uncovers disease-associated protein network alterations in AD. The identified AD-related network modules and their hub proteins generate new insights into the pathogenesis of sporadic AD. Our findings suggest new targets and biomarker candidates for AD diagnostic development and therapeutic intervention.
Akram A, Schmeidler J, Katsel P, Hof PR, Haroutunian V (2012) Association of ApoE and LRP mRNA levels with dementia and AD neuropathology. Neurobiol Aging 33: 628 e621–628 e614. Doi https://doi.org/10.1016/j.neurobiolaging.2011.04.010
Almaas E (2007) Biological impacts and context of network theory. J Exp Biol 210:1548–1558. https://doi.org/10.1242/jeb.003731
Andreev VP, Petyuk VA, Brewer HM, Karpievitch YV, Xie F, Clarke J, Camp D, Smith RD, Lieberman AP, Albin RL et al (2012) Label-free quantitative LC-MS proteomics of Alzheimer's disease and normally aged human brains. J Proteome Res 11:3053–3067. https://doi.org/10.1021/pr3001546
Babicki S, Arndt D, Marcu A, Liang Y, Grant JR, Maciejewski A, Wishart DS (2016) Heatmapper: web-enabled heat mapping for all. Nucleic Acids Res 44:W147–W153. https://doi.org/10.1093/nar/gkw419
Bamburg JR, Bernstein BW (2016) Actin dynamics and cofilin-actin rods in alzheimer disease. Cytoskeleton (Hoboken) 73:477–497. https://doi.org/10.1002/cm.21282
Bantscheff M, Schirle M, Sweetman G, Rick J, Kuster B (2007) Quantitative mass spectrometry in proteomics: a critical review. Anal Bioanal Chem 389:1017–1031. https://doi.org/10.1007/s00216-007-1486-6
Barabasi AL, Oltvai ZN (2004) Network biology: understanding the cell's functional organization. Nat Rev Genet 5:101–113. https://doi.org/10.1038/nrg1272
Bartzokis G (2011) Alzheimer's disease as homeostatic responses to age-related myelin breakdown. Neurobiol Aging 32:1341–1371. https://doi.org/10.1016/j.neurobiolaging.2009.08.007
Bomsztyk K, Denisenko O, Ostrowski J (2004) hnRNP K: one protein multiple processes. BioEssays 26:629–638. https://doi.org/10.1002/bies.20048
Borreca A, Gironi K, Amadoro G, Ammassari-Teule M (2016) Opposite dysregulation of fragile-X mental retardation protein and heteronuclear ribonucleoprotein C protein associates with enhanced APP translation in Alzheimer disease. Mol Neurobiol 53:3227–3234. https://doi.org/10.1007/s12035-015-9229-8
Braak H, Braak E (1991) Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol 82:239–259
Brohawn DG, O'Brien LC, Bennett JP Jr (2016) RNAseq analyses identify tumor necrosis factor-mediated inflammation as a major abnormality in ALS spinal cord. PLoS One 11:e0160520. https://doi.org/10.1371/journal.pone.0160520
Calafate S, Flavin W, Verstreken P, Moechars D (2016) Loss of Bin1 promotes the propagation of tau pathology. Cell Rep 17:931–940. https://doi.org/10.1016/j.celrep.2016.09.063
Calsolaro V, Edison P (2016) Neuroinflammation in Alzheimer's disease: current evidence and future directions. Alzheimers Dement 12:719–732. https://doi.org/10.1016/j.jalz.2016.02.010
Campanini B, Spyrakis F, Peracchi A, Mozzarelli A (2013) Serine racemase: a key player in neuron activity and in neuropathologies. Front Biosci (Landmark Ed) 18:1112–1128
Chapuis J, Flaig A, Grenier-Boley B, Eysert F, Pottiez V, Deloison G, Vandeputte A, Ayral AM, Mendes T, Desai S et al (2017) Genome-wide, high-content siRNA screening identifies the Alzheimer's genetic risk factor FERMT2 as a major modulator of APP metabolism. Acta Neuropathol 133:955–966. https://doi.org/10.1007/s00401-016-1652-z
Choi J, Levey AI, Weintraub ST, Rees HD, Gearing M, Chin LS, Li L (2004) Oxidative modifications and down-regulation of ubiquitin carboxyl-terminal hydrolase L1 associated with idiopathic Parkinson's and Alzheimer's diseases. J Biol Chem 279:13256–13264
Choi J, Rees HD, Weintraub ST, Levey AI, Chin LS, Li L (2005) Oxidative modifications and aggregation of cu,Zn-superoxide dismutase associated with Alzheimer and Parkinson diseases. J Biol Chem 280:11648–11655
Choi J, Sullards MC, Olzmann JA, Rees HD, Weintraub ST, Bostwick DE, Gearing M, Levey AI, Chin LS, Li L (2006) Oxidative damage of DJ-1 is linked to sporadic Parkinson and Alzheimer diseases. J Biol Chem 281:10816–10824
Colacurcio DJ, Nixon RA (2016) Disorders of lysosomal acidification-the emerging role of v-ATPase in aging and neurodegenerative disease. Ageing Res Rev 32:75–88. https://doi.org/10.1016/j.arr.2016.05.004
Craft GE, Chen A, Nairn AC (2013) Recent advances in quantitative neuroproteomics. Methods 61:186–218. https://doi.org/10.1016/j.ymeth.2013.04.008
Cuyvers E, Sleegers K (2016) Genetic variations underlying Alzheimer's disease: evidence from genome-wide association studies and beyond. Lancet Neurol 15:857–868. https://doi.org/10.1016/S1474-4422(16)00127-7
De Strooper B, Chavez Gutierrez L (2015) Learning by failing: ideas and concepts to tackle gamma-secretases in Alzheimer's disease and beyond. Annu Rev Pharmacol Toxicol 55:419–437. https://doi.org/10.1146/annurev-pharmtox-010814-124309
Donaldson JG, Jackson CL (2011) ARF family G proteins and their regulators: roles in membrane transport, development and disease. Nat Rev Mol Cell Biol 12:362–375. https://doi.org/10.1038/nrm3117
Dubey J, Ratnakaran N, Koushika SP (2015) Neurodegeneration and microtubule dynamics: death by a thousand cuts. Front Cell Neurosci 9:343. https://doi.org/10.3389/fncel.2015.00343
El Gaamouch F, Jing P, Xia J, Cai D (2016) Alzheimer's disease risk genes and lipid regulators. J Alzheimers Dis 53:15–29. https://doi.org/10.3233/JAD-160169
Fogel BL, Cho E, Wahnich A, Gao F, Becherel OJ, Wang X, Fike F, Chen L, Criscuolo C, De Michele G et al (2014) Mutation of senataxin alters disease-specific transcriptional networks in patients with ataxia with oculomotor apraxia type 2. Hum Mol genet 23: 4758-4769. Doi https://doi.org/10.1093/hmg/ddu190
Forner S, Baglietto-Vargas D, Martini AC, Trujillo-Estrada L, LaFerla FM (2017) Synaptic impairment in Alzheimer's disease: a dysregulated symphony. Trends Neurosci 40:347–357. https://doi.org/10.1016/j.tins.2017.04.002
Gearing M, Schneider JA, Rebeck GW, Hyman BT, Mirra SS (1995) Alzheimer's disease with and without coexisting Parkinson's disease changes: apolipoprotein E genotype and neuropathologic correlates. Neurology 45:1985–1990
Gygi SP, Corthals GL, Zhang Y, Rochon Y, Aebersold R (2000) Evaluation of two-dimensional gel electrophoresis-based proteome analysis technology. Proc Natl Acad Sci U S A 97:9390–9395. https://doi.org/10.1073/pnas.160270797
Heine SJ, Olive D, Gao JL, Murphy PM, Bukrinsky MI, Constant SL (2011) Cyclophilin a cooperates with MIP-2 to augment neutrophil migration. J Inflamm Res 4:93–104. https://doi.org/10.2147/JIR.S20733
Horvath S, Dong J (2008) Geometric interpretation of gene coexpression network analysis. PLoS Comput Biol 4:e1000117. https://doi.org/10.1371/journal.pcbi.1000117
Horvath S, Zhang B, Carlson M, Lu KV, Zhu S, Felciano RM, Laurance MF, Zhao W, Qi S, al CZ (2006) Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target. Proc Natl Acad Sci U S A 103:17402–17407. https://doi.org/10.1073/pnas.0608396103
Hutagalung AH, Novick PJ (2011) Role of Rab GTPases in membrane traffic and cell physiology. Physiol Rev 91:119–149. https://doi.org/10.1152/physrev.00059.2009
Karpievitch YV, Dabney AR, Smith RD (2012) Normalization and missing value imputation for label-free LC-MS analysis. BMC Bioinformatics 13 Suppl 16: S5. Doi https://doi.org/10.1186/1471-2105-13-S16-S5
Karran E, De Strooper B (2016) The amyloid cascade hypothesis: are we poised for success or failure? J Neurochem 139(Suppl 2):237–252. https://doi.org/10.1111/jnc.13632
Keck M, Androsova G, Gualtieri F, Walker A, von Ruden EL, Russmann V, Deeg CA, Hauck SM, Krause R, Potschka H (2017) A systems level analysis of epileptogenesis-associated proteome alterations. Neurobiol Dis 105:164–178. https://doi.org/10.1016/j.nbd.2017.05.017
Keller L, Xu W, Wang HX, Winblad B, Fratiglioni L, Graff C (2011) The obesity related gene, FTO, interacts with APOE, and is associated with Alzheimer's disease risk: a prospective cohort study. J Alzheimers Dis 23:461–469. https://doi.org/10.3233/JAD-2010-101068
Kiddle SJ, Sattlecker M, Proitsi P, Simmons A, Westman E, Bazenet C, Nelson SK, Williams S, Hodges A, Johnston C et al (2014) Candidate blood proteome markers of Alzheimer's disease onset and progression: a systematic review and replication study. J Alzheimers Dis 38:515–531. https://doi.org/10.3233/JAD-130380
Kilbride SM, Prehn JH (2013) Central roles of apoptotic proteins in mitochondrial function. Oncogene 32:2703–2711. https://doi.org/10.1038/onc.2012.348
Korolainen MA, Nyman TA, Aittokallio T, Pirttila T (2010) An update on clinical proteomics in Alzheimer's research. J Neurochem 112:1386–1414. https://doi.org/10.1111/j.1471-4159.2009.06558.x
Kovacs DM (2000) alpha2-macroglobulin in late-onset Alzheimer's disease. Exp Gerontol 35:473–479
Labbadia J, Morimoto RI (2015) The biology of proteostasis in aging and disease. Annu Rev Biochem 84:435–464. https://doi.org/10.1146/annurev-biochem-060614-033955
Langfelder P, Horvath S (2007) Eigengene networks for studying the relationships between co-expression modules. BMC Syst Biol 1:54. https://doi.org/10.1186/1752-0509-1-54
Langfelder P, Horvath S (2008) WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. https://doi.org/10.1186/1471-2105-9-559
Langfelder P, Mischel PS, Horvath S (2013) When is hub gene selection better than standard meta-analysis? PLoS One 8:e61505. https://doi.org/10.1371/journal.pone.0061505
Langfelder P, Zhang B, Horvath S (2008) Defining clusters from a hierarchical cluster tree: the dynamic tree cut package for R. Bioinformatics 24:719–720. https://doi.org/10.1093/bioinformatics/btm563
Lee SM, Sha D, Mohammed AA, Asress S, Glass JD, Chin LS, Li L (2013) Motor and sensory neuropathy due to myelin infolding and paranodal damage in a transgenic mouse model of Charcot-Marie-tooth disease type 1C. Hum Mol Genet 22:1755–1770. https://doi.org/10.1093/hmg/ddt022
Lepelletier FX, Mann DM, Robinson AC, Pinteaux E, Boutin H (2017) Early changes in extracellular matrix in Alzheimer's disease. Neuropathol Appl Neurobiol 43:167–182. https://doi.org/10.1111/nan.12295
Li NN, Tan EK, Chang XL, Mao XY, Zhang JH, Zhao DM, Liao Q, Yu WJ, Peng R (2013) Genetic association study between STK39 and CCDC62/HIP1R and Parkinson's disease. PLoS One 8:e79211. https://doi.org/10.1371/journal.pone.0079211
Lin WH, Chiu KC, Chang HM, Lee KC, Tai TY, Chuang LM (2001) Molecular scanning of the human sorbin and SH3-domain-containing-1 (SORBS1) gene: positive association of the T228A polymorphism with obesity and type 2 diabetes. Hum Mol Genet 10:1753–1760
Loos RJ, Yeo GS (2014) The bigger picture of FTO: the first GWAS-identified obesity gene. Nat Rev Endocrinol 10:51–61. https://doi.org/10.1038/nrendo.2013.227
McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR, Jr., Kawas CH, Klunk WE, Koroshetz WJ, Manly JJ, Mayeux R et al (2011) The diagnosis of dementia due to Alzheimer's disease: recommendations from the National Institute on Aging-Alzheimer's Association workgroups on diagnostic guidelines for Alzheimer's disease Alzheimers Dement 7: 263–269 Doi https://doi.org/10.1016/j.jalz.2011.03.005
Miller JA, Oldham MC, Geschwind DH (2008) A systems level analysis of transcriptional changes in Alzheimer's disease and normal aging. J Neurosci 28:1410–1420. https://doi.org/10.1523/JNEUROSCI.4098-07.2008
Miners JS, Clarke P, Love S (2017) Clusterin levels are increased in Alzheimer's disease and influence the regional distribution of Abeta. Brain Pathol 27:305–313. https://doi.org/10.1111/bpa.12392
Minjarez B, Calderon-Gonzalez KG, Rustarazo ML, Herrera-Aguirre ME, Labra-Barrios ML, Rincon-Limas DE, Del Pino MM, Mena R, Luna-Arias JP (2016) Identification of proteins that are differentially expressed in brains with Alzheimer's disease using iTRAQ labeling and tandem mass spectrometry. J Proteome 139:103–121. https://doi.org/10.1016/j.jprot.2016.03.022
Mirra SS, Heyman A, McKeel D, Sumi SM, Crain BJ, Brownlee LM, Vogel FS, Hughes JP, van Belle G, Berg L (1991) The consortium to establish a registry for Alzheimer's disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer's disease. Neurology 41:479–486
Montine TJ, Phelps CH, Beach TG, Bigio EH, Cairns NJ, Dickson DW, Duyckaerts C, Frosch MP, Masliah E, Mirra SS et al (2012) National Institute on Aging-Alzheimer's Association guidelines for the neuropathologic assessment of Alzheimer's disease: a practical approach. Acta Neuropathol 123:1–11. https://doi.org/10.1007/s00401-011-0910-3
Moya-Alvarado G, Gershoni-Emek N, Perlson E, Bronfman FC (2016) Neurodegeneration and Alzheimer's disease (AD). What can proteomics tell us about the Alzheimer's brain? Mol Cell Proteomics 15:409–425. https://doi.org/10.1074/mcp.R115.053330
Naudi A, Cabre R, Jove M, Ayala V, Gonzalo H, Portero-Otin M, Ferrer I, Pamplona R (2015) Lipidomics of human brain aging and Alzheimer's disease pathology. Int Rev Neurobiol 122:133–189. https://doi.org/10.1016/bs.irn.2015.05.008
Negyessy L, Nepusz T, Kocsis L, Bazso F (2006) Prediction of the main cortical areas and connections involved in the tactile function of the visual cortex by network analysis. Eur J Neurosci 23:1919–1930. https://doi.org/10.1111/j.1460-9568.2006.04678.x
Oldham MC, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath S, Geschwind DH (2008) Functional organization of the transcriptome in human brain. Nat Neurosci 11:1271–1282. https://doi.org/10.1038/nn.2207
Orsini F, De Blasio D, Zangari R, Zanier ER, De Simoni MG (2014) Versatility of the complement system in neuroinflammation, neurodegeneration and brain homeostasis. Front Cell Neurosci 8:380. https://doi.org/10.3389/fncel.2014.00380
Pore D, Gupta N (2015) The ezrin-radixin-moesin family of proteins in the regulation of B-cell immune response. Crit Rev Immunol 35:15–31
Pugazhenthi S, Qin L, Reddy PH (2017) Common neurodegenerative pathways in obesity, diabetes, and Alzheimer's disease. Biochim Biophys Acta 1863:1037–1045. https://doi.org/10.1016/j.bbadis.2016.04.017
Rajagopalan LE, Westmark CJ, Jarzembowski JA, Malter JS (1998) hnRNP C increases amyloid precursor protein (APP) production by stabilizing APP mRNA. Nucleic Acids Res 26:3418–3423
Ramoz N, Cai G, Reichert JG, Silverman JM, Buxbaum JD (2008) An analysis of candidate autism loci on chromosome 2q24-q33: evidence for association to the STK39 gene. Am J Med Genet B Neuropsychiatr Genet 147B: 1152–1158. Doi https://doi.org/10.1002/ajmg.b.30739
Reitz C, Tosto G, Mayeux R, Luchsinger JA, Group N-LNFS, Alzheimer's Disease Neuroimaging I (2012) Genetic variants in the fat and obesity associated (FTO) gene and risk of Alzheimer's disease. PLoS One 7: e50354. Doi https://doi.org/10.1371/journal.pone.0050354
Sandin M, Chawade A, Levander F (2015) Is label-free LC-MS/MS ready for biomarker discovery? Proteomics Clin Appl 9:289–294. https://doi.org/10.1002/prca.201400202
Scheltens P, Blennow K, Breteler MM, de Strooper B, Frisoni GB, Salloway S, Van der Flier WM (2016) Alzheimer's disease. Lancet 388:505–517. https://doi.org/10.1016/S0140-6736(15)01124-1
Shirasaki DI, Greiner ER, Al-Ramahi I, Gray M, Boontheung P, Geschwind DH, Botas J, Coppola G, Horvath S, Loo JA et al (2012) Network organization of the huntingtin proteomic interactome in mammalian brain. Neuron 75:41–57. https://doi.org/10.1016/j.neuron.2012.05.024
Shulman JM, Imboywa S, Giagtzoglou N, Powers MP, Hu Y, Devenport D, Chipendo P, Chibnik LB, Diamond A, Perrimon N et al (2014) Functional screening in drosophila identifies Alzheimer's disease susceptibility genes and implicates tau-mediated mechanisms. Hum Mol Genet 23:870–877. https://doi.org/10.1093/hmg/ddt478
Small SA, Petsko GA (2015) Retromer in Alzheimer disease, Parkinson disease and other neurological disorders. Nat Rev Neurosci 16:126–132. https://doi.org/10.1038/nrn3896
Soejitno A, Tjan A, Purwata TE (2015) Alzheimer's disease: lessons learned from Amyloidocentric clinical trials. CNS Drugs 29:487–502. https://doi.org/10.1007/s40263-015-0257-8
Stenmark H (2009) Rab GTPases as coordinators of vesicle traffic. Nat Rev Mol Cell Biol 10:513–525. https://doi.org/10.1038/nrm2728
Storey JD (2015) Qvalue: Q-value estimation for false discovery rate control
Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100:9440–9445. https://doi.org/10.1073/pnas.1530509100
Sulistio YA, Heese K (2016) The ubiquitin-proteasome system and molecular chaperone deregulation in Alzheimer's disease. Mol Neurobiol 53:905–931. https://doi.org/10.1007/s12035-014-9063-4
Tannu NS, Hemby SE (2006) Two-dimensional fluorescence difference gel electrophoresis for comparative proteomics profiling. Nat Protoc 1:1732–1742. https://doi.org/10.1038/nprot.2006.256
Udyavar AR, Hoeksema MD, Clark JE, Zou Y, Tang Z, Li Z, Li M, Chen H, Statnikov A, Shyr Y et al (2013) Co-expression network analysis identifies spleen tyrosine kinase (SYK) as a candidate oncogenic driver in a subset of small-cell lung cancer. BMC Syst Biol 7 Suppl 5: S1. Doi https://doi.org/10.1186/1752-0509-7-S5-S1
Vitvitsky VM, Garg SK, Keep RF, Albin RL, Banerjee R (2012) Na+ and K+ ion imbalances in Alzheimer's disease Biochim Biophys Acta 1822: 1671–1681 Doi https://doi.org/10.1016/j.bbadis.2012.07.004
Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH (2011) Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature 474:380–384. https://doi.org/10.1038/nature10110
Wang S, Yang F, Petyuk VA, Shukla AK, Monroe ME, Gritsenko MA, Rodland KD, Smith RD, Qian WJ, Gong CX et al (2017) Quantitative proteomics identifies altered O-GlcNAcylation of structural, synaptic and memory-associated proteins in Alzheimer's disease. J Pathol 243:78–88. https://doi.org/10.1002/path.4929
Wang Y, O'Connell JR, McArdle PF, Wade JB, Dorff SE, Shah SJ, Shi X, Pan L, Rampersaud E, Shen H et al (2009) From the cover: whole-genome association study identifies STK39 as a hypertension susceptibility gene. Proc Natl Acad Sci U S A 106:226–231. https://doi.org/10.1073/pnas.0808358106
Wisniewski JR (2017) Filter-aided sample preparation: the versatile and efficient method for proteomic analysis. Methods Enzymol 585:15–27. https://doi.org/10.1016/bs.mie.2016.09.013
Wisniewski JR, Rakus D (2014) Multi-enzyme digestion FASP and the 'Total protein Approach'-based absolute quantification of the Escherichia coli proteome. J Proteome 109:322–331. https://doi.org/10.1016/j.jprot.2014.07.012
Wisniewski JR, Zougman A, Nagaraj N, Mann M (2009) Universal sample preparation method for proteome analysis. Nat Methods 6:359–362. https://doi.org/10.1038/nmeth.1322
Ye H, Liu W (2015) Transcriptional networks implicated in human nonalcoholic fatty liver disease. Mol Gen Genomics 290:1793–1804. https://doi.org/10.1007/s00438-015-1037-3
Yerbury JJ, Ooi L, Dillin A, Saunders DN, Hatters DM, Beart PM, Cashman NR, Wilson MR, Ecroyd H (2016) Walking the tightrope: proteostasis and neurodegenerative disease. J Neurochem 137:489–505. https://doi.org/10.1111/jnc.13575
Yip AM, Horvath S (2007) Gene network interconnectedness and the generalized topological overlap measure. BMC Bioinformatics 8:22. https://doi.org/10.1186/1471-2105-8-22
Zenaro E, Pietronigro E, Della Bianca V, Piacentino G, Marongiu L, Budui S, Turano E, Rossi B, Angiari S, Dusi S et al (2015) Neutrophils promote Alzheimer's disease-like pathology and cognitive decline via LFA-1 integrin. Nat Med 21:880–886. https://doi.org/10.1038/nm.3913
Zhang B, Horvath S (2005) A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol 4: Article17. Doi https://doi.org/10.2202/1544-6115.1128
Zhang L, Liu YZ, Zeng Y, Zhu W, Zhao YC, Zhang JG, Zhu JQ, He H, Shen H, al TQ (2016) Network-based proteomic analysis for postmenopausal osteoporosis in Caucasian females. Proteomics 16:12–28. https://doi.org/10.1002/pmic.201500005
We are grateful to Emory Center for Neurodegenerative Disease Brain Bank donors and their family. This work was supported by National Institutes of Health (GM103613 and NS092343 to L.L. and NS093550 to L.S.C.); and pilot grant awards from Emory University Research Committee to L.L. and Atlanta Clinical and Translational Science Institute (ACTSI/URC) to L.S.C. Emory Center for Neurodegenerative Disease Brain Bank is supported in part by NIH grants P50 AG025688 and P30 NS055077.
Ethics approval and consent to participate
All research in this study was performed in accordance with the US National Institutions of Health guidelines for research involving human tissues and with the ethical standards and principles of the Declaration of Helsinki. All brain tissues were obtained from Emory Center for Neuro-degenerative Disease Brain Bank. Human postmortem brain tissues were acquired with Institutional Review Board (IRB) approval and informed consent from the subject or the family.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Demographic and neuropathological data of human AD patient and control cases. For each case, the age, gender, disease status, age at onset, disease duration, Braak stage, CERAD neuritic plaque score, frontal cortex neuritic plaque frequency, ApoE genotype, and postmortem interval (PMI) are provided. (XLSX 10 kb)
Table S2. Differential expression analysis of protein abundances in AD and control brains. List of all proteins with complete abundance data in AD and control brains is provided with their fold changes, P values, and q values. Differentially expressed proteins with significantly altered protein abundances (> 1.3-fold change; P < 0.05) in AD versus control are indicated in bold and also provided in separate tabs. (XLSX 660 kb)
Table S3. List of novel proteins with altered abundances in AD identified in the present study. The fold changes of protein abundances in AD versus control are provided with corresponding P values and q values. (XLSX 51 kb)
Table S4. Gene ontology (GO) term enrichment for differentially expressed proteins in AD. The enriched GO terms with associated P values (Benjamini-Hochberg FDR corrected) for biological processes, cellular compartments, and molecular functions are provided in separate tabs. (XLSX 35 kb)
Table S5. Protein co-expression network analysis by WGCNA. Network analysis of the entire proteomic data set from all AD and control cases identified 24 network modules, M1 to M24, coded by different colors according to the convention of WGCNA. Proteins that were not assigned to any module were coded by the color grey in M0. The complete list of proteins in each module and their module membership values (kME) are provided. (XLSX 724 kb)
Table S6. Gene ontology (GO) term enrichment for proteins in WGCNA modules. The enriched GO terms with associated P values (Benjamini-Hochberg FDR corrected) for biological processes, cellular components, and molecular functions are provided in separate tabs. (XLSX 147 kb)
About this article
Cite this article
Zhang, Q., Ma, C., Gearing, M. et al. Integrated proteomics and network analysis identifies protein hubs and network alterations in Alzheimer’s disease. acta neuropathol commun 6, 19 (2018). https://doi.org/10.1186/s40478-018-0524-2
- Alzheimer’s disease
- Brain proteome
- Mass spectrometry
- Quantitative proteomics
- Differential expression analysis
- Protein co-expression network analysis
- Network biology