Person‐specific differences in ubiquitin‐proteasome mediated proteostasis in human neurons

Abstract

BACKGROUND

Impairment of the ubiquitin‐proteasome system (UPS) has been implicated in abnormal protein accumulation in Alzheimer's disease. It remains unclear if genetic variation affects the intrinsic properties of neurons that render some individuals more vulnerable to UPS impairment.

METHODS

Induced pluripotent stem cell (iPSC)‐derived neurons were generated from over 50 genetically variant and highly characterized participants of cohorts of aging. Proteomic profiling, proteasome activity assays, and Western blotting were employed to examine neurons at baseline and in response to UPS perturbation.

RESULTS

Neurons with lower basal UPS activity were more vulnerable to tau accumulation following mild UPS inhibition. Chronic reduction in proteasome activity in human neurons induced compensatory elevation of regulatory proteins involved in proteostasis and several proteasome subunits.

DISCUSSION

These findings reveal that genetic variation influences basal UPS activity in human neurons and differentially sensitizes them to external factors perturbing the UPS, leading to the accumulation of aggregation‐prone proteins such as tau.

Highlights

Polygenic risk score for AD is associated with the ubiquitin‐proteasome system (UPS) in neurons. Basal proteasome activity correlates with aggregation‐prone protein levels in neurons. Genetic variation affects the response to proteasome inhibition in neurons. Neuronal proteasome perturbation induces an elevation in specific proteins involved in proteostasis. Low basal proteasome activity leads to enhanced tau accumulation with UPS challenge.

BACKGROUND

Regulated protein homeostasis is particularly essential for postmitotic, long-lived cells such as neurons. Human brain cells accumulate oxidative injury during aging, which can lead to widespread protein aggregation, increasing the demand for cellular protein degradation systems such as the ubiquitin-proteasome system (UPS) and the autophagy-lysosome pathway (ALP). With more than 800 proteins involved, the UPS degrades most water-soluble proteins in cells. The proteasome comprises a multiprotein catalytic core particle (20S) and multiprotein regulatory components (19S or 11S). The process of UPS starts with an enzymatic cascade that marks proteasome substrates with ubiquitin chains. Once recognized by the proteasome, substrates are deubiquitylated by deubiquitinases, unfolded by ATPases, and translocated into the 20S chamber for hydrolysis.

Impaired UPS activity has been reported in neurodegenerative diseases, including Alzheimer's disease (AD) and Parkinson's disease (PD). Abnormal protein accumulation is a common pathological feature of several neurodegenerative diseases, and reductions in normal protein turnover can contribute to the aggregation and accumulation of amyloid-β (Aβ), tau, and α-synuclein. Therefore, the accumulation of abnormal protein aggregates, including excessively phosphorylated microtubule-associated protein tau (p-tau) and α-synuclein, is likely mediated at least in part by the UPS. Aging is the leading risk factor for neurodegeneration, and studies show reduced UPS activity in aged cells and organs. Thus, disruption of the UPS may be an important link between aging, neurodegenerative disorders, and protein aggregation. However, the molecular mechanisms underlying why some individuals develop neurodegeneration while others of the same age do not are poorly understood. There is also evidence that accumulating toxic aggregates of p-tau, Aβ, and α-synuclein can disrupt proteasome function in transgenic animals and cells. Still, whether reduced UPS activity in neurons causes abnormal protein accumulation is unclear. This study provides evidence that in neurons, basal UPS activity is, in part, encoded within the genetic landscape of an individual and can influence the propensity for protein aggregation.

Human induced pluripotent stem cell (iPSC) -derived neurons (iNs) provide a highly controlled and reproducible model system for the study of disease biology. This system is particularly well-suited to analyze the initiating events in the disease that are directly downstream of genetic risk variants. We have generated iPSC lines from over 50 individuals in the Religious Order Study or Memory and Aging Project (ROSMAP) cohorts. These cohorts are composed of Catholic nuns, priests, and monks across the United States (ROS) and older persons recruited from communities around the Chicago metropolitan area (MAP). Upon enrollment, ROSMAP participants are cognitively unaffected and are primarily over the age of 65. Participants agree to take annual physiological examinations and cognitive function tests and consent to donate their organs, including their brain, after death. Brain tissue is assessed to quantify neuropathology, including Aβ plaques and tau tangles, and several layers of -omic profiles are acquired.

Studying iNs derived from these well-characterized ROSMAP participants that span the clinical and neuropathological spectrum of aging allows for examining the biological processes affected by genetic risk for the disease and elucidating the molecular pathways contributing to risk and resilience for disease. In this study, unbiased proteomic profiling of ROSMAP iNs led us to further interrogate basal UPS activity and tau accumulation in iNs across different genetic backgrounds. These studies revealed person-specific differences in UPS-mediated protein homeostasis in human neurons that were accompanied by differences in the accumulation of aggregation-prone proteins. We also confirm previous reports that accumulation of tau reduces proteasome activity and additionally report that modeling chronic reduced proteasome capacity with mild UPS inhibition induces accumulation of tau, suggesting a positive feedback loop in excitatory neurons that supports tau accumulation following an initiation event. Last, we show that neurons from different genetic backgrounds are differentially vulnerable to mild UPS insult, resulting in differences in tau accumulation. Taken together, we provide evidence for a genetically encoded risk of tau accumulation that impacts UPS activity in human neurons.

Research in context

• Systematic review: A systematic review was performed to examine the state of knowledge regarding genetic background effects on neuronal proteostasis.

• Interpretation: Our findings suggest that genetic variation influences basal ubiquitin-proteasome system (UPS) activity in human neurons and that this genetic variation differentially sensitizes neurons to external factors that can perturb the UPS resulting in the accumulation of aggregation-prone proteins such as tau.

• Future directions: Future directions will involve identifying the genetic variants affecting UPS vulnerability which will allow for the determination of a “UPS vulnerability score.” These studies would provide a foothold for identifying subtypes of neurodegenerative disease that may be especially responsive to interventions that boost proteostasis.

METHODS

iPSC-derived neurons

The iPSC from the ROSMAP cohort and the MAPT^P301L carrier were generated as described previously. Induced neurons were generated and cultured to Day 21 as described previously. An Institutional Review Board approved ROS and MAP studies of Rush University Medical Center. All participants signed informed consent, an Anatomical Gift Act, and a repository consent to allow their data and biospecimens to be repurposed. The iPSC lines were generated following Institutional Review Board (IRB) review and approval through Partners/BWH IRB (#2015P001676), and we certify that the study was performed in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and its later amendments. See Supplemental Table for details on human subjects' cognition and pathology.

Liquid chromatography-tandem mass spectrometry proteomic analysis

For each ROSMAP iN line, 1.2 million Day 4 iNs were plated per well of a six-well plate and were cultured for 21 days. On Day 21, iNs were washed twice with ice-cold Dulbecco's phosphate-buffered saline (DPBS) and flash-frozen in dry ice before being lysed in 100 μL lysis buffer. Liquid chromatography coupled to tandem mass spectrometry (MS) and the analyses were performed as described previously. Briefly, protein concentrations were determined by the bicinchoninic acid method. Tandem mass tag (TMT) labeling and high pH fractionation were performed. MS was conducted with a high-field asymmetric waveform ion mobility spectrometry Pro equipped with Orbitrap Eclipse (Thermo Scientific). Raw data were analyzed using the Proteome Discoverer Suite (version 2.3, Thermo Scientific). MS/MS spectra were searched against the UniProtKB human proteome database. Following spectral assignment, peptides were assembled into proteins. Post-quantification quality control procedures included normalization to a Global Internal Standard, removal of peptides with excessive missing, and imputation with the k-nearest neighbor algorithm. A ComBat algorithm was used to remove variance induced by the cell harvest batch. The final iN proteomics data file included 94,162 peptides with 11,508 unique proteins.

Polygenic risk score calculation

Genotype data were available for 2,067 ROS/MAP subjects between two batches: n_batch1 = 1,686 genotyped using the Affymetrix GeneChip 6.0 and n_batch2 = 381 genotyped using the Illumina OmniQuad Express platform. Details of raw genotype quality control (QC) have been previously described. Briefly, genotypes were first preprocessed using the TOPMed Imputation Server-recommended data preparation pipeline (). Following data preparation, each batch was imputed separately using the TOPMed Imputation Server (TOPMed reference r2), including Eagle (v2.4) for allelic phasing and Minimac4 (v1.5.7) for imputation. The resulting imputed data were then filtered for imputation quality (removing SNPs with r² < 0.8) and minor allele frequency (removing single nucleotide polymorphisms [SNPs] with MAF < 2.5e-3). All tri-allelic variants were removed before the imputed genotype data were merged and all overlapping SNPs mapped to rsIDs (dbSNP build 155). This resulted in a set of 9,329,439 high-quality, bi-allelic autosomal SNPs.

PRSice (v2.3.3) was used to calculate an Alzheimer's disease (AD) polygenic risk score (PRS) using the largest available genome-wide association study (GWAS) on AD and related dementias. To account for patterns of linkage disequilibrium (LD), we applied a 500Kb sliding window and r² threshold of 0.1 during variant clumping. We then set the SNP inclusion p-value threshold to 5e-8, including only genome-wide significant SNPs during variant thresholding. This resulted in 126 independent SNPs that were carried forward to the PRS calculation. An individual's PRS is a weighted sum of their allelic dosage at each of the 126 selected SNPs, weighted by the corresponding effect size estimate of that SNP estimated by the AD GWAS.

Gene set enrichment analyses

Gene set enrichment analyses (GSEA) were performed using the desktop GSEA application (). For GSEA of correlations, we generated rank files using the -log10(p-value)^∗(r-value) to factor both the significance as well as the magnitude and directionality of the effect. KEGG subset of canonical pathways (cp.kegg.v2022.1.Hs.symbol.gmt) was used. One thousand permutations were performed. Normalized Enrichment Scores (NES) reflect the amount of enrichment of each gene set in the top (+) or bottom (-) of the gene rank files that were normalized to account for gene set sizes.

Proteasome activity assay in iN protein lysates

Induced neurons vehicle-treated or treated with 5 nM bortezomib (2-, 8-, 24-, or 72 h) or 50 nM ONX-0914 (72 h) were washed with ice-cold DPBS on Day 21, followed by homogenization in 25 mM Hepes-KOH (pH 7.5), 5 mM MgCl₂, 10% glycerol, 0.1% NP-40, 1 mM dithiothreitol (DTT), 1 mM adenosine triphosphate (ATP), 0.1 mM phenylmethylsulfonyl fluoride (PMSF), 1 mM NaF, followed by sonication in a water bath sonicator (B2500A-MT, VWR) for 20 times for 10 s each, with 30 s pause in between pulses. Protein lysates were centrifuged at 10,000 × g for 10 min at 4°C. Samples normalized for protein concentration (5∼10 μg total protein) were loaded into a black-walled 96-well plate. Proteasome activity was monitored using Infinite 200 PRO plate reader (Tecan) with the i-control microplate reader software (Tecan). The reaction was started by adding 20 μM Suc-LLVY-amc (for chymotrypsin-like activity), or 20 μM Boc-LRR-amc (for trypsin-like activity), or 40 μM Ac-nLPnLD-amc (for caspase-like activity) in 50 mM Tris, 40 mM KCl, 5 mM MgCl₂, 0.05 mg/mL Bovine Serum Albumin (BSA), 1 mM ATP, 1 mM DTT, 1 mM NaF at the zero-time point. For background detection, 2 μM epoxomicin was added to the reactions with Suc-LLVY-amc; 100 μM MG132 was added to reactions with Boc-LRR-amc or Ac-nLPnLD-amc. The activity was assayed by measuring the fluorescence intensity for 1 h at 37°C in 1-min intervals (excitation 380 nm; emission 460 nm).

Western blot analysis

The cell lysates were separated by 4-12% Bis-Tris gels (Thermo Fisher), followed by transfer to nitrocellulose membrane using Criterion blotter system (Bio-Rad). Revert 700 total protein stain (LI-COR) was used to quantify total protein levels. Nitrocellulose membranes were incubated with Intercept blocking buffer (LI-COR) for 1 h at room temperature, followed by overnight incubation with mouse anti-Ubiquitin (Ubi-1) (1:500, Millipore), rabbit anti-PSMB5 (1:2000, Bethyl), rabbit anti-synaptophysin (1:10,000, Abcam), rabbit anti-HOMER1 (1:1,000, Synaptic Systems), chicken anti-MAP2 (1:20,000, Abcam), chicken anti-β-III tubulin (1:10,000, Novus Biologicals), rabbit anti-tau (K9JA, 1:1000, Dako), mouse anti-tau (Tau-5, 1:1000, Invitrogen), mouse anti-pT181-tau (1:1000, Invitrogen), mouse anti-pS202/T205-tau (AT8, 1:250, Invitrogen), rabbit anti-pS214-tau (1:1,000, Cell Signaling), rabbit anti-pT217-tau (1:1000, GeneTex), mouse anti-pS396-tau (1:1000, Cell Signaling), rabbit anti-pS404-tau (1:1000, Invitrogen), mouse anti-glyceraldehyde-3-phosphate dehydrogenase (GAPDH; 1:10,000, Proteintech). Membranes were washed three times with TBST and incubated with fluorescent dye-conjugated secondary antibodies (anti-mouse or anti-rabbit, 1:10,000, LI-COR) for 1 h at room temperature. Membranes were washed three times with TBST and then two times with TBS, followed by scanning using an Odyssey Infrared Imaging System (LI-COR).

Amyloid-β enzyme-linked immunosorbent assay

Quantification of Aβ in the conditioned media collected from the ROSMAP iN cultures was performed using a multiplex enzyme-linked immunosorbent assay (ELISA) for Aβ38, Aβ40, and Aβ42 (V-Plex Plus Aβ Panel, Mesoscale Discovery). Aβ37 was quantified as previously described.

Immunocytochemistry

Day 21 iPSC-derived neurons were washed with ice-cold PBS and fixed with 4% paraformaldehyde in PBS. Following fixation, the cells were permeabilized with 1% Triton X-100 in donkey serum (Jackson ImmunoResearch). The cells were probed with primary antibodies: chicken anti-β-tubulin III (1:500, Novus Biologicals), rabbit anti-Brn2 (1:100, Cell Signaling Technology), followed by secondary antibodies: anti-chicken Cy3 (1:2000, Jackson ImmunoResearch), anti-rabbit Cy5 (1:2000, Jackson ImmunoResearch). DAPI (1 μg/mL, Thermo Analysis) was used for nuclear staining. Images were acquired using an LSM710 confocal microscope (Zeiss).

Statistical analysis

Pearson correlation coefficients were used for examining the relationships between ROSMAP iN proteomics and (i) LOAD PRS, (ii) iN basal proteasome activity, and (iii) p-tau levels. Z-Scores were used to compare basal proteasome activity, PSMB5 (measured by WB), and pan-ubiquitinated protein across ROSMAP iNs. Unpaired t-tests were used to make the comparison between (i) NCI and AD, (ii) MAPT^+/+ and MAPT^P301L/+, (iii) DMSO and Bortezomib-treated, and (iv) DMSO and ONX-0914-treated.

RESULTS

Analyses of neurons from over 50 individuals implicate genetically encoded proteasome vulnerability in AD

To determine the pathways and molecular processes downstream of genetic risk for late-onset AD (LOAD) in human neurons, we performed unbiased proteomic profiling by TMT-MS of iNs from 53 ROSMAP participants (Figure , Supplemental Tables ). IPSC models capture genetic risk and resilience factors from the person from whom they were derived, and in this cohort a range of PRS for LOAD is represented (Figure , Supplemental Table ). We performed correlation analyses between proteomic profiles of iNs and LOAD PRS, followed by GSEA. Strikingly, the only KEGG (Kyoto Encyclopedia of Genes and Genomes) gene set negatively associated with LOAD PRS was “proteasome” (FDR = 0.017) (Figure , Supplemental Table ). Differential protein analysis comparing iNs from non-cognitively impaired (NCI) individuals (n = 35) and those from AD individuals (n = 18) revealed 760 differentially expressed proteins (DEPs) defined by an adjusted p-value < 0.05 and a logFC of > 0.25 or < −0.25 (Supplemental Table ). Pathway analysis of DEPs revealed associations with several processes, including the downregulation of proteins involved in proteostasis (Figure ). Examples of proteins driving these associations include proteasome components, ubiquitin ligases, and deubiquitinases (Figure ). In line with these findings, pseudobulk RNAseq data from human brain comparing AD and NCI reveals that proteasome components are downregulated in excitatory neurons in the AD brain (Figure ).

[IMAGE OMITTED. SEE PDF]

Genetic variation influences basal proteasome activity in human neurons

A key advantage of the iPSC system is that it captures human genetic variation in a well-controlled reductionist system, allowing for the study of biological domains impacted by known and unknown genetic risk and resilience variants. We, therefore, sought to determine whether neurons from these varied genetic backgrounds display differences in proteasome activity as suggested by our unbiased-omics analysis (Figure ). To determine if changes in protein levels of UPS components result in differential proteasome activity, we performed in vitro proteasome activity assays following cell lysis using the fluorogenic substrate Suc-LLVY-amc, which fluoresces when cleaved and detects chymotrypsin-like activity in protein lysates. To assess the proportion of cleavage due to the proteasome, parallel protein lysates receiving no treatment were compared to those receiving epoxomicin, a proteasome inhibitor targeting PSMB5 (the subunit responsible for chymotrypsin-like activity in the 20S proteasome). Of note, chymotrypsin-like activity is the rate-limiting catalytic site in proteasomal proteolysis. Across ROSMAP iNs, a spectrum of basal proteasome activity was observed (Figure , Supplemental Table ). Chymotrypsin-like activity significantly correlated to PSMB5 protein levels, supporting the validity of the LLVY cleavage assay (Figure ). The NGN2 direct induction protocol generates cultures with > 95% of cells expressing neuronal markers. Immunocytochemical staining of neurons (Figure ), proteomic profiles, and Western blotting of neuronal protein lysates (Supplemental Figure ) provide evidence that the differences in proteasome activity observed are not likely to be due to differences in neuronal differentiation efficiency. Rather, analyses of the correlation between TMT-MS profiles and proteasome activity revealed that the top three associated pathways were “oxidative phosphorylation,” “Parkinson's disease,” and “Alzheimer's disease,” each of which was elevated with lower proteasome activity (Figure , Supplemental Tables ). Each of the leading edge proteins driving the negative associations in Figure  is shown in the Venn diagram in Figure . This Venn diagram reveals the highly interconnected nature of these KEGG terms with regard to the specific proteins elevated with low proteasome activity (Figure ).

[IMAGE OMITTED. SEE PDF]

Basal proteasome activity correlates with the level of aggregation-prone proteins in iNs

The accumulation of aggregation-prone proteins is a common feature of several neurodegenerative disorders, including AD, PD, frontotemporal dementia (FTD), other tauopathies, amyotrophic lateral sclerosis (ALS), and Huntington's disease (HD). Here, we examined the relationship between basal proteasome activity and expression levels of aggregation-prone proteins relevant to these diseases, including amyloid-beta (Aβ), microtubule-associated protein tau (MAPT), α-synuclein (SNCA), superoxide dismutase [Cu-Zn] (SOD1), huntingtin (HTT), prion protein (PRNP), RNA-binding protein FUS (FUS), TAR DNA-binding protein 43 (TARDBP), and progranulin (GRN) in iNs. Proteasome activity was significantly associated with protein levels of several of these genes across ROSMAP iNs (Figure ). As proteasome activity increases, APP and Aβ, MAPT, SNCA, PRNP, and SOD1 are reduced (Figures , Supplemental Figure ). Intriguingly, the opposite is true for the RNA-binding proteins TDP-43 and FUS, which are elevated with higher proteasome activity (Figure , Supplemental Figure ). Examination of a correlation matrix of these proteins reveals strong associations among these different aggregation-prone proteins (Supplemental Figure ). While separate from the focus of this study, future experiments are warranted to determine whether there is a causal relationship between levels of these proteins in human neurons (for example, is there a causal link between levels of tau and TDP-43 in this set of iNs similar to what is observed in experimental models of TARDBP mutations).

[IMAGE OMITTED. SEE PDF]

Tau levels showed the strongest association with proteasome activity compared to each of the other aggregation-prone proteins analyzed by TMT-MS (Figure ). Posttranslational regulation of tau is complex, and different modifications of tau affect its function and aggregation state. Thus, we next performed a deeper analysis of tau to look at both the major forms of 50-60 kDa (MAJ) tau, as well as the high molecular weight (HMW) form of aggregated tau previously reported by us and others. In these Western blot analyses, we also examined different phospho-epitopes to determine whether distinct proteoforms of tau are differentially associated with particular proteasome components (Figures , see Supplemental Figure for representative Western blots for each epitope, Supplemental Table for all quantifications). Consistent with TMT-MS data, levels of the major form of total tau (measured with antibody K9JA) are reduced with increasing proteasome activity (Figure ). However, distinct associations are observed when specific phospho-epitopes are quantified, and differences are also observed when examining MAJ versus HMW forms of tau (Figure ). Intriguingly, levels of different phospho-epitopes of tau MAJ are differentially associated with proteasome subunit levels (Figure ). For example, levels of pS214 tau are strongly and inversely associated with several subunits of the 19S ATPases and chaperones (Figure ) but positively associated with 11S immunoproteasome subunits and proteasome assembly proteins (Figure ). Other phospho-tau epitopes, such as pT181 and pS202/T205, show positive correlations with various core proteasome components and negative associations with 11S and proteasome assembly proteins (Figure ).

Phosphorylated HMW-tau forms that associate with cognition also significantly and negatively correlate with the expression of multiple proteasome subunits, including the core catalytic subunits in the 20S proteasome (PSMB5 and PSMB6), the 19S ATPases (PSMC2 and PSMC4), and the deubiquitinases (PSMD7 and PSMD14) in the 19S proteasome (Figure ). The expression of proteasome assembly chaperone proteins (PSMG1, PSMG2, PSMG3, PSMG4, and POMP) positively correlates to the level of several phosphor-epitopes on HMW-tau (Figure ). Additionally, expressions of subunits in the 11S immunoproteasome (PSME1, PSME2, and PSME3) strongly positively correlate to phosphorylated HMW-tau (Figure ).

Accumulation of phosphorylated tau can induce a reduction in proteasome activity

Prior studies have demonstrated the causal role of p-tau accumulation in proteasome impairment in tau transgenic animals (Figure ). Studies of iPSC-derived neurons harboring FTD-associated MAPT mutations previously showed an accumulation of p-tau and dysregulation of the proteostasis network, including proteasome activity and lysosomal function. Here, we validated these findings using our iN system. In human iNs with mutant MAPT (MAPT^P301L/+), we observed increased p-tau (Figure ). In addition, we found that MAPT mutant neurons had reduced proteasome activity and a reduction in PSMB5 levels (Figure ). Interestingly, TMT-MS analyses revealed reduced levels of additional proteasome subunits (Figure , Supplemental Tables , S10), supporting prior findings that MAPT mutations lead to tau accumulation and proteasome disruption. Interestingly, gene sets involved in neuronal activity (KEGG long-term potentiation, KEGG long-term depression) are positively associated with MAPT mutation (Supplemental Table ), supporting the prior finding that proteasome inhibition prevents the turnover of synaptic proteins that mediate synaptic activity.

[IMAGE OMITTED. SEE PDF]

Establishing a neuronal system to model chronic, mild reductions in proteasome activity

Having observed that MAPT mutations can induce accumulation of tau which in turn leads to proteasome dysfunction, we next interrogated whether reductions in proteasome activity in human neurons would induce accumulation of tau. To this end, we aimed to identify conditions that could model a chronic, partial reduction in proteasome activity. We chose two proteasome inhibitors to examine: bortezomib, an Food and Drug Administration (FDA)-approved cancer treatment that binds reversibly to PSMB5, and ONX-0914, a selective inhibitor of the chymotrypsin-like subunit of the immunoproteasome (PSMB8) (Figure ). To explore all protein-level changes in human neurons with mild proteasome inhibition, we first examined the proteomic profile of iNs treated with 5 nM bortezomib for 72 h (Supplemental Table ). 5 nM was chosen based on findings that bortezomib inhibits proteasome-mediated intracellular proteolysis of long-lived proteins with a concentration of ∼0.1 μM that inhibits 50% of the proteolysis and has been shown to kill cells at 24 and 48 h with IC50s of 100 and 20 nM, respectively. At 72 h, neuronal morphology remained grossly intact (Supplemental Figure ), and the level of apoptosis regulators Bax and p53 were unaltered (Supplemental Figure ), while proteasome activity was reduced by 75% with an accumulation of ubiquitinated proteins (Figure , Supplemental Figure ). Using a cutoff of qval < 0.01 and absolute value of logFC > 0.2, 1,361 proteins were differentially expressed with bortezomib treatment for 72 h (Supplemental Table ). Multiple proteasome subunits, including PSMC4, PSMC5, and PSMD1, are significantly increased with bortezomib, suggesting a compensatory response to reduced proteasome activity in neurons. Interestingly, however, the expression of nuclear factor erythroid-derived 2-related factor 1 (NRF1), a well-established master transcription factor regulating the expression of proteasome subunits, was decreased by 10% in bortezomib-treated iNs compared to control after 72 h of low-dose proteasome inhibition (Supplemental Table ). Several deubiquitinases (DUBs), including PSMD7 and PSMD14 that form the minimal DUB-competent complex in the proteasome, are elevated in bortezomib-treated iNs (Figure ). Additionally, multiple ubiquitin-conjugating enzymes, including UBE2J1, whose expression in iNs inversely correlates to LOAD PRS (Supplemental Table ), are upregulated in bortezomib-treated iNs (Supplemental Table ). The increase of proteins that regulate ubiquitination indicates a neuronal response to the global accumulation of ubiquitinated proteins during proteasome inhibition. Intriguingly, HSP70 co-chaperone BCL-2-associated athanogene 3 (BAG3), a necessary component of chaperone-assisted selective macroautophagy, is among the top 5 proteins significantly upregulated in bortezomib-treated neurons compared to control (Figures ), supporting prior findings in other cellular models.

[IMAGE OMITTED. SEE PDF]

PSME1, PSME2, and PSME3 in the 11S immunoproteasome strongly correlate to the levels of p-tau in iNs (Figure ). To understand the impact of immunoproteasome inhibition on proteomic profile in human neurons, we treated iNs with 50 nM ONX-0914, an immunoproteasome-specific inhibitor, for 72 h. ONX-0914 reduced LLVY-amc hydrolysis in iN protein lysates by 50% (Figure ). Surprisingly, levels of pan-ubiquitinated protein were not significantly increased (Figure ), and only 15 proteins were significantly differentially expressed in ONX-0914-treated iNs (Figure , qval < 0.01 and absolute value of logFC > 0.2). Several of these 15 DEPs overlapped with DEPs from the bortezomib analysis, including POMP, PIR, and UBL5 (Figure ).

Genetic variation affects differential vulnerability to proteasome perturbation

To understand the consequence of the extended reduction of proteasome activity on p-tau homeostasis, we treated a set of iNs from eight genetic backgrounds that span the spectrum of proteasome activity with bortezomib at a low dose (5 nM) (Figure ). Proteasome activity in iNs was significantly reduced by 40% after 2 h of the treatment, regardless of basal proteasome activity, and it continued to decrease to approximately 25% of basal levels by 72 h (Figure ). Ubiquitinated protein began to accumulate after only 2 h of the treatment (Figure , Supplemental Figure ). HMW-pS214-tau began to accumulate between 8 and 24 h of proteasome inhibition (Figure ), while HMW-pS202/T205-tau began to accumulate after 24 h of proteasome inhibition (Figure ). Accumulation of HMW-pS214-tau before HMW-pS202/T205-tau supports the previously proposed hypothesis that tau phosphorylation at S214 is a priming event for further phosphorylation of tau at S202/T205 and other phosphorylation sites commonly found in aggregated tau. Intriguingly, iNs with lower basal proteasome activity (from individuals BR107 and BR33) showed a greater accumulation of HMW-p-tau than those with higher basal proteasome activity (from individuals BR09 and BR98) (Figure ). These results indicate that iNs with lower basal proteasome activity are more vulnerable to HMW-p-tau accumulation mediated by proteasome perturbation.

[IMAGE OMITTED. SEE PDF]

Treatment with immunoproteasome-specific ONX-0914 had only minor or no effects on tau accumulation, even though proteasome activity was reduced by approximately 50% (Supplemental Figure ). Thus, while p-tau levels are highly associated with immunoproteasome components (Figure ), the reduction of immunoproteasome activity does not dramatically affect p-tau levels.

DISCUSSION

This study reveals person-specific differences in UPS-mediated proteostasis in human neurons by integrating proteomic profiles, proteasome activity assay, and p-tau characterization in ROSMAP iNs. First, we identify the ubiquitin-proteasome system as a molecular pathway in AD. In ROSMAP iNs with a spectrum of basal proteasome activity, reduced proteasome components are associated with high LOAD PRS and tau accumulation. Second, we find that basal proteasome activity correlates with the level of aggregation-prone proteins in iNs. Third, we reveal the differential responsiveness of neurons in different genetic backgrounds to the proteasome inhibitor bortezomib. Last, we show that proteasome perturbation significantly increases specific proteins involved in proteostasis in iNs. Overall, our data support the hypothesis that individuals with low basal proteasome activity are more susceptible to proteasome perturbation and prone to protein accumulation (schematized in Figure ). Our datasets also implicate the crucial components involved in the UPS that may mediate differential susceptibility to abnormal protein accumulation in human neurons.

[IMAGE OMITTED. SEE PDF]

The data presented herein can be utilized to identify candidate proteins that may mediate the vulnerability of neurons to abnormal protein accumulation by comparing proteins in iNs that are (i) associated with LOAD PRS, (ii) correlated with basal proteasome activity, and (iii) differentially expressed under mild proteasome inhibition (Figure , Supplemental Tables). Multiple proteins in the UPS, including RNF145, ING3, and NBEAL1, are associated with both LOAD PRS and basal proteasome activity. Numerous proteins in the UPS that correlate with basal proteasome activity in iNs are also differentially expressed in iNs with proteasome perturbation (Supplemental Table ), including multiple proteasome subunits and DUBs. AMFR (an E3 ligase), DNAJB14 (a heat shock protein), SEC24D (involved in vesicle trafficking), SRPRB (a transmembrane GTPase on the ER membrane), and VSTM2L (involved in negative regulation of neuron apoptosis) each are differentially expressed in iNs with proteasome inhibition. Interestingly, their protein expression also correlates with both LOAD PRS and basal proteasome activity, suggesting their critical roles in mediating the susceptibility of neurons to abnormal protein accumulation.

Proteostasis in human neurons in the early stage of the disease

Proteostasis is especially crucial in post-mitotic, long-lived neurons, and the function of UPS and ALP, the two major pathways maintaining proteostasis in cells, plays a vital role in neuronal health. Much younger than neurons in aged brains, ROSMAP iNs studied here provide an experimental system for understanding cellular pathways mediating neuropathogenesis directly downstream of genetic risk factors. Here we show an association of LOAD genetic risk with reduced expression of proteins involved in proteostasis in human iNs (Figure ), indicating the genetic contribution to dysfunctional proteostasis pathways in LOAD. Several top hits involved in proteostasis pathways are linked to neurodegeneration. NBR1 (“neighbor of BRCA1 gene”) is an autophagy receptor that accumulates in Lewy bodies. UBE2J1, a ubiquitin-conjugating enzyme, is differentially expressed in AD. Additionally, GSEA suggests that the UPS is a significantly altered pathway in AD iNs, emphasizing the essential role of efficient proteasome-mediated protein turnover in AD. These results indicate that at least a subset of individuals who developed LOAD carried genetic risk factors associated with lower expression of proteins involved in proteostasis. More studies are warranted to identify the genetic drivers for the downregulation of proteins involved in UPS.

The impact of proteasome inhibition on neuronal proteomic profiles

In iNs treated with low-dose bortezomib for 72 h, we observe a dramatic increase in levels of multiple proteins associated with the cellular stress response, including HSPA6 and BAG3, whose elevated expression was also reported in aging human tissues. HSPA6, a not well-studied member of the heat shock 70-kDa proteins, is robustly induced and colocalized with cytoplasmic protein aggregates in human neurons under thermal stress or treated with MG132, a proteasome inhibitor. HSPA6 is the most highly elevated protein in bortezomib-treated neurons (Figure ) and is not present in the mouse or rat genome, underscoring the advantage of a controlled human model system for studying processes of complex neurobiological regulation. BAG3 regulates multiple cellular pathways and can mediate an increase in autophagic flux. Further, BAG3 has been identified as a master regulator in tau homeostasis. Its upregulation has been shown to reduce tau accumulation in primary neurons, while its downregulation exacerbated this effect. Highly elevated BAG3 levels in bortezomib-treated neurons suggest that neurons are responding to partial proteasome inhibition by enhancing the ALP, supporting cooperative crosstalk between the UPS and the ALP in human neurons. Additionally, in the bortezomib-treated neurons, highly elevated ZFAND5, a proteasome activator that binds ubiquitin conjugates and the proteasome, indicates a compensatory effect mediated by the components within the UPS. Interestingly, in both bortezomib- and ONX-0914-treated neurons, an elevation in PIR (pirin, a nuclear redox sensor) protein levels was observed, suggesting an elevation of oxidative stress induced by proteasome inhibition. The increased level of POMP (proteasome maturation protein), an essential factor facilitating the formation of the 20S proteasome, indicates a neuronal response to offset proteasome inhibition. In summary, this study aids in the identification of crucial elements in the proteostasis network (Figure ) and also enhances our understanding of the cellular response to the perturbation of proteostasis in human neurons.

ACKNOWLEDGMENTS

We thank the NeuroTechonology Studio at Brigham and Women's Hospital for providing Zeiss LSM710 Confocal instrument access and consultation on data acquisition and data analysis; iPSC NeuroHub at Brigham and Women's Hospital for technical assistance with iN differentiations; Alfred Goldberg and his lab members at Harvard Medical School for sharing techniques and insights on characterizations of the UPS; members of the Young-Pearse lab and Dennis Selkoe for their constructive input. The results published here are in whole or in part based on data obtained from the AD Knowledge Portal . This work was supported by an Alzheimer's Association Research Fellowship AARF-21-720471 and NIH grants R01AG055909, U01AG072572, U01AG061356, RF1NS117446.

CONFLICT OF INTEREST STATEMENT

The authors declare no conflicts of interest.

CONSENT STATEMENT

Informed consent was obtained by all human participants with proper IRB approvals.

DIVERSITY, EQUITY, AND INCLUSION STATEMENT

This study analyzed cells derived from deceased human subjects in the ROS and MAP cohorts. In these cohorts, race is self-declared and ethnicity is determined using genetic data. These data are provided in supplemental tables. Ethnicity and race distribution in this study match the overall statistics for the cohort. Ongoing studies are expanding the iPSC cohort to purposefully expand diversity in the iPSC lines available for study.

References

Reichmann D, Voth W, Jakob U. Maintaining a healthy proteome during oxidative stress. Mol Cell. 2018;69(2):203‐213. doi: [DOI: https://dx.doi.org/10.1016/j.molcel.2017.12.021]