# Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals

Pauline Favre, Melissa Pauling, Jacques Stout, Franz Hozer, Samuel Sarrazin, Christoph Abé, Martin Alda, Clara Alloza, Silvia Alonso-Lana, Ole A Andreassen, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Pauline Favre, Melissa Pauling, Jacques Stout, Franz Hozer, Samuel Sarrazin, et al. Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals. Neuropsychopharmacology, 2019, 44 (13), pp.2285 - 2293. 10.1038/s41386-019-0485-6. hal-05223805

HAL Id: hal-05223805

https://hal.science/hal-05223805v1

Submitted on 26 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for
updates

ARTICLE **OPEN**

# Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals

This article has been corrected since Advance Online Publication and a correction is also printed in this issue

Pauline Favre<sup>1,2</sup>, Melissa Pauling<sup>1,2</sup>, Jacques Stout<sup>1</sup>, Franz Hozer<sup>1,2,3,4</sup>, Samuel Sarrazin<sup>1,2,5,6</sup>, Christoph Abè<sup>7</sup>, Martin Alda<sup>1,8</sup>, Clara Alloza<sup>9,10</sup>, Silvia Alonso-Lana<sup>11,12</sup>, Ole A. Andreassen<sup>13,14</sup>, Bernhard T. Baune<sup>15,16</sup>, Francesco Benedetti<sup>17,18</sup>, Geraldo F. Busatto<sup>19,20</sup>, Erick J. Canales-Rodríguez<sup>11</sup>, Xavier Caseras<sup>21</sup>, Tiffany Moukbel Chaim-Avancini<sup>19,20</sup>, Christopher R. K. Ching<sup>22,23</sup>, Udo Dannlowski<sup>16</sup>, Michael Deppe<sup>24</sup>, Lisa T. Eyler<sup>25,26</sup>, Mar Fatjo-Vilas<sup>27,28</sup>, Sonya F. Foley<sup>29</sup>, Dominik Grotegerd<sup>16</sup>, Tomas Hajek<sup>8</sup>, Unn K. Haukvik<sup>13,14</sup>, Fleur M. Howells<sup>30,31</sup>, Neda Jahanshad<sup>22</sup>, Harald Kugel<sup>32</sup>, Trine V. Lagerberg<sup>13,14</sup>, Stephen M. Lawrie<sup>33</sup>, Julia O. Linke<sup>34,35</sup>, Andrew McIntosh<sup>33,36</sup>, Elisa M. T. Melloni<sup>37</sup>, Philip B. Mitchell<sup>38,39</sup>, Mircea Polosan<sup>40</sup>, Edith Pomarol-Clotet<sup>27,28</sup>, Jonathan Repple<sup>16</sup>, Gloria Roberts<sup>38,39</sup>, Annerine Roos<sup>41</sup>, Pedro G. P. Rosa<sup>19,20</sup>, Raymond Salvador<sup>28,29</sup>, Salvador Sarró<sup>28,29</sup>, Peter R. Schofield<sup>42,43</sup>, Mauricio H. Serpa<sup>17,18</sup>, Kang Sim<sup>44,45,46</sup>, Dan J. Stein<sup>41</sup>, Jess E. Sussmann<sup>33</sup>, Henk S. Temmingh<sup>30,47</sup>, Paul M. Thompson<sup>23</sup>, Norma Verdolini<sup>27,28,48</sup>, Eduard Vieta<sup>28,48</sup>, Michele Wessa<sup>34</sup>, Heather C. Whalley<sup>34</sup>, Marcus V. Zanetti<sup>18,19,49</sup>, Marion Leboyer<sup>1,50,51</sup>, Jean-François Mangin<sup>1</sup>, Chantal Henry<sup>52</sup>, Edouard Duchesnay<sup>1</sup> and Josselin Houenou<sup>1,2,50,51</sup>
for the ENIGMA Bipolar Disorder Working Group

Fronto-limbic white matter (WM) abnormalities are assumed to lie at the heart of the pathophysiology of bipolar disorder (BD); however, diffusion tensor imaging (DTI) studies have reported heterogeneous results and it is not clear how the clinical heterogeneity is related to the observed differences. This study aimed to identify WM abnormalities that differentiate patients with BD from healthy controls (HC) in the largest DTI dataset of patients with BD to date, collected via the ENIGMA network. We gathered individual tensor-derived regional metrics from 26 cohorts leading to a sample size of *N* = 3033 (1482 BD and 1551 HC). Mean fractional anisotropy (FA) from 43 regions of interest (ROI) and average whole-brain FA were entered into univariate mega- and meta-analyses to differentiate patients with BD from HC. Mega-analysis revealed significantly lower FA in patients with BD compared with HC in 29 regions, with the highest effect sizes observed within the corpus callosum (R² = 0.041, *P*<sub>corr</sub> < 0.001) and cingulum (right: R² = 0.041, left: R² = 0.040, *P*<sub>corr</sub> < 0.001). Lithium medication, later onset and short disease duration were related to higher FA along multiple ROIs. Results of the meta-analysis showed similar effects. We demonstrated widespread WM abnormalities in BD and highlighted that altered WM connectivity within the corpus callosum and the cingulum are strongly associated with BD. These brain abnormalities could represent a biomarker for use in the diagnosis of BD. Interactive three-dimensional visualization of the results is available at www.enigma-viewer.org.

Neuropsychopharmacology (2019) 44:2285–2293; https://doi.org/10.1038/s41386-019-0485-6

<footnote n="1">¹Neurospin, CEA, Université Paris-Saclay, Gif-sur-Yvette, France; ²INSERM Unit U955, Team 15, “Translational Psychiatry”, Créteil, France; ³Assistance Publique-Hôpitaux de Paris (AP-HP), Corentin-Celton Hospital, Department of Psychiatry, Issy-les-Moulineaux, France; ⁴Paris Descartes University, PRES Sorbonne Paris Cité, Paris, France; ⁵Pôle de psychiatrie, DHU PePSY, Hôpitaux Universitaires Mondor, Créteil, France; ⁶Bipol Falret, Fondation Falret, St Ouen, France; ⁷Clinical Neuroscience, Karolinska Institute, Stockholm, Sweden; ⁸Department of Psychiatry, Dalhousie University, Halifax, Canada; ⁹Division of Psychiatry, University of Edinburgh, Edinburgh, UK; ¹⁰Department of Child and Adolescent Psychiatry, IISGM, Hospital General Universitario Gregorio Marañón, Madrid, Spain; ¹¹FIDMAG Research Foundation, Barcelona, Spain; ¹²CIBERSAM, Barcelona, Spain; ¹³Department of Mental Health and Addiction, University of Oslo, Oslo, Norway; ¹⁴NORMENT K.G. Jebsen Centre for Psychosis Research, Oslo University Hospital, Oslo, Norway; ¹⁵Department of Psychiatry, University of Melbourne, Melbourne, Australia; ¹⁶Department of Psychiatry, University of Münster, Münster, Germany; ¹⁷Division of Neuroscience, San Raffaele Scientific Institute, Milano, Italy; ¹⁸University Vita-Salute San Raffaele, Milano, Italy; ¹⁹Laboratory of Psychiatric Neuroimaging (LIM-21), Departamento e Instituto de Psiquiatria, Hospital das Clínicas HCFMUSP, Faculdade de Medicina, Universidade de São Paulo, São Paulo, Brazil; ²⁰Centrer for Interdisciplinary Research on Applied Neurosciences (NAPNA), University of São Paulo, São Paulo, Brazil; ²¹MRC Centre for Neuropsychiatric Genetics and Genomics, Cardiff University, Cardiff, UK; ²²Interdepartmental Neuroscience Program, University of California, Los Angeles, CA, USA; ²³Imaging Genetics Center, Mark and Mary Stevens Neuroimaging and Informatics Institute, Keck School of Medicine of USC, University of Southern California, Marina del Rey, Los Angeles, CA, USA; ²⁴University of Münster, Department of Neurology, Münster, Germany; ²⁵Department of Psychiatry, University of California San Diego, La Jolla, CA, USA; ²⁶Mental Illness Research Education and Clinical Center, VA San Diego Healthcare System, La Jolla, CA, USA; ²⁷FIDMAG Research Foundation, Barcelona, Spain; ²⁸CIBERSAM, Barcelona, Spain; ²⁹Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, UK; ³⁰Department of Psychiatry and Mental Health, University of Cape Town, Cape Town, South Africa; ³¹Neuroscience Institute, University of Cape Town, Cape Town, South Africa; ³²University of Muenster, Institute of Clinical Radiology, Münster, Germany; ³³Department of Psychiatry, Royal Edinburgh Hospital, Edinburgh, UK; ³⁴Department of Clinical Psychology and Neuropsychology, Johannes Gutenberg-Universität Mainz, Mainz, Germany; ³⁵Emotion and Development Branch, National Institute of Mental Health, Bethesda, MD, USA; ³⁶Center for Cognitive Ageing and Cognitive Epidemiology, University of Edinburgh, Edinburgh, UK; ³⁷Psychiatry and Clinical Psychobiology, Division of Neuroscience, Scientific Institute and University Vita-Salute San Raffaele, Milano, Italy; ³⁸School of Psychiatry, University of New South Wales, Sydney, Australia; ³⁹Black Dog Institute, Sydney, Sydney, Australia; ⁴⁰Univ. Grenoble Alpes, Inserm, U1216, CHU Grenoble Alpes, Grenoble Institut Neurosciences, Grenoble, France; ⁴¹Dept of Psychiatry, SAMRC Unit on Risk & Resilience in Mental Disorders, University of Cape Town, Cape Town, South Africa; ⁴²Neuroscience Research Australia, Sydney, Australia; ⁴³School of Medical Sciences, University of New South Wales, Sydney, Australia; ⁴⁴West Region, Institute of Mental Health, Singapore, Singapore; ⁴⁵Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore; ⁴⁶Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore, Singapore; ⁴⁷Valkenberg Hospital, Cape Town, South Africa; ⁴⁸Institute of Neuroscience, Hospital Clinic, University of Barcelona, IDIBAPS, Barcelona, Spain; ⁴⁹Hospital Sírio-Libanés, São Paulo, Brazil; ⁵⁰Assistance Publique-Hôpitaux de Paris (AP-HP), CHU Mondor, Psychiatry Department, Créteil, France; ⁵¹Faculté de Médecine, Université Paris Est Créteil, Créteil, France and ⁵²Institut Pasteur, Unité Perception et Mémoire, Paris, France
Correspondence: Pauline Favre (pauline@favre-univ.fr)</footnote>
## INTRODUCTION

Bipolar disorder (BD) is a severe chronic mental illness that affects ~1% of the general population [1]. There is often a long period with inadequate treatment before the diagnosis is established [2]. Consequently, there is a great need to identify biomarkers of BD. A better understanding of the neurobiology of BD could ultimately help to refine the diagnosis and guide innovative interventions. Recent advances in magnetic resonance imaging (MRI) could help to achieve this goal.

Neural models of BD suggest a role of fronto-limbic dysconnectivity in the emergence of mood symptoms of BD [3, 4]. This model is mainly supported by results from functional MRI (fMRI) studies demonstrating that emotional instability in this disorder might be underpinned by abnormal connectivity between frontal and limbic regions [5, 6]. However, results from diffusion tensor imaging (DTI) studies, a technique that allows the exploration of structural connectivity in vivo, have highlighted far more extensive brain abnormalities in BD. Indeed, the first DTI studies identified alterations in limbic tracts [7–9], followed by numerous studies that reported WM alterations within non-limbic regions, such as the corpus callosum [10–15] and corona radiata [16]. Meta-analyses based on whole-brain data have revealed lower fractional anisotropy (FA), a metric derived from DTI known to be positively correlated with the directionality and coherence of white matter bundles [17], in patients with BD near the parahippocampal gyrus, subgenual cingulate cortex [18], temporo-parietal junction and cingulum [19].

Inconsistencies in the location of WM microstructure alterations may be related to limited sample sizes and diversity in methods to collect data from different populations and for DTI data analysis. Indeed, differences in sample characteristics such as age of onset, disease duration, psychotic features, and lithium treatment, all of which have been associated with WM features [12, 20–22], may have contributed to the inconsistency in previous findings. Consequently, large harmonized multi-center studies are required to improve the reliability of case-control findings.

The ENIGMA consortium presents a framework to identify generalizable biomarkers, by analyzing large samples with a harmonized processing pipeline—a strategy that has already identified widespread cortical alterations and specific subcortical volumetric abnormalities in patients with BD [23, 24]. Thus, we analyzed DTI data from the ENIGMA-BD working group with the objectives of (i) identifying reliable generalizable WM abnormalities in BD using mega- and meta- analytics; (ii) testing if clinical characteristics modulate WM microstructure using mega- analytics. Specifically, we expected more pronounced alterations (i.e., larger FA differences with respect to healthy controls) in WM microstructure in patients with a more severe course of illness, and a significant association with psychotropic medication.

## METHODS

### Samples

The ENIGMA-BD DTI working group, comprised of 26 cohorts spanning 12 countries, yielded a total of 3033 individuals (1551 healthy controls (HC) and 1482 patients with BD) included in this study. Demographic and clinical information from the whole sample is shown in Table 1; details of the contributing sites may be found in Table S1 and available clinical data for each site is provided in Table S2. Each cohort comprised a minimum of 12 subjects per group and a minimal ratio of patients to controls of 1:3, to allow for robust comparisons and meta-analysis. When needed, we randomly removed some subjects from a given group (mainly control subjects that were too numerous at 4 sites, except for one site that comprised too many patients in comparison to controls; for details, see Table S3). The current analysis includes data acquired until February 2018.

All participating sites obtained approval from their local ethics committees and all participants gave written informed consent. Participants younger than 18 or older than 65 as well as individuals with diffusion images with low quality after visual inspection (e.g., movement artifacts) were excluded from the analyses.

### Image processing

Acquisition parameters for each of the 26 sites are provided in Table S4. The pre-processing (i.e., eddy current and echo-planar corrections and tensor fitting) was performed at each site using harmonized analysis and quality control protocols from the ENIGMA consortium that have previously been applied in large-scale studies of schizophrenia [25]; recommended pipelines and procedures for the image analyses and quality control are provided online at the ENIGMA-DTI website (http://enigma.ini.usc.edu/protocols/dti-protocols/). After estimation of tensors, each site performed the image analysis and extracted the FA of each region of interest (ROI) (see description in Table S5) according to the ENIGMA-DTI protocol. The multi-subject JHU white matter parcellation atlas [26] was used to parcellate regions of interest from the ENIGMA template in MNI space. Mean FA from 43 regions of interest (ROI) as well as average whole-brain FA were then extracted for each participant across all cohorts.

### Mega-analysis

Our first aim was to identify WM microstructure differences between patients with BD and HC. We merged individual FA values of the 43 ROIs and Average FA (from each cohort) into one mega-analysis and entered them separately in a linear mixed model (using R software version 3.2.1. (R Core Team, 2015) and lme4 package [27]) including fixed effects for the diagnosis (patients vs. controls), age, sex, and random intercepts for each site:

$$FA \text{ ROI}_i = \text{Intercept} + \beta_1 * \text{Diagnosis} + \beta_2 * \text{Age} + \beta_3 * \text{Sex} + \text{random effect/site}$$

We used Bonferroni correction to control for multiple comparisons ($p < 0.05/44 = 0.0011$). We also assessed the influence of average FA (per subject) across the entire TBSS FA tract skeleton (including core and periphery FA [25]) on local FA differences observed in the first analysis by running the same models including average FA as a covariate.

We performed additional analyses to assess how age, sex, illness duration, age of onset, medication at the time of scan (lithium, antipsychotics, anticonvulsants, and antidepressants), illness severity, history of psychotic symptoms and type of BD (type I vs. type II) might have modulated the main effect of diagnosis. We tested the effect of age and sex by including age-by-diagnosis and sex-by-diagnosis interaction terms. We included medication and history of psychosis as dichotomous measures in the analyses (yes/no variables) and used the density of episodes as an index of illness severity (number of mood episodes/illness duration). Importantly, each analysis controlled for age and sex, so that associations with illness duration and the age of onset would not be confounded by global age differences.

Age, sex, and diagnosis were available for all participants, whereas the remaining variables were available for some sites only (see Table S2 for details of available data for each site).

### Meta-analysis

Given previous demonstrations of the usefulness of meta-analysis for multisite neuroimaging [28], we performed a meta-analysis to allow comparisons with previous ENIGMA studies and comparison across sites. Similarly to previous ENIGMA meta-analyses, we conducted a random-effects inverse-variance weighted meta-analysis (R, metaphore package), to combine Cohen's d effect size
<table><caption>**Table 1.** Descriptive statistics of sample</caption><thead><tr><th rowspan="2"></th><th colspan="3">Bipolar disorder</th><th colspan="3">Healthy controls</th><th rowspan="2">t/χ<sup>2</sup></th><th rowspan="2">p-value</th></tr><tr><th>N<sup>a</sup></th><th>mean/freq<sup>b</sup></th><th>std</th><th>N<sup>a</sup></th><th>mean/freq<sup>b</sup></th><th>std</th></tr></thead><tbody><tr><td>Mean age</td><td>1482</td><td>39.60</td><td>12.15</td><td>1551</td><td>35.08</td><td>12.10</td><td>10.11</td><td>&lt;0.001</td></tr><tr><td>Sex (% females)</td><td>1482</td><td>60.66%</td><td></td><td>1551</td><td>51.13%</td><td></td><td>25.77</td><td>&lt;0.001</td></tr><tr><td>Age of onset</td><td>1046</td><td>25.19</td><td>10.35</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Illness duration</td><td>1045</td><td>15.47</td><td>10.58</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Depression Score</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HDRS-17</td><td>115</td><td>8.63</td><td>8.29</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HDRS-21</td><td>285</td><td>6.42</td><td>8.06</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MADRS</td><td>230</td><td>9.45</td><td>9.79</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Number of depressive episodes</td><td>587</td><td>5.84</td><td>6.02</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mania Score (YMRS)</td><td>545</td><td>2.80</td><td>3.93</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Number of manic episodes</td><td>485</td><td>4.30</td><td>5.44</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Total number of major episodes</td><td>476</td><td>10.45</td><td>10.15</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Density of episodes<sup>c</sup></td><td>414</td><td>0.82</td><td>1.00</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>On medication</td><td>904</td><td>84.62%</td><td>36.09%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Antipsychotics</td><td>862</td><td>45.82%</td><td>49.85%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Antidepressants</td><td>903</td><td>31.34%</td><td>46.41%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anticonvulsants</td><td>812</td><td>39.29%</td><td>48.87%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lithium</td><td>824</td><td>42.48%</td><td>49.46%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>History of psychotic symptoms</td><td>908</td><td>54.30%</td><td>49.84%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lifetime alcohol abuse</td><td>272</td><td>8.09%</td><td>27.32%</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Onset time<sup>d</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Early</td><td>334</td><td>31.93%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Intermediate</td><td>534</td><td>51.05%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Late</td><td>178</td><td>17.02%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BD type</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BD1</td><td>432</td><td>77.01%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BD2</td><td>129</td><td>22.99%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mood phase</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Depressed</td><td>313</td><td>46.23%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Euthymic</td><td>336</td><td>49.63%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Hypomanic</td><td>2</td><td>0.30%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Manic</td><td>18</td><td>2.66%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Mixed</td><td>8</td><td>1.18%</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

of each of the 26 cohorts of the study, both for right and left tracts separately and for bilateral tracts (to allow comparison with other ENIGMA DTI working groups). We calculated the *I*² statistic to estimate the heterogeneity of the diagnostic effects across sites. This analysis was run following publicly available scripts on the ENIGMA-GitHub (https://github.com/ENIGMA-git).

## RESULTS

We included 1482 patients with BD and 1551 HC. The patients were significantly older than the controls (mean age BD = 39.6 years; mean age HC = 35.1 years; *t* = 10.11; *p* < 0.001) and comprised a higher proportion of females (60.7 vs. 51.1%; *χ*² = 25.77; *p* < 0.001). We included both age and sex as covariates in the mega- and meta-analyses, and tested for the age-by-diagnosis and sex-by-diagnosis interactions for further exploration of these effects.

## Mega-analysis

Linear mixed models revealed significantly lower FA in BD vs. HC along 29 out of 43 WM tracts and whole skeleton FA (see Table 2, Fig. 1). The largest effect sizes were found in the whole corpus callosum (CC) (*R*² = 0.0441; *P* < 1.0 × 10⁻²⁰), followed by the body (*R*² = 0.0368; *P* < 1.0 × 10⁻²⁰) and genu (*R*² = 0.0331; *P* < 1.0 × 10⁻²⁰) of the CC and the bilateral cinguli (right: *R*² = 0.0281; *P* < 1.0 × 10⁻²⁰; left: *R*² = 0.0269; *P* < 1.0 × 10⁻²⁰). Notably, we found lower FA in bilateral tracts, with the exception of the inferior fronto-occipital fasciculus, where significant difference was observed only in the right hemisphere. In a second analysis, with similar LMM but also covarying for average FA, we still observed lower FA in BD vs. HC across 19 tracts, meaning that the whole-brain average FA moderately influenced the results and that the effects were not exclusively driven by a global decrease in FA in patients (Table S6).

<footnote><sup>a</sup>Number of available data</footnote>

<footnote><sup>b</sup>Proportion calculated among the available data</footnote>

<footnote><sup>c</sup>Density = total episodes/illness duration</footnote>

<footnote><sup>d</sup>Early: < 18 years; 18 years < intermediate < 35 years; Late: > 35 years</footnote>
<table><caption>**Table 2.** Mega-analysis results: linear mixed model parameters sorted by effect size (descending order) for FA differences between bipolar patients and healthy controls after controlling for age and sex</caption><thead><tr><th>ROI</th><th>β</th><th>s.e.</th><th>t-value</th><th>P<sub>corr</sub> > |t|</th><th>R<sup>2</sup></th><th>[0.025</th><th>0.975]</th><th>Sign.</th></tr></thead><tbody><tr><td colspan="9"><strong>Projection fibers</strong></td></tr><tr><td>PTR.R</td><td>0.0098</td><td>0.0013</td><td>7.3542</td><td>1.09E-11</td><td>0.0176</td><td>0.0095</td><td>0.0281</td><td>***</td></tr><tr><td>PTR.L</td><td>0.0079</td><td>0.0013</td><td>5.9150</td><td>1.63E-07</td><td>0.0115</td><td>0.0051</td><td>0.0203</td><td>***</td></tr><tr><td>ACR.L</td><td>0.0065</td><td>0.0011</td><td>5.8177</td><td>2.91E-07</td><td>0.0110</td><td>0.0048</td><td>0.0197</td><td>***</td></tr><tr><td>CR.L</td><td>0.0046</td><td>0.0009</td><td>5.2310</td><td>7.94E-06</td><td>0.0088</td><td>0.0034</td><td>0.0168</td><td>***</td></tr><tr><td>ACR.R</td><td>0.0056</td><td>0.0011</td><td>4.9348</td><td>3.73E-05</td><td>0.0080</td><td>0.0029</td><td>0.0156</td><td>***</td></tr><tr><td>CR.R</td><td>0.0040</td><td>0.0009</td><td>4.6037</td><td>1.90E-04</td><td>0.0068</td><td>0.0022</td><td>0.0140</td><td>***</td></tr><tr><td>PCR.R</td><td>0.0040</td><td>0.0011</td><td>3.8079</td><td>6.29E-03</td><td>0.0048</td><td>0.0011</td><td>0.0110</td><td>**</td></tr><tr><td>ALIC.L</td><td>0.0040</td><td>0.0011</td><td>3.7958</td><td>6.61E-03</td><td>0.0047</td><td>0.0011</td><td>0.0108</td><td>**</td></tr><tr><td>ALIC.R</td><td>0.0038</td><td>0.0010</td><td>3.6899</td><td>1.00E-02</td><td>0.0044</td><td>0.0009</td><td>0.0104</td><td>*</td></tr><tr><td>PCR.L</td><td>0.0036</td><td>0.0010</td><td>3.4692</td><td>2.33E-02</td><td>0.0040</td><td>0.0007</td><td>0.0098</td><td>*</td></tr><tr><td>SCR.L</td><td>0.0028</td><td>0.0010</td><td>2.7731</td><td>2.46E-01</td><td>0.0026</td><td>0.0002</td><td>0.0075</td><td>NS</td></tr><tr><td>SCR.R</td><td>0.0022</td><td>0.0010</td><td>2.3003</td><td>9.46E-01</td><td>0.0017</td><td>0.0000</td><td>0.0060</td><td>NS</td></tr><tr><td>IC.L</td><td>0.0014</td><td>0.0008</td><td>1.6965</td><td>1.00E+00</td><td>0.0009</td><td>0.0000</td><td>0.0044</td><td>NS</td></tr><tr><td>RLIC.L</td><td>0.0014</td><td>0.0011</td><td>1.3106</td><td>1.00E+00</td><td>0.0006</td><td>0.0000</td><td>0.0036</td><td>NS</td></tr><tr><td>IC.R</td><td>0.0011</td><td>0.0008</td><td>1.3451</td><td>1.00E+00</td><td>0.0006</td><td>0.0000</td><td>0.0036</td><td>NS</td></tr><tr><td>CST.R</td><td>-0.0021</td><td>0.0016</td><td>-1.2829</td><td>1.00E+00</td><td>0.0005</td><td>0.0000</td><td>0.0035</td><td>NS</td></tr><tr><td>CST.L</td><td>-0.0017</td><td>0.0017</td><td>-1.0213</td><td>1.00E+00</td><td>0.0003</td><td>0.0000</td><td>0.0030</td><td>NS</td></tr><tr><td>RLIC.R</td><td>0.0009</td><td>0.0012</td><td>0.7712</td><td>1.00E+00</td><td>0.0002</td><td>0.0000</td><td>0.0025</td><td>NS</td></tr><tr><td>PLIC.L</td><td>-0.0007</td><td>0.0010</td><td>-0.7363</td><td>1.00E+00</td><td>0.0002</td><td>0.0000</td><td>0.0024</td><td>NS</td></tr><tr><td>PLIC.R</td><td>-0.0002</td><td>0.0010</td><td>-0.2474</td><td>1.00E+00</td><td>0.0000</td><td>0.0000</td><td>0.0018</td><td>NS</td></tr><tr><td colspan="9"><strong>Association fibers</strong></td></tr><tr><td>CGC.R</td><td>0.0136</td><td>0.0015</td><td>9.3757</td><td>0.00E+00</td><td>0.0281</td><td>0.0176</td><td>0.0408</td><td>***</td></tr><tr><td>CGC.L</td><td>0.0138</td><td>0.0015</td><td>9.1811</td><td>0.00E+00</td><td>0.0269</td><td>0.0166</td><td>0.0395</td><td>***</td></tr><tr><td>EC.L</td><td>0.0057</td><td>0.0009</td><td>6.0965</td><td>5.39E-08</td><td>0.0119</td><td>0.0054</td><td>0.0209</td><td>***</td></tr><tr><td>EC.R</td><td>0.0051</td><td>0.0009</td><td>5.6114</td><td>9.65E-07</td><td>0.0100</td><td>0.0042</td><td>0.0184</td><td>***</td></tr><tr><td>UNC.R</td><td>0.0103</td><td>0.0019</td><td>5.3636</td><td>3.87E-06</td><td>0.0096</td><td>0.0039</td><td>0.0178</td><td>***</td></tr><tr><td>UNC.L</td><td>0.0103</td><td>0.0020</td><td>5.1902</td><td>9.87E-06</td><td>0.0090</td><td>0.0035</td><td>0.0171</td><td>***</td></tr><tr><td>SS.L</td><td>0.0054</td><td>0.0012</td><td>4.6913</td><td>1.25E-04</td><td>0.0072</td><td>0.0024</td><td>0.0146</td><td>***</td></tr><tr><td>IFO.R</td><td>0.0080</td><td>0.0017</td><td>4.6490</td><td>1.53E-04</td><td>0.0071</td><td>0.0024</td><td>0.0144</td><td>***</td></tr><tr><td>SFO.R</td><td>0.0057</td><td>0.0013</td><td>4.2623</td><td>9.18E-04</td><td>0.0060</td><td>0.0017</td><td>0.0128</td><td>***</td></tr><tr><td>SFO.L</td><td>0.0053</td><td>0.0014</td><td>3.8511</td><td>5.29E-03</td><td>0.0049</td><td>0.0012</td><td>0.0112</td><td>**</td></tr><tr><td>FX.ST.R</td><td>0.0047</td><td>0.0014</td><td>3.3251</td><td>3.94E-02</td><td>0.0036</td><td>0.0006</td><td>0.0092</td><td>*</td></tr><tr><td>IFO.L</td><td>0.0059</td><td>0.0018</td><td>3.2375</td><td>5.36E-02</td><td>0.0035</td><td>0.0005</td><td>0.0090</td><td>NS</td></tr><tr><td>SS.R</td><td>0.0039</td><td>0.0012</td><td>3.1308</td><td>7.75E-02</td><td>0.0032</td><td>0.0004</td><td>0.0086</td><td>NS</td></tr><tr><td>FX.ST.L</td><td>0.0038</td><td>0.0013</td><td>2.9044</td><td>1.63E-01</td><td>0.0028</td><td>0.0003</td><td>0.0078</td><td>NS</td></tr><tr><td>CGH.R</td><td>0.0038</td><td>0.0019</td><td>2.0100</td><td>1.00E+00</td><td>0.0013</td><td>0.0000</td><td>0.0052</td><td>NS</td></tr><tr><td>CGH.L</td><td>0.0007</td><td>0.0018</td><td>0.3830</td><td>1.00E+00</td><td>0.0000</td><td>0.0000</td><td>0.0019</td><td>NS</td></tr><tr><td colspan="9"><strong>Commissural fibers</strong></td></tr><tr><td>CC</td><td>0.0123</td><td>0.0010</td><td>11.9431</td><td>0.00E+00</td><td>0.0441</td><td>0.0308</td><td>0.0594</td><td>***</td></tr><tr><td>BCC</td><td>0.0150</td><td>0.0014</td><td>10.7856</td><td>0.00E+00</td><td>0.0368</td><td>0.0247</td><td>0.0511</td><td>***</td></tr><tr><td>GCC</td><td>0.0123</td><td>0.0012</td><td>10.2936</td><td>0.00E+00</td><td>0.0331</td><td>0.0217</td><td>0.0468</td><td>***</td></tr><tr><td>SCC</td><td>0.0077</td><td>0.0009</td><td>8.1690</td><td>1.95E-14</td><td>0.0209</td><td>0.0119</td><td>0.0322</td><td>***</td></tr><tr><td>FX</td><td>0.0185</td><td>0.0025</td><td>7.4763</td><td>4.42E-12</td><td>0.0184</td><td>0.0101</td><td>0.0292</td><td>***</td></tr><tr><td>AverageFA</td><td>0.0025</td><td>0.0006</td><td>4.0044</td><td>2.80E-03</td><td>0.0050</td><td>0.0012</td><td>0.0113</td><td>**</td></tr></tbody></table>

ns not significant

*p<sub>corr</sub> < 0.05; **p<sub>corr</sub> < 0.01; ***p<sub>corr</sub> < 0.001

**Age and sex effects.** To examine differential effects of age and sex on group differences in FA values, we tested for age-by-diagnosis and sex-by-diagnosis interactions for each ROI. Results showed significant age-by-diagnosis interactions in bilateral superior corona radiata, the posterior limb of the internal capsule and left

cingulum, such that there was steeper apparent age-related decline in the HC than BD group in all but the cingulate gyrus portion of the cingulum, where the opposite was found (Table S7; Figure S1). We did not find any significant sex-by-diagnosis interaction (Table S8).
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Results of the mega-analysis. **a** Effect sizes of fractional anisotropy (FA) differences between patients with bipolar disorder (BD) and healthy controls projected on the 43 white matter (WM) tracts analyzed. **b** R squared (effect size) with confidence interval, sorted in increasing order of magnitude, for the regions showing significant differences between bipolar patients and healthy controls</figcaption></figure>

*Effects of clinical variables.* Within the BD group, we found a significant positive relationship of age at onset to FA in the right inferior fronto-occipital fasciculus (Table S9) and a negative association between illness duration and FA within the left cingulum (Table S10) (Fig. S2). In addition, we observed significantly lower FA in patients receiving vs. not receiving antipsychotics within the genu of the CC and in patients receiving vs. not receiving anticonvulsants within multiple ROIs (Figs. S3 and S4; Tables S11 and S12). In contrast, we found higher FA values in several regions among patients receiving vs. not receiving lithium (Fig. S5, Table S13).

We did not observe any significant relationships between FA and antidepressant medication, illness severity, history of psycho-tic symptoms, or BD subtype (I or II) (see Tables S14–S17).

Meta-analysis

Results from the meta-analysis revealed lower FA among 23 out of the 44 ROIs (43 tracts and the whole-brain skeleton) analyzed (Table 3, Fig. 2). Similarly to the mega-analysis, the results showed largest effect sizes for the whole CC ($d = -0.46$; $P = 7.86 \times 10^{-12}$), body of the CC ($d = -0.43$; $P = 5.41 \times 10^{-11}$), and left cingulum ($d = -0.39$; $P = 2.38 \times 10^{-8}$). Overall, the meta-analysis showed similar effects to the mega-analysis but was slightly less sensitive. The $t^2$ test indicates small to high heterogeneity across sites for all effect sizes ($t^2 = 0.002–69.24$). To allow comparison with other DTI studies of the ENIGMA consortium, we also conducted a meta-analysis based on bilateral tracts (i.e., 25 ROIs). We found significant decrease FA in patients with BD compared to HC along 15 fasciculi. Similarly, the higher effect sizes were observed

for the CC ($d = -0.46$; $P = 7.86 \times 10^{-12}$) and cingulum ($d = -0.39$; $P = 4.58 \times 10^{-8}$) (Figure S6, Table S18).

DISCUSSION

In the largest multi-center DTI study of BD to date, we found alterations of WM microstructure in patients with BD along multiple bundles, with strongest effects within the CC and the cingulum. FA was lower in patients in most ROIs, although effect sizes were small. Age, age of onset, illness duration as well as anticonvulsants and antipsychotic medications were associated with lower FA.

We collected individual data from 1482 patients and 1551 controls across 26 international cohorts, allowing a sample size considerably exceeding all prior DTI studies of BD. Unlike most studies that found localized WM alterations in BD, we identified widespread abnormalities (lower FA along 29 out of the 44 regions analyzed in the mega-analysis and 32 out of 44 ROIs in the meta-analysis). Similarly to results in the ENIGMA DTI schizophrenia project, this suggests a global profile of microstructural abnormalities in BD, which are however not specific to that disorder [25].

For both analyses (i.e., mega and meta), the largest effect sizes were observed within the CC and cingulum. This is consistent with a recent meta-analysis showing decreased FA within the CC, cingulum and the anterior superior longitudinal fasciculus in BD in comparison to controls [29]. The cingulum is a major pathway in the limbic system. Impairment of cingulum and uncinate structural integrity is in good agreement with previous models of altered fronto-limbic connectivity in BD [3, 30].
<table><caption>**Table 3.** Meta-analysis results: Cohen's *d* values, their *s.e.*, P-values and *I*² values (heterogeneity between sites) sorted by effect size (descending order) for FA differences between patients with bipolar disorder and healthy controls after controlling for age and sex</caption><thead><tr><th>ROI</th><th>Cohen's d</th><th>s.e.</th><th>p-value</th><th>I²</th><th>p-value (corr)</th><th>Sign.</th></tr></thead><tbody><tr><td colspan="7"><strong>Projection fibers</strong></td></tr><tr><td>ACR.L</td><td>-0.245</td><td>0.048</td><td>3.14E-07</td><td>25.562</td><td>8.86E-06</td><td>***</td></tr><tr><td>ACR.R</td><td>-0.217</td><td>0.051</td><td>2.07E-05</td><td>33.382</td><td>1.03E-03</td><td>**</td></tr><tr><td>CR.L</td><td>-0.202</td><td>0.053</td><td>1.32E-04</td><td>37.861</td><td>4.00E-03</td><td>**</td></tr><tr><td>CR.R</td><td>-0.180</td><td>0.056</td><td>1.26E-03</td><td>43.697</td><td>2.60E-02</td><td>*</td></tr><tr><td>ALIC.L</td><td>-0.158</td><td>0.056</td><td>4.41E-03</td><td>43.311</td><td>2.32E-01</td><td>NS</td></tr><tr><td>PCR.R</td><td>-0.152</td><td>0.043</td><td>4.26E-04</td><td>12.031</td><td>1.13E-01</td><td>NS</td></tr><tr><td>PCR.L</td><td>-0.136</td><td>0.055</td><td>1.37E-02</td><td>42.736</td><td>7.39E-01</td><td>NS</td></tr><tr><td>ALIC.R</td><td>-0.131</td><td>0.048</td><td>6.52E-03</td><td>26.872</td><td>2.78E-01</td><td>NS</td></tr><tr><td>SCR.L</td><td>-0.095</td><td>0.052</td><td>6.90E-02</td><td>36.840</td><td>1.00E+00</td><td>NS</td></tr><tr><td>SCR.R</td><td>-0.072</td><td>0.061</td><td>2.36E-01</td><td>53.242</td><td>1.00E+00</td><td>NS</td></tr><tr><td>IC.L</td><td>-0.070</td><td>0.055</td><td>2.01E-01</td><td>41.932</td><td>1.00E+00</td><td>NS</td></tr><tr><td>IC.R</td><td>-0.058</td><td>0.054</td><td>2.75E-01</td><td>39.437</td><td>1.00E+00</td><td>NS</td></tr><tr><td>RLIC.R</td><td>-0.047</td><td>0.055</td><td>3.90E-01</td><td>41.619</td><td>1.00E+00</td><td>NS</td></tr><tr><td>RLIC.L</td><td>-0.044</td><td>0.055</td><td>4.21E-01</td><td>41.598</td><td>1.00E+00</td><td>NS</td></tr><tr><td>CST.L</td><td>-0.012</td><td>0.065</td><td>8.49E-01</td><td>58.035</td><td>1.00E+00</td><td>NS</td></tr><tr><td>CST.R</td><td>0.012</td><td>0.062</td><td>8.50E-01</td><td>53.890</td><td>1.00E+00</td><td>NS</td></tr><tr><td>PLIC.L</td><td>0.031</td><td>0.053</td><td>5.60E-01</td><td>37.920</td><td>1.00E+00</td><td>NS</td></tr><tr><td>PLIC.R</td><td>0.034</td><td>0.044</td><td>4.35E-01</td><td>13.722</td><td>1.00E+00</td><td>NS</td></tr><tr><td colspan="7"><strong>Association fibers</strong></td></tr><tr><td>CGC.L</td><td>-0.391</td><td>0.062</td><td>2.99E-10</td><td>53.489</td><td>2.38E-08</td><td>***</td></tr><tr><td>CGC.R</td><td>-0.350</td><td>0.057</td><td>7.50E-10</td><td>45.173</td><td>2.00E-06</td><td>***</td></tr><tr><td>EC.L</td><td>-0.233</td><td>0.049</td><td>1.65E-06</td><td>27.203</td><td>1.06E-04</td><td>***</td></tr><tr><td>UNC.L</td><td>-0.231</td><td>0.052</td><td>8.16E-06</td><td>35.172</td><td>3.30E-04</td><td>***</td></tr><tr><td>SS.L</td><td>-0.220</td><td>0.044</td><td>6.96E-07</td><td>15.563</td><td>5.59E-05</td><td>***</td></tr><tr><td>UNC.R</td><td>-0.219</td><td>0.049</td><td>8.05E-06</td><td>28.252</td><td>1.30E-03</td><td>**</td></tr><tr><td>EC.R</td><td>-0.205</td><td>0.042</td><td>1.09E-06</td><td>8.287</td><td>2.58E-04</td><td>***</td></tr><tr><td>IFO.R</td><td>-0.180</td><td>0.047</td><td>1.19E-04</td><td>22.441</td><td>3.99E-03</td><td>**</td></tr><tr><td>IFO.L</td><td>-0.145</td><td>0.039</td><td>2.24E-04</td><td>0.000</td><td>1.41E-02</td><td>*</td></tr><tr><td>FX.ST.R</td><td>-0.145</td><td>0.045</td><td>1.43E-03</td><td>18.610</td><td>9.68E-02</td><td>NS</td></tr><tr><td>SFO.L</td><td>-0.144</td><td>0.051</td><td>4.67E-03</td><td>33.405</td><td>2.00E-01</td><td>NS</td></tr><tr><td>SS.R</td><td>-0.140</td><td>0.053</td><td>8.63E-03</td><td>39.037</td><td>3.17E-01</td><td>NS</td></tr><tr><td>FX.ST.L</td><td>-0.134</td><td>0.039</td><td>6.40E-04</td><td>0.002</td><td>5.12E-02</td><td>NS</td></tr><tr><td>SFO.R</td><td>-0.127</td><td>0.053</td><td>1.65E-02</td><td>38.320</td><td>7.26E-01</td><td>NS</td></tr><tr><td>CGH.R</td><td>-0.080</td><td>0.045</td><td>7.69E-02</td><td>18.794</td><td>1.00E+00</td><td>NS</td></tr><tr><td>CGH.L</td><td>-0.039</td><td>0.044</td><td>3.71E-01</td><td>14.618</td><td>1.00E+00</td><td>NS</td></tr><tr><td colspan="7"><strong>Commissural fibers</strong></td></tr><tr><td>CC</td><td>-0.462</td><td>0.055</td><td>5.08E-17</td><td>41.305</td><td>7.86E-12</td><td>***</td></tr><tr><td>BCC</td><td>-0.430</td><td>0.052</td><td>2.32E-16</td><td>35.479</td><td>5.41E-11</td><td>***</td></tr><tr><td>GCC</td><td>-0.373</td><td>0.066</td><td>1.78E-08</td><td>59.395</td><td>6.87E-06</td><td>***</td></tr><tr><td>SCC</td><td>-0.339</td><td>0.053</td><td>1.97E-10</td><td>37.906</td><td>5.66E-08</td><td>***</td></tr><tr><td>FX</td><td>-0.288</td><td>0.054</td><td>8.19E-08</td><td>39.029</td><td>7.84E-05</td><td>***</td></tr><tr><td>AverageFA</td><td>-0.260</td><td>0.076</td><td>5.69E-04</td><td>69.240</td><td>1.66E-01</td><td>NS</td></tr></tbody></table>

ns not significant

*p<sub>corr</sub> < 0.05; **p<sub>corr</sub> < 0.01; ***p<sub>corr</sub> < 0.001

mood switching [31]. Reduced FA within the CC was also reported in a meta-analysis of DTI studies in schizophrenia [25] and major depressive disorder [29], suggesting an overlapping involvement in both psychosis and affective disorders. Further studies are warranted to evaluate to what extent the CC is differentially affected in these disorders. Preliminary data suggest that disruption of interhemispheric connectivity is a disease marker rather than a vulnerability marker to BD [32]. Nonetheless, we identified extensive WM abnormalities suggesting that current pathophysiological models of BD are incomplete. Future models should not be limited to fronto-limbic networks, and should perhaps consider interhemispheric disconnectivity as a key feature of BD.

Importantly, the patient group was significantly older than the control group. Although we controlled for age in all analyses, it is possible that the linear models used are not fully accounting for the age-related variance [33]. However, the assessment of the effects of age revealed a significant interaction between age and diagnosis for only 4 ROIs out of the 43 analyzed. We found a significant increase in the effect of age in patients with BD for the left CGC only, while we found the reverse association for the bilateral SCR and the left PLIC, these effects were not anticipated and should be verified when replication samples become available.

We found that lithium intake was associated with higher FA in several tracts, as well as with global FA. Prior studies have suggested neuroprotective effects of lithium, on gray matter [23, 34–36] and white matter [37]. Higher FA associated with lithium use could reflect a direct influence of lithium on water diffusion or a beneficial effect on myelination [38], as suggested by the observation that lithium promotes myelin gene expression, morphological maturation, and remyelination in cultured oligodendrocytes via the Wnt/β-catenin and the Akt/CREB pathways [39]. In patients with BD, lithium may increase axial diffusivity in WM tracts also influenced by genetic variation in this pathway [22]. We also found lower FA in patients who received anticonvulsants in several tracts and average global FA. Further, patients who were on antipsychotic treatment showed lower FA within the genu of the CC. This is consistent with prior results suggesting a negative relationship between anticonvulsants, antipsychotics and cortical thickness or FA [23, 37]. However, it could be possible that the choice of the medication was driven by some patients’ particularities or unknown neurobiological characteristics, which are hard to assess with a cross-sectional design, leading to confounding by indication. Longitudinal clinical trials are needed to clarify this point.

We did not find significant differences between BD type I and type II. The power of prior meta-analyses of DTI studies has also been too low to perform this comparison. However, sensitivity analyses for these meta-analyses indicated that the sub-group of patients with BD I was driving the FA difference observed between patients with BD and HC [19, 29, 40]. Although we had enough power, the comparison of BD I vs. BD II did not replicate this result. Consistent with our results, however, ENIGMA analyses of T1-weighted anatomical MRI data of patients with BD did not yield any detectable differences between BD types [23, 24].

In sum, the multisite nature of the study is a strength that allowed us to detect small but significant differences. Our results seem to challenge the hypothesis of a precise localization for the WM alterations in BD. Indeed, we have highlighted extensive abnormalities, which do not seem to be specific to this psychiatric disorder. Lower FA across multiple bundles has already been consistently observed in studies of schizophrenia, with apparently higher effect sizes (e.g., [25]). Consequently, to build more precise neurobiological models of BD future studies should benefit from new advanced neuroimaging methods such as Neurite Orientation Dispersion and Density Imaging (NODDI) [41]. This recent processing model allows fine-grained measurement of the WM

In contrast, the role of the CC in pathophysiological models of BD is less straightforward. Disconnection in patients with BD with psychotic history has been suggested [12] but there is no clear evidence for the implication of the CC in emotion processing or
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Results of the meta-analysis. **a** Effect sizes for fractional anisotropy (FA) differences between patients with bipolar disorder (BD) and healthy controls projected on the 43 white matter (WM) tracts analyzed. **b** Cohen's *d* (effect size) sorted in increasing order of magnitude for significant differences between bipolar patients and healthy controls. Significant findings after Bonferroni correction are highlighted in blue. Error bars represent standard error</figcaption></figure>

microstructure, with physiological interpretation of the derived indices, and has already shown promising results in BD [42]. However, the large-scale application of such methods will only be possible with raw data sharing within international consortia. This will allow the application of advanced DTI models and whole-brain analyses, which are needed to better understand WM abnormalities observed in BD. Finally, longitudinal studies conducted in conjunction with advanced DTI protocols are essential to clarify the impact of pharmaceutical treatments on brain microstructure.

Some limitations are important to emphasize. We did not include other diffusion parameters in our analysis. Lower FA may represent abnormal fiber coherence but does not yield information on fiber density or myelination. The mean, radial and axial diffusivity measure would have added complementary information regarding the nature of WM alteration. However, we have focused on the most commonly used measure, which offers better comparability with prior studies. Also, most studies have highlighted a correlation between FA and these other measures, while their inclusion would have tripled the number of analyses. In addition, although we found "widespread" WM abnormalities in patients with BD, the robust ENIGMA DTI pipeline used to partition the ROIs involved only long and isolinear bundles. With this methodological approach (i.e., FSL TBSS), we cannot evaluate localized changes within the superficial WM, as have been previously observed in BD and schizophrenia [43]. Also, this methodological approach poorly reconstructs fiber crossings, which may have led to incomplete localization of group

differences. Further studies are warranted to identify more fine-
grained WM abnormalities in BD.

Importantly, retrospective multisite analyses have some limita-
tions. Differences in the acquisition parameters, magnet strength,
head coil and manufacturer provided software could have
impacted the results. However, we believe that our approach,
using a harmonized data processing pipeline, with a reliable
procedure, allows for the first time coordinated mega- and meta-
analyses with robust results.

Moreover, the effects of the covariates found here are only derived from post hoc analyses in cross-sectional studies with a somewhat limited representation of individuals with BD over age 50 (only 18% of the sample). Longitudinal studies would be more suitable to identify and predict the effect of age, illness duration/ severity and medication on WM microstructure in patients with BD. In addition, despite their importance, we were not able to test the relation between FA and other covariates, such as body mass index and frequent BD comorbidities (e.g., anxiety or substance use disorder). Too few sites had collected these measures to allow robust analyses. However, we believe that our sample is ecologically valid and captures the heterogeneity of BD.

With this unprecedented sample size, we found evidence for
widespread WM abnormalities in patients with BD and showed
differences in BD WM microstructure that were unobserved until
now. These results may inform future DTI studies with regard to
expected effect sizes, and the effects of several covariates and
clinical variables. We also highlighted that the CC and the
cingulum had the strongest decrease in FA in patients with BD. Despite growing evidence for altered structure of the CC in BD, its specific role in the pathophysiology of BD needs to be further integrated into neural models of BD.

## FUNDING AND DISCLOSURE

The researchers and studies included in this paper were funded by the: German Research Foundation (DFG, grant FOR2107 DA1151/5-1 and DA1151/5-2 to UD; SFB-TRR58, Projects C09 and Z02 to UD) and the Interdisciplinary Center for Clinical Research (IZKF) of the medical faculty of Münster (grant Dan3/012/17 to UD); German Research Foundation (SFB636/C6, WE3638/3-1); Oslo University Hospital, University of Oslo, Norwegian Research Council, South Eastern Norwegian Health Authorities; NIH Grant MH083968; VA Desert-Pacific Mental Illness Research Education and Clinical Center; Research Council of Norway (223273, 213837, 249711, 249795, 248238, 248778, and 262656), South East Norway Health Authority (2017-112) Kristian Gerhard Jebsen Stiftelsen (SKGJ-MED-008) and the European Community's Seventh Framework Program (FP7/2007-2013), grant agreement no. 602450 (IMAGEMEND); Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Brazil, 480370/2009-5); Australian National Medical and Health Research Council (Program Grant 1037196), and the Lansdowne Foundation; NIH grant U54 EB020403 from the BD2K Initiative, R01 MH116147, and P41 EB015922; CIBERSAM; NHG (SIG/12004) and SBIC (RP C-009) (KS); SAMRC (DJS); Canadian Institutes of Health Research (103703, 106469 and 64410); Nova Scotia Health Research Foundation, Dalhousie Clinical Research and Scholarship to T Hajek, NARSAD Young Investigator and Independent Investigator Awards (TH), JMAS SIM fellowship from the Royal College of Physicians of Edinburgh and an ESAT College Fellowship from the University of Edinburgh (HCW); Part of the Cardiff cohort was funded through a NARSAD Young Investigator Award (17319) (XC), and we are also grateful to the National Centre for Mental Health (NCMH) and the Bipolar Disorder Research Network for their support with recruitment; Spanish Ministry of Science, Innovation and Universities (PI15/00283) integrated into the Plan Nacional de I+D+I y cofinanciado por el ISCIII-Subdirección General de Evaluación y el Fondo Europeo de Desarrollo Regional (FEDER), CIBERSAM, and the Comissionat per a Universitats i Recerca del DIUE de la Generalitat de Catalunya to the Bipolar Disorders Group (2017 SGR 1365) and the project SLT006/17/00357, from PERIS 2016-2020 (Departament de Salut). CERCA Programme/Generalitat de Catalunya (EV); FAPESP-Brazil (#2009/14891-9, 2010/18672-7, 2012/23796-2 and 2013/03905-4), CNPq-Brazil (#478466/2009 and 480370/2009), the Wellcome Trust (UK) and the Brain & Behavior Research Foundation (NARSAD Independent Investigator Award (GFB); CRKC was supported by NIA T32AG058507; NIH/NIMH 5T32MH073526 and NIH grant U54EB020403 from the Big Data to Knowledge (BD2K) Program. CRKC has received partial research support from Biogen, Inc. (Boston, USA) for work unrelated to the topic of this manuscript; the Human Brain Project, funded from the European Union's Horizon 2020 Framework Program for Research and Innovation under the Specific Grant Agreements No. 785907 (SGA2) and No: 604102 (SGA1), and by the FRM DIC20161236445; IRMaGe MRI/Neurophysiology facility which was partly funded by the French program "Investissement d'Avenir" run by the "Agence Nationale pour la Recherche"; grant "Infrastructure d'avenir en Biologie Santé"—ANR-11-INBS-0006" and the Agence Nationale pour la Recherche (ANR-11-IDEX-0004 Labex BioPsy, ANR-10-COHO-10-01 psyCOH), Fondation pour la Recherche Médicale (Bioinformatique pour la biologie 2014) and the Fondation de l'Avenir (Recherche Médicale Appliquée 2014). PMT and NJ received research grant from Biogen, Inc., for research unrelated to this manuscript. Dr. Vieta has received grants and served as consultant, advisor or CME speaker for the following entities: AB-Biotics, Abbott, Allergan,

Angelini, Dainippon Sumitomo Pharma, Galenica, Janssen, Lundbeck, Novartis, Otsuka, Sage, Sanofi-Aventis, and Takeda. OAA has received speakers honorarium from Lundbeck, and is consultant for HealthLytix. DJS has received research grants and/or consultancy honoraria from Lundbeck and Sun. The remaining authors declare no competing interests.

## ADDITIONAL INFORMATION

**Supplementary Information** accompanies this paper at (https://doi.org/10.1038/s41386-019-0485-6).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## REFERENCES

1. Merikangas KR, Akiskal HS, Angst J, Greenberg PE, Hirschfeld R, Petukhova M, et al. Lifetime and 12-month prevalence of bipolar spectrum disorder in the National Comorbidity Survey replication. Arch Gen Psychiatry. 2007;64:543.

2. Goldberg JF, Ernst CL. Features associated with the delayed initiation of mood stabilizers at illness onset in bipolar disorder. J Clin Psychiatry. 2002;63:985-91.

3. Phillips ML, Swartz HA. A critical appraisal of neuroimaging studies of bipolar disorder: toward a new conceptualization of underlying neural circuitry and a road map for future research. Am J Psychiatry. 2014;171:829-43.

4. Strakowski SM, Adler CM, Almeida J, Altshuler LL, Blumberg HP, Chang KD, et al. The functional neuroanatomy of bipolar disorder: a consensus model. Bipolar Disord. 2012;14:313-25.

5. Favre P, Baciu M, Pichat C, Bougerol T, Polosan M. fMRI evidence for abnormal resting-state functional connectivity in euthymic bipolar patients. J Affect Disord. 2014;165:182-9.

6. Favre P, Polosan M, Pichat C, Bougerol T, Baciu M. Cerebral correlates of abnormal emotion conflict processing in euthymic bipolar patients: a functional MRI study. PLoS ONE. 2015;10:e0134961.

7. Houenou J, Wessa M, Douaud G, Leboyer M, Chanraud S, Perrin M, et al. Increased white matter connectivity in euthymic bipolar patients: diffusion tensor tractography between the subgenual cingulate and the amygdalo-hippocampal complex. Mol Psychiatry. 2007;12:1001-10.

8. Versace A, Almeida JRC, Hassel S, Walsh ND, Novelli M, Klein CR, et al. Elevated left and reduced right orbitomedial prefrontal fractional anisotropy in adults with bipolar disorder revealed by tract-based spatial statistics. Arch Gen Psychiatry. 2008;65:1041-52.

9. Wang F, Jackowski M, Kalmar JH, Chepenik LG, Tie K, Qiu M, et al. Abnormal anterior cingulum integrity in bipolar disorder determined through diffusion tensor imaging. Br J Psychiatry. 2008;193:126-9.

10. Wang F, Kalmar JH, Edmiston E, Chepenik LG, Bhagwagar Z, Spencer L, et al. Abnormal corpus callosum integrity in bipolar disorder: a diffusion tensor imaging study. Biol Psychiatry. 2008;64:730-33.

11. Leow A, Ajilore O, Zhan L, Arienzo D, GadElkarim J, Zhang A, et al. Impaired interhemispheric integration in bipolar disorder revealed with brain network analyses. Biol Psychiatry. 2013;73:183-93.

12. Sarrazin S, Poupon C, Linke J, Wessa M, Phillips M, Delavest M, et al. A multicenter tractography study of deep white matter tracts in bipolar I disorder: psychotic features and interhemispheric disconnectivity. JAMA Psychiatry. 2014;71:388.

13. Yurgelun-Todd DA, Silveri MM, Gruber SA, Rohan ML, Pimentel PJ. White matter abnormalities observed in bipolar disorder: a diffusion tensor imaging study. Bipolar Disord. 2007;9:504-12.

14. Roberts G, Wen W, Frankland A, Perich T, Holmes-Preston E, Levy F, et al. Interhemispheric white matter integrity in young people with bipolar disorder and at high genetic risk. Psychol Med. 2016;46:2385-96.

15. Repple J, Meinert S, Grotegerd D, Kugel H, Redlich R, Dohm K, et al. A voxel-based diffusion tensor imaging study in unipolar and bipolar depression. Bipolar Disord. 2017;19:23-31.

16. Cui L, Chen Z, Deng W, Huang X, Li M, Ma X, et al. Assessment of white matter abnormalities in paranoid schizophrenia and bipolar mania patients. Psychiatry Res Neuroimaging. 2011;194:347-53.

17. Basser PJ, Mattiello J, LeBihan D. MR diffusion tensor spectroscopy and imaging. Biophys J. 1994;66:259-67.

18. Vederine F-E, Wessa M, Leboyer M, Houenou J. A meta-analysis of whole-brain diffusion tensor imaging studies in bipolar disorder. Prog Neuropsychopharmacol Biol Psychiatry. 2011;35:1820-6.
19. Nortje G, Stein DJ, Radua J, Mataix-Cols D, Horn N. Systematic review and voxel-based meta-analysis of diffusion tensor imaging studies in bipolar disorder. *J Affect Disord.* 2013;150:192-200.

20. Zanetti MV, Jackowski MP, Versace A, Almeida JRC, Hassel S, Duran FbLS, et al. State-dependent microstructural white matter changes in bipolar I depression. *Eur Arch Psychiatry Clin Neurosci.* 2009;259:316-28.

21. Ha TH, Her JY, Kim JH, Chang JS, Cho HS, Ha K. Similarities and differences of white matter connectivity and water diffusivity in bipolar I and II disorder. *Neurosci Lett.* 2011;505:150-4.

22. Benedetti F, Bollettini I, Barberi I, Radaelli D, Poletti S, Locatelli C, et al. Lithium and GSK3-β promoter gene variants influence white matter microstructure in bipolar disorder. *Neuropsychopharmacology.* 2013;38:313.

23. Hibar D, Westlye L, Doan N, Jahanshad N, Cheung J, Ching C, et al. Cortical abnormalities in bipolar disorder: an MRI analysis of 6503 individuals from the ENIGMA Bipolar Disorder Working Group. *Mol Psychiatry.* 2018;23:932.

24. Hibar D, Westlye LT, van Erp TG, Rasmussen J, Leonardo CD, Faskowitz J, et al. Subcortical volumetric abnormalities in bipolar disorder. *Mol Psychiatry.* 2016;21:1710.

25. Kelly S, Jahanshad N, Zalesky A, Kochunov P, Agartz I, Alloza C, et al. Widespread white matter microstructural differences in schizophrenia across 4322 individuals: results from the ENIGMA Schizophrenia DTI Working Group. *Mol Psychiatry.* 2017;23:1261-9.

26. Mori S, Oishi K, Jiang H, Jiang L, Li X, Akhter K, et al. Stereotaxic white matter atlas based on diffusion tensor imaging in an ICBM template. *Neuroimage.* 2008;40:570-82.

27. Bates D, Maechler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. *Journal of Statistical Software*, 2015;67:1-48.

28. Kochunov P, Jahanshad N, Sprooten E, Nichols TE, Mandl RC, Almasy L, et al. Multi-site study of additive genetic effects on fractional anisotropy of cerebral white matter: comparing meta and mega-analytical approaches for data pooling. *Neuroimage.* 2014;95:136-50.

29. Wise T, Radua J, Nortje G, Cleare AJ, Young AH, Arnone D. Voxel-based meta-analytical evidence of structural disconnectivity in major depression and bipolar disorder. *Biol Psychiatry.* 2016;79:293-302.

30. Mahon K, Burdick KE, Szeszko PR. A role for white matter abnormalities in the pathophysiology of bipolar disorder. *Neurosci Biobehav Rev.* 2010;34:533-54.

31. Anderson LB, Paul LK, Brown WS. Emotional intelligence in agenesis of the corpus callosum. *Arch Clin Neuropsychol.* 2017;32:267-79.

32. Linke J, King AV, Poupon C, Hennerici MG, Gass A, Wessa M. Impaired anatomical connectivity and related executive functions: differentiating vulnerability and disease marker in bipolar disorder. *Biol Psychiatry.* 2013;74:908-16.

33. Miller GA, Chapman JP. Misunderstanding analysis of covariance. *J Abnorm Psychol.* 2001;110:40.

34. Hajek T, Bauer M, Pfennig A, Cullis J, Ploch J, O'Donovan C, et al. Large positive effect of lithium on prefrontal cortex N-acetylaspartate in patients with bipolar disorder: 2-centre study. *J Psychiatry Neurosci.* 2012;37:185.

35. Lyoo IK, Dager SR, Kim JE, Yoon SJ, Friedman SD, Dunner DL, et al. Lithium-induced gray matter volume increase as a neural correlate of treatment response in bipolar disorder: a longitudinal brain imaging study. *Neuropsychopharmacology.* 2010;35:1743-50.

36. Moore GJ, Cortese BM, Glitz DA, Zajac-Benitez C, Quiroz JA, Uhde TW, et al. A longitudinal study of the effects of lithium treatment on prefrontal and subgenual prefrontal gray matter volume in treatment-responsive bipolar disorder patients. *J Clin Psychiatry.* 2009;70:699-705.

37. Abramovic L, Boks MP, Vreeker A, Verkooijen S, van Bergen AH, Ophoff RA, et al. White matter disruptions in patients with bipolar disorder. *Eur Neuropsychopharmacol.* 2018;28:743-51.

38. Marlinge E, Bellivier F, Houenou J. White matter alterations in bipolar disorder: potential for drug discovery and development. *Bipolar Disord* 2014;16:97-112.

39. Meffre D, Massaad C, Grenier J. Lithium chloride stimulates PLP and MBP expression in oligodendrocytes via Wnt/β-catenin and Akt/CREB pathways. *Neuroscience.* 2015;284:962-71.

40. Yang C, Li L, Hu X, Luo Q, Kuang W, Lui S, et al. Psychoradiologic abnormalities of white matter in patients with bipolar disorder: diffusion tensor imaging studies using tract-based spatial statistics. *J Psychiatry Neurosci.* 2019;44:32.

41. Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC. NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. *Neuroimage.* 2012;61:1000-16.

42. Sarrazin S, Poupon C, Teillac A, Mangin J-F, Polosan M, Favre P, et al. Higher in vivo cortical intracellular volume fraction associated with lithium therapy in bipolar disorder: a multicenter NODDI study. *Psychother Psychosom.* 2019;88:171-6.

43. Ji E, Guevara P, Guevara M, Grigis A, Labra N, Sarrazin S, et al. Increased and decreased superficial white matter structural connectivity in schizophrenia and bipolar disorder. *Schizophr Bull.* 2019:sbz015. https://doi.org/10.1093/schbul/sbz015.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.

org/licenses/by/4.0/.

© The Author(s) 2019