Summary Rare variants enriched for functions in chromatin regulation and neuronal synapses have been linked to autism. How chromatin and DNA methylation interact with environmental exposures at synaptic genes in autism etiologies is currently unclear. Using whole genome bisulfite sequencing in brain tissue and a neuronal cell culture model carrying a 15q11.2-q13.3 maternal duplication, we find significant global DNA hypomethylation that is enriched over autism candidate genes and impacts gene expression. The cumulative effect of multiple chromosomal duplications and exposure to the pervasive persistent organic pollutant PCB 95 altered methylation of >1,000 genes. Hypomethylated genes were enriched for H2A.Z, increased maternal UBE3A in Dup15q corresponded to reduced levels of RING1B, and bivalently modified H2A.Z was altered by PCB 95 and duplication. These results demonstrate the compounding effects of genetic and environmental insults on the neuronal methylome that converge upon dysregulation of chromatin and synaptic genes.

Down Syndrome (DS) is the most common genetic cause of intellectual disability, in which an extra copy of human chromosome 21 (HSA21) affects regional DNA methylation profiles across the genome. Although DNA methylation has been previously examined at select regulatory regions across the genome in a variety of DS tissues and cells, differentially methylated regions (DMRs) have yet to be examined in an unbiased sequencing-based approach. Here, we present the first analysis of DMRs from whole genome bisulfite sequencing (WGBS) data of human DS and matched control brain, specifically frontal cortex. While no global differences in DNA methylation were observed, we identified 3,152 DS-DMRs across the entire genome, the majority of which were hypermethylated in DS. DS-DMRs were significantly enriched at CpG islands and de-enriched at specific gene body and regulatory regions. Functionally, the hypermethylated DS-DMRs were enriched for one-carbon metabolism, membrane transport, and glutamatergic synaptic signalling, while the hypomethylated DMRs were enriched for proline isomerization, glial immune response, and apoptosis. Furthermore, in a cross-tissue comparison to previous studies of DNA methylation from diverse DS tissues and reference epigenomes, hypermethylated DS-DMRs showed a strong cross-tissue concordance, while a more tissue-specific pattern was observed for the hypomethylated DS-DMRs. Overall, this approach highlights that low-coverage WGBS of clinical samples can identify epigenetic alterations to known biological pathways, which are potentially relevant to therapeutic treatments and include metabolic pathways. These results also provide new insights into the genome-wide effects of genetic alterations on DNA methylation profiles indicative of altered neurodevelopment and brain function.

Supplementary Figure 2: ASD DMRs heatmap by child outcome continuous measurements of cognition and autism severity versus potential cofounding variables.The plot shows a heatmap of ASD DMRs (y-axis) and the association of % methylation at each DMR with other measured variables. The first 5 child outcome variables on the x-axis includes 4 sub-categories of Mullen scores as well as composite score and autism severity score from the ADOS. Significant associations are red (p < 0.05). While ASD DMRs were highly associated with autism severity and to a lesser degree with early learning Mullen's scores, other potential confounding variables from MARBLES exhibited only rare associations with individual ASD DMRs.

Neonatal dried blood spots (NDBS) are a widely banked sample source that enable retrospective investigation into early-life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early-life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q &lt; 0.05) DMRs that distinguished DS from typical development (TD) and idiopathic developmental delay (DD) were identified. Machine learning feature selection refined these DMRs to 22 loci. The DS DMRs mapped to genes involved in neurodevelopment, metabolism, and transcriptional regulation. Based on comparisons to previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q &lt; 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q &lt; 0.05) enriched for blood-specific chromatin states. A ~ 28 kb block of hypermethylation was observed on chromosome 21 in the RUNX1 locus, which encodes a hematopoietic transcription factor whose binding motif was the most significantly enriched (q &lt; 0.05) overall and specifically within the hypomethylated DMRs. Finally, we also identified DMRs that distinguished DS NDBS based on the presence or absence of congenital heart disease (CHD). Together, these results not only demonstrate the utility of low-pass WGBS on NDBS samples for epigenome-wide association studies, but also provide new insights into the early-life mechanisms of epigenomic dysregulation resulting from trisomy 21.

Background Autism spectrum disorder (ASD) is a neurodevelopmental disorder with complex heritability and higher prevalence in males. The neonatal epigenome has the potential to reflect past interactions between genetic and environmental factors during early development and influence future health outcomes. Methods We performed whole-genome bisulfite sequencing of 152 umbilical cord blood samples from the MARBLES and EARLI high-familial risk prospective cohorts to identify an epigenomic signature of ASD at birth. Samples were split into discovery and replication sets and stratified by sex, and their DNA methylation profiles were tested for differentially methylated regions (DMRs) between ASD and typically developing control cord blood samples. DMRs were mapped to genes and assessed for enrichment in gene function, tissue expression, chromosome location, and overlap with prior ASD studies. DMR coordinates were tested for enrichment in chromatin states and transcription factor binding motifs. Results were compared between discovery and replication sets and between males and females. Results We identified DMRs stratified by sex that discriminated ASD from control cord blood samples in discovery and replication sets. At a region level, 7 DMRs in males and 31 DMRs in females replicated across two independent groups of subjects, while 537 DMR genes in males and 1762 DMR genes in females replicated by gene association. These DMR genes were significantly enriched for brain and embryonic expression, X chromosome location, and identification in prior epigenetic studies of ASD in post-mortem brain. In males and females, autosomal ASD DMRs were significantly enriched for promoter and bivalent chromatin states across most cell types, while sex differences were observed for X-linked ASD DMRs. Lastly, these DMRs identified in cord blood were significantly enriched for binding sites of methyl-sensitive transcription factors relevant to fetal brain development. Conclusions At birth, prior to the diagnosis of ASD, a distinct DNA methylation signature was detected in cord blood over regulatory regions and genes relevant to early fetal neurodevelopment. Differential cord methylation in ASD supports the developmental and sex-biased etiology of ASD and provides novel insights for early diagnosis and therapy.