|Home | About | Journals | Submit | Contact Us | Français|
Fifteen-year-old adolescents (N=109) in a longitudinal study of child development were recruited to examine differences in DNA methylation in relation to parent reports of adversity during the adolescents’ infancy and preschool periods. Microarray technology applied to 28,000 cytosine-guanine dinucleotide (CpG) sites within DNA derived from buccal epithelial cells showed differential methylation among adolescents whose parents reported high levels of stress during their children’s early lives. Maternal stressors in infancy and paternal stressors in the preschool years were most strongly predictive of differential methylation, and the patterning of such epigenetic marks varied by children’s gender. To the authors’ knowledge, this is the first report of prospective associations between adversities in early childhood and the epigenetic conformation of adolescents’ genomic DNA.
Experience “gets under the skin” early in life, and does so in ways that affect the course of human development. Heart disease, diabetes, obesity, depression, substance abuse, school success, premature mortality, disability at retirement, and accelerated aging and memory loss all have determinants in early life (Harkonmaki, et al., 2007). In particular, childhood psychosocial adversity affects multiple outcomes, including smoking, suicide, depression, obesity, illicit drug use, alcoholism, teen pregnancy, sexual risk behaviors and sexually transmitted diseases (Felitti, et al., 1998; Thomas, Hypponen, & Power, 2008). Although physical and emotional abuse and other dramatic adverse childhood events have health and developmental consequences (Heim, Newport, Mletzko, Miller, & Nemeroff, 2008), it is often the less memorable yet prevalent, co-occurring, day-to-day misfortunes of early childhood (e.g., impoverished parent-child interactions) that have a lasting influence on the subsequent life course (Kishiyama, Boyce, Jimenez, Perry, & Knight, 2009).
The biological plausibility of this scenario is well established. Experience plays a crucial role in the early phases of conception, prenatal, and post-natal periods of children's development (Swain, Lorberbaum, Kose, & Strathearn, 2007; Talge, Neal, & Glover, 2007). Sensitive periods in brain and biological development start in the prenatal period, reach a peak in the first few years of life, and continue at a declining rate throughout childhood and adolescence, such that the brain does not reach its "adult" state until the end of the second or beginning of the third decade of life (Fox, Levitt, & Nelson, 2010). Early sensory stimulation activates specific genes in different parts of the brain to differentiate neuronal function and establish sensory pathways. Sensory pathways, in turn, influence the development of neural pathways to other parts of the brain involved in coping, movement, language, cognition, and biological pathways, including the immune and hormone systems (McEwen, 2008). For example, key executive functions, regarding how an individual responds to social and emotional stimuli, develop in the prefrontal cortex from approximately age three to nine years, whereas neural connections to the prefrontal cortex, from centers in the midbrain that sense environmental threats, develop earlier (Gunnar & Nelson, 1994; Knudsen, Heckman, Cameron, & Shonkoff, 2006). Most importantly, early environments are mediated through relationships with primary caregivers. A prime example is infant attachment to a parental figure (usually the mother) in both higher primates and humans. Attachment, in turn, drives the development of neural pathways that help the developing child's brain become attuned to the immediate environment (Bowlby, 1969; Sroufe, 2005).
There has been substantial theoretical work on the multiple ecological factors associated with child development (Bronfenbrenner & Morris, 2006) and seminal empirical work on the cumulative risks these factors can pose (Rutter, 1979; Sameroff, 1998). Additional longitudinal studies have demonstrated the deleterious effects of accumulated risk on child outcomes (e.g., (Evans, Kim, Ting, Tesher, & Shannis, 2007)). Other researchers have focused on the child and adolescent sequelae of specific aspects of stressful family environments. For example, parental depression (Goodman & Gotlib, 1999; Halligan, Murray, Martins, & Cooper, 2007), marital conflict (Cummings & Davies, 2002), and parenting stress (Bayer, Hiscock, Ukoumunne, Price, & Wake, 2008) have been shown individually and jointly to predict child health and functioning. Research also suggests that the effects of children’s exposures to parents’ stress may vary depending on the timing of exposure (Essex, Klein, Cho, & Kalin, 2002; Lupien, McEwen, Gunnar, & Heim, 2009). The research base further leads to an expectation that mothers will contribute in different ways – and disproportionately at earlier stages – than fathers, although each may be expected to contribute to the child's development (Bornstein, 2002; Keller, Cummings, Peterson, & Davies, 2009; Ramchandani, et al., 2008). Furthermore, fathers and mothers may influence boys' and girls' development differently and to different degrees (Morrell & Murray, 2003).
In our own previous work, we have found that children’s exposure to parental stress factors, including parental depression, marital conflict, and family negative emotional climate, can contribute to children’s later mental health symptoms (e.g., (Essex, Klein, Cho, & Kraemer, 2003; Essex, Klein, Miech, & Smider, 2001)). We have shown that these effects sometimes vary as a function of child gender (Essex, et al., 2003; Essex, et al., 2001) and that fathers can make unique, significant contributions to children’s risk and well-being in the early years of development (Boyce, et al., 2005). We have also demonstrated the predictive utility of broader parental stress composites that incorporate the risk factors noted above plus others, including parenting stress and financial stress, for understanding children’s later mental health (Essex, et al., 2002; Essex, et al., 2006), as well as autonomic and adrenocortical reactivity to stress (Ellis, Essex, & Boyce, 2005).
Here, we seek to understand the ”biological embedding” of early experience (Hertzman, 1999). Biological embedding is said to occur when: experience gets under the skin and alters human biological processes; systematic differences in experience, under different nurturant conditions, lead to different bio-developmental states; the differences are stable and long-term; and these differences influence health, well-being, learning, or behavior over the life course. The underlying mechanisms of biological embedding have not yet been fully demonstrated, but interest is now turning to the prospect that epigenetic processes are involved.
It is now well established that information contained in the genome’s inherited nucleotide sequence is complemented by other, more modifiable information inherent in the chromatin structure, or epigenome. At a molecular level, the chromatin material that constitutes human (and all eukaryotic) chromosomes consists of strands of DNA wrapped around octamers of eight histone proteins, resembling beads on a string. Chromatin modifications, which include DNA methylation, histone acetylation, and a variety of other chemical ‘marks’, render the protein-encoding portions of a given gene more or less accessible to the molecular transcriptional apparatus that decodes the DNA sequence into messenger RNA and then specific protein products. The epigenetically mediated, differential expression of such protein products in turn guides or changes the function of the cell and the tissue or organ of which it is a part. Emerging epigenetic science is thus “the study of stable alterations in gene expression by non-genomic mechanisms, resulting in stable alterations in phenotypes” (Kramer, 2005). Although DNA methylation is just one of an array of epigenetic mechanisms, it is the one regarded as most stable and most representative of the broader epigenetic modification of a given locus.(Hochberg, et al.). For comprehensive reviews of this rapidly expanding field of study, readers are referred to: (Champagne & Mashoodh, 2009; Cole, 2009; Curley, Jensen, Mashoodh, & Champagne; Hochberg, et al.; Meaney, 2010).
Our working hypothesis in the current study was that early experience leaves a mark on the epigenome that, in turn, leads to stable changes in expression of genes critical for human development and health (Szyf, McGowan, & Meaney, 2008). It is known that the enzymatic machinery responsible for creating and maintaining DNA methylation is responsive to environmental exposures during both intra-uterine development and after birth in animals (e.g., (Jirtle & Skinner, 2007; Weaver, Cervoni, et al., 2004)). Most notably, epigenetic changes in expression of the gene encoding the glucocorticoid receptor (GR) are mediated, during an early sensitive period, by naturally occurring differences in maternal behavior, specifically the intensity and frequency of licking, grooming and arched-back nursing of rat pups (Meaney, 2010; Weaver, Diorio, Seckl, Szyf, & Meaney, 2004). Parallel evidence is now accumulating in humans (e.g., (Heijmans, et al., 2008; McGowan, et al., 2009; Oberlander, et al., 2008; Terry, et al., 2008)). The GR gene has shown increased methylation in human cord-blood DNA from newborns of depressed or anxious mothers (Oberlander, et al., 2008). In other work, the imprinted IGF2 gene has shown reduced methylation in blood DNA of sixty year-old individuals who had been prenatally exposed to famine (Heijmans, et al., 2008). It is speculated that epigenetic processes, including methylation, transduce early environmental signals into conditionally adaptive changes in metabolic, endocrine and neuroregulatory pathways. These changes, in turn, are responsible for systematic developmental biases toward distinct profiles of growth, metabolism, immune responsiveness, developmental pace, and behavior.
Epigenetic methylation has the capacity to modify gene expression, although it does not always do so. This capacity is nonetheless one piece of evidence making plausible the proposition that DNA methylation is important for human developmental plasticity (Gluckman, Hanson, Buklijas, Low, & Beedle, 2009). In this report, we ask a logical next question: are there lasting DNA methylation differences in the tissues of adolescent children who spent their infant or preschool years in environments of maternal and paternal stress? In other words, do typical differences in daily experience during the early years influence DNA methylation over developmental time?
The current study is based on a sub-sample of 109 children and their parents who participated in the Wisconsin Study of Families and Work (WSFW). The WSFW originally comprised 570 pregnant women and their partners who were recruited from prenatal clinics in and around Madison and Milwaukee, Wisconsin for a study of maternity leave and health outcomes. Eligible female participants were over the age of 18 years, in the second trimester of pregnancy, living with the baby’s biological father, and either employed or a full-time homemaker (see (Hyde, Klein, Essex, & Clark, 1995) for details). Of those eligible, 75% agreed to participate. Because the children were born over a 16-month period, they entered primary school in two cohorts (1996 and 1997). The present analyses are based on data from the 249 families in Cohort 1 who participated at the time of school entry. Parental stress data were obtained during the infancy and preschool periods, and child DNA samples were collected at age 15 years (M = 15.1 years, SD = 0.5). At age 15, 201 (81%) of Cohort 1 participants continued in the study; of these, 110 agreed to participate in the DNA collection and provided a sample sufficient for analysis. Methylation results for one participant were dropped as outliers. A total of 109 participants were thus included in analyses presented here.
At recruitment, the 109 mothers ranged in age from 20–41 years (Mdn = 29). Four percent of mothers had less than a high school education; 13% were high school graduates or equivalent; 31% had some post-secondary education or training; 38% were college graduates; and 14% had post-baccalaureate schooling or a degree; fathers’ education levels were comparable. Nearly all couples were married (95%); 44% were first-time parents; and 11% were ethnic minorities. Family incomes ranged from less than $10,000 per annum to $120,000 (Mdn = $45,000). There were no statistically significant differences between the demographic characteristics of the 109 families and the remainder of the 570 original families. There was a significant difference in the percentage of girls in the 109 sub-study families compared to that in the remainder (61% versus 49%, respectively; χ2 = 4.32, df = 1, p = .04). The 109 and remaining families did not significantly differ on measures of parental stress.
Measures of parental stress were obtained separately from mothers and fathers at multiple times during their children’s infancy and preschool periods. For each parent, stress scores were composites of five domains: 1) depression symptoms, 2) family expressed anger, 3) parenting stress, 4) role overload, and 5) financial stress. Except where noted, scores for the infancy period were averaged over the three assessments conducted at 1, 4, and 12 months post-delivery; scores for the preschool period were averaged over the two assessments conducted at children’s ages 3.5 and 4.5 years. Parental depression symptoms were assessed by the Center for Epidemiologic Studies-Depression scale (CES-D; (Radloff, 1977)). Across all assessments, αs were ≥ .85 for both mothers and fathers. Family expressed anger in infancy was assessed with the average of three items from the Partner Role Quality scale tapping overt marital conflict (e.g., concerned about “arguing or fighting”; (Barnett & Marshall, 1989)); for the preschool assessments, the marital conflict scores were combined, using principal components analysis (PCA), with scores from the Anger Expression Inventory (Spielberger, 1988) and the Negative subscale of the Family Expressiveness Questionnaire (Halberstadt, 1986). Across all assessments, for both mothers and fathers, αs were ≥ .78 for marital conflict, .68 for anger expression, and .85 for negative family expressiveness. The first component of the PCA for each preschool assessment accounted for more than 50% of the variance. Parenting stress was calculated (PCA weights) from scores on the Competence and Child Reinforces Parent subscales of the Parenting Stress Inventory (PSI; (Abidin, 1986)) and the average of three negative items (e.g., “I often feel angry with my child”) from the Childrearing Practices Report (CRPR; (Block, 1965)); in the infancy period, the CRPR was assessed only at 12 months. All αs were ≥ .76 for PSI Competence, .64 for PSI Child Reinforces Parent, and .62 (except at age 12 months, where α = .49) for parent-child negativity. The first component of each PCA accounted for more than 50% of the variance. Parental role overload was calculated (z scores, averaged; r range = .43 – .50, ps < .001) from the Role Restriction subscale of the PSI (Abidin, 1986) and the average of five items from the Role Overload scale (e.g., feeling “pulled apart by conflicting obligations”; (Barnett & Marshall, 1989)); in the infancy period, the Role Overload scale was assessed only at 12 months. For the two scales, all αs were ≥ .77. Financial stress was the average of four items (e.g., “how much difficulty making monthly payments”); all αs were ≥ .74. Finally, separately for the infancy and preschool periods and for mothers and fathers, scores for the five stress domains were combined using PCA. The first component of each of the four PCAs accounted for more than 50% of the variance, and all stress domains had factor loadings greater than .50 for each period and for each parent. Four parental stress variables were thus employed in subsequent analyses: mother and father stress scores for both the infancy and preschool periods. Due to parental non-participation at some time points, the sample size was somewhat different for the analysis of each specific parental stress variable: thus, for infancy, n = 106 mothers and 100 fathers, and for preschool, n = 106 mothers and 85 fathers. For each parent, the infancy and preschool stress scores were strongly correlated (r = .74 for mothers, .73 for fathers, ps < .001), whereas there were only very modest correlations between the mother- and father-report scores within and across the two developmental periods (r range = .26 – .34, ps < .001). To deal with skewness, log-transformed variables were used in analyses.
Before collection, participants rinsed their mouths with water twice for 15 s each. Buccal epithelial cells were collected using MasterAmpBuccal Swabs (Epicentre Biotechnologies), which were stored at −80°C prior to analysis. Genomic DNA was isolated using IsohelixBuccal DNA Isolation Kits (Cell Projects) and was purified and concentrated using DNA Clean & Concentrator (Zymo Research). Sample yield and purity were assessed spectrophotometrically using a NanoDrop ND-1000 (Thermo Scientific).
Quantitative DNA methylation measurements of purified genomic DNA were performed with the Infinium Human Methylation 27 BeadChip assay (Illumina). This microarray platform allows for the simultaneous measurement of the DNA methylation status within 27,578 cytosine-guanine dinucleotide (CpG) sites, covering more than 13,500 promoters, first exons of the coding region, or occasionally first intron of well-annotated genes, thereby representing more than half of all human genes. Each array also contains 254 CpG sites located in approximately 110 micro RNA loci. After random ordering of samples, 750ng of genomic DNA from each sample was bisulfite-converted using the EZ DNA Methylation Kit (Zymo Research). Bisulfite deaminates unmethylated cytosine to uracil, while methylated cytosine is protected from deamination. Bisulfite treatment thus converts epigenetic information to sequence-based information, which can be measured with methods similar to those used to distinguish single nucleotide polymorphisms. A sample of 160ng of bisulfite-converted DNA was whole-genome amplified, fragmented by an enzymatic process and hybridized to BeadChip arrays. Two oligonucleotide probes interrogated each CpG site, one probe with sequences targeting methylated DNA and the other containing sequences targeting unmethylated DNA. After extension with DNP-labeled and biotin-labeled dNTP, each array was stained with Cy5 labeled anti-DNP antibodies and Cy3 labeled streptavidin and scanned with the IlluminaiScanon a two-color channel to detect Cy3 labeled probes on the green channel and Cy5 labeled probes on the red channel. Using the IlluminaGenomeStudio software package, average beta values were calculated by dividing the methylated probe signal intensity by the sum of methylated and unmethylated probe signal intensities. Average beta values range from 0 (completely unmethylated) to 1 (fully methylated) and provide a quantitative readout of relative DNA methylation for each CpG site within the cell population interrogated. This method was highly reproducible, with replicates across runs achieving r > .99.
All normalization, calculations and statistical analysis was done using R 2.11.0 (http://www.R-project.org/). Average background intensity, measured by negative background probes present on the array, was subtracted from raw intensities to adjust for varying background signals across samples. Background adjustment was completed separately for raw data from the green and red channels to adjust for Cy3 and Cy5 differences. To control for batch effects across different sets of arrays, background adjusted raw data from both channels were quantile normalized separately (Bolstad, Irizarry, Astrand, & Speed, 2003). Average beta values were then recalculated per the following Illumina formula: background subtracted and quantile normalized intensities of methylated probes divided by the sum of normalized intensities from unmethylated and methylated probes. All CpG probes with a detection p > .05 were removed to ensure that only high-confidence probes were included in the subsequent analysis. The detection p-value assesses whether the signal intensity for the average beta value is significantly different from the negative probes, as well as consistent across all oligonucleotides interrogating a given CpG site on the array. CpG sites on the X-chromosome were removed to prevent confounding by gender. To increase analytic power, CpG sites whose average beta values across subjects were < 0.05 or > 0.95 were excluded, so that all sites would have between-sample differences in methylation > 5% of the average beta value.
This left a total of 18,231 CpG sites to be tested for associations with the four parental stress measures. Significance Analysis of Microarray (SAM, in the samr package of Bioconductor (Tibshirani, Chu, Hastie, & Narasimhan, 2010; Tusher, Tibshirani, & Chu, 2001)), using a ranked linear regression analyses, identified CpG sites with DNA methylation patterns associated with any of the stress variables. SAM assigns a score to each CpG site on the basis of differential DNA methylation relative to the standard deviation of repeated measurements. To account for multiple testing, the procedure uses permutations of repeated measures to estimate the false discovery rate (FDR), defined as the expected proportion of false positives among a ‘family’ of hypothesis tests. This allows the identification of a set of ‘candidate positives’ and is useful in statistical settings where priority is placed on the discernment of mostly true findings, rather than the strict avoidance of single false positives (Storey, 2011). The computed FDR thus measures the likelihood of a Type I error within a particular association, in a manner parallel to the p-value, but the background distribution to estimate the Type I error is generated from random assignment of stress variables to DNA methylation data. The q-value is the minimum FDR at which the test may be designated significant. In all cases, 1000 permutations were run and included FDRs up to 20%. Separate SAM analyses were run for the entire group of 109 subjects where stress data were available and for gender-specific subsamples. Reported results show differences in DNA methylation between the low and the high end of the regression line defined by regressing DNA methylation level on a parental stress measure.
Finally, the functional annotation tool from the DAVID bioinformatics resource (Huang da, Sherman, & Lempicki, 2009) was used to identify common biological processes and pathways, molecular functions, and cellular components, as defined by Gene Ontology (GO) criteria. For this analysis, and using the default, “medium stringency” setting of the DAVID instrument, the enrichment of GO terms in the target gene list was compared, using Fisher’s Exact test, to a background list containing all genes analyzed. The Functional Annotation Cluster tool aggregates enriched GO terms based on their co-occurrence in the same gene; the tool also computes an enrichment score, the minus log transformed geometric mean of p-values from GO terms within each cluster. We report the single cluster from the DAVID analysis of 139 methylation sites associated with maternal stress at infancy that had an enrichment score > 1.3, corresponding to a p of 0.05.
High levels of early adversity were defined as the upper 20% of PCA-derived scores for parent reported stress, for each reporter and time period. Among mothers reporting high stress during infancy, for example, at baby’s age 12 months the average CES-D score was 15 (with 43% ≥ the clinical threshold of 16 and 76% ≥ the clinical sub-threshold of 12), 67% reported marital conflict, 95% felt overloaded “occasionally” to “very often,” and 81% felt financially stressed. Among fathers reporting high stress in preschool, at child’s age 4½ years, the average CES-D score was 20 (with 59% ≥ 16 and 76% ≥ 12), 71% reported marital conflict, 88% felt overloaded “occasionally” to “very often,” and 65% felt financially stressed.
Raw data from the scanned microarrays are available to readers in the gene expression omnibus (GEO) database, a repository of high throughput gene expression data and hybridization arrays, which can be accessed under accession number GSE25892 at: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25892.
Differences in DNA methylation were significantly associated with both maternal and paternal stressors reported 10–15 years prior to the collection of the children’s DNA. Table 1 is a summary of these results for both the full sample and the child gender-specfic subsamples. Because of the large number of associations identified, the chance of reporting false positive associations is reduced by focusing only on CpG sites with FDRs < 20%. Sites are further differentiated as high confidence sites (HCSs) when FDRs were < 5% and medium confidence sites (MCSs) when FDRs were between 5% and 20%. Finally, the number of sites with > 5% difference in methylation between low and high stress conditions is indicated in parentheses in the Table, and the range of site-specific Spearman rho correlation coefficients is displayed for associations between stress and site methylation. Figure 1 is a visual representation of the same data for HCSs only. Readers interested in a more detailed rendition of these results are referred to a web site, http://www.cmmt.ubc.ca/research/investigators/kobor/lab/childstress, where supplementary tables list methylation status for all specific CpG loci associated with parental stress. The posted tables display, for each CpG site, the specific gene symbol, the site interrogated, the slope of the regression line relating parental stress to DNA methylation level, the range of relative DNA methylation, the percent differential DNA methylation between low and high stress children, Spearman rho correlation coefficients, and the FDR. Shaded rows identify sites with DNA methylation differences exceeding 5%.
In the full sample, associations of DNA methylation were found with two of the four stress variables. First, maternal stress in infancy bore high confidence associations, with increased DNA methylation in 139 CpG sites, all of which had high confidence FDRs of 0% and 37 of which had more than a 5% difference in DNA methylation between children exposed to low and high maternal stress in their first year of life (Suppl Table 1). Further, rank order correlations between maternal stress in infancy and increased methylation of CpG sites in adolescence ranged from .21 to .40. Second, although the association of maternal stress in infancy with DNA methylation did not extend into the preschool years, paternal stress during this period was associated with increased DNA methylation in 31 CpG sites, of which 2 were HCSs and 12 had methylation differences > 5% (Suppl Table 2). Rank order correlation magnitudes for paternal stress in the preschool period ranged from .34 to .45.
The microarray contained a limited number of CpG sites located in the gene regulatory regions for several candidate genes in which DNA sequence polymorphisms or differential DNA methylation have known associations with environmental stress or behavior. These included genes for the glucocorticoid receptor (NR3C1), dopamine receptor (DRD4), serotonin transporter (5 HTT), brain-derived neurotrophic factor (BDNF), catechol-O-methyltransferase (COMT), and dopamine transporter (DAT1). No significant associations were identified between DNA methylation in these genes and the measures of parental stressors; Spearman rho correlations ranged from −.39 to .26 but failed multiple testing correction (FDR > 20% in cases; detailed data are shown in Suppl Table 3).
Repeat analyses for male and female children produced patterns of association distinctive from but commensurate with those observed in the full sample. In the infancy period, for example, increased exposures to maternal stress were associated with differential DNA methylation in both genders, consistent with findings from the full sample. In contrast to the large number of sites identified in the full sample, however, far fewer emerged in gender-specific analyses. The fewer differentially methylated sites are likely due, at least in part, to reductions in statistical power among the gender sub-samples, but all identified associations had substantially stronger correlations. The 3 HCSs with increased DNA methylation identified in boys had rank order correlation coefficients in the range of .50 – .58 (Suppl Table 4), and the single HCS found in girls had rho = .50 (Suppl Table 5). Further, differential methylation among CpG sites in girls was related to both maternal and paternal stressors in infancy, but in opposing directions. That is, maternal stress in infancy was associated with increased methylation of one CpG site (Suppl Table 5), whereas paternal stress was linked to decreased methylation at three CpG sites, with the latter having rank order correlation coefficients in the range of −.45 – −.52 (Suppl Table 6). In boys, on the other hand, only maternal stress in infancy was associated with differential DNA methylation, in the direction of increased methylation among stress-exposed boys.
Gender-specific analyses of parental stress in the preschool period indicated that differential methylation was uniquely associated with exposure to paternal stress and was specific to girls. In these findings, analyses identified distinctive sets of CpG loci that were either more methylated or less methylated with increased exposure to paternal stress. Specifically, 320 CpG sites were found with increased paternal stress-related DNA methylation: 6 belonging to the high confidence group (rhos= .52 – .56) and 314 to the medium confidence group (rhos = .34 – .51; Suppl Table 7). These 320 CpG sites included 21 out of the 31 sites identified in the full sample, with rank order correlations becoming stronger in each case (Suppl Table 9). Gender subset analyses also identified 1066 CpG sites that were less methylated with higher paternal stress: 9 classified as HCSs (rhos = −.51 – −.61) and 1057 as MCSs (rhos = −.28 – −.50). Paralleling results for the full sample, there was no differential DNA methylation associated with exposure to maternal stress during the preschool period for either gender.
Taken together, results summarized in Table 1 and Figure 1 indicate that differences in DNA methylation were associated with three contrasting predictors of epigenetic modification: a) the gender of the parent reporting the stressors (mother v father), b) the developmental period in which the stressors were reported (infancy v preschool), and c) the gender of the child (male v female).
To further evaluate these associations, an additional level of stringency was next applied to the analyses. Specifically, we focused on CpG sites within sets of HCSs that had more than a 5% difference in absolute methylation from the lowest to the highest point on the regression ‘fit’ line (calculated as the slope × range of the stress variable). Although the 5% cut-off is arbitrary, it was reasoned that this added stringency would highlight sites with the most robust differentials in DNA methylation and therefore those most likely to have a possible functional impact. The numbers of robust (> 5% methylation difference) CpG loci are indicated in parentheses for each set of associations in Table 1, along with their correlation range. These ranges of rank order correlation coefficients were quite similar to those of the entire set, suggesting that applying the more stringent filter did not bias toward higher correlations but merely removed sites with less pronounced DNA methylation changes.
As noted above, of the 139 HCS CpG loci identified in analyses for exposure to maternal stress during the infancy period in the full sample, 37 had DNA methylation changes exceeding 5% (Suppl Table 1). We noted that, whereas almost all of the 37 CpG loci were from different genes, two belonged to the PGAM2 gene, which encodes phosphoglyceratemutase 2. These two associations are presented as scatterplots in Figure 2. As before, analysis by gender reduced the number of CpG sites associated with maternal stress in infancy. There were only 2 sites in the gender-specific analysis exceeding 5% methylation difference, but both were also identified for the full sample. Specifically, these sites were: a) a CpG site in the PKN1 gene promoter for boys, and b) a site in the C5orf21 promoter for girls. Scatterplots for these associations are presented in Figures 3A and 3B. As shown, the rank order correlations became stronger and their differences in DNA methylation more pronounced when analyzed by gender.
Filtering for greater robustness removed 2 of the 3 HCS CpG loci associated with paternal stress during infancy in girls. The remaining CpG site in the promoter of MGMT gene, which encodes O-6-methylguanine-DNA methyltransferase, had more than a 20% difference in methylation, the most pronounced difference in the study (Suppl Figure 1). Results for this site should be considered with caution, however, as our array contained 25 additional CpG sites in the MGMT gene, none of which were associated with parental stress.
Application of the filtering procedure to the analysis of paternal stress during the preschool period in the full sample showed that the 2 HCSs had less than 5% DNA methylation difference. For 12 of the 29 MCSs, however, the differential DNA methylation exceeded 5%, ranging up to 14% (Suppl Table 2). Of the 15 HCSs identified in gender-specific analyses as being either more or less methylated following paternal stress in girls, 8 had DNA methylation changes exceeding 5% (Suppl Tables 7 and 9). This relation between differential DNA methylation and preschool paternal stress in girls is illustrated in Figure 4. A CpG site in the TRDN gene promoter showed increased methylation in girls as paternal stress increased (Fig 4A), whereas a CpG site in the MKRN1 first coding exon was less methylated under the same conditions (Fig. 4B). Notably, the C19orf30 first coding exon contained 2 CpG sites that had more DNA methylation with increasing paternal stress in girls (Suppl Figure 2). One additional CpG in the C19orf30 gene belonged to the set of MCSs that were less methylated as paternal stress increased, while 5 others in the same gene were not associated with parental stress according to our criteria.
The above analyses focused principally on HCS sites, which had less than 5% FDRs, as a means of minimizing the likelihood of reporting false positive associations. Nonetheless, we also sought for ways to test both HCS and MCS loci for functional biological information. To do this, we searched the combined HCS and MCS datasets for increased and decreased DNA methylation of paternal stress at preschool for genes with more than one CpG site having DNA methylation change going in the same direction. This was motivated by our finding that some genes (e.g., PGAM2, c19orf30) had 2 associated CpG sites and by the underlying assumption that identification of several associated CpG sites in the same gene might reflect greater likelihood of a significant association. A large number of genes were identified using this approach (n=71 genes), prompting a reapplication of the 5% DNA methylation change screen, as described above (Suppl Table 10). This approach reduced the number of genes to 28, including several genes of potential functional interest. IGF2AS (3 out of 11 CpGs on the array with FDR < 20%; 1 out of 3 with > 5% change in DNA methylation), for example, is a gene that expresses a paternally imprinted antisense transcript of the insulin-like growth factor 2 gene. Another example, the TGFB2 gene (2 out of 2 CpGs on the array with FDR < 20%; 1 out of 2 with > 5% change in methylation) encodes a member of the transforming growth factor beta (TGFB) family of cytokines. A final example, the NNAT neuronatin gene (2 out of 7 CpGs on the array with FDR < 20%; both with > 5% change in methylation), is an imprinted gene encoding a proteolipid possibly involved in the regulation of ion channels during brain development.
Further, 2 genes in the reduced set of 28 had 3 CpG sites with > 5% change. These were the GNAS gene cluster (6 out of 29 CpGs with FDR < 20%; 3 out of the 6 with > 5% difference in DNA methylation, the latter shown in Figure 5) and the MGC33302 gene, also known as MFSD8/CLN7 (3 of 6 CpGs with FDR < 20%; all with > 5% change in methylation and all shown in Supp Figure 3). Although it is not clear whether the DNA methylation pattern measured in buccal epithelial cells is conserved in more central tissues, both of these loci are of potential interest for brain development and behavior. GNAS is an imprinted locus consisting of a complex arrangement of several protein-coding and non-coding transcripts. In mice, deletion of Nesp55, one of the protein-coding genes located in the GNAS locus and under the control of its imprinting region, influences behavioral reactivity to novel environments (Plagge, et al., 2005). Genetic variation in MFSD8/CLN7 is associated with the variant late-infantile onset neuronal ceroid lipofuscinoses, the most common neurodegenerative disorders of childhood (Siintola, et al., 2007).
In addition to gene-specific analyses, it is often informative to look at broader patterns of genes with DNA methylation differences to identify cellular or biological process commonalities. To do this, we focused on the set of stress-associated CpG sites containing only HCSs, as this minimized the chances of reporting commonalities of lesser significance. The web-based DAVID platform (see: http://david.abcc.ncifcrf.gov/) is a project of the National Institute of Allergy and Infectious Diseases that was designed and published to foster functional annotation analyses in microarray and proteomics studies (Huang da, et al., 2009). Because the algorithm is optimized for gene lists containing 100–1000 genes, the 139 HCSs that were associated in the present study with maternal stress in infancy were selected as the set of genes on which to perform the DAVID analysis and illustrate its potential utility. The functional annotation cluster that was identified using this approach had an enrichment score of 1.45, which is above the 1.3 value that corresponds to p<0.05 and represents the minimal score considered significant (Suppl Table 11). This cluster contained GO terms centered on biosynthetic and metabolic processes.
The study reported here had four principal findings of relevance to genetic research on child development. First, maternal stress during the infancy of a child—long suspected and more recently documented as a contributor to maladaptive developmental outcomes in children—was associated with a broad and substantial array of epigenetic marks assessed a decade and a half later, during children’s mid-adolescence. Interrogating nearly 28,000 CpG sites in the regulatory regions of 14,000 genes, a large number of DNA methylation sites were identified for which children whose mothers who were highly stressed in the child’s first year of life showed significantly higher levels of methylation than those less exposed, even after controlling for false discovery rates related to multiple comparisons. Such genomic sites with increased DNA methylation occurred both in genes with plausible linkages to family experiences of stress and adversity (e.g., NEUROG1, a bHLH transcription factor gene that is involved in neuronal differentiation and cell-type specification in distinct regions of the developing nervous system) and in genes with no such linkages (e.g., PGAM2.). To our knowledge, this is the first report of a long-term, prospective association between early childhood parental adversities and the epigenetic status of specific, regulatory sites in children’s genomic DNA.
Second, the presented data showed that experiences of adversity among both mothers and fathers bore reliable linkages to epigenetic profiles within specific developmental timeframes, with maternal stressors during the infancy period more potently predictive and paternal stressors during the preschool years more predictive. Fathers’ reports of high levels of adversity during their children’s preschool years were significantly associated with altered methylation status in a smaller, but still significant, number of regulatory loci. Findings therefore offer evidence for epigenomic correlates of biparental adversities, with sequenced, period-specific sensitivity to the stressors of one parent over the other.
Third, the patterning of epigenetic marks associated with early adversity varied by the gender of both parent and child, with fathers’ reported stressors more strongly associated with girls’ epigenetic modifications and mothers’ more related to those of both sexes. Such findings corroborate, at an epigenetic level of analysis, past observations of gender-specific influences of parenting on child development and well being. Father absence and non-participation in parenting has been linked, for example, to an earlier timing of puberty onset (Ellis, McFadyen-Ketchum, Dodge, Pettit, & Bates, 1999) and to the emergence of difficult temperamental characteristics in daughters, rather than in sons (Bezirganian & Cohen, 1992). The current study provides provisional evidence that sex-specific influences of parent-child relationships on maladaptive developmental processes may be subserved, in part, by epigenetic modifications capable of guiding such processes.
Fourth, the noted associations were identified epigenetically within buccal epithelial cells, a histological cell type not previously examined as a target tissue for social contextual, exposure-related epigenetic changes. Such DNA methylation differences in epithelial chromatin are likely not paralleled by similar marks in neural tissues, which are generally assumed to be the biological substrate for developmental events and processes. Nonetheless, the adversity-related epigenetic marks reported here constitute measurable epigenomic change that is associated with early life experience within proximal social contexts and has the potential for altering gene expression. Further, the epigenetic profiling of even readily accessible epithelial cells can advance understanding of issues such as the timing and persistence of early exposures’ biological correlates and the consistency of profiles across populations, geographies and time.
Taken together, these findings offer novel evidence for a biological embedding of early experience, or more specifically, the temporally remote correlates of early adverse experiences on the human epigenome and its regulatory role in the expression of specific genes, including genes that guide neurodevelopment. Both relative increases and decreases in promoter region methylation were detected in the genomes of adolescents whose parents had reported significant adversity in past years. Although DNA methylation is generally associated with a down-regulation in gene expression, the epigenetic control of differential transcription is almost certainly far more complex and involves a broad and diverse array of chromatin modifications (see, e.g., (Mehler, 2008)). As a consequence, discerning a coherent and functional ‘meaning’ of the reported findings lies beyond the present state of epigenetic science. We note, however, that to our knowledge these data constitute only the second report of altered DNA methylation in buccal epithelial cells associated with environmental exposure, following a recent study that reported DNA methylation marks in such cells from 5-year-old children associated with maternal smoking during the intra-uterine period (Breton, et al., 2009).
Importantly, our findings may be correctly viewed as an instantiation of ‘gene-environment interplay’ and the capacity for experience and genomic variation—allelic or epigenetic—to conjointly influence salient developmental endpoints (Gilbert & Epel, 2009; Rutter, Moffitt, & Caspi, 2006). One variety of such interplay is the moderation of early experiential effects by single nucleotide polymorphisms within genes affecting key neural circuitry and neurotransmission pathways (Caspi, et al., 2002). The most broadly recognized examples of epigenetic regulatory processes are those described within the caregiving behavior of the mother rat and within the dietary influences on coat color in mice. Studies by Meaney, Szyf and colleagues (Meaney & Szyf, 2005; Weaver, Cervoni, et al., 2004) have shown how phenotypic differences in the reactivity of the HPA system arise from the mother rat’s licking and grooming behavior in her pups’ first several post-natal days. Such behavior changes reactivity phenotypes by demethylating the binding site for transcription factor Egr1 in the enhancer region of the pups’ glucocorticoid receptor gene. Waterland and Jirtle (Waterland & Jirtle, 2003) supplemented the diets of pregnant, viable yellow Agouti mice with methyl-donors, such as folate, choline and betaine, and created dramatic differences in offspring coat color. Pups of mothers fed methyl-donor supplements had increased methylation of the Agouti allele that guides yellow fur development, resulting in suppressed gene expression and phenotypic reversion to a brown coat color. Far more limited examples of epigenetic processes in humans have included research by McGowan and colleagues (McGowan, et al., 2009; McGowan, et al., 2008) in the Meaney laboratory, revealing DNA methylation changes in promoter regions of the ribosomal RNA and glucocorticoid receptor genes among individuals with a history of childhood abuse and suicide. A recent report by van IJzendoorn and colleagues found higher levels of promoter methylation in the serotonin transporter gene among adult adoptees with unresolved trauma (van IJzendoorn, Caspers, Bakermans-Kranenburg, Beach, & Philibert, 2010). No prior work, however, has provided prospective human evidence of sustained epigenetic marks related to early life experience. The novelty of the current findings lies in the long-term association of family adversity over the first five years of life with epigenetic changes in a broad variety of methylation sites in adolescence and the identification of such changes in non-invasively accessible epithelial cells.
These findings, however novel, are also usefully commensurate with several past sets of observations in the developmental and biological sciences. A now expansive literature documents, with increasing methodological strengths, the epidemiologic associations between experiences of early adversity in the life a young family and their long term consequences for health and development (Hertzman & Boyce, 2010; Shonkoff, Boyce, & McEwen, 2009). Such literature includes evidence that adverse childhood events are associated with a variety of mental and physical health outcomes in mid-life (e.g., (Felitti, et al., 1998)), that adult chronic diseases are associated with nutritional and growth deficits in very early, even prenatal development (e.g., (Barker, 1990)), that maternal depression during infancy is associated with systematic differences in the regulation of children’s stress responsive biological systems (e.g., (Essex, et al., 2002)), and that growing up in conditions of poverty is linked to proinflammatory shifts in cytokine expression, which is related in turn to risk for chronic cardiovascular disease in adult life (e.g., (Miller, Maletic, & Raison, 2009)). A recent report relating early life stress to differences in both whole genome and candidate gene methylation in bonnet macaques extends such findings to the epigenome (Kinnally, et al., 2011), and the current study adds new evidence that chromatin structure may bear the ‘mark’ of early life stressors, even in human children. Epigenetic variation may thus constitute a biological ‘memory’ of early life experience, perhaps especially experiences of adversity, stress and misfortune.
The study findings also substantiate prior observations that, although maternal stressors are preeminent in their potential influence on the developmental outcomes of infants and young children, such influence involves both parents, with fathers’ effects on adaptation and risk becoming especially salient beyond infancy. Past research has shown, for example, that fathers’ supportive presence has key influences on biological processes, from pubertal timing (Ellis, 2004) to the regulation of stress reactivity (Roubinov & Luecken, 2010), and on developmental processes, from school readiness (Martin, Ryan, & Brooks-Gunn, 2010) to the acquisition of significant, presyndromal behavior problems (Boyce, et al., 2005; Lamb, 1997). Again, at an especially fundamental level of biological impact, current findings add to the body of evidence implicating biparental experiences of adversity and depression in the advent of maladaptive development and disorders of mental and physical health.
Increasingly, studies of human development also identify important gender differences in both trajectories of normative development and in the determinants of deviation from such trajectories (Essex, et al., 2003; Landman-Peeters, et al., 2008). These gender differences in developmental course may be traceable, in part, to the roles of mothers and fathers in the care of and attention to male and female children, a differential involvement that may also differ by socioeconomic status (Roopnarine, Fouts, Lamb, & Lewis-Elligan, 2005). Similarly, the findings of the current study suggest significant differences between boys and girls in the parental-adversity-associated differential methylation of gene promoter regions. Such differences occurred both in the magnitude of changes in DNA site methylation, as well as in the specific epigenetic sites involved.
Among the biological and methodological issues involved in studying chromatin structure are the profound epigenetic differences between cells of differing histological origin. Indeed, since all cells in the body contain the identical DNA sequence, it is the epigenetic differentiation of primordial, germ layer cells that allows and guides the embryogenesis of different tissues and organs. Thus, although the neural tissue epigenome is closest to the neurobiological processes that drive development, human epigenetic research is necessarily relegated to the examination of peripheral, accessible tissues, such as immune cells, gametes, and biopsy materials. A recent report by Talens et al, however, offers methodological assurances of meaningful variation in the human BEC epigenome, stability in such variation over years of time, and substantial overlap between BEC and leukocyte epigenetic marks (Talens, et al., 2010). The present study offers additional evidence that experience-related changes in the epigenome are detectable in the DNA of epithelial cells, which are readily accessible to sampling and study.
Interpretation of the results presented here should be informed by an awareness of specific weaknesses in the design and methods of the reported study. First, the current study sample comprised a subset of children from a larger, longitudinal study of child development and mental health. Although the examined subsample did not differ substantially in a variety of demographic descriptors, there remains the possibility that the children studied differed systematically from the larger cohort from which they were drawn. Second, while representative of the specific population from which participants were identified, findings may not be replicable within samples derived from more racially or socioeconomically diverse groups. Third, despite extensive efforts to minimize the likelihood of Type I errors, the search for differences among 28,000 potential methylation sites involved issues of multiple comparison, to say the least, and it is possible that false positive findings escaped efforts to preclude such errors. Fourth, the current study did not consider the epigenetic changes associated with prenatal versus infancy and preschool stress exposures. Replication and extension of the reported findings will thus be essential. Finally, useful extensions of the reported analyses would examine behavioral and psychological phenotypes that may be associated with both early parental stressors and differential methylation of gene promoters. Such analyses are now within reach of several research teams, and the co-authors of the current report look forward to pursuing linkages among early exposures, epigenetic marks and phenotypic variation in these and other, related data.
Even taking into account these limitations in design and methods, however, the reported research yielded singular observations linking parental adversity in early childhood to differences in chromatin structure a decade and a half later, in the mid-adolescence of the children studied. Such epigenetic differences may constitute a fundamental biological vestige of such early exposures, potentially altering the expression of genes affecting metabolic and physiologic pathways and changing trajectories of individual phenotypic development. As such, the observed epigenetic modifications could plausibly represent a ‘fundamental cause’ of the risk-enhancing trajectories leading to disease, disorder and troubled developmental outcomes.
This work was funded by the National Institute of Mental Health grants R01-MH044340, P50-MH069315, P50-MH084051, and R24-MH081797 and by research support to WTB, CH and MSK from the Canadian Institute for Advanced Research. MSK is a Scholar of the Mowafaghian Foundation through the Human Early Learning Partnership of the University of British Columbia. WTB holds the Sunny Hill Health Centre-BC Leadership Chair in Child Development. The authors wish to express their appreciation to the participating families who generously committed their time to the project and to Dr. Hunter Fraser (Stanford University), who provided important statistical assistance in the interpretation of data.
Marilyn J. Essex, University of Wisconsin School of Medicine and Public Health.
W. Thomas Boyce, University of British Columbia.
Clyde Hertzman, University of British Columbia.
Lucia L. Lam, University of British Columbia.
Jeffrey M. Armstrong, University of Wisconsin School of Medicine and Public Health.
Sarah M.A. Neumann, University of British Columbia.
Michael S. Kobor, University of British Columbia.