Bayesian Nonparametric Survival Regression for Optimizing Precision Dosing of Intravenous Busulfan in Allogeneic Stem Cell Transplantation Yanxun Xu† Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, USA Peter F. Thall Department of Biostatistics, U.T. M.D. Anderson Cancer Center, Houston, USA. William Hua Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, USA Borje S. Andersson Department of Stem Cell Transplantation and Cellular Therapy, U.T. M.D. Anderson Cancer Cen- ter, Houston, USA. Abstract. Allogeneic stem cell transplantation (allo-SCT) is now part of standard of care for acute leukemia (AL). To reduce toxicity of the pre-transplant conditioning regimen, intravenous busulfan is usually used as a preparative regimen for AL patients undergoing allo-SCT. Sys- temic busulfan exposure, characterized by the area under the plasma concentration versus time curve (AUC), is strongly associated with clinical outcome. An AUC that is too high is associated with severe toxicities, while an AUC that is too low carries increased risks of disease recurrence and failure to engraft. Consequently, an optimal AUC interval needs to be determined for ther- apeutic use. To address the possibility that busulfan pharmacokinetics and pharmacodynamics vary significantly with patient characteristics, we propose a tailored approach to determine op- timal covariate-specific AUC intervals. To estimate these personalized AUC intervals, we apply a flexible Bayesian nonparametric regression model based on a dependent Dirichlet process and Gaussian process, DDP-GP. Our analyses of a dataset of 151 patients identified optimal therapeutic intervals for AUC that varied substantively with age and whether the patient was in complete remission or had active disease at transplant. Extensive simulations to evaluate the DDP-GP model in similar settings showed that its performance compares favorably to alternative methods. We provide an R package, DDPGPSurv, that implements the DDP-GP model for a broad range of survival regression analyses. Keywords: Allogeneic stem cell transplantation; Bayesian nonparametrics; Personalized medicine; Survival regression. †Address for correspondence: Yanxun Xu, Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, MD, USA 21042. Email: [email protected].
23
Embed
Bayesian Nonparametric Survival Regression for Optimizing ...pfthall/main/JRSSC_2019_AUC-surv.pdfinferences in settings where the proportional hazards assumption, speci c parametric
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Bayesian Nonparametric Survival Regression for OptimizingPrecision Dosing of Intravenous Busulfan in Allogeneic StemCell Transplantation
Yanxun Xu†Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, USA
Peter F. Thall
Department of Biostatistics, U.T. M.D. Anderson Cancer Center, Houston, USA.
William Hua
Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, USA
Borje S. Andersson
Department of Stem Cell Transplantation and Cellular Therapy, U.T. M.D. Anderson Cancer Cen-
ter, Houston, USA.
Abstract. Allogeneic stem cell transplantation (allo-SCT) is now part of standard of care foracute leukemia (AL). To reduce toxicity of the pre-transplant conditioning regimen, intravenousbusulfan is usually used as a preparative regimen for AL patients undergoing allo-SCT. Sys-temic busulfan exposure, characterized by the area under the plasma concentration versus timecurve (AUC), is strongly associated with clinical outcome. An AUC that is too high is associatedwith severe toxicities, while an AUC that is too low carries increased risks of disease recurrenceand failure to engraft. Consequently, an optimal AUC interval needs to be determined for ther-apeutic use. To address the possibility that busulfan pharmacokinetics and pharmacodynamicsvary significantly with patient characteristics, we propose a tailored approach to determine op-timal covariate-specific AUC intervals. To estimate these personalized AUC intervals, we applya flexible Bayesian nonparametric regression model based on a dependent Dirichlet processand Gaussian process, DDP-GP. Our analyses of a dataset of 151 patients identified optimaltherapeutic intervals for AUC that varied substantively with age and whether the patient was incomplete remission or had active disease at transplant. Extensive simulations to evaluate theDDP-GP model in similar settings showed that its performance compares favorably to alternativemethods. We provide an R package, DDPGPSurv, that implements the DDP-GP model for a broadrange of survival regression analyses.
Figure 5: Estimated survival functions under the DDP-GP survival regression model for patients with
different CR status (Yes or No) and ages (30, 40, 50, 60). The patients are assigned AUC=5. The
dashed lines represent the point-wise 95% credible intervals for each survival curve.
are given in Figure 6. Our analyses confirm the existence, for each combination of CR status
and Age, of an optimal AUC range that yields higher expected survival times compared to an
AUC that is either below or above the optimal range. A very important inference is that these
optimal AUC ranges differ substantially between many of the (CR status, Age) combinations.
This has extremely important therapeutic implications when choosing an individual patient’s
targeted AUC. For example, the optimal AUC interval for a patient not in CR with Age=50
is 4.7 ± 0.47 = [4.23, 5.17] compared with the optimal interval 5.8 ± 0.58 = [5.22, 6.38] for
a patient in CR with Age=40. Since these intervals are disjoint, they suggest that these two
patients should have very different targeted AUC values to maximize their expected survival
times. The estimated mean survival times versus AUC under the alternative methods, PT,
TBP, RF, and BART, are included in the supplement Section 3. There are no meaningful
patterns we can observe in these figures.
In contrast with our inferences, (Bartelink et al., 2016) concluded that CR status has a
negligible effect on the optimal AUC However, the results reported by Bartelink et al. (2016)
were based on data from a large number of different medical centers, many different pretrans-
plant conditioning regimens were used, the PK-data were obtained from different laboratories,
with a very heterogeneous pediatric patient population having a large number of different di-
agnostic categories, including patients with malignant and non-malignant genetic disorders. In
contrast, our analyses are based on a much more homogeneous dataset. Our results indicate
that CR status is an important covariate, and that the optimal dose of AUC is higher for
patients who are in CR at transplant. Furthermore, the increased optimal AUC for patients
in CR at transplant versus patients not in CR is much larger in older patients, whereas these
16
3 4 5 6 7
34
56
7
CR=No, Age=30
AUC
Mea
n S
urvi
val (
log)
AUC=6.1
3 4 5 6 7
34
56
7
CR=No, Age=40
AUC
Mea
n S
urvi
val (
log)
AUC=5.5
3 4 5 6 7
34
56
7
CR=No, Age=50
AUCM
ean
Sur
viva
l (lo
g)
AUC=4.7
3 4 5 6 7
34
56
7
CR=No, Age=60
AUC
Mea
n S
urvi
val (
log)
AUC=4.1
3 4 5 6 7
46
810
12
CR=Yes, Age=30
AUC
Mea
n S
urvi
val (
log)
AUC=6.2
3 4 5 6 7
46
810
12
CR=Yes, Age=40
AUC
Mea
n S
urvi
val (
log)
AUC=5.8
3 4 5 6 7
46
810
12
CR=Yes, Age=50
AUC
Mea
n S
urvi
val (
log)
AUC=5.4
3 4 5 6 7
46
810
12
CR=Yes, Age=60
AUC
Mea
n S
urvi
val (
log)
AUC=5.1
Figure 6: Mean log survival time estimates under the DDP-GP model, as a function of AUC,for each of eight (CR status, Age) combinations. The gray area in each plot represents the95% credible interval for estimated mean survival, and the tick marks on the horizontal axis(rug plot) indicate the AUC values for patients in the data set. The red area represents theoptimal AUC range, defined as the estimated mean ±10%.
17
differences appear negligible in adolescents or young adults, similar to what was reported by
Bartelink, et al. (2016). Our results also demonstrate that, across all ages, mean survival time
for patients in CR is larger compared with those not CR.
To further illustrate how the optimal AUC ranges change with both CR and Age, we plotted
the optimal AUC ranges as Age is varied continuously, for CR=Yes and CR=No, in Figure
7. The negative association between optimal AUC and Age is clearly shown by this figure. It
also shows that, while CR status has virtually no effect on the optimal AUC interval for very
young patients with Age ≤ 28, the optimal AUC for patients in CR at transplant is increasingly
higher as Age increases, with the optimal intervals for CR = Yes versus CR = No becoming
completely disjoint for patients above 55 years of age. Thus, the lower portions of the curves
in Figure 7 for Age ≤ 28, agree with the conclusion of (Bartelink et al., 2016) for pediatric and
adolescent patients, while the higher portions for Age > 28, provide news insights. Again, this
demonstrates the importance of considering both CR status and Age when planning a targeted
AUC for a patient with a diagnosis of AML or MDS.
20 30 40 50 60
34
56
7
Age
Opt
imal
AU
C
CR=YesCR=No
Figure 7: Optimal AUC ranges versus age given CR status. The blue and red lines representthe optimal AUC for CR=Yes and No, respectively. The optimal AUC ranges are representedby the shaded regions above and below the optimal AUC.
7. Conclusions
We have proposed an extended Bayesian nonparametric DDP-GP model for survival regres-
sion having a generalized covariance structure, studied it by simulation, and applied it to
estimate personalized optimal dose intervals for IV busulfan in allo-SCT for AML/MDS. Our
simulations, constructed to mimic the dataset, show that the DDP-GP model provides more
18
accurate survival function estimates and optimal AUC range estimates compared with conven-
tional parametric to AFT models. Our analyses of the IV busulfan allo-SCT dataset identified
optimal AUC intervals, varying with the patient’s CR status and Age, that previously have
not been known for this treatment. Our results may have profound therapeutic implications,
since they provide a basis for personalized medicine by enabling physicians to prospectively
assign an optimized therapeutic target interval for each patient based on his/her CR status
and age.
More generally, we have developed an R package, DDPGPSurv, that implements the DDP-GP
model for a broad range of survival regression analyses. While the DDP-GP is more complex
than conventional survival regression models, its robustness and broad applicability make it
an attractive methodology for survival analysis. The DDP-GP based data analysis reported
here, while important in its own right, identified a nonlinear three-way interaction between age,
CR status, and AUC in their joint effect on survival time, as shown by Figures 6 and 7. This
pattern was identified despite the fact, noted above, that only the main effects were included in
the mean of the Gaussian process prior via the linear term β0 +β1Age+β2CR+β3AUC. This
is because the DDP-GP is essentially a mixture model, hence it can identify complex patterns
in the data that may be missed by conventional models. For the allo-SCT IV busulfan data,
this may be related to the multi-modality of the survival time distribution, seen in Figure 1.
This illustrates the practical advantage that, when applying the DDP-GP, one need not guess
or search for complex patterns in the linear term of the covariates, as is done routinely when
applying conventional survival regression models.
Acknowledgements
Peter Thall’s research was supported by NCI grant 5-R01-CA083932.
References
Andersson, B., Thall, P., Valdez, B., Milton, D., Al-Atrash, G., Chen, J., Gulbis, A., Chu,
D., Martinez, C., Parmar, S. et al. (2016) Fludarabine with pharmacokinetically guided
IV busulfan is superior to fixed-dose delivery in pretransplant conditioning of AML/MDS
patients. Bone Marrow Transplantation.
Andersson, B. S., Madden, T., Tran, H. T., Hu, W. W., Blume, K. G., Chow, D. S.-L., Cham-
plin, R. E. and Vaughan, W. P. (2000) Acute safety and pharmacokinetics of intravenous
busulfan when used with oral busulfan and cyclophosphamide as pretransplantation condi-
tioning therapy: a phase i study. Biology of Blood and Marrow Transplantation, 6, 548–554.
Andersson, B. S., Thall, P. F., Madden, T., Couriel, D., Wang, X., Tran, H. T., Anderlini, P.,
De Lima, M., Gajewski, J. and Champlin, R. E. (2002) Busulfan systemic exposure relative to
regimen-related toxicity and acute graft-versus-host disease: defining a therapeutic window
for iv bucy2 in chronic myelogenous leukemia. Biology of Blood and Marrow Transplantation,
8, 477–485.
19
Bartelink, I. H., Bredius, R. G., Belitser, S. V., Suttorp, M. M., Bierings, M., Knibbe, C. A.,
Egeler, M., Lankester, A. C., Egberts, A. C., Zwaveling, J. et al. (2009) Association between
busulfan exposure and outcome in children receiving intravenous busulfan before hematologic
stem cell transplantation. Biology of Blood and Marrow Transplantation, 15, 231–241.
Bartelink, I. H., Lalmohamed, A., van Reij, E. M., Dvorak, C. C., Savic, R. M., Zwaveling, J.,
Bredius, R. G., Egberts, A. C., Bierings, M., Kletzel, M. et al. (2016) Association of busulfan
exposure with survival and toxicity after haemopoietic cell transplantation in children and
young adults: a multicentre, retrospective cohort analysis. The Lancet Haematology, 3,
e526–e536.
Bredeson, C., LeRademacher, J., Kato, K., DiPersio, J. F., Agura, E., Devine, S. M., Appel-
baum, F. R., Tomblyn, M. R., Laport, G. G., Zhu, X. et al. (2013) Prospective cohort study
comparing intravenous busulfan to total body irradiation in hematopoietic cell transplanta-
tion. Blood, 122, 3871–3878.
Chipman, H. A., George, E. I. and McCulloch, R. E. (2012) BART: Bayesian additive regression
trees. Annals of Applied Statistics, 6, 266–298.
Copelan, E. A., Hamilton, B. K., Avalos, B., Ahn, K. W., Bolwell, B. J., Zhu, X., Aljurf, M.,
Van Besien, K., Bredeson, C. N., Cahn, J.-Y. et al. (2013) Better leukemia-free and overall
survival in aml in first remission following cyclophosphamide in combination with busulfan
compared to tbi. Blood, 122, 3863–3870.
De Iorio, M., Johnson, W. O., Muller, P. and Rosner, G. L. (2009) Bayesian nonparametric