Next Article in Journal
Dual Roles of Astrocyte-Derived Factors in Regulation of Blood-Brain Barrier Function after Brain Damage
Next Article in Special Issue
Metabolic and Molecular Basis of Sarcopenia: Implications in the Management of Urothelial Carcinoma
Previous Article in Journal
Deciduous DPSCs Ameliorate MPTP-Mediated Neurotoxicity, Sensorimotor Coordination and Olfactory Function in Parkinsonian Mice
Previous Article in Special Issue
Exercise Training Impacts Cardiac Mitochondrial Proteome Remodeling in Murine Urothelial Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 3
10.3390/ijms20030570
ijms-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

Transcriptional Analysis of Immunohistochemically Defined Subgroups of Non-Muscle-Invasive Papillary High-Grade Upper Tract Urothelial Carcinoma

by
Minsun Jung
1,
Jeong Hoon Lee
2,
Bohyun Kim
1,
Jeong Hwan Park
1,3 and
Kyung Chul Moon
1,4,*
1
Department of Pathology, Seoul National University College of Medicine, Seoul 03080, Korea
2
Seoul National University Biomedical Informatics (SNUBI), Division of Biomedical Informatics and Systems Biomedical Informatics National Core Research Center, Seoul National University College of Medicine, Seoul 03080, Korea
3
Department of Pathology, SMG-SNU Boramae Medical Center, Seoul 03080, Korea
4
Kidney Research Institute, Medical Research Center, Seoul National University College of Medicine, Seoul 03080, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(3), 570; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030570
Submission received: 27 December 2018 / Revised: 26 January 2019 / Accepted: 28 January 2019 / Published: 29 January 2019
(This article belongs to the Special Issue Molecular Research Efforts in Urothelial Carcinoma)
Browse Figures
Versions Notes

Abstract

:
Immunohistochemical (IHC) staining for CK5/6 and CK20 was reported to be correlated with the prognosis of early urothelial carcinoma in a way contrary to that of advanced tumors for unknown reasons. We aimed to characterize the gene expression profiles of subgroups of non-muscle-invasive papillary high-grade upper tract urothelial carcinoma (UTUC) classified by CK5/6 and CK20 expression levels: group 1 (CK5/6-high/CK20-low), group 2 (CK5/6-high/CK20-high), and group 3 (CK5/6-low/CK20-high). Expression of group 3 was predictive of worse prognosis of non-muscle-invasive papillary high-grade UTUC. Transcriptional analysis revealed 308 differentially expressed genes across the subgroups. Functional analyses of the genes identified cell adhesion as a common process differentially enriched in group 3 compared to the other groups, which could explain its high-risk phenotype. Late cell cycle/proliferation signatures were also enriched in group 3 and in some of the other groups, which may be used as a prognostic biomarker complementary to CK5/6 and CK20. Group 2, characterized by low levels of genes associated with mitogen-activated protein kinase and tumor necrosis factor signaling pathways, was hypothesized to represent the least cancerous subtype considering its normal urothelium-like IHC pattern. This study would facilitate the application of easily accessible prognostic biomarkers in practice.

Graphical Abstract

1. Introduction

Recent advances in genetic and bioinformatic technologies have increased our understanding of the molecular characteristics of urothelial carcinoma. Based on gene expression profiles, it was demonstrated that both muscle-invasive (MIBC) and non-muscle-invasive urinary bladder carcinoma (NMIBC) can be divided into intrinsic molecular subtypes [1,2,3,4,5,6]. Although the classifications of MIBC were named differently by respective research groups, they had significant overlap: common luminal and basal molecular subtypes were reported [7]. Luminal type tumors expressed terminal differentiation markers of the urothelium, such as KRT20, PPARG, GATA3, FOXA1, and UPKs [2,5]. The basal type was enriched with basal cytokeratin, squamous differentiation, and stem cell signature markers, including KRT5, KRT6, KRT14, and CD44 [2,3,5]. Molecular subtypes showed differences in prognosis, responsiveness to neoadjuvant chemotherapy, and targetable mutations; therefore, they would significantly influence MIBC treatment [5,8,9,10]. On the other hand, gene expression profiles of NMIBC showed 3 different clusters, named as class 1, class 2, and class 3, which differed in biologic signatures and in prognosis [6]. For example, class 1 tumors had high expression of early cell cycle genes and luminal type markers, class 2 was enriched with a late cell cycle signature and with the luminal type genes, and class 3 was characterized by high expression of basal type genes and long non-coding RNAs [6]. Despite elevated levels of luminal type markers, including KRT20, class 2 tumors showed poor prognosis, which indicated that class 2 may represent a molecular shift toward muscle-invasive carcinoma by taking the carcinoma in situ (CIS) pathway [6]. In another study on stage Ta urinary bladder carcinoma, cancers were segregated by copy number alteration profiles into two subtypes which varied in mammalian target of rapamycin complex 1 activity, metabolic pathways, and mutations [11].
Upper tract urothelial carcinoma (UTUC) is much less common than urinary bladder carcinoma, accounting for 5%–10% of total urothelial carcinoma, with its incidence estimated at 1/50,000 in the Western population [12]. UTUC showed high mortality rate with 5-year survival <50% for muscle invasive tumor and <10% for pT4 stage disease [12]. UTUC was revealed to have comparable types of mutations as urinary bladder carcinoma by targeted sequencing, but with different prevalence [13]. Comparison of public microarray data, including 12 UTUC and 20 urinary bladder carcinomas, demonstrated that both tumors were unidentifiable at the gene expression level, except for pT3 tumors [14]. It was also reported that UTUC was more enriched with luminal type genes compared with urinary bladder carcinoma [10,14]. However, the genomic landscape of UTUC has not been fully elucidated due to its rarity [15].
Several studies reported biomarkers of urothelial carcinoma using immunohistochemical (IHC) staining by reflecting its molecular frameworks, which were suggested to be relevant to prognosis and responsiveness to systematic treatment [16]. Luminal and basal type urinary bladder carcinoma were positive for luminal (CK20, GATA3, and uroplakin) and basal (CK5/6, CK14, and CD44) IHC markers, respectively [2,17]. In MIBC, positivity for CK5/6 and negativity for CK20 were correlated with poor survival [2]. In addition, high expression of CK5/6 and CK14 combined with low expression of FOXA1 and GATA3 in IHC staining was known to indicate a subtype with aggressive behavior or a basal/squamous-like (BASQ) subtype [18]. Interestingly, subgroups of non-muscle-invasive UTUC classified by CK5/6 and CK20 immunoprofiles showed contradicting prognoses to those of MIBC: CK5/6-low/CK20-high early UTUC had the worst outcome [19,20]. As a result, these IHC-defined subgroups of early UTUC supposedly have unique genetic features, the understanding of which would aid in comprehending the underlying pathobiology of non-muscle-invasive UTUC and facilitate the implementation of these IHC-based risk-stratifications in practice. Here, we analyzed the gene expression profiles of non-muscle-invasive papillary high-grade UTUC and revealed that CK5/6 and CK20 expression was not entirely relevant to molecular subtypes of this tumor. Instead, we postulated that cell adhesion and late cell cycle/proliferation processes were important to clinical behavior of subgroups defined by IHC staining for CK5/6 and CK20. In addition, tumors with CK5/6-high/CK20-high expression were suggested as a less cancerous subgroup.

2. Results

2.1. Subgroups of Non-Muscle-Invasive High-Grade Papillary UTUC Based on IHC Staining for CK5/6 and CK20

IHC staining for CK5/6, CK20, CK14, GATA3, FOXA1, CD44, and p53 was carried out for 50 non-muscle-invasive papillary UTUC tumors (Table S1). CK14 expression was almost absent in 98% of evaluable cases. By contrast, 95% and 76% of evaluable cases stained for GATA3 and FOXA1 in >80% of tumor cells, respectively. Correlations of IHC staining with pT stage, WHO grade, and the presence of CIS were analyzed. Among these characteristics, only WHO grade showed a significant association with IHC profiles: high WHO grade was associated with low CK5/6 (p = 0.004) and CD44 (p = 0.048), marginally with high CK20 (p = 0.083) and p53 (p = 0.051) expression, but not with CK14 (p = 1.000), GATA3 (p = 0.339), and FOXA1 (p = 0.778) staining (Figure 1A and Figure S1).
Further analysis was limited to high-grade tumors to remove grade-related genetic diversity [16]. We tried to classify non-muscle-invasive papillary high-grade UTUC using IHC staining for CK5/6 and CK20 as surrogate markers, which were known to be related to molecular subtypes and to the prognosis of UTUC and of urinary bladder carcinoma [19,20,21,22]. Fresh tissue with high (IHC score ≥ 6) or low (IHC score = 1) CK5/6 and CK20 expression was selected for RNA sequencing (RNA-seq). The same IHC criteria were marginally correlated with the progression-free survival (PFS) of the patients with non-muscle-invasive papillary high-grade UTUC (p = 0.071) retrieved from the published cohort [20]: tumors with CK5/6-low/CK20-high expression tend to show the worst PFS (Figure S2). Finally, three subgroups were established as follows: group 1, CK5/6-high/CK20-low; group 2, CK5/6-high/CK20-high; and group 3, CK5/6-low/CK20-high (Figure 1B). CK5/6 and CK20 were frequently stained in the whole layer of tumors in groups 1 and 3, respectively, without noticeable compartmentalization (Figure 1B). Although group 2 tumors were positive to both CK5/6 and CK20 in >50% of tumor cells, the expression of CK5/6 and CK20 was accentuated in basal and luminal cells, leaving out at least one cell layer of the basal and luminal portion, respectively, in all cases (Figure 1B). Clinicopathological data of the patients are summarized in Table 1. The median age of the patients was 69 years (range, 56–84) and the male/female sex ratio was 2:1. The tumors measured 3.6 ± 3.24 cm (mean ± standard deviation) in maximal diameter. Six (40%) patients were in stage pTa, and the other 9 (60%) patients were in pT1. CIS was observed in 5 cases (33.3%). There were no significant differences in clinicopathological parameters or in IHC profiles among the subgroups, except for CK5/6 and CK20 expression (Table 1). One sample that belonged to group 1 showed 20%–30% positivity for CK14. With the exception of this case, IHC staining for GATA3 and FOXA1 showed diffuse staining in all samples.

2.2. Identification of Differentially Expressed Genes (DEGs) between the Subgroups

Using p-value <0.01 and |fold change| >2 as the cut-offs, 132, 115, and 143 DEGs were identified between groups 1 and 2, groups 1 and 3, and groups 2 and 3, respectively (Figure 2A,B and Table S2). Group 1 had 58 upregulated and 74 downregulated genes compared to group 2. Group 3, relative to group 1 and 2, harbored increased levels of 57 and 58 genes and decreased levels of 52 and 91 genes, respectively. Consequently, a total of 308 DEGs across all comparisons were found. They were put into unsupervised hierarchical clustering that demonstrated the congregation of the subgroups defined by IHC staining (Figure 2C). Principal component analysis of the DEGs revealed that tumors with similar CK5/6 and CK20 expression profiles made clusters separated from one another in 2D space (Figure S3).

2.3. Functional Analysis of the Subgroups

Figure 3 shows the top 10 significant (false discovery date (FDR) < 0.05) results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of DEGs between the subgroups.

2.3.1. Group 3 was Enriched for Processes Related to Cellular Adhesion and Motility

Notably, DEGs between groups 3 and 1, or between groups 3 and 2 were enriched in cell adhesion and motility in GO analyses, including ‘cell adhesion’ (FDR < 0.001), ‘biological adhesion’ (FDR < 0.001), ‘cell–cell adhesion’ (FDR < 0.001), ‘regulation of cell adhesion’ (FDR = 0.015), and ‘cell junction’ (FDR = 0.021) between groups 3 and 1, and ‘cell adhesion’ (FDR < 0.001), ‘biological adhesion’ (FDR < 0.001), ‘cell–cell adhesion’ (FDR < 0.001), ‘cell–cell adhesion via plasma-membrane adhesion molecules’ (FDR = 0.028), ‘cell motility’ (FDR = 0.045), and ‘regulation of cell adhesion’ (FDR = 0.046) between groups 3 and 2. These processes involved several DEGs, for example, CASK, MYO10, SPINK5, FAT2, CLCA2, CELSR1, SHC1, LGALS8, COL7A1, WNT5A, LAMA2, CDH8, NRXN3, and LY6D that were downregulated, and WWC1, PKN1, FREM2, XBP1, CLDN4, PKP2, ANG, and PFN1 which were upregulated in group 3 compared to the other subgroups (Figure 4). Gene set enrichment analysis (GSEA) confirmed the alteration of cellular binding/junction/migration signatures in group 3, including diminished function of binding and enhanced function of cell migration (Figure 4). Furthermore, GO analysis of DEGs between groups 3 and 1 showed enrichment of junctional complexes, such as ‘adherens junction’ (FDR = 0.005)’ and ‘anchoring junction’ (FDR = 0.006). Finally, we tried to validate correlation of these DEGs related to cellular adhesion and motility with the prognoses of patients with urothelial carcinoma in The Cancer Genome Atlas (TCGA) database. As a result, expression levels of CASK (p = 0.024), LAMA2 (p = 0.040), LY6D (p = 0.018), and CLDN4 (p = 0.043) were significantly associated with the overall survival of patients with urinary bladder carcinoma (Figure 5). Finally, KEGG analysis of DEGs between groups 3 and 2 showed enrichment of ‘pathways in cancer’ (FDR = 0.001) and ‘PI3K-Akt signaling pathway’ (FDR = 0.010) (Figure 3C).

2.3.2. Group 2 Showed Downregulation of Mitogen-Activated Protein Kinase (MAPK) and Tumor Necrosis Factor (TNF) Signaling Pathways

Group 2 showed enrichment of MAPK process in GO analysis, compared to group 1 (‘MAPK cascade’, FDR = 0.044) and to group 3 (‘activation of MAPK activity’, FDR = 0.009; ‘regulation of MAP kinase activity’, FDR = 0.012; ‘MAPK cascade’, FDR = 0.018; ‘regulation of MAPK cascade’, FDR = 0.040). In addition, KEGG pathway analysis revealed that DEGs between groups 1 and 2 were enriched in TNF signaling pathway (FDR = 0.014) (Figure 3A). Most DEGs associated in these pathways were downregulated in group 2 tumors (Figure 6). GSEA also indicated group 2 tumors were less reactive to MAPK and TNF signaling cascades than group 1 tumors (Figure 6).

2.3.3. Expression of Biologic Signature Genes of the Subgroups

To further delineate the properties of the IHC-based subgroups, we evaluated the expression of known biologic markers (Figure 7) [1,3,6,23,24,25]. Each subgroup was enriched with different types of keratin: KRT5, KRT6, KRT14, and KRT15 in group 1, KRT18 and KRT20 in group 3, and all in group 2. The expression patterns of basal and luminal type markers clustered moderately with group 1 and group 3, respectively. Group 2 showed modest enrichment of both subtype markers. Urothelial differentiation markers, sharing some of the gene sets with those of luminal type, showed moderate overexpression in group 3. Gene expression of late cell cycle/proliferation signature overlapped substantially with group 3 and moderately with some of the other groups, which was also differentially enriched between groups 3 and 1 (‘cell proliferation’, FDR = 0.042) and between groups 2 and 1 (‘regulation of cell proliferation’, FDR = 0.012; ‘cell proliferation’, FDR = 0.025) in GO analyses. Consistently, the expression of genes related to the progression of early urothelial carcinoma [26] was higher in group 3 and 1 than in group 2 (Mann–Whitney U, p = 0.028), showing a huge parallel to that of late cell cycle/proliferation genes, in that it was high in group 3 and in two samples of group 1. The two group 1 tumors with elevated levels of cell cycle/proliferation genes tend to express KRT14 more than the others, and the tumor with the highest KRT14 level was the only one positive for CK14 by IHC staining. Expression levels of early cell cycle, epithelial-to-mesenchymal transition, and cancer stem cell markers did not significantly overlap with the subgroups based on IHC staining. Finally, independent molecular classifiers of MIBC and NMIBC were applied to the present tumors [3,6]. Although they yielded gene expression clusters similar to their original designs, they did not overlap with the present IHC-defined subgroups (Figure S4).

3. Discussion

In the present study, we performed transcriptional investigation on non-muscle-invasive papillary high-grade UTUC by comparing its subgroups defined by CK5/6 and CK20 expression. First, we showed that immunoreactions of low CK5/6, high CK20, low CD44, and high p53 were predictive of high WHO grade in non-muscle-invasive papillary UTUC [19,20,27]. The 4 IHC markers defining the BASQ subtype of MIBC, CK14-positive, CK5/6-positive, GATA3-negative, FOXA1-negative expression [18] were not readily applicable to non-muscle-invasive papillary UTUC because positivity for CK14 and negativity for GATA3 or for FOXA1 was rarely observed in this tumor. To minimize molecular divergence between low- and high-grade urothelial carcinoma [16], we excluded low-grade tumors and focused on high-grade tumors. Along with previous studies [19,21,28], CK5/6-low/CK20-high expression was marginally correlated with the poor PFS of non-muscle-invasive papillary high-grade UTUC, which was demonstrated in the patients from the published cohort [20]. Therefore, IHC staining of CK5/6 and CK20 was adapted as prognostic markers of non-muscle-invasive papillary high-grade UTUC, which could be relevant to its molecular classifications [2,5,19,20,22]. To ensure the protein expression of the RNA-seq samples, expression levels of CK5/6 and CK20 were determined based on IHC staining on adjacent tissue with strict cut-offs [29]. Consequently, mRNA levels of basal/intermediate and luminal keratins concurred with the IHC staining profiles of CK5/6 and CK20, overall [30].
CK5/6 and CK20 expression has been suggested to represent the differentiation status of urothelium; therefore, IHC staining for these proteins can be useful as a surrogate marker pertinent to molecular subtypes of urothelial carcinoma [2,5]. CK5/6 is expressed in basal and progenitor cells of the normal urothelium which is conserved in basal type urothelial carcinoma that supposedly obtains a basal cell-like undifferentiated phenotype [31,32]. On the other hand, high levels of CK20 expression in urothelial carcinoma, as normally observed in terminally differentiated urothelium, are suggested to indicate neoplasm in well-differentiated states [31]. These were well-correlated with prognostic differences of basal and luminal type MIBC [1,2,3,5,24,25,31,33,34]. Consistent with the correlation of keratin expression and basal/luminal types of MIBC, early urothelial carcinomas showing CK5/6-high/CK20-low and CK5/6-low/CK20-high immunoreactions were enriched for basal type and luminal type/urothelial differentiation signature genes, respectively, in this study and in the previous one [1]. However, the contrasting clinical outcomes associated with CK5/6 and CK20 expression, assessed at the protein or at the mRNA levels, were frequently reported in early urothelial carcinoma of the upper or lower tract [6,19,20,21] and, occasionally, in MIBC [22]. Furthermore, previous classifiers of urothelial carcinoma based on gene expression were assigned to the present cases in a way that was independent of their IHC staining [3,6]. Collectively, these results implied that distinct molecular frameworks, that would determine their respective behaviors, were inherent in the present subgroups rather than the basal vs. luminal paradigm.
Cell adhesion and motility, in agreement with other reports [1,33], was the distinct function that was differentially enriched between groups 3 and 1 and between groups 3 and 2. Of note, group 3 showed downregulation of several genes involved in major cell adhesion complexes, including adherens junction (CDH8, CLCA2) [35,36,37], desmosome (SPINK5) [38], focal adhesion (MYO10, SHC1) [39,40], and basement membrane and the interaction with it (LGALS8, LAMA2, COL7A1, ITGB4) [41,42], with tight junction (CLDN4) being the exception [43], which would promote migration and invasion of cancer cells. As suggested by Sjodahl et al. [1], the major high-risk subtype of NMIBC, ‘genomically unstable’, had low expression levels of genes engaged in major adhesion structures, such as adherens junction, desmosome, and gap junction, with the exception being tight junction (CLDNs), compared to the low-risk subtype, ‘urothelial-like A’. The ‘genomically unstable’ subtype showed IHC staining for CK5 and CK20, similar to that of group 3 in this study [1,33]. In line with this, functional enrichment analyses and GSEA supported the premise that the loss of non-apical cell-to-cell and cell-to-matrix adhesions is a major determinant of the high-risk subtype, group 3 [37]. In addition, the non-canonical Wnt/PCP pathway (FAT2, CELSR1, WNT5A), known to suppress early-staged cancers by regulating cell adhesion or migration, was downregulated in group 3 [44]. This was also substantiated by the results showing that the alteration of the genes observed in group 3 have been reported to promote invasion of and to be associated with the poor prognoses of urothelial carcinoma and other various malignancies, including downregulation of CLCA2 [36], SPINK5 [45], SHC1 [40], LGALS8 [42], COL7A1 [41], FAT2 [44], WNT5A [44], and NRXN3 [46], and upregulation of CLDN4 [43], PKN1 [47], FREM2 [48], XBP1 [49,50], PKP2 [51], ANG [52], and PFN1 [53]. We tried to validate the prognostic values of these genes in the independent TCGA urinary bladder carcinoma cohort. However, their prognostic consequences were inconsistent, and a few of those genes altered in group 3 even predicted longer survival in the TCGA database. This result, in conjunction with those of previous studies [35,41,44,53], suggested that the adhesion molecules may act differently in advanced disease, which accounted for 99.2% of the TCGA cohort. In treatment of non-muscle-invasive UTUC, adhesion molecules differentially expressed among the present subgroups provide candidates for targeted therapy [44,54], which warrants a deeper investigation.
Late cell cycle/proliferation signature genes were predominantly upregulated in most group 3 tumors, as well as in a few samples of the other subgroups, consistent with the molecular subtypes of NMIBC with poor prognosis suggested in previous reports, ‘genomically unstable’, ‘urothelial-like B’, and ‘squamous cell carcinoma-like’ subtypes, which shared elevated levels of these signature genes compared to the ‘urothelial-like A’ subtype [1,6,24,33]. High levels of late cell cycle/proliferation markers were also observed in aggressive subtypes of MIBC [6,24]. In addition, IHC staining for cyclin B1, as a late cell cycle marker, combined with CK5 expression and tumor histology in risk-stratified patients with T1 NMIBC [33,34]. Moreover, we showed that the enrichment pattern of late cell cycle/proliferation correlated closely with that of the NMIBC progression signature [26]. Therefore, it is feasible to postulate that the enhanced cell proliferation of group 3 would support its aggressive phenotype, which was possibly obtained by bypassing cell cycle checkpoints [1]. In addition, different molecular subtypes of NMIBC that showed diverse clinical behaviors, ‘urothelial-like A’, ‘urothelial-like B’, and ‘squamous cell carcinoma-like’, as suggested by Sjodahl et al., included tumors with CK5/6-high/CK20-low IHC staining [1,33,34], which was representative of group 1. Therefore, group 1 is suggested to contain prognostically heterogeneous tumors. Considering the remarkable overexpression of these signature genes observed in certain tumors of group 1, genes related to late cell cycle/proliferation are expected to be used as risk-stratification markers within group 1 tumors. Likewise, they would also be useful for group 2 tumors that showed modest discrepancies in the level of late cell cycle/progression signature within the subgroup. This seemed analogous to the ‘urothelial-like A’ and ‘urothelial-like B’ subtypes, as both subtypes could be positive for CK5/6 and CK20, similar to group 2 tumors [1,33]. Alternatively, IHC staining for CK14, observed only in one of the group 1 samples, was one of the characteristics of the ‘squamous cell carcinoma-like’ subtype [33]; therefore, CK14 immunostaining could be used for subclassifying group 1 tumors, even though it was rarely positive in non-muscle-invasive papillary high-grade UTUC.
The accentuated patterns of IHC staining for CK5/6 and CK20 observed in group 2 were reminiscent of the immunoprofiles that were associated with low-risk phenotypes of superficial papillary urothelial neoplasms which maintained stratified CK5/6 and CK20 expression of the normal urothelium [20,55]. Similar IHC staining for CK5/6 and CK20 was also reported in the low-risk ‘urothelial-like A’ subtype which was well-differentiated at the gene expression level [1,33,34]. In transcriptional investigation, group 2 had moderate enrichment for genes involved in urothelial differentiation as well as in basal and luminal type signatures. Moreover, most DEGs enriched in MAPK or TNF pathways were expressed at lower levels in group 2 than the other groups, except for SPRED1, DUSP7, and WNT5a, that are known to downregulate MAPK and TNF signaling cascades, respectively [56,57]. MAPK activation was shown to be an important mediator of cellular transformation induced by constitutively activated FGFR3, which have crucial roles in tumorigenesis of papillary urothelial carcinoma [58]. In addition, TNF, a major inflammatory cytokine, is known to participate in many steps of carcinogenesis and tumor progression [59]. As a result, we hypothesize that group 2 represents the least cancerous subtype which maintains molecular structures of the normal urothelium, as exemplified by expression of KRTs; the expression of CK5/6 and CK20 in this group is not entirely representative of the basal and luminal types, respectively [60]. It is still not clear if this subgroup is a precursor of progressed lesions [6], for example, group 1 or group 3, or a distinct subtype with different pathogenesis.
This study had some limitations. The number of the patients was constrained by the strict cut-off criteria of IHC staining. Due to the short follow-up duration of the patients subjected to RNA-seq, their exact outcome was unevaluable, even though it was validated in the independent patient cohort. Finally, cellular adhesion, late cell cycle/proliferation signature, MAPK, or TNF signaling pathway genes, suggested as potential biological keys in the present and in previous studies, need external validation for confirmation. However, this study had the merits of investigating the transcriptional characteristics of subgroups defined by routinely used IHC staining, CK5/6 and CK20, which would facilitate the application of the accumulating genetic/phenotypic information of urothelial carcinoma, in practice.
In conclusion, IHC staining for CK5/6 and CK20 classifies subgroups of non-muscle-invasive papillary high-grade UTUC which are prognostically relevant but molecularly different from previous gene expression subtypes established in MIBC or even in NMIBC. Group 2 tumors were postulated to represent a less cancerous subtype.

4. Materials and Methods

4.1. Sample Selection and Immunohistochemistry

Radical nephroureterectomy specimens of Korean patients with UTUC were collected from Seoul National University Hospital from January 2016 to June 2018. Harvested tumors were divided in two. One half was snap-frozen and stored for RNA-seq, and the other was made into a formalin-fixed paraffin embedded (FFPE) block. Protein expression of frozen tissue was evaluated based on the IHC staining of its counterpart FFPE block. IHC staining was done on 4 μm-thick slides with the BenchMark Autostainer (Ventana, Tucson, AZ) according to the manufacturer’s instructions and using the following antibodies: CK5/6 (1:100; D5/16 B4; Dako, Glostrup, Denmark), CK20 (1:50; Ks 20.8; Dako), CK14 (1:300; LL002; Cell Marque, Rocklin, CA), GATA3 (L50-823; 1:500; Cell Marque), FOXA1 (1:500; PA5-27157; Thermo Fisher, Waltham, MA), CD44 (1:100; 156-3C11; Thermo Fisher), and p53 (1:1000; DO7; Dako). Immunoreactivity was scored semiquantitatively from 1 (<10%) to 10 (≥90%) by a 10% scale. Grading was based on the 2004 WHO system [12]. Lymph node involvement or distant metastasis was not found at diagnosis in any patients. The patients did not receive any neoadjuvant treatments before surgery. Cases showing divergent differentiation were excluded because of their potential aberrant reaction to CK20 [22]. Comparison of IHC profiles was performed by Mann–Whitney U test or Kruskal–Wallis H test. Correlation of clinicopathological variables with IHC-defined subgroups was calculated by Fisher’s exact test. SPSS 25 software (IBM, Armonk, NY, USA) was used for statistical analysis, with a two-tailed p-value < 0.05 considered significant. This study was approved by the Institutional Review Board of Seoul National University Hospital (H-1810-148-983, 07 November 2018).

4.2. Validation of Prognosis Associated with Immunohistochemical Profiles of RNA-Seq Subgroups

Prognostic validation of the IHC criteria of RNA-seq subgroups was performed by Kaplan–Meier and log-rank tests of the PFS of the patients with non-muscle-invasive papillary high-grade UTUC who were retrieved from the previously published cohort [20]. Protein expression of CK5/6 and CK20 was re-evaluated semiquantitatively in accordance with the IHC criteria of IHC-defined subgroups, high (≥ 50%) or low (<10%), from the triplicate tissue microarray IHC slides made for the previous study [20]. From 93 patients with non-muscle-invasive papillary high-grade UTUC, 39 tumors met this IHC staining criteria, including 13 CK5/6-high/CK20-low, 4 CK5/6-high/Ck20-high, and 22 CK5/6-low/CK20-high tumors. The median follow-up period was 60 months (range, 1–207). Disease progression was defined by local recurrence except for urinary bladder or distant metastasis.

4.3. Transcription Analysis and Identification of Differentially Expressed Genes

mRNA was extracted using a punch with a plunger (3 mm in diameter) from the pure papillary tumor area (tumor cellularity >90%) of frozen tissue that was marked on a hematoxylin and eosin-stained frozen slide. RNA-seq was performed as previously described [61] using the Illumina HiSeq 2500 platform (Illumina, San Diego, CA, USA). Paired demultiplexed fastq files were generated, and initial quality control was performed using FastQC (Phred quality score >30) [62]. Adapter trimming was conducted using the Trimmomatic tool [63] followed by mapping to the human genome (UCSC hg19) using the HISAT2 and Bowtie2 [64]. Known transcripts were assembled with the StringTie tool [65]. Consequently, 21,780 genes were found across the samples. Quantile normalization was done on log2 (FPKM+1). To fit a linear model for each gene based on the three-group study design, the function of lmFit was used (limma) [66]. Then, the eBayes function was used to calculate the empirical Bayes moderated t-statistic.

4.4. Functional Enrichment Analysis

Functional enrichment of DEGs were found from GO [67] and KEGG pathway [68] database with FDR < 0.05 as a cut-off. The GSEA of all genes was used to interpret their functions from the public database [69]. Biological signatures were acquired from the respective publications [1,3,6,23,24,25]. The integrated score for the progression of NMIBC was calculated by subtracting the total value of COL4A3BP, MBNL2, NEK1, FABP4, and SKAP2 from that of KPNA2, BIRC5, UB2C, CDC25B, COL4A1, MSN, and COL18A1 [26].

4.5. TCGA Data Preprocessing and Survival Analysis

From the TCGA database, 399 samples of urinary bladder urothelial carcinoma were collected that contained mRNA expression and survival data. Most tumors were WHO high-grade (n = 375) and MIBC (n = 396). The median follow-up period was 1735 days, during which 109 patients died. Preprocessing of the mRNA expression data was performed in the R 3.4.3 statistical environment: mRNA expression level quantified by RSEM raw counts were normalized using the limma-voom, which applied linear modeling to voom-transformed read counts [70]. We determined the median value of gene expression as a cut-off point to divide patients into high- and low-expression groups. Kaplan–Meier and log-rank tests were applied for survival analysis.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/3/570/s1.

Author Contributions

Conceptualization, M.J., J.H.P, and K.C.M.; methodology, M.J. and J.H.L.; formal analysis, M.J. and J.H.L.; investigation, M.J. and K.C.M.; resources, M.J and B.K.; writing—original draft preparation, M.J. and J.H.L.; writing—review and editing, M.J. and K.C.M.; funding acquisition, K.C.M.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (grant number 2018R1D1A1B07045763).

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.

Abbreviations

MIBCmuscle-invasive bladder cancer
NMIBCnon-muscle-invasive bladder cancer
CIScarcinoma in situ
UTUCupper tract urothelial carcinoma
IHCimmunohistochemistry
BASQbasal/squamous-like
RNA-seqRNA sequencing
PFSprogression-free survival
DEGdifferentially expressed gene
FDRfalse discovery rate
GOgene ontology
KEGGKyoto Encyclopedia of Genes and Genomes
BPbiologic process
CCcellular component
MFmolecular function
GSEAgene set enrichment analysis
TCGAThe Cancer Genome Atlas
MAPKmitogen-activated protein kinase
TNFtumor necrosis factor

References

  1. Sjodahl, G.; Lauss, M.; Lovgren, K.; Chebil, G.; Gudjonsson, S.; Veerla, S.; Patschan, O.; Aine, M.; Ferno, M.; Ringner, M.; et al. A molecular taxonomy for urothelial carcinoma. Clin. Cancer. Res. 2012, 18, 3377–3386. [Google Scholar] [CrossRef] [PubMed]
  2. Choi, W.; Porten, S.; Kim, S.; Willis, D.; Plimack, E.R.; Hoffman-Censits, J.; Roth, B.; Cheng, T.; Tran, M.; Lee, I.L.; et al. Identification of distinct basal and luminal subtypes of muscle-invasive bladder cancer with different sensitivities to frontline chemotherapy. Cancer Cell 2014, 25, 152–165. [Google Scholar] [CrossRef] [PubMed]
  3. Damrauer, J.S.; Hoadley, K.A.; Chism, D.D.; Fan, C.; Tiganelli, C.J.; Wobker, S.E.; Yeh, J.J.; Milowsky, M.I.; Iyer, G.; Parker, J.S.; et al. Intrinsic subtypes of high-grade bladder cancer reflect the hallmarks of breast cancer biology. Proc. Natl. Acad. Sci. USA 2014, 111, 3110–3115. [Google Scholar] [CrossRef] [Green Version]
  4. The Cancer Genome Atlas Research Network. Comprehensive molecular characterization of urothelial bladder carcinoma. Nature 2014, 507, 315–322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Robertson, A.G.; Kim, J.; Al-Ahmadie, H.; Bellmunt, J.; Guo, G.; Cherniack, A.D.; Hinoue, T.; Laird, P.W.; Hoadley, K.A.; Akbani, R.; et al. Comprehensive molecular characterization of muscle-invasive bladder cancer. Cell 2017, 171, 540–556.e25. [Google Scholar] [CrossRef] [PubMed]
  6. Hedegaard, J.; Lamy, P.; Nordentoft, I.; Algaba, F.; Hoyer, S.; Ulhoi, B.P.; Vang, S.; Reinert, T.; Hermann, G.G.; Mogensen, K.; et al. Comprehensive transcriptional analysis of early-stage urothelial carcinoma. Cancer Cell 2016, 30, 27–42. [Google Scholar] [CrossRef] [PubMed]
  7. Aine, M.; Eriksson, P.; Liedberg, F.; Sjodahl, G.; Hoglund, M. Biological determinants of bladder cancer gene expression subtypes. Sci. Rep. 2015, 5, 10957. [Google Scholar] [CrossRef] [Green Version]
  8. Seiler, R.; Ashab, H.A.D.; Erho, N.; van Rhijn, B.W.G.; Winters, B.; Douglas, J.; Van Kessel, K.E.; Fransen van de Putte, E.E.; Sommerlad, M.; Wang, N.Q.; et al. Impact of molecular subtypes in muscle-invasive bladder cancer on predicting response and survival after neoadjuvant chemotherapy. Eur. Urol. 2017, 72, 544–554. [Google Scholar] [CrossRef]
  9. McConkey, D.J.; Choi, W.; Ochoa, A.; Siefker-Radtke, A.; Czerniak, B.; Dinney, C.P. Therapeutic opportunities in the intrinsic subtypes of muscle-invasive bladder cancer. Hematol. Oncol. Clin. N. Am. 2015, 29, 377–394. [Google Scholar] [CrossRef]
  10. McConkey, D.J.; Choi, W.; Shen, Y.; Lee, I.L.; Porten, S.; Matin, S.F.; Kamat, A.M.; Corn, P.; Millikan, R.E.; Dinney, C.; et al. A prognostic gene expression signature in the molecular classification of chemotherapy-naive urothelial cancer is predictive of clinical outcomes from neoadjuvant chemotherapy: A phase 2 trial of dose-dense methotrexate, vinblastine, doxorubicin, and cisplatin with bevacizumab in urothelial cancer. Eur. Urol. 2016, 69, 855–862. [Google Scholar] [CrossRef]
  11. Hurst, C.D.; Alder, O.; Platt, F.M.; Droop, A.; Stead, L.F.; Burns, J.E.; Burghel, G.J.; Jain, S.; Klimczak, L.J.; Lindsay, H.; et al. Genomic subtypes of non-invasive bladder cancer with distinct metabolic profile and female gender bias in KDM6A mutation frequency. Cancer Cell 2017, 32, 701–715.e07. [Google Scholar] [CrossRef] [PubMed]
  12. Roupret, M.; Babjuk, M.; Comperat, E.; Zigeuner, R.; Sylvester, R.J.; Burger, M.; Cowan, N.C.; Gontero, P.; Van Rhijn, B.W.G.; Mostafid, A.H.; et al. European association of urology guidelines on upper urinary tract urothelial carcinoma: 2017 update. Eur. Urol. 2018, 73, 111–122. [Google Scholar] [CrossRef] [PubMed]
  13. Sfakianos, J.P.; Cha, E.K.; Iyer, G.; Scott, S.N.; Zabor, E.C.; Shah, R.H.; Ren, Q.; Bagrodia, A.; Kim, P.H.; Hakimi, A.A.; et al. Genomic characterization of upper tract urothelial carcinoma. Eur. Urol. 2015, 68, 970–977. [Google Scholar] [CrossRef]
  14. Sanford, T.; Porten, S.; Meng, M.V. Molecular analysis of upper tract and bladder urothelial carcinoma: Results from a microarray comparison. PLoS ONE 2015, 10, e0137141. [Google Scholar] [CrossRef] [PubMed]
  15. Glaser, A.P.; Fantini, D.; Shilatifard, A.; Schaeffer, E.M.; Meeks, J.J. The evolving genomic landscape of urothelial carcinoma. Nat. Rev. Urol. 2017, 14, 215–229. [Google Scholar] [CrossRef] [PubMed]
  16. Audenet, F.; Attalla, K.; Sfakianos, J.P. The evolution of bladder cancer genomics: What have we learned and how can we use it? Urol. Oncol. 2018, 36, 313–320. [Google Scholar] [CrossRef] [PubMed]
  17. Dadhania, V.; Zhang, M.; Zhang, L.; Bondaruk, J.; Majewski, T.; Siefker-Radtke, A.; Guo, C.C.; Dinney, C.; Cogdell, D.E.; Zhang, S.; et al. Meta-analysis of the luminal and basal subtypes of bladder cancer and the identification of signature immunohistochemical markers for clinical use. EBioMedicine 2016, 12, 105–117. [Google Scholar] [CrossRef]
  18. Lerner, S.P.; McConkey, D.J.; Hoadley, K.A.; Chan, K.S.; Kim, W.Y.; Radvanyi, F.; Hoglund, M.; Real, F.X. Bladder cancer molecular taxonomy: Summary from a consensus meeting. Bladder Cancer 2016, 2, 37–47. [Google Scholar] [CrossRef]
  19. Sikic, D.; Keck, B.; Wach, S.; Taubert, H.; Wullich, B.; Goebell, P.J.; Kahlmeyer, A.; Olbert, P.; Isfort, P.; Nimphius, W.; et al. Immunohistochemiocal subtyping using CK20 and CK5 can identify urothelial carcinomas of the upper urinary tract with a poor prognosis. PLoS ONE 2017, 12, e0179602. [Google Scholar] [CrossRef]
  20. Jung, M.; Kim, B.; Moon, K.C. Immunohistochemistry of CK5/6, CD44, and CK20 as prognostic biomarkers of non-muscle-invasive papillary upper tract urothelial carcinoma. Histopathology 2019, 74, 483–493. [Google Scholar] [CrossRef]
  21. Breyer, J.; Wirtz, R.M.; Otto, W.; Erben, P.; Kriegmair, M.C.; Stoehr, R.; Eckstein, M.; Eidt, S.; Denzinger, S.; Burger, M.; et al. In stage pT1 non-muscle-invasive bladder cancer (NMIBC), high KRT20 and low KRT5 mRNA expression identify the luminal subtype and predict recurrence and survival. Virchows Arch. 2017, 470, 267–274. [Google Scholar] [CrossRef] [PubMed]
  22. Eckstein, M.; Wirtz, R.M.; Gross-Weege, M.; Breyer, J.; Otto, W.; Stoehr, R.; Sikic, D.; Keck, B.; Eidt, S.; Burger, M.; et al. mRNA-expression of KRT5 and KRT20 defines distinct prognostic subgroups of muscle-invasive urothelial bladder cancer correlating with histological variants. Int. J. Mol. Sci. 2018, 19, 3396. [Google Scholar] [CrossRef] [PubMed]
  23. Biton, A.; Bernard-Pierrot, I.; Lou, Y.; Krucker, C.; Chapeaublanc, E.; Rubio-Perez, C.; Lopez-Bigas, N.; Kamoun, A.; Neuzillet, Y.; Gestraud, P.; et al. Independent component analysis uncovers the landscape of the bladder tumor transcriptome and reveals insights into luminal and basal subtypes. Cell Rep. 2014, 9, 1235–1245. [Google Scholar] [CrossRef]
  24. Tan, T.Z.; Rouanne, M.; Tan, K.T.; Huang, R.Y.; Thiery, J.P. Molecular subtypes of urothelial bladder cancer: Results from a meta-cohort analysis of 2411 tumors. Eur. Urol. 2018. [Google Scholar] [CrossRef]
  25. Rebouissou, S.; Bernard-Pierrot, I.; de Reynies, A.; Lepage, M.L.; Krucker, C.; Chapeaublanc, E.; Herault, A.; Kamoun, A.; Caillault, A.; Letouze, E.; et al. EGFR as a potential therapeutic target for a subset of muscle-invasive bladder cancers presenting a basal-like phenotype. Sci. Transl. Med. 2014, 6, 244ra91. [Google Scholar] [CrossRef] [PubMed]
  26. Dyrskjot, L.; Reinert, T.; Novoradovsky, A.; Zuiverloon, T.C.; Beukers, W.; Zwarthoff, E.; Malats, N.; Real, F.X.; Segersten, U.; Malmstrom, P.U.; et al. Analysis of molecular intra-patient variation and delineation of a prognostic 12-gene signature in non-muscle invasive bladder cancer; technology transfer from microarrays to PCR. Br. J. Cancer 2012, 107, 1392–1398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kelsey, K.T.; Hirao, T.; Schned, A.; Hirao, S.; Devi-Ashok, T.; Nelson, H.H.; Andrew, A.; Karagas, M.R. A population-based study of immunohistochemical detection of p53 alteration in bladder cancer. Br. J. Cancer 2004, 90, 1572–1576. [Google Scholar] [CrossRef] [Green Version]
  28. Bertz, S.; Otto, W.; Denzinger, S.; Wieland, W.F.; Burger, M.; Stohr, R.; Link, S.; Hofstadter, F.; Hartmann, A. Combination of CK20 and Ki-67 immunostaining analysis predicts recurrence, progression, and cancer-specific survival in pT1 urothelial bladder cancer. Eur. Urol. 2014, 65, 218–226. [Google Scholar] [CrossRef]
  29. Sjodahl, G.; Eriksson, P.; Liedberg, F.; Hoglund, M. Molecular classification of urothelial carcinoma: Global mRNA classification versus tumour-cell phenotype classification. J. Pathol. 2017, 242, 113–125. [Google Scholar] [CrossRef]
  30. Schaafsma, H.E.; Ramaekers, F.C.; van Muijen, G.N.; Lane, E.B.; Leigh, I.M.; Robben, H.; Huijsmans, A.; Ooms, E.C.; Ruiter, D.J. Distribution of cytokeratin polypeptides in human transitional cell carcinomas, with special emphasis on changing expression patterns during tumor progression. Am. J. Pathol. 1990, 136, 329–343. [Google Scholar]
  31. Volkmer, J.P.; Sahoo, D.; Chin, R.K.; Ho, P.L.; Tang, C.; Kurtova, A.V.; Willingham, S.B.; Pazhanisamy, S.K.; Contreras-Trujillo, H.; Storm, T.A.; et al. Three differentiation states risk-stratify bladder cancer into distinct subtypes. Proc. Natl. Acad. Sci. USA 2012, 109, 2078–2083. [Google Scholar] [CrossRef] [Green Version]
  32. Ho, P.L.; Kurtova, A.; Chan, K.S. Normal and neoplastic urothelial stem cells: Getting to the root of the problem. Nat. Rev. Urol. 2012, 9, 583–594. [Google Scholar] [CrossRef] [PubMed]
  33. Sjodahl, G.; Lovgren, K.; Lauss, M.; Patschan, O.; Gudjonsson, S.; Chebil, G.; Aine, M.; Eriksson, P.; Mansson, W.; Lindgren, D.; et al. Toward a molecular pathologic classification of urothelial carcinoma. Am. J. Pathol. 2013, 183, 681–691. [Google Scholar] [CrossRef] [PubMed]
  34. Patschan, O.; Sjodahl, G.; Chebil, G.; Lovgren, K.; Lauss, M.; Gudjonsson, S.; Kollberg, P.; Eriksson, P.; Aine, M.; Mansson, W.; et al. A molecular pathologic framework for risk stratification of stage T1 urothelial carcinoma. Eur. Urol. 2015, 68, 824–832. [Google Scholar] [CrossRef] [PubMed]
  35. Bryan, R.T. Cell adhesion and urothelial bladder cancer: The role of cadherin switching and related phenomena. Philos. Trans. R. Soc. B Biol. Sci. 2015, 370, 20140042. [Google Scholar] [CrossRef] [PubMed]
  36. Ramena, G.; Yin, Y.; Yu, Y.; Walia, V.; Elble, R.C. CLCA2 interactor EVA1 is required for mammary epithelial cell differentiation. PLoS ONE 2016, 11, e0147489. [Google Scholar] [CrossRef]
  37. Koivunen, J.; Aaltonen, V.; Koskela, S.; Lehenkari, P.; Laato, M.; Peltonen, J. Protein kinase C alpha/beta inhibitor Go6976 promotes formation of cell junctions and inhibits invasion of urinary bladder carcinoma cells. Cancer Res. 2004, 64, 5693–5701. [Google Scholar] [CrossRef]
  38. Ishida-Yamamoto, A.; Igawa, S.; Kishibe, M.; Honma, M. Clinical and molecular implications of structural changes to desmosomes and corneodesmosomes. J. Dermatol. 2018, 45, 385–389. [Google Scholar] [CrossRef]
  39. Jacquemet, G.; Hamidi, H.; Ivaska, J. Filopodia in cell adhesion, 3D migration and cancer cell invasion. Curr. Opin. Cell Biol. 2015, 36, 23–31. [Google Scholar] [CrossRef] [Green Version]
  40. Zheng, Z.; Luan, X.; Zha, J.; Li, Z.; Wu, L.; Yan, Y.; Wang, H.; Hou, D.; Huang, L.; Huang, F.; et al. TNF-alpha inhibits the migration of oral squamous cancer cells mediated by miR-765-EMP3-p66Shc axis. Cell Signal 2017, 34, 102–109. [Google Scholar] [CrossRef] [PubMed]
  41. Brunner, A.; Tzankov, A. The role of structural extracellular matrix proteins in urothelial bladder cancer (review). Biomark. Insights 2007, 2, 418–427. [Google Scholar] [CrossRef] [PubMed]
  42. Kramer, M.W.; Waalkes, S.; Serth, J.; Hennenlotter, J.; Tezval, H.; Stenzl, A.; Kuczyk, M.A.; Merseburger, A.S. Decreased galectin-8 is a strong marker for recurrence in urothelial carcinoma of the bladder. Urol. Int. 2011, 87, 143–150. [Google Scholar] [CrossRef]
  43. Szekely, E.; Torzsok, P.; Riesz, P.; Korompay, A.; Fintha, A.; Szekely, T.; Lotz, G.; Nyirady, P.; Romics, I.; Timar, J.; et al. Expression of claudins and their prognostic significance in noninvasive urothelial neoplasms of the human urinary bladder. J. Histochem. Cytochem. 2011, 59, 932–941. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, Y. Wnt/Planar cell polarity signaling: A new paradigm for cancer therapy. Mol. Cancer Ther. 2009, 8, 2103–2109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhao, C.; Zou, H.; Zhang, J.; Wang, J.; Liu, H. An integrated methylation and gene expression microarray analysis reveals significant prognostic biomarkers in oral squamous cell carcinoma. Oncol. Rep. 2018, 40, 2637–2647. [Google Scholar] [CrossRef]
  46. Sun, H.T.; Cheng, S.X.; Tu, Y.; Li, X.H.; Zhang, S. FoxQ1 promotes glioma cells proliferation and migration by regulating NRXN3 expression. PLoS ONE 2013, 8, e55693. [Google Scholar] [CrossRef] [PubMed]
  47. O’Sullivan, A.G.; Mulvaney, E.P.; Kinsella, B.T. Regulation of protein kinase C-related kinase (PRK) signalling by the TPalpha and TPbeta isoforms of the human thromboxane A2 receptor: Implications for thromboxane- and androgen- dependent neoplastic and epigenetic responses in prostate cancer. Biochim. Biophys. Acta Mol. Basis Dis. 2017, 1863, 838–856. [Google Scholar] [CrossRef]
  48. Nagaishi, M.; Kim, Y.H.; Mittelbronn, M.; Giangaspero, F.; Paulus, W.; Brokinkel, B.; Vital, A.; Tanaka, Y.; Nakazato, Y.; Legras-Lachuer, C.; et al. Amplification of the STOML3, FREM2, and LHFP genes is associated with mesenchymal differentiation in gliosarcoma. Am. J. Pathol. 2012, 180, 1816–1823. [Google Scholar] [CrossRef]
  49. Fang, D.; Kitamura, H. Cancer stem cells and epithelial-mesenchymal transition in urothelial carcinoma: Possible pathways and potential therapeutic approaches. Int. J. Urol. 2018, 25, 7–17. [Google Scholar] [CrossRef]
  50. Sun, Y.; Jiang, F.; Pan, Y.; Chen, X.; Chen, J.; Wang, Y.; Zheng, X.; Zhang, J. XBP1 promotes tumor invasion and is associated with poor prognosis in oral squamous cell carcinoma. Oncol. Rep. 2018, 40, 988–998. [Google Scholar] [CrossRef]
  51. Takahashi, H.; Nakatsuji, H.; Takahashi, M.; Avirmed, S.; Fukawa, T.; Takemura, M.; Fukumori, T.; Kanayama, H. Up-regulation of plakophilin-2 and Down-regulation of plakophilin-3 are correlated with invasiveness in bladder cancer. Urology 2012, 79, 240.e1–e8. [Google Scholar] [CrossRef]
  52. Miyake, H.; Hara, I.; Yamanaka, K.; Gohji, K.; Arakawa, S.; Kamidono, S. Increased angiogenin expression in the tumor tissue and serum of urothelial carcinoma patients is related to disease progression and recurrence. Cancer 1999, 86, 316–324. [Google Scholar] [CrossRef] [Green Version]
  53. Frantzi, M.; Klimou, Z.; Makridakis, M.; Zoidakis, J.; Latosinska, A.; Borras, D.M.; Janssen, B.; Giannopoulou, I.; Lygirou, V.; Lazaris, A.C.; et al. Silencing of Profilin-1 suppresses cell adhesion and tumor growth via predicted alterations in integrin and Ca2+ signaling in T24M-based bladder cancer models. Oncotarget 2016, 7, 70750–70768. [Google Scholar] [CrossRef] [Green Version]
  54. Kowalski, M.; Entwistle, J.; Cizeau, J.; Niforos, D.; Loewen, S.; Chapman, W.; MacDonald, G.C. A Phase I study of an intravesically administered immunotoxin targeting EpCAM for the treatment of nonmuscle-invasive bladder cancer in BCGrefractory and BCG-intolerant patients. Drug Des. Dev. Ther. 2010, 4, 313–320. [Google Scholar] [CrossRef]
  55. Desai, S.; Lim, S.D.; Jimenez, R.E.; Chun, T.; Keane, T.E.; McKenney, J.K.; Zavala-Pompa, A.; Cohen, C.; Young, R.H.; Amin, M.B. Relationship of cytokeratin 20 and CD44 protein expression with WHO/ISUP grade in pTa and pT1 papillary urothelial neoplasia. Mod. Pathol. 2000, 13, 1315–1323. [Google Scholar] [CrossRef]
  56. Delire, B.; Starkel, P. The Ras/MAPK pathway and hepatocarcinoma: Pathogenesis and therapeutic implications. Eur. J. Clin. Investig. 2015, 45, 609–623. [Google Scholar] [CrossRef]
  57. Kim, J.H.; Park, S.; Chung, H.; Oh, S. Wnt5a attenuates the pathogenic effects of the Wnt/beta-catenin pathway in human retinal pigment epithelial cells via down-regulating beta-catenin and Snail. BMB Rep. 2015, 48, 525–530. [Google Scholar] [CrossRef]
  58. Wu, X.R. Urothelial tumorigenesis: A tale of divergent pathways. Nat. Rev. Cancer 2005, 5, 713–725. [Google Scholar] [CrossRef]
  59. Sethi, G.; Sung, B.; Aggarwal, B.B. TNF: A master switch for inflammation to cancer. Front. Biosci. 2008, 13, 5094–5107. [Google Scholar] [CrossRef]
  60. Van Batavia, J.; Yamany, T.; Molotkov, A.; Dan, H.; Mansukhani, M.; Batourina, E.; Schneider, K.; Oyon, D.; Dunlop, M.; Wu, X.R.; et al. Bladder cancers arise from distinct urothelial sub-populations. Nat. Cell Biol. 2014, 16, 982–991. [Google Scholar] [CrossRef]
  61. Martin, J.A.; Wang, Z. Next-generation transcriptome assembly. Nat. Rev. Genet. 2011, 12, 671–682. [Google Scholar] [CrossRef] [Green Version]
  62. Babraham Bioinformatics. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 18 December 2018).
  63. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed]
  65. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  67. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  68. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef]
  69. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  70. Law, C.W.; Chen, Y.; Shi, W.; Smyth, G.K. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. 2014, 15, R29. [Google Scholar] [CrossRef]
Ijms 20 00570 g001a 550Ijms 20 00570 g001b 550
Figure 1. IHC staining for CK5/6 and CK20 in non-muscle-invasive papillary upper tract urothelial carcinoma (UTUC). (A) CK5/6-low and CK20-high IHC staining is significantly related to high WHO grade of non-muscle-invasive papillary UTUC. Blue bars indicate the mean value ± standard deviation. (B) Representative images of IHC staining for CK5/6 and CK20 across subgroups of non-muscle-invasive papillary high-grade UTUC. Original magnification × 40.
Figure 1. IHC staining for CK5/6 and CK20 in non-muscle-invasive papillary upper tract urothelial carcinoma (UTUC). (A) CK5/6-low and CK20-high IHC staining is significantly related to high WHO grade of non-muscle-invasive papillary UTUC. Blue bars indicate the mean value ± standard deviation. (B) Representative images of IHC staining for CK5/6 and CK20 across subgroups of non-muscle-invasive papillary high-grade UTUC. Original magnification × 40.
Ijms 20 00570 g001aIjms 20 00570 g001b
Ijms 20 00570 g002 550
Figure 2. Identification of differentially expressed genes (DEGs) between the subgroups of non-muscle-invasive papillary high-grade UTUC. (A) Volcano plots of the DEGs between each subgroup. Genes that meet the threshold are in blue. (B) Venn diagram of all comparisons demonstrates 308 DEGs across the subgroups. (C) Unsupervised hierarchical clustering of all DEGs for all samples.
Figure 2. Identification of differentially expressed genes (DEGs) between the subgroups of non-muscle-invasive papillary high-grade UTUC. (A) Volcano plots of the DEGs between each subgroup. Genes that meet the threshold are in blue. (B) Venn diagram of all comparisons demonstrates 308 DEGs across the subgroups. (C) Unsupervised hierarchical clustering of all DEGs for all samples.
Ijms 20 00570 g002
Ijms 20 00570 g003a 550Ijms 20 00570 g003b 550
Figure 3. Top 10 significantly enriched biologic process (BP), cellular component (CC), and molecular function (MF) terms of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of each comparison. (A) Group 1 vs. Group 2, (B) Group 1 vs. Group 3, (C) Group 2 vs. Group 3. FDR, false discovery rate.
Figure 3. Top 10 significantly enriched biologic process (BP), cellular component (CC), and molecular function (MF) terms of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of each comparison. (A) Group 1 vs. Group 2, (B) Group 1 vs. Group 3, (C) Group 2 vs. Group 3. FDR, false discovery rate.
Ijms 20 00570 g003aIjms 20 00570 g003b
Ijms 20 00570 g004 550
Figure 4. Heatmaps of DEGs enriched in cellular adhesion/motility processes and related GSEA results. Group 3 was enriched with biologic themes of cellular adhesion and motility. ES, enrichment score.
Figure 4. Heatmaps of DEGs enriched in cellular adhesion/motility processes and related GSEA results. Group 3 was enriched with biologic themes of cellular adhesion and motility. ES, enrichment score.
Ijms 20 00570 g004
Ijms 20 00570 g005 550
Figure 5. Survival analyses of DEGs involved in cellular adhesion/motility in The Cancer Genome Atlas (TCGA) database. Four genes were significantly associated with the overall survival of the patients with urinary bladder carcinoma.
Figure 5. Survival analyses of DEGs involved in cellular adhesion/motility in The Cancer Genome Atlas (TCGA) database. Four genes were significantly associated with the overall survival of the patients with urinary bladder carcinoma.
Ijms 20 00570 g005
Ijms 20 00570 g006 550
Figure 6. Heatmaps of DEGs enriched in MAPK or TNF signaling pathways, along with related GSEA results. Group 2 was characterized by low levels of genes involved in MAPK or TNF signaling pathways, except for a few genes that play regulatory roles in these pathways. ES, enrichment score.
Figure 6. Heatmaps of DEGs enriched in MAPK or TNF signaling pathways, along with related GSEA results. Group 2 was characterized by low levels of genes involved in MAPK or TNF signaling pathways, except for a few genes that play regulatory roles in these pathways. ES, enrichment score.
Ijms 20 00570 g006
Ijms 20 00570 g007 550
Figure 7. Heatmaps of known biological signature genes extracted from previous reports, with samples arranged in each column and subgroups divided by black vertical lines. Progression (right bottom) was visualized according to the integrated score (red, high; blue, low). Uro differentiation, urothelial differentiation.
Figure 7. Heatmaps of known biological signature genes extracted from previous reports, with samples arranged in each column and subgroups divided by black vertical lines. Progression (right bottom) was visualized according to the integrated score (red, high; blue, low). Uro differentiation, urothelial differentiation.
Ijms 20 00570 g007
Table
Table 1. Clinicopathological characteristics of subgroups of non-muscle-invasive papillary high-grade UTUC and their IHC expression.
Table 1. Clinicopathological characteristics of subgroups of non-muscle-invasive papillary high-grade UTUC and their IHC expression.
VariablesGroup 1Group 2Group 3p-value
Number (%)5 (100)5 (100)5 (100)
Age 1.000
≥693 (60)3 (60)3 (60)
<692 (40)2 (40)2 (40)
Sex 1.000
Male3 (60)4 (80)3 (60)
Female2 (40)1 (20)2 (40)
Size 0.725
≥3.62 (40)1 (20)0 (0)
<3.63 (60)4 (80)5 (100)
Organ 0.301
Pelvis/calyx4 (80)1 (20)3 (60)
Ureter1 (20)4 (80)2 (40)
Stage 0.800
T14 (80)2 (40)3 (60)
Ta1 (20)3 (60)2 (40)
CIS 1.000
Present2 (40)1 (20)2 (40)
Absent3 (60)4 (80)3 (60)
IHC score (mean ± SD) 1
CK5/66.4 ± 2.417.4 ± 2.411 ± 0.000.001
CK201.0 ± 0.008.4 ± 1.5210 ± 0.00<0.001
CD446.0 ± 2.835.4 ± 3.213.0 ± 2.550.162
p537.5 ± 3.706.3 ± 2.224.8 ± 3.550.428
CK141.4 ± 0.891.0 ± 0.001.0 ± 0.001.000
GATA38.0 ± 3.9410.0 ± 0.0010.0 ± 0.000.286
FOXA16.4 ± 4.5110.0 ± 0.008.4 ± 1.950.066
1 Semiquantitative scores from 1 (<10%) to 10 (90%–100%) by a 10% scale. CIS, carcinoma in situ; IHC, immunohistochemistry.

Share and Cite

MDPI and ACS Style

Jung, M.; Lee, J.H.; Kim, B.; Park, J.H.; Moon, K.C. Transcriptional Analysis of Immunohistochemically Defined Subgroups of Non-Muscle-Invasive Papillary High-Grade Upper Tract Urothelial Carcinoma. Int. J. Mol. Sci. 2019, 20, 570. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030570

AMA Style

Jung M, Lee JH, Kim B, Park JH, Moon KC. Transcriptional Analysis of Immunohistochemically Defined Subgroups of Non-Muscle-Invasive Papillary High-Grade Upper Tract Urothelial Carcinoma. International Journal of Molecular Sciences. 2019; 20(3):570. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030570

Chicago/Turabian Style

Jung, Minsun, Jeong Hoon Lee, Bohyun Kim, Jeong Hwan Park, and Kyung Chul Moon. 2019. "Transcriptional Analysis of Immunohistochemically Defined Subgroups of Non-Muscle-Invasive Papillary High-Grade Upper Tract Urothelial Carcinoma" International Journal of Molecular Sciences 20, no. 3: 570. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030570

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