- 1Department of Neurology, Harvard Medical School, MassGeneral Institute for Neurodegenerative Disease, Massachusetts General Hospital, Charlestown, MA, United States
- 2Department of Pathology, Harvard Medical School, MassGeneral Institute for Neurodegenerative Disease, Massachusetts General Hospital, Charlestown, MA, United States
Alzheimer's Disease (AD) is associated with neuropathological changes, including aggregation of tau neurofibrillary tangles (NFTs) and amyloid-beta plaques. Mounting evidence indicates that vascular dysfunction also plays a key role in the pathogenesis and progression of AD, in part through endothelial dysfunction. Based on findings in animal models that tau pathology induces vascular abnormalities and cellular senescence, we hypothesized that tau pathology in the human AD brain leads to vascular senescence. To explore this hypothesis, we isolated intact microvessels from the dorsolateral prefrontal cortex (PFC, BA9) from 16 subjects with advanced Braak stages (Braak V/VI, B3) and 12 control subjects (Braak 0/I/II, B1), and quantified expression of 42 genes associated with senescence, cell adhesion, and various endothelial cell functions. Genes associated with endothelial senescence and leukocyte adhesion, including SERPINE1 (PAI-1), CXCL8 (IL8), CXCL1, CXCL2, ICAM-2, and TIE1, were significantly upregulated in B3 microvessels after adjusting for sex and cerebrovascular pathology. In particular, the senescence-associated secretory phenotype genes SERPINE1 and CXCL8 were upregulated by more than 2-fold in B3 microvessels after adjusting for sex, cerebrovascular pathology, and age at death. Protein quantification data from longitudinal plasma samples for a subset of 13 (n = 9 B3, n = 4 B1) subjects showed no significant differences in plasma senescence or adhesion-associated protein levels, suggesting that these changes were not associated with systemic vascular alterations. Future investigations of senescence biomarkers in both the peripheral and cortical vasculature could further elucidate links between tau pathology and vascular changes in human AD.
Introduction
Vascular dysfunction has become increasingly implicated in the pathogenesis of Alzheimer's Disease (AD) (1, 2). Cerebrovascular diseases including cerebral amyloid angiopathy (CAA) and atherosclerosis often co-present with AD (3–5), and reduced cerebral blood flow in the human AD brain is associated with both cognitive decline (6, 7) and tau accumulation (8–10). Proper function of the neurovascular unit, which is comprised primarily of endothelial cells, mural cells, and astrocytic endfeet, is essential for nutrient supply and protection from peripheral blood molecules via the blood-brain barrier (BBB). Impaired cerebral blood flow represents a key mechanism in neurodegeneration as it contributes to neuronal injury (11) and the accumulation of amyloid beta aggregates (12).
One of the neuropathological hallmarks of AD is the aggregation of phosphorylated tau protein into neurofibrillary tangles (NFTs). Several lines of research indicate that tau pathology may interact with endothelial cell changes to drive vascular impairment in AD. Senescence is a complex process in which cells irreversibly cease proliferating, change morphology, and present the senescence-associated secretory phenotype (SASP) (13). In neurons, cellular senescence and tau NFT formation have been reported to exhibit a positive feedback relationship (14, 15), and suppression of senescence prevents tau aggregation in vivo (15). Moreover, endothelial senescence is associated with vascular dysfunction via BBB disruption (16, 17), increased vascular stiffness (18, 19), and aberrant neurovascular coupling (10, 20–22).
The SASP profile includes pro-inflammatory cytokines, cell cycle regulators, and pro-angiogenic factors (13, 23, 24), which propagate senescence via autocrine and paracrine signaling (25, 26). This profile phenocopies, to some extent, changes we have observed in the microvasculature of tau-overexpressing transgenic (Tg4510) mice (27). We therefore examined whether endothelial senescence contributes to microvascular alterations in tau-related neurodegenerative disease. To investigate this hypothesis, we isolated intact microvessels from the dorsolateral prefrontal cortex (PFC, Brodmann area 9) of 16 AD and 12 control cases and measured expression of 42 genes related to senescence, cell adhesion, and general endothelial cell functions. Since many of the senescence-associated biomarkers are secreted, we also investigated whether changes in plasma protein levels could be detected in AD subjects in the years before death.
Materials and Methods
Subject Information
Twenty-eight cases were selected from the Neuropathology Core of the Massachusetts Alzheimer's Disease Research Center. Thirteen of the subjects also had longitudinal plasma sample collection over several years; these cases are indicated with an asterisk in the Plasma column in Table 1. All cases were assessed by a neuropathologist and scored for Alzheimer's Disease (AD) neuropathology burden according to the NIA-AA guidelines (28). Control cases (n = 12) were defined as subjects with a Braak neurofibrillary tangle (NFT) score of 0/I/II (B1) and were classified as low AD probability. One of the B1 subjects exhibited Lewy Body Disease with Primary Age-Related Tauopathy and one exhibited Progressive Supranuclear Palsy pathology, as labeled in Table 1. AD cases (n = 16) were defined as subjects with a Braak NFT score of V/VI (B3) indicating extensive tau pathology in the neocortex, with NIA-AA classification as intermediate or high AD probability. B3 subjects did not exhibit Lewy bodies, though 6 out of 16 had evidence of TDP-43 inclusions in the amygdala and hippocampus.
Cerebral Microvessel Isolation
We adapted this protocol for intact cortical microvessel isolation from Boulay et al. (29) and Hartz et al. (30), which is designed to protect the intercellular connections between endothelial cells and mural cells. All steps for microvessel isolation were performed on ice under RNAse-free conditions. Approximately 200 mg of frozen postmortem dorsolateral prefrontal cortex (PFC, Brodmann area 9) tissue was measured by a histologist for each case. Meninges were removed from the cortical surface and the cortical tissue was sliced into ~2 mm sections using a sterile razor blade in Hank's Balanced Salt Solution (HBSS) and HEPES buffer. The tissue sections were manually dissociated with a dounce homogenizer and centrifuged for 10 min at 2,000 g in 4°C. The pellet was resuspended in a dextran solution and centrifuged at 4,400 g for 15 min at 4°C. A myelin layer formed atop the supernatant, which was removed by inverting the centrifuge tube and carefully blotting the inside walls. The pellet was resuspended in the HBSS-HEPES buffer with RNAse-free bovine serum albumin (BSA) and filtered over a 20 μm mesh filter (Millipore) to collect microvessels. Microvessels were collected from the filter by gently stirring the filter in HBSS-HEPES-BSA buffer. The collected microvessels were washed and spun twice at 2,000 g for 5 min at 4°C. The supernatant was discarded and the pellet containing isolated microvessels was stored at −80°C prior to RNA isolation.
RNA Isolation and RT-qPCR
Frozen microvessel pellets or 25 mg of whole tissue was resuspended in RLT buffer (Qiagen) and sonicated for 20 pulses at 10% amplitude to lyse cells. The lysed cells were then centrifuged at 13,000 rpm for 3 min to pellet out residual cellular debris. RNA was extracted using the RNeasy Mini Kit (Qiagen) and was eluted in 30 μl of nuclease-free water. RNA was assessed using the NanoDrop spectrophotometer (ThermoFisher) and diluted to a standard concentration of 20 ng/μl for each sample. cDNA was synthesized using the RT2 HT First Strand kit (Qiagen) and was combined with RT2 SYBR Green Mastermix fluorescent dye (Qiagen). The cDNA reaction mix was loaded into a custom 96-well RT2 Profiler Array with pre-loaded primers. Each plate contained 42 genes of interest (Table 2), three housekeeping genes, a Human Genomic DNA Control (HGDC), a Reverse Transcription Control (RTC), and a Positive PCR Control (PPC). Two samples were run per plate. The plate was covered with a plastic seal to prevent sample evaporation and briefly spun at ~400 rpm to remove bubbles. The qPCR reactions were performed using the BioRad CFX96 Real-Time Detection System. The sequence began with a 10-min incubation at 95°C to activate the DNA polymerase enzyme. Fluorescence data collection then commenced with 40 cycles of alternating 15 s at 95°C and 60 s at 60°C.
The ΔΔCT relative quantification method was used to calculate gene expression (59). First, positive controls for reverse transcription and PCR (RTC and PPC) were analyzed for each sample according to manufacturer instructions. Gene Ct (cycle threshold) results were quality filtered to remove any data points with a Ct value greater than the sample Ct for the genomic DNA (HGDC) control primer. Any reaction with no fluorescence detection by 40 qPCR cycles was reported as not detected, and any gene with undetectable expression in more than 25% of samples in total was excluded. Three reference genes were included in the RT2-Profiler array: ACTB, GAPDH, and PGK1. The geometric mean of the reference gene Ct values was calculated for each subject and subtracted from the Ct value of each target gene, yielding ΔCT values. For each gene, the average ΔCT value in B1 subjects was then subtracted from each sample's ΔCT value to obtain ΔΔCT values. The relative quantification (RQ) value was calculated as 2−ΔΔCT, which was then log2-transformed to yield fold changes.
Cell-type specific genes were included to verify the presence of vascular cells in microvessel isolates and compared to bulk RNA measures from total cortex. Using the endothelial cell marker PECAM1 (CD31) and the vascular smooth muscle cell marker ACTA2, we confirmed these cells are significantly enriched in isolated microvessels compared to total cortex homogenate in 16 representative samples (Supplementary Figure 1A). Astrocytes (ALDH1L1) were not enriched in isolated microvessels compared with total cortex, and microglial (ITGAM) cell composition was not statistically different. There was also no statistical difference in cell composition between B1 vs. B3 microvessels when examining all n = 28 microvessel samples (Supplementary Figure 1B). Overall, endothelial cell and vascular smooth muscle cell expression was greater than astrocyte and microglia expression in isolated microvessels.
Plasma Proteomics
For 13 of the subjects (n = 4 B1, n = 9 B3), plasma was collected at two or three time points spaced at ~1-year intervals (Table 1, labeled with asterisks). Plasma was sent to Olink Proteomics for protein quantification using 92-plex antibody labeling with a proximity extension assay. Quantification results were normalized and log-2 transformed and were reported as Normalized Protein eXpression (NPX) values. Twenty-two of the biomarkers analyzed via RT-qPCR were also included in the Olink protein panels (Table 2). Quality control was applied by eliminating any data points with an NPX value below the manufacturer-specified protein limit of detection (LOD) threshold.
Statistical Analysis
All statistical analyses and data visualizations were performed with R (v4.0.0). A significance level was set at p < 0.05, and correction for multiple comparisons was applied where necessary using the Benjamini-Hochberg (BH) False Discovery Rate (FDR) (60). This was calculated with the p.adjust() function in R using the “stats” package (v4.0.0).
Univariate comparisons were conducted using Welch's two-way unpaired t-test using the t.test() function from the R “stats” package. To test the robustness of a given result, bootstrapped iterations were performed in which 75% (12/16) of B3 and 83% (10/12) of B1 subjects were randomly sampled 1,000 times, and the sample estimates were computed (61). The average fold change difference and proportion of significant iterations (p < 0.05) were reported to summarize bootstrapped iteration results. Regression models were implemented using the lm() function from the R “stats” package and results were reported using the tidy.lm() function from the “broom” R package (v0.5.6). ANOVA and Tukey's post-hoc HSD test were implemented using the aov() functions from the R “stats” package and the tukey_hsd() function from “rstatix” (v0.6.0). Note: Missing data were excluded on a gene-wise basis for t-test and regression analyses, though imputation with the highest Ct value per gene was also examined (62).
Principal component analysis (PCA) was implemented using the “FactoMineR” (v2.3) and “factoextra” (v1.0.7) packages in R. As PCA requires complete data, missing data were imputed on a gene-wise basis with the highest observed Ct value for the corresponding gene (62). Fold change data were centered to have a mean of zero before PCA was applied using the prcomp() function in FactoMineR. The proportion of variance per principal component (PC) as well as cumulative variance were examined with a Scree plot to determine how many components would be further analyzed.
Results
Senescence-Associated Genes Are Upregulated in B3 PFC Microvessels
We isolated intact microvessels from the dorsolateral prefrontal cortex (PFC, BA9) to examine gene expression changes related to the extent of AD-related tau pathology. Enrichment of vascular transcripts in microvessel preparations was first verified by RT-qPCR (see Methods, Supplementary Figure 1). The PFC was selected as a Braak III/IV (B2) region, meaning it does not generally exhibit tau NFT pathology in early-stage AD (63). Braak 0/I/II (B1) subjects are therefore unlikely to present any PFC tau pathology, while Braak V/VI (B3) subjects are likely to have a high tau NFT burden in the PFC (64).
To assess gene expression associated with senescence, cell adhesion, and endothelial cell function, we measured expression of 42 genes classified into eight hypothesis-driven functional groups (Table 2). These eight classes were proposed based on literature reviewed in Table 2. Of the 42 genes measured, CDKN2A and SELP were excluded due to undetectable expression in more than 25% of samples (defined as no fluorescence detection by 40 cycles or a Ct value above the Human Genomic DNA Contamination Ct value) (Supplementary Figure 2). For the remaining 40 genes, heatmap visualization of log2 fold changes showed an observable hotspot of senescence-associated gene upregulation in B3 PFC microvessels (Figure 1). The average functional group fold change value was calculated per subject (excluding the non-endothelial cell markers) which confirmed a 2.5-fold increase in senescence-associated gene mRNA in B3 microvessels (Figure 2A; p = 0.0030, BH-FDR = 0.0208, 95% CI: [0.5019, 2.1442]). Additionally, expression of cell adhesion genes was increased by 1.9-fold in B3 microvessels compared to B1 microvessels (Figure 2A; p = 0.0161, BH-FDR=0.0562, 95% CI: [0.1945, 1.7123]). Missing values were omitted from these mean fold change calculations, though very similar results were obtained by imputing missing values with the highest Ct value of the gene. The other functional groups examined did not exhibit differential expression in B1 vs. B3 PFC microvessels (Figure 2B).
Figure 1. Senescence-associated genes form a fold change hotspot in prefrontal cortex microvessels isolated from high-Braak (B3) vs. low-Braak (B1) samples. Differential gene expression between microvessels isolated from B3 vs. B1 samples, grouped by hypothesis-driven gene function groups. Values reflect log2-transformed fold change values relative to the B1 subject mean expression per gene. Relative quantification was performed such that the mean of B1 fold changes is zero for each gene. For this visualization, missing data were imputed with the ABC-Braak group mean fold change for the gene.
Figure 2. Senescence- and adhesion-associated genes are significantly upregulated in prefrontal cortex microvessels isolated from B3 vs. B1 subjects. (A) Senescence-associated genes are significantly upregulated in B3 cerebral microvessels with a log2 fold change of 1.323; Welch's unpaired t-test, p = 0.0030 (BH-corrected FDR 0.0208), 95% CI for B3 log2 fold change: 0.5019, 2.1442. Cell adhesion genes are also significantly upregulated with a log2 fold change of 0.953; Welch's unpaired t-test, p = 0.0161 (BH-FDR 0.0562), 95% CI for B3 log2 fold change: 0.1945, 1.7123. (B) Volcano plot showing t-test results. Dashed vertical lines denote magnitudes of 2-fold change and the dashed horizontal line denotes the cutoff for p < 0.05. (C) Stepwise multiple regression models were applied to senescence-associated and cell adhesion functional groups to measure the association between Braak stage and functional group average fold change. The first row indicates the independent variables upon which the average fold change was regressed per model for each functional group. Within each model, the ABC-Braak term regression coefficient (β1), p-value, and BH-FDR are reported. Significant ABC-Braak β1 terms are denoted with **p < 0.01, *p < 0.05.
Sex, cerebrovascular pathology, and age can contribute to AD neurodegeneration, both directly and via tau pathology (65). To adjust for these important and potentially confounding covariates, stepwise multiple regression models were implemented to compare senescence-associated and cell adhesion mean fold changes between B3 vs. B1 samples (Figure 2C). We created the binary variable “Cerebrovascular pathology” to indicate presence or absence of at least one of cerebral amyloid angiopathy (CAA), cerebrovascular disease (CVD), and infarcts in the subject's neuropathology report (Table 1). Senescence-associated gene expression was still significantly increased in B3 microvessels compared to B1 after adjusting for sex, cerebrovascular pathology, and age. Cell adhesion gene expression was significantly increased in B3 microvessels after adjusting for sex and cerebrovascular pathology but not age at death, indicating that microvascular adhesion gene expression may be related to age. We note there was no statistical difference in the age at death between B1 and B3 subjects (Supplementary Figure 3).
IL-8, PAI-1, and TIE1 Demonstrate Robust Upregulation in B3 PFC Microvessels Independent of Age, Sex, and Cerebrovascular Pathology
After observing the general upregulation of senescence-associated and cell adhesion functional group genes in B3 PFC microvessels, we next investigated changes in individual genes. Of the 40 genes that we measured in PFC microvessels which passed quality control filtering, 11 were significantly upregulated (with p < 0.05, BH-FDR <0.20): CXCL8 (IL8), SERPINE1 (PAI-1), CXCL1, CXCL2, CSF2, CDKN1A, ICAM-2, MAdCAM-1, SELE, TIE1, and DLL4 (Figure 3A). With the exception of CDKN1A, all of the senescence-associated genes in this group exhibited greater than a 2-fold average increase in mRNA expression in B3 microvessels. Further, these results were robust for CXCL8 (IL-8), SERPINE1 (PAI-1), CXCL1, CXCL2, TIE1, and ICAM-2, which were significantly upregulated (p < 0.05) in 62–94% of 1,000 bootstrapped iterations (Supplementary Figure 4). Missing values were excluded on a gene-wise basis for these comparisons, though these six genes were still significantly upregulated after imputing missing values with the highest Ct value of the gene.
Figure 3. Eleven genes related to senescence, adhesion, endothelial cell function, and VEGF/NOTCH pathways are significantly upregulated in PFC capillaries from B3 samples relative to B1 samples. (A) 11 genes show significant upregulation in B3 vs. B1 microvessels. Welch's unpaired t-test results are shown: **p < 0.01, *p < 0.05, all BH-FDR < 0.2. Plots are shaded by hypothesis-driven functional group. (B) Stepwise multiple regression models were applied to each of the 6 most-consistently upregulated genes to measure the association between Braak change and log2 fold change, adjusting for potential confounders. The first row indicates the independent variables upon which the log2 fold change was regressed per model for each gene. Within each model, the ABC-Braak term regression coefficient (β1), p-value, and BH-FDR are reported. Significant ABC-Braak β1 terms are marked in bold with **p < 0.01, *p < 0.05.
For each of these six robustly upregulated genes, individual linear models were fit to regress gene fold change on ABC-Braak score, with stepwise inclusion of sex, cerebrovascular pathology, and age as covariates. The increased expression in B3 samples remained significant for all six genes after adjusting for sex and cerebrovascular pathology, indicating their upregulation was independent of these factors (Figure 3B). However, with the inclusion of age at death in the regression model, only CXCL8 (IL-8), SERPINE1 (PAI-1), and TIE1 were still significantly upregulated in B3 microvessels, suggesting their upregulation was age-independent. For the other three genes (ICAM-2, CXCL1, CXCL2), elevated expression in B3 PFC microvessels may be associated with age.
Tau Pathology and Cerebrovascular Pathology Interact to Increase Microvessel MAdCAM-1 Expression
Approximately half of both the B1 and B3 subjects exhibited at least one form of cerebrovascular pathology (Table 1), which is consistent with population-based studies in elderly demented and non-demented individuals (66, 67). Given the documented association between tau pathology and cerebrovascular pathology (68–70), we investigated whether the two interact in association with microvessel gene expression. Two-way ANOVA revealed an interaction between tau pathology and cerebrovascular pathology in MAdCAM-1 expression in PFC microvessels, which was significant prior to adjusting for multiple comparisons (Supplementary Figure 5A). Post-hoc analysis via Tukey's HSD test revealed that MAdCAM-1 mRNA expression was upregulated by more than 16-fold (4 log2-fold) in B3 subjects with cerebrovascular pathology compared to all other subjects (Figure 4, Supplementary Figure 5B). Of note, the top three MAdCAM-1 expression values correspond to B3 subjects with concomitant cerebral amyloid angiopathy (CAA) and general cerebrovascular disease (CVD). However, we note that the two-way ANOVA was not significant for MAdCAM-1 when missing values were imputed with the highest MAdCAM-1 Ct value.
Figure 4. MAdCAM-1 upregulation demonstrates interaction between Braak stage and presence of cerebrovascular pathology. ANOVA with interaction demonstrated significant interaction between a high Braak score (B3) and presence of cerebrovascular pathology in MAdCAM1 expression out of all 40 genes quantified. Data shown in the boxplot reflect Tukey's HSD post-hoc test results, with **p < 0.01, *p < 0.05. A total of 6 samples did not have detectable expression, with an even distribution of missing data across the four groups. The number of samples with detected MAdCAM-1 expression per group are indicated at the bottom of the boxes.
Endothelial Senescence-Associated and Cell Adhesion Gene-Driven Principal Component Composite Score Distinguishes B3 vs. B1 Microvessels
To explore patterns of gene covariance in B1 vs. B3 cortical microvessels, principal component (PC) analysis was applied to centered fold change values for all 40 genes. The first two PCs collectively explained 51.2% of variance across the samples (Supplementary Figure 6). The top contributors to the first principal component (PC1), defined as those with more than one standard deviation above the mean contribution, included senescence-associated genes (IL6, CXCL8, CSF2) and cell adhesion genes (SELE and MAdCAM-1) (Figure 5A). As these genes are associated with endothelial dysfunction and/or senescence (23, 24, 31, 71–73), we hypothesized that PC-derived composite gene expression scores (a.k.a. PC composite score) would distinguish B3 from B1 microvessels. Indeed, the PC scores from the first and second principal components show partial separation of the B1 and B3 sample clusters, primarily along the PC1 axis (Figure 5B). This difference was statistically significant for the PC1 composite score, with B3 samples exhibiting significantly larger composite scores (Welch's unpaired t-test, p = 0.0064, 95% CI: [1.477, 7.956]; Figure 5C).
Figure 5. Principal component composite gene scores in B3 vs. B1 microvessels show variance driven by senescence and adhesion genes. (A) PC loadings are plotted for each gene in PC1 (x-axis) and PC2 (y-axis). The rightmost shaded region contains the genes considered to be top contributors to PC1, with their percent contribution at least one standard deviation above the mean contribution. Genes immediately adjacent to this top contributor subset are also labeled. Point and text colors correspond to each gene's hypothesis-driven functional group, listed on the right. (B) Samples are plotted by PC1 and PC2 composite gene fold change scores, with geometric encircling to outline B1 vs. B3 group boundaries. (C) Boxplots show PC1 composite gene fold change scores in B1 vs. B3 PFC microvessels. Welch's unpaired t-test was applied to compare B1 and B3 composite scores, with ** indicating p < 0.01. Composite gene fold change scores are significantly larger in B3 vs. B1 microvessels (p = 0.0064, 95% CI: [1.477, 7.956]).
Four of the five genes driving fold change composite scores—CXCL8 (IL8), CSF2, SELE, and MAdCAM-1—were identified as upregulated in B3 microvessels, with CXCL8 (IL8) upregulation particularly pronounced after adjusting for age, sex, and cerebrovascular pathology (Figure 3). By contrast, IL6 upregulation did not reach significance in B3 microvessels. Its strong contribution to the gene composite score that separated B1 vs. B3 microvessels may indicate that IL-6 is associated with ABC-Braak differences in the context of other senescence-associated and cell adhesion genes. Of note, IL6 expression differences were not driven by any other neuropathological variable, sex, APOE genotype, or cerebrovascular pathology (Supplementary Figure 7). Furthermore, PCA composite scores did not improve separation B1 vs. B3 microvessels based on these other covariates relative to ABC-Braak scores (Supplementary Figure 8).
NOS3 Antemortem Plasma Protein Expression Is Associated With Postmortem Cortical Microvessel Gene Expression
The cerebral vasculature interacts with the peripheral vasculature via plasma proteins and soluble blood factors (74). Therefore, we compared postmortem PFC microvessel gene expression with antemortem plasma protein expression in a subset of 13 subjects with both metrics available (Table 1, subjects with asterisks). In total, 22 proteins overlapped with cognate genes measured in our qPCR assay (Table 2), though SELP was not analyzed as it was excluded from qPCR analysis due to a high proportion of missing data. Protein expression was quantified as Normalized Protein eXpression (NPX), a log2-transformed unit that enables inter-subject comparison within a given protein. Of note, two ITGAM values were omitted as they were below the manufacturer-specified limit of detection (see Methods). Protein expression was measured longitudinally in two or three plasma samples per subject, though we focused here on the final plasma samples to minimize pre-mortem intervals.
Univariate regression revealed that NOS3 postmortem microvessel expression was significantly associated with antemortem plasma expression (Figure 6A; NPX β1 = 1.022, p = 0.0024, BH-FDR = 0.0511, R2 = 0.58). However, the pre-mortem interval ranged from 1 month to 6.7 years among samples. To account for this variability, we included pre-mortem interval as a covariate in a subsequent multiple regression. This still yielded a significant linear relationship between NOS3 plasma NPX and PFC microvessel gene fold change (NPX β1 = 1.070, p = 0.0031, BH-FDR = 0.0660). For visualization, the gene fold change values were adjusted for the partial residual of the model, which demonstrated a strong linear association between antemortem plasma NPX and postmortem gene expression (R2 = 0.94, Figure 6B).
Figure 6. NOS3 antemortem plasma expression is significantly associated with postmortem NOS3 gene expression in PFC capillaries of both B1 and B3 subjects. (A) Univariate regression of NOS3 plasma NPX (from final plasma sample per subject) vs. postmortem microvessel gene log2 fold change demonstrated a significant association between NPX and log2 fold change (β1 =1.022, p = 0.0024, BH-FDR = 0.0511). The shaded region indicates the 95% confidence interval of the regression slope. (B) After adjusting for years until death in multiple regression, there is still a significant association between NOS3 plasma NPX and postmortem capillary gene log2 fold change (β1 = 1.070, p = 0.0031, BH-FDR = 0.0660). For visualization, the adjusted log2 fold change values were calculated as the NPX β1 coefficient plus the partial residual associated between plasma NPX and years until death.
Plasma Senescence-Associated Protein Levels Are Not Associated With Postmortem Severity of AD-Related Neuropathological Changes
Six of the genes that were upregulated in B3 microvessels were also quantified in plasma as secreted proteins: IL8 (CXCL8), PAI-1 (SERPINE1), CXCL1, TIE1, ICAM-2, and SELE. We reasoned that the changes seen in the CNS microvessels might reflect systemic vascular factors, which might then be observed in the plasma; alternatively, CNS-only microvascular changes may be difficult to detect in the systemic circulation. We therefore investigated whether these biomarkers exhibited protein expression differences in B1 vs. B3 antemortem plasma. Univariate analysis comparing the protein NPX from the final plasma sample between B1 vs. B3 samples failed to yield any significant results (Figure 7A) as measured via t-test and linear regression. Similarly, multivariate analysis adjusting for sex, age at visit, and pre-mortem interval yielded no significant differences in protein NPX by Braak score. While the magnitude of NPX measurements cannot be directly compared between proteins, their relative changes from baseline can be compared to assess temporal stability. Plasma time points differed between subjects, with an average of 1.30 years (±0.38) between plasma samples per subject. To that end, PAI-1 levels exhibited the greatest fluctuations from baseline across all subjects, followed by CXCL1 (Figures 7B,C). By contrast, TIE1 plasma levels remained the most stable over time across all subjects.
Figure 7. B1 vs. B3 subject plasma expression of protein biomarkers corresponding to genes upregulated in B3 cortical PFC microvessels. (A) NPX values from each subject's last plasma sample for CXCL1, ICAM-2, IL8, PAI-1, SELE, and TIE1. B1 samples are shown in blue and B3 samples are shown in red. Welch's unpaired t-test p-values are shown above the boxes. (B) Plasma NPX values across three longitudinal samples per subject for CXCL1, ICAM-2, IL8, PAI-1, SELE, and TIE1. B1 samples are shown in blue and B3 samples are shown in red. Squares represent the ABC-Braak group mean NPX value at the corresponding plasma sample time point (which is an arbitrary unit, relative to each subject). (C) Relative plasma NPX change from baseline for CXCL1, ICAM-2, IL8, PAI-1, SELE, and TIE1. Values on the y-axis represent net change in plasma NPX from the first plasma sample, which is set to zero for all subjects for visualization. Squares represent the ABC-Braak group mean change in NPX at the corresponding plasma sample time point. Scales are fixed for all six proteins to compare relative temporal stability.
None of the other 15 plasma proteins assayed in this study showed statistically different plasma levels in B1 vs. B3 subjects. Further descriptive and exploratory analyses were carried out to examine effects of age and pre-mortem interval duration on plasma protein levels. In a multiple regression of plasma protein NPX on ABC-Braak, age at visit, and pre-mortem interval duration, IL6, IL8, and ITGAM exhibited negative slopes between pre-mortem interval vs. plasma NPX before adjusting for multiple comparisons (Table 3). This indicates that subjects with longer pre-mortem intervals generally had lower IL6, IL8, and ITGAM plasma protein levels. IL6 also exhibited a positive slope between age vs. plasma NPX, as did CXCL1, suggesting that IL6 and CXCL1 plasma secretion may increase with age.
Discussion
Tau pathology leads to senescence-associated transcriptomic changes (14, 15, 27) and subsequent vascular alterations in the mouse cortex, including transient capillary occlusion (27). Such intermittent episodes of capillary blockage arise from leukocyte adhesion in the endothelium, a phenomenon also reported in an APP-PS1 mutant mouse model of AD (75), as well as in chronic hypoperfusion (76) and traumatic brain injury (77, 78). This may be related to imbalanced capillary transit in AD, which can impair neuronal oxygenation even in the absence of gross reductions in cerebral blood flow (79, 80). Senescent endothelial cells secrete factors that attract peripheral leukocytes (81, 82) and upregulate surface leukocyte adhesion molecules (24, 73, 83, 84). Moreover, transcriptomic (85) and proteomic (86) analyses of human AD cortical microvessels have demonstrated upregulation of pathways related to senescence, including inflammation, leukocyte migration, and cell adhesion. However, this is the first study to directly examine the relationship between AD-related tau pathology and vascular senescence in both the cerebral vasculature and peripheral blood.
We report robust upregulation of the senescence- and leukocyte adhesion-associated genes CXCL8 (IL8), SERPINE1 (PAI-1), CXCL1, CXCL2, ICAM-2, and TIE1 in B3 cortical microvasculature. Interestingly, these genes are also involved in DNA damage response (DDR) signaling, a pathway that potently mediates cellular senescence (87–89). The pronounced upregulation of CXCL8 and SERPINE1 in B3 microvessels is notable, given their reported upregulation in senescent endothelial cells (23, 24, 90–92) and synergistic roles in leukocyte-endothelial adhesion (93). ICAM-2 also mediates leukocyte-endothelial adhesion in cerebral vasculature (94), and neurons upregulate ICAM-2 with the formation of tau NFTs (95). Tau burden and cerebrovascular pathology also appeared to additively interact with regard to microvascular expression of MAdCAM-1, a leukocyte adhesion molecule that is upregulated in the CNS following chronic vascular inflammation (96–98). Since MAdCAM-1 detection was variable across samples, future studies in a larger cohort are warranted to clarify the relationship between MAdCAM-1, tau pathology, and cerebrovascular pathologies.
One marker of cerebrovascular dysfunction is impaired blood-brain barrier (BBB) integrity, a phenomenon that has been linked to tau misfolding and aggregation (16, 99, 100). However, tau-independent BBB dysfunction is also observed with cognitive decline (101), possibly in relation to APOE ε4-associated BBB impairment in the hippocampus and medial temporal lobes (102). The BBB is comprised largely of endothelial cell tight junctions, supported by tight junction proteins such as OCLN, CLDN5, CDH5, and TJP1. We report here no change in gene expression of these four tight junction markers with tau progression in the AD cortical vasculature. This is a departure from previous reports of reduced cerebrovascular tight junction protein expression in AD (103) as well as the Parkinsonism–dementia complex of Guam (Guam PDC) tauopathy (104). However, Yamazaki et al. (17) showed that senescent endothelial cells exhibit aberrant tight junction protein localization without differences in overall tight junction protein expression, possibly reconciling this discrepancy. Furthermore, the assembly and localization of such tight junction proteins is largely influenced by post-transcriptional modifications (105).
Antemortem plasma protein expression did not associate with postmortem brain microvascular gene expression in this small sample, with the exception of NOS3. NOS3 is the endothelial source of nitric oxide (NO), a potent vasodilator that mediates vascular homeostasis and cerebral blood flow (106). While the association between gene and protein expression is generally tenuous, even within the same CNS cell type (107), future analysis in a larger cohort could determine the utility of NOS3 as a peripheral readout of cerebrovascular health. Of note, the temporal stability of a peripheral biomarker is an important consideration for its utility in studying the pathogenesis and prediction of AD. This is particularly relevant for PAI-1, which greatly fluctuated over time in both B1 and B3 plasma, possibly owing to its reported circadian variations (108) and cell surface binding after secretion (109). Furthermore, observed discrepancies in senescence-associated secretory phenotype (SASP) biomarker expression in cerebral vasculature vs. peripheral blood may be attributable to mismatches in transcriptional and translational upregulation.
One advantage of this study includes the preservation of the intercellular milieu in the neurovascular unit by isolating intact cortical microvessels. Additionally, the use of tau pathology scores, microvascular gene expression, and plasma protein levels from the same subjects enabled direct comparison across modalities. However, this study is limited by its small sample size, particularly in the plasma protein subgroup, with a bias toward AD samples. We included potential confounders such as sex, age, and cerebrovascular pathology in statistical models, though their roles in tau pathology and vascular dysfunction cannot be excluded. Additionally, AD commonly presents with other co-morbidities, including TDP-43 proteinopathy and Lewy body disease (110–112); future studies are warranted to compare tau-related vascular changes with and without such co-morbidities. Due to insufficient statistical power and a high degree of multicollinearity in this dataset, we did not adjust for APOE genotype nor for amyloid-beta burden, the effects of which also cannot be ruled out here. However, we note that the transcriptional changes shown in human AD cortical vasculature here are similar to those observed in the Tg4510 mouse model (15, 27), which exclusively over-expresses pathological tau. Future studies investigating individuals with high NFT burden but low amyloid accumulation, as seen in Primary Age Related Tauopathy (PART) (113), will be needed to confirm the specificity of these changes.
Looking forward, a larger sample size with evenly distributed pre-mortem intervals may reveal changes in the systemic vasculature that associate with cortical microvessel gene expression and tau pathology. Alternatively, rapidly advancing tau-PET neuroimaging represents a viable technique to track peripheral biomarker changes with tau accumulation in the AD brain. Tau tracer uptake correlates with the Braak stages of tau pathology progression in AD (114–116), and Ashton et al. (117) recently showed that tau-PET can be used to compare antemortem plasma protein expression and antemortem tau pathology. Real-time investigations of senescence-associated biomarker secretion in plasma and cerebral tau accumulation could better elucidate their relationship to vascular dysfunction in AD.
In conclusion, microvessels isolated from the human AD prefrontal cortex with extensive tau pathology upregulate genes involved in endothelial senescence and in recruiting leukocytes to the endothelium, which may contribute to AD-related vascular dysfunction and impaired cerebral blood flow. Future studies could identify peripheral biomarkers that are associated with vascular senescence and its relation to tau pathology in the human AD brain.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by Institutional Review Board of Partners Healthcare, Massachusetts General Hospital, Boston, MA. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
AB, SA, BH, and RB contributed to study concepts and experimental design. AB, MH, BC, and RB collected data. AB performed statistical analyses, with interpretation guided by input from BC, SD, and RB. MF provided access to biomaterials in the Massachusetts Alzheimer's Disease Research Center. AB drafted the manuscript, with critical revisions from BC, SD, BH, and RB for intellectual content. All co-authors reviewed the final manuscript and approved it for submission.
Funding
This work was supported by the NIH/NIA K99 AG061259 and by a New Vision Investigator Award. We would also like to thank the donors and their families for their contribution to the Massachusetts ADRC (NIH/NIA P30AG062421).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. Pia K. Webb for coordination of the longitudinal plasma sample collection and analysis. We thank Dr. Lori Chibnik for guidance and review of statistical analyses employed in this study. We also thank Patrick M. Dooley and Theresa R. Connors for identifying and procuring postmortem tissue used in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.575953/full#supplementary-material
References
1. de la Torre JC, Mussivan T. Can disturbed brain microcirculation cause Alzheimer's disease? Neurol Res. (1993) 15:146–53. doi: 10.1080/01616412.1993.11740127
2. de la Torre JC. The vascular hypothesis of Alzheimer's disease: bench to bedside and beyond. Neurodegener Dis. (2010) 7:116–21. doi: 10.1159/000285520
3. Roher AE, Debbins JP, Malek-Ahmadi M, Chen K, Pipe JG, Maze S, et al. Cerebral blood flow in Alzheimer's disease. Vasc Health Risk Manag. (2012) 8:599–611. doi: 10.2147/VHRM.S34874
4. Yarchoan M, Xie SX, Kling MA, Toledo JB, Wolk DA, Lee EB, et al. Cerebrovascular atherosclerosis correlates with Alzheimer pathology in neurodegenerative dementias. Brain. (2012) 135:3749–56. doi: 10.1093/brain/aws271
5. Toledo JB, Arnold SE, Raible K, Brettschneider J, Xie SX, Grossman M, et al. Contribution of cerebrovascular disease in autopsy confirmed neurodegenerative disease cases in the National Alzheimer's Coordinating Centre. Brain. (2013) 136:2697–706. doi: 10.1093/brain/awt188
6. Nagahama Y, Nabatame H, Okina T, Yamauchi H, Narita M, Fujimoto N, et al. Cerebral correlates of the progression rate of the cognitive decline in probable Alzheimer's disease. Eur Neurol. (2003) 50:1–9. doi: 10.1159/000070851
7. Benedictus MR, Leeuwis AE, Binnewijzend MAA, Kuijer JPA, Scheltens P, Barkhof F, et al. Lower cerebral blood flow is associated with faster cognitive decline in Alzheimer's disease. Eur Radiol. (2017) 27:1169–75. doi: 10.1007/s00330-016-4450-z
8. Bradley KM, O'Sullivan VT, Soper NDW, Nagy Z, King EM-F, Smith D, et al. Cerebral perfusion SPET correlated with Braak pathological stage in Alzheimer's disease. Brain. (2002) 125:1772–81. doi: 10.1093/brain/awf185
9. Binnewijzend MAA, Benedictus MR, Kuijer JPA, van der Flier WM, Teunissen CE, Prins ND, et al. Cerebral perfusion in the predementia stages of Alzheimer's disease. Eur Radiol. (2016) 26:506–14. doi: 10.1007/s00330-015-3834-9
10. Nedelska Z, Senjem ML, Przybelski SA, Lesnick TG, Lowe VJ, Boeve BF, et al. Regional cortical perfusion on arterial spin labeling MRI in dementia with Lewy bodies: associations with clinical severity, glucose metabolism and tau PET. NeuroImage Clin. (2018) 19:939–47. doi: 10.1016/j.nicl.2018.06.020
11. Matsumoto M, Hatakeyama T, Morimoto K, Yanagihara T. Cerebral blood flow and neuronal damage during progressive cerebral ischemia in gerbils. Stroke. (1990) 21:1470–7. doi: 10.1161/01.STR.21.10.1470
12. Sagare AP, Bell RD, Zlokovic BV. Neurovascular defects and faulty amyloid-β vascular clearance in Alzheimer's disease. J Alzheimer's Dis. (2012) 33:S87–S100. doi: 10.3233/JAD-2012-129037
13. Coppé J-P, Desprez P-Y, Krtolica A, Campisi J. The senescence-associated secretory phenotype: the dark side of tumor suppression. Annu Rev Pathol Mech Dis. (2010) 5:99–118. doi: 10.1146/annurev-pathol-121808-102144
14. Bussian TJ, Aziz A, Meyer CF, Swenson BL, van Deursen JM, Baker DJ. Clearance of senescent glial cells prevents tau-dependent pathology and cognitive decline. Nature. (2018) 562:578–82. doi: 10.1038/s41586-018-0543-y
15. Musi N, Valentine JM, Sickora KR, Baeuerle E, Thompson CS, Shen Q, et al. Tau protein aggregation is associated with cellular senescence in the brain. Aging Cell. (2018) 17:e12840. doi: 10.1111/acel.12840
16. Blair LJ, Frauen HD, Zhang B, Nordhues BA, Bijan S, Lin Y-C, et al. Tau depletion prevents progressive blood-brain barrier damage in a mouse model of tauopathy. Acta Neuropathol Commun. (2015) 3:8. doi: 10.1186/s40478-015-0186-2
17. Yamazaki Y, Baker DJ, Tachibana M, Liu C-C, van Deursen JM, Brott TG, et al. Vascular cell senescence contributes to blood–brain barrier breakdown. Stroke. (2016) 47:1068–77. doi: 10.1161/STROKEAHA.115.010835
18. Tian X-L, Li Y. Endothelial cell senescence and age-related vascular diseases. J Genet Genomics. (2014) 41:485–95. doi: 10.1016/j.jgg.2014.08.001
19. Merlini M, Wanner D, Nitsch RM. Tau pathology-dependent remodelling of cerebral arteries precedes Alzheimer's disease-related microvascular cerebral amyloid angiopathy. Acta Neuropathol. (2016) 131:737–52. doi: 10.1007/s00401-016-1560-2
20. Holmes HE, Colgan N, Ismail O, Ma D, Powell NM, O'Callaghan JM, et al. Imaging the accumulation and suppression of tau pathology using multiparametric MRI. Neurobiol Aging. (2016) 39:184–94. doi: 10.1016/j.neurobiolaging.2015.12.001
21. Qiu L, Ng G, Tan EK, Liao P, Kandiah N, Zeng L. Chronic cerebral hypoperfusion enhances Tau hyperphosphorylation and reduces autophagy in Alzheimer's disease mice. Sci Rep. (2016) 6:23964. doi: 10.1038/srep23964
22. Graves SI, Baker DJ. Implicating endothelial cell senescence to dysfunction in the ageing and diseased brain. Basic Clin Pharmacol Toxicol. (2020) 127:102–10. doi: 10.1111/bcpt.13403
23. Khan SY, Awad EM, Oszwald A, Mayr M, Yin X, Waltenberger B, et al. Premature senescence of endothelial cells upon chronic exposure to TNFα can be prevented by N-acetyl cysteine and plumericin. Sci Rep. (2017) 7:39501. doi: 10.1038/srep39501
24. Prattichizzo F, De Nigris V, Mancuso E, Spiga R, Giuliani A, Matacchione G, et al. Short-term sustained hyperglycaemia fosters an archetypal senescence-associated secretory phenotype in endothelial cells and macrophages. Redox Biol. (2018) 15:170–81. doi: 10.1016/j.redox.2017.12.001
25. Kuilman T, Michaloglou C, Vredeveld LCW, Douma S, van Doorn R, Desmet CJ, et al. Oncogene-induced senescence relayed by an interleukin-dependent inflammatory network. Cell. (2008) 133:1019–31. doi: 10.1016/j.cell.2008.03.039
26. Hubackova S, Krejcikova K, Bartek J, Hodny Z. IL1- and TGFβ-Nox4 signaling, oxidative stress and DNA damage response are shared features of replicative, oncogene-induced, and drug-induced paracrine ‘Bystander senescence.' Aging. (2012) 4:932–51. doi: 10.18632/aging.100520
27. Bennett RE, Robbins AB, Hu M, Cao X, Betensky RA, Clark T, et al. Tau induces blood vessel abnormalities and angiogenesis-related gene expression in P301L transgenic mice and human Alzheimer's disease. Proc Natl Acad Sci USA. (2018) 115:E1289–98. doi: 10.1073/pnas.1710329115
28. Hyman BT, Phelps CH, Beach TG, Bigio EH, Cairns NJ, Carrillo MC, et al. National Institute on Aging-Alzheimer's Association guidelines for the neuropathologic assessment of Alzheimer's disease. Alzheimers Dement. (2012) 8:1–13. doi: 10.1016/j.jalz.2011.10.007
29. Boulay AC, Saubaméa B, Declèves X, Cohen-Salmon M. Purification of mouse brain vessels. J Vis Exp. (2015) e53208. doi: 10.3791/53208
30. Hartz AMS, Schulz JA, Sokola BS, Edelmann SE, Shen AN, Rempe RG, et al. Isolation of cerebral capillaries from fresh human brain tissue. J Vis Exp. (2018) 57346. doi: 10.3791/57346
31. Kiss T, Nyúl-Tóth Á, Balasubramanian P, Tarantini S, Ahire C, DelFavero J, et al. Single-cell RNA sequencing identifies senescent cerebromicrovascular endothelial cells in the aged mouse brain. GeroScience. (2020) 42:429–44. doi: 10.1007/s11357-020-00177-1
32. Cotran R, Mayadas-Norton T. Endothelial adhesion molecules in health and disease. Pathol Biol. (1998) 46:164–70.
33. Reglero-Real N, Colom B, Bodkin JV, Nourshargh S. Endothelial cell junctional adhesion molecules. Arteriosc Thromb Vasc Biol. (2016) 36:2048–57. doi: 10.1161/ATVBAHA.116.307610
34. Zahr A, Alcaide P, Yang J, Jones A, Gregory M, dela Paz NG, et al. Endomucin prevents leukocyte–endothelial cell adhesion and has a critical role under resting and inflammatory conditions. Nat Commun. (2016) 7:10363. doi: 10.1038/ncomms10363
35. Piovella F, Nalli G, Malamani GD, Majolino I, Frassoni F, Sitar GM, et al. The ultrastructural localization of factor VIII-antigen in human platelets, megakaryocytes and endothelial cells utilizing a ferritin-labelled antibody. Br J Haematol. (1978) 39:209–13. doi: 10.1111/j.1365-2141.1978.tb01090.x
36. Nakhaei-Nejad M, Farhan M, Mojiri A, Jabbari H, Murray AG, Jahroudi N. Regulation of von willebrand factor gene in endothelial cells that are programmed to pluripotency and differentiated back to endothelial cells. Stem Cells. (2019) 37:542–54. doi: 10.1002/stem.2978
37. Partanen J, Armstrong E, Mäkelä TP, Korhonen J, Sandberg M, Renkonen R, et al. A novel endothelial cell surface receptor tyrosine kinase with extracellular epidermal growth factor homology domains. Mol Cell Biol. (1992) 12:1698–707. doi: 10.1128/MCB.12.4.1698
38. Dumont DJ, Yamaguchi TP, Conion RA, Rossant J, Breitman ML. Tek, a novel tyrosine kinase gene located on mouse chromosome 4, is expressed in endothelial cells and their presumptive precursors. Oncogene. (1992) 7:1471.
39. Maher F, Vannucci SJ, Simpson IA. Glucose transporter isoforms in brain: absence of GLUT3 from the blood—brain barrier. J Cereb Blood Flow Metab. (1993) 13:342–5. doi: 10.1038/jcbfm.1993.43
40. Zimmermann K, Opitz N, Dedio J, Renne C, Muller-Esterl W, Oess S. NOSTRIN: a protein modulating nitric oxide release and subcellular distribution of endothelial nitric oxide synthase. Proc Natl Acad Sci USA. (2002) 99:17167–72. doi: 10.1073/pnas.252345399
41. DeLisser HM, Newman PJ, Albelda SM. Molecular and functional aspects of PECAM-1/CD31. Immunol Today. (1994) 15:490–5. doi: 10.1016/0167-5699(94)90195-3
42. Guo X, Geng M, Du G. Glucose transporter 1, distribution in the brain and in neural disorders: its relationship with transport of neuroactive drugs through the blood-brain barrier. Biochem Genet. (2005) 43:175–87. doi: 10.1007/s10528-005-1510-5
43. Cuevas P, Gutierrez-Diaz JA, Reimers D, Dujovny M, Diaz FG, Ausman JI. Pericyte endothelial gap junctions in human cerebral capillaries. Anat Embryol. (1984) 170:155–9. doi: 10.1007/BF00319000
44. Wallez Y, Huber P. Endothelial adherens and tight junctions in vascular homeostasis, inflammation and angiogenesis. Biochim Biophys Acta Biomembr. (2008) 1778:794–809. doi: 10.1016/j.bbamem.2007.09.003
45. Dejana E, Orsenigo F, Molendini C, Baluk P, McDonald DM. Organization and signaling of endothelial cell-to-cell junctions in various regions of the blood and lymphatic vascular trees. Cell Tissue Res. (2009) 335:17–25. doi: 10.1007/s00441-008-0694-5
46. Thurston G, Kitajewski J. VEGF and Delta-Notch: interacting signalling pathways in tumour angiogenesis. Br J Cancer. (2008) 99:1204–9. doi: 10.1038/sj.bjc.6604484
47. Favara DM, Zois CE, Haider S, Pires E, Sheldon H, McCullagh J, et al. ADGRL4/ELTD1 silencing in endothelial cells induces ACLY and SLC25A1 and alters the cellular metabolic profile. Metabolites. (2019) 9:287. doi: 10.3390/metabo9120287
48. Tarkowski E, Issa R, Sjögren M, Wallin A, Blennow K, Tarkowski A, et al. Increased intrathecal levels of the angiogenic factors VEGF and TGF-beta in Alzheimer's disease and vascular dementia. Neurobiol Aging. (2002) 23:237–43. doi: 10.1016/S0197-4580(01)00285-8
49. Cho S-J, Park MH, Han C, Yoon K, Koh YH. VEGFR2 alteration in Alzheimer's disease. Sci Rep. (2017) 7:17713. doi: 10.1038/s41598-017-18042-1
50. Cho S-J, Yun S-M, Jo C, Jeong J, Park MH, Han C, et al. Altered expression of Notch1 in Alzheimer's disease. PLoS ONE. (2019) 14:e0224941. doi: 10.1371/journal.pone.0224941
51. Poss KD, Tonegawa S. Reduced stress defense in heme oxygenase 1-deficient cells. Proc Natl Acad Sci USA. (1997) 94:10925–30. doi: 10.1073/pnas.94.20.10925
53. Herz J, Strickland DK. LRP: a multifunctional scavenger and signaling receptor. J Clin Invest. (2001) 108:779–84. doi: 10.1172/JCI200113992
54. Lillis AP, Van Duyn LB, Murphy-Ullrich JE, Strickland DK. LDL receptor-related protein 1: unique tissue-specific functions revealed by selective gene knockout studies. Physiol Rev. (2008) 88:887–918. doi: 10.1152/physrev.00033.2007
55. Kanekiyo T, Bu G. The low-density lipoprotein receptor-related protein 1 and amyloid-Î2 clearance in Alzheimer's disease. Front Aging Neurosci. (2014) 6:93. doi: 10.3389/fnagi.2014.00093
56. Cahoy JD, Emery B, Kaushal A, Foo LC, Zamanian JL, Christopherson KS, et al. A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J Neurosci. (2008) 28:264–78. doi: 10.1523/JNEUROSCI.4178-07.2008
57. Yuan S-M. α-smooth muscle actin and ACTA2 gene expressions in vasculopathies. Rev Bras Cir Cardiovasc. (2015) 30:644–9. doi: 10.5935/1678-9741.20150081
58. Grabert K, Michoel T, Karavolos MH, Clohisey S, Baillie JK, Stevens MP, et al. Microglial brain region–dependent diversity and selective regional sensitivities to aging. Nat Neurosci. (2016) 19:504–16. doi: 10.1038/nn.4222
59. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.126
60. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. (1995) 57:289–300.
61. Kulesa A, Krzywinski M, Blainey P, Altman N. Sampling distributions and the bootstrap: the bootstrap can be used to assess uncertainty of sample estimates. Nat Methods. (2015) 12:477. doi: 10.1038/nmeth.3414
62. Carlyle BC, Kitchen RR, Kanyo JE, Voss EZ, Pletikos M, Sousa AMM, et al. A multiregional proteomic survey of the postnatal human brain. Nat Neurosci. (2017) 20:1787–95. doi: 10.1038/s41593-017-0011-2
63. Braak H, Del Tredici K. The preclinical phase of the pathological process underlying sporadic Alzheimer's disease. Brain. (2015) 138:2814–33. doi: 10.1093/brain/awv236
64. Bussire T, Gold G, Kvari E, Giannakopoulos P, Bouras C, Perl D, et al. Stereologic analysis of neurofibrillary tangle formation in prefrontal cortex area 9 in aging and Alzheimer's disease. Neuroscience. (2003) 117:577–92. doi: 10.1016/S0306-4522(02)00942-9
65. Vemuri P, Lesnick TG, Przybelski SA, Knopman DS, Lowe VJ, Graff-Radford J, et al. Age, vascular health, and Alzheimer disease biomarkers in an elderly sample. Ann Neurol. (2017) 82:706–18. doi: 10.1002/ana.25071
66. Fernando MS, Ince PG. Vascular pathologies and cognition in a population-based cohort of elderly people. J Neurol Sci. (2004) 226:13–7. doi: 10.1016/j.jns.2004.09.004
67. Schneider JA, Arvanitakis Z, Bang W, Bennett DA. Mixed brain pathologies account for most dementia cases in community-dwelling older persons. Neurology. (2007) 69:2197–204. doi: 10.1212/01.wnl.0000271090.28148.24
68. Williams S, Chalmers K, Wilcock GK, Love S. Relationship of neurofibrillary pathology to cerebral amyloid angiopathy in Alzheimer's disease. Neuropathol Appl Neurobiol. (2005) 31:414–21. doi: 10.1111/j.1365-2990.2005.00663.x
69. Castillo-Carranza DL, Nilson AN, Van Skike CE, Jahrling JB, Patel K, Garach P, et al. Cerebral microvascular accumulation of Tau oligomers in Alzheimer's disease and related tauopathies. Aging Dis. (2017) 8:257–66. doi: 10.14336/AD.2017.0112
70. Kim HJ, Park S, Cho H, Jang YK, San Lee J, Jang H, et al. Assessment of extent and role of Tau in subcortical vascular cognitive impairment using 18 F-AV1451 positron emission tomography imaging. JAMA Neurol. (2018) 75:999. doi: 10.1001/jamaneurol.2018.0975
71. Camós S, Mallolas J. Experimental models for assaying microvascular endothelial cell pathophysiology in stroke. Molecules. (2010) 15:9104–34. doi: 10.3390/molecules15129104
72. Wilhelmsen K, Mesa KR, Prakash A, Xu F, Hellman J. Activation of endothelial TLR2 by bacterial lipoprotein upregulates proteins specific for the neutrophil response. Innate Immun. (2012) 18:602–16. doi: 10.1177/1753425911429336
73. Coleman PR, Chang G, Hutas G, Grimshaw M, Vadas MA, Gamble JR. Age-associated stresses induce an anti-inflammatory senescent phenotype in endothelial cells. Aging. (2013) 5:913–24. doi: 10.18632/aging.100622
74. Zenaro E, Piacentino G, Constantin G. The blood-brain barrier in Alzheimer's disease. Neurobiol Dis. (2017) 107:41–56. doi: 10.1016/j.nbd.2016.07.007
75. Cruz Hernández JC, Bracko O, Kersbergen CJ, Muse V, Haft-Javaherian M, Berg M, et al. Neutrophil adhesion in brain capillaries reduces cortical blood flow and impairs memory function in Alzheimer's disease mouse models. Nat Neurosci. (2019) 22:413–20. doi: 10.1038/s41593-018-0329-4
76. Yata K, Nishimura Y, Unekawa M, Tomita Y, Suzuki N, Tanaka T, et al. In vivo imaging of the mouse neurovascular unit under chronic cerebral hypoperfusion. Stroke. (2014) 45:3698–703. doi: 10.1161/STROKEAHA.114.005891
77. Bullock R, Maxwell WL, Graham DI, Teasdale GM, Adams JH. Glial swelling following human cerebral contusion: an ultrastructural study. J Neurol Neurosurg Psychiatry. (1991) 54:427–34. doi: 10.1136/jnnp.54.5.427
78. Mazzoni MC, Schmid-Schonbein GW. Mechanisms and consequences of cell activation in the microcirculation. Cardiovasc Res. (1996) 32:709–19. doi: 10.1016/S0008-6363(96)00146-0
79. Jespersen SN, Ostergaard L. The roles of cerebral blood flow, capillary transit time heterogeneity, and oxygen tension in brain oxygenation and metabolism. J Cereb Blood Flow Metab. (2012) 32:264–77. doi: 10.1038/jcbfm.2011.153
80. Østergaard L, Aamand R, Gutiérrez-Jiménez E, Ho Y-CL, Blicher JU, Madsen SM, et al. The capillary dysfunction hypothesis of Alzheimer's disease. Neurobiol Aging. (2013) 34:1018–31. doi: 10.1016/j.neurobiolaging.2012.09.011
81. Davalos AR, Coppe J-P, Campisi J, Desprez P-Y. Senescent cells as a source of inflammatory factors for tumor progression. Cancer Metastasis Rev. (2010) 29:273–83. doi: 10.1007/s10555-010-9220-9
82. Lasry A, Ben-Neriah Y. Senescence-associated inflammatory responses: aging and cancer perspectives. Trends Immunol. (2015) 36:217–28. doi: 10.1016/j.it.2015.02.009
83. Gorgoulis VG, Pratsinis H, Zacharatos P, Demoliou C, Sigala F, Asimacopoulos PJ, et al. p53-Dependent ICAM-1 overexpression in senescent human cells identified in atherosclerotic lesions. Lab Investig. (2005) 85:502–11. doi: 10.1038/labinvest.3700241
84. Pietronigro E, Zenaro E, Constantin G. Imaging of leukocyte trafficking in Alzheimer's disease. Front Immunol. (2016) 7:33. doi: 10.3389/fimmu.2016.00033
85. Wang S, Qaisar U, Yin X, Grammas P. Gene expression profiling in Alzheimer's disease brain microvessels. J Alzheimer Dis. (2012) 31:193–205. doi: 10.3233/JAD-2012-120454
86. Grammas P, Ovase R. Inflammatory factors are elevated in brain microvessels in Alzheimer's disease. Neurobiol Aging. (2001) 22:837–42. doi: 10.1016/S0197-4580(01)00276-7
87. Di Leonardo A, Linke SP, Clarkin K, Wahl GM. DNA damage triggers a prolonged p53-dependent G1 arrest and long-term induction of Cip1 in normal human fibroblasts. Genes Dev. (1994) 8:2540–51. doi: 10.1101/gad.8.21.2540
88. Nakamura AJ, Chiang YJ, Hathcock KS, Horikawa I, Sedelnikova OA, Hodes RJ, et al. Both telomeric and non-telomeric DNA damage are determinants of mammalian cellular senescence. Epigenetics Chromatin. (2008) 1:6. doi: 10.1186/1756-8935-1-6
89. Bielak-Zmijewska A, Mosieniak G, Sikora E. Is DNA damage indispensable for stress-induced senescence? Mech Ageing Dev. (2018) 170:13–21. doi: 10.1016/j.mad.2017.08.004
90. Hampel B, Fortschegger K, Ressler S, Chang M, Unterluggauer H, Breitwieser A, et al. Increased expression of extracellular proteins as a hallmark of human endothelial cell in vitro senescence. Exp Gerontol. (2006) 41:474–81. doi: 10.1016/j.exger.2006.03.001
91. Kortlever RM, Higgins PJ, Bernards R. Plasminogen activator inhibitor-1 is a critical downstream target of p53 in the induction of replicative senescence. Nat Cell Biol. (2006) 8:877–84. doi: 10.1038/ncb1448
92. Vaughan DE, Rai R, Khan SS, Eren M, Ghosh AK. Plasminogen activator inhibitor-1 is a marker and a mediator of senescence. Arterioscler Thromb Vasc Biol. (2017) 37:1446–52. doi: 10.1161/ATVBAHA.117.309451
93. Marshall LJ, Ramdin LSP, Brooks T, DPhil PC, Shute JK. Plasminogen activator inhibitor-1 supports IL-8-mediated neutrophil transendothelial migration by inhibition of the constitutive shedding of endothelial IL-8/heparan sulfate/syndecan-1 complexes. J Immunol. (2003) 171:2057–65. doi: 10.4049/jimmunol.171.4.2057
94. Halai K, Whiteford J, Ma B, Nourshargh S, Woodfin A. ICAM-2 facilitates luminal interactions between neutrophils and endothelial cells in vivo. J Cell Sci. (2014) 127:620–9. doi: 10.1242/jcs.137463
95. Dunckley T, Beach TG, Ramsey KE, Grover A, Mastroeni D, Walker DG, et al. Gene expression correlates of neurofibrillary tangles in Alzheimer's disease. Neurobiol Aging. (2006) 27:1359–71. doi: 10.1016/j.neurobiolaging.2005.08.013
96. Hänninen A, Taylor C, Streeter PR, Stark LS, Sarte JM, Shizuru JA, et al. Vascular addressins are induced on islet vessels during insulitis in nonobese diabetic mice and are involved in lymphoid cell binding to islet endothelium. J Clin Invest. (1993) 92:2509–15. doi: 10.1172/JCI116859
97. Farkas S, Hornung M, Sattler C, Edtinger K, Steinbauer M, Anthuber M, et al. Blocking MAdCAM-1 in vivo reduces leukocyte extravasation and reverses chronic inflammation in experimental colitis. Int J Colorectal Dis. (2006) 21:71–8. doi: 10.1007/s00384-004-0709-y
98. Kuhbandner K, Hammer A, Haase S, Terbrack E, Hoffmann A, Schippers A, et al. MAdCAM-1-mediated intestinal lymphocyte homing is critical for the development of active experimental autoimmune encephalomyelitis. Front Immunol. (2019) 10:903. doi: 10.3389/fimmu.2019.00903
99. Kovac A, Zilkova M, Deli MA, Zilka N, Novak M. Human truncated tau is using a different mechanism from amyloid-β to damage the blood-brain barrier. J Alzheimer Dis. (2009) 18:897–906. doi: 10.3233/JAD-2009-1197
100. Majerova P, Michalicova A, Cente M, Hanes J, Vegh J, Kittel A, et al. Trafficking of immune cells across the blood-brain barrier is modulated by neurofibrillary pathology in tauopathies. PLoS ONE. (2019) 14:e0217216. doi: 10.1371/journal.pone.0217216
101. Nation DA, Sweeney MD, Montagne A, Sagare AP, D'Orazio LM, Pachicano M, et al. Blood–brain barrier breakdown is an early biomarker of human cognitive dysfunction. Nat Med. (2019) 25:270–6. doi: 10.1038/s41591-018-0297-y
102. Montagne A, Nation DA, Sagare AP, Barisano G, Sweeney MD, Chakhoyan A, et al. APOE4 leads to blood–brain barrier dysfunction predicting cognitive decline. Nature. (2020) 581:71–6. doi: 10.1038/s41586-020-2247-3
103. Yamazaki Y, Shinohara M, Shinohara M, Yamazaki A, Murray ME, Liesinger AM, et al. Selective loss of cortical endothelial tight junction proteins during Alzheimer's disease progression. Brain. (2019) 142:1077–92. doi: 10.1093/brain/awz011
104. Majerova P, Garruto RM, Kovac A. Cerebrovascular inflammation is associated with tau pathology in Guam parkinsonism dementia. J Neural Trans. (2018) 125:1013–25. doi: 10.1007/s00702-018-1883-3
105. Huang SF, Othman A, Koshkin A, Fischer S, Fischer D, Zamboni N, et al. Astrocyte glutathione maintains endothelial barrier stability. Redox Biol. (2020). 34:101576. doi: 10.1016/j.redox.2020.101576
106. Yang Y-M, Huang A, Kaley G, Sun D. eNOS uncoupling and endothelial dysfunction in aged vessels. Am J Physiol Circ Physiol. (2009) 297:H1829–36. doi: 10.1152/ajpheart.00230.2009
107. Sharma K, Schmitt S, Bergner CG, Tyanova S, Kannaiyan N, Manrique-Hoyos N, et al. Cell type– and brain region–resolved mouse brain proteome. Nat Neurosci. (2015) 18:1819–31. doi: 10.1038/nn.4160
108. Vaughan DE. PAI-1 and atherothrombosis. J Thromb Haemost. (2005) 3:1879–83. doi: 10.1111/j.1538-7836.2005.01420.x
109. Gerenu G, Martisova E, Ferrero H, Carracedo M, Rantamäki T, Ramirez MJ, et al. Modulation of BDNF cleavage by plasminogen-activator inhibitor-1 contributes to Alzheimer's neuropathology and cognitive deficits. Biochim Biophys Acta. (2017) 1863:991–1001. doi: 10.1016/j.bbadis.2017.01.023
110. Rahimi J, Kovacs GG. Prevalence of mixed pathologies in the aging brain. Alzheimer Res Ther. (2014) 6:82. doi: 10.1186/s13195-014-0082-1
111. James BD, Wilson RS, Boyle PA, Trojanowski JQ, Bennett DA, Schneider JA. TDP-43 stage, mixed pathologies, and clinical Alzheimer's-type dementia. Brain. (2016) 139:2983–93. doi: 10.1093/brain/aww224
112. Robinson JL, Lee EB, Xie SX, Rennert L, Suh E, Bredenberg C, et al. Neurodegenerative disease concomitant proteinopathies are prevalent, age-related and APOE4-associated. Brain. (2018) 141:2181–93. doi: 10.1093/brain/awy146
113. Crary JF, Trojanowski JQ, Schneider JA, Abisambra JF, Abner EL, Alafuzoff I, et al. Primary age-related tauopathy (PART): a common pathology associated with human aging. Acta Neuropathol. (2014) 128:755–66. doi: 10.1007/s00401-014-1349-0
114. Mielke MM, Hagen CE, Xu J, Chai X, Vemuri P, Lowe VJ, et al. Plasma phospho-tau181 increases with Alzheimer's disease clinical severity and is associated with tau- and amyloid-positron emission tomography. Alzheimer Dement. (2018) 14:989–97. doi: 10.1016/j.jalz.2018.02.013
115. Janelidze S, Mattsson N, Palmqvist S, Smith R, Beach TG, Serrano GE, et al. Plasma P-tau181 in Alzheimer's disease: relationship to other biomarkers, differential diagnosis, neuropathology and longitudinal progression to Alzheimer's dementia. Nat Med. (2020) 26:379–86. doi: 10.1038/s41591-020-0755-1
116. Thijssen EH, La Joie R, Wolf A, Strom A, Wang P, Iaccarino L, et al. Diagnostic value of plasma phosphorylated tau181 in Alzheimer's disease and frontotemporal lobar degeneration. Nat Med. (2020) 26:387–97. doi: 10.1038/s41591-020-0762-2
Keywords: Alzheimer's disease, tau pathology, neurofibrillary tangles, vascular dysfunction, endothelial senescence, gene expression, plasma biomarkers
Citation: Bryant AG, Hu M, Carlyle BC, Arnold SE, Frosch MP, Das S, Hyman BT and Bennett RE (2020) Cerebrovascular Senescence Is Associated With Tau Pathology in Alzheimer's Disease. Front. Neurol. 11:575953. doi: 10.3389/fneur.2020.575953
Received: 24 June 2020; Accepted: 12 August 2020;
Published: 16 September 2020.
Edited by:
Sonia Do Carmo, McGill University, CanadaReviewed by:
Andrej Kovac, Slovak Academy of Sciences (SAS), SlovakiaAntonio Giuliano Zippo, National Research Council, Italy
Copyright © 2020 Bryant, Hu, Carlyle, Arnold, Frosch, Das, Hyman and Bennett. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rachel E. Bennett, rebennett@mgh.harvard.edu