Childhood abuse is associated with methylation of multiple loci in adult DNA
BMC Medical Genomics volume 7, Article number: 13 (2014)
Childhood abuse is associated with increased adult disease risk, suggesting that processes acting over the long-term, such as epigenetic regulation of gene activity, may be involved. DNA methylation is a critical mechanism in epigenetic regulation. We aimed to establish whether childhood abuse was associated with adult DNA methylation profiles.
In 40 males from the 1958 British Birth Cohort we compared genome-wide promoter DNA methylation in blood taken at 45y for those with, versus those without, childhood abuse (n = 12 vs 28). We analysed the promoter methylation of over 20,000 genes and 489 microRNAs, using MeDIP (methylated DNA immunoprecipitation) in triplicate.
We found 997 differentially methylated gene promoters (311 hypermethylated and 686 hypomethylated) in association with childhood abuse and these promoters were enriched for genes involved in key cell signaling pathways related to transcriptional regulation and development. Using bisulfite-pyrosequencing, abuse-associated methylation (MeDIP) at the metalloproteinase gene, PM20D1, was validated and then replicated in an additional 27 males. Abuse-associated methylation was observed in 39 microRNAs; in 6 of these, the hypermethylated state was consistent with the hypomethylation of their downstream gene targets. Although distributed across the genome, the differentially methylated promoters associated with child abuse clustered in genome regions of at least one megabase. The observations for child abuse showed little overlap with methylation patterns associated with socioeconomic position.
Our observed genome-wide methylation profiles in adult DNA associated with childhood abuse justify the further exploration of epigenetic regulation as a mediating mechanism for long-term health outcomes.
Abuse in childhood, encompassing physical, sexual or emotional abuse, is a key component of a broader spectrum of child maltreatment . Life-long consequences of child abuse have been identified, including a greater risk of violence and delinquency, as well as adult depression and attempted suicide . Hazardous behaviors, such as smoking and alcoholism, have also been found to be associated with abuse in childhood [2–4] along with later disease risk factors, including obesity [1, 5], poorer immune function [1, 6–8] earlier menarche [9–11] and outcomes such as ischemic heart disease [6, 12, 13] and chronic obstructive lung disease [13, 14]. Explanations including biological mechanisms for long-term outcomes of child abuse have yet to be fully explored.
DNA methylation and histone modification play crucial roles in development, adaptation and response to environmental signals . Methylation of cytosine bases occurs at CpG sites and, in gene promoters, usually results in gene silencing, whereas loss of methylation is associated with activity. MicroRNAs that repress the expression of their often numerous target genes are also part of epigenetic regulation . MicroRNAs can down regulate key players in the epigenetic regulation machinery, but can also be silenced themselves by DNA methylation . Whilst epigenetic regulation, by definition, does not alter DNA sequence, DNA variants can influence methylation levels. However, DNA methylation associated with early adversity (prenatal famine) was found to be independent of that associated with genetic variation . Evidence to date suggests that stable changes in DNA methylation in the hippocampus of humans  and rats [20, 21] are triggered by maltreatment in early life.
Much DNA methylation is tissue specific  but most tissues are unavailable for population studies of living individuals. Given the multiple outcomes for childhood abuse, we hypothesize that DNA methylation associated with childhood abuse is system-wide . Several recent studies have supported the possibility of differential DNA methylation associations with social adversity in peripheral blood cells. For example, Borghol et al., demonstrated association of DNA methylation profiles with early life socioeconomic position in blood cells . Provencal et al., showed that differential maternal rearing is associated with differential DNA methylation profiles in both prefrontal cortex and blood T cells . Klengel et al., demonstrated childhood trauma-dependent DNA demethylation in functional glucocorticoid response elements of FKBP5 in blood cells . Mehta et al., have delineated recently DNA methylation signatures of child trauma and posttraumatic stress disorder in blood cells . Although blood cells turn over, they are derived from stem cells and progenitors that stay with us for a life long. Thus, it is plausible that a DNA methylation event in a stem cell population that is introduced in early life remains into adulthood.
We therefore aimed to establish whether childhood abuse is associated with adult gene promoter methylation in a genome-wide investigation of peripheral blood cells . We studied 40 adult males enrolled in the 1958 British Birth Cohort who have been found to have substantial variation in promoter methylation in over 6,000 genes, with a distinct methylation profile associated with socio-economic position . Those with childhood abuse in this cohort have been shown to have long-term associations with negative health outcomes, specifically, a greater prevalence of obesity among those who reported physical abuse in childhood .
All participants provided written consent and a blood sample for DNA analysis; ethical approval for a 45y biomedical survey and data analysis was given by the South-East Multi-Centre Research Ethics Committee(ref. 01/1/44) and the Joint UCL/UCLH Committees on the Ethnics of Human Research (Committee A) (ref. 08/H0714/40).
The selection of 40 adult males from the 1958 cohort  has been described previously  and are detailed in the Additional files. In brief, 17,638 participants were enrolled, all born in England, Scotland and Wales, during a single week in March 1958. At 45y, 4,177 males provided written consent and a blood sample for DNA analysis; ethical approval was given by the South-East Multi-Centre Research Ethics Committee. After exclusions (e.g. cancer or elevated C-reactive protein levels, immigrants), 3,362 white males were classified by socioeconomic position (SEP) and childhood abuse. Forty males were selected from extremes of SEP, including 12 who reported abuse (7 low and 5 high child SEP; 7 low and 5 high adult SEP). With exclusion of immigrants, the 1958 cohort shows little genetic population stratification .
Abuse was identified through participants' reports in a confidential questionnaire at 45y on the following experiences to age 16y:  "I was verbally abused by a parent";  "I suffered humiliation, ridicule, bullying or mental cruelty from a parent";  "I was physically abused by a parent –punched, kicked or hit or beaten with an object, or needed medical treatment";  "I was sexually abused by a parent". A report of any of these was scored as abuse. These questions were from the PATH Through Life Project including items derived from the Parental Bonding Instrument, the British National Survey of Health and Development and the US National Comorbidity Survey .
Measurement of relative DNA methylation levels
DNA sample preparation, methylated DNA immunoprecipitation (MeDIP) and microarray hybridization, scanning and data extraction were performed as described previously . Briefly, DNA was extracted from whole blood collected in EDTA at 45 years using an in-house, manual guanidine hydrochloride and ethanol precipitation method. DNA promoter methylation data from 20,533 genes and 489 microRNAs for the 40 participants were generated using MeDIP with an antibody that recognizes and binds 5-methylcytosine (DNA methylation) to isolate methylated DNA fragments. These fragments were then hybridized to custom-designed, high-density oligonucleotide microarrays, covering approximately 1000 bp upstream to 250 bp downstream at 100 bp spacing from the transcription start sites (TSS) in Ensembl (version 44). Microarray data files used in this study can be downloaded from the Gene Expression Omnibus (accession number: GSE31713). Three replicate microarrays were generated per individual and demonstrated adequate reproducibility . Both hierarchical clustering and principal components analysis applied to the 500 most variable probes across all microarrays showed that the three replicates clustered. Furthermore, >70% of the variance in these probes was explained by individual variation.
Microarray statistical analysis
The steps taken in the microarray statistical analysis are shown in Additional file 1: Figure S1 and justification for our approach is given in Additional files. Quality control involved generating MvA plots (i.e. plots of log(Cy5/Cy3) vs log(Cy5 × Cy3)) to identify those with severe dye biases or low signal. Microarrays deemed unacceptable were repeated, so no sample was excluded by quality control. Unsupervised clustering failed to identify batch effects related to hybridization date. Normalization of the final set of microarrays proceeded by computing log ratios of the bound (Cy5) and input (Cy3) microarray channel intensities for each microarray and then microarrays were normalized to one another using quantile normalization under the assumption that all samples have identical overall methylation levels. A probe was called differentially methylated if the modified t-statistic from ‘limma’  of Bioconductor  was significant (p < 0.05) and the log2 fold-difference of the mean group probe intensities was ≥0.25. A promoter was called differentially methylated if it contained a probe called differentially methylated and if it contained probes for which modified t-statistics were significantly higher or lower than the average probe on the microarray. Significance for the latter was determined by applying the Wilcoxon rank-sum test and then calculating a corresponding false discovery rate (FDR)  using the method of Benjamini and Hochberg . Promoters with FDR < 20% were called differentially methylated. This false discovery rate (FDR) was designed to test the chances of an overall false discovery among a series of related results. It is particularly useful for an exploratory analysis concerned with making general inferences from among a set of 'discoveries', rather than guarding against one or more individual false positives. The FDR threshold of 20% used here indicates that the expected proportion of promoters incorrectly called differentially methylated is around 20%. We find this threshold acceptable because this preliminary study is not meant to definitively characterise the epigenetic signatures of childhood abuse. In Figure 1, we present a heatmap showing probe methylation scores averaged across triplicate microarrays. Clustering was performed using Ward's hierarchical clustering algorithm with Pearson correlation distance as the distance metric.
All bioinformatic functional analysis was based on gene sets from GO , KEGG  and mSigDB . Enrichment for differential methylation was determined by applying the hypergeometric test to the overlap between known gene sets and those found in our study to be differentially methylated. FDR values were obtained by adjusting these significance levels over all gene sets and pathways considered. The differentially methylated gene set was then subjected to pathway analysis using Ingenuity Pathway Analysis software (http://ingenuity.com/products/pathways_analysis.html).
In assessing megabase regions of the genome, methylation patterns were obtained by computing the mean methylation score difference between abuse and non-abuse groups for each probe, generating a UCSC wiggle track file from these differences and then uploading it for display on the UCSC genome browser (http://genome.ucsc.edu/).
Validation and further methylation analysis
First, we validated the microarray calls, selecting 11 genes with the strongest methylation association with abuse (Additional file 2: Figure S2). Validation was performed using quantitative PCR (qPCR) of bound and input fractions of MeDIP with primers flanking the differentially methylated regions (Additional file 3: Table S1). Second, we validated two of these 11 genes, SLC17A3 and PM20D1, hypermethylated in association with abuse on MeDIP, by bisulfite pyrosequencing (in participants with sufficient DNA), as an independent method that measures methylation at specific sites . Next, bisulfite pyrosequencing analysis of PM20D1 was repeated on an additional 27 males selected using the same criteria as the original  group. Details of pyrosequencing conditions, including optimization of PCR amplification using 0, 50 and 100% methylation controls are provided in Additional files.
Cell type ratios in blood are known to fluctuate so certain methylation differences between individuals could be caused by different cell ratios, particularly in promoters of genes with cell-type specific methylation. To rule out this possibility in our analysis, we compared our results to published MeDIP , expression  and Illumina 450 K  profiles of purified blood cell types. In each published dataset, we identified all differentially methylated or expressed genes or probes (as appropriate) between all pairs of available blood cell type profiles and then compared those lists of differences to the list of differentially methylated genes or probes between the abused and non-abused individuals in our study. If variation in blood cell type ratios explains the methylation differences in our analysis, then we would expect to see at least one larger-than-expected intersection. In each case, however, hypergeometric tests failed to identify larger-than-expected intersections (p > 0.4 in each case). For the published MeDIP dataset , the microarray design used was similar to our design so we were able to construct a 1-1 mapping between over half of the probes across our respective designs. Probes were paired if they were closest and within 150 bp. Unfortunately, this MeDIP dataset only contained profiles for B and T cell purified cell types. We therefore expanded our analysis to include an expression dataset  with profiles for CD33+ (myeloid), CD34+, CD71+ (early erythroid), CD4+, CD8+, CD14+ (monocyte), CD19+ (B) and CD56 (natural killer) cells. We also included a recent Illumina 450 K dataset  with profiles for granulocytes, neutrophils, eosinophils, CD4+, CD8+, CD14+, CD19+ and CD56+ cells. For both these datasets, results were compared at the gene level.
Results and discussion
Physical, cognitive and emotional characteristics and biomarkers are listed for participants in Table 1. As expected, the abuse group showed more adverse characteristics than the non-abuse group, but differences did not reach conventional p-values in this small sample.
Hundreds of promoters are differentially methylated in association with child abuse
In total, 997 gene promoters were differentially methylated in association with childhood abuse, affecting 1141 different genes (Additional files). Of these promoters, 311 were hypermethylated and 686 were hypomethylated in abused compared to non-abused males. Figure 1 shows a heatmap depicting the relative methylation levels for all differentially methylated promoters and how they cluster within study participants. Even at more stringent thresholds (p < 0.01 and q < 0.05, see Methods), there were still 34 differentially methylated promoters corresponding to 58 different genes with similar proportions hypermethylated to hypomethylated. These cluster the study participants very similarly to the larger set of differentially methylated promoters (Additional file 4: Figure S3). To assess whether the broad methylation signature of childhood abuse was affected by the numerical imbalance of abused versus non-abused (N = 12 vs 28), we conducted a permutation analysis. We found that 997 differentially methylated promoters between abused and non-abused was larger than the number of differences associated with 82% of random partitions (410 of 500) of the participants with partition size ratios corresponding to 12 vs 28. To address any concern that the abuse associated methylation differences were reflecting differences in blood cell type ratios, we compared our results with recently published expression and methylation profiles of purified cell types [40–42]. We found no evidence of statistically significant overlaps between our results for abuse and cell-type specific methylation and expression patterns (p > 0.4, hypergeometric test; see Methods for details).
In 11 genes selected for validation, the direction of abuse associated methylation differences was confirmed using qPCR of bound and input MeDIP fractions (Additional file 2: Figure S2). We also confirmed abuse associated hypermethylation by pyrosequencing of sites in the promoter of SLC17A3 and the first exon of PM20D1 in the original samples (Figure 2A) and in an additional 27 males for PM20D1 (Figure 2B), and with SNP rs11540014 showing no association with methylation levels (data not shown). However, the associations in the promoter of SLC17A3 were not replicated in the additional 27 males.
Abuse-associated methylation clusters by biological function
Full results of functional analysis are given in Additional files. Differentially methylated gene promoters in abused males (1141 genes) were enriched in regulatory (169 genes) and developmental (230 genes) functions (Table 2). Central to both of these functions is the KEGG WNT signaling pathway; enriched for genes  for which promoters are hypomethylated in abused individuals, consistent with activation of this pathway in blood cells of the abuse group (Figure 3). No other KEGG pathway was enriched with differentially methylated genes at p <0.05 (uncorrected for multiple testing). Of the differentially methylated genes that perform some regulatory function, most (134 of 169) are hypomethylated in abused males. The regulation mainly affects transcription as indicated by enrichment of these genes in functional categories such as chromatin modification (28 genes), histone modification (11 genes) and transcription factor binding (35 genes). Similarly, most of the 230 developmental genes are hypomethylated in abused males (172 genes), best characterized by the general gene ontology category “multicellular organismal development” (163 genes). More specific subcategories do not show significant enrichment.
Differentially hypermethylated gene promoters in abused males are enriched in few functional categories. One of these, "cell surface receptor linked signal transduction", contains 125 genes with differentially methylated promoters of which 46 are hypermethylated in abused individuals. An Ingenuity functional analysis of the differentially methylated genes revealed similar molecular and cellular functions associated with transcriptional control (Additional file 5: Figure S4).
Abuse-associated methylation consistent with microRNA targeting
MicroRNA genes are, like DNA methylation, known to repress the expression of target genes. However, unlike an individual methylation mark which typically targets a single nearby gene, each microRNA is associated with a specific set of a few hundred target genes . We discovered an association of microRNA DNA hypermethylation with abuse. Of 489 microRNAs analysed, 39 were differentially methylated, of which 31 were hypermethylated in association with abuse. The target genes of six of these included a highly non-random proportion of genes with decreased promoter methylation in abused males (Table 3).
Abuse-associated hypomethylation and CpG density
DNA methylation in regions of relatively high CpG frequency, known as CpG islands, plays an important regulatory role in the otherwise CpG-depleted (≤40% of that expected) mammalian genome [44, 45]. In spite of the fact that MeDIP is known to enrich for methylation differences away from CpG islands , we observed unusually high CpG frequencies in promoters with reduced methylation levels in abused individuals. This frequency (0.86) is significantly higher than that observed in the average promoter (frequency = 0.42; p < 1.4 × 10-285) as well as promoters with increased methylation levels in abused individuals (frequency = 0.38; p < 4 × 10-138; Additional file 6: Figure S5). This frequency (0.86) is even higher than the 0.6 threshold used to define CpG islands.
Abuse-associated methylation clusters by genomic location
Differentially methylated DNA loci associated with early life environments tend to cluster in the genome [24, 47]. Chromosome-wide views of our data reveal megabase-sized regions significantly enriched for differentially methylated promoters (Figure 4). At the chromosomal level, chromosomes 16 and 17 were significantly enriched for hypomethylated promoters in abused individuals, whereas chromosome X was significantly enriched for hypermethylated promoters. At the megabase level, three regions were significantly enriched for differentially methylated promoters (p < 0.05). All were hypermethylated in abused individuals: chr1:246250000-247000000, chr14:100250000-101000000 and chr19:58500000-59250000 (genome assembly hg18), but only the regions on chromosomes 1 and 19 passed multiple testing correction with FDR below 0.2 (FDR < 0.006 and 0.0001, respectively; Figure 4D,E). The regions on chromosomes 14 and 19 each contain a cluster of microRNAs in which promoters account for all of their statistically significant site-specific differential methylation.
Clustering of differential promoter methylation, up to 2 Mb apart, was detectable across the entire genome (Additional file 7: Figure S6).
Socio-economic position (SEP) and abuse
Previously, we identified 1252 gene promoters associated with childhood SEP and 545 associated with adulthood SEP . Only 73 of 1252 (5.8%) and 19 of 545 (3.5%) gene promoters were also differentially methylated in association with childhood abuse. Just three (CTAGE5, GNG4, MYO1B) were differentially methylated in association with all three characteristics (childhood and adult SEP and childhood abuse). The association for PM20D1 was specific to abuse.
Blood DNA of 45y old males revealed differentially methylated gene promoters associated with abuse that occurred three decades earlier in childhood. There were several novel findings from our study. First, hundreds of specific promoter associations were uncovered, with approximately two-thirds hypomethylated in the abused group. Second, replication confirmed that hypermethylation in PM20D1 is associated with childhood abuse. Third, microRNA gene targets tended to be hypomethylated, particularly when the microRNA itself was hypermethylated. Fourth, differentially methylated genes were clustered in discrete functional pathways and in genomic locations. These findings support the hypothesis that the differences in DNA methylation we observed were non-random and reflect an organized biological process.
It is now known that genes act through functional and interacting pathways, so we adopted a genome-wide approach to DNA methylation analysis, recognizing that modest epigenetic changes in numerous genes could reset the function of gene networks having phenotypic effects. We found enrichment of differentially methylated promoters in the WNT signaling pathway complex with hypomethylation of 15 genes in the abused. Elsewhere this pathway complex has been found to play a key role in embryonic development and cellular proliferation , and is deregulated in some chronic health conditions such as obesity [49–51], diabetes [52–54], metabolic syndrome , cancer [56–59] and inflammatory processes [56–58]. Whilst recognizing that our findings do not provide evidence for causal links between child abuse and later outcomes, they nonetheless raise the prospect of mediation by epigenetic modifications.
Of particular note was hypermethylation of PM20D1 in association with abuse, given a previous study showing a variably methylated region at this metalloproteinase gene was hypermethylated in association with obesity . This association persisted over 10 years of follow-up in an elderly population. Interestingly, child abuse has been shown to be associated with adult obesity in the full 1958 cohort  and is suggested by our Table 1. It is perhaps surprising to note that both our association with childhood abuse and the association with obesity were observed in blood DNA when PM20D1 has its highest expression levels in the brain and lowest expression levels in blood . Given that it is highly conserved from yeast to human, it likely plays a key though little understood role in the cell. By contrast, SLC17A3 is like most of the genes differentially methylated in childhood abuse, most highly expressed in blood and a few specific brain regions (hypothalamus, prefrontal cortex, pituitary) . It appears conserved in fewer species, mainly the higher mammals, and the expressed protein acts as a voltage-driven transporter in blood. Given this basic role, it is likely essential at nearly all stages of life.
Further support for epigenetic regulation working through interacting pathways comes from the striking enrichment in the abuse group of hypermethylated microRNAs combined with hypomethylation of their respective gene targets across the genome. It implies that during typical development, active transcription of these microRNAs is combined with synergistic target methylation to create a double layer of repression of these target genes; a repression that is lifted in association with child abuse.
Intriguingly, hypomethylated gene promoters in abused individuals typically contained sequences with very high CpG frequency. Demethylation of such CpG-rich promoters in abused males suggests that abuse leads to increased activity of key basic cellular functions, such as gene regulation and development, as found in pathway analysis. Another genomic feature associated with abuse was the clustering of differential promoter methylation detectable across the entire genome, providing further evidence of genome-wide as well as gene-specific organization of epigenetic profiles.
Previously we observed genome-wide clustering in association with SEP, but importantly, the “methylation signature” for abuse differed, such that <10% of the differentially methylated regions overlapped with childhood SEP . Also, the differentially methylated genes were enriched in different functional pathways, notably, MAP kinase for SEP and WNT for child abuse. Further, the abuse associated differential methylation of microRNAs and their target genes was not seen for SEP. Whilst not ruling out generic associations with early life adversities, our findings suggest that different adversities are associated with different epigenetic changes to the genome.
Several methodological considerations arise here. First, reliable measurement of the frequency and severity of child abuse is not straightforward1. Child abuse was identified through participant’s report at 45y and was primarily emotional and physical abuse – only rarely sexual abuse. All measures have biases and inconsistencies yet retrospective reports are an accepted method of ascertainment in population studies . Furthermore, prospective identification of abuse is not feasible in large studies and likely to be unrepresentative. By contrast, retrospective self-report, used here, is feasible though it is likely to underestimate true levels of abuse. Second, given the scale of assessing methylation at all promoters, we could only study a small but selected sample. Whilst our study is imbalanced with respect to abuse (12 vs 28) it has the benefit of control for SEP. Third, we used DNA from whole blood to test our hypothesis, currently the only practical option for population based studies. We cannot know the extent to which our results relate to gene expression. Use of whole blood also raises the possibility that abuse-associated differences in B-to-T cell ratios might account for some of our observations. We have partly addressed this by noting that B-cell and T-cell expression and methylation profiles [40, 41] do not differ for many genes with abuse-associated methylation levels. Fourth, those abused in childhood might represent a distinct genetic group, but genetic differences alone are unlikely to account for all methylation differences observed here. Given the possibility of differences in epigenetic response due to genetic variation, future integrated studies of the epigenome and whole genome sequencing are an important next step. Fifth, our study is imbalanced including 28 controls compared to only 12 with childhood abuse resulting in reduced power to identify methylation differences. Nonetheless, this preliminary study was able to discover hundreds of differentially methylated promoters so future studies with better balance are likely discover many more. Finally, there is currently no ‘gold standard’ for measuring the methylome, yet MeDIP is a well-established genome-wide method that has been evaluated [46, 61–65] and we confirmed all the micro-array calls in the top 11 methylation differences. Current genome-wide methods are more complementary than interchangeable and each has its strengths and weaknesses. Our analyses included triplicate arrays and methylation differences were confirmed in selected genes using other gene-specific methods both here and previously . In using an analytic approach that was sensitive to subtle methylation associations across gene networks necessarily results in some false positives (for justification see Additional files). However, the non-random organization of methylation differences throughout the genome supports our main hypothesis that childhood abuse is associated with DNA methylation changes in adult blood.
In sum, the pattern of changes associated with child abuse detected in peripheral blood cells of 45 year-olds suggest that there is a system-wide readjustment of the epigenome to signals triggered by early life abuse. Our study does not demonstrate causality, nor can it demonstrate a temporal relationship between child abuse and DNA methylation levels in adulthood. It does, however, provide a justification for a range of studies addressing epigenetic responses to child abuse and their mediating role with later phenotypic outcomes.
Gilbert R, Widom CS, Browne K, Fergusson D, Webb E, Janson S: Burden and consequences of child maltreatment in high-income countries. Lancet. 2009, 373 (9657): 68-81. 10.1016/S0140-6736(08)61706-7.
Anda RF, Croft JB, Felitti VJ, Nordenberg D, Giles WH, Williamson DF, Giovino GA: Adverse childhood experiences and smoking during adolescence and adulthood. JAMA. 1999, 282 (17): 1652-1658. 10.1001/jama.282.17.1652.
Anda RF, Whitfield CL, Felitti VJ, Chapman D, Edwards VJ, Dube SR, Williamson DF: Adverse childhood experiences, alcoholic parents, and later risk of alcoholism and depression. Psychiatr Serv. 2002, 53 (8): 1001-1009.
Dube SR, Miller JW, Brown DW, Giles WH, Felitti VJ, Dong M, Anda RF: Adverse childhood experiences and the association with ever using alcohol and initiating alcohol use during adolescence. J Adolesc Health. 2006, 38 (4): 444-e441-410
Felitti VJ, Anda RF, Nordenberg D, Williamson DF, Spitz AM, Edwards V, Koss MP, Marks JS: Relationship of childhood abuse and household dysfunction to many of the leading causes of death in adults. The Adverse Childhood Experiences (ACE) study. Am J Prev Med. 1998, 14 (4): 245-258. 10.1016/S0749-3797(98)00017-8.
Dong M, Giles WH, Felitti VJ, Dube SR, Williams JE, Chapman DP, Anda RF: Insights into causal pathways for ischemic heart disease: adverse childhood experiences study. Circulation. 2004, 110 (13): 1761-1766. 10.1161/01.CIR.0000143074.54995.7F.
Lehman BJ, Taylor SE, Kiefe CI, Seeman TE: Relation of childhood socioeconomic status and family environment to adult metabolic functioning in the CARDIA study. Psychosom Med. 2005, 67 (6): 846-854. 10.1097/01.psy.0000188443.48405.eb.
Sumanen M, Koskenvuo M, Sillanmaki L, Mattila K: Childhood adversities experienced by working-aged coronary heart disease patients. J Psychosom Res. 2005, 59 (5): 331-335. 10.1016/j.jpsychores.2005.04.005.
Wise LA, Palmer JR, Rothman EF, Rosenberg L: Childhood abuse and early menarche: findings from the black women's health study. Am J Publ Health. 2009, 99 (Suppl 2): S460-S466.
Zabin LS, Emerson MR, Rowland DL: Childhood sexual abuse and early menarche: the direction of their relationship and its implications. J Adolesc Health. 2005, 36 (5): 393-400. 10.1016/j.jadohealth.2004.07.013.
Vigil JM, Geary DC, Byrd-Craven J: A life history assessment of early childhood sexual abuse in women. Dev Psychol. 2005, 41 (3): 553-561.
Fuller-Thomson E, Brennenstuhl S, Frank J: The association between childhood physical abuse and heart disease in adulthood: findings from a representative community sample. Child Abuse Negl. 2010, 34 (9): 689-698. 10.1016/j.chiabu.2010.02.005.
Goodwin RD, Stein MB: Association between childhood trauma and physical disorders among adults in the United States. Psychol Med. 2004, 34 (3): 509-520. 10.1017/S003329170300134X.
Springer KW: Childhood physical abuse and midlife physical health: testing a multi-pathway life course model. Soc Sci Med. 2009, 69 (1): 138-146. 10.1016/j.socscimed.2009.04.011.
Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33 (Suppl): 245-254.
Chuang JC, Jones PA: Epigenetics and microRNAs. Pediatr Res. 2007, 61 (5 Pt 2): 24R-29R.
Baer C, Claus R, Frenzel LP, Zucknick M, Park YJ, Gu L, Weichenhan D, Fischer M, Pallasch CP, Herpel E, Rehli M, Byrd JC, Wendtner CM, Plass C: Extensive promoter DNA hypermethylation and hypomethylation is associated with aberrant microRNA expression in chronic lymphocytic leukemia. Canc Res. 2012, 72 (15): 3775-3785. 10.1158/0008-5472.CAN-12-0803.
Tobi EW, Slagboom PE, van Dongen J, Kremer D, Stein AD, Putter H, Heijmans BT, Lumey LH: Prenatal famine and genetic variation are independently and additively associated with DNA methylation at regulatory loci within IGF2/H19. PloS One. 2012, 7 (5): e37933-10.1371/journal.pone.0037933.
McGowan PO, Sasaki A, D'Alessio AC, Dymov S, Labonte B, Szyf M, Turecki G, Meaney MJ: Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat Neurosci. 2009, 12 (3): 342-348. 10.1038/nn.2270.
Roth TL, Lubin FD, Funk AJ, Sweatt JD: Lasting epigenetic influence of early-life adversity on the BDNF gene. Biol Psychiatr. 2009, 65 (9): 760-769. 10.1016/j.biopsych.2008.11.028.
Weaver IC, Cervoni N, Champagne FA, D'Alessio AC, Sharma S, Seckl JR, Dymov S, Szyf M, Meaney MJ: Epigenetic programming by maternal behavior. Nat Neurosci. 2004, 7 (8): 847-854. 10.1038/nn1276.
Razin A, Szyf M: DNA methylation patterns. Formation and function. Biochim Biophys Acta. 1984, 782 (4): 331-342. 10.1016/0167-4781(84)90043-5.
Szyf M: How do environments talk to genes?. Nat Neurosci. 2013, 16 (1): 2-4.
Borghol N, Suderman M, McArdle W, Racine A, Hallett M, Pembrey M, Hertzman C, Power C, Szyf M: Associations with early-life socio-economic position in adult DNA methylation. Int J Epidemiol. 2012, 41 (1): 62-74. 10.1093/ije/dyr147.
Provencal N, Suderman MJ, Guillemin C, Massart R, Ruggiero A, Wang D, Bennett AJ, Pierre PJ, Friedman DP, Cote SM, Hallett M, Tremblay RE, Suomi SJ, Szyf M: The signature of maternal rearing in the methylome in rhesus macaque prefrontal cortex and T cells. J Neurosci Offic J Soc Neurosci. 2012, 32 (44): 15626-15642. 10.1523/JNEUROSCI.1470-12.2012.
Klengel T, Mehta D, Anacker C, Rex-Haffner M, Pruessner JC, Pariante CM, Pace TW, Mercer KB, Mayberg HS, Bradley B, Nemeroff CB, Holsboer F, Heim CM, Ressler KJ, Rein T, Binder EB: Allele-specific FKBP5 DNA demethylation mediates gene-childhood trauma interactions. Nat Neurosci. 2013, 16 (1): 33-41.
Mehta D, Klengel T, Conneely KN, Smith AK, Altmann A, Pace TW, Rex-Haffner M, Loeschner A, Gonik M, Mercer KB, Bradley B, Muller-Myhsok B, Ressler KJ, Binder EB: Childhood maltreatment is associated with distinct genomic and epigenetic profiles in posttraumatic stress disorder. Proc Natl Acad Sci U S A. 2013, 110 (20): 8302-8307. 10.1073/pnas.1217750110.
Thomas C, Hypponen E, Power C: Obesity and type 2 diabetes risk in midadult life: the role of childhood adversity. Pediatrics. 2008, 121 (5): e1240-e1249. 10.1542/peds.2007-2403.
Power C, Elliott J: Cohort profile: 1958 British birth cohort (National Child Development Study). Int J Epidemiol. 2006, 35 (1): 34-41.
Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678. 10.1038/nature05911.
Rosenman S, Rodgers B: Childhood adversity in an Australian population. Soc Psychiatry Psychiatr Epidemiol. 2004, 39 (9): 695-702. 10.1007/s00127-004-0802-0.
Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor, Volume 1. Edited by: Gentleman VC R, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Soric B: Statistical "discoveries" and effect-size estimation. J Am Stat Assoc. 1989, 84 (406): 608-610.
Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995, 57 (1): 289-300.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Colella S, Shen L, Baggerly KA, Issa JP, Krahe R: Sensitive and quantitative universal Pyrosequencing methylation analysis of CpG sites. Biotechniques. 2003, 35 (1): 146-150.
Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Graf S, Tomazou EM, Backdahl L, Johnson N, Herberth M, Howe KL, Jackson DK, Miretti MM, Fiegler H, Marioni JC, Birney E, Hubbard TJ, Carter NP, Tavare S, Beck S: An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res. 2008, 18 (9): 1518-1529. 10.1101/gr.077479.108.
Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004, 101 (16): 6062-6067. 10.1073/pnas.0400782101.
Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlen SE, Greco D, Soderhall C, Scheynius A, Kere J: Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PloS One. 2012, 7 (7): e41361-10.1371/journal.pone.0041361.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.
Bird AP: CpG-rich islands and the function of DNA methylation. Nature. 1986, 321 (6067): 209-213. 10.1038/321209a0.
Aissani B, Bernardi G: CpG islands: features and distribution in the genomes of vertebrates. Gene. 1991, 106 (2): 173-183. 10.1016/0378-1119(91)90197-J.
Nair SS, Coolen MW, Stirzaker C, Song JZ, Statham AL, Strbenac D, Robinson MW, Clark SJ: Comparison of methyl-DNA immunoprecipitation (MeDIP) and methyl-CpG binding domain (MBD) protein capture for genome-wide DNA methylation analysis reveal CpG sequence coverage bias. Epigenetics. 2011, 6 (1): 34-44. 10.4161/epi.6.1.13313.
McGowan PO, Suderman M, Sasaki A, Huang TC, Hallett M, Meaney MJ, Szyf M: Broad epigenetic signature of maternal care in the brain of adult rats. PloS One. 2011, 6 (2): e14739-10.1371/journal.pone.0014739.
Clevers H, van de Wetering M: TCF/LEF factor earn their wings. Trends Genet. 1997, 13 (12): 485-489. 10.1016/S0168-9525(97)01305-X.
Oh DY, Olefsky JM: Medicine. Wnt fans the flames in obesity. Science. 329 (5990): 397-398.
Ouchi N, Higuchi A, Ohashi K, Oshima Y, Gokce N, Shibata R, Akasaki Y, Shimono A, Walsh K: Sfrp5 is an anti-inflammatory adipokine that modulates metabolic dysfunction in obesity. Science. 2010, 329 (5990): 454-457. 10.1126/science.1188280.
Prestwich TC, Macdougald OA: Wnt/beta-catenin signaling in adipogenesis and metabolism. Curr Opin Cell Biol. 2007, 19 (6): 612-617. 10.1016/j.ceb.2007.09.014.
Hattersley AT: Prime suspect: the TCF7L2 gene and type 2 diabetes risk. J Clin Invest. 2007, 117 (8): 2077-2079. 10.1172/JCI33077.
Cauchi S, Froguel P: TCF7L2 genetic defect and type 2 diabetes. Curr Diab Rep. 2008, 8 (2): 149-155. 10.1007/s11892-008-0026-x.
Jin T: The WNT, signalling pathway and diabetes mellitus. Diabetologia. 2008, 51 (10): 1771-1780. 10.1007/s00125-008-1084-y.
Schinner S: Wnt-signalling and the metabolic syndrome. Horm Metab Res. 2009, 41 (2): 159-163. 10.1055/s-0028-1119408.
Pereira CP, Bachli EB, Schoedon G: The wnt pathway: a macrophage effector molecule that triggers inflammation. Curr Atheroscler Rep. 2009, 11 (3): 236-242. 10.1007/s11883-009-0036-4.
Schett G, Zwerina J, David JP: The role of Wnt proteins in arthritis. Nat Clin Pract Rheumatol. 2008, 4 (9): 473-480.
Takahashi-Yanaga F, Sasaguri T: The Wnt/beta-catenin signaling pathway as a target in drug discovery. J Pharmacol Sci. 2007, 104 (4): 293-302. 10.1254/jphs.CR0070024.
Waltzer L, Bienz M: The control of beta-catenin and TCF during embryonic development and cancer. Canc Metastasis Rev. 1999, 18 (2): 231-246.
Feinberg AP, Irizarry RA, Fradin D, Aryee MJ, Murakami P, Aspelund T, Eiriksdottir G, Harris TB, Launer L, Gudnason V, Fallin MD: Personalized epigenomic signatures that are stable over time and covary with body mass index. Sci Transl Med. 2010, 2 (49): 49ra67.
Irizarry RA, Ladd-Acosta C, Carvalho B, Wu H, Brandenburg SA, Jeddeloh JA, Wen B, Feinberg AP: Comprehensive high-throughput arrays for relative methylation (CHARM). Genome Res. 2008, 18 (5): 780-790. 10.1101/gr.7301508.
Jia J, Pekowska A, Jaeger S, Benoukraf T, Ferrier P, Spicuglia S: Assessing the efficiency and significance of Methylated DNA Immunoprecipitation (MeDIP) assays in using in vitro methylated genomic DNA. BMC Res Notes. 2010, 3: 240-10.1186/1756-0500-3-240.
Jin SG, Kadam S, Pfeifer GP: Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine. Nucleic Acids Res. 2010, 38 (11): e125-10.1093/nar/gkq223.
Rajendram R, Ferreira JC, Grafodatskaya D, Choufani S, Chiang T, Pu S, Butcher DT, Wodak SJ, Weksberg R: Assessment of methylation level prediction accuracy in methyl-DNA immunoprecipitation and sodium bisulfite based microarray platforms. Epigenetics. 2011, 6 (4): 410-415. 10.4161/epi.6.4.14763.
Robinson MD, Stirzaker C, Statham AL, Coolen MW, Song JZ, Nair SS, Strbenac D, Speed TP, Clark SJ: Evaluation of affinity-based genome-wide DNA methylation data: effects of CpG density, amplification bias, and copy number variation. Genome Res. 2010, 20 (12): 1719-1729. 10.1101/gr.110601.110.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/13/prepub
We dedicate this work to the memory of our recently deceased co-author Clyde Hertzman, who played a central role in conceiving of this study and inspiring its completion. Grants from the Canadian Institute of Health Research; MOP-42411 the Sackler program in Psychobiology and Epigenetics at McGill University (to M.S.); fellowship from the Canadian Institute of Health Research and the Fragile X Research Foundation of Canada (to N.B.); the UK Medical Research Council (MRC) (grant G0000934 to the clinical examination and DNA banking of the 1958 cohort); and funding received for the MRC Centre of Epidemiology for Child Health (Grant G0400546) (for SPP). Work undertaken at Great Ormond Street Hospital/University College London, Institute of Child Health is in part supported by funding from the Department of Health’s National Institute of Health Research (‘Biomedical Research Centres’ funding). Laboratory work in the Department of Social Medicine, University of Bristol was supported in part by a grant to MP from the Welcome Trust.
The authors declare that they have no competing interests.
Study was designed by CP, CH and MS. Participants were selected by CP and CH. Methylation analysis was completed by NB and JP. Bioinformatic analysis was performed by MS. The entire process was overseen by MS. JP, NB, CH, MP, SPP, CP, MS and MS all contributed to writing the manuscript. All authors read and approved the final manuscript.
Matthew Suderman, Nada Borghol, Jane J Pappas contributed equally to this work.
Electronic supplementary material
Additional file 2: Figure S2: Validation by qPCR. Eleven gene promoters identified by microarray as being differentially methylated were subjected to real-time PCR quantification of the enrichment by the MeDIP procedure. Results were normalized against a methylated luciferase gene-containing plasmid (control), which was added to every sample in equal quantity before MeDIP. The y-axis represents relative concentration levels generated by applying PCR to methylation-enriched DNA. Each real-time PCR reaction was performed in duplicate for all subjects. Shown are the averages per group. Error bars indicate the standard error of the mean. Above the chart are tracks of the regions with the microarray data. (The bars indicate the difference between the abuse and the non-abuse groups, bars descending from the physical map are regions that are more methylated in the abused than the non-abused group; lower tracks identify probes with the most statistically significant differences). Primers for each PCR are given in Additional file 3: Table S1. They were selected so that the forward primer (denote by ‘F’) binds to the left and the reverse primer (denoted by ‘R’) binds to the right of the most significantly different probe. In some cases, two sets of PCR primers were designed, denoted by ‘set1’ and ‘set2’. 85% of the eleven gene promoters show statistically significant PCR quantification differences (*: P<0.05; **: P<0.01), hence validating differences found by microarray. (ZIP 170 KB)
Additional file 4: Figure S3: Promoter methylation associated with childhood abuse. Heatmap showing MeDIP probe values from the 34 differentially methylated promoters (rows) across all 40 participants (columns) based on more stringent thresholds (q < 0.05 and p < 0.01, see Methods). Each promoter is represented by the probe most associated with childhood abuse. Blackened squares above the columns denote non-abuse males, white squares denote those with childhood abuse. Other covariates included are childhood and adulthood socio-economic position (white = low, gray = high). None appears to explain the main sample clusters. (TIFF 16 MB)
Additional file 5: Figure S4: Summary of functional analysis. Genes with hypermethylated or hypomethylated promoters in the abuse group were analysed by Ingenuity Pathway Analysis®. Gene categories enriched with this set of genes as well as enrichment p-values are listed. (TIFF 16 MB)
Additional file 6: Figure S5: CpG frequency in differentially methylated regions. Bars indicate average normalized CpG frequencies (observed/expected CpG frequency) of ‘all’ genomic regions profiled, regions ‘hypermethylated’ in abused individuals and regions ‘hypomethylated’ in abused individuals. Error bars depict standard deviation. The dashed line indicates the usual CpG frequency used to identify CpG islands. (TIFF 16 MB)
Additional file 7: Figure S6: Methylation dependencies across megabases. Shown are correlations of methylation differences from 500 kilobase regions at various distances apart. The level of clustering was quantified as the level of correlation between the differential methylation statistics within promoters at different distances apart. The solid grey region contains the 95% CI, and error bars contain the 95% CI for correlation values. (TIFF 16 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Suderman, M., Borghol, N., Pappas, J.J. et al. Childhood abuse is associated with methylation of multiple loci in adult DNA. BMC Med Genomics 7, 13 (2014). https://doi.org/10.1186/1755-8794-7-13
- Childhood abuse
- Early life environment
- DNA methylation