J.H. Smit

Genome-wide association for major depressive disorder: a possible role for the presynaptic protein piccolo

Major depressive disorder (MDD) is a common complex trait with enormous public health significance. As part of the Genetic Association Information Network initiative of the US Foundation for the National Institutes of Health, we conducted a genome-wide association study of 435 291 single nucleotide polymorphisms (SNPs) genotyped in 1738 MDD cases and 1802 controls selected to be at low liability for MDD. Of the top 200, 11 signals localized to a 167 kb region overlapping the gene piccolo (PCLO, whose protein product localizes to the cytomatrix of the presynaptic active zone and is important in monoaminergic neurotransmission in the brain) with P-values of 7.7 × 10−7 for rs2715148 and 1.2 × 10−6 for rs2522833. We undertook replication of SNPs in this region in five independent samples (6079 MDD independent cases and 5893 controls) but no SNP exceeded the replication significance threshold when all replication samples were analyzed together. However, there was heterogeneity in the replication samples, and secondary analysis of the original sample with the sample of greatest similarity yielded P=6.4 × 10−8 for the nonsynonymous SNP rs2522833 that gives rise to a serine to alanine substitution near a C2 calcium-binding domain of the PCLO protein. With the integrated replication effort, we present a specific hypothesis for further studies.

Meta-Analysis of Genome-Wide Association Data of Bipolar Disorder and Major Depressive Disorder

Meta-analysis of genome-wide association data of bipolar disorder and major depressive disorder

Genome-wide association study of Tourette's syndrome

Tourettes syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10−8); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10−6). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10−7 for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

The association between low vitamin D and depressive disorders

It has been hypothesized that hypovitaminosis D is associated with depression but epidemiological evidence is limited. We investigated the association between depressive disorders and related clinical characteristics with blood concentrations of 25-hydroxyvitamin D [25(OH)D] in a large cohort. The sample consisted of participants (aged 18–65 years) from the Netherlands Study of Depression and Anxiety (NESDA) with a current (N=1102) or remitted (N=790) depressive disorder (major depressive disorder, dysthymia) defined according to DSM-IV criteria, and healthy controls (N=494). Serum levels of 25(OH)D measured and analyzed in multivariate analyses adjusting for sociodemographics, sunlight, urbanization, lifestyle and health. Of the sample, 33.6% had deficient or insufficient serum 25(OH)D (<50 nmol l−1). As compared with controls, lower 25(OH)D levels were found in participants with current depression (P=0.001, Cohen’s d=0.21), particularly in those with the most severe symptoms (P=0.001, Cohen’s d=0.44). In currently depressed persons, 25(OH)D was inversely associated with symptom severity (β=−0.19, s.e.=0.07, P=0.003) suggesting a dose-response gradient, and with risk (relative risk=0.90, 95% confidence interval=0.82–0.99, P=0.03) of having a depressive disorders at 2-year follow-up. This large cohort study indicates that low levels of 25(OH)D were associated to the presence and severity of depressive disorder suggesting that hypovitaminosis D may represent an underlying biological vulnerability for depression. Future studies should elucidate whether—the highly prevalent—hypovitaminosis D could be cost-effectively treated as part of preventive or treatment interventions for depression.

Gene expression in major depressive disorder

The search for genetic variants underlying major depressive disorder (MDD) has not yet provided firm leads to its underlying molecular biology. A complementary approach is to study gene expression in relation to MDD. We measured gene expression in peripheral blood from 1848 subjects from The Netherlands Study of Depression and Anxiety. Subjects were divided into current MDD (N=882), remitted MDD (N=635) and control (N=331) groups. MDD status and gene expression were measured again 2 years later in 414 subjects. The strongest gene expression differences were between the current MDD and control groups (129 genes at false-discovery rate, FDR<0.1). Gene expression differences across MDD status were largely unrelated to antidepressant use, inflammatory status and blood cell counts. Genes associated with MDD were enriched for interleukin-6 (IL-6)-signaling and natural killer (NK) cell pathways. We identified 13 gene expression clusters with specific clusters enriched for genes involved in NK cell activation (downregulated in current MDD, FDR=5.8 × 10−5) and IL-6 pathways (upregulated in current MDD, FDR=3.2 × 10−3). Longitudinal analyses largely confirmed results observed in the cross-sectional data. Comparisons of gene expression results to the Psychiatric Genomics Consortium (PGC) MDD genome-wide association study results revealed overlap with DVL3. In conclusion, multiple gene expression associations with MDD were identified and suggest a measurable impact of current MDD state on gene expression. Identified genes and gene clusters are enriched with immune pathways previously associated with the etiology of MDD, in line with the immune suppression and immune activation hypothesis of MDD.

Further confirmation of the association between anxiety and CTNND2: replication in humans

The rat genome sequencing and mapping consortium found evidence for an association between the catenin‐δ2 gene (CTNND2) and anxious behaviour. We replicated these results in humans by carrying out a genetic association test in patients with panic disorder, social phobia, generalized anxiety disorder and/or agoraphobia (N = 1714) and controls (N = 4125). We further explored the association between CTNND2 and other psychiatric disorders based on publicly available genome‐wide association results. A gene‐based test showed that single nucleotide polymorphisms (SNPs) in CTNND2 have a significantly increased signal (P < 1e−5) and decreased P‐values. Single nucleotide polymorphism rs1012176 showed the strongest association with any anxiety disorder (odds ratio: 0.8128, SE = 0.063, P = 0.00099), but this effect was not significant after correction for multiple testing. In available genome‐wide association results from the Psychiatric Genomics Consortium we found that SNPs in CTNND2 collectively showed an increased signal for schizophrenia (P < 1e−5) and major depressive disorder (P < 1e−5), but not for bipolar disorder. These signals remained significant after correction for potential confounders. The association between CTNND2 and anxiety was not strong enough to be picked up in the current generation of human genome‐wide analyses, indicating the usefulness of and need for animal genetic studies to identify candidate genes for further study in human samples.

Gene expression in major depressive disorder: ERRATUM

Correction to: Molecular Psychiatry advance online publication, 26 May 2015; doi:10.1038/mp.2015.57 Following publication of the above article, the authors noticed that the Supplementary Figure legends were not published with the paper. The legends accompany this erratum. In addition, the Supplementary Tables were originally presented as PDF files.

Large normal-range TBP and ATXN7 CAG repeat lengths are associated with increased lifetime risk of depression

Depression is one of the most prevalent and debilitating psychiatric disorders worldwide. Recently, we showed that both relatively short and relatively long cytosine–adenine–guanine (CAG) repeats in the huntingtin gene (HTT) are associated with an increased risk of lifetime depression. However, to what extent the variations in CAG repeat length in the other eight polyglutamine disease-associated genes (PDAGs) are associated with depression is still unknown. We determined the CAG repeat sizes of ATXN1, ATXN2, ATXN3, CACNA1A, ATXN7, TBP, ATN1 and AR in two well-characterized Dutch cohorts—the Netherlands Study of Depression and Anxiety and the Netherlands Study of Depression in Older Persons—including 2165 depressed and 1058 non-depressed individuals—aged 18–93 years. The association between PDAG CAG repeat size and the risk for depression was assessed via binary logistic regression. We found that the odds ratio (OR) for lifetime depression was significantly higher for individuals with >10, compared with subjects with ≤10, CAG repeats in both ATXN7 alleles (OR=1.90, confidence interval (CI) 1.26–2.85). For TBP we found a similar association: A CAG repeat length exceeding the median in both alleles was associated with an increased risk for lifetime depression (OR=1.33, CI 1.00–1.76). In conclusion, we observed that carriers of either ATXN7 or TBP alleles with relatively large CAG repeat sizes in both alleles had a substantially increased risk of lifetime depression. Our findings provide critical evidence for the notion that repeat polymorphisms can act as complex genetic modifiers of depression.

Genetic loci influencing kidney function and chronic kidney disease.

Using genome-wide association, we identify common variants at 2p12–p13, 6q26, 17q23 and 19q13 associated with serum creatinine, a marker of kidney function (P = 10−10 to 10−15). Of these, rs10206899 (near NAT8, 2p12–p13) and rs4805834 (near SLC7A9, 19q13) were also associated with chronic kidney disease (P = 5.0 × 10−5 and P = 3.6 × 10−4, respectively). Our findings provide insight into metabolic, solute and drug-transport pathways underlying susceptibility to chronic kidney disease.