Introduction

Emotional instability is frequent in adolescence and experiencing hypomanic symptoms is a common adolescent phenomenon that is generally transient1 and might represent a developmental stage.2 This emotional instability may be marked and considered as subthreshold bipolarity (SBP) and vulnerability to mood disorders.3, 4 Lifetime prevalence of SBP has been found to be 5.7%,5 and a significant proportion (36–45%) of adolescents having SBP symptoms have been found to escalate to bipolar disorder (BD),6, 7 particularly when symptom load increases.8 Outcomes in full-syndrome depressive and anxiety disorders in young adulthood have also been reported.9

Recent evidence indicates that early white-matter (WM) abnormalities may have a significant part in the pathophysiology of BD, and represent an early marker of the disorder.10, 11, 12 Indeed, WM alterations denoted by lower fractional anisotropy (FA), a measure of water diffusion reflecting WM bundles’ coherence, have been reported in medication-naive adolescents with BD,11 or in pediatric samples involving adolescents.10 However, studies of very early stages of the disease, or studies of vulnerable subjects are scarce. A few diffusion tensor imaging (DTI) studies have investigated youth at familial risk for BD or depression, either healthy13, 14 or having comorbidities.15 Those studies have observed decreased FA in major WM tracts such as the cingulum, the corpus callosum, the uncinate fasciculus or the superior longitudinal fasciculus;10, 14, 15 one of these studies has suggested altered WM maturation during adolescence.13

In addition, BD and depression in adolescents have been variously associated with gray-matter (GM) reductions in prefrontal and cingulate cortices, and in the amygdala–hippocampus complex.16, 17, 18, 19, 20, 21, 22, 23, 24 Increased GM volume of the left superior temporal cortex has also been observed in BD adolescents.24

Although SBP might represent a very early stage of—or a vulnerability predisposition to—bipolar, depressive or anxiety disorders, so far no study has investigated brain correlates in adolescents displaying SBP symptoms. Thus, it is unknown whether such emotional instability in adolescents is associated with variations in the maturation of WM and GM.

We hypothesized that adolescents with SBP symptoms might exhibit differences in WM bundles and in brain regions implicated in mood disorders, such as frontal–limbic networks. We sought to investigate WM microstructure using DTI, and GM and WM volumes using T1-weighted magnetic resonance imaging in healthy adolescents with SBP symptoms. To this end, we used the IMAGEN database, a European database constituted across eight study centers in France, United Kingdom, Ireland and Germany, which includes 2223 adolescents recruited in middle schools around age 14 years.25

Participants and methods

The investigation was performed in accordance with the Declaration of Helsinki. The study was approved by the ethics committees in all participating institutions. Written informed assent and consent were obtained, respectively, from all adolescents and their parents after complete description of the study.

A detailed description of recruitment and assessment procedures, and exclusion and inclusion criteria has been published.25

Participants

After screening of the IMAGEN database, 42 adolescents with SBP symptoms and 126 controls were found eligible for this study. Clinical symptoms were self-reported by the participants using the Development and Well-Being Assessment (DAWBA),26 a computerized self-report assessment that generates diagnoses based on the DSM-IV (Statistical Manual of Mental Disorders, Fourth Edition) and on the ICD-10 (10th revision of the International Statistical Classification of Diseases), subsequently validated by a group of trained raters from IMAGEN consortium psychiatrists and psychologists. In the DAWBA, presence of symptoms is rated as ‘no’, ‘a little’ or ‘a lot’.

Subthreshold bipolarity was defined as a distinct period of abnormally and persistently elevated, expansive or irritable mood associated with one or more manic or hypomanic symptoms in number and duration below criteria threshold to qualify for a diagnosis of bipolar episode.3 The SBP group was constituted by those adolescents who had no identified diagnosis on any diagnostic categories and rated ‘a lot’ on at least one cardinal manic symptom (‘expansive mood’ or ‘unstable mood’) and on at least ‘accelerated speech’ and ‘decreased sleep’ using the adolescent DAWBA symptom assessment (www.dawba.info/Bipolar). These criteria included two out of three core symptoms for the diagnosis of mania,27, 28 expansive mood and pressure of speech that have been reported in 82% of subthreshold mania cases.29 In addition, we chose the ‘decreased need for sleep’ item that has been found in 94% of cases of mania30 and in 57% of subthreshold cases.29 We discarded the hyperactivity symptom of mania, as it has been reported to overlap with attention deficit-hyperactivity disorder.31 The participants’ distribution did not differ across centers (Table 1).

Table 1 Demographic, clinical and imaging characteristics in adolescents with subthreshold bipolar symptoms compared with controls

The control adolescents were selected using the DAWBA-generated probability levels of having a psychiatric diagnosis. We chose the adolescents whose probabilities to have any diagnosis were no more than 0.5%, who had no validated diagnosis, no depressive symptoms and whose rating of the selected manic symptoms was ‘no’ or ‘a little’. They were ‘matched’ for age and sex to the SBP group using a ‘genetic search algorithm’.32

The psychometric characterization was partly conducted in participants’ homes using the Psytools software (Delosis, London, UK). Self-ratings were retrieved for substance use (Alcohol Use Disorders Identification Test; European School Survey Project on Alcohol and Drugs),33, 34 and emotion and behavior (Strengths and Difficulties Questionnaire (SDQ).35

Five participants in the SBP group, and two in the control group reported alcohol abuse, and one SBP adolescent and no control reported cannabis dependence.

In addition, the adolescents were assessed with five subtests of the Wechsler Intelligence Scale for Children (WISC-IV) to estimate intelligence quotient.

Puberty status was self-assessed using the computerized Pubertal Development Status (PDS) questionnaire,36 and family history of mood disorders was assessed with the Family Interview for Genetic Studies (National Institute of Mental Health, Bethesda, MD, USA).25

Eight participants in the SBP group (19%) and 26 in the control group (21%) reported either history of depression (recurrent or single episode) or history of BD in first- or second-degree relatives (χ2 test=0, d.f.=1, P=0.99).

Magnetic resonance imaging data acquisition

T1-weighted and diffusion tensor images were obtained on 3 T scanners from a range of manufacturers (Siemens, Munich, Germany; Philips, Best, The Netherlands; General Electrics, Chalfont St Giles, UK; Bruker, Ettlingen, Germany). The imaging protocols comparability in the different scanners was ensured through a thorough standardization.25 All participants were instructed to close their eyes and keep as immobile as possible during image acquisition.

High-resolution T1-weighted images were acquired using a Magnetization Prepared Rapid Acquisition Gradient Echo sequence based on the ADNI protocol (http://www.loni.ucla.edu/ADNI/Cores/index.shtml). Scanning parameters were similar in all centers (sagittal slice plane; repetition time=2300 ms; echo time=2.8 ms; flip angle=8°; 256 × 256 × 160 matrix; voxel size: 1.1 × 1.1 × 1.1 mm). The diffusion tensor images were acquired using an Echo Planar imaging sequence (4 b=0 and 32 directions with b-value=1300 s/mm2; axial slice plane; echo time=104 ms; 128 × 128 × 60 matrix; voxel size 2.4 × 2.4 × 2.4 mm), adapted to tensor measurements (for example, FA, mean diffusivity (MD)) and tractography analysis.

Preprocessing of T1-weighted magnetic resonance imaging data

T1-weighted image preprocessing and analysis were performed with Statistical Parametric Mapping 8 software package (SPM8; Wellcome Department of Cognitive Neurology, London, UK; www.fil.ion.ucl.ac.uk/spm). Study-specific adolescent brain tissue templates (Tissue Probability Maps) were created from a subsample of 240 adolescents randomly drawn from the IMAGEN database (15 girls and 15 boys from each center, comparable for age) to take into account the age specificity of our sample. All IMAGEN T1-weighted images were then segmented using our customized Tissue Probability Maps and the unified segmentation model37 included in the new segment toolbox. The unified segmentation enables bias correction, spatial normalization and tissue segmentation into GM, WM and cerebrospinal fluid within the same generative model. Segmented images were modulated to preserve participant’s original tissue quantity after spatial normalization and resliced to 1.5 × 1.5 × 1.5 mm. Finally, the normalized, segmented and modulated GM images were smoothed using a Gaussian kernel of 10-mm full-width at half-maximum.

Preprocessing of diffusion data

Diffusion data preprocessing was performed using FMRIB Diffusion Toolbox in the FSL software (www.fmrib.ox.ac.uk/fsl) and consisted of affine registration to the first b=0 image for head motion and eddy currents correction, brain extraction using the Brain Extraction Tool (BET),38 and diffusion tensor fitting to obtain FA, MD, axial diffusivity (AD) and radial diffusivity (RD) images. Voxel-wise statistical analysis of the FA data was carried out using tract-based spatial statistics (TBSS) part of FSL.39 All participants’ FA data were aligned into a common space using the nonlinear registration tool FNIRT,40, 41 which uses a b-spline representation of the registration warp field.42 Next, the mean FA image was created and thinned to create a mean FA skeleton, which represents the centers of all tracts common to the group. This skeleton was then thresholded to FA>0.2 to keep only the main tracts. Each adolescent’s aligned FA, MD, AD and RD data were then projected onto this skeleton and the resulting data fed into voxel-wise cross-individual statistics.

Data quality control

Images were reconstructed and visually checked for major artifacts before additional processing. Eleven SBP and 19 control T1-weighted datasets were discarded because of low signal-to-noise ratio, motion artifacts and blink eye artifacts. Twelve SBP and 33 control DTI datasets were discarded because of head movement and poor tensor computation. After voxel-based morphometry preprocessing, one control T1-weighted dataset was discarded because of poor spatial normalization. After TBSS preprocessing, 4 SBP and 16 control DTI datasets were discarded because of poor spatial normalization. Finally, one SBP participant was discarded because of missing PDS assessment. Thus, 30 SBP adolescents and 106 controls were available for voxel-based morphometry, and 25 SBP and 77 controls were available for TBSS analysis.

Statistical analyses

Demographic and clinical data

Analyses were performed using R software (http://www.r-project.org/). Unpaired t-tests, analyses of covariance and χ2 tests were used for comparing continuous and categorical variables in between-group comparisons. Statistical significance was defined as P<0.05, two-tailed.

T1-weighted image analysis

Voxel-wise comparisons were performed using SPM8 within the framework of the general linear model. GM results were obtained within a frontal-limbic mask adapted from Nugent and co-workers,43 drawn using Wake Forest University PickAtlas,44 and including the cingulate, orbital, medial prefrontal and superior temporal cortices. Brain locations were reported as x, y and z coordinates in Montreal Neurologic Institute space. Between-group comparisons of T1-weighted images were performed on GM and WM tissues using a two-sample t-test with PDS, scanner type and total intracranial volume (TIV) as confounding covariates. Significance threshold was set at P<0.05, family-wise error corrected for multiple comparisons.

DTI data analysis

Voxel-wise group comparisons on FA, RD, AD and MD maps were tested in the framework of the general linear model using a randomization-based method (5000 permutations). We included PDS and DTI acquisition type as confounding covariates. Statistical thresholds were set at P<0.05 family-wise error corrected and threshold-free cluster enhancement corrected. The ICBM-DTI-81 WM labels atlas45 and the JHU tractography atlas46 were used within FSLView to locate the tracts that displayed significant differences between groups.

Results

The SBP participants had higher scores than the controls regarding pubertal status and all difficulties scores retrieved using the SDQ, particularly the total impact scale score that reflects overall social impairment.47 No other clinical difference was found between groups (Table 1).

Voxel-based morphometry analysis

No between-group difference was found in TIV and WM/TIV, but GM/TIV was significantly lower, and cerebrospinal fluid/TIV higher, in SBP (vs control) adolescents. On a regional basis, GM volume was significantly smaller in the SBP adolescents in the left anterior cingulate (ACC) (Table 2 and Figure 1).

Table 2 Differences in gray-matter volumes between subthreshold bipolar and healthy adolescents
Figure 1
figure 1

(a) Top, sagittal and axial views: T-maps of comparisons between subthreshold bipolar (SBP) adolescents and controls superimposed on the average T1-weighted magnetic resonance imaging (MRI) of all IMAGEN database participants. Blue colour indicates white-matter tracts where fractional anisotropy (FA) is lower; red indicates tracts where radial diffusivity (RD) is higher (P<0.05 family-wise error corrected) in SBP adolescents; green indicates regions with lower gray-matter volume in SBP adolescents (P<0.001, uncorrected); white indicates white-matter skeleton. FA and RD images are displayed using the ‘tbss_fill’ script, which allows better visualization of the regions with significant between-group differences. (b) Three-dimensional representation using the Anatomist software (http://brainvisa.info); statistical maps are projected onto a single IMAGEN participant brain mesh. The same color code as above is in use.

PowerPoint slide

There were no regions of larger GM volume in SBP adolescents. Voxel-based morphometry of WM images did not retrieve any significant cluster in any comparison.

TBSS analysis

Within the WM skeleton, SBP adolescents (vs controls) had lower FA values globally, while no between-group differences were found for MD, AD and RD global measures (Table 1).

Regionally, the SBP adolescents (vs controls) showed significantly lower FA in a large contiguous cluster (26 872 voxels, peak voxel at x=19, y=−51, z=26 Montreal Neurologic Institute coordinates) encompassing association tracts including the cingulum bundle, bilateral superior and inferior longitudinal, inferior fronto-occipital and uncinate fasciculi; projection fibers including the corticospinal tract, and anterior thalamic radiations; and commissural fibers including the genu, body and splenium of the corpus callosum (Table 3 and Figure 1).

Table 3 Differences in FA and RD between SBP and healthy adolescents

The SBP adolescents (vs controls) had significantly higher RD in two clusters (cluster 1: 1119 voxels, peak voxel at x=36, y=−56, z=27; cluster 2: 145 voxels, peak voxel at x=−32, y=−66, z=23 ) including posterior parts of the right cingulum, bilateral superior and inferior longitudinal fasciculi, bilateral inferior fronto-occipital fasciculi, anterior thalamic radiations and the forceps major (Table 3 and Figure 1).

There were no regions of significantly higher FA or lower RD in the SBP adolescents, and no regional between-groups difference was found in MD and AD.

In a slice-wise manner (Figure 2), the FA and RD profiles showed the distribution of the mean values slice-by-slice across the posterior-to-anterior axis of the brain for each group. On this one-dimensional summary, areas of lower FA and higher RD in SBP adolescents (vs controls) were specifically located in the posterior regions at the coronal slice coordinate y=−60. This coronal plane includes the posterior cingulate/precuneus region.

Figure 2
figure 2

Plots of mean fractional anisotropy (FA) (top) and mean radial diffusivity (RD) (bottom) along the posterior to anterior y axis (mm) of the white-matter skeleton. Red dots indicate values in subthreshold bipolar adolescents; blue dots indicate values in the controls; green bars indicate P-values; and green horizontal line the P<0.05 statistical threshold. Boxes indicate zoomed plots between slices 56 and 72.

PowerPoint slide

Discussion

This is the first study to have assessed WM microstructure and GM volumes in community adolescents with SBP. Our findings suggest widespread alterations in WM involving a number of tracts that continue to mature during adolescence. In addition, we observed a smaller left ACC GM volume in the SBP (vs control) adolescents.

Although most of SBP participants were not at familial risk for BD, it is remarkable that lower regional FA was detected in most of the tracts where decreased FA has been reported in previous studies of pediatric BD or depressed adolescent patients. The corpus callosum,14, 15, 48 the uncinate fasciculus,49 the cingulum bundle,14, 15 the inferior longitudinal fasciculus,15 the thalamic radiations, corticospinal tracts and the superior longitudinal fasciculus14 have all previously been reported to be altered in adolescents and also adults with BD.50

Several WM bundles with significant WM alterations in this study are involved in the circuitry implicated in emotional regulation. The cingulum bundle, for instance, a component of the circuitry subserving emotion, originates in various parts of the cingulate cortex, such as the rostral cingulate, and caudal cingulate, and projects rostrally to various areas implicated in mood disorders,51 including the DLPFC, the ACC and orbital–frontal cortex, and caudally to the ventral temporal lobe, retrosplenial cortex and parietal cortex.52 Also, WM changes were particularly marked in all parts of the corpus callosum. Most of the prefrontal fibers course through this bundle, notably the fibers originating in DLPFC and ACC, which course through the rostral part of the genu. Thus, WM alterations in those bundles might confer vulnerability to emotional dysregulation.

Among the other bundles with WM integrity alterations in the SBP group, the uncinate fasciculus is worth noting as it links the anterior temporal lobe and amygdala with the medial and orbital prefrontal cortices, and as such appears to be implicated in the regulation of behavior expression and emotional responses subserved by those structures.53

Adolescence is a critical period for WM and GM maturation, involving widespread reductions in cortical GM volume, increases in WM volume and FA and decreased RD and AD.54, 55, 56, 57, 58 WM maturation involves both myelination and axonal caliber increases, but the maturation mechanisms seem region-dependent and are still debated.56, 59, 60 Considerable variation has been reported in age-related FA changes across tracts, suggesting a hierarchical pattern of maturation in which some connections develop more slowly than others.56 It has been shown that WM development takes place along a posterorostral gradient, and that WM in posterior regions tends to develop earlier and faster.61 Our results indicate that WM integrity alterations in the SBP group were more distinct in the posterior cingulate regions, with lower FA and higher RD in the posterior regions, as denoted by the slice-wise analysis of the mean FA and RD values in each coronal slice along the posterorostral axis. RD decrease has been reported in postpubertal stage;62 therefore, pubertal status of the SBP group cannot account for the higher RD, as the SBP adolescents had higher pubertal status as compared with controls. Thus, the FA and RD changes in this study point out possibly reduced or delayed WM maturation, in the SBP adolescents, who might be lagging ‘behind’ the controls in a number of WM bundles, notably in regions that develop faster than others in early adolescence.

It is noteworthy that the posterior part of the cingulate bundle is adjacent to the PCC/precuneus region encompassing Brodmann areas 23/31, which also showed a trend for lower gray matter volume in the SBP adolescents. The posterior cingulate/precuneus has been identified as a prominent ‘hub’, central to resting connectivity networks for integration of neural processing; the dysfunction of such a region, which might develop with GM changes, has been hypothesized to contribute to psychiatric diseases.63

These WM findings were also concomitant to smaller regional GM volumes in SBP subjects in the ‘affective division’ of the ACC,64 and to a lesser extent in the medial prefrontal cortex, both regions where the cingulum bundle originates and projects. The ACC affective division is involved in assessing the salience of emotional and motivational information and the regulation of emotional responses, and is part of a medial prefrontal network involved in self-referential functions where GM decreases have been reported in early-onset recurrent major depression or BD.51 Of note, the age of peak cortical thickness in the cingulate has been found around 14 years,65 which is precisely the mean age of this sample.

Smaller GM volume in those prefrontal regions, associated with delayed WM maturation in projecting bundles, might contribute to inadequate modulation of subcortical structures, and lead to mood dysregulation, which could underlie vulnerability to BD. Indeed, a recent model postulates that reduced prefrontal modulation of subcortical and medial temporal structures within the anterior limbic network, associated with early WM abnormalities, contributes to further developmental abnormalities involving prefrontal connections with the amygdala, a structure supposed to be at the core of emotional dysregulation in BD.66, 67 Lower amygdala volume is the most replicated neuroanatomical finding in adolescents with BD.68 However, in line with our findings, no change in amygdala volume has been reported in adolescents at onset of BD,69, 70 nor in youth at familial risk for BD.71, 72

Thus, these findings in the SBP adolescents, showing both WM microstructure alterations in the frontal–limbic pathways and lower ACC/medial prefrontal GM volumes, support an altered interplay between WM and GM maturation in the GM regions and WM tracts involved in mood regulation, which might lead to later amygdala development abnormalities.

Limitations

The main limitation is the lack of known outcome of the participants. Indeed, the SBP adolescents had no prior psychiatry history, nor were they at more ‘familial risk’ of affective disorders than the controls. Thus, the SBP adolescents may have expressed extreme personality features in the range of normal mood instability. However, the inclusion symptoms were core symptoms of mania whose rates have been found to be high (82% for expansive mood and pressure of speech; 90% for unstable mood, 57% for sleep decrease) in adolescent SBP, and have proved to have a significant pathological meaning as regards outcomes.30 Furthermore, the SBP adolescents had significantly higher SDQ emotional and behavioral difficulties scores, which has been related with overall social impairment.47 Yet, given the outcomes of subjects with SBP,6, 9 it is difficult to interpret the meaning of the findings with regard to a specific disorder, as those early WM and GM alterations might reflect vulnerability to BD or to anxiety or depression as well. In addition, the risk of developing depression in adolescence has been related to pubertal status,73 which is consistent with the higher pubertal status in the present SBP group, although the groups were matched by sex and age.

The pathophysiological significance of altered directional diffusivities in relation to WM pathology is an ongoing area of research. DTI parameters could reflect several facets of WM, as axonal size, density, coherence and membrane structure, or lipid, protein and macromolecule content, or water compartmentalization. Previous studies have reported that lower FA might indicate decreased fiber organization and coherence, notably when associated with increased RD,74 or lack of myelin,75 and/or changes in other structural components of the axon and WM microenvironment. The absence of axial diffusivity change might reflect preserved axonal integrity, while radial diffusivity might inform on myelin sheaths.75 Thus, reduced FA with increased RD, and no axial diffusivity changes, as in the present SBP group, could indicate altered myelination without axonal loss. However, the histology and molecular changes in oligodendrocytes and myelin microenvironment reflected by DTI alterations cannot be tested in humans by noninvasive techniques.

In summary, our results suggest that SBP in adolescents is associated with structural connectivity alterations likely to affect the development of WM bundles and connected GM areas involved in emotion regulation.