Simultaneous estimation of bi-directional causal effects and heritable confounding from GWAS summary statistics

Identification and correction for collider bias in a genome-wide association study of diabetes-related heart failure

Type 2 diabetes (T2D) is a major risk factor for heart failure (HF) and has elevated incidence among individuals with HF. Since genetics and HF can independently influence T2D, collider bias may occur when T2D (i.e., collider) is controlled for by design or analysis. Thus, we conducted a genome-wide association study (GWAS) of diabetes-related HF with correction for collider bias. We first performed a GWAS of HF to identify genetic instrumental variables (GIVs) for HF and to enable bidirectional Mendelian randomization (MR) analysis between T2D and HF. We identified 61 genomic loci, significantly associated with all-cause HF in 114,275 individuals with HF and over 1.5 million controls of European ancestry. Using a two-sample bidirectional MR approach with 59 and 82 GIVs for HF and T2D, respectively, we estimated that T2D increased HF risk (odds ratio [OR] 1.07, 95% confidence interval [CI] 1.04-1.10), while HF also increased T2D risk (OR 1.60, 95% CI 1.36-1.88). Then we performed a GWAS of diabetes-related HF corrected for collider bias due to the study design of index cases. After removing the spurious association of TCF7L2 locus due to collider bias, we identified two genome-wide significant loci close to PITX2 (chromosome 4) and CDKN2B-AS1 (chromosome 9) associated with diabetes-related HF in the Million Veteran Program and replicated the associations in the UK Biobank. Our MR findings provide strong evidence that HF increases T2D risk. As a result, collider bias leads to spurious genetic associations of diabetes-related HF, which can be effectively corrected to identify true positive loci.

Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

Mendelian randomization: causal inference leveraging genetic data

Mendelian randomization (MR) leverages genetic information to examine the causal relationship between phenotypes allowing for the presence of unmeasured confounders. MR has been widely applied to unresolved questions in epidemiology, making use of summary statistics from genome-wide association studies on an increasing number of human traits. However, an understanding of essential concepts is necessary for the appropriate application and interpretation of MR. This review aims to provide a non-technical overview of MR and demonstrate its relevance to psychiatric research. We begin with the origins of MR and the reasons for its recent expansion, followed by an overview of its statistical methodology. We then describe the limitations of MR, and how these are being addressed by recent methodological advances. We showcase the practical use of MR in psychiatry through three illustrative examples - the connection between cannabis use and psychosis, the link between intelligence and schizophrenia, and the search for modifiable risk factors for depression. The review concludes with a discussion of the prospects of MR, focusing on the integration of multi-omics data and its extension to delineating complex causal networks.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

Estimating causality between smoking and abdominal obesity by Mendelian randomization

BACKGROUND AND AIMS

Smokers tend to have a lower body weight than non-smokers, but also more abdominal fat. It remains unclear whether or not the relationship between smoking and abdominal obesity is causal. Previous Mendelian randomization (MR) studies have investigated this relationship by relying upon a single genetic variant for smoking heaviness. This approach is sensitive to pleiotropic effects and may produce imprecise causal estimates. We aimed to estimate causality between smoking and abdominal obesity using multiple genetic instruments.

DESIGN

MR study using causal analysis using summary effect estimates (CAUSE) and latent heritable confounder MR (LHC-MR) methods that instrument smoking using genome-wide data, and also two-sample MR (2SMR) methods.

SETTING

Genome-wide association studies (GWAS) summary statistics from participants of European ancestry, obtained from the GWAS and Sequencing Consortium of Alcohol and Nicotine use (GSCAN), Genetic Investigation of Anthropometric Traits (GIANT) Consortium and the UK Biobank.

PARTICIPANTS

We used GWAS results for smoking initiation (n = 1 232 091), life-time smoking (n = 462 690) and smoking heaviness (n = 337 334) as exposure traits, and waist-hip ratio (WHR) and waist and hip circumferences (WC and HC) (n up to 697 734), with and without adjustment for body mass index (adjBMI), as outcome traits.

MEASUREMENTS

Smoking initiation, life-time smoking, smoking heaviness, WHR, WC, HC, WHRadjBMI, WCadjBMI and HCadjBMI.

FINDINGS

Both CAUSE and LHC-MR indicated a positive causal effect of smoking initiation on WHR (0.13 [95% confidence interval (CI) = 0.10, 0.16 and 0.49 (0.41, 0.57), respectively] and WHRadjBMI (0.07 (0.03, 0.10) and 0.31 (0.26, 0.37). Similarly, they indicated a positive causal effect of life-time smoking on WHR [0.35 (0.29, 0.41) and 0.44 (0.38, 0.51)] and WHRadjBMI [0.18 (0.13, 0.24) and 0.26 (0.20, 0.31)]. In follow-up analyses, smoking particularly increased visceral fat. There was no evidence of a mediating role by cortisol or sex hormones.

CONCLUSIONS

Smoking initiation and higher life-time smoking may lead to increased abdominal fat distribution. The increase in abdominal fat due to smoking is characterized by an increase in visceral fat. Thus, efforts to prevent and cease smoking can have the added benefit of reducing abdominal fat.

Breaking down causes, consequences, and mediating effects of telomere length variation on human health

BACKGROUND

Telomeres form repeated DNA sequences at the ends of chromosomes, which shorten with each cell division. Yet, factors modulating telomere attrition and the health consequences thereof are not fully understood. To address this, we leveraged data from 326,363 unrelated UK Biobank participants of European ancestry.

RESULTS

Using linear regression and bidirectional univariable and multivariable Mendelian randomization (MR), we elucidate the relationships between leukocyte telomere length (LTL) and 142 complex traits, including diseases, biomarkers, and lifestyle factors. We confirm that telomeres shorten with age and show a stronger decline in males than in females, with these factors contributing to the majority of the 5.4% of LTL variance explained by the phenome. MR reveals 23 traits modulating LTL. Smoking cessation and high educational attainment associate with longer LTL, while weekly alcohol intake, body mass index, urate levels, and female reproductive events, such as childbirth, associate with shorter LTL. We also identify 24 traits affected by LTL, with risk for cardiovascular, pulmonary, and some autoimmune diseases being increased by short LTL, while longer LTL increased risk for other autoimmune conditions and cancers. Through multivariable MR, we show that LTL may partially mediate the impact of educational attainment, body mass index, and female age at childbirth on proxied lifespan.

CONCLUSIONS

Our study sheds light on the modulators, consequences, and the mediatory role of telomeres, portraying an intricate relationship between LTL, diseases, lifestyle, and socio-economic factors.

Unravelling the complex causal effects of substance use behaviours on common diseases

BACKGROUND

Substance use behaviours (SUB) including smoking, alcohol consumption, and coffee intake are associated with many health outcomes. However, whether the health effects of SUB are causal remains controversial, especially for alcohol consumption and coffee intake.

METHODS

In this study, we assess 11 commonly used Mendelian Randomization (MR) methods by simulation and apply them to investigate the causal relationship between 7 SUB traits and health outcomes. We also combine stratified regression, genetic correlation, and MR analyses to investigate the dosage-dependent effects.

RESULTS

We show that smoking initiation has widespread risk effects on common diseases such as asthma, type 2 diabetes, and peripheral vascular disease. Alcohol consumption shows risk effects specifically on cardiovascular diseases, dyslipidemia, and hypertensive diseases. We find evidence of dosage-dependent effects of coffee and tea intake on common diseases (e.g., cardiovascular disease and osteoarthritis). We observe that the minor allele effect of rs4410790 (the top signal for tea intake level) is negative on heavy tea intake ( b ̂ G W A S = - 0.091 , s . e . = 0.007 , P = 4.90 × 10 - 35 ) but positive on moderate tea intake ( b ̂ G W A S = 0.034 , s . e . = 0.006 , P = 3.40 × 10 - 8 ) , compared to the non-tea-drinkers.

CONCLUSION

Our study reveals the complexity of the health effects of SUB and informs design for future studies aiming to dissect the causal relationships between behavioural traits and complex diseases.

Approaches to estimate bidirectional causal effects using Mendelian randomization with application to body mass index and fasting glucose

Mendelian randomization (MR) is an epidemiological framework using genetic variants as instrumental variables (IVs) to examine the causal effect of exposures on outcomes. Statistical methods based on unidirectional MR (UMR) are widely used to estimate the causal effects of exposures on outcomes in observational studies. To estimate the bidirectional causal effects between two phenotypes, investigators have naively applied UMR methods separately in each direction. However, bidirectional causal effects between two phenotypes create a feedback loop that biases the estimation when UMR methods are naively applied. To overcome this limitation, we proposed two novel approaches to estimate bidirectional causal effects using MR: BiRatio and BiLIML, which are extensions of the standard ratio, and limited information maximum likelihood (LIML) methods, respectively. We compared the performance of the two proposed methods with the naive application of UMR methods through extensive simulations of several scenarios involving varying numbers of strong and weak IVs. Our simulation results showed that when multiple strong IVs are used, the proposed methods provided accurate bidirectional causal effect estimation in terms of median absolute bias and relative median absolute bias. Furthermore, compared to the BiRatio method, the BiLIML method provided a more accurate estimation of causal effects when weak IVs were used. Therefore, based on our simulations, we concluded that the BiLIML should be used for bidirectional causal effect estimation. We applied the proposed methods to investigate the potential bidirectional relationship between obesity and diabetes using the data from the Multi-Ethnic Study of Atherosclerosis cohort. We used body mass index (BMI) and fasting glucose (FG) as measures of obesity and type 2 diabetes, respectively. Our results from the BiLIML method revealed the bidirectional causal relationship between BMI and FG in across all racial populations. Specifically, in the White/Caucasian population, a 1 kg/m2 increase in BMI increased FG by 0.70 mg/dL (95% confidence interval [CI]: 0.3517-1.0489; p = 8.43×10-5), and 1 mg/dL increase in FG increased BMI by 0.10 kg/m2 (95% CI: 0.0441-0.1640; p = 6.79×10-4). Our study provides novel findings and quantifies the effect sizes of the bidirectional causal relationship between BMI and FG. However, further studies are needed to understand the biological and functional mechanisms underlying the bidirectional pathway.

Bidirectional relationship between type 2 diabetes mellitus and coronary artery disease: Prospective cohort study and genetic analyses

BACKGROUND

While type 2 diabetes mellitus (T2DM) is considered a putative causal risk factor for coronary artery disease (CAD), the intrinsic link underlying T2DM and CAD is not fully understood. We aimed to highlight the importance of integrated care targeting both diseases by investigating the phenotypic and genetic relationships between T2DM and CAD.

METHODS

We evaluated phenotypic associations using data from the United Kingdom Biobank ( N = 472,050). We investigated genetic relationships by leveraging genomic data conducted in European ancestry for T2DM, with and without adjustment for body mass index (BMI) (T2DM: Ncase / Ncontrol = 74,124/824,006; T2DM adjusted for BMI [T2DM adj BMI]: Ncase / Ncontrol = 50,409/523,897) and for CAD ( Ncase / Ncontrol = 181,522/984,168). We performed additional analyses using genomic data conducted in multiancestry individuals for T2DM ( Ncase / Ncontrol = 180,834/1,159,055).

RESULTS

Observational analysis suggested a bidirectional relationship between T2DM and CAD (T2DM→CAD: hazard ratio [HR] = 2.12, 95% confidence interval [CI]: 2.01-2.24; CAD→T2DM: HR = 1.72, 95% CI: 1.63-1.81). A positive overall genetic correlation between T2DM and CAD was observed ( rg = 0.39, P = 1.43 × 10 -75 ), which was largely independent of BMI (T2DM adj BMI-CAD: rg = 0.31, P = 1.20 × 10 -36 ). This was corroborated by six local signals, among which 9p21.3 showed the strongest genetic correlation. Cross-trait meta-analysis replicated 101 previously reported loci and discovered six novel pleiotropic loci. Mendelian randomization analysis supported a bidirectional causal relationship (T2DM→CAD: odds ratio [OR] = 1.13, 95% CI: 1.11-1.16; CAD→T2DM: OR = 1.12, 95% CI: 1.07-1.18), which was confirmed in multiancestry individuals (T2DM→CAD: OR = 1.13, 95% CI: 1.10-1.16; CAD→T2DM: OR = 1.08, 95% CI: 1.04-1.13). This bidirectional relationship was significantly mediated by systolic blood pressure and intake of 3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitors, with mediation proportions of 54.1% (95% CI: 24.9-83.4%) and 90.4% (95% CI: 29.3-151.5%), respectively.

CONCLUSION

Our observational and genetic analyses demonstrated an intrinsic bidirectional relationship between T2DM and CAD and clarified the biological mechanisms underlying this relationship.

PheWAS-based clustering of Mendelian Randomisation instruments reveals distinct mechanism-specific causal effects between obesity and educational attainment

Mendelian Randomisation (MR) estimates causal effects between risk factors and complex outcomes using genetic instruments. Pleiotropy, heritable confounders, and heterogeneous causal effects violate MR assumptions and can lead to biases. To alleviate these, we propose an approach employing a Phenome-Wide association Clustering of the MR instruments (PWC-MR) and apply this method to revisit the surprisingly large apparent causal effect of body mass index (BMI) on educational attainment (EDU): [Formula: see text] = -0.19 [-0.22, -0.16]. First, we cluster 324 BMI-associated genetic instruments based on their association with 407 traits in the UK Biobank, which yields six distinct groups. Subsequent cluster-specific MR reveals heterogeneous causal effect estimates on EDU. A cluster enriched for socio-economic indicators yields the largest BMI-on-EDU causal effect estimate ([Formula: see text] = -0.49 [-0.56, -0.42]) whereas a cluster enriched for body-mass specific traits provides a more likely estimate ([Formula: see text] = -0.09 [-0.13, -0.05]). Follow-up analyses confirms these findings: within-sibling MR ([Formula: see text] = -0.05 [-0.09, -0.01]); MR for childhood BMI on EDU ([Formula: see text] = -0.03 [-0.06, -0.002]); step-wise multivariable MR ([Formula: see text] = -0.05 [-0.07, -0.02]) where socio-economic indicators are jointly modelled. Here we show how the in-depth examination of the BMI-EDU causal relationship demonstrates the utility of our PWC-MR approach in revealing distinct pleiotropic pathways and confounder mechanisms.

MRSL: a causal network pruning algorithm based on GWAS summary data

Causal discovery is a powerful tool to disclose underlying structures by analyzing purely observational data. Genetic variants can provide useful complementary information for structure learning. Recently, Mendelian randomization (MR) studies have provided abundant marginal causal relationships of traits. Here, we propose a causal network pruning algorithm MRSL (MR-based structure learning algorithm) based on these marginal causal relationships. MRSL combines the graph theory with multivariable MR to learn the conditional causal structure using only genome-wide association analyses (GWAS) summary statistics. Specifically, MRSL utilizes topological sorting to improve the precision of structure learning. It proposes MR-separation instead of d-separation and three candidates of sufficient separating set for MR-separation. The results of simulations revealed that MRSL had up to 2-fold higher F1 score and 100 times faster computing time than other eight competitive methods. Furthermore, we applied MRSL to 26 biomarkers and 44 International Classification of Diseases 10 (ICD10)-defined diseases using GWAS summary data from UK Biobank. The results cover most of the expected causal links that have biological interpretations and several new links supported by clinical case reports or previous observational literatures.

Subset scanning for multi-trait analysis using GWAS summary statistics

MOTIVATION

Multi-trait analysis has been shown to have greater statistical power than single-trait analysis. Most of the existing multi-trait analysis methods only work with a limited number of traits and usually prioritize high statistical power over identifying relevant traits, which heavily rely on domain knowledge.

RESULTS

To handle diseases and traits with obscure etiology, we developed TraitScan, a powerful and fast algorithm that identifies potential pleiotropic traits from a moderate or large number of traits (e.g. dozens to thousands) and tests the association between one genetic variant and the selected traits. TraitScan can handle either individual-level or summary-level GWAS data. We evaluated TraitScan using extensive simulations and found that it outperformed existing methods in terms of both testing power and trait selection when sparsity was low or modest. We then applied it to search for traits associated with Ewing Sarcoma, a rare bone tumor with peak onset in adolescence, among 754 traits in UK Biobank. Our analysis revealed a few promising traits worthy of further investigation, highlighting the use of TraitScan for more effective multi-trait analysis as biobanks emerge. We also extended TraitScan to search and test association with a polygenic risk score and genetically imputed gene expression.

AVAILABILITY AND IMPLEMENTATION

Our algorithm is implemented in an R package "TraitScan" available at https://github.com/RuiCao34/TraitScan.

Investigating the Shared Genetic Etiology Between Parkinson's Disease and Depression

Background

Depression is a common symptom in Parkinson's disease (PD), resulting from underlying neuropathological processes and psychological factors. However, the extent to which shared genetic risk factors contribute to the relationship between depression and PD is poorly understood.

Objective

To examine the effects of common genetic variants influencing the etiology of PD and depression risk at the genome-wide and local genomic regional level.

Methods

We comprehensively investigated the genetic relationship between PD and depression using genome-wide association studies data. First, we estimated the genetic correlation at the genome-wide level using linkage-disequilibrium score regression, followed by local genetic correlation analysis using the GWAS-pairwise method and functional annotation to identify genes that may jointly influence the risk for both traits. Also, we performed Latent Causal Variable, Latent Heritable Confounder Mendelian Randomization, and traditional Mendelian Randomization analyses to investigate the potential causal relationship.

Results

Although the genetic correlation between PD and depression was not statistically significant at the genome-wide level, GWAS-pairwise analyses identified 16 genomic segments associated with PD and depression, implicating nine genes. Further analyses revealed distinct patterns within individual genes, suggesting an intricate pattern. These genes involve various biological processes, including neurotransmitter regulation, senescence, and nucleo-cytoplasmic transport mechanisms. We did not observe genetic evidence of causality between PD and depression.

Conclusions

Our findings did not support a genome-wide genetic correlation or a causal association between both conditions. However, we identified genomic segments but identified genomic segments linked to distinct biological pathways influencing their etiology.Further research is needed to understand their functional consequences.

Distilling causality between physical activity traits and obesity via Mendelian randomization

BACKGROUND

Whether obesity is a cause or consequence of low physical activity levels and more sedentary time has not yet been fully elucidated. Better instrumental variables and a more thorough consideration of potential confounding variables that may influence the causal inference between physical activity and obesity are needed.

METHODS

Leveraging results from our recent genome-wide association study for leisure time moderate-to-vigorous intensity (MV) physical activity and screen time, we here disentangle the causal relationships between physical activity, sedentary behavior, education-defined by years of schooling-and body mass index (BMI), using multiple univariable and multivariable Mendelian Randomization (MR) approaches.

RESULTS

Univariable MR analyses suggest bidirectional causal effects of physical activity and sedentary behavior with BMI. However, multivariable MR analyses that take years of schooling into account suggest that more MV physical activity causes a lower BMI, and a higher BMI causes more screen time, but not vice versa. In addition, more years of schooling causes higher levels of MV physical activity, less screen time, and lower BMI.

CONCLUSIONS

In conclusion, our results highlight the beneficial effect of education on improved health and suggest that a more physically active lifestyle leads to lower BMI, while sedentary behavior is a consequence of higher BMI.

Loneliness and depression: bidirectional mendelian randomization analyses using data from three large genome-wide association studies

Major depression (MD) is a serious psychiatric illness afflicting nearly 5% of the world's population. A large correlational literature suggests that loneliness is a prospective risk factor for MD; correlational assocations of this nature may be confounded for a variety of reasons. This report uses Mendelian Randomization (MR) to examine potentially causal associations between loneliness and MD. We report on analyses using summary statistics from three large genome wide association studies (GWAS). MR analyses were conducted using three independent sources of GWAS summary statistics. In the first set of analyses, we used available summary statistics from an extant GWAS of loneliness to predict MD risk. We used two sources of outcome data: the Psychiatric Genomics Consortium (PGC) meta-analysis of MD (PGC-MD; N = 142,646) and the Million Veteran Program (MVP-MD; N = 250,215). Finally, we reversed analyses using data from the MVP and PGC samples to identify risk variants for MD and used loneliness outcome data from UK Biobank. We find robust evidence for a bidirectional causal relationship between loneliness and MD, including between loneliness, depression cases status, and a continuous measure of depressive symptoms. The estimates remained significant across several sensitivity analyses, including models that account for horizontal pleiotropy. This paper provides the first genetically-informed evidence that reducing loneliness may play a causal role in decreasing risk for depressive illness, and these findings support efforts to reduce loneliness in order to prevent or ameliorate MD. Discussion focuses on the public health significance of these findings, especially in light of the SARS-CoV-2 pandemic.

Deciphering genetic causality between inflammatory bowel disease and periodontitis through bi-directional two-sample Mendelian randomization

Inflammatory bowel disease (IBD) and periodontitis are reported to be closely associated; however, whether there is a causal association between them remains unclear. To explore the existence of this causality, this study applied a bidirectional two-sample Mendelian randomization (MR). The genetic variants were obtained from the summary statistics of genome-wide association studies of IBD, including its subtypes CD and UC, and periodontitis. 175, 148, 113, and six single-nucleotide polymorphisms were selected as instrumental variables for IBD, CD, UC, and periodontitis, respectively. In MR analysis, random-effects inverse-variance weighted was used as the primary method, and weighted median and MR Egger regression were applied as the complementary method. A series of sensitivity analyses were also conducted to ensure the reliability of the results. None of these analyses found a significant effect of genetically proxied IBD and its subtypes on periodontitis, and vice versa. Subsequent sensitivity analyses did not detect any horizontal pleiotropy and heterogeneity. Caution should be exerted when it comes to clinical relevance and further studies are needed to clarify the relationship between IBD and periodontitis.

Correction for Collider Bias in the Genome-wide Association Study of Diabetes-Related Heart Failure due to Bidirectional Relationship between Heart Failure and Type 2 Diabetes

Aims

Type 2 diabetes (T2D) is a major risk factor for heart failure (HF) across demographic groups. On the other hand, metabolic impairment, including elevated T2D incidence is a hallmark of HF pathophysiology. We investigated the bidirectional relationship between T2D and HF, and identified genetic associations with diabetes-related HF after correction for potential collider bias.

Methods

We performed a genome-wide association study (GWAS) of HF to identify genetic instrumental variables (GIVs) for HF, and to enable bidirectional Mendelian Randomization (MR) analysis between T2D and HF. Since genetics and HF can independently influence T2D, collider bias may occur when T2D (i.e., collider) is controlled for by design or analysis. Thus, we conducted GWAS of diabetes-related HF with correction for collider bias.

Results

We first identified 61 genomic loci, including 24 novel loci, significantly associated with all-cause HF in 114,275 HF cases and over 1.5 million controls of European ancestry. Combined with the summary statistics of a T2D GWAS, we obtained 59 and 82 GIVs for HF and T2D, respectively. Using a two-sample bidirectional MR approach, we estimated that T2D increased HF risk (OR 1.07, 95% CI 1.04-1.10), while HF also increased T2D risk (OR 1.60, 95% CI 1.36-1.88). Then we performed a GWAS of diabetes-related HF corrected for collider bias due to prevalent HF affecting incidence of T2D. After removing the spurious association of TCF7L2 locus due to collider bias, we identified two genome-wide significant loci close to PITX2 (chromosome 4) and CDKN2B-AS1 (chromosome 9) associated with diabetes-related HF in the Million Veteran Program, and replicated the associations in the UK Biobank study.

Conclusion

We identified novel HF-associated loci to enable bidirectional MR study of T2D and HF. Our MR findings support T2D as a HF risk factor and provide strong evidence that HF increases T2D risk. As a result, collider bias leads to spurious genetic associations of diabetes-related HF, which can be effectively corrected to identify true positive loci. Evaluation of collider bias should be a critical component when conducting GWAS of complex disease phenotypes such as diabetes-related cardiovascular complications.

Guidelines for performing Mendelian randomization investigations: update for summer 2023

This paper provides guidelines for performing Mendelian randomization investigations. It is aimed at practitioners seeking to undertake analyses and write up their findings, and at journal editors and reviewers seeking to assess Mendelian randomization manuscripts. The guidelines are divided into ten sections: motivation and scope, data sources, choice of genetic variants, variant harmonization, primary analysis, supplementary and sensitivity analyses (one section on robust statistical methods and one on other approaches), extensions and additional analyses, data presentation, and interpretation. These guidelines will be updated based on feedback from the community and advances in the field. Updates will be made periodically as needed, and at least every 24 months.

Understanding the relationship between loneliness, substance use traits and psychiatric disorders: A genetically informed approach

Loneliness is a common, yet distressing experience associated with adverse outcomes including substance use problems and psychiatric disorders. To what extent these associations reflect genetic correlations and causal relationships is currently unclear. We applied Genomic Structural Equation Modelling (GSEM) to dissect the genetic architecture between loneliness and psychiatric-behavioural traits. Included were summary statistics from 12 genome-wide association analyses, including loneliness and 11 psychiatric phenotypes (range N: 9,537 - 807,553). We first modelled latent genetic factors amongst the psychiatric traits to then investigate potential causal effects between loneliness and the identified latent factors, using multivariate genome-wide association analyses and bidirectional Mendelian randomization. We identified three latent genetic factors, encompassing neurodevelopmental/mood conditions, substance use traits and disorders with psychotic features. GSEM provided evidence of a unique association between loneliness and the neurodevelopmental/mood conditions latent factor. Mendelian randomization results were suggestive of bidirectional causal effects between loneliness and the neurodevelopmental/mood conditions factor. These results imply that a genetic predisposition to loneliness may elevate the risk of neurodevelopmental/mood conditions, and vice versa. However, results may reflect the difficulty of distiguishing between loneliness and neurodevelopmental/mood conditions, which present in similar ways. We suggest, overall, the importance of addressing loneliness in mental health prevention and policy.

Partner choice, confounding and trait convergence all contribute to phenotypic partner similarity

Partners are often similar in terms of their physical and behavioural traits, such as their education, political affiliation and height. However, it is currently unclear what exactly causes this similarity-partner choice, partner influence increasing similarity over time or confounding factors such as shared environment or indirect assortment. Here, we applied Mendelian randomization to the data of 51,664 couples in the UK Biobank and investigated partner similarity in 118 traits. We found evidence of partner choice for 64 traits, 40 of which had larger phenotypic correlation than causal effect. This suggests that confounders contribute to trait similarity, among which household income, overall health rating and education accounted for 29.8, 14.1 and 11.6% of correlations between partners, respectively. Finally, mediation analysis revealed that most causal associations between different traits in the two partners are indirect. In summary, our results show the mechanisms through which indirect assortment increases the observed partner similarity.

Genetic insights into the causal relationship between physical activity and cognitive functioning

Physical activity and cognitive functioning are strongly intertwined. However, the causal relationships underlying this association are still unclear. Physical activity can enhance brain functions, but healthy cognition may also promote engagement in physical activity. Here, we assessed the bidirectional relationships between physical activity and general cognitive functioning using Latent Heritable Confounder Mendelian Randomization (LHC-MR). Association data were drawn from two large-scale genome-wide association studies (UK Biobank and COGENT) on accelerometer-measured moderate, vigorous, and average physical activity (N = 91,084) and cognitive functioning (N = 257,841). After Bonferroni correction, we observed significant LHC-MR associations suggesting that increased fraction of both moderate (b = 0.32, CI95% = [0.17,0.47], P = 2.89e - 05) and vigorous physical activity (b = 0.22, CI95% = [0.06,0.37], P = 0.007) lead to increased cognitive functioning. In contrast, we found no evidence of a causal effect of average physical activity on cognitive functioning, and no evidence of a reverse causal effect (cognitive functioning on any physical activity measures). These findings provide new evidence supporting a beneficial role of moderate and vigorous physical activity (MVPA) on cognitive functioning.

Shared genetic architecture between attention-deficit/hyperactivity disorder and lifespan

There is evidence linking ADHD to a reduced life expectancy. The mortality rate in individuals with ADHD is twice that of the general population and it is associated with several factors, such as unhealthy lifestyle behaviors, social adversity, and mental health problems that may in turn increase mortality rates. Since ADHD and lifespan are heritable, we used data from genome-wide association studies (GWAS) of ADHD and parental lifespan, as proxy of individual lifespan, to estimate their genetic correlation, identify genetic loci jointly associated with both phenotypes and assess causality. We confirmed a negative genetic correlation between ADHD and parental lifespan (rg = -0.36, P = 1.41e-16). Nineteen independent loci were jointly associated with both ADHD and parental lifespan, with most of the alleles that increased the risk for ADHD being associated with shorter lifespan. Fifteen loci were novel for ADHD and two were already present in the original GWAS on parental lifespan. Mendelian randomization analyses pointed towards a negative causal effect of ADHD liability on lifespan (P = 1.54e-06; Beta = -0.07), although these results were not confirmed by all sensitivity analyses performed, and further evidence is required. The present study provides the first evidence of a common genetic background between ADHD and lifespan, which may play a role in the reported effect of ADHD on premature mortality risk. These results are consistent with previous epidemiological data describing reduced lifespan in mental disorders and support that ADHD is an important health condition that could negatively affect future life outcomes.

Role of the Gut-Brain Axis in the Shared Genetic Etiology Between Gastrointestinal Tract Diseases and Psychiatric Disorders: A Genome-Wide Pleiotropic Analysis

IMPORTANCE

Comorbidities and genetic correlations between gastrointestinal tract diseases and psychiatric disorders have been widely reported, with the gut-brain axis (GBA) hypothesized as a potential biological basis. However, the degree to which the shared genetic determinants are involved in these associations underlying the GBA is unclear.

OBJECTIVE

To investigate the shared genetic etiology between gastrointestinal tract diseases and psychiatric disorders and to identify shared genomic loci, genes, and pathways.

DESIGN, SETTING, AND PARTICIPANTS

This genome-wide pleiotropic association study using genome-wide association summary statistics from publicly available data sources was performed with various statistical genetic approaches to sequentially investigate the pleiotropic associations from genome-wide single-nucleotide variation (SNV; formerly single-nucleotide polymorphism [SNP]), and gene levels and biological pathways to disentangle the underlying shared genetic etiology between 4 gastrointestinal tract diseases (inflammatory bowel disease, irritable bowel syndrome, peptic ulcer disease, and gastroesophageal reflux disease) and 6 psychiatric disorders (schizophrenia, bipolar disorder, major depressive disorder, attention-deficit/hyperactivity disorder, posttraumatic stress disorder, and anorexia nervosa). Data were collected from March 10, 2021, to August 25, 2021, and analysis was performed from January 8 through May 30, 2022.

MAIN OUTCOMES AND MEASURES

The primary outcomes consisted of a list of genetic loci, genes, and pathways shared between gastrointestinal tract diseases and psychiatric disorders.

RESULTS

Extensive genetic correlations and genetic overlaps were found among 22 of 24 trait pairs. Pleiotropic analysis under a composite null hypothesis identified 2910 significant potential pleiotropic SNVs in 19 trait pairs, with 83 pleiotropic loci and 24 colocalized loci detected. Gene-based analysis found 158 unique candidate pleiotropic genes, which were highly enriched in certain GBA-related phenotypes and tissues, whereas pathway enrichment analysis further highlighted biological pathways primarily involving cell adhesion, synaptic structure and function, and immune cell differentiation. Several identified pleiotropic loci also shared causal variants with gut microbiomes. Mendelian randomization analysis further illustrated vertical pleiotropy across 8 pairwise traits. Notably, many pleiotropic loci were identified for multiple pairwise traits, such as 1q32.1 (INAVA), 19q13.33 (FUT2), 11q23.2 (NCAM1), and 1p32.3 (LRP8).

CONCLUSIONS AND RELEVANCE

These findings suggest that the pleiotropic genetic determinants between gastrointestinal tract diseases and psychiatric disorders are extensively distributed across the genome. These findings not only support the shared genetic basis underlying the GBA but also have important implications for intervention and treatment targets of these diseases simultaneously.

A Guide for Understanding and Designing Mendelian Randomization Studies in the Musculoskeletal Field

Mendelian randomization (MR) is an increasingly popular component of an epidemiologist's toolkit, used to provide evidence of a causal effect of one trait (an exposure, eg, body mass index [BMI]) on an outcome trait or disease (eg, osteoarthritis). Identifying these effects is important for understanding disease etiology and potentially identifying targets for therapeutic intervention. MR uses genetic variants as instrumental variables for the exposure, which should not be influenced by the outcome or confounding variables, overcoming key limitations of traditional epidemiological analyses. For MR to generate a valid estimate of effect, key assumptions must be met. In recent years, there has been a rapid rise in MR methods that aim to test, or are robust to violations of, these assumptions. In this review, we provide an overview of MR for a non-expert audience, including an explanation of these key assumptions and how they are often tested, to aid a better reading and understanding of the MR literature. We highlight some of these new methods and how they can be useful for specific methodological challenges in the musculoskeletal field, including for conditions or traits that share underlying biological pathways, such as bone and joint disease. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

No bidirectional relationship between depression and periodontitis: A genetic correlation and Mendelian randomization study

Background

Observational and in-vivo research suggested a bidirectional relationship between depression and periodontitis. We estimated the genetic correlation and examined directionality of causation.

Methods

The study used summary statistics from published genome wide association studies, with sample sizes ranging from 45,563 to 797,563 individuals of European ancestry. We performed linkage disequilibrium score regression (LDSC) to estimate global correlation and used Heritability Estimation from Summary Statistics (ρ-HESS) to further examine local genetic correlation. Latent Heritable Confounder Mendelian randomization (LHC-MR), Causal Analysis using Summary Effect estimates (CAUSE), and conventional MR approaches assessed bidirectional causation.

Results

LDSC observed only weak genetic correlation (r_g = 0.06, P-Value = 0.619) between depression and periodontitis. Analysis of local genetic correlation using ρ-HESS did not reveal loci of significant local genetic covariance. LHC-MR, CAUSE and conventional MR models provided no support for bidirectional causation between depression and periodontitis, with odds ratios ranging from 1.00 to 1.06 in either direction.

Conclusions

Results do not support shared heritability or a causal connection between depression and periodontitis.

Adjusting for collider bias in genetic association studies using instrumental variable methods

Genome‐wide association studies have provided many genetic markers that can be used as instrumental variables to adjust for confounding in epidemiological studies. Recently, the principle has been applied to other forms of bias in observational studies, especially collider bias that arises when conditioning or stratifying on a variable that is associated with the outcome of interest. An important case is in studies of disease progression and survival. Here, we clarify the links between the genetic instrumental variable methods proposed for this problem and the established methods of Mendelian randomisation developed to account for confounding. We highlight the critical importance of weak instrument bias in this context and describe a corrected weighted least‐squares procedure as a simple approach to reduce this bias. We illustrate the range of available methods on two data examples. The first, waist–hip ratio adjusted for body‐mass index, entails statistical adjustment for a quantitative trait. The second, smoking cessation, is a stratified analysis conditional on having initiated smoking. In both cases, we find little effect of collider bias on the primary association results, but this may propagate into more substantial effects on further analyses such as polygenic risk scoring and Mendelian randomisation.

Robust inference of bi-directional causal relationships in presence of correlated pleiotropy with GWAS summary data

To infer a causal relationship between two traits, several correlation-based causal direction (CD) methods have been proposed with the use of SNPs as instrumental variables (IVs) based on GWAS summary data for the two traits; however, none of the existing CD methods can deal with SNPs with correlated pleiotropy. Alternatively, reciprocal Mendelian randomization (MR) can be applied, which however may perform poorly in the presence of (unknown) invalid IVs, especially for bi-directional causal relationships. In this paper, first, we propose a CD method that performs better than existing CD methods regardless of the presence of correlated pleiotropy. Second, along with a simple but yet effective IV screening rule, we propose applying a closely related and state-of-the-art MR method in reciprocal MR, showing its almost identical performance to that of the new CD method when their model assumptions hold; however, if the modeling assumptions are violated, the new CD method is expected to better control type I errors. Notably bi-directional causal relationships impose some unique challenges beyond those for uni-directional ones, and thus requiring special treatments. For example, we point out for the first time several scenarios where a bi-directional relationship, but not a uni-directional one, can unexpectedly cause the violation of some weak modeling assumptions commonly required by many robust MR methods. We also offer some numerical support and a modeling justification for the application of our new methods (and more generally MR) to binary traits. Finally we applied the proposed methods to 12 risk factors and 4 common diseases, confirming mostly well-known uni-directional causal relationships, while identifying some novel and plausible bi-directional ones such as between body mass index and type 2 diabetes (T2D), and between diastolic blood pressure and stroke.

Inferring causal relationships between two traits based on observational data is one of the most important as well as challenging problems in scientific research. The potential impact of such an application on and beyond genetics/genomics is significant, such as in prioritizing molecular, clinical and behavioral targets for therapeutic and behavioral interventions. However, this problem, especially on a possibly bi-directional causal relationship between a pair of traits (in which the causal direction may be from either trait to the other, and both causal directions may be present at the same time), had been largely neglected until very recently. The increasing availability of large-scale GWAS summary data of various traits has popularized the development and application of Mendelian randomization (MR) methods for such a purpose. We point out some severe limitations with the current methods, mainly due to some new and unique challenges facing inference of bi-directional relationships as compared to that of uni-directional relationships that has been more commonly and exclusively considered in MR. By combining two basic ideas of bidirectional (or reciprocal) MR and Steiger’s correlation-based screening methods, we develop two new approaches based on constrained maximum likelihood (cML) and GWAS summary data to infer causal effects (as in typical MR) and SNP-trait correlations (as in Steiger’s method), called MR-cML and CD-cML respectively, demonstrating their similar effectiveness and more importantly their advantages over existing methods through extensive simulations, real data examples and statistical theory. In particular, our proposed two methods are robust to violations of all three valid IV assumptions, including presence of correlated pleiotropy.