Skip to main content

Knee pain trajectories over 18 months in non-Hispanic Black and non-Hispanic White adults with or at risk for knee osteoarthritis

Abstract

Background

Pain is the hallmark symptom of knee osteoarthritis (OA), and varies widely across individuals. Previous research has demonstrated both fluctuating and stable pain trajectories in knee OA using various time periods. Changes in pain assessed quarterly (i.e. 3-month intervals) in knee OA are relatively unknown. The current study aimed to investigate temporal variations in pain over a one and a half year period (18 months) based on quarterly characteristic pain assessments, and to examine differences in pain patterns by sociodemographic and baseline pain characteristics.

Methods

The sample included a prospective cohort of 188 participants (mean age 58 years; 63% female; 52% non-Hispanic Black) with or at risk for knee OA from an ongoing multisite investigation of ethnic/race group differences. Knee pain intensity was self-reported at baseline and quarterly over an18-month period. Baseline pain assessment also included frequency, duration, and total number of pain sites. Group-based trajectory modeling was used to identify distinct pain trajectories. Multinomial logistic regression was used to examine associations between sociodemographic characteristics, risk factors, and pain trajectory groups.

Results

Pain trajectories were relatively stable among a sample of adults with knee pain. Four distinct pain trajectories emerged in the overall sample, with the largest proportion of participants (35.1%) classified in the moderate-high pain group. There were significant relationships between age, education, income, ethnicity/race and trajectory group; with younger, less educated, lower income, and non-Hispanic Black participants had a greater representation in the highest pain trajectory group.

Conclusions

Pain remained stable across a one and a half-year period in adults with or at risk for knee osteoarthritis, based on quarterly assessments. Certain sociodemographic variables (e.g. ethnicity/race, education, income, age) may contribute to an increased risk of experiencing greater pain.

Peer Review reports

Background

Pain is the hallmark symptom of knee osteoarthritis (OA) and the primary reason people with knee OA seek medical care [1, 2]. Knee OA is among the leading causes of disability in older adults, with approximately 33 million adults in the U.S. diagnosed [3,4,5]. There is substantial variation in pain among individuals with or at risk for knee OA [6]. Multiple factors, in addition to radiographic disease progression [7, 8], influence the severity of knee OA-related pain. This is evidenced by prior research which identified distinct knee OA pain phenotypes based on neurophysiological and psychological indices [9,10,11]. There is also growing evidence for distinct pain trajectories within individuals with knee OA, which may help to inform clinical care. To date, the majority of studies have focused on annual assessments of pain to model pain progression, leaving critical questions regarding shorter-term changes, particularly quarterly changes in pain over one and a half years.

Recent scientific advances demonstrating broad heterogeneity in knee OA symptom and disease progression have challenged the historical perspective of OA being a slowly progressive disease characterized by increasing pain and articular joint degeneration [12,13,14,15,16,17,18,19]. A growing body of research examining pain and physical function trajectories in knee OA demonstrate pain progression is heterogeneous [16, 19,20,21,22,23,24]. The majority of studies have relied on annual assessments with follow-up ranging from 5 to 10 years, and have reported both progressive and non-progressive pain trajectories. For example, in a six-year longitudinal investigation of pain progression in persons with symptomatic knee OA, researchers found five pain trajectories ranging from “mild, non-progressive” to “severe, non-improving” [19]. Another study investigating pain progression annually over 5 years identified six unique pain trajectories ranging from “constant mild pain” to “constant severe pain”, with evidence for pain regression and progression [25]. One study identified five unique and stable pain trajectories based on annual assessments over a five-year period [16]. Despite these conflicting findings, it is evident that individuals with or at risk for knee OA vary in pain progression over time. It is likely that assessments taken more frequently might lead to more precise estimates of pain changes over time.

Evidence from a recent study assessing pain daily over a one-month time period, reported patients demonstrated relative stability in pain fluctuations [26]. However, the study was limited by the large percentage of participant drop out (approximately 50%) and by the use of descriptive statistics only. In another study investigating short-term fluctuations in OA-related pain over 26 weeks, researchers using latent class growth analysis, identified three distinct pain trajectories (i.e., “high, persistent pain”, “moderate, persistent pain”, and “low pain, improving”) [27]. However, this study was also limited by a high dropout rate (20%), and relatively short data collection period. While appealing, daily assessments are not feasible to capture changes in pain over a longer time span, which is critical for understanding the nature of symptom development. For example, researchers have noted that significant symptom progression can be seen in individuals with knee OA over a 1–2 year period [28,29,30,31,32,33]. Given that OA-related pain can be experienced as persistent or intermittent, with ongoing fluctuations in intensity, longitudinal studies using annual or bi-annual assessments may not fully capture clinical pain variability [18], and daily assessments place undue burden on patients. Therefore, identification of pain trajectories in knee OA based upon quarterly assessments may prove optimal by balancing the benefit to burden ratio [34].

Furthermore, despite evidence for associations between sociodemographic factors and pain trajectories [18, 34], little is known about the relationship between ethnicity/race and pain progression in knee OA. Previous research has demonstrated considerable health disparities in knee OA, with African Americans/non-Hispanic Blacks (NHBs) reporting greater overall levels of clinical pain and physical disability [21, 35,36,37], and experimental pain sensitivity [10, 37, 38]. Over a nine-year period of annual assessments, researchers found that African Americans (AAs) reported significantly higher pain at baseline and each measurement time-point compared to Whites (WHs) [21]. Furthermore, AAs reported more pain variation over time than WHs, yet this variability failed to reach a statistically significant difference. Overall, there was little year-to-year change in mean reported pain, except for the first year in which AAs showed a clinically significant decline in pain from baseline to 12 months [21]. Further examination of pain fluctuations using a more granular approach are warranted to better understand how pain progression influences disparities in knee OA. It is possible that greater baseline clinical pain predisposes individuals to a more severe pain trajectory, which if identified early through more frequent assessments, would improve understanding and help inform interventions to improve health outcomes.

The objectives of this study were to: (1) identify distinct pain trajectories of individuals with or at-risk for, symptomatic knee OA based on quarterly assessments over an 18-month period; 2) examine differences in sociodemographic factors, including ethnicity/race, across identified pain trajectories; and, 3) examine baseline pain characteristics among identified pain trajectory groups in NHB and NHW individuals, with or at-risk-for knee OA. We hypothesized that there would be heterogeneity (i.e. different patterns) of pain trajectories in the sample, and that these patterns would be associated with sociodemographic variables and clinical pain characteristics.

Methods

Study design

The current study is a secondary analysis of data from the Understanding Pain and Limitations in Osteoarthritic Disease-2 (UPLOAD-2), an ongoing, prospective multisite study that aims to elucidate the biopsychosocial mechanisms underlying ethnic/race group differences in knee pain between non-Hispanic Black (NHB) and non-Hispanic White (NHW) adults with or at risk for knee osteoarthritis (OA). The data set allows for the investigation of pain trajectories across 18 months using seven quarterly (3-month) pain assessments. The study was conducted at the University of Florida (UF) and the University of Alabama (UAB). Full study methods have been previously published [39, 40]. The study was approved by Institutional Review Boards at UF and UAB and was completed following the institutionally approved study protocol. Findings are presented in accordance with STROBE reporting guidelines.

Participants

Participants from the UPLOAD2 parent study were included in the current analysis if they had completed at least one quarterly health assessment and other measures of interest. This resulted in a sample of 188 individuals who were: (1) between 45 and 85 years of age, (2) self-identified as non-Hispanic Black (NHB) or non-Hispanic White (NHW), (3) having or being at risk for unilateral or bilateral symptomatic knee OA based upon American College of Rheumatology criteria [41], (4) able to complete multiple sessions, and (5) able to read and speak English [42]. All participants screened positive on a brief questionnaire that has previously shown 87% sensitivity and 92% specificity for radiographically confirmed symptomatic knee OA [43]. As the parent study was designed to evaluate progression of OA-related symptoms, a cohort of individuals were enrolled with a broad range of knee pain characteristics. Individuals were defined as at-risk of developing knee OA if they responded positively to the knee OA screening questions [43] but did not demonstrate radiographic joint degeneration (i.e., Kellgen- Lawrence grade ≥ 2) [44]. The broad inclusion criteria provides the opportunity to better understand factors associated with knee pain progression rather than OA pathophysiology itself. Exclusion criteria included: (1) prosthetic knee replacement or other clinically significant surgery to the arthritic knee; (2) heart disease, congestive heart failure, or history of acute myocardial infarction; (3) peripheral neuropathy; (4) systematic rheumatic disorders (e.g., rheumatoid arthritis, systemic lupus erythematosus, and fibromyalgia; (5) chronic daily opioid use; (6) neurological diseases such as Parkinson’s, multiple sclerosis, stroke with loss of sensory or motor function, or uncontrolled seizures; (7) greater pain in body sites other than knee; (8) hospitalization within preceding year for psychiatric illness; or (9) pregnant (positive human chorionic gonadotropin urine test) or breastfeeding. Participants were recruited from the community through posted fliers, radio and print media, orthopedic clinic recruitment and word-of-mouth referral, between August 2015 and May 2017. Participants provided written informed consent at the baseline visit, prior to study procedures. Following the baseline session, participants were informed that they would receive prompts via email and/or phone, to complete quarterly health assessments. Participants were compensated up to $300 for their involvement.

Procedures

Participant eligibility was determined through an initial telephone screening. One to two weeks following screening, enrolled participants completed an in-person health assessment session (HAS) which included informed consent, pain history questionnaires and a bilateral knee joint evaluation by the study’s rheumatologists or nurse practitioner. Participants were classified as either having, or being at risk for knee OA. After the examination, participants were asked to complete online questionnaires assessing self-reported baseline knee pain symptoms. Approximately 3 months following baseline data collection, participants completed a quarterly self-reported knee pain assessment via their preferred mode of responding (online, paper, telephone), which continued to be distributed each quarter over the 18-month period, for a total of seven assessment points.

Measures

Sociodemographic and health characteristics

Participants self-reported age, sex, ethnicity/race, highest education level, household income, during the telephone screening. BMI and Kellgren-Lawrence scores were assessed at baseline.

Baseline pain characteristics

Clinical pain was assessed at baseline (prior to quarterly pain assessments), across several domains important to understanding chronic pain trajectories [45, 46]. These included: (1) knee pain frequency (“On average, how many days per week do you experience pain in your knee?”); (2) knee pain intensity (worst, average, current); (3) knee pain duration (“For how long have you been experiencing pain in your knee?”); and, (4) total number of painful body sites. Knee pain frequency of more than 5 days per week was classified as “experiencing knee pain on most days”. Knee pain intensity was assessed using three items from the Graded Chronic Pain Scale (GCPS) [47], specific to knee pain (i.e., worst and average pain over the past 6 months, and current pain), rated on a 0 (“no pain”) to 10 (“pain as bad as could be”) numeric rating scale (NRS). These ratings were averaged and multiplied by 10 to generate a characteristic knee pain intensity score [47]. The NRS has been shown to be a reliable (ICC = 0.95) and valid (r = 0.94, visual analog scale) measure of OA pain [48, 49]. Knee pain duration was categorized as: (1) less than 6 months; (2) 6 months to 1 year; (3) 1 to 3 years; (4) 3 to 5 years; or, (5) more than 5 years. The total number of painful body sites was calculated as the sum of self-reported painful body sites reported. Participants could indicate if pain was experienced on the left, right, or both sides across 14 body sites. Items were coded as follows: hands (right or left = 1, both = 2); arms (right or left = 1, both = 2); neck (right or left = 1, both = 2); shoulders (right or left = 1, both = 2); head/face/jaw (right or left = 1, both = 2); chest (right or left = 1, both = 2); stomach (right or left = 1, both = 2); pelvis (right or left = 1, both = 2); upper back (right or left = 1, both = 2); lower back (right or left = 1, both = 2); knees (right or left = 1, both = 2); legs (right or left = 1, both = 2); feet (right or left = 1, both = 2), and/or 1 ‘other’ body region (free response; right or left = 1, both = 2, and summed to produce a total score with a range from 0 to 28. These four variables (pain frequency, pain intensity, pain duration, and pain sites) were summarized using descriptive statistics and compared across identified pain trajectory subgroups [45, 46].

The Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) [50], was administered at baseline to assess knee OA-related symptoms in the past 48 h. The three subscales of the WOMAC (i.e., pain, stiffness, and physical functioning) were summed to create a total symptom burden score ranging from 0 to 96 (i.e., WOMAC Global score), which was used for analysis. Higher scores indicate greater overall symptom burden. The WOMAC is a well-validated measure of pain and function in lower extremity OA [50, 51].

Quarterly pain assessment

Participants completed quarterly self-reported knee pain assessments every 3 months using the GCPS to assess pain intensity (“knee pain at its worst in the past week”; “knee pain on average in the past week”; and “knee pain right now”), each rated on a 0 (“no pain”) to 10 (“pain as bad as you can imagine”) NRS. NRS ratings were averaged to generate a knee pain intensity score ranging from 0 to 10 [47], which were used in to identify pain trajectories in the following analysis.

Statistical analysis

Analyses were performed in JMP Pro 14 and SAS 9.3 (SAS Institute Inc., Cary, NC). Continuous measures were summarized with means and standard deviations and categorical measures were summarized with percentages. To identify subgroups of individuals that have similar progressions (i.e. trajectories) of pain, group-based trajectory modeling (GBTM), using maximum likelihood estimation, was implemented with PROC TRAJ in SAS® software (SAS Institute, Cary, NC) [52,53,54]. GBTM was used to examine heterogeneity in clinical pain patterns; GBTM does not assume that all individuals in the population will follow a similar functional form of development (i.e., it does not assume one trajectory shape will fit all) [55]. Group assignment using GBTM takes a data-driven approach, with each individual clustered into the trajectory group to which they had the highest posterior probability of membership, rather than assigning individuals to groups based on traits or other subjective criteria [36]. It is important to note that these trajectory groups represent latent strata, in that individuals are following approximately the same development pattern. However, individuals do not “belong” to a group, but are assigned based on probability of membership; any given individual could have a developmental pattern that does not exactly follow the group-level trajectory.

GBTM was used to determine both the number of distinct trajectory groups and the shape of each trajectory (i.e. order of polynomial). Model-fitting followed a multi-stage iterative process. First, a one-group model was run with the highest order of polynomial (quartic), then regression terms for each order (zero-order [intercept], linear, quadratic, cubic, and quartic) were evaluated for statistical significance (p < 0.05). Then, additional groups were added (2-,3-,4-,5- groups) until best fitting model was determined. Bayesian information criteria (BIC) was used to identify the most parsimonious, best fitting model, (i.e. the model that has the best fit using the fewest number of trajectories) [56]. To compare model fits, Bayes factor was calculated, which is approximately 2 times the difference in BIC between 2 models (2 x [BIC more complex model – BIC simpler model]) [52]. A Bayes factor > 2 suggests positive evidence to support a meaningful change in BIC for a more complex model, with a Bayes factor ≥ 10 providing very strong evidence [52]. Additionally, to ensure stability of the model, too small group sizes of a given trajectory group are not recommended [54]. Because GBTM utilizes full information maximum likelihood estimation, participants with missing values are still included in the models. Missingness, at each time point, ranged from 29 to 55%, with higher rates at later time points. This pattern of missing data was similar between NHBs and NHWs. A majority of the sample (78.3%) had at least five time points with complete data.

The differences in age, sex, ethnicity/race and baseline pain characteristics among trajectory groups were examined using Chi-square analysis for categorical variables and analysis of variance (ANOVA) for continuous variables. Pairwise post-hoc multiple comparison between groups was estimated using Tukey’s HSD test. Welch’s correction for ANOVA was used if equal variance assumptions were not met. A multivariable nominal logistic regression was also conducted, including sociodemographic and risk factors, to determine independent predictors of trajectory group; effects summarized as odds ratio (OR) with 95% confidence intervals (95%CIs). Note, categories for some variables were combined for logistic regression analyses. P < 0.05 was considered statistically significant. All methods were carried out in accordance with relevant guidelines and regulations [57].

Results

Sample characteristics

The average age of participants was 58 (SD ± 7.8) years, with a majority of participants being female (63%). The sample consisted of approximately even representation of NHBs and NHWs. Nearly half of participants reported having knee pain most days and reported experiencing pain for > 5 years. NHB and NHW participants significantly differed on age (p = 0.005), highest education level (p = 0.03), income (p = 0.004), GCPS score (p < 0.001), and total WOMAC score (p < 0.001). See Table 1 for full sample characteristics.

Table 1 Descriptive Measures (n = 188), for Full Sample and Stratified by Ethnicity/race

Group-based trajectory modeling

Table 2 reports model fitting results for the GBTM analysis. In the overall sample, the best fitting model included 4 trajectory groups, all linear order. Trajectory groups are labeled according to characteristic knee pain intensity patterns (low, moderate-low, moderate-high, and high) (Fig. 1). The largest proportion of participants (35.1%) were in the moderate-high pain group, followed by 31.4% of participants in the moderate-low pain group. The low pain group had the smallest proportion of participants (14.1%), with the remaining 19.3% of participants in the high pain group. An overall pattern of worsening or improvement in reported pain across the 2-year study period was not observed in any of the four trajectory groups, with linear coefficients for each group the best-fitting model nor statistically significantly different from zero (0) (Table 2). Average posterior probabilities for trajectory groups were 80% or higher for all groups, with highest posterior probabilities high pain group: low (88%), moderate-low (87%), moderate-high (80%), and high (95%).

Table 2 Model Fitting for Group Based Trajectory Analysis
Fig. 1
figure 1

Pain trajectory groups in the overall sample over 7 quarterly (3-month intervals) assessments over an 18-month period with 95% confidence intervals. Note. Percentages indicate the proportion of the sample included in each trajectory group

Participant sociodemographic characteristics and trajectory groups

There were differences between trajectory groups for age (F(3,184) = 8.3, p < 0.001), ethnicity/race (χ 2 = 47.5, df = 3, p < 0.001), highest education level (χ 2 = 57.3, df = 3, p < 0.001), household income (χ 2 = 54.8, df = 3, p < 0.001), BMI (F(3,184) = 8.9, p < 0.001), and Kellgren-Lawrence (KL) score (χ 2 = 27.8, df = 3, p = 0.006) (Table 3). Participants in the high pain group (mean age = 54 years) were, on average, a decade younger than participants in the low pain group (mean age = 64 years). Non-Hispanic Black (NHB) participants were not represented in the low pain group and comprised three-quarters of the high pain group. Those with highest level of education of high school or less comprised the majority of the high pain (71.0%) and moderate-high pain (64.2%) groups. Similarly, those with lowest household income (<$20,000) comprise the vast majority of the high pain group (80%). BMI was significantly higher in high pain group compared to all other groups. Finally, the majority of the high pain group include participants with KL scores of 3 or higher, though, nearly 30% of the high pain group also included those with KL scores of 0 or 1. Interestingly, while nearly 60% of the low pain group included those with KL scores of 0 or 1, nearly one-third of the low pain group was comprised of individuals with KL scores of 3 or higher.

Table 3 Sample Characteristics by Pain Trajectory Group

In multivariable nominal logistic regression analysis (Table 4) with sociodemographic and risk factors, age (p = 0.031), ethnicity/race (p = 0.012), education (p = 0.003), household income (p < 0.001), and KL scores (p = 0.001) remained as independent predictors of pain trajectories; odds ratios used combined low and moderate-low pain group as reference. As age increased, participants had decreased likelihood of being in high pain group (OR = 0.90, 95%CI:0.82, 0.99). Compared to NHWs, NHBs were more likely to be in high pain (OR = 5.28, 95%CI:1.5, 18.66) and the moderate-high pain group (OR = 2.57, 95%CI:1.15, 5.71). For household income, those with incomes less than $20,000 were more likely to be in high pain group (OR = 8.50, 95%CI:1.62, 44.49) and moderate-high pain group (OR = 3.36, 95%CI:1.07, 10.51). BMI (p = 0.143) and gender (p = 0.337) were not statistically significant predictors in multivariable analysis.

Table 4 Results from Multivariable Nominal Logistic Regression

Baseline pain characteristics and trajectory groups

There were notable differences among pain trajectory groups in baseline pain characteristics. Specifically, knee pain frequency (χ 2 = 16.6, df = 3, p = 0.001), knee pain intensity (F(3,184) = 37.5, p < 0.001), and total number of pain sites (F(3,184) = 3.9, p = 0.011) differed significantly across pain trajectory groups. Three-quarters of participants in the high pain trajectory group reported experiencing knee pain on most days, compared to less than half of participants in each other pain group (Table 3). Participants in high pain groups also reported the greatest number of pain sites and had the highest knee pain intensity scores. Participants in the low pain trajectory group reported the fewest number of pain sites and had the lowest knee pain intensity scores (Table 3). There were similar trajectory group differences in total WOMAC scores (Table 3), with participants in the higher pain trajectories having higher baseline total WOMAC scores (F(3,182) = 33.8, p < 0.001).

Discussion

The overall intention of the current study was to investigate quarterly variations in pain over an 18-month period and to examine differences in pain patterns by sociodemographic, health, and baseline pain characteristics. Four distinct pain trajectories (i.e., low, moderate-low, moderate-high, and high) emerged within the overall sample. The majority of our participants fell within the moderate-low and moderate-high pain trajectory subgroups, based on quarterly pain assessments. The low pain group had the lowest representation across the sample. On average, identified trajectories remained stable within individuals over time, which is consistent with previous large studies using annual assessments [16, 23, 24, 58, 59].

Identified pain trajectory groups differed across sociodemographic and health variables, including age, ethnicity/race, education, income, BMI, and KL grade. Participants in the low pain trajectory group were older compared to the high pain group. This age finding is consistent with other pain trajectory studies conducted in the US [16], but is in contrast to studies conducted outside of the US, which show significant associations between increasing age and pain [59]. It is possible that older adults with more severe pain were excluded from the current study unintentionally due to the exclusion of individuals with comorbidities, resulting in a relatively healthy sample of older adults who were not experiencing high levels of pain. Another possible explanation is that younger participants are more likely to work and engage in activities which may result in greater self-reported pain. Importantly, Non-Hispanic Blacks (NHBs) were not represented in the low pain group and comprised the majority of the high pain group. Additionally, lower education, lower income, and BMI were also associated with higher pain trajectories.

An important consideration for this sample is that NHBs were, on average, younger and reported lower education and income compared to NHW participants. Despite this potential confounding, in multivariable analysis, age, ethnicity/race, education, and income remained as independent predictors of pain trajectories. In addition to sociodemographic factors, individuals with more significant radiographic joint disease (indicated by higher KL scores) were also more likely to be in higher pain trajectory groups. This finding is consistent with the study by Wesseling and colleagues (2015) indicating a higher KL score placed individuals at a greater risk for a moderate pain trajectory [60], but differs from the study by Collins et al. (2014) which demonstrated no differences in pain trajectories based on KL scores [16].

Pain trajectory groups also differed on baseline clinical pain characteristics such that the majority of individuals in the high pain trajectory group reported experiencing knee pain on most days, greater knee pain intensity at baseline, and more bodily pain sites compared to the other pain trajectory groups. This supports prior research which found comorbid musculoskeletal (i.e. hip) pain [60], and baseline pain severity [25], placed individuals with knee OA at greater risk of developing a moderate pain trajectory [60]. Consistent with prior research NHB participants reported higher levels of baseline clinical pain across trajectory groups compared to NHWs [21, 35].

The current findings regarding greater pain in NHBs with knee OA are supportive of prior research indicating that NHBs/African Americans have increased risk of experiencing greater OA-related pain over time [16, 21]. However, as noted in other research, ethnicity/race is a construct associated with environmental, social, and cultural experiences that impact a number of factors influencing health including multiple chronic pain conditions, access to health care and treatment, and management strategies [37, 61]. However, the current analyses are somewhat limited in the number of variables available to fully examine the complex array of factors contributing to environmental and sociodemographic influences on pain outcomes [62]. While the mechanisms underlying ethnic/race group differences have yet to be fully elucidated, research indicates central sensitization (e.g. hyperexcitability of the central nervous system), may place some individuals with knee OA at a greater risk for increased pain severity [10, 38, 63], and further investigation of factors contributing to observed ethnic/race differences in clinical and experimental pain are underway [10, 64,65,66].

To our knowledge this is among the first studies evaluating pain trajectories in OA using quarterly self-reports. Prior research of pain trajectories has relied on a few large samples, with annual assessments, that may not be fully representative of NHB and NHW adults with, or at-risk-for, knee OA. Investigating changes in knee pain quarterly across an 18-month time span provides a more nuanced appreciation for trajectory patterns of knee pain symptoms that may differ from those assessed annually, and can inform the development of more accurate prognosis pathways for precision care. One study utilizing quarterly pain assessments over 2 years in hip OA classified patients into five pain trajectories and showed symptom progression over time [67]. This contrasts with the current study’s findings and may be partly explained by the differences in symptom progression between hip and knee OA [68]. The lack of observed pain progression within identified trajectories could be due in part to adaptive behaviors, coping skills, and pharmacological pain management [69, 70]. Another explanation is that knee OA may have a sustained period of symptomatic stability followed by symptom escalation [16, 18, 21]. Overall, this has significant clinical relevance in terms of prevention, rehabilitation and symptom management, warranting further investigation.

Implications and future directions

The heterogeneity of pain identified in the current sample, and in prior studies, highlights the inter-individual differences in symptom progression over time. Understanding pain trajectories associated with knee OA can inform clinical care and best treatment practices. For example, individuals experiencing frequent, high intensity knee pain in conjunction with multisite pain may be at a higher risk of poor health outcomes, surgical complications, and may require multimodal pain management approaches [71]. Future efforts aimed at developing algorithms useful for identifying individuals at a higher risk for a greater symptom trajectory could inform the development of targeted interventions for individuals with knee OA pain. Future research examining associations between pain processing (e.g., peripheral and central sensitization), psychological factors and pain trajectories, will further our understanding of the mechanistic contributors to pain, and will help to identify other key indicators to inform clinical decision making. This could include studies into the relationship between pain self-management and temporal changes in symptoms to inform personalized medicine approaches which allow for prognosis-based treatment planning.

Strengths and limitations

The results of this study should be interpreted in light of its strengths and limitations. This was a cohort study with equal representation of NHB and NHW participants with, or at-risk-for, knee OA which allowed for important comparisons. We used group-based multi-trajectory modelling (GBMT), which allows for the possibility of distinct subgroups that can be identified within a population [55, 72]. Our findings add novel information regarding knee pain intensity trajectories across quarterly assessments over an 18-month period, thus reducing potential bias from symptom fluctuations captured at only one time point in a given year [73]. Furthermore, use of a validated measure, the Graded Chronic Pain Scale (GCPS) to assess pain quarterly provides a more robust characterization of symptoms targeted over the past week [47], rather than measures with a shorter symptom assessment period (e.g., WOMAC-48 h).

Several limitations must also be considered, including the proportion of missing data. Also, longitudinal evaluation of group means may obscure within-person variability, and misclassification is possible. Another limitation is our smaller sample (< 200), inadequate sample sizes can bias model solutions and can make it difficult identify smaller subgroups [74]. As quarterly assessments were completed online, we cannot be sure what influence participant environment had on pain assessments. However, participants were provided with the same instructions at each assessment point to reduce potential response bias. Additionally, the use of a single pain dimension (i.e., intensity), may not provide full information regarding symptom progression and health outcomes, despite being a well-validated instrument. Also, with the exception for exclusion if an individual was taking opioids on a daily basis, we did not control for pain medication in this study. However, previous research using the UPLOAD2 data set found that less than half of participants reported using pain medication, which did not differ significantly by ethnicity/race. And finally, our exclusion criteria may limit generalizability of our findings. Despite these limitations, this study improves understanding regarding pain trajectories in knee OA, providing the next level of evidence regarding symptom progression in knee OA by replicating findings from similar, larger studies using annual assessments in sample with more detailed, frequent assessments.

Conclusions

Our findings provide further evidence for stable, yet heterogeneous, pain trajectories in adults with, or at-risk-for, knee OA-related pain with greater temporal resolution. Four stable knee pain trajectories were identified in this sample, and differed significantly in respect to several sociodemographic factors including age, ethnicity/race, education, and income. NHBs who were on average, younger, less educated and reported lower income compared to NHWs were more likely to report high pain levels at the baseline and were categorized in the stable high pain trajectory over 18 months. Taken together, our findings further confirm the stability of pain progression in knee OA and improves our understanding of pain progression and highlights the importance of considering sociodemographic characteristics in understanding pain heterogeneity.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ANOVA:

Analysis of variance

BIC:

Bayesian information criterion

BMI:

Body mass index

CID:

Clinically important difference

GBTM:

Group-based trajectory modeling

GCPS:

Graded Chronic Pain Scale

NHB:

Non-Hispanic Black

NHW:

Non-Hispanic White

NRS:

Numerical rating scale

OA:

Osteoarthritis

WOMAC:

Western Ontario and McMaster Universities Osteoarthritis Index

References

  1. Lespasio MJ, Piuzzi NS, Husni ME, Muschler GF, Guarino A, Mont MA. Knee Osteoarthritis: A Primer. Perm J. 2017;21:16–183 2017/10/17. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5638628/pdf/16-183.pdf.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Hadler NM. Knee pain is the malady--not osteoarthritis. Ann Intern Med. 1992;116(7):598–9 1992/04/01.

    Article  CAS  PubMed  Google Scholar 

  3. Lawrence RC, Felson DT, Helmick CG, Arnold LM, Choi H, Deyo RA, et al. Estimates of the prevalence of arthritis and other rheumatic conditions in the United States. Part II. Arthritis Rheum. 2008;58(1):26–35 2008/01/01. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1002/art.23176.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Vos T, Flaxman AD, Naghavi M, Lozano R, Michaud C, Ezzati M, et al. Years lived with disability (YLDs) for 1160 sequelae of 289 diseases and injuries 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380(9859):2163–96 2012/12/19.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Hootman JM, Helmick CG, Brady TJ. A public health approach to addressing arthritis in older adults: the Most common cause of disability. Am J Public Health. 2012;102(3):426–33. https://doi.org/10.2105/AJPH.2011.300423.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Paradowski PT, Englund M, Lohmander LS, Roos EM. The effect of patient characteristics on variability in pain and function over two years in early knee osteoarthritis. Health Qual Life Outcomes. 2005;3(1):59 [cited 2021 Feb 11]. Available from: http://hqlo.biomedcentral.com/articles/10.1186/1477-7525-3-59.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Edwards RR, Cahalan C, Mensing G, Smith M, Haythornthwaite JA. Pain, catastrophizing, and depression in the rheumatic diseases. Nat Rev Rheumatol. 2011;7(4):216–24 2011/02/02.

    Article  CAS  PubMed  Google Scholar 

  8. Finan PH, Buenaver LF, Bounds SC, Hussain S, Park RJ, Haque UJ, et al. Discordance between pain and radiographic severity in knee osteoarthritis: findings from quantitative sensory testing of central sensitization. Arthritis Rheum. 2013;65(2):363–72 2012/09/11. Available from: https://www.ncbi.nlm.nih.gov/pubmed/22961435.

    Article  PubMed  Google Scholar 

  9. Cruz-Almeida Y, Cardoso J, Riley JL 3rd, Goodin B, King CD, Petrov M, et al. Physical performance and movement-evoked pain profiles in community-dwelling individuals at risk for knee osteoarthritis. Exp Gerontol. 2017;98:186–91 2017/08/27.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Cardoso JS, Riley JL 3rd, Glover T, Sibille KT, Bartley EJ, Goodin BR, et al. Experimental pain phenotyping in community-dwelling individuals with knee osteoarthritis. Pain. 2016;157(9):2104–14 Available from: https://www.ncbi.nlm.nih.gov/pubmed/27340911.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cruz-Almeida Y, King CD, Goodin BR, Sibille KT, Glover TL, Riley JL, et al. Psychological profiles and pain characteristics of older adults with knee osteoarthritis. Arthritis Care Res (Hoboken). 2013;65(11):1786–94. Available from: https://pubmed.ncbi.nlm.nih.gov/23861288. https://doi.org/10.1002/acr.22070.

    Article  Google Scholar 

  12. Emrani PS, Katz JN, Kessler CL, Reichmann WM, Wright EA, McAlindon TE, et al. Joint space narrowing and Kellgren–Lawrence progression in knee osteoarthritis: an analytic literature synthesis. Osteoarthr Cartil. 2008;16(8):873–82. Available from: http://www.sciencedirect.com/science/article/pii/S1063458407003925. https://doi.org/10.1016/j.joca.2007.12.004.

    Article  CAS  Google Scholar 

  13. Felson DT, Zhang Y, Hannan MT, Naimark A, Weissman BN, Aliabadi P, et al. The incidence and natural history of knee osteoarthritis in the elderly. The Framingham Osteoarthritis Study. Arthritis Rheum. 1995;38(10):1500–5 1995/10/01.

    Article  CAS  PubMed  Google Scholar 

  14. Peters TJ, Sanders C, Dieppe P, Donovan J. Factors associated with change in pain and disability over time: a community-based prospective observational study of hip and knee osteoarthritis. Br J Gen Pr. 2005;55(512):205–11 2005/04/06.

    Google Scholar 

  15. van Dijk GM, Dekker J, Veenhof C, van den Ende CHM. Course of functional status and pain in osteoarthritis of the hip or knee: a systematic review of the literature. Arthritis care res (Hoboken). 2006;55(5):779–85. Available from. https://doi.org/10.1002/art.22244.

    Article  Google Scholar 

  16. Collins JE, Katz JN, Dervan EE, Losina E. Trajectories and risk profiles of pain in persons with radiographic, symptomatic knee osteoarthritis: data from the osteoarthritis initiative. Osteoarthr Cartil. 2014;22(5):622–30 2014/03/26.

    Article  CAS  Google Scholar 

  17. Dieppe PA, Lohmander LS. Pathogenesis and management of pain in osteoarthritis. Lancet. 2005;365(9463):965–73. Available from: http://www.sciencedirect.com/science/article/pii/S0140673605710862. https://doi.org/10.1016/S0140-6736(05)71086-2.

    Article  CAS  PubMed  Google Scholar 

  18. Felson D, Niu J, Sack B, Aliabadi P, McCullough C, Nevitt MC. Progression of osteoarthritis as a state of inertia. Ann Rheum Dis. 2013;72(6):924 Available from: http://ard.bmj.com/content/72/6/924.abstract.

    Article  PubMed  Google Scholar 

  19. Nicholls E, Thomas E, van der Windt DA, Croft PR, Peat G. Pain trajectory groups in persons with, or at high risk of, knee osteoarthritis: findings from the Knee Clinical Assessment Study and the Osteoarthritis Initiative. Osteoarthr Cartil. 2014;22(12):2041–50 2014/10/12.

    Article  CAS  Google Scholar 

  20. White DK, Neogi T, Nguyen USDT, Niu J, Zhang Y. Trajectories of functional decline in knee osteoarthritis: The osteoarthritis initiative. Rheumatol (United Kingdom). 2016;55(5):801–8 [cited 2021 Feb 9]. Available from: /pmc/articles/PMC5009418/.

    Google Scholar 

  21. Vina ER, Ran D, Ashbeck EL, Kwoh CK. Natural history of pain and disability among African-Americans and Whites with or at risk for knee osteoarthritis: A longitudinal study. Osteoarthr Cartil. 2018;26(4):471–9 2018/02/07.

    Article  CAS  Google Scholar 

  22. Radojčić MR, Arden NK, Yang X, Strauss VY, Birrell F, Cooper C, et al. Pain trajectory defines knee osteoarthritis subgroups: a prospective observational study. Pain. 2020;161(12):2841–51 [cited 2021 Feb 9]. Available from: https://journals.lww.com/10.1097/j.pain.0000000000001975.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Dai Z, Lu N, Niu J, Felson DT, Zhang Y. Dietary Fiber Intake in Relation to Knee Pain Trajectory. Arthritis Care Res. 2017;69(9):1331–9 [cited 2021 Feb 11]. Available from: http://doi.wiley.com/10.1002/acr.23158.

    Article  CAS  Google Scholar 

  24. Schiphof D, Runhaar J, Waarsing JH, van Spil WE, van Middelkoop M, Bierma-Zeinstra SMA. The clinical and radiographic course of early knee and hip osteoarthritis over 10 years in CHECK (cohort hip and cohort knee). Osteoarthr Cartil. 2019;27(10):1491–500. https://doi.org/10.1016/j.joca.2019.06.002.

    Article  CAS  Google Scholar 

  25. Bastick AN, Wesseling J, Damen J, Verkleij SPJ, Emans PJ, Bindels PJE, et al. Defining knee pain trajectories in early symptomatic knee osteoarthritis in primary care: 5-year results from a nationwide prospective cohort study (CHECK). Br J Gen Pract. 2016;66(642):e32–9 [cited 2021 Feb 11]. Available from: https://bjgp.org/content/66/642/e32.

    Article  PubMed  Google Scholar 

  26. Trouvin A-P, Marty M, Goupille P, Perrot S. Determinants of daily pain trajectories and relationship with pain acceptability in hip and knee osteoarthritis. A national prospective cohort study on 886 patients. Jt Bone Spine. 2019;86:245–50. Available from. https://doi.org/10.1016/j.jbspin.2018.06.009.

    Article  Google Scholar 

  27. Mills K, Eyles JP, Martin MA, Hancock MJ, Hunter DJ. Exploratory study of 6-month pain trajectories in individuals with predominant patellofemoral osteoarthritis: A cohort study. J Orthop Sports Phys Ther. 2019;49(1):5–16 [cited 2021 Feb 9]. Available from: https://www.jospt.org/doi/abs/10.2519/jospt.2019.8354.

    Article  PubMed  Google Scholar 

  28. Alschuler KN, Molton IR, Jensen MP, Riddle DL. Prognostic value of coping strategies in a community-based sample of persons with chronic symptomatic knee osteoarthritis. Pain. 2013;154(12):2775–81 [cited 2020 Sep 11]. Available from: /pmc/articles/PMC4298486/?report=abstract.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Tanamas SK, Wluka AE, Davies-Tuck M, Wang Y, Strauss BJ, Proietto J, et al. Association of weight gain with incident knee pain, stiffness, and functional difficulties: A longitudinal study. Arthritis Care Res. 2013;65(1):34–43 [cited 2021 Feb 11]. Available from: https://pubmed.ncbi.nlm.nih.gov/22674832/.

    Article  Google Scholar 

  30. Mallen CD, Peat G, Thomas E, Lacey R, Croft P. Predicting poor functional outcome in community-dwelling older adults with knee pain: Prognostic value of generic indicators. Ann Rheum Dis. 2007;66(11):1456–61 [cited 2021 Feb 11]. Available from: https://pubmed.ncbi.nlm.nih.gov/17456527/.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Dunlop DD, Semanik P, Song J, Manheim LM, Shih V, Chang RW. Risk factors for functional decline in older adults with arthritis. Arthritis Rheum. 2005;52(4):1274–82 [cited 2021 Feb 11]. Available from: https://pubmed.ncbi.nlm.nih.gov/15818691/.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Belo JN, Berger MY, Koes BW. Bierma-Zeinstra SMA. Prognostic factors in adults with knee pain in general practice. Arthritis Rheum. 2009;61(2):143–51 [cited 2021 Jan 28]. Available from: http://doi.wiley.com/10.1002/art.24419.

    Article  CAS  PubMed  Google Scholar 

  33. Belo JN, Berger MY, Reijman M, Koes BW, Bierma-Zeinstra SMA. Prognostic factors of progression of osteoarthritis of the knee: A systematic review of observational studies [Internet]. Vol. 57, Arthritis Care and Research. Arthritis Rheum. 2007:13–26 [cited 2021 Feb 11]. Available from: https://pubmed.ncbi.nlm.nih.gov/17266080/.

  34. Previtali D, Andriolo L, Di Laura FG, Boffa A, Candrian C, Zaffagnini S, et al. Pain Trajectories in Knee Osteoarthritis—A Systematic Review and Best Evidence Synthesis on Pain Predictors. J Clin Med. 2020;9(9):2828 [cited 2020 Oct 5]. Available from: https://www.mdpi.com/2077-0383/9/9/2828.

    Article  PubMed Central  Google Scholar 

  35. Vaughn IA, Terry EL, Bartley EJ, Schaefer N, Fillingim RB. Racial-Ethnic Differences in Osteoarthritis Pain and Disability: A Meta-Analysis. J Pain. 2018:629-44. 2018/12/14.

  36. Booker S, Cardoso J, Cruz-Almeida Y, Sibille KT, Terry EL, Powell-Roach KL, et al. Movement-evoked pain, physical function, and perceived stress: An observational study of ethnic/racial differences in aging non-Hispanic Blacks and non-Hispanic Whites with knee osteoarthritis. Exp Gerontol. 2019;124:110622 [cited 2020 Sep 22]. Available from: https://pubmed.ncbi.nlm.nih.gov/31154005/.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Johnson AJ, Sibille KT, Cardoso J, Terry EL, Powell-Roach KL, Goodin B, et al. Patterns and correlates of self-management strategies for osteoarthritis related pain among older non-Hispanic Black and non-Hispanic White adults. Arthritis Care Res (Hoboken). 2020. Available from. https://doi.org/10.1002/acr.24396.

  38. Cruz-Almeida Y, Sibille KT, Goodin BR, Petrov ME, Bartley EJ, Riley JL 3rd, et al. Racial and ethnic differences in older adults with knee osteoarthritis. Arthritis Rheumatol. 2014;66(7):1800–10 2014/04/15. Available from: https://www.ncbi.nlm.nih.gov/pubmed/24729357.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Terry EL, Booker SQ, Cardoso JS, Sibille KT, Bartley EJ, Glover TL, et al. Neuropathic-Like Pain Symptoms in a Community-Dwelling Sample with or at Risk for Knee Osteoarthritis. Pain Med. 2019;21(1):125–37 2019/06/01.

    Article  PubMed Central  Google Scholar 

  40. Bartley EJ, Hossain NI, Gravlee CC, Sibille KT, Terry EL, Vaughn IA, et al. Race/ethnicity moderates the association between psychosocial resilience and movement-evoked pain in knee osteoarthritis. ACR Open Rheumatol. 2019;1(1):16–25. Available from. https://doi.org/10.1002/acr2.1002.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Altman R, Asch E, Bloch D, Bole G, Borenstein D, Brandt K, et al. Development of criteria for the classification and reporting of osteoarthritis. Classification of osteoarthritis of the knee. Diagnostic and Therapeutic Criteria Committee of the American Rheumatism Association. Arthritis Rheum. 1986;29(8):1039–49 1986/08/01.

    Article  CAS  PubMed  Google Scholar 

  42. Thompson KA, Bulls HW, Sibille KT, Bartley EJ, Glover TL, Terry EL, et al. Optimism and Psychological Resilience are Beneficially Associated With Measures of Clinical and Experimental Pain in Adults With or at Risk for Knee Osteoarthritis. Clin J Pain. 2018;34(12):1164–72 2018/07/24.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Roux CH, Saraux A, Mazieres B, Pouchot J, Morvan J, Fautrel B, et al. Screening for hip and knee osteoarthritis in the general population: predictive value of a questionnaire and prevalence estimates. Ann Rheum Dis. 2008;67(10):1406–11 2007/12/14.

    Article  CAS  PubMed  Google Scholar 

  44. Wright RW. Osteoarthritis Classification Scales: Interobserver Reliability and Arthroscopic Correlation. J Bone Jt Surg. 2014;96(14):1145–51 [cited 2021 Feb 11]. Available from: https://journals.lww.com/00004623-201407160-00001.

    Article  Google Scholar 

  45. Sibille KT, Chen H, Bartley EJ, Riley J, Glover TL, King CD, et al. Accelerated aging in adults with knee osteoarthritis pain: consideration for frequency, intensity, time, and total pain sites. Pain Reports. 2017;2(3):e591 Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5741297/.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Sibille KT, Steingrímsdóttir ÓA, Fillingim RB, Stubhaug A, Schirmer H, Chen H, et al. Investigating the burden of chronic pain: an inflammatory and metabolic composite. Pain Res Manag. 2016;2016:7657329. Craig KD, editor. P. Available from. https://doi.org/10.1155/2016/7657329.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Von Korff M, Ormel J, Keefe FJ, Dworkin SF. Grading the severity of chronic pain. Pain. 1992;50(2):133–49 1992/08/01.

    Article  Google Scholar 

  48. Alghadir AH, Anwer S, Iqbal A, Iqbal ZA. Test-retest reliability, validity, and minimum detectable change of visual analog, numerical rating, and verbal rating scales for measurement of osteoarthritic knee pain. J Pain Res. 2018;11:851–6 Available from: https://www.ncbi.nlm.nih.gov/pubmed/29731662.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Ferreira-Valente MA, Pais-Ribeiro JL, Jensen MP. Validity of four pain intensity rating scales. Pain. 2011;152(10):2399–404 2011/08/23.

    Article  PubMed  Google Scholar 

  50. Bellamy N, Buchanan WW, Goldsmith CH, Campbell J, Stitt LW. Validation study of WOMAC: a health status instrument for measuring clinically important patient relevant outcomes to antirheumatic drug therapy in patients with osteoarthritis of the hip or knee. J Rheumatol. 1988;15(12):1833–40 1988/12/01.

    CAS  PubMed  Google Scholar 

  51. Theiler R, Spielberger J, Bischoff HA, Bellamy N, Huber J, Kroesen S. Clinical evaluation of the WOMAC 3.0 OA Index in numeric rating scale format using a computerized touch screen version. Osteoarthr Cartil. 2002;10(6):479–81 Available from: https://pubmed.ncbi.nlm.nih.gov/12056851.

    Article  CAS  Google Scholar 

  52. Jones BL, Nagin DS, Roeder K. A SAS procedure based on mixture models for estimating developmental trajectories. Sociol Methods Res. 2001;29(3):374–93. https://doi.org/10.1177/0049124101029003005.

    Article  Google Scholar 

  53. Nagin DS. Analyzing developmental trajectories: a semiparametric, group-based approach. Psychol Methods. 1999;4(2):139–57. https://doi.org/10.1037/1082-989X.4.2.139.

    Article  Google Scholar 

  54. Jung T, Wickrama KAS. An introduction to latent class growth analysis and growth mixture modeling. Soc Personal Psychol Compass. 2008;2(1):302–17. https://doi.org/10.1111/j.1751-9004.2007.00054.x.

    Article  Google Scholar 

  55. Nagin DS, Odgers CL. Group-based trajectory modeling in clinical research. Annu Rev Clin Psychol. 2010;6:109–38 2010/03/03.

    Article  PubMed  Google Scholar 

  56. Vittinghoff E, McCulloch CE, Glidden DV, Shiboski SC. 3 - Linear and Non-Linear Regression Methods in Epidemiology and Biostatistics. In: Rao CR, Miller JP, Rao DC, editors. Essential Statistical Methods for Medical Statistics. Boston: North-Holland; 2011. p. 66–103. Available from: http://www.sciencedirect.com/science/article/pii/B9780444537379500062.

    Chapter  Google Scholar 

  57. van de Schoot R, Sijbrandij M, Winter SD, Depaoli S, Vermunt JK. The GRoLTS-Checklist: Guidelines for Reporting on Latent Trajectory Studies. Struct Equ Model A Multidiscip J. 2017;24(3):451–67 [cited 2021 Feb 9]. Available from: https://www.tandfonline.com/doi/full/10.1080/10705511.2016.1247646.

    Article  Google Scholar 

  58. Wieczorek M, Rotonda C, Guillemin F, Rat A-C. What have we learned from trajectory analysis of clinical outcomes in knee and hip osteoarthritis before surgery? Arthritis Care Res (Hoboken). 2019. Available from. https://doi.org/10.1002/acr.24069.

  59. Wieczorek M, Rotonda C, Coste J, Pouchot J, Saraux A, Guillemin F, et al. Trajectory analysis combining pain and physical function in individuals with knee and hip osteoarthritis: results from the French KHOALA cohort. Rheumatol. 2020:3488-98 2020/05/07.

  60. Wesseling J, Bastick AN, ten Wolde S, Kloppenburg M, Lafeber FPJG, Bierma-Zeinstra SMA, et al. Identifying trajectories of pain severity in early symptomatic knee osteoarthritis: a 5-year Followup of the cohort hip and cohort knee (CHECK) study. J Rheumatol. 2015;42(8):1470–7. Available from: http://europepmc.org/abstract/MED/26136492. https://doi.org/10.3899/jrheum.141036.

    Article  PubMed  Google Scholar 

  61. Lillis TA, Burns J, Aranda F, Burgess HJ, Purim-Shem-Tov YA, Bruehl S, et al. Race-related differences in acute pain complaints among inner-city women: the role of socioeconomic status. J Behav Med. 2019; Available from. https://doi.org/10.1007/s10865-019-00123-3.

  62. Hill CV, Perez-Stable EJ, Anderson NA, Bernard MA. The National Institute on Aging Health Disparities Research Framework. Ethn Dis. 2015;25(3):245–54 2015/12/18.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Arendt-Nielsen L, Nie H, Laursen MB, Laursen BS, Madeleine P, Simonsen OH, et al. Sensitization in patients with painful knee osteoarthritis. Pain. 2010;149(3):573–81 2010/04/27.

    Article  PubMed  Google Scholar 

  64. Reynolds Losin EA, Woo C-W, Medina NA, Andrews-Hanna JR, Eisenbarth H, Wager TD. Neural and sociocultural mediators of ethnic differences in pain. Available from; 2020. https://doi.org/10.1038/s41562-020-0819-8.

    Book  Google Scholar 

  65. Terry EL, Tanner JJ, Cardoso JS, Sibille KT, Lai S, Deshpande H, et al. Associations of pain catastrophizing with pain-related brain structure in individuals with or at risk for knee osteoarthritis: Sociodemographic considerations. Available from; 2020. https://doi.org/10.1007/s11682-020-00372-w.

    Book  Google Scholar 

  66. Mickle AM, Garvan C, Service C, Pop R, Marks J, Wu S, et al. Relationships Between Pain, Life Stress, Sociodemographics, and Cortisol: Contributions of Pain Intensity and Financial Satisfaction. Chronic Stress. 4:1–12 Available from: https://us.sagepub.com/en-.

  67. Verkleij SPJ, Hoekstra T, Rozendaal RM, Waarsing JH, Koes BW, Luijsterburg PAJ, et al. Defining discriminative pain trajectories in hip osteoarthritis over a 2-year time period. Ann Rheum Dis. 2012;71(9):1517 LP–1523 Available from: http://ard.bmj.com/content/71/9/1517.abstract.

    Article  Google Scholar 

  68. Dabare C, Le Marshall K, Leung A, Page CJ, Choong PF, Lim KK. Differences in presentation, progression and rates of arthroplasty between hip and knee osteoarthritis: Observations from an osteoarthritis cohort study-a clear role for conservative management. Int J Rheum Dis. 2017;20(10):1350–60 [cited 2021 Feb 12]. Available from: /pmc/articles/PMC5655735/.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Bastick AN, Wesseling J, Damen J, Verkleij SPJJ, Emans PJ, Bindels PJEE, et al. Defining knee pain trajectories in early symptomatic knee osteoarthritis in primary care: 5-year results from a nationwide prospective cohort study (CHECK). Br J Gen Pract. 2016;66(642):e32–9 2015/12/06. [cited 2020 Sep 11]. Available from: /pmc/articles/PMC4684033/?report=abstract.

    Article  PubMed  Google Scholar 

  70. Benyon K, Hill S, Zadurian N, Mallen C. Coping strategies and self-efficacy as predictors of outcome in osteoarthritis: a systematic review. Musculoskeletal Care. 2010;8(4):224–36 2010/10/22.

    Article  PubMed  Google Scholar 

  71. Kolasinski SL, Neogi T, Hochberg MC, Oatis C, Guyatt G, Block J, et al. 2019 American College of Rheumatology/Arthritis Foundation guideline for the Management of Osteoarthritis of the hand, hip, and knee. Arthritis Rheumatol. 2020;72(2):220–33. Available from. https://doi.org/10.1002/art.41142.

    Article  PubMed  Google Scholar 

  72. Nagin D. Group-based modeling of development: Harvard University Press; 2009.

    Google Scholar 

  73. Stone AA, Schwartz JE, Broderick JE, Shiffman SS. Variability of momentary pain predicts recall of weekly pain: A consequence of the peak (or salience) memory heuristic [Internet]. Vol. 31, Personality and Social Psychology Bulletin. Sage PublicationsSage CA: 2005;1340–6 [cited 2020 Oct 5]. Available from: http://journals.sagepub.com/doi/10.1177/0146167205275615

  74. van der Nest G, Lima Passos V, Candel MJJM, van Breukelen GJP. An overview of mixture modelling for latent evolutions in longitudinal data: modelling approaches, fit statistics and software. Adv Life Course Res. 2020;43:100323. https://doi.org/10.1016/j.alcr.2019.100323.

    Article  Google Scholar 

Download references

Funding

Research funding and support provided by NIH/NIA Grants R01AG054370 (KTS); R37AG033906 (RBF) and UF CTSA Grant UL1TR001427 and UAB CTSA Grant UL1TR001417 from the NIH Center for Advancing Translational Sciences; NIH/NIAMS K23AR076463–01 and NIH/NIMHD L60MD014496 (SQB); NIH/NINDS 4K22NS102334–03 and NIH/NIA 1P30AG059297 provided to the University of Florida and McKnight Brain Institute Career Development Award (ELT), minority supplement provided to the University of Florida (JSC). Data management was funded in part by NCATS grant UL1TR000064. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Author information

Authors and Affiliations

Authors

Contributions

AJJ made substantial contributions to the authorship of this manuscript by participating in the conception and design, analysis and interpretation of data, drafting the article, critically revising the article, provided final approval of the article, and provided administrative support. TV made substantial contributions to the authorship of this manuscript by participating in the conception and design, analysis and interpretation of data, drafting the article, critically revising the article, provided final approval of the article, and statistical expertise. SQB made substantial contributions to the authorship of this manuscript by participating in the conception and design, interpretation of the data, drafting the article, critically revising the article, provided final approval of the article, and obtained funding. JC made substantial contributions to the authorship of this manuscript by providing administrative and technical logistic support, collecting and assembling data, critically revising the article, and provided final approval of the article. ELT made substantial contributions to the authorship of this manuscript by participating in involvement in study implementation and data collection, interpretation of the data, critically revising the article, and provided final approval of the article. KPR made substantial contributions to the authorship of this manuscript by involvement in study implementation and data collection, critically revising the article, and provided final approval of the article. RS made substantial contributions to the authorship of this manuscript by participating in overall study design, data collection oversight, interpretation of data, critically revising the article, and provided final approval of the article. DK made substantial contributions to the authorship of this manuscript by involvement in study implementation and data collection, contributions related to conception and interpretation of data, critically revising the article, and provided final approval of the article. AA made substantial contributions to the authorship of this manuscript by providing administrative and technical logistic support, collecting and assembling data, critically revising the article, and provided final approval of the article. DTR made substantial contributions to the authorship of this manuscript by involvement in overall study design and data collection and analysis oversight, interpretation of the data, and final approval of the article. BRG made substantial contributions to the authorship of this manuscript by involvement in the overall study design and oversight, interpretation of data, critically revising the article, provided final approval of the article, and obtained funding. RBF made substantial contributions to the authorship of this manuscript by obtaining funding, responsibility for overall study design and oversight, participating in the conception and design, interpretation of data, critically revising the article, and provided final approval of the article. KTS made substantial contributions to the authorship of this manuscript by obtaining funding, responsibility for overall study design and oversight, participating in the conception and design, interpretation of data, drafting the article, critically revising the article, and provided final approval of the article. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Alisa J. Johnson.

Ethics declarations

Ethics approval and consent to participate

All study procedures were approved by the Institutional Review Boards (IRBs) at the University of Florida and the University of Alabama, Birmingham.

Consent for publication

All participants provided written informed consent prior to study involvement.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Johnson, A.J., Vasilopoulos, T., Booker, S.Q. et al. Knee pain trajectories over 18 months in non-Hispanic Black and non-Hispanic White adults with or at risk for knee osteoarthritis. BMC Musculoskelet Disord 22, 415 (2021). https://doi.org/10.1186/s12891-021-04284-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-021-04284-8

Keywords