pi :: Psychiatry Investigation

Cho, Lee, Woo, Choi, Greenwood, and Kelsoe: CDH13 and HCRTR2 May Be Associated with Hypersomnia Symptom of Bipolar Depression: A Genome-Wide Functional Enrichment Pathway Analysis

Brief Report

Psychiatry Investigation 2015;12(3):402-407.

Published online: July 06, 2015

DOI: https://doi.org/10.4306/pi.2015.12.3.402

CDH13 and HCRTR2 May Be Associated with Hypersomnia Symptom of Bipolar Depression: A Genome-Wide Functional Enrichment Pathway Analysis

Chul-Hyun Cho¹, Heon-Jeong Lee¹, Hyun Goo Woo², Ji-Hye Choi², Tiffany A. Greenwood³, John R. Kelsoe^3,⁴

¹Department of Psychiatry, Korea University College of Medicine, Seoul, Republic of Korea.

²Department of Physiology, Ajou University School of Medicine, Suwon, Republic of Korea.

³Department of Psychiatry, University of California San Diego, La Jolla, CA, USA.

⁴San Diego VA Healthcare System, San Diego, CA, USA.

Correspondence: Heon-Jeong Lee, MD, PhD. Department of Psychiatry, Korea University College of Medicine, 73 Inchon-ro, Seongbuk-gu, Seoul 136-705, Republic of Korea. Tel: +82-2-920-6721, Fax: +82-2-929-7679, leehjeong@korea.ac.kr

Received August 21, 2014 Revised November 17, 2014 Accepted November 17, 2014

(open-access, http://creativecommons.org/licenses/by-nc/3.0/):

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Although bipolar disorder is highly heritable, the identification of specific genetic variations is limited because of the complex traits underlying the disorder. We performed a genome-wide association study of bipolar disorder using a subphenotype that shows hypersomnia symptom during a major depressive episode. We investigated a total of 2,191 cases, 1,434 controls, and 703,012 single nucleotide polymorphisms (SNPs) in the merged samples obtained from the Translational Genomics Institute and the Genetic Association Information Network. The gene emerging as the most significant by statistical analysis was rs1553441 (odds ratio=0.4093; p=1.20×10^-5; Permuted p=6.0×10^-6). However, the 5×0^-8 threshold for statistical significance required in a genome-wide association study was not achieved. The functional enrichment pathway analysis showed significant enrichments in the adhesion, development-related, synaptic transmission-related, and cell recognition-related pathways. For further evaluation, each gene of the enriched pathways was reviewed and matched with genes that were suggested to be associated with psychiatric disorders by previous genetic studies. We found that the cadherin 13 and hypocretin (orexin) receptor 2 genes may be involved in the hypersomnia symptom during a major depressive episode of bipolar disorder.

Keywords: Genome-wide association study, Bipolar disorder, Hypersomnia, Functional enrichment pathway analysis, Bipolar depression

INTRODUCTION

Bipolar disorder (BD) patients suffer from both a low physical quality and low mental quality of life, even during an euthymic period.1 BD is known to be highly heritable; however, the identification of specific genetic variations resulted in limited findings.2 3 One reason for this low success rate in identifying genetic variations relies in the nature and complex traits of BD.4 As a consequence, dividing BD into subgroups according to clinical subphenotypes is a reliable alternative approach for future genetic studies in BD.5

Several studies suggested that dysregulation of the circadian rhythm is associated with the pathophysiology of BD.6 7 In general, BD patients have common specific circadian rhythm-related clinical manifestations such as: diurnal mood variation, periodicity of manic-depressive episodes, activity changes, alteration of the secretion of hormones and other endogenous substances, and sleep disturbances.8 9 10 Sleep disturbances may have a negative effect on emotional regulation.11 According to these results, it can be assumed that circadian rhythm dysregulation is very closely associated with BD.12

We hypothesized that a subphenotype group that shows hypersomnia symptom during a major depressive episode of bipolar disorder may represent a genetically distinct subtype of BD. We tested this hypothesis using a genome-wide association (GWA) analysis between subjects that showed hypersomnia symptom during a major depressive episode of bipolar disorder (hypersomnia symptom of bipolar depression, HBD) and subjects that did not show hypersomnia symptom during a major depressive episode of bipolar disorder (non-hypersomnia symptom of bipolar depression, NHBD). The subject samples were of European ancestry and were genotyped as part of the Genetic Association Information Network (GAIN) by the Bipolar Genome Study (BiGS).

METHODS

Subject ascertainment

Prior to genotyping, which was part of the BiGS, the unrelated bipolar I disorder subjects of European ancestry were selected from those collected by the National Institute of Mental Health Genetics Initiative for Bipolar Disorder. All subjects provided written informed consent in accordance to protocols from local institutional review boards. The subjects were interviewed using protocols from the Diagnostic Interview for Genetic Studies (DIGS).13 Information was obtained from family informants and medical records. The information was reviewed along with the interview by a panel of experienced clinicians to obtain a final best-estimate diagnosis. Control subjects were selected after they were certified through a National Institute of Mental Health-supported contract between Dr. Pablo Gejman and Knowledge Networks, Inc.14 All subjects donated a blood sample and were given medical questionnaires. The selected controls were matched for gender and ethnicity with the BD cases. Control subjects who had a history of BD, psychosis, or recurrent major depression were excluded from our study. The subject samples as part of the GAIN were obtained by Dr. Kelsoe who is a member of the BiGS.

Genotyping and cleaning

The first set of samples was genotyped at the Broad Institute, as part of GAIN, using the Affymetrix SNP Array 6.0 (Affymetrix; Santa Clara, CA, US) 1M SNP array. We obtained a total of 1,001 BD cases, 1,033 controls, and 724,067 single nucleotide polymorphisms (SNPs). These were available for analysis following an extensive quality control (QC) process.15 The QC process eliminated all individuals with >10% missing data, SNPs with poor allele clustering, duplicate errors, minor allele frequencies <0.05, and significant deviation from Hardy-Weinberg equilibrium at p<10^-6. The second set of samples was genotyped similarly to the first set of samples at the Translational Genomics Institute (TGEN) and underwent a comparable QC. From this set of samples we obtained 1,190 BD cases, 401 controls, and 728,187 SNPs available for analysis. An additional round of QC was performed on the merged samples from GAIN and TGEN. This merge resulted in a set of 703,012 SNPs that passed the imposed QC process.

Phenotypes

As part of the DIGS interview, bipolar I disorder subjects were queried as to whether they had a hypersomnia symptom during major depressive episodes. According to the answer, subjects were categorized either into the HBD or the NHBD group. Those who answered 'Unknown' to the question about hypersomnia symptom were categorized as the missing group. After filtering, there were 263 BD subjects in the HBD group and 112 subjects in the NHBD group.

Association analyses

To assess genetic factors contributing to HBD, we performed a genome-wide case-only analysis of HBD versus NHBD. This association analysis was performed using a logistic regression using PLINK16 with a covariance adjustment for sex and age. Adaptive permutations were performed to find the empirical significance of the results using PLINK.16

SNP imputation

Missing SNPs were imputed using the IMPUTE2 tool17 and the CEU panel of HapMap 3+1,000 Genomes Pilot haplotypes as a reference. The imputed SNPs were used for the functional enrichment pathway analysis.

Functional enrichment pathway analysis

Gene ontology analysis was performed on the gene sets harboring the identified SNPs (p<0.005) evaluated using the DAVID software.18 19 The enriched gene functions were identified from the HBD versus NHBD analysis using an enrichment score. To go further in our investigation, we made a list of genes of the enriched pathways and compared them with various genes that were suggested to be related to psychiatric disorders by previous genetic studies. We searched the available databases in PubMed (http://www.ncbi.nlm.nih.gov/pubmed/) and in the National Human Genome Research Institute (http://www.genome.gov/) for genes of enriched pathways. Some of the search terms included: the name of each gene of enriched pathways and "psychiatric disorder or depression or bipolar disorder or schizophrenia or circadian rhythm". After the review and matching, we selected genes from those of enriched pathways that were described as susceptible genes in psychiatric disorders by previous genetic studies.

RESULTS

Eighty associated SNPs with p<10^-4 were identified in this present genome-wide case-only analysis. We found that rs1553441 is the gene that appeared to be the most significant after statistical analysis (odds ratio=0.4093, p value of 1.20×10^-5, Permuted p=6.0×10^-6) (Supplementary Figure 1 in the online-only Data Supplement). This SNP is located within a region of the gene encoding the phospholipase D family, member 5 (PLD5) on chromosome 1q43. However, the most significant gene of this study did not reach the p value of 5×10^-8, which is the threshold for statistical significance in a GWA.16

To determine whether the associated SNPs showed enrichment in certain functional pathways, we performed a genome-wide functional enrichment pathway analysis using the imputed data. Gene ontology functions were identified for all genes that contained SNPs associated with a significance level of p<0.005. The functional enrichment analysis of HBD versus NHBD showed a significant enrichment of the adhesion, development-related, synaptic transmission-related, and cell recognition-related pathways (Supplementary Table 1 in the online-only Data Supplement). We can speculate that the altered functions of these pathways may contribute to the development of hypersomnia symptom during a major depressive episode of BP.

We aimed to investigate the genes of the enriched pathways in more detail. At first, we reviewed the results of previous genetic studies of psychiatric disorders found in the PubMed database (http://www.ncbi.nlm.nih.gov/pubmed/) and the database from the National Human Genome Research Institute (http://www.genome.gov/). To search the databases we used search terms such as: the name of each gene of enriched pathways and "psychiatric disorder or major depressive disorder or bipolar disorder or schizophrenia or circadian rhythm". In accordance with previously reported genetic findings, we investigated the significant genes of those identified from categories that showed significant enrichment by comparison to the results of previous genetic studies. We selected the genes of enriched pathways that matched those that have been suggested in the literature to be related to psychiatric disorders (Table 1). We ultimately found two matching genes: cadherin 13 (CDH13) and hypocretin (orexin) receptor 2 (HCRTR2). CDH13 is located in 16q24.2 and encodes a member of the cadherin superfamily. The encoded protein acts as a negative regulator of axon growth during neural differentiation.20 21 In the present pathway analysis, CDH13 was identified as a member of genes thought to be involved in cell adhesion, biological adhesion, cell-cell adhesion, cell projection organization, and cell motion. HCRTR2 is located in 6p12 and encodes a G protein-coupled receptor involved in regulation of feeding and sleep behavior.22 In the present pathway analysis, HCRTR2 was identified as a member of genes thought to be involved in cell-cell signaling, synaptic transmission, transmission of nerve impulse, and neurological system process.

DISCUSSION

We were not able to identify significant genes in a GWA study that evaluated HBD versus NHBD. As a consequence, we focused on the pathways identified in genome-wide functional enrichment pathway analysis that showed significant enriched function. Specifically, we evaluated the gene lists of these enriched pathways, and compared them with genes that were previously suggested in the literature to be related to psychiatric disorders. This procedure was found to be helpful in finding genes that are possibly related to HBD, although this is an unfamiliar method and the results do not always provide an explanation on how genes are related.

CDH13 is highly expressed in various brain regions such as the cerebral cortex, medulla, thalamus, and midbrain.20 It has been suggested that CDH13 is associated with various drug abuse-related phenotypes,23 24 25 26 27 28 especially with respect to comorbid depression and alcohol dependence.29 Disorders related to substance use, particularly alcohol dependence, show high comorbidity with BD.30 In addition, the presence of substance abuse disorders seriously affects the course and prognosis of BD.31 CDH13 expression also showed an association with adiponectin levels32 33 and depression.34 Previous studies investigated adult depressed patients that showed decreased circulating levels of adiponectin.35 36 37 One explanation for this finding could be that adiponectin is related to the Hypothalamic-Pituitary-Adrenal (HPA) axis. Several studies suggested a correlation between adiponectin and stress and development of various psychiatric disorders.38 Furthermore, TNF-α, IL-10, and IL-6, previously shown to be regulated by adiponectin, also showed a correlation with BD.39 40 TNF-α is suggested to produce an increased slow-wave sleep and a symptom of sleepiness.41 Orexin, also known as hypocretin-A and B, neuromodulatory peptides secreted from orexin neurons, signal through orexin receptor 1 and orexin receptor 2 (HCRTR1 and HCRTR2), which are G protein-coupled receptors.42 Orexin-containing neurons project out to various monoaminergic nuclei of the brain, including the locus coeruleus (norepinephrine), the raphe nuclei (5-hydroxytryptamine, 5-HT), and the ventral tegmental area (dopamine). Therefore, the orexin-containing neurons can modify the monoaminergic neurons.43 Furthermore, the presence of a functional positive feedback loop between the orexin and monoaminergic neurons has been suggested.44 Orexin is implicated in diverse functions such as feeding, drinking, the sleep-wake cycle, hormone secretion, autonomic function, drug addiction, and reward. To date, most research has focused on the role of orexin in depression compared to other mood disorders, and several studies suggested that orexin plays a significant role in depression.45 In particular, several studies have reported a role for orexin in sleep disturbance,46 47 reward system,48 feeding behavior,49 50 hippocampal neuronal plasticity,51 52 53 and monoamines.54 55 It is important to note that BP is an important portion of depressive symptoms. Moreover, atypical depressive symptoms, characterized by hypersomnia, hyperphagia (or weight gain), and leaden paralysis, are more common in patients with BP than unipolar depression.56 We speculate that representative symptoms of atypical depression, such as hypersomnia and hyperphagia, may be closely related to orexin's function. Moreover, a balanced effect of orexin's action on either the HCRTR1 or the HCRTR2 receptor is important in achieving an anti-depressant- or pro-depressant-like effect.57 Taken together, these previous findings suggest that hypersomnia symptoms during a major depressive episode of bipolar disorder are related to the functions of CDH13 and HCRTR2.

We have explored the genetic association of BD using the HBD subphenotype to focus on a more homogenous group of subjects that present similar clinical courses. Although the present GWA of HBD versus NHBD did not produce significant results, the functional enrichment pathway analysis identified significant enrichments of the adhesion, development-related, synaptic transmission-related, and cell recognition-related pathways. In addition, we found that CDH13 and HCRTR2 show a potential correlation with HBD and we used a matching process based on the results of previous genetic studies performed on psychiatric disorders. To date, only several studies have been reported about this correlation. In the future, more studies focusing on the correlation between these pathways or several genes and HBD are needed to validate the results reported in the present study.

Acknowledgments

Bipolar Genome Study Co-authors: John R. Kelsoe, Tiffany A. Greenwood, Caroline M. Nievergelt, Rebecca McKinney, Paul D. Shilling, Erin N. Smith - University of California, San Diego, CA, USA; Nicholas J. Schork, Cinnamon S. Bloss - Scripps Translational Science Institute, La Jolla, CA, USA; John I. Nurnberger, Jr., Howard J. Edenberg, Tatiana Foroud, Daniel L. Koller - Indiana University, Indianapolis, IN, USA; Elliot S. Gershon, Chunyu Liu, Judith A. Badner - University of Chicago, Chicago, IL, USA; William A. Scheftner - Rush University Medical Center, Chicago, IL, USA; William B. Lawson, Evaristus A. Nwulia, Maria Hipolito - Howard University, Washington, D.C., USA; William Coryell - University of Iowa, Iowa City, IA, USA; John Rice - Washington University, St. Louis, MO, USA; William Byerley - University of California, San Francisco, CA, USA; Francis J. McMahon, Thomas G. Schulze - National Institute of Mental Health Intramural Research Program, Bethesda, MD, USA; Wade H. Berrettini - University of Pennsylvania, Philadelphia, PA, USA; James B. Potash, Peter P. Zandi, Pamela B. Mahon - Johns Hopkins School of Medicine, Baltimore, MD, USA; Melvin G. McInnis, Sebastian Zöllner, Peng Zhang - University of Michigan, Ann Arbor, MI, USA; David W. Craig, Szabolics Szelinger - The Translational Genomics Research Institute, Phoenix, AZ, USA; Thomas B. Barrett - Portland Veterans Affairs Medical Center, Portland, OR, USA

This work was supported by grants from the National Institute of Mental Health (NIMH) and National Human Genome Research Institute (NHGRI) to JRK (MH078151, MH081804, and MH059567supplement), and by the Genetic Association Information Network (GAIN). HJL was supported by the Korea Health 21 R&D Project funded by the Ministry of Health & Welfare, Republic of Korea (HI11C1901) and the Future Environmental R&D grant funded by the Korea Environmental Industry and Technology Institute (No. RE201206020).

Supplementary Materials

The online-only Data Supplement is available with this article at http://dx.doi.org/10.4306/pi.2015.12.3.402.

Supplementary Figure 1

Genome-wide association results for hypersomnia bipolar depression. The Manhattan plot shows the results of susceptible loci relevant to hypersomnia bipolar depression. The chromosomal position is shown along the X-axis, whereas the -log (p value) for each single nucleotide polymorphism is shown along the Y-axis. The horizontal line indicates the p<10^-4 significance threshold. The arrows in the figure indicate the position of possible relevant genes above the set threshold.

pi-12-402-s001.pdf

Supplementary Table 1

Enrichment findings from the genome-wide functional enrichment pathway analysis of hypersomnia bipolar depression

pi-12-402-s002.pdf

References

1. Gutierrez-Rojas L, Gurpegui M, Ayuso-Mateos JL, Gutierrez-Ariza JA, Ruiz-Veguilla M, Jurado D. Quality of life in bipolar disorder patients: a comparison with a general population sample. Bipolar Disord 2008;10:625–634. PMID: 18657247.

2. Baum AE, Akula N, Cabanero M, Cardona I, Corona W, Klemens B, et al. A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder. Mol Psychiatry 2008;13:197–207. PMID: 17486107.

3. Smith EN, Koller DL, Panganiban C, Szelinger S, Zhang P, Badner JA, et al. Genome-wide association of bipolar disorder suggests an enrichment of replicable associations in regions near genes. PLoS Genet 2011;7:e1002134PMID: 21738484.

4. Risch N, Merikangas K. The future of genetic studies of complex human diseases. Science 1996;273:1516–1517. PMID: 8801636.

5. McQueen MB, Devlin B, Faraone SV, Nimgaonkar VL, Sklar P, Smoller JW, et al. Combined analysis from eleven linkage studies of bipolar disorder provides strong evidence of susceptibility loci on chromosomes 6q and 8q. Am J Hum Genet 2005;77:582–595. PMID: 16175504.

6. Kripke DF, Nievergelt CM, Joo E, Shekhtman T, Kelsoe JR. Circadian polymorphisms associated with affective disorders. J Circadian Rhythms 2009;7:2PMID: 19166596.

7. McCarthy MJ, Nievergelt CM, Kelsoe JR, Welsh DK. A survey of genomic studies supports association of circadian clock genes with bipolar disorder spectrum illnesses and lithium response. PLoS One 2012;7:e32091PMID: 22384149.

8. Kennedy SH, Kutcher SP, Ralevski E, Brown GM. Nocturnal melatonin and 24-hour 6-sulphatoxymelatonin levels in various phases of bipolar affective disorder. Psychiatry Res 1996;63:219–222. PMID: 8878318.

9. Kripke DF, Mullaney DJ, Atkinson M, Wolf S. Circadian rhythm disorders in manic-depressives. Biol Psychiatry 1978;13:335–351. PMID: 667233.

10. Leibenluft E, Albert PS, Rosenthal NE, Wehr TA. Relationship between sleep and mood in patients with rapid-cycling bipolar disorder. Psychiatry Res 1996;63:161–168. PMID: 8878312.

11. Harvey AG, Mullin BC, Hinshaw SP. Sleep and circadian rhythms in children and adolescents with bipolar disorder. Dev Psychopathol 2006;18:1147–1168. PMID: 17064432.

12. Mitterauer B. Clock genes, feedback loops and their possible role in the etiology of bipolar disorders: an integrative model. Med Hypotheses 2000;55:155–159. PMID: 10904433.

13. Nurnberger JI Jr, Blehar MC, Kaufmann CA, York-Cooler C, Simpson SG, Harkavy-Friedman J, et al. Diagnostic interview for genetic studies. Rationale, unique features, and training. NIMH Genetics Initiative. Arch Gen Psychiatry 1994;51:849–859. discussion 863-864PMID: 7944874.

14. Sanders AR, Levinson DF, Duan J, Dennis JM, Li R, Kendler KS, et al. The Internet-based MGS2 control sample: self report of mental illness. Am J Psychiatry 2010;167:854–865. PMID: 20516154.

15. Smith EN, Bloss CS, Badner JA, Barrett T, Belmonte PL, Berrettini W, et al. Genome-wide association study of bipolar disorder in European American and African American individuals. Mol Psychiatry 2009;14:755–763. PMID: 19488044.

16. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–575. PMID: 17701901.

17. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet 2009;5:e1000529PMID: 19543373.

18. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 2009;37:1–13. PMID: 19033363.

19. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009;4:44–57. PMID: 19131956.

20. Takeuchi T, Misaki A, Liang SB, Tachibana A, Hayashi N, Sonobe H, et al. Expression of T-cadherin (CDH13, H-Cadherin) in human brain and its characteristics as a negative growth regulator of epidermal growth factor in neuroblastoma cells. J Neurochem 2000;74:1489–1497. PMID: 10737605.

21. Redies C. Cadherins in the central nervous system. Prog Neurobiol 2000;61:611–648. PMID: 10775799.

22. Mignot E. A commentary on the neurobiology of the hypocretin/orexin system. Neuropsychopharmacology 2001;25(5 Suppl):S5–S13. PMID: 11682267.

23. Drgon T, Montoya I, Johnson C, Liu QR, Walther D, Hamer D, et al. Genome-wide association for nicotine dependence and smoking cessation success in NIH research volunteers. Mol Med 2009;15:21–27. PMID: 19009022.

24. Johnson C, Drgon T, Walther D, Uhl GR. Genomic regions identified by overlapping clusters of nominally-positive SNPs from genome-wide studies of alcohol and illegal substance dependence. PLoS One 2011;6:e19210PMID: 21818250.

25. Hart AB, Engelhardt BE, Wardle MC, Sokoloff G, Stephens M, de Wit H, et al. Genome-wide association study of d-amphetamine response in healthy volunteers identifies putative associations, including cadherin 13 (CDH13). PLoS One 2012;7:e42646PMID: 22952603.

26. Liu QR, Drgon T, Johnson C, Walther D, Hess J, Uhl GR. Addiction molecular genetics: 639,401 SNP whole genome association identifies many "cell adhesion" genes. Am J Med Genet B Neuropsychiatr Genet 2006;141B:918–925. PMID: 17099884.

27. Treutlein J, Cichon S, Ridinger M, Wodarz N, Soyka M, Zill P, et al. Genome-wide association study of alcohol dependence. Arch Gen Psychiatry 2009;66:773–784. PMID: 19581569.

28. Uhl GR, Drgon T, Johnson C, Walther D, David SP, Aveyard P, et al. Genome-wide association for smoking cessation success: participants in the Patch in Practice trial of nicotine replacement. Pharmacogenomics 2010;11:357–367. PMID: 20235792.

29. Edwards AC, Aliev F, Bierut LJ, Bucholz KK, Edenberg H, Hesselbrock V, et al. Genome-wide association study of comorbid depressive syndrome and alcohol dependence. Psychiatr Genet 2012;22:31–41. PMID: 22064162.

30. Merikangas KR, Akiskal HS, Angst J, Greenberg PE, Hirschfeld RM, Petukhova M, et al. Lifetime and 12-month prevalence of bipolar spectrum disorder in the National Comorbidity Survey replication. Arch Gen Psychiatry 2007;64:543–552. PMID: 17485606.

31. Salloum IM, Thase ME. Impact of substance abuse on the course and treatment of bipolar disorder. Bipolar Disord 2000;2:269–280. PMID: 11249805.

32. Jee SH, Sull JW, Lee JE, Shin C, Park J, Kimm H, et al. Adiponectin concentrations: a genome-wide association study. Am J Hum Genet 2010;87:545–552. PMID: 20887962.

33. Wu Y, Li Y, Lange EM, Croteau-Chonka DC, Kuzawa CW, McDade TW, et al. Genome-wide association study for adiponectin levels in Filipino women identifies CDH13 and a novel uncommon haplotype at KNG1-ADIPOQ. Hum Mol Genet 2010;19:4955–4964. PMID: 20876611.

34. Terracciano A, Tanaka T, Sutin AR, Sanna S, Deiana B, Lai S, et al. Genome-wide association scan of trait depression. Biol Psychiatry 2010;68:811–817. PMID: 20800221.

35. Leo R, Di Lorenzo G, Tesauro M, Cola C, Fortuna E, Zanasi M, et al. Decreased plasma adiponectin concentration in major depression. Neurosci Lett 2006;407:211–213. PMID: 16973279.

36. Narita K, Murata T, Hamada T, Takahashi T, Kosaka H, Sudo S, et al. Adiponectin multimer distribution, not absolute amount of plasma, correlates with depression severity in healthy elderly subjects. Prog Neuropsychopharmacol Biol Psychiatry 2008;32:124–127. PMID: 17761380.

37. Diniz BS, Teixeira AL, Campos AC, Miranda AS, Rocha NP, Talib LL, et al. Reduced serum levels of adiponectin in elderly patients with major depression. J Psychiatr Res 2012;46:1081–1085. PMID: 22633396.

38. Taylor VH, Macqueen GM. The role of adipokines in understanding the associations between obesity and depression. J Obes 2010;2010:pii: 748048

39. Modabbernia A, Taslimi S, Brietzke E, Ashrafi M. Cytokine alterations in bipolar disorder: a meta-analysis of 30 studies. Biol Psychiatry 2013;74:15–25. PMID: 23419545.

40. Fernandez-Real JM, Lopez-Bermejo A, Casamitjana R, Ricart W. Novel interactions of adiponectin with the endocrine system and inflammatory parameters. J Clin Endocrinol Metab 2003;88:2714–2718. PMID: 12788878.

41. Kapsimalis F, Basta M, Varouchakis G, Gourgoulianis K, Vgontzas A, Kryger M. Cytokines and pathological sleep. Sleep Med 2008;9:603–614. PMID: 18024171.

42. Sakurai T, Amemiya A, Ishii M, Matsuzaki I, Chemelli RM, Tanaka H, et al. Orexins and orexin receptors: a family of hypothalamic neuropeptides and G protein-coupled receptors that regulate feeding behavior. Cell 1998;92:573–585. PMID: 9491897.

43. Peyron C, Tighe DK, van den Pol AN, de Lecea L, Heller HC, Sutcliffe JG, et al. Neurons containing hypocretin (orexin) project to multiple neuronal systems. J Neurosci 1998;18:9996–10015. PMID: 9822755.

44. Eriksson KS, Sergeeva OA, Haas HL, Selbach O. Orexins/hypocretins and aminergic systems. Acta Physiol (Oxf) 2010;198:263–275. PMID: 19566795.

45. Nollet M, Leman S. Role of orexin in the pathophysiology of depression: potential for pharmacological intervention. CNS Drugs 2013;27:411–422. PMID: 23657787.

46. Allard JS, Tizabi Y, Shaffery JP, Manaye K. Effects of rapid eye movement sleep deprivation on hypocretin neurons in the hypothalamus of a rat model of depression. Neuropeptides 2007;41:329–337. PMID: 17590434.

47. Feng P, Vurbic D, Wu Z, Strohl KP. Brain orexins and wake regulation in rats exposed to maternal deprivation. Brain Res 2007;1154:163–172. PMID: 17466285.

48. Nocjar C, Zhang J, Feng P, Panksepp J. The social defeat animal model of depression shows diminished levels of orexin in mesocortical regions of the dopamine system, and of dynorphin and orexin in the hypothalamus. Neuroscience 2012;218:138–153. PMID: 22626650.

49. Lutter M, Krishnan V, Russo SJ, Jung S, McClung CA, Nestler EJ. Orexin signaling mediates the antidepressant-like effect of calorie restriction. J Neurosci 2008;28:3071–3075. PMID: 18354010.

50. Mahler SV, Smith RJ, Moorman DE, Sartor GC, Aston-Jones G. Multiple roles for orexin/hypocretin in addiction. Prog Brain Res 2012;198:79–121. PMID: 22813971.

51. Ito N, Yabe T, Gamo Y, Nagai T, Oikawa T, Yamada H, et al. I.c.v. administration of orexin-A induces an antidepressive-like effect through hippocampal cell proliferation. Neuroscience 2008;157:720–732. PMID: 18952152.

52. Howell OW, Doyle K, Goodman JH, Scharfman HE, Herzog H, Pringle A, et al. Neuropeptide Y stimulates neuronal precursor proliferation in the post-natal and adult dentate gyrus. J Neurochem 2005;93:560–570. PMID: 15836615.

53. Ito N, Yabe T, Nagai T, Oikawa T, Yamada H, Hanawa T. A possible mechanism underlying an antidepressive-like effect of Kososan, a Kampo medicine, via the hypothalamic orexinergic system in the stress-induced depression-like model mice. Biol Pharm Bull 2009;32:1716–1722. PMID: 19801833.

54. Adidharma W, Leach G, Yan L. Orexinergic signaling mediates light-induced neuronal activation in the dorsal raphe nucleus. Neuroscience 2012;220:201–207. PMID: 22710065.

55. Pail G, Huf W, Pjrek E, Winkler D, Willeit M, Praschak-Rieder N, et al. Bright-light therapy in the treatment of mood disorders. Neuropsychobiology 2011;64:152–162. PMID: 21811085.

56. Mitchell PB, Wilhelm K, Parker G, Austin MP, Rutgers P, Malhi GS. The clinical features of bipolar depression: a comparison with matched major depressive disorder patients. J Clin Psychiatry 2001;62:212–216. quiz 217PMID: 11305713.

57. Scott MM, Marcus JN, Pettersen A, Birnbaum SG, Mochizuki T, Scammell TE, et al. Hcrtr1 and 2 signaling differentially regulates depression-like behaviors. Behav Brain Res 2011;222:289–294. PMID: 21377495.

Table 1

Enrichment findings of the genome-wide functional enrichment pathway analysis in hypersomnia symptom of bipolar depression

A gene of interest is the result of a matching procedure based on previous genetic studies on psychiatric disorders. There were various genes identified as enriched genes of interest; however, we omitted all genes that did not overlap with the genetic findings from previously reported psychiatric genetic studies. Genes of annotation cluster 4 did not contain any genes that overlapped with previous studies, and was therefore omitted from this table. Supplementary Table 1 lists the full version of these results