Accelerometer-Derived Total Activity Counts, Bouted Minutes of Moderate to Vigorous Activity, and Insulin Resistance: NHANES 2003–2006
ORIGINAL RESEARCH — Volume 13 — October 20, 2016
William R. Boyer, MSH; Dana L. Wolff-Hughes, PhD; David R. Bassett, PhD; James R. Churilla, PhD; Eugene C. Fitzhugh, PhD
Suggested citation for this article: Boyer WR, Wolff-Hughes DL, Bassett DR, Churilla JR, Fitzhugh EC. Accelerometer-Derived Total Activity Counts, Bouted Minutes of Moderate to Vigorous Activity, and Insulin Resistance: NHANES 2003–2006. Prev Chronic Dis 2016;13:160159. DOI: http://dx.doi.org/10.5888/pcd13.160159.
The objective of this study was to compare the associations of accelerometer-derived total activity counts per day and minutes of bouted moderate to vigorous physical activity (MVPA) with insulin resistance.
The sample included 2,394 adults (aged ≥20 y) from the 2003–2006 National Health and Nutrition Examination Survey. Time spent in MVPA, measured by using 2 cutpoints (≥2,020 counts/min [MVPA2,020] and ≥760 counts/min [MVPA760]), was calculated for bouts of at least 8 to 10 minutes. Total activity counts per day reflects the total amount of activity across all intensities. Insulin resistance was measured via the homeostatic model assessment of insulin resistance (HOMA-IR) and the quantitative insulin sensitivity check index (QUICKI). Two nested regression models regressed HOMA-IR and QUICKI, respectively, on minutes of bouted MVPA and total activity counts per day. We used an adjusted Wald F statistic to illustrate strength of association.
After adjustment for covariates, total activity counts per day was more strongly associated with both HOMA-IR (adjusted Wald F = 36.83 , P < .001) and QUICKI (adjusted Wald F = 29.44, P < .001) compared with MVPA2,020 (HOMA-IR, adjusted Wald F = 4.00, P = .06; QUICKI, adjusted Wald F = 1.08, P = .31).Total activity counts per day was more strongly associated with both HOMA-IR (adjusted Wald F = 13.64, P < .001) and QUICKI (adjusted Wald F = 12.10, P < .001) compared with MVPA760 (HOMA-IR, adjusted Wald F = 1.13, P = .30; QUICKI, adjusted Wald F = 0.97, P = .33).
Our study indicated that total activity counts per day has stronger associations with insulin resistance compared with minutes of bouted MVPA. The most likely explanation is that total activity counts per day captures data on light physical activity and intermittent MVPA, both of which influence insulin resistance.
Insulin resistance is associated with an increased risk for cardiovascular disease (1), diabetes (2), and other metabolic disorders. An estimated 30% of the US adult population is insulin resistant (3). Investigating ways to prevent or delay insulin resistance, researchers have focused on the role of physical activity. Studies have highlighted the importance of both total volume (4) and intensity (5) of physical activity in relation to insulin resistance. One study (5) found a significant difference in insulin response after 7 days of exercise at 70% of maximal oxygen consumption but no difference in those exercising at 50% of maximal oxygen consumption. In contrast, another study (4) found lower levels of insulin resistance independent of intensity among people exercising at 170 minutes per week compared with people exercising 115 minutes per week. Results of cross-sectional studies also indicated that higher levels of moderate to vigorous physical activity (MVPA), both bouted (6) and nonbouted (7), were associated with lower levels of insulin resistance.
Wearable monitors that can assess physical activity on a minute-by-minute basis (8,9) have improved researchers’ ability to examine the relationship between physical activity and various chronic conditions such as insulin resistance (6,7,10). Methods of analyzing accelerometer data allow researchers to classify minutes spent in light physical activity (LPA), moderate physical activity (MPA), and vigorous physical activity (VPA); several cutpoints have been established to measure time spent in these intensities (11). However, recent studies show the limitations of these methods in accurately quantifying time (8,12); for example, these methods miss upwards of 50% of the time spent in MVPA. Because of these limitations, a new metric for quantifying accelerometer-derived data is being used: total activity counts per day (TAC/d). This metric incorporates data on physical activity across all intensities, frequencies, and durations (9).
Only 2 recent studies show that accelerometer-derived measures of TAC/d are associated with insulin resistance (13,14). A study published in 2014 (13) showed that TAC/d was significantly associated with measures of insulin sensitivity derived from the homeostatic model assessment of insulin resistance (HOMA-IR). Furthermore, the use of TAC/d eliminated the association between sedentary time and HOMA-IR. A study published in 2015 of a population-based sample (14) found that TAC/d was more strongly associated with multiple insulin resistance–related risk factors — fasting blood glucose, fasting insulin, and C-peptide — than was bouted MVPA.
Although both bouted and nonbouted MVPA are associated with insulin resistance, the aforementioned studies indicate that TAC/d should be considered when assessing the association between physical activity and insulin resistance. Furthermore, it is necessary measure the total volume of activity, because subsets of total volume, such as MVPA, may not fully illustrate the relationship between physical activity and insulin resistance (9,13). To date, no studies have investigated the associations between objectively measured physical activity and insulin resistance by measuring both MVPA and total volume of physical activity. The objective of this study was to determine which measure — MVPA or TAC/d —is more strongly associated with 2 markers of insulin resistance commonly used in epidemiological research — HOMA-IR and the quantitative insulin sensitivity check index (QUICKI) — in a representative sample of US adults. The results of this study will provide further insight into the relationship of TAC/d and bouted MVPA with insulin resistance and the use of these 2 measures in epidemiologic research.
This study used data from the 2003–2006 National Health and Nutrition Examination Survey (NHANES), a cross-sectional survey conducted by the National Center for Health Statistics that uses a complex, multistage sampling design to obtain a representative sample of the US population (15). NHANES was designed to provide national estimates of the health and nutritional status of noninstitutionalized US civilians aged 2 months or older. Our analytic sample (n = 2,394) met the following conditions: 1) they were adult men and women aged 20 years or older; 2) they attended a morning mobile examination center (MEC) after an overnight fast of 8 to 9 hours; 3) if female, they were not pregnant or lactating; 4) they had complete data on all the variables of interest; and 5) they had at least 4 days of accelerometer data with 10 or more hours per day of wear time.
Accelerometer-derived physical activity
Eligible participants were those who participated in the MEC component of NHANES, were ambulatory, and agreed to participate. Eligible participants were given an ActiGraph model 7164 accelerometer and were instructed to wear the device on the right hip for 7 consecutive days (16). Because the ActiGraph is not waterproof, participants were instructed to remove the device during activities in water, such as bathing or swimming (16). Participants were instructed to remove the device also while sleeping.
Accelerometer data were recorded in 1-minute epochs. Nonwear time was defined as 60 minutes or more of no accelerometer counts, with 2 minutes or fewer of limited movement between zero and 100 counts. Wear time was calculated by subtracting the total nonwear time from 24 hours. Participants who had at least 4 days with 10 hours or more hours of valid wear time per day were included in the analysis (16). Several cutpoints exist to denote minutes of MVPA (11). For our study, we used 2 cutpoints for bouted MVPA. One cutpoint (MVPA2,020) was defined as minutes consisting of 2,020 or more activity counts that were accumulated in bouts of at least 10 minutes (with 1–2 minutes allowed below the threshold) (17). A second cutpoint (MVPA760) was defined as minutes consisting of 760 or more activity counts that were accumulated in bouts of at least 10 minutes (with 1–2 minutes allowed below the threshold) (18). These cutpoints are consistently applied to NHANES data (11). TAC/d was created by summing up the total activity counts and dividing by the number of valid wear days (14).
Insulin resistance measures were HOMA-IR (19) and QUICKI (20), both of which are used in NHANES data analyses (21,22). The HOMA-IR equation was created as a surrogate marker of insulin resistance based on the negative-feedback loop between the liver and β cells (19). QUICKI is a measure of insulin sensitivity, with lower levels reflecting a greater degree of insulin resistance. The QUICKI method is an equation derived from fasting insulin and glucose levels (similar to HOMA-IR) and is based on the linear relationship of these variables with the gold-standard hyperinsulinemic–euglycemic clamp technique (20). Before calculating HOMA-IR and QUICKI, we corrected for differences between the 2005–2006 NHANES cycle and the 2003–2004 NHANES cycle in measuring insulin and glucose values (23). The assays used to measure insulin and glucose changed between 2003–2004 and 2005–2006. NHANES provides an equation to account for those differences, and it was applied before analyzing the data. Data on HOMA-IR were logarithmically transformed, because they were nonnormally distributed.
We included sociodemographic variables, wear-time variables, and markers of health and behavior as covariates. Sociodemographic variables were age (continuous variable in years), sex (male and female), race/ethnicity (non-Hispanic white, non-Hispanic black, Mexican American, and other), and education (
All analyses were conducted using SAS 9.4 (SAS Institute, Inc). We used the proc surveyreg procedure, which accounts for the complex sampling design and survey nonresponse inherent to NHANES (15). Fasting subsample weights were used, thus ensuring the proper subsample for all analyses. Weighted prevalence estimates were generated for all categorical variables. The weighted mean, standard error, and range were generated for all continuous variables.
Age-adjusted linear regressions were performed; these regressed TAC/d and MVPA separately on each measure of insulin resistance. Two nested linear regression models were created. Both physical activity metrics were simultaneously entered into the regression models to examine the difference in the magnitude of the association between each physical activity metric and insulin resistance. Model 1 included all covariates except waist circumference. Model 2 included all covariates from model 1 plus waist circumference. The Wald F statistic was used to indicate the strength of association for both physical activity metrics. A larger Wald F indicates a stronger association with the dependent variable (26). For the regression models using MVPA2,020 and TAC/d, multicollinearity was checked and not found to violate the variance inflation factor threshold (< 5) (27,28). In the regression models using MVPA760 and TAC/d the variance inflation factor statistic was 5.7, violating the threshold set for the analysis. A variance inflation factor of less than 5 is a conservative statistic, and other studies have used a value of less than 10 (29).
The mean age of our study sample was 48.3 years; 51.4% were women, 73.5% were non-Hispanic white, 32.1% had some college, 51.4% were nonsmokers, 26.9% had hypertension, and 6.5% had diabetes (Table 1). The age-adjusted linear regressions (Table 2) showed that MVPA2,020 was significantly associated with both HOMA-IR (β = −0.01, P < .001) and QUICKI (β = 0.0003, P < .001). MVPA760 was also significantly associated with both HOMA-IR (β = −0.003, P < .001) and QUICKI (β = 0.00001, P < .001). TAC/d was also significantly associated with both HOMA-IR (β = −0.000002, P < .001) and QUICKI (β = 0.00000004, P < .001).
The nested regression analyses that used both MVPA2,020 and TAC/d in the same regression model and adjusted for covariates (excluding waist circumference) showed that TAC/d was more strongly associated with HOMA-IR (adjusted Wald F = 58.97, P < .001) and QUICKI (adjusted Wald F = 46.71, P < .001) than was MVPA2,020 (HOMA-IR, adjusted Wald F = 0.62, P = .44; QUICKI, adjusted Wald F = 1.10, P = .30) (Table 3). After adjustment for waist circumference, TAC/d remained a stronger predictor of HOMA-IR (adjusted Wald F = 36.83, P < .001) and QUICKI (adjusted Wald F = 29.44, P < .001) compared with MVPA2,020 (HOMA-IR, adjusted Wald F = 4.00, P = .06; QUICKI, adjusted Wald F = 1.08, P = .31).
The nested regression analyses that used MVPA760 and TAC/d and adjusted for covariates (excluding waist circumference) found that TAC/d was more strongly associated with HOMA-IR (adjusted Wald F = 31.59, P < .001) and QUICKI (adjusted Wald F = 28.00, P < .001) than was MVPA760 (HOMA-IR, adjusted Wald F = 2.18, P = .15; QUICKI, adjusted Wald F = 1.93, P = .17) (Table 3). After adjustment for waist circumference, TAC/d remained a stronger predictor of HOMA-IR (adjusted Wald F = 13.64, P < .001) and QUICKI (adjusted Wald F = 12.10, P < .001) compared with MVPA760 (HOMA-IR, adjusted Wald F = 1.13, P = .30; QUICKI, adjusted Wald F = 0.97, P = .33).
The results of this study provide insight into the associations between objectively measured physical activity and insulin resistance. In separate models, both MVPA and TAC/d were significantly associated with insulin resistance. However, the combined models showed that TAC/d was more strongly associated with insulin resistance than was MVPA. Moreover, regardless of the cutpoints used, bouted minutes of MVPA were not significantly associated with insulin resistance in the fully adjusted models.
Our results, which used surrogate measures of insulin resistance, extend the results reported by investigators (14) who found that TAC/d, compared with bouted MVPA, was more strongly associated with insulin resistance–related cardiometabolic risk factors (ie, fasting glucose, fasting insulin, and C-peptide). Using data from the 2003–2006 NHANES, another study (13) found significant associations between TAC/d and HOMA-derived insulin sensitivity or β cell function, independent of inactive time. These results, as well as the results of our study, indicate that TAC/d is important to consider, independent of MVPA, when examining physical activity and insulin resistance.
Two possible explanations exist for why TAC/d is more strongly associated with insulin resistance compared with bouted MVPA. First, TAC/d is not limited by having to classify each minute according to categorical cutpoints that define sedentary time, LPA, MPA, or VPA; instead, activity counts per minute vary along a continuous scale (9,12). Second, TAC/d incorporates the intensity, frequency, and duration of activity, while bouted MVPA represents only a small subset of the overall volume of physical activity. Several studies highlight the independent contribution of LPA (30) and nonbouted MVPA (7,10) in improving insulin resistance profiles. The use of TAC/d accounts for the activity attributable to these other physical activity subcategories, thus strengthening the relationship between physical activity and insulin resistance. TAC/d is more strongly related to cardiometabolic risk than MVPA (14), and our study extends this finding to insulin resistance.
A major strength of this study is the external validity resulting from the use of a large, nationally representative sample. Another strength is the use of validated indices of insulin resistance that are highly correlated with the gold standard euglycemic–hyperinsulinemic clamp technique (19,20). This study also has several limitations. First is the use of bouted MVPA minutes using the 2,020 counts-per-minute cutpoint, which can underestimate time spent in MVPA (12). However, these cutpoints have consistently been used in NHANES analyses (11). Second, although the 760 counts-per-minute cutpoint captures more MVPA, it can also overestimate it (12). Third, physical activity may be underestimated by accelerometry because accelerometers cannot capture data on nonambulatory movements such as swimming and resistance training. Fourth, the counts used in this study were specific to the features of the ActiGraph accelerometer. Thus, these counts cannot be applied to other devices. Fifth, the metric TAC/d has no intuitive meaning (9): although it provides insight into the total volume of physical activity, it does not have the same kind of immediate meaning as the number of minutes, for example. However, population-referenced TAC/d percentiles standardized by age and sex were published in 2015 (31); these percentiles were used to characterize differences in physical activity levels across population subgroups in New York City (32). In the future, criterion-referenced standards for TAC/d may further our understanding of the relationship between physical activity and chronic disease.
In our study of a nationally representative sample of US adults, TAC/d, an index of total activity volume, was more strongly associated with measures of insulin resistance and insulin sensitivity compared with minutes of bouted MVPA using the cutpoints of 2,020 counts per minute and 760 counts per minute. We speculate that the stronger association resulted from the ability of TAC/d to capture counts associated with LPA and nonbouted MVPA, which are both independently associated with insulin resistance (7,10). Future analyses should examine total MVPA defined by the cutpoints used in our study and other common cutpoints to further understand the relationship between components of physical activity and insulin resistance. Finally, our study indicates that it is important to consider physical activity across the entire intensity spectrum when investigating physical activity and insulin resistance.
This research was not funded.
Corresponding Author: William Boyer, MSH, Department of Kinesiology, Recreation, and Sports Studies, University of Tennessee, 1914 Andy Holt Ave, 307 HPER Bldg, Knoxville, TN 37996-2700. Telephone: 904-502-3473. Email: email@example.com.
Author Affiliations: Dana L. Wolff-Hughes, Division of Cancer Control and Population Sciences, National Cancer Institute, Bethesda, Maryland; David R. Bassett, Eugene C. Fitzhugh, Department of Kinesiology, Recreation and Sport Studies, University of Tennessee, Knoxville, Tennessee; James R. Churilla, Department of Clinical and Applied Movement Sciences, University of North Florida, Jacksonville, Florida.
- Bonora E, Formentini G, Calcaterra F, Lombardi S, Marini F, Zenari L, et al. HOMA-estimated insulin resistance is an independent predictor of cardiovascular disease in type 2 diabetic subjects: prospective data from the Verona Diabetes Complications Study. Diabetes Care 2002;25(7):1135–41. CrossRef PubMed
- Bonora E, Kiechl S, Willeit J, Oberhollenzer F, Egger G, Meigs JB, et al. Population-based incidence rates and risk factors for type 2 diabetes in white individuals: the Bruneck study. Diabetes 2004;53(7):1782–9. CrossRef PubMed
- Ioannou GN, Bryson CL, Boyko EJ. Prevalence and trends of insulin resistance, impaired fasting glucose, and diabetes. J Diabetes Complications 2007;21(6):363–70. CrossRef PubMed
- Houmard JA, Tanner CJ, Slentz CA, Duscha BD, McCartney JS, Kraus WE. Effect of the volume and intensity of exercise training on insulin sensitivity. J Appl Physiol (1985) 2004;96(1):101–6.CrossRef PubMed
- Kang J, Robertson RJ, Hagberg JM, Kelley DE, Goss FL, DaSilva SG, et al. Effect of exercise intensity on glucose and insulin metabolism in obese individuals and obese NIDDM patients. Diabetes Care 1996;19(4):341–9. CrossRef PubMed
- Balkau B, Mhamdi L, Oppert JM, Nolan J, Golay A, Porcellati F, et al. Physical activity and insulin sensitivity: the RISC study. Diabetes 2008;57(10):2613–8. CrossRef PubMed
- Nelson RK, Horowitz JF, Holleman RG, Swartz AM, Strath SJ, Kriska AM, et al. Daily physical activity predicts degree of insulin resistance: a cross-sectional observational study using the 2003–2004 National Health and Nutrition Examination Survey. Int J Behav Nutr Phys Act 2013;10(1):10. CrossRef PubMed
- Chen KY, Bassett DR Jr. The technology of accelerometry-based activity monitors: current and future. Med Sci Sports Exerc 2005;37(11, Suppl):S490–500. CrossRef PubMed
- Bassett DR, Troiano RP, McClain JJ, Wolff DL. Accelerometer-based physical activity: total volume per day and standardized measures. Med Sci Sports Exerc 2015;47(4):833–8. CrossRefPubMed
- Ekelund U, Brage S, Griffin SJ, Wareham NJ; ProActive UK Research Group. Objectively measured moderate- and vigorous-intensity physical activity but not sedentary time predicts insulin resistance in high-risk individuals. Diabetes Care 2009;32(6):1081–6. CrossRef PubMed
- Tudor-Locke C, Camhi SM, Troiano RP. A catalog of rules, variables, and definitions applied to accelerometer data in the National Health and Nutrition Examination Survey, 2003–2006. Prev Chronic Dis 2012;9:E113. PubMed
- Crouter SE, DellaValle DM, Haas JD, Frongillo EA, Bassett DR. Validity of ActiGraph 2-regression model, Matthews cut-points, and NHANES cut-points for assessing free-living physical activity. J Phys Act Health 2013;10(4):504–14. CrossRef PubMed
- Maher C, Olds T, Mire E, Katzmarzyk PT. Reconsidering the sedentary behaviour paradigm. PLoS One 2014;9(1):e86403. CrossRef PubMed
- Wolff-Hughes DL, Fitzhugh EC, Bassett DR, Churilla JR. Total activity counts and bouted minutes of moderate-to-vigorous physical activity: relationships with cardiometabolic biomarkers using 2003–2006 NHANES. J Phys Act Health 2015;12(5):694–700. CrossRef PubMed
- National Center for Health Statistics. Centers for Disease Control and Prevention. NHANES analytic and reporting guidelines; 2006. http://www.cdc.gov/nchs/data/nhanes/nhanes_03_04/nhanes_analytic_guidelines_dec_2005.pdf. Accessed November 20, 2015.
- Troiano RP, Berrigan D, Dodd KW, Mâsse LC, Tilert T, McDowell M. Physical activity in the United States measured by accelerometer. Med Sci Sports Exerc 2008;40(1):181–8. CrossRefPubMed
- National Health and Nutrition Examination Survey (NHANES) 2003–2004 documentation, codebook, and frequencies. MEC exam component: physical activity monitor examination data. Atlanta (GA): Centers for Disease Control and Prevention; 2007. https://wwwn.cdc.gov/Nchs/Nhanes/2003-2004/PAXRAW_C.htm. Accessed November 15, 2015.
- Matthew CE. Calibration of accelerometer output for adults. Med Sci Sports Exerc 2005;37(11, Suppl):S512–22. CrossRef PubMed
- Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 1985;28(7):412–9. CrossRef PubMed
- Katz A, Nambi SS, Mather K, Baron AD, Follmann DA, Sullivan G, et al. Quantitative insulin sensitivity check index: a simple, accurate method for assessing insulin sensitivity in humans. J Clin Endocrinol Metab 2000;85(7):2402–10. CrossRef PubMed
- Cheng YJ, Gregg EW, De Rekeneire N, Williams DE, Imperatore G, Caspersen CJ, et al. Muscle-strengthening activity and its association with insulin sensitivity. Diabetes Care 2007;30(9):2264–70. CrossRef PubMed
- Boyer WR, Johnson TM, Fitzhugh EC, Richardson MR, Churilla JR. The associations between increasing degrees of homeostatic model assessment for insulin resistance and muscular strengthening activities among euglycaemic US adults. Diab Vasc Dis Res 2015;12(6):420–7. CrossRef PubMed
- National Health and Nutrition Examination Survey (NHANES) 2005–2006 documentation, codebook, and frequencies. Plasma fasting glucose and insulin. Atlanta (GA): Centers for Disease Control and Prevention; 2008. https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/GLU_D.htm. Accessed November 15, 2015.
- James PA, Oparil S, Carter BL, Cushman WC, Dennison-Himmelfarb C, Handler J, et al. 2014 evidence-based guideline for the management of high blood pressure in adults: report from the panel members appointed to the Eighth Joint National Committee (JNC 8). JAMA 2014;311(5):507–20. CrossRef PubMed
- American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care 2014;37(Suppl 1):S81–90. CrossRef PubMed
- Clogg CC, Petkova E, Haritou A. Statistical methods for comparing regression coefficients between models. Am J Sociol 1995;100(5):1261–93. CrossRef
- O’Brien R. A caution regarding rules of thumb for variance inflation factors. Qual Quant 2007;41(5):673–90. CrossRef
- Hair JF, Black WC, Babin BJ, Anderson RE. Multivariate data analysis. Saddle River (NJ): Prentice Hall; 2009.
- Belsley DA. Conditioning diagnostics: collinearity and weak data in regression. New York (NY): Wiley-Interscience; 1991.
- Herzig KH, Ahola R, Leppäluoto J, Jokelainen J, Jämsä T, Keinänen-Kiukaanniemi S. Light physical activity determined by a motion sensor decreases insulin resistance, improves lipid homeostasis and reduces visceral fat in high-risk subjects: PreDiabEx study RCT. Int J Obes 2014;38(8):1089–96. CrossRef PubMed
- Wolff-Hughes DL, Fitzhugh EC, Bassett DR, Churilla JR. Waist-worn actigraphy: population-referenced percentiles for total activity counts in U.S. adults. J Phys Act Health 2015;12(4):447–53.CrossRef PubMed
- Wolff-Hughes DL, Troiano RP, Boyer WR, Fitzhugh EC, McClain JJ. Use of population-referenced total activity counts percentiles to assess and classify physical activity of population groups. Prev Med 2016;87:35–40. CrossRef PubMed