Heterogeneity Mapping of Protein Expression in Tumors using Quantitative Immunofluorescence

Summary

Here we describe a method to quantify molecular heterogeneity in histological sections of tumor material using quantitative immunofluorescence, image analysis, and a statistical measure of heterogeneity. The method is intended for use in clinical biomarker development and analysis.

Abstract

Morphologic heterogeneity within an individual tumor is well-recognized by histopathologists in surgical practice. While this often takes the form of areas of distinct differentiation into recognized histological subtypes, or different pathological grade, often there are more subtle differences in phenotype which defy accurate classification (Figure 1). Ultimately, since morphology is dictated by the underlying molecular phenotype, areas with visible differences are likely to be accompanied by differences in the expression of proteins which orchestrate cellular function and behavior, and therefore, appearance. The significance of visible and invisible (molecular) heterogeneity for prognosis is unknown, but recent evidence suggests that, at least at the genetic level, heterogeneity exists in the primary tumor^1,2, and some of these sub-clones give rise to metastatic (and therefore lethal) disease.

Moreover, some proteins are measured as biomarkers because they are the targets of therapy (for instance ER and HER2 for tamoxifen and trastuzumab (Herceptin), respectively). If these proteins show variable expression within a tumor then therapeutic responses may also be variable. The widely used histopathologic scoring schemes for immunohistochemistry either ignore, or numerically homogenize the quantification of protein expression. Similarly, in destructive techniques, where the tumor samples are homogenized (such as gene expression profiling), quantitative information can be elucidated, but spatial information is lost. Genetic heterogeneity mapping approaches in pancreatic cancer have relied either on generation of a single cell suspension³, or on macrodissection⁴. A recent study has used quantum dots in order to map morphologic and molecular heterogeneity in prostate cancer tissue⁵, providing proof of principle that morphology and molecular mapping is feasible, but falling short of quantifying the heterogeneity. Since immunohistochemistry is, at best, only semi-quantitative and subject to intra- and inter-observer bias, more sensitive and quantitative methodologies are required in order to accurately map and quantify tissue heterogeneity in situ.

We have developed and applied an experimental and statistical methodology in order to systematically quantify the heterogeneity of protein expression in whole tissue sections of tumors, based on the Automated QUantitative Analysis (AQUA) system⁶. Tissue sections are labeled with specific antibodies directed against cytokeratins and targets of interest, coupled to fluorophore-labeled secondary antibodies. Slides are imaged using a whole-slide fluorescence scanner. Images are subdivided into hundreds to thousands of tiles, and each tile is then assigned an AQUA score which is a measure of protein concentration within the epithelial (tumor) component of the tissue. Heatmaps are generated to represent tissue expression of the proteins and a heterogeneity score assigned, using a statistical measure of heterogeneity originally used in ecology, based on the Simpson's biodiversity index⁷.

To date there have been no attempts to systematically map and quantify this variability in tandem with protein expression, in histological preparations. Here, we illustrate the first use of the method applied to ER and HER2 biomarker expression in ovarian cancer. Using this method paves the way for analyzing heterogeneity as an independent variable in studies of biomarker expression in translational studies, in order to establish the significance of heterogeneity in prognosis and prediction of responses to therapy.

Protocol

1. Tissue preparation

1. Microtomy
   1. Section formalin fixed paraffin embedded tumor blocks at 4 μm thickness using a rotary microtome.
   2. Place the sections onto positively charged microscopic slides.
   3. Incubate slides in a 37 °C oven overnight.
2. Storage
   1. Store sections in a -18 °C to -25 °C freezer, in order to prevent loss of antigenicity.

2. Immunofluorescence of tissue sections

1. Dewaxing and Rehydration
   1. Dewax slides in xylene twice for 5 min.
   2. Rehydrate slides in 99%, 99%, 80%, 50% of ethanol and running tap water for 2 min of each in order.
2. Antigen Retrieval
   1. Perform the heat induced antigen retrieval, using 0.15 mM Sodium Citrate, pH 6.0 buffer in a pressure cooker for 5 min in the microwave.
   2. Cool down slides for 20 min.
   3. Rinse slides in 0.05% PBST for 5 min on the rocker.
3. Blocking procedure
   1. Treat sections in 3% hydrogen peroxide for 10 min.
   2. Rinse slides in 0.05% PBST for 5 min on the rocker.
   3. Add the sections to Sequenza Immunostaining rack.
   4. Treat sections in Serum free protein block for 10 min.
4. Primary antibody incubation
   1. Incubate Primary antibody diluted to optimal dilution in Dako antibody diluent for 1 hour at room temperature.
   2. Rinse in 0.05% PBST three times for 5 min each.
5. Epithelial mask (Cytokeratin) incubation
   1. Incubate mouse anti-pan cytokeratin, diluted 1:50 in Dako antibody diluent overnight at 4 °C.
6. Epithelial mask visualization
   1. Prepare a 1:25 dilution of the goat anti-mouse Alexa555 secondary antibody in the pre-diluted Envision goat-rabbit HRP antibody solution. Incubate slides in the dark for 1.5 hrs at room temperature.
   2. Rinse in 0.05% PBST three times for 5min each.
7. Target visualization
   1. Transfer the slides from the Sequenza to the humidity chamber.
   2. Combine Target signal amplification diluent (HistoRx tub E) and the Cy5 Tyramide (HistoRx tube F) at 1:50 concentration. Vortex to mix thoroughly. Incubate slides in the dark for 10 min at room temperature.
   3. Rinse in 0.05% PBST three times for 5 min each.
   4. Dehydrate slides in 80% ethanol for 1 min.
   5. Air dry slides in the dark.
8. Counterstaining and coverslipping
   1. Apply DAPI mounting medium on the coverslips and place the coverslips over the tissue sections.
   2. Dry slides overnight in the dark.
   3. After the slides have dried, secure with nail polish to ensure long term preservation and keep in a 4 °C refrigerator.

3. Image capture using whole section scanning

1. Load up to five slides in the slide cassette of the Aperio ScanScope FL slide scanner.
2. For each slide, select the general region to image using the interface in the ScanScope Console. Select the region to encompass all the tissue on the slide, regardless of staining pattern or other observable aspects of the sample.
3. Using the ScanScope Console, for each channel (DAPI, Cy3 and Cy5) determine an optimal exposure as follows:
   1. Select a region of tissue to interrogate - this area should, ideally, be in the tumor region of tissue to provide the best representation.
   2. Using the ScanScope Console autoexposure feature, evaluate the exposure time suggested for each channel, along with image focus.
   3. Repeat over up to 3-5 regions of each sample to confirm stability of the value.
4. Select an additional region (as shown by a blue diamond on the ScanScope console image) away from the tissue, yet still within the coverslipped region of the slide, to define a background area where flat field correction images will be acquired for each filter.
5. Review these images prior to image acquisition. If the images appear to have high background or other image artifacts, the image region should be moved and new images acquired.
6. The acquired digital slide images are uploaded automatically into the Aperio Spectrum database.

4. AQUA automated image analysis

1. View the digital whole slide images in the Spectrum database and annotate tumor areas of interest in the context of the entire tissue and staining pattern(s).
2. Select the slide to analyze using AQUAnalysis from the digital slide list in the Spectrum database by double clicking the thumbnail image (alternatively, select the check box next to the image and then select 'View Images' from the menu at the top of the list). This automatically opens the ImageScope software and color merged image of the tissue sample.
3. Use this software to navigate the image using the tools provided.
4. Using the accompanying H&E (either the physical slide or an annotated image), annotate the region of interest on the fluorescent image. Multiple regions of interest, individual circled areas, can be selected on a single sample. There is also a 'negative' tool which is used to exclude areas within a selected region (i.e. damaged tissue, areas of poor histology, non-tumor areas, debris, etc.).
5. Save the annotation layer(s) on the image.
6. Return to the Spectrum main menu and select the check box next to the image that was just annotated. Select 'Analyze' from the menu strip across the top of the list.
7. In the Analysis window which will come up, select "HistoRx AQUA Analysis" from the pull-down menu.
8. If there are multiple layers of annotations, use the check boxes to select the appropriate layer(s) for analysis.
9. Press the 'Analyze' button when ready to begin.
10. At this point, a minimized console window will appear in the Windows taskbar. This will display the progress of the transfer of images from the Spectrum database to the local PC (the 'pull' operation).
11. Once the data is transferred, AQUAnalysis will launch. The user can then use all of the internal functions of the software to view, navigate and analyze images.
12. The annotated region of the image from the Spectrum database is broken into 512x512 pixel image tiles - each of which will be scored individually.
13. For this study, an unsupervised AQUA scoring algorithm, based on data clustering, was used⁸. In this algorithm, AQUA scores are generated as follows:
    1. The pan-cytokeratin image signal is thresholded above background and those pixels are used to define a tumor 'mask' which is then used for subsequent calculations. The high -expression pan-cytokeratin pixels also define the cytoplasmic/non-nuclear regions of the tumor cells.
    2. Pixels with high signal in the DAPI channel that are within the tumor mask are identified as nuclear in nature.
    3. The clustering algorithm also examines pixels which have low intensity for either DAPI staining or cytokeratin expression and attempts to partially attribute them to either the nuclear or cytoplasmic compartments, or background.
    4. Once all pixels are assigned to one of the tumor compartments or background, the intensity of signal from the target protein (ER or HER2) is then calculated within each of the compartments and normalized to the size of the compartments, thus producing AQUA scores.
    5. Images which did not have sufficient tumor representation were not scored (as defined by those with less than 5% of pixels representing tumor). These regions also produce a 'Final Result' value called 'Fail' and can thus be removed in subsequent statistical analyses.
    6. Additionally, if upon review, an operator identified fields which were not suitable for scoring due to sample or imaging artifacts (e.g. folded tissue, debris), those fields were redacted from scoring and produced a 'Final Result' called 'Fail'.
14. The results of AQUA scoring are tables of scores associated with each 512x512 pixel tile in the region of interest within a whole tissue section sample. These files, in .CSV format, can then be used for subsequent statistical analysis.

5. Statistical Analysis: All statistical analyses are performed in SPSS (IBM)

1. Where possible, for large operations, generate SPSS syntax code files (.sps) to perform data manipulations and analysis steps. Example syntax code is provided as supplementary files.
2. Read data files into SPSS. Remove unused variables in the raw data (e.g. instrument specific parameters, AQUA scores from non-utilized compartments).
3. Where regions have a 'Final Result' of 'Fail' (see above), mark these values as SYSMIS, which removes them from the numerical analyses.
4. Aggregate all data for all slides for each marker (ER or HER2) into a single master data set in SPSS format. Use this master file for all subsequent calculations.
5. Using an SPSS syntax file (one for each marker: ER or HER2 (see supplementary syntax files) data was further validated and Simpson's Indices of diversity calculated as follows:
   1. Generate a master analytical file, in order to prevent modification of the master data set previously generated.
   2. Validate the master data set against the raw data and images generated.
   3. Produce summary statistics for AQUA scores for each sample (count of fields, mean score, median score, minimum score, maximum score, and standard deviation of scores for sample).
   4. Check a random sample of results against the original raw data output files and against the raw image data scored.
6. To generate the heat map images (shown in Figures 2 and 3):
   1. 'Bin' the AQUA score data (using the 'Visual Binning' function in SPSS) into eight individual levels. For ER, nuclear scores were used, for HER2, cytoplasmic scores were used.
   2. Derive bins such that each bin represents an equal percentage of available cases from the entire data set.
   3. Use the SPSS function 'Split File' so that all subsequent output would be on a 'per slide/sample' level.
   4. Plot each region (corresponding to a 512x512 tile) by its x,y coordinates and color with a value corresponding to the bin into which it was allocated (see figures 2 and 3).
7. To generate heterogeneity scores, code was written to generate a Simpson's index of diversity using AQUA scores.
   1. The Simpson's index of diversity, as used here, is defined as:
      
      Where N is the total number of scores/fields used for a given sample and n is the number of scores/fields in a given group of standardized scores (produced as described below).
8. Use the master analytical file described above for data source. In the examples shown, all ER calculations used the nuclear AQUA score and HER2 calculations used the cytoplasmic AQUA score. Perform the calculations below for each sample/slide.
9. Since fluorescence data is subject to variance instability, such that error magnitude propagates with intensity, apply logarithmic (base 2) normalization to all AQUA scores.
10. Further transform the normalized (log transformed) data into z-(standard) scores.
    1. The Z-score transformation is calculated as follows: Z = (AQUA score - average of all scores for a sample)/standard deviation of scores for the sample. This translates the distribution of scores into a distribution of deviations around the central value of zero deviation.
11. Bin these Z-score values into six groups (<-2, -2 - -1, -1 - 0, 0 - 1, 1 - 2, >2).
12. Calculate the number of values allocated into each bin and sum for each bin. Use this value, along with the total number of fields, to calculate a value for each of the six bins used for the index.
13. Sum these values and subtract from one to produce the index of diversity.
14. The diversity index provides an indicator of the spread of scores around the mean for a particular sample. Thus, higher indices indicate a broader spread of scores, or increased heterogeneity of expression. Conversely, lower indices, indicate less variation in expression measurement and homogeneous expression in the sample. This is illustrated in Figures 2 and 3.

6. Representative Results:

Targeted therapy with relatively low toxicity agents, such as tamoxifen and trastuzumab, are not standard of care for ovarian cancer patients. There is good evidence that platinum-taxane first-line chemotherapy is superior to other chemotherapy regimens for ovarian cancer^9-13, and while 70-80% of patients initially respond, most will relapse and die of their disease. Measurement of biomarkers which might predict response to alternative therapies would be useful in selecting individualized therapy for each patient, and since chemotherapy exerts selection pressure on a tumor, surviving tumor cells may possess different characteristics and represent a subpopulation of the tumor before therapy; quantifying this change might therefore be useful. We applied our heterogeneity mapping approach to ER and HER2 expression in ovarian cancer specimens before and after chemotherapy, in order to measure quantitative changes in expression, and assess the relationship between biomarker expression before and after therapy. Both ER and HER2 demonstrate marked heterogeneity within individual tumors, and in some tumors, the Simpson's index changes after therapy from relatively high heterogeneity to low heterogeneity, represented by decreasing Simpson's index scores (see Figures 2 and 3). This supports the notion that clonal selection occurs in response to cytotoxic chemotherapy; more importantly, this may be exploited using targeted therapy against the residual population, in order to better personalize therapy for cancer patients.

Figure 1. Hematoxylin and eosin stained photomicrographs of different areas of the same ovarian cancer showing variable morphology. (A) x40 photomicrograph shows two adjacent areas of tumor with varying morphology. High power views (x200) reveal cells with cytoplasmic clearing and a solid pattern of growth in the upper part of the tumor (B), while the lower part of the tissue (C) has more homogeneous eosinophilic cytoplasm and a papillary growth pattern. This tumor would, however, be classified as serous papillary histological subtype for the purposes of further management.

Figure 2. Heterogeneity heatmaps of ER protein expression in cytokeratin-positive ovarian cancer tissue, before (upper panel) and after (lower panel) therapy.

Figure 3. Heterogeneity heatmaps of HER2 protein expression in cytokeratin-positive ovarian cancer tissue, before (upper panel) and after (lower panel) therapy.

Discussion

The method described herein permits quantification of molecular heterogeneity in standard formalin-fixed, paraffin-embedded histologic sections of tumor material. The method permits a value to be assigned to the degree of heterogeneity in a tissue section, so that this may then be taken into account as a variable in biomarker development and analysis. While in general it is preferable that tissue biomarkers are homogeneous with respect to expression, so that the assay is less susceptible to sampling error or bias, it might under some circumstances be useful to quantify heterogeneity as a parameter. For instance, as shown in the illustrative example for ER and HER2, common biomarkers show considerable heterogeneity of expression and it is currently unknown whether this represents an independent prognostic or predictive factor with respect to clinical outcome or response to therapy, respectively. Likewise, as illustrated, therapy itself might dynamically alter the constituent populations of cells, and therefore measurement of target enrichment might be useful in guiding treatment decisions.

However, the technique is limited by the same factors which limit traditional histopathological or biomarker (such as immunohistochemical) analyses. The result is dependent on the type, size, and quality of the tissue being analyzed, and therefore a single tissue section may not be representative of the entire tumor. Indeed, the Simpson's index may be unduly biased by size of the tissue (number of frames captured and AQUA scores measured). The immunogenicity of the epitopes being measured may be subject to uncontrollable pre-analytical factors which artifactually alter their expression across the tissue section (such as cold ischemia time, or the length of fixation, and edge artifact). This is particularly a problem for phospho-epitopes, which are notoriously labile^{14, 15}. Therefore the technique might be better suited to resection specimens rather than small biopsies, and performed on multiple sections rather than just one. Ultimately, 3D imaging techniques may offer an opportunity to better quantify this parameter.

The technique as presented may also be limited by biological factors, such as variability in expression of the cytokeratin masking epitope. This might be of particular concern in those cancers, including ovarian and breast cancer, which undergo epithelial-to-mesenchymal transition (EMT) or are enriched for non-epithelial components or stem cells¹⁶, as has recently been shown to occur in chemotherapy-treated ovarian cancer in vitro¹⁷, and in breast cancer patients in response to therapy¹⁸. This limitation may be overcome using alternative masking antibodies (or cocktails of antibodies), such as cytokeratin plus vimentin, which is upregulated in EMT¹⁶. In addition, since other components of the tissue (the microenvironment), such as fibroblasts or vascular endothelial cells are also increasingly being targeted by biological therapies, the technique may also be expanded to to score these components, using specific masks for compartments of interest, such as PECAM/ CD31 to stain for vascular endothelial cells.

Although the adaptation of the Simpson's index has been used as a measure of heterogeneity in the current protocol due to its simplicity, other measures of heterogeneity could be used, such as simple measurements of central tendency (mean, variance, median, etc). Additionally, the Simpson's index calculations could be modified to better suit a particular test population. The use of Z-score transformations allows the scoring to be more applicable for multiple markers since the heterogeneity is based on deviation around the mean for a data set. However, the overall range of scoring can be limited in this type of transformation. Some designs may be better suited to simply bin the actual AQUA scores for all samples of a particular marker set. The selection of the number of bins and cutoffs is also a limitation in the range of values which can result and may be optimized for alternative experimental designs.

We have used ER and HER2 in the current study as candidate biomarkers, since they are already used in the clinic, help answer relevant clinical questions with respect to the efficacy of targeted therapies, and are known to exhibit considerable heterogeneity and change in primary and distant disease¹⁹. However, the optimal marker of heterogeneity remains to be determined. Other possibilities might include cytogenetic markers of clonality, such as those previously used in interphase FISH studies²⁰. The approach requires further validation in large, well annotated clinical cohorts or clinical trials in order to further establish the relevance of this parameter in cancer biology.

Materials

Name	Company	Catalog Number	Comments
Name of the reagent	Company	Catalogue number	Comments
ScanScope FL	Aperio	ScanScope FL	Used as equipped from vendor
Spectrum	Aperio	Spectrum	Configuration and catalog numbers will vary on how local institutions set up their database and IT structure.
AQUAnalysis, in conjunction with AQUAjustment	HistoRx	V2.3	The AQUAjustment package allows AQUAnalysis to use images stored in Spectrum.
HER2	Dako	A0485	1 in 400 dilution
mouse anti-pan cytokeratin	Dako	M3515	1 in 50 dilution
Antibody diluent	Dako	S0809
Envision goat-rabbit HRP	Dako	K4003
Protein block	Dako	X0909
Goat anti-mouse Alexa555	Invitrogen	A21422
DAPI mounting medium	Invitrogen	P36931
Cy5 Tyramide	HistoRx	AQ-EMR1-00001
Sequenza Immunostaining Center	Thermo Scientific Shandon	73300001	Used here for bench-top immunofluorescence staining