Abstract
Intratumor heterogeneity in biologic properties and in relationships between various phenotypes may present a challenge for biologically targeted therapies. Understanding the relationships between different phenotypes in individual tumor types could help inform treatment selection. The goal of this study was to characterize spatial correlations of glucose metabolism, proliferation, and hypoxia in 2 histologic types of tumors. Methods: Twenty canine veterinary patients with spontaneously occurring sinonasal tumors (13 carcinomas and 7 sarcomas) were imaged with 18F-FDG, 18F-labeled 3′-deoxy-3′-fluorothymidine (18F-FLT), and 61Cu-labeled diacetyl-bis(N4-methylthiosemicarbazone) (61Cu-ATSM) PET/CT on 3 consecutive days. Precise positioning and immobilization techniques coupled with anesthesia enabled motionless scans with repeatable positioning. Standardized uptake values (SUVs) of gross sarcoma and carcinoma volumes were compared by use of Mann–Whitney U tests. Patient images were rigidly registered together, and intratumor tracer uptake distributions were compared. Voxel-based Spearman correlation coefficients were used to quantify intertracer correlations, and the correlation coefficients of sarcomas and carcinomas were compared. The relative overlap of the highest uptake volumes of the 3 tracers was quantified, and the values were compared for sarcomas and carcinomas. Results: Large degrees of heterogeneity in SUV measures and phenotype correlations were observed. Carcinoma and sarcoma tumors differed significantly in SUV measures, with carcinoma tumors having significantly higher 18F-FDG maximum SUVs than sarcoma tumors (11.1 vs. 5.0; P = 0.01) as well as higher 61Cu-ATSM mean SUVs (2.6 vs. 1.2; P = 0.02). Carcinomas had significantly higher population-averaged Spearman correlation coefficients than sarcomas in comparisons of 18F-FDG and 18F-FLT (0.80 vs. 0.61; P = 0.02), 18F-FLT and 61Cu-ATSM (0.83 vs. 0.38; P < 0.0001), and 18F-FDG and 61Cu-ATSM (0.82 vs. 0.69; P = 0.04). Additionally, the highest uptake volumes of the 3 tracers had significantly greater overlap in carcinomas than in sarcomas. Conclusion: The relationships of glucose metabolism, proliferation, and hypoxia were heterogeneous across different tumors, with carcinomas tending to have high correlations and sarcomas having low correlations. Consequently, canine carcinoma tumors are robust targets for therapies that target a single biologic property, whereas sarcoma tumors may not be well suited for such therapies. Histology-specific PET correlations have far-reaching implications for the robustness of biologic target definition.
Two types of biologic heterogeneity are often associated with malignant tumors: heterogeneity between tumors and heterogeneity within tumors. Biologic heterogeneity between tumors occurs when various degrees of biologic properties, such as oxygenation status or metastatic potential, are observed across different tumors; this heterogeneity causes variable clinical responses to standard-of-care therapies. Heterogeneity within tumors—or intratumor heterogeneity—occurs when a single tumor has regional variations in biologic—and even genetic—properties (1). Intratumor heterogeneity presents a challenge for advanced therapies in which specific biologic processes are targeted: certain tumor regions might express the biologic target, but other tumor regions might not. With increasing interest in therapies that target specific biologic mechanisms, measuring intratumor heterogeneity might help inform treatment selection and could influence staging and treatment response assessment. Functional imaging, such as PET, is a noninvasive method for determining how biologic properties are spatially distributed throughout a tumor.
Although measuring spatial patterns of a single phenotype within a tumor may help describe certain biologic aspects of the tumor, understanding the spatial relationships between different biologic properties provides greater insight into the underlying biologic status of the tumor. For example, Busk et al. found that the spatial concordance between tumor glucose metabolism and hypoxia varied with different histologic types of tumors, indicating that different histologic types of tumors have various degrees of dependence on aerobic glycolysis (2). Lohith et al. also found significant differences in the spatial relationships between glucose metabolism and hypoxia in lung adenocarcinoma and lung squamous cell carcinoma tumors, as assessed by PET (3). Understanding the spatial relationships between tumor phenotypes may help both researchers and clinicians in tasks such as biomarker development, therapy selection, and interpretation of responses to therapies.
Biologically conformal radiotherapy, or dose painting, is another area in which spatial concordance or discordance between different biologic properties within a tumor can have a large impact. In dose painting, a new paradigm in radiation therapy, a nonuniform dose is tailored to match the spatial distribution of the biologic properties of a tumor (4). However, it is uncertain which biologic property should be targeted. Tumors in which different biologic targets (e.g., cellular hypoxia and glucose metabolism) have similar spatial distributions would be robust targets for dose painting: targeting a single phenotype would result in a similar treatment plan to targeting other phenotypes. However, in tumors in which different biologic targets cluster in different regions of the tumor, dose painting treatment plans would be sensitive to the choice of biologic target. Measuring spatial distributions of different biologic targets could help identify histologic types of tumors that are most likely to respond to dose painting treatments.
The purposes of this study were to measure and compare the spatial distributions of 3 biologic properties in tumors—glucose metabolism, cellular proliferation, and hypoxia—in 2 histologic types of tumors. Surrogates of these properties were measured with PET radiotracers 18F-FDG, 18F-labeled 3′-deoxy-3′-fluorothymidine (18F-FLT), and 61Cu(II)-labeled diacetyl-bis(N4-methylthiosemicarbazone) (61Cu-ATSM), respectively. Canine veterinary patients with spontaneously occurring carcinoma and sarcoma sinonasal tumors were imaged. These patients were ideal subjects for this study because we were able to eliminate motion during imaging, ensure repeatable positioning, and accurately register images. Furthermore, dogs are excellent models for understanding human cancers (5). We compared the 18F-FDG, 18F-FLT, and 61Cu-ATSM PET uptake patterns of canine sarcoma and carcinoma tumors; our experiments included a comparison of uptake magnitudes and an analysis of the spatial relationships of the 3 tracer distributions.
MATERIALS AND METHODS
Patient Population
Twenty canine patients of different breeds were included in the study. Companion animal (pet) dogs were referred to the University of Wisconsin Veterinary Medical Teaching Hospital. The research protocol was approved by the Animal Care and Use Committee of the University of Wisconsin, and all canine owners signed a written informed consent form. Patients were diagnosed with CT, and biopsies for histopathologic evaluation were obtained. All dogs had nasal or paranasal sinus tumors of stages 1–3 (mean tumor volume, 56 cm3), with no evidence of distant metastases or intracranial extension. Histopathologic evaluation revealed that 12 dogs had adenocarcinomas, 6 had chondrosarcomas, 1 had a squamous cell carcinoma, and 1 had an osteosarcoma. A veterinarian monitored all imaging sessions.
Tracer Production
61Cu (half-life, 204.5 min) was produced via the 60Ni(d,n)61Cu reaction at the University of Wisconsin cyclotron. The 61Cu-ATSM tracer was synthesized via the process outlined by Avila-Rodriguez (6). Analytic high-performance liquid chromatography of 61Cu-ATSM revealed a radiochemical purity of greater than 95% and typical specific activities of 3.7 × 1010–11.1 × 1010 Bq/μmol (1–3 Ci/μmol). The 18F-FLT tracer was also produced at the University of Wisconsin cyclotron, and the 18F-FDG tracer was provided by IBA Molecular.
Imaging Procedure
Each canine patient underwent 18F-FDG, 18F-FLT, and 61Cu-ATSM PET/CT imaging on 3 consecutive days as part of a phase I/II clinical trial investigating dose painting. Scans were performed approximately 24 h apart.
The protocol for PET/CT imaging was as follows. Dogs fasted for at least 12 h before each scan. They were injected with 150–370 MBq (4–10 mCi) of tracer, depending on their size. After injection, the dogs were kept in a kennel to limit physical activity. The dogs were anesthetized during scans with an initial propofol bolus injection and then isoflurane plus 100% oxygen. The time between tracer injection and the start of 18F-FDG PET scans was 60 min; for 61Cu-ATSM scans, the interval was 3 h. 18F-FLT scans were acquired dynamically. For reproducible positioning between PET/CT scans, the dogs’ maxillae were positioned into custom dental molds fixed to the scanner couch, and their bodies were immobilized by vacuum mattresses (7).
All scans were performed on a Discovery VCT (GE Healthcare) PET/CT scanner. 18F-FDG and 61Cu-ATSM scans were 20-min 3-dimensional static acquisitions over a single 15-cm bed position. 18F-FLT scans were 90-min 3-dimensional dynamic acquisitions over a single 15-cm bed position. Emission data were attenuation-corrected and reconstructed with ordered-subset expectation maximization (2 iterations, 35 subsets, and 3-mm postfiltering). The image grid was 256 × 256 × 47, with voxel sizes of 2.0 × 2.0 × 3.3 mm. Activity measurements were converted to standardized uptake values (SUVs) for analysis. For 18F-FLT scans, SUVs were calculated by averaging the frames from 60 to 90 min.
Data Analysis
The magnitudes of 18F-FDG, 18F-FLT, and 61Cu-ATSM uptake in sarcoma (n = 7) and carcinoma (n = 13) tumors were compared. Veterinarians contoured gross tumor volumes (GTVs) on the basis of CT images. The maximum voxel SUV of the GTV (SUVmax), the mean SUV of a 1-cm3 sphere centered at the SUVmax (SUVpeak), and the mean SUV of the GTV (SUVmean) were calculated for each canine patient. The SUVmax, SUVpeak, and SUVmean of sarcomas and carcinomas were then compared, and 2-sided Mann–Whitney U tests were used to determine significant differences. Tests of significance were robust to outliers in the distribution. In addition, F tests were used to compare the variances of SUV measures in sarcomas and carcinomas.
For each patient, all CT images were cropped to include only the volume of interest and then rigidly registered in Amira (Visage Imaging Inc.) to a single reference CT image by use of mutual information; the resulting transformations were applied to their respective PET data. This process allowed for a voxel-based comparison of 18F-FDG, 18F-FLT, and 61Cu-ATSM uptake distributions in each tumor. Voxel-based Spearman rank correlation coefficients were calculated for each tracer combination: 18F-FLT–61Cu-ATSM, 18F-FDG–61Cu-ATSM, and 18F-FDG–18F-FLT. Spearman rank correlation coefficients were used instead of Pearson correlation coefficients to avoid assumptions of linearity and normality. Correlation coefficients were averaged for sarcomas and carcinomas with the Fisher transformation, and 2-sided t tests were used to determine significant differences between the 2 tumor histology z score distributions.
In addition to the voxel-based analysis, tracer uptake patterns in sarcoma and carcinoma tumors were compared by analyzing the degrees of overlap of tracer-avid uptake regions. Tracer-avid volumes were created for each tracer by applying a threshold to each tracer’s SUVmax within the GTV, resulting in 3 tracer-avid volumes per tumor. From these 3 volumes we calculated the percentage of overlap relative to the total threshold volume (i.e., the intersect relative to the union). Because the best threshold for segmenting tumor volumes is unknown, we applied several different thresholds, from 10% to 90% of the SUVmax in 10% increments. At each threshold level, the percentage of overlap was quantified. The overlap percentages were then averaged for both histologic types of tumors, and t tests were used to determine whether the degrees of overlap in the 2 histologic types of tumors were significantly different. To better visualize the histology-averaged tracer overlap, we created Venn diagrams (8) for each threshold level.
RESULTS
Magnitude of Tracer Uptake
Table 1 shows SUV measures averaged for sarcomas and carcinomas, with P values resulting from Mann–Whitney U tests. Carcinomas had significantly higher 18F-FDG SUVmax than sarcomas (11.1 vs. 5.0; P = 0.01) as well as higher 61Cu-ATSM SUVmean (2.6 vs. 1.2; P = 0.02). Differences in 18F-FDG SUVmean and SUVpeak approached significance (P = 0.057). For all SUV measures, carcinomas had higher average levels of uptake than sarcomas. Box plots of SUVmean and SUVmax measurements are shown in Figure 1. It is evident from Figure 1 that the interpatient range, or variance, in SUV measures was much larger for carcinomas than for sarcomas. This finding was confirmed by F tests, which showed the variances in the 2 histologic types of tumors to be significantly different (P < 0.01) in all cases except the 61Cu-ATSM SUVmax and SUVmean. The average tumor volumes in the 2 histologic types of tumors were not significantly different.
Spatial Correlations
Figure 2 shows voxel-based SUV scatterplots for 3 canine patients with carcinomas and 3 canine patients with sarcomas. For 4 of the 7 sarcoma tumors, the 18F-FLT–61Cu-ATSM scatterplots appeared to be bifurcated, with one arm of low hypoxia and various degrees of proliferation and another arm with low proliferation and various degrees of hypoxia. For the other 3 sarcoma tumors, the scatterplots were fan-shaped. On the other hand, for carcinoma tumors, correlations between phenotypes were generally high. Box plots of voxel-based Spearman correlation coefficients, separated by histologic types of tumors, are shown in Figure 3. Fisher-averaged correlation coefficients and t test results are shown in Table 2. Carcinomas had significantly higher intertracer correlations than sarcomas for all tracer comparisons. The greatest differences between the 2 histologic types of tumors were observed for the 18F-FLT–61Cu-ATSM correlations: the average correlation coefficients were 0.38 for sarcomas and 0.83 for carcinomas (P < 10−4).
Overlapping Volumes
Examples of tracer-avid volumes at a 70% threshold are shown for 12 tumors in Figure 4. Tracer-avid regions generally overlapped in carcinoma tumors but were spatially separated in sarcoma tumors. Figure 5 shows Venn diagram representations of intertracer overlap, averaged for carcinomas and sarcomas, at 9 SUV thresholds. For low thresholds (≤30%) and high thresholds (90%), no significant differences between the 2 histologic types of tumors were observed. However, for thresholds between 40% and 80%, significant differences between sarcomas and carcinomas were observed; t test results are shown in Table 3. No significant differences between sarcomas and carcinomas were observed for 18F-FDG–18F-FLT overlap at any threshold. For 61Cu-ATSM–18F-FDG and 61Cu-ATSM–18F-FLT overlap, differences between sarcomas and carcinomas were significant at thresholds of 40%–80%, with carcinomas having greater average overlap. The regions in which all 3 tracers overlapped were also significantly larger in carcinomas at thresholds of 40%–70%.
DISCUSSION
In comparing 18F-FDG, 18F-FLT, and 61Cu-ATSM uptake patterns in canine nasal tumors, we found that sarcoma tumors had significantly lower levels of 18F-FDG and 61Cu-ATSM uptake than carcinoma tumors. Two other studies—1 in dogs (9) and 1 in human lung cancer patients (10)—compared 18F-FDG uptake in sarcoma and carcinoma tumors and found similar results. Lower levels of 18F-FDG and 61Cu-ATSM uptake may be indicative of a less aggressive phenotype in sarcoma tumors and may partly explain why patients with sarcomas respond better to treatment (11). However, low levels of uptake may also be caused by a lower tumor cell density or poor vascularity, which are characteristic of chondrosarcoma tumors (12).
The present study was unique in its ability to establish accurate voxelwise phenotype correlations. As registration inaccuracies can significantly degrade voxelwise correlations between images (13), canine nasal tumors were chosen for the present study because the nondeformable bony anatomy surrounding the nasal cavity enables accurate image registration, mitigating such errors. Motion blurring can also make image comparison difficult by smoothing out heterogeneities and causing images to appear similar (14). However, because canine patients experienced no motion during scans as a result of anesthesia, differences between images were apparent, as was notably observed in patients with sarcomas. In addition, tumors were large enough to provide significant intratumor heterogeneity. Because scans were obtained on consecutive days, uncertainties associated with tumor deformation and cellular redistribution were minimal. These precisely controlled conditions allowed for an accurate assessment of phenotype correlations, which might not be reproducible in human studies.
Spatial correlations of 18F-FDG, 18F-FLT, and 61Cu-ATSM PET were heterogeneous across canine patients, with Spearman correlation coefficients ranging from 0.14 to 0.98. Carcinoma tumors had significantly higher correlations and levels of colocalization of tracers than sarcoma tumors. The greatest discordance between the 2 histologic types of tumors was noted for the 18F-FLT–61Cu-ATSM correlations. In 4 of the 7 patients with sarcomas, there were distinctive bifurcations in the 18F-FLT–61Cu-ATSM scatterplots (3 of these patients also had minor bifurcations in the 18F-FDG–61Cu-ATSM scatterplots), which have not been reported in the literature. The 2 branches, when mapped back into 3-dimensional image space, resulted in 2 spatially contiguous tumor regions. This finding suggests that significant phenotypic heterogeneity exists in sarcoma tumors, with distinct regions of proliferation, distinct regions of hypoxia, and various levels of metabolism.
The biologic underpinnings of the observed uptake patterns are unclear but may be related to different cellular responses to hypoxia. Assuming that high levels of 18F-FDG, 18F-FLT, and 61Cu-ATSM uptake indicate the colocalization of metabolism, proliferation, and hypoxia, the high intertracer correlations observed in carcinomas may be explained by increased cellular hypoxia–inducible factor 1 activity, which can activate metabolic enzymes and transporters and drive proliferation (15). This reasoning may be supported by a study by Kaira et al., who reported that human lung carcinoma patients had significantly higher levels of hypoxia–inducible factor 1α expression than lung sarcoma patients (10). Conversely, the low intertracer correlations observed in sarcomas may be due to a preference of sarcoma cells to undergo cell cycle arrest under hypoxic conditions, thereby lowering proliferation rates (16). However, because the SUV of a voxel reflects only an average of the many cellular processes occurring within its macroscopic volume, we cannot explain the trends with certainty. Interestingly, although sarcoma tumors had significantly more phenotypic heterogeneity than carcinoma tumors, canine nasal sarcoma tumors were reported to respond better to treatment than carcinoma tumors (11). In contrast, other studies showed that intratumor heterogeneity, although defined and measured differently, was associated with a poorer outcome (17).
Treatments designed to specifically target certain tumor phenotypes (e.g., dose painting) may result in different outcomes for carcinoma and sarcoma tumors. In carcinomas, in which proliferative, hypoxic, and metabolic phenotypes were colocalized, targeting a single phenotype would likely simultaneously target the other 2 phenotypes. Alternatively, in sarcomas, each tracer accumulated in different regions of the tumor, and targeting a single phenotype would likely miss tumor regions with high levels of expression of the other 2 phenotypes. It is also possible that in different histologic types of tumors, the relationships between various biologic properties and treatment resistance may be different. The latter was observed in a study involving a smaller subset of our patient population, in which we observed histology dependence when associating distributions of different pretreatment tracers with follow-up 18F-FDG distributions (18).
The present study is the first study to report comparisons of 3 tracer distributions in 2 histologic types of tumors. Results from previous studies comparing 2 tracers in 2 or more histologic types of tumors were in agreement with our results—namely, that spatial correlations of different biologic properties are histology-specific (2,3,19). Tumor types with poor correlations between metabolic activity and hypoxia include colon adenocarcinomas (20) and Lewis lung carcinomas (21) in rodents and lung squamous cell carcinomas (3) and sarcomas (19) in humans. On the other hand, high correlations have been found in patients with head and neck cancers (19) and lung adenocarcinomas (3). Hansen et al. compared 18F-FDG and 61Cu-ATSM distributions in 6 sarcoma and 3 carcinoma canine tumors and found moderate to strong voxel-based correlations but reported no comparisons of histologic types of tumors (22). In studies of 18F-FLT and hypoxia correlations in tumors, various degrees of correlation were found in human oropharyngeal tumors (14), high correlations were found in a 9L gliosarcoma rat model (23), and low correlations were found in non–small cell lung cancer xenografts in mice (24). For intratumor correlations of 18F-FDG and 18F-FLT, McCall et al. reported high correlations in oropharyngeal tumors (13), whereas Huang et al. reported low levels of 18F-FDG–18F-FLT overlap in non–small cell lung cancer xenografts (24). In other studies, 18F-FDG uptake was found to be strongly influenced by the proliferation rate (25). Coupled with our findings, these reports suggest that in different histologic types of tumors, the relationships of biologic properties are substantially different—although trends do appear to occur within individual histologic types.
In the analysis of overlapping volumes, thresholds of 90% and less than or equal to 30% did not result in significantly different degrees of overlap of the histologic types of tumors. This finding was an expected consequence of the use of SUV thresholds: high SUV thresholds result in very small volumes and therefore very little overlap of tracer-avid volumes; low thresholds create tracer-avid volumes that encompass most of the GTV, resulting in almost complete overlap. An additional drawback to the use of thresholds is their sensitivity to the SUVmax, which can be highly variable because of SUV uncertainties (26). On the other hand, voxel-based correlation coefficients are insensitive to the SUVmax but can be sensitive to registration, noise, and the reconstruction algorithm used (13).
There are uncertainties in the interpretations of the results of the present study, particularly for application to human cancers. Although humans and dogs have many biologic similarities—they are more similar in growth rate, vasculature, treatment response, and gene expression profile than humans and mice (5)—the use of anesthesia and 100% O2 during imaging may affect tracer uptake (27,28). 61Cu-ATSM also has a lower binding affinity for serum albumin in dogs than in humans, potentially influencing uptake (29). Furthermore, we grouped together all sarcoma tumors (6 chondrosarcomas and 1 osteosarcoma) and all carcinoma tumors (12 adenocarcinomas and 1 squamous cell carcinoma) for analysis because of their similarities in tracer uptake patterns. However, we recognize that genetic profiles can vary drastically in different sarcomas and carcinoma subtypes (30); therefore, phenotypic patterns may show different trends for different subtypes. Finally, there is uncertainty about the specificity of 61Cu-ATSM for hypoxia. Understanding of the mechanism of 61Cu-ATSM uptake has evolved over time (31), but this uptake was recently shown to be caused by impaired mitochondrial electron transport chain function—a consequence of cellular hypoxia (32). Nonetheless, some studies have reported cell line dependence in the specificity of 61Cu-ATSM for hypoxia (33). Likewise, for 18F-FDG and 18F-FLT uptake, no tracer is a perfect surrogate for a single biologic process: tracer uptake is influenced by confounding factors, such as cellular density and vascularity.
CONCLUSION
Using 61Cu-ATSM, 18F-FLT, and 18F-FDG PET imaging, we have shown that the interrelationships of hypoxia, proliferation, and metabolism are significantly different in canine sarcoma and carcinoma tumors. The findings suggest that intratumor phenotypic heterogeneity can vary significantly across different histologic types of tumors. Consequently, therapies that target a single biologic process may result in different outcomes for different histologic types of tumors.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. This work was funded by NIH grant R01 CA136927. No potential conflict of interest relevant to this article was reported.
Acknowledgments
We thank veterinary oncologist Lyndsay Kubicek for overseeing the imaging sessions and for her help in coordinating the study. We also thank members of the UW Image-Guided Therapy Group and the Veterinary Medical Teaching Hospital who assisted during PET/CT imaging and members of the UW Cyclotron Research Group for providing 18F-FLT and 61Cu-ATSM.
Footnotes
Published online Sep. 16, 2013.
- © 2013 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication March 4, 2013.
- Accepted for publication June 17, 2013.