Genome-wide association study reveals two new risk loci for bipolar disorder

A bipolar disorder-associated missense variant alters adenylyl cyclase 2 activity and promotes mania-like behavior

The single nucleotide polymorphism rs13166360, causing a substitution of valine (Val) 147 to leucine (Leu) in the adenylyl cyclase 2 (ADCY2), has previously been associated with bipolar disorder (BD). Here we show that the disease-associated ADCY2 missense mutation diminishes the enzyme´s capacity to generate the second messenger 3',5'-cylic adenosine monophosphate (cAMP) by altering its subcellular localization. We established mice specifically carrying the Val to Leu substitution using CRISPR/Cas9-based gene editing. Mice homozygous for the Leu variant display symptoms of a mania-like state accompanied by cognitive impairments. Mutant animals show additional characteristic signs of rodent mania models, i.e., they are hypersensitive to amphetamine, the observed mania-like behaviors are responsive to lithium treatment and the Val to Leu substitution results in a shifted excitatory/inhibitory synaptic balance towards more excitation. Exposure to chronic social defeat stress switches homozygous Leu variant carriers from a mania- to a depressive-like state, a transition which is reminiscent of the alternations characterizing the symptomatology in BD patients. Single-cell RNA-seq (scRNA-seq) revealed widespread Adcy2 mRNA expression in numerous hippocampal cell types. Differentially expressed genes particularly identified from glutamatergic CA1 neurons point towards ADCY2 variant-dependent alterations in multiple biological processes including cAMP-related signaling pathways. These results validate ADCY2 as a BD risk gene, provide insights into underlying disease mechanisms, and potentially open novel avenues for therapeutic intervention strategies.

Novel genetic loci of inhibitory control in ADHD and healthy children and genetic correlations with ADHD

Cumulative evidence has showed the deficits of inhibitory control in patients with attention deficit hyperactivity disorder (ADHD), which is considered as an endophenotype of ADHD. Genetic study of inhibitory control could advance gene discovery and further facilitate the understanding of ADHD genetic basis, but the studies were limited in both the general population and ADHD patients. To reveal genetic risk variants of inhibitory control and its potential genetic relationship with ADHD, we conducted genome-wide association studies (GWAS) on inhibitory control using three datasets, which included 783 and 957 ADHD patients and 1350 healthy children. Subsequently, we employed polygenic risk scores (PRS) to explore the association of inhibitory control with ADHD and related psychiatric disorders. Firstly, we identified three significant loci for inhibitory control in the healthy dataset, two loci in the case dataset, and one locus in the meta-analysis of three datasets. Besides, we found more risk genes and variants by applying transcriptome-wide association study (TWAS) and conditional FDR method. Then, we constructed a network by connecting the genes identified in our study, leading to the identification of several vital genes. Lastly, we identified a potential relationship between inhibitory control and ADHD and autism by PRS analysis and found the direct and mediated contribution of the identified genetic loci on ADHD symptoms by mediation analysis. In conclusion, we revealed some genetic risk variants associated with inhibitory control and elucidated the benefit of inhibitory control as an endophenotype, providing valuable insights into the mechanisms underlying ADHD.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Fine-mapping genomic loci refines bipolar disorder risk genes

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Abnormal patterns of sleep and waking behaviors are accompanied by neocortical oscillation disturbances in an Ank3 mouse model of epilepsy-bipolar disorder comorbidity

ANK3 is a leading bipolar disorder (BD) candidate gene in humans and provides a unique opportunity for studying epilepsy-BD comorbidity. Previous studies showed that deletion of Ank3-1b, a BD-associated variant of Ank3 in mice leads to increased firing threshold and diminished action potential dynamic range of parvalbumin (PV) interneurons and absence epilepsy, thus providing a biological mechanism linking epilepsy and BD. To explore the behavioral overlap of these disorders, we characterized behavioral patterns of Ank3-1b KO mice during overnight home-cage activity and examined network activity during these behaviors using paired video and EEG recordings. Since PV interneurons contribute to the generation of high-frequency gamma oscillations, we anticipated changes in the power of neocortical EEG signals in the gamma frequency range (> 25 Hz) during behavioral states related to human BD symptoms, including abnormal sleep, hyperactivity, and repetitive behaviors. Ank3-1b KO mice exhibited an overall increase in slow gamma (~25-45 Hz) power compared to controls, and slow gamma power correlated with seizure phenotype severity across behaviors. During sleep, increased slow gamma power correlated with decreased time spent in the rapid eye movement (REM) stage of sleep. Seizures were more common during REM sleep compared to non-REM (NREM) sleep. We also found that Ank3-1b KO mice were hyperactive and exhibited a repetitive behavior phenotype that co-occurred with increased slow gamma power. Our results identify a novel EEG biomarker associating Ank3 genetic variation with BD and epilepsy and suggest modulation of gamma oscillations as a potential therapeutic target.

Mapping microRNA expression quantitative trait loci in the prenatal human brain implicates miR-1908-5p expression in bipolar disorder and other brain-related traits

MicroRNA (miRNA) are small non-coding RNA involved in post-transcriptional gene regulation. Given their known involvement in early neurodevelopment processes, we here sought to identify common genetic variants associated with altered miRNA expression in the prenatal human brain. We performed small RNA sequencing on brain tissue from 112 genome-wide genotyped fetuses from the second trimester of gestation, identifying high-confidence (false discovery rate < 0.05) expression quantitative trait loci for 30 mature miRNA. Integrating our findings with genome-wide association study data for brain-related disorders, we implicate increased prenatal expression of miR-1908-5p as a risk mechanism for bipolar disorder and find that predicted mRNA targets of miR-1908-5p that are expressed in the fetal brain are enriched for common variant genetic association with the condition. Extending these analyses to other brain-related traits, we find that common genetic variation associated with increased miR-1908-5p expression in fetal brain is additionally associated with depressive symptoms, irritability, increased right cerebellum exterior volume and increased sleep duration in the general population. Our findings provide support to the view that altered miRNA expression can influence susceptibility to neuropsychiatric illness and suggest an early neurodevelopmental risk mechanism for bipolar disorder.

Emerging trends in gene and bipolar disorder research: a bibliometric analysis and network visualisation

This study aims to use a bibliometric technique to evaluate the scientific output of gene and bipolar disorder research. The search query related to gene and bipolar disorder from the Scopus database identified 1848 documents from 1951 to 2020. The growth in the publications increased since early 1990, peaked in 2011, and started to decline thereafter. High occurrence in author keywords suggests that some research topics, such as "polymorphism", "linkage" and "association study" have waned over time, whereas others, such as "DNA methylation," "circadian rhythm," "" and "meta-analysis," are now the emerging trends in gene and bipolar disorder research. The USA was the country with the highest production followed by the UK, Canada, Italy and Germany. The leading institutions were Cardiff University in the UK, the National Institute of Mental Health (NIMH) in the USA, King's College London in the UK and the University of California, San Diego in the USA. The leading journals publishing gene and bipolar literature were the American Journal of Medical Genetics Neuropsychiatric Genetics, Molecular Psychiatry and Psychiatric Genetics. The top authors in the number of publications were Craddock N, Serretti A and Rietschel M. According to the co-authorship network analysis of authors, the majority of the authors in the same clusters were closely linked together and originated from the same or neighbouring country. The findings of this study may be useful in identifying emerging topics for future research and promoting research collaboration in the field of genetic studies related to bipolar disorder.

Comparison of oxidative stress parameters, thiol-disulfide homeostasis, and pro-inflammatory cytokines levels in patients with bipolar disorder and their first-degree relatives

BACKGROUND

In this study, we aimed to compare the oxidative stress parameters, thiol-disulfide homeostasis, and plasma pro-inflammatory cytokines levels of patients with bipolar disorder (BD), BD patients' first-degree relatives (FDRs), and the healthy controls (HCs).

METHODS

Thirty-five patients with BD, 35 FDRs of BD, and 35 healthy controls (HCs) were included. The individuals' ages varied from 28 to 58, and the groups were well-matched in terms of age and gender. The total thiol (TT), native thiol (NT), disulfide (DIS), total oxidant status (TOS), total antioxidant status (TAS), IL-1β, IL-6, and TNF-α concentrations were measured from serum samples. The oxidative stress index (OSI) was calculated using mathematical formulas.

RESULTS

TOS was significantly higher in both patients and FDRs than HCs (p < 0.01 for all pairwise comparisons). OSI, DIS, oxidized thiol, and the ratio of thiol oxidation-reduction levels were significantly higher in both patients with BD and FDRs than HCs (p < 0.01 for all pairwise comparisons). TAS, TT, NT, and reduced thiol levels were significantly lower in both patients with BD and FDRs than HCs (p < 0.01 for all pairwise comparisons). IL-1β, IL-6, and TNF-α were significantly higher in both patients and FDRs than HCs (p < 0.01 for all pairwise comparisons).

LIMITATIONS

Small sample size.

CONCLUSIONS

Early diagnosis is important for treating of bipolar disorder. TT, NT, DIS, TOS, TAS, OSI, IL1-β, IL-6, and TNF-α can be used as potential biomarkers in the early diagnosis and intervention of BD. Furthermore, oxidative/antioxidative markers and plasma pro-inflammatory cytokine parameters may guide the determination of the disease's activity and response to treatment.

Palmitoylation controls the stability of 190 kDa ankyrin-G in dendritic spines and is regulated by ZDHHC8 and lithium

Introduction

AnkG, encoded by the ANK3 gene, is a multifunctional scaffold protein with complex isoform expression: the 480 and 270 kDa isoforms have roles at the axon initial segment and node of Ranvier, whereas the 190 kDa isoform (AnkG-190) has an emerging role in the dendritic shaft and spine heads. All isoforms of AnkG undergo palmitoylation, a post-translational modification regulating protein attachment to lipid membranes. However, palmitoylation of AnkG-190 has not been investigated in dendritic spines. The ANK3 gene and altered expression of AnkG proteins are associated with a variety of neuropsychiatric and neurodevelopmental disorders including bipolar disorder and are implicated in the lithium response, a commonly used mood stabilizer for bipolar disorder patients, although the precise mechanisms involved are unknown.

Result

Here, we showed that Cys70 palmitoylation stabilizes the localization of AnkG-190 in spine heads and at dendritic plasma membrane nanodomains. Mutation of Cys70 impairs AnkG-190 function in dendritic spines and alters PSD-95 scaffolding. Interestingly, we find that lithium reduces AnkG-190 palmitoylation thereby increasing its mobility in dendritic spines. Finally, we demonstrate that the palmitoyl acyl transferase ZDHHC8, but not ZDHHC5, increases AnkG-190 stability in spine heads and is inhibited by lithium.

Discussion

Together, our data reveal that palmitoylation is critical for AnkG-190 localization and function and a potential ZDHHC8/AnkG-190 mechanism linking AnkG-190 mobility to the neuronal effects of lithium.

Genetics of bipolar disorder: insights into its complex architecture and biology from common and rare variants

Bipolar disorder (BD) is a common mental disorder characterized by recurrent mood episodes, which causes major socioeconomic burdens globally. Though its disease pathogenesis is largely unknown, the high heritability of BD indicates strong contributions from genetic factors. In this review, we summarize the recent achievements in the genetics of BD, particularly those from genome-wide association study (GWAS) of common variants and next-generation sequencing analysis of rare variants. These include the identification of dozens of robust disease-associated loci, deepening of our understanding of the biology of BD, objective description of correlations with other psychiatric disorders and behavioral traits, formulation of methods for predicting disease risk and drug response, and the discovery of a single gene associated with bipolar disorder and schizophrenia spectrum with a large effect size. On the other hand, the findings to date have not yet made a clear contribution to the improvement of clinical psychiatry of BD. We overview the remaining challenges as well as possible paths to resolve them, referring to studies of other major neuropsychiatric disorders.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.

Gene Therapy for Neuropsychiatric Disorders: Potential Targets and Tools

Neuropsychiatric disorders that affect the central nervous system cause considerable pressures on the health care system and have a substantial economic burden on modern societies. The present treatments based on available drugs are mostly ineffective and often costly. The molecular process of neuropsychiatric disorders is closely connected to modifying the genetic structures inherited or caused by damage, toxic chemicals, and some current diseases. Gene therapy is presently an experimental concept for neurological disorders. Clinical applications endeavor to alleviate the symptoms, reduce disease progression, and repair defective genes. Implementing gene therapy in inherited and acquired neurological illnesses entails the integration of several scientific disciplines, including virology, neurology, neurosurgery, molecular genetics, and immunology. Genetic manipulation has the power to minimize or cure illness by inducing genetic alterations at endogenous loci. Gene therapy that involves treating the disease by deleting, silencing, or editing defective genes and delivering genetic material to produce therapeutic molecules has excellent potential as a novel approach for treating neuropsychiatric disorders. With the recent advances in gene selection and vector design quality in targeted treatments, gene therapy could be an effective approach. This review article will investigate and report the newest and the most critical molecules and factors in neuropsychiatric disorder gene therapy. Different genome editing techniques available will be evaluated, and the review will highlight preclinical research of genome editing for neuropsychiatric disorders while also evaluating current limitations and potential strategies to overcome genome editing advancements.

Early cognitive comorbidities before disease onset: A common symptom towards prevention of related brain diseases?

Brain diseases are very heterogeneous; however they also display multiple common risk factors and comorbidities. With a paucity of disease-modifying therapies, prevention became a health priority. Towards prevention, one strategy is to focus on similar symptoms of brain diseases occurring before disease onset. Cognitive deficits are a promising candidate as they occur across brain diseases before disease onset. Based on recent research, this review highlights the similarity of brain diseases and discusses how early cognitive deficits can be exploited to tackle disease prevention. After briefly introducing common risk factors, I review common comorbidities across brain diseases, with a focus on cognitive deficits before disease onset, reporting both experimental and clinical findings. Next, I describe network abnormalities associated with early cognitive deficits and discuss how these abnormalities can be targeted to prevent disease onset. A scenario on brain disease etiology with the idea that early cognitive deficits may constitute a common symptom of brain diseases is proposed.

Mapping brain gene coexpression in daytime transcriptomes unveils diurnal molecular networks and deciphers perturbation gene signatures

Brain tissue transcriptomes may be organized into gene coexpression networks, but their underlying biological drivers remain incompletely understood. Here, we undertook a large-scale transcriptomic study using 508 wild-type mouse striatal tissue samples dissected exclusively in the afternoons to define 38 highly reproducible gene coexpression modules. We found that 13 and 11 modules are enriched in cell-type and molecular complex markers, respectively. Importantly, 18 modules are highly enriched in daily rhythmically expressed genes that peak or trough with distinct temporal kinetics, revealing the underlying biology of striatal diurnal gene networks. Moreover, the diurnal coexpression networks are a dominant feature of daytime transcriptomes in the mouse cortex. We next employed the striatal coexpression modules to decipher the striatal transcriptomic signatures from Huntington's disease models and heterozygous null mice for 52 genes, uncovering novel functions for Prkcq and Kdm4b in oligodendrocyte differentiation and bipolar disorder-associated Trank1 in regulating anxiety-like behaviors and nocturnal locomotion.

Genetic and functional analyses implicate microRNA 499A in bipolar disorder development

Bipolar disorder (BD) is a complex mood disorder with a strong genetic component. Recent studies suggest that microRNAs contribute to psychiatric disorder development. In BD, specific candidate microRNAs have been implicated, in particular miR-137, miR-499a, miR-708, miR-1908 and miR-2113. The aim of the present study was to determine the contribution of these five microRNAs to BD development. For this purpose, we performed: (i) gene-based tests of the five microRNA coding genes, using data from a large genome-wide association study of BD; (ii) gene-set analyses of predicted, brain-expressed target genes of the five microRNAs; (iii) resequencing of the five microRNA coding genes in 960 BD patients and 960 controls and (iv) in silico and functional studies for selected variants. Gene-based tests revealed a significant association with BD for MIR499A, MIR708, MIR1908 and MIR2113. Gene-set analyses revealed a significant enrichment of BD associations in the brain-expressed target genes of miR-137 and miR-499a-5p. Resequencing identified 32 distinct rare variants (minor allele frequency < 1%), all of which showed a non-significant numerical overrepresentation in BD patients compared to controls (p = 0.214). Seven rare variants were identified in the predicted stem-loop sequences of MIR499A and MIR2113. These included rs142927919 in MIR2113 (p_nom = 0.331) and rs140486571 in MIR499A (p_nom = 0.297). In silico analyses predicted that rs140486571 might alter the miR-499a secondary structure. Functional analyses showed that rs140486571 significantly affects miR-499a processing and expression. Our results suggest that MIR499A dysregulation might contribute to BD development. Further research is warranted to elucidate the contribution of the MIR499A regulated network to BD susceptibility.

Correlations between peripheral levels of inflammatory mediators and frontolimbic structures in bipolar disorder: an exploratory analysis

BACKGROUND

Altered peripheral immune/inflammatory system and brain volumetric changes have been implicated in the pathophysiology of bipolar disorder (BD). This study aimed to evaluate how peripheral levels of cytokines are related to volumetric brain changes in euthymic patients with BD.

METHODS

Euthymic patients with BD (n = 21) and healthy controls (n = 22) were enrolled in this exploratory study. Blood samples were collected on the same day of clinical assessment and neuroimaging. Cytokines were measured through cytometric bead array method. Neuroimaging data were acquired using a sagittal three-dimensional magnetic resonance imaging T1-weighted fast field echo sequence and was processed using FreeSurfer.

RESULTS

Compared to controls, BD patients had significantly lower volumes in the cingulate, medial-orbitofrontal (MOF) and parahippocampal regions. We found a negative correlation between right MOF volume and interferon-gamma levels (β = -0.431, P = .049) and a positive correlation between interleukin-10 levels and left posterior cingulate volume (β = 0.457, P = .048).

CONCLUSION

Our results support the involvement of inflammatory pathways in structural brain changes in BD.

Genetic Polymorphisms in miR-137 and Its Target Genes, TCF4 and CACNA1C, Contribute to the Risk of Bipolar Disorder: A Preliminary Case-Control Study and Bioinformatics Analysis

Accumulating evidence has suggested that miR-137 and its target genes, CACNA1C, and TCF4, are amongst the most robustly implicated genes in psychiatric disorders. This preliminary study is aimed at investigating the effects of genetic variations in miR-137 (rs1625579A/C), TCF4 (rs1261084C/T), and CACNA1C (rs10774053A/G and rs10466907G/T) on BD susceptibility. We recruited 252 BD patients and 213 healthy subjects as the control group. Genotyping was performed using PCR-RFLP and ARMS-PCR methods. Enhanced risk of BD was found under the codominant homozygous, dominant, and allelic models of TCF4 rs1261084C/T, codominant homozygous and allelic models of CACNA1C rs10466907G/T polymorphisms, as well as codominant homozygous, dominant, recessive, and allelic models of the CACNA1C rs10774053A/G. Moreover, both TT/AG/GT/AA and TT/GG/GT/AC genotype combinations strongly increased the risk of BD in the participants. The bioinformatics analyses revealed that rs1261084C/T and rs10466907G/T created and disrupted binding sites of some miRNAs in the 3′-untranslated region of TCF4 and CACNA1C genes. In contrast, the rs10774053A/G created a new binding site for a major splicing factor and might have an effective role in the function of the CACNA1C protein. We have found that all the studied SNPs are positively associated with BD susceptibility. Replicated studies on different ethnicities are required to confirm these findings.

A Systematic Review of Genetic Polymorphisms Associated with Bipolar Disorder Comorbid to Substance Abuse

It is currently unknown which genetic polymorphisms are involved in substance use disorder (SUD) comorbid with bipolar disorder (BD). The research on polymorphisms in BD comorbid with SUD (BD + SUD) is summarized in this systematic review. We looked for case-control studies that genetically compared adults and adolescents with BD and SUD, healthy controls, and BD without SUD. PRISMA was used to create our protocol, which is PROSPERO-registered (identification: CRD4221270818). The following bibliographic databases were searched indefinitely until December 2021 to identify potentially relevant articles: PubMed, PsycINFO, Scopus, and Web of Science. This systematic review, after the qualitative analysis of the study selection, included 17 eligible articles. In the selected studies, 66 polymorphisms in 29 genes were investigated. The present work delivers a group of potentially valuable genetic polymorphisms associated with BD + SUD: rs11600996 (ARNTL), rs228642/rs228682/rs2640909 (PER3), PONQ192R (PON1), rs945032 (BDKRB2), rs1131339 (NR4A3), and rs6971 (TSPO). It is important to note that none of those findings have been confirmed by two or more studies; thus, we believe that all the polymorphisms identified in this review require additional evidence to be confirmed.

The molecular pathophysiology of mood disorders: From the analysis of single molecular layers to multi-omic integration

Next-generation sequencing now enables the rapid and affordable production of reliable biological data at multiple molecular levels, collectively referred to as "omics". To maximize the potential for discovery, computational biologists have created and adapted integrative multi-omic analytical methods. When applied to diseases with traceable pathophysiology such as cancer, these new algorithms and statistical approaches have enabled the discovery of clinically relevant molecular mechanisms and biomarkers. In contrast, these methods have been much less applied to the field of molecular psychiatry, although diagnostic and prognostic biomarkers are similarly needed. In the present review, we first briefly summarize main findings from two decades of studies that investigated single molecular processes in relation to mood disorders. Then, we conduct a systematic review of multi-omic strategies that have been proposed and used more recently. We also list databases and types of data available to researchers for future work. Finally, we present the newest methodologies that have been employed for multi-omics integration in other medical fields, and discuss their potential for molecular psychiatry studies.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

Background

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

Methods

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding–disrupting SNPs (hereafter referred to as “functional SNPs”) through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

Results

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding–disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

Conclusions

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01057-3.

Genetic Regulation of Vertebrate Forebrain Development by Homeobox Genes

Forebrain development in vertebrates is regulated by transcription factors encoded by homeobox, bHLH and forkhead gene families throughout the progressive and overlapping stages of neural induction and patterning, regional specification and generation of neurons and glia from central nervous system (CNS) progenitor cells. Moreover, cell fate decisions, differentiation and migration of these committed CNS progenitors are controlled by the gene regulatory networks that are regulated by various homeodomain-containing transcription factors, including but not limited to those of the Pax (paired), Nkx, Otx (orthodenticle), Gsx/Gsh (genetic screened), and Dlx (distal-less) homeobox gene families. This comprehensive review outlines the integral role of key homeobox transcription factors and their target genes on forebrain development, focused primarily on the telencephalon. Furthermore, links of these transcription factors to human diseases, such as neurodevelopmental disorders and brain tumors are provided.

Zika virus-induced TNF-α signaling dysregulates expression of neurologic genes associated with psychiatric disorders

BACKGROUND

Zika virus (ZIKV) is an emerging flavivirus of global concern. ZIKV infection of the central nervous system has been linked to a variety of clinical syndromes, including microcephaly in fetuses and rare but serious neurologic disease in adults. However, the potential for ZIKV to influence brain physiology and host behavior following apparently mild or subclinical infection is less well understood. Furthermore, though deficits in cognitive function are well-documented after recovery from neuroinvasive viral infection, the potential impact of ZIKV on other host behavioral domains has not been thoroughly explored.

METHODS

We used transcriptomic profiling, including unbiased gene ontology enrichment analysis, to assess the impact of ZIKV infection on gene expression in primary cortical neuron cultures. These studies were extended with molecular biological analysis of gene expression and inflammatory cytokine signaling. In vitro observations were further confirmed using established in vivo models of ZIKV infection in immunocompetent hosts.

RESULTS

Transcriptomic profiling of primary neuron cultures following ZIKV infection revealed altered expression of key genes associated with major psychiatric disorders, such as bipolar disorder and schizophrenia. Gene ontology enrichment analysis also revealed significant changes in gene expression associated with fundamental neurobiological processes, including neuronal development, neurotransmission, and others. These alterations to neurologic gene expression were also observed in the brain in vivo using several immunocompetent mouse models of ZIKV infection. Mechanistic studies identified TNF-α signaling via TNFR1 as a major regulatory mechanism controlling ZIKV-induced changes to neurologic gene expression.

CONCLUSIONS

Our studies reveal that cell-intrinsic innate immune responses to ZIKV infection profoundly shape neuronal transcriptional profiles, highlighting the need to further explore associations between ZIKV infection and disordered host behavioral states.

Genomic and neuroimaging approaches to bipolar disorder

BACKGROUND

To date, besides genome-wide association studies, a variety of other genetic analyses (e.g. polygenic risk scores, whole-exome sequencing and whole-genome sequencing) have been conducted, and a large amount of data has been gathered for investigating the involvement of common, rare and very rare types of DNA sequence variants in bipolar disorder. Also, non-invasive neuroimaging methods can be used to quantify changes in brain structure and function in patients with bipolar disorder.

AIMS

To provide a comprehensive assessment of genetic findings associated with bipolar disorder, based on the evaluation of different genomic approaches and neuroimaging studies.

METHOD

We conducted a PubMed search of all relevant literatures from the beginning to the present, by querying related search strings.

RESULTS

ANK3, CACNA1C, SYNE1, ODZ4 and TRANK1 are five genes that have been replicated as key gene candidates in bipolar disorder pathophysiology, through the investigated studies. The percentage of phenotypic variance explained by the identified variants is small (approximately 4.7%). Bipolar disorder polygenic risk scores are associated with other psychiatric phenotypes. The ENIGMA-BD studies show a replicable pattern of lower cortical thickness, altered white matter integrity and smaller subcortical volumes in bipolar disorder.

CONCLUSIONS

The low amount of explained phenotypic variance highlights the need for further large-scale investigations, especially among non-European populations, to achieve a more complete understanding of the genetic architecture of bipolar disorder and the missing heritability. Combining neuroimaging data with genetic data in large-scale studies might help researchers acquire a better knowledge of the engaged brain regions in bipolar disorder.

Teneurin4 dimer structures reveal a calcium‐stabilized compact conformation supporting homomeric trans‐interactions

Establishment of correct synaptic connections is a crucial step during neural circuitry formation. The Teneurin family of neuronal transmembrane proteins promotes cell–cell adhesion via homophilic and heterophilic interactions, and is required for synaptic partner matching in the visual and hippocampal systems in vertebrates. It remains unclear how individual Teneurins form macromolecular cis‐ and trans‐synaptic protein complexes. Here, we present a 2.7 Å cryo‐EM structure of the dimeric ectodomain of human Teneurin4. The structure reveals a compact conformation of the dimer, stabilized by interactions mediated by the C‐rich, YD‐shell, and ABD domains. A 1.5 Å crystal structure of the C‐rich domain shows three conserved calcium binding sites, and thermal unfolding assays and SAXS‐based rigid‐body modeling demonstrate that the compactness and stability of Teneurin4 dimers are calcium‐dependent. Teneurin4 dimers form a more extended conformation in conditions that lack calcium. Cellular assays reveal that the compact cis‐dimer is compatible with homomeric trans‐interactions. Together, these findings support a role for teneurins as a scaffold for macromolecular complex assembly and the establishment of cis‐ and trans‐synaptic interactions to construct functional neuronal circuits.

A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics

BACKGROUND

Recent genome-wide association studies (GWASs) have revealed the polygenic nature of psychiatric disorders and discovered a few of single-nucleotide polymorphisms (SNPs) associated with multiple psychiatric disorders. However, the extent and pattern of pleiotropy among distinct psychiatric disorders remain not completely clear.

METHODS

We analyzed 14 psychiatric disorders using summary statistics available from the largest GWASs by far. We first applied the cross-trait linkage disequilibrium score regression (LDSC) to estimate genetic correlation between disorders. Then, we performed a gene-based pleiotropy analysis by first aggregating a set of SNP-level associations into a single gene-level association signal using MAGMA. From a methodological perspective, we viewed the identification of pleiotropic associations across the entire genome as a high-dimensional problem of composite null hypothesis testing and utilized a novel method called PLACO for pleiotropy mapping. We ultimately implemented functional analysis for identified pleiotropic genes and used Mendelian randomization for detecting causal association between these disorders.

RESULTS

We confirmed extensive genetic correlation among psychiatric disorders, based on which these disorders can be grouped into three diverse categories. We detected a large number of pleiotropic genes including 5884 associations and 2424 unique genes and found that differentially expressed pleiotropic genes were significantly enriched in pancreas, liver, heart, and brain, and that the biological process of these genes was remarkably enriched in regulating neurodevelopment, neurogenesis, and neuron differentiation, offering substantial evidence supporting the validity of identified pleiotropic loci. We further demonstrated that among all the identified pleiotropic genes there were 342 unique ones linked with 6353 drugs with drug-gene interaction which can be classified into distinct types including inhibitor, agonist, blocker, antagonist, and modulator. We also revealed causal associations among psychiatric disorders, indicating that genetic overlap and causality commonly drove the observed co-existence of these disorders.

CONCLUSIONS

Our study is among the first large-scale effort to characterize gene-level pleiotropy among a greatly expanded set of psychiatric disorders and provides important insight into shared genetic etiology underlying these disorders. The findings would inform psychiatric nosology, identify potential neurobiological mechanisms predisposing to specific clinical presentations, and pave the way to effective drug targets for clinical treatment.

Astrocytes in Bipolar Disorder

Bipolar disorder (BD) is a complex group of neuropsychiatric disorders, typically comprising both manic and depressive episodes. The underlying neuropathology of BD is not established, but a consistent feature is progressive thinning of cortical grey matter (GM) and white matter (WM) in specific pathways, due to loss of subpopulations of neurons and astrocytes, with accompanying disturbance of connectivity. Dysregulation of astrocyte homeostatic functions are implicated in BD, notably regulation of glutamate, calcium signalling, circadian rhythms and metabolism. Furthermore, the beneficial therapeutic effects of the frontline treatments for BD are due at least in part to their positive actions on astrocytes, notably lithium, valproic acid (VPA) and carbamazepine (CBZ), as well as antidepressants and antipsychotics that are used in the management of this disorder. Treatments for BD are ineffective in a large proportion of cases, and astrocytes represent new therapeutic targets that can also serve as biomarkers of illness progression and treatment responsiveness in BD.

Effects of polygenic risk for suicide attempt and risky behavior on brain structure in young people with familial risk of bipolar disorder

Bipolar disorder (BD) is associated with a 20-30-fold increased suicide risk compared to the general population. First-degree relatives of BD patients show inflated rates of psychopathology including suicidal behaviors. As reliable biomarkers of suicide attempts (SA) are lacking, we examined associations between suicide-related polygenic risk scores (PRSs)-a quantitative index of genomic risk-and variability in brain structures implicated in SA. Participants (n = 206; aged 12-30 years) were unrelated individuals of European ancestry and comprised three groups: 41 BD cases, 96 BD relatives ("high risk"), and 69 controls. Genotyping employed PsychArray, followed by imputation. Three PRSs were computed using genome-wide association data for SA in BD (SA-in-BD), SA in major depressive disorder (SA-in-MDD) (Mullins et al., 2019, The American Journal of Psychiatry, 176(8), 651-660), and risky behavior (Karlsson Linnér et al., 2019, Nature Genetics, 51(2), 245-257). Structural magnetic resonance imaging processing employed FreeSurfer v5.3.0. General linear models were constructed using 32 regions-of-interest identified from suicide neuroimaging literature, with false-discovery-rate correction. SA-in-MDD and SA-in-BD PRSs negatively predicted parahippocampal thickness, with the latter association modified by group membership. SA-in-BD and Risky Behavior PRSs inversely predicted rostral and caudal anterior cingulate structure, respectively, with the latter effect driven by the "high risk" group. SA-in-MDD and SA-in-BD PRSs positively predicted cuneus structure, irrespective of group. This study demonstrated associations between PRSs for suicide-related phenotypes and structural variability in brain regions implicated in SA. Future exploration of extended PRSs, in conjunction with a range of biological, phenotypic, environmental, and experiential data in high risk populations, may inform predictive models for suicidal behaviors.

Association study of BDNF Val66Met gene polymorphism with bipolar disorder and lithium treatment response in Indian population

BACKGROUND

The association of the Val66Met (rs6265) polymorphism in the brain-derived neurotrophic factor (BDNF) gene with bipolar disorder (BD) and response to lithium treatment has been suggested, though inconsistently. The considerable diversity of allele frequency across different populations contributes to this. There is no data from South Asia till date. Hence, we examined the association of this polymorphism in BD cases from India, and its association with lithium treatment response.

METHODS

BD patients (N = 301) were recruited from the clinical services of National Institute of Mental Health and Neurosciences (NIMHANS), India. Lithium treatment response for 190 BD subjects was assessed using Alda scale by NIMH life charts. Patients with total score ⩾7 were defined as lithium responders (N = 115) and patients with score <7 were defined as lithium non-responders (N = 75). Healthy controls (N = 484) with no lifetime history of neuropsychiatric illness or a family history of mental illness were recruited as control set. Genotyping was performed by TaqMan genotyping assay.

RESULTS

Genotype and allele frequency of BDNF Val66Met SNP was significantly different (χ2 = 7.78, p = 0.02) in cases compared to controls, and the Val(G) allele was more frequent (χ2 = 7.08, p = 0.008) in BD patients. However, no significant difference is noted in genotype or allele frequencies of this polymorphism between the lithium responders and non-responders.

CONCLUSIONS

The Val(G) allele of BDNF Val66Met polymorphism is associated with risk of BD in this sample, but it is not related to response to lithium.

Physiological Roles of Mammalian Transmembrane Adenylyl Cyclase Isoforms

Adenylyl cyclases (ACs) catalyze the conversion of ATP to the ubiquitous second messenger cAMP. Mammals possess nine isoforms of transmembrane ACs, dubbed AC1-9, that serve as major effector enzymes of G protein-coupled receptors. The transmembrane ACs display varying expression patterns across tissues, giving potential for them having a wide array of physiologic roles. Cells express multiple AC isoforms, implying that ACs have redundant functions. Furthermore, all transmembrane ACs are activated by Gα_s so it was long assumed that all ACs are activated by Gα_s-coupled GPCRs. AC isoforms partition to different microdomains of the plasma membrane and form prearranged signaling complexes with specific GPCRs that contribute to cAMP signaling compartments. This compartmentation allows for a diversity of cellular and physiological responses by enabling unique signaling events to be triggered by different pools of cAMP. Isoform specific pharmacological activators or inhibitors are lacking for most ACs, making knockdown and overexpression the primary tools for examining the physiological roles of a given isoform. Much progress has been made in understanding the physiological effects mediated through individual transmembrane ACs. GPCR-AC-cAMP signaling pathways play significant roles in regulating functions of every cell and tissue, so understanding each AC isoform's role holds potential for uncovering new approaches for treating a vast array of pathophysiological conditions.

Genetic contributions to bipolar disorder: current status and future directions

Bipolar disorder (BD) is a highly heritable mental disorder and is estimated to affect about 50 million people worldwide. Our understanding of the genetic etiology of BD has greatly increased in recent years with advances in technology and methodology as well as the adoption of international consortiums and large population-based biobanks. It is clear that BD is also highly heterogeneous and polygenic and shows substantial genetic overlap with other psychiatric disorders. Genetic studies of BD suggest that the number of associated loci is expected to substantially increase in larger future studies and with it, improved genetic prediction of the disorder. Still, a number of challenges remain to fully characterize the genetic architecture of BD. First among these is the need to incorporate ancestrally-diverse samples to move research away from a Eurocentric bias that has the potential to exacerbate health disparities already seen in BD. Furthermore, incorporation of population biobanks, registry data, and electronic health records will be required to increase the sample size necessary for continued genetic discovery, while increased deep phenotyping is necessary to elucidate subtypes within BD. Lastly, the role of rare variation in BD remains to be determined. Meeting these challenges will enable improved identification of causal variants for the disorder and also allow for equitable future clinical applications of both genetic risk prediction and therapeutic interventions.

Identifying pleiotropic genes for complex phenotypes with summary statistics from a perspective of composite null hypothesis testing

Pleiotropy has important implication on genetic connection among complex phenotypes and facilitates our understanding of disease etiology. Genome-wide association studies provide an unprecedented opportunity to detect pleiotropic associations; however, efficient pleiotropy test methods are still lacking. We here consider pleiotropy identification from a methodological perspective of high-dimensional composite null hypothesis and propose a powerful gene-based method called MAIUP. MAIUP is constructed based on the traditional intersection-union test with two sets of independent P-values as input and follows a novel idea that was originally proposed under the high-dimensional mediation analysis framework. The key improvement of MAIUP is that it takes the composite null nature of pleiotropy test into account by fitting a three-component mixture null distribution, which can ultimately generate well-calibrated P-values for effective control of family-wise error rate and false discover rate. Another attractive advantage of MAIUP is its ability to effectively address the issue of overlapping subjects commonly encountered in association studies. Simulation studies demonstrate that compared with other methods, only MAIUP can maintain correct type I error control and has higher power across a wide range of scenarios. We apply MAIUP to detect shared associated genes among 14 psychiatric disorders with summary statistics and discover many new pleiotropic genes that are otherwise not identified if failing to account for the issue of composite null hypothesis testing. Functional and enrichment analyses offer additional evidence supporting the validity of these identified pleiotropic genes associated with psychiatric disorders. Overall, MAIUP represents an efficient method for pleiotropy identification.

SUPERGNOVA: local genetic correlation analysis reveals heterogeneous etiologic sharing of complex traits

Local genetic correlation quantifies the genetic similarity of complex traits in specific genomic regions. However, accurate estimation of local genetic correlation remains challenging, due to linkage disequilibrium in local genomic regions and sample overlap across studies. We introduce SUPERGNOVA, a statistical framework to estimate local genetic correlations using summary statistics from genome-wide association studies. We demonstrate that SUPERGNOVA outperforms existing methods through simulations and analyses of 30 complex traits. In particular, we show that the positive yet paradoxical genetic correlation between autism spectrum disorder and cognitive performance could be explained by two etiologically distinct genetic signatures with bidirectional local genetic correlations.

Kleine-Levin syndrome is associated with birth difficulties and genetic variants in the TRANK1 gene loci

Kleine-Levin syndrome (KLS) is a rare disorder characterized by severe episodic hypersomnia, with cognitive impairment accompanied by apathy or disinhibition. Pathophysiology is unknown, although imaging studies indicate decreased activity in hypothalamic/thalamic areas during episodes. Familial occurrence is increased, and risk is associated with reports of a difficult birth. We conducted a worldwide case-control genome-wide association study in 673 KLS cases collected over 14 y, and ethnically matched 15,341 control individuals. We found a strong genome-wide significant association (rs71947865, Odds Ratio [OR] = 1.48, P = 8.6 × 10^-9) within the 3'region of TRANK1 gene locus, previously associated with bipolar disorder and schizophrenia. Strikingly, KLS cases with rs71947865 variant had significantly increased reports of a difficult birth. As perinatal outcomes have dramatically improved over the last 40 y, we further stratified our sample by birth years and found that recent cases had a significantly reduced rs71947865 association. While the rs71947865 association did not replicate in the entire follow-up sample of 171 KLS cases, rs71947865 was significantly associated with KLS in the subset follow-up sample of 59 KLS cases who reported birth difficulties (OR = 1.54, P = 0.01). Genetic liability of KLS as explained by polygenic risk scores was increased (pseudo R ² = 0.15; P < 2.0 × 10^-22 at P = 0.5 threshold) in the follow-up sample. Pathway analysis of genetic associations identified enrichment of circadian regulation pathway genes in KLS cases. Our results suggest links between KLS, circadian regulation, and bipolar disorder, and indicate that the TRANK1 polymorphisms in conjunction with reported birth difficulties may predispose to KLS.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

Sex specific effect of gut microbiota on the risk of psychiatric disorders: A Mendelian randomisation study and PRS analysis using UK Biobank cohort

OBJECTIVE

The relationships between gut microbiota and brain-related diseases/traits remains not fully understood.

METHOD

A two-stage study was performed to investigate the relationships between gut microbiota and brain-related diseases/traits, and evaluate the potential sex specific effects of gut microbiota. In discovery stage, we systematically scanned the relationships between 515 brain-related diseases/traits and gut microbiota through two-sample Mendelian randomisation analysis. Using ∼500,000 individuals derived from the UK Biobank, polygenetic risk scoring (PRS) analysis was performed to validate the associations detected in discovery stage. To evaluate the potential sex-specific effect of gut microbiota on brain-related disorders, PRS analysis was conducted in female and male, respectively.

RESULTS

After systematically scanning diseases or traits, 41 of the 515 brain-related diseases/traits were identified to be associated with gut microbiota, such as Neuroticism score (P_2-MR = 0.0018), worrier/anxious feelings (P_2-MR = 0.0013), Suffer from 'nerves' (P_2-MR = 0.0062) and Nervous feelings (P_2-MR = 0.0158). 5 of 41 brain-related diseases or traits were successfully validated in UK Biobank, such as Neuroticism score (P_UK = 0.0024, P_UK-female = 0.0063, P_UK-male = 0.1142), Nervous feelings (P_UK = 0.0043, P_UK-female = 0.0115, P_UK-male = 0.1670) and Worrier/anxious feelings (P_UK = 0.0166, P_UK-female = 0.0196, P_UK-male = 0.2930).

CONCLUSION

Our results suggest that gut microbiota contributed more to brain-related diseases or traits in females than in males.Key pointsA two-stage study was performed to investigate the relationships between gut microbiota and brain-related diseases/traits.Using the individuals derived from the UK Biobank, polygenetic risk scoring analysis was performed to validate the associations detected in the discovery stage.Our results suggest that gut microbiota contributed more to brain-related diseases or traits in females than in males.

Mood Stabilizers in Psychiatric Disorders and Mechanisms Learnt from In Vitro Model Systems

Bipolar disorder (BD) and schizophrenia are psychiatric disorders that manifest unusual mental, behavioral, and emotional patterns leading to suffering and disability. These disorders span heterogeneous conditions with variable heredity and elusive pathophysiology. Mood stabilizers such as lithium and valproic acid (VPA) have been shown to be effective in BD and, to some extent in schizophrenia. This review highlights the efficacy of lithium and VPA treatment in several randomized, controlled human trials conducted in patients suffering from BD and schizophrenia. Furthermore, we also address the importance of using induced pluripotent stem cells (iPSCs) as a disease model for mirroring the disease's phenotypes. In BD, iPSC-derived neurons enabled finding an endophenotype of hyperexcitability with increased hyperpolarizations. Some of the disease phenotypes were significantly alleviated by lithium treatment. VPA studies have also reported rescuing the Wnt/β-catenin pathway and reducing activity. Another significant contribution of iPSC models can be attributed to studying the molecular etiologies of schizophrenia such as abnormal differentiation of patient-derived neural stem cells, decreased neuronal connectivity and neurite number, impaired synaptic function, and altered gene expression patterns. Overall, despite significant advances using these novel models, much more work remains to fully understand the mechanisms by which these disorders affect the patients' brains.

ANK3 Gene Polymorphism Rs10994336 Influences Executive Functions by Modulating Methylation in Patients With Bipolar Disorder

Background: A large body of evidence suggests that epigenetic modification including DNA methylation plays a critical role in BD's pathogenesis while the identification of methylation quantitative trait loci (meQTLs) shed light on the interpretation of the function of genetic variants in non-coding regions. The intronic single nucleotide polymorphism (SNP) rs10994336 within the ANK3 has emerged as one of the most replicated risk variants for bipolar disorder (BD) in genome-wide association studies. Whether rs10994336 functions as a meQTL to mediate the association between genotype and phenotype remains unclear.

Method: A total of 154 patients with BD and 181 healthy controls (HC) were recruited. The genotypes of rs10994336 and methylation levels of CpG sites within ANK3 were tested. Executive functions were assessed using a computerized version of the Wisconsin Card Sorting Test (WCST).

Results: Bipolar disorder patients with the risk-T allele of rs10994336 scored lower on tests of executive function compared to homozygous CC carriers, after controlling for age, gender, and education level. No significant difference was found in HC individuals. The risk-T allele is associated with a lower methylation level of CpG site cg02172182 in HC after multiple corrections and replicated in the BD group in the same direction. Further mediation analysis revealed that the cg02172182 methylation significantly mediated the association between the polymorphism rs10994336 and PE index of WCST in patients with BD.

Conclusion: Our study suggests that BD-related genetic variant rs10994336 in ANK3 impacts executive functions by modulating ANK3 methylation, supporting the theory that methylation acts as a mediator between genotype and phenotype.

Identify Consistent Cross-Modality Imaging Genetic Patterns via Discriminant Sparse Canonical Correlation Analysis

Sparse canonical correlation analysis (SCCA) is a bi-multivariate technique used in imaging genetics to identify complex multi-SNP-multi-QT associations. However, the traditional SCCA algorithm has been designed to seek a linear correlation between the SNP genotype and brain imaging phenotype, ignoring the discriminant similarity information between within-class subjects in brain imaging genetics association analysis. In addition, multi-modality brain imaging phenotypes are extracted from different perspectives and imaging markers from the same region consistently showing up in multimodalities may provide more insights for the mechanistic understanding of diseases. In this paper, a novel multi-modality discriminant SCCA algorithm (MD-SCCA) is proposed to overcome these limitations as well as to improve learning results by incorporating valuable discriminant similarity information into the SCCA algorithm. Specifically, we first extract the discriminant similarity information between within-class subjects by the sparse representation. Second, the discriminant similarity information is enforced within SCCA to construct a discriminant SCCA algorithm (D-SCCA). At last, the MD-SCCA algorithm is adopted to fully explore the relationships among different modalities of different subjects. In experiments, both synthetic dataset and real data from the Alzheimer's Disease Neuroimaging Initiative database are used to test the performance of our algorithm. The empirical results have demonstrated that the proposed algorithm not only produces improved cross-validation performances but also identifies consistent cross-modality imaging genetic biomarkers.

Identify Complex Imaging Genetic Patterns via Fusion Self-Expressive Network Analysis

In the brain imaging genetic studies, it is a challenging task to estimate the association between quantitative traits (QTs) extracted from neuroimaging data and genetic markers such as single-nucleotide polymorphisms (SNPs). Most of the existing association studies are based on the extensions of sparse canonical correlation analysis (SCCA) for the identification of complex bi-multivariate associations, which can take the specific structure and group information into consideration. However, they often take the original data as input without considering its underlying complex multi-subspace structure, which will deteriorate the performance of the following integrative analysis. Accordingly, in this paper, the self-expressive property is exploited for the reconstruction of the original data before the association analysis, which can well describe the similarity structure. Specifically, we first apply the within-class similarity information to construct self-expressive networks by sparse representation. Then, we use the fusion method to iteratively fuse the self-expressive networks from multi-modality brain phenotypes into one network. Finally, we calculate the imaging genetic association based on the fused self-expressive network. We conduct the experiments on both single-modality and multi-modality phenotype data. Related experimental results validate that our method can not only better estimate the potential association between genetic markers and quantitative traits but also identify consistent multi-modality imaging genetic biomarkers to guide the interpretation of Alzheimer's disease.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Candidate risk genes for bipolar disorder are highly conserved during evolution and highly interconnected

OBJECTIVES

Bipolar disorder (BPD) is a highly heritable psychiatric disorder whose genetic complexity and pathogenetic mechanisms are still being unraveled. The main goal of this work was to characterize BPD risk-gene candidates (identified by Nurnberger et al., JAMA Psychiatry 71:657, 2014, and Stahl et al., Nat. Genet. 51:793, 2019) with respect to their evolutionary conservation, associated phenotypes, and extent of gene-gene interactions.

METHODS

Database searches and BLAST were used to identify homologous counterparts of human BPD risk genes in C. elegans, zebrafish, and Drosophila. Phenotypes associated with the C. elegans genes were annotated and searched. With GeneMANIA, we characterized and quantified gene-gene interactions among members of the BPD gene set in comparison to randomly chosen gene sets of the same size.

RESULTS

BPD risk genes are highly conserved across species and are enriched for essential genes and genes associated with lethality and altered life span. They are significantly more interactive with each other in comparison to random genes. We identified syntenic blocks of risk genes, which provided potential insights into molecular pathways and co-morbidities associated with BPD including coronary disease, obesity, and decreased life expectancy.

CONCLUSIONS

BPD risk genes appear to be special in terms of their degree of conservation, interconnectedness, and pleiotropic effects that extend beyond a role in brain function. Key hub genes or pleiotropic regulatory components may represent attractive targets for future drug discovery.

Impaired blood-brain barrier in the microbiota-gut-brain axis: Potential role of bipolar susceptibility gene TRANK1

Bipolar disorder (BD) is a common psychiatric illness with high prevalence and disease burden. Accumulating susceptibility genes for BD have been identified in recent years. However, the exact functions of these genes remain largely unknown. Despite its high heritability, gene and environment interaction is commonly accepted as the major contributing factor to BD pathogenesis. Intestine microbiota is increasingly recognized as a critical environmental factor for human health and diseases via the microbiota‐gut‐brain axis. BD individuals showed altered diversity and compositions in the commensal microbiota. In addition to pro‐inflammatory factors, such as interleukin‐6 and tumour necrosis factor‐α, type 1 interferon signalling pathway is also modulated by specific intestinal bacterial strains. Disruption of the microbiota‐gut‐brain axis contributes to peripheral and central nervous system inflammation, which accounts for the BD aetiology. Administration of type 1 interferon can induce the expression of TRANK1, which is associated with elevated circulating biomarkers of the impaired blood‐brain barrier in BD patients. In this review, we focus on the influence of intestine microbiota on the expression of bipolar gene TRANK1 and propose that intestine microbiota‐dependent type 1 interferon signalling is sufficient to induce the over‐expression of TRANK1, consequently causing the compromise of BBB integrity and facilitating the entrance of inflammatory mediators into the brain. Activated neuroinflammation eventually contributes to the occurrence and development of BD. This review provides a new perspective on how gut microbiota participate in the pathogenesis of BD. Future studies are needed to validate these assumptions and develop new treatment targets for BD.

Independent replications and integrative analyses confirm TRANK1 as a susceptibility gene for bipolar disorder

Genetic analyses for bipolar disorder (BD) have achieved prominent success in Europeans in recent years, whereas its genetic basis in other populations remains relatively less understood. We herein report that the leading risk locus for BD in European genome-wide association studies (GWAS), the single-nucleotide polymorphism (SNP) rs9834970 near TRANK1 at 3p22 region, is also genome-wide significantly associated with BD in a meta-analysis of four independent East Asian samples including 5748 cases and 65,361 controls (p = 2.27 × 10^-8, odds ratio = 1.136). Expression quantitative trait loci (eQTL) analyses and summary data-based Mendelian randomization (SMR) analyses in multiple human brain samples suggest that lower TRANK1 mRNA expression is a principal BD risk factor explaining its genetic risk signals at 3p22. We also identified another SNP rs4789 in the 3' untranslated region (3'UTR) of TRANK1 showing stronger eQTL associations as well as genome-wide significant association with BD. Despite the relatively unclear neuronal function of TRANK1, our mRNA expression analyses in the human brains and in rat primary cortical neurons reveal that genes highly correlated with TRANK1 are significantly enriched in the biological processes related to dendritic spine, synaptic plasticity, axon guidance and circadian entrainment, and are also more likely to exhibit strong associations in psychiatric GWAS (e.g., the CACNA1C gene). Overall, our results support that TRANK1 is a potential BD risk gene. Further studies elucidating its roles in this illness are needed.

The functional polymorphisms linked with interleukin-1β gene expression are associated with bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a severe psychiatric illness attributable to multifactorial risk components (e.g. environmental stimuli, neuroinflammation, etc.), and genetic variations affecting these risk components are considered pivotal predisposing factors. The interleukin-1β (IL-1β) gene and its protein product have been repeatedly highlighted in the pathogenesis of BD. As functional polymorphisms and haplotypes linked with IL-1β mRNA expression have been reported, whether they are correlated with the risk of developing BD remains to be tested.

METHODS

To examine whether variations in the IL-1β gene locus confer genetic risk of BD, we recruited 930 BD patients and 912 healthy controls for the current study. All subjects were Han Chinese, and were age- and gender-matched. We tested seven functional single nucleotide polymorphisms (SNPs) spanning the IL-1β gene and one haplotype composed of three SNPs for their associations with risk of BD.

RESULTS

We found that the functional SNPs in the promoter region of IL-1β gene were significantly associated with risk of BD. The haplotype analyses further supported the involvement of IL-1β promoter SNPs in BD. The BD risk SNPs in our study have been previously reported to predict higher IL-1β levels in the brain and peripheral blood, which is consistent with the clinical observation of elevated IL-1β levels in the lymphocytes or peripheral blood of patients with BD compared with healthy subjects.

CONCLUSION

Our results support the contention that IL-1β is likely a risk gene for BD, and further investigations on this gene may promote our understanding and clinical management of this illness.

Bipolar multiplex families have an increased burden of common risk variants for psychiatric disorders

Multiplex families with a high prevalence of a psychiatric disorder are often examined to identify rare genetic variants with large effect sizes. In the present study, we analysed whether the risk for bipolar disorder (BD) in BD multiplex families is influenced by common genetic variants. Furthermore, we investigated whether this risk is conferred mainly by BD-specific risk variants or by variants also associated with the susceptibility to schizophrenia or major depression. In total, 395 individuals from 33 Andalusian BD multiplex families (166 BD, 78 major depressive disorder, 151 unaffected) as well as 438 subjects from an independent, BD case/control cohort (161 unrelated BD, 277 unrelated controls) were analysed. Polygenic risk scores (PRS) for BD, schizophrenia (SCZ), and major depression were calculated and compared between the cohorts. Both the familial BD cases and unaffected family members had higher PRS for all three psychiatric disorders than the independent controls, with BD and SCZ being significant after correction for multiple testing, suggesting a high baseline risk for several psychiatric disorders in the families. Moreover, familial BD cases showed significantly higher BD PRS than unaffected family members and unrelated BD cases. A plausible hypothesis is that, in multiplex families with a general increase in risk for psychiatric disease, BD development is attributable to a high burden of common variants that confer a specific risk for BD. The present analyses demonstrated that common genetic risk variants for psychiatric disorders are likely to contribute to the high incidence of affective psychiatric disorders in the multiplex families. However, the PRS explained only part of the observed phenotypic variance, and rare variants might have also contributed to disease development.

Whole Genome Interpretation for a Family of Five

Although best practices have emerged on how to analyse and interpret personal genomes, the utility of whole genome screening remains underdeveloped. A large amount of information can be gathered from various types of analyses via whole genome sequencing including pathogenicity screening, genetic risk scoring, fitness, nutrition, and pharmacogenomic analysis. We recognize different levels of confidence when assessing the validity of genetic markers and apply rigorous standards for evaluation of phenotype associations. We illustrate the application of this approach on a family of five. By applying analyses of whole genomes from different methodological perspectives, we are able to build a more comprehensive picture to assist decision making in preventative healthcare and well-being management. Our interpretation and reporting outputs provide input for a clinician to develop a healthcare plan for the individual, based on genetic and other healthcare data.

Prenatal testosterone triggers long-term behavioral changes in male zebra finches: unravelling the neurogenomic mechanisms

BACKGROUND

Maternal hormones, like testosterone, can strongly influence developing offspring, even generating long-term organizational effects on adult behavior; yet, the mechanisms facilitating these effects are still unclear. Here, we experimentally elevated prenatal testosterone in the eggs of zebra finches (Taeniopygia guttata) and measured male aggression in adulthood along with patterns of neural gene expression (RNA-seq) and DNA methylation (MethylC-Seq) in two socially relevant brain regions (hypothalamus and nucleus taenia of the amygdala). We used enrichment analyses and protein-protein interaction networks to find candidate processes and hub genes potentially affected by the treatment. We additionally identified differentially expressed genes that contained differentially methylated regions.

RESULTS

We found that males from testosterone-injected eggs displayed more aggressive behaviors compared to males from control eggs. Hundreds of genes were differentially expressed, particularly in the hypothalamus, including potential aggression-related hub genes (e.g., brain derived neurotrophic factor). There were also enriched processes with well-established links to aggressive phenotypes (e.g., somatostatin and glutamate signaling). Furthermore, several highly connected genes identified in protein-protein interaction networks also showed differential methylation, including adenylate cyclase 2 and proprotein convertase 2.

CONCLUSIONS

These results highlight genes and processes that may play an important role in mediating the effects of prenatal testosterone on long-term phenotypic outcomes, thereby providing insights into the molecular mechanisms that facilitate hormone-mediated maternal effects.

A linkage and exome study of multiplex families with bipolar disorder implicates rare coding variants of ANK3 and additional rare alleles at 10q11-q21

BACKGROUND

Bipolar disorder is a highly heritable psychiatric condition for which specific genetic factors remain largely unknown. In the present study, we used combined whole-exome sequencing and linkage analysis to identify risk loci and dissect the contribution of common and rare variants in families with a high density of illness.

METHODS

Overall, 117 participants from 15 Australian extended families with bipolar disorder (72 with affective disorder, including 50 with bipolar disorder type I or II, 13 with schizoaffective disorder-manic type and 9 with recurrent unipolar disorder) underwent whole-exome sequencing. We performed genome-wide linkage analysis using MERLIN and conditional linkage analysis using LAMP. We assessed the contribution of potentially functional rare variants using a genebased segregation test.

RESULTS

We identified a significant linkage peak on chromosome 10q11-q21 (maximal single nucleotide polymorphism = rs10761725; exponential logarithm of the odds [LODexp] = 3.03; empirical p = 0.046). The linkage interval spanned 36 protein-coding genes, including a gene associated with bipolar disorder, ankyrin 3 (ANK3). Conditional linkage analysis showed that common ANK3 risk variants previously identified in genome-wide association studies - or variants in linkage disequilibrium with those variants - did not explain the linkage signal (rs10994397 LOD = 0.63; rs9804190 LOD = 0.04). A family-based segregation test with 34 rare variants from 14 genes under the linkage interval suggested rare variant contributions of 3 brain-expressed genes: NRBF2 (p = 0.005), PCDH15 (p = 0.002) and ANK3 (p = 0.014).

LIMITATIONS

We did not examine non-coding variants, but they may explain the remaining linkage signal.

CONCLUSION

Combining family-based linkage analysis with next-generation sequencing data is effective for identifying putative disease genes and specific risk variants in complex disorders. We identified rare missense variants in ANK3, PCDH15 and NRBF2 that could confer disease risk, providing valuable targets for functional characterization.