Abstract
Background/Aim: Tyrosine kinases have crucial functions in cell signaling and proliferation. The phosphatidylinositol 3-kinase (PI3K) pathway is frequently deregulated in human cancer and is an essential regulator of cellular proliferation. We aimed to determine which tyrosine kinases contribute to resistance elicited by PI3K silencing and inhibition. Materials and Methods: To mimic catalytic inactivation of p110α/β, specific p110α (BYL719) and p110β (KIN193) inhibitors were used in addition to genetic knock-out in in vitro assays. Cell viability was assessed using crystal violet staining, whereas cellular transformation ability was analyzed by soft-agar growth assays. Results: Activated zeta chain of T-cell receptor-associated protein kinase 70 (ZAP70) generated resistance to PI3K inhibition. This resistance was via activation of the Janus kinase/signal transducer and activator of transcription 3 (JAK/STAT3) axis. We demonstrated that activated ZAP70 has a high transforming capability associated with the formation of malignant phenotype in untransformed cells and has the potential to be a tumor-initiating factor in cancer cells. Conclusion: ZAP70 may be a potent driver of proliferation and transformation in untransformed cells and is implicated in resistance to PI3K inhibitors in cancer cells.
Keywords: PI3K, cell signaling, tyrosine kinases, growth resistance, ZAP70
Protein phosphorylation causes activation of signal transduction pathways, which are crucial for many biological processes (1). The most well-characterized protein kinases in the human genome are protein serine/threonine kinases and protein tyrosine kinases (2). High levels of tyrosine phosphorylation correspond to enhanced proliferative state (3).
The phosphatidylinositol 3-kinase (PI3K) pathway is one of the most activated signaling pathways in human cancer. Upon activation of receptor tyrosine kinases (RTKs), recruitment of PI3K to the cell membrane occurs via coupling the Src homology-2 (SH2) domains of p85 with tyrosine phosphorylated residues of activated receptors (4). Oncogenic RAS and other tyrosine kinases can cause PI3K activation. This pathway is also a crucial regulator of cell proliferation, survival and motility (5). Although many findings support the idea that targeting the PI3K pathway is one of the most promising approaches in cancer therapy (4,6), inhibition of class I PI3K activity by using pharmacological inhibitors leads to retardation of cell proliferation rather than induction of apoptosis, bringing about a problem with generation of drug resistance in the longer run (7).
Previous studies have identified the contribution of several activated tyrosine or serine/threonine kinases to growth compensation upon PI3K inhibition (8-10). In addition to this, pharmacological inhibition of the PI3K pathway in cancer cells leads to up-regulation of receptor tyrosine kinases through feedback loops that circumvent baseline level of PI3K inactivation. These compensatory mechanisms include forkhead box O (FOXO) transcription factors along with mitogen-activated protein kinase (MAPK)-dependent pathways (11,12). Elucidation of the primary mechanisms of PI3K resistance might be a significant therapeutic approach for tumorigenesis (6). As an example, Pim-1 serine/threonine kinase (PIM1) regulates cell proliferation and survival that might lead to reduction of therapeutic drug efficacy by promoting the Janus kinases/signal transducer and activator of transcription proteins (JAK/STAT) signaling pathway in MYC proto-oncogene (MYC)-driven cancer (13-16), particularly in breast cancer.
Zeta chain of T-cell receptor associated protein kinase 70 (ZAP70) is a cytoplasmic tyrosine phosphoprotein. It belongs to spleen tyrosine kinase (SYK) family of nonreceptor tyrosine kinases (17). ZAP70 has a crucial role in T-cell receptor signaling activation and T-cell development (18). Based on structural studies, T-cell activation by antigens results in phosphorylation of LCK proto-oncogene, Src family tyrosine kinase (LCK) and binding of SH2 domains of ZAP70 to tyrosine phosphorylated immunoreceptor tyrosine-based activation motifs which leads to its recruitment to the membrane (19). Existing literature has also described non-immunological functions of ZAP70. ZAP70 was shown to promote migration and invasion of prostate cancer cell lines (20). Moreover, in mammalian oocytes and embryonic stem cells, ZAP70 has a role in maintaining stemness and differentiation by regulating the JAK–STAT3–MYC signaling axis (21). Additionally, ZAP70 was identified as a prognostic marker in colorectal cancer for radiation response (22). In this study, we wanted to identify potential tyrosine kinases that contribute to generating resistance to PI3K inhibition. ZAP70 was identified as a mediator of resistance to PI3K ablation in our genetic screen. We then aimed to understand the potential role of activated ZAP70 in tumor initiation, carcinogenesis and resistance to PI3K inhibition.
Materials and Methods
Generation of cell lines expressing activated tyrosine kinase and cell culture. p110αflox/flox; p110βflox/flox Mouse embryonic fibroblasts (MEFs) and SV40 large T-antigen-expressing human embryonic kidney 293 (HEK293T) cells, T47D and MCF7 luminal A type breast cancer cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Sigma–Aldrich, St Louis, MO, USA) with 4.5 g/l D-glucose and 1 mM L-glutamine (Sigma–Aldrich) supplemented with 8% fetal bovine serum (FBS, Biowest, Nuaillé, France), and 100 IU/ml penicillin and 100 μg/ml streptomycin. Human telomerase reverse transcriptase-immortalized retinal pigment epithelial-1 (RPE1-hTERT) cells were grown in DMEM/F-12, supplemented with 8% FBS, 1 mM L-glutamine and 100 IU/ml penicillin and 100 μg/ml streptomycin (Thermo Fisher Scientific, Waltham, MA, USA). The non-tumorigenic human epithelial breast cells MCF10A and human mammary epithelial cells (HMECs) were cultured in DMEM/F-12, supplemented with 4% FBS, penicillin (100 IU/ml), streptomycin (100 μg/ml), L-glutamine, 10 ng/ml epidermal growth factor (Sigma–Aldrich), 10 μg/ml insulin (Sigma–Aldrich) and 0.5 μg/ml hydrocortisone (Sigma–Aldrich). Jurkat cells were maintained in RPMI 1640 medium supplemented with 8% FBS, 1 mM L-glutamine, and 100 IU/ml penicillin/streptomycin. All cell lines were incubated at 37˚C with 5% CO2. The Jurkat and hTERT-RPE1 cells were purchased from the American Type Culture Collection (Manassas, VA, USA). All the other cell lines have previously been validated and characterized (23).
To identify a tyrosine kinase that potentially contributes to generating resistance to PI3K abrogation, we screened an activated tyrosine kinase library in MEFs which are genetically engineered with LoxP sites, inserted into the first exons of PIK3CA and PIK3CB (23,24). The LoxP sequences in these cells can be targeted by Cre recombinase or LacZ-expressing adenoviruses (Ad/Cre, Ad/LacZ; The University of Iowa, Viral Vector Core, Iowa City, IA, USA) for excision of the targeted exons of PIK3CA and PIK3CB to generate p110α/β double knock-out cells. The activated tyrosine kinase library has 73 open-reading frames (Table I). This pooled tyrosine kinase library consists of modified receptor and non-receptor tyrosine kinases. The activated tyrosine kinase library constructs have C-terminal dimerization tag; ETS variant transcription factor 6 (ETS6) dimerization domain (TEL) and are tagged with FLAG epitope (25). The tyrosine kinase pools were stably expressed in MEFs and screened using crystal violet growth assay. The TEL-RTK library was stably expressed in cells using the transfection protocol of Lipofectamine 2000 Reagent (Invitrogen, Carlsbad, CA, USA). To generate HEK293T cells transiently expressing the activated tyrosine kinase library, Lipofectamine 3000 Reagent (Invitrogen) was used. The pWzl green fluorescent protein-expressing (#12269; Addgene, Watertown, MA, USA) and myristoylated AKT serine/threonine kinase 1 (MYR-AKT1)-expressing (#15294; Addgene) MEFs were generated by retroviral transduction. The H-Ras G12V construct (#9051; Addgene) was stably expressed in MEFs and RPE1-hTERT cells using the protocol of Lipofectamine 2000 reagent (Invitrogen).
Cell viability assays and inhibitor treatment. The cells were seeded into 12-well or 6-well plates at a confluency of 30-40% and then treated with increasing concentrations of p110α-specific inhibitor BYL719 (0.3, 0.5 or 1 μM, Alpelisib; Selleck Chem, Houston, TX, USA) or p110β-specific inhibitor KIN193 (1 μM, AZD6482; Selleck Chem) or with both inhibitors combined (0.3, 0.5 or 1 μM) for approximately 7-10 days. The inhibitors were supplemented in reduced serum medium (4% FBS) to maintain physiologically relevant levels of growth factors. When control cells reached about 80% confluency, cells were fixed in 10% acetic acid and 10% ethanol for at least 24 hours at room temperature. Cells were then washed with 1× non-sterile phosphate-buffered saline and stained with 0.4% crystal violet and 20% ethanol for 1 hour at RT. Cells were washed with distilled water and air-dried. DNA intercalating dye was removed with 1 ml 10% acetic acid. Plates were incubated in destaining solution for 1 hour at room temperature. The optical density of the solution was measured at 595 nm.
Adenovirus/Cre recombinase transduction of MEFs. MEF cell lines were seeded in 6-well plates in duplicates at a density of 1×105 cells per well. The day after seeding, cells were treated with Ad/Cre at a multiplicity of infection of 75 in medium containing 2% FBS and incubated for 6-8 h. Then the medium was removed and growth medium was added into cells. The infection was repeated three times to facilitate near-complete excision of the first exons of endogenous PIK3CA and PIK3CB. The adenovirus/Cre-treated cells were seeded into 12-well plates at a density of 2×103 cells per well and cultured for 6 days for cellular viability assay.
Soft-agar assay. The soft-agar assay was performed in 6-well plates as triplicates. To form the bottom layer, 2.5 ml DMEM containing 2.25% low-melting agarose (Sigma–Aldrich) mixture was poured into each well and were solidified at room temperature for at least 30 minutes. Then 2×105 cells were resuspended in 9 ml DMEM solution and mixed with 3.2 ml 2.25% low-melting agarose. Two millimeters of this mixture was poured into each well and allowed to solidify. A total of 2-3 ml complete growth medium was added to each well and refreshed every 4-5 days. After 4 weeks, the colonies formed were stained with 0.005% crystal violet and 10% ethanol solution and counted under a microscope.
Antibodies and western blotting. The protein lysates were prepared from MEFs, RPE1-hTERT, MCF10A and HEK293T cell pellets. Then 20-30 μg of protein was resolved by sodium dodecyl sulphate–polyacrylamide gel electrophoresis and transferred to nitrocellulose membranes (Bio-Rad, Hercules, CA, USA). The membranes were blocked with 5% bovine serum albumin (Serva, Heidelberg, Germany) in TBS and then incubated with primary antibodies in TBS with 5% bovine serum albumin and 0.1% Tween. The following primary antibodies were used: PI3K p110α (C73F8; Cell Signaling, Danvers, MA, USA), PI3K p110β (C-8; Santa Cruz Biotechnology, Dallas, TX, USA), phospho-STAT3 (Tyr705) (Cell Signaling), ZAP70 (1E7.2) (Santa Cruz Biotechnology), phospho-ZAP70 (Tyr493/Syk (Tyr526) (Cell Signaling), STAT3 (Cell Signaling), phospho- AKT (Ser473) (Cell Signaling), phospho-AKT (Thr308) (D25E6) (Cell Signaling), phospho-p44/42 extracellular signal-regulated protein kinase 1/2 (ERK1/2) (Thr202/Tyr204) (Cell Signaling), phospho-p70 S6 kinase (Thr389) (Cell Signaling), S6 ribosomal protein (Cell Signaling), glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Cell Signaling), β-actin (clone AC74; Sigma–Aldrich) phospho-tyrosine 4G10 (Sigma–Aldrich).
In-silico analysis. Gene Expression Profiling Interactive Analysis (GEPIA; http://gepia2.cancer-pku.cn) is an interactive platform for analyzing RNA sequencing expression data of 9,736 tumors and 8,587 normal samples from The Cancer Genome Atlas (TCGA) and the Genotype Tissue Expression (GTEx) projects. Analysis of survival was conducted by selecting disease-free survival as the endpoint with a median cutoff of 50% for each gene of interest to generate Kaplan–Meier survival curves. In addition, survival curves were generated using KM plotter by restricting the analysis of progression-free survival by grade and chemotherapy treatment (26). University of Alabama Cancer web portal (UALCAN) database was used to depict the comparison of mRNA expression levels between tumor and normal TCGA datasets (27). The cBio Cancer Genomics Portal was used to indicate the alteration frequency of ZAP70 in various type of cancer using US National Cancer Institute (NCI)-60 cancer cell line studies (28).
Statistical analysis. Two-way analysis of variance and paired t-test were used for differential comparison between variables. All statistical analyses were performed in GraphPad Prism 8.0 (San Diego, CA, USA). p-Values less than 0.05 were considered significant.
Results
Library screen for identifying factors that generate resistance to PI3K inhibition. Our activated tyrosine kinase library comprises 73 open-reading frames of receptor and non-receptor tyrosine kinases which are activated via C-terminal dimerization domain of the TEL transcription factor (Table I) (25). We stably expressed activated tyrosine kinases in MEFs and generated p110α/β double knock-out cells. As a result of our screening, Pool 9, which included protein tyrosine kinase 2 beta (PTK2B), protein tyrosine kinase 6 (PTK6), ZAP70 and neurotrophic receptor tyrosine kinase 3 (NTRK3), partially restored growth capability upon dual loss of p110α/β (Figure 1A). In our screen, green fluorescent protein-expressing MEFs were used as a negative control, whereas MYR-AKT1-expressing MEFs served as a positive control as MYR-AKT1 localizes to the cell membrane and becomes constitutively activated (29). Then, we extracted proteins from MEFs that expressed FLAG-tagged activated tyrosine kinases in Pool 9 before and after Ad/Cre treatment and performed western blot analysis. Immunoblotting experiments confirmed the presence of immunoreactive bands compatible with the expected sizes of NTRK3, ZAP70 and PTK6, whereas PTK2B expression was not detected (Figure 1B). We also checked for the efficiency of PIK3CA and PIK3CB knock-out in the cells and found that p110α and p110β protein expression was abolished, confirming Ad/Cre targeting (Figure 1C).
TEL-ZAP70 promoted cell proliferation. Cellular growth assays using individual Pool 9 components demonstrated that standalone expression of activated ZAP70 enabled the MEFs to grow efficiently upon Ad/Cre treatment in comparison to those bearing other tyrosine kinases in the pool and to controls (Figure 2A). To understand the role of activated tyrosine kinases in MEFs, we stably transfected these cells with TEL-ZAP70 or TEL-PTK6, using MEFs transfected with MYR-AKT1 as a positive control and pBabe Neo (empty) as a negative control (Figure 2B). We conducted crystal violet cell viability assays. Our results showed that TEL-ZAP70 expression increased the rate of cellular growth to a level similar to MYR-AKT1-expressing MEFs. However, TEL-PTK6 expression in MEFs had no significant effect on cellular growth (Figure 2C).
Investigation of the growth-compensatory potential of activated ZAP70 upon PI3K knockout. In order to confirm the tyrosine kinase library screening, we followed two complementary approaches to inhibit PI3K signaling in MEFs. FirstIy, we took advantage of molecular genetics and induced genetic silencing of PIK3CA as well as PIK3CB, and secondly, we employed pharmacological inhibitors specific for p110α and -β isoforms. The immunoblot analysis demonstrated efficient p110α/β ablation in Ad/Cre-treated MEFs in comparison to controls (Figure 2D). The knockout of p110α/β caused growth retardation in control MEFs. Nevertheless, TEL-ZAP70 expression in MEFs partially restored the impairment of cellular growth imposed by p110α/β knockout, as did MYR-AKT1 expression in the positive control (Figure 2E). Of note, the compensation was more pronounced in TEL-ZAP70-expressing MEFs in comparison to MYR-AKT1-expressing MEFs.
For pharmacological inhibition, we used Food and Drug Administration-approved small-molecule inhibitor of PI3K-p110α, Alpelisib (BYL719), and a specific PI3K-p110β inhibitor, KIN193 (AZD6482) to selectively inhibit kinase activities of p110α and p110β respectively. In general, TEL-ZAP70 MEFs were growth-inhibited to a lesser degree than control MEFs in crystal violet growth assays. We particularly observed a statistically significant reduction of growth inhibition in response to the combination of BYL719 and KIN193 when the cells expressed activated ZAP70 (Figure 2F). These results suggest that activation of ZAP70-mediated signaling induces growth resistance to genetic or pharmacological inhibition of p110α/β.
TEL-ZAP70 promoted anchorage-independent growth in MEFs. MEFs are untransformed cells which are not capable of growing in an anchorage-independent manner or initiating tumors (30,31). To determine whether ZAP70 has the ability to transform cells, we performed soft-agar assays with our MEF lines that stably expressed activated ZAP70, as well as H-Ras G12V mutant, which constitutively activates H-Ras to enhance transforming ability in mesenchymal cells as a positive control (32,33).
In this experiment, expression of activated PTK6, which did not compensate for PI3K knockout previously, was used as an additional negative control (Figure 2A). Our results indicated that although activated PTK6- and pBabe Neo empty vector-expressing MEFs were not able to form colonies, ZAP70 overexpression triggered colony formation on soft agar (Figure 2G, right panel). Moreover, co-expression of TEL-ZAP70 and H-Ras G12V had an additive effect on colony formation, implicating a cooperation of ZAP70-mediated signaling with the MAPK pathway (Figure 2G).
ZAP70 expression leads to activation of the JAK/STAT3 pathway. Receptor and non-receptor tyrosine kinases might activate alternative signaling pathways besides PI3K via phosphorylation of downstream targets (34,35). We hypothesized that the JAK/STAT and MAPK pathways, which have been implicated in functional compensation of PI3K (10,16), might be activated upon ZAP70-induced partial resistance. To assess the molecular mechanisms that promote compensation of growth upon PI3K knockout in activated ZAP70-expressing MEFs, we conducted western blots to analyze activation of relevant pathways. According to our results, the level of STAT3 phosphorylation was elevated in TEL-ZAP70-expressing MEFs (Figure 2H). In addition, STAT3 phosphorylation was also elevated in MYR-AKT1-expressing MEFs, which can be explained by the essential role of STAT3 in PI3K-induced oncogenic transformation (36).
Additionally, a slight increase in ERK1/2 phosphorylation was observed (Figure 2H). However, the level of pS6K did not change. The lack of a change in the level of pS6K for ZAP70-expressing MEFs indicates a PI3K-independent mode of growth compensation. Taken together, these observations indicate that growth compensation mediated by activated ZAP70 might occur via phosphorylation/activation of JAK/STAT or MAPK signaling pathways.
To understand if endogenous ZAP70 expression is present in adherent cell lines, we conducted western blot analysis. To this end, Jurkat CD4+ T-cells and TEL-ZAP70-expressing T47D cells were used as a positive control. The immunoblot showed that all cell lines, except for p110α/β+/+ wild-type MEFs, expressed endogenous ZAP70 to various extents (Figure 3A). LCK tyrosine kinase was found to be an important mediator of phosphorylation-dependent activation of ZAP70 (37). To further illustrate if endogenously expressed ZAP70 is phosphorylated and therefore activated, we examined phosphorylation levels of ZAP70. We observed elevated levels of phospho-ZAP70 in TEL-ZAP70-expressing T47D cells, along with HEK293T and MCF7 cells (Figure 3B). Taken together, we detected expression of endogenous and possibly activated ZAP70 expression in several adherent cell lines, including HEK293T, MCF7 and T47D. Several studies showed LCK expression in several solid cancer types [reviewed in (37)], therefore it seems plausible that LCK phosphorylates and activates ZAP70 in a context-dependent manner.
Functional analysis of TEL-ZAP70 expression in alternative, non-transformed epithelial cell culture models. The morphological characteristics of MEFs are different from epithelial cells as they represent cells of mesenchymal origin. Since most carcinomas are derived from epithelial cells, we used non-transformed RPE1-hTERT cells and non-tumorigenic epithelial human breast cells (MCF10A) to understand the impact of activated ZAP70 on initiation of carcinogenesis. We generated cell lines with stable expression of activated ZAP70 along with positive and negative controls in MCF10A and RPE1-hTERT (Figure 4A). Next, we wanted to analyze the activity of ZAP70-mediated signaling pathways in these cells. Our results showed that the phospho-STAT3 level in both cell lines expressing activated ZAP70 was elevated (Figure 4B). Since activated ZAP70 had little effect on pERK and pS6K levels in our study, these results imply that the mechanistic target of rapamycin kinase (mTOR) and MAPK pathways might not be universally affected by ZAP70 activation. However, JAK/STAT signaling pathway appears to be consistently up-regulated by activated ZAP70 in a broad range of cellular models.
We next assessed whether activated ZAP70 had a positive effect on cell viability upon PI3K inhibition in RPE1-hTERT cells. We used TEL-ZAP70-expressing RPE1 cells along with controls (pBabe Neo and MYR-AKT1-expressing RPE1 cells). Crystal violet growth assays demonstrated that TEL-ZAP70 expression in RPE1-hTERT cells partially compensated growth upon BYL719-mediated inhibition of p110α (Figure 4C).
After determining the importance of activated ZAP70 for cellular proliferation in RPE1-hTERT cells, we wanted to see if TEL-ZAP70 had the ability to form colonies in anchorage-independent growth assays. We performed soft-agar growth assays and found that TEL-ZAP70-expressing RPE1-hTERT cells were able to form colonies (Figure 4D), although constitutively activated AKT was unable to induce colony formation. Taken together, expression of activated ZAP70 in non-transformed epithelial cells can increase their proliferative potential as well as their ability to transform into more malignant phenotypes.
Analyzing the involvement of ZAP70 in solid tumors in silico. ZAP70 activation has been associated with leukemia (38). We performed in silico analyses to determine if it might also be involved in initiation or progression of solid tumors. Analyses of the US National Cancer Institute’s NCI-60 cell lines studied using the cBioPortal platform revealed that, ZAP70 is moderately altered in renal cell carcinomas (9%), and these alterations mostly occur in the form of gene amplification (Figure 5).
Moreover, we investigated ZAP70 mRNA expression changes in tumor vs normal tissue for solid tumors. For this purpose, we used the University of Alabama cancer web portal in UALCAN and found that in many cancer types, such as renal clear-cell carcinoma, renal papillary-cell carcinoma, cervical squamous cell carcinoma, endocervical adenocarcinoma, esophageal carcinoma, prostate adenocarcinoma, stomach adenocarcinoma, uterine corpus endometrial carcinoma, as well as head and neck squamous carcinoma, ZAP70 mRNA levels were found to be higher in tumor tissues in comparison to matched controls (Figure 6).
We wished to expand our observations to other types of cancer. To this end, we performed overall survival analysis in various cancer types. The Kaplan–Meier plots showed that low expression of ZAP70 in gastric cancer correlated with better survival (Figure 7A). Additionally, the Kaplan–Meier curves generated using the GEPIA2 platform demonstrated that low ZAP70 expression correlated with better survival in renal clear-cell carcinoma, low-grade glioma, and uveal melanoma (Figure 7B-D). Next, to delineate whether ZAP70 has a differential effect on distinct stages of tumorigenesis, we plotted stage-specific Kaplan–Meier survival curves. Our results showed that low expression of ZAP70 in stage 2 and stage 3 gastric cancer correlated with better survival (Figure 8A and B). These in silico analyses support our findings implying a causative role for ZAP70 in various types of solid tumor. To further investigate involvement of ZAP70 in solid tumors, we plotted Kaplan–Meier survival curves in ovarian cancer using KM Plotter database. In survival analysis of patients with grade 1-2 ovarian cancer treated with paclitaxel/cisplatin chemotherapy, high ZAP70 expression levels significantly correlated with reduced progression-free survival (Figure 8C). In conclusion, ZAP70 might be a crucial prognostic marker for ovarian carcinogenesis and low ZAP70 expression correlated with better survival in patients treated with chemotherapy for early-stage ovarian cancer.
The role of ZAP70 overexpression in kidney cell lines. According to our in-silico analysis, renal cell carcinomas had the highest frequency of ZAP70 amplifications (Figure 5). We wanted to investigate the impact of ZAP70 overexpression in kidney-derived cell lines to consolidate our findings in previous cellular models. To that end, we used embryonic kidney-derived HEK293T cells and transiently expressed activated ZAP70 along with MYR-AKT1 as a positive control and the empty vector (pBabe Neo). Expression of TEL-ZAP70 dramatically enhanced STAT3 phosphorylation in HEK293T cells. However, levels of phospho-S6 kinase and phospho-ERK along with c-MYC did not change (Figure 9). These results imply that ZAP70 activity triggers activation of the JAK/STAT signaling pathway in epithelial as well as in mesenchymal cellular models of tumorigenesis. In conclusion, activated ZAP70 functions as a novel transforming factor in solid tumors and induces partial resistance to PI3K inhibition via phosphorylation and activation of STAT3 in the JAK/STAT signaling axis.
Discussion
In this study, we screened an activated tyrosine kinase library and identified ZAP70 as mediator of resistance upon simultaneous loss of p110α / p110β in MEFs (Figure 1).
ZAP70 plays critical roles in T-cell/B-cell development and T-cell receptor signaling as well as being involved in hematological malignancies (38). Our results describe a potentially novel role for ZAP70 in untransformed cell lines besides its immunological functions. Our cell viability assays showed that TEL-ZAP70 expression in MEFs significantly enhanced cellular proliferation in comparison to negative control pBABE NEO and positive control MYR-AKT1 (Figure 2). These analyses indicate that active ZAP70 can promote cellular proliferation in adherent cells.
Existing literature supports the notion that PI3K is a crucial regulator of cellular proliferation, and its impairment has been linked to severe growth defects in various cancer models (3,7). In order to abrogate PI3K signaling, we knocked out endogenous p110α/p110β in MEFs. For this, we utilized genetically engineered MEFs, where p110α and p110β can be effectively knocked out by expression of Ad-Cre. Although, active PTK6 expression was not able to restore growth in the p110α/p110β knockout system (Figure 2), TEL-ZAP70 expressing MEFs largely recovered growth upon p110α/p110β ablation. To abolish the catalytic activity of PI3K, we used clinically relevant p110α_specific Alpelisib and p110β-specific AZD6482 pharmacological inhibitors. Our results show that active ZAP70 expression can induce partial resistance to PI3K inhibition in various cellular models comprising both epithelial as well as mesenchymal cells (Figure 2 and Figure 4).
As proposed by Siveen et al., STAT3 can be activated by Src family kinases (39). Furthermore, Warmuth et al. suggested that the STAT3 signaling cascade is crucial in carcinogenesis (40). We performed western blot analysis to elucidate the mechanism of TEL-ZAP70-mediated activation of downstream pathway signaling upon PI3K ablation. Our results indicate that up-regulation of STAT3 phosphorylation on tyrosine 705 residue is a ubiquitous event in active ZAP70-expressing cell lines including RPE1-hTERT, MEFs, non-tumorigenic epithelial breast (MCF10A) and human embryonic kidney cells (HEK293T) cell lines (Figure 2, Figure 4 and Figure 9). Remarkably, p-ERK1/2 and p-S6 kinase p70 expression levels did not change in TEL-ZAP70 RPE1 cells. Additionally, dual-tyrosine phosphorylation of ERK in TEL-ZAP70 MEFs (Figure 2) and phosphorylation of S6K p70 levels in TEL-ZAP70 HEK293T (Figure 9) only slightly increased. Cha et al. suggested that the JAK/STAT3/c-MYC signaling axis is negatively regulated to modulate differentiation capability by ZAP70 in mouse embryonic stem cells (21). In contrast to their findings, our results suggested that STAT3 activation was positively correlated with ZAP70 expression (21). In addition, we found out that the c-MYC levels were not altered upon expression of active ZAP70 in TEL-ZAP70 expressing cells (Figure 4 and Figure 9). These contrasting observations might be explained by different genetic contexts in pluripotent stem cells versus differentiated or transformed mesenchymal/epithelial cell lines.
To further explore the function of ZAP70 in carcinogenesis, we investigated involvement of ZAP70 in anchorage-independent growth on soft agar in untransformed cell lines. The morphological characteristics of MEFs are different from epithelial cells as they represent cells of mesenchymal origin, whereas most carcinomas are of epithelial origin. Thus, we used non-transformed epithelial cells, MCF10A and RPE1-hTERT, to understand the impact of activated ZAP70 on initiation of carcinogenesis. We found that ZAP70-expressing MEF and RPE1 cells were capable of forming colonies on soft agar (Figure 2 and Figure 4). We also showed that ZAP70 had an additive effect in anchorage-independent growth along with H-Ras G12V (Figure 2). These findings imply that ZAP70 might have a function in tumor initiation. Although our research was mainly conducted on untransformed cell lines, it would be important to decipher the effects of ZAP70 in transformed cancer cells. Based on the literature, ZAP70 promotes cell migration and invasion of prostate cancer cell lines (20). Moreover, ZAP70 was identified as a prognostic marker in cervical squamous cell carcinoma (41), prostate adenocarcinoma (42) and colorectal cancer in response to radiation (22). Furthermore, Sun et al. indicated that ZAP70 may be a crucial regulator of metastasis in prostate cancer (42). Of note, Sadras et al. found that while ZAP70 has an oncogenic role, SYK, a member belonging to the same family, functions as a tumor suppressor in autoimmune diseases and B-cell malignancies (43).
According to our in-silico analysis, ZAP70 mRNA expression levels were found to be elevated in many solid tumor tissues in comparison to matched controls (Figure 6). Additionally, the Kaplan–Meier plots show that low ZAP70 expression in several cancer types, including gastric cancer, uveal melanoma, kidney clear cell carcinoma and low-grade glioma, was associated with better survival (Figure 7). Moreover, our survival analysis depicted that in patients with low-grade ovarian cancer who were treated with paclitaxel and cisplatin, low ZAP70 expression correlated with better survival (Figure 8). Consequently, high ZAP70 expression might have a tumorigenic/oncogenic function during the establishment of the disease and correlate with a worse prognosis.
In conclusion, our results implicate ZAP70 in tumor initiation as well as in resistance to PI3K inhibition. Although ZAP70 expression has been implicated in poor prognosis in hematological malignancies such as chronic lymphocytic leukemia (44), our results suggest that ZAP70 tyrosine kinase activation can promote cell proliferation and contribute to generation of resistance to PI3K inhibition. The induction of resistance to PI3K inhibition might be explained through activation of JAK/STAT; in particular, the JAK/STAT3 signaling axis. We believe that these findings are of strong clinical relevance as several PI3K inhibitors have already passed clinical trials and are being used in the clinic (45,46).
Conflicts of Interest
The Authors declare no competing conflicts of interest.
Authors’ Contributions
Onur Cizmecioglu: Conception and design, analysis and interpretation of data, and revision of the article. Melike Demir: Acquisition of data, analysis and interpretation of data, and drafting of the article.
Acknowledgements
This study was funded by The Scientific and Technological Research Council of Turkey (TÜBİTAK) Directorate of Science Fellowships and Grant Programs (BİDEB), project number 117C040.
References
1
Kondapalli L
,
Soltani K
&
Lacouture ME
. The promise of molecular targeted therapies: protein kinase inhibitors in the treatment of cutaneous malignancies. J Am Acad Dermatol.
53(2)
291
- 302
2005.
PMID:
16021125.
DOI:
10.1016/j.jaad.2005.02.011
3
Krause DS
&
Van Etten RA
. Tyrosine kinases as targets for cancer therapy. N Engl J Med.
353(2)
172
- 187
2005.
PMID:
16014887.
DOI:
10.1056/NEJMra044389
4
Liu P
,
Cheng H
,
Roberts TM
&
Zhao JJ
. Targeting the phosphoinositide 3-kinase pathway in cancer. Nat Rev Drug Discov.
8(8)
627
- 644
2009.
PMID:
19644473.
DOI:
10.1038/nrd2926
5
Vivanco I
&
Sawyers CL
. The phosphatidylinositol 3-Kinase AKT pathway in human cancer. Nat Rev Cancer.
2(7)
489
- 501
2002.
PMID:
12094235.
DOI:
10.1038/nrc839
6
Yang J
,
Nie J
,
Ma X
,
Wei Y
,
Peng Y
&
Wei X
. Targeting PI3K in cancer: mechanisms and advances in clinical trials. Mol Cancer.
18(1)
26
2019.
PMID:
30782187.
DOI:
10.1186/s12943-019-0954-x
8
Michmerhuizen NL
,
Leonard E
,
Kulkarni A
&
Brenner JC
. Differential compensation mechanisms define resistance to PI3K inhibitors in PIK3CA amplified HNSCC. Otorhinolaryngol Head Neck Surg.
1(2)
44
- 50
2016.
PMID:
28004037.
DOI:
10.15761/ohns.1000111
9
Fruman DA
,
Chiu H
,
Hopkins BD
,
Bagrodia S
,
Cantley LC
&
Abraham RT
. The PI3K pathway in human disease. Cell.
170(4)
605
- 635
2017.
PMID:
28802037.
DOI:
10.1016/j.cell.2017.07.029
10
Serra V
,
Scaltriti M
,
Prudkin L
,
Eichhorn PJ
,
Ibrahim YH
,
Chandarlapaty S
,
Markman B
,
Rodriguez O
,
Guzman M
,
Rodriguez S
,
Gili M
,
Russillo M
,
Parra JL
,
Singh S
,
Arribas J
,
Rosen N
&
Baselga J
. PI3K inhibition results in enhanced HER signaling and acquired ERK dependency in HER2-overexpressing breast cancer. Oncogene.
30(22)
2547
- 2557
2011.
PMID:
21278786.
DOI:
10.1038/onc.2010.626
11
Hanker AB
,
Kaklamani V
&
Arteaga CL
. Challenges for the clinical development of PI3K inhibitors: strategies to improve their impact in solid tumors. Cancer Discov.
9(4)
482
- 491
2019.
PMID:
30867161.
DOI:
10.1158/2159-8290.CD-18-1175
12
Muranen T
,
Selfors LM
,
Worster DT
,
Iwanicki MP
,
Song L
,
Morales FC
,
Gao S
,
Mills GB
&
Brugge JS
. Inhibition of PI3K/mTOR leads to adaptive resistance in matrix-attached cancer cells. Cancer Cell.
21(2)
227
- 239
2012.
PMID:
22340595.
DOI:
10.1016/j.ccr.2011.12.024
13
Nawijn MC
,
Alendar A
&
Berns A
. For better or for worse: the role of Pim oncogenes in tumorigenesis. Nat Rev Cancer.
11(1)
23
- 34
2011.
PMID:
21150935.
DOI:
10.1038/nrc2986
14
Song JH
,
Singh N
,
Luevano LA
,
Padi SKR
,
Okumura K
,
Olive V
,
Black SM
,
Warfel NA
,
Goodrich DW
&
Kraft AS
. Mechanisms behind resistance to PI3K inhibitor treatment induced by the PIM kinase. Mol Cancer Ther.
17(12)
2710
- 2721
2018.
PMID:
30190422.
DOI:
10.1158/1535-7163.MCT-18-0374
15
Le X
,
Antony R
,
Razavi P
,
Treacy DJ
,
Luo F
,
Ghandi M
,
Castel P
,
Scaltriti M
,
Baselga J
&
Garraway LA
. systematic functional characterization of resistance to PI3K inhibition in breast cancer. Cancer Discov.
6(10)
1134
- 1147
2016.
PMID:
27604488.
DOI:
10.1158/2159-8290.CD-16-0305
16
Chang M
,
Kanwar N
,
Feng E
,
Siu A
,
Liu X
,
Ma D
&
Jongstra J
. PIM kinase inhibitors downregulate STAT3(Tyr705) phosphorylation. Mol Cancer Ther.
9(9)
2478
- 2487
2010.
PMID:
20667852.
DOI:
10.1158/1535-7163.MCT-10-0321
17
Chan AC
,
Irving BA
,
Fraser JD
&
Weiss A
. The zeta chain is associated with a tyrosine kinase and upon T-cell antigen receptor stimulation associates with ZAP-70, a 70-kDa tyrosine phosphoprotein. Proc Natl Acad Sci U S A.
88(20)
9166
- 9170
1991.
PMID:
1717999.
DOI:
10.1073/pnas.88.20.9166
19
Kaur M
,
Singh M
&
Silakari O
. Insight into the therapeutic aspects of ‘Zeta-Chain Associated Protein Kinase 70 kDa’ inhibitors: a review. Cell Signal.
26(11)
2481
- 2492
2014.
PMID:
25049080.
DOI:
10.1016/j.cellsig.2014.06.017
20
Fu D
,
Liu B
,
Zang LE
&
Jiang H
. MiR-631/ZAP70: A novel axis in the migration and invasion of prostate cancer cells. Biochem Biophys Res Commun.
469(3)
345
- 351
2016.
PMID:
26620225.
DOI:
10.1016/j.bbrc.2015.11.093
21
Cha Y
,
Moon BH
,
Lee MO
,
Ahn HJ
,
Lee HJ
,
Lee KA
,
Fornace AJ Jr
,
Kim KS
,
Cha HJ
&
Park KS
. Zap70 functions to maintain stemness of mouse embryonic stem cells by negatively regulating Jak1/Stat3/c-Myc signaling. Stem Cells.
28(9)
1476
- 1486
2010.
PMID:
20641039.
DOI:
10.1002/stem.470
22
Huang MY
,
Wang JY
,
Chang HJ
,
Kuo CW
,
Tok TS
&
Lin SR
. CDC25A, VAV1, TP73, BRCA1 and ZAP70 gene overexpression correlates with radiation response in colorectal cancer. Oncol Rep.
25(5)
1297
- 1306
2011.
PMID:
21344162.
DOI:
10.3892/or.2011.1193
23
Cizmecioglu O
,
Ni J
,
Xie S
,
Zhao JJ
&
Roberts TM
. Rac1-mediated membrane raft localization of PI3K/p110β is required for its activation by GPCRs or PTEN loss. Elife.
5
e17635
2016.
PMID:
27700986.
DOI:
10.7554/eLife.17635
24
Yuzugullu H
,
Von T
,
Thorpe LM
,
Walker SR
,
Roberts TM
,
Frank DA
&
Zhao JJ
. NTRK2 activation cooperates with PTEN deficiency in T-ALL through activation of both the PI3K-AKT and JAK-STAT3 pathways. Cell Discov.
2
16030
2016.
PMID:
27672444.
DOI:
10.1038/celldisc.2016.30
25
Kuno Y
,
Abe A
,
Emi N
,
Iida M
,
Yokozawa T
,
Towatari M
,
Tanimoto M
&
Saito H
. Constitutive kinase activation of the TEL-Syk fusion gene in myelodysplastic syndrome with t(9;12)(q22;p12). Blood.
97(4)
1050
- 1055
2001.
PMID:
11159536.
DOI:
10.1182/blood.v97.4.1050
26
Győrffy B
. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J.
19
4101
- 4109
2021.
PMID:
34527184.
DOI:
10.1016/j.csbj.2021.07.014
27
Chandrashekar DS
,
Bashel B
,
Balasubramanya SAH
,
Creighton CJ
,
Ponce-Rodriguez I
,
Chakravarthi BVSK
&
Varambally S
. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia.
19(8)
649
- 658
2017.
PMID:
28732212.
DOI:
10.1016/j.neo.2017.05.002
28
Cerami E
,
Gao J
,
Dogrusoz U
,
Gross BE
,
Sumer SO
,
Aksoy BA
,
Jacobsen A
,
Byrne CJ
,
Heuer ML
,
Larsson E
,
Antipin Y
,
Reva B
,
Goldberg AP
,
Sander C
&
Schultz N
. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov.
2(5)
401
- 404
2012.
PMID:
22588877.
DOI:
10.1158/2159-8290.CD-12-0095
29
Zhao JJ
,
Gjoerup OV
,
Subramanian RR
,
Cheng Y
,
Chen W
,
Roberts TM
&
Hahn WC
. Human mammary epithelial cell transformation through the activation of phosphatidylinositol 3-kinase. Cancer Cell.
3(5)
483
- 495
2003.
PMID:
12781366.
DOI:
10.1016/s1535-6108(03)00088-6
30
Mori S
,
Chang JT
,
Andrechek ER
,
Matsumura N
,
Baba T
,
Yao G
,
Kim JW
,
Gatza M
,
Murphy S
&
Nevins JR
. Anchorage-independent cell growth signature identifies tumors with metastatic potential. Oncogene.
28(31)
2796
- 2805
2009.
PMID:
19483725.
DOI:
10.1038/onc.2009.139
31
Suzuki J
,
Sukezane T
,
Akagi T
,
Georgescu MM
,
Ohtani M
,
Inoue H
,
Jat PS
,
Goff SP
,
Hanafusa H
&
Shishido T
. Loss of c-abl facilitates anchorage-independent growth of p53- and RB- deficient primary mouse embryonic fibroblasts. Oncogene.
23(52)
8527
- 8534
2004.
PMID:
15378021.
DOI:
10.1038/sj.onc.1207894
32
Seeburg PH
,
Colby WW
,
Capon DJ
,
Goeddel DV
&
Levinson AD
. Biological properties of human c-Ha-ras1 genes mutated at codon 12. Nature.
312(5989)
71
- 75
1984.
PMID:
6092966.
DOI:
10.1038/312071a0
33
Colby WW
,
Hayflick JS
,
Clark SG
&
Levinson AD
. Biochemical characterization of polypeptides encoded by mutated human Ha-ras1 genes. Mol Cell Biol.
6(2)
730
- 734
1986.
PMID:
3537694.
DOI:
10.1128/mcb.6.2.730-734.1986
34
Bennasroune A
,
Gardin A
,
Aunis D
,
Crémel G
&
Hubert P
. Tyrosine kinase receptors as attractive targets of cancer therapy. Crit Rev Oncol Hematol.
50(1)
23
- 38
2004.
PMID:
15094157.
DOI:
10.1016/j.critrevonc.2003.08.004
35
Gocek E
,
Moulas AN
&
Studzinski GP
. Non-receptor protein tyrosine kinases signaling pathways in normal and cancer cells. Crit Rev Clin Lab Sci.
51(3)
125
- 137
2014.
PMID:
24446827.
DOI:
10.3109/10408363.2013.874403
36
Hart JR
,
Liao L
,
Yates JR 3rd
&
Vogt PK
. Essential role of Stat3 in PI3K-induced oncogenic transformation. Proc Natl Acad Sci U S A.
108(32)
13247
- 13252
2011.
PMID:
21788516.
DOI:
10.1073/pnas.1110486108
37
Bommhardt U
,
Schraven B
&
Simeoni L
. Beyond TCR signaling: emerging functions of lck in cancer and immunotherapy. Int J Mol Sci.
20(14)
3500
2019.
PMID:
31315298.
DOI:
10.3390/ijms20143500
39
Siveen KS
,
Prabhu KS
,
Achkar IW
,
Kuttikrishnan S
,
Shyam S
,
Khan AQ
,
Merhi M
,
Dermime S
&
Uddin S
. Role of non receptor tyrosine kinases in hematological malignances and its targeting by natural products. Mol Cancer.
17(1)
31
2018.
PMID:
29455667.
DOI:
10.1186/s12943-018-0788-y
40
Warmuth M
,
Damoiseaux R
,
Liu Y
,
Fabbro D
&
Gray N
. SRC family kinases: potential targets for the treatment of human cancer and leukemia. Curr Pharm Des.
9(25)
2043
- 2059
2003.
PMID:
14529415.
DOI:
10.2174/1381612033454126
41
Qin R
,
Cao L
,
Ye C
,
Wang J
&
Sun Z
. A novel prognostic prediction model based on seven immune-related RNAs for predicting overall survival of patients in early cervical squamous cell carcinoma. BMC Med Genomics.
14(1)
49
2021.
PMID:
33588862.
DOI:
10.1186/s12920-021-00885-3
42
Sun X
,
Wang L
,
Li H
,
Jin C
,
Yu Y
,
Hou L
,
Liu X
,
Yu Y
,
Yan R
&
Xue F
. Identification of microenvironment related potential biomarkers of biochemical recurrence at 3 years after prostatectomy in prostate adenocarcinoma. Aging (Albany NY).
13(12)
16024
- 16042
2021.
PMID:
34133324.
DOI:
10.18632/aging.203121
43
Sadras T
,
Martin M
,
Kume K
,
Robinson ME
,
Saravanakumar S
,
Lenz G
,
Chen Z
,
Song JY
,
Siddiqi T
,
Oksa L
,
Knapp AM
,
Cutler J
,
Cosgun KN
,
Klemm L
,
Ecker V
,
Winchester J
,
Ghergus D
,
Soulas-Sprauel P
,
Kiefer F
,
Heisterkamp N
,
Pandey A
,
Ngo V
,
Wang L
,
Jumaa H
,
Buchner M
,
Ruland J
,
Chan WC
,
Meffre E
,
Martin T
&
Müschen M
. Developmental partitioning of SYK and ZAP70 prevents autoimmunity and cancer. Mol Cell.
81(10)
2094
- 2111.e9
2021.
PMID:
33878293.
DOI:
10.1016/j.molcel.2021.03.043
44
Shvidel L
,
Bairey O
,
Tadmor T
,
Braester A
,
Ruchlemer R
,
Fineman R
,
Joffe E
,
Berrebi A
,
Polliack A
&
Israeli CLL Study Group
. Absolute lymphocyte count with extreme hyperleukocytosis does not have a prognostic impact in chronic lymphocytic leukemia. Anticancer Res.
35(5)
2861
- 2866
2015.
PMID:
25964568.
45
Juric D
,
Rodon J
,
Tabernero J
,
Janku F
,
Burris HA
,
Schellens JHM
,
Middleton MR
,
Berlin J
,
Schuler M
,
Gil-Martin M
,
Rugo HS
,
Seggewiss-Bernhardt R
,
Huang A
,
Bootle D
,
Demanse D
,
Blumenstein L
,
Coughlin C
,
Quadt C
&
Baselga J
. Phosphatidylinositol 3-kinase α-selective inhibition with alpelisib (BYL719) in PIK3CA-altered solid tumors: Results from the first-in-human study. J Clin Oncol.
36(13)
1291
- 1299
2018.
PMID:
29401002.
DOI:
10.1200/JCO.2017.72.7107
46
Engelman JA
. Targeting PI3K signalling in cancer: opportunities, challenges and limitations. Nat Rev Cancer.
9(8)
550
- 562
2009.
PMID:
19629070.
DOI:
10.1038/nrc2664