- Open Access
Longitudinal associations of DNA methylation and sleep in children: a meta-analysis
Clinical Epigenetics volume 14, Article number: 83 (2022)
Sleep is important for healthy functioning in children. Numerous genetic and environmental factors, from conception onwards, may influence this phenotype. Epigenetic mechanisms such as DNA methylation have been proposed to underlie variation in sleep or may be an early-life marker of sleep disturbances. We examined if DNA methylation at birth or in school age is associated with parent-reported and actigraphy-estimated sleep outcomes in children.
We meta-analysed epigenome-wide association study results. DNA methylation was measured from cord blood at birth in 11 cohorts and from peripheral blood in children (4–13 years) in 8 cohorts. Outcomes included parent-reported sleep duration, sleep initiation and fragmentation problems, and actigraphy-estimated sleep duration, sleep onset latency and wake-after-sleep-onset duration.
We found no associations between DNA methylation at birth and parent-reported sleep duration (n = 3658), initiation problems (n = 2504), or fragmentation (n = 1681) (p values above cut-off 4.0 × 10–8). Lower methylation at cg24815001 and cg02753354 at birth was associated with longer actigraphy-estimated sleep duration (p = 3.31 × 10–8, n = 577) and sleep onset latency (p = 8.8 × 10–9, n = 580), respectively. DNA methylation in childhood was not cross-sectionally associated with any sleep outcomes (n = 716–2539).
DNA methylation, at birth or in childhood, was not associated with parent-reported sleep. Associations observed with objectively measured sleep outcomes could be studied further if additional data sets become available.
Sleep is a key aspect of healthy functioning; disruptions in sleep as early as during childhood have been linked to a wide range of mental and physical health problems . Sleep is a complex, multifactorial phenotype, reflecting the influence of both genetic and environmental factors, beginning in utero. For example, prenatal exposure to maternal smoking , alcohol use  and depression [4, 5] during pregnancy have been linked to disrupted sleep in children. However, the biological factors underlying sleep development in childhood remain unclear.
Epigenetic processes have received much attention as a biological mechanism that may underlie differences in child development and mediate the effects of foetal or early-life exposures on later health. The most widely studied of epigenetic processes, DNA methylation (DNAm) of cytosine nucleotides followed by guanine (CpGs), has been linked to a variety of mental and physical health outcomes [6, 7], and to foetal exposures such as maternal smoking . It has been proposed that DNAm may also serve as a mechanism by which genetic and environmental exposures influence sleep . DNAm is also showing promise as a biological marker for disease prediction, early detection and risk stratification. This application is especially well suited for use in peripheral tissues, which are more readily available in humans but may not be causal for the phenotype of interest. For example, it is already possible to estimate a range of exposures, traits and health outcomes based on peripheral DNAm patterns alone (e.g. age, smoking, BMI) .
While associations between epigenetic changes and sleep have been previously assessed, these have mainly been investigated within the context of experimental studies of sleep deprivation or sleep disruption experienced during shift work [11,12,13]. Studies on DNAm in relation to sleep characteristics among children or adolescents are scarce. One epigenome-wide study utilizing a network-based approach identified a module of inter-correlated CpG sites in peripheral blood associated with actigraphy-estimated sleep duration among 10-year-olds (n = 188) . Another epigenome-wide pilot study reported an association between blood DNAm and diary-based sleep duration among 18- to 19-year-olds (n = 26) , and a third epigenome-wide study reported differences in buccal cell DNAm between 18- and 19-year-old monozygotic twins who were discordant for diurnal preference (n = 30) . One candidate region study reported an association between DNAm of long interspersed nuclear elements and metabolism-related genes in leukocytes, and actigraphy-estimated sleep duration and fragmentation among 14-year-olds (n = 351) . These prior studies, while encouraging as proof-of-concept, are limited due to small sample size and cross-sectional design, and vulnerable to chance findings due to single-cohort design.
Since in utero exposures and genetic factors are known to influence sleep outcomes later in life, we posit that DNAm at birth may be an early-life marker of sleep disturbances . Our primary aim was to investigate if DNAm (at birth) is prospectively associated with parent-reported sleep duration in children. As secondary aims, we investigated two other parent-reported sleep-related phenotypes (sleep fragmentation or initiation problems) and three actigraphy-estimated sleep measures (duration, sleep onset latency and wake-after-sleep-onset duration). In further analyses, we also investigated associations between childhood (4–13 years) DNAm and sleep characteristics. For the current meta-analysis, a total of 14 cohorts in the Pregnancy And Childhood Epigenetics (PACE) consortium conducted one or multiple epigenome-wide association study (EWAS) analyses to address the primary and secondary aims of our study (n = 3658 participants in primary analysis).
Fourteen North American and European cohorts, all members of the Pregnancy and Childhood Epigenetics Consortium (PACE) , participated in this meta-analysis. Eleven cohorts had data on newborn cord blood DNAm at birth and child sleep (Table 1): of these, ten cohorts had data on parent-reported sleep outcomes and three had data on actigraphy-estimated sleep outcomes. Eight cohorts had data on DNAm in peripheral blood in childhood and child sleep (Additional file 1: Table S1). Overlap within each cohort was large across data sets, i.e. same individuals formed the majority of the analytical sample in different cohort-level analyses of, for example, DNAm at birth and DNAm in childhood (within the same cohort): in total, five cohorts had data on DNAm both in cord blood and in peripheral blood in childhood (Additional file 1: Table S2). Individual cohorts are described in more detail in Additional file 2: Methods.
Exposure: DNA methylation at birth and in childhood
Cohorts collected samples of newborn blood from umbilical cord blood at birth (Table 1) and child blood using venepuncture at 4–13 years of age (Additional file 1: Table S1). DNAm was assessed with the Illumina® HumanMethylation450 (450 k) or the HumanMethylationEPIC (EPIC) BeadChip assay at Illumina or cohort-specific laboratories. Cohorts performed sample processing, quality control and normalization as described in Additional file 2: Methods.
We used normalized, untransformed DNAm beta values, ranging from 0 (completely unmethylated) to 1 (completely methylated), after trimming extreme outliers (3 × interquartile range from the quartile limit). Certain individuals were removed due to trimming of extreme values, resulting in different sample sizes across probes and cohorts. We excluded probes mapped to X or Y chromosomes, polymorphic CpGs which overlap with known single-nucleotide polymorphisms, probes with cohort-level call rate < 90%, control probes and cross-reactive probes (targeting repetitive sequences/co-hybridizing to alternate sequences) [20, 21].
Outcomes: parent-reported and actigraphy-estimated sleep measures
Parent-reported sleep duration was chosen as the primary outcome because of known associations with genetic and early-life environmental factors, comparability across studies, and large sample size (n = 3658, for other phenotypes n ≤ 2504) . Sleep duration was the parent-reported number of hours the child slept per night (e.g. "How many hours of sleep does your child get on most nights?") or the calculated difference between reported average time of falling asleep and waking up. This primary outcome was standardized within each cohort (mean = 0, SD = 1).
We also included parent-reported sleep initiation problems (e.g. "The child has difficulty getting to sleep at night") and sleep fragmentation problems (e.g. "The child wakes up more than twice per night") to consider multiple dimensions of sleep [22, 23]. These were used as secondary, binary outcomes.
Actigraphy-estimated sleep outcomes were chosen to objectively assess determined sleep characteristics. Sleep duration (total sleep time), sleep onset latency (time between lying down in bed and falling asleep) and wake-after-sleep-onset duration (time awake between falling asleep and final awakening) were recorded using accelerometery and averaged across the measurement period of 3–9 days. These outcomes were all standardized within each sample (mean = 0, SD = 1). Four cohorts used wrist-worn actigraphs to measure sleep (GENEActiv in Generation R and INMA; Actigraph GT3X in PROGRESS; Actiwatch AW7 in GLAKU). One cohort (CHOP) used armband accelerometers (SenseWear Armband): this cohort only had childhood DNAm, not newborn DNAm data, and thus only contributed to sensitivity analyses (see “Statistical analyses”).
Mean child age at parental assessment ranged between 4 and 12 years between studies. Mean age at actigraphy was 11–12 years in four cohorts and 4.7 years in one cohort (PROGRESS). Please see Table 1 and Additional file 1: Table S1 for cohort-level descriptives and Additional file 2: Methods for further details.
We adjusted for maternal smoking, maternal age, maternal education, child sex, child age at sleep assessment, cell counts, and surrogate variables for batch adjustment. Maternal smoking and education were coded and categorized according to data availability as explained in Additional file 2: Methods. Cell counts were estimated using the Houseman method  with the Bakulski reference panel  (cord blood samples) and Reinius reference panel  (childhood samples). In cord blood DNAm analyses, we additionally adjusted for gestational age at birth, and in childhood DNAm analyses, we adjusted for child age at venepuncture, if different from age at sleep assessment (Table S1). Covariates were chosen based on the previous literature to increase precision and to address potential confounders such as maternal smoking and socio-economic status [2, 8, 27, 28]. Cohort-level analysts were advised to adjust for ethnic differences in any multi-ethnic cohort by choosing the appropriate approach based on cohort characteristics and data availability: these cohort-specific covariates were country of birth (CHOP), self-reported ethnicity (Healthy Start, Viva), or principal components from genome-wide sequencing data (GLAKU, PREDO) (Additional file 2: Methods).
Each cohort-level EWAS was performed according to a predefined analysis plan. We used multivariate linear and logistic regression, for continuous (sleep duration and actigraphy-estimated outcomes) and binary (parent-reported sleep initiation and fragmentation) outcomes, respectively. All p values were two-sided. We excluded participants with missing data and siblings (one sibling per sibling pair, chosen at random). For more cohort-level information, please see Additional file 2: Methods.
Primary meta-analysis: newborn DNAm and parent-reported sleep duration
We examined associations between newborn cord blood DNAm and parent-reported sleep duration at 364,672 loci across 10 cohorts (max n = 3658). We combined 450 k and EPIC data, only including sites that are available on the 450 k array . We performed an inverse variance-weighted fixed-effects meta-analysis using R-3.6.1 (https://www.r-project.org) and METAL (release 2018–08-28) . To assess epigenome-wide statistical inflation, we calculated cohort- and meta-analysis-level genomic inflation factor lambdas (λ) and examined quantile–quantile plots (see Results).
Cohort-level results were meta-analysed at Erasmus MC. We verified the findings in shadow meta-analyses conducted independently at ISGlobal, using the meta-analytical tool GWAMA instead of METAL .
We reran the primary meta-analysis when including only (1) cohorts based in Europe, and (2) cohorts where sleep was assessed on school-aged children (i.e. mean age ≥ 7 years, corresponding to the highest primary school entrance age in the studied populations, http://uis.unesco.org/sites/default/files/documents/indicator-efa-official-entrance-age-to-primary-education.xlsx, accessed 2021/06/09), to increase precision at the cost of sample size.
Individual variations in cord and peripheral blood DNAm are only partly stable from birth to school age . To check if the timing of the DNAm measurement changed the findings, cohorts with available data (re-)ran the EWAS using DNAm data from blood samples collected at school age, and we then reran the meta-analysis on these cohort-level results.
As a final sensitivity analysis, we repeated the primary meta-analysis using a sample-size-weighted meta-analytical approach (also in METAL), which does not expect effect magnitude to be similar across cohorts, to complement the primary inverse variance-weighted approach.
U-shaped or other nonlinear associations were not tested for lack of a strong hypothesis.
Meta-analyses of secondary phenotypes
We ran five secondary meta-analyses that were otherwise similar to the primary meta-analysis but utilized alternative outcomes: parent-reported (1) sleep initiation and (2) fragmentation problems, and actigraphy-estimated (3) sleep duration, (4) sleep-onset latency, and (5) wake-after-sleep-onset duration.
Multiple testing correction and subthreshold findings
The DNAm of nearby sites is correlated, and a cut-off of 2.4 × 10–7 for multiple testing correction on 450 k-based EWASes has been recommended . In the current study, we conducted one primary and five secondary meta-analyses on six separate outcomes: the cut-off of epigenome-level hit was thus (2.4 × 10–7)/6 = 4.0 × 10–8. To confirm independence among the outcomes, we extracted eigenvalues from individual-level matrix of phenotype data using the meff function from poolr (Additional file 1: Figure S1).
Statistical significance cut-offs are inescapably somewhat arbitrary, and particularly in the case of borderline-significant findings, more information on effect magnitude and consistency across studies may be needed to balance type I and II error. We used a suggestive cut-off of p < 5.0 × 10−5 to select the top 25 subthreshold findings that came closest to statistical significance in our primary meta-analysis: this cut-off was selected in line with a recent meta-analysis on maternal anxiety and DNAm at birth . In supplementary tables, we report the site-level results of the primary meta-analysis and sensitivity analyses (as described above) for these sites.
Annotation and description of CpGs
Probes were annotated using meffil (hg19/b37) . Previously reported EWAS-based associations between CpGs of interest and child/adult phenotypes were derived from the EWAS catalogue (http://ewascatalog.org/ accessed 2021/06/09) . Function-related information was derived from Gene-Cards (https://www.genecards.org accessed 2021/06/09) and GWAS Catalog (https://www.ebi.ac.uk/gwas accessed 2021/06/09). For look-up analyses, single-nucleotide polymorphisms related to sleep duration in previous genome-wide-association studies [37,38,39,40,41] were identified using the Sleep Disorder Knowledge Portal (https://sleep.hugeamp.org, accessed 2021/06/23).
We then follow up the hits (none in the primary meta-analysis, two hits in the secondary meta-analyses) to report cross-tissue correspondence at these loci between DNAm in blood vs in the brain, the most relevant organ for sleep phenotype. We used previously published data on DNAm correlations across blood and brain tissue (http://epigenetics.essex.ac.uk/bloodbrain/ accessed 2021/06/09)  and circadian expression of genes of interest in human and mouse tissues (http://circadb.hogeneschlab.org/human, accessed 2021/06/21).
To test for functional enrichment of annotated genes, we performed gene ontology (GO) pathway analyses on genes annotated to CpGs with the lowest p values (p < 1.0 × 10–4) for each model separately (1 primary + 5 secondary meta-analyses), using missMethyl in R . The number of CpGs in these analyses was 44, 62, 51, 60, 43 and 85, for analyses on DNAm at birth and parent-rated sleep duration (primary meta-analysis), initiation problems, fragmentation, actigraphy-estimated sleep duration, sleep onset latency and wake-after-sleep-onset duration, respectively. Analyses were adjusted for the number of pathways tested based on the false discovery rate (FDR) .
We used the output of the primary and secondary meta-analyses to identify DMRs in newborn cord blood associated with parent-reported child sleep duration, initiation and fragmentation, actigraphy-estimated sleep duration, sleep-onset latency and wake-after-sleep-onset duration.
DMRs were analysed using two alternative approaches: DMRcate  and ipDMR . Briefly, DMRcate applies Gaussian kernel smoothing for t-statistics using a bandwidth lambda, and p values for regions are calculated based on the Satterthwaite method and corrected with FDR . On the other hand, ipDMR calculates p values for intervals bordered by two adjacent CpGs, performs the Benjamini–Hochberg (BH) procedure on the interval p values to select those significant intervals at a user-specified FDR threshold (seed threshold), joins all nearby significant intervals and CpGs, recalculates p values for each combined region using the original p values for all CpGs and then obtains FDR-adjusted p values for these regions . In DMRcate, alpha level for FDR as calculated on individual CpGs was 0.05; in ipDMR, FDR alpha level for initial selection of regions (seed) was also set at 0.05 , DMRcate lambda = 500 and C = 5 were set as previously recommended , otherwise default parameters for both approaches were applied.
Associations between newborn DNAm and child sleep were assessed in 11 cohorts, whose participants are described in Table 1. Of these 11 cohorts, 10 had data on parent-reported sleep outcomes and 3 had data on actigraphy-estimated outcomes. Two of these cohorts (Generation R, INMA) had both parent-reported and actigraphy-estimated sleep data: correlations across all different sleep outcomes within these cohorts are shown in Additional file 1: Figure S1.
Associations between childhood DNAm and child sleep were assessed in 8 cohorts, described in Additional file 1: Table S1. Five of these cohorts had DNAm data available among both newborns and in childhood and thus contributed to both analyses.
DNAm and parent-reported sleep duration
We did not observe associations between newborn DNAm and parent-reported child sleep duration in our primary meta-analysis among 3,658 children and a total of 364,672 CpGs, when correcting for multiple testing (cut-off p < 4.0 × 10–8) (Fig. 1A). We found no evidence of genomic inflation (λ = 1.01, quantile–quantile plot in Fig. 1B).
In sensitivity analyses, we repeated the primary meta-analysis among (1) European cohorts only, and (2) among cohorts with a mean age at sleep assessment ≥ 7 years (i.e. school-aged children), to harmonize the phenotype and reduce noise. We also repeated the analyses (3) with DNAm measured in childhood, rather than at birth, as our exposure variable, and reran the primary meta-analysis using (4) a sample-size-weighted, p value-based meta-analysis approach (that does not assume effect magnitude to be similar across cohorts) as an alternative to inverse variance-weighted fixed-effects meta-analysis. In line with the primary meta-analysis, none of these sensitivity analyses suggested any associations between DNAm and parent-reported sleep duration (Additional file 1: Table S3).
In Additional file 1: Table S3, we describe the subthreshold findings that came closest to statistical significance (i.e. the 25 CpG sites with a p value below the suggestive threshold of p < 5 × 10–5). For these subthreshold findings, we show the site-level results from the primary meta-analysis (of newborn DNAm and parent-reported sleep duration), and from analyses where we excluded non-European cohorts, excluded cohorts with child mean age < 7 years, used child (rather than newborn) DNAm as the exposure, or used actigraphy-estimated (rather than parent-reported) sleep duration as the outcome. Briefly, for only 5 out of 25 sites, all cohorts reported a consistent direction of effects in the primary meta-analysis: lower methylation of cg14340131 [annotated to HRAS] and higher methylation of cg01532396 [not annotated to any genes], cg04384689 [not annotated to any genes], cg10143030 [WDR43;SNORD92] and cg17853707 [ZNF91]. One of these sites, cg01532396 was the CpG that came closest to statistical significance in the meta-analysis (p = 2.22 × 10–6). We also investigated if any of these 25 subthreshold findings were located nearby (< 1Mbase) 172 single-nucleotide polymorphisms (SNP) that have previously been associated with sleep duration (https://sleep.hugeamp.org, accessed 2021/06/23). Only one CpG, cg24769432, was located within a 1Mbase region of a sleep-related SNP, rs4538155 (LINC01876 gene located in 2q24.1). Cohort-level results showed inconsistent directions of effect for associations between DNAm at this CpG and sleep duration. Further, none of the 25 subthreshold findings from the primary meta-analysis overlapped with any subthreshold findings (p < 5 × 10–5) from our other (secondary) meta-analyses of newborn DNAm and parent-reported sleep initiation or fragmentation problems.
DNAm and parent-reported sleep initiation and fragmentation problems
We observed no associations between newborn DNAm and parent-reported child sleep initiation (Additional file 1: Figure S2) or fragmentation (Additional file 1: Figure S3) problems among 2,504 and 1,681 children, respectively. We repeated the analyses using DNAm measured at childhood, rather than at birth: again, no associations were observed (Additional file 1: Figure S2–S3).
DNAm and actigraphy-estimated sleep outcomes
Newborn DNAm and actigraphy-estimated sleep data were available among 582 school-aged children from 3 cohorts. Additional file 1: Figure S4, S5 and S6 show meta-analytic associations between newborn DNAm and actigraphy-estimated sleep duration, sleep onset latency (time between lying down in bed and falling asleep) and wake-after-sleep-onset duration (time awake between falling asleep and final awakening), respectively.
We identified two statistically differentiated CpGs (p < 4.0 × 10–8) in analyses of newborn DNAm and child sleep, and none in the analyses of DNAm in childhood and child sleep. These two hits are described briefly below: for details, see Additional file 1: Table S4. First, lower cord blood DNAm at cg24815001 was associated with longer actigraphy-estimated sleep duration among 577 children (p = 3.31 × 10–8). The cohort-level direction of effects was negative in all three cohorts with data, and the meta-analytic effect estimate corresponded to a 0.69-SD-unit reduction in sleep duration z-scores per 10% increase in methylation (i.e. effect estimate − 6.9 per change from completely non-methylated (0) to completely methylated (1), standard error 1.2, I2 = 0). In our study, childhood DNAm at cg24815001 was not associated with actigraphy-based sleep duration (p = 0.84), nor was newborn cord blood DNAm at cg24815001 associated with parent-reported sleep duration (p = 0.96). We are not aware of previous EWAS studies linking cg24815001, an open-sea CpG in chromosome 7 with child (or adult) phenotypes (http://www.ewascatalog.org/?cpg=cg24815001, accessed 2021/06/09), nor is it annotated to any genes . At cg24815001, blood cell DNAm is only weakly correlated with DNAm in brain tissues (r < 0.06) (https://epigenetics.essex.ac.uk/bloodbrain/?probenameg=cg24815001, accessed 2021/06/09).
Second, lower cord blood DNAm at cg02753354 was associated with longer actigraphy-estimated sleep onset latency among 580 children (p = 8.8 × 10–9). However, heterogeneity between our three cohorts was considerable (I2 = 74.1). The results were driven by one cohort (Generation R): negative effect estimates were observed in Generation R (beta = − 11.9, p = 3.37 × 10–9, n = 257) and PROGRESS (beta = − 4.7, p = 0.16, n = 242), while INMA reported a positive effect estimate (beta = 7.7, p = 0.37, n = 81). Further, childhood DNAm of cg02753354 was not associated with sleep onset latency (p = 0.85, n = 712). Newborn DNAm of cg02753354 was not associated with parent-reported sleep initiation problems (p = 0.77, n = 1,712). This CpG in chromosome 19 is annotated to the ARHGAP45 gene, which has a circadian expression profile (http://circadb.hogeneschlab.org/human, last accessed 2021/06/22). ARHGAP45 codes for a precursor of the histocompatibility antigen HA-1 and a GTPase activator for the Rho-type GTPases (https://www.genecards.org/cgi-bin/carddisp.pl?gene=ARHGAP45&keywords=hmha1, last accessed 2021/06/09) and has been linked to endothelial integrity and immune cell maturation [48, 49]. In a previous study investigating the associations between DNA co-methylation modules and sleep in one of the cohorts included in this meta-analysis, namely Generation R, authors reported that a module containing a CpG located at ARHGAP27 was associated with actigraphy-assessed sleep duration . Additionally, a genome-wide association study based on UK Biobank and 23andMe data indicated that variants in ARHGAP27 were associated with self-reported sleep traits, including sleep duration . In previous epigenetic studies, lower cord blood DNAm at cg02753354 has been linked to higher birth weight  and gestational age at birth [52, 53], and a cross-sectional adult study suggested an association with ischaemic stroke . However, correlations between cg02753354 DNAm in blood vs another (perhaps more relevant) target tissue, the brain, are weak (r < 0.11) .
Gene ontology (GO) enrichment
We performed gene ontology (GO) pathway analyses on genes annotated to CpGs with the lowest p values (p < 1.0 × 10–4) (see Methods for details). GO term enrichment analyses did not suggest any GO terms to be over- or underrepresented within the gene set of interest in the primary or secondary meta-analyses.
Differentially methylated region (DMR) analyses
We used two alternative approaches (DMRcate and ipDMR) to identify DMRs in cord blood, with all six child sleep phenotypes. The results are shown in Additional file 1: Table S5. Neither approach identified any DMRs associated with parent-reported child sleep duration, initiation, or fragmentation.
DMRcate yielded 1 DMR including 6 CpGs (chromosome 11: start–end 2,292,890–2,293,048 [annotated to ASCL2]) associated with actigraphy-estimated sleep onset latency. On the other hand, ipDMR yielded 8 regions associated with actigraphy-estimated sleep phenotypes, 6 of which included at least two CpGs: sleep onset latency was associated with 3 DMRs (chr19:1,074,425–1,074,927 [ARHGAP45]; chr6:56,819,612–56,819,616 [BEND6;DST]; and chr17:46,669,566–46,669,645 [LOC404266;HOXB5]), and wake-after-sleep-onset duration with 3 DMRs (chr19:48,894,694–48,894,716 [KDELR1]; chr4:187,422,114–187,422,120 [not annotated to genes]; and chr6:33,048,254–33,048,287 [HLA-DPB1]). None of these DMRs identified using two alternative tools (DMRcate, ipDMR) showed any overlap.
This meta-analysis did not identify associations between DNAm at birth and parent-reported sleep duration among 3,658 children. Similarly, we found no associations between DNAm at birth and parent-reported sleep initiation or fragmentation problems.
Actigraphy can provide a more objective and potentially more sensitive measure of child sleep duration compared to parental reports [55, 56]. Actigraphy can also be used to measure sleep onset latency, an indicator of sleep initiation problems, and wake-after-sleep-onset duration, which can increase when the number of sleep fragments or the average time spent awake between fragments (or both) increases. However, collecting actigraphy data is often relatively cumbersome, compared to questionnaire-based data. The current meta-analysis brought together three cohorts that have both DNAm at birth and actigraphy-estimated child sleep data. Their results showed that lower DNAm at cg24815001 and cg02753354 at birth was associated with longer actigraphy-estimated sleep duration (n = 577) and longer actigraphy-estimated sleep onset latency (n = 580) in childhood, respectively. It is possible that actigraphy-estimated sleep assessments were more sensitive to picking up differences in child sleep that relate to DNAm than parental reports. However, due to the small number of participants and, in the case of cg02753354, heterogeneity between cohorts, we advise caution in accepting these findings as proof of underlying epigenetic signatures or mechanisms without further studies. Further, we identified some potentially differentially methylated regions in relation to actigraphy-estimated sleep onset latency and wake-after-sleep-onset duration, yet the small number of cohorts and the failure to replicate the findings using two alternative tools raise some doubt on whether these reflect meaningful differences in DNAm. Both individual CpG and region-level findings should be further studied in larger paediatric cohorts, should such large-scale data become available.
Our results do not provide consistent evidence that DNAm from blood at birth or in childhood would explain differences in child sleep, neither from parent-reported outcomes, nor actigraphy-based data. Alternative mechanisms that can affect inter-individual differences in child sleep characteristics include, for example, family bedtime practices, socio-economic disparities and genetic make-up . While the results of this large meta-analysis may be less vulnerable to false positives than the limited number of prior smaller, single-sample studies, it is also possible that discrepancies between our and previous studies reflect true biological phenomena. Most importantly, our study focussed on prospective associations with child sleep, while previous studies have mostly studied DNAm and sleep among adults (e.g. [11,12,13]). In a cross-sectional design, DNAm could reflect a biomarker of poor sleep, rather than a causal factor underlying variation in sleep. Further, sleep could be differentially restricted by external factors during childhood compared to adulthood, which in turn could affect the extent of observable epigenetically driven variance. It would be interesting to see if large-scale epigenome-wide studies that span to adolescence and adulthood could establish differences that emerge at this later stage.
DNAm at nearby sites is correlated, and there is disagreement over the optimal method of genome-wide correction. In the current study, however, the findings of the primary meta-analysis would have remained the same (in this case, null) if we had ignored the additional tests introduced by the secondary meta-analyses completely and only applied the more lenient Safari cut-off (2.4 × 10–7),  or if we had used the Bonferroni cut-off that assumes total independence of exposures and outcomes (0.05/364,672/6 = 2.3 × 10–8).
The current study has several strengths. This meta-analysis of previously unpublished EWAS results from a total of 14 cohorts is the largest epigenetic study on child sleep to date. It is extensive in terms of both sample size and the scope of the epigenome-wide rather than candidate site approach. The harmonized, predefined analysis plan increases comparability across studies, while the collaboration between many independent studies improves reliability and generalizability of the results. The variety of sleep outcomes provide a window into several important aspects of child sleep, both objective and subjective. Finally, the inclusion of DNAm data at both birth and in childhood accounts for the partial instability of DNAm across childhood  and can elucidate temporal relations beyond the scope of cross-sectional designs.
Our study also has limitations. First, while this is the largest study so far, it is possible we are still underpowered to identify subtle associations between DNAm and sleep. Further, differences between cohorts, relating for example to age at assessment and cultural and genetic differences between populations, can obscure associations that are specific to certain subgroups, and differences in how sleep was measured (e.g. phrasing of questionnaires) could introduce noise. Second, we combined the 450 k and EPIC array data sets, which could also introduce some noise; however, correlations between these arrays are high for DNAm in blood . Third, we used parent-reported sleep measures as primary outcomes because of comparability across studies, large sample size, and previously reported associations with genetic and early-life environmental factors. However, these measurements are prone to misclassification and may have contributed measurement error to our analyses . We attempted to overcome this limitation by investigating actigraphy-estimated sleep measures as more objective measures, when available. Fourth, we measured DNAm in peripheral blood and not in the brain. Beyond a potential mechanism, blood DNAm could be interpreted as an early-life marker of causal genetic, biological or environmental influences on sleep, thus potentially lending new insights into factors shaping sleep outcomes, rather than biological mechanism. Future research, including advanced causal inference tools and molecular research, would be necessary to establish whether identified hits in blood DNAm are the causal mediator or a proxy of other causal drivers.
Even though we found no consistent evidence of differential methylation related to child sleep, we encourage future research into sleep phenotypes that were beyond the scope of the current study. For example, clinical studies could identify epigenetic patterns associated with clinically diagnosed sleep disorders, which population-based studies may not capture. Future studies should also consider longitudinal associations between DNAm and child sleep. Repeated methylation assessment with concurrent sleep assessments would give insight into trajectories of changes in DNA methylation over time and their long-term effects on sleep. Further, epigenetic mechanisms such as histone modifications or effects limited to specific target tissues (e.g. neuroendocrine tissues) could be of interest, if large-scale analyses of these mechanisms become possible in the future.
In conclusion, we found no consistent evidence of an association between cord blood or peripheral blood DNAm and sleep among children. Larger studies or studies that focus to overcome some of our limitations could reveal subtle associations or confirm associations limited to objectively measured sleep outcomes.
Availability of data and materials
Site-level meta-analytical results are available at the EWAS Catalog [all site-level results from the primary and secondary meta-analyses will be publicly uploaded through EWAS Catalog, http://www.ewascatalog.org/, and link will be added upon acceptance or according to journal embargo policy]. For access to cohort-level data, requests can be sent directly to individual studies.
Meltzer LJ, Williamson AA, Mindell JA. Pediatric sleep health: it matters, and so does how we define it. Sleep Med Rev. 2021;57: 101425.
Stone KC, LaGasse LL, Lester BM, Shankaran S, Bada HS, Bauer CR, et al. Sleep problems in children with prenatal substance exposure: the maternal lifestyle study. Arch Pediatr Adolesc Med. 2010;164(5):452–6.
Pesonen AK, Räikkönen K, Matthews K, Heinonen K, Paavonen JE, Lahti J, et al. Prenatal origins of poor sleep in children. Sleep. 2009;32(8):1086–92.
El-Sheikh M, Kelly RJ, Bagley EJ, Wetter EK. Parental depressive symptoms and children’s sleep: the role of family conflict. J Child Psychol Psychiatr. 2012;53(7):806–14.
Toffol E, Lahti-Pulkkinen M, Lahti J, Lipsanen J, Heinonen K, Pesonen AK, et al. Maternal depressive symptoms during and after pregnancy are associated with poorer sleep quantity and quality and sleep disorders in 3.5-year-old offspring. Sleep Med. 2019;56:201–10.
Barker ED, Walton E, Cecil CAM. Annual Research Review: DNA methylation as a mediator in the association between risk exposure and child and adolescent psychopathology. J Child Psychol Psychiatr. 2018;59(4):303–22.
Breton CV, Marsit CJ, Faustman E, Nadeau K, Goodrich JM, Dolinoy DC, et al. Small-magnitude effect sizes in epigenetic end points are important in children’s environmental health studies: the children’s environmental health and disease prevention research center’s epigenetics working group. Environ Health Perspect. 2017;125(4):511–26.
Joubert BR, Felix JF, Yousefi P, Bakulski KM, Just AC, Breton C, et al. DNA methylation in newborns and maternal smoking in pregnancy: Genome-wide consortium meta-analysis. Am J Hum Genet. 2016;98(4):680–96.
Qureshi IA, Mehler MF. Epigenetics of sleep and chronobiology. Curr Neurol Neurosci Rep. 2014;14(3):432.
Hillary RF, Marioni RE. MethylDetectR: a software for methylation-based health profiling. Wellcome Open Res. 2021;5:283.
Gaine ME, Chatterjee S, Abel T. Sleep deprivation and the epigenome. Front Neural Circuits. 2018;12(February):1–10.
Lahtinen A, Puttonen S, Vanttola P, Viitasalo K, Sulkava S, Pervjakova N, et al. A distinctive DNA methylation pattern in insufficient sleep. Sci Rep. 2019;9(1):1–9.
White AJ, Kresovich JK, Xu Z, Sandler DP, Taylor JA. Shift work, DNA methylation and epigenetic age. Int J Epidemiol. 2019;48(5):1536–44.
Koopman-Verhoeff ME, Mulder RH, Saletin JM, Reiss I, van der Horst GTJ, Felix JF, et al. Genome-wide DNA methylation patterns associated with sleep and mental health in children: a population-based study. J Child Psychol Psychiatr. 2020;61(10):1061–9.
Huang H, Zhu Y, Eliot MN, Knopik VS, McGeary JE, Carskadon MA, et al. Combining human epigenetics and sleep studies in Caenorhabditis elegans: a cross-species approach for finding conserved genes regulating sleep. Sleep. 2017;40(6).
Wong CCY, Parsons MJ, Lester KJ, Burrage J, Eley TC, Mill J, et al. Epigenome-wide DNA methylation analysis of monozygotic twins discordant for diurnal preference. Twin Res Hum Genet. 2015;18(6):662–9.
Jansen EC, Dolinoy DC, O’Brien LM, Peterson KE, Chervin RD, Banker M, et al. Sleep duration and fragmentation in relation to leukocyte DNA methylation in adolescents. Sleep. 2019;42(9):1–9.
Nakahara K, Michikawa T, Morokuma S, Ogawa M, Kato K, Sanefuji M, et al. Association of maternal sleep before and during pregnancy with sleep and developmental problems in 1-year-old infants. Sci Rep. 2021;11(1):11834.
Felix JF, Joubert BR, Baccarelli AA, Sharp GC, Almqvist C, Annesi-Maesano I, et al. Cohort profile: pregnancy and childhood epigenetics (PACE) consortium. Int J Epidemiol. 2018;47:22–3.
Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9.
McCartney DL, Walker RM, Morris SW, McIntosh AM, Porteous DJ, Evans KL. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genomics Data. 2016;9:22–4.
Gregory AM, Sadeh A. Annual research review: sleep problems in childhood psychiatric disorders–a review of the latest science. J Child Psychol Psychiatr Mars. 2016;57(3):296–317.
Meltzer LJ, Biggs S, Reynolds A, Avis KT, Crabtree VM, Bevans KB. The children’s report of sleep patterns-sleepiness scale: a self-report measure for school-aged children. Sleep Med. 2012;13(4):385–9.
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13(1):1–6.
Bakulski KM, Feinberg JI, Andrews SV, Yang J, Brown SL, McKenney S, et al. DNA methylation of cord blood cell types: applications for mixed cell birth studies. Epigenetics. 2016;11(5):354–62.
Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlén SE, Greco D, et al. Differential DNA methylation in purified human blood cells: Implications for cell lineage and studies on disease susceptibility. PLoS ONE. 2012;7(7):e41361.
Bagley EJ, Kelly RJ, Buckhalt JA, El-Sheikh M. What keeps low-SES children from sleeping well: the role of presleep worries and sleep environment. Sleep Med. 2015;16(4):496–502.
McDade TW, Ryan CP, Jones MJ, Hoke MK, Borja J, Miller GE, et al. Genome-wide analysis of DNA methylation in relation to socioeconomic status during development and early adulthood. Am J Phys Anthropol. 2019;169(1):3–11.
Solomon O, MacIsaac J, Quach H, Tindula G, Kobor MS, Huen K, et al. Comparison of DNA methylation measured by Illumina 450K and EPIC BeadChips in blood of newborns and 14-year-old children. Epigenetics. 2018;13(6):655–64.
Willer CJ, Li Y, Abecasis GR. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–1.
Mägi R, Morris AP. GWAMA: software for genome-wide association meta-analysis. BMC Bioinform. 2010;11:288.
Mulder RH, Neumann A, Cecil CAM, Walton E, Houtepen LC, Simpkin AJ, et al. Epigenome-wide change and variation in DNA methylation in childhood: trajectories from birth to late adolescence. Hum Mol Genet. 2021;30(1):119–34.
Saffari A, Silver MJ, Zavattari P, Moi L, Columbano A, Meaburn EL, et al. Estimation of a significance threshold for epigenome-wide association studies. Genet Epidemiol. 2018;42:20–33.
Sammallahti S, Cortes Hidalgo AP, Tuominen S, Malmberg A, Mulder RH, Brunst KJ, et al. Maternal anxiety during pregnancy and newborn epigenome-wide DNA methylation. Mol Psychiatr. 2021;26(6):1832–45.
Min JL, Hemani G, Davey Smith G, Relton C, Suderman M. Meffil: Efficient normalization and analysis of very large DNA methylation datasets. Bioinformatics. 2018;34(23):3983–9.
Battram T, Yousefi P, Crawford G, Prince C, Babaei MS, Sharp G, et al. The EWAS Catalog: a database of epigenome-wide association studies. 2021.
Lane JM, Vlasac I, Anderson SG, Kyle SD, Dixon WG, Bechtold DA, et al. Genome-wide association analysis identifies novel loci for chronotype in 100,420 individuals from the UK Biobank. Nat Commun. 2016;7:1–10.
Lane JM, Liang J, Vlasac I, Anderson SG, Bechtold DA, Bowden J, et al. Genome-wide association analyses of sleep disturbance traits identify new loci and highlight shared genetics with neuropsychiatric and metabolic traits. Nat Genet. 2017;49(2):274–81.
Jones SE, Lane JM, Wood AR, van Hees VT, Tyrrell J, Beaumont RN, et al. Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat Commun. 2019;10(1):1–11.
Wang H, Lane JM, Jones SE, Dashti HS, Ollila HM, Wood AR, et al. Genome-wide association analysis of self-reported daytime sleepiness identifies 42 loci that suggest biological subtypes. Nat Commun. 2019;10(1):1–12.
Dashti HS, Jones SE, Wood AR, Lane JM, van Hees VT, Wang H, et al. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat Commun. 2019;10(1):1–12.
Hannon E, Lunnon K, Schalkwyk L, Mill J. Interindividual methylomic variation across blood, cortex, and cerebellum: Implications for epigenetic studies of neurological and neuropsychiatric phenotypes. Epigenetics. 2015;10(11):1024–32.
Phipson B, Maksimovic J, Oshlack A. MissMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, Lord RV, et al. De novo identification of differentially methylated regions in the human genome. Epigenet Chromatin. 2015;8(6):1–16.
Xu Z, Xie C, Taylor JA, Niu L. IpDMR: Identification of differentially methylated regions with interval P-values. Bioinformatics. 2021;37(5):711–3.
Mallik S, Odom GJ, Gao Z, Gomez L, Chen X, Wang L. An evaluation of supervised methods for identifying differentially methylated regions in Illumina methylation arrays. Brief Bioinform. 2019;20(6):2224–35.
Amado-Azevedo J, Reinhard NR, van Bezu J, van Nieuw Amerongen GP, van Hinsbergh VWM, Hordijk PL. The minor histocompatibility antigen 1 (HMHA1)/ArhGAP45 is a RacGAP and a novel regulator of endothelial integrity. Vascul Pharmacol févr. 2018;101:38–47.
He L, Valignat M, Zhang L, Gelard L, Zhang F, Le Guen V, et al. ARHGAP45 controls naïve T- and B-cell entry into lymph nodes and T-cell progenitor thymus seeding. EMBO Rep. 2021;22(4):e52196.
Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, et al. Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways. Nat Genet Mars. 2019;51(3):394–403.
Küpers LK, Monnereau C, Sharp GC, Yousefi P, Salas LA, Ghantous A, et al. Meta-analysis of epigenome-wide association studies in neonates reveals widespread differential DNA methylation associated with birthweight. Nat Commun. 2019;10(1):1–11.
Merid SK, Novoloaca A, Sharp GC, Küpers LK, Kho AT, Roy R, et al. Epigenome-wide meta-analysis of blood DNA methylation in newborns and children identifies numerous loci related to gestational age. Genome Med. 2020;12(1):1–17.
Bohlin J, Håberg SE, Magnus P, Reese SE, Gjessing HK, Magnus MC, et al. Prediction of gestational age based on genome-wide differentially methylated regions. Genome Biol. 2016;17(1):1–9.
Soriano-Tárraga C, Lazcano U, Giralt-Steinhauer E, Avellaneda-Gómez C, Ois Á, Rodríguez-Campello A, et al. Identification of 20 novel loci associated with ischaemic stroke. Epigenome-wide Assoc Study Epigenet. 2020;15(9):988–97.
Smith MT, McCrae CS, Cheung J, Martin JL, Harrod CG, Heald JL, et al. Use of actigraphy for the evaluation of sleep disorders and circadian rhythm sleep-wake disorders: an American academy of sleep medicine systematic review, meta-analysis, and GRADE assessment. J Clin Sleep Med. 2018;14(7):1209–30.
Galland BC, Short MA, Terrill P, Rigney G, Haszard JJ, Coussens S, et al. Establishing normal values for pediatric nighttime sleep measured by actigraphy: a systematic review and meta-analysis. Sleep. 2018;41(4):1–16.
Newton AT, Honaker SM, Reid GJ. Risk and protective factors and processes for behavioral sleep problems among preschool and early school-aged children: a systematic review. Sleep Med Rev. 2020;52: 101303.
Acknowledgements for each of the participating studies are listed in Additional file 2: Methods.
Statement of significance
This is the first meta-analysis of DNA methylation and child sleep. This study has several methodological strengths, including the harmonized epigenome-wide approach, large sample size, rich parent-reported and actigraphy-based sleep data, and the analysis of both prospective and cross-sectional associations between DNA methylation and sleep. We show that there are no consistent associations of DNA methylation measured from cord blood at birth, or DNA methylation measured from peripheral blood in childhood, and parent-reported sleep duration, sleep initiation problems or fragmentation. Interestingly, however, we identify associations between DNA methylation and actigraphy-based sleep duration and sleep onset latency. This suggests that these objectively measured sleep outcomes could be of particular interest, as we disentangle the basic biological mechanisms responsible for differences in sleep.
Funding for each of the participating studies is listed in Additional file 2: Methods.
Ethics approval and consent to participate
All studies acquired approval from local ethics committees and informed consent from participants (for cohort-level information, please see Additional file 2: Methods).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics of the participating cohorts in analyses of DNAm in childhood and child sleep outcomes. Table S2. Overlap between cord blood and peripheral blood in childhood DNA methylation analyses. Table S3. Site-specific results for the 25 CpGs that came closest to statistical significance (p<5.0×10-5) in the primary meta-analysis of DNAm at birth and parent-reported sleep duration in school age. Table S4. Secondary meta-analyses: significant associations between DNAm and child sleep (p<4.0×10-8). Table S5. Analyses of differentially methylated regions (DMRs) in cord blood at birth and child sleep. Figure S1. Correlations and independence of the six phenotypes of interest. Figure S2. DNAm and parent-reported child sleep initiation problems among school-aged children: Manhattan and quantile–quantile plots. Figure S3. DNAm and parent-reported child sleep fragmentation problems among school-aged children: Manhattan and quantile–quantile plots. Figure S4. DNAm and actigraphy-estimated child sleep duration among school-aged children: Manhattan and quantile–quantile plots. Figure S5. DNAm and actigraphy-estimated child sleep onset latency among school-aged children: Manhattan and quantile–quantile plots. Figure S6. DNAm and actigraphy-estimated child wake-after-sleep-onset duration among school-aged children: Manhattan and quantile–quantile plots.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sammallahti, S., Koopman-Verhoeff, M.E., Binter, AC. et al. Longitudinal associations of DNA methylation and sleep in children: a meta-analysis. Clin Epigenet 14, 83 (2022). https://doi.org/10.1186/s13148-022-01298-4
- Longitudinal studies