Abstract
Background Physical, psychological, and social aspects of health-related quality of life (HRQOL) among head and neck cancer (HNC) survivors may be more affected during the COVID-19 pandemic than before the pandemic. However, the impact is not yet understood well.
Methods Prospectively collected data from the NETherlands QUality of life and BIomedical Cohort study in HNC were used. All patients were diagnosed and treated before the COVID-19 pandemic. Patient reported outcome measures (PROMs) collected 24 and 36 months after treatment (M24 and M36) were compared between survivors who completed both assessments before the COVID-19 pandemic and those who completed M24 before but M36 during the pandemic. Personal, clinical, physical, psychological, social, and lifestyle characteristics of the survivors assessed at baseline or M24 were investigated as potential effect modifiers.
Results In total, 318 HNC survivors were included, of which 199 completed both M24 and M36 before the COVID-19 pandemic and 119 completed M24 before but M36 during the pandemic. Changes in HRQOL between 24 and 36 months follow-up did not differ between the two groups for any of the PROMs. However, in some subgroups of HNC survivors the COVID-19 pandemic negatively affected the course of HRQOL for several PROMs while it positively affected the course of HRQOL for other PROMs.
Conclusions The COVID-19 pandemic did not affect HRQOL in HNC survivors in general, but some subgroups were affected in a positive and others in a negative way.
Funding This work was supported by the Dutch Cancer Society [grant number VU 2013–5930] and the Dutch Cancer Society, Alpe Young Investigator Grant [grant number 12820].
Introduction
Head and neck cancer (HNC) and its treatment negatively affect the physical, psychological, and social aspects of health-related quality of life (HRQOL) not only before and during treatment but also in long-term survivors [1-6]. It may be that HRQOL of HNC survivors is even more affected in the setting of the COVID-19 pandemic which had an enormous health care, societal and economic impact, but this is not yet understood very well. Previous studies among cancer patients reported a negative impact of the COVID-19 pandemic on HRQOL related to physical [7, 8], psychological [7, 9-12], and social functioning [7, 10, 12]. Other studies found no effect of the COVID-19 pandemic on HRQOL [13, 14] or showed an improvement in physical, role, and social functioning [12] and less loneliness [10]. Several studies investigated possible effects of the COVID-19 pandemic on HRQOL among HNC patients specifically [15-19]. Gallo et al. reported a negative effect on physical, role and emotional functioning and no effect on cognitive and social functioning, global quality of life, or cancer related symptoms [15]. Hamilton et al. found a negative impact on cancer related symptoms [16]. Kirtane et al. conducted a qualitative study and found a negative impact of the COVID-19 pandemic on psychological distress and social isolation [17]. In contrast, Rodrigues-Oliveira et al. found no differences on symptoms of anxiety and depression between HNC patients during radiotherapy before or after the COVID-19 pandemic [18]. Büntzel et al. reported a negative effect on physical and psychological function, and isolation, and a positive effect on relations with family and nature [19].
An explanation for these different and sometimes contradictory findings on the impact of the COVID-19 pandemic among cancer patients might be that most investigators used cross-sectional study designs, and compared the results to reference values or historical cohorts [9-11, 13-19]. Only few studies investigated HRQOL longitudinally, in which part of the measurements were carried out before the COVID-19 pandemic and others during the pandemic [8, 12] or in which cancer patients were prospectively followed during the COVID-19 pandemic [7]. Another explanation for the various previous findings might be that the COVID-19 pandemic affected some patients more than others. Factors that might moderate the impact of the COVID-19 pandemic are age, sex, educational level, marital status, household composition, type and stage of cancer, treatment intent, physical health, job security, and pre-existing psychological problems [7, 9, 14]. Thus far, studies longitudinally investigating the effect of the COVID-19 pandemic on HRQOL among HNC survivors are lacking. The aim of this study is to prospectively investigate the impact of the COVID-19 pandemic on HRQOL among HNC survivors. The following two research questions are addressed: 1) Did the COVID-19 pandemic change the course of HRQOL over time (disease and tumor specific HRQOL, physical activity, symptoms of distress, anxiety, and depression, fear of cancer recurrence, and loneliness), 2) Can subgroups of HNC patients (in terms of personal, clinical, psychological, physical, social, lifestyle, and disease-related factors) be identified in which the course of HRQOL is positively or negatively affected as a consequence of the COVID-19 pandemic.
Methods
Study design and patients
Data and samples were used from the NETherlands QUality of life and BIomedical Cohort study in HNC (NET-QUBIC), a prospective cohort study among 739 HNC patients. Patients were recruited between March 2014 and June 2018. Patients were included before start of treatment (baseline) and data was collected at baseline (T0) and 3, 6, 12, 24, and 36 months (M3, M6, M12, M24 and M36) after end of treatment. Data was derived from an electronic case report form (eCRF) designed for NET-QUBIC, patient reported outcome measures (PROMs) and fieldwork assessments. Newly diagnosed HNC patients were included in NET-QUBIC if they were i) 18 years or older, ii) treated with curative intent for cancer of the oral cavity, oropharynx, hypopharynx, larynx, or unknown primary, and iii) able to write, read, and speak Dutch. Patients were excluded if they i) were unable to understand the questions or test instructions, ii) had severe psychiatric comorbidities (i.e. schizophrenia, Korsakoff’s syndrome, dementia), or iii) were unable to understand informed consent. In the current study, HNC survivors were included if they completed at least one PROM at M36. Furthermore, only the T0, M24 and M36 PROMs were used. PROMs completed before March 14th 2020 were considered ‘before COVID-19’ while PROMs completed on or after March 14th 2020 were considered ‘during COVID-19’. Consent procedures were approved by the Medical Ethical Committee of Amsterdam UMC location VUmc (METc VUmc 2013.301 (A2018.307)-NL45051.029.13)) and all participating hospitals and followed the Dutch Medical Research Involving Human Subjects Act. All included survivors signed informed consent.
Details of the key components of NET-QUBIC were published previously: the study population (including retention, attrition and potential selection bias), eCRF, the outcome assessment protocol, biobanking protocol, data management (collection and storage), and data and sample dissemination procedures [20, 21]. The STROBE guidelines were adhered when reporting the results of this study [22]. As this was a not pre-planned post-hoc analysis of the original NET-QUBIC, no formal sample size calculation for the current study has been performed.
Outcome measures
Disease specific and tumor specific HRQOL were measured with the European Organisation for Research and Treatment of Cancer (EORTC) Quality of Life Questionnaire (QLQ)-C30 and EORTC-QLQ-HN35, respectively. The EORTC QLQ-C30 is a 30-item questionnaire and consists of one global quality of life scale (QL), five functional scales (physical, role, emotional, cognitive and social functioning), three symptom scales (fatigue, nausea and vomiting, and pain) and six single item symptoms (dyspnoea, insomnia, appetite loss, constipation, diarrhoea and financial difficulties). A sum score (SumSC) is based on the five functional scales, the three symptom scales and five of the six single items (financial toxicity is not included) [23-25]. The EORTC-QLQ-HN35 consists of seven HNC specific symptoms (pain, swallowing, senses, speech, social eating, social contact, and sexuality) and ten single item symptoms (problems with teeth, dry mouth, sticky saliva, cough, opening the mouth wide, weight loss, weight gain, use of nutritional supplements, feeding tubes, and painkillers) [26]. QLQ-C30 and QLQ-HN35 Items were scored on a 4-point Likert scale and scores for each of the subscales range from 0 to 100, where higher scores on QL, SumSC and the functional scales indicate better HRQOL and functioning while higher scores on the symptom scales indicate more symptoms.
Physical activity was measured with the Physical Activity Scale for the Elderly (PASE). The PASE is a 13-item questionnaire measuring duration and frequency of leisure time, household and work-related physical activities [27, 28]. Subscale scores for each of the domains were calculated as well as a total score. Higher scores on the subscales and total score indicate more physical activity. The total activity score was also categorized as very poor, poor, fair, good, very good and excellent.
Distress and symptoms of anxiety and depression were measured using the Hospital Anxiety and Depression Scale (HADS). The HADS is a 14-item questionnaire measuring emotional distress and includes a total scale (HADS-T) and an anxiety (HADS-A) and depression (HADS-D) subscale [29]. Items were scored on a 4-point Likert scale and scores for each of the subscales range from 0 to 21, where higher scores indicate higher extent of distress, depression or anxiety symptoms.
Fear of cancer recurrence was measured with the Cancer Worry Scale (CWS). The CWS is an 8-item questionnaire measuring concerns about developing cancer or developing cancer again and the effect of these concerns on daily life [30]. Items were scored on a 4-point Likert scale, and a total scale score is calculated by summing all items, resulting in a total scale score ranging from 8 to 32. A higher score indicates higher extent of fear of cancer recurrence.
Loneliness was measured with the De Jong Gierveld Loneliness Scale (De Jong Gierveld). The De Jong Gierveld is an 11-item questionnaire measuring emotional and social loneliness [31, 32]. Items were scored on a 3-point Likert scale. Scores for emotional loneliness range from 0 to 6 and scores for social loneliness range from 0 to 5, where higher scores indicate higher loneliness. A total loneliness score was calculated by summing the scores of the two scales. The total loneliness score was also categorized as not lonely (score 0 to 2), moderate (score 3 to 8), severe (score 9 or 10) and very severe (score 11) [33].
Influencing factors
Data on personal, clinical, physical, psychological, social, and lifestyle characteristics were collected from eCRF data, PROMs and fieldwork assessments.
Personal factors (assessed at baseline) included age, sex, educational level (low, middle or high), living status (alone or cohabiting), marital status (married or not married), and personality. Personality was assessed by the extraversion subscale of the 60-item NEO Five Factor Inventory (NEO-FFI) questionnaire [34]. Items are scored on a 5-point Likert scale. The extraversion subscale consists of 12 item and ranges from 12 to 60, where a higher score indicates a higher level of extraversion.
Clinical factors included tumor location (oral cavity, oropharynx, hypopharynx, larynx or unknown primary), tumor stage (I/II or III/IV), treatment modality (single or multimodality treatment), World Health Organization (WHO) performance status (0, able to carry out all normal activity without restriction or ≥1, restricted in normal activities), comorbidity and cancer progression at 24 months follow-up. Comorbidity was assessed by the 27-item Adult Comorbidity Evaluation-27 (ACE-27) Index, which categorizes comorbidity as none, mild, moderate and severe [35]. Cancer progression status between end of treatment and M24 was categorized as residual disease, recurrence and/or second primary tumor at M24 or none of those at M24.
Physical impairments in instrumental activities in daily life was assessed by the Instrumental Activities of Daily Living (IADL) questionnaire [36].
Psychological characteristics included the presence of a major depressive disorder in the past year at M24, and presence of a lifetime major depressive disorder at M24. Presence of a major depressive disorder in the past year was assessed with the Composite International Diagnostic Interview (CIDI), which is based on DSM-IV criteria [37]. The CIDI was assessed yearly during the fieldwork assessment, i.e. at T0, M12 and 24. At T0, presence of a lifetime major depressive disorder was also assessed. A lifetime major depressive disorder was scored as present in case it was already present at baseline or when a major depressive disorder in the past year was diagnosed during the CIDI at M12 and/or M24.
Having paid work (yes or no) was assessed by the iMTA productivity cost questionnaire at baseline and at M24.
Lifestyle-related factors included excessive alcohol consumption (categorized as no or yes (at least 14 (women) or 21 (men) glasses of alcohol per week)), smoking behavior (categorized as current smoker or never smoker and former smoker), and body mass index (BMI) at 24 months.
Statistical analyses
Data is described by number and frequency in case of categorical variables and by mean and standard deviation (SD) in case of continuous variables. Differences in baseline characteristics between survivors who completed M36 before COVID-19 and survivors who completed M36 during COVID-19 as well as between survivors included in the current study and the other NET-QUBIC participants were analyzed by the chi-square test or the independent sample t-test, depending on the distribution of the variable. Difference in change from M24 to M436 of continuous PROMs between the groups were analyzed by linear mixed effects models, with fixed effect for group, follow-up measurement (M24 or M36) and their two-way interaction and a random intercept for subject. Differences in change of dichotomous or categorical PROMs were analyzed by generalized estimating equations (GEE), with a logit-link function (dichotomous PROMs) or cumulative logit-link function (categorical PROMs). Reported effect sizes, with corresponding 95% confidence intervals (CI) for the change between M24 and M36 include the difference between groups in change between M24 and M36 based on the estimated marginal means for continuous PROMs and odds ratios (OR) for dichotomous and categorical outcomes.
The modifying effect of potential influencing factors for the change between M24 and M36 was also analyzed using linear mixed effects model, including fixed effects for group, measurement and the potential effect modifier, all two-way interactions and the three-way interaction as well as a random intercept for subject. In case of multiple effect modifiers for one PROM, analyses were repeated after stratification for the factor with the lowest p-value (below 0.05) for the three-way interaction. Continuous effect modifiers were stratified by median split. Stratification was only done in case each stratum and each group (completely before COVID or partly during COVID) contained at least 20 survivors. All analyses were performed in SPSS version 28 (IBM Corp., Armonk, NY, USA). The two-sided significance level was set at 0.01, to account for the large number of PROMs analyzed.
Results
Study population
Of the 739 participants included in NET-QUBIC, 487 were known to be alive and included in NET-QUBIC at M24 of which 345 survivors completed at least one PROM questionnaire at the M36 assessment, 199 before the COVID-19 pandemic and 146 during the pandemic. Twenty-seven of these 146 survivors also filled in the questionnaires at M24 during the pandemic, and were excluded form analyses resulting in 119 survivors in the group with the M36 assessment during the pandemic (Figure 1). Baseline characteristics, stratified by group, are displayed in Table 1. There were no significant difference in characteristics between the groups. HNC survivors included in the current study (n=318) had lower comorbidity and less often a lifetime major depressive disorder at M24 compared to HNC survivors who were not included (n=142) in the current study (Supplementary Table 1).
The course of HRQOL in relation to the COVID-19 pandemic
There were no statistically significant differences in the change of PROMs between 24 and 36 months follow-up between survivors who completed all PROMs before COVID-19 and those who completed the M36 assessment during COVID-19 (Table 2 and Supplemental Tables 2 and 3).
For 31 of the 46 PROMs, one or more factors were identified that showed a different effect of the COVID-19 pandemic on the change of HRQOL between M24 and M36 (p<0.05; Supplemental Table 4). After stratification, the most important factors (with a p-value in at least one of the subgroups <0.01) were: gender, living status at baseline, treatment modality, comorbidity, cancer progression status at M24, presence of a lifetime major depressive disorder at M24, and BMI at M24 (Figure 2; Supplemental Table 5). Females deteriorated in PASE total score during COVID-19 while they improved before COVID-19 (categorical) (OR: 0.25, 95% CI: [0.093; 0.68], p=0.007). Survivors living alone at baseline improved on EORTC-HN35 painkiller use during COVID-19 while they deteriorated before COVID-19 (OR: 3.5, 95% CI [1.5; 8.1], p=0.003). HNC survivors treated by multimodality treatment who completed M36 during COVID-19 showed a deterioration on the EORTC-C30 loss of appetite scale whereas survivors who completed M36 before COVID-19 improved (effect size: −8.5, 95% CI [−13.7; −3.4], p=0.001). Survivors with residual disease, recurrence and/or second primary tumor at M24 who completed M36 during COVID-19 improved on the EORTC-C30 financial problems scale while survivors before COVID-19 deteriorated (effect size: 14.9, 95% CI: [5.4; 24.5], p=0.003). Survivors with moderate and severe comorbidity deteriorated on the EORTC-HN35 swallowing scale during COVID-19 while survivors before COVID-19 improved (effect size: −10.2, 95% CI [−16.5; −3.9], p=0.002). Survivors with a lifetime major depressive disorder present at M24 deteriorated in fear of recurrence during COVID-19 while they improved before COVID-19 (effect size: −3.0, 95% CI: [−5.1; − 0.84], p=0.007). Finally, survivors with a BMI at M24 above the median improved on emotional loneliness during COVID-19 while they deteriorated before COVID-19 (effect size: −0.72, 95% CI: − 1.2; −0.28], p=0.002).
Discussion
Over all, the COVID-19 pandemic did not influence the change in HRQOL between 24 and 36 months after treatment in HNC survivors. However, in some subgroups of HNC survivors the COVID-19 pandemic had a negative effect on some PROMs while in other groups there was a positive effect. The COVID-19 pandemic had a negative effect on females (worsening of physical activity), survivors treated with multimodality treatment (worsening of appetite loss), survivors with comorbidity (worsening of swallowing problems), survivors with a history of a major depressive disorder (worsening of fear of recurrence), and survivors with a high BMI (worsening of emotional loneliness). The COVID-19 pandemic had a positive effect on survivors living alone (decrease of painkiller use) and on survivors with disease progression (decrease of financial problems). Previous studies also reported that sex and pre-existing psychological problems moderate the impact of the COVID-19 pandemic [7, 9, 14]. Furthermore, these studies suggested that also age, educational level, marital status, household composition, type and stage of cancer, treatment intent, and physical health might moderate the effect of the COVID-19 pandemic but this was not confirmed by the current study [7, 9, 14].
An explanation for the negative and positive effects of the COVID-19 pandemic among subgroups of HNC survivors might be an altered and unequal access to follow-up or supportive care during the COVID-19 pandemic [7, 8, 17, 38, 39]. Among colorectal cancer survivors in follow-up care during the COVID-19 pandemic, role, emotional and social functioning, fatigue, sleep disturbance, and anxiety was worse in those survivors who had hospital visits canceled, postponed, or changed into digital care, compared with survivors without changes in their cancer care planning [13]. In contrast, a systematic review on studies investigating the effect of a reduction in follow-up frequency during the COVID-19 pandemic among breast cancer patients, showed no adverse effect on HRQOL (nor on survival) but improved cost-effectiveness of follow-up care. Four RCTs that investigated follow-up on-demand versus scheduled follow-up visits found no statistically significant differences in HRQOL [40]. What we learned from the COVID-19 pandemic is that supportive care services were capable to make significant changes in the provision of their care in a short period of time, and that eHealth was often used as integrated part of supportive care. However, it is known that some cancer survivors benefit more from eHealth than others [41]. HNC survivors may also have used peer support but a survey among health care professionals in the UK demonstrated that different types of peer support are available but that referral to peer support is complex and divers [42].
The key strength of this study is the prospective longitudinal research design. A limitation is that all HNC survivors were included and finished their primary treatment for HNC before the COVID-19 pandemic, and the results cannot be generalized to those who were diagnosed and treated during the COVID-19 pandemic. Also, for some variables, data were only available at baseline while their situation may have changed two years later (e.g. marital status). Moreover, we investigated many PROMs, and although we accounted for this by setting the significance level to 0.01, some of our results could have been significant by chance. Finally, we do not know whether participating HNC survivors were diagnosed with COVID-19. Cancer patients in general have higher odds to develop severe COVID-19 and to die of the consequences [43].
In conclusion, the course of HRQOL between 24 and 36 months after treatment in HNC survivors in general was not affected by the COVID-19 pandemic. However, the COVID-19 pandemic did change the course of some HRQOL domains or symptoms over time in some subgroups of HNC survivors. The development of personalized supportive care programs including regular care, eHealth, and peer support, tailored to the needs of the individual survivor may help to overcome disparities among HNC survivors.
Data Availability
The datasets geneThe dataset generated and analyzed (including a deindentified version of it) during the current study is not publicly available as the collection and integration of large amounts of personal, biological, genetic and diagnostic information precludes open access to the NET-QUBIC research data. Data are available from the corresponding author on reasonable request. On the NET-QUBIC website (www.kubusproject.nl) is described how NET-QUBIC data are made available for the research community. In short, a research proposal has to be submitted to the NET-QUBIC steering committee. The researcher has to be a member of the NET-QUBIC consortium, or a consortium member has to collaborate in the project. The NET-QUBIC steering committee will judge the proposal. The SPSS syntax of the analyses in this study are publically available on GitHub: https://github.com/b-lissenberg/COVID-and-cancer.
Funding
This work was supported by the Dutch Cancer Society [grant number VU 2013–5930] and the Dutch Cancer Society, Alpe Young Investigator Grant [grant number 12820]. The funding body had no role in the design of the study and collection, analysis, and interpretation of data nor in writing the manuscript.
Competing interests
BLW: none
FJ: received funding from the Dutch Cancer Society, Alpe Young Investigator Grant [grant number 12820]. All payments were made to the institution.
RJBdJ: none
FL: received funding from the Dutch Cancer Society [grant numbers VU2017-8288 and 11839]. All payments were made to the institution.
CRL: none
SFO: received funding from Celidex Therapeutics, payment or honoraria for lectures, presentations, speakers bureaus, manuscript writing or educational events from Merck and participated in the Data Safety Monitoring Board or Advisory Board of GenMab and Bristol Myers Squibb. All payments were made to the institution.
RPT: none
IMVdL: received funding from the Dutch Cancer Society (grant number VU 2013-5930). All payments were made to the institution.
Data availability
The dataset generated and analyzed (including a deindentified version of it) during the current study is not publicly available as the collection and integration of large amounts of personal, biological, genetic and diagnostic information precludes open access to the NET-QUBIC research data. Data are available from the corresponding author on reasonable request. On the NET-QUBIC website (www.kubusproject.nl) is described how NET-QUBIC data are made available for the research community. In short, a research proposal has to be submitted to the NET-QUBIC steering committee. The researcher has to be a member of the NET-QUBIC consortium, or a consortium member has to collaborate in the project. The NET-QUBIC steering committee will judge the proposal. The SPSS syntax of the analyses in this study are publically available on GitHub: https://github.com/b-lissenberg/COVID-and-cancer.