OpenSAFELY: A platform for analysing electronic health records designed for reproducible research.

Electronic health records (EHRs) and other administrative health data are increasingly used in research to generate evidence on the effectiveness, safety, and utilisation of medical products and services, and to inform public health guidance and policy. Reproducibility is a fundamental step for research credibility and promotes trust in evidence generated from EHRs. At present, ensuring research using EHRs is reproducible can be challenging for researchers. Research software platforms can provide technical solutions to enhance the reproducibility of research conducted using EHRs. In response to the COVID-19 pandemic, we developed the secure, transparent, analytic open-source software platform OpenSAFELY designed with reproducible research in mind. OpenSAFELY mitigates common barriers to reproducible research by: standardising key workflows around data preparation; removing barriers to code-sharing in secure analysis environments; enforcing public sharing of programming code and codelists; ensuring the same computational environment is used everywhere; integrating new and existing tools that encourage and enable the use of reproducible working practices; and providing an audit trail for all code that is run against the real data to increase transparency. This paper describes OpenSAFELY's reproducibility-by-design approach in detail.

Challenges in Estimating the Effectiveness of COVID-19 Vaccination Using Observational Data.

The COVID-19 vaccines were developed and rigorously evaluated in randomized trials during 2020. However, important questions, such as the magnitude and duration of protection, their effectiveness against new virus variants, and the effectiveness of booster vaccination, could not be answered by randomized trials and have therefore been addressed in observational studies. Analyses of observational data can be biased because of confounding and because of inadequate design that does not consider the evolution of the pandemic over time and the rapid uptake of vaccination. Emulating a hypothetical "target trial" using observational data assembled during vaccine rollouts can help manage such potential sources of bias. This article describes 2 approaches to target trial emulation. In the sequential approach, on each day, eligible persons who have not yet been vaccinated are matched to a vaccinated person. The single-trial approach sets a single baseline at the start of the rollout and considers vaccination as a time-varying variable. The nature of the confounding depends on the analysis strategy: Estimating "per-protocol" effects (accounting for vaccination of initially unvaccinated persons after baseline) may require adjustment for both baseline and "time-varying" confounders. These issues are illustrated by using observational data from 2 780 931 persons in the United Kingdom aged 70 years or older to estimate the effect of a first dose of a COVID-19 vaccine. Addressing the issues discussed in this article should help authors of observational studies provide robust evidence to guide clinical and policy decisions.

Describing the population experiencing COVID-19 vaccine breakthrough following second vaccination in England: a cohort study from OpenSAFELY.

BACKGROUND: While the vaccines against COVID-19 are highly effective, COVID-19 vaccine breakthrough is possible despite being fully vaccinated. With SARS-CoV-2 variants still circulating, describing the characteristics of individuals who have experienced COVID-19 vaccine breakthroughs could be hugely important in helping to determine who may be at greatest risk. METHODS: With the approval of NHS England, we conducted a retrospective cohort study using routine clinical data from the OpenSAFELY-TPP database of fully vaccinated individuals, linked to secondary care and death registry data and described the characteristics of those experiencing COVID-19 vaccine breakthroughs. RESULTS: As of 1st November 2021, a total of 15,501,550 individuals were identified as being fully vaccinated against COVID-19, with a median follow-up time of 149 days (IQR: ​107-179). From within this population, a total of 579,780 (<4%) individuals reported a positive SARS-CoV-2 test. For every 1000 years of patient follow-up time, the corresponding incidence rate (IR) was 98.06 (95% CI 97.93-98.19). There were 28,580 COVID-19-related hospital admissions, 1980 COVID-19-related critical care admissions and 6435 COVID-19-related deaths; corresponding IRs 4.77 (95% CI 4.74-4.80), 0.33 (95% CI 0.32-0.34) and 1.07 (95% CI 1.06-1.09), respectively. The highest rates of breakthrough COVID-19 were seen in those in care homes and in patients with chronic kidney disease, dialysis, transplant, haematological malignancy or who were immunocompromised. CONCLUSIONS: While the majority of COVID-19 vaccine breakthrough cases in England were mild, some differences in rates of breakthrough cases have been identified in several clinical groups. While it is important to note that these findings are simply descriptive and cannot be used to answer why certain groups have higher rates of COVID-19 breakthrough than others, the emergence of the Omicron variant of COVID-19 coupled with the number of positive SARS-CoV-2 tests still occurring is concerning and as numbers of fully vaccinated (and boosted) individuals increases and as follow-up time lengthens, so too will the number of COVID-19 breakthrough cases. Additional analyses, to assess vaccine waning and rates of breakthrough COVID-19 between different variants, aimed at identifying individuals at higher risk, are needed.

Using Bayesian networks to identify musculoskeletal symptoms influencing the risk of developing psoriatic arthritis in people with psoriasis.

OBJECTIVES: The aim of this study was to explore the use of Bayesian networks (BNs) to understand the relationships between musculoskeletal symptoms and the development of PsA in people with psoriasis. METHODS: Incident cases of psoriasis were identified for 1998 to 2015 from the UK Clinical Research Practice Datalink. Musculoskeletal symptoms (identified by Medcodes) were concatenated into primary groups, each made up of several subgroups. Baseline demographics for gender, age, BMI, psoriasis severity, alcohol use and smoking status were also extracted. Several BN structures were composed using a combination of expert knowledge and data-oriented modelling based on: (i) primary musculoskeletal symptom groups; (ii) musculoskeletal symptom subgroups and (iii) demographic variables. Predictive ability of the networks using the area under the receiver operating characteristic curve was calculated. RESULTS: Over one million musculoskeletal symptoms were extracted for the 90 189 incident cases of psoriasis identified, of which 1409 developed PsA. The BN analysis yielded direct relationships between gender, BMI, arthralgia, finger pain, fatigue, hand pain, hip pain, knee pain, swelling, back pain, myalgia and PsA. The best BN, achieved by using the more site-specific musculoskeletal symptom subgroups, was 76% accurate in predicting the development of PsA in a test set and had an area under the receiver operating characteristic curve of 0.73 (95% CI: 0.70, 0.75). CONCLUSION: The presented BN model may be a useful method to identify clusters of symptoms that predict the development of PsA with reasonable accuracy. Using a BN approach, we have shown that there are several symptoms which are predecessors of PsA, including fatigue, specific types of pain and swelling.

Gender stratified adjustment of the DAS28-CRP improves inter-score agreement with the DAS28-ESR in rheumatoid arthritis.

OBJECTIVES: To evaluate determinants of discordance between DAS28-ESR and DAS28-CRP and resulting impact on disease activity stratification in RA. METHODS: Paired DAS28-ESR and DAS28-CRP readings (n = 31 074) were obtained from the British Society for Rheumatology Biologics Register for RA. Factors influencing discordance between DAS28-ESR and DAS28-CRP were evaluated alongside the resulting effect on disease activity stratification. The impact of gender adjustment to the DAS28-CRP was evaluated. RESULTS: DAS28-CRP scores were ∼0.3 lower than DAS28-ESR overall, with greatest differences for women (-0.35) and patients over 50 years old (-0.34). Mean male DAS28-CRP scores were 0.15 less than corresponding DAS28-ESR scores. Discordance between DAS28-ESR and DAS28-CRP significantly impacted disease activity stratification at low disease activity and remission thresholds (32.0% and 66.6% concordance, respectively). Adjusting DAS28-CRP scores by gender significantly (P < 0.001) improved agreement with the DAS28-ESR. CONCLUSION: Discordance between DAS28-ESR and DAS28-CRP is greatest for women and patients over 50 years of age, and influences disease activity stratification. The proposed gender-adjusted DAS28-CRP improves inter-score agreement with DAS28-ESR, supporting more reliable disease activity stratification in treat-to-target approaches for RA.

Risk of type 2 diabetes and cardiovascular disease in an incident cohort of people with psoriatic arthritis: a population-based cohort study.

Objectives: To determine the risk of type 2 diabetes (T2D) and cardiovascular diseases in PsA patients compared with the general population and patients with psoriasis. Methods: Incident PsA patients aged 18-89 years were identified in the UK Clinical Practice Research Datalink between 1998 and 2014 and were matched (1:4 ratio) to a general population cohort and psoriasis cohort. The incidence of T2D, cerebrovascular disease, ischaemic heart disease and peripheral vascular disease (PVD) was calculated for each study cohort. Conditional Poisson regression was used to calculate adjusted relative risks. Results: We identified 6783 incident cases of PsA. The risk of T2D was significantly higher in the PsA cohort than in the general population and the psoriasis cohorts [adjusted relative risk 1.40 (CI95 1.15, 1.70) and adjusted relative risk 1.53 (CI95 1.19, 1.97), respectively]. The incidence of ischaemic heart disease, peripheral vascular disease and the three cardiovascular outcomes combined in the PsA cohort was significantly higher than in the general population. No significant differences in risk were observed between the PsA and psoriasis cohorts for any cardiovascular outcome. Conclusion: The development of T2D in an incident population of PsA is significantly higher than in psoriasis alone or in a general population, whereas the increased risk of cardiovascular disease in PsA and psoriasis is similar.

Risk of uveitis and inflammatory bowel disease in people with psoriatic arthritis: a population-based cohort study.

OBJECTIVES: To determine the risk of uveitis and inflammatory bowel disease (IBD) in patients with psoriatic arthritis (PsA) compared with the general population and patients with psoriasis. METHODS: A cohort study using data from the UK Clinical Practice Research Datalink between 1998 and 2014. Patients with incident PsA aged 18-89 years were identified and matched to a cohort of patients with psoriasis and a general population cohort. The incidence of uveitis, all IBD, Crohn's disease and ulcerative colitis was calculated for each study cohort and adjusted relative risks (RRadj) were calculated using conditional Poisson regression. RESULTS: 6783 incident cases of PsA were identified with a median age of 49 years. The risk of uveitis was significantly higher in the PsA cohort than in the general population and psoriasis cohorts (RRadj 3.55, 95% CI 2.21 to 5.70 and RRadj 2.13, 95% CI 1.40 to 3.24, respectively). A significant increase was observed for Crohn's disease (RRadj 2.96, 95% CI 1.46 to 6.00 and RRadj3.60, 95% CI 1.83 to 7.10) but not for ulcerative colitis (RRadj1.30, 95% CI 0.66 to 2.56 and RRadj0.98, 95% CI 0.50 to 1.92). CONCLUSIONS: In a primary care-based incidence cohort of patients with PsA, there were substantial risks of developing uveitis and/or Crohn's disease, but not ulcerative colitis, when compared with the general population and psoriasis controls.

Data Integration for the Assessment of Population Exposure to Ambient Air Pollution for Global Burden of Disease Assessment.

Air pollution is a leading global disease risk factor. Tracking progress (e.g., for Sustainable Development Goals) requires accurate, spatially resolved, routinely updated exposure estimates. A Bayesian hierarchical model was developed to estimate annual average fine particle (PM2.5) concentrations at 0.1° × 0.1° spatial resolution globally for 2010-2016. The model incorporated spatially varying relationships between 6003 ground measurements from 117 countries, satellite-based estimates, and other predictors. Model coefficients indicated larger contributions from satellite-based estimates in countries with low monitor density. Within and out-of-sample cross-validation indicated improved predictions of ground measurements compared to previous (Global Burden of Disease 2013) estimates (increased within-sample R2 from 0.64 to 0.91, reduced out-of-sample, global population-weighted root mean squared error from 23 µg/m3 to 12 µg/m3). In 2016, 95% of the world's population lived in areas where ambient PM2.5 levels exceeded the World Health Organization 10 µg/m3 (annual average) guideline; 58% resided in areas above the 35 µg/m3 Interim Target-1. Global population-weighted PM2.5 concentrations were 18% higher in 2016 (51.1 µg/m3) than in 2010 (43.2 µg/m3), reflecting in particular increases in populous South Asian countries and from Saharan dust transported to West Africa. Concentrations in China were high (2016 population-weighted mean: 56.4 µg/m3) but stable during this period.

Recovery of shape and size in a developing organ pair.

BACKGROUND: Paired organs in animals are largely bilaterally symmetric despite inherent noise in most biological processes. How is precise organ shape and size achieved during development despite this noise? Examining paired organ development is a challenge because it requires repeated quantification of two structures in parallel within living embryos. Here we combine bilateral quantification of morphology through time with asymmetric perturbations to study regulation of organ shape, size, and symmetry in developing organ pairs. RESULTS: We present quantitative live imaging tools to measure the shape and size of the developing inner ears on both the left and right side simultaneously over time. By quantifying variation between the left and right inner ear (intrinsic noise) and between different individuals (extrinsic noise), we find that initial variability decreases over time in normal development to achieve symmetry. Early asymmetry is increased by environmental stress, but symmetry is still recovered over subsequent developmental time. Using multiple unilateral perturbations including Fgf signaling and ultraviolet light, we find that growth can be adjusted to compensate for a range of initial size and shape differences. CONCLUSIONS: We propose that symmetry in developmental systems does not emerge through precise deterministic bilateral development, but rather through feedback mechanisms that adjust morphogenesis rates to account for variation. Developmental Dynamics 246:451-465, 2016. © 2017 Wiley Periodicals, Inc.

Interval between onset of psoriasis and psoriatic arthritis comparing the UK Clinical Practice Research Datalink with a hospital-based cohort.

Objectives: To describe the time interval between the onset of psoriasis and PsA in the UK primary care setting and compare with a large, well-classified secondary care cohort. Methods: Patients with PsA and/or psoriasis were identified in the UK Clinical Practice Research Datalink (CPRD). The secondary care cohort comprised patients from the Bath PsA longitudinal observational cohort study. For incident PsA patients in the CPRD who also had a record of psoriasis, the time interval between PsA diagnosis and first psoriasis record was calculated. Comparisons were made with the time interval between diagnoses in the Bath cohort. Results: There were 5272 eligible PsA patients in the CPRD and 815 in the Bath cohort. In both cohorts, the majority of patients (82.3 and 61.3%, respectively) had psoriasis before their PsA diagnosis or within the same calendar year (10.5 and 23.8%), with only a minority receiving their PsA diagnosis first (7.1 and 14.8%). Excluding those who presented with arthritis before psoriasis, the median time between diagnoses was 8 years [interquartile range (IQR) 2-15] in the CPRD and 7 years (IQR 0-20) in the Bath cohort. In the CPRD, 60.1 and 75.1% received their PsA diagnosis within 10 and 15 years of their psoriasis diagnosis, respectively; this was comparable with 57.2 and 67.7% in the Bath cohort. Conclusion: A similar distribution for the time interval between psoriasis and arthritis was observed in the CPRD and secondary care cohort. These data can inform screening strategies and support the validity of data from each cohort.

The safety of antivirals and neutralising monoclonal antibodies used in prehospital treatment of Covid-19.

OBJECTIVE: This proof-of-principle pharmacovigilance study used Electronic Health Record (EHR) data to examine the safety of sotrovimab, paxlovid and molnupiravir in prehospital treatment of Covid-19. METHOD: With NHS England approval, we conducted an observational cohort study using OpenSAFELY-TPP, a secure software-platform which executes analyses across EHRs for 24 million people in England. High-risk individuals with Covid-19 eligible for prehospital treatment were included. Adverse events (AEs) were categorised into events in the drug's Summary of Product Characteristics (SmPC), drug-reactions and immune-mediated. Cox models compared risk across treatments. A pre-pandemic record analysis was performed for comparative purposes. RESULTS: Between 2021-2023, 37,449 patients received sotrovimab, paxlovid or molnupiravir whilst 109,647 patients made up an eligible-but-untreated population. The 28-day rates of AEs were low: SmPC 0.34 per 1000 patient-years (95% CI 0.32-0.36); drug-reactions 0.01 (95% CI 0.01-0.02) and immune-mediated 0.03 (95% CI 0.03-0.04), and similar or lower than the pre-pandemic period. Compared with the eligible but untreated population, sotrovimab and paxlovid associated with a risk of SmPC AE [adjHR 1.36 (95% CI 1.15-1.62) and 1.28 (95% CI 1.05-1.55), respectively], whilst sotrovimab associated with a risk of drug-reactions [adjHR 2.95 (95% CI 1.56-5.55)] and immune-mediated events [adjHR 3.22 (95% CI 1.86-5.57)]. CONCLUSION: Sotrovimab, paxlovid and molnupiravir demonstrate acceptable safety profiles. Although the risk of AEs was greatest with sotrovimab, event rates were lower than comparative pre-pandemic period.

Impact of COVID-19 on recorded blood pressure screening and hypertension management in England: an analysis of monthly changes in the quality and outcomes framework indicators in OpenSAFELY.

BACKGROUND: The COVID-19 pandemic disrupted cardiovascular disease management in primary care in England. OBJECTIVE: To describe the impact of the pandemic on blood pressure screening and hypertension management based on a national quality of care scheme (Quality and Outcomes Framework, QOF) across key demographic, regional and clinical subgroups. METHODS: With NHS England approval, a population-based cohort study was conducted using OpenSAFELY-TPP on 25.2 million NHS patients registered at general practices (March 2019 to March 2023). We examined monthly changes in recorded blood pressure screening in the preceding 5 years in patients aged ≥45 years and recorded the hypertension prevalence and the percentage of patients treated to target (≤140/90 mmHg for patients aged ≤79 years and ≤150/90 mmHg for patients aged ≥80 years) in the preceding 12 months. RESULTS: The percentage of patients aged ≥45 years who had blood pressure screening recorded in the preceding 5 years decreased from 90% (March 2019) to 85% (March 2023). Recorded hypertension prevalence was relatively stable at 15% throughout the study period. The percentage of patients with a record of hypertension treated to target in the preceding 12 months reduced from a maximum of 71% (March 2020) to a minimum of 47% (February 2021) in patients aged ≤79 years and from 85% (March 2020) to a minimum of 58% (February 2021) in patients aged ≥80 years before recovery. Blood pressure screening rates in the preceding 5 years remained stable in older people, patients with recorded learning disability or care home status. CONCLUSIONS: The pandemic substantially disrupted hypertension management QOF indicators, which is likely attributable to general reductions of blood pressure measurement including screening. OpenSAFELY can be used to continuously monitor changes in national quality-of-care schemes to identify changes in key clinical subgroups early and support prioritisation of recovery from care disrupted by COVID-19.

Genetic control of organ shape and tissue polarity.

The mechanisms by which genes control organ shape are poorly understood. In principle, genes may control shape by modifying local rates and/or orientations of deformation. Distinguishing between these possibilities has been difficult because of interactions between patterns, orientations, and mechanical constraints during growth. Here we show how a combination of growth analysis, molecular genetics, and modelling can be used to dissect the factors contributing to shape. Using the Snapdragon (Antirrhinum) flower as an example, we show how shape development reflects local rates and orientations of tissue growth that vary spatially and temporally to form a dynamic growth field. This growth field is under the control of several dorsoventral genes that influence flower shape. The action of these genes can be modelled by assuming they modulate specified growth rates parallel or perpendicular to local orientations, established by a few key organisers of tissue polarity. Models in which dorsoventral genes only influence specified growth rates do not fully account for the observed growth fields and shapes. However, the data can be readily explained by a model in which dorsoventral genes also modify organisers of tissue polarity. In particular, genetic control of tissue polarity organisers at ventral petal junctions and distal boundaries allows both the shape and growth field of the flower to be accounted for in wild type and mutants. The results suggest that genetic control of tissue polarity organisers has played a key role in the development and evolution of shape.

Quantitative control of organ shape by combinatorial gene activity.

The development of organs with particular shapes, like wings or flowers, depends on regional activity of transcription factors and signalling molecules. However, the mechanisms that link these molecular activities to the morphogenetic events underlying shape are poorly understood. Here we describe a combination of experimental and computational approaches that address this problem, applying them to a group of genes controlling flower shape in the Snapdragon (Antirrhinum). Four transcription factors are known to play a key role in the control of floral shape and asymmetry in Snapdragon. We use quantitative shape analysis of mutants for these factors to define principal components underlying flower shape variation. We show that each transcription factor has a specific effect on the shape and size of regions within the flower, shifting the position of the flower in shape space. These shifts are further analysed by generating double mutants and lines that express some of the genes ectopically. By integrating these observations with known gene expression patterns and interactions, we arrive at a combinatorial scheme for how regional effects on shape are genetically controlled. We evaluate our scheme by incorporating the proposed interactions into a generative model, where the developing flower is treated as a material sheet that grows according to how genes modify local polarities and growth rates. The petal shapes generated by the model show a good quantitative match with those observed experimentally for each petal in numerous genotypes, thus validating the hypothesised scheme. This article therefore shows how complex shapes can be accounted for by combinatorial effects of transcription factors on regional growth properties. This finding has implications not only for how shapes develop but also for how they may have evolved through tinkering with transcription factors and their targets.

Factors associated with COVID-19 vaccine uptake in people with kidney disease: an OpenSAFELY cohort study.

OBJECTIVE: To characterise factors associated with COVID-19 vaccine uptake among people with kidney disease in England. DESIGN: Retrospective cohort study using the OpenSAFELY-TPP platform, performed with the approval of NHS England. SETTING: Individual-level routine clinical data from 24 million people across GPs in England using TPP software. Primary care data were linked directly with COVID-19 vaccine records up to 31 August 2022 and with renal replacement therapy (RRT) status via the UK Renal Registry (UKRR). PARTICIPANTS: A cohort of adults with stage 3-5 chronic kidney disease (CKD) or receiving RRT at the start of the COVID-19 vaccine roll-out was identified based on evidence of reduced estimated glomerular filtration rate (eGFR) or inclusion in the UKRR. MAIN OUTCOME MEASURES: Dose-specific vaccine coverage over time was determined from 1 December 2020 to 31 August 2022. Individual-level factors associated with receipt of a 3-dose or 4-dose vaccine series were explored via Cox proportional hazards models. RESULTS: 992 205 people with stage 3-5 CKD or receiving RRT were included. Cumulative vaccine coverage as of 31 August 2022 was 97.5%, 97.0% and 93.9% for doses 1, 2 and 3, respectively, and 81.9% for dose 4 among individuals with one or more indications for eligibility. Delayed 3-dose vaccine uptake was associated with younger age, minority ethnicity, social deprivation and severe mental illness-associations that were consistent across CKD severity subgroups, dialysis patients and kidney transplant recipients. Similar associations were observed for 4-dose uptake. CONCLUSION: Although high primary vaccine and booster dose coverage has been achieved among people with kidney disease in England, key disparities in vaccine uptake remain across clinical and demographic groups and 4-dose coverage is suboptimal. Targeted interventions are needed to identify barriers to vaccine uptake among under-vaccinated subgroups identified in the present study.

Comparative effectiveness of nirmatrelvir/ritonavir versus sotrovimab and molnupiravir for preventing severe COVID-19 outcomes in non-hospitalised high-risk patients during Omicron waves: observational cohort study using the OpenSAFELY platform.

Background: Timely evidence of the comparative effectiveness between COVID-19 therapies in real-world settings is needed to inform clinical care. This study aimed to compare the effectiveness of nirmatrelvir/ritonavir versus sotrovimab and molnupiravir in preventing severe COVID-19 outcomes in non-hospitalised high-risk COVID-19 adult patients during Omicron waves. Methods: With the approval of NHS England, we conducted a real-world cohort study using the OpenSAFELY-TPP platform. Patient-level primary care data were obtained from 24 million people in England and were securely linked with data on COVID-19 infection and therapeutics, hospital admission, and death, covering a period where both nirmatrelvir/ritonavir and sotrovimab were first-line treatment options in community settings (February 10, 2022-November 27, 2022). Molnupiravir (third-line option) was used as an exploratory comparator to nirmatrelvir/ritonavir, both of which were antivirals. Cox proportional hazards model stratified by area was used to compare the risk of 28-day COVID-19 related hospitalisation/death across treatment groups. Findings: A total of 9026 eligible patients treated with nirmatrelvir/ritonavir (n = 5704) and sotrovimab (n = 3322) were included in the main analysis. The mean age was 52.7 (SD = 14.9) years and 93% (8436/9026) had three or more COVID-19 vaccinations. Within 28 days after treatment initiation, 55/9026 (0.61%) COVID-19 related hospitalisations/deaths were observed (34/5704 [0.60%] treated with nirmatrelvir/ritonavir and 21/3322 [0.63%] with sotrovimab). After adjusting for demographics, high-risk cohort categories, vaccination status, calendar time, body mass index and other comorbidities, we observed no significant difference in outcome risk between nirmatrelvir/ritonavir and sotrovimab users (HR = 0.89, 95% CI: 0.48-1.63; P = 0.698). Results from propensity score weighted model also showed non-significant difference between treatment groups (HR = 0.82, 95% CI: 0.45-1.52; P = 0.535). The exploratory analysis comparing nirmatrelvir/ritonavir users with 1041 molnupiravir users (13/1041 [1.25%] COVID-19 related hospitalisations/deaths) showed an association in favour of nirmatrelvir/ritonavir (HR = 0.45, 95% CI: 0.22-0.94; P = 0.033). Interpretation: In routine care of non-hospitalised high-risk adult patients with COVID-19 in England, no substantial difference in the risk of severe COVID-19 outcomes was observed between those who received nirmatrelvir/ritonavir and sotrovimab between February and November 2022, when Omicron subvariants BA.2, BA.5, or BQ.1 were dominant. Funding: UK Research and Innovation, Wellcome Trust, UK Medical Research Council, National Institute for Health and Care Research, and Health Data Research UK.

Comparative effectiveness of sotrovimab and molnupiravir for preventing severe COVID-19 outcomes in patients on kidney replacement therapy: observational study using the OpenSAFELY-UKRR and SRR databases.

Background: Due to limited inclusion of patients on kidney replacement therapy (KRT) in clinical trials, the effectiveness of coronavirus disease 2019 (COVID-19) therapies in this population remains unclear. We sought to address this by comparing the effectiveness of sotrovimab against molnupiravir, two commonly used treatments for non-hospitalised KRT patients with COVID-19 in the UK. Methods: With the approval of National Health Service England, we used routine clinical data from 24 million patients in England within the OpenSAFELY-TPP platform linked to the UK Renal Registry (UKRR) to identify patients on KRT. A Cox proportional hazards model was used to estimate hazard ratios (HRs) of sotrovimab versus molnupiravir with regards to COVID-19-related hospitalisations or deaths in the subsequent 28 days. We also conducted a complementary analysis using data from the Scottish Renal Registry (SRR). Results: Among the 2367 kidney patients treated with sotrovimab (n = 1852) or molnupiravir (n = 515) between 16 December 2021 and 1 August 2022 in England, 38 cases (1.6%) of COVID-19-related hospitalisations/deaths were observed. Sotrovimab was associated with substantially lower outcome risk than molnupiravir {adjusted HR 0.35 [95% confidence interval (CI) 0.17-0.71]; P = .004}, with results remaining robust in multiple sensitivity analyses. In the SRR cohort, sotrovimab showed a trend toward lower outcome risk than molnupiravir [HR 0.39 (95% CI 0.13-1.21); P = .106]. In both datasets, sotrovimab had no evidence of an association with other hospitalisation/death compared with molnupiravir (HRs ranged from 0.73 to 1.29; P > .05). Conclusions: In routine care of non-hospitalised patients with COVID-19 on KRT, sotrovimab was associated with a lower risk of severe COVID-19 outcomes compared with molnupiravir during Omicron waves.

Changes in medication safety indicators in England throughout the covid-19 pandemic using OpenSAFELY: population based, retrospective cohort study of 57 million patients using federated analytics.

Objective: To implement complex, PINCER (pharmacist led information technology intervention) prescribing indicators, on a national scale with general practice data to describe the impact of the covid-19 pandemic on safe prescribing. Design: Population based, retrospective cohort study using federated analytics. Setting: Electronic general practice health record data from 56.8 million NHS patients by use of the OpenSAFELY platform, with the approval of the National Health Service (NHS) England. Participants: NHS patients (aged 18-120 years) who were alive and registered at a general practice that used TPP or EMIS computer systems and were recorded as at risk of at least one potentially hazardous PINCER indicator. Main outcome measure: Between 1 September 2019 and 1 September 2021, monthly trends and between practice variation for compliance with 13 PINCER indicators, as calculated on the first of every month, were reported. Prescriptions that do not adhere to these indicators are potentially hazardous and can cause gastrointestinal bleeds; are cautioned against in specific conditions (specifically heart failure, asthma, and chronic renal failure); or require blood test monitoring. The percentage for each indicator is formed of a numerator of patients deemed to be at risk of a potentially hazardous prescribing event and the denominator is of patients for which assessment of the indicator is clinically meaningful. Higher indicator percentages represent potentially poorer performance on medication safety. Results: The PINCER indicators were successfully implemented across general practice data for 56.8 million patient records from 6367 practices in OpenSAFELY. Hazardous prescribing remained largely unchanged during the covid-19 pandemic, with no evidence of increases in indicators of harm as captured by the PINCER indicators. The percentage of patients at risk of potentially hazardous prescribing, as defined by each PINCER indicator, at mean quarter 1 (Q1) 2020 (representing before the pandemic) ranged from 1.11% (age ≥65 years and non-steroidal anti-inflammatory drugs) to 36.20% (amiodarone and no thyroid function test), while Q1 2021 (representing after the pandemic) percentages ranged from 0.75% (age ≥65 years and non-steroidal anti-inflammatory drugs) to 39.23% (amiodarone and no thyroid function test). Transient delays occurred in blood test monitoring for some medications, particularly angiotensin-converting enzyme inhibitors (where blood monitoring worsened from a mean of 5.16% in Q1 2020 to 12.14% in Q1 2021, and began to recover in June 2021). All indicators substantially recovered by September 2021. We identified 1 813 058 patients (3.1%) at risk of at least one potentially hazardous prescribing event. Conclusion: NHS data from general practices can be analysed at national scale to generate insights into service delivery. Potentially hazardous prescribing was largely unaffected by the covid-19 pandemic in primary care health records in England.

Impact of the COVID-19 pandemic on antipsychotic prescribing in individuals with autism, dementia, learning disability, serious mental illness or living in a care home: a federated analysis of 59 million patients' primary care records in situ using OpenSAFELY.

BACKGROUND: The COVID-19 pandemic affected how care was delivered to vulnerable patients, such as those with dementia or learning disability. OBJECTIVE: To explore whether this affected antipsychotic prescribing in at-risk populations. METHODS: With the approval of NHS England, we completed a retrospective cohort study, using the OpenSAFELY platform to explore primary care data of 59 million patients. We identified patients in five at-risk groups: autism, dementia, learning disability, serious mental illness and care home residents. We calculated the monthly prevalence of antipsychotic prescribing in these groups, as well as the incidence of new prescriptions in each month. FINDINGS: The average monthly rate of antipsychotic prescribing increased in dementia from 82.75 patients prescribed an antipsychotic per 1000 patients (95% CI 82.30 to 83.19) in January-March 2019 to 90.1 (95% CI 89.68 to 90.60) in October-December 2021 and from 154.61 (95% CI 153.79 to 155.43) to 166.95 (95% CI 166.23 to 167.67) in care homes. There were notable spikes in the rate of new prescriptions issued to patients with dementia and in care homes. In learning disability and autism groups, the rate of prescribing per 1000 decreased from 122.97 (95% CI 122.29 to 123.66) to 119.29 (95% CI 118.68 to 119.91) and from 54.91 (95% CI 54.52 to 55.29) to 51.04 (95% CI 50.74 to 51.35), respectively. CONCLUSION AND IMPLICATIONS: We observed a spike in antipsychotic prescribing in the dementia and care home groups, which correlated with lockdowns and was likely due to prescribing of antipsychotics for palliative care. We observed gradual increases in antipsychotic use in dementia and care home patients and decreases in their use in patients with learning disability or autism.

Trends, variation, and clinical characteristics of recipients of antiviral drugs and neutralising monoclonal antibodies for covid-19 in community settings: retrospective, descriptive cohort study of 23.4 million people in OpenSAFELY.

Objective: To ascertain patient eligibility status and describe coverage of antiviral drugs and neutralising monoclonal antibodies (nMAB) as treatment for covid-19 in community settings in England. Design: Retrospective, descriptive cohort study, approved by NHS England. Setting: Routine clinical data from 23.4 million people linked to data on covid-19 infection and treatment, within the OpenSAFELY-TPP database. Participants: Outpatients with covid-19 at high risk of severe outcomes. Interventions: Nirmatrelvir/ritonavir (paxlovid), sotrovimab, molnupiravir, casirivimab/imdevimab, or remdesivir, used in the community by covid-19 medicine delivery units. Results: 93 870 outpatients with covid-19 were identified between 11 December 2021 and 28 April 2022 to be at high risk of severe outcomes and therefore potentially eligible for antiviral or nMAB treatment (or both). Of these patients, 19 040 (20%) received treatment (sotrovimab, 9660 (51%); molnupiravir, 4620 (24%); paxlovid, 4680 (25%); casirivimab/imdevimab, 50 (<1%); and remdesivir, 30 (<1%)). The proportion of patients treated increased from 9% (190/2220) in the first week of treatment availability to 29% (460/1600) in the latest week. The proportion treated varied by high risk group, being lowest in those with liver disease (16%; 95% confidence interval 15% to 17%); by treatment type, with sotrovimab favoured over molnupiravir and paxlovid in all but three high risk groups (Down's syndrome (35%; 30% to 39%), rare neurological conditions (45%; 43% to 47%), and immune deficiencies (48%; 47% to 50%)); by age, ranging from ≥80 years (13%; 12% to 14%) to 50-59 years (23%; 22% to 23%); by ethnic group, ranging from black (11%; 10% to 12%) to white (21%; 21% to 21%); by NHS region, ranging from 13% (12% to 14%) in Yorkshire and the Humber to 25% (24% to 25%) in the East of England); and by deprivation level, ranging from 15% (14% to 15%) in the most deprived areas to 23% (23% to 24%) in the least deprived areas. Groups that also had lower coverage included unvaccinated patients (7%; 6% to 9%), those with dementia (6%; 5% to 7%), and care home residents (6%; 6% to 7%). Conclusions: Using the OpenSAFELY platform, we were able to identify patients with covid-19 at high risk of severe outcomes who were potentially eligible to receive treatment and assess the coverage of these new treatments among these patients. In the context of a rapid deployment of a new service, the NHS analytical code used to determine eligibility could have been over-inclusive and some of the eligibility criteria not fully captured in healthcare data. However targeted activity might be needed to resolve apparent lower treatment coverage observed among certain groups, in particular (at present): different NHS regions, ethnic groups, people aged ≥80 years, those living in socioeconomically deprived areas, and care home residents.