Pleiotropy informed adaptive association test of multiple traits using genome-wide association study summary data

Genome-wide association study and polygenic score assessment of insulin resistance

Insulin resistance (IR) and beta cell dysfunction are the major drivers of type 2 diabetes (T2D). Genome-Wide Association Studies (GWAS) on IR have been predominantly conducted in European populations, while Middle Eastern populations remain largely underrepresented. We conducted a GWAS on the indices of IR (HOMA2-IR) and beta cell function (HOMA2-%B) in 6,217 non-diabetic individuals from the Qatar Biobank (QBB; Discovery cohort; n = 2170, Replication cohort; n = 4047) with and without body mass index (BMI) adjustment. We also developed polygenic scores (PGS) for HOMA2-IR and compared their performance with a previously derived PGS for HOMA-IR (PGS003470). We replicated 11 loci that have been previously associated with HOMA-IR and 24 loci that have been associated with HOMA-%B, at nominal statistical significance. We also identified a novel locus associated with beta cell function near VEGFC gene, tagged by rs61552983 (P = 4.38 × 10-8). Moreover, our best performing PGS (Q-PGS4; Adj R2 = 0.233 ± 0.014; P = 1.55 x 10-3) performed better than PGS003470 (Adj R2 = 0.194 ± 0.014; P = 5.45 x 10-2) in predicting HOMA2-IR in our dataset. This is the first GWAS on HOMA2 and the first GWAS conducted in the Middle East focusing on IR and beta cell function. Herein, we report a novel locus in VEGFC that is implicated in beta cell dysfunction. Inclusion of under-represented populations in GWAS has potentials to provide important insights into the genetic architecture of IR and beta cell function.

Identifying pleiotropic genes via the composite test amidst the complexity of polygenic traits

Identifying the causal relationship between genotype and phenotype is essential to expanding our understanding of the gene regulatory network spanning the molecular level to perceptible traits. A pleiotropic gene can act as a central hub in the network, influencing multiple outcomes. Identifying such a gene involves testing under a composite null hypothesis where the gene is associated with, at most, one trait. Traditional methods such as meta-analyses of top-hit $P$-values and sequential testing of multiple traits have been proposed, but these methods fail to consider the background of genome-wide signals. Since Huang's composite test produces uniformly distributed $P$-values for genome-wide variants under the composite null, we propose a gene-level pleiotropy test that entails combining the aforementioned method with the aggregated Cauchy association test. A polygenic trait involves multiple genes with different functions to co-regulate mechanisms. We show that polygenicity should be considered when identifying pleiotropic genes; otherwise, the associations polygenic traits initiate will give rise to false positives. In this study, we constructed gene-trait functional modules using the results of the proposed pleiotropy tests. Our analysis suite was implemented as an R package PGCtest. We demonstrated the proposed method with an application study of the Taiwan Biobank database and identified functional modules comprising specific genes and their co-regulated traits.

Dissecting the shared genetic architecture between Alzheimer's disease and frailty: a cross-trait meta-analyses of genome-wide association studies

Introduction: Frailty is the most common medical condition affecting the aging population, and its prevalence increases in the population aged 65 or more. Frailty is commonly diagnosed using the frailty index (FI) or frailty phenotype (FP) assessments. Observational studies have indicated the association of frailty with Alzheimer's disease (AD). However, the shared genetic and biological mechanism of these comorbidity has not been studied. Methods: To assess the genetic relationship between AD and frailty, we examined it at single nucleotide polymorphism (SNP), gene, and pathway levels. Results: Overall, 16 genome-wide significant loci (15 unique loci) (p meta-analysis < 5 × 10-8) and 22 genes (21 unique genes) were identified between AD and frailty using cross-trait meta-analysis. The 8 shared loci implicated 11 genes: CLRN1-AS1, CRHR1, FERMT2, GRK4, LINC01929, LRFN2, MADD, RP11-368P15.1, RP11-166N6.2, RNA5SP459, and ZNF652 between AD and FI, and 8 shared loci between AD and FFS implicated 11 genes: AFF3, C1QTNF4, CLEC16A, FAM180B, FBXL19, GRK4, LINC01104, MAD1L1, RGS12, ZDHHC5, and ZNF521. The loci 4p16.3 (GRK4) was identified in both meta-analyses. The colocalization analysis supported the results of our meta-analysis in these loci. The gene-based analysis revealed 80 genes between AD and frailty, and 4 genes were initially identified in our meta-analyses: C1QTNF4, CRHR1, MAD1L1, and RGS12. The pathway analysis showed enrichment for lipoprotein particle plasma, amyloid fibril formation, protein kinase regulator, and tau protein binding. Conclusion: Overall, our results provide new insights into the genetics of AD and frailty, suggesting the existence of non-causal shared genetic mechanisms between these conditions.

DNA methylation may partly explain psychotropic drug-induced metabolic side effects: results from a prospective 1-month observational study

BACKGROUND

Metabolic side effects of psychotropic medications are a major drawback to patients' successful treatment. Using an epigenome-wide approach, we aimed to investigate DNA methylation changes occurring secondary to psychotropic treatment and evaluate associations between 1-month metabolic changes and both baseline and 1-month changes in DNA methylation levels. Seventy-nine patients starting a weight gain inducing psychotropic treatment were selected from the PsyMetab study cohort. Epigenome-wide DNA methylation was measured at baseline and after 1 month of treatment, using the Illumina Methylation EPIC BeadChip.

RESULTS

A global methylation increase was noted after the first month of treatment, which was more pronounced (p < 2.2 × 10-16) in patients whose weight remained stable (< 2.5% weight increase). Epigenome-wide significant methylation changes (p < 9 × 10-8) were observed at 52 loci in the whole cohort. When restricting the analysis to patients who underwent important early weight gain (≥ 5% weight increase), one locus (cg12209987) showed a significant increase in methylation levels (p = 3.8 × 10-8), which was also associated with increased weight gain in the whole cohort (p = 0.004). Epigenome-wide association analyses failed to identify a significant link between metabolic changes and methylation data. Nevertheless, among the strongest associations, a potential causal effect of the baseline methylation level of cg11622362 on glycemia was revealed by a two-sample Mendelian randomization analysis (n = 3841 for instrument-exposure association; n = 314,916 for instrument-outcome association).

CONCLUSION

These findings provide new insights into the mechanisms of psychotropic drug-induced weight gain, revealing important epigenetic alterations upon treatment, some of which may play a mediatory role.

Genome-Wide Network Analysis of Above- and Below-Ground Co-growth in Populus euphratica

Tree growth is the consequence of developmental interactions between above- and below-ground compartments. However, a comprehensive view of the genetic architecture of growth as a cohesive whole is poorly understood. We propose a systems biology approach for mapping growth trajectories in genome-wide association studies viewing growth as a complex (phenotypic) system in which above- and below-ground components (or traits) interact with each other to mediate systems behavior. We further assume that trait-trait interactions are controlled by a genetic system composed of many different interactive genes and integrate the Lotka-Volterra predator-prey model to dissect phenotypic and genetic systems into pleiotropic and epistatic interaction components by which the detailed genetic mechanism of above- and below-ground co-growth can be charted. We apply the approach to analyze linkage mapping data of Populus euphratica, which is the only tree species that can grow in the desert, and characterize several loci that govern how above- and below-ground growth is cooperated or competed over development. We reconstruct multilayer and multiplex genetic interactome networks for the developmental trajectories of each trait and their developmental covariation. Many significant loci and epistatic effects detected can be annotated to candidate genes for growth and developmental processes. The results from our model may potentially be useful for marker-assisted selection and genetic editing in applied tree breeding programs. The model provides a general tool to characterize a complete picture of pleiotropic and epistatic genetic architecture in growth traits in forest trees and any other organisms.

Metabolic Syndrome: An Overview on Its Genetic Associations and Gene-Diet Interactions

Metabolic syndrome (MetS) is a cluster of cardiometabolic risk factors that includes central obesity, hyperglycemia, hypertension, and dyslipidemias and whose inter-related occurrence may increase the odds of developing type 2 diabetes and cardiovascular diseases. MetS has become one of the most studied conditions, nevertheless, due to its complex etiology, this has not been fully elucidated. Recent evidence describes that both genetic and environmental factors play an important role on its development. With the advent of genomic-wide association studies, single nucleotide polymorphisms (SNPs) have gained special importance. In this review, we present an update of the genetics surrounding MetS as a single entity as well as its corresponding risk factors, considering SNPs and gene-diet interactions related to cardiometabolic markers. In this study, we focus on the conceptual aspects, diagnostic criteria, as well as the role of genetics, particularly on SNPs and polygenic risk scores (PRS) for interindividual analysis. In addition, this review highlights future perspectives of personalized nutrition with regard to the approach of MetS and how individualized multiomics approaches could improve the current outlook.

A Meta-Analysis of the Genome-Wide Association Studies on Two Genetically Correlated Phenotypes Suggests Four New Risk Loci for Headaches

Headache is one of the commonest complaints that doctors need to address in clinical settings. The genetic mechanisms of different types of headache are not well understood while it has been suggested that self-reported headache and self-reported migraine were genetically correlated. In this study, we performed a meta-analysis of genome-wide association studies (GWAS) on the self-reported headache phenotype from the UK Biobank and the self-reported migraine phenotype from the 23andMe using the Unified Score-based Association Test (metaUSAT) software for genetically correlated phenotypes (N = 397,385). We identified 38 loci for headaches, of which 34 loci have been reported before and four loci were newly suggested. The LDL receptor related protein 1 (LRP1)-Signal Transducer and Activator of Transcription 6 (STAT6)-S hort chain D ehydrogenase/R eductase family 9C member 7 (SDR9C7) region in chromosome 12 was the most significantly associated locus with a leading p value of 1.24 × 10-62 of rs11172113. The One Cut homeobox 2 (ONECUT2) gene locus in chromosome 18 was the strongest signal among the four new loci with a p value of 1.29 × 10-9 of rs673939. Our study demonstrated that the genetically correlated phenotypes of self-reported headache and self-reported migraine can be meta-analysed together in theory and in practice to boost study power to identify more variants for headaches. This study has paved way for a large GWAS meta-analysis involving cohorts of different while genetically correlated headache phenotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-022-00078-7.

Insight into genetic, biological, and environmental determinants of sexual-dimorphism in type 2 diabetes and glucose-related traits

There is growing evidence that sex and gender differences play an important role in risk and pathophysiology of type 2 diabetes (T2D). Men develop T2D earlier than women, even though there is more obesity in young women than men. This difference in T2D prevalence is attenuated after the menopause. However, not all women are equally protected against T2D before the menopause, and gestational diabetes represents an important risk factor for future T2D. Biological mechanisms underlying sex and gender differences on T2D physiopathology are not yet fully understood. Sex hormones affect behavior and biological changes, and can have implications on lifestyle; thus, both sex-specific environmental and biological risk factors interact within a complex network to explain the differences in T2D risk and physiopathology in men and women. In addition, lifetime hormone fluctuations and body changes due to reproductive factors are generally more dramatic in women than men (ovarian cycle, pregnancy, and menopause). Progress in genetic studies and rodent models have significantly advanced our understanding of the biological pathways involved in the physiopathology of T2D. However, evidence of the sex-specific effects on genetic factors involved in T2D is still limited, and this gap of knowledge is even more important when investigating sex-specific differences during the life course. In this narrative review, we will focus on the current state of knowledge on the sex-specific effects of genetic factors associated with T2D over a lifetime, as well as the biological effects of these different hormonal stages on T2D risk. We will also discuss how biological insights from rodent models complement the genetic insights into the sex-dimorphism effects on T2D. Finally, we will suggest future directions to cover the knowledge gaps.

Genome-wide study on 72,298 individuals in Korean biobank data for 76 traits

Genome-wide association studies (GWAS) on diverse ancestry groups are lacking, resulting in deficits of genetic discoveries and polygenic scores. We conducted GWAS for 76 phenotypes in Korean biobank data, namely the Korean Genome and Epidemiology Study (KoGES) (n = 72,298). Our analysis discovered 2,242 associated loci, including 122 novel associations, many of which were replicated in Biobank Japan (BBJ) GWAS. We also applied several up-to-date methods for genetic association tests to increase the power, discovering additional associations that are not identified in simple case-control GWAS. We evaluated genetic pleiotropy to investigate genes associated with multiple traits. Following meta-analysis of 32 phenotypes between KoGES and BBJ, we further identified 379 novel associations and demonstrated the improved predictive performance of polygenic risk scores by using the meta-analysis results. The summary statistics of 76 KoGES GWAS phenotypes are publicly available, contributing to a better comprehension of the genetic architecture of the East Asian population.

Stratification of Type 2 Diabetes by Age of Diagnosis in the UK Biobank Reveals Subgroup-Specific Genetic Associations and Causal Risk Profiles

The pathogenesis of type 2 diabetes (T2D) might change with increasing age. Here, we used a stratification based on age of diagnosis to gain insight into the genetics and causal risk factors of T2D across different age-groups. We performed genome-wide association studies (GWAS) on T2D and T2D subgroups based on age of diagnosis (<50, 50-60, 60-70, and >70 years) (total of 24,986 cases). As control subjects, participants were at least 70 years of age at the end of follow-up without developing T2D (N =187,130). GWAS identified 208 independent lead single nucleotide polymorphism (SNPs) mapping to 69 loci associated with T2D (P < 1.0e-8). Among others, SNPs mapped to CDKN2B-AS1 and multiple independent SNPs mapped to TCF7L2 were more strongly associated with cases diagnosed after age 70 years than with cases diagnosed before age 50 years. Based on the different case groups, we performed two-sample Mendelian randomization. Most notably, we observed that of the investigated risk factors, the association between BMI and T2D attenuated with increasing age of diagnosis. Collectively, our results indicate that stratification of T2D based on age of diag-nosis reveals subgroup-specific genetics and causal determinants, supporting the hypothesis that the pathogenesis of T2D changes with increasing age.

Comparison of adaptive multiple phenotype association tests using summary statistics in genome-wide association studies

Genome-wide association studies have been successful mapping loci for individual phenotypes, but few studies have comprehensively interrogated evidence of shared genetic effects across multiple phenotypes simultaneously. Statistical methods have been proposed for analyzing multiple phenotypes using summary statistics, which enables studies of shared genetic effects while avoiding challenges associated with individual-level data sharing. Adaptive tests have been developed to maintain power against multiple alternative hypotheses because the most powerful single-alternative test depends on the underlying structure of the associations between the multiple phenotypes and a single nucleotide polymorphism (SNP). Here we compare the performance of six such adaptive tests: two adaptive sum of powered scores (aSPU) tests, the unified score association test (metaUSAT), the adaptive test in a mixed-models framework (mixAda) and two principal-component-based adaptive tests (PCAQ and PCO). Our simulations highlight practical challenges that arise when multivariate distributions of phenotypes do not satisfy assumptions of multivariate normality. Previous reports in this context focus on low minor allele count (MAC) and omit the aSPU test, which relies less than other methods on asymptotic and distributional assumptions. When these assumptions are not satisfied, particularly when MAC is low and/or phenotype covariance matrices are singular or nearly singular, aSPU better preserves type I error, sometimes at the cost of decreased power. We illustrate this trade-off with multiple phenotype analyses of six quantitative electrocardiogram traits in the Population Architecture using Genomics and Epidemiology (PAGE) study.

Lossless integration of multiple electronic health records for identifying pleiotropy using summary statistics

Increasingly, clinical phenotypes with matched genetic data from bio-bank linked electronic health records (EHRs) have been used for pleiotropy analyses. Thus far, pleiotropy analysis using individual-level EHR data has been limited to data from one site. However, it is desirable to integrate EHR data from multiple sites to improve the detection power and generalizability of the results. Due to privacy concerns, individual-level patients' data are not easily shared across institutions. As a result, we introduce Sum-Share, a method designed to efficiently integrate EHR and genetic data from multiple sites to perform pleiotropy analysis. Sum-Share requires only summary-level data and one round of communication from each site, yet it produces identical test statistics compared with that of pooled individual-level data. Consequently, Sum-Share can achieve lossless integration of multiple datasets. Using real EHR data from eMERGE, Sum-Share is able to identify 1734 potential pleiotropic SNPs for five cardiovascular diseases.