Pleiotropic Meta-Analysis of Cognition, Education, and Schizophrenia Differentiates Roles of Early Neurodevelopmental and Adult Synaptic Pathways

Cooperative Studies Program (CSP) #572: A Study of Serious Mental Illness in Veterans as a Pathway to personalized medicine in Schizophrenia and Bipolar Illness

BACKGROUND

Personalization of psychiatric treatment includes treatment of symptoms, cognition and functional deficits, suicide, and medical co-morbidities. VA Collaborative Study 572 examined a large sample of male and female veterans with schizophrenia (n=3,942) and with bipolar disorder (n=5,414) with phenotyping and genomic analyses. We present the results to date and future directions.

METHODS

All veterans received a structured diagnostic interview and assessments of suicidal ideation and behavior, PTSD, and health. Veterans with schizophrenia were assessed for negative symptoms and lifetime depression. All were assessed with a cognitive and functional capacity assessment. Data for genome wide association studies were collected. Controls came from the VA Million Veteran Program.

RESULTS

Suicidal ideation or behavior was present in 66%. Cognitive and functional deficits were consistent with previous studies. 40% of the veterans with schizophrenia had a lifetime major depressive episode and PTSD was present in over 30%. Polygenic risk score (PRS) analyses indicated that cognitive and functional deficits overlapped with PRS for cognition, education, and intelligence in the general population and PRS for suicidal ideation and behavior correlated with previous PRS for depression and suicidal ideation and behavior, as did the PRS for PTSD.

DISCUSSION

Results to date provide directions for personalization of treatment in SMI, veterans with SMI, and veterans in general. The results of the genomic analyses suggest that cognitive deficits in SMI may be associated with general population features. Upcoming genomic analyses will reexamine the issues above, as well as genomic factors associated with smoking, substance abuse, negative symptoms, and treatment response.

Genome-Wide Association Study of Motor Coordination

The ability to finely control our movement is key to achieving many of the educational milestones and life-skills we develop throughout our lives. Despite the centrality of coordination to early development, there is a vast gap in our understanding of the underlying biology. Like most complex traits, both genetics and environment influence motor coordination, however, the specific genes, early environmental risk factors and molecular pathways are unknown. Previous studies have shown that about 5% of school-age children experience unexplained difficulties with motor coordination. These children are said to have Developmental Coordination Disorder (DCD). For children with DCD, these motor coordination difficulties significantly impact their everyday life and learning. DCD is associated with poorer academic achievement, reduced quality of life, it can constrain career opportunities and increase the risk of mental health issues in adulthood. Despite the high prevalence of coordination difficulties, many children remain undiagnosed by healthcare professionals. Compounding under-diagnosis in the clinic, research into the etiology of DCD is severely underrepresented in the literature. Here we present the first genome-wide association study to examine the genetic basis of early motor coordination in the context of motor difficulties. Using data from the Avon Longitudinal Study of Parents and Children we generate a derived measure of motor coordination from four components of the Movement Assessment Battery for Children, providing an overall measure of coordination across the full range of ability. We perform the first genome-wide association analysis focused on motor coordination (N = 4542). No single nucleotide polymorphisms (SNPs) met the threshold for genome-wide significance, however, 59 SNPs showed suggestive associations. Three regions contained multiple suggestively associated SNPs, within five preliminary candidate genes: IQSEC1, LRCC1, SYNJ2B2, ADAM20, and ADAM21. Association to the gene IQSEC1 suggests a potential link to axon guidance and dendritic projection processes as a potential underlying mechanism of motor coordination difficulties. This represents an interesting potential mechanism, and whilst further validation is essential, it generates a direct window into the biology of motor coordination difficulties. This research has identified potential biological drivers of DCD, a first step towards understanding this common, yet neglected neurodevelopmental disorder.

Transcriptome-wide association analysis of brain structures yields insights into pleiotropy with complex neuropsychiatric traits

Structural variations of the human brain are heritable and highly polygenic traits, with hundreds of associated genes identified in recent genome-wide association studies (GWAS). Transcriptome-wide association studies (TWAS) can both prioritize these GWAS findings and also identify additional gene-trait associations. Here we perform cross-tissue TWAS analysis of 211 structural neuroimaging and discover 278 associated genes exceeding Bonferroni significance threshold of 1.04 × 10-8. The TWAS-significant genes for brain structures have been linked to a wide range of complex traits in different domains. Through TWAS gene-based polygenic risk scores (PRS) prediction, we find that TWAS PRS gains substantial power in association analysis compared to conventional variant-based GWAS PRS, and up to 6.97% of phenotypic variance (p-value = 7.56 × 10-31) can be explained in independent testing data sets. In conclusion, our study illustrates that TWAS can be a powerful supplement to traditional GWAS in imaging genetics studies for gene discovery-validation, genetic co-architecture analysis, and polygenic risk prediction.

Attachment Style Moderates Polygenic Risk for Posttraumatic Stress in United States Military Veterans: Results From the National Health and Resilience in Veterans Study

BACKGROUND

A polygenic risk score (PRS) derived from genome-wide association studies of posttraumatic stress disorder (PTSD) may inform risk for this disorder. To date, however, no known study has examined whether social environmental factors such as attachment style may moderate the relation between PRS and PTSD.

METHODS

We evaluated main and interactive effects of PRS and attachment style on PTSD symptoms in a nationally representative sample of trauma-exposed European-American U.S. military veterans (N = 2030). PRS was derived from a genome-wide association study of PTSD re-experiencing symptoms (N = 146,660) in the Million Veteran Program cohort. Using one-sample Mendelian randomization with data from the UK Biobank (N = 115,099), we evaluated the effects of re-experiencing PRS and attachment style on PTSD symptoms.

RESULTS

Higher re-experiencing PRS and secure attachment style were independently associated with PTSD symptoms. A significant PRS-by-attachment style interaction was also observed (β = -.11, p = .006), with a positive association between re-experiencing PRS and PTSD symptoms observed only among veterans with an insecure attachment style. One-sample Mendelian randomization analyses suggested that the association between PTSD symptoms and attachment style is bidirectional. PRS enrichment analyses revealed a significant interaction between attachment style and a variant mapping to the IGSF11 gene (rs151177743, p = 2.1 × 10^-7), which is implicated in regulating excitatory synaptic transmission and plasticity.

CONCLUSIONS

Attachment style may moderate polygenic risk for PTSD symptoms, and a novel locus implicated in synaptic transmission and plasticity may serve as a possible biological mediator of this association. These findings may help inform interpersonally oriented treatments for PTSD for individuals with high polygenic risk for this disorder.

Identifying loci with different allele frequencies among cases of eight psychiatric disorders using CC-GWAS

Psychiatric disorders are highly genetically correlated, but little research has been conducted on the genetic differences between disorders. We developed a new method (case-case genome-wide association study; CC-GWAS) to test for differences in allele frequency between cases of two disorders using summary statistics from the respective case-control GWAS, transcending current methods that require individual-level data. Simulations and analytical computations confirm that CC-GWAS is well powered with effective control of type I error. We applied CC-GWAS to publicly available summary statistics for schizophrenia, bipolar disorder, major depressive disorder and five other psychiatric disorders. CC-GWAS identified 196 independent case-case loci, including 72 CC-GWAS-specific loci that were not significant at the genome-wide level in the input case-control summary statistics; two of the CC-GWAS-specific loci implicate the genes KLF6 and KLF16 (from the Krüppel-like family of transcription factors), which have been linked to neurite outgrowth and axon regeneration. CC-GWAS loci replicated convincingly in applications to datasets with independent replication data.

A genome-wide association study identifies a gene network associated with paranoid schizophrenia and antipsychotics-induced tardive dyskinesia

In the present study we conducted a genome-wide association study (GWAS) in a cohort of 505 patients with paranoid schizophrenia (SCZ), of which 95 had tardive dyskinesia (TD), and 503 healthy controls. Using data generated by the PsychENCODE Consortium (PEC) and other bioinformatic databases, we revealed a gene network, implicated in neurodevelopment and brain function, associated with both these disorders. Almost all these genes are in gene or isoform co-expression PEC network modules important for the functioning of the brain; the activity of these networks is also altered in SCZ, bipolar disorder and autism spectrum disorders. The associated PEC network modules are enriched for gene ontology terms relevant to the brain development and function (CNS development, neuron development, axon ensheathment, synapse, synaptic vesicle cycle, and signaling receptor activity) and to the immune system (inflammatory response). Results of the present study suggest that orofacial and limbtruncal types of TD seem to share the molecular network with SCZ. Paranoid SCZ and abnormal involuntary movements that indicate the orofacial type of TD are associated with the same genomic loci on chromosomes 3p22.2, 8q21.13, and 13q14.2. The limbtruncal type of TD is associated with a locus on chromosome 3p13 where the best functional candidate is FOXP1, a high-confidence SCZ gene. The results of this study shed light on common pathogenic mechanisms for SCZ and TD, and indicate that the pathogenesis of the orofacial and limbtruncal types of TD might be driven by interacting genes implicated in neurodevelopment.

Gene-based association analysis identifies 190 genes affecting neuroticism

Neuroticism is a personality trait, which is an important risk factor for psychiatric disorders. Recent genome-wide studies reported about 600 genes potentially influencing neuroticism. Little is known about the mechanisms of their action. Here, we aimed to conduct a more detailed analysis of genes that can regulate the level of neuroticism. Using UK Biobank-based GWAS summary statistics, we performed a gene-based association analysis using four sets of within-gene variants, each set possessing specific protein-coding properties. To guard against the influence of strong GWAS signals outside the gene, we used a specially designed procedure called “polygene pruning”. As a result, we identified 190 genes associated with neuroticism due to the effect of within-gene variants rather than strong GWAS signals outside the gene. Thirty eight of these genes are new. Within all genes identified, we distinguished two slightly overlapping groups obtained from using protein-coding and non-coding variants. Many genes in the former group included potentially pathogenic variants. For some genes in the latter group, we found evidence of pleiotropy with gene expression. Using a bioinformatics analysis, we prioritized the neuroticism genes and showed that the genes that contribute to neuroticism through their within-gene variants are the most appropriate candidate genes.

Integrative analysis of genome-wide association studies identifies novel loci associated with neuropsychiatric disorders

Neuropsychiatric disorders, such as autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), schizophrenia (SCZ), bipolar disorder (BIP), and major depressive disorder (MDD) share common clinical presentations, suggesting etiologic overlap. A substantial proportion of SNP-based heritability for neuropsychiatric disorders is attributable to genetic components, and genome-wide association studies (GWASs) focusing on individual diseases have identified multiple genetic loci shared between these diseases. Here, we aimed at identifying novel genetic loci associated with individual neuropsychiatric diseases and genetic loci shared by neuropsychiatric diseases. We performed multi-trait joint analyses and meta-analysis across five neuropsychiatric disorders based on their summary statistics from the Psychiatric Genomics Consortium (PGC), and further carried out a replication study of ADHD among 2726 cases and 16299 controls in an independent pediatric cohort. In the multi-trait joint analyses, we found five novel genome-wide significant loci for ADHD, one novel locus for BIP, and ten novel loci for MDD. We further achieved modest replication in our independent pediatric dataset. We conducted fine-mapping and functional annotation through an integrative multi-omics approach and identified causal variants and potential target genes at each novel locus. Gene expression profile and gene-set enrichment analysis further suggested early developmental stage expression pattern and postsynaptic membrane compartment enrichment of candidate genes at the genome-wide significant loci of these neuropsychiatric disorders. Therefore, through a multi-omics approach, we identified novel genetic loci associated with the five neuropsychiatric disorders which may help to better understand the underlying molecular mechanism of neuropsychiatric diseases.

Insights into the Molecular Mechanisms of Liuwei Dihuang Decoction via Network Pharmacology

The traditional Chinese medicines (TCMs) have been used to treat diseases over a long history, but it is still a great challenge to uncover the underlying mechanisms for their therapeutic effects due to the complexity of their ingredients. Based on a novel network pharmacology-based approach, we explored in this study the potential therapeutic targets of Liuwei Dihuang (LWDH) decoction in its neuroendocrine immunomodulation (NIM) function. We not only collected the known targets of the compounds in LWDH but also predicted the targets for these compounds using the balanced substructure-drug-target network-based inference (bSDTNBI), which is a target prediction method based on network inferring developed by our laboratory. A "target-(pathway)-target" (TPT) network, in which targets of LWDH were connected by relevant pathways, was constructed and divided into several separate modules with strong internal connections. Then the target module that contributes the most to NIM function was determined through a contribution scoring algorithm. Finally, the targets with the highest contribution score to NIM-related diseases in this target module were recommended as potential therapeutic targets of LWDH.

Genome-Wide Association Studies of Schizophrenia and Bipolar Disorder in a Diverse Cohort of US Veterans

BACKGROUND

Schizophrenia (SCZ) and bipolar disorder (BIP) are debilitating neuropsychiatric disorders, collectively affecting 2% of the world's population. Recognizing the major impact of these psychiatric disorders on the psychosocial function of more than 200 000 US Veterans, the Department of Veterans Affairs (VA) recently completed genotyping of more than 8000 veterans with SCZ and BIP in the Cooperative Studies Program (CSP) #572.

METHODS

We performed genome-wide association studies (GWAS) in CSP #572 and benchmarked the predictive value of polygenic risk scores (PRS) constructed from published findings. We combined our results with available summary statistics from several recent GWAS, realizing the largest and most diverse studies of these disorders to date.

RESULTS

Our primary GWAS uncovered new associations between CHD7 variants and SCZ, and novel BIP associations with variants in Sortilin Related VPS10 Domain Containing Receptor 3 (SORCS3) and downstream of PCDH11X. Combining our results with published summary statistics for SCZ yielded 39 novel susceptibility loci including CRHR1, and we identified 10 additional findings for BIP (28 326 cases and 90 570 controls). PRS trained on published GWAS were significantly associated with case-control status among European American (P < 10-30) and African American (P < .0005) participants in CSP #572.

CONCLUSIONS

We have demonstrated that published findings for SCZ and BIP are robustly generalizable to a diverse cohort of US veterans. Leveraging available summary statistics from GWAS of global populations, we report 52 new susceptibility loci and improved fine-mapping resolution for dozens of previously reported associations.

The enigmatic fetal subplate compartment forms an early tangential cortical nexus and provides the framework for construction of cortical connectivity

The most prominent transient compartment of the primate fetal cortex is the deep, cell-sparse, synapse-containing subplate compartment (SPC). The developmental role of the SPC and its extraordinary size in humans remain enigmatic. This paper evaluates evidence on the development and connectivity of the SPC and discusses its role in the pathogenesis of neurodevelopmental disorders. A synthesis of data shows that the subplate becomes a prominent compartment by its expansion from the deep cortical plate (CP), appearing well-delineated on MR scans and forming a tangential nexus across the hemisphere, consisting of an extracellular matrix, randomly distributed postmigratory neurons, multiple branches of thalamic and long corticocortical axons. The SPC generates early spontaneous non-synaptic and synaptic activity and mediates cortical response upon thalamic stimulation. The subplate nexus provides large-scale interareal connectivity possibly underlying fMR resting-state activity, before corticocortical pathways are established. In late fetal phase, when synapses appear within the CP, transient the SPC coexists with permanent circuitry. The histogenetic role of the SPC is to provide interactive milieu and capacity for guidance, sorting, "waiting" and target selection of thalamocortical and corticocortical pathways. The new evolutionary role of the SPC and its remnant white matter neurons is linked to the increasing number of associative pathways in the human neocortex. These roles attributed to the SPC are regulated using a spatiotemporal gene expression during critical periods, when pathogenic factors may disturb vulnerable circuitry of the SPC, causing neurodevelopmental cognitive circuitry disorders.

The Genetics of Epilepsy

Epilepsy encompasses a group of heterogeneous brain diseases that affect more than 50 million people worldwide. Epilepsy may have discernible structural, infectious, metabolic, and immune etiologies; however, in most people with epilepsy, no obvious cause is identifiable. Based initially on family studies and later on advances in gene sequencing technologies and computational approaches, as well as the establishment of large collaborative initiatives, we now know that genetics plays a much greater role in epilepsy than was previously appreciated. Here, we review the progress in the field of epilepsy genetics and highlight molecular discoveries in the most important epilepsy groups, including those that have been long considered to have a nongenetic cause. We discuss where the field of epilepsy genetics is moving as it enters a new era in which the genetic architecture of common epilepsies is starting to be unraveled.

Polygenic Risk Scores Shed Light on the Relationship between Schizophrenia and Cognitive Functioning: Review and Meta-Analysis

Schizophrenia is a multifactorial disease associated with widespread cognitive impairment. Although cognitive deficits are one of the factors most strongly associated with functional impairment in schizophrenia (SZ), current treatment strategies hardly tackle these impairments. To develop more efficient treatment strategies in patients, a better understanding of their pathogenesis is needed. Recent progress in genetics, driven by large genome-wide association studies (GWAS) and the use of polygenic risk scores (PRS), has provided new insights about the genetic architecture of complex human traits, including cognition and SZ. Here, we review the recent findings examining the genetic links between SZ and cognitive functions in population-based samples as well as in participants with SZ. The performed meta-analysis showed a negative correlation between the polygenetic risk score of schizophrenia and global cognition (p < 0.001) when the samples rely on general and healthy participants, while no significant correlation was detected when the three studies devoted to schizophrenia patients were meta-analysed (p > 0.05). Our review and meta-analysis therefore argues against universal pleiotropy for schizophrenia alleles and cognition, since cognition in SZ patients would be underpinned by the same genetic factors than in the general population, and substantially independent of common variant liability to the disorder.

Genome-wide analysis identifies molecular systems and 149 genetic loci associated with income

Socioeconomic position (SEP) is a multi-dimensional construct reflecting (and influencing) multiple socio-cultural, physical, and environmental factors. In a sample of 286,301 participants from UK Biobank, we identify 30 (29 previously unreported) independent-loci associated with income. Using a method to meta-analyze data from genetically-correlated traits, we identify an additional 120 income-associated loci. These loci show clear evidence of functionality, with transcriptional differences identified across multiple cortical tissues, and links to GABAergic and serotonergic neurotransmission. By combining our genome wide association study on income with data from eQTL studies and chromatin interactions, 24 genes are prioritized for follow up, 18 of which were previously associated with intelligence. We identify intelligence as one of the likely causal, partly-heritable phenotypes that might bridge the gap between molecular genetic inheritance and phenotypic consequence in terms of income differences. These results indicate that, in modern era Great Britain, genetic effects contribute towards some of the observed socioeconomic inequalities.

Genome-wide analysis identifies molecular systems and 149 genetic loci associated with income

Socioeconomic position (SEP) is a multi-dimensional construct reflecting (and influencing) multiple socio-cultural, physical, and environmental factors. In a sample of 286,301 participants from UK Biobank, we identify 30 (29 previously unreported) independent-loci associated with income. Using a method to meta-analyze data from genetically-correlated traits, we identify an additional 120 income-associated loci. These loci show clear evidence of functionality, with transcriptional differences identified across multiple cortical tissues, and links to GABAergic and serotonergic neurotransmission. By combining our genome wide association study on income with data from eQTL studies and chromatin interactions, 24 genes are prioritized for follow up, 18 of which were previously associated with intelligence. We identify intelligence as one of the likely causal, partly-heritable phenotypes that might bridge the gap between molecular genetic inheritance and phenotypic consequence in terms of income differences. These results indicate that, in modern era Great Britain, genetic effects contribute towards some of the observed socioeconomic inequalities.

The genetic relationship between educational attainment and cognitive performance in major psychiatric disorders

Cognitive deficits are a core feature of psychiatric disorders like schizophrenia and bipolar disorder. Evidence supports a genome-wide polygenic score (GPS) for educational attainment (GPSEDU) can be used to explain variability in cognitive performance. We aimed to identify different cognitive domains associated with GPSEDU in a transdiagnostic clinical cohort of chronic psychiatric patients with known cognitive deficits. Bipolar and schizophrenia patients from the PsyCourse cohort (N = 730; 43% female) were used. Likewise, we tested whether GPSs for schizophrenia (GPSSZ) and bipolar disorder (GPSBD) were associated with cognitive outcomes. GPSEDU explained 1.5% of variance in the backward verbal digit span, 1.9% in the number of correctly recalled words of the Verbal Learning and Memory Test, and 1.1% in crystallized intelligence. These effects were robust to the influences of treatment and diagnosis. No significant associations between GPSSZ or GPSBD with cognitive outcomes were found. Furthermore, these risk scores did not confound the effect of GPSEDU on cognitive outcomes. GPSEDU explains a small fraction of cognitive performance in adults with psychiatric disorders, specifically for domains related to linguistic learning and working memory. Investigating such a proxy-phenotype longitudinally, could give intriguing insight into the disease course, highlighting at what time genes play a more influential role on cognitive performance. Better understanding the origin of these deficits might help identify those patients at risk for lower levels of functioning and poor social outcomes. Polygenic estimates may in the future be part of predictive models for more personalized interventions.