Monday, March 6, 2017

Human Brain Mapping



Abstract

The origin, structure, and function of the claustrum, as well as its role in neural computation, have remained a mystery since its discovery in the 17th century. Assessing the in vivo connectivity of the claustrum may bring forth useful insights with relevance to model the overall functionality of the claustrum itself. Using structural and diffusion tensor neuroimaging in N = 100 healthy subjects, we found that the claustrum has the highest connectivity in the brain by regional volume. Network theoretical analyses revealed that (a) the claustrum is a primary contributor to global brain network architecture, and that (b) significant connectivity dependencies exist between the claustrum, frontal lobe, and cingulate regions. These results illustrate that the claustrum is ideally located within the human central nervous system (CNS) connectome to serve as the putative “gate keeper” of neural information for consciousness awareness. Our findings support and underscore prior theoretical contributions about the involvement of the claustrum in higher cognitive function and its relevance in devastating neurological disease. Hum Brain Mapp 36:827–838, 2015. © 2014 Wiley Periodicals, Inc.

INTRODUCTION

Considerable mystery has surrounded the origin, structure, and function of the claustrum—a diminutive bilateral anatomical structure of the brain whose very name means “hidden away” [Crick and Koch, 2005]. The structure was identified in humans as early as 1672 as shown by the drawings of Thomas Willis [Bayer and Altman, 1991], though was first described (under the name “vormauer”) by Karl Friedrich Burdach in his seminal work, Von Baue und Leben des Gehirns in the early 19th century. Burdach himself, however, credited the discovery to the 1,786 drawings by FĂ©lix Vicq-d'Azyr . Despite this long history, knowledge of the claustrum, its organization, afferents, and efferents has been studied intermittently and relatively poorly explored until recently [Edelstein and Denaro, 2004; Smythies et al., 2014].
In animals, extensive neuronal projections have been noted between the claustrum and numerous cortical regions [Buchanan and Johnson, 2011; Crick and Koch, 2005; Minciacchi et al., 1995; Tanne-Gariepy et al., 2002; Wilhite et al., 1986] with white matter connections to multiple cortical layers [Carey et al., 1980] (although see Crick and Koch [2005]). Likewise, the claustrum exhibits many connections to subcortical structures [Arikuni and Kubota, 1985; Berke, 1960; Buchanan and Johnson, 2011; Crick and Koch, 2005; Salerno et al., 1984; Tanne-Gariepy et al., 2002]. The structure is present in all mammals with species such as tree shrews and cats having relatively large claustra by volume in contrast to monotremes, who apparently lack claustra altogether [Butler et al., 2002].
In electrophysiological studies, the claustrum has been shown to send signals to four thalamic nuclei: the medial geniculate nucleus, the lateralis posterior, centrum medianum, and ventralis lateralis [Chachich and Powell, 2004; Cortimiglia et al., 1987, 1991; Crescimanno et al., 1989; Spector et al., 1970; see also Sherk 2013, Chapter 5]. Excitotoxic lesioning in the rat has also shown evidence of frontostriatal connectivity as well [Grasby and Talk 2013]. Both the medial ectosylvian gyrus (A2) and the anterior ectosylvian gyrus (S2) also have been reported as transmitting information to the claustrum [Edelstein and Denaro, 2004; Hassmannova, 1977]. Moreover, the claustrum is topologically organized with frontal cortex being linked to the claustrum's anterior portion, the parietal cortex with its central and posterior parts, and the occipital and temporal cortices are linked to its posterior and inferior margins [Druga, 2014; Markowitsch et al., 1984]. Nontopographic projections to other parts of the same cortical area exist, but there is complete segregation between distinct cortical areas [Minciacchi et al., 1995]. Fiber bundles such as the corona radiata, uncinate fasciculi, and inferior occipitofrontal fasciculi project to the claustrum from the superior frontal, precentral, postcentral, superior parietal, and parietooccipital regions [Fernandez-Miranda et al., 2008a, b]. Only a relatively small number of claustral cells appear to project to the contralateral hemisphere [Markowitsch et al., 1984]. Opinions have differed, however, as to whether the claustrum receive inputs via brain stem and spinal afferents [Arikuni and Kubota, 1985; Edelstein and Denaro, 2004].
In humans, the volume of the claustrum is one quarter of 1% of the volume of the cerebral cortex [Crick and Koch, 2005; Edelstein and Denaro, 2004]. As illustrated in the Talairach and Tournoux [1988] human brain atlas, the claustrum extends 22 mm inferior-to-superior and 38 mm anterior-to-posterior. The right claustrum has an average volume of 828.83 mm3, while the left has a volume of 705.82 mm3 and lies approximately 1 mm from the insular cortex and about 1 mm from the putamen [Kapakin, 2011]. Bilaterally, it displays asymmetry in shape and anisotropy between the hemispheres [Cao et al., 2003] (as shown schematically in Fig. 1A). It projects downward and crosses the rhinalis fissure—though it may be interrupted by fibers from the uncinate fasciculus in some places—before extending to the lateral rhinencephalon. A thorough review of claustral gross anatomy can be found in Druga [2014]. Historically, the claustrum was often considered to be the innermost layer of the insula [Landau, 1919], though has also been associated with the basal ganglia as a potential pallial derivative [Puelles, 2014]. The claustrum contains both fusiform cells and pyramidal somata, which are indicative of cortical areas, but also contains subcortical cell types. Due to the variation in cell morphometry, the generally nonlayered structure of the claustrum cannot formally be pronounced as strictly cortical or subcortical [Mathur et al., 2009]. This might explain its only modest description in early as well as modern neuroscience texts and resources (e.g., http://neurolex.org/wiki/Category: Claustrum).
Figure 1.
(A) In the human brain, the claustrum is a thin band of cells located medial to the insular cortex and lateral to the putamen. (B) Three-dimensional models of the bilateral claustra were obtained from T1-weighted MPRAGE structural MRI image volumes. (C) White matter fiber tractography was performed and inter-regional connectivity was computed by determining the relative proportion of extracted fibers initiated or terminated within the boundaries of each anatomical parcel. Here, fibers linking the claustrum to other brain regions are illustrated in an example subject. (D) These measures are aggregated into an M × M connectivity matrix, shown here as an image in which black refers to no or relatively low inter-regional connectivity between region i and region j, varying to white reflecting a greater extent of region. (E) Regional morphometrics for the entire brain are represented as a circular “connectogram.” The measurements for each region are shown as “heat rings” including, from the outside inward, cortical thickness, regional volume, surface area, and connectivity density. The information contained in the computed connectivity matrix is used here to illustrate the pattern and strength of connections between brain regions. Line opacity is proportional to connection density, whereas color represents the average FA integrated along all pathways comprising that connection. Red is high average FA, green in medium average FA, and blue if low mean FA, according to the upper, middle, and lowest thirds of the FA distribution. Further information on connectogram construction and interpretation can be found in Irimia, Chambers et al. [2012b]. (F) The connectogram showing only those fiber pathways which eminate from or are connected to the left and/or right claustra. These fibers were subject to the graph theoretical analysis described in the body of the text.
Recent in vivo examinations in humans using diffusion tensor imaging (DTI) have shown the claustrum to possess long-range projections to many prominent Brodmann's areas (BA) [Milardi et al., 2013], including connections to visual cortex (BA 17, 18, 19, and 39), frontal areas (BA 8, 9, 10, 11, 12, and 34), as well as superior regions (BA 7, 5, 1/2/3, 4, 6, and 8) and language areas (e.g., BA 44, 45, 31, etc). Fernandez-Miranda et al. [2008a, b, 2012; Smythies et al., 2014] have traced the extensive, intricate claustrocortical system using both microsurgery in concert with DTI mapping. For their investigation, they examine these connections in subdivisions of the claustrum. They demonstrated that the claustrocortical fibers connect the dorsal claustrum with the superior frontal, precentral, postcentral, and posterior parietal cortices with a regular topographic organization. Additionally, they reiterate that the ventral portion of the claustrum is connected with the orbitofrontal, prefrontal, and temporal cortex and with the amygdala. These results confirm the patterns of connectivity previously observed in nonhuman species (see Supporting Information, Table 1).
Table 1. Volumes (mm3) and Connectivity of Extracted Brain Structures Sorted by Connection Density per cm3
 Both hemispheresRight hemisphereLeft hemisphere
 Volume [mm3]# of connections per cm3Volume [mm3]# of connections per cm3Volume [mm3]# of connections per cm3
StructureMeanSDMeanSDMeanSDMeanSDMeanSDMeanSD
Region
Clau654.6778.8834.20210.259694.9468.8630.7238.136614.4167.0137.68110.996
G_cingul-Post-ventral620.06169.0826.49312.445667.08177.6123.5769.853573.04146.4729.40914.038
G_subcallosal780.43317.1319.5539.747801.32253.1419.8378.508759.55370.4119.26910.882
S_suborbital830.12352.3517.9028.274542.69191.3521.8799.2561117.54214.5713.9254.469
G_front_inf-Orbital848.24223.0416.80810.023913.35237.3816.68110.177783.13187.2916.9349.916
S_temporal_transverse504.74126.5516.1038.148474.36125.1217.0928.224535.12121.1415.1157.991
S_pericallosal1575.07414.0815.474.6431788.36398.2114.2543.921361.77306.8716.6864.998
S_collat_transv_post621.63209.8613.7946.271727.31224.5213.1725.638515.94124.7214.4166.818
G_cingul-Post-dorsal1408.08299.413.3923.2361389.35300.9213.5682.9911426.81298.1913.2153.471
S_orbital_lateral609.08201.9413.2757.014642.71234.4912.7387.418575.45157.1613.8126.579
G_temp_sup SUP-G_T_TRANSV967.53256.0512.035.43868.14219.4412.9655.7241066.93252.3411.0954.974
S_interm_prim_Jensen637.81318.4311.7216.453727.84333.3511.0954.204547.78276.2512.3488.077
G_temp_sup-Plan_polar1684.57416.1511.1144.9461798.09426.710.1984.1871571.05374.1712.0315.472
G_occipital_sup2998.26640.6210.9734.793262.08694.4710.634.9962734.44449.4911.3164.575
S_orbital_med-olfact1256.38232.4110.2384.1281238.48186.3810.9524.5471274.28270.539.5243.542
G_cuneus2867.94492.319.4983.2642970.91505.319.4693.1822764.97458.739.5273.359
G_and_S_transv_frontopol2064.46577.448.9174.2462436.06518.87.9983.2831692.86349.469.8364.873
G_temp_sup-Plan_tempo1646.52446.238.6853.8741504.76389.79.784.0731788.27455.847.5893.34
S_occipital_ant1230.93386.328.5613.0091275.57407.458.0092.4981186.3360.499.1133.367
S_cingul-Marginalis1745.48336.558.4252.9051900.88338.848.0292.7291590.09253.148.823.034
G_front_inf-Triangul2637.3573.88.3713.8592602.16597.028.3544.2012672.45550.358.3883.505
G_and_S_paracentral2407.92420.448.3662.2712254.69375.698.4062.1922561.15408.368.3272.357
S_subparietal1816.81445.358.2693.0941917.39479.98.3413.0851716.23384.628.1973.117
G_rectus2110.24496.878.193.3431757.5298.178.7643.462462.98394.957.6173.135
G_and_S_frontomargin1945.5402.957.7864.0241818.95386.968.4344.072072.05379.847.1393.889
G_and_S_cingul-Mid-Post2647.99445.027.311.5942788.45472.087.0391.3922507.54367.947.5811.738
G_precuneus5798.45971.417.2342.6045798.64967.057.2982.6255798.27980.637.1692.595
G_front_inf-Opercular3186.14585.497.1352.4053060.85548.067.3022.3043311.42597.486.9682.502
S_oc_sup_and_transversal1991.96489.986.8682.9422162.72488.626.5482.8591821.21430.037.1883.003
S_oc-temp_lat1539.9394.266.7992.9051589.85358.646.5042.7071489.95422.847.0933.076
G_and_S_occipital_inf2666.47591.096.7612.72533.03552.686.282.2312799.91600.737.2423.035
S_collat_transv_ant1762.77451.416.7283.6361763.61419.996.7973.8241761.93482.926.6593.457
S_parieto_occipital2846.92581.036.6362.442918.64615.716.5742.6412775.19537.76.6982.231
Pole_occipital3699.11972.866.6272.8774466.21697.55.7252.152932.01476.667.5293.221
G_oc-temp_med-Lingual4787.69758.046.3821.9984787.5711.836.1741.8024787.88805.26.5912.166
G_and_S_cingul-Mid-Ant2878.24506.946.2091.8123032.81481.315.9671.6492723.67486.476.4511.939
G_parietal_sup5605.551153.546.1752.1345014.18877.726.5572.2776196.911094.465.7931.917
S_calcarine2948.55554.526.1161.6562897.59530.375.9361.4612999.5575.826.2961.819
G_postcentral3766.18719.696.1141.93595.87687.686.3291.9233936.5713.895.8991.862
S_temporal_inf2084.3547.485.972.5972055.57470.685.6232.1252113.03615.876.3172.966
G_temp_sup-Lateral5216.72991.975.5741.9654972.02905.545.6192.0675461.421018.335.531.866
S_oc_middle_and_Lunatus1444.58425.375.5612.5981461.36468.495.522.4171427.79379.015.6032.78
G_occipital_middle4776.68919.825.5372.0435102.78953.464.9531.8694450.58759.466.1212.052
G_oc-temp_med-Parahip3684.89719.735.5051.9573882.85747.55.4811.8953486.92635.085.532.027
G_and_S_subcentral2879.21603.965.4921.5942735.02548.095.7951.7133023.4625.165.191.411
Pole_temporal5539.3870.775.2532.4515493.25829.155.4762.5515585.35912.345.032.338
S_oc-temp_med_and_Lingual3151.12568.645.0071.6762919.38488.45.11.5663382.86550.44.9141.782
G_and_S_cingul-Ant4810.7788.34.7031.1495177.32754.954.631.0934444.09638.474.7761.204
G_precentral5998.47866.494.5341.1156019.14881.324.6761.2415977.79855.344.3920.959
G_orbital6264.03884.884.4221.5776514.48873.414.4731.6646013.58827.274.371.493
G_oc-temp_lat-fusifor4553.95791.074.351.5254476.53783.594.291.4374631.38794.854.411.613
S_precentral-sup-part2055.31544.844.0862.0732103.49626.224.2362.0072007.13447.073.9362.136
S_postcentral3541.86828.883.9451.5613247.51766.414.2651.4763836.21786.513.6251.586
S_central3438.17584.793.7291.8493417.99571.373.9222.0043458.35600.13.5351.667
S_orbital2294.95411.283.6922.1162354.25430.93.6232.3112235.64383.693.7611.911
S_precentral-inf-part2440.52579.993.6891.3652523.19605.963.5011.272357.85543.313.8761.435
S_intrapariet_and_P_trans4418.54718.413.6661.4364435.42686.324.0491.3824401.66752.223.2831.392
S_front_middle2797.12771.433.5191.8563235.95700.263.0851.5922358.29562.883.9532.003
S_front_inf3342.05733.613.5071.433138.64725.793.631.4183545.47686.623.3851.439
G_pariet_inf-Angular6348.841177.873.4891.33569491050.443.2311.1935748.68978.623.7471.422
G_pariet_inf-Supramar6123.031251.573.2940.8616010.841225.923.430.8766235.231272.943.1580.829
G_temporal_middle7423.521479.973.2160.9187667.881516.673.1020.9027179.171407.783.330.925
G_front_sup16504.562260.643.0380.88515933.012206.43.1370.89717076.112178.182.940.867
G_temporal_inf6709.91263.143.0011.1956570.021279.142.9631.2116849.771237.523.0391.185
S_front_sup4343.69871.012.7641.1954119.55769.282.8311.274567.82911.982.6971.118
G_front_middle9808.581799.062.480.8729376.611737.472.4520.87510240.561763.552.5090.872
S_temporal_sup9567.531470.562.1580.56810107.041449.781.9960.4759028.021287.252.320.609
Because of its cellular composition and wide-ranging connectivity, the claustrum has been proposed as the cornerstone of sensory integration—putatively receiving, assimilating, integrating, and channeling information throughout the brain from each sensory cortex [Edelstein and Denaro, 2004; Remedios et al., 2014]. Moreover, it has been posited that the claustrum assesses the congruence of information while binding multiple facets of perception into a whole, such as assigning motion, direction, and sound to a singular visual object [Naghavi et al., 2007]. Indeed, in a thorough series of reviews, Crick and Koch [2003, 2005] and Smythies et al. [2012a, b] have argued that the claustrum might form the foundation for the neural locus of consciousness. With little doubt, these important papers are principally responsible for the recent resurgence of interest in the claustrum, its structure, connectivity, and its function [Mathur, 2014]. The mechanisms of its role in brain networks, however, is not well understood especially in light of ongoing debates concerning claustral neurons as multisensory processors [Braak and Braak, 1982; Edelstein and Denaro, 2004; Remedios et al., 2010; Smith and Alloway, 2010].
Despite its consideration in a modest number of largely descriptive neuroimaging studies, however, the relative importance or strengths of the diverse and wide-ranging pattern of claustral connectivity have not, to date, been thoroughly explored using large human samples. Quantitative methods for the assessment of white matter networks offer an opportunity to explore the connectivity of the claustrum in vivo and to measure its relative influence on brain subnetworks. Indeed, graph theoretical analytic techniques have been widely applied to connectomics data from human neuroimaging studies to characterize the novel properties of structural and functional brain networks [Bassett and Bullmore, 2009; Bullmore and Sporns, 2009]. Graph-wide as well as nodal connectivity characteristics provide insights into the integrity and efficiency of information flow through a network and can shed light on the effects of injury or disease. They can also highlight the macrolevel importance of certain regions on the overall structure of a network and how those influence the connectivity of brain subnetworks.
With this in mind, we sought to quantitatively examine the relative importance of putative subnetworks, which involve the claustrum—assessed using DTI, structural neuroimaging, and graph theoretical analysis—in a cohort large enough to establish a reliable and comprehensive model of the full macroscale connectivity of this so often over-looked brain structure.

METHODS

Subjects and Data Acquisition

The study cohort included N = 100 healthy adult subjects (47 males and 53 females) with ages between 18.6 and 61.1 (mean and SD: 32.71 ± 11.6 years). Subjects provided informed written consent as required by the Declaration of Helsinki, U.S. 45 CFR 46, and neuroimage volume acquisition was conducted with the approval of the local ethics committees at the respective research institutions where data were acquired. Participants were recruited by advertisements in local newspapers and campus flyers. Subjects were all healthy and had no history of neurological or psychiatric illnesses. No participant had a current or past psychiatric diagnosis (including substance abuse) or was taking medications for any medical reasons. Additional exclusion criteria for all participants included left-handedness, hypertension, neurological illness, metal implants, and a history of head trauma with loss of consciousness for more than 5 min. Neuroimaging datasets were fully anonymized, and no linked coding or keys to subject identity were maintained. For these reasons, in compliance with the U.S. Health Insurance Portability and Accountability Act (HIPAA; http://www.hhs.gov/ocr/privacy), this study does not qualify as involving human subjects' materials.
Structural T1-weighted magnetic resonance imaging (MRI) and DTI volumes were acquired from each patient using a Siemens Trio Tim 3.0 Tesla whole body scanning system with a 12-channel head coil. For MRI, a Turbo MP-RAGE sequence (repetition time [TR] = 20 ms, echo time [TE] = 3 ms, flip angle = 25°, slice thickness = 1 mm, acquisition matrix = 256 × 256 × 256) was obtained. A DTI protocol having the following parameters was obtained: TR = 9.4 s, TE = 88 ms, flip angle = 90°, slice thickness = 2 mm, number of gradient directions = 68, acquisition matrix = 128 × 128 × 128. Two nondiffusion-weighted volumes were acquired for each subject (b-values: 0 and 1,000).

Neuroimage Data Processing

Data processing workflows were created using the Laboratory of Neuro Imaging (LONI) Pipeline Workflow Environment (http://pipeline.loni.ucla.edu; version 5.9.1). The LONI Pipeline is a graphical environment for construction, validation and execution of advanced neuroimaging data analysis protocols [Dinov et al., 2010]. It enables the creation of heterogeneous tool interactivity, automated data format conversion, allows grid utilization, facilitates data provenance, and provides a significant library of computational tools. It is built as a distributed grid computing environment and permits efficient tool integration, protocol validation, and broad resource distribution. Further details are available via the LONI Pipeline website.
A set of customized LONI Pipeline workflow was created to perform basic preprocessing of all subject data as a prelude to the application of regionally focused processing, atlas creation, and connectomics estimation. This initial workflow performed FreeSurfer 5.3 reconstruction [Dale et al., 1999; Fischl et al., 1999], stripped the skulls and other nonbrain tissues out of the image volumes, bias-corrected the T1-weighted structural images, performed eddy current correction on the DTI, and then reconstructed the diffusion image into 3D tracts in both .vtk and .trk file formats so that the white matter can be visually examined using commonly available tractography software. Specifically, DTI fractional anisotropy (FA) computation and fiber reconstruction was performed using the TrackVis Diffusion Toolkit (http://trackvis.org) via streamline tractography using a subject-specific B0 white-matter mask and a turning angle threshold of 35°. The workflow also carried out the spatial registration of the T1 and DTI images, and likewise applied the transformation to the cortical parcels of the FreeSurfer labelmap. Finally, connectivity matrices were computed to quantify the relative number of extracted DTI fibers and average FA values over them which connected each pair of cortical regions. These matrices are then combined with the cortical region morphometric analyses to generate whole-brain connectograms (see below). The LONI Pipeline workflows described herein, along with additional supporting workflows, are available as Supporting Information.

Identification and Segmentation of the Claustrum

While the claustrum has been assigned labelmap values in the context of the FreeSurfer implementation of the Destrieux et al. [2010] brain atlas, the version of FreeSurfer we used (version 5.3) does not, in fact, assign voxels to this structure. So, following the application of the initial image processing workflow, 3D Slicer (http://slicer.org) was used to manually delineate left and right claustrum masks in the sample of healthy subjects. For each subject, an axially resampled T1 anatomical was examined in a superior-to-inferior manner to identify the body of the claustrum. As described above and illustrated in Figure 1A, the claustrum is located between the external and internal capsules, medially from the insular cortex and laterally from the body of the putamen. Separate left and right binary masks of this region were manually drawn (by Carinna M. Torgerson (CMT)) in the axial plane and then edited in the sagittal plane for accuracy, before being visually reviewed by trained raters (by Andrei Irimia (AI) and John D. Van Horn (JDVH)), and the regional delineation refined as necessary.

Pooling and Applying Claustral Masks

Each subject's pair of claustral masks were then used as input into an “Atlas Creation” LONI Pipeline workflow, which spatially warped the claustrum label masks into Montreal Neurological Institute (MNI) Atlas space and, using FSL modules, spatially averaged them, before returning the average back into the native space of each individual subject via inverting the spatial warping transform. This averaging was performed to make claustral delineation consistent with the FreeSurfer mapping technique, which uses averaged data to generate accurate labelmaps. These newly applied masks of the claustrum (now in native space for all 100 subjects), the “Adding Labelmaps” workflow was used to spatially merge the left and right claustrum labels into the FreeSurfer parcellation and segmentation labelmap. At this stage, new connectivity matrices were generated for each subject based upon whole-brain labelmaps, which now included the spatially averaged and unwarped claustra, using the “Connectivity Calculation Workflow.”
Subsequently, the fully integrated labelmaps, which include the left and right claustra, were used as input into three different workflows. First, the “Whole Brain Plus Claustrum” workflow took the full labelmaps and generated connectograms depicting the whole, macroscale pattern of DTI connectivity across the complete collection of segmented brain regions. Finally, the “Only Claustrum” workflow used the connectivity matrices from the whole-brain plus claustrum analysis resulting in new files created for each subject depicting the specific connectivity patterns of the left and right claustrum in each subject. 1 Values of all connections originating or terminating in the insular cortex an putamen were set to zeroes to avoid misattribution due to proximity, as well as double-counting errors.

Connectogram Generation

Cortical parcellations were represented as a circular arrangement of 165 positioned elements representing the connectivity matrix of the left (quadrant IV of Fig. 1D) and right (quadrants I of Fig. 1D) cerebral hemispheres, and their jointly connecting fiber pathways (quadrants II and III), each positioned symmetrically with respect to the vertical axis. This graphical representation—termed a “connectogram”—has been described extensively elsewhere [Irimia et al., 2012a, b; Van Horn et al., 2012]. The connectogram is becoming widely used as a method for showing patterns of connectivity in a variety of research domains, as is demonstrated by the variety of topics mentioned on its Wikipedia page (http://en.wikipedia.org/wiki/Connectogram). Briefly, the connectograms represent cortical region labeling and morphometric measurements as color-coded rings and illustrate inter-regional connectivity as lines between lobar and regional segments [Irimia et al., 2012a, b]. Line opacity is proportional to connection density, while in these figures, specifically, red indicates high average FA (aggregated over fibers connecting those regions), green represents moderate FA, and blue represents low FA. FA measures the degree of preferred water diffusivity, ranging between zero (no preferred diffusion) and unity (strong directional diffusion) [Basser et al., 2000]. A complete list of parcellations, their abbreviations and associated Red-Green-Blue (RGB) codes are provided in the supplemental materials of Van Horn et al. [2012]. Parcellations were arranged within each lobe in the order of their location along the anteroposterior axis of the cortical surface associated with the published FreeSurfer normal population atlas [Destrieux et al., 2010].

Network Theoretical Analyses

To examine the network contributions of the claustrum, we used graph theoretical methods, which have been demonstrated to provide useful constructs in characterizing neuroimaging-derived brain networks [Sporns, 2012]. Graph theoretical metrics are useful for describing the properties of network architecture and have been successfully applied to patterns of inter-regional brain connectivity. Additionally, using a method for systematic removal of white matter connections in conjunction with graph theoretical analysis [Irimia and Van Horn, 2014], the left and right claustra were analyzed for their role in major brain networks (Fig. 1D). Such methods have provided a convenient approach for the computational modeling of the effects of brain injury [Van Horn et al., 2012] and for the assessment of nodal clustering and assortativity which can illustrate the presence of various brain sub-networks [van den Heuvel et al., 2008].
All connectivity metrics were determined using the Brain Connectivity Toolbox (https://sites.google.com/site/bctnet/). Specifically, the graph-wide metrics used here included network density, diameter, and assortativity, as well as the nodal metric betweeness centrality. To simplify interpretation, particularly short (<5 a="" accompanied="" analyses.="" analysis="" and="" as="" being="" by="" claustra="" claustrum="" close="" connections="" considered="" determine="" disregarded="" errors="" excision="" graph-wide="" i="" in="" included="" increase="" influential.="" insula="" leave-one-out="" likelihood="" manova="" metrics="" mm="" modeling="" most="" neighborhood="" not="" notably="" of="" our="" p="" performed.="" proximity="" putamen="" s="" statistical="" statistically="" strategy="" subnetwork="" the="" their="" then="" this="" to="" topological="" treatment="" type="" was="" were="" which="" with="" would="">

RESULTS

In our large sample of N = 100 healthy subjects, the claustrum was found to possess the highest density of fiber connections per unit volume out of all examined brain regions (Table 1; listed by name according to the FreeSurfer aparc.a2009 naming scheme, which can be found at https://surfer.nmr.mgh.harvard.edu/fswiki/FsTutorial/AnatomicalROI/FreeSurferColorLUT following the list of aparc.a2005 labels). The normalization by regional volume is commonly recommended practiced when comparing fiber tract densities [Bassett et al., 2011; Cahalane et al., 2012; Varkuti et al., 2011]. The complete pattern of inter-regional connectivity between all regions is illustrated as a connectogram in Figure 1E, while the specific pattern of claustral connectivity is depicted in Figure 1F. The claustra had their densest DTI fiber tract connections to frontal cortices, with more modest degrees of connectivity to parietal, temporal, and occipital lobes, respectively. Connections appeared the fewest involving limbic structures in agreement with the findings of Markowitsch et al. [1984]. Evidence was also found for fiber tracts between the claustra and the brain stem, in contrast to prior reports which hypothesized that the claustrum communicates with the brain stem and autonomic spinal neurons via intermediaries, such as the medial prefrontal cortex [Edelstein and Denaro, 2004; Hatam et al., 2013].
Our statistical analyses of the claustrum's topological neighborhood subnetwork, as measured using graph theoretical metrics with claustrum excision being considered as a level of treatment, were as follows.

Analysis of Graph-Level Metrics

The three specific network feature variables included in a MANOVA statistical model were assortativity, density, and graph diameter. We note that of the variety of available graph-based metrics, these were found to be only ones which were negligibly intercorrelated. But, importantly, they are particularly useful for assessing overall network structure and sensitivity to network alteration [Sporns, 2011]. These are defined as follows: (a) density is the ratio of present connections to all possible connections; (b) assortativity measures the correlation coefficient between the degrees of all nodes on two opposite ends of a link. A positive assortativity coefficient indicates that nodes tend to link to other nodes with the same or similar degree; while (c) graph diameter refers to the maximum network eccentricity, where eccentricity is the maximal shortest path length between a node and any other node. Under the MANOVA, it was found that the removal of the claustrum led to a statistically significant change in the network measure feature vector at the omnibus level, F(3,196) = 3.49, P ≤ 0.017 (FDR corrected). Systematic removal of individual network measures in the MANOVA was performed in order to identify their relative contribution to the overall model. The results, indicated by which variable had been removed, are as follows: assortativity—F(2,197) = 2.62; P ≤ 0.07; diameter—F(2,197) = 4.73, P ≤ 0.001; and density—F(2,197) = 1.87, P ≤ 0.16.

Analysis of Node-Level Metrics:

Betweenness centrality assesses the fraction of all shortest paths in the network which contain a given node. Of particular interest at the nodal level is the notion of the betweenness centrality of brain regions connecting to the claustrum and how betweenness centrality of nodes is altered when the connectivity to the claustrum is removed. Nodes with high values of betweenness centrality participate in a large number of shortest paths. Assessing this measure of a node's relative “position” in a network is often a more informative measure than measuring connectivity density alone [Sporns, 2011]. Measurement of betweenness centrality has become a useful strategy for understanding the elements of complex networks, notably in the brain [Mirzasoleiman and Jalili, 2011]. It was found that the systematic removal of claustral connections led to significant changes (P ≤ 0.05, FDR corrected) in the ANOVA main-effects models of nodal betweenness centrality involving the structures listed in Table 2.
Table 2. Influence of the Claustrum Removal on Nodal Betweenness Centrality (P < 0.05, FDR)
  With claustrumClaustrum removed   
StructureHemisphereMeanSDMeanSDFStudent's t-testP
G_front_inf-OpercularRight59.9817.4264.8915.574.41682.10160.0369
G_precentralRight37.9512.9444.179.383.89841.97440.0497
G_postcentralRight42.9815.8247.9015.235.00532.23730.0264
G_front_middleRight47.4711.9151.4914.954.42432.10340.0367
G_front_middleLeft45.9210.8949.4111.784.73872.17690.0312
G_postcentralLeft40.5214.2144.5713.344.29942.07350.0394
G_precentralLeft39.9211.2742.798.744.04882.01220.0456
G_front_inf-OpercularLeft55.6214.7665.2517.654.18532.04580.0421

DISCUSSION

As examined using tract-tracing and other methods [Berke, 1960; Markowitsch et al., 1984], claustral neuroanatomy is noteworthy across multiple species for receiving input from almost all regions of cortex and then directly projecting back to them [Buchanan and Johnson, 2011; Edelstein and Denaro, 2004; Kowianski et al., 1999; Smythies et al., 2012a, b]. Yet, while prior descriptive or small-sample human neuroimaging examinations of the claustrum's white matter architecture have been performed, to our knowledge, ours is the first population-level study to quantitatively examine the structural connectivity of the claustrum using in vivo brain mapping and graph analytic methods.
Here, using DTI tractography and graph theoretical analytics, we computationally verify in a large human sample (N = 100) that the human claustrum is widely connected as assessed by its contribution to graph-wide diameter. Dense connectivity exists to most cortical areas although the structure is only modestly connected to limbic and occipital regions. While many contemporary examinations of the claustrum subdivide the claustrum [Fernandez-Miranda et al., 2008a, b, 2012; Gattass et al., 2014; Smith and Alloway, 2014; Smith et al., 2012], we did not feel that any such divisions would be necessary for our study, as boundaries are generally drawn according to perceived function or previously reported cortical region connectivity. In contrast, this study sought to solidify understanding of the entire connectome of the claustrum and attempt to better comprehend how this connectivity shapes networks; theoretically, one could use our network analysis in the future to subdivide the region according to network participation. Careful and accurate delineation of the claustrum itself, however, is already difficult at standard imaging resolution, and therefore, research that seeks to understand single networks that the claustrum participates in ought to be carefully designed to achieve fine resolution.
In our results, the claustrum appears to play a central role in linking multiple disparate structural brain networks. Conversely, the nodal betweenness centrality of predominantly frontal regions connected to the claustrum change significantly more if their connections to the claustrum are systematically removed than any other region. Traditionally, betweenness centrality increases as the relative contributions of shorter pathways come to dominate a network. Increases in betweenness centrality following removal of claustrum connectivity indicate that the white matter networks associated with them are, more heavily influenced by shorter—presumably local—pathways. This illustrates that the claustrum occupies a unique, and presumably critical, “location” in the overall architecture of network connectivity in the brain. Given the traffic converging on this hub, damage to the claustrum may be difficult to compensate for, since efficiency would be so greatly compromised by distributing such a large number of connections over a larger area.

Indeed, disruption to corticoclaustral connectivity—forcing a greater dependence among local regional connections—appears to result in a variety of neurological symptoms (Supporting Information, Table 3). Studies from a variety of human neurological and psychiatric syndromes [Smythies et al., 2014] have identified altered claustral morphometry or have found disrupted patterns of white matter connectivity.

Disorders such as Wilson's Disease, for example, are noteworthy for their specificity of claustral involvement. Wilson's Disease is an autosomal recessive genetic disorder in which copper accumulates in brain tissues [Lorincz, 2010]. Copper deposition in the claustrum is notable using T1-weighted neuroimaging and is often considered a neurological hallmark of the disease, and can lead to general executive control difficulty, as well as various symptoms of memory dysfunction [Sener, 1993]. Other neurological conditions have been examined in the claustrum as well.

Negative correlations between anhedonia and metabolism in the claustrum have been shown in both patients with unipolar depression and bipolar disorder [Chen et al., 2011]. Claustral amyloid plaques accumulation has been implicated in the outcomes of Alzheimer's disease and aging [Fernandez-Miranda et al., 2008a, b; Morys et al., 1994]. Severity of delusions in schizophrenia is correlated with the reduction in left claustral volume [Cascella et al., 2011] and schizophrenia patients with hallucinations show signs of white matter excesses [Shapleske et al., 2002].

See Supporting Information for further summary of clinical syndromes associated with claustrum dysfunction.
Our identification of the strongest network dependencies existing between the claustrum, frontal, and cingulate cortices is a particularly noteworthy result. The simulated excision of the claustrum's representation in the overall network forces these regions to effectively bias their connectivity via more spatially local pathways. Such changes do not appear in other brain areas.

This finding supports previous findings of frontal lobe connectivity across many methodologies, including electrophysiology [Berke, 1960; Chachich and Powell, 2004; Cortimiglia et al., 1991; Grasby and Talk, 2013], anterograde, and retrograde tracers [Arikuni and Kubota, 1985; Bayer and Altman, 1991; Pearson et al., 1982; Sadowski et al., 1997], and neuroimaging [Fauvel et al., 2014; Milardi et al., 2013; Park et al., 2012] and highlights the importance of widening the scope of functional observations of the claustrum beyond the primary sensory and motor regions.

Some authors have suggested that the claustrum may act to synchronize cortical subnetworks, which are responsible for a variety of coordinated behaviors [Smith and Alloway, 2010] possibly serving to counterbalance spatially over-represented cortical regions [LeVay and Sherk, 1981; Minciacchi et al., 1995]. The claustrum may also be particularly useful in sensory cross-modal matching [Arnow et al., 2002] and possibly modulate the cortical neuron receptive field properties [Shima et al., 1996].

 Such putative roles emphasize the importance of inhibitory claustral connectivity in the regulation of conscious behaviors. Furthermore, it lends credence to the suggestions of Smythies et al. [2012a, b] that, unlike the thalamus, the claustrum may not be a strict multisensory processor, per se, but may be indispensable in organizing the information used by multisensory processors as relevant to the brain's executive functions.

In this way, the claustrum may well serve to filter signals about the relationships of all the thousands of sensorimotor inputs from the outside world to and from frontal and cingulate processing subnetworks.
Our study is, however, not without several limitations. Due to the paucity of large sample, in vivo human neuroimaging studies of the claustrum, we chose to focus more broadly on the macroscale structural connectomics properties of the structure, rather than examine factors which influence between-group differences.

For instance, we pooled male and female subjects together to get a picture of the connectivity of the average claustrum. We did not consider handedness or age-related effects. Handedness, gender, age, ethnicity, and a number of other phenotypic as well as genetic factors, will be very important for understanding the role of the claustrum in developmental, mental health, and age-related disorders.


Additionally, we made no attempt to identify subdivisions of the claustrum. Already exploring one of the smaller structures of the brain, it is arguable whether conventional imaging sequences provide sufficient spatial resolution or image contrast to permit the identification of tissue subclasses or reliable spatial landmarks thereof.

While earlier studies have suggested a compelling approach to parcellating claustral subdivisions using Magnetic Resonance (MR) [Pathak and Fernandez-Miranda, 2014], further validation and testing of such approaches are likely required before they can be of practical utility. Thus, here, we chose to focus on the claustrum as a single entity, which likely limits our ability to assign claustro-cortico projections with high specificity. Advances in ultrahigh field, high-resolution MRI will make such studies possible.

All in all, the high-resolution imaging of the claustrum and the exploration of phenomic variables is encouraged in future studies to determine how the structure and connectivity of the claustrum varies between genders, how it develops, and is altered due to aging and age-related disease.

While we have discussed how our results may support some of the functional findings of previous investigations, further research using functional imaging—such as fMRI—would be necessary to confirm that these functions are in fact correlated with these putative connections. Moreover, traditional functional imaging methods may not provide sufficient spatial resolution to accurately assess such a small structure; technological advancements may be necessary to obtain an accurate picture of the behavior of the claustrum. Hagmann et al.

 [2008] applied community detection or modularity analysis to demonstrate a close relationship between structural connections and functional connections. Therefore, an analysis of the structural connections of the claustrum may help generate hypotheses regarding its functional role in the networks to which it contributes. Despite these caveats, our results are in line with previous neuroimaging and non-neuroimaging findings.

Given the results of our large-scale analysis—and in concurrence with previous insightful research [Crick and Koch, 2005]—we concur with the hypothesis that the claustrum is likely more than a simple neural relay station in the overall context of inter-regional white matter networks. Its differential influence on brain region subnetworks typically associate with high cognitive activity, attention, and action suggest that the claustrum plays a central role in linking multiple sensory networks with those regions, which can interpret and take action on such information. In support of previous hypotheses on the role of consciousness [Smythies et al., 2012a, b], our results show that the claustrum's position in the macroscale network would be structurally capable of allowing salient sensory information to enter brain regions associated with conscious thought while keeping superfluous input out of awareness.

In conclusion, the claustrum embodies a highly connected yet still curious brain structure whose in vivo examination in humans using neuroimaging has remained challenging due to its small, thin, irregular shape, and location. Determining the functions of the claustrum in humans has formed the basis of recent empirical and theoretical exploration [Smythies et al., 2014]. Our study highlights the claustrum as a critical and “ideally located” component in brain connectomic architecture upon which frontal and cingulate regions appear particularly dependent.
These results have particular relevance for neurological disorders and diseases that affect the structure and connectivity of the claustrum and provide further support for the claustrum as a putative neural filter for conscious thought.

ACKNOWLEDGMENTS

The authors wish to acknowledge the artistic talents of the late Mr. Carlos Mena, his creativity, friendship, and, for this article, his contribution to the graphical illustration of claustral anatomy. Finally, the authors thank the dedicated staff of the Institute for Neuroimaging and Informatics based at the University of Southern California, Los Angeles, CA.

CONFLICT OF INTEREST STATEMENT

This research was conducted in absence of any commercial or financial relationships which could be construed as a potential conflict of interest.

CONTRIBUTIONS

CMT pre-processed all neuroimaging data, helped to prepare the manuscript and supplemental materials. AI conducted all graph theoretical connectivity analyses. SYMG helped to prepare figures and review the manuscript text for technical correctness. JDVH conceived of the examination, guided the necessary data processing and analyses, and helped to author the manuscript and supplemental materials. All authors had the opportunity to discuss the results, participated in writing/editing, and made comments on earlier versions of the manuscript.

Footnotes

  1. 1 Note: In our examination, due to issues of image resolution and the uncertainty in the literature over the precise intraregional boundaries for so doing, we purposefully did not attempt to segment the claustrum into putative substructures. However, see the chapter by Pathak and Fernandez-Miranda (2014), for a potential approach for subdividing the claustrum based on MRI data.

No comments:

Post a Comment