1. Introduction
Breast cancer is the most commonly diagnosed cancer in American women and is the second-leading cause of cancer-related deaths among women [1]. Breast cancers can be categorized into five separate subtypes: luminal A (ER+/PR+), luminal B (ER+/PR+/HER2 ± with high Ki-67), normal breast-like, HER2-positive (HER2-enriched), and triple-negative (ER−/PR−/HER2−; TNBC). Triple-negative and HER2-enriched subtypes are diagnosed in up to 30% of breast cancer cases and are considered more aggressive than other subtypes. As opposed to surgical resection and adjuvant chemotherapy, most patients with triple-negative and HER2-positive breast cancer now undergo systemic neoadjuvant chemotherapy as a standard of care, showing high pathological response and positive long-term outcome [2,3]. However, regardless of treatment modality, many breast cancer patients often relapse and present with recurrent breast cancer that currently remains incurable. Consequently, patients with recurrent breast cancer experience high rates of morbidity and mortality [4]. As such, there is an urgent need to identify new molecular targets that drive breast cancer growth and progression and then develop therapeutic intervention against these targets.
Janus kinase 2 (JAK2) is a non-receptor tyrosine kinase that is frequently amplified or aberrantly active in triple-negative and HER2-enriched breast cancers [5,6]. JAK2 functions as a central signaling hub that links oncogenic receptor tyrosine kinase signaling and interleukin receptor activity to Signal Transducer and Activator of Transcription 3 (STAT3) transcription factor [7]. JAK2 phosphorylates STAT3 on its Y705 residue, resulting in the dimerization and nuclear translocation of STAT3 and activation of STAT3 transcriptional activity, leading to the expression of STAT3 target genes that promote cell proliferation, differentiation, survival, and migration [8,9,10]. We have previously reported that JAK2–STAT3 can upregulate the expression of TWIST, iNOS, COX-2, and STAT1 [11,12,13,14]. Since JAK2–STAT3 signaling is frequently upregulated in aggressive breast cancers and the upregulation is correlated to poor overall clinical outcome, the pathway has been regarded as an important target of cancer drug development. Of note, three orally active JAK2 inhibitors, ruxolitinib, baricitinib, and fedratinib, have been approved by the FDA for rheumatoid arthritis and myelofibrosis [15,16,17]; they are being evaluated for human cancers in numerous clinical trials.
Tropomyosin receptor tyrosine kinase A (TrkA) is a neurotrophic receptor encoded by the NTRK1 gene and is implicated in numerous cancer types, in part, due to its propensity for forming oncogenic NTRK fusion genes that drive the malignant progression of some cancer types [18]. However, wild-type TrkA has also been shown to function as a potent oncogenic driver when overexpressed in cells. The binding of its cognate ligand β-Nerve Growth Factor (β-NGF) results in TrkA receptor dimerization and activation, allowing for TrkA to phosphorylate its substrates and activate pathways, such as RAS/MAPK and PI3 kinase, resulting in increased cell proliferation and survival [19]. Similar to the JAK2–STAT3 signaling, the TrkA pathway is implicated in breast cancer progression. The overexpression of TrkA confers increased cell proliferation, migration, and metastatic capacity of breast cancer cells along with concomitant reduction in apoptosis [20,21,22,23]. TrkA-overexpressing breast tumors were found to grow significantly faster in vivo with increased development of lung, liver, and brain metastases, suggesting that TrkA activation contributes to breast cancer metastasis [20,21,24]. Accordingly, TrkA, as well as other members of the TRK protein family, have become therapeutic targets for inhibitor design [18]. Two orally active TrkA inhibitors, Larotrectinib and Entrectinib, have been approved by the FDA for treating solid tumors bearing NTRK fusion genes [25,26].
Although previous reports show that the activation of TrkA by its cognate ligand, β-NGF, promotes the phosphorylation of STAT3 on its Y705 and S727 residues and promotes STAT3 transcriptional activity, the co-activation of JAK2–STAT3 and TrkA signaling pathways, and the mechanisms underlying their crosstalk, have never been investigated in any cell or tumor type [27,28]. In this study, we report for the first time that signaling crosstalk between the JAK2–STAT3 and TrkA pathways is mediated through TrkA–STAT3 protein–protein interactions and TrkA-mediated phosphorylation of STAT3 Y705, leading to STAT3 activation and increased expression of STAT3 target genes. Immunohistochemical (IHC) analyses of human breast cancer samples reveal a high expression of p-STAT3 (Y705) and p-TrkA (Y490) in the tissue, supporting the co-activation of JAK2–STAT3 and TrkA pathways. Survival analyses of publicly available datasets reveal dismal survival and metastasis-free survival in patients with co-activation of JAK2–STAT3 and TrkA. Finally, we show that TrkA interacts directly with and phosphorylates STAT3 on Y705 to promote oncogenic gene transcription, implicating a novel mechanism to activate STAT3-mediated breast cancer metastasis in TrkA-overexpressing tumors.
2. Results
2.1. JAK2–STAT3 and TrkA Pathways Are Significantly Co-Activated in Triple-Negative and HER2-Enriched Breast Cancers
To determine the extent to which JAK2–STAT3 and TrkA signaling pathways are co-activated in breast tumor samples, we performed IHC staining on 33 node-positive breast carcinomas across three major subtypes (triple-negative, HER2-enriched, and luminal subtypes) to detect phosphorylated TrkA (Y490) and STAT3 (Y705) proteins. Results showed that node-positive triple-negative breast cancer (TNBC) and HER2-positive breast tumors highly co-express p-TrkA and p-STAT3Of the nine TNBC samples, we found that 100% of them have high p-STAT3 and 77.8% of these samples concurrently express high p-TrkA. Similarly, 83% of HER2-positive breast tumors were highly stained for both p-TrkA and p-STAT3. Although over 80% of node-positive luminal samples showed high p-STAT3 levels, only 41% of these samples also highly expressed p-TrkA. These results indicated that concurrently activated JAK2–STAT3 and TrkA signaling is preferentially found in TNBC and HER2-positive breast cancer. We found a 61% (20/33) positivity rate for p-TrkA and p-STAT3 co-overexpression after examining each individual sample. Representative IHC images are shown in Figure 1A-bottom. As expected, p-TrkA was primarily detected on cell membrane (arrows), while p-STAT3 was found primarily within the nucleus (arrows). Since STAT3 phosphorylation by receptor tyrosine kinases and non-receptor tyrosine kinases occur in the cytoplasm, diffuse cytoplasmic staining was also observed. In a negative control experiment for IHC (no primary antibody, secondary antibody only), we did not observe any positive staining from a mammary fat pad xenograft of BT474-TtzmR, which is a breast cancer cell line that is resistant to trastuzumab. Then, we confirmed our IHC findings using immunofluorescent staining of a primary patient breast tumor and found that cells with p-TrkA were also positive for p-STAT3 (Figure 1B), supporting our hypothesis that these two proteins are co-expressed and co-localized in breast tumor tissues. Furthermore, we expanded our investigation of TrkA and STAT3 co-activation by performing a Western blot analysis using breast cancer cell lines of varying subtypes (Figure 1C). Taken together, our findings suggest that the co-activation of JAK2–STAT3 and TrkA pathways may preferentially occur in TNBC and HER2-enriched breast cancers.
To validate the above-mentioned results, we analyzed publicly available breast cancer databases to determine the extent to which target genes downstream of the JAK2–STAT3 and TrkA pathways are enriched in triple-negative and HER2-enriched breast cancers. To this end, we applied Gene Set Enrichment Analyses (GSEA) using a STAT3 activation signature, TrkA activation signature, or a combination of STAT3 and TrkA signatures to examine 1533 breast cancer patient samples from TCGA and GEO databases [29,30,31,32,33,34,35]. Here, we observed that TNBC and HER2-enriched breast cancer showed significantly higher enrichment of the STAT3 activation signature, but not the TrkA activation signature, when compared to luminal samples (Figure 1D,E). Importantly, the combined STAT3 and TrkA activation signature was significantly enriched in triple-negative and HER2-enriched breast cancers in comparison to luminal subtypes (Supplementary Figure S1A,B). Further analysis of the GEO datasets revealed significantly higher STAT3 activation scores in both triple-negative and HER2-enriched tumors when compared to luminal tumors (Figure 1F); TrkA activation is not significantly different between subtypes (Figure 1G). However, combined STAT3–TrkA activation signature is highest in TNBC (Figure 1H). Since luminal breast cancers can be divided into A (LumA) and B (LumB) subtypes while triple-negative breast cancers can be divided into claudin-low (CL; mesenchymal) and basal subtypes, we further analyzed the GEO dataset and found similar results (Supplementary Figure S1C–E). To determine whether our findings are relevant in other cancers, we performed GSEA using 592 colorectal cancer patient data as they frequently express oncogenic NTRK1 fusions [36] and found significant enrichment of STAT3 activation gene signature in patients with high NTRK1 mRNA levels (Supplementary Figure S2), suggesting that the signaling crosstalk may also be found in other malignant tissues. Taken together, the results in Figure 1 and Figures S1 and S2 demonstrated, for the first time, that JAK2–STAT3 and TrkA pathways are frequently co-activated in triple-negative and HER2-enriched breast cancers.
2.2. STAT3 and TrkA Proteins Directly Interact in Triple-Negative and HER2-Enriched Breast Cancers
The observation of co-overexpression of p-STAT3 and p-TrkA prompted us to examine if the co-expression leads to their physical interaction. Of note, these two proteins have never been reported to physically interact. For this, we conducted immunoprecipitation (IP) followed by Western blot using HEK293 cells transfected with flag-tagged STAT3 (STAT3-WT-Flag) and found an interaction between STAT3 and p-TrkA/total TrkA (Figure 2A). The interaction was confirmed in reciprocal IP using a p-TrkA antibody (Figure 2B). To determine whether the STAT3–TrkA interaction occurs in breast cancer cells, we immunoprecipitated endogenous p-TrkA (Y490) from MDA-MB-468 TNBC cells, and Western blot results confirm that endogenous STAT3 co-immunoprecipitates with p-TrkA (Figure 2C). A cell-free TrkA kinase assay followed by immunoprecipitation of p-TrkA reveals that recombinant STAT3 co-immunoprecipitates with p-TrkA, suggesting that these two proteins directly interact (Figure 2D). TrkA undergoes oncogenic fusions in certain cancer types but rarely in breast cancer [18]. According to our datamining of 8767 breast cancer samples (using cBioPortal), only one sample was found to express the TrkA fusion (1/8767). Nevertheless, we next determined whether TrkA fusions also interact with STAT3. To address this, we performed a cell-free kinase assay using GST-tagged Tropomyosin 3 (TPM3)–TrkA fusion protein and recombinant human STAT3, which is followed by immunoprecipitation using an anti-GST antibody. Western blot analysis reveals that STAT3 does not immunoprecipitate with TPM3–TrkA fusion protein (Figure 2E). Next, we asked whether the interaction between STAT3 and p-TrkA is dependent on the phosphorylation of STAT3 on its Y705 residue. To address this question, we compared constitutively active STAT3 (CA) with non-phosphorylation STAT3 mutant (Y705F) for their interaction with p-TrkA using IP-Western blot. The result showed that the STAT3-Y705F mutant lost the ability to interact with p-TrkA, indicating that the Y705 residue is required for the STAT3–TrkA interaction (Figure 2F). Interestingly, these data suggested that TrkA without Y490 phosphorylation can weakly interact with the STAT3-Y705F mutant (Figure 2F). To confirm whether TrkA kinase activity is critical for its interaction with STAT3, we immunoprecipitated flag-tagged STAT3 in HEK293 cells co-transfected with either TrkA-WT or TrkA-K538N (kinase-dead) [37]. Western blot analysis reveals that both TrkA-WT and TrkA-K538N (kinase-dead) mutant co-immunoprecipitate with STAT3, indicating that TrkA kinase activity is not critical for its interaction with STAT3 (Figure 2G).
To confirm the novel STAT3–TrkA interaction and the requirement of STAT3-Y705 for the interaction, we conducted immunofluorescent (IF) staining–confocal microscopy on HEK293 cells co-transfected with TrkA and either STAT3-WT-Flag or STAT3-Y705F-Flag plasmids. To induce phosphorylation of both STAT3 and TrkA, transfected cells were stimulated with human epidermal growth factor (hEGF) and human β-nerve growth factor (β-NGF), respectively. Results of these experiments (Figure 2H, top row) showed that p-TrkA co-localized with STAT3-WT-Flag on or close to the inner cell membrane, as indicated by the yellow merged signals (arrows) from p-TrkA (green) and STAT3-Flag (red). In contrast, p-TrkA failed to co-localize with STAT3-Y705F-Flag, which supports the IP–Western blot finding (Figure 2H, bottom row). Furthermore, we asked whether the STAT3–TrkA interaction can be detected in a triple-negative human breast cancer xenograft. For this, we performed IF staining and confocal microscopy to detect p-STAT3 and p-TrkA in in MDA-MB-231 mammary fat pad (MFP) tumors, and the result showed significant co-localization of both proteins as indicated by the yellow merged signal (arrows) (Figure 2I). Taken together, results in Figure 2 uncovered a novel physical interaction between STAT3 and TrkA proteins, and that their interaction is dependent on STAT3-Y705 phosphorylation.
2.3. TrkA Phosphorylates STAT3 at Y705
Since our data showed that phosphorylated TrkA can interact with STAT3 and the interaction is abolished when STAT3-Y705 is mutated, we hypothesized that TrkA phosphorylates STAT3 at its Y705 residue. While TrkA activation has been observed to promote the phosphorylation of STAT3 on Y705 and S727 residues [27,28] and the pharmacological inhibition of TrkA leads to reduction of STAT3 Y705 phosphorylation [38,39], TrkA phosphorylation of STAT3-Y705 has never been reported. To test this hypothesis, we first conducted a cell-free TrkA kinase assay using recombinant (rec) human TrkA and STAT3 proteins. Kinase reactions were analyzed using Western blot; results showed that TrkA phosphorylates STAT3 at Y705 (Figure 3A). To validate the results of cell-free kinase assay, we immunoprecipitated STAT3 and STAT3-Y705F from transfected HEK293 cells and subjected the immunoprecipitates to the cell-free TrkA kinase assay and Western blot. We found that immunoprecipitated STAT3, but not STAT3-Y705F, can be phosphorylated by TrkA (Figure 3B). We did not detect any signal from an antibody against p-STAT3 (S727), confirming that TrkA kinase specifically targets STAT3-Y705 for phosphorylation. To confirm that the induction of STAT3-Y705 phosphorylation is indeed due to TrkA kinase activity, we treated MDA-MB-231 TNBC cells with TrkA inhibitor Entrectinib. Our Western blot analysis shows significant reduction of p-TrkA (Y490) and concomitant reduction of p-STAT3 (Y705), supporting our hypothesis that TrkA kinase activity is required for STAT3 phosphorylation at Y705 (Figure 3C). To further validate TrkA phosphorylation of STAT3 at Y705, we used liquid chromatography/mass spectrometry (LC/MS) analysis to identify phosphopeptides within the recombinant human STAT3 protein that have been phosphorylated by recombinant TrkA kinase; the results showed the highest enrichment of p-STAT3 (Y705) peptide (Figure 3D; arrow). To determine if TrkA ligand, β-NGF, can induce STAT3 phosphorylation and nuclear transport, we conducted IF followed by confocal microscopy and we showed (Figure 3E) that β-NGF induced TrkA activation, as indicated by p-TrkA signals (green), and STAT3 activation as indicated by nuclear p-STAT3 (red). Together, these results show, for the first time, that TrkA phosphorylates STAT3 and STAT3-Y705 is targeted by TrkA.
2.4. TrkA Phosphorylation Promotes STAT3 Nuclear Import and Transcriptional Activity and Breast Cancer Stem Cells
STAT3 phosphorylation at Y705 leads to STAT3 nuclear translocation and activation of transcriptional activity. Previous studies have shown that STAT3-mediated transcriptional activity can be enhanced by the addition of TrkA ligand β-NGF [28] or overexpression of Trk fusion genes [27]; however, whether it is occurred directly through TrkA-mediated STAT3 activation has remained unclear. Here, we tested whether TrkA-mediated phosphorylation of STAT3 at Y705 activates the ability of STAT3 to under nuclear transport and activate a STAT3-targeted promoter. For this, we transfected HER2-enriched SKBR3 cells (with low endogenous TrkA and STAT3) with TrkA-WT or control vector plus pGAS-Luc (a luciferase reporter under the control of STAT3 response elements) and determined luciferase activity. As shown in Figure 4A, β-NGF stimulation of untransfected SKBR3 breast cancer cells following starvation induced significant upregulation of pGAS-Luc activity, suggesting that the activation of endogenous TrkA can induce STAT3 nuclear translocation and transcriptional activity. Accordingly, TrkA overexpression strongly activated pGAS-Luc activity indicating TrkA activates STAT3 transcriptional activity, which is in agreement with previous reports (Figure 4B) [27]. To determine whether TrkA ligand β-NGF induces pGAS-Luc and if TrkA and STAT3-WT co-induces pGAS-Luc activity, we co-transfected SKBR3 cells with STAT3-WT and pGAS-Luc plasmids (with or without the TrkA plasmid), serum-starved the cells, and stimulated cells with or without β-NGF, and determined luciferase activity. As shown in Figure 4C, β-NGF induced pGAS-Luc activity in the presence of STAT3 or TrkA-STAT3. Importantly, cells with TrkA-STAT3 co-transfection and β-NGF stimulation had the highest pGAS-Luc activity.
To determine if β-NGF-induced activation of STAT3 is attributed to TrkA kinase activity, we examined TrkA-WT versus TrkA-K538N (dominant-negative kinase-dead mutant) [37] in their ability to induce STAT3 Y705 phosphorylation and nuclear transport. HEK293 cells were co-transfected with STAT3-WT-Flag and either TrkA–WT (Figure 4D) or TrkA-K538N (Figure 4E), starved for serum for 16 h, stimulated with β-NGF, and examined for STAT3 Y705 phosphorylation and nuclear import using IF staining-confocal microscopy. As expected, neither p-TrkA nor p-STAT3 was detected in serum-starved cells. Stimulation with β-NGF activated TrkA, as indicated by the detection of p-TrkA (green signals) at the cell membrane (Figure 4D). Importantly, β-NGF induced p-STAT3, which was detected in the cytosol as well as in the nucleus, indicating that TrkA activation induced STAT3 Y705 phosphorylation and nuclear translocation (yellow merge signals, arrows; Figure 4D). In contrast, no phosphorylation events, in either TrkA or STAT3, were detected in cells overexpressing the kinase-dead TrkA-K538N mutant (Figure 4E).
Activated STAT3, a transcription factor, can induce the transcription of many target genes implicated in malignant progression, some of which are involved in cell proliferation, survival, invasion, and metastasis [40]. To further determine if STAT3 transcriptional activity can be activated through TrkA, we performed RT-PCR analysis of representative STAT3 target genes using SKBR3 cells overexpressing TrkA-WT or TrkA-K538N. Results showed that TrkA WT-overexpressing cells had significantly higher mRNA levels of SOX2 (Figure 4F) and MYC (Figure 4G), which are two known STAT3 target genes frequently associated with breast cancer cell stemness and aggressiveness [41,42]. Interestingly, cells overexpressing TrkA-WT also showed a modest increase in STAT3 mRNA expression, which is abrogated in TrkA-K538N expressing cells (Supplementary Figure S13A,B), suggesting that TrkA can regulate STAT3 expression as well as STAT3 transcriptional activity. Conversely, SOX2 or MYC mRNAs were significantly downregulated in cells expressing TrkA-K538N, nearly to the levels found in the vector controls, suggesting that TrkA kinase activity is essential for inducing STAT3 target gene transcription in breast cancer cells. To confirm that the increase in SOX2 and MYC gene expression is indeed due to TrkA kinase activity, we treated MDA-MB-468 TNBC cells, which have high levels of p-TrkA (Y490) (Figure 1C), with a TrkA inhibitor Entrectinib and performed RT-PCR for SOX2 and MYC. We found that while SOX2 expression was not altered upon TrkA inhibition (Figure 4H), MYC expression is significantly downregulated in the absence of TrkA activity (Figure 4I). It is important to note that treatment with Entrectinib does not impact NTRK1 mRNA (Supplementary Figure S13C); however, there is a significant reduction in STAT3 expression upon treatment with Entrectinib (Supplementary Figure S13D), which is in agreement with our results that TrkA kinase activity may also regulate STAT3 expression in addition to STAT3 activity (Supplementary Figure S13B). To confirm whether TrkA-mediated STAT3 activation can indeed promote breast cancer cell stemness, we transiently overexpressed TrkA in BT20 cells, which have relatively lower levels of activated TrkA and STAT3, with and without constitutively active STAT3 (STAT3-CA), and we performed flow cytometry to determine changes in the CD44high/CD24low population of cells. Relative to vector controls, TrkA-overexpressing BT20 cells had significantly higher percentage of CD44high/CD24low cells (4% vs. 10.3%, respectively) (Figure 4J). Compared to TrkA alone, the addition of STAT3-CA further enhanced the effect of TrkA in promoting the CD44high/CD24low population (10.3% vs. 23.3%), which supports the signaling crosstalk between JAK2/STAT3 and TrkA pathways (Figure 4K). Then, we performed mammosphere assays on BT20 cells transiently overexpressing TrkA, with and without STAT3-CA. Compared to vector controls, TrkA-overexpressing BT20 cells formed significantly higher number of mammospheres (Figure 4L). The addition of STAT3-CA enhanced the effect of TrkA on the mammosphere-forming ability of BT20 cells and yielded the highest number of mammospheres, which is reflected in the representative images (Figure 4M). To further validate the induction of cancer cell stemness in breast cancer cells overexpressing TrkA, we performed an aldehyde dehydrogenase (ALDH) assay in BT20 TNBC cells after transient transfection with TrkA, with and without STAT3-CA. The induction of ALDH activity is a marker for cancer cell stemness and may be predictive of poor clinical outcome in breast cancer [43]. Our ALDH assay results indicated that TrkA can increase ALDH activity in BT20 breast cancer cells when compared to vector controls (Figure 4N), although this did not reach statistical significance. However, in BT20 cells overexpressing TrkA, STAT3-CA significantly increased ALDH activity (Figure 4O). Taken together, these data suggest that TrkA-mediated phosphorylation and activation of STAT3 can induce STAT3 nuclear transport and the transcription of STAT3 target genes that promote breast cancer stem cells.
2.5. Co-Activation of JAK2–STAT3 and TrkA Pathways Is Correlated with Poor Overall and Bone Metastasis-Free Survival of Triple-Negative and HER2-Enriched Breast Cancers
We have shown that activated STAT3 and activated TrkA can be found in triple-negative and HER2-enriched breast cancers and that the combined JAK2–STAT3/TrkA activation signatures are significantly enriched in these two subtypes (Figure 1). Since triple-negative and HER2-positive breast cancers are more aggressive than the luminal subtype and have higher rates of metastasis, we examined GEO datasets to determine if JAK2–STAT3 and TrkA pathway activation, individually and jointly, is associated with overall metastasis-free survival in patients with these two tumor subtypes. In 166 TNBC patients, we found tumors with high STAT3 activation, high TrkA activation, and co-activation to correlate with a shortened time to develop overall metastasis (Figure 5A). In 96 HER2-enriched breast cancer patients, we observed that high STAT3 activation and high JAK2-STAT3/TrkA co-activation, but not high TrkA activation, correlate with a shortened time to develop overall metastasis (Figure 5B). Furthermore, in TNBC, we found that high TrkA activation and high co-activation, but not high STAT3 activation, were associated with a shortened time to develop bone metastasis (Figure 5C). In HER2-enriched breast cancer, high STAT3 activation and high co-activation, but not high TrkA activation, correlated with an increase potential to develop bone metastasis (Figure 5D). Together, these data indicate the novel prognostic value for the co-activation of JAK2–STAT3 and TrkA for metastatic triple-negative and HER2-enriched breast cancers.
2.6. Co-Activation of JAK2-STAT3 and TrkA Pathways Is Associated with a Shortened Time to Develop Brain and Lung Metastasis of Triple-Negative and HER2-Enriched Breast Cancers
Further survival analyses of triple-negative breast cancer datasets showed that high TrkA activation and high co-activation, but not STAT3 activation, were associated with shortened time to develop brain metastasis (Figure 6A). For HER2-enriched breast cancer (N = 96), we found that high STAT3 activation and high co-activation, but not TrkA activation, were associated with a shortened time to develop brain metastasis (Figure 6B). For lung metastasis-free survival of triple-negative breast cancer (N = 166), we observed that high STAT3 activation, high TrkA activation, and high co-activation, were associated with a high likelihood to develop lung metastasis (Figure 6C). For lung metastasis-free survival of HER2-enriched breast cancer, high STAT3 activation and high co-activation, but not TrkA activation, were associated a higher potential to develop lung metastasis (Figure 6D). Collectively, the results in Figure 6 indicate the novel prognostic value for the co-activation of JAK2–STAT3 and TrkA for brain- and lung-metastatic triple-negative and HER2-enriched breast cancers.
3. Discussion
In this study, we have made the following novel and impactful findings. (1) JAK2–STAT3 and TrkA pathways are significantly co-activated in triple-negative and HER2-positive human breast cancer specimens. (2) STAT3 and TrkA proteins interact directly as observed in cells and MFP xenograft tissues. (3) The STAT3–TrkA interaction is dependent on STAT3 Y705. (4) TrkA phosphorylates STAT3 at Y705 to induce its activation, nuclear import, and transcriptional activation of known STAT3 target genes, SOX2 and MYC. (5) Co-activation of TrkA and STAT3 promotes breast cancer stem cells. (6) JAK2–STAT3 and TrkA pathway co-activation is significantly higher in triple-negative and HER2-enriched breast cancers when compared to luminal subtypes of cancers. (7) High JAK2–STAT3 and TrkA co-activation in TNBC and HER2-enriched breast cancers is associated with shorter metastasis-free survival and significantly higher risks for lung, bone, and brain metastasis. Through these significant observations, our study uncovers TrkA-mediated STAT3 phosphorylation as a novel mechanism underlying the crosstalk between the JAK2–STAT3 and TrkA pathways in triple-negative and HER2-enriched breast cancers, and it identifies the novel prognostic value for the pathway co-activation.
While previous reports showed TrkA activation can lead to STAT3 phosphorylation [27,28] and that the inhibition of TrkA can abrogate STAT3 phosphorylation [38,39], the current study is the first to show that TrkA directly and physically interacts with STAT3 and that TrkA is a novel upstream kinase of STAT3. TrkA-mediated phosphorylation and activation of STAT3 represents a novel mechanism by which STAT3 is activated in human cancers. While JAK2 has been regarded as the key tyrosine kinase responsible for phosphorylating STAT3 at Y705, it has been reported that loss of JAK2 does not significantly affect STAT3 activation [44], suggesting that STAT3 phosphorylation and activation can be sustained independently of JAK2 activity. It has also been shown that JAK1 plays an important role in STAT3 activation in breast cancer [45]. Since all three kinases target Y705 of STAT3, it would be important to determine and compare the ability of these three STAT3-activating kinases in binding to STAT3 and phosphorylating STAT3.
Since STAT3 serves as a converging point of several oncogenic signaling pathways, including JAK1/2, EGFR, HER2, Src [40], and TrkA (current study) and STAT3 transcriptionally activates expression of many important cancer-related genes [40], it is critically important to completely inhibit STAT3 activity to effectively inhibit cancer growth and metastasis. In this regard, combined inhibition of multiple STAT3-activating kinases should be exploited. To facilitate the testing of these combination treatments, it is an important task to determine the extent to which these STAT3-activating kinases are co-activated in breast cancers. It is also imperative to examine if there is any cooperation among these kinases.
Given the observed co-activation of JAK2 and TrkA in triple-negative and HER2-enriched breast cancers, we speculate that co-targeting of JAK2 and TrkA would be a promising new therapeutic modality against breast cancer. Currently, there are two FDA-approved orally active TrkA inhibitors, entrectinib and larotrectinib. Entrectinib (Rozlytrek®) was approved in August 2019 for treatment of solid tumors that bore oncogenic NTRK1 fusions. Larotrectinib (Vitravki®) was approved in 2018 as a tissue-agnostic inhibitor for any solid tumors bearing NTRK1 fusions [25,26]. There are three FDA-approved JAK2 inhibitors, ruxolitinib, baricitinib, and fedratinib [15,16,17]. Ruxolitinib (also known as Jakafi®) was approved in 2011 for treatment of myelofibrosis and is currently in Phase 1 and 2 clinical trials for treating breast cancer patients. Baricitinib (Olumiant®) is approved for the treatment of rheumatoid arthritis and is not being evaluated for breast cancer. Fedratinib (Inrebic®) is approved for treating primary and secondary myelofibrosis; however, it has not entered clinical trials for breast cancer patients. Additional JAK1/2 inhibitors, pacritinib, gandotinib, lestaurtinib, and momelotinib, are in clinical trials for myelofibrosis, myeloid neoplasms, or relapsed lymphomas, but they have not yet been investigated clinically in breast cancer patients [46,47,48,49].
We also observed an increase in the transcription of STAT3 target genes, SOX2 and MYC, upon the overexpression of TrkA (Figure 4F,G). SOX2 is critical for maintaining the pluripotency of embryonic stem cells and has been credited with the induction and maintenance of breast cancer cell stemness [50,51]. The induction of breast cancer cell stemness is frequently implicated in breast cancer recurrence as well as therapeutic resistance [52]. Indeed, SOX2 expression is adequate to induce mammosphere formation and promotes tumor growth in a xenograft model of breast cancer [50]. In addition to promoting cell proliferation, MYC overexpression has also been implicated in breast cancer cell stemness. Yin et al. and other groups show that the overexpression of MYC in triple-negative breast cancer cells can induce the expression of ALDH1 protein, which is considered a stem cell marker in a number of solid tumors, resulting in higher rates of stem cell self-renewal [53,54]. In agreement with these findings, we show that TrkA-mediated induction breast cancer cell stemness and mammosphere-forming ability can be enhanced by STAT3, further supporting a cooperative signaling crosstalk between the TrkA and JAK2/STAT3 pathways (Figure 4L,M). Thus, TrkA-mediated activation of STAT3 can induce the expression of stemness genes to support breast cancer progression. Whether the downregulation of both JAK2–STAT3 and TrkA signaling can suppress breast cancer cell stemness remains to be examined.
In summary, our study reports that TrkA is a novel upstream activating kinase for STAT3, the TrkA and JAK2 pathways are concurrently activated in triple-negative and HER2-positive breast cancers, and the co-activation serves as a contributing factor to breast cancer stem cells and a new prognosis indicator for breast cancer metastasis to various organs. Together, our findings shed new light into the pathobiology of breast cancer with concurrent expression of TrkA and STAT3. Our observations also provide first evidence to support future testing of a novel combination therapy for metastatic triple-negative and HER2-enriched breast cancers through simultaneous inhibition of JAK2 and TrkA kinase activities.
4. Materials and Methods
4.1. Cell Culture, Cell Lines, and Reagents
All human cell lines were purchased from the American Type Culture Collection (Manassas, VA, USA). The breast cancer cell lines used in this study included triple-negative MDA-MB-231, MDB-MB-468, and BT20 cells, as well as LumB HER2-amplified SKBR3 cells. Cells were maintained in Corning DMEM (high glucose with sodium pyruvate; Corning Incorporated, Corning, NY, USA) or MEM (ThermoFisher Scientific, Waltham, MA, USA) supplemented with 10% heat-inactivated fetal bovine serum (Gibco, Gaithersburg, MD, USA) and 1% penicillin/streptomycin.
4.2. Immunohistochemistry (IHC) and Tumor Scoring
Breast cancer tissue microarray (BR10010e) was purchased from US Biomax (Rockville, MD, USA) and stained for p-TrkA (Y490) and p-STAT3 (Y705) as previously described [55,56]. IHC scoring was performed as previously described [55]. Briefly, histologic scores (H-scores) were computed from both percent positivity (A%, A = 1–100) and intensity (B = 0–3) using the following equation: H-score = A × B. For p-TrkA (Y490) and p-STAT3 (Y705) H-score, tumors were classified as low (H-score = 0–150) or high (H-score = 150–300).
4.3. Plasmids and Transfection
pCMV5-TrkA was generated by Raymond Birge (Addgene #15002, Watertown, MA, USA) [57]. pLEGFP-Y705F-STAT3 was made by George Stark (Addgene #71445) [58]. pLEGFP-STAT3-WT-Flag was generated by performing site-directed mutagenesis on the Y705F mutant plasmid to restore Y705 using Quik-Change II XL mutagenesis kit (Agilent Technologies, Santa Clara, CA, USA). Constitutively activate STAT3 (STAT3-CA) was generously gifted by Keping Xie at the University of Texas-MD Anderson Cancer Center [59]. pCMV5-TrkA-K538N kinase-dead mutant was generated by performing site-directed mutagenesis on the wild-type plasmid. Mutagenesis primers are described in Supplementary Table S1. All transfections were performed using XtremeGene HP transfection reagent (Roche, Basel, Switzerland) and OPTI-MEM I Reduced Serum Medium (Gibco) as per manufacturers’ protocols.
4.4. Metastasis-Free Survival, Gene Set Enrichment Analyses (GSEA), and Statistical Analysis
Publicly available breast tumor RNA-Seq expression profiles were retrieved from TCGA and GEO (GSE 2034, 2603, 5327, 12,276) with subtype and/or metastasis-free survival information from 166 triple-negative breast cancer patients and 96 HER2-enriched breast cancer patients. Median centering was used to generate STAT3 or TrkA activation signature comprised of known genes regulated by STAT3 [35] or TrkA (R-HSA-187037) [34] pathway activity. For Kaplan–Meier survival analyses, STAT3 or TrkA activation scores were calculated, and patients were stratified 50:50 (high vs. low activation) based on the median activation score and plotted using GraphPad Prism 8. The Log-rank test was used to determine significance. GSEA was performed by generating the Gene MatriX file (.gmx) using published gene signatures for STAT3 and/or TrkA activation [29,60,61]. Gene Cluster Text files (.gct) and Categorical Class files (.cls) were generated based on STAT3 and/or TrkA score in the GEO datasets. The number of permutations for GSEA was set to 1000 and the GEO gene list was used as the chip platform. TrkA fusion gene (TFG) mRNA was analyzed from TCGA breast cancer datasets. Results are represented as mean ± SD. Log-rank test and Tukey’s test were performed using GraphPad Prism 8 (GraphPad Software, San Diego, CA, USA).
4.5. Immunoprecipitation (IP) and Western Blot
For crude protein lysates, adherent cells were washed with cold 1× PBS before lysis with RIPA buffer (ThermoFisher Scientific) or SDS-free RIPA buffer. Lysis buffers were prepared with 1× Halt Protease/Phosphatase inhibitor (ThermoFisher Scientific). Lysates were sonicated and clarified by centrifugation before quantification with Bradford assay (Bio-Rad, Hercules, CA, USA). Flag IP experiments were performed as previously described [62]. To IP p-TrkA (Y490), anti-p-TrkA (Y490) antibody was conjugated to Dynabeads Protein G beads (Invitrogen, Carlsbad, CA, USA) as per manufacturer’s protocol before overnight incubation with crude protein lysates at 4 °C.
Western blot of IP products was performed as previously described [55]. All antibodies are described in Supplementary Table S2.
4.6. Immunofluorescence (IF) Staining and Confocal Microscopy
HEK293 cells were seeded at 1.5 × 104 cells/well into Nunc Lab-Tek II 8-well chamber slide (ThermoFisher) and transfected with the indicated plasmids. Transfected cells were starved overnight before stimulation with 100 ng/mL of β-NGF (R&D Systems, Minneapolis, MN, USA) or human EGF (Sigma Aldrich, St. Louis, MO, USA) for 20 min at 37 °C. Cells were fixed with 4% paraformaldehyde for 30 min at ambient temperature. Fixed cells were washed in 1× PBS before blocking (5% goat serum + 0.3% Triton X-100 in 1× PBS) for 1 h, followed by incubation with primary antibody overnight at 4 °C. Cells were incubated with a fluorophore-conjugated secondary antibody (see Supplementary Table S2) for 1 h at room temperature. Nuclei were stained with DAPI while mounting slides with a coverslip (Vector Labs, Burlingame, CA, USA). Slides were imaged using the Olympus FV1200 confocal microscope (Center Valley, PA, USA) or ImageXpress Micro Confocal System (Molecular Devices, San Jose, CA, USA).
For MFP tissue immunostaining, 5 µm frozen sections were used. All sections were fixed using cold acetone for 10 min. and blocked with 1× PBS/0.3% Tween/5% goat serum. Sections were stained using mouse anti-phospho-STAT3 (Y705) and rabbit anti-phospho-TrkA (Y490) simultaneously, in 1× PBS/0.3% Tween/5% goat serum overnight at 4 °C. Phospho-STAT3 (Y705) was detected with goat anti-mouse Alexa Fluor 488 (Invitrogen) while phospho-TrkA (Y490) was detected using goat anti-rabbit Alexa Fluor 568 (Invitrogen). Nuclei were labeled with DAPI, and slides were mounted with a cover slip (Vector Labs). Slides were imaged using the Olympus FV1200 confocal microscope, and Z-stack images were acquired and processed with FV10-ASW 4.2 viewer software.
4.7. TrkA Cell-Free Kinase Assay, Mass Spectrometry, and IP-Kinase Assay
Recombinant human STAT3 (Creative Biomart, Shirley, NY, USA) was incubated with enzymatically active recombinant human TrkA (Abcam, Cambridge, UK) or TPM3-TrkA (MyBioSource, San Diego, CA, USA) for 30 min at 30 °C in kinase reaction buffer (52 mM Tris, 43 mM MgCl2, 1 mg/mL BSA, pH 7.5) in the presence of 50 mM ATP (New England BioLabs, Ipswich, MA, USA). Kinase reactions were terminated by the addition of Laemmli sodium dodecyl sulfate (SDS) buffer and resolved on 8% SDS polyacrylamide gel. For LC/MS analyses, STAT3 bands were excised using sterile blades and submitted for liquid chromatography/mass spectrometry at the Massachusetts Institute of Technology Proteomics Core. For Coomassie staining, gels were fixed (50:40:10/methanol:water:acetic acid) for 1 h, followed by overnight staining with Coomassie R-250 (G-Biosciences, St. Louis, MO, USA) and de-staining until bands of interest were visible. Coomassie-stained gels were imaged using the Bio-Rad ChemiDoc. For IP-kinase assays, HEK293 cells were transfected with flag-tagged wild-type STAT3 or STAT3-Y705F and underwent Flag IP. Flag IP products were subjected to TrkA kinase assay as described and reaction products were analyzed by Western blot. For kinase assay-IP, TrkA cell-free kinase assay was repeated and recombinant active TrkA was immunoprecipitated using anti-TrkA (Y490) antibody, resolved using SDS-PAGE and examined for co-immunoprecipitation of recombinant STAT3 using Western blot.
4.8. Pharmacological Inhibition of TrkA
MDA-MB-231 cells were seeded on a six-well culture dish and allowed to adhere overnight. Cells were treated with 5 μM Entrectinib (Adooq BioScience, Irvine, CA, USA) or 1% DMSO vehicle overnight at 37 °C and immediately harvested for Western blot analysis.
4.9. STAT3 Luciferase Assay
Luciferase assays were performed using a dual-luciferase assay kit (Biotium, Fremont, CA, USA) and readouts were collected using the Molecular Devices iD3 plate reader (San Jose, CA, USA). Cells were co-transfected with a STAT3-responsive luciferase reporter (pGAS-luc), and either vector or TrkA (pCMV5-TrkA) plasmid for 28 h before serum starvation for 16 h. Then, cells were stimulated with 100 ng/mL β-NGF for 4 h before the addition of lysis buffer. Lysates underwent mechanical lysis by scraping followed by trituration with a pipette to ensure complete lysis. Firefly luciferase activity was measured as previously described [56,63]. Results are represented as mean ± SD. Student’s t-tests were performed using GraphPad Prism 8.
4.10. Quantitative RT-PCR
Total RNA was isolated from cells using QIAGEN RNEasy Mini Kit (Hilden, Germany). RT-PCR was performed as previously described using primers listed in Supplementary Table S3 [62,64].
4.11. Flow Cytometry, Mammosphere Assays and ALDH Activity Assay
Flow cytometry for CD44high/CD24low cells and mammosphere assays were performed as previously described [65]. For ALDH activity, BT20 cells were transfected for 24 h and seeded in mammosphere conditions for up five days. Mammospheres were then pelleted and subjected to ALDH activity assay using the Abcam colorimetric ALDH activity kit (Cambridge, UK).
5. Conclusions
The current study is the first to illustrate that JAK2–STAT3 and TrkA pathways are co-activated in metastatic triple-negative and HER2-enriched breast cancers. Additionally, we show that TrkA interacts with and phosphorylates STAT3 on its Y705 residue, inducing oncogenic gene transcription to support breast cancer cell stemness, and that loss of Y705 residue abrogates STAT3 interaction with TrkA. These findings suggest that JAK2–STAT3 and TrkA signaling pathways converge at STAT3, supporting the crosstalk between these two signaling cascades. Importantly, we show that JAK2–STAT3/TrkA pathway co-activation is a new prognostic indicator of distant metastasis to various organs, prompting a future study to examine combined inhibition of JAK2–STAT3 and TrkA pathways as a novel and viable combination therapy in metastatic triple-negative and HER2-enriched breast cancers.
Conceptualization, H.-W.L.; methodology, H.-W.L., A.T.R.; formal analysis, H.-W.L., A.T.R.; investigation, H.-W.L.; resources, H.-W.L.; writing—original draft preparation, A.T.R., H.-W.L.; writing—review and editing, H.-W.L., N.R.A., S.A.J., D.L.D., S.G.M., D.Z., G.L.W., A.A., C.J.W., Y.Y., A.T., M.D.C., J.R., G.J., R.S., P.S., J.L.; supervision, H.-W.L.; project administration, H.-W.L.; funding acquisition, H.-W.L., A.T.R. All authors have read and agreed to the published version of the manuscript.
This research was funded by NIH grants R01NS087169 (HWL), 1T32CA247819-01 (ATR), P30CA012197 (BP), 1R01CA228137-01A1 (HWL), as well as, DoD grants, W81XWH-17-1-0044 (HWL), W81XWH-19-1-0072 (HWL), W81XWH-19-1-0753 (HWL), and W81XWH-20-1-0044 (HWL; JL).
No new datasets were generated in this manuscript.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. STAT3 and TrkA signaling pathways are significantly co-activated in triple-negative and HER2-enriched breast cancers. (A) IHC staining of 33 node-positive breast cancer tumors showed strong positive staining for both activated STAT3 (p-Y705) and activated TrkA (p-Y490) in both TNBC and HER2-enriched samples. Representative IHC images from three different subtypes of tumors are shown at the bottom. Arrows indicate positive staining for p-TrkA (Y490) at the cell membrane and nuclear staining for p-STAT3 (Y705). In a negative control experiment for IHC (no primary antibody, secondary antibody only), a mammary fat pad xenograft of BT474-TtzmR, a breast cancer cell line that is resistant to trastuzumab, was used. Magnification, 20×. (B) Immunofluorescent staining of patient breast tumor tissue for p-TrkA (Y490) and p-STAT3 (Y705). Scale bars, 50 μm. (C) Western blot panel of breast cancer cell lines to determine TrkA and STAT3 activation as indicated by p-TrkA and p-STAT3, respectively. (D,E) GSEA analysis of GEO breast cancer datasets (GSE 2034, 2603, 5327, 12,276) comparing STAT3 and/or TrkA activation signatures in TNBC versus luminal breast cancers (D) and HER2-enriched versus luminal tumors (E). (F–H) Activation scores for STAT3 (F), TrkA (G), and STAT3-TrkA (H) pathways in three major subtypes of breast cancer. GEO datasets were used. Original Western blots can be found in Supplementary Figure S3.
Enlarge this image.Figure 2. STAT3 and TrkA directly interact in cells and breast tumor tissue. (A) IP–Western blot (WB) in HEK293 cells transfected with a STAT3-Flag plasmid. Flag antibody (Ab) was used in IP. (B) Reciprocal IP-WB using p-TrkA Ab for IP. (C) IP–WB in MDA-MB-468 cells using p-TrkA (Y490) Ab for IP. (D) IP–WB of recombinant human TrkA and STAT3 following a cell-free TrkA kinase assay. (E) IP–WB of recombinant human TPM3–TrkA and STAT3 following cell-free kinase assay. (F) IP–WB of lysates from HEK293 cells transfected with STAT3-CA-Flag or STAT3-Y705F-Flag plasmids. Flag Ab was used in IP. (G) IP–WB of lysates from HEK293 cells co-transfected with STAT3-WT-Flag and either TrkA-WT or TrkA-K538N (kinase-dead). Densitometry values are displayed below each Western blot. Original Western blots are found in Supplementary Figures S4–S8. (H) Representative IF–confocal microscopy of HEK293 cells co-transfected with a TrkA-WT plasmid plus either a STAT3-WT-Flag or STAT3-Y705F-Flag plasmid, and stimulated with β-NGF for 20 min. DAPI, nuclei (blue). p-TrkA, green. p-STAT3, red. Merged signals, arrows. Scale bar, 10 μm. (I) IF staining of a MDA-MB-231 mammary fad xenograft (MFP) tumor for p-STAT3 (green) and p-TrkA (red). DAPI, nuclei (blue). Merged signals, arrows. Scale bars are 10 μm, 5 μm in magnified images.
Enlarge this image.Figure 3. TrkA phosphorylates STAT3 at Y705. (A) Cell-free TrkA kinase assay followed by WB to detect p-STAT3 (Y705), total STAT3, total TrkA, and total phosphotyrosines (p-Tyr). Recombinant human STAT3 and TrkA proteins were used. Coomassie-stained gel serves as a loading control for recombinant proteins used in each condition (arrows). (B) Immunoprecipitated STAT3-WT-Flag or STAT3-Y705F-Flag proteins were subjected to the TrkA kinase assay, followed by WB to detect p-STAT3 (Y705 or S727). Recombinant proteins are indicated with arrows. (C) Western blot of MDA-MB-231 TNBC cells upon treatment with Entrectinib. Densitometry values are displayed below each Western blot. Original Western blots are found in Supplementary Figures S9–S11. (D) LC/MS analysis of recombinant human STAT3 after TrkA kinase to detect p-STAT3 (Y705) phosphopeptides. (E) Representative IF–confocal microscopy showed that β-NGF induced TrkA activation, as indicated by p-TrkA signals (green), and STAT3 activation, as indicated by nuclear p-STAT3 (red). Scale bars, 10 μm.
Enlarge this image.Figure 4. TrkA phosphorylation promotes STAT3 nuclear import and transcriptional activity, and breast cancer stem cells. (A) β-NGF stimulation of untransfected SKBR3 breast cancer cells enhances STAT3 transcriptional activity, as shown by luciferase promoter activity. (B) TrkA activates STAT3 transcriptional activity as shown by the luciferase promoter activity. SKBR3 cells transfected with TrkA-WT or control vector plus pGAS-Luc (a luciferase reporter under the control of STAT3 response elements) were examined for luciferase activity. (C) β-NGF induces pGAS-Luc activity in the presence of TrkA STAT3 or TrkA-STAT3 overexpression. SKBR3 cells were transfected with STAT3-WT plus pGAS-Luc plasmids and with or without the TrkA plasmid, serum-starved, stimulated cells with or without β-NGF for 20 min. and analyzed for luciferase activity. SKBR3 cells transfected with control vector, TrkA-WT, STAT3, or TrkA–STAT3 plasmid plus pGAS-Luc were serum-starved the cells, stimulated with β-NGF for 20 min, and examined for luciferase activity. (D,E) Representative IF-confocal microscopy showed β-NGF induced activation of STAT3 in the presence of TrkA (D), but not kinase-dead TrkA-K538N mutant (E). HEK293 cells were used. Merged signals, arrows. Scale bars (upper right), 10 μm. (F,G) RT-PCR showed that TrkA-WT but not TrkA-K538N induces expression of STAT3 target genes, SOX2 (F) and MYC (G). SKBR3 were used. (H,I) RT-PCR of SOX2 (H) and MYC (I) in MDA-MB-468 upon treatment with a TrkA inhibitor Entrectinib. (J,K) Flow cytometry of transfected BT20 cells showed that TrkA-WT increases the CD44high/CD24low population of breast cancer cells (J) and that this effect is enhanced upon the addition of constitutively active STAT3 (K). (L,M) TrkA-WT (L) enhances the mammosphere-forming ability of BT20 cells, which is further increased upon the addition of constitutively active STAT3 (M). Representative mammosphere images are presented below each graph. Scale bars, 500 μm. (N,O) Representative ALDH activity of BT20 mammospheres upon transfection with TrkA alone (N) or TrkA and constitutively active STAT3 (O). All experiments were done for at least three times. Student’s t-test were used to determine p-values; NS: not significant. Data represented as mean ± SD. Western blots and densitometry values for each transfection condition are displayed below each plot. Original Western blots are found in Supplementary Figures S12 and S14.
Enlarge this image.Figure 5. Co-activation of JAK2–STAT3 and TrkA pathways is correlated with poor overall and bone metastasis-free survival of triple-negative and HER2-enriched breast cancers. Using Kaplan–Meier analysis log-rank test, and the STAT3 or TrkA activation signatures, we analyzed patients whose expression profiles and survival data were obtained from GEO databases (GSE 2034, 2603, 5327, 12,276). (A) Overall metastasis-free survival in 166 triple-negative breast cancer. (B) Overall metastasis-free survival in 96 HER2 breast cancer. (C) Bone metastasis-free survival in 166 triple-negative breast cancer patients. (D) Bone metastasis-free survival in 96 HER2 breast cancer patients. Median survival times are indicated.
Enlarge this image.Figure 6. Co-activation of JAK2-STAT3 and TrkA pathways is associated with poor brain- and lung metastasis-free survival of triple-negative and HER2-enriched breast cancers. Using Kaplan–Meier analysis log-rank test, and the STAT3 or TrkA activation signatures, we analyzed patients whose expression profiles and survival data were obtained from GEO databases (GSE 2034, 2603, 5327, 12,276). (A) Brain metastasis-free survival in 166 triple-negative breast cancer patients. (B) Brain metastasis-free survival in HER2 breast cancer. (C) Lung metastasis-free survival in triple-negative breast cancer. (D) Lung metastasis-free survival in 96 HER2-enriched breast cancer patients. Median survival times are indicated.
Supplementary Materials
The following are available online at
References
1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin.; 2020; 70, pp. 7-30. [DOI: https://dx.doi.org/10.3322/caac.21590] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31912902]
2. de Azambuja, E.; Holmes, A.P.; Piccart-Gebhart, M.; Holmes, E.; Di Cosimo, S.; Swaby, R.F.; Untch, M.; Jackisch, C.; Lang, I.; Smith, I. et al. Lapatinib with trastuzumab for HER2-positive early breast cancer (NeoALTTO): Survival outcomes of a randomised, open-label, multicentre, phase 3 trial and their association with pathological complete response. Lancet Oncol.; 2014; 15, pp. 1137-1146. [DOI: https://dx.doi.org/10.1016/S1470-2045(14)70320-1]
3. Harbeck, N.; Gnant, M. Breast cancer. Lancet; 2017; 389, pp. 1134-1150. [DOI: https://dx.doi.org/10.1016/S0140-6736(16)31891-8]
4. Dent, R.; Valentini, A.; Hanna, W.; Rawlinson, E.; Rakovitch, E.; Sun, P.; Narod, S.A. Factors associated with breast cancer mortality after local recurrence. Curr. Oncol.; 2014; 21, pp. e418-e425. [DOI: https://dx.doi.org/10.3747/co.21.1563]
5. Balko, J.M.; Schwarz, L.J.; Luo, N.; Estrada, M.V.; Giltnane, J.M.; Davila-Gonzalez, D.; Wang, K.; Sanchez, V.; Dean, P.T.; Combs, S.E. et al. Triple-negative breast cancers with amplification of JAK2 at the 9p24 locus demonstrate JAK2-specific dependence. Sci. Transl. Med.; 2016; 8, 334ra353. [DOI: https://dx.doi.org/10.1126/scitranslmed.aad3001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27075627]
6. Duru, N.; Fan, M.; Candas, D.; Menaa, C.; Liu, H.C.; Nantajit, D.; Wen, Y.; Xiao, K.; Eldridge, A.; Chromy, B.A. et al. HER2-associated radioresistance of breast cancer stem cells isolated from HER2-negative breast cancer cells. Clin. Cancer Res.; 2012; 18, pp. 6634-6647. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-12-1436]
7. Yu, H.; Jove, R. The STATs of cancer--new molecular targets come of age. Nat. Rev. Cancer; 2004; 4, pp. 97-105. [DOI: https://dx.doi.org/10.1038/nrc1275] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14964307]
8. Teng, T.S.; Lin, B.; Manser, E.; Ng, D.C.; Cao, X. Stat3 promotes directional cell migration by regulating Rac1 activity via its activator betaPIX. J. Cell Sci.; 2009; 122, pp. 4150-4159. [DOI: https://dx.doi.org/10.1242/jcs.057109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19861492]
9. Xiong, H.; Zhang, Z.G.; Tian, X.Q.; Sun, D.F.; Liang, Q.C.; Zhang, Y.J.; Lu, R.; Chen, Y.X.; Fang, J.Y. Inhibition of JAK1, 2/STAT3 signaling induces apoptosis, cell cycle arrest, and reduces tumor cell invasion in colorectal cancer cells. Neoplasia; 2008; 10, pp. 287-297. [DOI: https://dx.doi.org/10.1593/neo.07971] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18320073]
10. Lo, H.W.; Cao, X.; Zhu, H.; Ali-Osman, F. Constitutively activated STAT3 frequently coexpresses with epidermal growth factor receptor in high-grade gliomas and targeting STAT3 sensitizes them to Iressa and alkylators. Clin. Cancer Res.; 2008; 14, pp. 6042-6054. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-07-4923]
11. Lo, H.W.; Cao, X.; Zhu, H.; Ali-Osman, F. Cyclooxygenase-2 is a novel transcriptional target of the nuclear EGFR-STAT3 and EGFRvIII-STAT3 signaling axes. Mol. Cancer Res.; 2010; 8, pp. 232-245. [DOI: https://dx.doi.org/10.1158/1541-7786.MCR-09-0391]
12. Lo, H.W.; Hsu, S.C.; Ali-Seyed, M.; Gunduz, M.; Xia, W.; Wei, Y.; Bartholomeusz, G.; Shih, J.Y.; Hung, M.C. Nuclear interaction of EGFR and STAT3 in the activation of the iNOS/NO pathway. Cancer Cell; 2005; 7, pp. 575-589. [DOI: https://dx.doi.org/10.1016/j.ccr.2005.05.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15950906]
13. Lo, H.W.; Hsu, S.C.; Xia, W.; Cao, X.; Shih, J.Y.; Wei, Y.; Abbruzzese, J.L.; Hortobagyi, G.N.; Hung, M.C. Epidermal growth factor receptor cooperates with signal transducer and activator of transcription 3 to induce epithelial-mesenchymal transition in cancer cells via up-regulation of TWIST gene expression. Cancer Res.; 2007; 67, pp. 9066-9076. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-07-0575]
14. Han, W.; Carpenter, R.L.; Cao, X.; Lo, H.W. STAT1 gene expression is enhanced by nuclear EGFR and HER2 via cooperation with STAT3. Mol. Carcinog.; 2013; 52, pp. 959-969. [DOI: https://dx.doi.org/10.1002/mc.21936]
15. Mesa, R.A.; Yasothan, U.; Kirkpatrick, P. Ruxolitinib. Nat. Rev. Drug Discov.; 2012; 11, pp. 103-104. [DOI: https://dx.doi.org/10.1038/nrd3652]
16. Mullard, A. FDA approves Eli Lilly’s baricitinib. Nat. Rev. Drug Discov.; 2018; 17, 460. [DOI: https://dx.doi.org/10.1038/nrd.2018.112] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29950712]
17. Fedratinib Becomes New Option in Myelofibrosis. Cancer Discov.; 2019; 9, 1332. [DOI: https://dx.doi.org/10.1158/2159-8290.CD-NB2019-102]
18. Lange, A.M.; Lo, H.W. Inhibiting TRK Proteins in Clinical Cancer Therapy. Cancers; 2018; 10, 105. [DOI: https://dx.doi.org/10.3390/cancers10040105]
19. Regua, A.T.; Doheny, D.; Arrigo, A.; Lo, H.W. Trk receptor tyrosine kinases in metastasis and cancer therapy. Discov. Med.; 2019; 28, pp. 195-203.
20. Lagadec, C.; Meignan, S.; Adriaenssens, E.; Foveau, B.; Vanhecke, E.; Romon, R.; Toillon, R.A.; Oxombre, B.; Hondermarck, H.; Le Bourhis, X. TrkA overexpression enhances growth and metastasis of breast cancer cells. Oncogene; 2009; 28, pp. 1960-1970. [DOI: https://dx.doi.org/10.1038/onc.2009.61]
21. Kyker-Snowman, K.; Hughes, R.M.; Yankaskas, C.L.; Cravero, K.; Karthikeyan, S.; Button, B.; Waters, I.; Rosen, D.M.; Dennison, L.; Hunter, N. et al. TrkA overexpression in non-tumorigenic human breast cell lines confers oncogenic and metastatic properties. Breast Cancer Res. Treat.; 2020; 179, pp. 631-642. [DOI: https://dx.doi.org/10.1007/s10549-019-05506-3]
22. Demont, Y.; Corbet, C.; Page, A.; Ataman-Onal, Y.; Choquet-Kastylevsky, G.; Fliniaux, I.; Le Bourhis, X.; Toillon, R.A.; Bradshaw, R.A.; Hondermarck, H. Pro-nerve growth factor induces autocrine stimulation of breast cancer cell invasion through tropomyosin-related kinase A (TrkA) and sortilin protein. J. Biol. Chem.; 2012; 287, pp. 1923-1931. [DOI: https://dx.doi.org/10.1074/jbc.M110.211714]
23. Descamps, S.; Toillon, R.A.; Adriaenssens, E.; Pawlowski, V.; Cool, S.M.; Nurcombe, V.; Le Bourhis, X.; Boilly, B.; Peyrat, J.P.; Hondermarck, H. Nerve growth factor stimulates proliferation and survival of human breast cancer cells through two distinct signaling pathways. J. Biol. Chem.; 2001; 276, pp. 17864-17870. [DOI: https://dx.doi.org/10.1074/jbc.M010499200]
24. Zhang, J.; Wang, L.S.; Ye, S.L.; Luo, P.; Wang, B.L. Blockage of tropomyosin receptor kinase a (TrkA) enhances chemo-sensitivity in breast cancer cells and inhibits metastasis in vivo. Int. J. Clin. Exp. Med.; 2015; 8, pp. 634-641.
25. Laetsch, T.W.; Hawkins, D.S. Larotrectinib for the treatment of TRK fusion solid tumors. Expert Rev. Anticancer Ther.; 2019; 19, pp. 1-10. [DOI: https://dx.doi.org/10.1080/14737140.2019.1538796]
26. Entrectinib OK’d for Cancers with NTRK Fusions, NSCLC. Cancer Discov.; 2019; 9, OF2. [DOI: https://dx.doi.org/10.1158/2159-8290.CD-NB2019-101] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31471291]
27. Miranda, C.; Fumagalli, T.; Anania, M.C.; Vizioli, M.G.; Pagliardini, S.; Pierotti, M.A.; Greco, A. Role of STAT3 in in vitro transformation triggered by TRK oncogenes. PLoS ONE; 2010; 5, e9446. [DOI: https://dx.doi.org/10.1371/journal.pone.0009446] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20209132]
28. Ng, Y.P.; Cheung, Z.H.; Ip, N.Y. STAT3 as a downstream mediator of Trk signaling and functions. J. Biol. Chem.; 2006; 281, pp. 15636-15644. [DOI: https://dx.doi.org/10.1074/jbc.M601863200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16611639]
29. Gatza, M.L.; Lucas, J.E.; Barry, W.T.; Kim, J.W.; Wang, Q.; Crawford, M.D.; Datto, M.B.; Kelley, M.; Mathey-Prevot, B.; Potti, A. et al. A pathway-based classification of human breast cancer. Proc. Natl. Acad. Sci. USA; 2010; 107, pp. 6994-6999. [DOI: https://dx.doi.org/10.1073/pnas.0912708107] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20335537]
30. Liberzon, A.; Birger, C.; Thorvaldsdottir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst.; 2015; 1, pp. 417-425. [DOI: https://dx.doi.org/10.1016/j.cels.2015.12.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26771021]
31. Huang, E.J.; Reichardt, L.F. Trk receptors: Roles in neuronal signal transduction. Annu. Rev. Biochem.; 2003; 72, pp. 609-642. [DOI: https://dx.doi.org/10.1146/annurev.biochem.72.121801.161629]
32. Friedman, W.J.; Greene, L.A. Neurotrophin signaling via Trks and p75. Exp. Cell Res.; 1999; 253, pp. 131-142. [DOI: https://dx.doi.org/10.1006/excr.1999.4705] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10579918]
33. Kaplan, D.R.; Miller, F.D. Neurotrophin signal transduction in the nervous system. Curr. Opin. Neurobiol.; 2000; 10, pp. 381-391. [DOI: https://dx.doi.org/10.1016/S0959-4388(00)00092-1]
34. Jassal, B.; Matthews, L.; Viteri, G.; Gong, C.; Lorente, P.; Fabregat, A.; Sidiropoulos, K.; Cook, J.; Gillespie, M.; Haw, R. et al. The reactome pathway knowledgebase. Nucleic Acids Res.; 2020; 48, pp. D498-D503. [DOI: https://dx.doi.org/10.1093/nar/gkz1031]
35. Dauer, D.J.; Ferraro, B.; Song, L.; Yu, B.; Mora, L.; Buettner, R.; Enkemann, S.; Jove, R.; Haura, E.B. Stat3 regulates genes common to both wound healing and cancer. Oncogene; 2005; 24, pp. 3397-3408. [DOI: https://dx.doi.org/10.1038/sj.onc.1208469] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15735721]
36. Ardini, E.; Bosotti, R.; Borgia, A.L.; De Ponti, C.; Somaschini, A.; Cammarota, R.; Amboldi, N.; Raddrizzani, L.; Milani, A.; Magnaghi, P. et al. The TPM3-NTRK1 rearrangement is a recurring event in colorectal carcinoma and is associated with tumor sensitivity to TRKA kinase inhibition. Mol. Oncol.; 2014; 8, pp. 1495-1507. [DOI: https://dx.doi.org/10.1016/j.molonc.2014.06.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24962792]
37. Eggert, A.; Ikegaki, N.; Liu, X.; Chou, T.T.; Lee, V.M.; Trojanowski, J.Q.; Brodeur, G.M. Molecular dissection of TrkA signal transduction pathways mediating differentiation in human neuroblastoma cells. Oncogene; 2000; 19, pp. 2043-2051. [DOI: https://dx.doi.org/10.1038/sj.onc.1203518] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10803465]
38. Ardini, E.; Menichincheri, M.; Banfi, P.; Bosotti, R.; De Ponti, C.; Pulci, R.; Ballinari, D.; Ciomei, M.; Texido, G.; Degrassi, A. et al. Entrectinib, a Pan-TRK, ROS1, and ALK Inhibitor with Activity in Multiple Molecularly Defined Cancer Indications. Mol. Cancer Ther.; 2016; 15, pp. 628-639. [DOI: https://dx.doi.org/10.1158/1535-7163.MCT-15-0758] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26939704]
39. Smith, K.M.; Fagan, P.C.; Pomari, E.; Germano, G.; Frasson, C.; Walsh, C.; Silverman, I.; Bonvini, P.; Li, G. Antitumor Activity of Entrectinib, a Pan-TRK, ROS1, and ALK Inhibitor, in ETV6-NTRK3-Positive Acute Myeloid Leukemia. Mol. Cancer Ther.; 2018; 17, pp. 455-463. [DOI: https://dx.doi.org/10.1158/1535-7163.MCT-17-0419] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29237803]
40. Carpenter, R.L.; Lo, H.W. STAT3 Target Genes Relevant to Human Cancers. Cancers; 2014; 6, pp. 897-925. [DOI: https://dx.doi.org/10.3390/cancers6020897]
41. Foshay, K.M.; Gallicano, G.I. Regulation of Sox2 by STAT3 initiates commitment to the neural precursor cell fate. Stem Cells Dev.; 2008; 17, pp. 269-278. [DOI: https://dx.doi.org/10.1089/scd.2007.0098] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18447642]
42. Kiuchi, N.; Nakajima, K.; Ichiba, M.; Fukada, T.; Narimatsu, M.; Mizuno, K.; Hibi, M.; Hirano, T. STAT3 is required for the gp130-mediated full activation of the c-myc gene. J. Exp. Med.; 1999; 189, pp. 63-73. [DOI: https://dx.doi.org/10.1084/jem.189.1.63] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9874564]
43. Ginestier, C.; Hur, M.H.; Charafe-Jauffret, E.; Monville, F.; Dutcher, J.; Brown, M.; Jacquemier, J.; Viens, P.; Kleer, C.G.; Liu, S. et al. ALDH1 is a marker of normal and malignant human mammary stem cells and a predictor of poor clinical outcome. Cell Stem Cell; 2007; 1, pp. 555-567. [DOI: https://dx.doi.org/10.1016/j.stem.2007.08.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18371393]
44. Sakamoto, K.; Triplett, A.A.; Schuler, L.A.; Wagner, K.U. Janus kinase 2 is required for the initiation but not maintenance of prolactin-induced mammary cancer. Oncogene; 2010; 29, pp. 5359-5369. [DOI: https://dx.doi.org/10.1038/onc.2010.274] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20639901]
45. Wehde, B.L.; Radler, P.D.; Shrestha, H.; Johnson, S.J.; Triplett, A.A.; Wagner, K.U. Janus Kinase 1 Plays a Critical Role in Mammary Cancer Progression. Cell Rep.; 2018; 25, pp. 2192-2207. [DOI: https://dx.doi.org/10.1016/j.celrep.2018.10.063] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30463015]
46. Hart, S.; Goh, K.C.; Novotny-Diermayr, V.; Tan, Y.C.; Madan, B.; Amalini, C.; Ong, L.C.; Kheng, B.; Cheong, A.; Zhou, J. et al. Pacritinib (SB1518), a JAK2/FLT3 inhibitor for the treatment of acute myeloid leukemia. Blood Cancer J.; 2011; 1, e44. [DOI: https://dx.doi.org/10.1038/bcj.2011.43] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22829080]
47. Berdeja, J.; Palandri, F.; Baer, M.R.; Quick, D.; Kiladjian, J.J.; Martinelli, G.; Verma, A.; Hamid, O.; Walgren, R.; Pitou, C. et al. Phase 2 study of gandotinib (LY2784544) in patients with myeloproliferative neoplasms. Leuk. Res.; 2018; 71, pp. 82-88. [DOI: https://dx.doi.org/10.1016/j.leukres.2018.06.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30025280]
48. Mascarenhas, J.; Baer, M.R.; Kessler, C.; Hexner, E.; Tremblay, D.; Price, L.; Sandy, L.; Weinberg, R.; Pahl, H.; Silverman, L.R. et al. Phase II trial of Lestaurtinib, a JAK2 inhibitor, in patients with myelofibrosis. Leuk. Lymphoma; 2019; 60, pp. 1343-1345. [DOI: https://dx.doi.org/10.1080/10428194.2018.1532509]
49. Verstovsek, S.; Courby, S.; Griesshammer, M.; Mesa, R.A.; Brachmann, C.B.; Kawashima, J.; Maltzman, J.D.; Shao, L.; Xin, Y.; Huang, D. et al. A phase 2 study of momelotinib, a potent JAK1 and JAK2 inhibitor, in patients with polycythemia vera or essential thrombocythemia. Leuk. Res.; 2017; 60, pp. 11-17. [DOI: https://dx.doi.org/10.1016/j.leukres.2017.05.002]
50. Leis, O.; Eguiara, A.; Lopez-Arribillaga, E.; Alberdi, M.J.; Hernandez-Garcia, S.; Elorriaga, K.; Pandiella, A.; Rezola, R.; Martin, A.G. Sox2 expression in breast tumours and activation in breast cancer stem cells. Oncogene; 2012; 31, pp. 1354-1365. [DOI: https://dx.doi.org/10.1038/onc.2011.338]
51. Wang, Z.; Oron, E.; Nelson, B.; Razis, S.; Ivanova, N. Distinct lineage specification roles for NANOG, OCT4, and SOX2 in human embryonic stem cells. Cell Stem Cell; 2012; 10, pp. 440-454. [DOI: https://dx.doi.org/10.1016/j.stem.2012.02.016]
52. Butti, R.; Gunasekaran, V.P.; Kumar, T.V.S.; Banerjee, P.; Kundu, G.C. Breast cancer stem cells: Biology and therapeutic implications. Int. J. Biochem. Cell Biol.; 2019; 107, pp. 38-52. [DOI: https://dx.doi.org/10.1016/j.biocel.2018.12.001]
53. Toledo-Guzman, M.E.; Hernandez, M.I.; Gomez-Gallegos, A.A.; Ortiz-Sanchez, E. ALDH as a Stem Cell Marker in Solid Tumors. Curr. Stem Cell Res. Ther.; 2019; 14, pp. 375-388. [DOI: https://dx.doi.org/10.2174/1574888X13666180810120012] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30095061]
54. Yin, S.; Cheryan, V.T.; Xu, L.; Rishi, A.K.; Reddy, K.B. Myc mediates cancer stem-like cells and EMT changes in triple negative breast cancers cells. PLoS ONE; 2017; 12, e0183578. [DOI: https://dx.doi.org/10.1371/journal.pone.0183578] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28817737]
55. Zhu, H.; Carpenter, R.L.; Han, W.; Lo, H.W. The GLI1 splice variant TGLI1 promotes glioblastoma angiogenesis and growth. Cancer Lett.; 2014; 343, pp. 51-61. [DOI: https://dx.doi.org/10.1016/j.canlet.2013.09.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24045042]
56. Carpenter, R.L.; Paw, I.; Dewhirst, M.W.; Lo, H.W. Akt phosphorylates and activates HSF-1 independent of heat shock, leading to Slug overexpression and epithelial-mesenchymal transition (EMT) of HER2-overexpressing breast cancer cells. Oncogene; 2015; 34, pp. 546-557. [DOI: https://dx.doi.org/10.1038/onc.2013.582]
57. Yano, H.; Cong, F.; Birge, R.B.; Goff, S.P.; Chao, M.V. Association of the Abl tyrosine kinase with the Trk nerve growth factor receptor. J. Neurosci. Res.; 2000; 59, pp. 356-364. [DOI: https://dx.doi.org/10.1002/(SICI)1097-4547(20000201)59:3<356::AID-JNR9>3.0.CO;2-G]
58. Dasgupta, M.; Unal, H.; Willard, B.; Yang, J.; Karnik, S.S.; Stark, G.R. Critical role for lysine 685 in gene expression mediated by transcription factor unphosphorylated STAT3. J. Biol. Chem.; 2014; 289, pp. 30763-30771. [DOI: https://dx.doi.org/10.1074/jbc.M114.603894]
59. Wei, D.; Le, X.; Zheng, L.; Wang, L.; Frey, J.A.; Gao, A.C.; Peng, Z.; Huang, S.; Xiong, H.Q.; Abbruzzese, J.L. et al. Stat3 activation regulates the expression of vascular endothelial growth factor and human pancreatic cancer angiogenesis and metastasis. Oncogene; 2003; 22, pp. 319-329. [DOI: https://dx.doi.org/10.1038/sj.onc.1206122]
60. Chang, J.T.; Gatza, M.L.; Lucas, J.E.; Barry, W.T.; Vaughn, P.; Nevins, J.R. SIGNATURE: A workbench for gene expression signature analysis. BMC Bioinform.; 2011; 12, 443. [DOI: https://dx.doi.org/10.1186/1471-2105-12-443]
61. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S. et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 15545-15550. [DOI: https://dx.doi.org/10.1073/pnas.0506580102] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16199517]
62. Sirkisoon, S.R.; Carpenter, R.L.; Rimkus, T.; Anderson, A.; Harrison, A.; Lange, A.M.; Jin, G.; Watabe, K.; Lo, H.W. Interaction between STAT3 and GLI1/tGLI1 oncogenic transcription factors promotes the aggressiveness of triple-negative breast cancers and HER2-enriched breast cancer. Oncogene; 2018; 37, pp. 2502-2514. [DOI: https://dx.doi.org/10.1038/s41388-018-0132-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29449694]
63. Lo, H.W.; Zhu, H.; Cao, X.; Aldrich, A.; Ali-Osman, F. A novel splice variant of GLI1 that promotes glioblastoma cell migration and invasion. Cancer Res.; 2009; 69, pp. 6790-6798. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-09-0886]
64. Sirkisoon, S.R.; Carpenter, R.L.; Rimkus, T.; Doheny, D.; Zhu, D.; Aguayo, N.R.; Xing, F.; Chan, M.; Ruiz, J.; Metheny-Barlow, L.J. et al. TGLI1 transcription factor mediates breast cancer brain metastasis via activating metastasis-initiating cancer stem cells and astrocytes in the tumor microenvironment. Oncogene; 2020; 39, pp. 64-78. [DOI: https://dx.doi.org/10.1038/s41388-019-0959-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31462709]
65. Doheny, D.; Sirkisoon, S.; Carpenter, R.L.; Aguayo, N.R.; Regua, A.T.; Anguelov, M.; Manore, S.G.; Arrigo, A.; Jalboush, S.A.; Wong, G.L. et al. Combined inhibition of JAK2-STAT3 and SMO-GLI1/tGLI1 pathways suppresses breast cancer stem cells, tumor growth, and metastasis. Oncogene; 2020; [DOI: https://dx.doi.org/10.1038/s41388-020-01454-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32929154]
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Aguayo, Noah R 1 ; Sara Abu Jalboush 1 ; Doheny, Daniel L 1
; Manore, Sara G 1
; Zhu, Dongqin 1 ; Wong, Grace L 1
; Austin, Arrigo 1 ; Wagner, Calvin J 1 ; Yang, Yu 1 ; Thomas, Alexandra 2 ; Chan, Michael D 3 ; Ruiz, Jimmy 4 ; Jin, Guangxu 5 ; Strowd, Roy 6 ; Sun, Peiqing 5 ; Lin, Jiayuh 7 ; Hui-Wen, Lo 8 1 Department of Cancer Biology, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
2 Department of Hematology and Oncology, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
3 Wake Forest Baptist Comprehensive Cancer Center, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
4 Department of Hematology and Oncology, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
5 Department of Cancer Biology, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
6 Wake Forest Baptist Comprehensive Cancer Center, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;
7 Department of Biochemistry and Molecular Biology, University of Maryland School of Medicine, Baltimore, MD 21201, USA;
8 Department of Cancer Biology, Wake Forest University School of Medicine, Winston-Salem, NC 27101, USA;