Qing Ji1, John O. Glass1, Elizabeth C. Duncan1, Amar Gajjar2, and Willburn E. Reddick1
1Diagnostic Imaging, St.Jude Children's Research Hospital, Memphis, TN, United States, 2Oncology, St.Jude Children's Research Hospital, Memphis, TN, United States
Synopsis
Global and regional cerebral
white matter changes induced by cranial
radiotherapy (CRT) were analyzed by comparing Pre-CRT and Post-CRT
diffusion tensor imaging (DTI) from 30 childhood medulloblastoma patients using
a structural connectivity network model and graph theory approaches. At the
global level, global network efficiency and character path were significantly
changed for the whole network. At the regional level, 17 of 82 network nodes
had significantly decreased local efficiencies, and 14 of those nodes also had
significantly decreased clustering coefficients. These findings suggest significant reduction in the microstructural
integrity immediately after CRT in this population.Introduction
Cranial radiotherapy
(CRT) is a radiotherapy technique critical to the treatment of primary or metastatic
brain tumors, which delivers a prescribed radiation dose uniformly across
the whole brain. Monitoring the brain tissue changes caused by CRT during and
after treatment is very important to understand the long-term effects induced
by CRT among the pediatric brain tumor survivors. Non-invasive diffusion tensor
imaging (DTI) not only offers an effective way to character tissue integrity by
fractional anisotropy (FA), but also makes it possible to study the brain as a structurally
connected network by fiber tractography. In this study, we report the results
of a graph theory analysis applied to a structural connectivity network model
to investigate the changes induced by CRT among pediatric brain tumor patients.
Method and Material
MR scans of 30 medulloblastoma patients (Age at exam 11.2+5.2 years) were used for this study. For each patient, two MR scans were performed, one was Pre-CRT and the other was Post-CRT. Each MR scan consists of anatomic 3D T1 weighted imaging and DTI (30 directions, 2 averages, b=700). For each scan, the anatomic imaging set was processed using Freesurfer [1] to obtain a re-sliced 3D isotropic (256 x 256 x 256, 1 mm resolution) T1 weighted image with 82 anatomic brain structures (7 sub-cortical and 34 cortical for each hemisphere). The eigenvectors and eigenvalues for the DTI data were first calculated using FSL FRMIB Toolbox [2] before being registered to the re-sliced T1 weighted image. To establish a reproducible network graph for each exam, probabilistic fiber tracking was then performed using FSL with 10,000 permutations using the 82 anatomic structures as graph nodes. If the connection between a pair of nodes was non-zero for 80% of all patient’s exams, the connection was considered valid. A total of 311 valid pairs of connections or edges were identified. The connection pathway between two nodes, which was the volume in image space that the connection fibers passed through, was extracted for each valid connection using a previously developed adaptation of the probabilistic fiber tracing technique [3]. The mean FA value of the connection pathway served as the quantitative measure for each edge. The graph network metrics were calculated at the global and nodal levels for each exam using Brain Connectome Toolbox [4] (http://site.google.com/a/brain-connectivity-toolbox.net/bct/Home). For each graph metric, the values of Pre-CRT and Post-CRT were compared across the patients using paired samples T-test.
Results
Figure 1 illustrated a typical pathway connecting two brain structures (left caudate and left rotralmiddlefrontal). Figure 2 demonstrated a typical 82 x 82 connection matrix with mean FA values of the connection pathway as the connection strength between nodes. In Table 1, the results of a comparison of Pre-CRT and Post-CRT for two global network metrics (global efficiency and character path) as well as the mean FA values were listed. We found that the two global network metrics were changed significantly in Post-CRT (p=0.007 for global efficiency and p=0.004 for character path) compared to Pre-CRT. Even though the mean FA value was also decreased in Post-CRT, it did not reach statistically significance (p=0.270). At the nodal level, we chose local efficiency and clustering coefficient as our network metrics. Among the 82 nodes, 17 were identified with significantly (p< 0.05) decreased local efficiencies and 14 of those nodes with significantly (p < 0.05) decreased clustering coefficients Post-CRT. These 17 nodes are listed in Table 2. It should be noted that the bilateral basal ganglia (caudate, pallidum and putamen) and thalamus were among those 17 structures.
Discussion and Conclusion
Changes in cerebral white matter due to CRT in pediatric brain tumor patients have been detected using a network model and graph metrics at both the global and regional level. Since the network model in this study combined tissue integrity measures with structural connection information, the network metrics were more sensitive for detecting the tissue response to CRT than the tissue properties alone.
Acknowledgements
This work was supported in part by Cancer Center Support (CORE) grant and R01CA90246 from the National Institutes of Health and by American Lebanese Syrian Associated Charities (ALSAC) References
[1] Fishl
B. et al, “Automatically
Parcellating the Human Cerebral Cortex”, Cerebral Cortex, 14,
pp. 11-22, 2004
[2]
Brehrens TE, et al. “Probabilistic
diffusion tractorgraphy with multiple fibre orientations: What can we gain?” Neuroimage, 34, pp. 144-55, 2007
[3] Ji Q.
et al, “Extraction of fiber Pathway using probabilistic fiber tracking”, OHBM, 2015
[4]
Rubinnov, M, et al, “Complex network
measures of brain connectivity: Uses and interpretations Original “, NeuroImage, 52, pp 1059-1069, 2010