Demographic and clinical data
Table 1 summarized the basic demographic and clinical characteristics. There were no statistically significant differences in age and gender between-group. The mean MMSE score of HCs was 28.59 ± 1.01 and the range was 27–30, the score of CSVD was 28.30 ±0.89 and the range was 27–30. The mean MoCA score of HCs was 27.62 ± 1.27 and the range was 26–30, the score of CSVD was 27.38 ± 1.05 and the range was 26–30. There were no significant differences in MMSE and MoCA scores between the two groups.
Table 1 Baseline data of subjects
|
Basic information
|
CSVD
(n=78)
|
HC
(n=49)
|
t/x2
|
P
|
Demographics
|
|
|
|
|
Age (year,x̄ ± s )
|
66.18 ± 1.42
|
66.53 ± 1.30
|
1.396
|
0.165
|
Male (n,%)
|
48(61.54)
|
29(59.18)
|
0.070
|
0.791
|
Education (year,x̄ ± s )
|
9.47 ± 1.94
|
9.88 ± 2.00
|
1.146
|
0.254
|
Cognitive function scale(score,x̄ ± s )
|
|
|
|
|
MMSE
|
28.30 ± 0.89
|
28.59 ± 1.01
|
1.696
|
0.092
|
MoCA
|
27.38 ± 1.05
|
27.62 ± 1.27
|
1.155
|
0.250
|
Behavioral performance
WM task accuracy and RT were analyzed between the two groups. The correct rate was significantly lower than the HC group(68.26±1.81 vs 83.5±3.39), shown in Figure 1B. Response time in CSVD group was significantly longer than that in HC group༈1756.48±42.8ms vs 1262.48±41.79ms༉, shown in Figure 1C, The difference was statistically significant༈all P<0.01༉, which suggested that WM behavior disorder existed in CSVD group.
Time-frequency analysis
Apply the short-time Fourier transform to obtain the time-frequency distribution of the EEGs of the characteristic channel during the WM task. During the correct trials, the time-frequency distribution of Fz in the HC group is shown in Fig. 2A. The diagram shows that the energy density of theta band is most obvious when executing tasks. During the correct trials, the time-frequency distribution of Fz in the CSVD group is shown in Fig. 2B. The diagram shows that the energy density of theta band is most obvious when executing tasks. Compared the power of different bands: in Fig. 2C, the power of theta band was higher than other bands in the HC group (theta/delta: t = 5.188, P < 0.001; theta/alpha: t = 5.522, P < 0.001; theta/beta: t = 48.37, P < 0.001; theta/gamma: t=53.06, P< 0.001), in Fig. 2D, the power of theta band was higher than other bands in the CSVD group (theta/delta: t = 4.782, P < 0.001; theta/alpha: t = 5.758, P < 0.001; theta/beta: t = 55.15, P < 0.001; theta/gamma: t=41.42, P< 0.001). Then, we compared the power of the five bands between the two groups (Fig. 2E). Compared with the HC group, the CSVD group exhibited significantly decreased theta band (t = 13.43, P < 0.001), beta band (t = 14.60, P < 0.001), gamma band (t = 3.763, P < 0.001) and alpha band (t = 2.691, P < 0.05). The delta band was no statistical difference (t = 1.65, P = 0.0.122) between the two groups. Fig. 2F shows the variation curve of whole brain average power density with frequency during the WM tasks for the HC group. Fig. 2G shows the variation curve of whole brain average power density with frequency during the WM tasks for the CSVD group. Fig. 2H shows the comparation curve of whole brain average power density with frequency during the WM tasks between the two groups. In Fig. 2I, we compared the 25 channel average power in theta band at each phase of WM task in the HC group, the power of encoding phase was lower than other phase (encoding/maintenance: t = 3.068, P < 0.01; encoding/retrieval: t = 4.009, P < 0.001). The power was no statistical difference between the maintenance phase and the retrieval phase in the HC group (t = 1.539, P = 0.125). In Fig. 2J, we compared the 25 channel average power in theta band at each phase of WM task in the CSVD group, the power of retrieval phase was higher than other phase (encoding/retrieval: t = 8.321, P < 0.01; maintenance/retrieval: t = 7.431, P < 0.001). The power was no statistical difference between the encoding phase and the maintenance phase in the CSVD group (t = 1.018, P = 0.309). Then, we compared the 25 channel average power in theta band at each phase between the two groups (Fig. 2K). The CSVD group exhibited significantly decreased encoding phase (t = 26.09, P < 0.001), maintenance phase (t = 24.52, P < 0.001) and retrieval phase (t = 26.73, P < 0.05) compared to that of the HC group.
Figure 3A, D, G show the energy topography of the theta band among the three phases in the HC group. Fig. 3B, E, H show the energy topography of the theta band among the three phases in the CSVD group. Both in the two groups, the topographic map shows the theta activity related to the frontal region, especially in the Fz channel in the midline area of the frontal region. We compared the power of the four brain regions between the two groups. In Fig. 3C, compared with the HC group, significantly decreased theta power during the encoding phase of the CSVD group in frontal region (t = 6.749, P < 0.001), parietal region (t = 3.923, P < 0.01), temporal region (t = 11.29, P < 0.001) and occipital region (t = 11.21, P < 0.01). From Fig. 3F, Compared with the HC group, significantly decreased theta power during the maintenance phase of the CSVD group in frontal region (t = 6.766, P < 0.001), parietal region (t = 3.067, P < 0.01) and temporal region (t = 12.23, P < 0.001). There was no statistical difference in occipital region (t = 3.557, P = 0.071). From Fig. 3I, Compared with the HC group, significantly decreased theta power during the retrieval phase of the CSVD group in frontal region (t = 6.78, P < 0.001), parietal region (t = 2.958, P < 0.05), temporal region (t = 10.15, P < 0.001) and occipital region (t = 21.38, P < 0.01) .
Functional connectivity across electrodes during WM tasks
The DTF value was used to evaluate the FC strength (theta band) during the WM task in the HC and CSVD groups. In Fig. 4, showed three phases connection matrix during the WM tasks of the two groups. Both in the two groups, DTF was mainly concentrated in the frontal region, especially in the Fz channel. However, Compaired to the HC group (Figure 4A, D, G), the FC strength of the CSVD group (Fig. 4B, E, H) was significantly weakened (the encoding phase, Fig. 4C, t = 14.220, P < 0.001; the maintenance phase, Fig. 4F, t = 7.765, P < 0.001); the retrieval phase, Fig. 4I, t = 13.030, P < 0.001). In addition, we compared the FC strength in theta band among three phases of WM task, both in the HC group (F = 0.28, df = 2, P = 0.752) and the CSVD group (F = 0.33, df = 2, P = 0.718) were no statistically difference. This indicates that it is the disease group rather than the phase of the WM task that affects the functional connectivity strength.
Functional connectivity distribution during the WM tasks
We use the topographic maps of DTF distribution to evaluate the spatial distribution of FC during the WM tasks. Consistent with the results of the previous study on neural oscillations, both in the two groups, the spatial distribution of DTF was mainly concentrated in the frontal region (Fig. 5A and B). The strength of FC in the CSVD group was significantly lower in the frontal region than the HC group (Figure 5C, t = 5.965, P <0.001). In addition, we further analyzed the FC of frontal region outflow and frontal region inflow. Figures 6A and B respectively show the average DTF matrices of the two groups in the four brain regions. The FC distribution mainly includes the frontal-frontal, frontal-parietal, and frontal-occipital in the HC groups (Fig. 6A). In the CSVD group (Fig. 6B), there were low strength and disordered connectivity among the four brain regions. From Fig. 6C, the FC strength of frontal outflow and the frontal inflow in the CSVD group were significantly lower than the HC group (t = 4.579, P <0.001; t = 5.402, P <0.001).
Correlation analysis
We calculated the Pearson correlation between DTF and behavioral indicators to further clarify the relationship between the strength of FC and behavioral performance in patients with CSVD. As shown in Figure 7A and B, there was a significant positive correlation between DTF and accuracy in both groups (CSVD: p <0.001, r = 0.966, HC༚p <0.001, r = 0.773). There was a significant negative correlation between DTF and reaction time in both groups (CSVD༚p <0.001, r = -0.763, HC༚p <0.001, r = -0.963).