Early childhood temperament predicts substance use in young adults.

Behavioral inhibition (BI) is an important early childhood marker of risk for later psychiatric problems. The current 20-year prospective, longitudinal study focused on individual differences in this early temperament and adolescent brain function. As adolescents, 83 participants initially identified in infancy with the temperament of BI were assessed using functional imaging to examine striatal responses to incentives. Five years later, as young adults, these participants provided self-report of their substance use. Our findings show that children's early temperament interacts with their striatal sensitivity to incentives in adolescence to predict their level of substance use in young adulthood. Those young adults who, as children, showed the highest levels of BI reported the greatest substance use if, as adolescents, they also exhibited striatal hypersensitivity to incentives. These longitudinal data delineate one developmental pathway involving early biology and brain mechanisms for substance use in young adulthood.

Introduction

Behavioral inhibition (BI) is a well-studied early childhood temperament, which is an important marker of risk for later psychiatric problems.[1, 2] Young children identified during infancy with this temperament show fearful responses to novelty, heightened negative affect, and in early childhood display heightened social reticence to unfamiliar peers.[3, 4, 5, 6] Furthermore, individuals characterized by BI are at increased risk for psychopathology, particularly anxiety.[7]

Recent neuroimaging studies report that adolescents with a history of childhood BI,[8, 9, 10] as well as adolescents with social anxiety,[11] show striatal hyper-activation not only in response negative monetary incentive cues, but also to positive incentive cues. This hypersensitivity, found across studies, suggests that BI involves perturbed activation in the striatum. This brain region is linked to motivated behavior,[12] risk-taking, and substance use,[13] and these findings resonate with data on abnormal striatum recruitment by non-drug incentives.[14] Indeed, although BI individuals are sensitive to threat,[3, 4, 5, 6] research suggests that BI individuals may be also more vulnerable to later substance use compared with typical peers.[1, 15, 16]

Substance use disorders impart considerable suffering, impairment and health-care costs, affecting physical and mental health, as well as school and job performance. Because predisposing factors can be traced back to childhood, early identification of risk phenotypes is a public health priority.[17, 18] Substance use emerges in adolescence and tends to peak in young adulthood,[19] when changes in social contexts, such as unstructured socializing, present new opportunities for experimentation and risky behaviors.[20, 21] Considerable work in this area focuses on children with disruptive behavior as an early-emerging risk factor.[22, 23] Yet, other individuals also face risk for substance use. In particular, anxious adolescents and young adults may turn to alcohol and illicit substances to cope with increasing stress from the new, challenging social environments of young adulthood.[24]

Complex interactions among factors identified early in life and those that arise later are thought to ultimately predict outcome. Indeed, prior research on temperament suggests that risk for negative outcomes and psychiatric problems is found primarily among BI children who also manifest unique electrophysiological and behavioral response profiles, related to attention and threat monitoring.[25, 26, 27, 28, 29] However, these prior studies mainly considered prediction of anxiety; no such prospective research has charted the ways in which early-life brain and behavioral measures predict substance use patterns over time. Advances in neuroscience, and neuroimaging in particular, highlight notable progress that has improved our understanding of the neurobiology of substance use and risk factors related to it.[30] Although research on anxiety implicates attention and associated amygdala–frontal dysfunction, research on substance use implicates the striatum. The present study examines how adolescent striatal hypersensitivity to incentives moderates the relation between early childhood BI and young adult substance use. To our knowledge, this is the first prospective longitudinal study to examine the interaction between early childhood temperament and adolescent brain function, as these two factors are independently associated with individual differences in substance use.

Approximately 10 years after being characterized with the temperament of BI in early childhood, adolescents completed the monetary incentive delay (MID) task,[12] while undergoing functional magnetic resonance imaging (fMRI). About 5 years after this neuroimaging experiment, participants completed self-report measures of substance use. This longitudinal design allowed us to examine whether substance use could be predicted from a life-course perspective involving the combined influence of temperament in early childhood and brain function in adolescence. Our a priori hypothesis was that striatal activation to monetary incentives would moderate the relation between early childhood BI and adult substance use, such that BI would be positively related to substance use among individuals displaying high striatal activation to incentives. This hypothesis is consistent with other work on BI and risk for psychiatric problems, such as anxiety, where biological indices related to anxiety have been shown to moderate outcome.[25, 28, 29]

Participants and methods

Participants and procedures

Participants were drawn from two cohorts participating in a longitudinal study of temperament. At the study's inception, 443 infants were screened at 4 months of age to assess their reactivity to novel auditory and visual stimuli. Videotapes of the screening procedure were coded for positive and negative affect and motor activity during the presentation of novel stimuli. A complete description of the stimuli, coding procedures and inter-rater reliability is presented elsewhere.[4, 5, 6] One hundred seventy-eight (86 male), primarily white infants from two parent, middle-class homes were selected to provide a wide range of temperamental reactivity to novelty measures. Specifically, 37% of the selected infants showed high negative/high motor reactivity, 29% showed high positive/high motor reactivity, and 34% showed low reactivity. These infants were followed across childhood and adolescence and into young adulthood as part of the larger longitudinal study.

Following the initial 4-month assessment participants were observed in the laboratory at 14 months, 24 months, 4 years and 7 years of age. Of the original selected sample, 159 mothers of child participants (81 male) completed a temperament questionnaire during at least one of the assessment points across infancy and childhood (n=139 at 14 months; n=133 at 24 months; n=133 at 4 years; and n=116 at 7 years). Those with temperament data at any time point were not significantly different by sex (all P values >0.07).

During adolescence, a total of 91 participants completed the MID task[12] while undergoing fMRI. Participants were excluded because of medication use during the time of scan (n=4), motion >3 mm (n=2) or technical difficulties (n=2). As a result 83 participants (39 male, M age=16.33, s.d.=2.79, Range: 10.74–21.00) had useable fMRI data. During young adulthood, 93 participants (40 male, M age=19.85, s.d.=0.97, Range: 18.10–21.18) completed substance use measures, including the customary drinking/drug use record (CDDR[31]), youth risk behavior surveillance (YRBS[32]), and Brief Young Adult Alcohol Consequences Questionnaire (B-YAACQ[33]). A total of 65 (30 male) participants with available fMRI and substance use data were included in the main analyses.

Temperament assessment

At 14 and 24 months of age infants were assessed in the laboratory and were presented with novel and unfamiliar objects and people.[5] Furthermore, maternal ratings of social fear were collected at these ages with the Toddler Behavior Assessment Questionnaire.[34] At ages 4 and 7 years, children's reticent behavior with unfamiliar peers was measured using Rubin's Play Observation Scale.[35] In addition, mothers rated shyness at these ages with the Colorado Child Temperament Inventory.[36] To create a BI composite measure, behavioral, and maternal report scores across the four testing waves were taken from the entire cohort (regardless of participation at adolescence or adulthood) and standardized. Z-scores were then averaged to create the composite measure of BI.

MID task

The MID task[12] included two runs of 56 contiguous 6 s trials, for 112 trials total. At the beginning of each trial a cue was presented for 250 ms. This was followed by a crosshair fixation delay of 2000-2500 ms. Subsequently, a target stimulus was presented for 160–260 ms. This was followed by feedback information (1650 ms). Three types of cues informed the incentive condition: circle cues (n=48) indicated a potential monetary gain, square cues (n=48) indicated a potential monetary loss, and triangle cues (n=16) indicated the absence of incentive. Information on magnitude level of the incentive value was also carried by the incentive cues in the form of lines within the circle or rectangle: single line for small incentive ($0.20; n=32), two lines for medium incentive ($1; n=32) or three lines for high incentive ($5; n=32). Participants were asked to press a button as quickly as possible in response to the target stimulus in order to either gain the anticipated reward or avoid the potential loss. After the target disappeared, participants were notified of both their current and cumulative dollar amount. The order of trial type was random within each run. Participants were told that they would receive a percentage of the dollar amount won. Participants completed a practice MID task in the scanner, providing an estimate of their reaction time. Task difficulty was standardized to a ∼66% success rate by adjusting target duration (five difficulty levels). This approach minimized the potential confound of large performance differences across participants.

fMRI acquisition

Scanning occurred in a Signa 3T scanner (General Electric, Waukesha, WI) for the first 68 participants. For the final 15 participants, scanning occurred in a GE 3T Excite HDX scanner. Both groups used the same GE head coil. The two scanning groups did not differ on the major demographic or research factors (all Ps >0.13), other than age at time of scanning (P<0.0001) and substance use (P<0.05). Thus, in order to account for any potential impacts of scanner differences we included scanner type as a covariate in our main analysis.

A Cedrus Lumina response box (Psychological Software Tools, San Pedro, CA) recorded behavioral data. Task stimuli were projected onto a screen at the foot of the scanner bed and viewed with mirrors mounted on the head coil. Foam padding constrained head movement.

Each brain volume consisted of 30 interleaved slices, 4 mm thick, acquired in the sagittal plane using a T2*-weighted echo-planar sequence with a repetition time of 2500 ms, echo time of 23 ms, and flip angle of 90°. Voxel dimension was 3.75 × 3.75 × 4.0 mm. Matrix size was 64 × 64, and field of view was 24 cm. To allow for signal stabilization, four acquisitions were obtained before task onset. A high-resolution structural image was also acquired for each subject using a T1-weighted standardized magnetization-prepared spoiled gradient-recalled echo sequence: 1241 mm slices, 8100 ms repetition time, 32 ms echo time, 15° flip angle, 256 × 256 matrix and 24 cm field of view.

fMRI data analysis

Analysis of Functional NeuroImages (AFNI) software was used for the analyses.[37] Error trials were removed from processing. Standard preprocessing of echo-planar data included slice time correction, motion correction and spatial smoothing with a 6 mm full-width half-maximum smoothing kernel. Signal deviations >2.5 s.d. from the mean were removed using an AFNI despiking algorithm applied on a voxel-wise basis and a bandpass filtering algorithm to remove cyclical fluctuations in signal (either >0.01 or <0.15 Hz) not temporally indicative of a hemodynamic response. We then normalized blood oxygen level dependent (BOLD) signal intensity to percentage signal change using each subject's voxel-wise time series mean as a baseline.

Preprocessed time series data for each individual were analyzed by multiple regression.[38] The regression model included 32 regressors, including both cues and feedback. Of these, there were eight regressors: six regressors modeling effects attributable to residual motion (using the motion correction factors in the x, y and z planes and in the yaw, pitch and roll dimensions), and two regressors modeling baseline and linear trends for each of the two runs. Regressors of interest included cues signaling trial type (for example, potential gains and losses) and were convolved with a gamma variate function that modeled a prototypical hemodynamic response.[39] Idealized signal time courses were estimated based on onset time of different event types during the task.

Contrasts of BOLD activation were created for each subject for cues signaling incentive magnitude of (1) large $ gain/loss versus no $, (2) medium $ gain/loss versus no $, and (3) small gain/loss $ versus no $. Our main analysis targeted the interaction between BI and striatal activation in predicting substance use. Therefore, we used a region of interest (ROI) approach focused on the caudate, putamen and nucleus accumbens. ROIs were defined by anatomical boundaries provided by AFNI after spatial normalization,[40] in line with previous work with this task.[8] Specifically, we extracted values from the ROIs to reflect the moderation of striatal activation in the association between BI and substance use. In this way we hoped to minimize potential biases in the analysis of the imaging data.[41]

One contrast value was generated per subject for each ROI. This approach was used to balance Type 1 errors associated with multi-voxel testing, which might accrue from multiple testing of specific brain regions or specific event types, with Type 2 errors associated with overly stringent statistical thresholds in between-subject imaging research.[42] Moreover, this approach also is consistent with prior findings for the MID task.[8] The individual contrast values were extracted from each anatomical ROI. Given that we were interested in the role of striatal activation to both positive and negative incentives in predicting substance use, we made an a priori decision to average across region (nucleus accumbens, caudate, putamen), incentive (gain, loss) and magnitude of incentive (small, medium, large) as well as across laterality (right, left). This generated a single striatal activation composite variable for each subject. With this approach, one index of striatal hypersentivity could serve as a moderator in the relation between BI and substance use. This approach also allowed us to test our hypothesis in one regression model. However, preliminary analyses using non-averaged variables indicated a similar pattern of results.

Substance use composite

During young adulthood, participants responded to self-report measures of substance use. First, participants completed the CDDR,[31] a structured interview that assesses recent and lifetime alcohol, drug and tobacco involvement in four domains: consumption, withdrawal characteristics, psychological/behavioral dependence and consequences. The CDDR has been found to have high reliability and validity for each domain for both abusing and community samples.[31] These include internal consistency (α=0.68–0.94), test–retest reliability (r=0.86 to 0.93) and inter-rater reliability (99% agreement). Validity scores include convergent validity (r=0.24–0.79), discriminant validity (distinguish substance abusing from non-abusing adolescents Ps >0.01) and criterion validity (k=0.51–0.84). Scoring of the CDDR included calculating whether participants ever used a substance (88% used alcohol, 53% used drugs), age of regular substance use (alcohol M=17.23, s.d.=1.51), age of onset for each substance (alcohol M=15.87, s.d.=2.47, drugs M=16.27, s.d.=1.92), heavy drinking (M=1.85 drinks per episode, s.d.=1.43), as well as time taken to achieve regular smoking (years taken to increase to weekly smoking M=1.41, s.d.=1.53; years taken to increase to daily smoking M=1.56, s.d.=1.36). Each of these scores was then standardized.

Second, participants completed items from the YRBS.[32] Items chosen from the larger YRBS assessed risk behaviors, including unsafe driving, aggression, tobacco use, alcohol consumption, drug use and risky sexual behavior. Test–retest reliability has been carried out for each item with kappas ranging from 23.6 to 90.5%, with a mean of 60.7%.[32] Data were scored according to standard YRBS procedures. Items with <10% endorsement of behavior were not analyzed. Items that were scored included questions pertaining to alcohol consumption, drug use and risky sexual behavior. Although many questions contain multiple response categories, standard YRBS reports dichotomize responses into ‘no risk' versus ‘at risk.' These dichotomized scores were then included in a principal components factor analysis carried out separately for each type of behavior. This analysis confirmed four factors that were saved as variables: alcohol consumption (loadings >0.67), drug use (loadings >0.74), risky sexual behavior (loadings >0.45) and substance-use-related sexual behavior (loadings >0.70).

Finally, participants completed the B-YAACQ,[33] a 24-item measure that asses alcohol-related consequences. Dichotomous items (yes/no) were summed for a total number of alcohol-related consequences experienced in the past year. This measure has demonstrated high internal consistency in previous research with college students (α=0.89).[33] A sum total score was calculated (M=4.92, s.d.=5.32) and standardized. Subsequently, the standardized scores of the B-YAACQ and CDDR, as well as the factors extracted from the YRBS, were averaged to create a substance use composite score.

Hierarchical regression analyses

As a preliminary step, a series of t-tests were used to examine sex-related differences on all predictor and outcome variables. Significant sex-related differences were found for BI, t(43.63)=2.22, P<0.03 with males (M=0.21, s.d.=0.90) displaying greater BI than females (M=−0.20, s.d.=0.50). No other significant differences were found, all Ps>0.37. Thus, sex was included as a covariate in the hierarchical regression analysis. Next, Pearson inter-correlations were carried out between predictor and outcome variables. A correlation between striatal activation and substance use was found at a trend level, P<0.08. No other significant correlations were found, all Ps>0.62.

Although we did not find a correlation between the predictor variables (that is, BI and striatal activation),[43] we examined our a priori hypothesis regarding the potential moderating role of striatal activation in the relation between BI and substance use. We conducted a hierarchical multiple regression analysis using the BI and striatal activation composites as continuous predictors and the substance use composite as an outcome. To reduce multicollinearity and aid in interpretation, mean-centered predictors were used. Next, the interaction term was computed as the product between the BI composite measure and striatal activation. Age at scan time, sex and scanner type were entered as covariates in the first step of the regression analysis. BI and striatal activation were also entered in the first step of the regression analysis as main effects. To test for significant moderation effects of striatal activation on the link between BI and substance use, the interaction product term between BI and striatal activation was entered in the second step of the regression analysis. Although the entire regression model was examined for significance, to test the moderation hypothesis, we examined whether the second step of the model significantly increased the variance explained by the model. Interactions were probed and plotted according to the guidelines by Aiken and West.[44] High and low levels of BI and striatal activation were computed as ±1 s.d.

Results

Means and s.d. of predictor (BI and striatal activation) and outcome (substance use) variables are presented in Table 1. We conducted a hierarchical multiple regression analysis with BI and a composite measure representing striatal activation during the MID task (Figure 1) as predictor variables. BI was a continuous composite score, which was comprised of behavioral and maternal report assessments conducted at 14 months, 24 months, 4 years and 7 years of age. In this analysis, a substance use composite was regressed on these predictors, with the key hypothesis testing the BI-by-brain function interaction. This interaction significantly improved the model, ΔR2=0.08, F(1,64)=6.03, P<0.017. The full model with the interaction term was significant F(6, 64)=2.26, P<0.05, f2=0.23, with a medium-to-large effect,[45] predicting 19% of the variance in adult substance use (Table 2 and Figure 2). To decompose this interaction, follow-up regressions were conducted. The findings indicate that within participants who exhibited high striatal activation, BI was positively related to substance use, β=0.40, t(64)=2.19, P<0.03. However, when participants showed low striatal activation, BI was not related to substance use, β=−0.28, t(64)=−1.46, P=0.15. In addition, within participants with high BI, striatal activation was related to substance use, β=0.47, t(64)=2.85, P<0.01. However, within participants with low BI, striatal activation was not related to substance use, β=−0.16, t(64)=−0.84, P=0.41.

Table 1

Means (and s.d.) of predictor and outcome variables in the present sample

Variable	M (s.d.)
Behavioral inhibition	−0.01 (0.73)
Striatal activation	0.08 (0.13)
Substance use composite	−0.03 (0.61)

Figure 1

Striatal activation during the MID task averaged across incentives versus neutral.

Table 2

Hierarchical multiple regression analysis predicting substance use in young adulthood

Variables by step	R2	β (t)
Step 1 (df 5/59)	0.11
Gender		−0.04 (−0.33)
Scanner		−0.18 (−1.26)
Age at scan time		−0.07 (−0.51)
BI		0.08 (0.57)
Striatal activation		0.19 (1.54)
Step 2 (df 6/58)	0.19*
Gender		−0.07 (−0.52)
Scanner		−0.17 (−1.25)
Age at scan time		−0.03 (−0.23)
BI		0.06 (0.49)
Striatal activation		0.16 (1.30)
BI × striatal activation		0.30 (2.46)*

Abbreviation: BI, behavioral inhibition.

*P<0.017.

Figure 2

Joint effect of BI and striatal activation on substance use.

Discussion

The present longitudinal study highlights the conditions under which a complex pattern of behavior, substance use, arises from bio-behavioral interactions that unfold in a developmental context. Early childhood BI interacted with adolescent striatal activation to predict substance use in young adults. Specifically, individuals who were behaviorally inhibited during early childhood and who exhibited high striatal activation in response to incentives as adolescents reported elevated substance use as young adults.

These data are consistent with research showing striatal hypersensitivity to incentives among BI individuals.[8, 9, 10] Our findings are also in line with previous studies, suggesting a relation between BI and substance use.[1, 15, 16] Our study extends this previous work by describing functional mechanisms that contribute to this developmental pathway.

During young adulthood, individuals often move away from the home environment and enter new social settings.[46] An early-life tendency to display social reticence may render some individuals vulnerable to the use of substances to enhance their ability to socialize.[47] Substance use may be viewed by these individuals as a means for social approach and ease the transition into novel social relationships.[15] Our findings suggest that this may be particularly salient for behaviorally inhibited individuals with heightened sensitivity to motivational cues, as manifested in hyper-sensitive striatal responses to positive and negative incentives.

Our hierarchical regression used a composite variable representing activation to both positive and negative incentives in different structures within the striatum (nucleus accumbens, caudate and putamen). We found that this average measure of striatal activation moderated the relation between early childhood BI and substance use. This finding is in line with previous work showing that striatal activation is sensitive to positive and negative incentive cues, substance use and BI.[8, 9, 10, 12, 48] Specifically, anticipation of increasing monetary reward has been found to selectively recruit nucleus accumbens,[12] while the caudate and putamen seem to respond more generally to salient events.[48] Furthermore, the results of the present study are in line with motivational accounts of substance use vulnerability, according to which individuals at risk for substance use show perturbations in striatal recruitment by cues predicting non-drug-related incentives.[49]

Although our findings indicate that striatal activation moderates the link between BI and substance use, our regression model did not reveal direct links between BI and substance use, nor between striatal activation and substance use. Previous studies have shown that the relation between BI and substance use is complex and involves moderating factors such as gender and risk-taking propensity.[15] Both high[1, 16] and low[50] BI may place individuals at increased risk for substance use. Finally, prior research on adverse outcomes, such as anxiety, suggests that neurocognitive and neurophysiological factors moderate the relations between BI and outcome. This work shows that BI predicts risk primarily among youth with biological profiles that also are associated with risk.[25, 26, 27, 28, 29] In the current study, this pattern of results emerged for striatal hyper-activity. Although the regression model did not reveal a significant relation between substance use and striatal activation, these variables were found to be related at a trend level in the present sample. This finding is in line with previous work on neural correlates associated with substance use.[13, 51] More importantly, among individuals with relatively high levels of striatal activation, early childhood BI did predict substance use.

Although previous work[8] has shown a significant difference between high and low BI individuals in striatal activation to incentives, we did not find a significant direct relation between BI and the striatal activation composite in the current study. This inconsistency may be a result of the larger variability of BI in the present sample, coupled with the complexity of associations between BI and reward function in prior studies. The lack of correlation in the present study, along with findings in prior studies, suggests that the relations between BI and striatal activation are complex. This is in line with previous studies,[9, 10] suggesting that increased striatal activation involves interactions between BI and various task conditions, such as incentive valence.

Few phenotypes are as well characterized in early childhood as BI. Children with this temperament exhibit a unique pattern of physiological reactivity and behavioral response to novelty. These findings, using early childhood temperament and adolescent brain function, are the first to trace the conditions under which a complex, multi-determined adult behavior is developed.