Association of Topoisomerase II (TOP2A) and Dual-Specificity Phosphatase 6 (DUSP6) Single Nucleotide Polymorphisms with Radiation Treatment Response and Prognosis of Lung Cancer in Han Chinese

Evaluation of the prognosis in patients with small-cell lung cancer treated by chemotherapy using tumor shrinkage rate-based radiomics

BACKGROUND

Small-cell lung cancer (SCLC) is a leading cause of cancer-related death. However, the prognostic value of the tumor shrinkage rate (TSR) after chemotherapy for SCLC is still unknown.

METHODS

We performed a retrospective analysis of 235 patients with SCLC. The TSR cutoff was determined based on receiver-operating characteristic curve analysis. The associations of TSR with progression-free survival (PFS) and overall survival (OS) were assessed using univariate and multivariate Cox proportional hazards models. Survival curves were obtained by the Kaplan-Meier method and compared using the log-rank test. Recurrence patterns after first-line treatment were summarized in a pie chart. A nomogram was constructed to validate the predictive role of the TSR in SCLC.

RESULTS

The TSR cutoff was identified to be  - 6.6%. Median PFS and OS were longer in the group with a TSR < -6.6% than in the group with a TSR ≥ - 6.6%. PFS and OS were also longer in patients with extensive SCLC when the TSR was < - 6.6% than when it was > - 6.6%. Brain metastasis-free survival was better in the group with a TSR < - 6.6%. There was a significant positive correlation between TSR and PFS. Furthermore, univariate and multivariate regression analyses showed that the TSR, patient age, and previous radiotherapy were independent prognostic factors for OS while TSR and M stage were independent prognostic factors for PFS.

CONCLUSIONS

The TSR may prove to be a good indicator of OS and PFS in patients receiving chemotherapy-based first-line treatment for SCLC.

Mortality burden of pre-treatment weight loss in patients with non-small-cell lung cancer: A systematic literature review and meta-analysis

Cachexia, with weight loss (WL) as a major component, is highly prevalent in patients with cancer and indicates a poor prognosis. The primary objective of this study was to conduct a meta-analysis to estimate the risk of mortality associated with cachexia (using established WL criteria prior to treatment initiation) in patients with non-small-cell lung cancer (NSCLC) in studies identified through a systematic literature review. The review was conducted according to PRISMA guidelines. Embase® and PubMed were searched to identify articles on survival outcomes in adult patients with NSCLC (any stage) and cachexia published in English between 1 January 2016 and 10 October 2021. Two independent reviewers screened titles, abstracts and full texts of identified records against predefined inclusion/exclusion criteria. Following a feasibility assessment, a meta-analysis evaluating the impact of cachexia, defined per the international consensus criteria (ICC), or of pre-treatment WL ≥ 5% without a specified time interval, on overall survival in patients with NSCLC was conducted using a random-effects model that included the identified studies as the base case. The impact of heterogeneity was evaluated through sensitivity and subgroup analyses. The standard measures of statistical heterogeneity were calculated. Of the 40 NSCLC publications identified in the review, 20 studies that used the ICC for cachexia or reported WL ≥ 5% and that performed multivariate analyses with hazard ratios (HRs) or Kaplan-Meier curves were included in the feasibility assessment. Of these, 16 studies (80%; n = 6225 patients; published 2016-2021) met the criteria for inclusion in the meta-analysis: 11 studies (69%) used the ICC and 5 studies (31%) used WL ≥ 5%. Combined criteria (ICC plus WL ≥ 5%) were associated with an 82% higher mortality risk versus no cachexia or WL < 5% (pooled HR [95% confidence interval, CI]: 1.82 [1.47, 2.25]). Although statistical heterogeneity was high (I2 = 88%), individual study HRs were directionally aligned with the pooled estimate, and there was considerable overlap in CIs across included studies. A subgroup analysis of studies using the ICC (HR [95% CI]: 2.26 [1.80, 2.83]) or WL ≥ 5% (HR [95% CI]: 1.28 [1.12, 1.46]) showed consistent findings. Assessments of methodological, clinical and statistical heterogeneity indicated that the meta-analysis was robust. Overall, this analysis found that ICC-defined cachexia or WL ≥ 5% was associated with inferior survival in patients with NSCLC. Routine assessment of both weight and weight changes in the oncology clinic may help identify patients with NSCLC at risk for worse survival, better inform clinical decision-making and assess eligibility for cachexia clinical trials.

Aberrant Expression and Prognostic Potential of IL-37 in Human Lung Adenocarcinoma

Interleukin-37 (IL-37) is a relatively new IL-1 family cytokine that, due to its immunoregulatory properties, has lately gained increasing attention in basic and translational biomedical research. Emerging evidence supports the implication of this protein in any human disorder in which immune homeostasis is compromised, including cancer. The aim of this study was to explore the prognostic and/or diagnostic potential of IL-37 and its receptor SIGIRR (single immunoglobulin IL-1-related receptor) in human tumors. We utilized a series of bioinformatics tools and -omics datasets to unravel possible associations of IL-37 and SIGIRR expression levels and genetic aberrations with tumor development, histopathological parameters, distribution of tumor-infiltrating immune cells, and survival rates of patients. Our data revealed that amongst the 17 human malignancies investigated, IL-37 exhibits higher expression levels in tumors of lung adenocarcinoma (LUAD). Moreover, the expression profiles of IL-37 and SIGIRR are associated with LUAD development and tumor stage, whereas their high mRNA levels are favorable prognostic factors for the overall survival of patients. What is more, IL-37 correlates positively with a LUAD-associated transcriptomic signature, and its nucleotide changes and expression levels are linked with distinct infiltration patterns of certain cell subsets known to control LUAD anti-tumor immune responses. Our data indicate the potential value of IL-37 and its receptor SIGIRR to serve as biomarkers and/or immune-checkpoint therapeutic targets for LUAD patients. Further, the data highlight the urgent need for further exploration of this cytokine and the underlying pathogenetic mechanisms to fully elucidate its implication in LUAD development and progression.

TOP2A correlates with poor prognosis and affects radioresistance of medulloblastoma

Radiotherapy remains the standard treatment for medulloblastoma (MB), and the radioresistance contributes to tumor recurrence and poor clinical outcomes. Nuclear DNA topoisomerase II-alpha (TOP2A) is a key catalytic enzyme that initiates DNA replication, and studies have shown that TOP2A is closely related to the therapeutic effects of radiation. In this study, we found that TOP2A was significantly upregulated in MB, and high expression of TOP2A related to poor prognosis of MB patients. Knockdown of TOP2A inhibited MB cell proliferation, migration, and invasion, whereas overexpression of TOP2A enhanced the proliferative and invasive ability of MB cells. Moreover, si-TOP2A transfection in combination with irradiation (IR) significantly reduced the tumorigenicity of MB cells, compared with those transfected with si-TOP2A alone. Cell survival curve analysis revealed that the survival fraction of MB cells was significantly reduced upon TOP2A downregulation and that si-TOP2A-transfected cells had decreased D₀, Dq, and SF₂ values, indicating that TOP2A knockdown suppresses the resistance to radiotherapy in MB cells. In addition, western blot analysis demonstrated that the activity of Wnt/β-catenin signaling pathway was inhibited after TOP2A downregulation alone or in combination with IR treatment, whereas overexpression of TOP2A exhibited the opposite effects. Gene set enrichment analysis also revealed that Wnt/β-catenin signaling pathway is enriched in TOP2A high-expression phenotypes. Collectively, these data indicate that high expression of TOP2A leads to poor prognosis of MB, and downregulation of TOP2A inhibits the malignant behaviour as well as the radioresistance of MB cells. The Wnt/β-catenin signaling pathway may be involved in the molecular mechanisms of TOP2A mediated reduced tumorigenicity and radioresistance of MB cells.

Comprehensive Bioinformatics Analysis of mRNA Expression Profiles and Identification of a miRNA-mRNA Network Associated with the Pathogenesis of Low-Grade Gliomas

Purpose

Low-grade glioma is the most common type of primary intracranial tumour, and the overall survival of patients with low-grade glioma (LGG) has shown no significant improvement over the past few decades. Therefore, it is crucial to understand the precise molecular mechanisms involved in the carcinogenesis of LGG.

Methods

To investigate the regulatory mechanisms of mRNA–miRNA networks related to LGG, in the present study, a comprehensive analysis of the genomic landscape between low-grade gliomas and normal brain tissues from the GEO and TCGA datasets was first conducted to identify differentially expressed genes (DEGs) and differentially expressed miRNAs in LGG. Following a series of analyses, including WGCNA, GO and KEGG analyses, PPI and key model analyses, and survival analysis of the DEGs with clinical phenotypes, the potential key genes were screened and identified, and the related miRNA–mRNA networks were subsequently constructed through miRWalk 3.0. Finally, the potential miRNA–mRNA networks were further validated in CGGA (Chinese Glioma Genome Atlas) datasets and clinical specimens by qRT-PCR.

Results

In our results, six hub genes, MELK, NCAPG, KIF4A, NUSAP1, CEP55, and TOP2A, were ultimately identified. Two regulatory pathways, miR-495-3p-TOP2A and miR-1224-3p-MELK, that regulate the pathogenesis of LGG were ultimately identified. Furthermore, the expression of miR-495-3p-TOP2A and miR-1224-3p-MELK in solid tissues was validated by qRT-PCR.

Conclusion

Our study identified hub genes and related miRNA–mRNA regulatory pathways that contribute to the carcinogenesis of LGG, which may help us reveal the mechanisms underlying the development of LGG.

Identification of LINC02310 as an enhancer in lung adenocarcinoma and investigation of its regulatory network via comprehensive analyses

BACKGROUND

Lung adenocarcinoma (LADC) is a major subtype of non-small cell lung cancer and has one of the highest mortality rates. An increasing number of long non-coding RNAs (LncRNAs) were reported to be associated with the occurrence and progression of LADC. Thus, it is necessary and reasonable to find new prognostic biomarkers for LADC among LncRNAs.

METHODS

Differential expression analysis, survival analysis, PCR experiments and clinical feature analysis were performed to screen out the LncRNA which was significantly related to LADC. Its role in LADC was verified by CCK-8 assay and colony. Furthermore, competing endogenous RNA (ceRNA) regulatory network construction, enrichment analysis and protein-protein interaction (PPI) network construction were performed to investigate the downstream regulatory network of the selected LncRNA.

RESULTS

A total of 2431 differentially expressed LncRNAs (DELncRNAs) and 2227 differentially expressed mRNAs (DEmRNAs) were from The Cancer Genome Atlas database. Survival analysis results indicated that lnc-YARS2-5, lnc-NPR3-2 and LINC02310 were significantly related to overall survival. Their overexpression indicated poor prognostic. PCR experiments and clinical feature analysis suggested that LINC02310 was significantly correlated with TNM-stage and T-stage. CCK-8 assay and colony formation assay demonstrated that LINC02310 acted as an enhancer in LADC. In addition, 3 targeted miRNAs of LINC02310 and 414 downstream DEmRNAs were predicted. The downstream DEmRNAs were then enriched in 405 Gene Ontology terms and 11 Kyoto Encyclopedia of Genes and Genomes pathways, which revealed their potential functions and mechanisms. The PPI network showed the interactions among the downstream DEmRNAs.

CONCLUSIONS

This study verified LINC02310 as an enhancer in LADC and performed comprehensive analyses on its downstream regulatory network, which might benefit LADC prognoses and therapies.

Identification of key differentially expressed mRNAs and microRNAs in non-small cell lung cancer using bioinformatics analysis

Non-small cell lung cancer (NSCLC) is a leading cause of mortality worldwide. However, the pathogenesis of NSCLC remains to be fully elucidated. Therefore, the present study aimed to explore the differential expression of mRNAs and microRNAs (miRNAs/miRs) in NSCLC and to determine how these RNA molecules interact with one another to affect disease progression. Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were identified from the GSE18842, GSE32863 and GSE29250 datasets downloaded from the Gene Expression Omnibus (GEO database). Functional and pathway enrichment analysis were performed based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. STRING, Cytoscape and MCODE were applied to construct a protein-protein interaction (PPI) network and to screen hub genes. The interactions between miRNAs and mRNAs were predicted using miRWalk 3.0 and a miRNA-mRNA regulatory network was constructed. The prognostic value of the identified hub genes was then evaluated via Kaplan-Meier survival analyses using datasets from The Cancer Genome Atlas. A total of 782 DEGs and 46 DEMs were identified from the 3 GEO datasets. The enriched pathways and functions of the DEGs and target genes of the DEMs included osteoclast differentiation, cell adhesion, response to a drug, plasma membrane, extracellular exosome and protein binding. A subnetwork composed of 11 genes was extracted from the PPI network and the genes in this subnetwork were mainly involved in the cell cycle, cell division and DNA replication. A miRNA-gene regulatory network was constructed with 247 miRNA-gene pairs based on 6 DEMs and 210 DEGs. Kaplan-Meier survival analysis indicated that the expression of ubiquitin E2 ligase C, cell division cycle protein 20, DNA topoisomerase IIα, aurora kinase A and B, cyclin B2, maternal embryonic leucine zipper kinase, slit guidance ligand 3, phosphoglucomutase 5, endomucin, cysteine dioxygenase type 1, dihydropyrimidinase-like 2, miR-130b, miR-1181 and miR-127 was significantly associated with overall survival of patients with lung adenocarcinoma. In the present study, a miRNA-mRNA regulatory network in NSCLC was established, which may provide future avenues for scientific exploration and therapeutic targeting of NSCLC.

Tumor shrinkage rate as a potential marker for the prediction of long-term outcome in advanced non-small cell lung cancer treated with first-line tyrosine kinase inhibitors

Context

Tyrosine kinase inhibitors (TKIs) targeting epidermal growth factor receptor (EGFR) play an indispensable role in the treatment of non-small cell lung cancer (NSCLC), leading to a survival major breakthrough, but there remains no uniform standard for predicting the efficacy of TKI therapy.

Aims

We retrospectively reviewed the use of EGFR-TKIs for advanced NSCLC between January 2009 and December 2017 in a hospital, which 169 patients who treated with first-line TKIs were enrolled.

Subjects and Methods

Multiple clinical factors, including histology, age, and sex, were analyzed. We calculated the tumor shrinkage rate (TSR) by measuring the longest diameters of the main mass by computed tomography (CT) before TKI therapy and the first CT after TKI therapy. We evaluated overall survival (OS) and progression-free survival (PFS) after first-line TKI therapy, and we assessed factors predicting survival using the Kaplan-Meier method.

Results

Eligible patients were sorted into higher (n = 83) and lower (n = 86) TSR groups according to the mean TSR of 0.49%. The 83 patients with a higher TSR had longer PFS and OS than those in the 86 patients with a lower TSR (14.83 vs. 8.40 months, P < 0.001, and 31.03 vs. 20.10 months, P < 0.001, respectively). Multivariate analyses revealed that TSR was an independent predictor of PFS and OS (PFS hazard ratio [HR]: 0.506, P < 0.001, and OS HR: 0.291, P < 0.001).

Conclusions

These cumulative data support that TSR may be an early predictor of the treatment efficacy in NSCLC with EGFR mutations treated with first-line TKIs.

Potential Functional Variants in DNA Repair Genes Are Associated with Efficacy and Toxicity of Radiotherapy in Patients with Non-Small-Cell Lung Cancer

Background

Lung cancer is one of the leading causes of cancer-related deaths. Radiotherapy, either alone or with chemotherapy, is still the primary treatment for patients with non-small-cell lung cancer (NSCLC). There are variations in how patients with NSCLC respond to radiotherapy and how toxic the therapy is. DNA repair gene polymorphisms are related to cancer development; however, their association with radiotherapy outcomes remains unknown. We hypothesized that gDNA repair gene variation could affect the efficacy and toxicity of radiotherapy in patients with NSCLC.

Methods

A total of 486 histologically confirmed patients with NSCLC were recruited from the Shengjing Hospital of China Medical University from July 2015 to September 2019. Eleven potentially functional single nucleotide polymorphisms (SNPs) in four DNA repair genes (XRCC1, XRCC2, XPD, and MSH2) were genotyped in these patients. A multiple factor logistic regression analysis was used to assess the association between these SNPs and the efficacy and toxicity of radiotherapy.

Results

Three SNPs, rs25487 (XRCC1), rs3218556 (XRCC2), and rs13181 (XPD), were all significantly associated with the efficacy of radiotherapy. The allele frequencies of the rs25487 CC genotype (OR = 0.457, 95% CI = 0.259–0.804, p=0.006) and the rs3218556 AG or AA genotypes (AG genotype: OR = 0.664, 95% CI = 0.442–0.999, p=0.049; AA genotype: OR = 0.380, 95% CI = 0.181–0.795, p=0.008) were both significantly higher in the response group than in the nonresponse group. For rs13181, the radiotherapy efficacy was associated with the heterozygous genotype GT (OR = 1.663, 95% CI = 1.057–2.614,p=0.027). Statistically significant associations between radiation-induced toxic reactions and rs25487 (XRCC1), rs3218556 (XRCC2), and rs13181 (XPD) were also observed. The rs13181GT genotype was associated with lower toxic reactions than the TT genotype (OR = 1.680, 95% CI = 1.035–2.728,p=0.035).

Conclusions

The variants rs25487 (XRCC1), rs3218556 (XRCC2), and rs13181 (XPD) all contribute to the efficacy and toxicity of radiotherapy in patients with NSCLC. Our findings may clarify the predictive value of DNA repair genes for prognosis in patients with NSCLC after radiotherapy. Further investigation of more genes and samples should be performed to confirm our findings.

Single nucleotide polymorphisms in the dual specificity phosphatase genes and risk of necrotizing enterocolitis in premature infant

BACKGROUND

Differences in the susceptibility of preterm infants to develop necrotizing enterocolitis (NEC) implicate potential genetic differences in response to the inflammatory stimuli leading to NEC. Dual specificity phosphatases (DUSPs) are a key suppressor pathway of the mitogen-activated protein kinase (MAPK) pro-inflammatory signaling pathway. We hypothesized that inherited single nucleotide polymorphisms (SNPs) in DUSP genes contribute to NEC susceptibility in premature infants.

METHODS

Patients admitted between 2010 and 2015 born at <  32 weeks GA and≤1,500 g BW with stage II+NEC (cases; n = 50) and age, weight-matched controls (n = 38) were included. Blood samples were collected for DNA isolation. Agena Mass Array assay was used to examine 31 SNPs in 9 different DUSP genes. Calculated minor allele frequencies (MAF) for cases and controls were compared using χ2 and logistic regression.

RESULTS

The presence of the rs704074 SNP was associated with a 48% decreased risk of developing NEC (OR 0.52; 95% CI 0.27- 1.01, p = 0.04). The odds of surgical NEC decreased by 78% (OR 0.22; 95% CI 0.06- 0.84, p = 0.027) for each copy of rs704074/G allele in patients with NEC.

CONCLUSION

In this small single-center pilot study, DUSP-6 SNP (rs704074) was associated with a lower risk of developing NEC and surgical NEC, the most severe form of NEC, in preterm infants.

Promoter polymorphisms of TOP2A and ERCC1 genes as predictive factors for chemotherapy in non-small cell lung cancer patients

Background

Topoisomerase 2‐alpha (TOP2A) is an enzyme that controls topologic changes in DNA during transcription and replication. ERCC1 is an enzyme that takes part in DNA repair processes. The purpose of this study was to assess the predictive role of particular single nucleotide polymorphisms (SNPs) in the promoter regions of TOP2A and ERCC1 genes in non‐small cell lung cancer patients (NSCLC) treated with chemotherapy.

Materials and methods

We enrolled 113 NSCLC patients treated in the first line with platinum‐based chemotherapy. Effectiveness was available for 71 patients. DNA was isolated from whole blood using the Qiamp DNA Blood Mini kit (Qiagen). We examined five SNPs: rs11615 (ERCC1), rs3212986 (ERCC1), rs13695 (TOP2A), rs34300454 (TOP2A), rs11540720 (TOP2A). Quantitative PCR using the TaqMan probe (ThermoFisher) was performed on a Eco Illumina Real‐Time PCR system device (Illumina Inc).

Results

Patients with the A/A genotype in rs11615 of the ERCC1 gene had significantly longer median progression free survival (PFS) (8.5 months; P = .0088). Patients with the C/C genotype in rs3212986 of the ERCC1 gene had longer median PFS (7 months; P = .05). Patients with the C/C genotype in rs34300454 of TOP2A gene had significantly higher median PFS (7.5 months; P = .0029). Carriers of the C/C genotype in rs34300454 of the TOP2A gene had significantly longer median OS (15.5 months; P = .0017). Patients with the A/A genotype in rs11615 of the ERCC1 gene had significantly higher risk of neutropenia (P = .0133).

Conclusions

Polymorphisms of the TOP2A and ERCC1 genes may be a predictive factor of toxicities and survival for chemotherapy in NSCLC patients.

Src Promotes Metastasis of Human Non-Small Cell Lung Cancer Cells through Fn14-Mediated NF-κB Signaling

Background

Src and Fn14 are implicated in the aggressiveness of non-small cell lung cancer (NSCLC) cells, yet the molecular mechanism is not fully understood.

Material/Methods

The proliferation, migration, and invasion of HCC827 cells with Src knockdown were examined in vitro. The expression of Fn14 and the activation of NF-κB signaling pathway in Src-silenced HCC827 cells were detected by western blot. The role of Fn14 in Src-regulated cell migration/invasion and activation of NF-κB signaling was investigated by overexpressing Fn14 in Src knockdown NSCLC cells. Furthermore, the pro-metastatic role of Src was validated in a NSCLC metastasis mouse model.

Results

Knockdown of Src inhibited the proliferation, migration, and invasion of HCC827 cells, which was associated with reduced levels of Fn14, p-IκBα, p-IKKβ, and nuclear NF-κB p65. Overexpression of Fn14 restored the potential of migration and invasion as well as the activation of NF-κB signaling in Src-silenced NSCLC cells. In addition, silencing of Src suppressed lung metastasis of HCC827 cells in mice, and inhibited the expression of Fn14 and nuclear translocation of NF-κB p65 in vivo.

Conclusions

The data demonstrated that the Src/Fn14/NF-κB axis plays a critical role in NSCLC metastasis.

Over-expression of TOP2A as a prognostic biomarker in patients with glioma

Topoisomerase (DNA) II alpha (TOP2A), an enzyme that controls and alters the topologic states of DNA during transcription, is aberrantly expressed in many cancers. However, few studies have investigated expression of TOP2A and its clinical significance in glioma. We retrieved six independent investigations from the Oncomine database and found that TOP2A is highly expressed in glioma tissues compared with corresponding normal controls. Similar results were also found in clinical specimens at the protein level. Immunohistochemical analysis indicated that TOP2A over expression was highly correlated with grade stage, KI67 positive percentage, IDH1 mutation, and age, but other clinical parameters such as sex distribution and tumor size were barely associated with high TOP2A gene expression. Meanwhile we used Prognos can to assess the prognostic value of TOP2A expression in glioma patients, and found that high expression was associated with poor prognosis of patients with glioma. Furthermore, we used the Gene-Cloud of Biotechnology Information (GCBI) bioinformatics platform predict the role of TOP2A in glioma. It was not only involved in DNA replication, chromosome condensation, and responses to DNA damage stimuli, but also promoted cancer cell mitotic cell cycle and apoptosis, and phosphatidylinositol-mediated signaling by regulating gene expression. By these approaches we demonstrate that TOP2A may be a reliable prognostic factor or therapeutic target in glioma.

High Expression of S100A6 Predicts Unfavorable Prognosis of Lung Squamous Cell Cancer

Background

S100 family of proteins is mainly involved in regulation of intracellular calcium homeostasis. Aberrant expression of S100 family members has been reported in many types of cancers. However, as a member of S100 family, the prognostic value of S100A6 for lung squamous cell carcinoma (SCC) has not been well-studied.

Material/Methods

Using immunohistochemistry, we investigated the expression of S100A6 in 177 patients with SCC and further divided the cohort into a high S100A6 expression group and a low S100A6 expression group. The chi-square test was applied to analyze the correlation between S100A6 expression and clinicopathological factors. Univariate analysis using the Kaplan-Meier method was performed to compare the difference in survival rates between the high S100A6 expression group and the low S100A6 expression group; multivariate analysis with Cox regression model was used to identify independent prognostic risk factors.

Results

In our experiment, we demonstrated that the expression of S100A6 was significantly associated with patient age and tumor differentiation. High-expression of S100A6 was shown to be substantially related to the unfavorable prognosis of SCC. Moreover, our results confirmed that S100A6 was an independent risk factor for SCC prognosis, and could predict unfavorable prognosis.

Conclusions

High-expression of S100A6 was identified as an independent unfavorable prognostic factor for SCC, suggesting that targeting S100A6 may result in the development of potential targeted drug for SCC.