A cultured approach to canine urothelial carcinoma: molecular characterization of five cell lines
© Shapiro et al. 2015
Received: 6 May 2015
Accepted: 26 August 2015
Published: 17 September 2015
Urothelial carcinoma (UC), also known as transitional cell carcinoma (TCC), of the bladder is the most common neoplasm affecting the canine urogenital system. To facilitate study of the disease in vitro, cell line models have been established from primary tumor biopsies. Their resemblance to the primary disease, however, has not been well defined. In the present study, we evaluated five canine UC cell lines via oligonucleotide array comparative genomic hybridization (oaCGH), fluorescence in situ hybridization (FISH), and gene expression analysis.
Comparison of genome wide DNA copy number profiles of the cell lines with primary biopsy specimens revealed redundancies in genomic aberrations, indicating that the cell lines retain the gross genomic architecture of primary tumors. As in the primary tumors, gain of canine chromosomes 13 and 36 and loss of chromosome 19 were among the most frequent aberrations evident in the cell lines. FISH analysis revealed chromosome structural aberrations, including tandem duplications, bi-armed chromosomes, and chromosome fusions, suggesting genome instability during neoplastic transformation. Gene expression profiling highlighted numerous differentially expressed genes, including many previously shown as dysregulated in primary canine UC and human bladder cancer. Pathway enrichment analysis emphasized pathways suspected to be at the crux of UC pathogenesis, including xenobiotic and lipid compound metabolism.
These data support valid use of the canine UC cell lines evaluated by confirming they provide an accurate and practical means to interrogate the UC at a molecular level. Moreover, the cell lines may provide a valuable model for furthering our understanding of aberrant metabolic pathways in UC development.
Urothelial carcinoma (UC), also referred to as transitional cell carcinoma, is the most common bladder neoplasm in the dog . While precise lifetime risk and incidence numbers in pet dogs are unknown, UC is estimated to affect more than 20,000 dogs each year in the United States . Due to the uncomfortable and potentially fatal consequences associated with a rapidly growing bladder tumor, the diagnosis of UC conveys a guarded prognosis and evokes considerable dog owner concern . Additionally, over 90 % of canine bladder tumors are invasive with metastatic potential; approximately 20 % of canine patients show overt metastases at the time of diagnosis and over 60 % at death . The devastating clinical course of disease, combined with a high rate of metastasis, emphasizes a need for early detection of the tumor. With improved detection and a better understanding of the disease, therapeutics may prove more efficacious, extend patient survival, and improve quality of life.
Previous studies have shown that the genomic landscape of canine UC is highly aberrant, with recurrent chromosome copy number aberrations affecting gene dosage on, most notably, canine chromosomes (CFA) 13, 19, and 36 . In addition, overwhelming breed predispositions suggest the risk for UC development is at least partially genetic. The inbred history of purebred dogs, resulting in reduced levels of genetic variation and enrichment in breed-associated loci, provides a biologic basis for breed susceptibility . Scottish terriers, for example, are diagnosed with UC 20 times more frequently than the average of all breeds, and Beagles, Shetland sheepdogs, Wire hair fox terriers, and West Highland white terriers are reported to be affected three to five times as frequently . The uniquely homogenous genetic population of purebred dogs makes the dog an excellent model for teasing out genetics that are potentially important in cancers of the highly heterogeneous human population. Previous studies also have found that canine UC closely resembles that of the human counterpart, in histopathological features, clinical behavior, and genomic aberrations, furthering the value of the dog to human medicine [2, 6, 4].
In addition to genetic risk factors, UC is associated with numerous environmental factors, including flea and tick dips, obesity, pesticides, and cyclophosphamide treatment [7–9]. While correlations of UC incidence with each of these risk factors have been found, mechanisms of carcinogenesis have yet to be elucidated. Problems in cyclic compound metabolism have been proposed, but until the mechanisms of UC development are elucidated, the provision of effective prevention and treatment strategies will remain challenging . The study of potential molecular drivers and therapeutic targets would be greatly expedited by a biologically appropriate in vitro model.
Cell lines provide a useful in vitro model for the study of disease. By facilitating characterization and manipulation of all molecular facets without harming the model animal, cell lines enable researchers to evaluate highly experimental ideas and therapies. Although canine UC cell lines have previously been established for such use , their resemblance to the primary tumor in the context of molecular behavior and genomic landscape has yet to be thoroughly evaluated. Affirming the cell lines recapitulate molecular characteristics of the primary tumor would increase their clinical predictive value and enhance their significance in translational studies.
In the present study, we investigated the molecular profile of five canine UC cell lines with oligonucleotide array comparative genomic hybridization (oaCGH) and fluorescence in situ hybridization (FISH). Clinical relevance of findings was validated by comparison of oaCGH profiles of cell lines with those of primary tumor samples, which showed conservation of major genomic aberrations in both sample cohorts. Subsequent gene expression profiling quantified mRNA from copy number neutral and aberrant regions of the genome. KEGG pathway and gene ontology (GO) analysis demonstrated that genes involved in lipid metabolism and cell cycle regulation were enriched among genes differentially expressed between neoplastic cell lines and healthy urothelium, highlighting a potential pathogenesis that is relevant to known risk factors. Our results affirm that the canine UC cell lines are genomically similar to the primary tumor, providing a relevant in vitro model for study of the molecular mechanisms of disease.
Characteristics of five canine UC cell lines
Cell Line ID
Average chromosome #
Specific structural changes
Control samples for expression analysis (healthy urothelium) were collected during necropsy of dogs with no clinical or histopathological signs of disease at North Carolina State University College of Veterinary Medicine. After collection, bladder tissue was stored in transport media (RPMI media supplemented with 10 % fetal bovine serum (FBS) until sample processing (<1 h). Samples were split in half: one half was fixed in 10 % neutral buffered formalin (NBF) before paraffin-embedding (FFPE) to enable pathologic evaluation, and the second half was used for urothelial cell isolation. FFPE samples were sectioned and stained with hematoxylin and eosin (H&E) for review by a board certified veterinary pathologist. All control samples were found to be histopathologically normal prior to use in subsequent protocols. Urothelial cells were isolated by scraping the bladder mucosa three times with a glass slide and rinsing with 1x phosphate buffered saline (PBS), as described previously . Urothelial scraping was performed in order to minimize contamination of urothelium by submucosal tissue layers. Cell scrapings were then analyzed under phase contrast to determine the proportion of epithelial cells to contaminating mesenchymal cells. Over 90 % of isolated cells were confirmed to be epithelial.
Cells from each of the five lines were cultured in DMEM/F12 without glutamine (Mediatech, Manassas, VA) and supplemented with 10 % FBS (Mediatech), 1 % Glutamax™ (Life Technologies, Carlsbad, CA), and 0.6 % Primocin™ (InvivoGen, San Diego, CA) to confluence, passaging when confluent. For each line, cells from the same flask were split into three aliquots and used simultaneously for chromosome preparation, DNA isolation, and RNA isolation. Four of the cell lines (K9TCC-PU-An, K9TCC-PU-In, K9TCC-PU-Mx, K9TCC-PU-Sh) were available as low passage (<p7) and so were harvested at an earlier (<p10) and later (>p16) passage to evaluate genomic evolution during culture. One cell line (K9TCC) was available only from p14 and so was harvested at p16 only.
DNA from each cell line ("test") was isolated using a DNeasy Kit according to manufacturer’s recommendations (Qiagen, Valencia, CA). Purified DNA was verified to be of high molecular weight and purity by agarose gel electrophoresis and spectrophotometry (Nanodrop™ 1000, Thermo Fisher Scientific, Wilmington, DE; 260:230 > 2.0 and 260:280 > 1.8). Since all five cell lines were female, a single female "reference" mix was generated by combining equimolar quantities of DNA extracted from ten different, healthy female dogs of various breeds (QIAmp DNA Midi Kit, Qiagen). Isolated DNA was labeled using the Genomic DNA Enzymatic Labeling Kit (Agilent Technologies, Santa Clara, CA) as described previously [4, 13]. Fluorescently labeled test and reference DNA samples were hybridized to Canine G3 Sureprint 180,000 feature oaCGH arrays (Agilent, AMADID 025522) for 40 h, as described previously [13, 14]. Arrays were scanned at 3 μm (Agilent, Model G2505C) and data extracted with Feature Extraction software (v10.9 Agilent) using the canFam2 genome build. CGH results were analyzed using the FASST2 algorithm in Nexus Copy Number (Biodiscovery, Hawthorne, CA). After aberration detection, cell lines and primary tumors were clustered based on aberrations using a complete linkage hierarchical clustering algorithm by Nexus Copy Number (Biodiscovery).
Genomic regions likely to be central to tumorigenesis were pinpointed by the Genomic Identification of Significant Targets in Cancer (GISTIC) algorithm in Nexus Copy Number, which identifies regions unlikely to be aberrant by chance when taking aberration frequency and amplitude into account . A G-score is computed to reflect the frequency and amplitude of the aberration, while a q-value suggests the likelihood (probability) of that aberration occurring by chance when looking at the overall aberration pattern throughout the genome. Lastly, aberrations in early and late passage cell lines were compared using the “Comparisons” tool in Nexus, using the early passage as baseline (p = 0.05, differential threshold = 1 %).
Primary tumor oaCGH data
Previously published oaCGH data from 31 primary canine UCs was compared to cell line aberrations . None of the cell lines used in the current study were derived from any of the 31 primary tumors. Using Nexus Copy Number, aberrations common to primary tumors and cell lines were compared using the “Comparisons” tool. Additionally, GISTIC analysis was repeated using primary tumor aberrations.
Fluorescence in situ hybridization (FISH)
Cells were harvested using conventional methods of colcemid-induced metaphase arrest (Life Technologies), hypotonic (KCl) treatment (Life Technologies), trypsinization (0.25 % trypsin, Mediatech), and fixation, as described previously . Fixed cell suspensions were dropped onto clean glass slides and aged for three days at room temperature prior to ethanol dehydration (70 %/90 %/100 %) and storage at −80 °C.
Regions chosen for FISH validation of oaCGH
Using Log2 ratios ascertained by oaCGH and also deduced from single locus-BAC probe probe enumeration, correlation coefficients were derived. To calculate the Log2 ratio, the average copy number of the aberrant locus was divided by the average copy number of the CFA11 neutral locus of each cell line. A Log2 of the resultant value was calculated and used as the FISH-derived Log2 for each copy number-aberrant locus. Additionally, a Mann–Whitney U test (two-tailed, p < 0.05) was performed using JMP Professional Statistical Software (v. 11, SAS, Cary, NC) to compare mean Log2 ratios from both methods.
RNA was extracted from cultured cells or control urothelial cells using the RNeasy Plus Mini Kit (Qiagen). RNA was immediately isolated from the healthy, fresh urothelial cells to minimize alterations in post-mortem and post-collection gene expression. RNA purity and integrity was assessed using spectrophotometry (Nanodrop) and a 2100 BioAnalyzer (Agilent), respectively. Samples with 260:230 and 260:280 > 2.0 and RIN > 8.0 were used for microarray analysis and qRT-PCR validation. RNA was labeled using the Quick Amp Labeling kit (Agilent). Purified RNA probes were hybridized to a one-color expression microarray (Agilent SurePrint G3 Canine 4x44k Expression Array). Arrays were scanned at 3 μm (Agilent, Model G2505C). Analysis was performed using the GeneSpring advanced platform (Agilent, v.11.5, 2011) and Nexus Expression (BioDiscovery, v.2.0, 2010). Prior to expression normalization, compromised probes were removed from analysis in GeneSpring, including features which were non-uniform, saturated, or population outliers. Normalized signal values were obtained by log transformation, followed by 75th percentile shift normalization and baseline shift to median of all samples. Normalized signal values (NS) were compared between RNA isolated from the cell lines and from urothelium of two histologically confirmed healthy dog bladders to determine the fold change in expression (NScell_line ÷ NScontrol_average = relative fold changegene). Significantly over- and underexpressed genes in demonstrated a ± 2-fold change in gene expression in the cell lines relative to the average of control bladders. The resulting gene expression profiles of cell lines and healthy control urothelium were clustered using an agglomerative, unsupervised hierarchical clustering algorithm in Nexus Expression (Biodiscovery).
Copy number aberrant genes analyzed by qRT-PCR
Canine UC cell lines display recurrent chromosome copy number and structural aberrations
Regions implicated by GISTIC analysis represent regions relevant to neoplasia
Compiled FISH and oaCGH Log2 data
Gene expression profiles of cell lines differ from those of normal urothelium
When profiled alongside healthy urothelium, the UC cell lines showed highly differential gene expression levels. Of 43,803 probes on the array, 8,688 probes, representing at least 3,303 Ensembl-annotated canine genes, displayed altered levels of expression with absolute fold change greater than 2.0 when normalized and compared to reference urothelium. Examples are shown in Table 4 and were selected based on the magnitude of fold change and potential relevance to tumorigenesis. Hierarchical expression clustering confirmed a vastly different expression landscape between tumors and control urothelium, each recognized as a separate outgroup based on DEGs (Additional file 2: Figure S1).
Using the DAVID tool for gene ontology (GO) and pathway analysis, DEG enrichment in functional pathways was evaluated. Among underexpressed genes (1859 total), the GO term “lipid metabolism” was the only GO term enriched (p = 0.03) in our dataset, with 189 genes underexpressed in the cell lines. Among overexpressed genes (1444 total), “DNA replication, recombination, and repair” was the most highly enriched term (p = 0.01), with 13 overexpressed genes. When evaluating all DEGs, “DNA replication and repair” was the single most commonly dysregulated functional ontology, represented by 10 % of all cell line DEGs.
KEGG pathway analysis in DAVID highlighted four pathways enriched for underexpressed genes: metabolism of xenobiotics by cytochrome p450, glutathione metabolism, drug metabolism, and fatty acid metabolism, resulting in a KEGG pathway clustering of these pathways with an enrichment score indicative of biological significance (Bonferroni = 0.05, enrichment score = 2.34). Additionally, 3.7 % of genes with five-fold or greater expression than normal urothelium (10 of 268) were involved in cell cycle regulatory pathways (p = 0.0004), as were 1.8 % (26 of 1444) of all upregulated genes (p = 0.00003). DAVID chromosome association analysis found 4.7 % of overexpressed genes (68 genes) were found on CFA 13 and 4.2 % on CFA 10 (61 genes). Similarly, 3.6 % (66 genes) of underexpressed genes were located on CFA 12 and 1.2 % on CFA 19 (23 genes).
Canine UC represents the most common genitourinary cancer in dogs and, due to the aggressive nature of the tumor, diagnosis conveys a guarded prognosis for the patient . To date, curative therapies are lacking for UC, which is likely due in part to delayed tumor diagnosis and the lack of knowledge about the molecular basis of tumor pathogenesis. It is known, however, that canine UC closely resembles prognostically unfavorable invasive human UC, making research in canine bladder cancer potentially valuable in improving the outlook for human UC patients [6, 2, 1, 11]. As a result, a valid in vitro cell line model would be advantageous to both species and could potentiate investigation of more precise molecular mechanisms involved in UC pathogenesis, including the development of therapies to target them. Previously, five canine UC cell lines were established for in vitro study of the disease, and four of them retain tumorigenicity in mice . With the potential to facilitate translational studies, care should be taken to determine that these cell lines provide a high clinical predictive value. Our study sought to evaluate: 1) the genomic similarity between primary UC and cell lines, 2) the genomic integrity of the cell lines over time, and 3) the gene expression landscape of canine UC cell lines.
Our results show that, although variation exists in the frequency of copy number aberrations between the cell lines and primary tumors, the most frequent copy number gains and losses are preserved in the cell lines, specifically loss of regions of CFA 19, and gain of regions of CFA 13 and 36 . Among primary tumors, 100 % of the canine UC cases possess at least one of these aberrations, 94 % possess two, and 68 % possess all three aberrations, emphasizing the diagnostic potential of copy number aberrations and a potentially critical role of these chromosomes in tumorigenesis. Similarly, all five cell lines had at least one of these aberrations (100 %), four had at least two (80 %), and two had all three aberrations (40 %). When oaCGH-called aberrations among the 31 primary tumors and five cell lines were clustered, no cell line outgroup was seen, suggesting the cell lines were more like primary tumors than one another other and demonstrating a lack of culture artifact in their genomic profiles. Thus, the cell lines maintain the valuable aberrations of the primary tumor and provide an in vitro tool with which to further investigate the significance of these aberrations in canine UC.
The relevance of the gains on CFA13 and 36 are further emphasized by GISTIC analysis. GISTIC analysis, which focuses on regions containing genes repeatedly implicated in neoplastic transformation and tumor growth, highlighted only seven regions aberrant in the cell lines: amplified regions on CFA10, 36, and 38 and gains on CFA13. In primary tumors, however, GISTIC analysis implicated numerous regions, including the majority of those seen in the cell lines. These results not only emphasize the importance of genes highlighted by GISTIC analysis, but may also suggest that a number of the mutations seen in the primary tumor may be passenger mutations unessential for tumor maintenance in the cell lines. Thus, we may be able to utilize the less genomically complex cell lines to hone in on driver mutations crucial for UC pathogenesis and progression.
In addition to major chromosome aberrations, chromothriptic-like events represent another manifestation of copy number change. First described in 2011, chromothriptic-like events represent small genomic regions comprised of numerous adjacent and alternating copy number changes [19, 20]. Chromothriptic-like events were previously found to be prevalent throughout the genome of primary canine UC, particularly on CFA36, 10, 16, 4 and 7 . Similarly, the cell lines show a high frequency and number of chromothriptic-like events throughout the genome, with 21 of 39 chromosomes demonstrating recurrent chromothriptic-like events in at least two cell lines. Among these are all of those observed at high frequency in primary tumors, except CFA16. The prevalence of chromothriptic-like events in the cell line genomes suggests inherent genome instability along with a massive chromosome shattering event. For example, tandem duplications observed on CFA10, 36, and 38, as well as translocations among these chromosomes (t(36;38)), may have arisen due to such chromothriptic events early in or leading to neoplastic transformation. Although arising from an unknown etiology, current hypotheses of chromothripsis in cancer implicate mutagenic exogenous compounds, which induce double-stranded DNA breaks throughout the genome, and chromosome missegregation [21, 22]. This type of chromosome shattering can lead to massive regional amplification, as well as losses and translocations, due to improper chromosome rejoining. Particularly when considering the potential etiology of risk-associated carcinogenic compounds, chromothripsis is a hugely interesting phenomenon in UC.
In addition to validating oaCGH-called aberrations, the ability to perform metaphase FISH on the cell lines allowed us to explore karyotypic reorganization occurring in canine UC, including recurrent tandem duplications responsible for CFA36 amplification and a CFA13 duplication and centromeric fusion often responsible for increases in CFA13. In our cell line analyses, copies of CFA10, 36, and 38 per cell were often too numerous to accurately count (>10-30 clearly distinguishable signals, particularly in K9TCC-Pu-Mx). Furthermore, in all cell lines, abnormally metacentric chromosomes were seen, suggesting prior mitotic dysfunction and/or breakage-fusion-bridge cycles and leading to the formation of new hybrid chromosomes.
The hyperdiploid karyotype of the cell lines, excluding K9TCC, further confirmed previous findings that canine UCs are often tetraploid (>50 % of UCs) . Since oaCGH is blind to ploidy, the correlation between oaCGH status and copy number requires knowledge of ploidy status, making a copy number neutral control critical to appropriate copy number interpretation. CFA11 was identified as copy number neutral using oaCGH values and FISH analysis; CGH values suggestive of neither a loss nor gain (Log2 ≈ 0 for all cell lines) and chromosome structure appeared normal. In polyploid cell lines, four copies of CFA11 were consistently observed, supporting a tetraploid genome. Similarly, when oaCGH implied a copy number loss, one to three copies were observed, showing a relative, but often not absolute, loss of genetic material. In the case of CFA19, only in diploid K9TCC did the chromosome count consistently reflect an absolute copy number loss (n < 2). Relative to a duplicated genome, however, two copies of CFA 19 (K9TCC-PU-An) reflected a loss rather than a normal diploid number of Log2 = 0.
Cells were repeatedly passaged (>10 passages) in order to determine the stability of the UC genomic profile over time. Little change was evident in oaCGH profiles or chromosome structure of cell lines after passaging cells from early to mid-late stage. Only three regions on CFA1, 2, and 9 were significantly different between early and late passages, representing less than 0.006 % of the entire genome in length. Although the region on CFA9 contains no known genes, the regions on CFA1 and 2 contain genes potentially relevant to the culture environment: insulin-like growth factor 2 receptor (IGF2R, CFA1) and a protocadherin cluster (PCDH, CFA2). IGF2R and PCDH are both potential adaptations to the culture environment, allowing improved utilization of media sugars and improving culture flask adhesion, respectively. Nevertheless, such few aberrations acquired over 10+ passages represent a very minor effect of repeated passage and cell culture. This characteristic of the cell lines allows experimentation and analysis over time without concern of culture-induced changes, a valuable feature for translational studies.
We should note that several factors might influence our results. One is that cell lines are ex vivo models. Thus, some aberrations detected in cell lines may be due to the altered culture environment. We were unable to analyze the primary tumor from which the cell lines were derived, precluding our ability to detect initial culture-related genomic changes. Secondly, since oaCGH is a cell population-based analysis, it is plausible that rare aberrations evaded detection by oaCGH, but would become more prominent with continued neoplastic cell proliferation. However, our cell lines sustained the major copy number changes seen in primary tumors, even after ten or more passages, leading us to conclude cell lines represent a good model for UC in vitro.
Factors incriminated in the development of UC include exposure to cyclic amines, including carcinogens in cigarette smoke, estrogen, and cyclophosphamide, as well as obesity [7, 6, 2, 1, 24–27]. Interestingly, DAVID analysis of genes underexpressed in the UC cell lines revealed GO enrichment of genes essential to the metabolism of such lipophilic cyclic compounds. Genes underexpressed in the UC cell lines were also enriched in metabolic KEGG pathways, most notably pathways involved in xenobiotic metabolism. Decreased expression of metabolic genes suggests that toxic compounds are not being handled properly by the bladder epithelium, potentially leading to a buildup of corrosive compounds in the urine and increased exposure of the urothelium to carcinogens. In particular, carboxylesterase-1 (CES-1; CFA2; fold change = −633.5) is involved in metabolism of organophosphate and pyrethroid insecticides, both associated with UC development [28–30]. Similarly, decreased to undetectable levels of urothelial CES1 are associated with UC in humans . In addition, CES enzymatic activity is highly variable among human individuals and may be associated with differences in canine breed susceptibilities [32, 33].
Of further interest, Cytochrome P450 2C19 (CYP2C19; CFA28; fold change = −8.1) is a potentially relevant metabolic gene central to the enriched cytochrome p450 KEGG pathway. CYP2C19 is one of the most important cytochrome p450 enzymes and is responsible for the metabolism of numerous cyclic amine xenobiotics, including the pyrethroid-derivatives of hydantoin present in common insecticides . Similar to CES1, genetic polymorphisms of CYP2C19 are associated with reduced metabolic capacity in humans, suggesting similar phenomena may occur among genetically distinct dog breeds [35, 36]. Our data corroborate previous speculations regarding the involvement of metabolic pathways in canine UC development and highlights a need for closer investigation.
Among genes overexpressed in the UC cell lines, GO terms were enriched for DNA replication, recombination, and repair. Similarly, KEGG analysis highlighted cell cycle regulator pathways, in which 4 % of genes with expression five-fold above that of normal urothelium were involved (1.8 % of all overexpressed genes). One of the most prominent overexpressed genes in our dataset and a key gene in DNA replication and cell cycle regulation--PTTG1 (CFA4; fold change=+57.5). As a result of its functions as a sister chromatid securin, microtubule nucleation regulator, and AKT activator, among others, PTTG1 overexpression potentiates tumor proliferation and invasion, as well as chromosome instability [37–39]. Interestingly, activating mutations in FGFR3, an upstream regulator of AKT activity, are implicated in constitutive activation of the AKT pathway commonly seen in human UC . PTTG1 has the potential to be a major driver of AKT activation in canine UC, providing a valuable therapeutic target and highlighting pathway dysregulation similarities in human and dog UC.
Our expression data also emphasized the potential impact of chromosome copy number on gene dosage and, therefore, gene expression. When analyzing genes located in regions of copy number aberration, it was shown that gene expression often varied directly with copy number change (r = 0.76). DAVID chromosome analysis showed enrichment of underexpressed genes on CFA19, which is frequently lost in UC cell lines. CFA 19 contains the gene for histamine N-methyltransferase (HNMT), a gene highly under expressed in our cell lines (fold change = −29.6). Although unreported in bladder cancer, reduced HNMT transcription has been reported in other human carcinomas and is posited to be involved in tumorigenesis . Decreased levels of HNMT lead to local increases of angiogenesis-promoting PTGS2 (COX2; CFA7). PTGS2 levels are increased in both canine and human UC, and this gene is similarly overexpressed in our cell line data (fold change = +76.9), providing a further mechanistic role for COX-2 inhibitors in UC treatment [42, 43]. Additionally, COX2 is a downstream effector of an activated AKT pathway, again suggesting a conservation of AKT overactivation in both canine and human UC development.
Gene expression microarray results
Lines differentially expressing
Average fold change
Lines with relative fold change ≥2 or ≤ −2
Average fold change
In, An, Mx, K9TCC
In, Sh, Mx, K9TCC
In, Sh, An, K9TCC
In, An, Mx, K9TCC
An, Mx, K9TCC
Despite strong correlations between copy number gain/loss and differential gene expression (Fig. 4), it is not known if increased genomic dosage directly leads to increased transcript levels. In fact, PABPC1, a gene located in the region of shared high frequency gain on HSA 8/CFA 13 in human and canine UC, was actually underexpressed (average fold change = −2.6) in the cell lines . However, PABPC1 is highly expressed, at both mRNA and protein levels, in normal urothelium, potentially obscuring subtle differences in expression while reducing their physiologic significance [45, 46]. Additionally, while human expression studies have shown relative overexpression of PABPC1 at the mRNA level, protein studies have shown a relative decrease in protein, espousing the need for protein evaluation alongside mRNA quantification . We cannot discount the importance of epigenetic factors in determining expression. In addition, disruption of gene promoters and enhancers, as may occur due to tandem duplication, deletion, and/or translocation, may alter predicted gene expression. Post-transcriptional gene regulation may also play a role due to possible alterations in silencing RNA. Regardless, it is reasonable to assume genomic copy number plays a role in determining corresponding mRNA levels.
The recurrence of specific genetic aberrations in UC emphasizes the importance of continued UC genomic research, particularly in predisposed dog breeds with reduced genetic variation and in association with prominent carcinogenic risk factors. A treatment found to be effective in canine cell lines, which then is effective in dogs, increases the clinical predictive value of the cell lines and, therefore, their value to biomedical research . The preservation of primary tumor aberrations in the cell line model provides an accessible and accurate means of performing genetic experiments in vitro, making the cell lines a valuable resource for translational UC research. Future research should involve utilizing the UC cell lines in a functional manner to elucidate mechanisms of UC pathogenesis.
This study was supported by the NCSU Cancer Genomics Fund (MB). Additional support for SS came from an NIH-T35 grant, a Triangle Community Foundation Hitching’s Young Investigator award, and the Comparative Biomedical Sciences DVM/PhD Program at NC State University College of Veterinary Medicine.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Mutsaers AJ, Widmer WR, Knapp DW. Canine transitional cell carcinoma. J Vet Intern Med. 2003;17(2):136–44.View ArticlePubMedGoogle Scholar
- Knapp DW, Ramos-Vara JA, Moore GE, Dhawan D, Bonney PL, Young KE. Urinary Bladder Cancer in Dogs, a Naturally Occurring Model for Cancer Biology and Drug Development. ILAR J. 2014;55(1):100–17.View ArticlePubMedGoogle Scholar
- Hamilton MJ, Sarcornrattana O, Illiopoulou M, Xie Y, Kitchell B. Questionnaire-based assessment of owner concerns and doctor responsiveness: 107 canine chemotherapy patients. J Small Animal Pract. 2012;53(11):627–33. doi:10.1111/j.1748-5827.2012.01269.x.View ArticleGoogle Scholar
- Shapiro SG, Raghunath S, Williams C, Motsinger-Reif AA, Cullen JM, Liu T et al. Canine urothelial carcinoma: genomically aberrant and comparatively relevant. Chromosome Res. 2015. doi:10.1007/s10577-015-9471-y
- Dobson JM. Breed-predispositions to cancer in pedigree dogs. ISRN Vet Sci. 2013;2013.Google Scholar
- Knapp D, Glickman N, DeNicola D, Bonney P, Lin T, Glickman L. Naturally-occurring canine transitional cell carcinoma of the urinary bladder: A relevant model of human invasive bladder cancer. Urol Oncol. 2000;5:47–59.View ArticlePubMedGoogle Scholar
- Glickman LT, Raghavan M, Knapp DW, Bonney PL, Dawson MH. Herbicide exposure and the risk of transitional cell carcinoma of the urinary bladder in Scottish Terriers. J Am Vet Med Assoc. 2004;224(8):1290–7.View ArticlePubMedGoogle Scholar
- Plotz PH, Klippel JH, Decker JL, Grauman D, Wolff B, Brown BC, et al. Bladder complications in patients receiving cyclophosphamide for systemic lupus erythematosus or rheumatoid arthritis. Ann Intern Med. 1979;91(2):221–3.View ArticlePubMedGoogle Scholar
- Glickman LT, Schofer FS, McKee LJ, Reif JS, Goldschmidt MH. Epidemiologic study of insecticide exposures, obesity, and risk of bladder cancer in household dogs. J Toxicol Environ Health. 1989;28(4):407–14. doi:10.1080/15287398909531360.View ArticlePubMedGoogle Scholar
- Zhang J, Wei S, Liu L, Nagana Gowda GA, Bonney P, Stewart J, et al. NMR-based metabolomics study of canine bladder cancer. Biochim Biophys Acta. 2012;1822(11):1807–14. doi:10.1016/j.bbadis.2012.08.001.View ArticlePubMedGoogle Scholar
- Dhawan D, Ramos-Vara JA, Stewart JC, Zheng R, Knapp DW. Canine invasive transitional cell carcinoma cell lines: in vitro tools to complement a relevant animal model of invasive urinary bladder cancer. Urol Oncol. 2009;27(3):284–92. doi:10.1016/j.urolonc.2008.02.015.View ArticlePubMedGoogle Scholar
- Harris Jr HW, Wade JB, Handler JS. Transepithelial water flow regulates apical membrane retrieval in antidiuretic hormone-stimulated toad urinary bladder. J Clin Invest. 1986;78(3):703–12. doi:10.1172/JCI112630.PubMed CentralView ArticlePubMedGoogle Scholar
- Poorman K, Borst L, Moroff S, Roy S, Labelle P, Motsinger-Reif A et al. Comparative cytogenetic characterization of primary canine melanocytic lesions using array CGH and fluorescence in situ hybridization. Chromosome Res. 2015. 23:171-186.Google Scholar
- Roode S, Rotroff D, Avery AC, Suter SE, Bienzle, D, Schiffman JD, et al. Genome-wide assessment of recurrent genomic imbalances in canine leukemia identifies evolutionarily conserved regions for subtype differentiation. Chromosome Res. 2015;on-line June 3rd.Google Scholar
- Beroukhim R, Getz G, Nghiemphu L, Barretina J, et al. Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. PNAS. 2007;104(50):20007–12.Google Scholar
- Breen M, Langford CF, Carter NP, Holmes NG, Dickens HF, Thomas R, et al. FISH mapping and identification of canine chromosomes. J Hered. 1999;90(1):27–30.View ArticlePubMedGoogle Scholar
- Breen M, Hitte C, Lorentzen T, Thomas R, Cadieu E, Sabacan L, et al. An integrated 4249 marker FISH/RH map of the canine genome. BMC Genomics. 2004;5(1):65.PubMed CentralView ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13. doi:10.1093/nar/gkn923.PubMed CentralView ArticleGoogle Scholar
- Stephens PJ, Greenman Cd, Fu B, Yang F, et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell. 2011;144(1):27–40.Google Scholar
- Cai H, Kumar N, Bagheri HC, von Mering C, Robinson MD, Baudis M. Chromothripsis-like patterns are recurring but heterogeneously distributed features in a survey of 22,347 cancer genome screens. BMC Genomics. 2014;15:82. doi:10.1186/1471-2164-15-82.PubMed CentralView ArticlePubMedGoogle Scholar
- Holland AJ, Cleveland DW. Losing balance: the origin and impact of aneuploidy in cancer. EMBO Rep. 2012;13(6):501–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Stephens PJ, Greenman CD, Fu B, Yang F, Bignell GR, Mudie LJ, et al. Massive Genomic Rearrangement Acquired in a Single Catastrophic Event during Cancer Development. Cell. 2011;144(1):27–40. doi:http://dx.doi.org/10.1016/j.cell.2010.11.055.PubMed CentralView ArticlePubMedGoogle Scholar
- Clemo FAS, DeNicola DB, Carlton WW, Morrison WB, Walker E. Flow Cytometric DNA Ploidy Analysis in Canine Transitional Cell Carcinoma of Urinary Bladders. Vet Pathol. 1994;31(2):207–15. doi:10.1177/030098589403100208.View ArticlePubMedGoogle Scholar
- Lawson M, Vasilaras A, De Vries A, Mactaggart P, Nicol D. Urological implications of cyclophosphamide and ifosfamide. Scand J Urol Nephrol. 2008;42(4):309–17. doi:10.1080/00365590701570953.View ArticlePubMedGoogle Scholar
- Zhang Y. Understanding the gender disparity in bladder cancer risk: the impact of sex hormones and liver on bladder susceptibility to carcinogens. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2013;31(4):287–304. doi:10.1080/10590501.2013.844755.View ArticlePubMedGoogle Scholar
- Talaska G. Aromatic amines and human urinary bladder cancer: exposure sources and epidemiology. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2003;21(1):29–43. doi:10.1081/gnc-120021372.View ArticlePubMedGoogle Scholar
- Brody JS. Transcriptome alterations induced by cigarette smoke. Int J Cancer. 2012;131(12):2754–62. doi:10.1002/ijc.27829.View ArticlePubMedGoogle Scholar
- Wheelock CE, Miller JL, Miller MJ, Phillips BM, Huntley SA, Gee SJ, et al. Use of carboxylesterase activity to remove pyrethroid-associated toxicity to ceriodaphnia dubia and hyalella azteca in toxicity identification evaluations. Environ Toxicol Chem/SETAC. 2006;25(4):973–84.View ArticleGoogle Scholar
- Nishi K, Huang H, Kamita SG, Kim I-H, Morisseau C, Hammock BD. Characterization of pyrethroid hydrolysis by the human liver carboxylesterases hCE-1 and hCE-2. Arch Biochem Biophys. 2006;445(1):115–23. doi:10.1016/j.abb.2005.11.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang D, Pearce RE, Wang X, Gaedigk R, Wan Y-JY, Yan B. Human Carboxylesterases HCE1 and HCE2: Ontogenic Expression, Inter-Individual Variability and Differential Hydrolysis of Oseltamivir, Aspirin. Deltamethrin Permethrin Biochem Pharmacol. 2009;77(2):238–47. doi:10.1016/j.bcp.2008.10.005.View ArticlePubMedGoogle Scholar
- Pontén F, Jirström K, Uhlen M. The Human Protein Atlas—a tool for pathology. J Pathol. 2008;216(4):387–93. doi:10.1002/path.2440.View ArticlePubMedGoogle Scholar
- Sanghani SP, Sanghani PC, Schiel MA, Bosron WF. Human Carboxylesterases: an update on CES1, CES2 and CES3. Protein Peptide Lett. 2009;16(10):1207-14.Google Scholar
- Merali Z, Ross S, Pare G. The pharmacogenetics of carboxylesterases: CES1 and CES2 genetic variants and their clinical effect. Drug Metabol Personalized Ther. 2014;29(3):143–51.Google Scholar
- Suzuki H, Kneller MB, Rock DA, Jones JP, Trager WF, Rettie AE. Active-site characteristics of CYP2C19 and CYP2C9 probed with hydantoin and barbiturate inhibitors. Arch Biochem Biophys. 2004;429(1):1–15. doi:http://dx.doi.org/10.1016/j.abb.2004.05.015.
- Tamer L, Ercan B, Ercan S, Ateş N, Ateş C, Öcal K, et al. CYP2C19 polymorphisms in patients with gastric and colorectal carcinoma. Int J Gastrointest Canc. 2006;37(1):1–5. doi:10.1385/IJGC:37:1:1.View ArticleGoogle Scholar
- Ruiter R, Bijl MJ, van Schaik RHN, Berns EMJJ, Hofman A, Coebergh J-WW, et al. CYP2C19*2 polymorphism is associated with increased survival in breast cancer patients using tamoxifen. Pharmacogenomics. 2010;11(10):1367–75. doi:10.2217/pgs.10.112.View ArticlePubMedGoogle Scholar
- Li H, Yin C, Zhang B, Sun Y, Shi L, Liu N, et al. PTTG1 promotes migration and invasion of human non-small cell lung cancer cells and is modulated by miR-186. Carcinogenesis. 2013;34(9):2145–55. doi:10.1093/carcin/bgt158.View ArticlePubMedGoogle Scholar
- Moreno-Mateos MA, Espina AG, Torres B, Gamez del Estal MM, Romero-Franco A, Rios RM, et al. PTTG1/securin modulates microtubule nucleation and cell migration. Mol Biol Cell. 2011;22(22):4302–11. doi:10.1091/mbc.E10-10-0838.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoon CH, Kim MJ, Lee H, Kim RK, Lim EJ, Yoo KC, et al. PTTG1 oncogene promotes tumor malignancy via epithelial to mesenchymal transition and expansion of cancer stem cell population. J Biol Chem. 2012;287(23):19516–27. doi:10.1074/jbc.M111.337428.PubMed CentralView ArticlePubMedGoogle Scholar
- Juanpere N, Agell L, Lorenzo M, de Muga S, López-Vilaró L, Murillo R, et al. Mutations in FGFR3 and PIK3CA, singly or combined with RAS and AKT1, are associated with AKT but not with MAPK pathway activation in urothelial bladder cancer. Hum Pathol. 2012;43(10):1573–82. doi:http://dx.doi.org/10.1016/j.humpath.2011.10.026.
- Kierska D, Szymanska H, Maslinski C. Decreased histamine content and metabolism in mammary cancer tissue from C3H mice. Agents Actions. 1992;37(3–4):227–31.View ArticlePubMedGoogle Scholar
- Knapp DW, Richardson RC, Chan TCK, Bottoms GD, Widmer WR, DeNicola DB, et al. Piroxicam Therapy in 34 Dogs With Transitional Cell Carcinoma of the Urinary Bladder. J Vet Intern Med. 1994;8(4):273–8. doi:10.1111/j.1939-1676.1994.tb03232.x.View ArticlePubMedGoogle Scholar
- Mohammed SI, Knapp DW, Bostwick DG, Foster RS, Khan KNM, Masferrer JL, et al. Expression of Cyclooxygenase-2 (COX-2) in Human Invasive Transitional Cell Carcinoma (TCC) of the Urinary Bladder. Cancer Res. 1999;59(22):5647–50.PubMedGoogle Scholar
- Sauter G, Carroll P, Moch H, Kallioniemi A, Kerschmann R, Narayan P, et al. c-myc copy number gains in bladder cancer detected by fluorescence in situ hybridization. Am J Pathol. 1995;146(5):1131–9.PubMed CentralPubMedGoogle Scholar
- Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A et al. Tissue-based map of the human proteome. Science. 2015;347(6220). doi:10.1126/science.1260419.
- Mutch DM, Berger A, Mansourian R, Rytz A, Roberts MA. The limit fold change model: a practical approach for selecting differentially expressed genes from microarray data. BMC Bioinformatics. 2002;3:17.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen R, Feng C, Xu Y. Cyclin-Dependent Kinase-Associated Protein Cks2 is Associated with Bladder Cancer Progression. J Int Med Res. 2011;39(2):533–40. doi:10.1177/147323001103900222.View ArticlePubMedGoogle Scholar