Quantifying the proportion of different cell types in the human cortex using DNA methylation profiles

BACKGROUND

Due to interindividual variation in the cellular composition of the human cortex, it is essential that covariates that capture these differences are included in epigenome-wide association studies using bulk tissue. As experimentally derived cell counts are often unavailable, computational solutions have been adopted to estimate the proportion of different cell types using DNA methylation data. Here, we validate and profile the use of an expanded reference DNA methylation dataset incorporating two neuronal and three glial cell subtypes for quantifying the cellular composition of the human cortex.

RESULTS

We tested eight reference panels containing different combinations of neuronal- and glial cell types and characterised their performance in deconvoluting cell proportions from computationally reconstructed or empirically derived human cortex DNA methylation data. Our analyses demonstrate that while these novel brain deconvolution models produce accurate estimates of cellular proportions from profiles generated on postnatal human cortex samples, they are not appropriate for the use in prenatal cortex or cerebellum tissue samples. Applying our models to an extensive collection of empirical datasets, we show that glial cells are twice as abundant as neuronal cells in the human cortex and identify significant associations between increased Alzheimer's disease neuropathology and the proportion of specific cell types including a decrease in NeuNNeg/SOX10Neg nuclei and an increase of NeuNNeg/SOX10Pos nuclei.

CONCLUSIONS

Our novel deconvolution models produce accurate estimates for cell proportions in the human cortex. These models are available as a resource to the community enabling the control of cellular heterogeneity in epigenetic studies of brain disorders performed on bulk cortex tissue.

InterpolatedXY: a two-step strategy to normalize DNA methylation microarray data avoiding sex bias

MOTIVATION

Data normalization is an essential step to reduce technical variation within and between arrays. Due to the different karyotypes and the effects of X chromosome inactivation, females and males exhibit distinct methylation patterns on sex chromosomes; thus, it poses a significant challenge to normalize sex chromosome data without introducing bias. Currently, existing methods do not provide unbiased solutions to normalize sex chromosome data, usually, they just process autosomal and sex chromosomes indiscriminately.

RESULTS

Here, we demonstrate that ignoring this sex difference will lead to introducing artificial sex bias, especially for thousands of autosomal CpGs. We present a novel two-step strategy (interpolatedXY) to address this issue, which is applicable to all quantile-based normalization methods. By this new strategy, the autosomal CpGs are first normalized independently by conventional methods, such as funnorm or dasen; then the corrected methylation values of sex chromosome-linked CpGs are estimated as the weighted average of their nearest neighbors on autosomes. The proposed two-step strategy can also be applied to other non-quantile-based normalization methods, as well as other array-based data types. Moreover, we propose a useful concept: the sex explained fraction of variance, to quantitatively measure the normalization effect.

AVAILABILITY AND IMPLEMENTATION

The proposed methods are available by calling the function 'adjustedDasen' or 'adjustedFunnorm' in the latest wateRmelon package (https://github.com/schalkwyk/wateRmelon), with methods compatible with all the major workflows, including minfi.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluation of nanopore sequencing for epigenetic epidemiology: a comparison with DNA methylation microarrays

Most epigenetic epidemiology to date has utilized microarrays to identify positions in the genome where variation in DNA methylation is associated with environmental exposures or disease. However, these profile less than 3% of DNA methylation sites in the human genome, potentially missing affected loci and preventing the discovery of disrupted biological pathways. Third generation sequencing technologies, including Nanopore sequencing, have the potential to revolutionise the generation of epigenetic data, not only by providing genuine genome-wide coverage but profiling epigenetic modifications direct from native DNA. Here we assess the viability of using Nanopore sequencing for epidemiology by performing a comparison with DNA methylation quantified using the most comprehensive microarray available, the Illumina EPIC array. We implemented a CRISPR-Cas9 targeted sequencing approach in concert with Nanopore sequencing to profile DNA methylation in three genomic regions to attempt to rediscover genomic positions that existing technologies have shown are differentially methylated in tobacco smokers. Using Nanopore sequencing reads, DNA methylation was quantified at 1779 CpGs across three regions, providing a finer resolution of DNA methylation patterns compared to the EPIC array. The correlation of estimated levels of DNA methylation between platforms was high. Furthermore, we identified 12 CpGs where hypomethylation was significantly associated with smoking status, including 10 within the AHRR gene. In summary, Nanopore sequencing is a valid option for identifying genomic loci where large differences in DNAm are associated with a phenotype and has the potential to advance our understanding of the role differential methylation plays in the aetiology of complex disease.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

DNA methylation-based sex classifier to predict sex and identify sex chromosome aneuploidy

Background

Sex is an important covariate of epigenome-wide association studies due to its strong influence on DNA methylation patterns across numerous genomic positions. Nevertheless, many samples on the Gene Expression Omnibus (GEO) frequently lack a sex annotation or are incorrectly labelled. Considering the influence that sex imposes on DNA methylation patterns, it is necessary to ensure that methods for filtering poor samples and checking of sex assignment are accurate and widely applicable.

Results

Here we presented a novel method to predict sex using only DNA methylation beta values, which can be readily applied to almost all DNA methylation datasets of different formats (raw IDATs or text files with only signal intensities) uploaded to GEO. We identified 4345 significantly (p<0.01) sex-associated CpG sites present on both 450K and EPIC arrays, and constructed a sex classifier based on the two first principal components of the DNA methylation data of sex-associated probes mapped on sex chromosomes. The proposed method is constructed using whole blood samples and exhibits good performance across a wide range of tissues. We further demonstrated that our method can be used to identify samples with sex chromosome aneuploidy, this function is validated by five Turner syndrome cases and one Klinefelter syndrome case.

Conclusions

This proposed sex classifier not only can be used for sex predictions but also applied to identify samples with sex chromosome aneuploidy, and it is freely and easily accessible by calling the ‘estimateSex’ function from the newest wateRmelon Bioconductor package (https://github.com/schalkwyk/wateRmelon).

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07675-2).

Recalibrating the epigenetic clock: implications for assessing biological age in the human cortex

Human DNA methylation data have been used to develop biomarkers of ageing, referred to as ‘epigenetic clocks’, which have been widely used to identify differences between chronological age and biological age in health and disease including neurodegeneration, dementia and other brain phenotypes. Existing DNA methylation clocks have been shown to be highly accurate in blood but are less precise when used in older samples or in tissue types not included in training the model, including brain. We aimed to develop a novel epigenetic clock that performs optimally in human cortex tissue and has the potential to identify phenotypes associated with biological ageing in the brain. We generated an extensive dataset of human cortex DNA methylation data spanning the life course (n = 1397, ages = 1 to 108 years). This dataset was split into ‘training’ and ‘testing’ samples (training: n = 1047; testing: n = 350). DNA methylation age estimators were derived using a transformed version of chronological age on DNA methylation at specific sites using elastic net regression, a supervised machine learning method. The cortical clock was subsequently validated in a novel independent human cortex dataset (n = 1221, ages = 41 to 104 years) and tested for specificity in a large whole blood dataset (n = 1175, ages = 28 to 98 years). We identified a set of 347 DNA methylation sites that, in combination, optimally predict age in the human cortex. The sum of DNA methylation levels at these sites weighted by their regression coefficients provide the cortical DNA methylation clock age estimate. The novel clock dramatically outperformed previously reported clocks in additional cortical datasets. Our findings suggest that previous associations between predicted DNA methylation age and neurodegenerative phenotypes might represent false positives resulting from clocks not robustly calibrated to the tissue being tested and for phenotypes that become manifest in older ages. The age distribution and tissue type of samples included in training datasets need to be considered when building and applying epigenetic clock algorithms to human epidemiological or disease cohorts.

The DNA methylome of human sperm is distinct from blood with little evidence for tissue-consistent obesity associations

Epidemiological research suggests that paternal obesity may increase the risk of fathering small for gestational age offspring. Studies in non-human mammals indicate that such associations could be mediated by DNA methylation changes in spermatozoa that influence offspring development in utero. Human obesity is associated with differential DNA methylation in peripheral blood. It is unclear, however, whether this differential DNA methylation is reflected in spermatozoa. We profiled genome-wide DNA methylation using the Illumina MethylationEPIC array in a cross-sectional study of matched human blood and sperm from lean (discovery n = 47; replication n = 21) and obese (n = 22) males to analyse tissue covariation of DNA methylation, and identify obesity-associated methylomic signatures. We found that DNA methylation signatures of human blood and spermatozoa are highly discordant, and methylation levels are correlated at only a minority of CpG sites (~1%). At the majority of these sites, DNA methylation appears to be influenced by genetic variation. Obesity-associated DNA methylation in blood was not generally reflected in spermatozoa, and obesity was not associated with altered covariation patterns or accelerated epigenetic ageing in the two tissues. However, one cross-tissue obesity-specific hypermethylated site (cg19357369; chr4:2429884; P = 8.95 × 10-8; 2% DNA methylation difference) was identified, warranting replication and further investigation. When compared to a wide range of human somatic tissue samples (n = 5,917), spermatozoa displayed differential DNA methylation across pathways enriched in transcriptional regulation. Overall, human sperm displays a unique DNA methylation profile that is highly discordant to, and practically uncorrelated with, that of matched peripheral blood. We observed that obesity was only nominally associated with differential DNA methylation in sperm, and therefore suggest that spermatozoal DNA methylation is an unlikely mediator of intergenerational effects of metabolic traits.

Systematic underestimation of the epigenetic clock and age acceleration in older subjects

Background

The Horvath epigenetic clock is widely used. It predicts age quite well from 353 CpG sites in the DNA methylation profile in unknown samples and has been used to calculate “age acceleration” in various tissues and environments.

Results

The model systematically underestimates age in tissues from older people. This is seen in all examined tissues but most strongly in the cerebellum and is consistently observed in multiple datasets. Age acceleration is thus age-dependent, and this can lead to spurious associations. The current literature includes examples of association tests with age acceleration calculated in a wide variety of ways.

Conclusions

The concept of an epigenetic clock is compelling, but caution should be taken in interpreting associations with age acceleration. Association tests of age acceleration should include age as a covariate.

DNA methylation aging clocks: challenges and recommendations

Epigenetic clocks comprise a set of CpG sites whose DNA methylation levels measure subject age. These clocks are acknowledged as a highly accurate molecular correlate of chronological age in humans and other vertebrates. Also, extensive research is aimed at their potential to quantify biological aging rates and test longevity or rejuvenating interventions. Here, we discuss key challenges to understand clock mechanisms and biomarker utility. This requires dissecting the drivers and regulators of age-related changes in single-cell, tissue- and disease-specific models, as well as exploring other epigenomic marks, longitudinal and diverse population studies, and non-human models. We also highlight important ethical issues in forensic age determination and predicting the trajectory of biological aging in an individual.

RNA sequencing of identical twins discordant for autism reveals blood-based signatures implicating immune and transcriptional dysregulation

Background

A gap exists in our mechanistic understanding of how genetic and environmental risk factors converge at the molecular level to result in the emergence of autism symptoms. We compared blood-based gene expression signatures in identical twins concordant and discordant for autism spectrum condition (ASC) to differentiate genetic and environmentally driven transcription differences, and establish convergent evidence for biological mechanisms involved in ASC.

Methods

Genome-wide gene expression data were generated using RNA-seq on whole blood samples taken from 16 pairs of monozygotic (MZ) twins and seven twin pair members (39 individuals in total), who had been assessed for ASC and autism traits at age 12. Differential expression (DE) analyses were performed between (a) affected and unaffected subjects (N = 36) and (b) within discordant ASC MZ twin pairs (total N = 11) to identify environmental-driven DE. Gene set enrichment and pathway testing was performed on DE gene lists. Finally, an integrative analysis using DNA methylation data aimed to identify genes with consistent evidence for altered regulation in cis.

Results

In the discordant twin analysis, three genes showed evidence for DE at FDR < 10%: IGHG4, EVI2A and SNORD15B. In the case-control analysis, four DE genes were identified at FDR < 10% including IGHG4, PRR13P5, DEPDC1B, and ZNF501. We find enrichment for DE of genes curated in the SFARI human gene database. Pathways showing evidence of enrichment included those related to immune cell signalling and immune response, transcriptional control and cell cycle/proliferation. Integrative methylomic and transcriptomic analysis identified a number of genes showing suggestive evidence for cis dysregulation.

Limitations

Identical twins stably discordant for ASC are rare, and as such the sample size was limited and constrained to the use of peripheral blood tissue for transcriptomic and methylomic profiling. Given these primary limitations, we focused on transcript-level analysis.

Conclusions

Using a cohort of ASC discordant and concordant MZ twins, we add to the growing body of transcriptomic-based evidence for an immune-based component in the molecular aetiology of ASC. Whilst the sample size was limited, the study demonstrates the utility of the discordant MZ twin design combined with multi-omics integration for maximising the potential to identify disease-associated molecular signals.

Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic autism in post-mortem human brain tissue

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

Bigmelon: tools for analysing large DNA methylation datasets

MOTIVATION

The datasets generated by DNA methylation analyses are getting bigger. With the release of the HumanMethylationEPIC micro-array and datasets containing thousands of samples, analyses of these large datasets using R are becoming impractical due to large memory requirements. As a result there is an increasing need for computationally efficient methodologies to perform meaningful analysis on high dimensional data.

RESULTS

Here we introduce the bigmelon R package, which provides a memory efficient workflow that enables users to perform the complex, large scale analyses required in epigenome wide association studies (EWAS) without the need for large RAM. Building on top of the CoreArray Genomic Data Structure file format and libraries packaged in the gdsfmt package, we provide a practical workflow that facilitates the reading-in, preprocessing, quality control and statistical analysis of DNA methylation data.We demonstrate the capabilities of the bigmelon package using a large dataset consisting of 1193 human blood samples from the Understanding Society: UK Household Longitudinal Study, assayed on the EPIC micro-array platform.

AVAILABILITY AND IMPLEMENTATION

The bigmelon package is available on Bioconductor (http://bioconductor.org/packages/bigmelon/). The Understanding Society dataset is available at https://www.understandingsociety.ac.uk/about/health/data upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Expression quantitative trait loci in the developing human brain and their enrichment in neuropsychiatric disorders

BACKGROUND

Genetic influences on gene expression in the human fetal brain plausibly impact upon a variety of postnatal brain-related traits, including susceptibility to neuropsychiatric disorders. However, to date, there have been no studies that have mapped genome-wide expression quantitative trait loci (eQTL) specifically in the human prenatal brain.

RESULTS

We performed deep RNA sequencing and genome-wide genotyping on a unique collection of 120 human brains from the second trimester of gestation to provide the first eQTL dataset derived exclusively from the human fetal brain. We identify high confidence cis-acting eQTL at the individual transcript as well as whole gene level, including many mapping to a common inversion polymorphism on chromosome 17q21. Fetal brain eQTL are enriched among risk variants for postnatal conditions including attention deficit hyperactivity disorder, schizophrenia, and bipolar disorder. We further identify changes in gene expression within the prenatal brain that potentially mediate risk for neuropsychiatric traits, including increased expression of C4A in association with genetic risk for schizophrenia, increased expression of LRRC57 in association with genetic risk for bipolar disorder, and altered expression of multiple genes within the chromosome 17q21 inversion in association with variants influencing the personality trait of neuroticism.

CONCLUSIONS

We have mapped eQTL operating in the human fetal brain, providing evidence that these confer risk to certain neuropsychiatric disorders, and identifying gene expression changes that potentially mediate susceptibility to these conditions.

Leveraging DNA-Methylation Quantitative-Trait Loci to Characterize the Relationship between Methylomic Variation, Gene Expression, and Complex Traits

Characterizing the complex relationship between genetic, epigenetic, and transcriptomic variation has the potential to increase understanding about the mechanisms underpinning health and disease phenotypes. We undertook a comprehensive analysis of common genetic variation on DNA methylation (DNAm) by using the Illumina EPIC array to profile samples from the UK Household Longitudinal study. We identified 12,689,548 significant DNA methylation quantitative trait loci (mQTL) associations (p < 6.52 × 10-14) occurring between 2,907,234 genetic variants and 93,268 DNAm sites, including a large number not identified by previous DNAm-profiling methods. We demonstrate the utility of these data for interpreting the functional consequences of common genetic variation associated with > 60 human traits by using summary-data-based Mendelian randomization (SMR) to identify 1,662 pleiotropic associations between 36 complex traits and 1,246 DNAm sites. We also use SMR to characterize the relationship between DNAm and gene expression and thereby identify 6,798 pleiotropic associations between 5,420 DNAm sites and the transcription of 1,702 genes. Our mQTL database and SMR results are available via a searchable online database as a resource to the research community.

A histone acetylome-wide association study of Alzheimer's disease identifies disease-associated H3K27ac differences in the entorhinal cortex

We quantified genome-wide patterns of lysine H3K27 acetylation (H3K27ac) in entorhinal cortex samples from Alzheimer's disease (AD) cases and matched controls using chromatin immunoprecipitation and highly parallel sequencing. We observed widespread acetylomic variation associated with AD neuropathology, identifying 4,162 differential peaks (false discovery rate < 0.05) between AD cases and controls. Differentially acetylated peaks were enriched in disease-related biological pathways and included regions annotated to genes involved in the progression of amyloid-β and tau pathology (for example, APP, PSEN1, PSEN2, and MAPT), as well as regions containing variants associated with sporadic late-onset AD. Partitioned heritability analysis highlighted a highly significant enrichment of AD risk variants in entorhinal cortex H3K27ac peak regions. AD-associated variable H3K27ac was associated with transcriptional variation at proximal genes including CR1, GPR22, KMO, PIM3, PSEN1, and RGCC. In addition to identifying molecular pathways associated with AD neuropathology, we present a framework for genome-wide studies of histone modifications in complex disease.

Schizophrenia-associated methylomic variation: molecular signatures of disease and polygenic risk burden across multiple brain regions

Genetic association studies provide evidence for a substantial polygenic component to schizophrenia, although the neurobiological mechanisms underlying the disorder remain largely undefined. Building on recent studies supporting a role for developmentally regulated epigenetic variation in the molecular aetiology of schizophrenia, this study aimed to identify epigenetic variation associated with both a diagnosis of schizophrenia and elevated polygenic risk burden for the disease across multiple brain regions. Genome-wide DNA methylation was quantified in 262 post-mortem brain samples, representing tissue from four brain regions (prefrontal cortex, striatum, hippocampus and cerebellum) from 41 schizophrenia patients and 47 controls. We identified multiple disease-associated and polygenic risk score-associated differentially methylated positions and regions, which are not enriched in genomic regions identified in genetic studies of schizophrenia and do not reflect direct genetic effects on DNA methylation. Our study represents the first analysis of epigenetic variation associated with schizophrenia across multiple brain regions and highlights the utility of polygenic risk scores for identifying molecular pathways associated with aetiological variation in complex disease.

5-hydroxymethylcytosine is highly dynamic across human fetal brain development

Background

Epigenetic processes play a key role in orchestrating transcriptional regulation during the development of the human central nervous system. We previously described dynamic changes in DNA methylation (5mC) occurring during human fetal brain development, but other epigenetic processes operating during this period have not been extensively explored. Of particular interest is DNA hydroxymethylation (5hmC), a modification that is enriched in the human brain and hypothesized to play an important role in neuronal function, learning and memory. In this study, we quantify 5hmC across the genome of 71 human fetal brain samples spanning 23 to 184 days post-conception.

Results

We identify widespread changes in 5hmC occurring during human brain development, notable sex-differences in 5hmC in the fetal brain, and interactions between 5mC and 5hmC at specific sites. Finally, we identify loci where 5hmC in the fetal brain is associated with genetic variation.

Conclusions

This study represents the first systematic analysis of dynamic changes in 5hmC across human neurodevelopment and highlights the potential importance of this modification in the human brain. A searchable database of our fetal brain 5hmC data is available as a resource to the research community at http://www.epigenomicslab.com/online-data-resources.

Electronic supplementary material

The online version of this article doi:(10.1186/s12864-017-4091-x) contains supplementary material, which is available to authorized users.

Paternal Age Alters Social Development in Offspring

OBJECTIVE

Advanced paternal age (APA) at conception has been linked with autism and schizophrenia in offspring, neurodevelopmental disorders that affect social functioning. The current study explored the effects of paternal age on social development in the general population.

METHOD

We used multilevel growth modeling to investigate APA effects on socioemotional development from early childhood until adolescence, as measured by the Strengths and Difficulties Questionnaire (SDQ) in the Twins Early Development Study (TEDS) sample. We also investigated genetic and environmental underpinnings of the paternal age effects on development, using the Additive genetics, Common environment, unique Environment (ACE) and gene-environment (GxE) models.

RESULTS

In the general population, both very young and advanced paternal ages were associated with altered trajectory of social development (intercept: p = .01; slope: p = .03). No other behavioral domain was affected by either young or advanced age at fatherhood, suggesting specificity of paternal age effects. Increased importance of genetic factors in social development was recorded in the offspring of older but not very young fathers, suggesting distinct underpinnings of the paternal age effects at these two extremes.

CONCLUSION

Our findings highlight that the APA-related deficits that lead to autism and schizophrenia are likely continuously distributed in the population.

Highly polygenic architecture of antidepressant treatment response: Comparative analysis of SSRI and NRI treatment in an animal model of depression

Response to antidepressant (AD) treatment may be a more polygenic trait than previously hypothesized, with many genetic variants interacting in yet unclear ways. In this study we used methods that can automatically learn to detect patterns of statistical regularity from a sparsely distributed signal across hippocampal transcriptome measurements in a large-scale animal pharmacogenomic study to uncover genomic variations associated with AD. The study used four inbred mouse strains of both sexes, two drug treatments, and a control group (escitalopram, nortriptyline, and saline). Multi-class and binary classification using Machine Learning (ML) and regularization algorithms using iterative and univariate feature selection methods, including InfoGain, mRMR, ANOVA, and Chi Square, were used to uncover genomic markers associated with AD response. Relevant genes were selected based on Jaccard distance and carried forward for gene-network analysis. Linear association methods uncovered only one gene associated with drug treatment response. The implementation of ML algorithms, together with feature reduction methods, revealed a set of 204 genes associated with SSRI and 241 genes associated with NRI response. Although only 10% of genes overlapped across the two drugs, network analysis shows that both drugs modulated the CREB pathway, through different molecular mechanisms. Through careful implementation and optimisations, the algorithms detected a weak signal used to predict whether an animal was treated with nortriptyline (77%) or escitalopram (67%) on an independent testing set. The results from this study indicate that the molecular signature of AD treatment may include a much broader range of genomic markers than previously hypothesized, suggesting that response to medication may be as complex as the pathology. The search for biomarkers of antidepressant treatment response could therefore consider a higher number of genetic markers and their interactions. Through predominately different molecular targets and mechanisms of action, the two drugs modulate the same Creb1 pathway which plays a key role in neurotrophic responses and in inflammatory processes. © 2016 The Authors. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics Published by Wiley Periodicals, Inc.

Transcriptome analysis of genes and gene networks involved in aggressive behavior in mouse and zebrafish

Despite moderate heritability estimates, the molecular architecture of aggressive behavior remains poorly characterized. This study compared gene expression profiles from a genetic mouse model of aggression with zebrafish, an animal model traditionally used to study aggression. A meta-analytic, cross-species approach was used to identify genomic variants associated with aggressive behavior. The Rankprod algorithm was used to evaluated mRNA differences from prefrontal cortex tissues of three sets of mouse lines (N = 18) selectively bred for low and high aggressive behavior (SAL/LAL, TA/TNA, and NC900/NC100). The same approach was used to evaluate mRNA differences in zebrafish (N = 12) exposed to aggressive or non-aggressive social encounters. Results were compared to uncover genes consistently implicated in aggression across both studies. Seventy-six genes were differentially expressed (PFP < 0.05) in aggressive compared to non-aggressive mice. Seventy genes were differentially expressed in zebrafish exposed to a fight encounter compared to isolated zebrafish. Seven genes (Fos, Dusp1, Hdac4, Ier2, Bdnf, Btg2, and Nr4a1) were differentially expressed across both species 5 of which belonging to a gene-network centred on the c-Fos gene hub. Network analysis revealed an association with the MAPK signaling cascade. In human studies HDAC4 haploinsufficiency is a key genetic mechanism associated with brachydactyly mental retardation syndrome (BDMR), which is associated with aggressive behaviors. Moreover, the HDAC4 receptor is a drug target for valproic acid, which is being employed as an effective pharmacological treatment for aggressive behavior in geriatric, psychiatric, and brain-injury patients. © 2016 Wiley Periodicals, Inc.

Effects of advanced paternal age on trajectories of social behavior in offspring

Our study is the first investigation of the effects of advanced paternal age (APA) on the developmental trajectory of social behavior in rodent offspring. Given the strong epidemiological association between APA and sexually dimorphic neurodevelopmental disorders that are characterized by abnormalities in social behavior (autism, schizophrenia), we assessed sociability in male and female inbred mice (C57BL/6J) across postnatal development (N = 104) in relation to paternal age. We found differences in early social behavior in both male and female offspring of older breeders, with differences in this social domain persisting into adulthood in males only. We showed that these social deficits were not present in the fathers of these offspring, confirming a de novo origin of an altered social trajectory in the offspring generation. Our results, highly novel in rodent research, support the epidemiological observations in humans and provide evidence for a causal link between APA, age-related changes in the paternal sperm DNA and neurodevelopmental disorders in their offspring.

Comparative mRNA analysis of behavioral and genetic mouse models of aggression

Mouse models of aggression have traditionally compared strains, most notably BALB/cJ and C57BL/6. However, these strains were not designed to study aggression despite differences in aggression-related traits and distinct reactivity to stress. This study evaluated expression of genes differentially regulated in a stress (behavioral) mouse model of aggression with those from a recent genetic mouse model aggression. The study used a discovery-replication design using two independent mRNA studies from mouse brain tissue. The discovery study identified strain (BALB/cJ and C57BL/6J) × stress (chronic mild stress or control) interactions. Probe sets differentially regulated in the discovery set were intersected with those uncovered in the replication study, which evaluated differences between high and low aggressive animals from three strains specifically bred to study aggression. Network analysis was conducted on overlapping genes uncovered across both studies. A significant overlap was found with the genetic mouse study sharing 1,916 probe sets with the stress model. Fifty-one probe sets were found to be strongly dysregulated across both studies mapping to 50 known genes. Network analysis revealed two plausible pathways including one centered on the UBC gene hub which encodes ubiquitin, a protein well-known for protein degradation, and another on P38 MAPK. Findings from this study support the stress model of aggression, which showed remarkable molecular overlap with a genetic model. The study uncovered a set of candidate genes including the Erg2 gene, which has previously been implicated in different psychopathologies. The gene networks uncovered points at a Redox pathway as potentially being implicated in aggressive related behaviors.

Tissue-specific patterns of allelically-skewed DNA methylation

While DNA methylation is usually thought to be symmetrical across both alleles, there are some notable exceptions. Genomic imprinting and X chromosome inactivation are two well-studied sources of allele-specific methylation (ASM), but recent research has indicated a more complex pattern in which genotypic variation can be associated with allelically-skewed DNA methylation in cis. Given the known heterogeneity of DNA methylation across tissues and cell types we explored inter- and intra-individual variation in ASM across several regions of the human brain and whole blood from multiple individuals. Consistent with previous studies, we find widespread ASM with > 4% of the ∼220,000 loci interrogated showing evidence of allelically-skewed DNA methylation. We identify ASM flanking known imprinted regions, and show that ASM sites are enriched in DNase I hypersensitivity sites and often located in an extended genomic context of intermediate DNA methylation. We also detect examples of genotype-driven ASM, some of which are tissue-specific. These findings contribute to our understanding of the nature of differential DNA methylation across tissues and have important implications for genetic studies of complex disease. As a resource to the community, ASM patterns across each of the tissues studied are available in a searchable online database: http://epigenetics.essex.ac.uk/ASMBrainBlood.

Methylation QTLs in the developing brain and their enrichment in schizophrenia risk loci

We characterized DNA methylation quantitative trait loci (mQTLs) in a large collection (n = 166) of human fetal brain samples spanning 56-166 d post-conception, identifying >16,000 fetal brain mQTLs. Fetal brain mQTLs were primarily cis-acting, enriched in regulatory chromatin domains and transcription factor binding sites, and showed substantial overlap with genetic variants that were also associated with gene expression in the brain. Using tissue from three distinct regions of the adult brain (prefrontal cortex, striatum and cerebellum), we found that most fetal brain mQTLs were developmentally stable, although a subset was characterized by fetal-specific effects. Fetal brain mQTLs were enriched amongst risk loci identified in a recent large-scale genome-wide association study (GWAS) of schizophrenia, a severe psychiatric disorder with a hypothesized neurodevelopmental component. Finally, we found that mQTLs can be used to refine GWAS loci through the identification of discrete sites of variable fetal brain methylation associated with schizophrenia risk variants.

Methylomic profiling of human brain tissue supports a neurodevelopmental origin for schizophrenia

BACKGROUND

Schizophrenia is a severe neuropsychiatric disorder that is hypothesized to result from disturbances ine arly brain development. There is mounting evidence to support a role for developmentally regulated epigenetic variation in the molecular etiology of the disorder. Here, we describe a systematic study of schizophrenia-associated methylomic variation in the adult brain and its relationship to changes in DNA methylation across human fetal brain development.

RESULTS

We profile methylomic variation in matched prefrontal cortex and cerebellum brain tissue from schizophrenia patients and controls, identifying disease-associated differential DNA methylation at multiple loci,particularly in the prefrontal cortex, and confirming these differences in an independent set of adult brain samples.Our data reveal discrete modules of co-methylated loci associated with schizophrenia that are enriched for genes involved in neurodevelopmental processes and include loci implicated by genetic studies of the disorder. Methylomic data from human fetal cortex samples, spanning 23 to 184 days post-conception, indicates that schizophrenia-associated differentially methylated positions are significantly enriched for loci at which DNA methylation is dynamically altered during human fetal brain development.

CONCLUSIONS

Our data support the hypothesis that schizophrenia has an important early neurodevelopmental component, and suggest that epigenetic mechanisms may mediate these effects.

Pervasive and opposing effects of Unpredictable Chronic Mild Stress (UCMS) on hippocampal gene expression in BALB/cJ and C57BL/6J mouse strains

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.

Intermediate DNA methylation is a conserved signature of genome regulation

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.

Methylomic trajectories across human fetal brain development

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.

Methylomic profiling implicates cortical deregulation of ANK1 in Alzheimer's disease

Alzheimer's disease (AD) is a chronic neurodegenerative disorder that is characterized by progressive neuropathology and cognitive decline. We performed a cross-tissue analysis of methylomic variation in AD using samples from four independent human post-mortem brain cohorts. We identified a differentially methylated region in the ankyrin 1 (ANK1) gene that was associated with neuropathology in the entorhinal cortex, a primary site of AD manifestation. This region was confirmed as being substantially hypermethylated in two other cortical regions (superior temporal gyrus and prefrontal cortex), but not in the cerebellum, a region largely protected from neurodegeneration in AD, or whole blood obtained pre-mortem from the same individuals. Neuropathology-associated ANK1 hypermethylation was subsequently confirmed in cortical samples from three independent brain cohorts. This study represents, to the best of our knowledge, the first epigenome-wide association study of AD employing a sequential replication design across multiple tissues and highlights the power of this approach for identifying methylomic variation associated with complex disease.

Alzheimer's disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci

Here, we leverage a unique 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 find that the level of methylation at 71 of the 415,848 interrogated CpGs is significantly associated with the burden of AD pathology, including CpGs in the ABCA7 and BIN1 regions, which harbor known AD susceptibility variants. We validate 11 of the differentially methylated regions in an independent set of 117 subjects. Further, we functionally validate these CpG associations and identify the nearby genes whose RNA expression is altered in AD: ANK1, CDH23, DIP2A, RHBDF2, RPL13, RNF34, SERPINF1 and SERPINF2. Our analyses suggest that these DNA methylation changes may have a role in the onset of AD since (1) they are seen in presymptomatic subjects and (2) six of the validated genes connect to a known AD susceptibility gene network.

Integrative mouse and human mRNA studies using WGCNA nominates novel candidate genes involved in the pathogenesis of major depressive disorder

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.

Mood-stabilizers differentially affect housekeeping gene expression in human cells

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.

The endogenous and reactive depression subtypes revisited: integrative animal and human studies implicate multiple distinct molecular mechanisms underlying major depressive disorder

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.

Transcriptomic changes in the frontal cortex associated with paternal age

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.

Putative transcriptomic biomarkers in the inflammatory cytokine pathway differentiate major depressive disorder patients from control subjects and bipolar disorder patients

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.

The inflammatory cytokines: molecular biomarkers for major depressive disorder?

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.

ATP-binding cassette sub-family F member 1 (ABCF1) is identified as a putative therapeutic target of escitalopram in the inflammatory cytokine pathway

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.

A data-driven approach to preprocessing Illumina 450K methylation array data

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.

Peripheral blood RNA gene expression profiling in patients with bacterial meningitis

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.

Pharmacoproteomic investigation into antidepressant response in two mouse inbred strains

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

Stochastic choice of allelic expression in human neural stem cells

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.

TCF4 (e2-2; ITF2): a schizophrenia-associated gene with pleiotropic effects on human disease

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.

DNA methylation at the Igf2/H19 imprinting control region is associated with cerebellum mass in outbred mice

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.

Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood

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.

Reduced anxiety and depression-like behaviours in the circadian period mutant mouse afterhours

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.

Peripheral blood RNA expression profiling in illicit methcathinone users reveals effect on immune system

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.

Hypothalamic gene expression profile indicates a reduction in G protein signaling in the Wfs1 mutant mice

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.