Skip to main content

Longitudinal associations of diurnal rest-activity rhythms with fatigue, insomnia, and health-related quality of life in survivors of colorectal cancer up to 5 years post-treatment

Abstract

Background

There is a growing population of survivors of colorectal cancer (CRC). Fatigue and insomnia are common symptoms after CRC, negatively influencing health-related quality of life (HRQoL). Besides increasing physical activity and decreasing sedentary behavior, the timing and patterns of physical activity and rest over the 24-h day (i.e. diurnal rest-activity rhythms) could also play a role in alleviating these symptoms and improving HRQoL. We investigated longitudinal associations of the diurnal rest-activity rhythm (RAR) with fatigue, insomnia, and HRQoL in survivors of CRC.

Methods

In a prospective cohort study among survivors of stage I-III CRC, 5 repeated measurements were performed from 6 weeks up to 5 years post-treatment. Parameters of RAR, including mesor, amplitude, acrophase, circadian quotient, dichotomy index, and 24-h autocorrelation coefficient, were assessed by a custom MATLAB program using data from tri-axial accelerometers worn on the upper thigh for 7 consecutive days. Fatigue, insomnia, and HRQoL were measured by validated questionnaires. Confounder-adjusted linear mixed models were applied to analyze longitudinal associations of RAR with fatigue, insomnia, and HRQoL from 6 weeks until 5 years post-treatment. Additionally, intra-individual and inter-individual associations over time were separated.

Results

Data were available from 289 survivors of CRC. All RAR parameters except for 24-h autocorrelation increased from 6 weeks to 6 months post-treatment, after which they remained relatively stable. A higher mesor, amplitude, circadian quotient, dichotomy index, and 24-h autocorrelation were statistically significantly associated with less fatigue and better HRQoL over time. A higher amplitude and circadian quotient were associated with lower insomnia. Most of these associations appeared driven by both within-person changes over time and between-person differences in RAR parameters. No significant associations were observed for acrophase.

Conclusions

In the first five years after CRC treatment, adhering to a generally more active (mesor) and consistent (24-h autocorrelation) RAR, with a pronounced peak activity (amplitude) and a marked difference between daytime and nighttime activity (dichotomy index) was found to be associated with lower fatigue, lower insomnia, and a better HRQoL. Future intervention studies are needed to investigate if restoring RAR among survivors of CRC could help to alleviate symptoms of fatigue and insomnia while enhancing their HRQoL.

Trial registration

EnCoRe study NL6904 (https://www.onderzoekmetmensen.nl/).

Background

The incidence of colorectal cancer (CRC) is steadily rising, posing an increasing challenge to public health [1]. At the same time, there has been a notable increase in survival rates, reflecting advancements in early detection due to national screening programs, improved treatments, and comprehensive care strategies, resulting in growing numbers of survivors of CRC [2]. Survivors of CRC often experience a decline in their health-related quality of life (HRQoL), due to symptoms caused by the cancer and/or its treatment [3]. Fatigue and insomnia (sleep problems) are among the most common and debilitating symptoms reported after CRC [4]. Up to two-thirds of survivors of CRC are still affected by one of these distressing symptoms up to three years after diagnosis, further exacerbating the physical and emotional burden on former patients with CRC [4, 5]. Therefore, it is important to alleviate symptoms of fatigue and insomnia during the CRC survivorship period to promote their HRQoL.

Partially shared underlying causes of fatigue and insomnia may include chrono-biological factors influencing daily rhythms, such as the timing, amplitude, and regularity of daily rhythms in physical activity [6, 7]. Multiple physiological processes such as sleep, immune function, and metabolism exhibit cyclic patterns approximating 24 h, referred to as circadian rhythms [8, 9]. These circadian rhythms are primarily entrained by light, through the central circadian clock located in the suprachiasmatic nucleus (SCN), but can also be affected by physical activity and dietary intake through peripheral clocks found in almost all organ systems [9, 10]. Previous research in survivors of CRC showed that less sedentary behavior and more physical activity were associated with lower fatigue and higher HRQoL in the first two years after the end of cancer treatment [11,12,13]. Besides the intensity and type of physical activity, the timing of physical activity and rest may also be important for perceived fatigue and HRQoL of this population. Appropriately timed diurnal rest-activity rhythms (RAR) may have the potential to restore synchronicity and improve alignment between central and peripheral clocks [10], thereby potentially also decreasing the risk of negative health outcomes such as fatigue and insomnia after cancer.

In previous research, wrist-worn actigraphy has been employed for measuring several parameters encompassing different aspects of RAR in patients with breast cancer and metastatic CRC [14,15,16,17]. These RAR parameters were related to mean activity levels (mesor), the contrast between peak- and mean activity levels (amplitude), the timing of the peak activity (acrophase), the difference between daytime and nighttime activity (dichotomy index), and consistency of RAR from one day to the next (24-h autocorrelation coefficient). In these patients, lower values for most of these parameters were observed compared to non-cancer adults, indicating a more disrupted RAR [15,16,17]. In turn, a more disrupted RAR has been associated with more fatigue, a lower physical and emotional well-being, and a lower quality of life in patients with breast cancer [15, 18, 19] and metastatic CRC [14, 20, 21]. However, longitudinal research on RAR in relation to patient-reported health outcomes among survivors of CRC is scarce. Therefore, we aimed to study longitudinal associations of RAR parameters with fatigue, insomnia, and HRQoL among survivors of stage I-III CRC, using five repeated measurements between 6 weeks and 5 years after CRC treatment. We hypothesized that less disrupted RAR, indicated by higher values for mesor, amplitude, dichotomy index, and 24-h autocorrelation coefficient, are associated with lower fatigue and insomnia, and a higher HRQoL after CRC treatment.

Methods

Study design and population

We used data of the Energy for Life after ColoRectal cancer (EnCoRe) study, an ongoing prospective cohort study of survivors of stage I-III CRC in the Netherlands (Netherlands Trial Register number: NL6904) [22]. From April 2012 onwards, all patients diagnosed with stage I-III CRC in three hospitals in the south of the Netherlands (Maastricht University Medical Center + , VieCuri Medical Center, and Zuyderland Medical Center) were eligible for inclusion. Exclusion criteria were stage IV CRC, age below 18 years, not able to understand and speak the Dutch language, a residential address outside of the Netherlands, or one or more comorbidities that could hinder successful participation (e.g., a cognitive disorder, or problems with visibility or hearing). Data were collected by trained dietitians visiting participants during 6 repeated home visits: one at diagnosis before the start of cancer treatment, and five after the end of cancer treatment (surgery, radiotherapy and/or chemotherapy) at 6 weeks, 6 months, 12 months, 24 months, and 60 months post-treatment. From October 2020 onwards, due to the COVID-19 pandemic, data were collected remotely via postal methods. All data collected between April 2012 and July 2018 were used for the present analyses. In addition, data collected until October 2021 were included for the 60-month post-treatment follow-up measurement. The Medical Ethics Committee of the University Hospital Maastricht and Maastricht University approved the study (METC 11–3-075). All patients provided written informed consent prior to participation, and the study was performed in accordance with the Declaration of Helsinki. A flow diagram describing recruitment and follow-up procedures of participants within the EnCoRe study can be found in Fig. 1.

Fig. 1
figure 1

Flow diagram of the inclusion of participants within the Energy for Life after ColoRectal cancer (EnCoRe) study from 2012 onwards and the number of post-treatment measurements available up to October 2021 included in the analyses of the current paper. 1Response rate = (persons with measurements)/(persons with measurements + persons lost to follow-up – persons died). The declining number of participants at subsequent time points is because not all participants included at diagnosis from April 2012 onwards had already reached all follow-up time points by July 2018 (6 weeks until 24 months post-treatment) or October 2021 (60 months post-treatment). 2Since the analyses in this paper were focused on diurnal rest-activity rhythms and fatigue, insomnia, and health-related quality of life (HRQoL) after colorectal cancer treatment, only post-treatment measurements with available data on diurnal rest-activity rhythms, fatigue, insomnia, or HRQoL, and covariates were included. 3Of the three persons without 6 weeks follow-up visits, one person did not have a 6 months follow-up visit before July 16th 2018. Of the six persons without 6 months follow-up visits, one person did not have a 12 months follow-up visit before July 16th 2018. 4All data collected from participants between April 2012 and July 2018 were used for the measurements from 6 weeks until 24 months post-treatment. In addition, these participants were followed over time and data was collected in October 2021 for the 60 months-treatment follow-up measurement

Diurnal rest-activity rhythm

Parameters of RAR were measured using the validated tri-axial MOX activity monitor (Maastricht Instruments B.V., Maastricht, The Netherlands) [23]. Participants wore the MOX on the anterior upper thigh ten cm above the knee for seven consecutive days (24 h/day) at every post-treatment time point, but not at diagnosis. The MOX measures raw acceleration data in three orthogonal sensor axes at a sampling rate of 25 Hz. These acceleration data were converted into activity counts in 1-min epochs (intervals), using the signal magnitude area [24]. The epochs were used for determining RAR parameters as described below. Monitor wear days with 24 h of wear time were considered valid (i.e. days with non-wear were excluded); only participants with ≥ 4 valid days including at least one weekend day were included in the analyses. Based on these criteria, 4.0% of accelerometer measurements were excluded.

RAR parameters included mesor, amplitude, acrophase, circadian quotient, dichotomy index, and 24-h autocorrelation. These measures were determined based on the activity counts at all post-treatment time points. These parameters were calculated using a custom-made MATLAB program (Maastricht Instruments B.V., Maastricht, The Netherlands), and are described in more detail below. Additionally, a schematic overview of these parameters is shown in Fig. 2. Each of the parameters, except for the 24-h autocorrelation, was calculated for every valid day, and afterwards averaged over all valid days available at each post-treatment time point. For the 24-h autocorrelation, one value for each post-treatment time point was obtained based on all days the accelerometer was worn at that specific time point. For obtaining the mesor, amplitude, and acrophase, the cosinor method, widely used within RAR studies, was used [25, 26].

Fig. 2
figure 2

Visualization of the operationalized diurnal rest-activity rhythm parameters: mesor, amplitude, acrophase, circadian quotient, dichotomy index, and 24-h autocorrelation coefficient, based on fictitious accelerometer data. A Example of a day the accelerometer was worn by a participant. Activity counts of the participant in each 1-min interval of a 24-h day were available. Based on how these activity counts were distributed across the 24-h day, the best cosine function was fitted according to the following formula: \((t)\hspace{0.17em}=\hspace{0.17em}M\hspace{0.17em}+\hspace{0.17em}Acos (2\uppi t/\mathrm{\rm T}\hspace{0.17em}+\hspace{0.17em}\upphi )\hspace{0.17em}+\hspace{0.17em}e(t)\). In this formula, M is the mesor, A is the amplitude, ϕ is the acrophase, Τ is the period, and e(t) is the error term. For each day the accelerometer was worn, a cosine function was fitted to determine the diurnal rest-activity parameters. Afterwards, these values were averaged across all available accelerometer days. B The dichotomy index was calculated as the percentage of in-bed (I) activity counts, expressed per 1-min interval, that were less than the median of out-of-bed (O) activity counts. The dichotomy index was calculated for each day the accelerometer was worn, and these values were averaged across all available days. C Overlap is shown between two different fictitious days the accelerometer was worn. The 24-h autocorrelation coefficient describes the average correlation between activity counts in 1-min intervals separated by 24 h. For the 24-h autocorrelation coefficient only one value was calculated, as the autocorrelation across all available accelerometer days. In this example, only day 1 and 2 are shown; however, activity counts separated by 24 h between day 2 and day 3, day 3 and 4 etc., are also included in the calculation

Mesor

The midline estimating statistic of the rhythm (mesor) is the 24-h rhythm-adjusted mean of the activity counts. Higher values for the mesor are indicative of increased activity over the 24-h day, whereas lower values indicate less activity [27].

Amplitude

The amplitude is the difference between the highest activity point of the cosinor curve (peak) and the mesor in activity counts. Higher values for amplitude indicate a larger contrast between average activity levels and peak activity [27].

Acrophase

The acrophase describes the clock time of the cosine weighted peak, and is expressed in decimal hours. Acrophase is often an appropriate estimate of the timing of a person’s 24h rhythm as it describes the timing of the cosine weighted peak [27].

Circadian quotient

The circadian quotient is calculated by dividing the amplitude by the mesor. The benefit of this measure is that, unlike unadjusted amplitude, it is correcting for average differences in activity. Higher values, i.e. higher amplitude relative to mesor, reflect a stronger RAR [16].

Dichotomy index

The dichotomy index (I < O) describes the percentage of in-bed (I) activity counts that are less than the median of out-of-bed (O) activity counts [28]. In this regard, wake and bed times were reported by participants in a structured 7-day sleep and dietary record. Lower values for I < O signify weaker RAR, whereas higher I < O values signal stronger RAR [28]. The I < O was calculated for each 24-h period ranging from 00:00 – 23:59h, on all valid wear days [28].

Twenty-four hour autocorrelation

The 24-h autocorrelation coefficient quantifies the regularity or consistency of the RAR from one day to the next, in other words between sequential days. More specifically, this coefficient describes the average correlation between activity (count) levels in 1-min epochs separated by 24 h. Higher and positive values for the 24-h autocorrelation coefficient indicate a more consistent RAR [16]. To have a representative number based on the full 7-day week at each post-treatment time point, the 24-h autocorrelation coefficient was only calculated in case of seven valid wear days. Due to non-wear on one or more days, 34% of accelerometer measurements were excluded for this parameter.

Patient-reported outcomes

Fatigue

The Checklist Individual Strength (CIS), a validated 20-item questionnaire, was used to acquire a multidimensional assessment of fatigue at all post-treatment time points [29]. The CIS consists of four subscales and a total score (ranging from 20 to 140) derived from summing the individual subscales. The subscales include subjective fatigue (ranging from 8 to 56), motivation- (ranging from 4 to 28), concentration- (ranging from 5 to 25), and activity-related fatigue (ranging from 3 to 21) [30]. Higher scores on all these scales indicate more fatigue. Participants’ fatigue was also measured using the fatigue symptom scale of the European Organization for the Research and Treatment of Cancer Quality of Life Questionnaire (EORTC QLQ-C30) [31]. The fatigue symptom scale consists of three items, and based on these items, a fatigue score was calculated which ranged from 0 to 100, with higher scores indicating more fatigue. The EORTC QLQ-C30 fatigue subscale is cancer-specific and has been validated in cancer populations [31].

Insomnia

The validated EORTC QLQ-C30 insomnia scale was used to assess the severity of sleep problems of participants at all post-treatment time points [31]. The insomnia scale consists of one item; “Have you had trouble sleeping in the past week?” with four response options from 1 “not at all” to 4 “very much”, based on which a score was calculated ranging from 0 to 100. Higher values indicate more sleep problems (insomnia). Recently, the one-item insomnia scale of the EORTC QLQ-C30 was found to be suitable for detecting sleep problems among patients with cancer and for testing associations of those sleep problems with other variables [32]. Additionally, sleep duration in hours was determined using self-reported wake and bed times by participants in a structured 7-day sleep and dietary record which was completed at all post-treatment time points on the same days as the accelerometer was worn. If missing, sleeping times were determined based on visual inspection of acceleration data (MOX activity monitor). Based on the reported sleeping times, the midpoint of sleep was calculated for every night of each participant. The average midpoint of sleep for each post-treatment time point was determined by averaging the midpoint of sleep across all available nights of every participant at that post-treatment time point.

Health-related quality of life

The widely used and well-validated EORTC QLQ-C30 was used to measure cancer-specific health-related quality of life (HRQoL) at all post-treatment time points [31]. The EORTC QLQ-C30 contains a global health/QoL scale, as well as subscales for physical, role, cognitive, emotional, and social functioning [31]. All scores were converted to a 0 to 100 scale, where higher scores indicated better QoL and functioning. In the current analyses, we only included global QoL and the physical functioning scale because these were hypothesized to be most likely related to RAR [14, 20].

Lifestyle, clinical, sociodemographic, and psychosocial factors

Clinical information, including cancer stage (I, II, III), type of treatment (including surgery, chemotherapy and/or radiotherapy), and tumor site (colon vs. rectum), were collected from medical records. Furthermore, self-reported data was retrieved on age, sex, educational level, employment status, smoking status, number of co-morbidities [33], and presence of a stoma. Trained dietitians performed anthropometric measurements at every post-treatment time point during home visits according to standardized procedures [34]. From October 2020 onwards, due to the COVID-19 pandemic, anthropometric data were self-reported. Participants’ body height and weight were measured in duplicate and then averaged to determine body mass index (BMI) in kg/m2. BMI was categorized using the World Health Organization (WHO) guidelines into underweight (BMI < 18.5 kg/m2), normal weight (18.5 ≤ BMI < 25 kg/m2), overweight (25 ≤ BMI < 30 kg/m2), or obesity (BMI ≥ 30 kg/m2). Dietary intake of participants, including alcohol consumption, was measured through 7-day food diaries collected at each post-treatment time point. Further, total time spent in total physical activity (all activities with an energy expenditure > 1.5 METs) was objectively determined at each post-treatment time point using the validated tri-axial MOX activity monitor. Similarly, prolonged sedentary behavior (both hours/day), i.e., time accumulated in sedentary bouts with a duration of at least 30 min, was determined at each post-treatment time point, as previously described in more detail [35]. Finally, the Hospital Anxiety and Depression Scale (HADS) was used to assess symptoms related to depression and anxiety at each post-treatment time point [36].

Statistical analyses

Descriptive analyses summarized characteristics of the study population, encompassing RAR parameters, fatigue, insomnia, HRQoL, as well as sociodemographic, psychosocial, lifestyle, and clinical variables at each post-treatment time point. Normally distributed quantitative variables were presented as mean (± SD), whereas non-normally distributed quantitative variables were presented as median (inter-quartile range, IQR). Categorical variables were described as frequency (percentage) across classes. For describing longitudinal changes over time in RAR parameters, fatigue, insomnia, and HRQoL, time was modeled as a categorical variable in linear mixed regression models, represented by dummy variables with the measurement at 6 months post-treatment selected as reference category (based on the observed pattern of changes over time). Furthermore, Pearson’s correlation coefficients were calculated to determine the correlation between RAR parameters, total physical activity, and prolonged sedentary time.

Longitudinal analyses were performed using linear mixed models, to analyze associations of RAR parameters with fatigue, insomnia, and HRQoL from 6 weeks to 60 months post-treatment. For comparability of results across parameters, all RAR parameters were standardized by dividing individual values by the mean standard deviation across all five post-treatment time points for each RAR parameter. Mesor, amplitude, circadian quotient, dichotomy index, and 24-h autocorrelation were modeled continuously. Acrophase was modeled categorically based on tertiles, as for this variable the linearity assumption was unrealistic. Longitudinal associations were adjusted for a priori-defined confounders, based on literature and causal reasoning. Fixed (time-invariant) confounders included age at enrollment (years), sex, neo-adjuvant radio- and/or chemotherapy (yes, no), adjuvant chemotherapy (yes, no), and education level (low, medium, high). Time-variant confounders, which were measured at all post-treatment time points, included number of comorbidities (0, 1, ≥ 2), BMI (kg/m2), stoma (yes, no), smoking (current, former, never), employment (employed, unemployed/retired), alcohol intake (g/day), and time since diagnosis (months). Associations for acrophase were additionally adjusted for the time-variant variable midpoint of sleep to account for potential differences in the timing of sleep. The use of random slopes was tested for all models with a likelihood-ratio test. These were added when the model improved statistically significantly. In addition, inter- and intra-individual associations were disaggregated by adding centered person-mean values to the model to estimate inter-individual associations (i.e., average differences between participants over time) and individual deviations from the person-mean value to estimate intra-individual associations (i.e., within-participant changes over time) [37].

To assess potential interaction between RAR parameters and sex, chemotherapy (yes/no), BMI (categorical at 6 weeks post-treatment), and time since diagnosis (months), in relation to fatigue, insomnia, and HRQoL, interaction terms were added into the linear mixed models. Statistical significance for the interaction term was set at P < 0.05.

As a sensitivity analysis, to obtain more insight into the possible direction of longitudinal associations, time-lag models were used, in which RAR parameters at earlier time points were coupled with fatigue, insomnia, and HRQoL outcomes at subsequent time points to mimic a more natural direction of associations. Two additional sensitivity analyses including additional adjustment for covariates were performed. First, we investigated associations of RAR with fatigue, insomnia, and HRQoL with additional adjustment for total physical activity and prolonged sedentary time. This was done to determine if observed associations could be (partly) attributed to a more (or less) active lifestyle in general, rather than 24h rhythms of activity and rest. Second, analyses were additionally adjusted for the total anxiety and depression score of participants at 6 weeks post-treatment, to investigate potential confounding effects. Finally, with regards to acrophase, we additionally investigated the phase difference as an exposure variable, as this represents an alternative method to look at RAR timing in relation to sleep timing [38]. The phase difference was calculated as the difference in clock-hours between the acrophase and the midpoint of sleep. All statistical analyses were performed using Stata 16.0 (StataCorp LLC) with statistical significance set at P < 0.05 (two-sided).

Results

Characteristics of the included study participants (n = 289) at 6 weeks post-treatment are presented in Table 1. About two-thirds of participants were men (67%), with an average age of 67 years (SD = 9) at enrollment. The majority of participants were colon cancer survivors (63%), whereas 37% were rectosigmoid or rectum cancer survivors. Approximately half of the participants were diagnosed with stage III CRC (44%), followed by stage I (31%), and stage II (25%). With regards to treatment, 89% of participants had surgery, 25% received radiotherapy (all except for one pre-operative), and 38% of participants received chemotherapy (18% of participants pre-operative and 29% post-operative, some received both). The average BMI of our population at 6 weeks post-treatment was 28 kg/m2 (SD = 5), and most participants (44%) had at least two co-morbidities. Median daily time spent in total physical activity (all activities with an energy expenditure > 1.5 METs) at 6 weeks post-treatment was 1.4 h (IQR: 1.0 – 1.9), whereas participants’ median daily time spent in prolonged sedentary behavior (in bouts with duration ≥ 30 min) was 4.8 h (IQR: 3.2 – 6.7) at 6 weeks post-treatment.

Table 1 Sociodemographic, lifestyle, and clinical characteristics of participants included in the Energy for Life after ColoRectal cancer (EnCoRe) study at 6 weeks after treatment

Descriptives of diurnal rest-activity rhythms, fatigue, insomnia, and HRQoL up to 60 months post CRC treatment

Descriptives of RAR parameters and fatigue, insomnia, and HRQoL from 6 weeks to 60 months post-treatment are presented in Table 2. At 6 weeks post-treatment, mesor, with higher values reflecting more activity during the 24-h day, was on average 3.71 (SD = 0.17). Similarly, amplitude, with higher values expressing a larger difference between peak- and mean activity levels, was 0.57 (SD = 0.18). Acrophase, indicating the clock-time of the peak activity, was 14:12h (IQR: 13:24h – 14:58h). Moreover, circadian quotient, which is the peak activity adjusted for mean activity with a higher value indicating a larger difference between peak and mean activity as the amplitude divided by the mesor, was 0.15 (SD = 0.04). Dichotomy index, with higher values reflecting a larger difference between daytime and nighttime activity, was 0.78 (SD = 0.21). Lastly, 24-h autocorrelation, with higher values indicating a more consistent rest-activity rhythm from one day to the next, was 0.16 (SD = 0.10). Linear mixed models showed that all RAR parameters, except for 24-h autocorrelation, increased statistically significantly from 6 weeks to 6 months post-treatment. The average absolute change in this time period was 0.03 (SD = 0.14), 0.03 (SD = 0.15), 0:01h (IQR: -0:26h – 0:39h), 0.01 (SD = 0.04), and 0.10 (SD = 0.23), for mesor, amplitude, acrophase, circadian quotient, and dichotomy index, respectively. Afterwards, from 6 to 60 months post-treatment, RAR parameters remained relatively stable, except for mesor where an increase at 60 months post-treatment was observed (0.05 (SD = 0.16)) (Table 2).

Table 2 Descriptive analyses of diurnal rest-activity rhythm parameters and fatigue, insomnia, and health-related quality of life in the study population of survivors of colorectal cancer from 6 weeks to 60 months post-treatment

All fatigue and insomnia scores significantly decreased, whereas physical functioning and global QoL significantly increased from 6 weeks to 6 months post-treatment. The decline in fatigue and insomnia scores, and the increase in physical functioning and global QoL appeared to continue up until 24 months post-treatment, although these changes were mostly non-significant (Table 2). Generally, moderate to strong positive correlations were observed among RAR parameters and with total physical activity, whereas moderate negative correlations were observed with sedentary behavior (Supplementary Fig. 1).

Longitudinal associations of diurnal rest-activity rhythm parameters with fatigue, insomnia, and HRQoL

Table 3 shows confounder-adjusted longitudinal associations of RAR parameters mesor, amplitude, acrophase, circadian quotient, dichotomy index, and 24-h autocorrelation with fatigue, insomnia, and HRQoL from 6 weeks to 60 months post-treatment, including overall, intra- and inter-individual associations.

Table 3 Longitudinal associations of diurnal rest-activity rhythm parameters with fatigue, insomnia, and health-related quality of life in the study population of survivors of colorectal cancer between 6 weeks and 60 months post-treatment

Fatigue

Higher values for mesor, amplitude, circadian quotient, and dichotomy index were statistically significantly associated with lower values for all fatigue outcomes (Table 3). For each SD increase in mesor (0.175), the CIS total score decreased on average by 5.3 points (95% CI: -6.9, -3.6), whereas each SD increase in amplitude (0.180) was associated with a decrease of 5.5 points (-7.1, -3.9) in the CIS total fatigue score. Similarly, each SD increase in circadian quotient (0.042) and dichotomy index (0.114) was significantly associated with decreases of 5.1 (-6.6, -3.6) and 3.1 points (-4.2, -2.1) in CIS total fatigue, respectively. Similar results were observed for the fatigue score as measured through the EORTC QLQ-C30 scale. Although for 24-h autocorrelation inverse associations were observed with all fatigue outcomes, only associations with subjective fatigue (CIS) (-1.3; -2.1, -0.4) and fatigue as measured by the EORTC subscale (-3.0; -4.6, -1.4) were statistically significant. No significant associations of the acrophase with fatigue were observed.

Most of the described associations with fatigue were driven by both intra-individual changes over time and inter-individual differences. However, generally, inter-individual associations were slightly stronger as compared to intra-individual associations. For example, inter-individual associations showed that participants whose average mesor was one SD (0.175) higher than other participants, reported lower levels of CIS total fatigue (β: -6.3; 95%CI: -9.1, -3.4), whereas for participants intra-individually increasing their mesor by one SD a smaller inverse association was found (-4.8; -6.8, -2.8).

Insomnia

Each SD increase in amplitude (0.180) and circadian quotient (0.042) was significantly associated with lower insomnia scores over time (-2.4; -4.3, -0.4 and -2.3; -4.2, -0.4, respectively). With regards to both amplitude and circadian quotient, this overall association was mainly driven by intra-individual changes over time (-2.9; -5.5, -0.2 and -2.5; -5.0, 0.1, respectively). For the RAR parameters mesor, acrophase, dichotomy index, and 24-h autocorrelation, no significant overall associations with insomnia were found (Table 3).

Physical functioning and global QoL

Higher values for mesor, amplitude, circadian quotient, dichotomy index, and 24-h autocorrelation were statistically significantly associated with better physical functioning and global QoL (Table 3). For each SD increase in mesor (0.175), physical functioning was on average 3.2 points (2.1, 4.4) higher and global QoL 3.7 points (2.5, 4.9) higher. Similarly, for each SD increase in amplitude (0.180), physical functioning was on average 4.7 points (3.7, 5.7) higher and global QoL 4.1 points (3.0, 5.2) higher. For each SD increase in circadian quotient (0.042), physical functioning was on average 4.6 points (3.6, 5.6) and global QoL 3.9 points (2.7, 5.0) higher. Similarly, for each SD increase in dichotomy index (0.114), physical functioning was on average 2.3 points (1.5, 3.1) higher and global QoL 1.7 points (0.7, 2.6) higher. Finally, each SD increase in 24-h autocorrelation (0.107) was associated with better physical functioning (4.1; 2.9, 5.4) and a higher global QoL (1.5; 0.1, 2.9). No significant associations of acrophase with either of these outcomes were observed.

Most of the described associations with physical functioning and global QoL appeared to be mainly driven by inter-individual differences over time, although also significant intra-individual associations were observed (Table 3). For example, inter-individual associations showed that participants whose average 24-h autocorrelation was one SD (0.107) higher than other participants, reported better physical functioning (5.7; 3.8, 7.7), whereas for participants intra-individually increasing their 24-h autocorrelation by one SD a smaller association was found (1.9; 0.5, 3.4).

Interaction and sensitivity analyses

For some of the described overall associations, significant interaction effects were observed with time since end of treatment. Therefore, in case of significant interaction, stratified effects for each post-treatment time point are visualized in Supplementary Fig. 2. In general, observed associations of RAR parameters with fatigue, insomnia, and HRQoL were strongest at the 6 weeks post-treatment time point, whereas for the RAR variable dichotomy index associations were strongest at 2–5 years post-treatment. No significant interaction effects were found for sex, chemotherapy, and BMI.

In the sensitivity analyses additionally adjusting for total physical activity and prolonged sedentary time, we found mostly similar associations as compared to the main results (Supplementary Fig. 3). Although slightly attenuated, most associations were in a similar direction and remained statistically significant, except for the associations of acrophase with fatigue outcomes, physical functioning, and global QoL, where associations were strengthened. Specifically, participants in higher tertiles, indicating a later clock-time for the cosine weighted peak, reported on average statistically significantly higher activity-related fatigue and lower physical functioning. Results from the sensitivity analysis additionally adjusting for anxiety and depression scores of participants at 6 weeks post-treatment were comparable to the results of the main analyses (Supplementary Fig. 4). In comparison with the results of the main analyses, the overall longitudinal associations of RAR parameters with fatigue, insomnia, and HRQoL as observed in the time-lag analyses were attenuated, and many betas were not statistically significant anymore (Supplementary Fig. 5). Finally, no significant associations were observed in the sensitivity analyses were we looked at the longitudinal association between the phase difference and fatigue, insomnia, and HRQoL (Supplementary Table 1).

Discussion

To our best knowledge, this study is the first to analyze how objectively assessed RAR are longitudinally associated with fatigue, insomnia, and HRQoL in survivors of CRC from 6 weeks up to 5 years post-treatment. In general, we observed that fatigue and insomnia decreased and HRQoL improved in survivors of CRC following treatment, suggesting recovery. However, some survivors still suffer from these complaints even after a prolonged period of time. Our main findings were that higher values for RAR parameters mesor, amplitude, circadian quotient, dichotomy index, and 24-h autocorrelation were significantly associated with lower fatigue and better physical functioning and global QoL over time. In addition, a higher amplitude and circadian quotient were significantly associated with less insomnia over time. For acrophase, no significant associations were observed with any of the outcomes. In general, the observed associations were in line with our hypothesis, highlighting that having a generally more active (mesor) and consistent (24-h autocorrelation) RAR, with a pronounced peak activity (amplitude) and a marked difference between daytime and nighttime activity (dichotomy index), was associated with less fatigue and a higher HRQoL in the five years after CRC treatment.

Our findings are in line with findings from prospective and cross-sectional studies in patients with metastatic CRC showing that a higher dichotomy index and 24-h autocorrelation coefficient were associated with lower fatigue and better physical functioning and global QoL [14, 20, 21]. Furthermore, we showed that also higher scores for other RAR parameters including mesor, amplitude, and circadian quotient were associated with lower fatigue and a better HRQoL. As far as we know, we are the first to report these associations in survivors of CRC up to 5 years post-treatment, although similar associations of mesor, amplitude, and circadian quotient with fatigue and HRQoL have been reported in breast cancer patients [15, 18, 19]. Our findings with regards to insomnia are in line with cross-sectional studies conducted in patients with breast and lung cancer, showing that a higher mesor, dichotomy index, and 24-h autocorrelation were associated with a better sleep quality, albeit not statistically significant in our study [39, 40]. In addition, we found statistically significant associations between a higher amplitude and circadian quotient and lower insomnia.

Previous research shows that survivors of CRC with more time spent in physical activity and less time spent sedentary reported on average lower levels of fatigue and a better HRQoL up to ten years post-diagnosis [11, 12, 41, 42]. This highlights that physical behavior could be an important lifestyle factor that may potentially empower patients with CRC to reduce their symptoms and subsequently improve their HRQoL. In the current study, we performed a sensitivity analysis where we additionally adjusted the associations between RAR parameters and fatigue, insomnia, and HRQoL for total physical activity and prolonged sedentary time, which showed high correlations with some of these parameters. The direction of observed associations was similar, with most associations remaining statistically significant. This finding suggests that independent of the total amount of physical activity or sedentary time itself, aspects related to RAR are also relevant for fatigue, insomnia, and HRQoL. A possible underlying mechanism that could explain the observed association may be that the timing of physical behavior could influence peripheral clocks which affect centrally (SCN) controlled circadian rhythms [9, 10]. If these circadian rhythms are disrupted or desynchronized, for example by not optimally timed physical activity, this may potentially influence fatigue, insomnia, and consequently HRQoL. In this regard, inflammation could play an important role as it is connected to both circadian rhythm disruption as well as patient-reported outcomes such as fatigue [43, 44]. However, also other common upstream influences affecting both RAR and patient-reported outcomes could play a role.

In our study, we observed significant intra-individual associations of mesor, amplitude, circadian quotient, dichotomy index, and 24-h autocorrelation with fatigue, physical functioning, and/or global QoL. These observed intra-individual associations indicate that an increase in these RAR parameters within participants after CRC treatment was significantly associated with decreased fatigue symptoms and/or increased HRQoL. Therefore, it would be interesting to investigate in an intervention study if patients with CRC with a tendency for a disrupted RAR can be assisted in restoring their RAR through targeting the relevant elements, and if this would also reduce their symptoms and increase their HRQoL. Most of our observed overall associations for the EORTC QLQ-C30 subscales appeared to be clinically relevant per 1 SD increase in RAR parameter, according to previously published thresholds to interpret observed changes [45, 46]. Largest intra-individual changes in RAR parameters were observed between 6 weeks and 6 months post-treatment, but these changes were still small compared to 1 SD. Nevertheless, the observed changes in RAR parameters still highlight the potential of interventions focusing on these parameters which may induce larger changes than we observed in our observational data. Ideally, these interventions should focus on the time-period shortly after CRC treatment, as we then observed the strongest associations (Supplementary Fig. 2).

A major strength of the current study is its prospective design with repeated and in-depth measurements of both RAR and patient-reported outcomes in a large group of survivors of CRC. Moreover, RAR of participants were objectively assessed using high-quality data of the 7-day worn tri-axial MOX activity monitor. Other strengths include the relatively high response rates at all post-treatment time points (up to 24 months > 90%, after 60 months post-treatment > 63%), small amounts of missing data as a result of intensive data collection, and the wide range of information collected on potential confounders. Finally, using linear mixed models enabled us to separate overall longitudinal associations into inter- and intra-individual associations, thereby providing valuable insights into the nature of observed relations over time. Despite the strengths of the current study, some limitations should be considered. This study is an observational study, which limits our ability to draw conclusions regarding causality, despite thorough adjustment for confounders in our analyses. Compared with results of the main analyses, overall longitudinal associations of RAR with fatigue, insomnia, and HRQoL were attenuated in the time-lag analyses, with many associations losing statistical significance. This indicates that observed associations are likely reciprocal. More specifically, RAR may have affected fatigue, insomnia, and HRQoL, but at the same time these outcomes could in turn also have altered RAR. Future randomized controlled trials should further investigate whether improving RAR can reduce fatigue and insomnia and improve HRQoL. An additional limitation could be the limited response rate at diagnosis to participate in the EnCoRe study (45%), which may have resulted in selection bias affecting the generalizability of our results. It is reasonable to think that individuals with worst RAR and most complaints did not participate, which may have led to an attenuation of our observed associations. In addition, because MOX activity monitor data was not collected before the start of treatment, we were unable to describe changes in activity habits and RAR parameters from pre-treatment to post-treatment. Finally, we cannot rule out the possibility of chance findings due to the high number of statistical tests performed. Nevertheless, the consistent significant and strength of longitudinal associations of RAR parameters with fatigue, insomnia, and HRQoL emphasize the importance of our findings.

Conclusions

Our findings show that in the first five years after CRC treatment, adhering to a generally more active and consistent RAR, with a pronounced peak activity and a marked difference between daytime and nighttime activity was found to be associated with lower fatigue, insomnia, and a better HRQoL among survivors of CRC. These findings provide valuable insights, which could help to fuel and design future intervention studies to investigate if restoring RAR has the potential to empower fatigued survivors of CRC to reduce fatigue and insomnia complaints and improve their HRQoL.

Availability of data and materials

Data described in the manuscript, code book, and analytic code will be made available upon request pending (e.g., application and approval, payment, other). Requests for data of the EnCoRe study can be sent to Dr. Martijn Bours, Department of Epidemiology, GROW School for Oncology and Reproduction, Maastricht University, the Netherlands (email: m.bours@maastrichtuniversity.nl).

Abbreviations

BMI:

Body mass index

CIS:

Checklist Individual Strength

CRC:

Colorectal cancer

EnCoRe:

Energy for Life after Colorectal Cancer

EORTC QLQ-C30:

European Organization for the Research and Treatment of Cancer Quality of Life Questionnaire

HADS:

Hospital Anxiety and Depression Scale

HRQoL:

Health-related quality of life

IQR:

Interquartile range

RAR:

Diurnal rest-activity rhythm

SCN:

Suprachiasmatic nucleus

SD:

Standard deviation

SQUASH:

Short Questionnaire to Assess Health-Enhancing Physical Activity

WCRF:

World Cancer Research Fund

WHO:

World Health Organization

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30.

    Article  PubMed  Google Scholar 

  3. Arndt V, Merx H, Stegmaier C, Ziegler H, Brenner H. Quality of life in patients with colorectal cancer 1 year after diagnosis compared with the general population: a population-based study. J Clin Oncol. 2004;22(23):4829–36.

    Article  PubMed  Google Scholar 

  4. O’Gorman C, Stack J, O’Ceilleachair A, Denieffe S, Gooney M, McKnight M, et al. Colorectal cancer survivors: an investigation of symptom burden and influencing factors. BMC Cancer. 2018;18:1–10.

    Article  Google Scholar 

  5. Howell D, Oliver T, Keller-Olaman S, Davidson J, Garland S, Samuels C, et al. Sleep disturbance in adults with cancer: a systematic review of evidence for best practices in assessment and management for clinical practice. Ann Oncol. 2014;25(4):791–800.

    Article  CAS  PubMed  Google Scholar 

  6. Innominato PF, Roche VP, Palesh OG, Ulusakarya A, Spiegel D, Lévi FA. The circadian timing system in clinical oncology. Ann Med. 2014;46(4):191–207.

    Article  CAS  PubMed  Google Scholar 

  7. Saligan LN, Olson K, Filler K, Larkin D, Cramp F, Sriram Y, et al. The biology of cancer-related fatigue: a review of the literature. Support Care Cancer. 2015;23:2461–78.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Koronowski KB, Sassone-Corsi P. Communicating clocks shape circadian homeostasis. Science. 2021;371(6530):eabd0951.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Lee Y. Roles of circadian clocks in cancer pathogenesis and treatment. Exp Mol Med. 2021;53(10):1529–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Hower IM, Harper SA, Buford TW. Circadian rhythms, exercise, and cardiovascular health. J Circadian Rhythms. 2018;16:7.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Kenkhuis MF, Van Roekel EH, Breedveld-Peters JJ, Breukink SO, Janssen-Heijnen ML, Keulen ET, et al. Longitudinal associations of sedentary behavior and physical activity with quality of life in colorectal cancer survivors. Med Sci Sports Exerc. 2021;53(11):2298.

    Article  PubMed  PubMed Central  Google Scholar 

  12. van Roekel EH, Duchâteau J, Bours M, Van Delden L, Breedveld-Peters J, Koole J, et al. Longitudinal associations of light-intensity physical activity with quality of life, functioning and fatigue after colorectal cancer. Qual Life Res. 2020;29:2987–98.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Van Blarigan EL, Meyerhardt JA. Role of physical activity and diet after colorectal cancer diagnosis. J Clin Oncol. 2015;33(16):1825.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Mormont M-C, Waterhouse J, Bleuzen P, Giacchetti S, Jami A, Bogdan A, et al. Marked 24-h rest/activity rhythms are associated with better quality of life, better response, and longer survival in patients with metastatic colorectal cancer and good performance status. Clin Cancer Res. 2000;6(8):3038–45.

    CAS  PubMed  Google Scholar 

  15. Berger AM, Wielgus K, Hertzog M, Fischer P, Farr L. Patterns of circadian activity rhythms and their relationships with fatigue and anxiety/depression in women treated with breast cancer adjuvant chemotherapy. Support Care Cancer. 2010;18:105–14.

    Article  PubMed  Google Scholar 

  16. Berger AM, Hertzog M, Geary CR, Fischer P, Farr L. Circadian rhythms, symptoms, physical functioning, and body mass index in breast cancer survivors. J Cancer Surviv. 2012;6:305–14.

    Article  PubMed  Google Scholar 

  17. Roveda E, Bruno E, Galasso L, Mulè A, Castelli L, Villarini A, et al. Rest-activity circadian rhythm in breast cancer survivors at 5 years after the primary diagnosis. Chronobiol Int. 2019;36(8):1156–65.

    Article  PubMed  Google Scholar 

  18. Sultan A, Choudhary V, Parganiha A. Worsening of rest-activity circadian rhythm and quality of life in female breast cancer patients along progression of chemotherapy cycles. Chronobiol Int. 2017;34(5):609–23.

    Article  PubMed  Google Scholar 

  19. Liu L, Rissling M, Neikrug A, Fiorentino L, Natarajan L, Faierman M, et al. Fatigue and circadian activity rhythms in breast cancer patients before and after chemotherapy: a controlled study. Fatigue. 2013;1(1–2):12–26.

    PubMed  PubMed Central  Google Scholar 

  20. Innominato PF, Focan C, Gorlia T, Moreau T, Garufi C, Waterhouse J, et al. Circadian rhythm in rest and activity: a biological correlate of quality of life and a predictor of survival in patients with metastatic colorectal cancer. Can Res. 2009;69(11):4700–7.

    Article  CAS  Google Scholar 

  21. Innominato PF, Komarzynski S, Palesh OG, Dallmann R, Bjarnason GA, Giacchetti S, et al. Circadian rest-activity rhythm as an objective biomarker of patient-reported outcomes in patients with advanced cancer. Cancer Med. 2018;7(9):4396–405.

    Article  PubMed  PubMed Central  Google Scholar 

  22. van Roekel EH, Bours MJ, de Brouwer CP, Ten Napel H, Sanduleanu S, Beets GL, et al. The applicability of the international classification of functioning, disability, and health to study lifestyle and quality of life of colorectal cancer survivors. Cancer Epidemiol Biomark Prev. 2014;23(7):1394–405.

    Article  Google Scholar 

  23. Annegarn J, Spruit MA, Uszko-Lencer NH, Vanbelle S, Savelberg HH, Schols AM, et al. Objective physical activity assessment in patients with chronic organ failure: a validation study of a new single-unit activity monitor. Arch Phys Med Rehabil. 2011;92(11):1852–7.e1.

    Article  PubMed  Google Scholar 

  24. Bijnens W, Aarts J, Stevens A, Ummels D, Meijer K. Optimization and validation of an adjustable activity classification algorithm for assessment of physical behavior in elderly. Sensors. 2019;19(24):5344.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Nelson W. Methods for cosinor-rhythmometry. Chronobiologia. 1979;6:305–23.

    CAS  PubMed  Google Scholar 

  26. Halberg F, Katinas GS. Chronobiologic glossary of the International Society for the Study of Biologic Rhythms. Int J Chronobiol. 1973;1(1):31–63.

    CAS  PubMed  Google Scholar 

  27. Lentz MJ. Time-series analysis-cosinor analysis: a special case. West J Nurs Res. 1990;12(3):408–12.

    Article  CAS  PubMed  Google Scholar 

  28. Natale V, Innominato PF, Boreggiani M, Tonetti L, Filardi M, Parganiha A, et al. The difference between in bed and out of bed activity as a behavioral marker of cancer patients: A comparative actigraphic study. Chronobiol Int. 2015;32(7):925–33.

    Article  PubMed  Google Scholar 

  29. Vercoulen JH, Swanink CM, Fennis JF, Galama JM, van der Meer JW, Bleijenberg G. Dimensional assessment of chronic fatigue syndrome. J Psychosom Res. 1994;38(5):383–92.

    Article  CAS  PubMed  Google Scholar 

  30. Servaes P, van der Werf S, Prins J, Verhagen S, Bleijenberg G. Fatigue in disease-free cancer patients compared with fatigue in patients with chronic fatigue syndrome. Support Care Cancer. 2001;9:11–7.

    Article  CAS  PubMed  Google Scholar 

  31. Aaronson NK, Ahmedzai S, Bergman B, Bullinger M, Cull A, Duez NJ, et al. The European Organization for Research and Treatment of Cancer QLQ-C30: a quality-of-life instrument for use in international clinical trials in oncology. J Natl Cancer Inst. 1993;85(5):365–76.

    Article  CAS  PubMed  Google Scholar 

  32. Hofmeister D, Schulte T, Hinz A. Sleep problems in cancer patients: a comparison between the Jenkins Sleep Scale and the single-item sleep scale of the EORTC QLQ-C30. Sleep Med. 2020;71:59–65.

    Article  PubMed  Google Scholar 

  33. Sangha O, Stucki G, Liang MH, Fossel AH, Katz JN. The Self-Administered Comorbidity Questionnaire: a new method to assess comorbidity for clinical and health services research. Arthritis Rheum. 2003;49(2):156–63.

    Article  PubMed  Google Scholar 

  34. Kenkhuis M-F, Van Roekel EH, Koole JL, Breedveld-Peters JJ, Breukink SO, Janssen-Heijnen ML, et al. Increases in adipose tissue and muscle function are longitudinally associated with better quality of life in colorectal cancer survivors. Sci Rep. 2021;11(1):12440.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Van Roekel EH, Winkler EA, Bours MJ, Lynch BM, Willems PJ, Meijer K, et al. Associations of sedentary time and patterns of sedentary time accumulation with health-related quality of life in colorectal cancer survivors. Prev Med Rep. 2016;4:262–9.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Zigmond AS, Snaith RP. The hospital anxiety and depression scale. Acta Psychiatr Scand. 1983;67(6):361–70.

    Article  CAS  PubMed  Google Scholar 

  37. Twisk JW, de Vente W. Hybrid models were found to be very elegant to disentangle longitudinal within-and between-subject relationships. J Clin Epidemiol. 2019;107:66–70.

    Article  PubMed  Google Scholar 

  38. Stone JE, Aubert XL, Maass H, Phillips AJ, Magee M, Howard ME, et al. Application of a limit-cycle oscillator model for prediction of circadian phase in rotating night shift workers. Sci Rep. 2019;9(1):1–12.

    Article  Google Scholar 

  39. Berger AM, Farr LA, Kuhn BR, Fischer P, Agrawal S. Values of sleep/wake, activity/rest, circadian rhythms, and fatigue prior to adjuvant breast cancer chemotherapy. J Pain Symptom Manage. 2007;33(4):398–409.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Chen H-M, Wu Y-C, Tsai C-M, Tzeng J-I, Lin C-C. Relationships of circadian rhythms and physical activity with objective sleep parameters in lung cancer patients. Cancer Nurs. 2015;38(3):215–23.

    Article  PubMed  Google Scholar 

  41. Lynch BM, van Roekel EH, Vallance JK. Physical activity and quality of life after colorectal cancer: overview of evidence and future directions. Expert Review of Quality of Life in Cancer Care. 2016;1(1):9–23.

    Article  Google Scholar 

  42. Eyl RE, Xie K, Koch-Gallenkamp L, Brenner H, Arndt V. Quality of life and physical activity in long-term (≥ 5 years post-diagnosis) colorectal cancer survivors-systematic review. Health Qual Life Outcomes. 2018;16:1–13.

    Article  Google Scholar 

  43. Wright KP Jr, Drake AL, Frey DJ, Fleshner M, Desouza CA, Gronfier C, et al. Influence of sleep deprivation and circadian misalignment on cortisol, inflammatory markers, and cytokine balance. Brain Behav Immun. 2015;47:24–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Bower JE, Lamkin DM. Inflammation and cancer-related fatigue: mechanisms, contributing factors, and treatment implications. Brain Behav Immun. 2013;30:S48–57.

    Article  CAS  PubMed  Google Scholar 

  45. Cocks K, King M, Velikova G, de Castro JG, St-James MM, Fayers P, et al. Evidence-based guidelines for interpreting change scores for the European Organisation for the Research and Treatment of Cancer Quality of Life Questionnaire Core 30. Eur J Cancer. 2012;48(11):1713–21.

    Article  CAS  PubMed  Google Scholar 

  46. Revicki D, Hays RD, Cella D, Sloan J. Recommended methods for determining responsiveness and minimally important differences for patient-reported outcomes. J Clin Epidemiol. 2008;61(2):102–9.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank all participants of the EnCoRe study and the health professionals in the three hospitals involved in the recruitment of participants of the study: Maastricht University Medical Centre+, VieCuri Medical Centre, and Zuyderland Medical Centre. We would also like to thank the MEMIC centre for data and information management for facilitating the logistic processes and data management of our study. Finally, we would like to thank the research dietitians and research assistant who are responsible for patient inclusion and follow-up, performing home visits, as well as data collection and processing.

Funding

Funding for SG_2021_076 was obtained from Wereld Kanker Onderzoek Fonds (WKOF), as part of the World Cancer Research Fund International grant programme. M.Y. Chong is supported by a grant from WKOF/WCRF: 2021/076. K.G. Frenken and E.H. van Roekel are supported by the Dutch Cancer Society (Grant No. 2021–1 / 13387). The EnCoRe study was also supported by the following other grants; Stichting Alpe d'Huzes within the research program “Leven met kanker” of the Dutch Cancer Society grants UM 2010–4867 and UM 2012–5653, by ERA-NET on Translational Cancer Research (TRANSCAN: Dutch Cancer Society (UM 2014–6877), by the WKOF/WCRF (WCRF: 2017/1619 and WCRF: 2016/1620), and by Kankeronderzoekfonds Limburg as part of Health Foundation Limburg grant 00005739. F.A.J.L. Scheer was supported in part by National Institutes of Health grants R01-HL140574 and R01-HL153969. All funders had no role in the design, analysis or writing of this article.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, MPW, MJLB, SJPME and EHvR; formal analysis, MYC and KGF; funding acquisition, MPW, MJLB, SJPME and EHvR; software, MPW; supervision, MPW, MJLB, SJPME and EHvR; writing–original draft, MYC and KGF; writing–review and editing, MPW, MJLB, SJPME, EHvR, AK, GKP, SOB, MJ, ETPK, WB, LMB, KM, FAJLS, KS, and JdV. All authors have read and agreed to the final version of the manuscript to be published.

Corresponding author

Correspondence to Marvin Y. Chong.

Ethics declarations

Ethics approval and consent to participate

This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects/patients were approved by the Medical Ethics Committee of the University Hospital Maastricht and Maastricht University (METC 11–3-075 | NL38786.068.11.). All patients provided written informed consent prior to participation.

Consent for publication

Not applicable.

Competing interests

Judith de Vos-Geelen has served as a consultant for Amgen, AstraZeneca, MSD, Pierre Fabre, and Servier, and has received institutional research funding from Servier, all outside the submitted work. F.A.J.L. Scheer served on the Board of Directors for the Sleep Research Society and has received consulting fees from the University of Alabama at Birmingham and Morehouse School of Medicine. F.A.J.L. Scheer interests were reviewed and managed by Brigham and Women’s Hospital and Partners HealthCare in accordance with their conflict of interest policies. F.A.J.L.S. consultancies are not related to the current work. The other authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Figure 1.

Figure comparing the Pearson correlation coefficients (n=1048) between diurnal rest-activity rhythms parameters, total physical activity, and prolonged sedentary time (time accumulated in sedentary bouts with a duration of at least 30 min, as measured by the MOX accelerometer) using data of all post-treatment time points. Abbreviations: CQ, circadian quotient; DI, dichotomy index; AC, 24-h autocorrelation; PA, physical activity; ST, sedentary time.

Additional file 2: Supplementary Figure 2.

Figure illustrating significant interaction effects of overall longitudinal associations between diurnal rest-activity rhythm parameters with fatigue, insomnia, and HRQoL with time since end of treatment. Stratified effects are shown separately for each post-treatment time point, and significantly different effects as compared to 6 weeks post-treatment are indicated with an asterisk. The P-value for trend indicates if the association between RAR parameters and outcomes differs significantly per 6 months increase in time since end of treatment. In case of no interaction, a horizontal line would be expected.

Additional file 3: Supplementary Figure 3.

Figure comparing the main results of the overall associations between diurnal rest-activity rhythms parameters and fatigue, insomnia, and HRQoL, and similar associations additionally adjusted for total physical activity and prolonged sedentary time (hours/day). Abbreviations: ME, mesor; AM, amplitude; AC2, acrophase tertile 2; AC3, acrophase tertile 3; CQ, circadian quotient; DI, dichotomy index; A24, 24-h autocorrelation.

Additional file 4: Supplementary Figure 4.

Figure comparing the main results of the overall associations between diurnal rest-activity rhythms parameters and fatigue, insomnia, and HRQoL, and similar associations additionally adjusted for anxiety and depression levels. Abbreviations: ME, mesor; AM, amplitude; AC2, acrophase tertile 2; AC3, acrophase tertile 3; CQ, circadian quotient; DI, dichotomy index; A24, 24-h autocorrelation.

Additional file 5: Supplementary Figure 5.

Figure comparing the main results of the overall associations between diurnal rest-activity rhythms parameters and fatigue, insomnia, and HRQoL, and similar associations from a time-lag model (exposure at one post-treatment time point combined with outcomes at the next post-treatment time point). Abbreviations: ME, mesor; AM, amplitude; AC2, acrophase tertile 2; AC3, acrophase tertile 3; CQ, circadian quotient; DI, dichotomy index; A24, 24-h autocorrelation.

Additional file 6.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chong, M.Y., Frenken, K.G., Eussen, S.J.P.M. et al. Longitudinal associations of diurnal rest-activity rhythms with fatigue, insomnia, and health-related quality of life in survivors of colorectal cancer up to 5 years post-treatment. Int J Behav Nutr Phys Act 21, 51 (2024). https://doi.org/10.1186/s12966-024-01601-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12966-024-01601-x

Keywords