Purpose:

Chemoradiotherapy for locally advanced cervical cancer fails in over a third of patients. Biomarkers with therapeutic implications are therefore needed. We investigated the relationship between an established prognostic marker, maximum standardized uptake value (SUVmax) on 18F-fluorodeoxyglucose positron emission tomography, and the inflammatory and immune state of cervical cancers.

Experimental Design:

An SUVmax most prognostic for freedom from progression (FFP) was identified and compared with known prognostic clinical variables in a cohort of 318 patients treated with definitive radiation with prospectively collected clinical data. Gene set enrichment analysis (GSEA) and CIBERSORT of whole-transcriptome data from 68 patients were used to identify biological pathways and immune cell subpopulations associated with high SUVmax. IHC using a tissue microarray (TMA, N = 82) was used to validate the CIBERSORT findings. The impact of macrophages on cervical cancer glucose metabolism was investigated in coculture experiments.

Results:

SUVmax <11.4 was most prognostic for FFP (P = 0.001). The GSEA showed that high SUVmax is associated with increased gene expression of inflammatory pathways, including JAK/STAT3 signaling. CIBERSORT and CD68 staining of the TMA showed high SUVmax tumors are characterized by a monocyte-predominant immune infiltrate. Coculture of cervical cancer cells with macrophages or macrophage-conditioned media altered glucose uptake, and IL6 and JAK/STAT3 signaling contribute to this effect.

Conclusions:

SUVmax is a prognostic marker in cervical cancer that is associated with activation of inflammatory pathways and tumor infiltration of myeloid-derived immune cells, particularly macrophages. Macrophages contribute to changes in cervical cancer glucose metabolism.

See related commentary by Williamson et al., p. 4136

Translational Relevance

18F-fluorodeoxyglocose pretreatment positron emission tomography (FDG-PET) is a valuable tool for evaluating cervical cancer, and maximum standardized uptake value (SUVmax) is a widely clinically reported metric. Previous series have shown SUVmax to be prognostic. In the current study, we confirm the prognostic value of SUVmax using prospectively collected clinical outcome data for 318 patients treated with standard-of-care radiotherapy. Using smaller cohorts, we identify gene signatures and immune cell populations associated with a high SUVmax using whole-transcriptome analysis and IHC analysis of a tissue microarray. We show high SUVmax is associated with inflammatory gene signatures (e.g., JAK/STAT3 signaling) and myeloid-derived immune infiltrates, particularly macrophages. We then demonstrate that coculture of macrophages with cervical cancer cells alters tumor JAK/STAT3 signaling, and this in turn increases glucose uptake. Therefore, SUVmax could be used to identify cancers susceptible to therapies targeting select inflammatory pathways and infiltration by tumor permissive immune subpopulations.

Despite significant advances in screening and prevention, cervical cancer remains a leading cause of cancer deaths in women worldwide (1). The standard treatment for locally advanced cervical cancer, definitive pelvic radiotherapy (RT) given with concurrent cisplatin chemotherapy, has not been improved upon for nearly 20 years, despite the fact that this therapy fails in approximately one third of patients (2–4). Treatments for recurrent and metastatic cervical cancer are limited, associated with significant morbidity, and rarely curative (5, 6). Therefore, innovative treatment strategies are needed.

A number of prognostic imaging markers have been developed for cervical cancer, including metrics derived from positron emission tomography with 18F-fluorodeoxyglucose (FDG-PET) and magnetic resonance imaging (MRI; refs. 7–11). We have previously demonstrated that persistent or progressive disease on posttreatment FDG-PET is a powerful negative prognostic marker, and that it is associated with PI3K/AKT pathway activation (8, 12). However, this is of limited value when making initial treatment decisions, thus a pretreatment marker is desirable. Pretreatment maximum standardized uptake value (SUVmax) of the primary tumor is particularly suited for this purpose. Other potentially more robust PET metrics have been proposed, but SUVmax is the most widely clinically reported, most widely studied, and we have previously reported its prognostic significance in patients with cervical cancer (7). A deeper knowledge of primary tumor biology associated with a high SUVmax would further our understanding of how high-risk cancers should be treated, and could transform FDG-PET signal from a prognostic marker to a marker predictive of response to specific therapies.

The purpose of this work was to identify targetable signaling pathways, genetic signatures, and immune cell populations enriched in cervical cancers with a high pretreatment SUVmax. We established an SUVmax cutoff most prognostic for freedom from progression (FFP) in a cohort of patients with cervical cancer treated with definitive RT using modern techniques to identify high-risk cancers, confirming our prior work that established SUVmax as prognostic in a new cohort of patients (7). We used whole-transcriptome RNA-seq data from a subcohort of patients to identify signaling pathways associated with high SUVmax using gene set enrichment analysis (GSEA). We likewise used the RNA-seq data to identify immune cell populations associated with high SUVmax using CIBERSORT, and confirmed these results using a tissue microarray (TMA). Finally, we performed benchtop studies to further explore the impact of macrophages on cervical cancer glucose metabolism and the signaling pathways identified in the GSEA.

Patient cohorts

SUVmax cohort

All studies involving patient records and samples were conducted in accordance with the Declaration of Helsinki. An SUVmax cutoff with the most prognostic value was determined in a cohort of patients with cervical cancer treated with definitive RT at our institution between 2006 and 2015. To be included, patients had to have completed the entire prescribed course of definitive-intent RT. FIGO stage I–IVa tumors (FIGO 2009 staging) and all histologic subtypes were allowed, patients with metastatic disease were excluded, all patients were >18 years old at the time of their diagnosis, and all patients underwent pretreatment FDG-PET/computed tomography (FDG-PET/CT) or FDG-PET/magnetic resonance imaging (FDG-PET/MRI). Details of our institutional practice for treating cervical cancer with RT have previously been reported (13). Briefly, patients were treated with external beam RT to the pelvis using intensity-modulated RT and to the para-aortic lymph nodes if they were involved. The primary cervical tumor was also treated with image-guided brachytherapy. Nearly all (90.1%) patients were treated with concurrent chemotherapy. This cohort of patients was analyzed with institutional review board (IRB) approval with waiver of consent (IRB No. 201601052).

Whole-transcriptome RNA-sequencing (RNA-seq) cohort

The second cohort of patients was used to identify overexpressed genes and pathways associated with high FDG uptake. This consisted of patients enrolled on a prospective tumor banking study with written informed consent (IRB No. 201105374). Overlap between this cohort and the SUVmax cohort is described in the results below and in Supplementary Fig. S1. Tumor samples with sufficient high-quality RNA for whole-transcriptome sequencing (RNA-seq), as defined by the criteria utilized for The Cancer Genome Atlas (TCGA), were included in this study (14). Tumors were not macrodissected and thus results are expected to include genes attributable to tumor and any infiltrating stromal or immune cells.

TMA IHC cohort

The third cohort of patients consisted of patients whose tumor samples were incorporated into a TMA, and were interpretable for IHC analysis. Overlap between this cohort and the SUVmax cohort is also described in the results below and in Supplementary Fig. S1. This was also done with a waiver of consent from the IRB for access to tissue specimens and clinical data (IRB No. 201807063). Cohort 3 was also composed of patients treated with definitive RT.

PET/CT imaging and analysis

FDG-PET was performed as part of routine staging evaluation of patients with newly diagnosed cervical cancer according to our institution's imaging protocols, which have been developed to maximize visualization of cervical tumors (15). Briefly, patients were required to fast for at least 4 hours prior to FDG injection, and blood glucose levels were required to be less than 200 mg/dL at the time of tracer injection (5/331 patients had a blood glucose over 200 mg/dL). PET/CT images were acquired approximately 60 minutes after FDG injection. Images were acquired on 1 of 3 PET/CT scanners: Siemens Biograph 2 (n = 107), Siemens Biograph 40-HD (n = 187), and Siemens mCT (n = 32), and on one PET/MRI scanner (Siemens mMR, n = 5). All images were reconstructed using standard iterative reconstruction algorithms.

As part of routine clinical practice, the SUVmax of the primary cervical tumor was measured by determining the maximum single-voxel SUV within a volume of interest encompassing the primary tumor. A representative image of a large cervical tumor with an automatically generated metabolic tumor volume is shown in Supplementary Fig. S2. The SUV is defined as: SUVmax = (tissue radioactivity concentration (mCi/mL) × patient weight (g))/injected dose (mCi).

Statistical analysis of clinical data

Given the significant overlap between the cohort of patients used in the SUV analysis and the cohort of patients used for the RNA-seq analysis, a comparison was made between patients included only in the SUV analysis cohort and patients included in both cohorts using the Wilcoxon two-sample test and Fisher exact test, as appropriate.

The optimal SUVmax for discriminating patients with good versus poor FFP was determined by outcome-oriented cutoff methodology utilizing an alpha-adjusted log-likelihood test (16). Once the optimal SUVmax was determined, the prognostic significance of this SUVmax for FFP, cancer-specific survival (CSS), and overall survival (OS) was evaluated with Kaplan–Meier/log-rank and Cox univariate and multivariate analysis. For FFP any recurrence was counted as an event and patients were censored at time of last follow-up. For CSS death from cervical cancer was counted as an event, and for OS, death from any cause was counted as an event. In either case, patients were censored at time of last follow-up or death.

A bootstrap Cox multivariate analysis was then performed to better illustrate the significance of SUVmax in the context of other known clinical prognostic markers. The bootstrap analysis was performed by taking 1,000 random samplings from the first cohort of patients, such that each sampling was half the size of the total cohort. A Cox multivariate analysis was then run on each sampling. The variables included in this analysis were SUVmax, FIGO stage, lymph node involvement by FDG-PET (none, pelvic, para-aortic), age, and histology (squamous vs. nonsquamous).

GSEA

We performed RNA-seq analysis on the patients in cohort 2. Specific details regarding this are included in the Supplementary Methods. The expressed genes were ranked according to fold change between patients with high SUVmax (≥11.41) and low SUVmax (<11.41) values and included to perform GSEA using default parameters (17–21). The H: Hallmark gene sets (n = 50) and the C2: curated gene sets (n = 4,762) from the Molecular Signatures Database (MSigDB) were used to discover altered biological pathways between patients with cervical cancer with an SUVmax above or below the cutoff point established in the SUVmax analysis (21).

Cell-type abundance analysis

Using CIBERSORT, we estimated the relative abundances of hematopoietic cells in the patients from cohort 2 (22). For signature genes, we used the leukocyte gene signature matrix, LM22, which contains 547 genes that distinguish 22 human hematopoietic cell phenotypes, including seven T-cell types, naïve and memory B cells, plasma cells, natural killer (NK) cells, and myeloid subsets. Relative fractions of B cells, CD4 T cells, CD8 T cells, neutrophils (M0, M1, and M2), macrophages, monocytes, dendritic cells, and NK cells were compared between patients with high and low SUVmax using the Student t test.

Prevalence of immune cell populations in cervical tumors using a tumor microarray

To evaluate the prevalence of various immune cells in human cervical tumors, we performed IHC for markers of macrophages and neutrophils on a TMA of human cervical cancer tumors. Details regarding this TMA have previously been published (23), and are outlined in more detail in the Supplementary Methods. Briefly, areas of cancer on formalin-fixed paraffin-embedded (FFPE) tumor specimens were identified and triplicate 1.5-mm cores were placed in a TMA. Sections 6 μm thick were then mounted onto polarized slides for IHC staining. Staining for macrophages was done on a Ventana automated stainer (Ventana Medical Systems). Staining was performed with anti-CD68 (0.4 μg/mL, monoclonal mouse, Ventana) and anti-CD163 (0.3 μg/mL, monoclonal mouse, Cell Marque). Neutrophil staining was performed manually, and slides were costained with CD66b (1:600, BD Biosciences) and neutrophil elastase (1:200, Novus).

Stained slides were digitally scanned and analyzed using the Aperio eSlide Manager platform (Leica Biosystems). Samples were considered either negative (≤1% of cells stained) or positive (>1% of cells) for all IHC stains used. Categorizing immune infiltrate as high versus low has previously been described, including for CD68 (24–26). Samples were scored according to immune cell infiltrate into the tumor. All tumor samples were scored by two independent blinded readers. Readers evaluated cores at 20× magnification and scored each core as either negative (≤1% of cells stained) or positive (>1% of cells stained). Agreement between readers was assessed with Cohen's kappa statistic. The scores from both readers for all available cores from each patient were averaged to generate a composite score rounded to the nearest whole number (<0.5 to 0 for negative, ≥0.5 to 1 for positive). This score was matched to the SUVmax for that tumor, and the SUVmax values of tumors that stained positive versus negative were compared using the Wilcoxon two-sample test.

Tissue culture and reagents

Two human cervical cancer cell lines obtained through the ATCC collection were evaluated: ME180 and SiHa. Experiments were performed on cell lines between passages 10 and 30. Mycoplasma testing was performed by the iPSC Center (ME180 and SiHa cells, February 20, 2019) at Washington University School of Medicine. All experiments were carried out under standard tissue culture conditions at 37°C and 5% CO2, with cells maintained in IMDM with 10% fetal bovine serum and 0.1 mg/mL gentamycin.

Macrophage coculture and conditioned media experiments

Human monocytic THP-1 cells were differentiated into M0, M1, and M2 macrophages using methods outlined in the Supplementary Methods. The generation of macrophage-conditioned media is also outlined in the Supplementary Methods.

In the coculture experiments, SiHa and ME180 cells were plated onto 6- or 12-well plates, and THP-1 monocytes were differentiated into M0, M1, or M1 macrophages in 6 or 12 Transwell inserts (membrane pore size of 0.4 μm, Corning). Cells were cocultured for 48 hours. For the conditioned media experiments, cervical cancer cells were maintained as above and directly exposed to conditioned media collected from different macrophages. Cervical cancer cells from both conditioned media and coculture experiments were harvested after 48 hours for Western blot analysis or a glucose uptake assay.

Western blotting

Activation of STAT3 by coculture with macrophages was confirmed using Western blotting with primary antibodies against total STAT3 (1:1,000, room temperature for 1 hour, Santa Cruz Biotech) and phosphorylated STAT3 (pSTAT3 (Tyr705), 1:2,000, 4°C overnight, Cell Signaling Technology). Blots were probed with horseradish peroxidase–conjugated anti-rabbit (1:3,000, Cell Signaling Technology) for 1 hour at room temperature. Blotting was also performed for actin (1:1,000; Santa Cruz Biotech) at room temperature for 1 hour. The STAT3 and actin primary antibodies were both conjugated to HRP. Whenever necessary, blots were stripped by incubating membrane with restore Western blot stripping buffer (Thermo scientific) with intermittent shaking and reprobed with required antibodies. For detection, Amersham ECL select (GE Healthcare) was used according to the manufacturer's protocol. Images were acquired using Chemidoc Imaging systems (BioRad).

Fluorescent 2-DG glucose uptake assay

A glucose uptake assay was performed using 2-[N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino]-2-deoxy-D-glucose (2-NBDG; Cayman) by flow-cytometric detection of fluorescence produced by the cells. Briefly, cervical cancer cells were plated in a 12-well plate and cocultured with macrophages or conditioned media collected from macrophages as described above. After 48 hours of exposure, cells were trypsinized, washed in PBS and stained by 2-NBDG. For 2-NBDG uptake in the experimental settings of α-IL6, cells were exposed with α-IL6 (500 μg/mL). After 48 hours of exposure, cells were washed three times with PBS to remove remaining α-IL6 from the culture media. Subsequently, cells were trypsinized, washed in PBS, and stained by 2-NBDG. Flow cytometry was performed using a MACSQuant Analyzer10 flow cytometer (Miltenyi Biotec). For each measurement, data from 10,000 single-cell events were collected using FACScan flow cytometer.

Cytokine array analysis

Cytokines were analyzed by human Discovery Assay using an array of 13 cytokines (Human Cytokine Array/Chemokine Array 13-Plex, Eve Technologies Corporation). For each sample, conditioned or culture media were collected from macrophages cocultured with ME180 cells as a supernatant and processed according to the company's specifications and requirements.

Patients

The cohort of patients used in the SUVmax analysis (cohort 1) consisted of 318 patients, the cohort of patients used for the RNA-seq analysis (cohort 2) consisted of 68 patients, and the cohort used for the TMA consisted of 82 patients. All but eight of the patients in cohort 2 were also included in cohort 1, and all but 11 of the patients in cohort 3 were also included in cohort 1. Patient characteristics are summarized in Table 1, and the overlap between the cohorts is summarized in Supplementary Fig. S1. We compared the characteristics of the patients included only in the SUVmax analysis with the patients included in both the SUVmax analysis and RNA-seq analysis, and also with the patients included in both the SUVmax analysis and the TMA analysis. There were no significant differences between these groups in terms of age at diagnosis, FIGO stage, lymph node involvement by FDG-PET, histology, or the primary tumor SUVmax (Table 1).

Table 1.

Patient characteristics of the three study cohorts (SUVmax, RNA-seq, and TMA).

SUV cohortRNA-seq CohortP valueTMA cohortP value
All patients 318 68 (57 in SUV cohort)  82 (67 in SUV cohort)  
Age at diagnosis   0.42  0.94 
Median (range) 50 (23–92) 53 (25–81)  52 (25–92)  
FIGO stage (%)   0.52  0.26 
125 (39.3%) 20 (29.4%)  29 (35.4%)  
II 119 (37.4%) 25 (36.8%)  28 (34.1%)  
III 69 (21.7%) 21 (30.9%)  22 (26.8%)  
IV 5 (1.6%) 2 (2.9%)  3 (3.7%)  
Lymph node involvement (%)   0.08  0.82 
None 146 (45.9%) 33 (48.5%)  40 (58.8%)  
Pelvic 122 (38.4%) 21 (30.9%)  27 (32.9%)  
Para-aortic 50 (15.7%) 14 (20.6%)  15 (18.3%)  
Histology   0.85  0.86 
Squamous 256 (80.5%) 55 (80.9%)  56 (80.5%)  
Nonsquamous 62 (19.5%) 13 (19.1%)  16 (19.5%)  
  Adenocarcinoma 42 (13.2%) 9 (13.2%)  10 (12.2%)  
  Small cell 8 (2.5%) 2 (2.9%)  2 (2.4%)  
  Adenosquamous 4 (1.3%) 2 (2.9%)  1 (1.2%)  
  Other 8 (2.5%)  3 (3.7%)  
SUVmax   0.65  0.19 
Mean (range) 14.6 (2.9–60.3) 17.0 (3.4–53.6)  15.4 (3.77–60.3)  
SUV cohortRNA-seq CohortP valueTMA cohortP value
All patients 318 68 (57 in SUV cohort)  82 (67 in SUV cohort)  
Age at diagnosis   0.42  0.94 
Median (range) 50 (23–92) 53 (25–81)  52 (25–92)  
FIGO stage (%)   0.52  0.26 
125 (39.3%) 20 (29.4%)  29 (35.4%)  
II 119 (37.4%) 25 (36.8%)  28 (34.1%)  
III 69 (21.7%) 21 (30.9%)  22 (26.8%)  
IV 5 (1.6%) 2 (2.9%)  3 (3.7%)  
Lymph node involvement (%)   0.08  0.82 
None 146 (45.9%) 33 (48.5%)  40 (58.8%)  
Pelvic 122 (38.4%) 21 (30.9%)  27 (32.9%)  
Para-aortic 50 (15.7%) 14 (20.6%)  15 (18.3%)  
Histology   0.85  0.86 
Squamous 256 (80.5%) 55 (80.9%)  56 (80.5%)  
Nonsquamous 62 (19.5%) 13 (19.1%)  16 (19.5%)  
  Adenocarcinoma 42 (13.2%) 9 (13.2%)  10 (12.2%)  
  Small cell 8 (2.5%) 2 (2.9%)  2 (2.4%)  
  Adenosquamous 4 (1.3%) 2 (2.9%)  1 (1.2%)  
  Other 8 (2.5%)  3 (3.7%)  
SUVmax   0.65  0.19 
Mean (range) 14.6 (2.9–60.3) 17.0 (3.4–53.6)  15.4 (3.77–60.3)  

Note: A substantial number of the patients in the RNA-seq and TMA cohorts were also included in the SUVmax analysis. The P values shown represent comparisons between patients included in both the RNA-seq and SUVmax analyses versus patients included only in the SUVmax analysis to determine if the RNA-seq subgroup was biasing the overall cohort (using the Fisher exact or Wilcoxon tests as appropriate). The same was done for the TMA subgroup.

Median follow-up was 3.2 years for cohort 1, 4.0 years for cohort 2, and 5.1 years for cohort 3. OS rates at 2 and 5 years were 76.1% and 66.8%, respectively, in cohort 1; 72.5% and 56.5%, respectively, in cohort 2; and 69.5% and 56.1%, respectively, in cohort 3. CSS rates at 2 and 5 years were 77.6% and 70.7%, respectively, in cohort 1; 64.7% at both 2 and 5 years in cohort 2; and 73.2% and 61.0%, respectively, in cohort 3.

Prognostic value of SUVmax

Using outcome-oriented cutoff methodology, an SUVmax < 11.41 was found to be most prognostic of FFP. This is similar to the SUVmax of 13.3 previously reported as being associated with the worst survival (7). This SUVmax cutoff was also prognostic of CSS, but not OS (Fig. 1A–C), and was prognostic for local (cervical) recurrence and any pelvic recurrence as well (Supplementary Fig. S3). On Cox univariate analysis, SUVmax (≥ vs. <11.41), FIGO stage, and lymph node involvement on PET (para-aortic vs. none) were found to be prognostic for FFP. On multivariate analysis for FFP, SUVmax, FIGO stage, and lymph node involvement were all significant (Table 2).

Figure 1.

Prognostic significance of SUVmax. An SUVmax cutoff of 11.4 was prognostic for FFP and cause-specific survival (A and B) but not OS. C, P values represent results of the log-rank test.

Figure 1.

Prognostic significance of SUVmax. An SUVmax cutoff of 11.4 was prognostic for FFP and cause-specific survival (A and B) but not OS. C, P values represent results of the log-rank test.

Close modal
Table 2.

Cox univariate and multivariate analyses for FFP.

UnivariateMultivariate
VariableHR (95% CI)PHR (95% CI)P
SUVmax (≥ vs. <11.41) 2.20 (1.37–3.54) 0.001 1.73 (1.06–2.81) 0.03 
FIGO stage 
 II vs. I 1.28 (0.80–2.06) 0.30 1.07 (0.66–1.73) 0.79 
 III vs. I 2.83 (1.76–4.55) <0.0001 1.98 (1.20–3.27) 0.008 
 IV vs. I 4.75 (1.67–13.5) 0.003 2.30 (0.77–6.88) 0.14 
 III/IV vs. I/II 2.59 (1.76–3.82) <0.0001 1.89 (1.26–2.86) 0.002 
Lymph nodes 
 Pelvic vs. none 1.30 (0.83–2.03) 0.25 1.27 (0.81–2.00) 0.30 
 Para-aortic vs. none 3.93 (2.47–6.24) <0.0001 3.03 (1.88–4.88) <0.0001 
Age 1.00 (0.98–1.01) 0.54   
Histology (squamous vs. nonsquamous) 0.99 (0.61–1.59) 0.96   
UnivariateMultivariate
VariableHR (95% CI)PHR (95% CI)P
SUVmax (≥ vs. <11.41) 2.20 (1.37–3.54) 0.001 1.73 (1.06–2.81) 0.03 
FIGO stage 
 II vs. I 1.28 (0.80–2.06) 0.30 1.07 (0.66–1.73) 0.79 
 III vs. I 2.83 (1.76–4.55) <0.0001 1.98 (1.20–3.27) 0.008 
 IV vs. I 4.75 (1.67–13.5) 0.003 2.30 (0.77–6.88) 0.14 
 III/IV vs. I/II 2.59 (1.76–3.82) <0.0001 1.89 (1.26–2.86) 0.002 
Lymph nodes 
 Pelvic vs. none 1.30 (0.83–2.03) 0.25 1.27 (0.81–2.00) 0.30 
 Para-aortic vs. none 3.93 (2.47–6.24) <0.0001 3.03 (1.88–4.88) <0.0001 
Age 1.00 (0.98–1.01) 0.54   
Histology (squamous vs. nonsquamous) 0.99 (0.61–1.59) 0.96   

We also performed a bootstrap multivariate analysis to further determine the relative prognostic significance of SUVmax (≥ vs. <11.41) in the context of other known clinical prognostic variables; namely, age at diagnosis, FIGO stage, nodal involvement on pretreatment FDG-PET/CT, and histology. Results summarized in Supplementary Table S1 demonstrate that in 1,000 random samplings from cohort 1, SUVmax proved to be a significant prognostic factor 63.3% of the time, FIGO stage 72.1% of the time with FIGO stage III portending a particularly poor prognosis, nodal involvement 96.5% of the time with para-aortic lymph node involvement being a powerful prognostic variable, and age and histology 12.6 and 12.0% of the time, respectively. Median hazard ratios and P values when each variable was included in the model are shown in Supplementary Table S1.

Gene sets associated with high SUVmax

Once the association of SUVmax and prognosis was established, we performed a GSEA on the RNA-seq data to determine gene sets associated with a high SUVmax. The GSEA was performed treating SUVmax as a dichotomous variable (≥ or <11.41). We focused on Hallmark gene sets, which represent well-defined biological states that display coherent expression, and have been generated using a computational methodology (27). The top three Hallmark gene sets associated with high SUVmax were for inflammatory response, TNFα signaling via NF-κΒ, and JAK/STAT3 signaling (Fig. 2). There was overlap between many of the genes in these gene sets (i.e., gene present in >1 gene set). Furthermore, many of these genes encode cytokines and chemokines associated with specific immune cell subpopulations, including cells within the myeloid lineage, and macrophages in particular (Fig. 2).

Figure 2.

Results of the GSEA. Considering the Hallmark pathways (well-defined biological states), the top three pathways significantly upregulated in tumors with a high SUVmax were Inflammatory Response (A), TNFα signaling via NF-κΒ (B), and IL6 JAK STAT3 signaling (C). All remained significant accounting for FDR, and the Inflammatory Response and NF-κΒ gene sets also had a family-wise error rate (FWER) < 0.05, though JAK/STAT3 signaling did not. These gene sets shared a number of genes, many of which help regulate immune cells (bolded), particularly myeloid-derived cells. The genes displayed were shared by at least two of the gene sets (D).

Figure 2.

Results of the GSEA. Considering the Hallmark pathways (well-defined biological states), the top three pathways significantly upregulated in tumors with a high SUVmax were Inflammatory Response (A), TNFα signaling via NF-κΒ (B), and IL6 JAK STAT3 signaling (C). All remained significant accounting for FDR, and the Inflammatory Response and NF-κΒ gene sets also had a family-wise error rate (FWER) < 0.05, though JAK/STAT3 signaling did not. These gene sets shared a number of genes, many of which help regulate immune cells (bolded), particularly myeloid-derived cells. The genes displayed were shared by at least two of the gene sets (D).

Close modal

Computational prediction of immune cell populations associated with high SUVmax

Given the association between high SUVmax and inflammatory pathway genes related to myeloid cell populations, we performed a deconvolution of the bulk RNA-seq data using the CIBERSORT algorithm. The relative prevalence of neutrophils, monocytes, and NK cells was significantly different between the high and low SUVmax groups, with higher proportions of neutrophils and monocytes and lower proportions of NK cells in the high SUVmax group (Fig. 3A). The proportion of macrophages almost met the threshold for significance (P = 0.066; Fig. 3A). Given the number of genes related to macrophages overexpressed in the high SUVmax tumors in the GSEA, we also examined the predicted prevalence of macrophage subgroups, and found that a greater proportion of M2-polarized macrophages were predicted in high SUVmax tumors (Fig. 3B).

Figure 3.

Results of the CIBERSORT analysis and the TMA. CIBERSORT identified neutrophils and monocytes as significantly more prevalent, and NK cells as significantly less prevalent in the high SUVmax tumors (A). Further subdividing macrophages, CIBERSORT identified M2-polarized macrophages as more prevalent in the high SUVmax tumors (B). We stained a TMA of human cervical cancer tumors for CD68 (a marker for macrophages) and CD66b and neutrophil elastase (markers for neutrophils). Example slides representing negative and positive staining (with a range of positive staining) are shown for both CD68 (C) and CD66b and neutrophil elastase costaining (D). There was a significant difference in the average SUVmax between samples that stained positive (n = 62) for CD68 versus those that were negative (n = 14, P = 0.006 by the Wilcoxon two-sample test; C). There was no significant difference in the average SUVmax in samples that were positive (n = 38) versus negative (n = 26) for neutrophils (P = 0.24; D).

Figure 3.

Results of the CIBERSORT analysis and the TMA. CIBERSORT identified neutrophils and monocytes as significantly more prevalent, and NK cells as significantly less prevalent in the high SUVmax tumors (A). Further subdividing macrophages, CIBERSORT identified M2-polarized macrophages as more prevalent in the high SUVmax tumors (B). We stained a TMA of human cervical cancer tumors for CD68 (a marker for macrophages) and CD66b and neutrophil elastase (markers for neutrophils). Example slides representing negative and positive staining (with a range of positive staining) are shown for both CD68 (C) and CD66b and neutrophil elastase costaining (D). There was a significant difference in the average SUVmax between samples that stained positive (n = 62) for CD68 versus those that were negative (n = 14, P = 0.006 by the Wilcoxon two-sample test; C). There was no significant difference in the average SUVmax in samples that were positive (n = 38) versus negative (n = 26) for neutrophils (P = 0.24; D).

Close modal

Prevalence of macrophages and neutrophils in a cervical cancer TMA

To validate the findings from the GSEA and CIBERSORT data suggesting that high SUVmax is associated with inflammatory pathways and immune cells in the myeloid lineage, we performed IHC staining on a TMA of human cervical cancer tumors to identify immune cell subpopulations. We specifically focused on surface markers for neutrophils and macrophages. CD68 expression was used as a general marker of macrophages, with good agreement between the two readers (kappa = 0.69). Tumor samples with minimal/negative CD68 staining (≤1%) showed a lower average SUVmax than samples that were positive (>1%) for CD68 (Fig. 3C). CD163 expression was used as a marker of M2-polarized macrophages. There was no significant difference in the SUVmax between samples that were negative for CD163 versus those that were positive (Supplementary Fig. S4). To evaluate neutrophil infiltrate, separate TMA slides were costained with CD66b and neutrophil elastase; there was good agreement between readers of neutrophil staining (kappa = 0.70). There was a nonsignificant positive correlation between SUVmax and neutrophil infiltrate in the tumor microenvironment (Fig. 3D).

M0 and M2 macrophages alter glucose metabolism in human cervical cancer cells through inflammatory signaling pathways

To study the relationship between macrophages, glucose uptake, and JAK/STAT3 signaling, human monocytic THP-1 cells were differentiated into macrophages using phorbol 12-myristate 13-acetate (PMA), then further polarized into M1 macrophages using interferon-γ (INFγ) and lipopolysaccharide (LPS), or M2 macrophages using IL4 and IL13. Human cervical cancer cells (SiHa and ME180) were then exposed to macrophage-conditioned media or were cocultured with macrophages in transwell plates.

Glucose metabolism was measured using fluorescent 2-NBDG by flow cytometry. Exposing ME180 and SiHa cells to either M0, M1, or M2 macrophage-conditioned media, or coculturing ME180 and SiHa cells with M0, M1, or M2 macrophages, resulted in a significant increase in 2-NBDG uptake by the cervical cancer cells. This was seen in both cell lines with all subtypes of macrophages (Fig. 4A). We therefore evaluated an array of 13 proinflammatory cytokines and chemokines that were secreted in this transwell system (Eve Technologies). This demonstrated an increase in a wide array of cytokines secreted by macrophages. Two cytokines, IL1β and IL6, showed consistent increases when macrophages were cocultured with cervical cancer cells, compared with the macrophages alone (Fig. 4B). Results from four cytokines are not shown (INFγ, IL2, IL4, and IL13) because either the signal for most conditions was too low to be detected (INFγ and IL2), or the signal from ME180 cells alone outweighed that from the macrophages by several orders of magnitude (IL4 and IL13).

Figure 4.

Changes in glucose metabolism, cytokine and chemokine expression, and JAK/STAT3 signaling seen with coculture of cervical cancer cells and macrophages and macrophage-conditioned media. Coculture of macrophages with ME180 and SiHa cells, or culturing ME180 and SiHa cells with macrophage-conditioned media, results in a significant increase in 2-NBDG signal measured by flow cytometry; this is most pronounced in the SiHa cells (A). Differentiated macrophages produce a number of cytokines; coculture with ME180 cells significantly increased expression of IL1β and IL6 compared with differentiated macrophages alone (B). Coculture with differentiated macrophages, or culturing with macrophage-conditioned media, induces STAT3 phosphorylation in ME180 cells, which have little STAT3 phosphorylation at baseline. This effect was not seen in SiHa cells, which show high baseline levels of STAT3 phosphorylation (C). SiHa cells cultured with IL6 NAB do show reduced STAT3 phosphorylation (D). Treatment of macrophages cocultured with cervical cancer cells with IL6 NAB alters 2-NBDG uptake in both ME180 and SiHa cells. This effect was most pronounced with SiHa cells, and coculture of either SiHa or ME180 cells with M2 macrophages (E). GSEA showed that the Hallmark Glycolysis pathway was significantly upregulated in tumors with high STAT3 expression (F). Experiments were carried out in triplicate (n = 3), with error bars representing standard deviation.

Figure 4.

Changes in glucose metabolism, cytokine and chemokine expression, and JAK/STAT3 signaling seen with coculture of cervical cancer cells and macrophages and macrophage-conditioned media. Coculture of macrophages with ME180 and SiHa cells, or culturing ME180 and SiHa cells with macrophage-conditioned media, results in a significant increase in 2-NBDG signal measured by flow cytometry; this is most pronounced in the SiHa cells (A). Differentiated macrophages produce a number of cytokines; coculture with ME180 cells significantly increased expression of IL1β and IL6 compared with differentiated macrophages alone (B). Coculture with differentiated macrophages, or culturing with macrophage-conditioned media, induces STAT3 phosphorylation in ME180 cells, which have little STAT3 phosphorylation at baseline. This effect was not seen in SiHa cells, which show high baseline levels of STAT3 phosphorylation (C). SiHa cells cultured with IL6 NAB do show reduced STAT3 phosphorylation (D). Treatment of macrophages cocultured with cervical cancer cells with IL6 NAB alters 2-NBDG uptake in both ME180 and SiHa cells. This effect was most pronounced with SiHa cells, and coculture of either SiHa or ME180 cells with M2 macrophages (E). GSEA showed that the Hallmark Glycolysis pathway was significantly upregulated in tumors with high STAT3 expression (F). Experiments were carried out in triplicate (n = 3), with error bars representing standard deviation.

Close modal

We further explored the impact of IL6 on glucose metabolism in ME180 and SiHa cells. IL6 is known to activate the JAK/STAT3 pathway, resulting in STAT3 phosphorylation. ME180 cells showed no baseline STAT3 phosphorylation, but showed STAT3 phosphorylation when cultured with macrophages or macrophage-conditioned media. SiHa cells showed baseline STAT3 phosphorylation, which was unchanged or possibly slightly diminished when cocultured with macrophages and macrophage-conditioned media (Fig. 4C). Treatment with neutralizing antibody to IL6 (500 ng/mL) resulted in a decrease in STAT3 expression in SiHa cells (D), as well as a decrease in glucose uptake, as measured by 2-NBDG when either SiHa or ME180 cells were cocultured with M2 macrophages (Fig. 4E).

Finally, to connect JAK/STAT3 signaling to glucose metabolism in human tumors, we repeated the GSEA, but ranking expressed genes according to STAT3 expression as opposed to SUVmax. The Hallmark Glycolysis pathway was significantly upregulated in tumors with high STAT3 expression [false discovery rate (FDR) q <0.0001, FWER P <0.0001).

In this work, we have made a connection between a high SUVmax derived from FDG-PET, an immune infiltrate characterized by myeloid-derived cells, particularly monocytes and M2 macrophages, and inflammatory signaling pathways. As high SUVmax is also prognostic, this could serve as an imaging biomarker to identify high-risk cancers with a specific inflammatory and immune profile, and guide use of novel targeted or immune therapies for these patients. To our knowledge, these are the first data to demonstrate a connection between tumor glucose metabolism, inflammatory signaling pathways, and immune cell infiltrate in cervical cancer.

Our data specifically suggest that a high SUVmax is associated with inflammatory pathways, including JAK/STAT3 and NF-κΒ signaling (Fig. 2), both of which have previously been implicated as important signaling pathways in cervical cancer. Activated STAT3 (pSTAT3) has been shown to be a poor prognostic marker in cervical squamous cell carcinoma (28, 29). Likewise, NF-κΒ signaling has been implicated in cervical cancer carcinogenesis, progression, and therapeutic resistance, though the exact role of NF-κΒ in cervical cancer appears to be complex (30, 31). In addition, many of the genes identified in these pathways in our cohort have been shown to have prognostic significance or important biological implications in cervical cancer, including IL6 (32).

The inflammatory pathway gene sets associated with high SUVmax on the GSEA share a number of genes, particularly cytokines and chemokines linked with myeloid-derived immune cells (Fig. 2). The association between SUVmax and myeloid-derived immune cells was further supported by our CIBERSORT and IHC analyses (Fig. 3). The CIBERSORT algorithm showed an association between SUVmax and neutrophils, monocytes, and M2-polarized macrophages, and a negative association with NK cells, and the IHC analysis of the TMA further supported the association between CD68-positive cells (i.e., monocytes and macrophages) and SUVmax. An association between neutrophils and SUV was not seen on the TMA, though this remains an active area of investigation (33). The negative association between SUVmax and NK cells seen in the CIBERSORT analysis is interesting given the known high glucose uptake of NK cells (34). It is possible that the tumor microenvironment is modulating expression of NK cell genes included in the LM22 signature used for CIBERSORT (35). This study therefore provides the important insight that the immune microenvironment and tumor glycolysis are interconnected.

We have further demonstrated a connection between glucose metabolism, macrophage infiltrate, and JAK/STAT3 signaling (Fig. 4). Specifically, we have shown that coculturing human cervical cancer cells with macrophages or growing cervical cancer cells with macrophage-conditioned media results in increased glucose uptake. Cytokine array profiling showed IL6 and IL1β concentrations were significantly increased in media from macrophages cocultured with ME180 cells, compared with media from the macrophages alone, suggesting an interplay between the cervical cancer cells and the macrophages that alters expression of these factors. Coculture with macrophages also increased STAT3 phosphorylation in ME180 cells, and treatment with neutralizing antibody to IL6 decreased glucose uptake in SiHa and ME180 cells. Finally, we demonstrated that the Hallmark Glycolysis gene set is significantly upregulated in human tumors with high STAT3 expression, once again utilizing GSEA.

In aggregate, these data demonstrate: (i) primary tumor SUVmax is associated with FFP and CSS in women with cervical cancer; (ii) high SUVmax is associated with inflammatory signaling pathways; (iii) high SUVmax is also associated with a myeloid-derived immune cell infiltrate, particularly macrophages. We propose that the relationship between cervical cancer glucose metabolism and macrophage immune infiltrate is mediated in part by cytokines and chemokines involved in JAK/STAT3 signaling. There is some other evidence in the literature that STAT3 can affect glucose metabolism through modifying hexokinase expression, FOXO1 and HIF1α signaling, and mitochondrial metabolism, though further study is needed to establish the exact mechanism of this connection in cervical cancer (36–38). Notably, we did not establish a connection between SUVmax and the Hallmark Hypoxia signature in our GSEA.

The clinical implications of inflammatory pathways and immune makeup in cervical cancers are actively being investigated and elucidated. Tumor-associated macrophages have been shown to be prevalent in cervical tumors, more so than regulatory T cells, suggesting they contribute to an immunosuppressive tumor microenvironment (24, 39, 40). Tumor-associated macrophages have also been associated with poor prognosis in cervical cancers (26, 41). JAK/STAT3 is also a targetable pathway, and ruxolitinib is clinically used in blood disorders, and is under investigation in a number of cancers (NCT02713386, NCT02876302, and NCT03120624; refs. 42, 43). Our data provide a link between a proinflammatory but immunosuppressive microenvironment in cervical cancer, and a routinely clinically measured imaging metric, the PET-derived SUVmax for FDG.

This study has important limitations. We have focused on the SUVmax of cervical tumors treated at a single institution with a robust PET imaging program over a recent period of time (since 2006). There are a number of potential critiques regarding using SUVmax as a quantitative marker. In this case, our PET data were derived from several different scanners with differing performance characteristics, and the discriminating value for high versus low SUVmax in this study reflects this specific data set. We have not validated the reproducibility of SUVmax in cervical cancer across a broad array of scanners, reconstruction methods, and different imaging centers, and it is therefore not clear how easily it could be adopted as a quantitative predictive biomarker across many institutions (44). Future studies could alternatively utilize SUVpeak, most typically defined as the average SUV within a small volume (e.g., 1 cm3) within the highest uptake region of the tumor, or an SUV ratio. These have been suggested as a more robust alternative to SUVmax less sensitive to differences in noise, scanner resolution, and reconstruction, and both have some evidence to support their use in cervical cancer (45–47). Although it has limitations, SUVmax is simple to measure, is one of the most widely reported FDG-PET metrics in standard clinical interpretations, and we have previously reported its prognostic significance in a separate cohort of patients with cervical cancer (7, 48, 49). Regardless, our data have outlined high FDG uptake as a poor prognostic factor that is correlated with the inflammatory state and immune infiltrate in cervical cancer tumors, even if our specific SUVmax cutoff may not hold true across a more diverse array of scanners and institutions. Because of its broad clinical availability, the connection between FDG uptake and inflammatory and immune state could have immediate clinical and therapeutic implications.

Future work should focus on further study and characterization of the interplay between immune infiltrate and cervical cancer metabolism, and on validating our connection between high SUVmax for FDG and a macrophage-predominant immune infiltrate. With regard to the connection between glucose metabolism and immune infiltrate, here we have broadly established an interplay between the presence of macrophages, or macrophage-conditioned growth medium, and cervical cancer glucose metabolism. We focused on factors associated with JAK/STAT3 signaling, and IL6 in particular, as our GSEA and cytokine array data suggested this as a potential candidate for mediating this effect. However, the interactions between macrophages and cervical cancer cells are certainly complex, and mediated by a host of other factors. We have also recently demonstrated a connection between FDG uptake, oxidative stress, and radioresistance in cervical cancer cell lines and animal models (50). Further work should explore the interplay between glucose utilization, oxidative stress, and the inflammatory and immune microenvironment in primary cervical tumors.

In sum, we have demonstrated that a high SUVmax for FDG is prognostic of outcomes in cervical cancer, and also associated with proinflammatory signaling pathways and a monocytic immune infiltrate. Ultimately, SUVmax may potentially be used as a surrogate marker for this microenvironment to identify cancers that can be treated with therapies targeting proinflammatory pathways (e.g., JAK/STAT3 signaling) or therapies that aim to reverse the infiltration of tumor-associated macrophages and other cells that suppress adaptive immunity.

J.M. Floberg reports grants from the Radiological Society of North America during the conduct of the study. B.A. Siegel reports grants from the NIH, ASTRO Resident Research Seed Grant, and RSNA Resident Research Grant during the conduct of the study; grants and personal fees from the American College of Radiology, Imaginab, Inc., and Blue Earth Diagnostics; and personal fees from the American Medical Foundation for Peer Review and Education, Avid Radiopharmaceuticals, Inc., Capella Imaging, LLC, Curium Pharma, Radiological Society of North America, and Siemens Healthineers outside the submitted work. S. Markovina reports grants from the NCI during the conduct of the study and grants and nonfinancial support from GlaskoSmithKline outside the submitted work. J.K. Schwarz reports grants from NIH, ASTRO, and RSNA during the conduct of the study; grants from Goldman Sachs Philanthropy Cancer Research Fund, Siteman Investment Program, 2019 AACR–Bristol Myers Squibb Mid-career Female Investigator Grant, NIH; and nonfinancial support from Calithera Biosciences outside the submitted work. No disclosures were reported by the other authors.

J.M. Floberg: Conceptualization, resources, data curation, formal analysis, funding acquisition, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. J. Zhang: Conceptualization, data curation, formal analysis, validation, visualization, methodology, writing–original draft, writing–review and editing. N. Muhammad: Investigation, visualization, methodology, writing–original draft, writing–review and editing. T.A. DeWees: Formal analysis, writing–original draft, writing–review and editing. M. Inkman: Formal analysis, writing–review and editing. K. Chen: Formal analysis, investigation, writing–review and editing. A.J. Lin: Formal analysis, investigation, methodology, writing–review and editing. R. Rashmi: Investigation, methodology, writing–review and editing. K. Jayachandran: Data curation, formal analysis, investigation, methodology, writing–review and editing. B.T. Edelson: Methodology, writing–review and editing. B.A. Siegel: Formal analysis, investigation, writing–review and editing. F. Dehdashti: Formal analysis, investigation, writing–review and editing. P.W. Grigsby: Conceptualization, resources, data curation, formal analysis, supervision, writing–review and editing. S. Markovina: Conceptualization, resources, supervision, investigation, methodology, writing–review and editing. J.K. Schwarz: Conceptualization, resources, supervision, funding acquisition, methodology, project administration, writing–review and editing.

We would like to thank Marina Platik and the Anatomic and Molecular Pathology Core Laboratory at Washington University in St. Louis for their assistance in staining the TMA. We would also like to acknowledge Songyang Wang, MD, PhD, and Min Tan from Stephanie Markovina's laboratory for their assistance with IHC and for supplying supplies, reagents, and expertise for investigating the JAK/STAT3 pathway. This work was supported in part by NIH R01CA181745 and the AACR–Bristol Myers Squibb Mid-career Female Investigator Award to J.K. Schwarz; by NIH K08CA237822 to S. Markovina, by NCI K22CA237839 to J. Zhang, and by ASTRO Resident Research Seed Grant 531448 and the RSNA Resident Research Grant to J.M. Floberg.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Jemal
A
,
Bray
F
,
Center
MM
,
Ferlay
J
,
Ward
E
,
Forman
D
. 
Global cancer statistics
.
CA Cancer J Clin
2011
;
61
:
69
90
.
2.
Eifel
PJ
,
Winter
K
,
Morris
M
,
Levenback
C
,
Grigsby
PW
,
Cooper
J
, et al
Pelvic irradiation with concurrent chemotherapy versus pelvic and para-aortic irradiation for high-risk cervical cancer: an update of radiation therapy oncology group trial (RTOG) 90-01
.
J Clin Oncol
2004
;
22
:
872
80
.
3.
Rose
PG
,
Ali
S
,
Watkins
E
,
Thigpen
JT
,
Deppe
G
,
Clarke-Pearson
DL
, et al
Long-term follow-up of a randomized trial comparing concurrent single agent cisplatin, cisplatin-based combination chemotherapy, or hydroxyurea during pelvic irradiation for locally advanced cervical cancer: a Gynecologic Oncology Group Study
.
J Clin Oncol
2007
;
25
:
2804
10
.
4.
Lanciano
R
,
Calkins
A
,
Bundy
BN
,
Parham
G
,
Lucci
JA
,
Moore
DH
, et al
Randomized comparison of weekly cisplatin or protracted venous infusion of fluorouracil in combination with pelvic radiation in advanced cervix cancer: a gynecologic oncology group study
.
J Clin Oncol
2005
;
23
:
8289
95
.
5.
Tewari
KS
,
Sill
MW
,
Penson
RT
,
Huang
H
,
Ramondetta
LM
,
Landrum
LM
, et al
Bevacizumab for advanced cervical cancer: final overall survival and adverse event analysis of a randomised, controlled, open-label, phase 3 trial (Gynecologic Oncology Group 240)
.
Lancet
2017
;
390
:
1654
63
.
6.
Schmidt
AM
,
Imesch
P
,
Fink
D
,
Egger
H
. 
Indications and long-term clinical outcomes in 282 patients with pelvic exenteration for advanced or recurrent cervical cancer
.
Gynecol Oncol
2012
;
125
:
604
9
.
7.
Kidd
EA
,
Siegel
BA
,
Dehdashti
F
,
Grigsby
PW
. 
The standardized uptake value for F-18 fluorodeoxyglucose is a sensitive predictive biomarker for cervical cancer treatment response and survival
.
Cancer
2007
;
110
:
1738
44
.
8.
Schwarz
JK
,
Siegel
BA
,
Dehdashti
F
,
Grigsby
PW
. 
Association of posttherapy positron emission tomography with tumor response and survival in cervical carcinoma
.
JAMA
2007
;
298
:
2289
95
.
9.
Kidd
EA
,
Siegel
BA
,
Dehdashti
F
,
Grigsby
PW
. 
Pelvic lymph node F-18 fluorodeoxyglucose uptake as a prognostic biomarker in newly diagnosed patients with locally advanced cervical cancer
.
Cancer
2010
;
116
:
1469
75
.
10.
Ho
JC
,
Allen
PK
,
Bhosale
PR
,
Rauch
GM
,
Fuller
CD
,
Mohamed
AS
, et al
Diffusion-weighted magnetic resonance imaging as a predictor of outcome in cervical cancer after chemoradiation
.
Int J Radiat Oncol Biol Phys
2017
;
97
:
546
53
.
11.
Halle
C
,
Andersen
E
,
Lando
M
,
Aarnes
EK
,
Hasvold
G
,
Holden
M
, et al
Hypoxia-induced gene expression in chemoradioresistant cervical cancer revealed by dynamic contrast-enhanced MRI
.
Cancer Res
2012
;
72
:
5285
95
.
12.
Schwarz
JK
,
Payton
JE
,
Rashmi
R
,
Xiang
T
,
Jia
Y
,
Huettner
P
, et al
Pathway-specific analysis of gene expression data identifies the PI3K/Akt pathway as a novel therapeutic target in cervical cancer
.
Clin Cancer Res
2012
;
18
:
1464
71
.
13.
Dyk
P
,
Jiang
N
,
Sun
B
,
DeWees
TA
,
Fowler
KJ
,
Narra
V
, et al
Cervical gross tumor volume dose predicts local control using magnetic resonance imaging/diffusion-weighted imaging-guided high-dose-rate and positron emission tomography/computed tomography-guided intensity modulated radiation therapy
.
Int J Radiat Oncol Biol Phys
2014
;
90
:
794
801
.
14.
Network CGAR, Medicine AECo, Services AB, Hospital BC, Medicine BCo, Hope BRIoCo
et al 
Integrated genomic and molecular characterization of cervical cancer
.
Nature
2017
;
543
:
378
84
.
15.
Wright
JD
,
Dehdashti
F
,
Herzog
TJ
,
Mutch
DG
,
Huettner
PC
,
Rader
JS
, et al
Preoperative lymph node staging of early-stage cervical carcinoma by [18F]-fluoro-2-deoxy-D-glucose-positron emission tomography
.
Cancer
2005
;
104
:
2484
91
.
16.
Contal
C
,
O'Quigley
J
. 
An application of changepoint methods in studying the effect of age on survival in breast cancer
.
Comput Stat Data Anal
1999
;
30
:
253
70
.
17.
Dobin
A
,
Gingeras
TR
. 
Mapping RNA-seq reads with STAR
.
Curr Protoc Bioinformatics
2015
;
51
:
11.4.1
9
.
18.
Liao
Y
,
Smyth
GK
,
Shi
W
. 
featureCounts: an efficient general purpose program for assigning sequence reads to genomic features
.
Bioinformatics
2014
;
30
:
923
30
.
19.
Trapnell
C
,
Williams
BA
,
Pertea
G
,
Mortazavi
A
,
Kwan
G
,
van Baren
MJ
, et al
Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation
.
Nat Biotechnol
2010
;
28
:
511
5
.
20.
White
NM
,
Zhao
SG
,
Zhang
J
,
Rozycki
EB
,
Dang
HX
,
McFadden
SD
, et al
Multi-institutional analysis shows that low PCAT-14 expression associates with poor outcomes in prostate cancer
.
Eur Urol
2017
;
71
:
257
66
.
21.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
,
Mukherjee
S
,
Ebert
BL
,
Gillette
MA
, et al
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci U S A
2005
;
102
:
15545
50
.
22.
Newman
AM
,
Liu
CL
,
Green
MR
,
Gentles
AJ
,
Feng
W
,
Xu
Y
, et al
Robust enumeration of cell subsets from tissue expression profiles
.
Nat Methods
2015
;
12
:
453
7
.
23.
Markovina
S
,
Wang
S
,
Henke
LE
,
Luke
CJ
,
Pak
SC
,
DeWees
T
, et al
Serum squamous cell carcinoma antigen as an early indicator of response during therapy of cervical cancer
.
Br J Cancer
2018
;
118
:
72
8
.
24.
Ring
KL
,
Yemelyanova
AV
,
Soliman
PT
,
Frumovitz
MM
,
Jazaeri
AA
. 
Potential immunotherapy targets in recurrent cervical cancer
.
Gynecol Oncol
2017
;
145
:
462
8
.
25.
Kamper
P
,
Bendix
K
,
Hamilton-Dutoit
S
,
Honoré
B
,
Nyengaard
JR
,
d'Amore
F
. 
Tumor-infiltrating macrophages correlate with adverse prognosis and Epstein-Barr virus status in classical Hodgkin's lymphoma
.
Haematologica
2011
;
96
:
269
76
.
26.
Kawachi
A
,
Yoshida
H
,
Kitano
S
,
Ino
Y
,
Kato
T
,
Hiraoka
N
. 
Tumor-associated CD204
.
Cancer Sci
2018
;
109
:
863
70
.
27.
Institute
B
. 
Gene set
enrichment analysis: MSigDB Collections
[database on the Internet]
; 
2017
.
Accessed March 15, 2019. Available from
: http://software.broadinstitute.org/gsea/msigdb/collections.jsp.
28.
Sobti
RC
,
Singh
N
,
Hussain
S
,
Suri
V
,
Bharti
AC
,
Das
BC
. 
Overexpression of STAT3 in HPV-mediated cervical cancer in a north Indian population
.
Mol Cell Biochem
2009
;
330
:
193
9
.
29.
Takemoto
S
,
Ushijima
K
,
Kawano
K
,
Yamaguchi
T
,
Terada
A
,
Fujiyoshi
N
, et al
Expression of activated signal transducer and activator of transcription-3 predicts poor prognosis in cervical squamous-cell carcinoma
.
Br J Cancer
2009
;
101
:
967
72
.
30.
Tilborghs
S
,
Corthouts
J
,
Verhoeven
Y
,
Arias
D
,
Rolfo
C
,
Trinh
XB
, et al
The role of nuclear factor-kappa B signaling in human cervical cancer
.
Crit Rev Oncol Hematol
2017
;
120
:
141
50
.
31.
DA Costa
RM
,
Bastos
MM
,
Medeiros
R
,
Oliveira
PA
. 
The NFκB signaling pathway in papillomavirus-induced lesions: friend or foe?
Anticancer Res
2016
;
36
:
2073
83
.
32.
Luan
S
,
An
Z
,
Bi
S
,
Chen
L
,
Fan
J
. 
Interleukin 6 receptor (IL-6R) was an independent prognostic factor in cervical cancer
.
Histol Histopathol
2018
;
33
:
269
76
.
33.
Wisdom
AJ
,
Hong
CS
,
Lin
AJ
,
Xiang
Y
,
Cooper
DE
,
Zhang
J
, et al
Neutrophils promote tumor resistance to radiation therapy
.
Proc Natl Acad Sci U S A.
2019
;
116
:
18584
9
.
34.
Chan
WK
,
Au
WY
,
Wong
CY
,
Liang
R
,
Leung
AY
,
Kwong
YL
, et al
Metabolic activity measured by F-18 FDG PET in natural killer-cell lymphoma compared to aggressive B- and T-cell lymphomas
.
Clin Nucl Med
2010
;
35
:
571
5
.
35.
Chang
WC
,
Li
CH
,
Chu
LH
,
Huang
PS
,
Sheu
BC
,
Huang
SC
. 
Regulatory T cells suppress natural killer cell immunity in patients with human cervical carcinoma
.
Int J Gynecol Cancer
2016
;
26
:
156
62
.
36.
Li
M
,
Jin
R
,
Wang
W
,
Zhang
T
,
Sang
J
,
Li
N
, et al
STAT3 regulates glycolysis via targeting hexokinase 2 in hepatocellular carcinoma cells
.
Oncotarget
2017
;
8
:
24777
84
.
37.
Zheng
M
,
Cao
MX
,
Yu
XH
,
Li
L
,
Wang
K
,
Wang
SS
, et al
STAT3 promotes invasion and aerobic glycolysis of human oral squamous cell carcinoma via inhibiting FoxO1
.
Front Oncol
2019
;
9
:
1175
.
38.
Poli
V
,
Camporeale
A
. 
STAT3-mediated metabolic reprograming in cellular transformation and implications for drug resistance
.
Front Oncol
2015
;
5
:
121
.
39.
Piersma
SJ
. 
Immunosuppressive tumor microenvironment in cervical cancer patients
.
Cancer Microenviron
2011
;
4
:
361
75
.
40.
Pedraza-Brindis
EJ
,
Sánchez-Reyes
K
,
Hernández-Flores
G
,
Bravo-Cuellar
A
,
Jave-Suárez
LF
,
Aguilar-Lemarroy
A
, et al
Culture supernatants of cervical cancer cells induce an M2 phenotypic profile in THP-1 macrophages
.
Cell Immunol
2016
;
310
:
42
52
.
41.
Zijlmans
HJ
,
Fleuren
GJ
,
Baelde
HJ
,
Eilers
PH
,
Kenter
GG
,
Gorter
A
. 
The absence of CCL2 expression in cervical carcinoma is associated with increased survival and loss of heterozygosity at 17q11.2
.
J Pathol
2006
;
208
:
507
17
.
42.
Stover
DG
,
Gil Del Alcazar
CR
,
Brock
J
,
Guo
H
,
Overmoyer
B
,
Balko
J
, et al
Phase II study of ruxolitinib, a selective JAK1/2 inhibitor, in patients with metastatic triple-negative breast cancer
.
NPJ Breast Cancer
2018
;
4
:
10
.
43.
Hurwitz
HI
,
Uppal
N
,
Wagner
SA
,
Bendell
JC
,
Beck
JT
,
Wade
SM
, et al
Randomized, double-blind, phase II study of ruxolitinib or placebo in combination with capecitabine in patients with metastatic pancreatic cancer for whom therapy with gemcitabine has failed
.
J Clin Oncol
2015
;
33
:
4039
47
.
44.
Velasquez
LM
,
Boellaard
R
,
Kollia
G
,
Hayes
W
,
Hoekstra
OS
,
Lammertsma
AA
, et al
Repeatability of 18F-FDG PET in a multicenter phase I study of patients with advanced gastrointestinal malignancies
.
J Nucl Med
2009
;
50
:
1646
54
.
45.
Vanderhoek
M
,
Perlman
SB
,
Jeraj
R
. 
Impact of the definition of peak standardized uptake value on quantification of treatment response
.
J Nucl Med
2012
;
53
:
4
11
.
46.
Schernberg
A
,
Reuze
S
,
Orlhac
F
,
Buvat
I
,
Dercle
L
,
Sun
R
, et al
A score combining baseline neutrophilia and primary tumor SUV
.
Eur J Nucl Med Mol Imaging
2018
;
45
:
187
95
.
47.
Kunos
C
,
Radivoyevitch
T
,
Abdul-Karim
FW
,
Faulhaber
P
. 
18F-fluoro-2-deoxy-D-glucose positron emission tomography standard uptake value ratio as an indicator of cervical cancer chemoradiation therapeutic response
.
Int J Gynecol Cancer
2011
;
21
:
1117
23
.
48.
Lodge
MA
,
Chaudhry
MA
,
Wahl
RL
. 
Noise considerations for PET quantification using maximum and peak standardized uptake value
.
J Nucl Med
2012
;
53
:
1041
7
.
49.
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
:
122S
50S
.
50.
Rashmi
R
,
Huang
X
,
Floberg
JM
,
Elhammali
AE
,
McCormick
ML
,
Patti
GJ
, et al
Radioresistant cervical cancers are sensitive to inhibition of glycolysis and redox metabolism
.
Cancer Res
2018
;
78
:
1392
403
.