Demographic characteristics
In total, 288 patients with normal cognitive function, SCD, MCI, AD, or SVaD were included. The characteristics of the samples are listed in Additional file 1. The mean age of participants at baseline was 72.53 ± 6.86 years, and their estimated brain age was 75.58 ± 4.79 years. The proportions of participants with a clinical diagnosis of MCI and dementia were 69.2% and 14.9%, respectively. Additionally, 83.0% showed cognitive impairment below the global CDR of 0.5.
Identification and characteristics of microbial community types by DMM
Two DMM clusters were identified in the microbiota at the phylum level (Figure 1). The cluster 1 group was characterized by six genera with predominant commensal bacteria: Firmicutes, Actinobacteriota, Proteobacteria, Bacteroidota, Verrucomicrobiota, and Euryarchaeota. The cluster 2 group was characterized by seven genera with predominant commensal bacteria: Firmicutes, Bacteroidota, Actinobacteriota, Proteobacteria, Verrucomicrobiota, Euryarchaeota, and Cyanobacteria. The microbial community was mainly composed of Firmicutes, Bacteroidota, Actinobacteriota, and Proteobacteria, which is consistent with the results of previous studies [35]. However, Bacteroidota and Proteobacteria were more abundant in the cluster 2 group than in the cluster 1 group.
Figure 1. Dirichlet multinomial mixtures (DMM) clustering using microbiome data
The cluster 1 group is characterized by six genera with predominant commensal bacteria: Firmicutes, Actinobacteriota, Proteobacteria, Bacteroidota, Verrucomicrobiota, and Euryarchaeota. The cluster 2 group is characterized by seven genera with predominant commensal bacteria: Firmicutes, Bacteroidota, Actinobacteriota, Proteobacteria, Verrucomicrobiota, Euryarchaeota, and Cyanobacteria. Bacteroidota,and Proteobacteria are more abundant in the cluster 2 group than in the cluster 1 group * p<0.05
The characteristics between the cluster groups did not differ significantly in terms of sex, the CDR score, diagnosis, or clinical assessment (Table 1). However, the cluster 2 group was older (p=0.020), had lower MMSE scores (p=0.038), higher CDR-SB scores (p=0.028), higher brain age_residual (p=0.026), and a higher Bacteroidetes/Firmicutes ratio (p<0.001) than the cluster 1 group did.
Table 1. Clinical characteristics of participants according to the data-driven cluster groups using microbiome data
|
Group 1 (N = 191)
|
Group 2 (N = 97)
|
X or T
|
P value
|
N, or mean
|
%, or SD
|
N, or mean
|
%, or SD
|
Age
|
71.84
|
7.06
|
73.84
|
6.30
|
-2.34
|
0.020
|
Sex (female)
|
132
|
68.40
|
71
|
72.40
|
0.51
|
0.502
|
CDR
|
|
|
|
|
4.28
|
0.233
|
0
|
6
|
3.10
|
14
|
1.00
|
|
|
0.5
|
160
|
82.90
|
75
|
76.50
|
|
|
1
|
23
|
11.90
|
19
|
19.40
|
|
|
2
|
4
|
2.10
|
3
|
3.10
|
|
|
Dx
|
|
|
|
|
1.52
|
0.468
|
SCD
|
11
|
5.70
|
8
|
8.20
|
|
|
MCI
|
139
|
72.00
|
64
|
65.30
|
|
|
Dementia
|
43
|
22.30
|
26
|
26.50
|
|
|
CDR-SB
|
2.25
|
0.15
|
2.86
|
0.23
|
4.89
|
0.028
|
MMSE
|
24.11
|
0.33
|
22.85
|
0.49
|
4.38
|
0.038
|
BMI
|
23.99
|
0.30
|
24.13
|
0.46
|
0.06
|
0.801
|
MNA
|
20.74
|
0.32
|
21.08
|
0.48
|
0.33
|
0.567
|
Brain age residual
|
-0.15
|
0.09
|
0.20
|
0.13
|
5.04
|
0.026
|
Bacteroidota/Firmicutes ratio
|
0.027
|
0.075
|
0.473
|
0.475
|
-9.188
|
< 0.001
|
Abbreviations: SD, standard deviation; CDR-SB, Clinical Dementia Rating Sum of Box; MMSE, Mini-Mental State Examination; BMI, body mass; MNA, Mini Nutritional Assessment; SCD, subjective cognitive decline; MCI, mild cognitive impairment
Elastic Net Predictions
The best model achieved an out-of-bag accuracy of 0.66 predicting groups (α=0, λ=68.9, p<0.001). The results for each individual factor are presented in Table 2, which shows the pooled p-values across the imputations. The cluster 1 group was associated with higher anxiety, lifetime alcohol consumption, and a higher likelihood of thyroid disorders, compared with the cluster 2 group. However, individuals with the cluster 1 group were younger, had lower brain age, lower global amyloid, higher overall cognitive function (working memory, spatial memory, word recall, and overall cognitive function), and higher hematocrit, hemoglobin, HDL, and vitamin B12 levels, compared with individuals with the cluster 2 group.
Table 2. Factors associated with cluster groups identified by elastic net regression
|
fcoeff average
|
fcoeff SD
|
P vales pooled
|
Compared to cluster 2, cluster 1 has
|
Chronological age
|
0.0010
|
1.77E-05
|
7.200E-06
|
lower age
|
Brain age_residual
|
0.0009
|
1.86E-06
|
3.520E-05
|
lower brain age
|
Global_amyloid_SUVR
|
0.0008
|
6.95E-05
|
0.001
|
lower amyloid
|
CDR-SB
|
0.0006
|
7.54E-06
|
0.019
|
better overall cognitive function
|
BAI
|
-0.0010
|
1.98E-06
|
2.090E-06
|
higher anxiety
|
Alcohol consumption
|
-0.0009
|
2.48E-06
|
1.080E-05
|
higher lifetime drinking
|
RCFT_DR_ z score
|
-0.0007
|
2.43E-05
|
0.002
|
higher spatial memory
|
RCFT_Copy_z score
|
-0.0006
|
6.43E-06
|
0.031
|
higher spatial memory
|
BNT_ z score
|
-0.0007
|
1.42E-05
|
0.007
|
higher naming function (word recall)
|
COWAT_Animal_z score
|
-0.0008
|
1.09E-05
|
0.001
|
higher working memory
|
Hematocrit
|
-0.0007
|
2.00E-05
|
0.002
|
higher hematocrit
|
Hemoglobin
|
-0.0007
|
2.06E-05
|
0.006
|
higher hemoglobin
|
HDL
|
-0.0007
|
3.31E-05
|
0.007
|
higher HDL
|
Vitamin B12
|
-0.0007
|
2.47E-05
|
0.009
|
higher Vitamin B12
|
History of thyroid disease
|
-0.0007
|
1.35E-06
|
0.014
|
higher likelihood of thyroid disorder
|
Abbreviations: SD, standard deviation; SUVR, standardized uptake value; CDR-SB, clinical dementia rating sum of box; BAI, Beck anxiety inventory; SCD, subjective cognitive decline; RCFT_DR, Rey complex figure delayed recall test; BNT, Boston naming test; COWAT, controlled oral word association test; HDL, High-Density Lipoprotein
Mediational analyses
We explored the mediating roles of brain age between group clusters and cognitive function (measured using MMSE and CDR-SB). For the mediation analysis, we divided the data into two clusters, coding cluster 1 as “0” and cluster 2 as “1.” Mediation analyses results are shown in Figure 2, Table 3, and Table 4. Table 3 presents the estimated regression coefficients for the clusters, brain age, and the MMSE score. As observed, individuals with the cluster 2 group exhibited greater brain age than did those with the cluster 1 group (ß =0.283, 95% confidence interval [CI]=0.044 to 0.523, p=0.021). Furthermore, when controlling for the cluster group, as a mediator, brain age significantly affected the MMSE score (ß=-1.244, 95% CI=-1.800 to -0.687, p<0.001). Although the direct effect was insignificant, the indirect effect of the cluster group on the MMSE score through brain age was significant (ß=-0.352, 95% CI=-0.690 to -0.015, p=0.041). The mediation analysis, adjusted for sex, age, and education, was consistent with the unadjusted analysis.
Figure 2. Mediation analysis among cluster group, brain age, and cognitive function measures
Individuals with the cluster 2 group exhibit greater brain age than do those with the cluster 1 group (ß =.436, 95% confidence interval [CI]=0.118 to 0.883, p=0.024). When controlling for the cluster group, as a mediator, brain age significantly affects the MMSE score (ß=-1.288, 95% CI=-1.793 to -0.783, p<0.001). Although the direct effect was insignificant, the indirect effect of the cluster group on the MMSE score through brain age is significant (ß=-0.370, 95% CI=-0.714 to -0.026, p=0.035). For the CDR-SB score outcome variable, individuals with the cluster 2 group exhibit greater brain age than do those with the cluster 1 group (ß=0.436, 95% CI=0.118 to 0.883, p=0.024). When controlling for the cluster group, as a mediator, brain age significantly affects the CDR-SB score (ß=0.425, 95% CI=0.189 to 0.661, p<0.001). Both direct and indirect effects of the cluster group on the CDR-SB score through brain age are not significant; however, the indirect effect shows the same tendency as the MMSE score outcome at the p value of 0.052 (direct: ß=0.301, 95% CI=-0.198 to 0.801, p=0.238; indirect: ß=0.122, 95% CI=-0.001 to -0.246, p=0.052, respectively). Sex, age, and education are adjusted for in the analysis.
Abbreviations: CDR-SB, Clinical Dementia Rating Sum of Box; MMSE, Mini-Mental State Examination; * p<0.05; **p<0.01
Table 3. Mediation analysis among cluster groups, brain age, and MMSE
|
Coeff (ß)
|
SE
|
T
|
P value
|
LLCI
|
ULCI
|
Unadjusted analysis
|
|
|
|
Outcome: Brain age residual
|
Cluster group
|
0.283
|
0.122
|
2.316
|
0.021*
|
0.044
|
0.523
|
Outcome: MMSE
|
Cluster group
|
-0.195
|
0.595
|
-0.328
|
0.743
|
-1.361
|
0.971
|
Brain age residual
|
-1.244
|
0.284
|
-4.382
|
<0.001*
|
-1.800
|
-0.687
|
Indirect effect
|
-0.352
|
0.172
|
-2.047
|
0.041*
|
-0.690
|
-0.015
|
Direct effect
|
-0.195
|
0.595
|
-0.328
|
0.743
|
-1.361
|
0.971
|
Total effect
|
-0.547
|
0.609
|
-0.899
|
0.368
|
-1.740
|
0.646
|
Adjusted analysis
|
|
|
|
Outcome: Brain age residual
|
Cluster group
|
0.436
|
0.193
|
2.262
|
0.024*
|
0.118
|
0.883
|
Sex
|
-0.004
|
0.127
|
-0.035
|
0.972
|
-0.254
|
0.245
|
Age
|
-0.002
|
0.009
|
-0.196
|
0.845
|
-0.019
|
0.015
|
Education
|
0.005
|
0.013
|
0.360
|
0.719
|
-0.020
|
0.029
|
Outcome: MMSE
|
Cluster group
|
-0.195
|
0.545
|
-0.358
|
0.720
|
-1.264
|
0.873
|
Brain age residual
|
-1.288
|
0.258
|
-5.000
|
<0.001*
|
-1.793
|
-0.783
|
Sex
|
-0.035
|
0.557
|
-0.063
|
0.950
|
-1.126
|
1.056
|
Age
|
0.047
|
0.038
|
1.245
|
0.213
|
-0.027
|
0.122
|
Education
|
0.427
|
0.055
|
7.777
|
<0.001*
|
0.319
|
0.534
|
Indirect effect
|
-0.370
|
0.176
|
-2.109
|
0.035*
|
-0.714
|
-0.026
|
Direct effect
|
-0.195
|
0.545
|
-0.358
|
0.720
|
-1.264
|
0.873
|
Total effect
|
-0.565
|
0.563
|
-1.004
|
0.315
|
-1.669
|
0.538
|
Abbreviations: MMSE, mini mental status examination; SE, standard error; LLCI, lower limit confidence interval; ULCI, upper limit confidence interval. * p<0.05
Table 4. Mediation analysis among cluster groups, brain age, and CDR-SB
|
Coeff (ß)
|
SE
|
T
|
P value
|
LLCI
|
ULCI
|
Unadjusted analysis
|
|
|
|
Outcome: Brain age residual
|
Cluster group
|
0.283
|
0.122
|
2.316
|
0.021*
|
0.044
|
0.523
|
Outcome: CDR-SB
|
Cluster group
|
0.288
|
0.252
|
1.140
|
0.254
|
-0.207
|
0.783
|
Brain age residual
|
0.425
|
0.120
|
3.524
|
<0.001*
|
0.188
|
0.661
|
Indirect effect
|
0.120
|
0.062
|
1.935
|
0.053
|
-0.002
|
-0.242
|
Direct effect
|
-0.195
|
0.595
|
-0.328
|
0.743
|
-1.361
|
0.971
|
Total effect
|
-0.547
|
0.609
|
-0.899
|
0.368
|
-1.740
|
0.646
|
Adjusted analysis
|
|
|
|
Outcome: Brain age residual
|
Cluster group
|
0.436
|
0.193
|
2.262
|
0.024*
|
0.118
|
0.883
|
Sex
|
-0.004
|
0.127
|
-0.035
|
0.972
|
-0.254
|
0.245
|
Age
|
-0.002
|
0.009
|
-0.196
|
0.845
|
-0.019
|
0.015
|
Education
|
0.005
|
0.013
|
0.360
|
0.719
|
-0.020
|
0.029
|
Outcome: CDR-SB
|
Cluster group
|
0.301
|
0.255
|
1.182
|
0.237
|
-0.198
|
0.800
|
Brain age residual
|
0.425
|
0.120
|
3.533
|
<0.001*
|
0.189
|
0.661
|
Sex
|
0.108
|
0.260
|
0.417
|
0.677
|
-0.401
|
0.618
|
Age
|
-0.006
|
0.018
|
-0.332
|
0.740
|
-0.041
|
0.029
|
Education
|
-0.014
|
0.026
|
-0.550
|
0.583
|
-0.064
|
0.036
|
Indirect effect
|
0.122
|
0.063
|
1.943
|
0.052
|
-0.001
|
0.246
|
Direct effect
|
0.301
|
0.255
|
1.181
|
0.238
|
-0.198
|
0.801
|
Total effect
|
0.423
|
0.258
|
1.641
|
0.101
|
-0.082
|
0.929
|
Abbreviations: CDR-SB, clinical dementia rating sum of box; SE, standard error; LLCI, lower limit confidence interval; ULCI, upper limit confidence interval. * p<0.05
For the CDR-SB score outcome variable in Table 4, individuals with the cluster 2 group exhibited greater brain age than did those with the cluster 1 group (ß=0.283, 95% CI=0.044 to 0.523, p=0.021). Furthermore, when controlling for the cluster group, as a mediator, brain age significantly affected the CDR-SB score (ß=0.425, 95% CI=0.188 to 0.661, p<0.001). In contrast, both direct and indirect effects of the cluster group on the CDR-SB score through brain age were not significant; however, the indirect effect showed the same tendency as the MMSE score outcome at the p value of 0.053 (direct: ß=-0.195, 95% CI=-1.361 to 0.971, p=0.743; indirect: ß=0.120, 95% CI=-0.002 to -0.242, p=0.053, respectively). The results of the mediation analysis, adjusted for sex, age, and education, was consistent with those of the unadjusted analysis.
Explorative analyses for mediation effects of plasma amyloid A40, A42, p-tau 217 and PSMD did not show any significant results (Additional file 1).