Introduction

Osteoarthritis (OA) is a chronic condition of the synovial joint, which includes the progressive degeneration of cartilage and the excess growth of bone, often leading to pain and functional impairment [1]. It is one of the leading causes of global disability, with adult prevalence rates reported between 8.5–22% for symptomatic radiographic knee OA [24] and 3.4–8.9% for symptomatic radiographic hip OA [3, 5, 6]. To determine risk factors for this disease, it is necessary to analyse previously collected data from longitudinal population cohorts. While there are some well-established risk factors for hip and knee OA, the relationship between physical activity (PA) and OA is inconsistent. This, in part, may be due to the heterogeneous definition of PA used in cohorts and the lack of differentiation between weight-bearing and non-weight-bearing activity.

PA is defined as any bodily movement that results in energy expenditure and is categorised by domains, including occupation, leisure time, daily living, and active travel. Assessment of PA in cohort studies of OA is usually captured by self-report questionnaires, typically including measures of frequency, intensity, duration, and type [7]. Much of the previous research on PA has been completed within the area of cardiovascular disease and obesity and many of the questionnaires and assessments were developed with such health outcomes in mind. Although observational OA-related cohorts have collected information on PA, the parameters measured differ between studies and the number of domains varies, making comparison and interpretation of results difficult. Equally as important is the lack of an available method to assess the degree of joint loading for different physical activities, which is a known risk factor for OA.

A recent systematic review of observational studies, subject to the above limitations, confirmed the association of injury, obesity, and occupational activity with knee and hip OA; however, increased volume and intensity of PA were found to be a risk factor for OA in four studies and protective in another [8]. Due to the heterogeneity in the definition of PA, meta-analysis was not possible for the risk of PA on OA.

To address the difficulties in comparing heterogeneous aggregate data within epidemiological research, an increasingly popular alternative to meta-analysis is individual participant data (IPD) meta-analysis, where the raw individual level data are used for statistical synthesis [9]. This method allows for using the combined power of multiple international cohort studies to address more complex research questions, such as the association between PA and OA [1018]. A key element of IPD meta-analysis is the harmonization of the variables, which requires a standard measure of both PA and OA from all the included cohorts. Unfortunately, population-based cohort studies have collected PA data using a variety of validated questionnaires and individual questions, which need to be harmonised to be included in any analysis.

To enable this, a number of issues need to be addressed: PA questionnaires and questions vary between cohorts; not all domains are available in all cohorts; some activities are attributes to varying domains (e.g., walking and cycling is included with sport and recreation in a number of cohorts, yet travel in others); duration, intensity, and frequency of PA are not addressed consistently between cohorts; and PA has not previously been assessed by degree of weight-bearing for use in studies of lower limb OA.

A consensus study, including an expert meeting and Delphi approach, was developed to address these issues. The aims of this study were to:

  1. 1.

    Determine the usability of a common metric (Metabolic equivalent of task [MET]) as a key method for harmonising PA assessments/questions between cohorts, and agree upon specific assumptions required to generate METs for each cohort (objectives 1 and 2).

  2. 2.

    Assess the available domains of PA and establish the appropriate assumptions needed to harmonise information between cohorts (objectives 3 and 4).

  3. 3.

    Evaluate the potential to use of a lower limb OA-specific PA measure within the cohorts, taking weight-bearing into consideration (objective 5).

  4. 4.

    Evaluate the use of national PA guidelines to determine the effect of meeting such recommendations on the association with lower limb OA (objective 6).

Methods

The process consisted of the following steps

PA expert committee

A multidisciplinary, geographically diverse expert committee was selected on the basis of individual expertise in PA; exercise medicine and OA and each were invited to participate in developing a PA variable for use in normal population-based cohort studies. The expert committee (n = 9 of the listed co-authors) met by video conference link in December 2014.

Expert consensus meeting

The steering group (consisting of authors LG, KL, and NA) conducted a systematic evaluation of international cohorts containing PA and OA data. Four key issues were identified by the steering group whilst establishing the methods to harmonise PA data between cohorts which had not been adequately addressed in the previous literature. Background topic information, the previous research where applicable and specific factors relating to individual cohorts were provided to experts before the meeting to assist informed decision-making. Each issue was presented and facilitated discussion undertaken until agreement was reached.

Delphi

If a question was raised within the expert meeting and no agreement could be made, due to either a requirement for further investigation or exploration of cohort data, then these items were addressed within a follow-up Delphi technique to obtain consensus. The Delphi technique, which is a structured process of anonymous iteration [19], consisted of an online questionnaire, which was sent to each expert, following the meeting. In accordance with the previous OARSI Delphi exercises to define OA diagnostic criteria [20], inclusion for measures within each round was based on having ≥60% of the votes. Evidenced-based information was provided to inform each item, and where applicable more in depth data from individual cohorts were provided. Unanimous decision-making was made within the first round.

The aims of the consensus study were addressed via six methodological objectives:

  1. 1.

    Determine the suitability of using METs as a key method for harmonising PA exposure variables between cohorts.

  2. 2.

    Determine methods for treating missing components of MET-min/week calculation.

  3. 3.

    Assess the domains of physical activity: how to treat missing domains.

  4. 4.

    Assess the domains of physical activity: the use of occupation as a PA domain in studies with OA as an outcome.

  5. 5.

    Evaluate the use of an OA-specific PA measure taking weight-bearing vs. non-weight-bearing activity into consideration.

  6. 6.

    Establish if thresholds based on national PA guidelines should be used to investigate the association of PA with OA.

Results

Determine the suitability of using METs as a key method for harmonising PA exposure variables between cohorts

Background information provided to experts

The steering group proposed using the MET to harmonise PA data across population cohorts, based on the results of the literature search and availability of measures within cohorts. PA data can be converted to MET-min/week using the 2011 Compendium of Physical Activities [21] if METs were not already reported. MET is a physiological measure which expresses the intensity and energy expenditure of an activity, and is defined as the ratio of the rate of energy expended during an activity to the rate of energy expended at rest (e.g., 1 MET is the rate of energy expenditure while at rest) [22]. Standard intensity thresholds have been established for MET values of activities with <3 as light, 3–5.9 as moderate, and ≥6 as vigorous [23]. MET-min/week is calculated by multiplying the MET value of an activity (from a standard Compendium of Physical Activities values [21] by the number of minutes per week an activity is done [24]). METs have been used for recommending activity in adults for chronic disease prevention and health promotion [23].

Points for expert discussion

Due to the variability between cohorts in the way PA exposure was reported and the wording of PA questions, the benefits, limitations, and feasibility of converting the questions from each cohort into MET-min/week as the standard unit were proposed as the first discussion point to be addressed within the expert consensus meeting.

Determine methods for treating missing components of MET-min/week calculation

Background information provided to experts

To convert PA exposure data to MET-min/week, the duration, frequency, and intensity of an activity are required. In certain cohorts within the current study, some of these parameters were not collected, and a number of assumptions needed to be made. The previous studies have adopted a standard bout time of 30 min when duration was missing from the physical activity question, but the validity of this assumption is unclear [25].

Points for expert discussion

Several common PA questionnaires do not include intensity, frequency, or duration of the activity, preventing a direct calculation of MET-min/week. A possible solution proposed by the steering group was to use a standard intensity MET value when calculating MET-min/week if the terminology included “low”, “moderate”, or “high”, and to use a standard of 30 min when duration was missing. These issues raised a requirement for additional investigation via a further Delphi study.

Further investigation recommended

For the Delphi exercise, median duration data for individual activities were prepared by the steering group from the cohorts where these data were available (Supplementary Appendix 1). These data were to be presented for agreement within the Delphi to ascertain agreement on the generalisability of the median durations for each activity. Due to some variation between standard activity times between countries, a proposal was made to use the median durations from the Hertfordshire Cohort Study ; a UK study of 3000 men and women born during the period 1931–1939 [15], for other UK cohorts with missing duration data, and similarly to use the durations from Johnston County Osteoarthritis Study: a US community-based, longitudinal study of approximately 3200 rural Caucasian and African American residents aged 45 and older [11], for US cohorts with missing data. The same method was proposed where frequency was identified as missing in one UK cohort (Supplementary Appendix 2).

Assess the domains of physical activity: how to treat missing domains

Background information provided to experts

In previous research, household activities have been included as an essential domain of PA to determine the minimum amount of PA associated with significantly lower risks of all-cause mortality [26, 27]. Household activities have also been specifically included within PA health guidelines [28]. The importance of household activities is also anticipated in determining the association of PA with lower limb OA, particularly due to the known risks of increased cumulative loads seen in such activities and knee OA [29]. Household activities are largely modifiable and are, therefore, an important consideration for any public health message involving physical activity.

Points for expert discussion

A significant number of the population-based cohorts did not ask questions regarding household activities and/or gardening, both of which have been identified as important contributors to daily activity loads. Experts were asked to consider the impact of excluding these domains from the primary physical activity variable, and if there were alternative actions to consider.

Assess the domains of physical activity: the use of occupation as a PA domain in studies with OA as an outcome

Background information provided to experts

Occupation, in particular manual labour occupations, is a well-established risk factor for OA [30, 31], and is, therefore, an important consideration in an analysis using OA as the primary outcome. The occupation section within the Compendium of Physical Activities [21] was designed for counting different tasks within one occupation over 1 day and is of limited use for providing an “average” MET value for a specific occupation (i.e., MET values of individual tasks which make up a complete occupation are given and are expected to be done for an entire hour). Occupation is the least modifiable domain of PA and approximately 50% of the consortium cohorts had no or limited occupation data available. It would, therefore, be difficult to calculate a valid exposure in MET-min/week from the available data. The potential for occupation activities to be overestimated is greater than for any other domain, particularly because the occupation MET value within the Compendium of Physical Activities assumes that a specific MET activity is completed for entire hour. To multiply this by hours, worked in 1 week would potentially overestimate weekly METs. The majority of consortium cohorts, which did have occupation data available, did not define tasks within occupation in anything less than 1 day.

Points for expert discussion

Experts were asked to discuss the limitations of excluding occupation from the overall physical activity variable specific to studies using osteoarthritis as an outcome. The preferred treatment of an osteoarthritis variable (e.g., MET-min/week and categories) was also explored.

Further investigation recommended

A method was required for categorisation into levels, defining occupations into manual and non-manual tasks. Due to the variety of occupation-related tasks available in each cohort, the occupation level reported in the Physical Activity Scale for the Elderly (PASE) questionnaire was used to facilitate the categorisation of occupation-related tasks. The PASE questionnaire is designed to assess the duration, frequency, exertion level, and amount of physical activity undertaken over a 7 day period [32]. The steering group matched the PASE occupation levels to the corresponding four levels of occupation, which were suggested in the expert consensus meeting. The selection of occupation-related tasks with each of the four levels of occupation (sedentary, light, light manual and heavy manual) was prepared to be presented within the follow-up Delphi exercise (Supplementary Appendix 3).

Evaluate the use of an OA-specific PA measure taking weight-bearing vs. non-weight-bearing activity into consideration

Background information provided to experts

OA is at least in part a mechanically driven disease [33]. The association between PA and risk of OA may, therefore, be dependent on joint loading and type of PA [34].

Points for expert discussion

There was a need to identify a method that could account for the weight-bearing component of PA. The fifth proposed discussion point for the meeting was, therefore, the potential to use previously established bone and joint loading questionnaires to quantify loading [35, 36].

Further investigation recommended

Informed by these decisions, each activity listed within the study cohorts was placed into corresponding loading categories of low, moderate, and high joint loading based on the degree of impact and loading (Supplementary Appendix 4). These were prepared by the steering group for presentation within a follow-up Delphi exercise.

Establish if thresholds based on national PA guidelines should be used to investigate the association of PA with OA

Background information provided to experts

It is important to ensure that a translatable public health message is provided; however, the relationship of meeting or not meeting PA guidelines and the risk of OA are unknown. The use of a threshold based on national guidelines such as those from the American College of Sports Medicine and the American Heart Association [37], the US Physical Activity Guidelines Advisory Committee (US PAGAC) report [23], and the U.K. Department of Health [38] (150 min/week of moderate-equivalent activity) was, therefore, proposed to investigate the effect of meeting current guidelines on the risk of OA.

An additional threshold was suggested to investigate the risk of inactivity on risk of OA; however, there is no global consensus on a defined threshold for inactivity. According to the US PAGAC report [23], <10 min/week of moderate activity is defined as ‘inactive’ and this is in comparison to the UK definition of ‘inactivity’ of <30 min/week moderate activity.

There is growing evidence that increasing steps per day provides health benefits [3942]. Although originally based on minimal evidence, the 10,000 steps a day guideline now provides a translatable and applicable PA recommendation [43].

Points for expert discussion

The final points proposed for discussion within the expert meeting were the use of current national guidelines as a threshold for PA in the investigation of the PA and OA. Also discussed was the use of either a previous arbitrary threshold for inactivity or a new data driven threshold, and the use of steps per day as an equivalent for the translatable outcome instead of 150 min of weekly moderate-equivalent activity.

Expert consensus meeting results

Within the expert meeting, each aim was discussed until consensus reached. Where new questions were raised, these were addressed within a follow-up Delphi exercise. There was unanimous agreement for every objective within the consensus meeting and follow-up Delphi. In summary, key agreements were made based upon: the use of MET-min/week as a method for harmonising PA variables between cohorts; defining methods for treating missing components of MET-min/week calculation, in particular the use of a value produced from comparable activities within a representative cohort; the exclusion of the domain of ‘occupation’ from total MET-min/week; and the need for a specific measure of ‘joint loading’ of an activity in studies of bone diseases, such as osteoarthritis.

Details of the decisions made within the consensus meeting and follow-up Delphi are shown in Tables 1 and 2.

Table 1 Decisions made within the expert consensus meeting
Table 2 Modified Delphi results (questions were based on evidence and supplementary information provided in Appendices 1–4)

Discussion

Our research describes a method to classify and harmonise PA data for epidemiology research in population-based OA-related cohorts, based on international consensus. The recommendations will allow the IPD meta-analysis of PA and incident lower limb OA to be undertaken using the most comprehensive PA data possible. In addition, it will be useful not only for future epidemiology research into OA that uses physical activity and IPD studies, but also to guide researchers planning to collect PA data for current or future epidemiological research.

Recommendations for combining data between cohorts, based on expert opinion in this study, are to use MET-min/week as a standardised measure of PA. In the instance of a missing parameter, such as duration, frequency, or intensity, the methods of assigning values using data derived from representative cohorts or the Compendium of Physical Activities have been agreed. The likely effects of occupation should be accounted for by categorisation of occupation type and stratification. The weight-bearing aspect within PA should be taken into consideration when using OA as an outcome.

Although national PA guidelines are an essential and translatable source of health prescription, they were not designed with OA in mind. Guidelines are required to address health and, therefore, recommendations relative to OA, an increasing global burden, are required to add to other disease areas. There are a number of domains of PA to consider when assessing the target for meeting national recommendations; these include leisure time, household/gardening, active travel, and occupational activities. These domains are particularly pertinent when considering the effect of PA on lower limb OA due to the weight-bearing nature of many activities within each and the difficulty is that there is currently no index to combine the physiological measure of METs and joint load.

An increasingly popular alternative to meeting 150 min/week of moderate-equivalent activity is the daily accumulation of 10,000 steps. Experts agreed that this could be a useful method of assessing PA against OA in the future. Recommendations from the expert consensus study suggest that a guideline based on the personalised optimum number of steps per day to reduce the risk of OA would be a valuable public health message.

A limitation of this study was that decisions had to be based on already data that were already available, because due to the requirement to use existing population-based cohorts to investigate the association between PA and OA. There are also known limitations for using METs as a measure of PA, particularly when making comparisons among a number of studies or populations. Studies have previously measured METs based on a varying numbers of activities from the total available from the domains of PA [44]. As a MET is the total volume of a given activity, which combines frequency, intensity, and duration, the effect of duration or intensity alone cannot be deciphered. IPD meta-analysis allows for the use of original raw data, so that all aspects of PA can be included and where no available in certain circumstances imputation methods can be applied. Likewise, the access to original data allows for data driven thresholds to derived and provides potential for observing individual parameters of PA, be it intensity or duration.

A further limitation to this study was the use of PA data from more than one decade ago to calculate median duration for those cohorts missing this data. Although PA levels are likely to have changed, since these data were collected the committee felt, this was still a more appropriate representation of duration than using a standard 30 min, which is likely to over or underestimate activity levels.

Our study provides the first expert consensus on the limitations of and the methods for harmonising PA data in population-based OA cohort studies. The application of these recommendations in future individual patient meta-analysis on PA and OA will provide a homogeneous way to assess PA in cohorts from around the world. It will also allow for quantifying the volume of PA and examine the shape of the dose–response curve for PA and OA as well as the ability to apply new thresholds for future national PA guidelines. These findings will also be useful for any study investigating PA and other long-term health outcomes in existing cohort data. The recommendations arising from this consensus study for the collection of PA data in normal population-based cohort studies are: the need for all parameters of a given activity (duration, frequency, and type/intensity) within a specified timeframe; PA measured in all domains of daily life (sport/leisure, household/gardening, active travel, and occupation); an occupation measure which can be used to calculate accurate MET-min/week value in addition to a manual labour (occupational tasks) in terms of OA risk; and a measure of joint loading for each reported activity.