Introduction
Hematopoietic stem cell transplantation (HSCT) is a potentially curative treatment for patients with hematological malignancies, and advances in transplantation techniques have markedly increased the number of long-term survivors 1 . As survival rates continue to improve, attention in nursing and survivorship care has shifted toward long-term functional outcomes and quality of life among HSCT survivors 2 . Return to work (RTW) is widely recognized as a key indicator of comprehensive rehabilitation, reflecting not only physical recovery but also psychosocial well-being, financial independence, and societal reintegration 3 . Beyond its economic benefits, RTW represents a critical step in “ returning to normalcy", strengthening survivors' self-identity, and serving as concrete evidence of recovery 3 .
Despite these benefits, employment remains a substantial challenge across the post-transplant trajectory. Previous studies have demonstrated that only 62.4% of survivors were employed one year post-transplant 4 , and up to 23% remained unemployed even after five years 5 . Importantly, these employment challenges are not confined to adult HSCT recipients; individuals who undergo transplantation during childhood or adolescence continue to experience vocational disadvantages well into adult life 6 , 7 . These findings highlight that the gap between clinical cure and functional recovery persists long after hospital discharge, driven by multifaceted barriers including treatment-related sequelae and persistent symptom burdens.
Among late effects after HSCT, fatigue is one of the most prevalent and persistent symptoms and is consistently identified as highly burdensome 8 . Fatigue compromises both cognitive and physical functioning essential for sustained employment 9 , and has been consistently associated with delayed or unsuccessful RTW during the first post-transplant year 4 , 5 . Chronic graft-versus-host disease (cGvHD) represents another complex barrier to vocational reintegration 10 - 13 . However, most prior studies have treated cGvHD simply as a binary variable (present vs. absent) 11 , 12 , overlooking the heterogeneity of its organ-specific manifestations. Crucially, binary classifications fail to capture the distinct functional impairments imposed by different organs. This limitation underscores the need for research that evaluates individual organ involvement to guide more precise, symptom-tailored survivorship interventions.
Existing evidence on treatment-related predictors of RTW remains inconsistent, making it difficult to accurately stratify risk of RTW outcomes and guide follow-up care. With respect to transplant modality, Morrison et al. (2016) did not identify HSCT type as a significant determinant of employment outcomes 4 , whereas Kurosawa et al. (2023) reported that transplant type significantly predicted RTW 12 . Similarly, Winterling et al. (2018) observed no significant association between total body irradiation (TBI) exposure and RTW 7 , while Bhatt et al. (2021) observed poorer employment outcomes among recipients of myeloablative conditioning that included TBI 14 . Such inconsistencies likely reflect heterogeneity in study populations and methodologies. This emphasizes the critical need for longitudinal approaches that do not view treatment factors in isolation but integrate them with dynamic, symptom-level predictors to better understand the complex survivorship trajectory.
Most existing research is constrained by cross-sectional designs that provide only a static snapshot of employment status. Such designs fail to capture the dynamic trajectory of RTW, overlooking critical patterns such as delayed workforce entry, recurrent sick leave, or later involuntary withdrawal from employment 5 , 7 , 12 . Moreover, few studies have evaluated the effects of organ-specific cGvHD 10 , leaving gaps in understanding how distinct symptom presentation differentially impede vocational functioning. To address these gaps, this study employed a retrospective longitudinal design to quantify the probability and timing of RTW, thereby accounting for variations in follow-up time and reflecting the dynamic nature of work resumption. We further evaluated the impact of symptom burden—specifically fatigue and organ-specific cGvHD—alongside treatment factors, to identify modifiable predictors of vocational recovery.
Materials and Methods
Study design and setting
This retrospective longitudinal cohort study investigated factors associated with RTW among adult HSCT recipients. The study was conducted at a tertiary comprehensive cancer center in Taiwan, which performs an average of 2–3 HSCT procedures per month. Data were collected from patients treated between January 1, 2016, and December 31, 2021. To allow for an adequate evaluation of long-term vocational outcomes and survivorship trajectories, patients were followed for a minimum of 12 months post-transplant. The study design and reporting adhered to the Strengthening the Reporting of Observational Studies in Epidemiology guidelines.
Participants
Eligible criteria included: (1) age ≥ 18 years; (2) receipt of HSCT during the study period; (3) history of paid employment prior to cancer diagnosis; and (4) survival with ≥ 12 months of follow-up after hospital discharge. Patients with incomplete demographic or clinical data were excluded. All consecutive eligible HSCT recipients were included to minimize selection bias and enhance representativeness. The final analytical cohort consisted of 2,334 longitudinal nursing assessment records derived from 80 HSCT recipients. During the follow-up period, 8 patients (10.0%) experienced disease relapse, and 6 patients (7.5%) died from treatment-related or disease-related complications; however, all longitudinal assessments prior to these events were included in the time-dependent analysis to maximize data utility. This extensive dataset allowed for a granular, time-dependent analysis of symptom trajectories and vocational outcomes.
Data source and variables
Data were extracted from institutional medical records and transplant follow-up databases, including routinely documented clinical assessments performed as part of survivorship care. Baseline variables were collected at the time of transplantation and included sociodemographic characteristics (age, sex, and educational attainment) and clinical variables, specifically primary diagnosis, conditioning regimen (stratified by anti-thymocyte globulin [ATG] or TBI), and HSCT type (autologous vs. allogeneic).
Time-varying variables were assessed at each follow-up visit to capture changes in clinical status over time. These included: (1) time since HSCT, measured in days post-transplant; (2) fatigue, which was assessed as part of routine clinical follow-up and reflects patients' functional status in survivorship care. Fatigue intensity was evaluated using a standard 0–10 Numeric Rating Scale; (3) organ-specific cGvHD manifestations, which were analyzed to capture the functional heterogeneity of cGvHD and its potential impact on vocational outcomes. Organ involvement was assessed for the skin, oral cavity, eyes, gastrointestinal tract, lungs, and joints/fascia; and (4) hepatic cGvHD, which was evaluated based on established clinical thresholds of serum bilirubin, aspartate aminotransferase (AST), and alanine aminotransferase (ALT) levels; (5) disease relapse status, which was monitored as a time-dependent covariate to adjust for its impact on vocational reintegration; and (6) a history of acute graft-versus-host disease (aGvHD), identified through longitudinal nursing records as documented episodes of aGvHD involving the skin, gastrointestinal tract, or liver.
RTW status was determined through routine nursing inquiries regarding employment status and was documented during follow-up visits as part of survivorship care.
Outcome measures
The primary outcome was time to RTW, defined as the interval from the date of HSCT to the first resumption of any form of paid employment. The probability and timing of RTW were estimated at key milestones (6, 12, and 24 months) using the Kaplan-Meier method. To account for varying follow-up durations, the incidence rate was calculated and expressed as the number of RTW events per 100 person-years of observation.
Statistical analyses
Descriptive statistics were used to summarize sociodemographic and clinical characteristics. As the Shapiro-Wilk test confirmed the non-normal distribution of the data, continuous variables were reported as medians with interquartile ranges (IQR), while categorical variables were presented as frequencies and percentages.
The probability of RTW over time was estimated using the Kaplan–Meier method, and the log-rank test was used to compare RTW probabilities between groups. To provide a clearer understanding of the data density, risk tables showing the number of patients at risk were included with the Kaplan–Meier curves.
Multivariable time-dependent Cox proportional hazards regression models were constructed to identify predictors of RTW. To minimize bias and adjust for potential confounders, variables such as age, sex, conditioning regimen, and disease relapse status were included in the models. To leverage the longitudinal nature of the dataset (2,334 assessments), time-varying variables—including fatigue and cGvHD status—were incorporated to capture their dynamic impact on vocational outcomes. For analysis, fatigue level was dichotomized according to the NCCN guideline definition, with scores of 4 or higher (≥ 4) considered high fatigue 15 .
Two separate models were developed: Model 1 (Entire Cohort) included all HSCT recipients to assess the impact of sociodemographic factors, clinical characteristics, and fatigue levels. Model 2 (Allogeneic Subgroup) was restricted to allogeneic HSCT recipients to evaluate the impact of organ-specific cGvHD involvement, while adjusting for other sociodemographic and clinical variables, including disease relapse status and history of aGvHD. For both models, the proportional hazards assumption was verified using Schoenfeld residuals to ensure the validity of the Cox regression.
Results were expressed as hazard ratios (HRs) with 95% confidence intervals (CIs) derived from the multivariable models. A two-tailed p < 0.05 was considered statistically significant. All analyses were conducted using IBM SPSS Statistics, version 26.0 (IBM Corp., Armonk, NY, USA).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of Koo Foundation Sun Yat-Sen Cancer Center (Approval No. 20230518A). The requirement for informed consent was waived because of the retrospective nature of the study and the use of de-identified data.
Results
Participant characteristics and data structure
The analytic cohort comprised 80 HSCT recipients. After excluding ineligible patients ( Figure 1 ), the final dataset comprised 2,334 longitudinal nursing assessment records. The median follow-up duration was 507 days (IQR 224–797), ranging from 8 to 1,505 days. During the study period, 8 patients (10.0%) experienced disease relapse, and 6 patients (7.5%) died.
The median age of participants was 53.5 years (IQR 46.8–58.2). The majority were male (62.5%) and had a college or higher level of education (50.0%). Clinically, half of the cohort was diagnosed with lymphoma (50.0%), followed by leukemia (26.3%). Slightly more than half (52.5%) underwent autologous HSCT, while 48.5% received allogeneic HSCT. Conditioning regimens containing ATG and TBI were used in 7.5% and 8.8% of patients, respectively ( Table 1 ).

Table 1.
Patient demographics and clinical characteristics (n=80)
Auto, autologous transplantation; Allo, allogeneic transplantation; ATG, anti-human thymocyte immunoglobulin; HSCT, hematopoietic stem cell transplantation; IQR, interquartile range; MM, multiple myeloma; SD standard deviation; TBI, total body irradiation
| Item | Median (IQR) /n | SD/% |
| Age | 53.5 (46.8–58.2) | 9.9 |
| Sex | ||
| Female | 30 | 37.5% |
| Male | 50 | 62.5% |
| Education level (n=78) | ||
| Elementary and high school | 15 | 19.2% |
| Senior high school | 24 | 30.8% |
| Graduate degree and above | 39 | 50.0% |
| Diagnosis | ||
| Leukemia | 21 | 26.3% |
| Lymphoma | 40 | 50.0% |
| MM | 11 | 13.8% |
| Solid Tumor & Others | 8 | 10.0% |
| Conditioning Regimen Containing ATG | ||
| Yes | 6 | 7.5% |
| No | 74 | 92.5% |
| TBI | ||
| Yes | 7 | 8.8% |
| No | 73 | 91.3% |
| Type of HSCT | ||
| Auto | 42 | 52.5% |
| Allo | 38 | 48.5% |
| Relapse | 8 | 10% |
| Mortality | 6 | 7.5% |
Incidence and incidence rate of return to work
During the follow-up period, 44 survivors (55.0%) successfully returned to work. The overall RTW incidence rate was 37.8 events per 100 person-years. The median time to first RTW was 312.5 days (IQR 144.0–567.2, range: 8–1,224 days). Based on Kaplan-Meier estimates, the probability of RTW for the entire cohort was 17.8% at 6 months, 34.1% at 12 months, and 53.3% at 24 months post-transplant, reaching an estimated 74.3% by the end of the follow-up period ( Table 2 and Figure 2 ).
Table 2.
Kaplan–Meier estimated probability of return to work among HSCT patients (n=80)
Auto, autologous transplantation; Allo, allogeneic transplantation; HSCT, hematopoietic stem cell transplantation; IQR, interquartile range; RTW, return to work *Estimated probabilities were derived from Kaplan–Meier survival analysis to account for censored observations and do not represent raw proportions (n/N).
| Time-to-event summary | ||
| Outcome Metric | Result | |
| Time to RTW (n=44 events) | ||
| Median (IQR), days | 312.5 (144.0–567.2) | |
| Range, days | 8–1,224 | |
| Kaplan–Meier estimates | ||
| Outcome Metric | Events (n) | Estimated probability* (%) |
| Total Cohort (n=80) | ||
| At 6 Months | 14 | 17.8 |
| At 1 Year | 26 | 34.1 |
| At 2 Years | 38 | 53.3 |
| At the end of follow-up | 44 | 74.3 |
| Allo (n=38) | ||
| At 6 Months | 4 | 10.8 |
| At 1 Year | 9 | 26.0 |
| At 2 Years | 17 | 52.8 |
| At the end of follow-up | 20 | 75.7 |
| Auto (n=42) | ||
| At 6 Months | 10 | 24.3 |
| At 1 Year | 17 | 41.4 |
| At 2 Years | 21 | 54.2 |
| At the end of follow-up | 24 | 68.0 |

Stratified by transplant type, autologous HSCT recipients had a relatively higher probability at 24 months (54.2% vs. 52.8%) and a faster rate of return (41.73 vs. 34.05 per 100 person-years) than allogeneic recipients ( Table 2 ), although the difference in the overall trajectory did not reach statistical significance (log-rank p =0.49). The Kaplan-Meier curves with the associated risk table illustrate the RTW trajectories between the two groups ( Figure 3 ), where autologous recipients appeared to achieve earlier vocational reintegration during the first year, though this trend was not sustained over the longer follow-up period.

Predictors of RTW in the entire cohort (Model 1)
After adjusting for sociodemographic and clinical characteristics, allogeneic HSCT was associated with a significantly lower likelihood of RTW (HR=0.17, 95% CI: 0.10–0.28, p < 0.001). Higher fatigue intensity (score ≥ 4) was also significantly associated with a lower likelihood of RTW (HR=0.67, 95% CI: 0.47–0.96, p =0.029).
Among the control variables, male sex (HR=3.87, p < 0.001), younger age (HR=0.95, p < 0.001), and a diagnosis of leukemia (HR=1.48, p =0.046) were significantly associated with a higher likelihood of RTW. Conversely, lymphoma (HR=0.49, p =0.010) and lower educational attainment (HR=0.57, p =0.022) were associated with delayed vocational recovery ( Table 3 ).
Table 3.
Hazard ratios for RTW among HSCT patients from multivariate Cox regression analysis based on 2,334 longitudinal assessments from 80 patients
Auto, autologous transplantation; Allo, allogeneic transplantation; ATG, anti-human thymocyte immunoglobulin; β, regression coefficient; CI, confidence interval; DF, degree of freedom; HSCT, hematopoietic stem cell transplantation; HR, hazard ratio; MM, Multiple myeloma; RTW, return to work; SE, standard error; TBI, total body irradiation; ref., reference Significance Levels: * p < 0.05; ** p < 0.001
| Items | β | SE | DF | p -value | HR | 95% CI of HR | |
| lower | upper | ||||||
| Age | -0.05 | 0.01 | 1 | < 0.001** | 0.95 | 0.93 | 0.97 |
| Sex | |||||||
| Female | ref. | ||||||
| Male | 1.35 | 0.17 | 1 | < 0.001** | 3.87 | 2.80 | 5.35 |
| Education Level | |||||||
| Elementary and high school | -0.56 | 0.24 | 1 | 0.022* | 0.57 | 0.36 | 0.92 |
| Senior high school | 0.03 | 0.17 | 1 | 0.840 | 1.03 | 0.75 | 1.43 |
| Graduate degree and above | ref. | ||||||
| Diagnosis | |||||||
| Leukemia | 0.39 | 0.20 | 1 | 0.046* | 1.48 | 1.01 | 2.17 |
| Lymphoma | -0.72 | 0.28 | 1 | 0.010* | 0.49 | 0.28 | 0.84 |
| MM | -0.68 | 0.36 | 1 | 0.065 | 0.51 | 0.25 | 1.04 |
| Solid Tumor & Others | ref. | ||||||
| Conditioning Regimen containing ATG | 0.27 | 0.20 | 1 | 0.172 | 1.31 | 0.89 | 1.94 |
| TBI | -0.27 | 0.23 | 1 | 0.242 | 0.76 | 0.49 | 1.20 |
| Type of HSCT | |||||||
| Auto | ref. | ||||||
| Allo | -1.78 | 0.27 | 1 | < 0.001** | 0.17 | 0.10 | 0.28 |
| Fatigue Level | |||||||
| Low-level (0-3) group | ref. | ||||||
| High-level (≥ 4) group | -0.40 | 0.18 | 1 | 0.029* | 0.67 | 0.47 | 0.96 |
| Relapse | -0.07 | 0.54 | 1 | 0.892 | 0.93 | 0.32 | 2.67 |
Impact of cGvHD and fatigue in allogeneic recipients (Model 2)
Among allogeneic HSCT recipients, organ-specific manifestations of cGvHD demonstrated differential associations with RTW. Involvement of joints/fascia (HR=0.39, 95% CI: 0.21–0.73, p =0.003) was significantly associated with delayed RTW. Although lung involvement showed a similar trend toward delayed RTW, it did not reach statistical significance after multivariable adjustment in this model (HR=0.44, 95% CI: 0.17–1.12, p =0.084). In contrast, oral (HR=1.67, 95% CI: 1.13–2.45, p =0.010) involvement was associated with a higher likelihood of RTW. After inclusion of organ-specific cGvHD variables, fatigue was no longer significantly associated with RTW ( p =0.287). Male sex (HR=3.40, p < 0.001) and younger age (HR=0.95, p < 0.001) remained significant predictors, while disease relapse (HR=0.34, p =0.100) and history of aGVHD (HR=0.74, p =0.460) did not show significant associations in this specific subgroup ( Table 4 ).
Table 4.
Hazard ratios for RTW among allo-HSCT patients with cGvHD and fatigue symptoms from multivariate Cox regression analysis based on 1,097 longitudinal assessments from 38 patients
aGvHD, acute graft-versus-host disease; allo-HSCT, allogeneic hematopoietic stem cell transplantation; β, regression coefficient; cGvHD, chronic graft-versus-host disease; CI, confidence interval; DF, degree of freedom; HR, hazard ratio; ref., reference; RTW, return to work; SE, standard error; GI, gastrointestinal Significance Levels: * p < 0.05; ** p < 0.001
| Items | β | SE | DF | p -value | HR | 95% CI of HR | |
| lower | upper | ||||||
| Age | -0.05 | 0.02 | 1 | < 0.001** | 0.95 | 0.92 | 0.98 |
| Sex | |||||||
| Female | ref. | ||||||
| Male | 1.22 | 0.27 | 1 | < 0.001** | 3.40 | 1.91 | 6.04 |
| Education Level | |||||||
| Elementary and high school | -0.68 | 0.65 | 1 | 0.300 | 0.51 | 0.14 | 1.82 |
| Senior high school | 0.11 | 0.27 | 1 | 0.680 | 1.12 | 0.66 | 1.89 |
| Graduate degree and above | ref. | ||||||
| Fatigue Level | |||||||
| Low-level (0-3) group | ref. | ||||||
| High-level (≥ 4) group | 0.31 | 0.29 | 1 | 0.287 | 1.37 | 0.77 | 2.44 |
| Organ-specific cGvHD involvement | |||||||
| Skin | 0.19 | 0.21 | 1 | 0.845 | 1.21 | 0.81 | 1.81 |
| Oral | 0.51 | 0.20 | 1 | 0.010* | 1.67 | 1.13 | 2.45 |
| Eyes | -0.20 | 0.27 | 1 | 0.489 | 0.83 | 0.48 | 1.42 |
| GI | 0.53 | 0.36 | 1 | 0.145 | 1.70 | 0.83 | 3.47 |
| Liver | 0.50 | 0.31 | 1 | 0.106 | 1.65 | 0.90 | 3.03 |
| Lung | -0.82 | 0.47 | 1 | 0.084 | 0.44 | 0.17 | 1.12 |
| Joint and Fascia | -0.94 | 0.32 | 1 | 0.003* | 0.39 | 0.21 | 0.73 |
| History of aGvHD | -0.31 | 0.42 | 1 | 0.460 | 0.74 | 0.33 | 1.66 |
| Disease relapse | -1.08 | 0.66 | 1 | 0.100 | 0.34 | 0.09 | 1.22 |
Discussion
This study provides longitudinal evidence that return to work after HSCT is not solely determined by treatment intensity or transplant type, but is largely influenced by functional symptom burden, particularly joint/fascia manifestations of cGvHD. While pulmonary involvement also showed a clinical trend toward delayed RTW, its statistical impact was less pronounced in our cohort after multivariable adjustment. By examining organ-specific cGvHD involvement over time, our findings extend prior work that treated cGvHD as a binary exposure and highlight the importance of organ-specific assessment in survivorship care.
In our cohort, 55% of HSCT survivors returned to work, with a median time to first RTW of 312.5 days (IQR 144.0–567.2) and an incidence rate of 37.8 events per 100 person-years. These findings reflect the gradual and heterogeneous recovery trajectory. Kaplan–Meier estimates indicated that the probability of RTW reached 74.3% by the end of follow-up, indicating that while initial reintegration is slow, a substantial proportion of survivors eventually resume employment over an extended period. This recovery timeline appears slower compared to U.S. cohorts, where 60–70% of survivors typically resume employment within one to two years 4 , 14 . Taken together, recovery in our cohort appears to proceed more slowly, potentially reflecting differences in transplant composition, symptom burden, and work-related cultural contexts.
Furthermore, when compared with other cancer populations in Taiwan, RTW among HSCT survivors appears to be delayed. Our findings show that only 34.1% of HSCT survivors returned to work within one year post-transplant, although the cumulative probability of RTW increased to 74.3% during extended follow-up. Previous studies have reported that approximately 70% of general cancer survivors returned to work within one year after diagnosis or treatment 16 , whereas 82.3% of cervical cancer survivors resumed employment within five years after diagnosis 17 . Collectively, these findings suggest that vocational recovery after HSCT is often prolonged and may require a longer recovery period than that observed in other cancer populations. The persistent symptom burden following transplantation may contribute to delayed work reintegration among HSCT survivors.
Regarding clinical factors, allogeneic transplantation emerged as a strong negative predictor of RTW. While Morrison et al. (2016) did not identify HSCT type as a significant determinant in a cross-sectional analysis 4 , our longitudinal findings highlight the time-dependent nature of this vocational gap. Clinically, work recovery is notably slower among allogeneic survivors—showing a slower rate of reintegration compared with autologous recipients, particularly during the first two years post-transplant. Notably, while disease relapse did not significantly impact the entire cohort, its influence appeared more pronounced among allogeneic recipients (HR=0.34, p =0.100), suggesting that the threat of relapse may represent a more direct barrier to work resumption in this subgroup, even if statistical significance was not reached due to the limited sample size. This prolonged trajectory is attributed to sustained immunosuppression, psychosocial distress, and, most critically, the burden of cGvHD 18 , necessitating a deeper examination of how specific symptoms impede recovery.
High levels of fatigue (score ≥ 4) were identified as a significant negative predictor of RTW in our entire cohort (Model 1), consistent with previous findings that fatigue is a significant barrier to RTW following HSCT 4 . However, a key finding of this study is that when organ-specific cGvHD symptoms were included in the allogeneic subgroup model, the independent effect of fatigue was no longer significant ( p =0.287). This likely reflects overlapping pathophysiology, as cGvHD often manifests with musculoskeletal pain, fascial stiffness, or pulmonary involvement—symptoms that exacerbate fatigue and functional limitation 19 , 20 . Our results suggest that specific physical impairments, particularly joints and fascia, may be the primary driver of physical impairment, thereby “ masking" the effect of general fatigue. While fatigue is a pervasive subjective experience, the structural and mechanical limitations imposed by fascial and joint cGvHD represent more direct barriers to the physical demands of employment. This emphasizes that managing the underlying cGvHD is essential to alleviating the symptom burden that hinders work reintegration.
We found that cGvHD symptoms involving the joints and fascia were significantly associated with lower RTW rates, while lung involvement showed a similar clinical trend. This aligns with Kurosawa et al. (2023), who identified respiratory and musculoskeletal issues as key determinants of sick leave and prolonged unemployment 12 . Possible explanations include that these objective functional limitations may be underreported or misinterpreted by patients as general fatigue. This further explains why organ-specific symptoms “ masked" the effect of fatigue in our multivariable model.
In contrast, oral manifestations of cGvHD in our cohort were associated with a higher likelihood of RTW. This finding differs from prior research, such as Eriksson et al. (2023), which linked overall symptom burden—particularly eating difficulties and appetite loss—to full-time sick leave. One plausible explanation is that oral symptoms, being among the earliest and most visible manifestations of cGvHD, are more likely to be promptly detected and proactively managed. Furthermore, oral involvements, unless severe, may pose fewer direct mechanical barriers to physical labor compared to musculoskeletal constraints. As a result, timely recognition and intervention, coupled with a lower impact on physical mobility, may help preserve the capacity necessary for sustained employment.
Strengths and limitations
One major strength of this study is the availability of extensive longitudinal follow-up data, which enabled a comprehensive analysis of RTW trajectories among HSCT survivors. The longitudinal design allowed for detailed evaluation of both the timing and dynamics of RTW, providing a more thorough understanding than what cross-sectional studies can offer. Another key strength is the use of standardized nurse-assessed symptom records. Unlike studies relying solely on patient self-reports, our approach minimizes recall bias and provides a rigorous clinical basis for evaluating functional barriers. Importantly, this study is among the first to specifically examine the impact of organ-specific cGvHD symptoms on RTW, thereby laying the groundwork for more targeted survivorship interventions. Additionally, the use of incidence rates and Kaplan-Meier estimates provided a standardized approach to capture RTW events across variable follow-up durations, offering a dynamic perspective that enhances comparability across clinical settings.
However, several limitations should be acknowledged. First, the relatively small sample size of the allogeneic subgroup (n=38) may have limited the statistical power to detect significant associations for certain variables, such as pulmonary involvement ( p =0.084), which showed a clinical trend toward delayed RTW but did not reach statistical significance. Second, the absence of National Institutes of Health (NIH)–defined cGvHD severity grading precluded the evaluation of global severity scores combining organ involvement numbers. Third, the lack of multidimensional fatigue assessment tools and psychological measures—such as depression, anxiety, or coping capacity—may have limited our ability to fully characterize their influence on work outcomes. Furthermore, the absence of detailed information on pre-HSCT employment characteristics, such as job intensity (full-time vs. part-time) and physical demands (manual labor vs. desk jobs), restricted our ability to analyze how different occupational requirements moderate the impact of cGvHD symptoms. In addition, although educational attainment was included, other socioeconomic factors, such as employer support, income, and family support, were not available for analysis. Collectively, these limitations suggest that RTW is influenced by complex, multidimensional factors—including disease severity, symptom burden, psychological status, and occupational context—that were not fully captured in this study.
Future research should adopt prospective, multicenter designs to incorporate more comprehensive assessments of cGvHD severity, psychological status, and workplace environmental factors, thereby enabling a more precise evaluation of their influence on employment outcomes. Such efforts will be critical to enhancing the generalizability of these findings and informing the development of tailored survivorship interventions to promote vocational recovery after HSCT.
Implications for practice
Based on these findings, a proactive and targeted clinical approach may support vocational recovery after HSCT. Given the substantial functional impact of musculoskeletal involvement, survivorship care may benefit from rigorous surveillance informed by established clinical frameworks. Specifically, early physical therapy referrals and regular joint range-of-motion assessments, alongside pulmonary function testing 21 , can facilitate early detection and timely management, thereby helping to preserve work-related physical capacity.
Furthermore, systematic organ-specific functional assessment is crucial, as recommended by the NIH Consensus Development Project on cGvHD 22 . Since nurse-led clinical assessments can identify subtle fascial or joint restrictions before they result in profound disability, integrating these evaluations into routine nursing workflows is essential. Emerging digital health tools, such as app-based monitoring systems, can enhance real-time observation of graft-versus-host manifestations and facilitate earlier intervention 23 .
Finally, integrated, multidisciplinary symptom management programs—including intensive physical and occupational therapy, pulmonary rehabilitation, oral care, and fatigue management—may play an important role in promoting functional recovery and facilitating vocational reintegration within survivorship care pathways 18 .
Conclusion
In conclusion, this longitudinal study demonstrated that more than half of HSCT survivors returned to work, with a median time to RTW of 312.5 days and Kaplan–Meier–estimated probability of 74.3% at extended follow-up. These reflect a gradual and heterogeneous recovery trajectory after transplantation, providing a clinically relevant benchmark for evaluating interventions. Allogeneic transplantation emerged as a strong negative predictor of RTW, highlighting the importance of early vocational counseling and individualized survivorship planning.
Moreover, organ-specific manifestations of cGvHD exerted distinct effects on employment outcomes: involvement of the joints and fascia significantly hindered recovery, whereas pulmonary involvement showed a similar restrictive trend. These results underscore the need for systematic, nurse-led clinical assessments to identify specific functional barriers and implement targeted rehabilitation strategies to facilitate successful work reintegration.
Conflicts of Interest
The authors declare no conflict of interest. Disclosure forms provided by the authors are available on the website.
Author Contributions
MKC designed the study, analyzed and interpreted the data, conducted the statistical analyses, and drafted the manuscript. YMH contributed to data collection, result interpretation, and manuscript revision. AMA contributed to data management, statistical analysis, methodology refinement, and interpretation of results. TDT and LWC contributed to data collection and manuscript revision. All authors reviewed and approved the final manuscript.
Disclosure of AI Use
Artificial intelligence tools were used exclusively for language polishing and sentence-level editing. The study design, data analysis, interpretation, and conclusions were developed entirely by the authors.
The authors thank the HSCT Functional Group at Koo Foundation Sun Yat-Sen Cancer Center for their support. Special appreciation is extended to Professor Shang-Liang Wu for data management and to Dr. Deborah Ng, Patrick Ng, and Professor Lih-Mih Chen for their critical reviews of the manuscript.
1. Passweg JR , Baldomero H , Bader P , Bonini C , Cesaro S , Dreger P , et al . Hematopoietic stem cell transplantation in Europe 2014: more than 40000 transplants annually . Bone Marrow Transplant . 2016 ; 51 : 786 - 92 . 26901709 10.1038/bmt.2016.20 PMC4895175
2. Gifford G , Gilroy N , Dyer G , Brice L , Kabir M , Greenwood M , et al . The experience of survival following allogeneic haematopoietic stem cell transplantation in New South Wales, Australia . Bone Marrow Transplant . 2016 ; 51 : 1361 - 8 . 27214081 10.1038/bmt.2016.135
3. Persoon S , Buffart LM , Chinapaw MJM , Nollet F , Frings-Dresen MH , Koning S , et al . Return to work experiences of patients treated with stem cell transplantation for a hematologic malignancy . Support Care Cancer . 2019 ; 27 : 2987 - 97 . 30591967 10.1007/s00520-018-4596-0 PMC6597586
4. Morrison EJ , Ehlers SL , Bronars CA , Patten CA , Brockman TA , Cerhan JR , et al . Employment Status as an Indicator of Recovery and Function One Year after Hematopoietic Stem Cell Transplantation . Biol Blood Marrow Transplant . 2016 ; 22 : 1690 - 5 . 27220264 10.1016/j.bbmt.2016.05.013
5. Tichelli A , Gerull S , Holbro A , Buser A , Nair G , Medingr M , et al . Inability to work and need for disability pension among long-term survivors of hematopoietic stem cell transplantation . Bone Marrow Transplant. . 2017 ; 52 : 1436 - 42 . 28650451 10.1038/bmt.2017.115
6. Bhatt NS , Brazauskas R , Tecca HR , Carreras J , Burns LJ , Phelan R , et al . Post-transplantation employment status of adult survivors of childhood allogeneic hematopoietic cell transplant: A report from the Center for International Blood and Marrow Transplant Research (CIBMTR) . Cancer . 2019 ; 125 : 144 - 52 . 30368771 10.1002/cncr.31781 PMC6310211
7. Winterling J , Johansson E , Wettergren L , Ljungman P , Alexanderson K . Occupational status among adult survivors following allogeneic stem cell transplantation in childhood . Eur J Cancer Care (Engl) . 2018 ; 27 : e12808 . 29315926 10.1111/ecc.12808
8. Ullrich CK , Baker KK , Carpenter PA , Flowers ME , Gooley T , Stevens S , et al . Fatigue in Hematopoietic Cell Transplantation Survivors: Correlates, Care Team Communication, and Patient-Identified Mitigation Strategies . Transplant Cell Ther . 2023 ; 29 : 200.e1 - 8 . 10.1016/j.jtct.2022.11.030 36494015
9. Alfano CM , Kent EE , Padgett LS , Grimes M , de Moor JS . Making Cancer Rehabilitation Services Work for Cancer Patients: Recommendations for Research and Practice to Improve Employment Outcomes . PM R . 2017 ; 9 : S398 - 406 . 28942911 10.1016/j.pmrj.2017.06.019 PMC5657535
10. Eriksson LV , Holmberg K , Lundh Hagelin C, Wengström Y , Bergkvist K , Winterling J . Symptom Burden and Recovery in the First Year After Allogeneic Hematopoietic Stem Cell Transplantation . Cancer Nurs . 2023 ; 46 : 77 - 85 . 35283470 10.1097/NCC.0000000000001077
11. Eriksson L , Wennman-Larsen A , Bergkvist K , Ljungman P , Winterling J . Important factors associated with sick leave after allogeneic haematopoietic stem cell transplantation-a 1-year prospective study . J Cancer Surviv . 2021 ; 15 : 933 - 41 . 33417181 10.1007/s11764-020-00986-5 PMC7790935
12. Kurosawa S , Yamaguchi T , Mori A , Matsuura T , Masuko M , Murata M , et al . Incidence and predictors of recurrent sick leave in survivors who returned to work after allogeneic hematopoietic cell transplantation . J Cancer Surviv . 2023 ; 17 : 781 - 94 . 36048313 10.1007/s11764-022-01250-8
13. Yu J , Bhatt V , Galvin J , Jariwala-Parikh K , Thiel E . Extended disability leave and related costs among employed patients with versus without graft-versus-host disease following hematopoietic stem cell transplantation . Support Care Cancer . 2025 ; 33 : 543 . 40471352 10.1007/s00520-025-09561-z PMC12141405
14. Bhatt NS , Brazauskas R , Salit RB , Syrijala K , Bo-Subait S , Tecca H , et al . Return to Work Among Young Adult Survivors of Allogeneic Hematopoietic Cell Transplantation in the United States . Transplant Cell Ther . 2021 ; 27 : 679.e1 - 8 . 10.1016/j.jtct.2021.04.013 PMC8425287 33895402
15. National Comprehensive Cancer Network . NCCN Clinical Practice Guidelines in Oncology: Cancer-Related Fatigue . Version 2. 2025 . https://www.nccn.org/professionals/physician_gls/pdf/fatigue.pdf [Accessed: 30 January 2026]
16. Chen WL , Chen YY , Wu WT , Ho CL , Wang CC . Life expectancy estimations and determinants of return to work among cancer survivors over a 7-year period . Sci Rep . 2021 ; 11 : 12858 . 34145368 10.1038/s41598-021-92306-9 PMC8213801
17. Sun YS , Chen WL , Wu WT , Wang CC . The Fact of Return to Work in Cervical Cancer Survivors and the Impact of Survival Rate: An 11-Year Follow-Up Study . Int J Environ Res Public Health . 2021 ; 18 : 10703 . 34682449 10.3390/ijerph182010703 PMC8535606
18. Salit RB , Lee SJ , Burns LJ , Shaw BE , Majhail NS , Bhatt NS , et al . Return-to-Work Guidelines and Programs for Post-Hematopoietic Cell Transplantation Survivors: An Initial Survey . Biol Blood Marrow Transplant . 2020 ; 26 : 1520 - 6 . 32360563 10.1016/j.bbmt.2020.04.022 PMC7546432
19. Perić Z , Basak G , Koenecke C , Moiseev I , Chauhan J , Asaithambi S , et al . Understanding the Needs and Lived Experiences of Patients With Graft-Versus-Host Disease: Real-World European Public Social Media Listening Study . JMIR Cancer . 2023 ; 9 : e42905 . 37948101 10.2196/42905 PMC10674148
20. Jagasia MH , Greinix HT , Arora M , Williams KM , Wolff D , Cowen EW , et al . National Institutes of Health Consensus Development Project on Criteria for Clinical Trials in Chronic Graft-Versus-Host Disease: I. The 2014 Diagnosis and Staging Working Group Report . Biol Blood Marrow Transplant . 2015 ; 21 : 389 - 401.e1 . 25529383 10.1016/j.bbmt.2014.12.001 PMC4329079
21. Tsai XC , Yao ZH , Chou SW , Liu YC , Chen CH , Yeh CH , et al . Taiwan consensus on pulmonary chronic graft-versus-host disease: A joint statement from the Taiwan Society of Blood and Marrow Transplantation (TBMT) and Taiwan Society of Pulmonary and Critical Care Medicine (TSPCCM) . Br J Haematol . 2025 ; 207 : 1815 - 31 . 40926299 10.1111/bjh.70074 PMC12624199
22. Kitko CL , Pidala J , Schoemans HM , Lawitschka A , Flowers ME , Cowen EW , et al . National Institutes of Health Consensus Development Project on Criteria for Clinical Trials in Chronic Graft-Versus-Host Disease: IIa. The 2020 Clinical Implementation and Early Diagnosis Working Group Report . Transplant Cell Ther . 2021 ; 27 : 545 - 57 . 33839317 10.1016/j.jtct.2021.03.033 PMC8803210
23. Redondo S , De Dios A, Gomis-Pastor M , Esquirol A , Aso O , Triquell M , et al . Feasibility of a New Model of Care for Allogeneic Stem Cell Transplantation Recipients Cacilitated by eHealth: The MY-Medula Pilot Study . Transplant Cell Ther . 2023 ; 29 : 385.e1 - 8 . 10.1016/j.jtct.2023.03.016 36948273