A framework for employing longitudinally collected multicenter electronic health records to stratify heterogeneous patient populations on disease history

Objective To facilitate patient disease subset and risk factor identification by constructing a pipeline which is generalizable, provides easily interpretable results, and allows replication by overcoming electronic health records (EHRs) batch effects. Material and Methods We used 1872 billing codes in EHRs of 102 880 patients from 12 healthcare systems. Using tools borrowed from single-cell omics, we mitigated center-specific batch effects and performed clustering to identify patients with highly similar medical history patterns across the various centers. Our visualization method (PheSpec) depicts the phenotypic profile of clusters, applies a novel filtering of noninformative codes (Ranked Scope Pervasion), and indicates the most distinguishing features. Results We observed 114 clinically meaningful profiles, for example, linking prostate hyperplasia with cancer and diabetes with cardiovascular problems and grouping pediatric developmental disorders. Our framework identified disease subsets, exemplified by 6 "other headache" clusters, where phenotypic profiles suggested different underlying mechanisms: migraine, convulsion, injury, eye problems, joint pain, and pituitary gland disorders. Phenotypic patterns replicated well, with high correlations of >= 0.75 to an average of 6 (2-8) of the 12 different cohorts, demonstrating the consistency with which our method discovers disease history profiles. Discussion Costly clinical research ventures should be based on solid hypotheses. We repurpose methods from single-cell omics to build these hypotheses from observational EHR data, distilling useful information from complex data. Conclusion We establish a generalizable pipeline for the identification and replication of clinically meaningful (sub)phenotypes from widely available high-dimensional billing codes. This approach overcomes datatype problems and produces comprehensive visualizations of validation-ready phenotypes.Molecular Epidemiolog

Psychiatric genomics: An update and an Agenda

The Psychiatric Genomics Consortium (PGC) is the largest consortium in the history of psychiatry. This global effort is dedicated to rapid progress and open science, and in the past decade it has delivered an increasing flow of new knowledge about the fundamental basis of common psychiatric disorders. The PGC has recently commenced a program of research designed to deliver “actionable” findings—genomic results that 1) reveal fundamental biology, 2) inform clinical practice, and 3) deliver new therapeutic targets. The central idea of the PGC is to convert the family history risk factor into biologically, clinically, and therapeutically meaningful insights. The emerging findings suggest that we are entering a phase of accelerated genetic discovery for multiple psychiatric disorders. These findings are likely to elucidate the genetic portions of these truly complex traits, and this knowledge can then be mined for its relevance for improved therapeutics and its impact on psychiatric practice within a precision medicine framework

Genome-wide association for major depressive disorder: a possible role for the presynaptic protein piccolo

Major depressive disorder (MDD) is a common complex trait with enormous public health significance. As part of the Genetic Association Information Network initiative of the US Foundation for the National Institutes of Health, we conducted a genome-wide association study of 435 291 single nucleotide polymorphisms (SNPs) genotyped in 1738 MDD cases and 1802 controls selected to be at low liability for MDD. Of the top 200, 11 signals localized to a 167 kb region overlapping the gene piccolo (PCLO, whose protein product localizes to the cytomatrix of the presynaptic active zone and is important in monoaminergic neurotransmission in the brain) with P-values of 7.7 × 1

Largest GWAS of PTSD (N=20 070) yields genetic overlap with schizophrenia and sex differences in heritability

The Psychiatric Genomics Consortium-Posttraumatic Stress Disorder group (PGC-PTSD) combined genome-wide case-control molecular genetic data across 11 multiethnic studies to quantify PTSD heritability, to examine potential shared genetic risk with schizophrenia, bipolar disorder, and major depressive disorder and to identify risk loci for PTSD. Examining 20 730 individuals, we report a molecular genetics-based heritability estimate (h 2 SNP) for European-American females of 29% that is similar to h 2 SNP for schizophrenia and is substantially higher than h 2 SNP in European-American males (estimate not distinguishable from zero). We found strong evidence of overlapping genetic risk between PTSD and schizophrenia along with more modest evidence of overlap with bipolar and major depressive disorder. No single-nucleotide polymorphisms (SNPs) exceeded genome-wide significance in the transethnic (overall) meta-analysis and we do not replicate previously reported associations. Still, SNP-level summary statistics made available here afford the best-available molecular genetic index of PTSD - for both European- and African-American individuals - and can be used in polygenic risk prediction and genetic correlation studies of diverse phenotypes. Publication of summary statistics for 1/410 000 African Americans contributes to the broader goal of increased ancestral diversity in genomic data resources. In sum, the results demonstrate genetic influences on the development of PTSD, identify shared genetic risk between PTSD and other psychiatric disorders and highlight the importance of multiethnic/racial samples. As has been the case with schizophrenia and other complex genetic disorders, larger sample sizes are needed to identify specific risk loci

A putative causal relationship between genetically determined female body shape and posttraumatic stress disorder

Background: The nature and underlying mechanisms of the observed increased vulnerability to posttraumatic stress disorder (PTSD) in women are unclear. Methods: We investigated the genetic overlap of PTSD with anthropometric traits and reproductive behaviors and functions in women. The analysis was conducted using female-specific summary statistics from large genome-wide association studies (GWAS) and a cohort of 3577 European American women (966 PTSD cases and 2611 trauma-exposed controls). We applied a high-resolution polygenic score approach and Mendelian randomization analysis to investigate genetic correlations and causal relationships. Results: We observed an inverse association of PTSD with genetically determined anthropometric traits related to body shape, independent of body mass index (BMI). The top association was related to BMI-adjusted waist circumference (WCadj; R = -0.079, P < 0.001, Q = 0.011). We estimated a relative decrease of 64.6% (95% confidence interval = 27.5-82.7) in the risk of PTSD per 1-SD increase in WCadj. MR-Egger regression intercept analysis showed no evidence of pleiotropic effects in this association (Ppleiotropy = 0.979). We also observed associations of genetically determined WCadj with age at first sexual intercourse and number of sexual partners (P = 0.013 and P < 0.001, respectively). Conclusions: There is a putative causal relationship between genetically determined female body shape and PTSD, which could be mediated by evolutionary mechanisms involved in human sexual behaviors

Serum magnesium and calcium levels in relation to ischemic stroke : Mendelian randomization study

ObjectiveTo determine whether serum magnesium and calcium concentrations are causally associated with ischemic stroke or any of its subtypes using the mendelian randomization approach.MethodsAnalyses were conducted using summary statistics data for 13 single-nucleotide polymorphisms robustly associated with serum magnesium (n = 6) or serum calcium (n = 7) concentrations. The corresponding data for ischemic stroke were obtained from the MEGASTROKE consortium (34,217 cases and 404,630 noncases).ResultsIn standard mendelian randomization analysis, the odds ratios for each 0.1 mmol/L (about 1 SD) increase in genetically predicted serum magnesium concentrations were 0.78 (95% confidence interval [CI] 0.69-0.89; p = 1.3 7 10-4) for all ischemic stroke, 0.63 (95% CI 0.50-0.80; p = 1.6 7 10-4) for cardioembolic stroke, and 0.60 (95% CI 0.44-0.82; p = 0.001) for large artery stroke; there was no association with small vessel stroke (odds ratio 0.90, 95% CI 0.67-1.20; p = 0.46). Only the association with cardioembolic stroke was robust in sensitivity analyses. There was no association of genetically predicted serum calcium concentrations with all ischemic stroke (per 0.5 mg/dL [about 1 SD] increase in serum calcium: odds ratio 1.03, 95% CI 0.88-1.21) or with any subtype.ConclusionsThis study found that genetically higher serum magnesium concentrations are associated with a reduced risk of cardioembolic stroke but found no significant association of genetically higher serum calcium concentrations with any ischemic stroke subtype

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.Peer reviewe

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg =-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness

Genome-wide association study identifies 30 Loci Associated with Bipolar Disorder

This paper is dedicated to the memory of Psychiatric Genomics Consortium (PGC) founding member and Bipolar disorder working group co-chair Pamela Sklar. We thank the participants who donated their time, experiences and DNA to this research, and to the clinical and scientific teams that worked with them. We are deeply indebted to the investigators who comprise the PGC. The views expressed are those of the authors and not necessarily those of any funding or regulatory body. Analyses were carried out on the NL Genetic Cluster Computer (http://www.geneticcluster.org ) hosted by SURFsara, and the Mount Sinai high performance computing cluster (http://hpc.mssm.edu).Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P<1x10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (GWS, p < 5x10-8) in the discovery GWAS were not GWS in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis 30 loci were GWS including 20 novel loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene-sets including regulation of insulin secretion and endocannabinoid signaling. BDI is strongly genetically correlated with schizophrenia, driven by psychosis, whereas BDII is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential new biological mechanisms for BD.This work was funded in part by the Brain and Behavior Research Foundation, Stanley Medical Research Institute, University of Michigan, Pritzker Neuropsychiatric Disorders Research Fund L.L.C., Marriot Foundation and the Mayo Clinic Center for Individualized Medicine, the NIMH Intramural Research Program; Canadian Institutes of Health Research; the UK Maudsley NHS Foundation Trust, NIHR, NRS, MRC, Wellcome Trust; European Research Council; German Ministry for Education and Research, German Research Foundation IZKF of Münster, Deutsche Forschungsgemeinschaft, ImmunoSensation, the Dr. Lisa-Oehler Foundation, University of Bonn; the Swiss National Science Foundation; French Foundation FondaMental and ANR; Spanish Ministerio de Economía, CIBERSAM, Industria y Competitividad, European Regional Development Fund (ERDF), Generalitat de Catalunya, EU Horizon 2020 Research and Innovation Programme; BBMRI-NL; South-East Norway Regional Health Authority and Mrs. Throne-Holst; Swedish Research Council, Stockholm County Council, Söderström Foundation; Lundbeck Foundation, Aarhus University; Australia NHMRC, NSW Ministry of Health, Janette M O'Neil and Betty C Lynch

The genetic architecture of the human cerebral cortex

INTRODUCTION The cerebral cortex underlies our complex cognitive capabilities. Variations in human cortical surface area and thickness are associated with neurological, psychological, and behavioral traits and can be measured in vivo by magnetic resonance imaging (MRI). Studies in model organisms have identified genes that influence cortical structure, but little is known about common genetic variants that affect human cortical structure. RATIONALE To identify genetic variants associated with human cortical structure at both global and regional levels, we conducted a genome-wide association meta-analysis of brain MRI data from 51,665 individuals across 60 cohorts. We analyzed the surface area and average thickness of the whole cortex and 34 cortical regions with known functional specializations. RESULTS We identified 306 nominally genome-wide significant loci (P < 5 × 10−8) associated with cortical structure in a discovery sample of 33,992 participants of European ancestry. Of the 299 loci for which replication data were available, 241 loci influencing surface area and 14 influencing thickness remained significant after replication, with 199 loci passing multiple testing correction (P < 8.3 × 10−10; 187 influencing surface area and 12 influencing thickness). Common genetic variants explained 34% (SE = 3%) of the variation in total surface area and 26% (SE = 2%) in average thickness; surface area and thickness showed a negative genetic correlation (rG = −0.32, SE = 0.05, P = 6.5 × 10−12), which suggests that genetic influences have opposing effects on surface area and thickness. Bioinformatic analyses showed that total surface area is influenced by genetic variants that alter gene regulatory activity in neural progenitor cells during fetal development. By contrast, average thickness is influenced by active regulatory elements in adult brain samples, which may reflect processes that occur after mid-fetal development, such as myelination, branching, or pruning. When considered together, these results support the radial unit hypothesis that different developmental mechanisms promote surface area expansion and increases in thickness. To identify specific genetic influences on individual cortical regions, we controlled for global measures (total surface area or average thickness) in the regional analyses. After multiple testing correction, we identified 175 loci that influence regional surface area and 10 that influence regional thickness. Loci that affect regional surface area cluster near genes involved in the Wnt signaling pathway, which is known to influence areal identity. We observed significant positive genetic correlations and evidence of bidirectional causation of total surface area with both general cognitive functioning and educational attainment. We found additional positive genetic correlations between total surface area and Parkinson’s disease but did not find evidence of causation. Negative genetic correlations were evident between total surface area and insomnia, attention deficit hyperactivity disorder, depressive symptoms, major depressive disorder, and neuroticism. CONCLUSION This large-scale collaborative work enhances our understanding of the genetic architecture of the human cerebral cortex and its regional patterning. The highly polygenic architecture of the cortex suggests that distinct genes are involved in the development of specific cortical areas. Moreover, we find evidence that brain structure is a key phenotype along the causal pathway that leads from genetic variation to differences in general cognitive function