Int J Biol Sci 2023; 19(5):1336-1351. doi:10.7150/ijbs.80737 This issue Cite
Research Paper
Hepatobiliary Center, The First Affiliated Hospital of Nanjing Medical University
# These authors contributed to the work equally and should be regarded as co-first authors.
Cholangiocarcinoma (CCA) is the second most common primary hepatic malignancy and associated with poor prognosis. Lack of therapeutic methods for CCA and insensitivity of targeted therapy and immunotherapy make its treatment challenging. NUF2, a component of Ndc80 kinetochore complex, is implicated in the initiation and development of multiple cancers. However, the role and mechanism of NUF2 in CCA is still unclear. In this research, we investigated the biological processes and underlying mechanisms of NUF2 in CCA. We discovered that the expression of NUF2 was upregulated in CCA and negatively correlated with prognosis. Changes in NUF2 levels had an impact on cell proliferation and migration. Moreover, NUF2 functioned as an oncogene to promote the progression of CCA through p38/MAPK signaling by inhibiting p62 binding of TFR1 and affecting its autophagic degradation. In addition, TFR1 promoted CCA progression and Kaplan-Meier analyses uncovered patients with high expression of TFR1 was associated with the poor survival. In conclusion, our study demonstrated that NUF2 promoted CCA progression by regulating TFR1 protein degradation, and the NUF2/TFR1/MAPK axis could be an excellent therapeutic target for CCA.
Keywords: Cholangiocarcinoma, NUF2, TFR1, MAPK, Autophagy
Cholangiocarcinoma (CCA) is a heterogeneous group of malignancies arising from varying locations within the biliary tree, which is associated with poor prognosis [1]. CCA is the second most frequent type of primary liver cancer, accounting for approximately 15% of all primary liver tumors and 3% of all gastrointestinal cancers [2]. Although in some countries CCA is a rare cancer, the morbidity and mortality are increasing exceptionally worldwide [3, 4]. The most common classification of CCA is intrahepatic cholangiocarcinoma (iCCA), perihilar cholangiocarcinoma (pCCA), and distant cholangiocarcinoma (dCCA) according to anatomical location, which also have different oncologic characters, clinical features, treatment strategies, and prognosis [5-7]. At present, a cure may only be achieved by radical resection or liver transplant after strict selection of early-stage disease; however, most patients present with advanced stage and lost the opportunity of surgery because of no obvious symptoms at the early stage [8-10]. Due to the scarcity of research on biomarkers and drug targets, CCA therapy is in a slow process. Hence, more effective biomarkers are urgent to be found for diagnosis and treatment [11-13].
The correct chromosome segregation during mitosis depends on the interaction between the spindle microtubules and the kinetochore, and it is well-recognized that aberrant chromosomal separation during mitosis is a common cause of carcinogenesis [14, 15]. NUF2, the component of Ndc80 kinetochore complex, is a crucial substance that maintains spindle microtubule-kinetochore connection throughout the metaphase of cell division [16, 17]. Previous studies showed that depletion of NUF2 resulted in a strong prometaphase block with an active spindle checkpoint [18]. Therefore, it is not surprising that the change of NUF2 expression leading to its dysfunction can promote tumor development. This has been verified in renal clear cell carcinoma [19], colorectal and gastric cancers [20]. However, there has been no research on the connection between NUF2 and CCA. So little is understood about the precise role, clinical implications, and underlying mechanism of NUF2 in CCA.
TFR1, the transferrin receptor, is an essential protein involved in the uptake of iron and regulation of cell growth [21]. Iron is necessary for DNA synthesis as well as other biological functions. Malignant cells frequently overexpress TFR1 due to its essential role in the pathophysiology of cancer cells, and this increased expression has been linked to a poor prognosis [22]. TFR1 is an appealing target for antibody-mediated therapy due to the high levels of TFR1 expression on malignant cells, as well as its extracellular accessibility, capacity for internalization, and crucial role in the pathophysiology of cancer cells. Although TFR1 has been used as a popular target for antibody-mediated cancer therapy over the years, interest in both targeting the receptor for delivery and using it as a direct anticancer agent is growing [23].
Autophagy is a conserved catabolic system that destroys extra or malfunctioning cell components through the lysosome [24]. It was reported that perturbations in autophagy have been found in many cancers, which have dual functions of promoting tumor and inhibiting tumor [25]. Autophagy is a lysosomal-mediated process including autophagosomes formation, autophagosome-lysosomal fusion, and autophagolysosomal degradation [26]. Autophagy consists of over 30 core autophagy-related proteins, including LC3B and p62. Double membrane vesicles are considered as a typical feature of autophagy. These vesicles could engulf the damaged proteins and organelles and fuse with lysosomal vesicles.
This study aimed to explore the functions and mechanisms of NUF2 in CCA. Our results demonstrated that high expression of NUF2 was closely related to poor prognosis, and NUF2 can act as an independent predictor for CCA. In addition, our in vitro experiments emphasized the role of NUF2/TFR1/MAPK axis in the progression of CCA, and can be used as a new target for the treatment of CCA in the future.
The human normal bile duct cell line HIBEpiC and cholangiocarcinoma cell lines, QBC939, RBE, HuCCT1 and HCCC9810 used in this study were purchased from the Chinese National Human Genome Center (Shanghai, China). Cells were cultured in DMEM culture medium supplemented with 10% fetal bovine serum (FBS)and 1% penicillin/streptomycin under 5% CO2 at 37 °C. NUF2 and TFR1 knockdown were used by siRNA. Cells were transfected with 50 nM siRNA against NUF2 (5'-GUGGCAUUAUCAACUUUAUTT-3', 5'-GCUACAACAAUCACUAAAUTT-3', 5'-CACCAAGAAUGAUCUUUAUTT-3'), TFR1 (5'-UAGAGAAUGCUGAUCUAGCUU-3', 5'-AACTTCAAGGTTTCTGCCAGC-3'), p62 (CUUCCGAAUCUACAUUAAATT) or scrambled control siRNA (Genepharma Co., Ltd, Shanghai, China) using lipofectamine 2000. The NUF2 overexpression lentiviruses and TFR1 overexpression plasmid were constructed by Cories Biotechnology Co., Ltd. (Nanjing, China). The cells were transfected according to the manufacturer's protocol. Knockdown and overexpression efficiency was evaluated by western blot.
A cohort of 38 pairs of CCA tissues and adjacent normal tissues were obtained from CCA patients who underwent curative resection in The First Affiliated Hospital of Nanjing Medicine University. Outdo Biotech Company (Shanghai, China) made the tissue microarray using samples from CCA patients who underwent surgery at The First Affiliated Hospital of Nanjing Medicine University between 2006 and 2017. The expression of NUF2 and TFR1 was scored by two independent pathologists and all differences that arose were resolved by discussion. The intensity was classified as negative (0), weak (1), moderate (2), or strong (3), and the density of positive cells in the target region was scored as follows: 0-25% (1), 26-50% (2), 51-75% (3), and 76-100% (4). The overall score was calculated by multiplying intensity and density. The expression level was categorized as low (0-4) and high (6-12). We showed representative IHC images for the scoring criteria in Figure S1. All the patients understood the procedure and provided written consent for the study.
Three pairs of QBC939/siNC and QBC939/siNUF2 were prepared for RNA-seq analysis by Personalbio (Shanghai, China). The fold change cut-off identified a significant difference in mRNA expression (P-value < 0.05 and |log2 FC| >1) between groups.
Cells were lysed with Nonidet P-40 Substitute. Protein complexes were captured using specific antibodies and IgG at 4 °C overnight and protein A+G magnetic beads were used to capture the antigen antibody complex the next day. Mixed at room temperature for two hours, discarded the supernatant, and obtained the antigen antibody complex for subsequent experiments. The protein complexes were separated by SDS-PAGE and detected with corresponding antibodies by western blot assays. Cut out candidate targets with obvious differences from the gel and send them to the company for mass spectrometry for protein identification.
Cells were plated in a 24-well plate with appropriate density and transfected with mRFP-GFP-LC3 adenovirus (Genepharma Co., Ltd, Shanghai, China) for 72h. Then we used a confocal microscope (Zeiss, LSM880) to observe cells. mRFP was utilized to label and track LC3. GFP expression was reduced, indicating that lysosomes and autophagosomes had combined to create autophagolysosomes, and GFP was quenched as a result of a pH adjustment, at which time only red fluorescence was seen. Following microscopic imaging, the red and green fluorescence pictures were combined; the yellow and red puncta were regarded as autophagosomes and autophagolysosomes, respectively. The quantity of yellow and red puncta was counted to determine the intensity of autophagy flux.
Autophagosomes were observed by transmission electron microscopy. 2.5% glutaraldehyde was used to chemically fix cells for 30 minutes at room temperature. The cells were then put into pyramid-tip moulds, dehydrated using a graduated series of ethanol, and polymerized for 72 hours at 60 °C. Sections that were cut to a thickness of 100 nm were gathered. The grids were stained for 20 minutes with 1% uranyl acetate, then left at room temperature for 10 minutes with Reynolds lead citrate. A Hitachi H-7650 transmission electron microscope (Hitachi High-Tech, Japan) operating at 80 kV was used to analyze the sections. Software called AMT-600 Image Capture Engine was used to take the pictures.
Fixed cells by 4% paraformaldehyde at 37 °C for 45 min and permeabilized by 0.3% Triton X-100 for 15 minutes. Then block cells with 5% donkey serum for 30 min and incubated them with specific primary antibodies at 4 °C overnight. We developed fluorescence by incubating with Cy3, 488 or FITC labeled goat anti-rabbit IgG for 2 h at room temperature. The nucleus was stained with DAPI. Fluorescence image were captured by a confocal fluorescence microscope.
All experiments were repeated at least three times. The results were presented as the mean ± standard error of the mean, unless stated otherwise. SPSS v24.0 (IBM, SPSS, Chicago, IL, USA) and Graphpad Prism 9 (GraphPad Software, La Jolla, CA, USA) were used to perform statistical analyses. Differences between the two groups were analyzed by Student's t-test. Moreover, Overall survival and disease-free survival were analyzed with Kaplan-Meier methods, and log-rank test was applied for comparison. A multivariate analysis was performed using Cox proportional hazard regression model. A p-value of <0.05 was considered significant: *p < 0.05, **p < 0.01, and ***p < 0.001.
More materials and methods were recorded in supplementary file.
Our previous RNA-seq found that NUF2 was upregulated in CCA tissues. The pan-cancer analysis of NUF2 in TCGA database showed that NUF2 was highly expressed in most tumors, including in CCA (Figure 1A). We next examined 38 pairs of CCA tissues in our center, and found that the mRNA expression of NUF2 was upregulated in CCA tissues compared to matched normal tissue samples (Figure 1B). Western blot also showed that the protein level of NUF2 was significantly increased in tumor tissues (Figure 1C). To analyze the relationship between NUF2 expression and prognosis, we performed tissue microarrays by IHC staining in 134 CCA patients. The results showed that NUF2 expression was higher in tumor tissues than matched normal bile duct tissues (Figure 1D-E). The ratio of lymphatic metastasis, tumor size>5cm, and vascular invasion was significantly higher in NUF2 high expression group than that in low expression group (Figure 1F-H). Kaplan-Meier survival analysis showed that patients with high expression of NUF2 had a shorter overall survival (OS) and disease-free survival (DFS) compared with low- NUF2 CCA patients (Figure 1I-J). These results indicated that NUF2 might be a critical regulator in CCA progression and functioned as an independent predictive factor of prognosis for CCA patients.
NUF2 was up regulated in CCA and predicated the poor prognosis. (A) Pan-cancer analysis of NUF2 and analysis of NUF2 mRNA level in CCA tissues and adjacent tissues using TCGA database. (B) qRT-PCR analysis of NUF2 expression in CCA tissues and adjacent tissues. (C) Western blot analysis of NUF2 expression in CCA tissues and adjacent tissues. (D) Immunohistochemical staining images of NUF2 in CCA tissues and adjacent tissues (Scale bar: 200 μm, 50 μm). (E) Statistical analysis of NUF2 expression according to the IHC assay. (F-H) Ratio of lymphatic metastasis, tumor size>5cm, and vascular invasion in high NUF2 group and low NUF2 group. (I) Overall survival and (J) Disease-free survival related to NUF2 expression by Kaplan-Meier survival curve analysis.
NUF2 promoted CCA cell proliferation and migration in vitro. (A) CCK8 assays, (B) clone formation, (C) EdU assays were performed to identify the proliferation ability after NUF2 knockdown or overexpression. (D) Wound healing and (E) transwell assays were performed to identify the migration ability after NUF2 knockdown or overexpression. (F) Flow cytometry analysis of the cell cycle of CCA cells. (G) The expression of N-cad and Vimentin was evaluated by IF staining. (H) Western blot analysis of the levels of vimentin, N-cadherin, E-cadherin, and snail in HuCCT1 (upper) and QBC939 (lower) cells.
We then explored the role of NUF2 in CCA progression in vitro by evaluating cell proliferation and migration. Compared with normal biliary epithelial cells, the expression of NUF2 was significantly higher in human CCA cell lines QBC939, RBE and HCCC-9810 (Figure S2A). We next chose QBC939 with the highest expression to knock down, and HuCCT1 with the lowest expression to overexpress. The transfected efficiency was validated by western blot (Figure S2B).
CCK8 assays results showed that knocking down NUF2 significantly inhibited the cells proliferation of CCA, while overexpression of NUF2 promoted cells proliferation (Figure 2A and Figure S3A). The clone formation assays demonstrated that NUF2 knockdown weakened the ability of clone formation, whilst NUF2 overexpression had the opposite situation (Figure 2B and Figure S3B). With EdU assays, the percentage of red cells which means that the cells in the proliferation phase decreased significantly after knocking down NUF2, while overexpression of NUF2 significantly increased the percentage of red cells (Figure 2C and Figure S3C). Furthermore, the results of wound healing and transwell assays showed that depletion of NUF2 notably attenuated the migration ability of QBC939 cells while overexpression of NUF2 significantly promoted the migration ability of HuCCT1 (Figure 2D-E and Figure S3D-E). Flow cytometry analysis showed that NUF2 knockdown significantly increased the proportion of cells in the G0/G1 phase and decreased the proportion of cells in the S and G2/M phase. And the percentage of cells in the S and G2/M phase was significantly higher in the NUF2 overexpression group, while the percentage of cells in theG0/G1 phase decreased significantly (Figure 2F and Figure S3F). Moreover, we found that knockdown of NUF2 promoted the expression of E-cadherin and inhibited the expression of N-cadherin, snail, vimentin (Figure 2G-H and Figure S3G). Overexpression of NUF2 had the opposite phenomenon. All these results indicated that NUF2 played an essential role in the cell proliferation, migration, and EMT of CCA.
To unlock the mechanism how NUF2 regulated CCA progression, we compared the transcriptome profiles of QBC939 cells with NUF2 knockdown and controls by RNA-seq analysis. A total of 361 differential genes were found, including 130 upregulated genes and 231 downregulated genes (Figure 3A). In KEGG analysis, the top 20 signaling pathways according to the p-value were represented in Figure 3B. The MAPK pathway contained the most differential genes and was the only one closely related to cell proliferation and migration in the RNA-seq result, suggesting that this pathway may be important in NUF2 regulating CCA progression. Further research was done to determine which of the three main MAPK signaling subgroups was responsible for the carcinogenic function of NUF2. Compared to control cells, all three subgroups were downregulated in NUF2 knockdown cells, especially in p-p38 and p-ERK (Figure 3C). Overexpression of NUF2 promoted the activations of p-p38 and p-ERK.
To understand the importance of the p-p38 or p-ERK in NUF2-mediated CCA progression, we performed rescue experiments using the p38 inhibitor SB203580 and ERK inhibitor U0126 (Figure S4). Regretfully, U0126 could not inhibit the tumor-promoting effect caused by NUF2 (Figure S5). However, SB203580 restored the CCA cells proliferation and migration caused by NUF2 overexpression. The CCK8, EdU assays and clone formation demonstrated that compared with the DMSO groups, the addition of p38 inhibitor significantly reduced the CCA cells proliferation (Figure 3D-F). Wound healing and transwell assays showed that the migration ability was weakened after adding p38 inhibitor (Figure 3G-H). These results implied that NUF2 promoted CCA progression maybe through regulating MAPK signaling pathway.
To further identify the potential mechanism of NUF2 regulating MAPK, a combined IP/MS approach was used to identify novel NUF2-mediated protein- protein interactions (Figure 4A). We found that TFR1 had the highest score among the differential proteins closely related to MAPK pathway in comparison of the CO-IP products by NUF2 antibody and IgG (Figure 4B). To confirm the MS data, we performed the CO-IP assay to verify the binding between two proteins (Figure 4C). Furthermore, confocal immunofluorescence microscopy captured a high degree of spatial concordance between endogenous NUF2 and TFR1 in QBC939 cells (Figure 4D). These findings demonstrated that NUF2 could interact with TFR1. Then we investigated whether the change of NUF2 protein expression could synchronously cause the change of TFR1 expression. The results showed that the protein level of TFR1 was upregulated in NUF2 overexpression group, while was down- regulated in NUF2 knockout group (Figure 4E). However, there was no obvious change on the mRNA level (Figure 4F). We then used public database to verify our experiments and the public data showed that NUF2 had no correlation with TFR1 at mRNA level (Figure S6). Furthermore, we also verified the changes of three MAPK subgroups after TFR1 knockdown or overexpression. We found that the difference of p-p38 in three subgroups was the most obvious, and no significant change in p-ERK and p-JNK (Figure 4G). Chirasani et al. showed that TFR1 can activate p38/MAPK pathway, thereby promote tumor progression by increasing intracellular ROS concentration [27]. In our results, overexpression NUF2 promoted accumulation of intracellular ROS, while knocking down NUF2 decreased the level of intracellular ROS (Figure 4H). These results indicated that NUF2-TFR1-p38 axis maybe play a major role in regulating CCA progression.
To explored the function of TFR1 in CCA, we knocked down and over expressed TFR1 in CCA cells. The efficiency was examined by western blot (Figure S7). We used CCK8, clone formation, EdU assays, wound healing and transwell assays to explore the capacity of proliferation and migration. The results showed that TFR1 knockdown impaired the CCA cells' capacity for proliferation and migration, while TFR1 overexpressing promoted the CCA cells proliferation and migration (Figure 5A-E).
Western blot showed that the protein level of NUF2 and TFR1 was significantly increased in tumor tissues (Figure 6A). Immunohistochemical staining of CCA tissues showed that TFR1 was upregulated in CCA tissues (Figure 6B) and it was positively correlated with NUF2 expression (r = 0.2080, p < 0.05) (Figure 6C). In addition, high expression of TFR1 was associated with lymph node metastasis, tumor size, and vascular invasion. People with high TFR1 expression had worse performance in these clinical indicators (Figure 6D-F). Kaplan-Meier analysis showed that CCA patients with high expression of TFR1 had shorter overall survival and disease-free survival (Figure 6G-H). The prognosis of people with both high expression of NUF2 and TFR1 was the worst, and people with both low expressions had the best prognosis (Figure 6I-J). These results indicated that TFR1 might be a critical regulator in CCA progression.
Next, we further explored how NUF2 regulated TFR1 protein expression. As we all know, there are two main ways of protein degradation: the ubiquitination-proteasome degradation pathway and lysosomal degradation pathway [28]. MG132 and chloroquine (CQ) are two mainstream inhibitors, which inhibit ubiquitination-mediated degradation and autophagy-regulated degradation respectively. Hence, in order to find out which way played the major role, QBC939 and HuCCT1cells were treated with MG132 or CQ. And we found that CQ restored the expression of TFR1 protein but not MG132 (Figure 7A-B). The results suggested that NUF2 regulated TFR1 degradation maybe through the autolysosome pathway.
To prove our hypothesis, we examined the expression of LC3B by confocal fluorescence microscopy. The results showed that knockdown of NUF2 induced the accumulation of LC3B while overexpression of NUF2 decreased the expression of LC3B (Figure 7C). Then after GFP-mRFP-LC3 lentivirus transfected in QBC939 and HuCCT1 cells, there were more yellow puncta which indicated an increased autophagic flux in siNUF2 group compared to siNC group while a significant decrease when overexpressing NUF2 (Figure 7D). Furthermore, we found that NUF2 knockdown cells contain abundant autophagic vacuoles while NUF2 overexpression was reversely viewed by transmission electron microscopy (TEM) (Figure 7E). We next measured the protein expression level of LC3B and p62 which are the most important autophagy-related markers. NUF2 overexpression potently downregulated LC3B protein level but increased p62 protein level in HuCCT1 while NUF2 knockdown reversed the situation (Figure 7F). These results implied NUF2 inhibited the autophagic flux in CCA cells.
Substantial pieces of evidence suggested that cargo receptors are essential for delivering substrates to the autophagosome for selective degradation [29]. Although p62 is the most well-known autophagy cargo receptor in autophagy, there are still some other cargo proteins that act as transporters. Next, we carried out a series of IP experiments on common autophagy cargo receptors to find out which was truly combined with TFR1. The results demonstrated that only p62 and TFR1 had obvious binding rather than other cargo receptors NBR1, NDP52, TOLLIP and OPTN (Figure 7H-I), and immunofluorescence experiment also verified the binding (Figure 7G). Next, we hypothesized that NUF2 affected autophagy by regulating TFR1 and NUF2 binding. We transfected the siNUF2 and siNC into QBC939 cells, added with autophagy inhibitor CQ as well as. Then TFR1 antibody were used for IP. As a result, we found that NUF2 knockdown group combined with more p62 in cells treated with CQ and overexpression of NUF2 led to consistent conclusions (Figure 7J). Importantly, si-p62 could rescue TFR1 protein level in NUF2 knockdown CCA cells (Figure 7K). These results indicated that NUF2 inhibited autophagy contributed to the accumulation of TFR1 in a p62-dependent manner.
NUF2 promoted CCA cell proliferation and migration via the p38/MAPK signaling pathway. (A) Hot map of 361 differential genes in the RNA-seq. (B) Top 20 signaling pathway maps enriched in the RNA Sequence. (C) Western blot analysis of NUF2, p-JNK, p-p38, p-ERK, JNK, p38, and ERK protein levels. (D) CCK8 assays, (E) EdU assays, and (F) clone formation showed SB203580 reverses the promotional effect of NUF2 overexpression on the proliferation of CCA cells. (G) Wound healing and (H) Transwell assays showed SB203580 reverses the promotional effect of NUF2 overexpression on the migration of CCA cells.
NUF2 interacted and co-localized with TFR1. (A) NUF2-associated proteins analyzed by MS. (B) Mass spectrogram of TFR1. (C) Co-IP analysis of NUF2 and TFR1. (D) Confocal immunofluorescence analysis of NUF2 and TFR1 expression. (E) Changes of TFR1 protein level after NUF2 knockdown or overexpression. (F) Changes of TFR1 mRNA level after NUF2 knockdown or overexpression. (G) Western blot analysis of TFR1, p-JNK, p-p38, p-ERK, JNK, p38, and ERK protein levels. (H) Intracellular ROS levels were evaluated.
TFR1 promoted CCA cell proliferation and migration in vitro. (A) CCK8 assays, (B) clone formation, and (C) EdU assays were performed to identify the proliferation ability. (D) Wound healing and (E) transwell assays were performed to identify the migration ability.
Next, we continued to verified the role of TFR1 in NUF2 regulating CCA progression. The clone formation, CCK8 and EdU assays were used to explore the proliferation ability of CCA cells (Figure 8A-C). Wound healing and transwell assays were performed to detect the effect on cell migration (Figure 8D-E). As expected, overexpression of TFR1 significantly reversed the tumor inhibitory effect of NUF2 knockdown on cell proliferation and migration, which suggested that TFR1 was involved in NUF2 induced CCA malignancy, at least partly through the regulation of TFR1. In addition, a xenograft model was used to investigate the role of NUF2/TRF1/p38 pathway in CCA in vivo. The result showed that overexpression of NUF2 promoted the tumor growth. However, the inhibition of TRF1 and p38 impaired the CCA growth caused by NUF2 overexpression (Figure 8F-H). In addition, IHC staining of xenografts showed that tissues from NUF2 overexpression group showed an increase Ki67 expression (Figure 8I). These results implied that the NUF2-TFR1-p38 axis played important roles in CCA progression.
Clinical significance of TFR1. (A) Western blot analysis of NUF2 and TFR1expression in CCA tissues and adjacent tissues. (B) Immunohistochemical staining images of TFR1 in CCA tissues and adjacent tissues. (C) Correlation analysis of NUF2 and TFR1 expression according to IHC staining. (D-F) Ratio of lymphatic metastasis, tumor size>5cm, and vascular invasion in high TFR1 group and low TFR1 group. (G) Overall survival and (H) Disease-free survival related to TFR1 expression by Kaplan-Meier survival curve analysis. (I) Overall survival and (J) Disease-free survival were further analyzed by dividing patients into subgroups with the double high expression, single high expression, and dual low expression of NUF2 and TFR1.
NUF2 is an oncogene that reported in clear cell renal cell carcinoma [30], melanoma [31], hepatocellular carcinoma [32], and breast cancer [33]. More importantly, in a phase I clinical trial, NUF2 peptide vaccination has been used to treat castration resistant prostate cancer (CRPC), and achieved encouraging results [34]. Together, these findings showed that NUF2 played a unique function in the development and growth of tumors, making it a viable and efficient target for molecule-targeted therapy in a variety of malignancies. Regretfully, there was no study on the role of NUF2 in CCA. In our research, genomics and proteomics data implied that the mRNA and protein abundances of NUF2 were high in the tumor tissues of CCA, and patients with higher NUF2 expression had a poorer prognosis. Univariate and multivariate Cox regression analyses revealed that NUF2 was an independent predictive factor for the prognosis of CCA patients and this was consistent with the results of others' researches [32, 35, 36]. Further results showed that NUF2 promoted the proliferation and migration of CCA cells in vitro, as well as in vivo.
Up to now, most studies have shown that NUF2 mainly affected tumorigenesis by regulating cell cycle and DNA replication [36, 37], which was consistent with the function of NUF2 itself. Nevertheless, the direct targets that mediated NUF2's biological actions in malignancies, particularly in CCA, were not well understood. In this study, we performed RNAseq and mass spectrometry analysis and identified TFR1 and MAPK which served as cancer promoting factors to promote the development of various tumors [38-42]. Therefore, we first explored the role of NUF2/TFR1/MAPK axis in CCA. Our rescue experiments proved that NUF2 played an important role in the proliferation and migration of CCA cells through TFR1/MAPK. Our results showed that NUF2 could interact with TFR1, and control the downstream pathway by changing its protein expression. Compared with the control group, after NUF2 knockdown, TFR1 protein expression decreased, cell proliferation and migration ability weakened. Mechanistically, NUF2 inhibited the autophagic degradation of TFR1 by preventing the binding of TFR1 and p62, leading to the accumulation of TFR1 in cells and promoting the occurrence and development of tumors. As a result, our research advanced knowledge of the precise molecular processes underpinning NUF2-mediated CCA development.
Knockdown of NUF2 promoted autophagy. (A) Cells with NUF2 knockdown or overexpression were treated with CQ (50μM). (B) Cells with NUF2 knockdown or overexpression were treated with MG-132 (10μM). (C) Immunofluorescence confirmed the levels of LC3B involved in autophagy NUF2 knockdown or overexpression in HuCCT1 (upper) and QBC939 (lower) cells. (D) mCherry-GFP-LC3 analysis showed Autophagic flux in CCA cells. (E) Typical transmission electron microscopic images of endogenic autophagic double membrane structures in CCA cells. (F) The protein levels of LC3B and p62 were detected by western blot analysis. (G) Confocal immunofluorescence analysis of TFR1 and p62 expression. (H) IP analysis of the binding of TFR1 and cargo receptors. (I) IP analysis of p62 and TFR1. (J) CO-IP demonstrated the combining ability of TFR1 and p62 in NUF2 knockdown or overexpression cells. (K) Western blot showed that si-p62 can rescue TFR1 protein level in NUF2 knockdown cells.
TFR1 reversed the effect of NUF2 on CCA cells. (A) EdU assays, (B) CCK8 assays, and (C) Clone formation showed TFR1 overexpression reversed the inhibited effect of NUF2 knockdown on the proliferation of CCA cells. (D) Wound healing and (E) transwell assays showed TFR1 overexpression reversed the inhibited effect of NUF2 knockdown on the migration of CCA cells. (F-G) Volume and weight of xenograft tumor (n = 6). (H) Typical images of subcutaneous xenograft tumor. (I) Ki67, NUF2, LC3B, and p-p38 expression in different groups of xenografts.
Schematic representation showing the proposed mechanism of NUF2.
TFR1 is a type II transmembrane glycoprotein which regulates iron import and is highly expressed in proliferating exuberantly cells, such as tumor cells [40, 43, 44]. TFR1 is abnormally expressed in many human cancers. As a cancer promoting protein, it was reported to induce cancer cell proliferation by protecting cancer cells against natural killer (NK) cells [45] and inhibiting apoptosis [46], as well as regulating NF-KB [47], MAPK [48], AKT [49], and JAK/STAT pathways [50]. These results explain why NUF2 could promote the proliferation and migration of CCA via TFR1 in our study. In addition, TFR1 functioned as a universal prognosis predictor [51]. Increased expression of TFR1 correlated with advanced stage and poorer prognosis in multiple cancers [52, 53]. It was reported that TFR1 could be degraded by ubiquitin proteasome [54] or ubiquitinated and degraded in lysosome [55, 56]. According to our experiments, we found that TFR1 was mainly degraded by autophagic lysosome pathway in CCA. The autophagy-lysosomal system is crucial for intracellular protein degradation, nutrition cycling, scavenging, and stability maintenance [57]. Although TFR1 has been reported to be degraded by autophagy, the specific mechanism was still unknown. We studied for the first time how NUF2 played a role in the autophagic degradation of TFR1 through p62, and uncovered the molecular mechanism of NUF2 promoting the progression of cholangiocarcinoma through TFR1.
Overactivation of MAPK signaling is a common event in multiple cancers. MAPK pathway contains three major subtypes, including ERK, JNK, and p38. Among them, p38 kinase plays a key role in the occurrence and development of CCA by promoting proliferation, invasion, inflammation and angiogenesis [58, 59]. In our research, we found that NUF2 activated p38/MAPK by increasing TFR1 protein expression. ERK and JNK did not seem to play an obvious role in the TFR1/MAPK pathway. This was consistent with other's work [60]. As a transmembrane glycoprotein involved in Iron transport, the essential roles of TFR1 in various biological events, such as oxygen delivery, electron transport, and enzymatic reactions, depend on its redox activity [61-63]. TFR1 can promoted mitochondrial respiration and ROS production, both of which are key players in tumor cell growth and survival [64]. In addition, p-p38 is a major ROS target and ROS can affect expression of p-p38 in cholangiocarcinoma [65]. It has been reported that TFR1 can activate p38/MAPK pathway, thereby promote tumor progression by increasing intracellular ROS concentration [27].
In general, we found a new effect of NUF2 on the proliferation and metastasis of CCA cells, revealing how NUF2-TFR1 interaction impacted the stability of TFR1 to activate p38/MAPK signaling, which can work as an ideal therapeutic target in future. Consequently, targeting NUF2 or NUF2-TFR1 interaction may provide a promising therapeutic approach for CCA.
CCA: cholangiocarcinoma; HCC: hepatocellular carcinoma; iCCA: intrahepatic cholangiocarcinoma; pCCA: perihilar cholangiocarcinoma; dCCA: distant cholangiocarcinoma; TMA: tissue microarray; (RT)PCR: real-time reverse-transcription; siRNA: small interfering RNA; shRNA: short hairpin RNA; IP: immunoprecipitation; TEM: transmission electron microscopy; IF: Immunofluorescence; EMT: epithelial mesenchymal transformation; TCGA: The Cancer Genome Atlas; MAPK: mitogen activated protein kinase; IHC: immunochemistry staining.
Supplementary methods, figures and table.
This work was supported by Jiangsu Provincial Key Research and Development Program (BE2016789).
The studies were reviewed and approved by the Ethics Committee of The First Affiliated Hospital of Nanjing Medical University.
JS, WJ, JC, and CL designed the study. JS performed most of the experiments and analyzed data. TZ, XX, YW, RC, XS, SF and SL assisted in completing some experiments. JS produced the initial draft of the manuscript. XL, CL and YZ supervised the research and review the manuscript. YC, YW, and JW provided human CCA tissue. XL obtained funding for research.
The authors have declared that no competing interest exists.
1. Banales JM, Marin JJG, Lamarca A, Rodrigues PM, Khan SA, Roberts LR. et al. Cholangiocarcinoma 2020: the next horizon in mechanisms and management. Nat Rev Gastroenterol Hepatol. 2020;17:557-88
2. Banales JM, Cardinale V, Carpino G, Marzioni M, Andersen JB, Invernizzi P. et al. Expert consensus document: Cholangiocarcinoma: current knowledge and future perspectives consensus statement from the European Network for the Study of Cholangiocarcinoma (ENS-CCA). Nat Rev Gastroenterol Hepatol. 2016;13:261-80
3. Fabris L, Andersen JB, Fouassier L. Intrahepatic cholangiocarcinoma: A single-cell resolution unraveling the complexity of the tumor microenvironment. J Hepatol. 2020;73:1007-9
4. Rizvi S, Khan SA, Hallemeier CL, Kelley RK, Gores GJ. Cholangiocarcinoma - evolving concepts and therapeutic strategies. Nature reviews Clinical oncology. 2018 15
5. Chen T, Li K, Liu Z, Liu J, Wang Y, Sun R. et al. WDR5 facilitates EMT and metastasis of CCA by increasing HIF-1α accumulation in Myc-dependent and independent pathways. Mol Ther. 2021;29:2134-50
6. Li Z, Liu J, Chen T, Sun R, Liu Z, Qiu B. et al. HMGA1-TRIP13 axis promotes stemness and epithelial mesenchymal transition of perihilar cholangiocarcinoma in a positive feedback loop dependent on c-Myc. Journal of Experimental & Clinical Cancer Research: CR. 2021;40:86
7. Liu Z, Sun R, Zhang X, Qiu B, Chen T, Li Z. et al. Transcription factor 7 promotes the progression of perihilar cholangiocarcinoma by inducing the transcription of c-Myc and FOS-like antigen 1. EBioMedicine. 2019;45:181-91
8. Primrose JN, Fox RP, Palmer DH, Malik HZ, Prasad R, Mirza D. et al. Capecitabine compared with observation in resected biliary tract cancer (BILCAP): a randomised, controlled, multicentre, phase 3 study. Lancet Oncol. 2019;20:663-73
9. Abou-Alfa GK, Macarulla T, Javle MM, Kelley RK, Lubner SJ, Adeva J. et al. Ivosidenib in IDH1-mutant, chemotherapy-refractory cholangiocarcinoma (ClarIDHy): a multicentre, randomised, double-blind, placebo-controlled, phase 3 study. Lancet Oncol. 2020;21:796-807
10. Abou-Alfa GK, Sahai V, Hollebecque A, Vaccaro G, Melisi D, Al-Rajabi R. et al. Pemigatinib for previously treated, locally advanced or metastatic cholangiocarcinoma: a multicentre, open-label, phase 2 study. Lancet Oncol. 2020;21:671-84
11. Diggs LP, Ruf B, Ma C, Heinrich B, Cui L, Zhang Q. et al. CD40-mediated immune cell activation enhances response to anti-PD-1 in murine intrahepatic cholangiocarcinoma. J Hepatol. 2021;74:1145-54
12. O'Rourke CJ, Munoz-Garrido P, Andersen JB. Molecular Targets in Cholangiocarcinoma. Hepatology (Baltimore, Md). 2021;73(Suppl 1):62-74
13. Silverman IM, Hollebecque A, Friboulet L, Owens S, Newton RC, Zhen H. et al. Clinicogenomic Analysis of -Rearranged Cholangiocarcinoma Identifies Correlates of Response and Mechanisms of Resistance to Pemigatinib. Cancer Discov. 2021;11:326-39
14. Cheerambathur DK, Gassmann R, Cook B, Oegema K, Desai A. Crosstalk between microtubule attachment complexes ensures accurate chromosome segregation. Science. 2013;342:1239-42
15. Foley EA, Kapoor TM. Microtubule attachment and spindle assembly checkpoint signalling at the kinetochore. Nat Rev Mol Cell Biol. 2013;14:25-37
16. Nabetani A, Koujin T, Tsutsumi C, Haraguchi T, Hiraoka Y. A conserved protein, Nuf2, is implicated in connecting the centromere to the spindle during chromosome segregation: a link between the kinetochore function and the spindle checkpoint. Chromosoma. 2001;110:322-34
17. DeLuca JG, Howell BJ, Canman JC, Hickey JM, Fang G, Salmon ED. Nuf2 and Hec1 are required for retention of the checkpoint proteins Mad1 and Mad2 to kinetochores. Curr Biol. 2003;13:2103-9
18. DeLuca JG, Moree B, Hickey JM, Kilmartin JV, Salmon ED. hNuf2 inhibition blocks stable kinetochore-microtubule attachment and induces mitotic cell death in HeLa cells. J Cell Biol. 2002;159:549-55
19. Xie X, Lin J, Fan X, Zhong Y, Chen Y, Liu K. et al. LncRNA CDKN2B-AS1 stabilized by IGF2BP3 drives the malignancy of renal clear cell carcinoma through epigenetically activating NUF2 transcription. Cell Death & Disease. 2021;12:201
20. Kaneko N, Miura K, Gu Z, Karasawa H, Ohnuma S, Sasaki H. et al. siRNA-mediated knockdown against CDCA1 and KNTC2, both frequently overexpressed in colorectal and gastric cancers, suppresses cell proliferation and induces apoptosis. Biochem Biophys Res Commun. 2009;390:1235-40
21. Neckers LM, Trepel JB. Transferrin receptor expression and the control of cell growth. Cancer Invest. 1986;4:461-70
22. Fan K, Cao C, Pan Y, Lu D, Yang D, Feng J. et al. Magnetoferritin nanoparticles for targeting and visualizing tumour tissues. Nat Nanotechnol. 2012;7:459-64
23. Candelaria PV, Leoh LS, Penichet ML, Daniels-Wells TR. Antibodies Targeting the Transferrin Receptor 1 (TfR1) as Direct Anti-cancer Agents. Frontiers In Immunology. 2021;12:607692
24. Mukhopadhyay S, Panda PK, Sinha N, Das DN, Bhutia SK. Autophagy and apoptosis: where do they meet? Apoptosis. 2014;19:555-66
25. Cao Y, Luo Y, Zou J, Ouyang J, Cai Z, Zeng X. et al. Autophagy and its role in gastric cancer. Clin Chim Acta. 2019;489:10-20
26. Andrejeva G, Gowan S, Lin G, Wong Te Fong A-CL, Shamsaei E, Parkes HG. et al. phosphatidylcholine synthesis is required for autophagosome membrane formation and maintenance during autophagy. Autophagy. 2020;16:1044-60
27. Chirasani SR, Markovic DS, Synowitz M, Eichler SA, Wisniewski P, Kaminska B. et al. Transferrin-receptor-mediated iron accumulation controls proliferation and glutamate release in glioma cells. J Mol Med (Berl). 2009;87:153-67
28. Dikic I. Proteasomal and Autophagic Degradation Systems. Annu Rev Biochem. 2017;86:193-224
29. Liu J, Wu X, Wang H, Wei J, Wu Q, Wang X. et al. HFE inhibits type I IFNs signaling by targeting the SQSTM1-mediated MAVS autophagic degradation. Autophagy. 2021;17:1962-77
30. Lin J, Chen X, Yu H, Min S, Chen Y, Li Z. et al. NUF2 Drives Clear Cell Renal Cell Carcinoma by Activating HMGA2 Transcription through KDM2A-mediated H3K36me2 Demethylation. International Journal of Biological Sciences. 2022;18:3621-35
31. Wang D, Chen J, Li B, Jiang Q, Liu L, Xia Z. et al. A noncoding regulatory RNA Gm31932 induces cell cycle arrest and differentiation in melanoma via the miR-344d-3-5p/Prc1 (and Nuf2) axis. Cell Death & Disease. 2022;13:314
32. Wang Y, Tan PY, Handoko YA, Sekar K, Shi M, Xie C. et al. NUF2 is a valuable prognostic biomarker to predict early recurrence of hepatocellular carcinoma after surgical resection. International journal of cancer. 2019;145:662-70
33. Lv S, Xu W, Zhang Y, Zhang J, Dong X. NUF2 as an anticancer therapeutic target and prognostic factor in breast cancer. International journal of oncology. 2020;57:1358-67
34. Obara W, Sato F, Takeda K, Kato R, Kato Y, Kanehira M. et al. Phase I clinical trial of cell division associated 1 (CDCA1) peptide vaccination for castration resistant prostate cancer. Cancer Sci. 2017;108:1452-7
35. Chen M, Li S, Liang Y, Zhang Y, Luo D, Wang W. Integrative Multi-Omics Analysis of Identified NUF2 as a Candidate Oncogene Correlates with Poor Prognosis and Immune Infiltration in Non-Small Cell Lung Cancer. Front Oncol. 2021;11:656509
36. Jiang X, Jiang Y, Luo S, Sekar K, Koh CKT, Deivasigamani A. et al. Correlation of NUF2 Overexpression with Poorer Patient Survival in Multiple Cancers. Cancer Res Treat. 2021;53:944-61
37. Hu P, Chen X, Sun J, Bie P, Zhang L-D. siRNA-mediated knockdown against NUF2 suppresses pancreatic cancer proliferation in vitro and in vivo. Bioscience Reports. 2015 35
38. Corte-Rodríguez M, Blanco-González E, Bettmer J, Montes-Bayón M. Quantitative Analysis of Transferrin Receptor 1 (TfR1) in Individual Breast Cancer Cells by Means of Labeled Antibodies and Elemental (ICP-MS) Detection. Anal Chem. 2019;91:15532-8
39. Arumov A, Liyanage PY, Trabolsi A, Roberts ER, Li L, Ferreira BCLB. et al. Optimized Doxorubicin Chemotherapy for Diffuse Large B-cell Lymphoma Exploits Nanocarrier Delivery to Transferrin Receptors. Cancer Research. 2021;81:763-75
40. Wagner EF, Nebreda AR. Signal integration by JNK and p38 MAPK pathways in cancer development. Nat Rev Cancer. 2009;9:537-49
41. Zang WJ, Hu YL, Qian CY, Feng Y, Liu JZ, Yang JL. et al. HDAC4 promotes the growth and metastasis of gastric cancer via autophagic degradation of MEKK3. Br J Cancer. 2022;127:237-48
42. Ge QY, Chen J, Li GX, Tan XL, Song J, Ning D. et al. GRAMD4 inhibits tumour metastasis by recruiting the E3 ligase ITCH to target TAK1 for degradation in hepatocellular carcinoma. Clin Transl Med. 2021;11:e635
43. Kawabata H. Transferrin and transferrin receptors update. Free Radic Biol Med. 2019;133:46-54
44. Gammella E, Buratti P, Cairo G, Recalcati S. The transferrin receptor: the cellular iron gate. Metallomics. 2017;9:1367-75
45. Jiang X-P, Elliott RL. Decreased Iron in Cancer Cells and Their Microenvironment Improves Cytolysis of Breast Cancer Cells by Natural Killer Cells. Anticancer Res. 2017;37:2297-305
46. Pham CG, Bubici C, Zazzeroni F, Papa S, Jones J, Alvarez K. et al. Ferritin heavy chain upregulation by NF-kappaB inhibits TNFalpha-induced apoptosis by suppressing reactive oxygen species. Cell. 2004;119:529-42
47. Kenneth NS, Mudie S, Naron S, Rocha S. TfR1 interacts with the IKK complex and is involved in IKK-NF-κB signalling. Biochem J. 2013;449:275-84
48. Senyilmaz D, Virtue S, Xu X, Tan CY, Griffin JL, Miller AK. et al. Regulation of mitochondrial morphology and function by stearoylation of TFR1. Nature. 2015;525:124-8
49. Cheng X, Fan K, Wang L, Ying X, Sanders AJ, Guo T. et al. TfR1 binding with H-ferritin nanocarrier achieves prognostic diagnosis and enhances the therapeutic efficacy in clinical gastric cancer. Cell Death & Disease. 2020;11:92
50. Cui C, Cheng X, Yan L, Ding H, Guan X, Zhang W. et al. Downregulation of TfR1 promotes progression of colorectal cancer via the JAK/STAT pathway. Cancer Manag Res. 2019;11:6323-41
51. Essaghir A, Demoulin J-B. A minimal connected network of transcription factors regulated in human tumors and its application to the quest for universal cancer biomarkers. PloS One. 2012;7:e39666
52. Chan KT, Choi MY, Lai KKY, Tan W, Tung LN, Lam HY. et al. Overexpression of transferrin receptor CD71 and its tumorigenic properties in esophageal squamous cell carcinoma. Oncol Rep. 2014;31:1296-304
53. Habashy HO, Powe DG, Staka CM, Rakha EA, Ball G, Green AR. et al. Transferrin receptor (CD71) is a marker of poor prognosis in breast cancer and can predict response to tamoxifen. Breast Cancer Res Treat. 2010;119:283-93
54. Wu Y, Jiao H, Yue Y, He K, Jin Y, Zhang J. et al. Ubiquitin ligase E3 HUWE1/MULE targets transferrin receptor for degradation and suppresses ferroptosis in acute liver injury. Cell Death Differ. 2022
55. Fujita H, Iwabu Y, Tokunaga K, Tanaka Y. Membrane-associated RING-CH (MARCH) 8 mediates the ubiquitination and lysosomal degradation of the transferrin receptor. J Cell Sci. 2013;126:2798-809
56. Xiong Q, Li X, Li W, Chen G, Xiao H, Li P. et al. WDR45 Mutation Impairs the Autophagic Degradation of Transferrin Receptor and Promotes Ferroptosis. Front Mol Biosci. 2021;8:645831
57. Xie C, Li N, Wang H, He C, Hu Y, Peng C. et al. Inhibition of autophagy aggravates DNA damage response and gastric tumorigenesis via Rad51 ubiquitination in response to infection. Gut Microbes. 2020;11:1567-89
58. Ling S, Xie H, Yang F, Shan Q, Dai H, Zhuo J. et al. Metformin potentiates the effect of arsenic trioxide suppressing intrahepatic cholangiocarcinoma: roles of p38 MAPK, ERK3, and mTORC1. J Hematol Oncol. 2017;10:59
59. Platanias LC. Map kinase signaling pathways and hematologic malignancies. Blood. 2003;101:4667-79
60. Shirakihara T, Yamaguchi H, Kondo T, Yashiro M, Sakai R. Transferrin receptor 1 promotes the fibroblast growth factor receptor-mediated oncogenic potential of diffused-type gastric cancer. Oncogene. 2022;41:2587-96
61. Wilson MT, Reeder BJ. Oxygen-binding haem proteins. Exp Physiol. 2008;93:128-32
62. Rouault TA, Tong W-H. Iron-sulphur cluster biogenesis and mitochondrial iron homeostasis. Nat Rev Mol Cell Biol. 2005;6:345-51
63. Costas M, Mehn MP, Jensen MP, Que L. Dioxygen activation at mononuclear nonheme iron active sites: enzymes, models, and intermediates. Chem Rev. 2004;104:939-86
64. Wu Y, Ma Z, Mai X, Liu X, Li P, Qi X. et al. Identification of a Novel Inhibitor of TfR1 from Designed and Synthesized Muriceidine A Derivatives. Antioxidants (Basel, Switzerland). 2022 11
65. Thanee M, Loilome W, Techasen A, Sugihara E, Okazaki S, Abe S. et al. CD44 variant-dependent redox status regulation in liver fluke-associated cholangiocarcinoma: A target for cholangiocarcinoma treatment. Cancer Sci. 2016 107
Corresponding authors: Prof. Xiangcheng Li Address: 300# Guangzhou Road, Nanjing, Jiangsu Province, China, 210000. E-mail: drxcliedu.cn Dr. Changxian Li Address: 300# Guangzhou Road, Nanjing, Jiangsu Province, China, 210000. E-mail: doclicx20com
Received 2022-11-11
Accepted 2023-1-27
Published 2023-2-21