ORIGINAL RESEARCH article

Front. Psychiatry, 09 January 2014

Sec. Neuroimaging

Volume 4 - 2013 | https://doi.org/10.3389/fpsyt.2013.00186

Asymmetry of the Structural Brain Connectome in Healthy Older Adults

  • 1. Department of Neurosciences, Medical University of South Carolina, Charleston, SC, USA

  • 2. Department of Psychology, University of South Carolina, Columbia, SC, USA

  • 3. Department of Communication Sciences and Disorders, University of South Carolina, Columbia, SC, USA

Abstract

Background: It is now possible to map neural connections in vivo across the whole brain (i.e., the brain connectome). This is a promising development in neuroscience since many health and disease processes are believed to arise from the architecture of neural networks.

Objective: To describe the normal range of hemispheric asymmetry in structural connectivity in healthy older adults.

Materials and Methods: We obtained high-resolution structural magnetic resonance images (MRI) from 17 healthy older adults. For each subject, the brain connectome was reconstructed by parcelating the probabilistic map of gray matter into anatomically defined regions of interested (ROIs). White matter fiber tractography was reconstructed from diffusion tensor imaging and streamlines connecting gray matter ROIs were computed. Asymmetry indices were calculated regarding ROI connectivity (representing the sum of connectivity weight of each cortical ROI) and for regional white matter links. All asymmetry measures were compared to a normal distribution with mean = 0 through one-sample t-tests.

Results: Leftward cortical ROI asymmetry was observed in medial temporal, dorsolateral frontal, and occipital regions. Rightward cortical ROI asymmetry was observed in middle temporal and orbito-frontal regions. Link-wise asymmetry revealed stronger connections in the left hemisphere between the medial temporal, anterior, and posterior peri-Sylvian and occipito-temporal regions. Rightward link asymmetry was observed in lateral temporal, parietal, and dorsolateral frontal connections.

Conclusion: We postulate that asymmetry of specific connections may be related to functional hemispheric organization. This study may provide reference for future studies evaluating the architecture of the connectome in health and disease processes in older individuals.

Introduction

With recent advances in structural neuroimaging, it is now possible to track medium and large-scale pathways of white matter fibers and to construct the map of neural connectivity across the entire brain (the brain connectome) (, ). The brain connectome constitutes a promising development in neuropsychiatry since many physiological and pathological processes are believed to affect the architecture of neural networks (). For example, normal cognitive development spanning from childhood to senior years is traditionally believed to be associated with maturation of neural networks supporting cognitive domains such as attention (), memory (), language (), and executive function (). Similarly, neurological diseases including epilepsy and dementia are associated with pathological rearrangements or impoverishment of normal networks at a systems level (, ). Furthermore, psychiatric diseases such as schizophrenia (), bipolar disorder (), and addiction () are related to reinforcements of pathological networks.

In order to better understand how different biological processes affect the normal connectome, it is important to accurately characterize the connectome organization in healthy individuals. Specifically, it is important to define the degree of individual variability and hemispheric asymmetry that can be expected from the normal population.

Since hemispheric asymmetry is abnormal in neuropsychiatric diseases such as schizophrenia (), bipolar disorder (, ), and depression (), in this study we aimed to describe the patterns of hemispheric asymmetry and individual variability of the brain connectome obtained from a cohort of healthy senior individuals. We employed high-resolution magnetic resonance images (MRI) to reconstruct structural brain connectivity based upon white matter pathways linking anatomically defined cortical regions of interest (ROIs). We aimed to describe the group-wise distribution of asymmetries involving cortical connectivity and white matter pathways.

Materials and Methods

Subjects

We studied 17 right-handed healthy subjects (mean age 53 years, SD = 7 years, range = 40–76) who were recruited from the local community. None of the subjects had a history of neurological, psychiatric, or chronic medical illnesses. All patients signed an informed consent to participate in this study. The Institutional Review Board at the University of South Carolina approved this study.

Image acquisition

All subjects underwent MRI scanning at a 3T Siemens Trio equipped with a 12-channel head coil located at the University of South Carolina, yielding: (1) T1-weighed images (3D MP-RAGE, TR = 2250 ms, TE = 4.15 ms, 256 × 256 matrix, 256 mm × 256 mm FOV, parallel imaging GRAPPA = 2, 80 reference lines, TA = 377 s); and (2) Diffusion weighted images (dMRI)-EPI scan (30-directions with b = 1000 s/mm2 and b = 2000 s/mm2, TR = 6100 ms, TE = 101 ms, 82 × 82 matrix, 222 mm × 222 mm FOV, parallel imaging GRAPPA = 2, 80, 45 contiguous 2.7 mm axial slices, TA = 390 s).

Image pre-processing

The construction of the connectome involved two parallels pre-processing steps, namely, the segmentation of the cerebral cortex into multiple anatomical ROIs and reconstruction of white matter fibers. These steps are explained below:

Segmentation of the cerebral cortex

T1-weighted MR images were converted into NIfTI format utilizing the dcm2nii tool from the MRIcron software package (). Images in native space were non-linearly normalized into standard MNI space using the Clinical Toolbox () employing unified segmentation-normalization routines as part of the software Statistical Parametric Mapping (SPM8). Of note, the Clinical Toolbox was particularly designed to accurately quantify tissue volumes in seniors and older adults ().

This step yielded probabilistic maps of gray and white matter in MNI space.

Next, the base b = 0 T2-weighted dMRI volume was linearly transformed to standard space utilizing a boundary-based registration approach (). This step was performed using FMRIB’s Linear Image Registration Tool (FLIRT), as part of FMRIB Software Library (FSL) (). The registration parameters were then used to transform the Automated Anatomical Labeling (AAL) atlas () and the white and gray matter probabilistic maps onto the dMRI space. Once in dMRI space, a map of cortical regions segmented according to AAL was obtained by overlaying the registered AAL atlas onto the registered probabilistic gray matter map. The intersection between these images (including only voxels with a probability greater than 50% of being gray matter) represented the segmented cortical map. A list describing all ROIs used in this study can be observed in Table A1 in Appendix.

White matter fiber reconstruction

Extraction of diffusion gradients was performed with dcm2nii (). The dMRI volumes were aligned to the b = 0 dMRI image using the FSL FLIRT tool (). In diffusion space, whole brain tractography was reconstructed with the software Diffusion Toolkit () according to the following parameters: (1) angle threshold = 45°, (2) inclusion mask derived from the average of diffusion weighted signal and from the white matter probabilistic map registered to dMRI space, (3) FACT propagation algorithm, (4) spline filter.

Calculation of the connectome

From each patient, the path of each tractography streamline was assessed. All streamlines were seeded in white matter, and the end-points of each streamline were computed. Streamlines with end-points within ROIs were counted as links between these ROIs. Streamlines with end-points outside ROIs were discarded. After all streamlines were assessed, the result was a weighted connectivity matrix A, where the entry Aij represented the number of streamlines connecting regions i and j (i.e., the weighted link between i and j). Note that only direct links between regions i and j were included in the link Aij.

Finally, each link weight was corrected based on the surface of the connected ROIs and the distance between the ROIs, as proposed by Hagmann et al., where the link weight is inversely proportional to the sum of the linked ROI surfaces and fiber length, in order to account for tractography bias related to size of connected ROIs and distance traveled by the streamline ().

Connectome asymmetry

Asymmetry was evaluated for ROIs and for links

Regions of interested asymmetry was calculated by assessing the sum of link weights connecting an ROI, in comparison with the homologous ROI in the contralateral hemisphere. The connectivity of each ROI was computed without discrimination regarding the opposite end of the streamline. For example, when computing the connectivity of the left hippocampus, all possible connections of the left hippocampus were computed, including connections to ROIs in the ipsilateral and contralateral hemispheres. Connections of the left hippocampus to itself or to the contralateral hippocampus were discarded.

For each ROI, cortical connectivity asymmetry was calculated according to the following asymmetry index (AI) = (R-L)/[(R + L)/2]; where L represents the connectivity of the ROI in the left hemisphere and R represents the connectivity of the ROI in the homologous ROI in the right hemisphere. A one-sample t-test was performed to evaluate whether the distribution of ROI asymmetries across all subjects was statistically different than a distribution with mean = 0. This step was performed for each ROI. Note that ROI asymmetry represents the sum of connections to an ROI, irrespective to which other ROI is being linked to that ROI. Thus, ROI asymmetry should be interpreted in the context of link-wise asymmetry.

Link-wise hemispheric asymmetry was calculated by assessing the difference in weight for each link, according to the same asymmetry index AI = (R-L)/[(R + L)/2]; where L represents the link between ROIs within the left hemisphere and R represents the link between the homologous ROIs within the right hemisphere. Note that inter-hemispheric connections were excluded from this calculation. Cerebellar links were also excluded. Only supratentorial links within the same hemisphere were assessed. A one-sample t-test was performed to evaluate whether the distribution of asymmetries for each link across subjects was statistically different than a distribution with mean = 0.

Results

The average connectivity matrix is demonstrated in Figure 1. The ROIs are numbered in accordance with the glossary from Table A1 in Appendix. Briefly, regions 1–45 represent ROIs located in the left hemisphere, whilst regions 46–90 represent ROIs located in the right hemisphere. As such, within the average connectivity matrix with 90 × 90 entries, the upper left quadrant demonstrates links between ROIs in the left hemisphere, and the lower right quadrant demonstrates links between ROIs within the right hemisphere. The right upper and the left lower quadrants of the connectivity matrix represent reciprocal left to right connections. Since connections are not directed (i.e., the strength of connectivity from region i to j is the same as the connectivity from region j to i), the matrix is symmetrical along its main diagonal. As expected, links within the same hemisphere exhibited a higher weight compared with connectivity links between different hemispheres.

Figure 1

The average link-wise hemispheric asymmetry is also demonstrated in Figure 1. This matrix only includes 45 × 45 entries, where each cell represents the asymmetry index for each link. In order to avoid false estimations of asymmetry in links that were tracked only in a few subjects, we only included links that were tracked in greater than 75% of subjects (i.e., at least 13/17 subjects). Each matrix cell represents the asymmetry index for the connection between the ROI listed in the column and the ipsilateral ROI listed in the row. The ROIs are numbered in accordance with Table A1 in Appendix (from 1 to 45, regardless of side).

The distribution of ROI asymmetries is shown in Figure 2. The box-plot demonstrates the range of asymmetries across all subjects for each ROI. As expected, there were no regions of extreme asymmetry. A one-sample t-test revealed that among all possible 45 ROIs, a significant leftward asymmetry was noted on the inferior occipital region, fusiform gyrus, lingual gyrus, parahippocampal gyrus, amygdala, inferior frontal operculum, superior frontal gyrus, and mid cingulate gyrus. Conversely, a rightward asymmetry was noted for the middle temporal gyrus and superior frontal orbital region. These results are shown in Table 1.

Table 1

RegionAsymmetry indexpT
LEFTWARD ASYMMETRY
Occipital inf−0.03810.0023−3.6133
Fusiform−0.04860.0065−3.1304
Lingual−0.03760.0067−3.1154
Parahippocampal−0.03200.0073−3.0732
Amygdala−0.02370.0241−2.4916
Frontal inf oper−0.02570.0361−2.2873
Frontal sup−0.02800.0464−2.1587
Cingulum mid−0.02320.0488−2.1321
RIGHTWARD ASYMMETRY
Temporal mid0.03680.01002.9201
Frontal sup orb0.04450.02832.4112

Asymmetry of cortical ROI global connectivity.

Only ROIs with a t-test p value less than 0.05 are shown.

Figure 2

Link-wise asymmetry was observed toward both hemispheres. A significant leftward asymmetry was noted on the following reciprocal connections: amydgala to parahippocampal gyrus; middle to inferior occipital regions; fusiform to lingual gyri; fusiform to occipital inferior gyrus; insula to inferior frontal opercular region; fusiform to parahippocampal gyrus; precuneus to lingual gyrus; angular to superior parietal region; precuneus to mid cinculate gyrus. Conversely, a significant rightward asymmetry was noted on the following reciprocal connections: superior temporal gyrus to rolandic opercular region; middle temporal gyrus to middle occipital region; inferior to superior parietal regions; insula to inferior frontal triangularis region; precuneus to superior parietal region. These results are summarized in Table 2.

Table 2

Link betweenAsymmetry indexPT
LEFTWARD ASYMMETRY
AmygdalaParahippocampal−0.84110.0011−3.9764
Occipital infOccipital mid−0.81390.0080−3.0298
FusiformLingual−0.54960.0087−2.9857
FusiformOccipital inf−0.72230.0113−2.8617
InsulaFrontal inf oper−0.77260.0156−2.7069
FusiformParahippocampal−0.46850.0239−2.4952
PrecuneusLingual−0.66170.0406−2.2281
AngularParietal sup−0.63780.0421−2.2092
PrecuneusCingulum mid−0.53400.0453−2.1715
RIGHTWARD ASYMMETRY
Temporal supRolandic oper0.64820.00403.3592
Temporal midOccipital mid0.64240.00982.9310
Parietal infParietal sup0.69730.01132.8605
InsulaFrontal inf tri0.49120.01292.7959
PrecuneusParietal sup0.68940.01572.7016

Asymmetry of cortical links.

Only links with a t-test p value less than 0.05 are shown.

These results are demonstrated in Figure 1, which shows the distribution of T scores and p values for all links, represented in 45 × 45 matrices, where each matrix cell demonstrates the T value or the p value for the asymmetry distribution for the connection between the row ROI and the column ROI. Figure 3 demonstrates the anatomical distribution of links with an absolute mean asymmetry higher than 0.5. The three-dimensional anatomical reconstruction of regional links was defined based on an in-house developed atlas of anatomical connectivity involving all 90 ROIs used in this study. The location of travel of all streamlines connecting each possible pair of ROIs was defined based on the spatial distribution of center points (or “centroid”) of serial transverse sections across the white matter streamlines corresponding to each link, through an in-house modified version of the methods described by Garyfallidis et al (). The centroids for each link were connected to define the main pathway of streamline travel. This step was repeated for all possible links. The utility of this approach was exclusively for anatomical visualization of the results, as demonstrated in Figure 3.

Figure 3

Discussion

In this study, we investigated the individual variability of hemispheric asymmetry in structural connectivity in healthy senior individuals. By reconstructing the structural connectome from each subject, and by assessing the distribution of asymmetry indices related to cortical ROI connectivity and link weight, we observed that, within a sample composed of healthy older adults, there was mild but noticeable hemispheric asymmetry in structural connectivity.

We observed a more prominent leftward asymmetry in cortical connectivity, i.e., a larger number of regions demonstrated a higher degree of connectivity in the left hemisphere. Specifically, cortical ROIs located in the occipital lobe, medial temporal, dorsolateral and medial frontal lobe and cingulate exhibited a higher weight of connectivity in the left hemisphere. Conversely, fewer cortical regions demonstrated a rightward asymmetry, with the middle temporal gyrus and the orbito-frontal regions exhibiting a significantly higher cortical connectivity on the right hemisphere.

In turn, we also observed regional hemispheric asymmetry in relationship with the strength of connectivity between specific ROIs. In accordance with the previous observation about global cortical connectivity, a leftward asymmetry was also more commonly observed among links. Interestingly, a leftward asymmetry was observed on peri-Sylvian and medial temporal – occipital regions. Conversely, a rightward asymmetry was noted on lateral temporal – frontal – occipital and parietal links. While did not test the relationship between these links and cognitive performance, it is possible to speculate that some lateral asymmetry may be related to functional specialization of some of these connections. For instance, leftward asymmetry may be associated with dominant hemisphere language processing involving verbal memory (parahippocampal – amydala connections) (), phonological processing (insula – frontal operculum connections) (), and semantic retrieval () (angular gyrus – superior parietal connections). Conversely, a rightward asymmetry may be observed in networks associated with visual-spatial processing () (intraparietal connections, precuneus-parietal connections).

It should be noted that we adopted a liberal statistical threshold (i.e., the level of statistical significance from the one-sample t-tests was not corrected based on multiple comparisons). We adopted a liberal threshold since we expected that the degree of asymmetry exhibited by our population would be mild, and given the number of multiple comparisons, a more stringent threshold would preclude the evaluation of the locations with a higher degree of asymmetry. Nonetheless, given our sample size, it is possible that some of our observed asymmetries may be related to sample bias and may constitute false positives. For this reason, we recommend the interpretation of the results from this manuscript in this context. The asymmetry index is possibly a better representation of the magnitude of asymmetry, rather than over-emphasizing the importance of links or ROIs with p < 0.05.

The results reported in this study should help the contextual evaluation of other connectome studies applied to health and disease. We demonstrated the range of asymmetry in a small cohort of normal older adults with the purpose of providing a reference for future studies evaluating processes that affect neural network organization. Thus, future studies should also be interpreted with special attention to specific characteristics of the demographics from the population studied.

Statements

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Appendix

Table A1

NumberSideROINumberSideROI
1LeftPrecentral46RightPrecentral
2LeftFrontal sup47RightFrontal sup
3LeftFrontal sup orb48RightFrontal sup orb
4LeftFrontal mid49RightFrontal mid
5LeftFrontal mid orb50RightFrontal mid orb
6LeftFrontal inf oper51RightFrontal inf oper
7LeftFrontal inf tri52RightFrontal inf tri
8LeftFrontal inf orb53RightFrontal inf orb
9LeftRolandic oper54RightRolandic oper
10LeftSupp motor area55RightSupp motor area
11LeftOlfactory56RightOlfactory
12LeftFrontal sup medial57RightFrontal sup medial
13LeftFrontal mid orb58RightFrontal mid orb
14LeftRectus59RightRectus
15LeftInsula60RightInsula
16LeftCingulum ant61RightCingulum ant
17LeftCingulum mid62RightCingulum mid
18LeftCingulum post63RightCingulum post
19LeftHippocampus64RightHippocampus
20LeftParahippocampal65RightParahippocampal
21LeftAmygdala66RightAmygdala
22LeftCalcarine67RightCalcarine
23LeftCuneus68RightCuneus
24LeftLingual69RightLingual
25LeftOccipital sup70RightOccipital sup
26LeftOccipital mid71RightOccipital mid
27LeftOccipital inf72RightOccipital Inf
28LeftFusiform73RightFusiform
29LeftPostcentral74RightPostcentral
30LeftParietal sup75RightParietal sup
31LeftParietal inf76RightParietal inf
32LeftSupramarginal77RightSupramarginal
33LeftAngular78RightAngular
34LeftPrecuneus79RightPrecuneus
35LeftParacentral lobule80RightParacentral lobule
36LeftCaudate81RightCaudate
37LeftPutamen82RightPutamen
38LeftPallidum83RightPallidum
39LeftThalamus84RightThalamus
40LeftHeschl85RightHeschl
41LeftTemporal sup86RightTemporal sup
42LeftTemporal pole sup87RightTemporal pole sup
43LeftTemporal mid88RightTemporal mid
44LeftTemporal pole mid89RightTemporal pole mid
45LeftTemporal inf90RightTemporal inf

Legends for ROI numbering.

Summary

Keywords

asymmetry, structural connectome, connectivity, diffusion tensor imaging, aging

Citation

Bonilha L, Nesland T, Rorden C and Fridriksson J (2014) Asymmetry of the Structural Brain Connectome in Healthy Older Adults. Front. Psychiatry 4:186. doi: 10.3389/fpsyt.2013.00186

Received

01 October 2013

Accepted

23 December 2013

Published

09 January 2014

Volume

4 - 2013

Edited by

Qiyong Gong, West China Hospital of Sichuan University, China

Reviewed by

James Gee, University of Pennsylvania, USA; Su Lui, West China Hospital of Sichuan University, China

Copyright

*Correspondence: Leonardo Bonilha, Department of Neurology, Comprehensive Epilepsy Center, Medical University of South Carolina, 96 Jonathan Lucas Street, Suite 301 CSB, Charleston, SC 29425, USA e-mail:

This article was submitted to Neuropsychiatric Imaging and Stimulation, a section of the journal Frontiers in Psychiatry.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics