Search

Differences between sexes contribute to variation in the levels of fasting glucose and insulin. Epidemiological studies established a higher prevalence of impaired fasting glucose in men and impaired glucose tolerance in women, however, the genetic component underlying this phenomenon is not established. We assess sex-dimorphic (73,089/50,404 women and 67,506/47,806 men) and sex-combined (151,188/105,056 individuals) fasting glucose/fasting insulin genetic effects via genome-wide association study meta-analyses in individuals of European descent without diabetes. Here we report sex dimorphism in allelic effects on fasting insulin at IRS1 and ZNF12 loci, the latter showing higher RNA expression in whole blood in women compared to men. We also observe sex-homogeneous effects on fasting glucose at seven novel loci. Fasting insulin in women shows stronger genetic correlations than in men with waist-to-hip ratio and anorexia nervosa. Furthermore, waist-to-hip ratio is causally related to insulin resistance in women, but not in men. These results position dissection of metabolic and glycemic health sex dimorphism as a steppingstone for understanding differences in genetic effects between women and men in related phenotypes.

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

While nearly all common genomic variants associated with schizophrenia have no known function, one corresponds to a missense variant associated with change in efficiency of a metal ion transporter, ZIP8, coded by SLC39A8. This variant has been linked to a range of phenotypes and is believed to be under recent selection pressure, but its impact on health is poorly understood. We sought to understand phenotypic implications of this variant in a large genomic biobank using an unbiased phenome-wide approach. Specifically, we generated 50 topics based on diagnostic codes using latent Dirichlet allocation, and examined them for association with the risk variant. Then, any significant topics were further characterized by examining association with individual diagnostic codes contributing to the topic. Among 50 topics, 1 was associated at an experiment-wide significance threshold (beta = 0.003, uncorrected p = 0.00049), comprising predominantly brain-related codes, including intracranial hemorrhage, cerebrovascular disease, and delirium/dementia. These results suggest that a functional variant previously associated with schizophrenia risk also increases liability to cerebrovascular disease. They further illustrate the utility of a topic-based approach to phenome-wide association.

Event-related potential (ERP) components have been used to assess cognitive functions in patients with psychotic illness. Evidence suggests that among patients with psychosis there is a distinct heritable neurophysiologic phenotypic subtype captured by impairments across a range of ERP measures. In this study, we investigated the genetic basis of this "globally impaired" ERP cluster and its relationship to psychosis and cognitive abilities. We applied K-means clustering to six ERP measures to re-derive the globally impaired (n = 60) and the non-globally impaired ERP clusters (n = 323) in a sample of cases with schizophrenia (SCZ = 136) or bipolar disorder (BPD = 121) and healthy controls (n = 126). We used genome-wide association study (GWAS) results for SCZ, BPD, college completion, and childhood intelligence as the discovery datasets to derive polygenic risk scores (PRS) in our study sample and tested their associations with globally impaired ERP. We conducted mediation analyses to estimate the proportion of each PRS effect on severity of psychotic symptoms that is mediated through membership in the globally impaired ERP. Individuals with globally impaired ERP had significantly higher PANSS-positive scores (β = 3.95, P = 0.005). The SCZ-PRS was nominally associated with globally impaired ERP (unadjusted P = 0.01; R 2  = 3.07%). We also found a significant positive association between the college-PRS and globally impaired ERP (FDR-corrected P = 0.004; R 2  = 6.15%). The effect of college-PRS on PANSS positivity was almost entirely (97.1%) mediated through globally impaired ERP. These results suggest that the globally impaired ERP phenotype may represent some aspects of brain physiology on the path between genetic influences on educational attainment and psychotic symptoms.

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N=78 199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5 × 10 -8 ), we found suggestive evidence (P<5 × 10 -6 ) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.

Variants at microRNA-137 (MIR137), one of the most strongly associated schizophrenia risk loci identified to date, have been associated with poorer cognitive performance. As microRNA-137 is known to regulate the expression of ~1900 other genes, including several that are independently associated with schizophrenia, we tested whether this gene set was also associated with variation in cognitive performance. Our analysis was based on an empirically derived list of genes whose expression was altered by manipulation of MIR137 expression. This list was cross-referenced with genome-wide schizophrenia association data to construct individual polygenic scores. We then tested, in a sample of 808 patients and 192 controls, whether these risk scores were associated with altered performance on cognitive functions known to be affected in schizophrenia. A subgroup of healthy participants also underwent functional imaging during memory (n=108) and face processing tasks (n=83). Increased polygenic risk within the empirically derived miR-137 regulated gene score was associated with significantly lower performance on intelligence quotient, working memory and episodic memory. These effects were observed most clearly at a polygenic threshold of P=0.05, although significant results were observed at all three thresholds analyzed. This association was found independently for the gene set as a whole, excluding the schizophrenia-associated MIR137 SNP itself. Analysis of the spatial working memory fMRI task further suggested that increased risk score (thresholded at P=10 -5 ) was significantly associated with increased activation of the right inferior occipital gyrus. In conclusion, these data are consistent with emerging evidence that MIR137 associated risk for schizophrenia may relate to its broader downstream genetic effects.

We performed a genome-wide association study of 6447 bipolar disorder (BD) cases and 12 639 controls from the International Cohort Collection for Bipolar Disorder (ICCBD). Meta-analysis was performed with prior results from the Psychiatric Genomics Consortium Bipolar Disorder Working Group for a combined sample of 13 902 cases and 19 279 controls. We identified eight genome-wide significant, associated regions, including a novel associated region on chromosome 10 (rs10884920; P=3.28 × 10 -8 ) that includes the brain-enriched cytoskeleton protein adducin 3 (ADD3), a non-coding RNA, and a neuropeptide-specific aminopeptidase P (XPNPEP1). Our large sample size allowed us to test the heritability and genetic correlation of BD subtypes and investigate their genetic overlap with schizophrenia and major depressive disorder. We found a significant difference in heritability of the two most common forms of BD (BD I SNP-h 2 =0.35; BD II SNP-h 2 =0.25; P=0.02). The genetic correlation between BD I and BD II was 0.78, whereas the genetic correlation was 0.97 when BD cohorts containing both types were compared. In addition, we demonstrated a significantly greater load of polygenic risk alleles for schizophrenia and BD in patients with BD I compared with patients with BD II, and a greater load of schizophrenia risk alleles in patients with the bipolar type of schizoaffective disorder compared with patients with either BD I or BD II. These results point to a partial difference in the genetic architecture of BD subtypes as currently defined.

The ability to predict psychiatric readmission would facilitate the development of interventions to reduce this risk, a major driver of psychiatric health-care costs. The symptoms or characteristics of illness course necessary to develop reliable predictors are not available in coded billing data, but may be present in narrative electronic health record (EHR) discharge summaries. We identified a cohort of individuals admitted to a psychiatric inpatient unit between 1994 and 2012 with a principal diagnosis of major depressive disorder, and extracted inpatient psychiatric discharge narrative notes. Using these data, we trained a 75-topic Latent Dirichlet Allocation (LDA) model, a form of natural language processing, which identifies groups of words associated with topics discussed in a document collection. The cohort was randomly split to derive a training (70%) and testing (30%) data set, and we trained separate support vector machine models for baseline clinical features alone, baseline features plus common individual words and the above plus topics identified from the 75-topic LDA model. Of 4687 patients with inpatient discharge summaries, 470 were readmitted within 30 days. The 75-topic LDA model included topics linked to psychiatric symptoms (suicide, severe depression, anxiety, trauma, eating/weight and panic) and major depressive disorder comorbidities (infection, postpartum, brain tumor, diarrhea and pulmonary disease). By including LDA topics, prediction of readmission, as measured by area under receiver-operating characteristic curves in the testing data set, was improved from baseline (area under the curve 0.618) to baseline+1000 words (0.682) to baseline+75 topics (0.784). Inclusion of topics derived from narrative notes allows more accurate discrimination of individuals at high risk for psychiatric readmission in this cohort. Topic modeling and related approaches offer the potential to improve prediction using EHRs, if generalizability can be established in other clinical cohorts.

Multiple studies have examined the risk of prenatal antidepressant exposure and risk for autism spectrum disorder (ASD) or attention-deficit hyperactivity disorder (ADHD), with inconsistent results. Precisely estimating such risk, if any, is of great importance in light of the need to balance such risk with the benefit of depression and anxiety treatment. We developed a method to integrate data from multiple New England health systems, matching offspring and maternal health data in electronic health records to characterize diagnoses and medication exposure. Children with ASD or ADHD were matched 1:3 with children without neurodevelopmental disorders. Association between maternal antidepressant exposure and ASD or ADHD liability was examined using logistic regression, adjusting for potential sociodemographic and psychiatric confounding variables. In new cohorts of 1245 ASD cases and 1701 ADHD cases, along with age-, sex- and socioeconomic status matched controls, neither disorder was significantly associated with prenatal antidepressant exposure in crude or adjusted models (adjusted odds ratio 0.90, 95% confidence interval 0.50-1.54 for ASD; 0.97, 95% confidence interval 0.53-1.69 for ADHD). Pre-pregnancy antidepressant exposure significantly increased risk for both disorders. These results suggest that prior reports of association between prenatal antidepressant exposure and neurodevelopmental disease are likely to represent a false-positive finding, which may arise in part through confounding by indication. They further demonstrate the potential to integrate data across electronic health records studies spanning multiple health systems to enable efficient pharmacovigilance investigation., Competing Interests: Dr Perlis has received consulting fees or served on scientific advisory boards for Genomind, Healthrageous, Pamlab, Perfect Health, Pfizer, Proteus Biomedical, Psybrain and RID Ventures, and received patent fees/royalties from Concordant Rater Systems (now UBC/Medco). The remaining authors declare no conflict of interest.

Variability in response to drug use is common and heritable, suggesting that genome-wide pharmacogenomics studies may help explain the 'missing heritability' of complex traits. Here, we describe four independent analyses in 33 781 participants of European ancestry from 10 cohorts that were designed to identify genetic variants modifying the effects of drugs on QT interval duration (QT). Each analysis cross-sectionally examined four therapeutic classes: thiazide diuretics (prevalence of use=13.0%), tri/tetracyclic antidepressants (2.6%), sulfonylurea hypoglycemic agents (2.9%) and QT-prolonging drugs as classified by the University of Arizona Center for Education and Research on Therapeutics (4.4%). Drug-gene interactions were estimated using covariable-adjusted linear regression and results were combined with fixed-effects meta-analysis. Although drug-single-nucleotide polymorphism (SNP) interactions were biologically plausible and variables were well-measured, findings from the four cross-sectional meta-analyses were null (Pinteraction>5.0 × 10(-8)). Simulations suggested that additional efforts, including longitudinal modeling to increase statistical power, are likely needed to identify potentially important pharmacogenomic effects.

Major depressive disorder (MDD) is a common psychiatric illness characterized by low mood and loss of interest in pleasurable activities. Despite years of effort, recent genome-wide association studies (GWAS) have identified few susceptibility variants or genes that are robustly associated with MDD. Standard single-SNP (single nucleotide polymorphism)-based GWAS analysis typically has limited power to deal with the extensive heterogeneity and substantial polygenic contribution of individually weak genetic effects underlying the pathogenesis of MDD. Here, we report an alternative, gene-set-based association analysis of MDD in an effort to identify groups of biologically related genetic variants that are involved in the same molecular function or cellular processes and exhibit a significant level of aggregated association with MDD. In particular, we used a text-mining-based data analysis to prioritize candidate gene sets implicated in MDD and conducted a multi-locus association analysis to look for enriched signals of nominally associated MDD susceptibility loci within each of the gene sets. Our primary analysis is based on the meta-analysis of three large MDD GWAS data sets (total N=4346 cases and 4430 controls). After correction for multiple testing, we found that genes involved in glutamatergic synaptic neurotransmission were significantly associated with MDD (set-based association P=6.9 × 10(-4)). This result is consistent with previous studies that support a role of the glutamatergic system in synaptic plasticity and MDD and support the potential utility of targeting glutamatergic neurotransmission in the treatment of MDD.

Our understanding of the genetics of type 2 diabetes mellitus (T2DM) has changed, in part owing to implementation of genome-wide association studies as a method for unraveling the genetic architecture of complex traits. These studies enable a global search throughout the nuclear genome for variants that are associated with specific phenotypes. Currently, single nucleotide polymorphisms in about 24 different genetic loci have been associated with T2DM. Most of these genetic loci are associated with the insulin secretion pathway rather than insulin resistance. Study design, heritability differences and the intrinsic properties of in vivo insulin resistance measures might partially explain why only a few loci associated with insulin resistance have been detected through genome-wide association approaches. Despite the success of these approaches at detecting loci associated with T2DM, currently known associations explain only a small amount of the genetic variance involved in the disease. Compared with previous studies, larger cohorts might be needed to identify variants of smaller effect sizes and lower allele frequencies. Finally, the current list of genetic loci that are related to T2DM does not seem to offer greater predictive value in determining diabetes risk than do commonly used phenotypic risk factors and family history.