Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Prognostic value and peripheral immunologic correlates of early FDG PET response imaging in a phase II trial of risk-adaptive chemoradiation for unresectable non-small cell lung cancer

View ORCID ProfileStephen R. Bowen, Daniel S. Hippe, Hannah M. Thomas, Balukrishna Sasidharan, Paul D. Lampe, Christina S. Baik, Keith D. Eaton, Sylvia Lee, Renato G. Martins, Rafael Santana-Davila, Delphine Chen, Paul E. Kinahan, Robert S. Miyaoka, Hubert J. Vesselle, A. McGarry Houghton, Ramesh Rengan, Jing Zeng
doi: https://doi.org/10.1101/2021.02.22.21251919
Stephen R. Bowen
1Department of Radiation Oncology, University of Washington School of Medicine, Seattle WA
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stephen R. Bowen
Daniel S. Hippe
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hannah M. Thomas
3Department of Radiation Oncology, Christian Medical College Vellore, India
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Balukrishna Sasidharan
3Department of Radiation Oncology, Christian Medical College Vellore, India
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul D. Lampe
4Human Biology Division, Fred Hutchinson Cancer Research Center, Seattle, WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christina S. Baik
5Division of Medical Oncology, Department of Medicine, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Keith D. Eaton
5Division of Medical Oncology, Department of Medicine, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sylvia Lee
5Division of Medical Oncology, Department of Medicine, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Renato G. Martins
5Division of Medical Oncology, Department of Medicine, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Rafael Santana-Davila
5Division of Medical Oncology, Department of Medicine, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Delphine Chen
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul E. Kinahan
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert S. Miyaoka
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hubert J. Vesselle
2Department of Radiology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. McGarry Houghton
4Human Biology Division, Fred Hutchinson Cancer Research Center, Seattle, WA
6Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ramesh Rengan
1Department of Radiation Oncology, University of Washington School of Medicine, Seattle WA
6Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jing Zeng
1Department of Radiation Oncology, University of Washington School of Medicine, Seattle WA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: jzeng13@uw.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Introduction We hypothesized that FDG PET imaging during chemoradiation for unresectable non-small cell lung cancer (NSCLC) is prognostic for survival, and that tumor PET response is correlated with peripheral T-cell function.

Methods 45 patients with AJCCv7 stage IIB-IIIB NSCLC enrolled on the phase II FLARE-RT trial and received platinum-doublet chemotherapy concurrent with 6 weeks of radiation (NCT02773238). FDG PET imaging was performed prior to treatment start and after 24 Gy of radiation (week 3). PET response status was prospectively defined by multifactorial radiologic interpretation. PET responders received 60 Gy in 30 fractions, while non-responders received concomitant boosts to 74 Gy in 30 fractions. Peripheral blood was drawn synchronously with PET imaging, from which germline DNA sequencing, T-cell receptor (TCR) sequencing, and plasma cytokines analysis were performed.

Results Median follow-up was 18.8 months, 1-year overall survival (OS) 82%, 1-year progression-free survival (PFS) 53%, and 1-year locoregional control (LRC) 88%. Higher mid-treatment PET total lesion glycolysis was detrimental to OS (1-yr 87% vs. 63%, p<0.001), PFS (1-yr 60% vs. 26%, p=0.044) and LRC (1-yr 94% vs. 65%, p=0.012), even after adjustment for clinical/treatment factors. Higher PD-L1 tumor proportion score (TPS) was correlated with PET response (p=0.017): 6/6 patients with high PD-L1 (TPS>50%) were classified as PET responders, while 4/5 patients classified as PET non-responders had negative PD-L1 (TPS<1%). Higher TCR richness and clone distribution slope was associated with improved OS (p=0.018-0.035); clone distribution slope was correlated with PET response (p=0.031). Germline DNA alterations in immunologic pathways had an outsized effect on OS and PET response; of the top 30 SNPs ranked by association with PET response status (p<0.016), a plurality (13/30) came from immunologic pathways.

Conclusions Mid-chemoradiation PET imaging is prognostic for survival; PET response may be linked to tumor and peripheral T-cell biomarkers.

Introduction

The current standard of care for unresectable locally advanced non-small cell lung cancer (LA-NSCLC) is concurrent chemoradiation plus a year of adjuvant durvalumab, a programmed death-ligand 1 (PD-L1) immune checkpoint inhibitor, which improved survival over chemoradiation alone1,2. However, patient outcomes remain suboptimal even for well-selected clinical trial eligible patients, with an overall survival at 2 years of 66.3%1. Clinical trials seek to improve outcomes through intensification of therapy, both in terms of additional immune modulating therapy as well as localized radiation dose escalation. Beyond durvalumab, other PD-1/PD-L1 pathway checkpoint inhibitors are also being tested as consolidation therapy after chemoradiation, concurrent with chemoradiation, and neoadjuvant to chemoradiation3. As treatment intensifies, patient tolerance is increasingly challenged, with 29.1% of patients experiencing serious adverse events in the PACIFIC trial and 15.4% of patients discontinuing treatment due to adverse events. Although durvalumab improved 2-year overall survival (OS) from 55.6% to 66.3% in the PACIFIC trial, only a minority of the patients that received durvalumab derived a survival benefit from the additional treatment. A prognostic test that could predict patient outcome would help better select patients who need additional therapies, and spare toxicity in patients who do not.

FDG PET imaging during chemoradiation for NSCLC has been correlated with treatment response and survival4-8. Van Elmpt et al. found that for patients whose mean tumor standardized uptake values (SUV) decreased by >15% on mid-treatment PET, 2-year OS was 92% compared to 33% for patients with a decrease in mean SUV <15%9. Tumor FDG PET avidity has also been correlated with PD-L1 expression, with multiple clinical series suggesting SUVmax may be associated with PD-L1 positivity, and that SUVmax could be a potential predictive marker of response to anti-PD-1 therapy in NSCLC patients10-12.

In addition to predicting for survival, FDG PET avidity may also be predictive for local recurrences13,14, both on the pre-chemoradiation scan and the mid-chemoradiation scan. Decades of trials to overcome local recurrence risk by radiation treatment intensification in unresectable NSCLC have yielded limited success. Results of RTOG 0617 showing inferior outcomes with 74 Gy versus 60 Gy during concurrent chemoradiation was against expectations and reset the standard of care radiation dose to around 60 Gy15. However, local control remains a problem in unresectable NSCLC with locoregional failures of around 20-30% by 2-years and rising to 50% by 3-5 years1,16,17. There are multiple theories regarding the lack of benefit with 74 Gy versus 60 Gy in RTOG 0617, including higher dose to normal tissues being detrimental (esophagitis grade and heart dose were both highly correlated with survival)15. Since uniform dose escalation in an unselected group of patients was detrimental to survival, multiple clinical trials are underway to test other strategies for dose escalation, to improve cancer control and survival. Some trials have tested giving all patients a radiation boost based on pre-treatment PET scans, targeting the most FDG-avid regions with a higher dose of radiation18,19, whereas other trials dose-escalate all patients based on mid-treatment PET20. Since some patients have durable cancer control after 60 Gy of radiation, not all patients may benefit from dose escalation. We designed a phase II trial FLARE-RT (NCT02773238) where select patients undergo dose escalation based on mid-treatment response as assessed by PET, and only non-responders undergo radiation dose escalation (74 Gy) for the second half of chemoradiation21.

We hypothesized that early FDG PET response imaging during chemoradiation has prognostic value in patients treated on the FLARE-RT trial, and that a robust tumor PET response is linked to peripheral T-cell function, given the clinical and preclinical data that a robust tumor response to radiation requires an intact immune system, including the presence of CD8+ T-cells22,23. The purpose of this investigation was to (i) validate the prognostic value of early FDG PET response imaging, (ii) evaluate peripheral immunologic correlates of PET imaging, and (iii) identify complementary markers of tumor biology (PET response) and disease burden (PET total lesion glycolysis). Identification of biomarkers of early treatment response/resistance can support clinical management decisions, and inform the design of next-generation biomarker-guided and risk-adaptive clinical trials to improve outcomes in patients with LA-NSCLC.

Materials and Methods

Clinical Trial Protocol

The FLARE-RT phase II trial (NCT02773238) schema is shown in Figure 1. This research was IRB approved and conducted in accord with the ethical standards of our institution. Patients with histopathologically confirmed unresectable AJCCv7 stage IIB-IIIB non-small lung cancer and ECOG performance status 0-1 were enrolled from 2016-2020. Baseline FDG PET imaging was followed by initiation of definitive chemoradiation planned for 60 Gy in 30 fractions. Chemotherapy consisted of a platinum doublet per physician choice (Table 1) and started at the same time as radiation. Consolidation immunotherapy was allowed post-chemoradiation once it became standard of care. During the third week of chemoradiation (after 24 Gy delivered), FDG PET imaging was performed to assess early treatment response. Patients prospectively classified as responders continued to receive a total dose of 60 Gy in 30 fractions, while those classified as non-responders received a concomitant dose boost over the final 15 fractions to a total dose of 74 Gy in 30 fractions, as previously described21. The precision dose boost was spatially localized and redistributed to conform to residual intratumoral FDG avidity. Peripheral blood draws were completed synchronously with PET imaging at baseline, week 3 mid-treatment, and 3-month post-treatment time points. The primary trial endpoint was 2-year OS, while secondary endpoints included 1-year PFS, 1-year LRC, and pulmonary toxicity (CTCAE v4 grade 2+ pneumonitis).

Figure 1.
  • Download figure
  • Open in new tab
Figure 1.

FLARE-RT phase II trial schema of risk-adaptive chemoradiation for patients with unresectable locally advanced non-small cell lung cancer.

View this table:
  • View inline
  • View popup
Table 1.

Patient characteristics (n=45).

PET/CT Response Assessment & Residual Total Lesion Glycolysis

Week 3 PET response status was prospectively defined by multifactorial radiologic interpretation. Quantitative assessments were based on interval changes in standardized uptake value (SUV) metrics (SUVmax, SUVmean, SUVpeak), metabolic tumor volume (MTV), and total lesion glycolysis (TLG) using semi-automatic gradient-based segmentation of the primary tumor and adjacent involved lymph nodes. PET responders presented with greater than 20% decrease in at least one metric of FDG avidity (SUVmax, SUVmean, SUVpeak) and at least one metric of FDG volumetric extent (MTV, TLG), along with clinical radiographic evidence of interval response to treatment. Metabolic tumor volume contours were defined by a commercially validated segmentation algorithm (PET Edge, MIM Software, Cleveland, OH) that achieved improved inter-observer agreement compared to manual contouring and reduced sensitivity to image reconstruction compared to fixed threshold contouring24,25. Prospective FLARE protocol-defined mid-chemoradiation PET response status showed substantial agreement to retrospective mid-chemoradiation PET response assessment by PERCIST 1.026,27 (Kappa 0.72 [0.51-0.93]). Residual total lesion glycolysis on week 3 PET (TLGmidtx [g]) was calculated as the product of MTV [mL] and SUVmean [g/mL]. Given the skewed distribution of TLGmidtx across patients, we dichotomized the upper tail of the distribution from the bulk, which was achieved at the 80th percentile corresponding to a threshold of TLGmidtx = 250 g. Dichotomized TLGmidtx defined patient risk strata in survival analyses.

Peripheral Blood Genomic & Immunologic Assays

Sequencing of peripheral germline DNA was conducted using the Infinium Global Screening Array (Illumina Inc. San Diego, CA), which included sample QC, library prep+QC, cluster optimization, and sequencing steps. Sequencing libraries were prepared with a Covaris LE220 system and PerkinElmer Sciclone NGSx workstation. Library QC was ensured with two Agilent Bioanalyzer 2100, an Agilent 2200 TapeStation, a Life Technologies Qubit 2.0 fluorometer and an ABI StepOne Real-Time PCR system. Genotypes with GenCall scores > 0.15 were considered of sufficient quality for inclusion. Genome sequencing reports were filtered to include single nucleotide polymorphism (SNP) genotypes from a list of predefined relevant candidate genes (see Supplemental Table 1) belonging to different pathway families: 96 DNA repair genes28,29, 53 immunology genes30-32, 38 oncology genes33,34, and 27 lung biology genes35,36.

View this table:
  • View inline
  • View popup
Supplemental Table 1.

Predefined Candidate Genes for 4 Pathways: DNA Repair, Immunology, Oncology, and Lung Biology

Molecular T-cell receptor (TCR) beta chain CDR3 sequencing was carried out on the ImmunoSEQ platform (Adaptive Biotechnologies, Snohomish, WA) with a survey sampling depth of 120,000 T-cells. We characterized the immune repertoire through diversity metrics, which capture the number of unique T-cell receptors (richness) and the relative abundance of each receptor (clonality / evenness). Richness was averaged over the following measures due to high collinearity (Spearman r > 0.93): iChao137, Efron Thisted Estimator38, and Daley Smith Estimator39. Clonality / evenness measures included the following: Pielou Evenness40, Simpson Clonality and Evenness41, and Clone Distribution Slope42. In general, higher richness and higher clone distribution slope indicate greater TCR diversity with the presence of more rare clones.

Exploratory single immune cell functional assays were performed for 4 patients on the IsoLight system using IsoCode Human Adaptive Immune chips (Isoplexis, Branford, CT), which can measure 25+ cytokine secretions on a single cell base of approximately 1000 individual cells per chip. We restricted our analysis to abundance of polyfunctional CD8+ T-cells, following stimulation and secretion of at least 2 effector cytokines (e.g. MIP1β, IFNγ). Lastly, detectable plasma concentrations from an immunologic panel of 43 cytokines were measured via a combination of the Luminex 200 microbead system (Luminex, Austin TX) and standard ELISA.

Statistical Analysis

Overall survival (OS), progression-free survival (PFS), and distant metastatic-free (DM) survival rates were estimated by Kaplan-Meier. Locoregional control (LRC) and pneumonitis (PNM) rates were adjusted for distant progression and death as competing risks. Differences between risk strata were evaluated with log rank or Gray’s testing. Univariate and bivariate hazard ratios (HR) were estimated by Cox regression or Fine & Gray competing risk regression.

Unsupervised learning of peripheral germline DNA sequencing SNP genotypes was performed with hierarchical clustering and summarized with cumulative frequency distributions of rank ordered SNPs for each genetic pathway family (DNA repair, immunology, oncology, lung biology). Genotype associations to PET response status were characterized by logistic regression with Firth penalization for estimating odds ratios (OR). Genotype associations to OS and PFS were characterized by Cox regression with Firth penalization for estimating hazard ratios (HR). Genotype associations to TLGmidtx were characterized by Spearman rank correlation.

Associations between PD-L1 tumor proportion score (TPS) and PET response status were summarized by logistic regression with Firth penalization, Spearman rank correlation, and Fisher exact testing. PD-L1 TPS associations with TLGmidtx were estimated by Spearman rank correlation, while associations to histology, stage, and presence of driver mutation were estimated by Fisher exact testing.

TCR diversity metric trends over pre-treatment, mid-treatment, and post-treatment timepoints were estimated by Spearman rank correlation and pairwise changes by Wilcoxon signed rank testing. TCR diversity metric associations to OS and PFS were assessed by Permutation testing of Harrell’s c-index. TCR diversity metric associations to PET response class were assessed by groupwise Wilcoxon rank sum testing. TCR diversity metric associations to TLGmidtx were estimated by Spearman rank correlation.

Unsupervised learning of plasma cytokines was performed with hierarchical clustering and summarized with a 2-dimensional dendrogram heatmap of cytokine concentration z-scores with rugs for OS, PFS, PET response, and TLGmidtx classes. Associations of cytokines with PET response status were assessed by Fisher exact testing, associations with OS and PFS were estimated by Cox regression with Firth penalization, and associations with TLGmidtx were estimated by Spearman rank correlation. All statistical analyses and data visualization were carried out in OriginPro 2020b (OriginLab, Northampton, MA) and R version 4.0.3 (R Statistical Language, Vienna, Austria).

Throughout, two-sided tests were used with p < 0.05 considered statistically significant, without adjustment for the number of comparisons.

Results

Patient Characteristics

Patient characteristics are summarized in Table 1. Forty-nine patients were enrolled but 4 patients did not initiate any clinical trial therapy (3 were found to have metastatic disease on repeat baseline imaging, and 1 patient chose to have surgery instead). In the intention-to-treat population (n=45) with median age of 63 years (range 34-78 years), most patients enrolled on the FLARE-RT trial had unresectable, stage III, multi-station N2 positive NSCLC. Approximately half of patients received consolidation immune checkpoint inhibitor therapy when the standard of care changed mid-point during this trial (Table 1), and about half the patients were treated with proton beam radiation (Table 1) and half with photon intensity-modulated radiation therapy (IMRT).

Clinical Outcomes

Figure 2 summarizes clinical outcomes of the FLARE-RT phase II trial from the time of consent. The median time between consent and start of chemoradiation was 11 days (range: 0 - 24 days). After a median follow-up of 18.8 months (range 3.0-49.9 months) in the intention-to-treat (ITT) population (n=45), 1-year OS was 82%, 1-year PFS was 53%, 1-year LRC was 88%, 1-year DM was 61%, and 6-month pneumonitis (PNM) cumulative incidence was 25%. Two of 45 patients terminated chemoradiation early due to progressive disease during chemoradiation, and switched to systemic therapy alone, but were included in the ITT analysis. Patients who received durvalumab immune checkpoint inhibitor therapy had numerically higher cumulative incidence of grade 2+ pneumonitis versus those who did not (1-yr PNM 52% vs. 26%, Gray p=0.20). Receipt of consolidation durvalumab according to the PACIFIC regimen resulted in numerically higher PFS but did not reach statistical significance in our trial cohort (median PFS 19.4 months vs. 10.9 months, log rank p=0.45).

Figure 2.
  • Download figure
  • Open in new tab
Figure 2.

FLARE-RT clinical trial outcomes: overall survival (OS, A), progression-free survival (PFS, B), PFS stratified by receipt of consolidation immune checkpoint inhibitor therapy (ICI, C), distant metastatic-free survival (DM, D), locoregional control (LRC, E), CTCAE v4 grade 2 or higher pneumonitis (PNM, F). KM: Kaplan-Meier, CR: competing risk (distant progression / death).

Week 3 Mid-treatment PET Risk Stratification

Figure 3A-C shows that residual total lesion glycolysis (TLGmidtx) on week 3 PET, dichotomized at the 80th percentile (250 g), was significantly associated with a detriment in OS (log rank p < 0.001), PFS (log rank p = 0.044), and LRC (Gray p = 0.012). In the setting of FLARE-RT trial protocol PET response-adaptive dose escalation, PET response status was not associated with outcomes (p > 0.25). However, Figure 3D captures a highly significant interaction between PET response status and TLGmidtx that modulated PFS (log rank p < 0.00001), in which the highest risk of rapid disease progression was in PET non-responders with TLGmidtx > 250 g. Of the 9 patients with TLGmidtx > 250 g, 6/9 were classified as PET responders and received standard chemoradiation without PET-guided dose escalation. Forest plots in Figure 3E-F summarize the associations of TLGmidtx with risk of death or risk of disease progression, which were maintained after bivariate adjustment with clinical and treatment factors, including radiation target volume (PTVpretx) and PET response status.

Figure 3.
  • Download figure
  • Open in new tab
Figure 3.

Kaplan-Meier overall survival (OS), competing risk-adjusted locoregional control (LRC), and Kaplan-Meier progression-free survival (PFS) stratified by week 3 mid-treatment FDG PET total lesion glycolysis (TLGmidtx) (A-C), along with interaction of total lesion glycolysis with PET response status for PFS (D). Forest plots of total lesion glycolysis OS and PFS univariate and bivariate hazard ratios, adjusted for individual clinical and treatment factors (E-F). PET-NR: PET non-responder; PET-R: PET responder.

Peripheral Germline DNA Sequencing

Figure 4 depicts the cumulative frequency of SNP genotypes of different genetic pathway families (DNA repair, immunology, lung biology, oncology) in 19 patients rank ordered by their effect size on unsupervised hierarchical cluster membership (A), risk of death (B), PET response class membership (C), and residual total lesion glycolysis correlation (D). A larger frequency of SNP gene alterations in the immunology pathways belong to the main hierarchical cluster and have significant association to overall survival (Figure 4A-B), relative to other pathway families which reside closer to the diagonal reference line of equal frequency contribution. Differences in gene alteration frequency between FLARE-RT PET non-responder and responder subgroups were observed across several pathway families, particularly immunology and oncology (Figure 4C). By contrast, there was no preferential correlation of specific pathway families with total lesion glycolysis, which all have cumulative frequency distributions that closely follow the diagonal reference line (Figure 4D). Of the top 30 SNPs ranked by association with PET response status (p < 0.016), a plurality (13/30) came from immunologic pathways, while none of these same SNPs were associated with total lesion glycolysis (p > 0.11). Immunologic pathway genetic differences in JAK1 associated with PET response status and PFS (p=0.033-0.040), while differences in STAT1 (p=0.001-0.030) and IFNγ (p=0.002-0.060) associated with PET response status and hierarchical cluster membership. Other associations with PET response status included germline genetic alterations in VEGFC (p=0.002), ALK (p=0.018) and EGFR (p=0.047). Associations with OS included immuno-oncogene alterations in ALK (p=0.012), PIK3CA (p=0.013) and TNF (p=0.041). None of these gene alterations were associated with total lesion glycolysis (p>0.38). Without adjustments for multiple comparisons, our results are hypothesis generating rather than conclusive.

Figure 4.
  • Download figure
  • Open in new tab
Figure 4.

Peripheral DNA microarray cumulative frequency of single nucleotide polymorphism (SNP) gene alterations by pathways between unsupervised hierarchical clusters (A), differential risk of death (B), PET response groups (C), and PET total lesion glycolysis (TLGmidtx) (D). The diagonal reference line represents equal frequency contributions from all SNPs. The frequency of gene alterations in immunologic pathways has an outsized effect on risk of death and hierarchical cluster membership relative to other pathways. PET response status shows highly significant association to gene alterations across several pathways, while PET TLGmidtx correlations to gene alterations are more randomly distributed without linkages to specific pathways. OR: odds ratio, HR: hazard ratio, r: Spearman rank correlation.

PET Response vs. Tumor PD-L1 Expression

PET response status was correlated with PD-L1 tumor proportion score (TPS) in 20 patients: 6/6 patients with high PD-L1 TPS (≥50%) were PET responders and 6/7 with moderate PD-L1 TPS (1-49%) were PET responders. By contrast, 4/5 patients classified as PET non-responders had negative PD-L1 TPS (<1%). The association between PET response and PD-L1 TPS was statistically significant as a trend across 3 expression levels (<1%, 1-49%, ≥50%, Firth p=0.017), rank correlation (Spearman r = 0.54, p=0.014) and categorical frequency (<1% vs. ≥1%, Fisher p=0.031). The association to PET response status was not confounded by availability of PD-L1 testing results: 15/20 were PET responders when PD-L1 was available, compared to 14/25 PET responders when unavailable (Fisher p=0.22). Tumor PD-L1 expression was neither correlated to mid-treatment total lesion glycolysis (Spearman r = 0.04, p=0.88) nor associated to histology, stage, or presence of a driver mutation (p>0.28).

TCR Diversity Metrics & CD8+ T-cell Polyfunctionality

TCR sequencing was performed at pre-treatment (n=15), mid-treatment (n=15), and post-treatment (n=10) time points. No patients received immunotherapy pre-chemoradiation treatment or mid-treatment, and 10/10 patients were receiving immunotherapy (durvalumab) at the post-treatment time point. TCR average richness declined significantly across time points (Figure 5A, p=0.001). Pairwise changes were most significant between pre-treatment and post-treatment time points (Wilcoxon signed rank p=0.037). Higher pre-treatment TCR richness (Figure 5B, Permutation p=0.018), higher pre-treatment clone distribution slope (Figure 5D, Permutation p=0.035), and smaller decline in clone distribution slope (Permutation p=0.050) were associated with improved OS. Pre-treatment TCR clone distribution slope, pre- and mid-treatment clonality, and mid-treatment evenness were correlated with PET response status (Wilcoxon rank-sum p=0.031-0.048). None of the TCR diversity metrics were strongly correlated with mid-treatment residual total lesion glycolysis (median Spearman |r| = 0.17 [0.04-0.52]). Patients classified as mid-treatment PET responders had higher percentage of peripheral polyfunctional CD8+ T-cells compared to mid-treatment PET non-responders (Figure 5C, 0.7-1.3% vs. 0-0.1% expressing MIP1β and IFNγ). There was no association between total lesion glycolysis and peripheral CD8+ T-cell functionality.

Figure 5.
  • Download figure
  • Open in new tab
Figure 5.

T-cell receptor (TCRβ CDR3) richness boxplots over pre/mid/post-treatment timepoints (A), TCR richness and clone distribution slope boxplots grouped by overall survival events with Kaplan-Meier insets (B,D), and CD8+ T-cell polyfunctionality for pairs of PET responders and PET non-responders (C).

Plasma Cytokine Levels

Figure 6 displays the mid-treatment peripheral blood plasma cytokine hierarchical clustering dendrogram with heatmap scaled by z-scores of cytokine concentration levels and rugs for OS, PFS, and PET response status. While the 23 patients were not distinctly clustered by cytokine levels, hierarchical clusters of cytokines included the following: (i) MIP1β, IFNγ, TNFR1, TNFR2, TNFα; (ii) VEGF, TGFβ1. Among these select cytokines, lower mid-treatment TNFR1 relative to baseline was associated with worse PFS (HR 0.43, p=0.005) and higher mid-treatment TNFα relative to baseline was associated with worse OS (HR 1.59, p=0.05). Higher mid-treatment TGFβ1 showed a trend towards correlation with PET response (Fisher p=0.065).

Figure 6.
  • Download figure
  • Open in new tab
Figure 6.

Mid-treatment peripheral blood plasma cytokines hierarchical cluster dendrogram heatmap with rugs for overall survival (OS), progression-free survival (PFS), and FLARE PET response status.

Discussion

We conducted a phase II trial testing selective radiation dose escalation in patients classified as non-responders on mid-chemoradiation PET/CT (FLARE-RT trial), with the objective that this would improve overall survival due to better local control in patients at highest risk of local recurrence. All of the patients received concurrent platinum doublet chemotherapy, and approximately half the patients received durvalumab post-chemoradiation in accordance to the PACIFIC trial results that changed the standard of care for this population2. Although our trial population is not identical to the PACIFIC trial population (which only included patients who recovered from toxicity of chemoradiation) nor the patients in RTOG 0617 (none received durvalumab post-chemoradiation)15, our results compare favorably against both trials: our 1-year OS was 82% (vs 83.1% in PACIFIC and 80.0% in RTOG 0617), our 1-year PFS was 53% (vs 55.9% in PACIFIC and 49.2% in RTOG 0617), and our 1-year LRC was 88% (vs 83.7% in RTOG 0617, and 12.6% lung recurrence plus 6.5% nodal recurrence in PACIFIC). We saw an increase in both PFS and pneumonitis with durvalumab, although neither were statistically significant in our relatively small sample size. In our trial, disease recurrence was driven primarily by distant metastatic progression, with distant metastasis being a site of first recurrence in 19/23 (83%) patients. This suggests improvement in systemic disease control is needed in a subset of patients to further improve patient outcomes.

We found that mid-chemoradiation PET imaging independently predicts for OS and PFS. Multiple prior series have found changes in SUVmax during chemoradiation to be prognostic for survival 9,43. In our trial, being a PET responder versus non-responder no longer conferred a difference in OS or PFS, perhaps due to the fact we dose-escalated patients based on PET changes during chemoradiation, mitigating the prognostic power of the PET response. However, similar to other published reports4,9,10, we did find that mid-treatment TLG independently predicted for OS and PFS, potentially serving as a marker of disease burden, even after adjustments for clinical and treatment factors. Patients classified as PET non-responders on the mid-treatment scan who also presented with a high residual TLG had the worst outcomes. This strategy could potentially help risk-stratify patients for further treatment intensification.

Our findings suggest a correlation between mid-treatment PET response to chemoradiation and peripheral immunologic status, although our results are only hypothesis generating and not conclusive. Mid-treatment PET response status was associated to PD-L1 expression, with PD-L1 positive tumors likely to be PET responders and PD-L1 negative tumors likely to be PET non-responders. Tumor response to radiation treatment requires an intact immune system with functional T-cells22, and radiation has been shown to induce a pro-inflammatory tumor microenvironment via the induction of an immunogenic cell death44,45. Since PD-L1 expression plays a major role in suppressing adaptive immunity, it is possible that radiation treatment could produce a pro-inflammatory environment that helps overcome PD-L1 induced immune suppression.

Our analysis of germline DNA SNPs found differences in gene alteration frequency between PET responders and non-responder, particularly in immunology pathways, which also correlated with survival. The JAK-STAT signaling pathway plays critical roles in cytokine receptors and can modulate the polarization of T helper cells. A recent report by Shahamatdar et al.46 analyzing germline variants in the TCGA cohort demonstrated that host genetics are associated with phenotypes that describe the immune component of the tumor microenvironment; they found one SNP associated with the amount of infiltrating follicular helper T cells; and 23 candidate genes, some of which are involved in cytokine-mediated signaling. Our patient cohort supports prior findings that patients with greater TCR richness in the peripheral blood had improved survival compared with patients with less clonal richness47-49. TCR richness declined with chemoradiation although a smaller mid-treatment decline in TCR diversity was associated with improved survival. Pre-treatment TCR clone distribution slope, among other diversity metrics, was correlated with PET response status but none of the TCR diversity metrics were correlated with mid-treatment residual total lesion glycolysis. This suggests that PET-response status may be a marker of cancer/patient biology, whereas TLG may be an independent measure of cancer disease burden. Single cell functional assays were performed on the peripheral blood for a small subset of patients, suggesting mid-treatment PET responders had higher percentages of peripheral polyfunctional CD8+ T-cells compared to non-responders. We also saw that select plasma cytokines are associated with OS and PFS and correlated with PET response, including TNFα, TNFR1, TGFβ, and MIP1β.

In conclusion, in a prospective phase 2 trial of response-adaptive radiation dose escalation for patients with unresectable NSCLC, we confirm that mid-chemoradiation PET imaging is prognostic for survival and found that PET response status may be linked with peripheral T-cell function. The combination of PET response imaging and peripheral blood biomarkers could be used to guide further clinical trials of treatment intensification by identifying patients at highest risk of treatment failure.

Data Availability

Data will be made publicly available and upon request at the conclusion of the clinical trial.

Acknowledgments

This work was funded by NIH / NCI R01CA204301. We thank Priya Vissamraju and Christina Lo for coordination of the FLARE-RT clinical trial. We acknowledge the efforts of nuclear medicine and radiation oncology technologists during PET/CT acquisitions, radiation therapy planning, and image-guided radiation therapy delivery.

Footnotes

  • Financial Support: NIH/NCI R01CA204301

References

  1. 1.↵
    Antonia SJ, Villegas A, Daniel D, et al. Overall Survival with Durvalumab after Chemoradiotherapy in Stage III NSCLC. The New England journal of medicine 2018; 379(24): 2342–50.
    OpenUrlCrossRefPubMed
  2. 2.↵
    Antonia SJ, Villegas A, Daniel D, et al. Durvalumab after Chemoradiotherapy in Stage III Non-Small-Cell Lung Cancer. The New England journal of medicine 2017; 377(20): 1919–29.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Zeng J, Bowen SR. Treatment Intensification in Locally Advanced/Unresectable NSCLC Through Combined Modality Treatment and Precision Dose Escalation. Seminars in Radiation Oncology 2021; In press.
  4. 4.↵
    Gensheimer MF, Hong JC, Chang-Halpenny C, et al. Mid-radiotherapy PET/CT for prognostication and detection of early progression in patients with stage III non-small cell lung cancer. Radiotherapy and oncology: journal of the European Society for Therapeutic Radiology and Oncology 2017; 125(2): 338–43.
    OpenUrl
  5. 5.
    Yossi S, Krhili S, Muratet JP, Septans AL, Campion L, Denis F. Early assessment of metabolic response by 18F-FDG PET during concomitant radiochemotherapy of non-small cell lung carcinoma is associated with survival: a retrospective single-center study. Clin Nucl Med 2015; 40(4): e215–21.
    OpenUrlCrossRefPubMed
  6. 6.
    Huang W, Zhou T, Ma L, et al. Standard uptake value and metabolic tumor volume of (1)(8)F-FDG PET/CT predict short-term outcome early in the course of chemoradiotherapy in advanced non-small cell lung cancer. European journal of nuclear medicine and molecular imaging 2011; 38(9): 1628–35.
    OpenUrlCrossRefPubMed
  7. 7.
    Kong FS, Li L, Wang W, et al. Greater reduction in mid-treatment FDG-PET volume may be associated with worse survival in non-small cell lung cancer. Radiotherapy and oncology: journal of the European Society for Therapeutic Radiology and Oncology 2019; 132: 241–9.
    OpenUrl
  8. 8.↵
    Brodin NP, Tome WA, Abraham T, Ohri N. (18)F-Fluorodeoxyglucose PET in Locally Advanced Non-small Cell Lung Cancer: From Predicting Outcomes to Guiding Therapy. PET Clin 2020; 15(1): 55–63.
    OpenUrlCrossRef
  9. 9.↵
    Kim N, Kim JS, Geol Lee C. Predictive value of interim 18F-FDG-PET in patients with non-small cell lung cancer treated with definitive radiation therapy. PloS one 2020; 15(7): e0236350.
    OpenUrl
  10. 10.↵
    Chen HHW, Su WC, Guo HR, Lee BF, Chiu NT. Prognostic value of volumetric metabolic parameter changes determined by during and after radiotherapy-based (18) F-FDG PET/CT in stage III non-small cell lung cancer. The Kaohsiung journal of medical sciences 2019; 35(3): 151–9.
    OpenUrl
  11. 11.
    Chen HH, Lee BF, Su WC, et al. The increment in standardized uptake value determined using dual-phase 18F-FDG PET is a promising prognostic factor in non-small-cell lung cancer. European journal of nuclear medicine and molecular imaging 2013; 40(10): 1478–85.
    OpenUrlPubMed
  12. 12.↵
    Chen HH, Chiu NT, Su WC, Guo HR, Lee BF. Prognostic value of whole-body total lesion glycolysis at pretreatment FDG PET/CT in non-small cell lung cancer. Radiology 2012; 264(2): 559–66.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Ohri N, Piperdi B, Garg MK, et al. Pre-treatment FDG-PET predicts the site of in-field progression following concurrent chemoradiotherapy for stage III non-small cell lung cancer. Lung cancer 2015; 87(1): 23–7.
    OpenUrl
  14. 14.↵
    Na F, Wang J, Li C, Deng L, Xue J, Lu Y. Primary tumor standardized uptake value measured on F18-Fluorodeoxyglucose positron emission tomography is of prediction value for survival and local control in non-small-cell lung cancer receiving radiotherapy: meta-analysis. Journal of thoracic oncology: official publication of the International Association for the Study of Lung Cancer 2014; 9(6): 834–42.
    OpenUrl
  15. 15.↵
    Bradley JD, Paulus R, Komaki R, et al. Standard-dose versus high-dose conformal radiotherapy with concurrent and consolidation carboplatin plus paclitaxel with or without cetuximab for patients with stage IIIA or IIIB non-small-cell lung cancer (RTOG 0617): a randomised, two-by-two factorial phase 3 study. The Lancet Oncology 2015; 16(2): 187–99.
    OpenUrlCrossRefPubMed
  16. 16.↵
    Bradley JD, Paulus R, Graham MV, et al. Phase II trial of postoperative adjuvant paclitaxel/carboplatin and thoracic radiotherapy in resected stage II and IIIA non-small-cell lung cancer: Promising long-term results of the radiation therapy oncology group - RTOG 9705. Journal of Clinical Oncology 2005; 23(15): 3480–7.
    OpenUrlAbstract/FREE Full Text
  17. 17.
    Machtay M, Bae K, Movsas B, et al. Higher biologically effective dose of radiotherapy is associated with improved outcomes for locally advanced non-small cell lung carcinoma treated with chemoradiation: an analysis of the Radiation Therapy Oncology Group. International journal of radiation oncology, biology, physics 2012; 82(1): 425–34.
    OpenUrlCrossRefPubMed
  18. 18.↵
    van Diessen J, De Ruysscher D, Sonke JJ, et al. The acute and late toxicity results of a randomized phase II dose-escalation trial in non-small cell lung cancer (PET-boost trial). Radiotherapy and oncology: journal of the European Society for Therapeutic Radiology and Oncology 2019; 131: 166–73.
    OpenUrl
  19. 19.↵
    Raman S, Bissonnette JP, Warner A, et al. Rationale and Protocol for a Canadian Multicenter Phase II Randomized Trial Assessing Selective Metabolically Adaptive Radiation Dose Escalation in Locally Advanced Non-small-cell Lung Cancer (NCT02788461). Clinical lung cancer 2018; 19(5): e699–e703.
    OpenUrl
  20. 20.↵
    Kong FM, Ten Haken RK, Schipper M, et al. Effect of Midtreatment PET/CT-Adapted Radiation Therapy With Concurrent Chemotherapy in Patients With Locally Advanced Non-Small-Cell Lung Cancer: A Phase 2 Clinical Trial. JAMA oncology 2017; 3(10): 1358–65.
    OpenUrl
  21. 21.↵
    Lee E, Zeng J, Miyaoka RS, et al. Functional lung avoidance and response-adaptive escalation (FLARE) RT: Multimodality plan dosimetry of a precision radiation oncology strategy. Medical physics 2017; 44(7): 3418–29.
    OpenUrl
  22. 22.↵
    Lee Y, Auh SL, Wang Y, et al. Therapeutic effects of ablative radiation on local tumor require CD8+ T cells: changing strategies for cancer treatment. Blood 2009; 114(3): 589–95.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Demaria S, Formenti SC. Role of T lymphocytes in tumor response to radiotherapy. Frontiers in oncology 2012; 2: 95.
    OpenUrl
  24. 24.↵
    Aristophanous M, Penney BC, Pelizzari CA. The development and testing of a digital PET phantom for the evaluation of tumor volume segmentation techniques. Medical physics 2008; 35(7): 3331–42.
    OpenUrlPubMed
  25. 25.↵
    Werner-Wasik M, Nelson AD, Choi W, et al. What is the best way to contour lung tumors on PET scans? Multiobserver validation of a gradient-based method using a NSCLC digital PET phantom. International journal of radiation oncology, biology, physics 2012; 82(3): 1164–71.
    OpenUrlCrossRefPubMed
  26. 26.↵
    O JH, Lodge MA, Wahl RL. Practical PERCIST: A Simplified Guide to PET Response Criteria in Solid Tumors 1.0. Radiology 2016; 280(2): 576–84.
    OpenUrl
  27. 27.↵
    Wahl RL, Jacene H, Kasamon Y, Lodge MA. From RECIST to PERCIST: Evolving Considerations for PET response criteria in solid tumors. J Nucl Med 2009; 50 Suppl 1: 122S–50S.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    Haricharan S, Bainbridge MN, Scheet P, Brown PH. Somatic mutation load of estrogen receptor-positive breast tumors predicts overall survival: an analysis of genome sequence data. Breast cancer research and treatment 2014; 146(1): 211–20.
    OpenUrlCrossRefPubMed
  29. 29.↵
    Wood RD, Mitchell M, Sgouros J, Lindahl T. Human DNA repair genes. Science 2001; 291(5507): 1284–9.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    Hwang S, Kwon AY, Jeong JY, et al. Immune gene signatures for predicting durable clinical benefit of anti-PD-1 immunotherapy in patients with non-small cell lung cancer. Scientific reports 2020; 10(1): 643.
    OpenUrl
  31. 31.
    Patel SJ, Sanjana NE, Kishton RJ, et al. Identification of essential genes for cancer immunotherapy. Nature 2017; 548(7669): 537–42.
    OpenUrlCrossRefPubMed
  32. 32.↵
    Hugo W, Zaretsky JM, Sun L, et al. Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic Melanoma. Cell 2016; 165(1): 35–44.
    OpenUrlCrossRefPubMed
  33. 33.↵
    Chen HY, Yu SL, Chen CH, et al. A five-gene signature and clinical outcome in non-small-cell lung cancer. The New England journal of medicine 2007; 356(1): 11–20.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    Zhu CQ, Ding K, Strumpf D, et al. Prognostic and predictive gene signature for adjuvant chemotherapy in resected non-small-cell lung cancer. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 2010; 28(29): 4417–24.
    OpenUrl
  35. 35.↵
    Noth I, Zhang Y, Ma SF, et al. Genetic variants associated with idiopathic pulmonary fibrosis susceptibility and mortality: a genome-wide association study. The Lancet Respiratory medicine 2013; 1(4): 309–17.
    OpenUrl
  36. 36.↵
    Fingerlin TE, Murphy E, Zhang W, et al. Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis. Nature genetics 2013; 45(6): 613–20.
    OpenUrlCrossRefPubMed
  37. 37.↵
    Chiu CH, Wang YT, Walther BA, Chao A. An improved nonparametric lower bound of species richness via a modified good-turing frequency formula. Biometrics 2014; 70(3): 671–82.
    OpenUrlCrossRefPubMed
  38. 38.↵
    EFRON B, THISTED R. Estimating the number of unseen species: How many words did Shakespeare know? Biometrika 1976; 63(3): 435–47.
    OpenUrlCrossRefWeb of Science
  39. 39.↵
    Daley T, Smith AD. Predicting the molecular complexity of sequencing libraries. Nat Methods 2013; 10(4): 325–7.
    OpenUrlCrossRefPubMed
  40. 40.↵
    Pielou EC. The measurement of diversity in different types of biological collections. Journal of Theoretical Biology 1966; 13: 131–44.
    OpenUrlCrossRefWeb of Science
  41. 41.↵
    Simpson EH. Measurement of Diversity. Nature 1949; 163(4148): 688-.
    OpenUrlCrossRefWeb of Science
  42. 42.↵
    DeWolf S, Grinshpun B, Savage T, et al. Quantifying size and diversity of the human T cell alloresponse. JCI Insight 2018; 3(15).
  43. 43.↵
    van Elmpt W, Ollers M, Dingemans AM, Lambin P, De Ruysscher D. Response assessment using 18F-FDG PET early in the course of radiotherapy correlates with survival in advanced-stage non-small cell lung cancer. Journal of nuclear medicine: official publication, Society of Nuclear Medicine 2012; 53(10): 1514–20.
    OpenUrl
  44. 44.↵
    Weichselbaum RR, Liang H, Deng L, Fu YX. Radiotherapy and immunotherapy: a beneficial liaison? Nature reviews Clinical oncology 2017; 14(6): 365–79.
    OpenUrl
  45. 45.↵
    Golden EB, Apetoh L. Radiotherapy and immunogenic cell death. Semin Radiat Oncol 2015; 25(1): 11–7.
    OpenUrlCrossRefPubMed
  46. 46.↵
    Shahamatdar S, He MX, Reyna MA, et al. Germline Features Associated with Immune Infiltration in Solid Tumors. Cell reports 2020; 30(9): 2900–8 e4.
    OpenUrl
  47. 47.↵
    Kansy BA, Shayan G, Jie HB, et al. T cell receptor richness in peripheral blood increases after cetuximab therapy and correlates with therapeutic response. Oncoimmunology 2018; 7(11): e1494112.
    OpenUrl
  48. 48.
    Postow MA, Manuel M, Wong P, et al. Peripheral T cell receptor diversity is associated with clinical outcomes following ipilimumab treatment in metastatic melanoma. Journal for immunotherapy of cancer 2015; 3: 23.
    OpenUrl
  49. 49.↵
    Charles J, Mouret S, Challende I, et al. T-cell receptor diversity as a prognostic biomarker in melanoma patients. Pigment cell & melanoma research 2020; 33(4): 612–24.
    OpenUrl
Back to top
PreviousNext
Posted February 23, 2021.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Prognostic value and peripheral immunologic correlates of early FDG PET response imaging in a phase II trial of risk-adaptive chemoradiation for unresectable non-small cell lung cancer
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Prognostic value and peripheral immunologic correlates of early FDG PET response imaging in a phase II trial of risk-adaptive chemoradiation for unresectable non-small cell lung cancer
Stephen R. Bowen, Daniel S. Hippe, Hannah M. Thomas, Balukrishna Sasidharan, Paul D. Lampe, Christina S. Baik, Keith D. Eaton, Sylvia Lee, Renato G. Martins, Rafael Santana-Davila, Delphine Chen, Paul E. Kinahan, Robert S. Miyaoka, Hubert J. Vesselle, A. McGarry Houghton, Ramesh Rengan, Jing Zeng
medRxiv 2021.02.22.21251919; doi: https://doi.org/10.1101/2021.02.22.21251919
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Prognostic value and peripheral immunologic correlates of early FDG PET response imaging in a phase II trial of risk-adaptive chemoradiation for unresectable non-small cell lung cancer
Stephen R. Bowen, Daniel S. Hippe, Hannah M. Thomas, Balukrishna Sasidharan, Paul D. Lampe, Christina S. Baik, Keith D. Eaton, Sylvia Lee, Renato G. Martins, Rafael Santana-Davila, Delphine Chen, Paul E. Kinahan, Robert S. Miyaoka, Hubert J. Vesselle, A. McGarry Houghton, Ramesh Rengan, Jing Zeng
medRxiv 2021.02.22.21251919; doi: https://doi.org/10.1101/2021.02.22.21251919

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Oncology
Subject Areas
All Articles
  • Addiction Medicine (175)
  • Allergy and Immunology (421)
  • Anesthesia (97)
  • Cardiovascular Medicine (901)
  • Dentistry and Oral Medicine (171)
  • Dermatology (102)
  • Emergency Medicine (257)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (407)
  • Epidemiology (8796)
  • Forensic Medicine (4)
  • Gastroenterology (405)
  • Genetic and Genomic Medicine (1864)
  • Geriatric Medicine (179)
  • Health Economics (388)
  • Health Informatics (1292)
  • Health Policy (644)
  • Health Systems and Quality Improvement (493)
  • Hematology (207)
  • HIV/AIDS (396)
  • Infectious Diseases (except HIV/AIDS) (10577)
  • Intensive Care and Critical Care Medicine (564)
  • Medical Education (193)
  • Medical Ethics (52)
  • Nephrology (218)
  • Neurology (1762)
  • Nursing (104)
  • Nutrition (267)
  • Obstetrics and Gynecology (343)
  • Occupational and Environmental Health (461)
  • Oncology (966)
  • Ophthalmology (283)
  • Orthopedics (107)
  • Otolaryngology (177)
  • Pain Medicine (118)
  • Palliative Medicine (43)
  • Pathology (265)
  • Pediatrics (557)
  • Pharmacology and Therapeutics (266)
  • Primary Care Research (221)
  • Psychiatry and Clinical Psychology (1846)
  • Public and Global Health (3992)
  • Radiology and Imaging (655)
  • Rehabilitation Medicine and Physical Therapy (344)
  • Respiratory Medicine (536)
  • Rheumatology (216)
  • Sexual and Reproductive Health (178)
  • Sports Medicine (166)
  • Surgery (197)
  • Toxicology (37)
  • Transplantation (107)
  • Urology (80)