To analyze the association of α-syn load and distribution with Aβ and tau pathology in AD, we analyzed immunohistochemical stainings of up to 28 brain regions per case in 72 AD patients (Table 1). The cohort had a mean age at death of 72.8 years (± 11.5 years standard deviation). 56% of the subjects were female. Most of the cases had a Braak and Braak stage VI and a Thal phase 5, corresponding to an advanced stage of AD. For 66 cases, information about the ApoE status was available, revealing at least one ApoE4 allele in 58% of the subjects.
The deposit covered area was automatically quantified by random forest classifiers in 1016 regions in α-syn stainings, 1292 regions in tau stainings, and 1098 regions in Aβ stainings. By thresholding, AD patients were assigned to αSyn-, comprising 41% of the cases, and αSyn+, including 59% of the cases (Figs. 1 and 2). The latter were further divided in three α-syn distribution patterns (Fig. 1, Fig. 2, Table 2, Supplementary Table 6): αSyn + A, comprising around one third of the α-syn positive cases with an almost exclusive amygdala-entorhinal α-syn load; αSyn + B, including around 12% of the α-syn positive cases and characterized by a high brainstem α-syn load without cortical spread and a low amygdala involvement; αSyn + C, comprising around half of the α-syn positive cases and presenting with at least focal cortical α-syn deposits together with the highest amygdala-entorhinal and a relatively high brainstem α-syn load. All groups and subgroups were evaluated regarding their Aβ and tau load, revealing distinct loads in different brain regions.
Table 1
Demographic and neuropathological overview of α-syn groups in Alzheimer's disease
| | Available n | All | αSyn- | αSyn+ | Statistic, p-value |
|---|
n (%) | 72 | 72 (100%) | 29 (40%) | 43 (60%) | |
Sex (female:male) | 72 | 40:32 | 13:16 | 27:16 | χ 2 = 1.6, p = 0.21 |
Age at death [years] | 71 | 72.8 ± 11.5 | 73.7 ± 10.5 | 72.2 ± 12.0 | U = 561.5, p = 0.6 |
Braak and Braak (IV:V:VI) | 72 | 8:14:50 | 1:7:21 | 7:7:29 | χ 2 = 3.2, p = 0.21 |
Thal phase (3:4:5) | 69a | 2:7:60 | 1:2:25 | 1:5:35 | χ 2 = 3.4, p = 0.5 |
TDP43 (neg:pos) | 49 | 25:24 | 14:7 | 11:17 | χ 2 = 2.6, p = 0.11 |
ApoE4 allele (neg:pos) | 66 | 28:38 | 15:12 | 13:26 | χ 2 = 2.4, p = 0.12 |
Age at death presented as mean ± first standard deviation; U two-sided Mann-Whitney U test; χ 2 chi-squared test. neg negative, pos positive. aThe three missing cases have a Thal phase ≥ 3
Table 2
Demographic and neuropathological overview of α-syn distribution subgroups
| | Avail. n | αSyn- | αSyn + A (amygdala pred.) | αSyn + B (brainstem pred.) | αSyn + C (cortical α-syn) | Statistic, p-value |
|---|
n (%) | 71 | 29 (41%) | 15 (21%) | 5 (7%) | 22 (31%) | |
Sex (female:male) | 39:32 | 13:16 | 11:4 | 2:3 | 13:9 | χ 2 = 3.9, p = 0.28 |
Age at death | 70 | 73.7 ± 10.5 | 76.5 ± 9.2 | 76.6 ± 5.6 | 68.5 ± 13.6 | K = 4.07, p = 0.25 |
Braak and Braak (IV:V:VI) | 71 | 1:7:21 | 0:1:14 | 2:2:1 | 5:4:13 | χ 2 = 15.5, p = 0.016b |
Thal phase (3:4:5) | 68a | 1:2:25 | 0:0:15 | 0:2:3 | 1:3:16 | χ 2 = 14.6, p = 0.27 |
TDP43 (neg:pos) | 49 | 14:7 | 5:6 | 3:2 | 3:9 | χ 2 = 5.6, p = 0.13 |
ApoE4 allele (neg:pos) | 65 | 15:12 | 6:7 | 1:4 | 5:15 | χ 2 = 5.5, p = 0.14 |
| Significant p-values were labeled in bold. Age at death presented as mean ± first standard deviation; K Kruskal-Wallis test; χ 2 chi-squared test; Avail. available, neg negative, pos positive, pred. predominant. aThe three missing cases have a Thal phase ≥ 3. bBraak and Braak staging was significantly different between αSyn- and αSyn + B (χ 2 = 8.7, p = 0.013), and αSyn + A and αSyn + B (χ 2 = 11.5, p = 0.003) |
α-syn load and distribution in AD
The αSyn- and αSyn + cases showed a comparable age distribution (independent t-test: p = 0.6). In the αSyn + group, there were comparatively more female than male subjects while the αSyn- group had a slight male predominance, although the difference was not significant (p = 0.21). There was no significant difference regarding the Braak and Braak staging (p = 0.21) or Thal phase (p = 0.5) between groups. There were proportionally more TDP43 positive cases in the αSyn + group, however, not significantly (p = 0.11). There were also more cases carrying at least one ApoE4 allele in proportion to non-carriers in the αSyn + group than in the αSyn- group but also not significantly (p = 0.12). Thus, there might be a female sex, TDP43, and ApoE4 bias in the α-syn positive group, even without reaching significance. However, it is not clear if this association is causally related or a limitation of the available cohort.
By definition, the αSyn + cases showed a higher α-syn load than αSyn- cases. Performing multiple linear regression correcting for the specific region name, age, and sex, there was a significant difference between αSyn + and αSyn- cases in cortical, subcortical, hippocampal, amygdala-entorhinal, and brainstem region clusters (Table 3), confirming the split into these two groups. The αSyn + cases showed the highest median α-syn load in the amygdala-entorhinal area, followed by the brainstem and hippocampal region, and low coverage in subcortical areas (Fig. 2). The α-syn load in cortical regions was low in the median but showed a large variability and thereby reached the highest values of covered area in single subjects. These findings suggest a region-dependent predestination for α-syn deposits in AD with a special focus on the amygdala in many cases, but a broad inter-patient variability.
As additional control analyses, we conducted multiple linear regression without correction for sex and age or with ApoE4 as an additional control factor. Both models showed significantly higher α-syn load in the αSyn + group in all region clusters. Furthermore, applying linear mixed-effects models with correction for age, sex, and a random factor for subject ID, only the difference in the amygdala-entorhinal region remained significant, indicating a strong difference in the amygdala (Supplementary Table 2). Comparing the α-syn load of αSyn + with αSyn- cases in 28 brain regions separately under correction of age and sex, there was a significantly higher α-syn covered area in the substantia nigra (p = 0.005), amygdala (p = 0.005), entorhinal cortex (p = 0.023), and olfactory bulb (p = 0.023), suggesting these regions as a focus of α-syn co-pathology in AD (Supplementary Table 3). Other brain regions, e.g., the hippocampus and insula cortex, are also affected. However, probably due to the small absolute numbers, the p-values were not statistically significant for other brain regions.
Table 3
Comparison of the α-syn covered area between αSyn + vs. αSyn- cases with multiple linear regression controlling for age and sex and correction for false discovery rate
Region cluster | n (αSyn-) | n (αSyn+) | Median [IQR] [%] of αSyn- cases | Median [IQR] [%] of αSyn + cases | β, p-value (age, sex corrected) |
|---|
cortical | 73 | 348 | 0.001 [0.0004; 0.003] | 0.09 [0.008; 0.48] | β = 0.009, p = 0.004 |
subcortical | 22 | 97 | 0.001 [0.0001; 0.003] | 0.06 [0.015; 0.35] | β = 0.003, p = 0.026 |
hippocampal | 30 | 176 | 0.0006 [0.0002; 0.002] | 0.17 [0.017; 0.61] | β = 0.005, p = 0.004 |
Amygdala-entorhinal | 54 | 70 | 0.0013 [0.0002; 0.006] | 0.67 [0.30; 1.4] | β = 0.013, p = < 0.001 |
brainstem | 36 | 72 | 0.0025 [0.0007; 0.006] | 0.19 [0.014; 0.72] | β = 0.005, p = 0.004 |
Significant p-values and the highest median α-syn covered area of the αSyn + group were labeled in bold. ID subject ID, IQR interquartile range
Comparing the sex distribution in the α-syn positive subgroups with the αSyn- cases, there was a female preponderance in αSyn + A, however, without reaching significance (p = 0.28). There was a trend towards younger age at death in αSyn + C with a mean age of 68.5 years (± 13.6 years standard deviation) in comparison with 76.5 ± 9.2 years in αSyn + A, 76.6 ± 5.6 years in αSyn + B, and 73.7 ± 10.5 years in αSyn- cases. Although this finding did not reach significance (p = 0.25) and there was a broad variability between cases, this observation suggests a negative association between cortically spread α-syn pathology in AD and survival. The Braak and Braak staging distribution was shifted towards lower Braak and Braak stages in αSyn + B, which reached significance when comparing αSyn- and αSyn + B (p = 0.013), as well as between αSyn + A and αSyn + B (p = 0.003). There was no significant difference regarding the Thal phases (p = 0.27). Where TDP43 information was available, two thirds of the αSyn- cases were also TDP43 negative while three quarters of the αSyn + C subgroup were TDP43 positive. Performing a chi-squared test over these groups, there was also no significant difference (p = 0.13). Regarding the presence or absence of the ApoE4 allele, 75% of the αSyn + C cases had at least one ApoE4 allele while it was more balanced in αSyn- and αSyn + A cases, although these group comparisons did not reach significance in a chi-squared test (p = 0.14).
To confirm that the α-syn distribution subgroups vary in their α-syn distributions, we applied multiple linear regression controlling for specific region names, sex, and age. Detailed results are presented in Fig. 2 and Supplementary Tables 6 and 9. In pairwise tests, all groups are significantly different from each other in their α-syn load across cortical regions, with the highest α-syn load in αSyn + C and, after a large gap, αSyn + B in second place. αSyn + C and to a lesser extent αSyn + B show significantly higher subcortical α-syn load than αSyn- cases. αSyn + C significantly shows the highest hippocampal and amygdala-entorhinal α-syn load, much higher than the actual amygdala-entorhinal predominant α-syn subgroup αSyn + A. αSyn + B and αSyn + C show higher brainstem α-syn loads than subgroup αSyn + A. Within αSyn + A, the highest α-syn load is in the amygdala and lower in other brain regions. αSyn + B shows the highest α-syn levels in the brainstem with low values in other brain regions, affirming its definition. Interestingly, αSyn + C manifests with an α-syn amygdala predominance next to high deposit loads in some cortical regions, and often a lower, but still high amount in other brain regions. The high deposit load in the amygdala in αSyn + C suggests a general α-syn sensitivity of the amygdala in AD, independent of the exact α-syn distribution type. In total, the identified distribution patterns propose the presence of distinct pathological α-syn accumulation features with overlaps, e.g., in the amygdala.
Tau load in relation to α-syn distribution
According to the inclusion criteria of Braak and Braak stage ≥ IV, all AD cases showed marked tau pathology. The most affected area was the amygdala-entorhinal region, followed by the hippocampal region and the cortical region in third place (Table 4). There was a low tau covered area in the brainstem and subcortical areas. To examine potential associations between tau and α-syn loads, we compared the tau covered area of αSyn- vs. αSyn + cases with multiple linear regression, correcting for the specific region name, age, and sex. Interestingly, there was no significant effect of α-syn presence on tau load in any brain region cluster (Fig. 3, Table 4), suggesting independent accumulation of α-syn and (AT8-) hyperphosphorylated tau.
As additional control analyses, we conducted multiple linear regression without correction for sex and age or with additional correction for ApoE4. Following the previous analysis, there was no significant difference regarding the tau load between αSyn- and αSyn + groups in all region clusters. Furthermore, applying linear mixed-effects models with correction for age, sex, and a random factor for subject ID also yielded no significant difference (Supplementary Table 2). Comparing the tau load of αSyn + vs. αSyn- cases in 28 brain regions separately under correction of age and sex, there was no significant difference (Supplementary Table 4). These findings support a theory of tau accumulation that is independent from α-syn deposits.
Table 4
Comparison of the tau covered area between αSyn + vs. αSyn- cases with multiple linear regression controlling for age and sex and correction for false discovery rate
Region cluster | n (αSyn-) | n (αSyn+) | Median [IQR] [%] of αSyn- cases | Median [IQR] [%] of αSyn + cases | β, p-value (age, sex corrected) |
|---|
cortical | 216 | 355 | 14.4 [7.7; 23] | 13.2 [5.3; 23.1] | β=-0.012, p = 0.19 |
subcortical | 39 | 72 | 1.5 [0.6; 4] | 1.7 [0.6; 5.6] | β=-0.003, p = 0.72 |
hippocampal | 133 | 222 | 19.8 [13.1; 27.5] | 18.9 [10.4; 27.3] | β=-0.017, p = 0.19 |
Amygdala-entorhinal | 42 | 63 | 26.3 [18.3; 33] | 19.8 [12; 30.6] | β=-0.043, p = 0.12 |
brainstem | 35 | 55 | 2.4 [1.7; 3.6] | 1.9 [01.3; 3] | β=-0.002, p = 0.72 |
The highest median tau covered areas of αSyn + and αSyn- groups were labeled in bold. ID subject ID, IQR interquartile range
To examine whether tau distribution varies between α-syn positive subgroups, we performed multiple linear regression controlling for specific region names, age, and sex (Fig. 3, Supplementary Tables 7 and 9). After FDR correction, there was a significantly decreased tau load in αSyn + B compared to αSyn- cases in cortical (p < 0.001), hippocampal (p < 0.001), and amygdala-entorhinal regions (p = 0.021). Furthermore, but only with correction for age and sex, there was a significantly higher tau load in αSyn + A than in αSyn- cases across cortical regions (p = 0.004), indicating a positive association between α-syn in the amygdala and cortical tau accumulation. On the other hand, there was a significantly lower tau load in αSyn + C than in αSyn- cases across cortical regions (p = 0.022), which was also only evident when controlling for age and sex, suggesting a relatively lower cortical tau load at death when cortical α-syn load is apparent. These findings were comparable with additional statistical correction for ApoE4 carriage, except for the lower tau load of αSyn + B in the amygdala-entorhinal region.
Aβ load in relation to α-syn distribution
The analyzed AD cases showed marked Aβ pathology, predominantly corresponding to Thal phase 5 (Table 1). The most affected areas were the parietal, frontal, and temporal cortices, followed by the amygdala, hippocampal, subcortical, and brainstem areas, which were impacted to a markedly lesser extent (Fig. 4, Table 5, refer to Supplementary Table 5 for results per region). To examine potential associations between Aβ and α-syn loads, we compared the Aβ covered area of αSyn- vs. αSyn + cases with multiple linear regression, correcting for the specific region name, age, and sex. There was a significantly higher Aβ load in cortical brain regions of αSyn + cases (Fig. 4, Table 5), suggesting an association of cortical Aβ with α-syn load.
As control analyses, we conducted multiple linear regression without correction for sex and age. Again, there was a significant difference regarding the Aβ load between αSyn- and αSyn + groups in cortical regions (Supplementary Table 2). Additionally, there were significantly higher Aβ covered areas in subcortical and hippocampal regions, suggesting a positive association between Aβ and α-syn across regions. Supplementing the multiple linear regression model with ApoE4 next to sex, age, and region name, the cortical Aβ load showed a trend but was not significantly different (p = 0.077). In a further control analysis, applying linear mixed-effects models with correction for age, sex, and a random factor for subject ID, there was also no significant difference, probably due to overcorrection (Supplementary Table 2). Regarding the 28 brain regions separately, the Aβ load was higher in the αSyn + vs. αSyn- group in the occipital sulcus, the insula cortex and the parahippocampal gyrus, however, these effects did not stay significant after FDR correction or after correction for age and sex (Supplementary Table 5). Thus, the increase of the Aβ load in αSyn + AD cases becomes particularly apparent when multiple regions are considered in one analysis, it is mostly evident in cortical areas, and the effect is partly explained by ApoE4 carriage.
Table 5
Comparison of the Aβ covered area between αSyn + vs. αSyn- groups with multiple linear regression controlling for age and sex and correction for false discovery rate
Region cluster | n (αSyn-) | n (αSyn+) | Median [IQR] [%] of αSyn- cases | Median [IQR] [%] of αSyn + cases | β, p-value (age, sex corrected) |
|---|
cortical | 149 | 288 | 3.9 [2.2; 6] | 5.4 [3.0; 9.7] | β = 0.017, p = 0.003 |
subcortical | 20 | 93 | 0.5 [0.2; 0.9] | 1.5 [0.4; 3.3] | β = 0.012, p = 0.09 |
hippocampal | 98 | 226 | 0.9 [0.2; 1.8] | 1.2 [0.3; 3.5] | β = 0.007, p = 0.09 |
Amygdala-entorhinal | 16 | 31 | 2.4 [1.5; 2.9] | 1.5 [0.9; 3.3] | β=-0.002, p = 0.86 |
brainstem | 26 | 32 | 0.7 [0.3; 1.3] | 0.7 [0.3; 1.6] | β<-0.001, p = 0.86 |
Significant p-values and the highest median Aβ covered area of αSyn- and αSyn + groups were labeled in bold. ID subject ID, IQR interquartile range
To examine whether the increased Aβ load can be attributed to specific α-syn positive subgroups, we applied multiple linear regression controlling for region names, age, and sex (Fig. 4, Supplementary Table 9). After FDR correction, there was a significantly increased Aβ load in αSyn + A compared to αSyn- cases across cortical regions (p = 0.037) and subcortical regions (p = 0.048). Additionally, there was a significantly increased Aβ load in αSyn + C compared to αSyn- (p = 0.01) across cortical regions, suggesting that the finding described above of more cortical Aβ in αSyn + cases is mainly driven by α-syn subgroups αSyn + A and αSyn + C. With additional correction for ApoE4 carriage, there was a significantly higher Aβ load in the cortical regions of the αSyn + A vs. αSyn- (p = 0.023) and αSyn + B (p = 0.0024) and in the hippocampal region of αSyn + A vs. αSyn + C (p = 0.010), supporting the notion of a particularly higher Aβ load in αSyn + A.
α-syn co-pathology in relation to age, sex, and ApoE genotype
We examined the association of α-syn co-pathology in AD with age at death, sex, and ApoE status (Fig. 5). In detail, we applied multiple linear regression with α-syn covered area as the target variable and sex as a predictor variable across region clusters, controlling for the specific region names. There was a higher α-syn load in cortical regions in female vs. male cases (β=-0.0049, p = 0.038), which did not remain significant after FDR correction (p = 0.19), suggesting a slight trend towards higher cortical α-syn load in female subjects. The results were comparable after additionally correcting for age. The α-syn load did not differ between female and male cases in other brain regions.
To examine the association with age, we defined three age groups: <65 years at death (< 65), 65 to < 75 years (65–75), 75 years or older (≥ 75). Thereby, it should be noted that all cases pertain to advanced stages of AD. We applied multiple linear regression with α-syn covered area as the target variable and age group as a predictor variable across region clusters, controlling for specific region names. Before FDR correction, there was a significantly higher cortical α-syn load in 65–75 aged AD patients compared with < 65 cases (β = 0.007, p = 0.033). This result was not significant after FDR correction or correction for sex. Interestingly, there was a significantly lower α-syn load in the hippocampal region in 65–75 aged AD patients compared with < 65 cases (β=-0.0033, p = 0.030), which was also significant after correction for sex but not after FDR correction. The amygdala-entorhinal α-syn load was significantly lower in ≥ 75 aged patients compared with 65–75 cases (β=-0.0078, p = 0.033), which was also significant after correction for sex but not after FDR correction. In total, these results suggest that α-syn co-pathology in general appears independent from patient age, but a higher hippocampal and amygdala-entorhinal α-syn load might be associated with a younger age at death to a certain extent. Another explanation could be that younger patients with initiated protein deposition cascades can accumulate higher α-syn loads in the hippocampus and amygdala until death, maybe due to fewer life-limiting comorbidities. However, this trend was not reflected in cortical regions.
In order to evaluate the association of α-syn load in AD with the ApoE genotype, we compared AD cases with at least one ApoE4 allele to cases without ApoE4. Again, multiple linear regression was applied with α-syn load as the target variable and ApoE status as the predictor variable across region clusters, controlling for the specific region names. The cortical α-syn load of ApoE4 carriers was significantly lower (β=-0.0055, p = 0.034) but did not remain significant after FDR correction or correction for sex and age. On the other hand, there was a significantly higher α-syn load in hippocampal (β = 0.0045, p = 0.0019) and amygdala-entorhinal regions (β = 0.0066, p = 0.038) of ApoE4 carriers which was significant after correction for age and sex but only the difference in the hippocampal regions stayed significant after FDR correction (p = 0.009 without and p = 0.016 with correction for age and sex). These results suggest ApoE4 as a risk factor for higher hippocampal and putatively amygdala-entorhinal α-syn load, which in turn might be associated with a younger age at death.
Aβ and tau load in relation to age, sex, and ApoE genotype
Additionally, we investigated the relation of age, sex, and ApoE genotype regarding tau and Aβ load (Supplementary Fig. 2 and Fig. 3). We applied multiple linear regression with tau or Aβ covered area as the target variable and sex, age, or ApoE as a predictor variable across region clusters, controlling for the specific region names and FDR-correction. There was a significantly higher tau load in male patients in the hippocampal (β = 0.029, p = 0.009) and amygdala-entorhinal regions (β = 0.06, p = 0.009), which was also significant after correction for age (p = 0.022, respectively). Conversely, the Aβ load was significantly higher in female patients in cortical (β=-0.022, p < 0.001) and hippocampal regions (β=-0.013, p < 0.001), which was significant after correction for age. These findings suggest a sex imbalance towards tau in male and Aβ in female cases.
Regarding different age groups, all with advanced disease stages, there was a significantly lower cortical tau load in the oldest group (≥ 75 years at death) than in the younger age groups, < 65 (β=-0.024, p < 0.001) and 65–75 (β=-0.024, p = 0.038). Both findings were significant after correction for sex. In line with this observation, there was a significantly higher Aβ load in the youngest age group, < 65 years, than in 65–75 years old patients in hippocampal regions (β=-0-016, p < 0.001) and in brainstem regions in 65–75 (β=-0.015, p = 0.007) and ≥ 75 years cases (β=-0.007, p = 0.0027). The findings remained significant after correction for sex and suggest a higher deposit load in younger AD cases at death.
Concerning the presence of at least one ApoE4 allele, there was no significant association with tau covered areas but with further age and sex correction, there was a significantly decreased tau load in ApoE4 carriers in cortical (β=-0.020, p = 0.03), hippocampal (β=-0.026, p = 0.03), and amygdala-entorhinal regions (β=-0.043, p = 0.049). Regarding Aβ, there was a higher Aβ load in cortical regions of ApoE4 carriers (β = 0.019, p = 0.001), also significant after age and sex correction. This finding coincides with the high Aβ load in the αSyn + C cases with a relatively high proportion of ApoE4 carriers. ApoE4 might be related to disseminated α-syn deposition and to a higher cortical Aβ load with a speculative causal relationship.