• VRd is superior to Rd in veterans with myeloma, and the evidence of benefit is strongest in those with the highest levels of frailty.

  • More intensive treatment of a patient who is frail with myeloma may be a more effective treatment of frailty itself.

Although randomized controlled trial data suggest that the more intensive triplet bortezomib-lenalidomide-dexamethasone (VRd) is superior to the less intensive doublet lenalidomide-dexamethasone (Rd) in patients newly diagnosed with multiple myeloma (MM), guidelines have historically recommended Rd over VRd for patients who are frail and may not tolerate a triplet. We identified 2573 patients (median age, 69.7 years) newly diagnosed with MM who were initiated on VRd (990) or Rd (1583) in the national US Veterans Affairs health care System from 2004 to 2020. We measured frailty using the Veterans Affairs Frailty Index. To reduce imbalance in confounding, we matched patients for MM stage and 1:1 based on a propensity score. Patients who were moderate-severely frail had a higher prevalence of stage III MM and myeloma-related frailty deficits than patients who were not frail. VRd vs Rd was associated with lower mortality (hazard ratio [HR], 0.81; 95% confidence interval [CI], 0.70-0.94) in the overall matched population. Patients who were moderate-severely frail demonstrated the strongest association (HR 0.74; 95% CI, 0.56-0.97), whereas the association weakened in those who were mildly frail (HR, 0.80; 95% CI, 0.61-1.05) and nonfrail (HR, 0.86; 95% CI, 0.67-1.10). VRd vs Rd was associated with a modestly higher incidence of hospitalizations in the overall population, but this association weakened in patients who were moderate-severely frail. Our findings confirm the benefit of VRd over Rd in US veterans and further suggest that this benefit is strongest in patients with the highest levels of frailty, arguing that more intensive treatment of myeloma may be more effective treatment of frailty itself.

The optimal intensity of treatment in adults who are frail with multiple myeloma (MM) is often uncertain, given the higher risk of harm in patients who are frail on standard, intensive regimens.1,2 Randomized controlled trials (RCTs) directly or indirectly exclude patients who are frail, limiting the generalizability of results to a large portion of patients treated in practice.3,4 Although analyses from RCTs and observational studies consistently show that frailty is associated with higher risks of treatment toxicity, treatment discontinuation, and mortality,5 this evidence does little to inform the key clinical question, that is, whether to administer to patients who are frail, in practice, the same intensity of treatments evaluated in the younger, fitter patients in RCTs.

An illustrative example is Southwest Oncology Group (SWOG) S0777, the seminal randomized phase 3 trial that demonstrated an overall survival benefit associated with a more intensive triplet (bortezomib-lenalidomide-dexamethasone or VRd) vs a less intensive doublet (lenalidomide-dexamethasone or Rd) regimen in patients newly diagnosed with MM not intended for immediate hematopoietic stem cell transplant.6 These findings solidified the recommendation to favor initiation on the triplet VRd over the doublet Rd in most patients. However, many experts and guidelines have historically recommended Rd over VRd for patients who are frail and may not be able to tolerate a triplet.7,8 To date, there has been no evidence suggesting that patients who are frail treated with VRd experience either lesser benefits and/or greater harms than those treated with Rd. Although there are newer triplet regimens that are preferred in older patients (eg, daratumumab-lenalidomide-dexamethasone9), examining the effects of VRd vs Rd provides an opportunity to examine whether patients who are frail, in practice, fare worse on a more intensive initial regimen.

We designed a retrospective cohort study in US veterans, a population that is even more underrepresented in MM clinical trials than the general population, to evaluate for differences in survival associated with VRd vs Rd in newly diagnosed MM. We hypothesized that, similar to the results of SWOG S0777, we would observe a survival benefit associated with VRd over Rd in the overall population. Furthermore, we hypothesized that, similar to assumptions behind expert recommendations for clinical practice, the survival benefit associated with VRd would diminish in patients with high levels of frailty than in relatively fit patients, whereas the risk of hospitalization would increase (reflecting higher toxicity).

Study design and population

Emulating SWOG S0777,6 we designed a retrospective cohort study using a new-user design that uses data from the Veterans Affairs (VA) Corporate Data Warehouse (CDW), which collects clinical, billing, and electronic health record (EHR) information from veterans treated in VA facilities throughout the United States.10 We also used data from the national VA Cancer Registry, which collects information on cancer diagnosis and treatment compiled and submitted by local cancer registry staff at each of the 132 VA Medical Centers that diagnose and treat veterans with cancer.11,12 Using a selection algorithm validated against clinician chart review,13 we included all veterans newly diagnosed with MM who initiated treatment in VA Medical centers with either VRd or Rd from years 2004 (date of the first instance of using Rd) to 2020 (see supplemental Methods regarding our selection algorithm). We excluded veterans who received transplant within the first four months of treatment initiation. Although we could not measure a clinician’s intention to transplant at the time of initiation of VRd or Rd, as was measured in SWOG S0777,6 we measured the receipt of transplant at any point during follow-up. This study was approved by the VA Boston Health Care System Institutional Review Board. The study was conducted in accordance with the Declaration of Helsinki.

Frailty, MM stage, and covariates

We measured frailty using the Veterans Affairs Frailty Index (VA-FI), an electronic frailty index developed by geriatricians that carries both predictive and construct validity, including when validated against clinical measures of frailty assessed by geriatricians and the International Myeloma Working Group Frailty Score (supplemental Table 2).14-16 In brief, the VA-FI measures frailty based on the deficit-accumulation approach, which is one of the most widely-studied models of frailty.17 A total of 31 aging-related health deficits were measured from the EHR data and administrative claims, using diagnostic and procedural codes spanning the domains of morbidity (14 deficits), function (8), cognition (3), sensory (3), and other (3). These health deficits were measured within a 3-year assessment period leading up to the date of treatment initiation. A score was calculated, for each patient, as the proportion of possible health deficits that are present in an individual, ranging from 0 to 1, in which higher values indicate more severe frailty.18 Based on validated cut-points,14,15,19-22 veterans were classified as nonfrail (VA-FI < 0.2), mildly frail (VA-FI = 0.2-0.3), or moderate-severely frail (VA-FI ≥ 0.3).

We measured MM stage from both laboratory and EHR data within the CDW (supplemental Methods). We further measured pretreatment covariates spanning sociodemographic variables (age, sex, race, and rurality of the home address of the veteran23), MM risk24 (hemoglobin, platelets, calcium, and creatinine, measured within 90 days before initial treatment and imputing missing values as a function of all other covariates using multivariable regression), comorbidity burden (Charlson Comorbidity Index25), polypharmacy (having ≥ 5 chronic medications), and health care use (number of unplanned hospitalizations, emergency department visits, and skilled nursing facility stays in the year before initial treatment).

Outcomes

Our primary outcome was overall survival, measured using vital status information in the VA CDW (which provides 98.3% sensitivity and 97.6% exact agreement against the US National Death Index).26 Our secondary outcome was unplanned hospitalizations at the VA Health care System with a duration of stay ≥1 day.27 For mortality, veterans were followed up through 5 years (mimicking the SWOG S0777 median follow-up time) until death or the end of the study period, 31 December 2021, after which their data were censored. For hospitalizations, veterans were followed up for 1 year until death, the end of study, or until their last inpatient/outpatient record in the CDW, after which their data were censored.

Analysis

We compared baseline characteristics between VRd and Rd treatment groups. To reduce confounding within frailty categories, we matched patients based on the MM stage and estimated a propensity score (PS) using regularized logistic regression that modeled the probability of being initially treated with VRd vs Rd as a function of all covariates. We did not find differences in survival over increasing calendar years of treatment initiation (supplemental Figure 1), so we did not include initial treatment year as a covariate in the PS model. We performed a 1:1 nearest-neighbor matching with a caliper of 1% of the PS.

For our primary analysis, we evaluated the association between VRd vs Rd and overall survival within each matched frailty category and in the combined overall population, first by using Kaplan-Meier analyses to evaluate differences in survival between VRd and Rd groups and next by using Cox proportional hazards regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs). We assessed for effect heterogeneity across frailty categories.28 For our secondary analysis, we evaluated the association between VRd vs Rd and the rates of unplanned hospitalizations. We conducted sensitivity analyses in multivariable Cox regression models, adjusting for any covariates that showed evidence of remaining imbalance (standardized mean difference [SMD] > 0.1) across those treated with VRd vs Rd after matching.29,30 We also repeated analyses with landmarking follow-up to begin on day 15 after the date of the first treatment with either VRd or Rd, accounting for potential immortal time bias.31 All analyses were consistent with an intent-to-treat approach, analyzing veterans according to their initial treatment assignment. Analyses were conducted in R 4.1.2 and Stata 17.0. We followed the STROBE guidelines in the design and reporting of observational studies.32 

Characteristics of the study population

We identified 2573 patients with newly diagnosed MM (median age, 69.7 years; interquartile range, 74.2-75.8) who were initiated either on VRd (990) or Rd (1583) from 2004 to 2020 (Figure 1). After matching those within each frailty category based on the stage and PS, all covariates were adequately balanced between the VRd and Rd groups in the combined overall population and within each frailty category, except for minimal residual imbalance in prior emergency department use in patients who were nonfrail (SMD = 0.103) and prior skilled nursing facility use in patients who were moderate-severely frail (SMD == 0.111; Table 1; supplemental Tables 4-6). Compared with patients who were nonfrail, patients who were moderate-severely frail were older (median age, 67.6 years vs 71.1 years) and had a higher prevalence of stage III MM (19.5% vs 32.9%) and myeloma-related health deficits, such as anemia (47.2% vs 85%), kidney disease (11.2% vs 46.8%), and bone disease or pathologic skeletal fractures (4.6% vs 18.5%). A total of 170 patients (21.6%) in the VRd group and 134 patients (17%) in the Rd group underwent a transplant at some point in follow-up. In patients who were moderate-severely frail, 15 (8.7%) in the VRd group and 14 (8.1%) in the Rd group underwent transplants.

Figure 1.

Flow diagram of the selection of the study population.

Figure 1.

Flow diagram of the selection of the study population.

Close modal
Table 1.

Characteristics of veterans with newly diagnosed MM initially treated with either VRd or Rd before and after matching based on stage and PS within each frailty category

Rd before matching
(N = 1583)
VRd before matching
(N = 990)
SMD before matchingRd after matching (N = 788)VRd after matching (N = 788)SMD after matching
Age, median (IQR), y 70.6 (65.1-78.1) 68.7 (63.2-72.6) 0.399 68.3 (62.9-74.3) 69.7 (64.4-73.5) 0.041 
Race/ethnicity, n (%)   0.27   0.034 
Non-Hispanic White 1035 (65.4) 526 (53.1)  469 (59.5) 467 (59.3)  
Black 362 (22.9) 334 (33.7)  215 (27.3) 223 (28.3)  
Hispanic 70 (4.4) 59 (6.0)  42 (5.3) 42 (5.3)  
Other 116 (7.3) 71 (7.2)  62 (7.9) 56 (7.1)  
Frailty, n (%)   0.043   < 0.001 
Nonfrail (VA-FI < 0.2) 760 (48.0) 463 (46.8)  389 (49.4) 389 (49.4)  
Mild frailty (VA-FI = 0.2-0.3) 439 (27.7) 294 (29.7)  226 (28.7) 226 (28.7)  
Moderate-severe frailty (VA- FI ≥ 0.3) 384 (24.3) 233 (23.5)  173 (22.0) 173 (22.0)  
Charlson Comorbidity Index, median (IQR) 1.00 (0.00-2.00) 1.00 (0.00-2.00) 0.063 1.00 (0.00-2.00) 1.00 (0.00-2.00) 0.010 
Myeloma stage, n (%)   0.361   < 0.001 
Stage I 349 (22.0) 91 (9.2)  76 (9.6) 76 (9.6)  
Stage II 350 (22.1) 241 (24.3)  200 (25.4) 200 (25.4)  
Stage III 527 (33.3) 392 (39.6)  318 (40.4) 318 (40.4)  
Missing 357 (22.6) 266 (26.9)  194 (24.6) 194 (24.6)  
Pretreatment laboratory levels, median (IQR)       
Creatinine (mg/dL) 1.20 (0.97-1.50) 1.18 (0.96-1.51) 0.013 1.15 (0.94-1.46) 1.15 (0.94-1.50) 0.026 
Calcium (mg/dL) 9.10 (8.70-9.60) 9.20 (8.80-9.70) 0.168 9.20 (8.80-9.66) 9.20 (8.80-9.70) 0.002 
Hemoglobin (g/dL) 11.10 (9.70-12.50) 10.90 (9.30-12.47) 0.066 11.01 (9.60-12.62) 11.10 (9.50-12.60) 0.010 
Platelets (n per mL) 205.00 (158.50-250.00) 202.00 (160.00-252.75) 0.014 203.00 (158.00-250.00) 204.42 (162.00-252.00) 0.029 
Rurality of veteran, n (%)   0.099   0.030 
Rural 561 (35.4) 306 (30.9)  263 (33.4) 261 (33.1)  
Urban 1012 (63.9) 679 (68.6)  518 (65.7) 522 (66.2)  
Unknown 10 (0.6) 5 (0.5)  7 (0.9) 5 (0.6)  
Polypharmacy (≥5 medications), n (%) 522 (33.0) 333 (33.6) 0.014 265 (33.6) 264 (33.5) 0.003 
Healthcare utilization in the year before treatment, median (IQR)       
Any hospital utilization 0.00 (0.00-1.00) 0.00 (0.00-1.00) 0.017 0.00 (0.00-1.00) 0.00 (0.00-1.00) 0.040 
Any ED utilization 0.00 (0.00-1.00) 1.00 (0.00-2.00) 0.309 0.00 (0.00-2.00) 1.00 (0.00-2.00) 0.060 
Any NH utilization 0.00 (0.00-0.00) 0.00 (0.00-0.00) 0.096 0.00 (0.00-0.00) 0.00 (0.00-0.00) 0.019 
Rd before matching
(N = 1583)
VRd before matching
(N = 990)
SMD before matchingRd after matching (N = 788)VRd after matching (N = 788)SMD after matching
Age, median (IQR), y 70.6 (65.1-78.1) 68.7 (63.2-72.6) 0.399 68.3 (62.9-74.3) 69.7 (64.4-73.5) 0.041 
Race/ethnicity, n (%)   0.27   0.034 
Non-Hispanic White 1035 (65.4) 526 (53.1)  469 (59.5) 467 (59.3)  
Black 362 (22.9) 334 (33.7)  215 (27.3) 223 (28.3)  
Hispanic 70 (4.4) 59 (6.0)  42 (5.3) 42 (5.3)  
Other 116 (7.3) 71 (7.2)  62 (7.9) 56 (7.1)  
Frailty, n (%)   0.043   < 0.001 
Nonfrail (VA-FI < 0.2) 760 (48.0) 463 (46.8)  389 (49.4) 389 (49.4)  
Mild frailty (VA-FI = 0.2-0.3) 439 (27.7) 294 (29.7)  226 (28.7) 226 (28.7)  
Moderate-severe frailty (VA- FI ≥ 0.3) 384 (24.3) 233 (23.5)  173 (22.0) 173 (22.0)  
Charlson Comorbidity Index, median (IQR) 1.00 (0.00-2.00) 1.00 (0.00-2.00) 0.063 1.00 (0.00-2.00) 1.00 (0.00-2.00) 0.010 
Myeloma stage, n (%)   0.361   < 0.001 
Stage I 349 (22.0) 91 (9.2)  76 (9.6) 76 (9.6)  
Stage II 350 (22.1) 241 (24.3)  200 (25.4) 200 (25.4)  
Stage III 527 (33.3) 392 (39.6)  318 (40.4) 318 (40.4)  
Missing 357 (22.6) 266 (26.9)  194 (24.6) 194 (24.6)  
Pretreatment laboratory levels, median (IQR)       
Creatinine (mg/dL) 1.20 (0.97-1.50) 1.18 (0.96-1.51) 0.013 1.15 (0.94-1.46) 1.15 (0.94-1.50) 0.026 
Calcium (mg/dL) 9.10 (8.70-9.60) 9.20 (8.80-9.70) 0.168 9.20 (8.80-9.66) 9.20 (8.80-9.70) 0.002 
Hemoglobin (g/dL) 11.10 (9.70-12.50) 10.90 (9.30-12.47) 0.066 11.01 (9.60-12.62) 11.10 (9.50-12.60) 0.010 
Platelets (n per mL) 205.00 (158.50-250.00) 202.00 (160.00-252.75) 0.014 203.00 (158.00-250.00) 204.42 (162.00-252.00) 0.029 
Rurality of veteran, n (%)   0.099   0.030 
Rural 561 (35.4) 306 (30.9)  263 (33.4) 261 (33.1)  
Urban 1012 (63.9) 679 (68.6)  518 (65.7) 522 (66.2)  
Unknown 10 (0.6) 5 (0.5)  7 (0.9) 5 (0.6)  
Polypharmacy (≥5 medications), n (%) 522 (33.0) 333 (33.6) 0.014 265 (33.6) 264 (33.5) 0.003 
Healthcare utilization in the year before treatment, median (IQR)       
Any hospital utilization 0.00 (0.00-1.00) 0.00 (0.00-1.00) 0.017 0.00 (0.00-1.00) 0.00 (0.00-1.00) 0.040 
Any ED utilization 0.00 (0.00-1.00) 1.00 (0.00-2.00) 0.309 0.00 (0.00-2.00) 1.00 (0.00-2.00) 0.060 
Any NH utilization 0.00 (0.00-0.00) 0.00 (0.00-0.00) 0.096 0.00 (0.00-0.00) 0.00 (0.00-0.00) 0.019 

ED, emergency department; IQR, interquartile range; NH, nursing home utilization (including postacute care/rehabilitation).

VRd vs Rd and survival: overall and based on frailty severity

A total of 1213 patients died during the study period, with a median follow-up time of 3.1 years. In the combined overall population, after matching within each frailty category based on the stage and PS, VRd was associated with improved survival vs Rd (log-rank P = .007; Figure 2) and lower mortality (HR, 0.81; 95% CI, 0.70-0.94;). When viewing each matched frailty category separately, patients who were moderate-severely frail demonstrated the strongest association between VRd vs Rd and lower mortality (HR, 0.74; 95% CI, 0.56-0.97; Figures 3A and 4), whereas the association was weaker in those who were mildly frail (HR, 0.80; 95% CI, 0.61-1.05; Figures 3B and 4) and nonfrail (HR, 0.86; 95% CI, 0.67-1.10; Figures 3C and 4). There was no evidence of the effect of heterogeneity across frailty categories (P = .73). In our exploratory analysis stratifying based on both frailty and stage, the largest magnitude of benefit associated with VRd vs Rd was seen in patients who were moderate-severely frail with stage III disease (HR, 0.57; 95% CI, 0.36-0.90; supplemental Figure 2).

Figure 2.

Survival of VRd vs. Rd in combined overall population. Kaplan-Meier analyses of patients with newly diagnosed MM, evaluating differences in survival between those initially treated with VRd (orange) and those initially treated with Rd (blue), after matching based on stage and PS within each frailty category and combining these for the overall population.

Figure 2.

Survival of VRd vs. Rd in combined overall population. Kaplan-Meier analyses of patients with newly diagnosed MM, evaluating differences in survival between those initially treated with VRd (orange) and those initially treated with Rd (blue), after matching based on stage and PS within each frailty category and combining these for the overall population.

Close modal
Figure 3.

Survival of VRd vs. Rd within each frailty category. Kaplan-Meier analyses of patients who were (A) moderate-severely frail, (B) mildly frail, and (C) nonfrail with newly diagnosed MM, evaluating differences in survival between those initially treated with VRd (orange) and those initially treated with Rd (blue), after matching for stage and PS within each frailty category.

Figure 3.

Survival of VRd vs. Rd within each frailty category. Kaplan-Meier analyses of patients who were (A) moderate-severely frail, (B) mildly frail, and (C) nonfrail with newly diagnosed MM, evaluating differences in survival between those initially treated with VRd (orange) and those initially treated with Rd (blue), after matching for stage and PS within each frailty category.

Close modal
Figure 4.

Mortality of VRd vs. Rd within each frailty category and in combined overall population. Forest plot of HRs and 95% CIs from Cox proportional hazards models evaluating association of VRd vs Rd with mortality, based on frailty severity and in the combined overall population after matching or stage and PS within frailty categories.

Figure 4.

Mortality of VRd vs. Rd within each frailty category and in combined overall population. Forest plot of HRs and 95% CIs from Cox proportional hazards models evaluating association of VRd vs Rd with mortality, based on frailty severity and in the combined overall population after matching or stage and PS within frailty categories.

Close modal

VRd vs Rd and hospitalizations: overall and based on frailty severity

There were 2440 unplanned VA hospitalizations in the year after treatment initiation. In the combined overall population after matching within each frailty category based on stage and PS, veterans initiated on VRd had a modestly higher incidence of hospitalizations than those initiated on Rd (VRd, 1.44 hospitalizations per 5 person-years [95% CI, 1.34-1.54]; Rd, 1.18 hospitalizations per 5 person-years [95% CI, 1.10-1.27]); incidence rate ratio [IRR], 1.22; 95% CI, 1.03-1.34; Table 2). When viewing each frailty category separately, there was evidence that the association between VRd vs Rd and increased incidence of hospitalizations was present in patients who were nonfrail (IRR, 1.24; 95% CI, 1.06-1.46) and mildly frail (IRR, 1.27; 95% CI, 1.06-1.54), but this association weakened in patients who were moderate-severely frail (IRR, 1.14; 95% CI, 0.96-1.36).

Table 2.

Incidence rate of unplanned hospitalizations in patients newly diagnosed with MM treated with VRd vs Rd, based on frailty severity and in combined overall population after matching for stage and PS within frailty categories

PopulationIncidence rate per 5 person-years (95% CI)Incidence rate ratio (95% CI)
VRdRd
Moderate-severe frailty (VA-FI ≥ 0.3) 1.36 (1.21-1.54) 1.19 (1.05-1.35) 1.14 (0.96-1.36) 
Mild frailty (VA-FI = 0.2-0.3) 1.44 (1.26-1.65) 1.13 (0.99-1.29) 1.27 (1.06-1.54) 
Nonfrail (VA-FI < 0.2) 1.51 (1.34-1.69) 1.21 (1.08-1.35) 1.24 (1.06-1.46) 
Overall 1.44 (1.34-1.54) 1.18 (1.1-1.27) 1.22 (1.10-1.34) 
PopulationIncidence rate per 5 person-years (95% CI)Incidence rate ratio (95% CI)
VRdRd
Moderate-severe frailty (VA-FI ≥ 0.3) 1.36 (1.21-1.54) 1.19 (1.05-1.35) 1.14 (0.96-1.36) 
Mild frailty (VA-FI = 0.2-0.3) 1.44 (1.26-1.65) 1.13 (0.99-1.29) 1.27 (1.06-1.54) 
Nonfrail (VA-FI < 0.2) 1.51 (1.34-1.69) 1.21 (1.08-1.35) 1.24 (1.06-1.46) 
Overall 1.44 (1.34-1.54) 1.18 (1.1-1.27) 1.22 (1.10-1.34) 

Sensitivity analyses

Results were similar in our sensitivity analyses, which showed that VRd vs Rd was associated with a reduction in mortality in the combined overall population and that patients who were moderate-severely frail demonstrated the strongest evidence of mortality reduction (supplemental Table 7).

Similar to the results of the SWOG S0777 trial, we found in US veterans newly diagnosed and treated for MM a mortality benefit associated with the more intensive triplet VRd over the less intensive doublet Rd. However, the evidence of benefit was strongest in veterans with the highest levels of frailty, contrary to our hypothesis that the benefit would diminish. Moreover, although there was a modestly higher incidence of hospitalizations associated with VRd vs Rd, this association was diminished, rather than strengthened, among veterans with the highest levels of frailty. Our findings suggest that, on average, veterans with the highest levels of frailty derive an even greater net benefit from VRd over Rd than fit veterans, challenging historical expert opinion-based recommendations to consider Rd over VRd in patients who are frail.7,8 

To date, to our knowledge, our study is one of the largest designed to evaluate benefits and harms of more intensive vs less intensive therapy in patients with high levels of frailty newly diagnosed with and treated for MM. Within the SWOG S0777 trial itself, the authors found that the overall mortality benefit conferred by VRd was maintained in subgroup analyses stratified based on the intent to undergo transplant and actual receipt of transplant.6,33 In a post hoc analysis recently presented as a conference abstract, the magnitude of mortality benefit was weaker in patients aged ≥65 years (n = 269) vs patients aged <65 years (n = 202).34 However, frailty was not measured in SWOG S0777, and only 14% of trial participants had an Eastern Cooperative Oncology Group performance status > 1, making it difficult to compare results across the fitter trial population and our frailer veteran population. Investigators have used a simplified frailty score in post hoc analyses of other RCTs to find that the benefit of more intensive MM regimens is generally preserved in frail subgroups.35-37 Similar to most RCTs, however, exclusion criteria of these trials restricted the enrollment of patients with poorer performance status and multimorbidity, limiting both the number of patients who are frail available for subgroup analyses and the generalizability of results (ie, patients who are frail as defined in the trial may not correspond to those treated in practice).35,36,38 Differences in frailty measures and the categorization of frailty severity may also hinder comparisons of outcomes between different trials as well as between trials and real-world studies.39 

Although observational studies other than ours have investigated outcomes of VRd and Rd in real-world clinical populations with newly diagnosed MM, many in descriptive analyses do not allow for valid head-to-head comparisons.40-43 A conference abstract reporting on a cohort study comparing VRd with Rd suggested that VRd showed survival benefit in patients age <75 years but not in patients age ≥75 years; however, these analyses were not adjusted for any confounders, with substantial missingness in covariates.44 Better-powered multivariable analyses of large clinical and administrative databases have compared triplet and doublet regimens, but they either lacked data on survival45 or did not control for important confounders such as MM stage and/or frailty.46,47 Our analysis leveraged several features that reduce bias in observational studies examining treatments. Emulating a specific target trial in SWOG S0777, we used a new user, active comparator design with validated, highly specific selection algorithms to reduce misclassification, selection, and immortal time biases.48,49 

We also leveraged both the structured and unstructured data within VA to control for a wide range of pretreatment confounders. Multiple pretreatment factors in our cohort could have driven the individual decisions to initiate treatment with either VRd or Rd, including MM aggressiveness, concerns regarding triplet regimens used in patients who were frail or regarding bortezomib-specific side effects in patients with multiple comorbidities, and logistical challenges/patient preferences concerning travel to oncology clinics to receive bortezomib injections (vs oral-only Rd). We controlled for frailty by stratifying based on the levels of severity and balanced these within each level across VRd vs Rd treatment groups and for many other pretreatment covariates, including MM stage, MM-specific laboratory variables (eg, creatinine), rurality, and comorbidity burden. Moreover, evolving treatment practices and greater awareness of trial data (with the main analysis of SWOG S0777 published in 2017) could have influenced treatment decisions over the time period of our study, but we did not find any association between the calendar year of treatment initiation and survival. Although residual unmeasured confounding could still remain, as with any observational study, the control for measured confounding and the other study features noted above mitigate the common threats to internal validity of observational analyses while maintaining good external validity in our frail veteran population.

Our findings call for viewing frailty not merely as a measure of treatment tolerability but more broadly as a summary measure of a patient’s health status that reflects both cancer and noncancer contributors. The electronic VA-FI we used to measure frailty is based on one of the most widely-accepted constructs of frailty in deficit accumulation,50 which carries an advantage over other frailty constructs used in MM51 in its ability to evaluate the individual deficits contributing to the overall frailty score. When looking at the frailty deficits in our veteran population, myeloma-related deficits and aggressive (stage III) MM were more prevalent in patients who were moderate-severely frail than in those who were nonfrail. Moreover, these myeloma-related deficits were more prevalent in patients who were frail in this study than in those identified through the VA-FI in a more general veteran population.15,52 MM and its complications likely drive this higher prevalence of frailty, and as frailty increases, MM and its complications contribute disproportionately more to frailty than other, nononcologic frailty deficits. Therefore, more intensive (eg, VRd) vs less intensive (eg, Rd) treatment for a veteran who is frail and having MM may actually deliver more effective treatment of frailty itself. Further supporting this interpretation is our exploratory analysis that found that the frailest patients with the most aggressive disease (stage III MM) derived the greatest mortality benefit associated with VRd compared with all other frailty/stage subgroups.

However, the added benefit of a more intensive regimen must outweigh its added risk of toxicity, which increases with frailty.2,53,54 Indeed, the recommendation to consider Rd in patients who are frail stems from the concern that the added toxicity of triplets in such patients outweighs their treatment benefit. We did observe a slightly higher incidence of unplanned hospitalizations associated with VRd thanwith Rd, consistent with the higher rate of adverse events associated with VRd in SWOG S0777,6 but this association weakened in patients who were moderate-severely frail. If it is true that MM and its complications are underlying a large portion of the frailty of a veteran, then more effective MM treatment may lower the risk of myeloma-associated hospitalizations to balance out the higher risk of treatment toxicity–associated hospitalizations. Other harms must be considered in association with VRd vs Rd, including more trips to the clinic to receive bortezomib injections that could be burdensome to patients who are frail and considered as not worth the survival benefit. Moreover, patients who are frail may not accept a survival benefit if treatment with VRd leads to severe peripheral neuropathy, especially if the neuropathy is disabling. Future prospective studies should measure peripheral neuropathy and other frailty deficits using validated and responsive measures to determine what aspects of frailty may be changing in response to VRd. Nonetheless, our findings argue that for patients with high levels of frailty who prioritize prolongation of survival, the more intensive triplet VRd is on average superior to the doublet Rd and should be recommended over the doublet, not vice versa.

There are limitations to our study. First, we did not measure the dose of lenalidomide and bortezomib used at treatment initiation, nor did we measure postprotocol regimens used after VRd or Rd. Our goal in this study was to compare the triplet vs the doublet regimens to address the question of whether patients who are frail fare less well on 3 drugs rather than 2, and our intention-to-treat analysis mirrored SWOG S0777 in primarily evaluating the effect of the initial assignment to VRd or Rd. Second, unmeasured confounding, for example, with respect to cytogenetics, could be present in our study. However, we did measure disease aggressiveness via the main MM staging systems and achieved a balance across VRd vs Rd with respect to stage and other measures of MM risk. Moreover, although we modeled our population after SWOG S0777, because patients initiating therapy were not eligible for immediate stem cell transplantation, we could not measure and, thus, balance at treatment initiation the intention of an oncologist to prescribe a transplant at some point in the future. Although, in the overall combined population, a modestly higher percentage (21.6%) of patients in the VRd group underwent transplantation at some point in our follow-up period compared with the Rd group (17%), this could have been in part due to subsequent improvement in health and disease status with more effective initial therapy rather than differences in baseline confounding. Of note, transplantation rates during the follow-up period were nearly equal across VRd and Rd in patients who were moderate-severely frail. Lastly, the competing risk of mortality could have affected our measurement of the incidence rates of unplanned hospitalizations.

In conclusion, our findings not only confirm the net benefit of VRd over Rd in US veterans newly diagnosed with MM but also suggest that this benefit is strongest in patients with the highest levels of frailty, countering historical recommendations to consider doublets in this population. Although there are newer triplet regimens that are preferred in older patients, our findings argue that frailty is not merely a measure of treatment tolerability and a patient who is frail with MM should be considered treatable cause of their frailty, wherein more intensive treatment may be more effective. Rigorous observational studies of real-world practice should be used to determine how the benefits and harms of novel regimens tested in RCTs translate to frail, underrepresented populations in practice.

This work was supported by the VA Office of Research and Development, Cooperative Studies Program (N.R.F, N.V.D., and M.B.); HCSRN-OAICs AGING Initiative Pilot (NIAR33-AG057806) (C.D. and N.R.F.); VA Career Development Award IK2CX002218 (C.D.); VA Merit Review Award 1I01BX001584 (N.C.M.); National Institutes of Health grants P01-155258-07 and P50-100707 (N.C.M.); American Heart Association 857078 (N.R.F.); and NIA K24-AG073527 (D.H.K.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The views expressed are those of the authors and do not represent the views of VA or the US Government.

Contribution: C.D., J.L., D.H.K., N.C.M., N.R.F., and J.A.D. conceptualized and designed the study; J.L., J.B., J.C., and N.R.F. performed data curation and quality control; C.D., J.L., J.B., J.C., and N.R.F. performed data analysis; and all authors contributed to the interpretation of results and manuscript writing.

Conflict-of-interest disclosure: N.C.M. reported receiving personal fees from Bristol Myers Squibb, Janssen, Amgen, Takeda, OncoPep, AbbVie, Karyopharm, Novartis, Legend, Raqia, Adaptive Biotechnolology, and Pfizer outside submitted work; and IP licensed to OncoPep and held stocks in C4 Therapeutics. The remaining authors declare no competing financial interests.

Correspondence: Clark DuMontier, Division of Aging Dana-Farber Cancer Institute, Division of Population Sciences Assistant Professor of Medicine, Geriatrician and Clinical Investigator VA Boston Healthcare System, New England GRECC Brigham and Women’s Hospital, Harvard Medical School, 150 S Huntington Ave, Boston, MA 02130; e-mail: cdumontier@bwh.harvard.edu.

1.
Grant
SJ
,
Freeman
CL
,
Rosko
AE
.
Treatment of older adult or frail patients with multiple myeloma
.
Hematology Am Soc Hematol Educ Program
.
2021
;
2021
(
1
):
46
-
54
.
2.
Palumbo
A
,
Bringhen
S
,
Mateos
MV
, et al
.
Geriatric assessment predicts survival and toxicities in elderly myeloma patients: an International Myeloma Working Group report
.
Blood
.
2015
;
125
(
13
):
2068
-
2074
.
3.
White
P
,
Cooley
A
,
Bal
S
, et al
.
Disparities in clinical trial participation among older adults with multiple myeloma in the United States
.
J Geriatric Oncol
.
2022
;
13
(
8
):
1241
-
1243
.
4.
Ludmir
EB
,
Mainwaring
W
,
Lin
TA
, et al
.
Factors associated with age disparities among Cancer Clinical Trial participants
.
JAMA Oncol
.
2019
;
5
(
12
):
1769
-
1773
.
5.
Cook
G
,
Larocca
A
,
Facon
T
,
Zweegman
S
,
Engelhardt
M
.
Defining the vulnerable patient with myeloma-a frailty position paper of the European Myeloma Network
.
Leukemia
.
2020
;
34
(
9
):
2285
-
2294
.
6.
Durie
BGM
,
Hoering
A
,
Abidi
MH
, et al
.
Bortezomib with lenalidomide and dexamethasone versus lenalidomide and dexamethasone alone in patients with newly diagnosed myeloma without intent for immediate autologous stem-cell transplant (SWOG S0777): a randomised, open-label, phase 3 trial
.
Lancet
.
2017
;
389
(
10068
):
519
-
527
.
7.
Mina
R
,
Bringhen
S
,
Wildes
TM
,
Zweegman
S
,
Rosko
AE
.
Approach to the older adult with multiple myeloma
.
Am Soc Clin Oncol Educ Book
.
2019
;
39
:
500
-
518
.
8.
Rajkumar
SV
.
Multiple myeloma: overview of management
. UpToDate. 2022. Updated. Accessed 30 December 2022. https://www.uptodate.com/contents/multiple-myeloma-overview-of-management.
9.
Facon
T
,
Kumar
S
,
Plesner
T
, et al
.
Daratumumab plus lenalidomide and dexamethasone for untreated myeloma
.
N Engl J Med
.
2019
;
380
(
22
):
2104
-
2115
.
10.
Price
LE
,
Shea
K
,
Gephart
S
.
The Veterans Affairs's Corporate Data Warehouse: uses and implications for nursing research and practice
.
Nurs Adm Q
.
2015
;
39
(
4
):
311
-
318
.
11.
Klein
M
,
Scaria
G
,
Ganti
AK
.
Utilization of the Veterans Affairs Central Cancer Registry to evaluate lung cancer outcomes
.
Semin Oncol
.
2019
;
46
(
4-5
):
321
-
326
.
12.
Zullig
LL
,
Jazowski
SA
,
Chawla
N
, et al
.
Summary of Veterans Health Administration Cancer Data Sources
.
J Registry Manag
.
2019
;
46
(
3
):
76
-
83
.
13.
La
J
,
DuMontier
C
,
Hassan
H
, et al
.
Validation of algorithms to select patients with multiple myeloma and patients initiating myeloma treatment in the National Veterans Affairs Healthcare System
.
Pharmacoepidemiol Drug Saf
.
2023
;
32
(
5
):
558
-
566
.
14.
Orkaby
AR
,
Nussbaum
L
,
Ho
YL
, et al
.
The burden of frailty among U.S. veterans and its association with mortality, 2002-2012
.
J Gerontol A Biol Sci Med Sci
.
2019
;
74
(
8
):
1257
-
1264
.
15.
Cheng
D
,
DuMontier
C
,
Yildirim
C
, et al
.
Updating and validating the U.S. Veterans Affairs Frailty Index: transitioning from ICD-9 to ICD-10
.
J Gerontol A Biol Sci Med Sci
.
2021
;
76
(
7
):
1318
-
1325
.
16.
American Geriatrics Society 2023 Annual Scientific Meeting
.
J Am Geriatr Soc
.
2023
;
71
(
S1
):
S1
-
S398
.
17.
Rockwood
K
,
Mitnitski
A
.
Frailty defined by deficit accumulation and geriatric medicine defined by frailty
.
Clin Geriatr Med
.
2011
;
27
(
1
):
17
-
26
.
18.
Searle
SD
,
Mitnitski
A
,
Gahbauer
EA
,
Gill
TM
,
Rockwood
K
.
A standard procedure for creating a frailty index
.
BMC Geriatr
.
2008
;
8
(
1
):
24
.
19.
Armstrong
JJ
,
Andrew
MK
,
Mitnitski
A
,
Launer
LJ
,
White
LR
,
Rockwood
K
.
Social vulnerability and survival across levels of frailty in the Honolulu-Asia Aging Study
.
Age Ageing
.
2015
;
44
(
4
):
709
-
712
.
20.
Blodgett
J
,
Theou
O
,
Kirkland
S
,
Andreou
P
,
Rockwood
K
.
Frailty in NHANES: comparing the frailty index and phenotype
.
Arch Gerontol Geriatr
.
2015
;
60
(
3
):
464
-
470
.
21.
Hoover
M
,
Rotermann
M
,
Sanmartin
C
,
Bernier
J
.
Validation of an index to estimate the prevalence of frailty among community-dwelling seniors
.
Health Rep
.
2013
;
24
(
9
):
10
-
17
.
22.
Pajewski
NM
,
Williamson
JD
,
Applegate
WB
, et al
.
Characterizing frailty status in the systolic blood pressure intervention trial
.
J Gerontol A Biol Sci Med Sci
.
2016
;
71
(
5
):
649
-
655
.
23.
Veterans Affairs Office of Rural Health
.
Rural veterans
. Accessed 20 January 2023. https://www.ruralhealth.va.gov/aboutus/ruralvets.asp.
24.
Kyle
RA
,
Gertz
MA
,
Witzig
TE
, et al
.
Review of 1027 patients with newly diagnosed multiple myeloma
.
Mayo Clin Proc
.
2003
;
78
(
1
):
21
-
33
.
25.
Gasparini
A
.
comorbidity: an R package for computing comorbidity scores
.
J Open Source Softw
.
2018
;
3
(
23
):
648
.
26.
Sohn
MW
,
Arnold
N
,
Maynard
C
,
Hynes
DM
.
Accuracy and completeness of mortality data in the Department of Veterans Affairs
.
Popul Health Metr
.
2006
;
4
:
2
.
27.
Vincent
BM
,
Wiitala
WL
,
Burns
JA
,
Iwashyna
TJ
,
Prescott
HC
.
Using Veterans Affairs Corporate Data Warehouse to identify 30-day hospital readmissions
.
Health Serv Outcomes Res Methodol
.
2018
;
18
(
3
):
143
-
154
.
28.
Fleiss
JL
.
The statistical basis of meta-analysis
.
Stat Methods Med Res
.
1993
;
2
(
2
):
121
-
145
.
29.
Stuart
EA
,
Lee
BK
,
Leacy
FP
.
Prognostic score-based balance measures can be a useful diagnostic for propensity score methods in comparative effectiveness research
.
J Clin Epidemiol
.
2013
;
66
(
8 suppl
):
S84
-
S90.e1
.
30.
Nguyen
TL
,
Collins
GS
,
Spence
J
, et al
.
Double-adjustment in propensity score matching analysis: choosing a threshold for considering residual imbalance
.
BMC Med Res Methodol
.
2017
;
17
(
1
):
78
.
31.
Mi
X
,
Hammill
BG
,
Curtis
LH
,
Greiner
MA
,
Setoguchi
S
.
Impact of immortal person-time and time scale in comparative effectiveness research for medical devices: a case for implantable cardioverter-defibrillators
.
J Clin Epidemiol
.
2013
;
66
(
8 suppl
):
S138
-
S144
.
32.
von Elm
E
,
Altman
DG
,
Egger
M
, et al
.
The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies
.
Lancet
.
2007
;
370
(
9596
):
1453
-
1457
.
33.
Durie
BGM
,
Hoering
A
,
Sexton
R
, et al
.
Longer term follow-up of the randomized phase III trial SWOG S0777: bortezomib, lenalidomide and dexamethasone vs. lenalidomide and dexamethasone in patients (Pts) with previously untreated multiple myeloma without an intent for immediate autologous stem cell transplant (ASCT)
.
Blood Cancer J
.
2020
;
10
(
5
):
53
.
34.
Durie
BGM
,
Lam
A
,
Ngo
M
,
Pei
H
,
Ammann
EM
.
Post-hoc analysis of efficacy and safety in the SWOG S0777 Trial stratified by age
.
Blood
.
2022
;
140
(
suppl 1
):
10009
-
10010
.
35.
Mateos
MV
,
Dimopoulos
MA
,
Cavo
M
, et al
.
Daratumumab plus bortezomib, melphalan, and prednisone versus bortezomib, melphalan, and prednisone in transplant-ineligible newly diagnosed multiple myeloma: Frailty Subgroup Analysis of ALCYONE
.
Clin Lymphoma Myeloma Leuk
.
2021
;
21
(
11
):
785
-
798
.
36.
Facon
T
,
Cook
G
,
Usmani
SZ
, et al
.
Daratumumab plus lenalidomide and dexamethasone in transplant-ineligible newly diagnosed multiple myeloma: frailty subgroup analysis of MAIA
.
Leukemia
.
2022
;
36
(
4
):
1066
-
1077
.
37.
Schjesvold
F
,
Bringhen
S
,
G Richardson
P
, et al
.
Isatuximab plus pomalidomide and dexamethasone in frail patients with relapsed/refractory multiple myeloma: ICARIA-MM subgroup analysis
.
Am J Hematol
.
2021
;
96
(
11
):
E423
-
e427
.
38.
Frieden
TR
.
Evidence for health decision making - beyond randomized, controlled trials
.
N Engl J Med
.
2017
;
377
(
5
):
465
-
475
.
39.
Mian
H
,
McCurdy
A
,
Giri
S
, et al
.
The prevalence and outcomes of frail older adults in clinical trials in multiple myeloma: a systematic review
.
Blood Cancer J
.
2023
;
13
(
1
):
6
.
40.
Zamagni
E
,
Dhanasiri
S
,
Ghale
A
,
Moore
A
,
Roussel
M
.
Real-world analysis of patient characteristics, treatment outcomes, and healthcare resource utilization across Europe in patients with newly diagnosed multiple myeloma ineligible for stem cell transplantation who received lenalidomide- or bortezomib-based regimens
.
Leuk Lymphoma
.
2021
;
62
(
10
):
2492
-
2501
.
41.
Medhekar
R
,
Ran
T
,
Fu
AZ
,
Patel
S
,
Kaila
S
.
Real-world patient characteristics and treatment outcomes among nontransplanted multiple myeloma patients who received Bortezomib in combination with Lenalidomide and Dexamethasone as first line of therapy in the United States
.
BMC Cancer
.
2022
;
22
(
1
):
901
.
42.
Joseph
NS
,
Kaufman
JL
,
Dhodapkar
MV
, et al
.
Long-term follow-up results of lenalidomide, bortezomib, and dexamethasone induction therapy and risk-adapted maintenance approach in newly diagnosed multiple myeloma
.
J Clin Oncol
.
2020
;
38
(
17
):
1928
-
1937
.
43.
Cejalvo
MJ
,
Bustamante
G
,
González
E
, et al
.
Treatment patterns and outcomes in real-world transplant-ineligible patients newly diagnosed with multiple myeloma
.
Ann Hematol
.
2021
;
100
(
7
):
1769
-
1778
.
44.
Decaux
O
,
Perrot
A
,
Frenzel
L
, et al
.
RD vs. VRD: what is the first-line, real-life management of multiple myeloma patients ineligible for stem cell auto transplantation in the Emmy cohort?
.
Blood
.
2021
;
138
(
suppl 1
):
3766
.
45.
Chari
A
,
Parikh
K
,
Ni
Q
,
Abouzaid
S
.
Treatment patterns and clinical and economic outcomes in patients with newly diagnosed multiple myeloma treated with lenalidomide- and/or bortezomib-containing regimens without stem cell transplant in a real-world setting
.
Clin Lymphoma Myeloma Leuk
.
2019
;
19
(
10
):
645
-
655
.
46.
He
J
,
Schmerold
L
,
Van Rampelbergh
R
, et al
.
Treatment pattern and outcomes in newly diagnosed multiple myeloma patients who did not receive autologous stem cell transplantation: a real-world observational study : treatment pattern and outcomes in patients with multiple myeloma
.
Adv Ther
.
2021
;
38
(
1
):
640
-
659
.
47.
Jimenez-Zepeda
VH
,
Venner
C
,
McCurdy
A
, et al
.
Real-world outcomes with bortezomib-containing regimens and lenalidomide plus dexamethasone for the treatment of transplant-ineligible multiple myeloma: a multi-institutional report from the Canadian Myeloma Research Group database
.
Br J Haematol
.
2021
;
193
(
3
):
532
-
541
.
48.
Lund
JL
,
Richardson
DB
,
Stürmer
T
.
The active comparator, new user study design in pharmacoepidemiology: historical foundations and contemporary application
.
Curr Epidemiol Rep
.
2015
;
2
(
4
):
221
-
228
.
49.
Yoshida
K
,
Solomon
DH
,
Kim
SC
.
Active-comparator design and new-user design in observational studies
.
Nat Rev Rheumatol
.
2015
;
11
(
7
):
437
-
441
.
50.
Rockwood
K
,
Mitnitski
A
.
Frailty in relation to the accumulation of deficits
.
J Gerontol A Biol Sci Med Sci
.
2007
;
62
(
7
):
722
-
727
.
51.
Giri
S
,
Williams
G
,
Rosko
A
, et al
.
Simplified frailty assessment tools: are we really capturing frailty or something else?
.
Leukemia
.
2020
;
34
(
7
):
1967
-
1969
.
52.
DuMontier
C
,
Fillmore
NR
,
Yildirim
C
, et al
.
Contemporary analysis of electronic frailty measurement in older adults with multiple myeloma treated in the National US Veterans Affairs Healthcare System
.
Cancers (Basel)
.
2021
;
13
(
12
):
3053
.
53.
Hurria
A
,
Mohile
S
,
Gajra
A
, et al
.
Validation of a prediction tool for chemotherapy toxicity in older adults with cancer
.
J Clin Oncol
.
2016
;
34
(
20
):
2366
-
2371
.
54.
Zweegman
S
,
Engelhardt
M
,
Larocca
A
;
EHA SWG on ‘Aging and Hematology’
.
Elderly patients with multiple myeloma: towards a frailty approach?
.
Curr Opin Oncol
.
2017
;
29
(
5
):
315
-
321
.

Author notes

C.D. and J.L. are joint first authors.

J.A.D. and N.R.F. are joint senior authors.

The derived data generated in this research are available, as permitted by VA policy, on reasonable request from the corresponding author, Clark DuMontier (cdumontier@bwh.harvard.edu).

The full-text version of this article contains a data supplement.