Darina Czamara

Genome-Wide Association Study Identifies Novel Restless Legs Syndrome Susceptibility Loci on 2p14 and 16q12.1

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.

EPIBLASTER-fast exhaustive two-locus epistasis detection strategy using graphical processing units

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 Pearsons 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.

Restless Legs Syndrome-associated intronic common variant in Meis1 alters enhancer function in the developing telencephalon

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.

Genome-Wide Association Study of Antidepressant Treatment-Emergent Suicidal Ideation

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.

IRGM variants and susceptibility to inflammatory bowel disease in the German population.

Background & Aims Genome-wide association studies identified the autophagy gene IRGM to be strongly associated with Crohns 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.

Replication and meta-analysis of TMEM132D gene variants in panic disorder.

A recent genome-wide association study in patients with panic disorder (PD) identified a risk haplotype consisting of two single-nucleotide polymorphisms (SNPs) (rs7309727 and rs11060369) located in intron 3 of TMEM132D to be associated with PD in three independent samples. Now we report a subsequent confirmation study using five additional PD case–control samples (n=1670 cases and n=2266 controls) assembled as part of the Panic Disorder International Consortium (PanIC) study for a total of 2678 cases and 3262 controls in the analysis. In the new independent samples of European ancestry (EA), the association of rs7309727 and the risk haplotype rs7309727–rs11060369 was, indeed, replicated, with the strongest signal coming from patients with primary PD, that is, patients without major psychiatric comorbidities (n=1038 cases and n=2411 controls). This finding was paralleled by the results of the meta-analysis across all samples, in which the risk haplotype and rs7309727 reached P-levels of P=1.4e−8 and P=1.1e−8, respectively, when restricting the samples to individuals of EA with primary PD. In the Japanese sample no associations with PD could be found. The present results support the initial finding that TMEM132D gene contributes to genetic susceptibility for PD in individuals of EA. Our results also indicate that patient ascertainment and genetic background could be important sources of heterogeneity modifying this association signal in different populations.

MEIS1 and BTBD9: genetic association with restless leg syndrome in end stage renal disease

Background Restless legs syndrome (RLS) is a sleep related movement disorder that occurs both in an idiopathic form and in symptomatic varieties. RLS is a frequent and distressing comorbidity in end stage renal disease (ESRD). For idiopathic RLS (iRLS), genetic risk factors have been identified, but their role in RLS in ESRD has not been investigated yet. Therefore, a case–control association study of these variants in ESRD patients was performed. Methods The study genotyped 10 iRLS associated variants at four loci encompassing the genes MEIS1, BTBD9, MAP2K5/SKOR1, and PTPRD, in two independent case–control samples from Germany and Greece using multiplex PCR and MALDI-TOF (matrix assisted laser desorption/ionisation time-of-flight) mass spectrometry. Statistical analysis was performed as logistic regression with age and gender as covariates. For the combined analysis a Cochran–Mantel–Haenszel test was applied. Results The study included 200 RLS-positive and 443 RLS-negative ESRD patients in the German sample, and 141 and 393 patients, respectively, in the Greek sample. In the German sample, variants in MEIS1 and BTBD9 were associated with RLS in ESRD (Pnom≤0.004, ORs 1.52 and 1.55), whereas, in the Greek sample, there was a trend for association to MAP2K5/SKOR1 and BTBD9 (Pnom≤0.08, ORs 1.41 and 1.33). In the combined analysis including all samples, BTBD9 was associated after correction for multiple testing (Pcorrected=0.0013, OR 1.47). Conclusions This is the first demonstration of a genetic influence on RLS in ESRD patients with BTBD9 being significantly associated. The extent of the genetic predisposition could vary between different subgroups of RLS in ESRD.

Genetic Differences in the Immediate Transcriptome Response to Stress Predict Risk-Related Brain Function and Psychiatric Disorders

Summary 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. Video Abstract

MEIS1 and BTBD9 - genetic association with RLS in end-stage renal disease

Background Restless legs syndrome (RLS) is a sleep related movement disorder that occurs both in an idiopathic form and in symptomatic varieties. RLS is a frequent and distressing comorbidity in end stage renal disease (ESRD). For idiopathic RLS (iRLS), genetic risk factors have been identified, but their role in RLS in ESRD has not been investigated yet. Therefore, a case–control association study of these variants in ESRD patients was performed. Methods The study genotyped 10 iRLS associated variants at four loci encompassing the genes MEIS1, BTBD9, MAP2K5/SKOR1, and PTPRD, in two independent case–control samples from Germany and Greece using multiplex PCR and MALDI-TOF (matrix assisted laser desorption/ionisation time-of-flight) mass spectrometry. Statistical analysis was performed as logistic regression with age and gender as covariates. For the combined analysis a Cochran–Mantel–Haenszel test was applied. Results The study included 200 RLS-positive and 443 RLS-negative ESRD patients in the German sample, and 141 and 393 patients, respectively, in the Greek sample. In the German sample, variants in MEIS1 and BTBD9 were associated with RLS in ESRD (Pnom≤0.004, ORs 1.52 and 1.55), whereas, in the Greek sample, there was a trend for association to MAP2K5/SKOR1 and BTBD9 (Pnom≤0.08, ORs 1.41 and 1.33). In the combined analysis including all samples, BTBD9 was associated after correction for multiple testing (Pcorrected=0.0013, OR 1.47). Conclusions This is the first demonstration of a genetic influence on RLS in ESRD patients with BTBD9 being significantly associated. The extent of the genetic predisposition could vary between different subgroups of RLS in ESRD.

Allelic differences between Europeans and Chinese for CREB1 SNPs and their implications in gene expression regulation, hippocampal structure and function, and bipolar disorder susceptibility

Bipolar disorder (BD) is a polygenic disorder that shares substantial genetic risk factors with major depressive disorder (MDD). Genetic analyses have reported numerous BD susceptibility genes, while some variants, such as single-nucleotide polymorphisms (SNPs) in CACNA1C have been successfully replicated, many others have not and subsequently their effects on the intermediate phenotypes cannot be verified. Here, we studied the MDD-related gene CREB1 in a set of independent BD sample groups of European ancestry (a total of 64 888 subjects) and identified multiple SNPs significantly associated with BD (the most significant being SNP rs6785[A], P=6.32 × 10−5, odds ratio (OR)=1.090). Risk SNPs were then subjected to further analyses in healthy Europeans for intermediate phenotypes of BD, including hippocampal volume, hippocampal function and cognitive performance. Our results showed that the risk SNPs were significantly associated with hippocampal volume and hippocampal function, with the risk alleles showing a decreased hippocampal volume and diminished activation of the left hippocampus, adding further evidence for their involvement in BD susceptibility. We also found the risk SNPs were strongly associated with CREB1 expression in lymphoblastoid cells (P<0.005) and the prefrontal cortex (P<1.0 × 10−6). Remarkably, population genetic analysis indicated that CREB1 displayed striking differences in allele frequencies between continental populations, and the risk alleles were completely absent in East Asian populations. We demonstrated that the regional prevalence of the CREB1 risk alleles in Europeans is likely caused by genetic hitchhiking due to natural selection acting on a nearby gene. Our results suggest that differential population histories due to natural selection on regional populations may lead to genetic heterogeneity of susceptibility to complex diseases, such as BD, and explain inconsistencies in detecting the genetic markers of these diseases among different ethnic populations.