Using compositional data analysis to explore accumulation of sedentary behavior, physical activity and youth health.

PURPOSE: The study aimed to describe youth time-use compositions, focusing on time spent in shorter and longer bouts of sedentary behavior and physical activity (PA), and to examine associations of these time-use compositions with cardiometabolic biomarkers.
METHODS: Accelerometer and cardiometabolic biomarker data from 2 Australian studies involving youths 7-13 years old were pooled (complete cases with accelerometry and adiposity marker data, n = 782). A 9-component time-use composition was formed using compositional data analysis: time in shorter and longer bouts of sedentary behavior; time in shorter and longer bouts of light-, moderate-, or vigorous-intensity PA; and "other time" (i.e., non-wear/sleep). Shorter and longer bouts of sedentary time were defined as <5 min and ≥5 min, respectively. Shorter bouts of light-, moderate-, and vigorous-intensity PA were defined as <1 min; longer bouts were defined as ≥1 min. Regression models examined associations between overall time-use composition and cardiometabolic biomarkers. Then, associations were derived between ratios of longer activity patterns relative to shorter activity patterns, and of each intensity level relative to the other intensity levels and "other time", and cardiometabolic biomarkers.
RESULTS: Confounder-adjusted models showed that the overall time-use composition was associated with adiposity, blood pressure, lipids, and the summary score. Specifically, more time in longer bouts of light-intensity PA relative to shorter bouts of light-intensity PA was significantly associated with greater body mass index z-score (zBMI) (β = 1.79; SE = 0.68) and waist circumference (β = 18.35, SE = 4.78). When each activity intensity was considered relative to all higher intensities and "other time", more time in light- and vigorous-intensity PA, and less time in sedentary behavior and moderate-intensity PA, were associated with lower waist circumference.
CONCLUSION: Accumulating PA, particularly light-intensity PA, in frequent short bursts may be more beneficial for limiting adiposity compared to accumulating the same amount of PA at these intensities in longer bouts.

Introduction

An emerging body of work links both physical activity (PA) and sedentary behavior (SED) to health across the lifespan,1, 2, 3 with evidence from studies on adults suggesting that both patterns of accumulation and total volume performed are important.4, 5, 6 Evidence concerning when these movement behaviors should start in youth is important for developing recommendations over the life course. Consequently, there is increasing interest in understanding the way in which daily activity is accumulated, including the timing, duration, and frequency of bouts of PA and SED, and how such patterns may be associated with cardiometabolic health in youth. There is a need to establish at an early age how activity patterns may be associated with cardiometabolic biomarkers, such as obesity and blood lipids, given that they can track into later life.

Evidence from large cohort studies and a meta-analysis suggest that, for adults, high-intensity PA in bouts lasting ≥10 min is associated with a lower mortality risk. Both total SED and longer average SED bout durations are associated with higher mortality risk, but breaking up SED with light-intensity PA (LPA) helps improve adiposity and postprandial glycaemia. Whilst evidence in adults suggests that different PA and SED accumulation patterns play an important role in cardiometabolic health, a systematic review of studies in youth yielded findings that were inconsistent and mostly limited to measures of adiposity. Much remains unknown about whether there may be consequences for cardiometabolic risk, especially in younger populations, depending on whether SED and PA are accumulated sporadically (in shorter bouts) or whether they are accumulated in a more sustained pattern (in longer bouts). Such information is key for the development of targeted interventions that will benefit health outcomes and that will inform the evolution of 24-h movement guidelines.

From a 24-h behavior perspective, sleep and waking movement behaviors occur on an activity continuum that ranges from low- (e.g., sleep, sitting) to high-intensity (e.g., jumping, skipping) activity. A change in time spent in 1 activity intensity (e.g., moderate-to-vigorous-intensity PA (MVPA)) will consequently result in a change in at least one other intensity (e.g., LPA; sleep) in the 24-h period., Thus, considering overall daily time-use composition of total time spent in various activity intensities (i.e., volumes) can yield a different perspective than considering each intensity in isolation. This same principle can be applied to patterns of accumulation, such as time spent in bouts of varying duration for each intensity.

Using compositional data analysis, the relative distribution of time-use between different activities can be examined, classifying activities by their total volume of intensity (e.g., total SED), activity accumulation patterns (e.g., time in longer vs. shorter bouts), other attributes (e.g., behavioral context), or multiple attributes simultaneously., These respective elements of time-use can then be modelled simultaneously and tested for their combined associations with biomarkers of cardiometabolic health.10, 11, 12 This compositional data analysis approach is more advanced than other types of analyses in that it allows for the simultaneous consideration of components that sum to a whole, without statistical problems such as collinearity, which is often encountered in some traditional approaches., To date, the application of compositional data analysis to the study of time-use and health has mostly been applied to adults, rather than to youth. Furthermore, to the best of our knowledge, there have been no studies, in either adults or youth, that have simultaneously considered activity intensity and bout duration components across the activity spectrum. Although 1 previous study has investigated bout duration alongside total volume using a compositional approach and found that engaging in shorter vs. longer SED bouts yielded some benefits in limiting adiposity, it did not investigate other intensity patterns and health outcomes beyond adiposity. Consequently, our study aimed to (1) describe daily time-use composition of youth 7‒13 years old, including time spent in longer and shorter SED and PA bouts, and (2) examine associations between time-use compositions and cardiometabolic risk biomarkers.

Materials and methods

Participant information

Data from 2 cluster-randomized, school-based trials—Transform-Us! (Australian clinical trials registration number12609000715279 (19/08/2009), ISRCTN83725066 (30/06/2010)) and LOOK (ACTRN12615000066583 (23/01/2015))—were pooled for the purpose of this study. Baseline data from 2010 from Transform-Us! and data from 2009 (time-point 5, i.e., 4 years post-randomization) from LOOK were used., Ethical approval was provided by the Deakin University Human Research Ethics Committee (Transform-Us!: EC 2009-141) and by the Australian Capital Territory Health Human Research Ethics Committee (LOOK: ETH.9/05.687). Parents provided written informed consent for their children to participate in each assessment (e.g., accelerometry, anthropometry, lipids). Data were obtained from 1219 participants, but the initial analytic sample included only those participants who had valid accelerometry data (n = 843; 69% of the original sample). Details of each study, and comparisons between studies, are reported elsewhere.

Data collection and measures

Cardiometabolic biomarkers

Seven cardiometabolic measures and related blood biomarkers were collected. Height (cm), weight (kg), and waist circumference (WC; cm) were collected using standardized procedures. Body mass index (BMI; kg/m2) was calculated and converted to z-value, standardized for age and sex, using the World Health Organization Child Growth Standards (zBMI). Systolic blood pressure (SBP) and diastolic blood pressure measurements were taken while participants were in a rested, seated, or supine position. Overnight fasted blood samples were collected for the assessment of high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (lipids)., A continuous combined cardiometabolic risk score (CMR score) was derived using the z-value = (value − mean)/SD for WC, SBP, diastolic blood pressure, LDL-C, HDL-C, and triglycerides, as previously described., Since the spread of observed values and prevalence of deviating values for individual risk factors in a sample of reasonably healthy children are likely to be low compared to adults, the use of this CMR score is appropriate. The calculation of the CMR score is in line with previous calculations used in research on this age group.,

Accelerometry

Participants wore an ActiGraph (ActiGraph, Pensacola, FL, USA) accelerometer on the right hip during waking hours for 8 consecutive days. The output of the different ActiGraph models (GT1M in LOOK and GT3X in Transform-Us!) have been shown to be comparable and can be pooled. ActiGraph data collected in 5-s (LOOK) and 15-s (Transform-Us!) epochs were downloaded and (if applicable) converted to 15-s epochs using ActiLife software (ActiGraph). These were then processed using a customized Excel (Microsoft, Redmond, WA, USA) macro. Sustained periods of ≥20 min of consecutive zeros were used to define non-wear time., Using validated, age-specific cut-points,, total time was subdivided into SED, LPA, moderate-intensity PA (MPA), vigorous-intensity PA (VPA), and “other time” (i.e., non-wear/sleep). These cut-points were developed in a laboratory study that showed that they could be used to distinguish differing levels of PA intensity as well as sedentary time in children. An overview of the age-specific, 15-s epoch cut-points used is provided in Supplementary Table 1. The “other time” category could not be separated out into sleep and non-wear time because sleep data were not collected in either study. Total volumes of SED, LPA, MPA, and VPA were separated into duration of longer and shorter bouts. The thresholds for defining longer and shorter bouts were based on previous literature and preliminary exploration of accumulation patterns in this sample., For SED, longer bouts were defined as ≥5 min and shorter bouts as <5 min. For LPA, MPA, and VPA, longer bouts were defined as ≥1 min and shorter bouts as <1 min. Data related to accelerometer wear time were totaled per day and then averaged across all valid days (i.e., weekdays with ≥8 h and weekend days with ≥7 h of accelerometer wear time). Only participants with ≥4 valid days of accelerometer wear time were included in analyses, with no weekend-day requirement.

Covariates

Participant age was included as a continuous variable, and sex was included as a binary variable. Socioeconomic status (SES) of school location was transformed into deciles based on the national Index of Relative Socio-Economic Disadvantage of the Socio-Economic Indexes for Areas and was categorized as low (1–2), mid (3–8), or high (9 –10) SES. The data source (LOOK or Transform-Us!) was coded as a categorical covariate. We did not adjust the models for the LOOK intervention arm because no differences were observed between arms.

Statistical analyses

Analyses were performed in STATA Version 15.0 (StataCorp, College Station, TX, USA) and RStudio Version 1.4.453 (R Version 3.6.3; The R Foundation for Statistical Computing, Vienna, Austria), using the “compositions” (“acomp” framework), zCompositions (“MultLN”) and “stats (“lm” function) packages.30, 31, 32 For the purpose of our analysis, the proportion of the 24-h day spent in different intensities was normalized for each participant so that their sum equaled 1. A time-use composition of 9 distinct components (time in longer and shorter bouts of SED, LPA, MPA, and VPA, and “other time”) was formed using compositional data analysis. Most activities were ubiquitous, except for the VPA longer bouts component, which was monitored as 0 for 10 participants (1%). A limitation of compositional data analysis is that 0 values cannot be included in log ratios because dividing by 0 or taking the logarithm of 0 are undefined mathematical operations. Accordingly, these 0 were treated as VPA that occurs at a level less than the detection limit (1 epoch over the monitoring period) and thus were replaced using the multiplicative lognormal univariate replacement method of the zCompositions package (detection limit set at the lowest observed value in the sample (0. 143)).

The proportion of time spent in different bouts (longer vs. shorter) of SED, LPA, MPA, VPA, and “other time” were reported using standard descriptive statistics (i.e., median, interquartile range, range) as well as the compositional mean. The compositional mean, or center, is the vector of geometric means of its parts, rescaled to sum up to 24 h, and is coherent with the inter-dependent nature of compositional data.,,

Linear regression models were used to test associations between the time-use composition and the cardiometabolic health biomarkers. All models used robust standard errors to account for clustering by school. Adjusted models controlled for age, sex, SES, and data source (i.e., LOOK or Transform-Us!). The nine-part time-use composition was modelled as a set of 8 isometric log-ratio (ilr) coordinates.34, 35, 36 Values estimated by such models are invariant with regards to the orthonormal basis selected. Accordingly, a basis (i.e., sequential binary partition) consisting of 8 ilr coordinates reflecting comparisons of time spent in longer vs. shorter bouts (ilr coordinates 2 (SED), 4 (LPA), 6 (MPA), and 8 (VPA)), as well as comparisons of volumes of each intensity (SED, LPA, MPA, and VPA) vs. more intense activities and/or “other time” (ilr coordinates 1, 3, 5, and 7), was created. A detailed overview of this basis is displayed in the sequential binary partition and sign matrix in Supplementary Table 2.

The 8 ilr coordinates were modelled simultaneously. Thus, estimates of effects of each bout duration are controlled for other bout durations, the relative volumes of each intensity, and “other time”. Vice versa, effects of relative volumes of each intensity were controlled for “other time”, the other activity intensities, and their accumulation pattern in longer vs. shorter bouts. This sequential binary partition was chosen to allow direct comparisons between time in longer and shorter bouts within each intensity. This is novel within the PA field, where typically ratios between one-versus-remaining behaviors are investigated when using compositional data analysis (e.g.,,).

Global F tests of these 8 ilr coordinates, which reflect the overall effect of the time-use composition, were reported. In addition, individual regression coefficients for each of the 8 ilr coordinates were reported. Because the change in outcome depends on which other compositional parts change to compensate, plots with estimated outcomes were created. These provide a more interpretable indication of relationships between activity accumulation pattern (more time in longer vs. shorter bouts) and cardiometabolic biomarkers compared with the estimated regression parameters for individual ilr coordinates. These plots included the estimated mean cardiometabolic health biomarker, with 95% confidence intervals (95%CIs). In these estimations, the overall time within the intensity of interest was kept constant, while the ratio between longer and shorter behaviors was alternated. Times in all other components were set at the mean observed values. Means were estimated for girls and boys separately, keeping age (mean = 10.5 years), SES (mid) and data source (LOOK) constant (Table 1). Estimations were only conducted for the range of activity levels occurring in the observed data, rounded to the nearest 5 min for plotting purposes. That is, data were not extrapolated beyond the observed range for the exposure variables (Table 2).

Table 1

Participant characteristics (subset adiposity, n = 782; mean ± SD or %).

Demographic characteristics	Whole sample	Girl (n = 431)	Boy (n = 351)
Age (year)	10.5 ± 1.7	10.3 ± 1.7	10.6 ± 1.7
Gender (% female)	55	100	0
SES (%)a
High	3	3	3
Mid	36	36	35
Low	61	61	62
Data source (% in LOOK)	53	50	56
Cardiometabolic biomarkers
zBMI	0.5 ± 1.1	0.5 ± 1.1	0.6 ± 1.2
WC (cm)	64.1 ± 9.0	63.3 ± 9.0	65.1 ± 8.9
Systolic blood pressure (mm Hg)b	106.5 ± 10.3	105.6 ± 10.4	107.5 ± 10.0
Diastolic blood pressure (mm Hg)b	61.0 ± 7.5	61.0 ± 7.4	61.0 ± 7.6
HDL-C (mmol/L)c	1.5 ± 0.3	1.4 ± 0.3	1.5± 0.3
LDL-C (mmol/L)c	2.5 ± 0.7	2.6 ± 0.7	2.5 ± 0.7
Triglycerides (mmol/L)d	0.8 ± 0.4	0.9 ± 0.4	0.8 ± 0.4
CMR scored	0.0 ± 3.4	0.2 ± 3.2	-0.1 ± 3.7

Abbreviations: CMR score = cardiometabolic risk score; HDL-C = high-density lipoprotein cholesterol; LDL-C = low-density lipoprotein cholesterol; SES = socioeconomic status; WC = waist circumference; zBMI = body mass index converted to the World Health Organization Child Growth Standards age and sex standardized z-values.

SES of school location was transformed into deciles based on the Index of Relative Socio-Economic Disadvantage of the national Socio-Economic Indexes for Areas and was categorized as low (1–2), mid (3–8), or high (9–10) SES.

n = 637 (55% girls).

n = 525 (55% girls).

n = 404 (54% girls).

Table 2

Standard and compositional descriptive statistics of the proportion of time spent in longer and shorter bouts of SED, LPA, MPA, VPA, and “other time” (subset adiposity, n = 782).

	Standard descriptive statistics (min/day)				Compositional descriptive statistics
	Q1	Median	Q3	Range	Mean (%)	Mean (min/24h)
SED
Longer bouts	122.8	160.2	208.1	28.0–512.2	12	169.6
Shorter bouts	243.6	265.2	293.1	162.7–465.7	19	268.0
LPA
Longer bouts	87.9	101.4	119.3	44.0–197.3	7	103.9
Shorter bouts	114.8	126.4	136.9	79.6–189.7	9	126.0
MPA
Longer bouts	5.5	7.7	10.6	1.0–40.2	<1	8.6
Shorter bouts	27.2	36.5	47.1	6.4–82.5	3	37.4
VPA
Longer bouts	1.9	4.1	7.5	0.1–49.1	<1	5.8
Shorter bouts	8.5	13.5	19.6	1.5–43.6	1	14.9
“Other time”	614.1	668.2	718.9	301.8–889.2	46	664.3

Notes: Longer bouts of SED were defined as ≥5-min bouts. Shorter SED was defined as bouts of <5 min. Similarly, longer PA bouts (i.e., LPA, MPA, and VPA) were defined as ≥1-min bouts and shorter PA bouts as the time accumulated in <1-min bouts. “Other time” consisted of non-wear/sleep.

Abbreviations: LPA = light-intensity physical activity; MPA = moderate-intensity physical activity; Q1 = the first interquartile; Q3 = the third interquartile; SED = sedentary behavior; VPA = vigorous-intensity physical activity.

Results

Participant characteristics

Of the 843 participants with valid accelerometry data, 61 (7%) had missing data for BMI, WC, and/or covariates, and were therefore excluded from data analysis. A total of 782 participants were included in the analysis (subset adiposity). Next, subsets with complete data were created for blood pressure (n = 637 (76% of the sample with valid accelerometry)), lipids (n = 525 (62%)), and CMR score analysis (n = 404 (48%)). These smaller analytic samples were attributable to fewer children participating in blood pressure and blood assessments. The included sample for each subset was compared with the excluded participants, and no major differences were observed (comparisons were previously published). Demographic characteristics and cardiometabolic biomarker values for girls vs. boys were comparable, and thus the data were pooled for analysis (Table 1).

Participants were on average 10 years old (55% girls), and most (61%) attended schools located in low SES areas relative to the Australian average (Table 1). Standard and compositional descriptive statistics of the proportion of time spent in longer and shorter bouts of SED, LPA, MPA, and VPA are presented in Table 2. Approximately 12% and 19% of the total time per 24-h day was spent in longer and shorter SED bouts, respectively. The proportions of time spent in longer and shorter bouts of LPA, MPA, and VPA ranged from <1% to 9%.

Associations between time-use composition and cardiometabolic biomarkers

Results from the linear regression analysis modelling the relationship between the time-use compositions, including the 8 ilr coordinates, and cardiometabolic biomarkers are provided in Table 3 and Supplementary Table 3. Table 3 presents the global F test and the parameters corresponding to time in longer vs. shorter SED, LPA, MPA, and VPA (ilr coordinates 2, 4, 6, and 8). Supplementary Table 3 presents the parameters corresponding to relative volumes of SED, LPA, MPA, and VPA (ilr coordinates 1, 3, 5, and 7). These were obtained within the same time-use composition models as the parameters presented in Table 3. While all results are presented in Table 3 and Supplementary Table 3, the following text only focuses on the overall time-use composition models and the individual ilr regression coefficients within the confounder-adjusted models, which showed evidence of an association between the time-use composition and the cardiometabolic biomarkers (global F test p < 0.05; specific values reported in Table 3).

Table 3

Associations of the ratio of time in longer vs. shorter bouts, controlled for all other intensities, with cardiometabolic biomarkers in children aged 7–13 years.

		SED (ilr 2)		LPA (ilr 4)		MPA (ilr 6)		VPA (ilr 8)		Overall
Outcome	Model	β (SE)	p	β (SE)	p	β (SE)	p	β (SE)	p	F(df)	p	Adj. R2
zBMI	Unadj.	–0.28 (0.48)	0.558	1.43 (0.65)	0.027	–0.53 (0.27)	0.050	0.23 (0.16)	0.146	6.74(8,773)	<0.001	0.06
(n = 782)	Adj.	–0.33 (0.48)	0.489	1.79 (0.68)	0.009	–0.01 (0.30)	0.983	0.14 (0.16)	0.395	5.78(13,768)	<0.001	0.07
WC (cm)	Unadj.	–8.83 (3.40)	0.010	27.26 (4.62)	<0.001	1.36 (1.94)	0.481	0.70 (1.12)	0.530	30.41(8,773)	<0.001	0.23
(n = 782)	Adj.	–5.99 (3.37)	0.076	18.35 (4.78)	<0.001	–2.61 (2.12)	0.220	1.91 (1.12)	0.087	22.61(13,768)	<0.001	0.26
SBP (mm Hg)	Unadj.	–4.04 (4.28)	0.347	18.91 (5.84)	0.001	12.26 (2.54)	<0.001	–3.75 (1.52)	0.014	22.28(8,628)	<0.001	0.21
(n = 637)	Adj.	–1.40 (4.19)	0.738	8.61 (5.88)	0.144	4.42 (2.76)	0.110	–1.52 (1.50)	0.313	18.93(13,623)	<0.001	0.27
DBP (mm Hg)	Unadj.	0.16 (3.52)	0.964	–1.08 (4.81)	0.823	1.20 (2.09)	0.566	0.59 (1.25)	0.635	0.92(8,628)	0.497	0.00
(n = 637)	Adj.	0.14 (3.58)	0.969	–0.71 (5.03)	0.888	1.61 (2.36)	0.495	0.77 (1.28)	0.549	0.90(13,623)	0.550	0.00
HDL-C (mmol/L)	Unadj.	–0.07 (0.16)	0.661	–0.30 (0.21)	0.151	–0.28 (0.09)	0.001	0.00 (0.05)	0.999	14.17(8,516)	<0.001	0.17
(n = 525)	Adj.	–0.08 (0.16)	0.593	–0.14 (0.22)	0.535	–0.09 (0.09)	0.349	–0.04 (0.05)	0.470	10.70(13,511)	<0.001	0.19
LDL-C (mmol/L)	Unadj.	–0.03 (0.35)	0.948	0.27 (0.48)	0.576	–0.52 (0.20)	0.008	0.21 (0.11)	0.063	2.08(8,516)	0.036	0.02
(n = 525)	Adj.	–0.08 (0.36)	0.820	0.40 (0.52)	0.437	–0.51 (0.22)	0.021	0.21 (0.12)	0.074	1.42(13,511)	0.146	0.01
Triglycerides (mmol/L)	Unadj.	0.26 (0.19)	0.180	0.07 (0.26)	0.791	0.18 (0.11)	0.093	0.09 (0.06)	0.155	8.13(8,516)	<0.001	0.10
(n = 525)	Adj.	0.17 (0.20)	0.379	0.16 (0.28)	0.564	0.08 (0.12)	0.481	0.08 (0.06)	0.203	6.80(13,511)	<0.001	0.13
CMR	Unadj.	–0.51 (1.76)	0.770	5.67 (2.40)	0.018	1.82 (1.02)	0.075	0.64 (0.63)	0.313	18.70(8,395)	<0.001	0.26
(n = 404)	Adj.	0.16 (1.77)	0.704	2.78 (2.53)	0.272	–0.33 (1.12)	0.766	1.19 (0.64)	0.062	13.70(13,390)	<0.001	0.29

Notes: Bold values denote statistical significance at the p < 0.05 level. Unadjusted and adjusted models accounted for school clustering using robust standard error and included effects of the time-use composition tested as a series of ilr coordinates. Adjusted models additionally adjusted for data source (LOOK or Transform-Us!), age, sex and SES. All coordinates represent ilr of simplex coefficients of time in longer vs. shorter bouts of 1 intensity, keeping the total time within that intensity constant. The β-estimates evaluate associations between an outcome and an increase in time accumulated in longer vs. shorter bouts. Longer bouts of SED were defined as bouts ≥5 min. Shorter SED was defined as bouts of <5 min. Similarly, longer PA bouts (i.e., LPA, MPA, and VPA) were defined as bouts ≥1 min, and shorter PA bouts were defined as the time accumulated in bouts of <1 min. “Other time” consisted of non-wear/sleep.

Abbreviations: Adj. = adjusted; CMR score = cardiometabolic risk score; DBP = diastolic blood pressure; df = degree of freedom; HDL-C = high-density lipoprotein cholesterol; ilr = isometric log ratio; LDL-C = low-density lipoprotein cholesterol; LPA = light-intensity physical activity; MPA = moderate-intensity physical activity; SBP = systolic blood pressure; SED = sedentary behavior; Unadj. = Unadjusted; VPA = vigorous-intensity physical activity; WC = waist circumference; zBMI = body mass index converted to the World Health Organization Child Growth Standards age and sex standardized z-values.

The unadjusted and adjusted models showed evidence of an association of the overall time-use composition with 6 of 8 cardiometabolic health biomarkers: zBMI, WC, SBP, HDL-C, triglycerides, and CMR score (global F test p < 0.05; specific values reported in Table 3). A significant association of the time-use composition with LDL-C was seen in the unadjusted model (global F test p = 0.036), but was attenuated with adjustment for confounders (global F-test p = 0.146). Closer inspection of the individual ilr coordinates when the test of the full composition was significant (zBMI, WC, SBP, HDL-C, triglycerides, and CMR score) (Table 3) showed that a higher ratio of time in longer relative to shorter bouts of LPA (Long LPA: Short LPA) was associated with higher zBMI (β = 1.79; SE = 0.68; p = 0.009) and WC (β = 18.35; SE = 4.78; p < 0.001). Thus, higher levels of continuous LPA relative to sporadic LPA was associated with higher estimated zBMI and WC. The associations for other individual ilr coordinates, across the other significant, adjusted time-use models (global F test p < 0.05), were smaller and not statistically significant (all p ≥ 0.05; specific values reported in Table 3). The unadjusted models revealed some further associations that were attenuated in size and no longer statistically significant after adjustment for confounders.

More time in LPA relative to MPA, VPA, and “other time” was associated with lower WC (β = ‒9.10; SE = 3.75; p = 0.016). More time in VPA relative to “other time” was also associated with lower WC (β = ‒2.87; SE = 1.16; p = 0.014). In addition, more time in SED vs. all other behaviors, and more time in MPA relative to VPA and “other time”, was associated with higher WC (β = 7.60; SE = 3.34; p = 0.023 and β = 4.08; SE = 1.84; p = 0.027; respectively). After accounting for confounders, no evidence of an association between any of the other individual ilr coordinates with zBMI, WC, SBP, HDL-C, triglycerides, or CMR score, was observed (Supplementary Table 3).

Estimated cardiometabolic biomarkers based on different ratios of time spent in longer vs. shorter bouts

Estimated mean zBMI and WC for different ratios of longer and shorter LPA are visually presented in Fig. 1. Only these individual ilr coordinates indicated significance (p < 0.05) within the adjusted models, with an observed overall significance between the time-use composition and a cardiometabolic health marker (global F test p < 0.05). All other estimations are visualized in the Supplementary Figs. 1‒8.

Fig. 1

(A) Estimated zBMI and (B) WC for different ratios of longer and shorter LPA bouts. Estimated adiposity markers for different ratios of longer and shorter LPA bouts. Plots include estimated mean zBMI and WC with 95% confidence intervals. The overall time within the intensity of interest was kept constant, while the ratio between longer and shorter behaviors was alternated. Times in all other components were set at the mean observed values. Blue stands for boys, red stands for girls, and grey stands for overlapping 95% confidence intervals. LPA = light-intensity physical activity; WC = waist circumference; zBMI = body mass index converted to the World Health Organization Child Growth Standards age and sex standardized z-values.

Fig. 1A, 1B shows that, relative to the existing compositional mean, performing an additional 10 min of LPA in 1 or more continuous longer bouts rather than in shorter bouts was consistent with approximately 0.2 higher mean zBMI (0.7 vs. 0.9, respectively) and 2.2 cm higher mean WC (65.3 cm vs. 67.5 cm, respectively) in girls. Similar estimations were predicted for boys (Fig. 1).

Discussion

This study examined daily time-use compositions of youth 7–13 years old and their associations with cardiometabolic biomarkers. The study is novel in that it not only factored in activity intensity (SED, LPA, MPA, and VPA) but also determined whether these activity intensities were accumulated in longer bouts (≥5-min SED and ≥1-min PA) or shorter bouts (<5-min SED and <1-min PA) within a 24-h composition. In particular, the novel sequential binary partition chosen was advantageous for the comparison of relative time spent in longer vs. shorter bouts. Although this provided a slightly complicated view of the relative volume of time spent in each activity intensity versus others, it makes the current research an original contribution to activity pattern research in youth.

The results demonstrated that the overall time-use composition was associated with most cardiometabolic biomarkers, with one of the important components of the composition being LPA bout duration, particularly in relation to adiposity. Results indicated that the same amount of LPA accumulated with a more continuous accumulation pattern (i.e., more time in longer rather than shorter bouts) was linked with poorer WC and zBMI. While the cross-sectional design limits causal inference, it is possible that performing a specific amount of LPA frequently in shorter bursts rather than in longer bouts at one time may be beneficial for limiting adiposity.

The results also suggested that different relative time-use of total volumes (i.e., long and short bouts combined) of SED, LPA, MPA, and VPA were important components of the time-use composition. Specifically, more SED (less PA and “other time”) and more MPA (less VPA and “other time”) were associated with higher average WC, and more VPA (less “other time”) and, somewhat unexpectedly, more LPA (less MPA, VPA, and “other time”) were associated with lower WC. Effect sizes were in a similar direction for other health markers, but were not statistically significant. Taken collectively, the findings across the different ilr parameters suggest that the most favorable cardiometabolic biomarker levels were observed with higher high-intensity PA (VPA) and lower sedentary time. These findings are consistent with previous research using similar methods (e.g., 15-s epoch accelerometer data and compositional data analysis,) and confirm that participating in high-intensity PA, including MPA and VPA, may be important for youth health.

Comparisons to previous literature

To the best of our knowledge, no research has used compositional data analysis to examine bout durations across the activity spectrum, including SED and PA bouts, to assess associations with a range of cardiometabolic outcomes. One study by Gába and colleagues used a similar approach, except they only investigated SED bout durations alongside total volumes of other intensities and their associations with adiposity. Their results suggested that adiposity could be improved by replacing middle-length SED bouts (defined as 10–29 min) with shorter SED bouts (1–9 min), yet they found no associations for replacements with long SED bouts (≥30 min duration). This contrasts with our study, which did not find significant associations for a particular (shorter or longer) type of SED bouts. These contrasting results may be explained by methodological differences, such as the sequential binary partition chosen and the different thresholds for longer vs. shorter bouts. Because we identified no other studies that have used compositional analysis to assess bout durations in youth, it is difficult to further compare our findings, or the findings of Gába and colleagues, with additional research.

Nevertheless, a few studies have investigated youth activity patterns and associations with cardiometabolic biomarkers, using different analytical methods. For example, Holman and colleagues used logistic regression models to evaluate all MVPA (i.e., total volume), and MVPA accumulated in short bouts (i.e., <5 min, <10 min) or long bouts (i.e., ≥5 min, ≥10 min) among participants who were 6‒19 years old and did not observe differences in high cardiometabolic risk vs. normal cardiometabolic risk. Despite some methodological differences (e.g., 60-s epoch vs. 15-s epoch; different thresholds for bout durations), their finding are in line with the findings in our study, which also failed to detect evidence to suggest that certain MPA and/or VPA bouts are more strongly associated with health than others.

Our findings are also partially in line with a study by Aadland et al., who compared shorter and longer bouts of PA in youth using a multivariate pattern analysis approach. Whereas our study only found significant effects of bout length for LPA, the Aadland et al. study found that short intermittent bursts of PA of all intensities—including MPA and VPA—were favorable for youth cardiometabolic health. Aadland et al. explained their findings through the existence of an interrelation between LPA and VPA, with short bouts of LPA being positively related to VPA and vice versa. Whilst this finding was not specifically observed in our sample, the similar LPA observations across these 2 studies suggest that the accumulation of LPA may be relevant for youth health and therefore warrants further investigation.

Practical and public health implications

Our findings suggest that activity patterns may play a role in children's cardiometabolic health—and particularly adiposity—beyond total volumes of PA and SED. Our results were more supportive of encouraging children to accumulate PA through facilitating their natural sporadic accumulation patterns, rather than by trying to alter these patterns by accruing PA in longer bouts. One important caveat in the practical application of this type of evidence is that statistical models of effects of patterns at a fixed volume are hypothetical scenarios of what might be expected. In practical, real-life terms, changes in volume may well occur alongside any intervention that targets accumulation patterns because these 2 characteristics of PA are intrinsically linked. Although the main focus of our study was on PA and SED bouts, and not necessarily on total volumes, the results do suggest that it is not just bouts that are important for health. Specifically, more relative time in total VPA (less “other time”) and LPA (less MPA, VPA, and “other time”) was associated with lower WC, and more relative time in SED (less PA and “other time”) and MPA (less VPA and “other time”) was associated with higher WC. While further evidence is required to determine the impact of longer activity bouts on children's health, the present findings support the recent removal of a minimum threshold of 10 min sustained MVPA in the US guidelines; instead, the focus should be on total volume of MVPA.

Strengths and limitations

The main strengths of our study are the large sample size and the objective measurement of activity patterns and cardiometabolic biomarkers. The main limitations in most previous PA and SED accelerometer and health studies are the combined assessment of intensities, adjustment for other intensities and the consequent potential collinearity issues arising from these factors (particularly when investigating SED with the inclusion of PA adjustments). This issue was overcome in our study through the use of compositional data analysis, which allows for handling of multiple PA and SED patterns within a joint statistical model.

Our study also had some limitations. First, cross-sectional data were pooled. Thus, the estimated differences in cardiometabolic biomarkers cannot be directly interpreted as an effect of time reallocation from one component to another. Second, our study was not designed to explain the possible mechanisms by which significant effects of bout duration were seen for adiposity but not for blood biomarkers (e.g., lipids, CMR score). However, we might speculate that the participants’ age range and their limited cumulative exposure to unhealthy lifestyle behaviors, such as extensive sitting, may have been contributing factors. In addition, there was a higher power for detection of potential associations with adiposity vs. other cardiometabolic biomarkers, which was due to the adiposity subset being a larger subset compared to those who consented to blood collection. Longitudinal studies with larger samples with data on risk factors other than adiposity are needed in order to investigate the long-term health effects of continuous and sporadic activity patterns. Third, the measurement protocol used in our study assumes that habitual patterns were captured, yet it only truly reflects the capture of patterns for the past 4–7 days. Fourth, as in much of the current literature, the data used in our study came from behaviors classified by waist-worn accelerometers processed by applying thresholds to epoch data rather than being measured by posture-based devices. This measurement approach has acceptable validity for capturing total volume of PA but has limited validity for measuring accumulation patterns such as SED bout durations. The accelerometer data used in our study was not examined as raw data, was not collected 24 h/day and did not include heart rate data. Consequently, this may have failed to capture some aspects of these youths’ highly intermittent accumulation patterns, their physiological responses related to certain accumulation patterns and potentially important differences between sleep and non-wear time. Because physiology and sleep are important factors in youth health,, future studies should consider using continuous observations with 24-h wear protocols. Further improvements could come from using posture-based devices (such as activPAL, http://www.palt.com/) to measure SED; ideally, heart rate could also be monitored, thus capturing physiological responses. Fifth, although models were adjusted for age, it was not possible to adjust for puberty because this was not assessed in the Transform-Us! study. It should be noted that the Transform-Us! participants were younger than the LOOK participants. Biological differences in maturity may have influenced results, but chronological age is not a strong indicator for maturity stage. In addition, no data on ethnicity and/or race were collected in our study. Future studies could incorporate these characteristics. Finally, because our study was exploratory and a large number of tests were performed without significance adjustment, there is an increased likelihood of false discovery due to multiple testing. Nevertheless, because behavioral research in youth is an understudied area, a deliberate decision was made not to adjust for multiple comparisons., Although the results presented here focus on statistically significant findings (using p values, which can be problematic), the full results, including effect sizes, standard errors, and plots, are available in the tables included in the Supplementary Appendix, which allows for interpretation with and without p values.

Conclusion

In summary, this study was the first, to our knowledge, to investigate associations between patterns of PA and SED and a range of cardiometabolic health markers using compositional data analysis in youth. The main findings were that the time-use composition, which specifically included time in longer and shorter bouts of SED and PA, was significantly associated with most cardiometabolic biomarkers. Specifically, the accumulation of the same amount of LPA in a more sporadic rather than sustained manner (i.e., shorter vs. longer bouts) was associated with lower zBMI and WC. While the findings are only cross-sectional, they suggest that accumulating LPA frequently in short bursts may result in less adiposity than engaging in the same amount of LPA in longer bouts less regularly. However, the preliminary findings of this study and the extant literature should be corroborated (or refuted) with evidence from other samples, collected using alternative device-based measurements (especially those having high validity for accumulation patterns and assessment of SED), and that have measurements of sleep and longitudinal designs. Additional research is warranted in order to develop targeted interventions that benefit health outcomes and to establish whether existing movement guidelines should be refined to include recommendations relating to specific activity patterns and intensities.