Predicting time to prostate cancer recurrence based on joint models for non-linear longitudinal biomarkers and event...

Biological markers that are both sensitive and specific for tumour regrowth or metastasis are increasingly becoming available and routinely monitored during the regular follow-up of patients treated for cancer. Obtained by a simple blood test, these markers provide an inexpensive non-invasive means for the early detection of recurrence (or progression). Currently, the longitudinal behaviour of the marker is viewed as an indicator of early disease progression, and is applied by a physician in making clinical decisions. One marker that has been studied for use in both population screening for early disease and for detection of recurrence in prostate cancer patients is PSA. The elevation of PSA levels is known to precede clinically detectable recurrence by 2 to 5 years, and current clinical practice often relies partially on multiple recent rises in PSA to trigger a change in treatment. However, the longitudinal trajectory for individual markers is often non-linear; in many cases there is a decline immediately following radiation therapy or surgery, a plateau during remission, followed by an exponential rise following the recurrence of the cancer. The aim of this article is to determine the multiple aspects of the longitudinal PSA biomarker trajectory that can be most sensitive for predicting time to clinical recurrence. Joint Bayesian models for the longitudinal measures and event times are utilized based on non-linear hierarchical models, implied by unknown change-points, for the longitudinal trajectories, and a Cox proportional hazard model for progression times, with functionals of the longitudinal parameters as covariates in the Cox model. Using Markov chain Monte Carlo sampling schemes, the joint model is fit to longitudinal PSA measures from 676 patients treated at Massachusetts General Hospital between the years 1988 and 1995 with follow-up to 1999. Based on these data, predictive schemes for detecting cancer recurrence in new patients based on their longitudinal trajectory are derived.