Next Article in Journal
Pediatric Hepatocellular Adenomas: What Is Known and What Is New?
Previous Article in Journal
The Impact of the COVID-19 Pandemic on Cancer Mortality in Pennsylvania: A Retrospective Study with Geospatial Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes

1
Labcorp Oncology, Durham, NC 27560, USA
2
Division of Gynecologic Oncology, Duke Cancer Institute, Durham, NC 27710, USA
3
Division of Medicine, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
4
Fred Hutchinson Cancer Center, Seattle, WA 98109, USA
5
Mantech International, Virginia Beach, VA 23452, USA
6
OmniSeq, Inc., Buffalo, NY 14203, USA
7
Wake Forest Baptist Comprehensive Cancer Center, Winston-Salem, NC 27157, USA
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(19), 4789; https://doi.org/10.3390/cancers15194789
Submission received: 11 August 2023 / Revised: 18 September 2023 / Accepted: 25 September 2023 / Published: 29 September 2023
(This article belongs to the Section Cancer Therapy)

Abstract

:

Simple Summary

We compared RNA next-generation sequencing (RNA-seq) to standard immunohistochemistry (IHC) for PD-L1 expression measurement and associations with pembrolizumab immunotherapy outcomes in NSCLC patient tumors. RNA-seq and IHC PD-L1 score interpretation agreed for 80% of patients, and an RNA-seq “high” cutoff that accurately separated IHC high versus low or negative expression was identified. However, RNA-seq could not discern PD-L1 IHC from negative expression due to the limited sensitivity of IHC as a reference test. High PD-L1 expression by RNA-seq alone and in combination with IHC high or low status was associated with better pembrolizumab outcomes in NSCLC patients than IHC alone.

Abstract

Programmed cell death ligand (PD-L1) expression by immunohistochemistry (IHC) lacks sensitivity for pembrolizumab immunotherapy selection in non-small cell lung cancer (NSCLC), particularly for tumors with low expression. We retrospectively evaluated transcriptomic PD-L1 by mRNA next-generation sequencing (RNA-seq). In an unselected NSCLC patient cohort (n = 3168) tested during standard care (2017–2021), PD-L1 IHC and RNA-seq demonstrated moderate concordance, with 80% agreement overall. Most discordant cases were either low or negative for PD-L1 expression by IHC but high by RNA-seq. RNA-seq accurately discriminated PD-L1 IHC high from low tumors by receiver operator curve (ROC) analysis but could not distinguish PD-L1 IHC low from negative tumors. In a separate pembrolizumab monotherapy cohort (n = 102), NSCLC tumors classified as PD-L1 high versus not high by RNA-seq had significantly improved response, progression-free survival, and overall survival as an individual measure and in combination with IHC high or low status. PD-L1 IHC status (high or low) trended toward but had no significant associations with improved outcomes. Conventional PD-L1 IHC testing has inherent limitations, making it an imperfect reference standard for evaluating novel testing technologies. RNA-seq offers an objective PD-L1 measure that could represent a complementary method to IHC to improve NSCLC patient selection for immunotherapy.

1. Introduction

Immune checkpoint inhibitors have improved outcomes over chemotherapy for patients with non-small cell lung cancer (NSCLC) [1,2]. Checkpoint inhibitors are a class of immunotherapy drugs designed to enhance the body’s natural immune response against cancer cells. They work by targeting specific proteins, known as immune checkpoints, that regulate the immune system’s response to prevent it from attacking normal cells. When these checkpoints are inhibited, it allows the immune system to recognize and attack cancer cells more effectively. One of the most well-known immune checkpoints is the programmed cell death protein 1 (PD-1) and its ligand, programmed death-ligand 1 (PD-L1). When PD-1 on the surface of T-cells binds with PD-L1 on the surface of cancer cells, it sends a signal that suppresses the immune response, allowing cancer cells to evade detection and destruction by the immune system and leading to PD-L1 protein expression [3,4,5]. Currently, eligible NSCLC patients receive pembrolizumab anti-PD-1-containing regimens as a new frontline treatment standard based on the expression of PD-L1 [6].
The immunohistochemistry (IHC) 22C3 FDA-approved companion diagnostic assay for evaluation of PD-L1 protein expression in NSCLC tissue biopsy specimens is indicated for selecting patients for pembrolizumab immunotherapy [7]. While IHC 22C3 measurement technology is intended to capture the full range of possible PD-L1 expression across tumors, the test lacks high sensitivity, misclassifies tumors, and may negatively impact therapeutic efficacy [8]. Problems with the IHC 22C3 assay arise from challenges to pathologist interpretation from limited NSCLC core biopsy specimens; PD-L1 expression can be highly heterogeneous within a given tumor, between primary and metastatic sites, and for different specimen types [9]. This introduces complexity to calculating an accurate tumor proportion score (TPS), which is the percentage of tumor cells subjectively showing complete or partial linear membranous staining, in accordance with FDA-approved pathologist scoring guidelines for IHC 22C3 testing [10,11,12,13,14].
The IHC 22C3 assay also renders an incomplete picture of PD-L1 in NSCLC because it excludes counting expression on immune cells, which are regulated by a different mechanism with an independent role in immunotherapy response. High tumoral expression of PD-L1 reflects epigenetic dysregulation of the PD-L1 gene and is associated with poor immune cell infiltration, desmoplastic stroma, and mesenchymal histologic features. In contrast, high immune cell expression of PD-L1 is more common and reflects IFN-γ-induced adaptive regulation, accompanied by increased tumor-infiltrating lymphocytes and effector T-cells [15]. Capturing PD-L1 expression on both tumor and immune cells in the landmark KEYNOTE-001 pembrolizumab trial using IHC was initially evaluated based on the 22C3 assay prototype used to enroll patients. The trial showed PD-L1 IHC 22C3 expression in at least 50% of tumor cells was correlated with improved response, which led to the initial FDA approval of pembrolizumab monotherapy in the NSCLC second-line treatment setting. Post-trial, the PD-L1 IHC 22C3 assay was modified to exclude immune cell staining from 22C3 companion diagnostic scoring validation. PD-L1 TPS ≥ 50% was selected as the IHC 22C3 companion diagnostic cutoff from among highly similar results for four scoring methods because it was the simplest approach to implement in clinical practice, with the investigators concluding PD-L1 expression in NSCLC is tumor-specific [16,17].
Since KEYNOTE-001, however, the 50% TPS cutoff for high PD-L1 expression by IHC 22C3 has proved problematic. Most notably, the KEYNOTE-042 trial initially showed NSCLC patients with any positive PD-L1 IHC staining (≥1% TPS) had a survival advantage for frontline pembrolizumab over chemotherapy (16.7 vs. 12.1 months), followed by the FDA’s lowering of the cutoff from ≥50% TPS to ≥1% TPS on the label indication. However, a later pre-specified subgroup analysis showed KEYNOTE-042 patients with low PD-L1 IHC expression had no significant survival benefit compared with patients who received chemotherapy (13.4 vs. 12.1 months) [18,19,20,21]. As a result, professional practice guidelines now recommend NSCLC patients with high (≥50% TPS) PD-L1 IHC expression receive frontline pembrolizumab monotherapy, while patients with low (≥1% TPS) or negative (<1% TPS) PD-L1 IHC expression receive combination pembrolizumab plus chemotherapy, in contradiction to the FDA label indication [22].
The PD-L1 IHC SP142 companion diagnostic assay for atezolizumab in NSCLC, which does include tumor and immune cell staining [23,24], has also struggled to demonstrate concordance and predictive value. A comparison of the Ventana SP142, Ventana SP263, and Dako 22C3 IHC-validated assays showed the percentage of PD-L1-positive cells was consistently higher among immune cells than tumor cells across all three tests, but assay concordance was lower for immune cell vs. tumoral expression measurement [25]. Other studies have also shown both greater variability, subjectivity, and poor inter-pathologist agreement in scoring PD-L1 expression on immune cells compared with tumor cells when using the percentage of tumor with immune cell infiltrate [14,26,27]. Despite these inconsistencies, PD-L1 expression on either tumor or immune cells has been associated with an improved response to checkpoint blockade. In the POPLAR trial, second-line, NSCLC patients with PD-L1 IHC SP142 expression ≥ 1% in either tumor or immune cells had significantly better overall survival when treated with atezolizumab monotherapy versus docetaxel (median 15.5 months versus 9.2 months, p = 0.005) [28]. Problematically, however, the later IMpower110 NSCLC frontline atezolizumab trial compared IHC SP142 to IHC 22C3 in pre-specified exploratory analyses and showed the SP142 assay was less sensitive for staining both tumor and immune cell PD-L1 [29].
Continued advancement in precision oncology clinical care requires a re-evaluation of testing methods for established biomarkers, especially when accumulating data demonstrates the imperfect nature of existing approaches like IHC. Compared with IHC, next-generation sequencing (NGS) objectively measures PD-L1 mRNA (RNA-seq) by digitally counting target molecules, enabling more sensitive and precise transcriptome quantification and a continuum of measurements across a large dynamic range of expression. RNA-seq methods have also previously been shown to be positively correlated with IHC testing for PD-L1 [25,30,31,32,33].
IHC and RNA-seq are distinct techniques with their own advantages and limitations. RNA-seq is more sensitive; however, standard NGS does not discern the detection of PD-L1 mRNA in tumor versus stroma. RNA-seq does not quantitatively assign PD-L1 expression levels to specific cell types, although the clinical value of making this distinction using current IHC techniques also appears to be variable, with NSCLC tumors that harbor one and not the other potentially representing different patient subgroups [15,25,34]. In addition to these divergent detection methods, underlying biological mechanisms can also result in differences in PD-L1 expression, including mRNA stability, which can affect the availability of translation into proteins, translational regulation factors, post-translational modifications, and protein degradation. Understanding these mechanisms is also crucial for the accurate interpretation of PD-L1 expression data.
In the absence of a single, perfect technology and with the growing use of NGS in NSCLC patient standard care, we posit there may be a role for measuring PD-L1 expression by both IHC and RNA-seq. Here, we sought to address unanswered questions about the sensitivity, potential clinical utility, and real-world application of RNA-seq to improve the selection of NSCLC patients for immunotherapy.

2. Methods

2.1. Patient Cohorts

Approval for this study was obtained from an independent institutional review board (IRB), the Western Copernicus Group (www.wcgclinical.com (accessed on 14 September 2023)), under protocol #1340120. NSCLC patients treated between 2017 and 2021 at over 300 oncology practices in the United States whose FFPE tumor tissue biopsy samples underwent testing in a commercial laboratory by clinically validated comprehensive immune profiling (CIP) were considered for inclusion. Cases with non-NSCLC histology failed RNA-seq testing, or failed/inconclusive IHC results were excluded. The final set of NSCLC patients (n = 3270) was divided into 2 groups—an unselected cohort with unknown treatment status (n = 3168) to evaluate assay concordance and explore alternative RNA-seq cut-offs, and a cohort of pembrolizumab monotherapy-treated patients (n = 102) to evaluate outcomes associations with results for each assay. The pembrolizumab outcomes cohort included patients with stage IIIB/IV disease who were confirmed negative for ALK and EGFR mutations, received at least one dose of pembrolizumab monotherapy post-CIP test and had at least 90 days of follow-up (excluding death) post-first dose.

2.2. PD-L1 Expression Testing and Clinical Reporting Cutoffs for IHC and RNA-Seq

PD-L1 expression by IHC and RNA-seq was assessed for the same tissue from each patient as part of a larger comprehensive immune profiling (CIP) assay ordered alongside comprehensive genomic profiling (CGP) as part of standard patient care. For CIP, formalin-fixed paraffin-embedded (FFPE) tissue specimens (3–5 unstained positively charged slides cut at 5 µm), plus one hematoxylin and eosin-stained slide, were submitted along with the corresponding pathology report for each order. Histology confirmation, tumor enrichment, and estimation of the minimum required 20% tumor purity for RNA-seq were determined through microscopic review and cell selection by a board-certified pathologist. PD-L1 IHC protein expression was measured with the FDA-approved Dako PD-L1 IHC 22C3 pharmDx assay (Agilent, Santa Clara, CA) on Autostainer Link 48 (Agilent, Santa Clara, CA) and scored as the percentage of viable tumor cells showing % membrane staining at any intensity as a tumor proportion score (TPS) following FDA scoring guidelines [7]. TPS scores were interpreted as ≥50% (high), 1–49% (low), or <1% (negative) per NCCN guideline recommendations for pembrolizumab regimen treatment selection [22].
PD-L1 by RNA-seq was semi-quantitatively measured by analytically validated amplicon-based targeted RNA next-generation sequencing as part of a 64-transcript digital gene expression panel relating to the anticancer immune response as previously described [35]. RNA was extracted from FFPE using the truXTRAC FFPE extraction kit (Covaris, Inc., Woburn, MA, USA), with some modifications to the instructions by the manufacturer. After purification, RNA was dissolved in 50 µL of water, and the yield was measured through the Quant-iT RNA HS assay (Thermo Fisher Scientific, Waltham, MA, USA), as per the manufacturer’s recommendation. For library preparation, the pre-defined titer of 10 ng RNA was referred to as the acceptance criterion. Torrent Suite’s plugin immuneResponseRNA (v5.2.0.0) was used for the absolute reading of the RNA sequence. The RNA expression of 64 different genes was measured. Transcript abundance was normalized to internal housekeeping gene profiles and ranked (0–100 percentile) in a standardized manner to a reference population of 735 tumors spanning 35 histologies as previously described [35]. In clinical testing, a percentile rank ≥75 is interpreted and reported as the cutoff for “high” versus “not high” expression for each of the 64 RNA-seq genes in the CIP panel, including PD-L1. The clinical cutoff for reporting was determined based on the rank-score distributions observed in assay validation. RNA-seq expression results for immune genes that are drug targets in clinical trials are currently reported for research use only.

2.3. Clinical Data and Pembrolizumab Outcomes Measures

For all patients, clinicopathological data included age, sex, histology (squamous or non-squamous), specimen type (primary or metastatic), and anatomic site of metastasis (adrenal gland, bone, brain, liver, lymph node, pleura, soft tissue, and others). Additional clinical data for pembrolizumab-treated patients included smoking history (ever/never), pembrolizumab line of therapy (first, second, subsequent), weighted performance status, calculated from Eastern Cooperative Group (ECOG) scores (0–5) ascertained at post-test patient encounters between date of CIP report until date of last dose of pembrolizumab or death, response status (complete, partial, stable, or progressive disease), progression-free survival (PFS), calculated as the number of months between first dose of pembrolizumab and progression date, and overall survival (OS), calculated as the number of months between the first dose of pembrolizumab and date of death.

2.4. Statistical Analysis

PD-L1 assay concordance between IHC and RNA-seq was assessed by independent samples Kruskal-Wallis median test and proportion analysis (chi-square) using clinical cutoffs for IHC tumor proportion score (TPS) as high (≥50%), low (1–49%), or negative (<1%), or for RNA-seq percentile rank as high (≥75) or not high (<75). Receiver operator characteristic (ROC) analysis was also used in the concordance analysis cohort to explore the accuracy (sensitivity and specificity) of potential alternative cutoffs for RNA-seq using IHC as the gold reference standard.
In the pembrolizumab-treated cohort, we assessed the predictive value of individual and combined PD-L1 marker status by RNA-seq and IHC for patient treatment outcomes. Proportion analysis (chi-square) was used to assess for differences in treatment response by marker status, and Kaplan-Meier analysis assessed associations between marker status and survival. Cox proportional hazards multiple regression models were used to evaluate the predictive value of PD-L1 by marker status for PFS and OS, adjusting for potential covariates, including age, sex, performance status (weighted ECOG), histology, smoking history, and pembrolizumab line of treatment. All data were analyzed using SPSS Statistics for Windows, version 26.0 (IBM Corp., Armonk, NY, USA).

3. Results

3.1. Concordance between IHC 22C3 and RNA-Seq for PD-L1 Expression in NSCLC

Although IHC is scored on a 0–100 scale for TPS, results were significantly left-skewed toward cases with low (1–49% TPS) expression, with the largest proportion of cases (37.7%) also falling in this group (Figure 1A). Most cases were also classified as “not high” by RNA-seq (n = 2036, 64.3%), falling below the ≥75 percentile rank cutoff used in clinical reporting (Figure 1B). However, a greater number of cases were high by RNA-seq compared with IHC (35.7% vs. 30.8%, respectively).
Overall, IHC and RNA-seq demonstrated moderate concordance for PD-L1 expression. As seen in Figure 1C,D, the medians and distributions of IHC TPS results by RNA-seq percentile rank clinical cutoff group (high or not high) were significantly different from one another, as were the medians and distributions of RNA-seq percentile rank results by IHC TPS clinical cutoff group (high, low, or negative) (Kruskal-Wallis test with significance values adjusted by Bonferroni correction for multiple tests; p < 0.001).
Furthermore, as seen in Figure 1E, 2539 (80%) of PD-L1 IHC and RNA-seq results were concordant with the clinical cutoff group. The proportion of NSCLC patients in each IHC cutoff group (high, low, or negative) was significantly associated with the proportion of patients in the corresponding RNA-seq cutoff group (high or not high), as was the reverse scenario (Chi-Square < 0.001 for each scenario). Most tumors (n = 739, 76%) that were high by IHC were also high by RNA-seq. Similarly, most tumors that were low (n = 925, 78%) or negative (n = 875, 88%) by IHC were also not high by RNA-seq.
Approximately 20% (629/3168) of all NSCLC patient tumors in the cohort had discordant findings for PD-L1 expression by IHC and RNA-seq. The largest number of discordant cases had either low (268) or negative (125) tumoral expression by IHC but were determined to be high expressers by RNA-seq. The remaining discordant cases (236) were high by IHC but not high by RNA-seq (Figure 1E). Discordant PD-L1 expression between IHC and RNA-seq is also not readily discernable under the microscope, as illustrated by example photomicrographs (Figure 1F).

3.2. Alternative PD-L1 RNA-Seq Cutoffs and Assay Performance

Given the known suboptimal IHC test sensitivity for PD-L1, particularly at low levels of expression, we explored alternative PD-L1 RNA-seq cutoffs that might better discern IHC results as the standard reference test than the 75th percentile rank historically used in clinical reporting. Receiver operator characteristic (ROC) analyses were performed using five sets of NSCLC patient tumors with two different PD-L1 IHC TPS cutoff groups in each set: IHC high versus negative (Figure 2A); IHC low versus negative (Figure 2B); IHC high versus low (Figure 2C); IHC high + low versus negative (Figure 2D); and IHC high (TPS ≥ 50) versus low + negative (Figure 2E). A summary Youden index score (sensitivity + 1-specificity) was calculated for each ROC curve, with an index score ≥ 50% considered the empirical benchmark of an efficacious diagnostic test [36].
RNA-seq was best able to differentiate IHC-high versus negative tumors at a 69th percentile rank cutoff for PD-L1, demonstrating high sensitivity and specificity (Figure 2A). However, a more pressing and clinically relevant problem is that oncologists need to differentiate between IHC-high versus low tumors to support the selection of pembrolizumab monotherapy versus combination pembrolizumab plus chemotherapy. For this scenario, RNA-seq accurately classified patients in each group at a 76th percentile rank cutoff (AUC = 0.84, Youden index = 53; p < 0.0001) (Figure 2B). Similarly, RNA-seq accurately classified IHC-high versus low + negative tumors (AUC = 0.87, Youden index = 58; p < 0.0001) (Figure 2C). RNA-seq, however, could not discern IHC low versus negative tumors (Figure 2D) or IHC high + low versus negative tumors (Figure 2E), lacking sensitivity and specificity for both scenarios. Given the high sensitivity and specificity of a 76th percentile RNA-seq rank cutoff identified in ROC analysis for discerning the most clinically relevant group of PD-L1 expression cases (IHC high versus low) and its proximity to the 75th percentile rank cutoff already in use for clinically reporting a high result, subsequent analyses in the study were based on a 75th percentile rank cutoff for PD-L1 by RNA-seq.

3.3. Associations with Pembrolizumab Monotherapy Outcomes

The potential predictive value of PD-L1 by RNA-seq for pembrolizumab monotherapy response, progression-free survival (PFS), and overall survival (OS) was explored using the 75th percentile rank cutoff in a cohort of NSCLC patients who received pembrolizumab monotherapy as their first immunotherapy following testing (n = 102) (Table 1). Given recommendations from professional practice guidelines to select NSCLC patients with high PD-L1 expression for pembrolizumab, this real-world treatment cohort was biased for PD-L1 high-expressing tumors (82.3% of patients). These patients were also more likely to have a smoking history (89.2%), non-squamous histology (81.4%), a weighted ECOG performance status between 1 and 2 (61.8%), and have received their first course of pembrolizumab in the frontline setting (84.3%). Very few patients in the pembrolizumab-treated cohort received prior (non-EGFR/ALK) targeted therapy (2.9%) or immunotherapy (2%).
Of 102 pembrolizumab-treated patients, 74 (72.5%) had required scan data and measurable disease at baseline and were retrospectively evaluable for pembrolizumab response based on RECIST v1.1.
As shown in Figure 3, as an individual measure, patients with high PD-L1 expression by RNA-seq had a significantly better response rate (RR) for pembrolizumab than those whose tumors were below the 75th percentile rank cutoff for high expression (65% vs. 18%, p < 0.001). In contrast, we observed no significant differences in response rates for pembrolizumab monotherapy between IHC high (≥50% TPS) versus low (<1% TPS) expression (54% vs. 38%, p = 0.306).
Patients with combined PD-L1 “double high” status by both RNA-seq and IHC had the best response rate of any group (66%). Patients with “double low” (IHC low + RNA not high) status fared better than patients who were high by IHC but not high by RNA-seq (25% vs. 14%, respectively).
Comparing response rates for combined versus single marker status, patients that were high by RNA-seq + low by IHC had better response rates than patients who were high by IHC alone (60% vs. 54%). Patients who were low by IHC and not high by RNA-seq had worse response rates than patients who were low by IHC alone (25% vs. 38%).
Kaplan-Meier analysis was used to evaluate associations between PD-L1 marker status for IHC and RNA-seq with median progression-free survival (PFS) and overall survival (OS) for pembrolizumab-treated NSCLC patients in the validation cohort. There were no significant differences between median PFS or OS for patients with PD-L1 IHC high (≥50% TPS) versus low (1–49% TPS) tumors (Figure 4A,B). In contrast, patients with RNA-seq high (≥75 rank) versus not high (<75 rank) expression had significantly improved median PFS (19.2 versus 4.0 months, p < 0.001) (Figure 4C) and median OS (24.3 vs. 5.4 months, p = 0.035) (Figure 4D).
For combined PD-L1 marker status, patients with PD-L1 RNA-high tumors had significantly better PFS than patients with tumors that were not high by PD-L1 RNA, regardless of PD-L1 IHC status (Figure 4E). Specifically, patients with PD-L1 “double high” results had the best PFS (18 months), while patients who were high by RNA-seq and low by IHC had the next best PFS (14.1 months). However, patients with RNA-seq results that were not high had significantly worse PFS, even when IHC results were high (3.2 months) and low positive (5.7 months). Combined PD-L1 marker status did not achieve statistical significance for OS; however, results followed the same trend as PFS, with RNA-seq high vs. not high status demonstrating better survival, irrespective of PD-L1 IHC results (Figure 4F).
Cox proportional hazards regression models were used to assess the predictive value of IHC and RNA-seq PD-L1 marker status for pembrolizumab PFS and OS, adjusting for potential covariates—age, sex, histologic subtype (squamous vs. non-squamous), smoking history, ECOG performance status, and line of treatment (Figure 5). Performance status was a significant covariate (p < 0.001) in all multiple regression analyses. As an individual measure, only RNA-seq high (compared with not high) status was significantly associated with a reduction in the likelihood of death for both PFS (HR = 0.62; 95% CI = 0.42–0.82, p < 0.001) and OS (HR = 0.45; 95% CI = 0.24–0.84, p = 0.05).
For combined marker status, patients with PD-L1 “double high” results had the most significant reduction in risk of death for PFS compared with IHC low + RNA-seq not high status (HR = 0.32; CI = 0.14–0.77, p = 0.01). Combined marker status was not significant for OS; however, results were directionally aligned with PD-L1 RNA marker status; when RNA results were high, patients had better HR for OS, and when RNA was not high, patients had worse HR for OS, regardless of IHC status.

4. Discussion

Higher levels of PD-L1 expression have been correlated with increased clinical benefit from immunotherapy compared with chemotherapy, but current IHC methods for PD-L1 measurement are imprecise for predictive use in NSCLC. Intratumoral heterogeneity, error-prone, subjective scoring methods, and limited tissue samples present challenges for reliable quantification of PD-L1 by IHC in NSCLC and lead to a lack of clinical test sensitivity [8,9,37,38]. These challenges are of particular concern to treating oncologists, as most NSCLC patient tumors do not have clear negative (<1% TPS) or high (≥50% TPS) PD-L1 expression but rather harbor low expression (1–49%). PD-L1 low expression is the subgroup that is also likely to be misclassified for PD-L1 status by IHC, contributing to unclear clinical efficacy for these patients [13,39].
Given the problems with current IHC testing, our study sought to evaluate the potential clinical utility of transcriptomic PD-L1 in a cohort of NSCLC patient tumors simultaneously tested by standard IHC and by RNA-seq. RNA-seq demonstrated reasonable overall concordance with IHC for PD-L1 scores and their interpretations using standard clinical reporting cutoffs for each test (Figure 1C,D), with 80% of cases demonstrating concordant results (Figure 1E). Using ROC analysis, our data also showed PD-L1 by RNA-seq accurately discerns IHC high (≥50% TPS) from IHC low (1–49% TPS) cases at a 76th percentile rank score (Figure 2B,C), a cutoff that is reasonably equivalent to the 75th percentile rank score already used in RNA-seq clinical reporting.
Our intent was to use IHC as the gold standard to define cutoffs for high, moderate, and low PD-L1 expression by RNA-seq. Problematically, however, we observed a significantly skewed distribution of TPS scores among PD-L1 IHC-low tumors toward 1% in a possible 1–49% range (Figure 1A). This reflects IHC 22C3 companion diagnostic scoring guidelines for pathologists, which call for calculation of the percentage of viable tumor cells showing linear membranous staining of any intensity and may ultimately drive over-representation of patients with very low PD-L1 expression by IHC, as well as some of the observed discord between IHC and RNA-seq in our study. As would be expected, a high proportion (78%) of PD-L1 IHC-low cases were also classified as not high for PD-L1 by RNA-seq (Figure 1E). Thus, our attempts to identify a cutoff within the RNA-seq “not high” group that would accurately discern IHC low (1–49% TPS) from IHC negative (<1% TPS) tumors failed (Figure 2D,E). This underscores the problems arising from the compulsory use of IHC as a reference standard and ultimately limits our evaluation of the clinical utility of RNA-seq as an individual measure beyond “high” versus “not high” status.
As individual measures, PD-L1 RNA-seq high versus low status was associated with improved pembrolizumab outcomes for all measures—response, progression-free survival (PFS), and overall survival (OS). We found similar trends but no significant associations with outcomes for PD-L1 IHC high or low versus negative status. Perhaps of greatest import, the evaluation of outcome associations for combined marker status demonstrated clear potential in every scenario for RNA-seq to further stratify tumors with high or low PD-L1 expression by IHC that are either more or less likely to respond and survive. Given that it has been shown in prospective clinical trials that PD-L1 expression by IHC in either tumor or immune cells is predictive of response to checkpoint blockade in NSCLC [28,40,41], our RNA-seq findings may be evidence that it is unnecessary to discern between tumor and stroma and that better assay sensitivity by NGS is more critical.
Of note, we found that NSCLC tumors for patients that were simultaneously high for PD-L1 expression by IHC but not by RNA-seq had a worse response rate (14%) than those with low expression by IHC alone (38%) or with “double low” expression by IHC and RNA-seq (25%). Even when PD-L1 is expressed, the absence of tumor-infiltrating lymphocytes in the tumor microenvironment, required to attack cancer cells, will limit pembrolizumab efficacy in these “cold” tumors. Another potential explanation in these cases, however, may be that IHC detected high levels of stable but non-functional PD-L1 protein and that IHC does not perfectly correlate with PD-L1′s functional activity [42].
There are also several limitations to our outcome analyses. Increasing PD-L1 IHC expression trended toward improved NSCLC pembrolizumab response and survival; however, our statistically negative results deviate from positive clinical trial findings and are to be viewed with caution. First, our real-world patient data are biased toward high PD-L1 IHC expression (82.3% of patients), as these patients are the most often treated by pembrolizumab monotherapy. In turn, a limited number of patients in our analysis harbored low PD-L1 expression by IHC (n = 18; 17.7%), challenging statistical comparisons. This stands in contrast to the number and proportion of patients with RNA-seq high (n = 69; 67.7%) and not high (n = 33, 32.3%) status, which is not as disparate. Unlike prospective trials, our real-world retrospective study was not powered to evaluate outcomes by assay technology for different levels of PD-L1 expression. This consideration is particularly relevant to the strength of our combined marker status findings for outcome associations.

5. Conclusions and Future Directions

Our findings support additional study of RNA-seq as a complementary testing strategy for standard IHC for PD-L1 expression. Combined measurement of PD-L1 by IHC and RNA-seq may better discern patients who are more or less likely to benefit from checkpoint inhibition, a particular problem for the largest subgroup of NSCLC patients with low PD-L1 expression by IHC and the least certainty of response to checkpoint inhibition. While it is possible that the development of new, more sensitive IHC technology can improve upon the current state of IHC testing for PD-L1, all biomarkers are susceptible to detection method bias. Furthermore, while it’s true that standard RNA sequencing does not inherently discriminate between tumor and stromal cells, advanced computational methods can help deconvolute the data to estimate the contribution of different cell types to the overall expression. These methods use known gene expression signatures associated with cell type-specific markers to infer their presence in a sample. This can provide insights into the contribution of stromal cells to PD-L1 expression.
Perhaps more importantly, growing evidence substantiates that the use of any one single marker for immunotherapy selection is a suboptimal strategy. Determining the likelihood of immunotherapy response requires measuring many more anti-cancer immune cycle responses and resistance markers in the tumor microenvironment, including those related to neoantigens, antigen presentation, DNA repair, oncogenic pathways, and their mediators [43]. Broad, comprehensive immune profiling by RNA-seq that includes PD-L1 in addition to expression for many other immune genes can be run in parallel with comprehensive genomic profiling to measure the growing spectrum of biomarkers needed for treatment selection in NSCLC. For example, in our own work, we have integrated tumor immunogenic and cell proliferation signatures with PD-L1 IHC and TMB testing and have shown that tumors that are at least moderately immunogenic have significantly better survival for checkpoint inhibition compared with highly proliferative tumors with weak immunogenicity. These signatures also identified patients who may benefit from checkpoint inhibition but were unlikely to respond based on standard PD-L1 IHC and TMB marker status alone [44]. Other investigations have also contributed new knowledge toward evolving the future use of RNA-seq signatures and multi-marker approaches to better assess responses to immunotherapy [45,46,47,48].
Advancing multi-marker approaches toward standard clinical use requires the inclusion of comprehensive immune profiling in early-stage diagnostic and drug development strategies using NGS-based technologies to overcome the technical challenges of IHC-based methods.

Author Contributions

Conceptualization, M.K.N., S.P. and J.M.C.; methodology, M.K.N., S.P. and J.M.C.; formal analysis, M.K.N. and Y.H.L.; investigation, M.K.N., R.A.P., G.K.D., L.D., P.D., S.Z., N.M., E.S., S.R., S.P. and J.M.C.; data curation, Y.H.L. and N.M.; writing—original draft preparation, M.K.N. and J.M.C.; writing—review and editing, R.A.P., G.K.D., L.D., P.D., E.S. and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Western Institutional Review Board (IRB) Copernicus Group (WCG protocol # 1340120), approved 9/2/2022.

Informed Consent Statement

The WCG IRB determined this study met requirements for a waiver of informed consent under 45 CFR 46 116(f) (2018 Requirements) and 45 CFR 46.116(d) (Pre-2018 Requirements).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

Mary K. Nesline is/was an employee of Labcorp and declares the following relationships: stock ownership from Labcorp. Rebecca A. Previs is/was an employee of Labcorp and declares the following relationships: stock ownership from Labcorp. She has served on advisory boards for Myriad Genetics and Natera. Grace K. Dy declares that she has no conflict of interest that might be relevant to the contents of this manuscript. Lei Deng declares that he has no conflict of interest that might be relevant to the contents of this manuscript. Yong Hee Lee is/was an employee of OmniSeq. Paul DePietro is/was an employee of OmniSeq and declares the following relationships: stock ownership from Labcorp. Shengle Zhang is/was an employee of OmniSeq and declares the following relationships: stock ownership from Labcorp. Nathan Meyers is/was an employee of OmniSeq. Eric Severson is/was an employee of Labcorp and declares the following relationships: stock ownership from Labcorp. Shakti Ramkissoon is/was an employee of Labcorp and declares the following relationships: stock ownership from Labcorp. Sarabjot Pabla is/was an employee of OmniSeq and declares the following relationships: stock ownership from Labcorp. Jeffrey M. Conroy is/was an employee of OmniSeq and declares the following relationships: stock ownership from Labcorp.

References

  1. Remon, J.; Passiglia, F.; Ahn, M.J.; Barlesi, F.; Forde, P.M.; Garon, E.B.; Gettinger, S.; Goldberg, S.B.; Herbst, R.S.; Horn, L.; et al. Immune Checkpoint Inhibitors in Thoracic Malignancies: Review of the Existing Evidence by an IASLC Expert Panel and Recommendations. J. Thorac. Oncol. 2020, 15, 914–947. [Google Scholar] [CrossRef] [PubMed]
  2. Shields, M.D.; Marin-Acevedo, J.A.; Pellini, B. Immunotherapy for Advanced Non–Small Cell Lung Cancer: A Decade of Progress. Am. Soc. Clin. Oncol. Educ. Book 2021, 41, e105–e127. [Google Scholar] [CrossRef] [PubMed]
  3. Postow, M.A.; Callahan, M.K.; Wolchok, J.D. Immune Checkpoint Blockade in Cancer Therapy. J. Clin. Oncol. 2015, 33, 1974–1982. [Google Scholar] [CrossRef] [PubMed]
  4. Topalian, S.L.; Drake, C.G.; Pardoll, D.M. Immune Checkpoint Blockade: A Common Denominator Approach to Cancer Therapy. Cancer Cell 2015, 27, 450–461. [Google Scholar] [CrossRef] [PubMed]
  5. Ribas, A.; Wolchok, J.D. Cancer immunotherapy using checkpoint blockade. Science 2018, 359, 1350–1355. [Google Scholar] [CrossRef]
  6. Stenehjem, D.; Lubinga, S.; Betts, K.A.; Tang, W.; Jenkins, M.; Yuan, Y.; Hartman, J.; Rao, S.; Lam, J.; Waterhouse, D. Treatment patterns in patients with metastatic non-small-cell lung cancer in the era of immunotherapy. Futur. Oncol. 2021, 17, 2940–2949. [Google Scholar] [CrossRef]
  7. Dako. PD-L1 IHC 22C3 pharmDx: Non-Small Cell Lung Cancer [Interpretation Manual]; Dako: Santa Clara, CA, USA, 2021; Available online: https://www.agilent.com/cs/library/usermanuals/public/29158_pd-l1-ihc-22C3-pharmdx-nsclc-interpretation-manual.pdf (accessed on 10 May 2023).
  8. Ben Dori, S.; Aizic, A.; Sabo, E.; Hershkovitz, D. Spatial heterogeneity of PD-L1 expression and the risk for misclassification of PD-L1 immunohistochemistry in non-small cell lung cancer. Lung Cancer 2020, 147, 91–98. [Google Scholar] [CrossRef]
  9. Hwang, D.M.; Albaqer, T.; Santiago, R.C.; Weiss, J.; Tanguay, J.; Cabanero, M.; Leung, Y.; Pal, P.; Khan, Z.; Lau, S.C.M.; et al. Prevalence and Heterogeneity of PD-L1 Expression by 22C3 Assay in Routine Population-Based and Reflexive Clinical Testing in Lung Cancer. J. Thorac. Oncol. 2021, 16, 1490–1500. [Google Scholar] [CrossRef]
  10. Wang, Y.; Wu, J.; Deng, J.; She, Y.; Chen, C. The detection value of PD-L1 expression in biopsy specimens and surgical resection specimens in non-small cell lung cancer: A meta-analysis. J. Thorac. Dis. 2021, 13, 4301–4310. [Google Scholar] [CrossRef]
  11. Xu, H.; Chen, X.; Lin, D.; Zhang, J.; Li, C.; Zhang, D.; Zhang, X. Conformance assessment of PD-L1 expression between primary tumour and nodal metastases in non-small-cell lung cancer. Onco. Targets Ther. 2019, 12, 11541–11547. [Google Scholar] [CrossRef]
  12. Moutafi, M.K.; Tao, W.; Huang, R.; Haberberger, J.; Alexander, B.; Ramkissoon, S.; Ross, J.S.; Syrigos, K.; Wei, W.; Pusztai, L.; et al. Comparison of programmed death-ligand 1 protein expression between primary and metastatic lesions in patients with lung cancer. J. Immunother. Cancer 2021, 9, e002230. [Google Scholar] [CrossRef]
  13. Naso, J.R.; Banyi, N.; Al-Hashami, Z.; Zhu, J.; Wang, G.; Ionescu, D.N.; Ho, C. Discordance in PD-L1 scores on repeat testing of non-small cell lung carcinomas. Cancer Treat. Res. Commun. 2021, 27, 100353. [Google Scholar] [CrossRef]
  14. McLaughlin, J.; Han, G.; Schalper, K.A.; Carvajal-Hausdorf, D.; Pelekanou, V.; Rehman, J.; Velcheti, V.; Herbst, R.; LoRusso, P.; Rimm, D.L. Quantitative Assessment of the Heterogeneity of PD-L1 Expression in Non–Small-Cell Lung Cancer. JAMA Oncol. 2016, 2, 46. [Google Scholar] [CrossRef] [PubMed]
  15. Kowanetz, M.; Zou, W.; Gettinger, S.N.; Koeppen, H.; Kockx, M.; Schmid, P.; Kadel Iii, E.E.; Wistuba, I.; Chaft, J.; Rizvi, N.A.; et al. Differential regulation of PD-L1 expression by immune and tumor cells in NSCLC and the response to treatment with atezolizumab (anti-PD-L1). Proc. Natl. Acad. Sci. USA 2018, 115, E10119–E10126. [Google Scholar] [CrossRef] [PubMed]
  16. Dolled-Filhart, M.; Roach, C.; Toland, G.; Stanforth, D.; Jansson, M.; Lubiniecki, G.M.; Ponto, G.; Emancipator, K. Development of a companion diagnostic for pembrolizumab in non-small cell lung cancer using immunohistochemistry for programmed death ligand-1. Arch. Pathol. Lab. Med. 2016, 140, 1243–1249. [Google Scholar] [CrossRef]
  17. Roach, C.; Zhang, N.; Corigliano, E.; Jansson, M.; Toland, G.; Ponto, G.; Dolled-Filhart, M.; Emancipator, K.; Stanforth, D.; Kulangara, K. Development of a Companion Diagnostic PD-L1 Immunohistochemistry Assay for Pembrolizumab Therapy in Non–Small-cell Lung Cancer. Appl. Immunohistochem. Mol. Morphol. 2016, 24, 392–397. [Google Scholar] [CrossRef] [PubMed]
  18. Isaacs, J.; Clarke, J.; Ready, N. Keynote 42: Pembrolizumab, PD-L1, and where to draw the line. Ann. Transl. Med. 2020, 8, 517. [Google Scholar] [CrossRef]
  19. Mok, T.S.K.; Wu, Y.; Kudaba, I.; Kowalski, D.M.; Cho, B.C.; Turna, H.Z.; Castro, G.; Srimuninnimit, V.; Laktionov, K.K.; Bondarenko, I.; et al. Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): A randomised, open-label, controlled, phase 3 trial. Lancet 2019, 6736, 1819–1830. [Google Scholar] [CrossRef] [PubMed]
  20. Morgensztern, D. KEYNOTE-042 and the role for single agent pembrolizumab in patients with PD-L1 tumor proportion score 1-49. J. Thorac. Dis. 2019, 11 (Suppl. 15), S1963–S1965. [Google Scholar] [CrossRef]
  21. Shah, A.T.; Neal, J.W. The role of Anti-PD-1/PD-L1 monotherapy as first-line treatment of metastatic NSCLC without targetable mutations and PD-L1 TPS 1–49%. Precis. Cancer Med. 2021, 4, 9. [Google Scholar] [CrossRef]
  22. National Comprehensive Cancer Network. NCCN Clinical Practice Guidelines in Oncology (NCCN Guidelines®) Non-Small Cell Lung Cancer, Version 3.2022; National Comprehensive Cancer Network: Fort Washington, PA, USA, 2022; Available online: https://www.nccn.org/professionals/physician_gls/default.aspx (accessed on 10 May 2023).
  23. Genentech. Tecentriq® (Atezolizumab) [Package Insert]; U.S. Food and Drug Administration: South San Francisco, CA, USA, 2022. Available online: https://www.accessdata.fda.gov/scripts/cder/daf/index.cfm?event=overview.process&ApplNo=761034 (accessed on 10 May 2023).
  24. FDA Approves Atezolizumab for First-Line Treatment of Metastatic NSCLC with High PD-L1 Expression. Available online: https://www.fda.gov/drugs/resources-information-approved-drugs/fda-approves-atezolizumab-first-line-treatmentmetastatic-nsclc-high-pd-l1-expression (accessed on 10 May 2023).
  25. Tsimafeyeu, I.; Imyanitov, E.; Zavalishina, L.; Raskin, G.; Povilaitite, P.; Savelov, N.; Kharitonova, E.; Rumyantsev, A.; Pugach, I.; Andreeva, Y.; et al. Agreement between PDL1 immunohistochemistry assays and polymerase chain reaction in non-small cell lung cancer: CLOVER comparison study. Sci. Rep. 2020, 10, 3928. [Google Scholar] [CrossRef] [PubMed]
  26. Scheel, A.H.; Dietel, M.; Heukamp, L.C.; Jöhrens, K.; Kirchner, T.; Reu, S.; Rüschoff, J.; Schildhaus, H.-U.; Schirmacher, P.; Tiemann, M.; et al. Harmonized PD-L1 immunohistochemistry for pulmonary squamous-cell and adenocarcinomas. Mod. Pathol. 2016, 29, 1165–1172. [Google Scholar] [CrossRef] [PubMed]
  27. Rehman, J.A.; Han, G.; Carvajal-Hausdorf, D.E.; Wasserman, B.E.; Pelekanou, V.; Mani, N.L.; McLaughlin, J.; Schalper, K.A.; Rimm, D.L. Quantitative and pathologist-read comparison of the heterogeneity of programmed death-ligand 1 (PD-L1) expression in non-small cell lung cancer. Mod. Pathol. 2017, 30, 340–349. [Google Scholar] [CrossRef] [PubMed]
  28. Fehrenbacher, L.; Spira, A.; Ballinger, M.; Kowanetz, M.; Vansteenkiste, J.; Mazieres, J.; Park, K.; Smith, D.; Artal-Cortes, A.; Lewanski, C.; et al. Atezolizumab versus docetaxel for patients with previously treated non-small-cell lung cancer (POPLAR): A multicentre, open-label, phase 2 randomised controlled trial. Lancet 2016, 387, 1837–1846. [Google Scholar] [CrossRef]
  29. Herbst, R.S.; Giaccone, G.; de Marinis, F.; Reinmuth, N.; Vergnenegre, A.; Barrios, C.H.; Morise, M.; Felip, E.; Andric, Z.; Geater, S.; et al. Atezolizumab for First-Line Treatment of PD-L1–Selected Patients with NSCLC. N. Engl. J. Med. 2020, 383, 1328–1339. [Google Scholar] [CrossRef]
  30. Erber, R.; Stöhr, R.; Herlein, S.; Giedl, C.; Rieker, R.J.; Fuchs, F.; Ficker, J.H.; Hartmann, A.; Veltrup, E.; Wirtz, R.M.; et al. Comparison of PD-L1 mRNA Expression Measured with the CheckPoint Typer® Assay with PD-L1 Protein Expression Assessed with Immunohistochemistry in Non-small Cell Lung Cancer. Anticancer Res. 2017, 37, 6771–6778. [Google Scholar] [CrossRef]
  31. Conroy, J.M.; Pabla, S.; Nesline, M.K.; Glenn, S.T.; Papanicolau-Sengos, A.; Burgher, B.; Andreas, J.; Giamo, V.; Wang, Y.; Lenzo, F.L.; et al. Next generation sequencing of PD-L1 for predicting response to immune checkpoint inhibitors. J. Immunother. Cancer 2019, 7, 18. [Google Scholar] [CrossRef]
  32. Ma, W.; De Dios, I.; Antzoulatos, S.; Estella, J.; Albitar, M. Measuring PD-L1 mRNA expression using targeted next generation sequencing and correlation with interferon pathway activation genes in lung and colorectal cancers. J. Clin. Oncol. 2020, 38 (Suppl. 15), e21534. [Google Scholar] [CrossRef]
  33. Imyanitov, E.N.; Ivantsov, A.O.; Tsimafeyeu, I.V. Harmonization of Molecular Testing for Non-Small Cell Lung Cancer: Emphasis on PD-L1. Front. Oncol. 2020, 10, 549198. [Google Scholar] [CrossRef]
  34. Theelen, W.S.M.E.; Kuilman, T.; Schulze, K.; Zou, W.; Krijgsman, O.; Peters, D.D.G.C.; Cornelissen, S.; Monkhorst, K.; Sarma, P.; Sumiyoshi, T.; et al. Absence of PD-L1 expression on tumor cells in the context of an activated immune infiltrate may indicate impaired IFNγ signaling in non-small cell lung cancer. PLoS ONE 2019, 14, e0216864. [Google Scholar] [CrossRef]
  35. Conroy, J.M.; Pabla, S.; Glenn, S.T.; Burgher, B.; Nesline, M.; Papanicolau-Sengos, A.; Andreas, J.; Giamo, V.; Lenzo, F.L.; Hyland, F.C.L.; et al. Analytical Validation of a Next-Generation Sequencing Assay to Monitor Immune Responses in Solid Tumors. J. Mol. Diagn. 2018, 20, 95–109. [Google Scholar] [CrossRef] [PubMed]
  36. Youden, W.J. Index for rating diagnostic tests. Cancer 1950, 3, 32–35. [Google Scholar] [CrossRef] [PubMed]
  37. Lantuejoul, S.; Sound-Tsao, M.; Cooper, W.A.; Girard, N.; Hirsch, F.R.; Roden, A.C.; Lopez-Rios, F.; Jain, D.; Chou, T.Y.; Motoi, N.; et al. PD-L1 Testing for Lung Cancer in 2019: Perspective from the IASLC Pathology Committee. J. Thorac. Oncol. 2020, 15, 499–519. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, M.; Wang, S.; Trapani, J.A.; Neeson, P.J. Challenges of PD-L1 testing in non-small cell lung cancer and beyond. J. Thorac. Dis. 2020, 12, 4541–4548. [Google Scholar] [CrossRef] [PubMed]
  39. Fundytus, A.; Booth, C.M.; Tannock, I.F. How low can you go? PD-L1 expression as a biomarker in trials of cancer immunotherapy. Ann. Oncol. 2021, 32, 833–836. [Google Scholar] [CrossRef]
  40. Peters, S.; Gettinger, S.; Johnson, M.L.; Jänne, P.A.; Garassino, M.C.; Christoph, D.; Toh, C.K.; Rizvi, N.A.; Chaft, J.E.; Carcereny Costa, E.; et al. Phase II trial of atezolizumab as first-line or subsequent therapy for patients with programmed death-ligand 1-selected advanced non-small-cell lung cancer (BIRCH). J. Clin. Oncol. 2017, 35, 2781–2789. [Google Scholar] [CrossRef]
  41. Rittmeyer, A.; Barlesi, F.; Waterkamp, D.; Park, K.; Ciardiello, F.; von Pawel, J.; Gadgeel, S.M.; Hida, T.; Kowalski, D.M.; Dols, M.C.; et al. Atezolizumab versus docetaxel in patients with previously treated non-small-cell lung cancer (OAK): A phase 3, open-label, multicentre randomised controlled trial. Lancet 2017, 389, 255–265. [Google Scholar] [CrossRef]
  42. Taube, J.M.; Klein, A.; Brahmer, J.R.; Xu, H.; Pan, X.; Kim, J.H.; Chen, L.; Pardoll, D.M.; Topalian, S.L.; Anders, R.A. Association of PD-1, PD-1 ligands, and other features of the tumor immune microenvironment with response to anti-PD-1 therapy. Clin. Cancer Res. 2014, 20, 5064–5074. [Google Scholar] [CrossRef]
  43. Sholl, L.M. Biomarkers of response to checkpoint inhibitors beyond PD-L1 in lung cancer. Mod. Pathol. 2022, 35 (Suppl. 1), 66–74. [Google Scholar] [CrossRef]
  44. Pabla, S.; Seager, R.J.; Van Roey, E.; Gao, S.; Hoefer, C.; Nesline, M.K.; DePietro, P.; Burgher, B.; Andreas, J.; Giamo, V.; et al. Integration of tumor inflammation, cell proliferation, and traditional biomarkers improves prediction of immunotherapy resistance and response. Biomark. Res. 2021, 9, 56. [Google Scholar] [CrossRef]
  45. Ayers, M.; Nebozhyn, M.; Cristescu, R.; McClanahan, T.K.; Perini, R.; Rubin, E.; Cheng, J.D.; Kaufman, D.R.; Loboda, A. Molecular profiling of cohorts of tumor samples to guide clinical development of pembrolizumab as monotherapy. Clin. Cancer Res. 2019, 25, 1564–1573. [Google Scholar] [CrossRef] [PubMed]
  46. Hwang, S.; Kwon, A.Y.; Jeong, J.Y.; Kim, S.; Kang, H.; Park, J.; Kim, J.H.; Han, O.J.; Lim, S.M.; An, H.J. Immune gene signatures for predicting durable clinical benefit of anti-PD-1 immunotherapy in patients with non-small cell lung cancer. Sci. Rep. 2020, 10, 643. [Google Scholar] [CrossRef] [PubMed]
  47. Sorokin, M.; Ignatev, K.; Poddubskaya, E.; Vladimirova, U.; Gaifullin, N.; Lantsov, D.; Garazha, A.; Allina, D.; Suntsova, M.; Barbara, V.; et al. RNA sequencing in comparison to immunohistochemistry for measuring cancer biomarkers in breast cancer and lung cancer specimens. Biomedicines 2020, 8, 114. [Google Scholar] [CrossRef] [PubMed]
  48. Poma, A.M.; Bruno, R.; Pietrini, I.; Alì, G.; Pasquini, G.; Proietti, A.; Vasile, E.; Cappelli, S.; Chella, A.; Fontanini, G. Biomarkers and gene signatures to predict durable response to pembrolizumab in non-small cell lung cancer. Cancers 2021, 13, 3828. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Concordance of IHC and RNA-seq for PD-L1 expression, NSCLC (n = 3168). (A,B) Distribution IHC tumor proportion scores (TPS) and RNA-seq ranks by clinical reporting cutoffs. (C,D) Median IHC TPS scores by RNA-seq rank clinical reporting cutoff and medial RNA-seq median rank scores by IHC TPS clinical reporting cutoff, Kruskal-Wallis test with significance values adjusted by Bonferroni correction for multiple tests (p < 0.001). (E) The proportion of patients in each IHC clinical reporting cutoff group was significantly associated with the proportion of patients in the corresponding RNA clinical reporting cutoff group (Chi-Square < 0.001). (F) Photomicrographs depict representative PD-L1 IHC staining for each of the subgroups defined as high, low or negative TPS scores within each of the RNA percentile rank clinical reporting categories (high or not high). Scale bars are 100 μm.
Figure 1. Concordance of IHC and RNA-seq for PD-L1 expression, NSCLC (n = 3168). (A,B) Distribution IHC tumor proportion scores (TPS) and RNA-seq ranks by clinical reporting cutoffs. (C,D) Median IHC TPS scores by RNA-seq rank clinical reporting cutoff and medial RNA-seq median rank scores by IHC TPS clinical reporting cutoff, Kruskal-Wallis test with significance values adjusted by Bonferroni correction for multiple tests (p < 0.001). (E) The proportion of patients in each IHC clinical reporting cutoff group was significantly associated with the proportion of patients in the corresponding RNA clinical reporting cutoff group (Chi-Square < 0.001). (F) Photomicrographs depict representative PD-L1 IHC staining for each of the subgroups defined as high, low or negative TPS scores within each of the RNA percentile rank clinical reporting categories (high or not high). Scale bars are 100 μm.
Cancers 15 04789 g001
Figure 2. Receiver Operating Characteristic (ROC) analysis for alternate RNA-seq cutoff values for PD-L1 using IHC clinical cutoffs as the reference standard, NSCLC concordance cohort (n = 3168). (A) IHC high vs. IHC negative; (B) IHC High vs. IHC low; (C) IHC High vs. IHC Low + IHC Negative; (D) IHC High + IHC Low vs. IHC Negative, and; (E) IHC Low vs. IHC Negative.
Figure 2. Receiver Operating Characteristic (ROC) analysis for alternate RNA-seq cutoff values for PD-L1 using IHC clinical cutoffs as the reference standard, NSCLC concordance cohort (n = 3168). (A) IHC high vs. IHC negative; (B) IHC High vs. IHC low; (C) IHC High vs. IHC Low + IHC Negative; (D) IHC High + IHC Low vs. IHC Negative, and; (E) IHC Low vs. IHC Negative.
Cancers 15 04789 g002
Figure 3. Pembrolizumab response rates (RR) by PD-L1 marker status (n = 74).
Figure 3. Pembrolizumab response rates (RR) by PD-L1 marker status (n = 74).
Cancers 15 04789 g003
Figure 4. Kaplan-Meier survival curves for median progression free survival (PFS) and median overall survival (OS) by individual or combined PD-L1 IHC and RNA-seq marker status. (A,B) IHC high vs IHC low; (C,D) RNA high vs RNA not high, and; (E,F) IHC + RNA-seq combined status.
Figure 4. Kaplan-Meier survival curves for median progression free survival (PFS) and median overall survival (OS) by individual or combined PD-L1 IHC and RNA-seq marker status. (A,B) IHC high vs IHC low; (C,D) RNA high vs RNA not high, and; (E,F) IHC + RNA-seq combined status.
Cancers 15 04789 g004
Figure 5. Forest plot of hazard ratios (HR) for predicting progression free survival (PFS) and overall survival (OS) for pembrolizumab by individual and combined PD-L1 marker status for IHC and RNA-seq. Cox proportional regression analyses adjusted for potential covariates: age, sex, histologic subtype, smoking history, weighted performance status, and pembrolizumab line of treatment.
Figure 5. Forest plot of hazard ratios (HR) for predicting progression free survival (PFS) and overall survival (OS) for pembrolizumab by individual and combined PD-L1 marker status for IHC and RNA-seq. Cox proportional regression analyses adjusted for potential covariates: age, sex, histologic subtype, smoking history, weighted performance status, and pembrolizumab line of treatment.
Cancers 15 04789 g005
Table 1. Pembrolizumab outcomes cohort by PD-L1 RNA-seq expression status (n = 102).
Table 1. Pembrolizumab outcomes cohort by PD-L1 RNA-seq expression status (n = 102).
RNA High
(≥75 Rank) n = 69 (67.7)
RNA Not High
(<75 Rank) n = 33 (32.3)
Total (n = 102)
PD-L1 IHC statusHigh (≥50% TPS)64 (92.8)20 (60.6)84 (82.3)
Low (1–49% TPS)5 (7.2)13 (29.4)18 (17.7)
Age (mean, years)707070
SexFemale39 (56.5)20 (60.6)59 (57.8)
Male30 (43.5)13 (39.4)43 (42.2)
Ever smoker (yes)61 (88.4)30 (90.9)91 (89.2)
HistologyNon-Squamous61 (88.4)22 (66.7)83 (81.4)
Squamous8 (11.6)11 (33.3)19 (18.6)
Tissue SitePrimary34 (49.3)21 (63.8)55 (53.9)
Metastatic35 (50.7)12 (36.4)47 (46.1)
Performance
status (weighted)
<122 (31.9)3 (9.1)25 (24.5)
1 to <239 (56.5)24 (72.7)63 (61.8)
2 to <46 (8.7)5 (15.2)11 (10.8)
Brain metastases (yes)11 (15.9)4 (12.1)15 (14.7)
Pembrolizumab line of treatment159 (85.5)27 (81.8)86 (84.3)
28 (11.6)4 (12.1)12 (11.8)
≥32 (2.9)2 (6.1)4 (3.9)
Prior treatment history (yes)Chemotherapy8 (11.6)5 (15.2)13 (12.7)
Targeted therapy1 (1.4)2 (6.0)3 (2.9)
Immunotherapy2 (2.8)0 (0.0)2(2.0)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nesline, M.K.; Previs, R.A.; Dy, G.K.; Deng, L.; Lee, Y.H.; DePietro, P.; Zhang, S.; Meyers, N.; Severson, E.; Ramkissoon, S.; et al. PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes. Cancers 2023, 15, 4789. https://doi.org/10.3390/cancers15194789

AMA Style

Nesline MK, Previs RA, Dy GK, Deng L, Lee YH, DePietro P, Zhang S, Meyers N, Severson E, Ramkissoon S, et al. PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes. Cancers. 2023; 15(19):4789. https://doi.org/10.3390/cancers15194789

Chicago/Turabian Style

Nesline, Mary K., Rebecca A. Previs, Grace K. Dy, Lei Deng, Yong Hee Lee, Paul DePietro, Shengle Zhang, Nathan Meyers, Eric Severson, Shakti Ramkissoon, and et al. 2023. "PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes" Cancers 15, no. 19: 4789. https://doi.org/10.3390/cancers15194789

APA Style

Nesline, M. K., Previs, R. A., Dy, G. K., Deng, L., Lee, Y. H., DePietro, P., Zhang, S., Meyers, N., Severson, E., Ramkissoon, S., Pabla, S., & Conroy, J. M. (2023). PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes. Cancers, 15(19), 4789. https://doi.org/10.3390/cancers15194789

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop