Purpose:

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.

Experimental Design:

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.

Results:

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.

Conclusions:

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.

Translational Relevance

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.

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.

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).

The weight of each TGFβ gene (Supplementary Fig. S4) was calculated as follows:
Likewise, the weight of each alt-EJ gene (Supplementary Fig. S4) was calculated as:

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

Taking this into account, the TGFβ and alt-EJ weighted expression scores were calculated for each tumor as the sum of the expression of the genes from each signature multiplied by their factors:
The βAltw score conveys in one value the relative expression of both signatures in each tumor and is computed as follows:

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.

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.

Figure 1.

Evaluation of TGFβ and alt-EJ gene expression signatures in HNSC tumors. A, Schematic illustration of 15 PDX and 22 primary HNSC tumors collected for NanoString or immunostaining assays. B, Unsupervised hierarchical clustering based on the expression (high, yellow, low, blue) of TGFβ (pink) and alt-EJ (green) signature genes using the NanoString custom platform and RNA extracted from the 37 HNSC samples. HPV positive (purple) or negative or not determined (gray) are indicated with the tissue origin from PDX (red) or primary patient specimens (light blue) in the bottom bar.

Figure 1.

Evaluation of TGFβ and alt-EJ gene expression signatures in HNSC tumors. A, Schematic illustration of 15 PDX and 22 primary HNSC tumors collected for NanoString or immunostaining assays. B, Unsupervised hierarchical clustering based on the expression (high, yellow, low, blue) of TGFβ (pink) and alt-EJ (green) signature genes using the NanoString custom platform and RNA extracted from the 37 HNSC samples. HPV positive (purple) or negative or not determined (gray) are indicated with the tissue origin from PDX (red) or primary patient specimens (light blue) in the bottom bar.

Close modal

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).

Figure 2.

TGFβ signaling status associates with DNA repair proficiency. A, Percentage of pSMAD2-positive cells in HNSC explants (n = 37). Representative images of a pSMAD2-low sample (*) and a pSMAD2-high sample (+). B, Percentage of HNSC cells with five or more 53BP1 foci at 5 hours after 5 Gy irradiation (n = 19). Representative images of a HNSC sample with few cells with residual 53BP1 foci (*) and a sample with high residual 53BP1 foci (+). C, Correlation of percent nuclear pSMAD2-positive cells with the expression of TGFβ signature genes calculated from NanoString data in the same sample. D, Percentage of 53BP1 foci positive cells after irradiation correlated with the expression of the alt-EJ signature genes from the same HNSC sample. E, Anticorrelation of percent nuclear pSMAD2-positive cells with percent 53BP1-positive cells after irradiation in the same panel of HNSC explants. F, Anticorrelation of TGFβ signature with alt-EJ signature genes’ expression in the same panel of HNSC explants. G, Percentage of 53BP1 foci positive cells are positively correlated with alt-EJ score from olaparib-treated HNSC explants (n = 19). H, Anticorrelation of TGFβ signature score with percentage of 53BP1 foci positive cells from olaparib-treated HNSC explants. Spearman correlations were performed for r and P values. Trend lines (blue) from linear regression are shown. Purple bars or dots are HPV-positive samples.

Figure 2.

TGFβ signaling status associates with DNA repair proficiency. A, Percentage of pSMAD2-positive cells in HNSC explants (n = 37). Representative images of a pSMAD2-low sample (*) and a pSMAD2-high sample (+). B, Percentage of HNSC cells with five or more 53BP1 foci at 5 hours after 5 Gy irradiation (n = 19). Representative images of a HNSC sample with few cells with residual 53BP1 foci (*) and a sample with high residual 53BP1 foci (+). C, Correlation of percent nuclear pSMAD2-positive cells with the expression of TGFβ signature genes calculated from NanoString data in the same sample. D, Percentage of 53BP1 foci positive cells after irradiation correlated with the expression of the alt-EJ signature genes from the same HNSC sample. E, Anticorrelation of percent nuclear pSMAD2-positive cells with percent 53BP1-positive cells after irradiation in the same panel of HNSC explants. F, Anticorrelation of TGFβ signature with alt-EJ signature genes’ expression in the same panel of HNSC explants. G, Percentage of 53BP1 foci positive cells are positively correlated with alt-EJ score from olaparib-treated HNSC explants (n = 19). H, Anticorrelation of TGFβ signature score with percentage of 53BP1 foci positive cells from olaparib-treated HNSC explants. Spearman correlations were performed for r and P values. Trend lines (blue) from linear regression are shown. Purple bars or dots are HPV-positive samples.

Close modal

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).

Figure 3.

Association of TGFβ and alt-EJ genes with functional biological readouts and gene-to-gene correlations. Volcano plots showing the statistical significance (y axis, −log10P value) of the correlation (x axis, rs) between the expression of the individual TGFβ (pink) and alt-EJ (green) signature genes and the percentage of pSMAD2-positive cells (A) or the percentage of cells with 53BP1 foci (B) in the panel of HNSC tumor explants. C, Gene correlation matrix showing the rp between the expression of each pair of TGFβ (pink) and alt-EJ (green) signature genes in the HNSC tumor explants (yellow = negative rp; blue = positive rp). Genes are displayed in first principal component order.

Figure 3.

Association of TGFβ and alt-EJ genes with functional biological readouts and gene-to-gene correlations. Volcano plots showing the statistical significance (y axis, −log10P value) of the correlation (x axis, rs) between the expression of the individual TGFβ (pink) and alt-EJ (green) signature genes and the percentage of pSMAD2-positive cells (A) or the percentage of cells with 53BP1 foci (B) in the panel of HNSC tumor explants. C, Gene correlation matrix showing the rp between the expression of each pair of TGFβ (pink) and alt-EJ (green) signature genes in the HNSC tumor explants (yellow = negative rp; blue = positive rp). Genes are displayed in first principal component order.

Close modal

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).

Figure 4.

Refinement of the TGFβ and alt-EJ signatures by weighting the functional relevance of each gene. Relative weight of the TGFβ signature genes (A) and the alt-EJ signature genes (B) according to the strength of their association with the frequency of pSMAD2-positive cells (column 1, rs for TGFβ genes and −rs for alt-EJ genes), the strength of their association with the frequency of 53BP1-positive cells (column 2, −rs for TGFβ genes and rs for alt-EJ genes), and their centrality degree (column 3) in the HNSC tumor explants. The fourth column represents the mean of the other three columns. The dot size indicates the absolute value, and the color indicates its direction (blue, positive; red, negative).

Figure 4.

Refinement of the TGFβ and alt-EJ signatures by weighting the functional relevance of each gene. Relative weight of the TGFβ signature genes (A) and the alt-EJ signature genes (B) according to the strength of their association with the frequency of pSMAD2-positive cells (column 1, rs for TGFβ genes and −rs for alt-EJ genes), the strength of their association with the frequency of 53BP1-positive cells (column 2, −rs for TGFβ genes and rs for alt-EJ genes), and their centrality degree (column 3) in the HNSC tumor explants. The fourth column represents the mean of the other three columns. The dot size indicates the absolute value, and the color indicates its direction (blue, positive; red, negative).

Close modal

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.

Figure 5.

Comparison of βAlt and βAltw as predictors of pancancer OS after genotoxic treatment. Cox regression analyses of 500 bootstrapping sets generated from TCGA-pancancer dataset. A, HR and 95% CI of the top versus bottom βAlt tertile of patients calculated for each sample set using a Cox model. B, Kaplan–Meier OS curves of the top (blue) versus bottom (red) βAlt tertile in the TCGA-pancancer dataset. P values were calculated with log-rank tests here and in C. C, HR and 95% CI of the top versus bottom βAltw tertile of patients calculated using a Cox model. D, Kaplan–Meier OS curves of the top (blue) versus bottom (red) βAltw tertile in TCGA-pancancer dataset. E, Comparison of HR for βAlt and βAltw for the 500 test sets (paired t test, ****, P < 0. 0001).

Figure 5.

Comparison of βAlt and βAltw as predictors of pancancer OS after genotoxic treatment. Cox regression analyses of 500 bootstrapping sets generated from TCGA-pancancer dataset. A, HR and 95% CI of the top versus bottom βAlt tertile of patients calculated for each sample set using a Cox model. B, Kaplan–Meier OS curves of the top (blue) versus bottom (red) βAlt tertile in the TCGA-pancancer dataset. P values were calculated with log-rank tests here and in C. C, HR and 95% CI of the top versus bottom βAltw tertile of patients calculated using a Cox model. D, Kaplan–Meier OS curves of the top (blue) versus bottom (red) βAltw tertile in TCGA-pancancer dataset. E, Comparison of HR for βAlt and βAltw for the 500 test sets (paired t test, ****, P < 0. 0001).

Close modal

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.

Figure 6.

The βAltw score predicts clinical outcomes after genotoxic therapy in independent HNSC and ovarian cancer datasets. A, βAltw top (blue) tertile is associated with better OS (log-rank, P = 0.019) compared with bottom (red) tertile of patients with HNSC from TCGA. B, βAltw top (blue) tertile is associated with better cancer-specific survival (log-rank, P = 0.015) compared with bottom (red) tertile of HPV-negative patients with oral squamous carcinoma (GSE41613). C, βAltw top (blue) tertile is associated with better OS (log-rank, P = 0.036) compared with bottom (red) tertile of naïve patients with stage II–III high-grade ovarian carcinoma treated with adjuvant platinum chemotherapy (GSE26712). D, The same population as in C classified by optimal (dark blue and dark red) versus suboptimal (light blue and orange) debulking status. βAltw top (light blue) tertile of suboptimally debulked patients is associated with better OS (log-rank, P = 0. 0017) compared with patients in the bottom (orange) tertile. In all plots, P values were calculated with log-rank tests.

Figure 6.

The βAltw score predicts clinical outcomes after genotoxic therapy in independent HNSC and ovarian cancer datasets. A, βAltw top (blue) tertile is associated with better OS (log-rank, P = 0.019) compared with bottom (red) tertile of patients with HNSC from TCGA. B, βAltw top (blue) tertile is associated with better cancer-specific survival (log-rank, P = 0.015) compared with bottom (red) tertile of HPV-negative patients with oral squamous carcinoma (GSE41613). C, βAltw top (blue) tertile is associated with better OS (log-rank, P = 0.036) compared with bottom (red) tertile of naïve patients with stage II–III high-grade ovarian carcinoma treated with adjuvant platinum chemotherapy (GSE26712). D, The same population as in C classified by optimal (dark blue and dark red) versus suboptimal (light blue and orange) debulking status. βAltw top (light blue) tertile of suboptimally debulked patients is associated with better OS (log-rank, P = 0. 0017) compared with patients in the bottom (orange) tertile. In all plots, P values were calculated with log-rank tests.

Close modal

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.

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.

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.

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.

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.

1.
Kirshner
J
,
Jobling
MF
,
Pajares
MJ
,
Ravani
SA
,
Glick
AB
,
Lavin
MJ
, et al
.
Inhibition of transforming growth factor-beta1 signaling attenuates ataxia telangiectasia mutated activity in response to genotoxic stress
.
Cancer Res
2006
;
66
:
10861
9
.
2.
Martinez-Ruiz
H
,
Illa-Bochaca
I
,
Omene
C
,
Hanniford
D
,
Liu
Q
,
Hernando
E
, et al
.
A TGFbeta-miR-182-BRCA1 axis controls the mammary differentiation hierarchy
.
Sci Signal
2016
;
9
:
ra118
.
3.
Liu
Q
,
Ma
L
,
Jones
T
,
Palomero
L
,
Pujana
MA
,
Martinez-Ruiz
H
, et al
.
Subjugation of TGFb Signaling by human papilloma virus in head and neck squamous cell carcinoma shifts DNA repair from homologous recombination to alternative end joining
.
Clin Cancer Res
2018
;
24
:
6001
14
.
4.
Chen
J
,
Shukla
V
,
Farci
P
,
Andricovich
J
,
Jogunoori
W
,
Kwong
LN
, et al
.
Loss of the transforming growth factor-β effector β2-Spectrin promotes genomic instability
.
Hepatology
2017
;
65
:
678
93
.
5.
Kim
MR
,
Lee
J
,
An
YS
,
Jin
YB
,
Park
IC
,
Chung
E
, et al
.
TGFbeta1 protects cells from gamma-IR by enhancing the activity of the NHEJ repair pathway
.
Mol C Res
2015
;
13
:
319
29
.
6.
Liu
Q
,
Lopez
K
,
Murnane
J
,
Humphrey
T
,
Barcellos-Hoff
MH
.
Misrepair in context: TGFbeta regulation of DNA repair
.
Front Oncol
2019
;
9
:
799
.
7.
Ceccaldi
R
,
Rondinelli
B
,
D'Andrea
AD
.
Repair pathway choices and consequences at the double-strand break
.
Trends Cell Biol
2016
;
26
:
52
64
.
8.
Iliakis
G
,
Murmann
T
,
Soni
A
.
Alternative end-joining repair pathways are the ultimate backup for abrogated classical non-homologous end-joining and homologous recombination repair: implications for the formation of chromosome translocations
.
Mutat Res Genet Toxicol Environ Mutagen
2015
;
793
:
166
75
.
9.
Sallmyr
A
,
Tomkinson
AE
.
Repair of DNA double-strand breaks by mammalian alternative end-joining pathways
.
J Biol Chem
2018
;
293
:
10536
46
.
10.
Sfeir
A
,
Symington
LS
.
Microhomology-mediated end joining: a back-up survival mechanism or dedicated pathway?
Trends Biochem Sci
2015
;
40
:
701
14
.
11.
Wood
RD
,
Doublie
S
.
DNA polymerase theta (POLQ), double-strand break repair, and cancer
.
DNA Repair
2016
;
44
:
22
32
.
12.
Bouquet
SF
,
Pal
A
,
Pilones
KA
,
Demaria
S
,
Hann
B
,
Akhurst
RJ
, et al
.
Transforming growth factor b1 inhibition increases the radiosensitivity of breast cancer cells in vitro and promotes tumor control by radiation in vivo
.
Clin Cancer Res
2011
;
17
:
6754
65
.
13.
Hardee
ME
,
Marciscano
AE
,
Medina-Ramirez
CM
,
Zagzag
D
,
Narayana
A
,
Lonning
SM
, et al
.
Resistance of glioblastoma-initiating cells to radiation mediated by the tumor microenvironment can be abolished by inhibiting transforming growth factor-beta
.
Cancer Res
2012
;
72
:
4119
29
.
14.
Zhang
M
,
Herion
TW
,
Timke
C
,
Han
N
,
Hauser
K
,
Weber
KJ
, et al
.
Trimodal glioblastoma treatment consisting of concurrent radiotherapy, temozolomide, and the novel TGF-β receptor I kinase inhibitor LY2109761
.
Neoplasia
2011
;
13
:
537
49
.
15.
Du
S
,
Bouquet
F
,
Lo
C-H
,
Pellicciotta
I
,
Bolourchi
S
,
Parry
R
, et al
.
Attenuation of the DNA damage response by TGFβ inhibitors enhances radiation sensitivity of NSCLC cells in vitro and in vivo
.
Int J Radiat Oncol Biol Phys
2014
;
91
:
91
9
.
16.
Liu
Q
,
Palomero
L
,
Moore
J
,
Guix
I
,
Espín
R
,
Aytés
A
, et al
.
Loss of TGFβ signaling increases alternative end-joining DNA repair that sensitizes to genotoxic therapies across cancer types
.
Sci Transl Med
2021
;
13
:
eabc4465
.
17.
Ang
KK
,
Harris
J
,
Wheeler
R
,
Weber
R
,
Rosenthal
DI
,
Nguyen-Tan
PF
, et al
.
Human papillomavirus and survival of patients with oropharyngeal cancer
.
N Engl J Med
2010
;
363
:
24
35
.
18.
Nees
M
,
Geoghegan
JM
,
Munson
P
,
Prabhu
V
,
Liu
Y
,
Androphy
E
, et al
.
Human papillomavirus type 16 E6 and E7 proteins inhibit differentiation-dependent expression of transforming growth factor-beta 2 in cervical keratinocytes
.
Cancer Res
2000
;
60
:
4289
98
.
19.
Andarawewa
KL
,
Erickson
AC
,
Chou
WS
,
Costes
SV
,
Gascard
P
,
Mott
JD
, et al
.
Ionizing radiation predisposes nonmalignant human mammary epithelial cells to undergo transforming growth factor b induced epithelial to mesenchymal transition
.
Cancer Res
2007
;
67
:
8662
70
.
20.
Bennardo
N
,
Cheng
A
,
Huang
N
,
Stark
JM
.
Alternative-NHEJ is a mechanistically distinct pathway of mammalian chromosome break repair
.
PLoS Genet
2008
;
4
:
e1000110
.
21.
DiBiase
SJ
,
Zeng
ZC
,
Chen
R
,
Hyslop
T
,
Curran
WJ
Jr
,
Iliakis
G
.
DNA-dependent protein kinase stimulates an independently active, nonhomologous, end-joining apparatus
.
Cancer Res
2000
;
60
:
1245
53
.
22.
Howard
SM
,
Yanez
DA
,
Stark
JM
.
DNA damage response factors from diverse pathways, including DNA crosslink repair, mediate alternative end joining
.
PLoS Genet
2015
;
11
:
e1004943
.
23.
Geiss
GK
,
Bumgarner
RE
,
Birditt
B
,
Dahl
T
,
Dowidar
N
,
Dunaway
DL
, et al
.
Direct multiplexed measurement of gene expression with color-coded probe pairs
.
Nat Biotechnol
2008
;
26
:
317
25
.
24.
Gu
Z
,
Eils
R
,
Schlesner
M
.
Complex heatmaps reveal patterns and correlations in multidimensional genomic data
.
Bioinformatics
2016
;
32
:
2847
9
.
25.
Sidaway
P
.
Glioblastoma subtypes revisited
.
Nat Rev Clin Oncol
2017
;
14
:
587
-.
26.
Cancer Genome Atlas Research Network
.
Integrated genomic characterization of pancreatic ductal adenocarcinoma
.
Cancer Cell
2017
;
32
:
185
203
.
27.
Wang
S
,
Liu
X
.
The UCSCXenaTools R package: a toolkit for accessing genomics data from UCSC Xena platform, from cancer multi-omics to single-cell RNA-seq
.
J Open Source Softw
2019
;
4
:
1627
.
28.
Eschrich
S
,
Zhang
H
,
Zhao
H
,
Boulware
D
,
Lee
JH
,
Bloom
G
, et al
.
Systems biology modeling of the radiation sensitivity network: a biomarker discovery platform
.
Int J Radiat Oncol Biol Phys
2009
;
75
:
497
505
.
29.
Liu
J
,
Lichtenberg
T
,
Hoadley
KA
,
Poisson
LM
,
Lazar
AJ
,
Cherniack
AD
, et al
.
An integrated TCGA Pan-Cancer clinical data resource to drive high-quality survival outcome analytics
.
Cell
2018
;
173
:
400
16
.
30.
Bonome
T
,
Levine
DA
,
Shih
J
,
Randonovich
M
,
Pise-Masison
CA
,
Bogomolniy
F
, et al
.
A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer
.
Cancer Res
2008
;
68
:
5478
86
.
31.
Lohavanichbutr
P
,
Méndez
E
,
Holsinger
FC
,
Rue
TC
,
Zhang
Y
,
Houck
J
, et al
.
A 13-gene signature prognostic of HPV-Negative OSCC: discovery and external validation
.
Clin Cancer Res
2013
;
19
:
1197
203
.
32.
Chen
EG
,
Wang
P
,
Lou
H
,
Wang
Y
,
Yan
H
,
Bi
L
, et al
.
A robust gene expression-based prognostic risk score predicts overall survival of lung adenocarcinoma patients
.
Oncotarget
2018
;
9
:
6862
71
.
33.
Veldman-Jones
MH
,
Brant
R
,
Rooney
C
,
Geh
C
,
Emery
H
,
Harbron
CG
, et al
.
Evaluating robustness and sensitivity of the NanoString technologies nCounter platform to enable multiplexed gene expression analysis of clinical samples
.
Cancer Res
2015
;
75
:
2587
93
.
34.
Asaithamby
A
,
Hu
B
,
Chen
DJ
.
Unrepaired clustered DNA lesions induce chromosome breakage in human cells
.
Proc Natl Acad Sci U S A
2011
;
108
:
8293
8
.
35.
Castellanos-Martin
A
,
Castillo-Lluva
S
,
Saez-Freire Mdel
M
,
Blanco-Gomez
A
,
Hontecillas-Prieto
L
,
Patino-Alonso
C
, et al
.
Unraveling heterogeneous susceptibility and the evolution of breast cancer using a systems biology approach
.
Genome Biol
2015
;
16
:
40
.
36.
Riester
M
,
Wei
W
,
Waldron
L
,
Culhane
AC
,
Trippa
L
,
Oliva
E
, et al
.
Risk prediction for late-stage ovarian cancer by meta-analysis of 1525 patient samples
.
J Natl Cancer Inst
2014
;
106
:
dju048
.
37.
Piskol
R
,
Huw
L
,
Sergin
I
,
Kljin
C
,
Modrusan
Z
,
Kim
D
, et al
.
A clinically applicable gene-expression classifier reveals intrinsic and extrinsic contributions to consensus molecular subtypes in primary and metastatic colon cancer
.
Clin Cancer Res
2019
;
25
:
4431
42
.
38.
Wichmann
G
,
Rosolowski
M
,
Krohn
K
,
Kreuz
M
,
Boehm
A
,
Reiche
A
, et al
.
The role of HPV RNA transcription, immune response-related gene expression and disruptive TP53 mutations in diagnostic and prognostic profiling of head and neck cancer
.
Int J Cancer
2015
;
137
:
2846
57
.
39.
Leemans
CR
,
Snijders
PJF
,
Brakenhoff
RH
.
The molecular landscape of head and neck cancer
.
Nat Rev Cancer
2018
;
18
:
269
82
.
40.
Mehanna
H
,
Robinson
M
,
Hartley
A
,
Kong
A
,
Foran
B
,
Fulton-Lieuw
T
, et al
.
Radiotherapy plus cisplatin or cetuximab in low-risk human papillomavirus-positive oropharyngeal cancer (De-ESCALaTE HPV): an open-label randomised controlled phase 3 trial
.
Lancet
2019
;
393
:
51
60
.
41.
Matsumoto
K
,
Nishimura
M
,
Onoe
T
,
Sakai
H
,
Urakawa
Y
,
Onda
T
, et al
.
PARP inhibitors for BRCA wild type ovarian cancer; gene alterations, homologous recombination deficiency and combination therapy
.
Jpn J Clin Oncol
2019
;
49
:
703
7
.
42.
Mirza
MR
,
Monk
BJ
,
Herrstedt
J
,
Oza
AM
,
Mahner
S
,
Redondo
A
, et al
.
Niraparib maintenance therapy in platinum-sensitive, recurrent ovarian cancer
.
N Engl J Med
2016
;
375
:
2154
64
.
43.
Helleday
T
.
Polθ inhibitors unchained
.
Nat Cancer
2021
;
2
:
581
3
.

Supplementary data