Analysis of the intra- and intertumoral heterogeneity of hypoxia in pancreatic cancer patients receiving the nitroimidazole tracer pimonidazole.

BACKGROUND: Hypoxia is thought to be an adverse feature of pancreatic cancer, but direct measurement in patients is technically challenging. To address this, we characterised the intra/interpatient heterogeneity of hypoxia in surgical specimens from patients who received the 2-nitroimidazole tracer pimonidazole pre-operatively.
METHODS: Pimondazole was given intravenously 16-20 h before pancreatectomy, and the extent and intratumoral heterogeneity of hypoxia determined by image analysis applied to multiple tissue blocks stained by immunohistochemistry. Intra/interpatient heterogeneity was estimated by variance component analysis.
RESULTS: Pimonidazole staining was analysed in 10 tumours. The extent of labelling varied amongst patients (0-26%), with a broader range of hypoxia in the epithelial (1-39%) compared with the stromal (1-13%) compartments. Variance component analysis demonstrated greater inter- than intrapatient variability of hypoxia, and that multiple (4-5) tumour sections are required to provide a consistent evaluation of its extent in individual tumours.
CONCLUSIONS: There is significant intra- and intertumoral heterogeneity of hypoxia in pancreatic cancers, and these do not appear to be generally more hypoxic than other cancer types. This study establishes the feasibility to assess hypoxia in pancreatic cancer patients using pimonidazole, but questions the reliability of measurements made using a single tissue section.

Pancreatic ductal adenocarcinoma has the poorest treatment outcome of the common cancers, and current projections estimate that it will become the second leading cause of cancer death in western countries after lung cancer by 2030 (Rahib ). Early whole-genome and targeted sequencing studies have partially revealed the extent of genetic complexity inherent in this disease, suggesting that highly specific targeting of molecular aberrations unique to epithelial malignant cells is unlikely to result in sustained benefit for large groups of patients (Jones ; Alexandrov ). Contributing additional complexity is the tumour stroma—a heterogenous collection of fibroblasts, immunomodulatory and vascular elements embedded in a biochemically diverse extracellular matrix (ECM) believed to be of particular relevance to pancreatic ductal adenocarcinoma (PDAC) tumorigenesis (Neesse ).

The dense desmoplastic stroma that characterises PDAC has been proposed as a determinant of biological aggression, promoting invasion and metastasis as well as resistance to standard therapies (Olive ; Provenzano ). Several mechanisms have been suggested, including compromised tumour perfusion that limits drug penetration; the promotion of an acidotic, hypoxic microenvironment; and the activation of tumour-promoting signalling by ECM components such as collagen-1, the matrix metalloproteinases and TGF-β (Armstrong ; He ; Olive ; Neesse ).

The development of microregional hypoxia in solid tumours is well described with higher levels of hypoxia being associated with poor patient outcome due to inferior treatment response and more aggressive tumour behaviour (Fyles ; Brown and Wilson, 2004; Nordsmark ; Rischin ; Vaupel, 2008; Overgaard, 2011). PDAC has been historically considered unusually hypoxic, based on poor perfusion by radiological contrast agents ((Fusaroli ; Matsubara ) and Eppendorf electrode probe pO2 measurements made intraoperatively in a series of seven patients with resectable pancreatic cancers (Koong ). Still, although preclinical data support tumoral hypoxia as a relevant factor in the aggressive biology of PDAC (Buchler ; Spivak-Kroizman ), this remains to be clinically validated.

Using orthotopically grown patient-derived xenograft models of PDAC, we observed a wide range in hypoxia, and a striking correlation between higher levels of hypoxia and increased cellular proliferation, tumour growth rate and greater metastatic potential (Chang ). These results support the concept that hypoxia is associated with biological aggression of PDAC. To establish the clinical relevance of these findings, we initiated the PIMO-PANC trial (NCT01248637), which administers the 2-nitroimidazole hypoxia tracer pimonidazole to patients pre-operatively. Pimonidazole undergoes bioreductive metabolism under low oxygen conditions to form stable adducts that can be detected in tissue sections by immunohistochemistry (Raleigh ; Varia ) (Bussink ). This ongoing study has a target accrual of 100 patients, and is powered to assess the prognostic importance of hypoxia in early stage pancreatic cancer. In the present paper, we describe an interim analysis of the intra/interpatient heterogeneity of hypoxia, and the importance of sampling error due to intratumoral heterogeneity when assessing tumour hypoxia using histological markers. Our results suggest that this is a larger problem than is often appreciated.

Materials and methods

Materials

Pimonidazole hydrochloride (Hypoxyprobe-1) was obtained from Natural Pharmacia International (NPI) Inc, Belmont, MA, USA and has IND (Investigational New Drug) status for use for the clinical evaluation of hypoxia. Immunohistochemical (IHC) detection of stable pimonidazole protein adducts was performed using (mouse) monoclonal IgG1, Hypoxyprobe MAb1 (Raleigh ; Raleigh ), biotin blocking kit (Vector Lab, Burlington, ON, Canada), streptavidin biotin detection (Signet Pathology System, Deham, MA, USA), Nova Red substrate (Vector Lab), Mayer's haematoxylin counterstaining and Permount mounting medium.

Patients

Approval was obtained from the University Health Network Research Ethics Board. The eligible patient cohort included those scheduled for surgery with a presumptive diagnosis of pancreatic adenocarcinoma. Male and female patients, over 18 years of age and of all ethnicities, were included. Written informed consent was obtained before patient registration.

Drug administration

Enroled patients received one dose of intravenous pimonidazole (0.5 gm m−2) over a minimum of 30 min on the day before surgery. Pimonidazole was administered on an out-patient basis under medical supervision with infusion being completed 16–20 h before surgery. Manufacturer's recommendations are for tumour sampling to take place at 16–24 h post pimonidazole administration ensuring minimal levels of circulating pimonidazole at the time of tumour harvesting. Pimonidazole labelling is therefore reflective of in vivo hypoxia status of tumour and not of hypoxia sustained during tumour harvesting.

Tissue processing

Resected tumours were processed for paraffin embedding and histological examination according to standard pathological practice, yielding multiple paraffin blocks from the primary tumour and regional lymph nodes. This material was reviewed (ND, SS), and blocks representative of primary tumour were selected. Serial sections were stained for haematoxylin and eosin (H&E) and pimonidazole (monoclonal IgG1 antibody, Hypoxyprobe MAb1). The IHC protocol was followed as described previously and outlined below (Raleigh ). Stained sections were digitised using an Aperio ScanScope XT (Aperio Technologies, Vista, CA, USA) whole slide scanner at × 20 magnification (resolution 0.5 μM per pixel).

Pimonidazole IHC

Unless specified, all steps were performed at ambient temperature with PBS washes between incubations. Paraffin sections were dried in a 60 °C oven overnight, dewaxed and endogenous peroxidase blocked by incubation in 3% hydrogen peroxide for 10 min. Following microwave antigen retrieval, endogenous biotin was blocked (Vector blocking kit). Slides were incubated overnight in pimonidazole mouse monoclonal antibody followed by biotinylated anti-mouse IgG, streptavidin biotin (30 min each) and Nova Red substrate (5 min), then counterstained in Mayer's haematoxylin.

Histological tumour analysis

Image analysis included quantitation of (1) proportion of epithelial and stromal tumour compartments and (2) hypoxic fraction in whole tumour (i.e., in combined epithelial and stromal regions) as well as in individual compartments. ‘Hypoxic tumour areas' were defined as regions demonstrating any positive immunostaining for pimonidazole. Images were managed on Spectrum, and analysis was completed using Genie, Aperio's pattern recognition software. This image analysis platform differentiates different histological features through computer-based learning algorithms by running repetitive training sequences using examples of different cellular morphologies to define ‘classes.'

H&E sections were first reviewed to identify appropriate tumour regions for analysis, excluding large areas of necrosis, atrophic pancreas and non-neoplastic tissues. A region of interest was manually selected for each pimonidazole IHC image based on H&E. The Genie classes ‘epithelium,' ‘stroma' and ‘other' were defined, with ‘other' including necrosis, luminal spaces and any ‘non-tumour' tissue to be excluded from analysis. On the basis of an observer-set threshold/training, a unique Genie classifier was developed for each patient's tumour in order to differentiate epithelial from stromal tumour compartments. Aperio's Positive Pixel v9 algorithm was applied to quantify hypoxic percentages in whole tumour with non-tumour tissue (including necrotic and luminal spaces) defined as ‘other' being excluded and epithelial and stromal compartments separately. Microscopic scoring analysis was also conducted independently by one of the authors (JS), assigning each section a score (in 5% increments) representing per cent hypoxia in the epithelial tumour compartment.

Statistical analysis

On the basis of our previous studies of the intratumoral heterogeneity of hypoxia, we estimated that an analysis of five sections per tumour in 10 patient tumours would provide a reasonable initial estimate of hypoxia levels in resected pancreatic cancers (Thrall ; Iakovlev ; Maseide ; Pintilie ). Hypoxia levels were measured in multiple sections of the same specimen (each from a different tissue block). Variance of hypoxia within and across tumours was calculated using variance component analysis. The intraclass correlation coefficient (ICC), or the per cent variation attributable to that between tumours, was also calculated from simulations. The ICC is a measure of reliability ranging between 0 and 1. A value of 0.85 or higher is considered indicative of good reliability reflecting the fact that the variability between tumours explains most of the total variability. In contrast, low values for ICC reflect poor reliability suggesting that a large proportion of the total variance is due to the variance within a tumour. To estimate the incremental benefit of analysis of multiple sections per patient, we calculated the ICC for the average value of hypoxia for the scenarios when 1–6 sections per patient were to be analysed.

The relationship between hypoxia in epithelial and stromal compartments was assessed using the Spearman correlation coefficient.

Results

Of the first 16 patients who received pimonidazole and proceeded to surgery, 5 patients were intraoperatively found to have advanced disease precluding definitive resection. In these patients with locally advanced or metastatic disease, biopsies taken to confirm/establish histological diagnosis were evaluated qualitatively for pimonidazole staining. Regions of hypoxia, as indicated by positive pimonidazole staining, were observed in the liver and lymph node metastases (Supplementary Fig 1). The group of 11 resected tumours contained 1 common bile duct cancer, leaving 10 PDACs that are the subject of this study. Relevant clinicopathologic characteristics of these patients/tumours are outlined in Supplementary Table 1. The majority (8 out of 10) had a primary tumour located in the head of the pancreas, and (9 out of 10) had lymph node involvement. Four out of the 10 patients completed neoadjuvant treatment before surgery—3 received 2 months of gemcitabine and erlotinib on a perioperative clinical trial (NCT00733746) and 1 concurrent (5-FU) chemoradiotherapy.

Quantitative estimations of the extent of stroma and hypoxia were made in the 10 resected PDACs by image analysis conducted on multiple sections per patient tumour, with each section cut from different blocks representing geographically distinct tumour regions. Number of sections analysed was 3 (in 1 patient), 5 (2 patients), 6 (5 patients) and 7 (2 patients). There was significant heterogeneity of the proportion of tumour stroma both within and across individual tumours with a range (across patients) of 38–78% (per sectional tumour area analysed); the largest intrapatient range was 41–83% (pt 008) (Figure 1).

Figure 1

Heterogeneity of extent of stroma across patient tumours. The table outlines stromal content across 10 resected PDACs represented as mean value across all tumour sections analysed. Included histology are representative images of tumours with high (001: 78%) and moderate (017: 42%) stromal content; scale bar, 500 μM.

The pattern of pimonidazole staining was cord like as previously described as characteristic of nitroimidazole markers (Raleigh ; Raleigh ; Evans ). Positive staining was observed in nuclear and cytoplasmic cellular regions in both epithelial and stromal tumour compartments. Pimonidazole staining in tumour was heterogeneous across and within patient tumours, with adjacent normal pancreas demonstrating no pimonidazole staining consistent with an absence of severe tissue hypoxia (Figure 2). There was no evidence of cytotoxic treatment effect or necrosis in tumours resected from patients receiving pre-operative treatment. Further, pimonidazole staining in treated tumours was similar to that observed in untreated tumours with respect to staining patterns and intensity.

Figure 2

Representative histology (H&E (columns 2 and 4) and pimonidazole immunohistochemistry (IHC) (columns 3 and 5)) for 10 resected tumours included in heterogeneity analysis. (H&E for individual patients are labelled with patient study ID (001 to 017) with corresponding pimonidazole IHC in column to right). Note absence of pimonidazole staining in adjacent normal pancreatic tissue (column 1).

The mean hypoxic percentage (defined as positive pimonidazole staining per tumoral area analysed) in whole tumours (including both epithelial and stromal tumour compartments) was 9.5% with a range of 1–26% (Figure 3). Higher levels of hypoxia were measured in epithelial (mean 19% range: 1–39%)) compared with stromal (mean 5% range: 1–13%) regions (Figure 4) with a significant correlation between the two (Spearman correlation coefficient r=0.75; P=0.011). Positive pimonidazole staining was noted in all the 4 tumours resected from patients who had received pre-operative treatment, 2 at low level (1% in epithelial tumour) and 2 at more moderate levels (15 and 12%, respectively, in epithelial tumour).

Figure 3

Hypoxic percentage (HP) in whole tumour (i.e., including both epithelial and stromal tumour components) defined as positive pimonidazole staining per whole-tumour section analysed. Each dot represents the HP measured in a single section and the bar the mean value from analysis of all sections.

Figure 4

Hypoxia in epithelial and stromal tumour compartments. The table includes mean and range hypoxic percentages in whole tumour, and in the epithelial and stromal tumour compartments. Dot plots show the individual values and the means for the epithelial (A) and stromal (B) compartments.

Variance component analysis demonstrated greater inter- than intrapatient variability of hypoxia as indicated by positive pimonidazole staining (Table 1a). Within a given patient tumour, there was greater heterogeneity of pimonidazole staining within the stromal (in comparison with the epithelial) compartment (as indicated by the lower ICC values in tumour stroma of 0.50–0.91 vs 0.79–0.97 in epithelial tumour), whereas the greater variance between patient tumours was in the epithelial compartment (Table 1b). These data suggest that if tumour sampling is adequate to account for heterogeneity of pimonidazole staining in the epithelial component, it should be feasible to resolve biological differences related to severe hypoxia between tumours.

Table 1a

Results of variance component analysis of heterogeneity of hypoxia in resected tumours. Intra- and interpatient variance of hypoxia (indicated by HP or hypoxic percentage) in whole tumour and epithelial and stromal compartments

	Intrapatient variance	Interpatient variance	Ratio of intra/interpatient variance
Hypoxia in whole tumour	0.66	1.13	0.37
Hypoxia in epithelial tumour	0.49	1.83	0.21
Hypoxia in stromal tumour	0.74	0.74	0.5

Table 1b

Intraclass coefficient (ICC) analysis estimating reliability of measures of HP made in 1–10 tumour sections per patient tumour

Sections analysed per patient tumour	ICC epithelial	ICC stroma
1	0.79	0.5
2	0.88	0.67
3	0.92	0.75
4	0.94	0.8
5	0.95	0.83
6	0.96	0.86
7	0.96	0.88
8	0.97	0.89
9	0.97	0.9
10	0.97	0.91

Analysis demonstrates increased reliability of estimates with increasing number of sections analysed, with analysis of 2 and 6 sections being required to provide an ICC greater than 0.85 in epithelial and stromal tumour, respectively.

The average of repeat measures of a particular tumour feature provides a more accurate description of that particular characteristic than a single measurement. The greater the number of repeats, the greater the accuracy of the estimate becomes as the variance is reduced. We considered estimates of hypoxia in the epithelial tumour compartment provided by measurements made in 1–5 sections per patient tumour. As expected, there was a decrease in variance from 1 to 3 repeat measures made per tumour with further reduction as more sections were analysed. Consistent with our initial hypothesis, analysis of 4–5 sections per patient tumour appears to be able to provide some discriminatory power to differentiate among patients based on extent of pimonidazole staining (Figure 5).

Figure 5

Results of simulations estimating proportion of hypoxic epithelial tumour (indicated by positive pimonidazole staining in epithelial tumour cells) obtained when analysing 1–5 sections per patient tumour. In all graphs, x axis indicates patient ID (in rank order of increasing levels of hypoxia) and y axis indicates hypoxic fraction or proportion (log transformed) with 95% confidence interval. As anticipated, the variance in estimates of hypoxia decrease (indicated by the narrowing of the confidence intervals) as the number of sections analysed increases.

Discussion

Although a number of previous studies have assessed hypoxia based on the expression of surrogate markers such as HIF-1α and carbonic anhydrase IX (CAIX) in pancreatectomy specimens (Couvelard ; Sun ; Hoffmann ; Cheng ), to our knowledge this is the first to use an extrinsic hypoxia tracer administered to the patient pre-operatively, and the first systematic study of sampling error related to its intratumoral heterogeneity. The accuracy of histology-based assessments of molecular markers relies on the extent of tumour sampling, which in turn should be guided by the heterogeneity of marker tissue expression (Pintilie ). We therefore planned our analysis to include sections from multiple distinct tumour blocks. Defining ‘hypoxic tissue' as that demonstrating positive pimonidazole uptake, we observed significant inter- and intratumoral heterogeneity in both epithelial and stromal tumour compartments. Variance component analysis and calculation of intraclass correlation (ICC) coefficient values helped determine the number of repeat measures required to accurately estimate hypoxia. Given that an ICC value of >0.85 is generally considered acceptable with respect to the robustness of measures obtained, we noted that analysis of a minimum of three sections per tumour (ICC: 0.92) was required to estimate hypoxia in the epithelial tumour compartment. Consistent with the higher variance of hypoxia in stroma, analysis of a minimum of 6 sections (ICC: 0.86) appear to be required to estimate stromal hypoxia (Table 1b).

These results raise concerns about the robustness of the assessment of tumoral hypoxia through the use of core biopsies and tissue microarrays in PDAC. They also substantiate prior work demonstrating that the accurate assessment of heterogenous (hypoxia) markers such as CAIX requires analysis of multiple core biopsies obtained systematically from geographically distant tumour locations (either from distinct tumour blocks or different biopsies) (Thrall ; Iakovlev ; Maseide ); strategies that tend to be more extensive than those currently employed (Couvelard ; Sun ). Therefore, there is a trade-off between minimising sampling error due to intratumoral heterogeneity of a histological marker, and the increased costs of analysing multiple tissue sections.

In contrast to a prevalent belief that PDAC are unusually hypoxic, we observed that a larger proportion of each tumour examined did not stain positively for pimonidazole. The most ‘hypoxic' tumours exhibited positive staining in 20–30% of epithelial tumour areas. Furthermore, 2 out of 10 tumours showed minimal levels of pimonidazole positive staining in epithelial regions (Figure 4). Our series also included a small number of patients who had received neoadjuvant chemo/radiotherapy before surgical resection. The range of pimonidazole-detectable hypoxia was similar in the tumours from patients who had received pre-operative treatment (range of epithelial hypoxia: 1–15%) to that of the larger group. However, it is not possible in this small series to comment on the impact of treatment on the extent of tumoral hypoxia and the potential of systemic therapy and/or radiation to modulate tumoral oxygenation status is beyond the scope of the current study but worth exploring in subsequent trials.

Clinical hypoxia research has been complicated by the diverse methods and variable thresholds of oxygenation considered, as reviewed extensively (Olive ; Bussink ; Wilson and Hay, 2011). In the context of PDAC, the only specific study of hypoxia utilised Eppendorf electrode probes, often considered a ‘gold standard' given its direct assessment of oxygenation status. Although the levels of tissue hypoxia observed in the present study are similar to those reported in nitroimidazole-based studies of other tumour types, they are much less severe than those reported using polarographic measurements made in a similar patient population (Brizel ; Koong ; Kaanders ; Nordsmark ; Chang ). This lack of concordance has been previously observed, and attributed to differences in technique resolution, necrosis artefact and tumour heterogeneity among others (Bussink ; Nordsmark ; Jankovic ; Nordsmark ). Equally important, the bioreductive activation and binding of nitroimidazole markers like pimonidazole requires very low levels of oxygen, and therefore identifies severely hypoxic tissue at 02 thresholds <1–5 mm Hg. However, moderate levels of hypoxia (pO2 of 10–20 mm Hg) might also be clinically relevant, given the upregulation of HIF1α and related pathways at that threshold (Hockel and Vaupel, 2001a, 2001b).

Although several hypoxia studies have been completed using nitroimidazoles in human tumours, few have attempted quantitative IHC image analysis as we have described. A small number have used quantitative immunofluorescence (IF) to study relationships amongst hypoxia, vasculature and proliferation (Kaanders ; Hoogsteen ; Evans ), but IHC studies to date have focused on manual, qualitative scoring schemes (Carnell ; Goethals ; Hoskin ). The potential advantages of IHC over IF include the ability to characterize hypoxia in the context of tissue morphology highly relevant for PDAC given its complex histology. Although most IHC based tissue biomarkers are evaluated manually, there is a significant literature describing inter-observer and inter-laboratory variability in the assessment of histology markers (Polley ). Therefore, there is interest in moving towards more quantitative analysis of tissue-based biomarkers with a goal of reducing variability and increasing throughput (Keller ; O'Hurley ). However, visual evaluation by an expert pathologist continues to be the standard approach for most IHC markers, and we noted a high concordance with automated analysis of hypoxia in the epithelial tumour component (Spearman's correlation r=0.91; P=0.0003; Supplementary Figure 2).

The important conclusions from this study are that the extent of severe hypoxia in PDAC can be assessed by preoperative administration of pimonidazole, and that it is highly heterogeneous. The analysis of multiple whole sections from the same tumour shows considerable intratumoral heterogeneity, which cautions against the use of single sections or biopsy samples to assess hypoxia in PDAC. The variance of hypoxia within tumours is less than that between individual patients, supporting the feasibility to stratify patients using this approach provided that the tumour is adequately sampled.