Demographic, genetic and molecular biomarkers, and neuropsychological data
There were no significant differences in age, gender, education, levels of thirty-eight serum lipid metabolites, or any candidate genotypes, except apolipoprotein E (APOE) genotypes within any groups of participants. Significant cognitive decline as signed by lower MMSE scores and higher ADAS-cog scores, gradually decreased Aβ level, and increased Tau and p-Tau levels were identified in ADS individuals compared to the CN. In terms of cognitive scores or CSF biomarkers, there were no significant differences between early MCI (EMCI) and late (LMCI) groups. No genotypes deviated from the Hardy-Weinberg equilibrium with all p values above 0.05. More details of demographic, lipid pathway-based genotypes, CSF biomarkers, and global cognition are displayed in Table 1 and Table S1 and S2.
Table 1
Demographic data, lipid pathway-based genotypes, cerebrospinal fluid biomarkers, and global cognitive performance across the AD spectrum.
|
CN
|
EMCI
|
LMCI
|
AD
|
P
values
|
Items
|
(n = 51)
|
(n = 26)
|
(n = 26)
|
(n = 21)
|
Age (years)
|
74.08 ± 5.79
|
70.04 ± 6.87
|
70.81 ± 7.14
|
71.81 ± 7.77
|
0.051
|
Gender (F/M)
|
30/21
|
14/12
|
11/15
|
9/12
|
0.447*
|
Education (years)
|
16.31 ± 2.59
|
15.27 ± 2.51
|
16.31 ± 2.51
|
15.14 ± 2.76
|
0.157
|
Multiple protective genes
|
CLU T status (TC + TT/CC)
|
38/13
|
17/9
|
18/8
|
13/8
|
0.710*
|
LDLR A status (AG + AA/GG)
|
38/13
|
16/10
|
15/11
|
11/10
|
0.240*
|
LRP1 T status (TC + TT/CC)
|
15/36
|
8/18
|
5/21
|
10/11
|
0.214*
|
PICALM A status (AG + AA/GG)
|
27/24
|
15/11
|
16/10
|
11/10
|
0.884*
|
Multiple risk genes
|
APOE ε4 status (+/−)
|
15/36
|
14/12
|
12/14
|
15/6
|
0.008*
|
SORL1 G status (TG + GG/TT)
|
19/32
|
12/14
|
12/14
|
6/15
|
0.548*
|
CETP A status (AG + AA/GG)
|
46/5
|
24/2
|
23/3
|
19/2
|
0.974*
|
ABCA1 G status (AG + GG/AA)
|
45/6
|
25/1
|
21/5
|
19/2
|
0.369*
|
BIN1 C status (TC + CC/TT)
|
27/24
|
17/9
|
15/11
|
11/10
|
0.739*
|
Cerebrospinal fluid biomarkers
|
Aβ (pg/ml)
|
192.79 ± 50.17bc
|
183.61 ± 50.66d
|
168.77 ± 50.80
|
140.40 ± 43.59
|
0.001
|
Tau (pg/ml)
|
68.53 ± 34.14c
|
79.32 ± 51.89d
|
86.01 ± 52.19e
|
129.29 ± 61.42
|
< 0.001
|
pTau (pg/ml)
|
34.18 ± 16.58bc
|
39.60 ± 24.73d
|
48.42 ± 33.50
|
55.23 ± 26.13
|
0.005
|
Global cognitive performance
|
MMSE
|
28.84 ± 1.16abc
|
27.92 ± 2.13d
|
27.73 ± 1.54e
|
22.67 ± 2.50
|
< 0.001
|
ADAS-Cog
|
10.76 ± 6.53abc
|
14.19 ± 6.58d
|
16.96 ± 5.32e
|
35.81 ± 8.99
|
< 0.001
|
Notes: *, p values were obtained using a Chi-square test; other p values were obtained from a one-way ANOVA. Unless indicated, data are presented as means ± standard deviation. Post hoc analyses were used with least significance difference correction (p < 0.05): a: statistical difference detected between CN group and EMCI group; b: statistical difference was detected between CN group and LMCI group; c: statistical difference was detected between CN group and AD group; d: statistical difference was detected between EMCI group and AD group; e: statistical difference was detected between LMCI group and AD group. Abbreviations: CN, cognitively normal; EMCI, early mild cognitive impairment; LMCI, late mild cognitive impairment; AD, Alzheimer’s disease; M/F, male/female; CLU, clusterin; LDLR, low density lipoprotein receptor; LRP1, low density lipoprotein receptor-related protein 1; PICALM, phosphatidylinositol-binding clathrin assembly protein; APOE, apolipoprotein E; SORL1, sortilin-related receptor 1; CETP, cholesterol ester transfer protein; ABCA1, ATP-binding cassette transporter A1; BIN1, bridging integrator 1; Aβ, amyloid-1 to 42 peptide; Tau, total tau; pTau, tau phosphorylated at the threonine 181 position; MMSE, mini-mental state examination; ADAS-Cog, Alzheimer’s Disease Assessment Scale-Cognitive Subscale. |
Relationships among polygenic scores, lipid metabolites, CSF biomarkers and cognitive performance
First, binomial nonlinear connections were discovered between five cholesterol metabolism related biomarkers in the serum and ADAS-cog score, including serum total cholesterol (SERUM_C), esterified cholesterol (ESTC), free cholesterol (FREEC), phosphatidylcholine (PC), and sphingomyelins (SM). Besides, TOTPG (total phosphoglycerides), total choline (TOTCHO), and S_HDL_P (total lipids, phospholipids, total cholesterol, cholesterol esters, free cholesterol and triglycerides in small HDL particles) levels represented a correlation trend with ADAS-cog or MMSE scores. Then, three cholesterol metabolism related markers including SERUM_C, ESTC and SM were also linearly correlated with Tau but not Aβ and p-Tau levels of CSF in the ADS. In addition, linear regression between polygenic scores and CSF biomarkers disclosed that genetic risk scores (GRS) was significantly correlated with Aβ, Tau, and even p-Tau levels, while genetic protective score (GPS) was not correlated with any of CSF biomarkers. It is interesting that relative risk scores (RRS = GRS - GPS) was significantly influenced the Aβ, Tau but not p-Tau levels in the ADS. Of note, the correlations between GRS_n (GRS without APOE), RRS_n (RRS without APOE) and CSF markers were not found, so the graphs were not presented here. All the corresponding graphs above were plotted in Fig. 1. Meanwhile, regression analyses revealed that CSF biomarkers (including Aβ, Tau and p-Tau) could significantly impact global cognitive performance in a nonlinear manner (Fig. S1).
Dynamic trajectory of large-scale brain network roles across the AD spectrum
To explore the dynamic trajectory both within- and between RSNs in ADS patients, pairwise functional connections (correlations) were extracted within- and between ten predefined large-scale functional brain networks: auditory network (AUD), cingulo-opercular network (CON), dorsal attention network (DAN), DMN, fronto-parietal network (FPN), SAN, sensory network (SMN), subcortical network (SUB), ventral attention network (VAN), and visual network (VIS), as derived from the brain atlas of Power et al. [33]. By mapping the group mean within-network connectivity (WNC) and between-network connectivity (BNC) to a 2D parameter space, the mean functional role of 10 RSNs was qualitatively described across the ADS (Fig. 2A-D). From the means of individual WNC and BNC values (depicted by horizontal and vertical dotted lines in Fig. 2F, detailed information provided in Supplemental Materials), the RSNs from the CN group were consequently classified into four network roles: cohesive connector (SAN, DAN, SMN, and SUB), cohesive provincial (VIS), incohesive connector (AUD, FPN and CON), or incohesive provincial (DMN and VAN) (Fig. 2A). In addition to DMN and VAN, which exhibited both weaker cohesive connector and cohesive provincial roles, the other eight networks in the ADS represented divergent network roles compared to those of the CN group. Specifically, SAN, DAN, and AUD represented incohesive provincial and connector networks in patients with EMCI, LMCI and AD, respectively, the converse of that observed in CN subjects (Fig. 2A-D). The graphs visually demonstrated how the network roles of these large-scale RSNs dynamically changed with severity of disease (Fig. 2E). Interestingly, the strengths of WNC and BNC exhibited a dynamically weakened trend, except for SUB, as disease progressed through the ADS, indicating that spatiotemporal patterns of large-scale RSNs represent a rebounding network mode rather than cascading network failure, as described previously [16].
Group-level Comparison Of Network Connectivity In Ad Spectrum Individuals
We next tested differences in WNC and BNC in terms of large-scale RSNs among the four groups. Firstly, we obtained distinctive WNC and BNC matrices of the 10 RSNs for the four groups (Fig. 2G). Clearly, five RSNs (DAN, FPN, SAN, VAN, and VIS) exhibited significantly differential WNC among the disease spectrum (Fig. 2H). Although VAN and DMN were found to be incohesive provincial networks in four groups (Fig. 2A-D), VAN exhibited significantly lower WNC and BNC across the ADS (Fig. 2A-D). Similarly, five RSNs, including DAN, SAN, SMN, CON, and AUD, were found to have incohesive connector roles and provincial networks that had lower connectivity in the ADS relative to CN subjects. It was noted that the FPN shifted from an incohesive connector to incohesive provincial network while the VIS changed from cohesive connector to a cohesive provincial network, both representing lower connectivity in the ADS relative to CN subjects. In addition, the SUB displayed more cohesive connector and cohesive provincial networking, having the greatest connectivity in ADS individuals compared with CN subjects. Furthermore, ADS patients also showed significantly differential one-versus-all-other-network connectivity in the DAN, FPN, SAN, VAN, and VIS networks compared with controls (Fig. 2H).
In addition, pairwise BNC was calculated as the mean connectivity between each pair of RSN. Connectivity profiles of patients with EMCI, LMCI and AD were compared with controls. Figure 2H demonstrates that pairwise BNC was significantly different for ADS patients in the following pairs: AUD-VAN, AUD-VIS, CON-DAN, CON-VAN, DMN-FPN, DMN-VAN, DAN-FPN, DAN-SMN, FPN-VAN, SUB-VIS and VAN-VIS. Furthermore, post hoc analysis indicated that the source of the significant differences in these pairwise BNC groups was at the large-scale network level. Specifically, ADS patients were characterized by continuous hypoconnectivity and dynamically hyperconnected links among the ten predefined RSNs as disease progressed (Fig. 2I and 2J). These original alterations of large-scale networks may initially reproduce those spatiotemporal pattern discrepancies, accounting for proposed molecular pathophysiological mechanisms at the distributed network level.
Correlation patterns of large-scale network connectivity with CSF biomarkers and cognitive performance in the AD spectrum
To explore the potential relationship between the dynamic trajectory of connectivity of the RSNs and CSF biomarkers or cognitive performance, a new method of combination network analysis and CCA was utilized. Recent studies have demonstrated that CCA, a powerful multivariate approach that seeks to identify clusters of maximal correlation between two groups of variables, can detect associations between structural or functional connectivity and other phenotypic measures [33, 34]. Using this method, we demonstrated that large-scale brain network abnormalities were significantly correlated with phenotypic variations and molecular biomarkers in the ADS. In the first step, univariate correlation was used to test the composition of the clinical CCA mode with each of the two clinical variables (MMSE and ADAS-cog). We observed that clinical CCA mode was highly correlated with MMSE score (r = 0.78, p < 0.0001) and ADAS-cog score (r = 0.78, p < 0.0001) (Fig. 3A). Similarly, we identified that CSF CCA mode was highly correlated with levels of Tau (r = 0.72, p < 0.0001) and pTau (r = 0.68, p < 0.0001), and moderately correlated with Aβ42 (r = 0.55, p < 0.0001) (Fig. 3B). As shown in Fig. 3C, the network CCA mode was significantly correlated with 55 original network variables (Table S4). In total, the CCA model of network was significantly correlated with clinical variate CCA (Fig. 3D, r = 0.93, p < 0.0001) and CSF CCA (Fig. 3E, r = 0.95, p < 0.0001) modes, respectively.
Correlation patterns of lipid pathway-based genetic variants and lipoproteins with large-scale network connectivity in AD spectrum patients
Similarly, we also performed CCA to ascertain the association of brain network connectivity measures with accumulated lipid-related genetic scores and lipoproteins in the blood of ADS patients. We firstly tested univariate correlations for each of the 3 gene variables and 38 serum lipid variables in order to better understand the composition of gene CCA and serum lipid CCA modes. We found that first gene CCA mode was highly correlated with GRS (r = 1, p < 0.0001), GPS (r = 1, p < 0.0001) and RRS (r = 1, p < 0.0001) (Fig. 4A). As shown in Fig. 4B, serum lipid CCA mode was significantly correlated with all 38 original serum lipid variables (Table S5). The results of third pair CCA mode of network and gene variate were again highly significant (Fig. 4C, r = 0.97, p < 0.0001), as was fourth pair CCA mode of network and serum lipid variable (Fig. 4D, r = 0.82, p < 0.0001).
In order to determine the potential for APOE ε4 genotype to alter the association between lipid metabolism-related genes and dynamic changes in RSNs, we constructed a second gene CCA mode and found that second order gene CCA mode was highly correlated with GRS_n (r = 1, P < 0.0001) and RRS_n scores (r = 1, P < 0.0001) after removing the APOE ε4 genotype (Fig. S2A). Fifth pair CCA mode of the network and three gene score variables where removed APOE ε4 OR values was also significantly correlated (Fig. S2C, r = 0.94, p < 0.0001).
Post hoc analysis revealed the potential of distinctive lipid-related genetic scores and lipoproteins on large-scale network connectivity, CSF biomarkers, and cognitive performance
To determine the direction and magnitude of these associations between network CCA mode and a single variate of a clinical indicator, we conducted post hoc correlation analysis. As illustrated in Fig. 5A, nineteen lipoproteins were mostly associated with increased network connectivity and seven with decreased network connectivity within- and between- ten predefined RSNs. Furthermore, GRS was positively associated with increased network connectivity within the SUB and negatively associated with altered network connectivity between CON-VAN, DAN-VAN, FPN-VAN, SAN-VAN, SUB-VAN, VAN-VIS, and AUD-VAN, while GPS was negatively associated with decreased network connectivity between the FPN and SUB. Similarly, RRS was positively associated with decreased network connectivity between DMN-SUB, DAN-SUB, FPN-SUB, CON-DMN, AUD-CON, and negatively associated with decreased network connectivity between AUD-VAN, and DAN-VAN, whereas RRS was associated with increased network connectivity within the SUB. In addition, MMSE was negatively correlated with decreased network connectivity between DMN and SAN, while ADAS-cog and Tau were mostly associated with increased network connectivity between SAN-SUB, and DMN-SAN, as shown in Fig. 5B. It is noteworthy that GRS was only associated with decreased network connectivity between SAN-VAN after removing the effects of the APOE ε4 genotype (Fig. S3). Detailed information for these correlation coefficients (r) and p values are described in Table S6.
SVM analyses identified potential lipid-associated imaging biomarker for AD spectrum
After post-hoc analysis, we performed correlation analysis, and found that there were six functional connections significantly correlated to all nineteen lipoproteins and sixteen functional connections related to all three gene scores (Table S7). Then, these twenty-two features were used for classification. The SVM model revealed that the lipid-associated twenty-two functional connections represented higher capacity to discriminate disease spectrum (AUC between 0.82–0.92), as shown in Fig. 6.