Distinct employment interference profiles in patients with breast cancer prior to and for 12 months following surgery.

PURPOSE: To identify subgroups of female breast cancer patients with distinct self-reported employment interference (EI) profiles and determine which demographic, clinical, and symptom characteristics, and quality of life outcomes were associated with subgroup membership.
METHODS: Women with breast cancer (n = 385) were assessed for changes in EI over ten times, from prior to, through 12 months after breast cancer surgery. Latent profile analysis (LPA) was used to identify subgroups of patients with distinct EI profiles.
RESULTS: Three distinct EI profiles (i.e., None - 26.2% (n = 101), Low - 42.6% (n = 164), High - 31.2% (n = 120)) were identified. Compared to the None and Low groups, patients in the High group were more likely to be younger. Higher proportions in the High group were non-White, pre-menopausal prior to surgery, had more advanced stage disease, had received an axillary lymph node dissection, had received neoadjuvant chemotherapy, had received adjuvant chemotherapy, and had a re-excision or mastectomy on the affected breast within 6 months after surgery. In addition, these patients had lower quality of life scores. Compared to the None group, the High group had higher levels of trait and state anxiety, depressive symptoms, fatigue and sleep disturbance and lower levels of cognitive function.
CONCLUSIONS: This study provides new knowledge regarding EI profiles among women in the year following breast cancer surgery. The non-modifiable risk factors (e.g., younger age, being non-White, having more advanced stage disease) can inform current screening procedures. The potentially modifiable risk factors can be used to develop interventions to improve employment outcomes of breast cancer patients.

Introduction

In 2020, over 270,000 new cases of female breast cancer were diagnosed in the United States [1]. A breast cancer diagnosis and its treatment have significant short- and long- term impacts on patients’ employment [2]. Depending on an individual’s situation, breast cancer patients report varying levels of absenteeism from work, reduced work ability, limited work performance and career progression, and at times, termination of employment [2, 3]. In addition, compared to men, women are in relatively less stable occupations, as well as in lower income and more casual employment, that results in more acute employment interference (EI) from their cancer [4]. The ability to return to work and to satisfactorily engage in work contributes to continued insurance benefits, financial well-being, sense of identity and accomplishment, and a return to a resemblance of normal [5].

The impacts of cancer on employment is a global issue [2]. In a Swedish cross sectional study of 756 women with breast cancer, 56% were still on sick leave at 4 to 6 weeks after surgery [6]. In a French study [7], over 80% of breast cancer survivors who were in paid work prior to their diagnosis had a median sick leave duration of 10.8 months. Older age, lower educational level, receipt of chemotherapy (CTX) or radiotherapy, lymphoedema, and insufficient support from employers or colleagues were associated with a delay in return to work [7]. For women who returned to work, diminished performance was associated with pain [8], anxiety [9-13], fatigue [9-13], depression [9-13], and cognitive dysfunction [10, 13, 14]. These limitations affected productivity, duration of the work day, and/or changes in work roles and responsibilities [5].

In the survivorship literature, most studies have focused on how the aforementioned factors affect return to work as a binary outcome [2, 15]. This binary outcome captures only extreme consequences, namely: no loss or complete loss of employment. In contrast, an assessment of the overall extent of employment interference (EI) has the advantage of evaluating the spectrum of potential negative outcomes. This approach may provide a more accurate picture of the impact of cancer on work and employment. Prior to our recent publication [16], no studies had examined patients’ subjective experiences with the extent of EI, as well as changes in their experiences over time and associated predictors. In our recent longitudinal study [16], we used hierarchical linear modeling (HLM) to evaluate for inter-individual differences in self-reported EI in 387 women with breast cancer over 12 months. In this study [16], factors associated with higher levels of self-reported EI were younger age, a lower annual household income, higher pain intensity, higher sleep disturbance, and the receipt of an axillary lymph node dissection (ALND), a sentinel lymph node biopsy (SLNB), adjuvant CTX, complementary or alternative therapies, or a re-excision or mastectomy.

While our previous study [16] provides some insights into the characteristics associated with EI in patients with breast cancer, the statistical approach used did not allow for the identification of subgroups of patients with distinct EI profiles. The use of person-centered analytic approaches, like latent profile analysis (LPA), allows for the identification of these subgroups, as well as for the identification of demographic and clinical characteristics associated with subgroup membership. This type of analysis can provide information on high risk patients who can be targeted for appropriate interventions. Therefore, the purposes of this study were to use LPA to identify subgroups of breast cancer patients with distinct self-reported EI profiles and to determine which demographic, clinical, and symptom characteristics, as well as quality of life (QOL) outcomes, were associated with each subgroup. In order to be able to determine risk factors for EI, all of the patients in the sample were included regardless of their employment status at enrollment.

Methods

Patients and settings

The methods used in this analysis are described in detail in previous publications [17-20] and are summarized in this methods section. This descriptive, longitudinal analysis is part of a larger National Cancer Institute-funded study that evaluated for neuropathic pain and lymphedema in women who underwent breast cancer surgery [17-19]. Patients were recruited from Breast Care Centers located in a Comprehensive Cancer Center, two public hospitals, and four community practices in Northern California. Patients were eligible to participate if they were women > 18 years of age who would undergo breast cancer surgery on one breast; were able to read, write, and understand English; agreed to participate; and gave written informed consent. Patients were excluded if they were having breast cancer surgery on both breasts and/or had distant metastases at the time of diagnosis.

Study procedures

The study was approved by the Committee on Human Research at the University of California, San Francisco and by the Institutional Review Board at each of the study sites. During the patient’s preoperative visit, a clinical staff member explained the study to the patient, determined her willingness to participate, and introduced her to the research nurse. The research nurse met with the women, determined eligibility, and obtained written informed consent prior to surgery. After obtaining consent, patients completed the enrollment questionnaires an average of 4 days prior to surgery and at one, two, three, four, five, six, eight, 10, and 12 months after surgery. The research nurse met with the patients in the Clinical Research Center or in their homes. Patients’ medical records were reviewed for disease and treatment information.

Instruments

At enrollment, patients completed a demographic questionnaire, the Karnofsky Performance Status (KPS) scale [21], and the Self-Administered Comorbidity Questionnaire (SCQ) [22]. The KPS scale, that ranged from 40 (disabled, requires special care and assistance) to 100 (normal, no complaints; no evidence of recurrence) was used to evaluate functional status [21]. The SCQ evaluates the occurrence of, treatment for, and impact of 13 common medical conditions. The total SCQ score ranges from 0 to 39 [22].

For this study, EI was assessed using a single item from the Quality of Life Scale-Patient Version (QOL-PV) [23], namely, “To what degree has your illness or treatment interfered with your employment?”. Patients were asked to rate this item using a 0 (no problem) to 10 (severe problem) numeric rating scale (NRS). No gold standard exists to measure employment outcomes [2]. Our approach allowed us to assess patients’ subjective perceptions of the overall impact of cancer and its treatment on their employment [16]. Another advantage of the evaluation of self-reported EI with a single item rather than return to work (yes/no) is that the latter outcome is rather simplistic and does not take into account external factors that are unrelated to cancer but may influence employment, such as aging and changes in outlook on life [5, 24]. The phrasing of the EI item has the advantage of having the patient focus on how their employment was disrupted specifically by cancer and its treatment.

The Spielberger State-Trait Anxiety Inventories (STAI-S, STAI-T) were used to assess the patient’s transitory emotional response to a stressful situation and her predisposition to anxiety, respectively. Scores for each scale are summed and can range from 20 to 80. Higher scores indicate greater anxiety. Cut-off scores of > 31.8 and > 32.2 indicate high levels of trait and state anxiety, respectively [25]. Both inventories have well-established validity and reliability [26, 27]. In this study, Cronbach’s alphas for the STAI-T and STAI-S were .88 and .95, respectively.

The 20-item Center for Epidemiological Studies Depression Scale (CES-D) was used to evaluate the major symptoms in the clinical syndrome of depression. Scores can range from 0 to 60, with scores of > 16 indicating the need for individuals to seek clinical evaluation. The CES-D has well-established validity and reliability [28-30]. In this study, its Cronbach’s alpha was 0.90.

The 18-item Lee Fatigue Scale (LFS) was used to assess physical fatigue and energy [31]. Each item was rated on a 0 to 10 NRS. Total fatigue and energy scores were calculated as the mean of the 13 fatigue and the five energy items. Higher scores indicate greater fatigue severity and higher levels of energy. Patients were asked to rate each item based on how they felt “right now”. Cut-off scores of > 4.4 and < 4.8 indicate high levels of fatigue and low levels of energy, respectively [32]. The LFS has well established validity and reliability [31]. In this study, Cronbach’s alphas for fatigue and energy scales were .96 and .93, respectively.

The 16-item Attentional Function Index (AFI) was used to measure attentional function [33]. Each item was rated on a 0 to 10 NRS. A higher mean score indicates greater capacity to direct attention [33, 34]. Scores are grouped into categories of attentional function (i.e., < 5.0 low function, 5.0 to 7.5 moderate function, > 7.5 high function) [35]. The AFI has well established validity and reliability [34, 36]. In this study, its Cronbach’s alpha was .95.

The 21-item General Sleep Disturbance Scale (GSDS) was used to assess the quality of sleep in the past week. Each item was rated on a 0 (never) to 7 (everyday) NRS. The GSDS total score is the sum of 21 items that can range from 0 (no disturbance) to 147 (extreme sleep disturbance). A GSDS total score of > 43 indicates a significant level of sleep disturbance [37]. The GSDS has well established validity and reliability [38, 39]. In this study, its Cronbach’s alpha was .86.

The occurrence of breast pain prior to surgery was determined by asking the question “Are you experiencing pain in your affected breast?” If women responded yes, they rated their average and worst pain using a 0 (no pain) to 10 (worst imaginable pain) NRS.

The 41-item QOL-PV was used to assess four domains of QOL (i.e., physical well-being, psychological well-being, social well-being, spiritual well-being). Items are rated on a 0 to 10 NRS. Mean subscale and total scores were calculated. Higher scores indicate better QOL. The QOL-PV has well established validity and reliability [40, 41]. Cronbach’s alphas for the QOL-PV physical well-being, psychological well-being, social well-being, and spiritual well-being subscales, as well as the total QOL scale, were: .80, .86, .80, .63, and .86, respectively.

Data analysis

Descriptive statistics and frequency distributions were computed for sample characteristics, symptom severity scores, and QOL scores using SPSS version 27 (IBM Corporation, Armonk, NY). As previously described [42], unconditional LPA was used to identify the profiles of self-reported EI that characterized unobserved subgroups (i.e., latent classes) of patients over the 12 months of the study. First, all of the patients who reported a zero for the EI item across the 10 assessments were categorized in the None group. Then, we identified subgroups of patients based on their profiles of means across the 10 assessments for the EI item from the QOL-PV. In order to incorporate the expected correlations among the repeated measures, we included covariance among the EI scores that were up to four occasions apart (i.e., a covariance structure with a lag of four). In this way, we retained the within-person correlation among the self-reported EI scores, while we focused on the patterns of means that distinguished among the latent classes. We limited the covariance structure to a lag of four to accommodate the expected reduction in correlation that would be introduced by decreased stability in the EI ratings as the separation of months increased and to reduce model complexity.

Estimation was carried out with full information maximum likelihood with standard errors and a Chi-square test that are robust to non-normality and non-independence of observations (“estimator = MLR”). Model fit was evaluated to identify the best solution that characterized the observed latent class structure with the Bayesian Information Criterion (BIC), the Vuong-Lo-Mendell-Rubin (VLMR) likelihood ratio test the K vs. K− 1 model, entropy, and latent class percentages that were large enough to be reliable (i.e., likely to replicate in new samples) [43]. Missing data were accommodated with the use of the Expectation-Maximization algorithm [44]. Mixture models, like LPA, are known to produce solutions at local maxima. Therefore, our models were fit with from 800 to 1600 random starts. This approach ensured that the estimated model was replicated many times and was not due to a local maximum. Estimation was done with Mplus Version 7.2 [45].

After identifying the latent class solution that best fit the data, differences among EI groups, in demographic and clinical characteristics, symptom scores, and QOL scores, obtained at enrollment, were evaluated using analyses of variance, Chi-square, and Kruskal-Wallis analyses. A p-value of < 0.05 was considered statistically significant. Post-hoc contrasts were done using a Bonferroni corrected p-value of < 0.017 (0.05/3 pairwise comparisons).

Results

Latent profile analysis

Data from 385 patients with breast cancer were used in the LPA. As shown in Fig. 1, 26.2% of patients (n = 101) did not report any EI for all of the assessments and were named the None class. Among the remaining 284 patients, two distinct latent classes were identified. A two-class model was selected because its BIC was lower than the BIC for the 1-class solution (Table 1). The VLMR was statistically significant for the 2-class solution, indicating that two classes fit the data better than one class. In addition, the VLMR was not significant for the 3-class solution, indicating that too many classes were extracted.

Fig. 1

Employment interference trajectories for patients in each of the latent classes

Table 1

Employment Interference: Latent Profile Solutions and Fit Indices for One- Through Three-Class Solutions

Model	LL	AIC	BIC	Entropy	VLMR
1 Class	− 5115.38	10,318.75	10,479.30	n/a	n/a
2 Class	− 4937.42	9982.84	10,179.88	.92	355.91+
3 Class	− 4885.99	9899.98	10,133.52	.87	102.85ns

+p < .0001

Abbreviations: AIC Akaike’s Information Criterion, BIC Bayesian Information Criterion, LL log-likelihood, n/a not applicable for one class, ns not significant, VLMR Vuong-Lo-Mendell-Rubin likelihood ratio test for the K vs. K-1 model

As shown in Fig. 1, the largest proportion of the patients were classified in the Low class (n = 164, 42.6%). This group reported a mean EI score of 2.4 (+ 2.9) at enrollment, that increased slightly at month 1 and then decreased over the 12 months of the study. The second class that consisted of 31.2% (n = 120) of the women had a mean EI score of 6.2 (+ 3.5) at enrollment and was named the High class. This group’s EI scores were in the moderate range at enrollment and increased and remained in the moderate to high range in the 12 months following surgery.

Differences in demographic and clinical characteristics

As shown in Table 2, no differences were found among the three groups in body mass index, marital status, living arrangements, number of hours worked per week (for those who were employed), time since diagnosis, type of surgery, receipt of SLNB, receipt of breast reconstruction at the time of surgery, or pain in the affected breast. In addition, no differences were found among the three groups in the receipt of radiation therapy, hormonal therapy, or breast reconstruction in the 12 months following surgery.

Table 2

Differences in demographic and clinical characteristics among the employment interference latent classes at enrollment

Characteristic	None (0)26.2% (n = 101)	Low (1)42.6% (n = 164)	High (2)31.2% (n = 120)	Statistics
	Mean (SD)	Mean (SD)	Mean (SD)
Age (years)	62.4 (11.4)	54.8 (10.3)	48.7 (9.5)	F = 48.04, p < 0.001 0 > 1 > 2
Education (years)	15.0 (2.6)	16.1 (2.5)	15.9 (2.7)	F = 5.46, p = 0.005 0 < 1 and 2
Self-Administered Comorbidity Questionnaire score	4.8 (2.9)	3.8 (2.8)	4.3 (2.7)	F = 3.90, p = 0.021 0 > 1
Karnofsky Performance Status score	95.7 (8.0)	93.2 (10.7)	91.3 (11.0)	F = 5.17, p = 0.006 0 > 2
Body mass index (kg/m2)	27.5 (6.7)	26.1 (5.6)	27.4 (6.2)	F = 2.46, p = 0.086
	% (n)	% (n)	% (n)
Ethnicity – Non-White	30.7 (31)	26.2 (43)	47.5 (57)	X2 = 14.66, p = 0.001 0 and 1 < 2
Married/partnered (% yes)	41.6 (42)	38.4 (63)	45.0 (54)	X2 = 1.24, p = 0.537
Lives alone (% yes)	27.7 (28)	22.6 (37)	20.8 (25)	X2 = 1.56, p = 0.459
Currently employed (% yes)	16.8 (17)	68.9 (113)	48.3 (58)	X2 = 67.85, p < 0.001 0 < 1 and 2, 1 > 2
Number of work hours per week (% yes)				X2 = 1.43, p = 0.838
1–35 h/week	38.9 (7)	36.3 (41)	44.8 (26)
36–44 h/week	44.4 (8)	48.7 (55)	39.7 (23)
45h hours/week	16.7 (3)	15.0 (17)	15.5 (9)
Annual household income				KW, p < 0.001 0 < 1, 1 > 2
< $30,000a	23.1 (18)	10.3 (15)	34.4 (33)
$30,000 to <$70,000	28.2 (22)	24.0 (35)	25.0 (24)
$70,000 to <$100,000	14.1 (11)	20.5 (30)	10.4 (10)
> $100,000	34.6 (27)	45.2 (66)	30.2 (29)
Gone through menopause prior to surgery (% yes)	79.2 (80)	65.9 (108)	50.8 (61)	X2 = 19.50, p < 0.001 0 and 1 > 2
Days since cancer diagnosis (mean (SD))	57.3 (63.9)	67.7 (83.7)	84.8 (77.0)	KW, p = 0.079
Days since cancer diagnosis (median)	39.0	36.0	40.50
Stage of disease				KW, p < 0.001 0 and 1 < 2
Stage 0	16.8 (17)	23.8 (39)	14.2 (17)
Stage I	53.5 (54)	37.2 (61)	26.7 (32)
Stage IIa and IIb	25.7 (26)	32.9 (54)	45.8 (55)
Stage IIIa, IIIb, IIIc, and IVa	4.0 (4)	6.1 (10)	13.3 (16)
Type of surgery				X2 = 6.84, p = 0.033 No significant post hoc contrasts
Breast conservation	85.1 (86)	74.4 (122)	85.0 (102)
Mastectomy	14.9 (15)	25.6 (42)	15.0 (18)
Sentinel lymph node biopsy (% yes)	84.2 (85)	84.8 (139)	79.2 (95)	X2 = 1.69, p = 0.430
Axillary lymph node dissection (% yes)	23.8 (24)	31.7 (52)	54.2 (65)	X2 = 24.82, p < 0.001 0 and 1 < 2
Underwent reconstruction at the time of surgery (% yes)	14.9 (15)	25.0 (41)	23.3 (28)	X2 = 4.01, p = 0.135
Received neoadjuvant therapy (% yes)	11.9 (12)	17.7 (29)	30.8 (37)	X2 = 13.37, p = 0.001 0 and 1 < 2
Pain in the affected breast prior to surgery (% yes)	20.8 (21)	27.4 (45)	31.7 (38)	X2 = 3.32, p = 0.191
Received radiation therapy during the 12 months following surgery (% yes)	68.3 (69)	70.1 (115)	75.0 (90)	X2 = 1.35, p = 0.510
Received chemotherapy during the 12 months following surgery (% yes)	26.7 (27)	29.9 (49)	45.8 (55)	X2 = 11.10, p = 0.004 0 and 1 < 2
Received hormonal therapy during the 12 months following surgery (% yes)	64.4 (65)	61.0 (100)	55.8 (67)	X2 = 1.73, p = 0.422
Received complementary therapy in the 12 months following surgery (% yes)	19.8 (20)	39.6 (65)	36.7 (44)	X2 = 11.82, p = 0.003 0 < 1 and 2
Had breast reconstruction in the 12 months following surgery (% yes)	7.9 (8)	14.6 (24)	11.7 (14)	X2 = 2.69, p = 0.260
Had re-excision or mastectomy on the affected breast in the six months following surgery (% yes)	24.8 (25)	25.6 (42)	41.7 (50)	X2 = 10.50, p = 0.005 0 and 1 < 2

areference group

Abbreviations: kg kilograms, KW Kruskal-Wallis, m meters squared, SD standard deviation

Compared to the other two groups, patients in the High group were more likely to be younger, non-White, pre-menopausal prior to surgery, had more advanced stage disease, and had received an ALND, neoadjuvant CTX, adjuvant CTX, and a re-excision or mastectomy on the affected breast within 12 months after surgery. Compared to the Low group, patients in the High group were less likely to be employed and had a lower annual household income prior to surgery. Compared to the None group, patients in the Low group were more likely to be younger, had a lower comorbidity score, and had a higher annual household income. Compared to the None group, patients in the other two groups had more years of education, were more likely to be employed prior to surgery, and more likely to use complementary therapy in the 12 months following surgery.

Differences in symptom severity scores

As shown in Table 3, based on their responses at enrollment, no differences were found among the three groups in energy or pain scores prior to surgery. Compared to the None group, the High group had higher scores for trait and state anxiety, depressive symptoms, fatigue, and sleep disturbance and lower scores for attentional function. Compared to the None group, the Low group had higher fatigue scores at enrollment.

Table 3

Differences in symptom severity scores among the employment interference latent classes at enrollment

Symptoma	None (0)26.2% (n = 101)	Low (1)42.6% (n = 164)	High (2)31.2% (n = 120)	Statistics
	Mean (SD)	Mean (SD)	Mean (SD)
Trait anxiety score (> 31.8)	33.7 (8.6)	35.1 (8.4)	36.7 (9.5)	F = 3.20, p = 0.042 0 < 2
State anxiety score (> 32.2)	38.7 (13.4)	41.7 (13.0)	43.3 (13.0)	F = 3.52, p = 0.031 0 < 2
Center of Epidemiologic Studies-Depression Scale score (> 16.0)	11.2 (9.3)	13.8 (9.2)	15.3 (10.4)	F = 5.01, p = 0.007 0 < 2
Lee Fatigue Scale – fatigue score (> 4.4)	2.3 (2.0)	3.2 (2.2)	3.6 (2.6)	F = 9.32, p < 0.001 0 < 1 and 2
Lee Fatigue Scale – energy score (< 4.8)	5.1 (2.8)	4.9 (2.2)	4.7 (2.5)	F = 0.58, p = 0.558
Attentional Function Index score (< 5.0 is low, 5.0 to 7.5 is moderate, > 7.5 is high)	7.1 (2.0)	6.6 (1.8)	6.1 (2.0)	F = 7.06, p = 0.001 0 > 2
General Sleep Disturbance Scale score (> 43.0)	43.7 (21.9)	48.0 (21.3)	51.8 (20.4)	F = 4.06, p = 0.018 0 < 2
Average pain intensity score in the breast prior to surgery	0.5 (1.5)	0.5 (1.1)	0.7 (1.5)	KW, p = 0.293
Worst pain intensity score in the breast prior to surgery	0.6 (1.6)	0.9 (1.7)	1.2 (2.3)	KW = 0.200

aClinically meaningful cutpoints for symptom severity are in parentheses

Abbreviations: KW Kruskal-Wallis, SD standard deviation

Differences in QOL scores

As shown in Table 4, based on their responses at enrollment, significant differences were found among the three groups in psychological well-being, social well-being, and total QOL scores (None > Low > High). Compared to the other two groups, patients in the High group had lower physical well-being scores. Compared to the Low group, the High group had higher spiritual well-being scores.

Table 4

Differences in quality of life scores among the employment interference latent classes at enrollment

Quality of Life	None (0)26.2% (n = 101)	Low (1)42.6% (n = 164)	High (2)31.2% (n = 120)	Statistics
	Mean (SD)	Mean (SD)	Mean (SD)
Physical well-being	8.5 (1.4)	8.2 (1.4)	7.3 (1.8)	F = 19.05, p < 0.001 0 and 1 > 2
Psychological well being	6.7 (1.8)	5.7 (1.7)	5.1 (1.7)	F = 20.39, p < 0.001 0 > 1 > 2
Social well-being	8.3 (1.3)	7.2 (1.6)	5.5 (2.0)	F = 69.95, p < 0.001 0 > 1 > 2
Spiritual well-being	5.6 (1.8)	5.5 (1.7)	6.1 (2.0)	F = 3.84, p = 0.022 1 < 2
Total quality of life score	7.2 (1.2)	6.5 (1.2)	5.8 (1.4)	F = 29.45, p < 0.001 0 > 1 > 2

Abbreviation: SD standard deviation

Discussion

This longitudinal study is the first to use LPA to identify subgroups of patients with distinct self-reported EI profiles in a large sample of women who were assessed prior to and for 12 months following breast cancer surgery. Cancer- and treatment-related EI is an important issue for women, with 42.6% reporting relatively low levels of EI and 31.2% reporting relatively high levels of EI over the 12 months after surgery. While direct comparisons are difficult because of differences in study measures, our findings are consistent with previous research that highlighted the high prevalence of EI in these patients [2, 5].

Rather than only selecting women who were employed at the pre-operative assessment, our approach of including all women regardless of current employment status allowed us to determine risk factors for EI in the entire breast cancer patient cohort. There are two advantages of using this approach. First, unemployment at the pre-surgical assessment is not necessarily permanent. That is, women who were unemployed could still have the capacity and desire to engage in paid work over the course and after the completion of their treatment. Second, the risk factors identified in this study apply to all women receiving surgery and are applicable for the identification of high risk patients in the clinic. While only 48.8% of these women (n = 188) reported “currently working for pay” at enrollment, the assessments of EI identified that 73.8% of our patients (n = 284) experienced some level of EI over the 12-month period. This discrepancy may be due to changes in women’s employment status before and following surgery. While these findings support our approach of including the entire sample of women in our analysis, they suggest that in future studies, investigators should evaluate employment status, intent to seek employment, and EI at multiple time points following surgery.

Compared to our previous HLM analysis of EI [16], the current LPA identified additional risk factors in those women who experienced high levels of EI (see Additional file 1). In terms of additional non-modifiable risk factors, compared to women with no EI, patients with high levels of EI had a higher level of education; were more likely to be non-White and pre-menopausal prior to surgery; had more advanced stage disease; and had received neoadjuvant CTX. These risk factors are known to impact employment [46]. While findings from a meta-analysis of previous studies found that lower levels of education were associated with an increased level of unemployment following a cancer diagnosis [46], we found that compared to patients with no EI, women in the other two EI groups had more years of education. This inconsistent finding may be related to the older age of the women in the No EI group and that a higher proportion of the patients in this group were retired or working part-time. In contrast, given that the women in the other two groups were younger and better educated, they may have had more demanding jobs and a higher expectation to return to work [5, 46]. Consistent with the literature, a non-White background [3, 47], receipt of neoadjuvant CTX [46], and having more advanced disease [46] were associated with worse employment outcomes in patients with breast cancer.

Consistent with the findings from a systematic review [48], our study highlights a number of potentially modifiable risk factors associated with EI, including a lower functional status, lower levels of energy and cognitive function, as well as higher levels of trait and state anxiety, depressive symptoms, sleep disturbance, and fatigue. While the clear association between a higher symptom burden and poorer work outcomes is not surprising, our findings provide additional evidence of the need to systematically incorporate symptom management interventions in return to work plans. According to a Cochrane review [49], while current return to work interventions may have one or a combination of physical, psychosocial and vocation components, they do not consistently incorporate a symptom management focus. Our findings suggest that the management of the modifiable risk factors including unrelieved symptoms in these interventions may facilitate a patient’s ability to return to work with lower levels of EI. Recently, Alfano and colleagues called for a personalized, tailored approach to address factors that affect cancer survivors’ ability to work during and/or after treatment [50]. To enhance the personalization of interventions, additional research is warranted to determine the most appropriate combination of interventions, based on individual risk factors [50].

Consistent with findings from a systematic review [51], lower levels of EI were associated with better physical, psychological, and social well-being. Of note, fewer studies have examined the relationships between spiritual well-being and work outcomes in cancer survivors [51]. Similar to our previous report on the relationships between financial toxicity and quality of life [20], patients with high levels of EI reported higher levels of spiritual well-being. This finding may be partially explained by the higher proportion of Non-white (47.5%) patients in the High group. Non-white patients are more likely to report higher levels of spiritual well-being [52], especially among African Americans who have a religious and church affiliation [53].

This study has some limitations. First, because no single gold standard is available to measure employment outcomes, EI was measured using a single item. This approach did not allow for a comprehensive assessment of all relevant aspects of work ability and performance (e.g., work satisfaction, time off work, work quality). While our single item is valid and reliable, future studies should perform a more detailed evaluation of EI. In addition, this study was conducted in Breast Care Centers in the United States, which limits the generalizability of the findings to other countries where women may have different experiences with EI. This study was conducted pre-COVID-19. Future studies should investigate how EI profiles and associated factors were or are being impacted by the pandemic.

In conclusion, this longitudinal study with 12 months of follow-up provides new knowledge about distinct subgroups of women with varying levels of EI during and following breast cancer surgery. The non-modifiable and potentially modifiable risk factors associated with EI can assist clinicians to identify high risk patients. The risk factors associated with the High and Low groups can be incorporated into current screening procedures and/or treatment algorithms that would facilitate referrals for multidisciplinary return-to-work and symptom management interventions.

Additional file 1.