Brain responses to anticipating and receiving beer: Comparing light, at-risk, and dependent alcohol users

Background Impaired brain processing of alcohol-related rewards has been suggested to play a central role in alcohol use disorder. Yet, evidence remains inconsistent, and mainly originates from studies in which participants passively observe alcohol cues or taste alcohol. Here we designed a protocol in which beer consumption was predicted by incentive cues and contingent on instrumental action, closer to real life situations. We predicted that anticipating and receiving beer (compared with water) would elicit activity in the brain reward network, and that this activity would correlate with drinking level across participants. Methods The sample consisted of 150 beer-drinking males, aged 18-25 years. Three groups were defined based on AUDIT scores: light drinkers (n=40), at-risk drinkers (n=63), and dependent drinkers (n=47). fMRI measures were obtained while participants engaged in the Beer Incentive Delay task involving beer- and water-predicting cues, followed by real sips of beer or water. Results During anticipation, outcome notification and delivery of beer compared with water, higher activity was found in a reward-related brain network including the medial prefrontal cortex, orbitofrontal cortex and amygdala. Yet, no activity was observed in the striatum, and no differences were found between the groups. Conclusions Our results reveal that anticipating, obtaining and tasting beer activates parts of the brain reward network, but that these brain responses do not differentiate between different drinking levels. We speculate that other factors, such as cognitive control or sensitivity to social context, may be more discriminant predictors of drinking behaviour in young adults.


Introduction
Excessive alcohol use has been associated with risky behaviour and increased mortality (Gmel et al., 2011;Hingson et al., 2009), resulting in a high economic and disease burden worldwide (Rehm et al., 2009). Improving prevention and treatment of alcohol use disorders (AUD) requires a better understanding of the underlying mechanisms. At the brain level, disrupted reward processing has been considered as one of the key mechanisms contributing to AUD, and more generally to addictive behaviours (Charlet et al., 2013;Galandra et al., 2017;Hommer et al., 2011;Luijten et al., 2017;Volkow and Morales, 2015;Wiers et al., 2007). Most studies that have probed the reactivity of the brain reward network in AUD have used monetary rewards, showing altered response patterns in the striatum and medial prefrontal cortex (although the direction of these effects remains inconsistent, see (Galandra et al., 2017;Huys et al., 2016a) for recent reviews). In contrast, fewer studies have investigated the processing of alcohol-related rewards in AUD. This is important in the light of recent literature arguing that the brain reward network responds differently to addiction versus non-addiction related rewards (Leyton and Vezina, 2013;Sescousse et al., 2013).
Previous studies investigating alcohol-related reward processing have mostly focused on visual alcohol cue-reactivity as well as alcohol tasting. Cues that have been repeatedly paired with alcohol use elicit wanting responses resulting from the activation of the brain reward network, including the ventral striatum (VS), amygdala, medial prefrontal cortex (mPFC) and orbitofrontal cortex (OFC) (Chase et al., 2011;Heinz et al., 2009;Kuhn and Gallinat, 2011;Schacht et al., 2013;Schacht et al., 2011). These responses are proposed to be over-sensitized in substance use disorders, triggering exaggerated motivation and approach behaviour towards the substance of abuse (Robinson and Berridge, 2008;Volkow et al., 2010). However, as revealed by a recent meta-analysis, reward-related brain responses to alcohol cues do not seem to distinguish non-dependent from dependent alcohol users (Schacht et al., 2013), challenging the importance of the wanting component in explaining AUD. Yet, the studies included in this meta-analysis have a few limitations. First, most of them have examined brain responses to alcohol pictures that are not followed by actual alcohol consumption and thus presumably do not carry the same incentive value as in real life. In order to investigate whether conditioned alcohol cues elicit exaggerated brain responses in AUD, it is important to ensure that these cues are predictive of real alcohol consumption. Second, it is noteworthy that many cue-reactivity studies in dependent alcohol users did not include a control group, or had relatively low sample sizes (see Table 1 in (Schacht et al., 2013)), thereby limiting the ability to detect group differences.
More recently, a few studies have examined brain responses to tasting alcohol, investigating the hedonic or 'liking' properties of alcohol (Claus et al., 2011;Filbey et al., 2008a;Filbey et al., 2008b;Korucuoglu et al., 2016;Oberlin et al., 2016;Oberlin et al., 2013). These studies have revealed that, compared with soft drinks, alcohol elicits activity in the reward-related brain network, in particular in the VS and mPFC, among heavy drinking (young) adults. Importantly, moderate to strong correlations were found between activation in these regions and the severity of alcohol use problems (Claus et al., 2011;Filbey et al., 2008a;Filbey et al., 2008b), suggesting that reward-related brain responses to the taste of alcohol may represent a neurobiological marker of AUD. Yet, these studies have employed passive tasting paradigms in which alcohol was administered in a fully predictable manner and independently of any instrumental action. This may provide an incomplete picture of alcohol reward processing in AUD, given that reward-related brain activity is heavily dependent on both the unpredictability of the reward (Delgado et al., 2005) and the requirement for instrumental action (Elliott et al., 2004;Zink et al., 2004). This is all the more important as brain activity in response to drugs of abuse is thought to depend on whether these drugs are received passively or following contingent action (Jacobs et al., 2003).
In this study, we aimed to investigate whether problematic alcohol use is associated with abnormal reward-related brain responses to anticipating, obtaining and tasting beer, while addressing the limitations of previous studies. To this aim, we used fMRI combined with a novel task design inspired by the Monetary Incentive Delay Task (Knutson and Greer, 2008). This task, that we refer to as the Beer Incentive Delay task, used alcohol instead of monetary rewards. Specifically, participants were exposed to abstract cues predicting the delivery of either beer or water (latter being used as a control condition) and had to react fast enough to a visual target in order to receive the predicted reward in the mouth via a tube. Brain responses reflecting motivation or "wanting" were measured during the period preceding the motor action (cue + delay), while brain responses reflecting pleasure or "liking" were measured at the time of reward outcome notification and reward delivery, separately.
Importantly, we recruited a large cohort of 150 participants spanning the whole spectrum from light to dependent drinkers. In order to validate our novel task design, we first tested whether beer anticipation, outcome notification and delivery would elicit higher brain responses in the reward-related network compared with water. Then we hypothesized that these brain responses, in particular in the VS, would differentiate light, at-risk and dependent alcohol drinkers.

Material and methods
Participants Participants were recruited via flyers distributed throughout the Radboud University campus and Nijmegen city, as well as via online advertisement including the University recruitment website. Potential participants completed an online screening to assess their eligibility to participate (see detailed flow-chart in Supplementary Figure 1). Inclusion criteria were: 1) age 18-25, 2) being male, and 3) drinking beer. Exclusion criteria were MRI contraindications and a history of brain injury. Participants were categorized into three groups: light, at-risk and dependent drinkers, based on Alcohol Use Disorders Identification Test (AUDIT) scores (Saunders et al., 1993). Participants with an AUDIT score < 8 were considered light drinkers, those with an AUDIT score between 8 and 15 were considered at-risk drinkers, and those with an AUDIT score > 15 were considered dependent drinkers, in line with AUDIT scoring criteria (Saunders et al., 1993). Moreover, dependent drinkers had to drink more than 22 alcoholic drinks per week (Saunders et al., 1994), and meet DSM-IV criteria for alcohol dependence as assessed with a semi-structured interview (MINI) administered by a psychologist in training (Sheehan et al., 1997;Van Vliet and De Beurs, 2007). The light and at-risk drinkers all consumed less than 22 alcoholic drinks a week. All participants participated voluntary, gave written consent and received a financial compensation of 50 euros (with an additional 10 euros for the MINI interview for the dependent drinkers). The study was approved by the regional ethical committee CMO-Arnhem-Nijmegen (# 2014/043).
The initial sample consisted of 165 individuals. Seven individuals were incorrectly included, as they did not meet the combined requirement of drinking less than 22 drinks/week with an AUDIT score between 0 and 15, or drinking more than 22 drinks/week with an AUDIT score >15. The data from these 7 participants were discarded before performing any data analysis. Five participants further dropped out during data collection.
Finally, the data from one participant was missing because of technical problems with the pumps delivering the drinks, and data from two participants were excluded due to excessive head motion in the scanner (>3mm). Characteristics of the final sample (n=150) are presented in Table 1. The light drinkers (n=40), at-risk drinkers (n=63), and dependent drinkers (n=47) were matched for age and education. The groups differed in alcohol consumption levels as well as smoking status, but not in drinking desire, measured with the Desire for Alcohol Questionnaire (Love et al., 1998).
Procedure Following the online screening, participants completed two behavioural sessions in a Bar-lab, followed by a separate fMRI session (data from the Bar-lab will be reported elsewhere, see Supplementary Table 1 for a complete overview of the data collected within the larger context of this study). FMRI data acquisition took place between 4:00 and 10:00 pm, coinciding with typical drinking hours. Participants were asked to abstain from drinking alcohol in the 24 hours preceding testing, as verified using a breath-analyzer. Participants performed two tasks in the scanner, including the Beer Incentive Delay (BID) task in which participants could earn sips of beer and water (see below). Participants consumed a glass of water before scanning to homogenize the level of thirst across participants. After scanning, a breath-analyzer was again used to check whether the blood-alcohol-levels (BACs) were below the legal .05 limit before participants were allowed to leave.
Beer Incentive Delay task We used a modified version of the Monetary Incentive Delay task (Knutson and Greer, 2008;Knutson et al., 2000), in which the rewards were 3 mL sips of either beer or water (Figure 1). The sips were delivered using two StepDos 03RC fluid pumps with tubes that were placed in the participants' mouth. In the anticipation phase, participants first saw a cue informing them about the opportunity to earn beer (yellow triangle) or water (blue square), followed by a variable delay materialized by a fixation cross.
Then a visual target appeared, and participants were instructed to respond to it as fast as possible using a button press. If the response was fast enough, positive feedback was provided in the form of a green tick (outcome notification phase), followed by the drink delivery in the mouth and then swallowing (delivery phase). When the response was too slow, a red cross was presented, ending that trial. Both the beer and water conditions consisted of 30 pseudo-randomized trials. Ten practice trials preceded the task. Reaction times during the practice trials were used to tailor task difficulty to each individual (by adjusting the time limit for responding to the target), which was further continuously adjusted online to ensure an overall success rate of ~66% in each condition (Knutson et al., 2000).
The task duration was approximately 20 min.
Behavioural analyses Reaction times on successful trials were analyzed using a mixed-ANOVA design, with Drink (beer/water) as a within-subject factor, and Group (light/atrisk/dependent drinkers) as a between-subject factor. Liking of the beer and water was assessed at the end of the fMRI session using Likert-scales ranging from 1 to 10 with the questions 'How much did you like the beer/water?'. Liking ratings were analyzed using the same ANOVA design as for reaction times. fMRI data acquisition Imaging was conducted on a PRISMA(Fit) 3T Siemens scanner, using a 32-channel head coil. Blood oxygen level-dependent (BOLD) sensitive functional images were acquired with a whole brain T2*-weighted sequence using multi-echo echoplanar imaging (EPI) (35 axial slices, matrix 64x64, voxel size = 3.5x3.5x3.0 mm, repetition time = 2250 ms, echo times = [9.4 18.8 28.2 37.6 ms], flip angle = 90°). The BOLD data acquisition sequence was updated during the course of the study, due to the discovery of MRI noise artefacts. The sequence parameters remained identical, except for the slice order which changed from ascending to interleaved. We took some measures in our analyses to 1) remove the artefacts, and 2) model the change in scanning sequence halfway through the study (see below). A high-resolution T1 scan was acquired in each participant (192 sagittal slices, field of view 256 mm, voxel size=1.0x1.0x1.0 mm, repetition time= 300 ms, echo time 3.03 ms).

fMRI data analyses
Pre-processing steps were conducted in SPM8 (www.fil.ion.ucl.ac.uk/spm). For each volume, the 4 echo images were combined into a single one, weighing all echoes equally. Standard pre-processing steps were performed on the functional data: realignment to the first image of the time series, co-registration to the structural image, normalization to MNI space based on the segmentation and normalization of the structural image, and spatial smoothing with an 8 mm Gaussian kernel. In addition, two cleaning methods were incorporated into the pipeline to ensure optimal removal of artefacts and thorough de-noising of the data: 1) a Principal Component Analyses (PCA) to filter out slice-specific noise components (Viviani et al., 2005) before pre-processing, and 2) an independent component analysis (ICA)-based automatic removal of motion artifacts using FSL (www.fmrib.ox.ac.uk/fsl) after pre-processing (ICA-AROMA; (Pruim et al., 2015a;Pruim et al., 2015b). This pipeline has previously been found to be efficient to take care of the MRI noise artefacts identified in the first half of our data (Nieuwhof et al., 2017).
After pre-processing, the data were modelled using a general linear model. The anticipation phase was modelled with a boxcar function as the combination of the cue and delay periods (duration 3.5 to 9.5 s). The outcome notification phase was modelled with separate regressors for correct and incorrect responses using stick-functions. The delivery phase was modelled with a boxcar function as the combination of the drink and swallow periods for correct trials (duration 7 s). The beer and water conditions were modelled separately. Six motion parameters were included and a temporal high-pass filter with a cutoff of 128s was applied to remove low-frequency noise. For each task phase (anticipation, outcome notification, and delivery), contrast images were calculated for beer>water and then entered in second-level analyses.
In order to validate the task, we first examined the brain activity elicited by the beer>water contrast across all participants using one-sample T-tests, separately for each task phase. We used the same procedure to further examine the brain activity elicited by each drink separately, using first-level Beta images for the beer or water condition.
To examine group differences, we performed separate one-way ANOVAs for each task phase, with Group as a between-subject factor (light/at-risk/dependent drinkers).
Additionally, we performed a regression analysis across all participants to identify brain regions in which activity elicited by the beer>water contrast would scale with a continuous measure of drinking level. We computed this continuous measure based on a combination of AUDIT and weekly drinking scores. Specifically, we used a principal component analysis (using 'PCA' in MATLAB) that reduced the correlation between these scores while retaining most of their information (Janssen et al., 2017;Jolliffe, 2002). We selected the first principal component as a composite measure of drinking level, that explained 96.0% of the common variance across the AUDIT and weekly drinking scores. The scanning sequence (before/after discovery of artefacts) was added as a binary covariate of no interest in all fMRI analyses. All T-maps were thresholded with a voxel-level uncorrected p<.001, combined with a clusterlevel family-wise error (FWE) corrected p<.05, accounting for multiple comparisons across 1 0 level FWE corrected p<.05 across the whole brain (since cluster-level correction is not available for F-maps in SPM8).
Given our a priori hypothesis about the ventral striatum (VS), region-of-interest (ROI) analyses were performed using an anatomical mask of the VS (Murray et al., 2008). Percent signal change for the beer and water conditions was extracted in this ROI using the rfxplot toolbox (Glascher, 2009) and analyzed using frequentist ANOVAs in SPSS, as well as by Bayesian statistics using JASP (Wagenmakers et al., 2017).

Behavioural results
The average success rate of 65% (SD: ± 6) approached the intended 66%. Average reaction times are shown in results suggest that participants across the three groups liked the water more than the beer in the BID task.
Imaging results All analyses referred to in this section can be accessed at https://neurovault.org/collections/TOHDTVIQ. First, we examined brain responses to beer compared with water, in the three different phases of the task and across all participants ( Figure 2). These analyses revealed that during the anticipation phase, the right medial PFC The beer versus water contrast did not elicit the expected reward-related activations in the striatum. In addition, liking ratings revealed that water was rated as more pleasurable than beer. We thus reasoned that the lack of striatal activity in the beer>water contrast might reflect the fact that the water condition elicits comparable or even higher striatal activity compared with the beer condition. To test this hypothesis, we examined brain activation patterns for the beer and water conditions separately. We found that, across the three phases of the task, the beer and water conditions indeed recruit very similar brain regions including the VS and putamen, insula, DLPFC, ACC, and somatosensory cortex (Figure 3, for other foci, see Supplementary Table 3). This observation confirms that the reward-related brain network is activated in this task, but to a similar extent in the beer and water conditions, thereby explaining why their direct contrast does not produce the expected activation in the striatum.
Then, we examined whether the groups differed in their brain responses to the beer versus water conditions in any of the three phases of the task. In contrast to our hypothesis, we did not observe any significant differences surviving multiple comparisons across the whole brain between light drinkers, at-risk drinkers, and dependent drinkers.
In order to further examine individual differences, we performed a regression analysis using our composite measure of drinking level (see Methods) as a regressor, across all participants. In line with the results of the above group analysis, this regression analysis did not reveal any brain activity in the beer>water contrast scaling with drinking level, in any of the three phases of the task.
Finally, we performed an ROI analysis restricted to an anatomical mask of the VS.
Specifically, we extracted the percent signal change for the beer and water conditions for the various phases of the task and examined potential group differences for the beer>water contrast using a one-way ANOVA ( Figure 4). Again, the results showed no group differences during the anticipation (F (2,147) =.955, p=.387), outcome notification (F (2,147) =.511, p=.601) and delivery phases (F (2,147) =.097, p=.908). In order to assess whether this lack of significant group differences can be interpreted as evidence for the null hypothesis (H0, no group difference), we performed Bayesian analyses with default flat priors in JASP (Wagenmakers et al., 2017). The Bayes factor quantifying the relative evidence in favour of H0 over H1 (significant difference between groups) was BF01=6.557, thus providing moderate evidence for the null hypothesis of no group difference. Similarly, the Bayes factors for the outcome notification and delivery phases were BF01=9.710 and BF01=13.679, respectively, indicating strong evidence for the null hypothesis of no group difference in terms of VS activation.

Discussion
Across groups, our results revealed increased brain activity in reward-related brain areas during the anticipation, outcome notification and delivery of beer compared with water, suggesting that our novel task design is well-suited to examine the processing of alcoholrelated rewards. Yet, in contrast to our hypotheses, no brain activity was found in the VS in the beer vs water comparison, and no group differences were observed between light, at-risk and dependent drinkers in any of the phases of the task. Whilst disrupted reward processing has been implicated as a core component of substance use disorders, brain responses to anticipating and receiving beer did not differentiate drinkers with different levels of alcohol use in our study. If replicated, our findings would suggest that individual differences in alcohol use may not be associated with critical differences in the processing of alcoholrelated rewards in the brain.
For validation purposes, we first examined whether the contrast of beer vs water elicited the expected activations in reward-related brain areas. Across the various phases of the task, we found such activations in the mPFC, ACC, OFC and amygdala, all part of the reward brain network (Diekhof et al., 2008;O'Doherty, 2004). More specifically, during the anticipation phase, we found higher activity for the beer vs water condition in the mPFC, ACC and OFC, areas known to be involved in the perception of craving-related stimuli In contrast to our hypothesis, we did not find activation differences between the beer and water conditions in the VS. Further analyses revealed that this was due to both beer and water conditions activating the VS to a similar extent. Unexpectedly, we also found that beerpredicting cues elicited similar reaction times (i.e. motivation) as water-predicting cues, and that liking ratings for the beer condition were lower than for the water condition. This might reflect the fact that drinking small sips of beer through a tube feels different and less pleasurable than drinking from a glass. Interestingly, our results are in line with previous reports showing that water is similar to other caloric beverages in terms subjective liking and wanting (Wegman et al., 2018), and activates a large brain network including reward-related regions (de Araujo et al., 2003). These observations call for the use of artificial saliva as a more neutral control condition (Veldhuizen et al., 2007), as well as a better evaluation of the hedonic properties of beer delivered through a tube in future fMRI studies. Importantly though, the lack of VS activity when comparing the beer and water conditions across all participants does not prevent us from examining group differences in this area, as it could be that VS activity is only present in individuals with relatively high drinking levels.
However, in contrast with this hypothesis, we found no differences between light, atrisk and dependent drinkers when examining whole-brain brain responses to beer vs water, in any of the phases of the BID task. In the VS, Bayesian statistics further provided moderate to strong evidence in favour of the null hypothesis of no group difference. While these findings appear to be at odds with major reward-related addiction theories (Bjork et al., 2012;Blum et al., 2000;Robinson and Berridge, 2001), our results are in line with a previous metaanalysis showing no differences in brain responses to visual alcohol cues between dependent and non-dependent drinkers (Schacht et al., 2013). Below we propose several explanations for the lack of group differences in our study.
First, our participants were relatively young (M age =21.45) and none of them was seeking treatment for their alcohol use, which means that our dependent group might represent a mild form of AUD. Moreover, at this age, the cumulative exposure to alcohol is still limited both in time and quantity, compared with older AUD populations. These specificities might play a role in the absence of reward-related brain abnormalities in our study, which might only arise in older, treatment-seeking AUD populations. If so, that would suggest that such abnormalities are a consequence of alcohol use rather than a predisposing factor. In support of this interpretation, a review of the literature indicates that most studies reporting positive correlations between level of use and reward-related brain activation included dependent in-treatment individuals (Hommer et al., 2011). Another feature of our study that could explain the lack of group differences is the use of abstract cues in the task.
Indeed, it has been suggested that sensitized brain responses to incentive cues in addiction would only be observed when using explicit addiction-related cues, such as alcohol pictures (Leyton and Vezina, 2013). In contrast, abstract and non-familiar cues like geometric shapes might lead to blunted reward-related brain responses. Future studies should directly test this prediction.
Alternatively, it may be that brain responses to alcohol-related rewards simply . For young adults such as those included in the present study, drinking alcohol may only be rewarding when it is accompanied by social interaction. Therefore, future research may further examine this social aspect of alcohol use in relation to reward processing.
While the relatively large sample size is a major strength of the current study, some limitations have to be acknowledged. First, the ecological validity of our task can be questioned, as individuals were receiving sips of beer through a tube, while lying down in the MRI scanner. This is obviously a different experience than having a full drink from a glass in a more relaxing setting. Yet, to date, this is the closest way to examine brain responses to the taste of beer. Second, this study is cross-sectional, and longitudinal data is needed to examine whether brain responses to beer can predict future alcohol use. Eventually, such data will provide insight into the transition from alcohol use to AUD or resilience for developing AUD.
To conclude, in this group of young adults, brain responses to the anticipation and consumption of beer were not related to individual differences in the level of alcohol use, thereby challenging the role of alcohol-related reward processing in explaining AUD. We suggest that future studies focus on the role of cognitive control and sensitivity to social context as potentially more discriminant predictors of alcohol use in young drinkers.  i  j  t  e  n  ,  M  .  ,  M  a  c  h  i  e  l  s  e  n  ,  M  .  W  .  ,  V  e  l  t  m  a  n  ,  D  .  J  .  ,  H  e  s  t  e  r  ,  R  .  ,  d  e  H  a  a  n  ,  L  .  ,  F  r  a  n  k  e  n  ,  I  .  H  .  ,  2  0  1  4  .  S  y  s  t  e  m  a  t  i  c   r  e  v  i  e  w  o  f  E  R  P  a  n  d  f  M  R  I  s  t  u  d  i  e  s  i  n  v  e  s  t  i  g  a  t  i  n  g  i  n  h  i  b  i  t  o  r  y  c  o  n  t  r  o  l  a  n  d  e  r  r  o  r  p  r  o  c  e  s  s  i  n  g  i  n  p  e  o  p  l  e   w  i  t  h  s  u  b  s  t  a  n  c  e  d  e  p  e  n  d  e  n  c  e  a  n  d  b  e  h  a  v  i  o  u  r  a  l  a  d  d  i  c  t  i  o  n  s  .  J  o  u  r  n  a  l  o  f  P  s  y  c  h  i  a  t  r  y  &  N  e  u  r  o  s  c  i  e  n  c  e   3  9  ,  1  4  9  -1  6  9  .   L  u  i  j  t  e  n  ,  M  .  ,  S  c  h  e  l  l  e  k  e  n  s  ,  A  .  F  .  ,  K  u  h  n  ,  S  .  ,  M  a  c  h  i  e  l  s  e  ,  M  .  W  .  ,  S  e  s  c  o  u  s  s  e  ,  G  .  ,  2  0  1  7  .  D  i  s  r  u  p  t  i  o  n  o  f  R  e  w  a  r I  s  t  r  i  a  t  a  l   r  e  s  p  o  n  s  e  t  o  a  l  c  o  h  o  l  c  u  e  s  :  a  h  i  e  r  a  r  c  h  i  c  a  l  l  i  n  e  a  r  m  o  d  e  l  i  n  g  a  p  p  r  o  a  c  h  .  N  e  u  r  o  i  m  a  g  e  5  6  ,  6  1  -6  8  .   S  c  h  u  l  t  z  ,  W  .  ,  2  0  1  6  .  D  o  p  a  m  i  n  e  r  e  w  a  r  d  p  r  e  d  i  c  t  i  o  n  e  r  r  o  r  c  o  d  i  n  g  .  D  i  a  l  o  g  u  e  s  i  n  C  l  i  n  i  c  a  l  N  e  u  r  o  s  c  i  e  n  c  e  1  8  ,   2  3  -3  2  .   S  e  b  o  l  d  ,  M  .  ,  N  e  b  e  ,  S  .  ,  G  a  r  b  u  s  o  w  ,  M  .  ,  G  u  g  g  e  n  m  o  s  ,  M  .  ,  S  c  h  a  d  ,  D  .  J  .  ,  B  e  c  k  ,  A  .  ,  K  u  i  t  u  n  e  n  -P  a  u        Box height represents the interquartile range (IQR), black lines represent the median, crosses represent the mean, and whiskers represent the largest and smallest values no further than 1.5*IQR. Single data points are values located outside the whiskers. Frequentists statistics show no significant differences between groups, while Bayesian statistics provide moderate to strong evidence in favour of no group differences. Note: the scales are different between the figures due to differences in the amplitude of the BOLD response for the various phases of the task.