Genome-wide association studies (GWAS), although efficient to detect genes involved in complex diseases, are not designed to measure the real effect of the genes. This is illustrated here by the example of IL2RA in multiple sclerosis (MS). Association between IL2RA and MS is clearly established, although the functional variation is still unknown: the effect of IL2RA might be better described by several SNPs than by a single one. This study investigates whether a pair of SNPs better explains the observed linkage and association data than a single SNP. In total, 522 trio families and 244 affected sib-pairs were typed for 26 IL2RA SNPs. For each SNP and pairs of SNPs, the phased genotypes of patients and controls were compared to determine the SNP set offering the best risk discrimination. Consistency between the genotype risks provided by the retained set and the identical by descent allele sharing in affected sib-pairs was assessed. After controlling for multiple testing, the set of SNPs rs2256774 and rs3118470, provides the best discrimination between the case and control genotype distributions (P-corrected=0.009). The relative risk between the least and most at-risk genotypes is 3.54 with a 95% confidence interval of [2.14–5.94]. Furthermore, the linkage information provided by the allele sharing between affected sibs is consistent with the retained set (P=0.80) but rejects the SNP reported in the literature (P=0.006). Establishing a valid modeling of a disease gene is essential to test its potential interaction with other genes and to reconstruct the pathophysiological pathways.

Dyslexia is a developmental disorder characterised by extensive difficulties in the acquisition of reading or spelling. Genetic influence is estimated at 50–70%. However, the link between genetic variants and phenotypic deficits is largely unknown. Our aim was to investigate a role of genetic variants of FOXP2, a prominent speech and language gene, in dyslexia using imaging genetics. This technique combines functional magnetic resonance imaging (fMRI) and genetics to investigate relevance of genetic variants on brain activation. To our knowledge, this represents the first usage of fMRI-based imaging genetics in dyslexia. In an initial case/control study (n=245) for prioritisation of FOXP2 polymorphisms for later use in imaging genetics, nine SNPs were selected. A non-synonymously coding mutation involved in verbal dyspraxia was also investigated. SNP rs12533005 showed nominally significant association with dyslexia (genotype GG odds ratio recessive model=2.1 (95% confidence interval 1.1–3.9), P=0.016). A correlated SNP was associated with altered expression of FOXP2 in vivo in human hippocampal tissue. Therefore, influence of the rs12533005-G risk variant on brain activity was studied. fMRI revealed a significant main effect for the factor ‘genetic risk' in a temporo-parietal area involved in phonological processing as well as a significant interaction effect between the factors ‘disorder' and ‘genetic risk' in activation of inferior frontal brain areas. Hence, our data may hint at a role of FOXP2 genetic variants in dyslexia-specific brain activation and demonstrate use of imaging genetics in dyslexia research.

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.

Alcohol dependence (AD) is an important contributory factor to the global burden of disease. The etiology of AD involves both environmental and genetic factors, and the disorder has a heritability of around 50%. The aim of the present study was to identify susceptibility genes for AD by performing a genome-wide association study (GWAS). The sample comprised 1,333 male in-patients with severe DSM-IV AD and 2,168 controls. These included 487 patients and 1,358 controls from a previous GWAS study by our group. All individuals were of German descent. Single marker tests and a polygenic score based analysis to assess the combined contribution of multiple markers with small effects were performed. The SNP rs1789891, which is located between the ADH1B and ADH1C genes, achieved genome-wide significance (p=1.27E–8; OR=1.46). Other markers from this region were also associated with AD, and conditional analyses indicated that these made a partially independent contribution. The SNP rs1789891 is in complete linkage disequilibrium with the functional Arg272Gln variant (p=1.24E–7, OR=1.31) of the ADH1C gene, which has been reported to modify the rate of ethanol oxidation to acetaldehyde in vitro. A polygenic score based approach produced a significant result (p=9.66E–9). This is the first GWAS of AD to provide genome-wide significant support for the role of the ADH gene cluster and to suggest a polygenic component to the etiology of AD. The latter result suggests that many more AD susceptibility genes still await identification.

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−5; 0.76 [0.67–0.87]) and of rs7720838 (p = 6.91×10−4; 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−7 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

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.

Background

Given that genome wide association studies (GWAS) of psychiatric disorders have identified only a small number of convincingly associated variants, there is interest in seeking additional evidence for associated variants using tests of gene-gene interaction. Comprehensive pair-wise SNP-SNP interaction analysis is computationally intensive and the penalty for multiple testing is severe given the number of interactions possible. Aiming to minimize these statistical and computational burdens, we have explored approaches to prioritise SNPs for interaction analyses.

Methods

Primary interaction analyses were performed using the Wellcome Trust Case Control Consortium Bipolar Disorder GWAS (1868 cases, 2938 controls). Replication analyses were performed using the Genetic Association Information Network BD dataset (1001 cases, 1033 controls). SNPs were prioritized for interaction analysis that showed evidence for association that surpassed a number of nominally significant thresholds, are within genome-wide significant genes, or are within genes that are functionally related.

Results

For no set of prioritized SNPs did we obtain evidence to support the hypothesis that the selection strategy identified pairs of variants that were enriched for true (statistical) interactions.

Conclusions

SNPs prioritized according to a number of criteria do not have a raised prior probability for significant interaction that is detectable in samples of this size. As is now widely accepted for single SNP analysis, we argue the use of significance levels reflecting only the number of tests performed does not offer an appropriate degree of protection against the potential for GWAS studies to generate an enormous number of false positive interactions.

Dyslexia affects 5–10% of school-aged children and is therefore one of the most common learning disorders. Research on auditory event related potentials (AERP), particularly the mismatch negativity (MMN) component, has revealed anomalies in individuals with dyslexia to speech stimuli. Furthermore, candidate genes for this disorder were found through molecular genetic studies. A current challenge for dyslexia research is to understand the interaction between molecular genetics and brain function, and to promote the identification of relevant endophenotypes for dyslexia. The present study examines MMN, a neurophysiological correlate of speech perception, and its potential as an endophenotype for dyslexia in three groups of children. The first group of children was clinically diagnosed with dyslexia, whereas the second group of children was comprised of their siblings who had average reading and spelling skills and were therefore “unaffected” despite having a genetic risk for dyslexia. The third group consisted of control children who were not related to the other groups and were also unaffected. In total, 225 children were included in the study. All children showed clear MMN activity to/da/−/ba/contrasts that could be separated into three distinct MMN components. Whilst the first two MMN components did not differentiate the groups, the late MMN component (300–700 ms) revealed significant group differences. The mean area of the late MMN was attenuated in both the dyslexic children and their unaffected siblings in comparison to the control children. This finding is indicative of analogous alterations of neurophysiological processes in children with dyslexia and those with a genetic risk for dyslexia, without a manifestation of the disorder. The present results therefore further suggest that the late MMN might be a potential endophenotype for dyslexia.

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 likely contributes to MD. We present data from a genome-wide association study revealing a neuron-specific neutral amino acid transporter (SLC6A15) as a novel 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 novel pathophysiological mechanism for MD that may be accessible to drug targeting.

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.

Inspired by the localization, on 15q21.2 of the CYP19A1 gene in the linkage region of speech and language disorders, and a rare translocation in a dyslexic individual that was brought to our attention, we conducted a series of studies on the properties of CYP19A1 as a candidate gene for dyslexia and related conditions. The aromatase enzyme is a member of the cytochrome P450 super family, and it serves several key functions: it catalyzes the conversion of androgens into estrogens; during early mammalian development it controls the differentiation of specific brain areas (e.g. local estrogen synthesis in the hippocampus regulates synaptic plasticity and axonal growth); it is involved in sexual differentiation of the brain; and in songbirds and teleost fishes, it regulates vocalization. Our results suggest that variations in CYP19A1 are associated with dyslexia as a categorical trait and with quantitative measures of language and speech, such as reading, vocabulary, phonological processing and oral motor skills. Variations near the vicinity of its brain promoter region altered transcription factor binding, suggesting a regulatory role in CYP19A1 expression. CYP19A1 expression in human brain correlated with the expression of dyslexia susceptibility genes such as DYX1C1 and ROBO1. Aromatase-deficient mice displayed increased cortical neuronal density and occasional cortical heterotopias, also observed in Robo1−/− mice and human dyslexic brains, respectively. An aromatase inhibitor reduced dendritic growth in cultured rat neurons. From this broad set of evidence, we propose CYP19A1 as a candidate gene for human cognitive functions implicated in reading, speech and language.

Electronic supplementary material

The online version of this article (doi:10.1007/s10519-012-9532-3) contains supplementary material, which is available to authorized users.

The level of body iron storage and the erythropoietic need for iron are indicated by the serum levels of ferritin and soluble transferrin receptor (sTfR), respectively. A meta-analysis of five genome-wide association studies on sTfR and ferritin revealed novel association to the PCSK7 and TMPRSS6 loci for sTfR and the HFE locus for both parameters. The PCSK7 association was the most significant (rs236918, P = 1.1 × 10E−27) suggesting that proprotein convertase 7, the gene product of PCSK7, may be involved in sTfR generation and/or iron homeostasis. Conditioning the sTfR analyses on transferrin saturation abolished the HFE signal and substantially diminished the TMPRSS6 signal while the PCSK7 association was unaffected, suggesting that the former may be mediated by transferrin saturation whereas the PCSK7-associated effect on sTfR generation appears to be more direct.

Background

Suicidal behaviour can be conceptualised as a continuum from suicidal ideation, to suicidal attempts to completed suicide. In this study we identify genes contributing to suicidal behaviour in the depression study RADIANT.

Methodology/Principal Findings

A quantitative suicidality score was composed of two items from the SCAN interview. In addition, the 251 depression cases with a history of serious suicide attempts were classified to form a discrete trait. The quantitative trait was correlated with younger onset of depression and number of episodes of depression, but not with gender. A genome-wide association study of 2,023 depression cases was performed to identify genes that may contribute to suicidal behaviour. Two Munich depression studies were used as replication cohorts to test the most strongly associated SNPs. No SNP was associated at genome-wide significance level. For the quantitative trait, evidence of association was detected at GFRA1, a receptor for the neurotrophin GDRA (p = 2e-06). For the discrete trait of suicide attempt, SNPs in KIAA1244 and RGS18 attained p-values of <5e-6. None of these SNPs showed evidence for replication in the additional cohorts tested. Candidate gene analysis provided some support for a polymorphism in NTRK2, which was previously associated with suicidality.

Conclusions/Significance

This study provides a genome-wide assessment of possible genetic contribution to suicidal behaviour in depression but indicates a genetic architecture of multiple genes with small effects. Large cohorts will be required to dissect this further.

Motivation: In recent years, numerous genome-wide association studies have been conducted to identify genetic makeup that explains phenotypic differences observed in human population. Analytical tests on single loci are readily available and embedded in common genome analysis software toolset. The search for significant epistasis (gene–gene interactions) still poses as a computational challenge for modern day computing systems, due to the large number of hypotheses that have to be tested.

Results: In this article, we present an approach to epistasis detection by exhaustive testing of all possible SNP pairs. The search strategy based on the Hilbert–Schmidt Independence Criterion can help delineate various forms of statistical dependence between the genetic markers and the phenotype. The actual implementation of this search is done on the highly parallelized architecture available on graphics processing units rendering the completion of the full search feasible within a day.

Availability:The program is available at http://www.mpipsykl.mpg.de/epigpuhsic/.

Contact: tony@mpipsykl.mpg.de

Motivation: High-throughput-sequencing (HTS) technologies are the method of choice for screening the human genome for rare sequence variants causing susceptibility to complex diseases. Unfortunately, preparation of samples for a large number of individuals is still very cost- and labor intensive. Thus, recently, screens for rare sequence variants were carried out in samples of pooled DNA, in which equimolar amounts of DNA from multiple individuals are mixed prior to sequencing with HTS. The resulting sequence data, however, poses a bioinformatics challenge: the discrimination of sequencing errors from real sequence variants present at a low frequency in the DNA pool.

Results: Our method vipR uses data from multiple DNA pools in order to compensate for differences in sequencing error rates along the sequenced region. More precisely, instead of aiming at discriminating sequence variants from sequencing errors, vipR identifies sequence positions that exhibit significantly different minor allele frequencies in at least two DNA pools using the Skellam distribution. The performance of vipR was compared with three other models on data from a targeted resequencing study of the TMEM132D locus in 600 individuals distributed over four DNA pools. Performance of the methods was computed on SNPs that were also genotyped individually using a MALDI-TOF technique. On a set of 82 sequence variants, vipR achieved an average sensitivity of 0.80 at an average specificity of 0.92, thus outperforming the reference methods by at least 0.17 in specificity at comparable sensitivity.

Availability: The code of vipR is freely available via: http://sourceforge.net/projects/htsvipr/

Contact: altmann@mpipsykl.mpg.de

Galanin (GAL) is an estrogen-inducible neuropeptide, highly expressed in brain regions reported to be involved in regulation of mood and anxiety. GAL possibly has a direct modulatory effect on hypothalamic–pituitary–adrenal (HPA)-axis regulation. Recent data from pharmacological and genetic studies indicate a significant function of GAL in stress-related disorders. By using a tag SNP approach covering the locus encoding preprogalanin (PPGAL), earlier findings of female-specific associations of polymorphisms in this locus with panic disorder were expanded to a larger sample of 268 outpatients with anxiety disorders (ADs). Within a larger sample of 541 inpatients with major depressive disorder (MDD), we then tested associations of one PPGAL tag SNP with specific depression symptom clusters and HPA-axis activity assessed by the combined dexamethasone-suppression/CRH-stimulation test both at inpatient admission and discharge (n=298). Gender specificity as well as dependence of the association on levels of circulating estrogens was analyzed. Genotyping revealed high linkage disequilibrium in the promoter area of the PPGAL gene, which includes several estrogen-response elements. Confirming earlier results, rs948854, tagging this promoter region, was associated with more severe anxiety pathology in female AD patients, but not in males. In premenopausal female MDD patients, the same allele of rs948854 was associated with more severe vegetative but not cognitive depressive symptoms at discharge and worse treatment response on antidepressant medication. Furthermore, this allele was associated with higher HPA-axis activity at admission. No significant case–control associations could be observed. However, because of power limitations of both patient samples, small effects cannot be excluded. The reported associations in independent samples of AD and MDD support an estrogen-dependent function of GAL in pathophysiology of anxiety and depression, affecting response to antidepressant treatment.

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.

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.

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (p=5.12 × 10−9, OR 1.23 [1.150-1.324]) in a genome-wide association study of 2,748 migraineurs from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis p-value of 1.60 × 10−11 (OR 1.18 [1.127 – 1.244]). rs1835740 is located between the astrocyte elevated gene 1 (MTDH/AEG-1) and plasma glutamate carboxypeptidase (PGCP). In an expression quantitative trait study in lymphoblastoid cell lines transcript levels of the MTDH/AEG-1 were found to have a significant correlation to rs1835740. Our data establish rs1835740 as the first genetic risk factor for migraine.

Single-nucleotide polymorphisms (SNPs) in the FKBP5, GRIK4, and HTR2A genes have been shown to be associated with response to citalopram treatment in the STAR*D sample, but only associations with FKBP5 have so far been tested in the Munich Antidepressant Response Signature (MARS) project. Response and remission of depressive symptoms after 5 weeks of antidepressant treatment were tested against 82 GRIK4 and 37 HTR2A SNPs. Association analysis was conducted in about 300 depressed patients from the MARS project, 10% of whom had bipolar disorder. The most predictive SNPs from these two genes and rs1360780 in FKBP5 were then genotyped in a total of 387 German depressed in-patients to analyze potential additive and interactive effects of these variants. We could not replicate previous findings of the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study in our sample. Although not statistically significant, the effect for the best GRIK4 SNP of STAR*D (rs1954787, p=0.076, pcorrected=0.98) seemed to be in the same direction. On the other hand, the nominally significant association with the top HTR2A SNPs of STAR*D (rs7997012, allelic, p=0.043, pcorrected=0.62) was with the opposite risk allele. The GRIK4 SNP (rs12800734, genotypic, p=0.0019, pcorrected=0.12) and the HTR2A SNP (rs17288723, genotypic, p=0.0011, pcorrected=0.02), which showed the strongest association with remission in our sample, had not been reported previously. Associations across all genetic markers within the GRIK4 (genotypic, p=0.022) or HTR2A (genotypic, p=0.012) locus using the Fisher's product method (FPM) were also significant. In all 374 patients, the best predictive model included a main effect for GRIK4 rs12800734 and two significant interactions between GRIK4 rs12800734 and FKBP5 rs1360780, and GRIK4 rs12800734 and HTR2A rs17288723. This three SNP model explained 13.1% of the variance for remission after 5 weeks (p=0.00051 for the model). Analyzing a sub-sample of 194 patients, plasma ACTH (p=0.002) and cortisol (p=0.021) responses of rs12800734 GG (GRIK4) carriers, who also showed favorable treatment response, were significantly lower in the second combined dexamethasone (dex)/corticotrophin-releasing hormone (CRH) test before discharge compared with the other two genotype groups. Despite large differences in ethnicity and design compared with the STAR*D study, our results from the MARS study further support both independent and interactive involvement of GRIK4, HTR2A and FKBP5 in antidepressant treatment response.

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.

The QT interval, a measure of cardiac repolarization, predisposes to ventricular arrhythmias and sudden cardiac death (SCD) when prolonged or shortened. A common variant in NOS1AP is known to influence repolarization. We analyze genome-wide data from five population-based cohorts (ARIC, KORA, SardiNIA, GenNOVA and HNR) with a total of 15,842 individuals of European ancestry, to confirm the NOS1AP association and identify nine additional loci at P < 5 × 10−8. Four loci map near the monogenic long-QT syndrome genes KCNQ1, KCNH2, SCN5A and KCNJ2. Two other loci include ATP1B1 and PLN, genes with established electrophysiological function, whereas three map to RNF207, near LITAF and within NDRG4-GINS3-SETD6-CNOT1, respectively, all of which have not previously been implicated in cardiac electrophysiology. These results, together with an accompanying paper from the QTGEN consortium, identify new candidate genes for ventricular arrhythmias and SCD.