Modeling the relationship between fluorodeoxyglucose uptake and tumor radioresistance as a function of the tumor microenvironment.
As the most commonly used functional imaging modality in oncology, [sup.18]F-labeled fluorodeoxyglucose positron emission tomography (FDG-PET) has played a valuable role since first being synthesized in 1978 . Most malignant tumors exhibit an elevated glucose uptake and FDG-PET, as a close glucose chemical analogy, results in an image of glucose uptake in the patient, offering unique information for cancer detection, staging, target definition, and response monitoring [2, 3].
As a significant predictor of prognosis in radiation therapy (RT), FDG uptake seems to reflect an increased radioresistance; however, this is poorly understood. Many clinical outcome analyses have verified that high uptake of FDG in a tumor is correlated with increased local failure and shorter survival for many tumor sites, as summarized in several metaanalyses [4-7]. Therefore, FDG-avid regions in a tumor are recognized as a possible target for dose escalation to compensate for the radioresistance [8, 9]. Recently, utilizing a novel meta-analysis tool, we showed that FDG-avid head and neck tumors require about 20% more doses to equalize the local control rate with FDG nonavid tumors , although tumor size confounded that analysis to an unknown extent.
Enhanced glycolysis of tumor cells is certainly related to hypoxia, because hypoxic cells produce energy (in the form of ATP molecules) through glycolysis, without oxygen. However, it is also known that tumor cells can show increased glycolysis even in the presence of oxygen, compared to normal cells (the "Warberg effect") [11,12]. The oncolytic appetite for glycolysis is thought to be caused by a number of genetic or possibly epigenetic changes that drive malignancy [13,14].
Many studies have been carried out to correlate FDG uptake with various physiological parameters, such as hypoxia, proliferation, blood flow, histology, and differentiation, utilizing FDG-PET and immunohistochemical methods [15-20]. However, although several studies have shown the relationship between the FDG uptake and hypoxia or proliferation, the underlying mechanism of FDG uptake in a tumor is still unclear.
In this study, we do not try to resolve the detailed mechanism of FDG uptake. However, we do test various assumptions correlating FDG uptake (and presumably glucose consumption) with local cell microenvironmental conditions. The model is an attempt to incorporate known radiobiological effects that have been established as being important to understand radiotherapy treatment response, such as varying access to oxygen and glucose as well as the basic mathematical features of tumors, including variable growth fractions and cell loss factors. The key starting point of the model is that there is a limited amount of chemical resources for each tumor subvolume and that this level of resources is assumed to remain constant over a course of radiotherapy. We therefore used the model to determine the assumptions relating FDG uptake to the underlying cell compartments that best fits the observed correspondence between FDG uptake and reduced local control.
2.1. State-Driven Tumor Response Model. To explore the potential relationship between FDG-PET uptake and classical radiobiological mechanisms, a previously developed state-based tumor response model was used . In the mechanistic model, a tumor was assumed to be comprised of many small tumorlets of a PET-voxel-comparable size (4 x 4 x 4 [mm.sup.3]). Each tumorlet is comprised of three subpopulations of cells based on the level of proliferation, hypoxia, and cell loss, which is thought to be related to the available amount of oxygen and glucose, as shown by Kiran et al. . Figure 1 shows the three compartments of a tumorlet in a typical tumor microenvironment. Proliferation was assumed to take place only in the P-compartment, wherein cells have access to glucose and oxygen, whereas cell loss in the absence of treatment only takes place in the H-compartment, wherein cells have no access to either oxygen or glucose. The intermediate cells in the I-compartment, which have access to glucose but not oxygen, were assumed to be metabolically active but do not proliferate.
After RT begins, damaged cells (as deterministically calculated by linear-quadratic model) become doomed according to a compartment-specific radiosensitivity, leading to mitotic cell death when doomed cells are attempting to proliferate in the P-compartment. As doomed cells die, metabolically active but hypoxic cells in the 7-compartment receive oxygen and move into the P-compartment, causing reoxygenation. The model assumes a locally-constant blood supply over the course of RT. The transition of cells between compartments is determined dynamically, with the P-compartment always taking as many cells from the I-compartment as can be supported.
As a discrete-time simulation algorithm, the model keeps updating the number of cells in each compartment in a small time step. In each time step, a fraction of cells in the P-compartment proliferate, cell loss takes place in the H-compartment, the doomed cells in the P-compartment die probabilistically following mitosis, and the cells are re-compartmentalized based on the capacity of each compartment.
At the time of RT fraction, a fraction of viable cells in each compartment becomes doomed according to compartment-specific radiosensitivity value.
The model can be used to evaluate clinically important phenomena, including fraction size effects, the reoxygenation effect, repopulation effects, and tumor regression. Details of the mathematical model are available elsewhere .
2.2. Hypothetical FDG Uptake Patterns. We consider different relative uptakes of FDG within the three compartments. The subpopulations of tumor cells in the model were distinguished based on oxygen and glucose availability : only the P- and I-compartments were thought to be associated with FDG uptake given the model assumptions. Based on correlation studies between FDG and physiological factors, three different potential relationships between FDG uptake and cell subpopulations are considered here: (1) FDG uptake is proportional to the total number of metabolically viable cells (pattern I); (2) FDG uptake is associated mainly with proliferating cells with a reduced contribution from intermediate cells (pattern II); and (3) the FDG uptake is associated mainly with the intermediate cells with a reduced contribution from proliferating cells (pattern III), all shown in Figure 2. A reduced contribution from extremely hypoxic (starving and dying) cells in the H-compartment was additionally tested (uptake pattern IV).
2.3. Model Simulation. The model simulations were performed for a tumorlet treated with a standard RT regime of 2 Gy/fx (5 fx/week). The size of tumorlet was set to be 64 [mm.sup.3] (based on typical PET-voxel size) with 6.4 x [10.sup.7] cells, assuming the tumor cell density of [10.sup.6] [mm.sup.-3]. Only a small subset of cells is known to have stem-cell-like property and 1% of viable cells were assumed to be clonogenic cells in the model .
Based on previous work , relevant parameter values for head and neck squamous cell carcinoma (HNSCC) were used, including the radiosensitivity of the P-compartment ([[alpha].sub.p] = 0.382 [Gy.sup.-1], [alpha]/[beta] = 6.63 Gy). Hypoxic cells in the I- and H-compartments are considered to be only in the G0/G1-phase and the OER values for the I- and H-compartments were assumed to be 2 and 1.37, respectively, considering the lower OER of the G0/G1-phase and reduced repair capability of chronically hypoxic/highly stressed cells in the H-compartment .
The initial distribution of cells in each compartment is determined based on the presumed growth fraction (GF) and cell loss factor (CLF). As the GF increases, more cells are in proliferation and the P-compartment becomes larger. Given a higher CLF, more cells are starving and dying, resulting in a larger H-compartment. To consider a wide range of possible microenvironmental conditions, all possible ranges of initial cell distributions were simulated, based on various combinations of GFs and CLFs (1315 combinations). The parameter values for different compartments and the ranges of GF and CLF used in the simulation are summarized in Table 1.
2.4. Correlation between FDG Uptake and [TD.sub.50]. For each assumed FDG uptake pattern, the FDG uptake values were quantified for all possible ranges of initial distributions of cells in each compartment. Then, simulations were performed to determine the tumor dose for 50% control ([TD.sub.50]) for all the possible ranges of initial cell distributions. The [TD.sub.50] value was estimated from the Poisson formalism based on the total number of clonogenic cells in the tumorlet. Finally, we correlated the resulting FDG uptake values with the corresponding [TD.sub.50] values and determined the linear correlation coefficient for comparison purposes.
Estimated [TD.sub.50] values for all possible ranges of initial conditions are shown in Figure 3(a). [TD.sub.50] increased as the growth fraction decreased, due to a larger fraction of cells in the hypoxic compartments. [TD.sub.50] decreased as the cell loss factor increased but the dependency of [TD.sub.50] on the cell loss factor was much less significant compared to dependence on the growth fraction. The FDG uptake value was also quantified for all possible initial conditions based on the hypothesized FDG uptake patterns, as shown in Figure 3(b) for the uptake pattern I.
Hypothetical correlations between FDG uptake levels and corresponding [TD.sub.50] values are shown in Figure 4, for the four different hypothetical FDG uptake patterns. For the first assumed pattern (pattern I), where the FDG uptake was assumed to be proportional to the total number of viable cells, only a weak positive correlation was observed between FDG uptake and [TD.sub.50] value with a coefficient of determination ([R.sup.2]) of 0.38. When the uptake pattern was assumed to be mainly associated with cell proliferation with a minor contribution from intermediate cells (pattern II), a significant negative correlation existed between FDG uptake and [TD.sub.50] ([R.sup.2] = 0.76). For the hypothesis that metabolically-viable hypoxic cells in the I-compartment are avid for FDG uptake (pattern III), a strong positive correlation resulted ([R.sup.2] = 0.85). Inclusion of a reduced contribution from extremely hypoxic/starving cells in the H-compartment (pattern IV) yielded almost the same result as uptake pattern III, but with a slightly higher value of coefficient of determination ([R.sup.2] = 0.86), because the cells in H-compartment are either dying or reoxygenated into the I-compartment during the course of RT.
To better understand the potential relationship between FDG uptake and growth fraction (GF), the correlation between FDG uptake (for pattern IV) and [TD.sub.50] was evaluated for a fixed cell loss factor of 0.9, as shown in Figure 5. For the assumed uptake pattern, the FDG uptake was inversely correlated with the median GF in each bin.
Potential causes for the observed clinical correlation between FDG and radioresistance was explored using a mathematical model, in which classical radiobiological mechanisms were incorporated. Several different FDG uptake patterns were explored. Among the assumed FDG uptake patterns, when the metabolically viable hypoxic compartment (the I-compartment) was assumed to dominate the FDG uptake (uptake pattern III or IV), a significant positive correlation between FDG uptake and the required [TD.sub.50] was observed ([R.sup.2] = 0.85 or 0.86, resp.), implying that the increased cellular radioresistance due to chronic hypoxia may be the cause of the clinically observed increase in tumor radioresistance.
We tested for any potential role of starving and dying cells (in the H-compartment). When the total number of cells in each hypoxic compartment (I- or H-compartment) was correlated with the required [TD.sub.50] values, a strong positive correlation was observed for the I-compartment ([R.sup.2] = 0.86) but a weak negative correlation for the H-compartment ([R.sup.2] = 0.38), as shown in Figure 6.
This implies that the extremely hypoxic cells in the H-compartment are not the cause of an increased radioresistance. Mathematically, this is due to the lower OER value of the H-compartment used for the model simulation ([OER.sub.H] = 1.37), compared to the OER of the I-compartment ([OER.sub.I] = 2). To better understand this point, higher values of [OER.sub.H] ([OER.sub.H] = 2 or 3) were simulated (results not shown). Although the significance of the correlation was slightly reduced as the [OER.sub.H] increased; quantitatively similar relationships were observed. Even when the [OER.sub.H] was assumed to be as high as 3, which is certainly unreasonable radiobiologically, a strong correlation between the required [TD.sub.50] and the intermediate hypoxia in the I-compartment was still present ([R.sup.2] = 0.83).
In this work, it was assumed that there is no proliferation in the I-compartment for simplicity. Although it is known that some cells can proliferate even in hypoxic condition , the effect was evaluated to be not that significant due to the trade-off between increased proliferation and increased mitotic cell death, as shown in the previous work .
The ratio of [TD.sub.50] values between high versus low FDG uptake tumorlets, separated by median normalized FDG uptake, was estimated to be about 1.24 for the most likely uptake pattern, pattern III. This implies that high FDG uptake tumorlets require about 24% extra dose, which agrees with the clinically estimated range of the extra dose (~20%) for FDG-avid head and neck cancers . However, this estimate might be an over-estimate of the effect, since it assumes maximum heterogeneity of CLF and GF parameters, as shown in Table 1.
The results suggest that the total number of metabolically viable hypoxic cells (in I-compartment of the model) is a deterministic factor in tumor response and this subpopulation might be associated with the voxel-by-voxel correlation studies performed by Pugachev et al.  and Rajendran et al. , in which the microregional relationship between FDG and hypoxia was observed. Also the result supports the study of Wouters and Brown, in which the importance of the cells at intermediate oxygen level was emphasized .
In this work, the potential relationship between FDG-PET uptake and classical radiobiological mechanisms was explored using a mathematical framework based on a law of conservation of chemical resources. Several different FDG uptake patterns were hypothesized and the estimated FDG uptake values were correlated with tumor dose for 50% control ([TD.sub.50]) using the state-driven tumor response model. The model could generate a correlation between FDG uptake and an increase in tumor radioresistance when it was assumed that cells receiving glucose but not oxygen take up more glucose than cells that are well-oxygenated. Establishing this correspondence further, and understanding its limitations, will require appropriate FDG-PET clinical datasets.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This research was supported by Enid A. Haupt Medical Physics Endowed Chair fund at Memorial Sloan-Kettering Cancer Center, as well as the Nuclear Science and Engineering Institute of the University of Missouri and the Department of Radiation Oncology of the Washington University School of Medicine.
 T Ido, C.-N. Wan, and V. Casella, "Labeled 2-deoxy-D-glucose analogs. [sup.18]F-labeled 2-deoxy-2-fluoro-D-glucose, 2-deoxy-2-fluoro-D-mannose and [sup.14]C-deoxy-2-fluoro-D-glucose," Journal of Labelled Compounds and Radiopharmaceuticals, vol. 14, no. 2, pp. 175-183, 1978.
 M. H. M. Janssen, H. J. W. L. Aerts, M. C. Ollers et al., "Tumor delineation based on time-activity curve differences assessed with dynamic fluorodeoxyglucose positron emission tomography-computed tomography in rectal cancer patients," International Journal of Radiation Oncology Biology Physics, vol. 73, no. 2, pp. 456-465, 2009.
 S. Ben-Haim and P. Ell, "18F-FDG PET and PET/CT in the evaluation of cancer treatment response," Journal of Nuclear Medicine, vol. 50, no. 1, pp. 88-99, 2009.
 P. Xie, M. Li, H. Zhao, X. Sun, Z. Fu, and J. Yu, "F-FDG PET or PET-CT to evaluate prognosis for head and neck cancer: a metaanalysis," Journal of Cancer Research and Clinical Oncology, vol. 137, no. 7, pp. 1085-1093, 2011.
 H. Ghooshkhanei, G. Treglia, G. Sabouri, R. Davoodi, and R. Sadeghi, "Risk stratification and prognosis determination using 18F-FDG PET imaging in endometrial cancer patients: a systematic review and meta-analysis," Gynecologic Oncology, vol. 132, no. 3, pp. 669-676, 2014.
 K. Pak, G. J. Cheon, H. Y. Nam et al., "Prognostic value of metabolic tumor volume and total lesion glycolysis in head and neck cancer: a systematic review and meta-analysis," Journal of Nuclear Medicine, vol. 55, no. 6, pp. 884-890, 2014.
 Z. Wang, J.-Q. Chen, J.-L. Liu, X.-G. Qin, and Y. Huang, "FDG-PET in diagnosis, staging and prognosis of pancreatic carcinoma: a meta-analysis," World Journal of Gastroenterology, vol. 19, no. 29, pp. 4808-4817, 2013.
 M. Feng, F.-M. Kong, M. Gross, S. Fernando, J. A. Hayman, and R. K. Ten Haken, "Using fluorodeoxyglucose positron emission tomography to assess tumor volume during radiotherapy for non-small-cell lung cancer and its potential impact on adaptive dose escalation and normal tissue sparing," International Journal of Radiation Oncology Biology Physics, vol. 73, no. 4, pp. 1228-1234, 2009.
 D. de Ruysscher, S. Wanders, A. Minken et al., "Effects of radiotherapy planning with a dedicated combined PET-CT-simulator of patients with non-small cell lung cancer on dose limiting normal tissues and radiation dose-escalation: a planning study," Radiotherapy and Oncology, vol. 77, no. 1, pp. 5-10, 2005.
 J. Jeong, J. S. Setton, N. Y. Lee, J. H. Oh, and J. O. Deasy, "Estimate of the impact of FDG-avidity on the dose required for head and neck radiotherapy local control," Radiotherapy and Oncology, vol. 111, no. 3, pp. 340-347, 2014.
 O. Warburg, "On the origin of cancer cells," Science, vol. 123, no. 3191, pp. 309-314, 1956.
 R. A. Gatenby and R. J. Gillies, "Why do cancers have high aerobic glycolysis?" Nature Reviews Cancer, vol. 4, no. 11, pp. 891-899, 2004.
 J.-W. Kim and C. V. Dang, "Cancer's molecular sweet tooth and the warburg effect," Cancer Research, vol. 66, no. 18, pp. 8927-8930, 2006.
 M. G. V. Heiden, L. C. Cantley, and C. B. Thompson, "Understanding the warburg effect: the metabolic requirements of cell proliferation," Science, vol. 324, no. 5930, pp. 1029-1033, 2009.
 J. G. Rajendran, D. A. Mankoff, F. O'Sullivan et al., "Hypoxia and glucose metabolism in malignant tumors: evaluation by [18F]Fluoromisonidazole and [18F]fluorodeoxyglucosepositron emission tomography imaging," Clinical Cancer Research, vol. 10, no. 7, pp. 2245-2252, 2004.
 A. Pugachev, S. Ruan, S. Carlin et al., "Dependence of FDG uptake on tumor microenvironment," International Journal of Radiation Oncology Biology Physics, vol. 62, no. 2, pp. 545-553, 2005.
 M. Zimny, B. Gagel, E. DiMartino et al., "FDG--a marker of tumour hypoxia? A comparison with [18F] fluoromisonidazole and p[O.sub.2]-polarography in metastatic head and neck cancer," European Journal of Nuclear Medicine and Molecular Imaging, vol. 33, no. 12, pp. 1426-1431, 2006.
 I. Buchmann, U. Haberkorn, I. Schmidtmann et al., "Influence of cell proportions and proliferation rates on FDG uptake in squamous-cell esophageal carcinoma: a PET study," Cancer Biotherapy and Radiopharmaceuticals, vol. 23, no. 2, pp. 172-180, 2008.
 R. A. Dierckx and C. van de Wiele, "FDG uptake, a surrogate of tumour hypoxia?" European Journal of Nuclear Medicine and Molecular Imaging, vol. 35, no. 8, pp. 1544-1549, 2008.
 H. Vesselle, A. Salskov, E. Turcotte et al, "Relationship between non-small cell lung cancer FDG uptake at PET, tumor histology, and Ki-67 proliferation index," Journal of Thoracic Oncology, vol. 3, no. 9, pp. 971-978, 2008.
 J. Jeong, K. I. Shoghi, and J. O. Deasy, "Modelling the interplay between hypoxia and proliferation in radiotherapy tumour response," Physics in Medicine and Biology, vol. 58, no. 14, pp. 4897-4919, 2013.
 K. L. Kiran, D. Jayachandran, and S. Lakshminarayanan, "Mathematical modelling of avascular tumour growth based on diffusion of nutrients and its validation," The Canadian Journal of Chemical Engineering, vol. 87, no. 5, pp. 732-740, 2009.
 C. Hemmings, "The elaboration of a critical framework for understanding cancer: the cancer stem cell hypothesis," Pathology, vol. 42, no. 2, pp. 105-112, 2010.
 N. Chan, M. Koritzinsky, H. Zhao et al., "Chronic hypoxia decreases synthesis of homologous recombination proteins to offset chemoresistance and radioresistance," Cancer Research, vol. 68, no. 2, pp. 605-614, 2008.
 K. I. E. M. Wijffels, H. A. M. Marres, J. P W. Peters, P F. J. W. Rijken, A. J. van der Kogel, and J. H. A. M. Kaanders, "Tumour cell proliferation under hypoxic conditions in human head and neck squamous cell carcinomas," Oral Oncology, vol. 44, no. 4, pp. 335-344, 2008.
 B. G. Wouters and J. M. Brown, "Cells at intermediate oxygen levels can be more important than the "hypoxic fraction" in determining tumor response to fractionated radiotherapy," Radiation Research, vol. 147, no. 5, pp. 541-550, 1997
Jeho Jeong and Joseph O. Deasy
The Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA
Correspondence should be addressed to Joseph O. Deasy; email@example.com
Received 4 July 2014; Accepted 25 August 2014; Published 8 September 2014
Academic Editor: Iuliana Toma-Dasu
Table 1: Parameter values for different compartments and the ranges of GF and CLF used in the model simulation. Compartment P I H Oxygen 1 2 1.37 enhancement ratio (OER) Glucose uptake Varies depending on the hypothetical FDG level per cell uptake pattern (a) Proliferation 50-100% (b) -- -- fraction Cell loss Mitotic cell -- Cell loss due mechanism and death due to to starvation rate irradiation (cell loss (survival rate halftime: 2 of each progeny days) after mitosis: 0.3) Growth fraction 0.01 to 1/(2 + (GF) range CLF) [0.01 applied step] Cell loss 0.03 to 0.99 factor (CLF) [0.03 step] range applied (a) The relative uptake ratios in P:I: H are 1:1:0, 3:1:0, 1:3:0, and 2:5:2 for patterns I, II, III, and IV, respectively, as shown in Figure 2. (b) Depends on the fullness of the compartment: 50% of proliferation was assumed when the compartment is full, and 100% of proliferation was assumed when the number of cells is less than half of the capacity of the P-compartment.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Jeong, Jeho; Deasy, Joseph O.|
|Publication:||Computational and Mathematical Methods in Medicine|
|Date:||Jan 1, 2014|
|Previous Article:||Screening for cervical cancer using automated analysis of PAP-smears.|
|Next Article:||Dosimetric algorithm to reproduce isodose curves obtained from a LINAC.|