Research Papers:

Copy number variations of circulating, cell-free DNA in urothelial carcinoma of the bladder patients treated with radical cystectomy: a prospective study

PDF |  HTML  |  Supplementary Files  |  How to cite

Oncotarget. 2017; 8:56398-56407. https://doi.org/10.18632/oncotarget.17657

Metrics: PDF 2012 views  |   HTML 3106 views  |   ?  

Armin Soave, Felix K.H. Chun, Timo Hillebrand, Michael Rink, Lars Weisbach, Bettina Steinbach, Margit Fisch, Klaus Pantel and Heidi Schwarzenbach _


Armin Soave1, Felix K.-H. Chun1, Timo Hillebrand2, Michael Rink1, Lars Weisbach1, Bettina Steinbach3, Margit Fisch1, Klaus Pantel3 and Heidi Schwarzenbach3

1Department of Urology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany

2AJ Innuscreen GmbH, Berlin, Germany

3Department of Tumor Biology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany

Correspondence to:

Heidi Schwarzenbach, email: [email protected]

Keywords: urothelial carcinoma of the bladder, copy number variation, multiplex ligation-dependent probe amplification (MLPA), cell-free DNA (cfDNA), chromosomal region

Received: March 21, 2017     Accepted: April 26, 2017     Published: May 07, 2017


The aim of the present study was to establish a rapid profiling method using multiplex ligation-dependent probe amplification (MLPA) and characterize copy number variations (CNV) in circulating, cell-free DNA (cfDNA) in 85 urothelial carcinoma of the bladder (UCB) patients treated with radical cystectomy (RC). MLPA was tested for the use of cfDNA extracted from serum and plasma by various commercial extraction kits. Eighteen probes served as reference to control denaturation, ligation and amplification efficiency. MLPA was exclusively suitable for cfDNA extracted from serum. Serum from 72 patients (84.7%) could be analyzed. Thirty-five patients (48.6%) had presence of CNV in cfDNA. The median CNV count in patients with presence of CNV was 2. Predominantly, CNV were located in the genes CDH1, ZFHX3, RIPK2 and PTEN in 15 patients (20.8%), 12 patients (16.7%), 9 patients (12.5%) and 7 patients (9.7%), respectively. CNV in TSG1, RAD21, KIAA0196, ANXA7 and TMPRSS2 were associated with presence of variant UCB histology (p = 0.029, 0.029, 0.029, 0.029, 0.043, respectively). Furthermore, CNV in miR-15a, CDH1 and ZFHX3 were associated with presence of incidental prostate cancer (p = 0.023, 0.003, 0.025, respectively). Patients with CNV in KLF5, ZFHX3 and CDH1 had reduced cancer-specific survival, compared to patients without CNV in these genes (pairwise p = 0.028, 0.026, 0.044, respectively). MLPA represents an efficient method for the detection of CNV among numerous genes on various chromosomal regions. CNV in specific genes seem to be associated with aggressive UCB biologic features and presence of incidental prostate cancer, and may have a negative impact on cancer-specific survival.


Urothelial carcinoma of the bladder (UCB) is the second most common genitourinary malignancy and a fatal disease, accounting for 74,000 new cases and 16,000 deaths in the USA in 2015 [1]. For patients with muscle-invasive and recurrent high-risk non-muscle invasive UCB, radical cystectomy (RC) with bilateral pelvic lymphadenectomy represents the golden standard surgical treatment [2]. Despite major improvements in surgical techniques, imaging, perioperative management and systemic chemotherapy, outcomes have remained stable over the past decades [3], and a relevant number of patients will experience disease progression within two years after RC [4]. Various clinico-pathologic UCB features and biomarkers predicting disease recurrence and progression following RC have been introduced [5]; however, none has succeeded in daily clinical practice. Thus, there is still an urgent need for new biomarkers allowing accurate prediction of the true UCB tumor biology, helping to select patients who might best benefit from multimodal treatments and emerging targeted therapy.

Tumor cells and healthy cells release their DNA into the circulatory system. Tumor cell derived circulating cell-free DNA (cfDNA) represents an encouraging blood-based biomarker in various malignancies [68] including UCB [9]. UCB is a heterogeneous disease with complex underlying genomic alterations [10], which can be detected in cfDNA during tumor growth and disease progression [11, 12]. Real-time extraction of cfDNA from blood plasma or serum offers the promising opportunity to reveal the molecular UCB biology and course of the disease. In this regard, copy number variations (CNV) comprising DNA amplifications and deletions are a prominent source of genetic variations in cfDNA. Multiplex ligation-dependent probe amplification (MLPA) is a semi-quantitative technique for determining the relative CNV of multiple tumor suppressor genes and oncogenes in a single multiplex PCR-based analysis. To characterize CNV in cfDNA by an easy and rapid method, the aim of the present study was to establish MLPA for the use of cfDNA, together with a data analyses software custom-developed for this assay. In a single reaction, MLPA allows analyzing CNV in 43 chromosomal regions containing 37 genes.


Evaluation of the cfDNA extraction from serum and plasma

To test whether MLPA assay is suitable for the analyses of cfDNA in UCB patients, we compared the extractions of cfDNA from plasma and serum using DNA extraction kits from different companies (Qiagen, Macherey Nagel and Analytik Jena). Although cfDNA was additionally precipitated with carrier cfDNA from plasma or serum, we did not get any results on CNV by MLPA using the kits from Qiagen and Macherey Nagel. Exclusively, the PME free-circulating DNA Extraction kit (PME kit, Analytik Jena) was able to provide reliable and solid data on CNV, but only with serum. Adequate amounts of cfDNA could even be precipitated from serum by the PME kit without carrier DNA. With the NanoDrop spectrometer, we measured average concentrations of 700 ng and 100 ng in 1 ml serum and plasma, respectively, indicating that the PME kit is rather suitable for serum than plasma and the failure of the performance of MLPA using cfDNA extracted from plasma with the PME kit. Due to the supplementation of carrier DNA, the amounts of cfDNA extracted from serum or plasma with the kits from Qiagen and Macherey Nagel could not unambiguously determined. For all further experiments with the MLPA assay, we used the PME kit and at least 50 ng cfDNA (in 5 μl solution buffer) extracted from serum. Supplementary Figure 1 shows exemplary analyzable and unanalyzable electropherograms derived from capillary electrophoresis using the PME kit and the QiAmp DNA Blood Mini kit, respectively.

Clinico-pathologic UCB characteristics according to the CNV status

Using the cfDNA and genomic DNA (negative control) extracted from serum and the corresponding leukocytes of each patient, respectively, by the PME kit, we analyzed 43 chromosomal regions containing 37 genes for CNV by the MLPA assay. The CNV status was defined as the presence of DNA amplifications (DNA gains) or DNA deletions (DNA losses). As expected from normal, wild type (wt) DNA, genomic leukocyte DNA was analyzable in all 85 UCB patients and did not show any CNV, thus no genomic variations. Serum cfDNA samples were analyzable in 72 patients (85.0%) out of the 85 UCB patients. Due to the low serum DNA amounts, 13 serum samples could not be analyzed. In total, 35 out of the 72 analyzable patients (48.6%) had presence of CNV in cfDNA. The median CNV count in patients with presence of CNV was 2 (mean: 2.4; standard deviation: 1.4; range: 1–6). Most CNV were located in the genes CDH1, ZFHX3 (both copy number gains in chromosome 16), RIPK2 (copy number losses in chromosome 8) and PTEN (copy number losses in chromosome 10) in 15 patients (20.8%), 12 patients (16.7%), 9 patients (12.5%) and 7 patients (9.7%), respectively (Supplementary Table 2). Figure 1 shows a box plot of an exemplary CNV analysis in a serum sample. The CNV status was not associated with any clinico-pathologic UCB feature (Table 1). However, CNV in TSG1, RAD21, KIAA0196, ANXA7 and TMPRSS2 were associated with presence of squamous and non-squamous cell differentiation variant UCB histology (all p ≤ 0.043; Table 2). In addition, CNV in miR-15a, CDH1 and ZFHX3 were associated with presence of incidental prostate cancer in the RC specimen (all p-values ≤ 0.025; Table 2).


Figure 1: The box plot shows data of an exemplary serum sample (as calculated by Coffalyser.Net software). The DNA probes are arranged by chromosomal locations. The target-specific probes have a blue, orange and white background (left), whereas the reference probes have a grey background (right). Red point indicates a significant decrease in CNV, whereas yellow/orange points are ambiguous. A reduced copy number is clearly detected in PTEN. The data were calculated by intra- and inter-sample comparisons. Intra-sample normalization was performed by dividing the fluorescence signal of each target-specific probe by the signal of every single reference probe in this probe. The median of all these ratios of this probe is the normalization constant. Subsequently, inter-sample comparison was performed by dividing the normalization constant of each probe of this sample by the average normalization constant of all 72 reference (leukocyte) samples.

Table 1: Descriptive characteristics of 72 urothelial carcinoma of the bladder patients treated with radical cystectomy and bilateral lymphadenectomy stratified by the copy number variations status

(n =72)

CNV status
(n = 37)

CNV status
(n = 35)


Age [year; median (IQR)]

70.4 (60.5; 74.8)

70.2 (64.9; 75.9)

70.1 (56.5; 73.4)


Gender (n; %)


57 (79.2)

27 (73.0)

30 (85.7)



15 (20.8)

10 (27.0)

5 (14.3)

Clinical tumor stage (n; %)


cTa, cTis

4 (5.6)

3 (8.1)

1 (2.9)


20 (27.8)

11 (29.7)

9 (25.7)


45 (62.5)

21 (56.8)

24 (68.6)


3 (4.2)

2 (5.4)

1 (2.9)

Clinical tumor grade (n; %)



1 (1.4)

0 (0)

1 (2.9)


6 (8.3)

3 (8.1)

3 (8.6)


65 (90.3)

34 (91.9)

31 (88.6)

Intravesical chemo- and/or immunotherapy prior to RC (n; %)



56 (77.8)

28 (75.7)

28 (80.0)


16 (22.2)

9 (24.3)

7 (20.0)

Number of TURB prior to RC [median (IQR)]

1 (1; 2)

1 (1; 3)

1 (1; 2)


Days between last TURB and RC [median (IQR)]

39 (27; 61)

49 (28; 79)

36 (25; 50)


Pathologic tumor stage (n; %)


pT0, pTa, pTis

21 (29.2)

8 (21.6)

13 (37.1)


5 (6.9)

4 (10.8)

1 (2.9)


20 (27.8)

10 (27.0)

10 (28.6)


15 (20.8)

7 (18.9)

8 (22.9)


11 (15.3)

8 (21.6)

3 (8.6)

Combined tumor stage (n; %)


Localized (pT≤2)

46 (63.9)

22 (59.5)

24 (68.6)

Advanced (pT3-4)

26 (36.1)

15 (40.9)

11 (31.4)

Combined disease stage (n; %)


≤ pT2 and pN0

43 (59.7)

21 (56.8)

22 (62.9)

≥ pT3 or pN1-3

29 (40.3)

16 (43.2)

13 (37.1)

Pathologic tumor grade (n; %)



2 (2.8)

6 (16.2)

7 (20.0)


59 (81.9)

31 (83.8)

28 (80.0)

Concomitant carcinoma in situ (n; %)



50 (69.4)

26 (70.3)

24 (68.6)


22 (30.6)

11 (29.7)

11 (31.4)

Lymphovascular invasion (n; %)



53 (73.6)

26 (70.3)

27 (77.1)


19 (26.4)

11 (29.7)

8 (22.9)

Micro-vessel invasion (n; %)



65 (90.3)

33 (89.2)

32 (91.4)


7 (9.7)

4 (10.8)

3 (8.6)

Lymph node status (n; %)



54 (75.0)

25 (67.6)

29 (82.9)


18 (25.0)

12 (32.4)

6 (17.1)

Number of lymph nodes removed [median, (IQR)]

13 (9; 19)

13 (7; 20)

15 (9; 19)


Soft tissue surgical margin status (n; %)



61 (84.7)

29 (78.4)

32 (91.4)


11 (15.3)

8 (21.6)

3 (8.6)

Urothelial carcinoma histology (n; %)


Pure UCB

53 (73.6)

29 (78.4)

24 (68.6)

Presence of squamous cell differentiation

9 (12.5)

4 (10.8)

5 (14.3)

Presence of non-squamous cell differentiation

10 (13.9)

4 (10.8)

6 (17.1)

Presence of incidental prostate cancer in the RC specimen (n; %)



35 (48.6)

22 (59.5)

13 (37.1)


37 (51.4)

15 (40.5)

22 (62.9)

Adjuvant chemotherapy (n; %)


Not administered

53 (73.6)

26 (70.3)

27 (77.1)


19 (26.4)

11 (29.7)

8 (22.9)

Adjuvant chemotherapy regimen (n; %)



11 (57.9)

7 (63.6)

4 (50.0)


8 (42.1)

4 (36.4)

4 (50.0)

CNV status negative: no DNA amplifications or deletions.

CNV status positive: presence of DNA amplifications or deletions.

* = Fisher’s Exact test.

** = Pearson χ2 test.

*** = Mann-Whitney-U test.

Abbreviations: CNV = copy number variations; IQR = interquartile range; RC = radical cystectomy; TURB = transurethral resection of the bladder.

Table 2: Urothelial carcinoma histology and presence of incidental prostate cancer in the radical cystectomy specimen of 72 urothelial carcinoma of the bladder patients treated with radical cystectomy and bilateral lymphadenectomy stratified by copy number variations of selected genes


Outcomes according to the CNV status

The median follow-up of cancer survivors was 16 months (IQR: 4; 28). Actuarial two-year recurrence-free, cancer-specific and overall survival estimates were 67% ± 7% (standard error), 92% ± 4% and 88% ± 5% respectively.

In Kaplan-Meier analyses there was no difference in recurrence-free, cancer-specific and overall survival according to the CNV status in cfDNA (pairwise p = 0.409, 0.419 and 0.477, respectively; Figure not shown). However, patients with copy number gains in KLF5, ZFHX3 and CDH1 had reduced cancer-specific survival, compared to patients without CNV in these genes ZFHX3 and CDH1, respectively (pairwise p ≤ 0.044; Figure 2). In univariable logistic regression analysis, copy number gains in KLF5 (hazard ratio (HR): 3.2; 95% confidence interval (CI): 1.2–9.1; p = 0.025) and ZFHX3 (HR: 2.5; 95%CI: 1.1–6.0; p = 0.038) were associated with cancer-specific mortality.


Figure 2: Kaplan-Meier plots of cancer-specific survival stratified by CNV in ZFHX3 (A), KLF5 (B) and CDH1 (C) in 72 UCB patients treated with RC and bilateral pelvic lymphadenectomy. Top curves show UCB patients with no CNV (no genomic aberrations), and bottom curves show patients with CNV comprising DNA amplifications.


In the current study, we tested whether the MLPA assay is suitable for the analyses of cfDNA in blood of UCB patients. In contrast to the high amounts of intact genomic DNA derived from cells, plasma/serum frequently delivers insufficient amounts of cfDNA that is additionally fragmented. These limitations of getting adequate amounts of qualitatively good cfDNA and the additional prevalence of wild-type cfDNA over tumor cfDNA in blood indicates the difficulties to establish reliable multiplex genetic analyses with tumor derived cfDNA. With a specific extraction method, i.e., the PME kit and provided that at least 50 ng of cfDNA can be recovered from this kit, we found CNV in serum cfDNA from nearly each second UCB patient. The reliability of the CNV detection was examined by a repetitive analysis of some serum cfDNA samples and by the fact that genomic leukocyte DNA did not show any CNV in all 85 UCB patients.

Currently, it still remains a matter of debates, whether serum or plasma is the optimal source for cfDNA analyses [13]. However, within the framework of our study, our findings indicate that serum is only suitable for the MLPA assay. The lack of getting data with plasma is not based on the cfDNA amounts and may rather be due to anticoagulants required for the preparation of plasma, such as heparin, acid citrate dextrose (ACD) or EDTA. These agents may affect the MLPA [14]. In addition, preanalytical variables [13] may also affect this assay. However, their description is beyond the scope of the present study. For the extraction of cfDNA from plasma and serum, we used various established commercial extraction kits, which are commonly applied by numerous research groups. Notably, only the PME cfDNA Extraction kit was suitable for CNV analyses with the MLPA assay. Possibly, the components of the other extraction kits may have an influence on the performance of MLPA. In this context, it has also been demonstrated that the NucleoSpin Plasma XS Kit from Macherey-Nagel [15] is superior to the QIAamp DNA blood mini kit [1618] in terms of yield, purity and efficiency of small DNA fragment retrieval. Likewise, the NucleoSpin Plasma XS Kit did also not provide data with the MLPA assay. Therefore, further studies are needed to shed more light on the important issue of differences among various cfDNA extraction methods from serum.

To carry out the MLPA assay using serum cfDNA, it should also be considered that tumor-derived cfDNA is diluted by wild type cfDNA in the blood of cancer patients that may camouflage the detection of genetic alterations in tumor-derived cfDNA. Therefore, the establishment of a reliable detection method of CNV in cfDNA is still challenging. Although the low prevalence of tumor-derived cfDNA is a drawback of carrying out such analyses [6, 19] we could efficiently detect CNV in cfDNA. MLPA overcomes even the limitations of in situ hybridization with its limited resolution of longer than 20 kb DNA molecules which is not suitable for detection of exon-length CNV. Furthermore, the analysis of in situ hybridization assays is labor-extensive, and cannot be scaled to high-throughput and -multiplex testing [20].

In respect to our data, the majority of CNV were located in oncogenes and tumor suppressor genes, including CDH1, ZFHX3, RIPK2 and PTEN. These genes play a major role in cancer progression in various malignancies, including bladder cancer [2124]. For the first time, our findings suggest that CNV in these genes may be of particular relevance for the understanding of the UCB biology. We found that almost half of patients had presence of incidental prostate cancer in the RC specimen, which is in line to findings of previous studies [25]. The presence of incidental prostate cancer was associated with CNV in miR-15a, CDH1 and ZFHX3. We detected copy number gains in CDH1 and ZFHX3 on chromosome 16, and copy number losses in miR-15a on chromosome 13. Similarly, using comparative genomic hybridization, other authors have suggested that gains on chromosome 16 and losses on chromosome 13 may be common findings in non-metastatic and metastatic prostate cancer [26]. CDH1 is a cadherin cell adhesion molecule that is involved in epithelial-mesenchymal transition [21]. Using RNA-Seq and real-time qPCR, CDH1 was shown to be up-regulated in recurrent muscle-invasive cisplatin-resistant UCB tissue compared with adjacent non-tumor tissue [27]. Interestingly, genomic gains in the transcription factor ZFHX3 could also be found in circulating tumor cells from prostate cancer patients [28]. Moreover, the putative tumor suppressor miR-15a was described to be homozygously deleted in prostate cancer cell lines and xenografts [29]. Expression of this microRNA inhibited cell proliferation, promoted apoptosis of cancer cells, and suppressed tumorigenicity of diverse cancer types, both in vitro and in vivo [30].

Transformation from pure UCB to variant UCB histology has been suggested being a loss of differentiation [31]. Still, genetic and epigenetic modifications during the transformation from pure UCB to variant UCB histology remain mainly unresolved [31, 32]. In this regard, we found CNV in specific genes, such as TSG1, RAD21, KIAA0196, ANXA7 and TMPRSS2, to be associated with the presence of variant UCB histology at RC. RAD21 encodes a key component of the cohesin complex, which is essential for chromosome segregation, and RAD21 deregulation may impact survival in breast cancer [33]. TSG1 and KIAA0196 are poorly characterized genes. In prostate cancer, KIAA0196 is amplified and associated with poor prognosis [34]. The cancer-specific expression of ANXA7, a GTPase, has been described as a diagnostic marker of cancer and a potential target for cancer treatment. Cross talk of ANXA7 with PTEN and EGF receptor led to constitutive activation of PI3K-AKT signaling, a central pathway of tumor cell survival and proliferation [35]. To sum up, our findings suggest that genetic instability in variant histology may result in CNV in specific genes. However, further studies with larger patient cohorts with variant UCB histology are needed to verify this hypothesis.

We found that the CNV status was not associated with outcomes. However, patients with CNV in specific genes, such as ZFHX3, KLF5 and CDH1, had reduced cancer-specific survival compared to patients without CNV in these genes. Patients with copy number gains in KLF5 and ZFHX3 were at a 3.2-fold and 2.5-fold increased risk of cancer-specific mortality, respectively. Previously, it was reported that the transcription factor KLF5 is involved in tumorigenesis of UCB, by promoting cell proliferation, migration and angiogenesis [3638]. Thus, KLF5 could become a promising therapeutic target molecule for UCB. Recently, in an established xenograft mouse model of colon cancer, the drug ML264 efficiently inhibited growth of the tumor within 5 days of treatment by inhibiting the expression of KLF5 and EGR1, a transcriptional activator of KLF5 [39]. Moreover, in non-muscle invasive UCB treated with TURBT, the tumor suppressive transcription factor ZFHX3 seems to be an independent predictor for disease recurrence [22].

To our knowledge this is the first study investigating CNV in cfDNA using MLPA in UCB patients treated with RC, however, it is not devoid of some limitations. First and foremost, overall sample size and follow-up data are limited. Therefore, we cannot exclude that findings may be different in larger patient cohorts with extended follow-up. In addition, the limited number of events, that could also be caused by the prevalence of normal wild type cfDNA in blood of UCB patients, impeded multivariable analysis of risk factors for disease recurrence and survival. Nevertheless, our study remains the first and currently largest using MLPA in UCB patients treated with RC. Approximately every second patient had an incidental prostate cancer in the RC specimens. Thus, CNV in cfDNA detected in our study may be derived from prostate or bladder cancer cells, as well as from both. Further investigations on CNV in the primary tumor together with CNV in cfDNA are necessary to shed light on the source of cfDNA, whether it stems from the primary tumor, circulating tumor cells or micrometastases. A further limitation of our study is the fact that no single analyses on each chromosomal region by real-time Taqman PCR have been performed due to the insufficient serum cfDNA amounts extracted from our patient cohort. However, these analyses concerns future studies, since the main focus of the present study was to establish a quick technique without excessive statistical efforts to reliably analyze multiply CNV in cfDNA.

MLPA represents a simple and efficient method for the detection of CNV among numerous genes on various chromosomal regions. Prior to RC, approximately half of UCB patients harbor CNV in different tumor suppressor genes and oncogenes. CNV in specific genes are associated with aggressive UCB biologic features and presence of incidental prostate cancer in the RC specimen. In addition, CNV in specific genes may have a negative impact on cancer-specific survival. The inclusion of MLPA in future studies is recommended to validate our findings in larger patient cohorts.


Patient cohort

After written informed consent, we prospectively enrolled 85 UCB patients treated with RC and bilateral pelvic lymphadenectomy without neoadjuvant chemotherapy at the University Medical Center Hamburg-Eppendorf between 2011 and 2014. Recurrent Ta, T1, or carcinoma in situ (CIS), refractory to transurethral resection of the bladder (TURBT) with or without intravesical immunotherapy or chemotherapy, or muscle invasive UCB were indications for RC. Preoperative staging consisted of computed tomography (CT) of the thorax and abdomen/pelvis, and bone scan and cranium imaging when clinically indicated. Exclusion criteria included metastatic disease at preoperative staging, a history of any other malignancy, previous systemic chemotherapy or radiation, as well as incomplete clinico-pathologic or follow-up data. In total, 19 patients (26.4%) received adjuvant chemotherapy based on tumor stage, overall health status, renal function and patients’ desire. Adjuvant chemotherapy was consistently platinum-based and generally started within 90 days after RC. The study was approved by the local ethics committee (No. PV3962).

Pathological evaluation

The complete surgical RC specimen was inked, and multiple sections were obtained from the bladder and the tumor in addition to the regional lymph nodes and ureters. Tumor stage and nodal status were assessed according to the tumor, lymph node and metastasis (TNM) system. Tumor grade was assessed according to the 1998 World Health Organization (WHO) grading system [40]. Concomitant CIS was defined as the presence of CIS in conjunction with another tumor other than CIS alone. Lymphovascular invasion (LVI) was defined as the unequivocal presence of tumor cells within an endothelium-lined space without underlying muscular walls [41]. Micro-vascular invasion (MVI) was defined as the presence of tumor cells within a vessel with a vascular wall and red blood cells in the lumen [42]. A positive soft tissue surgical margin (STSM) was defined as the presence of tumor at inked areas of soft tissue on the RC specimen [43]. Presence of variant UCB histology was defined as the presence of UCB combined with any variant histology. Variant UCB histologies were classified corresponding to the WHO Classification of Tumors [44]. Incidental prostate cancer was defined as presence of prostate cancer in the RC specimens [25].

cfDNA extraction

Preoperative blood samples were usually collected on the day prior to RC at a median of 39 days [interquartile range (IQR): 27; 61] after the preceding TURB. Serum and plasma were prepared from 6 ml whole blood. cfDNA was extracted from serum and plasma using diverse DNA extraction kits (i.e., QiAmp DNA Blood Mini kit, Qiagen, Hilden, Germany; QiAmp Circulating Nucleic Acid kit, Qiagen; NucleoSpin Plasma XS kit, Macherey Nagel, Düren, Germany; PME free-circulating DNA Extraction kit, Analytik Jena, Germany). cfDNA was extracted from 2 ml serum or plasma as well as leukocytes (reference) from 6 ml EDTA blood, and performed according to the manufacturer´s instructions. Quantification and quality of the extracted cfDNA were determined spectrophotometrically using the NanoDrop Spectrometer ND-1000 (Thermo Fisher Scientific, Wilmington, DE, USA).

MLPA assay

CNV were determined using 5 μl (50 ng) cfDNA, 5 μL (100 ng) leukocyte (reference) DNA from 85 BCa patients and the SALSA MLPA probemix X049-A1 kit (MRC Holland, Amsterdam, The Netherlands). This kit contains a probe mix of 43 sequences of 37 genes to be analyzed and 22 reference genes (Supplementary Table 1). According to the manufacturer’s instructions, the MLPA probes were hybridized and ligated to denaturated serum, plasma and leukocyte DNA. During the subsequent PCR, all ligated probes were amplified simultaneously using the same PCR primer pair, of which one PCR primer was fluorescently labelled. Fragment separation was done by capillary electrophoresis on an automated ABI 3130 DNA analyzer (Applied Biosystems, Freiburg, Germany), yielding a specific electropherogram.

Data normalization

Data normalization was carried out by Coffalyser.Net analysis software (wwww.mlpa.com). It consists of 2 steps: intra- and intersample normalization. For intrasample normalization, within each sample, each probe peak was compared with the peaks of the reference probes. Reference probes located on various chromosomes detect sequences that are expected to have a normal copy number (CN) in all samples. The determined relative probe signals were then used for intersample normalization (Supplementary Table 1). Final probe ratios were determined by comparing the relative probe peak in the cfDNA sample of interest with those of all leukocyte DNA samples. Leukocyte DNA samples are expected to have a normal CN for both the reference and target probe. To avoid false positive data due to the quality and quantity of the serum cfDNA, only unambiguous values were used (Figure 1), and PCR was repeated.

Follow-up regimen

For the first year following RC, patients were generally seen every three months, from the second to fifth years every six months, and annually thereafter. Follow-up comprised a history, serum chemistry evaluation and physical examination. Diagnostic imaging of the abdomen including the urinary tract (e.g. ultrasonography and/or intravenous urography, CT of the abdomen/pelvis with intravenous contrast) and chest radiography were performed at least annually or when clinically indicated. Further radiographic evaluations (i.e., bone or brain scans, magnetic resonance imaging) were conducted at the discretion of the treating physician when clinically indicated.

Disease recurrence was defined as local failure in the operative site, regional lymph nodes, or distant metastasis. Upper tract urothelial carcinoma was not considered as disease recurrence but metachronous tumor. Cancer-specific mortality was defined as death from UCB. Overall mortality was defined as death from any cause. The cause of death was determined by the treating physician, by chart review corroborated by death certificates, or by death certificates alone [45]. Perioperative mortality (i.e., death within 30 days of surgery) was censored at time of death for bladder cancer-specific survival analyses.

Statistical analyses

The co-primary endpoints of the present study were disease recurrence, cancer-specific and overall mortality according to CNV in multiple tumor suppressor genes and oncogenes. The indicator variables (i.e., CNV) were analyzed as categorical variables. Associations between categorical variables were assessed using the Fisher exact and χ2-test. Differences in continuous variables were analyzed using the Mann-Whitney-U test (two categories) and the Kruskal-Wallis test (three or more categories). Recurrence-free, cancer-specific and overall survival probabilities were estimated using the Kaplan-Meier method and differences between groups were assessed using the Log rank statistic. Univariable Cox regression models assessed time to disease recurrence, cancer-specific and overall mortality. All tests are two-sided and a p-value of < 0.05 was set to be statistically significant. All analyses were performed with SPSS 20 (SPSS Inc., IBM Corp., Armonk, NY).


We thank Dr. Lan Kluwe for her assistance with the MLPA data evaluation by the Coffalyser.Net analysis software.


T.H., director of J Innuscreen GmbH, provided the PME free-circulating DNA Extraction kits free of charge. H.S. received the SALSA MLPA probemix X049-A1 kits from MRC Holland, Amsterdam, The Netherlands, free of charge. No potential conflicts of interest were disclosed by the other authors.


1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2015. CA Cancer J Clin. 2015; 65:5–29.

2. Witjes JA, Compérat E, Cowan NC, De Santis M, Gakis G, Lebret T, Ribal MJ, Van der Heijden AG, Sherif A, and European Association of Urology. EAU guidelines on muscle-invasive and metastatic bladder cancer: summary of the 2013 guidelines. Eur Urol. 2014; 65:778–92.

3. Stein JP, Lieskovsky G, Cote R, Groshen S, Feng AC, Boyd S, Skinner E, Bochner B, Thangathurai D, Mikhail M, Raghavan D, Skinner DG. Radical cystectomy in the treatment of invasive bladder cancer: long-term results in 1,054 patients. J Clin Oncol. 2001; 19:666–75.

4. Rink M, Lee DJ, Kent M, Xylinas E, Fritsche HM, Babjuk M, Brisuda A, Hansen J, Green DA, Aziz A, Cha EK, Novara G, Chun FK, et al, and Bladder Cancer Research Consortium. Predictors of cancer-specific mortality after disease recurrence following radical cystectomy. BJU Int. 2013 (3 Pt B); 111:E30–36.

5. Rink M, Cha EK, Green D, Hansen J, Robinson BD, Lotan Y, Sagalowsky AI, Chun FK, Karakiewicz PI, Fisch M, Scherr DS, Shariat SF. Biomolecular predictors of urothelial cancer behavior and treatment outcomes. Curr Urol Rep. 2012; 13:122–35.

6. Schwarzenbach H, Hoon DS, Pantel K. Cell-free nucleic acids as biomarkers in cancer patients. Nat Rev Cancer. 2011; 11:426–37.

7. Schwarzenbach H, Chun FK, Isbarn H, Huland H, Pantel K. Genomic profiling of cell-free DNA in blood and bone marrow of prostate cancer patients. J Cancer Res Clin Oncol. 2011; 137:811–19.

8. Schwarzenbach H. Circulating nucleic acids and protease activities in blood of tumor patients. Expert Opin Biol Ther. 2012 (Suppl 1); 12:S163–69.

9. Birkenkamp-Demtröder K, Nordentoft I, Christensen E, Høyer S, Reinert T, Vang S, Borre M, Agerbæk M, Jensen JB, Ørntoft TF, Dyrskjøt L. Genomic Alterations in Liquid Biopsies from Patients with Bladder Cancer. Eur Urol. 2016; 70:75–82.

10. Cancer Genome Atlas Research Network. Comprehensive molecular characterization of urothelial bladder carcinoma. Nature. 2014; 507:315–22.

11. Ralla B, Stephan C, Meller S, Dietrich D, Kristiansen G, Jung K. Nucleic acid-based biomarkers in body fluids of patients with urologic malignancies. Crit Rev Clin Lab Sci. 2014; 51:200–31.

12. Abbosh PH, Rosenberg JE, Plimack ER. Circulating biomarkers to guide systemic therapy for urothelial carcinoma. Urol Oncol. 2016; 34:502–09.

13. Bronkhorst AJ, Aucamp J, Pretorius PJ. Cell-free DNA: preanalytical variables. Clin Chim Acta. 2015; 450:243–53.

14. Xue X, Teare MD, Holen I, Zhu YM, Woll PJ. Optimizing the yield and utility of circulating cell-free DNA from plasma and serum. Clin Chim Acta. 2009; 404:100–04.

15. Kirsch C, Weickmann S, Schmidt B, Fleischhacker M. An improved method for the isolation of free-circulating plasma DNA and cell-free DNA from other body fluids. Ann N Y Acad Sci. 2008; 1137:135–39.

16. Legler TJ, Liu Z, Mavrou A, Finning K, Hromadnikova I, Galbiati S, Meaney C, Hultén MA, Crea F, Olsson ML, Maddocks DG, Huang D, Fisher SA, et al. Workshop report on the extraction of foetal DNA from maternal plasma. Prenat Diagn. 2007; 27:824–29.

17. Su YH, Wang M, Block TM, Landt O, Botezatu I, Serdyuk O, Lichtenstein A, Melkonyan H, Tomei LD, Umansky S. Transrenal DNA as a diagnostic tool: important technical notes. Ann N Y Acad Sci. 2004; 1022:81–89.

18. Müller I, Beeger C, Alix-Panabières C, Rebillard X, Pantel K, Schwarzenbach H. Identification of loss of heterozygosity on circulating free DNA in peripheral blood of prostate cancer patients: potential and technical improvements. Clin Chem. 2008; 54:688–96.

19. Diaz LA Jr, Bardelli A. Liquid biopsies: genotyping circulating tumor DNA. J Clin Oncol. 2014; 32:579–86.

20. Kondrashova O, Love CJ, Lunke S, Hsu AL, Waring PM, Taylor GR, and Australian Ovarian Cancer Study (AOCS) Group. High-Throughput Amplicon-Based Copy Number Detection of 11 Genes in Formalin-Fixed Paraffin-Embedded Ovarian Tumour Samples by MLPA-Seq. PLoS One. 2015; 10:e0143006.

21. Bryan RT, Tselepis C. Cadherin switching and bladder cancer. J Urol. 2010; 184:423–31.

22. Kawaguchi M, Hara N, Bilim V, Koike H, Suzuki M, Kim TS, Gao N, Dong Y, Zhang S, Fujinawa Y, Yamamoto O, Ito H, Tomita Y, et al. A diagnostic marker for superficial urothelial bladder carcinoma: lack of nuclear ATBF1 (ZFHX3) by immunohistochemistry suggests malignant progression. BMC Cancer. 2016; 16:805.

23. Guirado M, Gil H, Saenz-Lopez P, Reinboth J, Garrido F, Cozar JM, Ruiz-Cabello F, Carretero R. Association between C13ORF31, NOD2, RIPK2 and TLR10 polymorphisms and urothelial bladder cancer. Hum Immunol. 2012; 73:668–72.

24. Sfakianos JP, Lin Gellert L, Maschino A, Gotto GT, Kim PH, Al-Ahmadie H, Bochner BH. The role of PTEN tumor suppressor pathway staining in carcinoma in situ of the bladder. Urol Oncol. 2014; 32:657–62.

25. Bruins HM, Djaladat H, Ahmadi H, Sherrod A, Cai J, Miranda G, Skinner EC, Daneshmand S. Incidental prostate cancer in patients with bladder urothelial carcinoma: comprehensive analysis of 1,476 radical cystoprostatectomy specimens. J Urol. 2013; 190:1704–09.

26. Zitzelsberger H, Engert D, Walch A, Kulka U, Aubele M, Höfler H, Bauchinger M, Werner M. Chromosomal changes during development and progression of prostate adenocarcinomas. Br J Cancer. 2001; 84:202–08.

27. Zhang S, Liu Y, Liu Z, Zhang C, Cao H, Ye Y, Wang S, Zhang Y, Xiao S, Yang P, Li J, Bai Z. Transcriptome profiling of a multiple recurrent muscle-invasive urothelial carcinoma of the bladder by deep sequencing. PLoS One. 2014; 9:e91466.

28. Gupta S, Li J, Kemeny G, Bitting RL, Beaver J, Somarelli J, Ware KE, Gregory S, Armstrong AJ. Whole genomic copy number alterations in circulating tumor cells from men with abiraterone or enzalutamide resistant metastatic castration-resistant prostate cancer. Clin Cancer Res. 2017; 23:1346–57.

29. Porkka KP, Ogg EL, Saramäki OR, Vessella RL, Pukkila H, Lähdesmäki H, van Weerden WM, Wolf M, Kallioniemi OP, Jenster G, Visakorpi T. The miR-15a-miR-16-1 locus is homozygously deleted in a subset of prostate cancers. Genes Chromosomes Cancer. 2011; 50:499–509.

30. Huang E, Liu R, Chu Y. miRNA-15a/16: as tumor suppressors and more. Future Oncol. 2015; 11:2351–63.

31. Soave A, Schmidt S, Dahlem R, Minner S, Engel O, Kluth LA, John LM, Hansen J, Schmid M, Sauter G, Shariat SF, Fisch M, Rink M. Does the extent of variant histology affect oncological outcomes in patients with urothelial carcinoma of the bladder treated with radical cystectomy? Urol Oncol. 2015; 33:21.e1–21.e9.

32. Xylinas E, Rink M, Robinson BD, Lotan Y, Babjuk M, Brisuda A, Green DA, Kluth LA, Pycha A, Fradet Y, Faison T, Lee RK, Karakiewicz PI, et al. Impact of histological variants on oncological outcomes of patients with urothelial carcinoma of the bladder treated with radical cystectomy. Eur J Cancer. 2013; 49:1889–97.

33. Yan M, Xu H, Waddell N, Shield-Artin K, Haviv I, McKay MJ, Fox SB, and kConFab authors. Enhanced RAD21 cohesin expression confers poor prognosis in BRCA2 and BRCAX, but not BRCA1 familial breast cancers. Breast Cancer Res. 2012; 14:R69.

34. Jalava SE, Porkka KP, Rauhala HE, Isotalo J, Tammela TL, Visakorpi T. TCEB1 promotes invasion of prostate cancer cells. Int J Cancer. 2009; 124:95–102.

35. Leighton X, Eidelman O, Jozwik C, Pollard HB, Srivastava M. ANXA7-GTPase as Tumor Suppressor: Mechanisms and Therapeutic Opportunities. Methods Mol Biol. 2017; 1513:23–35.

36. Du C, Gao Y, Xu S, Jia J, Huang Z, Fan J, Wang X, He D, Guo P. KLF5 promotes cell migration by up-regulating FYN in bladder cancer cells. FEBS Lett. 2016; 590:408–18.

37. Gao Y, Wu K, Chen Y, Zhou J, Du C, Shi Q, Xu S, Jia J, Tang X, Li F, Hui K, He D, Guo P. Beyond proliferation: KLF5 promotes angiogenesis of bladder cancer through directly regulating VEGFA transcription. Oncotarget. 2015; 6:43791–805. https://doi.org/10.18632/oncotarget.6101.

38. Chen C, Benjamin MS, Sun X, Otto KB, Guo P, Dong XY, Bao Y, Zhou Z, Cheng X, Simons JW, Dong JT. KLF5 promotes cell proliferation and tumorigenesis through gene regulation and the TSU-Pr1 human bladder cancer cell line. Int J Cancer. 2006; 118:1346–1355.

39. Ruiz de Sabando A, Wang C, He Y, García-Barros M, Kim J, Shroyer KR, Bannister TD, Yang VW, Bialkowska AB. ML264, A Novel Small-Molecule Compound That Potently Inhibits Growth of Colorectal Cancer. Mol Cancer Ther. 2016; 15:72–83.

40. Epstein JI, Amin MB, Reuter VR, Mostofi FK. The World Health Organization/International Society of Urological Pathology consensus classification of urothelial (transitional cell) neoplasms of the urinary bladder. Bladder Consensus Conference Committee. Am J Surg Pathol. 1998; 22:1435–48.

41. Lotan Y, Gupta A, Shariat SF, Palapattu GS, Vazina A, Karakiewicz PI, Bastian PJ, Rogers CG, Amiel G, Perotte P, Schoenberg MP, Lerner SP, Sagalowsky AI. Lymphovascular invasion is independently associated with overall survival, cause-specific survival, and local and distant recurrence in patients with negative lymph nodes at radical cystectomy. J Clin Oncol. 2005; 23:6533–39.

42. Hong SK, Kwak C, Jeon HG, Lee E, Lee SE. Do vascular, lymphatic, and perineural invasion have prognostic implications for bladder cancer after radical cystectomy? Urology. 2005; 65:697–702.

43. Novara G, Svatek RS, Karakiewicz PI, Skinner E, Ficarra V, Fradet Y, Lotan Y, Isbarn H, Capitanio U, Bastian PJ, Kassouf W, Fritsche HM, Izawa JI, et al. Soft tissue surgical margin status is a powerful predictor of outcomes after radical cystectomy: a multicenter study of more than 4,400 patients. J Urol. 2010; 183:2165–70.

44. Eble J, Sauter G, Epstein JI, Sesterhenn IA (editors). Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs. World Health Organization Classification of Tumours. IARCPress. Lyon. 2004.

45. Rink M, Fajkovic H, Cha EK, Gupta A, Karakiewicz PI, Chun FK, Lotan Y, Shariat SF. Death certificates are valid for the determination of cause of death in patients with upper and lower tract urothelial carcinoma. Eur Urol. 2012; 61:854–55.

Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 17657