All subjects were predominantly Russian from Moscow and the Moscow region. The enrolled cohort included 56 HTT mutation carriers (27 men, 29 women), a mean age 32.5 ± 1.1 years with a CAG repeat length of 43.9 ± 0.5 and 41 age-matched healthy subjects (15 men, 26 women), mean age of 32.8 ± 1.7 years. The group HTT mutations carriers comprised 45 individuals with PreHD (19 men, 26 women), a mean age of 31.4 ± 1.1 years with a CAG repeat length of 43.2 ± 0.5; and 11 individuals with EMHD (8 men and 3 women), a mean age of 36.9 ± 2.9 years with a CAG repeat length of 46.8 ± 1.5. The disease burden score (DBS) giving a measure of toxic load was calculated for each pre-HD participant using the formula: (CAG repeat length — 35.5) × age (Tabrizi et al., 2012).

“Preclinical” in this article denotes the absence of unequivocal motor signs as assessed with the Unified Huntington’s Disease Rating Scale (UHDRS) (Huntington Study Group, 1996), established during a subject’s most recent visit to the outpatient clinic of the Research Center of Neurology (Moscow). That visit took place no longer than one week before the EEG and neuropsychological examinations. All participants earlier underwent predictive DNA testing. Mutation carriers were defined as persons having more than 36 CAG repeats in one allele of the HTT gene. PreHD was defined as individuals with HTT CAG repeat expansions >36 and a the diagnostic confidence level (DCL) < 2 of the Total Motor Score (TMS) as part of the UHDRS, whereas EMHD was defined as individuals with an HTT CAG repeat expansion >36 and a DCL = 4. All examined EMHD patients had >40 HTT CAG repeat.

All patients with EMHD were at stage 2 of the Huntington’s Disease Integrated Staging System (HD-ISS) (Tabrizi et al., 2022). Among the 45 PreHD individuals who underwent EEG examinations, 40 subjects were at stage ≤1 of the HD-ISS. Five individuals with PreHD did not meet the HD-ISS criteria due to having a CAG expansion in the range of 37–39 (three subjects had 39, one had 38, and one had 37 CAG repeats in the HTT gene). Among the 19 PreHD individuals who underwent both fMRI and EEG examinations, 17 were at stage ≤1 of the HD-ISS. Two HD patient relatives with PreHD did not meet the HD-ISS criteria due to having 39 CAG repeats. Control analysis, excluding HTT mutation carriers with 37–39 CAG repeats, confirmed all the received results.

Exclusion criteria were neurological and psychiatric comorbidity, any medication, as well as the presence of clinical HD based on an unequivocal motor diagnosis on the UHDRS.

Healthy controls were individuals with no relatives having HD or other neurodegenerative diseases. Controls underwent neurological examination and cognitive screening, and were found to be free of dementia and other medical, psychiatric and neurological conditions, including cerebrovascular diseases, hypertension, epilepsy and endogenous disorders. Exclusion criteria were a history of neurological and psychiatric diseases and any kind of memory impairment. Pre-HD subjects and HC also underwent a neuropsychological battery that included the following tests: the Mini Mental State Examination (MMSE) (Folstein et al., 1975), the modified pictures memory test (van der Hiele et al., 2007), the Controlled Oral Word Association Test (FAS) (Benton and Hamsher, 1989), Luria memory words test (LMWT) (Luria, 1966),and the Spielberger state–trait anxiety inventory (STAI) (Spielberger, 1983), the Total Functional Capacity (TFC) from the UHDRS, Symbol Digit Modalities Test (SDMT) (Smith, 1982), Stroop Word Reading Test (Stroop, 1935).

In addition to other examinations, fMRI rsFC was examined in a subgroup of 29 HTT mutation carriers (14 men, 15 women), a mean age 35.2 ± 1.4 years with a CAG repeat length of 44.1 ± 0.7 and 25 age-matched healthy subjects (8 men, 17 women), a mean age 38.0 ± 1.9 years. The group HTT mutations carriers included 19 PreHD individuals (7 men, 12 women), a mean age 34.8 ± 1.4 years with CAG repeat length of 43.0 ± 0.6; and 10 individuals with EMHD (7 men and 3 women), a mean age of 36.0 ± 3.0 years with a CAG repeat length of 46.3 ± 1.5.

Structural images were acquired using a T1-weighted MPRAGE sequence: TR = 1,900 ms, TE = 2.47 ms; FOV = 256 × 256 mm2; flip angle = 10° slice thickness 1.0 mm; interslice distance 1 mm; number of slices = 176.

Functional scans were obtained at rest using T2_-weighted EPI sequence: TR = 1,500 ms, TE = 30 ms, flip angle 70°, slice thickness 2 mm, FOV = 190 mm, FoV phase 100.0%. The subjects were instructed to relax as much as possible, to lie quietly with their eyes closed (to exclude stimulation of the visual system) and not to think about anything in particular.

For analysis of rsFC signals in fMRI images we used CONN, which is a MATLAB-based open source toolbox (Functional Connectivity SPM Toolbox 2017, McGovern Institute for Brain Research, Massachusetts Institute of Technology) (Whitfield-Gabrieli and Nieto-Castanon, 2012).

CONN toolbox version 18b in conjunction with the SPM 12 software package (Wellcome Department of Cognitive Neurology, London, United Kingdom, www.fil.ion.ucl.ac.uk/spm) was used to perform all preprocessing steps. Functional images were slice-time corrected, realigned (motion corrected), and coregistered to their respective T1-weighted anatomical image. Images were then normalized to the Montreal Neurological Institute (MNI) standard space and spatially smoothed with an 8-mm Gaussian filter.

Denoising methods was then applied to minimize the impact of artifactual sources of signal variability. This included band-pass filtering (0.01–0.1 Hz), scrubbing (volumes showing displacement larger than the 97 × th percentile were censored), regressing out of the first 10 principal components (aCompCor) calculated within the maps of white matter (five components) and cerebrospinal fluid (five components) and regressing out of 24 head motion parameters, including linear and rotational indices, their temporal derivatives and their squared values.

The CONN-18.b toolbox was used to obtain a linear measure of functional connectivity based on bivariate correlation and bivariate regression coefficients between seed areas for ROI-to-ROI analysis (Whitfield-Gabrieli and Nieto-Castanon, 2012). ROIs of the whole brain were drawn from the template provided by CONN (conn/rois/atlas.nii). For the purpose of the analysis, BOLD signal time courses were converted to normally distributed scores with Fisher’s transformation, which allows for the use of second-level general linear Model analysis.

We performed a region of interest (ROI) analysis whereby ROIs were anatomically defined using the FSL Harvard-Oxford maximum probability cortical atlas, with bilateral regions divided into left and right hemispheres (164 ROIs).

Differences in demographic scores between the groups (HTT mutation carriers vs. HC; PreHD vs. HC; EMHD vs. HC; PreHD vs. EMHD) were tested using analysis of variance (ANOVA) for continuous variables (age, education), and the Mann–Whitney U test for categorical variables (sex). EEG parameters from each group were tested for normal distribution by the Wilk–Shapiro test, and in no cases were the data skewed.

The significance of the differences between the log-transformed EEG parameters was estimated using repeated measures of ANOVA in the general linear model (GLM), with group (HTT mutation carriers vs. HC; PreHD vs. HC; EMHD vs. HC; PreHD vs. EMHD as a between-subjects factor, and bands as a within-subject factor; age and gender were included in the analysis as covariates. Separate estimations using repeated measure ANOVA were performed on the EEG, firstly for the traditional broad bands, and subsequently for the 1-Hz subbands of theta and alpha (4–13 Hz). Post-hoc comparisons for between-subject effects and within-subject effects were analyzed using the Duncan test, and the level of significance was set to p < 0.05.

During fMRI rsFC analysis, our primary focus was on differences between two groups, HTT mutation carriers and HC, using t-tests, considering the smaller sample size compared to the EEG analysis. As the groups of the HTT mutations and HC had nonsignificant differences in gender and age, we did not include these factors in the analysis.

We compared ROI-to-ROI rsFCs between the HTT mutation carriers vs. HC using two-tailed t-tests. The resulting statistical maps were set with p < 0.05 false discovery rate (FDR) corrected.

In addition, we used one-sample t-test to assess rsFCs in the HC group. A similar analysis was conducted using one-sample t-test to evaluate rsFCs in the group of HTT mutation carries. In all analysis, significance was determined using two-tailed t-tests, with a significance level set at p < 0.05 FDR corrected. The analysis specifically was focused on rsFCs within circuits where a significant difference between HTT mutation carriers and HC had already been found.

To investigate associations between CAG repeat numbers in HTT gene, DBS, EEG characteristics and psychometric measures, Pearson correlation analysis was used whenever the data followed a normal distribution and Spearman rank correlations were calculated in other cases. The results were adjusted for multiple comparisons by applying FDR corrections at p < 0.05 using the Benjamini and Hochberg (1995). We conducted a correlation analysis, including only the relative power of the standard EEG bands (delta, theta, alpha. beta1, beta2) and the 1-Hz ranges (specifically 7–8 Hz, 8–9 Hz and the difference between 7–8 and 4–5 Hz), as these frequencies have shown significant differences in PreHD, EMHD and HC. An uncorrected significance level of p < 0.05 was considered to indicate a tendency.

Regression analyses between fMRI ROI-ROI connectivity strength and EEG alpha relative power were conducted directly in the Conn Toolbox using generalized PsychoPhyisiological Interaction analysis. The results were assessed utilizing an analysis with an FDR correction at a significance level of p < 0.05. This correction took into account the total number of pairwise correlations run in our ROI-ROI analysis (164×163/2), and computes an FDR correction accordingly. We examined only the interrelation between alpha relative power and fMRI rsFC, as the alpha rhythm has been previously identified as the most affected in HD.

The ANOVA results showed a significant interaction between band and groups (F[4,352] = 2.79, p = 0.026). Post-hoc comparison showed that the alpha relative power was decreased (p = 0.038) and the beta2 relative power was increased (p = 0.018) in HTT mutation carriers as compared to HC individuals (Figure 1A).

ANOVA revealed a robust effect of group, HTT mutation carriers vs. HC individuals, on the relative power of the 1 Hz spectral frequency bands in the range 4–13 Hz (F(Alper et al., 2006; Zhang et al., 2023) = 11.56, p = 0.001). The ANOVA results also showed a Group x Band interaction (F[8, 704] = 2,07, p=,037). In the group of individuals carrying HTT mutations, the relative power of the 7–8 Hz, 8–9 Hz and 9–10 Hz frequency bands was reduced compared to that in the HC individuals (p = 0.016, p = 0.008 and p = 0.017 for the 7–8 Hz, 8–9 Hz and 9–10 Hz frequency bands respectively) (Figure 1B).

We performed a further analysis of the differences in EEG spectral power separately for the PreHD, EMHD, and HC groups. A significant interaction effect was observed between the group and band factors on the relative power of traditional broad EEG bands (F[8,364] = 5.0 p=,00001) (Figure 2A). Post hoc comparison indicated that there was no significant difference between the PreHD and HC groups. Post hoc comparison showed that in the EMHD group compared to the HC group, the delta relative power was increased (p = 0.001), the alpha relative power was reduced (p = 0.00002) and the beta2 relative power was increased (p = 0.0004). In the EMHD group delta power was higher (p = 0.03), alpha power was lower (p = 0.0002), and beta2 power was higher (p = 0.01) than in the PreHD group. The effect size of the difference between the EMHD patients and the HC individuals was large for the delta relative power (Hedge’s g = 0.91), alpha relative power (Hedge’s g = 1.00), and beta2 relative power (Hedge’s g = 1.16).

ANOVA yielded significant effect of group on the relative power (preHD vs. healthy controls) of the 1 Hz spectral frequency bands in the range 4–13 Hz (F[2.91] = 18.32, p < 0.000001). A significant interaction effect was observed between the group and band factors (F[16.728] = 2.79, p = 0.0002). Post hoc comparisons showed that in the PreHD compared to the healthy controls, the relative power was significantly reduced in the 7–8 Hz (p = 0.03) and in the 8–9 Hz (p = 0.02) frequency bands (Figure 3). The effect size of the difference between the PreHD individuals and the HC individuals was significant in the 7–8 Hz (Hedge’s g = 0.76) and 8–9 Hz (Hedge’s g = 0.63) frequency ranges.

The patients with EMHD, in comparison to healthy controls, exhibited a significant reduction in the relative power within the 7–8 Hz (p = 0.01), 8–9 Hz (p = 0.0006), and 9–10 Hz (p = 0.0003), 10–11 Hz (p = 0.0002) (Figure 2B). The patients with EMHD, in comparison to PreHD individuals, showed a significant reduction in the relative power within the 8–9 Hz (p = 0.04), 9–10 Hz (p = 0.01), and 10–11 Hz (p = 0.003) ranges (Figure 2B). The effect size of the difference between the EMHD patients and the HC individuals was substantial in the 7–8 Hz (Hedge’s g = 1.22), 8–9 Hz (Hedge’s g = 1.67), 9–10 Hz (Hedge’s g = 1.105), and the 10–11 Hz (Hedge’s g = 0.85) frequency ranges.

An analysis of the correlation between the number of CAG repeats in the HTT gene and EEG parameters revealed a positive correlation between the number of CAG repeats and the delta and theta relative power (r = 0.46, p = 0.001; p-FDR = 0.005; r = 0.41, p = 0.002, p-FDR = 0.008 for delta and theta, respectively). The alpha relative power was inversely correlated with the number of CAG repeats (r = −0.43, p = 0.001, p-FDR = 0.005). The value for the correlation of DBS with the alpha relative power was r = −0.28, p = 0.04, p-FDR > 0.05.

The number of CAG repeats was inversely correlated with the relative power of the alpha subband at 8–9 Hz (r = −0.36, p = 0.009, p-FDR = 0.021) The relative power of the 7–8 Hz subband was not associated with the number of CAG repeats, but the difference in relative power of the 7–8 Hz and 4–5 Hz frequency subbands (7–8 Hz% minus 4–5 Hz%) showed significant negative correlation with the number of CAG repeats (r = −0.54, p = 0.00004, p-FDR = 0.0006). A significant negative correlation was also revealed between the DBS and the relative power in 8–9 Hz range (r = −0.38, p = 0.005, p-FDR = 0.013) as well as the difference in the relative power of the 7–8 Hz and 4–5 Hz frequency subbands (7–8 Hz% minus 4–5 Hz%) (r = −0.41, p = 0.003, p-FDR = 0.01).

We found significant correlation between the number of CAG repeats in HTT and age: r = −0.34, p = 0.012. To control for the potential influence of age on the correlations between CAG and EEG, we performed partial correlation analysis. The results have shown a significant partial correlation between CAG repeats and the relative power of the alpha and delta bands (r = −0.4, p = 0.003 and r = 0.41, p = 0.002 for alpha and delta respectively), as well as with the power of 1-Hz range 8–9 Hz (r = −0.42, p = 0.001), and the difference in relative power between the 7–8 Hz and 4–5 Hz subbands (r = −0.54, p = 0.0001.

In HTT mutation carriers, word fluency (FAS) scores showed inverse correlation with the delta (r = −0.35, p = 0.01) relative power and positive correlation with the alpha relative power (r = 0.37, p = 0.007). Analysis of correlation between spectral power of 1-Hz EEG frequency bands and psychometric characteristics revealed a tendency toward positive correlation between verbal fluency and relative power of the alpha subband at 8–9 Hz (r = 0.31, p = 0.025), and a tendency toward negative correlation between verbal fluency and relative power of the theta subband 4–5 Hz (r = −0.30, p = 0.03). The relative power of the 7–8 Hz subband correlated with total correct on memory working test (r = 0.44, p = 0.008). However, none of the correlations between EEG characteristics and psychometric scores survived FDR correction for multiple comparisons. No significant correlation was observed between relative power of any of the examined frequency bands and neuropsychological characteristics in the healthy control group.

In addition to all other examinations, fMRI rsFC was assessed in a subgroup of 29 HTT mutation carriers with a CAG repeat length of 44.1 ± 0.7 and 25 age-matched healthy subjects. Their demographic and psychometric characteristics of the participants who underwent fMRI evaluation were similar to those from the main groups (Table 2).

The ROI-to-ROI significant fMRI rsFC differences between the HTT mutation carriers and healthy controls are presented in Figure 3 and in Table 3. The positive rsFC in the following circuits of the frontal cortex was higher in HC individuals than in the HTT mutation carriers: the left superior frontal gyrus (SFG l) – the right rostral prefrontal cortex (RPFC r, salience network); the left superior frontal gyrus (SFG l) the left prefrontal cortex (LPFC l); the left superior frontal gyrus (SFG l) the anterior cingulate cortex (ACC, salience network) and the right frontal pole (FPr) –the left superior temporal gyrus (pSTG l). The positive rsFC was also higher in several fronto-cerebellar circuits in HC individuals than in the HTT mutation carriers.

The healthy controls had stronger negative rsFC than the HTT mutation carriers in the following interhemispheric circuits: the right hippocampus (Hippocampus r)– the left superior lateral occipital cortex (sLOC l); the right superior parietal lobule – the left posterior temporal fusiform cortex (pTFusC l); and the right superior lateral occipital cortex (sLOC r) – the left lingual gyrus (LG l). The negative rsFC was also higher in HC individuals than in the HTT mutation carriers in several cortico-cerebellar circuits, between the cerebellum and the following posterior cortical regions: sLOC r, SPL r, right posterior middle temporal gyrus (pMTG r).

The HTT mutation carriers also had stronger negative rsFC than the healthy controls between the right putamen and the contralateral frontal cortex (MidFG l).

The effect size of the difference in fMRI rsFC between the HTT mutation carriers and the HC was substantial, exceeding 0.8 (Table 3).

The PreHC individuals had higher positive rsFC than the patients with EMHD in the following circuits: the left rostral prefrontal cortex (RPFC l, salience network) – the left posterior superior temporal gyrus (pSTG l, language network (beta = 0.26, T = 3.63, p-FDR = 0.048); the left middle temporal gyrus (pMTG l) and the left rostral prefrontal cortex (RPFC l) (beta = 0.29, T = 3.8, p-FDR = 0.034).

The positive rsFC between the left postcentral gyrus (PostCG l) and left sensorimotor lateral cortex was lower in the PreHD individuals than in the patients with EMHD ((beta = −0.38, T = −4.67, p-FDR = 0.012).

The PreHC individuals had stronger negative rsFC than the patients with EMHD in the following circuits: the left lingual gyrus (LG l) and the left posterior superior temporal gyrus (pSTG l) (beta = −0.24, T = -4.31, p-FDR = 0.03). Negative rsFC was also higher in PreHD individuals than in EMHD patients in several cortico-cerebellar circuits: the left cerebellum 45 (cereb45 l) – the left posterior middle temporal gyrus (pMTG l) (beta = −0.31, T = -5.23, p-FDR = 0.0027); the vermis45 – the left posterior temporal gyrus (pSTG l) (beta = −0.28, T = -4.87, p-FDR = 0.007); the vermis45 – left angular gyrus (AG l)) (beta = −0.26, T = -3.87, p-FDR = 0.034); the vermis45 – left supramarginal gyrus (SMG l)) (beta = −0.24, T = -3.78, p-FDR = 0.034); the vermis45 – left middle temporal gyrus (pMTG l)) (beta = −0.24, T = -3.68, p-FDR = 0.034).

Regression analysis between EEG alpha relative power and fMRI was conducted on HTT mutation carriers and HC subjects who underwent both examinations. The subgroup included 28 HTT mutation carrier (14 men, 14 women), a mean age 35.3 ± 1.4 years with a CAG repeat length of 44.2 ± 0.7 and 17 HC subjects (4 men, 13 women), a mean age 39.9 ± 2.2 years. The group HTT mutations carriers included 18 PreHD individuals (7 men, 11 women), a mean age 34.9 ± 1.5 years with CAG repeat length of 43.0 ± 0.6; and 10 individuals with EMHD (7 men and 3 women), a mean age of 36.0 ± 3.0 years with a CAG repeat length of 46.3 ± 1.5.

In the overall sample that included both HTT mutation carriers and healthy control participants, reduced alpha relative power was associated with increased connectivity in several brain circuits (Table 4, Figure 4). The rsFCs of the frontal cortex included the MPFC.DMN-related networks and the putamen-frontal network. The occipital cortex circuits comprised rsFC between the right inferior lateral occipital cortex (iLOC r) and the right temporo-occipital fusiform cortex (TOFusC r). The networks of the temporal cortex included the rsFC between the right superior temporal cortex (pSTG r) and the right parahippocampal cortex (PaHC), as well as the interhemispheric connections involving the left hippocampus with the right posterior middle temporal gyrus (pMTG) and the right posterior inferior temporal gyrus (pITG). Additionally, they involved connectivity between the right pITG and the PaHC l and cerebellum.

We found that in the group of HTT mutation carriers, including individuals with PreHD and EMHD, there was a significant decrease in the relative power of both the alpha band and narrow-band activity specifically within the range of 7–10 Hz, when compared to the HC group. In the comparison between individuals with PreHD and HC, the only significant alteration observed was a decrease in the relative power of EEG within a narrow frequency range at the theta-alpha border and slow alpha activity (7–9 Hz), with Hedge’s g effect sizes of 0.76 and 0.63 for the 7–8 Hz and 8–9 Hz, respectively. These effect sizes, while not extremely large, still indicate a notable and statistically meaningful distinction between the two groups within these specific frequency bands. These results suggest that these frequency ranges may be particularly relevant for understanding the differences in neural activity between individuals with PreHD and those without the condition.

In the EMHD patients, apart from the more pronounced decrease in power within the7-10 Hz, a significant reduction in power within the classic alpha band was also detected in the present study. The difference in alpha relative power between EMHD and HC individuals was substantial, with an effect size of 1.0. Previous studies have consistently reported a reduction in alpha power in individuals with HD (de Tommaso et al., 2003; Bellotti et al., 2004; van der Hiele et al., 2007; Hunter et al., 2010; Painold et al., 2010, 2011; Leuchter et al., 2017; Piano et al., 2017; Odish et al., 2018; Delussi et al., 2020). The reduction in alpha power observed in the EMHD group may reflect a dysregulation of cortical excitability due to the alteration of inhibitory GABAergic mechanisms and may also be associated with excitotoxicity (Palva and Palva, 2007).

Previous studies have reported a decrease in power at the theta-alpha border in individuals with preHD and EMHD (Painold et al., 2011; Ponomareva et al., 2014; Hawellek et al., 2022).

This frequency range, 7–9 Hz, has been associated with thalamo-cortical mechanisms involved in cognitive functions. A decrease in 7–8 Hz activity has been observed during the retrieval of semantic memories (Slotnick et al., 2002), whereas this activity typically increases during a state of deep relaxation and reverie (Sokhadze et al., 2008). These results suggest a potential link between the decline in 7–8 Hz power and cortical disinhibition in preHD and HD patients.

Furthermore, a decrease in power within this frequency range has also been observed during sleep in preHD and EMHD patients (Lazar et al., 2015). Researchers have suggested that this decrease may underlie deficits in neural restoration following wakefulness. The hypothesis linking 7–8 Hz activity to the processes of HD development was supported by evidence that the use of the huntingtin-lowering oligonucleotide tominersen led to a partial recovery of activity within this frequency range (Hawellek et al., 2022).

In the EMHD group, we found an increase in beta2 activity. Previous EEG studies in HD, including the research by Bylsma et al. (1994) identified an increase in beta activity, while contrasting results were found in other studies (Painold et al., 2011; Piano et al., 2017). The discrepancies in the alteration of the beta band may be related to the heterogeneity of HD among patients and variations in the EEG analysis methods used. In our study, as well as in the research of Bylsma et al. (1994), normalization was applied to the analysis of EEG spectral power. Normalization minimizes interindividual variability in absolute power, which arises from factors such as skull thickness and impedance, and allows for a better comparison of an individual’s power distribution across frequency bands. However, a reduction in the spectral power of the alpha band found in EMHD can result in an increase in the relative spectral power of other frequency ranges, especially in the beta bands, which typically have lower power. Results obtained from the R6/2 mouse model of HD show an increase in beta activity during the disease progression (Fisher et al., 2013, 2016). Further research is needed to clarify whether an increase in beta activity may indeed occur in some patients during specific stages of HD progression.

We observed a correlation between an increase in delta power and a decrease in alpha power, and poorer performance on cognitive tests, particularly verbal fluency. We also revealed a positive correlation between the relative power of the 7–8 Hz sub-band relative power and performance on memory working test. The correlations between EEG characteristics and psychometric scores did not survive FDR correction for multiple comparisons and should be considered as tendencies.

These findings align with prior studies that have shown a connection between a reduction in alpha and increase in delta power and cognitive decline in individuals with HD (Bylsma et al., 1994; Painold et al., 2010, 2011; Ponomareva et al., 2014; Odish et al., 2018). Changes in spectral power in individuals with HD have been found to correlate with performance on cognitive tests, including the Mini Mental Status Examination, FAS verbal fluency test, Stroop Word Reading Test, and the Symbol Digit Modalities Test.

In the present study a positive correlation was found between the length of the CAG repeats in the HTT gene and the relative power of delta and theta frequencies in HTT mutation carriers. On the other hand, the CAG repeat length was negatively correlated with the alpha band and its subbands. Additionally, a negative correlation was found between the difference in relative power of the 7–8 Hz and 4–5 Hz and the number of CAG repeats. These results validate our prior findings on a larger sample size (Ponomareva et al., 2014). Previous research has also revealed a correlation between EEG characteristics, such as the relative delta anterior–posterior gradient, and the number of CAG repeats (Nguyen et al., 2010, for review, Hunter et al., 2010; Piano et al., 2017, for review).

Numerous studies of PreHD and EMHD have employed data-driven approaches to explore whole-brain connectivity and these studies have consistently demonstrated altered functional connectivity patterns (Dumas et al., 2013; Harrington et al., 2015; Poudel et al., 2015; Sánchez-Castañeda et al., 2017; Espinoza et al., 2018; Pini et al., 2020; Polosecki et al., 2020). The results from these studies have shown a range of differences, which can be attributed to the variation in the cohorts studied and the methodologies used in the investigations.

Evidence from fMRI studies indicates the presence of distributed anticorrelated networks in the brain (Fox et al., 2005; Martinez-Gutierrez et al., 2022). Task-positive networks, such as the frontoparietal and salience networks, are involved in cognitive control and attention modulation, becoming more active during external goal-directed tasks. On the other hand, the DMN, crucial for memory function, shows reduced activity during external goal-directed tasks (Sheline and Raichle, 2013). In healthy adults, rsFC between task-positive and task-negative networks is typically characterized by negative correlation.

The current study identified significant changes within cortico-cortical, putamen-cortical, cortico-hippocampal, and cortico-cerebellar networks in individuals carrying HTT mutations, when compared to healthy controls. Altered connectivity was seen in networks that show both negative and positive connectivity in healthy individuals. The effect size for the difference in fMRI rsFC in these circuits between HTT mutation carriers and the HC was notably large, exceeding 0.8.

HD is primarily characterized by the progressive degeneration of the striatum, with the putamen being one of the key regions affected. This neurodegeneration is largely attributed to the loss of GABAergic medium spiny neurons (Spampanato et al., 2008). We found a significant difference in rsFC between the right putamen and the left middle frontal gyrus when comparing HTT mutation carriers to healthy controls. In HTT mutation carriers this circuit displayed negative rsFC, while connectivity values in HC individuals were not significantly different from zero. Previous studies demonstrated reduced rsFC in fronto-striatal circuits in preHD and HD patients (Poudel et al., 2015; Espinoza et al., 2018; Nair et al., 2022). However, Kronenburger et al. (2019) found increased striatum-prefrontal connectivity in premanifest HD subjects. The authors suggested that enhanced rsFC may imply a compensatory mechanism, where additional cortical regions are recruited to subserve functions that have been impaired due to HD pathology.

Neurodegeneration in the frontal lobe is a prominent feature of HD (Aylward et al., 1998; Matsui et al., 2014). Prior studies have reported aberrant rsFC patterns in the frontal lobe of HD patients (Thiruvady et al., 2007; Harrington et al., 2015; Langley et al., 2021). Reduced rsFC was found within the bilateral superior frontal gyrus between ventral attentional center and the frontal area (Zhang et al., 2023, for review). However, other studies have found an increase in frontal resting-state functional connectivity (rsFC) and enhanced frontal-posterior rsFC (Werner et al., 2014; Harrington et al., 2015). The current study showed that in HTT mutation carriers, compared to healthy controls, a significant decrease in positive rsFC is observed in several cortical networks, particularly within the frontal lobes. Specifically, decreased rsFC was found between the left superior frontal lobe and the regions within salience networks.

We have revealed aberrant cortico-cerebellar rsFC in HTT mutation carriers compared to HC individuals. In accordance with prior research, the findings of the present study identified aberrant connectivity between the cerebellum and the cortical regions implicated in executive functions, including the prefrontal and parietal cortices (Wolf et al., 2015; Sarappa et al., 2017; Franklin et al., 2021). These findings align with the disturbance of executive function in HD, supporting the notion of involvement of the cerebellum in the nonmotor symptoms of the disease.

In the current study, we found diminished negative rsFC in individuals with preHD compared to the healthy controls, particularly in interhemispheric networks. Specifically, we observed reduced rsFC between the right superior parietal lobule (SPL) and the left posterior fusiform cortex (pTFusC), between the right lateral occipital cortex (sLOCr) and the left lingual gyrus (LGl), and between the right hippocampus and the left sLOC. Previous studies have suggested the importance of abnormal interhemispheric connectivity in the pathophysiology of HD (Harrington et al., 2015).

Our findings indicate a prominent prevalence of neural networks showing an inverse correlation between fMRI connectivity patterns and EEG alpha relative power in HTT mutation carriers and HC individuals. This finding aligns with the results of previous studies that have reported a negative association between the functional connectivity of cortical circuits and alpha power (Scheeringa et al., 2012; Timmermann et al., 2023).

In particular, we found a significant relationship between alpha power and the rsFC of circuit located in frontal regions and the putamen. Studies show that the putamen is involved in modulating alpha rhythm (Hughes and Crunelli, 2005). Early dysfunction and degeneration of the putamen in HD may influence the alpha rhythm and may underlie the association of fMRI rsFC and alpha power. Our results are consistent with studies suggesting that the impact on the cortico-striato-thalamo-cortical circuitry is an important factor in the abnormal reduction of intrinsic alpha oscillations observed throughout the cortex in HD (Hughes and Crunelli, 2005; Alper et al., 2006). Furthermore, the degeneration of GABAergic projections from the putamen to the frontal cortex can impair inhibitory function and contribute to the dysregulation of cortical rhythms (Smith et al., 1998).

The present study’s results indicated that the rsFC between the MPFC, a central hub within the DMN, and the SFG r was associated with a decrease in alpha rhythm. Previous studies have reported that the MPFC and the SFGr are involved in the pathophysiology of HD (Wolf et al., 2012; Pini et al., 2020). Accumulating evidence suggests that the MPFC and the SFG r work together as part of a cognitive control network (Mathewson et al., 2014). The MPFC provides contextual information and top-down signals, while the SFG r exerts cognitive control to suppress resting alpha rhythms across multiple cortical areas.

We also found significant association between a decrease of EEG alpha power and an increase of rsFC between the right inferior lateral occipital cortex (iLOC) and the right temporo-occipital fusiform cortex (TOFusC). It was consistently shown that lower alpha power indicates higher cortical excitability and lower perceptual threshold (Foxe and Snyder, 2011). Our results suggest that pathological process in HD may disrupt normal communication between visual areas, and that the reduction in alpha power could indicate compensatory increases in connectivity in response to disease-related changes.

Additionally, we found an association of the alpha power with the interhemispheric circuits between the left hippocampus with the pMTG r and the pITG r, as well as between the right pITG, the left PaHC and the cerebellum. Prior research has shown that hippocampal dysfunction worsens in a linear manner during the development of HD. Furthermore, studies have suggested that abnormalities in the hippocampus could serve as a marker for ongoing cognitive decline in individuals with HD (Harris et al., 2019). Our results indicate that the alteration of hippocampal networks may be a factor contributing to the reduction in alpha power in HD.

The precise mechanisms underlying the alteration of the alpha rhythm in HD are still not fully understood. The findings of the current study indicate that alteration in several cortico-cortical networks, including the putamen-frontal, DMN-related and cortico-hippocampal circuits, can contribute to the reduction in alpha power in HD.

The limitations of the present study include the relatively small number of participants, who underwent fMRI examination, especially in the EMHD group. This suggests that the results on the association of EEG and rsFC in HD are preliminary. Nevertheless, the findings on the interdependence of the EEG and fMRI characteristics can contribute to an understanding of the underlying mechanisms behind alpha rhythm abnormalities in HTT mutation carriers. Another limitation is the nonsimultaneous registration of EEG and fMRI. However, as in the resting state of these characteristics are relatively stable, separate EEG and fMRI recordings across experimental sessions are sufficient (Scrivener, 2021).

In conclusion, the study identified significant differences in EEG power between individuals with PreHD and healthy controls, with a notable decrease in power within a specific frequency range at the theta-alpha border and slow alpha activity (7–9 Hz). In EMHD patients, a reduction in classical alpha band power was observed compared to HC, while a decrease in relative power within the 7–9 Hz frequency range effectively distinguishes preHD individuals from healthy controls.

Furthermore, fMRI analysis revealed alterations in functional connectivity within various brain networks, particularly within the frontal lobe, putamen-cortical, and cortico-cerebellar networks, in individuals carrying the HTT mutation compared to healthy controls.

The analysis of the relationship between EEG and fMRI rsFC revealed a correlation between a decrease in alpha power, observed in individuals with EMHD, and a reduction in inhibitory control in cortico-cortical networks. Such networks include putamen-cortical, DMN-related and cortico-hippocampal circuits. These alterations in functional connectivity may contribute to the observed EEG changes in HD.

Overall, the findings suggest that EEG and fMRI provide valuable information about pathophysiological mechanisms of HD progression during preclinical and early manifest stages of the disease. A decrease in inhibitory control within the putamen-cortical, DMN-related and cortico-hippocampal circuits, accompanied by a reduction in alpha and theta-alpha border oscillatory activity, could potentially contribute to cognitive decline in HD.