ABSTRACT
Alzheimer’s disease (AD) is a public health crisis that grows as populations age. Hallmarks of this neurodegenerative disease include aggregation of beta-amyloid peptides and hyperphosphorylated tau proteins in the brain. Variants of the APOE gene are the greatest known risk factors for sporadic AD. As emerging players in AD pathophysiology, extracellular vesicles (EVs) contain proteins, lipids, and RNAs and are involved in disposal of cellular toxins and intercellular communication. AD-related changes in the molecular composition of EVs may contribute to pathophysiology and lend insights into disease mechanisms. We recently adapted a method for separation of brain-derived EVs (bdEVs) from post-mortem tissues. Using this method, we isolated bdEVs from AD patients with different APOE genotypes (n=23) and controls (n=7). bdEVs were counted, sized, and subjected to parallel small RNA sequencing and proteomic analysis. Numerous bdEV-associated RNAs and proteins correlated with AD pathology and APOE genotype. Some of the identified entities have been implicated previously in important AD-related pathways, including amyloid processing, neurodegeneration, and metabolic functions. These findings provide further evidence that bdEVs and their molecular cargo modulate development and progression of AD.
INTRODUCTION
Alzheimer’s disease (AD, (1)) is a public health crisis (2, 3) that continues to grow as the population ages, demanding new insights into pathophysiology. Mechanisms of pathogenesis are manifold and incompletely understood. Synaptic connections are lost through neurodegeneration. Extracellular deposition of β-amyloid (Aβ—fragments of the amyloid precursor protein, APP) occurs in the form of neuritic plaques (NPs), while intracellular aggregation of hyperphosphorylated tau produces neurofibrillary tangles (NFTs) (4–8). Similar to misfolded prions, pathogenic Aβ and tau can template misfolding and protein modifications at recipient sites (9–11). Wide distribution and high density of NPs and NFTs correspond with greater cognitive impairment. Other hallmarks of AD include vascular amyloidosis, gliosis, and blood-brain barrier disruption (12).
Numerous genetic risk factors for AD have been identified. Mutations in genes such as presenilin 1 (PSEN1) and APP contribute to early-onset, familial AD (FAD). In sporadic AD, variants of the apolipoprotein E (APOE) gene (APOE) are the greatest known genetic risk factor (13–16). A central player in lipid homeostasis, APOE assists in intercellular lipid transfer, binding both cholesterol and low-density lipoprotein receptors in addition to binding Aβ. Three major APOE isoforms exist, encoded by alleles ε2, ε3, and ε4, and produce proteins that differ by only 1-2 amino acids but have different binding partner interactions. Whereas APOE2 is protective against AD relative to APOE3, APOE4 is associated with increased risk, particularly in women and certain ethnic populations (17). In Caucasians, the risk of late-onset AD is higher for APOE4 homozygosity (12-fold) and heterozygosity (2-3-fold), with some indication of greater risk for early-onset disease (18). A majority of AD patients in the US have at least one ε4 allele. The critical role of APOE in central nervous system (CNS) lipid metabolism is not limited to AD. The ε4 allele has been associated with increased cholesterol, ceramide, and sphingomyelin in brains of HIV dementia patients (19), a difference not due to inflammation (activated microglia) or astrogliosis. Furthermore, interactions of biological sex (AD disproportionately affects women (20)) and APOE4 have been noted in AD (17,21,22) and AD models (23, 24).
Extracellular vesicles (EVs) comprise a diversity of lipid bilayer membrane-delimited particles that dispose cellular toxins and mediate intercellular communication (25–27). Roles for EVs have recently been recognized in AD pathophysiology (28–30). EVs participate in neurodegeneration in AD, Parkinson’s, and prion diseases in part by spreading misfolded proteins (31–35). In AD, Aβ and tau have been found in or on EVs in AD models and patients (32–34,36–42). EVs are also found “at the scene of the crime”: in amyloid plaques (36, 43). EV cargo betray AD and mild cognitive impairment (MCI) (37,39,44–48). Furthermore, some EVs contribute to amyloid clearance by glia (49–51), suggesting that EVs from healthy cells block AD pathology. However, context matters: in certain models, reducing EV release may diminish pathology (52, 53). Differential expression of various microRNAs was observed in brain tissue (54–56) and serum of AD patients (57), suggesting that the RNA cargo of EVs may contribute to AD pathogenesis.
The study of EVs and APOE, like the study of EVs in AD in general, is just beginning (58–61), but the rapidly growing literature shows that EVs and APOE contribute to AD via multiple mechanisms. At the same time, APOE variants act in lipid homeostasis and trafficking, processes that also govern EV release and uptake. APOE+ EVs have recently been characterized in cerebrospinal fluid (CSF) and cell culture, with changes in fluids of AD and MCI patients (62). In the pigment cell amyloid model system (63), APOE was associated with EVs in the multivesicular body, where it nucleated formation of cytotoxic amyloid aggregates (59, 64) before release from the cell. APOE is prominently sorted into EVs, particularly from neurons and astrocytes (58). APOE shuttling may increase upon exposure to Aβ peptides (58). Consistent with an influence on endosome-origin EV biogenesis, APOE is a well-known determinant of transmission of certain enveloped viruses, a role exerted in the cell of origin, not the recipient cell (65). Studies in the periphery have touched on the role of EVs in the context of APOE and atherosclerosis (e.g., (66–68)), but without finding direct connections.
In this study, we separated and characterized EVs from human AD and control brain tissue. Most studies of EVs in AD to date focus on EV functions or EV-associated biomarkers in biofluids or cell culture supernatant. However, EVs are most likely to act locally, in the tissue of origin. We therefore obtained brain tissue of AD cases with different APOE genotypes as well as control brain tissue. We used our modification of a rigorous method for separation of tissue EVs (69, 70) to obtain and characterize bdEVs in accordance with the recommendations of the Minimal Information for Studies of EVs (MISEV) and related initiatives of the International Society for Extracellular Vesicles (ISEV) (26,71–73). RNA and protein profiling revealed multiple molecular entities associated with Alzheimer’s disease and, to a lesser extent, with APOE genotype in AD.
METHODS
Tissue collection, processing, and approvals
Human brain tissues were obtained from Johns Hopkins Alzheimer’s Disease Research Center. All collections were approved by the Johns Hopkins University Institutional Review Board. Following external examination and weighing of the autopsy brain, the right cerebral hemisphere was cut in coronal slabs and frozen on prechilled metal plates and stored at –80°C. A mixture of temporal cortex and inferior parietal cortex tissues from 23 Alzheimer’s disease (AD) patients and 7 non-AD controls were included in the study (Table 1). AD patients were diagnosed according to Braak (74) and CERAD (75, 76) criteria. For APOE genotyping, genomic DNA was extracted from tissue using the DNeasy kit (Qiagen). Genotyping was performed using standard procedures (77).
Human cortex tissues used in this study
Separation of extracellular vesicles from brain tissue
EVs were separated from brain tissues using our published protocol (70). Before extraction, a small (∼50 mg) piece of tissue was stored at −80°C for later protein and RNA extraction. The remaining frozen tissue was weighed and briefly sliced on dry ice and then incubated in 75 U/ml collagenase type 3 (Worthington #CLS-3, S8P18814) in Hibernate-E solution for 20 min at 37°C. PhosSTOP and Complete Protease Inhibitor (SigmaAldrich PS/PI 4906837001/11697498001) solution were added to stop digestion. The dissociated tissue was spun at 300 × g for 10 min at 4°C. Supernatant was transferred to a fresh tube and spun at 2000 × g for 15 min at 4°C. Supernatant was filtered through a 0.22 μm filter (Millipore Sigma, SLGS033SS) to deplete tissue debris and spun at 10,000 × g for 30 min at 4°C (Thermo Fisher swinging-bucket rotor model AH-650, k-factor 53, acceleration and deceleration settings of 9, using Ultra-Clear tubes with 5 ml capacity). The pellet was resuspended in 100 μl PBS. This resuspension is termed the “10,000 × g pellet” (10K). The 10,000 × g supernatant was concentrated by 100 kilodalton (kDa) MWCO protein concentrator (Thermo Fisher 88524) from 5 ml to 0.5 ml and loaded onto a size-exclusion chromatography column (qEV original, IZON Science SP1-USD, Christchurch, New Zealand) and eluted by PBS. 0.5 ml fractions were collected. The first 3 ml (Fractions 1-6) of the eluate were discarded as the void volume and subsequent 0.5 mL fractions were collected. For the purposes of this study, a total of 2 ml eluate (Fractions 7-10) were pooled and ultracentrifuged for 70 min at 110,000 × g (average) at 4°C (swinging-bucket rotor model TH-641, Thermo Fisher, k factor 114 at max speed, acceleration and deceleration settings of 9, using thinwall polypropylene tubes with 13.2 ml capacity). Supernatant was removed, and the pellet was resuspended in 100 μl PBS as the purified EV fraction. Fractions were stored at -80°C.
Brain homogenate preparation for protein and RNA
For protein extraction, brain homogenates (BH) were prepared by grinding tissue in cold PBS containing PI/PS with a handheld homogenizer (Kontes Pellet Pestle Motor) for 10 sec. RIPA lysis buffer (Cell Signaling Technology 9806) was added and the mixture was sonicated using an ultrasonic ice bath at 20 kHz for 4 × 20 sec, with a 10 sec interval between each sonication. Homogenate was rotated at 4°C for 2 h and spun 15 min at 14,000 × g at 4°C. Supernatant was transferred to tubes and stored at -80°C.
For RNA extraction, Trizol (Thermo Fisher 15596018) was added to frozen brain tissue and homogenized with Lysing Matrix D beads (MP Biomedicals 116913100) using a benchtop homogenizer at speed setting of 4.0 m/s for 20 sec (FastPrep-24, MP Biomedicals). After homogenization, the supernatant was collected and subjected to RNA extraction by miRNeasy Mini Kit (Qiagen 217004) according to the manufacturer’s instructions.
Nano-flow analysis
EV concentration and size distribution were estimated using the nanoFCM flow nanoAnalyzer (NanoFCM Co.) per manufacturer’s instructions. The instrument was calibrated for concentration and size using 200 nm polystyrene beads and a silica nanosphere cocktail (provided by NanoFCM as pre-mixed silica beads with diameters of 68, 91, 113, and 151 nm), respectively. EV and 10K preparations were resuspended in 50 μl of PBS, diluted as needed (typically 1:200 dilution for EV and 1:50 dilution for 10K fractions), and loaded and recorded for 1 minute. Using a calibration curve, the flow rate and side scatter intensity were converted into particle numbers.
RNA extraction and quality control
RNA was extracted by miRNeasy Mini Kit (Qiagen 217004) according to the manufacturer’s instructions. The small RNA profiles of samples were analyzed using capillary electrophoresis by RNA 6000 Pico Kit (Agilent Technologies 5067-1513) on a Fragment Analyzer (Advanced Analytical). Total RNA and small RNA from BH were analyzed using capillary electrophoresis by RNA 6000 Nano Kit (Agilent Technologies 5067-1511) and RNA 6000 Pico Kit.
Small RNA sequencing
EV RNA was concentrated to 6 µl using a Savant SpeedVac Vacuum concentrator. Small RNA libraries were constructed from 50 ng of RNA extracted from brain homogenate or 5 µl of RNA from bdEVs using the Ion Total RNA-Seq Kit V2 (Life Technologies 4475936). Barcodes were attached using the Ion Xpress™ RNA-Seq Barcode 1-16 Kit (Life Technologies 4471250) according to the manufacturer’s protocol and as previously published (78). The yield and size distribution of the small RNA libraries (96 nt to 250 nt) were assessed using the Agilent 2100 Bioanalyzer™ instrument with DNA 1000 chip (Agilent Technologies 5067-1504). Multiplexed libraries were equally pooled to 45 pM and prepared for deep sequencing using the Ion Chef system (Life Technologies 4484177) and sequenced on the Ion Torrent S5™ using Ion™ 540 chips (Life Technologies A27765).
RNA sequencing data analysis
Original BAM files were converted into FASTQ format using picard tools (SamToFastq command). Reads shorter than 15 nt were removed from the raw FASTQ data using cutadapt software v1.18. The size-selected reads were aligned to human reference transcriptomes using bowtie software (1 mismatch tolerance) in a sequential manner. Specifically, reads were first mapped to rRNA, tRNA, RN7S, snRNA, snoRNA, scaRNA, VT-RNA, Y-RNA as well as the mitochondrial genome. All reads which did not map to the above RNA species were aligned to human miRNA references (miRBase 22 release). The remaining reads were further aligned to protein-coding mRNAs and long non-coding RNA (lncRNA) references (GENCODE Release 29).
The numbers of reads mapped to each RNA type were extracted using eXpress software based on a previous publication (79). To compare the RNA biotype distribution, reads of different RNA biotypes were normalized as reads per million total reads (RPM). Differential gene expression was quantified using R/Bioconductor packages DESeq or DESeq2 as described (80, 81).
Mass spectrometry
Samples were resuspended in 1 X RIPA buffer (20mM Tris-HCl pH7.5, 150mM NaCl, 1mM Na2EDTA, 1mM EGTA, 1% NP-40, 1% SDS, 2.5mM sodium pyrophosphate, 1mM β-glycerophosphate, 1mM Na3VO4, 1μg/ml leupeptin) and protease inhibitors and incubated on ice for 5 min. The samples were sonicated for 15 min in an ice water bath before centrifugation at 14,000 × g at 4°C for 10 min. The supernatant was collected and assessed for protein concentration using the micro BCA assay (Thermo Fisher Scientific 23235). 3 μg of brain homogenate and 1.5 µg of 10K pellet and EV samples were buffer-exchanged prior to mass spectrometry to remove detergent. Proteins were resuspended in 8M Urea, 50 mM Tris pH=8.3. 1 µL of TCEP (tris [2-carboxyethyl] phosphine hydrochloride, 200 mM solution in water) was then added to the samples and incubated for 4 hours at 21°C in a ThermoMixer (Eppendorf AG). 4 µL of 1M IAA (iodoacetamide in water) was then added and samples were incubated in the dark at 21°C. 800 µL of 50 mM Tris (pH 8.3) and 1 μg trypsin were then added to samples prior to overnight incubation at 37°C. 10 µL of 10% trifluoroacetic acid (TFA) was added to each sample to acidify. Samples were cleaned using stage-tips preparations using 3 plugs of Empore polystyrenedivinylbenzene (SBD-XC) copolymer disks (Sigma Aldrich, MO, USA) for solid phase extraction. Peptides were reconstituted in 0.1% formic acid and 2% acetonitrile and loaded onto a trap column (C18 PepMap 100 μm i.d. × 2 cm trapping column, Thermo Fisher Scientific) at 5 µL/min for 6 min using a Thermo Scientific UltiMate 3000 RSLCnano system and washed for 6 min before switching the precolumn in line with the analytical column (BEH C18, 1.7 μm, 130 Å and 75 μm ID × 25 cm, Waters). Separation of peptides was performed at 45°C, 250 nL/min using a linear ACN gradient of buffer A (water with 0.1% formic acid, 2% ACN) and buffer B (water with 0.1% formic acid, 80% ACN), starting from 2% buffer B to 13% B in 6 min and then to 33% B over 70 min followed by 50% B at 80 min. The gradient was then increased from 50% B to 95% B for 5 min and maintained at 95% B for 1 min. The column was then equilibrated for 4 min in water with 0.1% formic acid, 2% ACN. Data were collected on a Q Exactive HF (Thermo Fisher Scientific) in Data Dependent Acquisition mode using m/z 350–1500 as MS scan range at 60 000 resolution. HCD MS/MS spectra were collected for the 7 most intense ions per MS scan at 60 000 resolution with a normalized collision energy of 28% and an isolation window of 1.4 m/z. Dynamic exclusion parameters were set as follows: exclude isotope on, duration 30 s, and peptide match preferred. Other instrument parameters for the Orbitrap were MS maximum injection time 30 ms with AGC target 3 × 106, MSMS for a maximum injection time of 110 ms with AGT target of 1 × 105.
Proteomics data analysis
Human protein sequences (last modified date: 16 May 2019) were downloaded from the Uniprot database and used as the database for the search engine. Common Repository of Adventitious Proteins (CRAP) was used as the potential lab contaminant database. Protein identification was performed using the proteomics search engine Andromeda built in to Maxquant V 1.16.0. Trypsin with a maximum of two missed cleavages was used as the cleavage enzyme. Carbamidomethyl of cysteine was set as fixed modification and oxidation of methionine was set as variable modification. The Percolator results were set to reflect a maximum of 1% false discovery rate (FDR). The Label Free quantification was done with match between runs using a match window of 0.7 min. Large label free quantification (LFQ) ratios were stabilized to reduce the sensitivity for outliers. For human datasets, data scaling was done using the cyclic loess method, and scaled data were visualized with a PCA plot. For differential abundance analysis, nested factorial design was set up for the analysis, where each subtype of the disease was nested within the main disease category and contrasts for the main categories were computed by averaging the subtypes.
The protein interaction and cluster protein function prediction was done by Protein-Protein Interaction Networks Functional Enrichment Analysis (STRING) (82). Kyoto Encyclopedia of Genes and Genomes (KEGG) (83) was used to enrich pathway involvement of identified proteins. Statistical significance of enrichment was determined by the tools mentioned above. Only nominally significant categories (false discovery rate (FDR) < 0.05) were included for analysis.
Electrochemiluminescence-linked (ECL) immunoassay for total tau and phosphorylated tau detection
Total tau and phosphorylated tau at threonine 231 (phosTau) were measured in BH, 10K, and EVs using an ECL-immunoassay (Meso Scale Discovery K15121D) according to the manufacturer’s instructions. In brief, 10K and EV samples were diluted 1:10 while BH was diluted 1:100 in 2% block buffer (2% BSA in tris buffer) containing 0.5% triton X-100. Samples were incubated for one hour on the plate. After washing the plate, the SULFO-TAG Anti-Total Tau Antibody was added and incubated with the plate for one hour. After washing, MSD Read Buffer was added, and the plate was read immediately with a Quick plus SQ 120 MM instrument. Data analysis was done on MSD DISCOVERY WORKBENCH software version 2.0.
Statistical analysis
Statistical significance of particle count, particle: protein ratio, size distribution, RNA biotype differences, and tau and phosTau level differences between AD and control groups and between either of two APOE genotype groups were determined by two-tailed Welch’s t-test.
Data availability
Nucleic acid sequencing data have been deposited with the Gene Expression Omnibus, accession GSE159541. Proteomics data files are available on request. We have submitted all relevant details of our experiments to the EV-TRACK knowledgebase (EV-TRACK ID: EV200126) (73).
RESULTS
Separation of EVs from AD and control brain tissue
Following the protocol illustrated in Figure 1A, we separated brain-derived extracellular vesicles (bdEVs) from control (n=7) and AD (n=23) individuals with different APOE genotypes (Table 1). A small amount (∼50 mg) of each tissue was set aside to produce brain homogenate (BH) to assess protein and RNA profiles of the source material. After enzymatic digestion of the remaining tissue and initial filtering, 10,000 x g ultracentrifuged pellets were collected and termed “10K” as an intermediate product of EV separation. The 10K supernatant was then separated by size exclusion chromatography (SEC) (70) and concentrated into a more pure EV preparation. Fractions were processed for characterization, including proteomics and small RNA profiling.
Alzheimer’s disease (AD) and control brain tissue-derived EV (bdEV) enrichment and characterization. (A) Workflow for bdEV enrichment, small RNA sequencing and proteomics. Following digestion, centrifugation, and filtration steps, 10,000 x g pellets from AD and control brain tissue (BH) (as indicated in Table 1) were collected and defined as the 10K fraction. Size-exclusion chromatography (SEC) was applied to 10,000 x g supernatants to enrich EVs. RNA and proteins from BH, 10K, and EVs were then isolated and subjected to small RNA sequencing and mass spectrometry. (B) Left: Particle concentrations of 10K and EV fractions of AD and control samples were measured by NFCM. Particle concentration for each group was normalized by tissue mass (per 100 mg). Right: Ratio of particles to protein (particles/µg). Protein concentrations of 10K and EV fractions were measured by BCA protein assay. Data are mean +/- SD. ns: no significant difference (p > 0.05) between AD and control by two-tailed Welch’s t-test. (C) Size distributions of 10K and EV fractions (AD and control) were measured by NFCM and calculated as particles in a specific size bin versus total detected particles in each sample (percentage). Data are presented as mean +/- SD. *p ≤ 0.05, ****p ≤ 0.0001 by two-tailed Welch’s t-test. (D) Principal component analysis (PCA) based on extracellular vesicle (EV) marker expression per proteomics analysis. EV markers used for PCA: CD81, CD9, FLOT1, FLOT2, RAB1A, RAB7A, TUBA1B, TUBB4B, ANXA2, ANXA5, ANXA6, ACTN1, GAPDH. Data are presented as mean +/- SD.
Basic EV characterization
Particle concentration per 100 mg tissue input and particle size distribution were determined by nano-flow cytometry measurement (NFCM). No significant particle yield difference was detected between AD and control brain-derived 10K and EV fractions (Figure 1B left).
Particle:protein ratio was also calculated to evaluate EV purity. This ratio was similar between AD and control group (Figure 1B right). However, the EV fraction had a higher particle:protein ratio compared with 10K, as expected and consistent with greater protein contamination of the 10K pellet (Figure 1B right). Fractions from AD and control groups also had similar size distributions. However, more small particles in the approximately 45-50 nm diameter range were observed in AD samples compared with controls (Figure 1C). No significant particle count differences or size distribution shifts were observed for samples with different APOE genotypes (ε2/3, ε3/3, ε3/4/, and ε4/4; Figure S1A-B). We examined the expression levels of 13 proteins that are commonly reported to be associated with EVs (CD81, CD9, FLOT1, FLOT2, RAB1A, RAB7A, TUBA1B, TUBB4B, ANXA2, ANXA5, ANXA6, ACTN1, GAPDH). These proteins were not found to be significantly differentially abundant between AD and controls. Principal component analysis (PCA) also showed different EV marker expression patterns between BH, 10K, and EV fractions, but not between AD and controls (Figure 1D).
Particle yield and size distribution of bdEVs from AD samples with different APOE genotypes. Particle concentration (A) and size distribution (B) of 10K and EV fractions from AD patients classified by APOE genotype as measured by NFCM. Data are presented as mean +/- SD. ns: no significant difference (p > 0.05) between the two APOE genotype groups by two-tailed Welch’s t-test.
Assessment of RNA biotype distribution
Small RNA sequencing of BH, 10K, and EV fractions yielded 3.2M (± 1M), 2.2M (± 0.9M), and 0.75M (± 0.36M) reads, respectively (M = million, 1 x 10^6). After adapter clipping and removing reads shorter than 15 nt, 71.28% (± 4.4%) of BH, 73.79% (±6.8%) of 10K, and 41.83% (± 9.0%) of EV reads mapped to the human genome (hg38). Reads mapped to various RNA biotypes were normalized to reads per million mapped reads (RPM) (Fig S2A-C). Overall, no major differences in biotype distribution were observed between AD and control conditions, although AD EVs showed a slight decrease in mitochondrial RNA (mtRNA) (Figure S2C). We then determined enrichment of ncRNA biotypes in EVs relative to tissues and assessed whether the enrichment pattern differed between AD and control brains. To do this, we normalized the RPM of each ncRNA biotype to its average RPM in brain tissues (Figure S2D). Consistent with earlier reports (70,84–86), miRNA and snoRNA were relatively underrepresented in EVs compared with brain homogenates. The same was observed for mtRNA and, to a lesser extent, for snRNA. Y-RNA was equally abundant in BH, 10K, and EVs. tRNA was enriched in 10K and was equally abundant in EVs and BH. RN7S (SRP-RNA) was significantly enriched in the purified EV fraction of AD samples, but not of controls. Vault RNAs were strongly enriched in EVs, although this may have been due to co-purification of extracellular Vault particles with EVs, as reported earlier (87). Interestingly, both protein-coding (mRNAs) and long non-coding RNAs (lncRNAs) were enriched in EVs compared with BH and 10K, and these reads mapped to thousands of different genes. However, because of our size selection and library preparation methods, these reads likely represent the RNA degradome and not full-length RNAs.
Small RNA biotypes in BH, 10K, and EV fractions from AD and control patients. (A) The log2 reads per million mapped reads (Log2RPM) of small RNA biotypes in BH (top), 10K (middle), and EVs (bottom) of AD and control. Data are presented as mean +/- SD. **p ≤ 0.01 by two-tailed Welch’s t-test. (B) Normalized RNA class read counts in 10K and EVs were scaled to the normalized RNA class read counts in BH. Data are presented as mean +/-SD. *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001 by two-tailed Welch’s t-test.
EVs are associated with differentially expressed miRNAs related to AD pathology
Examining individual miRNAs, 26 miRNAs (10K) and 16 miRNAs (EVs) differed significantly (FDR < 0.05) between AD and controls (Figure 2A). Of these, two miRNAs in 10K and nine miRNAs in EVs differed by more than two-fold (Figure S3A-B). Table 2 provides a summary of these dysregulated miRNAs, indicating up/downregulation and involvement in AD as reported in previous studies, since up/downregulation of some EV-associated miRNAs was consistent with reports of brain tissue expression (54–56). By APOE genotype, focusing on the presumed highest-risk ε4/4 and the lowest-risk ε2/3 carriers in our study, nine miRNAs in 10K differed more than two-fold (with p value < 0.05), while only two miRNAs (miR-379-5p and miR-199a-5p) differed in EVs (Figure 2B, S4A-B, and Table 3). However, none of the same miRNAs were found in the comparisons of AD versus controls and APOE ε4/4 versus ε2/3.
Normalized count of miRNAs differentially expressed between AD and control with fold change > 2 in 10K (A) and EVs (B) of patients classified by APOE genotype. *p value < 0.05 by DeSeq2 differential gene expression analysis based on the negative binomial distribution.
bdEV miRNAs with differential expression in AD. (A) Volcano plots showing 10K (left) and EV (right) miRNA log2 fold changes (Log2FC) and p values (pval) for AD versus control. (B) Volcano plots showing 10K (left) and EV (right) miRNA log2 fold changes (Log2FC) and p values (pval) for APOE ε4/4 and APOE ε2/3 carriers. (A) and (B): Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05, and red: FDR < 0.05. (C) miRNA log2 fold changes (Log2FC) for AD and control for BH were plotted against 10K (upper left) and EVs (upper right). Venn diagrams of miRNAs differentially expressed between AD and controls in 10K (bottom left) versus BH, and EVs versus BH (bottom right). (D) The normalized abundance of eight miRNAs dysregulated in AD and control in common in BH, 10K, and EV fractions. Data are presented as mean +/- SD. (E) miRNA log2 fold changes (Log2FC) between APOE ε4/4 and APOE ε2/3 in BH were plotted against 10K (upper left) and EVs (upper right). Venn diagrams of miRNAs differentially expressed between APOE ε4/4 and APOE ε2/3 in 10K (bottom left) versus BH, and EVs versus BH (bottom right). (C) and (E): Dashed lines indicate log2FC of one (up or down). Colored dots indicate miRNAs with p values < 0.05 in BH (blue), in 10K/EV (red), or in both BH and 10K/EVs (purple), or unchanged transcripts (grey). Dot size represents the mean normalized abundance of individual miRNAs. (F) The normalized abundance of two miRNAs dysregulated in APOE ε4/4 vs APOE ε2/3 in BH, 10K, and EV fractions. Data are presented as mean +/- SD.
Significantly different miRNAs (adjust P value < 0.05) with at least twofold change between AD and control in EVs and 10K. Up/downregulation of these EV miRNAs was compared with brain miRNA change in three earlier studies. Studies that have identified involvement of these miRNAs in AD are listed, as well as the potential target genes of the miRNAs.
Significantly different miRNAs (adjusted P value < 0.05) with twofold change between APOE ε4/4 and ε2/3 carrier in EVs and 10K. Up/downregulation of these EV miRNAs was compared with brain miRNA change in three earlier studies. Studies that have identified involvement of these miRNAs in AD are listed, as well as potential target genes of the miRNAs.
To assess whether AD-associated miRNA differences in the extracellular space reflect overall changes in brain tissue, we compared miRNA fold-changes of 10K and EVs with BH. We highlighted miRNAs that were differentially abundant only in BH (blue dots, Figure 2C), only in 10K or EV (red dots), or in both EVs and BH (purple dots). Around 30% of differentially abundant EV miRNAs reflected differences that could also be observed in BH (18 miRNAs in 10K and 13 in EV, Figure 2C). miRNAs that showed robust differential abundance (FDR < 0.05 and log2 fold-change > 2) were found within this subset of miRNAs. Eight miRNAs showed consistent changes across BH, 10K, and EVs (Figure 2D). Analyzing the higher- and lower-risk APOE genotypes as above, more than 90 miRNAs were differentially abundant in BH (blue dots), while relatively few miRNAs differed in EVs (red dots). Of the differentially abundant miRNAs, 14 showed consistent change between 10K and BH, but only three between EVs and BH (Figure 2E). Curiously, the red blood cell-specific miR-451a was consistently differentially abundant across all fractions (Figure 2F). Whether this is a coincidental contaminant or indicative of some biological difference may require further investigation. miR-328-2p, in contrast, was downmodulated in BH and 10K, but upmodulated in the EV fraction (Figure 2F).
Differential expression of other non-coding RNA biotypes in brain tissue and EVs related to AD pathology
Since incorporation into EVs of non-miRNA non-coding RNAs can be modulated by external stimuli imposed on cells (86, 88), and since these ncRNAs may also contribute to intercellular communication (89, 90), we examined individual snoRNAs, snRNAs, Y-RNAs, and tRNAs. snoRNAs were relatively depleted from EVs (Figure S2D), and we did not find any significant differences in snoRNA content of BH, 10K or EVs in AD versus controls (data not shown). Additionally, no differences in snRNA were found for 10K or EVs in AD versus controls, although a significant upregulation of RNU4-2 in AD brain homogenates was observed (Figure S5A). We also observed slight differences in Y-RNA content of AD versus control EVs (Figure S5B). Compared with non-AD controls, 10K fractions of AD patients had higher levels of Y4-RNA (p value < 0.05), while EVs contained lower levels of Y1-RNA (p value < 0.05).
Normalized count of miRNAs differentially expressed in APOE ε4/4 versus APOE ε2/3 carriers with fold change > 2 in 10K (A) and EVs (B). *p value < 0.05 by DeSeq2 differential gene expression analysis based on the negative binomial distribution.
Differential expression of RNU and Y-RNAs related to AD pathology. (A) RNU4 expression in AD versus control in BH. Data are presented as mean +/- SD. (B) Y-RNAs differentially expressed in AD versus control in BH, 10K, and EVs. Data are presented as mean +/- SD. *p ≤ 0.05 as determined by DeSeq2 differential gene expression analysis based on the negative binomial distribution.
Interestingly, we found many differences in tRNA content of EVs in AD versus controls (Figure 3A). Most of these differences were subtle, however, as none of the tRNAs were enriched more than twofold. In EVs, we observed changes in various isodecoders of tRNA-Cys-GCA and tRNA-Gly-GCC. Additionally, we observed changes in tRNA-Val, tRNA-Arg, tRNA-Pro, and different tRNA-Ala isoacceptors. In contrast, no clear differences in tRNA content were found in 10K. Next, we compared the tRNA content of 10K and EVs from APOE ε4/4 versus ε/3 carriers (Figure 3B). The fold-differences in tRNA dependent on APOE genotype were larger than those dependent on AD status, although fewer significant tRNAs were found. Two tRNA-Gly-CCC isodecoders were enriched in 10K from APOE ε2/3 carriers, while no significant differences were observed in EVs (as defined by FDR < 0.05).
bdEV tRNAs with differential expression in AD. (A) Volcano plots showing 10K (left) and EV (right) tRNA log2 fold changes (Log2FC) and p values (pval) for AD versus control. (B) Volcano plots showing 10K (left) and EV (right) tRNA log2 fold changes (Log2FC) and p values (pval) for APOE ε4/4 versus APOE ε2/3. (A) and (B): Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05, and red: FDR < 0.05. (C) tRNA log2 fold changes (Log2FC) for AD and control for BH were plotted against 10K (left) and EVs (right). (D) tRNA log2 fold changes (Log2FC) between APOE ε4/4 and APOE ε2/3 carriers for BH were plotted against 10K (left) and EVs (right). (C) and (D): Dashed lines indicate log2FC of one (up or down). Colored dots indicate tRNAs with p values < 0.05 in BH (blue), in 10K/EV (red), or in both BH and 10K/EVs (purple), or unchanged transcripts (grey). Dot size represents the mean normalized abundance of individual tRNAs.
We subsequently compared the differences in EV-associated tRNA with differences in brain homogenates (Figure 3C). Like the miRNAs, most of the differences were observed only in the brain homogenates, or in EVs, while a minor number of tRNAs differed in both. Comparing APOE genotype-dependent differences in the tRNA content of EVs and BH, we observed a great difference in BH, or in EVs, but only a small number of tRNAs that differed in both (Figure 3D).
Proteomics and regulatory pathways
Analysis of data collected by label-free mass spectrometry identified 1675, 801, and 1197 proteins in BH, 10K, and EVs (Supplemental Table 1), respectively. Most proteins were identified in both AD and control in BH (73.6%), 10K (66.5%), and EVs (66.6%) (Figure 4A, Figure S6A). Enrichment of common and AD “unique” proteins in BH, 10K, and EVs according to presumed regulatory pathways as identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database was observed for pathways directly or indirectly related to neurodegenerative diseases (Figure 4B, Figure S6B). Metabolic pathways (e.g., carbon metabolism, glycolysis, and oxidative phosphorylation) were widely enriched for common and AD-unique proteins. Proteins with known involvement in Alzheimer’s disease, Parkinson’s disease, and Huntington’s disease were enriched in AD-unique proteins in 10K and EVs (Figure 4B) and common proteins in BH (Figure S6B). Principal component analysis (PCA) of the EV proteome content showed a separation of AD and control groups (Figure 4C, right). In contrast, no such separation was observed for 10K (Figure 4C, left) or BH (Figure S6C).
Proteins related to AD pathology in brain tissues. (A) Venn diagram of proteins identified in BH of AD and control. (B) The top 10 pathways ranked by FDR-corrected p value of BH proteins according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). (C) Principal component analysis (PCA) based on proteome content of BH. (D) Volcano plots showing BH protein log2 fold changes (Log2FC) and p values (pval) for AD versus control. Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05, and red: FDR < 0.05.
Number of proteins identified in each type of sample/fraction
bdEV protein content reflects differences between AD and control brain tissues. (A) Venn diagrams of proteins identified in 10K and EVs from AD and control patients. (B) Top 10 pathways ranked by FDR-corrected p value of 10K and EV proteins according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). (C) Principal component analysis (PCA) based on proteome content of 10K (left) and EVs (right). (D) Volcano plots showing 10K (left) and EV (right) protein log2 fold changes (Log2FC) and p values (pval) for AD versus control. Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05, and red: FDR < 0.05. (E) Expression level of proteins differentially expressed between AD and control with fold change > 2 (non-adjusted p value < 0.05) in 10K (left) and EVs (right). Data are presented as mean +/- SD. (F) String protein interaction network analysis indicated that 18 proteins (out of 31 proteins showing differential expression (non-adjusted p value < 0.05) between AD and control in EVs) were enriched with known high confidence (0.7 on a scale of 0-1) protein-protein interactions. Protein clusters are indicated with different colors based on predicted functions.
Differentially abundant proteins: AD versus control
Label-free quantitation (LFQ) was used to identify up- and downregulated proteins in AD (Figure 4D and Figure S6D). More proteins were up- or downregulated in AD in the EV fraction (Figure 4D, right) than in 10K (Figure 4D left) or BH (Figure S6D). EVs have increased potential to indicate the difference between AD and controls as compared with 10K and BH. Three proteins (10K) and 11 proteins (EVs) differed by more than two-fold between AD and control (Figure 4E, S7A-B). Examining proteins that were differentially abundant in EVs by Protein-Protein Interaction Networks Functional Enrichment Analysis (STRING), 18 out of 31 had high protein to protein interaction confidence scores (0.7 on a scale of 0-1), participating in AD-related processes such as neurodegeneration, neurotrophin signaling, oxidative regulation, and metabolic regulation (Figure 4F). Prominent among these was microtubule-associated protein tau (MAPT), which forms neurofibrillary tangles in AD brains. To validate these results, we next measured the concentration of total tau and tau with phosphorylated threonine 231 (phosTau) in lysed BH, 10K, and EVs by ECL-immunoassay. Normalized to total protein input, both total tau and phosTau were significantly increased in bdEVs of AD versus controls, but not in the 10K fraction (Figure 5A). In BH, only phosTau was significantly increased in AD compared with controls (Figure 5A). Results were similar when data were normalized by starting amount of tissue (Figure S8). We then assessed the possible correlations of tau or phosTau in BH and 10K or EV fractions. No strong correlation was observed for total tau levels (BH vs 10K or EVs; Figure 5B) or for phosTau in BH and 10K. However, phosTau levels in BH and EVs were positively correlated (R=0.48, p=0.007; Figure 5B).
Expression level of proteins differentially expressed between AD (n=23) and control (n=7) with fold change > 2 in 10K (A) and EVs (B). Data are presented as mean +/- SD. *p value < 0.05 by Nested Factorial design differential abundance analysis.
Higher total tau and phosphorylated tau protein levels in bdEVs from AD patients vs controls. The level of total tau and phosphorylated tau protein per 100 mg tissue input in 10K and EVs. Data are presented as mean +/- SD. ns: no significant difference (p > 0.05), **p ≤ 0.01 by two-tailed Welch’s t-test.
Higher total tau and phosphorylated tau protein levels in bdEVs from AD patients compared with controls. (A) The level of total tau and phosphorylated tau protein per 100 µg protein input in BH, 10K, and EVs. Data are presented as mean +/- SD. ns: no significant difference (p > 0.05), *p ≤ 0.05, ***p ≤ 0.001, ****p ≤ 0.0001 by two-tailed Welch’s t-test. (B) Scattered plot of tau and phosphorylated tau protein levels between BH with 10K (left), and BH with EVs (right) in AD and controls. Trendlines are shown in black, and the grey area depicts 95% confidence intervals. Pearson correlation coefficient (R) and significance (p) are shown based on n = 30 AD and control samples.
Proteomics: APOE ε4/4 and APOE ε2/3
The analyses above were next restricted to APOE ε4/4 and ε2/3 genotypes (AD samples only). Most proteins were detected in more than 50% of individuals in each group (i.e., n > 3 in each; Figure 6A, Figure S9A). Most proteins were found in common in the ε4/4 and ε2/3 groups (71.4% in 10K, 70.3% in EVs, and 81.5% in BH). KEGG pathway analyses of APOE ε4/4-“unique” proteins in 10K and EVs corresponded with metabolism-related pathways (carbon, glutathione, and galactose metabolism) (Figure 6B). KEGG pathway analysis of APOE ε2/3-“unique” proteins in 10K and EVs did not return enriched pathways (data not shown), although numerous pathways were enriched in BH of APOE ε2/3 carriers (Figure S9B). Moreover, both ε4/4 and ε2/3 contained enriched proteins in BH that were involved in Alzheimer’s disease, Huntington’s disease, and oxidative phosphorylation (Figure S9B).
Protein related to APOE genotypes in brain tissues. (A) Venn diagram of proteins identified in BH from APOE ε4/4 and APOE ε2/3 carriers. (B) Volcano plots showing BH protein log2 fold changes (Log2FC) and p-value (pval) for APOE ε4/4 and APOE ε2/3 carriers. Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05, and red: FDR < 0.05. (C) The pathways significant for APOE ε4/4 and APOE ε2/3 unique proteins in BH according to the Kyoto Encyclopedia of Genes and Genomes (KEGG).
bdEV proteins with differential expression in APOE ε4/4 and APOE ε2/3 carriers. (A) Venn diagrams of proteins identified in 10K and EVs from APOE ε4/4 and APOE ε2/3 carriers. (B) Significant pathways of APOE ε4/4 10K and EV proteins according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). (C) Volcano plots showing 10K (left) and EV (right) protein log2 fold changes (Log2FC) and non-adjusted p values (pval) between APOE ε4/4 and APOE ε2/3 carriers. Thresholds for two-fold change and p value < 0.05 are indicated by dashed lines. Significant changes are indicated with different colors. Grey: non-significant (Not Sig), black: non-adjusted p value < 0.05. (D) Expression level of proteins differentially expressed between APOE ε4/4 and APOE ε2/3 in 10K (left) and EVs (right). Data are presented as mean +/- SD.
The fold changes and p values of individual proteins upregulated and downregulated in ε4/4 vs ε2/3 are shown in Figure 6C (10K and EVs) and Figure S9C (BH). The expression levels of dysregulated proteins between ε4/4 vs ε2/3 are shown in Figure 6D and Figure S10A-B. As for AD pathology overall, more differentially expressed proteins were identified in 10K and EVs than in BH. The potential functions of these proteins in neurodegenerative disease progression are summarized in Tables 4 and 5.
Expression level of proteins differentially expressed in APOE ε4/4 versus APOE ε2/3 for 10K (A) and EVs (B). Data are presented as mean +/- SD.
The potential function of significantly different proteins (P value < 0.05) between APOE ε4/4 and ε2/3 carriers in EVs
The potential function of significantly different proteins (P value < 0.05) between APOE ε4/4 and ε2/3 carrier in 10K
DISCUSSION
The roles of EVs in regulating CNS diseases have been inferred predominantly from studies of in vitro models and biofluid EVs, with growing but still limited study of tissue EVs. Here, we compared protein and small RNA contents of brain homogenate with those of a “10K” pelleted extracellular fraction and a purified EV fraction of control and late-stage AD brain, including several APOE genotypes. Although the concentrations of recovered particles did not differ significantly with AD pathology, several bdEV miRNAs, tRNAs, and proteins were dysregulated not only between AD and controls, but also between the APOE genotypes representing the most “distant” risk groups in our study, the ε4/4 and ε2/3 carriers. Proteome differences between AD and controls were most pronounced for EVs, suggesting that EV proteins may have the best biomarker potential, at least for late-stage AD. The dysregulated molecules identified in our study, especially those involved in aging and neurodegeneration pathways, may be involved in CNS disease mechanisms and constitute new biomarkers for disease monitoring.
Does AD or APOE genotype affect EV production?
Consistent with previous studies on bdEVs from human cortex (91, 92), we observed a similar bdEV recovery for AD and control. However, since bdEVs are released as a mixture of various subtypes and from diverse cells, this finding does not exclude the possibility that biogenesis of specific EV subtypes may be affected by AD. For example, MHC class I bdEVs were reportedly enhanced in preclinical AD patients compared with controls and late-stage AD (92). Also, the APOE ε4 allele was reported to be involved in neuronal, endosomal, and lysosomal system dysfunction (93). In another study, control brain tissues from human and mouse had reduced EV counts in the presence of the APOE ε4 allele (61). However, in our late-stage AD samples (Braak stage 5-6, CERAD B-C), no particle concentration difference was associated with the APOE ε4 allele or other variants. Consistent with previous findings, this suggests that effects of APOE alleles, if present, may affect EVs and their effects on brain function in a measurable manner only before late-stage disease (13, 94). Further studies of APOE will be needed in young brain and during early disease stages.
Altered bdEV miRNA expression in AD
RNA sequencing revealed several differentially expressed (DE) miRNAs in AD bdEVs with potential involvement in AD. Importantly, these DE miRNAs largely overlapped with results of earlier studies on miRNAs in AD brain tissues (54–56). For example, miRNA-132-3p was downregulated in AD EVs. It was previously found to be downregulated in late-onset AD, negatively correlated with Braak stages and formation of tau tangles in the prefrontal cortex (55). miR-132 promotes neurite outgrowth and synapse formation (95, 96), with expression driven by the neurotrophin-responsive transcription factor CREB (97). Additionally, miR-132 knockout mice express increased levels of phosphorylated tau (98) and display memory deficits (99). In line with downregulated miR-132-5p, we found that tau levels were increased in AD EVs. Two other miRNAs downregulated in AD EVs (miR-485-3p and miR-338-5p) target beta-secretase 1 (BACE1) (100, 101). Together with γ-secretase, BACE1 cleaves APP into aggregate-forming Aβ peptides (102). Downregulated miR-129-5p and miR-432-5p regulate neuroinflammation (103, 104), and miR-129-3p targets fragile X mental retardation 1 (FMR1), which controls migration of cortical neurons (105, 106). Although our study predominantly highlights downregulated miRNAs, upregulated miRNAs have also been reported in AD bdEVs (107). We found at least three upregulated miRNAs in common with that study: miR-483-5p, miR-152-3p, and miR-30a-5p (107). Beyond these validations, differences between the studies may be explained by use of different brain regions for bdEV isolation, as well as differences in separation methods, possibly affecting purity and recovery of EV subtypes. Overall, we found various dysregulated miRNAs in AD EVs that appear to be involved in AD hallmarks such as APP and tau processing, neuronal function, memory formation, and neuroinflammation.
bdEV miRNAs and APOE genotype
In contrast with AD versus control comparisons, EV miRNA profile differences by APOE were not as pronounced. Nevertheless, several EV-associated miRNAs differed. In the 10K fraction of APOE ε4/4 carriers, miR-29b-5p was downregulated. This miRNA targets BACE, which contributes to Aβ aggregation (Hebert et al PNAS 2008). Additionally, we observed downregulation of miR-379-5p and miR-410-3p, which have been previously implicated in anxiety behavior (108), miR-23a-5p, which is linked to neuronal cell loss (109), and miR-483-5p, which is involved in neuronal development (110). Our findings are limited, however, since only late-stage carriers were examined, with late Braak stage and similar CNS pathology. miRNAs are quite possibly involved in APOE genotype interactions much earlier in AD, perhaps even in young, putatively healthy individuals. More research is thus needed to address the involvement of these miRNAs in AD pathogenesis.
AD-dependent differences in additional bdEV-associated ncRNAs
AD-dependent differences in bdEV-associated ncRNAs such as tRNAs and Y-RNAs were also apparent. These ncRNAs are highly abundant in EVs from various sources (84–86,111–113), but their roles in health and disease remain largely unexplored. We found various AD-dependent differences in tRNA content of EVs, as well as APOE-dependent differences in tRNAs in the 10K fraction. Changes in cellular tRNA abundance have been implicated in neurological disorders (114, 115), and dysregulation of cellular tRNA levels affects protein expression in cancer (116). Interestingly, tRNA fragments released with EVs may modulate gene expression in recipient cells (117). Small changes in bdEV-associated Y-RNAs (Y1 and Y4) were also observed. Y-RNA is abundantly detected in EVs from various biological sources (118). Cellular Y-RNA serves as a scaffold for RNA-binding proteins (119). In AD patients, dysregulation of Y3-RNA binding to enhancer protein HuD (ELAVL4) was shown to cause alternative splicing in neurons (120). Additionally, EV-associated Y4-RNA activated endosomal RNA sensor TLR7 in monocytes (90). In brain tissue of AD patients, we observed a significant difference in U4 RNA, which has not been previously reported in AD. As part of the spliceosome complex, U4 is involved in pre-mRNA splicing (121). In the prefrontal cortex of AD brains, snRNA U1 was found to be upregulated and to accumulate in tau plaques (122, 123). Upregulation of U1 was additionally implicated in aberrant mRNA splicing and altered APP expression (122). Our findings encourage further investigations into the role of extracellular non-miRNA ncRNAs in AD and other neurological conditions.
Implications of bdEV-associated proteins in AD
bdEV-associated tau proteins
EVs have been widely reported to propagate misfolded protein MAPT (tau) in neurodegenerative disease (42,124–126). Consistent with other bdEV studies (42,91,125), we found a higher level of tau and phosTau protein in EVs of AD brain compared with controls, adding more evidence to EV-mediated tau propagation. However, we did not see consistent associations with APOE genotype, possibly because our samples sizes were small, and possibly because we examined brains from individuals with late-stage disease. Further study is needed, including study of younger and early-stage individuals with different genotypes, to decide if EVs act as a protective factor to clear tau from the CNS to the periphery, or as a neurotoxic factor to spread tau between neurons and glia, or perhaps some as-yet incompletely understood balance between the two extremes.
bdEV-associated antioxidant proteins
Upregulated members of the antioxidant peroxiredoxin (PRDX) family in AD EVs include PRDX1, PRDX2, and PRDX6. Oxidative damage is often implicated in the progression of neurodegenerative diseases. We and our collaborators previously found that PRDXs were highly abundant in EVs derived from both human induced pluripotent stem cells (iPSCs) and mesenchymal stem cells (MSCs), and that these proteins help to reduce oxidative stress in induced senescent cells (127). Although the protective role of PRDX family members in AD (128–130) and the association of PRDX with EVs (127,131,132) have been reported before, this is the first study to show PRDXs upregulating in tissue EVs from AD brain. We also found another antioxidant defense protein, deglucase DJ1 (PARK7), enriched in APOE ε4/4 compared with ε2/3 samples. PARK protects cells against hydrogen peroxide-induced death (133, 134) upon oxidative stress. We hypothesize that this factor may spread via EVs as a neuroprotective response in patients with the risk factor ε4.
bdEV-associated mitochondrial and metabolism pathway proteins
Mounting evidence suggests that mitochondrial and metabolism dysfunction contribute to AD by compromising the energy supply (135–138). We found numerous bdEV-associated proteins that contribute to mitochondrial and metabolic pathways. The ATPase enzyme member valosin-containing protein/p97 (VCP) was upregulated in AD bdEVs. VCP is thought to bring about mitophagy impairment, leading to mitochondrial dysfunction in AD (139, 140). Higher levels of tricarboxylic acid cycle (TCA) protein mitochondrial aconitate hydratase (ACO2) were previously found in multiple brain regions of AD patients (141) and here appeared to be enriched in bdEVs of APOE ε4/4 patients compared with ε2/3. We also identified changes in bdEV proteins that are involved in glycolytic and carbon metabolism, including transketolase (TKT) (142, 143), alpha-enolase (ENO1) (144, 145), and fructose-bisphosphate aldolase A (ALDOA) (146). Further exploration of protein modification and enzymatic activity of these proteins is now indicated to unveil mechanisms of energy metabolism deficits in AD.
bdEV-associated proteins in neurogenesis
Impairment of pathways involved in neurogenesis/differentiation are thought to be early clinical events in aging and AD development (136,147–150). Among neurogenesis-related proteins, the 14-3-3 protein family of ubiquitous phosphoserine/threonine–binding proteins is highly abundant in brain and implicated in nervous system development and maintenance (151–153). Our study shows dysregulation of 14-3-3 isoforms ε (YWHAE), ζ (YWHAZ), and γ (YWHAG) in the 10K and EV fractions. 14-3-3 proteins indirectly regulate activation of Rho family GTPases (154). Interestingly, several proteins in the Rho-GTPase cycle were also dysregulated in our bdEVs, including Ras-related C3 botulinum toxin substrate 1/3 (RAC1/3), and Rho GDP-dissociation inhibitors 1 and 2 (ARHGDIA, GDI2). We also found that contactin-1 (CNTN1) and neurofascin (NFASC), two proteins involved in axonal guidance and neuron projection development (155–158), were downregulated in AD bdEVs, which may reflect neurodegeneration. In addition, some neuroprotective proteins were also found to be downregulated in bdEVs from APOE ε4/4 patients compared to ε2/3, such as γ-enolase (ENO2) (159–161) and brevican core protein (BCAN) (162–164).
In summary, this is the first study thoroughly examining bdEV small RNA and protein content from clinical AD patients with different APOE genotypes. Our results identified the composition of RNAs (including miRNAs, tRNAs, and Y-RNAs) and proteins in bdEVs that varied with AD pathology and APOE genotypes. These bdEV-transported molecules may play critical roles in modulating AD progression and thus have potential for exploitation as biomarkers and therapeutic targets.
Data Availability
Nucleic acid sequencing data have been deposited with the Gene Expression Omnibus, accession GSE159541. We have submitted all relevant details of our experiments to the EV-TRACK knowledgebase (EV-TRACK ID: EV200126). Any other data are available on reasonable request.
https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE159541
CONFLICTS OF INTEREST
The authors report no conflicts of interest.
ACKNOWLEDGMENTS
This work was supported in part by grants from the US National Institutes of Health: AI144997 (to KWW, with support for TAPD), MH118164 and AG057430 (to VM and KWW), by UG3CA241694 (to KWW), supported by the NIH Common Fund, through the Office of Strategic Coordination/Office of the NIH Director, and by the National Health and Medical Research Council of Australia (GNT1132604 to AFH). The Johns Hopkins Alzheimer’s Disease Research Center provided tissues and is funded by P30 AG066507 (PI: M. Albert). VM and KWW gratefully acknowledge support from the Richman Family Precision Medicine Center of Excellence in Alzheimer’s Disease including helpful comments and advice from founder and director Constantine Lyketsos. Thanks to: Kenneth Pienta, Johns Hopkins University School of Medicine, for access to the nanoFCM flow nanoAnalyzer platform; Mitchell Science Writing for manuscript editing and citation formatting; and the La Trobe University Comprehensive Proteomics Platform. We also thank members of the Witwer and Retrovirus Laboratories, Johns Hopkins University School of Medicine, and various members of the International Society for Extracellular Vesicles for valuable discussions and support.
Footnotes
Email addresses: Yiyao Huang yiyaohuang{at}jhu.edu, Tom A. P. Driedonks tdriedo1{at}jhu.edu, Lesley Cheng l.cheng{at}latrobe.edu.au, Andrey Turchinovich andrey.turchinovich{at}sciberg.com, Harinda Rajapaksha K.Rajapaksha{at}latrobe.edu.au, Tanina Arab tarab1{at}jhu.edu, Bonita H. Powell bpowel15{at}jhu.edu, Olga Pletniková olgaplet{at}buffalo.edu, Javier Redding redding54{at}jhmi.edu, Juan C. Troncoso troncoso{at}jhmi.edu, Laura Vella laura.vella{at}florey.edu.au, Lei Zheng nfyyzhenglei{at}smu.edu.cn, Andrew F. Hill Andrew.Hill{at}latrobe.edu.au, Vasiliki Mahairaki vmachai1{at}jhmi.edu, Kenneth W. Witwer kwitwer1{at}jhmi.edu
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.
- 6.
- 7.
- 8.↵
- 9.↵
- 10.
- 11.↵
- 12.↵
- 13.↵
- 14.
- 15.
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.
- 30.↵
- 31.↵
- 32.↵
- 33.
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.
- 39.↵
- 40.
- 41.
- 42.↵
- 43.↵
- 44.↵
- 45.
- 46.
- 47.
- 48.↵
- 49.↵
- 50.
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵
- 120.↵
- 121.↵
- 122.↵
- 123.↵
- 124.↵
- 125.↵
- 126.↵
- 127.↵
- 128.↵
- 129.
- 130.↵
- 131.↵
- 132.↵
- 133.↵
- 134.↵
- 135.↵
- 136.↵
- 137.
- 138.↵
- 139.↵
- 140.↵
- 141.↵
- 142.↵
- 143.↵
- 144.↵
- 145.↵
- 146.↵
- 147.↵
- 148.
- 149.
- 150.↵
- 151.↵
- 152.
- 153.↵
- 154.↵
- 155.↵
- 156.
- 157.
- 158.↵
- 159.↵
- 160.
- 161.↵
- 162.↵
- 163.
- 164.↵