Understanding the intricacies of lethal prostate cancer poses specific challenges due to difficulties in accurate modeling of metastasis in vivo. Here we show that NPK EYFP mice (for Nkx3.1 CreERT2/+ ; Pten flox/flox ; Kras LSL-G12D/+ ; R26R-CAG-LSL-EYFP/+) develop prostate cancer with a high penetrance of metastasis to bone, thereby enabling detection and tracking of bone metastasis in vivo and ex vivo. Transcriptomic and whole-exome analyses of bone metastasis from these mice revealed distinct molecular profiles conserved between human and mouse and specific patterns of subclonal branching from the primary tumor. Integrating bulk and single-cell transcriptomic data from mouse and human datasets with functional studies in vivo unravels a unique MYC/RAS co-activation signature associated with prostate cancer metastasis. Finally, we identify a gene signature with prognostic value for time to metastasis and predictive of treatment response in human patients undergoing androgen receptor therapy across clinical cohorts, thus uncovering conserved mechanisms of metastasis with potential translational significance.