UroMark—a urinary biomarker assay for the detection of bladder cancer
- Andrew Feber1Email authorView ORCID ID profile,
- Pawan Dhami1,
- Liqin Dong1,
- Patricia de Winter2,
- Wei Shen Tan2,
- Mónica Martínez-Fernández3, 10,
- Dirk S. Paul1,
- Antony Hynes-Allen2,
- Sheida Rezaee1,
- Pratik Gurung2, 5,
- Simon Rodney2,
- Ahmed Mehmood2,
- Felipe Villacampa9, 10,
- Federico de la Rosa9, 10,
- Charles Jameson4,
- Kar Keung Cheng6,
- Maurice P. Zeegers6, 7,
- Richard T. Bryan5,
- Nicholas D. James8,
- Jesus M. Paramio3, 10,
- Alex Freeman4,
- Stephan Beck1 and
- John D. Kelly1, 2
© The Author(s). 2017
Received: 2 November 2016
Accepted: 12 December 2016
Published: 31 January 2017
Bladder cancer (BC) is one of the most common cancers in the western world and ranks as the most expensive to manage, due to the need for cystoscopic examination. BC shows frequent changes in DNA methylation, and several studies have shown the potential utility of urinary biomarkers by detecting epigenetic alterations in voided urine. The aim of this study is to develop a targeted bisulfite next-generation sequencing assay to diagnose BC from urine with high sensitivity and specificity.
We defined a 150 CpG loci biomarker panel from a cohort of 86 muscle-invasive bladder cancers and 30 normal urothelium. Based on this panel, we developed the UroMark assay, a next-generation bisulphite sequencing assay and analysis pipeline for the detection of bladder cancer from urinary sediment DNA. The 150 loci UroMark assay was validated in an independent cohort (n = 274, non-cancer (n = 167) and bladder cancer (n = 107)) voided urine samples with an AUC of 97%. The UroMark classifier sensitivity of 98%, specificity of 97% and NPV of 97% for the detection of primary BC was compared to non-BC urine.
Epigenetic urinary biomarkers for detection of BC have the potential to revolutionise the management of this disease. In this proof of concept study, we show the development and utility of a novel high-throughput, next-generation sequencing-based biomarker for the detection of BC-specific epigenetic alterations in urine.
KeywordsBladder cancer Epigenetics Urine Next generation sequencing Diagnostic
Bladder cancer represents one of the most common malignancies in the western world, ranking 8th in incidence and ranks 13th in terms of cancer mortality worldwide . Cystoscopy is the gold standard test for the detection of bladder cancer, but it is operator-dependent with a sensitivity of 90–97% [2–6]. However, cystoscopy is an invasive procedure requiring clinic or hospital attendance and poses a small but significant risk of infection . Although the sensitivity of cystoscopy is less than absolute, patient perception of the test is such that the performance of an alternative (such as a non-invasive) test should have a sensitivity of 95% or greater .
Bladder cancer carries a significant health economic burden in the UK, with the management of bladder cancer costing in excess of £55M/year [9, 10]. A significant proportion of that cost is due to the need for cystoscopy to rule out the presence of cancer. Over 110,000 cystoscopies are performed each year in the UK for patients presenting with haematuria, and a similar number for surveillance cystoscopies are performed for known non-muscle-invasive bladder cancer patients. However, given that only 10% of haematuria patients undergoing cystoscopy will have a diagnosis of bladder cancer, a non-invasive assay which can rule out the presence of cancer with a high degree of certainty will not only reduce the economic burden of cystoscopy, but also minimise the requirement for this invasive procedure in the majority of patients without cancer [9, 10]. Although several commercial assays have FDA approval for the detection of bladder cancer from urine, none are approved as stand-alone tests to replace cystoscopy . This includes urine cytology, which is frequently used as a diagnostic aid in conjunction with cystoscopy, but as with the other commercial assays, has low sensitivity to detect cancer other than high-grade disease, and carcinoma in situ thus cannot replace cystoscopy [12, 13].
Changes in DNA methylation play a key role in malignant transformation, leading to the silencing of tumour suppressor genes and overexpression of oncogenes . Despite its plasticity, DNA methylation is ontogenically relatively stable, a property which can be exploited to develop diagnostic assays resulting in an active area of research in the field of urinary-based biomarkers for the non-invasive detection of bladder cancer [15–18]. To date, DNA methylation biomarker panels have contained a relatively small number of loci, in part due to technological limitations and the requirement to retain diagnostic specificity [11–19]. Although these panels have shown promise , they have in general not reached the required sensitivity to replace cystoscopy with an inherent weakness being the limited number of targets included to maintain specificity [19–24]. More recently, assays based on mutation and methylation targets have shown high sensitivity but remain to be validated [20, 25].
Emerging techniques that utilise next-generation DNA sequencing (NGS) hold particular promise for the development of highly sensitive epigenetic biomarker panels. For example, the microdroplet-based PCR amplification of bisulfite-converted DNA followed by NGS of the amplified target loci (termed RainDrop BS-Seq) enables the sensitive, specific and simultaneous amplification of up to 4000 bisulfite-converted target loci . We and others have shown the utility of this approach to validate epigenetic alterations in a range of tissues [27–29]. In this proof of concept study, we describe the development of the UroMark assay, which uses high-throughput targeted bisulphite sequencing of urinary sediment cell DNA to provide a read-out of presence or absence of bladder cancer. This assay, in which a large comprehensive panel shows diagnostic precision, achieving high sensitivity and specificity for disease detection, would represent a potential paradigm shift in the diagnosis and surveillance of bladder cancer.
Patient characteristics of primary tissues used in (cohort 1) discovery and (cohort 2) validation
Cohort 1 (N = 116)
Cohort 2 (N = 199)
For target validation, cohort 2 (n = 199) was an independent dataset obtained from The Cancer Genome Atlas (TCGA) (https://tcga-data.nci.nih.gov/tcga/dataAccessMatrix.htm?mode=ApplyFilter&showMatrix=true&diseaseType=BLCA&tumorNormal=TN&tumorNormal=T&tumorNormal=NT&platformType=2&platformType=42). This contains data for 144 muscle-invasive bladder cancers and 20 normal urothelium samples. We supplemented the TCGA data with a further 35 methylomes generated from non-muscle-invasive disease, representative of low-grade disease from the Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT) (Madrid) (cohort 2, Table 1).
Patient characteristics of urine samples used in the assessment of the UroMark assay
Cohort 3 (N = 86)
Cohort 4 (N = 188)
The studies were conducted under the following ethics approvals: For primary tissue obtained at UCLH (10/H1306/42, 15/YH0311), BCPP (06/MRE04/65) and Madrid (CEIC 10/50); and urinary validation: 06/Q0104/57, 10/H1306/42 and 15/YH0311.
Genome-wide methylation profiling
Five hundred nanograms of DNA was bisulfite-converted and hybridised to the Infinium 450K Human Methylation array (Illumina) and processed in accordance with the manufacturer’s recommendations. DNA bisulfite conversion was carried out using the EZ DNA Methylation Kit (Zymo Research) as per manufacturer’s instructions. The R statistical software (version 3.1.2 ) was used for the subsequent data analysis. The ChAMP analysis pipeline was used to extract and analyse data from iDat files . Samples were normalised using BMIQ [33, 34]. Raw β-values (methylation value) were subjected to a stringent quality control analysis as follows: samples showing reduced coverage and probes containing SNPs were removed, and only probes with detection levels above background across all samples were retained (detection P < 0.01).
Panel identification, selection and classification
The UroMark panel was defined using pre-set criteria in the training cohort (Table 1) as follows: in order for probes to be considered as potential biomarker candidates, they had to show no or very low methylation (β <10%) in normal urothelium, blood and non-cancer urine samples and methylation (β) of >50% in bladder cancer. A <10% cut was selected from the analysis of beta values from fully unmethylated control DNA. The probes passing this filter were subsequently used to generate a classifier using a random forest model.
The random forest classification model was selected as it has been shown to be effective with a limited number of predictors (i.e. number of loci being compared) in comparison with the number of training points (i.e. the number of samples within the training cohort) . The random forest model was implemented through the Bioconductor package CAReT (version 6.0-24) . To define a robust classifier, the following steps were implemented: (A) A random selection of 80% of cases and controls was selected as training cohort, with the remaining 20% retained to form a test set; (B) a random forest model was developed for each training cohort; (C) a predicted class (cancer or normal) was generated for each corresponding test cohort; (D) for each iteration, the model and area under the curve (AUC) value were noted; (E) steps A–D were repeated 100 times, and the optimal model was determined by comparison of each AUC value. The optimal model was fixed and applied to all subsequent analysis.
Comparative testing of UroMark assay and classifiers based on the best performing 3, 5 and 10 loci panels was performed by a simple logistic regression model to calculate probabilities of each combination. The performance of individual loci is shown in Additional file 1: Table S4. The best performing combination of 3, 5 and 10 probes were selected based on a false positive rate of <10% in the training cohort. Area under the ROC (receiver operating characteristic) curve for the UroMark assay and best performing 3, 5 and 10 marker panels were calculated using the Bioconductor package pROC (version 1.7). Loci were involved in OTX1, ONECUT2, ZNF154, TBX2 and ZIC4.
RainDance microdroplet PCR of urinary DNA
RainDrop BS-seq was performed as previously described [27, 28]. Primers were designed for targeted regions in Additional file 2: Table S1. For microdroplet PCR, 7.20 μL of bisulfite-treated urinary DNA was added to 4.70 μL of 10× High-Fidelity Buffer (Invitrogen), 1.80 μL of 50 mM MgSO4 (Invitrogen), 1.62 μL of 10 mM dNTP solution mix (NEB), 3.60 μL of 4 mol L−1 betaine solution (Sigma-Aldrich), 3.60 μL of droplet stabiliser (RainDance Technologies), 1.80 μL of 100% dimethyl sulfoxide (Sigma-Aldrich) and 0.72 μL of 5 U/μL Platinum Taq Polymerase High-Fidelity (Invitrogen), to a total volume of 25 μL. The sample plate was sealed using an ALPS 50V microplate heat sealer (Thermo Scientific). The bisulfite-treated genomic DNA template mix was then applied to a fully automated ThunderStorm system (RainDance Technologies) following the manufacturer’s instructions. In brief, primer panel droplets (MethylSeq Solution, RainDance Technologies) were dispensed to a microfluidic chip. The DNA template mix was converted into droplets within the microfluidic chip. The primer pair droplets and template droplets were then paired together in a 1:1 ratio. The paired droplets passed through an electric field inducing the discrete droplets to coalesce into a single PCR droplet (26 pL); approximately one million PCR droplets are collected per sample.
The PCR droplets were processed in a PTC-225 thermocycler (MJ Research) as follows: 94 °C for 2 min; 55 cycles of 94 °C for 30 s, 54 °C for 45 s and 68 °C for 80 s; followed by 68 °C for 10 min; 4 °C until further processing. The ramp rate was set to 1 °C per second. Following PCR amplification, 70 μL of droplet destabilizer (RainDance Technologies) were added to each sample to break the PCR droplet emulsion and release the amplicons contained within the droplets. The solution was mixed well and incubated for 15 min at room temperature. The samples were purified using Agencourt AMPure XP magnetic beads (Beckman Coulter) following the manufacturer’s protocol. For each sample, 234 μL of beads were used. The samples were eluted from magnetic beads in 40 μL of EB Buffer. The integrity and concentration (fragment range 120–300 bp) of purified amplicon DNA were assessed using a High Sensitivity DNA Kit (Agilent Technologies) on a 2100 Bioanalyzer (Agilent Technologies).
To prepare the samples for high-throughput DNA sequencing, Illumina adapter sequences and unique barcodes were introduced through an additional PCR step. Fifteen nanograms of purified amplified DNA were added to 3.25 μL of 10× High-Fidelity Buffer, 0.88 μL of 50 mM MgSO4, 0.88 μL of 10 mM dNTP solution mix, 2.50 μL of 4 mol L−1 betaine solution, 1.25 μL of 100% dimethyl sulfoxide, 0.50 μl of 5 U/μL Platinum Taq Polymerase High-Fidelity and 2.5 μL of 5 μM PCR primers, to a total volume of 25 μL. All primer sequences are provided in Additional file 3: Table S2.
The samples were amplified as follows: 94 °C for 2 min; 10 cycles of 94 °C for 30 s, 56 °C for 45 s and 68 °C for 60 s; followed by 68 °C for 10 min; 4 °C until further processing. DNA was purified using a MinElute PCR Purification Kit (QIAGEN) according to the manufacturer’s protocol. Purified DNA was eluted in 10 μL of EB Buffer. The samples were quantified (fragment range 100–400 bp and 220–500 bp) using a DNA 1000 Kit on a 2100 Bioanalyzer. Fifty nanograms of each sample were subsequently pooled. The resulting sequencing library was quantified using a Qubit dsDNA BR Assay Kit.
High-throughput DNA sequencing
The pooled sequencing library (12 pmol L−1) and custom sequencing primers (0.5 μmol L−1) were applied to a MiSeq 300-cycle PE consumable cartridge (Illumina) according to the manufacturer’s protocol. The DNA sequences of the custom sequencing primers are provided in Additional file 4: Table S3. The sequencing was performed on a MiSeq DNA sequencer (Illumina) using 150 bp paired end reads.
Data and statistical analyses
Sequencing adapters were trimmed from the raw sequencing reads using the fastq-mcf tool of ea-utils v1.1.2-537. Trimmed sequencing data were mapped to an in silico bisulfite-converted human reference genome (GRCh37) using Bismark v0.7.12 . Methylation information was extracted using the methylation extractor tool of Bismark v0.7.12 . Targeted DNA sequencing analyses were performed using the R package TEQC v3.2.0 .
Genome-wide DNA methylation profiling to define bladder cancer-specific loci
To derive a bladder cancer-specific DNA methylation signature, which would allow classification of independent samples, we used a random forest model which resulted in a signature consisting of 150 CpG loci (Fig. 1b) which on the test set, resulted in a cross-validated sensitivity of 100% and specificity of 100% for the detection of cancer (Additional file 5: Figure S1A, B).
Validation of the detection panel
Developing the UroMark assay as a high-throughput NGS assay
The UroMark assay is designed as a high-throughput microdroplet-based PCR amplification system using RainDrop BS-seq [27, 28]. The encapsulation of distinct PCR reactions in microdroplets combined with NGS allows the targeted bisulfite sequencing of a large number of unique regions in parallel from limited substrate [27, 28]. We have previously validated this technology and have shown it to be highly correlated with the Infinium Human 450K DNA methylation array and also shown its utility with low template input [27, 28].
We designed a bisulfite-converted primer library to determine the methylation state of the 150 selected genomic loci. Primers were designed to interrogate both Watson and Crick strands independently where possible. Bisulfite-treated urinary DNA was subsequently used as a template for the microdroplet-based PCR amplification reaction with a RainDance ThunderStorm system.
To validate the UroMark assay using RainDrop BS-seq, we tested a second independent cohort of 188 cases (cohort 4). DNA from urinary sediment cells was obtained from 55 patients with bladder cancer, and 133 patients confirmed to be cancer-free on cystoscopy and upper track imaging (Table 2, cohort 4). All samples analysed had >30 ng of DNA; NGS data was analysed as described [27, 28]. The fraction of aligned sequencing reads mapping to the target amplicons ranged between 94.5–98.7% across the sample cohort. All amplicons of the panel amplified across the sample cohort. The average sequence coverage across all samples was 1254-fold (range 123–2673).
Combining all urine samples processed using the UroMark (n = 274), including non-cancer (n = 167) and bladder cancer (n = 107) (cohort 3 and cohort 4), the UroMark assay predicted the presence of bladder cancer with an AUC of 97% and a NPV of 97% (Fig. 4b).
Comparison with small panels
To understand how the performance of the UroMark assay was compared to the best performing single markers and small-marker panels, combinations of the top performing three, five and ten biomarkers were identified from the training cohort. Any individual marker was positive based on a simple methylation cut-off generated for each locus referenced on the methylation value of normal urothelium. The best performing markers were combined into panels of three, five or ten loci along with the loci involved in previously published in regression-based predictive classifiers developed to explore the potential for these ‘oligo’ panels [19, 21–23].
Test performance characteristics of small panels of markers (3, 5 and 10) and single markers to detect bladder cancer compared to the UroMark assay
Top individual marker
This study reports the development and proof of principle testing of high-throughput target bisulphite sequencing assay, UroMark, to interrogate cancer epigenetic alterations in urinary sediment. A non-invasive test for the detection of bladder cancer has the potential to revolutionise the diagnostic pathway, for both haematuria investigations and improved surveillance strategies for patients with established disease. In proof of concept testing of a cohort of cancer and control urine samples, we show that the high sensitivity and specificity obtained with the UroMark assay has performance characteristics which are similar to cystoscopy.
DNA methylation patterns are highly cancer cell-specific, and the ontogenic stability of these epigenetic events makes DNA methylation an ideal biomarker for the detection and diagnosis of disease. Changes in global DNA methylation patterns are a common feature of neoplastic transformation and a frequent event in bladder cancer [16–18, 38]. Previous studies have shown that changes in DNA methylation status of bladder cancer, both non-muscle-invasive and muscle-invasive bladder cancers as well as normal urothelium are reflected in the methylation status of urinary sediment cells and as such, could be a useful diagnostic marker [16, 17, 38]. Although methylation-based detection assays (alone or in combination with somatic mutation) show promise with sensitivity to detect bladder cancer between 65–98%, they have not as yet progressed into clinical practise [8, 17, 19, 21–24]. The number of loci that can be included within a panel using traditional technology to detect methylation and the low volume of substrate DNA that can be extracted from urinary sediment for large numbers of candidates have until recently been a limiting factor for assay development. Furthermore, the reliance on uniform methylation (low inter-tumour variability) as well as the effect of intra-tumour heterogeneity indicate that the performance of ‘oligo panel’ assays will be limited across a wide spectrum of stage and grade . Novel technologies, combining next-generation bisulfite sequencing with large-scale multiplex PCR, overcome these issues allowing the interrogation of a large panel of epigenetic biomarkers from a single sample, and as such, we believe this represents a paradigm shift in development strategy [27, 28].
In order to annotate the epigenetic alterations involved in bladder cancer and to define a biomarker panel, we have carried out the largest unbiased genome-wide DNA methylation screens of bladder cancer to date. Although our initial discovery panel was predominately high-grade disease, it is of note that the majority or >98% of the alterations present in high-grade bladder cancer were confirmed in low-grade disease. Using stringent criteria, we have defined a detection panel of 150 CpG loci which is relatively large but was necessary to detect all of the 260 bladder cancers included in the assay development cohort based on individual methylation expression and representing the spectrum of stage and grade. We believe the large panel and customised random forest analysis pipeline to retain specificity and sensitivity can overcome shortcomings of traditional biomarker panels which are constrained by technology. The use of NGS, particularly in cancer diagnosis, is becoming routine, and consistent with this, we demonstrate the potential for a large-scale highly multiplexed next-generation assay, the performance characteristics being highly sensitive and specific, can be achieved with a NPV comparable to cystoscopy.
The current assay was developed to answer the specific question around the primary diagnosis of bladder cancer in patients with haematuria and has yet to be tested in the recurrence setting. Urinary biomarker assays for the detection of recurrent bladder cancer have generally fared less well than in the primary diagnosis setting with sensitivities ranging from 56–80%. This low sensitivity for the detection of recurrent disease is likely due to low urinary concentration of tumour cells harbouring the small number specific alterations analysed or the analytical sensitivity of the assays used. The large panel of loci utilised in the UroMark panel in combination with the analytical sensitivity potential achieved with next generation sequencing may also allow the UroMark assay to compare favourably with cystoscopy in the recurrent setting.
The next stage will be the robust testing (Phase III biomarker road map ) of this assay in two MRC-funded trials (NCT02676180 and NCT0278428) which are currently recruiting across multiple sites in the UK. In the current development studies, we used a mixed cohort of tumours and non-cancer controls from various sources. The objective of the Phase III studies is to determine the NPV of UroMark for the detection of bladder cancer in a population of patients referred for investigation of haematuria.
In this proof of concept study, we show the potential utility of a highly multiplex bisulphite sequencing assay for the detection of bladder cancer from urinary sediment. The use of a non-invasive assay which rules out the presence of cancer with a high degree of certainty has the potential to revolutionise the treatment of bladder cancer.
The authors would like to acknowledge the UCL genomics facility at the Institute of Child Health and UCL Cancer Institute.
AF and JDK are supported by the UCL/UCLH Comprehensive Biomedical Research Centre, MRC (MR/M025411/1), Orchid and the Rosetrees Trust. DP and SB are funded by the EU-FP7 Project BLUEPRINT (282510) and the Wellcome Trust (99148). JMP is funded by the following: MINECO grant SAF2012-34378 and SAF2015-66015-R, FEDER cofounded CB/16/00228; Comunidad Autónoma de Madrid grant S2010/BMD-2470 (Oncocycle Program); and AES grant ISCIII-RETIC RD12/0036/0009. MMF is supported by an EMBO fellowship (EMBO ASTF 81-2014/Award). WST and SR are supported by the Urology Foundation.
Availability of data and materials
Upon acceptance of this manuscript, the dataset supporting the conclusions of this article will be made available through EGA.
AFe designed the study, interpreted the data and prepared the manuscript. PD, LD, PdW, WTS, AH, SR, PG, SR and AM prepared, processed and analysed the samples. MMF, FV, FR, KKC, MPZ, RTB,NDJ and LMP provided the samples. CJ and AFr provided pathological review. DSP and SB provided technical expertise and assistance with the study design. JDK designed the study, interpreted the data and prepared the manuscript. All authors have read and approved the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The studies were conducted under the following ethics approvals: For primary tissue obtained at UCLH (10/H1306/42, 15/YH0311), BCPP (06/MRE04/65) and Madrid (CEIC 10/50); and urinary validation: 06/Q0104/57, 10/H1306/42 and 15/YH0311.
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.
- Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–E86.View ArticlePubMedGoogle Scholar
- Khadra MH, Pickard RS, Charlton M, Powell PH, Neal DE. A prospective analysis of 1,930 patients with hematuria to evaluate current diagnostic practice. J Urol. 2000;163:524–7.View ArticlePubMedGoogle Scholar
- Jocham D, Stepp H, Waidelich R. Photodynamic diagnosis in urology: state-of-the-art. Eur Urol. 2008;53:1138–48.View ArticlePubMedGoogle Scholar
- Denzinger S, Burger M, Walter B, Knuechel R, Roessler W, Wieland WF, et al. Clinically relevant reduction in risk of recurrence of superficial bladder cancer using 5-aminolevulinic acid-induced fluorescence diagnosis: 8-year results of prospective randomized study. Urology. 2007;69:675–9.View ArticlePubMedGoogle Scholar
- Zaak D, Kriegmair M, Stepp H, Stepp H, Baumgartner R, Oberneder R, et al. Endoscopic detection of transitional cell carcinoma with 5-aminolevulinic acid: results of 1012 fluorescence endoscopies. Urology. 2001;57:690–4.View ArticlePubMedGoogle Scholar
- Schlake A, Crispen PL, Cap AP, Atkinson T, Davenport D, Preston DM. NMP-22, urinary cytology, and cystoscopy: a 1 year comparison study. Can J Urol. 2012;19:6345–50.PubMedGoogle Scholar
- Burke DM, Shackley DC, O'Reilly PH. The community-based morbidity of flexible cystoscopy. BJU Int. 2002;89:347–9.View ArticlePubMedGoogle Scholar
- Yossepowitch O, Herr HW, Donat SM. Use of urinary biomarkers for bladder cancer surveillance: patient perspectives. J Urol. 2007;177:1277–82.View ArticlePubMedGoogle Scholar
- Burns MB, Lackey L, Carpenter MA, Rathore A, Land AM, Leonard B, et al. APOBEC3B is an enzymatic source of mutation in breast cancer. Nature. 2013;494:366–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelly JD, Fawcett DP, Goldberg LC. Assessment and management of non-visible haematuria in primary care. BMJ. 2009;338:a3021.View ArticlePubMedGoogle Scholar
- Tilki D, Burger M, Dalbagni G, Grossman HB, Hakenberg OW, Palou J, et al. Urine markers for detection and surveillance of non-muscle-invasive bladder cancer. Eur Urol. 2011;60:484–92.View ArticlePubMedGoogle Scholar
- Lotan Y, Roehrborn CG. Sensitivity and specificity of commonly available bladder tumor markers versus cytology: results of a comprehensive literature review and meta-analyses. Urology. 2003;61:109–18.View ArticlePubMedGoogle Scholar
- van Rhijn BW, van der Poel HG, van der Kwast TH. Urine markers for bladder cancer surveillance: a systematic review. Eur Urol. 2005;47:736–48.View ArticlePubMedGoogle Scholar
- Kulis M, Esteller M. DNA methylation and cancer. Adv Genet. 2010;70:27–56.PubMedGoogle Scholar
- Kandimalla R, Van Tilborg AA, Zwarthoff EC. DNA methylation-based biomarkers in bladder cancer. Nat Rev Urol. 2013;10:327–35.View ArticlePubMedGoogle Scholar
- Beukers W, Hercegovac A, Vermeij M, Kandimalla R, Blok AC, van der Aa MM, et al. Hypermethylation of the polycomb group target gene PCDH7 in bladder tumors from patients of all ages. J Urol. 2013;190:311–6.View ArticlePubMedGoogle Scholar
- Kandimalla R, Masius R, Beukers W, Bangma CH, Orntoft TF, Dyrskjot L, et al. A 3-plex methylation assay combined with the FGFR3 mutation assay sensitively detects recurrent bladder cancer in voided urine. Clin Cancer Res. 2013;19:4760–9.View ArticlePubMedGoogle Scholar
- Su SF, de Castro Abreu AL, Chihara Y, Tsai Y, Andreu-Vieyra C, Daneshmand S, et al. A panel of three markers hyper- and hypomethylated in urine sediments accurately predicts bladder cancer recurrence. Clin Cancer Res. 2014;20:1978–89.View ArticlePubMedGoogle Scholar
- van Kessel KE, Van Neste L, Lurkin I, Zwarthoff EC, Van Criekinge W. Evaluation of an epigenetic profile for the detection of bladder cancer in patients with hematuria. J Urol. 2016;195:601–7.View ArticlePubMedGoogle Scholar
- Frantzi M, Van Kessel KE, Zwarthoff EC, Marquez M, Rava M, Malats N, et al. Development and validation of urine-based peptide biomarker panels for detecting bladder cancer in a multi-center study. Clin Cancer Res. 2016;22:OF1–OF10.View ArticleGoogle Scholar
- Costa VL, Henrique R, Danielsen SA, Duarte-Pereira S, Eknaes M, Skotheim RI, et al. Three epigenetic biomarkers, GDF15, TMEFF2, and VIM, accurately predict bladder cancer from DNA-based analyses of urine samples. Clin Cancer Res. 2010;16:5842–51.View ArticlePubMedGoogle Scholar
- Fantony JJ, Abern MR, Gopalakrishna A, Owusu R, Jack Tay K, Lance RS, et al. Multi-institutional external validation of urinary TWIST1 and NID2 methylation as a diagnostic test for bladder cancer. Urol Oncol. 2015;33:387 e1-6.View ArticlePubMedGoogle Scholar
- Renard I, Joniau S, van Cleynenbreugel B, Collette C, Naome C, Vlassenbroeck I, et al. Identification and validation of the methylated TWIST1 and NID2 genes through real-time methylation-specific polymerase chain reaction assays for the noninvasive detection of primary bladder cancer in urine samples. Eur Urol. 2010;58:96–104.View ArticlePubMedGoogle Scholar
- Ward DG, Baxter L, Gordon NS, Ott S, Savage RS, Beggs AD, et al. Multiplex PCR and next generation sequencing for the non-invasive detection of bladder cancer. PLoS One. 2016;11:e0149756.View ArticlePubMedPubMed CentralGoogle Scholar
- Dahmcke CM, Steven KE, Larsen LK, Poulsen AL, Abdul-Al A, Dahl C, et al. A prospective blinded evaluation of urine-DNA testing for detection of urothelial bladder carcinoma in patients with gross hematuria. Eur Urol. 2016;70(6):916–19.View ArticlePubMedGoogle Scholar
- Tewhey R, Warner JB, Nakano M, Libby B, Medkova M, David PH, et al. Microdroplet-based PCR enrichment for large-scale targeted sequencing. Nat Biotechnol. 2009;27:1025–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Paul DS, Guilhamon P, Karpathakis A, Butcher LM, Thirlwell C, Feber A, et al. Assessment of RainDrop BS-seq as a method for large-scale, targeted bisulfite sequencing. Epigenetics. 2014;9(5):678–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Guilhamon P, Eskandarpour M, Halai D, Wilson GA, Feber A, Teschendorff AE, et al. Meta-analysis of IDH-mutant cancers identifies EBF1 as an interaction partner for TET2. Nat Commun. 2013;4:2166.View ArticlePubMedPubMed CentralGoogle Scholar
- Komori HK, LaMere SA, Torkamani A, Hart GT, Kotsopoulos S, Warner J, et al. Application of microdroplet PCR for large-scale targeted bisulfite sequencing. Genome Res. 2011;21:1738–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowe R, Rakyan VK. Marmal-aid—a database for Infinium HumanMethylation450. BMC Bioinformatics. 2013;14:359.View ArticlePubMedPubMed CentralGoogle Scholar
- Tan WS, Rodney S, Lamb B, Feneley M, Kelly J. Management of non-muscle invasive bladder cancer: a comprehensive analysis of guidelines from the United States. Eur Asia Cancer Treat Rev. 2016;47:22–31.View ArticleGoogle Scholar
- Kiss B, Schneider S, Thalmann GN, Roth B. Is thermochemotherapy with the Synergo system a viable treatment option in patients with recurrent non-muscle-invasive bladder cancer? Int J Urol. 2015;22:158–62.View ArticlePubMedGoogle Scholar
- Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics. 2014;30:428–30.View ArticlePubMedGoogle Scholar
- Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29:189–96.View ArticlePubMedGoogle Scholar
- Kuhn M. Building predictive models in R using the caret package. J Stat Software. 2008;28(5):1–26.View ArticleGoogle Scholar
- Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27:1571–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Hummel M, Bonnin S, Lowy E, Roma G. TEQC: an R package for quality control in target capture experiments. Bioinformatics. 2011;27:1316–7.View ArticlePubMedGoogle Scholar
- Kandimalla R, van Tilborg AA, Zwarthoff EC. DNA methylation-based biomarkers in bladder cancer. Nat Rev Urol. 2013;10(6):327–35.View ArticlePubMedGoogle Scholar
- Gerlinger M, Catto JW, Orntoft TF, Real FX, Zwarthoff EC, Swanton C. Intratumour heterogeneity in urologic cancers: from molecular evidence to clinical implications. Eur Urol. 2015;67:729–37.View ArticlePubMedGoogle Scholar
- CRUK. CRUK Biomarker Roadmap. http://www.cancerresearchuk.org/prod_consump/groups/cr_common/@fre/@fun/documents/generalcontent/cr_027486.pdf.