Next Article in Journal
A Profound Basic Characterization of eIFs in Gliomas: Identifying eIF3I and 4H as Potential Novel Target Candidates in Glioma Therapy
Next Article in Special Issue
MicroRNAs Possibly Involved in the Development of Bone Metastasis in Clear-Cell Renal Cell Carcinoma
Previous Article in Journal
Abundance of Microvascular Endothelial Cells Is Associated with Response to Chemotherapy and Prognosis in Colorectal Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 6
10.3390/cancers13061480
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

miRNA Expression Characterizes Histological Subtypes and Metastasis in Penile Squamous Cell Carcinoma

by
Hiresh Ayoubian
1,
Joana Heinzelmann
1,2,
Sebastian Hölters
1,†,
Oybek Khalmurzaev
1,3,
Alexey Pryalukhin
4,‡,
Philine Loertzer
1,5,
Julia Heinzelbecker
1,
Stefan Lohse
6,
Carol Geppert
7,
Hagen Loertzer
5,
Heiko Wunderlich
8,
Rainer M. Bohle
4,
Michael Stöckle
1,
Vsevolod Borisovich Matveev
3,
Arndt Hartmann
7 and
Kerstin Junker
1,*
1
Department of Urology and Paediatric Urology, Saarland University, 66421 Homburg, Germany
2
Department of Ophthalmology, Martin-Luther-University Halle-Wittenberg, 06120 Halle/Saale, Germany
3
Department of Urology, Federal State Budgetary Institution “N.N. Blokhin National Medical Research Center of Oncology” of the Ministry of Health of the Russian Federation, 115478 Moscow, Russia
4
Institute of Pathology, Saarland University Medical Centre, 66421 Homburg, Germany
5
Westpfalz-Klinikum, Clinic of Urology and Paediatric Urology, 67655 Kaiserslautern, Germany
6
Institute of Virology, Saarland University, 66421 Homburg, Germany
7
Institute of Pathology, University Erlangen-Nuremberg, 91054 Erlangen, Germany
8
St. Georg Klinikum, Clinic of Urology and Paediatric Urology, 99817 Eisenach, Germany
*
Author to whom correspondence should be addressed.
†
Current affiliation: Ingenieurbüro für Gesundheitswesen GmbH, 04107 Leipzig, Germany.
‡
Current affiliation: Department of Pathology, Landesklinikum Wiener Neustadt, 2700 Wien, Austria.
Cancers 2021, 13(6), 1480; https://doi.org/10.3390/cancers13061480
Submission received: 2 February 2021 / Revised: 12 March 2021 / Accepted: 18 March 2021 / Published: 23 March 2021
(This article belongs to the Special Issue Non-Coding RNAs in Urological Malignancies: Relevance for Tumorigenesis and Biomarker Development)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Penile squamous cell carcinoma (PSCC) is the most common type of penile cancer (PeCa) and is associated with human papillomavirus (HPV) in about 50% of cases. It is of high clinical impact to identify patients who are at high risk of metastasis and are likely to benefit from adjuvant therapies. Today, valid prognostic biomarkers are scarce in penile cancer. In the present study, we attempted to identify miRNAs involved in tumor development and metastasis in distinct histological subtypes with or without HPV infection in PSCC patients. We confirmed that specific miRNAs could serve as potential diagnostic and prognostic markers in single PSCC subtypes and are associated with HPV infection.

Abstract

Although microRNAs are described as promising biomarkers in many tumor types, little is known about their role in PSCC. Thus, we attempted to identify miRNAs involved in tumor development and metastasis in distinct histological subtypes considering the impact of HPV infection. In a first step, microarray analyses were performed on RNA from formalin-fixed, paraffin-embedded tumor (22), and normal (8) tissue samples. Microarray data were validated for selected miRNAs by qRT-PCR on an enlarged cohort, including 27 tumor and 18 normal tissues. We found 876 significantly differentially expressed miRNAs (p ≤ 0.01) between HPV-positive and HPV-negative tumor samples by microarray analysis. Although no significant differences were detected between normal and tumor tissue in the whole cohort, specific expression patterns occurred in distinct histological subtypes, such as HPV-negative usual PSCC (95 differentially expressed miRNAs, p ≤ 0.05) and HPV-positive basaloid/warty subtypes (247 differentially expressed miRNAs, p ≤ 0.05). Selected miRNAs were confirmed by qRT-PCR. Furthermore, microarray data revealed 118 miRNAs (p ≤ 0.01) that were significantly differentially expressed in metastatic versus non-metastatic usual PSCC. The lower expression levels for miR-137 and miR-328-3p in metastatic usual PSCC were validated by qRT-PCR. The results of this study confirmed that specific miRNAs could serve as potential diagnostic and prognostic markers in single PSCC subtypes and are associated with HPV-dependent pathways.
Keywords:
PSCC; miRNA; HPV; microarray

1. Introduction

The incidence of penile cancer (PeCa) is low at 0.3 to 1.0 per 100,000 in Western Europe and the United States. The most common PeCa type is penile squamous cell carcinoma (PSCC) [1,2,3]. A recent systematic review and meta-analysis of the published literature found that human papillomavirus (HPV) DNA is detectable in 50.8% of penile cancer cases, of which HPV-16 is the most prevalent viral subtype, accounting for 63% of the HPV-attributable cases, followed by HPV-6 (8.1%), HPV-18 (6.9%), and, to a lesser extent, HPV-11, HPV-45, HPV-33, HPV-52, and HPV-31, among others [4]. As such, the 2016 WHO classification categorizes PSCC into HPV-related (such as basaloid, warty, and warty–basaloid) and non-HPV-related subtypes, which include the usual histological subtype as the most common [5,6]. Previous studies have explained two different mechanisms for oncogenesis in PSCC. First, the HPV-independent mechanism is linked to p53 mutation and shows nuclear p53 expression in atypical basal keratinocytes without P16INK4a expression [7,8,9]. P16INK4a downregulation through overexpression of the BMI-1 polycomb gene and inactivation of P16INK4a through DNA methylation has also been described in PSCC [10,11]. Second, the HPV-dependent mechanism involves E6 and E7 viral oncoproteins binding to p53 and retinoblastoma gene protein Rb. These host regulatory proteins disrupt DNA repair, growth arrest, apoptosis via E6 [12], and inactivation of hypo-phosphorylated Rb protein via E7, resulting in nuclear accumulation of cyclin-dependent kinase inhibitor P16INK4a [10,13,14,15].
In contrast to the well-known role of HPV in tumor development, the impact of HPV infection on the prognosis of PSCC remains unclear and is still under discussion [16,17,18].
It is of high clinical impact to identify patients who are at a high risk of metastasis and are likely to benefit from invasive therapies such as early inguinal lymph node dissection or adjuvant systemic chemotherapy. We and others proved regional lymph node involvement and lymph vessel invasion as independent prognostic factors in PSCC [19,20,21]. However, these parameters are not sufficient for individual therapy selection to improve treatment in high-risk patients and to prevent adverse effects in low-risk patients. Therefore, prognostic biomarkers are urgently needed.
miRNAs are conserved, short, single-stranded, non-coding RNAs of 19 to 25 nucleotides involved in post-transcriptional regulation of target genes by binding to the 3′ untranslated regions (3′UTRs) of the messenger RNA (mRNA) [22,23]. It has been shown that miRNAs play a crucial role in a broad range of biological processes, such as the proliferation, metabolism, apoptosis, and metastasis of tumors [24,25,26]. Changes in miRNA expression have been reported in a variety of human tumors and can potentially be employed as clinical biomarkers of diagnosis, prognosis, metastasis, and therapy resistance [27,28]. It was previously revealed that HPV infection via viral oncoproteins E6 and E7 altered the expression of cellular oncogenic and tumor suppressor miRNAs by modulating the c-Myc, p53, and E2F genes [29,30,31]. While it can be concluded that HPV infection has an effective impact on many biological processes, particularly through its interaction with the miRNA network, the impact of HPV infection on miRNA expression in PSCC remains unclear [16,17,18]. Furthermore, little is known about the role of miRNAs in tumor development and progression in this tumor type in general, and even less in histological subtypes. Therefore, miRNA expression analyses were performed in the present study to identify deregulated miRNAs associated with tumorigenesis and metastasis in distinct histological subtypes of PSCC considering HPV status.

2. Materials and Methods

2.1. Sample Collection and Patient Data

The patient samples were obtained from urology departments in Germany and Russia as described in our previous cohort study [19] (Table 1 and Table S1). Thirty samples, including 22 tumor tissue and 8 normal tissue samples, were used for screening by microarray analysis. For validation by PCR, 45 samples (27 tumor and 18 normal tissue samples) were analyzed, including those from the screening cohort. The hematoxylin and eosin (H&E)-stained tissue sections were assessed by two pathologists to confirm histologic subtypes as previously described [19]. Written informed consent was obtained from patients in most cases. For those individuals from whom we could not obtain informed consent, the data were analyzed anonymously. A regional ethics committee approved this investigation (Saarland ethical committee votes: 42/17, 220/19).

2.2. DNA and RNA Isolation

Total DNA and RNA (including small RNA) were isolated from formalin-fixed, paraffin-embedded (FFPE) tissue sections, prepared by punches from defined tumor-affected or normal areas, using AllPrep® DNA/RNA FFPE Kits (Qiagen, Hilden Germany) according to the manufacturer’s instructions. Detection of HPV infection was performed by HPV polymerase chain reaction (PCR) for GP5+/6+ and Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) primers using the LightCycler® 1.5 instrument (Roche Diagnostics GmbH) as previously described [19,32]. The concentration and purity of the RNA was assessed using the ratio of absorbance levels at 230, 260, and 280 nm, determined using a NanoDrop™ ND-1000 Spectrophotometer (ThermoFisher Scientific, Dreieich, Germany).

2.3. miRNA Microarray and Data Analysis

MicroRNA labelling and hybridization were performed using the miRNA Complete Labelling and Hybridization Kit (Agilent, Waldbronn, Germany) according to the manufacturer’s instructions. Briefly, 100 ng of isolated RNA from PSCC and normal tissue was dephosphorylated, labelled with pCp-Cy3, and hybridized to the Agilent human miRNA microarray. Microarray data analysis was performed using Qlucore Omics Explorer 3.5 (Qlucore AB, Lund, Sweden). In brief, microarrays were scanned using the Agilent Microarray Scanner, and expression values were extracted using the Agilent Feature Extraction (FE) software. Unpaired Student t-tests were performed using Qlucore Omics Explorer. The differential expression of miRNA between groups was considered statistically significant when p-values were less than 0.05, false discovery rates (q) were less than 0.3, and a 2-fold change in miRNA expression was detected. The total gene signal was normalized to the 75th percentile of signal intensity. Unsupervised hierarchical clustering and principal component analysis (PCA) were also conducted using Qlucore Omics Explorer.

2.4. Real-Time RT-PCR Quantification of miRNAs

To validate the miRNA microarray data, reverse transcription PCR (qRT-PCR) was performed using TaqMan® MicroRNA Reverse Transcription Kits or TaqMan® Advanced miRNA cDNA Synthesis Kits, followed by a polymerase chain reaction (PCR) step using specific TaqMan® miRNA primers and TaqMan® Fast Advanced Master Mix, according to the manufacturer’s protocol. PCRs were run in triplicate using LightCycler® 480 (Roche diagnostics Deutschland GmbH, Mannheim, Germany). miR-191-5p (ThermoFisher) was used as an endogenous control to normalize miRNA input in the real-time PCR (qRT-PCR) assay. The relative quantitative method of 2−ΔΔCT was used to measure changes in selected miRNAs. The ΔCt value was calculated by subtracting the Ct values of miR-191-5p from the Ct values of the miRNA of interest. Deregulation of miRNA was considered validated if there was at least a 2-fold expression difference between each group and the deregulation was in the same direction as in the microarray data. Ct values above 35 were considered as not expressing miRNA. The nonparametric Mann–Whitney U test was used to analyze the qRT-PCR data (IBM® SPSS® Statistics version 25). A p-value of less than 0.05 was considered a statistically significant threshold for upregulated or downregulated miRNA.

3. Results

3.1. Microarray Analysis

To determine the possible association between HPV infection and miRNA expression and to test whether miRNAs could serve as diagnostic or prognostic markers, 12 HPV-negative and 10 HPV-positive tumor samples of different histological subtypes, as well as 8 non-tumor tissue samples related to both HPV-positive and HPV-negative tumors, were subjected to miRNA microarray analysis. The comparison of HPV-positive and HPV-negative cases revealed significant differences in miRNA expression profiles (Figure 1a). In total, 876 miRNAs showed significantly different expression between HPV-positive and HPV-negative tumors (p ≤ 0.01, q ≤ 0.027), including 868 downregulated and 8 upregulated miRNAs in HPV-positive PSCC (Table S2). However, we could not find statistically significant differences in miRNA expression profiles between HPV-related and non-HPV-related histological subtypes (Figure S1) (q ≥ 0.81). There was also a statistically significant difference between non-tumor tissues based on HPV status (p ≤ 0.05, q ≤ 0.37), resulting in the upregulation of 24 miRNAs and the downregulation of 12 miRNAs in HPV-positive cases (Figure 1b). Surprisingly, we did not find significant differences in miRNA expression between tumor and normal samples in general (q ≥ 0.99) (Figure S2). From PCA and hierarchical clustering of the miRNA expression profiles of usual, warty, basaloid, and warty–basaloid PSCC subtypes, we found distinct miRNA profiles of basaloid tumors, and we found the usual subtype to be closer to the warty and warty–basaloid subtypes (F ≤ 0.01; Figure 2).
Based on these differences, we further investigated miRNA expression in the context of histological subtypes considering HPV status. A comparison of the profiles of the HPV-negative usual subtype and normal tissues revealed 95 differentially expressed miRNAs, including 91 that were upregulated and 4 that were downregulated in tumor tissue with a 2-fold change (p ≤ 0.05, q ≤ 0.23; Figure 3 and Table S3). In this comparison, miR-4512 had the highest expression, with a 5.1-fold increase, and miR-99a-5p had the lowest expression, with an 8-fold decrease, relative to normal tissue. To identify deregulated miRNAs involved in metastasis within this subgroup, the miRNA expression profiles in metastatic versus non-metastatic tissues were then analyzed (Figure 4 and Table S4). They were found to be statistically significantly different (p ≤ 0.01, q ≤ 0.08), with a total of 118 significantly deregulated miRNAs, of which 28 were upregulated and 90 were downregulated in metastatic tumors. Here, miR-4498 had the highest expression, with a 5.4-fold increase, whereas miR-137 (with a 3.7-fold decrease) and miR-301b (with a 4-fold decrease) had the lowest expression in usual metastatic tumors.
Due to the low frequency of basaloid tumors and the histological similarities between warty, basaloid, and warty–basaloid tumors, these three PSCC subtypes were combined into a single group. A comparison of this HPV-positive group to the related normal tissues revealed 247 miRNAs that were differently expressed (p ≤ 0.05, q ≤ 0.17). The majority of the miRNAs (n = 245) were downregulated, while only two were upregulated (Figure 5). Here, miR-181d-5p had the highest expression, with a 3.41-fold increase, whereas miR-211-5p had the lowest expression, with a 16.96-fold decrease in tumor tissue (Table S5). The analysis of microarray data did not reveal statistically significantly deregulated miRNAs between metastatic and non-metastatic tissues in these subtypes (q ≥ 0.79; Figure S3).

3.2. Validation of Differentially Expressed miRNAs by qRT-PCR

Based on the results of the microarray analysis (p- and q-values and fold change) and the potential biological functions in the context of tumorigenesis and metastasis, we selected eight miRNAs for validation using TaqMan® RT-qPCR from an enlarged cohort, including samples used for the microarray analyses. We focused on analysis in specific subtypes as described before. The significantly lower expression level for miR-99a-5p (p = 0.017) in the HPV-negative usual subtype compared to normal tissue samples was validated (Figure 6a). No statistically significant differences were found for either miR-125b-2-3p or miR-105-5p in this subgroup comparison (Figure S4), despite significant differences in microarray analysis. Furthermore, significantly lower expression levels for miR-137 (p = 0.004) and miR-328-3p (p = 0.032) in metastatic usual subtypes compared to non-metastatic tumors were also confirmed (Figure 6b,c). No significant difference was found for miR-509-5p in this subgroup comparison (Figure S4). Finally, we validated the higher expression of miR-181d-5p (p = 0.002) (Figure 6d), and no or very low expression of miR-211-5p was found in the HPV-positive basaloid, warty, and warty–basaloid subtypes compared to normal tissue (Table 2). Due to the non-significant results of microarray analysis when comparing the metastatic and non-metastatic tumors, we did not perform further qRT-PCR experiments.

4. Discussion

Currently, no biomarkers are available for the differential diagnosis, prognostic evaluation, or therapy response prediction of PSCC. However, this is a clinical need since individual prognostic evaluation and, therefore, treatment decisions concerning prevention of or invasive therapy for metastatic dissemination are still a challenge in this rare but often aggressive disease.
HPV infection is one of the most common risk factors in PSCC [1]. The prognostic role of HPV, however, is still debated [19,33,34,35,36]. The current WHO classification of PSCC is now based on clinicopathologic uniqueness and HPV status, and it assigns histological subtypes to HPV-related and non-HPV-related PSCC [5]. It should be considered that there is great diversity among histopathological subtypes in the two major groups concerning aggressiveness and outcome. Whereas verrucous and pseudohyperplastic PSCC are well differentiated and have an excellent prognosis, pseudoglandular and sarcomatoid PSCC are aggressive tumors with poor outcomes, and all are non-HPV-related PSCC. Usual PSCC represents the most frequent non-HPV-related subtype with heterogeneous outcome. In HPV-related PSCC, the basaloid subtype represents a poorly differentiated tumor type with a mainly worse disease course. Therefore, individual histopathological subtypes should be considered in clinical and molecular studies in addition to the HPV-status.
In the present study, we investigated miRNA expression profiles considering histopathological subtypes and HPV status. Microarray analysis showed that 876 miRNAs are differentially expressed at statistically significant levels in HPV-positive relative to HPV-negative PSCC. This suggests that HPV infection may indeed affect the expression of miRNAs in PSCC, resulting in distinct molecular pathways, as shown previously [36,37,38]. Interestingly, a positive HPV status in normal tissue was also associated with changes in miRNA expression. This could explain why statistically significant differences in miRNA expression profiles between tumor and normal tissues in general were not identified. In contrast to the real HPV status, a comparison of HPV-related (subtypes with predominant but not exclusive HPV infection based on histological subtypes) and non-HPV-related (subtypes with predominantly negative HPV status with only some HPV infection based on histological subtypes) PSCC did not reveal statistically significant changes in miRNA expression as defined by the current WHO classification system [5]. In concordance with published data, we found that 23% of HPV-negative PSCC are classified as HPV-related based on histological subtype. These data suggest that real HPV infection, not classification in HPV-related groups, is the relevant biological event associated with miRNA expression.
As our previous findings support the assertion that histological subtypes are more relevant than HPV status for prognosis [19], miRNA expression in different histological subtypes was further analyzed. Since PSCC is a rare tumor, little is known about miRNA alterations in distinct histological PSCC subtypes. The results of the present study confirm that histological subtypes are characterized by distinct miRNA expression profiles, especially between basaloid and usual PSCC. Therefore, the miRNA expression in these two main histological subtypes was further separately investigated. In contrast to the whole cohort, statistically significant differences in miRNA expression were identified between normal and tumor tissues in the HPV-negative usual and HPV-positive basaloid, warty, and warty–basaloid PSCC subtypes. To validate the microarray data, quantification via qRT-PCR was performed for selected miRNAs. We confirmed the lower expression of miR-99a-5p in tumor tissues compared to normal tissues in the usual PSCC subtype, supporting data on the tumor-suppressive role of this miRNA in other tumor types [39,40,41].
In the miRNA expression profile comparison of the warty, basaloid, and warty–basaloid group versus normal tissue, miR-211-5p and miR-181d-5p were the most downregulated and upregulated miRNAs, respectively. Interestingly, miR-211-5p expression in this group was not detected at all in tumor tissues based on the cut-off criteria in qRT-PCR analysis. These results are in line with the microarray data, which found a 16-fold downregulation in miR-211-5p in this group compared to normal samples. These data are consistent with a previous investigation of small RNA sequencing that found a 300-fold lower expression of miR-211-5p in cancerous penile tissues relative to adjacent non-cancerous tissues [42]. miR-211-5p has been reported to play a crucial role as a tumor suppressor in various cancer types [43,44,45,46]
The results of the qRT-PCR analysis also confirmed the overexpression of miR-181d-5p. The upregulation of miR-181d-5p has been found to induce hepatocarcinogenesis [47]. This miRNA packed in exosomes secreted by cancer associated fibroblasts enhances the epithelial–mesenchymal transition (EMT) and aggressiveness of breast cancer by targeting caudal-related homeobox 2 (CDX2) [48]. CDX2 loss or reduced expression has been associated with advanced tumor stage, including metastasis, as well as poor prognosis in colorectal cancer [49,50]. It is therefore possible that miR-181d-5p overexpression directly or indirectly reduces the expression of the CDX2 protein in PSCC, enforcing tumorigenesis.
We were then interested in finding putative prognostic miRNAs. We could not detect miRNA expression differences in association with metastasis in the whole cohort (data not shown). Therefore, we analyzed altered miRNAs possibly involved in metastasis by comparing metastatic and non-metastatic HPV-negative usual tumors. We found a clear pattern that characterized metastatic primary tumors by microarray analysis. The lower expression of miR-137 and miR-328-3p was confirmed by qRT-PCR. It has been shown that miR-137, as a tumor suppressor, is downregulated in many cancer types [51,52,53]. miR-137 overexpression inhibits tumor growth and increases chemosensitivity to paclitaxel and cisplatin in lung cancer cell lines [54]. In this regard, it would be interesting to investigate the putative role of this miRNA for chemotherapy response in PSCC, too. Downregulation of miR-328-3p has been linked to its role as a tumor metastasis suppressor by targeting matrix metalloproteases [55,56,57]. Zhang et al. also showed downregulation of miR-328-3p in penile cancer tissues relative to adjacent non-cancerous tissues by Next-Generation Small RNA Sequencing [42]. MiRmiR-328-3p inhibits metastasis in colorectal cancer via inactivation of the PI3K/Akt signaling pathway [58] and reduces invasion as well as EMT in liver cancer via by targeting endoplasmatic reticulum metallo protease 1 (ERMP1) to inhibit AKT phosphorylation [59].
We also investigated differential miRNA expression between metastatic and non-metastatic tumors in combined basaloid, warty, and warty–basaloid subtypes. However, the heterogeneity in this group was greatly elevated, as reflected by the high false discovery rate (q ≥ 0.79), perhaps due to the small sample size and heterogeneity. Therefore, the results were not sufficient for validation using qRT-PCR in this study, and this will be explored in the future with a larger sample cohort. Our data confirmed that distinct miRNAs could serve as prognostic biomarkers in primary tumor tissues to predict metastasis and could therefore support decision-making for local lymph node treatment, as well as adjuvant therapy.
It is accepted that the tumor microenvironment plays an important role in tumor development and progression. In this regard, immune checkpoints are of high interest not only as prognostic markers, but also as therapeutic targets. High programmed death-ligand 1 (PD-L1) expression is associated with poor prognosis in a variety of solid tumors [60,61,62,63]. Forty-eight percent of PSCC cases are PD-L1-positive; these are primarily associated with HPV-negative subtypes and poor survival [64]. PD-L1 expression is associated with lymph node metastasis in PSCC patients [64,65]. The correlation between miRNA expression—including that of miR-138-5p, miR-15a, miR-15b, miR-16, miR-34a, miR-142-5p, and miR-200—and PD-L1 is well known [60,66,67,68,69]. We found by microarray analysis that miR-138-5p was upregulated in usual PSCC. miR-138-5p is a tumor suppressor that targets PD-L1 [70,71,72]. Therefore, it would be interesting to investigate the potential role of miR-138-5p as a biomarker for therapy selection using checkpoint inhibitors in PSCC patients in the future.
In general, few data on miRNA expression in PSCC are published. Higher expression levels of miR-223-3p, miR-107, and miR-21 [36,42,73], and lower expression levels of miR-1, miR-101, and miR-204, have been reported as potential biomarkers in patients with PSCC [74]. However, we did not find significant expression differences for the mentioned miRNAs, which might be due to several reasons, including technical differences. We used punches from defined tumor and normal areas obtained from FFPE blocks to isolate RNA. miRNA expression was analyzed by microarrays instead of NGS with different bioinformatics tools. In addition, possibly more importantly, we investigated the altered miRNAs in respect to histological subtypes separately. Further studies should also include long noncoding RNAs as another important class of non-coding RNAs regulating the complex network of RNA translation.

5. Conclusions

This is the first study, to our knowledge, investigating miRNA expression in PSCC considering histological subtypes, HPV infection, and metastatic status in a single study. The identified miRNAs could serve as potential diagnostic and prognostic markers for single PSCC subtypes. Our data support the assertion that PSCC subtypes are characterized by distinct molecular alterations and should therefore be investigated separately, as already shown by clinical data [19]. It seems that HPV infection induces distinct miRNA expression in PSCC. The influence and degree of altered miRNAs associated with HPV infection in the different PSCC histological subtypes requires additional research in a larger cohort.

Supplementary Materials

The following data are available online at https://www.mdpi.com/2072-6694/13/6/1480/s1: Figure S1: Unsupervised hierarchical clustering of differentially expressed miRNAs in HPV-related versus non-HPV-related PSCC, Figure S2: Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in tumor versus normal tissues, Figure S3: Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in metastatic versus non- metastatic (HPV-positive) basaloid, warty, and warty basaloid histological subtype, Figure S4: miRNAs without significant differences as determined by TaqMan® RT-qPCR, Table S1: Summary of clinical and histopathological data of the PSCC samples used for microarray analysis and qRT-PCR validation, Table S2: Significantly differentially expressed miRNAs in HPV-positive versus HPV-negative PSCC, Table S3: Significantly differentially expressed miRNAs in usual (HPV-negative) PSCC versus normal samples, Table S4: Significantly differentially expressed miRNAs in metastatic versus non-metastatic usual (HPV-negative) PSCC, Table S5: Significantly differentially expressed miRNAs in combined basaloid and warty basaloid (HPV-positive) PSCC versus normal tissue samples (HPV-positive).

Author Contributions

Conceptualization K.J. and M.S.; methodology H.A., J.H. (Joana Heinzelmann), S.H., O.K. and S.L.; experimental analysis H.A., J.H. (Joana Heinzelmann), S.H.; data and sample acquisition O.K., A.P., H.L., P.L., H.W., V.B.M., R.M.B., A.H., C.G.; writing—original draft preparation, H.A. and K.J.; writing—review and editing, K.J., M.S., S.L., J.H. (Joana Heinzelmann), J.H. (Julia Heinzelbecker); supervision K.J.; funding acquisition J.H. (Joana Heinzelmann), K.J., M.S. All authors have read and agreed to the published version of the manuscript.

Funding

Medical Faculty: Saarland University: HOMFOR project initiation funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Saarland Ethics Committee (220/19, 14/02/2019).

Informed Consent Statement

Informed consent was obtained from subjects involved in the study. If informed consent could not be obtained data were analyzed anonymously in accordance with ethical guidelines.

Data Availability Statement

The data presented in this study are available in supplementary material tables and on request from the corresponding author.

Acknowledgments

The authors thank Alexander Vogt and Maria Link for excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bleeker, M.C.G.; Heideman, D.A.M.; Snijders, P.J.F.; Horenblas, S.; Dillner, J.; Meijer, C.J.L.M. Penile cancer: Epidemiology, pathogenesis and prevention. World J. Urol. 2009, 27, 141–150. [Google Scholar] [CrossRef]
  2. Jemal, A.; Siegel, R.; Ward, E.; Murray, T.; Xu, J.; Thun, M.J. Cancer Statistics. CA A Cancer J. Clin. 2007, 57, 43–66. [Google Scholar] [CrossRef]
  3. Parkin, D.M.; Muir, C.S. Cancer Incidence in Five Continents. Comparability and quality of data. IARC Sci. Publ. 1992, 120, 45–173. [Google Scholar]
  4. Olesen, T.B.; Sand, F.L.; Rasmussen, C.L.; Albieri, V.; Toft, B.G.; Norrild, B.; Munk, C.; Kjær, S.K. Prevalence of human papillomavirus DNA and p16INK4a in penile cancer and penile intraepithelial neoplasia: A systematic review and meta-analysis. Lancet Oncol. 2019, 20, 145–158. [Google Scholar] [CrossRef]
  5. Moch, H.; Cubilla, A.L.; Humphrey, P.A.; Reuter, V.E.; Ulbright, T.M. The 2016 WHO Classification of Tumours of the Urinary System and Male Genital Organs—Part A: Renal, Penile, and Testicular Tumours. Eur. Urol. 2016, 70, 93–105. [Google Scholar] [CrossRef]
  6. Miralles-Guri, C.; Bruni, L.; Cubilla, A.L.; Castellsagué, X.; Bosch, F.X.; De Sanjosé, S. Human papillomavirus prevalence and type distribution in penile carcinoma. J. Clin. Pathol. 2009, 62, 870–878. [Google Scholar] [CrossRef]
  7. Chaux, A.; Pfannl, R.; Rodríguez, I.M.; Barreto, J.E.; Velazquez, E.F.; Lezcano, C.; Piris, A.; Netto, G.J.; Cubilla, A.L. Distinctive Immunohistochemical Profile of Penile Intraepithelial Lesions. Am. J. Surg. Pathol. 2011, 35, 553–562. [Google Scholar] [CrossRef]
  8. Pinto, A.P.; Miron, A.; Yassin, Y.; Monte, N.; Woo, T.Y.C.; Mehra, K.K.; Medeiros, F.; Crum, C.P. Differentiated vulvar intraepithelial neoplasia contains Tp53 mutations and is genetically linked to vulvar squamous cell carcinoma. Mod. Pathol. 2010, 23, 404–412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Rocha, R.M.; Ignácio, J.A.; Jordán, J.; Carraro, D.M.; Lisboa, B.; Lopes, A.; Carvalho, K.C.; Da Cunha, I.W.; Cubilla, A.; Guimarães, G.C.; et al. A clinical, pathologic, and molecular study of p53 and murine double minute 2 in penile carcinogenesis and its relation to prognosis. Hum. Pathol. 2012, 43, 481–488. [Google Scholar] [CrossRef]
  10. Ferreux, E.; Lont, A.P.; Horenblas, S.; Gallee, M.P.W.; Raaphorst, F.M.; Doeberitz, M.V.K.; Meijer, C.J.L.M.; Snijders, P.J.F. Evidence for at least three alternative mechanisms targeting the p16INK4A/cyclin D/Rb pathway in penile carcinoma, one of which is mediated by high-risk human papillomavirus. J. Pathol. 2003, 201, 109–118. [Google Scholar] [CrossRef]
  11. Emmanuel, A.; Nettleton, J.; Watkin, N.; Berney, D.M. The molecular pathogenesis of penile carcinoma—current developments and understanding. Virchows Archiv. 2019, 475, 397–405. [Google Scholar] [CrossRef]
  12. Vanthoor, J.; Vos, G.; Albersen, M. Penile cancer: Potential target for immunotherapy? World J. Urol. 2020, 1–7. [Google Scholar] [CrossRef]
  13. Mannweiler, S.; Sygulla, S.; Winter, E.; Regauer, S. Two major pathways of penile carcinogenesis: HPV-induced penile cancers overexpress p16ink4a, HPV-negative cancers associated with dermatoses express p53, but lack p16ink4a overexpression. J. Am. Acad. Dermatol. 2013, 69, 73–81. [Google Scholar] [CrossRef]
  14. Afonso, L.A.; Carestiato, F.N.; Ornellas, A.A.; Ornellas, P.; Rocha, W.M.; Cordeiro, T.I.; Lisboa, D.C.; Alves, G.B.; Cavalcanti, S.M.B. Human papillomavirus, Epstein-Barr virus, and methylation status of p16ink4ain penile cancer. J. Med Virol. 2017, 89, 1837–1843. [Google Scholar] [CrossRef]
  15. Gage, J.R.; Meyers, C.; O Wettstein, F. The E7 proteins of the nononcogenic human papillomavirus type 6b (HPV-6b) and of the oncogenic HPV-16 differ in retinoblastoma protein binding and other properties. J. Virol. 1990, 64, 723–730. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, F.; Li, B.; Xie, X. The roles and clinical significance of microRNAs in cervical cancer. Histol. Histopathol. 2015, 31. [Google Scholar]
  17. Zheng, Z.-M.; Wang, X. Regulation of cellular miRNA expression by human papillomaviruses. Biochim. Biophys. Acta (BBA) Bioenerg. 2011, 1809, 668–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Snoek, B.C.; Babion, I.; Koppers-Lalic, D.; Pegtel, D.M.; Steenbergen, R.D. Altered microRNA processing proteins in HPV-induced cancers. Curr. Opin. Virol. 2019, 39, 23–32. [Google Scholar] [CrossRef]
  19. Hölters, S.; Khalmurzaev, O.; Pryalukhin, A.; Loertzer, P.; Janssen, M.; Heinzelbecker, J.; Ueberdiek, S.; Pfuhl, T.; Smola, S.; Agaimy, A.; et al. Challenging the prognostic impact of the new WHO and TNM classifications with special emphasis on HPV status in penile carcinoma. Virchows Archiv. 2019, 475, 211–221. [Google Scholar] [CrossRef]
  20. Sun, M.; Djajadiningrat, R.S.; Alnajjar, H.M.; Trinh, Q.-D.; Graafland, N.M.; Watkin, N.; Karakiewicz, P.I.; Horenblas, S. Development and external validation of a prognostic tool for prediction of cancer-specific mortality after complete loco-regional pathological staging for squamous cell carcinoma of the penis. BJU Int. 2014, 116, 734–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Paner, G.P.; Stadler, W.M.; Hansel, D.E.; Montironi, R.; Lin, D.W.; Amin, M.B. Updates in the Eighth Edition of the Tumor-Node-Metastasis Staging Classification for Urologic Cancers. Eur. Urol. 2018, 73, 560–569. [Google Scholar] [CrossRef]
  22. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350–355. [Google Scholar] [CrossRef] [PubMed]
  23. Bartel, D.P. MicroRNAs. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  24. Esquela-Kerscher, A.; Slack, F.J. Oncomirs—microRNAs with a role in cancer. Nat. Rev. Cancer 2006, 6, 259–269. [Google Scholar] [CrossRef]
  25. He, L.; Hannon, G.J. MicroRNAs: Small RNAs with a big role in gene regulation. Nat. Rev. Genet. 2004, 5, 522–531. [Google Scholar] [CrossRef]
  26. Bala, S.; Marcos, M.; Szabo, G. Emerging role of microRNAs in liver diseases. World J. Gastroenterol. 2009, 15, 5633–5640. [Google Scholar] [CrossRef]
  27. Lee, K.T.-W.; Tan, J.-K.; Lam, A.K.-Y.; Gan, S.-Y. MicroRNAs serving as potential biomarkers and therapeutic targets in nasopharyngeal carcinoma: A critical review. Crit. Rev. Oncol. 2016, 103, 1–9. [Google Scholar] [CrossRef] [Green Version]
  28. Cho, W.C. MicroRNAs: Potential biomarkers for cancer diagnosis, prognosis and targets for therapy. Int. J. Biochem. Cell Biol. 2010, 42, 1273–1281. [Google Scholar] [CrossRef] [PubMed]
  29. Wald, A.I.; Hoskins, E.E.; Wells, S.I.; Ferris, R.L.; Khan, S.A. Human papillomavirus alters microRNA profiles in squamous cell carcinoma of the head and neck (SCCHN) cell lines. Head Neck 2011, 33, 504–512. [Google Scholar] [CrossRef] [Green Version]
  30. Gao, D.; Zhang, Y.; Zhu, M.; Liu, S.; Wang, X. miRNA Expression Profiles of HPV-Infected Patients with Cervical Cancer in the Uyghur Population in China. PLoS ONE 2016, 11, e0164701. [Google Scholar] [CrossRef] [Green Version]
  31. Santos, J.M.; Da Silva, S.P.; Costa, N.R.; Gil Da Costa, R.M.; Medeiros, R. The Role of MicroRNAs in the Metastatic Process of High-Risk HPV-Induced Cancers. Cancers 2018, 10, 493. [Google Scholar] [CrossRef] [Green Version]
  32. Husman, A.-M.D.R.; Walboomers, J.M.M.; Brule, A.J.C.V.D.; Meijer, C.J.L.M.; Snijders, P.J.F. The use of general primers GP5 and GP6 elongated at their 3′ ends with adjacent highly conserved sequences improves human papillomavirus detection by PCR. J. Gen. Virol. 1995, 76, 1057–1062. [Google Scholar] [CrossRef] [PubMed]
  33. Djajadiningrat, R.S.; Jordanova, E.S.; Kroon, B.K.; Van Werkhoven, E.; De Jong, J.; Pronk, D.T.; Snijders, P.J.; Horenblas, S.; Heideman, D.A. Human Papillomavirus Prevalence in Invasive Penile Cancer and Association with Clinical Outcome. J. Urol. 2015, 193, 526–531. [Google Scholar] [CrossRef] [PubMed]
  34. Gregoire, L.; Cubilla, A.L.; Reuter, V.E.; Haas, G.P.; Lancaster, W.D. Preferential Association of Human Papillomavirus With High-Grade Histologic Variants of Penile-Invasive Squamous Cell Carcinoma. J. Natl. Cancer Inst. 1995, 87, 1705–1709. [Google Scholar] [CrossRef] [PubMed]
  35. Lont, A.P.; Kroon, B.K.; Horenblas, S.; Gallee, M.P.W.; Berkhof, J.; Meijer, C.J.L.M.; Snijders, P.J.F. Presence of high-risk human papillomavirus DNA in penile carcinoma predicts favorable outcome in survival. Int. J. Cancer 2006, 119, 1078–1081. [Google Scholar] [CrossRef] [PubMed]
  36. Kuasne, H.; Barros-Filho, M.C.; Busso-Lopes, A.; Marchi, F.A.; Pinheiro, M.; Muñoz, J.J.M.; Scapulatempo-Neto, C.; Faria, E.F.; Guimarães, G.C.; Lopes, A.; et al. Integrative miRNA and mRNA analysis in penile carcinomas reveals markers and pathways with potential clinical impact. Oncotarget 2017, 8, 15294–15306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Barzon, L.; Cappellesso, R.; Peta, E.; Militello, V.; Sinigaglia, A.; Fassan, M.; Simonato, F.; Guzzardo, V.; Ventura, L.; Blandamura, S.; et al. Profiling of Expression of Human Papillomavirus–Related Cancer miRNAs in Penile Squamous Cell Carcinomas. Am. J. Pathol. 2014, 184, 3376–3383. [Google Scholar] [CrossRef] [PubMed]
  38. Peta, E.; Cappellesso, R.; Masi, G.; Sinigaglia, A.; Trevisan, M.; Grassi, A.; Di Camillo, B.; Vassarotto, E.; Fassina, A.; Palù, G.; et al. Down-regulation of microRNA-146a is associated with high-risk human papillomavirus infection and epidermal growth factor receptor overexpression in penile squamous cell carcinoma. Hum. Pathol. 2017, 61, 33–40. [Google Scholar] [CrossRef] [PubMed]
  39. Tsai, T.-F.; Lin, J.-F.; Chou, K.-Y.; Lin, Y.-C.; Chen, H.-E.; Hwang, T.I.-S. miR-99a-5p acts as tumor suppressor via targeting to mTOR and enhances RAD001-induced apoptosis in human urinary bladder urothelial carcinoma cells. OncoTargets Ther. 2018, 11, 239–252. [Google Scholar] [CrossRef] [Green Version]
  40. Chen, Y.-T.; Yao, J.-N.; Qin, Y.-T.; Hu, K.; Wu, F.; Fang, Y.-Y. Biological role and clinical value of miR-99a-5p in head and neck squamous cell carcinoma (HNSCC): A bioinformatics-based study. FEBS Open Bio 2018, 8, 1280–1298. [Google Scholar] [CrossRef]
  41. Saito, R.; Maruyama, S.; Kawaguchi, Y.; Akaike, H.; Shimizu, H.; Furuya, S.; Kawaida, H.; Ichikawa, D. miR-99a-5p as Possible Diagnostic and Prognostic Marker in Patients With Gastric Cancer. J. Surg. Res. 2020, 250, 193–199. [Google Scholar] [CrossRef] [PubMed]
  42. Zhang, L.; Wei, P.; Shen, X.; Zhang, Y.; Xu, B.; Zhou, J.; Fan, S.; Hao, Z.; Shi, H.; Zhang, X.; et al. MicroRNA Expression Profile in Penile Cancer Revealed by Next-Generation Small RNA Sequencing. PLoS ONE 2015, 10, e0131336. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, W.; Liu, Z.; Zhang, X.; Liu, J.; Gui, J.; Cui, M.; Li, Y. miR-211-5p is down-regulated and a prognostic marker in bladder cancer. J. Gene Med. 2020, 22, 1–10. [Google Scholar] [CrossRef]
  44. Wang, K.; Jin, W.; Jin, P.; Fei, X.; Wang, X.; Chen, X. miR-211-5p Suppresses Metastatic Behavior by Targeting SNAI1 in Renal Cancer. Mol. Cancer Res. 2017, 15, 448–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Quan, J.; Pan, X.; He, T.; Lin, C.; Lai, Y.; Chen, P.; Zhang, Z.; Yang, S.; Wang, T.; Lai, Y. Tumor suppressor miR-211-5p is associated with cellular migration, proliferation and apoptosis in renal cell carcinoma. Exp. Ther. Med. 2018, 15, 4019–4028. [Google Scholar] [CrossRef]
  46. Chen, L.-L.; Zhang, Z.-J.; Yi, Z.-B.; Li, J.-J. MicroRNA-211-5p suppresses tumour cell proliferation, invasion, migration and metastasis in triple-negative breast cancer by directly targeting SETBP. Br. J. Cancer 2017, 117, 78–88. [Google Scholar] [CrossRef]
  47. Wang, B.; Hsu, S.-H.; Majumder, S.K.; Kutay, H.; Huang, W.; Jacob, S.T.; Ghoshal, K. TGFβ-mediated upregulation of hepatic miR-181b promotes hepatocarcinogenesis by targeting TIMP. Oncogene 2009, 29, 1787–1797. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, H.; Wei, H.; Wang, J.; Li, L.; Chen, A.; Li, Z. MicroRNA-181d-5p-Containing Exosomes Derived from CAFs Promote EMT by Regulating CDX2/HOXA5 in Breast Cancer. Mol. Ther. Nucleic Acids 2020, 19, 654–667. [Google Scholar] [CrossRef] [PubMed]
  49. Graule, J.; Uth, K.; Fischer, E.; Centeno, I.; Galván, J.A.; Eichmann, M.; Rau, T.T.; Langer, R.; Dawson, H.; Nitsche, U.; et al. CDX2 in colorectal cancer is an independent prognostic factor and regulated by promoter methylation and histone deacetylation in tumors of the serrated pathway. Clin. Epigenetics 2018, 10, 120. [Google Scholar] [CrossRef] [Green Version]
  50. Zhang, B.Y.; Jones, J.C.; Briggler, A.M.; Hubbard, J.M.; Kipp, B.R.; Sargent, D.J.; Dixon, J.G.; Grothey, A. Lack of Caudal-Type Homeobox Transcription Factor 2 Expression as a Prognostic Biomarker in Metastatic Colorectal Cancer. Clin. Color. Cancer 2017, 16, 124–128. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, L.-L.; Lu, S.-X.; Li, M.; Li, L.-Z.; Fu, J.; Hu, W.; Yang, Y.-Z.; Luo, R.-Z.; Zhang, C.Z.; Yun, J.-P. FoxD3-regulated microRNA-137 suppresses tumour growth and metastasis in human hepatocellular carcinoma by targeting AKT. Oncotarget 2014, 5, 5113–5124. [Google Scholar] [CrossRef] [Green Version]
  52. Zhu, L.-L.; Chen, W.; Wang, J.; Gan, T.; Wang, Y.-P.; Yang, J.-L. Expression Difference of MicroRNAs in Colorectal Adenoma and Colorectal Cancer. Sichuan Da Xue Xue Bao Yi Xue Ban J. Sichuan Univ. Med Sci. Ed. 2015, 46, 851–855. [Google Scholar]
  53. Chen, L.; Wang, X.; Wang, H.; Li, Y.; Yan, W.; Han, L.; Zhang, K.; Zhang, J.; Wang, Y.; Feng, Y.; et al. miR-137 is frequently down-regulated in glioblastoma and is a negative regulator of Cox. Eur. J. Cancer 2012, 48, 3104–3111. [Google Scholar] [CrossRef] [PubMed]
  54. Shen, H.; Wang, L.; Ge, X.; Jiang, C.-F.; Shi, Z.-M.; Li, D.-M.; Liu, W.-T.; Yu, X.; Shu, Y.-Q. MicroRNA-137 inhibits tumor growth and sensitizes chemosensitivity to paclitaxel and cisplatin in lung cancer. Oncotarget 2016, 7, 20728–20742. [Google Scholar] [CrossRef]
  55. Zhang, M.; Zhang, J.; Zhou, Q. Elevated expression of microRNA-328-3p suppresses aggressive malignant behaviors via targeting matrix metalloprotease 16 in osteosarcoma. OncoTargets Ther. 2019, 12, 2063–2070. [Google Scholar] [CrossRef] [Green Version]
  56. Yang, S.-F.; Lee, W.-J.; Tan, P.; Tang, C.-H.; Hsiao, M.; Hsieh, F.-K.; Chien, M.-H. Upregulation of miR-328 and inhibition of CREB-DNA-binding activity are critical for resveratrol-mediated suppression of matrix metalloproteinase-2 and subsequent metastatic ability in human osteosarcomas. Oncotarget 2014, 6, 2736–2753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Li, J.-Z.; Li, J.; Liu, B.-Z. MicroRNA-328-3p inhibits malignant progression of hepatocellular carcinoma by regulating MMP-9 level. Eur. Rev. Med Pharmacol. Sci. 2019, 23, 9331–9340. [Google Scholar]
  58. Pan, S.; Ren, F.; Li, L.; Liu, D.; Li, Y.; Wang, A.; Li, W.; Dong, Y.; Guo, W. MiR-328-3p inhibits cell proliferation and metastasis in colorectal cancer by targeting Girdin and inhibiting the PI3K/Akt signaling pathway. Exp. Cell Res. 2020, 390, 111939. [Google Scholar] [CrossRef]
  59. Lu, H.; Hu, J.; Li, J.; Lu, W.; Deng, X.; Wang, X. miR-328-3p overexpression attenuates the malignant proliferation and invasion of liver cancer via targeting Endoplasmic Reticulum Metallo Protease 1 to inhibit AKT phosphorylation. Ann. Transl. Med. 2020, 8, 754. [Google Scholar] [CrossRef] [PubMed]
  60. Wang, Q.; Lin, W.; Tang, X.; Li, S.; Guo, L.; Lin, Y.; Kwok, H.F. The Roles of microRNAs in Regulating the Expression of PD-1/PD-L1 Immune Checkpoint. Int. J. Mol. Sci. 2017, 18, 2540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Zhou, C.; Tang, J.; Sun, H.; Zheng, X.; Li, Z.; Sun, T.; Li, J.; Wang, S.; Zhou, X.; Sun, H.; et al. PD-L1 expression as poor prognostic factor in patients with non-squamous non-small cell lung cancer. Oncotarget 2017, 8, 58457–58468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Lin, Y.-M.; Sung, W.-W.; Hsieh, M.-J.; Tsai, S.-C.; Lai, H.-W.; Yang, S.-M.; Shen, K.-H.; Chen, M.-K.; Lee, H.; Yeh, K.-T.; et al. High PD-L1 Expression Correlates with Metastasis and Poor Prognosis in Oral Squamous Cell Carcinoma. PLoS ONE 2015, 10, e0142656. [Google Scholar] [CrossRef] [Green Version]
  63. Gu, X.; Dong, M.; Liu, Z.; Mi, Y.; Yang, J.; Zhang, Z.; Liu, K.; Jiang, L.; Zhang, Y.; Dong, S.; et al. Elevated PD-L1 expression predicts poor survival outcomes in patients with cervical cancer. Cancer Cell Int. 2019, 19, 146. [Google Scholar] [CrossRef] [PubMed]
  64. Ottenhof, S.R.; Djajadiningrat, R.S.; De Jong, J.; Thygesen, H.H.; Horenblas, S.; Jordanova, E.S. Expression of Programmed Death Ligand 1 in Penile Cancer is of Prognostic Value and Associated with HPV Status. J. Urol. 2017, 197, 690–697. [Google Scholar] [CrossRef]
  65. Udager, A.M.; Liu, T.-Y.; Skala, S.L.; Magers, M.J.; McDaniel, A.S.; Spratt, D.E.; Feng, F.Y.; Siddiqui, J.; Cao, X.; Fields, K.L.; et al. Frequent PD-L1 expression in primary and metastatic penile squamous cell carcinoma: Potential opportunities for immunotherapeutic approaches. Ann. Oncol. 2016, 27, 1706–1712. [Google Scholar] [CrossRef] [Green Version]
  66. Kao, S.C.; Cheng, Y.Y.; Williams, M.; Kirschner, M.B.; Madore, J.; Lum, T.; Sarun, K.H.; Linton, A.; McCaughan, B.; Klebe, S.; et al. Tumor Suppressor microRNAs Contribute to the Regulation of PD-L1 Expression in Malignant Pleural Mesothelioma. J. Thorac. Oncol. 2017, 12, 1421–1433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Wang, X.; Li, J.; Dong, K.; Lin, F.; Long, M.; Ouyang, Y.; Wei, J.; Chen, X.; Weng, Y.; He, T.; et al. Tumor suppressor miR-34a targets PD-L1 and functions as a potential immunotherapeutic target in acute myeloid leukemia. Cell. Signal. 2015, 27, 443–452. [Google Scholar] [CrossRef]
  68. Jia, L.; Xi, Q.; Wang, H.; Zhang, Z.; Liu, H.; Cheng, Y.; Guo, X.; Zhang, J.; Zhang, Q.; Zhang, L.; et al. miR-142-5p regulates tumor cell PD-L1 expression and enhances anti-tumor immunity. Biochem. Biophys. Res. Commun. 2017, 488, 425–431. [Google Scholar] [CrossRef]
  69. Chen, L.; Gibbons, D.L.; Goswami, S.; Cortez, M.A.; Ahn, Y.-H.; Byers, L.A.; Zhang, X.; Yi, X.; Dwyer, D.; Lin, W.; et al. Metastasis is regulated via microRNA-200/ZEB1 axis control of tumour cell PD-L1 expression and intratumoral immunosuppression. Nat. Commun. 2014, 5, 5241. [Google Scholar] [CrossRef]
  70. Sha, H.-H.; Wang, D.-D.; Chen, D.; Liu, S.-W.; Wang, Z.; Yan, D.-L.; Dong, S.-C.; Feng, J.-F. MiR-138: A promising therapeutic target for cancer. Tumor Biol. 2017, 39. [Google Scholar] [CrossRef] [Green Version]
  71. Yeh, M.; Oh, C.S.; Yoo, J.Y.; Kaur, B.; Lee, T.J. Pivotal role of microRNA-138 in human cancers. Am. J. Cancer Res 2019, 9, 1118–1126. [Google Scholar] [PubMed]
  72. Song, N.; Li, P.; Song, P.; Li, Y.; Zhou, S.; Su, Q.; Li, X.; Yu, Y.; Li, P.; Feng, M.; et al. MicroRNA-138-5p Suppresses Non-small Cell Lung Cancer Cells by Targeting PD-L1/PD-1 to Regulate Tumor Microenvironment. Front. Cell Dev. Biol. 2020, 8, 540. [Google Scholar] [CrossRef] [PubMed]
  73. Pinho, J.D.; Silva, G.E.B.; Júnior, A.A.L.T.; Belfort, M.R.D.C.; Macedo, J.M.; Da Cunha, I.W.; Quintana, L.G.; Calixto, J.D.R.R.; Nogueira, L.R.; Coelho, R.W.P.; et al. MIR-107, MIR-223-3P and MIR-21-5P Reveals Potential Biomarkers in Penile Cancer. Asian Pac. J. Cancer Prev. 2020, 21, 391–397. [Google Scholar] [CrossRef] [PubMed]
  74. Hartz, J.M.; Engelmann, D.; Fürst, K.; Marquardt, S.; Spitschak, A.; Goody, D.; Protzel, C.; Hakenberg, O.W.; Pützer, B.M. Integrated Loss of miR-1/miR-101/miR-204 Discriminates Metastatic from Nonmetastatic Penile Carcinomas and Can Predict Patient Outcome. J. Urol. 2016, 196, 570–578. [Google Scholar] [CrossRef] [PubMed]
Cancers 13 01480 g001 550
Figure 1. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in HPV-positive versus HPV-negative (a) PSCC tumor tissues (p ≤ 0.01; q ≤ 0.027; fold change ≥ 2) and (b) normal tissues (p ≤ 0.05; q ≤ 0.37; fold change ≥ 2).
Figure 1. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in HPV-positive versus HPV-negative (a) PSCC tumor tissues (p ≤ 0.01; q ≤ 0.027; fold change ≥ 2) and (b) normal tissues (p ≤ 0.05; q ≤ 0.37; fold change ≥ 2).
Cancers 13 01480 g001
Cancers 13 01480 g002 550
Figure 2. Multiple group comparison of differentially expressed miRNAs in basaloid, warty–basaloid, warty, and usual histological subtypes using the F-test (q ≤ 0.01). (a) Heatmap; (b) principal component analysis plot.
Figure 2. Multiple group comparison of differentially expressed miRNAs in basaloid, warty–basaloid, warty, and usual histological subtypes using the F-test (q ≤ 0.01). (a) Heatmap; (b) principal component analysis plot.
Cancers 13 01480 g002
Cancers 13 01480 g003 550
Figure 3. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in usual (HPV-negative) tumors versus related normal tissue samples (p ≤ 0.05; q ≤ 0.23; fold change ≥ 2).
Figure 3. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in usual (HPV-negative) tumors versus related normal tissue samples (p ≤ 0.05; q ≤ 0.23; fold change ≥ 2).
Cancers 13 01480 g003
Cancers 13 01480 g004 550
Figure 4. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in metastatic versus non-metastatic (HPV-negative) usual PSCC (p ≤ 0.01; q ≤ 0.08; fold change ≥ 2).
Figure 4. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in metastatic versus non-metastatic (HPV-negative) usual PSCC (p ≤ 0.01; q ≤ 0.08; fold change ≥ 2).
Cancers 13 01480 g004
Cancers 13 01480 g005 550
Figure 5. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in basaloid, warty, and warty–basaloid (HPV-positive) PSCC versus related normal tissue samples (p ≤ 0.05; q ≤ 0.17; fold change ≥ 2).
Figure 5. Heat map after unsupervised hierarchical clustering of differentially expressed miRNAs in basaloid, warty, and warty–basaloid (HPV-positive) PSCC versus related normal tissue samples (p ≤ 0.05; q ≤ 0.17; fold change ≥ 2).
Cancers 13 01480 g005
Cancers 13 01480 g006 550
Figure 6. Validation of differentially expressed miRNAs in different PSCC histological subtypes as determined by TaqMan® RT-qPCR (p ≤ 0.05). (a) Mean ΔCt values of (HPV-negative) usual PSCC and normal tissue samples; (b,c) mean ΔCt values of metastatic versus non-metastatic usual (HPV-negative) tumors; (d) mean ΔCt values of warty–basaloid (HPV-positive) PSCC versus normal tissue samples. Lower ΔCt implies higher expression level. Mann–Whitney U test was used to determine p-values, FC: Fold change.
Figure 6. Validation of differentially expressed miRNAs in different PSCC histological subtypes as determined by TaqMan® RT-qPCR (p ≤ 0.05). (a) Mean ΔCt values of (HPV-negative) usual PSCC and normal tissue samples; (b,c) mean ΔCt values of metastatic versus non-metastatic usual (HPV-negative) tumors; (d) mean ΔCt values of warty–basaloid (HPV-positive) PSCC versus normal tissue samples. Lower ΔCt implies higher expression level. Mann–Whitney U test was used to determine p-values, FC: Fold change.
Cancers 13 01480 g006
Table
Table 1. Summary of the clinical and histopathological data of the penile squamous cell carcinoma (PSCC) samples. n/a: Not available. In addition to lymph node metastases, visceral metastases were found in liver and lung. HPV: Human papillomavirus.
Table 1. Summary of the clinical and histopathological data of the penile squamous cell carcinoma (PSCC) samples. n/a: Not available. In addition to lymph node metastases, visceral metastases were found in liver and lung. HPV: Human papillomavirus.
PATIENT CHARACTERISTICS N%
HISTOLOGICAL SUBTYPEusual1851.4
warty–basaloid617.1
basaloid822.9
warty12.9
clear cell12.9
verrucous12.9
PRIMARY TUMORpT1a822.9
pT1b411.4
pT21131.4
pT31028.6
n/a25.7
REGIONAL LYMPH NODESN01748.6
N1617.1
N212.9
N3925.7
n/a25.7
HISTOLOGIC GRADEG125.7
G21645.7
G31542.9
n/a25.7
METASTASIS STATUSNon-metastatic1851.4
Metastatic1748.6
HPV STATUSHPV+1748.6
HPV-1851.4
Table
Table 2. Summary of miRNAs validated by microarray and RT-qPCR in different histological subtypes.
Table 2. Summary of miRNAs validated by microarray and RT-qPCR in different histological subtypes.
miRNAsFold Change in Microarray,
p-Value		Fold Change in RT-qPCR,
p-Value		HPV StatusComparisonHistological Subtype
miR-99a-5p−8.0
(p = 0.037)		−3.5
(p = 0.017)		HPV-negativeTumor vs. NormalUsual
miR-137−3.7
(p = 0.004)		−8.54
(p = 0.004)		HPV-negativeMetastatic vs. Non-metastatic
miR-328-3p−2.7
(p = 0.007)		−1.98
(p = 0.032)		HPV-negativeMetastatic vs. Non-metastatic
miR-181d-5p+3.4
(p = 0.004)		+2.01
(p = 0.002)		HPV-positiveTumor vs. NormalBasaloid, Warty, and Warty–basaloid
miR-211-5p−16.9
(p = 0.046)		Not definedHPV-positiveTumor vs. Normal
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ayoubian, H.; Heinzelmann, J.; Hölters, S.; Khalmurzaev, O.; Pryalukhin, A.; Loertzer, P.; Heinzelbecker, J.; Lohse, S.; Geppert, C.; Loertzer, H.; et al. miRNA Expression Characterizes Histological Subtypes and Metastasis in Penile Squamous Cell Carcinoma. Cancers 2021, 13, 1480. https://doi.org/10.3390/cancers13061480

AMA Style

Ayoubian H, Heinzelmann J, Hölters S, Khalmurzaev O, Pryalukhin A, Loertzer P, Heinzelbecker J, Lohse S, Geppert C, Loertzer H, et al. miRNA Expression Characterizes Histological Subtypes and Metastasis in Penile Squamous Cell Carcinoma. Cancers. 2021; 13(6):1480. https://doi.org/10.3390/cancers13061480

Chicago/Turabian Style

Ayoubian, Hiresh, Joana Heinzelmann, Sebastian Hölters, Oybek Khalmurzaev, Alexey Pryalukhin, Philine Loertzer, Julia Heinzelbecker, Stefan Lohse, Carol Geppert, Hagen Loertzer, and et al. 2021. "miRNA Expression Characterizes Histological Subtypes and Metastasis in Penile Squamous Cell Carcinoma" Cancers 13, no. 6: 1480. https://doi.org/10.3390/cancers13061480

APA Style

Ayoubian, H., Heinzelmann, J., Hölters, S., Khalmurzaev, O., Pryalukhin, A., Loertzer, P., Heinzelbecker, J., Lohse, S., Geppert, C., Loertzer, H., Wunderlich, H., Bohle, R. M., Stöckle, M., Matveev, V. B., Hartmann, A., & Junker, K. (2021). miRNA Expression Characterizes Histological Subtypes and Metastasis in Penile Squamous Cell Carcinoma. Cancers, 13(6), 1480. https://doi.org/10.3390/cancers13061480

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop