Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

Network signatures based on gene pair expression ratios improve classification and the analysis of muscle-invasive urothelial cancer

Urothelial cancer (UC) is highly recurrent and can progress from non-invasive (NMIUC) to a more aggressive muscle-invasive (MIUC) subtype that invades the muscle tissue layer of the bladder. We present a proof of principle study that network-based features of gene pairs can be used to improve classifier performance and the functional analysis of urothelial cancer gene expression data. In the first step of our procedure each individual sample of a UC gene expression dataset is inflated by gene pair expression ratios that are defined based on a given network structure. In the second step an elastic net feature selection procedure for network-based signatures is applied to discriminate between NMIUC and MIUC samples. We performed a repeated random subsampling cross validation in three independent datasets. The network signatures were characterized by a functional enrichment analysis and studied for the enrichment of known cancer genes. We observed that the network-based gene signatures from meta collections of proteinprotein interaction (PPI) databases such as CPDB and the PPI databases HPRD and BioGrid improved the classification performance compared to single gene based signatures. The network based signatures that were derived from PPI databases showed a prominent enrichment of cancer genes (e.g., TP53, TRIM27 and HNRNPA2Bl). We provide a novel integrative approach for large-scale gene expression analysis for the identification and development of novel diagnostical targets in bladder cancer. Further, our method allowed to link cancer gene associations to network-based expression signatures that are not observed in gene-based expression signatures.

Bcl-2 predicts response to neoadjuvant chemotherapy and is overexpressed in lymph node metastases of urothelial cancer of the bladder

PURPOSE To assess whether Bcl-2, an inhibitor of the apoptotic cascade, can predict response to neoadjuvant chemotherapy in patients with urothelial cancer of the bladder (UCB). METHODS Bcl-2 expression was analyzed in 2 different tissue microarrays (TMAs). One TMA was constructed of primary tumors and their corresponding lymph node (LN) metastases from 152 patients with chemotherapy-naive UCB treated by cystectomy and pelvic lymphadenectomy (chemotherapy-naive TMA cohort). The other TMA was constructed of tumor samples obtained from 55 patients with UCB before neoadjuvant chemotherapy (transurethral resection of the bladder cancer) and after cystectomy with pelvic lymphadenectomy (residual primary tumor [ypT+], n = 38); residual LN metastases [ypN+], n = 24) (prechemotherapy/postchemotherapy TMA cohort). Bcl-2 overexpression was defined as 10% or more cancer cells showing cytoplasmic immunoreactivity. RESULTS In both TMA cohorts, Bcl-2 overexpression was significantly (P<0.05) more frequent in LN metastases than in primary tumors (chemotherapy-naive TMA group: 18/148 [12%] in primary tumors vs. 39/143 [27%] in metastases; postchemotherapy TMA: ypT+7/35 [20%] vs. ypN+11/19 [58%]). In the neoadjuvant setting, patients with Bcl-2 overexpression in transurethral resection of the bladder cancer specimens showed significantly (P = 0.04) higher ypT stages and less regression in their cystectomy specimens than did the control group, and only one-eighth (13%) had complete tumor regression (ypT0 ypN0). In survival analyses, only histopathological parameters added significant prognostic information. CONCLUSIONS Bcl-2 overexpression in chemotherapy-naive primary bladder cancer is related to poor chemotherapy response and might help to select likely nonresponders.

The impact of biomarkers in multivariate algorithms for bladder cancer diagnosis in patients with hematuria

BACKGROUND: We appraised 23 biomarkers previously associated with urothelial cancer in a case-control study. Our aim was to determine whether single biomarkers and/or multivariate algorithms significantly improved on the predictive power of an algorithm based on demographics for prediction of urothelial cancer in patients presenting with hematuria. METHODS: Twenty-two biomarkers in urine and carcinoembryonic antigen (CEA) in serum were evaluated using enzyme-linked immunosorbent assays (ELISAs) and biochip array technology in 2 patient cohorts: 80 patients with urothelial cancer, and 77 controls with confounding pathologies. We used Forward Wald binary logistic regression analyses to create algorithms based on demographic variables designated prior predicted probability (PPP) and multivariate algorithms, which included PPP as a single variable. Areas under the curve (AUC) were determined after receiver-operator characteristic (ROC) analysis for single biomarkers and algorithms. RESULTS: After univariate analysis, 9 biomarkers were differentially expressed (t test; P

Comprehensive profiling and localisation of the matrix metalloproteinases in urothelial carcinoma

The matrix metalloproteinases (MMPs) are endopeptidases which break down the extracellular matrix and regulate cytokine and growth factor activity. Several MMPs have been implicated in the promotion of invasion and metastasis in a broad range of tumours including urothelial carcinoma. In this study, RNA from 132 normal bladder and urothelial carcinoma specimens was profiled for each of the 24 human MMPs, the four endogenous tissue inhibitors of MMPs (TIMPs) and several key growth factors and their receptors using quantitative real time RT-PCR. Laser capture microdissection (LCM) of RNA from 22 tumour and 11 normal frozen sections was performed allowing accurate RNA extraction from either stromal or epithelial compartments. This study confirms the over expression in bladder tumour tissue of well-documented MMPs and highlights a range of MMPs which have not previously been implicated in the development of urothelial cancer. In summary, MMP-2, MT1-MMP and the previously unreported MMP-28 were very highly expressed in tumour samples while MMPs 1, 7, 9, 11, 15, 19 and 23 were highly expressed. There was a significant positive correlation between transcript expression and tumour grade for MMPs 1, 2, 8, 10, 11, 12, 13, 14, 15 and 28 (P < 0.001). At the same confidence interval, TIMP-1 and TIMP-3 also correlated with increasing tumour grade. LCM revealed that most highly expressed MMPs are located primarily within the stromal compartment except MMP-13 which localised to the epithelial compartment. This work forms the basis for further functional studies, which will help to confirm the MMPs as potential diagnostic and therapeutic targets in early bladder cancer.

MMP-2 and MMP-9 in lymph-node-positive bladder cancer

Matrix metalloproteinases (MMP), particularly MMP-2 and MMP-9, participate in tumour progression and metastasis in various cancers. Their significance in urothelial cancer of the bladder (UCB) is unclear. Expression analysis of MMP-2 and MMP-9 in tissue microarrays (TMA) constructed of corresponding samples from histopathological normal urothelium, tumour centre and invasion front of primary tumours and lymph-node (LN) metastases might help to elucidate their relevance in UCB.

Occupational and non-occupational risk factors in bladder cancer patients in an industrialized area located in former east-Germany

Purpose: Several occupational carcinogens are metabolized by polymorphic enzymes. The distribution of the polymorphic enzymes N-acetyltransferase 2 (NAT2; substrates: aromatic amines), glutathione S-transferase M1 (GSTM1; substrates: e.g., reactive metabolites of polycyclic aromatic hydrocarbons), and glutathione S-transferase T1 (GSTT1; substrates: small molecules with 1 - 2 carbon atoms) were investigated. Material and Methods: At the urological department in Lutherstadt Wittenberg, 136 patients with a histologically proven transitional cell cancer of the urinary bladder were investigated for all occupations performed for more than 6 months. Several occupational and non-occupational risk factors were asked. The genotypes of NAT2, GSTM1, and GSTT1 were determined from leucocyte DNA by PCR. Results: Compared to the general population in Middle Europe, the percentage of GSTT1 negative persons (22.1%) was ordinary; the percentage of slow acetylators (59.6%) was in the upper normal range, while the percentage of GSTM1 negative persons (58.8%) was elevated in the entire group. Shifts in the distribution of the genotypes were observed in subgroups who had been exposed to asbestos (6/6 GSTM1 negative, 5/6 slow acetylators), rubber manufacturing (8/10 GSTM1 negative), and chlorinated solvents (9/15 GSTM1 negative). Conclusions: The overrepresentation of GSTM1 negative bladder cancer patients also in this industrialized area and more pronounced in several occupationally exposed subgroups points to an impact of the GSTM1 negative genotype in bladder carcinogenesis.

The null allele of GSTM1 does not affect susceptibility to solar keratoses in the Australian white population

Solar keratoses (SKs) are induced by exposure to UV radiation and are capable of undergoing transformation to squamous cell carcinoma (SCC).1 The two main factors influencing the occurrence of SK are the sensitivity of the skin to sunlight and the total duration of solar exposure. These factors are responsible for the high incidence of SK in Australia. Although the influence of genetic factors is not defined, there is evidence that the gene encoding the enzyme, glutathione S-transferase, may be implicated in cancer predisposition and therefore SK. Glutathione S-transferase Mu-1 (GSTM1) is an isoenzyme involved in the detoxification of carcinogens. The GSTM1 protein is completely absent in approximately 50% of white persons. This absence is caused by a homozygous gene deletion on chromosome 1p resulting in a null genotype.2 Katoh3 showed that the frequency of the GSTM1 null genotype was significantly higher in 85 patients with urothelial cancer (61.2%; p < 0.05), suggesting that the null genotype may increase cancer susceptibility. This finding was supported by Lafuente et al.4 who found evidence that persons who lack the GSTM1 gene have approximately twice the chance of experiencing malignant melanoma. Further research in the United Kingdom found that patients with two or more skin tumors of different types, basal cell carcinoma (BCC) and SCC, had a significantly higher frequency of GSTM1 null genotypes than controls (71%; p = 0.033). However the GSTM1 genotype in patients with only SCC was not excessive in this population.5 Persons residing in northern Australia have the highest incidence of nonmelanoma skin cancer (SCC and BCC) in the world6 and receive far greater solar exposure than persons residing in the United Kingdom. It is possible that the GSTM1 null genotype may affect susceptibility to SK, which may act as SCC precursors, in Australians exposed to these high levels of solar radiation.

Occurrence of urinary tract tumors in miners highly exposed to dinitrotoluene

Between 1984 and 1997, six cases of urothelial cancer and 14 cases of renal cell cancer occurred in a group of 500 underground mining workers in the copper-mining industry of the former German Democratic Republic, with high exposures to explosives containing technical dinitrotoluene. Exposure durations ranged from 7 to 37 years, and latency periods ranged from 21 to 46 years. The incidences of both urothelial and renal cell tumors in this group were much higher than anticipated on the basis of the cancer registers of the German Democratic Republic by factors of 4.5 and 14.3, respectively. The cancer cases and a representative group of 183 formerly dinitrotoluene- exposed miners of this local industry were interviewed for their working history and grouped into four exposure categories. This categorization of the 14 renal cell tumor cases revealed no dose-dependency concerning explosives in any of the four exposure categories and was similar to that of the representative group of employees, whereas the urothelial tumor cases were predominantly confined to the high-exposure categories. Furthermore, all identified tumor patients were genotyped by polymerase chain reaction, using lymphocyte DNA, regarding their genetic status of the polymorphic xenobiotic metabolizing enzymes, including the N-acetyltransferase 2 and the glutathione-S-transferases M1 and T1. This genotyping revealed remarkable distributions only for the urothelial tumor cases, who were exclusively identified as 'slow acetylators.' This points to the possibility of human carcinogenicity of dinitrotoluene, with regard to the urothelium as the target tissue.

Nephrotoxicity and nephrocarcinogenicity of dinitrotoluene : new aspects to be considered

Technical dinitrotoluene (DNT) is a mixture of 2,4- and 2,6-DNT. In humans, industrial or environmental exposure can occur orally, by inhalation, or by skin contact. The classification of DNT as an 'animal carcinogen' is based on the formation of malignant tumors in kidneys, liver, and mammary glands of rats and mice. Clear signs of toxic nephropathy were found in rats dosed with DNT, and the concept was derived of an interrelation between renal toxicity and carcinogenicity. Recent data point to the carcinogenicity of DNT on the urinary tract of exposed humans. Between 1984 and 1997, 6 cases of urothelial cancer and 14 cases of renal cell cancer were diagnosed in a group of 500 underground mining workers in the copper mining industry of the former GDR and having high exposures to explosives containing technical DNT. The incidences of both urothelial and renal cell tumors in this group were 4.5 and 14.3 times higher, respectively, than anticipated on the basis of the cancer registers of the GDR. The genotyping of all identified tumor patients for the polymorphic enzymes NAT2, GSTM1, and GSTT1 identified the urothelial tumor cases as exclusively 'slow acetylates'. A group of 161 miners highly exposed to DNT was investigated for signs of subclinical renal damage. The exposures were categorized semi-quantitatively into 'low', 'medium', 'high', and 'very high'. A straight dose-dependence of the excretion of urinary biomarker proteins with the ranking of exposure was seen. Biomarker excretion (alpha1-microglobulin, glutathione S-transferases alpha and pi) indicated that DNT-induced damage was directed toward the tubular system. New data on DNT-exposed humans appear consistent with the concept of cancer initiation by DNT isomers and the subsequent promotion of renal carcinogenesis by selective damage to the proximal tubule. The differential pathways of metabolic activation of DNT appear to apply to the proximal tubule of the kidney and to the urothelium of the renal pelvis and lower urinary tract as target tissues of carcinogenicity.

Collectives of diagnostic biomarkers identify high-risk subpopulations of hematuria patients: exploiting heterogeneity in large-scale biomarker data

Background: Ineffective risk stratification can delay diagnosis of serious disease in patients with hematuria. We applied a systems biology approach to analyze clinical, demographic and biomarker measurements (n = 29) collected from 157 hematuric patients: 80 urothelial cancer (UC) and 77 controls with confounding pathologies.

Methods: On the basis of biomarkers, we conducted agglomerative hierarchical clustering to identify patient and biomarker clusters. We then explored the relationship between the patient clusters and clinical characteristics using Chi-square analyses. We determined classification errors and areas under the receiver operating curve of Random Forest Classifiers (RFC) for patient subpopulations using the biomarker clusters to reduce the dimensionality of the data.

Results: Agglomerative clustering identified five patient clusters and seven biomarker clusters. Final diagnoses categories were non-randomly distributed across the five patient clusters. In addition, two of the patient clusters were enriched with patients with ‘low cancer-risk’ characteristics. The biomarkers which contributed to the diagnostic classifiers for these two patient clusters were similar. In contrast, three of the patient clusters were significantly enriched with patients harboring ‘high cancer-risk” characteristics including proteinuria, aggressive pathological stage and grade, and malignant cytology. Patients in these three clusters included controls, that is, patients with other serious disease and patients with cancers other than UC. Biomarkers which contributed to the diagnostic classifiers for the largest ‘high cancer- risk’ cluster were different than those contributing to the classifiers for the ‘low cancer-risk’ clusters. Biomarkers which contributed to subpopulations that were split according to smoking status, gender and medication were different.

Conclusions: The systems biology approach applied in this study allowed the hematuric patients to cluster naturally on the basis of the heterogeneity within their biomarker data, into five distinct risk subpopulations. Our findings highlight an approach with the promise to unlock the potential of biomarkers. This will be especially valuable in the field of diagnostic bladder cancer where biomarkers are urgently required. Clinicians could interpret risk classification scores in the context of clinical parameters at the time of triage. This could reduce cystoscopies and enable priority diagnosis of aggressive diseases, leading to improved patient outcomes at reduced costs. © 2013 Emmert-Streib et al; licensee BioMed Central Ltd.

Urinary thrombomodulin levels were significantly higher following occupational exposure to chemicals, in the presence of dipstick protein, but not in the presence of dipstick blood

Currently, there are no biomarkers which can identify patients with an increased risk of developing urothelial cancer as a result of occupational chemical exposure. The aim of this study was to evaluate the relationships between final diagnosis and 22 biomarkers measured in urine, serum and plasma collected from 156 hematuric patients. Fourteen of the 80 patients (17.5%) with urothelial cancer and 13/76 (17.1%) of the controls were deemed to have a history of chemical exposure. We applied Fisher's exact tests to explore associations between chemical exposure and final diagnosis, and tumor stage and grade, where applicable; ANOVA and t-test to compare age across patients with and without chemical exposure; and Zelen's exact test to evaluate relationships across final diagnosis, chemical exposure and smoking. Following pre-selection of biomarkers using Lasso, we identified biomarkers with differential levels across patients with and without chemical exposure using Welch's t-test. Using a one-sided t-test and considering multiple testing using FDR, we observed that TM levels in urine were significantly higher in samples from patients with a history of chemical exposure regardless of their diagnosis as control or urothelial cancer (one-sided t-test, p_UC = 0.014 and p_CTL = 0.043); in the presence of dipstick protein and when urinary pH levels ≤ 6 (p = 0.003), but not in the presence of dipstick blood (p = 0.115). Urothelial cancer patients with a history of chemical exposure were significantly younger (64.1 years) than those without chemical exposure (70.2 years) (one-sided t-test p-value = 0.012); and their tumors were higher grade (Fisher's exact test; p = 0.008). There was a strong association between a history of chemical exposure and smoking in urothelial cancer patients (Zelen's exact test; p = 0.025). Elevated urinary thrombomodulin levels could have the potential to identify chemical exposure in hematuric patients at high risk of developing urothelial cancer.

Cytostatic in vitro Effects of DTCM-Glutarimide on Bladder Carcinoma Cells

Bladder cancer is a common malignancy worldwide. Despite the increased use of cisplatin-based combination therapy, the outcomes for patients with advanced disease remain poor. Recently, altered activation of the PI3K/Akt/mTOR pathway has been associated with reduced patient survival and advanced stage of bladder cancer, making its upstream or downstream components attractive targets for therapeutic intervention. In the present study, we showed that treatment with DTCM-glutaramide, a piperidine that targets PDK1, results in reduced proliferation, diminished cell migration and G1 arrest in 5637 and T24 bladder carcinoma cells. Conversely, no apoptosis, necrosis or autophagy were detected after treatment, suggesting that reduced cell numbers in vitro are a result of diminished proliferation rather than cell death. Furthermore previous exposure to 10 mu g/ml DTCM-glutarimide sensitized both cell lines to ionizing radiation. Although more studies are needed to corroborate our findings, our results indicate that PDK1 may be useful as a therapeutic target to prevent progression and abnormal tissue dissemination of urothelial carcinomas.