Device-assessed sleep and physical activity in individuals recovering from a hospital admission for COVID-19: a multicentre study.

BACKGROUND: The number of individuals recovering from severe COVID-19 is increasing rapidly. However, little is known about physical behaviours that make up the 24-h cycle within these individuals. This study aimed to describe physical behaviours following hospital admission for COVID-19 at eight months post-discharge including associations with acute illness severity and ongoing symptoms.
METHODS: One thousand seventy-seven patients with COVID-19 discharged from hospital between March and November 2020 were recruited. Using a 14-day wear protocol, wrist-worn accelerometers were sent to participants after a five-month follow-up assessment. Acute illness severity was assessed by the WHO clinical progression scale, and the severity of ongoing symptoms was assessed using four previously reported data-driven clinical recovery clusters. Two existing control populations of office workers and individuals with type 2 diabetes were comparators.
RESULTS: Valid accelerometer data from 253 women and 462 men were included. Women engaged in a mean ± SD of 14.9 ± 14.7 min/day of moderate-to-vigorous physical activity (MVPA), with 12.1 ± 1.7 h/day spent inactive and 7.2 ± 1.1 h/day asleep. The values for men were 21.0 ± 22.3 and 12.6 ± 1.7 h /day and 6.9 ± 1.1 h/day, respectively. Over 60% of women and men did not have any days containing a 30-min bout of MVPA. Variability in sleep timing was approximately 2 h in men and women. More severe acute illness was associated with lower total activity and MVPA in recovery. The very severe recovery cluster was associated with fewer days/week containing continuous bouts of MVPA, longer total sleep time, and higher variability in sleep timing. Patients post-hospitalisation with COVID-19 had lower levels of physical activity, greater sleep variability, and lower sleep efficiency than a similarly aged cohort of office workers or those with type 2 diabetes.
CONCLUSIONS: Those recovering from a hospital admission for COVID-19 have low levels of physical activity and disrupted patterns of sleep several months after discharge. Our comparative cohorts indicate that the long-term impact of COVID-19 on physical behaviours is significant.

Background

There have been over 330 million confirmed cases of COVID-19 and over 5.5 million deaths [1]. Of the 15.6 million cases in the UK, there have been over 170,000 deaths and > 660,000 patients admitted to hospital [2]. As mortality improves, the number of post-hospitalisation survivors of COVID-19 is increasing. In some studies, more than 70% have not fully recovered by five months after discharge and have a substantial mental and physical health burden [3]. Given this, the pressing need for research to inform and support rehabilitation post-hospitalisation with COVID-19 is evident.

Emerging evidence suggests that physical activity [4, 5], good quality sleep, and regular sleep patterns [5] are associated with lower odds of being admitted to hospital or dying with COVID-19. This may occur through a reduction in chronic inflammation [6, 7] and lower cardiometabolic risk factors, features associated with an increased risk of COVID-19 [8], and/or through enhanced immunity [9]. The ongoing burden of symptoms associated with poor recovery [3] may have a detrimental impact on physical activity and sleep behaviours in post-hospitalisation survivors of COVID-19.

The post-hospitalisation COVID-19 (PHOSP-COVID) study is a large prospective multicentre follow-up study with the aim of understanding and improving long-term health outcomes following COVID-19 (https://phosp.org). Cluster analysis identified four recovery phenotypes relating to the severity of physical, mental, and cognitive health impairments an average of five months post-hospitalisation with COVID-19 [3].

The aim of this study was to describe accelerometer-assessed physical behaviours in patients post-hospitalisation with COVID-19 and to understand whether there are differences in physical behaviours associated with acute illness severity or the four recovery clusters. Physical behaviours within the PHOSP-COVID cohort were also described relative to a cohort of office workers and a cohort of adults with type 2 diabetes.

Methods

Participants and methods

PHOSP-COVID is a prospective longitudinal cohort study recruiting patients aged ≥ 18 years who were discharged from 83 National Health Service (NHS) hospitals across England, Northern Ireland, Scotland, and Wales following admission to a medical assessment or ward for confirmed or clinician-diagnosed COVID-19. Participants were excluded if they: i) had a confirmed diagnosis of a pathogen unrelated to the objectives of this study, ii) attended an accident and emergency department but were not admitted, iii) had another life-limiting illness with life expectancy less than six months. Clinical data collected as part of standard care was obtained from all patients discharged from participating hospitals following participant consent. An additional research visit was offered to all patients from discharge at 3 months (range 2–7 months). The 1077 participants who attended a dedicated research visit within one-year post-discharge were eligible for inclusion [3]. The included participants were assessed at a median five-months (range 2–7 months) between March and November 2020. Recruitment of participants continued until April 2022.

All study participants provided written informed consent. The study was approved by the Leeds West Research Ethics Committee (20/YH/0225) and is registered on the ISRCTN Registry (ISRCTN10980107).

Accelerometer data collection

Physical behaviours (i.e., physical activity and sleep) were assessed using the GENEActiv accelerometer (GENEActiv Original, ActivInsights, Kimbolton, UK). The monitors were initialised to record triaxial accelerations for 21 days at 30 Hz, with participants being asked to wear the monitor for 14 days.

Where possible, participants received the monitor and instructions by post within one month of their PHOSP-COVID research visit (2–7 months from discharge). Participants were instructed to start wearing the monitor on their non-dominant wrist immediately upon receiving it and to wear it 24 h/day. Participants were asked to return their monitors in a prepaid envelope after the 14-day assessment period.

Accelerometer data processing

Accelerometer files were processed with R-package GGIR version 2.2–0 (http://cran.r-project.org) [10]. Participants were excluded if they had < 3 days of valid wear (defined as ≥ 16 h/day). Signal processing in GGIR includes autocalibration using local gravity as a reference; detection of non-wear; calculation of the average magnitude of dynamic acceleration corrected for gravity (Euclidean Norm minus 1 g with negative values rounded up to zero, ENMO), averaged over 5-s epochs and expressed in milli-gravitational units (mg). Non-wear was imputed using the default setting, that is, invalid data were imputed by the average at similar time-points on different days of the week. Participants were excluded if post-calibration error was > 0.01 g (10 mg), they had < 3 days of valid wear (defined as ≥ 16 h per day), or if wear data were not present for each 15-min period of the 24-h cycle. Definitions of the physical behaviour variables are shown in the Supplementary material.

Physical behaviour characteristics included average acceleration in mg (a proxy for physical activity volume), moderate-to-vigorous intensity physical activity (MVPA) in minutes accumulated in ≥ 1-min bouts (> 100 mg) [11], light-intensity activity in minutes (between 40 and 100 mg) [12], inactive time (a proxy for sedentary time) in minutes (< 40 mg) [12], and intensity of the most active continuous 30 and 10 min/day (mg).

The percentage of participants meeting the guidelines of 150 min per week of MVPA [13] was calculated based on their average daily MVPA (i.e., if > 21.4 min/day) and used in analyses as a binary variable.

Weekly physical activity characteristics included the number of days/week with 10- and 30-min continuous MVPA. Weekly variables were limited to participants with at least seven valid days of data.

Sleep characteristics included total sleep time (hours, time between sleep onset and wake time minus any time awake), sleep efficiency (%), and mid-sleep variability (within-person standard deviation of mid-sleep time). Sleep mid-point variability describes how variable people are in the timing of their sleep.

Disease exposures

Acute illness severity was determined by the WHO clinical progression scale [14], defined as: i) 3–4 = no continuous supplemental oxygen needed, ii) 5 = continuous supplemental oxygen only, iii) 6 = continuous positive airway pressure ventilation (CPAP), bi-level positive airway pressure (BIPAP) or high flow nasal oxygen, iv) 7–9 = invasive mechanical ventilation (IMV) or extra-corporeal membrane oxygenation (ECMO).

The severity of ongoing symptoms after discharge was categorised on clusters derived previously within PHOSP-COVID where unsupervised machine learning using data from a battery of patient-reported outcomes and physical tests, identified four recovery clusters described as: 1) Very severe mental and physical health impairment, 2) Severe mental and physical health impairment, 3) Moderate mental and physical health impairment with pronounced cognitive impairment, 4) Mild mental and physical health impairment [3]. This outcome is referred to as recovery clusters.

Data for the severity of ongoing symptoms were missing in 27% of included participants due to missing data on the patient-reported outcomes and physical tests used for the cluster analysis [3].

Covariates

Data on sex, age at admission, ethnicity, number of chronic diseases, body mass index (BMI), and deprivation were included in this study (Supplementary material).

Comparative cohorts

In response to the lack of baseline data for the PHOSP-COVID cohort, the cohort was compared to accelerometer data collected in a cohort of office workers [15] and individuals with type 2 diabetes [16, 17] for descriptive purposes. Details of the comparative cohorts are described in the Supplementary material. Accelerometer data from all three cohorts were processed using identical methods.

Data selection

For this analysis, only those with valid accelerometer data were included. Of the 1077 participants in the PHOSP-COVID dataset, postal addresses were available to the research team for 853 participants within one month of follow-up clinical and research data collection visits allowing accelerometers to be posted, of which 796 were returned, with 715 providing valid data (Supplementary Figure S1).

Statistical analysis

Differences in physical activity and sleep across acute illness severity and recovery clusters were assessed using generalised linear models. Continuous variables were analysed using a normal distribution with an identity link. Model selection was informed by the Akaike Information Criterion. Although some physical activity variables displayed non-parametric distributions, adjusted model fit was not meaningfully improved using different distribution or log links once covariates were added. Physical activity bout data was analysed using a Poisson distribution with a log-linear link as count data. Binary logistic regression was used to investigate the odds of not meeting 150 min of MVPA per week across acute illness severity and recovery clusters and reported as odds ratios (95% CI).

Data were adjusted for age at admission, sex, ethnicity, deprivation, number of comorbidities, season of data collection, number of wear days (activity outcomes) or wear nights (sleep outcomes). Interactions between sex and acute illness severity/recovery cluster were included to determine whether differences in physical activity or sleep variables across acute illness severity or recovery clusters varied by sex. Data are reported as sex-stratified marginal means (95% CI) derived from this model.

Generalised linear models were also used to examine associations between the variables that made up the ongoing severity cluster definitions (breathlessness, fatigue, anxiety, depression, post-traumatic stress disorder, physical performance, and cognition—Supplementary material) and physical activity and sleep characteristics. Variables were standardised and analysed as continuous variables. After generating the main effect for each exposure, cluster variable by sex interactions were added to the models and significant interactions were further stratified by sex. Data are reported as beta-coefficients (95% CI).

A sensitivity analysis was conducted, removing healthcare workers to examine whether healthcare work status had an impact on sleep variables.

Data were analysed using SPSS (version 26.0). A p-value of < 0.05 was considered statistically significant.

Results

Of the 1077 participants included in the PHOSP-COVID dataset, 715 (253 women, 462 men) had valid accelerometer and acute illness severity data. Of these, 521 (185 women, 336 men) also had data for the severity of ongoing symptoms. Participant characteristics are displayed in Table 1. 151 (32.7%) men and 54 (21.3%) women received invasive mechanical ventilation (IMV) during the acute illness (WHO class 7–9), and 172 (37.2%) men and 66 (26.1%) women were classified within the very severe recovery cluster. Participants’ characteristics with accelerometer data compared to those without are displayed in Supplementary Table S1. The proportion of the population within the different classifications and clusters of disease severity were similar in those with complete and missing accelerometer data. However, those with complete data were older (59 vs. 55 years), with a higher proportion from White ethnicities (69.8% vs. 58.3%) and the least deprived quintile (20.3% vs. 14.6%).

Table 1

Participant characteristics

		Women (n = 253)		Men (n = 462)
Categorical variables		Count	Column %	Count	Column%
WHO disease severity class (acute COVID-19 severity)a	Class 3–4	72	28.5%	67	14.5%
	Class 5	89	35.2%	166	35.9%
	Class 6	38	15.0%	78	16.9%
	Class 7–9	54	21.3%	151	32.7%
Recovery cluster	Cluster 4: Mild	38	15.0%	48	10.4%
	Cluster 3: Moderate	50	19.8%	58	12.6%
	Cluster 2: Severe	31	12.3%	58	12.6%
	Cluster 1: Very Severe	66	26.1%	172	37.2%
	Missing	68	26.9%	126	27.3%
Comorbidities	No comorbidity	71	28.1%	134	29.0%
	1 comorbidity	53	20.9%	92	19.9%
	2 + comorbidities	129	51.0%	236	51.1%
Ethnicity	White	175	69.2%	324	70.4%
	South Asian	32	12.6%	60	13.0%
	Black	28	11.1%	21	4.6%
	Other	11	4.3%	32	7.0%
	Missing	7	2.8%	23	5.0%
IMD (quintile)b	1 (most deprived)	37	14.7%	98	21.2%
	2	59	23.4%	107	23.2%
	3	55	21.8%	72	15.6%
	4	46	18.3%	95	20.6%
	5 (least deprived)	55	21.8%	90	19.5%
Number of days where a 10-min moderate-intensity bout of physical activity was undertaken	0	139	55.8%	199	44.2%
	1	48	19.3%	82	18.2%
	2	28	11.2%	58	12.9%
	3 +	34	13.7%	111	24.7%
Number of days where a 30-min bout of moderate-intensity physical activity was undertaken	0	191	76.7%	285	63.3%
	1	33	13.3%	68	15.1%
	2	14	5.6%	45	10.0%
	3 +	11	4.4%	52	11.6%
Participants meeting > 150 min of moderate-to-vigorous intensity physical activity per week		68	26.9%	172	37.2%
Continuous variables		Mean	SD	Mean	SD
Age (Years)		58	14	60	12
BMI (kg/m2)c		32.2	7.9	30.4	6.1
Physical activity volume (mg)		19.3	6.6	19.6	7.2
Time spent in moderate or vigorous intensity physical activity (minutes/day)		14.9	14.8	21.1	22.3
Time spent in light intensity physical activity (minutes/day)		149.8	54.9	138.8	51.4
Time spent inactive (hours/day)		12.1	1.7	12.6	1.7
Intensity of the most active 10 min (mg)		61.0	39.0	73.6	64.0
Intensity of the most active 30 min (mg)		40.4	25.4	49.6	47.8
Total sleep time (hours/day)		7.2	1.1	6.9	1.1
Sleep Efficiency (%)		85.7	5.2	85.0	5.7
Variability in sleep midpoint (SD in minutes)		120.5	89.5	113.4	84.9
Number of valid days		14	2	14	2
Number of valid nights		13	3	13	3

a WHO clinical progression scale: 3–4 = no continuous supplemental oxygen needed, 5 = continuous supplemental oxygen only, 6 = continuous positive airway pressure ventilation (CPAP), bi-level positive airway pressure (BIPAP) or high flow nasal oxygen, 7–9 = invasive mechanical ventilation (IMV) or extra-corporeal membrane oxygenation (ECMO)

b IMD = Index of Multiple Deprivation

c BMI = body mass index

The summary variables from the accelerometer data for women and men are displayed in Table 1. The median time from discharge to accelerometer wear was 245 days [IQR 178–276 days]. The median time from the PHOSP-COVID research visit to accelerometer wear was 65 days [IQR 11–93 days]. Accelerometer data were available for a mean of 14 valid days. Women engaged in a mean ± SD of 14.9 ± 14.7 min/day of MVPA, with 12.1 ± 1.7 h/day spent inactive and 7.2 ± 1.1 h/day asleep. The same values for men were 21.0 ± 22.3 and 12.6 ± 1.7 h/day and 6.9 ± 1.1 h/day, respectively. Variability in sleep midpoint was ~ 2 h in men and women. Over 60% of both women and men did not have any days in a week that contained a 30-min bout of MVPA, e.g., walking, with most women (56%) also not having any days with a bout of 10-min of MVPA (Table 1).

Associations with disease severity

Across acute illness severity, those who had the most severe acute illness had ~ 1–2 mg lower volume of physical activity (p = 0.045) and less time spent in MVPA (p = 0.032) (Table 2). Women who received IMV undertook the lowest levels of MVPA [13.7 min/day; 95% CI 7.3, 20.2] (Table 2). Women and men with the most severe acute disease were 3.38 (95% CI 1.29, 8.85) and 2.17 (95% CI 1.06, 4.45) times more likely, respectively, to not meet physical activity recommendations for health compared to those with the least severe disease (Supplementary Figure S2). The pattern of number of days/week with continuous bouts of MVPA was similar across acute illness severity (Fig. 1).

Table 2

Physical activity and sleep variables across WHO classes of acute illness severity

	Women				Men
Physical activity variables	Class 3–4	Class 5	Class 6	Class 7–9	Class 3–4	Class 5	Class 6	Class 7–9	P for class	P for class x sex
Physical activity volume (mg)	21.3 (19.1, 23.3)	22.3 (20.4, 24.23)	19.9 (17.4, 22.4)	19.8 (17.6, 22.0)	21.9 (19.8, 24.0)	22.6 (21.0, 24.2)	23.3 (21.2, 25.4)	21.0 (19.3, 22.8)	0.045	0.248
Moderate to vigorous intensity physical activity (minutes/day)	18.7 (12.9, 24.6)	19.0 (13.6, 24.5)	16.0 (8.7, 23.2)	13.7 (7.3, 20.2)	25.1(19.1, 31.1)	28.5 (23.8, 33.2)	27.0 (21.1, 33.0)	21.6 (16.6, 26.7)	0.032	0.783
Light intensity physical activity (minutes/day)	167.8 (152.0, 183.5)	171.0 (156.2, 185.7)	150.2 (130.7, 169.6)	158.2 (140.7, 175.7)	154.0 (137.8, 170.2)	153.8 (141.3, 166.3)	160.7 (144.6, 176.9)	149.1 (135.5, 162.8)	0.386	0.150
Inactivity (hours/day)	12.0 (11.5, 12.5)	12.2 (11.8, 12.7)	12.3 (11.7, 12.9)	12.5 (12.0, 13.0)	12.7 (12.2, 13.3)	12.5 (12.1, 12.9)	12.5 (12.0, 13.0)	13.0 (12.5, 13.4)	0.201	0.569
Intensity of the most active continuous 30 min (mg)	54.6 (42.1, 67.2)	55.9 (44.2, 67.5)	47.8 (32.2, 63.3)	51.7 (38.0, 65.3)	62.8 (50.0, 75.7)	66.5 (56.5, 76.5)	62.9 (50.2, 75.7)	63.1(52.2, 73.9)	0.659	0.934
Intensity of the most active continuous 10 min (mg)	76.8 (59.6, 93.9)	85.3 (69.3, 101.2)	71.4 (50.1, 92.6)	74.7 (56.0, 93.5)	94.7 (77.1, 112.3)	98.8 (85.2, 112.5)	96.6 (79.1, 114.1)	88.5 (73.6, 103.4)	0.334	0.831
Sleep variables
Total sleep time (hours/day)	7.0 (6.6, 7.3)	6.8 (6.5, 7.2)	7.0 (6.6, 7.4)	6.7 (6.4, 7.1)	6.6 (6.2, 6.9)	6.5 (6.3, 6.8)	6.5 (6.2, 6.8)	6.5 (6.2, 6.8)	0.514	0.771
Sleep efficiency (%)	86.8 (85.1, 88.4)	86.1 (84.6, 87.7)	85.9 (83.8, 88.0)	85.8 (84.0, 87.6)	85.9 (84.3, 87.7)	85.4 (84.0, 86.7)	84.8 (83.1, 86.5)	85.5 (84.0, 86.9)	0.524	0.951
Sleep midpoint variability (minutes)	91.8 (66.1, 117.6)	92.7 (68.7, 116.8)	88.1 (56.4, 120.0)	113.7 (85.4, 141.9)	95.6 (69.1, 122.0)	91.9 (71.5, 112.3)	106.9 (80.6, 133.2)	85.7 (63.5, 108.0)	0.854	0.151

Data reported as marginal mean (95% CI). Adjusted for age, sex, ethnicity, deprivation, number of comorbidities, season of data collection and number of wear days (physical activity variables) or wear nights (sleep variables). Models included 675 participants with complete cluster and covariate data

Fig. 1

Proportion of participants within each WHO class of COVID-19 severity undertaking continuous bouts of physical activity. Data display the proportion within each WHO class achieving 0, 1, 2 and 3–7 days per week with a bout of 10 min (top panel) or 30 min (bottom) at least moderate-intensity physical activity. Adjusted for age, sex, ethnicity, deprivation, number of comorbidities, season of data collection, and number of wear days. 10-min bouts: p for difference by cluster = 0.132, p for sex x cluster = 0.395. 30-min bouts: p for difference by cluster = 0.217, p for sex x cluster = 0.128

Across recovery clusters, there was no difference in daily volume of physical activity and time spent in light activity or MVPA (Table 3). However, time spent inactive was greater in men than women (p = 0.013), with men in the very severe recovery cluster spending the most time inactive [13.2 h/day; 95% CI 12.6, 13.7]. Men in the severe and the very severe recovery clusters were also 2.52 (95% CI 1.19, 5.36) and 3.48 (95% CI 1.41, 8.59) times more likely, respectively, to not meet physical activity recommendations for health compared to those in the less severe recovery clusters (Supplementary Figure S3).

Table 3

Physical activity and sleep variables across four recovery clusters

	Women				Men
Physical activity variables	Cluster 4: Mild	Cluster 3: Moderate	Cluster 2: Severe	Cluster 1: Very Severe	Cluster 4: Mild	Cluster 3: Moderate	Cluster 2: Severe	Cluster 1: Very Severe	P for class	P for class x sex
Physical activity volume (mg)	21.5 (19.5, 23.6)	22.2 (19.5, 24.9)	19.3 (16.9, 21.8)	20.6 (18.1, 23.1)	22.3 (20.6, 24.1)	22.3 (20.2, 24.5)	22.7 (20.5, 24.9)	21.1 (18.6, 23.5)	0.395	0.293
Moderate to vigorous intensity physical activity (minutes/day)	19.2 (13.0, 25.4)	18. (10.3, 26.4)	16.0 (8.8, 23.3)	16.8 (9.2, 24.3)	29.3 (24.1, 34.5)	26.7 (20.3, 33.1)	24.0 (17.3, 30.6)	21.2 (13.8, 28.6)	0.156	0.756
Light intensity physical activity (minutes/day)	164.8 (149.6, 180.1)	168.2 (148.4, 187.9)	145.2 (127.4, 163.0)	156.3 (137.9, 174.7)	148.1 (135.4, 160.8)	153.2 (137.5, 168.9)	158.3 (141.9, 174.8)	139.3 (121.2, 157.3)	0.357	0.065
Inactivity (hours/day)	12.2 (11.7, 12.7)	11.9 (11.3, 12.6)	12.5 (12.0, 13.1)	12.5 (11.9, 13.1)	12.8 (12.4, 13.2)	12.8 (12.2, 13.3)	12.1 (11.5, 12.6)	13.2 (12.6, 13.7)	0.131	0.013
Intensity of the most active continuous 30 min (mg)	56.2 (41.9, 70.5)	59.6 (41.2, 78.1)	53.0 (36.3, 69.8)	45.2 (27.9, 62.4)	70.4 (58.4, 82.4)	69.29 (54.5, 84.1)	69.4 (53.8, 85.0)	55.9 (38.6, 73.2)	0.170	0.954
Intensity of the most active continuous 10 min (mg)	86.1 (66.6, 105.6)	84.6 (59.4, 109.7)	78.2 (55.3, 101.0)	67.0 (43.5, 90.56)	104.6 (88.3, 121.0)	97.43 (77.3, 117.6)	95.6 (74.3, 116.9)	80.5 (57.0, 104.1)	0.079	0.982
Sleep variables
Total sleep time (hours/day)	7.0 (6.6, 7.3)	6.8 (6.4, 7.2)	7.2 (6.8, 7.6)	6.8 (6.4, 7.3)	6.5 (6.2, 6.8)	6.4 (6.1, 6.8)	6.9 (6.5, 7.2)	6.5 (6.1, 6.9)	0.039	0.960
Sleep efficiency (%)	85.9 (84.2, 87.6)	85.1(82.9, 87.3)	87.3 (85.2, 89.5)	86.7 (84.6, 88.7)	85.6 (84.2, 87.1)	84.6 (82.8, 86.3)	85.2 (83.4, 87.1)	85.1 (83.1, 87.2)	0.412	0.534
Sleep midpoint variability (minutes)	92.4 (66.3, 87.6)	95.3 (61.1, 87.3)	74.5 (43.7, 89.3)	112.0 (80.3, 88.7)	87.4 (65.6, 87.1)	97.8 (70.7, 86.3)	78.9 (50.4, 87.1)	105.5 (74.3, 87.2)	0.070	0.956

Data reported as marginal mean (95% CI). Adjusted for age, sex, ethnicity, deprivation, number of comorbidities, season of data collection and number of wear days (physical activity variables) or wear nights (sleep variables). Models included 490 participants with complete cluster and covariate data

There was also a notable difference in the number of days/week on which longer bouts of physical activity were undertaken. In the very severe recovery cluster, over 80% of women and men did not undertake a bout of MVPA lasting 30 min on any day of the week, with over 60% not undertaking a bout lasting 10 min on any day (Fig. 2). The frequency of longer bouts of MVPA was substantially higher in those within the mild recovery cluster; over 20% in the mild cluster undertook at least a 10-min bout of MVPA on at least 3 days/week (Fig. 2). Men and women in the severe recovery cluster had the longest total sleep time (p = 0.039), with the very severe recovery cluster having the greatest sleep midpoint variability, although this did not reach significance (p = 0.070) (Table 3).

Fig. 2

Proportion of participants within each recovery cluster undertaking continuous bouts of physical activity. Data display the proportion within each cluster achieving 0, 1, 2 and 3–7 days per week with a bout of 10 min (top panel) or 30 min (bottom) at least moderate-intensity physical activity. Adjusted for age, sex, ethnicity, deprivation, number of comorbidities, season of data collection, and number of wear days. 10-min bouts: p for difference by cluster < 0.001, p for sex x cluster = 0.672. 30-min bouts: p for difference by cluster = 0.021, p for sex x cluster = 0.262

Sleep variables across acute illness severity and recovery clusters showed a similar pattern when removing healthcare workers (N = 98) (Supplementary Table S2).

The associations of recovery cluster variables with physical activity and sleep characteristics are shown in Fig. 3 (data shown in Supplementary Table S3). Lower severity of symptoms, except for cognition and anxiety, were positively associated with physical activity (p < 0.05). More severe breathlessness, fatigue, anxiety, depression, and post-traumatic stress disorder were all associated with greater sleep midpoint variability (p < 0.01). More severe depression was also associated with lower sleep efficiency (p = 0.039). Associations of physical performance with MVPA and intensity of the most active continuous 30/10 min were stronger in men than women (p for interaction < 0.05) (Supplementary Table S3 and S4).

Fig. 3

Forest plots of associations between patient-reported outcomes with physical activity and sleep characteristics. Data are shown as beta-coefficients (95% CI) showing the difference per SD in the exposure. Adjusted for age, sex, ethnicity, deprivation, number of comorbidities, season of data collection, and number of wear days (physical activity variables) or nights (sleep variables). FACIT = Functional Assessment of Chronic Illness Therapy; SPPB = short physical performance battery; MoCA = Montreal Cognitive Assessment; GAD-7 = General Anxiety Disorder 7 Questionnaire; PHQ-9 = Patient Health Questionnaire – 9; PTSD = post-traumatic stress disorder; PCL-5 = Post Traumatic Stress Disorder Checklist

Comparison to cohorts of office workers and those with type 2 diabetes

Characteristics of the SMART Work and Life (SWL) (n = 232) and Chronotype of Patients with Type 2 Diabetes and Effect on Glycaemic Control (CODEC) (n = 685) comparator cohorts are shown in Supplementary Table S5 with differences in physical behaviours shown in Supplementary Table S6. Mean age was within five years of the PHOSP-COVID cohort for both comparator cohorts, while the CODEC cohort was well-matched for key characteristics including sex, multimorbidity status and BMI. Overall, activity was higher in the SWL and CODEC cohorts compared to PHOSP-COVID, with notably higher activity in the SWL cohort. The differences in activity volume of 1.1 mg in CODEC and 3.2 mg in SWL. Those in PHOSP-COVID spent ~ 17 fewer minutes in light-intensity activity than those in CODEC (p = 0.004) and ~ 16 fewer minutes in MVPA (p < 0.001) than those in SWL, with men also spending more time inactive compared to CODEC (p = 0.020). The intensity of the most active 30/10 min was also lowest in PHOSP-COVID (p < 0.001). The frequency of continuous bouts of MVPA/week in PHOSP-COVID was similar to CODEC, but notably lower compared to SWL (p < 0.001) (Supplementary Figure S4). The variability in sleep midpoint was at least three times greater in PHOSP-COVID compared to SWL and CODEC (p < 0.001), with sleep efficiency also lower (p < 0.001) being 83.5% [95%CI 82.3, 84.7] in women and 82.7% [95%CI 81.6, 83.9] in men, but ≥ 86.5% in SWL and CODEC. Total sleep time was similar across the three cohorts (p > 0.092).

Discussion

Women and men recovering from a hospital admission for COVID-19 had low levels of physical activity and high variability in sleep timing. More severe acute illness was associated with approximately 1–2 mg lower volume of physical activity (approximating 500–1000 fewer steps per day) [18] and fewer minutes accumulated in MVPA, whereas more severe recovery clusters were associated with a low frequency of continuous sessions of physical activity per week, longer total sleep time and high variability in sleep timing. Further, the physical activity and sleep profile were notably worse in the PHOSP-COVID cohort than a similarly aged office worker comparator cohort [15], with the 3.2 mg difference in activity volume approximating 1600 fewer steps per day [18]. Relative to our well-matched comparator group with type 2 diabetes [17], differences in activity volume approximated 550 fewer steps per day [18], approximately 17 fewer minutes spent in light-intensity activity, with three times higher variability in sleep timing and lower sleep efficiency.

Among COVID-19 sufferers, sleep disturbance is one of the most commonly reported symptoms, irrespective of acute illness severity, and is highly prevalent following hospital discharge [19]. Irregular sleep patterns as observed in this cohort, independent of total sleep time, can lead to circadian disruption which is a risk factor for metabolic syndrome, obesity, dyslipidemia, and diabetes [20].

Being inactive, defined as not meeting the physical activity guidelines of 150 min of MVPA per week [21], is a risk factor for acute COVID-19 severity, with those who are inactive being 2.26 times more likely to be admitted to hospital, 1.73 times more likely to need intensive care, and 2.49 times more likely to die [4]. Given this, and that habitual physical activity is generally fairly stable [22], it is possible that the lower physical activity in those with more severe acute COVID-19 may reflect their activity levels prior to infection with COVID-19.

The recovery clusters in PHOSP-COVID are not closely associated with acute illness severity [3], consistent with previous research [23]. Thus, the differences observed in continuous bouts of physical activity and sleep across recovery clusters likely reflect the participants’ current mental and physical health impairment. This suggests that the ability to sustain a 10- or 30-min session of activity without resting is compromised in those with more severe ongoing symptoms. Further, the sleep routine appears to be disrupted across all participants within the dataset, but particularly in those in the more severe cluster. Both behaviours are associated with the multiple impairments that characterise the more severe clusters [3], including anxiety and depression [24-26], fatigue [27], and physical function [28]. These data suggest that rehabilitation pathways that have been set up to manage recovery from COVID-19 should focus on the spectrum of behaviours that encompass the 24-h period, including facilitating a return to normal patterns of physical activity, including being able to undertake longer bouts of physical activity, along with focusing on addressing sleep disruption issues.

Strengths and limitations

Key strengths of this study are its size, the comprehensively phenotyped multicentre cohort with novel clinical phenotypes, and accelerometer-assessed physical behaviours at scale. Irrespective, the study has several limitations. Notably, it was not possible to obtain measures of physical behaviours for the participants before they were infected with COVID-19. To account for this, we compared the data to a similarly aged cohort of office workers [15] and a cohort of adults with type 2 diabetes who were well-matched on key characteristics including sex, multimorbidity status, and BMI [17]. However, the data were collected on the comparator cohorts prior to the pandemic; we acknowledge that patterns of physical behaviours may also have been impacted due to the COVID-19 restrictions that have been imposed in the UK (and worldwide) to limit the spread of the virus. Variability in sleep timing and sleep efficiency were the main differences between the PHOSP-COVID and comparator cohorts in the present study. We have previously shown that these sleep-related variables did not differ before and during COVID-19 restrictions, suggesting the differences observed are unlikely due to differences in the measurement period [29]. Due to missing patient-reported outcome data within the PHOSP-COVID cohort, a cluster assignment was not derived for all participants. Finally, over 66% of the cohort had valid accelerometer data but those with data tended to be older, from less deprived communities, and with a lower proportion from ethnic minority communities. Therefore, the data presented may not be generalizable to all those recovering from a hospital admission for COVID-19. Although we used data from a well-phenotyped cohort of patients recovering from a hospital admission for COVID-19, it is not possible to disentangle to what extent these results are specific to COVID-19 or reflect recovery from acute illness requiring hospitalisation more generally.

Conclusions

Survivors of a hospital admission for COVID-19 have low levels of physical activity and significantly disrupted patterns of sleep several months after discharge. Acute illness severity was associated with lower total and moderate-to-vigorous activity following discharge, whereas more severe recovery clusters were associated with substantially fewer bouts of continuous physical activity and greater variability in sleep timing. Without modification, these behaviours are likely to result in further future disease. Physical activity, particularly sustained continuous bouts, and variability in sleep timing are potential treatable traits for survivors of COVID-19.

Additional file 1. PHOSP-COVID Collaborative Group.

Additional file 2. Supplementary material.