Search

We conducted a genome-wide association study (GWAS) on multiple sclerosis (MS) susceptibility in German cohorts with 4888 cases and 10,395 controls. In addition to associations within the major histocompatibility complex (MHC) region, 15 non-MHC loci reached genome-wide significance. Four of these loci are novel MS susceptibility loci. They map to the genes L3MBTL3, MAZ, ERG, and SHMT1. The lead variant at SHMT1 was replicated in an independent Sardinian cohort. Products of the genes L3MBTL3, MAZ, and ERG play important roles in immune cell regulation. SHMT1 encodes a serine hydroxymethyltransferase catalyzing the transfer of a carbon unit to the folate cycle. This reaction is required for regulation of methylation homeostasis, which is important for establishment and maintenance of epigenetic signatures. Our GWAS approach in a defined population with limited genetic substructure detected associations not found in larger, more heterogeneous cohorts, thus providing new clues regarding MS pathogenesis.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, 19 associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biological pathways.

Depression risk is exacerbated by genetic factors and stress exposure; however, the biological mechanisms through which these factors interact to confer depression risk are poorly understood. One putative biological mechanism implicates variability in the ability of cortisol, released in response to stress, to trigger a cascade of adaptive genomic and non-genomic processes through glucocorticoid receptor (GR) activation. Here, we demonstrate that common genetic variants in long-range enhancer elements modulate the immediate transcriptional response to GR activation in human blood cells. These functional genetic variants increase risk for depression and co-heritable psychiatric disorders. Moreover, these risk variants are associated with inappropriate amygdala reactivity, a transdiagnostic psychiatric endophenotype and an important stress hormone response trigger. Network modeling and animal experiments suggest that these genetic differences in GR-induced transcriptional activation may mediate the risk for depression and other psychiatric disorders by altering a network of functionally related stress-sensitive genes in blood and brain., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

The efflux pump P-glycoprotein (P-gp), a gene product of the ABCB1 gene, plays a pivotal role in the transfer of various molecules across the blood-brain barrier. P-gp protects the brain by selectively extruding its substrates, including certain antidepressive drugs, thereby limiting their uptake into the brain. Uhr et al. [2008] first showed that ABCB1 variants predicted the remission to antidepressants with P-gp substrate properties in patients suffering from major depression (MD). Other studies investigating the influence of ABCB1 polymorphisms on antidepressant treatment response produced inconclusive results. In this meta-analysis, we systematically summarized 16 pharmacogenetic studies focused on the association of ABCB1 variants and antidepressant treatment outcome in patients with MD (overall n = 2695). We investigated the association of treatment outcome and six ABCB1 single nucleotide polymorphisms (SNPs): rs2032583, rs2235015, rs2235040, rs1045642, rs2032582, rs1128503. We stratified for admission status, ethnicity, and prescription of concomitant medication. SNP rs2032583 showed a nominally significant association across all studies (P = 0.035, SNP was studied in a total of 2,037 patients) and a significant Bonferroni-corrected association among inpatients (P = 1.5 × 10(-05) , n = 485). Also SNP rs2235015 was significantly associated with antidepressant treatment outcome withstanding Bonferroni correction (P = 3.0 × 10(-04) ) among inpatients in a smaller subsample (n = 195). There were no significant associations of the other SNPs tested with antidepressant treatment outcome. Future pharmacogenetic association studies should focus on the role of the ABCB1 SNP rs2032583 in antidepressant outcome prediction., (© 2015 Wiley Periodicals, Inc.)

Genomic copy number variants (CNVs) have been implicated in multiple psychiatric disorders, but not much is known about their influence on anxiety disorders specifically. Using next-generation sequencing (NGS) and two additional array-based genotyping approaches, we detected CNVs in a mouse model consisting of two inbred mouse lines showing high (HAB) and low (LAB) anxiety-related behavior, respectively. An influence of CNVs on gene expression in the central (CeA) and basolateral (BLA) amygdala, paraventricular nucleus (PVN), and cingulate cortex (Cg) was shown by a two-proportion Z-test (p = 1.6 x 10-31), with a positive correlation in the CeA (p = 0.0062), PVN (p = 0.0046) and Cg (p = 0.0114), indicating a contribution of CNVs to the genetic predisposition to trait anxiety in the specific context of HAB/LAB mice. In order to confirm anxiety-relevant CNVs and corresponding genes in a second mouse model, we further examined CD-1 outbred mice. We revealed the distribution of CNVs by genotyping 64 CD 1 individuals using a high-density genotyping array (Jackson Laboratory). 78 genes within those CNVs were identified to show nominally significant association (48 genes), or a statistical trend in their association (30 genes) with the time animals spent on the open arms of the elevated plus-maze (EPM). Fifteen of them were considered promising candidate genes of anxiety-related behavior as we could show a significant overlap (permutation test, p = 0.0051) with genes within HAB/LAB CNVs. Thus, here we provide what is to our knowledge the first extensive catalogue of CNVs in CD-1 mice and potential corresponding candidate genes linked to anxiety-related behavior in mice.

Background: Obesity is strongly associated with major depressive disorder (MDD) and various other diseases. Genome-wide association studies have identified multiple risk loci robustly associated with body mass index (BMI). In this study, we aimed to investigate whether a genetic risk score (GRS) combining multiple BMI risk loci might have utility in prediction of obesity in patients with MDD., Methods: Linear and logistic regression models were conducted to predict BMI and obesity, respectively, in three independent large case-control studies of major depression (Radiant, GSK-Munich, PsyCoLaus). The analyses were first performed in the whole sample and then separately in depressed cases and controls. An unweighted GRS was calculated by summation of the number of risk alleles. A weighted GRS was calculated as the sum of risk alleles at each locus multiplied by their effect sizes. Receiver operating characteristic (ROC) analysis was used to compare the discriminatory ability of predictors of obesity., Results: In the discovery phase, a total of 2,521 participants (1,895 depressed patients and 626 controls) were included from the Radiant study. Both unweighted and weighted GRS were highly associated with BMI (P < 0.001) but explained only a modest amount of variance. Adding 'traditional' risk factors to GRS significantly improved the predictive ability with the area under the curve (AUC) in the ROC analysis, increasing from 0.58 to 0.66 (95% CI, 0.62-0.68; χ(2) = 27.68; P < 0.0001). Although there was no formal evidence of interaction between depression status and GRS, there was further improvement in AUC in the ROC analysis when depression status was added to the model (AUC = 0.71; 95% CI, 0.68-0.73; χ(2) = 28.64; P <0.0001). We further found that the GRS accounted for more variance of BMI in depressed patients than in healthy controls. Again, GRS discriminated obesity better in depressed patients compared to healthy controls. We later replicated these analyses in two independent samples (GSK-Munich and PsyCoLaus) and found similar results., Conclusions: A GRS proved to be a highly significant predictor of obesity in people with MDD but accounted for only modest amount of variance. Nevertheless, as more risk loci are identified, combining a GRS approach with information on non-genetic risk factors could become a useful strategy in identifying MDD patients at higher risk of developing obesity.

Dyslexia is one of the most common childhood disorders with a prevalence of around 5-10% in school-age children. Although an important genetic component is known to have a role in the aetiology of dyslexia, we are far from understanding the molecular mechanisms leading to the disorder. Several candidate genes have been implicated in dyslexia, including DYX1C1, DCDC2, KIAA0319, and the MRPL19/C2ORF3 locus, each with reports of both positive and no replications. We generated a European cross-linguistic sample of school-age children - the NeuroDys cohort - that includes more than 900 individuals with dyslexia, sampled with homogenous inclusion criteria across eight European countries, and a comparable number of controls. Here, we describe association analysis of the dyslexia candidate genes/locus in the NeuroDys cohort. We performed both case-control and quantitative association analyses of single markers and haplotypes previously reported to be dyslexia-associated. Although we observed association signals in samples from single countries, we did not find any marker or haplotype that was significantly associated with either case-control status or quantitative measurements of word-reading or spelling in the meta-analysis of all eight countries combined. Like in other neurocognitive disorders, our findings underline the need for larger sample sizes to validate possibly weak genetic effects.

Genome-wide association studies (GWAS) identified the MEIS1 locus for Restless Legs Syndrome (RLS), but causal single nucleotide polymorphisms (SNPs) and their functional relevance remain unknown. This locus contains a large number of highly conserved noncoding regions (HCNRs) potentially functioning as cis-regulatory modules. We analyzed these HCNRs for allele-dependent enhancer activity in zebrafish and mice and found that the risk allele of the lead SNP rs12469063 reduces enhancer activity in the Meis1 expression domain of the murine embryonic ganglionic eminences (GE). CREB1 binds this enhancer and rs12469063 affects its binding in vitro. In addition, MEIS1 target genes suggest a role in the specification of neuronal progenitors in the GE, and heterozygous Meis1-deficient mice exhibit hyperactivity, resembling the RLS phenotype. Thus, in vivo and in vitro analysis of a common SNP with small effect size showed allele-dependent function in the prospective basal ganglia representing the first neurodevelopmental region implicated in RLS.

Approximately 20% of individuals with Parkinson's disease (PD) report a positive family history. Yet, a large portion of causal and disease-modifying variants is still unknown. We used exome sequencing in two affected individuals from a family with late-onset familial PD followed by frequency assessment in 975 PD cases and 1014 ethnically-matched controls and linkage analysis to identify potentially causal variants. Based on the predicted penetrance and the frequencies, a variant in PLXNA4 proved to be the best candidate and PLXNA4 was screened for additional variants in 862 PD cases and 940 controls, revealing an excess of rare non-synonymous coding variants in PLXNA4 in individuals with PD. Although we cannot conclude that the variant in PLXNA4 is indeed the causative variant, these findings are interesting in the light of a surfacing role of axonal guidance mechanisms in neurodegenerative disorders but, at the same time, highlight the difficulties encountered in the study of rare variants identified by next-generation sequencing in diseases with autosomal dominant or complex patterns of inheritance.

Objectives: DMBT is an antibacterial pattern recognition and scavenger receptor. In this study, we analyzed the role of DMBT1 single nucleotide polymorphisms (SNPs) regarding inflammatory bowel disease (IBD) susceptibility and examined their functional impact on transcription factor binding and downstream gene expression., Methods: Seven SNPs in the DMBT1 gene region were analyzed in 2073 individuals including 818 Crohn's disease (CD) patients and 972 healthy controls in two independent case-control panels. Comprehensive epistasis analyses for the known CD susceptibility genes NOD2, IL23R and IL27 were performed. The influence of IL23R variants on DMBT1 expression was analyzed. Functional analysis included siRNA transfection, quantitative PCR, western blot, electrophoretic mobility shift and luciferase assays., Results: IL-22 induces DMBT1 protein expression in intestinal epithelial cells dependent on STAT3, ATF-2 and CREB1. IL-22 expression-modulating, CD risk-associated IL23R variants influence DMBT1 expression in CD patients and DMBT1 levels are increased in the inflamed intestinal mucosa of CD patients. Several DMBT1 SNPs were associated with CD susceptibility. SNP rs2981804 was most strongly associated with CD in the combined panel (p = 3.0 × 10(-7), OR 1.42; 95% CI 1.24-1.63). All haplotype groups tested showed highly significant associations with CD (including omnibus P-values as low as 6.1 × 10(-18)). The most strongly CD risk-associated, non-coding DMBT1 SNP rs2981804 modifies the DNA binding sites for the transcription factors CREB1 and ATF-2 and the respective genomic region comprising rs2981804 is able to act as a transcriptional regulator in vitro. Intestinal DMBT1 expression is decreased in CD patients carrying the rs2981804 CD risk allele., Conclusion: We identified novel associations of DMBT1 variants with CD susceptibility and discovered a novel functional role of rs2981804 in regulating DMBT1 expression. Our data suggest an important role of DMBT1 in CD pathogenesis.

Most psychiatric disorders are moderately to highly heritable. The degree to which genetic variation is unique to individual disorders or shared across disorders is unclear. To examine shared genetic etiology, we use genome-wide genotype data from the Psychiatric Genomics Consortium (PGC) for cases and controls in schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorders (ASD) and attention-deficit/hyperactivity disorder (ADHD). We apply univariate and bivariate methods for the estimation of genetic variation within and covariation between disorders. SNPs explained 17-29% of the variance in liability. The genetic correlation calculated using common SNPs was high between schizophrenia and bipolar disorder (0.68 ± 0.04 s.e.), moderate between schizophrenia and major depressive disorder (0.43 ± 0.06 s.e.), bipolar disorder and major depressive disorder (0.47 ± 0.06 s.e.), and ADHD and major depressive disorder (0.32 ± 0.07 s.e.), low between schizophrenia and ASD (0.16 ± 0.06 s.e.) and non-significant for other pairs of disorders as well as between psychiatric disorders and the negative control of Crohn's disease. This empirical evidence of shared genetic etiology for psychiatric disorders can inform nosology and encourages the investigation of common pathophysiologies for related disorders.

Background: Data from clinical studies and results from animal models suggest an involvement of the neurotrophin system in the pathology of depression and antidepressant treatment response. Genetic variations within the genes coding for the brain-derived neurotrophic factor (BDNF) and its key receptor Trkb (NTRK2) may therefore influence the response to antidepressant treatment., Methods: We performed a single and multi-marker association study with antidepressant treatment outcome in 398 depressed Caucasian inpatients participating in the Munich Antidepressant Response Signature (MARS) project. Two Caucasian replication samples (N = 249 and N = 247) were investigated, resulting in a total number of 894 patients. 18 tagging SNPs in the BDNF gene region and 64 tagging SNPs in the NTRK2 gene region were genotyped in the discovery sample; 16 nominally associated SNPs were tested in two replication samples., Results: In the discovery analysis, 7 BDNF SNPs and 9 NTRK2 SNPs were nominally associated with treatment response. Three NTRK2 SNPs (rs10868223, rs1659412 and rs11140778) also showed associations in at least one replication sample and in the combined sample with the same direction of effects (Pcorr  = .018, Pcorr  = .015 and Pcorr  = .004, respectively). We observed an across-gene BDNF-NTRK2 SNP interaction for rs4923468 and rs1387926. No robust interaction of associated SNPs was found in an analysis of BDNF serum protein levels as a predictor for treatment outcome in a subset of 93 patients., Conclusions/limitations: Although not all associations in the discovery analysis could be unambiguously replicated, the findings of the present study identified single nucleotide variations in the BDNF and NTRK2 genes that might be involved in antidepressant treatment outcome and that have not been previously reported in this context. These new variants need further validation in future association studies.

Attention-deficit/hyperactivity disorder (ADHD) and dyslexia belong to the most common neuro-behavioral childhood disorders with prevalences of around 5% in school-aged children. It is estimated that 20-60% of individuals affected with ADHD also present with learning disorders. We investigated the comorbidity between ADHD symptoms and reading/spelling and math difficulties in two on-going population-based birth cohort studies. Children with ADHD symptoms were at significantly higher risk of also showing reading/spelling difficulties or disorder (Odds Ratio (OR) = 2.80, p = 6.59×10⁻¹³) as compared to children without ADHD symptoms. For math difficulties the association was similar (OR = 2.55, p = 3.63×10⁻⁰⁴). Our results strengthen the hypothesis that ADHD and learning disorders are comorbid and share, at least partially, the same underlying process. Up to date, it is not clear, on which exact functional processes this comorbidity is based.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Restless legs syndrome (RLS) is a common multifactorial disease. Some genetic risk factors have been identified. RLS susceptibility also has been related to iron. We therefore asked whether known iron-related genes are candidates for association with RLS and, vice versa, whether known RLS-associated loci influence iron parameters in serum. RLS/control samples (n = 954/1814 in the discovery step, 735/736 in replication 1, and 736/735 in replication 2) were tested for association with SNPs located within 4 Mb intervals surrounding each gene from a list of 111 iron-related genes using a discovery threshold of P = 5 × 10(-4). Two population cohorts (KORA F3 and F4 with together n = 3447) were tested for association of six known RLS loci with iron, ferritin, transferrin, transferrin-saturation, and soluble transferrin receptor. Results were negative. None of the candidate SNPs at the iron-related gene loci was confirmed significantly. An intronic SNP, rs2576036, of KATNAL2 at 18q21.1 was significant in the first (P = 0.00085) but not in the second replication step (joint nominal P-value = 0.044). Especially, rs1800652 (C282Y) in the HFE gene did not associate with RLS. Moreover, SNPs at the known RLS loci did not significantly affect serum iron parameters in the KORA cohorts. In conclusion, the correlation between RLS and iron parameters in serum may be weaker than assumed. Moreover, in a general power analysis, we show that genetic effects are diluted if they are transmitted via an intermediate trait to an end-phenotype. Sample size formulas are provided for small effect sizes.

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

Background & Aims: Genome-wide association studies identified the autophagy gene IRGM to be strongly associated with Crohn's disease (CD) but its impact in ulcerative colitis (UC), its phenotypic effects and potential epistatic interactions with other IBD susceptibility genes are less clear which we therefore analyzed in this study., Methodology/principal Findings: Genomic DNA from 2060 individuals including 817 CD patients, 283 UC patients, and 961 healthy, unrelated controls (all of Caucasian origin) was analyzed for six IRGM single nucleotide polymorphisms (SNPs) (rs13371189, rs10065172 = p.Leu105Leu, rs4958847, rs1000113, rs11747270, rs931058). In all patients, a detailed genotype-phenotype analysis and testing for epistasis with the three major CD susceptibility genes NOD2, IL23R and ATG16L1 were performed. Our analysis revealed an association of the IRGM SNPs rs13371189 (p = 0.02, OR 1.31 [95% CI 1.05-1.65]), rs10065172 = p.Leu105Leu (p = 0.016, OR 1.33 [95% CI 1.06-1.66]) and rs1000113 (p = 0.047, OR 1.27 [95% CI 1.01-1.61]) with CD susceptibility. There was linkage disequilibrium between these three IRGM SNPs. In UC, several IRGM haplotypes were weakly associated with UC susceptibility (p<0.05). Genotype-phenotype analysis revealed no significant associations with a specific IBD phenotype or ileal CD involvement. There was evidence for weak gene-gene-interaction between several SNPs of the autophagy genes IRGM and ATG16L1 (p<0.05), which, however, did not remain significant after Bonferroni correction., Conclusions/significance: Our results confirm IRGM as susceptibility gene for CD in the German population, supporting a role for the autophagy genes IRGM and ATG16L1 in the pathogenesis of CD.

Background: The hypothalamic-pituitary-adrenal (HPA) axis is essential to control physiological stress responses in mammals. Its dysfunction is related to several mental disorders, including anxiety and depression. The aim of this study was to identify genetic loci underlying the endocrine regulation of the HPA axis., Method: High (HAB) and low (LAB) anxiety-related behaviour mice were established by selective inbreeding of outbred CD-1 mice to model extremes in trait anxiety. Additionally, HAB vs. LAB mice exhibit comorbid characteristics including a differential corticosterone response upon stress exposure. We crossbred HAB and LAB lines to create F1 and F2 offspring. To identify the contribution of the endocrine phenotypes to the total phenotypic variance, we examined multiple behavioural paradigms together with corticosterone secretion-based phenotypes in F2 mice by principal component analysis. Further, to pinpoint the genomic loci of the quantitative trait of the HPA axis stress response, we conducted genome-wide multipoint oligogenic linkage analyses based on Bayesian Markov chain Monte Carlo approach as well as parametric linkage in three-generation pedigrees, followed by a two-dimensional scan for epistasis and association analysis in freely segregating F2 mice using 267 single-nucleotide polymorphisms (SNPs), which were identified to consistently differ between HAB and LAB mice as genetic markers., Results: HPA axis reactivity measurements and behavioural phenotypes were represented by independent principal components and demonstrated no correlation. Based on this finding, we identified one single quantitative trait locus (QTL) on chromosome 3 showing a very strong evidence for linkage (2ln (L-score) > 10, LOD > 23) and significant association (lowest Bonferroni adjusted p < 10-28) to the neuroendocrine stress response. The location of the linkage peak was estimated at 42.3 cM (95% confidence interval: 41.3 - 43.3 cM) and was shown to be in epistasis (p-adjusted < 0.004) with the locus at 35.3 cM on the same chromosome. The QTL harbours genes involved in steroid synthesis and cardiovascular effects., Conclusion: The very prominent effect on stress-induced corticosterone secretion of the genomic locus on chromosome 3 and its involvement in epistasis highlights the critical role of this specific locus in the regulation of the HPA axis.

Recent genome-wide association studies (GWAS) and metaanalyses revealed genetic associations for ANK3 (ankyrin 3) and CACNA1C (alpha 1C subunit of the L-type voltage gated calcium channel) with bipolar disorder (BPD). Several findings from clinical, epidemiological, and genetic studies point towards a common biological background of BPD and major depressive disorder (MDD). We were interested whether this also applies for ANK3 and CACNA1C and tested associations of single nucleotide polymorphisms (SNPs) in these genes with MDD in two Caucasian case-control samples. Sample 1 (Munich Antidepressant Response Signature Project/MARS - MDD) consisted of 720 depressed inpatients and 542 psychiatric healthy controls. Sample 2 (unipolar recurrent depression (URD)) consisted of 827 patients with URD and 860 psychiatric healthy controls. After stringent quality control we analyzed 262 SNPs (sample 1) and 504 SNPs (sample 2) and imputed further 5771 SNPs (sample 1) and 5534 SNPs (sample 2) from Hapmap Phase 2 data in the ANK3 and CACNA1C gene regions. Additionally, a metaanalysis of both samples was performed. Several SNPs in both genes were nominally associated with MDD with the highest association in the 3'-region of ANK3 (rs10994143, nominal p = 3.3*10(-4)) in the metaanalysis of both samples. None of these results remained significant after correction for multiple testing. No association of MDD with SNPs previously reported in BPD studies could be detected. By analyzing the LD-structure, our highest associated SNPs could not be linked to the SNPs previously reported in BPD. Regarding ANK3 and CACNA1C, our findings do not support a strong genetic link between BPD and MDD for these two genes., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Objectives: Glutamate has been implicated in the pathophysiology and treatment of mood disorders possibly by affecting the regulation of the hypothalamus-pituitary-adrenocortical (HPA) axis. Growing evidence suggests an important role of the metabotropic glutamate receptor 1 (mGlu1) in depression-related phenotypes. To test whether these findings can also be supported by human genetics data, we explored polymorphisms within the metabotropic glutamate receptor 1 gene (GRM1) for their association with unipolar depression (UPD) as well as with biological phenotypes of this disorder., Methods: We first tested the association of 43 tag-SNPs covering the GRM1 locus with UPD in 350 patients and 370 matched controls. We then investigated the effects of the associated SNPs on hippocampal glutamate levels estimated using ¹H-MR-spectroscopy (¹H-MRS) and on endocrine measures from the combined dexamethasone-suppression/CRH stimulation (dex/CRH) test., Results: Within the GRM1 locus, 22 SNPs showed nominally significant association with UPD, of which 6 withstood corrections for multiple testing (rs2268666 with best allelic p=7.0×10⁻⁵). Supportive evidence for an association with UPD was gained from a second independent sample with 904 patients and 1012 controls. Furthermore, patients homozygous for the non-risk genotypes showed reduced hippocampal glutamate levels as measured by ¹H-MRS, a more pronounced normalization of HPA-axis hyperactivity as well as a better antidepressant treatment outcome., Conclusions: These results suggest that the combination of genetic and biological markers may allow to subgroup patients into etiopathogenetically more relevant subcategories which could guide clinicians in their antidepressant treatment choices., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Emergence of suicidal ideation (TESI) during treatment with antidepressants in major depression led to a black box warning. We performed a genome-wide association study to identify genetic markers, which increase the risk for this serious side effect. TESI was evaluated in depressed in-patients (N=397) and defined by an emergence of suicidal thoughts during hospitalization without suicidal thoughts at admission using the suicide item (3) of the Hamilton Depression Rating Scale. Genotype distribution of 405.383 single-nucleotide polymorphisms (SNPs) in patients with TESI (N=32/8.1%) was compared to patients without increase in suicidal ideation (N=329/82.9%) and to a subgroup never reported suicidal ideation (N=79/19.9%). Top results were analyzed in an independent sample (N=501). None variant reached genome-wide significance, the best associated SNP was rs1630535 (p-value=1.3 × 10(-7)). The top 79 SNPs could be analyzed in an independent sample, and 14 variants showed nominal significant association with the same risk allele in the replication sample. A discriminant analysis classifying patients using these 79 SNPs revealed a 91% probability to classify TESI vs non-TESI cases correctly in the replication sample. Although our data need to be interpreted carefully owing to the small numbers in both cohorts, they suggest that a combination of genetic markers might indeed be used to identify patients at risk for TESI.

Background: Genome-wide association studies identified a PTGER4 expression-modulating region on chromosome 5p13.1 as Crohn's disease (CD) susceptibility region. The study aim was to test this association in a large cohort of patients with inflammatory bowel disease (IBD) and to elucidate genotypic and phenotypic interactions with other IBD genes., Methodology/principal Findings: A total of 7073 patients and controls were genotyped: 844 CD and 471 patients with ulcerative colitis and 1488 controls were analyzed for the single nucleotide polymorphisms (SNPs) rs4495224 and rs7720838 on chromosome 5p13.1. The study included two replication cohorts of North American (CD: n = 684; controls: n = 1440) and of German origin (CD: n = 1098; controls: n = 1048). Genotype-phenotype, epistasis and transcription factor binding analyses were performed. In the discovery cohort, an association of rs4495224 (p = 4.10×10⁻⁵; 0.76 [0.67-0.87]) and of rs7720838 (p = 6.91×10⁻⁴; 0.81 [0.71-0.91]) with susceptibility to CD was demonstrated. These associations were confirmed in both replication cohorts. In silico analysis predicted rs4495224 and rs7720838 as essential parts of binding sites for the transcription factors NF-κB and XBP1 with higher binding scores for carriers of the CD risk alleles, providing an explanation of how these SNPs might contribute to increased PTGER4 expression. There was no association of the PTGER4 SNPs with IBD phenotypes. Epistasis detected between 5p13.1 and ATG16L1 for CD susceptibility in the discovery cohort (p = 5.99×10⁻⁷ for rs7720838 and rs2241880) could not be replicated in both replication cohorts arguing against a major role of this gene-gene interaction in the susceptibility to CD., Conclusions/significance: We confirmed 5p13.1 as a major CD susceptibility locus and demonstrate by in silico analysis rs4495224 and rs7720838 as part of binding sites for NF-κB and XBP1. Further functional studies are necessary to confirm the results of our in silico analysis and to analyze if changes in PTGER4 expression modulate CD susceptibility.

Background: Genome-wide association studies identified PTPN2 (protein tyrosine phosphatase, non-receptor type 2) as susceptibility gene for inflammatory bowel diseases (IBD). However, the exact role of PTPN2 in Crohn's disease (CD) and ulcerative colitis (UC) and its phenotypic effect are unclear. We therefore performed a detailed genotype-phenotype and epistasis analysis of PTPN2 gene variants., Methodology/principal Findings: Genomic DNA from 2131 individuals of Caucasian origin (905 patients with CD, 318 patients with UC, and 908 healthy, unrelated controls) was analyzed for two SNPs in the PTPN2 region (rs2542151, rs7234029) for which associations with IBD were found in previous studies in other cohorts. Our analysis revealed a significant association of PTPN2 SNP rs2542151 with both susceptibility to CD (p = 1.95×10⁻⁵; OR 1.49 [1.34-1.79]) and UC (p = 3.87×10⁻², OR 1.31 [1.02-1.68]). Moreover, PTPN2 SNP rs7234029 demonstrated a significant association with susceptibility to CD (p = 1.30×10⁻³; OR 1.35 [1.13-1.62]) and a trend towards association with UC (p = 7.53×10⁻²; OR 1.26 [0.98-1.62]). Genotype-phenotype analysis revealed an association of PTPN2 SNP rs7234029 with a stricturing disease phenotype (B2) in CD patients (p = 6.62×10⁻³). Epistasis analysis showed weak epistasis between the ATG16L1 SNP rs2241879 and PTPN2 SNP rs2542151 (p = 0.024) in CD and between ATG16L1 SNP rs4663396 and PTPN2 SNP rs7234029 (p = 4.68×10⁻³) in UC. There was no evidence of epistasis between PTPN2 and NOD2 and PTPN2 and IL23R. In silico analysis revealed that the SNP rs7234029 modulates potentially the binding sites of several transcription factors involved in inflammation including GATA-3, NF-κB, C/EBP, and E4BP4., Conclusions/significance: Our data confirm the association of PTPN2 variants with susceptibility to both CD and UC, suggesting a common disease pathomechanism for these diseases. Given recent evidence that PTPN2 regulates autophagosome formation in intestinal epithelial cells, the potential link between PTPN2 and ATG16L1 should be further investigated.

Background: IL12B encodes the p40 subunit of IL-12, which is also part of IL-23. Recent genome-wide association studies identified IL12B and IL23R as susceptibility genes for inflammatory bowel disease (IBD). However, the phenotypic effects and potential gene-gene interactions of IL12B variants are largely unknown., Methodology/principal Findings: We analyzed IL12B gene variants regarding association with Crohn's disease (CD) and ulcerative colitis (UC). Genomic DNA from 2196 individuals including 913 CD patients, 318 UC patients and 965 healthy, unrelated controls was analyzed for four SNPs in the IL12B gene region (rs3212227, rs17860508, rs10045431, rs6887695). Our analysis revealed an association of the IL12B SNP rs6887695 with susceptibility to IBD (p = 0.035; OR 1.15 [95% CI 1.01-1.31] including a trend for rs6887695 for association with CD (OR 1.41; [0.99-1.31], p = 0.066) and UC (OR 1.18 [0.97-1.43], p = 0.092). CD patients, who were homozygous C/C carriers of this SNP, had significantly more often non-stricturing, non-penetrating disease than carriers of the G allele (p = 6.8×10(-5); OR = 2.84, 95% CI 1.66-4.84), while C/C homozygous UC patients had less often extensive colitis than G allele carriers (p = 0.029; OR = 0.36, 95% CI 0.14-0.92). In silico analysis predicted stronger binding of the minor C allele of rs6887695 to the transcription factor RORα which is involved in Th17 differentiation. Differences regarding the binding to the major and minor allele sequence of rs6887695 were also predicted for the transcription factors HSF1, HSF2, MZF1 and Oct-1. Epistasis analysis revealed weak epistasis of the IL12B SNP rs6887695 with several SNPs (rs11889341, rs7574865, rs7568275, rs8179673, rs10181656, rs7582694) in the STAT4 gene which encodes the major IL-12 downstream transcription factor STAT4 (p<0.05) but there was no epistasis between IL23R and IL12B variants., Conclusions/significance: The IL12B SNP rs6887695 modulates the susceptibility and the phenotype of IBD, although the effect on IBD susceptibilty is less pronounced than that of IL23R gene variants.

In 2008, van den Oord et al. identified in a genome-wide association study the MAM domain containing glycosylphosphatidylinositol anchor 2 gene (MDGA2 or MAMDC1) as a new candidate for neuroticism. In addition to the replication attempt of this association, we further investigated the role of MDGA2 with respect to harm avoidance (HA), a personality trait highly related to neuroticism. In a sample of mentally healthy volunteers (n=541) and depressed patients (n=199), neuroticism and HA were assessed with Eysenck's Revised Personality Questionnaire and Cloninger's Tridimensional Personality Questionnaire. Genotypic information (Illumina Bead Chip HumanHap300) of 100 single nucleotide polymorphisms (SNPs) located in the MDGA2 gene (±5 kb) was available, and additional four SNPs for replication were imputed. We were able to replicate the association between MDGA2 and neuroticism for the strongest SNPs of the genome-wide association study. It could further be shown that volunteers homozygous for the T-allele of SNP rs2416054 showed higher scores in the HA4 subscale 'Fatigability and asthenia' (Pnominal=0.0006), remaining significant after correction for multiple testing (Pwy-corrected=0.045). The same SNP also showed an association with HA4 in the patients' sample (Pnominal=0.03). Our finding provides further support for a link between variants in the MDGA2 gene and specific neuroticism-related phenotypes.

Background: The pregnane X receptor (PXR/NR1I2) is an important regulator of xenobiotic metabolism and intestinal integrity. However, there are controversial studies on the role of PXR/NR1I2 in inflammatory bowel disease (IBD). We therefore initiated the largest analysis to date on PXR/NR1I2 gene variants in IBD patients., Methods: Genomic DNA from 2823 individuals of Caucasian origin including 859 patients with Crohn's disease (CD), 464 patients with ulcerative colitis (UC), and 1500 healthy, unrelated controls was analyzed for eight PXR/NR1I2 single nucleotide polymorphisms (SNPs) (rs12721602 (-25564), rs3814055 (-25385), rs1523128 (-24756), rs1523127 (-24381), rs45610735 = p.Gly36Arg (+106), rs6785049 (+7635), rs2276707 (+8055), and rs3814057 (+11156)). In addition, detailed haplotype and genotype-phenotype analyses were performed., Results: The PXR/NR1I2 SNP rs2276707 was weakly associated with UC susceptibility (P = 0.01; odds ratio [OR] 1.27 [1.06-1.52]). None of the other PXR/NR1I2 SNPs were associated with UC or CD susceptibility. However, several rare PXR/NR1I2 haplotypes were highly associated with CD susceptibility. In CD, the strongest disease association was found for a haplotype consisting of the SNPs rs12721602-rs3814055-rs1523128-rs1523127-rs12721607-rs6785049-rs2276707-rs3814057 (omnibus P-value: 6.50 × 10(-15)) which was found in two separate cohorts (cohort I = discovery cohort: CD: n = 492, controls: n = 793; P = 4.51 × 10(-17); Bonferroni corrected: P = 1.27 × 10(-15); cohort II = replication cohort: CD: n = 367, controls: n = 707; P = 7.12 × 10(-4); P(corr) = 1.99 × 10(-2))., Conclusions: Several PXR/NR1I2 haplotypes contribute to CD susceptibility, suggesting a role for PXR in the IBD pathogenesis of a certain patient subcohort. Given the accumulating evidence for an important role of PXR in intestinal inflammation, further analyses are required to investigate the functional and pharmacogenetic implications of these PXR/NR1I2 gene variants in IBD., (Copyright © 2010 Crohn's & Colitis Foundation of America, Inc.)

Restless legs syndrome (RLS) is a sensorimotor disorder with an age-dependent prevalence of up to 10% in the general population above 65 years of age. Affected individuals suffer from uncomfortable sensations and an urge to move in the lower limbs that occurs mainly in resting situations during the evening or at night. Moving the legs or walking leads to an improvement of symptoms. Concomitantly, patients report sleep disturbances with consequences such as reduced daytime functioning. We conducted a genome-wide association study (GWA) for RLS in 922 cases and 1,526 controls (using 301,406 SNPs) followed by a replication of 76 candidate SNPs in 3,935 cases and 5,754 controls, all of European ancestry. Herein, we identified six RLS susceptibility loci of genome-wide significance, two of them novel: an intergenic region on chromosome 2p14 (rs6747972, P = 9.03 × 10(-11), OR = 1.23) and a locus on 16q12.1 (rs3104767, P = 9.4 × 10(-19), OR = 1.35) in a linkage disequilibrium block of 140 kb containing the 5'-end of TOX3 and the adjacent non-coding RNA BC034767., Competing Interests: The authors have declared that no competing interests exist.

Background: The carcinoembryonic antigen-related cell adhesion molecule 6 (CEACAM6) acts as a receptor for adherent-invasive E. coli (AIEC) and its ileal expression is increased in patients with Crohn's disease (CD). Given its contribution to the pathogenesis of CD, we aimed to investigate the role of genetic variants in the CEACAM6 region in patients with inflammatory bowel diseases (IBD)., Methodology: In this study, a total of 2,683 genomic DNA samples (including DNA from 858 CD patients, 475 patients with ulcerative colitis (UC), and 1,350 healthy, unrelated controls) was analyzed for eight CEACAM6 SNPs (rs10415946, rs1805223 = p.Pro42Pro, rs4803507, rs4803508, rs11548735 = p.Gly239Val, rs7246116 = pHis260His, rs2701, rs10416839). In addition, a detailed haplotype analysis and genotype-phenotype analysis were performed. Overall, our genotype analysis did not reveal any significant association of the investigated CEACAM6 SNPs and haplotypes with CD or UC susceptibility, although certain CEACAM6 SNPs modulated CEACAM6 expression in intestinal epithelial cell lines. Despite its function as receptor of AIEC in ileal CD, we found no association of the CEACAM6 SNPs with ileal or ileocolonic CD. Moreover, there was no evidence of epistasis between the analyzed CEACAM6 variants and the main CD-associated NOD2, IL23R and ATG16L1 variants., Conclusions: This study represents the first detailed analysis of CEACAM6 variants in IBD patients. Despite its important role in bacterial attachment in ileal CD, we could not demonstrate a role for CEACAM6 variants in IBD susceptibility or regarding an ileal CD phenotype. Further functional studies are required to analyze if these gene variants modulate ileal bacterial attachment.

Major depression (MD) is one of the most prevalent psychiatric disorders and a leading cause of loss in work productivity. A combination of genetic and environmental risk factors probably contributes to MD. We present data from a genome-wide association study revealing a neuron-specific neutral amino acid transporter (SLC6A15) as a susceptibility gene for MD. Risk allele carrier status in humans and chronic stress in mice were associated with a downregulation of the expression of this gene in the hippocampus, a brain region implicated in the pathophysiology of MD. The same polymorphisms also showed associations with alterations in hippocampal volume and neuronal integrity. Thus, decreased SLC6A15 expression, due to genetic or environmental factors, might alter neuronal circuits related to the susceptibility for MD. Our convergent data from human genetics, expression studies, brain imaging, and animal models suggest a pathophysiological mechanism for MD that may be accessible to drug targeting., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Detection of epistatic interaction between loci has been postulated to provide a more in-depth understanding of the complex biological and biochemical pathways underlying human diseases. Studying the interaction between two loci is the natural progression following traditional and well-established single locus analysis. However, the added costs and time duration required for the computation involved have thus far deterred researchers from pursuing a genome-wide analysis of epistasis. In this paper, we propose a method allowing such analysis to be conducted very rapidly. The method, dubbed EPIBLASTER, is applicable to case-control studies and consists of a two-step process in which the difference in Pearson's correlation coefficients is computed between controls and cases across all possible SNP pairs as an indication of significant interaction warranting further analysis. For the subset of interactions deemed potentially significant, a second-stage analysis is performed using the likelihood ratio test from the logistic regression to obtain the P-value for the estimated coefficients of the individual effects and the interaction term. The algorithm is implemented using the parallel computational capability of commercially available graphical processing units to greatly reduce the computation time involved. In the current setup and example data sets (211 cases, 222 controls, 299468 SNPs; and 601 cases, 825 controls, 291095 SNPs), this coefficient evaluation stage can be completed in roughly 1 day. Our method allows for exhaustive and rapid detection of significant SNP pair interactions without imposing significant marginal effects of the single loci involved in the pair.

Background: Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients., Methodology/principal Findings: Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn's disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction., Conclusions/significance: Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.

It has been repeatedly shown that mismatch negativity (MMN), an event related potential measurement, reveals differences between dyslexic children and age-matched controls. MMN reflects the automatic detection of deviance between a stream of incoming sounds presented to the passive listener, and deficits in MMN (i.e. attenuated amplitudes) have been especially reported in dyslexia for detecting differences between speech sounds (e.g./ba/vs./da/). We performed an association analysis in 200 dyslexic children. This analysis focused on two MMN components, an early MMN (188-300 ms) and a late MMN (300-710 ms), and the dyslexia candidate genes KIAA0319 and DCDC2 on chromosome 6. Additionally, we imputed rare variants located in this region based on the 1000 genomes project. We identified four rare variants that were significantly associated with the late MMN component. For three of these variants, which were in high LD to each other, genotyping confirmed the association signal. Our results point to an association between late MMN and rare variants in a candidate gene region for dyslexia.

Background: The aims were to analyze two novel NOD2 variants (rs2066843 and rs2076756) in a large cohort of patients with inflammatory bowel disease and to elucidate phenotypic consequences., Methodology/principal Findings: Genomic DNA from 2700 Caucasians including 812 patients with Crohn's disease (CD), 442 patients with ulcerative colitis (UC), and 1446 healthy controls was analyzed for the NOD2 SNPs rs2066843 and rs2076756 and the three main CD-associated NOD2 variants p.Arg702Trp (rs2066844), p.Gly908Arg (rs2066847), and p.Leu1007fsX1008 (rs2066847). Haplotype and genotype-phenotype analyses were performed. The SNPs rs2066843 (p = 3.01×10(-5), OR 1.48, [95% CI 1.23-1.78]) and rs2076756 (p = 4.01×10(-6); OR 1.54, [95% CI 1.28-1.86]) were significantly associated with CD but not with UC susceptibility. Haplotype analysis revealed a number of significant associations with CD susceptibility with omnibus p values <10(-10). The SNPs rs2066843 and rs2076756 were in linkage disequilibrium with each other and with the three main CD-associated NOD2 mutations (D'>0.9). However, in CD, SNPs rs2066843 and rs2076756 were more frequently observed than the other three common NOD2 mutations (minor allele frequencies for rs2066843 and rs2076756: 0.390 and 0.380, respectively). In CD patients homozygous for these novel NOD2 variants, genotype-phenotype analysis revealed higher rates of a penetrating phenotype (rs2076756: p = 0.015) and fistulas (rs2076756: p = 0.015) and significant associations with CD-related surgery (rs2076756: p = 0.003; rs2066843: p = 0.015). However, in multivariate analysis only disease localization (p<2×10(-16)) and behaviour (p = 0.02) were significantly associated with the need for surgery., Conclusion/significance: The NOD2 variants rs2066843 and rs2076756 are novel and common CD susceptibility gene variants.