Adaptation of the targeted capture Methyl-Seq platform for the mouse genome identifies novel tissue-specific DNA methylation patterns of genes involved in neurodevelopment

DNA Methylation and Demethylation Underlie the Sex Difference in Estrogen Receptor Alpha in the Arcuate Nucleus

INTRODUCTION

Neurons expressing estrogen receptor (ER) ɑ in the arcuate (ARC) and ventromedial (VMH) nuclei of the hypothalamus sex-specifically control energy homeostasis, sexual behavior, and bone density. Females have more ERɑ neurons in the VMH and ARC than males, and the sex difference in the VMH is eliminated by neonatal treatment with testosterone or a DNA methylation inhibitor.

OBJECTIVE

Here, we tested the roles of testosterone and DNA methylation/demethylation in development of ERɑ in the ARC.

METHODS

ERɑ was examined at birth and weaning in mice that received vehicle or testosterone subcutaneously, and vehicle or DNA methyltransferase inhibitor intracerebroventricularly, as neonates. To examine effects of DNA demethylation on the ERɑ cell number in the ARC, mice were treated neonatally with small interfering RNAs against ten-eleven translocase enzymes. The methylation status of the ERɑ gene (Esr1) was determined in the ARC and VMH using pyrosequencing of bisulfite-converted DNA.

RESULTS

A sex difference in ERɑ in the ARC, favoring females, developed between birth and weaning and was due to programming effects of testosterone. Neonatal inhibition of DNA methylation decreased ERɑ in the ARC of females, and an inhibition of <underline>de</underline>methylation increased ERɑ in the ARC of males. The promoter region of Esr1 exhibited a small sex difference in percent of total methylation in the ARC (females > males) that was opposite to that in the VMH (males > females).

CONCLUSION

DNA methylation and demethylation regulate ERɑ cell number in the ARC, and methylation correlates with activation of Esr1 in this region.

DNA methylome perturbations: an epigenetic basis for the emergingly heritable neurodevelopmental abnormalities associated with maternal smoking and maternal nicotine exposure†

Maternal smoking during pregnancy is associated with an ensemble of neurodevelopmental consequences in children and therefore constitutes a pressing public health concern. Adding to this burden, contemporary epidemiological and especially animal model research suggests that grandmaternal smoking is similarly associated with neurodevelopmental abnormalities in grandchildren, indicative of intergenerational transmission of the neurodevelopmental impacts of maternal smoking. Probing the mechanistic bases of neurodevelopmental anomalies in the children of maternal smokers and the intergenerational transmission thereof, emerging research intimates that epigenetic changes, namely DNA methylome perturbations, are key factors. Altogether, these findings warrant future research to fully elucidate the etiology of neurodevelopmental impairments in the children and grandchildren of maternal smokers and underscore the clear potential thereof to benefit public health by informing the development and implementation of preventative measures, prophylactics, and treatments. To this end, the present review aims to encapsulate the burgeoning evidence linking maternal smoking to intergenerational epigenetic inheritance of neurodevelopmental abnormalities, to identify the strengths and weaknesses thereof, and to highlight areas of emphasis for future human and animal model research therein.

Developmental nicotine exposure engenders intergenerational downregulation and aberrant posttranslational modification of cardinal epigenetic factors in the frontal cortices, striata, and hippocampi of adolescent mice

BACKGROUND

Maternal smoking of traditional or electronic cigarettes during pregnancy, which constitutes developmental nicotine exposure (DNE), heightens the risk of neurodevelopmental disorders including ADHD, autism, and schizophrenia in children. Modeling the intergenerationally transmissible impacts of smoking during pregnancy, we previously demonstrated that both the first- and second-generation adolescent offspring of nicotine-exposed female mice exhibit enhanced nicotine preference, hyperactivity and risk-taking behaviors, aberrant rhythmicity of home cage activity, nicotinic acetylcholine receptor and dopamine transporter dysfunction, impaired furin-mediated proBDNF proteolysis, hypocorticosteronemia-related glucocorticoid receptor hypoactivity, and global DNA hypomethylation in the frontal cortices and striata. This ensemble of multigenerational DNE-induced behavioral, neuropharmacological, neurotrophic, neuroendocrine, and DNA methylomic anomalies recapitulates the pathosymptomatology of neurodevelopmental disorders such as ADHD, autism, and schizophrenia. Further probing the epigenetic bases of DNE-induced multigenerational phenotypic aberrations, the present study examined the expression and phosphorylation of key epigenetic factors via an array of immunoblot experiments.

RESULTS

Data indicate that DNE confers intergenerational deficits in corticostriatal DNA methyltransferase 3A (DNMT3A) expression accompanied by downregulation of methyl-CpG-binding protein 2 (MeCP2) and histone deacetylase 2 (HDAC2) in the frontal cortices and hippocampi, while the expression of ten-eleven translocase methylcytosine dioxygenase 2 (TET2) is unaltered. Moreover, DNE evokes multigenerational abnormalities in HDAC2 (Ser394) but not MeCP2 (Ser421) phosphorylation in the frontal cortices, striata, and hippocampi.

CONCLUSIONS

In light of the extensive gene regulatory roles of DNMT3A, MeCP2, and HDAC2, the findings of this study that DNE elicits downregulation and aberrant posttranslational modification of these factors in both first- and second-generation DNE mice suggest that epigenetic perturbations may constitute a mechanistic hub for the intergenerational transmission of DNE-induced neurodevelopmental disorder-like phenotypes.

Comprehensive analysis of DNA methylation and gene expression in orally tolerized T cells

T cell anergy is known to be a crucial mechanism for various types of immune tolerance, including oral tolerance. The expression of several anergy-specific genes was reportedly up-regulated in anergic T cells, and played important roles in the cells. However, how the genes were up-regulated has not been understood. In this study, we comprehensively analyzed the altered gene expression and DNA methylation status in T cells tolerized by oral antigen in vivo. Our results showed that many genes were significantly up-regulated in the orally tolerized T cells, and most of the genes found in this study have not been reported previously as anergy related genes; for example, ribosomal protein L41 (FC = 3.54E06, p = 3.70E-09: Fisher's exact test; the same applies hereinafter) and CD52 (FC = 2.18E05, p = 3.44E-06). Furthermore, we showed that the DNA methylation statuses of many genes; for example, enoyl-coenzyme A delta isomerase 3 (FC = 3.62E-01, p = 3.01E-02) and leucine zipper protein 1 (FC = 4.80E-01, p = 3.25E-02), including the ones distinctly expressed in tolerized T cells; for example, latexin (FC = 3.85E03, p = 4.06E-02 for expression; FC = 7.75E-01, p = 4.13E-01 for DNA methylation) and small nuclear ribonucleoprotein polypeptide F (FC = 3.12E04, p = 4.46E-04 for expression; FC = 8.56E-01, p = 5.15E-01 for DNA methylation), changed during tolerization, suggesting that the distinct expression of some genes was epigenetically regulated in the tolerized T cells. This study would contribute to providing a novel clue to the fine understanding of the mechanism for T cell anergy and oral tolerance.

A Rat Methyl-Seq Platform to Identify Epigenetic Changes Associated with Stress Exposure

As genomes of a wider variety of animals become available, there is an increasing need for tools that can capture dynamic epigenetic changes in these animal models. The rat is one particular model animal where an epigenetic tool can complement many pharmacological and behavioral studies to provide insightful mechanistic information. To this end, we adapted the SureSelect Target Capture System (referred to as Methyl-Seq) for the rat, which can assess DNA methylation levels across the rat genome. The rat design targeted promoters, CpG islands, island shores, and GC-rich regions from all RefSeq genes. To implement the platform on a rat experiment, male Sprague Dawley rats were exposed to chronic variable stress for 3 weeks, after which blood samples were collected for genomic DNA extraction. Methyl-Seq libraries were constructed from the rat DNA samples by shearing, adapter ligation, target enrichment, bisulfite conversion, and multiplexing. Libraries were sequenced on a next-generation sequencing platform and the sequenced reads were analyzed to identify DMRs between DNA of stressed and unstressed rats. Top candidate DMRs were independently validated by bisulfite pyrosequencing to confirm the robustness of the platform. Results demonstrate that the rat Methyl-Seq platform is a useful epigenetic tool that can capture methylation changes induced by exposure to stress.

Age-specific changes in genome-wide methylation enrich for Foxa2 and estrogen receptor alpha binding sites

The role of DNA methylation patterns in complex phenotypes remains unclear. To explore this question, we adapted our methods for rare variant analysis to characterize genome-wide murine DNA hybridization array to investigate methylation at CpG islands, shores, and regulatory elements. We have applied this platform to compare age and tissue- specific methylation differences in the brain and spleen of young and aged mice. As expected from prior studies, there are clear global differences in organ-specific, but not age-specific, methylation due mostly to changes at repetitive elements. Surprisingly, out of 200,000 loci there were only 946 differentially methylated cytosines (DMCs) between young and old samples (529 hypermethylated, 417 hypomethylated in aged mice) compared to thousands of tissue-specific DMCs. Hypermethylated loci were clustered around the promoter region of Sfi1, exon 2 of Slc11a2, Drg1, Esr1 and Foxa2 transcription factor binding sites. In particular, there were 75 hypermethylated Foxa2 binding sites across a 2.7 Mb region of chromosome 11. Hypomethylated loci were clustered around Mid1, Isoc2b and genome-wide loci with binding sites for Foxa2 and Esr1, which are known to play important roles in development and aging. These data suggest discreet tissue-independent methylation changes associated with aging processes such as cell division (Sfi1, Mid1), energy production (Drg1, Isoc2b) and cell death (Foxa2, Esr1).

Chronic social stress induces DNA methylation changes at an evolutionary conserved intergenic region in chromosome X

Chronic stress resulting from prolonged exposure to negative life events increases the risk of mood and anxiety disorders. Although chronic stress can change gene expression relevant for behavior, molecular regulators of this change have not been fully determined. One process that could play a role is DNA methylation, an epigenetic process whereby a methyl group is added onto nucleotides, predominantly cytosine in the CpG context, and which can be induced by chronic stress. It is unknown to what extent chronic social defeat, a model of human social stress, influences DNA methylation patterns across the genome. Our study addressed this question by using a targeted-capture approach called Methyl-Seq to investigate DNA methylation patterns of the dentate gyrus at putative regulatory regions across the mouse genome from mice exposed to 14 days of social defeat. Findings were replicated in independent cohorts by bisulfite-pyrosequencing. Two differentially methylated regions (DMRs) were identified. One DMR was located at intron 9 of Drosha, and it showed reduced methylation in stressed mice. This observation replicated in one of two independent cohorts. A second DMR was identified at an intergenic region of chromosome X, and methylation in this region was increased in stressed mice. This methylation difference replicated in two independent cohorts and in Major Depressive Disorder (MDD) postmortem brains. These results highlight a region not previously known to be differentially methylated by chronic social defeat stress and which may be involved in MDD.

"Same difference": comprehensive evaluation of four DNA methylation measurement platforms

Background

DNA methylation in CpG context is fundamental to the epigenetic regulation of gene expression in higher eukaryotes. Changes in methylation patterns are implicated in many diseases, cellular differentiation, imprinting, and other biological processes. Techniques that enrich for biologically relevant genomic regions with high CpG content are desired, since, depending on the size of an organism’s methylome, the depth of sequencing required to cover all CpGs can be prohibitively expensive. Currently, restriction enzyme-based reduced representation bisulfite sequencing and its modified protocols are widely used to study methylation differences. Recently, Agilent Technologies, Roche NimbleGen, and Illumina have ventured to both reduce sequencing costs and capture CpGs of known biological relevance by marketing in-solution custom-capture hybridization platforms. We aimed to evaluate the similarities and differences of these four methods considering each platform targets approximately 10–13% of the human methylome.

Results

Overall, the regions covered per platform were as expected: targeted capture-based methods covered > 95% of their designed regions, whereas the restriction enzyme-based method covered > 70% of the expected fragments. While the total number of CpG loci shared by all methods was low, ~ 24% of any platform, the methylation levels of CpGs covered by all platforms were concordant. Annotation of CpG loci with genomic features revealed roughly the same proportions of feature annotations across the four platforms. Targeted capture methods comprise similar types and coverage of annotations and, relative to the targeted methods, the restriction enzyme method covers fewer promoters (~ 9%), CpG shores (~ 8%) and unannotated loci (~ 11%).

Conclusions

Although all methods are largely consistent in terms of covered CpG loci, the commercially available capture methods result in covering nearly all CpG sites in their target regions with few off-target loci and covering similar proportions of annotated CpG loci, the restriction-based enrichment results in more off-target and unannotated CpG loci. Quality of DNA is very important for restriction-based enrichment and starting material can be low. Conversely, quality of the starting material is less important for capture methods, and at least twice the amount of starting material is required. Pricing is marginally less for restriction-based enrichment, and the number of samples that can be prepared is not restricted to the number of capture reactions a kit supports. However, the advantage of capture libraries is the ability to custom design areas of interest. The choice of the technique would be decided by the number of samples, the quality and quantity of DNA available and the biological areas of interest since comparable data are obtained from all platforms.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0190-4) contains supplementary material, which is available to authorized users.

A comprehensive review of genetic and epigenetic mechanisms that regulate BDNF expression and function with relevance to major depressive disorder

Major depressive disorder (MDD) is a mood disorder that affects behavior and impairs cognition. A gene potentially important to this disorder is the brain derived neurotrophic factor (BDNF) as it is involved in processes controlling neuroplasticity. Various mechanisms exist to regulate BDNF's expression level, subcellular localization, and sorting to appropriate secretory pathways. Alterations to these processes by genetic factors and negative stressors can dysregulate its expression, with possible implications for MDD. Here, we review the mechanisms governing the regulation of BDNF expression, and discuss how disease-associated single nucleotide polymorphisms (SNPs) can alter these mechanisms, and influence MDD. As negative stressors increase the likelihood of MDD, we will also discuss the impact of these stressors on BDNF expression, the cellular effect of such a change, and its impact on behavior in animal models of stress. We will also describe epigenetic processes that mediate this change in BDNF expression. Similarities in BDNF expression between animal models of stress and those in MDD will be highlighted. We will also contrast epigenetic patterns at the BDNF locus between animal models of stress, and MDD patients, and address limitations to current clinical studies. Future work should focus on validating current genetic and epigenetic findings in tightly controlled clinical studies. Regions outside of BDNF promoters should also be explored, as should other epigenetic marks, to improve identification of biomarkers for MDD.

Analysis of DNA modifications in aging research

As geroscience research extends into the role of epigenetics in aging and age-related disease, researchers are being confronted with unfamiliar molecular techniques and data analysis methods that can be difficult to integrate into their work. In this review, we focus on the analysis of DNA modifications, namely cytosine methylation and hydroxymethylation, through next-generation sequencing methods. While older techniques for modification analysis performed relative quantitation across regions of the genome or examined average genome levels, these analyses lack the desired specificity, rigor, and genomic coverage to firmly establish the nature of genomic methylation patterns and their response to aging. With recent methodological advances, such as whole genome bisulfite sequencing (WGBS), bisulfite oligonucleotide capture sequencing (BOCS), and bisulfite amplicon sequencing (BSAS), cytosine modifications can now be readily analyzed with base-specific, absolute quantitation at both cytosine-guanine dinucleotide (CG) and non-CG sites throughout the genome or within specific regions of interest by next-generation sequencing. Additional advances, such as oxidative bisulfite conversion to differentiate methylation from hydroxymethylation and analysis of limited input/single-cells, have great promise for continuing to expand epigenomic capabilities. This review provides a background on DNA modifications, the current state-of-the-art for sequencing methods, bioinformatics tools for converting these large data sets into biological insights, and perspectives on future directions for the field.

Sexually divergent DNA methylation patterns with hippocampal aging

DNA methylation is a central regulator of genome function, and altered methylation patterns are indicative of biological aging and mortality. Age‐related cellular, biochemical, and molecular changes in the hippocampus lead to cognitive impairments and greater vulnerability to neurodegenerative disease that varies between the sexes. The role of hippocampal epigenomic changes with aging in these processes is unknown as no genome‐wide analyses of age‐related methylation changes have considered the factor of sex in a controlled animal model. High‐depth, genome‐wide bisulfite sequencing of young (3 month) and old (24 month) male and female mouse hippocampus revealed that while total genomic methylation amounts did not change with aging, specific sites in CG and non‐CG (CH) contexts demonstrated age‐related increases or decreases in methylation that were predominantly sexually divergent. Differential methylation with age for both CG and CH sites was enriched in intergenic and intronic regions and under‐represented in promoters, CG islands, and specific enhancer regions in both sexes, suggesting that certain genomic elements are especially labile with aging, even if the exact genomic loci altered are predominantly sex‐specific. Lifelong sex differences in autosomal methylation at CG and CH sites were also observed. The lack of genome‐wide hypomethylation, sexually divergent aging response, and autosomal sex differences at CG sites was confirmed in human data. These data reveal sex as a previously unappreciated central factor of hippocampal epigenomic changes with aging. In total, these data demonstrate an intricate regulation of DNA methylation with aging by sex, cytosine context, genomic location, and methylation level.

Genome-wide Methyl-Seq analysis of blood-brain targets of glucocorticoid exposure

Chronic exposure to glucocorticoids (GCs) can lead to psychiatric complications through epigenetic mechanisms such as DNA methylation (DNAm). We sought to determine whether epigenetic changes in a peripheral tissue can serve as a surrogate for those in a relatively inaccessible tissue such as the brain. DNA extracted from the hippocampus and blood of mice treated with GCs or vehicle solution was assayed using a genome-wide DNAm platform (Methyl-Seq) to identify differentially methylated regions (DMRs) induced by GC treatment. We observed that ∼70% of the DMRs in both tissues lost methylation following GC treatment. Of the 3,095 DMRs that mapped to the same genes in both tissues, 1,853 DMRs underwent DNAm changes in the same direction. Interestingly, only 209 DMRs (<7%) overlapped in genomic coordinates between the 2 tissues, suggesting tissue-specific differences in GC-targeted loci. Pathway analysis showed that the DMR-associated genes were members of pathways involved in metabolism, immune function, and neurodevelopment. Also, changes in cell type composition of blood and brain were examined by fluorescence-activated cell sorting. Separation of the cortex into neuronal and non-neuronal fractions and the leukocytes into T-cells, B-cells, and neutrophils showed that GC-induced methylation changes primarily occurred in neurons and T-cells, with the blood tissue also undergoing a shift in the proportion of constituent cell types while the proportion of neurons and glia in the brain remained stable. From the current pilot study, we found that despite tissue-specific epigenetic changes and cellular heterogeneity, blood can serve as a surrogate for GC-induced changes in the brain.

Reduced DNA methylation of FKBP5 in Cushing's syndrome

FKBP5 encodes a co-chaperone of HSP90 protein that regulates intracellular glucocorticoid receptor sensitivity. When it is bound to the glucocorticoid receptor complex, cortisol binds with lower affinity to glucocorticoid receptor. Cushing's syndrome is associated with memory deficits, smaller hippocampal volumes, and wide range of cognitive impairments. We aimed at evaluating blood DNA methylation of FKBP5 and its relationship with memory and hippocampal volumes in Cushing's syndrome patients. Polymorphism rs1360780 in FKBP5 has also been assessed to determine whether genetic variations can also govern CpG methylation. Thirty-two Cushing's syndrome patients and 32 matched controls underwent memory tests, 3-Tesla MRI of the brain, and DNA extraction from total leukocytes. DNA samples were bisulfite treated, PCR amplified, and pyrosequenced to assess a total of 41CpG-dinucleotides in the introns 1, 2, 5, and 7 of FKBP5. Significantly lower intronic FKBP5 DNA methylation in CS patients compared to controls was observed in ten CpG-dinucleotides. DNA methylation at these CpGs correlated with left and right HV (Intron-2-Region-2-CpG-3: LHV, r = 0.73, p = 0.02; RHV, r = 0.58, p = 0.03). Cured and active CS patients showed both lower methylation of intron 2 (92.37, 91.8, and 93.34 %, respectively, p = 0.03 for both) and of intron 7 (77.08, 73.74, and 79.71 %, respectively, p = 0.02 and p < 0.01) than controls. Twenty-two subjects had the CC genotype, 34 had the TC genotype, and eight had the TT genotype. Lower average DNA methylation in intron 7 was observed in the TT subjects compared to CC (72.5vs. 79.5 %, p = 0.02) and to TC (72.5 vs. 79.0 %, p = 0.03). Our data demonstrate, for the first time, a reduction of intronic DNA methylation of FKBP5 in CS patients.

Comparison of Methyl-capture Sequencing vs. Infinium 450K methylation array for methylome analysis in clinical samples

Interindividual variability in the epigenome has gained tremendous attention for its potential in pathophysiological investigation, disease diagnosis, and evaluation of clinical intervention. DNA methylation is the most studied epigenetic mark in epigenome-wide association studies (EWAS) as it can be detected from limited starting material. Infinium 450K methylation array is the most popular platform for high-throughput profiling of this mark in clinical samples, as it is cost-effective and requires small amounts of DNA. However, this method suffers from low genome coverage and errors introduced by probe cross-hybridization. Whole-genome bisulfite sequencing can overcome these limitations but elevates the costs tremendously. Methyl-Capture Sequencing (MC Seq) is an attractive intermediate solution to increase the methylome coverage in large sample sets. Here we first demonstrate that MC Seq can be employed using DNA amounts comparable to the amounts used for Infinium 450K. Second, to provide guidance when choosing between the 2 platforms for EWAS, we evaluate and compare MC Seq and Infinium 450K in terms of coverage, technical variation, and concordance of methylation calls in clinical samples. Last, since the focus in EWAS is to study interindividual variation, we demonstrate the utility of MC Seq in studying interindividual variation in subjects from different ethnicities.

Developmental and environmental influences on physiology and behavior--2014 Alan N. Epstein Research Award

Environmental factors acting during development of an individual may influence future health and disease susceptibility. Stressors, including altered diet, psychosocial stress, and immune challenge, during gestation can have negative consequences on the intrauterine environment and increase disease susceptibility of the developing fetus. The long-term effects on offspring have been observed in humans and include greater susceptibility to psychiatric disease, such as depression and anxiety disorders, and adverse metabolic conditions including obesity, diabetes and cardiovascular disease. Studies in my laboratory use rodent models and incorporate a multilevel approach to determine the behavioral, physiological, and neurobiological correlates of disease development as a consequence of early life stressors. The road I took in developing this research program was a rather circuitous one and navigating that path would not have been possible without the many mentors, colleagues, fellows and students who provided critical support. Although my name appears on the plaque of the Alan N. Epstein Research Award, I share this with all those I had the privilege of working with along that road, as briefly summarized in this article.