Design:
This study was part of a multi-site, longitudinal cohort study of Canadian children with newly diagnosed JIA (BBOP; Biologically-based Outcome Predictors in JIA). BBOP aimed to identify the inter-relationships of biologic, environmental, and lifestyle factors as predictors of childhood arthritis outcomes. 15 We analyzed BBOP data to characterize the trajectory of PA levels in JIA.
Participants had a clinical assessment and completed all study questionnaires at baseline, 12 and 24 months. Where developmentally appropriate, children completed their own self-report questionnaires. Otherwise, they were completed by a parent.
In compliance with the Helsinki Declaration, the study was approved by the Biomedical Research Ethics Board, University of Saskatchewan: #07-86 and by the research ethics boards at each of the other 10 participating sites. Parents and, as appropriate, children provided written informed consent; children provided assent, as applicable.
BBOP Recruitment:
Children with a new diagnosis of JIA were enrolled at 11 Canadian pediatric rheumatology centers between December 2007 and December 2012. Enrollment criteria included: (i) consenting participants who met International League of Associations for Rheumatology (ILAR) JIA classification criteria 16 and (ii) were diagnosed within six months of symptom onset. As physical activity assessment tools were applicable only for participants aged 6 to 16 years, younger BBOP participants were not included. The cohort comprised participants from each of the seven ILAR JIA subtypes. The aim was to recruit sufficient numbers for each category rather than to strive to achieve a typical JIA subset distribution. To achieve this, only participants with polyarthritis or systemic JIA (the least prevalent subtypes) were eligible during the first six months of enrollment after which any JIA subtype was eligible. Participants who were unable to communicate in English or French were excluded from the study.
Reference population
The reference group for PA measures for this BBOP study analysis were the 154 Canadian children (82 females and 72 males) who participated in the University of Saskatchewan's Pediatric Bone Mineral Accrual Study (1991-97, entry age 8 to 15 years). 17
Measurements:
Physical Activity
Participants completed the Physical Activity Questionnaire (PAQ) at all three study visits. 1819 The PAQ is a self-administered, 7-day physical activity recall questionnaire that assesses participation in different physical activities, as well as activity during physical education class, lunch break, recess, after school, in the evenings and on weekends. The Physical Activity Questionnaire for Children (PAQ-C) has been validated in children 5 to 12 years of age20 and was completed by school age children > 6 years and ≤ 12 years of age in this study. 18 The Physical Activity Questionnaire for Adolescents (PAQ-A) is recommended for adolescents and was completed by children ≥ 13 years of age.21 The PAQ-C/A was completed at visits during the school year and through the summer break. Both measures have been validated in healthy children and in children with chronic disease. 2122
The PAQ-A and PAQ-C comprise 8 and 9 PA questionnaire items, respectively. Each item is scored between 1 (low PA) and 5 (high PA); a mean score of all items constitutes the overall PAQ score. Questionnaire items 9 (PAQ-A) and 10 (PAQ-C), which ask participants whether anything prevented them from doing normal physical activities, are only used to gauge whether the responses represent the individual’s typical PA levels and are not included in the calculation of the PAQ-C/A score.
Clinical Assessment
The attending pediatric rheumatologist recorded the total number of joints with active arthritis (0-71 joints) and provided a physician global assessment (PGA) of participant disease activity measured on a 0–10 visual analog scale (VAS) where 0 = no activity and 10 = maximum activity. Participants were asked to complete a parent/patient global assessment of well-being, measured on a 10-cm VAS where 0 = very well and 10 = very poor. If participants required bloodwork as part of their routine clinical care, the erythrocyte sedimentation rate (ESR) was recorded and normalized to a 0-10 scale as outlined by Consolaro et al. 23 The PGA of disease activity, parent/patient global assessment of well-being, active joint count, and ESR were used to calculate the Juvenile Arthritis Disease Activity Score (JADAS). 24 The JADAS-71 was calculated as the simple sum of the scores of its 4 components, which yields a global score of 0–101.
Participants completed the Child Health Assessment Questionnaire (CHAQ) 25 as a parent-reported measure of their functional impairment. As part of the CHAQ, they were asked to rate their pain in the past week on a 10-cm visual analogue scale (VAS) from 0 (no pain) to 10 (very severe pain). 2627
Socioeconomic information
Participants completed a questionnaire at study enrollment that ascertained place of residence (as determined by the Canadian postal code forward sortation area), population of the community in which the participant resided, parental marital status, parental education, parental occupation, and household income.
Psychosocial Stress
Three self-administered questionnaires were completed at all study visits to assess psychosocial stressors.
The Juvenile Arthritis Quality of Life Questionnaire (JAQQ) measures physical and psychosocial functioning, specifically rating difficulty with functions due to arthritis or its treatment. 28 The JAQQ is a self-administered questionnaire that is applicable to all age groups and JIA subtypes. It consists of 74 items distributed across 4 dimensions. The psychosocial function dimension was used for this study and includes 22 items that are scored in terms of severity from 1 (none of the time) to 7 (all of the time). The mean of the 5 highest scoring items was calculated to give the JAQQ-psychosocial score. Mean scores range from 1 to 7, with 7 indicating the worst function.
The Children’s Hassles Scale measures the frequency and impact of daily hassles, defined as irritating and/or distressing demands that to some degree characterize everyday interactions with the environment. 29 It consists of 25 items scored from 0 (this did not happen) to 3 (yes this happened, and it made you feel very sad). Mean scores range from 0 to 3, with 3 indicating the highest impact of daily hassles.
The Stressful Life Events Checklist (SLEC) measures the occurrence of stressful life events in the prior year. 30 Two versions of the checklist were used: (i) a 40-item checklist for children < 12 years with total scores ranging from 0 (no stressors occurred in the past year) to 40 (all stressors occurred in the past year) and (ii) a 47-item checklist for adolescents ≥ 12 with total scores ranging from 0-47.
Examples of the JAQQ, Hassles Scale and SLEC questionnaire items are shown in Supplementary Table 1.
Statistical Analysis:
Statistical analyses were performed with SPSS Statistics Professional, version 23, R, version 3.2.2., and MATLAB R2019a. The data had 21% missing values, which were imputed by Amelia II, R package which imputes missing values employing an expectation maximization algorithm. 31,32 Age- and sex-specific z-scores for the PAQ-C/A were generated for 8 to 15-year-old participants using normative data from the Saskatchewan Pediatric Bone Mineral Accrual Study. 33Mann–Whitney U test was used to examine the differences between girls’ and boys’ PA scores. The Kruskal-Wallis test was used to assess the difference between PA scores of the seven JIA subtypes. The associations between PAQ-C/A and clinical disease activity (JADAS-71), function (CHAQ), SES and measures of psychosocial stress (JAQQ psychosocial domain, SLEC, and Hassles scale) were assessed at baseline using Spearman correlations. A linear mixed model for panel data was used to estimate the association between PA over time with age, sex, and disease activity measures including number of active joints, ESR, c-reactive protein (CRP), and JADAS-71 (fixed effect). The measures of psychosocial stress and function were considered as variables with random effects. The best model was selected using a theoretical likelihood ratio test.