Noise-augmented directional clustering of genetic association data identifies distinct mechanisms underlying obesity

Exploring the shared genetic landscape of diabetes and cardiovascular disease: findings and future implications

Diabetes is a rapidly growing global health concern projected to affect one in eight adults by 2045, which translates to roughly 783 million people. The profound metabolic alterations often present in dysglycaemia significantly increase the risk of cardiovascular complications. While genetic susceptibility plays a crucial role in diabetes and its vascular complications, identifying genes and molecular mechanisms that influence both diseases simultaneously has proven challenging. A key reason for this challenge is the pathophysiological heterogeneity underlying these diseases, with multiple processes contributing to different forms of diabetes and specific cardiovascular complications. This molecular heterogeneity has limited the effectiveness of large-scale genome-wide association studies (GWAS) in identifying shared underlying mechanisms. Additionally, our limited knowledge of the causal genes, cell types and disease-relevant states through which GWAS signals operate has hindered the discovery of common molecular pathways. This review highlights recent advances in genetic epidemiology, including studies of causal associations that have uncovered genetic and molecular factors influencing both dysglycaemia and cardiovascular complications. We explore how disease subtyping approaches can be critical in pinpointing the unique molecular signatures underlying both diabetes and cardiovascular complications. Finally, we address critical research gaps and future opportunities to advance our understanding of both diseases and translate these discoveries into tangible benefits for patient care and population health.

Distinct pathway-based effects of blood pressure and body mass index on cardiovascular traits: comparison of novel Mendelian randomization approaches

BACKGROUND

Mendelian randomization (MR) leverages trait associated genetic variants as instrumental variables (IVs) to determine causal relationships in epidemiology. However, genetic IVs for complex traits are typically highly heterogeneous and, at a molecular level, exert effects on different biological processes. Exploration of the biological underpinnings of such heterogeneity can enhance our understanding of disease mechanisms and inform therapeutic strategies. Here, we introduce a new approach to instrument partitioning based on enrichment of Mendelian disease categories (pathway-partitioned) and compare it to an existing method based on genetic colocalization in contrasting tissues (tissue-partitioned).

METHODS

We employed individual- and summary-level MR methodologies using SNPs grouped by pathway informed by proximity to Mendelian disease genes affecting the renal system or vasculature (for blood pressure (BP)), or mental health and metabolic disorders (for body mass index (BMI)). We compared the causal effects of pathway-partitioned SNPs on cardiometabolic outcomes with those derived using tissue-partitioned SNPs informed by colocalization with gene expression in kidney, artery (BP), or adipose and brain tissues (BMI). Additionally, we assessed the likelihood that estimates observed for partitioned exposures could emerge by chance using random SNP sampling.

RESULTS

Our pathway-partitioned findings suggest the causal relationship between systolic BP and heart disease is predominantly driven by vessel over renal pathways. The stronger effect attributed to kidney over artery tissue in our tissue-partitioned MR hints at a multifaceted interplay between pathways in the disease aetiology. We consistently identified a dominant role for vessel (pathway) and artery (tissue) driving the negative directional effect of diastolic BP on left ventricular stroke volume and positive directional effect of systolic BP on type 2 diabetes. We also found when dissecting the BMI pathway contribution to atrial fibrillation that metabolic-pathway and brain-tissue IVs predominantly drove the causal effects relative to mental health and adipose in pathway- and tissue-partitioned MR analyses, respectively.

CONCLUSIONS

This study presents a novel approach to dissecting heterogeneity in MR by integrating clinical phenotypes associated with Mendelian disease. Our findings emphasize the importance of understanding pathway-/tissue-specific contributions to complex exposures when interpreting causal relationships in MR. Importantly, we advocate caution and robust validation when interpreting pathway-partitioned effect size differences.

Genetic subtyping of obesity reveals biological insights into the uncoupling of adiposity from its cardiometabolic comorbidities

Obesity is a highly heterogeneous disease that cannot be captured by one single adiposity trait. Here, we performed a multi-trait analysis to study obesity in the context of its common cardiometabolic comorbidities, acknowledging that not all individuals with obesity suffer from cardiometabolic comorbidities and that not all those with normal weight clinically present without them. We leveraged individual-level genotype-phenotype data of 452,768 individuals from the UK Biobank and designed uncoupling phenotypes that are continuous and range from high adiposity with a healthy cardiometabolic profile to low adiposity with an unhealthy cardiometabolic profile. Genome-wide association analyses of these uncoupling phenotypes identified 266 independent variants across 205 genomic loci where the adiposity-increasing allele is also associated with a lower cardiometabolic risk. Consistent with the individual variant effects, a genetic score (GRSuncoupling) that aggregates the uncoupling effects of the 266 variants was associated with lower risk of cardiometabolic disorders, including dyslipidemias (OR=0.92, P=1.4×10-89), type 2 diabetes (OR=0.94, P=6×10-21), and ischemic heart disease (OR=0.96, P=7×10-11), despite a higher risk of obesity (OR=1.16, P=4×10-108), which is in sharp contrast to the association profile observed for the adiposity score (GRSBFP). Nevertheless, a higher GRSuncoupling score was also associated with a higher risk of other, mostly weight-bearing disorders, to the same extent as the GRSBFP. The 266 variants clustered into eight subsets, each representing a genetic subtype of obesity with a distinct cardiometabolic risk profile, characterized by specific underlying pathways. Association of GRSuncoupling and GRSBFP with levels of 2,920 proteins in plasma found 208 proteins to be associated with both scores. The majority (85%) of these overlapping GRS-protein associations were directionally consistent, suggesting adiposity-driven effects. In contrast, levels of 32 (15%) proteins (e.g. IGFBP1, IGFBP2, LDLR, SHBG, MSTN) had opposite directional effects between GRSBFP and GRSuncoupling, suggesting that cardiometabolic health, and not adiposity, associated with their levels. Follow-up analyses provide further support for adipose tissue expandability, insulin secretion and beta-cell function, beiging of white adipose tissue, inflammation and fibrosis. They also highlight mechanisms not previously implicated in uncoupling, such as hepatic lipid accumulation, hepatic control of glucose homeostasis, and skeletal muscle growth and function. Taken together, our findings contribute new insights into the mechanisms that uncouple adiposity from its cardiometabolic comorbidities and illuminate some of the heterogeneity of obesity, which is critical for advancing precision medicine.

Using clustering of genetic variants in Mendelian randomization to interrogate the causal pathways underlying multimorbidity from a common risk factor

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.

Roles of naïve CD4+ T cells and their differentiated subtypes in lung adenocarcinoma and underlying potential regulatory pathways

BACKGROUND

Naïve CD4+ T cells and their differentiated counterparts play a significant regulatory role in the tumor immune microenvironment, yet their effects on lung adenocarcinoma (LUAD) are not fully understood.

METHODS

We utilized Mendelian randomization to assess the causal association between naïve CD4+ T cells and LUAD. Employing a modified single-sample Gene Set Enrichment Analysis (ssGSEA) algorithm with The Cancer Genome Atlas (TCGA) database, we determined the infiltration levels of naïve CD4+ T cells and their differentiation subtypes and investigated their correlation with clinical characteristics. Potential regulatory pathways of T helper cells were identified through Mantel tests and Kyoto Encyclopedia of Genes and Genomes (KEGG) database enrichment analysis.

RESULTS

Mendelian randomization analysis revealed an inhibitory effect of naïve CD4+ T cells on LUAD (false discovery rate < 0.05), which was corroborated by observational experiments using TCGA database. Specifically, T helper cell type 2 demonstrated a promotive effect on LUAD in terms of overall, disease-free, and progression-free survival (p < 0.05). Moreover, regulatory T cells exhibited a protective effect on LUAD in terms of disease-specific survival (p < 0.01). Concurrently, we explored the overall impact of naïve CD4+ T cell differentiation subtypes on LUAD, revealing upregulation in pathways such as neutrophil degranulation, MAPK family signaling pathways, and platelet activation, signaling, and aggregation.

CONCLUSION

Naïve CD4+ T cells and their differentiated counterparts play essential regulatory roles in the tumor immune microenvironment, demonstrating bidirectionality in their effects.Thus, elucidating the mechanisms and developing novel cell differentiation-inducing agents will benefit anti-cancer therapy.

Systemic proteome adaptions to 7-day complete caloric restriction in humans

Surviving long periods without food has shaped human evolution. In ancient and modern societies, prolonged fasting was/is practiced by billions of people globally for religious purposes, used to treat diseases such as epilepsy, and recently gained popularity as weight loss intervention, but we still have a very limited understanding of the systemic adaptions in humans to extreme caloric restriction of different durations. Here we show that a 7-day water-only fast leads to an average weight loss of 5.7 kg (±0.8 kg) among 12 volunteers (5 women, 7 men). We demonstrate nine distinct proteomic response profiles, with systemic changes evident only after 3 days of complete calorie restriction based on in-depth characterization of the temporal trajectories of ~3,000 plasma proteins measured before, daily during, and after fasting. The multi-organ response to complete caloric restriction shows distinct effects of fasting duration and weight loss and is remarkably conserved across volunteers with >1,000 significantly responding proteins. The fasting signature is strongly enriched for extracellular matrix proteins from various body sites, demonstrating profound non-metabolic adaptions, including extreme changes in the brain-specific extracellular matrix protein tenascin-R. Using proteogenomic approaches, we estimate the health consequences for 212 proteins that change during fasting across ~500 outcomes and identified putative beneficial (SWAP70 and rheumatoid arthritis or HYOU1 and heart disease), as well as adverse effects. Our results advance our understanding of prolonged fasting in humans beyond a merely energy-centric adaptions towards a systemic response that can inform targeted therapeutic modulation.

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.

A Bayesian approach to Mendelian randomization using summary statistics in the univariable and multivariable settings with correlated pleiotropy

Mendelian randomization uses genetic variants as instrumental variables to make causal inferences on the effect of an exposure on an outcome. Due to the recent abundance of high-powered genome-wide association studies, many putative causal exposures of interest have large numbers of independent genetic variants with which they associate, each representing a potential instrument for use in a Mendelian randomization analysis. Such polygenic analyses increase the power of the study design to detect causal effects; however, they also increase the potential for bias due to instrument invalidity. Recent attention has been given to dealing with bias caused by correlated pleiotropy, which results from violation of the "instrument strength independent of direct effect" assumption. Although methods have been proposed that can account for this bias, a number of restrictive conditions remain in many commonly used techniques. In this paper, we propose a Bayesian framework for Mendelian randomization that provides valid causal inference under very general settings. We propose the methods MR-Horse and MVMR-Horse, which can be performed without access to individual-level data, using only summary statistics of the type commonly published by genome-wide association studies, and can account for both correlated and uncorrelated pleiotropy. In simulation studies, we show that the approach retains type I error rates below nominal levels even in high-pleiotropy scenarios. We demonstrate the proposed approaches in applied examples in both univariable and multivariable settings, some with very weak instruments.

Mendelian Randomization as a Tool for Cardiovascular Research: A Review

Importance

Mendelian randomization (MR) is a statistical approach that has become increasingly popular in the field of cardiovascular disease research. It offers a way to infer potentially causal relationships between risk factors and outcomes using observational data, which is particularly important in cases where randomized clinical trials are not feasible or ethical. With the growing availability of large genetic data sets, MR has become a powerful and accessible tool for studying the risk factors for cardiovascular disease.

Observations

MR uses genetic variation associated with modifiable exposures or risk factors to mitigate biases that affect traditional observational study designs. The approach uses genetic variants that are randomly assigned at conception as proxies for exposure to a risk factor, mimicking a randomized clinical trial. By comparing the outcomes of individuals with different genetic variants, researchers may draw causal inferences about the effects of specific risk factors on cardiovascular disease, provided assumptions are met that address (1) the association between each genetic variant and risk factor and (2) the association of the genetic variants with confounders and (3) that the association between each genetic variant and the outcome only occurs through the risk factor. Like other observational designs, MR has limitations, which include weak instruments that are not strongly associated with the exposure of interest, linkage disequilibrium where genetic instruments influence the outcome via correlated rather than direct effects, overestimated genetic associations, and selection and survival biases. In addition, many genetic databases and MR studies primarily include populations genetically similar to European reference populations; improved diversity of participants in these databases and studies is critically needed.

Conclusions and Relevance

This review provides an overview of MR methodology, including assumptions, strengths, and limitations. Several important applications of MR in cardiovascular disease research are highlighted, including the identification of drug targets, evaluation of potential cardiovascular risk factors, as well as emerging methodology. Overall, while MR alone can never prove a causal relationship beyond reasonable doubt, MR offers a rigorous approach for investigating possible causal relationships in observational data and has the potential to transform our understanding of the etiology and treatment of cardiovascular disease.

Using genetic association data to guide drug discovery and development: Review of methods and applications

Evidence on the validity of drug targets from randomized trials is reliable but typically expensive and slow to obtain. In contrast, evidence from conventional observational epidemiological studies is less reliable because of the potential for bias from confounding and reverse causation. Mendelian randomization is a quasi-experimental approach analogous to a randomized trial that exploits naturally occurring randomization in the transmission of genetic variants. In Mendelian randomization, genetic variants that can be regarded as proxies for an intervention on the proposed drug target are leveraged as instrumental variables to investigate potential effects on biomarkers and disease outcomes in large-scale observational datasets. This approach can be implemented rapidly for a range of drug targets to provide evidence on their effects and thus inform on their priority for further investigation. In this review, we present statistical methods and their applications to showcase the diverse opportunities for applying Mendelian randomization in guiding clinical development efforts, thus enabling interventions to target the right mechanism in the right population group at the right time. These methods can inform investigators on the mechanisms underlying drug effects, their related biomarkers, implications for the timing of interventions, and the population subgroups that stand to gain the most benefit. Most methods can be implemented with publicly available data on summarized genetic associations with traits and diseases, meaning that the only major limitations to their usage are the availability of appropriately powered studies for the exposure and outcome and the existence of a suitable genetic proxy for the proposed intervention.

Therapeutic potential of IL6R blockade for the treatment of sepsis and sepsis-related death: A Mendelian randomisation study

BACKGROUND

Sepsis is characterised by dysregulated, life-threatening immune responses, which are thought to be driven by cytokines such as interleukin 6 (IL-6). Genetic variants in IL6R known to down-regulate IL-6 signalling are associated with improved Coronavirus Disease 2019 (COVID-19) outcomes, a finding later confirmed in randomised trials of IL-6 receptor antagonists (IL6RAs). We hypothesised that blockade of IL6R could also improve outcomes in sepsis.

METHODS AND FINDINGS

We performed a Mendelian randomisation (MR) analysis using single nucleotide polymorphisms (SNPs) in and near IL6R to evaluate the likely causal effects of IL6R blockade on sepsis (primary outcome), sepsis severity, other infections, and COVID-19 (secondary outcomes). We weighted SNPs by their effect on CRP and combined results across them in inverse variance weighted meta-analysis, proxying the effect of IL6RA. Our outcomes were measured in UK Biobank, FinnGen, the COVID-19 Host Genetics Initiative (HGI), and the GenOSept and GainS consortium. We performed several sensitivity analyses to test assumptions of our methods, including utilising variants around CRP and gp130 in a similar analysis. In the UK Biobank cohort (N = 486,484, including 11,643 with sepsis), IL6R blockade was associated with a decreased risk of our primary outcome, sepsis (odds ratio (OR) = 0.80; 95% confidence interval (CI) 0.66 to 0.96, per unit of natural log-transformed CRP decrease). The size of this effect increased with severity, with larger effects on 28-day sepsis mortality (OR = 0.74; 95% CI 0.47 to 1.15); critical care admission with sepsis (OR = 0.48, 95% CI 0.30 to 0.78) and critical care death with sepsis (OR = 0.37, 95% CI 0.14 to 0.98). Similar associations were seen with severe respiratory infection: OR for pneumonia in critical care 0.69 (95% CI 0.49 to 0.97) and for sepsis survival in critical care (OR = 0.22; 95% CI 0.04 to 1.31) in the GainS and GenOSept consortium, although this result had a large degree of imprecision. We also confirm the previously reported protective effect of IL6R blockade on severe COVID-19 (OR = 0.69, 95% CI 0.57 to 0.84) in the COVID-19 HGI, which was of similar magnitude to that seen in sepsis. Sensitivity analyses did not alter our primary results. These results are subject to the limitations and assumptions of MR, which in this case reflects interpretation of these SNP effects as causally acting through blockade of IL6R, and reflect lifetime exposure to IL6R blockade, rather than the effect of therapeutic IL6R blockade.

CONCLUSIONS

IL6R blockade is causally associated with reduced incidence of sepsis. Similar but imprecisely estimated results supported a causal effect also on sepsis related mortality and critical care admission with sepsis. These effects are comparable in size to the effect seen in severe COVID-19, where IL-6 receptor antagonists were shown to improve survival. These data suggest that a randomised trial of IL-6 receptor antagonists in sepsis should be considered.