MiR-204 silencing in intraepithelial to invasive cutaneous squamous cell carcinoma progression
© The Author(s). 2016
Received: 25 November 2015
Accepted: 29 June 2016
Published: 25 July 2016
Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer and frequently progresses from an actinic keratosis (AK), a sun-induced keratinocyte intraepithelial neoplasia (KIN). Epigenetic mechanisms involved in the phenomenon of progression from AK to cSCC remain to be elicited.
Expression of microRNAs in sun-exposed skin, AK and cSCC was analysed by Agilent microarrays. DNA methylation of miR-204 promoter was determined by bisulphite treatment and pyrosequencing. Identification of miR-204 targets and pathways was accomplished in HaCat cells. Immunofluorescence and immunohistochemistry were used to analyze STAT3 activation and PTPN11 expression in human biopsies.
cSCCs display a marked downregulation of miR-204 expression when compared to AK. DNA methylation of miR-204 promoter was identified as one of the repressive mechanisms that accounts for miR-204 silencing in cSCC. In HaCaT cells miR-204 inhibits STAT3 and favours the MAPK signaling pathway, likely acting through PTPN11, a nuclear tyrosine phosphatase that is a direct miR-204 target. In non-peritumoral AK lesions, activated STAT3, as detected by pY705-STAT3 immunofluorescence, is retained in the membrane and cytoplasm compartments, whereas AK lesions adjacent to cSCCs display activated STAT3 in the nuclei.
Our data suggest that miR-204 may act as a “rheostat” that controls the signalling towards the MAPK pathway or the STAT3 pathway in the progression from AK to cSCC.
KeywordsmiR-204 Cutaneous squamous cell carcinoma Actinic keratosis Sun-induced keratinocyte intraepithelial neoplasia DNA methylation MAPK STAT3
Cutaneous squamous cell carcinoma (cSCC) is an epidermal keratinocyte derived skin malignant tumor with an incidence of 16 per 100,000 people in Europe . Chronic solar radiation exposure is the most important risk factor for cSCC development since cSCC preferentially arises from precancerous lesions such as UV-induced actinic keratosis (AK) . AK is a keratinocyte intraepithelial neoplasm (KIN or in situ cSCC) that may progress to invasive cSCC at a rate estimated between 0.025 and 16 % for an individual lesion per year , and patients with multiple AKs have an annual risk of developing invasive cSCC ranging from 0.15 to 80 % . While most cSCCs are cured by surgery and/or radiotherapy, a subset of tumors may show an aggressive clinical behaviour, with increased local recurrences and metastases . Patients with metastatic cSCC have a disease-specific survival at 1 year of 44–56 % .
The genetic basis for the progression from actinic keratosis to cSCC is not completely understood [4, 7–10]. TP53 mutations are an early phenomenon already seen in 40 % of in situ cSCC . NOTCH inactivating mutations have recently been found in ~75 % of cSCC  but information about their presence in AK is lacking. Although increased genomic instability is found in cSCC, AKs already harbour cytogenetic alterations at several loci . We previously demonstrated that EGFR and MYC gains are already frequent in AKs and not significantly increased in invasive cSCC [14, 15]. Thus, further genomic and epigenetic studies are needed to better understand the progression from AK to cSCC.
MicroRNAs (miRNAs) are a class of small non-coding RNAs that can regulate gene expression of a very broad set of targets including loci in the protein-coding region of mRNAs, 5′ UTRs, intronic and intergenic transcripts and other non–protein-coding RNAs . MiRNAs are proposed to regulate 60 % of all protein-coding genes in humans . Every miRNA can have multifunctional roles as either repressors or enhancers of their targets expression by diverse pathways .
Aberrant miRNA expression is linked to cancer development and progression and affects several processes including proliferation, apoptosis, differentiation and invasiveness. In cSCC most of the altered miRNAs are downregulated (miR-125b, miR-34a, miR-214, miR-124, miR-361, miR-193b, miR-365a, miR-20a, miR-199a) [19–25] whereas only a small number of miRNAs have been found to be upregulated and act as onco-miRNAs, being involved in angiogenesis, colony formation, migration and invasion, and metastatic spread (miR-365, miR-9, miR-184, miR-21, miR-31, miR-135b, miR-205, let-7a) [25–34]. However, most of these miRNAs have been identified by comparing miRNA expression in cSCC versus healthy skin, and only few studies were aimed to compare miRNA expression in AKs and cSCC [19, 23, 31]. Here we explore the miRNA expression signature of cSCC and AK to identify miRNAs implicated in the development and evolution from AK to cSCC.
All samples were obtained at the Department of Dermatology of Hospital del Mar (Barcelona, Spain) and immediately embedded in OCT and frozen in liquid nitrogen.
Four millimeter punch biopsies were taken from moderately differentiated cSCCs (20 lesions in sun-exposed areas from 20 patients). Five KIN3 actinic keratoses and 5 control chronically sun-exposed (CSE) skin samples were obtained by shave biopsy from 10 different patients. Clinicopathological data of the patients can be found in Additional file 1: Figure S1. Samples were macrodissected from frozen blocks after examination of hematoxylin-eosin stained sections such that all samples displayed a minimum 70 % enrichment for tumor or dysplastic cells. RNA quality and concentration was measured using the Agilent Bioanalyser. Only samples with RIN ≥6 (Bioanalyzer 2100, Agilent Technologies) were included in the study. RNA (10 ng) was used as template for whole transcriptome amplification and cDNA synthesis. The microarray technology used was Agilent Human miRNA microarrays (V2, G4470B, Agilent Technologies, Sanger Database version 10.1). Data normalization was performed with Feature Extraction Software (Agilent Technologies). Extracted intensities were background corrected using the normexp method with an offset of 50 . Background corrected log2-transformed intensities were normalized using invariants normalization to make data from all arrays comparable . Microarray probes were collapsed to miRNAs by taking the median intensity of the respective probes per miRNA. For determining differentially regulated miRNAs moderated t-tests were applied using Limma (Linear Models for Microarray Data) . MiRNAs with FDR adjusted p-value <0.05 and additionally a fold change exceeding 1.2 in absolute value were selected as the relevant ones. All statistical analyses were performed with the Bioconductor project in the R statistical environment .
DNA methylation analysis by pyrosequencing
Bisulphite conversion of 500 ng of each DNA sample was performed with EZ DNA Methylation-Gold Kit (Zymo Research, Milan Italy) according to the manufacturer’s recommendations. PCR for miR-204 promoters was performed with 1 μl of bisulphite converted DNA under standard conditions with biotinylated primers using an annealing temperature of 60 °C. Primer sequences are: miR-204 Forward (5′-[Btn]TGGTTTTTTTTTAAATTAAGTTAGTAAAGT-3′); miR-204 Reverse (5′- ACAACCTACACAAAACAACCTATAATC-3′); miR-204 Sequencing (5′-ACAACCTACTCCAAAAT-3′) and were designed with PyroMark Assay Design 2.0. PCR products were observed on 2 % agarose gels before pyrosequencing analysis. Reactions were performed in a PyroMark Q96 System version 2.0.6 (Qiagen, Milan, Italy) and the methylation values of the CpG dinucleotides were obtained using Pyro Q-CpG 1.0.9 (Qiagen, Milan, Italy).
Cell culture and expression arrays
HaCaT cells were cultured in DMEM medium (Sigma) supplemented with 10 % fetal bovine serum in standard conditions. SCC13 and SCC12 cells were cultured in DMEM:F12 medium (Sigma) supplemented with 10 % fetal bovine serum, penicillin–streptomycin (100 mg/ml, Invitrogen), adenine (180 mM, Sigma), human recombinant EGF (10 ng/ml Peprotech), insulin (5 mg/ml, Sigma), cholera toxin (0.1 nM, Gentaur) and hydrocortisone (0.5 mg/ml, Sigma). Human epidermal keratinocytes (HEK, Life Technologies) were cultured in KSFM medium (Thermo Fisher). To restore the expression of DNA-methylated miR-204, HaCaT cells were treated with the DNA demethylating agent 5-aza-2′-deoxycytidine (Sigma) at 10 μM for 120 h. To generate retroviral stocks for stable expression of miR-204 antagomir, 293T cells were transfected with Fugene (Roche) with miRZIP-204 anti-miR-204 construct (shRNA targeting miR-204 microRNA, cloned in pGreenPuro HIV based lentivector, BioCat), or with the corresponding pGreenPuro scrambled hairpin negative control (BioCat). HaCaT cells were transduced with lentiviral supernatants in the presence of Polybrene (4 μg/ml; Sigma). Two days after the infection, control and miR-204 antagomir expressing HaCaT cells were FACS-sorted. Total RNA extraction from GFP-positive cells was performed with Genelute mammalian total RNA kit (Sigma). Samples were then processed following the protocols provided by Ambion (Ambion WT Expression, Protocol Number 4425209; Revised B 05/2009) and Affymetrix (Affymetrix WT Terminal Labeling and Hybridization User Manual, Protocol Number 702808 and Expression, Wash, Stain and Scan User Manual Protocol Number 702731). The expression platform used was Human Gene 2.0 ST (Affymetrix®).
Transient expression of miRIDIAN microRNA human hsa-miR-204-5p-hairpin inhibitor (Dharmacon), pre-miR™ miRNA-204-5p precursor (Life Technologies) or the corresponding controls (miRIDIAN microRNA hairpin inhibitor negative control #1, Dharmacon, and pre-miR™ miRNA precursor negative control #1, Life Technologies, respectively) at 50 nM was performed with DHARMAFECT 4 as described earlier . Twenty four hours later, cells were treated with 10 ng/ml EGF, IL6 or FGF9 (Peprotech) for the times indicated, in the absence or the presence of 50 μM NSC87877 (Sigma).
RNA extraction and quantitative reverse transcription PCR (qRT-PCR)
RNA was extracted from frozen biopsies with the mirVana™ miRNA Isolation Kit (Applied Biosystems). qRT-PCR for miR-204 transcripts was performed with TaqMan® microRNA assays (Applied Biosystems), using RNU48 levels to normalize the data. Real-time PCR reactions were performed in triplicate in 384-well plates and were run in a QuantStudio™ 12K Flex Real-Time PCR System. The 2-(ΔΔCt) method  was used to determine relative quantitative levels of miR-204. qRT-PCR to analyze gene expression was performed using SYBR Green PCR master mix (Applied Biosystems, Life Technologies) and GAPDH expression was used as endogenous control for normalization purposes. Primers used for qRT-PCR were: PTPN11 (5′- AAAGGGGAGAGCAATGACGG-3′; 5′- GGGGCTGCTTGAGTTGTAGT-3′); MFAP3 (5′-CCCGCTGTGCTTTTGTTGTA-3′; 5′-ATCCGAGTGGGAACTTGCTG-3′); FGF9 (5′-GCAGTCACGGACTTGGATCA-3′; 5′-TCCAGAATGCCAAATCGGCT-3′); FGF2 (5′- CCCCAGAAAACCCGAGCGA-3′; 5′- CGTCCGCTAATCTGGCACC-3′); CBL (5′- ATTCTCCATGGCCCCACAAG-3′; 5′- GAATATGGCCGGTCTGGAGG-3′); MAP3K2 (5′- GACCAACTGGAGATTGGGCA-3′; 5′- GGTCTCAGGACTATCGGGGT-3′); C1QBP (5′- TTTGCATCTGCACGTGTTCG-3′; 5′- CAGGAGCTGCCGGAAAGG-3′); PDP2 (5′- CTTGCAAACACGCCTACTGT-3′; 5′- GCACAAAAAGCTGCACCAGT-3′); IL6 (5′-AGCCACTCACCTCTTCAGAAC-3′; 5′-GTGCCTCTTTGCTGCTTTCAG-3′); GAPDH (5′-AGTCAGCCGCATCTTCTTTTG-3′; 5′-AAATCCGTTGACTCCGACCTT-3′).
Protein extraction and Western blot
These protocols were performed following standard techniques. STAT3 (C-20), ERK2 (C-14) and PTPN11 (SH-PTP2 B1) antibodies were purchased from Santa Cruz Biotech. Activated MAP kinase (diphosphorylated ERK1/2) antibody was from Sigma. Lamin B1, Actin, phospho-Y705-STAT3 and phospho-S727-STAT3 antibodies were from Abcam. HRP-conjugated secondary antibodies were purchased from Dako. Western blot quantification was performed using Image J densitometry software. Intensity of individual bands was normalized to Actin signal and referred to control conditions.
Transfection and luciferase assay for PTPN11 3′ UTR construct
Oligonucleotides corresponding to the PTPN11 3′UTR region containing the miR-204 binding site (5′-TCGAGTTGCAATTGTCAAGTGTTTTGTTGTAGCTTAGTATCCATAAGGGAAAC TTAGACTATAGACATAGATCTGC-3′; 5′-GGCCGCAGATCTATGTCTATAGTCTAAGTT TCCCTTATGGATACTAAGCTACAACAAAACACTTGACAATTGCAAC-3′) were annealed and inserted in the psiCHECK2 plasmid (Promega), immediately downstream from the stop codon of Renilla luciferase, and the sequence was confirmed by sequencing. For the luciferase assay, the transfection was performed using Fugene HD (Roche). Once ready, HaCaT cells were transfected with the reporter plasmid along with miR-204 antagomir, miR-204 mimic or with their corresponding hairpin inhibitor or precursor negative controls at 50 nM. Renilla and Firefly luciferase activities were measured consecutively using a dual-luciferase assay (Promega) 24 h after transfection. Renilla luciferase activity was normalized to Firefly luciferase activity expressed from the psiCHECK2 plasmid.
Immunofluorescence and immunohistochemistry
Immunofluorescence was performed on 4 μm-thick tissue sections of 14 non-peritumoral AKs, 11 cSCC and their corresponding adjacent AK samples. Heat-induced antigen retrieval was performed in 10 mM citrate (pH 6) for 15 min in a pressure cooker. The slides were then incubated with STAT3 (phospho-Y705) antibody (Abcam) for 12 h. After washing, the Envision + antibody reagent was applied (Dako). Reactions were developed with the TSA™ Plus Cyanine3/Fluorescein System (PerkinElmer). Nuclei were counterstained with Hoesch and examined with a fluorescence microscope (Olympus BX61). Cases where pY705-STAT3 signal was detected in more than 20 % of cells were considered as positive. PTPN11 immunohistochemistry was performed by antigenic retrieval in 10 mM citrate (pH 6) for 15 min in a pressure cooker. After primary antibody incubation, the Envision + System-HRP antibody reagent was applied (Dako) and reactions were developed using diaminobenzidine (DAB, Dako).
Identification of miRNAs differentially expressed in actinic keratosis (AK) and cutaneous squamous cell carcinoma (cSCC)
MiR-204 expression is frequently silenced by DNA methylation in cSCC
Identification of miR-204 targets in non-transformed keratinocytes
PTPN11 (also called SHP2) is a member of the protein tyrosine phosphatase (PTP) family that directly affects STAT activity in A431 cell nuclei , while FAM117B is a protein-coding gene of unknown function expressed in the nucleus of basal keratinocytes . Because the FGF-STAT3 pathway was identified to be deregulated in miR-204-depleted HaCaT we hypothesized that the PTPN11 phosphatase could be directly targeted by miR-204 to modulate STAT3 signaling. To investigate this possibility, HaCaT cells were transiently transfected with miR-204 antagomir (hairpin inhibitor), miR-204 mimic (pre-miRNA-204 precursor) or with the corresponding antagomir or precursor negative controls (Fig. 3c). According to the gene expression data, miR-204 antagomir transfected cells displayed lower PTPN11 mRNA levels than their control antagomir counterparts, whereas the expression of the miR-204 mimic enhanced PTPN11 levels (Additional file 1: Figure S3A). MiRNAs can silence cytoplasmic mRNA translation either by triggering an endonuclease cleavage, by promoting translation repression or by accelerating mRNA decapping . Since miRNAs can interfere with mRNA translation in a manner independent of changes at the mRNA level  we tested miR-204 impact on PTPN11 protein expression. Contrary to the effect observed on PTPN11 mRNA levels, transfection of miR-204 antagomir in HaCaT cells up-regulated PTPN11 protein levels, whereas miR-204 mimic resulted in downregulation of PTPN11 expression (Fig. 3d). Next, we performed a reporter assay in which the Renilla luciferase reporter gene was under control of the PTPN11 3′ UTR. We observed that sequestration of endogenous miR-204 with the specific antagomir increased the expression of the Renilla reporter, while the expression of the miR-204 mimic reduced the expression of the reporter, confirming that PTPN11 may be a primary target of miR-204 through its 3′ UTR (Fig. 3e).
Finally we validated the microarray data by qRT-PCR analysis of some differentially expressed genes in HaCaT cells transfected with miR-204 antagomir or miR-204 mimic (Fig. 3f). The fact that miR-204 inhibition represses while miR-204 induces the expression of the validated genes suggests that miR-204 might indirectly regulate their expression. We also tested some cellular functions identified by IPA as modulated by miR-204 expression. Cell cycle status assessed by BrdU staining to detect DNA synthesis revealed that miR-204 depletion in HaCaT cells impairs BrdU incorporation while cells transfected with miR-204 mimic displayed an enhancement in BrdU staining (Additional file 1: Figure S3B). In contrast, cell migration enhancement was observed in HaCaT cells transfected with miR-204 mimic, whereas no effect was observed in those transfected with the antagomir (Additional file 1: Figure S3C). Therefore gene expression and migration modulation by miR-204 were more effectively observed when this miRNA was overexpressed, whereas miR-204 depletion mostly resulted in mild gene expression changes and unaffected cell motility. This could be explained by the relatively low endogenous miR-204 levels due to partial DNA methylation of the miR-204 locus in HaCaT cells (see above), despite being a non-tumorigenic cell line. Since interfering endogenous miR-204 levels in HaCaT cells would likely result in subtle gene expression changes, from then on the analysis of miR-204-mediated changes was carried out by miR-204 overexpression.
MiR-204 affects the expression of EGF-modulated genes and regulates STAT3 signaling
Next we tested whether miR-204 overexpression impinged on EGF-induced gene expression. To this end, HaCaT cells transfected with control or miR-204 mimics were stimulated with EGF for 30 min or 24 h. MiR-204 mimic-transfected cells displayed higher basal and induced mRNA levels of the tested genes at short time points. However, mRNA levels of these genes were further increased in control cells at 24 h, whereas the levels of these mRNAs dropped in miR-204 mimic cells (Fig. 4b). These results indicated that miR-204 expression modulates EGF signaling intensity and timing, with an initial activation followed by a downregulation at longer time periods.
The miR-204 target PTPN11 can modulate STAT3 [53–55]. Accordingly, the kinetic of PTPN11 expression after the stimulation with the cytokine closely resembled that of total STAT3 and, as observed for STAT3, enhanced PTPN11 expression was stronger in miR-204 mimic cells than control cells at these time points (Fig. 5a).
STAT3 can be activated in a biphasic pattern, with a first wave of activation within 1 h of cytokine treatment and a sustained wave in which STAT3 remains active for many hours . One of the cytokines able to induce sustained STAT3 activation is IL6 . We analysed MAPK and STAT3 pathways after 24 h of treatment with EGF, FGF9 and IL6. At this time point pERK1/2 levels were higher in miR-204 mimic cells stimulated with EGF, FGF9 and IL6 (Fig. 5b). In contrast, control cells displayed enhanced STAT3 expression by FGF9 and IL6 and sustained Y705-STAT3 phosphorylation by IL6 that were not observed in miR-204 overexpressing cells. It has been reported that Y705-STAT3 phosphorylation can be negatively modulated by ERK induced phosphorylation of S727-STAT3 . However, similar levels of S727 phosphorylation could be observed regardless miR-204 expression. PTPN11 levels were higher in untreated cells as well as upon FGF9 and IL6 treatments in control cells when compared to similar conditions in miR-204 overexpressing cells at this time point. All these data indicated that miR-204 induced STAT3-mediated signaling at short time points but resulted in impairment in STAT3 expression and activation triggered by different cytokines signaling through this pathway at longer time points.
PTPN11 has been previously reported to regulate the STAT3 signaling pathway in different cellular contexts. To directly address its role in our experimental conditions, HaCaT cells were stimulated for 24 h with EGF, FGF9 and IL6 in the presence of NSC87877, that specifically inhibits PTPN11 , and activated and total STAT3 levels were analysed by western blot. Control cells displayed sustained Y705-STAT3 phosphorylation and enhanced STAT3 expression by IL6 that were not observed in NSC87877 treated cells (Fig. 5c). This result parallels that of the effect of IL6 in miR-204 overexpressing HaCaT cells, supporting the notion that PTPN11 mediates, at least in part, some of the miR-204 regulatory effects on STAT3 signaling pathway.
Differential subcellular localization of pY705-STAT3 in non-peritumoral and cSCC-adjacent AKs
MiRNAs alterations are involved in the development of several skin cancers . We have identified miR-204, previously found to be down-regulated in head and neck squamous cell carcinoma (HNSCC) [62, 63] and cSCC , as the only miRNA to be significantly downregulated in cSCC when compared to AK. MiR-204 is located at 9q21.1-q22-3, a cancer-associated genomic region exhibiting a high frequency of loss of heterozigosity (LOH) in HNSCC, but rarely lost in cSCC. DNA methylation represents the most studied epigenetic regulatory mechanism of modifying gene expression in both coding and non-coding genes such as microRNAs [64–67]. We here report hypermethylation of the miR-204 promoter in AKs and in a high proportion of cSCCs. Since this genomic region also displays STAT3-binding sites, and phosphorylated STAT3 suppresses miR-204 expression in pulmonary arterial hypertension , STAT3 activation by alterations in proteins modulating this pathway could account as an alternative mechanism for miR-204 dowregulation in cSCC . Therefore, different DNA-methylation dependent and independent mechanisms could mediate miR-204 repression in cSCC.
Several studies in epithelial cell lines have shown that members of the miR-200 family are repressed during epithelial-to-mesenchymal transition (EMT) . Some in vitro studies with hepatic, renal, ovarian and breast cancer cell lines have shown that loss of miR-204 results in activation of an EMT, leading to cancer cell migration and invasion [71, 72]. We have previously described that EMT is common in metastatic cSCC . Although we have not detected a relationship between miR-204 and EMT, the demonstration of miR-204 silencing in cSCC samples regardless of their biological behaviour makes unlikely the possibility that miR-204 plays a relevant role in this process.
All miRNAs have complex multifunctional and cell context-dependent roles as either repressors or enhancers . Accordingly, we observed a positive short-term induction of the STAT3 pathway by miR-204, whereas in the long-term miR-204 downregulates STAT3 pathway, while favouring MAPK pathway activation. Conversely, low levels of miR-204 in cSCC would favour STAT3 pathway activation, which is also mediated by IL6 . In pulmonary hypertension an auto-regulatory positive feedback has been reported, in which STAT3 activation might perpetuate miR-204 downregulation .
The MAPK/ERK pathway is the most important cell survival pathway in non-tumorigenic keratinocytes, and is triggered by EGF. Forced expression of EGFR in immortalized keratinocytes (HaCaT cells) has been associated with enhanced EGFR activation but not with Y705-STAT3 phosphorylation . In contrast, malignant progression of the keratinocyte towards SCC is associated with negative MAPK regulation and EGFR-induced STAT3 activation [50–52], which drives expression of multiple antiapoptotic molecules [50, 52, 75–78]. Accordingly, STAT3 expression was found to be consistently increased in cSCC [59, 60, 79, 80]. Immunofluorescence detection of pY705-STAT3 signal suggests an inverse correlation between the pY705-STAT3 subcelular localization and miR-204 expression in AK and cSCC. Non-peritumoral distant AKs showed phosphorylated STAT3 excluded from the nuclei, whereas AKs adjacent to cSCC displayed nuclear pY705-STAT3 signal when their corresponding tumors were also positive for this marker. Nuclear p705-STAT3 signal in adjacent AKs and cSCC [59, 60] (our own data) suggests that miR-204 already suffers a downregulation in some areas of the AK that will progress to cSCC. On the contrary, pY705-STAT3 signal at the plasma membrane/cytoplasm in some of the non-peritumoral AKs suggests that certain degree of STAT3 activation is already present in AK, but miR-204 expression would inhibit the expression of proteins regulating effective pY705-STAT3 translocation to the nucleus, where these proteins act as transcription factors.
MicroRNA-204 may act as a “rheostat” that controls the signaling towards the MAPK pathway or the STAT3 pathway. MiR-204 would favour MAPK pathway and make the AK cell insensitive to the induction of the STAT3 pathway despite the presence of procarcinogenic inflammatory mediators. In most cSCCs, miR-204 promoter DNA methylation accounts for miR-204 silencing, which would result in STAT3 activation and translocation to the nucleus, a common feature found in 80 % of cSCCs.
cSCC, cutaneous squamous cell carcinoma; AK, actinic keratosis; KIN, keratinocyte intraepidermal neoplasm; CSE skin, chronically sun-exposed skin; HNSCC, head and neck squamous cell carcinoma; EMT, epithelial-to-mesenchymal transition; microRNA, miRNA; IPA, Ingenuity Pathway Analysis
The authors would like to thank than Sergi Mojal, Lara Nonell, Sarah Bonnin, Anna Ferrer and Tania Lobato for their technical support.
This study was supported by grants PI15/00236, PI041728, PI10/00785, FIS PI11/02070, RD09/0076/00036, RD06/0020/0040, and PS09/00973 from Fondo de Investigación Sanitaria (FIS), Instituto de Salud Carlos III FEDER, SAF2010-16089, Ministerio de Economía y Competitividad, Spain, from the “Xarxa de Bancs de tumors sponsored by Pla Director d’Oncologia de Catalunya (XBTC)”. IHM is an investigator at the Miguel Servet program (Instituto de Salud Carlos III). ADL is funded by the ISCIII through a research contract Río Hortega (CM14/00067). ME is an ICREA Research Professor.
Availability of data and materials
Additional data are available as Supplementary information.
AT, BE, AD-L, EH-R, JS, ME, RMP and IH-M conceived the experiments and analysed data. RS, AD-L, EA, IH-M carried out experiments. All authors were involved in writing the paper and had final approval of the submitted and published versions.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Approval to conduct this study was obtained from the Ethics Committees from the Institution (Number 2010/3913/I), in accordance with the guidelines of the Helsinki Declaration of 1975, revised in 1983.
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.
- Stern RS, Lunder EJ. Risk of squamous cell carcinoma and methoxsalen (psoralen) and UV-A radiation (PUVA). A meta-analysis. Arch Dermatol. 1998;134(12):1582–5.View ArticlePubMedGoogle Scholar
- Mittelbronn MA, Mullins DL, Ramos-Caro FA, Flowers FP. Frequency of pre-existing actinic keratosis in cutaneous squamous cell carcinoma. Int J Dermatol. 1998;37(9):677–81.View ArticlePubMedGoogle Scholar
- Glogau RG. The risk of progression to invasive disease. J Am Acad Dermatol. 2000;42(1 Pt 2):23–4.View ArticlePubMedGoogle Scholar
- Ratushny V, Gober MD, Hick R, Ridky TW, Seykora JT. From keratinocyte to cancer: the pathogenesis and modeling of cutaneous squamous cell carcinoma. J Clin Invest. 2012;122(2):464–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Brantsch KD, Meisner C, Schonfisch B, Trilling B, Wehner-Caroli J, Rocken M, et al. Analysis of risk factors determining prognosis of cutaneous squamous-cell carcinoma: a prospective study. Lancet Oncol. 2008;9(8):713–20.View ArticlePubMedGoogle Scholar
- Wells III JL, Shirai K. Systemic therapy for squamous cell carcinoma of the skin in organ transplant recipients. Am J Clin Oncol. 2012;35(5):498–503.View ArticlePubMedGoogle Scholar
- Lambert SR, Mladkova N, Gulati A, Hamoudi R, Purdie K, Cerio R, et al. Key differences identified between actinic keratosis and cutaneous squamous cell carcinoma by transcriptome profiling. Br J Cancer. 2014;110(2):520–9.View ArticlePubMedGoogle Scholar
- Einspahr JG, Calvert V, Alberts DS, Curiel-Lewandrowski C, Warneke J, Krouse R, et al. Functional protein pathway activation mapping of the progression of normal skin to squamous cell carcinoma. Cancer Prev Res (Phila). 2012;5(3):403–13.View ArticleGoogle Scholar
- Ruini C, Witkowski AM, Cesinaro A, De Teixeira CN, Pellacani G. From actinic keratosis to squamous cell carcinoma: evidence of morphologic and biologic progression. J Am Acad Dermatol. 2015;72(1 Suppl):S8–S10.View ArticlePubMedGoogle Scholar
- Chisholm C, Greene Jr JF. Progression from atypical/dysplastic intraepidermal proliferations and carcinoma in situ to invasive tumors: a pathway based on current knowledge. Am J Dermatopathol. 2011;33(8):803–10.View ArticlePubMedGoogle Scholar
- Campbell C, Quinn AG, Ro YS, Angus B, Rees JL. p53 mutations are common and early events that precede tumor invasion in squamous cell neoplasia of the skin. J Invest Dermatol. 1993;100(6):746–8.View ArticlePubMedGoogle Scholar
- Wang NJ, Sanborn Z, Arnett KL, Bayston LJ, Liao W, Proby CM, et al. Loss-of-function mutations in Notch receptors in cutaneous and lung squamous cell carcinoma. Proc Natl Acad Sci U S A. 2011;108(43):17761–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Rehman I, Takata M, Wu YY, Rees JL. Genetic change in actinic keratoses. Oncogene. 1996;12(12):2483–90.PubMedGoogle Scholar
- Toll A, Salgado R, Yebenes M, Martin-Ezquerra G, Gilaberte M, Baro T, et al. MYC gene numerical aberrations in actinic keratosis and cutaneous squamous cell carcinoma. Br J Dermatol. 2009;161(5):1112–8.View ArticlePubMedGoogle Scholar
- Toll A, Salgado R, Yebenes M, Martin-Ezquerra G, Gilaberte M, Baro T, et al. Epidermal growth factor receptor gene numerical aberrations are frequent events in actinic keratoses and invasive cutaneous squamous cell carcinomas. Exp Dermatol. 2010;19(2):151–3.View ArticlePubMedGoogle Scholar
- Londin E, Loher P, Telonis AG, Quann K, Clark P, Jing Y, et al. Analysis of 13 cell types reveals evidence for the expression of numerous novel primate- and tissue-specific microRNAs. Proc Natl Acad Sci U S A. 2015;112(10):E1106–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH, Cuppen E. Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005;120(1):21–4.View ArticlePubMedGoogle Scholar
- Maute RL, Dalla-Favera R, Basso K. RNAs with multiple personalities. Wiley Interdiscip Rev RNA. 2014;5(1):1–13.View ArticlePubMedGoogle Scholar
- Gastaldi C, Bertero T, Xu N, Bourget-Ponzio I, Lebrigand K, Fourre S, et al. miR-193b/365a cluster controls progression of epidermal squamous cell carcinoma. Carcinogenesis. 2014;35(5):1110–20.View ArticlePubMedGoogle Scholar
- Kanitz A, Imig J, Dziunycz PJ, Primorac A, Galgano A, Hofbauer GF, et al. The expression levels of microRNA-361-5p and its target VEGFA are inversely correlated in human cutaneous squamous cell carcinoma. PLoS One. 2012;7(11), e49568.View ArticlePubMedPubMed CentralGoogle Scholar
- Lefort K, Brooks Y, Ostano P, Cario-Andre M, Calpini V, Guinea-Viniegra J, et al. A miR-34a-SIRT6 axis in the squamous cell differentiation network. EMBO J. 2013;32(16):2248–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang SH, Zhou JD, He QY, Yin ZQ, Cao K, Luo CQ. MiR-199a inhibits the ability of proliferation and migration by regulating CD44-Ezrin signaling in cutaneous squamous cell carcinoma cells. Int J Clin Exp Pathol. 2014;7(10):7131–41.PubMedPubMed CentralGoogle Scholar
- Xu N, Zhang L, Meisgen F, Harada M, Heilborn J, Homey B, et al. MicroRNA-125b down-regulates matrix metallopeptidase 13 and inhibits cutaneous squamous cell carcinoma cell proliferation, migration, and invasion. J Biol Chem. 2012;287(35):29899–908.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamane K, Jinnin M, Etoh T, Kobayashi Y, Shimozono N, Fukushima S, et al. Down-regulation of miR-124/-214 in cutaneous squamous cell carcinoma mediates abnormal cell proliferation via the induction of ERK. J Mol Med (Berl). 2013;91(1):69–81.View ArticleGoogle Scholar
- Sand M, Skrygan M, Georgas D, Sand D, Hahn SA, Gambichler T, et al. Microarray analysis of microRNA expression in cutaneous squamous cell carcinoma. J Dermatol Sci. 2012;68(3):119–26.View ArticlePubMedGoogle Scholar
- Bhandari A, Gordon W, Dizon D, Hopkin AS, Gordon E, Yu Z, et al. The Grainyhead transcription factor Grhl3/Get1 suppresses miR-21 expression and tumorigenesis in skin: modulation of the miR-21 target MSH2 by RNA-binding protein DND1. Oncogene. 2013;32(12):1497–507.View ArticlePubMedGoogle Scholar
- Bruegger C, Kempf W, Spoerri I, Arnold AW, Itin PH, Burger B. MicroRNA expression differs in cutaneous squamous cell carcinomas and healthy skin of immunocompetent individuals. Exp Dermatol. 2013;22(6):426–8.View ArticlePubMedGoogle Scholar
- Dziunycz P, Iotzova-Weiss G, Eloranta JJ, Lauchli S, Hafner J, French LE, et al. Squamous cell carcinoma of the skin shows a distinct microRNA profile modulated by UV radiation. J Invest Dermatol. 2010;130(11):2686–9.View ArticlePubMedGoogle Scholar
- Olasz EB, Seline LN, Schock AM, Duncan NE, Lopez A, Lazar J, et al. MicroRNA-135b regulates leucine zipper tumor suppressor 1 in cutaneous squamous cell carcinoma. PLoS One. 2015;10(5), e0125412.View ArticlePubMedPubMed CentralGoogle Scholar
- Syed DN, Lall RK, Mukhtar H. MicroRNAs and photocarcinogenesis. Photochem Photobiol. 2015;91(1):173–87.View ArticlePubMedGoogle Scholar
- Wang A, Landen NX, Meisgen F, Lohcharoenkal W, Stahle M, Sonkoly E, et al. MicroRNA-31 is overexpressed in cutaneous squamous cell carcinoma and regulates cell motility and colony formation ability of tumor cells. PLoS One. 2014;9(7), e103206.View ArticlePubMedPubMed CentralGoogle Scholar
- White RA, Neiman JM, Reddi A, Han G, Birlea S, Mitra D, et al. Epithelial stem cell mutations that promote squamous cell carcinoma metastasis. J Clin Invest. 2013;123(10):4390–404.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou M, Liu W, Ma S, Cao H, Peng X, Guo L, et al. A novel onco-miR-365 induces cutaneous squamous cell carcinoma. Carcinogenesis. 2013;34(7):1653–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou M, Zhou L, Zheng L, Guo L, Wang Y, Liu H, et al. miR-365 promotes cutaneous squamous cell carcinoma (CSCC) through targeting nuclear factor I/B (NFIB). PLoS One. 2014;9(6):e100620.View ArticlePubMedPubMed CentralGoogle Scholar
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, et al. A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007;23(20):2700–7.View ArticlePubMedGoogle Scholar
- Pradervand S, Weber J, Thomas J, Bueno M, Wirapati P, Lefort K, et al. Impact of normalization on miRNA microarray expression profiling. RNA. 2009;15(3):493–501.View ArticlePubMedPubMed CentralGoogle Scholar
- Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.PubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.View ArticlePubMedPubMed CentralGoogle Scholar
- Martinez-Romero C, Rooman I, Skoudy A, Guerra C, Molero X, Gonzalez A, et al. The epigenetic regulators Bmi1 and Ring1B are differentially regulated in pancreatitis and pancreatic ductal adenocarcinoma. J Pathol. 2009;219(2):205–13.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- Wang FE, Zhang C, Maminishkis A, Dong L, Zhi C, Li R, et al. MicroRNA-204/211 alters epithelial physiology. FASEB J. 2010;24(5):1552–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Ying Z, Li Y, Wu J, Zhu X, Yang Y, Tian H, et al. Loss of miR-204 expression enhances glioma migration and stem cell-like phenotype. Cancer Res. 2013;73(2):990–9.View ArticlePubMedGoogle Scholar
- Jansson MD, Lund AH. MicroRNA and cancer. Mol Oncol. 2012;6(6):590–610.View ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120(1):15–20.View ArticlePubMedGoogle Scholar
- Wu TR, Hong YK, Wang XD, Ling MY, Dragoi AM, Chung AS, et al. SHP-2 is a dual-specificity phosphatase involved in Stat1 dephosphorylation at both tyrosine and serine residues in nuclei. J Biol Chem. 2002;277(49):47572–80.View ArticlePubMedGoogle Scholar
- Asplund A, Gry BM, Sundquist C, Stromberg S, Edlund K, Ostman A, et al. Expression profiling of microdissected cell populations selected from basal cells in normal epidermis and basal cell carcinoma. Br J Dermatol. 2008;158(3):527–38.View ArticlePubMedGoogle Scholar
- Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006;20(5):515–24.View ArticlePubMedGoogle Scholar
- Jost M, Huggett TM, Kari C, Boise LH, Rodeck U. Epidermal growth factor receptor-dependent control of keratinocyte survival and Bcl-xL expression through a MEK-dependent pathway. J Biol Chem. 2001;276(9):6320–6.View ArticlePubMedGoogle Scholar
- Raj D, Brash DE, Grossman D. Keratinocyte apoptosis in epidermal development and disease. J Invest Dermatol. 2006;126(2):243–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Bowman T, Garcia R, Turkson J, Jove R. STATs in oncogenesis. Oncogene. 2000;19(21):2474–88.View ArticlePubMedGoogle Scholar
- Quadros MR, Peruzzi F, Kari C, Rodeck U. Complex regulation of signal transducers and activators of transcription 3 activation in normal and malignant keratinocytes. Cancer Res. 2004;64(11):3934–9.View ArticlePubMedGoogle Scholar
- Rubin GJ, Tweardy DJ, Melhem MF. Asynchronous modulation of transforming growth factor alpha and epidermal growth factor receptor protein expression in progression of premalignant lesions to head and neck squamous cell carcinoma. Clin Cancer Res. 1998;4(1):13–20.Google Scholar
- Blechacz BR, Smoot RL, Bronk SF, Werneburg NW, Sirica AE, Gores GJ. Sorafenib inhibits signal transducer and activator of transcription-3 signaling in cholangiocarcinoma cells by activating the phosphatase shatterproof 2. Hepatology. 2009;50(6):1861–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Usenko T, Chan G, Torlakovic E, Klingmuller U, Neel BG. Leukemogenic Ptpn11 allele causes defective erythropoiesis in mice. PLoS One. 2014;9(10), e109682.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu JH, Goswami R, Cai X, Exum ST, Huang X, Zhang L, et al. Regulation of the platelet-derived growth factor receptor-beta by G protein-coupled receptor kinase-5 in vascular smooth muscle cells involves the phosphatase Shp2. J Biol Chem. 2006;281(49):37758–72.View ArticlePubMedGoogle Scholar
- Wang Y, van Boxel-Dezaire AH, Cheon H, Yang J, Stark GR. STAT3 activation in response to IL-6 is prolonged by the binding of IL-6 receptor to EGF receptor. Proc Natl Acad Sci U S A. 2013;110(42):16975–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung J, Uchida E, Grammer TC, Blenis J. STAT3 serine phosphorylation by ERK-dependent and -independent pathways negatively modulates its tyrosine phosphorylation. Mol Cell Biol. 1997;17(11):6508–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Sung SS, Yip ML, Lawrence HR, Ren Y, Guida WC, et al. Discovery of a novel shp2 protein tyrosine phosphatase inhibitor. Mol Pharmacol. 2006;70(2):562–70.View ArticlePubMedGoogle Scholar
- Bito T, Sumita N, Ashida M, Budiyanto A, Ueda M, Ichihashi M, et al. Inhibition of epidermal growth factor receptor and PI3K/Akt signaling suppresses cell proliferation and survival through regulation of Stat3 activation in human cutaneous squamous cell carcinoma. J Skin Cancer. 2011;2011:874571.View ArticlePubMedGoogle Scholar
- Suiqing C, Min Z, Lirong C. Overexpression of phosphorylated-STAT3 correlated with the invasion and metastasis of cutaneous squamous cell carcinoma. J Dermatol. 2005;32(5):354–60.View ArticlePubMedGoogle Scholar
- Sand M, Sand D, Altmeyer P, Bechara FG. MicroRNA in non-melanoma skin cancer. Cancer Biomark. 2012;11(6):253–7.PubMedGoogle Scholar
- Lee Y, Yang X, Huang Y, Fan H, Zhang Q, Wu Y, et al. Network modeling identifies molecular functions targeted by miR-204 to suppress head and neck tumor metastasis. PLoS Comput Biol. 2010;6(4), e1000730.View ArticlePubMedPubMed CentralGoogle Scholar
- Nohata N, Hanazawa T, Kinoshita T, Okamoto Y, Seki N. MicroRNAs function as tumor suppressors or oncogenes: aberrant expression of microRNAs in head and neck squamous cell carcinoma. Auris Nasus Larynx. 2013;40(2):143–9.View ArticlePubMedGoogle Scholar
- Leveille N, Melo CA, Rooijers K, Diaz-Lagares A, Melo SA, Korkmaz G, et al. Genome-wide profiling of p53-regulated enhancer RNAs uncovers a subset of enhancers controlled by a lncRNA. Nat Commun. 2015;6:6520.View ArticlePubMedPubMed CentralGoogle Scholar
- Quintavalle C, Mangani D, Roscigno G, Romano G, Diaz-Lagares A, Iaboni M, et al. MiR-221/222 target the DNA methyltransferase MGMT in glioma cells. PLoS One. 2013;8(9), e74466.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandoval J, Diaz-Lagares A, Salgado R, Servitje O, Climent F, Ortiz-Romero PL, et al. MicroRNA expression profiling and DNA methylation signature for deregulated microRNA in cutaneous T-cell lymphoma. J Invest Dermatol. 2015;135(4):1128–37.View ArticlePubMedGoogle Scholar
- Simo-Riudalbas L, Diaz-Lagares A, Gatto S, Gagliardi M, Crujeiras AB, Matarazzo MR, et al. Genome-wide DNA methylation analysis identifies novel hypomethylated non-pericentromeric genes with potential clinical implications in ICF syndrome. PLoS One. 2015;10(7), e0132517.View ArticlePubMedPubMed CentralGoogle Scholar
- Courboulin A, Paulin R, Giguere NJ, Saksouk N, Perreault T, Meloche J, et al. Role for miR-204 in human pulmonary arterial hypertension. J Exp Med. 2011;208(3):535–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishio E, Miura Y, Kawaguchi M, Morita A. Nuclear translocation of ATBF1 is a potential prognostic marker for skin cancer. Acta Dermatovenerol Croat. 2012;20(4):239–45.PubMedGoogle Scholar
- Korpal M, Kang Y. The emerging role of miR-200 family of microRNAs in epithelial-mesenchymal transition and cancer metastasis. RNA Biol. 2008;5(3):115–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Imam JS, Plyler JR, Bansal H, Prajapati S, Bansal S, Rebeles J, et al. Genomic loss of tumor suppressor miRNA-204 promotes cancer cell migration and invasion by activating AKT/mTOR/Rac1 signaling and actin reorganization. PLoS One. 2012;7(12), e52397.View ArticlePubMedPubMed CentralGoogle Scholar
- Qiu YH, Wei YP, Shen NJ, Wang ZC, Kan T, Yu WL, et al. miR-204 inhibits epithelial to mesenchymal transition by targeting slug in intrahepatic cholangiocarcinoma cells. Cell Physiol Biochem. 2013;32(5):1331–41.View ArticlePubMedGoogle Scholar
- Toll A, Masferrer E, Hernandez-Ruiz ME, Ferrandiz-Pulido C, Yebenes M, Jaka A, et al. Epithelial to mesenchymal transition markers are associated with an increased metastatic risk in primary cutaneous squamous cell carcinomas but are attenuated in lymph node metastases. J Dermatol Sci. 2013;72(2):93–102.View ArticlePubMedGoogle Scholar
- Bertero T, Lu Y, Annis S, Hale A, Bhat B, Saggar R, et al. Systems-level regulation of microRNA networks by miR-130/301 promotes pulmonary hypertension. J Clin Invest. 2014;124(8):3514–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Buettner R, Mora LB, Jove R. Activated STAT signaling in human tumors provides novel molecular targets for therapeutic intervention. Clin Cancer Res. 2002;8(4):945–54.PubMedGoogle Scholar
- Chan KS, Carbajal S, Kiguchi K, Clifford J, Sano S, DiGiovanni J. Epidermal growth factor receptor-mediated activation of Stat3 during multistage skin carcinogenesis. Cancer Res. 2004;64(7):2382–9.View ArticlePubMedGoogle Scholar
- Sumita N, Bito T, Nakajima K, Nishigori C. Stat3 activation is required for cell proliferation and tumorigenesis but not for cell viability in cutaneous squamous cell carcinoma cell lines. Exp Dermatol. 2006;15(4):291–9.View ArticlePubMedGoogle Scholar
- Yin W, Cheepala S, Roberts JN, Syson-Chan K, DiGiovanni J, Clifford JL. Active Stat3 is required for survival of human squamous cell carcinoma cells in serum-free conditions. Mol Cancer. 2006;5:15.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen SY, Takeuchi S, Moroi Y, Hayashida S, Kido M, Chen SJ, et al. Overexpression of phosphorylated-ATF2 and STAT3 in cutaneous squamous cell carcinoma, Bowen’s disease and basal cell carcinoma. J Dermatol Sci. 2008;51(3):210–5.View ArticlePubMedGoogle Scholar
- Zhao L, Li W, Marshall C, Griffin T, Hanson M, Hick R, et al. Srcasm inhibits Fyn-induced cutaneous carcinogenesis with modulation of Notch1 and p53. Cancer Res. 2009;69(24):9439–47.View ArticlePubMedPubMed CentralGoogle Scholar