CRY2 Is Associated with Depression

Abstract

Background

Abnormalities in the circadian clockwork often characterize patients with major depressive and bipolar disorders. Circadian clock genes are targets of interest in these patients. CRY2 is a circadian gene that participates in regulation of the evening oscillator. This is of interest in mood disorders where a lack of switch from evening to morning oscillators has been postulated.

Principal Findings

We observed a marked diurnal variation in human CRY2 mRNA levels from peripheral blood mononuclear cells and a significant up-regulation (P = 0.020) following one-night total sleep deprivation, a known antidepressant. In depressed bipolar patients, levels of CRY2 mRNA were decreased (P = 0.029) and a complete lack of increase was observed following sleep deprivation. To investigate a possible genetic contribution, we undertook SNP genotyping of the CRY2 gene in two independent population-based samples from Sweden (118 cases and 1011 controls) and Finland (86 cases and 1096 controls). The CRY2 gene was significantly associated with winter depression in both samples (haplotype analysis in Swedish and Finnish samples: OR = 1.8, P = 0.0059 and OR = 1.8, P = 0.00044, respectively).

Conclusions

We propose that a CRY2 locus is associated with vulnerability for depression, and that mechanisms of action involve dysregulation of CRY2 expression.

About the Authors:

Catharina Lavebratt

* E-mail: [email protected]

Affiliation: Neurogenetics Unit, Department of Molecular Medicine and Surgery, Karolinska Institutet at Karolinska University Hospital Solna, Stockholm, Sweden

Louise K. Sjöholm

Affiliation: Neurogenetics Unit, Department of Molecular Medicine and Surgery, Karolinska Institutet at Karolinska University Hospital Solna, Stockholm, Sweden

Pia Soronen

Affiliation: Public Health Genomics Unit, Department of Chronic Disease Prevention, National Institute for Health and Welfare, Helsinki, Finland

Tiina Paunio

Affiliation: Public Health Genomics Unit, Department of Chronic Disease Prevention, National Institute for Health and Welfare, Helsinki, Finland

Marquis P. Vawter

Affiliation: Functional Genomics Laboratory, Department of Psychiatry and Human Behavior, University of California Irvine School of Medicine, Irvine, California, United States of America

William E. Bunney

Affiliation: Functional Genomics Laboratory, Department of Psychiatry and Human Behavior, University of California Irvine School of Medicine, Irvine, California, United States of America

Rolf Adolfsson

Affiliation: Division of Psychiatry, Department of Clinical Sciences, University of Umeå, Umeå, Sweden

Yvonne Forsell

Affiliation: Department of Public Health Science, Karolinska Institutet at Karolinska University Hospital Solna, Stockholm, Sweden

Joseph C. Wu

Affiliation: Functional Genomics Laboratory, Department of Psychiatry and Human Behavior, University of California Irvine School of Medicine, Irvine, California, United States of America

John R. Kelsoe

Affiliation: Department of Psychiatry, University of California San Diego School of Medicine, La Jolla, California, United States of America

Timo Partonen

Affiliation: Mood, Depression and Suicidal Behaviour Unit, Department of Mental Health and Substance Abuse Services, National Institute for Health and Welfare, Helsinki, Finland

Martin Schalling

Affiliation: Neurogenetics Unit, Department of Molecular Medicine and Surgery, Karolinska Institutet at Karolinska University Hospital Solna, Stockholm, Sweden

Introduction

Rhythms that approximate the 24-hour day-night cycle or light-dark transitions are called circadian. Abnormalities of the circadian pacemaker system is often seen in mood disorders (i.e. major depressive and bipolar disorders) as evidenced with sleep stage abnormalities and the clinical efficacy of sleep deprivation [1] as well as with the therapeutic mechanisms of lithium [2] that is prescribed primarily for bipolar disorder. Approximately one tenth of all mood disorders follow a seasonal pattern and area hence categorized as seasonal affective disorder (SAD) [3]. Season-bound mood episodes may occur in both depressive and bipolar disorders and emerge in any season, but the most common type is winter depression, a condition in which major depressive episodes routinely occur in the wintertime and remit the following summer [4]. In one study as much as 93% of the winter depression cases had a diagnosis of bipolar disorder [5] although other studies tend to show a greater proportion of unipolar recurrent major depression among winter depression cases [4]. Depressive episodes in winter depression are highly recurrent and appear to be clearly endogenous as there is no psychosocial factor that would account for their onset. For bipolar type 1 disorder, the heritability estimate is very high [6]. Therefore, winter depression, whether part of major depressive disorder or bipolar disorder, as well as bipolar type 1 disorder, with or without a seasonal pattern, provide excellent models for studying the molecular mechanisms of mood disorder [7], [8].

Advances and delays in phase (the location within a cycle at a particular time) and reduced amplitudes (intensities of the oscillations) have been reported in adult patients [9]–[11], while children with winter depression have circadian rhythms that are accurate in phase but low in amplitude [12]. This agrees with the hypothesis that not only phase shifts but also amplitude attenuations [13] contribute to the pathogenesis, and suggests that there are circadian clockwork abnormalities having relevance to the daily reset and synchronization. This is further supported by findings of the seasonal changes in sensitivity to light exposure in winter depression, these patients having supersensitivity to light in terms of melatonin suppression during winter [14], abnormal melatonin levels in patients with seasonal or non-seasonal depressive disorder [15], and abnormalities in circadian alignments in patients with bipolar disorder [16]. Therefore, the internal misalignment (i.e. the sleep-wake cycle is no longer in phase with the circadian rhythms) may account for the pathogenesis of mood disorders in general. Light exposures and sleep manipulations in patients with mood disorder are tools for exploration and elucidation of the mechanisms driving the circadian and seasonal clockworks.

CRY proteins [17] differ from many transcription factors that take part in the circadian clockwork, since they have no PAS domains. This suggests that they are unique in the core of the circadian clockwork where they act [18]. Both CRY1 and CRY2 operate in the retina and non-visual light detection pathways in a manner that is important for the internal alignment [19]–[24]. Of these two, CRY2 is highly expressed in the brain in particular and has a dose-dependent inhibitory effect on the activated ARNTL, whereas both CRY2 and CRY1 repress all four combinations of the ARNTL (ARNTL2) – CLOCK (NPAS2) protein heterodimers that act as transcriptional activators in the core of the circadian clockwork [25]. These basic findings give a rationale for the CRY2 gene as a target of high interest and relevance in our study.

Here, we assessed CRY2 gene expression in eight healthy volunteers for 48 hours and in 13 patients with bipolar disorder before and after a one-night sleep deprivation, and we report that CRY2 mRNA levels are reduced and unresponsive to sleep deprivation in depressed patients with bipolar disorder. To determine whether CRY2 genetic variation is associated with depression, we analyzed circadian clock gene variants in two separate population-based samples, a Swedish sample and a Finnish sample, and report that the CRY2 gene is associated with winter depression.

Results

CRY2 Expression Analysis

CRY2 mRNA levels displayed marked diurnal variation in healthy controls (n = 8). In addition, total sleep deprivation induced a 2.0-fold increase in CRY2 mRNA levels (P = 0.020) in the controls (Figure 1). The control gene GAPDH mRNA levels did not show diurnal variation and were not changed by the effect of sleep deprivation (P = 0.66).

[Figure omitted. See PDF.]

Figure 1. CRY2 gene expression in healthy volunteers and bipolar disorder patients.

CRY2 mRNA expression was not changed by sleep deprivation in bipolar disorder patients (triangle) and was significantly decreased compared to healthy volunteers (square) during sleep deprivation. Blood draw 1 (x-axis) was at 7 p.m., then there was a blood draw every 6 hours. Sleep deprivation was at time points #6, #7, #8 and #9. Bars indicate mean and error bars indicate SEM.

https://doi.org/10.1371/journal.pone.0009407.g001

In samples from patients in a depressive state of bipolar disorder (n = 13), sleep deprivation did not induce any increase in CRY2 mRNA levels. CRY2 mRNA expression was significantly decreased in the samples from patients as compared with the samples from healthy controls during sleep deprivation (P = 0.041 for the diagnosis and sleep deprivation interaction using ANOVA; P = 0.029 for the difference between patients and controls during sleep deprivation in post hoc test). GAPDH showed similar mRNA levels in the samples from controls to those in the bipolar disorder samples (P = 0.12).

Genetic Association Analysis

Analyses were performed on two samples of clinical cases with a diagnosis of winter depression (DSM-IV) and healthy screened controls (Swedish sample: 118 cases and 1011 controls, Finnish sample: 86 cases and 1096 controls). CRY2 allele frequency was significantly associated with winter depression in the Swedish sample, as shown with three SNPs (rs10838524 risk allele A, rs10838527 risk allele G, and rs3824872 risk allele A) having OR = 1.6, OR = 2.1 and OR = 1.8 (P = 0.0017, P = 0.00074 and 0.00070, respectively). The association to rs10838524 was confirmed in the Finnish sample (OR = 1.7, P = 0.0020) (Table 1) and indicated that homozygosity of the minor allele (overall minor allele frequency of about 0.5) was increased in the cases in both population samples (Table 2). The associated SNPs (rs10838527 and rs3824872) were present in LD block 1 (rs7123390, rs10838527, rs3824872) in both samples (Figures 2A and 2B).

[Figure omitted. See PDF.]

Figure 2. The suggested vulnerability locus in CRY2.

LD structure of the SNPs analyzed in A) the Swedish sample, and B) the Finnish sample. C) Scheme of the suggested location of the functional variant in the Finnish sample, in the Swedish sample, and the common locus relevant under the assumption of the same functional variant in the two samples. An open end of locus box (seen for the Swedish sample) indicates that this is not the locus border. The numbers in the squares of A and B) represent the pair-wise D' value, empty squares stand for D' = 1. Pink-red color indicates a pair-wise LOD ≥2 with redness proportional to D'. Haplotype blocks are formed if 95% of comparisons are “strong LD”, i.e. the 95% CI of D' is within [0.7–0.98].

https://doi.org/10.1371/journal.pone.0009407.g002

[Figure omitted. See PDF.]

Table 1. SNP allele frequency association analysis.

https://doi.org/10.1371/journal.pone.0009407.t001

[Figure omitted. See PDF.]

Table 2. Genotype association analysis of the SNPs with suggestive allele frequency association.

https://doi.org/10.1371/journal.pone.0009407.t002

Four haplotypes were identified and a significant difference in distribution of haplotypes between the cases and controls was found for the Swedish sample (χ2 = 8.7, P = 0.034) and the Finnish sample (χ2 = 14.6, P = 0.0022; see Supplemental data Table S1). A haplotype of the risk alleles GGA was more frequent in the Swedish cases compared with controls (OR = 1.7, P = 0.012) and even more so when adding the rs10838524 risk allele to the haplotype to be AGGA (OR = 1.8, P = 0.0059; Table 3). The block 1 (rs7123390, rs10838527, rs3824872) haplotype GAC which in the Swedish sample has the protective alleles from rs10838527 and rs3824872 was a risk haplotype in the Finnish sample (OR = 1.8, P = 0.00010), and similarly so when adding the rs10838524 allele G to the haplotype to be GGAC (OR = 1.8, P = 0.00044; Table 3). The Finnish data suggest a vulnerability locus increasing the risk for winter depression upstream rs10838527. This was further supported by the fact that changing only the rs7123390 allele in the Finnish risk haplotype (GAC) turned it into a protective haplotype (AAC) in the Finnish sample (OR = 0.61, P = 0.012), and the fact that the AGAC was non-existent of in the samples. Considering only the Swedish data the vulnerability locus could reside on either side of, but not exactly at rs7123390. If assuming the same functional variation (although probably opposite risk alleles in Swedes and Finns), the Swedish and the Finnish samples together hence suggest the vulnerability locus to reside in-between rs7123390 and rs10838527 or upstream rs7123390 (Figure 2C).

[Figure omitted. See PDF.]

Table 3. Haplotype association analysis.

https://doi.org/10.1371/journal.pone.0009407.t003

Discussion

Our key findings herein are that CRY2 gene variation and expression levels are associated with depression. CRY2 mRNA levels are lowered in blood mononuclear cells from depressed patients with bipolar disorder after total sleep deprivation in comparison to healthy controls, and CRY2 gene variation was associated with winter depression in both Swedish and Finnish patients. Differences in CRY2 risk haplotypes were observed between the Swedes and Finns. Overall there are genetic differences between these populations [26]. The fact that Swedes and Finns have different histories in terms of going through genetic bottle-necks that could have affected selection in operation at the CRY2 locus with regard to depression vulnerability, would possibly explain the differences in CRY2 risk haplotypes seen between the two samples, and explain a putative allelic heterogeneity (different vulnerability alleles in Swedes and Finns) of the functional CRY2 polymorphism. Assuming the same functional variant in the Swedish and the Finnish samples, the two samples would assist in narrowing the interval of that functional variant (Fig. 2C). The risk haplotypes spanned from CRY2 intron 1 to 3′UTR, and pointed at a location of a potential function variation somewhere upstream exon 12 in 3′UTR (rs10838527) from the Finnish data, with a distance from the intron 1 marker rs10838524 to rs10838527 corresponding to ∼33 kb (NCBI build 130). The vulnerability locus from the Swedish data overlapped the Finnish locus, but for at the intron 7 marker rs7123390, and extended downstream rs10838527.

Human CRY2 mRNA levels have been demonstrated to undergo circadian oscillation in fibroblasts [27] and in hematopoietic stem cells [28]. Our results herein add on these findings and demonstrate the effect of sleep deprivation on the circadian oscillations of CRY2 mRNA in human peripheral blood mononuclear cells. Our data indicate that depression in bipolar disorder is related to lowered levels of CRY2 mRNA. Sleep deprivation led to CRY2 mRNA increase in controls whereas depressed bipolar individuals were non-responsive.

CRY2 has been purified from human cells and its properties have been described [29], but the understanding of CRY2 regulation and CRY2 functions still need elucidation. Earlier studies have demonstrated that CRY2 transcription is driven by the canonical circadian circuit but through a unique mechanism of action not yet characterized [30]–[32]. CRY2 is unique among the canonical genes of the circadian clock, since experimental findings of its effects on the circadian clockwork violate the predictions from a theory on the roles of PER1, PER2, CRY1 and CRY2 [33]–[36]. Deletion of Cry2 gene lengthened the circadian period by approximately 48 min (from 23.7 to 24.5 hours) [37] and in fact strengthened circadian amplitudes and restored the lost rhythms in Per2 mutant mice [38], [39]. These findings fit in the view that activation of Per1 and Per2 genes occurs in the morning, whereas activation of Cry1 and Cry2 genes occurs in the evening [40]. Abnormalities in CRY2 regulation would therefore leave the morning oscillator intact, which agrees with findings in patients with winter depression [41] and those with bipolar disorder [42]. However, disruption of CRY2 regulation may compromise the evening oscillator and the subsequent proper oscillator reactions to light-dark transitions [43]. Cry2 null mutants are less decelerated, or are more accelerated, by light exposure than Cry1 null mutants [36]. The former phenotype does match with that seen in winter depression [44] and, to some extent, in bipolar type 1 disorder [45].

The action of melatonin on Cry1/Cry2 expression is hypothesized to form the basis of refractory reactions to stimulation with light [40], and the subsequent subsensitive or supersensitive responses to light exposure. Phase control through the melatonin-guided interval of Per1/Per2 to Cry1/Cry2 expression peaks may have relevance to mood disorders, as the internal alignment and the melatonin signal (amplitude and phase) are abnormal in patients with bipolar disorder and those with winter depression in particular [15], [16], [46]. Our findings now suggest that a dawn component [47], [48] and a dusk component (herein), i.e. PER2 and CRY2 variants respectively, are affected in winter depression. Such influence may well contribute to phase angle differences [49] and entrainment errors [11] that have been found in patients with winter depression.

To sum up, the phenotype of the Cry2 knockout mice together with mechanistic data are in line with the finding in this report suggesting that CRY2 has a role in winter depression.

Our current findings add to earlier gene expression studies in blood that have reported a number of potential peripheral biomarkers for bipolar disorder [50], [51], and earlier gene expression findings in brain indicating that clock genes are implicated in a model of bipolar disorder, the first report being [52], followed by [53] for CRY2. Circadian genes with a polymorphism previously reported to be associated to human bipolar disorder include CLOCK [54], and NR1D1, ARNTL and PER3 [55]–[57]. Moreover, ARNTL, RORA, RORB and RXRG were associated with bipolar disorder in a meta-analysis integrating data from genome-wide association studies and human and animal model expression studies [51]. PER2, NPAS2 and ARNTL genetic variation have been found associated with winter depression [47], [48], [58], and recent experimental data demonstrated these three proteins to regulate transcription of MAOA and to have subsequent influence on depressive behavior [59].

Other circadian genes were tested for association to winter depression in parallel with CRY2. The CRY2 finding reached significance even after correction for multiple testing considering other circadian gene SNPs tested for association to winter depression in parallel with CRY2. The corrected threshold for significance corresponded to a nominal P-value of 0.001 and was calculated applying a Bonferroni correction considering the partial LD between markers [60], [61]. As a novel finding we report here that CRY2 expression is reduced and non-responsive to sleep deprivation in depressed patients with bipolar disorder, suggesting that CRY2 plays a role in bipolar disorder.

There are some limitations in our work. A limitation in the CRY2 expression analysis was the few time points for the cases, since they do not allow analysis of the CRY2 rhythm and possible phase shift among the cases. However, irrespectively of whether there was a phase shift or not among the cases, CRY2 mRNA levels were significantly lower in cases compared to controls after sleep deprivation, which suggests that CRY2 levels is a trait marker that is non-responsive to the antidepressant sleep deprivation. Though, we cannot exclude that the CRY2 mRNA levels in the patients were influenced by medication. All the patients were on an SSRI together with a mood stabilizer. Lithium is known to affect the circadian clock through inhibition of GSK3 [2], and there are reports of SSRI influencing the circadian rhythms [62]. Since all controls and 10 of the patients were studied during the summer, the season of the examinations does not likely explain the CRY2 mRNA level difference between patients and controls. A limitation in the genetic analysis was the number of patients. However, the cases and the controls in both the Swedish and the Finnish materials represented ethnically homogeneous populations, and the controls were extensively scored to exclude any mental illness, reducing bias due to the ethnic variation and dilution of genetic effects due to disorder heterogeneity among the controls.

In conclusion, we show data suggesting that variation in CRY2 links to depression. CRY2 mRNA expression levels were lowered in PBMCs from depressed patients with bipolar disorder. In two ethnically homogeneous population-based samples, CRY2 SNPs were associated with winter depression. Although this study contains a replication of the genetic findings, it remains limited due to the lack of replication of the expression findings, thus warranting further studies.

Methods

CRY2 Expression Study

Informed consent was obtained from each participant using an approved University of California Institutional Review Board (IRB) protocol.

Thirteen patients (10 men, 3 women, out of which 11 with European, one with Asian and one with African decent) with bipolar type 1 disorder according to DSM-IV criteria, aged 40.2 years on average (SD = 13.4, range = 18–57), were enrolled in a sleep deprivation study at University of California Irvine (UCI) and San Diego (UCSD). During the study, all the patients were depressed, with the score on the 21-item Hamilton Depression Rating Scale [63] being 17.3 (SD = 6.4) on average. Eight healthy volunteers (4 women, 4 men, all with European decent), aged 23.6 years on average (SD = 5.9, range = 19–34 years), served as controls. Both the patients and the controls were hospitalized at the sleep research center at University of California at Irvine Medical Center (UCIMC) for 48 hours and deprived of sleep for 21 hours after an overnight stay (Figure 2). The season for the study of the patients was summer (n = 10) and winter (n = 3), and all controls were studied in the summer. All the patients were on an SSRI and either lithium (n = 8), valproate (n = 1) or lamotrigine (n = 4). None of the controls was on medication.

For the controls, venous blood samples were drawn at 9 different times, beginning at 7 p.m. and every 6 hours thereafter. For the patients, the blood was drawn at 1 p.m. on the days before and after sleep deprivation, corresponding to the time points #4 and #8 for the control samples. The patient blood was drawn in standard Vacutainer tubes without additive (10 ml, Becton Dickinson, Franklin Lakes, NJ, USA). From each control the blood was collected into both a standard acid citrate dextrose (ACD) tube (Becton Dickinson, Franklin Lakes, NJ, USA) as well as a Tempus Blood RNA Tube (Applied Biosystems, Foster City, CA, USA) in the same blood draw, yielding the total of 25–35 ml of whole blood. To check for potential difference in RNA degradation between control and bipolar disorder samples, the level of the house-keeping transcript GAPDH (glyceraldehyde-3-phosphate dehydrogenase) was analyzed with qPCR (data in the results section). To further rule out possible bias in comparison between the patients and controls due to type of collection tube, we compared CRY2 and GAPDH mRNA levels using qPCR between the Tempus tube and the ACD tube from the controls collected at the same blood draw. The CRY2 and GAPDH mRNA levels were not different between tube types (data not shown).

Peripheral Blood Mononuclear Cells (PBMC) Isolation and RNA Extraction

All samples from both the patients and the controls were prepared similarly. Within a few minutes in room temperature after blood draw, the whole blood samples were layered onto Ficoll (Amersham Biosciences, Piscataway, NJ, USA), and PBMC were separated by density gradient centrifugation at 2500 rpm at room temperature for 20 min. The resulting ‘buffy’ coat at the interface was added to 5 ml phosphate buffered saline (PBS) at pH of 7.4 (Invitrogen, Carlsbad, CA, USA) and cell counts were taken. Cells were centrifuged at 1000 rpm for 10 min at room temperature. The resulting pellets were resuspended in 1 ml Trizol and stored at minus 80°C. Total RNA was extracted using the standard Trizol isolation protocol (Invitrogen, Carlsbad, CA, USA). The RNA was resuspended in 50 µL DEPC water, cleaned by passing over silica-based mini-spin columns (Qiagen RNeasy PlusMini Kit, Valencia, CA, USA) and analyzed for quality and quantity on a 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA) and concentration was adjusted to 1 µg/µl.

Real Time Quantitative PCR (qPCR)

The total RNA (1 µg RNA) from both the patients and the controls were identically synthesized into complementary DNA (cDNA) using Oligo d(T)16 primer and TaqMan Reverse Transcription Reagents and (Applied Biosystems, Foster City, CA, USA) at 25°C for 10 min, 48°C for 30 min and 95°C for 5 min. The quantitative PCR (qPCR) was performed on an ABI 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) with triplicates for each cDNA, using Robot Biomek3000 (Beckman Coulter, Fullerton, CA, USA), with all samples run in the same plates with same normalization procedures. The reaction was performed in 12.5 µl consisting of 6.25 µl 2xSYBR Green Master Mix (Applied Biosystems, Foster City, CA, USA); 2.5 pmol of each primer; 2 µl 1∶10 dilution of cDNA template (corresponding approximately to 4 ng RNA). The thermal cycling profiles were 50°C for 2 min (incubation), 95°C for 10 min (activation), 45 cycles at 95°C for 15 sec (denaturation) and 60°C for 1 min (annealing/extension), then by dissociation step as 95°C for 15 sec, 60°C for 15 sec, and 95°C for 15 sec. The CRY2 primer sequences spanned exon 6–7 junction within target sequence (Affymetrix probeset 3329058, transcript ID 3329029) and were: forward 5′-CCTACCTGCGCTTTGGTTGT-3′, reverse 5′-TGCTGTTCCGCTTCACCTTT-3′. The primers were pre-tested using brain cDNA, genomic DNA, no temple control (NTC) and RT minus on the ABI 7900HT. The results ensured good amplification of cDNA, and that neither residual genomic DNA nor primer-dimer contributed to measurements.

Genetic Association Study

The local ethic committees (National Public Health Institute, Karolinska Institutet, and University of Umeå) approved the study protocol, and all the participants signed an informed consent after the protocol had been fully explained.

Patients were recruited from outpatient services at which senior specialists in psychiatry, familiar with our research program, were working. All patients met the DSM-IV diagnostic criteria for major depressive disorder with the seasonal (winter) pattern [64]. The consensus diagnosis of two independent psychiatrists was required for inclusion for each patient. Patients were from Sweden (118 patients: 13.6% men, 86.4% women) and Finland (86 patients: 29.1% men, 70.9% women), both being of local origin.

Controls matched for ethnicity and nationality (1011 Swedes: 39.2% men, 60.8% women; and 1096 Finns: 29.0% men, 71.0% women) had no current or past psychotic or mood disorder as assessed with interviews or self-rated questionnaires and were representative of the Swedish and Finnish adult populations. These population-based studies have been described more in detail earlier [65]–[67]. The Swedish cases and controls lived in Västerbotten and Stockholm areas. The Swedish population has no strong internal genetic borders [24]. The Finnish cases lived in the capital area (Helsinki and its surroundings), with inhabitants from all the regions in Finland, and the controls were collected from the nation-wide health examination study.

Single-Nucleotide Polymorphism (SNP) Selection and Genotyping

Haplotype-tagging SNPs covering the variation in CRY2 were selected using the HapMap database [68], applying the cut-off values of 0.8 for r2 and of 0.1 for the minor allele frequency (MAF). The genotyping was performed using the SEQUENOM iPLEX Application with the MassARRAY System (Sequenom, Inc., San Diego, CA, USA). To control for quality, 2.5% of the samples were genotyped in duplicates. The genotypes of these replicas were all in agreement.

The four CRY2 SNPs genotyped fulfilled the criteria for SNP association test including a success rate of at least 90% of the reactions, and the Hardy-Weinberg equilibrium (HWE) with P>0.05 among the controls. The DNA samples from 10 Swedish cases and 53 controls, and 5 Finnish cases and 24 controls were excluded from analysis, since genotyping of them failed in more than 20% of a separate set of 111 SNP assays.

Statistical Analyses

For the RNA analysis, the cycle threshold (Ct) was determined approximately in the middle of exponential phase of the amplification, and generally default condition on the software SDS 2.3 was used. The average Ct value was accepted if the CV was lower than 2% and SD was lower than 0.39. The Ct values for CRY2 and the house-keeping gene GAPDH were normalized using the housekeeping genes C1ORF82 and TFG. The selection of C1ORF82 and TFG and the normalization method was based on the GENORM [69]. Within each of the case and the control groups, correlation for CRY2 and GAPDH expression levels to age was tested and showed no significant correlation. Among controls, the CRY2 and GAPDH expression levels were compared between pre-sleep deprivation time points (pooled time points from #1 to #4) and sleep deprived time points (pooled time points from #6 to #9), where the individual was used as a factor in a mixed model ANOVA.

CRY2 and GAPDH expression levels comparison between bipolar patients (time point 1 p.m.) and controls (pooled time points from #1 to #4, and from #6 to #9) was done using ANOVA with the covariates gender, age and sleep deprivation (yes/no). Significant effect in ANOVA (p<0.05) was followed-up with post hoc test.

The SNPs were analyzed for allele frequency differences between the cases and controls in each population sample using logistic regression. Gender was used as a covariate, since women were overrepresented among Swedish cases compared with controls. To obtain empirical significance values, permutation tests with up to 10,000 permutations were performed. There was a similarity between the gender corrected effect size and the effect size when no gender was corrected for, especially in the Finnish sample (Table 1). Hence, the gender was not generally controlled for in the subsequent analyses. However, suggestive genotypic association findings in the Swedish sample were verified using logistic regression with the gender as a covariate, and displayed a very minor gender effect. The power was ≥0.8 to detect allele frequency differences at odds ratio (OR) of ≥1.8 in the Swedish sample and at OR of ≥2.0 in the Finnish sample when the MAF was of 0.1 to 0.3 and the non-corrected P was of <0.05 (http://pngu.mgh.harvard.edu/~purcell/gpc/cc2.html). Those SNPs whose allele frequency difference was indicative (P<0.05) were analyzed for genotype association. Calculations were performed using the PLINK program, version 1.04 [70].

The linkage disequilibrium (LD) measure D' was calculated between the SNPs using the Haploview program, version 3.2 [71], which applies an approach similar to the partition-ligation-expectation-maximization algorithm [72]. Haplotype blocks were constructed using the LD block parameters [73], and the D' confidence interval algorithm in the Haploview program. Tests for haplotype frequency difference between the cases and controls were calculated for the haplotype blocks harboring the SNPs having P of <0.05 for allele association using the Haploview program. The SNP rs10838524 in CRY2 was included in haplotype analysis even though it was not included in any block but was nearby and associated as a single SNP with winter depression.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

All the detected haplotypes of the CRY2 block.

https://doi.org/10.1371/journal.pone.0009407.s001

(0.04 MB DOC)

Author Contributions

Conceived and designed the experiments: CL TP MS. Performed the experiments: LS PS MV RA. Analyzed the data: CL LS PS TP MV WEB YF JW JK TP MS. Contributed reagents/materials/analysis tools: CL TP MV WEB RA YF JW JK TP MS. Wrote the paper: CL LS PS TP MV WEB RA YF JW JK TP MS.

Citation: Lavebratt C, Sjöholm LK, Soronen P, Paunio T, Vawter MP, Bunney WE, et al. (2010) CRY2 Is Associated with Depression. PLoS ONE 5(2): e9407. https://doi.org/10.1371/journal.pone.0009407

References

1. Germain A, Kupfer DJ (2008) Circadian rhythm disturbances in depression. Hum Psychopharmacol 23: 571–585.A. GermainDJ Kupfer2008Circadian rhythm disturbances in depression.Hum Psychopharmacol23571585

2. Yin L, Wang J, Klein PS, Lazar MA (2006) Nuclear receptor Rev-erbalpha is a critical lithium-sensitive component of the circadian clock. Science 311: 1002–1005.L. YinJ. WangPS KleinMA Lazar2006Nuclear receptor Rev-erbalpha is a critical lithium-sensitive component of the circadian clock.Science31110021005

3. Faedda GL, Tondo L, Teicher MH, Baldessarini RJ, Gelbard HA, et al. (1993) Seasonal mood disorders: patterns of seasonal recurrence in mania and depression. Arch Gen Psychiatry 50: 17–23.GL FaeddaL. TondoMH TeicherRJ BaldessariniHA Gelbard1993Seasonal mood disorders: patterns of seasonal recurrence in mania and depression.Arch Gen Psychiatry501723

4. Partonen T, Lönnqvist J (1998) Seasonal affective disorder. Lancet 352: 1369–1374.T. PartonenJ. Lönnqvist1998Seasonal affective disorder.Lancet35213691374

5. Rosenthal NE, Sack DA, Gillin JC, Lewy AJ, Goodwin FK, et al. (1984) Seasonal affective disorder: a description of the syndrome and preliminary findings with light therapy. Arch Gen Psychiatry 41: 72–80.NE RosenthalDA SackJC GillinAJ LewyFK Goodwin1984Seasonal affective disorder: a description of the syndrome and preliminary findings with light therapy.Arch Gen Psychiatry417280

6. Kieseppä T, Partonen T, Haukka J, Kaprio J, Lönnqvist J (2004) High concordance of bipolar I disorder in a nationwide sample of twins. Am J Psychiatry 161: 1814–1821.T. KieseppäT. PartonenJ. HaukkaJ. KaprioJ. Lönnqvist2004High concordance of bipolar I disorder in a nationwide sample of twins.Am J Psychiatry16118141821

7. Bunney WE, Bunney BG (2000) Molecular clock genes in man and lower animals: possible implications for circadian abnormalities in depression. Neuropsychopharmacology 22: 335–345.WE BunneyBG Bunney2000Molecular clock genes in man and lower animals: possible implications for circadian abnormalities in depression.Neuropsychopharmacology22335345

8. McClung CA (2007) Circadian genes, rhythms and the biology of mood disorders. Pharmacol Ther 114: 222–232.CA McClung2007Circadian genes, rhythms and the biology of mood disorders.Pharmacol Ther114222232

9. Lewy AJ, Sack RL, Miller LS, Hoban TM (1987) Antidepressant and circadian phase-shifting effects of light. Science 235: 352–354.AJ LewyRL SackLS MillerTM Hoban1987Antidepressant and circadian phase-shifting effects of light.Science235352354

10. Souêtre E, Salvati E, Belugou JL, Pringuey D, Candito M, et al. (1989) Circadian rhythms in depression and recovery: evidence for blunted amplitude as the main chronobiological abnormality. Psychiatry Res 28: 263–278.E. SouêtreE. SalvatiJL BelugouD. PringueyM. Candito1989Circadian rhythms in depression and recovery: evidence for blunted amplitude as the main chronobiological abnormality.Psychiatry Res28263278

11. Teicher MH, Glod CA, Magnus E, Harper D, Benson G, et al. (1997) Circadian rest-activity disturbances in seasonal affective disorder. Arch Gen Psychiatry 54: 124–130.MH TeicherCA GlodE. MagnusD. HarperG. Benson1997Circadian rest-activity disturbances in seasonal affective disorder.Arch Gen Psychiatry54124130

12. Glod CA, Teicher MH, Polcari A, McGreenery CE, Ito Y (1997) Circadian rest-activity disturbances in children with seasonal affective disorder. J Am Acad Child Adolesc Psychiatry 36: 188–195.CA GlodMH TeicherA. PolcariCE McGreeneryY. Ito1997Circadian rest-activity disturbances in children with seasonal affective disorder.J Am Acad Child Adolesc Psychiatry36188195

13. Czeisler CA, Kronauer RE, Mooney JJ, Anderson JL, Allan JS (1987) Biologic rhythm disorders, depression, and phototherapy: a new hypothesis. Psychiatr Clin North Am 10: 687–709.CA CzeislerRE KronauerJJ MooneyJL AndersonJS Allan1987Biologic rhythm disorders, depression, and phototherapy: a new hypothesis.Psychiatr Clin North Am10687709

14. Thompson C, Stinson D, Smith A (1990) Seasonal affective disorder and season-dependent abnormalities of melatonin suppression by light. Lancet 336: 703–706.C. ThompsonD. StinsonA. Smith1990Seasonal affective disorder and season-dependent abnormalities of melatonin suppression by light.Lancet336703706

15. Srinivasan V, Pandi-Perumal SR, Trakht I, Spence DW, Hardeland R (2009) Pathophysiology of depression: role of sleep and the melatonergic system. Psychiatry Res 165: 201–214.V. SrinivasanSR Pandi-PerumalI. TrakhtDW SpenceR. Hardeland2009Pathophysiology of depression: role of sleep and the melatonergic system.Psychiatry Res165201214

16. Lamont EW, Legault-Coutu D, Cermakian N, Boivin DB (2007) The role of circadian clock genes in mental disorders. Dialogues Clin Neurosci 9: 333–342.EW LamontD. Legault-CoutuN. CermakianDB Boivin2007The role of circadian clock genes in mental disorders.Dialogues Clin Neurosci9333342

17. Hsu DS, Zhao X, Zhao S, Kazantsev A, Wang RP, et al. (1996) Putative human blue-light photoreceptors hCRY1 and hCRY2 are flavoproteins. Biochemistry 35: 13871–13877.DS HsuX. ZhaoS. ZhaoA. KazantsevRP Wang1996Putative human blue-light photoreceptors hCRY1 and hCRY2 are flavoproteins.Biochemistry351387113877

18. Kume K, Zylka MJ, Sriram S, Shearman LP, Weaver DR, et al. (1999) mCRY1 and mCRY2 are essential components of the negative limb of the circadian clock feedback loop. Cell 98: 193–205.K. KumeMJ ZylkaS. SriramLP ShearmanDR Weaver1999mCRY1 and mCRY2 are essential components of the negative limb of the circadian clock feedback loop.Cell98193205

19. Miyamoto Y, Sancar A (1998) Vitamin B2-based blue-light photoreceptors in the retinohypothalamic tract as the photoactive pigments for setting the circadian clock in mammals. Proc Natl Acad Sci U S A 95: 6097–6102.Y. MiyamotoA. Sancar1998Vitamin B2-based blue-light photoreceptors in the retinohypothalamic tract as the photoactive pigments for setting the circadian clock in mammals.Proc Natl Acad Sci U S A9560976102

20. Hoang N, Schleicher E, Kacprzak S, Bouly JP, Picot M, et al. (2008) Human and Drosophila cryptochromes are light activated by flavin photoreduction in living cells. PLoS Biol 6: e160.N. HoangE. SchleicherS. KacprzakJP BoulyM. Picot2008Human and Drosophila cryptochromes are light activated by flavin photoreduction in living cells.PLoS Biol6e160

21. Cameron MA, Barnard AR, Hut RA, Bonnefont X, van der Horst GT, et al. (2008) Electroretinography of wild-type and Cry mutant mice reveals circadian tuning of photopic and mesopic retinal responses. J Biol Rhythms 23: 489–501.MA CameronAR BarnardRA HutX. BonnefontGT van der Horst2008Electroretinography of wild-type and Cry mutant mice reveals circadian tuning of photopic and mesopic retinal responses.J Biol Rhythms23489501

22. Buijs RM, Wortel J, Van Heerikhuize JJ, Feenstra MG, Ter Horst GJ, et al. (1999) Anatomical and functional demonstration of a multisynaptic suprachiasmatic nucleus adrenal (cortex) pathway. Eur J Neurosci 11: 1535–1544.RM BuijsJ. WortelJJ Van HeerikhuizeMG FeenstraGJ Ter Horst1999Anatomical and functional demonstration of a multisynaptic suprachiasmatic nucleus adrenal (cortex) pathway.Eur J Neurosci1115351544

23. de la Iglesia HO, Schwartz WJ (2002) A subpopulation of efferent neurons in the mouse suprachiasmatic nucleus is also light responsive. Neuroreport 13: 857–860.HO de la IglesiaWJ Schwartz2002A subpopulation of efferent neurons in the mouse suprachiasmatic nucleus is also light responsive.Neuroreport13857860

24. Luo AH, Aston-Jones G (2009) Circuit projection from suprachiasmatic nucleus to ventral tegmental area: a novel circadian output pathway. Eur J Neurosci 29: 748–760.AH LuoG. Aston-Jones2009Circuit projection from suprachiasmatic nucleus to ventral tegmental area: a novel circadian output pathway.Eur J Neurosci29748760

25. Dardente H, Fortier EE, Martineau V, Cermakian N (2007) Cryptochromes impair phosphorylation of transcriptional activators in the clock: a general mechanism for circadian repression. Biochem J 402: 525–536.H. DardenteEE FortierV. MartineauN. Cermakian2007Cryptochromes impair phosphorylation of transcriptional activators in the clock: a general mechanism for circadian repression.Biochem J402525536

26. Lappalainen T, Hannelius U, Salmela E, von Döbeln U, Lindgren CM, et al. (2009) Population structure in contemporary Sweden – a Y-chromosomal and mitochondrial DNA analysis. Ann Hum Genet 73: 61–73.T. LappalainenU. HanneliusE. SalmelaU. von DöbelnCM Lindgren2009Population structure in contemporary Sweden – a Y-chromosomal and mitochondrial DNA analysis.Ann Hum Genet736173

27. Miyazaki K, Nagase T, Mesaki M, Narukawa J, Ohara O, et al. (2004) Phosphorylation of clock protein PER1 regulates its circadian degradation in normal human fibroblasts. Biochem J 380: 95–103.K. MiyazakiT. NagaseM. MesakiJ. NarukawaO. Ohara2004Phosphorylation of clock protein PER1 regulates its circadian degradation in normal human fibroblasts.Biochem J38095103

28. Tsinkalovsky O, Smaaland R, Rosenlund B, Sothern RB, Hirt A, et al. (2007) Circadian variations in clock gene expression of human bone marrow CD34+ cells. J Biol Rhythms 22: 140–150.O. TsinkalovskyR. SmaalandB. RosenlundRB SothernA. Hirt2007Circadian variations in clock gene expression of human bone marrow CD34+ cells.J Biol Rhythms22140150

29. Ozgur S, Sancar A (2003) Purification and properties of human blue-light photoreceptor cryptochrome 2. Biochemistry 42: 2926–2932.S. OzgurA. Sancar2003Purification and properties of human blue-light photoreceptor cryptochrome 2.Biochemistry4229262932

30. Ukai-Tadenuma M, Kasukawa T, Ueda HR (2008) Proof-by-synthesis of the transcriptional logic of mammalian circadian clocks. Nat Cell Biol 10: 1154–1163.M. Ukai-TadenumaT. KasukawaHR Ueda2008Proof-by-synthesis of the transcriptional logic of mammalian circadian clocks.Nat Cell Biol1011541163

31. Mirsky HP, Liu AC, Welsh DK, Kay SA, Doyle FJ 3rd (2009) A model of the cell-autonomous mammalian circadian clock. Proc Natl Acad Sci U S A 106: 11107–11112.HP MirskyAC LiuDK WelshSA KayFJ Doyle 3rd2009A model of the cell-autonomous mammalian circadian clock.Proc Natl Acad Sci U S A1061110711112

32. Zhang EE, Liu AC, Hirota T, Miraglia LJ, Welch G, et al. (2009) A genome-wide RNAi screen for modifiers of the circadian clock in human cells. Cell 139: 199–210.EE ZhangAC LiuT. HirotaLJ MiragliaG. Welch2009A genome-wide RNAi screen for modifiers of the circadian clock in human cells.Cell139199210

33. Thresher RJ, Vitaterna MH, Miyamoto Y, Kazantsev A, Hsu DS, et al. (1998) Role of mouse cryptochrome blue-light photoreceptor in circadian photoresponses. Science 282: 1490–1494.RJ ThresherMH VitaternaY. MiyamotoA. KazantsevDS Hsu1998Role of mouse cryptochrome blue-light photoreceptor in circadian photoresponses.Science28214901494

34. van der Horst GT, Muijtjens M, Kobayashi K, Takano R, Kanno S, et al. (1999) Mammalian Cry1 and Cry2 are essential for maintenance of circadian rhythms. Nature 398: 627–630.GT van der HorstM. MuijtjensK. KobayashiR. TakanoS. Kanno1999Mammalian Cry1 and Cry2 are essential for maintenance of circadian rhythms.Nature398627630

35. Vitaterna MH, Selby CP, Todo T, Niwa H, Thompson C, et al. (1999) Differential regulation of mammalian period genes and circadian rhythmicity by cryptochromes 1 and 2. Proc Natl Acad Sci U S A 96: 12114–12119.MH VitaternaCP SelbyT. TodoH. NiwaC. Thompson1999Differential regulation of mammalian period genes and circadian rhythmicity by cryptochromes 1 and 2.Proc Natl Acad Sci U S A961211412119

36. Spoelstra K, Daan S (2008) Effects of constant light on circadian rhythmicity in mice lacking functional cry genes: dissimilar from per mutants. J Comp Physiol A 194: 235–242.K. SpoelstraS. Daan2008Effects of constant light on circadian rhythmicity in mice lacking functional cry genes: dissimilar from per mutants.J Comp Physiol A194235242

37. Thompson CL, Sancar A (2002) Photolyase/cryptochrome blue-light photoreceptors use photon energy to repair DNA and reset the circadian clock. Oncogene 21: 9043–9056.CL ThompsonA. Sancar2002Photolyase/cryptochrome blue-light photoreceptors use photon energy to repair DNA and reset the circadian clock.Oncogene2190439056

38. Liu AC, Welsh DK, Ko CH, Tran HG, Zhang EE, et al. (2007) Intercellular coupling confers robustness against mutations in the SCN circadian clock network. Cell 129: 605–616.AC LiuDK WelshCH KoHG TranEE Zhang2007Intercellular coupling confers robustness against mutations in the SCN circadian clock network.Cell129605616

39. Oster H, Yasui A, van der Horst GT, Albrecht U (2002) Disruption of mCry2 restores circadian rhythmicity in mPer2 mutant mice. Genes Dev 16: 2633–2638.H. OsterA. YasuiGT van der HorstU. Albrecht2002Disruption of mCry2 restores circadian rhythmicity in mPer2 mutant mice.Genes Dev1626332638

40. Lincoln GA, Andersson H, Hazlerigg D (2003) Clock genes and the long-term regulation of prolactin secretion: evidence for a photoperiod/circannual timer in the pars tuberalis. J Neuroendocrinol 15: 390–397.GA LincolnH. AnderssonD. Hazlerigg2003Clock genes and the long-term regulation of prolactin secretion: evidence for a photoperiod/circannual timer in the pars tuberalis.J Neuroendocrinol15390397

41. Koorengevel KM, Beersma DG, den Boer JA, van den Hoofdakker RH (2002) A forced desynchrony study of circadian pacemaker characteristics in seasonal affective disorder. J Biol Rhythms 17: 463–475.KM KoorengevelDG BeersmaJA den BoerRH van den Hoofdakker2002A forced desynchrony study of circadian pacemaker characteristics in seasonal affective disorder.J Biol Rhythms17463475

42. Elsass P, Mellerup ET, Rafaelsen OJ, Theilgaard A (1979) Lithium effects on time estimation and mood in manic-melancholic patients: a study of diurnal variations. Acta Psychiatr Scand 60: 263–271.P. ElsassET MellerupOJ RafaelsenA. Theilgaard1979Lithium effects on time estimation and mood in manic-melancholic patients: a study of diurnal variations.Acta Psychiatr Scand60263271

43. Daan S, Albrecht U, van der Horst GT, Illnerová H, Roenneberg T, et al. (2001) Assembling a clock for all seasons: are there M and E oscillators in the genes? J Biol Rhythms 16: 105–116.S. DaanU. AlbrechtGT van der HorstH. IllnerováT. Roenneberg2001Assembling a clock for all seasons: are there M and E oscillators in the genes?J Biol Rhythms16105116

44. Thompson C, Childs PA, Martin NJ, Rodin I, Smythe PJ (1997) Effects of morning phototherapy on circadian markers in seasonal affective disorder. Br J Psychiatry 170: 431–435.C. ThompsonPA ChildsNJ MartinI. RodinPJ Smythe1997Effects of morning phototherapy on circadian markers in seasonal affective disorder.Br J Psychiatry170431435

45. Nurnberger JI Jr, Adkins S, Lahiri DK, Mayeda A, Hu K, et al. (2000) Melatonin suppression by light in euthymic bipolar and unipolar patients. Arch Gen Psychiatry 57: 572–579.JI Nurnberger JrS. AdkinsDK LahiriA. MayedaK. Hu2000Melatonin suppression by light in euthymic bipolar and unipolar patients.Arch Gen Psychiatry57572579

46. Wehr TA, Duncan WC Jr, Sher L, Aeschbach D, Schwartz PJ, et al. (2001) A circadian signal of change of season in patients with seasonal affective disorder. Arch Gen Psychiatry 58: 1108–1114.TA WehrWC Duncan JrL. SherD. AeschbachPJ Schwartz2001A circadian signal of change of season in patients with seasonal affective disorder.Arch Gen Psychiatry5811081114

47. Partonen T, Treutlein J, Alpman A, Frank J, Johansson C, et al. (2007) Three circadian clock genes Per2, Arntl, and Npas2 contribute to winter depression. Ann Med 39: 229–238.T. PartonenJ. TreutleinA. AlpmanJ. FrankC. Johansson2007Three circadian clock genes Per2, Arntl, and Npas2 contribute to winter depression.Ann Med39229238

48. Lavebratt C, Sjöholm LK, Partonen T, Schalling M, Forsell Y (2009) PER2 variantion is associated with depression vulnerability. Am J Med Genet B Neuropsychiatr Genet Aug 19: C. LavebrattLK SjöholmT. PartonenM. SchallingY. Forsell2009PER2 variantion is associated with depression vulnerability.Am J Med Genet B Neuropsychiatr Genet Aug19

49. Lewy AJ, Lefler BJ, Emens JS, Bauer VK (2006) The circadian basis of winter depression. Proc Natl Acad Sci U S A 103: 7414–7419.AJ LewyBJ LeflerJS EmensVK Bauer2006The circadian basis of winter depression.Proc Natl Acad Sci U S A10374147419

50. Kanazawa T, Chana G, Glatt SJ, Mizuno H, Masliah E, et al. (2008) The utility of SELENBP1 gene expression as a biomarker for major psychotic disorders: replication in schizophrenia and extension to bipolar disorder with psychosis. Am J Med Genet B Neuropsychiatr Genet 147B: 686–689.T. KanazawaG. ChanaSJ GlattH. MizunoE. Masliah2008The utility of SELENBP1 gene expression as a biomarker for major psychotic disorders: replication in schizophrenia and extension to bipolar disorder with psychosis.Am J Med Genet B Neuropsychiatr Genet147B686689

51. Le-Niculescu H, Patel SD, Bhat M, Kuczenski R, Faraone SV, et al. (2009) Convergent functional genomics of genome-wide association data for bipolar disorder: comprehensive identification of candidate genes, pathways and mechanisms. Am J Med Genet B Neuropsychiatr Genet 150B: 155–181.H. Le-NiculescuSD PatelM. BhatR. KuczenskiSV Faraone2009Convergent functional genomics of genome-wide association data for bipolar disorder: comprehensive identification of candidate genes, pathways and mechanisms.Am J Med Genet B Neuropsychiatr Genet150B155181

52. Niculescu AB 3rd, Segal DS, Kuczenski R, Barrett T, Hauger RL, et al. (2000) Identifying a series of candidate genes for mania and psychosis: a convergent functional genomics approach. Physiol Genomics 4: 83–91.AB Niculescu 3rdDS SegalR. KuczenskiT. BarrettRL Hauger2000Identifying a series of candidate genes for mania and psychosis: a convergent functional genomics approach.Physiol Genomics48391

53. Ogden CA, Rich ME, Schork NJ, Paulus MP, Geyer MA, et al. (2004) Candidate genes, pathways and mechanisms for bipolar (manic-depressive) and related disorders: an expanded convergent functional genomics approach. Mol Psychiatry 9: 1007–1029.CA OgdenME RichNJ SchorkMP PaulusMA Geyer2004Candidate genes, pathways and mechanisms for bipolar (manic-depressive) and related disorders: an expanded convergent functional genomics approach.Mol Psychiatry910071029

54. Shi J, Wittke-Thompson JK, Badner JA, Hattori E, Potash JB, et al. (2008) Clock genes may influence bipolar disorder susceptibility and dysfunctional circadian rhythm. Am J Med Genet Part B Neuropsychiatr Genet 147B: 1047–1055.J. ShiJK Wittke-ThompsonJA BadnerE. HattoriJB Potash2008Clock genes may influence bipolar disorder susceptibility and dysfunctional circadian rhythm.Am J Med Genet Part B Neuropsychiatr Genet147B10471055

55. Mansour HA, Wood J, Logue T, Chowdari KV, Dayal M, et al. (2006) Association study of eight circadian genes with bipolar I disorder, schizoaffective disorder and schizophrenia. Genes Brain Behav 5: 150–157.HA MansourJ. WoodT. LogueKV ChowdariM. Dayal2006Association study of eight circadian genes with bipolar I disorder, schizoaffective disorder and schizophrenia.Genes Brain Behav5150157

56. Nievergelt CM, Kripke DF, Barrett TB, Burg E, Remick RA, et al. (2006) Suggestive evidence for association of the circadian genes PERIOD3 and ARNTL with bipolar disorder. Am J Med Genet. B Neuropsychiatr Genet 141: 234–241.CM NievergeltDF KripkeTB BarrettE. BurgRA Remick2006Suggestive evidence for association of the circadian genes PERIOD3 and ARNTL with bipolar disorder. Am J Med Genet.B Neuropsychiatr Genet141234241

57. Kripke DF, Nievergelt CM, Joo EJ, Shekhtman T, Kelsoe JR (2009) Circadian polymorphisms associated with affective disorders. J Circadian Rhythm 7: 2.DF KripkeCM NievergeltEJ JooT. ShekhtmanJR Kelsoe2009Circadian polymorphisms associated with affective disorders.J Circadian Rhythm72

58. Johansson C, Willeit M, Smedh C, Ekholm J, Paunio T, et al. (2003) Circadian clock-related polymorphisms in seasonal affective disorder and their relevance to diurnal preference. Neuropsychopharmacology 28: 734–739.C. JohanssonM. WilleitC. SmedhJ. EkholmT. Paunio2003Circadian clock-related polymorphisms in seasonal affective disorder and their relevance to diurnal preference.Neuropsychopharmacology28734739

59. Hampp G, Ripperger JA, Houben T, Schmutz I, Blex C, et al. (2008) Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood. Curr Biol 18: 678–683.G. HamppJA RippergerT. HoubenI. SchmutzC. Blex2008Regulation of monoamine oxidase A by circadian-clock components implies clock influence on mood.Curr Biol18678683

60. Nyholt DR (2004) A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet 74: 765–769.DR Nyholt2004A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other.Am J Hum Genet74765769

61. Gao X, Starmer J, Martin ER (2008) A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol 32: 361–369.X. GaoJ. StarmerER Martin2008A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms.Genet Epidemiol32361369

62. Sprouse J, Braselton J, Reynolds L (2006) Fluoxetine modulates the circadian biological clock via phase advances of suprachiasmatic nucleus neuronal firing. Biol Psychiatry 60: 896–899.J. SprouseJ. BraseltonL. Reynolds2006Fluoxetine modulates the circadian biological clock via phase advances of suprachiasmatic nucleus neuronal firing.Biol Psychiatry60896899

63. Hamilton M (1960) A rating scale for depression. J Neurol Neurosurg Psychiatry 23: 56–62.M. Hamilton1960A rating scale for depression.J Neurol Neurosurg Psychiatry235662

64. American Psychiatric Association (1994) Diagnostic and Statistical Manual of Mental Disorders, 4th ed. Washington, DC: American Psychiatric Press. 886 p.American Psychiatric Association1994Diagnostic and Statistical Manual of Mental Disorders, 4th ed.Washington, DCAmerican Psychiatric Press886

65. Lundberg I, Damström Thakker K, Hällström T, Forsell Y (2005) Determinants of non-participation, and the effects of non-participation on potential cause-effect relationships, in the PART study on mental disorders. Soc Psychiatry Psychiatr Epidemiol 40: 475–483.I. LundbergK. Damström ThakkerT. HällströmY. Forsell2005Determinants of non-participation, and the effects of non-participation on potential cause-effect relationships, in the PART study on mental disorders.Soc Psychiatry Psychiatr Epidemiol40475483

66. Pirkola SP, Isometsä E, Suvisaari J, Aro H, Joukamaa M, et al. (2005) DSM-IV mood-, anxiety- and alcohol use disorders and their comorbidity in the Finnish general population–results from the Health 2000 Study. Soc Psychiatry Psychiatr Epidemiol 40: 1–10.SP PirkolaE. IsometsäJ. SuvisaariH. AroM. Joukamaa2005DSM-IV mood-, anxiety- and alcohol use disorders and their comorbidity in the Finnish general population–results from the Health 2000 Study.Soc Psychiatry Psychiatr Epidemiol40110

67. Perälä J, Suvisaari J, Saarni SI, Kuoppasalmi K, Isometsä E, et al. (2007) Lifetime prevalence of psychotic and bipolar I disorders in a general population. Arch Gen Psychiatry 64: 19–28.J. PeräläJ. SuvisaariSI SaarniK. KuoppasalmiE. Isometsä2007Lifetime prevalence of psychotic and bipolar I disorders in a general population.Arch Gen Psychiatry641928

68. The International HapMap Consortium (2005) A haplotype map of the human genome. Nature 437: 1299–1320.The International HapMap Consortium2005A haplotype map of the human genome.Nature43712991320

69. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: research0034.1–0034.11.J. VandesompeleK. De PreterF. PattynB. PoppeN. Van Roy2002Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.Genome Biol3research0034.10034.11

70. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

71. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.JC BarrettB. FryJ. MallerMJ Daly2005Haploview: analysis and visualization of LD and haplotype maps.Bioinformatics21263265

72. Qin ZS, Niu T, Liu JS (2002) Partition-ligation-expectation-maximization algorithm for haplotype inference with single-nucleotide polymorphisms. Am J Hum Genet 71: 1242–1247.ZS QinT. NiuJS Liu2002Partition-ligation-expectation-maximization algorithm for haplotype inference with single-nucleotide polymorphisms.Am J Hum Genet7112421247

73. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229.SB GabrielSF SchaffnerH. NguyenJM MooreJ. Roy2002The structure of haplotype blocks in the human genome.Science29622252229