Skip to main content

Elevated expression of MKRN3 in squamous cell carcinoma of the head and neck and its clinical significance

Abstract

Background

Squamous cell carcinoma of the head and neck (SCCHN) is one of the most common types of cancer that cause a substantial number of cancer-related deaths. Our previous study has revealed that makorin ring finger protein 3 (MKRN3) may act as a key regulator of the SCCHN tumorigenesis; however, its specific role in SCCHN progression has not been reported.

Methods

The Cancer Genome Atlas (TCGA) data analysis and quantitative polymerase chain reaction (qPCR) were used to quantify the MKRN3 mRNA expression levels in SCCHN; immunohistochemical staining or immunoblotting analyses were performed to detect MKRN3 protein expression. Kaplan–Meier plotter was used to assess the prognostic values of MKRN3 in terms of overall survival and disease-free survival. The expression differences based on various clinicopathological features were evaluated using subgroup analysis and forest map analysis. The regulatory mechanism of MKRN3 was further investigated using gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Subsequently, STRING was used to perform a co-expression and enrichment analysis for MKRN3. Homologous modeling, molecular docking, and western blot analyses were performed to investigate the relationship between MKRN3 and its potential target gene P53.

Results

MKRN3 was ectopically expressed between cancerous and noncancerous SCCHN tissues, and its expression level was tightly associated with high T classifications as well as advanced clinical stages. qPCR analysis revealed that MKRN3 was upregulated in the SCCHN cell line. Moreover, Kaplan–Meier and Cox regression analyses indicated that SCCHN patients with high MKRN3 expression had poorer prognosis and that MKRN3 was a potential prognostic marker for SCCHN. Using gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses, we determined that MKRN3 may be involved in the regulation of synthesis and metabolism and cell growth, death and motility, as well as cancer pathways associated with SCCHN progression. Mechanism investigation further revealed that P53, a potential target of MKRN3, may be involved in the SCCHN tumorigenesis mediated by MKRN3.

Conclusions

We performed a comprehensive evaluation of the clinical significance of MKRN3 and explored its underlying mechanisms. We concluded that MKRN3 represents a valuable predictive biomarker and potential therapeutic target in SCCHN.

Background

Squamous cell carcinoma of the head and neck (SCCHN) is a histopathological diagnosis that encompasses epithelial malignancies arising from the paranasal sinuses, nasal cavity, oral cavity, pharynx, and larynx [1]. Despite the advances in oncology treatment strategies, limited improvement in the 5-year survival rate of SCCHN patients has been made in the recent decades [2]. The unfavorable outcome has been primarily attributed to late diagnosis, loco-regional recurrences, and cervical lymph node metastasis [3]. Thus, it is important to investigate the pathogenesis of SCCHN and identify molecular biomarkers with prognostic significance in order to enable optimized therapeutic strategies and prolong patient survival when feasible.

In our previous study, we identified MKRN3 as a candidate regulator of SCCHN tumorigenesis using prediction algorithms [4]. MKRN3, an imprinted gene located on the long arm of chromosome 15 (Prader–Willi region), encodes makorin ring finger protein 3, which is involved in the processes of gene transcription and ubiquitination [5, 6]. The MKRN3 protein has four zinc finger domains including three C3H1 motifs and one C3H4 ring finger with presumed E3 ubiquitin ligase activity [7]. It has been initially shown that MKRN3 deficiency causes central precocious puberty in humans [8]. Thus, the biological functions of MKRN3 were further investigated. Currently, MKRN3 is considered a novel imprinted gene involved in the progression of osteosarcoma [9] and non–small cell lung cancer [10]. However, the specific role of MKRN3 in SCCHN is yet to be determined, particularly with respect to its association with clinical outcomes.

Therefore, our study aimed to investigate the clinical relevance of MKRN3 expression in SCCHN. In the present work, the expression of MKRN3 in SCCHN tissue samples and cell lines has been investigated for the first time. Further experiments were performed to assess whether MKRN3 expression is correlated with clinicopathological parameters and to gain insight into the biological pathways and mechanisms regulated by MKRN3 that are involved in SCCHN pathogenesis.

Materials and methods

Data acquisition and processing

A workflow chart of this study is shown in Fig. 1. The expression profiles and clinical information of the 522 SCCHN and 44 adjacent noncancerous epithelial samples were obtained from The Cancer Genome Atlas (TCGA) database (Table 1). The MKRN3 expression levels and overall survival are shown in Additional file 6: Table S6; the main clinical and pathological parameters of SCCHN patients are summarized in Table 1. A P-value < 0.05 was set as the significance threshold for differential expression and clinical outcome. The immunohistochemistry (IHC) data was downloaded from Protein Atlas (https://www.proteinatlas.org). IHC evaluation was based on the staining intensity (0–3) and degree (0–4). The MKRN3 protein expression was classified into two groups: high set (score: 4–7) and low set (score: 0–3). Five fields per IHC slides were randomly selected for evaluation and all IHC were independently and blindly assessed as well as scored by investigators.

Fig. 1
figure1

Study flowchart

Table 1 Clinical characteristics of TCGA SCCHNs patients

Cell culture

SCCHN cell lines were kindly provided by Dr. Joseph Califano (University of California, San Diego, USA) and Dr. Zhuo G. Chen (Emory University Winship Cancer Institute,USA) or purchased from ATCC as previously described [11,12,13]. Human-derived dysplastic oral keratinocytes (DOK), served as a normal cell line, were grown in RPMI 1640 medium (Hyclone, Logan, UT). Tu686 cells were cultured in Dulbecco’s modified Eagle medium (DMEM)/F12 medium (Hyclone). Fadu and CAL27 cells were maintained in DMEM basic medium, whereas JHU011 and Tca8113 cells were kept in RPMI 1640 medium. The culture media were supplemented with 10% fetal bovine serum (FBS) (Gibco, NYC, New York, NY), and all cells were cultured in a humidified incubator at 37 °C and 5% CO2.

SCCHN patient samples

In this study, a total of 25 SCCHN tissues and 23 adjacent normal tissues were collected from patients in Xiangya Hospital, China from September 2015 to December 2017. All samples were collected for quantitative real-time PCR. Additionally, clinicopathologic information of patients that was collected included age at diagnosis, T status, lymph node metastasis, clinical stage, and histology grade. Ages of the 25 SCCHN patients ranged between 38 and 69 years (mean, 56.7 years), and more detailed information are shown in Additional file 5: Table S5. Exclusion criteria for this study were as follows: age < 18 years or age > 75 years; histologically unconfirmed SCCHN; and incomplete clinicopathological parameters. We also excluded patients who had comorbidities and underwent other radiotherapy and chemotherapy strategy.

The study was approved by the Research Ethics Committee of Xiangya Hospital, Central South University, Changsha, China, and all samples were used for analysis under written informed consent from the patients.

Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR)

The total RNA was isolated from SCCHN tissues and cell lines using the TRIzol reagents (Invitrogen, Carlsbad, CA, USA), and cDNA was reverse transcribed using the All-in-One™ mRNA cDNA synthesis Kit (GeneCopoeia, Rockville, MD, USA) following the manufacturers’ protocols. Relative MKRN3 expression values were calculated using the 2CT method and normalized using the GAPDH expression levels [14, 15]. The primers used are listed in Additional file 1: Table S1.

Western blotting assay

Total cell protein was lysed and extracted and separated using 8–12% SDS-PAGE gels and transferred onto PVDF membranes (Millipore, Bedford, MA, USA). Next, membranes were blocked with 5% skimmed milk and incubated with primary antibodies at 4 °C overnight. GAPDH was used as the loading control. The protein expression was imaged using Image Lab 4.1 (Bio-Rad, Hercules, CA, USA) with enhanced chemiluminescence reagents. The relevant antibody information are listed in Additional file 2: Table S2.

Functional enrichment analyses

Protein–protein interaction analysis (PPI), gene ontology (GO) biological process, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were employed in STRING (http://string-db.org) online database to identify over-represented GO terms in biological processes as well as over-represented KEGG pathway terms. The SWISS-MODEL (https://swissmodel.expasy.org) [16], ZDOCK server (https://zdock.umassmed.edu) [17] and PDBePISA Browser (https://www.ebi.ac.uk) [18] were used for homologous modeling, molecular docking of protein and calculate the Minimum Free Energy. For this analysis, a false discovery rate (FDR) < 0.05 and − log FDR > 1.301 were considered to indicate statistical significance.

Statistical analysis

All data were analyzed using SPSS version 22.0 (IBM Corp., Armonk, NY, USA). Results are presented as the mean ± standard deviation. Statistical differences between two groups were determined using Student’s t-test (for equal variances) or Mann–Whitney U test (for unequal variances). Additionally, survival curves were plotted using the Kaplan–Meier method and compared using the log-rank test. P-values < 0.05 were considered statistically significant.

Results

Elevated expression of MKRN3 in SCCHN

TCGA sequencing data were used to delineate differential expression scatterplots and paired difference analyses. As a result, MKRN3 level was markedly higher in SCCHN samples than in corresponding adjacent noncancerous tissues (Fig. 2a; P < 0.01). This was supported by 40 paired cases (Fig. 2b; P < 0.01). Next, the relative expression of MKRN3 was quantified in the 25 SCCHN samples and in 23 samples of adjacent epithelium using quantitative polymerase chain reaction. Our data revealed that the MKRN3 expression was upregulated in most SCCHN cases compared to those observed in the adjacent tissues (Fig. 2c; P < 0.01). Further analysis also showed that the MKRN3 expression level was greater in SCCHN than in DOK cells (Fig. 2d). Finally, we collected the IHC data for MKRN3 expression in SCCHN specimens to investigate whether the MKRN3 protein level was altered. As shown in Fig. 2e, f, the density and intensity of MKRN3 expression in SCCHN tumor tissues were significantly increased compared to those of para-carcinoma tissues. Together, these results suggest that MKRN3 expression is elevated in SCCHN.

Fig. 2
figure2

Elevated expression of MKRN3 in SCCHN. a, b Differential expression of MKRN3 in SCCHN tissue samples according to TCGA data (a) and paired difference analyses (b). c MKRN3 mRNA expression levels in patients with normal tumor tissues. d Relative expression level of MKRN3 mRNA in six different SCCHN cell lines and an immortalized non-malignant cell line derived from oral mucosa (DOK). The gene expression was normalized to GAPDH. e, f Representative immunohistochemical staining shows the expression of the MKRN3 protein in normal (e) and SCCHN tissues (f) according to TCGA (200× , 400×). Data are presented as the mean ± standard deviation. P-values were calculated using the Student’s t-test. *P < 0.05; **P < 0.01; ***P < 0.001. TCGA, The Cancer Genome Atlas

Correlation between clinicopathological features and MKRN3 expression in SCCHN

As summarized in Table 2, a high expression of MKRN3 was closely associated with smoking, advanced clinical stages, and high T classifications in patients with SCCHN (Table 2; all P < 0.05). These data establish the oncogenic role of MKRN3 in SCCHN.

Table 2 Correlations between the expression of MKRN3 and clinicopathological parameters

A high expression level of MKRN3 predicts a worse prognosis in SCCHN patients

The survival analysis indicated that patients with high MKRN3 expression levels had a worse prognosis than those with low MKRN3 expression levels in terms of overall survival (520 cases; Fig. 3a; P < 0.01) and disease-free survival (389 cases; Fig. 3b; P < 0.05). Furthermore, the subgroup analysis revealed that high MKRN3 expression was an unfavorable factor for the prognosis of SCCHN patients at stages III + IV (Fig. 4a; P < 0.001) and T3 + 4 (Fig. 4b; P < 0.05), but not in the N (Fig. 4c, P  > 0.05) or M classification (Fig. 4d, P > 0.05). Univariate Cox regression analyses revealed that age, sex, clinical stage, lymph node metastasis, T classification, and status of MKRN3 expression were significantly associated with overall survival (Table 3; all P < 0.05). Additionally, clinical stage, metastasis, T classification, and status of MKRN3 expression were significantly associated with disease-free survival (Table 4; all P < 0.05). Nevertheless, the multivariate Cox regression analyses showed that metastasis, T classification, and MKRN3 expression level were determined to be independent factors with prognostic value for the overall survival (Table 3) and disease-free survival (Table 4) of SCCHN patients. Additionally, the forest plot analysis of overall survival showed statistically significant associations between age, clinical stage, lymph node metastasis, T classification metastasis, and MKRN3 expression and the outcome of SCCHN patients (Fig. 5; all P < 0.05). Collectively, these findings indicate that MKRN3 represents a valuable biomarker in the surveillance and prognosis in SCCHN patients.

Fig. 3
figure3

High MKRN3 expression levels predicts a worse prognosis in patients with SCCHN. Kaplan–Meier survival analysis of MKRN3 expression in terms of overall survival (a) and disease-free survival (b) in SCCHN patients

Fig. 4
figure4

Correlation between MKRN3 expression and clinicopathologic characteristics. a Subgroup analysis of clinical stage (stages I + II vs. stages III + IV). b T classification (T1 + 2 vs. T3 + T4). c N classification (N0 vs. N +). d M classification (M0 vs. M1). P-values were calculated using the Mann–Whitney U test

Table 3 Cox model analysis of overall survival
Table 4 Cox model analysis of disease-free survival
Fig. 5
figure5

Forest map analysis of MKRN3 expression and clinicopathologic characteristics in SCCHN. P-values were calculated using the Mann–Whitney U test. *P < 0.05; **P < 0.01; ***P < 0.001

Functional analysis of MKRN3

The functional enrichment clustering of MKRN3 showed that a total of 63 categories from GO biological process, such as DNA synthesis and repair, cellular response to tumor microenvironment, regulation of cellular cycle, and translation were identified as important for cancer development (Fig. 6a; Additional file 3: Table S3). KEGG analysis based on these nodes revealed that 17 pathways in cancer, RNA synthesis and metabolism, and cell growth, death, and motility that were significantly enhanced during SCCHN progression were regulated by MKRN3 (Fig. 6b; Additional file 4: Table S4). These factors were closely related to the occurrence and development of cancer, suggesting that MKRN3 expression exhibits a strong relationship with the progress of SCCHN. Further PPI analysis of MKRN3 illustrated that there were 31 nodes based on a combined score ≥ 0.7 in the STRING analysis, and that P53 might be a direct target gene of MKRN3 (Fig. 7a). Therefore, we used homologous modeling and molecular docking for MKRN3 and P53. As shown in Fig. 7b, the ring finger domain of MKRN3 was supposed to form like a “Goldfish”-like shape and tightly bind with the P53, and the Minimum Free Energy value were − 126.6 kcal/mol. Subsequently, western blot analysis showed that MKRN3 substantially repressed the expression of P53 protein (Fig. 7c). Therefore, our data suggest that P53 might be involved in the MKRN3-mediated SCCHN tumorigenesis.

Fig. 6
figure6

Functional enrichment analysis of MKRN3 in SCCHN. a Enriched KEGG biological pathways. b Enriched GO terms in the “biological process” category with an FDR < 0.02. The columns indicate different significance levels, and the curve indicates the number of genes. FDR, false discovery rate; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes

Fig. 7
figure7

P53 might be a target gene of MKRN3. a STRING analysis in the protein–protein interaction of MKRN3. Only the proteins with more than one interaction are displayed. b The homologous modeling and molecular docking with MKRN3 and P53. Red and Blue cartoon represent MKRN3 and P53, respectively. The rectangle highlights the interacted domain. c Relative expression of P53 protein in Tcal8113 and JHU011 cells that transfected with MKRN3 cDNA and normalized for GAPDH

Discussion

In this study, we found that elevated MKRN3 expression was correlated with the aggressive tumor characteristics in SCCHN patients, making MKRN3 an independent prognostic predictor for the survival of SCCHN patients. Further functional analysis of MKRN3 provided insight into the biological pathways and mechanisms involved in SCCHN pathogenesis that are regulated by MKRN3. Our data provide evidence that MKRN3 plays an important role in tumor progression and may serve as a critical promoter of SCCHN.

Previous studies investigating MKRN3 has focused on its role as a ubiquitin E3 ligase during puberty initiation [5]. In 2013, mutations with loss of function in the MKRN3 gene were identified from human families with central precocious puberty using whole-exome sequencing [19]. MKRN3 is a maternally imprinted gene located in the Prader–Willi syndrome critical region (chromosome15q11eq13), and only subjects who inherit the mutation from their father develop central precocious puberty [20]. Although MKRN3 is postulated to be an inhibitor of gonadotropin-releasing hormone secretion [21]; the molecular mechanism through which MKRN3 influences the gonadotropin-releasing hormone network remains unclear. In recent years, the correlation between the change in hormonal regulation and cancer has been investigated. Studies have shown that early age at menarche may act as a cancer promoter in breast carcinomas [22], whereas in other types of cancers, the change in hormonal status can inhibit cancer cells invasion [23]. To our knowledge, there was only one report that has evaluated the prognostic value of sex hormone receptor expression in 50 patients with laryngeal squamous cell carcinoma, whose physiological changes occurred due to rapid development during puberty [24]. However, there are few studies on the correlation between MKRN3 expression and tumorigenesis, especially in SCCHN. Thus far, MKRN3 has only been found as an oncogene associated with gastric cancer [25] and imprinted genes in the process of human osteosarcoma [9]. Nevertheless, no research has provided more pervasive evidence to strengthen the links between MKRN3 and tumorigenesis.

With the rapid development of whole-genome sequencing and tumor databases, large-scale global gene expression profiling and database mining becomes more convenient for identifying a potential correlation between gene expression profiles and overall survival in a variety of malignancies [26], including SCCHN [27,28,29,30]. The clinical significance of MKRN3, especially its prognostic value in SCCHN, was a key highlight of our current investigation. Herein, we provide evidence that MKRN3 plays a key role in SCCHN progression. Our findings demonstrated that the MKRN3 expression level was markedly increased in SCCHN samples compared to the corresponding adjacent noncancerous tissues. Most importantly, Kaplan–Meyer survival and Cox regression analyses based on the expression level of MKRN3 strengthened the notion that MKRN3 represents a valuable prognosis biomarker with predictive potential in patients with SCCHN. Further studies are required to test whether MKRN3 expression in other solid carcinomas with large patient numbers to broaden its clinical significance.

Consistently, through GO and KEGG functional enrichment analyses, we found that a comprehensive molecular mechanism of MKRN3 action in SCCHN included interferon gamma–mediated signaling pathway, hypoxia-inducible factor 1 signaling pathway, DNA/RNA synthesis and metabolism, and cell cycle regulation. These terms were closely related to the occurrence and development of cancer [31,32,33,34,35]. and further mechanism investigation revealed that P53, a tumor suppressor gene [36], might be a direct target gene of MKRN3. thus, suggesting that MKRN3 is an oncogene and could represent a novel targeted therapeutic strategy for treating SCCHN.

In conclusion, the present study revealed that MKRN3 was upregulated in SCCHN tissues, and its expression may represent a potential marker for prognostic evaluation of SCCHN. However, our study has limitations because the results obtained from bioinformatics analysis are insufficient and need to be confirmed via functional experimental and mechanistic exploration. Furthermore, only one sample out of 522 SCCHN exhibited distant metastases. Therefore, further investigation is required to determine whether MKRN3 may represent an intriguing and novel therapeutic target in SCCHN.

Conclusion

In this study, we comprehensively analyzed the expression and prognostic values of MKRN3 in SCCHN. We found that MKRN3 possesses a diagnostic value in SCCHN progress. Molecular mechanisms provided important clues for developing novel therapeutic targets in SCCHN and P53 may represent a potential target gene. Overall, MKRN3 might have applications as a prognosis biomarker with predictive potential and a novel therapeutic target in SCCHN patients. Nevertheless, these findings need to be confirmed in the future.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

SCCHN:

Squamous cell carcinoma of the head and neck

MKRN3:

Makorin ring finger protein 3

TCGA:

The Cancer Genome Atlas

IHC:

Immunohistochemistry

DMEM:

Dulbecco’s modified Eagle medium

GAPDH:

Glyceraldehyde-3-phosphate dehydrogenase

SDS-PAGE:

Sodium dodecyl sulfatepolyacrylamide gel electrophoresis

PPI:

Protein–protein interaction analysis

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

References

  1. 1.

    Chow LQM. Head and neck cancer. N Engl J Med. 2020;382(1):60–72.

    CAS  Article  Google Scholar 

  2. 2.

    Miller KD, Siegel RL, Lin CC, Mariotto AB, Kramer JL, Rowland JH, Stein KD, Alteri R, Jemal A. Cancer treatment and survivorship statistics, 2016. CA Cancer J Clin. 2016;66(4):271–89.

    Article  Google Scholar 

  3. 3.

    Argiris A, Karamouzis MV, Raben D, Ferris RL. Head and neck cancer. Lancet. 2008;371(9625):1695–709.

    CAS  Article  Google Scholar 

  4. 4.

    Zhang S, Li G, Liu C, Lu S, Jing Q, Chen X, Zheng H, Ma H, Zhang D, Ren S, et al. miR-30e-5p represses angiogenesis and metastasis by directly targeting AEG-1 in squamous cell carcinoma of the head and neck. Cancer Sci. 2020;111(2):356–68.

    CAS  Article  Google Scholar 

  5. 5.

    Liu H, Kong X, Chen F. Mkrn3 functions as a novel ubiquitin E3 ligase to inhibit Nptx1 during puberty initiation. Oncotarget. 2017;8(49):85102–9.

    Article  Google Scholar 

  6. 6.

    Abreu AP, Toro CA, Song YB, Navarro VM, Bosch MA, Eren A, Liang JN, Carroll RS, Latronico AC, Rønnekleiv OK, et al. MKRN3 inhibits the reproductive axis through actions in kisspeptin-expressing neurons. J Clin Invest. 2020;130(8):4486–500.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Bulcao Macedo D, Nahime Brito V, Latronico AC. New causes of central precocious puberty: the role of genetic factors. Neuroendocrinology. 2014;100(1):1–8.

    Article  Google Scholar 

  8. 8.

    Abreu AP, Dauber A, Macedo DB, Noel SD, Brito VN, Gill JC, Cukier P, Thompson IR, Navarro VM, Gagliardi PC, et al. Central precocious puberty caused by mutations in the imprinted gene MKRN3. N Engl J Med. 2013;368(26):2467–75.

    CAS  Article  Google Scholar 

  9. 9.

    Li Y, Meng G, Guo QN. Changes in genomic imprinting and gene expression associated with transformation in a model of human osteosarcoma. Exp Mol Pathol. 2008;84(3):234–9.

    CAS  Article  Google Scholar 

  10. 10.

    Li K, Zheng X, Tang H, Zang YS, Zeng C, Liu X, Shen Y, Pang Y, Wang S, Xie F, et al. E3 ligase MKRN3 is a tumor suppressor regulating PABPC1 ubiquitination in non-small cell lung cancer. J Exp Med. 2021. https://doi.org/10.1084/jem.20210151.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Qin Y, Wang J, Zhu G, Li G, Tan H, Chen C, Pi L, She L, Chen X, Wei M, et al. CCL18 promotes the metastasis of squamous cell carcinoma of the head and neck through MTDH-NF-kappaB signalling pathway. J Cell Mol Med. 2019;23(4):2689–701.

    CAS  Article  Google Scholar 

  12. 12.

    Jing Q, Li G, Chen X, Liu C, Lu S, Zheng H, Ma H, Qin Y, Zhang D, Zhang S, et al. Wnt3a promotes radioresistance via autophagy in squamous cell carcinoma of the head and neck. J Cell Mol Med. 2019;23(7):4711–22.

    CAS  Article  Google Scholar 

  13. 13.

    Li G, Liu Y, Su Z, Ren S, Zhu G, Tian Y, Qiu Y. MicroRNA-324-3p regulates nasopharyngeal carcinoma radioresistance by directly targeting WNT2B. Eur J Cancer. 2013;49(11):2596–607.

    CAS  Article  Google Scholar 

  14. 14.

    Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.

    CAS  Article  Google Scholar 

  15. 15.

    Zhang S, He Y, Liu C, Li G, Lu S, Jing Q, Chen X, Ma H, Zhang D, Wang Y, et al. miR-93-5p enhances migration and invasion by targeting RGMB in squamous cell carcinoma of the head and neck. J Cancer. 2020;11(13):3871–81.

    CAS  Article  Google Scholar 

  16. 16.

    Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296-w303.

    CAS  Article  Google Scholar 

  17. 17.

    Pierce BG, Wiehe K, Hwang H, Kim BH, Vreven T, Weng Z. ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinformatics. 2014;30(12):1771–3.

    CAS  Article  Google Scholar 

  18. 18.

    Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372(3):774–97.

    CAS  Article  Google Scholar 

  19. 19.

    Soriano-Guillén L, Argente J. Central precocious puberty, functional and tumor-related. Best Pract Res Clin Endocrinol Metab. 2019;33(3): 101262.

    Article  Google Scholar 

  20. 20.

    Abreu AP, Macedo DB, Brito VN, Kaiser UB, Latronico AC. A new pathway in the control of the initiation of puberty: the MKRN3 gene. J Mol Endocrinol. 2015;54(3):R131-139.

    CAS  Article  Google Scholar 

  21. 21.

    Collaborative Group on Hormonal Factors in Breast Cancer. Menarche, menopause, and breast cancer risk: individual participant meta-analysis, including 118 964 women with breast cancer from 117 epidemiological studies. Lancet Oncol. 2012;13(11):1141–51.

    Article  Google Scholar 

  22. 22.

    Day FR, Thompson DJ, Helgason H, Chasman DI, Finucane H, Sulem P, Ruth KS, Whalen S, Sarkar AK, Albrecht E, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49(6):834–41.

    CAS  Article  Google Scholar 

  23. 23.

    Yager JD. Mechanisms of estrogen carcinogenesis: the role of E2/E1-quinone metabolites suggests new approaches to preventive intervention—a review. Steroids. 2015;99(Pt A):56–60.

    CAS  Article  Google Scholar 

  24. 24.

    Atef A, El-Rashidy MA, Elzayat S, Kabel AM. The prognostic value of sex hormone receptors expression in laryngeal carcinoma. Tissue Cell. 2019;57:84–9.

    CAS  Article  Google Scholar 

  25. 25.

    Xu Y, Sun JY, Jin YF, Yu H. PCAT6 participates in the development of gastric cancer through endogenously competition with microRNA-30. Eur Rev Med Pharmacol Sci. 2018;22(16):5206–13.

    CAS  PubMed  Google Scholar 

  26. 26.

    Zhou L, Huang W, Yu HF, Feng YJ, Teng X. Exploring TCGA database for identification of potential prognostic genes in stomach adenocarcinoma. Cancer Cell Int. 2020;20:264.

    Article  Google Scholar 

  27. 27.

    Wang J, Chen X, Tian Y, Zhu G, Qin Y, Chen X, Pi L, Wei M, Liu G, Li Z, et al. Six-gene signature for predicting survival in patients with head and neck squamous cell carcinoma. Aging. 2020;12(1):767–83.

    CAS  Article  Google Scholar 

  28. 28.

    Wang X, Yin Z, Zhao Y, He M, Dong C, Zhong M. Identifying potential prognostic biomarkers in head and neck cancer based on the analysis of microRNA expression profiles in TCGA database. Mol Med Rep. 2020;21(3):1647–57.

    CAS  PubMed  Google Scholar 

  29. 29.

    He Y, Ji P, Li Y, Wang R, Ma H, Yuan H. Genetic variants were associated with the prognosis of head and neck squamous carcinoma. Front Oncol. 2020;10:372.

    Article  Google Scholar 

  30. 30.

    Tao W, Tian G, Xu S, Li J, Zhang Z, Li J. NAT10 as a potential prognostic biomarker and therapeutic target for HNSCC. Cancer Cell Int. 2021;21(1):413.

    CAS  Article  Google Scholar 

  31. 31.

    Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, Dimitriadoy S, Liu DL, Kantheti HS, Saghafinia S, et al. Oncogenic signaling pathways in the cancer genome atlas. Cell. 2018;173(2):321-337.e310.

    CAS  Article  Google Scholar 

  32. 32.

    Hiroi M, Mori K, Sekine K, Sakaeda Y, Shimada J, Ohmori Y. Mechanisms of resistance to interferon-gamma-mediated cell growth arrest in human oral squamous carcinoma cells. J Biol Chem. 2009;284(37):24869–80.

    CAS  Article  Google Scholar 

  33. 33.

    Lehmann S, Te Boekhorst V, Odenthal J, Bianchi R, van Helvert S, Ikenberg K, Ilina O, Stoma S, Xandry J, Jiang L, et al. Hypoxia induces a HIF-1-dependent transition from collective-to-amoeboid dissemination in epithelial cancer cells. Curr Biol. 2017;27(3):392–400.

    CAS  Article  Google Scholar 

  34. 34.

    De Schutter H, Kimpe M, Isebaert S, Nuyts S. A systematic assessment of radiation dose enhancement by 5-Aza-2ʹ-deoxycytidine and histone deacetylase inhibitors in head-and-neck squamous cell carcinoma. Int J Radiat Oncol Biol Phys. 2009;73(3):904–12.

    Article  Google Scholar 

  35. 35.

    Zhang C, Zhang B, Meng D, Ge C. Comprehensive analysis of DNA methylation and gene expression profiles in cholangiocarcinoma. Cancer Cell Int. 2019;19:352.

    Article  Google Scholar 

  36. 36.

    Jiang L, Kon N, Li T, Wang SJ, Su T, Hibshoosh H, Baer R, Gu W. Ferroptosis as a p53-mediated activity during tumour suppression. Nature. 2015;520(7545):57–62.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (nos. 81773243 81874133, 82103641 and 82103631), the National Key Research and Development Program of China (2016YFC0902003), and the Natural Science Foundation of Hunan Province (nos. 2018JJ2630 and 2017JJ3488).

Author information

Affiliations

Authors

Contributions

SZ: development of methodology and acquisition, analysis and manuscript writing. CL and GL: acquisition and interpretation of data (provided TCGA data sets, acquired SCCHN cells, etc.). YL, XW and YQ: conception, design, study supervision and editing the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xingwei Wang or Yuanzheng Qiu.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Research Ethics Committee of Xiangya Hospital, Central South University, Changsha, China, and all samples were used for analysis under written informed consent from the patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. The sequence of primers.

Additional file 2

: Table S2. The antibody information.

Additional file 3

: Table S3. The GO biological process.

Additional file 4

: Table S4. The KEGG pathways.

Additional file 5

: Table S5. The clinical characteristics of 25 SCCHN patients.

Additional file 6

: Table S6. The data of MKRN3 expression and overall survival.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, S., Liu, C., Li, G. et al. Elevated expression of MKRN3 in squamous cell carcinoma of the head and neck and its clinical significance. Cancer Cell Int 21, 557 (2021). https://doi.org/10.1186/s12935-021-02271-6

Download citation

Keywords

  • Squamous cell carcinoma of the head and neck
  • MKRN3
  • P53
  • Tumorigenesis
  • Prognostic factor