Mueller, Ryan

Abstract Introduction: Urine cell-free DNA (cfDNA) is a promising ultra-noninvasive analyte for genitourinary cancers (GU) cancers. In this study, we investigated urine cfDNA fragmentomic features and evaluated their ability to predict the tumor tissue of origin. Methods: A total of 204 GU cancer patients and 34 healthy adults were enrolled in this prospective study. Preoperative urine was collected from 89 bladder cancer (BC) patients (64% muscle-invasive) undergoing cystectomy, 65 renal cell carcinoma (RCC) patients undergoing nephrectomy and 50 metastatic prostate cancer (mPC) patients. Urine cfDNA was isolated and sequenced at 5x genome-wide coverage on a NovaSeq S4 flow cell. To identify malignant tissue of origin, we analyzed somatic copy number alterations (CNAs), genome-wide fragment length ratios, and 4-mer end motifs. CNAs and tumor fraction (TFx) were quantified using ichorCNA in 1-Mb windows, while chromosome arm-level fragment count z-scores and fragment length ratios were calculated in 5-Mb bins across the genome. Relative frequencies of all 256 possible 4-mer end motifs were calculated for each sample. Individual logistic regression models were developed based on each fragmentomic feature using leave-one-out cross validation (LOOCV). Outputs from these models were combined into multi-feature XGBoost models. Cohort-specific end-motif enrichments were identified using Kruskal-Wallis tests with FDR correction. Motifs were then hierarchically clustered based on their cohort-wise median frequencies to define groups of motifs with shared fragmentation patterns. Results: Across urine cfDNA in the GU cohort, BC showed the highest mean TFx at 10.0%, followed by mPC at 3.6% and RCC at 2.7%. A machine learning model based on 4-mer end motifs achieved a mean ROC AUC of 0.933 with classification accuracies of 63% for healthy, 87% for BC, 72% for mPC, and 83% for RCC. This model outperformed models based on ichorCNA (AUC = 0.767), chromosome arm fragment-count z-scores (AUC = 0.827), and fragment length ratios (AUC = 0.674). Combining end motifs with other fragmentomic features did not improve prediction accuracy. Multi-feature models incorporating end motifs with ichorCNA, chromosome arm-level z-scores, or fragment-length ratios yielded mean AUCs between 0.908-0.924, lower than the end-motif only model. Kruskal-Wallis tests identified 179 of 256 end motifs (79%) as significantly different across healthy, BC, mPC, and RCC. Hierarchical clustering of cohort-level median motifs revealed distinct end motif patterns: Motifs beginning with CT were enriched in RCC; AAAA was strongly enriched in BC; and motifs beginning with CA, AC, and AG were enriched in mPC. Conclusion: Urine cfDNA 4-mer end motifs capture GU cancer-specific biological patterns that enable accurate tissue-of-origin prediction, underscoring their promise as a diagnostic tool. Citation Format: Jessica Linford, Pradeep S. Chauhan, Irfan Alahi, Yohan Kim, Arpit Panda, Nathan Colon, Ryan Mueller, Faridi Qaium, Eric H. Kim, Melissa A. Reimers, Zachary L. Smith, Woodson W. Smelser, Fabrice Lucien-Matteoni, Aadel A. Chaudhuri. Urine cfDNA 4-mer end-motif signatures outperform other fragmentomic features for tissue-of-origin classification in genitourinary cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 2590.

Abstract Introduction: We previously developed a urine liquid biopsy assay to measure urine tumor DNA (utDNA) levels in bladder cancer (BC) patients. Here, we extended this approach significantly by layering in cell-free DNA fragmentomic and cell-free RNA analysis, and extending the method to the full spectrum of genitourinary (GU) cancers, to enable both GU cancer detection and tumor tissue of origin prediction. Methods: A total of 155 GU cancer patients and 34 healthy adults were enrolled in this prospective study. Preoperative urine was collected from 94 bladder cancer (BC) patients, including 65 with muscle-invasive disease, all undergoing cystectomy. Additionally, urine samples were collected from 18 renal cell carcinoma (RCC) patients undergoing nephrectomy and 42 metastatic prostate cancer (mPC) patients. Urine cfDNA and cfRNA were isolated and sequenced, with cfRNA sequenced to 100 million reads per sample on an Illumina NovaSeq X and cfDNA sequenced at 5x genome-wide coverage on a Novoseq S4 flow cell. To identify malignant tissue of origin, we analyzed somatic copy number alterations (CNAs), genome-wide fragment length ratios, and 4-mer end motifs. Tumor fraction (TFx) was estimated using ichorCNA, while DELFI was employed for fragment length ratios in 5 megabase bins across the genome and chromosome arm-level z-scores. Individual logistic regression models were developed based on each fragmentomic feature, and validated using a training and held-out testing approach. Outputs from these models were combined into a single comprehensive model, also trained and held-out tested using the same approach. Urine cfRNA transcript isoform abundance was quantified using Salmon and normalized to Transcripts Per Million (TPM) to account for sequencing depth variations. Results: Bladder cancer (BC) showed the highest mean tumor fraction (TFx, 10.5%), followed by RCC (3.8%) and mPC (3.4%). Among 94 BC patients, 42% achieved a pCR, while 58% had residual disease and did not. Patients without pCR had significantly higher urine TFx compared to those with pCR (14.7% vs. 4.7%, p=0.0008). A machine learning model based on genome-wide fragment length ratios achieved an ROC AUC of 0.74, while models using genome-wide CNAs and end motifs achieved AUCs of 0.81 and 0.88, respectively. Combining all features into a single model yielded an AUC of 0.92, with classification accuracies of 89% for bladder cancer, 80% for RCC, and 60% for prostate cancer. Additionally, urine cfRNA analysis using the top 1, 000 coding transcripts achieved an AUC of 0.95, with classification accuracies of 100% for mPC, 87% for BC, and 78% for RCC. Conclusion: Fragmentomic features from urine cfDNA and transcriptomic features from urine cfRNA demonstrate strong potential for predicting tumor tissue-of-origin in GU cancers, underscoring the promise of non-invasive diagnosis across the spectrum of GU malignancies through urine liquid biopsy. Citation Format: Pradeep Singh Chauhan, Irfan Alahi, Jessica Linford, Arpit Panda, Nathan Colon, Ryan Mueller, Faridi Qaium, Peter K. Harris, John Sheng, Eric H. Kim, Melissa A. Reimers, Zachary L. Smith, Woodson W. Smelser, Aadel A. Chaudhuri. Predicting genitourinary cancer tissue-of-origin using urine cell-free DNA fragmentomics and cell-free RNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 3258.

Abstract Purpose: Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor signaling inhibitors (ARSI) is often lethal. Liquid biopsy biomarkers for this deadly form of disease remain under investigation, and underpinning mechanisms remain ill-understood. Experimental Design: We applied targeted cell-free DNA (cfDNA) sequencing to 126 patients with mCRPC from three academic cancer centers and separately performed genome-wide cfDNA methylation sequencing on 43 plasma samples collected prior to the initiation of first-line ARSI treatment. To analyze the genome-wide sequencing data, we performed nucleosome positioning and differential methylated region analysis. We additionally analyzed single-cell and bulk RNA sequencing data from 14 and 80 patients with mCRPC, respectively, to develop and validate a stem-like signature, which we inferred from cfDNA. Results: Targeted cfDNA sequencing detected AR/enhancer alterations prior to first-line ARSIs that correlated with significantly worse progression-free survival (P = 0.01; HR = 2.12) and overall survival (P = 0.02; HR = 2.48). Plasma methylome analysis revealed that AR/enhancer lethal mCRPC patients have significantly higher promoter-level hypomethylation than AR/enhancer wild-type mCRPC patients (P < 0.0001). Moreover, gene ontology and CytoTRACE analysis of nucleosomally more accessible transcription factors in cfDNA revealed enrichment for stemness-associated transcription factors in patients with lethal mCRPC. The resulting stemness signature was then validated in a completely held-out cohort of 80 patients with mCRPC profiled by tumor RNA sequencing. Conclusions: We analyzed a total of 220 patients with mCRPC, validated the importance of cell-free AR/enhancer alterations as a prognostic biomarker in lethal mCRPC, and showed that the underlying mechanism for lethality involves reprogramming developmental states toward increased stemness. See related commentary by Nawfal et al., p. 7

Abstract Background: Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor (AR)-targeted agents is often lethal. Unfortunately, biomarkers for this deadly disease remain under investigation, and underpinning mechanisms are ill-understood. Here, we applied epigenomic approaches to identify cell-free (cf) DNA features associated with these lethal AR-resistant mCRPC patients. Methods: Plasma from 99 mCRPC patients was collected from two independent institutions, either prior to first-line AR-targeted therapy (n = 63) or during treatment (n = 36). We applied EnhanceAR-Seq to detect alterations in the AR locus including the upstream enhancer. We then applied genome-wide cfDNA epigenomics to 43 pre-treatment plasma samples and performed both fragmentomics and methylation sequencing to delineate the biology of AR/enhancer-altered lethal mCRPC. Finally, we identified stemness genes by applying CytoTRACE to external single-cell RNA-seq data from mCRPC, and applied these to cfDNA via promoter-level methylation analysis. We then assessed these stem-associated genes in tumor bulk RNA-seq data from an independent mCRPC cohort. Results: AR/enhancer alterations were detected in 44% of pre-treatment cfDNA samples and correlated with significantly worse progression-free survival (PFS) (HR = 2.12, P = 0.01) and overall survival (OS) (HR = 2.48; P = 0.02). AR/enhancer alterations were detected in 19% of on-treatment cfDNA samples and were also associated with profoundly worse PFS (HR = 15.38, P = 0.0003) and OS (HR = 15.53, P = 0.002). Further, genome-wide cfDNA fragmentomic analysis revealed that binding sites for developmentally relevant transcription factors such as FOXA1, NKX-3 and HOXB13 were significantly more accessible in AR/enhancer-altered lethal mCRPC patients. Gene set enrichment analysis of the 20 most accessible transcription factors revealed several pathways associated with stem cell development. Interestingly, the 10 most stem-associated genes obtained using CytoTRACE applied to an external cohort of single-cell RNA-seq data from mCRPC (He et al. 2021) were promoter-hypomethylated in the plasma cfDNA of AR/enhancer-altered lethal mCRPC patients (P = 0.005). Further, metagene analysis of these stem-associated genes revealed worse PFS by Kaplan-Meier analysis (P = 0.04), which was confirmed by multivariate cox proportional hazards modeling. Our stemness metagene further stratified overall survival (P = 0.03) when applied to an external cohort of 80 mCRPC patients from Abida et al. 2019. Conclusions: Alterations in AR including the upstream enhancer are remarkably prognostic in mCRPC with pre-treatment predictive potential. The mechanism can be inferred via cfDNA epigenomic analysis and appears to involve increased tumor cell stemness in AR-altered lethal mCRPC. Citation Format: Pradeep Singh Chauhan, Irfan Alahi, Savar Sinha, Ryan Mueller, Alexander L. Shiang, Jace Webster, Ha X. Dang, Debanjan Saha, Lilli Greiner, Breanna Yang, Elisa M. Ledet, Ramandeep K. Babbra, Wenjia Feng, Peter K. Harris, Faridi Qaium, Ellen B. Jaeger, Patrick J. Miller, Sydney A. Caputo, Oliver A. Sartor, Russell K. Pachynski, Christopher A. Maher, Aadel A. Chaudhuri. Epigenomic analysis of plasma cell-free DNA identifies stemness features associated with worse prognosis in AR-altered lethal metastatic castration resistant prostate cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 983.