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.

To Adjust or Not to Adjust? When a "Confounder" Is Only Measured After Exposure.

Advice regarding the analysis of observational studies of exposure effects usually is against adjustment for factors that occur after the exposure, as they may be caused by the exposure (or mediate the effect of exposure on outcome), so potentially leading to collider stratification bias. However, such factors could also be caused by unmeasured confounding factors, in which case adjusting for them will also remove some of the bias due to confounding. We derive expressions for collider stratification bias when conditioning and confounding bias when not conditioning on the mediator, in the presence of unmeasured confounding (assuming that all associations are linear and there are no interactions). Using simulations, we show that generally neither the conditioned nor the unconditioned estimate is unbiased, and the trade-off between them depends on the magnitude of the effect of the exposure that is mediated relative to the effect of the unmeasured confounders and their relations with the mediator. We illustrate the use of the bias expressions via three examples: neuroticism and mortality (adjusting for the mediator appears the least biased option), glycated hemoglobin levels and systolic blood pressure (adjusting gives smaller bias), and literacy in primary school pupils (not adjusting gives smaller bias). Our formulae and simulations can inform quantitative bias analysis as well as analysis strategies for observational studies in which there is a potential for unmeasured confounding.

Evaluating the relationship between circulating lipoprotein lipids and apolipoproteins with risk of coronary heart disease: A multivariable Mendelian randomisation analysis.

BACKGROUND: Circulating lipoprotein lipids cause coronary heart disease (CHD). However, the precise way in which one or more lipoprotein lipid-related entities account for this relationship remains unclear. Using genetic instruments for lipoprotein lipid traits implemented through multivariable Mendelian randomisation (MR), we sought to compare their causal roles in the aetiology of CHD. METHODS AND FINDINGS: We conducted a genome-wide association study (GWAS) of circulating non-fasted lipoprotein lipid traits in the UK Biobank (UKBB) for low-density lipoprotein (LDL) cholesterol, triglycerides, and apolipoprotein B to identify lipid-associated single nucleotide polymorphisms (SNPs). Using data from CARDIoGRAMplusC4D for CHD (consisting of 60,801 cases and 123,504 controls), we performed univariable and multivariable MR analyses. Similar GWAS and MR analyses were conducted for high-density lipoprotein (HDL) cholesterol and apolipoprotein A-I. The GWAS of lipids and apolipoproteins in the UKBB included between 393,193 and 441,016 individuals in whom the mean age was 56.9 y (range 39-73 y) and of whom 54.2% were women. The mean (standard deviation) lipid concentrations were LDL cholesterol 3.57 (0.87) mmol/L and HDL cholesterol 1.45 (0.38) mmol/L, and the median triglycerides was 1.50 (IQR = 1.11) mmol/L. The mean (standard deviation) values for apolipoproteins B and A-I were 1.03 (0.24) g/L and 1.54 (0.27) g/L, respectively. The GWAS identified multiple independent SNPs associated at P < 5 × 10-8 for LDL cholesterol (220), apolipoprotein B (n = 255), triglycerides (440), HDL cholesterol (534), and apolipoprotein A-I (440). Between 56%-93% of SNPs identified for each lipid trait had not been previously reported in large-scale GWASs. Almost half (46%) of these SNPs were associated at P < 5 × 10-8 with more than one lipid-related trait. Assessed individually using MR, LDL cholesterol (odds ratio [OR] 1.66 per 1-standard-deviation-higher trait; 95% CI: 1.49-1.86; P < 0.001), triglycerides (OR 1.34; 95% CI: 1.25-1.44; P < 0.001) and apolipoprotein B (OR 1.73; 95% CI: 1.56-1.91; P < 0.001) had effect estimates consistent with a higher risk of CHD. In multivariable MR, only apolipoprotein B (OR 1.92; 95% CI: 1.31-2.81; P < 0.001) retained a robust effect, with the estimate for LDL cholesterol (OR 0.85; 95% CI: 0.57-1.27; P = 0.44) reversing and that of triglycerides (OR 1.12; 95% CI: 1.02-1.23; P = 0.01) becoming weaker. Individual MR analyses showed a 1-standard-deviation-higher HDL cholesterol (OR 0.80; 95% CI: 0.75-0.86; P < 0.001) and apolipoprotein A-I (OR 0.83; 95% CI: 0.77-0.89; P < 0.001) to lower the risk of CHD, but these effect estimates attenuated substantially to the null on accounting for apolipoprotein B. A limitation is that, owing to the nature of lipoprotein metabolism, measures related to the composition of lipoprotein particles are highly correlated, creating a challenge in making exclusive interpretations on causation of individual components. CONCLUSIONS: These findings suggest that apolipoprotein B is the predominant trait that accounts for the aetiological relationship of lipoprotein lipids with risk of CHD.

Causal Associations of Adiposity and Body Fat Distribution With Coronary Heart Disease, Stroke Subtypes, and Type 2 Diabetes Mellitus: A Mendelian Randomization Analysis.

BACKGROUND: The implications of different adiposity measures on cardiovascular disease etiology remain unclear. In this article, we quantify and contrast causal associations of central adiposity (waist-to-hip ratio adjusted for body mass index [WHRadjBMI]) and general adiposity (body mass index [BMI]) with cardiometabolic disease. METHODS: Ninety-seven independent single-nucleotide polymorphisms for BMI and 49 single-nucleotide polymorphisms for WHRadjBMI were used to conduct Mendelian randomization analyses in 14 prospective studies supplemented with coronary heart disease (CHD) data from CARDIoGRAMplusC4D (Coronary Artery Disease Genome-wide Replication and Meta-analysis [CARDIoGRAM] plus The Coronary Artery Disease [C4D] Genetics; combined total 66 842 cases), stroke from METASTROKE (12 389 ischemic stroke cases), type 2 diabetes mellitus from DIAGRAM (Diabetes Genetics Replication and Meta-analysis; 34 840 cases), and lipids from GLGC (Global Lipids Genetic Consortium; 213 500 participants) consortia. Primary outcomes were CHD, type 2 diabetes mellitus, and major stroke subtypes; secondary analyses included 18 cardiometabolic traits. RESULTS: Each one standard deviation (SD) higher WHRadjBMI (1 SD≈0.08 U) associated with a 48% excess risk of CHD (odds ratio [OR] for CHD, 1.48; 95% confidence interval [CI], 1.28-1.71), similar to findings for BMI (1 SD≈4.6 kg/m2; OR for CHD, 1.36; 95% CI, 1.22-1.52). Only WHRadjBMI increased risk of ischemic stroke (OR, 1.32; 95% CI, 1.03-1.70). For type 2 diabetes mellitus, both measures had large effects: OR, 1.82 (95% CI, 1.38-2.42) and OR, 1.98 (95% CI, 1.41-2.78) per 1 SD higher WHRadjBMI and BMI, respectively. Both WHRadjBMI and BMI were associated with higher left ventricular hypertrophy, glycemic traits, interleukin 6, and circulating lipids. WHRadjBMI was also associated with higher carotid intima-media thickness (39%; 95% CI, 9%-77% per 1 SD). CONCLUSIONS: Both general and central adiposity have causal effects on CHD and type 2 diabetes mellitus. Central adiposity may have a stronger effect on stroke risk. Future estimates of the burden of adiposity on health should include measures of central and general adiposity.

Reading comprehension difficulties in children with rolandic epilepsy.

AIM: Difficulties in reading comprehension can arise from either word reading or listening comprehension difficulties, or a combination of the two. We sought to determine whether children with rolandic epilepsy had poor reading comprehension relative to typically developing comparison children, and whether such difficulties were associated with word reading and/or general language comprehension difficulties. METHOD: In this cross-sectional study, children with rolandic epilepsy (n=25; 16 males, 9 females; mean age 9y 1mo, SD 1y 7mo) and a comparison group (n=39; 25 males, 14 females; mean age 9y 1mo, SD 1y 3mo) completed assessments of reading comprehension, listening comprehension, word/non-word reading, speech articulation, and Non-verbal IQ. RESULTS: Reading comprehension and word reading were worse in children with rolandic epilepsy (F1,61 =6.89, p=0.011, ηp2=0.10 and F1,61 =6.84, p=0.011, ηp2=0.10 respectively), with listening comprehension being marginal (F1,61 =3.81, p=0.055, ηp2=0.06). Word reading and listening comprehension made large and independent contributions to reading comprehension, explaining 70% of the variance. INTERPRETATION: Children with rolandic epilepsy may be at risk of reading comprehension difficulties. Thorough assessment of individual children is required to ascertain whether the difficulties lie with decoding text, or with general comprehension skills, or both. WHAT THIS PAPER ADDS: Children with rolandic epilepsy may be at risk of poor reading comprehension. This was related to poor word reading, poor listening comprehension, or both. Reading comprehension interventions should be tailored to the profile of difficulties.

Correcting the Standard Errors of 2-Stage Residual Inclusion Estimators for Mendelian Randomization Studies.

Mendelian randomization studies use genotypes as instrumental variables to test for and estimate the causal effects of modifiable risk factors on outcomes. Two-stage residual inclusion (TSRI) estimators have been used when researchers are willing to make parametric assumptions. However, researchers are currently reporting uncorrected or heteroscedasticity-robust standard errors for these estimates. We compared several different forms of the standard error for linear and logistic TSRI estimates in simulations and in real-data examples. Among others, we consider standard errors modified from the approach of Newey (1987), Terza (2016), and bootstrapping. In our simulations Newey, Terza, bootstrap, and corrected 2-stage least squares (in the linear case) standard errors gave the best results in terms of coverage and type I error. In the real-data examples, the Newey standard errors were 0.5% and 2% larger than the unadjusted standard errors for the linear and logistic TSRI estimators, respectively. We show that TSRI estimators with modified standard errors have correct type I error under the null. Researchers should report TSRI estimates with modified standard errors instead of reporting unadjusted or heteroscedasticity-robust standard errors.

Harnessing publicly available genetic data to prioritize lipid modifying therapeutic targets for prevention of coronary heart disease based on dysglycemic risk.

Therapeutic interventions that lower LDL-cholesterol effectively reduce the risk of coronary artery disease (CAD). However, statins, the most widely prescribed LDL-cholesterol lowering drugs, increase diabetes risk. We used genome-wide association study (GWAS) data in the public domain to investigate the relationship of LDL-C and diabetes and identify loci encoding potential drug targets for LDL-cholesterol modification without causing dysglycemia. We obtained summary-level GWAS data for LDL-C from GLGC, glycemic traits from MAGIC, diabetes from DIAGRAM and CAD from CARDIoGRAMplusC4D consortia. Mendelian randomization analyses identified a one standard deviation (SD) increase in LDL-C caused an increased risk of CAD (odds ratio [OR] 1.63 (95 % confidence interval [CI] 1.55, 1.71), which was not influenced by removing SNPs associated with diabetes. LDL-C/CAD-associated SNPs showed consistent effect directions (binomial P = 6.85 × 10(-5)). Conversely, a 1-SD increase in LDL-C was causally protective of diabetes (OR 0.86; 95 % CI 0.81, 0.91), however LDL-cholesterol/diabetes-associated SNPs did not show consistent effect directions (binomial P = 0.15). HMGCR, our positive control, associated with LDL-C, CAD and a glycemic composite (derived from GWAS meta-analysis of four glycemic traits and diabetes). In contrast, PCSK9, APOB, LPA, CETP, PLG, NPC1L1 and ALDH2 were identified as "druggable" loci that alter LDL-C and risk of CAD without displaying associations with dysglycemia. In conclusion, LDL-C increases the risk of CAD and the relationship is independent of any association of LDL-C with diabetes. Loci that encode targets of emerging LDL-C lowering drugs do not associate with dysglycemia, and this provides provisional evidence that new LDL-C lowering drugs (such as PCSK9 inhibitors) may not influence risk of diabetes.

Mendelian randomization of blood lipids for coronary heart disease.

AIMS: To investigate the causal role of high-density lipoprotein cholesterol (HDL-C) and triglycerides in coronary heart disease (CHD) using multiple instrumental variables for Mendelian randomization. METHODS AND RESULTS: We developed weighted allele scores based on single nucleotide polymorphisms (SNPs) with established associations with HDL-C, triglycerides, and low-density lipoprotein cholesterol (LDL-C). For each trait, we constructed two scores. The first was unrestricted, including all independent SNPs associated with the lipid trait identified from a prior meta-analysis (threshold P < 2 × 10(-6)); and the second a restricted score, filtered to remove any SNPs also associated with either of the other two lipid traits at P ≤ 0.01. Mendelian randomization meta-analyses were conducted in 17 studies including 62,199 participants and 12,099 CHD events. Both the unrestricted and restricted allele scores for LDL-C (42 and 19 SNPs, respectively) associated with CHD. For HDL-C, the unrestricted allele score (48 SNPs) was associated with CHD (OR: 0.53; 95% CI: 0.40, 0.70), per 1 mmol/L higher HDL-C, but neither the restricted allele score (19 SNPs; OR: 0.91; 95% CI: 0.42, 1.98) nor the unrestricted HDL-C allele score adjusted for triglycerides, LDL-C, or statin use (OR: 0.81; 95% CI: 0.44, 1.46) showed a robust association. For triglycerides, the unrestricted allele score (67 SNPs) and the restricted allele score (27 SNPs) were both associated with CHD (OR: 1.62; 95% CI: 1.24, 2.11 and 1.61; 95% CI: 1.00, 2.59, respectively) per 1-log unit increment. However, the unrestricted triglyceride score adjusted for HDL-C, LDL-C, and statin use gave an OR for CHD of 1.01 (95% CI: 0.59, 1.75). CONCLUSION: The genetic findings support a causal effect of triglycerides on CHD risk, but a causal role for HDL-C, though possible, remains less certain.

Genetic Evidence for Causal Relationships Between Maternal Obesity-Related Traits and Birth Weight.

IMPORTANCE: Neonates born to overweight or obese women are larger and at higher risk of birth complications. Many maternal obesity-related traits are observationally associated with birth weight, but the causal nature of these associations is uncertain. OBJECTIVE: To test for genetic evidence of causal associations of maternal body mass index (BMI) and related traits with birth weight. DESIGN, SETTING, AND PARTICIPANTS: Mendelian randomization to test whether maternal BMI and obesity-related traits are potentially causally related to offspring birth weight. Data from 30,487 women in 18 studies were analyzed. Participants were of European ancestry from population- or community-based studies in Europe, North America, or Australia and were part of the Early Growth Genetics Consortium. Live, term, singleton offspring born between 1929 and 2013 were included. EXPOSURES: Genetic scores for BMI, fasting glucose level, type 2 diabetes, systolic blood pressure (SBP), triglyceride level, high-density lipoprotein cholesterol (HDL-C) level, vitamin D status, and adiponectin level. MAIN OUTCOME AND MEASURE: Offspring birth weight from 18 studies. RESULTS: Among the 30,487 newborns the mean birth weight in the various cohorts ranged from 3325 g to 3679 g. The maternal genetic score for BMI was associated with a 2-g (95% CI, 0 to 3 g) higher offspring birth weight per maternal BMI-raising allele (P = .008). The maternal genetic scores for fasting glucose and SBP were also associated with birth weight with effect sizes of 8 g (95% CI, 6 to 10 g) per glucose-raising allele (P = 7 × 10(-14)) and -4 g (95% CI, -6 to -2 g) per SBP-raising allele (P = 1×10(-5)), respectively. A 1-SD ( ≈ 4 points) genetically higher maternal BMI was associated with a 55-g higher offspring birth weight (95% CI, 17 to 93 g). A 1-SD ( ≈ 7.2 mg/dL) genetically higher maternal fasting glucose concentration was associated with 114-g higher offspring birth weight (95% CI, 80 to 147 g). However, a 1-SD ( ≈ 10 mm Hg) genetically higher maternal SBP was associated with a 208-g lower offspring birth weight (95% CI, -394 to -21 g). For BMI and fasting glucose, genetic associations were consistent with the observational associations, but for systolic blood pressure, the genetic and observational associations were in opposite directions. CONCLUSIONS AND RELEVANCE: In this mendelian randomization study, genetically elevated maternal BMI and blood glucose levels were potentially causally associated with higher offspring birth weight, whereas genetically elevated maternal SBP was potentially causally related to lower birth weight. If replicated, these findings may have implications for counseling and managing pregnancies to avoid adverse weight-related birth outcomes.

Effects of BMI, fat mass, and lean mass on asthma in childhood: a Mendelian randomization study.

BACKGROUND: Observational studies have reported associations between body mass index (BMI) and asthma, but confounding and reverse causality remain plausible explanations. We aim to investigate evidence for a causal effect of BMI on asthma using a Mendelian randomization approach. METHODS AND FINDINGS: We used Mendelian randomization to investigate causal effects of BMI, fat mass, and lean mass on current asthma at age 7½ y in the Avon Longitudinal Study of Parents and Children (ALSPAC). A weighted allele score based on 32 independent BMI-related single nucleotide polymorphisms (SNPs) was derived from external data, and associations with BMI, fat mass, lean mass, and asthma were estimated. We derived instrumental variable (IV) estimates of causal risk ratios (RRs). 4,835 children had available data on BMI-associated SNPs, asthma, and BMI. The weighted allele score was strongly associated with BMI, fat mass, and lean mass (all p-values<0.001) and with childhood asthma (RR 2.56, 95% CI 1.38-4.76 per unit score, p = 0.003). The estimated causal RR for the effect of BMI on asthma was 1.55 (95% CI 1.16-2.07) per kg/m2, p = 0.003. This effect appeared stronger for non-atopic (1.90, 95% CI 1.19-3.03) than for atopic asthma (1.37, 95% CI 0.89-2.11) though there was little evidence of heterogeneity (p = 0.31). The estimated causal RRs for the effects of fat mass and lean mass on asthma were 1.41 (95% CI 1.11-1.79) per 0.5 kg and 2.25 (95% CI 1.23-4.11) per kg, respectively. The possibility of genetic pleiotropy could not be discounted completely; however, additional IV analyses using FTO variant rs1558902 and the other BMI-related SNPs separately provided similar causal effects with wider confidence intervals. Loss of follow-up was unlikely to bias the estimated effects. CONCLUSIONS: Higher BMI increases the risk of asthma in mid-childhood. Higher BMI may have contributed to the increase in asthma risk toward the end of the 20th century. Please see later in the article for the Editors' Summary.

Lack of identification in semiparametric instrumental variable models with binary outcomes.

A parameter in a statistical model is identified if its value can be uniquely determined from the distribution of the observable data. We consider the context of an instrumental variable analysis with a binary outcome for estimating a causal risk ratio. The semiparametric generalized method of moments and structural mean model frameworks use estimating equations for parameter estimation. In this paper, we demonstrate that lack of identification can occur in either of these frameworks, especially if the instrument is weak. In particular, the estimating equations may have no solution or multiple solutions. We investigate the relationship between the strength of the instrument and the proportion of simulated data sets for which there is a unique solution of the estimating equations. We see that this proportion does not appear to depend greatly on the sample size, particularly for weak instruments (ρ(2) ≤ 0.01). Poor identification was observed in a considerable proportion of simulated data sets for instruments explaining up to 10% of the variance in the exposure with sample sizes up to 1 million. In an applied example considering the causal effect of body mass index (weight (kg)/height (m)(2)) on the probability of early menarche, estimates and standard errors from an automated optimization routine were misleading.

Severity of bias of a simple estimator of the causal odds ratio in Mendelian randomization studies.

Mendelian randomization studies estimate causal effects using genetic variants as instruments. Instrumental variable methods are straightforward for linear models, but epidemiologists often use odds ratios to quantify effects. Also, odds ratios are often the quantities reported in meta-analyses. Many applications of Mendelian randomization dichotomize genotype and estimate the population causal log odds ratio for unit increase in exposure by dividing the genotype-disease log odds ratio by the difference in mean exposure between genotypes. This 'Wald-type' estimator is biased even in large samples, but whether the magnitude of bias is of practical importance is unclear. We study the large-sample bias of this estimator in a simple model with a continuous normally distributed exposure, a single unobserved confounder that is not an effect modifier, and interpretable parameters. We focus on parameter values that reflect scenarios in which we apply Mendelian randomization, including realistic values for the degree of confounding and strength of the causal effect. We evaluate this estimator and the causal odds ratio using numerical integration and obtain approximate analytic expressions to check results and gain insight. A small simulation study examines finite sample bias and mild violations of the normality assumption. For our simple data-generating model, we find that the Wald estimator is asymptotically biased with a bias of around 10% in fairly typical Mendelian randomization scenarios but which can be larger in more extreme situations. Recently developed methods such as structural mean models require fewer untestable assumptions and we recommend their use when the individual-level data they require are available. The Wald-type estimator may retain a role as an approximate method for meta-analysis based on summary data.

Comparative effectiveness of BNT162b2 versus mRNA-1273 covid-19 vaccine boosting in England: matched cohort study in OpenSAFELY-TPP.

OBJECTIVE: To compare the effectiveness of the BNT162b2 mRNA (Pfizer-BioNTech) and mRNA-1273 (Moderna) covid-19 vaccines during the booster programme in England. DESIGN: Matched cohort study, emulating a comparative effectiveness trial. SETTING: Linked primary care, hospital, and covid-19 surveillance records available within the OpenSAFELY-TPP research platform, covering a period when the SARS-CoV-2 delta and omicron variants were dominant. PARTICIPANTS: 3 237 918 adults who received a booster dose of either vaccine between 29 October 2021 and 25 February 2022 as part of the national booster programme in England and who received a primary course of BNT162b2 or ChAdOx1. INTERVENTION: Vaccination with either BNT162b2 or mRNA-1273 as a booster vaccine dose. MAIN OUTCOME MEASURES: Recorded SARS-CoV-2 positive test, covid-19 related hospital admission, covid-19 related death, and non-covid-19 related death at 20 weeks after receipt of the booster dose. RESULTS: 1 618 959 people were matched in each vaccine group, contributing a total 64 546 391 person weeks of follow-up. The 20 week risks per 1000 for a positive SARS-CoV-2 test were 164.2 (95% confidence interval 163.3 to 165.1) for BNT162b2 and 159.9 (159.0 to 160.8) for mRNA-1273; the hazard ratio comparing mRNA-1273 with BNT162b2 was 0.95 (95% confidence interval 0.95 to 0.96). The 20 week risks per 1000 for hospital admission with covid-19 were 0.75 (0.71 to 0.79) for BNT162b2 and 0.65 (0.61 to 0.69) for mRNA-1273; the hazard ratio was 0.89 (0.82 to 0.95). Covid-19 related deaths were rare: the 20 week risks per 1000 were 0.028 (0.021 to 0.037) for BNT162b2 and 0.024 (0.018 to 0.033) for mRNA-1273; hazard ratio 0.83 (0.58 to 1.19). Comparative effectiveness was generally similar within subgroups defined by the primary course vaccine brand, age, previous SARS-CoV-2 infection, and clinical vulnerability. Relative benefit was similar when vaccines were compared separately in the delta and omicron variant eras. CONCLUSIONS: This matched observational study of adults estimated a modest benefit of booster vaccination with mRNA-1273 compared with BNT162b2 in preventing positive SARS-CoV-2 tests and hospital admission with covid-19 20 weeks after vaccination, during a period of delta followed by omicron variant dominance.

The effect of elevated body mass index on ischemic heart disease risk: causal estimates from a Mendelian randomisation approach.

BACKGROUND: Adiposity, assessed as elevated body mass index (BMI), is associated with increased risk of ischemic heart disease (IHD); however, whether this is causal is unknown. We tested the hypothesis that positive observational associations between BMI and IHD are causal. METHODS AND FINDINGS: In 75,627 individuals taken from two population-based and one case-control study in Copenhagen, we measured BMI, ascertained 11,056 IHD events, and genotyped FTO(rs9939609), MC4R(rs17782313), and TMEM18(rs6548238). Using genotypes as a combined allele score in instrumental variable analyses, the causal odds ratio (OR) between BMI and IHD was estimated and compared with observational estimates. The allele score-BMI and the allele score-IHD associations used to estimate the causal OR were also calculated individually. In observational analyses the OR for IHD was 1.26 (95% CI 1.19-1.34) for every 4 kg/m(2) increase in BMI. A one-unit allele score increase associated with a 0.28 kg/m(2) (95 CI% 0.20-0.36) increase in BMI and an OR for IHD of 1.03 (95% CI 1.01-1.05) (corresponding to an average 1.68 kg/m(2) BMI increase and 18% increase in the odds of IHD for those carrying all six BMI increasing alleles). In instrumental variable analysis using the same allele score the causal IHD OR for a 4 kg/m(2) increase in BMI was 1.52 (95% CI 1.12-2.05). CONCLUSIONS: For every 4 kg/m(2) increase in BMI, observational estimates suggested a 26% increase in odds for IHD while causal estimates suggested a 52% increase. These data add evidence to support a causal link between increased BMI and IHD risk, though the mechanism may ultimately be through intermediate factors like hypertension, dyslipidemia, and type 2 diabetes. This work has important policy implications for public health, given the continuous nature of the BMI-IHD association and the modifiable nature of BMI. This analysis demonstrates the value of observational studies and their ability to provide unbiased results through inclusion of genetic data avoiding confounding, reverse causation, and bias.

Waning effectiveness of BNT162b2 and ChAdOx1 covid-19 vaccines over six months since second dose: OpenSAFELY cohort study using linked electronic health records.

OBJECTIVE: To estimate waning of covid-19 vaccine effectiveness over six months after second dose. DESIGN: Cohort study, approved by NHS England. SETTING: Linked primary care, hospital, and covid-19 records within the OpenSAFELY-TPP database. PARTICIPANTS: Adults without previous SARS-CoV-2 infection were eligible, excluding care home residents and healthcare professionals. EXPOSURES: People who had received two doses of BNT162b2 or ChAdOx1 (administered during the national vaccine rollout) were compared with unvaccinated people during six consecutive comparison periods, each of four weeks. MAIN OUTCOME MEASURES: Adjusted hazard ratios for covid-19 related hospital admission, covid-19 related death, positive SARS-CoV-2 test, and non-covid-19 related death comparing vaccinated with unvaccinated people. Waning vaccine effectiveness was quantified as ratios of adjusted hazard ratios per four week period, separately for subgroups aged ≥65 years, 18-64 years and clinically vulnerable, 40-64 years, and 18-39 years. RESULTS: 1 951 866 and 3 219 349 eligible adults received two doses of BNT162b2 and ChAdOx1, respectively, and 2 422 980 remained unvaccinated. Waning of vaccine effectiveness was estimated to be similar across outcomes and vaccine brands. In the ≥65 years subgroup, ratios of adjusted hazard ratios for covid-19 related hospital admission, covid-19 related death, and positive SARS-CoV-2 test ranged from 1.19 (95% confidence interval 1.14 to 1.24)to 1.34 (1.09 to 1.64) per four weeks. Despite waning vaccine effectiveness, rates of covid-19 related hospital admission and death were substantially lower among vaccinated than unvaccinated adults up to 26 weeks after the second dose, with estimated vaccine effectiveness ≥80% for BNT162b2, and ≥75% for ChAdOx1. By weeks 23-26, rates of positive SARS-CoV-2 test in vaccinated people were similar to or higher than in unvaccinated people (adjusted hazard ratios up to 1.72 (1.11 to 2.68) for BNT162b2 and 1.86 (1.79 to 1.93) for ChAdOx1). CONCLUSIONS: The rate at which estimated vaccine effectiveness waned was consistent for covid-19 related hospital admission, covid-19 related death, and positive SARS-CoV-2 test and was similar across subgroups defined by age and clinical vulnerability. If sustained to outcomes of infection with the omicron variant and to booster vaccination, these findings will facilitate scheduling of booster vaccination.

Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses.

In this paper, the authors describe different instrumental variable (IV) estimators of causal risk ratios and odds ratios with particular attention to methods that can handle continuously measured exposures. The authors present this discussion in the context of a Mendelian randomization analysis of the effect of body mass index (BMI; weight (kg)/height (m)(2)) on the risk of asthma at age 7 years (Avon Longitudinal Study of Parents and Children, 1991-1992). The authors show that the multiplicative structural mean model (MSMM) and the multiplicative generalized method of moments (MGMM) estimator produce identical estimates of the causal risk ratio. In the example, MSMM and MGMM estimates suggested an inverse relation between BMI and asthma but other IV estimates suggested a positive relation, although all estimates had wide confidence intervals. An interaction between the associations of BMI and fat mass and obesity-associated (FTO) genotype with asthma explained the different directions of the different estimates, and a simulation study supported the observation that MSMM/MGMM estimators are negatively correlated with the other estimators when such an interaction is present. The authors conclude that point estimates from various IV methods can differ in practical applications. Based on the theoretical properties of the estimators, structural mean models make weaker assumptions than other IV estimators and can therefore be expected to be consistent in a wider range of situations.

Using genetic loci to understand the relationship between adiposity and psychological distress: a Mendelian Randomization study in the Copenhagen General Population Study of 53,221 adults.

OBJECTIVE: We used genetic variants that are robustly associated with adiposity to examine the causal association of adiposity with psychological distress. METHODS: We examined the association of adiposity with psychological distress in a large (N = 53,221) general population cohort of 20- to 99-year-old adults from Copenhagen, Denmark. Psychological distress was assessed using four questions that asked about: feeling stressed; not accomplishing very much; wanting to give up; and regular use of antidepressants/sedatives. We used the genetic loci FTO rs9939609 and MC4R rs17782313 as instrumental variables for adiposity quantified by body mass index (BMI) and waist to hip ratio (WHR). RESULTS: In conventional multivariable analyses, BMI and WHR were positively associated with distress. For example, the odds ratio of reporting not accomplishing for each additional standard deviation increase for BMI was 1.11 (95% CI: 1.09, 1.13) and for WHR was 1.10 (95% CI: 1.08, 1.13) in the fully adjusted analyses. In contrast, instrumental variable analyses showed an inverse association of adiposity on distress; corresponding odds ratio in instrumental variable analyses was 0.64 (95% CI: 0.46, 0.89) for BMI and 0.49 (95% CI: 0.25, 0.94) for WHR (P-values for difference between the two approaches both = 0.001). CONCLUSION: The inverse associations of adiposity and psychological distress when genetic variants are used as instrumental variables could be explained by biological pathways linking adiposity and distress. The positive associations of adiposity with distress in multivariable analyses might be explained by residual confounding or reverse causality.

The role of combinatorial health technologies in supporting older people with long-term conditions: Responsibilisation or co-management of healthcare?

Neoliberalism, austerity and health responsibilisation are increasingly informing policies and practices designed to encourage older patients to take responsibility for the management of their own healthcare. Combined with an ageing population, novel ways to address the increasing healthcare needs of older people have become a priority, with the emergence in recent years of new models of integrated care enhanced by combinatorial health technologies (CHTs). This paper presents qualitative findings from the evaluation of one programme, the Lancashire and Cumbria Innovation Alliance (LCIA) Test Bed, a programme funded by NHS England and conducted in England between 2016 and 2018. Drawing on data from patients, family carers, and staff members involved in the programme, this paper explores the extent to which CHTs, as part of the LCIA Test Bed programme, contributed to health responsibilisation amongst older people with complex health conditions. Through this programme, we find that relationships between patients, family carers and healthcare professionals combined to create a sense of reassurance and shared responsibility for all parties. Our findings suggest the need for a more nuanced approach to responsibilisation and self-management for older people living with complex health conditions. By focusing on co-management - and recognising the potential of CHTs to facilitate this approach - there is potential to increase patient confidence in managing their health condition, reduce carer burden, and enhance clinician satisfaction in their work roles. While neoliberal agendas are focused on self-management and self-responsibility of one's own health care, with technology as a facilitator of this, our findings suggest that the successful use of CHTs for older people with complex health conditions may instead be rooted in co-management. This paper argues that co-management may be a more successful model of care for patients, carers and clinicians.

Collider bias undermines our understanding of COVID-19 disease risk and severity.

Numerous observational studies have attempted to identify risk factors for infection with SARS-CoV-2 and COVID-19 disease outcomes. Studies have used datasets sampled from patients admitted to hospital, people tested for active infection, or people who volunteered to participate. Here, we highlight the challenge of interpreting observational evidence from such non-representative samples. Collider bias can induce associations between two or more variables which affect the likelihood of an individual being sampled, distorting associations between these variables in the sample. Analysing UK Biobank data, compared to the wider cohort the participants tested for COVID-19 were highly selected for a range of genetic, behavioural, cardiovascular, demographic, and anthropometric traits. We discuss the mechanisms inducing these problems, and approaches that could help mitigate them. While collider bias should be explored in existing studies, the optimal way to mitigate the problem is to use appropriate sampling strategies at the study design stage.