Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7.

SARS-CoV-2 lineage B.1.1.7, a variant that was first detected in the UK in September 20201, has spread to multiple countries worldwide. Several studies have established that B.1.1.7 is more transmissible than pre-existing variants, but have not identified whether it leads to any change in disease severity2. Here we analyse a dataset that links 2,245,263 positive SARS-CoV-2 community tests and 17,452 deaths associated with COVID-19 in England from 1 November 2020 to 14 February 2021. For 1,146,534 (51%) of these tests, the presence or absence of B.1.1.7 can be identified because mutations in this lineage prevent PCR amplification of the spike (S) gene target (known as S gene target failure (SGTF)1). On the basis of 4,945 deaths with known SGTF status, we estimate that the hazard of death associated with SGTF is 55% (95% confidence interval, 39-72%) higher than in cases without SGTF after adjustment for age, sex, ethnicity, deprivation, residence in a care home, the local authority of residence and test date. This corresponds to the absolute risk of death for a 55-69-year-old man increasing from 0.6% to 0.9% (95% confidence interval, 0.8-1.0%) within 28 days of a positive test in the community. Correcting for misclassification of SGTF and missingness in SGTF status, we estimate that the hazard of death associated with B.1.1.7 is 61% (42-82%) higher than with pre-existing variants. Our analysis suggests that B.1.1.7 is not only more transmissible than pre-existing SARS-CoV-2 variants, but may also cause more severe illness.

Prediction Under Interventions: Evaluation of Counterfactual Performance Using Longitudinal Observational Data.

Predictions under interventions are estimates of what a person's risk of an outcome would be if they were to follow a particular treatment strategy, given their individual characteristics. Such predictions can give important input to medical decision-making. However, evaluating the predictive performance of interventional predictions is challenging. Standard ways of evaluating predictive performance do not apply when using observational data, because prediction under interventions involves obtaining predictions of the outcome under conditions that are different from those that are observed for a subset of individuals in the validation dataset. This work describes methods for evaluating counterfactual performance of predictions under interventions for time-to-event outcomes. This means we aim to assess how well predictions would match the validation data if all individuals had followed the treatment strategy under which predictions are made. We focus on counterfactual performance evaluation using longitudinal observational data, and under treatment strategies that involve sustaining a particular treatment regime over time. We introduce an estimation approach using artificial censoring and inverse probability weighting that involves creating a validation dataset mimicking the treatment strategy under which predictions are made. We extend measures of calibration, discrimination (c-index and cumulative/dynamic AUCt) and overall prediction error (Brier score) to allow assessment of counterfactual performance. The methods are evaluated using a simulation study, including scenarios in which the methods should detect poor performance. Applying our methods in the context of liver transplantation shows that our procedure allows quantification of the performance of predictions supporting crucial decisions on organ allocation.

Effectiveness of mRNA COVID-19 Vaccines as First Booster Doses in England: An Observational Study in OpenSAFELY-TPP.

BACKGROUND: The UK delivered its first "booster" COVID-19 vaccine doses in September 2021, initially to individuals at high risk of severe disease, then to all adults. The BNT162b2 Pfizer-BioNTech vaccine was used initially, then also Moderna mRNA-1273. METHODS: With the approval of the National Health Service England, we used routine clinical data to estimate the effectiveness of boosting with BNT162b2 or mRNA-1273 compared with no boosting in eligible adults who had received two primary course vaccine doses. We matched each booster recipient with an unboosted control on factors relating to booster priority status and prior COVID-19 immunization. We adjusted for additional factors in Cox models, estimating hazard ratios up to 182 days (6 months) following booster dose. We estimated hazard ratios overall and within the following periods: 1-14, 15-42, 43-69, 70-97, 98-126, 127-152, and 155-182 days. Outcomes included a positive SARS-CoV-2 test, COVID-19 hospitalization, COVID-19 death, non-COVID-19 death, and fracture. RESULTS: We matched 8,198,643 booster recipients with unboosted controls. Adjusted hazard ratios over 6-month follow-up were: positive SARS-CoV-2 test 0.75 (0.74, 0.75); COVID-19 hospitalization 0.30 (0.29, 0.31); COVID-19 death 0.11 (0.10, 0.14); non-COVID-19 death 0.22 (0.21, 0.23); and fracture 0.77 (0.75, 0.78). Estimated effectiveness of booster vaccines against severe COVID-19-related outcomes peaked during the first 3 months following the booster dose. By 6 months, the cumulative incidence of positive SARS-CoV-2 test was higher in boosted than unboosted individuals. CONCLUSIONS: We estimate that COVID-19 booster vaccination, compared with no booster vaccination, provided substantial protection against COVID-19 hospitalization and COVID-19 death but only limited protection against positive SARS-CoV-2 test. Lower rates of fracture in boosted than unboosted individuals may suggest unmeasured confounding. Observational studies should report estimated vaccine effectiveness against nontarget and negative control outcomes.

Diagnosis of human leptospirosis: systematic review and meta-analysis of the diagnostic accuracy of the Leptospira microscopic agglutination test, PCR targeting Lfb1, and IgM ELISA to Leptospira fainei serovar Hurstbridge.

BACKGROUND: Leptospirosis is an underdiagnosed infectious disease with non-specific clinical presentation that requires laboratory confirmation for diagnosis. The serologic reference standard remains the microscopic agglutination test (MAT) on paired serum samples. However, reported estimates of MAT's sensitivity vary. We evaluated the accuracy of four index tests, MAT on paired samples as well as alternative standards for leptospirosis diagnosis: MAT on single acute-phase samples, polymerase chain reaction (PCR) with the target gene Lfb1, and ELISA IgM with Leptospira fainei serovar Hurstbridge as an antigen. METHODS: We performed a systematic review of studies reporting results of leptospirosis diagnostic tests. We searched eight electronic databases and selected studies that tested human blood samples and compared index tests with blood culture and/or PCR and/or MAT (comparator tests). For MAT selection criteria we defined a threshold for single acute-phase samples according to a national classification of leptospirosis endemicity. We used a Bayesian random-effect meta-analysis to estimate the sensitivity and specificity of MAT in single acute-phase and paired samples separately, and assessed risk of bias using the Quality Assessment of Studies of Diagnostic Accuracy Approach- 2 (QUADAS-2) tool. RESULTS: For the MAT accuracy evaluation, 15 studies were included, 11 with single acute-phase serum, and 12 with paired sera. Two included studies used PCR targeting the Lfb1 gene, and one included study used IgM ELISA with Leptospira fainei serovar Hurstbridge as antigen. For MAT in single acute-phase samples, the pooled sensitivity and specificity were 14% (95% credible interval [CrI] 3-38%) and 86% (95% CrI 59-96%), respectively, and the predicted sensitivity and specificity were 14% (95% CrI 0-90%) and 86% (95% CrI 9-100%). Among paired MAT samples, the pooled sensitivity and specificity were 68% (95% CrI 32-92%) and 75% (95% CrI 45-93%) respectively, and the predicted sensitivity and specificity were 69% (95% CrI 2-100%) and 75% (2-100%). CONCLUSIONS: Based on our analysis, the accuracy of MAT in paired samples was not high, but it remains the reference standard until a more accurate diagnostic test is developed. Future studies that include larger numbers of participants with paired samples will improve the certainty of accuracy estimates.

Measurement error models with zero inflation and multiple sources of zeros, with applications to hard zeros.

We consider measurement error models for two variables observed repeatedly and subject to measurement error. One variable is continuous, while the other variable is a mixture of continuous and zero measurements. This second variable has two sources of zeros. The first source is episodic zeros, wherein some of the measurements for an individual may be zero and others positive. The second source is hard zeros, i.e., some individuals will always report zero. An example is the consumption of alcohol from alcoholic beverages: some individuals consume alcoholic beverages episodically, while others never consume alcoholic beverages. However, with a small number of repeat measurements from individuals, it is not possible to determine those who are episodic zeros and those who are hard zeros. We develop a new measurement error model for this problem, and use Bayesian methods to fit it. Simulations and data analyses are used to illustrate our methods. Extensions to parametric models and survival analysis are discussed briefly.

Emulated trial investigating effects of multiple treatments: estimating combined effects of mucoactive nebulisers in cystic fibrosis using registry data.

INTRODUCTION: People with cystic fibrosis (CF) are often on multiple long-term treatments, including mucoactive nebulisers. In the UK, the most common mucoactive nebuliser is dornase alfa (DNase). A common therapeutic approach for people already on DNase is to add hypertonic saline (HS). The effects of DNase and HS used alone have been studied in randomised trials, but their effects in combination have not. This study investigates whether, for people already prescribed DNase, adding HS has additional benefit for lung function or use of intravenous antibiotics. METHODS: Using UK CF Registry data from 2007 to 2018, we emulated a target trial. We included people aged 6 years and over who were prescribed DNase without HS for 2 years. We investigated the effects of combinations of DNase and HS over 5 years of follow-up. Inverse-probability-of-treatment weighting was used to control confounding. The period predated triple combination CF transmembrane conductance regulator modulators in routine care. RESULTS: 4498 individuals were included. At baseline, average age and forced expiratory volume in 1 s (FEV1%) predicted were 21.1 years and 69.7 respectively. During first year of follow-up, 3799 individuals were prescribed DNase alone; 426 added HS; 57 switched to HS alone and 216 were prescribed neither. We found no evidence that adding HS improved FEV1% at 1-5 years, or use of intravenous antibiotics at 1-4 years, compared with DNase alone. CONCLUSION: For individuals with CF prescribed DNase, we found no evidence that adding HS had an effect on FEV1% or prescription of intravenous antibiotics. Our study illustrates the emulated target trial approach using CF Registry data.

Causal inference in survival analysis using longitudinal observational data: Sequential trials and marginal structural models.

Longitudinal observational data on patients can be used to investigate causal effects of time-varying treatments on time-to-event outcomes. Several methods have been developed for estimating such effects by controlling for the time-dependent confounding that typically occurs. The most commonly used is marginal structural models (MSM) estimated using inverse probability of treatment weights (IPTW) (MSM-IPTW). An alternative, the sequential trials approach, is increasingly popular, and involves creating a sequence of "trials" from new time origins and comparing treatment initiators and non-initiators. Individuals are censored when they deviate from their treatment assignment at the start of each "trial" (initiator or noninitiator), which is accounted for using inverse probability of censoring weights. The analysis uses data combined across trials. We show that the sequential trials approach can estimate the parameters of a particular MSM. The causal estimand that we focus on is the marginal risk difference between the sustained treatment strategies of "always treat" vs "never treat." We compare how the sequential trials approach and MSM-IPTW estimate this estimand, and discuss their assumptions and how data are used differently. The performance of the two approaches is compared in a simulation study. The sequential trials approach, which tends to involve less extreme weights than MSM-IPTW, results in greater efficiency for estimating the marginal risk difference at most follow-up times, but this can, in certain scenarios, be reversed at later time points and relies on modelling assumptions. We apply the methods to longitudinal observational data from the UK Cystic Fibrosis Registry to estimate the effect of dornase alfa on survival.

Estimating sensitivity and specificity of diagnostic tests using latent class models that account for conditional dependence between tests: a simulation study.

BACKGROUND: Latent class models are increasingly used to estimate the sensitivity and specificity of diagnostic tests in the absence of a gold standard, and are commonly fitted using Bayesian methods. These models allow us to account for 'conditional dependence' between two or more diagnostic tests, meaning that the results from tests are correlated even after conditioning on the person's true disease status. The challenge is that it is not always clear to researchers whether conditional dependence exists between tests and whether it exists in all or just some latent classes. Despite the increasingly widespread use of latent class models to estimate diagnostic test accuracy, the impact of the conditional dependence structure chosen on the estimates of sensitivity and specificity remains poorly investigated. METHODS: A simulation study and a reanalysis of a published case study are used to highlight the impact of the conditional dependence structure chosen on estimates of sensitivity and specificity. We describe and implement three latent class random-effect models with differing conditional dependence structures, as well as a conditional independence model and a model that assumes perfect test accuracy. We assess the bias and coverage of each model in estimating sensitivity and specificity across different data generating mechanisms. RESULTS: The findings highlight that assuming conditional independence between tests within a latent class, where conditional dependence exists, results in biased estimates of sensitivity and specificity and poor coverage. The simulations also reiterate the substantial bias in estimates of sensitivity and specificity when incorrectly assuming a reference test is perfect. The motivating example of tests for Melioidosis highlights these biases in practice with important differences found in estimated test accuracy under different model choices. CONCLUSIONS: We have illustrated that misspecification of the conditional dependence structure leads to biased estimates of sensitivity and specificity when there is a correlation between tests. Due to the minimal loss in precision seen by using a more general model, we recommend accounting for conditional dependence even if researchers are unsure of its presence or it is only expected at minimal levels.

Diagnostic performance of the IMMY cryptococcal antigen lateral flow assay on serum and cerebrospinal fluid for diagnosis of cryptococcosis in HIV-negative patients: a systematic review.

BACKGROUND: The incidence of cryptococcosis amongst HIV-negative persons is increasing. Whilst the excellent performance of the CrAg testing in people living with HIV is well described, the diagnostic performance of the CrAg LFA has not been systematically evaluated in HIV-negative cohorts on serum or cerebrospinal fluid. METHODS: We performed a systematic review to characterise the diagnostic performance of IMMY CrAg® LFA in HIV-negative populations on serum and cerebrospinal fluid. A systematic electronic search was performed using Medline, Embase, Global Health, CENTRAL, WoS Science Citation Index, SCOPUS, Africa-Wide Information, LILACS and WHO Global Health Library. Studies were screened and data extracted from eligible studies by two independent reviewers. A fixed effect meta-analysis was used to estimate the diagnostic sensitivity and specificity. RESULTS: Of 447 records assessed for eligibility, nine studies met our inclusion criteria, including 528 participants overall. Amongst eight studies that evaluated the diagnostic performance of the IMMY CrAg® LFA on serum, the pooled median sensitivity was 96% (95% Credible Interval (CrI) 68-100%) with a pooled specificity estimate of 96% (95%CrI 84-100%). Amongst six studies which evaluated the diagnostic performance of IMMY CrAg® LFA on CSF, the pooled median sensitivity was 99% (95%CrI 95-100%) with a pooled specificity median of 99% (95%CrI 95-100%). CONCLUSIONS: This review demonstrates a high pooled sensitivity and specificity for the IMMY CrAg® LFA in HIV-negative populations, in keeping with findings in HIV-positive individuals. The review was limited by the small number of studies. Further studies using IMMY CrAg® LFA in HIV-negative populations would help to better determine the diagnostic value of this test.

Accuracy of the direct agglutination test for diagnosis of visceral leishmaniasis: a systematic review and meta-analysis.

BACKGROUND: Parasitological investigation of bone marrow, splenic or lymph node aspirations is the gold standard for the diagnosis of visceral leishmaniasis (VL). However, this invasive test requires skilled clinical and laboratory staff and adequate facilities, and sensitivity varies depending on the tissue used. The direct agglutination test (DAT) is a serological test that does not need specialised staff, with just minimal training required. While previous meta-analysis has shown DAT to have high sensitivity and specificity when using parasitology as the reference test for diagnosis, meta-analysis of DAT compared to other diagnostic techniques, such as PCR and ELISA, that are increasingly used in clinical and research settings, has not been done. METHODS: We conducted a systematic review to determine the diagnostic performance of DAT compared to all available tests for the laboratory diagnosis of human VL. We searched electronic databases including Medline, Embase, Global Health, Scopus, WoS Science Citation Index, Wiley Cochrane Central Register of Controlled Trials, Africa-Wide Information, LILACS and WHO Global Index. Three independent reviewers screened reports and extracted data from eligible studies. A meta-analysis estimated the diagnostic sensitivity and specificity of DAT. RESULTS: Of 987 titles screened, 358 were selected for full data extraction and 78 were included in the analysis, reporting on 32,822 participants from 19 countries. Studies included were conducted between 1987-2020. Meta-analysis of studies using serum and DAT compared to any other test showed pooled sensitivity of 95% (95%CrI 90-98%) and pooled specificity of 95% (95%CrI 88-98%). Results were similar for freeze-dried DAT and liquid DAT when analysed separately. Sensitivity was lower for HIV-positive patients (90%, CrI 59-98%) and specificity was lower for symptomatic patients (70%, CrI 43-89%). When comparing different geographical regions, the lowest median sensitivity (89%, CrI 67-97%) was in Western Asia (five studies). CONCLUSIONS: This systematic review and meta-analysis demonstrates high estimated pooled sensitivity and specificity of DAT for diagnosis of VL, although sensitivity and specificity were lower for different patient groups and geographical locations. This review highlights the lack of standardisation of DAT methods and preparations, and the lack of data from some important geographical locations. Future well-reported studies could provide better evidence to inform test implementation for different patient populations and use cases. PROSPERO REGISTRATION: CRD42021240830.

Extracorporeal Membrane Oxygenation for Severe Acute Respiratory Distress Syndrome Associated with COVID-19: An Emulated Target Trial Analysis.

Rationale: Whether patients with coronavirus disease (COVID-19) may benefit from extracorporeal membrane oxygenation (ECMO) compared with conventional invasive mechanical ventilation (IMV) remains unknown. Objectives: To estimate the effect of ECMO on 90-day mortality versus IMV only. Methods: Among 4,244 critically ill adult patients with COVID-19 included in a multicenter cohort study, we emulated a target trial comparing the treatment strategies of initiating ECMO versus no ECMO within 7 days of IMV in patients with severe acute respiratory distress syndrome (PaO2/FiO2 < 80 or PaCO2 ⩾ 60 mm Hg). We controlled for confounding using a multivariable Cox model on the basis of predefined variables. Measurements and Main Results: A total of 1,235 patients met the full eligibility criteria for the emulated trial, among whom 164 patients initiated ECMO. The ECMO strategy had a higher survival probability on Day 7 from the onset of eligibility criteria (87% vs. 83%; risk difference, 4%; 95% confidence interval, 0-9%), which decreased during follow-up (survival on Day 90: 63% vs. 65%; risk difference, -2%; 95% confidence interval, -10 to 5%). However, ECMO was associated with higher survival when performed in high-volume ECMO centers or in regions where a specific ECMO network organization was set up to handle high demand and when initiated within the first 4 days of IMV and in patients who are profoundly hypoxemic. Conclusions: In an emulated trial on the basis of a nationwide COVID-19 cohort, we found differential survival over time of an ECMO compared with a no-ECMO strategy. However, ECMO was consistently associated with better outcomes when performed in high-volume centers and regions with ECMO capacities specifically organized to handle high demand.

Estimating distribution of length of stay in a multi-state model conditional on the pathway, with an application to patients hospitalised with Covid-19.

Multi-state models are used to describe how individuals transition through different states over time. The distribution of the time spent in different states, referred to as 'length of stay', is often of interest. Methods for estimating expected length of stay in a given state are well established. The focus of this paper is on the distribution of the time spent in different states conditional on the complete pathway taken through the states, which we call 'conditional length of stay'. This work is motivated by questions about length of stay in hospital wards and intensive care units among patients hospitalised due to Covid-19. Conditional length of stay estimates are useful as a way of summarising individuals' transitions through the multi-state model, and also as inputs to mathematical models used in planning hospital capacity requirements. We describe non-parametric methods for estimating conditional length of stay distributions in a multi-state model in the presence of censoring, including conditional expected length of stay (CELOS). Methods are described for an illness-death model and then for the more complex motivating example. The methods are assessed using a simulation study and shown to give unbiased estimates of CELOS, whereas naive estimates of CELOS based on empirical averages are biased in the presence of censoring. The methods are applied to estimate conditional length of stay distributions for individuals hospitalised due to Covid-19 in the UK, using data on 42,980 individuals hospitalised from March to July 2020 from the COVID19 Clinical Information Network.

Using Negative Control Outcomes and Difference-in-Differences Analysis to Estimate Treatment Effects in an Entirely Treated Cohort: The Effect of Ivacaftor in Cystic Fibrosis.

When an entire cohort of patients receives a treatment, it is difficult to estimate the treatment effect in the treated because there are no directly comparable untreated patients. Attempts can be made to find a suitable control group (e.g., historical controls), but underlying differences between the treated and untreated can result in bias. Here we show how negative control outcomes combined with difference-in-differences analysis can be used to assess bias in treatment effect estimates and obtain unbiased estimates under certain assumptions. Causal diagrams and potential outcomes are used to explain the methods and assumptions. We apply the methods to UK Cystic Fibrosis Registry data to investigate the effect of ivacaftor, introduced in 2012 for a subset of the cystic fibrosis population with a particular genotype, on lung function and annual rate (days/year) of receiving intravenous (IV) antibiotics (i.e., IV days). We consider 2 negative control outcomes: outcomes measured in the pre-ivacaftor period and outcomes among persons ineligible for ivacaftor because of their genotype. Ivacaftor was found to improve lung function in year 1 (an approximately 6.5-percentage-point increase in ppFEV1), was associated with reduced lung function decline (an approximately 0.5-percentage-point decrease in annual ppFEV1 decline, though confidence intervals included 0), and reduced the annual rate of IV days (approximately 60% over 3 years).

Lung function in children with cystic fibrosis in the USA and UK: a comparative longitudinal analysis of national registry data.

RATIONALE: A previous analysis found significantly higher lung function in the US paediatric cystic fibrosis (CF) population compared with the UK with this difference apparently decreasing in adolescence and adulthood. However, the cross-sectional nature of the study makes it hard to interpret these results. OBJECTIVES: To compare longitudinal trajectories of lung function in children with CF between the USA and UK and to explore reasons for any differences. METHODS: We used mixed effects regression analysis to model lung function trajectories in the study populations. Using descriptive statistics, we compared early growth and nutrition (height, weight, body mass index), infections (Pseudomonas aeruginosa, Staphylococcus aureus) and treatments (rhDnase, hypertonic saline, inhaled antibiotics). RESULTS: We included 9463 children from the USA and 3055 children from the UK with homozygous F508del genotype. Lung function was higher in the USA than in the UK when first measured at age six and remained higher throughout childhood. We did not find important differences in early growth and nutrition, or P.aeruginosa infection. Prescription of rhDNase and hypertonic saline was more common in the USA. Inhaled antibiotics were prescribed at similar levels in both countries, but Tobramycin was prescribed more in the USA and colistin in the UK. S. aureus infection was more common in the USA than the UK. CONCLUSIONS: Children with CF and homozygous F508del genotype in the USA had better lung function than UK children. These differences do not appear to be explained by early growth or nutrition, but differences in the use of early treatments need further investigation.

Projecting the impact of triple CFTR modulator therapy on intravenous antibiotic requirements in cystic fibrosis using patient registry data combined with treatment effects from randomised trials.

BACKGROUND: Cystic fibrosis (CF) is a life-threatening genetic disease, affecting around 10 500 people in the UK. Precision medicines have been developed to treat specific CF-gene mutations. The newest, elexacaftor/tezacaftor/ivacaftor (ELEX/TEZ/IVA), has been found to be highly effective in randomised controlled trials (RCTs) and became available to a large proportion of UK CF patients in 2020. Understanding the potential health economic impacts of ELEX/TEZ/IVA is vital to planning service provision. METHODS: We combined observational UK CF Registry data with RCT results to project the impact of ELEX/TEZ/IVA on total days of intravenous (IV) antibiotic treatment at a population level. Registry data from 2015 to 2017 were used to develop prediction models for IV days over a 1-year period using several predictors, and to estimate 1-year population total IV days based on standards of care pre-ELEX/TEZ/IVA. We considered two approaches to imposing the impact of ELEX/TEZ/IVA on projected outcomes using effect estimates from RCTs: approach 1 based on effect estimates on FEV1% and approach 2 based on effect estimates on exacerbation rate. RESULTS: ELEX/TEZ/IVA is expected to result in significant reductions in population-level requirements for IV antibiotics of 16.1% (~17 800 days) using approach 1 and 43.6% (~39 500 days) using approach 2. The two approaches require different assumptions. Increased understanding of the mechanisms through which ELEX/TEZ/IVA acts on these outcomes would enable further refinements to our projections. CONCLUSIONS: This work contributes to increased understanding of the changing healthcare needs of people with CF and illustrates how Registry data can be used in combination with RCT evidence to estimate population-level treatment impacts.

Mediation of the total effect of cystic fibrosis-related diabetes on mortality: A UK Cystic Fibrosis Registry cohort study.

AIM: To investigate whether the effect of cystic fibrosis-related diabetes (CFRD) on the composite outcome of mortality or transplant could act through lung function, pulmonary exacerbations and/or nutritional status. METHODS: A retrospective cohort of adult cystic fibrosis (CF) patients who had not been diagnosed with CFRD were identified from the UK Cystic Fibrosis Registry (n = 2750). Rate of death or transplant was compared between patients who did and did not develop CFRD (with insulin use) during follow-up using Poisson regression, separately by sex. Causal mediation methods were used to investigate whether lung function, pulmonary exacerbations and nutritional status lie on the causal pathway between insulin-treated CFRD and mortality/transplant. RESULTS: At all ages, the mortality/transplant rate was higher in both men and women diagnosed with CFRD. Pulmonary exacerbations were the strongest mediator of the effect of CFRD on mortality/transplant, with an estimated 15% [95% CI: 7%, 28%] of the effect at 2 years post-CFRD diagnosis attributed to exacerbations, growing to 24% [95% CI: 9%, 46%] at 4 years post-diagnosis. Neither lung function nor nutritional status were found to be significant mediators of this effect. Estimates were similar but with wider confidence intervals in a cohort that additionally included people with CFRD but not using insulin. CONCLUSION: There is evidence that pulmonary exacerbations mediate the effect of CFRD on mortality but, as they are estimated to mediate less than one-quarter of the total effect, the mechanism through which CFRD influences survival may involve other factors.

Diagnostic accuracy of multiplex respiratory pathogen panels for influenza or respiratory syncytial virus infections: systematic review and meta-analysis.

Respiratory syncytial virus (RSV) and influenza viruses are important global causes of morbidity and mortality. We evaluated the diagnostic accuracy of the Luminex NxTAG respiratory pathogen panels (RPPs)™ (index) against other RPPs (comparator) for detection of RSV and influenza viruses. Studies comparing human clinical respiratory samples tested with the index and at least one comparator test were included. A random-effect latent class meta-analysis was performed to assess the specificity and sensitivity of the index test for RSV and influenza. Risk of bias was assessed using the QUADAS-2 tool and certainty of evidence using GRADE. Ten studies were included. For RSV, predicted sensitivity was 99% (95% credible interval [CrI] 96-100%) and specificity 100% (95% CrI 98-100%). For influenza A and B, predicted sensitivity was 97% (95% CrI 89-100) and 98% (95% CrI 88-100) respectively; specificity 100% (95% CrI 99-100) and 100% (95% CrI 99-100), respectively. Evidence was low certainty. Although index sensitivity and specificity were excellent, comparators' performance varied. Further research with clear patient recruitment strategies could ascertain performance across different populations.Protocol Registration: Prospero CRD42021272062.

Survival estimates in European cystic fibrosis patients and the impact of socioeconomic factors: a retrospective registry cohort study.

BACKGROUND: Median survival for cystic fibrosis (CF) patients in Europe is unknown and is likely to be influenced by socioeconomic factors. Using the European CF Society Patient Registry (ECFSPR), median survival estimates were obtained for CF patients across Europe and the impact of socioeconomic status on survival was examined. METHODS: CF subjects known to be alive and in the ECFSPR between 2010 and 2014 were included. Survival curves were estimated using the Kaplan-Meier method. Differences in the survival curves were assessed using the log-rank test. Cox regression was used to estimate the association between socioeconomic factors and the age-specific hazard of death, with adjustment for sex, age at diagnosis, CF transmembrane conductance regulator (CFTR) genotype and transplant status. RESULTS: The final analysis included 13 countries with 31 987 subjects (135 833 person-years of follow-up) and 1435 deaths. Median survival age for these patients in the ECFSPR was 51.7 (95% CI 50.0-53.4) years. After adjusting for potential confounders age at diagnosis, sex, CFTR genotype and transplant status, there remained strong evidence of an association between socioeconomic factors and mortality (p<0.001). Countries in the highest third of healthcare spending had a 46% lower hazard of mortality (HR 0.54, 95% CI 0.45-0.64) than countries in the lowest third of healthcare spending. CONCLUSIONS: Median survival for patients with CF in Europe is comparable to that reported in other jurisdictions and differs by socioeconomic factors.

Using generalized linear models to implement g-estimation for survival data with time-varying confounding.

Using data from observational studies to estimate the causal effect of a time-varying exposure, repeatedly measured over time, on an outcome of interest requires careful adjustment for confounding. Standard regression adjustment for observed time-varying confounders is unsuitable, as it can eliminate part of the causal effect and induce bias. Inverse probability weighting, g-computation, and g-estimation have been proposed as being more suitable methods. G-estimation has some advantages over the other two methods, but until recently there has been a lack of flexible g-estimation methods for a survival time outcome. The recently proposed Structural Nested Cumulative Survival Time Model (SNCSTM) is such a method. Efficient estimation of the parameters of this model required bespoke software. In this article we show how the SNCSTM can be fitted efficiently via g-estimation using standard software for fitting generalised linear models. The ability to implement g-estimation for a survival outcome using standard statistical software greatly increases the potential uptake of this method. We illustrate the use of this method of fitting the SNCSTM by reanalyzing data from the UK Cystic Fibrosis Registry, and provide example R code to facilitate the use of this approach by other researchers.

Importance of patient bed pathways and length of stay differences in predicting COVID-19 hospital bed occupancy in England.

BACKGROUND: Predicting bed occupancy for hospitalised patients with COVID-19 requires understanding of length of stay (LoS) in particular bed types. LoS can vary depending on the patient's "bed pathway" - the sequence of transfers of individual patients between bed types during a hospital stay. In this study, we characterise these pathways, and their impact on predicted hospital bed occupancy. METHODS: We obtained data from University College Hospital (UCH) and the ISARIC4C COVID-19 Clinical Information Network (CO-CIN) on hospitalised patients with COVID-19 who required care in general ward or critical care (CC) beds to determine possible bed pathways and LoS. We developed a discrete-time model to examine the implications of using either bed pathways or only average LoS by bed type to forecast bed occupancy. We compared model-predicted bed occupancy to publicly available bed occupancy data on COVID-19 in England between March and August 2020. RESULTS: In both the UCH and CO-CIN datasets, 82% of hospitalised patients with COVID-19 only received care in general ward beds. We identified four other bed pathways, present in both datasets: "Ward, CC, Ward", "Ward, CC", "CC" and "CC, Ward". Mean LoS varied by bed type, pathway, and dataset, between 1.78 and 13.53 days. For UCH, we found that using bed pathways improved the accuracy of bed occupancy predictions, while only using an average LoS for each bed type underestimated true bed occupancy. However, using the CO-CIN LoS dataset we were not able to replicate past data on bed occupancy in England, suggesting regional LoS heterogeneities. CONCLUSIONS: We identified five bed pathways, with substantial variation in LoS by bed type, pathway, and geography. This might be caused by local differences in patient characteristics, clinical care strategies, or resource availability, and suggests that national LoS averages may not be appropriate for local forecasts of bed occupancy for COVID-19. TRIAL REGISTRATION: The ISARIC WHO CCP-UK study ISRCTN66726260 was retrospectively registered on 21/04/2020 and designated an Urgent Public Health Research Study by NIHR.