Skip to main content

METHODS article

Front. Surg., 29 July 2021
Sec. Neurosurgery

Constrained-Spherical Deconvolution Tractography in the Evaluation of the Corticospinal Tract in Glioma Surgery

  • 1Juha Hernesniemi International Neurosurgery Center, Zhengzhou University People's Hospital (Henan Provincial People's Hospital), Zhengzhou, China
  • 2Department of Radiology, People's Hospital of Zhengzhou University (Henan Provincial People's Hospital), Zhengzhou, China
  • 3Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, ON, Canada
  • 4Multi Scale Robotics Laboratory, ETH Zurich, Zurich, Switzerland
  • 5Department of Neurosurgery, Center for Neuromodulation, Mount Sinai Health System, New York, NY, United States

Introduction: Tractography has demonstrated utility for surgical resection in the setting of primary brain tumors involving eloquent white matter (WM) pathways.

Methods: Twelve patients with glioma in or near eloquent motor areas were analyzed. The motor status was recorded before and after surgery. Two different tractography approaches were used to generate the motor corticospinal tract (CST): Constrained spherical deconvolution probabilistic tractography (CSD-Prob) and single tensor deterministic tractography (Tens-DET). To define the degree of disruption of the CST after surgical resection of the tumor, we calculated the percentage of the CST affected by surgical resection, which was then correlated with the postoperative motor status. Moreover, the fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) of the CST generated by the CSD-Prob and the Tens-DET was measured and compared between the ipsilesional and contralesional side.

Results: The CST was identified in all patients and its trajectory was displaced by the tumor. Only the CSD-Prob approach showed the CST with the characteristic fan-like projections from the precentral gyrus to the brainstem. Disruption of the CST was identified in 6/6 with postoperative motor deficit by CSD-Prob approach and in 5/6 in the Tens-DET. The degree of disruption was significantly associated with the motor deficit with the CSD-Prob approach (rho = −0.88, p = 0.021). However, with the Tens-DET approach the CST disruption did not show significant association with the motor function (rho = −0.27, p = 0.6). There was a significant decrease in FA (p = 0.006) and a significant increase in MD (p = 0.0004) and RD (p = 0.005) on the ipsilesional CST compared with the contralesional CST only with the CSD-Prob approach.

Conclusion: CSD-Prob accurately represented the known anatomy of the CST and provided a meaningful estimate of microstructural changes of the CST affected by the tumor and its macrostructural damage after surgery. Newer surgical planning stations should include advanced models and algorithms of tractography in order to obtain more meaningful reconstructions of the WM pathways during glioma surgery.

Introduction

Surgical resection is part of the multidisciplinary treatment of gliomas and the main goal is to achieve a safe gross total resection (1). When these tumors are located near or within eloquent brain regions, the surgical resection poses significant risk of damaging those regions. In these cases, mapping the eloquent white matter (WM) pathways related with the tumor has been considered useful to plan the surgical approaches (2). Diffusion-weighted image (DWI) is a magnetic resonance image (MRI) sequence that is based on the movement of the water molecules in every voxel of brain tissue. In the white matter, due to the increased hindrance to diffusion across the axonal membranes, the movement of the water molecules detected by DWI follow the directions of the major axonal bundles. Tractography is a technique that models the white matter pathways in the brain using DWI data. Tractography has demonstrated to be useful for surgical resection in the setting of primary brain tumors in eloquent areas, as long as the WM pathways generated are anatomically accurate (3).

The single-tensor model is the most commonly used to estimate WM fiber orientations from DWI data. This model defines a single predominant orientation of the movement of water molecules in each voxel. Likewise, the deterministic algorithm is the tracking method based on the principal diffusion direction modeled voxel by voxel (46). The combination of the single-tensor model with the deterministic algorithm is one of the first techniques used for tractography and is currently the most commonly used in surgical planning stations. The effect of intra-axial tumors on the WM (edema, infiltration, displacement, etc.) may affect the ability of less contemporary tractography techniques such as single-tensor deterministic to provide accurate reconstructions of the WM pathways (7). Ideally, the streamlines generated by tractography, which is the 3D representation of the tract of interest, should accurately define the trajectory and roughly approximate the cross-sectional area of the real underlying WM pathway. This way, a reliable plan for safe resection can be tailored based on the relationship between the tumor and the WM pathways.

Meanwhile, high-order tractography methods have been developed in the past 20 years. These approaches have demonstrated to improve the anatomical accuracy and the validity of the results (8). Constrained-spherical deconvolution (CSD) tractography has been previously evaluated in patients with intra-axial brain tumors. This approach provided valuable quantitative and qualitative information that has not been obtained with single-tensor approaches to aid preoperative planning for these patients (9). When comparing deterministic and probabilistic tractography algorithms, high-order tractography (q-ball probabilistic) techniques demonstrated better sensitivity and specificity than tensor-based approaches to locate the functional motor pathways using intraoperative stimulation as the gold standard (10). These data also evidenced the very poor sensitivity of the current tensor-based approaches in finding the ventral aspects of the motor corticospinal tract. The effect of edema and tumor mass effect further reduced the sensitivity of these approaches.

Considering that the aim of tractography for brain tumor surgery is to obtain meaningful reconstructions of the major WM pathways, there is still controversy on the advantage of advanced approaches in terms of its utility in clinical practice. Clinically meaningful reconstructions are urgently needed as image-guided minimally invasive approaches (laser ablation, focused ultrasound, radiosurgery, etc.) are becoming more popular in the management of tumors in eloquent areas or deeper locations (11, 12). The aim of this study is to investigate the utility of high-order tractography approaches in patients with intra-axial tumors and their clinical outcomes. To this aim, we analyzed patients operated for gliomas near or within motor areas. We compared the single-tensor deterministic approach (Tens-DET), which is the available approach in surgical planning stations with the CSD probabilistic approach (CSD-Prob) as a high-order more advanced tractography approach.

Methods

Patient Selection

Fifty-three patients underwent glioma resection surgery with subsequent pathological confirmation according to the 2016 WHO brain tumor classification criteria (13) at the Henan Provincial People's Hospital between January 2017 and January 2019. From these 53 patients, we selected patients based on the following criteria: tumor resection surgery, tumors within or near motor eloquent areas, and mild motor deficit (4/5) at the clinical presentation that would suggest displacement rather than damage of the motor corticospinal tract (CST). Additionally, we excluded patients if they had significant postoperative brain shift or edema that would preclude accurate co-registration of the preoperative and postoperative acquisitions. We recorded the clinical and demographic data, extent of resection, and permanent postoperative motor deficits (no recovery after last follow-up ~2 months after surgery). The motor deficit was evaluated with the medical research council scale ranging from 1 (no contraction) to 5 (normal muscle strength). Gross total resection was defined as resection of more than 90% of the tumor volume. The local ethics board approved the study and consent was provided by all patients.

Data Acquisition, Pre-processing, and Registration

The MRI data were acquired using a 3.0T whole-body MRI scanner (Siemens electronic medical system). The 3-dimensional preoperative and postoperative (within 2 days after surgery) T1-weighted images with gadolinium were acquired with the following imaging parameters: repetition time = 2,300 ms, echo time = 2.3 ms, inversion time = 900 ms, slice thickness = 1 mm, flip angle = 8, number of excitations = 1, field of view = 240 * 240 mm2, voxel sizes = 0.9 * 0.9 * 1.0 mm3, and acquisition time = 421 s. Preoperative high-angular resolution diffusion images (HARDI) were acquired with a thickness of 4.0 mm that covers the whole brain with the following parameters: repetition time = 3,000 ms, echo time = 69 ms, flip angle = 180, b-values = 1,000 s/m2, field of view = 220 * 220 mm2, voxel sizes = 1.7 * 1.7 * 4.0 mm3, number of excitations = 4, contiguous slices = 25, acquisition time = 257 s, and 64 gradient directions. The data was pre-processed using MRtrix and FSL tools (14, 15). Eddy-current, bias-field, and subject movement were corrected (16). The response function and the spherical deconvolution were estimated, and fiber orientation distributions (FOD) were obtained. Finally, linear registration was utilized to co-register the preoperative (T1-weighted and diffusion images) with the postoperative acquisitions. Using the postoperative acquisitions, we contoured the resection cavity, and this volume was later used to calculate the disruption of the CST by the surgical resection.

Tractography Analysis

The CST was generated using two different approaches: the CSD-probabilistic tractography (CSD-Prob) and the single tensor deterministic tractography (Tens-DET). We defined three regions of interest (ROIs) based on the neuroanatomy of the CST as follows (17): Seed ROI: The precentral gyrus. Inclusion ROIs (10 mm): one sphere at the posterior half of the posterior limb of the internal capsule and one at the posterolateral third of the cerebral peduncle (Figure 1). The fiber tracking was performed in the ipsilesional and contralesional side with 10,000 streamlines, cut-off: 0.15 (FOD amplitude for CSD-Prob and FA for Tens-DET), step size: 0.5 × voxelsize, angle threshold: 45, and seeding density 103 × the number of selected streamlines. To determine the microstructural changes of the CST as a result of the proximity of the tumor, we calculated a series of diffusion scalars including fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) along the ipsilesional and contralesional CST. To this aim, we used the map of the motor CST to sample the DTI scalars along the CST generated by the Tens-DET and the CSD-Prob. These parameters are indirect measures of white matter structural integrity (14). Finally, the degree of surgical damage of the CST, henceforth referred to as “tract disruption,” was defined as follows: We first generated the ipsilesional corticospinal tract (CST) of each patient with a fixed arbitrary number of streamlines (10,000), which are the 3D representation of the trajectory of a given WM bundle (5). We used the postoperative images to segment the tumor resection cavity. Then, we subtracted the CST streamlines passing through the cavity of the tumor using MRtrix tools. We then defined the disruption of the CST as the relative percentage of streamlines excluded by the resection cavity. Finally, we used this percentage to correlate with the motor status.

FIGURE 1
www.frontiersin.org

Figure 1. Study outline. First, diffusion weighted images were preprocessed with subject-motion, Eddy current distortion, and bias field correction. Fiber orientation distributions were calculated with constrained-spherical deconvolution and the directionally-encoded colors were obtained. The precentral gyrus was delineated using the T1 and defined as the seed region of interest (ROI). The other inclusion ROIs were selected in the directionally-encoded colors sequence based on the described trajectory of the motor CST. The posterior half of the posterior limb of the internal capsule and the posterolateral third of the cerebral peduncle (white arrows).

Statistical Analyses

The diffusion scalars calculated from the CST (FA, MD, RD, and AD) are values between 0 to 1 and were compared between the ipsilesional and contralesional side using the Wilcoxon test. The correlation analysis of the CST disruption with the postoperative motor status was performed using the Spearman correlation test. We used the SPSS software package 22.0 (IBM Corp, Armonk, NY) and the significance level was set at p < 0.05.

Results

We included 12 patients, six men and eight women with a mean age of 50 years (detailed clinical and demographic data are shown in Table 1). Pathological results yielded WHO grade II glioma in four patients and WHO IV glioblastoma (GBM) in eight patients. Preoperatively, six patients were neurologically intact and six presented with mild motor deficit (grade 4/5). Gross total resection was achieved in 9 patients. The CST was identified in all patients with the two tractography approaches with different patterns of displacement (Figure 2). The CSD-Prob approach showed the CST with the characteristic fan-like configuration from the dorsal and ventral precentral gyrus to the brainstem. The Tens-DET approach showed a straight trajectory from the dorsal precentral gyrus to the brainstem missing the ventral part of the precentral gyrus. This approach also underestimated the width of the CST and showed more spurious streamlines as compared with the CSD-Prob approach (Figure 3). Among the six patients with postoperative motor deficit, the CSD-Prob approach revealed CST disruption in different degree, and this was strongly associated with the degree of motor deficit in these patients (rho = −0.88, p = 0.021). The Tens-DET approach on the other hand, showed disruption of the CST in five of the six patients with motor deficit without significant association in the motor status (rho = −0.27, p = 0.6) (Table 1, Figure 4). Finally, only the CSD-Prob identified a significant decrease in fractional anisotropy (p = 0.0006) and an increase in mean and radial diffusivity (p = 0.004, 0.005) of the CST between the ipsilesional and the contralesional hemisphere (Figure 5). There was no significant difference in the axial diffusivity with the CSD-Prob approach and with none of the diffusion scalars with the Tens-DET approach.

TABLE 1
www.frontiersin.org

Table 1. Clinical and demographics data of the patients.

FIGURE 2
www.frontiersin.org

Figure 2. Motor CST displacement. (A), Axial, coronal, and 3D view of a medial displacement of the CST. (B), Axial, sagittal, and 3D view of a posterior displacement of the CST. (C), Axial, sagittal, and 3D view of an anterior displacement of the CST.

FIGURE 3
www.frontiersin.org

Figure 3. Comparison of the CSD-Prob and Tens-DET approaches in 5 subjects. The left panels show the 3D reconstruction of the CST by the CSD-Prob with the fan-like configuration of the cortical projections of the CST. The right panels show the CST by the Tens-DET approach with a straight trajectory to the dorsal part of the precentral gyrus, absence of the ventral part of the motor CST, and more prevalence of spurious fibers than the CSD-Prob (arrows).

FIGURE 4
www.frontiersin.org

Figure 4. Correlation analyses show that the degree of motor deficit had a strong negative correlation with the CST disruption by the tumor resection. This significant association was found with the CSD-Prob approach and not with the Tens-DET approach likely due to the missing parts of the ventral part of the precentral gyrus.

FIGURE 5
www.frontiersin.org

Figure 5. Boxplots showing statistically significant differences in fractional anisotropy, mean diffusivity, and radial diffusivity between the ipsilesional and contralesional side only with the CSD-Prob approach. **P ≤ 0.01.

Discussion

The CST was generated by two different tractography techniques in the setting of brain tumor and peritumoral edema. The CSD-Prob approach has demonstrated accuracy in the representation of the CST as previously demonstrated (9). The CSD-Prob also showed a meaningful reconstruction of the CST since there was a consistent degree of CST disruption significantly correlated with postoperative motor deficit. The Tens-DET approach showed limited accuracy to generate the CST with more spurious streamlines than the CSD-Prob approach. Most of these spurious streamlines originated at the level of the brainstem and continued propagating to the contralateral side following different bundles including contralateral CST, cerebellum, and optic radiations. This may be related with a more prevalence of crossing fibers in the brainstem and the difficulty of the Tens-DET to resolve crossing fibers. Modifications in the tracking parameters of the Tens-DET approach could have minimized the presence of these spurious streamlines but the comparison of different tracking parameters in the Tens-DET approaches is beyond the scope of this article. In three of the six cases with postoperative motor deficit, the disruption of the CST defined by the Tens-DET approach was not consistent with the degree of motor deficit, likely due to the missing parts of the ventral precentral gyrus. Also, only the CSD-Prob approach identified a significant decrease in FA and a significant increase in the MD and RD of the CST on the ipsilesional side. The changes in FA, MD, and RD may be interpreted as an alteration of the microstructure in the white matter close to the tumor due to displacement of the axons or increase in free tissue water. These changes are related with tumor infiltration or the effect of the peritumoral edema on the WM tracts (9).

Surgery for gliomas in eloquent areas can be challenging when the extent of resection is compromised by brain function. Patients with mild motor deficit, who usually respond to steroid treatment, may have WM pathways intact and the precise location of these pathways is desirable in order to protect them during surgery. In addition to special techniques such as awake craniotomy (18) or neurophysiological monitorization, the use of tractography is a critical tool for the surgeon to minimize the risk of postoperative deficits during glioma surgery (2, 11, 12). The single-tensor deterministic technique, which is currently the most frequently used algorithm in surgical planning stations, prevents accurate preoperative identification of WM pathways. For example, Kinoshita et al., showed that the tractography reconstruction of the pyramidal tract was not accurate to estimate the cross-sectional area of the pathway despite being “anatomically” correct (19). Other authors on the other hand have found that tractography can accurately predict the position of the WM pathways and correlate with clinical outcomes (20, 21). These discrepancies are associated with major differences in algorithms, tracking settings, software, and data acquisition. For a WM pathway to be meaningful for surgical planning, it should represent the real anatomical trajectory and should closely estimate the cross-sectional diameter. This information would help, not only to understand the relationship between the tumor and the pathway, but also in the trade-off between extent of resection and the risk of postoperative deficit.

Reconstruction of WM pathways using tractography has several considerations to take into account. Usually, the ROIs for fiber tracking are identified based on neuroanatomy. The more ROIs, the more anatomical accuracy is achieved. Because of the distortion of the normal anatomy produced by the tumor, the use of directionally-encoded color (DEC) images would aid in the placement of the ROIs (22, 23). This image provides directional information of the white matter pathways: blue (superior-inferior direction), red (left-right direction), and green (anterior-posterior direction). By placing ROIs in blue areas of the theoretical location of the CST, it is possible to trace the trajectory of the CST by coupling anatomical and directional information. During tracking, single-tensor deterministic approaches use the main fiber orientation in the voxels to generate the pathway and this limits the anatomical accuracy because of the high prevalence of crossing fibers in the WM (24). Another problem of the single-tensor techniques is the inclusion of a number of false positive streamlines depending on the FA threshold established in the software. The present work and others have demonstrated that the fractional anisotropy is affected in the WM close to the tumor because of infiltration or edema. Therefore, false negatives will also limit the results of single-tensor tractography approaches.

CSD probabilistic tractography approaches are not free of limitations, as the sensitivity increases, the specificity of the WM pathway decreases as a substantial amount of false positive streamlines start to appear in the tracking process (25, 26). Current evidence has demonstrated that high-order tractography techniques may improve the false negative rate but have the problem of more false positive streamlines (25). However, the use of anatomical landmarks as waypoints, increases the accuracy in the tractography reconstruction (27, 28). Furthermore, when using probabilistic algorithms the use of thresholding would also improve the accuracy of the results (26). In this work, we used three ROIs using the directional information provided by the directionally-encoded colors sequence. This may help to minimize biases, which would be required to obtain meaningful reconstructions for surgical planning. Despite using clinical acquisitions, our results demonstrated that the CSD-Prob tractography reconstructions closely resemble the known anatomy of the motor CST (8, 29) in the presence of significant edema or displacement of the tracts. The disruption of the CST also showed strong correlation with postoperative motor outcomes. Due to the low sample size and other limitations, our results have to be interpreted with caution. However, we anticipate that the inclusion of high-order tractography approaches in newer surgical planning stations will allow to better characterize its utility for minimally invasive surgical techniques, including laser ablation (30), focused ultrasound, and radiosurgery.

The main limitations of our work include the small sample size and the retrospective nature of this analysis, which prevents a more comprehensive comparison of all tractography techniques available. The aim of this study was to demonstrate that the current methods of tractography used in surgical planning stations are not as reliable as other high-order tractography techniques available. Another limitation is the use of clinical acquisition instead of research acquisitions. Clinical acquisitions are usually limited by timing and patient comfort and therefore subsequent limitations in the acquisition parameters. Our DWI acquisition included highly anisotropic voxels and evidence suggest that this voxel reconstruction could add bias to the tractography reconstructions (31). Also, the diffusion scalars, particularly FA in regions containing crossing fibers, could be underestimated with acquisitions using anisotropic voxels (32). Therefore, it is recommended to use voxels with isotropic resolution in DWI acquisitions to improve the tractography results (33). Another consideration in the use of tractography for surgical planning is the practical meaning of the streamline reconstructions. The absolute number of streamlines lack biological meaning since the streamlines only represent a 3D model of the pathway. Also, the use of the streamlines as a measure of pathway density may be biased by the number of false positive connections (34). In this work, we used the percentage of streamline reduction as a semi-quantitative measure of tract disruption. Some evidence showed that the density of streamlines could be considered a marker of the structural integrity of the underlying WM pathway (35, 36). Streamline tractography has also been used to validate functional connectivity and to explain biological processes such as aging (37, 38).

Conclusion

CSD-Prob accurately represented the known anatomy of the motor CST and provided a meaningful estimate of microstructural changes of the CST affected by the tumor and its macrostructural damage after surgery. Compared with CSD probabilistic approaches, single-tensor deterministic approaches seemed to underestimate the cross-sectional area of the CST in regions close to the tumor. Newer surgical planning stations should include advance tractography approaches in order to obtain more meaningful reconstructions of the WM pathways during glioma surgery.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Henan Provincial People's Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ZS acquired and analyzed the data and wrote the manuscript. JY and YS acquired the data. XB and MW designed the research and acquired the data. CS analyzed the data and wrote the manuscript. JH supervised the research, designed the research, and wrote the manuscript. BN designed the research and wrote the manuscript. AZ and JA-C designed the research, acquired and analyzed the data, supervised the research, wrote the manuscript, and contributed the analytical tools. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the Heidi Demetriades Foundation and the ETH Zurich Foundation to AZ and by the Henan Provincial People's Hospital Outstanding Talents Funding Grant Project.

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.

Publisher's Note

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.

References

1. Sanai N, Berger MS. Glioma extent of resection and its impact on patient outcome. Neurosurgery. (2008) 62:753–64; discussion 264-6. doi: 10.1227/01.neu.0000318159.21731.cf

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Bello L, Castellano A, Fava E, Casaceli G, Riva M, Scotti G, et al. Intraoperative use of diffusion tensor imaging fiber tractography and subcortical mapping for resection of gliomas: technical considerations. Neurosurg Focus. (2010) 28:E6. doi: 10.3171/2009.12.FOCUS09240

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fernandez-Miranda JC, Pathak S, Engh J, Jarbo K, Verstynen T, Yeh FC, et al. High-definition fiber tractography of the human brain: neuroanatomical validation and neurosurgical applications. Neurosurgery. (2012) 71:430–53. doi: 10.1227/NEU.0b013e3182592faa

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Basser PJ, Jones DK. Diffusion-tensor MRI: theory, experimental design and data analysis - a technical review. NMR Biomed. (2002) 15:456–67. doi: 10.1002/nbm.783

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Mori S, van Zijl PCM. Fiber tracking: principles and strategies - a technical review. NMR Biomed. (2002) 15:468–80. doi: 10.1002/nbm.781

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Tournier, J.-D, Mori S, Leemans A. Diffusion tensor imaging and beyond. Magn Reson Med. (2011) 65:1532–56. doi: 10.1002/mrm.22924

CrossRef Full Text | Google Scholar

7. Kuhnt D, Bauer MHA, Egger J, Richter M, Kapur T, Sommer J, et al. Fiber tractography based on diffusion tensor imaging compared with high-angular-resolution diffusion imaging with compressed sensing: initial experience. Neurosurgery. (2013) 72(Suppl. 1):165–75. doi: 10.1227/NEU.0b013e318270d9fb

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Farquharson S, Tournier, J.-D., Calamante F, Fabinyi G, Schneider-Kolsky M, et al. White matter fiber tractography: why we need to move beyond DTI. J Neurosurg. (2013) 118:1367–77. doi: 10.3171/2013.2.JNS121294

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Mormina E, Longo M, Arrigo A, Alafaci C, Tomasello F, Calamuneri A, et al. MRI tractography of corticospinal tract and arcuate fasciculus in high-grade gliomas performed by constrained spherical deconvolution: qualitative and quantitative analysis. AJNR Am J Neuroradiol. (2015) 36:1853–8. doi: 10.3174/ajnr.A4368

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Bucci M, Mandelli ML, Berman JI, Amirbekian B, Nguyen C, Berger MS, et al. Quantifying diffusion MRI tractography of the corticospinal tract in brain tumors with deterministic and probabilistic methods. NeuroImage Clin. (2013) 3:361–8. doi: 10.1016/j.nicl.2013.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Ganau L, Ligarotti GKI, Ganau M. Predicting complexity of tumor removal and postoperative outcome in patients with high-grade gliomas. Neurosurg Rev. (2018) 41:371–3. doi: 10.1007/s10143-017-0921-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Ganau M, Ligarotti GK, Apostolopoulos V. Real-time intraoperative ultrasound in brain surgery: neuronavigation and use of contrast-enhanced image fusion. Quant Imaging Med Surg. (2019) 9:350–8. doi: 10.21037/qims.2019.03.06

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Louis DN, Perry A, Reifenberger G, von Deimling A, Figarella-Branger D, Cavenee WK, et al. The 2016 World Health Organization classification of tumors of the central nervous system: a summary. Acta Neuropathol. (2016) 131:803–20. doi: 10.1007/s00401-016-1545-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Jenkinson M, Beckmann CF, Behrens TEJ, Woolrich MW, Smith SM. FSL. Neuroimage. (2012) 62:782–90. doi: 10.1016/j.neuroimage.2011.09.015

CrossRef Full Text | Google Scholar

15. Tournier JD, Calamante F, Connelly A. MRtrix: diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol. (2012) 22:53–66. doi: 10.1002/ima.22005

CrossRef Full Text | Google Scholar

16. Andersson JLR, Sotiropoulos SN. An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage. (2016) 125:1063–78. doi: 10.1016/j.neuroimage.2015.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Mori S, Wakana S, Zijl PCM, van, Nagae-Poetscher LM. MRI Atlas of Human White Matter. London: Elsevier (2005) 10–20.p.

Google Scholar

18. Rivera-Rivera PA, Rios-Lago M, Sanchez-Casarrubios S, Salazar O, Yus M, González-Hidalgo M, et al. Cortical plasticity catalyzed by prehabilitation enables extensive resection of brain tumors in eloquent areas. J Neurosurg. (2017) 126:1323–33. doi: 10.3171/2016.2.JNS152485

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kinoshita M, Yamada K, Hashimoto N, Kato A, Izumoto S, Baba T, et al. Fiber-tracking does not accurately estimate size of fiber bundle in pathological condition: initial neurosurgical experience using neuronavigation and subcortical white matter stimulation. Neuroimage. (2005) 25:424–9. doi: 10.1016/j.neuroimage.2004.07.076

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chen X, Weigel D, Ganslandt O, Buchfelder M, Nimsky C. Prediction of visual field deficits by diffusion tensor imaging in temporal lobe epilepsy surgery. Neuroimage. (2009) 45:286–97. doi: 10.1016/j.neuroimage.2008.11.038

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Caverzasi E, Hervey-Jumper SL, Jordan KM, Lobach IV, Li J, Panara V, et al. Identifying preoperative language tracts and predicting postoperative functional recovery using HARDI q-ball fiber tractography in patients with gliomas. J Neurosurg. (2016) 125:33–45. doi: 10.3171/2015.6.JNS142203

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sedrak M, Gorgulho A, Bari A, Behnke E, Frew A, Gevorkyan I, et al. Diffusion tensor imaging (DTI) and colored fractional anisotropy (FA) mapping of the subthalamic nucleus (STN) and the globus pallidus interna (GPi). Acta Neurochir. (2010) 152:2079–84. doi: 10.1007/s00701-010-0813-4

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Sedrak M, Gorgulho A, Frew A, Behnke E, DeSalles A, Pouratian N. Diffusion tensor imaging and colored fractional anisotropy mapping of the ventralis intermedius nucleus of the thalamus. Neurosurgery. (2011) 69:1124–9; discussion 1129-30. doi: 10.1227/NEU.0b013e3182296a42

CrossRef Full Text | Google Scholar

24. Jeurissen B, Leemans A, Tournier J.-D, Jones DK, Sijbers J. Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging. Hum Brain Mapp. (2013) 34:2747–66. doi: 10.1002/hbm.22099

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Maier-Hein KH, Neher PF, Houde, J.-C., Côté M.-A., et al. The challenge of mapping the human connectome based on diffusion tractography. Nat Commun. (2017) 8:1349. doi: 10.1038/s41467-017-01285-x

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sarwar T, Ramamohanarao K, Zalesky A. Mapping connectomes with diffusion MRI: deterministic or probabilistic tractography? Magn Reson Med. (2019) 81:1368–84. doi: 10.1002/mrm.27471

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Calabrese E, Badea A, Cofer G, Qi Y, Johnson GA. A diffusion MRI tractography connectome of the mouse brain and comparison with neuronal tracer data. Cereb Cortex. (2015) 25:4628–37. doi: 10.1093/cercor/bhv121

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Schilling KG, Petit L, Rheault F, Remedios S, Pierpaoli C, Anderson AW, et al. Brain connections derived from diffusion MRI tractography can be highly anatomically accurate-if we know where white matter pathways start, where they end, and where they do not go. Brain Struct Funct. (2020) 225:2387–402. doi: 10.1007/s00429-020-02129-z

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Chenot Q, Tzourio-Mazoyer N, Rheault F, Descoteaux M, Crivello F, Zago L, et al. A population-based atlas of the human pyramidal tract in 410 healthy participants. Brain Struct Funct. (2019) 224:599–612. doi: 10.1007/s00429-018-1798-7

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Avecillas-Chasin JM, Atik A, Mohammadi AM, Barnett GH. Laser thermal therapy in the management of high-grade gliomas. Int J Hyperth. (2020) 37:44–52. doi: 10.1080/02656736.2020.1767807

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Neher P, Stieltjes B, Wolf I, Meinzer H-P, Fritzsche KH. Analysis of tractography biases introduced by anisotropic voxels. ISMRM Volume. (2013) 21:1.

Google Scholar

32. Oouchi H, Yamada K, Sakai K, Kizu O, Kubota T, Ito H, et al. Diffusion anisotropy measurement of brain white matter is affected by voxel size: underestimation occurs in areas with crossing fibers. Am J Neuroradiol. (2007) 28:1102–6. doi: 10.3174/ajnr.A0488

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jones DK, Knösche TR, Turner R. White matter integrity, fiber count, and other fallacies: the do's and don'ts of diffusion MRI. Neuroimage. (2013) 73:239–54. doi: 10.1016/j.neuroimage.2012.06.081

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Girard G, Whittingstall K, Deriche R, Descoteaux M. Towards quantitative connectivity analysis: reducing tractography biases. Neuroimage. (2014) 98:266–78. doi: 10.1016/j.neuroimage.2014.04.074

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Catani M, Mesulam MM, Jakobsen E, Malik F, Martersteck A, Wieneke C, et al. A novel frontal pathway underlies verbal fluency in primary progressive aphasia. Brain. (2013) 136:2619–28. doi: 10.1093/brain/awt163

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Forkel SJ, Thiebaut de Schotten M, Dell'Acqua F, Kalra L, Murphy DGM, Williams SCR, et al. Anatomical predictors of aphasia recovery: a tractography study of bilateral perisylvian language networks. Brain. (2014) 137:2027–39. doi: 10.1093/brain/awu113

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Honey CJ, Sporns O, Cammoun L, Gigandet X, Thiran JP, Meuli R, et al. Predicting human resting-state functional connectivity from structural connectivity. Proc Natl Acad Sci USA. (2009) 106:2035–40. doi: 10.1073/pnas.0811168106

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Betzel RF, Byrge L, He Y, Goñi J, Zuo X.-N, Sporns O, et al. Changes in structural and functional connectivity among resting-state networks across the human lifespan. Neuroimage. (2014) 102:345–57. doi: 10.1016/j.neuroimage.2014.07.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: glioma, tractography, white matter, surgical resection, corticospinal tract

Citation: Sheng Z, Yu J, Chen Z, Sun Y, Bu X, Wang M, Sarica C, Hernesniemi J, Nelson BJ, Zemmar A and Avecillas-Chasin JM (2021) Constrained-Spherical Deconvolution Tractography in the Evaluation of the Corticospinal Tract in Glioma Surgery. Front. Surg. 8:646465. doi: 10.3389/fsurg.2021.646465

Received: 26 December 2020; Accepted: 24 May 2021;
Published: 29 July 2021.

Edited by:

Ben Jeurissen, University of Antwerp, Belgium

Reviewed by:

John A. Thompson, University of Colorado, United States
Alberto De Luca, University Medical Center Utrecht, Netherlands

Copyright © 2021 Sheng, Yu, Chen, Sun, Bu, Wang, Sarica, Hernesniemi, Nelson, Zemmar and Avecillas-Chasin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Josue M. Avecillas-Chasin, josueavecillas@hotmail.com; Ajmal Zemmar, ajmal.zemmar@gmail.com

These authors have contributed equally to this work and share senior authorship

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.