Joint associations between objectively measured physical activity volume and intensity with body-fatness. The Fenland Study =========================================================================================================================== * Tim Lindsay * Katrien Wijndaele * Kate Westgate * Paddy Dempsey * Tessa Strain * Emanuella De Lucia Rolfe * Nita G Forouhi * Simon Griffin * Nick J Wareham * Søren Brage ## Abstract **Background/Objectives** Physical activity energy expenditure (PAEE) represents the total volume of all physical activity. This can be accumulated as different underlying intensity profiles. Although volume and intensity have been studied in isolation, less is known about their joint association with health. We examined this association with body-fatness in a population-based sample of middle-aged British women and men. **Methods** 6148 women and 5320 men from the Fenland study with objectively-measured physical activity from individually calibrated combined heart rate and movement sensing and DXA-derived body-fat percentage (BF%) were included in the analyses. We used linear and compositional isocaloric substitution analysis to examine associations of PAEE and its intensity composition with body-fatness. Sex-stratified models were adjusted for socio-economic and dietary covariates. **Results** PAEE was inversely associated with body-fatness in women (beta=-0.16 (95%CI: −0.17; −0.15) BF% per kJ·day-1·kg-1) and men (beta=−0.09 (95%CI: −0.10; −0.08) BF% per kJ·day-1·kg-1). Intensity composition was significantly associated with body-fatness, beyond that of PAEE; the reallocation of energy to vigorous physical activity (>6 METs) from other intensities was associated with less body-fatness, whereas light activity (1.5-3 METs) was positively associated. However, light activity was the main driver of overall PAEE volume, and the relative importance of intensity was marginal compared to that of volume; the difference between PAEE in tertile 1 and 2 in women was associated with 3 percentage-point lower BF%. Higher vigorous physical activity in the same group to the maximum observed value was associated with 1 percentage-point lower BF%. **Conclusions** In this large, population-based cohort study with objective measures, PAEE was inversely associated with body-fatness. Beyond the PAEE association, greater levels of intense activity were also associated with lower body-fatness. This contribution was marginal relative to PAEE. These findings support current guidelines for obesity prevention which emphasise moving more over the specific intensity or duration of that activity. Keywords * Physical activity * body fatness * adiposity * compositional data * isocaloric * isotemporal * adults ## Introduction Physical activity energy expenditure (PAEE) is the most variable component of total energy expenditure (TEE) and is within voluntary control in free-living settings. Physical activity is also characterised by its intensity, which is usually expressed as the metabolic equivalent of task (MET), with 1 MET corresponding to resting metabolic rate. As energy expenditure is a function of intensity and time, it follows that many different behavioural patterns, or intensity profiles, can underpin any given total volume of PAEE. For example, large amounts of time spent in light physical activity (LPA), or a lesser amount of time spent in moderate physical activity (MPA), or an even shorter period of time spent in vigorous physical activity (VPA) could all lead to the same PAEE. Several studies have shown that physical activity is associated with a range of metabolic conditions, including obesity or hard end-points such as mortality1–5. Typically, overall volume and intensity of physical activity (PA) have been examined in parallel, and less is known about the role of intensity beyond that of volume. In order to advance this understanding, these two dimensions of physical activity need to be examined in an integrated model. Recently, a number of isotemporal substitution studies have sought to examine the associations between the reallocation of time in one type of behaviour for another. For example, time spent sedentary for time in moderate-to-vigorous physical activity (MVPA)4, 6–8. However, such models do not account for the inherent increase in PAEE that occurs when a fixed amount of time is reallocated to a more intense behaviour. Simultaneous examination of volume and intensity therefore requires a methodological adaptation to enable integrated analysis. Examples of comparable integrated analyses exist in nutritional epidemiology, where the combined role of total energy intake and macronutrient energy composition are studied simultaneously in the general population9, 10. In essence, these studies evaluate whether it is the total number of calories, or their source or exchanges between the different sources, which drives associations with health. In contrast, little is known about the role of the energy expenditure composition of physical activity in determining body-fatness, beyond that of total PAEE. Investigating this in free-living individuals requires valid energy expenditure estimates across the full spectrum of intensity and well-characterised body-fatness in population studies, large enough that different combinations of high and low PA volumes with different intensity profiles naturally occur. Here, we aimed to evaluate the integrated joint associations between objectively measured PA volume, intensity and body-fatness in a large sample of free-living adults. ## Methods ### Participants The Fenland Study is an ongoing population-based observational study of 12 435 young and middle-aged adults, the research methods for which have been previously described (DOI: [https://doi.org/10.1186/ISRCTN72077169](https://doi.org/10.1186/ISRCTN72077169))11. Briefly, participants born between 1950 and 1975 were recruited between 2005 and 2015 from general practice lists in the East of England, UK. Exclusion criteria were pregnancy, physician-diagnosed diabetes, inability to walk unaided, psychosis, and terminal illness. All participants provided written informed consent and the study was approved by the local ethics committee (NRES Committee – East of England Cambridge Central) and performed in accordance with the Declaration of Helsinki. ### Measures #### Anthropometry and other clinical measures Participants attended a clinical research facility after an overnight fast. All measurements were taken by trained research staff following standardised procedures. Height (cm) was measured with a rigid, wall mounted stadiometer (SECA 240; Seca, Birmingham, UK) and weight (kg) was measured in light clothing with calibrated scales (TANITA model BC-418 MA; Tanita, Tokyo, Japan). Total body fat mass (FM, in grams) was determined by dual-energy X-ray absorptiometry using a Lunar Prodigy Advanced with the enCORE™ software version 14.10.022 (GE Medical Healthcare, Hatfield, UK) as validated against the gold-standard 4-compartment method12. Standard imaging and positioning protocols were applied. Briefly, the system was first calibrated using a spine phantom made of calcium hydroxyapatite, embedded in a lucite block. The coefficient of variation for scanning precision was 2% for total fat mass (30 consecutive scans). To scan participants, they were positioned lying supine within the scanning area of the DXA bed. Anatomical regional boundaries were demarcated and corrected if necessary. For the present analysis, we excluded those with missing DXA data (n=566) but included DXA-scanned participants with medical implants, amputees and minor scanning artefacts (n=135). For participants who were too large to fit within the scanning area (n=108), the symmetry method was used during image processing, e.g. an unscanned left arm was assumed to match the right arm. The primary outcome variable was body fat percentage (BF%) and the secondary outcome was fat mass index (FMI, calculated as fat mass divided by height squared). #### Physical activity assessment Physical activity was measured objectively by fitting participants with a combined heart rate and uniaxial movement sensor (Actiheart, CamNtech, Papworth, UK), attached to the chest with standard ECG electrodes13. Heart rate was individually calibrated using a treadmill test as previously described14. At the end of the clinical visit, participants were asked to wear the sensor, initialised to collect data at 1-min resolution, for the following 6 days, and to return the monitor by freepost. Participants were advised that the device was waterproof and should be worn continuously, including during showering, water-based activities, and sleeping, whilst continuing with their usual activities. It could be removed to change electrodes, spare sets of which were provided. Heart rate data were pre-processed15, individually calibrated14, and combined with acceleration to estimate instantaneous PAEE (intensity) according to methods previously described and validated16, 17. Energy spent at each PA intensity level was calculated as a fraction of PAEE by dividing energy expended at that level by total PAEE. We did this for multiple intensity categories in high resolution (every 0.25 MET) and also grouped into the broader categories of sedentary behaviour/sleep (SS: <1.5 METs), light physical activity (LPA: 1.5-3 METs), moderate physical activity (MPA: 3-6 METs), and vigorous physical activity (VPA: >6 METs), where 1 MET = 71 J·min-1·kg-1 (∼3.5 ml O2·min-1·kg-1). As a sensitivity analysis, we redefined the non-sedentary categories as LPA (1.5-4 METs), MPA (4-7 METs), and VPA (>7 METs). For individuals registering no MPA or VPA time, these analytical fractions were replaced with a value of 0.0001 below the lowest recorded non-zero value for the population, since some analyses require log-transformation of these exposure variables. Participants were excluded from the present analysis if they had failed to wear their sensor for at least 72 hours overall (n=319), or at least 8 hours cumulative wear for each quadrant of the day (3am to 9am, 9am to 3pm, 3pm to 9pm, 9pm to 3am); these latter criteria ensure behavioural information is available for all parts of the day from at least two different days. Furthermore, activity records were excluded if they did not measure 0 m·s-2 (zero movement) at some point during the monitoring period as this indicates a technical problem with the acceleration sensor (n=76). #### Covariates Demographic, lifestyle and health variables were collected using self-report. These included age, sex, marital status (single, married, widowed/separated/divorced), education (compulsory, further – A-level/apprenticeship/sub-degree level, higher – degree level or above), household income level (<£20,000, £20,000-£40,000, >£40,000), ethnicity (White, South Asian, Black, East Asian, Others and unknown), location (Cambridge, Ely, Wisbech), and smoking status (never, ex-smoker, current). Participants with missing socio-demographic covariate data were retained for analysis; missing data in categorical variables were coded as a separate category. Habitual diet over the previous year was self-reported using a validated food frequency questionnaire (FFQ)18, from which estimated total energy intake (kJ·day-1) from carbohydrates, protein, fats, and alcohol were derived. Participants with no dietary information were excluded from the analysis (n=6). Plasma vitamin C was measured as a biomarker of fruit and vegetable intake and indicator of overall diet quality19. Participants with plasma vitamin C levels below the assay detection threshold were coded at the minimum detectable level (n=45). Participants missing this measure were included in analysis by imputation from age, sex, self-reported fruit and vegetable intake and vitamin supplementation (n=263). ### Data analysis All analyses were stratified by sex owing to fundamentally different adiposity patterns between men and women. For sample characteristics, means and standard deviations are reported for continuous variables and proportions are reported for categorical variables. All analyses were undertaken in Stata 14 (StataCorp, Texas). #### Associations between PAEE and intensity of PA on body-fatness Multivariable linear regression was used to examine the associations between body-fatness outcomes (BF% and FMI) and the PA exposure variables, with all models controlling for the potential confounding variables described above. First, we examined the association between total PAEE as a single exposure and both adiposity outcomes. We then explored the relative association of PA intensity through a series of sequential linear regressions for the fraction of PAEE spent above a given intensity threshold (1-11 METs), controlling for overall PAEE. Beta coefficients were graphically represented to estimate the difference in BF% associated with a 1% difference in PAEE above a given MET threshold. #### Isocaloric substitution analysis To assess the combined association between the relative PA intensities and body-fatness, we simulated the effect of reallocating energy expenditure between the four broader intensity categories on body adiposity, using two types of models. First, we used linear isocaloric substitution analysis as per methods previously described20, except with percentage of PAEE substituted, rather than time. In brief, each iteration of the model simultaneously included total PAEE and its fractions expended in the intensity categories, systematically excluding (dropping) the fraction of energy expended in one intensity category. Regression coefficients (95% CI) of the intensity components included in these models provide an estimate of the change in the outcome variable when reallocating energy expended in the dropped intensity to another by percentage of total PAEE, while decreasing energy expended in the dropped intensity category by the same percentage. We checked appropriateness of the linear approximation of these relationships by visual inspection of scatterplots for each exposure-outcome, controlled (residualised) for the other covariates in the model. As an alternative approach to assess the combined association between the relative PA intensities and body-fatness, we modelled the composition of PAEE using a compositional data analysis approach according to methods previously described4, 7. Fractional energy spent at each intensity of PA was transformed into sets of isometric log-ratios (ILRs). Each set of ILRs was regressed against the outcome of interest, controlling for total PAEE and confounders to produce a set of beta coefficients. These beta coefficients were then exponentiated (back transformed) and scaled with the geometric mean of PAEE in the relevant energy intensity for the population stratum of interest (men, women, tertile of PAEE). These results were then plotted against reallocations of PAEE accumulated at the given intensity in kJ·day-1·kg-1, whilst holding total PAEE constant. To provide an equivalent of linear isocaloric substitution using compositional analysis, pairwise reallocations were also modelled6. For all models, variance inflation factors were checked to assess potential issues of collinearity, and distribution of residuals were visually inspected for approximation to normality. ## Results ### Participant characteristics A total of 11 468 participants (6148 women, 5320 men) were included in the analyses (Table 1). Overall, women and men were of similar age, but women had lower body mass index (BMI) and higher FMI and BF% than men. View this table: [Table 1:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T1) Table 1: Participant characteristics. The Fenland study 2005-2015. Women also accumulated lower levels of total PAEE, with a mean (SD) of 50 (20) compared to 59 (23) kJ·day-1·kg-1 in men. Women and men accumulated PAEE differently, with LPA, MPA, and VPA respectively constituting 59%, 27%, and 3% of women’s total PAEE on average; the comparable figures for men were 51%, 35%, and 7%. ### Associations between PAEE and intensity of PA on body-fatness Higher levels of PAEE were linearly associated with lower BF% and FMI in both sexes (Figure 1); the adjusted beta coefficients (95%CI) being −0.16 (−0.17; −0.15) %BF per kJ·day-1·kg-1 in women and −0.09 (−0.10; −0.08) %BF per kJ·day-1·kg-1 in men. For FMI, corresponding beta coefficients were −0.07 (−0.08; −0.07) kg·m-2 per kJ·day-1·kg-1 in women and −0.03 (−0.04; −0.03) kg·m-2 per kJ·day-1·kg-1 in men. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F1) Figure 1: Sex-stratified binscatter of PAEE (kJ/day/kg) and body fat percentage and fat mass index, adjusted for all socio-demographic and dietary covariates. Each bin represents the mean values of 5% of the sex-stratified cohort. Women = grey diamonds, men = black circles. The relative intensity composition of the accumulated PAEE was also a significant factor in the association with body-fatness. Figure 2 shows the linear beta coefficients for the association with fraction of PAEE accumulated above increasingly higher intensity thresholds, whilst controlling for total PAEE; the association was steeper at higher intensities as indicated by progressively larger negative beta coefficients. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F2) Figure 2: Sex-stratified plot of the beta coefficients from 11 separate, sequential, multivariable linear regressions of the fraction of PAEE spent above each intensity (MET) threshold. Women = light grey diamonds, men = dark grey circles. Error bars are 95% confidence intervals. ### Isocaloric substitution analysis Table 2 shows the linear isocaloric substitution analyses for the percentage of energy expended at SS, LPA, MPA, or VPA, stratified by sex and adjusted for all dietary and socio-demographic covariates. These integrated analyses model the effect of overall PAEE on body-fatness, while simultaneously simulating the effect of substituting PAEE spent in any one intensity for another. In both sexes, higher levels of PAEE were associated with lower BF%, independent of the underlying intensity distribution. This association was stronger in women than men; each additional 1 kJ·day-1·kg-1 of active energy expended was associated with a 0.12 percentage-points lower body fat in women, compared to 0.05 percentage-points in men. For a given PAEE volume, the isocaloric substitutions of LPA for energy expended in MPA or VPA, or MPA for energy expended in VPA, were associated with a significantly lower body-fatness. The estimated association was stronger in women than in men (Women: LPA to VPA −0.26 %BF per 1 percentage-point PAEE, MPA to VPA −0.24 %BF per 1 percentage-point PAEE; Men: LPA to VPA −0.17 %BF per 1 percentage-point PAEE, MPA to VPA −0.12 %BF per 1 percentage-point PAEE). These results were not materially different when modelling FMI instead of BF% (Supplementary Table 1). Redefining the intensity categories to LPA (1.5 – 4 METs), MPA (4 – 7 METs), and VPA (>7 METs), produced statistically similar results, with a trend towards increased beta-coefficients for LPA to VPA substitutions (Supplementary Table 2). View this table: [Table 2:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T2) Table 2: Isocaloric substitution of physical activity energy expenditure and body fat percentage. To further examine the relative contribution of intensity and volume, we stratified analyses by PAEE tertile. This stratification revealed a larger linear PAEE beta-coefficient in first tertile women compared to second and third (−0.27, −0.11, and −0.12). Universally across all PAEE tertiles, reallocation of energy to VPA from any other intensity was associated with a significantly lower body fat percentage (Supplementary Table 3). Using compositional data analysis of the integrated volume-intensity association with body-fatness, overall PAEE remained significantly inversely associated with fatness in both women and men (Supplementary Table 4). The four compositional intensity reallocation curves for each sex are shown in Figure 3, with the accompanying boxplots indicating the relative size of the respective energy reservoirs. These curves show men and women had a similar estimated difference in BF% resulting from an isocaloric 1 kJ·day-1·kg-1 reallocation of active energy to a specific activity intensity, drawing proportionately from the three other intensities in the composition. Reallocation to VPA was significantly associated with lower body-fatness, whereas reallocation to LPA was significantly associated with higher body-fatness, and reallocation to MPA was not significant. However, reallocation to SS differed between men and women; in men, it was associated with lower body-fatness, whereas in women this association was non-significant. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F3) Figure 3: Estimated difference in body fat percentage associated with the isocaloric reallocation of PAEE to different intensities and box plots of the distribution of the PAEE composition, stratified by sex. The top panel shows the relative estimated difference in body fat percentage associated with an isocaloric reallocation of energy proportionately from all behaviours to the intensity of interest, as modelled by compositional data analysis. The origin (x=0,y=0) represents no change in the intensity composition of PAEE from the mean composition of the group of interest (women and men). The bottom panel illustrates the relative size of each reservoir of energy across women and men. Group PAEE values are mean (SD). SS = 0-1.5 METs, LPA = 1.5-3 METs, MPA = 3-6 METs, VPA >6 METs. To model the specific reallocation of energy from one specific intensity to another, we also conducted pairwise compositional analysis. This analysis is presented in Figure S1, alongside a graphical depiction of the linear substitution model. The regression results underpinning these curves can be found in Supplementary Table 5. To further explore the relative importance of intensity and volume, we conducted a sensitivity analysis of the compositional analysis, with intensity categories redefined as above. This analysis showed a slightly greater magnitude of association for MPA energy in women, but the association was non-significant in men as before. In contrast, the estimated difference in BF% per kJ·day-1·kg-1 active energy reallocated to VPA was greater (Figure S2, Supplementary Table 6). The relative importance of volume and composition is highlighted in Figure 4 (coefficients in Supplementary Table 7), which stratified the relationship between compositional intensity and body fat percentage by tertile of PAEE. On average, body fat percentage in the bottom vs top tertile of PAEE was 41% vs 34% in women and 30% vs 25% in men. The estimated difference in body fatness from reallocating PAEE to VPA was preserved across tertiles but was comparatively much smaller than the estimated effect of the absolute PAEE level itself. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F4) Figure 4: Estimated difference in body fat percentage associated with the isocaloric reallocation of PAEE to different intensities, stratified by sex and tertile of PAEE, as modelled by compositional data analysis. The origin (x=0,y=0) represents no change in the intensity composition of PAEE from the mean composition of the group of interest. Comparing groups by tertile at the intercept highlights the associated difference in body fat percentage attributable to varying levels of PAEE. Group PAEE values are mean (range). SS = 0-1.5 METs, LPA = 1.5-3 METs, MPA = 3-6 METs, VPA >6 METs. For comparison, equivalent isotemporal associations between physical activity intensity and body composition are shown in Supplementary Figures S3 and S4. The regressions underpinning these graphs are found in Supplementary Tables 8, 9 and 10. These models estimate the effect of intensity reallocation in the time domain and are not controlled for overall PAEE volume. ## Discussion In this large population-based cohort study of 11 468 free-living women and men with individually-calibrated, objectively-assessed physical activity and DXA-derived body composition, we show that overall physical activity volume is strongly inversely associated with body-fatness. Additionally, our integrated volume-and-intensity analyses demonstrate that higher proportions of more intense physical activity are associated with lower body-fatness, compared to when a similar volume of activity is accumulated at a lower intensity. The association between overall PAEE and body-fatness is in keeping with other population studies that highlight the important role that PA plays in maintaining a healthy body weight21–23. Our study confirms this association across the middle years of adult life. Uniquely, we have also shown the relative importance of PAEE compared to intensity; this is particularly prominent in our analysis of intensity composition, stratified by PAEE tertile. These results show that although more vigorous intensity activity is associated with lower body-fatness, the relative importance of intensity is smaller compared to the overarching contribution of PAEE volume. Indeed, all isocaloric associations must be considered in the context of overall PAEE. This is particularly important when considering the relative contribution of LPA, which was associated with higher body-fatness. In isolation, this result may appear paradoxical but when interpreted as a nested exposure within overall PAEE volume, it is not. Given the contribution of LPA to overall PAEE, especially in women11, and the challenges that some people face in being active, for example due to advancing age, ailments or other factors24, the value of LPA should not be discounted. Indeed, assuming associations are causal, our results suggest that were a woman to increase her PAEE from the 1st to the 2nd tertile, an average increase of 17 kJ·day-1·kg-1, her BF% would be roughly 3 percentage points lower. This contrasts to approximately 1 percentage point lower BF% if her PAEE remained fixed, but she accumulated an additional 30% of that PAEE from VPA (the maximum percentage observed within this group). If she were to increase both volume and VPA as above, the combined estimated benefit would be 4.5 percentage points. This is somewhat in contrast to the isotemporal substitution results, which may be perceived to have more direct relevance to the type of public health messaging which focuses on people’s time budget. These generally show that more time spent in higher intensity categories is associated with lower body fat. However, in the isotemporal analysis, the reallocation of time to higher intensity activity is associated with an inherent increase in volume. Thus, it is not possible to disentangle the estimated effects of volume and intensity. In addition, isotemporal analysis results also do not account for intensity variations within each intensity category, e.g. in our present analyses a 30-min activity at 3.1 METs is treated the same as a 30-min activity at 5.9 METs. In the integrated volume-and-intensity model in the energy space, this problem is effectively avoided as these two activities will result in substantively different PAEE levels as captured in the volume variable. This analysis acknowledges the nested nature of the exposures and aims to provide a simple answer to the primary question of whether overall physical activity is related to the outcome of interest, and secondarily what sub-dimensions of physical activity may or may not be associated beyond that. Apart from the more biologically intuitive specification of such isocaloric models, the public health relevance of these models can be interpreted by recasting the issue of potential behaviour change away from time to considering life as a series of daily tasks. For example, a common daily task for most working people is their commute to work. Whilst desirable, it may not be possible for everyone to switch to greener and more active modes of transport for the entire journey, but it would be possible for most car drivers and people who use public transport to incorporate more walking by parking further away or getting off the bus one stop earlier. This would increase activity volume but would also increase pressure on the time budget unless also accompanied by an increase in intensity. Our results indicate clear benefits to body-fatness of both such behaviour changes. Were it not initially possible to increase activity volume this way, small upward adjustments to the intensity by which set tasks are undertaken, for example by walking a set distance a bit faster, would still yield some health benefits, and importantly, it would also create a time surplus. Over time this may be sufficient to also accommodate more significant volume changes for even greater health benefits. Conversely, any small changes in, for example, the built environment, which result in subtle intensity decreases of set physical tasks or indeed the abolition of some tasks altogether, may have negative consequences for population levels of obesity if the associations we report here are causal. Our observations from this large cross-sectional study of UK adults are supported by evidence from trials seeking to delineate the contribution of activity intensity independent of PAEE. In a systematic review and meta-analysis of 28 intervention studies, Keating et al concluded that calorically matched high-intensity interval training and moderate-intensity continuous training provided similar benefits for reducing body fat, despite limitations of small sample sizes and incomplete control for confounding caused by variations in diet and physical activity outside of the training sessions that were part of the intervention25. A large population study such as ours provides evidence for the habitual physical activity perspective and has the advantage of being able to identify real-life activity profiles that associate with lower body-fatness, whilst controlling statistically for a range of confounders including habitual diet. Our results, at least in the time domain, are also consistent with previous observations in older British adults employing a similar exposure estimation technique, in which a positive association for sedentary time and inverse associations for both LPA and MVPA on body-fatness were reported, with a stronger magnitude of association for MVPA21. Comparable associations were also reported in a small sample of older high-risk Spanish adults with metabolic syndrome using wrist accelerometry-derived physical activity and DXA-derived body-fatness as estimated using linear isotemporal substitution analysis26. Likewise, an examination of BMI and waist circumference outcomes in US adults and the role of time spent in different physical activity intensity categories, modelled with compositional analysis, demonstrated an inverse cross-sectional association between time in MVPA and those outcomes4. Further, a recent longitudinal compositional study in elderly, Central European women, demonstrated that the longitudinal reallocation of time from MVPA to sedentary behaviour was associated with increased BMI and body-fatness27. Our findings also line-up with recently reported associations between activity and all-cause mortality in the UK Biobank cohort5. These findings suggest that the health benefits of higher intensity of activity, while significant, were secondary to that of volume. Taken together, these findings support the recently released WHO guidelines, which, above all else, emphasise the notion that every move counts28. Our study has several strengths including the large population size, the strength of the objective assessment methods for the exposures and outcomes, and the ability to control for a range of potential confounders. A further strength is the similarity of results across a range of modelling techniques. It is limited by its cross-sectional nature and further longitudinal studies and randomised controlled trials are needed before definitive statements about causality can be made. ## Conclusion In this population-based study of objectively measured physical activity and body-fatness, our integrated analysis of activity volume and intensity show that total PAEE and the pattern of accumulation of PAEE were both significantly associated with body-fatness outcomes independent of one another. At similar levels of PAEE, a greater proportion of energy expended at a higher intensity is associated with lower body-fatness. However, this association is secondary in order of magnitude to that of overall volume. ## Data Availability Contact the MRC Epidemiology Unit Cambridge with any data request. ## Funding The Fenland study was funded by the Medical Research Council and the Wellcome Trust. The current work was supported by the Medical Research Council (S.B., K.Wi., T.S., K.We, P.C.D., grant number MC\_UU\_12015/3, MC\_UU\_00006/4), (S.G., grant number MC\_UU\_12015/4, MC\_UU\_00006/6), (N.J.W., grant number MC\_UU\_12015/1, MC\_UU\_00006/1), (N.G.F., grant number MC\_UU_12015/5, MC_UU_00006/3); the National Institute of Health Research Cambridge (NIHR) Biomedical Research Centre (K.We., E.D.L.R., S.B., N.G.F., and N.J.W., grant number IS-BRC-1215-20014); National Health and Medical Research Council of Australia Research Fellowship (P.C.D., grant number 1142685), and the Cambridge Trust and St Catharine’s College (T.L.). The funders had no role in the design, analysis or writing of this article. ## Authors’ contributions The authors contributed to the present manuscript as follows: Idea for analysis (T.L., S.B.); acquisition, analysis of raw physical activity data (K.We., S.B.); acquisition, analysis of raw anthropometry and DXA data (E.D.L.R.) epidemiological data analysis (T.L.); drafting of the manuscript (T.L.); revising work critically for important intellectual content (all authors); approval of the final version before submission (all authors). Chief Investigator (N.J.W.) and Principal Investigators (N.G.F., S.G., S.B.) of the Fenland Study. ## Competing interests The authors declare no competing financial interests. ## Supplement Legend View this table: [Supplementary Table 1:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T3) Supplementary Table 1: Isocaloric substitution of physical activity energy expenditure and fat mass index. View this table: [Supplementary Table 2:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T4) Supplementary Table 2: Isocaloric substitution of physical activity energy expenditure and body fat percentage using alternative intensity thresholds (MPA > 4 METs, VPA > 7METs). View this table: [Supplementary Table 3:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T5) Supplementary Table 3: Isocaloric substitution of physical activity energy expenditure and body fat percentage by tertile of PAEE. View this table: [Supplementary Table 4:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T6) Supplementary Table 4: Relationship between isocaloric z1 ILR coordinate and body fat percentage. View this table: [Supplementary Table 5:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T7) Supplementary Table 5: Relationship between isocaloric pairwise ILR coordinates and body fat percentage. View this table: [Supplementary Table 6:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T8) Supplementary Table 6: Relationship between isocaloric z1 ILR coordinate and body fat percentage (MPA > 4MET, VPA > 7METs). View this table: [Supplementary Table 7:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T9) Supplementary Table 7: Relationship between isocaloric z1 ILR coordinate and body fat percentage by tertile of PAEE. View this table: [Supplementary Table 8:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T10) Supplementary Table 8: Isotemporal substitution of physical activity and body fat percentage. View this table: [Supplementary Table 9:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T11) Supplementary Table 9: Relationship between isotemporal z1 ILR coordinate and body fat percentage. View this table: [Supplementary Table 10:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/T12) Supplementary Table 10: Relationship between isotemporal pairwise ILR coordinates and body fat percentage. ![Figure S1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F5.medium.gif) [Figure S1:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F5) Figure S1: Pairwise and linear isocaloric reallocation of PAEE from one intensity to another. The top panel shows the results of pairwise compositional analysis, whereas the bottom panel shows linear substitution analysis. Both models estimate the difference in body fat percentage per 1% of PAEE reallocated. Group PAEE values are mean (SD). SS = 0-1.5 METs, LPA = 1.5-3 METs, MPA = 3-6 METs, VPA >6 METs. ![Figure S2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F6.medium.gif) [Figure S2:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F6) Figure S2: Estimated difference in body fat % associated with the isocaloric reallocation of PAEE to different intensities and box plots of the distribution of the PAEE composition, stratified by sex. Intensity thresholds redefined as SS = 0-1.5 METs, LPA = 1.5-4 METs, MPA = 4-7 METs, VPA >7 METs. The top panel shows the relative estimated difference in body fat percentage associated with an isocaloric reallocation of energy proportionately from all behaviours to the intensity of interest, as modelled by compositional data analysis. The origin (x=0,y=0) represents no change in the intensity composition of PAEE from the mean composition of the group of interest (women and men). The bottom panel illustrates the relative size of each reservoir of energy across women and men. Group PAEE values are mean (SD). ![Figure S3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F7.medium.gif) [Figure S3:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F7) Figure S3: Estimated difference in body fat percentage associated with the isotemporal reallocation of time to different intensities, and box plots of the distribution of time by intensities, stratified by sex. The top panel shows the relative estimated difference in body fat percentage associated with the reallocation of time proportionately from all behaviours to the intensity of interest, as modelled by compositional data analysis. The origin (x=0,y=0) represents no change in the intensity composition of time from the mean composition of the group of interest (women and men). The bottom panel illustrates the relative size of each reservoir of time across women and men. SS = 0-1.5 METs, LPA = 1.5-3 METs, MPA = 3-6 METs, VPA >6 METs. ![Figure S4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/28/2021.03.26.21254409/F8.medium.gif) [Figure S4:](http://medrxiv.org/content/early/2021/03/28/2021.03.26.21254409/F8) Figure S4: Pairwise and linear isotemporal reallocation of time from one intensity to another. The top panel shows the results of pairwise compositional analysis, whereas the bottom panel shows linear substitution analysis. Both models estimate the percentage difference in body fat % per minute per day reallocated. SS = 0-1.5 METs, LPA = 1.5-3 METs, MPA = 3-6 METs, VPA >6 METs. ## Acknowledgements We are grateful to the Fenland Study participants for their willingness and time to take part. The Fenland Study also has a dedicated Patient and Public Involvement panel, who provided input on the acceptability of the study protocols and the means by which participant data confidentiality is ensured. We thank all members of the following teams at the MRC Epidemiology Unit responsible for practical aspects of the study: Study Coordination, Field Epidemiology, Anthropometry Team, Physical Activity Technical Team, IT, Data Management, and Statistics. We would like to specifically acknowledge Stefanie Hollidge and Richard Powell for their assistance in the processing of physical activity and DXA data for this study. ## Footnotes * Competing interests: The authors declare no competing financial interests. * Received March 26, 2021. * Revision received March 26, 2021. * Accepted March 28, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Smith AD, Crippa A, Woodcock J, Brage S. Physical activity and incident type 2 diabetes mellitus: a systematic review and dose–response meta-analysis of prospective cohort studies. Diabetologia 2016; 59: 2527–2545. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-016-4079-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27747395&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) 2. 2.Moore SC, Lee I-M, Weiderpass E, Campbell PT, Sampson JN, Kitahara CM et al. Association of Leisure-Time Physical Activity With Risk of 26 Types of Cancer in 1.44 Million Adults. Jama Intern Med 2016; 176: 816. 3. 3.Pearce M, Strain T, Kim Y, Sharp SJ, Westgate K, Wijndaele K, et al. Estimating physical activity from self-reported behaviours in large-scale population studies using network harmonisation: findings from UK Biobank and associations with disease outcomes. Int J Behav Nutr Phy 2020; 17: 40. 4. 4.Chastin SFM, Palarea-Albaladejo J, Dontje ML, Skelton DA. Combined Effects of Time Spent in Physical Activity, Sedentary Behaviors and Sleep on Obesity and Cardio-Metabolic Health Markers: A Novel Compositional Data Analysis Approach. Plos One 2015; 10: e0139984. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0139984&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26461112&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) 5. 5.Strain T, Wijndaele K, Dempsey PC, Sharp SJ, Pearce M, Jeon J et al. Wearable-device-measured physical activity and future health risk. Nat Med 2020;: 1–7. 6. 6.Dumuid D, Pedišić Ž, Stanford TE, Martín-Fernández J-A, Hron K, Maher CA et al. The compositional isotemporal substitution model: A method for estimating changes in a health outcome for reallocation of time between sleep, physical activity and sedentary behaviour. Stat Methods Med Res 2017; 28: 846–857. 7. 7.Dumuid D, Stanford TE, Martin-Fernández J-A, Pedišić Ž, Maher CA, Lewis LK et al. Compositional data analysis for physical activity, sedentary time and sleep research. Stat Methods Med Res 2017; 27: 3726–3738. 8. 8.Wijndaele K, White T, Andersen LB, Bugge A, Kolle E, Northstone K, et al. Substituting prolonged sedentary time and cardiovascular risk in children and youth: a meta-analysis within the International Children’s Accelerometry database (ICAD). Int J Behav Nutr Phy 2019; 16: 96. 9. 9.Li SX, Imamura F, Schulze MB, Zheng J, Ye Z, Agudo A et al. Interplay between genetic predisposition, macronutrient intake and type 2 diabetes incidence: analysis within EPIC-InterAct across eight European countries. Diabetologia 2018; 61: 1325–1332. 10. 10.Willett WC, Howe GR, Kushi LH. Adjustment for total energy intake in epidemiologic studies. Am J Clin Nutrition 1997; 65: 1220S–1228S. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czoxMDoiNjUvNC8xMjIwUyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzI4LzIwMjEuMDMuMjYuMjEyNTQ0MDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 11. 11.Lindsay T, Westgate K, Wijndaele K, Hollidge S, Kerrison N, Forouhi N, et al. Descriptive epidemiology of physical activity energy expenditure in UK adults (The Fenland study). Int J Behav Nutr Phy 2019; 16: 126. 12. 12.Watson LPE, Venables MC, Murgatroyd PR. An Investigation Into the Differences in Bone Density and Body Composition Measurements Between 2 GE Lunar Densitometers and Their Comparison to a 4-Component Model. J Clin Densitom 2017; 20: 498–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jocd.2017.06.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28756995&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) 13. 13.Brage S, Brage N, Franks PW, Ekelund U, Wareham NJ. Reliability and validity of the combined heart rate and movement sensor Actiheart. Eur J Clin Nutr 2005; 59: 1602118. 14. 14.Brage S, Ekelund U, Brage N, Hennings MA, Froberg K, Franks PW et al. Hierarchy of individual calibration levels for heart rate and accelerometry to measure physical activity. J Appl Physiol 2007; 103: 682–692. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/japplphysiol.00092.2006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17463305&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248459600038&link_type=ISI) 15. 15.Stegle O, Fallert SV, MacKay DJ, Brage S. Gaussian Process Robust Regression for Noisy Heart Rate Data. Ieee T Bio-med Eng 2008; 55: 2143–2151. 16. 16.Brage S, Brage N, Franks PW, Ekelund U, Wong M-Y, Andersen LB et al. Branched equation modeling of simultaneous accelerometry and heart rate monitoring improves estimate of directly measured physical activity energy expenditure. J Appl Physiol 2004; 96: 343–351. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/japplphysiol.00703.2003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12972441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186984400046&link_type=ISI) 17. 17.Brage S, Westgate K, Franks PW, Stegle O, Wright A, Ekelund U et al. Estimation of Free-Living Energy Expenditure by Heart Rate and Movement Sensing: A Doubly-Labelled Water Study. Plos One 2015; 10: e0137206. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0137206&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26349056&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) 18. 18.Bingham SA, Welch AA, McTaggart A, Mulligan AA, Runswick SA, Luben R et al. Nutritional methods in the European Prospective Investigation of Cancer in Norfolk. Public Health Nutr 2001; 4: 847–858. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1079/PHN2000102&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11415493&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173798200015&link_type=ISI) 19. 19.Bingham S. Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers. Int J Epidemiol 1997; 26: 137S – 151. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/26.suppl_1.S137&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9126542&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WQ90700015&link_type=ISI) 20. 20.Mekary RA, Willett WC, Hu FB, Ding EL. Isotemporal Substitution Paradigm for Physical Activity Epidemiology and Weight Change. Am J Epidemiol 2009; 170: 519–527. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwp163&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19584129&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268586800015&link_type=ISI) 21. 21.Bann D, Kuh D, Wills AK, Adams J, Brage S, Cooper R et al. Physical Activity Across Adulthood in Relation to Fat and Lean Body Mass in Early Old Age: Findings From the Medical Research Council National Survey of Health and Development, 1946–2010. Am J Epidemiol 2014; 179: 1197–1207. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwu033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24722997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F28%2F2021.03.26.21254409.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000336094200007&link_type=ISI) 22. 22.Brage S, Lindsay T, Venables M, Wijndaele K, Westgate K, Collins D et al. Descriptive epidemiology of energy expenditure in the UK: findings from the National Diet and Nutrition Survey 2008-15. Int J Epidemiol 2020; 49: 1007–1021. 23. 23.Ekelund U, Brage S, Franks PW, Hennings S, Emms S, Wong M-Y et al. Physical activity energy expenditure predicts changes in body composition in middle-aged healthy whites: effect modification by age. Am J Clin Nutrition 2005; 81: 964–969. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI4MS81Lzk2NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAzLzI4LzIwMjEuMDMuMjYuMjEyNTQ0MDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 24. 24.Bredland EL, Söderström S, Vik K. Challenges and motivators to physical activity faced by retired men when ageing: a qualitative study. Bmc Public Health 2018; 18: 627. 25. 25.Keating SE, Johnson NA, Mielke GI, Coombes JS. A systematic review and meta-analysis of interval training versus moderate-intensity continuous training on body adiposity. Obes Rev 2017; 18: 943–964. 26. 26.Galmes-Panades AM, Konieczna J, Varela-Mato V, Abete I, Babio N, Fiol M et al. Changes in physical activity, sedentary behaviour and body composition: longitudinal analysis in the PREDIMED)-Plus trial. undefined 2020. doi:rs.3.rs-44115/v1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=rs.3.rs-44115/v1&link_type=DOI) 27. 27.Pelclová J, Štefelová N, Dumuid D, Pedišić Ž, Hron K, Gába A et al. Are longitudinal reallocations of time between movement behaviours associated with adiposity among elderly women? A compositional isotemporal substitution analysis. Int J Obesity 2020; 44: 1–8. 28. 28.Bull FC, Al-Ansari SS, Biddle S, Borodulin K, Buman MP, Cardon G et al. World Health Organization 2020 guidelines on physical activity and sedentary behaviour. Brit J Sport Med 2020; 54: 1451–1462. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiYmpzcG9ydHMiO3M6NToicmVzaWQiO3M6MTA6IjU0LzI0LzE0NTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wMy8yOC8yMDIxLjAzLjI2LjIxMjU0NDA5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)