We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
How do we estimate survival? External validation of a tool for survival estimation in patients with metastatic bone disease—decision analysis and comparison of three international patient populations
BMC Cancervolume 15, Article number: 424 (2015)
We recently developed a clinical decision support tool, capable of estimating the likelihood of survival at 3 and 12 months following surgery for patients with operable skeletal metastases. After making it publicly available on www.PATHFx.org, we attempted to externally validate it using independent, international data.
We collected data from patients treated at 13 Italian orthopaedic oncology referral centers between 2010 and 2013, then applied to PATHFx, which generated a probability of survival at three and 12-months for each patient. We assessed accuracy using the area under the receiver-operating characteristic curve (AUC), clinical utility using Decision Curve Analysis (DCA), and compared the Italian patient data to the training set (United States) and first external validation set (Scandinavia).
The Italian dataset contained 287 records with at least 12 months follow-up information. The AUCs for the three-month and 12-month estimates was 0.80 and 0.77, respectively. There were missing data, including the surgeon’s estimate of survival that was missing in the majority of records. Physiologically, Italian patients were similar to patients in the training and first validation sets. However notable differences were observed in the proportion of those surviving three and 12-months, suggesting differences in referral patterns and perhaps indications for surgery.
PATHFx was successfully validated in an Italian dataset containing missing data. This study demonstrates its broad applicability to European patients, even in centers with differing treatment philosophies from those previously studied.
Estimating survival in patients with skeletal metastases is important to set patient and physician expectations, as well as to guide surgical decision making [1–4]. During the preoperative evaluation, survival estimates can help surgeons carefully avoid under- or overtreatment of the disease by identifying which patients are likely to benefit from surgery but also whether a more durable implant may be necessary [5, 6]. Though most physicians are able to derive subjective survival estimates, they are generally inaccurate, and treating surgeons may be uncomfortable recording them in the medical record, or communicating them directly to patients .
With this in mind, we developed a Bayesian Belief Network capable of estimating three and twelve month survival in patients undergoing surgery for skeletal metastases . It is intended to guide, not replace, good clinical judgment and uses prognostic variables previously demonstrated to risk stratify these patients, including the oncologic diagnosis [8, 9], the extent of disease , the patient’s performance status , and basic laboratory assessments . In addition to delivering the likelihood of survival, PATHFx also estimates the quality of evidence supporting that estimation, which can be used by the treating surgeon to qualify each estimate. Following its development, we ensured that the tool was suitable for the clinical setting by performing decision curve analysis [13, 14], externally validated it using Scandinavian registry data , and made it publicly available to the international community, without charge, on www.PATHFx.org. However, the success of this tool depends on its performance in a variety of cultures, patient populations and institutions that may have differing treatment philosophies from those previously studied. In addition, though PATHFx was designed using the records of patients with metastases of the appendicular and axial skeleton, the first external validation set lacked patients treated for axial lesions. As such, additional validation studies are needed that include patients with both appendicular and axial metastases.
The purpose of this study was to (1) externally validate the PATHFx tool in an Italian patient population by evaluating accuracy by ROC analysis and clinical utility using DCA, and (2) compare the distributions of patients to both the training set (U.S.) and first external validation (Scandinavian) datasets, respectively.
The Italian Society of Orthopaedic and Traumatology (SIOT) established the Bone Metastasis Study Group in order to study patients with bone metastases and improve treatment. In the current study we retrospectively reviewed the records of 287 patients from 2010 to 2013 treated at one of thirteen Italian referral centers. Each record contained the 17 demographic and clinical variables, required of the PATHFx models. Survival was defined as the time elapsed from the date of surgery to the date of death or last follow-up. All records had sufficient follow-up to establish 12-month survival. This study received local ethical approval from the Università Campus Bio-Medico di Roma (Prot:15/13 19 June 2013). Informed consent was not required prior to using de-identified registry data.
Though data was collected from multiple Italian centers, the indications for surgery were standardized. In general, patients with metastatic disease of the extremities were offered surgery to prevent or treat a pathologic fracture, according to the Mirel’s criteria . The surgical indications for spine metastasis were: intractable pain, the onset of neurological deficits, caused by the compression of the myeloradicular structures of the neoplastic mass or by a pathological fracture of the vertebra, the mechanical instability of the spinal segment affected by the metastasis that causes a disabling mechanical pain and/or a neurological deficit and a failure of the previous therapy .
The PATHFx models are Bayesian Belief Networks comprised of ten prognostic features . These include: age at the time of surgery, sex, indication for surgery (impending or completed pathologic fracture), number of bone metastases (solitary or multiple), surgeon’s estimate of survival (postoperatively, in months), presence or absence of visceral metastases, presence or absence of lymph node metastases, preoperative hemoglobin concentration (g/dL, on admission to the hospital, prior to transfusion, if applicable), absolute lymphocyte count (K/μL), and the patient’s primary oncologic diagnosis, classified into one of three groups as previously described . For example, lung, gastric, and hepatocellular carcinoma and melanoma were assigned to Group 1; sarcomas and other carcinomas, Group 2; and breast, prostate, renal cell, and thyroid carcinoma, multiple myeloma, and malignant lymphoma, Group 3.
The definitions used for this study were similar to those previously described . Briefly, an impending pathologic fracture was one in which the degree of bone and/or cortical disruption warranted, in the opinion of the treating surgeon, prophylactic surgical stabilization to prevent fracture. Lesions that resulted in a change in bone length, alignment, rotation, or loss of height as determined by imaging, were considered completed pathologic fractures. Biopsy-proven and/or clinically obvious metastases to organs within the chest, abdomen or brain were considered visceral metastases. Only biopsy-proven metastases to the lymph nodes were considered indicative of lymph node involvement.
Using commercially available software (FasterAnalytics, DecisionQ Corp., Washington, DC, USA), we applied data contained in the Italian validation set to the to PATHFx, which estimated the likelihood of postoperative survival at three and 12 months, for each record. We then performed Receiver-operating characteristic (ROC) curve analysis and calculated the area under the ROC curve (AUC) as a measure of accuracy. The models were used “as-is” and were not re-fit or otherwise improved using either the Scandinavian or Italian validation sets. Bayesian Belief Networks retain functionality in the presence of missing data so no other imputation methods were employed. Validation was considered successful if the AUC was greater than 0.70 and was determined a priori. We chose this threshold because the authors consider it to be the lowest acceptable limit, however, Decision Curve Analysis (DCA)  was performed to determine whether the models should be used clinically.
The characteristics of the Italian set were compared to those of the training set and first external validation set. Continuous variables were tabulated and presented as mean (standard deviation), median (interquartile range) and categorical variables as number (%) (Table 1). The distribution of each continuous variable was compared with the normal distribution using the Shapiro-Wilk test. Equality of variance for continuous variables was determined using the Brown-Forsythe and Levene test. Statistical differences between continuous variables versus the bivariate outcome variables were evaluated using the Mann–Whitney U-test and the post hoc Tukey-Kramer assessment. Categorical variables were also tabulated and associations compared using Fisher’s exact test or chi-square analysis, depending on the number of expected values in the contingency matrix. A two-tailed α of 0.05 was considered statistically significant. We used JMP® Version 9.0.2 (SAS Institute, Inc, Cary, NC, USA) and R© Version 3.0.2 (R Foundation for Statistical Computing, Vienna, Austria) for all statistical estimations.
Two-hundred eighty seven (287) records had adequate follow-up information to establish survival at 3 and 12 months postoperatively and thus comprised the validation set. None of these records were excluded.
PATHFx correctly classified three-month survival in 253 of 287 (88 %) patients, and 12-month survival in 199 of 287 (69 %) patients. On ROC curve analysis, the AUCs were 0.80 and 0.77, respectively, for the three and 12-month survival, respectively.
Decision analysis revealed that PATHFx should be used, rather than assume all patients or no patients would ultimately survive longer than 12 months. However, since 93 % of Italian patients survived longer than three months, DCA indicated that outcomes may be better if orthopaedic surgeons assumed all patients would survive three months, rather than use the three month model (Fig. 1).
As expected, the demographic and clinical features of patients in the validation set differed from those of patients in the U.S. training set. Several features differed significantly (p < 0.05) including, presence of visceral and lymph node metastases, number of bone metastases, and three and 12-month survival. Nonsignificant differences were observed in age at surgery, gender, preoperative hemoglobin concentration, absolute lymphocyte count, oncologic diagnosis grouping, pathologic fracture status, ECOG performance status, and the surgeon’s estimate of survival. When compared to the Scandinavian set, most features differed significantly (p < 0.05) with the exception of gender, preoperative hemoglobin concentration, absolute lymphocyte count and the presence of visceral metastases. Most features in the validation set had some degree of missing data, also summarized in Tables 1 and 2. Notable features included the surgeon’s estimate of survival (missing in 87 %), absolute lymphocyte count (missing in 23 %), and ECOG performance (missing in 20 %), all of which are important first- or second-degree predictors of survival in the PATHFx tool.
We successfully externally validated PATHFx in an Italian patient population including patients with both axial and appendicular metastases. In doing so, we confirm the model’s ability to estimate the likelihood of survival at two time points useful for orthopaedic surgical decision-making. This is the second external validation study and demonstrates the model is also generalizable to the Italian patient population.
When one considers the goals of treating patients with skeletal metastases are to relieve pain and to restore function for the maximum amount of time, careful estimates of survival, such as those provided by PATHFx, are necessary to avoid over- or undertreatment of the disease. For example, if a surgeon considers nonoperative treatment, a very low probability of survival at three months may support this decision. By extension, if a surgeon were to consider using a less invasive and less durable implant such as an intramedullary nail, longer estimates of survival such as 3–12 months would support this decision. Conversely, estimates of survival greater than one year may support the decision to use a more durable implant such as a prosthesis in the case of extremity tumors, or more complicated spine procedures including vertebrectomies and combined anterior and posterior techniques. In fact, this study represents the first external validation of PATHFx in patients with axial (n = 34, 12 %), as well as appendicular skeletal metastases. This is important because although the training set contained 33 (18 %) spine patients, the Scandinavian external validation set contained only patients with extremity metastases.
Though there are several prognostic scoring systems designed for spine patients [10, 11, 18], none provide the surgeon with an estimation of the likelihood of survival at three and 12 months, which the authors consider to be useful for surgical decision-making. In addition, a recent analysis of seven prognostic tools demonstrated the Modified Bauer method  to be most reliable . PATHFx codifies the presence of visceral metastases, number of skeletal metastases and diagnosis grouping, which are all used by the modified Bauer method, and may explain why both models function accurately, in this setting. Nevertheless, it may be important to consider neurologic impairment in patients with spine metastases, as recommended by Tokuhashi . However, ECOG performance status which is also used by PATHFx may be an acceptable surrogate, since it is prognostic in patients with both axial  and appendicular metastases  notwithstanding the obvious differences in impairment due to neurologic as opposed to end-stage metastatic bone involvement. Still, we recognize the importance of a tool useful in the treatment of all patients with skeletal metastases—not simply those with spine or appendicular involvement. As such, the performance or applicability of PATHFx in the present validation set is encouraging.
PATHFx performed well in the Italian patient population, despite significant differences when compared to the training and previous validation sets (Tables 1 and 2). Importantly, 93 % of Italian patients survived longer than 3 months, which is much higher than either of the two previously studied groups. This may represent a key difference in the Italian patients, or more likely treatment philosophy and patient selection when compared to the U.S. and Scandinavian centers.
In addition, 63 % of Italian patients survived more than 12 months. This is twice as many as in the Scandinavian validation set, and nearly twice that observed in the training set. This, too may indicate key differences in patient selection and is surprising since there was a similar proportions of patients with pathologic fractures (p = 0.08), more-favorable diagnosis group (Group 3) (p = 0.42) and good performance status (ECOG 0,1,2) (p = 0.39) when compared to the training set. Still this may be explained by referral patterns among the Italian centers. Italian oncologists typically refer patients with excellent prognoses for orthopaedic consultation. In patients with more extensive disease and less favorable prognoses, however, surgery may be deemed unsuitable in the eyes of the oncologist, which obviates the need for an orthopaedic opinion. However, this practice may exclude patients that may benefit from less invasive stabilization or palliative procedures [20–22].
Nearly half of Italian patients included in this study presented with a solitary skeletal metastasis. This was unexpected, given that this proportion is much higher than both the training and previous validation sets, and could represent more effective disease surveillance practices than those in Scandinavia or the U.S. However, given the differences in referral patterns discussed above, it is more likely that Italian patients with less favorable prognoses—especially in the setting of impending pathologic fractures—were not referred for surgical management.
The accuracy for the three and 12-month models was 0.80 and 0.77, respectively. When compared to the original cross-validation AUCs of 0.86 and 0.83 , this represents a non-trivial, but acceptable 0.06-point degradation in model accuracy and is similar to that observed following external validation in the Scandinavian set (0.79 and 0.76, respectively) . Still by maintaining accuracy in differing patient populations, we believe PATHFx is sufficiently robust, and DCA suggests it may be used clinically, while undergoing additional external validation in more diverse patient populations.
The PATHFx models were designed to help surgeons avoid overtreatment or undertreatment of skeletal metastases. Previous work demonstrated that the models were suitable for clinical use, and that overly optimistic or pessimistic estimates generated by PATHFx were of unequal clinical significance . This is perhaps most important in the three month model that was designed to help surgeons identify which patients may benefit from a surgical or nonsurgical course of treatment. The present study demonstrated 34 (12 %) of records were misclassified by the three-month model. Of these, survival was overestimated in 13 (5 %) records, representing the maximum number of potentially unnecessary surgeries performed at the end of life. However, this estimate should be considered the theoretical maximum, since it likely includes patients who met surgical criteria and died of complications unrelated to the progression of disease. These results are more accurate than those observed in the Scandinavian set in which three month survival was overestimated in 15 % of records  and may be due to the larger proportion of Italian patients who survived greater than three months. If we consider that between 6 and 23 % of patients die within six weeks of surgery [23–25], then the clinical impact of such overestimates may fall within the acceptable norm.
Though one may consider an AUC of 0.8 for the three-month model to be sufficiently accurate, decision analysis helps illustrate the clinical impact of applying the model to a population in which virtually every patient referred for orthopaedic management of metastatic bone disease survives three months. Following DCA, we observe that at threshold probabilities (the point at which surgeons become indecisive about whether to offer surgery) less than 15 %, the model is equivalent to one in which all patients are expected to survive greater than three months. At thresholds >90 %, the three-month model should result in better outcomes. However, at thresholds between 15 and 90 %, an Italian orthopaedic surgeon is better off treating patients as if all will survive more than three-months, rather than use the three-month model. In the latter case, an erroneous underestimate may prompt the surgeon to withhold surgery from one in ten patients in whom it was otherwise indicated.
By extension, the 12-month model was designed to support decisions surrounding the type of procedure, as well as implant durability required for each patient. Of the 88 records misclassified by the 12-month model, survival was underestimated in 44 (15 %) cases. This represents the theoretical maximum proportion of patients at risk for implant failure if a less durable implant were used. This is higher than that observed in the Scandinavian validation set , in which 12-month survival was underestimated in 7.6 % of records. Though long term follow-up data were not available for this study beyond 12 months, we expect the proportion of patients surviving greater than 24 and 36 months, to decrease considerably. This trend has been observed previously [6, 8, 15], and further decreases the theoretical number of implants at risk, over time.
One of the most salient features of PATHFx is the ability to function in the presence of missing data. This attribute is particularly important when considering the surgeon’s estimate of survival was missing in 87 % of Italian and 100 % of Scandinavian records. The models maintained their accuracy because BBNs encode the information contained within the surgeon’s estimate in terms of shared, probabilistic relationships with other features, allowing one to “export” palliative expertise into settings where it may not exist. Though caution should be used when entering the surgeon’s estimate, those who are unsure of their estimate—or experience level—may simply leave it blank. Doing so will maintain accuracy of the model, while not introducing undue bias.
This study has several limitations. First, we developed PATHFx using the records of patients who underwent orthopaedic surgery for their skeletal metastases. Thus, it may not be applicable to all patients with metastatic disease, especially those in who are treated non-operatively. Second, similar to the previous Scandinavian external validation set, the Italian patient population was relatively homogeneous. However, we sought to obtain a representative sampling of Italian patients by collecting data from thirteen centers. Additionally, it is possible that PATHFx may become more accurate, by including other features potentially associated with survival in this patient population such as alkaline phosphatase , N-telopeptide [27, 28], and C-Reactive protein , or the degree of neurologic impairment as suggested by Tokuhashi. In addition, the time points chosen for PATHFx (three and 12-months) were initially chosen by two of the authors (JAF and JHH) because they are useful for orthopaedic surgical decision-making. Based on a recent study of practice patterns , other time points such as one-month and six-month survival may be needed in addition to the three and 12-month estimates to help surgeons decide on an operative strategy. One and six-month models are currently under development and would allow for a direct comparison with existing tools to estimate survival in patients with axial metastases, such as the Tokuhashi method . Next, the degree of experience required by surgeons to provide useful, as opposed to confounding, surgeon’s estimates is under further study, as the present study is too small to derive any meaningful information. Finally, PATHFx is a clinical decision support tool and should not supplant good clinical judgment by the treating surgeon and clinical team. Palliative surgery, by definition, can be appropriate even in patients with very short life expectancies, and low estimates of survival generated by any prognostic tool should not be used to deny these types of interventions if otherwise clinically indicated.
In conclusion, we successfully validated PATHFx using an Italian dataset containing patients with axial and appendicular skeletal metastases. This is the second external validation study and demonstrates that the tool is suited for clinical use in Italy. However, the three-month model should be used with caution in an Italian population, wherein nearly all (93 %) of patients referred for orthopaedic management of skeletal metastases are likely to survive longer than three months. Prospective, multicenter validation is necessary to confirm utility in other diverse patient populations and clinical settings, over time, which will provide an opportunity to assess whether the addition of newer, potentially prognostic variables could increase accuracy.
Wedin R, Bauer HC, Wersäll P. Failures after operation for skeletal metastatic lesions of long bones. Clin Orthop Relat Res. 1999:128–39.
Bauer HCF. Controversies in the surgical management of skeletal metastases. J Bone Joint Surg (Br). 2005;87:608–17.
Wedin R, Bauer H. Surgical treatment of skeletal metastatic lesions of the proximal femur: endoprosthesis or reconstruction nail? J Bone Joint Surg Br. 2005;87:1653–7.
Quinn RH, Randall RL, Benevenia J, Berven SH, Raskin KA. Contemporary management of metastatic bone disease: tips and tools of the trade for general practitioners. Instr Course Lect. 2014;63:431–41.
Steensma M, Boland PJ, Morris CD, Athanasian E, Healey JH. Endoprosthetic treatment is more durable for pathologic proximal femur fractures. Clin Orthop Relat Res. 2012;470:920–6.
Forsberg JA, Wedin R, Bauer H. Which implant is best after failed treatment for pathologic femur fractures? Clin Orthop Relat Res. 2013;471:735–40.
Lamont EB, Christakis NA. Prognostic disclosure to patients with cancer near the end of life. Ann Intern Med. 2001;134:1096–105.
Forsberg JA, Eberhardt J, Boland PJ, Wedin R, Healey JH. Estimating survival in patients with operable skeletal metastases: an application of a bayesian belief network. PLoS One. 2011;6, e19956.
Katagiri H, Takahashi M, Wakai K, Sugiura H, Kataoka T, Nakanishi K. Prognostic factors and a scoring system for patients with skeletal metastasis. J Bone Joint Surg (Br). 2005;87:698–703.
Bauer H, Wedin R. Survival after surgery for spinal and extremity metastases: prognostication in 241 patients. Acta Orthop. 1995;66:143–6.
Tokuhashi Y, Matsuzaki H, Oda H, Oshima M, Ryu J. A revised scoring system for preoperative evaluation of metastatic spine tumor prognosis. Spine. 2005;30:2186–91.
Hansen BH, Keller J, Laitinen M, Berg P, Skjeldal S, Trovik C, et al. The Scandinavian sarcoma group skeletal metastasis register. Survival after surgery for bone metastases in the pelvis and extremities. Acta Orthop Scand Suppl. 2004;75:11–5.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. 2006;26:565–74.
Forsberg JA, Sjoberg D, Chen Q-R, Vickers A, Healey JH. Treating metastatic disease: which survival model is best suited for the clinic? Clin Orthop Relat Res. 2013;471:843–50.
Forsberg JA, Wedin R, Bauer HCF, Hansen BH, Laitinen M, Trovik CS, et al. External validation of the Bayesian estimated tools for survival (BETS) models in patients with surgically treated skeletal metastases. BMC Cancer. 2012;12:493.
Mirels H. The classic: metastatic disease in long bones a proposed scoring system for diagnosing impending pathologic fractures. Clin Orthop Relat Res. 2003;415:S4.
Gasbarrini A, Boriani S, Capanna R, Casadei R, Di Martino A, Silvia Spinelli M, et al. Italian orthopaedic society bone metastasis study group: management of patients with metastasis to the vertebrae: recommendations from the Italian orthopaedic society (SIOT) bone metastasis study group. Expert Rev Anticancer Ther. 2014;14:143–50.
Hirabayashi H, Ebara S, Kinoshita T, Yuzawa Y, Nakamura I, Takahashi J, et al. Clinical outcome and survival after palliative surgery for spinal metastases: palliative surgery in spinal metastases. Cancer. 2003;97:476–84.
Wibmer C, Leithner A, Hofmann G, Clar H, Kapitan M, Berghold A, et al. Survival analysis of 254 patients after manifestation of spinal metastases. Spine. 2011;36:1977–86.
Yazawa Y, Frassica FJ, Chao EY, Pritchard DJ, Sim FH, Shives TC. Metastatic bone disease. A study of the surgical treatment of 166 pathologic humeral and femoral fractures. Clin Orthop Relat Res. 1990:213–19.
Nielsen OS, Munro AJ, Tannock IF. Bone metastases: pathophysiology and management policy. J Clin Oncol. 1991;9:509–24.
Perez CA, Bradfield JS, Morgan HC. Management of pathologic fractures. Cancer. 1972;29:684–93.
Clohisy DR, Le CT, Cheng EY, Dykes DC, Thompson RC. Evaluation of the feasibility of and results of measuring health-status changes in patients undergoing surgical treatment for skeletal metastases. J Orthop Res. 2000;18:1–9.
Quinn R, Drenga J. Perioperative morbidity and mortality after reconstruction for metastatic tumors of the proximal femur and acetabulum. J Arthroplast. 2006;21:227–32.
Bollen L, van der Linden YM, Pondaag W, Fiocco M, Pattynama BPM, Marijnen CAM, et al. Prognostic factors associated with survival in patients with symptomatic spinal bone metastases: a retrospective cohort study of 1,043 patients. Neuro-Oncology. 2014;16:991–8.
Halabi S, Lin C-Y, Kelly WK, Fizazi KS, Moul JW, Kaplan EB, et al. Updated prognostic model for predicting overall survival in first-line chemotherapy for patients with metastatic castration-resistant prostate cancer. J Clin Oncol. 2014;32:671–7.
Brown JE, Thomson CS, Ellis SP, Gutcher SA, Purohit OP, Coleman RE. Bone resorption predicts for skeletal complications in metastatic bone disease. Br J Cancer. 2003;89:2031–7.
Brown JE, Cook RJ, Major P, Lipton A, Saad F, Smith M, et al. Bone turnover markers as predictors of skeletal complications in prostate cancer, lung cancer, and other solid tumors. J Natl Cancer Inst. 2005;97:59–69.
Bataille R, Boccadoro M, Klein B, Durie B, Pileri A. C-reactive protein and beta-2 microglobulin produce a simple and powerful myeloma staging system. Blood. 1992;80:733–7.
Steensma M, Healey JH. Trends in the surgical treatment of pathologic proximal femur fractures among musculoskeletal tumor society members. Clin Orthop Relat Res. 2013;471:2000–6.
We thank the surgeons and staff at each of the 13 Italian centers involved in this study.
Role of the funding source
The Società Italiana di Ortopedia e Traumatologia Bone Metastasis Study Group facilitated data collection from the thirteen participating centers. No outside funding was used to support the analysis, preparation or submission of this manuscript.
The authors declare that they have no competing interests. PATHFx was developed using publicly available software (FasterAnalytics, DecisionQ Inc. Washington DC, USA) in collaboration with the corresponding author (JAF), and the web-enabled user interface is a product of DecisionQ Inc.
Conception and design: JAF, AP, MSS, RW, JHH. Acquisition of data: MSS, AP, VI, PD, PR, GM, AG, RB, RP, FF, AL, ADM, FN, FC, MAR, DAC, VD, RC. Analysis and interpretation of data: JAF, AP, MSS, VI, PAD, PR, GM, AG, RB, RP, FF, AL, ADM, FN, FC, MAR, DAC, VD, RC, JHH. Literature searches: JAF, MSS, AP. Drafting of manuscript: JAF, MSS, RW. Critical revision: JAF, MSS, RW, JHH. Statistical expertise: JAF. Obtaining funding: JAF, AP, MSS, VI, PAD, PR, GM, AG, RB, RP, FF, AL, ADM, FN, FC, MAR, DAC, VD, RC. Supervision: JAF, AP, RW. All authors read and approved the final manuscript.