Pleiotropy-Based Decomposition of Genetic Risk Scores: Association and Interaction Analysis for Type 2 Diabetes and CAD

Defining type 2 diabetes polygenic risk scores through colocalization and network-based clustering of metabolic trait genetic associations

BACKGROUND

Type 2 diabetes (T2D) is a heterogeneous and polygenic disease. Previous studies have leveraged the highly polygenic and pleiotropic nature of T2D variants to partition the heterogeneity of T2D, in order to stratify patient risk and gain mechanistic insight. We expanded on these approaches by performing colocalization across GWAS traits while assessing the causality and directionality of genetic associations.

METHODS

We applied colocalization between T2D and 20 related metabolic traits, across 243 loci, to obtain inferences of shared casual variants. Network-based unsupervised hierarchical clustering was performed on variant-trait associations. Partitioned polygenic risk scores (PRSs) were generated for each cluster using T2D summary statistics and validated in 21,742 individuals with T2D from 3 cohorts. Inferences of directionality and causality were obtained by applying Mendelian randomization Steiger's Z-test and further validated in a pediatric cohort without diabetes (aged 9-12 years old, n = 3866).

RESULTS

We identified 146 T2D loci that colocalized with at least one metabolic trait locus. T2D variants within these loci were grouped into 5 clusters. The clusters corresponded to the following pathways: obesity, lipodystrophic insulin resistance, liver and lipid metabolism, hepatic glucose metabolism, and beta-cell dysfunction. We observed heterogeneity in associations between PRSs and metabolic measures across clusters. For instance, the lipodystrophic insulin resistance (Beta - 0.08 SD, 95% CI [- 0.10-0.07], p = 6.50 × 10-32) and beta-cell dysfunction (Beta - 0.10 SD, 95% CI [- 0.12, - 0.08], p = 1.46 × 10-47) PRSs were associated to lower BMI. Mendelian randomization Steiger analysis indicated that increased T2D risk in these pathways was causally associated to lower BMI. However, the obesity PRS was conversely associated with increased BMI (Beta 0.08 SD, 95% CI 0.06-0.10, p = 8.0 × 10-33). Analyses within a pediatric cohort supported this finding. Additionally, the lipodystrophic insulin resistance PRS was associated with a higher odds of chronic kidney disease (OR 1.29, 95% CI 1.02-1.62, p = 0.03).

CONCLUSIONS

We successfully partitioned T2D genetic variants into phenotypic pathways using a colocalization first approach. Partitioned PRSs were associated to unique metabolic and clinical outcomes indicating successful partitioning of disease heterogeneity. Our work expands on previous approaches by providing stronger inferences of shared causal variants, causality, and directionality of GWAS variant-trait associations.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis

BACKGROUND

Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D.

METHODS

We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test.

RESULTS

Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis.

CONCLUSIONS

These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

A scalable approach to characterize pleiotropy across thousands of human diseases and complex traits using GWAS summary statistics

Genome-wide association studies (GWASs) across thousands of traits have revealed the pervasive pleiotropy of trait-associated genetic variants. While methods have been proposed to characterize pleiotropic components across groups of phenotypes, scaling these approaches to ultra-large-scale biobanks has been challenging. Here, we propose FactorGo, a scalable variational factor analysis model to identify and characterize pleiotropic components using biobank GWAS summary data. In extensive simulations, we observe that FactorGo outperforms the state-of-the-art (model-free) approach tSVD in capturing latent pleiotropic factors across phenotypes while maintaining a similar computational cost. We apply FactorGo to estimate 100 latent pleiotropic factors from GWAS summary data of 2,483 phenotypes measured in European-ancestry Pan-UK BioBank individuals (N = 420,531). Next, we find that factors from FactorGo are more enriched with relevant tissue-specific annotations than those identified by tSVD (p = 2.58E-10) and validate our approach by recapitulating brain-specific enrichment for BMI and the height-related connection between reproductive system and muscular-skeletal growth. Finally, our analyses suggest shared etiologies between rheumatoid arthritis and periodontal condition in addition to alkaline phosphatase as a candidate prognostic biomarker for prostate cancer. Overall, FactorGo improves our biological understanding of shared etiologies across thousands of GWASs.

The triglyceride-glucose index is associated with atherosclerosis in patients with symptomatic coronary artery disease, regardless of diabetes mellitus and hyperlipidaemia

BACKGROUND

Diabetes and hyperlipidaemia are both risk factors for coronary artery disease, and both are associated with a high triglyceride-glucose index (TyG index). The TyG index has been presented as a marker of insulin resistance (IR). Its utility in predicting and detecting cardiovascular disease has been reported. However, few studies have found it to be a helpful marker of atherosclerosis in patients with symptomatic coronary artery disease (CAD). The purpose of this study was to demonstrate that the TyG index can serve as a valuable marker for predicting coronary and carotid atherosclerosis in symptomatic CAD patients, regardless of diabetes mellitus and hyperlipidaemia.

METHODS

This study included 1516 patients with symptomatic CAD who underwent both coronary artery angiography and carotid Doppler ultrasound in the Department of Cardiology at Tianjin Union Medical Center from January 2016 to December 2022. The TyG index was determined using the Ln formula. The population was further grouped and analysed according to the presence or absence of diabetes and hyperlipidaemia. The Gensini score and carotid intima-media thickness were calculated or measured, and the patients were divided into four groups according to TyG index quartile to examine the relationship between the TyG index and coronary or carotid artery lesions in symptomatic CAD patients.

RESULTS

In symptomatic CAD patients, the TyG index showed a significant positive correlation with both coronary lesions and carotid plaques. After adjusting for sex, age, smoking, BMI, hypertension, diabetes, and the use of antilipemic and antidiabetic agents, the risk of developing coronary lesions and carotid plaques increased across the baseline TyG index. Compared with the lowest quartile of the TyG index, the highest quartile (quartile 4) was associated with a greater incidence of coronary heart disease [OR = 2.55 (95% CI 1.61, 4.03)] and carotid atherosclerotic plaque [OR = 2.31 (95% CI 1.27, 4.20)] (P < 0.05). Furthermore, when compared to the fasting blood glucose (FBG) or triglyceride (TG) level, the TyG index had a greater area under the ROC curve for predicting coronary lesions and carotid plaques. The subgroup analysis demonstrated the TyG index to be an equally effective predictor of coronary and carotid artery disease, regardless of diabetes and hyperlipidaemia.

CONCLUSION

The TyG index is a useful marker for predicting coronary and carotid atherosclerosis in patients with symptomatic CAD, regardless of diabetes mellitus and hyperlipidaemia. The TyG index is of higher value for the identification of both coronary and carotid atherosclerotic plaques than the FBG or TG level alone.

Lipid-Associated GWAS Loci Predict Antiatherogenic Effects of Rosuvastatin in Patients with Coronary Artery Disease

We have shown that lipid-associated loci discovered by genome-wide association studies (GWAS) have pleiotropic effects on lipid metabolism, carotid intima-media thickness (CIMT), and CAD risk. Here, we investigated the impact of lipid-associated GWAS loci on the efficacy of rosuvastatin therapy in terms of changes in plasma lipid levels and CIMT. The study comprised 116 CAD patients with hypercholesterolemia. CIMT, total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG) were measured at baseline and after 6 and 12 months of follow-up, respectively. Genotyping of fifteen lipid-associated GWAS loci was performed by the MassArray-4 System. Linear regression analysis adjusted for sex, age, body mass index, and rosuvastatin dose was used to estimate the phenotypic effects of polymorphisms, and p-values were calculated through adaptive permutation tests by the PLINK software, v1.9. Over one-year rosuvastatin therapy, a decrease in CIMT was linked to rs1689800, rs4846914, rs12328675, rs55730499, rs9987289, rs11220463, rs16942887, and rs881844 polymorphisms (Pperm < 0.05). TC change was associated with rs55730499, rs11220463, and rs6065906; LDL-C change was linked to the rs55730499, rs1689800, and rs16942887 polymorphisms; and TG change was linked to polymorphisms rs838880 and rs1883025 (Pperm < 0.05). In conclusion, polymorphisms rs1689800, rs55730499, rs11220463, and rs16942887 were found to be predictive markers for multiple antiatherogenic effects of rosuvastatin in CAD patients.

Polygenic Risk Score in African populations: progress and challenges

Polygenic Risk Score (PRS) analysis is a method that predicts the genetic risk of an individual towards targeted traits. Even when there are no significant markers, it gives evidence of a genetic effect beyond the results of Genome-Wide Association Studies (GWAS). Moreover, it selects  single nucleotide polymorphisms (SNPs) that  contribute to the disease with low effect size  making it more precise at individual level risk prediction. PRS  analysis addresses the shortfall of GWAS by taking into account the SNPs/alleles with  low effect size but play an indispensable role to the observed phenotypic/trait variance.  PRS analysis has  applications that investigate the genetic basis of several traits, which includes rare diseases. However, the accuracy of PRS analysis depends on the genomic data of the underlying population. For instance, several studies  show   that obtaining higher prediction power of PRS analysis is challenging for non-Europeans. In this manuscript, we review the conventional PRS methods and their application to sub-Saharan African communities. We conclude that  lack of sufficient GWAS data and tools is  the limiting factor of applying PRS analysis to sub-Saharan populations.   We recommend developing Africa-specific PRS methods and tools for estimating and analyzing  African population data   for clinical  evaluation of PRSs of interest and predicting  rare diseases.

Polygenic Risk Score in African populations: progress and challenges

Polygenic Risk Score (PRS) analysis is a method that predicts the genetic risk of an individual towards targeted traits. Even when there are no significant markers, it gives evidence of a genetic effect beyond the results of Genome-Wide Association Studies (GWAS). Moreover, it selects  single nucleotide polymorphisms (SNPs) that  contribute to the disease with low effect size  making it more precise at individual level risk prediction. PRS  analysis addresses the shortfall of GWAS by taking into account the SNPs/alleles with  low effect size but play an indispensable role to the observed phenotypic/trait variance.  PRS analysis has  applications that investigate the genetic basis of several traits, which includes rare diseases. However, the accuracy of PRS analysis depends on the genomic data of the underlying population. For instance, several studies  show   that obtaining higher prediction power of PRS analysis is challenging for non-Europeans. In this manuscript, we review the conventional PRS methods and their application to sub-Saharan African communities. We conclude that  lack of sufficient GWAS data and tools is  the limiting factor of applying PRS analysis to sub-Saharan populations.   We recommend developing Africa-specific PRS methods and tools for estimating and analyzing  African population data   for clinical  evaluation of PRSs of interest and predicting  rare diseases.

A scalable variational approach to characterize pleiotropic components across thousands of human diseases and complex traits using GWAS summary statistics

Genome-wide association studies (GWAS) across thousands of traits have revealed the pervasive pleiotropy of trait-associated genetic variants. While methods have been proposed to characterize pleiotropic components across groups of phenotypes, scaling these approaches to ultra large-scale biobanks has been challenging. Here, we propose FactorGo, a scalable variational factor analysis model to identify and characterize pleiotropic components using biobank GWAS summary data. In extensive simulations, we observe that FactorGo outperforms the state-of-the-art (model-free) approach tSVD in capturing latent pleiotropic factors across phenotypes, while maintaining a similar computational cost. We apply FactorGo to estimate 100 latent pleiotropic factors from GWAS summary data of 2,483 phenotypes measured in European-ancestry Pan-UK BioBank individuals (N=420,531). Next, we find that factors from FactorGo are more enriched with relevant tissue-specific annotations than those identified by tSVD (P=2.58E-10), and validate our approach by recapitulating brain-specific enrichment for BMI and the height-related connection between reproductive system and muscular-skeletal growth. Finally, our analyses suggest novel shared etiologies between rheumatoid arthritis and periodontal condition, in addition to alkaline phosphatase as a candidate prognostic biomarker for prostate cancer. Overall, FactorGo improves our biological understanding of shared etiologies across thousands of GWAS.

Do gene-environment interactions have implications for the precision prevention of type 2 diabetes?

The past decades have seen a rapid global rise in the incidence of type 2 diabetes. This surge has been driven by diabetogenic environmental changes that may act together with a genetic predisposition to type 2 diabetes. It is possible that there is a synergistic gene-environment interaction, where the effects of the diabetogenic environment depend on the genetic predisposition to type 2 diabetes. Randomised trials have shown that it is possible to delay, or even prevent the development of type 2 diabetes in individuals at elevated risk through behavioural modification, focusing on weight loss, physical activity and diet. There is wide heterogeneity between individuals regarding the effectiveness of these interventions, which could, in part, be due to genetic differences. However, the studies of gene-environment interactions performed thus far suggest that behavioural modifications appear equally effective in reducing the incidence of type 2 diabetes from the stage of impaired glucose tolerance, regardless of the known underlying genetic predisposition. Recent studies suggest that there may be several subtypes of type 2 diabetes, which give new opportunities for gaining insight into gene-environment interactions. At present, the role of gene-environment interactions in the development of type 2 diabetes remains unclear. With many puzzle pieces missing in the general picture of type 2 diabetes development, the available evidence of gene-environment interactions is not ready for translation to individualised type 2 diabetes prevention based on genetic profiling.

Chronic neuropsychiatric sequelae of SARS-CoV-2: Protocol and methods from the Alzheimer's Association Global Consortium

Introduction

Coronavirus disease 2019 (COVID-19) has caused >3.5 million deaths worldwide and affected >160 million people. At least twice as many have been infected but remained asymptomatic or minimally symptomatic. COVID-19 includes central nervous system manifestations mediated by inflammation and cerebrovascular, anoxic, and/or viral neurotoxicity mechanisms. More than one third of patients with COVID-19 develop neurologic problems during the acute phase of the illness, including loss of sense of smell or taste, seizures, and stroke. Damage or functional changes to the brain may result in chronic sequelae. The risk of incident cognitive and neuropsychiatric complications appears independent from the severity of the original pulmonary illness. It behooves the scientific and medical community to attempt to understand the molecular and/or systemic factors linking COVID-19 to neurologic illness, both short and long term.

Methods

This article describes what is known so far in terms of links among COVID-19, the brain, neurological symptoms, and Alzheimer's disease (AD) and related dementias. We focus on risk factors and possible molecular, inflammatory, and viral mechanisms underlying neurological injury. We also provide a comprehensive description of the Alzheimer's Association Consortium on Chronic Neuropsychiatric Sequelae of SARS-CoV-2 infection (CNS SC2) harmonized methodology to address these questions using a worldwide network of researchers and institutions.

Results

Successful harmonization of designs and methods was achieved through a consensus process initially fragmented by specific interest groups (epidemiology, clinical assessments, cognitive evaluation, biomarkers, and neuroimaging). Conclusions from subcommittees were presented to the whole group and discussed extensively. Presently data collection is ongoing at 19 sites in 12 countries representing Asia, Africa, the Americas, and Europe.

Discussion

The Alzheimer's Association Global Consortium harmonized methodology is proposed as a model to study long-term neurocognitive sequelae of SARS-CoV-2 infection.

Key Points

The following review describes what is known so far in terms of molecular and epidemiological links among COVID-19, the brain, neurological symptoms, and AD and related dementias (ADRD)The primary objective of this large-scale collaboration is to clarify the pathogenesis of ADRD and to advance our understanding of the impact of a neurotropic virus on the long-term risk of cognitive decline and other CNS sequelae. No available evidence supports the notion that cognitive impairment after SARS-CoV-2 infection is a form of dementia (ADRD or otherwise). The longitudinal methodologies espoused by the consortium are intended to provide data to answer this question as clearly as possible controlling for possible confounders. Our specific hypothesis is that SARS-CoV-2 triggers ADRD-like pathology following the extended olfactory cortical network (EOCN) in older individuals with specific genetic susceptibility.The proposed harmonization strategies and flexible study designs offer the possibility to include large samples of under-represented racial and ethnic groups, creating a rich set of harmonized cohorts for future studies of the pathophysiology, determinants, long-term consequences, and trends in cognitive aging, ADRD, and vascular disease.We provide a framework for current and future studies to be carried out within the Consortium. and offers a "green paper" to the research community with a very broad, global base of support, on tools suitable for low- and middle-income countries aimed to compare and combine future longitudinal data on the topic.The Consortium proposes a combination of design and statistical methods as a means of approaching causal inference of the COVID-19 neuropsychiatric sequelae. We expect that deep phenotyping of neuropsychiatric sequelae may provide a series of candidate syndromes with phenomenological and biological characterization that can be further explored. By generating high-quality harmonized data across sites we aim to capture both descriptive and, where possible, causal associations.

Shared components of heritability across genetically correlated traits

Most disease-associated genetic variants are pleiotropic, affecting multiple genetically correlated traits. Their pleiotropic associations can be mechanistically informative: if many variants have similar patterns of association, they may act via similar pleiotropic mechanisms, forming a shared component of heritability. We developed pleiotropic decomposition regression (PDR) to identify shared components and their underlying genetic variants. We validated PDR on simulated data and identified limitations of existing methods in recovering the true components. We applied PDR to three clusters of five to six traits genetically correlated with coronary artery disease (CAD), asthma, and type II diabetes (T2D), producing biologically interpretable components. For CAD, PDR identified components related to BMI, hypertension, and cholesterol, and it clarified the relationship among these highly correlated risk factors. We assigned variants to components, calculated their posterior-mean effect sizes, and performed out-of-sample validation. Our posterior-mean effect sizes pool statistical power across traits and substantially boost the correlation (r2) between true and estimated effect sizes (compared with the original summary statistics) by 94% and 70% for asthma and T2D out of sample, respectively, and by a predicted 300% for CAD.

Genome-wide pleiotropy analysis identifies novel blood pressure variants and improves its polygenic risk scores

Systolic and diastolic blood pressure (S/DBP) are highly correlated modifiable risk factors for cardiovascular disease (CVD). We report here a bidirectional Mendelian Randomization (MR) and horizontal pleiotropy analysis of S/DBP summary statistics from the UK Biobank (UKB)-International Consortium for Blood Pressure (ICBP) (UKB-ICBP) BP genome-wide association study and construct a composite genetic risk score (GRS) by including pleiotropic variants. The composite GRS captures greater (1.11-3.26 fold) heritability for BP traits and increases (1.09- and 2.01-fold) Nagelkerke's R² for hypertension and CVD. We replicated 118 novel BP horizontal pleiotropic variants including 18 novel BP loci using summary statistics from the Million Veteran Program (MVP) study. An additional 219 novel BP signals and 40 novel loci were identified after a meta-analysis of the UKB-ICBP and MVP summary statistics but without further independent replication. Our study provides further insight into BP regulation and provides a novel way to construct a GRS by including pleiotropic variants for other complex diseases.

Pharmacogenetics of Bronchodilator Response: Future Directions

PURPOSE OF REVIEW

Several genome-wide association studies (GWASs) of bronchodilator response (BDR) to albuterol have been published over the past decade. This review describes current knowledge gaps, including pharmacogenetic studies of albuterol response in minority populations, effect modification of pharmacogenetic associations by age, and relevance of BDR phenotype characterization to pharmacogenetic findings. New approaches, such as leveraging additional "omics" data to focus pharmacogenetic interrogation, as well as developing polygenic risk scores in asthma treatment responses, are also discussed.

RECENT FINDINGS

Recent pharmacogenetic studies of albuterol response in minority populations have identified genetic polymorphisms in loci (DNAH5, NFKB1, PLCB1, ADAMTS3, COX18, and PRKG1), that are associated with BDR. Additional studies are needed to replicate these findings. Modification of the pharmacogenetic associations for SPATS2L and ASB3 polymorphisms by age has also been published. Evidence from metabolomic and epigenomic studies of BDR may point to new pharmacogenetic targets. Lastly, a polygenic risk score for response to albuterol has been developed but requires validation in additional cohorts. In order to expand our knowledge of pharmacogenetics of BDR, additional studies in minority populations are needed. Consideration of effect modification by age and leverage of other "omics" data beyond genomics may also help uncover novel pharmacogenetic loci for use in precision medicine for asthma treatment.

A cross-population atlas of genetic associations for 220 human phenotypes

Current genome-wide association studies do not yet capture sufficient diversity in populations and scope of phenotypes. To expand an atlas of genetic associations in non-European populations, we conducted 220 deep-phenotype genome-wide association studies (diseases, biomarkers and medication usage) in BioBank Japan (n = 179,000), by incorporating past medical history and text-mining of electronic medical records. Meta-analyses with the UK Biobank and FinnGen (n_total = 628,000) identified ~5,000 new loci, which improved the resolution of the genomic map of human traits. This atlas elucidated the landscape of pleiotropy as represented by the major histocompatibility complex locus, where we conducted HLA fine-mapping. Finally, we performed statistical decomposition of matrices of phenome-wide summary statistics, and identified latent genetic components, which pinpointed responsible variants and biological mechanisms underlying current disease classifications across populations. The decomposed components enabled genetically informed subtyping of similar diseases (for example, allergic diseases). Our study suggests a potential avenue for hypothesis-free re-investigation of human diseases through genetics.