Introduction

Alcohol use disorders (AUDs) are associated with adverse physical, psychological and social consequences, leading to 3.3 million deaths each year1. While AUD patients’ profile of cognitive impairment2,3,4,5 has been associated with the specific susceptibility of frontal cortex6, neuroimaging studies highlighted an extensive damage fitting with the diffuse brain hypothesis7. Voxel-based morphometry (VBM) studies have indeed shown, in AUDs, an atrophic pattern involving lateral prefrontal cortex, anterior (ACC) and posterior cingulate cortex, insular-opercular cortex, thalamus, hippocampus and striatum8,9. This evidence supports the hypothesis that the cognitive alterations observed in addictions involve a cortico-striatal-thalamic circuit including the key-nodes of the executive and salience networks10, with the latter mediating the switch between rest and controlled cognition and behaviour10,11,12.

Since abstinence appears to reverse behavioural and neuro-structural damage13, treatment protocols would benefit from a detailed characterization of the neural correlates of cognitive impairment(s). Previous studies have addressed a relationship between atrophy and performance in different domains, e.g. executive functioning14,15, cognitive control16,17 and memory18. However, their inconsistent results14,15, and the lack of comprehensive neuropsychological assessments, highlight the need of further evidence on the neural bases of cognitive impairment(s) in AUDs.

We pursued this goal at the neuro-structural level, by coupling VBM with an extensive neuro-cognitive assessment. Unlike previous studies we adopted a multivariate approach to investigate, in AUDs, defective cognitive domains transcending specific tasks. We expected impaired performance to be predicted by GM density in regions within the executive or salience networks.

Results

Cognitive impairments in alcoholic patients

There was no significant group difference concerning age, education or nicotine consumption (Table 1). Compared with controls, patients displayed significantly worse performance in the ENB global score, and in different sub-scores: immediate recall, interference memory-10”, TMT-A and B, overlapping pictures and clock-drawing (Table 2a). TMT performance is indexed by the total time to completion19, i.e. shorter time corresponds to better performance. Age was negatively related to performance in TMT-A (i.e. positively correlated with response time; r = 0.46, p = 0.003), and both immediate and delayed recall (both r = −0.31, p = 0.046). Except for the clock-drawing task, even when controlling for age or education via ANCOVA we confirmed, in patients vs. controls, worse performance in all these tasks (Table 2b,c). A correction for multiple comparisons based on False Discovery Rate (FDR) confirmed the latter evidence, with the only exception of the immediate recall score which showed a marginal trend (p-corrected = 0.058).

Table 1 Demographics and alcohol use variables.
Table 2 Neuro-cognitive performance.

This evidence was refined by a multivariate analysis of neuro-cognitive data. Based on the Kaiser-Guttman criterion, considering components with eigenvalues > 120, the initial dataset of 15 ENB2 scores was reduced to 6 components explaining 74.89% of the variance (Supplementary Table S1). These components involve different domains such as visual-constructional abilities, verbal learning, basic-level and high-level executive processes, and language (Table 3). The sixth component, associated with the digit span and cognitive estimation tasks, might reflect the efficiency of processes related to cognitive estimation21. A strongly significant group difference was found in the factor score of the third component (F(1,39) = 11.58, p = 0.002). The latter (henceforth “basic-level executive component”) reflects attentional (TMT-A) and working-memory (interference-memory tasks) performance. TMT-A provides the strongest contribution to this component (r = −0.78, p < 0.001), followed by interference-memory-10” (r = 0.71, p < 0.001) and interference-memory-30” (r = 0.69, p < 0.001). Marginal evidence for impaired performance in patients (F(1,39) = 3.74, p = 0.060) was also found in the fourth component (“high-level executive component”), associated with TMT-B, overlapping pictures and abstraction tasks. Therefore, only the basic-level executive component and TMT-A response time (its main contributing variable) were considered in subsequent analyses.

Table 3 Principal component analysis of neuro-cognitive data.

Grey matter atrophy in alcoholic patients

VBM results highlighted a distributed pattern of GM density reduction in alcoholic patients vs. controls (Table 4a; Fig. 1A; Supplementary Table S2). GM atrophy involved both the dorsal (superior medial gyrus and dorsal sector of anterior cingulate cortex (dACC)) and ventral (rectus gyrus) sectors of the medial prefrontal cortex, alongside the pars triangularis of the right inferior frontal gyrus. GM density was significantly reduced, in patients, also in the rolandic operculum and posterior insular cortex, bilaterally but with a right-sided prevalence, as well as in the posterior temporal cortex (superior and middle temporal gyri) bilaterally. Also the left postcentral gyrus and the medial parietal cortex (precuneus and posterior cingulate cortex) displayed a significant GM atrophy, alongside the ventral striatum, caudal thalamus and left amygdala.

Table 4 Neuro-structural correlates of executive impairment in AUDs.
Figure 1
figure 1

Grey-matter density reduction in AUDs. The brain regions in which grey matter density was (A) significantly reduced in AUD patients vs. controls; (B) positively correlated with executive performance; (C) negatively correlated with TMT-A response time (p < 0.025 corrected for multiple comparisons).

Correlation between grey matter density and executive performance

We first assessed a relationship between voxel-wise GM density and the degree of cognitive impairment regardless of group. We found a positive correlation between the scores of the impaired executive component and GM density in fronto-insular and fronto-basal structures, alongside the cerebellum (Table 4b; Fig. 1B; Supplementary Table S3). While the anterior insula was involved bilaterally, in the left hemisphere we observed a wider cluster extending from the lateral orbital gyrus to the posterior insula. Executive performance was also related to GM density in the left amygdala and ventral striatum, in a right-hemispheric posterior insular cluster extending from the rolandic operculum to the superior temporal and supramarginal gyri, as well as in the left cerebellum.

Group differences in the correlation between grey matter density and executive performance (interaction analysis)

We found no significant interactive effect of group on the relationship between GM density and executive performance.

Correlation between grey matter density and TMT-A response time

Across both groups, TMT-A response time was negatively correlated with GM density in several brain structures (Table 4c; Fig. 1C; Supplementary Table S4), i.e. dACC alongside anterior insula (extending into inferior frontal cortex) and posterior insula (extending into the supramarginal gyrus) bilaterally. An analogous relationship was also observed in the right fusiform gyrus, precuneus, occipital cortex and right amygdala.

Group differences in the correlation between grey matter density and TMT-A response time (interaction analysis)

We found no significant interactive effect of group on the relationship between GM density and TMT-A response time.

Common neuro-structural effects of AUDs and correlation with executive performance

A conjunction-analysis highlighted the spatial overlap between the voxels associated with the neuro-structural and cognitive impairments observed in alcoholic patients (Figs 2 and 3). This analysis unveiled a clear correspondence between the regions where GM density is a) reduced in patients vs. controls, and b) related to executive performance (Table 4d; Figs 2A and 3; Supplementary Table S5). Such overlap was found bilaterally in a cluster encompassing the anterior-middle insular cortex, as well as in the medial temporal cortex and ventral striatum. We found further common effects in the left inferior frontal cortex (from pars opercularis to pars orbitalis) and right posterior insular cortex (from the rolandic operculum to the superior temporal and supramarginal gyri).

Figure 2
figure 2

Correlation between grey-matter density and executive performance in AUDs. The brain regions in which grey matter density was (A) both significantly reduced in AUD patients vs. controls and positively correlated with executive performance; (B) both significantly reduced in AUD patients vs. controls and negatively correlated with TMT-A response time (p < 0.025 corrected for multiple comparisons). The scatterplots in panel (A) additionally depict the relationship between executive performance and average grey matter density in the left and right fronto-insular cortex, either in healthy controls (HC), alcoholic patients (PT) or both.

Figure 3
figure 3

Common neuro-structural effects of AUDs and correlation with executive performance. The brain regions showing specific vs. common effects of AUDs, executive performance or TMT-A response time (p < 0.025 corrected for multiple comparisons).

We then extended this approach to the TMT-A task, to investigate a significant overlap between the regions in which GM density was both related to response time and reduced in patients vs. controls (Table 4e; Figs 2B and 3; Supplementary Table S6). Such overlap was found in the dorsal and ventral sectors of the medial prefrontal cortex (dACC and rectus gyrus, respectively), bilateral insula (extending rostrally into the left inferior frontal and caudally into the supramarginal and superior temporal gyri), as well as limbic and basal structures (particularly amygdala and ventral striatum).

Regions-of-Interest statistical analyses

Homogeneity-of-slopes models confirmed the lack of qualitative group differences in the relationship between GM density and either executive performance or TMT-A response time. The representative scatterplots depicted in Fig. 2A indeed show comparable slopes for the regression of executive performance on GM density in the left (F(1) = 0.28, p = 0.598) and right (F(1) = 0.52, p = 0.474) anterior insular clusters which also displayed a significant atrophy in patients.

When considering all the clusters showing GM atrophy in alcoholic patients (Fig. 1A and Table 4a), a strongly significant multiple regression model (F(1,38) = 11.32, p < 0.0005) showed that 34% of variability in the executive component is explained by average GM density in the dACC (β = 0.94, t(38) = 4.44, p < 0.0001, observed power = 0.99; adjusted r2 = 0.34). As expected, constraining such approach to the three clusters in which GM density was also related to executive performance (as highlighted by a conjunction-analysis; Fig. 2A and Table 4d) increased the explained amount of its variance (F(1,38) = 29.43, p < 0.00001; adjusted r2 = 0.42). In this case, the only retained predictor was GM density in the left anterior insular cluster.

Salience network structural integrity and executive impairments

Based both on our initial hypotheses, and on VBM evidence of structural impairment in the key nodes of the salience network (Figs 1 and 2), we aimed to assess the degree to which such damage accounts for the observed cognitive impairment. To this purpose, we assessed a relationship between executive performance and GM density in the voxels which, besides showing overlapping effects of interest (atrophy in patients, correlation with executive performance and TMT-A response time), were additionally included in the salience network (as highlighted by Neurosynth meta-analytic evidence; see 4.7). Such overlap was found in the anterior insula and frontal operculum bilaterally (left xyz: −38, 18, −10; right xyz: 45 12 −4), alongside right posterior insula (xyz: 50 −30 20) (Fig. 4). Average GM density in this set of regions was significantly correlated with the ENB2 global score (r = 0.33, p = 0.036), and particularly with executive performance (r = 0.64, p < 0.0001) (Fig. 4). The latter relationship appeared to be more strongly driven by TMT-A (r = −0.66, p < 0.0001) than by the other significantly related behavioural measures, i.e. working-memory (interference-memory-10”) and overlapping pictures (both r = 0.38, p = 0.014). Homogeneity-of-slopes models confirmed the lack of significant group differences in the relationship between average GM density in these regions and executive (F(1) = 0.26, p = 0.610) or TMT-A (F(1) = 0.29, p = 0.596) performance.

Figure 4
figure 4

Salience network and executive impairment in AUDs. The overlap between the salience network (red) and the brain regions showing common effects of AUDs and either executive performance (blue) or TMT-A response time (green). The scatterplots depict the significant relationship between average grey matter density in the overlapping voxels (white colour) and either executive performance (left) or TMT-A response time (right).

Overall, GM density in the salience network voxels highlighted by this analysis accounted for 43% and 40% of variability in, respectively, TMT-A and executive performance in the whole sample (p < 0.00001).

Discussion

We investigated the neural bases of cognitive impairment(s) in AUDs, by coupling a comprehensive neuro-cognitive assessment with neuro-structural VBM evidence of regional GM atrophy.

In patients, a global cognitive impairment was mainly driven by abnormal performance in tasks tapping working-memory (interference memory), visuomotor processing speed and attention (TMT-A), as well as divided attention, switching and mental flexibility (TMT-B; Overlapping pictures). The most defective tasks, i.e. TMT-A and interference-memory, clustered in a strongly impaired basic-level executive domain involving attention and working-memory. Within a global cognitive impairment22, our data thus confirm the prominence of an executive disorder in AUDs23. It is still debated, however, whether this deficit reflects a specific susceptibility of frontal regions (“frontal lobe” hypothesis6), or a diffuse pattern involving other cortical or subcortical structures (“diffuse brain hypothesis”7). We used VBM to distinguish between these hypotheses, by providing a detailed characterization of the neuro-structural bases of cognitive impairment in AUDs.

In line with previous data9, VBM results highlighted a diffuse pattern of GM decrease, in patients, along the subcortical (ventral striatum, thalamus, hippocampus and amygdala,) and cortical (ventromedial and posterior dorsomedial frontal cortex) components of the meso-cortico-limbic pathway (Fig. 1A). Based on functional magnetic resonance imaging (fMRI) evidence for its role in adaptive behavioural learning24,25, an impairment of this pathway has been suggested to underpin the development and maintenance of addiction, via negative reinforcement26. GM atrophy was also found in the dACC and insular cortex27, i.e. the main nodes of the so-called “salience network”10. This evidence of widespread brain damage fits with the notion that the cognitive and behavioural alterations observed in substance-use disorders, including AUDs, might reflect functional imbalances within a cortico-striato-thalamo-cortical regulatory loop underlying executive control and self-regulation10.

This circuit underpins another crucial component of behavioural learning28, i.e. regulating the switch between rest and effortful cognitive activity based on the salience of external stimuli with respect to behavioural goals. This mechanism is driven by the salience network, involving both cortical (dACC and anterior insula) and subcortical (basal ganglia, thalamus and amygdala) nodes of the loop10,12. Their functional interaction is central to detect salient (i.e., motivationally relevant) stimuli, and to enable them to modulate cognition and behaviour via executive processes of response inhibition and selection associated with the lateral prefrontal cortex10. This neural mechanism facilitates the access to attention and working-memory resources, by activating the dorsal attentional and fronto-parietal executive control networks, when relevant stimuli are detected12. The vicious circle linking executive and behavioural impairments might then be further reinforced, in AUDs, by abnormal processing of salient stimuli, due to the structural impairment of the key nodes of the salience network. To test this hypothesis, we assessed the spatial overlap between the regions displaying atrophy in patients and those in which GM density relates to impaired cognitive performance.

We first assessed the relationship between GM density and performance in the executive domain showing the greatest impairment in alcoholic patients. The lack of significant task-by-group interactions suggests that no qualitative difference exists, between patients and controls, in this relationship. Instead, correlational analyses highlighted evidence of quantitative group differences in this respect (Figs 1B and 2A): in the whole sample, executive performance was positively correlated with GM density in different clusters encompassing fronto-insular and medial temporal cortex bilaterally, alongside right posterior insular cortex, left cerebellum and ventral striatum. While showing a clear segregation between patients and controls, the scatterplots displayed in Fig. 2A indeed highlight, in the two groups, comparable slopes in the relationship between executive performance and GM density in the left and right anterior insula.

In line with behavioural evidence, TMT-A response time provided the strongest contribution in explaining the relationship between brain atrophy and executive dysfunction in AUDs. Moreover, the lack of a significant task-by-group interaction confirmed that such relationship reflects quantitative, rather than qualitative, differences along a continuum from normal to impaired conditions. In the case of TMT-A this relationship involves a general slowing of performance associated with grey matter loss in regions largely overlapping those previously described for executive performance, with the additional contribution of the dACC (Figs 1C4). The involvement of this structure supports the relationship between executive performance, particularly in the TMT-A, and GM density in the salience network10. Figure 4 (white colour) indeed shows the overlap between a meta-analytic map of this network and the bilateral anterior and right posterior insular clusters showing both atrophy in alcoholic patients and a significant correlation between GM density and executive performance. In the whole sample, average GM density in these voxels accounted for almost half of variability in TMT-A response time (adjusted r2 = 43%) and executive performance (adjusted r2 = 40%).

Overall, these data support an interpretation of AUD patients’ executive impairment in terms of defective access to attention and working-memory resources, due to the structural damage of a salience-detection neural mechanism in charge of activating the executive network12. In particular, the insular cortex is considered to link perceptual, cognitive and autonomic information, by relaying salience signals related to body states, generated by viscero-autonomic sensors and transmitted by thalamic nuclei29. An interoceptive awareness of salient stimuli involves the posterior insular sector of the network, which we also found to be both related to executive performance in our sample, and structurally impaired in patients. These data might thus reflect a domain-general role of the insular salience node in enhancing access to computational resources. Instead, the specific association between TMT-A response time and GM density in the dACC sector of the salience network (Fig. 4, yellow colour) may reflect the prominent role of this region in conflict monitoring and response selection30.

Within a diffuse pattern of GM atrophy, our data thus highlight an impairment of the salience network as a specific neuro-structural correlate of executive dysfunction in AUDs, likely related to previous evidence of defective functional activation of the executive network31. The present evidence of a more severe impairment of the insular and dACC nodes of the salience network, compared with the lateral prefrontal nodes of the executive network, fits with a recent meta-analysis of previous VBM studies on AUDs9. The prominent damage of insular and anterior cingulate cortex likely reflects the susceptibility of so-called von Economo neurons, localized in these regions32, to the neurotoxic effects of alcohol27. Importantly, there are limitations to interpreting the functional significance of neuro-structural data. However, the present results complement previous evidence at the functional level, i.e. abnormal connectivity between the insular and dACC nodes of the salience network highlighted by studies based both on fMRI33 and arterial spin labelling34.

A limitation of this study is represented by the lack of in-depth measures of executive functioning and cognitive control, such as response inhibition; in the trade-off between comprehensiveness and specificity, however, opting for a broad neuropsychological assessment allowed us to highlight the selective impairment of an executive domain transcending specific attentional and working-memory tasks. This novel evidence contributes to a lively debate between two opposite views rooted in the focal vs. global nature of the cognitive and neural dysfunctions in AUDs, and that our data seem to reconcile. The cognitive impairment in different tasks, associated with functional alterations in multiple brain structures, might indeed reflect the impairment of a specific executive domain, due to the structural damage of a well-established salience network activating the executive control network when relevant stimuli are detected. Another limitation is represented by a relatively small sample; by focusing on few cognitive variables showing the strongest impairment in patients, however, we identified a strongly significant relationship between its severity and the extent of GM atrophy (Fig. 4). Finally, we report only cross-sectional evidence which will require further support from future longitudinal studies.

Identifying a specific neural substrate of executive impairment in AUDs has several implications. First, we observed quantitative, rather than qualitative, group differences in the relationship between GM density and domain- or task-specific executive performance. This evidence might seem to challenge the notion of adaptive neural mechanisms supporting cognitive performance via compensatory brain regions35. Such neuroadaptation, however, might become detectable at the structural level only after longer abstinence periods. Second, the salience-detection mechanism discussed above might also influence decision-making, and thus abstinence vs. relapse, by underpinning the interplay between a “reflexive” limbic system driving automatic behaviours and a “reflective” fronto-striatal network activating executive control23. A defective switching mechanism between these systems may thus bias decision-making processes towards bottom-up impulsive signals, at the expenses of top-down goal-driven attentional resources required to exert behavioural control over alcohol search and consumption36,37. This hypothesis fits with previous morphometric evidence showing that AUD patients’ increased impulsivity correlates with the degree of grey matter atrophy in the anterior insular and cingulate nodes of the salience network11.

The present results highlight several directions for future research. Growing evidence shows that long-term abstinence can reverse both behavioural5,38 and neuro-structural13,14,39 alterations. Moreover, preliminary evidence highlights positive effects of cognitive remediation40 and neurostimulation41 on patients’ cognitive performance and craving. However, treatment efficacy depends on the integrity of executive functions42, and the presence of cognitive impairment requires adapting management strategies based on individual profiles43. Therefore, both the design of treatment protocols, and the assessment of their effectiveness, require sensitive benchmark metrics of executive functioning. Moreover, the efficacy of neurostimulation protocols depends on the choice of the target area, which in turn will reflect previous evidence of a significant relationship between its structural or functional properties and cognitive performance. By showing a restricted set of tasks which are particularly sensitive to AUD patients’ executive impairment, and its neural correlates, our results might thus help tailor remediation or neurostimulation treatment protocols to target specific brain networks and their associated cognitive functions.

Methods

See Supplementary methods for additional details.

Participants

Twenty-three alcoholic patients (9 females; mean age: 45.69 years ± 7.82) and 18 healthy control subjects (8 females; mean age: 44.83 years ± 8.86) participated in the study. Patients were interviewed to determine their drinking history, including the amount, type and lifetime duration of alcohol usage. Alcohol consumption was calculated as the average number of standard units of alcohol (UA) per day (one UA: 330 ml beer, 125 ml wine, or 40 ml hard liquor, corresponding to 12 g of ethanol) (see Table 1 for demographic, as well as nicotine and alcohol use, variables).

Inclusion criteria for patients were: (1) age between 20 and 60 years; (2) a diagnosis of alcohol dependence according to DSM-V criteria. Exclusion criteria for both groups were: (1) presence/history of neurological/psychiatric disorders other than AUDs, or any comorbid disorder except for nicotine dependence; (2) family history of neurological/psychiatric disorders; (3) current use of any psychotropic substance/medication; (4) past brain injury or loss of consciousness; (5) major medical disorders (e.g. kidney or liver diseases, sever diabetes and/or malnutrition); (6) inability to undergo the neuropsychological assessment; (7) contraindications to magnetic resonance imaging (MRI). Controls were excluded in case of presence/history of alcohol abuse. Patients joined the protocol after being detoxified for at least 10 days, via medically supported standard treatments. However, they had ceased benzodiazepine treatment at least 8 days before scanning. Controls were at least abstinent 10 days before scanning. All participants provided written informed consent to the experimental procedure, which was approved by the Ethical Committee of ICS Maugeri (Pavia, Italy). The investigation was conducted in accordance with the latest version of the Declaration of Helsinki.

Neuro-cognitive assessment

All participants underwent a neuro-cognitive evaluation based on the Brief neuropsychological examination44, a well-validated battery for the Italian population including tasks for different cognitive domains: attention (trail making (TMT-A and TMT-B)), memory (digit span, immediate and delayed prose memory), working-memory (10- and 30-seconds interference-memory), executive functions (TMT-B, cognitive estimation, abstract reasoning, phonemic fluency, clock drawing and overlapping pictures), as well as perceptive and praxis skills (praxis abilities, spontaneous drawing and copy drawing task). The battery results in a score for every task, alongside an overall score of global cognitive status.

Statistical analysis of neuro-cognitive data

For each task, we checked the normality of the score distribution across the whole sample. Based on this assessment, we then examined age and group effects by means of parametric or non-parametric two-sample and correlation tests. For the tasks showing a significant effect of age or education, we run an Analysis of covariance (ANCOVA) to assess group differences on cognitive performance after removing their effect. We applied a primary statistical threshold of p < 0.05, one-tailed due to a priori hypotheses of cognitive impairment in AUDs7,45, and then performed a correction for multiple comparisons based on FDR.

We investigated superordinate cognitive domains, transcending specific tasks, in which performance was impaired in patients. After assessing the suitability of the correlation matrix (Keiser-Meyer-Olkin Measure of Sampling Adequacy = 0.61; Bartlett’s test of sphericity <0.001; Supplementary Table S7) we performed a principal component analysis on the 15 ENB2 raw scores. Due to the ambiguity of the scree plot (Supplementary Fig. 1), we used the Kaiser-Guttman criterion to determine the number of components to be retained (i.e., components with eigenvalue >1). An orthogonal rotation (Varimax) was used to facilitate the interpretation of components20 (Supplementary Tables S8,S9). To investigate group differences in cognitive performance, we used an ANOVA (with Bonferroni correction for multiple comparisons) on the resulting factor scores for each subject/component.

MRI data acquisition

We used a 3 Tesla General Electrics Discovery scanner to collect a high-resolution 3D T1-weighted IR-prepared FSPGR (BRAVO) brain scan acquired along the AC-PC plane (152 slices, FOV = 24 cm, reconstruction matrix = 256 × 256, slice thickness = 1 mm). A T2-weighted image was also collected for diagnostic purposes.

VBM data pre-processing and whole-brain statistical analysis

We performed image pre-processing and statistical analyses using SPM12 (http://www.fil.ion.ucl.ac.uk/spm) and the CAT12 toolbox (http://www.neuro.uni-jena.de/cat/). The pre-processing included a correction for bias-field inhomogeneities, spatial normalization using the DARTEL algorithm46, segmentation into GM, white matter (WM) and cerebrospinal fluid (CSF)47, and smoothing with an 8 mm gaussian kernel.

Statistical analyses included: (a) a two-sample t-test, to assess a decrease of GM density in alcoholic patients vs. controls; (b) multiple regressions, to assess a relationship between GM density and performance in the domain/task displaying the greatest impairment in patients; and (c) full factorial models (two-sample t-test plus a behavioural covariate), to assess group differences in the relationship between GM density and performance (i.e., a significantly different regression slope in patients vs. controls). For (b) and (c), in separate analyses we modelled either the factor score of the basic-level executive component, or response-time of the TMT-A task to examine its contribution to the component. Multiple regressions and full factorial models highlighted, respectively, quantitative or qualitative group differences in the relationship between GM density and domain/task performance. We modelled age to remove its potentially confounding effect, and applied an internal GM threshold of 0.15 to prevent artefacts on the GM-WM border due to voxel misclassification. We used conjunction-null analyses48 to assess the predicted anatomical overlap between the regions in which GM density was both reduced in patients vs. controls, and related to executive or TMT-A performance.

Since the above analyses involved two behavioural measures, we adjusted our primary statistical threshold to p < 0.025 corrected for multiple comparisons with FDR (as implemented in SPM12) at the voxel or cluster level. We applied threshold-free cluster enhancement (TFCE49) with 5000 permutations per contrast and correction for multiple comparisons. This approach has been shown to increase the sensitivity of VBM findings50.

Regions-of-Interest statistical analyses

We investigated whether, and to what extent, the pattern of cognitive impairment observed in alcoholic patients is explained by the degree of regional GM atrophy. We used the SPM toolbox Marsbar (http://marsbar.sourceforge.net/) to create binary masks of the clusters displaying the effects reported above, i.e. a) GM atrophy in patients vs. controls; b) common effect of GM atrophy in patients and correlation with basic-level executive performance. Using the SPM toolbox REX (http://web.mit.edu/swg/software.htm), the average GM density in these regions was extracted for each subject, and entered in offline analyses. Namely, we used average GM density in the observed clusters as simultaneous predictors of a multiple regression model, to assess their global and relative efficacy for predicting executive performance.

Meta-analytic evidence of an overlap with the salience network

We used the Neurosynth toolbox (http://neurosynth.org) to produce a meta-analytic map of the salience network (see Supplementary Methods), and then examine its spatial overlap with the regions displaying common effects of interest (atrophy in patients, correlation with executive performance and TMT-A response time).

We used the Marsbar toolbox, as described above, to create spatial maps corresponding to the conjunction of our effects of interest, i.e. inclusion in the salience network alongside significant GM atrophy in AUDs, correlation with executive performance and with TMT-A response time. Then, to evaluate the extent to which the morphometric properties of the resulting regions account for cognitive performance, we replicated the procedure described above to extract average GM density from the commonly involved voxels, for offline multiple regressions.