Abstract
Loss of TGFβ signaling increases error-prone alternative end-joining (alt-EJ) DNA repair. We previously translated this mechanistic relationship as TGFβ and alt-EJ gene expression signatures, which we showed are anticorrelated across cancer types. A score representing anticorrelation, βAlt, predicts patient outcome in response to genotoxic therapy. Here we sought to verify this biology in live specimens and additional datasets.
Human head and neck squamous carcinoma (HNSC) explants were treated in vitro to test whether the signatures report TGFβ signaling, indicated by SMAD2 phosphorylation, and unrepaired DNA damage, indicated by persistent 53BP1 foci after irradiation or olaparib. A custom NanoString assay was implemented to analyze the signatures’ expression in explants. Each signature gene was then weighted by its association with functional responses to define a modified score, βAltw, that was retested for association with response to genotoxic therapies in independent datasets.
Most genes in each signature were positively correlated with the expected biological response in tumor explants. Anticorrelation of TGFβ and alt-EJ signatures measured by NanoString was confirmed in explants. βAltw was significantly (P < 0.001) better than βAlt in predicting overall survival in response to genotoxic therapy in The Cancer Genome Atlas (TCGA) pancancer patients and in independent HNSC and ovarian cancer patient datasets.
Association of the TGFβ and alt-EJ signatures with their biological response validates TGFβ competency as a key mediator of DNA repair that can be readily assayed by gene expression. The predictive value of βAltw supports its development to assist in clinical decision making.
The anticorrelation of TGFβ and error-prone alternative end-joining transcriptomic signatures represents a mechanistic relationship and indicates key biological processes that provide clinical insight. The βAltw score, or a similar means to assess this biology, may serve as a predictive biomarker for patients receiving genotoxic therapy in either head and neck squamous carcinoma or ovarian cancer. The clinical utility of these signatures needs to be further validated in a prospective clinical trial to determine whether the βAltw score can provide sufficient predictive information to stratify and help guide patient management. If so, this mechanism-based score could enable more personalized cancer therapy to assist in clinical decision making.
Introduction
Of its myriad roles, TGFβ maintenance of genomic stability is among the least studied. TGFβ regulates the expression or function of key DNA repair proteins encoded by ATM (ataxia telangiectasia mutated; ref. 1), BRCA1 (breast cancer 1 gene; refs. 2, 3), FANCD2 (Fanconi anemia complementation group D2; ref. 4), and LIG4 (DNA ligase 4; ref. 5), which are necessary for maintenance of genomic integrity (reviewed in ref. 6). Inhibition or loss of TGFβ signaling decreases canonical DNA repair by homologous recombination (HRR) and non-homologous end-joining (NHEJ), which sensitizes cells to genotoxic treatments (1, 3, 5). Faulty DNA repair is a hallmark of cancer, and specific repair defects can provide the basis for response to precision therapies (7). Failure of HRR or NHEJ can increase the use of less-efficient repair by alternative end-joining (alt-EJ, also called microhomology-mediated end-joining; refs. 8–11). Because inhibition of TGFβ signaling increases sensitivity to DNA damage by radiation or chemotherapy in preclinical models of breast, brain, lung, and head and neck cancer (3, 12–15), identifying TGFβ signaling defects in cancer may present a specific therapeutic opportunity (16).
Patients with head and neck squamous cell carcinoma (HNSC) associated with human papilloma virus (HPV) etiology exhibit a striking sensitivity to standard genotoxic therapy with cisplatin and radiotherapy (17). In addition to degrading p53 and RB proteins, HPV also blocks TGFβ by targeting its receptors and signal transduction (18). We have previously demonstrated that loss of TGFβ signaling competency in HPV-positive cancer is the mechanism by which canonical DNA double-strand break (DSB) repair shifts from HRR to alt-EJ. This repair pathway choice is recapitulated by blocking TGFβ signaling in HPV-negative cells (3). Consistent with compromised DNA repair upon loss of functional TGFβ signaling, cases of HPV-negative HNSC from The Cancer Genome Atlas (TCGA) characterized by low expression of TGFβ target genes were also associated with better overall survival (OS) following standard-of-care chemoradiation compared with those with high expression of TGFβ target genes (3). Similar regulation of the DNA damage response by TGFβ was also evident in glioblastoma cell lines and glioblastoma specimens in TCGA (16).
We translated the mechanistic relationship between TGFβ and DNA repair into functional gene expression signatures consisting of TGFβ targets and genes necessary for alt-EJ (16). The TGFβ signature consists of 50 genes that were reciprocally regulated in MCF10A treated for 7 days with TGFβ or a small-molecule inhibitor of the TGFβ type I receptor kinase (3, 19). The 36-gene alt-EJ competency signature was curated from the literature and a screen of DNA damage repair gene using the EJ2GFP reporter (16, 20, 21). To classify patients according to the relationship between their TGFβ and alt-EJ transcriptional profiles, we defined a score, βAlt, to convey the relative expression of these signatures in each tumor. A high βAlt score, indicative of high alt-EJ and low TGFβ gene expression, correlates with specific mutational signatures, genomic instability characteristics, and better patient outcome in response to genotoxic treatment (16).
Given that both gene signatures were derived from in vitro studies of human cell lines, the nontransformed breast cell line MCF10A in the case of TGFβ target genes (19) and U2OS osteosarcoma cells for alt-EJ components (22), here we sought to validate that each signature reflects its respective biology in cancer specimens. We designed a custom NanoString panel to test each gene's functional association with biological response that could facilitate application to retrospective analysis of archival specimens. Through functional assays in HNSC tumor explants, targeted gene expression analyses, and integrative data modeling we demonstrate the biological coherence of the signatures, and subsequently defined an optimized score, termed βAltw, that was validated in independent datasets. The results confirm that TGFβ regulation of DNA damage repair is an important biomarker of outcome after genotoxic treatment with radiotherapy or platinum chemotherapy. The developed tools and methods can be implemented to prospectively predict patient response.
Materials and Methods
Explant cultures
HNSC tumor tissue and ovarian cancer specimen were collected from patients who consented in writing under the U.S. Common Rule and was reviewed and approved by the University of California, San Francisco (UCSF, San Francisco, CA) Institutional Review Board. Establishment of HNSC patient-derived xenograft (PDX) and tumor collection were reviewed and approved by the UCSF Institutional Animal Care and Use Committee. Explants were established on floating rafts from HNSC PDX from immunocompromised mice (n = 15) and patient primary tumors (n = 22) collected during surgery as described previously (3). Characteristics of patients donating primary tissues are shown in Supplementary Table S1. Specimens were kept in DMEM and transported on ice. Tissues were dissected and one portion was frozen in liquid nitrogen for RNA extraction, and the remainder was portioned into approximately 1 mm3 fragments for explant cultures. Some explants were irradiated with 5 Gy at 48 hours; others were treated at 24 hours with 10 μm olaparib for 24 hours. All were harvested 5 hours after treatment. Specimens were embedded in OCT (Sigma-Aldrich) and frozen on ethanol/dry ice bath. Blocks were kept at −80°C before cryosectioning.
Immunofluorescence and microscopy
Staining and visualization of 53BP1 DNA damage foci, indicative of unrepaired DNA damage, was assessed in treated explants and pSMAD2 was assayed in untreated explants as described previously (3). In brief, cryosections were fixed with 4% paraformaldehyde for 15 minutes and then permeabilized with 0.5% TritonX-100 followed by blocking with the supernatant of 0.5% casein stirred for 1 h in PBS. Sections were incubated with antibodies to 53BP1 (Bethyl, catalog no. A700-011, 1:500) or SMAD2 phosphorylated on serine 465/467 at (Cell Signaling Technology, catalog no. 3108, 1:200) at 4°C overnight in a humidified chamber. After three rinses with PBS, sections were incubated for 1 hour with secondary donkey anti-rabbit IgG (Alexa Fluor 488/555, Invitrogen) or donkey anti-mouse IgG (Alexa Fluor 488/555, Invitrogen). Cell nuclei were counterstained with 4′,6-diamidino-2-phenylindole, dihydrochloride (DAPI). Slides were mounted in Vectashield (Sigma). A 40× objective with 0.95 numerical aperture was used on a Zeiss Axiovert equipped with epifluorescence. In-home developed macros in the opensource platform Fiji-ImageJ (NIH, Bethesda, MA) were used for image analyses of 8-bit images for each channel of fluorescence. The DAPI channel was used to generate the region of interest and five or more images per sample were randomly taken based on nuclear dye alone. At least 100 cells were analyzed for each sample. For analyses of radiation induced foci, spontaneous foci from sham-treated controls were subtracted unless otherwise noted.
NanoString assay
A custom NanoString panel was used consisting of 50 genes induced by chronic TGFβ, 36 genes necessary for execution of alt-EJ (16), and 12 housekeeping genes. Total RNA was extracted from samples using the TRIzol reagent (Invitrogen) and the miRNAeasy Mini Kit (Qiagen) for frozen HNSC or ovarian cancer samples. Formalin-fixed paraffin-embedded ovarian cancer specimens were reviewed and approved by the Institutional Review Board at L’Hospitalet del Llobregat under the Declaration of Helsinki. RNA was prepared from formalin-fixed paraffin-embedded ovarian cancer sections using the RNeasy FFPE Kit (Qiagen) according to manufacturer's instructions. A total of 250 ng of total RNA from each sample were hybridized following the manufacturer's protocol (NanoString Technologies). Gene expression was quantified using the standard nCounter methodology with multiplexed color-coded probe pairs (23). The raw expression data were processed and normalized using the nSolver software (NanoString). Normalized counts for HNSC and ovarian specimens were separately log2 transformed and mean-centered per gene by converting into z-scores.
Bioinformatic analyses
Gene expression heatmaps were constructed with unsupervised hierarchical clustering using the R package ComplexHeatmap (24). Euclidean distance was used as the similarity metric and the Ward.D2 method as the between-cluster distance metric. The gene correlation matrix was created by computing the Pearson correlation coefficient (rp) between the expression of every pair of TGFβ and alt-EJ signature genes, using the R package corrplot (https://www.rdocumentation.org/packages/corrplot). Genes were displayed on the basis of the weight of their contribution to the first principal component of the gene expression profiles.
Weighted gene coexpression networks were built using the R package ggraph (https://ggraph.data-imaginist.com). Each gene was represented as a node and was colored by its signature. An edge between two genes corresponded to the expression correlation (rp); and rp > 0.007 was used as a cutoff for edge generation. Layout force-directed algorithm from Fruchterman-Reingold was applied for network construction. To identify central genes in the network, the weighted centrality degree of each gene was calculated considering the number and weight of the edges connecting to any other gene of the same signature.
The βAltw score
To compute βAltw, signature genes were weighted using biological data generated from the HNSC explants based on: (i) a centrality degree within the corresponding signature calculated on the basis of the gene coexpression network analysis after rescaling it into a 0–0.5 range; (ii) the strength of the positive correlation (Spearman correlation coefficient, rs) between expression and the corresponding biologic measurement (frequency of cells with pSMAD2-positive nuclei for the genes from the TGFβ signature or cells with five or more 53BP1 foci indicative of unrepaired DNA damage for the genes from the alt-EJ signature); and (iii) the strength of the negative correlation (rs multiplied by −1) between expression and the other biologic measurement (53BP1 for the genes from the TGFβ signature and pSMAD2 for the genes from the alt-EJ signature).
Next, a factor (Supplementary Fig. S4) was assigned to each gene based on its weight with the following formula:
If weight genei > 0: Factor genei = 1 + weight genei
If weight genei ≤ 0: Factor genei = 1
Datasets
Gene expression data of TCGA-pancancer cohort were downloaded from the Genomic Data Commons portal from the file EBPlusPlusAdjustPANCAN_IlluminaHiSeq_RNASeqV2.geneExp.tsv. The downloaded gene expression values were trimmed mean of M values normalized, log2 transformed and mean-centered per gene by converting them into z-scores. Primary solid tumor samples were analyzed. Glioblastoma samples categorized as “neural” were excluded (25), as were mislabeled pancreatic cancer samples (26). OS, tumor stage, and age from all TCGA patients whose standard of care would include genotoxic radiotherapy and/or chemotherapy based on their cancer type and stage (n = 4597) were obtained from the dataset “TCGA-CDR-SupplementalTableS1.xlsx” from Genomic Data Commons in November 2020 (16).
Gene expression of TCGA-HNSC cohort was downloaded from the dataset “TCGA.HNSC.sampleMap/HiSeqV2” of the University of California Santa Cruz Xena platform using the R package UCSCXenaTools in November 2020, excluding those patients whose primary curative treatment had been surgery so that the remaining (n = 419) were likely to have received genotoxic treatment with radiotherapy and/or chemotherapy. Gene expression had been measured by RNA sequencing (RNA-seq) with the platform IlluminaHiSeq_RNASeqV2 and values had been RSEM normalized and log2(x+1) transformed. The downloaded gene expression values were mean centered per gene by converting them into z-scores for primary tumor samples. OS, tumor stage, and age from TCGA-HNSC patients were obtained from the dataset “TCGA-CDR-SupplementalTableS1.xlsx” from Genomic Data Commons in November 2020. Information about HPV status from the TCGA-HNSC tumors was downloaded from cBioPortal in November 2020 from the project "Head and Neck Squamous Cell Carcinoma (TCGA, PanCancer Atlas)." TCGA-HNSC patient treatment information was downloaded from the University of California Santa Cruz Xena platform using the R package UCSCXenaTools (27).
The GSE41613 HNSC dataset (n = 97) and the GSE26712 ovarian cancer dataset (n = 185) were downloaded from the Gene Expression Omnibus in March and April 2021, respectively, using the R package GEOquery. In GSE41613, gene expression had been measured with the platform GPL570 and normalized into log2 gcRMA signal. In GSE26712, gene expression had been measured with the platform GPL96 and normalized into RMA signal value. For genes with multiple probes, the average expression of the probes was calculated. The downloaded gene expression values were mean-centered per gene by converting them into z-scores.
Gene expression data from the NCI-60 cell lines (n = 60) were downloaded from cBioPortal in April 2021 and the surviving fraction after 2 Gy (SF2) of the cell lines was obtained from the literature based on reported clonogenic assays (28).
Statistical analyses
Kaplan–Meier survival curves, as defined in TCGA (29), GEO GSE26712 (30), and GSE41613 (31), were generated for βAltw score tertile 1 versus 3 via R package “survminer” and the multivariable Cox regressions were performed using the R package “survival” as indicated. Adjusted HR with corresponding 95% confidence interval (CI) were reported. The Cox proportional hazards assumption was assessed graphically with the “coxph” function via the Schoenfeld residuals, and there was no evidence of PH violation. All statistical comparisons were two sided and considered as statistically significant at P < 0.05.
To compare the performance of the βAltw score and the original βAlt in terms of predicting OS after genotoxic treatment, we followed a similar methodology as described previously (32). TCGA-pancancer dataset of patients whose standard of care would include genotoxic radiotherapy and/or chemotherapy (n = 4,597) was randomly split into 500 surrogate datasets using the bootstrapping resampling method with the R package “boot” to calculate the original βAlt score and the βAltw score. The OS HR calculated according to each score for the top and bottom tertile. The resulting HR obtained with the βAltw score and the original βAlt were compared using a paired t test.
Code availability
Data availability
All data in the article that are not from open-access datasets are available from the corresponding author upon request.
Results
TGFβ and alt-EJ signatures associate with biological readouts of TGFβ signaling and DNA repair
Gene expression of TGFβ and alt-EJ signatures was assessed by a targeted NanoString custom panel in 15 HNSC PDXs and 22 HNSC patients’ primary tumors from which we then established explants (Fig. 1A). Unsupervised hierarchical clustering of the HNSC tumors and their transcriptomic phenotype clustered genes from both signatures into two major clades characterized by low TGFβ and high alt-EJ or high TGFβ and low alt-EJ gene expression (Fig. 1B), reproducing the anticorrelation we reported using RNA-seq data from TCGA (16). As reported therein, most HPV-positive samples were in the dendrogram arm with low TGFβ and high alt-EJ.
NanoString assays are well suited to retrospective studies (33). To test how the panel reports on the TGFβ and alt-EJ pathways in archival specimens, RNA was extracted from 18 snap-frozen and 22 formalin-fixed paraffin embedded tissue sections of ovarian cancer. Unsupervised hierarchical clustering of the NanoString signatures revealed two clades characterized by low TGFβ and high alt-EJ or high TGFβ and low alt-EJ genes expression (Supplementary Fig. S1A). Duplicates (indicated by red bars) were adjacent, demonstrating the reproducibility of gene expression measured by NanoString.
We next sought to validate the biological significance of the TGFβ signature by assessing pathway activity by measuring phosphorylated SMAD2 (pSMAD2; Fig. 2A). DNA repair proficiency was assessed 5 hours after irradiation with 5 Gy by quantifying persistent 53BP1 foci (Fig. 2B; ref. 34). The mean expression of the TGFβ signature genes was significantly correlated with the percentage of pSMAD2-positive cells (rs = 0.45, P = 0.0067), supporting that the signature is indicative of TGFβ signaling competency (Fig. 2C). Likewise, the mean expression of the alt-EJ signature was positively correlated with unrepaired DNA damage (rs = 0.51, P = 0.03), supporting the alt-EJ signature as one of less efficient DNA repair (Fig. 2D).
The percentage of pSMAD2- and 53BP1-positive cells were negatively correlated across all specimens (Fig. 2E; rs = −0.5, P = 0.037), confirming the correlation between robust TGFβ signaling and competent DNA damage repair (3). Consistent with our prior study, HPV-positive specimens had fewer pSMAD2-positive cells and more cells with 53BP1 foci indicative of unrepaired DNA damage and the TGFβ and alt-EJ signatures were negatively correlated (rs = −0.41, P = 0.013; Fig. 2F).
Olaparib, an FDA-approved PARP inhibitor, can cause DSB in cells that are HRR deficient. Our prior work indicated that TGFβ loss or inhibition decreases HRR and increases sensitivity to PARP inhibition (3). To further test this relationship, we treated explants with olaparib for 24 hours and examined 53BP1 foci after 5 hours of recovery. As for irradiated explants, unrepaired DNA damage marked by 53BP1 foci in olaparib-treated explants was positively correlated with alt-EJ signature expression (rs = 0.59, P = 0.01; Fig. 2G). 53BP1 foci were anticorrelated with mean expression of the TGFβ signature (rs = −0.61, P = 0.005) and the percent pSMAD2-positive cells (rs = −0.65, P = 0.004; Fig. 2H). Thus, tumors with low TGFβ signaling indicated by pSMAD2 expression or TGFβ signature expression have more unrepaired damage following PARP inhibition.
To evaluate the relevance of each TGFβ and alt-EJ gene, we next assessed the expression of each gene per sample as a function of the percentage of pSMAD2- and 53BP1-positive cells (Fig. 3A and B). Most TGFβ genes were positively correlated with the frequency of pSMAD2-positive cells (37/50) and negatively correlated with 53BP1-positive cells (36/50), whereas most alt-EJ genes were positively correlated with residual DNA damage marked by 53BP1 (30/36) and negatively correlated with pSMAD2 (31/36).
Next, the Pearson correlation coefficient (rp) was computed between expression values of each pair of genes to construct a correlation matrix (Fig. 3C). Endorsing the premise that each signature embodies a distinct functional pathway, most genes from the same signature were highly correlated. Conversely, gene pairs between the signatures were found to be commonly anticorrelated, consistent with TGFβ suppression of alt-EJ genes (16). We then constructed TGFβ and alt-EJ gene correlation networks in the HNSC specimens to identify “hub” genes (i.e., genes with a relatively high number of positive correlations, expected to be functionally more relevant within the signature) and calculated the weighted centrality degree of each component (Supplementary Fig. S2A and S2B). Similar networks were computed from the NanoString analysis of the ovarian cancer specimens (Supplementary Fig. S2C and S2D). Notably, functional ranks of the genes based on their association with pSMAD2 and 53BP1 cellular readouts were similar to their ranking by network centrality (Fig. 4).
Association between gene ranks and clinical outcomes
The previously defined βAlt score calculated using single-sample gene set enrichment analysis of the TGFβ and alt-EJ signatures reports their anticorrelation and is associated with response to genotoxic cancer therapy (16). A high βAlt score represented specimens in which the expression of TGFβ target genes was low and expression of alt-EJ genes was high, whereas a low βAlt score indicated the opposite. To test whether the above evidence of differences in the strength of association between biology and specific genes in each signature improved the translational relevance of the signatures, we calculated βAlt score using only the top-ranked 15 genes or the bottom-ranked 15 genes from each signature for patients with HNSC from TCGA, excluding those whose primary curative treatment had been surgery (n = 419). Survival was significantly associated with a βAlt score calculated using only the top 15 genes (P = 0.0093; HR = 0.62; 95% CI, 0.43–0.89; Supplementary Fig. S3A), whereas one calculated using the bottom 15 genes was not (P = 0.49; HR = 0.89; 95% CI, 0.63–1.24; Supplementary Fig. S3B). This suggests that the genes with the greatest weight based on the biology and centrality exhibited by the HNSC tumor explants are also the most clinically relevant.
Given that 11 TGFβ signature genes and 4 alt-EJ signature genes were negatively weighted (Fig. 4), we examined whether excluding these genes would increase βAlt prognostic power compared to the original score. Notably, both the original (P = 0.00099; HR = 0.53; 95% CI, 0.37–0.78; Supplementary Fig. S3C) and shortened βAlt scores (P = 0.0036; HR = 0.59; 95% CI, 0.41–0.84; Supplementary Fig. S3D) were significantly associated with HNSC patient OS (P = 0.0009; HR = 0.53; 95% CI, 0.37–0.78 versus P = 0.0036; HR = 0.59; 95% CI, 0.41–0.84). Hence, we concluded that removing low weighted genes does not improve βAlt prognostic capacity.
Development of the βAltw score and validation of its predictive value
We next developed a model, termed βAltw, based on the sum of gene expression levels weighted by their estimated functional relevance (Materials and Methods). To compare the performance of the βAltw score to the original βAlt score we used pancancer patients in TCGA whose standard of care would include radiotherapy and/or genotoxic chemotherapy based on their cancer type and stage (n = 4,597; ref. 16) and the bootstrap resampling method to split the dataset into 500 surrogate datasets (35). Patients in the top tertile had significantly longer OS compared with those in the bottom tertile using either the original βAlt score (Fig. 5A and B) or the βAltw score (Fig. 5C and D). However, the HR (0.61; 95% CI, 0.54–0.68) for the βAltw score of the resampled patient tertiles were significantly lower (P < 0.0001) than those calculated using the original βAlt score (HR = 0.64; 95% CI, 0.57–0.73), indicating superiority of the βAltw score as a biomarker of prognosis in response to genotoxic treatment.
The βAltw score is mechanistically based on TGFβ control of the response to DNA damage (16). To confirm this, we evaluated βAltw prediction of the radiosensitivity of NCI-60 pancancer cell lines (n = 60). As expected, the weighted TGFβ and alt-EJ signatures were significantly anticorrelated (rs = −0.51, P < 0.0001; Supplementary Fig. S5A). Radiation sensitivity, measured as the surviving fraction after exposure to 2 Gy (SF2), was significantly correlated with βAltw (rs = −0.36, P = 0.0046; Supplementary Fig. S5B), which supports the functional validity of βAltw.
To further test the predictive power of βAltw, we reanalyzed TCGA-HNSC patients (n = 419). As shown for the original βAlt (Supplementary Fig. S3A), patients with a high βAltw score had significantly better OS compared with those with a low βAltw (P = 0.019; HR = 0.65; 95% CI, 0.45–0.94; Fig. 6A). Although in a multivariable Cox regression adjusted for HPV status, age, and stage, the association of the βAltw score with OS compared as a function of tertile was not significant (P = 0.1099; HR = 0.697; 95% CI, 0.45–1.08; Supplementary Table S2), it was significant as a continuous variable (P = 0.0250; HR = 0.995; 95% CI, 0.99–0.99). To eliminate the potential impact of HPV or tumor location as confounding variables, we analyzed HPV-negative patients with oral squamous carcinoma from the GSE41613 dataset (n = 97). Patients from the high βAltw tertile had significantly better cancer-specific survival than those from the low βAltw tertile (P = 0.015; HR = 0.30; 95% CI, 0.11–0.84; Fig. 6B). Multivariable Cox regression analysis including age and stage maintained βAltw statistical significance (P = 0.0393; HR = 0.33; 95% CI, 0.11–0.95; Supplementary Table S2). These analyses support the predictive capacity of the βAltw score in patients with HNSC.
Prior analysis of βAlt in TCGA showed a significant association with outcomes of patients with ovarian cancer in which standard of care is genotoxic chemotherapy (16); thus we sought to test βAltw using ovarian cancer data from the GSE26712 dataset of patients with naïve stage II–III high-grade ovarian cancer treated with adjuvant platinum-based chemotherapy (n = 185). Patients from the high βAltw tertile experience significantly longer OS than those from the low βAltw tertile (P = 0.036; HR = 0.64; 95% CI, 0.42–0.97; Fig. 6C). In this dataset, lack of additional stage and age information precluded multivariable Cox regression analysis.
Patients with advanced stage ovarian cancer are first treated by extensive debulking surgery, followed by genotoxic platinum chemotherapy combined with paclitaxel. Patients who are optimally debulked have a substantially improved survival compared with patients who are left with bulky residual disease (36). Hence, we conducted a subset analysis of this dataset to determine whether the βAltw score is equally predictive of outcome as a function of debulking status. Comparison of the OS of the top and bottom βAltw tertile from optimally and suboptimally debulked patients showed that those with a high βAltw had substantially better outcomes (P = 0.0017; Fig. 6D) in response to standard-of-care treatment with platinum chemotherapy. Notably, suboptimally debulked patients with a high βAltw had an equivalent OS as patients who were optimally debulked, underscoring the prognostic significance of the βAltw score in patients treated with genotoxic therapies.
Discussion
The induction of DNA damage by radiotherapy or genotoxic chemotherapy is arguably the most widely deployed cancer treatment approach. However, chemoresistance/radioresistance constitutes a major obstacle to effective or personalized treatment. Thus, identifying novel predictive biomarkers is imperative to guide therapeutic decisions and improve cancer patients’ survival. In this study, we showed that TGFβ and alt-EJ reciprocal gene expression signatures, which are predictive of patient outcome in response to genotoxic therapy (16), are significantly correlated with their respective biology in HNSC explants. The TGFβ signaling signature was associated with the frequency of pSMAD2-positive cells whereas the alt-EJ signature was correlated with unrepaired DNA damage. These data formed the basis for weighting the signatures into a single score, termed βAltw, which was clinically validated in independent HNSC and ovarian cancer patient datasets.
Several insights were gained from this research. First, the strong correlation between biological response in the human HNSC explants treated in vitro together with the NanoString analysis of TGFβ and alt-EJ genes demonstrate that the signatures accurately report TGFβ signaling and functional responses to DNA damage, which further supports the extensive control of TGFβ signaling in DNA repair choice (16).
Second, the original βAlt biomarker used whole-transcriptome profiling, where the requirements for high-quality RNA or sufficient bulk tumor can impede widespread clinical adoption (37). Implementation of the custom NanoString panel is a cost-effective alternative to whole-genome expression profiling and suited for analysis of archival tissues, allowing assessment of large patient cohorts for which only formalin-fixed, paraffin-embedded tissue is available. The technical advantages of this method were substantiated by the reproducibility of the gene expression in archival ovarian cancer replicates.
Third, comparing the performance of predicting pancancer patient survival after genotoxic treatment between the original βAlt and the βAltw demonstrates the superiority of βAltw and that it is a robust predictor of patient prognosis in response to genotoxic treatment, as was confirmed in independent HNSC and ovarian cancer datasets. In HNSC, appropriate therapy decision and stratification remain a major challenge (38). HPV positivity in patients with oropharyngeal carcinoma is a potent stratification factor that is associated with better survival due to increased sensitivity to radiation and platinum chemotherapy (17, 39). HPV-positive cancer lack TGFβ responsiveness, which impairs DNA repair by HRR and increases the use of error-prone alt-EJ (3). Inhibiting TGFβ in HPV-negative HNSC also increases alt-EJ (3). However, HPV is not the only means by which cancers become TGFβ incompetent, which was substantiated by TCGA pancancer anticorrelation of these signatures (16). Here, we show that βAltw predictive capacity is maintained in HPV-negative patients with HNSC independently of clinical characteristics including age and stage. Efforts are increasing to deintensify standard cancer treatments in HPV-positive patients (40); likewise, identification of HPV-negative patients resistant to genotoxic therapy could be used to enrich selection for clinical trials of intensified therapy. Thus, if further validated in prospective clinical trials, this biomarker could provide predictive and prognostic information to personalize therapies to the individual patient's likelihood of response.
Fourth, the validation of the βAltw score in patients with ovarian cancer treated with genotoxic therapy highlights fundamental tumor biology that has direct implications for prognosis. Ovarian cancer is the most lethal gynecologic malignancy in women worldwide (41). Ovarian cancer is sensitive to platinum-based chemotherapy and the current standard is carboplatin and paclitaxel in the first-line setting. Interestingly, βAltw predicts that loss of TGFβ signaling overcomes the survival disadvantage shown for suboptimally debulked patients (30).
Moreover, the biology of alt-EJ points to additional therapeutic choices. Maintenance with PARP inhibitors after platinum-based therapy is approved for patients with ovarian cancer with BRCA1 or BRCA2 mutations but also provides significant benefit in patients with platinum-sensitive, relapsed, high-grade serous ovarian cancer (42). Repair by alt-EJ relies on polymerase θ (encoded by POLQ), for which clinically viable inhibitors were recently identified (43). Cells in which TGFβ signaling is blocked are sensitized to PARP inhibition and silencing POLQ increases response (3). Thus, targeting cancers using alt-EJ by inhibiting PARP or polymerase θ selectively kills cancer cells via synthetic lethality and spares TGFβ competent normal cells. Importantly, although TGFβ regulates HRR utilization via regulation of BRCA1, HRR deficiency is not required for alt-EJ to increase upon TGFβ inhibition (16). Hence, the anticorrelation of TGFβ and alt-EJ reported by βAlt is a clinically actionable basis by which to stratify patients.
This study has limitations. Although we validated the βAltw score in published datasets, this biomarker has not been tested prospectively in a clinical trial, which is a necessary step before adoption of the score as a prognostic tool in the clinical setting. Also, the description of genotoxic treatments given to patients in some public datasets, such as TCGA, is not annotated in individual detail, requiring inferences about patient treatments based on the standard of care of each cancer type and stage. In addition, gene expression levels do not necessarily reflect gene function, given that many proteins are not regulated at the transcriptional level. Finally, comparison of outcomes as a function of βAltw score tertile was used as an unbiased approach, but further studies are needed to determine the optimal cutoff for clinical application.
In summary, the TGFβ and alt-EJ transcriptomic signatures represent functional biological processes, and their anticorrelation provides important clinical insight. The βAltw score, or a similar means to assess this biology, may serve as a predictive biomarker for patients receiving genotoxic radiation or chemotherapy. The clinical utility of these signatures needs to be validated in a prospective clinical trial to determine whether the βAlt score can provide sufficient predictive information to stratify and help guide patient management in either HNSC or ovarian cancer. If so, this mechanism-based score could assist in clinical decision making and enable more personalized cancer therapy for patients.
Nonetheless, the biological validation of TGFβ and alt-EJ signatures, introduction of a custom gene set in a platform that can be used for retrospective analysis of existing specimens, and the development of the βAltw score support the further investigation of this biology to understand patient response to genotoxic therapies.
Authors' Disclosures
Q. Liu reports grants from the National Natural Science Foundation of China (grant no. 82073007) during the conduct of the study; in addition, Q. Liu has a patent entitled “DNA Damage Deficits in Cancer Cells” pending. M.A. Pujana reports a patent entitled “DNA Damage Deficits in Cancer Cells” pending. J. Piulats reports other support from Clovis and AstraZeneca and grants and other support from MSD and Pfizer during the conduct of the study, as well as grants and other support from BeiGene, BMS, Janssen, and Sanofi and other support from Astellas outside the submitted work. S.S. Yom reports grants from EMD Serono and Bristol-Myers Squibb outside the submitted work. A. Ashworth reports personal fees, non-financial support, and other support from Tango Therapeutics, Kyttaro, and Yingli; non-financial support and other support from Azkarra Therapeutics, Ovibio, Phoenix MD, Ambagon, Earli, Bluestar, Prolynx, Gladiator, Cura, and Circle; personal fees from Genentech, GSK, and Genvivo; grants and other support from SPARC; and grants from AstraZeneca outside the submitted work; and A. Ashworth holds patents on the use of PARP inhibitors held jointly with AstraZeneca from which he has benefitted financially (and may do so in the future). M.H. Barcellos-Hoff reports grants from NIH during the conduct of the study; personal fees and non-financial support from Pathway Innovations, Inc; grants, personal fees, and non-financial support from Genentech, Inc.; personal fees from EMD-Serano; and grants and personal fees from Varian outside the submitted work; in addition, M.H. Barcellos-Hoff has a patent entitled “DNA damage Deficit in Cancer Cells” pending. No disclosures were reported by the other authors.
Authors' Contributions
I. Guix: Data curation, software, formal analysis, writing–original draft, writing–review and editing. Q. Liu: Data curation, formal analysis, writing–original draft, writing–review and editing. M.A. Pujana: Data curation, supervision, writing–review and editing. P. Ha: Resources, writing–review and editing. J. Piulats: Conceptualization, writing–review and editing. I. Linares: Writing–review and editing. F. Guedea: Writing–review and editing. J.-H. Mao: Formal analysis, supervision, writing–review and editing. A. Lazar: Data curation, formal analysis, supervision, writing–original draft, writing–review and editing. J. Chapman: Resources, writing–review and editing. S.S. Yom: Writing–review and editing. A. Ashworth: Funding acquisition, methodology, writing–review and editing. M.H. Barcellos-Hoff: Conceptualization, resources, writing–original draft, writing–review and editing.
Acknowledgments
The authors would like to thank William Chou, Trevor Jones, and Colin Foster at UCSF for technical support. The results presented here are partly based on data generated by TCGA Research Network (https://www.cancer.gov/tcga), and we would like to express our gratitude to TCGA consortia and their coordinators for the data provision and clinical information used in this study.
This research was supported by funding from the NIH (R01 CA239235) and the UCSF Resource Allocation Program to M.H. Barcellos-Hoff and A. Ashworth. Tissue processing and immunostaining was performed using the Helen Diller Family Comprehensive Cancer Center Pathology Shared Resource, supported by the NCI of the NIH under Award Number P30CA082103. The National Natural Science Foundation of China (grant no. 82073007) to Q. Liu, the Generalitat de Catalunya (SGR 2017-449, CERCA program to IDIBELL) and Carlos III Institute of Health (ISCIII), funded by FEDER funds (PI21/01306) to M.A. Pujana.
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.