The Brain Hierarchical Atlas: partitioning the brain into subnetworks embedding both structural and functional connections

I Diez, P Bonifazi, I Escudero, B Mateos, MA Munoz, S Stramaglia, and JM Cortes. The Brain Hierarchical Atlas. 22nd Annual Meeting of the Organization for Human Brain Mapping , OHBM 2016 [pdf]

Background. Most used brain atlases are purely anatomical or structural [1-4], and others are purely functional, like the well-known Brodmann atlas based on lesion studies or data-driven methods [5-6]. Although obtaining suitable brain partitions (or atlases) has been studied intensively (see [7] and references therein), to the best of our knowledge, the first brain partition that accounts for subnetworks that are relevant to both structure and function is the Brain Hierarchical Atlas (BHA), published in [8] and presented here.

Methods. This work was approved by the Ethics Committee at the Cruces University Hospital; all the methods were carried out in accordance to approved guidelines (Principal Investigator: Jesus M Cortes). A population of N=12 (6 males) healthy subjects, aged between 24 and 46 (33.5 ± 8.7), provided information consent forms before the magnetic resonance imaging session. For all the participants, we acquired for each subject both diffusion tensor imaging and resting functional magnetic resonance imaging. Our results have been fully validated with a population of N=30 subjects from the WU-Minn Human Connectome Project [9]. Large-scale brain networks were obtained after magnetic resonance imaging (MRI). Whilst the T1 was mainly used for preprocessing and coregistration, diffusion tensor imaging provided structural networks (SN) of gray- matter areas connected by white-matter tracts, and variations in the blood-oxygenation-level-dependent T2* signal provided functional networks (FN). Both FN and SN were obtained at a resolution of 2514 brain regions, making FN and SC to be represented by (2514*2514) matrices for each of the subjects. For details on how we built these networks, see [8]. Applying a standard hierarchical agglomerative clustering algorithm to FN and SN, we compared at the moduli level FN and SN by exploiting their hierarchical modular organization (figure 1). This approach provided a hierarchical tree or dendrogram in which nodes were progressively merged together into moduli following a nested hierarchy of “vicinity” (which reflects distance in the correlation for FN and distance in fiber number for SN). More specifically, we compared FN and SN by employing the template of hierarchical modular organization derived from FN to visualize SC, and vice versa. We searched for the best common partition shared by structure and function by maximizing the ‘cross- modularity’ index X [8], an index which is large for a given brain partition if the corresponding Newman’s modularities of the two matrices FN and SN are large and there is also a large within-module similarity between both divisions (i.e., a large fraction of existing intra-module links are shared by both networks). Thus, a large cross-modularity value indicates that both FN and SN are highly modular and, at the same time, the moduli are internally wired in a similar way.

Results. After maximization of the cross-modularity X, we found that the partition into 20 moduli derived from FN (as portrayed in figure 1) was optimal, although similar quality partitions can be obtained in the range of M from 10 to 30 [8]. The Brain Hierarchical Atlas (BHA), composed by these 20 moduli or ‘subnetworks’ maximizing X is illustrated in figure 1. We studied the similarity between the BHA and other existing brain partitions (figure 2): (i) The automated anatomical labeling (AAL) [2]; (ii) The resting state networks (RSNs) [5]; and (iii), The Brodmann atlas (included in the software http://www.mricro.com). By comparing with AAL, we found that different anatomical areas can be included in a single BHA subnetwork, clearly highlighting that a BHA subnetwork simultaneously can recruit distinct brain circuits. Conversely, the same area of the AAL might be included in several BHA subnetworks, highlighting the anatomical overlap of the latter. Similarly, we compared BHA with the well-studied RSNs; two of them displayed a strong overlap (>0.5) with two distinct BHA subnetworks, the “Sensory motor“ with the subnetwork 3 and the “Medial visual” with the subnetwork 4, while the other RSNs overlap with more than one BHA subnetworks. Finally, when comparing with the Brodmann atlas, we found that the Brodmann area number 18 was well-characterized by the subnetwork 4 and that the Brodmann area number 20 matched to the subnetwork 18.

Conclusions. The BHA (i) represents a distinct brain partition from those previously described in the literature; (ii) incorporate distinct both structural and functional brain regions into a single operative subnetwork/unit; and (iii) can overlap and share both anatomical and functional brain regions. In the same way as alterations in resting state networks have been reported in several brain pathologies and diseases, we expect that the use of the BHA, with simultaneous focus on structure and function, might help also in diagnosing disease.

References:

[1] Lancaster J, Woldorff M, Parsons L, Liotti M, Freitas C, Rainey L, Kochunov P, Nickerson D, Mikiten S, and Fox P (2000), ‘Automated talairach atlas labels for functional brain mapping’, Human Brain Mapping, vol. 10, pp. 120-131

[2] Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, Mazoyer B, and Joliot M (2002), ‘Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain’, Neuroimage, vol. 15, pp. 273-289.

[3] Eickhoff S, Stephan K, Mohlberg H, Grefkes C, Fink G, Amunts K, and Zilles K (2005), ‘A new spm toolbox for combining probabilistic cytoarchitectonic maps and functional imaging data’, Neuroimage, vol. 25, pp. 1325-1335.

[4] Desikan R, Segonne F, Fischl B, Quinn B, Dickerson B, Blacker D, Buckner R, Dale A, Maguire R, Hyman B, Albert M, and Killiany R. (2006), ‘An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest’, Neuroimage, vol. 31, pp. 968-980.

[5] Beckmann CF, DeLuca M, Devlin JT and, Smith SM (2005), ‘Investigations into resting-state connectivity using independent component analysis’. Philos. Trans. R. Soc. Lond. B. Biol. Sci., vol. 360, pp. 1001-1013

[6] Craddock R, James G, Holtzheimer P, Hu X, and Mayberg H (2012), ‘A whole brain fmri atlas generated via spatially constrained spectral clustering’, Human Brain Mapping, vol. 33, pp. 1914-1928.

[7] Craddock R, Jbabdi S, Yan C, Vogelstein J, Castellanos F, Martino AD, Kelly C, Heberlein K, Colcombe S, and Milham M (2013), ‘Imaging human connectomes at the macroscale’, Nature Methods, vol. 10, pp. 524- 539.

[8] Diez I, Bonifazi P, Escudero I, Mateos B, Munoz MA, Stramaglia S, and Cortes JM (2015), ‘A novel brain partition highlights the modular skeleton shared by structure and function’, Scientific Reports, vol. 5, pp. 10532. The atlas can be downloaded at http://www.nitrc.org/projects/biocr_hcatlas/

[9] Van Essen DC et. al. (2013), ‘The WU-Minn Human Connectome Project: An overview’, Neuroimage, vol. 80, pp. 62-79

You may also like...