Kristel R. van Eijk

Aging effects on DNA methylation modules in human brain and blood tissue

BackgroundSeveral recent studies reported aging effects on DNA methylation levels of individual CpG dinucleotides. But it is not yet known whether aging-related consensus modules, in the form of clusters of correlated CpG markers, can be found that are present in multiple human tissues. Such a module could facilitate the understanding of aging effects on multiple tissues.ResultsWe therefore employed weighted correlation network analysis of 2,442 Illumina DNA methylation arrays from brain and blood tissues, which enabled the identification of an age-related co-methylation module. Module preservation analysis confirmed that this module can also be found in diverse independent data sets. Biological evaluation showed that module membership is associated with Polycomb group target occupancy counts, CpG island status and autosomal chromosome location. Functional enrichment analysis revealed that the aging-related consensus module comprises genes that are involved in nervous system development, neuron differentiation and neurogenesis, and that it contains promoter CpGs of genes known to be down-regulated in early Alzheimers disease. A comparison with a standard, non-module based meta-analysis revealed that selecting CpGs based on module membership leads to significantly increased gene ontology enrichment, thus demonstrating that studying aging effects via consensus network analysis enhances the biological insights gained.ConclusionsOverall, our analysis revealed a robustly defined age-related co-methylation module that is present in multiple human tissues, including blood and brain. We conclude that blood is a promising surrogate for brain tissue when studying the effects of age on DNA methylation profiles.

Using genome-wide pathway analysis to unravel the etiology of complex diseases

Several genome‐wide association studies (GWAS) have been published on various complex diseases. Although, new loci are found to be associated with these diseases, still only very little of the genetic risk for these diseases can be explained. As GWAS are still underpowered to find small main effects, and gene‐gene interactions are likely to play a role, the data might currently not be analyzed to its full potential. In this study, we evaluated alternative methods to study GWAS data. Instead of focusing on the single nucleotide polymorphisms (SNPs) with the highest statistical significance, we took advantage of prior biological information and tried to detect overrepresented pathways in the GWAS data. We evaluated whether pathway classification analysis can help prioritize the biological pathways most likely to be involved in the disease etiology. In this study, we present the various benefits and limitations of pathway‐classification tools in analyzing GWAS data. We show multiple differences in outcome between pathway tools analyzing the same dataset. Furthermore, analyzing randomly selected SNPs always results in significantly overrepresented pathways, large pathways have a higher chance of becoming statistically significant and the bioinformatics tools used in this study are biased toward detecting well‐defined pathways. As an example, we analyzed data from two GWAS on type 2 diabetes (T2D): the Diabetes Genetics Initiative (DGI) and the Wellcome Trust Case Control Consortium (WTCCC). Occasionally the results from the DGI and the WTCCC GWAS showed concordance in overrepresented pathways, but discordance in the corresponding genes. Thus, incorporating gene networks and pathway classification tools into the analysis can point toward significantly overrepresented molecular pathways, which cannot be picked up using traditional single‐locus analyses. However, the limitations discussed in this study, need to be addressed before these methods can be widely used. Genet. Epidemiol. 33:419–431, 2009.

Genetic analysis of DNA methylation and gene expression levels in whole blood of healthy human subjects

BackgroundThe predominant model for regulation of gene expression through DNA methylation is an inverse association in which increased methylation results in decreased gene expression levels. However, recent studies suggest that the relationship between genetic variation, DNA methylation and expression is more complex.ResultsSystems genetic approaches for examining relationships between gene expression and methylation array data were used to find both negative and positive associations between these levels. A weighted correlation network analysis revealed that i) both transcriptome and methylome are organized in modules, ii) co-expression modules are generally not preserved in the methylation data and vice-versa, and iii) highly significant correlations exist between co-expression and co-methylation modules, suggesting the existence of factors that affect expression and methylation of different modules (i.e., trans effects at the level of modules). We observed that methylation probes associated with expression in cis were more likely to be located outside CpG islands, whereas specificity for CpG island shores was present when methylation, associated with expression, was under local genetic control. A structural equation model based analysis found strong support in particular for a traditional causal model in which gene expression is regulated by genetic variation via DNA methylation instead of gene expression affecting DNA methylation levels.ConclusionsOur results provide new insights into the complex mechanisms between genetic markers, epigenetic mechanisms and gene expression. We find strong support for the classical model of genetic variants regulating methylation, which in turn regulates gene expression. Moreover we show that, although the methylation and expression modules differ, they are highly correlated.

NEK1 variants confer susceptibility to amyotrophic lateral sclerosis

To identify genetic factors contributing to amyotrophic lateral sclerosis (ALS), we conducted whole-exome analyses of 1,022 index familial ALS (FALS) cases and 7,315 controls. In a new screening strategy, we performed gene-burden analyses trained with established ALS genes and identified a significant association between loss-of-function (LOF) NEK1 variants and FALS risk. Independently, autozygosity mapping for an isolated community in the Netherlands identified a NEK1 p.Arg261His variant as a candidate risk factor. Replication analyses of sporadic ALS (SALS) cases and independent control cohorts confirmed significant disease association for both p.Arg261His (10,589 samples analyzed) and NEK1 LOF variants (3,362 samples analyzed). In total, we observed NEK1 risk variants in nearly 3% of ALS cases. NEK1 has been linked to several cellular functions, including cilia formation, DNA-damage response, microtubule stability, neuronal morphology and axonal polarity. Our results provide new and important insights into ALS etiopathogenesis and genetic etiology.

Expression QTL analysis of top loci from GWAS meta-analysis highlights additional schizophrenia candidate genes

There is genetic evidence that schizophrenia is a polygenic disorder with a large number of loci of small effect on disease susceptibility. Genome-wide association studies (GWASs) of schizophrenia have had limited success, with the best finding at the MHC locus at chromosome 6p. A recent effort of the Psychiatric GWAS consortium (PGC) yielded five novel loci for schizophrenia. In this study, we aim to highlight additional schizophrenia susceptibility loci from the PGC study by combining the top association findings from the discovery stage (9394 schizophrenia cases and 12 462 controls) with expression QTLs (eQTLs) and differential gene expression in whole blood of schizophrenia patients and controls. We examined the 6192 single-nucleotide polymorphisms (SNPs) with significance threshold at P<0.001. eQTLs were calculated for these SNPs in a sample of healthy controls (n=437). The transcripts significantly regulated by the top SNPs from the GWAS meta-analysis were subsequently tested for differential expression in an independent set of schizophrenia cases and controls (n=202). After correction for multiple testing, the eQTL analysis yielded 40 significant cis-acting effects of the SNPs. Seven of these transcripts show differential expression between cases and controls. Of these, the effect of three genes (RNF5, TRIM26 and HLA-DRB3) coincided with the direction expected from meta-analysis findings and were all located within the MHC region. Our results identify new genes of interest and highlight again the involvement of the MHC region in schizophrenia susceptibility.

The effect of childhood maltreatment and cannabis use on adult psychotic symptoms is modified by the COMT Val158Met polymorphism

BACKGROUND Cannabis use and childhood maltreatment are independent risk factors for the development of psychotic symptoms. These factors have been found to interact in some but not all studies. One of the reasons may be that childhood maltreatment and cannabis primarily induce psychotic symptoms in genetically susceptible individuals. In this context, an extensively studied psychosis vulnerability gene is catechol-methyl-transferase (COMT). Therefore, we aimed to examine whether the COMT Val(158)Met polymorphism (rs4680) moderates the interaction between childhood maltreatment and cannabis use on psychotic symptoms in the general population. METHOD The discovery sample consisted of 918 individuals from a cross-sectional study. For replication we used an independent sample of 339 individuals from the general population. RESULTS A significant three-way interaction was found between childhood maltreatment, cannabis use, and the COMT genotype (rs4680) in the discovery sample (P=0.006). Val-homozygous individuals displayed increased psychotic experiences after exposure to both cannabis use and childhood maltreatment compared to Met-heterozygous and Met-homozygous individuals. Supportive evidence was found in the replication sample with similar effect and direction even though the results did not reach statistical significance (P=0.25). CONCLUSIONS These findings suggest that a functional polymorphism in the COMT gene may moderate the interaction between childhood maltreatment and cannabis use on psychotic experiences in the general population. In conclusion, the COMT Val(158)Met polymorphism may constitute a genetic risk factor for psychotic symptoms in the context of combined exposure to childhood maltreatment and cannabis use.

Identification of schizophrenia-associated loci by combining DNA methylation and gene expression data from whole blood.

Emerging evidence suggests that schizophrenia (SZ) susceptibility involves variation at genetic, epigenetic and transcriptome levels. We describe an integrated approach that leverages DNA methylation and gene expression data to prioritize genetic variation involved in disease. DNA methylation levels were obtained from whole blood of 260 SZ patients and 250 unaffected controls of which a subset with gene expression data was available. By assessing DNA methylation and gene expression in cases and controls, we identified 432 CpG sites with differential methylation levels that are associated with differential gene expression. We hypothesized that genetic factors involved in these methylation levels may be associated with the genetic risk of SZ susceptibility. To test this hypothesis, we used results from the Psychiatric Genomics Consortium SZ genome-wide association study (GWAS). We observe an enrichment of SZ-associated SNPs in the mQTLs of which the associated CpG site is also correlated with differential gene expression in SZ. While this enrichment was already apparent when using nominal significant thresholds, enrichment was even more pronounced when applying more stringent significance levels. One locus, previously identified as susceptibility locus in a SZ GWAS, involves SNP rs11191514:C>T, which regulates DNA methylation of calcium homeostasis modulator 1 that is also associated with differential gene expression in patients. Overall, our results suggest that epigenetic variation plays an important role in SZ susceptibility and that the integration of analyses of genetic, epigenetic and gene expression profiles may be a biologically meaningful approach for identifying disease susceptibility loci, even when using whole blood data in studies of brain-related disorders.

A Common Variant in ERBB4 Regulates GABA Concentrations in Human Cerebrospinal Fluid

The neuregulin 1 (NRG1) receptor ErbB4 is involved in the development of cortical inhibitory GABAergic circuits and NRG1-ErbB4 signaling has been implicated in schizophrenia (SCZ). A magnetic resonance spectroscopy (1H-MRS) study has demonstrated that a single-nucleotide polymorphism in ERBB4, rs7598440, influences human cortical GABA concentrations. Other work has highlighted the significant impact of this genetic variant on expression of ERBB4 in the hippocampus and dorsolateral prefrontal cortex in human post mortem tissue. Our aim was to examine the association of rs7598440 with cerebrospinal fluid (CSF) GABA levels in healthy volunteers (n=155). We detected a significant dose-dependent association of the rs7598440 genotype with CSF GABA levels (G-allele standardized β=−0.23; 95% CIs: −0.39 to −0.07; P=0.0066). GABA concentrations were highest in A homozygous, intermediate in heterozygous, and lowest in G homozygous subjects. When excluding subjects on psychotropic medication (three subjects using antidepressants), the results did not change (G-allele standardized β=−0.23; 95% CIs: −0.40 to −0.07; P=0.0051). The explained variance in CSF GABA by rs7598440 in our model is 5.2% (P=0.004). The directionality of our findings agrees with the aforementioned 1H-MRS and gene expression studies. Our observation therefore strengthens the evidence that the A-allele of rs7598440 in ERBB4 is associated with increased GABA concentrations in the human central nervous system (CNS). To our knowledge, our finding constitutes the first confirmation that CSF can be used to study genotype–phenotype correlations of GABA levels in the CNS. Such quantitative genetic analyses may be extrapolated to other CSF constituents relevant to SCZ in future studies.

DRD2 Schizophrenia-Risk Allele Is Associated With Impaired Striatal Functioning in Unaffected Siblings of Schizophrenia Patients

A recent Genome-Wide Association Study showed that the rs2514218 single nucleotide polymorphism (SNP) in close proximity to dopamine receptor D2 is strongly associated with schizophrenia. Further, an in silico experiment showed that rs2514218 has a cis expression quantitative trait locus effect in the basal ganglia. To date, however, the functional consequence of this SNP is unknown. Here, we used functional Magnetic resonance imaging to investigate the impact of this risk allele on striatal activation during proactive and reactive response inhibition in 45 unaffected siblings of schizophrenia patients. We included siblings to circumvent the illness specific confounds affecting striatal functioning independent from gene effects. Behavioral analyses revealed no differences between the carriers (n= 21) and noncarriers (n= 24). Risk allele carriers showed a diminished striatal response to increasing proactive inhibitory control demands, whereas overall level of striatal activation in carriers was elevated compared to noncarriers. Finally, risk allele carriers showed a blunted striatal response during successful reactive inhibition compared to the noncarriers. These data are consistent with earlier reports showing similar deficits in schizophrenia patients, and point to a failure to flexibly engage the striatum in response to contextual cues. This is the first study to demonstrate an association between impaired striatal functioning and the rs2514218 polymorphism. We take our findings to indicate that striatal functioning is impaired in carriers of the DRD2 risk allele, likely due to dopamine dysregulation at the DRD2 location.

Serotonin 2B receptor slows disease progression and prevents degeneration of spinal cord mononuclear phagocytes in amyotrophic lateral sclerosis

Microglia are the resident mononuclear phagocytes of the central nervous system and have been implicated in the pathogenesis of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS). During neurodegeneration, microglial activation is accompanied by infiltration of circulating monocytes, leading to production of multiple inflammatory mediators in the spinal cord. Degenerative alterations in mononuclear phagocytes are commonly observed during neurodegenerative diseases, yet little is known concerning the mechanisms leading to their degeneration, or the consequences on disease progression. Here we observed that the serotonin 2B receptor (5-HT2B), a serotonin receptor expressed in microglia, is upregulated in the spinal cord of three different transgenic mouse models of ALS. In mutant SOD1 mice, this upregulation was restricted to cells positive for CD11b, a marker of mononuclear phagocytes. Ablation of 5-HT2B receptor in transgenic ALS mice expressing mutant SOD1 resulted in increased degeneration of mononuclear phagocytes, as evidenced by fragmentation of Iba1-positive cellular processes. This was accompanied by decreased expression of key neuroinflammatory genes but also loss of expression of homeostatic microglial genes. Importantly, the dramatic effect of 5-HT2B receptor ablation on mononuclear phagocytes was associated with acceleration of disease progression. To determine the translational relevance of these results, we studied polymorphisms in the human HTR2B gene, which encodes the 5-HT2B receptor, in a large cohort of ALS patients. In this cohort, the C allele of SNP rs10199752 in HTR2B was associated with longer survival. Moreover, patients carrying one copy of the C allele of SNP rs10199752 showed increased 5-HT2B mRNA in spinal cord and displayed less pronounced degeneration of Iba1 positive cells than patients carrying two copies of the more common A allele. Thus, the 5-HT2B receptor limits degeneration of spinal cord mononuclear phagocytes, most likely microglia, and slows disease progression in ALS. Targeting this receptor might be therapeutically useful.