Developmentally dynamic changes in DNA methylation in the human pancreas.

Development of the human pancreas requires the precise temporal control of gene expression via epigenetic mechanisms and the binding of key transcription factors. We quantified genome-wide patterns of DNA methylation in human fetal pancreatic samples from donors aged 6 to 21 post-conception weeks. We found dramatic changes in DNA methylation across pancreas development, with > 21% of sites characterized as developmental differentially methylated positions (dDMPs) including many annotated to genes associated with monogenic diabetes. An analysis of DNA methylation in postnatal pancreas tissue showed that the dramatic temporal changes in DNA methylation occurring in the developing pancreas are largely limited to the prenatal period. Significant differences in DNA methylation were observed between males and females at a number of autosomal sites, with a small proportion of sites showing sex-specific DNA methylation trajectories across pancreas development. Pancreas dDMPs were not distributed equally across the genome and were depleted in regulatory domains characterized by open chromatin and the binding of known pancreatic development transcription factors. Finally, we compared our pancreas dDMPs to previous findings from the human brain, identifying evidence for tissue-specific developmental changes in DNA methylation. This study represents the first systematic exploration of DNA methylation patterns during human fetal pancreas development and confirms the prenatal period as a time of major epigenomic plasticity.

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.

Uncertainty quantification of reference-based cellular deconvolution algorithms.

The majority of epigenetic epidemiology studies to date have generated genome-wide profiles from bulk tissues (e.g., whole blood) however these are vulnerable to confounding from variation in cellular composition. Proxies for cellular composition can be mathematically derived from the bulk tissue profiles using a deconvolution algorithm; however, there is no method to assess the validity of these estimates for a dataset where the true cellular proportions are unknown. In this study, we describe, validate and characterize a sample level accuracy metric for derived cellular heterogeneity variables. The CETYGO score captures the deviation between a sample's DNA methylation profile and its expected profile given the estimated cellular proportions and cell type reference profiles. We demonstrate that the CETYGO score consistently distinguishes inaccurate and incomplete deconvolutions when applied to reconstructed whole blood profiles. By applying our novel metric to >6,300 empirical whole blood profiles, we find that estimating accurate cellular composition is influenced by both technical and biological variation. In particular, we show that when using a common reference panel for whole blood, less accurate estimates are generated for females, neonates, older individuals and smokers. Our results highlight the utility of a metric to assess the accuracy of cellular deconvolution, and describe how it can enhance studies of DNA methylation that are reliant on statistical proxies for cellular heterogeneity. To facilitate incorporating our methodology into existing pipelines, we have made it freely available as an R package (https://github.com/ds420/CETYGO).

DNA methylation signatures of Alzheimer's disease neuropathology in the cortex are primarily driven by variation in non-neuronal cell-types.

Alzheimer's disease (AD) is a chronic neurodegenerative disease characterized by the progressive accumulation of amyloid-beta and neurofibrillary tangles of tau in the neocortex. We profiled DNA methylation in two regions of the cortex from 631 donors, performing an epigenome-wide association study of multiple measures of AD neuropathology. We meta-analyzed our results with those from previous studies of DNA methylation in AD cortex (total n = 2013 donors), identifying 334 cortical differentially methylated positions (DMPs) associated with AD pathology including methylomic variation at loci not previously implicated in dementia. We subsequently profiled DNA methylation in NeuN+ (neuronal-enriched), SOX10+ (oligodendrocyte-enriched) and NeuN-/SOX10- (microglia- and astrocyte-enriched) nuclei, finding that the majority of DMPs identified in 'bulk' cortex tissue reflect DNA methylation differences occurring in non-neuronal cells. Our study highlights the power of utilizing multiple measures of neuropathology to identify epigenetic signatures of AD and the importance of characterizing disease-associated variation in purified cell-types.

Functional characterization of the schizophrenia associated gene AS3MT identifies a role in neuronal development.

Genome-wide association studies (GWAS) have identified multiple genomic regions associated with schizophrenia, although many variants reside in noncoding regions characterized by high linkage disequilibrium (LD) making the elucidation of molecular mechanisms challenging. A genomic region on chromosome 10q24 has been consistently associated with schizophrenia with risk attributed to the AS3MT gene. Although AS3MT is hypothesized to play a role in neuronal development and differentiation, work to fully understand the function of this gene has been limited. In this study we explored the function of AS3MT using a neuronal cell line (SH-SY5Y). We confirm previous findings of isoform specific expression of AS3MT during SH-SY5Y differentiation toward neuronal fates. Using CRISPR-Cas9 gene editing we generated AS3MT knockout SH-SY5Y cell lines and used RNA-seq to identify significant changes in gene expression in pathways associated with neuronal development, inflammation, extracellular matrix formation, and RNA processing, including dysregulation of other genes strongly implicated in schizophrenia. We did not observe any morphological changes in cell size and neurite length following neuronal differentiation and MAP2 immunocytochemistry. These results provide novel insights into the potential role of AS3MT in brain development and identify pathways through which genetic variation in this region may confer risk for schizophrenia.

Genome-wide characterization of mitochondrial DNA methylation in human brain.

Background: There is growing interest in the role of DNA methylation in regulating the transcription of mitochondrial genes, particularly in brain disorders characterized by mitochondrial dysfunction. Here, we present a novel approach to interrogate the mitochondrial DNA methylome at single base resolution using targeted bisulfite sequencing. We applied this method to investigate mitochondrial DNA methylation patterns in post-mortem superior temporal gyrus and cerebellum brain tissue from seven human donors. Results: We show that mitochondrial DNA methylation patterns are relatively low but conserved, with peaks in DNA methylation at several sites, such as within the D-LOOP and the genes MT-ND2, MT-ATP6, MT-ND4, MT-ND5 and MT-ND6, predominantly in a non-CpG context. The elevated DNA methylation we observe in the D-LOOP we validate using pyrosequencing. We identify loci that show differential DNA methylation patterns associated with age, sex and brain region. Finally, we replicate previously reported differentially methylated regions between brain regions from a methylated DNA immunoprecipitation sequencing study. Conclusions: We have annotated patterns of DNA methylation at single base resolution across the mitochondrial genome in human brain samples. Looking to the future this approach could be utilized to investigate the role of mitochondrial epigenetic mechanisms in disorders that display mitochondrial dysfunction.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing.

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

Sites of active gene regulation in the prenatal frontal cortex and their role in neuropsychiatric disorders.

Common genetic variation appears to largely influence risk for neuropsychiatric disorders through effects on gene regulation. It is therefore possible to shed light on the biology of these conditions by testing for enrichment of associated genetic variation within regulatory genomic regions operating in specific tissues or cell types. Here, we have used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) to map open chromatin (an index of active regulatory genomic regions) in bulk tissue, NeuN+ and NeuN- nuclei from the prenatal human frontal cortex, and tested enrichment of single-nucleotide polymorphism (SNP) heritability for five neuropsychiatric disorders (autism spectrum disorder, attention deficit hyperactivity disorder [ADHD], bipolar disorder, major depressive disorder, and schizophrenia) within these regions. We observed significant enrichment of SNP heritability for ADHD, major depressive disorder, and schizophrenia within open chromatin regions (OCRs) mapped in bulk fetal frontal cortex, and for all five tested neuropsychiatric conditions when we restricted these sites to those overlapping histone modifications indicative of enhancers (H3K4me1) or promoters (H3K4me3) in fetal brain. SNP heritability for neuropsychiatric disorders was significantly enriched in OCRs identified in fetal frontal cortex NeuN- as well as NeuN+ nuclei overlapping fetal brain H3K4me1 or H3K4me3 sites. We additionally demonstrate the utility of our mapped OCRs for prioritizing potentially functional SNPs at genome-wide significant risk loci for neuropsychiatric disorders. Our data provide evidence for an early neurodevelopmental component to a range of neuropsychiatric conditions and highlight an important role for regulatory genomic regions active within both NeuN+ and NeuN- cells of the prenatal brain.

Novel epigenetic clock for fetal brain development predicts prenatal age for cellular stem cell models and derived neurons.

Induced pluripotent stem cells (iPSCs) and their differentiated neurons (iPSC-neurons) are a widely used cellular model in the research of the central nervous system. However, it is unknown how well they capture age-associated processes, particularly given that pluripotent cells are only present during the earliest stages of mammalian development. Epigenetic clocks utilize coordinated age-associated changes in DNA methylation to make predictions that correlate strongly with chronological age. It has been shown that the induction of pluripotency rejuvenates predicted epigenetic age. As existing clocks are not optimized for the study of brain development, we developed the fetal brain clock (FBC), a bespoke epigenetic clock trained in human prenatal brain samples in order to investigate more precisely the epigenetic age of iPSCs and iPSC-neurons. The FBC was tested in two independent validation cohorts across a total of 194 samples, confirming that the FBC outperforms other established epigenetic clocks in fetal brain cohorts. We applied the FBC to DNA methylation data from iPSCs and embryonic stem cells and their derived neuronal precursor cells and neurons, finding that these cell types are epigenetically characterized as having an early fetal age. Furthermore, while differentiation from iPSCs to neurons significantly increases epigenetic age, iPSC-neurons are still predicted as being fetal. Together our findings reiterate the need to better understand the limitations of existing epigenetic clocks for answering biological research questions and highlight a limitation of iPSC-neurons as a cellular model of age-related diseases.

Attenuated Induction of the Unfolded Protein Response in Adult Human Primary Astrocytes in Response to Recurrent Low Glucose.

Aims/hypothesis: Recurrent hypoglycaemia (RH) is a major side-effect of intensive insulin therapy for people with diabetes. Changes in hypoglycaemia sensing by the brain contribute to the development of impaired counterregulatory responses to and awareness of hypoglycaemia. Little is known about the intrinsic changes in human astrocytes in response to acute and recurrent low glucose (RLG) exposure. Methods: Human primary astrocytes (HPA) were exposed to zero, one, three or four bouts of low glucose (0.1 mmol/l) for three hours per day for four days to mimic RH. On the fourth day, DNA and RNA were collected. Differential gene expression and ontology analyses were performed using DESeq2 and GOseq, respectively. DNA methylation was assessed using the Infinium MethylationEPIC BeadChip platform. Results: 24 differentially expressed genes (DEGs) were detected (after correction for multiple comparisons). One bout of low glucose exposure had the largest effect on gene expression. Pathway analyses revealed that endoplasmic-reticulum (ER) stress-related genes such as HSPA5, XBP1, and MANF, involved in the unfolded protein response (UPR), were all significantly increased following low glucose (LG) exposure, which was diminished following RLG. There was little correlation between differentially methylated positions and changes in gene expression yet the number of bouts of LG exposure produced distinct methylation signatures. Conclusions/interpretation: These data suggest that exposure of human astrocytes to transient LG triggers activation of genes involved in the UPR linked to endoplasmic reticulum (ER) stress. Following RLG, the activation of UPR related genes was diminished, suggesting attenuated ER stress. This may be a consequence of a successful metabolic adaptation, as previously reported, that better preserves intracellular energy levels and a reduced necessity for the UPR.

DNA methylation meta-analysis reveals cellular alterations in psychosis and markers of treatment-resistant schizophrenia.

We performed a systematic analysis of blood DNA methylation profiles from 4483 participants from seven independent cohorts identifying differentially methylated positions (DMPs) associated with psychosis, schizophrenia, and treatment-resistant schizophrenia. Psychosis cases were characterized by significant differences in measures of blood cell proportions and elevated smoking exposure derived from the DNA methylation data, with the largest differences seen in treatment-resistant schizophrenia patients. We implemented a stringent pipeline to meta-analyze epigenome-wide association study (EWAS) results across datasets, identifying 95 DMPs associated with psychosis and 1048 DMPs associated with schizophrenia, with evidence of colocalization to regions nominated by genetic association studies of disease. Many schizophrenia-associated DNA methylation differences were only present in patients with treatment-resistant schizophrenia, potentially reflecting exposure to the atypical antipsychotic clozapine. Our results highlight how DNA methylation data can be leveraged to identify physiological (e.g., differential cell counts) and environmental (e.g., smoking) factors associated with psychosis and molecular biomarkers of treatment-resistant schizophrenia.

Genome-wide DNA methylation meta-analysis in the brains of suicide completers.

Suicide is the second leading cause of death globally among young people representing a significant global health burden. Although the molecular correlates of suicide remains poorly understood, it has been hypothesised that epigenomic processes may play a role. The objective of this study was to identify suicide-associated DNA methylation changes in the human brain by utilising previously published and unpublished methylomic datasets. We analysed prefrontal cortex (PFC, n = 211) and cerebellum (CER, n = 114) DNA methylation profiles from suicide completers and non-psychiatric, sudden-death controls, meta-analysing data from independent cohorts for each brain region separately. We report evidence for altered DNA methylation at several genetic loci in suicide cases compared to controls in both brain regions with suicide-associated differentially methylated positions enriched among functional pathways relevant to psychiatric phenotypes and suicidality, including nervous system development (PFC) and regulation of long-term synaptic depression (CER). In addition, we examined the functional consequences of variable DNA methylation within a PFC suicide-associated differentially methylated region (PSORS1C3 DMR) using a dual luciferase assay and examined expression of nearby genes. DNA methylation within this region was associated with decreased expression of firefly luciferase but was not associated with expression of nearby genes, PSORS1C3 and POU5F1. Our data suggest that suicide is associated with DNA methylation, offering novel insights into the molecular pathology associated with suicidality.

Genome-wide methylomic analysis in individuals with HNF1B intragenic mutation and 17q12 microdeletion.

Heterozygous mutation of the transcription factor HNF1B is the most common cause of monogenetic developmental renal disease. Disease-associated mutations fall into two categories: HNF1B intragenic mutations and a 1.3 Mb deletion at chromosome 17q12. An increase in neurodevelopmental disorders has been observed in individuals harbouring the 17q12 deletion but not in patients with HNF1B coding mutations.Previous investigations have concentrated on identifying a genetic cause for the increase in behavioural problems seen in 17q12 deletion carriers. We have taken the alternative approach of investigating the DNA methylation profile of these two HNF1B genotype groups along with controls matched for age, gender and diabetes status using the Illumina 450K DNA methylation array (total sample n = 60).We identified a number of differentially methylated probes (DMPs) that were associated with HNF1B-associated disease and passed our stringent experiment-wide significance threshold. These associations were largely driven by the deletion patients and the majority of the significant probes mapped to the 17q12 deletion locus. The observed changes in DNA methylation at this locus were not randomly dispersed and occurred in clusters, suggesting a regulatory mechanism reacting to haploinsufficiency across the entire deleted region.Along with these deletion-specific changes in DNA methylation, we also identified a shared DNA methylation signature in both mutation and deletion patient groups indicating that haploinsufficiency of HNF1B impacts on the methylome of a number of genes, giving further insight to the role of HNF1B.

Application of a novel molecular method to age free-living wild Bechstein's bats.

The age profile of populations fundamentally affects their conservation status. Yet, age is frequently difficult to assess in wild animals. Here, we assessed the use of DNA methylation of homologous genes to establish the age structure of a rare and elusive wild mammal: the Bechstein's bat (Myotis bechsteinii). We collected 62 wing punches from individuals whose ages were known as a result of a long-term banding study. DNA methylation was measured at seven CpG sites from three genes, which have previously shown age-associated changes in humans and laboratory mice. All CpG sites from the tested genes showed a significant relationship between DNA methylation and age, both individually and in combination (multiple linear regression R2  = 0.58, p < 0.001). Despite slight approximation around estimates, the approach is sufficiently precise to place animals into practically useful age cohorts. This method is of considerable practical benefit as it can reliably age individual bats. It is also much faster than traditional capture-mark-recapture techniques, with the potential to collect information on the age structure of an entire colony from a single sampling session to better inform conservation actions for Bechstein's bats. By identifying three genes where DNA methylation correlates with age across distantly related species, this study also suggests that the technique can potentially be applied across a wide range of mammals.

Mapping epigenetic changes to the host cell genome induced by Burkholderia pseudomallei reveals pathogen-specific and pathogen-generic signatures of infection.

The potential for epigenetic changes in host cells following microbial infection has been widely suggested, but few examples have been reported. We assessed genome-wide patterns of DNA methylation in human macrophage-like U937 cells following infection with Burkholderia pseudomallei, an intracellular bacterial pathogen and the causative agent of human melioidosis. Our analyses revealed significant changes in host cell DNA methylation, at multiple CpG sites in the host cell genome, following infection. Infection induced differentially methylated probes (iDMPs) showing the greatest changes in DNA methylation were found to be in the vicinity of genes involved in inflammatory responses, intracellular signalling, apoptosis and pathogen-induced signalling. A comparison of our data with reported methylome changes in cells infected with M. tuberculosis revealed commonality of differentially methylated genes, including genes involved in T cell responses (BCL11B, FOXO1, KIF13B, PAWR, SOX4, SYK), actin cytoskeleton organisation (ACTR3, CDC42BPA, DTNBP1, FERMT2, PRKCZ, RAC1), and cytokine production (FOXP1, IRF8, MR1). Overall our findings show that pathogenic-specific and pathogen-common changes in the methylome occur following infection.

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.

Epigenome-Wide DNA Methylation Analysis of Monozygotic Twins Discordant for Diurnal Preference.

Diurnal preference is an individual's preference for daily activities and sleep timing and is strongly correlated with the underlying circadian clock and the sleep-wake cycle validating its use as an indirect circadian measure in humans. Recent research has implicated DNA methylation as a mechanism involved in the regulation of the circadian clock system in humans and other mammals. In order to evaluate the extent of epigenetic differences associated with diurnal preference, we examined genome-wide patterns of DNA methylation in DNA from monozygotic (MZ) twin-pairs discordant for diurnal preference. MZ twins were selected from a longitudinal twin study designed to investigate the interplay of genetic and environmental factors in the development of emotional and behavioral difficulties. Fifteen pairs of MZ twins were identified in which one member scored considerably higher on the Horne-Ostberg Morningness-Eveningness Questionnaire (MEQ) than the other. Genome-wide DNA methylation patterns were assessed in twins' buccal cell DNA using the Illumina Infinium HumanMethylation450 BeadChips. Quality control and data pre-processing was undertaken using the wateRmelon package. Differentially methylated probes (DMPs) were identified using an analysis strategy taking into account both the significance and the magnitude of DNA methylation differences. Our data indicate that DNA methylation differences are detectable in MZ twins discordant for diurnal preference. Moreover, downstream gene ontology (GO) enrichment analysis on the top-ranked diurnal preference associated DMPs revealed significant enrichment of pathways that have been previously associated with circadian rhythm regulation, including cell adhesion processes and calcium ion binding.

Methylomic analysis of monozygotic twins discordant for childhood psychotic symptoms.

Childhood psychotic symptoms are associated with increased rates of schizophrenia, other psychiatric disorders, and suicide attempts in adulthood; thus, elucidating early risk indicators is crucial to target prevention efforts. There is considerable discordance for psychotic symptoms between monozygotic twins, indicating that child-specific non-genetic factors must be involved. Epigenetic processes may constitute one of these factors and have not yet been investigated in relation to childhood psychotic symptoms. Therefore, this study explored whether differences in DNA methylation at age 10 were associated with monozygotic twin discordance for psychotic symptoms at age 12. The Environmental Risk (E-Risk) Longitudinal Twin Study cohort of 2,232 children (1,116 twin pairs) was assessed for age-12 psychotic symptoms and 24 monozygotic twin pairs discordant for symptoms were identified for methylomic comparison. Children provided buccal samples at ages 5 and 10. DNA was bisulfite modified and DNA methylation was quantified using the Infinium HumanMethylation450 array. Differentially methylated positions (DMPs) associated with psychotic symptoms were subsequently tested in post-mortem prefrontal cortex tissue from adult schizophrenia patients and age-matched controls. Site-specific DNA methylation differences were observed at age 10 between monozygotic twins discordant for age-12 psychotic symptoms. Similar DMPs were not found at age 5. The top-ranked psychosis-associated DMP (cg23933044), located in the promoter of the C5ORF42 gene, was also hypomethylated in post-mortem prefrontal cortex brain tissue from schizophrenia patients compared to unaffected controls. These data tentatively suggest that epigenetic variation in peripheral tissue is associated with childhood psychotic symptoms and may indicate susceptibility to schizophrenia and other mental health problems.

Genome-wide methylomic analysis of monozygotic twins discordant for adolescent depression.

BACKGROUND: Adolescent depression is a common neuropsychiatric disorder that often continues into adulthood and is associated with a wide range of poor outcomes including suicide. Although numerous studies have looked at genetic markers associated with depression, the role of epigenetic variation remains relatively unexplored. METHODS: Monozygotic (MZ) twins were selected from an adolescent twin study designed to investigate the interplay of genetic and environmental factors in the development of emotional and behavioral difficulties. There were 18 pairs of MZ twins identified in which one member scored consistently higher (group mean within the clinically significant range) on self-rated depression than the other. We assessed genome-wide patterns of DNA methylation in twin buccal cell DNA using the Infinium HumanMethylation450 BeadChip from Illumina. Quality control and data preprocessing was undertaken using the wateRmelon package. Differentially methylated probes (DMPs) were identified using an analysis strategy taking into account both the significance and the magnitude of DNA methylation differences. The top differentially methylated DMP was successfully validated by bisulfite-pyrosequencing, and identified DMPs were tested in postmortem brain samples obtained from patients with major depressive disorder (n = 14) and matched control subjects (n = 15). RESULTS: Two reproducible depression-associated DMPs were identified, including the top-ranked DMP that was located within STK32C, which encodes a serine/threonine kinase, of unknown function. CONCLUSIONS: Our data indicate that DNA methylation differences are apparent in MZ twins discordant for adolescent depression and that some of the disease-associated variation observed in buccal cell DNA is mirrored in adult brain tissue obtained from individuals with clinical depression.