- Open Access
Integration of BOLD-fMRI and DTI into radiation treatment planning for high-grade gliomas located near the primary motor cortexes and corticospinal tracts
Radiation Oncology volume 10, Article number: 64 (2015)
The main objective of this study was to evaluate the efficacy of integrating the blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) and diffusion tensor imaging (DTI) data into radiation treatment planning for high-grade gliomas located near the primary motor cortexes (PMCs) and corticospinal tracts (CSTs).
A total of 20 patients with high-grade gliomas adjacent to PMCs and CSTs between 2012 and 2014 were recruited. The bilateral PMCs and CSTs were located in the normal regions without any overlapping with target volume of the lesions. BOLD-fMRI, DTI and conventional MRI were performed on patients (Karnofsky performance score ≥ 70) before radical radiotherapy treatment. Four different imaging studies were conducted in each patient: a planning computed tomography (CT), an anatomical MRI, a DTI and a BOLD-fMRI. For each case, three treatment plans (3DCRT, IMRT and IMRT_PMC&CST) were developed by 3 different physicists using the Pinnacle planning system.
Our study has shown that there was no significant difference between the 3DCRT and IMRT plans in terms of dose homogeneity, but IMRT displayed better planning target volume (PTV) dose conformity. In addition, we have found that the Dmax and Dmean to the ipsilateral and contralateral PMC and CST regions were considerably decreased in IMRT_PMC&CST group (p < 0.001).
In conclusion, integration of BOLD-fMRI and DTI into radiation treatment planning is feasible and beneficial. With the assistance of the above-described techniques, the bilateral PMCs and CSTs adjacent to the target volume could be clearly marked as OARs and spared during treatment.
Gliomas, which contain oligodendroglia, astrocytic and ependymal lesions are the most common primary intracranial tumors. High-grade gliomas, which make up 35 to 45% of all newly diagnosed primary brain tumors worldwide, have a very poor prognosis . The three-dimensional conformal radiotherapy (3DCRT) has been considered as the standard therapy for patients with high-grade gliomas and intensity-modulated radiotherapy (IMRT) is becoming increasingly used to improve dose conformity and spare critical normal tissues. However, the risk of radiation-induced brain injury increases with the increase of radiation dose [2-4]. The strenuous endeavor has been made to diminish radiation complications.
With the assistance of conventional magnetic resonance imaging (MRI) and planning computed tomography (CT) data, many critical intracranial structures, such as lens, optic nerves and optic chiasm are well demarcated. However, it is difficult to accurately locate eloquent cortices and fiber connections in the white matter of the brain by routine neuroimaging. Excessive irradiation of eloquent cortices and white matter fiber tracts is unavoidable. Blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) and diffusion tensor imaging (DTI) have recently been used to identify the primary motor cortexes (PMCs) and corticospinal tracts (CSTs). These imaging techniques have been implemented in modern neuronavigation systems and used to guide the surgical removal of critically located intracranial lesions [5,6]. The purpose of our study was to evaluate whether the incorporation of BOLD-fMRI and DTI data into the 3D treatment planning process could spare the healthy brain and sensitive parts of the brain from high doses of radiation.
Ethical Approval was obtained from the General Hospital of Ningxia Medical University Review Board and written informed consent was obtained from patients. The study was conducted with strict adherence to the Declaration of Helsinki Principles.
A total of 20 patients with high-grade gliomas adjacent to PMCs and CSTs between May 2012 and February 2014 were recruited from the General Hospital of Ningxia Medical University, China. Eleven male patients and 9 female patients aged from 24 to 66 year-old were enrolled in this study.
The glioma tissues in our study included 14 astrocytomas (WHO Grade III) and 6 glioblastomas (WHO Grade IV). The bilateral PMCs and CSTs were located in the normal regions without any overlapping with target volume of the lesions. BOLD-fMRI, DTI and conventional MRI were performed on patients (Karnofsky performance score ≥ 70) before radical radiotherapy treatment.
Data acquisition and analysis
Four different imaging studies were conducted in each patient: a planning CT for radiosurgery treatment and target tracking during radiation therapy treatment delivering; an anatomical MRI to deliver a complete set of morphological MR data; a DTI to provide white matter tractography and a BOLD- fMRI to provide brain activation maps. Axial CT images (3-mm slice thickness) were taken by a wide-bore Siemens Somatom Sensation Open CT scanner (Siemens, Germany). MRI volumes were acquired using a Signal HDx 3.0 T MRI scanner (General Electric Company, USA).
BOLD fMRI data were obtained using fat-saturated single-shot gradient echo planar imaging (EPI) (TE = 35 ms, TR = 3,000 ms, acquisition matrix = 64 × 64 pixels, FOV = 240 mm × 240 mm2, flip angle = 90°, NEX = 1, 3 mm thickness). A block design paradigm (5 cycles, 30 sec on and 30 sec off) was utilized. Functional areas relevant to each treatment region were probed by Somatosensory tasks (finger tapping with audio cue). Following the acquisition of the functional data, gadolinium-enhanced high-resolution images were acquired (TR/TE = 450/14, flip angle = 90°, matrix = 256 × 256, FOV = 240 mm × 240 mm2, and slice thickness 3 mm skip 0 mm). After the images were taken, data were transferred to the Matlab workstation for analysis. The DTI data acquisition sequence was a spin echo-echo planar imaging (SE-EPI) sequence with TR = 10,000 ms, TE = 98.8 ms, acquisition matrix = 128 × 128 pixels; FOV = 240 mm × 240 mm2; slice thickness = 3.0 mm. Diffusion-weighted imaging with b factor of 1,000 mm2/s was taken along 25 noncollinear directions. The acquisition time of DTI sequence was 280 seconds. DTI data was analyzed on-line by the advantage workstation of the MR scanner (AW 4.4). For CST analysis, a seed region of interest (ROI) and a target ROI were placed on the posterior limb of the internal capsule and pons (anterior blue portion on the color map). Fiber tracking employed fractional anisotropy (FA) threshold of 0.2 and a tract angular change of 30°. The color-coded FA maps were merged with the anatomical MRI images. Regions of interest were drawn on the fused FA maps.
The fused fMRI activation maps and the white matter tracts overlaid on the anatomical MRI volume were exported as separate grayscale dicom images and loaded onto Pinnacle planning system software version 9.2 (Philips Medical Systems, Netherlands). The anatomical MRI images were registered with the CT volume for each patient. Figure 1 shows the anatomical MRI images were registered with the corresponding axial CT planes for a glioma case.
The target and organs at risk (OARs), i.e., optic nerves, optic chiasm and brain stem were precisely described using CT/anatomical MRI images. Both eyes were protected to avoid beam damage during treatment planning. The PMCs and the CSTs situated near the target were defined by a radiologist and a neurosurgeon, using the tractography images and the fused activation maps. Gross tumor volume (GTV) was described as the operative cavity with any remaining contrast-enhancing tissue on T1-weighted magnetic resonance imaging or as unresected enhancing tumor. The initial clinical target volume (CTV1) was defined as the T2 hyper intensity area (edema) with a 20 mm expansion. An initial planning target volume (PTV1) was created by adding a 30 mm expansion to the CTV1 to account for setup uncertainties. A second clinical target volume (CTV2) was defined as the contrast enhancement region in T1 with an additional 25 mm margin. A PTV2 was generated by adding a 30 mm expansion to the CTV2.
For each case, three treatment plans were developed by 3 different physicists using the Pinnacle planning system. The first physicist created conventional 3DCRT plans (3DCRT) and the PMCs and the CSTs situated near the target were not taken in account by the physicist. The target and standard morphological OARs were considered in this plan. The second physicist developed IMRT plans (IMRT) and the PMCs and the CSTs situated near the target were not considered by the physicist. The third physicist developed IMRT plans (IMRT_PMC&CST), and the PMCs and the CSTs situated near the target were considered. Figure 2 shows the axial isodose distribution of a patient with high-grade glioma.
The treatment plans met the requirement that at least 95% of the PTV receives the prescribed dose. Cumulative doses to the lenses, optic nerves, optic chiasm, and brainstem were limited to a maximum dose of 54 Gy for the last three structures and as low as practically achievable for the former. For conventional 3DCRT treatment, the prescribed dose was 50 Gy to the PTV1, immediately followed by 10 Gy to the PTV2, with a total cumulative dose of 50 Gy to the PTV1 and 60 Gy to the PTV 2 both at 2 Gy per fraction. For IMRT plans, the prescribed dose was 50 Gy to the PTV1 and 60 Gy to the PTV2, which were delivered concurrently over 30 daily fractions, with a fractional dose of 2 Gy to the PTV2.
Comparison criteria for the radiation treatment plans
The dose volume histograms (DVH) data were obtained from each patient. The dose coverage was analyzed according to the mean dose (Dmean), maximum dose (Dmax), conformity index (CI) and homogeneity index (HI). The CI was defined as follows : CI = VRI/PTV, where VRI represents the volume covered by the prescription dose. A CI value of 1.0 indicates that the volume of the prescription isodose surface is equal to that of the PTV. The HI was defined as follows : HI = (D2 − D98)/D50, where Dx% represents the dose delivered to x% of the PTV. Lower HI values indicate a more homogeneous target dose. OARs (e.g., brainstem, optic chiasm, optic nerves, and lenses) and PMCs and CSTs were compared based on the values of Dmax and Dmean.
The comparison of parameters between different plans was analyzed by the paired two-tailed Student t test. Differences were considered statistically significant at p < 0.05.
Comparison of target volume coverage and OAR sparing between 3DCRT and IMRT
Parameters related to dose coverage planning for 3DCRT and IMRT are presented in Tables 1 and 2. The results indicated that there was no significant difference between the 3DCRT and IMRT plans in terms of dose homogeneity, but IMRT displayed better PTV dose conformity. Regarding the comparison of PTV1 Dmax, PTV1 Dmean, PTV2 Dmax and PTV2 Dmean, there was no significant difference between the 3DCRT and IMRT plans. The dosimetric details of brainstem, optic chiasm, optic nerves, and lenses revealed no significant differences between the two plans and all of these organs were strictly maintained within the dose limitations. The Dmax and Dmean of PMCs and CSTs were observed in both 3DCRT and IMRT plans (Table 3); however no significant difference was found between the two plans.
Comparison of target volume coverage and OAR sparing between IMRT and IMRT_PMC&CST
Treatment plan parameters are shown in Tables 4 and 5. According to the data presented, both PTV1 and PTV2 (Dmax, Dmean, CI and HI) were analyzed and showed no significant differences between two groups. The Dmax and Dmean to these conventional OARs (e.g., brainstem, optic chiasm, optic nerves, and lenses) showed no significant differences between the two plans. The Dmax to the ipsilateral and contralateral PMC and CST regions was considerably decreased by 28.7%, 24.5%, 20.2% and 37.6%, respectively. The Dmean to the ipsilateral and contralateral PMC and CST regions was considerably decreased by 27.8%, 30.4%, 23.1% and 33.4%, respectively (Table 6).
Radiation therapy is commonly applied to the brain tumors due to its ability to control cell growth; however, radiation therapy can have detrimental effects on the central nervous system causing neurological complications. The response of cerebral tissue to radiation can lead to the deficits in neural functions [9,10]. The extent of neurologic deficit is associated with the location and size of radiation-induced brain injury [11,12]. Efforts dedicated to the precise division of brain lesions have been made to reduce the risk of neurological complications caused by the radiation therapy.
The DTI and BOLD-fMRI have recently been used to identify the white-matter pathways and functional structures of the brain. In our previous study, we proposed a clinically feasible protocol of integrating BOLD-fMRI and DTI to optimize the extent of resection involving the cortical motor areas and subcortical white matter tracts in patients with brain gliomas. Those information helped neurosurgeons resected the maximum amount of tumor while still preserving the most critical cortices of the brain, thus resulting in enhanced postoperative quality of life for patients . The incorporation of this information for radiosurgery planning has also been suggested. Liu et al. has reported a novel method to integrate the fMRI brain activation map with treatment planning for stereotactic radiosurgery (SRS). Direct irradiation of the eloquent cortices was avoided by multiple radiation arcs or static radiation beams in SRS planning, and the average dose reduction to the eloquent cortices was 32% . In addition, it has been reported that the risk of radiation-induced neuropathy was minimized by the integration of tractography of the brain white matter with DTI into radiation treatment planning of radiosurgery using Gamma Knife . Moreover, Pantelis et al. has demonstrated that critical structures of brain could be marked and spared with the aid of the integration of BOLD-fMRI and DTI into CyberKnife stereotactic radiosurgery .
In this study, BOLD-fMRI and DTI were used to localize the bilateral PMCs and CSTs and the information obtained from these two technologies were integrated into radiation treatment planning. The first part of our study (3DCRT versus IMRT) indicated that there was no significant reduction in the dose to bilateral PMCs and CSTs between the 3DCRT and IMRT plans. The critical structures adjacent to the target volume marked as OARs can be better spared during the IMRT planning process due to a steep dose gradient and a high conformity [17,18]. The second part of our study (IMRT versus IMRT_ PMC&CST), has shown that a significant reduction in the dose to bilateral PMCs and CSTs regions can be achieved without compromising the coverage of planning target volume and the limiting dose to these conventional OARs.
Sparing of the bilateral PMCs and CSTs does not represent any significant breakthrough in the treatment of brain tumors, but we have demonstrated that it is feasible to reduce the irradiation of critical structures adjacent to the target volume. The development of the most appropriate IMRT plan for the patient could be achieved by the identification of the important functional structures of the brain tissues proximal to the tumors. Sparing these vital functional structures is important to maintain quality of life, even in those patients with restricted life expectancy. The current study has shown that the DTI examination and MR Spectroscopy are valuable tools to differentiate the postoperative recurrent glioma from the radiation injury for patients with a glioma [19,20].
We are currently investigating paradigms for Broca’s and Wernicke’s areas (speech center), Broadmann-17 functional structures (visual center) and optic tracts. The application of the fMRI in low grade cases has also been validated.
In conclusion, integration of BOLD-fMRI and DTI into radiation treatment planning is feasible and beneficial. With the assistance of the above-described techniques, the structures adjacent to suspicious cancerous lesions could be clearly marked as OARs and spared during treatment. However, a wider investigation and the longer-term clinical follow up are required to further validate the effect of the integration of BOLD-fMRI and DTI on sparing normal tissues.
Blood oxygen level dependent functional magnetic resonance imaging
Diffusion tensor imaging
Three-dimensional conformal radiation treatment
Primary motor cortexes
Organs at risk
Intensity-modulated radiation therapy
- IMRT_ PMC&CST:
Intensity-modulated radiotherapy planning with PMC and CST information
World Health Organization
Magnetic resonance imaging
Spin echo-echo planar imaging
Region of interest
Gross tumor volume
Initial clinical target volume
Initial planning target volume
Second clinical target volume
Second planning target volume
Dose volume histograms
The maximum dose
Chang J, Narayana A. Functional MRI for radiotherapy of gliomas. Technol Cancer Res Treat. 2010;9(4):347–58.
Bleehen NM, Stenning SP. A Medical Research Council trial of two radiotherapy doses in the treatment of grades 3 and 4 astrocytoma. The Medical Research Council Brain Tumour Working Party. Br J Cancer. 1991;64(4):769–74.
Butler JM, Rapp SR, Shaw EG. Managing the cognitive effects of brain tumor radiation therapy. Curr Treat Options Oncol. 2006;7(6):517–23.
Roman DD, Sperduto PW. Neuropsychological effects of cranial radiation: current knowledge and future directions. Int J Radiat Oncol Biol Phys. 1995;31(4):983–98.
Kamada K, Todo T, Masutani Y, Aoki S, Ino K, Takano T, et al. Combined use of tractography-integrated functional neuronavigation and direct fiber simulation. J Neurosurg. 2005;102(4):664–72.
Pirotte B, Voordecker P, Neuroschl C, Baleriaux D, Wikler D, Metens T, et al. Combination of functional magnetic resonance imaging-guided neuronavigation and intraoperative cortical brain mapping improves targeting of motor cortex stimulation in neuropathic pain. Neurosurgery. 2008;62(6 Suppl 3):941–56.
Feuvret L, Noël G, Mazeron JJ, Bey P. Conformity index: a review. Int J Radiat Oncol Biol Phys. 2006;64(2):333–42.
Hodapp N. The ICRU report 83: prescribing, recording and reporting photon-beam intensity-modulated radiation therapy (IMRT). Strahlenther Onkol. 2012;188(1):97–9.
Keime-Guibert F, Napolitano M, Delattre JY. Neurological complications of radiotherapy and chemotherapy. J Neurol. 1998;245(11):695–708.
DeAngelis LM, Delattre J-Y, Posner JB. Radiation-induced dementia in patients cured of brain metastases. Neurology. 1989;39(6):789–96.
Flickinger JC, Kondziolka D, Lunsford LD, Kassam A, Phuong LK, Liscak R, et al. Development of a model to predict permanent symptomatic postradiosurgery. Int J Radiat Oncol Biol Phys. 2000;46(5):1143–8.
Voges J, Treuer H, Sturm V, Büchner C, Lehrke R, Kocher M, et al. Risk analysis of linear accelerator radiosurgery. Int J Radiat Oncol Biol Phys. 1996;36(5):1055–63.
Jia XX, Yu Y, Wang XD, Ma H, Zhang QH, Huang XY, et al. fMRI-driven DTT assessment of corticospinal tracts prior to cortex resection. Can J Neurol Sci. 2013;40(4):558–63.
Liu WC, Schulder M, Narra V, Kalnin AJ, Cathcart C, Jacobs A, et al. Functional magnetic resonance imaging aided radiation treatment planning. Med Phys. 2000;27(7):1563–72.
Maruyama K, Kamada K, Shin M, Itoh D, Aoki S, Masutani Y, et al. Integration of threedimensional corticospinal tractography into treatment planning for gamma knife surgery. J Neurosurg. 2005;102(4):673–7.
Pantelis E, Papadakis N, Verigos K, Stathochristopoulou I, Antypas C, Lekas L, et al. Integration of functional MRI, white matter tractography in stereotactic radiosurgery clinical practice. Int J Radiat Oncol Biol Phys. 2010;78(1):257–67.
Moretto F, Rampino M, Munoz F, Ruo Redda MG, Reali A, Balcet V, et al. Conventional 2D (2DRT) and 3D conformal radiotherapy (3DCRT) versus intensity-modulated radiotherapy (IMRT) for nasopharyngeal cancer treatment. Radiol Med. 2014;119(8):634–41.
Arbea L, Ramos LI, Martinez-Monge R, Moreno M, Aristu J. Intensity-modulated radiation therapy (IMRT) vs. 3D conformal radiotherapy (3DCRT) in locally advanced rectal cancer (LARC): dosimetric comparison and clinical implications. Radiat Oncol. 2010;5:17.
Xu JL, Li YL, Lian JM, Dou SW, Yan FS, Wu H, et al. Distinction between postoperative recurrent glioma and radiation injury using MR diffusion tensor imaging. Neuroradiology. 2010;52(12):1193–9.
Rabinov JD, Lee PL, Barker FG, Louis DN, Harsh GR, Cosgrove GR, et al. In vivo 3-T MR spectroscopy in the distinction of recurrent glioma versus radiation effects: initial experience. Radiology. 2002;225(3):871–9.
This work was supported by the Natural Science Foundation of China (Grant No. 81260373) and the Natural Science Foundation of Ningxia (Grant No. NZ11269).
The authors declare that they have no competing interests.
XDW, HCX, XYH and XXJ conceived the study and draft the manuscript. MLW and XDW performed the statistical analysis and drafted the manuscript. HM, YHG, XSX, YLG and HH participated in the analysis of the treatment of radiotherapy plans and draft this manuscript. YX edited this manuscript. All authors read and approved the final manuscript.
Minglei Wang and Xiaodong Wang contributed equally to this work.
About this article
Cite this article
Wang, M., Ma, H., Wang, X. et al. Integration of BOLD-fMRI and DTI into radiation treatment planning for high-grade gliomas located near the primary motor cortexes and corticospinal tracts. Radiat Oncol 10, 64 (2015) doi:10.1186/s13014-015-0364-1
- High-grade gliomas
- Blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI)
- Diffusion tensor imaging (DTI)
- Three-dimensional conformal radiation treatment (3DCRT)
- Intensity-modulated radiation therapy (IMRT)
- Radiation treatment planning