Digital quantitative assessment of PD-L1 using digital spatial profiling.

The assessment of programmed death 1 ligand 1 (PD-L1) expression by Immunohistochemistry (IHC) is the US Food and Drug Administration (FDA)-approved predictive marker to select responders to checkpoint blockade anti-PD-1/PD-L1 axis immunotherapies. Different PD-L1 immunohistochemistry (IHC) assays use different antibodies and different scoring methods in tumor cells and immune cells. Multiple studies have compared the performance of these assays with variable results. Here, we investigate an alternative method for assessment of PD-L1 using a new technology known as digital spatial profiling. We use a previously described standardization tissue microarray (TMA) to assess the accuracy of the method and compare digital spatial profiler (DSP) to each FDA-approved PD-L1 assays, one LDT assay and three quantitative fluorescence assays. The standardized cell line Index tissue microarray contains 10 isogenic cells lines in triplicates expressing various ranges of PD-L1. The dynamic range of PD-L1 digital counts was measured in the ten cell lines on the Index TMA using the GeoMx DSP assay and read on the nCounter platform. The digital method shows very high correlation with immunohistochemistry scored with quantitative software and with quantitative fluorescence. High correlation of PD-L1 digital DSP counts were seen between rows on the same Index TMA. Finally, experiments from two Index TMAs showed reproducibility of DSP counts were independent of variable slide storage time over a three-week period after antibody labeling but before collection of cleaved tags. In summary, DSP appears to have quantitative potential comparable to quantitative immunohistochemistry. It is possible that this technology could be used as a PD-L1 protein measurement system for companion diagnostic testing for immune therapy.

Introduction

Assessment of Programmed cell death ligand 1 (PD-L1) expression has been evaluated as a predictive diagnostic test to identify patients that may benefit from immune checkpoint blockade therapies such as anti-PD-L1 (atezolizumab and durvalumab) and anti-PD-1 (nivolumab and pembrolizumab) [1-4]. For each drug, a unique diagnostic assay has been developed. Four drug specific tests are approved by the US Food and Drug Administration (FDA) as either companion and/or complementary diagnostic assaysusing four different antibodies and four proprietary protocols to evaluate the PD-L1 expression by chromogenic immunohistochemistry (IHC). Each assay is read subjectively by a pathologist with a range of cutoffs defining expression positivity and cell specific expression either in tumor cell (TC) and/or immune cells (IC) [5-7]. Harmonization studies have been done to compare the assays within specific tumor types [8-15]. Despite strong agreement in TC PD-L1 IHC scoring among experienced pathologists, there is a poor concordance in IC PD-L1 IHC scoring or low PD-L1 scores [8, 12, 14–17]. This diversity of assays and scoring systems has led to confusion amongst both oncologists and pathologist.

A possible solution to this issue would be an objective, reproducible, and accurate measurement of PD-L1 in tissue sections. Digital Spatial Profiling (DSP), has potential to provide a platform technology for this sort of measurement [18]. In its current form, DSP is performed on the GeoMx platform and not amenable to usage as a companion diagnostic test. But the underlying technology could be reconfigured to make a closed system, spatially informed measurement of a biomarker. DSP technology is based on a UV cleavable DNA tag conjugated to an antibody. Spatial information is conveyed by selection of a defined area in which the UV laser directed. This region of interest can be a molecularly defined spatial compartment, for example epithelial cells defined by cytokeratin expression or T-lymphocytes defined by CD3. The region of interest serves as a guide for the UV laser, which cleaves off the DNA tags which are collected and subsequently counted using the Nanostring Barcode system [18]. Thus, there can be digital counts of antibody/antigen interactions in tumor cells, or in immune cells, as defined by a molecular compartmentalization tool. Although a different, CLIA lab friendly, platform would be required, we can envision this underlying technology being used as a digital method for cell specific measurement of PD-L1.

As a proof of concept and first step to determine the digital accuracy and reproducibility of DSP, we used a standardized cell line Index tissue microarray (TMA) with 10 isogenic cells lines expressing negative, weak, intermediate and high levels of PD-L1 protein in triplicate to objectively compare PD-L1 counts to quantitative PD-L1 expression with both chromogenic and fluorescent IHC. We also test various slide storage conditions and reproducibility over time. Future studies have been and will be done on actual tumor tissue, and are separately submitted for publication.

Materials and Methods

Standardized cell line index TMA and study design

Standardized cell line index TMA map with 3-fold redundancy consists of a 10 isogenic cell line array expressing negative, low, intermediate, and high PD-L1 protein from Horizon Dx (Figure 1a). The study plan was to stain three standardized cell line index TMA’s with a pre-commercial version of the Human Immune Oncology panel, a 45 unique antibody panel conjugated with ultraviolet light photocleavable linker to unique oligonucleotide readout on the nCounter platform for a quantitative measure (Supplemental Table 1). The TMA’s were also stained with a nuclear stain (SYTO 13) as a morphology marker for visualization on the GeoMx platform (Figure 1b) [18]. For first two experiments, slides were stained and collected in the same week, whereas, for the third experiment, the slide was stained, and each row was collected one week apart (Figure 2). For storage in between collections, slides were mounted with ~200ul aqueous mounting media (SouthernBiotech, 0100–01) to cover the cell pellets and stored in the dark at room temperature. To unmount the slides, the slides were place in TBST for 1hr until coverslips easily slid of the cell pellets.

Figure 1:

Index TMA layout for PD-L1 quantification.

(A) Schematic map along with ten standardized cell line codes and (B) representative image of the stained standardized TMA and 300um diameter ROIs as white circles.

Figure 2:

Overview of Experimental Plan.

Three standardized TMA were stained with Human IO panel and nuclear stain. First two experiments were stained and collected in the same week. Third experiment was stained and each row was collected one week apart.

Sample preparation for protein digital spatial profiling

TMA slides were stained using a protocol previously described [19]. Briefly, formalin-fixed paraffin-embedded (FFPE) 5µm sections were deparaffinized at 60°C for 30 minutes, then incubated thrice in CitriSolv (Decon Labs, Inc., 1601) for 5 min. Rehydration was performed with 100% ethanol twice for 10 min, 95% ethanol twice for 10 min and dH2O twice for 5 min. Antigen retrieval was performed with 1× citrate buffer pH 6.0 (Sigma-Aldrich, C9999) in a lid covered metal jar placed into a pressure cooker (BioSB, BSB 7008) on high pressure and temperature setting for 15 min. The metal jar was removed from the pressure cooker and cooled at room temperature for 25 min. Slides were washed with 1XTBST (Cell Signaling Technology, 9997) five times for 2 min and blocked with buffer W (NanoString, Seattle WA) for 1 h. Tissue sections were then stained with the Beta Human IO Protein panel and a nuclear stain as a morphology markers in Buffer W overnight at 4°C. Slides were washed with 1XTBST thrice for 10 min, fixed with 4% PFA (Thermo Scientific, 28906) for 30 min and washed again with 1XTBST twice for 5 min. DNA was counterstained with 500 nM SYTO 13 (NanoString, Seattle WA) diluted in 1XTBS (Cell Signaling Technology, 12498) for 15 min. Excess SYTO 13 was removed with one wash of 1XTBS for 5 min, and each slide was processed in an automated manner on the GeoMx DSP platform. SYTO 13 DNA dye was used as a visualization marker to identify regions of interest (ROIs). Per section, one ROI of 300um was selected and photocleaved oligonucleotides were released upon ultraviolet radiation exposure (Figure 1b). Released photocleaved oligonucleotides were collected via a microcapillary tube inspiration robotic system and transferred into a 96-microwell plate.

Hybridization assay

Hybridization of photocleaved oligonucleotides to optical fluorescent barcodes or GeoMx Hyb codes was performed using the nCounter Protein HybCode (NanoString) as directed by the manufacture’s protocol. In short, the photocleaved oligonucleotides were dried down and resuspended in 10μl of DEPC water. For the hybridization reaction 3μl of the reconstituted oligonucleotides, 5μl of ddH20 and 8μl of the appropriate GeoMx HybCode was combined and hybridized at 67°C for 18h with a 70°C heated lid on a thermocycler. Post-hybridization, 3μl per hybridization is pooled and processed using the nCounter PrepStation and Digital Analyzer according to the manufacturer’s instructions (NanoString, Seattle WA).

nCounter data analyses

Digital raw counts from barcodes corresponding to protein probes were first normalized with internal spike-in controls (ERCCs) to account for system variation, and then normalized to the endogenous control Histone H3. For signal to noise ratio, ERCC normalized digital counts were divided by geometric mean of all isotype controls (Mouse IgG1, Mouse IgG2a and Rabbit IgG).

Quantitative immunofluorescence

TMA slides were stained using a protocol described in our previous study [15]. Primary antibodies for PD-L1 E1L3N and SP142 were incubated overnight at 4 °C and clone SP263 was incubated for 20 minutes at 37 °C. Then, slides were incubated in rabbit EnVision reagent for 1 h at RT and Cy5-Tyramide was used to amplify target signal. Finally, TMAs were stained with 1:250 4,6-diamidino-2-phenylindole (DAPI) for 10 min at RT and mounted with Prolong Gold antifade mounting reagent. Image analysis was performed using AQUA method of QIF (NavigateBP), which generates a score by dividing the sum of target pixel intensities by the area of all nuclei compartment (DAPI staining). Three independent experiments were performed for each antibody.

Immunohistochemistry and digital analyses

Automated systems were used for different chromogenic PD-L1 assays, using protocols specified by corresponding manufacturer. Per the FDA-approved assays, 22C3 and 28–8 were stained with the Dako Autostainer Link 8 Instrument, and SP263 and SP142 were used on the Benchmark Ultra from Ventana Medical Systems, Inc. The LDT E1L3N assay was performed on a Leica Bond Autostainer. After chromogenic staining using DAB (3,3’-Diaminobenzidine), slides were scanned on the Aperio ScanScope XT platform. PD-L1 expression was quantified using the open source software QuPath [20]. An optimized algorithm was used for cell segmentation based on the size of the nucleus and cell expansion, and for DAB intensity quantification of PD-L1 expression based on DAB optical density (OD) mean. The settings were adjusted to avoid false positive detection. Results were shown as percentage of PD-L1+ cells or as OD of the chromogenic staining divided by mm2. Ten independent experiments were performed for each assay, including two slides per experiment.

Statistical analyses

Pearson’s correlation (R2) was computed for regression within GeoMx DSP or among GeoMx DSP, quantitative immunofluorescence (QIF) and IHC DAB results. All data sets were analyzed and plotted using GraphPad Prism v7.0 software for Windows (GraphPad Software, Inc., La Jolla, CA).

Results

To compare the DSP technology to conventional methods of assessment of protein on slides, we assessed PD-L1 expression compared to QIF and IHC DAB using standardized cell line index TMA. An average of six cores per cell line were evaluated for PD-L1 counts (three cores per cell line/two experiments) by DSP. This was compared to QIF quantification using three different PD-L1 antibodies on an average of 9 cores per cell line. Concordance of PD-L1 measurement between QIF scores obtained with three antibodies E1L3N, SP142, and SP263 and DSP counts was high, with coefficients >0.91 (Figure 3a–c). Five conventional IHC assays based on chromogenic visualization were performed and an average of 60 cores per cell line were measured for each assay (three cores per cell line/ten experiments/two slides per experiment). Four of the IHC assays are those currently FDA approved for companion or complementary testing for immunotherapy including assays named SP142, SP263, 22C3 and 28–8. A commonly used LDT assay based on the E1L3N antibody was also tested [15]. In each case, agreement with DSP counts for PD-L1 expression was excellent with coefficients ranging from 0.90 to 0.94 (Figure 3d–h).

Figure 3:

Correlation of Log2 transformed PD-L1 data from GeoMx DSP, QIF and IHC DAB.

Regression of PD-L1 protein expression by DSP with (A-C) QIF assay performed using 3 antibodies (E1L3N, SP142, and SP263) and (D-H) IHC DAB assay performed using 5 antibodies (E1L3N, SP142, SP263, 22C3 and 28–8). Each dot represents average of 2 TMAs (GeoMx DSP), 3 TMAs (QIF) and 20 TMAs (IHC DAB) with 3 pellet per cell clone in one TMA.

In order to assess the PD-L1 expression dynamic range between IHC and DSP, the Index TMA was used. The dynamic range of PD-L1 expression on the Index TMA was evaluated either as signal to noise ratio or PD-L1 counts. Signal to noise ratio in cell lines 8063 and 11198 was below 3 for all three rows in both the Index TMAs. Thus, confirming them as negative control. As expected, signal to noise ratio was highest in the highest expressing cell lines with one line averaging nearly 10,000 counts and a signal to noise ratio of nearly 100 (Figure 4 and Supp Figure 1).

Figure 4:

Reproducibility across rows and experiments stained and collected in the same week.

Bar graph show (A) signal to noise ratio and (B) PD-L1 counts for three rows from same experiment. Regression on Log2 scale for PD-L1 counts (C) between two rows from same experiment and (D) average PD-L1 counts of three rows between two experiments.

Correlation of PD-L1 counts in ten cell lines on the same index TMA were also assessed for reproducibility between rows and experiments on different days. The concordance was very high between rows (R2=0.97, figure 4c) and experiments on different days (average of three cores per cell line) from two Index TMAs (R2=0.99, figure 4d). Similar associations were observed for DSP counts of all targets from same cell line across two rows from same Index TMA and two experiments (average of three cores per cell line) from two Index TMAs, with coefficient 0.98 (Supp Figure 2a,b).

Finally, to test the robustness and reproducibility of GeoMx DSP count upon slide storage after antibody labeling but before cleavage and sipping of cleaved tags, we stained a standardized cell line index TMA with 45 markers and collected cleaved tags one week apart from each row for three consecutive weeks (Figure 2). For ten cell lines collected in three consecutive weeks, the distribution of un-normalized Histone H3 and average isotype control counts were consistent whereas the raw PD-L1 counts showed dynamic range of the PD-L1 expression in each week as expected (Figure 5a). Similar to previous finding described above, signal to noise ratio and PD-L1 counts were consistent and unchanged by longer storage prior to tag cleavage and collection (Figure 4b,c). Regression analysis of PD-L1 counts after progressively longer storage showed excellent concordance between three rows of ten cell lines collected on three consecutive weeks, with coefficients ranging from 0.97 to 0.98. Excellent correlations were observed for DSP counts of all targets from same cell line between rows of ten cell lines collected one week apart from each other from same Index TMA, with coefficient 0.97 (Supp Figure 2c,d).

Figure 5:

Reproducibility across rows collected one week apart from same staining.

(A) Distribution of PD-L1, Histone H3 and average IC raw counts from all cell pellet collected one week apart. Bar graph show (B) signal to noise ratio and (C) PD-L1 counts for three rows. Regression on Log2 scale for PD-L1 counts (D-E) between two rows collected on different weeks.

Discussion

Expression of PD-L1 protein by IHC has been utilized to predict response to the anti-PD-1/PD-L1 axis immunotherapies in more than 15 cancer types [21]. There are currently four FDA approved PD-L1 IHC assays as companion and/or complementary diagnostic tests, with different antibodies, cut points, scoring systems and cell specific expression associated with each test resulting in confusion in pathology labs for patients and their physicians. Since many studies have showed variable results with current assays [8-15], we have begun efforts to generate more quantitative and reproducible results. While DSP is not ready for the clinic, the technology behind the platform shows the kind of robustness and reproducibility that is typical for laboratory medicine pathology tests (generally better than IHC).

DSP instrument uses DMD mirrors to perform the spatially resolved profiling. This technology allows selection of regions of different shape or size, including regular (i.e. grid pattern) or irregular (i.e. areas defined by cell borders or proximity to specific tissue type) boundaries. The current platform on which the DSP assay is delivered has a number of limitations, representing the early state of the technology. First, profiling with DSP is expensive and it is inconvenient to fit the assay into an 8-hour day over current pathology IHC laboratory assays, but DSP has the ability to get high-multiplexed, quantitative and reproducible information. Second, DSP system is challenging for profiling every cell in the tissue at single cell resolution. Although imaging based methods have ability to generate multiplexed information on every cell in the tissue, DSP offers high-plex profiling capabilities in ROIs and tissue compartments. Therefore, imaging based methods and DSP are complementary to each other. Third, DSP technology aligns well with Laser Capture Microdissection technologies. However, DSP offers advantage of being high-throughput and automated. Lastly, DSP assay also requires manual or molecular selection of compartments for analysis which would need to be more accurately defined prior to clinical translation. The results of this study allow protocol modifications that may facilitate the development of this technology by increasing operator convenience. A major limitation of this study includes the usage of isogenic cell line for assessing PD-L1 sensitivity and reproducibility instead of tissue. Future studies will assess tissue where expression of PD-L1 may be different or possess post-translational modifications in tumors that might not be reflected in these cell lines. However, in testing analytic validity, this approach has advantages including elimination of tissue heterogeneity from the experiment in the comparison of platform reproducibility.

In summary, we were able to objectively quantity the PD-L1 expression in this standardized cell line Index TMA and report the DSP assay to be highly reproducible and independent of slide storage. While not yet “packaged” for the pathology lab, we believe the underlying DSP technology is truly digital and has the potential to accurately measure protein expression on a slide in a “closed box” system. We look forward to next generation DSP devices that may have the capacity to automate, standardize and quantify protein within spatially defined compartments in an objective manner.