Introduction
Obesity is the fifth leading cause of death globally and is one of the leading causes of disability (Di Angelantonio et al., Reference Di Angelantonio, Bhupathiraju, Wormser, Gao, Kaptoge, de Gonzalez and Hu2016; GBD 2015 Obesity Collaborators et al., Reference Afshin, Forouzanfar, Reitsma, Sur, Estep and Murray2017; Nyberg et al., Reference Nyberg, Batty, Pentti, Virtanen, Alfredsson, Fransson and Kivimäki2018). It is an even greater problem among people with severe mental illness (SMI) including bipolar disorders (BD). Based on meta-analysis of 120 studies, the pooled point prevalence of overweight and obesity in SMI was 60% and people with SMI had on average three times greater odds of obesity than the general population (Afzal et al., Reference Afzal, Siddiqi, Ahmad, Afsheen, Aslam, Ali and Zavala2021). Aside from the impact of obesity on general health, obesity is commonly associated with structural brain alterations (Dekkers, Jansen, & Lamb, Reference Dekkers, Jansen and Lamb2019; Fernández-Andújar, Morales-García, & García-Casares, Reference Fernández-Andújar, Morales-García and García-Casares2021; García-García et al., Reference García-García, Michaud, Dadar, Zeighami, Neseliler, Collins and Dagher2019; Janowitz et al., Reference Janowitz, Wittfeld, Terock, Freyberger, Hegenscheid, Volzke and Grabe2015; Willette & Kapogiannis, Reference Willette and Kapogiannis2015), and with an increased risk of cognitive impairment and dementia (Beydoun, Beydoun, & Wang, Reference Beydoun, Beydoun and Wang2008; Pedditzi, Peters, & Beckett, Reference Pedditzi, Peters and Beckett2016; Singh-Manoux et al., Reference Singh-Manoux, Dugravot, Shipley, Brunner, Elbaz, Sabia and Kivimaki2018; Tang et al., Reference Tang, Zhao, Lu, Zhang, Zhang, Xin and Stehouwer2021). These issues may be particularly relevant in individuals who already have an increased risk of brain alterations (Hibar et al., Reference Hibar, Westlye, Doan, Jahanshad, Cheung, Ching and Andreassen2018), cognitive impairment (Bora, Yucel, & Pantelis, Reference Bora, Yucel and Pantelis2009), and obesity (Vancampfort et al., Reference Vancampfort, Stubbs, Mitchell, De Hert, Wampers, Ward and Correll2015), such as people with BD. Individuals with BD and comorbid obesity face very specific challenges which require dedicated management and research efforts. However, very few studies up to date have investigated how the presence of obesity interacts with brain and cognitive changes in major psychiatric disorders.
Studying people with SMI and obesity could help identify preventable/treatable risk factors for neurostructural alterations, which may be associated with currently intractable psychiatric outcomes, including cognitive impairment. Indeed, we and others have shown that BD complicated by obesity-related metabolic alterations, specifically diabetes, is associated with lower psychosocial functioning (Hajek et al., Reference Hajek, Slaney, Garnham, Ruzickova, Passmore and Alda2005), higher rates of rapid cycling (Hajek et al., Reference Hajek, Hahn, Slaney, Garnham, Green, Ruzickova and Alda2008), poor treatment response (Calkin et al., Reference Calkin, Ruzickova, Uher, Hajek, Slaney, Garnham and Alda2015), and worse psychiatric outcomes (Calkin et al., Reference Calkin, van de Velde, Ruzickova, Slaney, Garnham, Hajek and Alda2009). It could also explain why people with the same psychiatric diagnosis differ so markedly in their neurobiological/clinical outcomes, thus moving toward individualized medicine. Furthermore, weight gain is a common side effect of many psychiatric medications, which are also frequently associated with changes in brain structure. Therefore, we need to better understand the role obesity plays in the links between psychiatric disorders or medications and brain structure.
Previous studies in relatively small (76–112 participants) and highly selected groups [i.e. people with first episode of mania (Bond et al., Reference Bond, Ha, Lang, Su, Torres, Honer and Yatham2014, Reference Bond, Lang, Noronha, Kunz, Torres, Su and Yatham2011, Reference Bond, Su, Honer, Dhanoa, Batres, Lee and Yatham2019), adolescent BD participants (Islam, Metcalfe, MacIntosh, Korczak, & Goldstein, Reference Islam, Metcalfe, MacIntosh, Korczak and Goldstein2018), or offspring of people with BD (Mansur et al., Reference Mansur, McIntyre, Cao, Lee, Japiassú, Chen and Lin2018)] have suggested that obesity may be associated with brain alterations in BD, possibly with a stronger effect size or with some regional specificity compared to controls (Bond et al., Reference Bond, Ha, Lang, Su, Torres, Honer and Yatham2014, Reference Bond, Lang, Noronha, Kunz, Torres, Su and Yatham2011, Reference Bond, Su, Honer, Dhanoa, Batres, Lee and Yatham2019; Islam et al., Reference Islam, Metcalfe, MacIntosh, Korczak and Goldstein2018; Mansur et al., Reference Mansur, McIntyre, Cao, Lee, Japiassú, Chen and Lin2018). In a study including 2735 individuals (McWhinney et al., Reference McWhinney, Abé, Alda, Benedetti, Bøen and Del Mar Bonnin2021a), we demonstrated that some of the most replicated subcortical brain alterations in BD, including larger ventricles, were to a large extent (up to 47%) mediated by obesity and that both BD and obesity were associated with similar subcortical alterations.
We need larger studies in more generalizable samples to better understand how the brain correlates of obesity map onto the cortical alterations in BD. To this end, we investigated the association between BD, medications, obesity, and neurostructural measures in a large, highly generalizable, multicenter sample from the ENIGMA-BD Working Group.
Methods
Participating sites
The ENIGMA-BD Working Group aims to improve replication and generalizability of neuroimaging studies of BD by combining existing, independently collected neuroimaging samples of BD from around the world (Ching et al., Reference Ching, Hibar, Gurholt, Nunes, Thomopoulos and Abé2022; Hibar et al., Reference Hibar, Westlye, Doan, Jahanshad, Cheung, Ching and Andreassen2018, Reference Hibar, Westlye, van Erp, Rasmussen, Leonardo, Faskowitz and Andreassen2016; McWhinney et al., Reference McWhinney, Abé, Alda, Benedetti, Bøen and Del Mar Bonnin2021a, Reference McWhinney, Abé, Alda, Benedetti, Bøen and Del Mar Bonnin2022a; Nunes et al., Reference Nunes, Schnack, Ching, Agartz, Akudjedu and Alda2020). Seventeen independently collected ENIGMA-BD samples from 13 countries on six continents contributed individual-level structural MRI data, medication information, specifically medications used at the time of scanning for the following medication categories (lithium, first-, second-generation antipsychotics, anticonvulsants, antidepressants), and body mass index (BMI) values from a total of 1231 individuals with BD and 1601 healthy controls. Table 1 shows all participant characteristics.
Demographic, diagnostic, and treatment characteristics of sample
a There were no missing age or BMI values. We used the Welch two-sample t test (unequal variance assumed), which relies on a Welch–Satterthwaite degrees of freedom adjustment, resulting in varying degrees of freedom.
Online Supplementary Tables S1 and S2 list the demographic/clinical details for each cohort. Online Supplementary Table S3 provides the diagnostic instruments used to obtain diagnosis and clinical information. Online Supplementary Table S4 lists exclusion criteria for study enrolment. Briefly, all studies used standard diagnostic instruments, including Structured Clinical Interview for DSM Disorders (SCID; N = 12), Mini International Neuropsychiatric Interview (MINI; N = 2) and Diagnostic Interview for Genetic Studies (DIGS; N = 1). Most studies (N = 10) included both bipolar I (BDI) and bipolar II (BDII) disorders, six studies included only BDI and one study included only BDII participants. Substance abuse was an exclusion criterion in nine studies. Most studies did not exclude comorbidities, other than substance abuse. Consequently, the sample is a broad, ecologically valid, and generalizable representation of BD. All participating sites received approval from local ethics committees, and all participants provided written informed consent. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.
Data acquisition and parcellation
High-resolution T1-weighted brain anatomical MRI scans were acquired at each site; see online Supplementary Table S5 for details of scan acquisition. All groups used the same ENIGMA-standardized analytical protocol, including visual and statistical quality assessment, as documented at: http://enigma.ini.usc.edu/protocols/imaging-protocols/. These protocols are standardized across the consortium, are open-source, and freely available online, in order to foster open science/replication/reproducibility. They were applied in the previous publications by our group (Hibar et al., Reference Hibar, Westlye, Doan, Jahanshad, Cheung, Ching and Andreassen2018; Nunes et al., Reference Nunes, Schnack, Ching, Agartz, Akudjedu and Alda2020) and more broadly in large-scale ENIGMA studies of major depression, schizophrenia, attention deficit hyperactivity disorder (ADHD), obsessive compulsive disorder (OCD), post traumatic stress disorder (PTSD), epilepsy, and autism (Thompson et al., Reference Thompson, Jahanshad, Ching, Salminen, Thomopoulos and Bright2020).
Briefly, using the freely available and extensively validated FreeSurfer software, we performed parcellations of 34 cortical regions, per hemisphere (left and right), based on the Desikan–Killiany atlas. All segmented regions were used as target regions of interest (ROIs) for analysis. We also computed total intracranial volume (ICV) to standardize surface area estimates. Visual quality controls were performed on an ROI level aided by the ENIGMA-standardized visual inspection guide including pass/fail parcellation examples. In addition, we generated diagnostic histogram plots for each site and outliers which deviated from the site mean for each structure at >3 standard deviations were flagged for further review. All ROIs failing quality inspection were withheld from subsequent analyses, see online Supplementary Table S6. Previous analyses from the ENIGMA-BD Working Group showed that scanner field strength, voxel volume, and the version of FreeSurfer used for parcellation did not significantly influence the effect size estimates.
Statistical modeling
In this mega-analysis, we used linear mixed modeling (package nlme version 3.1-152 in R version 4.1.1) with individual subject cortical thickness or cortical surface area as dependent variables and with both BMI and group (participants with BD or healthy controls) as predictors. In each case, age, sex, and hemisphere (left or right) were also included as fixed predictors. Total ICV was included as a covariate in models of cortical surface area. Models also included a random effect of hemisphere within participants and a random effect of data collection site.
We created one model per region, with each model including both hemispheres and all of the covariates described above. We used BMI as a continuous variable, which captures more variability between participants, increases sensitivity, and was the preferred approach in most previous studies (Dekkers et al., Reference Dekkers, Jansen and Lamb2019). BMI was normally distributed (online Supplementary Fig. S1). We checked the normality of model residuals using QQ plots and tested for multicollinearity using the variance inflation factor (VIF, shown in online Supplementary Table S7) of all predictor variables included in modeling. Variance in regional volumes was comparable between groups.
In post hoc analyses among individuals with BD, we separately explored the statistical effects of commonly prescribed medications. As the rates of monotherapy were low in this sample (see Table 1), we studied the association between number of jointly used medication classes (zero through three, including anticonvulsants, antipsychotics, and antidepressants) and BMI or cortical thickness or surface area. In the same model, we separately estimated the effects of current Li treatment. We used the same covariates and random-effect structure as described above. Interactions between BMI and either the number of medication classes or Li prescription were included where significant. The partial effect of the number of medication classes while adjusting for BMI was also compared with its effect without adjusting for BMI, but with all other covariates and random effects remaining. The a priori decision to analyze the effects of Li separately was motivated by the fact that statistical effects of Li on brain measures, which tend to be positive, may cancel the statistical effects of other medications such as antipsychotics, anticonvulsants, which tend to be negative (Hajek et al., Reference Hajek, Cullis, Novak, Kopecek, Hoschl, Blagdon and Alda2012a; Hajek, Kopecek, Hoschl, & Alda, Reference Hajek, Kopecek, Hoschl and Alda2012b).
We adjusted all p values for multiple comparisons using false discovery rate (FDR), with adjusted p values reported, at α = 0.05. We calculated effect sizes for between-group differences (partial d), and associations between BMI and ROI volumes (partial r), together with their 95% confidence intervals (CIs), using model coefficients and their standard error (s.e.) (Nakagawa & Cuthill, Reference Nakagawa and Cuthill2007).
Mediation analysis
We tested whether the variance in regional thickness that was associated with the number of jointly used medication classes (zero through three, including anticonvulsants, antipsychotics, and antidepressants; direct path) remained significant after also accounting for variance associated with BMI (indirect path) in individuals with BD. The number of medication classes was modeled as the associated variable, BMI as the mediating variable, and cortical thickness was the dependent variable. We modeled the direct effect of the number of jointly used medication classes on thickness, in comparison with the indirect effect of this association through BMI as a mediator, corrected for age, sex, prescription of Li, and random effects. To test this, we built 5000 bootstrapped models using random selection with replacement. This method non-parametrically identified the 95% CI for effect sizes. The bootstrap CI, which did not include zero, indicated a significant indirect effect. Simulation research indicates that the bootstrap method is more robust to non-normality and has better type I error control than the Sobel test (Hayes, Reference Hayes2009). Nevertheless, for methodological consistency, we also applied the Sobel test to investigate whether accounting for BMI significantly mitigated group-related differences in thickness. All of these analyses were performed in R version (4.1.1).
These analyses were applied only to regions which met the criteria for mediation: (1) the number of medication classes was a significant predictor of the ROI thickness, and (2) the number of medication classes was a significant predictor of the mediator (BMI), and when modelled jointly, (3) BMI was a significant predictor of the thickness, and (4) the strength of the coefficient of the previously significant independent variable (number of medication classes) was reduced. These criteria applied to the fusiform gyrus.
Results
Regional morphometric differences by diagnosis and BMI
When modeled jointly, numerous regions showed significant partial effects of either BMI, diagnosis, or both (Fig. 1, Table 2). Participants with BD showed significantly thinner cortex relative to controls in all regions except for the entorhinal cortex and temporal pole. Higher BMI was associated with thinner cortex in nine of the same regions as BD, and it was uniquely associated with thinner entorhinal cortex. Surface area did not significantly differ between groups in any region, while higher BMI was associated with larger surface area in the isthmus of the cingulate gyrus (online Supplementary Table S8).
Standardized coefficients for group differences (blue) and BMI effects (red) in predicting the cortical thickness of each region. Significant effects are shown using a filled marker (FDR-adjusted p < 0.05).
Results of multiple regression analysis in cortical thickness, including effect sizes for between-group differences (Cohen's d), 95% confidence interval, BMI effect sizes (part r), and their FDR-adjusted p values
Significance is shown using asterisks (*p < 0.05)
BMI and group significantly interacted in lateral occipital cortical thickness (online Supplementary Fig. S2), with control participants showing a significant negative association between BMI and cortical thickness [t(2391) = −3.03, p = 0.002], while no significant association was seen in those with BD [t(2391) = 0.31, p = 0.757]. There was no interaction between BMI and sex for any of the regions or any of the measures. The full list of interactions is shown in online Supplementary Table S9.
Medications, clinical variables, BMI, and brain structure
In individuals with BD, higher BMI was associated with greater number of jointly used medication classes (i.e. anticonvulsant, antipsychotic, and/or antidepressant medications) per participant [t(1100) = 4.89, p < 0.001], but not with Li treatment [t(736) = −0.42, p = 0.676].
The number of medication classes was significantly associated with smaller cortical thickness in 22 of 34 regions (64.7%), and in nearly all instances, these associations remained significant when controlling for the effects of BMI (online Supplementary Table S10). Exceptions included the pars opercularis, superior temporal gyrus, and supramarginal gyrus. There was a significant interaction between BMI and the number of medication classes in the isthmus cingulate [t(1858) = −4.34, p = 0.001], with progressively steeper associations between BMI and thickness in those with more medications (see online Supplementary Fig. S3).
There was an interaction between current Li-use and BMI in 13 of the ROIs (38.2%), including caudal and rostral anterior cingulate, medial orbitofrontal gyrus, postcentral gyrus, pars opercularis, pars triangularis, rostral and caudal middle frontal, superior frontal, superior temporal, supramarginal, frontal pole, and insula, such that people who were prescribed Li at the time of scanning showed stronger negative association between BMI and cortical thickness than individuals with BD who were not treated with Li (online Supplementary Table S10 and Fig. S4). Amongst the remaining regions, which did not show interaction between BMI and Li, Li was positively associated with cortical thickness even when controlling for the negative effect of BMI in nine regions, including cuneus, precuneus, inferior parietal, lateral occipital, lingual, paracentral, precentral, pericalcarine, and superior parietal gyri.
BMI, Li treatment, and the number of medication classes showed negligible multicollinearity (VIF < 1.004). BMI was not significantly associated with illness duration, history of psychotic symptoms, diagnostic subtype, or mood state (online Supplementary Table S11).
Mediating effect of medications
Only the fusiform gyrus met the criteria for investigating whether BMI mediates the relationship between the number of medication classes and cortical thickness (online Supplementary Table S10). Specifically, there was a significant indirect effect of the number of medication classes on lower fusiform gyrus thickness through BMI (Est = −0.015, 95% CI −0.022 to −0.008), with 34.6% mediation (Z = 3.46, p < 0.001, see Fig. 2).
The effect of medication classes and BMI on cortical thickness. Path (c) represents the direct effect, while (a) through (b) represent the indirect path through BMI, and (c′) represents the adjusted direct effect after accounting for BMI. We show standardized coefficients along with their 95% CI derived from bootstrapping. Significant effects (95% CI that excludes zero) are marked by asterisks. In all models, we controlled for the covariates age, sex, Li treatment, and data collection site, while paths b, c, and c’ were additionally adjusted for a random effect of hemisphere.
Discussion
In this study of 2832 individuals, we found substantial overlap between regions associated with BMI and BD. Specifically, with exception of a single ROI (entorhinal cortex), all of the regions which were negatively associated with BMI were also negatively associated with the diagnosis of BD. In contrast, only a single ROI (isthmus of the cingulate gyrus) showed association between BMI and surface area, which was positive. Importantly, about a third of the negative association between number of psychiatric medications and cortical thickness in the fusiform gyrus was partially mediated by the association between number of medications and higher BMI. Furthermore, the association between BMI and cortical thickness in isthmus cingulate became steeper with higher number of medication classes used jointly. In contrast, significant negative associations between BMI and cortical thickness in isthmus cingulate and rostral anterior cingulate became non-significant when modelled jointly with a significant positive effect of Li on cortical thickness. There was an interaction between Li and BMI, such that the brain correlates of BMI were more pronounced in individuals with v. without current Li treatment in a number of frontal regions, as well as anterior cingulate. The statistical effect of BMI on brain structure was linear in all regions, thus it would be most pronounced in people with obesity, but also manifest in overweight individuals.
The cortical correlates of obesity in this study closely replicated findings from previous large-scale studies, which also reported negative associations between obesity and cortical thickness or volume, especially in caudal and rostral anterior cingulate, entorhinal cortex, and several frontal lobe regions (Janowitz et al., Reference Janowitz, Wittfeld, Terock, Freyberger, Hegenscheid, Volzke and Grabe2015; McWhinney et al., Reference McWhinney, Brosch, Calhoun, Crespo-Facorro, Crossley, Dannlowski and Hajek2022b; Opel et al., Reference Opel, Thalamuthu, Milaneschi, Grotegerd, Flint, Leenings and Dannlowski2021), but with much fewer correlates in surface area (McWhinney et al., Reference McWhinney, Brosch, Calhoun, Crespo-Facorro, Crossley, Dannlowski and Hajek2022b; Opel et al., Reference Opel, Thalamuthu, Milaneschi, Grotegerd, Flint, Leenings and Dannlowski2021). These are some of the same regions which are consistently associated with BD, but also with other major psychiatric disorders (Hibar et al., Reference Hibar, Westlye, Doan, Jahanshad, Cheung, Ching and Andreassen2018; van Erp et al., Reference van Erp, Walton, Hibar, Schmaal, Jiang, Glahn and Turner2018). Interestingly, in all of these regions both BMI and BD showed partial association with cortical thickness, when controlling for the other factor. In other words, brain changes in cingulate and frontal regions will be greater in people with both obesity and BD than in those with either condition alone.
It is unclear whether these regions are particularly susceptible to obesity, or whether their changes predispose individuals to obesity. The mechanisms through which brain structure could predispose to obesity involve alterations in the reward system (Opel et al., Reference Opel, Redlich, Grotegerd, Dohm, Haupenthal, Heindel and Dannlowski2015), cue triggered learning and Pavlovian conditioning to hedonic food (Meyer, Risbrough, Liang, & Boutelle, Reference Meyer, Risbrough, Liang and Boutelle2015), and in appetitive behavior (Löscher, Brandt, & Ebert, Reference Löscher, Brandt and Ebert2003; Malkova, Mishkin, Suomi, & Bachevalier, Reference Malkova, Mishkin, Suomi and Bachevalier2010). On the other hand, obesity could affect brain structure through a range of mechanisms, including among others systemic inflammation, oxidative stress, insulin resistance/diabetes, hypertension or dyslipidemia (Cox et al., Reference Cox, Lyall, Ritchie, Bastin, Harris, Buchanan and Deary2019; Goldstein et al., Reference Goldstein, Baune, Bond, Chen, Eyler, Fagiolini and Fiedorowicz2020; Hajek et al., Reference Hajek, Calkin, Blagdon, Slaney, Uher and Alda2014; Parimisetty et al., Reference Parimisetty, Dorsemans, Awada, Ravanan, Diotel and Lefebvre d'Hellencourt2016; Van Gaal, Mertens, & De Block, Reference Van Gaal, Mertens and De Block2006; Willette & Kapogiannis, Reference Willette and Kapogiannis2015; Wisse, Reference Wisse2004). Speculating about how the regions described above could predispose to obesity would be post-hoc and inconclusive. However, previous studies investigating associations between genetic risk for obesity and brain structure reported involvement of a much smaller set of regions, i.e. surface area of lateral occipital lobe (Opel et al., Reference Opel, Thalamuthu, Milaneschi, Grotegerd, Flint, Leenings and Dannlowski2021) or medial prefrontal volume (Opel et al., Reference Opel, Redlich, Kaehler, Grotegerd, Dohm, Heindel and Dannlowski2017) than those reported here and in other previous studies of association between obesity and brain structure (Janowitz et al., Reference Janowitz, Wittfeld, Terock, Freyberger, Hegenscheid, Volzke and Grabe2015; McWhinney et al., Reference McWhinney, Brosch, Calhoun, Crespo-Facorro, Crossley, Dannlowski and Hajek2022b; Opel et al., Reference Opel, Thalamuthu, Milaneschi, Grotegerd, Flint, Leenings and Dannlowski2021). Therefore, brain correlates of obesity most likely include both causes and consequences of obesity. Considering the greater extent of associations between obesity and brain structure relative to the cih the much smaller extent of associations between genetic risk for obesity and brain structure, most brain changes likely represent consequences of obesity.
An especially interesting and relevant question is the role of obesity in brain effects of psychiatric medications. In keeping with other studies, we found that antipsychotics and anticonvulsants were negatively associated with brain structure (Andreasen, Liu, Ziebell, Vora, & Ho, Reference Andreasen, Liu, Ziebell, Vora and Ho2013; Fleisher et al., Reference Fleisher, Truran, Mai, Langbaum, Aisen, Cummings and Tariot2011; Fusar-Poli et al., Reference Fusar-Poli, Smieskova, Kempton, Ho, Andreasen and Borgwardt2013; Hibar et al., Reference Hibar, Westlye, van Erp, Rasmussen, Leonardo, Faskowitz and Andreassen2016, Reference Hibar, Westlye, Doan, Jahanshad, Cheung, Ching and Andreassen2018; Tariot et al., Reference Tariot, Schneider, Cummings, Thomas, Raman, Jakimovich and Aisen2011; Van Gestel et al., Reference Van Gestel, Franke, Petite, Slaney, Garnham, Helmick and Hajek2019) and positively with BMI (Mitchell et al., Reference Mitchell, Vancampfort, Sweers, van Winkel, Yu and De Hert2013; Tek et al., Reference Tek, Kucukgoncu, Guloksuz, Woods, Srihari and Annamalai2016). Considering the negative association between obesity and brain structure, perhaps the negative associations between medications and brain structure could be mediated by obesity (Joober, Schmitz, Malla, Sengupta, & Karma, Reference Joober, Schmitz, Malla, Sengupta and Karma2006; McWhinney et al., Reference McWhinney, Abé, Alda, Benedetti, Bøen and Del Mar Bonnin2021a, Reference McWhinney, Kolenic, Franke, Fialova, Knytl, Matejka and Hajek2021b). Across most regions, the association between number of medications and cortical thickness remained significant when we controlled for BMI, which is in keeping with another study (Jorgensen et al., Reference Jorgensen, Nesvag, Nerland, Morch-Johnsen, Westlye, Lange and Agartz2017). In a single region, the fusiform gyrus, about a third of the association between medications and cortical thickness was related to the indirect association between medications and BMI and between BMI and brain structure. In a single region, isthmus cingulate, the association between BMI and cortical thickness became steeper with growing numbers of medications used jointly. So, while the obesitogenic effects of medications may mediate or moderate the negative association between medications and brain structure in some regions, for the most part, the negative statistical effect of medications on brain structure was independent from their obesitogenic effects.
A separate question is the interplay between the putative neuroprotective effects of Li and its impact on weight. Interestingly, in this study, Li remained positively associated with cortical thickness across numerous regions, even when we controlled for negative associations between BMI and regional thickness. Conversely, a significant negative association between BMI and cortical thickness became non-significant in the isthmus cingulate and rostral anterior cingulate when modelled jointly with the significant positive effect of Li on cortical thickness. In other words, BMI did not cancel the positive association between Li and cortical thickness, while Li did cancel the negative association between BMI and cortical thickness for some regions.
Interestingly, in some regions, people who were using Li at the time of scanning showed stronger negative association between BMI and cortical thickness than BD individuals not treated with Li. These findings were quite robust and replicated across approximately one-third of regions. It is possible that regions which are negatively associated with BD will not show negative effects of additional variables, such as BMI, unless the impact of the illness is mitigated by, for example, Li. Indeed, we found this interaction mostly in regions which were more strongly associated with diagnosis than with BMI, including medial orbitofrontal gyrus as well as pars opercularis and triangularis of the inferior frontal gyrus. In addition, this interaction among individuals with BD could have decreased the apparent effect of BMI in these regions. Indeed, in the whole group, BMI was not associated with thickness of these regions. This is also in keeping with greater negative association between BMI and cortical thickness in controls than in people with BD.
Since we do not understand the origins of brain alterations in BD, it is highly relevant to study variables which could be associated with the structure of the brain, such as BMI. On the theoretical level, such studies could help explain the differences in brain measures among people with the same diagnosis, as they may differ in BMI. They could also explain similarities in brain measures across people with BD or schizophrenia, as they share a high risk of obesity. On the clinical level, comorbid obesity in people with major psychiatric disorders is associated with poor functioning, greater risk of chronicity, disability and suicide, poor treatment response, and functional deterioration (Bora, Akdede, & Alptekin, Reference Bora, Akdede and Alptekin2017; Calkin et al., Reference Calkin, Ruzickova, Uher, Hajek, Slaney, Garnham and Alda2015, Reference Calkin, van de Velde, Ruzickova, Slaney, Garnham, Hajek and Alda2009; Hajek et al., Reference Hajek, Hahn, Slaney, Garnham, Green, Ruzickova and Alda2008, Reference Hajek, Slaney, Garnham, Ruzickova, Passmore and Alda2005; McIntyre et al., Reference McIntyre, Danilewitz, Liauw, Kemp, Nguyen, Kahn and Taylor2010; Salvi et al., Reference Salvi, Salvo, Korčáková, Torriero, Aragno, Kolenič and Hajek2020). The associations between obesity and brain structure might help explain these links and provide new treatment options for some of these currently difficult to treat outcomes. Lastly, similar studies could help identify risk factors for neuroimaging outcomes, which may provide new opportunities for prevention or treatment of brain alterations with dietary/lifestyle medication or surgical interventions focused on weight management (Mansur et al., Reference Mansur, Zugman, Ahmed, Cha, Subramaniapillai, Lee and McIntyre2017; Mueller et al., Reference Mueller, Möller, Horstmann, Busse, Lepsien, Blüher and Pleger2015; Shan et al., Reference Shan, Li, Liu, Nie, Yin, Zhang and Shan2019; Tuulari et al., Reference Tuulari, Karlsson, Antikainen, Hirvonen, Pham, Salminen and Nummenmaa2016).
This study benefits from several unique advantages. The large sample size (2832 individuals) allowed us to test for interactions among relevant factors, which could not be conclusively studied in smaller, less powered studies. The multi-site nature of this study, with data from 17 sites in 13 countries, ensured highly generalizable representation of BD from around the world. We focused on overweight/obesity, which is a highly prevalent, but understudied factor in relation to brain structure in BD and which could also provide important insights into negative associations between obesitogenic medications and brain structure. In addition to novel findings on the interplay between BMI, BD, or medications and brain structure, we provide several replications of previous findings of associations between obesity or BD and specific regions of interest.
This study has the following limitations. The cross-sectional nature of our study does not allow us to discern the direction of the association, as brain alterations may predate or result from obesity. More detailed markers beyond BMI were not broadly available throughout the ENIGMA-BD Working Group. At the same time, BMI is much easier to acquire and is by far the most frequently used measure (García-García et al., Reference García-García, Michaud, Dadar, Zeighami, Neseliler, Collins and Dagher2019; Willette & Kapogiannis, Reference Willette and Kapogiannis2015), thus allowing for a more direct comparison with previous studies. Due to confidentiality reasons related to legacy datasets, we could not access raw, whole-brain data and could not use methods, such as voxel-based morphometry. Aside from using ENIGMA-standardized processing methods, we also addressed any differences between scanners statistically by using mixed models and including site as a random factor in all analyses. While there are other approaches, this is still by far the most utilized and accepted method for dealing with site effects (McWhinney et al., Reference McWhinney, Abé, Alda, Benedetti, Bøen and Del Mar Bonnin2021a; Thompson et al., Reference Thompson, Jahanshad, Ching, Salminen, Thomopoulos and Bright2020). Information about medications was limited to current usage at the time of scan. The study was not designed to comprehensively test the effects of medication, which would require a randomized controlled design. Therefore, the medication findings should be interpreted with caution, as medication prescriptions in clinical practice are not random. Also, medication details were limited to current prescription, without any measures of duration, dosage, compliance, previous medication exposure, treatment response or symptom levels at the time of prescription, so we cannot address the effects of these factors. The basic ENIGMA covariates, which are available across sites, did not contain any measures of cognitive/psychosocial functioning. Therefore, we could not evaluate the structure/function links. Fat content near the MRI coil may lead to slight signal intensity changes, but the vast majority of individuals in this study were normal weight to overweight. Psychiatric and other medical comorbidities, which might not be available for all the patients enrolled, may influence the interplay between BMI, BD, and neuroimaging findings. Finally, using other neuroimaging modalities could provide further insights into the mechanisms of the BMI effect. Last but not least, caution is needed when interpreting mediation analyses in observational studies.
To conclude, we confirmed consistent associations between higher BMI and lower cortical thickness across the cerebral mantle, in regions which were also associated with BD. There were few or no correlates of either condition with cortical surface area. In most regions number of medications remained associated with lower cortical thickness regardless of BMI, but there were also instances of mediation and moderation of associations between number of medication classes and cortical thickness by BMI. In terms of Li treatment, either the positive association between Li and cortical thickness was present regardless of BMI or people treated with Li showed steeper association between BMI and cortical thickness. All in all, BMI is important for understanding the neuroanatomical alterations in BD and the neurostructural correlates of psychiatric medications. We need prospective studies to investigate whether obesity is a modifiable/preventable risk factor for brain alterations in BD and whether the obesity-related negative psychiatric outcomes are related to obesity-related brain alterations.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291723000223.
Financial support
We gratefully acknowledge the following contributions and research funding sources that made this study possible: PT and CRKC of the Marina del Rey studies were supported by NIH grant U54 EB020403 from the Big Data to Knowledge (BD2K) Program; CRKC also acknowledges NIA T32AG058507 and partial research support from Biogen, Inc. (Boston, USA) for work unrelated to the topic of this manuscript. The St. Göran study was supported by grants from the Swedish Research Council (2018-02653), the Swedish foundation for Strategic Research (KF10-0039), the Swedish Brain foundation, and the Swedish Federal Government under the LUA/ALF agreement (ALF 20170019, ALFGBG-716801). This work is also part of the German multicenter consortium ‘Neurobiology of Affective Disorders. A translational perspective on brain structure and function’, funded by the German Research Foundation (Deutsche Forschungsgemeinschaft DFG; Forschungsgruppe/Research Unit FOR2107). Principal investigators (PIs) with respective areas of responsibility in the FOR2107 consortium are: Work Package WP1, FOR2107cohort and brain imaging: TK (speaker FOR2107; DFG grant numbers KI 588/14-1, KI 588/14-2), UD (co-speaker FOR2107; DA 1151/5-1, DA 1151/5-2), AK (KR 3822/5-1, KR 3822/7-2), IN (NE 2254/1-1 and NE 2254/2-1), CK (KO 4291/3-1). Further support from the German sites were provided by MNC and FOR2107-Muenster: This work was funded by the German Research Foundation (SFB-TRR58, Project C09 to UD) and the Interdisciplinary Center for Clinical Research (IZKF) of the medical faculty of Münster (grant Dan3/012/17 to UD and grant SEED11/18 to NO); FOR2107-Muenster: This work was supported by grants from the Interdisciplinary Center for Clinical Research (IZKF) of the medical faculty of Münster (grant MzH 3/020/20 to TH) and the German Research Foundation (DFG grants HA7070/2-2, HA7070/3, HA7070/4 to TH). The NUIG sample was supported by the Health Research Board (HRA_POR/2011/100). The Medellin studies (GIPSI) were supported by the PRISMA UNION TEMPORAL (UNIVERSIDAD DE ANTIOQUIA / HOSPITAL SAN VICENTE FUNDACIÓN), Colciencias-INVITACIÓN 990 de 3 de agosto de 2017, Codigo 99059634. The San Raffaele site was supported by the Italian Ministry of Health RF-2011-02350980 project. This research was also supported by the Irish Research Council (IRC) Postgraduate Scholarship, Ireland awarded to LN and to GM, and by the Health Research Board (HRA-POR-324) awarded to DMC. We thank the participants and the support of the Welcome-Trust HRB Clinical Research Facility and the Centre for Advanced Medical Imaging, St. James Hospital, Dublin, Ireland. The NUIG sample was supported by the Health Research Board (HRA_POR/2011/100). JS and RTK received support from the William K. Warren Foundation National Institute of Mental Health (R21MH113871); JS also acknowledges the National Institute of General Medical Sciences (P20GM121312). This study was also funded by EU-FP7-HEALTH-222963 ‘MOODIN- FLAME’ and EU-FP7-PEOPLE-286334 ‘PSYCHAID’. The Barcelona group would like to thank CIBERSAM (EPC) and the Instituto de Salud Carlos III (PI18/00877 and PI19/00394) for their support. This work was supported by the Singapore Bioimaging Consortium (RP C009/2006) research grant awarded to KS. The CIAM group (FMH – PI) was supported by the University Research Committee, University of Cape Town, and South African funding bodies National Research Foundation and Medical Research Council; DJS from CIAM was supported by the SAMRC. The Sydney studies were supported by the Australian National Health and Medical Research Council (NHMRC) Program Grant 1037196, Project Grants 1063960 and 1066177, the Lansdowne Foundation, Good Talk, and Keith Pettigrew Family; as well as the Janette Mary O'Neil Research Fellowship to JMF. The study was also supported by NIMH grant number: R01 MH090553(to RAO). Funding for the Oslo-Malt cohort was provided by the South Eastern Norway Regional Health Authority (2015-078), the Ebbe Frøland foundation, and a research grant from Throne-Holst. EV acknowledges the support of the Spanish Ministry of Science and Innovation (PI15/00283, PI18/00805) integrated into the Plan Nacional de I + D + I and co-financed by the ISCIII-Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER), the Instituto de Salud Carlos III, the CIBER of Mental Health (CIBERSAM), the Secretaria d'Universitats i Recerca del Departament d'Economia i Coneixement (2017 SGR 1365), the CERCA Programme, and the Departament de Salut de la Generalitat de Catalunya for the PERIS grant SLT006/17/00357. Lastly, this study was supported by the Canadian Institutes of Health Research (103703, 106469, 142255, 180449, and 487045), Nova Scotia Health Research Foundation, Dalhousie Clinical Research Scholarship to TH, Brain & Behavior Research Foundation (formerly NARSAD), and 2007 Young Investigator and 2015 Independent Investigator Awards to TH. PMT and CRKC received a grant from Biogen, Inc., for research unrelated to this manuscript. DJS has received research grants and/or consultancy honoraria from Lundbeck and Sun. LNY has received speaking/consulting fees and/or research grants from Abbvie, Alkermes, Allergan, AstraZeneca, CANMAT, CIHR, Dainippon Sumitomo Pharma, Janssen, Lundbeck, Otsuka, Sunovion, and Teva. TE received speaker's honoraria from Lundbeck and Janssen Cilag. EV has received grants and served as consultant, advisor, or CME speaker for the following entities (unrelated to the present work): AB-Biotics, Abbott, Allergan, Angelini, Dainippon Sumitomo Pharma, Ferrer, Gedeon Richter, Janssen, Lundbeck, Otsuka, Sage, Sanofi-Aventis, and Takeda. The authors declare no conflicts of interest.
Open access
