- Micro report
- Open Access
Postmortem gene expression profiles in the habenulae of suicides: implication of endothelial dysfunction in the neurovascular system
Molecular Brain volume 15, Article number: 48 (2022)
The habenula (Hb) is an epithalamic structure that links multiple forebrain areas with the mid/hindbrain monoaminergic systems. As an anti-reward center, it has been implicated in the etiology of various neuropsychiatric disorders, particularly those associated with dysregulated reward circuitry. In this regard, Hb has been proposed as a therapeutic target for treatment-resistant depression associated with a higher risk of suicide. Therefore, we aimed to gain insight into the molecular signatures of the Hb in association with suicide in individuals with major depression. Postmortem gene expression analysis identified 251 differentially expressed genes (DEGs) in the Hb tissue of suicides in comparison with Hb tissues from neurotypical individuals. Subsequent bioinformatic analyses using single-cell transcriptome data from the mouse Hb showed that the levels of a subset of endothelial cell-enriched genes encoding cell–cell junctional complex and plasma membrane-associated proteins, as well as the levels of their putative upstream transcriptional regulators, were significantly affected in suicides. Although our findings are based on a limited number of samples, the present study suggests a potential association of endothelial dysfunction in the Hb with depression and suicidal behavior.
The habenula (Hb) is an epithalamic structure located at the dorsomedial posterior end of the thalamus in mammals. It serves as a neuroanatomical bridge, linking multiple forebrain structures such as the frontal cortex, basal forebrain, and hypothalamus, with monoaminergic nuclei in the mid- and hindbrain. The Hb has been recently recognized for its role in encoding aversive stimuli as an anti-reward center, and its implications in neuropsychiatric disorders associated with dysregulated reward circuitry, such as depressive disorders, schizophrenia, and substance use disorders [1, 2]. Notably, Hb has been proposed as a therapeutic target for treatment-resistant major depressive disorder with a higher risk of suicide [3, 4], and its altered connectivity has been associated with suicidal ideation . Our previous study also revealed that several genes involved in cholinergic signaling exhibit decreased mRNA expression levels in Hb tissue of suicides . To gain systemic insights into the molecular signatures of human Hb underlying suicide with major depression, the present study examined genome-wide gene expression profiles in the Hb from suicide completers compared to unaffected controls.
For this purpose, we compared postmortem gene expression patterns between groups using a microarray platform. All subjects were male Caucasians, and groups (suicides and controls, n = 10 for each group) were matched for age, postmortem interval (PMI), and pH; Age: 44.20 ± 3.35 years (mean ± SE) for controls and 40.50 ± 4.35 for suicides; PMI: 26.45 ± 4.95 h for controls and 36.10 ± 4.32 for suicides; pH: 6.671 ± 0.047 for controls and 6.427 ± 0.095 for suicides (see Additional file 1 for experimental procedures and Additional file 2: Table S1 for more detailed subject information). Among 20,085 detected genes, RNA expression levels of 796 were significantly different between groups (p < 0.05), and 251 genes at the criteria of p < 0.05, and > 1.2-fold change (FC) were selected as differentially expressed genes (DEGs) for subsequent bioinformatic analyses (Additional file 2: Table S2). ARAP3, PHLDB2, PAQR5, PMAIP1, and GPR116 were the top-five upregulated DEGs, and NPNT, HSPA1A/B, SLC22A4, PAX3, and UBC were the most downregulated genes (Fig. 1a). We then compared our DEGs with disease-associated genes curated in the PsyGeNET database  and found that 28 DEGs are implicated in several psychiatric disorders (Fig. 1b). Gene-disease associations clearly show that these genes are associated with depressive disorder, bipolar disorder, substance abuse (alcohol, cannabis, or cocaine), and schizophrenia (Fig. 1c).
Next, we categorized the DEGs by cell type enrichment in their expression using single-cell level gene expression profiles obtained from murine Hb . Interestingly, cell type enrichment analysis raised a possibility that endothelial cells (ECs) constituting brain microvasculature could be one of the most affected cell types in the suicidal Hb; 59.3% of DEGs with cell-type-restricted expression were EC-enriched genes in the mouse Hb region (Fig. 1d and Additional file 2: Table S3). Among 155 DEGs found in the murine Hb dataset, 47 DEGs were highly enriched or even exclusively expressed in ECs, and most of them, except HSPA1A and HSPB1, exhibited significantly increased expression levels in the suicidal Hb (Fig. 1e, upper panel). The EC-enriched DEGs included genes encoding key components of cell–cell junctions such as CDH5, ESAM, LSR, and OCLN, along with those encoding plasma membrane-associated proteins such as ABCG2, ADCY4, ATP10A, CD34, FLT1, ICAM2, KDR, PTPRB, SEMA3G, SLC2A1, SLC38A5, and TIE1, implying possible alterations in endothelial functions, particularly involving blood–brain barrier (BBB) permeability . Protein–protein interaction (PPI) analysis also supported this notion. A PPI network was constructed to elucidate functional interactions between DEGs (221 nodes and 272 edges), and subsequent subcluster network analysis proposed the top-30 genes consisting of 20 EC-enriched genes and 2 genes with relatively higher expression levels in ECs, which also included genes linked with BBB permeability (Fig. 1f).
In good agreement with cell-type enrichment and PPI analyses for the identified DEGs, gene ontology (GO) analysis also suggested endothelial dysfunction. Among the featured GO biological processes, three terms such as “maintenance of BBB,” “positive regulation of vasculature development” and “endothelial cell differentiation,” are intricately linked with the roles of ECs in controlling BBB and/or microvasculature in brain tissues (Fig. 1G and Additional file 2: Table S4). Concerted regulation of EC-enriched genes suggested that a subset of upstream transcription factors (TFs) required for EC differentiation and function may underlie suicide-associated alterations in Hb gene expression profiles. To test this idea, we conducted gene-TF co-expression analysis using the ARCHS4 database (https://maayanlab.cloud/archs4/) and found that putative target genes of several EC-specific TFs were significantly enriched among the DEGs (Fig. 1h and Additional file 2: Table S5). Among the top-five co-expressed TFs, BCL3B, SOX18 and EPAS1 mRNA expression levels were significantly higher in the suicides than controls (Additional file 2: Table S2). It should be noted that the SOXF subgroup of the SOX (SRY-box transcription factor) family such as SOX7, SOX17, and SOX18 is well known as EC-enriched TFs required for vasculogenesis  and plays a major role in controlling BBB permeability through reciprocal crosstalk with canonical Wnt/β-catenin signaling . In addition, another EC-enriched TF, EPAS1 (endothelial PAS domain protein 1) was predicted to be activated in the suicidal Hb (Fig. 1e); its mRNA expression is known to be inducible in response to hypoxic/ischemic stimuli in brain tissues . One may suppose that the altered expression of EC-specific genes may represent different EC proportions between groups. However, it should be noted that the mRNA expression levels of well-established EC markers such as CLDN5, ICAM1, VCAM1, PTPRC, FLT4, and TEK were not significantly different between the groups (Additional file 2: Table S6). Therefore, it is more reasonable that the differential mRNA expression of a subset of EC-enriched genes may represent altered microvasculature functions in the Hb.
Several postmortem studies have implicated brain endothelial dysfunction in discrete brain areas including the prefrontal cortex, nucleus accumbens, and amygdala in the etiology of certain types of depressive disorders [13, 14]; altered gene expression profiles responsible for impaired BBB resistance and enhanced angiogenic signaling have been particularly linked to depressive disorders . As the Hb is one of the brain regions where ECs are abundant , it could represent a major site of action of EC-enriched and psychiatric disease-associated genes such as CD34, CDH5, FN1, HSPA1A, KDR, TGM2, and VWF (Fig. 1b and 1e). Therefore, it is plausible that DEGs promoting BBB permeability and vasculature development in the Hb may underlie depressive phenotypes, including suicidal behavior.
In conclusion, the present study identified a set of genes whose mRNA expression is significantly altered in the Hb of suicides, and suggests that endothelial dysfunction, presumably linked to Hb microvasculatures and BBB integrity, may play a role in the etiology of suicide associated with major depression. Although the present study is primarily based on transcriptome profiles of bulk Hb tissues with a limited number of samples, we provide putative molecular and cellular targets, implying possible importance of Hb ECs in mood regulation and suicidal behavior. The functional consequences of EC dysfunction on neuroglial properties in the Hb, as well as emotional abnormalities, should be further delineated.
Availability of data and materials
The datasets used and analyzed in the present study are available from the corresponding authors upon reasonable request. Raw microarray data have been submitted to the Gene Expression Omnibus (GEO) repository (Accession number: GSE199536).
Differentially expressed gene
Oligodendrocyte lineage cell
GO biological process
Hu H, Cui Y, Yang Y. Circuits and functions of the lateral habenula in health and in disease. Nat Rev Neurosci. 2020;21(5):277–95.
Metzger M, Souza R, Lima LB, Bueno D, Gonçalves L, Sego C, et al. Habenular connections with the dopaminergic and serotonergic system and their role in stress-related psychiatric disorders. Eur J Neurosci. 2021;53(1):65–88.
Sartorius A, Henn FA. Deep brain stimulation of the lateral habenula in treatment resistant major depression. Med Hypotheses. 2007;69(6):1305–8.
Zanos P, Gould TD. Mechanisms of ketamine action as an antidepressant. Mol Psychiatry. 2018;23(4):801–11.
Qiao D, Zhang A, Sun N, Yang C, Li J, Zhao T, et al. Altered static and dynamic functional connectivity of habenula associated with suicidal ideation in first-episode, drug-naïve patients with major depressive disorder. Front Psychiatry. 2020;11:608197.
Han S, Yang SH, Kim JY, Mo S, Yang E, Song KM, et al. Down-regulation of cholinergic signaling in the habenula induces anhedonia-like behavior. Sci Rep. 2017;7(1):900.
Gutiérrez-Sacristán A, Grosdidier S, Valverde O, Torrens M, Bravo À, Piñero J, et al. PsyGeNET: a knowledge platform on psychiatric disorders and their genes. Bioinformatics. 2015;31(18):3075–7.
Hashikawa Y, Hashikawa K, Rossi MA, Basiri ML, Liu Y, Johnston NL, et al. Transcriptional and spatial resolution of cell types in the mammalian habenula. Neuron. 2020;106(5):743-758.e5.
Munji RN, Soung AL, Weiner GA, Sohet F, Semple BD, Trivedi A, et al. Profiling the mouse brain endothelial transcriptome in health and disease models reveals a core blood–brain barrier dysfunction module. Nat Neurosci. 2019;22(11):1892–902.
Lilly AJ, Lacaud G, Kouskoff V. SOXF transcription factors in cardiovascular development. Semin Cell Dev Biol. 2017;63:50–7.
Corada M, Orsenigo F, Bhat GP, Conze LL, Breviario F, Cunha SI, et al. Fine-tuning of Sox17 and canonical Wnt coordinates the permeability properties of the blood–brain barrier. Circ Res. 2019;124(4):511–25.
Chavez JC, Baranova O, Lin J, Pichiule P. The transcriptional activator hypoxia inducible factor 2 (HIF-2/EPAS-1) regulates the oxygen-dependent expression of erythropoietin in cortical astrocytes. J Neurosci. 2006;26(37):9471–81.
Empana JP, Boutouyrie P, Lemogne C, Jouven X, van Sloten TT. Microvascular contribution to late-onset depression: mechanisms, current evidence, association with other brain diseases, and therapeutic perspectives. Biol Psychiatry. 2021;90(4):214–25.
Zeng D, He S, Ma C, Wen Y, Song W, Xu Q, et al. Network-based approach to identify molecular signatures in the brains of depressed suicides. Psychiatry Res. 2020;294: 113513.
Le Foll B, French L. Transcriptomic characterization of the human habenula highlights drug metabolism and the neuroimmune system. Front Neurosci. 2018;12:742.
We appreciate Dr. Gustavo Turecki (Douglas Hospital Research Center, McGill University, Montreal, Canada) for providing habenular RNA samples and personal communications for the results. We would like to thank Editage (www.editage.co.kr) for English language editing.
This work was supported by National Research Foundation of Korea (NRF) grants funded by the Korea Government Ministry of Science and ICT (NRF-2019M3C7A1032764 to GHS and NRF-2017M3C7A1079692 to HK), and by the Korea University Research Grants to GHS and HK.
Ethics approval and consent to participate
Postmortem habenular tissue was obtained from the Douglas-Bell Canada Brain Bank (Douglas Mental Health University Institute, Verdun, Quebec) under the guidelines and approval of the Research Ethics Committee of the Douglas Hospital Research Centre. This study was approved by the Korea University Anam Hospital Institutional Review Board (Approval No. ED13217).
Consent for publication
All authors have approved the contents of this manuscript and have provided consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Materials and methods including subject information, RNA isolation, microarray analysis, and subsequent bioinformatic analyses.
: Table S1. Subject information; Table S2. Full list of the main 251 DEGs; Table S3. Cell-type enrichment analysis using single-cell level gene expression profiles obtained from murine habenula; Table S4. Top-10 GO-BP terms and related genes; Table S5. Top-10 ARCHS4 co-expressed TFs and putative target DEGs; Table S6. List of EC markers that did not exhibit significant changes in mRNA expression levels.
About this article
Cite this article
Kim, H.J., Yoo, H., Kim, J.Y. et al. Postmortem gene expression profiles in the habenulae of suicides: implication of endothelial dysfunction in the neurovascular system. Mol Brain 15, 48 (2022). https://doi.org/10.1186/s13041-022-00934-7
- Gene expression