OBJECTIVE

Intensive health behavior and lifestyle treatment (IHBLT) is recommended for children aged 6–18 years with obesity. The objective was to evaluate the effectiveness of Fit Together, a health care and parks and recreation partnership to deliver IHBLT.

METHODS

A randomized controlled trial was conducted from 2018 to 2021. Youths (aged 5–17 years) with obesity were recruited from primary care clinics and randomized to a waitlist control or Fit Together (ie, clinical obesity care plus group-based lifestyle sessions at a local recreation center). Primary outcomes, child body mass index relative to the 95th percentile (BMIp95) and submaximal heart rate, were collected at baseline and 6 months. Generalized estimating equation models were used to assess changes in primary outcomes for those affected and not affected by COVID-19 study disruptions.

RESULTS

Participants (n = 255) had a mean (SD) age of 10.0 (3.0) years, were 39% Hispanic, and were 38% non-Hispanic Black. Intervention youths not affected by COVID-19 disruptions experienced a significant decrease in BMIp95 (β = −3.05; 95% confidence interval [CI], −5.08 to −1.01) compared with controls. There was no difference in BMIp95 between intervention and control youths affected by COVID-19 disruptions (β = −3.25; 95% CI, −7.98 to 1.48). For the entire cohort, intervention youths had a significant decrease in BMIp95 compared with control youths (β = −3.32; 95% CI, −5.69 to −0.96). Submaximal heart rate was only available for the nondisrupted group, but there was no difference between intervention and control youths (β = −7.18; 95% CI, −16.12 to 1.76).

CONCLUSION

Effective child obesity treatment can be implemented in local communities through a partnership between clinical practices and parks and recreation providers. Future research will explore this model in combination with newer, more effective obesity treatments.

What’s Known on This Subject:

The AAP Clinical Practice guideline recommends intensive health behavior and lifestyle treatment for pediatric obesity. Access to treatment is low, particularly for marginalized populations. A major barrier for health care settings is meeting the at least 26-hour threshold needed to demonstrate treatment effectiveness.

What This Study Adds:

This randomized controlled trial demonstrated the feasibility and effectiveness of a partnership between health care and parks and recreation, as compared with usual care, on body mass index among diverse youths aged 5-17 years with obesity.

Pediatric obesity is a complex, multifactorial disease with genetic, socioeconomic, and environmental influences.1 As of 2020, 14.7 million United States youths aged 2–19 years were affected by obesity,2,3 and the current obesity prevalence is almost certainly higher as there was an observed doubling of the rate of youth body mass index (BMI) increase during the COVID-19 pandemic.4 Obesity is unlikely to resolve without treatment and leads to significant health and social consequences.5–7 The 2023 American Academy of Pediatrics (AAP) Clinical Practices Guideline for the Evaluation and Treatment of Children and Adolescents with Obesity outlines safe and effective treatment options, including intensive health behavior and lifestyle treatment (IHBLT), pharmacotherapy, and metabolic and bariatric surgery.8 IHBLT is the cornerstone of treatment, providing intensive nutrition and physical activity opportunities for children and their families (ie, ≥26 hours of in-person contact over 3–12 months).9 However, IHBLT is limited by poor access, particularly in underserved communities, and the limited capacity of clinics alone to deliver the recommended treatment intensity.10 Additionally, little is known about the effectiveness of IHBLT among low-income and racial and ethnic minority groups.

One proposed solution has been for health care systems to partner with community organizations to increase capacity and improve access to IHBLT.11 Public parks and recreation (P&R) centers provide low-cost or free access to physical activity and are accessible to diverse populations. Access to P&R facilities is associated with higher levels of physical activity among children and adults.12,13 We have previously developed a model of pediatric obesity treatment called Fit Together that pairs health care organizations with local P&R providers to deliver IHBLT.14,15 In the Fit Together model, health care partners screen for obesity, provide counseling, and treat comorbidities while P&R partners provide IHBLT in the community setting. In pilot testing, Fit Together has demonstrated feasibility in engaging a diverse population of youths, the ability to meet the 26-hour threshold, and effectiveness in improving self-reported physical activity and quality of life.15,16 Based on these preliminary data, we sought to test Fit Together in a randomized controlled trial (Hearts & Parks). Therefore, the aim of the current study was to compare the Fit Together model with clinical care alone to determine the effectiveness of Fit Together on youth BMI and cardiorespiratory fitness among a diverse sample of youths with obesity.

Hearts & Parks was a 2-arm, randomized controlled crossover trial that aimed to evaluate the effectiveness of the Fit Together model of youth obesity treatment in comparison with routine primary care obesity management (ClinicalTrials.gov ID: NCT03339440). A detailed description of the design and methods is available elsewhere.17 This study and all procedures were approved by the Duke University Health System Institutional Review Board (Pro00086684). The Consolidated Standards of Reporting Trials reporting guidelines were used in this study (Supplementary File 1).18 

Participants were recruited from four primary care clinics in Durham, North Carolina, from 2018 to 2020. The eligibility criteria included (1) age 5–17 years and BMI at least at the 95th percentile for age and sex at the time of referral, (2) youth and caregiver ability to speak/read English or Spanish, and (3) caregiver ownership of a device capable of sending and receiving text messages. The exclusion criteria included (1) families living more than 20 miles from the community program, (2) endogenous or genetic causes of obesity, (3) taking medication that causes weight gain or loss, (4) participation in an obesity treatment program within the previous 12 months, or (5) either the caregiver or child having a significant health problem that would limit participation. During well-child visits, primary care physicians (PCPs) were alerted to potentially eligible participants by a best-practice alert (BPA) in the electronic health record. The BPA was programmed specifically for the study and would pop up if a patient met the inclusion criteria. PCPs then had the option to refer the patient to the study (on family consent), opt the patient out (eg, if they had severe obesity that would prevent them from participating or urgent comorbidities), or decline on behalf of the family. Trained study members confirmed eligibility, contacted families regarding participation, and proceeded with consent/assent before the baseline study visit.

Fit Together is a 6-month intervention that includes clinical visits with a pediatric clinic and lifestyle sessions at a local P&R center. In Hearts & Parks, the Fit Together intervention included clinical visits at the Duke Healthy Lifestyles pediatric obesity treatment clinic and access to Bull City Fit, a community program delivered at a Durham P&R center.14 At Healthy Lifestyles, participants received care from a team of medical providers, registered dietitians, physical therapists and a licensed counselor. Providers worked with patients and families to set and monitor health-related goals and monitor patient safety. Patients were encouraged to visit Healthy Lifestyles once a month. Visits with each provider lasted 30 minutes, and patients often saw multiple providers on the same day. Participants were also encouraged to attend Bull City Fit at least twice a week. Bull City Fit was offered 6 days a week and, at each session, there were opportunities for unstructured and structured physical activity. A weekly cooking class and nutrition education session was also offered. To promote family engagement, at least 1 adult household member was required to attend and participate with the child, and siblings were encouraged to join.15 

The waitlist control group received nutrition and activity counseling from their PCP during standard well-child visits and at additional follow-up visits at the discretion of the PCP. All PCPs at the 4 practices received training on motivational interviewing and written materials on nutrition and activity counseling. Families also received a non–obesity-related literacy intervention consisting of a monthly $15 gift card to a locally owned bookstore. At the end of the 6-month waiting period, they were invited to participate in the Fit Together program.

The target sample size was 270 participants based on previous research showing that a meaningful difference in mean 6-month BMI z score is −0.10.7 This sample size calculation accounted for a 30% loss to follow-up rate at 6 months and was computed using a 2-sided t test at an alpha significance level of 0.05.

Participants were randomized to the intervention or comparison group in a 1∶1 allocation ratio. Randomization was conducted using the Research Electronic Data Capture system by the study biostatistician. Research coordinators notified patients of their assigned group following the baseline visit. Given the nature of the intervention, neither patients nor providers could be blinded to the treatment assignment.

Data were collected at baseline and at 6 months at the clinic or community site. The study was powered to detect differences in 2 primary outcomes: youth BMI and cardiorespiratory fitness.17 

Height was assessed using a Seca stadiometer (Seca, Hamburg, Germany), and weight was assessed using a Tanita body composition scale (Tanita, Tokyo, Japan) using standard protocols. Measures were taken in triplicate and averaged. BMI was calculated based on age- and sex-specific reference values from the Centers for Disease Control and Prevention growth charts.19,20 Because of the known challenges in assessing changes in BMI among youths with obesity, height and weight measures were converted to BMI relative to the 95th percentile (BMIp95).21 Additionally, the following classes of obesity were calculated: class I obesity (BMI ≥95th percentile or BMI ≥30), class II obesity (BMI ≥120% of the 95th percentile or BMI ≥35), and class III obesity (BMI ≥140% of the 95th percentile or BMI ≥40).

Cardiorespiratory fitness was assessed using a submaximal 3-minute step test, a valid and reliable test of cardiorespiratory fitness in youths.22 A baseline heart rate (HR) measure was recorded before the start of the test using a pulse oximeter. Using a metronome, participants were then asked to step up and down on a platform for 3 minutes at a pace of 96 beats per minute. Another HR measure was taken immediately following the completion of the stepping protocol as a measure of peak HR (ie, 3-minute HR). Then the participant was directed to sit in a chair for 1 minute, where their HR was taken again to assess submaximal HR (ie, 4-minute HR). The 4-minute HR values were included in the analysis.

For intervention participants, attendance at clinic visits and Bull City Fit sessions were recorded to assess if participants were able to achieve the recommended 26 hours of treatment. Each visit to Healthy Lifestyles was counted as 1 engagement hour, and each Bull City Fit visit was counted as 2 engagement hours. Total engagement was calculated as the sum of all hours during the first 6 months of the study. Engagement hours were classified into the following categories: 0, less than 5, 5–25, or more than 26 hours.

At baseline, caregivers completed a demographic questionnaire that included child sex, race, ethnicity, annual family income, caregiver employment, food insecurity, and language spoken at home.

On March 15, 2020, all in-person study visits were halted because of COVID-19 mitigation measures. At this point, enrollment was concluded; however, to complete follow-up study visits, a remote data-collection protocol was established. Digital scales and stadiometers were mailed to participants’ homes. Study visits were conducted through direct video observation of child height and weight measures. The 3-minute step test was discontinued because of the inability to complete it remotely. The intervention was adapted to an online format, in which Healthy Lifestyles visits were conducted by telehealth, and Bull City Fit sessions were moved to an online platform, which included synchronous Zoom sessions or asynchronous recordings.

Descriptive statistics were computed to summarize sample characteristics at baseline. Means and SDs were calculated for all continuous variables and frequencies (%) for categorical variables. Differences among the groups at baseline were not statistically tested.23 To account for the COVID-19 disruption, the cohort was divided into 2 groups: (1) “not affected” and (2) “affected” for the primary analysis. The not affected group included participants with a 6-month visit where at least 1 of the primary outcomes was measured before March 15, 2020 (ie, not affected by protocol disruptions). Participants who did not have a 6-month visit were also considered not affected if their baseline visit was before September 15, 2019 (ie, 6 months before the COVID-19 shutdown). The remaining participants were considered affected owing to some exposure to COVID-19 adaptations of the study protocol during the active study period.

For each group (ie, affected and not affected), generalized estimating equation models were used to assess the effectiveness of the intervention on the primary outcomes (ie, BMIp95 and submaximal HR). An unstructured covariance structure was used in all models to account for correlation between repeated measures. Within-family correlations were weak and thus were not accounted for in the models (sensitivity analyses not shown). The models for the primary analyses included a “time*intervention” term to test for an intervention effect. Additionally, the model for the BMIp95 outcome was adjusted for race and ethnicity, and the model for the 4-minute HR outcome was adjusted for race and ethnicity, age, and sex. Sensitivity analyses were conducted with all participants regardless of COVID-19 disruptions, and additional models were also tested for subgroup differences by race, age (<12 years vs ≥12 years), and sex. The results of subgroup analyses should be interpreted with caution because the study was not sufficiently powered for these analyses. All analyses were conducted using an intent-to-treat approach. Descriptive statistics were used to describe participant engagement with the intervention. All analyses were performed using SAS version 9.4.

The flow of participant recruitment is shown in Figure 1. The final sample included 255 youths with a mean (SD) age of 10 (3.0) years. Nearly half of the participants were male (47.8%), and the majority were either Hispanic (38.8%) or non-Hispanic Black (38.4%). More than half of the participants had class I obesity (54.5%). Full baseline youth and family-level characteristics are shown in Table 1. Slightly fewer than half (112/255, 44%) were affected by COVID-19 disruptions to study protocols. Youth and family-level characteristics of participants not affected by and affected by COVID-19 disruptions are shown in Supplementary File 2.

FIGURE 1.

Flow of participation in the Hearts & Parks trial. ITT, intention to treat.

FIGURE 1.

Flow of participation in the Hearts & Parks trial. ITT, intention to treat.

Close modal
TABLE 1.

Baseline Participant Demographics

Control (n = 124)Intervention (n = 131)Total (N = 255)
Age, mean (SD), years 10.2 (3.0) 9.8 (3.0) 10.0 (3.0) 
Sex (female), n (%) 67 (54.0) 66 (50.4) 133 (52.2) 
Race and ethnicity, n (%) 
 Non-Hispanic white 21 (16.9) 21 (16.0) 42 (16.5) 
 Non-Hispanic Black 45 (36.3) 53 (40.5) 98 (38.4) 
 Hispanic 53 (42.7) 46 (35.1) 99 (38.8) 
 Othera 5 (4.0) 11 (8.4) 16 (6.3) 
BMI category, n (%) 
 Overweightb 7 (5.6) 13 (9.9) 20 (7.8) 
 Class I obesity 75 (60.5) 64 (48.9) 139 (54.5) 
 Class II obesity 33 (26.6) 42 (32.1) 75 (29.4) 
 Class III obesity 9 (7.3) 12 (9.2) 21 (8.2) 
Annual family income 
 ≤$19 999, n (%) 33 (30.3) 35 (31.5) 68 (30.9) 
 $20 000–$74 999, n (%) 51 (46.8) 45 (40.5) 96 (43.6) 
 ≥$75 000, n (%) 25 (22.9) 31 (27.9) 56 (25.5) 
 Missing, n 15 20 35 
Participating caregiver 
 Mother, n (%) 103 (87.3) 108 (90.8) 211 (89.0) 
 Father, n (%) 12 (10.2) 8 (6.7) 20 (8.4) 
 Other, n (%) 3 (2.5) 3 (2.5) 6 (2.5) 
 Missing, n 12 18 
Participating caregiver education 
 Less than a high school degree, n (%) 25 (21.4) 27 (22.7) 52 (22.0) 
 High school degree or GED, n (%) 25 (21.4) 24 (20.2) 49 (20.8) 
 Some college or technical school, n (%) 36 (30.8) 33 (27.7) 69 (29.2) 
 A 4-y college degree or higher, n (%) 31 (26.5) 35 (29.4) 66 (28.0) 
 Missing, n 12 19 
Participating caregiver employment status 
 Employed for wages, n (%) 64 (55.2) 81 (69.8) 145 (62.5) 
 Out of work, n (%) 10 (8.6) 7 (6.0) 17 (7.3) 
 Homemaker, n (%) 27 (23.3) 19 (16.4) 46 (19.8) 
 Other: student, retired, or self-employed, n (%) 15 (12.9) 9 (7.7) 24 (10.4) 
 Missing, n 15 23 
Food insecure (yes),c n (%) 46 (39.7) 37 (31.9) 83 (35.8) 
Eligible for free and reduced-price lunch (yes),d n (%) 70 (60.3) 72 (59.0) 142 (59.7) 
Language spoken at home 
 English, n (%) 66 (56.4) 77 (65.3) 143 (60.9) 
 Spanish, n (%) 22 (18.8) 21 (17.8) 43 (18.3) 
 Other, n (%) 0 (0.0) 2 (1.7) 2 (0.9) 
 Multiple, n (%) 29 (24.8) 18 (15.3) 47 (20.0) 
 Missing, n 13 20 
Control (n = 124)Intervention (n = 131)Total (N = 255)
Age, mean (SD), years 10.2 (3.0) 9.8 (3.0) 10.0 (3.0) 
Sex (female), n (%) 67 (54.0) 66 (50.4) 133 (52.2) 
Race and ethnicity, n (%) 
 Non-Hispanic white 21 (16.9) 21 (16.0) 42 (16.5) 
 Non-Hispanic Black 45 (36.3) 53 (40.5) 98 (38.4) 
 Hispanic 53 (42.7) 46 (35.1) 99 (38.8) 
 Othera 5 (4.0) 11 (8.4) 16 (6.3) 
BMI category, n (%) 
 Overweightb 7 (5.6) 13 (9.9) 20 (7.8) 
 Class I obesity 75 (60.5) 64 (48.9) 139 (54.5) 
 Class II obesity 33 (26.6) 42 (32.1) 75 (29.4) 
 Class III obesity 9 (7.3) 12 (9.2) 21 (8.2) 
Annual family income 
 ≤$19 999, n (%) 33 (30.3) 35 (31.5) 68 (30.9) 
 $20 000–$74 999, n (%) 51 (46.8) 45 (40.5) 96 (43.6) 
 ≥$75 000, n (%) 25 (22.9) 31 (27.9) 56 (25.5) 
 Missing, n 15 20 35 
Participating caregiver 
 Mother, n (%) 103 (87.3) 108 (90.8) 211 (89.0) 
 Father, n (%) 12 (10.2) 8 (6.7) 20 (8.4) 
 Other, n (%) 3 (2.5) 3 (2.5) 6 (2.5) 
 Missing, n 12 18 
Participating caregiver education 
 Less than a high school degree, n (%) 25 (21.4) 27 (22.7) 52 (22.0) 
 High school degree or GED, n (%) 25 (21.4) 24 (20.2) 49 (20.8) 
 Some college or technical school, n (%) 36 (30.8) 33 (27.7) 69 (29.2) 
 A 4-y college degree or higher, n (%) 31 (26.5) 35 (29.4) 66 (28.0) 
 Missing, n 12 19 
Participating caregiver employment status 
 Employed for wages, n (%) 64 (55.2) 81 (69.8) 145 (62.5) 
 Out of work, n (%) 10 (8.6) 7 (6.0) 17 (7.3) 
 Homemaker, n (%) 27 (23.3) 19 (16.4) 46 (19.8) 
 Other: student, retired, or self-employed, n (%) 15 (12.9) 9 (7.7) 24 (10.4) 
 Missing, n 15 23 
Food insecure (yes),c n (%) 46 (39.7) 37 (31.9) 83 (35.8) 
Eligible for free and reduced-price lunch (yes),d n (%) 70 (60.3) 72 (59.0) 142 (59.7) 
Language spoken at home 
 English, n (%) 66 (56.4) 77 (65.3) 143 (60.9) 
 Spanish, n (%) 22 (18.8) 21 (17.8) 43 (18.3) 
 Other, n (%) 0 (0.0) 2 (1.7) 2 (0.9) 
 Multiple, n (%) 29 (24.8) 18 (15.3) 47 (20.0) 
 Missing, n 13 20 

Abbreviations: BMI, body mass index; GED, general education degree.

a

“Other” included: American Indian or Alaskan Native, Native Hawaiian or Pacific Islander, or multiple races.

b

Some participants had obesity at eligibility screening but lost weight between screening and baseline.

c

A total of 116 control participants and 116 intervention participants reported food security status.

d

A total of 116 control participants and 122 intervention participants reported free and reduced-price lunch status.

The effect of Fit Together on BMIp95 and fitness by COVID-19 affected status is shown in Table 2. For youths not affected by COVID-19 disruptions, intervention youths experienced a significant decrease in BMIp95 (β = −3.05; 95% confidence interval [CI], −5.08 to −1.01) as compared with control. However, we did not observe a significant BMIp95 difference between intervention and control youths affected by COVID-19 disruptions (β = −3.25; 95% CI, −7.98 to 1.48). In a sensitivity analysis comparing the entire cohort, intervention youths had a significant decrease in BMIp95 compared with control youths (β = −3.32; 95% CI, −5.69 to −0.96). Because of the discontinuation of the fitness assessment, comparisons are only available for the not affected group. There was no significant difference between intervention and control youths in 4-minute HR (β = −7.18; 95% CI, −16.12 to 1.76).

TABLE 2.

Effect of the Fit Together Model on BMIp95 and Fitness of Youths Aged 5–17 Years With Obesity by COVID-19 Affected Status

Baseline6 Months
nMean (SD)nMean (SD)Adjusted Estimated
Mean Differences of
Differences (95% CI)a,b,c
P Value
BMIp95 
 Not affected 
  Control 68 116.3 (14.2) 48 116.8 (13.1) −3.05 (−5.08 to −1.01) 0.003 
  Intervention 75 117.1 (14.7) 50 116.1 (15.5) – – 
 Affected 
  Control 56 116.3 (13.8) 40 121.9 (16.8) −3.25 (−7.98 to 1.48) 0.178 
  Intervention 56 119.9 (15.9) 28 123.0 (18.5) – – 
 All participants 
  Control 124 116.3 (14.0) 88 119.1 (15.0) −3.32 (−5.69 to −0.96) 0.006 
  Intervention 131 118.3 (15.3) 78 118.6 (16.9) – – 
Fitness: 4-minute heart rate 
 Not affected 
  Control 35 106.5 (19.4) 38 112.3 (16.1) −7.18 (−16.12 to 1.76) 0.115 
  Intervention 49 102.4 (18.9) 44 102.5 (16.1) – – 
 Affectedd 
  Control 43 107.9 (13.7) – – – 
  Intervention 50 99.1 (18.9) – – – 
 All participantsd 
  Control 78 107.3 (16.4) 38 112.3 (16.1) – – 
  Intervention 99 100.7 (18.9) 44 102.5 (16.1) – – 
Baseline6 Months
nMean (SD)nMean (SD)Adjusted Estimated
Mean Differences of
Differences (95% CI)a,b,c
P Value
BMIp95 
 Not affected 
  Control 68 116.3 (14.2) 48 116.8 (13.1) −3.05 (−5.08 to −1.01) 0.003 
  Intervention 75 117.1 (14.7) 50 116.1 (15.5) – – 
 Affected 
  Control 56 116.3 (13.8) 40 121.9 (16.8) −3.25 (−7.98 to 1.48) 0.178 
  Intervention 56 119.9 (15.9) 28 123.0 (18.5) – – 
 All participants 
  Control 124 116.3 (14.0) 88 119.1 (15.0) −3.32 (−5.69 to −0.96) 0.006 
  Intervention 131 118.3 (15.3) 78 118.6 (16.9) – – 
Fitness: 4-minute heart rate 
 Not affected 
  Control 35 106.5 (19.4) 38 112.3 (16.1) −7.18 (−16.12 to 1.76) 0.115 
  Intervention 49 102.4 (18.9) 44 102.5 (16.1) – – 
 Affectedd 
  Control 43 107.9 (13.7) – – – 
  Intervention 50 99.1 (18.9) – – – 
 All participantsd 
  Control 78 107.3 (16.4) 38 112.3 (16.1) – – 
  Intervention 99 100.7 (18.9) 44 102.5 (16.1) – – 

Abbreviations: BMIp95, body mass index relative to the 95th percentile for age and sex; CI, confidence interval.

a

BMIp95 models were adjusted for race and ethnicity.

b

The fitness model was adjusted for age, sex, and race and ethnicity.

c

Only participants who had measures at both baseline and follow-up contributed to the estimation of the adjusted models.

d

The fitness measure was discontinued during the COVID-19 pandemic because of an inability to assess virtually; thus, no 6-month outcome data are available for participants affected by COVID-19.

In exploratory subgroup analyses, there were no significant between-group differences by age category (ie, <12 years vs ≥12 years) and sex (ie, female vs male) for either BMIp95 or 4-minute HR (Table 3). For youths affected by COVID-19 disruptions, there was a significant difference such that non-Hispanic Black youths had decreased BMIp95 (β = −8.58; 95% CI, −15.42 to −1.72).

TABLE 3.

Effect of the Fit Together Model on BMIp95 and Fitness of Youths With Obesity by Age, Sex, Race and Ethnicity, and COVID-19 Affected Status

Not Affected by COVID-19 Protocol DisruptionAffected by COVID-19 Protocol Disruption
naAdjusted Estimated Mean Differences of Differences (95% CI)P ValuenaAdjusted Estimated Mean Differences of Differences (95% CI)P Value
BMIp95 
 Age 
  <12 years 102 −2.83 (−5.25 to −0.42) 0.022 90 −3.53 (−8.84 to 1.77) 0.192 
  ≥12 years 41 −3.43 (−7.06 to 0.21) 0.065 22 −0.82 (−10.96 to 9.33) 0.874 
  ≥12 years (ref: <12 years) 143 −0.59 (−4.94 to 3.75) 0.789 112 2.72 (−9.00 to 14.43) 0.650 
 Sex 
  Female 79 −1.84 (−4.42 to 0.74) 0.162 54 −6.90 (−12.89 to −0.91) 0.024 
  Male 64 −4.83 (−7.91 to −1.74) 0.002 58 −0.99 (−7.90 to 5.92) 0.779 
  Male (ref: female) 143 −2.99 (−7.00 to 1.03) 0.145 112 5.91 (−3.18 to 15.00) 0.203 
 Race and ethnicity 
  NH white 29 −6.19 (−10.33 to −2.05) 0.003 13 0.33 (−8.45 to 9.10) 0.942 
  NH Black 46 −2.39 (−5.71 to 0.93) 0.159 52 −8.58 (−15.45 to −1.72) 0.014 
  Hispanic 57 −1.77 (−4.85 to 1.31) 0.260 42 2.94 (−4.46 to 10.35) 0.436 
  NH Black (ref: NH white) 75 3.80 (−1.51 to 9.11) 0.161 65 −8.91 (−20.06 to 2.23) 0.117 
  Hispanic (ref: NH white) 86 4.41 (−0.75 to 9.57) 0.094 55 2.61 (−8.87 to 14.10) 0.655 
  Hispanic (ref: NH Black) 103 0.62 (−3.92 to 5.15) 0.790 94 11.53 (1.43 to 21.63) 0.025 
Fitness: 4-minute heart rateb 
 Age 
  <12 years 82 −10.51 (−20.57 to −0.44) 0.041 – – – 
  ≥12 years 36 1.93 (−14.63 to 18.50) 0.819 – – – 
  ≥12 years (ref: <12 years) 118 12.44 (−6.32 to 31.21) 0.194 – – – 
 Sex 
  Female 67 −12.51 (−23.75 to −1.28) 0.029 – – – 
  Male 51 3.25 (−10.05 to 16.55) 0.632 – – – 
  Male (ref: female) 118 15.76 (−1.53 to 33.05) 0.074 – – – 
 Race and ethnicity 
  NH white 25 1.24 (−15.62 to 18.09) 0.886 – – – 
  NH Black 35 −6.24 (−23.24 to 10.76) 0.472 – – – 
  Hispanic 50 −12.23 (−26.71 to 2.24) 0.098 – – – 
  NH Black (ref: NH white) 60 −7.48 (−31.81 to 16.85) 0.547 – – – 
  Hispanic (ref: NH white) 75 −13.47 (−35.98 to 9.04) 0.241 – – – 
  Hispanic (ref: NH Black) 85 −5.99 (−28.05 to 16.07) 0.595 – – – 
Not Affected by COVID-19 Protocol DisruptionAffected by COVID-19 Protocol Disruption
naAdjusted Estimated Mean Differences of Differences (95% CI)P ValuenaAdjusted Estimated Mean Differences of Differences (95% CI)P Value
BMIp95 
 Age 
  <12 years 102 −2.83 (−5.25 to −0.42) 0.022 90 −3.53 (−8.84 to 1.77) 0.192 
  ≥12 years 41 −3.43 (−7.06 to 0.21) 0.065 22 −0.82 (−10.96 to 9.33) 0.874 
  ≥12 years (ref: <12 years) 143 −0.59 (−4.94 to 3.75) 0.789 112 2.72 (−9.00 to 14.43) 0.650 
 Sex 
  Female 79 −1.84 (−4.42 to 0.74) 0.162 54 −6.90 (−12.89 to −0.91) 0.024 
  Male 64 −4.83 (−7.91 to −1.74) 0.002 58 −0.99 (−7.90 to 5.92) 0.779 
  Male (ref: female) 143 −2.99 (−7.00 to 1.03) 0.145 112 5.91 (−3.18 to 15.00) 0.203 
 Race and ethnicity 
  NH white 29 −6.19 (−10.33 to −2.05) 0.003 13 0.33 (−8.45 to 9.10) 0.942 
  NH Black 46 −2.39 (−5.71 to 0.93) 0.159 52 −8.58 (−15.45 to −1.72) 0.014 
  Hispanic 57 −1.77 (−4.85 to 1.31) 0.260 42 2.94 (−4.46 to 10.35) 0.436 
  NH Black (ref: NH white) 75 3.80 (−1.51 to 9.11) 0.161 65 −8.91 (−20.06 to 2.23) 0.117 
  Hispanic (ref: NH white) 86 4.41 (−0.75 to 9.57) 0.094 55 2.61 (−8.87 to 14.10) 0.655 
  Hispanic (ref: NH Black) 103 0.62 (−3.92 to 5.15) 0.790 94 11.53 (1.43 to 21.63) 0.025 
Fitness: 4-minute heart rateb 
 Age 
  <12 years 82 −10.51 (−20.57 to −0.44) 0.041 – – – 
  ≥12 years 36 1.93 (−14.63 to 18.50) 0.819 – – – 
  ≥12 years (ref: <12 years) 118 12.44 (−6.32 to 31.21) 0.194 – – – 
 Sex 
  Female 67 −12.51 (−23.75 to −1.28) 0.029 – – – 
  Male 51 3.25 (−10.05 to 16.55) 0.632 – – – 
  Male (ref: female) 118 15.76 (−1.53 to 33.05) 0.074 – – – 
 Race and ethnicity 
  NH white 25 1.24 (−15.62 to 18.09) 0.886 – – – 
  NH Black 35 −6.24 (−23.24 to 10.76) 0.472 – – – 
  Hispanic 50 −12.23 (−26.71 to 2.24) 0.098 – – – 
  NH Black (ref: NH white) 60 −7.48 (−31.81 to 16.85) 0.547 – – – 
  Hispanic (ref: NH white) 75 −13.47 (−35.98 to 9.04) 0.241 – – – 
  Hispanic (ref: NH Black) 85 −5.99 (−28.05 to 16.07) 0.595 – – – 

Abbreviations: BMIp95, body mass index relative to the 95th percentile for age and sex; CI, confidence interval; NH, non-Hispanic; ref, reference.

a

The number of participants used in the analysis.

b

The fitness measure was discontinued during the COVID-19 pandemic because of an inability to assess virtually; thus, no 6-month outcome data are available for participants affected by COVID-19.

There were also several significant within-group differences for youths not affected by COVID-19 disruptions. Younger (β = −2.83; 95% CI, −5.25 to −0.42), male (β = −4.83; 95% CI, −7.91 to −1.74), and non-Hispanic white (β = −6.19; 95% CI, −10.33 to −2.05) intervention youths had significant decreases in BMIp95 compared with control youths. Female intervention youths had a significant improvement in 4-minute HR compared with control female youths (β = −12.51; 95% CI, −23.75 to −1.28).

Mean (SD) engagement hours for intervention participants was 8.9 (15.5) hours (range, 0–80 hours). Participants not affected by COVID-19 disruptions had a higher mean number of engagement hours compared with those affected by COVID-19 disruptions (11.5 ± 17.1 vs 5.4 ± 12.5 hours, respectively). In total, 10% (13/131) of participants received a high-intensity dose (ie, ≥26 hours) and nearly one-third of participants (30%, 39/131) received a moderate-intensity dose (ie, 5–25 hours). The hours of engagement for intervention participants are shown in Figure 2. No adverse events related to the study protocol were reported.

FIGURE 2.

Hours of engagement for intervention participants.

FIGURE 2.

Hours of engagement for intervention participants.

Close modal

This study found that Fit Together, a pediatric obesity treatment intervention delivering IHBLT through a clinic-community partnership, led to improved weight outcomes in youths with obesity; however, no effect on cardiorespiratory fitness was detected. In addition, Fit Together was able to reach a diverse group of youths (>80% from marginalized racial and ethnic groups). The findings highlight the effectiveness of this clinic-community model and will inform future pediatric obesity medicine practice and research.

Recommendations from the AAP and US Preventive Services Task Force support the use of IHBLT to treat pediatric obesity.8 In meta-analyses examining the effectiveness of behavioral treatment programs, trials with at least 26 contact hours demonstrated an increased intervention effect on BMI (mean difference, − 1.4) compared with trials with less than 26 hours (mean difference, − 0.3).24 However, few interventions have been able to deliver the recommended intensity, even in clinical programs designed to treat pediatric obesity.25 Findings from this study show that the Fit Together model provides an avenue to implement IHBLT and achieve BMI reductions. Of note, few youths in this study met the at least 26-hour IHBLT recommendation, suggesting a need to further investigate the dose-response relationship and heterogeneity in treatment response within IHBLT interventions.

Findings in this study are similar to those of previous work demonstrating the effectiveness of engaging community organizations in the treatment of youth obesity. For instance, the Mind, Exercise, Nutrition, Do it (MEND) program engaged youths with obesity and their families in community-based physical activity and education, which yielded significant decreases in youth BMI in several trials.26,27 MEND was adapted to the Healthy Weight and Your Child program that has been widely disseminated through YMCA organizations.28,29 This suggests that the Fit Together model could potentially be scaled and sustained within community settings.

Although findings were in the expected direction, there was no significant effect of Fit Together on youth cardiorespiratory fitness. The most likely reason is insufficient power because of discontinuation of the fitness assessment during the COVID-19 pandemic. Additionally, it is likely that youths who experienced reductions in BMIp95 also made dietary changes that contributed to reductions in weight, which may explain the reduction in BMI without changes in fitness. In subgroup analyses, there was a significant difference in fitness among girls not affected by COVID-19 disruptions, which may be attributed to improved physical activity in a population that traditionally has lower physical activity; however, this should be interpreted with caution.30,31 Given the cardiometabolic benefits of improved fitness, even in the absence of weight loss, it is necessary for future obesity treatment interventions to continue to explore the effect on youth fitness.32 

Prior pediatric obesity treatment studies have often failed to reach priority populations, including marginalized racial and ethnic groups and low-wealth populations that experience a greater obesity burden.33,34 In the current study, we were able to engage youths and families from predominantly priority populations, which may be attributed to the access provided by partnering with Durham P&R. However, subgroup analyses only revealed a significant effect among non-Hispanic white youths, suggesting a need to further understand the effect of this intervention among marginalized racial and ethnic groups.

This study had several limitations. First, COVID-19 safety measures affected both data collection and intervention delivery. Assessment of youth fitness was discontinued because of difficulties administering the protocol virtually, leading to missing data for 1 of the co-primary outcomes. The intervention was shifted to an online format, and although both synchronous and asynchronous options were available, engagement with the online program was low and likely did not match the intensity of in-person session. Additionally, although the Healthy Lifestyles–P&R partnership was feasible in this study, it may not be generalizable to primary care. Healthy Lifestyles is a specialized tertiary care obesity treatment clinic affiliated with an academic medical institution. However, other studies suggest that this model can be implemented in primary care settings.14 

A partnership among health care and P&R providers demonstrated effectiveness in improving BMI within a program designed to deliver IHBLT for youths with obesity. Future research should seek to further understand the heterogeneity of treatment response with respect to BMI reduction and other key health outcomes as well as avenues for dissemination and the sustainability of this approach.

Drs Armstrong and Skinner conceptualized and designed the study, coordinated and supervised data collection, oversaw analyses, interpreted results, and reviewed and revised the manuscript. Dr Neshteruk oversaw analyses, interpreted results, drafted the initial manuscript, and reviewed and revised the manuscript. Drs. Li, Kraus, Shah, Story, Zucker, and Perrin and Mr Jones conceptualized and designed the study and critically reviewed the manuscript for important intellectual content. Ms Zizzi and Dr Windom were responsible for data collection and reviewed the manuscript for important intellectual content. Dr Wagner interpreted results, drafted the initial manuscript, and reviewed and revised the manuscript. Dr Hong led the analyses and reviewed and revised the manuscript. Ms Truong and Mr Burrows carried out the analyses and reviewed and revised the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

CONFLICT OF INTEREST DISCLOSURES: The authors have no conflicts of interest to disclose.

FUNDING: This study was supported by funding from the American Heart Association (17SFRN33670990, 17SFRN33671003, 17SFRN33700117, and 17SFRN33700155), the National Institutes of Health (R24-DK110492 and 5R01-HL127009), and the Duke Endowment. Additionally, Dr Neshteruk was supported by the National Institutes of Health (5K12HL138030). None of the funding sources had a role in the design, execution, analysis, or interpretation of the data of this study.

ClinicalTrials.gov (NCT03339440)

COMPANION PAPER: A companion to this article can be found online at www.pediatrics.org/cgi/doi/10.1542/peds.2024-069092.

We would like to thank all the youths and families who participated in this project. Additionally, we would like to acknowledge Emily D’Agostino, Jessilyn Dunn, Janna Howard, Neha Pagidipati, Eric Peterson, and Charlene Wong for their contributions and insights.

AAP

American Academy of Pediatrics

BMI

body mass index

BMIp95

body mass index relative to the 95th percentile for age and sex

BPA

best-practice alert

CI

confidence interval

HR

heart rate

IHBLT

intensive health behavior and lifestyle treatment

MEND

Mind, Exercise, Nutrition, Do it

P&R

parks and recreation

PCP

primary care physician

1
Jastreboff
AM
,
Kotz
CM
,
Kahan
S
,
Kelly
AS
,
Heymsfield
SB
.
Obesity as a disease: the Obesity Society 2018 position statement
.
Obesity (Silver Spring)
.
2019
;
27
(
1
):
7
9
. PubMed doi: 10.1002/oby.22378
2
Stierman
B
,
Afful
J
,
Carroll
MD
, et al
.
National Health and Nutrition Examination Survey 2017–March 2020 prepandemic data files development of files and prevalence estimates for selected health outcomes
.
Natl Health Stat Report.
2021
;(
158
). PubMed doi: 10.15620/cdc:106273
3
Stierman
B
,
Ogden
CL
,
Yanovski
JA
,
Martin
CB
,
Sarafrazi
N
,
Hales
CM
.
Changes in adiposity among children and adolescents in the United States, 1999-2006 to 2011-2018
.
Am J Clin Nutr
.
2021
;
114
(
4
):
1495
1504
. PubMed doi: 10.1093/ajcn/nqab237
4
Lange
SJ
,
Kompaniyets
L
,
Freedman
DS
, et al;
DNP3
.
Longitudinal trends in body mass index before and during the COVID-19 pandemic among persons aged 2–19 years - United States, 2018–2020
.
MMWR Morb Mortal Wkly Rep
.
2021
;
70
(
37
):
1278
1283
. PubMed doi: 10.15585/mmwr.mm7037a3
5
Umer
A
,
Kelley
GA
,
Cottrell
LE
,
Giacobbi
P
Jr
,
Innes
KE
,
Lilly
CL
.
Childhood obesity and adult cardiovascular disease risk factors: a systematic review with meta-analysis
.
BMC Public Health
.
2017
;
17
(
1
):
683
. PubMed doi: 10.1186/s12889-017-4691-z
6
Moradi
M
,
Mozaffari
H
,
Askari
M
,
Azadbakht
L
.
Association between overweight/obesity with depression, anxiety, low self-esteem, and body dissatisfaction in children and adolescents: a systematic review and meta-analysis of observational studies
.
Crit Rev Food Sci Nutr
.
2022
;
62
(
2
):
555
570
. PubMed doi: 10.1080/10408398.2020.1823813
7
Skinner
AC
,
Perrin
EM
,
Moss
LA
,
Skelton
JA
.
Cardiometabolic risks and severity of obesity in children and young adults
.
N Engl J Med
.
2015
;
373
(
14
):
1307
1317
. PubMed doi: 10.1056/NEJMoa1502821
8
Hampl
SE
,
Hassink
SG
,
Skinner
AC
, et al
.
Clinical practice guideline for the evaluation and treatment of children and adolescents with obesity
.
Pediatrics
.
2023
;
151
(
2
):
e2022060640
. PubMed doi: 10.1542/peds.2022-060640
9
Nicholson
WK
,
Silverstein
M
,
Wong
JB
, et al;
US Preventive Services Task Force
.
Interventions for high body mass index in children and adolescents: US Preventive Services Task Force recommendation statement
.
JAMA
.
2024
;
332
(
3
):
226
232
. PubMed doi: 10.1001/jama.2024.11146
10
Haemer
M
,
Cluett
S
,
Hassink
SG
, et al
.
Building capacity for childhood obesity prevention and treatment in the medical community: call to action
.
Pediatrics
.
2011
;
128
Suppl 2(
Suppl 2
):
S71
7
. PubMed doi: 10.1542/peds.2011-0480G
11
Dietz
WH
,
Solomon
LS
,
Pronk
N
, et al
.
An integrated framework for the prevention and treatment of obesity and its related chronic diseases
.
Health Aff (Millwood)
.
2015
;
34
(
9
):
1456
1463
. PubMed doi: 10.1377/hlthaff.2015.0371
12
Gordon-Larsen
P
,
Nelson
MC
,
Page
P
,
Popkin
BM
.
Inequality in the built environment underlies key health disparities in physical activity and obesity
.
Pediatrics
.
2006
;
117
(
2
):
417
424
. PubMed doi: 10.1542/peds.2005-0058
13
Roemmich
JN
,
Epstein
LH
,
Raja
S
,
Yin
L
,
Robinson
J
,
Winiewicz
D
.
Association of access to parks and recreational facilities with the physical activity of young children
.
Prev Med
.
2006
;
43
(
6
):
437
441
. PubMed doi: 10.1016/j.ypmed.2006.07.007
14
Neshteruk
CD
,
Skinner
AC
,
Counts
J
, et al
.
Translating knowledge into action for child obesity treatment in partnership with Parks and Recreation: study protocol for a hybrid type II trial
.
Implement Sci
.
2023
;
18
(
1
):
6
. PubMed doi: 10.1186/s13012-023-01264-5
15
Hoffman
J
,
Frerichs
L
,
Story
M
, et al
.
An integrated clinic-community partnership for child obesity treatment: a randomized pilot trial
.
Pediatrics
.
2018
;
141
(
1
):
e20171444
. PubMed doi: 10.1542/peds.2017-1444
16
Andrews
M
,
Sawyer
C
,
Frerichs
L
, et al
.
Feasibility of a clinic-community partnership to treat childhood obesity
.
Clin Pediatr (Phila)
.
2018
;
57
(
7
):
783
791
. PubMed doi: 10.1177/0009922817734359
17
Armstrong
SC
,
Windom
M
,
Bihlmeyer
NA
, et al
.
Rationale and design of “Hearts & Parks”: study protocol for a pragmatic randomized clinical trial of an integrated clinic-community intervention to treat pediatric obesity
.
BMC Pediatr
.
2020
;
20
(
1
):
308
. PubMed doi: 10.1186/s12887-020-02190-x
18
Schulz
KF
,
Altman
DG
,
Moher
D
;
CONSORT Group
.
CONSORT 2010 statement: updated guidelines for reporting parallel group randomised trials
.
BMJ
.
2010
;
340
(
mar23 1
):
c332
. PubMed doi: 10.1136/bmj.c332
19
Ogden
CL
,
Kuczmarski
RJ
,
Flegal
KM
, et al
.
Centers for Disease Control and Prevention 2000 growth charts for the United States: improvements to the 1977 National Center for Health Statistics version
.
Pediatrics
.
2002
;
109
(
1
):
45
60
. PubMed doi: 10.1542/peds.109.1.45
20
Hales
CM
,
Freedman
DS
,
Akinbami
L
,
Wei
R
,
Ogden
CL
.
Evaluation of alternative body mass index (BMI) metrics to monitor weight status in children and adolescents with extremely high BMI using CDC BMI-for-age growth charts
.
Vital Health Stat 1
.
2022
;(
197
):
1
42
.
21
Freedman
DS
,
Berenson
GS
.
Tracking of BMI z scores for severe obesity
.
Pediatrics
.
2017
;
140
(
3
):
e20171072
. PubMed doi: 10.1542/peds.2017-1072
22
Jankowski
M
,
Niedzielska
A
,
Brzezinski
M
,
Drabik
J
.
Cardiorespiratory fitness in children: a simple screening test for population studies
.
Pediatr Cardiol
.
2015
;
36
(
1
):
27
32
. PubMed doi: 10.1007/s00246-014-0960-0
23
de Boer
MR
,
Waterlander
WE
,
Kuijper
LDJ
,
Steenhuis
IHM
,
Twisk
JWR
.
Testing for baseline differences in randomized controlled trials: an unhealthy research behavior that is hard to eradicate
.
Int J Behav Nutr Phys Act
.
2015
;
12
(
1
):
4
. PubMed doi: 10.1186/s12966-015-0162-z
24
O’Connor
E
,
Evans
C
,
Henninger
M
,
Redmond
N
,
Senger
C
,
Thomas
R
. Interventions for High Body Mass Index in Children and Adolescents: An Evidence Update for the US Preventive Services Task Force. Agency for Healthcare Research and Quality;
2024
.
25
Kirk
S
,
Armstrong
S
,
King
E
, et al
.
Establishment of the Pediatric Obesity Weight Evaluation Registry: a national research collaborative for identifying the optimal assessment and treatment of pediatric obesity
.
Child Obes
.
2017
;
13
(
1
):
9
17
. PubMed doi: 10.1089/chi.2016.0060
26
Butte
NF
,
Hoelscher
DM
,
Barlow
SE
, et al
.
Efficacy of a community- versus primary care-centered program for childhood obesity: TX CORD RCT
.
Obesity (Silver Spring)
.
2017
;
25
(
9
):
1584
1593
. PubMed doi: 10.1002/oby.21929
27
Sacher
PM
,
Kolotourou
M
,
Chadwick
PM
, et al
.
Randomized controlled trial of the MEND program: a family-based community intervention for childhood obesity
.
Obesity (Silver Spring).
2010
;
18
(S1)(
suppl 1
):
S62
S68
. PubMed doi: 10.1038/oby.2009.433
28
Liebhart
JL
,
Hassink
SG
,
Hohman
KH
,
Reilly
E
,
Lindros
J
,
Womack
SJ
;
The Healthy Weight and Your Child Collaborative
.
A Pilot Project to Improve Obesity-related Primary Care and Cultivate Clinic-Community Linkages for Pediatric Weight Management
.
American Academy of Pediatrics
;
2019
.
29
Fiechtner
L
,
Perkins
M
,
Biggs
V
, et al
.
Comparative effectiveness of clinical and community-based approaches to healthy weight
.
Pediatrics
.
2021
;
148
(
4
):
e2021050405
. PubMed doi: 10.1542/peds.2021-050405
30
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
188
. PubMed doi: 10.1249/mss.0b013e31815a51b3
31
Armstrong
S
,
Wong
CA
,
Perrin
E
,
Page
S
,
Sibley
L
,
Skinner
A
.
Association of physical activity with income, race/ethnicity, and sex among adolescents and young adults in the United States: findings from the National Health and Nutrition Examination Survey, 2007-2016
.
JAMA Pediatr
.
2018
;
172
(
8
):
732
740
. PubMed doi: 10.1001/jamapediatrics.2018.1273
32
Nyström
CD
,
Henriksson
P
,
Martínez-Vizcaíno
V
, et al
.
Does cardiorespiratory fitness attenuate the adverse effects of severe/morbid obesity on cardiometabolic risk and insulin resistance in children? A pooled analysis
.
Diabetes Care
.
2017
;
40
(
11
):
1580
1587
. PubMed doi: 10.2337/dc17-1334
33
Caprio
S
,
Daniels
SR
,
Drewnowski
A
, et al
.
Influence of race, ethnicity, and culture on childhood obesity: implications for prevention and treatment: a consensus statement of Shaping America’s Health and the Obesity Society
.
Diabetes Care
.
2008
;
31
(
11
):
2211
2221
. PubMed doi: 10.2337/dc08-9024
34
Skelton
JA
,
Irby
MB
,
Beech
BM
,
Rhodes
SD
.
Attrition and family participation in obesity treatment programs: clinicians’ perceptions
.
Acad Pediatr
.
2012
;
12
(
5
):
420
428
. PubMed doi: 10.1016/j.acap.2012.05.001

Supplementary data