Network-Level Structural Abnormalities of Cerebral Cortex in Type 1 Diabetes Mellitus

Abstract

Type 1 diabetes mellitus (T1DM) usually begins in childhood and adolescence and causes lifelong damage to several major organs including the brain. Despite increasing evidence of T1DM-induced structural deficits in cortical regions implicated in higher cognitive and emotional functions, little is known whether and how the structural connectivity between these regions is altered in the T1DM brain. Using inter-regional covariance of cortical thickness measurements from high-resolution T1-weighted magnetic resonance data, we examined the topological organizations of cortical structural networks in 81 T1DM patients and 38 healthy subjects. We found a relative absence of hierarchically high-level hubs in the prefrontal lobe of T1DM patients, which suggests ineffective top-down control of the prefrontal cortex in T1DM. Furthermore, inter-network connections between the strategic/executive control system and systems subserving other cortical functions including language and mnemonic/emotional processing were also less integrated in T1DM patients than in healthy individuals. The current results provide structural evidence for T1DM-related dysfunctional cortical organization, which specifically underlie the top-down cognitive control of language, memory, and emotion.

About the Authors:

In Kyoon Lyoo

Contributed equally to this work with: In Kyoon Lyoo, Sujung Yoon

* E-mail: [email protected]

Affiliation: Ewha Brain Institute & College of Pharmacy, Graduate School of Pharmaceutical Sciences, Ewha W. University, Seoul, South Korea

Sujung Yoon

Contributed equally to this work with: In Kyoon Lyoo, Sujung Yoon

Affiliation: Department of Psychiatry, Catholic University of Korea College of Medicine, Seoul, South Korea

Perry F. Renshaw

Affiliation: The Brain Institute, The University of Utah, Salt Lake City, Utah, United States of America

Jaeuk Hwang

Affiliation: Department of Psychiatry, Soonchunhyang University College of Medicine, Seoul, South Korea

Sujin Bae

Affiliation: The Brain Institute, The University of Utah, Salt Lake City, Utah, United States of America

Gail Musen

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Psychiatry, Harvard Medical School, Boston, Massachusetts, United States of America

Jieun E. Kim

Affiliation: Department of Brain and Cognitive Sciences, Ewha W. University Graduate School, Seoul, South Korea

Nicolas Bolo

Affiliation: Department of Psychiatry, Harvard Medical School, Boston, Massachusetts, United States of America

Hyeonseok S. Jeong

Affiliation: Interdisciplinary Program in Brain Science, Seoul National University College of Natural Sciences, Seoul, South Korea

Donald C. Simonson

Affiliation: Department of Internal Medicine, Brigham and Women’s Hospital, Boston, Massachusetts, United States of America

Sun Hea Lee

Affiliation: Ewha Brain Institute & College of Pharmacy, Graduate School of Pharmaceutical Sciences, Ewha W. University, Seoul, South Korea

Katie Weinger

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Psychiatry, Harvard Medical School, Boston, Massachusetts, United States of America

Jiyoung J. Jung

Affiliation: Interdisciplinary Program in Brain Science, Seoul National University College of Natural Sciences, Seoul, South Korea

Christopher M. Ryan

Affiliation: Department of Psychiatry, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, United States of America

Yera Choi

Affiliation: Interdisciplinary Program in Brain Science, Seoul National University College of Natural Sciences, Seoul, South Korea

Alan M. Jacobson

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Psychiatry, Winthrop University Hospital, Mineola, New York, United States of America

Introduction

The central nervous system (CNS) serves as one of the representative end-organ targets of metabolic insult related to type 1 diabetes mellitus (T1DM) [1]. End-organ damage to the brain by T1DM involves changes in synaptic plasticity as well as direct damage to neurons [1]. Given the dynamic nature of the brain regional organization in response to neuronal damage [2]–[4], disease processes of diabetes may elicit rearrangement of the intra- and inter-regional brain architectures. Furthermore, because T1DM may begin in early childhood [5] when the brain undergoes profound developmental changes such as neuronal proliferation, synaptic pruning, and myelination [6]–[8], the potential for brain structural and functional reorganization may increase in response to metabolic insult by T1DM.

To render the human brain capable of exerting its characteristic functions, it is essential not only that individual neuroanatomical areas responsible for specific function are operational but also that the complex and hierarchal networks between corresponding areas are comparably efficient [9]. Recent advances in image analysis techniques have enabled us to identify and evaluate the complex organization of brain regional connections using graph theoretical network analysis [9], [10]. Specifically, structural network analysis using interregional covariance of cortical thickness measurements can provide important information on structural interactive connectivity of the human brain [11]–[14]. Disrupted structural network organization has increasingly been reported in aged population [14], as well as in patients with schizophrenia [15], dementia [16], [17], epilepsy [18], and multiple sclerosis [13]. However, to the best of our knowledge, no studies thus far have examined the inter-regional network efficiencies of the T1DM brain, despite well-recognized T1DM-related end-organ damages to the brain [1]. In the present study, we sought to determine the differences in topological attributes of underlying cortical structural networks between young adults with T1DM and healthy adults. To investigate broad network architectures with emphasis on their functional significance, we divided the cerebral cortex into the intrinsic structural sub-network systems that mainly subserve cortical functions including strategic/executive control, language, mnemonic/emotional processing, sensorimotor, and visual functions, based on prior knowledge of structural and functional modularity of regional connectivity [11], [19]. In light of the crucial role of the prefrontal regions in top-down cognitive control of other brain functions including language, mnemonic/emotional processing, and sensorimotor function [20]–[23], we expected reduced inter-network efficiencies between intrinsic structural sub-network systems in T1DM patients relative to control subjects.

Materials and Methods

Ethics statement

The study protocol was approved by the institutional committee on human subjects of the Joslin Diabetes Center, Boston, MA, USA, the Brain Imaging Center of the McLean Hospital, Boston, MA, USA, and the Seoul National University Hospital, Seoul, South Korea. All participants were given a complete explanation of the study details and then provided written informed consent to participate in the study.

Participants

Adult T1DM patients between the ages of 25 and 40 and age- and sex-matched healthy control subjects were recruited from the Joslin Diabetes Center. The sample has been described in detail elsewhere [24]. Among subjects of the original cohort, 81 T1DM patients and 38 healthy subjects who had never experienced a depressive episode were included in the present study.

Acquisition of high-resolution of MR images was performed at the Brain Imaging Center of the McLean Hospital. Subjects were excluded if they had major neurological or medical disorders, if they had ever been diagnosed with psychosis, schizophrenia, bipolar disorder, attention-deficit hyperactivity disorder, or cocaine, heroin, or alcohol dependence, as assessed using the Structural Clinical Interview for DSM-IV [25], or if they had any contraindications to magnetic resonance (MR) imaging. Serious diabetic complications including clinically significant diabetic nephropathy, painful or symptomatic neuropathy, proliferative retinopathy requiring laser treatment, or gastroparesis were additional exclusion criteria for the T1DM group. Diabetic characteristics were obtained from the medical records and laboratory tests. Lifetime average hemoglobin A1C (HbA1C) was measured using the average value of HbA1C levels, grouped and time-weighted every 4 years for entire duration of the disease [24], which indicates the level of lifetime glycemic control. Among 57 T1DM patients whose information regarding insulin therapy is available, 39 T1DM patients received rapid- or short-acting insulin analogues while 33 did intermediate- or long-acting insulin analogues. Several patients were treated with one or more classes of insulin analogues. Twelve patients were managed using the insulin pump. Demographic and clinical characteristics of study participants are demonstrated in Table 1.

[Figure omitted. See PDF.]

Table 1. Group characteristics of patients with T1DM and control subjects.

https://doi.org/10.1371/journal.pone.0071304.t001

MR image acquisition and measurement of cortical thickness

Brain MR images were obtained using a 1.5 Tesla GE whole body imaging system (Horizon LX, GE Medical systems, Milwaukee, Wisconsin, USA) at the Brain Imaging Center of the McLean Hospital. Coronal T1 weighted images were produced using a three–dimensional spoiled gradient echo pulse sequence (124 slices, slice thickness  = 1.5 mm, echo time [TE]  = 5 ms, repetition time [TR]  = 35 ms, matrix  = 256×192; field of view [FOV]  = 24 cm, flip angle [FA]  = 45°, number of excitations [NEX]  = 1). Axial T2 weighted images (TE  = 80 ms, TR = 3,000 ms, 256×192 matrix; FOV  = 20 cm, FA = 90°, NEX  = 0.5, slice thickness/gap  = 3/0 mm) and Fluid Attenuated Inversion Recovery axial images (TE  = 133 ms, TR  = 9,000 ms, inversion time  = 2,200 ms, matrix  = 256×192, FOV  = 20 cm, FA = 90°, NEX  = 1, slice thickness/gap  = 5/2 mm) were obtained to screen for potential brain structural abnormalities.

Cortical thickness measurements were conducted using the FreeSurfer Tools (http://surfer.nmr.mgh.harvard.edu/) and the method implemented in this study has been previously described in detail elsewhere [24], [26]–[28]. In brief, initial surface was obtained for further reconstructing the pial surface using deformable surface algorithm [29] after the segmentation of white matter [30], tessellation of gray/white matter boundaries [30], and automatic correction of topological defects [31]. All surfaces of each individual image data were visually inspected for accuracy of spatial registration and gray/white matter segmentation and were manually corrected, as recommended by the guideline. Cortical thickness at each vertex was defined as the shortest distances from the pial to the corresponding gray/white surface. This measurement for the thickness of the cerebral cortex has been validated using histological [32] and manual [32], [33] estimations with sub-millimeter precision. Reliability for this measurement has also been established across field strength, scanner upgrade, and manufacturer [34].

Inter-regional correlations and construction of structural cortical networks

Cerebral cortices on MR images were parcellated into 32 gyral-based regions for each hemisphere, using an automated labeling system that is distributed in the FreeSurfer Tool (Figure S1 in File S1) [35]. Cortical thickness of each vertex belonging to the specific parcellated region was averaged to define cortical thickness of that region.

Parcellated regions were further categorized into five segregated intrinsic structural network systems subserving different neurobehavioral functions including strategic/executive control, language, mnemonic/emotional processing, sensorimotor, and visual functions based on prior shared knowledge on structural and functional clusters (Figure S1 in File S1) [11], [19]. The structural network system responsible for strategic/executive control was composed of 14 parcellated regions mainly located in the dorsolateral prefrontal and cingulate cortices. The structural network system subserving speech and language-related function comprised 14 parcellated regions mainly located in the auditory cortex as well as Broca's and Wernicke's areas. The structural network system supporting mnemonic/emotional processing was composed of 18 parcellated regions that are primarily located in medial prefrontal and paralimbic cortices. The structural network system responsible for sensorimotor function was composed of 10 parcellated regions primarily placed in the primary or association cortices for sensory or motor function. The structural network system supporting visual function was composed of 8 parcellated regions which are located in the occipital cortex.

Partial correlation analyses of cortical thickness in each possible pair of 64 parcellated regions after adjusting for age and sex were used for the measurement of structural associations between regions. Interregional connection matrices were made for each of the T1DM and control groups by calculating the partial correlation coefficients for all subjects in each study group, representing the specific anatomical connections between each of the 2016 [ = (64×63)/2] possible pairs of parcellated regions. Each connection matrix was thresholded into a binarized matrix with a fixed sparsity (S), which was defined as the total number of edges (connections) in the graph divided by the maximum possible number of edges, 2016 [14], [16]. This thresholding method enabled us to exclusively examine the alterations of topological organization based on connection matrices comprising same number of edges in each group [16], [36], [37]. A wide range of sparsity levels (0.07≤S≤0.40) were applied to threshold the connection matrix with an incremental interval of 0.01 and the topological organizations of resultant matrices at each sparsity level were examined for both study groups [14], [16]. Brain network analyses were performed with an open source Matlab toolbox (http://www.brain-connectivity-toolbox.net) [10].

Measurement of global and local efficiencies of whole-brain structural cortical networks

Segregation and integration of whole-brain structural network reflect the efficiency of connections within densely interconnected regions and throughout remote brain regions, respectively [9]. The local efficiency (Eloc), a common measure of segregation of cortical network, specifically quantifies how efficiently the network transfers information at the clustering level [10]. The global efficiency (Eglob), a common measure of integration of cortical network, quantifies how efficiently the network exchanges the information at the global level [10].

Measurement of inter-network organizations of intrinsic structural network systems

The length of path reflecting potential routes of information between regions could represent an estimate for efficiency of network integration [10], [38]. The average shortest path length (L) between all pairs of nodes in the network has been commonly used to measure the efficiency of a network [10], [38]. Structural sub-network systems subserving strategic/executive control, language, mnemonic/emotional processing, and sensorimotor function, which are known to be related to T1DM pathology, were at the focus of the present study.

Based on evidence of the prefrontal regions' vital role in top-down control for other brain function [20]–[23], inter-network efficiencies between structural network system for strategic/executive control and other intrinsic network systems were also measured by calculating the average inverse shortest path length (1/L) with the connection submatrix for corresponding network systems: strategic/executive control and language, strategic/executive control and mnemonic/emotional processing, and strategic/executive control and sensorimotor function.

Network hub identification and analysis of hierarchical organization of whole-brain structural network

Among various measures of centrality [10], we used two measures to identify the regional hubs of structural network: degree (K) and betweenness-centrality (b). The degree of a specific region (i), Ki, is defined as the number of edges (connections) to other regions throughout the network [10], [12]. The betweenness centrality of a specific region (i), bi, is defined as the number of shortest paths between any two regions that pass this specific region [10]. For hub identification, degree and the betweenness-centrality values were calculated at a specific sparsity threshold of 0.23. This threshold level was selected to ensure that all cortical parcellated regions of both study groups were included in the brain network to minimize the number of false-positive connections [16]. The regions with the normalized betweenness-centrality (Bi) greater than 1.5 and the degree higher than the mean degree of whole-brain networks are considered as regional hubs of structural network [14], [16], [39]. The hierarchical organization has a feature characterized by having many hubs with high degree (total connectivity) but low clustering (local connectivity). The clustering coefficient (C) of a specific region (i), Ci, is defined as the ratio of the number of connections of that region with its nearest regions (neighbors) to the maximum number of possible connections in the network [10]. The hierarchical structure of the network was measured by the β coefficient of the power-law relationship between the clustering coefficient (C) and degree (k): C = k–β [15], [40], [41]. We further classified hubs with lower clustering, which have clustering coefficients lower than the mean values of whole-brain structural networks, as hubs at a hierarchically higher level [9], [15].

As supplementary analyses, modularity of the structural cortical network (m) was also examined [11]. The process of modularity optimization found the optimal community structure, which is a subdivision of the structural cortical network into the non-overlapping groups of nodes, for achieving the maximum network modularity. We also defined intra-module connectivity of each node as the z score of the total number of degrees of a specific node (i) within a module (Zi) [42], [43]. Inter-modular connectivity was measured by the participation coefficient (Pi) [42], [43].

Statistical analyses

Group differences in demographic characteristics were analyzed using the independent t-test or chi-square test.

Differences in topological parameters for structural cortical network between the T1DM and control groups were examined using a non-parametric permutation test [44]. The regional cortical thickness values for each subject were randomly assigned to one of two groups and the partial correlation matrices for each randomized group were recomputed. Binarized matrices were then obtained using the same sparsity threshold as in the real brain network. The network parameters were calculated at each sparsity and between-group differences in all parameters were also obtained. The randomization steps were repeated 1000 times to sample the permutation distribution of differences in network parameters and 95% points of each distribution were used as critical values for a one-tailed test of null hypothesis with a probability of type 1 error of 0.05.

Results

Global and local efficiencies of whole-brain structural networks

We first assessed the global and local efficiencies of whole-brain structural networks, and compare the results between the T1DM and control groups. The global efficiency values (Eglob) and the local efficiency (Eloc) of whole-brain structural networks of T1DM patients were low at a wide range of sparsity levels relative to control subjects (inset graphs in Figure 1). Statistical analyses demonstrated that the T1DM group had significantly lower Eglob than the control group in the sparsity range of 0.21≤S≤0.24. Statistically significant T1DM-related lower Eloc values were also observed in sparsity range of 0.16≤S≤0.23 and 0.37≤S≤0.40 (Figure 1).

[Figure omitted. See PDF.]

Figure 1. Between-group differences in global efficiency, local efficiency, and hierarchical organization of whole-brain structural networks.

The graphs showed the differences in network parameters between the T1DM and control groups (blue line). The mean values and 95% of confidence interval of the null distribution of between-group differences obtained from 1000 permutation tests at each sparsity level were represented as gray circles and error bars, respectively. Asterisks indicate significant differences in parameters between the T1DM and control groups at P<0.05. Inset graphs show the global efficiency (Eglob), local efficiency (Eloc), and hierarchical organization (β) of whole-brain structural networks in T1DM (red line) and control (gray line) subjects as functions of sparsity thresholds. Abbreviations: T1DM, type 1 diabetes mellitus.

https://doi.org/10.1371/journal.pone.0071304.g001

Inter-network efficiencies between intrinsic structural sub -network systems

In order to examine whether the hierarchical cognitive control of the prefrontal cortex was altered in T1DM patients, we investigated inter-network efficiencies between structural network systems for strategic/executive control and other major brain functions including language, mnemonic/emotional processing, and sensorimotor function.

Average inverse shortest path length (1/L) between intrinsic structural sub-network systems was compared between groups. Statistical analyses showed that average inverse shortest path length between the strategic/executive control system and language system was lower in T1DM patients than in control subjects at 0.08≤S≤0.09, 0.12≤S≤0.14, and 0.23≤S≤0.25 (Figure 2A). T1DM patients also showed lower values of average inverse shortest path length between the strategic/executive control system and mnemonic/emotional processing system relative to control subjects in the sparsity range of 0.11, 0.14, 0.16≤S≤0.24, and 0.37≤S≤0.40 (Figure 2B). However, no between-group difference was observed in average inverse shortest path length between the strategic/executive control system and sensorimotor system (Figure 2C).

[Figure omitted. See PDF.]

Figure 2. Connections and between-group differences in inter-network efficiencies.

Connections of cortical structural network for strategic/executive control with other networks mediating language (A), mnemonic/emotional processing (B), and sensorimotor function (C) and between-group differences in inter-network efficiencies are presented in the left and right panels, respectively. Brain templates in figures demonstrate cortical parcellated regions for corresponding intrinsic cortical structural sub-network systems and inter-network connections at the sparsity threshold of 0.23. Red arrows in graphs indicate the sparsity of 0.23 that whole-brain structural networks of both control and T1DM groups included all 64 connected brain regions. Hub regions shown in Figure 3 are marked as larger white circles with the radius in proportion of the value of Bi. The graphs showed the differences in average inverse shortest path length of submatrix for each corresponding intrinsic cortical structural network system between the T1DM and control groups (blue line). The mean values and 95% of confidence interval of the null distribution of between-group differences in parameters obtained from 1000 permutation tests at each sparsity level were represented as gray circles and error bars, respectively. Asterisks indicate significant differences in average inverse shortest path length between the T1DM and control groups at P<0.05. Inset graphs show average inverse shortest path length of submatrix representing inter-network efficiency, was plotted as a function of sparsity thresholds in T1DM (red line) and control (gray line) subjects. Abbreviations: T1DM, type 1 diabetes mellitus.

https://doi.org/10.1371/journal.pone.0071304.g002

Distributions of structural hubs and hierarchical organization of whole-brain structural networks

Sixty-four gyral-based parcellated regions were defined to generate a correlation matrix (Figure S1 in File S1). Graph theoretical properties of each parcellated cortical region including degree (Ki), clustering coefficient (Ci), and normalized betweenness-centrality (Bi) values for both T1DM and control groups are demonstrated in Figure S2 in File S1.

As shown in Figure 3, the regional distribution of hubs and inter-hub structural connections in T1DM patients was considerably different from that in control subjects. The most prominent between-group differences were observed in the prefrontal and paralimbic regions that are known to be responsible for strategic/executive function and mnemonic/emotional processing, respectively. In particular, while healthy individuals had several hubs at a hierarchically higher level in the dorsolateral and paralimbic regions (red circles of Figure 3A), T1DM patients did not show this regional pattern of hierarchically high-level hubs in the same regions. The dorsolateral prefrontal hubs of T1DM patients were highly clustering, and thus had a hierarchically lower level (blue circles of Figure 3B). In addition, the degree and betweenness-centrality measures of the prefrontal hubs were reduced in the T1DM group as compared with the control group (Figures S2 and S3 in File S1).

[Figure omitted. See PDF.]

Figure 3. Location of hubs and inter-hub structural connections in whole-brain structural networks.

Regions (brain templates of panels A for control and B T1DM subjects) in orange, green, blue, yellow, and light purple colors represent each intrinsic cortical structural sub-network system subserving strategic/executive control, language, mnemonic/emotional processing, sensorimotor, and visual functions, respectively. Figures depict hub regions and significant inter-hub structural connections of control (A) and T1DM (B) subjects at the sparsity threshold of 0.23. A given region was identified as a hub of whole-brain structural networks if its normalized betweenness-centrality (Bi) was greater than 1.5 and its degree (Ki) was above the network mean at the sparsity threshold of 0.23. Red circles denote the hub regions with low clustering (less than average clustering of whole-brain structural networks) indicating hubs at a higher position in the hierarchical organization. Blue circles denote the hub regions with high clustering (greater than average clustering of whole-brain structural networks) indicating hubs at a lower position in the hierarchical organization. The radius of circles is in proportion of the value of Bi of the region at sparsity threshold of 0.23. Abbreviations: T1DM, type 1 diabetes mellitus; R, right; L, left; SFC, superior frontal cortex; rMFC, rostral middle frontal cortex; cMFC, caudal middle frontal cortex; ICC, isthmus cingulate cortex; IFCc, inferior frontal cortex- pars opercularis; IPC, inferior parietal cortex; supM, supramarginal cortex; STC, superior temporal cortex; MTC, middle temporal cortex; rACC, rostral anterior cingulate cortex; ITC, inferior temporal cortex; Fus, fusiform cortex; postCen, postcentral cortex; preCun, precuneus cortex; Cun, cuneus cortex.

https://doi.org/10.1371/journal.pone.0071304.g003

Modular organization of structural cortical network in the T1DM and control groups was presented in Figure S4 in File S1. Similar findings for the group differences in regional distribution of nodes with higher inter-module connectivity were also observed in these modular analyses. Four connector hubs (defined as regions of Pi >0.05 and Zi >1.0) [44], which were located on the right dorsolateral prefrontal cortex, the right anterior cingulate cortex, and the left temporal cortex, were identified in the structural cortical network of the control group (Figure S5 in File S2). However, the connector hubs in the T1DM group were mainly distributed in the right temporal cortex not in the prefrontal cortex (Figure S5 in File S2).

Furthermore, the β values, which indicate the level of hierarchy in whole-brain structural networks, were also found to be larger in control subjects than in T1DM patients in a range of the sparsity (0.13≤S≤0.21, 0.24, 0.28) suggesting a more hierarchical organization of networks in the control group (Figure 1).

Discussion

Our data provide novel evidence of T1DM-related alterations in the topological organization of structural brain networks, which was assessed by using cortical thickness data from high-resolution T1-weighted MR imaging. We found that hierarchically high-level prefrontal hubs were relatively absent in T1DM patients, and that the levels of the hierarchical organization of whole-brain structural networks were also lower in T1DM patients than in control subjects. By utilizing the division of brain regions into segregated intrinsic sub-network systems with functional significance, we also found that the structural connections between the prefrontal regions exerting top-down cognitive controls and other regions subserving complex brain functions including language and mnemonic/emotional processing were likely to be less integrated in T1DM patients relative to control subjects.

Previous studies have reported that the normal human brain demonstrates significant levels of hierarchical organization, as it is characterized by optimal connections for top-down control with minimum wiring costs [40], [41], both in anatomical and functional networks [9], [15], [41]. Our finding that control subjects had hierarchically higher positioned hubs predominantly in the prefrontal regions is consistent with previous literature [15]. In light of the prefrontal regional deficits frequently observed in T1DM patients [45]–[50] along with the dorsolateral prefrontal thickness reductions in the present T1DM sample (Figure S6 in File S2) [24], loss of hierarchical organization in these regions could be anticipated in T1DM patients.

The absence of hierarchically organized hubs in key regions, such as the prefrontal regions, suggests that the contribution of these regions to the network efficiency may be reduced in T1DM patients. Our inter-network analysis results (Figure 2) corroborate this speculation. Along with the similar profile of global and local efficiencies of whole brain networks between the normal and T1DM brains, this finding also suggests that T1DM-related network disorganization, which would reflect inefficient connections between brain areas, may occur in a region-specific manner.

A recent functional study, which assessed resting-state brain networks in young adults with T1DM, has demonstrated similar findings to ours [51]. The connectivity of networks related to ventral attention, language, and fronto-parietal circuits was decreased in T1DM patients relative to control subjects. In contrast, the connectivity of sensorimotor and visual networks appeared to be preserved in T1DM patients.

Although sample characteristics (old age vs. young age) and measurements for network construction (white matter connectivity vs. cortical thickness) might be different from ours, it should be noted that the structural white matter network for the type 2 diabetes mellitus (T2DM) brain was also altered in terms of decreased global and local efficiency [52]. Furthermore, T2DM-related whole-brain white matter network abnormalities were associated with slowing of information process [52]. Future studies would be necessary to examine the neurocognitive as well as clinical relevance of T1DM-related cortical network disorganization for example, the altered top-down control role of the prefrontal cortex.

As T1DM often starts during childhood and adolescence [5], the current findings could be understood from a developmental perspective. Emerging evidence suggests that regions that subserve more primary functions including sensory or motor processing typically mature in infancy and early childhood [53]. However, brain regions responsible for more complex functions may require additional intrinsic and extrinsic inputs for maturation, and thereby mature at much later ages from childhood through adolescence [6], [8], [53]–[55]. In the context of this heterochronous developmental trajectory, T1DM and related metabolic insults may typically influence the brain regions predominantly implicated in complex cognitive functions rather than more primary sensorimotor functions. Moreover, childhood and adolescence is a critical period during which patients may be more vulnerable to diabetes-related metabolic disturbances than adulthood [5].

Considering a recent longitudinal observation suggesting minimal age-related brain volume loss in adolescents with T1DM [56], the current alterations in structural connections of the frontal cortex with other brain regions may reflect the effects of T1DM-related metabolic insult on the brain over a longer period of time beyond adolescence. For example, several cellular mechanisms including hyperglycemia-induced mitochondrial dysfunction and changes in synaptic plasticity [57], [58] could also have interactions with developmental reorganization of cortex in response to T1DM.

It is unclear why the structural connections subserving prefrontal controls on sensorimotor cortices were likely to be retained relative to those on other functions. Subclinical involvements of the retina or peripheral nerves in T1DM, which could occur even in the early stage of disease [59], [60], might precipitate the brain activities to compensate for the altered information conveyed by the peripheral nervous systems [51].

As suggested by previous diffusion tensor imaging studies that could provide anatomical information on the connections between brain regions, white matter connectivity in a wide range of brain areas may be affected in T1DM patients relative to control subjects [61]–[63]. Although findings of specific regional vulnerability of white matter to T1DM-related metabolic insult may be rather inconsistent, partly because of age span of study participants as well as implemented methodologies, studies demonstrated not only fronto-temporal but also parieto-occipital regional involvements related to T1DM [61]–[63]. The present with a sample of young adults with T1DM provides the descriptive findings of relatively preserved hubs in the parieto-occipital regions in T1DM. The functional relevance of regionally altered hierarchy of the T1DM brain should be further investigated in future studies using multimodal imaging analyses.

T1DM may lead to several neurobehavioral abnormalities, as the CNS manifestations of T1DM include impaired performances on several cognitive domains as well as a wide range of emotional problems such as depression and anxiety [5], [64]–[66]. An increasing number of neuroimaging studies suggest that the frontal and temporal brain regions may be the neuroanatomical correlates of these neurobehavioral abnormalities [24], [45]–[50], [67]–[69], in that these regions, which are known to have high insulin-receptor density and require high levels of cerebral energy, are particularly vulnerable to abnormal brain glucose metabolism [70]–[72].

The current findings suggest that altered inter-network connectivity between these key regions may play an important role in T1DM-associated neurobehavioral abnormalities, in addition to T1DM-related volume or thickness reductions with a similar regional pattern. In contrast, network efficiencies between brain regions implicated in the earliest developed functions such as sensorimotor function were preserved in T1DM patients.

One might assume that the metabolic load including history of glycemic control might affect the extent of alterations of structural cortical network in T1DM patients as shown in dose-dependent effects of hyperglycemia on neurochemical disturbances in the prefrontal cortex of T1DM patients [46]. Although our subgroup analyses using a non-parametric permutation test could not find the differences in global and local efficiencies of whole-brain structural cortical network between the good and poor glycemic control T1DM subgroups (Figure S7 in File S2) partly because of a relatively small sample size of each subgroup, the potential effects of disease load including glycemic control levels should be examined in the future studies with a sufficient sample size of T1DM. Furthermore, a longitudinal follow up of the T1DM sample to evaluate the progression of illness would also be important in providing valuable information on the dose-dependent relationship between T1DM-related clinical characteristics and alterations in structural cortical network.

Several limitations should be taken into account in interpreting the results. First, it remains to be determined whether the inter-regional correlation patterns of cortical thickness could reflect the alterations in intrinsic network systems for specific brain functions. However, growing evidence suggests that shared morphological organization may occur in functionally synchronous brain regions due to mutually trophic or neuroplastic influences [15], [19]. Recent structural imaging studies using diffusion tensor imaging also validated that the correlations of cortical thickness could indirectly indicate the anatomical connectivity of white matter tracts between regions [73]. For example, strong inter-regional correlations of cortical thickness were observed particularly where regions are axonally connected [73].

It should also be noted that anatomical connection matrices based on cortical thickness data that were used in the present study, were constructed by measuring inter-regional correlation coefficients over a group of subjects. In contrast to individual levels of connection matrices yielded from functional imaging data, we could not determine the network information at individual levels as well as the effects of inter-individual clinical variances on network organization. Furthermore, future studies using different cortical parcellation maps would be needed to replicate the current results based on Desikan-Killiany cortical atlas [35].

Although evidence based on histological and neuroimaging studies has indicated that cortical thickness may reflect composite estimations of size, density, and arrangement of cortical neurons, glia, or nerve fibers [74]–[76], a pending consensus on cellular substrates of changes in cortical thickness should also be considered in inferring the altered patterns of the anatomical connections matrix based on cortical thickness.

An Increasing body of literature has provided evidence that topological changes in brain networks may reflect disease progression or plastic changes in response to external stimuli [2]–[4]. Our data indicate that T1DM, occurring in childhood, may influence the normal development of topological organizations specifically responsible for higher complex functions such as strategic/executive control and mnemonic/emotional processing. Specifically, structural organization underlying top-down cognitive control of language, memory, and emotion by the prefrontal cortex may be affected in T1DM patients. These findings of T1DM-related network disorganization also provide a clue to the brain-based mechanisms of common long-term T1DM complications including cognitive decline or emotional disturbances.

Supporting Information

[Figure omitted. See PDF.]

File S1.

Supporting figures. Figure S1. Cortical parcellation maps for construction of cortical structural network. The entire cerebral cortex was segmented into 64 cortical regions according to a prior surface parcellation based on the Desikan-Killiany Atlas for the construction of cortical structural network correlation matrix. Regions in orange, green, blue, yellow, and light purple colors represent each cortical structural sub-network exerting strategic/executive control, language, mnemonic/emotional processing, sensorimotor, and visual functions, respectively. Abbreviations: R, right; L, left; SFC, superior frontal cortex; rMFC, rostral middle frontal cortex; cMFC, caudal middle frontal cortex; IFCa, inferior frontal cortex- pars orbitalis; IFCb, inferior frontal cortex- pars triangularis; IFCc, inferior frontal cortex, pars opercularis; LOFC, lateral orbitofrontal cortex; MOFC, medial orbitofrontal cortex; preCen, precentral cortex; paraCen, paracentral cortex; Fpol, frontal pole; STC, superior temporal cortex; MTC, middle temporal cortex; ITC, inferior temporal cortex; TTC, transverse temporal cortex; Fus, fusiform cortex; Ento, entorhinal cortex; paraHi, parahippocampal cortex; Tpol, temporal pole; SPC, superior parietal cortex; IPC, inferior parietal cortex; supM, supramarginal cortex; postCen, postcentral cortex; preCun, precuneus cortex; rACC, rostral anterior cingulate cortex; cACC, caudal anterior cingulate cortex; PCC, posterior cingulate cortex; ICC, isthmus cingulate cortex; LOC, lateral occipital cortex; Cun, Cuneus; PeriCal, pericalcarine cortex; Lin, lingual cortex. Figure S2. Graph theoretical properties of parcellated cortical regions. Major brain functions subserved by parcellated regions are represented in distinct colors (orange, strategic/executive control; green, language function; blue, mnemonic/emotional processing; yellow, sensorimotor function; light purple, visual function, respectively). Abbreviations corresponding to the names of parcellated cortical regions are demonstrated in Figure S1. A. Comparison in the degree (Ki) of each parcellated cortical region between control and T1DM subjects. The length of each bar indicates the degree of each parcellated cortical region. Degree is defined as the number of edges (connections) to other regions throughout the network at the sparsity threshold of 0.23. Regions were sorted by degree in descending order for each group. The values for degree of hierarchically high hub regions are presented as red bars while those of hierarchically low hub regions are presented as blue bars. Note that the high-degree regions of control subjects were likely to be at a hierarchically high level (low-clustering). In contrast, the high-degree regions of T1DM subjects were likely to be located at a hierarchically lower level (high-clustering). B. Comparison in the clustering coefficient (Ci) of each parcellated cortical region between control and T1DM subjects. The length of each bar indicates the clustering coefficient of each parcellated cortical region. The clustering coefficient is defined as the ratio of the number of connections of that region with its nearest regions (neighbors) to the maximum number of possible connections in the network at the sparsity threshold of 0.23. Regions were sorted by clustering coefficient in descending order for each group. The values for degree of hierarchically high hub regions are presented as red bars while those of hierarchically low hub regions are presented as blue bars. C. Comparison in the normalized betweenness-centrality (Bi) of each parcellated cortical region between control and T1DM subjects. The length of each bar indicates the normalized betweenness-centrality of each parcellated cortical region. Betweenness-centrality is defined as the number of shortest paths between any two regions that pass this specific region at the sparsity threshold of 0.23. Regions were sorted by normalized betweenness-centrality in descending order for each group. The values for degree of hierarchically high hub regions are presented as red bars while those of hierarchically low hub regions are presented as blue bars. Note that the regions with high centrality were likely to be at a hierarchically high level (low-clustering) in control subjects. In contrast, the regions with high centrality were likely to be located at a hierarchically lower level (high-clustering) in T1DM patients. Figure S3. Topological properties of hubs in whole-brain structural networks of control (A) and T1DM subjects (B). Bar graphs depict the degrees (Ki) and clustering coefficients (Ci) of each hub region in control (left) and T1DM subjects (right) at sparsity threshold of 0.23. Red dotted lines of bar graphs indicate the mean clustering coefficients of whole-brain network of control (Cp = 0.55, left) and T1DM (Cp  = 0.49, right) subjects. Abbreviations: T1DM, type 1 diabetes mellitus; R, right; L, left; SFC, superior frontal cortex; rMFC, rostral middle frontal cortex; cMFC, caudal middle frontal cortex; ICC, isthmus cingulate cortex; IFCc, inferior frontal cortex- pars opercularis; IPC, inferior parietal cortex; supM, supramarginal cortex; STC, superior temporal cortex; MTC, middle temporal cortex; rACC, rostral anterior cingulate cortex; ITC, inferior temporal cortex; Fus, fusiform cortex; postCen, postcentral cortex; preCun, precuneus cortex; Cun, cuneus cortex. Figure S4. Modular organization of the structural cortical network of the T1DM and control groups. The binarized matrix at a sparsity threshold of 0.23 for the control and T1DM groups were presented. Four modules were detected in the control group (left column). Module I contains 22 cortical regions including RSFC, RcMFC, LcMFC, RpreCen, LpreCen, RparaCen, LparaCen, RSPC, LSPC, LIPC, RsupraM, RpostCen, LpostCen, RpreCun, LpreCun, LcACC, RSTC, LSTC, RTCC, LEnto, LTPol, and RCun (yellow regions in the left column). Module II contains 22 cortical regions including RrMFC, LIFCc, RIFCb, RIFCa, RLOFC, RMOFC, LMOFC, RFpol, LFPol, RIPC, RrACC, LrACC, RcACC, RPCC, LPCC, RICC, LICC, RFus, LTTC, RparaHi, LparaHi, and RLOC (green regions in the left column). Module III contains 14 cortical regions including LSFC, LrMFC, RIFCc, LIFCb, LIFCa, LLOFC, LsupraM, RMTC, LMTC, RITC, LITC, LFus, REnto, and RTpol (blue regions in the left column). Module IV contains 6 cortical regions including LLOC, Lcun, RperiCal, LperiCal, RLin, and LLin (light yellow regions in the left column). The T1DM structural cortical network consists of 4 modules (right column). Module I contains 20 cortical regions including RSFC, LSFC, RrMFC, LrMFC, RcMFC, LcMFC, RIFCb, LIFCb, RpreCen, LpreCen, RparaCen, LparaCen, RSPC, LSPC, RIPC, LIPC, RsupraM, LsupraM, LpostCen, and RpreCun (yellow regions in the right column). Module II contains 20 cortical regions including LIFCc, LIFCa, RLOFC, LLOFC, RMOFC, LMOFC, RFpol, LFpol, RrACC, LrACC, RcACC, LcACC, RPCC, LPCC, RICC, LICC, RTPol, RparaHi, LparaHi, and RCun (green regions in the right column). Module III contains 17 cortical regions including RIFCc, RIFCa, LpreCun, RSTC, LSTC, RMTC, LMTC, RITC, LITC, RFus, LFus, RTTC, LTTC, REnto, LEnto, LTPol, and RLOC (blue regions in the right column). Module IV contains 7 cortical regions including RpostCen, LLOC, LCun, RperiCal, LperiCal, RLin, and LLin (light yellow regions in the right column). Note that the control group had more dorsolateral prefrontal regions of higher inter-module connectivity (Pi >0.6)(red circles, LSFC, LrMFC, RrMFC) than the T1DM group (red circle. LrMFC). See also Figure S5. Abbreviations: T1DM, type 1 diabetes mellitus. Abbreviations corresponding to the names of parcellated cortical regions are demonstrated in Figure S1.

https://doi.org/10.1371/journal.pone.0071304.s001

(DOC)

File S2.

Supporting figures. Figure S5. Inter-module connectivity (participation coefficient, Pi) and intra-module connectivity (z scores of within-module degree, Zi) of each cortical region for the control (left column) and T1DM (right column) groups. Regions with blue arrows indicate those with participation coefficients greater than 0.6 suggesting the strong inter-module connectivity. Regions with red arrows indicate those with z scores of within-module degree greater than 1. The node with Pi >0.05 was defined as a connector node while that with Pi ≤0.05 was defined as a provincial node (Meunier et al, 2009). The node with Zi >1 was defined as a hub (Meunier et al, 2009). According to this criteria, 4 connector hubs including RrMFC, RcMFC, RrACC, and LITC were identified in the structural cortical network of the control group. T1DM cortical network also had 4 connector hubs which were located on RrACC, RSTC, RMTC, and RFus. Abbreviations: T1DM, type 1 diabetes mellitus.Abbreviations corresponding to the names of parcellated cortical regions are demonstrated in Figure S1. Figure S6. Group difference in cortical thickness of each parcellated region between control and T1DM subjects. Group differences in cortical thickness of each parcellated region are represented as z scores using group mean thickness and standard deviation values of the control group (gray bars). Regions were sorted by the z scores of thickness difference in ascending order. A negative z score indicates the number of standard deviation units thinner than expected thickness of control subjects. A reduction in dorsolateral prefrontal cortical thickness (left superior frontal, left caudal middle frontal, and right superior cortices, regions with z scores less than −0.70) was observed in T1DM patients relative to control subjects. Major brain functions subserved by each parcellated regions are represented in separate colors (orange, strategic/executive control; green, language function; blue, mnemonic/emotional processing; yellow, sensorimotor function; light purple, visual function, respectively). Abbreviations corresponding to the name of parcellated cortical regions are demonstrated in Figure S1. Figure S7. Global efficiency and local efficiency of whole-brain structural networks in the T1DM subgroups according to lifetime glycemic control levels. A non-parametric permutation test with 1000 repetitions showed no statistically significant differences in mean values of Eglob as well as local efficiency (Eloc) across sparsity between the good and poor glycemic control T1DM groups. 1Lifetime average HbA1C level less than 7% (n = 18, mean lifetime average HbA1C = 6.55%). 2Lifetime average HbA1C level 7% or greater (n = 63, mean lifetime average HbA1C = 8.40%).

https://doi.org/10.1371/journal.pone.0071304.s002

(DOC)

Author Contributions

Conceived and designed the experiments: IKL PFR GM AMJ. Performed the experiments: IKL PFR GM NB KW AMJ. Analyzed the data: IKL SY PFR JH SB GM JEK HSJ SHL KW JJJ YC AMJ. Contributed reagents/materials/analysis tools: IKL PFR AMJ. Wrote the paper: IKL SY PFR JH SB GM JEK NB HSJ DCS SHL KW JJJ CMR YC AMJ. Supervised study: IKL PFR AMJ. Critically revised the manuscript: IKL SY PFR JH SB GM JEK NB HSJ DCS SHL KW JJJ CMR YC AMJ.

Citation: Lyoo IK, Yoon S, Renshaw PF, Hwang J, Bae S, Musen G, et al. (2013) Network-Level Structural Abnormalities of Cerebral Cortex in Type 1 Diabetes Mellitus. PLoS ONE 8(8): e71304. https://doi.org/10.1371/journal.pone.0071304

References

1. Klein JP, Waxman SG (2003) The brain in diabetes: molecular changes in neurons and their implications for end-organ damage. Lancet Neurol 2: 548–554.

2. Ghosh A, Sydekum E, Haiss F, Peduzzi S, Zorner B, et al. (2009) Functional and anatomical reorganization of the sensory-motor cortex after incomplete spinal cord injury in adult rats. J Neurosci 29: 12210–12219.

3. Shu N, Liu Y, Li J, Li Y, Yu C, et al. (2009) Altered anatomical network in early blindness revealed by diffusion tensor tractography. PLoS One 4: e7228.

4. Yogarajah M, Focke NK, Bonelli SB, Thompson P, Vollmar C, et al. (2010) The structural plasticity of white matter networks following anterior temporal lobe resection. Brain 133: 2348–2364.

5. Biessels GJ, Deary IJ, Ryan CM (2008) Cognition and diabetes: a lifespan perspective. Lancet Neurol 7: 184–190.

6. Giedd JN, Rapoport JL (2010) Structural MRI of pediatric brain development: what have we learned and where are we going? Neuron 67: 728–734.

7. Power JD, Fair DA, Schlaggar BL, Petersen SE (2010) The development of human functional brain networks. Neuron 67: 735–748.

8. Toga AW, Thompson PM, Sowell ER (2006) Mapping brain maturation. Trends Neurosci 29: 148–159.

9. Bullmore E, Sporns O (2009) Complex brain networks: graph theoretical analysis of structural and functional systems. Nat Rev Neurosci 10: 186–198.

10. Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52: 1059–1069.

11. Chen ZJ, He Y, Rosa-Neto P, Germann J, Evans AC (2008) Revealing modular architecture of human brain structural networks by using cortical thickness from MRI. Cereb Cortex 18: 2374–2381.

12. He Y, Chen ZJ, Evans AC (2007) Small-world anatomical networks in the human brain revealed by cortical thickness from MRI. Cereb Cortex 17: 2407–2419.

13. He Y, Dagher A, Chen Z, Charil A, Zijdenbos A, et al. (2009) Impaired small-world efficiency in structural cortical networks in multiple sclerosis associated with white matter lesion load. Brain 132: 3366–3379.

14. Zhu W, Wen W, He Y, Xia A, Anstey KJ, et al. (2012) Changing topological patterns in normal aging using large-scale structural networks. Neurobiol Aging 33: 899–913.

15. Bassett DS, Bullmore E, Verchinski BA, Mattay VS, Weinberger DR, et al. (2008) Hierarchical organization of human cortical networks in health and schizophrenia. J Neurosci 28: 9239–9248.

16. He Y, Chen Z, Evans A (2008) Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer's disease. J Neurosci 28: 4756–4766.

17. Seeley WW, Crawford RK, Zhou J, Miller BL, Greicius MD (2009) Neurodegenerative diseases target large-scale human brain networks. Neuron 62: 42–52.

18. Raj A, Mueller SG, Young K, Laxer KD, Weiner M (2010) Network-level analysis of cortical thickness of the epileptic brain. Neuroimage 52: 1302–1313.

19. Zielinski BA, Gennatas ED, Zhou J, Seeley WW (2010) Network-level structural covariance in the developing brain. Proc Natl Acad Sci U S A 107: 18191–18196.

20. Miyashita Y (2004) Cognitive memory: cellular and network machineries and their top–down control. Science 306: 435–440.

21. Miller BT, D'Esposito M (2005) Searching for “the top” in top-down control. Neuron 48: 535–538.

22. Del Arco A, Mora F (2009) Neurotransmitters and prefrontal cortex-limbic system interactions: implications for plasticity and psychiatric disorders. J Neural Transm 116: 941–952.

23. Kouneiher F, Charron S, Koechlin E (2009) Motivation and cognitive control in the human prefrontal cortex. Nat Neurosci 12: 939–945.

24. Lyoo IK, Yoon S, Jacobson AM, Hwang J, Musen G, et al. (2012) Prefrontal cortical deficits in type 1 diabetes mellitus: brain correlates of comorbid depression. Arch Gen Psychiatry 69: 1267–1276.

25. First MB, Gibbon M, Spitzer RL, Williams JBW (1996) User's Guide for the Structured Clinical Interview for DSM-IV Axis I Disorders, Research Version 2.0. New York: Biometrics Research Department, New York State Psychiatric Institute.

26. Lyoo IK, Sung YH, Dager SR, Friedman SD, Lee JY, et al. (2006) Regional cerebral cortical thinning in bipolar disorder. Bipolar Disord 8: 65–74.

27. Lyoo IK, Kim JE, Yoon SJ, Hwang J, Bae S, et al. (2011) The neurobiological role of the dorsolateral prefrontal cortex in recovery from trauma. Longitudinal brain imaging study among survivors of the South Korean subway disaster. Arch Gen Psychiatry 68: 701–713.

28. Mollica RF, Lyoo IK, Chernoff MC, Bui HX, Lavelle J, et al. (2009) Brain structural abnormalities and mental health sequelae in South Vietnamese ex-political detainees who survived traumatic head injury and torture. Arch Gen Psychiatry 66: 1221–1232.

29. Fischl B, Dale AM (2000) Measuring the thickness of the human cerebral cortex from magnetic resonance images. Proc Natl Acad Sci U S A 97: 11050–11055.

30. Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9: 179–194.

31. Fischl B, Liu A, Dale AM (2001) Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. IEEE Trans Med Imaging 20: 70–80.

32. Rosas HD, Liu AK, Hersch S, Glessner M, Ferrante RJ, et al. (2002) Regional and progressive thinning of the cortical ribbon in Huntington's disease. Neurology 58: 695–701.

33. Kuperberg GR, Broome MR, McGuire PK, David AS, Eddy M, et al. (2003) Regionally localized thinning of the cerebral cortex in schizophrenia. Arch Gen Psychiatry 60: 878–888.

34. Han X, Jovicich J, Salat D, van der Kouwe A, Quinn B, et al. (2006) Reliability of MRI-derived measurements of human cerebral cortical thickness: the effects of field strength, scanner upgrade and manufacturer. Neuroimage 32: 180–194.

35. Desikan RS, Segonne F, Fischl B, Quinn BT, Dickerson BC, et al. (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31: 968–980.

36. Achard S, Bullmore E (2007) Efficiency and cost of economical brain functional networks. PLoS Comput Biol 3: e17.

37. Stam CJ, Jones BF, Nolte G, Breakspear M, Scheltens P (2007) Small-world networks and functional connectivity in Alzheimer's disease. Cereb Cortex 17: 92–99.

38. Watts DJ, Strogatz SH (1998) Collective dynamics of 'small-world' networks. Nature 393: 440–442.

39. Yan C, Gong G, Wang J, Wang D, Liu Det al (2011) Sex- and Brain Size-Related Small-World Structural Cortical Networks in Young Adults: A DTI Tractography Study. Cereb Cortex 17: 92–99.

40. Ravasz E, Barabasi AL (2003) Hierarchical organization in complex networks. Phys Rev E Stat Nonlin Soft Matter Phys 67: 026112.

41. Supekar K, Musen M, Menon V (2009) Development of large-scale functional brain networks in children. PLoS Biol 7: e1000157.

42. Chen ZJ, He Y, Rosa-Neto P, Gong G, Evans AC (2011) Age-related alterations in the modular organization of structural cortical network by using cortical thickness from MRI. Neuroimage 56: 235–245.

43. Meunier D, Achard S, Morcom A, Bullmore E (2009) Age-related changes in modular organization of human brain functional networks. Neuroimage 44: 715–723.

44. Bullmore ET, Suckling J, Overmeyer S, Rabe-Hesketh S, Taylor E, et al. (1999) Global, voxel, and cluster tests, by theory and permutation, for a difference between two groups of structural MR images of the brain. IEEE Trans Med Imaging 18: 32–42.

45. Musen G (2008) Cognition and brain imaging in type 1 diabetes. Curr Diab Rep 8: 132–137.

46. Lyoo IK, Yoon SJ, Musen G, Simonson DC, Weinger K, et al. (2009) Altered prefrontal glutamate-glutamine-gamma-aminobutyric acid levels and relation to low cognitive performance and depressive symptoms in type 1 diabetes mellitus. Arch Gen Psychiatry 66: 878–887.

47. Pell GS, Lin A, Wellard RM, Werther GA, Cameron FJ, et al. (2012) Age-related loss of brain volume and T2 relaxation time in youth with type 1 diabetes. Diabetes Care 35: 513–519.

48. Northam EA, Rankins D, Lin A, Wellard RM, Pell GS, et al. (2009) Central nervous system function in youth with type 1 diabetes 12 years after disease onset. Diabetes Care 32: 445–450.

49. Wessels AM, Simsek S, Remijnse PL, Veltman DJ, Biessels GJ, et al. (2006) Voxel-based morphometry demonstrates reduced grey matter density on brain MRI in patients with diabetic retinopathy. Diabetologia 49: 2474–2480.

50. Ferguson SC, Blane A, Perros P, McCrimmon RJ, Best JJet al (2003) Cognitive ability and brain structure in type 1 diabetes: relation to microangiopathy and preceding severe hypoglycemia. Diabetes 52: 149–156.

51. van Duinkerken E, Schoonheim MM, Sanz-Arigita EJ, IJzerman RG, Moll AC, et al. (2012) Resting-state brain networks in type 1 diabetic patients with and without microangiopathy and their relation to cognitive functions and disease variables. Diabetes 61: 1814–1821.

52. Reijmer YD, Leemans A, Brundel M, Kappelle LJ, Biessels GJ (2013) Disruption of the cerebral white matter network is related to slowing of information processing speed in patients with type 2 diabetes. Diabetes 62: 2112–2115.

53. Rapoport JL, Gogtay N (2008) Brain neuroplasticity in healthy, hyperactive and psychotic children: insights from neuroimaging. Neuropsychopharmacology 33: 181–197.

54. Luna B, Sweeney JA (2004) The emergence of collaborative brain function: FMRI studies of the development of response inhibition. Ann N Y Acad Sci 1021: 296–309.

55. Luna B, Thulborn KR, Munoz DP, Merriam EP, Garver KE, et al. (2001) Maturation of widely distributed brain function subserves cognitive development. Neuroimage 13: 786–793.

56. Perantie DC, Koller JM, Weaver PM, Lugar HM, Black KJ, et al. (2011) Prospectively determined impact of type 1 diabetes on brain volume during development. Diabetes 60: 3006–3014.

57. Brownlee M (2005) The pathobiology of diabetic complications: a unifying mechanism. Diabetes 54: 1615–1625.

58. Kamal A, Biessels GJ, Duis SE, Gispen WH (2000) Learning and hippocampal synaptic plasticity in streptozotocin-diabetic rats: interaction of diabetes and ageing. Diabetologia 43: 500–506.

59. van Dijk HW, Verbraak FD, Kok PH, Garvin MK, Sonka M, et al. (2010) Decreased retinal ganglion cell layer thickness in patients with type 1 diabetes. Investigative ophthalmology & visual science 51: 3660–3665.

60. Almeida S, Riddell MC, Cafarelli E (2008) Slower conduction velocity and motor unit discharge frequency are associated with muscle fatigue during isometric exercise in type 1 diabetes mellitus. Muscle & nerve 37: 231–240.

61. Antenor-Dorsey JA, Meyer E, Rutlin J, Perantie DC, White NH, et al. (2013) White matter microstructural integrity in youth with type 1 diabetes. Diabetes 62: 581–589.

62. Franc DT, Kodl CT, Mueller BA, Muetzel RL, Lim KO, et al. (2011) High connectivity between reduced cortical thickness and disrupted white matter tracts in long-standing type 1 diabetes. Diabetes 60: 315–319.

63. Kodl CT, Franc DT, Rao JP, Anderson FS, Thomas W, et al. (2008) Diffusion tensor imaging identifies deficits in white matter microstructure in subjects with type 1 diabetes that correlate with reduced neurocognitive function. Diabetes 57: 3083–3089.

64. Barnard KD, Skinner TC, Peveler R (2006) The prevalence of co-morbid depression in adults with Type 1 diabetes: systematic literature review. Diabet Med 23: 445–448.

65. Brands AM, Biessels GJ, de Haan EH, Kappelle LJ, Kessels RP (2005) The effects of type 1 diabetes on cognitive performance: a meta-analysis. Diabetes Care 28: 726–735.

66. McIntyre RS, Kenna HA, Nguyen HT, Law CW, Sultan F, et al. (2010) Brain volume abnormalities and neurocognitive deficits in diabetes mellitus: points of pathophysiological commonality with mood disorders? Adv Ther 27: 63–80.

67. Lobnig BM, Kromeke O, Optenhostert-Porst C, Wolf OT (2006) Hippocampal volume and cognitive performance in long-standing Type 1 diabetic patients without macrovascular complications. Diabet Med 23: 32–39.

68. Musen G, Lyoo IK, Sparks CR, Weinger K, Hwang J, et al. (2006) Effects of type 1 diabetes on gray matter density as measured by voxel-based morphometry. Diabetes 55: 326–333.

69. Perantie DC, Wu J, Koller JM, Lim A, Warren SL, et al. (2007) Regional brain volume differences associated with hyperglycemia and severe hypoglycemia in youth with type 1 diabetes. Diabetes Care 30: 2331–2337.

70. Northam EA, Rankins D, Cameron FJ (2006) Therapy insight: the impact of type 1 diabetes on brain development and function. Nat Clin Pract Neurol 2: 78–86.

71. Siesjo BK (1988) Hypoglycemia, brain metabolism, and brain damage. Diabetes Metab Rev 4: 113–144.

72. Amiel SA (1997) Hypoglycaemia in diabetes mellitus – otecting the brain. Diabetologia 40 Suppl 2S62–68.

73. Lerch JP, Worsley K, Shaw WP, Greenstein DK, Lenroot RK, et al. (2006) Mapping anatomical correlations across cerebral cortex (MACACC) using cortical thickness from MRI. Neuroimage 31: 993–1003.

74. Narr KL, Toga AW, Szeszko P, Thompson PM, Woods RPet al (2005) Cortical thinning in cingulate and occipital cortices in first episode schizophrenia. Biol Psychiatry 58: 32–40.

75. Rauch SL, Wright CI, Martis B, Busa E, McMullin KG, et al. (2004) A magnetic resonance imaging study of cortical thickness in animal phobia. Biol Psychiatry 55: 946–952.

76. Shaw P, Lerch J, Greenstein D, Sharp W, Clasen Let al (2006) Longitudinal mapping of cortical thickness and clinical outcome in children and adolescents with attention-deficit/hyperactivity disorder. Arch Gen Psychiatry 63: 540–549.