- Open Access
A multiplex urinary immunoassay for bladder cancer detection: analysis of a Japanese cohort
Journal of Translational Medicine volume 14, Article number: 287 (2016)
Bladder cancer (BCa) is among the most commonly diagnosed malignancies worldwide, and due the high rate of post-operative disease recurrence, it is one of the most prevalent in many countries. The development of non-invasive molecular assays that can accurately detect and monitor BCa would be a major advance, benefiting both patients and healthcare systems. We have previously identified a urinary protein biomarker panel that is being developed for application in at-risk patient cohorts. Here, we investigated the potential utility of the multiplex assay in a Japanese cohort.
The Japanese study cohort collected from urology clinics at two institutions was comprised of a total of 288 subjects. The protein biomarker panel (IL8, MMP9, MMP10, ANG, APOE, SDC1, A1AT, PAI1, CA9, VEGFA) was monitored in voided urine samples collected prior to cystoscopy using a custom multiplex ELISA assay. The diagnostic performance of the biomarker panel was assessed using receiver operator curves, predictive modeling and descriptive statistics.
Urinary biomarker concentrations were significantly elevated in cases versus controls, and in cases with high-grade and muscle-invasive tumors. The AUC for the 10-biomarker assay was 0.892 (95 % confidence interval 0.850–0.934), with an overall diagnostic sensitivity specificity of 0.85 and 0.81, respectively. A predictive model trained on the larger institutional cohort correctly identified 99 % of the cases from the second institution.
Urinary levels of a 10-biomarker panel enabled discrimination of patients with BCa. The multiplex urinary diagnostic assay has the potential to be developed for the non-invasive detection of BCa in at-risk Japanese patients.
Bladder cancer (BCa) is a major burden for both patients and health care systems worldwide . At initial presentation, the majority of cases are non-muscle invasive tumors, which are not immediately life-threatening. However, more than 70 % of patients with operable BCa will have a recurrence during the first 2 years after diagnosis and the subsequent lesions can progress to being muscle invasive and metastatic . Thus, once treated, BCa patients are under continual surveillance with routine cystoscopy examinations for early detection of new BCa development. Due to the prolonged natural history of BCa, plus the prolonged and invasive nature of follow-up and treatment strategies, it is one of the most expensive malignancies to manage on a per-patient basis . Although Japan and other Asian countries have a somewhat lower incidence of BCa than Europe and the Americas, the clinico-pathological characteristics of the disease are similar. In Japan, around 12,000 patients are newly diagnosed and 5000 patients die from the disease annually. While the overall incidence trend has been relatively stable over the past decade [4, 5], because the onset is most often seen in elderly patients, it is becoming a major social issue in an aging Japanese society.
In a series of previous studies, we have identified panels of protein biomarkers that are significantly associated with BCa [6, 7]. Through several studies [8–10] we have refined the target composition, and the immunoassay approach , in order to develop a viable non-invasive assay for BCa detection. In this study, we tested the potential utility of the multiplex immunoassay for the detection of BCa in a Japanese cohort. Evaluation of 288 subjects obtained from two independent institutes confirmed a significant association of the tested biomarkers with the presence of BCa. The multiplex assay achieved a strong overall diagnostic performance achieving 85 % sensitivity and 81 % specificity (AUC 0.892). This retrospective phase II study confirms the potential of using urinary protein biomarker signatures for the non-invasive detection of BCa, and suggests that the described multiplex immunoassay could aid in Japanese urology patient management.
Patients and specimen processing
Under Western Institutional Review Board approval (IRB #Rosser 2014-1), previously collected and banked voided urine samples shipped to University of Hawaii Cancer Center for analysis. Originally, voided urine samples were collected prior to cystoscopy, and stored frozen until analysis. The study cohort (Table 1) was comprised of subjects from two independent institutions; the Kyoto cohort consisted of 80 subjects with newly diagnosed BCa, and 67 subjects with no active BCa (controls) including subjects with voiding symptoms, renal or prostate cancer, urinary tract infections or urolithiasis. The Nara cohort consisted of 131 subjects with newly diagnosed BCa. For the bladder cancer case group, histological confirmation of urothelial carcinoma, including grade and stage was defined from excised tissue. Frozen aliquots of urine samples were thawed and protein content was measured using a Pierce 660-nm Protein Assay Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA). The concentration of urinary creatinine was measured using a commercially available enzymatic assay (Cat#KGE005 R&D Systems Inc., Minneapolis, MN, USA) according to the manufacturer’s instructions. Normalization to creatinine is the standard means of expressing analytes in urine to minimize the impact of patient hydration status .
The concentrations of the 10 proteins (IL8, MMP9, MMP10, ANG, APOE, SDC1, A1AT, PAI1, CA9 and VEGFA) were monitored using a custom multiplex immunoassay using MULTI-ARRAY® technology (Meso Scale Diagnostics, LLC, Rockville, MD). The multiplex assay, using MULTI-SPOT® plates, is based on a proprietary combination of electrochemiluminescence detection and patterned arrays [13–15]. Final monoclonal antibody pairs (capture and detection) were selected based on sensitivity, specificity, physical properties, and recognition of native protein, as described previously . In order to ensure detection across the range of protein concentrations, 7 of the 10 assays (IL8, MMP9, MMP10, APOE, PAI1, CA9 and VEGFA) were multiplexed in one well, and the other 3 assays (ANG, SDC1, A1AT) in a separate well. Urine samples were diluted fourfold for the 7-plex assay and 200-fold for the 3-plex assay. A seven point standard curve across the 4 log dynamic range of the assays was included in the current assay design. Urine samples were handled on ice and diluted with MSD Assay Diluent 37 . Samples and standards (50 μl) were added to the plate and incubated for 2 h at RT. After washing, SULFO-TAG conjugated detection antibody mix (25 μl) was added and incubated for 2 h at RT. After washing, the ensuing immunoassay complex was incubated with MSD Read Buffer (150 μl) and electrochemiluminescence measured on the QuickPlex® SQ 120 (MSD) instrument. Standard curves were constructed using MSD Discovery Workbench® 4.0, which allows for the selection of multiple non-linear and linear equations to fit the standard curve. Optimal curve fits were determined by visual graph evaluation and comparison of akaike’s information criteria (AIC) values .
Wilcoxon rank sum tests were used to determine the association between each biomarker and BCa. Nonparametric receiver operating characteristic (ROC) curves were generated to plot assay sensitivity against the false-positive rate (1-specificity). The relative ability of each biomarker to indicate BCa was evaluated by calculating the area under the curve (AUC), and AUCs were compared by Chi-square test. The sensitivity and specificity of each biomarker and their combinations were estimated at the optimal cutoff value defined by the Youden index . Furthermore, the formula for calculating the combined score of prediction was −3.3600 + (0.4889)*normalized IL8 + (−2.9371)*normalized MMP9 + (0.2112)*normalized A1AT + (−0.0470)*normalized ANG + (−0.5617)*normalized VEGFA + (1.7788)*normalized CA9 + (−0.1412)*normalized MMP10 + (−3.0268)*normalized APOE + (47.5054)*normalized PAI1 + (7.4465)*normalized SDC1 > 0.80385). To assess the independent association between biomarkers and BCa, we used logistic regression analysis with BCa status (yes vs. no) as the response variable and biomarker concentrations as explanatory variables. Multiple-level logistic regression was used to evaluate biomarker association with tumor grade and muscle-invasive disease. The all-subset method was used to evaluate the predictive value of each possible combination of biomarkers, and the Bayesian information criterion (BIC) was used to compare models. The BIC, a widely used criterion in model selection, balances the model likelihood and the number of biomarkers included in the model . The Bootstrap method (using 1000 Bootstrap samples) was used  to select the most efficient and stable predictive model. Sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated from ROC data. Statistical significance in this study was set at P < 0.05 and all reported P values were 2-sided. All analyses were performed using SAS software version 9.4.
Clinical, pathologic and demographic characteristics of the 288 subjects (183 BCa, 96 benign controls, and 41 healthy volunteers) comprising the study cohort are listed in Table 1. Of the 211 bladder cancer cases, 170 were classified as non-muscle invasive bladder cancer (NMIBC; stages Ta, Tis, T1), and 41 were muscle invasive bladder cancer (MIBC; stage ≥ T2). Eighty-seven cases were reported as low-grade carcinoma and 98 cases as high-grade. Gender distribution (BCa incidence is higher in men), and slightly older median age for patients in the BCa group reflect typical BCa incidence statistics .
Urinary biomarker levels and assay performance
Biomarker concentration data were normalized to urinary creatinine  to adjust for urine volume variations and to be consistent with our previous studies [8–11]. Urinary concentrations of all 10 biomarkers were significantly elevated in patients with bladder cancer compared with controls (Table 2). The concentrations of 9 of the 10 biomarkers (except SDC1) were also significantly elevated in patients with high-grade BCa relative to low-grade BCa, and MIBC relative to NMIBC (Table 2). The diagnostic performance of each individual biomarker was analyzed using nonparametric ROC analyses . As observed previously, many of the individual biomarkers achieved respectable diagnostic performance values, with SDC1 and SERPINA1 both recording AUC values >0.80. Table 3 provides AUC and corresponding sensitivity, specificity, PPV, and NPV values for all biomarkers tested. The optimal test was derived by combination of the 10 biomarkers into a diagnostic panel. In line with our previous multi-institutional studies with cohorts from the USA and Europe [9, 10], the 10-biomarker assay achieved an AUC of 0.8925 in the Japanese cohort (Fig. 1), with an overall diagnostic sensitivity of 0.895 and specificity of 0.806 (Table 3). As expected, sensitivity and specificity values of the combined assay (Table 4) were higher for high-grade (90.3/85.1 %) and for MIBC bladder cancer (85.4/95.5 %), but the overall performance was comparable. Although data was incomplete, voided urine cytology data indicated sensitivity values of ~45 % for those subjects tested in this cohort.
The overall classification data was obtained through analysis of data from the total cohort, but given that subjects were combined from two independent institutes, we could also investigate the ability of a predictive model derived from one institute to identify cases in the second institute. The Kyoto cohort (147 subjects) had a balance of 80 cases and 67 controls, so this was used to build a predictive model. The best logistic regression model for the 10-biomarker panel correctly identified 56 of the 80 Kyoto cases. Although this was less accurate than the overall classification analyses described above, when applied to the Nara cohort the Kyoto-derived predictive model correctly identified 130 of the 131 Nara cases as harboring BCa.
Bladder cancer is among the most common worldwide malignancies, and due to the high recurrence rate, it is also one of the most burdensome for both patients and healthcare systems. Current guidelines for BCa diagnosis recommend cystoscopy coupled with voided urine cytology (VUC), but cystoscopy is an invasive procedure associated with discomfort and anxiety, infection and trauma [21, 22]. Non-invasive VUC is used routinely as an adjunct to cystoscopy, but the subjective assay can be impaired by inter-observer variability and has poor sensitivity, especially for low-grade and low-stage tumors [23, 24]. Urine-based, quantitative assays that can non-invasively detect BCa would improve the rapid diagnosis of BCa and help to avoid unnecessary invasive cystoscopy and biopsy procedures, however, while a number of such assays have been proposed, to date, such assays have lacked adequate accuracy to replace VUC, or to support or guide cystoscopy. The lack of accuracy of current urinary tests may be due to the fact that they typically rely on measuring a single biomarker. A shift from single biomarker assays to multiplex molecular signatures that reflect the molecular pathways of BCa provides an opportunity to develop assays with clinical utility across the breadth of disease states. The advent of high-throughput technologies has facilitated the discovery of more complex molecular signatures with diagnostic or prognostic potential, and a number of gene expression signature assays are now being incorporated into clinical practice [25, 26].
Through proteomic profiling of naturally micturated urine, we were able to identify a set of urinary proteins that were associated with the presence of BCa [6, 7]. One advantage of profiling the urine component directly is the ability to compare samples collected from subjects with non-malignant conditions, a situation that is not feasible using surgically excised solid tissue samples. The candidate biomarkers were refined through independent validation studies [27, 28], and a 10-biomarker panel was established for development. The diagnostic performance of the 10-protein assay was then tested in two independent cohorts, one  comprised of 127 subjects (diagnostic sensitivity of 92 %, specificity of 97 %), and one  comprised of 308 patients obtained from multiple institutions in the USA and Europe (sensitivity 74 %, specificity 90 %). Taken together, these studies illustrate the reproducibility of the diagnostic protein panel in terms of a robust biomarker panel that achieved similar performance data across multiple, independent cohorts. Based on these results, we investigated the feasibility of developing a multiplex assay that could accurately and simultaneously monitor the diagnostic biomarkers in an efficient format for potential clinical application. A custom-designed multiplex assay using MULTI-ARRAY® technology (meso scale diagnostics, LLC) was constructed, and the analytical performance was compared with data obtained from individual ELISA assays directed at each of the same 10 urinary proteins. The multiplex assay achieved excellent concordance, and improved detection range and technical sensitivity . The diagnostic performance of the assay was confirmed in an independent cohort of 200 subjects (sensitivity 85 %, specificity 81 %). The integration of multiplex molecular signatures into a single assay is beneficial for clinical translation through reduced sample volume, decreased processing time, low cost analysis and low reagent consumption, and the MULTI-ARRAY® assay can be readily implemented into a CLIA certified laboratory setting using existing instrumentation and skillsets.
Here, we extended the evaluation of the multiplex BCa diagnostic assay to a phase II, multicenter study with a cohort comprised of 288 Japanese subjects. For overall classification, the multiplex assay achieved a sensitivity of 84.8 % and specificity of 80.6 % (AUC of 0.892). Biomarker expression and assay performance increased with increasing tumor stage and grade. The results in the Japanese cohort were very much in line with our previous studies evaluating cohorts from the USA and Europe [9–11], but as we refine and optimize the technical aspects of the assay and analyze additional samples, we expect to be able to further improve assay performance. While subtle differences in composition of optimal urinary diagnostic protein panels may exist in cohorts comprised of different ethnic populations, given the translation of an assay optimized in diverse cohorts to a Japanese population, these data are very encouraging, and further validate the concept that a multiplex panel of urinary protein biomarkers can assist the noninvasive diagnosis of BCa.
We recognize that the study has a number of limitations. The cohort was comprised of >50 % cases, but disease prevalence is typically considerably lower in routine urologic practice. While it was important to initially test enough cases to achieve statistical significance, it will be necessary to perform additional studies that reflect urology clinic presentation to assess predictive value of the assay. It will also be necessary to include more diverse controls that may be under-represented in our study cohort. To address these issues, we have launched a multi-institute, prospective clinical trial, which will assess the multiplex diagnostic assay in subjects with gross hematuria, microscopic hematuria and history of BCa on tumor surveillance. Such a study would minimize selection bias, better represent urological disease prevalence, and evaluate potential confounding comorbidities in the study population. Furthermore, in this study, we focused only on subjects with primary BCa, however, we have previously reported that the multiplex assay is also accurate for the detection of recurrent BCa (sensitivity 79 %, specificity 88 %), outperforming the Urovysion cytogenetic assay and VUC . The inclusion of Japanese patients on routine surveillance after primary BCa treatment will be an additional goal and will enable evaluation of potential prognostic utility of the assay. Finally, we have initiated the development of BCa diagnostic nomograms that incorporate biomarker data with relevant clinical information (e.g., age, sex, race, and tobacco history) in US cohort studies, and this can also be extended to future Japanese cohort studies.
Bladder cancer is a common neoplastic disease with a high rate of recurrence and progression, and the recurrence phenomenon makes it one of the most prevalent cancers worldwide. The development of robust non-invasive, urine-based assay for the detection of BCa is clinically urgent. In this study, we have been able to successfully translate a multiplex diagnostic assay derived from diverse cohort studies to the analysis of a Japanese population. The diagnostic assay achieved encouraging performance values and will be the focus of ongoing studies to investigate the potential added value of the multiplex assay if integrated into clinical decision making.
receiver operator curves
area under the curve
plasminogen activator inhibitor-1
carbonic anhydrase 9
vascular endothelial growth factor A
meso scale diagnostics
akaike’s information criteria
bayesian information criterion
positive predictive value
negative predictive value
non-muscle invasive bladder cancer
muscle invasive bladder cancer
voided urine cytology
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Millan-Rodriguez F, Chechile-Toniolo G, Salvador-Bayarri J, Palou J, Algaba F, Vicente-Rodriguez J. Primary superficial bladder cancer risk groups according to progression, mortality and recurrence. J Urol. 2000;164(3 Pt 1):680–4.
Botteman MF, Pashos CL, Redaelli A, Laskin B, Hauser R. The health economics of bladder cancer: a comprehensive review of the published literature. Pharmacoeconomics. 2003;21(18):1315–30.
Kakehi Y, Hirao Y, Kim WJ, Ozono S, Masumori N, Miyanaga N, Nasu Y, Yokomizo A. Bladder cancer working group report. Jpn J Clin Oncol. 2010;40(Supplement 1):i57–64.
Zaitsu M, Toyokawa S, Tonooka A, Nakamura F, Takeuchi T, Homma Y, Kobayashi Y. Sex differences in bladder cancer pathology and survival: analysis of a population-based cancer registry. Cancer Med. 2015;4(3):363–70.
Kreunin P, Zhao J, Rosser CJ, Urquidi V, Lubman DM, Goodison S. Bladder cancer associated glycoprotein signatures revealed by urinary proteomic profiling. J Proteom Res. 2007;6(7):2631–9.
Yang N, Feng S, Shedden K, Xie X, Liu Y, Rosser CJ, Lubman DM, Goodison S. Urinary glycoprotein biomarker discovery for bladder cancer detection using LC-MS/MS and label-free quantification. Clin Cancer Res. 2011;17(10):3349–59.
Goodison S, Chang M, Dai Y, Urquidi V, Rosser CJ. A multi-analyte assay for the non-invasive detection of bladder cancer. PLoS ONE. 2012;7:e47469.
Rosser CJ, Ross S, Chang M, Dai Y, Mengual L, Zhang G, Kim J, Urquidi V, Alcaraz A, Goodison S. Multiplex protein signature for the detection of bladder cancer in voided urine samples. J Urol. 2013;190(6):2257–62.
Chen LM, Chang M, Dai Y, Chai KX, Dyrskjøt L, Sanchez-Carbayo M, Szarvas T, Zwarthoff EC, Lokeshwar V, Jeronimo C, Parker AS, Ross S, Borre M, Orntoft TF, Young PR, Urquidi V, Goodison S, Rosser CJ. External validation of a multiplex urinary protein panel for the detection of bladder cancer in a multicenter cohort. Cancer Epidemiol Biomark Prev. 2014;23(9):1804–12.
Shimizu Y, Furuya H, Bryant GP, Chan O, Dai Y, Thornquist MD, Goodison S, Rosser CJ. A multiplex immunoassay for the non-invasive detection of bladder cancer. J Transl Med. 2016;14(1):1.
Reid CN, Stevenson M, Abogunrin F, Ruddock MW, Emmert-Streib F, Lamont JV, et al. Standardization of diagnostic biomarker concentrations in urine: the hematuria caveat. PLoS ONE. 2012;7:e53354.
O’Shannessy DJ, Somers EB, Palmer LM, Thiel RP, Oberoi P, Heath R, Marcucci L. Serum folate receptor alpha, mesothelin and megakaryocyte potentiating factor in ovarian cancer: association to disease stage and grade and comparison to CA125 and HE4. J Ovarian Res. 2013;6(1):29.
Zhang Q, Davis KJ, Hoffmann D, Vaidya VS, Brown RP, Goering PL. Urinary biomarkers track the progression of nephropathy in hypertensive and obese rats. Biomark Med. 2014;8:85–94.
Eastman PS, Manning WC, Qureshi F, Haney D, Cavet G, Alexander C, Hesterberg LK. Characterization of a multiplex, 12-biomarker test for rheumatoid arthritis. J Pharm Biomed Anal. 2012;70:415–24.
Motulsky H, Christopoulos H. Fitting models to biological data using linear and nonlinear regression: a practical guide to curve fitting. Oxford: Oxford University Press; 2004. p. 17.
Fluss R, Faraggi D, Reiser B. Estimation of the Youden Index and its associated cutoff point. Biom J. 2005;47:458–72.
Leon LF, Tsai CL. Assessment of model adequacy for Markov regression time series models. Biometrics. 1998;54:1165–75.
Austin PC, Tu JV. Bootstrap methods for developing predictive models. Am Stat. 2004;58:131–7.
Pepe MS, Feng Z, Janes H, Bossuyt PM, Potter JD. Pivotal evaluation of the accuracy of a biomarker used for classification or prediction: standards for study design. J Natl Cancer Inst. 2008;100:1432–8.
Avritscher EB, Cooksley CD, Grossman HB, Sabichi AL, Hamblin L, Dinney CP, et al. Clinical model of lifetime cost of treating bladder cancer and associated complications. Urology. 2006;68:549–53.
Stein M, Lubetkin D, Taub HC, Skinner WK, Haberman J, Kreutzer ER. The effects of intraurethrallidocaine anesthetic and patient anxiety on pain perception during cystoscopy. J Urol. 1994;151:1518–21.
Nakamura K, Kasraeian A, Iczkowski KA, Chang M, Pendleton J, Anai S, et al. Utility of serial urinary cytology in the initial evaluation of the patient with microscopic hematuria. BMC Urol. 2009;9:12.
Kumar A, Kumar R, Gupta NP. Comparison of NMP22 Bladderchek test and urine cytology for the detection of recurrent bladder cancer. Jpn J Clin Oncol. 2006;36:172–5.
Nguyen B, Cusumano PG, Deck K, Kerlin D, Garcia AA, Barone JL, Rivera E, Yao K, de Snoo FA, van den Akker J, Stork-Sloots L, Generali D. Comparison of molecular subtyping with blueprint, mammaprint, and targetprint to local clinical subtyping in breast cancer patients. Ann Surg Oncol. 2012;19:3257–63.
Abramowicz M, Zuccotti G, Pflomm JM. A stool DNA test (Cologuard) for colorectal cancer screening. JAMA. 2014;312(23):2566.
Urquidi V, Chang M, Dai Y, Kim J, Wolfson ED, Goodison S, Rosser CJ. IL8 as a urinary biomarker for the detection of bladder cancer. BMC Urol. 2012;12:12.
Urquidi V, Kim J, Chang M, Dai Y, Rosser CJ, Goodison S. CCL18 in a multiplex urine-based assay for the detection of bladder cancer. PLoS ONE. 2012;7(5):e37797.
Rosser CJ, Chang M, Dai Y, Ross S, Mengual L, Alcaraz A, Goodison S. Urinary protein biomarker panel for the detection of recurrent bladder cancer. Cancer Epidemiol Biomarkers Prev. 2014;23(7):1340–5.
OO, YM, KY, MM, SO, KF and KT provided samples reviewed ELISA results and drafted and reviewed manuscript. YS and HF organized the project and collated all data. YD organized and analyzed data. SG participated in study design and helped to draft the manuscript. CJR conceived the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
We thank the individuals that consented to be included in the study cohort for this project.
The following authors declare that they have no competing interests. Osamu Ogawa, Yoshiyuki Matsui, Makito Miyake, Sayuri Ohnishi, Kiyohide Fujimoto, Yunfeng Dai, Yoshiko Shimizu, Hideki Furuya. The following authors declare that they have a competing interest. Steve Goodison affiliated with Nonagen Bioscience Corp. Charles J Rosser affiliated with Nonagen Bioscience Corp.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article.
Consent for publication
The manuscript does not contain any individual person’s data in any form.
Ethics approval and consent to participate
This study was approved by the Western Institutional Review Board approval (IRB #Rosser 2014-1).
This work was supported by research grants NIH/NCI R44 CA173921 (SG), Weinman Foundation Fund (CJR), 1 R01 CA198887 (CJR) and 5P30CA0717890-6071 (CJR).