Predicting development of adolescent drinking behaviour from whole brain structure at 14 years of age
- University Medical Center Hamburg-Eppendorf, Germany
- Max Planck Institute for Human Development, Germany
- Heidelberg University, Germany
- Trinity College Dublin, Ireland
- University Medical Centre Hamburg-Eppendorf, Germany
- King’s College London, United Kingdom
- University of Mannheim, Germany
- Université Paris-Saclay, France
- University of Vermont, United States
- University of Nottingham, United Kingdom
- Charité – Universitätsmedizin Berlin, Germany
- Physikalisch-Technische Bundesanstalt (PTB), Germany
- Institut National de la Santé et de la Recherche Médicale, INSERM Unit 1000 “Neuroimaging & Psychiatry”, University ParisSud, University Paris Descartes, France
- Holland Bloorview Kids Rehabilitation Hospital, Canada
- University of Toronto, Canada
- University Medical Centre Göttingen, Germany
- Technische Universität Dresden, Germany
Abstract
Adolescence is a common time for initiation of alcohol use and development of alcohol use disorders. The present study investigates neuroanatomical predictors for trajectories of future alcohol use based on a novel voxel-wise whole-brain structural equation modeling framework. In 1814 healthy adolescents of the IMAGEN sample, the Alcohol Use Disorder Identification Test (AUDIT) was acquired at three measurement occasions across five years. Based on a two-part latent growth curve model, we conducted whole-brain analyses on structural MRI data at age 14, predicting change in alcohol use score over time. Higher grey-matter volumes in the caudate nucleus and the left cerebellum at age 14 years were predictive of stronger increase in alcohol use score over 5 years. The study is the first to demonstrate the feasibility of running separate voxel-wise structural equation models thereby opening new avenues for data analysis in brain imaging.
https://doi.org/10.7554/eLife.44056.001eLife digest
Puberty is a time of transformation. Physical changes in the body occur alongside changes in personality and behaviour. Compared to children, adolescents tend to be risk-takers and novelty-seekers. They crave new sensations and experiences, as well as social interaction with their peers. It is around puberty that many people try alcohol for the first time. But it is not clear why people differ in their drinking habits, and why a small minority of young adults go on to become dependent on alcohol.
Part of the answer may lie in changes in the brain. Differences in the size and structure of brain regions contribute to differences in behaviour between individuals. During adolescence, the brain undergoes extensive re-modelling. It forms new connections, while also pruning away connections that are unused. Could differences in brain structure at puberty lead to differences in alcohol consumption in early adulthood?
Kühn et al. scanned the brains of about 1,800 healthy adolescents at the age of 14 and then again at 19 (within the context of the IMAGEN study). At three time points, the teenagers also filled in questionnaires about their use of alcohol. Two areas of the brain – the caudate nucleus and the left cerebellum – were larger at age 14 in teenagers who would increase their alcohol consumption by age 19. The larger the areas at age 14, the bigger the increase in alcohol consumption over time. Notably, there was no relationship between the size of either brain area at the age of 14 and how much alcohol the individuals drank at the same age.
These results may help us to understand why some young adults develop harmful drinking habits, whereas most do not. The findings are part of a large and complex picture. Other factors, such as social influences, also shape alcohol consumption. However, the findings of Kühn et al. suggest that differences in brain structure may make some individuals more likely to increase how much alcohol they drink than others. Understanding these biological differences could help researchers to develop measures to prevent addiction in young adults.
https://doi.org/10.7554/eLife.44056.002Introduction
Adolescence is a critically vulnerable time for the development of alcohol drinking habits that may lead to considerable consequences later in life including the development of alcohol addiction. Importantly, the period of adolescence coincides with substantial behavioural changes together with structural and functional brain development. Cortico-striatal regions play an important role in the regulation of behaviour and might therefore play a role in progress and maintenance of habits such as drinking (Heinz, 2002). In particular, it has been proposed that drug addiction involves dysfunctions of brain circuitry related to the neurotransmitter dopamine that lead to alterations in both impulsive and compulsive behaviour (Koob and Kreek, 2007). Since early exposure to drugs may alter brain development during adolescence this may set the stage for cognitive problems in adulthood, which translate into behavioral consequences throughout life (Jacobus and Tapert, 2013). Hence, it is of particular importance to predict the acceleration of alcohol use as early as possible during adolescence in order to intervene timely.
Nees and colleagues reported that reward-related brain activation aided in the prediction of early-onset drinking in adolescents at age 14 years in a subset of the data used in the present paper (Nees et al., 2012). In a functional imaging study on youth of age 12–14 years, prior to initiation of alcohol use, it was found that teens classified as transitioning to heavy alcohol use by age 18 had less blood oxygen level dependent (BOLD) activation in frontal, temporal, and parietal cortices in a response inhibition task (Norman et al., 2011). These reports support the notion that brain differences present early during adolescence may leave certain youth vulnerable to addictive behaviors. In an earlier publication on the same data set, we focussed on the prediction of changes in alcohol-related problems between age 14 and 16 years based on gyrification of the orbitofrontal cortex (Kühn et al., 2016). We argued that it is important to use behavioral difference scores when aiming at predicting prospective behavior rather than absolute measures at the prospective time point only. In applying difference scores rather than absolute measures, all available information can be taken into account to evaluate potentially problematic behaviour.
Within the present study, we want to go beyond previous prediction attempts by including a trajectory of alcohol-related behavior over a longer period of time, from age 14 to 16 and 19 years of age. Obtaining three measurement occasions enables the modelling of change on a latent level within a structural equation modelling (SEM) framework. We set out to investigate brain structural predictors at age 14 of the latent change trajectory of the alcohol use score on a voxel-based whole-brain basis. Methodologically this goes beyond previous studies with an SEM approach including brain data, since commonly only extracted regions of interests or global brain measures (such as intracranial volume, white matter hyperintensities) have been employed in an SEM context (Ritchie et al., 2015; Kievit et al., 2014). The innovation of this paper’s approach lies in analysing the whole brain on voxel level. In doing so, we are able to on the one hand extract possible predictors for change in alcohol consumption on a fine-grade level (i.e. voxel), and on the other hand to exploratory analyse the whole brain in search of predictive regions. This approach offers the major advantage of giving consideration to both, specific structures of the brain on a rather fine-graded level, and the broad, exploratory search for generally influential areas.
Results
We started the analyses with estimating the two-part latent growth curve model on the clinical data only containing an intercept and a linear growth factor for both alcohol use vs. non-use as well as for the alcohol use score. Table 1 shows the classification of the AUDIT-scores concerning severity of use (Figure 1).
Figure 1
Preparation of AUDIT Sum-Scores for two-part latent growth mixture model.
Upper row: data from original scale (Sum Score), zeros are shown in black and indicate non-drinking individuals. Middle row: Transformation of data into consumer and non-consumer without fine-grading of alcohol use scores. Bottom row: Alcohol use score (AUDIT Sum-Score) for individuals who drink at all. Note that to enhance readability of the figure, sum-scales (upper and bottom row) are truncated at a score of 20.
Table 1
Severity of alcohol use at three measurement occasions according to AUDIT.*
https://doi.org/10.7554/eLife.44056.004| Baseline n** = 1794 (100%) | Follow-up 1 n = 1439 (100%) | Follow-up 2 n = 1284 (100%) | |
|---|---|---|---|
| No use at all | 855 (47.7%) | 255 (17.7%) | 95 (7.4%) |
| Unproblematic use | 872 (48.7%) | 961 (66.7%) | 823 (64.1%) |
| Medium level of alcohol problems | 64 (3.6%) | 218 (15.3%) | 329 (10.7%) |
| High level of alcohol problems | 2 (0.1%) | 4 (0.3%) | 28 (2.2%) |
| Indicating dependence | 1 (0.1%) | 1 (0.1%) | 9 (0.9%) |
-
*Note: Categorization is based on the interpretation guideline of the World Health Organization: Cut-offs scores are: 0–7 = unproblematic use, 8–15: simple advice focused on the reduction of hazardous drinking warranted, 16–19: brief counseling and continued monitoring warranted, above 20: further diagnostic for alcohol dependence strongly warranted.
**Note: 20 individuals had missing data, in total adding up to 1814.
In the discrete part of the model, a linear growth model demonstrated better fit to alcohol use vs. non-use than an intercept-only model, Δχ2 (Koob and Kreek, 2007) (3, N = 1814)=653.63, p<0.001. Inclusion of a quadratic growth factor did not improve model fit, Δχ2 (Koob and Kreek, 2007) (1, N = 1814)=−0.74, p=0.390, therefore we refrained from a quadratic growth factor.
In the continuous part of the model, a linear growth model likewise demonstrated a better fit to alcohol use scores than an intercept-only model, Δχ2 (Koob and Kreek, 2007)(3, N = 1814)=860.58, p<0.001. Inclusion of a quadratic growth factor did not improve model fit, Δχ2 (Koob and Kreek, 2007)(1, N = 1814)=−0.58, p=0.450. Therefore, we accepted the model with intercept and a linear growth factor on the clinical data as our working model. We then added the nuisance covariates age, sex, scanner and total brain volume to the model by predicting the slope of the continuous part of the model (no matter whether the regression paths were significant or not), since it is common practice in neuroimaging studies to control for these nuisance variables. The final model on the clinical data, not yet including the brain data (since this varied for each voxel of the brain) demonstrated an acceptable model fit, χ2 = 444, df = 65, RMSEA = 0.057, CFI = 0.785, SRMR = 0.065.
Concerning the latent intercept and slope for both parts of the model, intercepts were significantly different from zero and variances of the intercepts for both parts of the model were significant, suggesting significant interindividual heterogeneity around the estimated mean level of alcohol drinking use vs. non-use and the alcohol use score at age 14 (for estimates see Table 2). The covariance between the intercepts of the two parts of the model was 0.124, p<0.001, indicating that adolescents with a higher propensity to engage in alcohol drinking also engaged in it more frequently and vice versa.
Table 2
Estimated parameters in probability of use vs. non-use and alcohol use score with nuisance variables on the clinical data (not yet including brain data)
https://doi.org/10.7554/eLife.44056.005| Intercept | Slope | |||
|---|---|---|---|---|
| Estimate | SE | Estimate | SE | |
| Part 1: Prevalence of alcohol drinking (use vs. non-use)=discrete part of the model | ||||
| Mean | 0.568** | 0.011 | 0.188** | 0.006 |
| Variance | 0.090** | 0.009 | 0.024** | 0.004 |
| Part 2: Alcohol use score of AUDIT = continuous part of the model | ||||
| Mean | 0.693** | 0.037 | 0.498 | 0.642 |
| Variance | 0.618** | 0.087 | 0.218** | 0.046 |
| Regression onto Part two slope | ||||
| Sex | −0.183** | 0.046 | ||
| Age | −0.000 | 0.000 | ||
| TBV | 0.000* | 0.000 | ||
| Site_London | 0.410* | 0.163 | ||
| Site_Nottingham | 0.368* | 0.161 | ||
| Site_Dublin | 0.517* | 0.167 | ||
| Site_Berlin | 0.091 | 0.170 | ||
| Site_Hamburg | 0.122 | 0.162 | ||
| Site_Mannheim | 0.038 | 0.163 | ||
| Site_Paris | 0.079 | 0.163 | ||
| Site_Dresden | −0.044 | 0.163 | ||
| Covariances | ||||
| Covariance between intercept and slope in Part 1 | −0.033** | 0.005 | ||
| Covariance between intercept and slope in Part 2 | −0.078 | 0.050 | ||
| Covariance between the intercepts of Part 1 and Part 2 | 0.124** | 0.012 | ||
-
*p < 0.05 , **p<0.001, SE = standard error, TBV = total brain volume
Turning to the growth parameters, for the continuous part of the model the estimated mean of the slope was not significantly different from zero, indicating that, on average, no change in drinking habits emerged over time. However, the variance was significantly different from zero, indicating interindividual differences in change of drinking behaviour between participants. For the discrete part of the model both the mean and the variance of the slope were significantly different from zero, indicating change on average as well as interindividually. The positive mean of the slope indicated an increasing propensity for drinking across time. Intercept and slope covaried significantly (−0.033, p<0.001) within the discrete part of the model, however, numerically the correlation coefficient was so small that we do not think the association necessitates in-depth interpretation. No significant correlation emerged between intercept and slope for the continuous part of the model, indicating that a relation between alcohol use score at age 14 and change in this behavior was not captured in our model.
Based on this two-part latent growth curve model on the clinical data, we entered the brain data and conducted a whole-brain analysis on grey matter probability maps at age 14 years predicting change in alcohol use score, that is the latent slope in the continuous part of the model, over time from each voxel in the brain (Figure 2). Within the model an increase in alcohol use score is reflected as a positive and a decrease in alcohol use score as a negative latent slope mean. We found a positive association within bilateral caudate nucleus (around −14, 24, 7 and 17, 24, 6) and left cerebellum (Lobule VIII/IX, around −16,–53, −56), where higher grey matter volume predicts greater change in alcohol use scores (p<0.001, cluster >100 voxels, Figure 3). We repeated the same analysis on log-transformed data, to address the problem of skewness in the data with almost identical results. The same analysis on white matter probability maps did not result in any significant clusters. Neither did repeating the same analyses when investigating a regression path between white or grey matter voxels and the change in alcohol use vs. non-use. Moreover, we conducted the same whole-brain analysis while predicting the latent intercept of the continuous part of the model, which reflects how high individuals score on AUDIT on average. Here also no significant clusters emerged, neither for white nor grey matter maps.
Figure 2
Two-part latent growth mixture model.
c = continuous, d = discrete, BL = baseline, FU = follow up, I = intercept, S = slope, TBV = total brain volume, MRI site was not a single indicator as depicted for reasons of simplicity, but consisted of 9–1 separate indicators dummy coding the different scanners used.
Figure 3
Brain regions showing a significant regression path from brain voxel to the latent slope of alcohol use score increase over time.
The higher the grey matter volume the larger the slope increase.
Discussion
The goal of the present study was to unravel structural brain predictors at age 14 of the trajectory of alcohol use scores over the course of adolescence, namely between the age of 14 and 19 years. For this reason, we used a two-part latent growth curve model since it decomposes the semicontinuous outcome measure into a dichotomous use vs. non-use and a continuous alcohol use score part. The mean and variance of the intercepts for both parts of the model were significant, and the covariance between both intercepts was significant, indicating that adolescents with a higher propensity to engage in alcohol drinking also engaged in it more frequently and vice versa. For the continuous part of the model, the estimated mean of the slope was not significantly different, indicating that, on average, no change in drinking habits emerged. However, the variance was significant, indicating interindividual differences in change of drinking behaviour between participants. For the discrete part of the model, both the mean and the variance of the slope were significantly different from zero.
Since we were most interested in early brain-predictors of the intraindividual changes of alcohol use scores we computed a separate SEM for each brain voxel acquired at baseline (age 14 years). To obtain a brain map we plotted the resulting statistics of the regression path between brain and changes in alcohol use score back into brain space where we observed that higher grey matter volume in bilateral caudate nucleus and in left cerebellum was associated with a stronger increase (slope) in alcohol use scores. No associations were observed between grey matter brain data and the slope or intercept of the dichotomous use vs. non-use or the intercept of the continuous part of the model, nor in the white matter of the brain.
Neurodevelopmental changes in caudate and cerebellum
The direction of the association, namely between increase in alcohol use scores and higher brain volume is remarkable, since it may reflect a disturbance of brain development during adolescence. Caudate and cerebellum have been described as undergoing changes over the lifetime that resemble an inverted U-shape in gray matter volume that peaks during adolescence (Durston et al., 2001; Brenhouse and Andersen, 2011). However, on a longitudinal data set the changes of caudate over the course of adolescence were the smallest in comparison to other subcortical brain structures with caudate, putamen and nucleus accumbens peaking at earlier ages than amygdala and hippocampus (Goddings et al., 2014). In tendency females seem to show a peak around age 11 years and males seem to show a decrease in caudate volume over adolescence (Goddings et al., 2014; Brain Development Cooperative Group, 2012). The pattern of changes in the cerebellum seem less clear with one study reporting decreases across adolescence (Ostby et al., 2009) and another showing an inverse U-shape pattern with a peak around 15 years of age (Brain Development Cooperative Group, 2012). However, the literature on age related changes in brain volume during adolescence does not help to solve the question whether our observed effects reflect a deceleration of pruning in bilateral caudate and cerebellum or an overproduction that surpasses the normal overproduction of synapses.
Caudate nucleus and cerebellum as a predictor of alcohol problems
The present finding of a predictive value of bilateral caudate volume for the trajectory of alcohol use during adolescence fits nicely to previous studies in search of brain structural predictors of drinking.
One study on 40 adolescents using FreeSurfer showed that at baseline (12–17 years of age) participants who transitioned into heavy drinking after 3 years showed smaller left cingulate, pars triangularis, and rostral anterior cingulate volume, and less right cerebellar white matter volumes (Squeglia et al., 2014).
Other studies focussed on activation of the striatum, for example during an fMRI reward paradigm, reporting increases in brain activation in caudate nucleus at reward delivery in participants who had their first drink during puberty as compared to those who started after puberty (Boecker-Schlier et al., 2017). In another study focussing on brain function, 43 postpubertal participants (aged 18–21 years) were followed over a year and classified into moderate and heavy drinkers or transitioners who started drinking heavily. In a cue-reactivity paradigm showing alcohol stimuli, transitioners showed higher activation in bilateral caudate nucleus, orbitofrontal cortex medial frontal cortex/anterior cingulate and left insula (Dager et al., 2014). However, the two latter results have not been corrected for grey matter volume, and it could well be that the stronger brain activation observed in caudate (among others) was actually driven by higher grey matter volume (to correct for these structural effects dedicated toolboxes have been developed; Casanova et al., 2007).
Of note, it is very interesting that the cluster in caudate nucleus and cerebellum found to be predictive of the future changes in alcohol use scores were not related to the mean level of alcohol use scores in each individual, since we observe no significant regions in the analysis where the latent intercept of the continuous part of the model predicts the respective brain voxels. This indicates that the grey matter volume in caudate nucleus and cerebellum has a value in predicting changes, but not the average alcohol use in general at the timepoint analysed. On the basis of the present analyses, one might argue that preventive measures at this age should focus on the development of alcohol use rather than on drinking habits at the age of 14, as the results suggest that specific characteristics of those brain regions prepare the ground for future alcohol consumption.
Prediction of the start of alcohol drinking
The fact that the brain-based prediction of the latent slope of the dichotomous part of the model showed no significant resulting clusters indicates that, using the presented methods, brain data is not indicative for the prognosis when individuals start drinking. In this sense, the present results could be regarded as in line with a previous publication from the IMAGEN data set (Whelan et al., 2014) in which it was shown that brain information (at age 14 years) added only moderately in the prediction of binge drinking at age 16 years. Moreover, rather unspecific brain measures have been added in the analyses, namely overall regional grey matter volume and the ratio of grey and white matter volume. In a different data set (Squeglia et al., 2017), likewise only diffuse regions added to the prediction of initiation of alcohol use during adolescence. This later observation makes the prediction of the changes in alcohol use scores over time from very distinct brain regions the more remarkable. While we learned from previous studies that brain characteristics might not be the key factor in explaining the start of alcohol consumption in adolescence, we learn from the present study that brain structural characteristics are of relevance when considering the development of alcohol use in adolescence.
Caudate nucleus and psychiatric disease
Interestingly, a recent meta-analysis on brain imaging studies focussing on brain structural alterations across psychiatric disorders has revealed consistency in increases within bilateral striatum when comparing psychiatric patients to controls (Goodkind et al., 2015). From this meta-analysis on cross-sectional data, it is unclear whether these striatal increases in psychiatric patients were already present during adolescence or whether they occurred around disease onset or over the course of the disease. But it is interesting that the direction of the effect and the localisation bear resemblance to the results of our study, although we focussed on trajectories of alcohol use which were, for most of the participants, far from the actual diagnosis of alcohol addiction.
Novel methodological approach: whole-brain structural equation modeling
To our knowledge up to now structural equation modeling on brain imaging data has been conducted solely based on data derived from regions of interest (ROIs) (Kievit et al., 2014; Kühn et al., 2017; McArdle et al., 2004; Raz et al., 2005). However, this has the strong disadvantage that the obtained results are restricted to the regions selected for the analysis at hand. The present study demonstrates the feasibility of running separate structural equation models for each and every voxel of the brain and therefore plot the voxel-wise results of structural equation models back into brain space. This approach is not restricted to growth curve models but can be applied to all models in an SEM context. It can for example also be applied to measurement models, in order to relate brain structure or function not only to a separate task performance measure but rather to the latent factor representing the shared variance of a set of different performance measures from the same cognitive domain. This offers a new avenue of structural equation modeling on neuroimaging data in an unbiased, comprehensive way.
Conclusion
The present study revealed structural brain predictors (at 14 years of age) of the trajectory of alcohol use scores between the age of 14 and 19 years. A two-part latent growth curve model was utilized to decompose the semicontinuous AUDIT outcome measure into a dichotomous use vs. non-use and a continuous alcohol use scale part. We predicted the slope of use vs. non-use and of alcohol use scores by voxel-wise grey and white matter probability maps at baseline. To obtain brain maps as a result, we plotted the statistics of the regression path between brain and slope back into brain space. We observed that higher grey matter volume in bilateral caudate nucleus and in left cerebellum at age 14 years was associated with a stronger increase in alcohol use scores. This finding fits well to previous studies pointing at an association between increases in striatum and psychiatric disease. Potentially this is due to neurodevelopmental interindividual differences since adolescence is a period of brain structural in- and decreases. Our finding may reflect a deceleration of pruning or an overproduction that surpasses the normal developmental overproduction of synapses. Future research with repeated neuroimaging measurements is needed to solve this neurodevelopmental question.
Materials and methods
Participants
We used data of 1794 healthy 14-year-old adolescents (mean age = 14.4, SD = 0.45 years; 54% males) who were recruited within the scope of the IMAGEN project, a European multi-centre genetic-neuroimaging study in adolescence (Schumann et al., 2010). The selection of the participants was based on the fact that structural imaging data at age of 14 years was present. At the time of analyses reported here (age 14 years), retest data at age 16–17 years was present for 1439 participants (mean age = 16.6, SD = 0.64 years; 55% males) and at age 19 years for 1284 participants (mean age = 19.0, SD = 0.77 years; 53% males). Written informed consent was obtained from all participants as well as from their legal guardians. The adolescents were recruited from secondary schools. The study was approved by all local ethics committees separately (in Germany this was accomplished by the medical ethics committee of the University of Heidelberg, reference number: 2007-024N-MA) and approved by the head teachers of the respective schools. Participants with a medical condition or neurological disorders were excluded. All participating subjects were assessed by means of self-rating and two external ratings (by their parents and a psychiatrist specialized in pediatrics) based on ICD-10 as well as DSM-IV (The Development and Well-Being Assessment Interview, DAWBA; Goodman et al., 2000).
Questionnaire
Request a detailed protocolWe administered the Alcohol Use Disorder Identification Test (AUDIT, Babor and Higgins-Biddle, 2001) at Baseline (age 14 years), Follow-up 1 (age 16–17 years) and Follow-up 2 (age 19 years) to identify alcohol use. We computed the total score by adding the scores of all 10 items.
Scanning procedure
Request a detailed protocolStructural MRI was performed on 3 Tesla scanners from three manufacturers (Siemens: five sites; Philips: two sites; and General Electric: two sites). The details of the entire MR protocol are described elsewhere (Schumann et al., 2010). In this study, we used the T1-weighted images. These high-resolution anatomical MRIs were obtained using a three-dimensional magnetization prepared gradient-echo (MPRAGE) sequence based on the ADNI protocol (http://adni.loni.usc.edu/methods/documents/mri-protocols/; modified for the IMAGEN study to give a 1.1 × 1.1 × 1.1 mm3 voxel size).
Voxel-based morphometry
Request a detailed protocolFor the present report, structural MR data of 2072 adolescents were available. We excluded all participants where the image quality was suboptimal, most likely due to movement. The visual quality control was carried out by 10 independent raters.
Structural data acquired at age 14 years was preprocessed by means of the VBM8 toolbox (http://dbm.neuro.uni-jena.de/vbm.html) and SPM8 (http://www.fil.ion.ucl.ac.uk/spm) with default parameters. The VBM8 toolbox involves bias correction, tissue classification and affine registration. The affine registered grey matter (GM) and white matter (WM) segmentations were used to build a customized DARTEL (diffeomorphic anatomical registration through exponentiated lie algebra) template. Then warped GM and WM segments were created. Modulation was applied in order to preserve the volume of a particular tissue within a voxel by multiplying voxel values in the segmented images by the Jacobian determinants derived from the spatial normalization step. In effect, the analysis of modulated data tests for regional differences in the absolute amount (volume) of GM/WM. Images were smoothed with a FWHM (full-width at half maximum) kernel of 8 mm.
Structural equation modelling
Request a detailed protocolAnalyses were conducted within a structural equation modeling (SEM) framework using MPlus and R. We implemented a two-part latent growth curve model (Muthen, 2001; Olsen and Schafer, 2001) since the AUDIT scores were zero inflated. In a two-part latent growth curve model, zeros are valid values with its own meaning and not just proxies for missingness. Information that is contained by zeros and the specific values of the non-zeros is qualitatively different and might even be differentially influenced by covariates (cf., Olsen and Schafer, 2001). In a two-part latent growth curve model, the presence or absence of a behavior and, if present, the manifestation of a specific behavior can be modeled simultaneously in one model. In our study, the original distribution of the alcohol use variable (AUDIT-score) was decomposed into two parts (see Figure 1). Then, each was modeled by separate, but correlated, growth functions (see Figure 2). For the discrete part of the model scores of zero were separated from the rest of the distribution by creating a binary indicator variable that distinguished any positive alcohol use score (=1) from nouse (=0) (lower part of Figure 2). For the continuous part of the model, the continuous indicator variables representing the AUDIT score, given that it was above zero, were used (upper part of Figure 2). In this latter part of the model, substance non-use within each time point was treated as missing data, following standard assumptions of data missing at random (MAR; Little and Rubin, 1987). In that way participants who did not drink alcohol throughout the study contributed little information to the growth parameter estimates, but all information to alcohol use was used to estimate the growth parameters. We used Maximum Likelihood estimator for our analyses.
We controlled for the effects of age at baseline, sex, total brain volume (TBV) and site (by dummy coding eight of all nine neuroimaging sites, for simplicity we represent this by only one manifest variable in Figure 2) onto the slope of the continuous part of the growth model. As criteria for model fit we report Root Mean Square Error of Approximation (RMSEA), Comparative Fit Index (CFI), and Standardized Root Mean Square Residual (SRMR). Values of the CFI above 0.90 denote a well-fitting model, whereas for the RMSEA and the SRMR values less than 0.08 may be interpreted as acceptable model-fit.
Since our main interest was, which voxels of the brain from grey matter and white matter maps at age 14 years predict the continuous slope representing the increase of alcohol drinking over time, the final SEM models were conducted including a brain variable. For this purpose, we ran separate two-part growth mixture models for each and every voxel of the brain and repeated this analyses for the grey and white matter maps separately. This approach is exploratory as no specific areas of the brain are extracted (as would be the case in the more traditional analyses of ROI). While it bears the risk of inflated type-II-errors, it also enables the investigation of links between brain structures and behavior that have not been strictly established yet. For each SEM that was generated with a specific voxel as covariate, we plotted the estimate/p-value of the regression path from brain to the slope back into the brain, to be able to plot the resulting image using R (the scripts can be obtained from the first author, more flexible functions that may be helpful in conducting similar analyses can be found in the 'neuropointillist' R package, https://github.com/IBIC/neuropointillist, Madhyastha et al., 2018). This approach is not limited to growth curve models but can easily be extended to all types of models that can be implemented in MPlus). Finally, we thresholded the resulting maps with a threshold of p<0.001 and a cluster threshold of k > 100. We repeated the same analyses while predicting the discrete slope from voxels of the grey and white matter maps, although our main interest was on brain-wise predictors of the slope derived from the continuous part of the model.
Data availability
This study uses human brain data which cannot be completely de-identified. Moreover, it was not part of the written consent of the participants for the data to be publicly shared. Researchers may access the dataset through a request to the IMAGEN consortium: https://imagen-europe.com/resources/imagen-project-proposal/.
References
-
Developmental trajectories during adolescence in males and females: a cross-species understanding of underlying brain changesNeuroscience & Biobehavioral Reviews 35:1687–1703.https://doi.org/10.1016/j.neubiorev.2011.04.013
-
Anatomical MRI of the developing human brain: what have we learned?Journal of the American Academy of Child and Adolescent Psychiatry 40:1012–1020.https://doi.org/10.1097/00004583-200109000-00009
-
The development and Well-Being assessment: description and initial validation of an integrated assessment of child and adolescent psychopathologyJournal of Child Psychology and Psychiatry 41:645–655.https://doi.org/10.1111/j.1469-7610.2000.tb02345.x
-
Neurotoxic effects of alcohol in adolescenceAnnual Review of Clinical Psychology 9:703–721.https://doi.org/10.1146/annurev-clinpsy-050212-185610
-
Stress, dysregulation of drug reward pathways, and the transition to drug dependenceAmerican Journal of Psychiatry 164:1149–1159.https://doi.org/10.1176/appi.ajp.2007.05030503
-
Current methods and limitations for longitudinal fMRI analysis across developmentDevelopmental Cognitive Neuroscience 33:118–128.https://doi.org/10.1016/j.dcn.2017.11.006
-
Structural modeling of dynamic changes in memory and brain structure using longitudinal data from the normative aging studyThe Journals of Gerontology Series B: Psychological Sciences and Social Sciences 59:P294–P304.https://doi.org/10.1093/geronb/59.6.P294
-
Neural activation during inhibition predicts initiation of substance use in adolescenceDrug and Alcohol Dependence 119:216–223.https://doi.org/10.1016/j.drugalcdep.2011.06.019
-
A Two-Part Random-Effects model for semicontinuous longitudinal dataJournal of the American Statistical Association 96:730–745.https://doi.org/10.1198/016214501753168389
-
The IMAGEN study: reinforcement-related behaviour in normal brain function and psychopathologyMolecular Psychiatry 15:1128–1139.https://doi.org/10.1038/mp.2010.4
-
Brain volume reductions in adolescent heavy drinkersDevelopmental Cognitive Neuroscience 9:117–125.https://doi.org/10.1016/j.dcn.2014.02.005
-
Neural predictors of initiating alcohol use during adolescenceAmerican Journal of Psychiatry 174:172–185.https://doi.org/10.1176/appi.ajp.2016.15121587
Article and author information
Author details
Dimitri Papadopoulos Orfanos
Ulman Lindenberger
Funding
Horizon 2020 Framework Programme
- Simone Kühn
H2020 European Research Council
- Simone Kühn
Seventh Framework Programme
- Simone Kühn
Bundesministerium für Bildung und Forschung
- Simone Kühn
Max-Planck-Gesellschaft (Open-access funding)
- Simone Kühn
Jacobs Foundation (JRF 2016-2018)
- Simone Kühn
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work received support from the following sources: the European Union-funded FP6 Integrated Project IMAGEN (Reinforcement-related behaviour in normal brain function and psychopathology) (LSHM-CT- 2007–037286), the Horizon 2020 funded ERC Advanced Grant ‘STRATIFY’ (Brain network based stratification of reinforcement-related disorders) (695313), ERANID (Understanding the Interplay between Cultural, Biological and Subjective Factors in Drug Use Pathways) (PR-ST-0416–10004), BRIDGET (JPND: BRain Imaging, cognition Dementia and next generation GEnomics) (MR/N027558/1), the FP7 projects IMAGEMEND(602450; IMAging GEnetics for MENtal Disorders) and MATRICS (603016), the Innovative Medicine Initiative Project EU-AIMS (115300–2), the Medical Research Council Grant 'c-VEDA’ (Consortium on Vulnerability to Externalizing Disorders and Addictions) (MR/N000390/1), the Swedish Research Council FORMAS, the Medical Research Council, the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London, the Bundesministeriumfür Bildung und Forschung (BMBF grants 01GS08152; 01EV0711; eMED SysAlc01Z × 1311A; Forschungsnetz AERIAL 01EE1406A, 01EE1406B), the Deutsche Forschungsgemeinschaft (DFG grants SM 80/7–2, SFB 940/2), the Medical Research Foundation and Medical research council (grant MR/R00465X/1). Further support was provided by grants from: ANR (project AF12-NEUR0008-01 - WM2NA, and ANR-12-SAMA-0004), the Fondation de France, the Fondation pour la Recherche Médicale, the Mission Interministérielle de Lutte-contre-les-Drogues-et-les-Conduites-Addictives (MILDECA), the Assistance-Publique-Hôpitaux-de-Paris and INSERM (interface grant), Paris Sud University IDEX 2012; the National Institutes of Health, Science Foundation Ireland (16/ERCD/3797), U.S.A. (Axon, Testosterone and Mental Health during Adolescence; RO1 MH085772-01A1), and by NIH Consortium grant U54 EB020403, supported by a cross-NIH alliance that funds Big Data to Knowledge Centres of Excellence. SK has been funded by a Heisenberg grant from the German Science Foundation (DFG KU 3322/1–1, SFB936/C7), the European Union (ERC-2016-StG-Self-Control-677804) and the Jacobs Foundation (JRF 2016–2018). FN has been funded by a Heisenberg grant from the German Science Foundation (NE 1383/14–1).
Ethics
Human subjects: Written informed consent was obtained from all participants as well as from their legal guardians. The adolescents were recruited from secondary schools. The study was approved by the local ethics committees (in Germany this was accomplished by the medical ethics committee of the University of Heidelberg, reference number: 2007-024N-MA) and approved by the head teachers of the respective schools.
Version history
- Received: November 30, 2018
- Accepted: May 7, 2019
- Version of Record published: July 2, 2019 (version 1)
Copyright
© 2019, Kühn et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,684
- views
-
- 243
- downloads
-
- 19
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Predicting development of adolescent drinking behaviour from whole brain structure at 14 years of age
eLife 8:e44056.
https://doi.org/10.7554/eLife.44056
Further reading
-
- Neuroscience
The central tendency bias, or contraction bias, is a phenomenon where the judgment of the magnitude of items held in working memory appears to be biased toward the average of past observations. It is assumed to be an optimal strategy by the brain and commonly thought of as an expression of the brain’s ability to learn the statistical structure of sensory input. On the other hand, recency biases such as serial dependence are also commonly observed and are thought to reflect the content of working memory. Recent results from an auditory delayed comparison task in rats suggest that both biases may be more related than previously thought: when the posterior parietal cortex (PPC) was silenced, both short-term and contraction biases were reduced. By proposing a model of the circuit that may be involved in generating the behavior, we show that a volatile working memory content susceptible to shifting to the past sensory experience – producing short-term sensory history biases – naturally leads to contraction bias. The errors, occurring at the level of individual trials, are sampled from the full distribution of the stimuli and are not due to a gradual shift of the memory toward the sensory distribution’s mean. Our results are consistent with a broad set of behavioral findings and provide predictions of performance across different stimulus distributions and timings, delay intervals, as well as neuronal dynamics in putative working memory areas. Finally, we validate our model by performing a set of human psychophysics experiments of an auditory parametric working memory task.
-
- Neuroscience
The brain regulates food intake in response to internal energy demands and food availability. However, can internal energy storage influence the type of memory that is formed? We show that the duration of starvation determines whether Drosophila melanogaster forms appetitive short-term or longer-lasting intermediate memories. The internal glycogen storage in the muscles and adipose tissue influences how intensely sucrose-associated information is stored. Insulin-like signaling in octopaminergic reward neurons integrates internal energy storage into memory formation. Octopamine, in turn, suppresses the formation of long-term memory. Octopamine is not required for short-term memory because octopamine-deficient mutants can form appetitive short-term memory for sucrose and to other nutrients depending on the internal energy status. The reduced positive reinforcing effect of sucrose at high internal glycogen levels, combined with the increased stability of food-related memories due to prolonged periods of starvation, could lead to increased food intake.
