26
|
Malki K, Mineur YS, Tosto MG, Campbell J, Karia P, Jumabhoy I, Sluyter F, Crusio WE, Schalkwyk LC. Pervasive and opposing effects of Unpredictable Chronic Mild Stress (UCMS) on hippocampal gene expression in BALB/cJ and C57BL/6J mouse strains. BMC Genomics 2015; 16:262. [PMID: 25879669 PMCID: PMC4412144 DOI: 10.1186/s12864-015-1431-6] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2014] [Accepted: 03/03/2015] [Indexed: 01/03/2023] Open
Abstract
Background BALB/cJ is a strain susceptible to stress and extremely susceptible to a defective hedonic impact in response to chronic stressors. The strain offers much promise as an animal model for the study of stress related disorders. We present a comparative hippocampal gene expression study on the effects of unpredictable chronic mild stress on BALB/cJ and C57BL/6J mice. Affymetrix MOE 430 was used to measure hippocampal gene expression from 16 animals of two different strains (BALB/cJ and C57BL/6J) of both sexes and subjected to either unpredictable chronic mild stress (UCMS) or no stress. Differences were statistically evaluated through supervised and unsupervised linear modelling and using Weighted Gene Coexpression Network Analysis (WGCNA). In order to gain further understanding into mechanisms related to stress response, we cross-validated our results with a parallel study from the GENDEP project using WGCNA in a meta-analysis design. Results The effects of UCMS are visible through Principal Component Analysis which highlights the stress sensitivity of the BALB/cJ strain. A number of genes and gene networks related to stress response were uncovered including the Creb1 gene. WGCNA and pathway analysis revealed a gene network centered on Nfkb1. Results from the meta-analysis revealed a highly significant gene pathway centred on the Ubiquitin C (Ubc) gene. All pathways uncovered are associated with inflammation and immune response. Conclusions The study investigated the molecular mechanisms underlying the response to adverse environment in an animal model using a GxE design. Stress-related differences were visible at the genomic level through PCA analysis highlighting the high sensitivity of BALB/cJ animals to environmental stressors. Several candidate genes and gene networks reported are associated with inflammation and neurogenesis and could serve to inform candidate gene selection in human studies and provide additional insight into the pathology of Major Depressive Disorder. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1431-6) contains supplementary material, which is available to authorized users.
Collapse
|
27
|
Elliott G, Hong C, Xing X, Zhou X, Li D, Coarfa C, Bell RJ, Maire CL, Ligon KL, Sigaroudinia M, Gascard P, Tlsty TD, Harris RA, Schalkwyk LC, Bilenky M, Mill J, Farnham PJ, Kellis M, Marra MA, Milosavljevic A, Hirst M, Stormo GD, Wang T, Costello JF. Intermediate DNA methylation is a conserved signature of genome regulation. Nat Commun 2015; 6:6363. [PMID: 25691127 PMCID: PMC4333717 DOI: 10.1038/ncomms7363] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2014] [Accepted: 01/23/2015] [Indexed: 01/06/2023] Open
Abstract
The role of intermediate methylation states in DNA is unclear. Here, to comprehensively identify regions of intermediate methylation and their quantitative relationship with gene activity, we apply integrative and comparative epigenomics to 25 human primary cell and tissue samples. We report 18,452 intermediate methylation regions located near 36% of genes and enriched at enhancers, exons and DNase I hypersensitivity sites. Intermediate methylation regions average 57% methylation, are predominantly allele-independent and are conserved across individuals and between mouse and human, suggesting a conserved function. These regions have an intermediate level of active chromatin marks and their associated genes have intermediate transcriptional activity. Exonic intermediate methylation correlates with exon inclusion at a level between that of fully methylated and unmethylated exons, highlighting gene context-dependent functions. We conclude that intermediate DNA methylation is a conserved signature of gene regulation and exon usage.
Collapse
|
28
|
Spiers H, Hannon E, Schalkwyk LC, Smith R, Wong CCY, O'Donovan MC, Bray NJ, Mill J. Methylomic trajectories across human fetal brain development. Genome Res 2015; 25:338-52. [PMID: 25650246 PMCID: PMC4352878 DOI: 10.1101/gr.180273.114] [Citation(s) in RCA: 195] [Impact Index Per Article: 21.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
Abstract
Epigenetic processes play a key role in orchestrating transcriptional regulation during development. The importance of DNA methylation in fetal brain development is highlighted by the dynamic expression of de novo DNA methyltransferases during the perinatal period and neurodevelopmental deficits associated with mutations in the methyl-CpG binding protein 2 (MECP2) gene. However, our knowledge about the temporal changes to the epigenome during fetal brain development has, to date, been limited. We quantified genome-wide patterns of DNA methylation at ∼ 400,000 sites in 179 human fetal brain samples (100 male, 79 female) spanning 23 to 184 d post-conception. We identified highly significant changes in DNA methylation across fetal brain development at >7% of sites, with an enrichment of loci becoming hypomethylated with fetal age. Sites associated with developmental changes in DNA methylation during fetal brain development were significantly underrepresented in promoter regulatory regions but significantly overrepresented in regions flanking CpG islands (shores and shelves) and gene bodies. Highly significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small number of regions showing sex-specific DNA methylation trajectories across brain development. Weighted gene comethylation network analysis (WGCNA) revealed discrete modules of comethylated loci associated with fetal age that are significantly enriched for genes involved in neurodevelopmental processes. This is, to our knowledge, the most extensive study of DNA methylation across human fetal brain development to date, confirming the prenatal period as a time of considerable epigenomic plasticity.
Collapse
|
29
|
De Jager PL, Srivastava G, Lunnon K, Burgess J, Schalkwyk LC, Yu L, Eaton ML, Keenan BT, Ernst J, McCabe C, Tang A, Raj T, Replogle J, Brodeur W, Gabriel S, Chai HS, Younkin C, Younkin SG, Zou F, Szyf M, Epstein CB, Schneider JA, Bernstein BE, Meissner A, Ertekin-Taner N, Chibnik LB, Kellis M, Mill J, Bennett DA. Alzheimer's disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci. Nat Neurosci 2014; 17:1156-63. [PMID: 25129075 PMCID: PMC4292795 DOI: 10.1038/nn.3786] [Citation(s) in RCA: 646] [Impact Index Per Article: 64.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2014] [Accepted: 07/16/2014] [Indexed: 02/07/2023]
Abstract
We used a collection of 708 prospectively collected autopsied brains to assess the methylation state of the brain's DNA in relation to Alzheimer's disease (AD). We found that the level of methylation at 71 of the 415,848 interrogated CpGs was significantly associated with the burden of AD pathology, including CpGs in the ABCA7 and BIN1 regions, which harbor known AD susceptibility variants. We validated 11 of the differentially methylated regions in an independent set of 117 subjects. Furthermore, we functionally validated these CpG associations and identified the nearby genes whose RNA expression was altered in AD: ANK1, CDH23, DIP2A, RHBDF2, RPL13, SERPINF1 and SERPINF2. Our analyses suggest that these DNA methylation changes may have a role in the onset of AD given that we observed them in presymptomatic subjects and that six of the validated genes connect to a known AD susceptibility gene network.
Collapse
|
30
|
Malki K, Tosto MG, Jumabhoy I, Lourdusamy A, Sluyter F, Craig I, Uher R, McGuffin P, Schalkwyk LC. Integrative mouse and human mRNA studies using WGCNA nominates novel candidate genes involved in the pathogenesis of major depressive disorder. Pharmacogenomics 2014; 14:1979-90. [PMID: 24279853 DOI: 10.2217/pgs.13.154] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
AIM This study aims to identify novel genes associated with major depressive disorder and pharmacological treatment response using animal and human mRNA studies. MATERIALS & METHODS Weighted gene coexpression network analysis was used to uncover genes associated with stress factors in mice and to inform mRNA probe set selection in a post-mortem study of depression. RESULTS A total of 171 genes were found to be differentially regulated in response to both early and late stress protocols in a mouse study. Ten human genes, orthologous to mouse genes differentially expressed by stress, were also found to be dysregulated in depressed cases in a human post-mortem brain study from the Stanley Foundation Brain Collection. CONCLUSION Several novel genes associated with depression were uncovered, including NOVA1 and USP9X. Moreover, we found further evidence in support of hippocampal neurogenesis and peripheral inflammation in major depressive disorder.
Collapse
|
31
|
Powell TR, Powell-Smith G, Haddley K, Mcguffin P, Quinn J, Schalkwyk LC, Farmer AE, D'Souza UM. Mood-stabilizers differentially affect housekeeping gene expression in human cells. Int J Methods Psychiatr Res 2014; 23:279-88. [PMID: 24677680 PMCID: PMC6878232 DOI: 10.1002/mpr.1435] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/12/2013] [Revised: 06/12/2013] [Accepted: 07/15/2013] [Indexed: 11/08/2022] Open
Abstract
Recent studies have revealed that antidepressants affect the expression of constitutively expressed "housekeeping genes" commonly used as normalizing reference genes in quantitative polymerase chain reaction (qPCR) experiments. There has yet to be an investigation however on the effects of mood-stabilizers on housekeeping gene stability. The current study utilized lymphoblastoid cell lines (LCLs) derived from patients with mood disorders to investigate the effects of a range of doses of lithium (0, 1, 2 and 5 mM) and sodium valproate (0, 0.06, 0.03 and 0.6 mM) on the stability of 12 housekeeping genes. RNA was extracted from LCLs and qPCR was used to generate cycle threshold (Ct ) values which were input into RefFinder analyses. The study revealed drug-specific effects on housekeeping gene stability. The most stable housekeeping genes in LCLs treated: acutely with sodium valproate were ACTB and RPL13A; acutely with lithium were GAPDH and ATP5B; chronically with lithium were ATP5B and CYC1. The stability of GAPDH and B2M were particularly affected by duration of lithium treatment. The study adds to a growing literature that the selection of appropriate housekeeping genes is important for the accurate normalization of target gene expression in experiments investigating the molecular effects of mood disorder pharmacotherapies.
Collapse
|
32
|
Malki K, Keers R, Tosto MG, Lourdusamy A, Carboni L, Domenici E, Uher R, McGuffin P, Schalkwyk LC. The endogenous and reactive depression subtypes revisited: integrative animal and human studies implicate multiple distinct molecular mechanisms underlying major depressive disorder. BMC Med 2014; 12:73. [PMID: 24886127 PMCID: PMC4046519 DOI: 10.1186/1741-7015-12-73] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/13/2013] [Accepted: 04/10/2014] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND Traditional diagnoses of major depressive disorder (MDD) suggested that the presence or absence of stress prior to onset results in either 'reactive' or 'endogenous' subtypes of the disorder, respectively. Several lines of research suggest that the biological underpinnings of 'reactive' or 'endogenous' subtypes may also differ, resulting in differential response to treatment. We investigated this hypothesis by comparing the gene-expression profiles of three animal models of 'reactive' and 'endogenous' depression. We then translated these findings to clinical samples using a human post-mortem mRNA study. METHODS Affymetrix mouse whole-genome oligonucleotide arrays were used to measure gene expression from hippocampal tissues of 144 mice from the Genome-based Therapeutic Drugs for Depression (GENDEP) project. The study used four inbred mouse strains and two depressogenic 'stress' protocols (maternal separation and Unpredictable Chronic Mild Stress) to model 'reactive' depression. Stress-related mRNA differences in mouse were compared with a parallel mRNA study using Flinders Sensitive and Resistant rat lines as a model of 'endogenous' depression. Convergent genes differentially expressed across the animal studies were used to inform candidate gene selection in a human mRNA post-mortem case control study from the Stanley Brain Consortium. RESULTS In the mouse 'reactive' model, the expression of 350 genes changed in response to early stresses and 370 in response to late stresses. A minimal genetic overlap (less than 8.8%) was detected in response to both stress protocols, but 30% of these genes (21) were also differentially regulated in the 'endogenous' rat study. This overlap is significantly greater than expected by chance. The VAMP-2 gene, differentially expressed across the rodent studies, was also significantly altered in the human study after correcting for multiple testing. CONCLUSIONS Our results suggest that 'endogenous' and 'reactive' subtypes of depression are associated with largely distinct changes in gene-expression. However, they also suggest that the molecular signature of 'reactive' depression caused by early stressors differs considerably from that of 'reactive' depression caused by late stressors. A small set of genes was consistently dysregulated across each paradigm and in post-mortem brain tissue of depressed patients suggesting a final common pathway to the disorder. These genes included the VAMP-2 gene, which has previously been associated with Axis-I disorders including MDD, bipolar depression, schizophrenia and with antidepressant treatment response. We also discuss the implications of our findings for disease classification, personalized medicine and case-control studies of MDD.
Collapse
|
33
|
Smith RG, Fernandes C, Kember R, Schalkwyk LC, Buxbaum J, Reichenberg A, Mill J. Transcriptomic changes in the frontal cortex associated with paternal age. Mol Autism 2014; 5:24. [PMID: 24655730 PMCID: PMC3998024 DOI: 10.1186/2040-2392-5-24] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2013] [Accepted: 03/10/2014] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND Advanced paternal age is robustly associated with several human neuropsychiatric disorders, particularly autism. The precise mechanism(s) mediating the paternal age effect are not known, but they are thought to involve the accumulation of de novo (epi)genomic alterations. In this study we investigate differences in the frontal cortex transcriptome in a mouse model of advanced paternal age. FINDINGS Transcriptomic profiling was undertaken for medial prefrontal cortex tissue dissected from the male offspring of young fathers (2 month old, 4 sires, n = 16 offspring) and old fathers (10 month old, 6 sires, n = 16 offspring) in a mouse model of advancing paternal age. We found a number of differentially expressed genes in the offspring of older fathers, many previously implicated in the aetiology of autism. Pathway analysis highlighted significant enrichment for changes in functional networks involved in inflammation and inflammatory disease, which are also implicated in autism. CONCLUSIONS We observed widespread alterations to the transcriptome associated with advanced paternal age with an enrichment of genes associated with inflammation, an interesting observation given previous evidence linking the immune system to several neuropsychiatric disorders including autism.
Collapse
|
34
|
Powell TR, McGuffin P, D'Souza UM, Cohen-Woods S, Hosang GM, Martin C, Matthews K, Day RK, Farmer AE, Tansey KE, Schalkwyk LC. Putative transcriptomic biomarkers in the inflammatory cytokine pathway differentiate major depressive disorder patients from control subjects and bipolar disorder patients. PLoS One 2014; 9:e91076. [PMID: 24618828 PMCID: PMC3949789 DOI: 10.1371/journal.pone.0091076] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2013] [Accepted: 02/06/2014] [Indexed: 12/12/2022] Open
Abstract
Mood disorders consist of two etiologically related, but distinctly treated illnesses, major depressive disorder (MDD) and bipolar disorder (BPD). These disorders share similarities in their clinical presentation, and thus show high rates of misdiagnosis. Recent research has revealed significant transcriptional differences within the inflammatory cytokine pathway between MDD patients and controls, and between BPD patients and controls, suggesting this pathway may possess important biomarker properties. This exploratory study attempts to identify disorder-specific transcriptional biomarkers within the inflammatory cytokine pathway, which can distinguish between control subjects, MDD patients and BPD patients. This is achieved using RNA extracted from subject blood and applying synthesized complementary DNA to quantitative PCR arrays containing primers for 87 inflammation-related genes. Initially, we use ANOVA to test for transcriptional differences in a 'discovery cohort' (total n = 90) and then we use t-tests to assess the reliability of any identified transcriptional differences in a 'validation cohort' (total n = 35). The two most robust and reliable biomarkers identified across both the discovery and validation cohort were Chemokine (C-C motif) ligand 24 (CCL24) which was consistently transcribed higher amongst MDD patients relative to controls and BPD patients, and C-C chemokine receptor type 6 (CCR6) which was consistently more lowly transcribed amongst MDD patients relative to controls. Results detailed here provide preliminary evidence that transcriptional measures within inflammation-related genes might be useful in aiding clinical diagnostic decision-making processes. Future research should aim to replicate findings detailed in this exploratory study in a larger medication-free sample and examine whether identified biomarkers could be used prospectively to aid clinical diagnosis.
Collapse
|
35
|
Martin C, Tansey KE, Schalkwyk LC, Powell TR. The inflammatory cytokines: molecular biomarkers for major depressive disorder? Biomark Med 2014; 9:169-80. [PMID: 24524646 DOI: 10.2217/bmm.14.29] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Cytokines are pleotropic cell signaling proteins that, in addition to their role as inflammatory mediators, also affect neurotransmitter systems, brain functionality and mood. Here we explore the potential utility of cytokine biomarkers for major depressive disorder. Specifically, we explore how genetic, transcriptomic and proteomic information relating to the cytokines might act as biomarkers, aiding clinical diagnosis and treatment selection processes. We advise future studies to investigate whether cytokine biomarkers might differentiate major depressive disorder patients from other patient groups with overlapping clinical characteristics. Furthermore, we invite future pharmacogenetic studies to investigate whether early antidepressant-induced changes to cytokine mRNA or protein levels precede behavioral changes and act as longer-term predictors of clinical antidepressant response.
Collapse
|
36
|
Tarantino LM, Schalkwyk LC. Introduction to mammalian genome special issue: genetics of behavior. Mamm Genome 2014. [DOI: 10.1007/s00335-013-9496-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
37
|
Powell TR, Tansey KE, Breen G, Farmer AE, Craig IW, Uher R, McGuffin P, D'Souza UM, Schalkwyk LC. ATP-binding cassette sub-family F member 1 (ABCF1) is identified as a putative therapeutic target of escitalopram in the inflammatory cytokine pathway. J Psychopharmacol 2013; 27:609-15. [PMID: 23719290 DOI: 10.1177/0269881113490329] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
The inflammatory cytokine pathway may be a potential therapeutic target for major depressive disorder (MDD). Previous reports suggest that antidepressants have anti-inflammatory properties and can cause a reduction in proinflammatory cytokines. Recent evidence suggests this might be mediated at the level of the transcriptome. The current study investigated the transcription of 86 genes in the inflammatory cytokine pathway both at baseline and after eight weeks of escitalopram treatment in MDD patients who were either clinical responders (n=25) or non-responders (n=21), using a subset of samples in the Genome-Based Therapeutic Drugs for Depression project (GENDEP). Changes in expression between baseline and eight weeks of treatment were assessed using two-tailed t-tests. To establish if any significant expression changes related to clinical response, the magnitude of the relative expression change between baseline and eight weeks of treatment was established and binary logistic regressions were used to compare differences between responders and non-responders. ATP-binding cassette sub-family F member 1 (ABCF1), a translational regulator of the inflammatory cytokine pathway showed a significant increase in expression after escitalopram treatment which was significantly greater in responders compared to non-responders, suggesting that ABCF1 may play a role in mediating antidepressant response.
Collapse
|
38
|
Smith RG, Reichenberg A, Kember RL, Buxbaum JD, Schalkwyk LC, Fernandes C, Mill J. Advanced paternal age is associated with altered DNA methylation at brain-expressed imprinted loci in inbred mice: implications for neuropsychiatric disease. Mol Psychiatry 2013; 18:635-6. [PMID: 22733127 DOI: 10.1038/mp.2012.88] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
39
|
Pidsley R, Y Wong CC, Volta M, Lunnon K, Mill J, Schalkwyk LC. A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics 2013; 14:293. [PMID: 23631413 PMCID: PMC3769145 DOI: 10.1186/1471-2164-14-293] [Citation(s) in RCA: 705] [Impact Index Per Article: 64.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2013] [Accepted: 04/04/2013] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND As the most stable and experimentally accessible epigenetic mark, DNA methylation is of great interest to the research community. The landscape of DNA methylation across tissues, through development and in disease pathogenesis is not yet well characterized. Thus there is a need for rapid and cost effective methods for assessing genome-wide levels of DNA methylation. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a very useful addition to the available methods for DNA methylation analysis but its complex design, incorporating two different assay methods, requires careful consideration. Accordingly, several normalization schemes have been published. We have taken advantage of known DNA methylation patterns associated with genomic imprinting and X-chromosome inactivation (XCI), in addition to the performance of SNP genotyping assays present on the array, to derive three independent metrics which we use to test alternative schemes of correction and normalization. These metrics also have potential utility as quality scores for datasets. RESULTS The standard index of DNA methylation at any specific CpG site is β = M/(M + U + 100) where M and U are methylated and unmethylated signal intensities, respectively. Betas (βs) calculated from raw signal intensities (the default GenomeStudio behavior) perform well, but using 11 methylomic datasets we demonstrate that quantile normalization methods produce marked improvement, even in highly consistent data, by all three metrics. The commonly used procedure of normalizing betas is inferior to the separate normalization of M and U, and it is also advantageous to normalize Type I and Type II assays separately. More elaborate manipulation of quantiles proves to be counterproductive. CONCLUSIONS Careful selection of preprocessing steps can minimize variance and thus improve statistical power, especially for the detection of the small absolute DNA methylation changes likely associated with complex disease phenotypes. For the convenience of the research community we have created a user-friendly R software package called wateRmelon, downloadable from bioConductor, compatible with the existing methylumi, minfi and IMA packages, that allows others to utilize the same normalization methods and data quality tests on 450K data.
Collapse
|
40
|
Lill M, Kõks S, Soomets U, Schalkwyk LC, Fernandes C, Lutsar I, Taba P. Peripheral blood RNA gene expression profiling in patients with bacterial meningitis. Front Neurosci 2013; 7:33. [PMID: 23515576 PMCID: PMC3600829 DOI: 10.3389/fnins.2013.00033] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2013] [Accepted: 02/26/2013] [Indexed: 11/13/2022] Open
Abstract
OBJECTIVES The aim of present study was to find genetic pathways activated during infection with bacterial meningitis (BM) and potentially influencing the course of the infection using genome-wide RNA expression profiling combined with pathway analysis and functional annotation of the differential transcription. METHODS We analyzed 21 patients with BM hospitalized in 2008. The control group consisted of 18 healthy subjects. The RNA was extracted from whole blood, globin mRNA was depleted and gene expression profiling was performed using GeneChip Human Gene 1.0 ST Arrays which can assess the transcription of 28,869 genes. Gene expression profile data were analyzed using Bioconductor packages and Bayesian modeling. Functional annotation of the enriched gene sets was used to define the altered genetic networks. We also analyzed whether gene expression profiles depend on the clinical course and outcome. In order to verify the microarray results, the expression levels of ten functionally relevant genes with high statistical significance (CD177, IL1R2, IL18R1, IL18RAP, OLFM4, TLR5, CPA3, FCER1A, IL5RA, and IL7R) were confirmed by quantitative real-time (qRT) PCR. RESULTS There were 8569 genes displaying differential expression at a significance level of p < 0.05. Following False Discovery Rate (FDR) correction, a total of 5500 genes remained significant at a p-value of < 0.01. Quantitative RT-PCR confirmed the differential expression in 10 selected genes. Functional annotation and network analysis indicated that most of the genes were related to activation of humoral and cellular immune responses (enrichment score 43). Those changes were found in both adults and in children with BM compared to the healthy controls. The gene expression profiles did not significantly depend on the clinical outcome, but there was a strong influence of the specific type of pathogen underlying BM. CONCLUSION This study demonstrates that there is a very strong activation of immune response at the transcriptional level during BM and that the type of pathogen influences this transcriptional activation.
Collapse
|
41
|
Malki K, Campbell J, Davies M, Keers R, Uher R, Ward M, Paya-Cano J, Aitchinson KJ, Binder E, Sluyter F, Kuhn K, Selzer S, Craig I, McGuffin P, Schalkwyk LC. Pharmacoproteomic investigation into antidepressant response in two mouse inbred strains. Proteomics 2013; 12:2355-65. [PMID: 22696452 DOI: 10.1002/pmic.201100306] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In this study, we present a pharmacoproteomic investigation of response to antidepressants two inbred strains. Our aim was to uncover molecular mechanisms underlying antidepressant action and identify new biomarkers to determine therapeutic response to two antidepressants with proven efficacy in the treatment of depression but divergent mechanisms of action. Mice were treated with the pro-noradrenergic drug nortriptyline, the pro-serotonergic drug escitalopram or saline. Quantitative proteomic analyses were undertaken on hippocampal tissue from a study design that used two inbred mouse strains, two depressogenic protocols and a control condition, (maternal separation, chronic mild stress, control), two antidepressant drugs and two dosing protocols. The proteomic analysis was aimed at the identification of specific drug-response markers. Complementary approaches, 2DE and isobaric tandem mass tagging (TMT), were applied to the selected experimental groups. To investigate the relationship between proteomic profiles, depressogenic protocols and drug response, 2DE and TMT data sets were analysed using multivariate methods. The results highlighted significant strain- and stress-related differences across both 2DE and TMT data sets and identified the three gene products involved in serotonergic (PXBD5, YHWAB, SLC25A4) and one in noradrenergic antidepressant action (PXBD6).
Collapse
|
42
|
Jeffries AR, Perfect LW, Ledderose J, Schalkwyk LC, Bray NJ, Mill J, Price J. Stochastic choice of allelic expression in human neural stem cells. Stem Cells 2013; 30:1938-47. [PMID: 22714879 DOI: 10.1002/stem.1155] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Monoallelic gene expression, such as genomic imprinting, is well described. Less well-characterized are genes undergoing stochastic monoallelic expression (MA), where specific clones of cells express just one allele at a given locus. We performed genome-wide allelic expression assessment of human clonal neural stem cells derived from cerebral cortex, striatum, and spinal cord, each with differing genotypes. We assayed three separate clonal lines from each donor, distinguishing stochastic MA from genotypic effects. Roughly 2% of genes showed evidence for autosomal MA, and in about half of these, allelic expression was stochastic between different clones. Many of these loci were known neurodevelopmental genes, such as OTX2 and OLIG2. Monoallelic genes also showed increased levels of DNA methylation compared to hypomethylated biallelic loci. Identified monoallelic gene loci showed altered chromatin signatures in fetal brain, suggesting an in vivo correlate of this phenomenon. We conclude that stochastic allelic expression is prevalent in neural stem cells, providing clonal diversity to developing tissues such as the human brain.
Collapse
|
43
|
Navarrete K, Pedroso I, De Jong S, Stefansson H, Steinberg S, Stefansson K, Ophoff RA, Schalkwyk LC, Collier DA. TCF4 (e2-2; ITF2): a schizophrenia-associated gene with pleiotropic effects on human disease. Am J Med Genet B Neuropsychiatr Genet 2013; 162B:1-16. [PMID: 23129290 DOI: 10.1002/ajmg.b.32109] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/06/2012] [Accepted: 09/27/2012] [Indexed: 12/22/2022]
Abstract
Common SNPs in the transcription factor 4 (TCF4; ITF2, E2-2, SEF-2) gene, which encodes a basic Helix-Loop-Helix (bHLH) transcription factor, are associated with schizophrenia, conferring a small increase in risk. Other common SNPs in the gene are associated with the common eye disorder Fuch's corneal dystrophy, while rare, mostly de novo inactivating mutations cause Pitt-Hopkins syndrome. In this review, we present a systematic bioinformatics and literature review of the genomics, biological function and interactome of TCF4 in the context of schizophrenia. The TCF4 gene is present in all vertebrates, and although protein length varies, there is high conservation of primary sequence, including the DNA binding domain. Humans have a unique leucine-rich nuclear export signal. There are two main isoforms (A and B), as well as complex splicing generating many possible N-terminal amino acid sequences. TCF4 is highly expressed in the brain, where plays a role in neurodevelopment, interacting with class II bHLH transcription factors Math1, HASH1, and neuroD2. The Ca(2+) sensor protein calmodulin interacts with the DNA binding domain of TCF4, inhibiting transcriptional activation. It is also the target of microRNAs, including mir137, which is implicated in schizophrenia. The schizophrenia-associated SNPs are in linkage disequilibrium with common variants within putative DNA regulatory elements, suggesting that regulation of expression may underlie association with schizophrenia. Combined gene co-expression analyses and curated protein-protein interaction data provide a network involving TCF4 and other putative schizophrenia susceptibility genes. These findings suggest new opportunities for understanding the molecular basis of schizophrenia and other mental disorders.
Collapse
|
44
|
Pidsley R, Fernandes C, Viana J, Paya-Cano JL, Liu L, Smith RG, Schalkwyk LC, Mill J. DNA methylation at the Igf2/H19 imprinting control region is associated with cerebellum mass in outbred mice. Mol Brain 2012; 5:42. [PMID: 23216893 PMCID: PMC3541153 DOI: 10.1186/1756-6606-5-42] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2012] [Accepted: 12/02/2012] [Indexed: 11/10/2022] Open
Abstract
Background Insulin-like growth factor 2 (Igf2) is a paternally expressed imprinted gene regulating fetal growth, playing an integral role in the development of many tissues including the brain. The parent-of-origin specific expression of Igf2 is largely controlled by allele-specific DNA methylation at CTCF-binding sites in the imprinting control region (ICR), located immediately upstream of the neighboring H19 gene. Previously we reported evidence of a negative correlation between DNA methylation in this region and cerebellum weight in humans. Results We quantified cerebellar DNA methylation across all four CTCF binding sites spanning the murine Igf2/H19 ICR in an outbred population of Heterogeneous Stock (HS) mice (n = 48). We observe that DNA methylation at the second and third CTCF binding sites in the Igf2/H19 ICR shows a negative relationship with cerebellar mass, reflecting the association observed in human post-mortem cerebellum tissue. Conclusions Given the important role of the cerebellum in motor control and cognition, and the link between structural cerebellar abnormalities and neuropsychiatric phenotypes, the identification of epigenetic factors associated with cerebellum growth and development may provide important insights about the etiology of psychiatric disorders.
Collapse
|
45
|
Davies MN, Volta M, Pidsley R, Lunnon K, Dixit A, Lovestone S, Coarfa C, Harris RA, Milosavljevic A, Troakes C, Al-Sarraj S, Dobson R, Schalkwyk LC, Mill J. Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood. Genome Biol 2012; 13:R43. [PMID: 22703893 PMCID: PMC3446315 DOI: 10.1186/gb-2012-13-6-r43] [Citation(s) in RCA: 499] [Impact Index Per Article: 41.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2012] [Revised: 05/21/2012] [Accepted: 06/15/2012] [Indexed: 01/13/2023] Open
Abstract
Background Dynamic changes to the epigenome play a critical role in establishing and maintaining cellular phenotype during differentiation, but little is known about the normal methylomic differences that occur between functionally distinct areas of the brain. We characterized intra- and inter-individual methylomic variation across whole blood and multiple regions of the brain from multiple donors. Results Distinct tissue-specific patterns of DNA methylation were identified, with a highly significant over-representation of tissue-specific differentially methylated regions (TS-DMRs) observed at intragenic CpG islands and low CG density promoters. A large proportion of TS-DMRs were located near genes that are differentially expressed across brain regions. TS-DMRs were significantly enriched near genes involved in functional pathways related to neurodevelopment and neuronal differentiation, including BDNF, BMP4, CACNA1A, CACA1AF, EOMES, NGFR, NUMBL, PCDH9, SLIT1, SLITRK1 and SHANK3. Although between-tissue variation in DNA methylation was found to greatly exceed between-individual differences within any one tissue, we found that some inter-individual variation was reflected across brain and blood, indicating that peripheral tissues may have some utility in epidemiological studies of complex neurobiological phenotypes. Conclusions This study reinforces the importance of DNA methylation in regulating cellular phenotype across tissues, and highlights genomic patterns of epigenetic variation across functionally distinct regions of the brain, providing a resource for the epigenetics and neuroscience research communities.
Collapse
|
46
|
Keers R, Pedroso I, Breen G, Aitchison KJ, Nolan PM, Cichon S, Nöthen MM, Rietschel M, Schalkwyk LC, Fernandes C. Reduced anxiety and depression-like behaviours in the circadian period mutant mouse afterhours. PLoS One 2012; 7:e38263. [PMID: 22719873 PMCID: PMC3376117 DOI: 10.1371/journal.pone.0038263] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2012] [Accepted: 05/04/2012] [Indexed: 01/27/2023] Open
Abstract
Background Disruption of the circadian rhythm is a key feature of bipolar disorder. Variation in genes encoding components of the molecular circadian clock has been associated with increased risk of the disorder in clinical populations. Similarly in animal models, disruption of the circadian clock can result in altered mood and anxiety which resemble features of human mania; including hyperactivity, reduced anxiety and reduced depression-like behaviour. One such mutant, after hours (Afh), an ENU-derived mutant with a mutation in a recently identified circadian clock gene Fbxl3, results in a disturbed (long) circadian rhythm of approximately 27 hours. Methodology Anxiety, exploratory and depression-like behaviours were evaluated in Afh mice using the open-field, elevated plus maze, light-dark box, holeboard and forced swim test. To further validate findings for human mania, polymorphisms in the human homologue of FBXL3, genotyped by three genome wide case control studies, were tested for association with bipolar disorder. Principal Findings Afh mice showed reduced anxiety- and depression-like behaviour in all of the behavioural tests employed, and some evidence of increased locomotor activity in some tests. An analysis of three separate human data sets revealed a gene wide association between variation in FBXL3 and bipolar disorder (P = 0.009). Conclusions Our results are consistent with previous studies of mutants with extended circadian periods and suggest that disruption of FBXL3 is associated with mania-like behaviours in both mice and humans.
Collapse
|
47
|
|
48
|
Sikk K, Kõks S, Soomets U, Schalkwyk LC, Fernandes C, Haldre S, Aquilonius SM, Taba P. Peripheral blood RNA expression profiling in illicit methcathinone users reveals effect on immune system. Front Genet 2012; 2:42. [PMID: 22303338 PMCID: PMC3268596 DOI: 10.3389/fgene.2011.00042] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2011] [Accepted: 06/22/2011] [Indexed: 02/05/2023] Open
Abstract
Methcathinone (ephedrone) is relatively easily accessible for abuse. Its users develop an extrapyramidal syndrome and it is not known if this is caused by methcathinone itself, by side-ingredients (manganese), or both. In the present study we aimed to clarify molecular mechanisms underlying this condition. We used microarrays to analyze whole-genome gene expression patterns of peripheral blood from 20 methcathinone users and 20 matched controls. Gene expression profile data were analyzed by Bayesian modeling and functional annotation. Of 28,869 genes on the microarrays, 326 showed statistically significant differential expression with FDR adjusted p-values below 0.05. Quantitative real-time PCR confirmed differential expression for the most of the genes selected for validation. Functional annotation and network analysis indicated activation of a gene network that included immunological disease, cellular movement, and cardiovascular disease functions (enrichment score 42). As HIV and HCV infections were confounding factors, we performed additional stratification of subjects. A similar functional activation of the “immunological disease” category was evident when we compared subjects according to injection status (past versus current users, balanced for HIV and HCV infection). However, this difference was not large therefore the major effect was related to the HIV status of the subjects. Mn–methcathinone abusers have blood RNA expression patterns that mostly reflect their HIV and HCV infections.
Collapse
|
49
|
Kõks S, Soomets U, Plaas M, Terasmaa A, Noormets K, Tillmann V, Vasar E, Fernandes C, Schalkwyk LC. Hypothalamic gene expression profile indicates a reduction in G protein signaling in the Wfs1 mutant mice. Physiol Genomics 2011; 43:1351-8. [PMID: 22028430 DOI: 10.1152/physiolgenomics.00117.2011] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The Wfs1 gene codes for a protein with unknown function, but deficiency in this protein results in a range of neuropsychiatric and neuroendocrine syndromes. In the present study we aimed to find the functional networks influenced by Wfs1 in the hypothalamus. We performed gene expression profiling (Mouse Gene 1.0 ST Arrays) in Wfs1-deficient mice; 305 genes were differentially expressed with nominal P value<0.01. FDR (false discovery rate)-adjusted P values were significant (0.007) only for two genes: C4b (t=9.66) and Wfs1 (t=-9.03). However, several genes related to G protein signaling were very close to the FDR-adjusted significance level, such as Rgs4 (regulator of G protein signaling 4) that was downregulated (-0.34, t=-5.4) in Wfs1-deficient mice. Changes in Rgs4 and C4b expression were confirmed by QRT-PCR. In humans, Rgs4 is in the locus for bipolar disease (BPD), and its expression is downregulated in BPD. C4b is a gene related to the neurodegenerative diseases. Functional analysis including the entire data set revealed significant alterations in the canonical pathway "G protein-coupled receptor signaling." The gene expression profile in the hypothalami of the Wfs1 mutant mice was significantly similar to the profiles of following biological functions: psychological disorders, bipolar disorder, mood disorder. In conclusion, hypothalamic gene expression profile resembles with some molecular pathways functionally related to the clinical syndromes in the Wolfram syndrome patients.
Collapse
|
50
|
Dempster EL, Pidsley R, Schalkwyk LC, Owens S, Georgiades A, Kane F, Kalidindi S, Picchioni M, Kravariti E, Toulopoulou T, Murray RM, Mill J. Disease-associated epigenetic changes in monozygotic twins discordant for schizophrenia and bipolar disorder. Hum Mol Genet 2011; 20:4786-96. [PMID: 21908516 PMCID: PMC3221539 DOI: 10.1093/hmg/ddr416] [Citation(s) in RCA: 312] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Studies of the major psychoses, schizophrenia (SZ) and bipolar disorder (BD), have traditionally focused on genetic and environmental risk factors, although more recent work has highlighted an additional role for epigenetic processes in mediating susceptibility. Since monozygotic (MZ) twins share a common DNA sequence, their study represents an ideal design for investigating the contribution of epigenetic factors to disease etiology. We performed a genome-wide analysis of DNA methylation on peripheral blood DNA samples obtained from a unique sample of MZ twin pairs discordant for major psychosis. Numerous loci demonstrated disease-associated DNA methylation differences between twins discordant for SZ and BD individually, and together as a combined major psychosis group. Pathway analysis of our top loci highlighted a significant enrichment of epigenetic changes in biological networks and pathways directly relevant to psychiatric disorder and neurodevelopment. The top psychosis-associated, differentially methylated region, significantly hypomethylated in affected twins, was located in the promoter of ST6GALNAC1 overlapping a previously reported rare genomic duplication observed in SZ. The mean DNA methylation difference at this locus was 6%, but there was considerable heterogeneity between families, with some twin pairs showing a 20% difference in methylation. We subsequently assessed this region in an independent sample of postmortem brain tissue from affected individuals and controls, finding marked hypomethylation (>25%) in a subset of psychosis patients. Overall, our data provide further evidence to support a role for DNA methylation differences in mediating phenotypic differences between MZ twins and in the etiology of both SZ and BD.
Collapse
|