Proteomic characterization of esophageal squamous cell carcinoma response to immunotherapy reveals potential therapeutic strategy and predictive biomarkers.

Immunotherapy is the first-line therapy for esophageal squamous cell carcinoma (ESCC), yet many patients do not respond due to drug resistance and the lack of reliable predictive markers. We collected 73 ESCC patients (including discovery cohort and validation cohort) without immune thrombocytopenia and undergoing anti-PD1 immunotherapy. Proteomic and phosphoproteomic analysis of 73 ESCC treatment-naive samples by mass spectrometry-based label-free quantification were applied to explore the potential resistant and sensitive mechanisms, and identify predictive markers of ESCC immunotherapy. Comparative analysis found the pathways related to immune and mitochondrial functions were associated with ESCC immunotherapy sensitivity; while platelet activation bioprocess showed negative correlation with CD8+ T cells and related to ESCC immunotherapy non-sensitivity. Finally, we identified 10 ESCC immunotherapy predictive biomarkers with high accuracy (≥ 0.90) to predict the immunotherapeutic response, which was validated in the independent cohort.

Longitudinal plasma proteome profiling reveals the diversity of biomarkers for diagnosis and cetuximab therapy response of colorectal cancer.

Cetuximab therapy is the major treatment for colorectal cancer (CRC), but drug resistance limits its effectiveness. Here, we perform longitudinal and deep proteomic profiling of 641 plasma samples originated from 147 CRC patients (CRCs) undergoing cetuximab therapy with multi-course treatment, and 90 healthy controls (HCs). COL12A1, THBS2, S100A8, and S100A9 are screened as potential proteins to distinguish CRCs from HCs both in plasma and tissue validation cohorts. We identify the potential biomarkers (RRAS2, MMP8, FBLN1, RPTOR, and IMPDH2) for the initial response prediction. In a longitudinal setting, we identify two clusters with distinct fluctuations and construct the model with high accuracy to predict the longitudinal response, further validated in the independent cohort. This study reveals the heterogeneity of different biomarkers for tumor diagnosis, the initial and longitudinal response prediction respectively in the first course and multi-course cetuximab treatment, may ultimately be useful in monitoring and intervention strategies for CRC.

Proteomic characterization of the colorectal cancer response to chemoradiation and targeted therapies reveals potential therapeutic strategies.

Chemoradiation and targeted therapies are the major treatments for colorectal cancer (CRC); however, molecular properties associated with therapy resistance are incompletely characterized. Here, we profile the proteome of 254 tumor tissues from patients with CRC undergoing chemotherapy, chemoradiation, or chemotherapy combined with targeted therapy. Proteome-based classification reveals four subtypes featured with distinct biological and therapeutic characteristics. The integrative analysis of CRC cell lines and clinical samples indicates that immune regulation is significantly associated with drug sensitivity. HSF1 can increase DNA damage repair and cell cycle, thus inducing resistance to radiation, while high expression of HDAC6 is negatively associated with response of cetuximab. Furthermore, we develop prognostic models with high accuracy to predict the therapeutic response, further validated by parallel reaction monitoring (PRM) assay in an independent validation cohort. This study provides a rich resource for investigating the mechanisms and indicators of chemoradiation and targeted therapy in CRC.

Proteogenomics of different urothelial bladder cancer stages reveals distinct molecular features for papillary cancer and carcinoma in situ.

The progression of urothelial bladder cancer (UC) is a complicated multi-step process. We perform a comprehensive multi-omics analysis of 448 samples from 190 UC patients, covering the whole spectrum of disease stages and grades. Proteogenomic integration analysis indicates the mutations of HRAS regulated mTOR signaling to form urothelial papilloma rather than papillary urothelial cancer (PUC). DNA damage is a key signaling pathway in the progression of carcinoma in situ (CIS) and related to APOBEC signature. Glucolipid metabolism increase and lower immune cell infiltration are associated with PUC compared to CIS. Proteomic analysis distinguishes the origins of invasive tumors (PUC-derived and CIS-derived), related to distinct clinical prognosis and molecular features. Additionally, loss of RBPMS, associated with CIS-derived tumors, is validated to increase the activity of AP-1 and promote metastasis. This study reveals the characteristics of two distinct branches (PUC and CIS) of UC progression and may eventually benefit clinical practice.

Plasma proteomic profiling discovers molecular features associated with upper tract urothelial carcinoma.

Upper tract urothelial carcinoma (UTUC) is often diagnosed late and exhibits poor prognosis. Limited data are available on potential non-invasive biomarkers for disease monitoring. Here, we investigate the proteomic profile of plasma in 362 UTUC patients and 239 healthy controls. We present an integrated tissue-plasma proteomic approach to infer the signature proteins for identifying patients with muscle-invasive UTUC. We discover a protein panel that reflects lymph node metastasis, which is of interest in identifying UTUC patients with high risk and poor prognosis. We also identify a ten-protein classifier and establish a progression clock predicting progression-free survival of UTUC patients. Finally, we further validate the signature proteins by parallel reaction monitoring assay in an independent cohort. Collectively, this study portrays the plasma proteomic landscape of a UTUC cohort and provides a valuable resource for further biological and diagnostic research in UTUC.

Modulation of cellular metabolism by protein crotonylation regulates pancreatic cancer progression.

Protein lysine crotonylation has been recently identified as a vital posttranslational modification in cellular processes, particularly through the modification of histones. We show that lysine crotonylation is an important modification of the cytoplastic and mitochondria proteins. Enzymes in glycolysis, the tricarboxylic acid (TCA) cycle, fatty acid metabolism, glutamine metabolism, glutathione metabolism, the urea cycle, one-carbon metabolism, and mitochondrial fusion/fission dynamics are found to be extensively crotonylated in pancreatic cancer cells. This modulation is mainly controlled by a pair of crotonylation writers and erasers including CBP/p300, HDAC1, and HDAC3. The dynamic crotonylation of metabolic enzymes is involved in metabolism regulation, which is linked with tumor progression. Interestingly, the activation of MTHFD1 by decrotonylation at Lys354 and Lys553 promotes the development of pancreatic cancer by increasing resistance to ferroptosis. Our study suggests that crotonylation represents a metabolic regulatory mechanism in pancreatic cancer progression.

Proteomic characteristics reveal the signatures and the risks of T1 colorectal cancer metastasis to lymph nodes.

The presence of lymph node metastasis (LNM) affects treatment strategy decisions in T1NxM0 colorectal cancer (CRC), but the currently used clinicopathological-based risk stratification cannot predict LNM accurately. In this study, we detected proteins in formalin-fixed paraffin-embedded (FFPE) tumor samples from 143 LNM-negative and 78 LNM-positive patients with T1 CRC and revealed changes in molecular and biological pathways by label-free liquid chromatography tandem mass spectrometry (LC-MS/MS) and established classifiers for predicting LNM in T1 CRC. An effective 55-proteins prediction model was built by machine learning and validated in a training cohort (N=132) and two validation cohorts (VC1, N=42; VC2, N=47), achieved an impressive AUC of 1.00 in the training cohort, 0.96 in VC1 and 0.93 in VC2, respectively. We further built a simplified classifier with nine proteins, and achieved an AUC of 0.824. The simplified classifier was performed excellently in two external validation cohorts. The expression patterns of 13 proteins were confirmed by immunohistochemistry, and the IHC score of five proteins was used to build an IHC predict model with an AUC of 0.825. RHOT2 silence significantly enhanced migration and invasion of colon cancer cells. Our study explored the mechanism of metastasis in T1 CRC and can be used to facilitate the individualized prediction of LNM in patients with T1 CRC, which may provide a guidance for clinical practice in T1 CRC.

Most patients with early-stage colorectal cancer can be treated with a minimally invasive procedure. Surgeons use a flexible tool to remove precancerous or cancerous cells, cutting the risk of death from colorectal cancer in half. But a small number of early-stage colorectal cancer patients are at risk of their cancer spreading to the lymph nodes. These patients need more extensive surgery. Clinicians use risk stratification tools to decide which patients need more extensive surgery. Unfortunately, the existing risk stratification tools are not very accurate. The current approach, which analyzes colon tissue for cancerous changes, classifies 70% to 80% of early-stage colorectal cancer patients as high risk for cancer spread. But only about 8% to 16% of patients in the high risk group have lymph node metastasis. As a result, many patients undergo unnecessary, invasive surgery. Zhuang, Zhuang, Chen, Qin, et al. developed a more accurate way to predict which patients are at risk of lymph node metastasis using proteins. In the experiments, the team analyzed the proteins in tumor samples from 143 patients with early colorectal cancer who did not have lymph node metastases and 78 patients with metastases. Zhuang et al. then used machine learning to build a prediction tool that used 55 proteins to identify patients at risk of metastases. The new approach was more accurate than existing tools and simplified versions with only nine or five proteins also performed better than existing tools. This work provides preliminary evidence that protein-based models using as few as five proteins can more accurately identify which patients are at risk of metastasis. These models may reduce the number of patients who undergo unnecessary invasive surgery. The experiments also identified potential targets for therapies to prevent or treat lymph metastases. For example, they showed that low levels of the RHOT2 protein predict metastasis.

Comprehensive proteogenomic characterization of early duodenal cancer reveals the carcinogenesis tracks of different subtypes.

The subtypes of duodenal cancer (DC) are complicated and the carcinogenesis process is not well characterized. We present comprehensive characterization of 438 samples from 156 DC patients, covering 2 major and 5 rare subtypes. Proteogenomics reveals LYN amplification at the chromosome 8q gain functioned in the transmit from intraepithelial neoplasia phase to infiltration tumor phase via MAPK signaling, and illustrates the DST mutation improves mTOR signaling in the duodenal adenocarcinoma stage. Proteome-based analysis elucidates stage-specific molecular characterizations and carcinogenesis tracks, and defines the cancer-driving waves of the adenocarcinoma and Brunner's gland subtypes. The drug-targetable alanyl-tRNA synthetase (AARS1) in the high tumor mutation burden/immune infiltration is significantly enhanced in DC progression, and catalyzes the lysine-alanylation of poly-ADP-ribose polymerases (PARP1), which decreases the apoptosis of cancer cells, eventually promoting cell proliferation and tumorigenesis. We assess the proteogenomic landscape of early DC, and provide insights into the molecular features corresponding therapeutic targets.

Integrative proteogenomic characterization of early esophageal cancer.

Esophageal squamous cell carcinoma (ESCC) is malignant while the carcinogenesis is still unclear. Here, we perform a comprehensive multi-omics analysis of 786 trace-tumor-samples from 154 ESCC patients, covering 9 histopathological stages and 3 phases. Proteogenomics elucidates cancer-driving waves in ESCC progression, and reveals the molecular characterization of alcohol drinking habit associated signatures. We discover chromosome 3q gain functions in the transmit from nontumor to intraepithelial neoplasia phases, and find TP53 mutation enhances DNA replication in intraepithelial neoplasia phase. The mutations of AKAP9 and MCAF1 upregulate glycolysis and Wnt signaling, respectively, in advanced-stage ESCC phase. Six major tracks related to different clinical features during ESCC progression are identified, which is validated by an independent cohort with another 256 samples. Hyperphosphorylated phosphoglycerate kinase 1 (PGK1, S203) is considered as a drug target in ESCC progression. This study provides insight into the understanding of ESCC molecular mechanism and the development of therapeutic targets.

Proteogenomic characterization of cholangiocarcinoma.

BACKGROUND AND AIMS: Cholangiocarcinoma (CCA) is a highly heterogeneous cancer with limited understanding and few effective therapeutic approaches. We aimed at providing a proteogenomic CCA characterization to inform biological processes and treatment vulnerabilities. APPROACH AND RESULTS: Integrative genomic analysis with functional validation uncovered biological perturbations downstream of driver events including DPCR1 , RBM47 mutations, SH3BGRL2 copy number alterations, and FGFR2 fusions in CCA. Proteomic clustering identified three subtypes with distinct clinical outcomes, molecular features, and potential therapeutics. Phosphoproteomics characterized targetable kinases in CCA, suggesting strategies for effective treatment with CDK and MAPK inhibitors. Patients with CCA with HBV infection showed increased antigen processing and presentation (APC) and T cell infiltration, conferring a favorable prognosis compared with those without HBV infection. The characterization of extrahepatic CCA recommended the feasible application of vascular endothelial-derived growth factor inhibitors. Multiomics profiling presented distinctive molecular characteristics of the large bile duct and the small bile duct of intrahepatic CCA. The immune landscape further revealed diverse tumor immune microenvironments, suggesting immune subtypes C1 and C5 might benefit from immune checkpoint therapy. TCN1 was identified as a potential CCA prognostic biomarker, promoting cell growth by enhancing vitamin B12 metabolism. CONCLUSIONS: We characterized the proteogenomic landscape of 217 CCAs with 197 paired normal adjacent tissues and identified their subtypes and potential therapeutic targets. The multiomics analyses with other databases and some functional validations have indicated strategies regarding the clinical, biological, and therapeutic approaches to the management of CCA.

Proteogenomic insights into the biology and treatment of pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a devastating disease with poor prognosis. Proteogenomic characterization and integrative proteomic analysis provide a functional context to annotate genomic abnormalities with prognostic value. METHODS: We performed an integrated multi-omics analysis, including whole-exome sequencing, RNA-seq, proteomic, and phosphoproteomic analysis of 217 PDAC tumors with paired non-tumor adjacent tissues. In vivo functional experiments were performed to further illustrate the biological events related to PDAC tumorigenesis and progression. RESULTS: A comprehensive proteogenomic landscape revealed that TP53 mutations upregulated the CDK4-mediated cell proliferation process and led to poor prognosis in younger patients. Integrative multi-omics analysis illustrated the proteomic and phosphoproteomic alteration led by genomic alterations such as KRAS mutations and ADAM9 amplification of PDAC tumorigenesis. Proteogenomic analysis combined with in vivo experiments revealed that the higher amplification frequency of ADAM9 (8p11.22) could drive PDAC metastasis, though downregulating adhesion junction and upregulating WNT signaling pathway. Proteome-based stratification of PDAC revealed three subtypes (S-I, S-II, and S-III) related to different clinical and molecular features. Immune clustering defined a metabolic tumor subset that harbored FH amplicons led to better prognosis. Functional experiments revealed the role of FH in altering tumor glycolysis and in impacting PDAC tumor microenvironments. Experiments utilizing both in vivo and in vitro assay proved that loss of HOGA1 promoted the tumor growth via activating LARP7-CDK1 pathway. CONCLUSIONS: This proteogenomic dataset provided a valuable resource for researchers and clinicians seeking for better understanding and treatment of PDAC.

Integrated proteomic and transcriptomic landscape of macrophages in mouse tissues.

Macrophages are involved in tissue homeostasis and are critical for innate immune responses, yet distinct macrophage populations in different tissues exhibit diverse gene expression patterns and biological processes. While tissue-specific macrophage epigenomic and transcriptomic profiles have been reported, proteomes of different macrophage populations remain poorly characterized. Here we use mass spectrometry and bulk RNA sequencing to assess the proteomic and transcriptomic patterns, respectively, of 10 primary macrophage populations from seven mouse tissues, bone marrow-derived macrophages and the cell line RAW264.7. The results show distinct proteomic landscape and protein copy numbers between tissue-resident and recruited macrophages. Construction of a hierarchical regulatory network finds cell-type-specific transcription factors of macrophages serving as hubs for denoting tissue and functional identity of individual macrophage subsets. Finally, Il18 is validated to be essential in distinguishing molecular signatures and cellular function features between tissue-resident and recruited macrophages in the lung and liver. In summary, these deposited datasets and our open proteome server ( http://macrophage.mouseprotein.cn ) integrating all information will provide a valuable resource for future functional and mechanistic studies of mouse macrophages.

Integrated proteogenomic characterization across major histological types of pituitary neuroendocrine tumors.

Pituitary neuroendocrine tumor (PitNET) is one of the most common intracranial tumors. Due to its extensive tumor heterogeneity and the lack of high-quality tissues for biomarker discovery, the causative molecular mechanisms are far from being fully defined. Therefore, more studies are needed to improve the current clinicopathological classification system, and advanced treatment strategies such as targeted therapy and immunotherapy are yet to be explored. Here, we performed the largest integrative genomics, transcriptomics, proteomics, and phosphoproteomics analysis reported to date for a cohort of 200 PitNET patients. Genomics data indicate that GNAS copy number gain can serve as a reliable diagnostic marker for hyperproliferation of the PIT1 lineage. Proteomics-based classification of PitNETs identified 7 clusters, among which, tumors overexpressing epithelial-mesenchymal transition (EMT) markers clustered into a more invasive subgroup. Further analysis identified potential therapeutic targets, including CDK6, TWIST1, EGFR, and VEGFR2, for different clusters. Immune subtyping to explore the potential for application of immunotherapy in PitNET identified an association between alterations in the JAK1-STAT1-PDL1 axis and immune exhaustion, and between changes in the JAK3-STAT6-FOS/JUN axis and immune infiltration. These identified molecular markers and alternations in various clusters/subtypes were further confirmed in an independent cohort of 750 PitNET patients. This proteogenomic analysis across traditional histological boundaries improves our current understanding of PitNET pathophysiology and suggests novel therapeutic targets and strategies.

Proteomic characterization of gastric cancer response to chemotherapy and targeted therapy reveals new therapeutic strategies.

Chemotherapy and targeted therapy are the major treatments for gastric cancer (GC), but drug resistance limits its effectiveness. Here, we profile the proteome of 206 tumor tissues from patients with GC undergoing either chemotherapy or anti-HER2-based therapy. Proteome-based classification reveals four subtypes (G-I-G-IV) related to different clinical and molecular features. MSI-sig high GC patients benefit from docetaxel combination treatment, accompanied by anticancer immune response. Further study reveals patients with high T cell receptor signaling respond to anti-HER2-based therapy; while activation of extracellular matrix/PI3K-AKT pathway impair anti-tumor effect of trastuzumab. We observe CTSE functions as a cell intrinsic enhancer of chemosensitivity of docetaxel, whereas TKTL1 functions as an attenuator. Finally, we develop prognostic models with high accuracy to predict therapeutic response, further validated in an independent validation cohort. This study provides a rich resource for investigating the mechanisms and indicators of chemotherapy and targeted therapy in GC.

Proteomic analysis reveals key differences between squamous cell carcinomas and adenocarcinomas across multiple tissues.

Squamous cell carcinoma (SCC) and adenocarcinoma (AC) are two main histological subtypes of solid cancer; however, SCCs are derived from different organs with similar morphologies, and it is challenging to distinguish the origin of metastatic SCCs. Here we report a deep proteomic analysis of 333 SCCs of 17 organs and 69 ACs of 7 organs. Proteomic comparison between SCCs and ACs identifies distinguishable pivotal pathways and molecules in those pathways play consistent adverse or opposite prognostic roles in ACs and SCCs. A comparison between common and rare SCCs highlights lipid metabolism may reinforce the malignancy of rare SCCs. Proteomic clusters reveal anatomical features, and kinase-transcription factor networks indicate differential SCC characteristics, while immune subtyping reveals diverse tumor microenvironments across and within diagnoses and identified potential druggable targets. Furthermore, tumor-specific proteins provide candidates with differentially diagnostic values. This proteomics architecture represents a public resource for researchers seeking a better understanding of SCCs and ACs.

Integrated proteogenomic characterization of urothelial carcinoma of the bladder.

BACKGROUND: Urothelial carcinoma (UC) is the most common pathological type of bladder cancer, a malignant tumor. However, an integrated multi-omics analysis of the Chinese UC patient cohort is lacking. METHODS: We performed an integrated multi-omics analysis, including whole-exome sequencing, RNA-seq, proteomic, and phosphoproteomic analysis of 116 Chinese UC patients, comprising 45 non-muscle-invasive bladder cancer patients (NMIBCs) and 71 muscle-invasive bladder cancer patients (MIBCs). RESULT: Proteogenomic integration analysis indicated that SND1 and CDK5 amplifications on chromosome 7q were associated with the activation of STAT3, which was relevant to tumor proliferation. Chromosome 5p gain in NMIBC patients was a high-risk factor, through modulating actin cytoskeleton implicating in tumor cells invasion. Phosphoproteomic analysis of tumors and morphologically normal human urothelium produced UC-associated activated kinases, including CDK1 and PRKDC. Proteomic analysis identified three groups, U-I, U-II, and U-III, reflecting distinct clinical prognosis and molecular signatures. Immune subtypes of UC tumors revealed a complex immune landscape and suggested the amplification of TRAF2 related to the increased expression of PD-L1. Additionally, increased GARS, related to subtype U-II, was validated to promote pentose phosphate pathway by inhibiting activities of PGK1 and PKM2. CONCLUSIONS: This study provides a valuable resource for researchers and clinicians to further identify molecular pathogenesis and therapeutic opportunities in urothelial carcinoma of the bladder.

Trans-acting non-synonymous variant of FOXA1 predisposes to hepatocellular carcinoma through modulating FOXA1-ERα transcriptional program and may have undergone natural selection.

Interplay of pioneer transcription factor forkhead box A1 (FOXA1) and estrogen receptor has been implicated in sexual dimorphism in hepatocellular carcinoma (HCC), but etiological relevance of its polymorphism was unknown. In the case control study (1152 patients versus1242 controls), we observed significant increase in HCC susceptibility in hepatitis B virus carriers associated with a non-synonymous Thr83Ala variant of FOXA1 (odds ratio [OR], 1.28; 95% confidence interval [CI], 1.11-1.48, for Ala83-containing genotype, after validation in an independent population with 933 patients versus 1030 controls), a tightly linked (CGC)5/6or7 repeat polymorphism at its promoter (OR 1.32; 95% CI 1.10-1.60, for (CGC)6or7-repeat-containing genotype), and their combined haplotype (OR 1.50; 95% CI 1.24-1.81, for (CGC)6or7-Ala83 haplotype). The susceptible FOXA1-Ala83 impairs its interaction with ERα, attenuates transactivation toward some of their dual target genes, such as type 1 iodothyronine deiodinase, UDP glucuronosyltransferase 2 family, polypeptide B17 and sodium/taurocholate cotransporting polypeptide, but correlates with strengthened cellular expression of α-fetoprotein (AFP) and elevated AFP serum concentration in HCC patients (n = 1096). The susceptible FOXA1 cis-variant with (CGC)6or7 repeat strengthens the binding to transcription factor early growth response 1 and enhances promoter activity and gene expression. Evolutionary population genetics analyses with public datasets reveal significant population differentiation and unique haplotype structure of the derived protective FOXA1-Thr83 and suggest that it may have undergone positive natural selection in Chinese population. These findings epidemiologically highlight the functional significance of FOXA1-ERα transcriptional program and regulatory network in liver cancer development.

Proteomics-based screening of the target proteins associated with antidepressant-like effect and mechanism of Saikosaponin A.

Depression is a commonly occurring neuropsychiatric disease with an increasing incidence rate. Saikosaponin A (SA), a major bioactive component extracted from Radix Bupleuri, possesses anti-malignant cell proliferation, anti-inflammation, anti-oxidation and liver protective effects. However, few studies have investigated SA's antidepressant effects and pharmacological mechanisms of action. Our study aimed to explore the anti-depression effect of SA and screen the target proteins regulated by SA in a rat model of chronic unpredictable mild stress (CUMS)-induced depression. Results showed that 8-week CUMS combined with separation could successfully produce depressive-like behaviours and cause a decrease of dopamine (DA) in rat hippocampus, and 4-week administration of SA could relieve CUMS rats' depressive symptoms and up-regulated DA content. There were 15 kinds of significant differentially expressed proteins that were detected not only between the control and CUMS groups, but also between the CUMS and SA treatment groups. Proline-rich transmembrane protein 2 (PRRT2) was down-regulated by CUMS while up-regulated by SA. These findings reveal that SA may exert antidepressant effects by up-regulating the expression level of PRRT2 and increasing DA content in hippocampus. The identification of these 15 differentially expressed proteins, including PRRT2, provides further insight into the treatment mechanism of SA for depression.

Transketolase (TKT) activity and nuclear localization promote hepatocellular carcinoma in a metabolic and a non-metabolic manner.

BACKGROUND: Metabolic reprogramming is one of the hallmarks of cancer cells. The pentose phosphate pathway (PPP), a branch of glycolysis, is an important metabolic pathway for the survival and biosynthesis of cancer cells. Transketolase (TKT) is a key enzyme in the non-oxidative phase of PPP. The mechanistic details of TKT in hepatocellular carcinoma (HCC) development remain unclear. METHODS: TKT level and subcellular location were examined in HCC cell lines and tissue samples. We established the TKT overexpression and knocking-down stable cells in HCC cell lines. Proliferation, migration, viability and enzyme activity assays in vitro, tumor growth and metastasis assays in vivo were employed to test the effects of TKT on HCC development. GFP-tagged TKT truncations and mutants were used to locate the nuclear localization sequence (NLSs) of TKT. Cross-linking co-IP/MS was applied to identify the interaction proteins of nuclear TKT. RESULTS: We showed that TKT increased the proliferation and migration of HCC cells, as well as the viability under oxidative stress in vitro and accelerated the growth and metastasis of HCC cells in vivo. We found as a key enzyme of PPP, TKT could promote the proliferation, cell cycle, migration and viability by regulating the metabolic flux. Moreover, it was firstly reported that unlike other key enzymes in PPP, TKT showed a strong nuclear localization in HCC cells. We found not only high TKT expression, but also its nuclear localization was a prediction for poor prognosis of HCC patients. We further identified the nuclear localization sequences (NLS) for TKT and demonstrated the NLS mutations decreased the pro-tumor function of TKT independent of the enzyme activity. Cross-linking Co-IP/MS showed that nuclear TKT interacted with kinases and transcriptional coregulators such as EGFR and MAPK3, which are associated with cell activation or stress response processes. EGF treatment significantly increased the viability and proliferation of HCC cells in the enzyme-inactivating mutation TKT-D155A overexpression cells but not in the NLS-D155A double mutant group, which could be blocked by EGFR inhibitor erlotinib treatment. CONCLUSIONS: Our research suggests that in addition to the metabolic manner, TKT can promote the development of HCC in a non-metabolic manner via its nuclear localization and EGFR pathway.