ABSTRACT
Mendelian randomization (MR) is an epidemiological approach that utilizes genetic variants as instrumental variables to estimate the causal effect of an exposure on a health outcome. This paper investigates an MR scenario in which genetic variants aggregate into clusters that identify heterogeneous causal effects. Such variant clusters are likely to emerge if they affect the exposure and outcome via distinct biological pathways. In the multi-outcome MR framework, where a shared exposure causally impacts several disease outcomes simultaneously, these variant clusters can provide insights into the common disease-causing mechanisms underpinning the co-occurrence of multiple long-term conditions, a phenomenon known as multimorbidity. To identify such variant clusters, we adapt the general method of agglomerative hierarchical clustering to multi-sample summary-data MR setup, enabling cluster detection based on variant-specific ratio estimates. Particularly, we tailor the method for multi-outcome MR to aid in elucidating the causal pathways through which a common risk factor contributes to multiple morbidities. We show in simulations that our "MR-AHC" method detects clusters with high accuracy, outperforming the existing methods. We apply the method to investigate the causal effects of high body fat percentage on type 2 diabetes and osteoarthritis, uncovering interconnected cellular processes underlying this multimorbid disease pair.
ABSTRACT
Self-reported shorter/longer sleep duration, insomnia, and evening preference are associated with hyperglycaemia in observational analyses, with similar observations in small studies using accelerometer-derived sleep traits. Mendelian randomization (MR) studies support an effect of self-reported insomnia, but not others, on glycated haemoglobin (HbA1c). To explore potential effects, we used MR methods to assess effects of accelerometer-derived sleep traits (duration, mid-point least active 5-h, mid-point most active 10-h, sleep fragmentation, and efficiency) on HbA1c/glucose in European adults from the UK Biobank (UKB) (n = 73,797) and the MAGIC consortium (n = 146,806). Cross-trait linkage disequilibrium score regression was applied to determine genetic correlations across accelerometer-derived, self-reported sleep traits, and HbA1c/glucose. We found no causal effect of any accelerometer-derived sleep trait on HbA1c or glucose. Similar MR results for self-reported sleep traits in the UKB sub-sample with accelerometer-derived measures suggested our results were not explained by selection bias. Phenotypic and genetic correlation analyses suggested complex relationships between self-reported and accelerometer-derived traits indicating that they may reflect different types of exposure. These findings suggested accelerometer-derived sleep traits do not affect HbA1c. Accelerometer-derived measures of sleep duration and quality might not simply be 'objective' measures of self-reported sleep duration and insomnia, but rather captured different sleep characteristics.
Subject(s)
Accelerometry , Blood Glucose , Glycated Hemoglobin , Mendelian Randomization Analysis , Sleep , Humans , Glycated Hemoglobin/analysis , Glycated Hemoglobin/metabolism , Sleep/genetics , Sleep/physiology , Blood Glucose/analysis , Male , Female , Middle Aged , Adult , Self Report , Aged , Sleep Initiation and Maintenance Disorders/geneticsABSTRACT
AIMS/HYPOTHESIS: Older adults are under-represented in trials, meaning the benefits and risks of glucose-lowering agents in this age group are unclear. The aim of this study was to assess the safety and effectiveness of sodium-glucose cotransporter 2 inhibitors (SGLT2i) in people with type 2 diabetes aged over 70 years using causal analysis. METHODS: Hospital-linked UK primary care data (Clinical Practice Research Datalink, 2013-2020) were used to compare adverse events and effectiveness in individuals initiating SGLT2i compared with dipeptidyl peptidase-4 inhibitors (DPP4i). Analysis was age-stratified: <70 years (SGLT2i n=66,810, DPP4i n=76,172), ≥70 years (SGLT2i n=10,419, DPP4i n=33,434). Outcomes were assessed using the instrumental variable causal inference method and prescriber preference as the instrument. RESULTS: Risk of diabetic ketoacidosis was increased with SGLT2i in those aged ≥70 (incidence rate ratio compared with DPP4i: 3.82 [95% CI 1.12, 13.03]), but not in those aged <70 (1.12 [0.41, 3.04]). However, incidence rates with SGLT2i in those ≥70 was low (29.6 [29.5, 29.7]) per 10,000 person-years. SGLT2i were associated with similarly increased risk of genital infection in both age groups (incidence rate ratio in those <70: 2.27 [2.03, 2.53]; ≥70: 2.16 [1.77, 2.63]). There was no evidence of an increased risk of volume depletion, poor micturition control, urinary frequency, falls or amputation with SGLT2i in either age group. In those ≥70, HbA1c reduction was similar between SGLT2i and DPP4i (-0.3 mmol/mol [-1.6, 1.1], -0.02% [0.1, 0.1]), but in those <70, SGLT2i were more effective (-4 mmol/mol [4.8, -3.1], -0.4% [-0.4, -0.3]). CONCLUSIONS/INTERPRETATION: Causal analysis suggests SGLT2i are effective in adults aged ≥70 years, but increase risk for genital infections and diabetic ketoacidosis. Our study extends RCT evidence to older adults with type 2 diabetes.
ABSTRACT
The solute carrier organic anion transporter family member 1B1 (SLCO1B1) encodes the organic anion-transporting polypeptide 1B1 (OATP1B1 protein) that transports statins to liver cells. Common genetic variants in SLCO1B1, such as *5, cause altered systemic exposure to statins and therefore affect statin outcomes, with potential pharmacogenetic applications; yet, evidence is inconclusive. We studied common and rare SLCO1B1 variants in up to 64,000 patients from UK Biobank prescribed simvastatin or atorvastatin, combining whole-exome sequencing data with up to 25-year routine clinical records. We studied 51 predicted gain/loss-of-function variants affecting OATP1B1. Both SLCO1B1*5 alone and the SLCO1B1*15 haplotype increased LDL during treatment (beta*5 = 0.08 mmol/L, p = 6 × 10-8; beta*15 = 0.03 mmol/L, p = 3 × 10-4), as did the likelihood of discontinuing statin prescriptions (hazard ratio*5 = 1.12, p = 0.04; HR*15 = 1.05, p = 0.04). SLCO1B1*15 and SLCO1B1*20 increased the risk of General Practice (GP)-diagnosed muscle symptoms (HR*15 = 1.22, p = 0.003; HR*20 = 1.25, p = 0.01). We estimated that genotype-guided prescribing could potentially prevent 18% and 10% of GP-diagnosed muscle symptoms experienced by statin patients, with *15 and *20, respectively. The remaining common variants were not individually significant. Rare variants in SLCO1B1 increased LDL in statin users by up to 1.05 mmol/L, but replication is needed. We conclude that genotype-guided treatment could reduce GP-diagnosed muscle symptoms in statin patients; incorporating further SLCO1B1 variants into clinical prediction scores could improve LDL control and decrease adverse events, including discontinuation.
Subject(s)
Exome Sequencing , Hydroxymethylglutaryl-CoA Reductase Inhibitors , Liver-Specific Organic Anion Transporter 1 , Aged , Female , Humans , Male , Middle Aged , Atorvastatin/therapeutic use , Exome Sequencing/methods , Hydroxymethylglutaryl-CoA Reductase Inhibitors/therapeutic use , Hydroxymethylglutaryl-CoA Reductase Inhibitors/adverse effects , Liver-Specific Organic Anion Transporter 1/genetics , Polymorphism, Single Nucleotide , Simvastatin/therapeutic use , Treatment Outcome , UK Biobank , United KingdomABSTRACT
Pharmacogenetic variants are associated with clinical outcomes during Calcium Channel Blocker (CCB) treatment, yet whether the effects are modified by genetically predicted clinical risk factors is unknown. We analyzed 32,000 UK Biobank participants treated with dihydropiridine CCBs (mean 5.9 years), including 23 pharmacogenetic variants, and calculated polygenic scores for systolic and diastolic blood pressures, body fat mass, and other patient characteristics. Outcomes included treatment discontinuation and heart failure. Pharmacogenetic variant rs10898815-A (NUMA1) increased discontinuation rates, highest in those with high polygenic scores for fat mass. The RYR3 variant rs877087 T-allele alone modestly increased heart failure risks versus non-carriers (HR:1.13, p = 0.02); in patients with high polygenic scores for fat mass, lean mass, and lipoprotein A, risks were substantially elevated (HR:1.55, p = 4 × 10-5). Incorporating polygenic scores for adiposity and lipoprotein A may improve risk estimates of key clinical outcomes in CCB treatment such as treatment discontinuation and heart failure, compared to pharmacogenetic variants alone.
Subject(s)
Cardiovascular Diseases , Heart Failure , Hypertension , Humans , Calcium Channel Blockers/therapeutic use , Antihypertensive Agents/therapeutic use , Hypertension/drug therapy , Pharmacogenomic Variants , Cardiovascular Diseases/chemically induced , Risk Factors , Heart Failure/drug therapy , Heart Disease Risk Factors , Lipoprotein(a)/therapeutic useABSTRACT
BACKGROUND: The causal relationship between maternal smoking in pregnancy and reduced offspring birth weight is well established and is likely due to impaired placental function. However, observational studies have given conflicting results on the association between smoking and placental weight. We aimed to estimate the causal effect of newly pregnant mothers quitting smoking on their placental weight at the time of delivery. METHODS: We used one-sample Mendelian randomization, drawing data from the Avon Longitudinal Study of Parents and Children (ALSPAC) (N = 690 to 804) and the Norwegian Mother, Father and Child Cohort Study (MoBa) (N = 4267 to 4606). The sample size depends on the smoking definition used for different analyses. The analysis was performed in pre-pregnancy smokers only, due to the specific role of the single-nucleotide polymorphism (SNP) rs1051730 (CHRNA5 - CHRNA3 - CHRNB4) in affecting smoking cessation but not initiation. RESULTS: Fixed effect meta-analysis showed a 182 g [95%CI: 29,335] higher placental weight for pre-pregnancy smoking mothers who continued smoking at the beginning of pregnancy, compared with those who stopped smoking. Using the number of cigarettes smoked per day in the first trimester as the exposure, the causal effect on placental weight was 11 g [95%CI: 1,21] per cigarette per day. Similarly, smoking at the end of pregnancy was causally associated with higher placental weight. Using the residuals of birth weight regressed on placental weight as the outcome, we showed evidence of lower offspring birth weight relative to the placental weight, both for continuing smoking at the start of pregnancy as well as continuing smoking throughout pregnancy (change in z-score birth weight adjusted for z-score placental weight: -0.8 [95%CI: -1.6,-0.1]). CONCLUSION: Our results suggest that continued smoking during pregnancy causes higher placental weights.
Subject(s)
Mendelian Randomization Analysis , Placenta , Female , Humans , Pregnancy , Birth Weight/genetics , Cohort Studies , Longitudinal Studies , Smoking/adverse effectsABSTRACT
Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.
Subject(s)
Body Mass Index , Coronary Artery Disease , Diabetes Mellitus, Type 2 , Glucagon-Like Peptide-1 Receptor , Mendelian Randomization Analysis , Phenotype , Humans , Glucagon-Like Peptide-1 Receptor/genetics , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/drug therapy , Coronary Artery Disease/genetics , Coronary Artery Disease/drug therapy , Polymorphism, Single Nucleotide , Genetic Predisposition to DiseaseABSTRACT
AIMS/HYPOTHESIS: A precision medicine approach in type 2 diabetes could enhance targeting specific glucose-lowering therapies to individual patients most likely to benefit. We aimed to use the recently developed Bayesian causal forest (BCF) method to develop and validate an individualised treatment selection algorithm for two major type 2 diabetes drug classes, sodium-glucose cotransporter 2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP1-RA). METHODS: We designed a predictive algorithm using BCF to estimate individual-level conditional average treatment effects for 12-month glycaemic outcome (HbA1c) between SGLT2i and GLP1-RA, based on routine clinical features of 46,394 people with type 2 diabetes in primary care in England (Clinical Practice Research Datalink; 27,319 for model development, 19,075 for hold-out validation), with additional external validation in 2252 people with type 2 diabetes from Scotland (SCI-Diabetes [Tayside & Fife]). Differences in glycaemic outcome with GLP1-RA by sex seen in clinical data were replicated in clinical trial data (HARMONY programme: liraglutide [n=389] and albiglutide [n=1682]). As secondary outcomes, we evaluated the impacts of targeting therapy based on glycaemic response on weight change, tolerability and longer-term risk of new-onset microvascular complications, macrovascular complications and adverse kidney events. RESULTS: Model development identified marked heterogeneity in glycaemic response, with 4787 (17.5%) of the development cohort having a predicted HbA1c benefit >3 mmol/mol (>0.3%) with SGLT2i over GLP1-RA and 5551 (20.3%) having a predicted HbA1c benefit >3 mmol/mol with GLP1-RA over SGLT2i. Calibration was good in hold-back validation, and external validation in an independent Scottish dataset identified clear differences in glycaemic outcomes between those predicted to benefit from each therapy. Sex, with women markedly more responsive to GLP1-RA, was identified as a major treatment effect modifier in both the UK observational datasets and in clinical trial data: HARMONY-7 liraglutide (GLP1-RA): 4.4 mmol/mol (95% credible interval [95% CrI] 2.2, 6.3) (0.4% [95% CrI 0.2, 0.6]) greater response in women than men. Targeting the two therapies based on predicted glycaemic response was also associated with improvements in short-term tolerability and long-term risk of new-onset microvascular complications. CONCLUSIONS/INTERPRETATION: Precision medicine approaches can facilitate effective individualised treatment choice between SGLT2i and GLP1-RA therapies, and the use of routinely collected clinical features for treatment selection could support low-cost deployment in many countries.
Subject(s)
Diabetes Mellitus, Type 2 , Sodium-Glucose Transporter 2 Inhibitors , Male , Humans , Female , Diabetes Mellitus, Type 2/complications , Sodium-Glucose Transporter 2 Inhibitors/therapeutic use , Sodium-Glucose Transporter 2 Inhibitors/pharmacology , Hypoglycemic Agents/adverse effects , Glucagon-Like Peptide-1 Receptor Agonists , Liraglutide/therapeutic use , Bayes Theorem , Glucose , Phenotype , Glucagon-Like Peptide-1 ReceptorABSTRACT
BACKGROUND: Diabetes (regardless of type) and obesity are associated with a range of musculoskeletal disorders. The causal mechanisms driving these associations are unknown for many upper limb pathologies. We used genetic techniques to test the causal link between glycemia, obesity and musculoskeletal conditions. METHODS: In the UK Biobank's unrelated European cohort (N = 379â708) we performed mendelian randomisation (MR) analyses to test for a causal effect of long-term high glycaemia and adiposity on four musculoskeletal pathologies: frozen shoulder, Dupuytren's disease, carpal tunnel syndrome and trigger finger. We also performed single-gene MR using rare variants in the GCK gene. RESULTS: Using MR, we found evidence that long-term high glycaemia has a causal role in the aetiology of upper limb conditions. A 10-mmol/mol increase in genetically predicted haemoglobin A1C (HbA1c) was associated with frozen shoulder: odds ratio (OR) = 1.50 [95% confidence interval (CI), 1.20-1.88], Dupuytren's disease: OR = 1.17 (95% CI, 1.01-1.35), trigger finger: OR = 1.30 (95% CI, 1.09-1.55) and carpal tunnel syndrome: OR = 1.20 (95% CI, 1.09-1.33). Carriers of GCK mutations have increased odds of frozen shoulder: OR = 7.16 (95% CI, 2.93-17.51) and carpal tunnel syndrome: OR = 2.86 (95% CI, 1.50-5.44) but not Dupuytren's disease or trigger finger. We found evidence that an increase in genetically predicted body mass index (BMI) of 5 kg/m2 was associated with carpal tunnel syndrome: OR = 1.13 (95% CI, 1.10-1.16) and associated negatively with Dupuytren's disease: OR = 0.94 (95% CI, 0.90-0.98), but no evidence of association with frozen shoulder or trigger finger. Trigger finger (OR 1.96 (95% CI, 1.42-2.69) P = 3.6e-05) and carpal tunnel syndrome [OR 1.63 (95% CI, 1.36-1.95) P = 8.5e-08] are associated with genetically predicted unfavourable adiposity increase of one standard deviation of body fat. CONCLUSIONS: Our study consistently demonstrates a causal role of long-term high glycaemia in the aetiology of upper limb musculoskeletal conditions. Clinicians treating diabetes patients should be aware of these complications in clinic, specifically those managing the care of GCK mutation carriers. Upper limb musculoskeletal conditions should be considered diabetes complications.
Subject(s)
Bursitis , Carpal Tunnel Syndrome , Diabetes Mellitus , Dupuytren Contracture , Hyperglycemia , Musculoskeletal Diseases , Trigger Finger Disorder , Humans , Dupuytren Contracture/epidemiology , Dupuytren Contracture/genetics , Dupuytren Contracture/complications , Carpal Tunnel Syndrome/epidemiology , Carpal Tunnel Syndrome/genetics , Carpal Tunnel Syndrome/complications , Trigger Finger Disorder/complications , Hyperglycemia/complications , Hyperglycemia/epidemiology , Hyperglycemia/genetics , Upper Extremity , Musculoskeletal Diseases/complications , Risk Factors , Bursitis/complications , Obesity/complications , Obesity/epidemiology , Obesity/geneticsABSTRACT
BACKGROUND: The handling of missing data is a challenge for inference and regression modelling. A particular challenge is dealing with missing predictor information, particularly when trying to build and make predictions from models for use in clinical practice. METHODS: We utilise a flexible Bayesian approach for handling missing predictor information in regression models. This provides practitioners with full posterior predictive distributions for both the missing predictor information (conditional on the observed predictors) and the outcome-of-interest. We apply this approach to a previously proposed counterfactual treatment selection model for type 2 diabetes second-line therapies. Our approach combines a regression model and a Dirichlet process mixture model (DPMM), where the former defines the treatment selection model, and the latter provides a flexible way to model the joint distribution of the predictors. RESULTS: We show that DPMMs can model complex relationships between predictor variables and can provide powerful means of fitting models to incomplete data (under missing-completely-at-random and missing-at-random assumptions). This framework ensures that the posterior distribution for the parameters and the conditional average treatment effect estimates automatically reflect the additional uncertainties associated with missing data due to the hierarchical model structure. We also demonstrate that in the presence of multiple missing predictors, the DPMM model can be used to explore which variable(s), if collected, could provide the most additional information about the likely outcome. CONCLUSIONS: When developing clinical prediction models, DPMMs offer a flexible way to model complex covariate structures and handle missing predictor information. DPMM-based counterfactual prediction models can also provide additional information to support clinical decision-making, including allowing predictions with appropriate uncertainty to be made for individuals with incomplete predictor data.
Subject(s)
Diabetes Mellitus, Type 2 , Humans , Bayes Theorem , Diabetes Mellitus, Type 2/drug therapy , Clinical Decision-Making , UncertaintyABSTRACT
BACKGROUND: Mental health conditions represent one of the major groups of non-transmissible diseases. Physical activity (PA) and sedentary time (ST) have been shown to affect mental health outcomes in opposite directions. In this study, we use accelerometery-derived measures of PA and ST from the UK Biobank (UKB) and depression, anxiety and well-being data from the UKB mental health questionnaire as well as published summary statistics to explore the causal associations between these phenotypes. METHODS: We used MRlap to test if objectively measured PA and ST associate with mental health outcomes using UKB data and summary statistics from published genome-wide association studies. We also tested for bidirectional associations. We performed sex stratified as well as sensitivity analyses. RESULTS: Genetically instrumented higher PA was associated with lower odds of depression (OR = 0.92; 95% CI: 0.88, 0.97) and depression severity (beta = - 0.11; 95% CI: - 0.18, - 0.04), Genetically instrumented higher ST was associated higher odds of anxiety (OR = 2.59; 95% CI: 1.10, 4.60). PA was associated with higher well-being (beta = 0.11, 95% CI: 0.04; 0.18) and ST with lower well-being (beta = - 0.18; 95% CI: - 0.32, - 0.03). Similar findings were observed when stratifying by sex. There was evidence for a bidirectional relationship, with higher genetic liability to depression associated with lower PA (beta = - 0.25, 95% CI: - 0.42; - 0.08) and higher well-being associated with higher PA (beta = 0.15; 95% CI: 0.05, 0.25). CONCLUSIONS: We have demonstrated the bidirectional effects of both PA and ST on a range of mental health outcomes using objectively measured predictors and MR methods for causal inference. Our findings support a causal role for PA and ST in the development of mental health problems and in affecting well-being.
Subject(s)
Depression , Sedentary Behavior , Humans , Depression/epidemiology , Genome-Wide Association Study , Anxiety/epidemiology , Exercise , Mendelian Randomization Analysis/methodsABSTRACT
BACKGROUND: Major depressive disorder (MDD) has a significant impact on global burden of disease. Complications in clinical management can occur when response to pharmacological modalities is considered inadequate and symptoms persist (treatment-resistant depression (TRD)). We aim to investigate inflammation, proxied by C-reactive protein (CRP) levels, and body mass index (BMI) as putative causal risk factors for depression and subsequent treatment resistance, leveraging genetic information to avoid confounding via Mendelian randomisation (MR). METHODS: We used the European UK Biobank subcohort ([Formula: see text]), the mental health questionnaire (MHQ) and clinical records. For treatment resistance, a previously curated phenotype based on general practitioner (GP) records and prescription data was employed. We applied univariable and multivariable MR models to genetically predict the exposures and assess their causal contribution to a range of depression outcomes. We used a range of univariable, multivariable and mediation MR models techniques to address our research question with maximum rigour. In addition, we developed a novel statistical procedure to apply pleiotropy-robust multivariable MR to one sample data and employed a Bayesian bootstrap procedure to accurately quantify estimate uncertainty in mediation analysis which outperforms standard approaches in sparse binary outcomes. Given the flexibility of the one-sample design, we evaluated age and sex as moderators of the effects. RESULTS: In univariable MR models, genetically predicted BMI was positively associated with depression outcomes, including MDD ([Formula: see text] ([Formula: see text] CI): 0.133(0.072, 0.205)) and TRD (0.347(0.002, 0.682)), with a larger magnitude in females and with age acting as a moderator of the effect of BMI on severity of depression (0.22(0.050, 0.389)). Multivariable MR analyses suggested an independent causal effect of BMI on TRD not through CRP (0.395(0.004, 0.732)). Our mediation analyses suggested that the effect of CRP on severity of depression was partly mediated by BMI. Individuals with TRD ([Formula: see text]) observationally had higher CRP and BMI compared with individuals with MDD alone and healthy controls. DISCUSSION: Our work supports the assertion that BMI exerts a causal effect on a range of clinical and questionnaire-based depression phenotypes, with the effect being stronger in females and in younger individuals. We show that this effect is independent of inflammation proxied by CRP levels as the effects of CRP do not persist when jointly estimated with BMI. This is consistent with previous evidence suggesting that overweight contributed to depression even in the absence of any metabolic consequences. It appears that BMI exerts an effect on TRD that persists when we account for BMI influencing MDD.
Subject(s)
Depressive Disorder, Major , Female , Humans , Body Mass Index , Depressive Disorder, Major/epidemiology , Depressive Disorder, Major/genetics , Bayes Theorem , Depression/epidemiology , Depression/genetics , Inflammation/geneticsABSTRACT
Background: The causal relationship between maternal smoking in pregnancy and reduced offspring birth weight is well established and is likely due to impaired placental function. However, observational studies have given conflicting results on the association between smoking and placental weight. We aimed to estimate the causal effect of newly pregnant mothers quitting smoking on their placental weight at the time of delivery. Methods: We used one-sample Mendelian randomization, drawing data from the Avon Longitudinal Study of Parents and Children (ALSPAC) (up to N = 805) and the Norwegian Mother, Father and Child Cohort Study (MoBa) (up to N = 4475). The analysis was performed in pre-pregnancy smokers only, due to the specific role of the genetic instrument SNP rs1051730 (CHRNA5 - CHRNA3 - CHRNB4) in affecting smoking cessation but not initiation. Results: Fixed effect meta-analysis showed a 175 g [95%CI: 16, 334] higher placental weight for pre-pregnancy smoking mothers who continued smoking at the beginning of pregnancy, compared with those who stopped smoking. Using the number of cigarettes smoked per day in the first trimester as the exposure, the causal estimate was a 12 g [95%CI: 2,22] higher placental weight per cigarette per day. Results were similar when the smoking exposures were measured at the end of pregnancy. Using the residuals of birth weight regressed on placental weight as the outcome, we showed weak evidence of lower offspring birth weight relative to the placental weight for continuing smoking. Conclusion: Our results suggest that continued smoking during pregnancy causes higher placental weights.
ABSTRACT
BACKGROUND: The benefit of physical activity (PA) for increasing longevity is well-established, however, the impact of diurnal timing of PA on mortality remains poorly understood. We aimed to derive circadian PA patterns and investigate their associations with all-cause mortality. METHODS: We used 24 h PA time series from 96,351 UK Biobank participants aged between 42 and 79 years at accelerometry in 2013-2015. Functional principal component analysis (fPCA) was applied to obtain circadian PA patterns. Using multivariable Cox proportional hazard models, we related the loading scores of these fPCs to estimate risk of mortality. RESULTS: During 6.9 years of follow-up, 2,850 deaths occurred. Four distinct fPCs accounted for 96% of the variation of the accelerometry data. Using a loading score of zero (i.e., average overall PA during the day) as the reference, a fPC1 score of + 2 (high overall PA) was inversely associated with mortality (Hazard ratio, HR = 0.91; 95% CI: 0.84-0.99), whereas a score of -2 (low overall PA) was associated with higher mortality (1.69; 95% CI: 1.57-1.81; p for non-linearity < 0.001). Significant inverse linear associations with mortality were observed for engaging in midday PA instead of early and late PA (fPC3) (HR for a 1-unit increase 0.88; 95% CI: 0.83-0.93). In contrast, midday and nocturnal PA instead of early and evening PA (fPC4) were positively associated with mortality (HR for a 1-unit increase 1.16; 95% CI: 1.08-1.25). CONCLUSION: Our results suggest that it is less important during which daytime hours one is active but rather, to engage in some level of elevated PA for longevity.
Subject(s)
Accelerometry , Biological Specimen Banks , Humans , Adult , Middle Aged , Aged , Exercise , United KingdomABSTRACT
Previous studies have linked higher body mass index (BMI) to lower subjective well-being in adult European ancestry populations. However, our understanding of these relationships across different populations is limited. Here, we investigated the association between BMI and well-being in people of (a) East Asian and (b) European ancestry in the China Kadoorie Biobank (CKB) and UK Biobank (UKB), respectively. Mendelian randomisation (MR) methods were used to test the relationship between BMI with (a) health satisfaction and (b) life satisfaction. One-sample MR enabled us to test effects in men and women separately and to test the role of cultural contexts by stratifying our analyses by urban and rural home location in both China and the UK. Further, we implemented a control function method to test the linearity of the BMI-well-being relationship. We found evidence of different associations between BMI and well-being in individuals of East Asian versus European ancestry. For example, a genetically instrumented higher BMI tentatively associated with higher health satisfaction in people of East Asian ancestry, especially in females (ß: 0.041, 95% CI: 0.002, 0.081). In contrast, there was a robust inverse association between higher genetically instrumented BMI and health satisfaction in all European ancestry UKB participants (ß: -0.183, 95% CI: -0.200, -0.165, Pdifference < 1.00E-15). We also showed the importance of considering non-linear relationships in the MR framework by providing evidence of non-linear relationships between BMI and health and life satisfaction. Overall, our study suggests potential setting-specific causality in the relationship between BMI and subjective well-being, with robust differences observed between East Asians and Europeans when considering very similar outcomes. We highlight the importance of (a) considering potential non-linear relationships in causal analyses and (b) testing causal relationships in different populations, as the casual nature of relationships, especially relationships influenced by social processes, may be setting-specific.
Subject(s)
Body Mass Index , East Asian People , European People , Health Status , Adult , Female , Humans , Male , ChinaABSTRACT
Telomere attrition is one of biological aging hallmarks and may be intervened to target multiple aging-related diseases, including Alzheimer's disease and Alzheimer's disease related dementias (AD/ADRD). The objective of this study was to assess associations of leukocyte telomere length (TL) with AD/ADRD and early markers of AD/ADRD, including cognitive performance and brain magnetic resonance imaging (MRI) phenotypes. Data from European-ancestry participants in the UK Biobank (n = 435,046) were used to evaluate whether mid-life leukocyte TL is associated with incident AD/ADRD over a mean follow-up of 12.2 years. In a subsample without AD/ADRD and with brain imaging data (n = 43,390), we associated TL with brain MRI phenotypes related to AD or vascular dementia pathology. Longer TL was associated with a lower risk of incident AD/ADRD (adjusted Hazard Ratio [aHR] per SD = 0.93, 95% CI 0.90-0.96, p = 3.37 × 10-7 ). Longer TL also was associated with better cognitive performance in specific cognitive domains, larger hippocampus volume, lower total volume of white matter hyperintensities, and higher fractional anisotropy and lower mean diffusivity in the fornix. In conclusion, longer TL is inversely associated with AD/ADRD, cognitive impairment, and brain structural lesions toward the development of AD/ADRD. However, the relationships between genetically determined TL and the outcomes above were not statistically significant based on the results from Mendelian randomization analysis results. Our findings add to the literature of prioritizing risk for AD/ADRD. The causality needs to be ascertained in mechanistic studies.
Subject(s)
Alzheimer Disease , Humans , Alzheimer Disease/genetics , Alzheimer Disease/pathology , Mendelian Randomization Analysis , Biological Specimen Banks , Leukocytes , Telomere/genetics , Telomere/pathology , United KingdomABSTRACT
Many partial identification problems can be characterized by the optimal value of a function over a set where both the function and set need to be estimated by empirical data. Despite some progress for convex problems, statistical inference in this general setting remains to be developed. To address this, we derive an asymptotically valid confidence interval for the optimal value through an appropriate relaxation of the estimated set. We then apply this general result to the problem of selection bias in population-based cohort studies. We show that existing sensitivity analyses, which are often conservative and difficult to implement, can be formulated in our framework and made significantly more informative via auxiliary information on the population. We conduct a simulation study to evaluate the finite sample performance of our inference procedure, and conclude with a substantive motivating example on the causal effect of education on income in the highly selected UK Biobank cohort. We demonstrate that our method can produce informative bounds using plausible population-level auxiliary constraints. We implement this method in the [Formula: see text] package [Formula: see text].
ABSTRACT
Multivariable Mendelian randomisation (MVMR) is an instrumental variable technique that generalises the MR framework for multiple exposures. Framed as a regression problem, it is subject to the pitfall of multicollinearity. The bias and efficiency of MVMR estimates thus depends heavily on the correlation of exposures. Dimensionality reduction techniques such as principal component analysis (PCA) provide transformations of all the included variables that are effectively uncorrelated. We propose the use of sparse PCA (sPCA) algorithms that create principal components of subsets of the exposures with the aim of providing more interpretable and reliable MR estimates. The approach consists of three steps. We first apply a sparse dimension reduction method and transform the variant-exposure summary statistics to principal components. We then choose a subset of the principal components based on data-driven cutoffs, and estimate their strength as instruments with an adjusted F-statistic. Finally, we perform MR with these transformed exposures. This pipeline is demonstrated in a simulation study of highly correlated exposures and an applied example using summary data from a genome-wide association study of 97 highly correlated lipid metabolites. As a positive control, we tested the causal associations of the transformed exposures on coronary heart disease (CHD). Compared to the conventional inverse-variance weighted MVMR method and a weak instrument robust MVMR method (MR GRAPPLE), sparse component analysis achieved a superior balance of sparsity and biologically insightful grouping of the lipid traits.
Subject(s)
Coronary Disease , Genome-Wide Association Study , Humans , Mendelian Randomization Analysis/methods , Causality , LipidsABSTRACT
BACKGROUND: Extensive evidence links higher body mass index (BMI) to higher odds of depression in people of European ancestry. However, our understanding of the relationship across different settings and ancestries is limited. Here, we test the relationship between body composition and depression in people of East Asian ancestry. METHODS: Multiple Mendelian randomisation (MR) methods were used to test the relationship between (a) BMI and (b) waist-hip ratio (WHR) with depression. Firstly, we performed two-sample MR using genetic summary statistics from a recent genome-wide association study (GWAS) of depression (with 15,771 cases and 178,777 controls) in people of East Asian ancestry. We selected 838 single nucleotide polymorphisms (SNPs) correlated with BMI and 263 SNPs correlated with WHR as genetic instrumental variables to estimate the causal effect of BMI and WHR on depression using the inverse-variance weighted (IVW) method. We repeated these analyses stratifying by home location status: China versus UK or USA. Secondly, we performed one-sample MR in the China Kadoorie Biobank (CKB) in 100,377 participants. This allowed us to test the relationship separately in (a) males and females and (b) urban and rural dwellers. We also examined (c) the linearity of the BMI-depression relationship. RESULTS: Both MR analyses provided evidence that higher BMI was associated with lower odds of depression. For example, a genetically-instrumented 1-SD higher BMI in the CKB was associated with lower odds of depressive symptoms [OR: 0.77, 95% CI: 0.63, 0.95]. There was evidence of differences according to place of residence. Using the IVW method, higher BMI was associated with lower odds of depression in people of East Asian ancestry living in China but there was no evidence for an association in people of East Asian ancestry living in the USA or UK. Furthermore, higher genetic BMI was associated with differential effects in urban and rural dwellers within China. CONCLUSIONS: This study provides the first MR evidence for an inverse relationship between BMI and depression in people of East Asian ancestry. This contrasts with previous findings in European populations and therefore the public health response to obesity and depression is likely to need to differ based on sociocultural factors for example, ancestry and place of residence. This highlights the importance of setting-specific causality when using genetic causal inference approaches and data from diverse populations to test hypotheses. This is especially important when the relationship tested is not purely biological and may involve sociocultural factors.
Subject(s)
Body Composition , Depression , East Asian People , Genome-Wide Association Study , Female , Humans , Male , Body Composition/genetics , Body Mass Index , Depression/epidemiology , Depression/genetics , Mendelian Randomization Analysis , Obesity/genetics , Polymorphism, Single Nucleotide/genetics , ChinaABSTRACT
BACKGROUND: Mendelian randomization (MR) leverages genetic data as an instrumental variable to provide estimates for the causal effect of an exposure X on a health outcome Y that is robust to confounding. Unfortunately, horizontal pleiotropy-the direct association of a genetic variant with multiple phenotypes-is highly prevalent and can easily render a genetic variant an invalid instrument. METHODS: Building on existing work, we propose a simple method for leveraging sex-specific genetic associations to perform weak and pleiotropy-robust MR analysis. This is achieved by constructing an MR estimator in which pleiotropy is perfectly removed by cancellation, while placing it within the powerful machinery of the robust adjusted profile score (MR-RAPS) method. Pleiotropy cancellation has the attractive property that it removes heterogeneity and therefore justifies a statistically efficient fixed effects model. We extend the method from the typical two-sample summary-data MR setting to the one-sample setting by adapting the technique of Collider-Correction. Simulation studies and applied examples are used to assess how the sex-stratified MR-RAPS estimator performs against other common approaches. RESULTS: The sex-stratified MR-RAPS method is shown to be robust to pleiotropy even in cases where all genetic variants violated the standard Instrument Strength Independent of Direct Effect assumption. In some cases where the strength of the pleiotropic effect additionally varied by sex (and so perfect cancellation was not achieved), over-dispersed MR-RAPS implementations can still consistently estimate the true causal effect. In applied analyses, we investigate the causal effect of waist-hip ratio (WHR), an important marker of central obesity, on a range of downstream traits. While the conventional approaches suggested paradoxical links between WHR and height and body mass index, the sex-stratified approach obtained a more realistic null effect. Nonzero effects were also detected for systolic and diastolic blood pressure as well as high-density and low-density lipoprotein cholesterol. DISCUSSION: We provide a simple but attractive method for weak and pleiotropy robust causal estimation of sexually dimorphic traits on downstream outcomes, by combining several existing approaches in a novel fashion.