INTRODUCTION
Repair of double-strand breaks (DSBs), the most toxic and mutagenic form of DNA damage, involves multiple proteins, but the meiotic recombination 11 homolog (MRE11) complex with RAD50 adenosine triphosphatase (ATPase) and NBS1 phospho-binding protein (MRN) complex plays a central initiating and orchestrating role (
1,
2). At DSBs, the variant histone H2AX, which is inserted within ~2 to 25% of nucleosome core particles throughout the genome, is constitutively phosphorylated on tyrosine-142 (Y142) by the adenosine triphosphate (ATP)–dependent chromatin remodeling complex Williams syndrome transcription factor (WSTF) (
3). In response to DNA damage, serine-139 undergoes phosphorylation (pS139, also known as γH2AX) on both sides of a DNA DSB by phosphatidylinositol 3-kinase (PI3K)–related kinases (
4), thus providing a critical signal and platform for recruitment of DNA damage repair (DDR) proteins (
5). Early DNA damage signaling is dominated by H2AX dually phosphorylated on Y142 and S139 (pH2AX). Yet, around 2 hours after DNA damage, a switch to the γH2AX state of S139 phosphorylation is required to complete repair (
3). The transcriptional coactivator eyes absent (EYA) phosphatase facilitates Y142 dephosphorylation (
6). MRN is recruited to DSB sites via multiple independent processes. Many cell-based assays show that γH2AX plays a pivotal role in MRN recruitment to chromatin through association of NBS1-MDC1-H2AX (
7–
10). The NBS1 FHA domain interacts with phosphorylated mediator of DNA damage checkpoint protein 1 (MDC1), which, in turn, binds to γH2AX, establishing a mode of MRN recruitment to DNA damage sites (
11,
12). Rad17, through its association with NBS1, also indirectly recruits MRE11 to DSBs at an early stage, independent of MDC1 (
13). Evolutionarily conserved complement component 1 Q subcomponent-binding protein (C1QBP) interacts with MRE11-RAD50 to form a complex without NBS1, stabilizing MRE11 but preventing its DNA binding (
14). Thus, multiple processes may enrich MRE11 at DSBs, including some not involving MRN complex. Overall, these data show that timely MRE11 recruitment dictates DSB repair initiation and efficiency as critical for cell survival and genome stability (
2).
Growth factor receptor–bound protein 2 (GRB2) is a metazoan adapter protein essential for receptor tyrosine kinase (RTK)–induced RAS/mitogen-activated protein kinase (MAPK) activation (
15,
16). Functionally, GRB2 acts in early signaling complexes (ESCs): Within seconds after RTK activation, a GRB2-SOS complex is recruited by the receptor to the membrane for RAS activation (
17). GRB2 is so essential for RTK-induced RAS activation that GRB2-knockout (KO) mice are embryonically lethal due to defective endodermal cell differentiation and epiblast formation (
18). Although the nuclear localization of GRB2 was first reported over 20 years ago (
19,
20), no function was assigned. Based primarily upon inferential evidence, a GRB2-PTEN signaling axis was recently proposed to act in DNA repair and genomic stability with possible implications in DDR pathways in HeLa and 293T cells (
21). Yet, a substantial presence and role of GRB2 in the nucleus and its specific molecular and cellular functions remained uncertain. We therefore set out to define and test the molecular basis for GRB2 as a bona fide DDR protein, including possible nuclear activities, as a foundation for future GRB2 research.
Here, investigation of GRB2 in multiple human cells and tissues provides comprehensive and systematic evidence of its nuclear localization. Whereas we found that GRB2 depletion affects DNA repair, our direct measurements found limited correlation between PTEN and GRB2 in DNA repair. Instead, our results define a previously unidentified nucleoplasmic GRB2-MRE11 (GM) complex as an MRE11 pool possibly distinct from its canonical MRN complex (
2). Moreover, GRB2 preferentially promotes homology-directed repair (HDR) and suppresses alternative end joining (Alt-EJ). Furthermore, GRB2-KO was synthetically lethal to the loss of poly(adenosine 5′-diphosphate–ribose) polymerase (PARP) function, and GRB2-mediated MRE11 recruitment to phosphorylated H2AX was indispensable for timely HDR. At DNA damage sites, MRE11 release enabling efficient HDR was regulated by GRB2 ubiquitination on lysine-109 (K109). Release was controlled by the E3 ubiquitin ligase retinoblastoma binding protein 6 (RBBP6), which acts in maintaining common fragile site stability (
22). Mutant cell lines with GRB2-KO, reconstitution of GM interaction disrupting GRB2 mutants, and RBBP6 depletion showed HDR defects. This GM molecular and cellular axis for efficient HDR was further supported by data from patients with breast cancer. In our The Cancer Genome Atlas (TCGA) analyses, high
GRB2 expression showed worse survival only in HDR-proficient patients with high
MRE11 expression. Provocatively, immunohistochemical (IHC) analysis of human normal and breast cancer tissues supported correlation between cancer progression stage and nuclear GRB2 (nGRB2), where high levels of nGRB2 occur in late-stage patients. These findings suggest that
GRB2 and
MRE11 coexpression levels merit testing as a prognostic biomarker in HDR-proficient patients, paving the way to identify patient groups without
BRCA mutations who may favorably respond to PARP inhibitor (PARPi). Overall, we find an unexpected GRB2 function in timely and robust recruitment and ubiquitination-regulated release of MRE11 that promotes HDR and suppresses Alt-EJ, suggesting that the GM complex acts in maintaining genome integrity.
DISCUSSION
DNA repair and genome fidelity critically depend on timely recruitment and regulated release or handoffs of damage excision enzymes, as initially shown for APE1 in base excision repair (BER) (
40). For DSBs, our collective findings unveil GRB2 adapter complex for recruitment and release of MRE11 for efficient HDR. GRB2’s cytoplasmic role in proliferative RAS/MAPK/ERK (extracellular signal–regulated kinase) kinase pathway activation is thus unexpectedly complemented by a nuclear role in targeting MRE11 for DSB repair, as logically associated with replication and proliferative stress. The GRB2-SH2 domain binds MRE11 through a unique binding interface tailored to GRB2 DDR function. As RPA70 binds to a conserved OB1 fold in the nSH3 domain of GRB2, which occurs in many DNA repair proteins [e.g., BRCA2, BLM-complex proteins, and ligase-1 (
41)], further studies may test and uncover GRB2 as an adapter for recruitment of other DNA repair proteins. Currently, GRB2 targeting of MRE11 and RPA to γH2AX for efficient HDR has apparent analogies to recruitment of NEIL1 glycosylase to oxidation-susceptible open chromatin sites for efficient BER, as revealed by the notable correspondence between NEIL1 occupancy and mutation rates along the genome (
42).
GRB2-KO cells show an immediate measurable downregulation of DSB repair activities except for Alt-EJ. This could partly reflect loss of RAS/MAPK signaling in GRB2-KO cells, but reconstitution with separation-of-function GRB2 mutant suggests otherwise. Also, we found that the adapter protein XRCC1, which forms a repair complex with MRE11 and POLQ that supports Alt-EJ activity (
35,
36), is up-regulated in GRB2-depleted cells. Thus, the interplay between nuclear and cytosolic GRB2 redistribution and of GRB2-MRE11 promotion of HDR versus XRCC1-MRE11 enabling Alt-EJ merit future investigation. Notably,
K109RGRB2 retains all cytoplasmic GRB2 functions, but it cannot bind to and recruit MRE11 to the DNA damage site. In GRB2-KO cells, elevated Alt-EJ is returned to normal levels, and NHEJ down-regulation is fully rescued by
K109RGRB2 reconstitution.
K109RGRB2 partially rescues SSA. As EXO1 can potentially substitute for MRE11 in DSB repair, EXO1 may enable resection for SSA (
43,
44).
K109RGRB2, which is defective in MRE11 binding, and
K109AGRB2, which is defective in MRE11 release, both failed to rescue HDR, highlighting the role of GM complex recruitment and release for HDR.
In summary, GRB2 nuclear functions revealed here provide mechanistic insight into MRE11 recruitment to DSB sites with functional consequences for HDR and human cancer. Our TCGA database analysis showed that GRB2 affects survival of patients with both HDR proficiency and high MRE11 expression. Furthermore BRCA-proficient patients with high MRE11 and low GRB2 are a noncanonical HDR-deficient group with cancer progression vulnerability. Both GRB2-KD and GRB2-KO cells are HDR deficient, which causes synthetic lethality to PARPi-treated cells. Overall, our results unveil GRB2’s key role in DDR for efficiently recruiting and releasing MRE11 to initiate HDR rather than Alt-EJ with implications for cancer biology and for targeting synthetic lethality in precision medicine anticancer efforts.
MATERIALS AND METHODS
Reagents
Antibodies against p-ERK1/2 (4370S), ERK1/2 (4695S), p-Akt (S473) (4060S), Akt (4685S), ubiquitin (3936S), ubiquitin (3933S), cleaved PARP (Asp214) (5625S), cleaved caspase-3 (Asp175) (9664S), cleaved caspase-7 (Asp198) (8438S), MRE11 (4847S), RPA70 (2267S), H2AX (7631S), p-H2AX (S139) (20E3) (9718S), RAD50 (3427S), NBS1(14956S), XRCC1 (2735S), histone H3 (4499S), and β-actin (3700S) were ordered from Cell Signaling. Antibodies against GRB2 (C-23) (sc-255), glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (sc-47724), E-cadherin (sc-21791), lamin A (sc-71481), CBLL1 (sc-517157), RING2 (sc-101109), RBBP6 (M56) (sc-9962), GST (sc-138), GFP (sc-9996), and SOS1 (sc-10803) were purchased from Santa Cruz Biotechnology. Antibodies against γH2AX (S139) (05-636) and p-H2AX (Y142) (07-1590) were ordered from MilliporeSigma. Antibodies against MRE11 (ab214), NBS1 (ab181729), and RAD50 (ab124682) were purchased from Abcam. Antibodies against Strep-tag (A00626), red fluorescent protein (RFP)–tag (A00682), and Flag-tag (A00187) were ordered from GenScript. Antibody against Rad51 (GTX100469) was from GeneTex. Antibody against CtIP (61141) was from Active Motif. Antibody against BrdU (347580) was from BD Biosciences, and antibody against RBBP6 (NBP1-49535) was from Novus Biologicals. CBLL1 shRNA1 (V2LHS_157953), CBLL1 shRNA2 (V3LHS_319113), RING2 shRNA1 (V2LHS_188269), RING2 shRNA2 (V3LHS_640474), RBBP6 shRNA1 (V2LHS_255063), RBBP6 shRNA2 (V3LHS_353550), MRE11 shRNA1 (V2LHS_202468), MRE11 shRNA2 (V2LHS_202464), PTEN shRNA1 (V2LHS_231477), PTEN shRNA2 (V2LHS_192536), and GRB2 shRNA (V2LHS_137365) were ordered from Dharmacon. Flag-MRE11, Flag-RAD50, and Flag-H2AX plasmids were purchased from OriGene. Mutants H2AX S139A, H2AX Y142, GRB2 K109R, GRB2 K109A, and GRB2 DF/AA were generated by site-directed mutagenesis. GFP-MRE11 plasmid was a gift from S. Jackson (University of Cambridge, UK). NBS1-GFP plasmid was a gift from J. Lukas (University of Copenhagen, Denmark). GFP-RPA70 plasmid was provided by M. S. Wold (University of Iowa, USA). Flag-XRCC1 and Flag-POLQ plasmids were provided by S. Mitra (Houston Methodist, USA). Phosphorylated EGFR synthetic peptide FLPVPE(pY)INQSVPKR, MRE11 peptide PWVNYQDGNLN, tyrosine-phosphorylated MRE11 peptide PWVN(pY)QDGNLN, phosphorylated H2AX (pH2AX) peptide PSGGKKATQA(pS)QEY, phosphorylated H2AX (pH2AX) peptide PSGGKKATQASQE(pY), dual-phosphorylated H2AX (pH2AX) peptide PSGGKKATQA(pS)QE(pY), and GRB2 OB1 peptide KYDFKATADDELSFKRG were purchased from GenScript. GFP-trap, Flag agarose beads, and Strep-Tactin beads were purchased from Chromotek (gta-10), MilliporeSigma (A2220), and MilliporeSigma (71592), respectively.
Cell culture
Human embryonic kidney (HEK) 293T, NIH3T3, and A431 cells [purchased from the American Type Culture Collection (ATCC)] and MEFs and H2AX-KO MEFs (gifts from A. Nussenzweig, Center for Cancer Research, National Cancer Institute) were maintained in Dulbecco’s modified Eagle’s high-glucose medium. HeLa cells (ATCC) were maintained in RPMI 1640 medium, and HAP1 cells (Horizon Discovery) were maintained in Iscove’s modified Dulbecco’s medium. U2OS cells (ATCC) were maintained in McCoy’s 5A medium. Media were supplemented with 10% (v/v) fetal bovine serum (FBS) and 1% antibiotic/antimycotic (Lonza) in a humidified incubator with 10% CO
2. Stable cells containing Strep-tagged GRB2 or GRB2 mutants were produced as described previously (
45). The I-SceI reporter assay (a gift from J. Stark, City of Hope) was performed as described previously (
33). The sulforhodamine B cytotoxicity and colony formation assays were performed as described previously (
46).
CRISPR-Cas9–mediated GRB2-KO
GRB2 was knocked out by using the CRISPR-Cas9 system according to the protocol described previously (
47,
48). The guide RNA (gRNA) targeting GRB2 was designed using an online tool available at
http://crispr.mit.edu. DNA primers (forward primer, 5′-CACCGGAGCCGGAAGTCTTCCTC-3′; reverse primer, 5′-AAACGAGGAAGTACTTCCCGGCTCC-3′) for the GRB2 gRNA and reverse complement sequence plus adapters needed for ligation were synthesized from IDT and ligated into the LentiCRISPR v2 (Addgene #52961) or LentiCRISPR v2-Blast (Addgene #83480). Correct insertion of GRB2 gRNA was confirmed by sequencing the constructs. Lentivirus particles were generated in HEK293T cells. Cells were infected with lentivirus and then selected with puromycin or blasticidin. Single clones were picked up, and GRB2-KO was verified by Western blotting. Mixture pools of different clones were used for experiments.
Western blots
Cells were seeded onto 10-cm dishes and grown for at least 24 hours before experiments. Cells were lysed with radioimmunoprecipitation assay (RIPA) buffer [20 mM tris-HCl (pH 7.5), 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% IGEPAL, 1% sodium deoxycholate, 2.5 mM sodium pyrophosphate, 1 mM β-glycerophosphate, 1 mM Na3VO4, and leupeptin (1 μg/ml)] supplemented with Protease Inhibitor Cocktail Set III (EMD Millipore) to obtain total proteins. Cell fractionation was done using Subcellular Protein Fractionation Kit for Cultured Cells (Thermo Fisher Scientific) according to the manufacturer’s instructions.
Protein expression and purification
Expression and purification of GRB2 from bacteria have been described previously (
31). The K109R and K109A mutant GRB2 was generated by site-directed mutagenesis and purified in a manner similar to the WT. Human MRE11 core (residues 1 to 411) purification was described previously (
30). The generation of individual GRB2 domains was as described before (
49). The purification of RPA70N (1-120) was performed according to a published procedure (
50).
Immunofluorescence
Designated cells were grown on coverslips, fixed with the addition of 4% (w/v) paraformaldehyde (pH 8.0) with or without pre-extraction, and washed four times with phosphate-buffered saline (PBS; pH 8.0). After permeabilization with 0.5% Triton X-100 on ice for 5 min, cells were washed three times with PBS and incubated with blocking buffer [PBS, 3% bovine serum albumin (BSA), 5% FBS, and 0.5% Triton X-100] for 2 hours at room temperature or overnight at 4°C. Following a further three washes with PBS, cells were incubated with primary antibody overnight in PBS, 3% BSA, and 0.5% Triton X-100. Cells were then washed five or six times with PBS and incubated with the fluorescence-conjugated secondary antibody for 2 to 3 hours. Following another PBS wash, coverslips were mounted onto a slide with mounting medium (0.1% p-phenylenediamine and 75% glycerol in PBS at pH 7.5 to 8.0). For pre-extraction before fixation, cells were treated with cytoskeleton buffer [10 mM Pipes (pH 6.8), 100 mM NaCl, 300 mM sucrose, 3 mM MgCl2, 1 mM EGTA, and 0.5% Triton X-100] for 5 min on ice and followed by washing with stripping buffer [10 mM tris-HCl (pH 7.4), 10 mM NaCl, 3 mM MgCl2, 1% Tween 20, and 0.25% sodium deoxycholate] for 5 min on ice. Cells were imaged using either a Leica SP5 II or Zeiss LSM710 confocal microscope.
UV laser micro-irradiation
UV-LMI was performed as described previously (
51): Cells were grown on 35-mm glass-bottom dishes (Matsunami), micro-irradiated with a MicroPoint ablation system (Photonics Instruments) with 60% laser output, and visualized with a Nikon Eclipse TE2000-U inverted microscope.
Alkaline comet assay
Untreated and IR-treated cells embedded in low–melting point agarose were layered over a slide with a thin film of agarose. Cells were then lysed overnight at 4°C in lysis buffer [10 mM tris-HCl (pH 10), 2.5 M NaCl, 100 mM EDTA, 5% dimethyl sulfoxide (DMSO), and 1% Triton X-100]. After equilibration in electrophoresis buffer [300 mM NaOH, 1 mM EDTA (pH > 13)] for 40 min at 4°C, electrophoresis was performed to identify DNA damage by the migration of small fragments of broken DNA away from the nucleus, forming a “tail” to the dense circular comet head of undamaged DNA. Nuclear DNA was visualized by staining with SYBR Gold solution, and images were captured using a Leica SP5 II microscope. The tail moment of each cell was analyzed by OpenComet software.
Annexin V/PI staining assay
To determine the percentage of olaparib-induced apoptotic cells, an annexin V/PI staining assay was used according to the standard protocol (Life Technologies). Briefly, cells were cultured and treated with DMSO or 50 μM olaparib for 72 hours. Then, cells were trypsinized, washed, and resuspended in annexin-binding buffer. Cells were labeled by adding annexin V–Alexa Fluor 647 and PI to each sample. After incubation for 15 min at room temperature in the dark, samples were analyzed on a flow cytometer (FACSCanto II, BD Biosciences) for the detection of annexin V– and PI-positive subpopulations. DMSO-treated cells were used as a control, and each experiment was performed in triplicate. Further data analysis was performed with FlowJo V10 software.
Microscale thermophoresis
MST was measured using the Monolith NT.115 system (Nano Temper). Proteins were fluorescently labeled with Atto 488 according to the manufacturer’s protocol. Labeling efficiency was determined to be 1:1 (protein:dye) by measuring the absorbance at 280 and 488 nm. A solution of peptides or proteins in 0.01 M Hepes (pH 7.4), 0.15 M NaCl, and 0.005% (v/v) Surfactant P20 was serially diluted, typically from about 100 μM to 30 nM, in the presence of 100 nM labeled protein. The samples were loaded into silica capillaries (Polymicro Technologies) after incubation at room temperature for 15 min. Measurements were performed at 22°C using 20% light-emitting diode power and 40% infrared laser power. Measurements were also carried out using 20 and 60% infrared laser power for comparison. Data analyses were performed using Nano Temper Analysis software using the Kd curve fitting function. Raw data were exported and fitting curves were generated using Prism 8 (GraphPad Software) for presentation.
Molecular docking
We used the available GRB2 structure from the PDB (PDB ID, 1GRI). The structure was prepared for docking by eliminating water and cofactors; a monomeric unit was used. The GRB2 coordinates were prepared and minimized using Protein Preparation Wizard (Schrödinger Suite 2019-2). We aimed the molecular docking to the K109 site. We docked ligands (in this specific experimental peptide) with length ≤20 Å with cubic box dimensions of 15 Å. We also generated a mutant of GRB2, K109R, using the function “simply mutate” in Coot model-building software. In this structure, we aimed the docking to R109. We used the peptide docking utility implemented in the Schrödinger Suite; the peptides were added as a sequence in the specific section “Define peptide to dock.” We used the following peptide-derived sequences for MRE11, EGFR, and γH2AX: Mre11-PWVNYQDGN, EGFR-NPVYHNQPL, and H2AX-SQEY, respectively. For the docking experiments, we used Glide (Schrödinger Release 2019-2). The best docking results were selected and compared in superimposed sessions in the PyMOL molecular graphics system (PyMOL 2.0; Schrödinger).
IHC and staining quantification
Breast cancer tumor tissue array slides (BC081116d; Biomax), normal human tissue array slides (MNO961; Biomax), and the paraffin-embedded tissue sections of different mouse tissues were collected for determination of GRB2 distribution. Sections were dewaxed and rehydrated following a standard dewaxing protocol. Then, the samples were exposed to 10 mM citric acid buffer (pH 6.0 for 20 min at 105°C) for heat antigen retrieval and exposed to 3% H2O2 for 12 min to block endogenous peroxidase activity. Subsequently, the samples were blocked with goat serum for 1 hour and incubated with anti-GRB2 primary antibody (1:100) overnight. The VECTASTAIN ABC HRP Kit (peroxidase, rabbit immunoglobulin G) and DAB Substrate Kit (peroxidase, horseradish peroxidase) were used to develop color, followed by nuclear counterstaining with hematoxylin. Last, the slides were mounted and visualized. Representative regions were selected and photographed. Two pathologists were tasked with evaluating IHC staining, and the H-score method was applied for calculating the staining score of each sample. In brief, each sample was assessed by staining intensity (negative, 0; weak, 1; moderate, 2; strong, 3) and staining extent (0 to 100%), with H-score = [1 × (staining extent of weak staining) + 2 × (staining extent of moderate staining) + 3 × (staining extent of strong staining)]. Three samples were used to quantify H-score for each tissue, and the average of H-score and SD for all the cases were calculated and presented.
Mass spectrometry
To identify proteins associated with GRB2 in the nucleus, GRB2 was immunoprecipitated with Strep-Tactin beads (MilliporeSigma #GE28-9355-99) from the nuclear extracts of 293T cells expressing Strep-GRB2 or Strep-Ctrl. Three experimental samples and three matched negative controls were prepared independently for MS. Using the MyriMatch and Sequest search engines, the MS spectra were matched to the human protein database. Scaffold 3 was used to analyze the fold enrichment values of each sample compared to the matched negative control sample. Proteins enriched at least 1.5-fold with a P value of less than 0.05 (Fisher test) as computed by Scaffold 3 were selected and listed. The proteins in the list were used to conduct the pathway and network analysis by BINGO in Cytoscape 3.7.2. Hypergeometric tests were performed and corrected by Benjamini and Hochberg false discovery rate correction. A P value of less than 0.05 was considered significant.
Statistical analysis
All statistical analyses were performed using GraphPad Prism 8 using multiple t test function.
TCGA data analysis
We retrieved the gene expression and clinical patient data to estimate the Kaplan-Meier survival curves from TCGA using the TCGA-Assembler suite available at
https://github.com/compgenome365/TCGA-Assembler-2. For gene expression analyses between tumor and matched controls, we selected datasets with at least 10 matched controls. Data were plotted with the R packages ggpllot2 and ggpubr. Survival curves and hazard ratios were obtained with the survminer, survival, and dplyr R packages. Signature 3 mutations were obtained from
https://synapse.org/#!Synapse:syn11801497.
Acknowledgments
We thank J. Chen for UV-LMI instrument access, K. Schlacher for discussion, and L. S. Warden and T. M. Link for protein purifications. Editorial assistance was provided by MD Anderson’s Editing Services, Research Medical Library. Funding: This research was supported by the University Cancer Foundation via the Institutional Research Grant program at The University of Texas MD Anderson Cancer Center (to Z.A.); by Cancer Prevention Research Institute of Texas (CPRIT) grants RP180813 and RP130397; and by NIH grants R01 CA200231, P01 CA092584, R35 CA220430, and 1S10OD012304-01. J.A.T.’s efforts are also supported by a Robert A. Welch Chemistry Chair. Author contributions: Z.A. and J.A.T. directed the research; Z.Y., J.A.T., and Z.A. conceived the research plan; Z.Y., S.X., Y.S., and Z.A. performed the experiments; Z.Y., S.X., Y.S., J.A.T., and Z.A. designed experiments and analyzed results; Y.S. and A.B. contributed to MS and TCGA data analysis; D.M. and C.-L.T. performed computational docking and analysis; A.S., Q.S., G.P., P.G.L., D.E.J., and B.W. provided critical reagents and resources; and Z.Y., S.X., J.A.T., and Z.A. wrote the paper. Competing interests: J.A.T., Z.A., and Z.Y. are inventors on a patent application related to this work filed by The University of Texas MD Anderson Cancer Center. The authors declare no other competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.