Microstructurally Informed Subject-Specific Parcellation of the Corpus Callosum using Axonal Water Fraction

doi:10.21203/rs.3.rs-3645723/v1

Download PDF

Short Report

Microstructurally Informed Subject-Specific Parcellation of the Corpus Callosum using Axonal Water Fraction

https://doi.org/10.21203/rs.3.rs-3645723/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

The corpus callosum (CC) is the most important interhemispheric white matter (WM) structure composed of several anatomically and functionally distinct WM tracts. Resolving these tracts is a challenge since the callosum appears relatively homogenous in conventional structural imaging. Commonly used callosal parcellation methods such as the Hofer/Frahm scheme rely on rigid geometric guidelines to separate the substructures that are limited to consider individual variation. Here we present a novel subject-specific and microstructurally-informed method for callosal parcellation based on axonal water fraction (ƒ) known as a diffusion metric reflective of axon caliber and density. We studied 30 healthy subjects from the Human Connectome Project (HCP) dataset with multi-shell diffusion MRI. The biophysical parameter ƒ was derived from compartment-specific WM modeling. Inflection points were identified where there were concavity changes in ƒ across the CC to delineate callosal subregions. We observed relatively higher ƒ in anterior and posterior areas consisting of a greater number of small diameter fibers and lower ƒ in posterior body areas of the CC consisting of a greater number of large diameter fibers. Based on degree of change in ƒ along the callosum, seven callosal subregions can be consistently delineated for each individual. We observe that ƒ can capture differences in underlying tissue microstructures and seven subregions can be identified across CC. Therefore, this method provides microstructurally informed callosal parcellation in a subject-specific way, allowing for more accurate analysis in the corpus callosum.

corpus callosum

parcellation

axonal water fraction

diffusion MRI

microstructure

The corpus callosum (CC) is the largest commissural fiber bundle with more than 200 million axons, connecting left and right hemispheres of the brain. It is critical to interhemispheric communication and the global transfer of information across the brain. The CC is important in neurological disease and commonly involves conditions such as traumatic brain injury (Levin et al. 2000), demyelinating disorders (Evangelou et al. 2000; Ge et al. 2004) and brain tumor extension (Ho et al. 2013). Also, it has been implicated in a variety of learning, behavioral and affective disorders such as dyslexia (von Plessen et al. 2002), schizophrenia (Narr et al. 2002), depression (Lacerda et al. 2005) and autism (Prigge et al. 2013).

Commonly, the CC is divided into 3–5 arbitrary subregions from anterior to posterior to facilitate anatomic and imaging studies because different subregions connect specific cortical regions, serve different purposes, and have different propensities for disease. Unfortunately, there is a lack of clear boundaries that allow in vivo imaging parcellation of the CC since the structure appears relatively homogenous on conventional structural imaging. Most common research approaches rely on arbitrary, gross geometry to separate the CC into parts: for example, the Witelson scheme (Witelson 1989) defines subdivisions at 1/3, 1/2, 2/3, and 4/5 along the total callosal length anterior to posterior that is based on a mix of non-human primate and human datasets; Hofer and Frahm (Hofer and Frahm 2006), on the other hand, suggested subdivisions at 1/6, 1/2, 2/3, and 3/4 along the callosal length basing their system on presumed cortical projections derived from DTI-based fiber tractography (Basser and Pierpaoli 1996). However, such schemes ignore any individual variation in callosal structure which is clearly present even at a macroscopic level. Moreover, differences in fiber composition across the CC are observed on microscopy of human autopsy specimens, revealing mainly higher density of small and mid-diameter axons (< 2 µm in diameter) at anteriorly, large axons (> 2 µm in diameter) in the posterior mid-body, and mixed-size axons most posteriorly (Aboitiz et al. 1992).

Recent advances in compartment-specific white matter (WM) modeling of multi-shell diffusion MRI provide promising imaging markers that reflect underlying tissue microstructure (Fieremans et al. 2011; Novikov et al. 2018). In particular, axonal water fraction (ƒ) representing the volume of intra-axonal water relative to total intra and extra-axonal water volume, is known to reflect axon caliber and axon density based on animal and human studies with pathologic correlate (Jelescu et al. 2016; Margoni et al. 2019; Barazany et al. 2009). Distributions of varying axon caliber and density are observed in different sectors of the CC using electron microscopy (Aboitiz et al. 1992). Thus, we hypothesize that ƒ is able to capture differences in fiber composition as they vary along the anteroposterior extent of the midline CC. In this study, we present a novel subject-specific method for callosal parcellation based on the biophysical parameter ƒ and compare the results against callosal subregions as defined by the commonly used Hofer and Frahm scheme (Hofer and Frahm 2006).

Subjects

This study includes 30 healthy subjects (age range, 22–35 years; 15 males) with multi-shell diffusion images from the Human Connectome Project (HCP) (Van Essen et al. 2013), an open-access multi-center dataset with high-quality 3T MR images. Institutional review board approval and participants’ informed consent were obtained at the participating institutions. Subjects have no documented history of mental illness, neurological disorder, or physical illness.

Diffusion MRI Data

MRI data were acquired on two 3T Connectome scanners (Skyra, Siemens, Erlangen, Germany) with a 32-channel head coil. Diffusion imaging was performed with b-value = 1000, 2000, 3000 s/mm² along 90 diffusion encoding directions for each, using multiband (factor of three) spin-echo EPI with 6/8 partial fourier. For this study, we used only b-values up to 2000 s/mm² since high b-values are typically thought to distort quadratic fitting of diffusion kurtosis imaging (DKI) (Jensen et al. 2005). Other imaging parameters are: FOV = 210 mm × 210 mm, resolution = 1.25 × 1.25 × 1.25 mm³, matrix = 168 × 144, 111 slices, TR/TE = 5520/89.5 ms, bandwidth = 1488 Hz/pixel.

Image Processing and Analyses

Diffusion MRI Processing

We used datasets that were preprocessed including b0 image intensity normalization, echo planar imaging distortion correction (FSL’s function, ‘topup’), eddy current and motion correction (FSL’s function, ‘eddy’), and gradient nonlinearity correction (Glasser et al. 2013).

Axonal water fraction, ƒ, is calculated by ƒ = K_max/(K_max+3), where K_max is the maximum kurtosis overall diffusion directions (Fieremans et al. 2011) using in-house image processing software.

Corpus Callosum Parcellation and Evaluation

Parcellation of the CC followed the following procedure (Fig. 1): 1) obtaining callosal masks in the midsagittal plane using the JHU ICBM-DTI-81 WM labels atlas (Mori et al. 2008); 2) generating the callosal centerline by using FSL’s ‘tbss_skeleton’ commend (Lee et al. 1994); 3) plotting mean ƒ value of voxels perpendicular to the local centerline along the CC; and 4) identification of inflection points where there are concavity changes in ƒ (i.e., where the second derivative of ƒ = 0) to delineate callosal subregions based on the highest gradient in ƒ across the CC.

The results were compared with the Hofer and Frahm’s geometric partitioning scheme (Hofer and Frahm 2006) placing subdivisions at 1/6, 1/2, 2/3 and 3/4 along the callosal length for each subject. In addition, we employed whole-brain fiber tractography using MRtrix3 (Brain Research Institute, Melbourne, Australia) (Tournier et al. 2019) to map the cortical projections of the CC subregions derived from our method. Tracking parameters included 10 million streamlines with second-order integration over fiber orientation distributions (iFOD2), a step size of 0.6 mm, and a minimum length of track of 6.25 mm. Among whole-brain fiber tracts, only streamlines that traversed each callosal subregion were selected. The projected cortical regions were referred to by Brodmann area number (Brodmann 2006).

Figure 2A shows that ƒ is relatively higher in anterior and posterior regions and lower in the posterior body region of the CC in representative subjects. Six inflection points of ƒ are consistently found leading to the delineation of seven subregions (Fig. 2B). On average, borders between subregions are identified at 1/8, 1/3, 1/2, 2/3, 5/7 and 4/5 of the total callosal length going from anterior to posterior. As shown in Fig. 2B, our borders (black arrowhead) do differ from Hofer and Frahm’s geometric borders (top white bar) for all individuals. The pattern of ƒ is consistent throughout all 30 healthy subjects and the value of ƒ ranges from 0.38 to 0.81 in the CC (Fig. S1).

Figure 3 shows cortical projections of the callosal fiber bundles from the seven subregions, including prefrontal (A1: Brodmann 10–11), frontal (A2: Brodmann 9), supplementary motor (A3: Brodmann 8), premotor (A4: Brodmann 6), motor (A5: Brodmann 4), sensory (A6: Brodmann 1–3), and parietal (A7: Brodmann 7) combined with temporal and occipital regions, which are consistent with previous reports (Fabri et al. 2014; Park et al. 2008). The callosal parietal, temporal and occipital fiber bundles could not be separated.

We present a novel callosal parcellation method using a biophysical parameter ƒ derived from advanced WM modeling based on multi-shell diffusion MRI. This study shows that ƒ does indeed capture differences in underlying tissue microstructures that are likely reflective of axon caliber and density across the CC ventrodorsally. We observe relatively higher ƒ in the anterior regions that contain about 72% of thinner, lightly myelinated fibers with small diameters of 0.2–1 µm in the total fiber population, and lower ƒ in the posterior mid-body regions that contain a great number of thicker, heavily myelinated fibers with relatively larger diameter (> 3 µm) as reported using electron microscopy by Aboitiz, et al (Aboitiz et al. 1992).

Current commonly used geometric parcellation schemes relying on fixed partitioning ratios such as the Witelson and Hofer/Frahm methods are not sensitive to individual variability of underlying tissue compositions. Our method is consistently able to separate WM bundles that track to either side of the marginal sulcus which separates the paracentral lobule from the precuneus, anatomically and functionally distinct areas. This would not always be expected to be the case using simple geometric parcellation methods across all individuals; for example, S5 (Fig. 2B) shows that the Hofer/Frahm scheme combines portions of the sensory sector (blue) and the precuneus (purple) while our method clearly separates them. In addition, our method is able to distinguish specific callosal subregions associated with frontal (A2: Brodmann 9), supplementary motor (A3: Brodmann 8) and premotor (A4: Brodmann 6) areas that are previously conglomerated in the Hofer/Frahm scheme. Interestingly, Fig. 3 shows the supplementary motor area (SMA) (A3; yellow) to be the most variable in terms of the cross-sectional area represented in the CC among these callosal regions. The SMA is known to be variable across individuals without good predictors of who may suffer SMA syndrome after frontal lobe surgery or insult (Ribas 2010; Baker et al. 2018). Being able to accurately parse cross-hemisphere tracts that contribute to specific anatomic and functional pathways could help guide future exploration such as the role of supplementary motor and premotor tracts in motor function and dysfunction.

Some limitations of the proposed method include being incumbent on adequate spatial resolution as is available in the HCP dataset, particularly in individuals who may have a narrow callosal isthmus (region between the body and splenium) which theoretically may lead to partial volume effects. In addition, the biophysical white matter modeling (Fieremans et al. 2011) used here to determine ƒ makes the assumption of highly aligned tracts. While this is probably a reasonable assumption in the CC as it is a structure with high transverse directionality, it is possible that some of the inability of the approach to resolve parietal, temporal and occipital white matter bundles separately within the posterior callosum is due to the presence of crossing fibers. Finally, this is a study of healthy adults; pathology of the CC may make parcellation difficult, though this in and of itself could serve as important information.

We present a novel callosal parcellation method based on ƒ, reflective of biophysical factors of the underlying callosal microstructure. This method provides microstructurally-informed callosal subregions in a subject-specific way, allowing for more biologically-based analysis of regions of the CC.

corpus callosum

axonal water fraction

HCP

Human Connectome Project

DKI

diffusion kurtosis imaging

iFOD2

second-order integration over fiber orientation distributions

SMA

supplementary motor area

Funding: This work was supported in part by grant funding from the National Institute of Health (NIH): R01 NS119767-01A1, R01 NS039135-11, R56 NS119767. This work was also performed under the rubric of the Center for Advanced Imaging Innovation and Research (CAI2R), an NIBIB Biomedical Technology Resource Center (NIH P41 EB017183).

Competing Interests:

The authors have no relevant financial or non-financial interests to disclose.

Author Contributions:

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Sohae Chung. The first draft of the manuscript was written by Sohae Chung and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Data Availability:

The datasets used in the current study are available in the HCP repository (https://www.humanconnectome.org/study/hcp-young-adult/data-releases).

Aboitiz F, Scheibel AB, Fisher RS, Zaidel E (1992) Fiber composition of the human corpus callosum. Brain Res 598 (1-2):143-153
Baker CM, Burks JD, Briggs RG, Smitherman AD, Glenn CA, Conner AK, Wu DH, Sughrue ME (2018) The crossed frontal aslant tract: A possible pathway involved in the recovery of supplementary motor area syndrome. Brain Behav 8 (3):e00926. doi:10.1002/brb3.926
Barazany D, Basser PJ, Assaf Y (2009) In vivo measurement of axon diameter distribution in the corpus callosum of rat brain. Brain 132 (Pt 5):1210-1220. doi:10.1093/brain/awp042
Basser PJ, Pierpaoli C (1996) Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J Magn Reson B 111 (3):209-219
Brodmann K (2006) Brodmann’s localisation in the cerebral cortex (trans: Garey LJ). Springer, New York
Evangelou N, Konz D, Esiri MM, Smith S, Palace J, Matthews PM (2000) Regional axonal loss in the corpus callosum correlates with cerebral white matter lesion volume and distribution in multiple sclerosis. Brain 123 ( Pt 9):1845-1849
Fabri M, Pierpaoli C, Barbaresi P, Polonara G (2014) Functional topography of the corpus callosum investigated by DTI and fMRI. World J Radiol 6 (12):895-906. doi:10.4329/wjr.v6.i12.895
Fieremans E, Jensen JH, Helpern JA (2011) White matter characterization with diffusional kurtosis imaging. Neuroimage 58 (1):177-188. doi:10.1016/j.neuroimage.2011.06.006
Ge Y, Law M, Johnson G, Herbert J, Babb JS, Mannon LJ, Grossman RI (2004) Preferential occult injury of corpus callosum in multiple sclerosis measured by diffusion tensor imaging. J Magn Reson Imaging 20 (1):1-7. doi:10.1002/jmri.20083
Glasser MF, Sotiropoulos SN, Wilson JA, Coalson TS, Fischl B, Andersson JL, Xu J, Jbabdi S, Webster M, Polimeni JR, Van Essen DC, Jenkinson M, Consortium WU-MH (2013) The minimal preprocessing pipelines for the Human Connectome Project. Neuroimage 80:105-124. doi:10.1016/j.neuroimage.2013.04.127
Ho ML, Moonis G, Ginat DT, Eisenberg RL (2013) Lesions of the corpus callosum. AJR Am J Roentgenol 200 (1):W1-16. doi:10.2214/AJR.11.8080
Hofer S, Frahm J (2006) Topography of the human corpus callosum revisited--comprehensive fiber tractography using diffusion tensor magnetic resonance imaging. Neuroimage 32 (3):989-994. doi:10.1016/j.neuroimage.2006.05.044
Jelescu IO, Zurek M, Winters KV, Veraart J, Rajaratnam A, Kim NS, Babb JS, Shepherd TM, Novikov DS, Kim SG, Fieremans E (2016) In vivo quantification of demyelination and recovery using compartment-specific diffusion MRI metrics validated by electron microscopy. Neuroimage 132:104-114. doi:10.1016/j.neuroimage.2016.02.004
Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K (2005) Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn Reson Med 53 (6):1432-1440. doi:10.1002/mrm.20508
Lacerda AL, Brambilla P, Sassi RB, Nicoletti MA, Mallinger AG, Frank E, Kupfer DJ, Keshavan MS, Soares JC (2005) Anatomical MRI study of corpus callosum in unipolar depression. J Psychiatr Res 39 (4):347-354. doi:10.1016/j.jpsychires.2004.10.004
Lee TC, Kashyap RL, Chu CN (1994) Building skeleton models via 3-D medial surface axis thinning algorithms. CVGIP: Graphical models and image processing 56 (6):462-478
Levin HS, Benavidez DA, Verger-Maestre K, Perachio N, Song J, Mendelsohn DB, Fletcher JM (2000) Reduction of corpus callosum growth after severe traumatic brain injury in children. Neurology 54 (3):647-653
Margoni M, Petracca M, Schiavi S, Fabian M, Miller A, Lublin FD, Inglese M (2019) Axonal water fraction as marker of white matter injury in primary-progressive multiple sclerosis: a longitudinal study. Eur J Neurol 26 (8):1068-1074. doi:10.1111/ene.13937
Mori S, Oishi K, Jiang H, Jiang L, Li X, Akhter K, Hua K, Faria AV, Mahmood A, Woods R, Toga AW, Pike GB, Neto PR, Evans A, Zhang J, Huang H, Miller MI, van Zijl P, Mazziotta J (2008) Stereotaxic white matter atlas based on diffusion tensor imaging in an ICBM template. Neuroimage 40 (2):570-582. doi:10.1016/j.neuroimage.2007.12.035
Narr KL, Cannon TD, Woods RP, Thompson PM, Kim S, Asunction D, van Erp TG, Poutanen VP, Huttunen M, Lonnqvist J, Standerksjold-Nordenstam CG, Kaprio J, Mazziotta JC, Toga AW (2002) Genetic contributions to altered callosal morphology in schizophrenia. J Neurosci 22 (9):3720-3729. doi:20026309
Novikov DS, Veraart J, Jelescu IO, Fieremans E (2018) Rotationally-invariant mapping of scalar and orientational metrics of neuronal microstructure with diffusion MRI. Neuroimage 174:518-538. doi:10.1016/j.neuroimage.2018.03.006
Park HJ, Kim JJ, Lee SK, Seok JH, Chun J, Kim DI, Lee JD (2008) Corpus callosal connection mapping using cortical gray matter parcellation and DT-MRI. Hum Brain Mapp 29 (5):503-516. doi:10.1002/hbm.20314
Prigge MB, Lange N, Bigler ED, Merkley TL, Neeley ES, Abildskov TJ, Froehlich AL, Nielsen JA, Cooperrider JR, Cariello AN, Ravichandran C, Alexander AL, Lainhart JE (2013) Corpus Callosum Area in Children and Adults with Autism. Res Autism Spectr Disord 7 (2):221-234. doi:10.1016/j.rasd.2012.09.007
Ribas GC (2010) The cerebral sulci and gyri. Neurosurg Focus 28 (2):E2. doi:10.3171/2009.11.FOCUS09245
Tournier JD, Smith R, Raffelt D, Tabbara R, Dhollander T, Pietsch M, Christiaens D, Jeurissen B, Yeh CH, Connelly A (2019) MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. Neuroimage 202:116137. doi:10.1016/j.neuroimage.2019.116137
Van Essen DC, Smith SM, Barch DM, Behrens TE, Yacoub E, Ugurbil K, Consortium WU-MH (2013) The WU-Minn human connectome project: an overview. Neuroimage 80:62-79. doi:10.1016/j.neuroimage.2013.05.041
von Plessen K, Lundervold A, Duta N, Heiervang E, Klauschen F, Smievoll AI, Ersland L, Hugdahl K (2002) Less developed corpus callosum in dyslexic subjects--a structural MRI study. Neuropsychologia 40 (7):1035-1044
Witelson SF (1989) Hand and sex differences in the isthmus and genu of the human corpus callosum. A postmortem morphological study. Brain 112 ( Pt 3):799-835

No competing interests reported.

FigureS1.docx

Download PDF

Editorial decision: Revision requested
29 Dec, 2023
Reviews received at journal
26 Dec, 2023
Reviewers agreed at journal
06 Dec, 2023
Reviewers invited by journal
04 Dec, 2023
Editor assigned by journal
03 Dec, 2023
Submission checks completed at journal
22 Nov, 2023
First submitted to journal
21 Nov, 2023

You are reading this latest preprint version

Microstructurally Informed Subject-Specific Parcellation of the Corpus Callosum using Axonal Water Fraction

Status:

Version 1

Abstract

Figures

INTRODUCTION

MATERIALS AND METHODS

Subjects

Diffusion MRI Data

Image Processing and Analyses

Diffusion MRI Processing

Corpus Callosum Parcellation and Evaluation

RESULTS

DISCUSSION

CONCLUSIONS

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1