Diffusion MRI and Novel Texture Analysis in Osteosarcoma Xenotransplants Predicts Response to Anti-Checkpoint Therapy

Combinations of targeted drugs have been employed to treat sarcomas, however, response rates have not improved notably, therefore emphasizing the need for novel treatments. In addition, imaging approaches to assess therapeutic response is lacking, as currently measurable indices, such as volume and/or diameter, do not accurately correlate with changes in tumor biology. In this study, quantitative and profound analyses of magnetic resonance imaging (MRI) were developed to evaluate these as imaging biomarkers for MK1775 and Gem in an osteosarcoma xenotransplant model at early time-points following treatment. Notably, we showed that Gem and Gem+MK1775 groups had significantly inhibited tumor growth by day 4, which was presaged by elevations in mean ADC by 24 hours post treatment. Significant differences were also observed at later time points for the Gem+MK1775 combination and MK1775 therapy. ADC distribution and entropy (randomness of ADC values) were also elevated by 24 hours following therapy. Immunohistochemistry demonstrated that these treatment-related increases in ADC correlated with apoptosis and observed cell condensations (dense- and exploded bodies). These findings underline the role of ADC as a quantitative imaging biomarker for therapy-induced response and show promising clinical relevance in the sarcoma patient population.


Introduction
Sarcomas are malignancies of mesenchymal origin that develop within the supporting connective tissues of the body. Although rare in adults (i.e. 1% of all adult cancers), sarcomas are prevalent in children and young adolescents and account for 15% of all childhood cancers [1]. Sarcomas are classified as soft tissue cancers (myoleio, myxoid, Ewing's, lipo), gastrointestinal stromal tumors (GIST) and primary bone cancers (osteosarcomas), with the latter being the focus of this work.
Osteosarcoma, the most common type of malignant bone cancer, has a peak incidence in the second decade of life and affects more males than females [2,3]. Osteosarcoma arises within bone and may, with time spread to the lungs, making the need for early detection and systemic chemotherapy treatment crucial to maximize survival. Similar to other cancers, the current mainstay of osteosarcoma therapy is surgical resection of the primary lesion and chemotherapy. The development of novel therapeutics for sarcomas is hampered by the low incidence and the inter-and intratumoral heterogeneity. Hence, although combinations of novel targeted drugs have been investigated and used, cure rates have only modestly improved in the last decades [4]. Thus, to improve survival and quality of life in sarcoma patients, there is a pressing need to develop novel therapies that exhibit high efficacy.
Previously, we have evaluated the effectiveness of a novel Wee-1 cell cycle checkpoint inhibitor, MK1775, in combination with gemcitabine for its potential to treat osteosarcomas [5]. In a subsequent study, we observed that MK1775 led to cell death in all sarcoma cells and tumor explants independent of their p53 status and showed strong synergistic interaction with gemcitabine [6]. Furthermore, we demonstrated that in a patient-derived osteosarcoma xenotransplant mouse model, MK1775 alone, and in combination with gemcitabine, caused marked terminal differentiation, apoptotic cell death and increased DNA damage by day 14. Because of the aforementioned challenges to measure response by traditional imaging modalities, we have investigated alternative methods that can predict clinical response to MK1775. In this current study, we investigated diffusion MRI at 7 T to evaluate the therapeutic effects of MK1775 and gemcitabine in patient-derived xenotransplant mouse models, particularly at early time-points.
MRI is the modality of choice for evaluating sarcomas and preand post gadolinium enhanced T 1 -and T 2 -weighted multiplanar imaging is currently the clinical standard-of-care (SOC). Due to its excellent tissue contrast and ability to provide anatomical detail in multiple planes, MRI is used in SOC to evaluate measures of response, such as size, necrosis and presence of contrast enhancement. However, it is difficult to correlate changes in these parameters with response to treatment and patient survival in sarcoma [7]. This is likely due to the inherent heterogeneous nature of these tumors. Furthermore, therapy-induced volume changes often occur later and cannot be reliably used as a measure of response during the early course of treatment. Therefore, the need for novel modalities that can assess functional and biologic information more accurately and at an earlier stage in the treatment is of utmost importance in order to affect patient care.
As such, more sophisticated MR approaches including dynamic contrast enhanced MRI (DCE-MRI) [8][9][10] and 1 H magnetic resonance spectroscopy (MRS) [11,12] have been used to predict therapeutic response in a variety of cancers, including sarcoma. In a recent study Guo et al showed that DCE-MRI could be an early imaging indicator of osteosarcoma treatment outcomes; viz. histologic response, event-free survival (EFS) and overall survival) [10]. An enduring challenge in these studies, though, is standardization and reproducibility of DCE-MRI [13], which generally requires an accurate choice of input function [14]. Nonetheless, this highlights the potential role of functional MRmodalities as biomarkers of response for osteosarcoma.
An appealing alternative MRI technique is diffusion-weighted MRI (DWI), which generally has a lower test-retest variance [15]. DWI provides information about the extent and direction of restricted water motion. If data are acquired with different diffusion-weightings (b-values), water mobility in biological tissue can be quantified to yield Apparent Diffusion Coefficients (ADC). It has been shown in numerous studies that the ADC is strongly affected by temperature, extra-cellular space as well as tortuosity, cellularity and integrity of cell membranes [8,16,17]. By quantifying ADC, DWI may therefore serve as a quantitative biomarker for determining lesion cellularity prior to treatment and subsequent changes in cellularity (cytotoxic or vasogenic edema) or loss of cell membrane integrity that may occur in response to successful therapy. Several previous studies have demonstrated increases in ADC values at early time-points following treatment correlated to cytotoxicity [18][19][20], while others have shown early decreases following treatment of anti-angiogenic therapy to reduce vasogenic edema [21,22]. Responses of ADC to therapy in sarcoma have been inconsistent [23,24] and we propose that such studies are limited by evaluating the entire tumor as a region of interest (ROI). This is especially true in sarcomas, which demonstrate significant intratumoral heterogeneity and thus variations in response to therapy and clinical behavior [25].
Thus, to establish a sensitive MRI-based approach for evaluating early chemotherapeutic response, we acquired T 2weighted Fast Spin-Echo (FSE) datasets for quantification of tumor volume and diffusion-weighted sequences (DW-FSE) to assess ADC properties within the tumors and potential therapy-induced alterations therein. To complement mean ADC calculations, we also quantified incremental ADC map pixel fractions, area under the difference (post response) ADC curves and the lack of symmetry in the ADC distribution (i.e. skewness and kurtosis) prior to and after, treatment. Describing the randomness in ADC values, ADC entropy also was determined to investigate potential therapy-induced differences in tumor texture. As a final corroboration of the in vivo MR data, histological staining with cleaved Caspase 3 (CC3) for apoptosis, cH2AX for DNA damage and standard Hematoxylin & Eosin (H&E) was performed and quantified for tissue obtained by day 4.

Experimental Setup
Animal procedures and xenotransplant protocols were approved by the Institutional Animal Care and Use Committee and Institutional Review Board (University of South Florida). As previously described, fresh tissue was obtained from a chemonaïve, 52 year old osteosarcoma patient at the time of initial core biopsy of a distal femur osteosarcoma, which was metastatic to the lungs at the time of presentation [6]. The patient provided written informed consent. Patient tumor was implanted subcutaneously into both flanks of 16 six-week-old female SHO/SCID athymic mice (Charles River Laboratory). At a tumor volume of 500 mm 3 , mice were randomly assigned to the following treatment groups with N = 4 and thus 8 tumors per group: (1) control; (2) MK1775 (30 mg/kg. p.o. twice daily on days 1, 3, 8 and 10); (3) gemcitabine (100 mg/kg, i.p. once daily on days 1, 3, 8, 10); (4) MK1775+ gemcitabine (Combo) in the above mentioned doses. Drug doses and schedules were chosen based on the investigator brochure [6,26,27].
Prior to imaging, mice were anesthetized using 2% isoflurane and restrained in a mouse cradle. Respiratory function was monitored using the SAII system (Small Animal Instruments, Inc.) and temperature control was achieved by a gas unit, set to maintain a body temperature of 3761uC.

Imaging Protocol and Analysis
MRI was performed at baseline and 24 hours following each drug administration, and on day 14. All data were acquired using a 7-T horizontal magnet (ASR 310, Agilent Technologies, Inc.) equipped with nested 205/120/HDS gradient insert and a bore size of 310 mm. Using a 72-mm quadrature birdcage coil (Agilent Technologies, Inc.), axial T 2 -weighted fast spin-echo (FSE) sequences were acquired (TE/TR = 60/1403 ms) with slice thickness of 1.5 mm and in-plane resolution of 136 mm over 6 minutes. Applying identical slice plane and spatial resolution, diffusion weighted (DW) FSE sequences using four b-values = 50, 500, 1000 and 2000 and TE/TR = 36/1881 ms also were acquired over 12 minutes.
Image reconstruction and volumetric analyses were performed in VnmrJ (Agilent Technologies, Inc.) while ADC analysis used inhouse routines in MATLAB. Manually drawn regions of interest (ROIs) encompassing the tumors were applied to the T 2 -and diffusion-weighted datasets for tumor volume and ADC. Tumor ADC was calculated using the following equation: S(b) = S(0)6e 2ADC*b where b is the diffusion sensitizing factor, S(b) and S(0) the signal intensity with and without diffusion weighting, respectively. Percentage change in mean ADC was calculated according to: DADC [%] = (ADC t /ADC 0 )6100, where ADC 0 denotes the baseline ADC value. Similarly, the percent change in tumor volume was calculated as DV [%] = (V t /V 0 )6100 with V 0 being the tumor volume obtained pre-treatment. ADC maps were generated by nonlinear least squares regression of a mono-exponential to the experimental signal intensity for all four b-values and.
To compare across ADC values in a given tumor, an incremental lower ADC limit was set at 3610 25 mm 2 /s, and by 512 increments, all values above the lower threshold with an upper limit of 5610 23 were quantified and plotted to compare pre-and post treatment for all time points and mice [28]. Therapy-induced alterations in ADC values were further characterized by area under the curve (AUC) analysis, which was performed by subtracting the pre-treatment ADC pixel fraction curve from the post-treatment curve for each animal. Changes in skewness and kurtosis also were quantified and are displayed in histograms showing the difference between baseline values pre-treatment and 24 hours following treatment.

Texture Analysis
Visual analysis demonstrated intra-tumoral change and heterogeneity in the diffusion-weighted datasets, which appeared to correlate with treatment. To quantify these variations and subtract image texture features (i.e. distinct regions or repeating patterns of pixelated behavior), entropy filtering was applied to the ADC maps [29] [30]. Entropy filter is the implementation of Shannon entropy; a texture-based statistical measure of the variation in the histogram distribution, which reflects the predictability of intensity values within a given ROI. The region was manually delineated to encompass the entire tumor and a set of consecutive masks was generated, thus forming a volume of interest (VOI).
The entropy of the ADC maps was computed using the values within ADC map VOI. A custom entropy algorithm was implemented in MATLAB. Due to size limitations of the ADC maps (1286128 pixels per horizontal slice), entropy computation was performed globally over the entire VOI and the ADC values were used to form a histogram vector, which was subdivided into four discrete bins. For each tumor, a single entropy coefficient was computed and used as a measure of variability within that VOI.
Intensity variations within individual VOIs were visualized using 3D graphs of entropy coefficients. Instead of computing one entropy score for the entire tumor, entropy was computed for each pixel based on the ADC values of its neighboring pixels (969 window), which were also within the tumor VOI. Each pixel was then color-coded based on the value of its local entropy coefficient.

Histological Evaluation
Employing the identical experimental setup, animal experiments were performed in two sets of mice. The first batch (N = 16), focused on the cytotoxic effects of MK-1775 and Gem in osteosarcoma xenografts through day 14, with this as the immunohistochemical endpoint [6]. Quantitative analysis of these datasets demonstrated significant ADC increases by 2 to 4 days following therapy. Therefore, the current study (N = 16) focused on these early times-points with immunohistochemistry performed on day 4. Histological processing and staining was performed as previously described [6].
H&E stained histology slides of xenograft tumors, cH2AX and CC3 were scanned using the Aperio TM (Vista, CA) ScanScope XT with a 2006/0.8NA objective lens at a rate of 5 minutes per slide via Basler tri-linear-array. Histological pattern recognition was conducted using the GenieH (Aperio, Vista, CA, USA) software platform to segment and classify necrosis, viable tumor and other non-target tissues (i.e. adipose, muscle and skin). In addition, TissueStudioH (Definiens, Munich, Germany) image analysis platform was used to segment, classify and extract size and intensity features for single cells in the viable tumor and necrotic regions of each image.
Quantification of CC3 stained exploded bodies was performed by spectrally unmixing images into hematoxylin and DAB channels using the TRIO image acquisition and analysis system (Caliper; Boston, MA). Each channel was rerun through the segmentation and classification algorithm described earlier and the number of exploded bodies, with, and without CC3 staining was recorded.

Statistical Analysis
Statistical analysis was performed using Graphpad Prism software (Graphpad, San Diego, CA, USA) and one-way analysis of variance (ANOVA) followed by the Tukey test for comparison of mean values. A confidence interval of 95% was chosen and thus statistical significance was pre-determined at p,0.05. For each measured parameter (i.e. tumor volume, mean ADC values, ADC-AUC, entropy and histological scoring), the averaged mean value for each of the four cohorts on a specific day was compared to one another and the results are presented in Table 1. Figure 1A shows that high resolution FSE datasets allow for the delineation of the osteosarcoma tumors for all animals at all imaging time points. The borders of each tumor were distinctly visible making the lesions easy to distinguish from surrounding tissue and thereby facilitating segmentation for volume and diffusion properties. In addition, these images also displayed heterogeneity and microstructural details within the tumors, particularly with an increase in tumor size. Figure 1A further displays the difference in tumor sizes between the various groups. Figure 1B shows the time course of tumor volume quantifications, which showed that the Gem and Combo groups appeared to exhibit decreased tumor growth relative to the MK1775 and control tumors. Both Gem and Combo groups had significantly lower tumor volumes on day 4, but by day 14 only the Combo group was significantly lower than the untreated controls (Table 1A). The MK1775 group grew similar to the controls up until day 11, but appear to have decreased tumor volumes by day 14. However, this can be explained by having only 2 remaining quantifiable imaging acquisitions from this group on day 14 and this possible decrease cannot be statistically verified.

Quantification of ADC, Skewness and Kurtosis
To investigate the response of apparent water diffusion across the tumors to therapy, ADC maps were generated for all animals. Following established protocols, animals were treated on days 1, 3, 8 and 10. Focusing on the early time points, Figure 2A demonstrates a notable increase in ADC by 24 hours following the first and second treatment (i.e. day 2 and 4, respectively) in the Gem and Combo groups while the controls and MK1775 showed that the ADC was unaltered. To quantify this effect, each animal was used as its own control, and the mean ADC was estimated, then group averaged and displayed as percent change from the baseline value with standard deviations. As shown in Figure 2B, the ADC of the control group and MK1775 did not change significantly over the course of the study. The Gem group had a dramatic and significant increase in ADC evident by day 2, which persisted until day 9, and increased slightly following the last dose on day 10. Similar to the Gem cohort, the combo group demonstrated a significant increase in ADC following initiation of therapy, but in contrast to Gem, the ADC remained elevated until 5 days following the last treatment, at which time these tumors had started to re-grow (cf. Figure 1B). Statistically, significant differences in ADC were established for the controls compared to Gem and Combo on day 2, as well as controls and Combo, and MK1775 and Combo, on day 9 (Table 1B).
ADC histograms were generated for each mouse prior to and following therapy and using a similar approach as Galons et al, these were used to generate cumulative pixel fraction plots [28]. Subtracting the pre-from the post-therapy pixel fraction plots is a sensitive measure of movement to higher ADC values and the area under this resulting curve (AUC) was calculated and plotted as the average value for each treatment group as a function of time ( Figure 2C). These results were in agreement with both volumetric and mean ADC analysis and demonstrated significantly larger shifts in AUC for the Gem and Combo groups than MK1775 and controls at most time points (Table 1C). Statistical significance was reached between control and Combo on day 2 and 9, and control and Gem on day 4.
To further characterize diffusion properties within the tumors, additional descriptive statistical measures of the ADC distribution; viz. skewness and kurtosis, were determined. Skewness and kurtosis describe the lack of symmetry and the breadth of the ADC pixel distribution, respectively, and alterations therein are direct indicators of ADC changes across tumors. Figures 3A-D show Table 1. P-values obtained from a one-way analysis of variance (ANOVA) followed by the Tukey test for comparison of mean values between the stated groups on the specified time points for (A) tumor volumes, (B) percent change in mean ADC, (C) ADC-AUC analysis and (D) CC3 and cH2AX staining.  the raw ADC histograms of the pixel distributions pre-and 24 hours post the first treatment. Whereas the histogram for a representative control animal showed similar ADC pixel distributions prior to and following treatment, it is evident that the treated animals showed alterations within 24 hours after the first drug administration. In particular, the ADC distribution post-treatment appears to be markedly skewed and have a broader peak for the Combo, and to a lesser extent the Gem and MK1775 groups. The ADC histograms thus correlate to the incremental ADC pixel fraction plots, which demonstrate a shift of ADCs towards higher values in Gem and Combo treated animals compared to MK1775 and controls (Figures 3 E-H). Quantitatively, the relative difference in skewness and kurtosis for each group 24 hours following treatment is shown in Figures 3I  and 3J. In agreement with the mean ADC and AUC data at the same time point ( Figure 3K & L), all three treated groups, and Combo in particular, demonstrated negative skewness and kurtosis following therapy while the control group showed a minor positive increase. Although this negative trend was evident for all three treated groups, statistical significance could only be established between controls and Combo (p = 0.018 and p = 0.043, skewness and kurtosis respectively) at 24 hours.

Tumor Microstructure and Entropy
Entropy is a texture-based statistical measure of the randomness in an ADC histogram such that a large variation in ADC is associated with higher entropy values. Presented in Figure 4A, post treatment ADC entropy plots on day 2 demonstrate evident differences between the various groups. The most apparent distinction is between the controls and the Combo and Gem groups, with the controls showing large regions of substantially lower entropy values (blue color). In addition, the entropy distributions varied considerably for all treated groups, and especially so in Gem and Combo in particular, compared to controls. In agreement with this visual assessment, the quantitative data demonstrated in Figure 4B shows a significantly larger change in ADC entropy with treatment for the Gem (34.2%) group compared to both MK1775 (11.6%, p = 0.023) and controls (10.6%, p = 0.022). While the Combo group also displayed a 30.2% increase in ADC-entropy compared to pre-treatment values, statistical significance could not be established.

Quantitative Immunohistochemistry (CC3, cH2AX and H&E)
For the evaluation of drug-mediated cytotoxic and apoptotic effects, immunohistochemical analyses were performed with cleaved caspase 3 (CC3) antibody and phosphorylated histone (cH2AX). As a marker for apoptotic cell death, CC3 activity was significantly higher on day 4 in the Gem and Combo groups compared to controls but were not elevated in the MK1775 group ( Figure 5A and B, Table 1D). The cH2AX staining appears elevated in the GEM and Combo groups ( Figure 5C), however these differences are not statistically significant.
Standard H&E sections also were analyzed, however, no visual or quantitative differences with regards to potential bone formation, necrosis or average cell area were detected at day 4 ( Figure S1). Similarly, quantitative analysis of cellularity did not statistically differ between groups at this time-point, however, a trend was observed correlating treatment response with decreased cellularity ( Figure S1B). Interestingly, as determined by microscopic analysis of the H&E stained sections, we noted a varying amount of small dark bodies that seemed to be particularly increased in the Gem and Combo groups ( Figure 6A). Verifying these observations, a quantification of these dense units demonstrated a higher percentage in all treated groups, but Gem in particular (p = 0.038), compared to controls ( Figure 6B). To further evaluate these unknown structures, the fractions of these bodies that stained positively for CC3 were determined. Although all three treated groups contained more apoptotic dense bodies than the controls, statistical significance could not be established between any of the groups.
We also observed dark dense bodies with exceedingly large and bulgy masses that were increased in the Gem and Combo groups compared to controls ( Figure 6C). Classified as ''exploded dense bodies'', quantification of these structures demonstrated significantly higher levels in the Gem (p = 0.048) and Combo (p = 0.004), than in controls ( Figure 6D). The Combo group also showed significantly higher burden of these exploded units compared to the MK1775 treated mice (p = 0.032). To evaluate whether these exploded bodies stained positively for CC3 and thus comprised the apoptotic portion of these dense units demonstrated in Figure 6C, spectral unmixing image analysis was used to deconvolve and register the two signals and determined that 97% of these structures were fully apoptotic ( Figure 6E & F). No intra-group differences were detected.

Discussion
In this study, we used diffusion MRI to evaluate the response of osteosarcoma to Gem and MK1775. However, the focus of the current study was to evaluate quantitative approaches to assess biologically relevant therapy-induced responses at earlier timepoints than conventional volumetric analysis. Specifically, we investigated multiple descriptors of ADC parameters and compared these to subsequent changes in tumor volumes and tumor microstructure. Significant effects of therapy on tumor volume were observed by day 4 (72 hours after first dose) and the volume changes correlated with ADC. Significant increases in mean tumor ADC and ADC distribution could be detected within 24 hours following initiation of treatment, and these presaged subsequent volume responses. The Combo group had elevated ADC and tumor volume responses that were sustained at later time points relative to the GEM treatment alone. Quantitatively, these observations are in agreement with previous findings wherein ADC increased shortly after initiation of therapy in both humans and animal models [8,17,20,[31][32][33]. These rapid increases in ADC following treatment have been ascribed to reduced cellularity associated with cytotoxicity, as well as impaired cell membrane integrity associated with necrosis, both of which allow for more unimpeded diffusion [31,34,35]. Another contributing factor to elevated ADC values may be cell shrinkage due to increased interstitial osmolality or apoptosis [36,37]. It should be noted, however, that some therapies can result in a decrease in ADC, through reduction of vascular permeability by anti-VEGF therapies, which leads to reduced vasogenic edema and reduced interstitial osmolality [38]. Hence, interpretation of ADC data requires accurate knowledge of the biology of the tumor and the anti-tumor therapy mechanism of action. Regardless of mechanism, increased ADC is interpreted to be synonymous with increased interstitial space, and vice-versa.
In this study, histological analysis did not demonstrate significant therapy-induced necrosis by day 4, however, a trend of reduced cellularity as well as a larger range in cell sizes for the treated groups was observed, which may be explained by the presence of the unusually smaller dense and exploded bodies. More important, these findings correlated with the increases in mean ADC. For all animals, this increase in ADC also was associated with a negatively skewed ADC distribution as represented by the skewness and kurtosis. Consistently, the Combo group had more exploded bodies and greater perturbations in skewness and kurtosis of the ADC histogram compared to GEM or MK1775 alone.
Skewness and kurtosis are statistical measures of the asymmetry and breadth of a distribution, respectively. Theoretically, if the entire distribution moves to higher values, the mean ADC should yield the greatest statistically significant changes, as was the case in this current study. Previously, assessments of ADC skewness and kurtosis have been applied for various tumors; however, its usefulness thus far has been limited to the brain [39][40][41]. Capturing a different microenvironment for diffusivity, ADC histograms provide information that may be masked by mean ADC values [20], and should therefore complement standard ADC analysis. Further investigating ADC distribution, we found incremental ADC pixel fraction plots and ADC-AUC analysis informative in establishing therapy response. Providing quantitative information about ADC distribution, AUC analysis comparing pre-and post-treatment plots, showed statistically significant differences between groups, which also correlated with volume, mean ADC, skewness and kurtosis.
Perhaps more intriguing, ADC-entropy measurements and 3dimensional entropy plots demonstrated the intra-tumoral ADC variability and correlated with tumor growth, mean ADC, ADC distribution, AUC analysis and histology. Describing the variability of ADC values within a certain ROI, ADC entropy calculations were first reported in the literature by Benedict et al in 2007 and was therein recognized as a sensitive marker of overall brain damage that strongly correlated with cognitive impairment [29]. More recently, a study by Kierans et al. demonstrated that ADC entropy showed significantly greater accuracy than mean ADC for characterization of malignancy of adnexal lesions [30]. In this Diffusion MRI Texture Predicts Response in Sarcoma PLOS ONE | www.plosone.org current study, ADC entropy was able to predict treatment response following cancer therapy at earlier time-points than tumor volume changes. To our knowledge, this application of ADC-entropy has previously not been investigated nor reported.
The final segment of our approach evaluated potential mechanisms responsible for these observed therapy-induced changes in ADC and ADC distribution. On day 4 following therapy, we found that CC3 staining was significantly elevated in the Gem and Combo groups relative to controls and MK1775, indicating that apoptosis was a mechanism of rapid treatment response. This correlates well with the peak increases in ADC that occurred on day 4 for these two groups.
Equally compelling, was a significant increase in condensed cell bodies in the Gem group relative to controls. A subset of these dense bodies had an ''exploded'' appearance, and were highly correlated with all three treatment groups. As determined by quantitative image analysis, 98% percent of these dense exploded bodies stained fully for CC3 and, thus, apoptosis most likely is the cell-death mechanism leading to the formation of these condensed bodies. Since the presence of CC3 is a temporal event in the pathway to apoptotic cell death, the inability to detect CC3 in all dense bodies does not rule out apoptosis as being the sole mechanism. However, it is possible that apoptosis occurs in only a subset of cells undergoing cell condensation.
Interestingly, while MK1775 did not induce significant apoptosis or elevated ADC values at day 4, our studies have shown that when administered for a prolonged time period (i.e. 10 days or longer), there is a cumulative therapeutic effect of MK1775 as a mono-therapy. In addition, histological differences indicating differentiation (osteogenesis) and DNA damage are clearly evident by day 14. Therefore, our results suggest that ADC is a predictive biomarker for acute cytotoxic effects, however, not as sensitive to predict slower, cumulative processes that occur at later time points.
In conclusion, therapeutic effects of Gem, administered solely or in combination with MK1775, induced significant tumor growth control by day 4, which was presaged by significant elevations in mean tumor ADC, ADC distribution and entropy by 24 hours following therapy. Treatment related increases in ADC, CC3 staining and observed cell condensations indicated that apoptosis was a significant component of the response to treatment regimens that included Gem and that the increased ADC values in this study correlated with the apoptotic response. At early time points, the combination GEM+MK1775 showed significantly increased responses relative to GEM or MK1775 alone by multiple measures. As expected, if tumor control occurs more slowly, e.g. with MK1775 monotherapy, there are no significant changes in ADC or ADC distribution, although there is a slight trend to higher values. We interpret this observation to indicate that therapy-induced edema can be systemically absorbed faster than it is created. These results show promising clinical relevance in the sarcoma patient population as agents that can induce a rapid