Bladder Cancer (BLCa) inter-patient heterogeneity is considered the primary cause of tumor reoccurrence and treatment failure, suggesting that BLCa patients could benefit from a more personalized treatment approach. Patient-derived organoids (PDOs) have been successfully used as a functional model for predicting drug response in different cancer types. In our study, we established BLCa PDO cultures from different BLCa stages. BLCa PDOs preserve the histological and molecular heterogeneity of the parental tumors, including their multiclonal genetic landscapes. BLCa PDOs consistently share key genetic alterations detected in parental tumors, mirroring tumor evolution in longitudinal sampling. Our drug screening pipeline was implemented using BLCa PDOs, testing both standard-of-care and additional FDA-approved compounds for other solid tumors. Integrative analysis of drug response profiles with matched PDO genomic analysis was used to determine enrichment thresholds for candidate markers of therapy resistance and sensitivity. By assessing the clinical history of longitudinally sampled cases, the clonal evolution of the disease could be determined and matched with drug response profiles. In conclusion, we have developed a clinically relevant pipeline for drug response profile assessment and discovery of candidate markers of therapy resistance.