Trends for opioid prescribing and the impact of the COVID-19 pandemic in patients with rheumatic and musculoskeletal diseases between 2006 and 2021.

OBJECTIVE: To investigate opioid prescribing trends and assess the impact of the COVID-19 pandemic on opioid prescribing in rheumatic and musculoskeletal diseases (RMDs). METHODS: Adult patients with RA, PsA, axial spondyloarthritis (AxSpA), SLE, OA and FM with opioid prescriptions between 1 January 2006 and 31 August 2021 without cancer in UK primary care were included. Age- and gender-standardized yearly rates of new and prevalent opioid users were calculated between 2006 and 2021. For prevalent users, monthly measures of mean morphine milligram equivalents (MME)/day were calculated between 2006 and 2021. To assess the impact of the pandemic, we fitted regression models to the monthly number of prevalent opioid users between January 2015 and August 2021. The time coefficient reflects the trend pre-pandemic and the interaction term coefficient represents the change in the trend during the pandemic. RESULTS: The study included 1 313 519 RMD patients. New opioid users for RA, PsA and FM increased from 2.6, 1.0 and 3.4/10 000 persons in 2006 to 4.5, 1.8 and 8.7, respectively, in 2018 or 2019. This was followed by a fall to 2.4, 1.2 and 5.9, respectively, in 2021. Prevalent opioid users for all RMDs increased from 2006 but plateaued or dropped beyond 2018, with a 4.5-fold increase in FM between 2006 and 2021. In this period, MME/day increased for all RMDs, with the highest for FM (≥35). During COVID-19 lockdowns, RA, PsA and FM showed significant changes in the trend of prevalent opioid users. The trend for FM increased pre-pandemic and started decreasing during the pandemic. CONCLUSION: The plateauing or decreasing trend of opioid users for RMDs after 2018 may reflect the efforts to tackle rising opioid prescribing in the UK. The pandemic led to fewer people on opioids for most RMDs, providing reassurance that there was no sudden increase in opioid prescribing during the pandemic.

The epidemiology of psoriatic arthritis in the UK: a health intelligence analysis of UK Primary Care Electronic Health Records 1991-2020.

OBJECTIVES: Epidemiological estimates of psoriatic arthritis (PsA) underpin the provision of healthcare, research, and the work of government, charities and patient organizations. Methodological problems impacting prior estimates include small sample sizes, incomplete case ascertainment, and representativeness. We developed a statistical modelling strategy to provide contemporary prevalence and incidence estimates of PsA from 1991 to 2020 in the UK. METHODS: Data from Clinical Practice Research Datalink (CPRD) were used to identify cases of PsA between 1st January 1991 and 31st December 2020. To optimize ascertainment, we identified cases of Definite PsA (≥1 Read code for PsA) and Probable PsA (satisfied a bespoke algorithm). Standardized annual rates were calculated using Bayesian multilevel regression with post-stratification to account for systematic differences between CPRD data and the UK population, based on age, sex, socioeconomic status and region of residence. RESULTS: A total of 26293 recorded PsA cases (all definitions) were identified within the study window (77.9% Definite PsA). Between 1991 and 2020 the standardized prevalence of PsA increased twelve-fold from 0.03 to 0.37. The standardized incidence of PsA per 100,000 person years increased from 8.97 in 1991 to 15.08 in 2020, an almost 2-fold increase. Over time, rates were similar between the sexes, and across socioeconomic status. Rates were strongly associated with age, and consistently highest in Northern Ireland. CONCLUSION: The prevalence and incidence of PsA recorded in primary care has increased over the last three decades. The modelling strategy presented can be used to provide contemporary prevalence estimates for musculoskeletal disease using routinely collected primary care data.

Modelling and classifying joint trajectories of self-reported mood and pain in a large cohort study.

It is well-known that mood and pain interact with each other, however individual-level variability in this relationship has been less well quantified than overall associations between low mood and pain. Here, we leverage the possibilities presented by mobile health data, in particular the "Cloudy with a Chance of Pain" study, which collected longitudinal data from the residents of the UK with chronic pain conditions. Participants used an App to record self-reported measures of factors including mood, pain and sleep quality. The richness of these data allows us to perform model-based clustering of the data as a mixture of Markov processes. Through this analysis we discover four endotypes with distinct patterns of co-evolution of mood and pain over time. The differences between endotypes are sufficiently large to play a role in clinical hypothesis generation for personalised treatments of comorbid pain and low mood.

"Take up to eight tablets per day": Incorporating free-text medication instructions into a transparent and reproducible process for preparing drug exposure data for pharmacoepidemiology.

PURPOSE: Routinely collected prescription data provides drug exposure information for pharmacoepidemiology, informing start/stop dates and dosage. Prescribing information includes structured data and unstructured free-text instructions, which can include inherent variability, such as "one to two tablets up to four times a day". Preparing drug exposure data from raw prescriptions to a research ready dataset is rarely fully reported, yet assumptions have considerable implications for pharmacoepidemiology. This may have bigger consequences for "pro re nata" (PRN) drugs. Our aim was, using a worked example of opioids and fracture risk, to examine the impact of incorporating narrative prescribing instructions and subsequent drug preparation assumptions on adverse event rates. METHODS: R-packages for extracting free-text medication prescription instructions in a structured form (doseminer) and an algorithm for transparently processing drug exposure information (drugprepr) were developed. Clinical Practice Research Datalink GOLD was used to define a cohort of adult new opioid users without prior cancer. A retrospective cohort study was performed using data between January 1, 2017 and July 31, 2018. We tested the impact of varying drug preparation assumptions by estimating the risk of opioids on fracture risk using Cox proportional hazards models. RESULTS: During the study window, 60 394 patients were identified with 190 754 opioid prescriptions. Free-text prescribing instruction variability, where there was flexibility in the number of tablets to be administered, was present in 42% prescriptions. Variations in the decisions made during preparing raw data for analysis led to marked differences impacting the event number (n = 303-415) and person years of drug exposure (5619-9832). The distribution of hazard ratios as a function of the decisions ranged from 2.71 (95% CI: 2.31, 3.18) to 3.24 (2.76, 3.82). CONCLUSIONS: Assumptions made during the drug preparation process, especially for those with variability in prescription instructions, can impact results of subsequent risk estimates. The developed R packages can improve transparency related to drug preparation assumptions, in line with best practice advocated by international pharmacoepidemiology guidelines.

Sleep Disturbance and Quality of Life in Rheumatoid Arthritis: Prospective mHealth Study.

BACKGROUND: Sleep disturbances and poor health-related quality of life (HRQoL) are common in people with rheumatoid arthritis (RA). Sleep disturbances, such as less total sleep time, more waking periods after sleep onset, and higher levels of nonrestorative sleep, may be a driver of HRQoL. However, understanding whether these sleep disturbances reduce HRQoL has, to date, been challenging because of the need to collect complex time-varying data at high resolution. Such data collection is now made possible by the widespread availability and use of mobile health (mHealth) technologies. OBJECTIVE: This mHealth study aimed to test whether sleep disturbance (both absolute values and variability) causes poor HRQoL. METHODS: The quality of life, sleep, and RA study was a prospective mHealth study of adults with RA. Participants completed a baseline questionnaire, wore a triaxial accelerometer for 30 days to objectively assess sleep, and provided daily reports via a smartphone app that assessed sleep (Consensus Sleep Diary), pain, fatigue, mood, and other symptoms. Participants completed the World Health Organization Quality of Life-Brief (WHOQoL-BREF) questionnaire every 10 days. Multilevel modeling tested the relationship between sleep variables and the WHOQoL-BREF domains (physical, psychological, environmental, and social). RESULTS: Of the 268 recruited participants, 254 were included in the analysis. Across all WHOQoL-BREF domains, participants' scores were lower than the population average. Consensus Sleep Diary sleep parameters predicted the WHOQoL-BREF domain scores. For example, for each hour increase in the total time asleep physical domain scores increased by 1.11 points (ß=1.11, 95% CI 0.07-2.15) and social domain scores increased by 1.65 points. These associations were not explained by sociodemographic and lifestyle factors, disease activity, medication use, anxiety levels, sleep quality, or clinical sleep disorders. However, these changes were attenuated and no longer significant when pain, fatigue, and mood were included in the model. Increased variability in total time asleep was associated with poorer physical and psychological domain scores, independent of all covariates. There was no association between actigraphy-measured sleep and WHOQoL-BREF. CONCLUSIONS: Optimizing total sleep time, increasing sleep efficiency, decreasing sleep onset latency, and reducing variability in total sleep time could improve HRQoL in people with RA.

A semi-parametric mixed models for longitudinally measured fasting blood sugar level of adult diabetic patients.

BACKGROUND: At the diabetic clinic of Jimma University Specialized Hospital, health professionals provide regular follow-up to help people with diabetes live long and relatively healthy lives. Based on patient condition, they also provide interventions in the form of counselling to promote a healthy diet and physical activity and prescribing medicines. The main purpose of this study is to estimate the rate of change of fasting blood sugar (FBS) profile experienced by patients over time. The change may help to assess the effectiveness of interventions taken by the clinic to regulate FBS level, where rates of change close to zero over time may indicate the interventions are good regulating the level. METHODS: In the analysis of longitudinal data, the mean profile is often estimated by parametric linear mixed effects model. However, the individual and mean profile plots of FBS level for diabetic patients are nonlinear and imposing parametric models may be too restrictive and yield unsatisfactory results. We propose a semi-parametric mixed model, in particular using spline smoothing to efficiently analyze a longitudinal measured fasting blood sugar level of adult diabetic patients accounting for correlation between observations through random effects. RESULTS: The semi-parametric mixed models had better fit than the linear mixed models for various variance structures of subject-specific random effects. The study revealed that the rate of change in FBS level in diabetic patients, due to the clinic interventions, does not continue as a steady pace but changes with time and weight of patients. CONCLUSIONS: The proposed method can help a physician in clinical monitoring of diabetic patients and to assess the effect of intervention packages, such as healthy diet, physical activity and prescribed medicines, because individualized curve may be obtained to follow patient-specific FBS level trends.