- Short report
- Open Access
Case-control meta-analysis of blood DNA methylation and autism spectrum disorder
Molecular Autism volume 9, Article number: 40 (2018)
Several reports have suggested a role for epigenetic mechanisms in ASD etiology. Epigenome-wide association studies (EWAS) in autism spectrum disorder (ASD) may shed light on particular biological mechanisms. However, studies of ASD cases versus controls have been limited by post-mortem timing and severely small sample sizes. Reports from in-life sampling of blood or saliva have also been very limited in sample size and/or genomic coverage. We present the largest case-control EWAS for ASD to date, combining data from population-based case-control and case-sibling pair studies.
DNA from 968 blood samples from children in the Study to Explore Early Development (SEED 1) was used to generate epigenome-wide array DNA methylation (DNAm) data at 485,512 CpG sites for 453 cases and 515 controls, using the Illumina 450K Beadchip. The Simons Simplex Collection (SSC) provided 450K array DNAm data on an additional 343 cases and their unaffected siblings. We performed EWAS meta-analysis across results from the two data sets, with adjustment for sex and surrogate variables that reflect major sources of biological variation and technical confounding such as cell type, batch, and ancestry. We compared top EWAS results to those from a previous brain-based analysis. We also tested for enrichment of ASD EWAS CpGs for being targets of meQTL associations using available SNP genotype data in the SEED sample.
In this meta-analysis of blood-based DNA from 796 cases and 858 controls, no single CpG met a Bonferroni discovery threshold of p < 1.12 × 10− 7. Seven CpGs showed differences at p < 1 × 10− 5 and 48 at 1 × 10− 4. Of the top 7, 5 showed brain-based ASD associations as well, often with larger effect sizes, and the top 48 overall showed modest concordance (r = 0.31) in direction of effect with cerebellum samples. Finally, we observed suggestive evidence for enrichment of CpG sites controlled by SNPs (meQTL targets) among the EWAS CpG hits, which was consistent across EWAS and meQTL discovery p value thresholds.
No single CpG site showed a large enough DNAm difference between cases and controls to achieve epigenome-wide significance in this sample size. However, our results suggest the potential to observe disease associations from blood-based samples. Among the seven sites achieving suggestive statistical significance, we observed consistent, and stronger, effects at the same sites among brain samples. Discovery-oriented EWAS for ASD using blood samples will likely need even larger samples and unified genetic data to further understand DNAm differences in ASD.
The etiology of autism spectrum disorder (ASD) may involve epigenetic mechanisms. Indirect evidence supporting this hypothesis comes from the observation that children with Rett, Fragile X, and Angelman syndromes often show impaired communication and exhibit repetitive behaviors [1,2,3], two core domains affected in autism. All three of these syndromes are caused by epigenetic defects [4,5,6,7,8,9]. Additional evidence stems from genetic studies of rare variation in non-syndromic forms of ASD. Although these studies have primarily identified private variants associated with ASD, there is now strong evidence that the ASD-associated variants converge upon three biological pathways, one of which is chromatin remodeling [10,11,12]. Finally, there is direct evidence from case-control postmortem brain studies supporting epigenetic involvement in ASD. Several candidate gene-based studies have shown altered epigenetic states associated with autism [13,14,15,16,17,18]. Genome-scale screens have identified changes in DNA methylation (DNAm) at specific CpG sites [19, 20] as well as global changes in non-CpG methylation levels  in postmortem cerebral tissue from individuals with ASD relative to controls. Studies of cerebral cortex tissue has revealed genomic spreading of histone H3 lysine 4 methylation and histone H3 lysine 27 acetylation marks, away from the promoter region, among a subset of individuals with ASD compared to controls [22, 23].
Examination of the affected tissue, i.e., the brain, can provide important insights into potential mechanisms of disease etiology; however, there are considerable limitations with these types of studies. They suffer from severely small sample sizes, have historically had low genomic coverage, and often lack comprehensive unified clinical, demographic, and genomic data. Importantly, they are based on autopsy-derived tissue and do not reflect epigenetic marks in a living individual, are not at optimal developmental timing, and may be influenced by life experiences and cause of death. To overcome these barriers, complementary, large population-based autism epigenetic studies using accessible tissues, such as blood, from living individuals are needed. DNAm from peripheral tissues is useful to analyze for its potential to mimic DNAm signatures from the brain. But outside of brain concordance, it can also be useful for its ability reflect consequences of ASD, provided that sample collection is performed following ASD onset. Moreover, peripheral tissue DNAm signatures related to ASD may result from environmental exposures that simultaneously contribute to ASD. To date, three genome-scale epigenetic studies of autism in accessible peripheral tissues have been performed. One study of peripheral blood from 50 monozygotic twin pairs, including 6 pairs discordant for ASD at age 15, examined DNAm at over 27,000 CpG sites in promoter regions. The authors found suggestive evidence for epigenetic alterations associated with ASD and associated traits within families . Similarly, an investigation of DNAm at CpG island regions in lymphoblastoid cell lines, obtained from 7 twin pairs including 3 discordant for ASD, found ASD-related DNAm changes at the RORA gene . Both of these studies were limited by the small number of samples examined, lack of genome-scale coverage, and specific focus on twin pairs with a lack of extension to the general population. Ectoderm cell lineage-derived buccal cells, obtained from 47 ASD cases and 48 controls born to mothers aged 35 and older, have also shown suggestive epigenetic alterations associated with ASD . While suggestive, it is unclear how these buccal-based epigenetic findings relate to a population sample and in a larger number of individuals. Thus, more research in accessible tissues from larger population-based, non-familial, samples is needed.
Here, we overcome previous limitations and perform the largest epigenome-scale examination of DNAm, to date, among two large US case-control studies of autism: the Study to Explore Early Development, phase I (SEED I) and the Simons Simplex Collection (SSC). Both measured DNAm at over 450,000 loci in childhood blood samples from either population-based cases and controls (SEED I) or discordant sibling pairs (SSC). Meta-analysis across both sets included 796 ASD cases and 858 controls. In addition to CpG-specific differential DNAm, we explored the set of blood-derived differentially methylated sites for their concordance in post-mortem brain tissue and their enrichment for genetically controlled CpG sites.
Study to Explore Early Development (SEED)
The Study to Explore Early Development is a multi-site case-control study with population-based ascertainment. In SEED phase 1, a total of 3899 families were recruited across 6 study sites (California, Colorado, Georgia, Maryland, North Carolina, and Pennsylvania) and classified into 3 groups according to a child’s diagnosis: an autism spectrum disorder (ASD) group, a general population control group (POP), and a (non-ASD) developmental delay group. Details regarding participant recruitment, biospecimen collection, and final outcome classification have been previously described [27, 28]. Briefly, eligible children were born in one of the catchment areas between September 1, 2003, and August 31, 2006, which corresponded to being aged 2–5 years at the time of SEED phase I enrollment, resided in the same catchment area at the time of initial contact, and were required to live with a knowledgeable caregiver who could communicate in English (or in English or Spanish in California or Colorado) . Biospecimens were collected when the children were between the ages of 3 and 5 years. Children with possible ASD and DDs were ascertained through multiple sources providing services for children with developmental disorders including hospitals, individual providers, clinics, and education and intervention programs. Parents with a child with an ASD or DD diagnosis could also contact the study directly to enroll. General population controls were ascertained through random sampling of vital records in the catchment areas . This provides a more diverse segment of the population than solely recruiting participants from autism clinics.
Primary caregivers completed the Social Communications Questionnaire (SCQ) , a screener for autism spectrum disorder, during the study invitation phone call. Children with an SCQ score below 11 and without a previous ASD diagnosis were asked to participate in a general developmental evaluation in the clinic using the Mullen Scale of Early Learning (MSEL) . If the SCQ score was above 11, the child had previously received an ASD diagnosis, or a clinician suspected ASD during the clinic visit, the child additionally completed a full ASD evaluation that included the Autism Diagnostic Observation Schedule (ADOS) [31,32,33] and the Autism Diagnostic Interview Revised (ADI-R) [34, 35]. ASD was confirmed based on scores on the ADI-R and ADOS, as described in detail elsewhere . Institutional review boards at each study site and at the Centers for Disease Control and Prevention (CDC) approved the SEED study. Informed consent was obtained from all enrolled participants. For this study, we measured methylation among a subset of SEED phase 1 individuals (n = 980) with genome-wide genotyping data, a complete caregiver interview, an ASD or POP classification, and a sufficient amount of DNA available for methylation measurements.
Simons Simplex Collection (SSC) A complete description of the SSC, which enrolled and collected biospecimens from children and adolescents aged 4–18 years, can be found elsewhere . Briefly, a geneticist and a clinical psychologist were appointed as co-principal investigators at each site. Probands were evaluated with a battery of diagnostic measures, including the Autism Diagnostic Interview – Revised (ADI-R)  and the Autism Diagnostic Observation Schedule (ADOS) . Other instruments provided additional measures of the core features of autism, as well as of intellectual ability (verbal and nonverbal), adaptive behavior, emotional and behavior problems, motor function, and language. A description of instruments employed can be found at https://sfari.org/ssc-instruments. A comprehensive family medical history was obtained that included the proband’s prenatal and perinatal history, developmental milestones, immunizations, medications, dietary supplements, and common behavioral treatments. Emphasis was placed on common “comorbidities” including gastrointestinal complaints, sleep irregularities, and seizures. In addition, questions were asked about genetic, autoimmune, and psychiatric disorders in members of the extended family. Probands were excluded who were younger than 4 years of age or older than 18. Probands were also excluded for conditions that might compromise the validity of diagnostic instruments, such as nonverbal mental age below 18 months, severe neurological deficits, birth trauma, perinatal complications, or genetic evidence of fragile X or Down syndromes. A complete description of exclusion/inclusion criteria can be found at http://sfari.org. Measures of adaptive function, behavior-emotional problems, and symptoms of autism were examined in parents and siblings as well as probands.
Thus, the SSC represents a unique, well-described sample of able children and adolescents with relatively severe ASD, as indicated by ADI-R and ADOS Calibrated Severity Scores .
Reliability of data
To maximize the consistency of clinical observations across sites, each clinician was trained in administration of the ADOS and ADI-R to achieve research reliability as judged by expert clinicians. Most clinicians who had not previously received research training required 4–6 months of practice. Videotapes of interviews were exchanged to ensure that reliability requirements were met and maintained throughout the study. Error rates were very low, averaging less than 0.50 errors/1000 data points. Most errors could be corrected immediately, resulting in an unusually clean data set for a multisite study of this size. During each visit, a blood sample was collected from each study participant and DNA was extracted from blood cells, while plasma was stored for future use.
DNA methylation data quality control (QC) and processing
For the SEED samples, genomic DNA was isolated from 980 whole blood samples using the QIAsymphony midi kit (Qiagen). For each sample, 500 ng of DNA was bisulfite treated using the 96-well EZ DNA methylation kit (Zymo Research). Samples were randomized within and across plates, and across two main processing dates to minimize batch effects, and run on the Illumina HumanMethylation450 BeadChip. Background correction and dye-bias equalization was performed using the function preprocessNoob() [39, 40] in the minfi R package . We included 12 cross-plate duplicates for quality control purposes; pairwise correlation metrics for the duplicate samples ranged from 0.990 to 0.997 with a mean correlation value equal to 0.995. Samples were removed if they had low overall intensity (median unmethylated or methylated signal < 11) or had a detection p value > 0.01 in more than 1% of probes (N = 7), or if reported sex did not match predicted sex generated via the minfi function getSex() (N = 3). Probes were removed if they had a detection p value > 0.01 in more than 10% of samples (n = 702) and then if they had been previously identified as being ambiguously mapped  (n = 29,146). Following QC, the analytic data included DNAm for 455,664 sites on 970 samples. We further removed two samples who were missing a final outcome classification, leaving a total of 453 cases and 515 controls used for association analyses.
For the SSC samples, 500 ng of human genomic DNA was sodium bisulfite-treated for cytosine to thymine conversion using the EZ DNA Methylation Gold kit (Zymo Research). A total of 728 samples (from 364 families) were randomized within and across plates to minimize batch effects and run on the Illumina HumanMethylation450 BeadChip. Additional details have been previously described . Similar quality control procedures as used for the SEED samples were used for the SSC samples. After background correction and dye-bias equalization, samples were removed for low overall intensity (median unmethylated or methylated signal < 11) or for detection p value > 0.01 in more than 1% of probes (N = 42). Probes were removed if they had a detection p value > 0.01 in more than 10% of samples (n = 483) and then if they had been previously identified as being ambiguously mapped (n = 29,213). These steps resulted in an analytic data set with 455,816 sites on 686 samples, consisting of 343 proband-sibling pairs.
Finally, for all SEED and SSC samples, we estimated cell type proportions for six different cell types (granulocytes, monocytes, CD4 T cells, CD8 T cells, B cells, and natural killer cells) using the estimateCellCounts() function in the minfi R package. Estimation incorporated reference data from 60 samples generated from 6 healthy adult men .
Genotype data quality control and processing
Whole genome genotyping data was available for 943 of the 970 SEED 1 samples which passed DNAm quality control steps. After genotype measurement using the Illumina HumanOmni1-Quad BeadChip, standard quality control measures were applied, including removal of samples with < 95% SNP call rate, sex discrepancies, relatedness (Pi-hat > 0.2), or excess hetero- or homozygosity, and removal of markers with < 98.5% call rate, or that were monomorphic. Phasing was performed using SHAPEIT  followed by SNP imputation via the IMPUTE2 software , using 1000 Genomes Project samples as reference. Principal components to account for ancestry were determined via the EigenStrat program .
Epigenome-wide association testing and meta-analysis
For the SEED data, we used linear regression modeling of the M value (the ratio of methylated to total signal determined at every probe in every sample)  as a dependent variable and ASD status, sex, and surrogate variables (SVs) (described below) as independent variables. We implemented this model using the lmFit() function in the limma R package , separately for each of the DNAm probes that passed QC. For the SSC data, we implemented a generalized estimating equation (GEE) model using the gee() function in the gee R package  to account for the correlation inherent to the familial structure in the data. We used a fixed correlation structure of 0.5 for each sibling pair, and regressed M value onto ASD status, sex, and SVs.
To account for sources of technical and biological variability in our association analyses, we estimated surrogate variables (SVs)  in the cleaned SEED and SSC dataset to include as covariates in our downstream analyses. SVs have been shown to capture and adjust for differences related to batch effects and cell type proportions across samples in a wide variety of simulated settings , and to remove the effects of unwanted sources of technical and biological variation . In order to explicitly address the strong confounding effect of sex resulting from the high degree of male bias in ASD diagnosis, we removed sex chromosomes, where DNAm values strongly correlate to sex, before SV estimation, and included sex along with ASD status in the model used for SV estimation. We then used a data-driven procedure individually in the SEED and SSC data to select the number of SVs to include in the association models.
First, to examine the relationship between each SV and known sources of technical and biological variation, we estimated the association between each estimated SV and cell type composition, principal components of genetic ancestry, and processing batch. We then generated a visual representation of the degree of association with these variables using a heat map (Additional file 1: Figures S1a, S2a).
We next examined the influence of iterative inclusion of SVs as adjustment variables in our association regression models. To do this, we first ran a case-control association model with adjustment for the strongest estimated SV , then progressively included the next strongest SV in the analysis and continued this procedure until all estimated SVs were included. For each model, we recorded the inflation factor, or lambda, calculated via the estlambda() function from the GenABEL R package , and visualized the relationship between number of SVs adjusted for and lambda values (Additional file 1: Figures S1b, S2b). We chose the number of SVs to include in the model by considering both the number of SVs at which the estimated lambda values began to plateau and where the known potential confounders appeared to be captured by one or more SVs. We chose to include 19 SVs in the SEED association analysis and 14 in the SSC analysis.
After completing each association analysis, we then performed a meta-analysis using the METAL software  on the 445,068 probes that were present in both the SEED and SSC cleaned datasets. Our approach weighted individual study effects by sample size and also took into account the direction of effect. We also computed the false discovery rate (FDR) using the Benjamini-Hochberg method . We also determined statistical power for this meta-analysis a priori using an estimation method specifically designed for epigenome-wide association studies .
Comparison of blood epigenome-wide association studies (EWAS) hits to brain-based DNAm
We sought to compare the consistency of top EWAS results from the blood-based meta-analysis to our previous analysis of post-mortem brain samples from ASD cases and controls . These data consist of DNAm from three brain regions: cerebellum, prefrontal cortex, and temporal cortex. For the CpG sites reaching suggestive levels of significance (p < 1 × 10− 4) in the meta-analysis, we computed mean differences between cases and controls in each of these three brain regions. We then computed Pearson correlations and quadrant count ratios between the blood effect sizes and three lists of brain effect sizes. We computed quadrant count ratios as the sum of concordant effect sizes (both positive or both negative) minus the sum of discordant effect sizes, all divided by the total number of effect sizes being compared.
Methylation quantitative trait loci (meQTL) query and meQTL target enrichment test
We were interested in exploring the propensity of CpG sites that reached a level of suggestive significance in the EWAS meta-analysis to be significantly associated with nearby SNPs. We used joint DNAm and genotype data to define SNPs associated DNAm, sometimes referred to as “methylation quantitative trait loci (meQTLs)”, and the CpG sites under genetic control, or “meQTL targets”. We then tested for enrichment of meQTL targets in the top ranked CpG sites from the meta-analysis.
In lieu of applying our SV selection method (see “Epigenome-wide association testing and meta-analysis”) to every SNP-CpG association test in the meQTL query, we conducted separate meQTL queries in each processing batch of the SEED data (NBatch1 = 606; NBatch2 = 362). In each batch, we first used a data-driven procedure we have described in detail previously  to select three key parameters for the meQTL query: the SNP minor allele frequency threshold for inclusion, the CpG variability threshold for inclusion, and the maximum distance between SNP and CpG site to be considered for analysis. Briefly, this procedure selects parameters to ensure 80% power to detect a 5% DNAm difference with each addition of the minor allele, at a Bonferroni-defined significance threshold. We then performed the meQTL query in each batch using the MatrixEQTL R package, adjusting for sex, the first five principal components to account for genetic ancestry, and the first two principal components derived from the cell composition estimates. We then defined SNP-CpG association pairs based on results that can gain 100% power in the parameter survey (“permissive”), 90% power (“intermediate”), and 80% power (“stringent”). If a SNP-CpG association pair was significant at a designated threshold in each batch, the CpG site was labeled a meQTL target under that threshold for the downstream enrichment analysis.
We tested for enrichment of meQTL targets in ASD-associated CpG sites. We examined this using two ASD EWAS meta-analysis p value thresholds (p < 1 × 10− 3 and p < 1 × 10− 4) and the three meQTL p value thresholds. In each enrichment test, we accounted for the two main features of CpG sites likely to affect results: the degree of variability in DNAm at that CpG and the number of SNPs in the boundary considered. To do this, we binned each CpG site by decile according to these factors. For each EWAS/meQTL threshold scenario, we compared the proportion of meQTL targets among ASD-related CpGs to a null distribution of randomly selected CpGs, equal in count to the number of ASD-associated sites, matched on the same variability and nearby SNP decile. We defined a fold enrichment statistic as the count of meQTL targets in the ASD-associated CpGs divided by the mean proportion of meQTL targets from the null set, and an enrichment p value as the number of null CpG sets with a count of meQTL targets that was equal to or exceeded the count in the ASD-associated CpG list.
ASD EWAS meta-analysis in blood
We performed a meta-analysis over the 445,068 probes that were present in both the SEED (Table 1) and SSC (Additional file 2: Table S1) cleaned datasets. Figure 1a shows the range of p values and effect sizes detected in our meta-analysis. No CpG sites reached a Bonferroni level of significance, and effect sizes were modest (1.12 × 10− 7; Fig. 1a). The genomic inflation factor (λ) was 1.03, with a slight separation from expectation at the tail (Fig. 1b). A total of 48 CpG sites met or exceeded a p value < 1 × 10− 4, and 7 CpG sites (Table 2) reached a significance level of p < 1 × 10− 5. We have provided a full list of summary statistics for both SEED and SSC for all 445,068 probes (Additional file 3: Table S2). Based on our analytic sample size, we had 80% power to detect a 3.8% DNAm difference between cases and controls at a Bonferonni level of significance.
Consistency of blood EWAS hits in brain
We considered the consistency of signal for the 48 blood-based CpGs with suggestive significance (p value < 1 × 10− 4), among results from three different brain regions with data available from our previous analysis of post-mortem brain samples and ASD . The cerebellum exhibited a moderate degree of concordance in effect size and direction (r = 0.31; QCR = 0.33), although prefrontal cortex (r = 0.02; QCR = 0.125) and temporal cortex (r = − 0.10; QCR = − 0.125) showed only minimal concordance (Additional file 4: Table S3, Additional file 5: Figure S3). When considering the seven CpG sites with more stringent blood-based p values < 1 × 10− 5, the direction of effect was consistent for at least five of these seven in all three brain region results, with typically larger effect sizes (Table 3). The CpG site with the largest effect size in blood (cg09671955) displayed consistent, and larger, effect sizes in all three brain regions (Table 2, Additional file 4: Table S3).
meQTL target enrichment test
When considering all CpGs associated with ASD at a liberal p < 1 × 10− 3 EWAS threshold, we found meQTL target enrichment at nominal significance (penrichment = 0.041), but not at a significance threshold corrected for the six tests performed overall (Table 4). All other combinations of EWAS and meQTL p values tested displayed similar suggestive levels of significance (0.089 ≤ penrichment ≤ 0.243) and a consistent direction of effect towards enrichment. Also, tests conducted for CpGs meeting the more stringent EWAS p value threshold (1 × 10− 4) displayed a consistently greater effect size than their corresponding tests from the more liberal EWAS threshold.
We report results from a large study investigating the relationship between ASD and DNAm. A case-control meta-analysis of peripheral blood samples from the Study to Explore Early Development and the Simons Simplex Collection revealed that none of the 455,068 CpG sites investigated were associated with ASD at a genome-wide significant threshold. However, 48 CpG sites reached suggestive significance levels at p < 1 × 10− 4, including 7 CpGs at p < 1 × 10− 5. Associations with ASD at these sites display moderate concordance with post-mortem brain sample results from a previous study and display suggestive evidence for enrichment of SNP-controlled CpG sites, or meQTL targets.
Blood-based DNAm signatures of ASD may reflect causal mechanisms that mimic signatures seen in the brain. However, case-control sampling, by design occurring after disease onset, can also reflect consequences of ASD or of environmental exposures. Nonetheless, given the potential involvement of epigenetic mechanisms in ASD [4,5,6,7,8,9,10,11,12], and the availability of blood samples from existing studies, this was an important project to pursue. In contrast to our null findings in blood samples, previous work using brain samples has shown specific DNAm to be associated with ASD [15, 19, 20]. The difference in tissue type—specifically that brain-based signatures better reflect ASD etiology—might explain the inconsistent results across these studies. Nonetheless, a previous study of ASD and DNAm in peripheral blood  using the earlier 27K Illumina array reported numerous differentially methylated sites. Our study, despite a much higher sample size, did not observe associations at these same sites. However, the previous study’s reported differentially methylated sites were based on a ranking that integrated both degree of statistical significance and effect size; none of their single CpG associations achieved statistical significance at a Bonferroni correction level. Also, the previous study did not explicitly account for potential confounding by cell type proportions or address potential batch effects in the analytic pipeline, as this was not yet commonplace in EWAS pipelines at the time of their report. We used a rigorous, data-driven method to account for these factors and control the genome-wide type I error rate. It is also possible that true differentially methylated positions for ASD exist in blood but have smaller effects sizes than we were powered to detect. Recent EWAS discoveries have indeed found replicable very small effect sizes, particularly in environmental health .
Ideally, tissue-relevant samples from earliest in utero development would provide the greatest mechanistic insights. However, this is not typically available and our study is based instead on a case-control design using blood-based samples. Recent examples of blood-based epigenetic work in brain-based disorders have shown utility, particularly when focused on integration of genetic and epigenetic variation or when inferring biological pathways rather than discovery of specific CpGs [57, 59,60,61,62,63]. For example, DNAm under genetic control (meQTL targets) can inform genetic associations observed for ASD. Further, consideration of aggregated sets of CpGs associated with ASD, rather than single sites, can elucidate pathways of interest . In our own analysis, the most differential blood-based CpGs had consistent effect sizes and directions, although weaker, with brain-based results, particularly for cerebellum. Blood-based CpGs were also moderately enriched (though not to a statistically significant extent) for meQTL targets. These results suggest that blood DNAm can be reflective of DNAm in affected tissues and suggest genetic control of DNAm as a mechanism for this occurrence, at least in an ASD context. More precise evidence is needed, but given the easy accessibility of blood for DNAm measurements versus brain , the utility of blood-based DNAm research in ASD is worthy of additional consideration.
In summary, our results suggest a potential role of genetic factors in contributing to DNAm differences in ASD. This is consistent with a recently reported meta-analysis combining these results with an additional European case control sample. In that report, while ASD itself was not associated with DNAm at a genome-wide significance threshold, polygenic risk scores for ASD were associated with DNAm at particular CpGs . In contrast, we show enrichment for meQTL targets (DNAm at least partially controlled by SNP variation) among CpGs with at least marginal association with ASD. We further characterize these ASD DNAm associations and compare to brain-based results. Together, our work and the report by Hannon et al. suggest more investigation of the relationship between genetic risk, epigenetic marks, and ASD is warranted in larger samples and across multiple tissues and designs. To this end, we have provided our full summary statistics and meta-analysis results. The need for greater sample sizes mimics the initial stages of genetic variation discovery in ASD, for which large mega-analyses are starting to pay dividends .
Autism spectrum disorder
Epigenome-wide association study
Methylation quantitative trait loci
Study to Explore Early Development
Single nucleotide polymorphism
Simons Simplex Collection
Jedele KB. The overlapping spectrum of Rett and Angelman syndromes: a clinical review. Semin Pediatr Neurol. 2007;14:108–17.
Mount RH, Charman T, Hastings RP, Reilly S, Cass H. Features of autism in Rett syndrome and severe mental retardation. J Autism Dev Disord. 2003;33:435–42.
Budimirovic DB, Kaufmann WE. What can we learn about autism from studying fragile X syndrome? Dev Neurosci. 2011;33:379–94.
Amir RE, Van den Veyver IB, Wan M, Tran CQ, Francke U, Zoghbi HY. Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. Nat Genet. 1999;23:185–8.
Wan M, Lee SS, Zhang X, Houwink-Manville I, Song HR, Amir RE, et al. Rett syndrome and beyond: recurrent spontaneous and familial MECP2 mutations at CpG hotspots. Am J Hum Genet. 1999;65:1520–9.
Sutcliffe JS, Nakao M, Christian S, Orstavik KH, Tommerup N, Ledbetter DH, et al. Deletions of a differentially methylated CpG island at the SNRPN gene define a putative imprinting control region. Nat Genet. 1994;8:52–8.
Kishino T, Lalande M, Wagstaff J. UBE3A/E6-AP mutations cause Angelman syndrome. Nat Genet. 1997;15:70–3.
Bell MV, Hirst MC, Nakahori Y, MacKinnon RN, Roche A, Flint TJ, et al. Physical mapping across the fragile X: hypermethylation and clinical expression of the fragile X syndrome. Cell. 1991;64:861–6.
Vincent A, Heitz D, Petit C, Kretz C, Oberlé I, Mandel JL. Abnormal pattern detected in fragile-X patients by pulsed-field gel electrophoresis. Nature. 1991;349:624–6.
Sanders SJ, He X, Willsey AJ, Ercan-Sencicek AG, Samocha KE, Cicek AE, et al. Insights into autism spectrum disorder genomic architecture and biology from 71 risk loci. Neuron. 2015;87:1215–33.
Pinto D, Delaby E, Merico D, Barbosa M, Merikangas A, Klei L, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94:677–94.
Krumm N, O’Roak BJ, Shendure J, Eichler EE. A de novo convergence of autism genetics and molecular neuroscience. Trends Neurosci. 2014;37:95–105.
Nagarajan RP, Hogart AR, Gwye Y, Martin MR, LaSalle JM. Reduced MeCP2 expression is frequent in autism frontal cortex and correlates with aberrant MECP2 promoter methylation. Epigenetics. 2006;1:e1–11.
Gregory SG, Connelly JJ, Towers AJ, Johnson J, Biscocho D, Markunas CA, et al. Genomic and epigenetic evidence for oxytocin receptor deficiency in autism. BMC Med. 2009;7:62.
James SJ, Shpyleva S, Melnyk S, Pavliv O, Pogribny IP. Elevated 5-hydroxymethylcytosine in the Engrailed-2 (EN-2) promoter is associated with increased gene expression and decreased MeCP2 binding in autism cerebellum. Transl Psychiatry. 2014;4:e460.
James SJ, Shpyleva S, Melnyk S, Pavliv O, Pogribny IP. Complex epigenetic regulation of engrailed-2 (EN-2) homeobox gene in the autism cerebellum. Transl Psychiatry. 2013;3:e232.
Zhu L, Wang X, Li X-L, Towers A, Cao X, Wang P, et al. Epigenetic dysregulation of SHANK3 in brain tissues from individuals with autism spectrum disorders. Hum Mol Genet. 2014;23:1563–78.
Mor M, Nardone S, Sams DS, Elliott E. Hypomethylation of miR-142 promoter and upregulation of microRNAs that target the oxytocin receptor gene in the autism prefrontal cortex. Mol Autism. 2015;6:46.
Ladd-Acosta C, Hansen KD, Briem E, Fallin MD, Kaufmann WE, Feinberg AP. Common DNA methylation alterations in multiple brain regions in autism. Mol Psychiatry. 2014;19:862–71.
Nardone S, Sams DS, Reuveni E, Getselter D, Oron O, Karpuj M, et al. DNA methylation analysis of the autistic brain reveals multiple dysregulated biological pathways. Transl Psychiatry. 2014;4:e433.
Ellis SE, Gupta S, Moes A, West AB, Arking DE. Exaggerated CpH methylation in the autism-affected brain. Mol Autism. 2017;8:6.
Shulha HP, Cheung I, Whittle C, Wang J, Virgil D, Lin CL, et al. Epigenetic signatures of autism: trimethylated H3K4 landscapes in prefrontal neurons. Arch Gen Psychiatry. 2012;69:314–24.
Sun W, Poschmann J, Cruz-Herrera Del Rosario R, Parikshak NN, Hajan HS, Kumar V, et al. Histone acetylome-wide association study of autism spectrum disorder. Cell 2016;167:1385–1397.e11.
Wong CCY, Meaburn EL, Ronald A, Price TS, Jeffries AR, Schalkwyk LC, et al. Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits. Mol Psychiatry. 2014;19:495–503.
Nguyen A, Rauch TA, Pfeifer GP, Hu VW. Global methylation profiling of lymphoblastoid cell lines reveals epigenetic contributions to autism spectrum disorders and a novel autism candidate gene, RORA, whose protein product is reduced in autistic brain. FASEB J Off Publ Fed Am Soc Exp Biol. 2010;24:3036–51.
Berko ER, Suzuki M, Beren F, Lemetre C, Alaimo CM, Calder RB, et al. Mosaic epigenetic dysregulation of ectodermal cells in autism spectrum disorder. PLoS Genet. 2014;10:e1004402.
Schendel DE, Diguiseppi C, Croen LA, Fallin MD, Reed PL, Schieve LA, et al. The Study to Explore Early Development (SEED): a multisite epidemiologic study of autism by the Centers for Autism and Developmental Disabilities Research and Epidemiology (CADDRE) network. J Autism Dev Disord. 2012;42:2121–40.
Wiggins LD, Levy SE, Daniels J, Schieve L, Croen LA, DiGuiseppi C, et al. Autism spectrum disorder symptoms among children enrolled in the Study to Explore Early Development (SEED). J Autism Dev Disord. 2015;45:3183–94.
Rutter M, Bailey A, Lord C. The social communication questionnaire: manual. Western Psychological Services. West Psychol Serv. 2003;
Mullen EM. Mullen Scales of Early Learning. Am Guid Serv Inc. 1995;
Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, et al. The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000;30:205–23.
Gotham K, Risi S, Pickles A, Lord C. The Autism Diagnostic Observation Schedule: revised algorithms for improved diagnostic validity. J Autism Dev Disord. 2007;37:613–27.
Lord C, Rutter M, DeLavore P, Risi S. Autism Diagnostic Observation Schedule-WPS (ADOS-WPS). West Psychol Serv. 1999;
Lord C, Rutter M, Le Couteur A. Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994;24:659–85.
Rutter M, Le Couteur A, Lord C. ADI-R: Autism diagnostic interview-revised: manual. Western Psychological Services. West Psychol Serv. 2003;
Wiggins LD, Piazza V, Robins DL. Comparison of a broad-based screen versus disorder-specific screen in detecting young children with an autism spectrum disorder. Autism Int J Res Pract. 2014;18:76–84.
Fischbach GD, Lord C. The Simons Simplex Collection: a resource for identification of autism genetic risk factors. Neuron. 2010;68:192–5.
Gotham K, Pickles A, Lord C. Standardizing ADOS scores for a measure of severity in autism spectrum disorders. J Autism Dev Disord. 2009;39:693–705.
Fortin J-P, Labbe A, Lemire M, Zanke BW, Hudson TJ, Fertig EJ, et al. Functional normalization of 450k methylation array data improves replication in large cancer studies. Genome Biol. 2014;15(12):503.
Triche TJ, Weisenberger DJ, Van Den Berg D, Laird PW, Siegmund KD. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucleic Acids Res. 2013;41:e90.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinforma Oxf Engl. 2014;30:1363–9.
Chen Y, 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:203–9.
Alisch RS, Barwick BG, Chopra P, Myrick LK, Satten GA, Conneely KN, et al. Age-associated DNA methylation in pediatric populations. Genome Res. 2012;22:623–32.
Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlén S-E, 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:e41361.
Delaneau O, Zagury J-F, Marchini J. Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods. 2013;10:5–6.
Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44:955–9.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–9.
Du P, Zhang X, Huang C-C, Jafari N, Kibbe WA, Hou L, et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics. 2010;11:587.
limma. Bioconductor. http://bioconductor.org/packages/limma/. Accessed 21 Sep 2016.
Author VJCP to R by TL (versions 3 13 and 4 4) and BR (version 4 13) N that maintainers are not available to give advice on using a package they did not. gee: Generalized Estimation Equation Solver. 2015. https://CRAN.R-project.org/package=gee. Accessed 19 Jun 2018.
Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007;3:1724–35.
McGregor K, Bernatsky S, Colmegna I, Hudson M, Pastinen T, Labbe A, et al. An evaluation of methods correcting for cell-type heterogeneity in DNA methylation studies. Genome Biol. 2016;17:84.
GenABEL.pdf. http://www.genabel.org/packages/GenABEL. Accessed 15 Mar 2017.
Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinforma Oxf Engl. 2010;26:2190–1.
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:289–300.
Tsai P-C, Bell JT. Power and sample size estimation for epigenome-wide association scans to detect differential DNA methylation. Int J Epidemiol. 2015;
Andrews SV, Ellis SE, Bakulski KM, Sheppard B, Croen LA, Hertz-Picciotto I, et al. Cross-tissue integration of genetic and epigenetic data offers insight into autism spectrum disorder. Nat Commun. 2017;8:1011.
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:511–26.
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:1024–32.
Davies MN, Volta M, Pidsley R, Lunnon K, Dixit A, Lovestone S, et al. Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood. Genome Biol. 2012;13:R43.
Bakulski KM, Halladay A, Hu VW, Mill J, Fallin MD. Epigenetic research in neuropsychiatric disorders: the “tissue issue”. Curr Behav Neurosci Rep. 2016;3:264–74.
Montano C, Taub MA, Jaffe A, Briem E, Feinberg JI, Trygvadottir R, et al. Association of DNA methylation differences with schizophrenia in an epigenome-wide association study. JAMA Psychiatry. 2016;73:506–14.
Edgar RD, Jones MJ, Meaney MJ, Turecki G, Kobor MS. BECon: a tool for interpreting DNA methylation findings from blood in the context of brain. Transl Psychiatry. 2017;7:e1187.
Hannon E, Schendel D, Ladd-Acosta C, Grove J, iPSYCH-Broad ASD Group, Hansen CS, et al. Elevated polygenic burden for autism is associated with differential DNA methylation at birth. Genome Med. 2018;10:19.
Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol Autism. 2017;8:21.
Johnson SB, Whitney G, McAuliffe M, Wang H, McCreedy E, Rozenblit L, et al. Using global unique identifiers to link autism collections. J Am Med Inform Assoc JAMIA. 2010;17:689–95.
This project was supported by Centers for Disease Control and Prevention (CDC) Cooperative Agreements announced under the following RFAs: 01086, 02199, DD11-002, DD06-003, DD04-001, and DD09-002. The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention. The DNA methylation assays were supported by Autism Speaks Award #7659 and the genotype assays were supported by NIEHS (R01ES019001; R01ES017646). S. Andrews was supported by the Burroughs-Wellcome Trust training grant: Maryland, Genetics, Epidemiology and Medicine (MD-GEM). The SSC was supported by Simons Foundation (SFARI) award and NIH grant MH089606, both awarded to S.T. Warren.
Availability of data and materials
The SEED I data generated and analyzed for this study are not publicly available due to lack of explicit consent for such sharing in the written informed consents for SEED sites, according to the CDC IRB that governs the SEED network. The SEED data generated and analyzed for this study are not publicly available due to confidentiality and informed consent requirements. We have submitted the raw 450K methylation data from the SSC samples to the National Database for Autism Research (NDAR; ndar.nih.gov), which can be found under collection #300.
Ethics approval and consent to participate
This study was approved by the institutional review boards at each SEED site: SEED 1 recruitment was approved by the IRBs of each recruitment site: Institutional Review Board (IRB)-C, CDC Human Research Protection Office; Kaiser Foundation Research Institute (KFRI) Kaiser Permanente Northern California IRB, Colorado Multiple IRB, Emory University IRB, Georgia Department of Public Health IRB, Maryland Department of Health and Mental Hygiene IRB, Johns Hopkins Bloomberg School of Public Health Review Board, University of North Carolina IRB and Office of Human Research Ethics, IRB of The Children’s Hospital of Philadelphia, and IRB of the University of Pennsylvania. All enrolled families provided written consent for participation. This methylation substudy was approved as an amendment of the Johns Hopkins Institutional Review Board (IRB) approval. For participants from the Simons simplex collection, parents consented and children assented as required by each local institutional review board, which included a coalition of clinics located at Michigan, Yale, Emory, Columbia, Vanderbilt, McGill Washington, and Harvard Universities (Children’s Hospital of Boston), and at the Universities of Washington, Illinois (Chicago), Missouri, UCLA, and the Baylor College of Medicine. To protect the privacy of participants, Global Unique Identifiers (GUIDs) were constructed from personal information using an algorithm devised in collaboration with scientists at the NIH . Each clinic retained personal identifiers on site and transmitted de-identified GUIDs to a central database, as described previously .
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.
Figures S1-S2. Depiction of surrogate variable selection process for SEED (S1) and SSC (S2). Panel A: Heatmap indicating degree of association with known potential technical variables or confounders with estimated surrogate variables. Panel B: Inflation factor (lambda) calculated for progressively including surrogate variables in association models. The number of surrogate variables to include in the ultimate association testing model was to determine to be that which properly controlled the inflation factor and adequately captured known technical variables or confounders. See “Methods” for additional explanation. (PDF 19 kb)
Table S1. Demographic characteristics for samples in the SSC (S2) dataset. (XLSX 11 kb)
Table S2. Full summary statistics and meta-analysis results for all 445,608 CpG sites that were present in both the cleaned SEED and SSC datasets. (CSV 40707 kb)
Table S3. Concordance between suggestively associated (p < 1 × 10− 4) CpG sites in peripheral blood and their corresponding effect sizes in three brain regions. (CSV 3 kb)
Figures S3. Quadrant plots depicting concordance in effect sizes between suggestively associated (p < 1 × 10− 4) CpG sites in peripheral blood and three brain regions. A) Prefrontal cortex B) Temporal Cortex C) Cerebellum. Points in red indicate those sites with p < 1 × 10− 5 in peripheral blood. (PNG 21 kb)
About this article
- DNA methylation
- Autism spectrum disorders
- Peripheral blood
- Study to Explore Early Development
- Simons Simplex Collection