Abstract
Keywords
ABBREVIATIONS
Cytokeratin Array comparative genomic hybridization Urothelial carcinoma Fetal bovine serum Normal human urothelial cells 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide Senescence-associated, acidic β-galactosidase Quantitative polymerase chain reaction Human papilloma virus
INTRODUCTION
Cell culture models of normal urothelial cells are valuable for studying urothelial physiology and differentiation, for investigating the function of genetic and epigenetic changes found in urothelial carcinoma and as controls in the evaluation of the tumor specificity of novel cancer drugs. Primary cultures of urothelial cells, usually established from healthy ureters removed during tumor nephrectomy, provide one valuable model [1]. Importantly, these cultures can be induced to differentiate into urothelium-like epithelia by several protocols [2, 3]. However, they have a finite lifespan, display large interindividual variation and are not broadly available. Thus, several immortalized cell lines have been generated from such cultures by introducing viral oncogenes (SV-HUC1, UROtsa) or telomerase (TERT-NHUC) [4, 5]. Whereas viral oncogenes like the SV40 largeT-antigen elicit partial transformation by interfering with cell cycle control and checkpoints by inactivating RB1 and p53, urothelial cell lines immortalized by telomerase retain proper cell cycle control and are generally genetically stable, but, over time, they lose the ability to differentiate, which is associated with changes in the expression of important chromatin regulator proteins [6, 7].
We obtained the cell line HBLAK from a company supplying primary cultures of normal urothelial cells, who noted that one batch of these cells did not senesce as usual after several passages. Having confirmed this observation, we characterized this cell line described as “spontaneously immortalized” by the distributor in detail, investigating its expression of urothelial markers, its ability to differentiate, its karyotype and typical alterations associated with urothelial transformation, as well as its suitability as a normal control when investigating new therapeutic compounds. Although this cell line originates from ostensibly normal bladder tissue and shares many characteristics with normal urothelial cells, we discovered some genetic changes thought to contribute to immortalization of urothelial cells. Taken together, our results suggest that HBLAK cells retain many phenotypic characteristics of normal urothelial cells, but also harbor a limited number of genetic changes known to contribute to immortalization of normal and cancerous urothelial cells.
MATERIALS AND METHODS
Cell culture and induction of differentiation
The human urothelial cell line HBLAK, spontaneously outgrown from a standard lot of normal human uroepithelial cells, was kindly provided by CELLnTEC (Switzerland), from where the cell line is commercially available. This culture was established from morphologically inconspicuous primary urothelium taken from the urinary bladder of an 80 years old male patient who underwent surgery for a disease other than genitourinary cancer. The patient was not diagnosed for bladder cancer at the time of the surgery. Due to ethical guidelines no further follow-up data is available for the donor. HBLAK cells were cultured according to the manufacturer’s recommendation in serum-free CnT-Prime Epithelial Culture Medium (CELLnTEC) at 37°C and with 5% CO2.
To study their differentiation capacity the cells were seeded subconfluently into 6-well plates and subsequently treated according to various published protocols [2, 3]. In one protocol, EGFR inhibition and PPARγ activation were achieved by treatment with 1 μM PD153035 (Merck, Germany) and 1 μM Troglitazone (Cayman, USA) in CnT-Prime Epithelial Culture Medium (ready to use and supplemented with EGF, CELLnTEC) or in KFSM medium without supplemented EGF (Life Technologies, Germany) for seven days, during which the medium was changed twice. For some experiments 5% FBS (Biochrom, Germany) was added to cells in CnT-Prime Medium. For a second, calcium-based protocol near-confluent cultures were maintained in CnT-Prime Medium supplemented with CaCl2 to a final concentration of 2 mM, with additional 5% FBS in some experiments, for a period of 10–14 days. The medium was changed every 3 days.
For comparison, primary cultures of normal urothelial cells (NHUC) were established from healthy ureters removed during tumor nephrectomy. These samples were collected with informed consent of the patients and their use was approved by the Ethics Committee of the Medical Faculty of the Heinrich-Heine-University, study number 1788. The cells were cultured as published previously [8] in KFSM medium supplemented with 5 ng/ml EGF and 50 μg/ml bovine pituitary extract (Life Technologies).
TERT-NHUC cells, kindly provided by Prof. M.A. Knowles, University of Leeds, UK, were cultured as previously described [9].
Cell proliferation analysis
For calculating cell doubling time, cell viability was determined over a period of four days by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide dye reduction assay (MTT, Sigma Aldrich, Germany). IC50 concentrations for Cisplatin (Accord Healthcare, Germany) were determined by treatment of different passages of HBLAK, different passages of primary NHUC cells and triplicate analysis of TERT-NHUC cells in their respective standard media at various concentrations for 72 h.
Senescent cells were identified by acidic β-galactosidase staining. Following fixation for 5 minutes in 2% formaldehyde and 0.2% glutaraldehyde and washing, cells were stained with fresh SA-β-Gal staining solution (1 mg/ml X-Gal, 150 mM NaCl, 2 mM MgCl2, 5 mM K3Fe(CN)6, 5 mM K4Fe(CN)6) for 4 h at 37°C. Images were taken using the NIS-Elements software with a Nikon Eclipse TE2000-S microscope (Nikon, Germany).
Measurements of urothelial markers
Expression of urothelial markers was determined by qRT-PCR or immunofluorescence. RNA was extracted by the RNeasy Mini Kit as recommended by the manufacturer (Qiagen, Germany). One μg RNA was reverse transcribed into cDNA using the QuantiTect Reverse Transcription Kit (Qiagen), but with an extended incubation time of 30 min at 42 °C. qRT-PCR was performed with QuantiTect SYBR Green RT-PCR Kit (Qiagen) according to the manufacturer’s instructions. Expression of Cytokeratins (CK) 14 and 20, Uroplakin 2 (UPK2), p53, TERT and p16INK4A mRNAs was determined using primers as listed in Table 1. The housekeeping gene TATA-box binding protein (TBP) was used as a reference gene. Real time qPCR was performed using initial denaturation at 95°C for 15 min and 45 cycles of amplification including denaturation at 94°C for 15 s, annealing for 30 s (temperatures see Table 1) and elongation at 72°C for 30 s. Reactions were performed with the LightCycler 96 platform (Roche, Germany).
For immunofluorescence staining, cells were seeded onto glass slides and allowed to grow for 24 h. After washing with PBS, cells were fixed with 4% formaldehyde in PBS for 30 min and permeabilized by incubation with 0.5% Triton X-100 in PBS for 3 min. Slides were washed, blocked with 1% BSA in PBS for 30 min. Primary antibodies (CK5 1:100, ab-53121; CK14 1:50, ab-7800; CK7 1:200, ab68459; all Abcam, UK; CK7 1:50, MSK032 Zytomed Systems, Germany) and Alexa 488-conjugated or Alexa 594-conjugated secondary antibodies (1:500; Life Technologies) were applied for 1 h each. Nuclei were stained by DAPI (1 μg/ml; Roche) for 5 min. Images were taken using a Nikon Eclipse 400 microscope with the NIS Elements D2-30 software.
Karyotyping, DNA extraction, STR profiling, HPV testing and arrayCGH analysis
STR profiling was performed by standard DNA fingerprint analysis. Cytogenetic characterization was performed as previously described [10]. The karyotype was designated according to the International System for Human Cytogenetic Nomenclature [11]. Schematic presentations of the karyotypic findings were generated by the cytogenetic data analysis system (CyDAS, Germany) [12].
DNA was extracted by a standard proteinase K/phenol/chloroform extraction technique. Comparative genomic hybridization (CGH) was performed using an array of 60mer DNA-oligonucleotide probes with a median overall probe spacing of 13 kb (SurePrint G3 Human CGH Microarray 180 K, Agilent Technologies, USA). Array CGH (aCGH) was done and evaluated essentially according to the protocol of the manufacturer as described elsewhere [10].
HPV testing was performed by means of the digene® HC2 HPV DNA Test (Qiagen) according to the manufacturer’s protocol.
Mutation analysis by sequencing
For p53 mutation analysis, cDNA was amplified by PCR to yield a 1100 bp amplicon containing exons 1-10 of TA-p53 as described [10] (primers # 15 and 16). Mutations in the telomerase reverse transcriptase (TERT) promoter region were detected following PCR amplification. PCR products were purified using the DNA Clean and Concentrator Kit (Zymo Research, USA) and analyzed by Sanger sequencing. Deletion of p16INK4A was detected by amplification of exon 3 from genomic DNA. Primer sequences are detailed in Table 1. Mutation analysis for FGFR3 was performed as described [13].
Library preparation, next generation exome sequencing and data analysis
The exome NGS library was generated from HBLAK high quality genomic DNA with the Ion AmpliSeq™ Exome RDY Kit (Life Technologies, Germany). Amplification and adapter ligation were performed according to the manufacturer’s protocol (MAN0010084 Rev.C, Life Technologies). The library was barcoded during the ligation reaction using the Ion Xpress™ Barcode Adapter Kit (Life Technologies). Library concentration was determined by Qubit measurement (Life Technologies) and library fragment length was analyzed using the 2100 Bioanalyzer (Agilent Technologies, Germany). The library was diluted to a final concentration of 8 pM. Emulsion-PCR and enrichment were performed on the Ion OneTouch™ 2 System (Life Technologies) using the Ion PI™ Hi-Q™ OT2 200 Kit following the manufacturer’s protocol (MAN0010857 Rev.B, Life Technologies). Sequencing was performed on the Ion ProtonTM System (Life Technologies) using Ion PI™ v2 chip with the Ion PI™ Hi-Q™ Sequencing 200 Kit according to the manufacturer’s protocol MAN0010947 Rev.B, Life Technologies).
Amplicon sequences were aligned to the human reference genome GRCh37 (hg19) in the target region of the AmpliSeq™ exome using the Torrent SuiteTM software 4.4.3 with the implemented TMAP algorithm (Life Technologies). Subsequently, variants were detected using the Torrent Variant Caller Plugin v4.4.3.3 and the predefined parameter set ‘somatic low stringency’. Detailed annotations to each detected variant were collected in the databases of the NCBI Reference Sequence Database (RefSeq, version from 2015/03/22), Exome Sequencing Project (ESP6500), 1000 Genomes Project (1000 g), ClinVar (2015/06/29) and dbSNP142 using the open source software tool ANNOVAR [14]. Further information was added by the implemented prediction tools MutationTaster [15–19]. The Integrative Genomics Viewer (IGV 2.3) was used to visualize the read alignments and to check possible errors due to e.g. homopolymeric stretches [20]. Further information on filtering parameters and according numbers of variant is given in supplementary file 1, sheet 1.
Western blotting
Western blot analysis was used to determine p53 and p21 expression subsequent to 24 h of genotoxic stress by treatment with 5 μM/ 10 μM Cisplatin, 100 μM Etoposide, 1 μM Doxorubicin or irradiation with UV light (40 J/m2). Total protein was extracted by lysing the cells in a buffer containing 150 mM NaCl, 1% Triton X-100, 0.5% deoxycholate, 1% Nonidet P-40, 0.1% SDS, 1 mM EDTA, 50 mM Tris (pH 7.6), protease inhibitor cocktail (P-8340, Sigma-Aldrich) and phosphatase inhibitor cocktail (P-0044, Sigma Aldrich) for 30 min on ice. Protein concentration was determined by BCA protein assay (Pierce, Germany). Equal amounts of protein were separated in SDS-PAGE gels and transferred to PVDF membranes (Millipore, Germany). The membrane was blocked with 5% non-fat milk in TBS-T (150 mM NaCl, 10 mM Tris, pH 7.4, and 0.1% Tween-20). Primary antibodies (p53 Ab-6, 1:1000, OP43-100UG, Oncogene Science, USA; p21 Sx118, 1:1000, 556430, BD Biosciences, Germany) were applied overnight at 4°C. α-Tubulin was used as a loading control (1:10000, ab4074, Abcam). After washing, the membrane was incubated with the suitable horseradish peroxidase-conjugated secondary antibody (Santa Cruz, USA) for 1 h and exposed using Super SignalTM West Femto (Thermo Fisher Scientific, Germany) or WesternBright Quantum Kit (Biozym, Germany) according to the manufacturer’s recommendations. Signal intensities were recorded using a LI-COR C-Digit scanner (LI-COR, Germany).
RESULTS
Phenotype of HBLAK cells
HBLAK cells were received from the company and then regularly expanded at a 1:4 ratio for at least 50 further passages. No evidence of diminishing proliferation or morphological transformation was noted. In particular, we did not observe indications of a senescent phenotype during continuous passaging (Fig. 1a). Population doubling time was estimated at about 30 h (Fig. 1b).
Similar to primary cultures of normal human urothelial cells in low-calcium, serum-free KFSM medium (NHUC; Fig. 1c left panel), HBLAK cells cultured as monolayers in low-calcium, serum-free CnT-Prime medium displayed a proliferating and non-differentiated epithelial phenotype with prominent cell-cell contacts (Fig. 1c right panel). To study their differentiation state, RNA expression of urothelial differentiation markers in HBLAK cells was compared to that of three independent cultures of primary uroepithelial cells. CK14, a marker of precursor basal urothelial cells [21] was significantly more strongly expressed in HBLAK cells than in the NHUC cultures (Fig. 1d). Markers of differentiated urothelial cells such as UPK2 and CK20 were variably expressed among the three NHUC cultures, but undetectable in HBLAK. Analysis of protein expression by immunofluorescence revealed that all HBLAK cells expressed CK5 (Fig. 1e, right panel), a marker of basal and intermediate urothelial cells. CK14 was also broadly expressed, but not in all cells (Fig. 1e, left panel). Co-staining displayed that CK14 negative cells were still positive for CK5 (Fig. 1f). HBLAK cells stained also positive for CK7, a common marker of urothelial cells (Fig. 1g). Like CK5, CK7 was expressed in all HBLAK cells, but the staining appeared less homogeneous in intensity (Fig. 1h). Taken together, HBLAK cells phenotypically resemble primary urothelial cell cultures, but consist more homogeneously of basal-like cells compared to NHUC cultures.
Differentiation capacity of HBLAK cells
NHUC cells can be induced to differentiate by two established protocols [2, 3]. In one protocol, simultaneous inhibition of epidermal growth factor receptor (EGFR) signaling with PD153035 and activation of PPARγ by its cognate agonist Troglitazone led to significant induction of the differentiation markers UPK2 and CK20 after one week (Fig. 2a), whereas expression of the precursor marker CK14 was concomitantly lost. In the second protocol, differentiation was induced by increased calcium concentration over a period of two weeks. In addition to induction of terminal differentiation markers, visible stratification of the cultures was observed with this protocol (Fig. 2b bottom). However, induction of differentiation markers appeared to be not as pronounced as in the first protocol, but could be augmented by addition of 5% FBS, which also diminished expression of squamous differentiation markers [3] (data not shown).
Although differentiation markers CK20 and UPK2 were not induced in HBLAK cells by treatment with the EGFR inhibitor and the PPARγ agonist, if the differentiation experiment was performed in the standard CnT-Prime culture medium, which is supplemented with EGF by the manufacturer (Fig. 2c), changes in the morphology of HBLAK cells were similar compared to those of NHUC (Fig. 2d center left). Induction of UPK2 and phenotypical changes could be achieved by further addition of 5% FBS (Fig. 2d center right). In a different low-calcium medium (KFSM), where EGF supplementation is optional, the addition of PD153035 and Troglitazone resulted in a significant induction of CK20 and UPK2 in the absence of EGF (Fig. 2c).
Treatment of HBLAK cells by the high-calcium protocol altered the morphological appearance of the monolayer (Fig. 2d bottom left), but no significant stratification was observed, unlike in NHUC cells (Fig. 2b bottom). In HBLAK cells treated with high calcium, CK14 expression increased, but neither UPK2 nor CK20 became induced (Fig. 2c). Addition of 5% FBS led to a more urothelial-like expression pattern (Fig. 2d bottom right) with diminished CK14 and induction of UPK2, albeit not CK20.
Genetic changes in HBLAK cells
For further characterization of the HBLAK cell line a STR profile was generated as a reference for future analyses and in order to prove that the cell line is unique (Table 2).
Standard karyotyping revealed a hyperdiploid karyotype with trisomies of chromosomes 5 and 20 in all metaphases; two subclones were characterized by additional gains of chromosomes 7 or 9, respectively. No structural chromosomal changes were discernible. The formal description of the karyotype is 48-50, XY, +5, +20[9]/49, XY, +5, +9, +20[5]/50, XY,+5, +7, +20, +20[4]. These findings were confirmed by array-CGH, which also revealed further smaller gains and deletions (Table 3). Notable among these were several distinct deletions on chromosome 9, which was overall gained (i.e. trisomic in many metaphases). The arrayCGH data are best explained by the assumption that all copies of chromosome 9 lack a small 9p21.3 fragment, whereas only one copy contains additional deletions at 9p22.2, 9q21.3 and 9q33-4. The homozygous deletion at 9p21.3 affected a 60 kb region including part of the CDKN2A locus, the MTAP gene and part of the non-coding RNA ANRIL gene. Deletion of CDKN2A/p16INK4 was further confirmed by PCR for CDKN2A exon 3 (Fig. 3a) and for p16 mRNA on samples obtained from different passage numbers. By real time qRT-PCR p16 mRNA expression was undetectable compared to controls (data not shown).
The commonly gained chromosome 5 contains the telomerase reverse transcriptase (TERT) locus at 5p15.33. Mutations in the TERT promoter region leading to elevated telomerase expression are highly prevalent in urothelial carcinoma [22]. Indeed, sequencing of the TERT promoter region in HBLAK revealed the most common C228T mutation (Fig. 3b). TERT mRNA expression was determined by real time qRT-PCR in different passages of HBLAK cells and other cell lines with wildtype (HEK293, 253J) or mutant TERT promoter (VM-CUB1) as well as in TERT-NHUC cells [22] (Fig. 3c). In HBLAK cells, across all passages tested, TERT mRNA was well detectable, other than in normal urothelial cells (NHUC), but was lower than in UC cell lines and in TERT-NHUC cells with engineered telomerase overexpression.
Gains of chromosome 20 are observed in urothelial tumors, but also quite consistently in NHUC cell lines immortalized by HPV oncogenes [23]. We therefore ascertained the lack of HPV DNA in HBLAK by a clinical routine assay.
As p53 inactivation is observed in normal cells immortalized by viral oncogenes and also in many cancerous urothelial cell lines [24, 25], we investigated its status in HBLAK. Expression of the mRNA encoding full-length p53 was well detectable, but not as strongly as in the UC cell lines J82 and VM-CUB1 with mutant p53 (Fig. 3d). Sequencing revealed a wild-type sequence heterozygous for the well-characterized codon 72 G⟶C (Arg⟶Pro) polymorphism. Typical of wild-type p53, the protein was weakly expressed under standard growth conditions (Fig. 4a).
For further characterization of the genetic background of HBLAK cells next generation exome sequencing was performed. Exome sequencing confirmed the deletion of CDKN2A/p16INK4A exon 3 and the lack of p53 mutations. In brief, after filtering, as described in the supplementary file sheet 1, 1104 variants affecting 967 genes were detected compared to the hg19 reference genome (supplementary file, sheet 2). Type and total quantity of variants across the genome are summarized in Table 4. Quantity and type of variants per chromosome are displayed in Fig. 3e. On every chromosome, nonsynonymous single-nucleotide variants (SNV) accounted for the majority of differences, while frameshift and nonframeshift variants were infrequent. Chromosome 1 harbored the most nonsynonymous SNVs, chromosomes 13 and 21 were the least affected. No mutations were detected in genes commonly mutated in urothelial carcinoma (supplementary file, sheets 2 and 3). In particular, FGFR3, the most commonly mutated gene in non-invasive UC and especially in low-grade papillary tumors neither harbored point mutations nor displayed loss of exons indicative of gene fusions, as confirmed by an independent SNaPshot assay.
Sensitivity of HBLAK cells to cytotoxic treatment
Finally, we investigated whether HBLAK cells retained normal responses to cytotoxic treatment. Western blot analysis revealed induction of p53 expression following a 24 h treatment with cytotoxic compounds or UV-irradiation. Induction of p21 expression was also detectable (Fig. 4a).
HBLAK cells were also included in dose response analyses for various pharmaceutical compounds. Exemplarily, we illustrate the sensitivity of HBLAK cells towards treatment with cisplatin for 72 h in Fig. 4b. IC50 concentrations were about 5 μM, higher than those of a primary non-immortalized NHUC, but lower than those of TERT-NHUC cells.
DISCUSSION
Previous reports on artificially immortalized uroepithelial cells have demonstrated altered expression of genes particularly involved in the regulation of differentiation, self-renewal and tumorigenesis, limiting their suitability as models to study these processes [6]. Thus, we characterized the utility and limitations of the new HBLAK cell line, reported as “spontaneously immortalized”, as a cell culture model for normal uroepithelial cells.
Data presented here and in previous studies [26, 27] demonstrate that HBLAK cells behave more like normal primary (NHUC) or telomerase-immortalized urothelial cells (TERT-NHUC) than like urothelial carcinoma (UC) cell lines in drug sensitivity assays and can serve as a convenient control in a similar fashion as TERT-NHUC. One reason for this behavior should be that, in contrast to most urothelial carcinoma cells, HBLAK cells express wildtype p53, which can be induced by treatment with cytotoxic compounds and retain more normal checkpoint signaling. Accordingly, p21 induction can be observed.
The expression pattern of cytokeratins in HBLAK cells corresponds to that of basal or early intermediate urothelial cells, which express CK5 and less regularly CK14, but never CK20 [28]. Accordingly, UPK2 expression, which is weak in intermediate cells and maximal in terminally differentiated umbrella cells, was absent. Importantly, differentiation markers are more homogeneous in HBLAK than in primary NHUC cultures. The homogeneous and reproducible differentiation state of the cell line, in addition to its unlimited propagation, represents another advantage in several applications.
Like primary NHUC cells, HBLAK cells can be induced to differentiate. However, their ability to differentiate is more strongly dependent on the particular protocol than that of NHUC cultures and they do not stratify well. Accordingly, treatment with PD153035 and Troglitazone [2, 3] induced the urothelial differentiation marker UPK2 only in medium without EGF. Likewise, induction of differentiation by high calcium strictly required addition of FBS, which in NHUC enhances calcium-induced urothelial differentiation, but is not as obligatory. Moreover, even after inclusion of FBS in the protocol HBLAK cells retained considerable CK14 expression, indicating incomplete differentiation or a shift towards squamous differentiation, which occurs also when NHUC differentiate in culture under suboptimal conditions [2, 30]. Similar results have been reported for a spontaneously immortalized keratinocyte cell line named HaCaT, which can be induced to differentiate accompanied by the according expression changes of early differentiation markers, but which lost its ability for stratification [31].
A crucial question is how HBLAK became immortalized. Experimentally, NHUC, like many other cell types, can be immortalized by transformation with SV40 largeT-antigen [4] or human papillomavirus (HPV16) E7 [23, 33], which bind and inactivate p53 and RB1 (T-antigen) or RB1 (E7) [24]. Inactivation of p53 and RB1 is also common in invasive UC and accordingly, in many UC cell lines, but essentially absent in papillary UC. Immortalization of NHUCs by DNA virus oncogenes is accompanied by accumulation of additional genetic mutations and restored telomerase activity. Notably, immortalization of NHUC by HPV E7 resulted reproducibly in gain of chromosome 20q, which was also found in HBLAK cells [23]. In a few cell lines derived by HPV E7 transfer, an additional amplification event occurred involving the 20q13.2 region. Other than chromosome 20q gain per se, this amplification was associated with genomic instability, despite retention of wild-type p53 [23]. Interestingly, both 20q gain and 20q13.2 amplification have been reported for various cancer types [33] including UC [34–36], where amplifications of chromosome 20q generally appear to be associated with immortalization, oncogenic transformation and malignant progression. Whereas high level amplifications correlated with tumor progression in breast cancer [37], low-level amplifications appeared to be rather associated with early stage ovarian carcinoma [38]. Similarly, 20q gain also occurred as an early event in normal prostatic cells immortalized by TERT overexpression, where a number of oncogenes in the 20q13.2 region cooperate to support immortalization [39]. Gains at 20q thus appear to represent an early event in the evolution of a malignant phenotype [39]. Taken together, these findings suggest that the gain chromosome 20 in HBLAK cells is not incidental, but contributes to the immortalization of the cell line. However, in HBLAK cells this gain is unrelated to HPV oncogenes as the cell line tested negative for HPV.
Likely, one mechanism underlying spontaneous immortalization of HBLAK cells is telomerase reactivation. In UC, activation of telomerase is a common early event and is typically caused by mutations in the TERT gene promoter [22, 40]. Further, gain of chromosome 5p15.33, harboring the TERT gene, was found to be significantly associated with disease progression in UC [41]. Intriguingly, E7 transformed NHUC cells also gained chromosome 5, which may explain their telomerase overexpression. Accordingly, NHUC can be immortalized by exogenous TERT overexpression [7, 9], like several other cell types [39, 43]. Gain of chromosome 5, including 5p15.33, represents the second consistent chromosomal change in the HBLAK cell line. Moreover, sequencing of the TERT promoter in HBLAK revealed the most common TERT promoter mutation in urothelial carcinoma and TERT mRNA was moderately overexpressed compared to normal urothelial cells, albeit not as strongly as in UC cell lines.
In addition to these obvious chromosomal alterations, further chromosomal alterations in HBLAK cells were not seen across all metaphases or were revealed only by array-CGH. In particular, metaphases contained two or three copies of chromosome 9 with various alterations. All deletions detected in HBLAK can also occur in UC, especially in papillary tumors [44]. Loss of chromosome 9 is the most common chromosomal aberration in UC and the most consistently deleted region is 9p21.3 encompassing CDKN2A encoding the cell cycle inhibitor p16INK4. As loss of p16INK4 contributes to senescence bypass, the homozygous deletion at 9p21.3 removing the 3’-end of CDKN2A in HBLAK cells is likely a third factor permitting their immortalization. Thus, the TERT mutation in conjunction with gain of chromosomes 5 and chromosome 20q and with loss of p16INK4 may concomitantly contribute to immortalization of the HBLAK cell line. Of note, since the original tissue is not available, we do not know whether these changes were already present in the original tissue or arose under culture conditions. However, a number of previous reports on a spontaneously immortalized keratinocyte cell line (HaCaT), immortalization experiments with human mammary epithelial cells (HMEC) and normal prostatic cells suggest that normal epithelial cells can spontaneously acquire the above mentioned alterations in the course of immortalization and that these alterations are sufficient to immortalize normal, non-cancerous cells [31, 49]. Notably, these studies also demonstrate clearly that cells with this limited extent of genetic changes can retain an otherwise normal epithelial phenotype and do not necessarily undergo transformation towards a malignant phenotype. Cells with such genetic alterations therefore do not necessarily originate from cancer cells, but may be susceptible to transformation by further genetic changes.
This interpretation is further supported by the lack of further mutations in genes associated with UC documented by whole exome sequencing and exemplified by the results on two other genes frequently mutated in UC, but unaltered in HBLAK. FGFR3 is the most commonly mutated gene in non-invasive UC and especially in low-grade papillary tumors [50], whereas p53 mutations are more common in invasive UC. In HBLAK, p53 protein size and expression, basal and following DNA damage, fit to a wild-type protein. While MDM2 amplifications do occur in papillary urothelial carcinoma [51], none was found in HBLAK. Similarly, neither arrayCGH nor exome sequencing revealed any relevant alterations in the FGFR3 gene.
Taken together, HBLAK cells harbor some genetic alterations contributing to immortalization of urothelial carcinoma cells (TERT mutation, p16INK4A loss), but lack additional genetic changes characteristic for this tumor type. Notably, we did not obtain evidence for genetic changes characteristic of invasive UC such as p53 mutation and general chromosomal instability. Also a more detailed analysis of the genetic background by exome sequencing confirmed that HBLAK cells harbor only a limited number of genetic changes and, especially, that genes found to be frequently mutated in comprehensive genomic studies of urothelial carcinoma are unaltered in HBLAK. Although neither further primary tissue nor clinical follow-up of the donor is available, our results better support the interpretation that HBLAK cells represent a preneoplastic state of urothelial cells, rather than being derived from an undetected urothelial carcinoma.
In conclusion, HBLAK cells retained, on one hand, many characteristics of normal urothelial cells, including response to cytotoxic agents, allowing its application in studies on tumor-specific effects of novel antineoplastic agents without the limitations by interindividual heterogeneity, limited lifespan and irregular supply of primary cultures. On the other hand, no immortalized epithelial cell line can be truly normal. Interestingly, the immortalization of HBLAK cells appears to be brought about by specific genetic alterations that are thought to lead to immortalization during urothelial carcinogenesis. However, the cells do not contain genetic alterations typical of invasive UC and not even an FGFR3 mutation, suggesting that they represent a preneoplastic stage of urothelial carcinoma. In this regard, they may provide a unique in vitro carcinogenesis model to study factors involved in the stepwise progression of preneoplastic, immortalized urothelial cells towards cancer, alike to the HaCaT cell line used in research on squamous cell carcinoma progression [48]. Although we do not have access to donor tissue in order to ascertain whether the observed genetic changes were already present in the original tissue, none of our results from various analyses indicates that any of the observed genetic alterations was acquired in later passages in culture. We also consider it unlikely that three independent genetic changes favoring immortalization were acquired and fixed in a generally genetically stable cell culture. Accordingly, we propose that immortalized clones of urothelial cells may exist in vivo that contain genetic changes analogous to those experimentally introduced in order to immortalize NHUC.
CONFLICT OF INTEREST
The authors declare that they have no conflict of interest. In particular, they did not receive financial support by the CELLnTEC Company and the company did not exert any influence on the results or conclusions of this study.
Footnotes
ACKNOWLEDGMENTS
We are very grateful to Dr. Peter Girling, CELLnTEC, Bern, Switzerland for providing HBLAK cells and donor information. We gratefully acknowledge help with cell culturing by Marija Lenzen. We thank Dr. Tanja Arent for performing the STR profiling, Dr. Robert Stöhr for the additional FGRF3 mutation analysis and Ursula Grolik for HPV testing. M.P. and G.N. were financially supported by a grant from the Deutsche Forschungsgemeinschaft, E.K. is a postdoctoral fellow of the Bodossaki foundation, and M.A.S. is supported by a fellowship of the Duesseldorf School of Oncology (funded by the Comprehensive Cancer Centre Duesseldorf/ Deutsche Krebshilfe and the Medical Faculty of the Heinrich-Heine-University Duesseldorf).
