Next Article in Journal
Regulation of Eosinophilia in Asthma—New Therapeutic Approaches for Asthma Treatment
Next Article in Special Issue
Use of Cytokine Mix-, Imiquimod-, and Serum-Induced Monoculture and Lipopolysaccharide- and Interferon Gamma-Treated Co-Culture to Establish In Vitro Psoriasis-like Inflammation Models
Previous Article in Journal
Primary Human Cardiomyocytes and Cardiofibroblasts Treated with Sera from Myocarditis Patients Exhibit an Increased Iron Demand and Complex Changes in the Gene Expression
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 10
Issue 4
10.3390/cells10040814
cells-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

Penile Cancer-Derived Cells Molecularly Characterized as Models to Guide Targeted Therapies

by
Hellen Kuasne
1,2,†,
Luisa Matos do Canto
1,†,
Mads Malik Aagaard
1,
Juan Jose Moyano Muñoz
2,
Camille De Jamblinne
3,
Fabio Albuquerque Marchi
2,
Cristovam Scapulatempo-Neto
4,5,
Eliney Ferreira Faria
4,6,
Ademar Lopes
7,
Sébastien Carréno
3 and
Silvia Regina Rogatto
1,*
1
Department of Clinical Genetics, University Hospital of Southern Denmark-Vejle, Institute of Regional Health Research, University of Southern Denmark, 7100 Vejle, Denmark
2
International Research Center—A.C.Camargo Cancer Center, São Paulo 01508-010, Brazil
3
Institute for Research in Immunology and Cancer and Département de Pathologie et de Biologie Cellulaire, Université de Montréal, Montréal, QC H3C 3J7, Canada
4
Molecular Oncology Research Center, Barretos Cancer Hospital, Barretos 14784-400, Brazil
5
Diagnósticos da América—DASA, Barueri, São Paulo 06455-010, Brazil
6
Uro-oncology and Robotic Surgery, Hospital Felicio Rocho, Belo Horizonte 30110-934, Brazil
7
Pelvic Surgery Department, A.C.Camargo Cancer Center, São Paulo 01508-010, Brazil
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this paper.
Cells 2021, 10(4), 814; https://doi.org/10.3390/cells10040814
Submission received: 25 February 2021 / Revised: 30 March 2021 / Accepted: 2 April 2021 / Published: 6 April 2021
(This article belongs to the Special Issue Cell-Based Models of Diseases for Drug Discovery)
Browse Figures
Versions Notes

Abstract

:
Penile cancer (PeCa) is a common disease in poor and developing countries, showing high morbidity rates. Despite the recent progress in understanding the molecular events involved in PeCa, the lack of well-characterized in vitro models precludes new advances in anticancer drug development. Here we describe the establishment of five human primary penile cancer-derived cell cultures, including two epithelial and three cancer-associated fibroblast (CAF) cells. Using high-throughput genomic approaches, we found that the epithelial PeCa derived- cells recapitulate the molecular alterations of their primary tumors and present the same deregulated signaling pathways. The differentially expressed genes and proteins identified are components of key oncogenic pathways, including EGFR and PI3K/AKT/mTOR. We showed that epithelial PeCa derived cells presented a good response to cisplatin, a common therapeutic approach used in PeCa patients. The growth of a PeCa-derived cell overexpressing EGFR was inhibited by EGFR inhibitors (cetuximab, gefitinib, and erlotinib). We also identified CAF signature markers in three PeCa-derived cells with fibroblast-like morphology, indicating that those cells are suitable models for PeCa microenvironment studies. We thus demonstrate the utility of PeCa cell models to dissect mechanisms that promote penile carcinogenesis, which are useful models to evaluate therapeutic approaches for the disease.

1. Introduction

Penile cancer (PeCa) is an aggressive and mutilating disease that presents a high incidence in poor and developing countries [1,2]. According to two population-based cancer registry surveys, the survival of PeCa patients has not improved in Europe or the United States in the last three decades [3,4]. Although ongoing preclinical studies show promising results for more personalized therapy, the lack of improvement in survival is probably due to delayed diagnosis and lack of advances in curative standardized treatment options (reviewed in [5]).
Genetic and epigenetic alterations associated with PeCa development and progression have pointed out potential therapeutic targets (reviewed in [6,7]). In the last few years, we have reported molecular markers and deregulated pathways that are potentially helpful for discovering such targets in PeCa [8,9,10,11]. However, the scarcity of in vitro or animal models emerges as an obstacle that hampers the establishment of new and efficient therapeutic strategies for PeCa [12].
Cell cultures are valuable models for evaluating the molecular mechanisms underlying tumor initiation, progression, and drug response. Few cell lines derived from primary PeCa have been reported, and none of them are commercially available [13,14,15,16,17,18,19,20]. Two studies established and characterized cell lines derived from primary PeCa and corresponding lymph node metastasis [17,20]. We also reported a comprehensive characterization of a cell culture and xenograft derived from a verrucous tumor, which accounts for 2–8% of PeCa cases [10]. From 21 fresh PeCa samples, Chen et al. [18] established one cell line derived from a metastatic lymph node presenting a deleterious TP53 mutation. This cell line was sensitive to cisplatin (commonly used in advanced PeCa) and epirubicin. The authors suggested that epirubicin was an effective chemotherapeutic agent, though it is not commonly used in PeCa treatment. Zhou et al. [19] established a panel of five PeCa cell lines sensitive to cisplatin. They reported that cell lines presenting EGFR DNA amplification and/or protein overexpression were resistant to anti-EGFR therapies. The authors presumed that HRAS and PI3KCA mutations might be related to resistance to anti-EGFR therapy. These studies have demonstrated that tumor-derived cell lines are good models for testing therapeutic agents and investigating drug resistance mechanisms, drug discovery, and targeted treatments.
Cells in the tumor microenvironment also play a crucial role in cancer progression and sensitivity to therapy. Studies in urological cancer types have drawn attention to the influence of stromal-epithelial interactions on tumor growth, invasion, and immune response [21,22,23]. Among stromal cells, cancer-associated fibroblasts (CAFs) modulate cancer metastasis through several mechanisms [24].
In the present study, we established and characterized five penile cancer-derived cells from cancer tissues (two epithelial and three CAFs). We evaluated the morphology of these cells and their ability to proliferate, migrate, and invade. Different -omics approaches were applied to characterize these derived PeCa cells molecularly. We showed that the tumor epithelial cells retained the genetic features of primary tissues. The response to cisplatin and EGFR inhibitors was also investigated. Overall, our results showed that these newly established cells could be used in pre-clinical assays to investigate drug response in PeCa.

2. Materials and Methods

Specimens derived from primary PeCa obtained after surgery from patients naïve of treatment were cultured under sterile conditions. The study was conducted following the Declaration of Helsinki and approved by the Human Research Ethics Committees of the A.C.Camargo Cancer Center (Protocol 1230/2009) and Barretos Cancer Hospital (Protocol 363-2010), São Paulo, Brazil. All subjects provided informed consent.
The clinical and pathological characteristics of the patients are described in Table 1. Five PeCa derived cells (cell 2 to cell 6) were successfully established. Cells 2, 3, and 6 were derived from mixed usual-sarcomatoid, verrucous, and basaloid PeCa subtypes, respectively. Cells 4 and 5 were derived from the usual subtype. The PeCa from patients 2, 5, and 6 showed a high clinical stage (III and IV). Only the PeCa from patient 2 was positive for the human papillomavirus (HPV16). A xenograft model derived from cell 3 was previously published by our group [10]. Translatomic and reverse-phase protein arrays (RPPA) for each tumor-derived cell (cells 2, 3, 4, 5, and 6) were compared with the foreskin cell line obtained from a healthy individual (cell 1), kindly donated by Dr. Silvya Stuchi Maria-Engler, Clinical Chemistry and Toxicology Department, University of São Paulo, SP, BR.

2.1. Establishment of Penile Cancer-Derived Cells

Histopathological analyses confirmed the presence of tumor cells in the samples used for culture. Five PeCa samples were dissociated and plated as previously described [10]. Briefly, minced tumor fragments were seeded in 25 cm2 culture flasks containing 3:1 KSFM (keratinocyte serum-free medium)–DMEM/F12 (Dulbecco’s modified Eagle medium/nutrient mixture F-12) (GIBCO, Carlsbad, CA, USA) added with 2.5% fetal bovine serum (FBS) (HYCLONE, Waltham, MA, USA), 30 µg/mL of bovine pituitary extract (BPE) (GIBCO, Carlsbad, CA, USA), 0.2 ng/mL of epithelial growth factor (EGF) (GIBCO, Carlsbad, CA, USA), antibiotics (100 IU/mL penicillin G and 100 mg/mL streptomycin) (Sigma-Aldrich, St. Louis, MO, USA), and 25 µg/mL of Fungizone® Antimycotic (GIBCO, Carlsbad, CA, USA). After confluence, cells were treated with 0.05% trypsin/0.02% EDTA (Sigma-Aldrich, St. Louis, MO, USA) and replicated for at least 10 passages (P10).

2.2. Morphological Characterization and Functional Assays

The morphology of PeCa cells was evaluated by phase-contrast microscopy (Nikon TE2000, Amsterdam, Netherlands) and immunofluorescence by using Texas Red: actin/phalloidin (Thermo Fisher Scientific, Waltham, MA, USA), FITC (fluorescein isothiocyanate): tubulin (Thermo Fisher Scientific, Waltham, MA, USA), and DAPI (4′,6-diamidino-2-phenylindole): nucleus (Vector Laboratories, Burlingame, CA, USA).
The doubling time was calculated by seeding the cells at a density of 1.5 × 105 cells in a 6 cm2 plate. Cell count was performed using Trypan blue (GIBCO, Carlsbad, CA, USA) dye exclusion in a Neubauer chamber every 24 h for eight days (the medium was replaced every three days). The doubling time was determined from the growth curves using the data from three independent experiments, each with three technical replicates.
Migration and invasion assays were conducted in transwell (Costar, Corning, New York, NY, USA) and Matrigel™-coated transwell plates (BD Biosciences, San Jose, CA, USA), respectively, as previously described [9]. Briefly, 5 × 105 cells suspended in 100 µL of serum-free medium were added to the upper chamber of a transwell plate. The lower chamber contained 2.5% FBS (HYCLONE, Waltham, MA, USA), 30 µg/mL BPE (GIBCO, Carlsbad, CA, USA), and 0.2 ng/mL EGF (GIBCO, Carlsbad, CA, USA) supplemented medium, which acted as a chemoattractant, for 24 h at 37 °C. Cells incubated with serum-free medium only were used as controls. After fixation, the cells that migrated across the membrane were stained with hematoxylin-eosin. The number of cells that migrated and/or invaded were counted and scored from 0 to 4 (score 0: 0–5%, score 1: 5–25%, score 2: 25–50%, score 3: 50–75%, and score 4: 75–100%).

2.3. Molecular Profiling of Penile Cancer-Derived Cells

Genomic DNA and mRNA were isolated from all PeCa cells at passage 10. DNA was extracted using DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA, USA). The human HPV genotyping was performed using the Linear Array HPV Genotyping Test (Roche Molecular Diagnostics, Branchburg, NJ, USA).
Copy number alterations (CNAs) and copy-neutral loss of heterozygosity (cnLOH) were assessed using the CytoScan HD platform (Thermo Fisher Scientific, Waltham, MA, USA). Data were analyzed with the Chromosome Analysis Suite 3.0 (ChAS) software (Thermo Fisher Scientific, Waltham, MA, USA, v.4.0), considering at least 25 markers for losses/mosaic losses, 50 markers for gains/mosaic gains, and cnLOHs with a minimum of 5 Mb, as previously described [25,26]. A distogram showing dissimilarity between samples (matched primary tumor and cell culture) was performed using the R package Rawcopy [27].
The mutational profiling of primary tumors and their derived cells was investigated by targeted next-generation sequencing (tNGS) of 105 cancer-related genes (SureSelectXT Custom Panel, Agilent, Santa Clara, CA, USA). The DNA libraries of all coding exons, intron–exon boundaries, and 3′ and 5′ UTR of the genes were prepared using SureSelectQXT Library Prep Kit (Agilent, Santa Clara, CA, USA). The DNA sequencing was carried out using 2 × 150-bp paired-end technology on the Next-Seq500 Illumina platform (Illumina, San Diego, CA, USA). All sequenced data were aligned to UCSC hg19. The initial analysis was performed as previously described [26]. Variants were filtered out based on their classification (benign—B or likely benign—LB, according to the ACMG—American College of Medical Genetics and Genomics or ClinVar), and frequency (>0.1 in the GnomAD v2.1.1—the Genome Aggregation Database [28], and/or in the ABRaOM—Online Archive of Brazilian Mutations [29] datasets). We excluded variants in homopolymer regions.
Actively translating ribosomes and their associated mRNAs were isolated, as described by Roffé et al. [30]. mRNA was extracted from polysome fractions using TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA), according to the manufacturer’s recommendations. RNA quality and quantity assessment were conducted using an RNA 6000 Nano labchip (Bioanalyzer, Agilent, Santa Clara, CA, USA) and Nanodrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), respectively. Translatomic profiling was investigated with the Clariom™ D Assay platform (Thermo Fisher Scientific, Waltham, MA, USA). Labeling, hybridization, and washing followed the manufacturer’s recommendations. Normalization, quality control, and analysis were assessed by the Transcriptome Analysis Console (TAC software Thermo Fisher Scientific, Waltham, MA, USA, v.4.0). The expression levels of PeCa-derived cells were compared to cell 1 (normal foreskin). Only probes with adjusted p-value < 0.001 and |fold change| > 4 were considered significant. Experiments were performed in duplicate. Expression levels of genes considered markers of CAFs or related to their activation were evaluated [31,32,33].
A high-throughput antibody-microarray-based technique (reverse-phase protein array—RPPA) was used to evaluate the protein expression profile in all PeCa-derived cells compared to cell 1. The RPPA experiments were performed in the RPPA core facility of MD Anderson Cancer Center, Houston, Texas, USA. Protein lysates were extracted according to the RPPA protocol (https://www.mdanderson.org/research/research-resources/core-facilities/functional-proteomics-rppa-core.html, accessed on 4 March 2021). Graphics representing RPPA results were build using the RPPApipe tool [34].
Gene set enrichment analysis was performed on the translatomic data using the single-sample gene set enrichment analysis (ssGSEA) tool (http://genepattern.broadinstitute.org, accessed on 6 April 2021). Annotated gene sets were obtained from the Reactome and oncogenic signature sub-collection of the Molecular Signature Database (MSigDB) (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp, accessed on 6 April 2021). Ingenuity Pathway Analysis (IPA v01.12) (QIAGEN Inc, Germantown, MD, USA) was used to predict the top transcriptional regulators that are either activated or inhibited in PeCa-derived cells (upstream regulator analysis).
The gene expression levels of one study on penile cancer (GSE57955) publicly available [35] were assessed to identify the EGFR expression signatures (Reactome signaling by EGFR and KEEG ERBB signaling pathways). Heatmaps were generated using the Morpheus web tool (https://software.broadinstitute.org/morpheus/. accessed on 6 April 2021).

3. Results

3.1. Penile Cancer Derived-Cells Present Diverse Morphology, Proliferation, Migratory, and Invasive Characteristics

Two (cells 2 and 3) out of five PeCa-derived cells presented epithelial morphology (Figure 1A,B) whereas cells 4, 5, and 6 presented fibroblast-like morphology. Cells 2 and 3 showed positive expression of cytokeratin (Figure S1A) and epiplakin (VHL-PPK1) (Figure S2). Cell 3 presented both cytokeratin and vimentin expression (Figure S1A). Although cell 2 was derived from an HPV16 case, the positivity was lost after five passages in vitro.
Cells presenting the fibroblast-like morphology (Figure 1A,B) showed positive expression of vimentin in both immunocytochemistry and RPPA experiments, among other markers related to CAFs (Figures S1A and S2).
The doubling time was calculated for all cell cultures (Figure S1B). Cells 2, 3, 4, and 6 presented an average doubling time of 23.1 h, whereas cell 5 had a doubling time of 27.7 h (Figure S1B).
The ability of all cell lines to migrate and invade was evaluated using transwell assays. Cell 4 presented the highest average score of migration (score 4), followed by cells 3, 6 (score 3), and 2 (score 2.5), with cell 5 showing the lowest score (score 0.5) (Figure 1C,D). The invasive potential was high in cells 3, 4, and 6 (2 to 4 scores) and low in cell 2 (0 to 1 score), whereas cell 5 did not invade (0 score) (Figure 1D and Figure S1C). The ability to migrate and invade was not related to the morphology of the different cell lines.

3.2. Penile Cancer-Derived Cells Recapitulate the Molecular Profile of PeCa Primary Tissues

Genomic copy number alterations found in primary tumors were compared with their PeCa-derived cells (Cytoscan HD platform, Thermo Fisher Scientific, Waltham, MA, USA), revealing a high similarity level. The cells presented at least 80% matching single nucleotide polymorphisms with the primary tumor, confirming their parental origin (Figure 2A). In this analysis, cells 3 and 6 presented 100% genomic similarities with their matched primary tumors. A lower similarity was observed between the tumor and its derived cell 5 (~80%). Supplementary Table S1 and Figure S3 show the genomic alterations identified in primary tumors and their derived cells.
Overall, the primary tumors and their derived cells presented a lower number of chromosomal imbalances, with the highest number of altered regions in cell 3 (44) and the lowest in cell 4 (8). Cells with epithelial morphology shared a higher number of common alterations with the matched primary tissues. Primary tumor 2 and its derived cell 2 presented 24 chromosomal imbalances, including a homozygous loss of 9p21.3 (CDKN2A and CDKN2B genes). Cell 3 and its primary tumor 3 presented 28 genomic alterations in common (Table S1), including gains of 4q12 (PDGFRA), an amplification of 11p15.5 (H19 and IGF2 genes), and a loss of 6p25.3 (DUSP22). A loss of 8p, encompassing the DLC1 gene, was identified only in cell 3.
Although CAFs shared common altered regions with their primary tumors, the alterations were detected in a lower frequency compared to epithelial cells. Three chromosomal imbalances (gains of 14q32.33 and losses of 11p15.4 and 5q23) were shared by tumor 4 and its derived cell 4. Primary tumor 5 and its derived cell 5 presented 13 CNAs, including gains of 3p12.1, 14q32.33, and 17q21.31 and losses of 5q23.1, 6p25.3, and 12p11.21. Cell 6 presented gains of 3q26.31 (TNFSF10), 14q32.33, 16p11.2, and 17q21.31, and losses of 5q35.3, 6p25.3, 8p11.22, 8p22, 12p13.33, and 19p12. Tumor 6 and its derived cell 6 presented the highest number of cnLOH (Figure S3).
We also detected similar genomic alterations in more than one derived tumor cell (Table S1). Among them, cells 2 and 3 shared a loss of 9p21.3 (CDKN2A and CDKN2B), all PeCa derived cells presented gains of 14q32.33, and cells 5 and 6 presented gains of 17q21.31 (KANSL1). A detailed description of the genomic alterations detected in our cells paired with their respective primary tumors is presented in Table S1 and Figure S3.
Targeted next-generation sequencing was performed in all PeCa-derived cells and three matching primary tumors (patients 2, 5, and 6) using a custom panel composed of 105 cancer-related genes. After filtering, we identified variants in all but one cell line (cell 3). A total of 18 single nucleotide variants - SNVs (nonsynonymous, stop gains, splice site, three and five prime untranslated regions) and two indels mapping to 17 genes were detected. The variants (missense, frameshift, and truncation) identified in the primary tumors and their matched derived cells are represented in Figure 2B and Table 2. All variants found in tumor 2 were also observed in cell 2 (two variants of PIK3CA and one of FGFR1). Tumor 5 and its derived cell 5 presented variants in nine different genes, among them ALK, BRCA2, NOTCH1, and PIK3CA. One variant in a CpG island, in the promoter region of the APC gene, was common to tumor 6 and its derived cell 6. Two variants were found exclusively in cell 2 and not in its tumor (MMP1 and STAT3), whereas two variants were found only in the tumor from which cell 6 was derived (ERBB2 and MLH3). The gene most frequently altered was PIK3CA (two variants in cell/tumor 2, one in cell/tumor 5). Table 2 summarizes the tNGS results found in the cells compared with the primary tumors.

3.3. Identification of Potential Therapeutic Targets in PeCa Cells Using Translatomic, Pathways, and Protein Analysis

To evaluate the actively translated mRNAs, we performed a polysome profiling analysis using a sucrose gradient, which involves the separation of mRNAs into fractions according to the number of bound ribosomes. The mRNA profile from the polysomal fraction of each cell was investigated using the Clariom™ D Assay platform (ThermoFisher, USA) (Figure 3A,B). A heatmap containing the top 250 differentially expressed genes (DEG) showed three clusters. The first was composed of cell 3 (epithelial morphology), the second cluster of cell 2 (epithelial morphology), and the third cluster of cells 4, 5, and 6 (fibroblast-like morphology—CAFs). Penile cancer-derived cells 4, 5, and 6 presented overexpression of markers related to a CAF signature (MMP2, REAB3B, COL6A1, COL6A2, CTSK, THY1, PDGFRA, DCN, and fibroblast activation ACTA2) [31,32] (Figure 3B). The fibroblast activation protein α gene (FAP) did not show increased expression. The top 15 DEG of cells 2 to 6 compared with cell 1 (normal foreskin) and the number of DEGs for each comparison are shown in Table 3.
A single-sample gene set enrichment analysis (ssGSEA) was performed, and the main pathways identified as dysregulated in the PeCa cells are depicted in Figure 3C,D. Using oncogenic signatures (Figure 3C) and the Reactome (Figure 3D) sub-collection of the MSigDB, we identified that EGFR, PI3K, and mTOR pathways were dysregulated in cell 2. Cell 2 also presented a high score for several cell junction signatures (Figure 3D). Cell 3 expressed genes related to the regulation of apoptosis (BCL2L1 and BCL2) and pathways related to SRC and ERK activation. Cells 4, 5, and 6 presented pathways related to the CAF phenotype, including YAP and MYC signatures (Figure 3C,D).
We also evaluated the protein profile of 304 pre-selected antibodies involved in pathways potentially dysregulated in cancer (Figure S2). Figure 3E shows the main altered proteins or phosphorylated isotypes in PeCa cells. Potential therapeutic targets were identified among the dysregulated proteins identified in the PeCa cells. Cell 2 presented increased expression of EGFR, HER3, and VEGFR2, which are described as targets for therapy and tested in clinical trials (NCT01728233—dacomitinib). Increased EGFR expression was observed in cell 2 (Figure 3F). In addition, several proteins downstream to the EGFR receptor were dysregulated (Figure 3E). Cell 2 also presented a decreased expression of p16INK4A, which correlates with the genomic deletion of the CDKN2A gene detected by CNA analysis (Figure S3A).
High levels of phosphorylated AKT were observed in cells 3 and 5. The epithelial cells marker, epiplakin (VHL-PPK1), was highly expressed in cells 2 and 3. A weak expression of epiplakin was observed in cell 4, whereas no expression of this marker was found in cells 5 or 6. Vimentin, a marker of mesenchymal cells, was expressed in cells 5 and 6 (Figure 3E and Figure S2).
The expression levels of PI3K/AKT/mTOR, EGFR, ERBB2, TP53, and CDKN2A proteins found in our cells were previously reported in PeCa (Table 4). In addition to EGFR and CDKN2A described above, we also detected alterations in targetable tyrosine kinase receptors, including ERBB2 and ERBB3, and several proteins downstream to this pathway, such as AKT, S6, and mTOR (Table 4, Figure 3E).

3.4. Identification of Potential Therapeutic Targets for PeCa and Chemo-Sensitivity Assays

We performed chemo-sensitivity assays in epithelial PeCa cells 2 and 3. Using the Ingenuity Pathway Analysis software (Qiagen, Valencia, CA, USA), we identified potential drugs that target the differentially expressed kinases and proteins that showed higher expression (Table 5). The upstream regulator analysis of cell 2 revealed several genes predicted to be regulated by EGFR (Figure 4A,B). EGFR was a core molecule in both mRNA (Figure 4A) and protein analysis in cell 2 (Figure 4B).
Based on the evidence that EGFR was dysregulated in cell 2 at mRNA and protein (RPPA and Western blot) levels (Figure 3E,F), we tested the chemo-sensitivity of this cell using anti-EGFR inhibitors (cetuximab, gefitinib, and erlotinib). Cell 3 was used as a negative control since it did not show EGFR overexpression. We also tested cell viability in response to cisplatin (commonly used in the treatment of advanced PeCa) in cells 2 and 3. The IC50 and dose-response curve results showed that these cell lines were sensitive to cisplatin. However, only cell 2 (EGFR overexpression) was sensitive to anti-EGFR drugs (Figure 4C).
Having shown that EGFR inhibition can block the proliferation of PeCa cells in vitro, we evaluated whether PeCa primary tissues overexpress the EGFR gene. Samples that presented overexpression of the EGFR signature most likely have activation of the EGFR pathway. We accessed EGFR mRNA-related signatures (MySigDB) in a set of 36 primary PeCa tissues previously published by our group (GSE57955) [35]. We identified a cluster of samples (~30% of the cases) showing EGFR mRNA signature overexpression (Figure 5A,B). Several genes from the EGFR signature (downstream to EGFR) were positively correlated with EGFR expression (Figure 5C).

4. Discussion

An in vitro culture of tumor cells is a valuable model for drug development and preclinical drug testing in oncology. At least in part, these cells maintained the molecular alterations of the parental tumor and thus can be used in functional studies [56,57]. However, establishing cancer cell lines from fresh tumor tissues represents a technical challenge [58]. These models also present certain limitations (such as the lack of tumor heterogeneity, the absence of components of the tumor microenvironment, and genotypic and phenotypic drift during culture) that must be taken into consideration [59].
The development of in vitro models of PeCa is hampered by the low incidence of this tumor type. In this study, we derived PeCa cells from fresh penile primary tumors. Two of the established PeCa-derived cells presented a typical polygonal epithelial cell morphology (cells 2 and 3), and three presented fibroblast-like morphology (cells 4, 5, and 6). These CAFs are components of the tumor microenvironment playing crucial roles in tumor progression and treatment response [60,61]. Interesting approaches were developed by Di Donato et al [23]., where the prostate tumor growth was stimulated by androgen in co-cultures using CAFs derived from patients with positive AR (androgen receptor) and cell lines that were either AR positive or negative [23]. Bladder cancer progression, migration, and invasion were also shown to be triggered by factors secreted by CAFs [22]. In both cases, inhibiting CAF signaling sufficed to decrease tumor aggressiveness [22,23]. Although several studies indicated the influence of the microenvironment in response to chemotherapy, we selected only cells showing epithelial features (cells 2 and 3) to evaluate the treatment response since the selected drugs target the tumor cells and not the microenvironment.
The primary tumor from patient 2 was positive for high-risk HPV, whereas the derived cell 2 was HPV negative. Attempts to reproduce HPV replication in standard cell culture have been unsuccessful, mostly because the replication process is linked to the differentiation of keratinocytes, and it is challenging to recreate the stratified structure of the epithelium in vitro [12,62].
The primary tumor that generated cell 2 was a mixed tumor (usual with few sarcomatoid differentiation areas). Sarcomatoid PeCas are rare and aggressive tumors (1–2% of PeCa) associated with metastasis and poor prognosis [63,64]. Using phase-contrast microscopy, we exclusively detected epithelial cells without spindle cells, expected in the sarcomatoid component (Figure 1A,B). Indeed, microarray and RPPA analysis showed the expression of several epithelial markers in cell 2, including KRT5, KRT17, CDH1, DSG3 (Table 3), and epiplakin (VHL-PPK1) (Figure S2). These results suggest that only the usual component of this mixed PeCa was selected in vitro.
Very few studies have evaluated the mRNA expression profile of PeCa cell lines [18,19]. To our knowledge, we were the first to evaluate the expression levels of ribosome-associated mRNAs (translatomic analysis) in PeCa-derived cells. Since these mRNAs are associated with several ribosomes, they are likely actively translated [65]. Cell 2 was the only PeCa-derived cell that presented high EGFR overexpression at gene and protein levels. Gene set enrichment analysis (Figure 3C,D) and upstream regulator analysis (Figure 4A,B) also showed that several genes downstream to the EGFR pathway and EGFR interactors were dysregulated. Silva Amancio et al. [42] reported increased EGFR protein expression (3+: 67 of 139 cases) and its association with cancer recurrence and perineural invasion. Chaux et al. [66] found EGFR overexpression in 44% of the PeCas evaluated. These studies suggested that targeting the EGFR pathway would be an effective therapeutic strategy to treat a subset of cases. In contrast, Zhou et al. [19] identified EGFR protein overexpression (Western blot) in four of five PeCa cell lines and resistance to erlotinib and afatinib. Only one of these cell lines showed a minor sensitivity to cetuximab. The authors suggested that HRAS and PI3KCA alterations might be related to anti-EGFR therapy resistance. A chemo-sensitivity assay with a panel of drugs targeting EGFR (cetuximab, gefitinib, and erlotinib) showed that cell 2 was sensitive to anti-EGFR drugs (Figure 4C).
Anti-EGFR agents have been used to treat PeCa patients, mainly as a salvage treatment after first-line chemotherapy failure [67,68,69]. Di Lorenzo et al. [67] reported that 50% of 28 patients (24 treated with cetuximab) were sensitive to anti-EGFR monoclonal antibodies [67]. PeCa patients treated with nimotuzumab after chemotherapy failure showed clinical response or stable disease [69]. Necchi et al. [70] reported that dacomitinib (pan-HER inhibitor) was active and well-tolerated in patients with advanced PeCa and may represent an option when combined chemotherapy cannot be administered. The evaluation of downstream effectors of EGFR signaling and how these mutations can affect therapy response over time should be taken into account. Cell 2 presented a good initial response against EGFR inhibitors, but the occurrence of the two PIK3CA pathogenic variants could be translated into the later acquisition of resistance to treatment, as described in metastatic colorectal cancer [71]. However, it was not possible to confirm this assumption because the patient was only treated by surgery and lost to follow-up. The effect of PIK3CA mutation and resistance to EGFR inhibitors in PeCa should be better investigated. Other receptors such as ERRB3, phospho-ERBB3, and VEGFR2 were also overexpressed in this cell line and represent potential therapeutic opportunities.
Several potentially targetable dysregulated oncogenic pathways have been described in PeCa [9,38,72,73,74], including PI3K/Akt/mTOR, EGFR, and ERRB2, among others. Recently, clinical trials testing immunotherapy in PeCa patients have shown promising results [75]. Positive expression of PD-L1 was found in 48–60% of cases [75], indicating that a significant proportion of these patients do not respond to immune-checkpoint inhibitors. Genetic studies have shown that the PI3K/Akt/mTOR pathway is significantly altered in PeCa cases [9,38,73,74]. Here, cells 2 and 5 presented alterations involving the PI3K/Akt/mTOR signaling pathway axis. We identified several mechanisms involved in the dysregulation of this pathway in PeCa cells, including PIK3CA mutation and increased expression levels of mTOR, RICTOR, phosphorylated S6 (cell 2), and phosphorylated AKT (cell 5), among others. The PI3K/AKT/mTOR signaling pathway has been associated with cell proliferation, migration, metabolism, and survival [76]. This pathway was dysregulated in cell 5, suggesting that the microenvironment cells can also be a potential target for cancer therapy [77].
Previously we described a similar genomic profile of cell 3 at passages 1, 5, and 10, thus showing genomic stability after cell culture [10]. Here, we used translatomic and protein expression analyses to better characterize this cell line. The enrichment analysis using DEGs revealed increased expression levels of Hedgehog pathway genes. Hedgehog signaling pathway inhibitors have been developed for cancer treatment [78]. Cell 3 also presented the dysregulation of other targetable pathways such as SRC and ERK (Figure 3C,D), widely investigated in several tumor types. This PeCa-derived cell presented a complex pattern of gene alteration that affected several oncogenic pathways. Thus, combined therapies targeting one or more dysregulated pathways are promising alternatives to treat PeCa, and clinical trials should be encouraged.
We showed that cell 2 and cell 3 were sensitive to cisplatin, but only cell 2 (EGFR overexpression) was sensitive to anti-EGFR drugs (Figure 4C). As expected, cell 3 did not respond to EGFR inhibitors. The EGFR status is a crucial step to be investigated before treating patients with EGFR inhibitors. The in-depth molecular characterization used in our PeCa-derived cells can assist the selection of drugs to be tested in representative models of PeCa.
The genomic analysis performed in five penile cancer-derived cells revealed a high level of similarities with the primary tumors. The most relevant differentially expressed genes and proteins are components of critical oncogenic pathways supporting previous studies published in the literature. We identified proteins with the potential to be targeted in PeCa patients, including EGFR, VEGFR2, PIK3/AKT/mTOR, MAPK, and SRC. Considering that our PeCa-derived cells closely recapitulated the molecular features of their primary tumors and presented the same dysregulated pathways, they are excellent models to be used in preclinical tests. Moreover, the cancer-associated fibroblasts can be used for microenvironment-related studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cells10040814/s1, Table S1: Genomic alterations shared in primary tumors and their derived PeCa cells, Figure S1: Characteristics of the penile cancer-derived cells; Figure S2: Heatmap of the reverse phase protein array (RPPA) showing the protein expression levels (in alphabetical order); Figure S3: Ideogram of primary tumors (thin lines in pink) and derived cell lines (thin lines in blue) showing the copy number alterations (CytoScan HD, Affymetrix).

Author Contributions

Conceptualization, H.K., S.C., and S.R.R.; methodology, H.K., L.M.d.C., J.J.M.M., C.D.J., M.M.A., and F.A.M.; software, M.M.A. and F.A.M.; formal analysis, H.K., L.M.d.C., C.S.-N., and F.A.M.; resources, C.S.-N., E.F.F., A.L., S.C., and S.R.R.; writing—original draft preparation, H.K. and S.R.R.; writing—review and editing, all authors; supervision, S.R.R.; funding acquisition, S.R.R., and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received financial support from the National Institute of Science and Technology in Oncogenomics (São Paulo Research Foundation–FAPESP: #2008/57887-9 and the National Council for Scientific and Technological Development–CNPq: #573589/08-9) and the Research Council of Lillebaelt Hospital, Denmark. H.K. received scholarship awards from FAPESP (#2013/03667-6 and #2015/25373-0).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Human Research Ethics Committees of the A.C.Camargo Cancer Center (Protocol 1230/2009) and Barretos Cancer Hospital (Protocol 363-2010), São Paulo, Brazil.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to acknowledge Barretos Cancer Hospital and A.C.Camargo Cancer Center, SP, Brazil, for providing human specimens. Our special thanks to the patients who agreed to participate in this study and Silvya Stuchi Maria-Engler (Clinical Chemistry and Toxicology Department, University of São Paulo, SP, Brazil) for the donation of the foreskin cell line. We are grateful to Annabeth H Peterson, Glaucia Hajj, and Carlos Eduardo Fonseca Alves for their support during the development of our study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Christodoulidou, M.; Sahdev, V.; Houssein, S.; Muneer, A. Epidemiology of penile cancer. Curr. Probl. Cancer 2015, 39, 126–136. [Google Scholar] [CrossRef] [PubMed]
  2. Douglawi, A.; Masterson, T.A. Updates on the epidemiology and risk factors for penile cancer. Transl. Androl. Urol. 2017, 6, 785–790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Verhoeven, R.; Janssen-Heijnen, M.; Saum, K.; Zanetti, R.; Caldarella, A.; Holleczek, B.; Brewster, D.; Hakulinen, T.; Horenblas, S.; Brenner, H.; et al. Population-based survival of penile cancer patients in Europe and the United States of America: No improvement since 1990. Eur. J. Cancer 2013, 49, 1414–1421. [Google Scholar] [CrossRef] [PubMed]
  4. Trama, A.; Foschi, R.; Larrañaga, N.; Sant, M.; Fuentes-Raspall, R.; Serraino, D.; Tavilla, A.; Van Eycken, L.; Nicolai, N.; Hackl, M.; et al. Survival of male genital cancers (prostate, testis and penis) in Europe 1999–2007: Results from the EUROCARE-5 study. Eur. J. Cancer 2015, 51, 2206–2216. [Google Scholar] [CrossRef] [PubMed]
  5. Thomas, A.; Necchi, A.; Muneer, A.; Tobias-Machado, M.; Tran, A.T.H.; Van Rompuy, A.-S.; Spiess, P.E.; Albersen, M. Penile cancer. Nat. Rev. Dis. Prim. 2021, 7, 1–24. [Google Scholar] [CrossRef]
  6. Kuasne, H.; Marchi, F.A.; Rogatto, S.R.; Cólus, I.M.D.S. Epigenetic Mechanisms in Penile Carcinoma. Int. J. Mol. Sci. 2013, 14, 10791–10808. [Google Scholar] [CrossRef] [Green Version]
  7. Thomas, A.; Vanthoor, J.; Vos, G.; Tsaur, I.; Albersen, M. Risk factors and molecular characterization of penile cancer. Curr. Opin. Urol. 2020, 30, 202–207. [Google Scholar] [CrossRef]
  8. Busso-Lopes, A.F.; Marchi, F.A.; Kuasne, H.; Scapulatempo-Neto, C.; Trindade-Filho, J.C.S.; De Jesus, C.M.N.; Lopes, A.; Guimarães, G.C.; Rogatto, S.R. Genomic Profiling of Human Penile Carcinoma Predicts Worse Prognosis and Survival. Cancer Prev. Res. 2015, 8, 149–156. [Google Scholar] [CrossRef] [Green Version]
  9. 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] [Green Version]
  10. Muñoz, J.J.A.M.; Drigo, S.A.; Kuasne, H.; Villacis, R.A.R.; Marchi, F.A.; Domingues, M.A.C.; Lopes, A.; Santos, T.G.; Rogatto, S.R. A comprehensive characterization of cell cultures and xenografts derived from a human verrucous penile carcinoma. Tumor Biol. 2016, 37, 11375–11384. [Google Scholar] [CrossRef]
  11. Muñoz, J.J.; Drigo, S.A.; Barros-Filho, M.C.; Marchi, F.A.; Scapulatempo-Neto, C.; Pessoa, G.S.; Guimarães, G.C.; Filho, J.C.S.T.; Lopes, A.; Arruda, M.A.; et al. Down-Regulation of SLC8A1 as a Putative Apoptosis Evasion Mechanism by Modulation of Calcium Levels in Penile Carcinoma. J. Urol. 2015, 194, 245–251. [Google Scholar] [CrossRef]
  12. Medeiros-Fonseca, B.; Cubilla, A.; Brito, H.; Martins, T.; Medeiros, R.; Oliveira, P.; Gil da Costa, R. Experimental Models for Studying HPV-Positive and HPV-Negative Penile Cancer: New Tools for An Old Disease. Cancers 2021, 13, 460. [Google Scholar] [CrossRef] [PubMed]
  13. Yamane, I.; Tsuda, T. Establishment of a Cell Line in vitro from the Lesion of a Clinical Case of Penis Cancroid. Tohoku J. Exp. Med. 1966, 88, 9–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ishikawa, S.; Kanoh, S.; Nemoto, S. Establishment of a cell line (TSUS-1) derived from a human squamous cell carcinoma of the penis. Hinyokika kiyo. Acta Urol. Jpn. 1983, 29, 373–376. [Google Scholar]
  15. Gentile, G.; Giraldo, G.; Stabile, M.; Beth-Giraldo, E.; Lonardo, F.; Kyalwazi, S.K.; Perone, L.; Ventruto, V. Cytogenetic study of a cell line of human penile cancer. Ann. de Génétique 1987, 30, 164–169. [Google Scholar]
  16. Tsukamoto, T. Establishment and characterization of a cell line (KU-8) from squamous cell carcinoma of the penis. Keio J. Med. 1989, 38, 277–293. [Google Scholar] [CrossRef] [PubMed]
  17. Naumann, C.M.; Sperveslage, J.; Hamann, M.F.; Leuschner, I.; Weder, L.; Al-Najar, A.A.; Lemke, J.; Sipos, B.; Jünemann, K.-P.; Kalthoff, H. Establishment and Characterization of Primary Cell Lines of Squamous Cell Carcinoma of the Penis and its Metastasis. J. Urol. 2012, 187, 2236–2242. [Google Scholar] [CrossRef]
  18. Chen, J.; Yao, K.; Li, Z.; Deng, C.; Wang, L.; Yu, X.; Liang, P.; Xie, Q.; Chen, P.; Qin, Z.; et al. Establishment and characterization of a penile cancer cell line, penl1, with a deleterious TP53 mutation as a paradigm of HPV-negative penile carcinogenesis. Oncotarget 2016, 7, 51687–51698. [Google Scholar] [CrossRef] [Green Version]
  19. Zhou, Q.-H.; Deng, C.-Z.; Li, Z.-S.; Chen, J.-P.; Yao, K.; Huang, K.-B.; Liu, T.-Y.; Liu, Z.-W.; Qin, Z.-K.; Zhou, F.-J.; et al. Molecular characterization and integrative genomic analysis of a panel of newly established penile cancer cell lines. Cell Death Dis. 2018, 9, 1–14. [Google Scholar] [CrossRef] [Green Version]
  20. Fenner, F.; Goody, D.; Protzel, C.; Erbersdobler, A.; Richter, C.; Hartz, J.M.; Naumann, C.M.; Kalthoff, H.; Herchenröder, O.; Hakenberg, O.W.; et al. E2F1 Signalling is Predictive of Chemoresistance and Lymphogenic Metastasis in Penile Cancer: A Pilot Functional Study Reveals New Prognostic Biomarkers. Eur. Urol. Focus 2018, 4, 599–607. [Google Scholar] [CrossRef]
  21. Wang, L.; Saci, A.; Szabo, P.M.; Chasalow, S.D.; Castillo-Martin, M.; Domingo-Domenech, J.; Siefker-Radtke, A.; Sharma, P.; Sfakianos, J.P.; Gong, Y.; et al. EMT- and stroma-related gene expression and resistance to PD-1 blockade in urothelial cancer. Nat. Commun. 2018, 9, 3503. [Google Scholar] [CrossRef]
  22. Goulet, C.R.; Champagne, A.; Bernard, G.; Vandal, D.; Chabaud, S.; Pouliot, F.; Bolduc, S. Cancer-associated fibroblasts induce epithelial–mesenchymal transition of bladder cancer cells through paracrine IL-6 signalling. BMC Cancer 2019, 19, 1–13. [Google Scholar] [CrossRef] [Green Version]
  23. Di Donato, M.; Zamagni, A.; Galasso, G.; Di Zazzo, E.; Giovannelli, P.; Barone, M.V.; Zanoni, M.; Gunelli, R.; Costantini, M.; Auricchio, F.; et al. The androgen receptor/filamin a complex as a target in prostate cancer microenvironment. Cell Death Dis. 2021, 12, 1–17. [Google Scholar] [CrossRef]
  24. Sahai, E.; Astsaturov, I.; Cukierman, E.; DeNardo, D.G.; Egeblad, M.; Evans, R.M.; Fearon, D.; Greten, F.R.; Hingorani, S.R.; Hunter, T.; et al. A framework for advancing our understanding of cancer-associated fibroblasts. Nat. Rev. Cancer 2020, 20, 174–186. [Google Scholar] [CrossRef] [Green Version]
  25. Villacis, R.A.R.; Basso, T.R.; Canto, L.M.; Nóbrega, A.F.; Achatz, M.I.; Rogatto, S.R. Germline large genomic alterations on 7q in patients with multiple primary cancers. Sci. Rep. 2017, 7, srep41677. [Google Scholar] [CrossRef]
  26. Canto, L.M.D.; Larsen, S.J.; Kupper, B.E.C.; Begnami, M.D.F.D.S.; Scapulatempo-Neto, C.; Petersen, A.H.; Aagaard, M.M.; Baumbach, J.; Aguiar, S.J.; Rogatto, S.R. Increased Levels of Genomic Instability and Mutations in Homologous Recombination Genes in Locally Advanced Rectal Carcinomas. Front. Oncol. 2019, 9, 395. [Google Scholar] [CrossRef] [PubMed]
  27. Mayrhofer, M.; Viklund, B.; Isaksson, A. Rawcopy: Improved copy number analysis with Affymetrix arrays. Sci. Rep. 2016, 6, 36158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Karczewski, K.J.; Genome Aggregation Database Consortium; Francioli, L.C.; Tiao, G.; Cummings, B.B.; Alfoeldi, J.; Wang, Q.S.; Collins, R.L.; Laricchia, K.M.; Ganna, A.; et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020, 581, 434–443. [Google Scholar] [CrossRef]
  29. Naslavsky, M.S.; Yamamoto, G.L.; De Almeida, T.F.; Ezquina, S.A.M.; Sunaga, D.Y.; Pho, N.; Bozoklian, D.; Sandberg, T.O.M.; Brito, L.A.; Lazar, M.; et al. Exomic variants of an elderly cohort of Brazilians in the ABraOM database. Hum. Mutat. 2017, 38, 751–763. [Google Scholar] [CrossRef] [PubMed]
  30. Roffé, M.; Hajj, G.N.; Azevedo, H.F.; Alves, V.S.; Castilho, B.A. IMPACT Is a Developmentally Regulated Protein in Neurons That Opposes the Eukaryotic Initiation Factor 2α Kinase GCN2 in the modulation of Neurite Outgrowth*. J. Biol. Chem. 2013, 288, 10860–10869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Liu, T.; Han, C.; Wang, S.; Fang, P.; Ma, Z.; Xu, L.; Yin, R. Cancer-associated fibroblasts: An emerging target of anti-cancer immunotherapy. J. Hematol. Oncol. 2019, 12, 1–15. [Google Scholar] [CrossRef]
  32. Kalluri, R.; Zeisberg, M. Fibroblasts in cancer. Nat. Rev. Cancer 2006, 6, 392–401. [Google Scholar] [CrossRef]
  33. Nurmik, M.; Ullmann, P.; Rodriguez, F.; Haan, S.; Letellier, E. In search of definitions: Cancer-associated fibroblasts and their markers. Int. J. Cancer 2020, 146, 895–905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Eichner, J.; Heubach, Y.; Ruff, M.; Kohlhof, H.; Strobl, S.; Mayer, B.; Pawlak, M.; Templin, M.F.; Zell, A. RPPApipe: A pipeline for the analysis of reverse-phase protein array data. Biosystems 2014, 122, 19–24. [Google Scholar] [CrossRef] [PubMed]
  35. Kuasne, H.; Cólus, I.M.D.S.; Busso, A.F.; Hernandez-Vargas, H.; Barros-Filho, M.C.; Marchi, F.A.; Scapulatempo-Neto, C.; Faria, E.F.; Lopes, A.; Guimarães, G.C.; et al. Genome-wide methylation and transcriptome analysis in penile carcinoma: Uncovering new molecular markers. Clin. Epigenetics 2015, 7, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Stankiewicz, E.; Prowse, D.M.; Ng, M.; Cuzick, J.; Mesher, D.; Hiscock, F.; Lu, Y.-J.; Watkin, N.; Corbishley, C.; Lam, W.; et al. Alternative HER/PTEN/Akt Pathway Activation in HPV Positive and Negative Penile Carcinomas. PLoS ONE 2011, 6, e17517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Chaux, A.; Munari, E.; Cubilla, A.L.; Hicks, J.; Lecksell, K.; Burnett, A.L.; Netto, G.J. Immunohistochemical expression of the mammalian target of rapamycin pathway in penile squamous cell carcinomas: A tissue microarray study of 112 cases. Histopathology 2014, 64, 863–871. [Google Scholar] [CrossRef] [PubMed]
  38. Azizi, M.; Tang, D.H.; Verduzco, D.; Peyton, C.C.; Chipollini, J.; Yuan, Z.; Schaible, B.J.; Zhou, J.-M.; Johnstone, P.A.; Giuliano, A.; et al. Impact of PI3K-AKT-mTOR Signaling Pathway Up-regulation on Prognosis of Penile Squamous-Cell Carcinoma: Results From a Tissue Microarray Study and Review of the Literature. Clin. Genitourin. Cancer 2019, 17, e80–e91. [Google Scholar] [CrossRef] [PubMed]
  39. Faraj, S.F.; Chaux, A.; Gonzalez-Roibon, N.; Munari, E.; Cubilla, A.L.; Shih, I.-M.; Netto, G.J. Immunohistochemical expression of ARID1A in penile squamous cell carcinomas: A tissue microarray study of 112 cases. Hum. Pathol. 2015, 46, 761–766. [Google Scholar] [CrossRef] [Green Version]
  40. Gunia, S.; Erbersdobler, A.; Hakenberg, O.W.; Koch, S.; May, M. C-MET is expressed in the majority of penile squamous cell carcinomas and correlates with polysomy-7 but is not associated with MET oncogene amplification, pertinent histopathologic parameters, or with cancer-specific survival. Pathol.-Res. Pract. 2013, 209, 215–220. [Google Scholar] [CrossRef]
  41. Arya, M.; Thrasivoulou, C.; Henrique, R.; Millar, M.; Hamblin, R.; Davda, R.; Aare, K.; Masters, J.R.; Thomson, C.; Muneer, A.; et al. Targets of Wnt/ß-Catenin Transcription in Penile Carcinoma. PLoS ONE 2015, 10, e0124395. [Google Scholar] [CrossRef]
  42. Amancio, A.M.T.D.S.; Da Cunha, I.W.; Neves, J.I.; Quetz, J.D.S.; Carraro, D.M.; Rocha, R.M.; Zequi, S.C.; Cubilla, A.L.; Da Fonseca, F.P.; Lopes, A.; et al. Epidermal growth factor receptor as an adverse survival predictor in squamous cell carcinoma of the penis. Hum. Pathol. 2017, 61, 97–104. [Google Scholar] [CrossRef] [PubMed]
  43. Dorff, T.B.; Schuckman, A.K.; Schwartz, R.; Rashad, S.; Bulbul, A.; Cai, J.; Pinski, J.; Ma, Y.; Danenberg, K.; Skinner, E.; et al. Epidermal Growth Factor Receptor, Excision-Repair Cross-Complementation Group 1 Protein, and Thymidylate Synthase Expression in Penile Cancer. Clin. Genitourin. Cancer 2016, 14, 450–456.e1. [Google Scholar] [CrossRef] [PubMed]
  44. Ferrandiz-Pulido, C.; Masferrer, E.; Toll, A.; Hernandez-Losa, J.; Mojal, S.; Pujol, R.M.; Cajal, S.R.Y.; De Torres, I.; Garcia-Patos, V. mTOR Signaling Pathway in Penile Squamous Cell Carcinoma: pmTOR and peIF4E Over Expression Correlate with Aggressive Tumor Behavior. J. Urol. 2013, 190, 2288–2295. [Google Scholar] [CrossRef]
  45. Steinestel, J.; Al Ghazal, A.; Arndt, A.; Schnoeller, T.J.; Schrader, A.J.; Moeller, P.; Steinestel, K. The role of histologic subtype, p16(INK4a) expression, and presence of human papillomavirus DNA in penile squamous cell carcinoma. BMC Cancer 2015, 15, 220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Cubilla, A.L.; Lloveras, B.; Alejo, M.; Clavero, O.; Chaux, A.; Kasamatsu, E.; Monfulleda, N.; Tous, S.; Alemany, L.; Klaustermeier, J.; et al. Value of p16INK4a in the Pathology of Invasive Penile Squamous Cell Carcinomas. Am. J. Surg. Pathol. 2011, 35, 253–261. [Google Scholar] [CrossRef]
  47. Tang, D.H.; Clark, P.E.; Giannico, G.; Hameed, O.; Chang, S.S.; Gellert, L.L. Lack of P16 ink4a Over Expression in Penile Squamous Cell Carcinoma is Associated with Recurrence after Lymph Node Dissection. J. Urol. 2015, 193, 519–525. [Google Scholar] [CrossRef]
  48. 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]
  49. Gunia, S.; Kakies, C.; Erbersdobler, A.; Hakenberg, O.W.; Koch, S.; May, M. Expression of p53, p21 and cyclin D1 in penile cancer: p53 predicts poor prognosis. J. Clin. Pathol. 2011, 65, 232–236. [Google Scholar] [CrossRef]
  50. 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]
  51. Deng, C.; Li, Z.; Guo, S.; Chen, P.; Chen, X.; Zhou, Q.; Chen, J.; Yu, X.; Wu, X.; Ma, W.; et al. Tumor PD-L1 expression is correlated with increased TILs and poor prognosis in penile squamous cell carcinoma. OncoImmunology 2017, 6, e1269047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. 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] [PubMed] [Green Version]
  53. Cocks, M.; Taheri, D.; Ball, M.W.; Bezerra, S.M.; Rodriguez, M.D.C.; Ricardo, B.F.; Bivalacqua, T.J.; Sharma, R.B.; Meeker, A.; Chaux, A.; et al. Immune-checkpoint status in penile squamous cell carcinoma: A North American cohort. Hum. Pathol. 2017, 59, 55–61. [Google Scholar] [CrossRef] [PubMed]
  54. 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] [PubMed]
  55. Termini, L.; Fregnani, J.H.; Boccardo, E.; Da Costa, W.H.; Longatto-Filho, A.; Andreoli, M.A.; Costa, M.C.; Lopes, A.; Da Cunha, I.W.; Soares, F.A.; et al. SOD2 immunoexpression predicts lymph node metastasis in penile cancer. BMC Clin. Pathol. 2015, 15, 3. [Google Scholar] [CrossRef] [Green Version]
  56. Lippert, T.H.; Ruoff, H.-J.; Volm, M. Current Status of Methods to Assess Cancer Drug Resistance. Int. J. Med Sci. 2011, 8, 245–253. [Google Scholar] [CrossRef] [Green Version]
  57. Garnett, M.J.; McDermott, U. The evolving role of cancer cell line-based screens to define the impact of cancer genomes on drug response. Curr. Opin. Genet. Dev. 2014, 24, 114–119. [Google Scholar] [CrossRef] [Green Version]
  58. Cheung, P.F.-Y.; Yip, C.W.; Ng, L.W.-C.; Lo, K.W.; Wong, N.; Choy, K.W.; Chow, C.; Chan, K.F.; Cheung, T.T.; Poon, R.T.-P.; et al. Establishment and characterization of a novel primary hepatocellular carcinoma cell line with metastatic ability in vivo. Cancer Cell Int. 2014, 14, 103. [Google Scholar] [CrossRef] [Green Version]
  59. Ibarrola-Villava, M.; Cervantes, A.; Bardelli, A. Corrigendum to Preclinical models for precision oncology. BBACAN 1870/2 (2018) 239–246. Biochim. et Biophys. Acta (BBA)-Bioenerg. 2019, 1872, 188292. [Google Scholar] [CrossRef]
  60. Truffi, M.; Sorrentino, L.; Corsi, F. Fibroblasts in the Tumor Microenvironment. Adv. Exp. Med. Biol. 2020, 1234, 15–29. [Google Scholar] [CrossRef]
  61. Miyai, Y.; Esaki, N.; Takahashi, M.; Enomoto, A. Cancer-associated fibroblasts that restrain cancer progression: Hypotheses and perspectives. Cancer Sci. 2020, 111, 1047–1057. [Google Scholar] [CrossRef] [Green Version]
  62. Burd, E.M. Human Papillomavirus and Cervical Cancer. Clin. Microbiol. Rev. 2003, 16, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Guimarães, G.C.; Cunha, I.W.; Soares, F.A.; Lopes, A.; Torres, J.; Chaux, A.; Velazquez, E.F.; Ayala, G.; Cubilla, A.L. Penile Squamous Cell Carcinoma Clinicopathological Features, Nodal Metastasis and Outcome in 333 Cases. J. Urol. 2009, 182, 528–534. [Google Scholar] [CrossRef] [PubMed]
  64. Sanchez, D.F.; Soares, F.; Alvarado-Cabrero, I.; Cañete, S.; Fernández-Nestosa, M.J.; Rodríguez, I.M.; Barreto, J.; Cubilla, A.L. Pathological factors, behavior, and histological prognostic risk groups in subtypes of penile squamous cell carcinomas (SCC). Semin. Diagn. Pathol. 2015, 32, 222–231. [Google Scholar] [CrossRef]
  65. King, H.A.; Gerber, A.P. Translatome profiling: Methods for genome-scale analysis of mRNA translation. Briefings Funct. Genom. 2014, 15, 22–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Chaux, A.; Munari, E.; Katz, B.; Sharma, R.; Lecksell, K.; Cubilla, A.L.; Burnett, A.L.; Netto, G.J. The epidermal growth factor receptor is frequently overexpressed in penile squamous cell carcinomas: A tissue microarray and digital image analysis study of 112 cases. Hum. Pathol. 2013, 44, 2690–2695. [Google Scholar] [CrossRef]
  67. Di Lorenzo, G.; Buonerba, C.; Ferro, M.; Calderoni, G.; Bozza, G.; Federico, P.; Tedesco, B.; Ruggieri, V.; Aieta, M.; Federico, P. The epidermal growth factor receptors as biological targets in penile cancer. Expert Opin. Biol. Ther. 2014, 15, 473–476. [Google Scholar] [CrossRef] [Green Version]
  68. Gupta, S.; Sonpavde, G. Emerging Systemic Therapies for the Management of Penile Cancer. Urol. Clin. N. Am. 2016, 43, 481–491. [Google Scholar] [CrossRef]
  69. Huang, K.-B.; Liu, R.-Y.; Peng, Q.-H.; Li, Z.-S.; Jiang, L.-J.; Guo, S.-J.; Zhou, Q.-H.; Liu, T.-Y.; Deng, C.-Z.; Yao, K.; et al. EGFR mono-antibody salvage therapy for locally advanced and distant metastatic penile cancer: Clinical outcomes and genetic analysis. Urol. Oncol. Semin. Orig. Investig. 2019, 37, 71–77. [Google Scholar] [CrossRef] [PubMed]
  70. Necchi, A.; Vullo, S.L.; Perrone, F.; Raggi, D.; Giannatempo, P.; Calareso, G.; Nicolai, N.; Piva, L.; Biasoni, D.; Catanzaro, M.; et al. First-line therapy with dacomitinib, an orally available pan-HER tyrosine kinase inhibitor, for locally advanced or metastatic penile squamous cell carcinoma: Results of an open-label, single-arm, single-centre, phase 2 study. BJU Int. 2017, 121, 348–356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Parseghian, C.M.; Napolitano, S.; Loree, J.M.; Kopetz, S. Mechanisms of Innate and Acquired Resistance to Anti-EGFR Therapy: A Review of Current Knowledge with a Focus on Rechallenge Therapies. Clin. Cancer Res. 2019, 25, 6899–6908. [Google Scholar] [CrossRef]
  72. Kuasne, H.; Barros-Filho, M.C.; Marchi, F.A.; Drigo, S.A.; Scapulatempo-Neto, C.; Faria, E.F.; Rogatto, S.R. Nuclear loss and cytoplasmic expression of androgen receptor in penile carcinomas: Role as a driver event and as a prognosis factor. Virchows Archiv für Pathol. Anat. und Physiol. und für Klin. Med. 2018, 473, 607–614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Jacob, J.M.; Ferry, E.K.; Gay, L.M.; Elvin, J.A.; Vergilio, J.-A.; Ramkissoon, S.; Severson, E.; Necchi, A.; Killian, J.K.; Ali, S.M.; et al. Comparative Genomic Profiling of Refractory and Metastatic Penile and Nonpenile Cutaneous Squamous Cell Carcinoma: Implications for Selection of Systemic Therapy. J. Urol. 2019, 201, 541–548. [Google Scholar] [CrossRef] [PubMed]
  74. Chahoud, J.; Pickering, C.R.; Pettaway, C.A. Genetics and penile cancer. Curr. Opin. Urol. 2019, 29, 364–370. [Google Scholar] [CrossRef] [PubMed]
  75. McGregor, B.; Sonpavde, G. Immunotherapy for advanced penile cancer—Rationale and potential. Nat. Rev. Urol. 2018, 15, 721–723. [Google Scholar] [CrossRef]
  76. Fruman, D.A.; Chiu, H.; Hopkins, B.D.; Bagrodia, S.; Cantley, L.C.; Abraham, R.T. The PI3K Pathway in Human Disease. Cell 2017, 170, 605–635. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Huang, T.-X.; Guan, X.-Y.; Fu, L. Therapeutic targeting of the crosstalk between cancer-associated fibroblasts and cancer stem cells. Am. J. Cancer Res. 2019, 9, 1889–1904. [Google Scholar]
  78. Skoda, A.M.; Simovic, D.; Karin, V.; Kardum, V.; Vranic, S.; Serman, L. The role of the Hedgehog signaling pathway in cancer: A comprehensive review. Bosn. J. Basic Med. Sci. 2018, 18, 8–20. [Google Scholar] [CrossRef]
Cells 10 00814 g001 550
Figure 1. Morphological features, migration, and invasion characteristics of five penile cancer-derived cells. Cells were cultivated until passage P10. (A) Morphology of epithelial (cells 2 and 3) and fibroblast-like (cells 4, 5, and 6) cells by phase-contrast microscopy and (B) immunofluorescence (Texas Red: actin /phalloidin; FITC (fluorescein isothiocyanate): tubulin; and DAPI (4′,6-diamidino-2-phenylindole): nucleus) (10X magnification) (Nikon TE2000). Transwell migration and invasion assay (C,D) images and (E) quantification. The fields were evaluated in 4X magnification and were scored from 0 to 4 according to the percentage of cells that migrated and/or invaded (score 0: 0–5%, score 1: 5–25%, score 2: 25–50%, score 3: 50–75%, and score 4: 75–100%). Cells 2, 3, 4, and 6 presented high migratory capacity (scores 2.5 to 4), whereas cell 5 had a low score (score 0.5). Invasion capacity was observed in cells 3, 4, and 6. Chemotaxis was induced by adding 2.5% FBS (fetal bovine serum), BPE (bovine pituitary extract) (30ug/mL), and EGF (epithelial growth factor) (0.2 ng/mL) to the media, which is represented by (+). (−): control (media without chemo-attractants). Bars indicate the average of the scoring obtained in three independent experiments.
Figure 1. Morphological features, migration, and invasion characteristics of five penile cancer-derived cells. Cells were cultivated until passage P10. (A) Morphology of epithelial (cells 2 and 3) and fibroblast-like (cells 4, 5, and 6) cells by phase-contrast microscopy and (B) immunofluorescence (Texas Red: actin /phalloidin; FITC (fluorescein isothiocyanate): tubulin; and DAPI (4′,6-diamidino-2-phenylindole): nucleus) (10X magnification) (Nikon TE2000). Transwell migration and invasion assay (C,D) images and (E) quantification. The fields were evaluated in 4X magnification and were scored from 0 to 4 according to the percentage of cells that migrated and/or invaded (score 0: 0–5%, score 1: 5–25%, score 2: 25–50%, score 3: 50–75%, and score 4: 75–100%). Cells 2, 3, 4, and 6 presented high migratory capacity (scores 2.5 to 4), whereas cell 5 had a low score (score 0.5). Invasion capacity was observed in cells 3, 4, and 6. Chemotaxis was induced by adding 2.5% FBS (fetal bovine serum), BPE (bovine pituitary extract) (30ug/mL), and EGF (epithelial growth factor) (0.2 ng/mL) to the media, which is represented by (+). (−): control (media without chemo-attractants). Bars indicate the average of the scoring obtained in three independent experiments.
Cells 10 00814 g001
Cells 10 00814 g002 550
Figure 2. Genomic comparison of penile cancer (PeCa) samples and their derived cells (A) Sample identity distogram showing the degree of similarity among penile cancer-derived cells and their primary tumor tissue. (B) Variants identified by targeted next-generation sequencing in PeCa and derived PeCa cells. DNA from PeCa 3 and 4 were not available for this analysis. Oncoprint was generated by the cBioportal Oncoprinter tool (http://www.cbioportal.org/oncoprinter, accessed on 6 April 2021).
Figure 2. Genomic comparison of penile cancer (PeCa) samples and their derived cells (A) Sample identity distogram showing the degree of similarity among penile cancer-derived cells and their primary tumor tissue. (B) Variants identified by targeted next-generation sequencing in PeCa and derived PeCa cells. DNA from PeCa 3 and 4 were not available for this analysis. Oncoprint was generated by the cBioportal Oncoprinter tool (http://www.cbioportal.org/oncoprinter, accessed on 6 April 2021).
Cells 10 00814 g002
Cells 10 00814 g003 550
Figure 3. Molecular profile of penile cancer (PeCa)-derived cells. (A) Polysomal profiling of cells 2, 3, 4, 5, and 6. Polyribosome (polysome) fractionation was performed using a sucrose density gradient (5–50%). (B) Translatomic profiling (Clariom™ D Assay platform, Thermo Fisher Scientific, Waltham, MA, USA) was performed in duplicate for each derived PeCa cell culture and cell 1 (normal foreskin). A heatmap shows three clusters with probes differentially expressed compared to cell 1. The first cluster is composed of cell 3 replicates (tumor epithelial cells), cluster 2 is composed of cell 2 replicates (tumor epithelial cells), and the third one of cells 4, 5, and 6 (cancer-associated fibroblasts—CAFs). Blue: down-expressed probes, red: overexpressed probes. Genes indicated in the right of the heatmap are known markers of CAFs. (C,D) Gene set enrichment analysis was performed by ssGSEA (single-sample gene set enrichment analysis: https://www.genepattern.org/modules/docs/ssGSEAProjection/4, accessed on 6 April 2021). Oncogenic signature and Reactome sub-collection came from the Molecular Signature Database (MySigDB) (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp, accessed on 6 April 2021). Red: highest score in the ssGSEA, blue: lowest score. UP: up-regulated genes; DN: downregulated genes; V1: version 1. (E) Heatmap of reverse phase protein array (RPPA) representing proteins or phosphorylated isotypes in PeCa cells. The protein or phosphorylated isotype expression levels in tumor cells was compared to cell 1 (fold change—FC). The FC levels of each biomarker of the five penile cancer-derived cells are indicated in colors (red for the highest and blue for the lowest). (F) Western blot showing EGFR overexpression in cell 2 compared to the other PeCa cells.
Figure 3. Molecular profile of penile cancer (PeCa)-derived cells. (A) Polysomal profiling of cells 2, 3, 4, 5, and 6. Polyribosome (polysome) fractionation was performed using a sucrose density gradient (5–50%). (B) Translatomic profiling (Clariom™ D Assay platform, Thermo Fisher Scientific, Waltham, MA, USA) was performed in duplicate for each derived PeCa cell culture and cell 1 (normal foreskin). A heatmap shows three clusters with probes differentially expressed compared to cell 1. The first cluster is composed of cell 3 replicates (tumor epithelial cells), cluster 2 is composed of cell 2 replicates (tumor epithelial cells), and the third one of cells 4, 5, and 6 (cancer-associated fibroblasts—CAFs). Blue: down-expressed probes, red: overexpressed probes. Genes indicated in the right of the heatmap are known markers of CAFs. (C,D) Gene set enrichment analysis was performed by ssGSEA (single-sample gene set enrichment analysis: https://www.genepattern.org/modules/docs/ssGSEAProjection/4, accessed on 6 April 2021). Oncogenic signature and Reactome sub-collection came from the Molecular Signature Database (MySigDB) (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp, accessed on 6 April 2021). Red: highest score in the ssGSEA, blue: lowest score. UP: up-regulated genes; DN: downregulated genes; V1: version 1. (E) Heatmap of reverse phase protein array (RPPA) representing proteins or phosphorylated isotypes in PeCa cells. The protein or phosphorylated isotype expression levels in tumor cells was compared to cell 1 (fold change—FC). The FC levels of each biomarker of the five penile cancer-derived cells are indicated in colors (red for the highest and blue for the lowest). (F) Western blot showing EGFR overexpression in cell 2 compared to the other PeCa cells.
Cells 10 00814 g003
Cells 10 00814 g004 550
Figure 4. Enrichment analysis of direct and indirect mRNA interactions and proteins dysregulated in penile cancer (PeCa) cells show EGFR as a core molecule. Sensitivity assays to EGFR inhibitors are depicted for cells 2 and 3. Upstream regulator analysis (IPA, Ingenuity Pathway Analysis) indicates direct (solid lines) and indirect (dashed lines) interaction of mRNA (A) and proteins, which is exemplified in PeCa-derived cell 2. (B) EGFR is a core molecule in both mRNA and protein analyses. (C) Chemosensitivity assays were performed in cell 2 and cell 3 for the identification of IC50 and dose-response curves. Cells were treated with cisplatin and EGFR inhibitors (cetuximab, gefitinib, and erlotinib). Cell 2 and cell 3 were sensitive to cisplatin. Only cell 2 (EGFR overexpression) was sensitive to EGFR inhibitors.
Figure 4. Enrichment analysis of direct and indirect mRNA interactions and proteins dysregulated in penile cancer (PeCa) cells show EGFR as a core molecule. Sensitivity assays to EGFR inhibitors are depicted for cells 2 and 3. Upstream regulator analysis (IPA, Ingenuity Pathway Analysis) indicates direct (solid lines) and indirect (dashed lines) interaction of mRNA (A) and proteins, which is exemplified in PeCa-derived cell 2. (B) EGFR is a core molecule in both mRNA and protein analyses. (C) Chemosensitivity assays were performed in cell 2 and cell 3 for the identification of IC50 and dose-response curves. Cells were treated with cisplatin and EGFR inhibitors (cetuximab, gefitinib, and erlotinib). Cell 2 and cell 3 were sensitive to cisplatin. Only cell 2 (EGFR overexpression) was sensitive to EGFR inhibitors.
Cells 10 00814 g004
Cells 10 00814 g005 550
Figure 5. Approximately 30% of penile cancer overexpresses genes related to the EGFR signaling pathway. Heatmap illustrating the Reactome signaling by EGFR (50 genes) (A) and KEGG ERBB signaling pathway (87 genes) (B) signature of 36 PeCa samples (GSE57955). Rows represent normalized mRNA expression for a single gene (Agilent microarray, Agilent, Santa Clara, CA, USA), and each column represents one of the 36 previously analyzed PeCa samples. All tumor samples were normalized with a pool of normal samples, and the values were ranked by their expression levels (blue: down-expressed genes, red: overexpressed genes). The red cluster (A,B) presents overexpression of several genes belonging to the EGFR signature. (C) Heatmap representing a matrix of the Pearson correlation of genes from the Reactome signaling by EGFR in the 36 PeCa samples (GSE57955). EGFR (first column) was positively correlated with several genes of the PIK3 and MAPK pathways. A significant P-value of Pearson correlation is represented by (*) only for the EGFR gene (first line).
Figure 5. Approximately 30% of penile cancer overexpresses genes related to the EGFR signaling pathway. Heatmap illustrating the Reactome signaling by EGFR (50 genes) (A) and KEGG ERBB signaling pathway (87 genes) (B) signature of 36 PeCa samples (GSE57955). Rows represent normalized mRNA expression for a single gene (Agilent microarray, Agilent, Santa Clara, CA, USA), and each column represents one of the 36 previously analyzed PeCa samples. All tumor samples were normalized with a pool of normal samples, and the values were ranked by their expression levels (blue: down-expressed genes, red: overexpressed genes). The red cluster (A,B) presents overexpression of several genes belonging to the EGFR signature. (C) Heatmap representing a matrix of the Pearson correlation of genes from the Reactome signaling by EGFR in the 36 PeCa samples (GSE57955). EGFR (first column) was positively correlated with several genes of the PIK3 and MAPK pathways. A significant P-value of Pearson correlation is represented by (*) only for the EGFR gene (first line).
Cells 10 00814 g005
Table
Table 1. Clinical and pathological characteristics of penile cancer patients.
Table 1. Clinical and pathological characteristics of penile cancer patients.
Patient/
Cell Line		Age (Years)Histological SubtypeSurgery TypeHPVTNMPerineural InvasionClinical StageFollow-up (Months)
2/Cell 285Usual + SarcomatoidPartial penectomyHPV16T3N0M0YesIII120
3/Cell 371VerrucousPartial penectomyNegativeT1N0M0NoI138
4/Cell 443UsualPartial penectomyNegativeT1N0M0YesI102
5/Cell 557UsualTotal penectomyNegativeT4N0M1YesIV30
6/Cell 670BasaloidPartial penectomyNegativeT2N1M2NoIV18
HPV—human papillomavirus; TNM—TNM Staging System (T: Tumor size, N: regional lymph nodes, M: Distant metastasis).
Table
Table 2. Variants identified by targeted next-generation sequencing in the penile cancer-derived cell cultures.
Table 2. Variants identified by targeted next-generation sequencing in the penile cancer-derived cell cultures.
IDGeneClassificationChr:
Location 		Type of AlterationdbSNPTranscriptBase Change
Cell 2/tumor 2PIK3CAP3:17,8936,082Missensers121913273NM_006218.4c.1624G > A
Cell 2/tumor 2PIK3CAP3:178,936,093Missensers121913275NM_006218.4c.1635G > C
Cell 2/tumor 2FGFR1VUS8:38,270,403Otherrs1364534792NM_023110.3c.*743dupA
Cell 2MMP1VUS11:102,667,445Missense NM_002421.4c.575G > C
Cell 2STAT3P17:40,486,045LoF NM_139276.2c.820C > T
Cell 3-
Cell 4RAD50VUS5:131,953,850Missensers143189763NM_005732.4c.3253A > G
Cell 4MMP1VUS11:102,668,717LoFrs139018071NM_002421.4c.105+2T > C
Cell 4FLT3LP13:28,589,804Missensers903856095NM_004119.3c.2576G > A
Cell 5/tumor 5COL11A1LP1:103,380,339Missense NM_001854.4c.3845G > T
Cell 5/tumor 5ALKVUS2:29,455,260Missense NM_004304.5c.2542G > A
Cell 5/tumor 5PIK3CAP3:178,952,085Missensers121913279NM_006218.4c.3140A > G
Cell 5/tumor 5FBXW7LP4:153,247,367Missensers747241612NM_001349798.2c.1435C > G
Cell 5/tumor 5NOTCH1P9:139,413,064LoF NM_017617.5c.1078G > T
Cell 5/tumor 5BRCA2P13:32,968,951LoFrs80359212NM_000059.3c.9382C > T
Cell 5/tumor 5IGF1RVUS15:99507206Other NM_000875.5c.*6535T > G
Cell 5/tumor 5NF1VUS17:29,702,854Otherrs909909591NM_001042492.3c.*1683_*1685delGAA
Cell 5/tumor 5SMAD4LP18:48,575,116Missense NM_005359.6c.310C > T
Cell 6/tumor 6Near APCVUS5:112,043,188upstream transcript variantrs1554060178NC_000005.10112707490:C:G
Tumor 6MLH3VUS14:75,483,802Missense NM_001040108.1c.4345C > T
Tumor 6ERBB2LP17:37,866,662Missense NM_004448.3c.829G > T
Chr: chromosome. c. represents the coding sequence position. LoF: Loss of function. P: Pathogenic. LP: Likely pathogenic. VUS: variant of uncertain significance. dbSNP: The Single Nucleotide Polymorphism Database
Table
Table 3. Top 15 differentially expressed genes in penile cancer cells compared with the normal reference cell (cell 1). Translatomic profiling was performed using the microarray platform Clariom™ D Assay platform (Thermo Fisher Scientific, Waltham, MA, USA) after polysomal mRNA enrichment. The number within brackets represents differentially expressed genes after enrichment analysis, indicating genes more likely to be translated into proteins.
Table 3. Top 15 differentially expressed genes in penile cancer cells compared with the normal reference cell (cell 1). Translatomic profiling was performed using the microarray platform Clariom™ D Assay platform (Thermo Fisher Scientific, Waltham, MA, USA) after polysomal mRNA enrichment. The number within brackets represents differentially expressed genes after enrichment analysis, indicating genes more likely to be translated into proteins.
Cell 2 (564)Cell 3 (1199)Cell 4 (262)Cell 5 (205)Cell 6 (163)
Top 15 Overexpressed Genes
GeneFCGeneFCGeneFCGeneFCGeneFC
S100A854897.5OSR11865.57SFRP21936.2IL24305.3TMEM176B2376.2
COL17A135757.4FOXG11259.5TRPA1964.7APOD101.4TMEM176A1686.2
KRT525856.7PRKAR2B619.61LINC01436458.9SFRP297.6APOD935.5
S100A925339.7HOXD8333.72TMEM176B288.4STMN283.6SLC14A1873.7
LCN223704.2LGALSL172.59TMEM176A252.5EGR179.7XGY2446.0
CDH118954.7PORCN119.5XGY2247.8XG66.9XG374.6
LAMC217269.9NRN186.36SAT1113.8SNORD50A59.6CCND2343.2
MAL215941.8TBX1574.35XG100.1ZBTB1651.1FAM105A213.8
GJB214969.1PRPF3963.2IL3398.9TNFRSF2144.0SAA1163.9
FXYD313924.8MSI241.57APOD98.2DUXAP1040.9STMN2160.7
KRT1711118.8SLC44A140.69OSR188.5USP5337.3AQP1158.7
DSG310259.7STX636.88PDGFRL72.3LINC0129631.6IL13RA2155.7
FGFBP110054.2CUL334.26CCND265.9RGCC25.6DUSP6147.2
PI310052.7S100A732.87CLU62.3SFRP122.5FMO3138.3
S100A29269.6RTCB29.47FGF758.3EIF3A21.5ITGA8138.2
Top 15 Down-Expressed Genes
GeneFCGeneFCGeneFCGeneFCGeneFC
BGN–15958.7CCL2–24130.3GFRA1–1175.8GFRA1–1156.6GFRA1–738.2
GREM1–14798.9CLDN11–9744.1ADAM12–1155.9ADAM12–615.7CPA3–522.3
UCHL1–13551.2UCHL1–9705.2IGFBP3–556.8POSTN–585.0POSTN–468.1
MGST1–13191.8LUM–9606.2PSG5–416.9CPA3–422.7CXCL12–169.9
THY1–11417.6SULF1–7813.0CPA3–411.1UCP2–309.4PLPP4–147.7
MXRA8–8181.8HIST1H2BM–7623.5PLPP4–378.8UCHL1–255.8PSG2–119.4
LOXL1–8004.0MXRA8–6464.9CNN1–365.3ACKR3–209.7PSG1–106.9
ENG–7918.5PSG5–6399.8CLDN11–275.9CNN1–189.1HIST1H3F–96.9
FBLN5–6160.9MGST1–6082.7PLPPR4–245.8COL1A1–133.3BEX1–95.9
MFAP4–3745.0MRPL20–4740.9UCP2–243.7HIST1H3F–114.1PSG8–66.6
ADGRA2–3508.0FARP1–4690.7LBH–220.6KRT7–79.1SHOX–62.4
COL1A2–2748.2THY1–4498.7PSG8–167.9MYBL2–73.4UCP2–55.3
GFRA1–2705.4TGFBR1–4386.3SLC7A5–162.6LYPD6B–63.9DAPK1–49.0
F2RL2–2520.2HIST2H4B–4137.9PSG11–160.0SEL1L3–53.0KCND3–45.9
MYADM–2340.8CSRP1–4129.1MEST–138.6KCND3–44.6COLEC12–45.2
FC: fold change.
Table
Table 4. Summary of the studies in penile cancer that evaluated the expression (immunohistochemistry or immunofluorescence) of the same proteins presented in our RPPA (reverse-phase protein arrays) analysis.
Table 4. Summary of the studies in penile cancer that evaluated the expression (immunohistochemistry or immunofluorescence) of the same proteins presented in our RPPA (reverse-phase protein arrays) analysis.
ProteinNumber of CasesExpressionReferenceCell
23456
Akt1 148↑Stankiewicz et al., 2011 [36]↓↑-↑-
pAkt 112↑Chaux et al., 2014 [37]-↑-↑↓
pAkt 148↑Stankiewicz et al., 2011 [36]-↑-↑↓
pAkt57↑Azizi et al., 2019 [38]-↑-↑↓
ARID1A112↑Faraj et al., 2015 [39]↑↑↓↑↑
c-MET92↑Gunia et al., 2013 [40]-----
c-MYC141↓Arya et al., 2015 [41]----↓
cyclinD1141↓Arya et al., 2015 [41]-----
EGFR139↑Silva Amancio et al., 2017 [42]↑--↑-
EGFR52↑Dorff et al., 2016 [43]↑--↑-
EI2F13↑Fenner et al., 2018 [20]↑↑-↑-
eIF4E67↑Ferrandiz-Pulido et al., 2013 [44]↑↑--↓
peIF4E67↑Ferrandiz-Pulido et al., 2013 [44]-----
HER2148↑Stankiewicz et al., 2011 [36]↑↑---
HER3148↑Stankiewicz et al., 2011 [36]↑↑---
HER4148↑Stankiewicz et al., 2011 [36]-----
p1658↑Steinestel et al., 2015 [45]↓↓-↓↓
p16202↑Cubilla et al., 2011 [46]↓↓-↓↓
p16119↑Tang et al., 2015 [47]↓↓-↓↓
p16123↑Mannweiler et al., 2013 [48]↓↓-↓↓
p4E-BP167↑Ferrandiz-Pulido et al., 2013 [44]↑↑---
p53123↑Mannweiler et al., 2013 [48]-----
p53110↑Gunia et al., 2013 [49]-----
p53297↑Rocha et al., 2012 [50]-----
PDL1116↑Deng et al., 2017 [51]-↑---
PDL137↑Udager et al., 2016 [52]-↑---
PDL152↑Cocks et al., 2017 [53]-↑---
PDL1213↑Ottenhof et al., 2017 [54]-↑---
pmTOR67↑Ferrandiz-Pulido et al., 2013 [44]↑---↑
pmTOR112↑Chaux et al., 2014 [37]↑---↑
PTEN112↓Chaux et al., 2014 [37]-----
PTEN148↓Stankiewicz et al., 2011 [36]↓↑-↓↓
PTEN57↓Azizi et al., 2019 [38]↓↑-↓↓
pS657↑Azizi et al., 2019 [38]↑--↑-
SOD2125↑Termini et al., 2015 [55]↑↑↑↑-
↑: increased expression; ↓: decreased expression; -: no change in expression levels.
Table
Table 5. Drugs that target receptor tyrosine kinase genes and proteins. The tyrosine kinases showed high expression levels (fold change > 2) in penile cancer-derived cells 2 and 3. Data were generated using the Ingenuity Pathway Analysis software (Qiagen, Valencia, CA, USA).
Table 5. Drugs that target receptor tyrosine kinase genes and proteins. The tyrosine kinases showed high expression levels (fold change > 2) in penile cancer-derived cells 2 and 3. Data were generated using the Ingenuity Pathway Analysis software (Qiagen, Valencia, CA, USA).
Cell LineGene Symbol amRNA FCLocationDrug(s) b
mRNA
Cell 2FRK144.007NucleusRegorafenib
Cell 2ERBB3 c121.938Plasma membraneAfatinib/cetuximab/osimertinib/erlotinib
Cell 2MST1R54.569Plasma membraneCrizotinib/erlotinib/gefitinib/pazopanib
Cell 2PRKDC43.713NucleusCC-115/MSC2490484A/panulisib
Cell 2EGFRc38.586Plasma membraneGefitinib/erlotinib/afatinib/cetuximab
Cell 2MET30.91Plasma membraneCrizotinib/cabozantinib/ABT-700/altiratinib
Cell 3SRC28.443CytoplasmDasatinib/bosutinib/blinatumomab/AZD0424
Cell 2DDR124.59Plasma membraneBlinatumomab/dacomitinib/dasatinib/imatinib
RPPA
Cell 2KDR4.729Plasma membrane5-azacytidine/sorafenib/sAEE788/apatinib
Cell 2RPS6KA13.775CytoplasmPMD-026
Cell 2MAP2K13.493CytoplasmARRY-424704/AS703988/binimetinib
Cell 2MTOR2.925NucleusABI-009/apitolisib/sirolimus/tacrolimus/everolimus
Cell 2SRC2.792CytoplasmAZD0424/dasatinib/bosutinib/ponatinib
Cell 3CDK12.637NucleusAlvocidib/dinaciclib/milciclib/riviciclib
Cell 2MKNK12.62CytoplasmBAY1143269/dacomitinib/ETC-1907206/tomivosertib
Cell 3AKT12.386CytoplasmBAY1125976/capivasertib/patasertib/miransertib
Cell 2BRD42.241NucleusAZD5153/BI 894999/PLX2853/PLX51107
Cell 2CDK12.051NucleusAlvocidib/dinaciclib/milciclib/riviciclib
a Only kinases were listed. b A maximum of four drugs was listed. c EGFR and ERBB3 presented gene and protein overexpression. Only the fold change (FC) found in the mRNA analysis is listed.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kuasne, H.; Canto, L.M.d.; Aagaard, M.M.; Muñoz, J.J.M.; Jamblinne, C.D.; Marchi, F.A.; Scapulatempo-Neto, C.; Faria, E.F.; Lopes, A.; Carréno, S.; et al. Penile Cancer-Derived Cells Molecularly Characterized as Models to Guide Targeted Therapies. Cells 2021, 10, 814. https://doi.org/10.3390/cells10040814

AMA Style

Kuasne H, Canto LMd, Aagaard MM, Muñoz JJM, Jamblinne CD, Marchi FA, Scapulatempo-Neto C, Faria EF, Lopes A, Carréno S, et al. Penile Cancer-Derived Cells Molecularly Characterized as Models to Guide Targeted Therapies. Cells. 2021; 10(4):814. https://doi.org/10.3390/cells10040814

Chicago/Turabian Style

Kuasne, Hellen, Luisa Matos do Canto, Mads Malik Aagaard, Juan Jose Moyano Muñoz, Camille De Jamblinne, Fabio Albuquerque Marchi, Cristovam Scapulatempo-Neto, Eliney Ferreira Faria, Ademar Lopes, Sébastien Carréno, and et al. 2021. "Penile Cancer-Derived Cells Molecularly Characterized as Models to Guide Targeted Therapies" Cells 10, no. 4: 814. https://doi.org/10.3390/cells10040814

APA Style

Kuasne, H., Canto, L. M. d., Aagaard, M. M., Muñoz, J. J. M., Jamblinne, C. D., Marchi, F. A., Scapulatempo-Neto, C., Faria, E. F., Lopes, A., Carréno, S., & Rogatto, S. R. (2021). Penile Cancer-Derived Cells Molecularly Characterized as Models to Guide Targeted Therapies. Cells, 10(4), 814. https://doi.org/10.3390/cells10040814

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