Transcription factor FoxM1 is the downstream target of c-Myc and contributes to the development of prostate cancer
© The Author(s). 2018
Received: 15 November 2017
Accepted: 27 February 2018
Published: 20 March 2018
Prostate cancer is a common malignancy and the second leading cause of cancer death in men. Elevated expression of the transcription factor FoxM1 and c-Myc has been identified in prostate cancer. However, the potential mechanism of elevated FoxM1 and c-Myc to the development of prostate cancer has not been identified.
In this report, the mRNA level of FoxM1 and c-Myc was detected in 30 prostate cancer and para-cancer tissues. Then, we detected the expression level of FoxM1 by real-time PCR and Western blot after disturbance of the expression level of c-Myc in PC-3 cells. Whether c-Myc could bind to FoxM1 promoter was identified by ChIP assay. Finally, the migratory, invasive, and proliferative abilities in FoxM1 overexpressing and silencing PC-3 cells were detected by wound healing, transwell assay, CCK-8 assays, and Ki-67 protein level.
We found that the expression level of FoxM1 and c-Myc were both increased in prostate cancer samples compared with para-cancer samples. The expression level of FoxM1 was changed consistent with the protein level of c-Myc. ChIP assay detected the direct binding of c-Myc in FoxM1 gene promoter. Lastly, overexpression of FoxM1 increased the migratory, invasive, and proliferative abilities of PC-3 cells, and its downregulation significantly decreased the migratory, invasive, and proliferative abilities.
In conclusion, FoxM1 was significantly increased in prostate cancer samples, and it could regulate the proliferative and invasive ability of prostate cancer cells which might be a new target for prostate cancer. Besides, c-Myc could regulate the expression level of FoxM1 by directly binding to its gene promoter.
Prostate cancer is one of the common malignant tumors in western countries . At present, the incidence of prostate cancer ranks the second highest of all malignancies in the world and the mortality rate ranks fifth , although the incidence and mortality of prostate cancer in China is much lower than other developed countries. In recent years, the number of persons who diagnosed with prostate cancer in China has obviously increased with the aging of the population and the change of people’s eating habits. . So, studying the molecular mechanism of prostate cancer is important to its identification and treatment.
The dysregulation of tumor suppressor genes and oncogenes in cancer cells is the most common event in the pathogenesis of prostate cancer . These changes could be the result and at the same time be the cause of tumor. Studying the function of those genes is still important to clarify the mechanism of the pathogenesis of tumor. FoxM1 is an important member of the FOX gene family, which is only specifically expressed in the proliferative stage cells . It can also promote the transition from cell cycle G1 phase to S phase and G2 phase to M phase by inducing the expression of cell cycle-related factors [6, 7]. So, it is suggested that FoxM1 gene could regulate the cell proliferation ability. FoxM1 could also regulate cell differentiation and apoptosis [8–10]. Recent studies have confirmed that abnormal expression of FoxM1 is highly associated with the generation and progression of various human malignancies including breast cancer, liver cancer, lung cancer, and bladder cancer [11–14]. Its expression level is often associated with the degree of malignancy which indicated its potential to be used as a new target for clinical cancer therapy. However, the expression level of FoxM1 in prostate cancer and its clinical significance in the progression of prostate cancer have not been fully studied.
c-Myc, a member of the MYC gene family, is involved in regulating a lot of biological activities . From now on, c-Myc has been found overexpressed in colon cancer, breast cancer, lung cancer, and prostate cancer [16–19]. It is a bidirectional regulatory gene in promoting cell proliferation and inducing cell apoptosis. c-Myc could promote proliferative ability of cells in the stimulation of hematopoietic growth factor. In the absence of hematopoietic growth factor, c-Myc could mediate apoptosis. c-Myc gene also belongs to the apoptosis regulatory gene . In addition, c-Myc could also regulate cell cycle and cell metabolism . In the quiescent cells, the expression level of c-Myc is very low. It accumulates as the initial response gene and maintains at high level throughout the cell cycle when stimulated with growth factors. Once the expression level of c-Myc is out of control, it can lead to tumorigenesis. c-Myc regulates its target genes by directly binding or recruiting histone modification enzyme to the gene’s promoter region [22, 23]. The predicted c-Myc binding sites in FoxM1’s promoter suggested its regulatory role for FoxM1 expression.
In this study, the increased mRNA level of FoxM1 and c-Myc in prostate cancer tissues was further confirmed. Then, we mostly focused on the mechanism of the increased FoxM1 expression level to the pathogenesis of prostate cancer and the regulatory factors of FoxM1. According to ChIP assay and real-time PCR assay, we proved FoxM1 was regulated by c-Myc. Further experiments showed that FoxM1 could promote the proliferative and invasive ability of prostate cancer cells. Altogether, the increased expression of FoxM1 in prostate tissues might be caused by c-Myc overexpression, and FoxM1 promoted the development of prostate cancer.
All tissues were obtained from patients who underwent surgical procedures for prostate cancer treatment or diagnoses in Ningbo No.2 Hospital. In total, 30 prostate cancer tissues and 30 para-cancer tissues were obtained for experiments in this study. All the tissues were collected for RNA isolation.
Human PC-3 cells were purchased from American Type Culture Collection. All the cells were cultured in 1640 media supplemented with 10% fetal bovine serum in a 5% CO2 humidified atmosphere at 37 °C.
RNA isolation and real-time PCR
After PC-3 cells were transfected with c-Myc overexpression plasmid and small interference RNAs (the small interference RNA of c-Myc was synthesized by RiboBio), the cells were collected by centrifugation and discarded the supernatant. According to the operation instructions, the RNA was extracted by Trizol reagent. After extraction, the RNA was quantified and purified by ultraviolet spectrophotometer. The integrity of RNA was identified by 0.8% agarose gel electrophoresis. Then, 1 μg of total RNA was reversed by ABI reverse transcriptional kit according to the manufacturer’s instruction. Reverse transcription reaction procedure was as follows: 25 °C for 10 min, 37 °C for 120 min, and 85 °C for 5 min. Lastly, the real-time PCR amplification scheme was carried out according to the manufacturer’s instructions of ROCH. The primers used to detect the expression of FoxM1 and c-Myc were as follows: FoxM1 forward: GGAGCAGCGACAGGTTAAGG, reverse: GTTGATGGCGAATTGTATCATGG c-Myc forward: GGCTCCTGGCAAAAGGTCA; reverse: CTGCGTAGTTGTGCTGATGT.
After PC-3 cells were transfected with c-Myc and FoxM1 overexpression plasmid and small interference RNAs, cells were collected by centrifugation to discard the supernatant. The cells are collected by centrifugation and discarded. The cell deposition was re-suspended in cell lysis buffer (50 mM Tris, pH 7.4, 0.5% NP-40 and 0.01% SDS) containing protease inhibitors. After 5 min of boiling in SDS-loading buffer, the supernatant was separated by sodium dodecyl sulfate (SDS)-polyacrylamide gel. Then, the protein in the gel was transferred to the PVDF membrane. In order to prevent the non-specific protein and protein interactions, the PVDF membrane was blocked at room temperature in Tris saline Tween (TST) buffer containing 5% non-fat milk. The antibodies were diluted in TST buffer (c-Myc, 1:4000 (Santa Cruz,); FoxM1, 1:2000 (Santa Cruz,); β-Actin, 1:5000 (Santa Cruz)); PVDF membranes were maintained TST buffer containing diluted antibody and shaked at 4 °C overnight. The PVDF membrane was washed with TST buffer the next day, and then, a secondary antibody with peroxidase conjugation was incubated at room temperature for 1 h. Finally, ECL solution was prepared in a dark room, and exposure time was measured according to fluorescence intensity.
PC-3 cells were grown into 96-well cell culture plates for 5 × 105 cells each well. After cells were transfected with c-Myc and FoxM1 overexpression plasmid and small interference RNAs, PC-3 cells were harvested in medium 24 and 48 h after normal culture for proliferation detection. Cell proliferation ability was detected according to the manufacturer’s instruction of CCK8 assay kit. Absorbance was detected using a microplate reader at the wavelength of 450 nm.
Wound healing assay
PC-3 cells were transfected with FoxM1 overexpression plasmid and small interference RNA on the second day of the plate. After 6 h transfection, the sterile 10-μl white tip was used to draw a straight line from top to bottom at the bottom of a 3.5-cm dish and replaced with serum-free medium. At the same time, pictures were taken as 0 h in × 40 microscope. Following culture for 24 and 48 h, pictures were taken as 24 h in × 40 microscope.
The invasive ability of PC-3 cells was detected by using Matrigel transmembrane invasion assay. The transwell chambers were coated with Matrigel. PC-3 cells were transfected either with FoxM1 overexpressing plasmid or with FoxM1 siRNA on the second day of the plate. After transfection for 24 h, PC-3 cells were harvested for calculating and plating into the upper chamber and bottom wells were filled with complete medium. After incubation for 24 h, scraping with a cotton swab was used to remove cells from the upper surface of the filter. The invaded cells were stained with crystal violet solution after being fixed with methanol. Lastly, the number of penetrated cells of five randomly selected fields were counted, then the mean number was calculated.
ChIP assays were conducted according to the manufacturer (Upstate) of ChIP assay kit. For quantitative ChIP analysis, real-time PCR was carried out by applying SYBR Green PCR Master Mix according to the instructions of the Takara manufacturer. Then, 1:4, 1:16, 1:16, and 1:25 diluted input samples were used to make the relative quantitative standard curves. All PCR results from non-precipitate samples were transferred into their respective input standard curve to normalize the number of cells and primers efficiency. The procedure of real-time PCR is applied according to the method described in the report. Δ Ct value in the first uses the following formula: Δ Ct = Ct (sample) − Ct (input), Δ Δ Ct (Δ Ct sample − Δ Ct negative control), and fold difference (2 (− Δ Δ Ct treatment)/2 (− Δ Δ Ct control) was derived. Primers used in ChIP assay are as follows: F1: ATCCAGCCTTTGCACGCCTGA; R1: AGAGGTTAATGATACTTTGAGT; F2: TGGAGACGGTGTCTCCCTCTGT; R2: AGCGAAACCGCTTCTCTACTA.
All the results in this report were showed as the mean value ± standard deviation (SD). Student’s t test was used for all the statistical analyses. Each experiment in this research was repeated at least in triplicate. Each p value less than 0.05 is thought to be with great significance.
The expression level of FoxM1 and c-Myc was increased in prostate cancer tissues
FoxM1 was regulated by c-Myc
c-Myc bound to FoxM1 gene promoter
FoxM1 regulated the proliferative ability of PC-3 cells
FoxM1 regulated the invasive and migratory ability of PC-3 cells
Abnormal proliferative and invasive abilities are the phenotypes of cancer cells; we found the increased proliferative and invasive abilities in prostate cancers were possibly caused by FoxM1. In this report, we found increased mRNA level of FoxM1 was caused by c-Myc which was a new mechanism in prostate cancer.
We found the expression level of FoxM1 was increased in prostate cancer tissues which was consistent with the reports of the other cancers. FoxM1 was an important regulator to the proliferative and invasive ability of prostate cancer cells which contributed to the development of prostate cancer. As an important transcription factor, FoxM1 could regulate a wide range of other biological processes, including differentiation, metabolism, and apoptosis. In the study of FoxM1 in gastric cancer, FoxM1 could inhibit cell apoptosis and promote proliferation of gastric cells . FoxM1 also played a vital role in promoting angiogenesis and lymph node metastasis by regulating the expression of MMP in gastric cancer . Besides, the researchers found FoxM1 could participate in regulating the telomerase activity in gastric cancer cells . In hepatocellular carcinoma, FoxM1 could regulate cell cycle arrest, and the expression of cyclin B1, p27, and cyclin D1 was all changed in FoxM1 knock out MHCC-97H cells . So, studying the other functions of FoxM1 in prostate cancer was important to fully understand its role in the progression of prostate cancer.
c-Myc could regulate the expression level of FoxM1 by binding to its promoter region which was a newfound mechanism in prostate cancer. Besides, its regulatory mechanism includes the indirect regulatory mechanism mediated by other proteins . For example, c-Myc could regulate its target genes by interacting with transformation/transcription domain-associated protein (TRRAP) and other histone acetyltransferase . The specific mechanism of c-Myc in regulating FoxM1 by directly binding or being assisted by other proteins still needs further investigation. Except transcription activation, the other studies found that c-Myc possesses transcriptional inhibition activity . It depends on the sequence nearby c-Myc binding sites. Both the increased expression of FoxM1 and c-Myc would promote the development of prostate cancer, and they might contribute to this disease by different mechanisms. In this work, c-Myc upregulated the expression of FoxM1 to promote the cancer development which might be used as a target in clinical treatment.
The regulatory mechanism of FoxM1 was complex; c-Myc was one of the important factors to the increased expression level of FoxM1 in PC-3 cells. FoxM1 was regulated by a wide variety of transcription factors, enzyme, and microRNAs. From the previous report, TNF-α could induce the upregulation of FoxM1 in hepatocellular cancers . In breast cancer, FoxM1 was the downstream target of HER2, then contributes to the pathogenesis of breast cancer . In cervical cancer, E6 oncoprotein augmented the expression of FoxM1 through MZF1/NKX2 . All these FoxM1 regulatory mechanisms in prostate cancer need further investigation.
In conclusion, we found increased expression of FoxM1 in prostate cancer samples compared with normal prostate samples. FoxM1 contributed to the pathogenesis of prostate cancer by regulating prostate cancer cells’ proliferative, invasive, and migratory abilities. It might be used as a novel target for clinical treatment of prostate cancer. Besides, we found a new regulatory mechanism of FoxM1 in prostate cancer that was c-Myc could regulate the expression of FoxM1 by directly binding to its promoter.
This research was supported by the hospital technology foundation of the Ningbo No.2 Hospital.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
XR and HP designed and initiated this study. HP, YZ, and WW carried out the experiments of cell culture. SS performed RT-PCR analyses in prostate tissues. XR and HP wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The Ethics Committee of the Ningbo No.2 Hospital approved this study. The analysis was performed in accordance with the ethical standards of the hospital and the tenets of the Declaration of Helsinki/Declaration of Istanbul. The patient reported in this study provided written consent.
Consent for publication
The consent for publication of the individual patient data from the patient has been obtained by Ningbo No.2 Hospital.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Roshandel G, Boreiri M, Sadjadi A, Malekzadeh R. A diversity of cancer incidence and mortality in West Asian populations. Ann Glob Health. 2014;80:346–57.View ArticlePubMedGoogle Scholar
- Duong Q, Hill CL Jr, Janitz AE, Campbell JE. Trends in lung and bronchus, prostate, female breast, and colon and rectum cancers incidence and mortality in Oklahoma and the United States from 1999 to 2012. J Okla State Med Assoc. 2016;109:347–53.PubMedPubMed CentralGoogle Scholar
- Zhao R, Cheng G, Wang B, Qin C, Liu Y, Pan Y, Wang J, Hua L, Zhu W, Wang Z. BMI and serum lipid parameters predict increasing risk and aggressive prostate cancer in Chinese people. Oncotarget. 2017;8:66051–60.PubMedPubMed CentralGoogle Scholar
- Turner DP, Watson DK. ETS transcription factors: oncogenes and tumor suppressor genes as therapeutic targets for prostate cancer. Expert Rev Anticancer Ther. 2008;8:33–42.View ArticlePubMedGoogle Scholar
- Wang X, Quail E, Hung NJ, Tan Y, Ye H, Costa RH. Increased levels of forkhead box M1B transcription factor in transgenic mouse hepatocytes prevent age-related proliferation defects in regenerating liver. Proc Natl Acad Sci U S A. 2001;98:11468–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Wei Y, Sun Q, Zhao L, Wu J, Chen X, Wang Y, Zang W, Zhao G. LncRNA UCA1-miR-507-FOXM1 axis is involved in cell proliferation, invasion and G0/G1 cell cycle arrest in melanoma. Med Oncol. 2016;33:88.View ArticlePubMedGoogle Scholar
- Kelleher FC, O'Sullivan H. FOXM1 in sarcoma: role in cell cycle, pluripotency genes and stem cell pathways. Oncotarget. 2016;7:42792–804.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Liu Y, Ni H, Ding C, Zhang X, Zhang Z. FoxM1 overexpression promotes cell proliferation and migration and inhibits apoptosis in hypopharyngeal squamous cell carcinoma resulting in poor clinical prognosis. Int J Oncol. 2017;51:1045–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Wen N, Wang Y, Wen L, Zhao SH, Ai ZH, Wang Y, Wu B, Lu HX, Yang H, Liu WC, Li Y. Overexpression of FOXM1 predicts poor prognosis and promotes cancer cell proliferation, migration and invasion in epithelial ovarian cancer. J Transl Med. 2014;12:134.View ArticlePubMedPubMed CentralGoogle Scholar
- Ustiyan V, Wert SE, Ikegami M, Wang IC, Kalin TV, Whitsett JA, Kalinichenko VV. Foxm1 transcription factor is critical for proliferation and differentiation of Clara cells during development of conducting airways. Dev Biol. 2012;370:198–212.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu CP, Yu S, Shi L, Wang S, Li ZX, Wang YH, Sun CJ, Liang J. FoxM1 promotes epithelial-mesenchymal transition of hepatocellular carcinoma by targeting Snai1. Mol Med Rep. 2017;16:5181–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Inoguchi S, Seki N, Chiyomaru T, Ishihara T, Matsushita R, Mataki H, Itesako T, Tatarano S, Yoshino H, Goto Y, Nishikawa R, Nakagawa M, Enokida H. Tumour-suppressive microRNA-24-1 inhibits cancer cell proliferation through targeting FOXM1 in bladder cancer. FEBS Lett. 2014;588:3170–9.View ArticlePubMedGoogle Scholar
- Wei P, Zhang N, Wang Y, Li D, Wang L, Sun X, Shen C, Yang Y, Zhou X, Du X. FOXM1 promotes lung adenocarcinoma invasion and metastasis by upregulating SNAIL. Int J Biol Sci. 2015;11:186–98.View ArticlePubMedPubMed CentralGoogle Scholar
- Abdeljaoued S, Bettaieb I, Nasri M, Adouni O, Goucha A, El Amine O, Boussen H, Rahal K, Gamoudi A. Overexpression of FOXM1 is a potential prognostic marker in male breast cancer. Oncol Res Treat. 2017;40:167–72.View ArticlePubMedGoogle Scholar
- Aprelikova O, Chen K, El Touny LH, Brignatz-Guittard C, Han J, Qiu T, Yang HH, Lee MP, Zhu M, Green JE. The epigenetic modifier JMJD6 is amplified in mammary tumors and cooperates with c-Myc to enhance cellular transformation, tumor progression, and metastasis. Clin Epigenetics. 2016;8:38.View ArticlePubMedPubMed CentralGoogle Scholar
- Kugimiya N, Nishimoto A, Hosoyama T, Ueno K, Enoki T, Li TS, Hamano K. The c-MYC-ABCB5 axis plays a pivotal role in 5-fluorouracil resistance in human colon cancer cells. J Cell Mol Med. 2015;19:1569–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Jia Y, Zhao S, Zhang X, Wang X, Han X, Wang Y, Ma M, Shi J, Liu L. BIN1 reverses PD-L1-mediated immune escape by inactivating the c-MYC and EGFR/MAPK signaling pathways in non-small cell lung cancer. Oncogene. 2017;Google Scholar
- Sadeghi S, Hojati Z, Tabatabaeian H. Cooverexpression of EpCAM and c-Myc genes in malignant breast tumours. J Genet. 2017;96:109–18.View ArticlePubMedGoogle Scholar
- Li E, Liu L, Li F, Luo L, Zhao S, Wang J, Kang R, Luo J and Zhao Z. PSCA promotes prostate cancer proliferation and cell-cycle progression by up-regulating c-Myc. Prostate. 2017;77(16):1563–1572.Google Scholar
- Shortt J, Johnstone RW. Oncogenes in cell survival and cell death. Cold Spring Harb Perspect Biol. 2012;4Google Scholar
- Polioudakis D, Bhinge AA, Killion PJ, Lee BK, Abell NS, Iyer VR. A Myc-microRNA network promotes exit from quiescence by suppressing the interferon response and cell-cycle arrest genes. Nucleic Acids Res. 2013;41:2239–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Iraci N, Gherardi S, Gamble LD, Wood KM, Perini G, Lunec J, Tweddle DA. p53 is a direct transcriptional target of MYCN in neuroblastoma. Cancer Res. 2010;70:1377–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Tikhanovich I, Zhao J, Bridges B, Kumer S, Roberts B, Weinman SA. Arginine methylation regulates c-Myc-dependent transcription by altering promoter recruitment of the acetyltransferase p300. J Biol Chem. 2017;292:13333–44.View ArticlePubMedGoogle Scholar
- Okada K, Fujiwara Y, Takahashi T, Nakamura Y, Takiguchi S, Nakajima K, Miyata H, Yamasaki M, Kurokawa Y, Mori M, Doki Y. Overexpression of forkhead box M1 transcription factor (FOXM1) is a potential prognostic marker and enhances chemoresistance for docetaxel in gastric cancer. Ann Surg Oncol. 2013;20:1035–43.View ArticlePubMedGoogle Scholar
- Yu J, Wang X, Li Y, Tang B. Tanshinone IIA suppresses gastric cancer cell proliferation and migration by downregulation of FOXM1. Oncol Rep. 2017;37:1394–400.View ArticlePubMedPubMed CentralGoogle Scholar
- Zeng J, Wang L, Li Q, Li W, Bjorkholm M, Jia J, Xu D. FoxM1 is up-regulated in gastric cancer and its inhibition leads to cellular senescence, partially dependent on p27 kip1. J Pathol. 2009;218:419–27.View ArticlePubMedGoogle Scholar
- Wu QF, Liu C, Tai MH, Liu D, Lei L, Wang RT, Tian M, Lu Y. Knockdown of FoxM1 by siRNA interference decreases cell proliferation, induces cell cycle arrest and inhibits cell invasion in MHCC-97H cells in vitro. Acta Pharmacol Sin. 2010;31:361–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Dang CV. c-Myc target genes involved in cell growth, apoptosis, and metabolism. Mol Cell Biol. 1999;19:1–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu X, Tesfai J, Evrard YA, Dent SY, Martinez E. c-Myc transformation domain recruits the human STAGA complex and requires TRRAP and GCN5 acetylase activity for transcription activation. J Biol Chem. 2003;278:20405–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Xia L, Mo P, Huang W, Zhang L, Wang Y, Zhu H, Tian D, Liu J, Chen Z, Zhang Y, Chen Z, Hu H, Fan D, Nie Y, Wu K. The TNF-alpha/ROS/HIF-1-induced upregulation of FoxMI expression promotes HCC proliferation and resistance to apoptosis. Carcinogenesis. 2012;33:2250–9.View ArticlePubMedGoogle Scholar
- Francis RE, Myatt SS, Krol J, Hartman J, Peck B, McGovern UB, Wang J, Guest SK, Filipovic A, Gojis O, Palmieri C, Peston D, Shousha S, Yu Q, Sicinski P, Coombes RC, Lam EW. FoxM1 is a downstream target and marker of HER2 overexpression in breast cancer. Int J Oncol. 2009;35:57–68.PubMedPubMed CentralGoogle Scholar
- Chen PM, Cheng YW, Wang YC, Wu TC, Chen CY, Lee H. Up-regulation of FOXM1 by E6 oncoprotein through the MZF1/NKX2-1 axis is required for human papillomavirus-associated tumorigenesis. Neoplasia. 2014;16:961–71.View ArticlePubMedPubMed CentralGoogle Scholar