- Primary research
- Open Access
NAT10 as a potential prognostic biomarker and therapeutic target for HNSCC
Cancer Cell International volume 21, Article number: 413 (2021)
Increasing evidence has demonstrated the critical roles of mRNA modification regulators on multiple types of cancers. However, it is still poorly known about the prognostic and therapeutic value of mRNA modification regulators in HNSCC.
The gene expression profile of 36 mRNA modification regulators and their corresponding clinical data were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Stepwise regression in R with both directions was used to construct a model for the prognosis of HNSCC. Univariate Cox regression survival analysis was performed to identify the most significant risk gene. Gene set enrichment analysis (GSEA) was applied to determine the cancer-associated pathways with NAT10. Immunohistochemistry (IHC) staining was performed to evaluate the expression of NAT10 in formalin fixed paraffin-embedded (FFPE) samples of HNSCC. Univariate and multivariate Cox regression survival analysis performed to identify the independent risk factors associated with the OS of patients with HNSCC. HNSCC cell lines (Cal-27, FaDu, and Detroit-562) were transfected with short interfering RNA (siRNA) targeting NAT10 or treated with Remodelin, a small-molecule inhibitor of NAT10. Knockdown efficiency of siRNA was assessed by quantitative real-time PCR (qRT-PCR) and western blotting. In addition, CCK-8 assay, scratch assay and transwell assay were used to examine the proliferation, migration, and invasion abilities of the three HNSCC cell lines after NAT10 was inhibited genetically and pharmaceutically. Cell cycle and cell apoptosis assays were performed by flow cytometry. Finally, the therapeutic value of Remodelin in HNSCC was evaluated via a patient-derived xenograft (PDX) model. The statistical analysis was performed with SPSS 23.0.
A risk prediction model containing 10 mRNA modification regulators was constructed and showed prognostic value in HNSCC. NAT10 was further identified as a key risk gene and independent prognostic factor in TCGA HNSCC dataset. The GSEA analysis suggested that high NAT10 expression was associated with MYC, E2F, G2M checkpoint, mTORC1, DNA repair and oxidative phosphorylation pathways. NAT10 protein expression was significantly up-regulated in tumour cells compared to normal epithelial cells in FFPE samples and increased NAT10 protein expression was correlated with poor overall survival of 267 HNSCC patients. Genetic depletion of NAT10 using siRNA or chemical inhibition of NAT10 using Remodelin resulted in reduced cell proliferation, migration and invasion abilities in Cal-27, FaDu and Detroit-562 cells. Knockdown of NAT10 using siRNA significantly increased cell cycle arrest in S/G2-phase. Remodelin significantly inhibited tumour growth and tumour cell proliferation in the PDX model of HNSCC.
NAT10 could be a potential prognostic marker and a therapeutic target for HNSCC.
Head and neck squamous cell carcinoma (HNSCC), arising from the lip, oral cavity, oropharynx, hypopharynx and larynx, is the eighth most frequent cancer worldwide, accounting for almost 800,000 new cases and approximately 390,000 cancer-related deaths in 2020 . Despite significant advances in HNSCC treatments including surgery, chemotherapy, radiotherapy, targeted therapy and immunotherapy, the 5-year survival rate and long-term prognosis of these patients have not dramatically improved in recent decades [2, 3]. HNSCC is a remarkably heterogeneous malignant tumour characterized by biological and clinical heterogeneity . Currently, the generally used clinical classification system for HNSCC patients is TNM staging, however the patients with the same anatomical site and TNM stage differ in prognosis and clinical outcomes . Tobacco and alcohol consumption are the classical risk factors for HNSCC . Recently, human papilloma virus (HPV) infection has become the most recognized key factor in HNSCC . In addition, numerous molecular and clinical risk factors have been investigated with limited clinical utility . Thus, new prognostic biomarkers, therapeutic targets and treatment strategies for patients with HNSCC are urgently needed.
It is well accepted that gene expression is regulated in a complex manner through genetic and epigenetic control . Genetic abnormalities, including gene mutation, deletion, amplification, or chromosomal translocation, and abnormal epigenetic alterations, such as DNA or histone modification, may lead to the development of human diseases, including cancers [10, 11]. Recent advances in mRNA modification have opened an entirely new layer of gene expression regulation in human cancer [12, 13]. Growing evidence indicates that mRNA modifications and associated regulators contribute to the initiation, progression, and drug response of cancers [12, 14]. Only a few types of mRNA modifications have been revealed, such as N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 3-methylcytosine (m3C), pseudouridine (Ψ) and N4-acetylcytidine (ac4C) [15,16,17,18]. The turnover of each mRNA modification is regulated by related enzymes called regulators, such as METTL3–METTL14 as the methyltransferase for m6A, FTO as the demethylase for m6A and YTHDC2 as the m6A reader . Currently, the roles of mRNA modification regulators in cancer is still poorly understood, especially in HNSCC. Moreover, the prognostic and therapeutic potential of mRNA modification regulators for HNSCC remain unclear.
In this study, we analyzed 36 mRNA modification regulators and constructed a risk prediction model containing 10 potential prognostic risk signatures in The Cancer Genome Atlas (TCGA) HNSCC dataset. Of the candidate genes in these potential prognostic risk signatures, NAT10 was the most promising prognostic risk gene in the TCGA and GEO datasets. The HNSCC dataset from TCGA revealed that the highly expressed NAT10 gene is an independent prognostic factor of HNSCC. IHC analyses revealed high protein expression of NAT10 in HNSCC FFPE samples, which indicated a poor overall survival rate in 267 HNSCC patients. In addition, NAT10 plays an essential role in regulating the proliferation, migration, invasion and cell cycle of tumour cells, thus promoting tumour growth and progression. Remodelin, a small-molecule inhibitor of NAT10, significantly suppressed the growth of HNSCC in a patient-derived xenograft (PDX) model, indicating that targeting NAT10 may be a new therapeutic strategy for HNSCC patients.
Dataset acquisition and screening of mRNA modification regulators
TCGA (https://www.cancer.gov/tcga.) is a public database containing clinical information and RNA-seq results of HNSCC patients, which is widely used for cancer research . Our data were downloaded from TCGA through the FireBrowse (http://firebrowse.org/?cohort=HNSC&download_dialog=true#) website. Data were obtained for 566 samples (520 were primary HNSCC samples, two were metastatic samples and 44 were normal samples) with clinical and RNA-seq information. In addition, Gene Expression Omnibus (GEO) contains many datasets that could be used for data mining . The GSE65858 dataset contains 270 samples with complete clinical and mRNA expression data was downloaded. RNA-seq data of 36 mRNA modification regulators were selected and visualized in a heatmap by R software (version 3.6.3). Tumour stage, age, clinical T, N and M stage, gender, pathological grade and modification type were considered for heatmap clustering.
Construction of a risk prediction model for HNSCC
Stepwise regression in R with both directions in the mRNA modification regulator gene list based on the TCGA HNSCC cohort was used to construct a model for the prognosis of HNSCC, and 10 risk genes associated with the prognosis of HNSCC were screened. Pearson correlation analysis was performed among the 10 genes. The risk score of each gene was calculated by multiplying the expression of the gene and its coefficient as reported previously . In addition, the risk score for each patient was defined as the sum of each gene’s score. Patients were divided into a high-risk group and a low-risk group based on the median value of their risk score. The expression levels of the 10 risk genes in the high-risk group and low-risk group of HNSCC as well as normal samples were visualized. The distributions of the risk scores and overall survival status were visualized, and the Kaplan–Meier method was applied to generate overall survival (OS) curves for the high- and low-risk groups. The log-rank test was used to evaluate statistical significance. To identify genes correlated with OS, univariate Cox regression survival analysis was performed, the hazard ratio and P values were calculated, and a forest plot was used to display the results.
Gene set enrichment analysis (GSEA) of NAT10
GSEA is used to explore the consistency and variability of phenotypes to find common or different biofunctions . GSEA software (version 4.0.3, https://www.gsea-msigdb.org/gsea/downloads.jsp) was applied to identify the pathways associated with NAT10 in HNSCC. The 520 samples were divided into high and low expression groups based on the expression of NAT10, the median expression level of NAT10 was set as the cut-off, and there were 260 samples in each group. The files needed to run GSEA were made according to the help information of the GSEA website, and the predefined “hallmark gene sets” were used for analyses. The normalized enrichment score (NES) was obtained by running GSEA, and the cut-off for significance was defined as follows: a nominal P-value of < 0.05 and FDR q-value of < 0.05.
H&E and IHC staining to detect the expression of NAT10 in 267 HNSCC samples
HNSCC samples of 267 patients (215 patients diagnosed with oral squamous carcinoma and 52 diagnosed with oropharyngeal squamous carcinoma) were obtained from Shanghai 9th People’s Hospital with approval from the ethics committee and informed consent was obtained from the patients. Haematoxylin and eosin (H&E) and the immunohistochemistry (IHC) staining of NAT10 (NAT10 antibody, Ab194297, Abcam, 1:1000 dilution) of samples were performed. For HE and IHC staining, the procedure followed a protocol reported previously . The expression of NAT10 in tumours was estimated by two professional pathologists and defined as follows: high expression: cytoplasm or nucleus was stained dark brown; low expression: cytoplasm or nucleus was stained light brown, and no expression: cytoplasm or nucleus was not stained.
Survival analysis based on 267 HNSCC patients
The gender, age, tumour site, grade, tumour stage, HPV status, tobacco and alcohol use, and primary or recurrent tumour status of 267 patients were collected. Phone calls were used to follow up the survival information of patients. Univariate and multivariate Cox regression survival analysis was performed to identify the independent risk factors associated with the OS of patients. Hazard ratios and p values were calculated by SPSS 23.0. For overall survival analyses, the log-rank test was performed, and a p value of < 0.05 was considered statistically significant.
Construction of short interfering RNA (siRNA) and transfection of HNSCC cell lines
The siRNAs targeting NAT10 were constructed (Hanbio Biotechnology Co., Ltd., Shanghai). The sequences were as follows: si-NAT10-1 (5ʹ-GCACCACUGCUGAGAAUAATT-3ʹ), si-NAT10-2 (5ʹ-GCUCCUCAAGUUCUGGAAATT-3ʹ), si-NAT10-3 (5ʹ-GCAUGGACCUCUCUGAAUATT-3ʹ) and si-Scramble (5ʹ-UAUUCAGAGAGGUCCAUGCTT-3ʹ). 5-Carboxyfluorescein (FAM) was applied to trace siRNA in transfected Cal-27, FaDu and Detroit-562 cells. Cal-27 cells were purchased from the American Type Culture Collection (ATCC). FaDu and Detroit-562 cells were purchased from the National Collection of Authenticated Cell Cultures. Cell lines were confirmed to be free of contamination by short tandem repeat analysis. To verify the efficiency of NAT10 downregulation by siRNA, tumour cells were seeded in 6-well plates. 40 nM siRNAs and 20 µL siRNA Fit solution were added to each well. After culture for 6–17 h, the medium was removed and replaced with fresh complete medium. The cells were collected at 48 h after transfection for reverse transcription polymerase chain reaction and western blotting.
Quantitative real-time PCR (qRT-PCR)
Total RNA of the siRNA- and si-scramble-transfected cells was obtained, ReverTra Ace® qPCR RT Kit (TOYOBO) was used for reverse transcription. SYBR Green Real-time RCR Master Mix (TOYOBO) was applied for quantitative real-time PCR analysis according to the instructions. The primers were designed as follows: NAT10(F: GCAGCCACAAACATTCGCTA; R: AGGAGGATGACCACTAGCCC). The primers for the standard internal gene were designed as follows: GAPDH(F: TCAAGGCTGAGAACGGGAAG; R: TCGCCCCACTTGATTTTGGA); The primers of the top 10 ranked genes in MYC targets were as follows: LDHA (F: TAGCAGATTTGGCAGAGAGTAT; R: CAAGGAACACTAAGGAAGACAT); MYC (F: CGGAAACGACGAGAACAGT; R: CATAGGTGATTGCTCAGGACAT); HSPD1 (F: CACTCGTCTTGAATAGGCTAAA; R: AATCCCTCTTCTCCAAACACT); CDC20 (F: ATTCCTTCCCTGCCAGAC; R: GCCAGTACATTCCCAGAACTC); RRP9 (F: CGTGAAGGTGTGGAATGT; R: GTCCTGGTGTCCGAAGA); ODC1 (F: ATGTGAATGATGGCGTCTAT; R: ACTTCTCATCTGGTTTAGGTCTC); NOP16 (F: CTCGGGACCTCATTGACTAT; R: GTCACTCCACCTCCATCTTC); NOLC1 (F: AAAGGCGGCAGTGGTAGT; R: GGAAATGTGTTAGGGGTCTGA); KPNA2 (F: GCCCGTCTTCACAGATTC; R: GCTCCACATTGACCTCTATTC); PHB (F: GTAGGGGAAGGGACTCATT; R: GCAGTGTGATGTTGACATTCT); b-actin (F: AAGGTGACAGCAGTCGGTT; R: TGTGTGGACTTGGGAGAGG). The primers were provided by BioSune Biotechnology Co., Ltd., Shanghai. The reaction procedures were as follows: 95 °C, 180 s; 94 °C, 15 s; 60 °C, 30 s; 72 °C, 30 s. The melting curve procedures were 95 °C, 10 s; 65 °C, 60 s; 97 °C, 1 s.
Western blotting was applied to verify the expression of NAT10 in tumour cells at the protein level after transfection with siRNA. Total protein was extracted and the concentration was measured. Twenty micrograms of the protein of each sample were loaded for sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and subsequently transferred to polyvinylidene fluoride (PVDF) membranes. Then, the membranes were blocked with 5% skim milk and TBST was used for membranes washing. Primary antibody (NAT10 antibody, Ab194297, Abcam, 1:1000 dilution; GAPDH antibody: 5174s, CST) was incubated in 4 °C overnight, followed by incubation with the horseradish peroxidase-conjugated (HRP) secondary antibodies for 1 h at room temperature. Finally, the detection of the immunoreactive proteins were visualized by enhanced chemiluminescence (ECL, Millipore) according to the manufacturer's instructions.
In vitro studies
Three cell lines mentioned above were used for in vitro studies. For the cell proliferation assay, cells were pretreated with NAT10 siRNA or NAT10 inhibitor Remodelin hydrobromide  (Product ID: S7641, Selleck) at a concentration range of 5–60 μM. A total of 3000 cells were seeded in 96-well plates, and Cell Counting Kit-8 (CCK-8, Dojindo, Japan) was diluted in DMEM (1:10) to measure the optical density (OD) value at 450 nm on a multimode microplate reader (SpectraMax i3, Molecular Devices).
To determine the migration and invasion ability of the tumour cells, the three cell lines were transfected with siRNA or treated with Remodelin for 48 h. Cells were seeded into 6-well plates, and a scratch assay was performed with 10-μL pipette tips. Serum-free Dulbecco’s modified Eagle’s medium (DMEM) was added. After culturing for 24 h, pictures were taken, and ImageJ was used to calculate the area without cells. In addition, 105 cells were seeded in the apical chamber of the transwell in 24-well cell culture plates with serum-free DMEM, and DMEM containing 20% FBS was added to the basolateral chamber of the transwell. After culturing for 48 h, 4% paraformaldehyde (PFA) was used for fixation, and crystal violet staining was performed. For the invasion assay, the apical chamber of the transwell was precoated with Matrigel. The other procedures were the same as mentioned above.
Flow cytometry was used to explore the changes in the cell cycle and cell apoptosis after transfection with NAT10 siRNA for 48 h, and scrambled siRNA was used as a control. For the cell cycle assay, cells were fixed with 70% alcohol overnight and then stained with mixed stain buffer (50 μg/mL propidium iodide (PI), 100 μg/mL RNase A, 0.2% Triton X-100) at 4 °C for 30 min. For the apoptosis assay, cells were collected and stained according to the manufacturer’s instructions for the FITC Annexin V apoptosis detection kit (Product ID: 556547, BD Pharmingen). BD FACSCalibur was used to detect positive events, and FlowJo (version 10.6.2) was used for data analysis.
Evaluation of the therapeutic value of Remodelin via the PDX model
BALB/c nude mice (6 weeks, male, approximately 20 g) were obtained from Shanghai Sippr-BK Laboratory Animal Co. Ltd. and maintained in specific-pathogen-free facilities at Shanghai Ninth People’s Hospital. All experimental procedures were approved by the Laboratory Animal Care and Use Committees of the hospital. PDX model was constructed as we reported previously . Each mouse was anesthetized with inhalation of 2–4% isoflurane/O2 for less than 30 min, then tumour tissues from patients were cut into 20–30 mm3 pieces and seeded in the flanks of mice. When the tumour volume reached 1500–2000 mm3, the mice were sacrificed. Tumour tissues were cut into pieces and seeded again for passaging. The expression of NAT10 (NAT10 antibody, Ab194297, Abcam, 1:1000 dilution) and Ki-67 (Ki-67 antibody, Ab927742, Abcam, 1:1000 dilution) in the PDX model were evaluated by IHC. Third-passage mice with tumour volumes of 150–250 mm3 were used to evaluate the therapeutic value of the NAT10 inhibitor (Remodelin). Remodelin was dissolved in dimethyl sulfoxide (DMSO) and diluted in sterile Tween 80 and 45% 2-hydroxypropyl-b-cyclodextrin solution . The final formulation was 20% DMSO, 65% (45% 2-hydroxypropyl-b-cyclodextrin solution) and 15% Tween 80. The inhibitor and corresponding vehicle were administered intragastrically daily at 100 mg/kg/day. Tumour volume (TV) was recorded every 3–5 days, and mice were sacrificed after therapy for 4 weeks. The tumour growth inhibition rate (TGI) was calculated as
SPSS 23.0 was used for statistical analyses, and GraphPad Prism 9.0.2 was applied to generate figures. For survival analyses, the log-rank test was performed, and a P-value of < 0.05 was considered to indicate a significant difference.
In vitro and in vivo results are presented as means ± SD (standard deviation). A t test was used to compare two groups, and one-way ANOVA was used to compare three or more groups. Two-way ANOVA was used to evaluate the therapeutic value of Remodelin, and P-values < 0.05 were considered to indicate statistical significance.
The constructed risk prediction model showed prognostic value in HNSCC
To construct a model to predict the prognosis of patients, 36 modification regulators of mRNA were selected and analyzed (Additional file 1: Table S1). The workflow chart of the study is shown in Fig. 1a. The mRNA levels of these regulators in the TCGA cohort and their corresponding clinical information were visualized (Fig. 1b). Stepwise regression in R with both direction methods was applied, and a risk prediction model containing 10 modification regulators was constructed. The risk score of each patient was calculated. Patients were defined as the high-risk and low-risk groups based on their risk scores, and shorter survival time was observed with increasing risk score (Fig. 1c). In addition, a heatmap was generated to visualize the expression of the 10 risk genes (Fig. 1c). The mRNA levels of the 10 genes among the high-risk, low-risk and normal groups were also displayed (Fig. 1d). The OS analysis illustrated that the high-risk group had a lower survival probability than the low-risk group (Fig. 1e).
NAT10 was a key risk gene and independent prognostic factor of HNSCC
To further identify the most essential risk gene of the model, univariate Cox analysis was performed on TCGA and GEO datasets, and the hazard ratio of each gene was calculated (RBM15B was not found in the GEO dataset). The results indicated that NAT10 was the most significant gene with a higher hazard ratio among the 10 modification regulators (Fig. 2a, b). In addition, the expression of NAT10 was upregulated in HNSCC samples than in normal samples (Fig. 2c), and a lower OS probability was observed in the NAT10 high expression group of HNSCC patients (Fig. 2d). GSEA was performed to explore the enriched pathways associated with the high expression of NAT10. As shown in Fig. 2e, the gene sets were enriched in the MYC, E2F, G2M checkpoint, mTORC1, DNA repair and oxidative phosphorylation pathways. These pathways are associated with tumour progression, which indicates the essential role of NAT10 in HNSCC.
High protein expression of NAT10 indicated worse OS in HNSCC patients
To validate the prognostic role of NAT10 at the protein level, H&E and IHC staining of NAT10 from 267 HNSCC clinical FFPE samples were performed, and representative images of negative, low and high staining of NAT10 are shown (Fig. 3a). NAT10 was mainly located in the nucleus and more highly expressed in the tumour cells than that in normal epithelial cells (Fig. 3b). Based on the expression of NAT10, patients were divided into three groups (NAT10 negative, NAT10 low and NAT10 high group), and the clinical pathological parameters in these groups of patients were archived and compared with each other. We found that age (≥ 65), tumour site (oral cavity) and tumour type (recurrent) were strongly associated with NAT10 high group (Table 1). Univariate Cox regression survival analysis was performed to identify the risk factors for overall survival, and the results illustrated that NAT10 high group, age (≥ 65), tumour stage (III and IV), HPV-negative status, and recurrent tumour type were risk factors for OS (Additional file 2: Table S2). Low survival probability was found in patients with high expression of NAT10 (Fig. 3c), age ≥ 65 (Fig. 3d), stage III and IV disease (Fig. 3e), HPV-negative status (Fig. 3f) and recurrent tumours (Fig. 3g). Multivariate Cox regression survival analysis was further performed among the five mentioned risk factors above, and the results revealed that high expression of NAT10, age ≥ 65 and recurrent tumours were independent risk factors for HNSCC patients’ OS (Table 2).
Downregulation of NAT10 inhibited the proliferation ability of three HNSCC cell line cells
Three siRNAs of NAT10 and si-Scramble were constructed, and their knockdown efficacy in three HNSCC cell lines (Cal-27, FaDu and Detroit-562) was verified by qRT-PCR and western blotting (Additional file 3: Fig. S1). To reveal the role of NAT10 in cell proliferation ability, the three cell lines were treated with two siRNAs of NAT10, and si-Scramble was administered as a control. Reduced proliferation ability was observed after treating the cell lines with siRNAs (Fig. 4a). In addition, Remodelin as a chemical inhibitor of NAT10 was applied at different concentrations and the CCK-8 assay indicated that Remodelin inhibited the proliferation of HNSCC cells. The half-maximal inhibitory concentrations (IC50s) of Remodelin ranged from 14.38 to 32.77 µM (Fig. 4b). Remodelin significantly inhibited the proliferation ability, and proliferation rates were decreased as the concentration were increased (Fig. 4c). Cell cycles were arrested in the S and G2 stages by NAT10 siRNA (Fig. 4d). However, there was no difference in cell apoptosis between the si-Scramble and si-NAT10 groups (Fig. 4e).
NAT10 inhibition decreased the migration and invasion abilities of HNSCC cell lines
Tumour cell migration and invasion are essential for cancer progression. A cell scratch assay was used to evaluate the horizontal migration ability. Our results demonstrated that the migration ability was decreased by treatment of Remodelin and NAT10 siRNA (Fig. 5a–f). Transwell chambers were used in the cell migration and invasion assays. For the cell invasion assay, the apical chamber was precoated with Matrigel. Cell lines were treated with Remodelin or NAT10 siRNA. Approximately 105 cells were seeded in the apical chamber of the transwell. Crystal violet staining revealed that the migration ability and invasion ability of HNSCC cells were decreased after NAT10 was inhibited (Fig. 5g).
Remodelin inhibited the growth of HNSCC in PDX model
To evaluate the therapeutic value of Remodelin in HNSCC in vivo, PDX models were constructed. Remodelin was given daily by intragastric administration for 4 weeks. Tumour volume was recorded and TGI was calculated. As a result, Remodelin significantly inhibited the growth of HNSCC PDXs in vivo. The TGI was 55.7%, and the mice’s body weights were not affected (Fig. 6a, b). In addition, the tumour weight was significantly decreased in the treatment group compared with the vehicle group (Fig. 6c, d). IHC staining were performed to detect the expression of NAT10 and Ki-67 in tumours from patients (Fig. 6e), the PDX model without Remodelin administration (Fig. 6f), and the PDX model with Remodelin administration (Fig. 6g). As shown in the pictures, lower expression of Ki-67 was observed in the Remodelin treatment group, which indicated that Remodelin inhibited tumour proliferation ability in vivo.
Remodelin regulated the expression of MYC and LDHA
MYC targets is the most enriched pathway in GSEA results, qRT-PCR was performed to evaluate the expression of the top 10 enriched genes of MYC targets in three cell lines after treated with Remodelin for 48 h (Fig. 7a). Among the top 10 genes, MYC was significantly downregulated by Remodelin, but the expression of LDHA was upregulated in all three cell lines (Fig. 7b, c). The results indicated MYC and LDHA regulation was associated with the inhibition of NAT10.
It is clear that the mRNA modifications and related regulators play critical roles in tumorigenesis and progression of various cancers . In this study, we firstly applied bioinformatic methods to construct a risk prediction model containing 10 mRNA modification regulators related to six reported distinct forms of mRNA modifications. Moreover, survival analysis of this prediction model showed that prognosis of the high-risk group had a lower OS survival probability than the low-risk group. Based on univariate Cox analysis performed on TCGA and GEO datasets, NAT10 was identified as an independent essential gene in this risk model. In addition, we further identified NAT10 as an independent prognostic factor associated with the OS of HNSCC patients in our clinical cohort. Furthermore, the cell proliferation, migration and invasion abilities of three HNSCC cell lines were inhibited by genetic depletion of NAT10 using siRNA or chemical inhibitor, Remodelin. Knockdown of NAT10 using siRNA increased cell cycle arrest in S/G2-phase. Finally, the therapeutic value of Remodelin was evaluated in a HNSCC PDX model and found to inhibit the growth of tumours.
It is of great significance to identify risk models and biomarkers for predicting the prognosis of patients and the selection of appropriate clinical strategies . In our study, patients with HNSCC could be stratified into a high-risk cohort with worse OS and a low-risk cohort with better OS based on our risk prediction model in TCGA dataset. Univariate Cox analysis further revealed that NAT10 was the most significant gene with the highest hazard ratio among these model signatures in TCGA and GEO datasets. NAT10 (N-acetyltransferase 10) is a histone acetyltransferase associated with many biological processes mainly distributed in the nucleolus [28,29,30]. With acetyltransferase and RNA binding activities, NAT10 is the only human enzyme to catalyze ac4C production in eukaryotic RNA . In-depth studies demonstrated that ac4C as an mRNA modification plays a key role in controlling telomerase activity, the synthesis of rRNA, the repair of DNA damage and the stabilization of mRNA. It is also related to the development, progression, and prognosis of various types of human diseases, especially cancer. [16, 31, 32]. As the only ac4C modification producer, NAT10 is significantly associated with the prognosis of several malignant tumours and potently promotes the tumour progression [30, 33,34,35,36,37,38,39]. High expression of NAT10 in hepatocellular carcinoma [30, 37, 39] and colorectal cancer  is associated with a worse prognosis in patients. However, few studies have assessed the role of NAT10 in HNSCC. Based on public data from the TCGA and 267 HNSCC samples in our cohort, we identified NAT10 as a potential prognostic biomarker in HNSCC patients. In this study, NAT10 was more highly expressed in HNSCC tissues than in normal tissues, and that upregulated expression of NAT10 indicated worse OS of HNSCC. In hepatocellular carcinoma, the tumorigenic mechanism involves the redistribution of NAT10 from nucleoli to karyotheca, the cytoplasm and the cell membrane, and NAT10 can upregulate mutant p53 (a tumour suppressor gene) levels [30, 37, 39]. In addition, laryngeal cancer cell migration was reported to be promoted by the MYCT1/NAT10 axis . Likewise, MYC targets had the highest normalized enrichment score (NES) among all the gene sets in our GSEA results, and qRT-PCR results of the top 10 ranked genes of MYC targets indicated Remodelin downregulated the expression of MYC and upregulated the expression of LDHA in HNSCC cell lines.
Increasing evidence indicated that NAT10 might be a potential therapeutic target for cancer treatment [40,41,42,43]. Therefore, we analyzed the therapeutic potential of NAT10 in HNSCC and observed that the proliferation, migration and invasion abilities of HNSCC cell lines were decreased after NAT10 was inhibited by siRNA and Remodelin in this study. Besides, our result showed that the cell cycle of the three HNSCC cell lines was arrested in the G2/S stage after downregulation of NAT10 by siRNA. Interestingly, our GSEA results also revealed gene sets in the NAT10 high expression group involved in the G2/M checkpoint and associated with DNA repair. NAT10-mediated acetylation of MORC2 could force the cell to pass through the G2 checkpoint and confer resistance to DNA-damaging treatments in a breast cancer study . Furthermore, our in vivo study yielded preclinical data relevant to the treatment of HNSCC in the clinic. Specifically, Remodelin inhibited the growth of a HNSCC-PDX model with high expression of NAT10. Remodelin is a small-molecule compound explored as potent inhibitor of NAT10 and firstly reported that targeting NAT10 might be an alternative strategy for the treatment of laminopathies and aging . More recently, Remodelin showed a reverse effect on cell proliferation, cell invasion, and migration of several tumours [37, 43]. Although the research of Remodelin is still at a very early stage, our results indicate that NAT10 may be a promising drug target. Moreover, inhibition of NAT10 might be a new therapeutic strategy for HNSCC. We anticipate that Remodelin could be examined as a drug for HNSCC in clinical trials in the very near future.
Tumour cell proliferation, migration, invasion, and G2/S stage promotion were proven to be regulated by NAT10, thus aggravating the progression of HNSCC in the present study. However, one limitation of our study is that the other GSEA results need to be verified in further studies. For example, the specific role of NAT10 in E2F target, and mTORC1 signalling in HNSCC needs in-depth research. In addition, oxidative phosphorylation is essential for tumour growth  and it is necessary to explore whether NAT10 contributes to the metabolism of tumours. Moreover, targeting NAT10 alone or in combination with other drugs may be a new strategy for tumour treatment in further investigations.
Above all, this study reported the construction of a risk prediction model and the identification of NAT10 as a critical tumour growth-favoring factor in HNSCC. High expression of NAT10 was associated with an unfavorable prognosis in HNSCC patients in the TCGA cohort and was verified in our clinical cohort, which included 267 patients. NAT10 promoted tumour growth by regulating the migration, invasion and cell cycle of tumour cells. Remodelin, an inhibitor of NAT10, significantly suppressed the growth of HNSCC in a PDX model, indicating that Remodelin may be a new candidate drug for HNSCC treatment.
In conclusion, NAT10 was a key gene in HNSCC development and progression. As verified by IHC results of 267 HNSCC samples, NAT10 is associated with the poor prognosis in HNSCC patients. In addition, in vitro and in vivo studies proved that NAT10 could potentially be the therapeutic target in HNSCC. The migration, invasion and cell cycle promotion of the tumour cells may be the mechanisms by which NAT10 promotes tumour progression.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in GEO (http://www.ncbi.nlm.nih.gov/geo) and TCGA through firebrowse (http://firebrowse.org/?cohort=HNSC&download_dialog=true#) website.
American Type Culture Collection
Cell Counting Kit-8
Dulbecco’s Modified Eagle Medium
Formalin fixed paraffin-embedded
Gene Expression Omnibus
Gene set enrichment analysis
Hematoxylin and eosin
Head and neck squamous cell carcinoma
Human papilloma virus
Half maximal inhibitory concentration
Microrchidia family CW-type zinc finger 2
MYC target 1
Normalized enrichment score
Quantitative real-time polymerase chain reaction
The Cancer Genome Atlas
Tumour growth inhibition rate
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Chow LQM. Head and neck cancer. N Engl J Med. 2020;382(1):60–72.
Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. 2020;6(1):92.
Leemans CR, Snijders PJF, Brakenhoff RH. The molecular landscape of head and neck cancer. Nat Rev Cancer. 2018;18(5):269–82.
Hammerman PS, Hayes DN, Grandis JR. Therapeutic insights from genomic studies of head and neck squamous cell carcinomas. Cancer Discov. 2015;5(3):239–44.
Cramer JD, Burtness B, Le QT, Ferris RL. The changing therapeutic landscape of head and neck cancer. Nat Rev Clin Oncol. 2019;16(11):669–83.
Lydiatt WM, Patel SG, O’Sullivan B, Brandwein MS, Ridge JA, Migliacci JC, Loomis AM, Shah JP. Head and neck cancers-major changes in the American Joint Committee on cancer eighth edition cancer staging manual. CA Cancer J Clin. 2017;67(2):122–37.
Sepiashvili L, Bruce JP, Huang SH, O’Sullivan B, Liu FF, Kislinger T. Novel insights into head and neck cancer using next-generation “omic” technologies. Cancer Res. 2015;75(3):480–6.
Deng X, Su R, Weng H, Huang H, Li Z, Chen J. RNA N(6)-methyladenosine modification in cancers: current status and perspectives. Cell Res. 2018;28(5):507–17.
Lu Y, Chan YT, Tan HY, Li S, Wang N, Feng Y. Epigenetic regulation in human cancer: the potential role of epi-drug in cancer therapy. Mol Cancer. 2020;19(1):79.
Jones PA, Issa JP, Baylin S. Targeting the cancer epigenome for therapy. Nat Rev Genet. 2016;17(10):630–41.
Boriack-Sjodin PA, Ribich S, Copeland RA. RNA-modifying proteins as anticancer drug targets. Nat Rev Drug Discov. 2018;17(6):435–53.
Barbieri I, Kouzarides T. Role of RNA modifications in cancer. Nat Rev Cancer. 2020;20(6):303–22.
Chen XY, Zhang J, Zhu JS. The role of m(6)A RNA methylation in human cancer. Mol Cancer. 2019;18(1):103.
Zaccara S, Ries RJ, Jaffrey SR. Reading, writing and erasing mRNA methylation. Nat Rev Mol Cell Biol. 2019;20(10):608–24.
Jin G, Xu M, Zou M, Duan S. The processing, gene regulation, biological functions, and clinical relevance of N4-acetylcytidine on RNA: a systematic review. Mol Ther Nucleic Acids. 2020;20:13–24.
Davalos V, Blanco S, Esteller M. SnapShot: messenger RNA modifications. Cell. 2018;174(2):498-498.e491.
Wiener D, Schwartz S. The epitranscriptome beyond m(6)A. Nat Rev Genet. 2021;22(2):119–31.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42.
Hutter C, Zenklusen JC. The cancer genome atlas: creating lasting value beyond its data. Cell. 2018;173(2):283–5.
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Zhao X, Cui L. Development and validation of a mA RNA methylation regulators-based signature for predicting the prognosis of head and neck squamous cell carcinoma. Am J Cancer Res. 2019;9(10):2156–69.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Zhou R, Shi C, Tao W, Li J, Wu J, Han Y, Yang G, Gu Z, Xu S, Wang Y, et al. Analysis of mucosal melanoma whole-genome landscapes reveals clinically relevant genomic aberrations. Clin Cancer Res. 2019;25(12):3548–60.
Larrieu D, Britton S, Demir M, Rodriguez R, Jackson SP. Chemical inhibition of NAT10 corrects defects of laminopathic cells. Science. 2014;344(6183):527–32.
Balmus G, Larrieu D, Barros AC, Collins C, Abrudan M, Demir M, Geisler NJ, Lelliott CJ, White JK, Karp NA, et al. Targeting of NAT10 enhances healthspan in a mouse model of human accelerated aging syndrome. Nat Commun. 2018;9(1):1700.
Joyner MJ, Paneth N. Promises, promises, and precision medicine. J Clin Investig. 2019;129(3):946–8.
Shen Q, Zheng X, McNutt MA, Guang L, Sun Y, Wang J, Gong Y, Hou L, Zhang B. NAT10, a nucleolar protein, localizes to the midbody and regulates cytokinesis and acetylation of microtubules. Exp Cell Res. 2009;315(10):1653–67.
Zhang H, Hou W, Wang HL, Liu HJ, Jia XY, Zheng XZ, Zou YX, Li X, Hou L, McNutt MA, et al. GSK-3β-regulated N-acetyltransferase 10 is involved in colorectal cancer invasion. Clin Cancer Res. 2014;20(17):4717–29.
Zhang X, Liu J, Yan S, Huang K, Bai Y, Zheng S. High expression of N-acetyltransferase 10: a novel independent prognostic marker of worse outcome in patients with hepatocellular carcinoma. Int J Clin Exp Pathol. 2015;8(11):14765–71.
Arango D, Sturgill D, Alhusaini N, Dillman AA, Sweet TJ, Hanson G, Hosogane M, Sinclair WR, Nanan KK, Mandler MD, et al. Acetylation of cytidine in mRNA promotes translation efficiency. Cell. 2018;175(7):1872-1886.e1824.
Karthiya R, Wasil SM, Khandelia P. Emerging role of N4-acetylcytidine modification of RNA in gene regulation and cellular functions. Mol Biol Rep. 2020;47(11):9189–99.
Liang P, Hu R, Liu Z, Miao M, Jiang H, Li C. NAT10 upregulation indicates a poor prognosis in acute myeloid leukemia. Curr Probl Cancer. 2020;44(2):100491.
Zi J, Han Q, Gu S, McGrath M, Kane S, Song C, Ge Z. Targeting NAT10 induces apoptosis associated with enhancing endoplasmic reticulum stress in acute myeloid leukemia cells. Front Oncol. 2020;10:598107.
Zhang ZX, Zhang WN, Sun YY, Li YH, Xu ZM, Fu WN. CREB promotes laryngeal cancer cell migration via MYCT1/NAT10 axis. Onco Targets Ther. 2018;11:1323–31.
Liu HY, Liu YY, Yang F, Zhang L, Zhang FL, Hu X, Shao ZM, Li DQ. Acetylation of MORC2 by NAT10 regulates cell-cycle checkpoint control and resistance to DNA-damaging chemotherapy and radiotherapy in breast cancer. Nucleic Acids Res. 2020;48(7):3638–56.
Ma R, Chen J, Jiang S, Lin S, Zhang X, Liang X. Up regulation of NAT10 promotes metastasis of hepatocellular carcinoma cells through epithelial-to-mesenchymal transition. Am J Transl Res. 2016;8(10):4215–23.
Zhang H, Hou W, Wang HL, Liu HJ, Jia XY, Zheng XZ, Zou YX, Li X, Hou L, McNutt MA, et al. GSK-3beta-regulated N-acetyltransferase 10 is involved in colorectal cancer invasion. Clin Cancer Res. 2014;20(17):4717–29.
Li Q, Liu X, Jin K, Lu M, Zhang C, Du X, Xing B. NAT10 is upregulated in hepatocellular carcinoma and enhances mutant p53 activity. BMC Cancer. 2017;17(1):605.
Dalhat MH, Altayb HN, Khan MI, Choudhry H. Structural insights of human N-acetyltransferase 10 and identification of its potential novel inhibitors. Sci Rep. 2021;11(1):6051.
Zhang X, Chen J, Jiang S, He S, Bai Y, Zhu L, Ma R, Liang X. N-acetyltransferase 10 enhances doxorubicin resistance in human hepatocellular carcinoma cell lines by promoting the epithelial-to-mesenchymal transition. Oxid Med Cell Longev. 2019;2019:7561879.
Wu J, Zhu H, Wu J, Chen W, Guan X. Inhibition of N-acetyltransferase 10 using remodelin attenuates doxorubicin resistance by reversing the epithelial–mesenchymal transition in breast cancer. Am J Transl Res. 2018;10(1):256–64.
Oh TI, Lee YM, Lim BO, Lim JH. Inhibition of NAT10 suppresses melanogenesis and melanoma growth by attenuating microphthalmia-associated transcription factor (MITF) expression. Int J Mol Sci. 2017;18(9):1924.
Frederick M, Skinner HD, Kazi SA, Sikora AG, Sandulache VC. High expression of oxidative phosphorylation genes predicts improved survival in squamous cell carcinomas of the head and neck and lung. Sci Rep. 2020;10(1):6380.
The authors acknowledge The Ninth People’s Hospital affiliated to the Shanghai Jiao Tong University School of Medicine for financial support.
This work was supported by CAMS Innovation Fund for Medical Sciences (CIFMS, Project No. 2019-I2M-5-037), National Natural Science Foundation of China (Project No. 81872187), Emerging Frontier Technology Joint Research Project (Project No. SHDC12018104), Shanghai Clinical Research Center for Oral Diseases (Project No. 19MC1910600), Shanghai Municipal Key Clinical Specialty (Project No. shslczdzk01601).
Ethics approval and consent to participate
All animal experiments were approved by and performed in accordance with the guidelines of the Ninth People’s Hospital Affiliated to Shanghai Jiao Tong University School of Medicine (SH9H-2021-A012-SB).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The 36 modification regulators of mRNA.
Univariate survival analyses of factors contributing to OS in Chinese HNSCC patients.
The expression of NAT10 was inhibited at the RNA and protein levels by siRNA. (a) qRT-PCR was used to verify the expression of NAT10. After transfection with three siRNAs, the expression of NAT10 was downregulated in the three cell lines. (b) Western blotting was applied to verify the expression of NAT10 at the protein level. NAT10 was strongly expressed in cells transfected with si-Scramble, while NAT10 was expressed at low levels in cells transfected with si- NAT10-1 and si- NAT10-2.
About this article
Cite this article
Tao, W., Tian, G., Xu, S. et al. NAT10 as a potential prognostic biomarker and therapeutic target for HNSCC. Cancer Cell Int 21, 413 (2021). https://doi.org/10.1186/s12935-021-02124-2
- mRNA modification regulator
- Therapeutic target