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.

Health outcomes are frequently shaped by difficult to dissect inter-relationships between biological, behavioral, social and environmental factors. DNA methylation patterns reflect such multivariate intersections, providing a rich source of novel biomarkers and insight into disease etiologies. Recent advances in whole-genome bisulfite sequencing enable investigation of DNA methylation over all genomic CpGs, but existing bioinformatic approaches lack accessible system-level tools. Here, we develop the R package Comethyl, for weighted gene correlation network analysis of user-defined genomic regions that generates modules of comethylated regions, which are then tested for correlations with multivariate sample traits. First, regions are defined by CpG genomic location or regulatory annotation and filtered based on CpG count, sequencing depth and variability. Next, correlation networks are used to find modules of interconnected nodes using methylation values within the selected regions. Each module containing multiple comethylated regions is reduced in complexity to a single eigennode value, which is then tested for correlations with experimental metadata. Comethyl has the ability to cover the noncoding regulatory regions of the genome with high relevance to interpretation of genome-wide association studies and integration with other types of epigenomic data. We demonstrate the utility of Comethyl on a dataset of male cord blood samples from newborns later diagnosed with autism spectrum disorder (ASD) versus typical development. Comethyl successfully identified an ASD-associated module containing regions mapped to genes enriched for brain glial functions. Comethyl is expected to be useful in uncovering the multivariate nature of health disparities for a variety of common disorders. Comethyl is available at github.com/cemordaunt/comethyl with complete documentation and example analyses.

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables 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 < 0.05) differentially methylated regions (DMRs) that distinguished DS from typical development and idiopathic developmental delay 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 with previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q < 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q < 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 < 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.

Rhythmic oscillations of physiological processes depend on integrating the circadian clock and diurnal environment. DNA methylation is epigenetically responsive to daily rhythms, as a subset of CpG dinucleotides in brain exhibit diurnal rhythmic methylation. Here, we show a major genetic effect on rhythmic methylation in a mouse Snord116 deletion model of the imprinted disorder Prader-Willi syndrome (PWS). More than 23,000 diurnally rhythmic CpGs are identified in wild-type cortex, with nearly all lost or phase-shifted in PWS. Circadian dysregulation of a second imprinted Snord cluster at the Temple/Kagami-Ogata syndrome locus is observed at the level of methylation, transcription, and chromatin, providing mechanistic evidence of cross-talk. Genes identified by diurnal epigenetic changes in PWS mice overlapped rhythmic and PWS-specific genes in human brain and are enriched for PWS-relevant phenotypes and pathways. These results support the proposed evolutionary relationship between imprinting and sleep, and suggest possible chronotherapy in the treatment of PWS and related disorders.

Wilson disease (WD) is caused by mutations in the copper transporter ATP7B, leading to copper accumulation in the liver and brain. Excess copper inhibits S-adenosyl-L-homocysteine hydrolase, leading to variable WD phenotypes from widespread alterations in DNA methylation and gene expression. Previously, we demonstrated that maternal choline supplementation in the Jackson toxic milk (tx-j) mouse model of WD corrected higher thioredoxin 1 (TNX1) transcript levels in fetal liver. Here, we investigated the effect of maternal choline supplementation on genome-wide DNA methylation patterns in tx-j fetal liver by whole-genome bisulfite sequencing (WGBS). Tx-j Atp7b genotype-dependent differences in DNA methylation were corrected by choline for genes including, but not exclusive to, oxidative stress pathways. To examine phenotypic effects of postnatal choline supplementation, tx-j mice were randomized to one of six treatment groups: with or without maternal and/or continued choline supplementation, and with or without copper chelation with penicillamine (PCA) treatment. Hepatic transcript levels of TXN1 and peroxiredoxin 1 (Prdx1) were significantly higher in mice receiving maternal and continued choline with or without PCA treatment compared to untreated mice. A WGBS comparison of human WD liver and tx-j mouse liver demonstrated a significant overlap of differentially methylated genes associated with ATP7B deficiency. Further, eight genes in the thioredoxin (TXN) pathway were differentially methylated in human WD liver samples. In summary, Atp7b deficiency and choline supplementation have a genome-wide impact, including on TXN system-related genes, in tx-j mice. These findings could explain the variability of WD phenotype and suggest new complementary treatment options for WD.

DNA methylation acts at the interface of genetic and environmental factors relevant for autism spectrum disorder (ASD). Placenta, normally discarded at birth, is a potentially rich source of DNA methylation patterns predictive of ASD in the child. Here, we performed whole methylome analyses of placentas from a prospective study MARBLES (Markers of Autism Risk in Babies-Learning Early Signs) of high-risk pregnancies. A total of 400 differentially methylated regions (DMRs) discriminated placentas stored from children later diagnosed with ASD compared to typically developing controls. These ASD DMRs were significantly enriched at promoters, mapped to 596 genes functionally enriched in neuronal development, and overlapped genetic ASD risk. ASD DMRs at CYP2E1 and IRS2 reached genome-wide significance, replicated by pyrosequencing and correlated with expression differences in brain. Methylation at CYP2E1 associated with both ASD diagnosis and genotype within the DMR. In contrast, methylation at IRS2 was unaffected by within DMR genotype but modified by preconceptional maternal prenatal vitamin use. This study therefore identified two potentially useful early epigenetic markers for ASD in placenta.

Background

Autism spectrum disorder (ASD) is a neurodevelopmental disorder that affects more than 1% of children in the USA. ASD risk is thought to arise from both genetic and environmental factors, with the perinatal period as a critical window. Understanding early transcriptional changes in ASD would assist in clarifying disease pathogenesis and identifying biomarkers. However, little is known about umbilical cord blood gene expression profiles in babies later diagnosed with ASD compared to non-typically developing and non-ASD (Non-TD) or typically developing (TD) children.

Methods

Genome-wide transcript levels were measured by Affymetrix Human Gene 2.0 array in RNA from cord blood samples from both the Markers of Autism Risk in Babies-Learning Early Signs (MARBLES) and the Early Autism Risk Longitudinal Investigation (EARLI) high-risk pregnancy cohorts that enroll younger siblings of a child previously diagnosed with ASD. Younger siblings were diagnosed based on assessments at 36 months, and 59 ASD, 92 Non-TD, and 120 TD subjects were included. Using both differential expression analysis and weighted gene correlation network analysis, gene expression between ASD and TD, and between Non-TD and TD, was compared within each study and via meta-analysis.

Results

While cord blood gene expression differences comparing either ASD or Non-TD to TD did not reach genome-wide significance, 172 genes were nominally differentially expressed between ASD and TD cord blood (log₂(fold change) > 0.1, p < 0.01). These genes were significantly enriched for functions in xenobiotic metabolism, chromatin regulation, and systemic lupus erythematosus (FDR q < 0.05). In contrast, 66 genes were nominally differentially expressed between Non-TD and TD, including 8 genes that were also differentially expressed in ASD. Gene coexpression modules were significantly correlated with demographic factors and cell type proportions.

Limitations

ASD-associated gene expression differences identified in this study are subtle, as cord blood is not the main affected tissue, it is composed of many cell types, and ASD is a heterogeneous disorder.

Conclusions

This is the first study to identify gene expression differences in cord blood specific to ASD through a meta-analysis across two prospective pregnancy cohorts. The enriched gene pathways support involvement of environmental, immune, and epigenetic mechanisms in ASD etiology.

Background

Wilson disease (WD) is an autosomal recessive disease caused by mutations in ATP7B encoding a copper transporter. Consequent copper accumulation results in a variable WD clinical phenotype involving hepatic, neurologic, and psychiatric symptoms, without clear genotype-phenotype correlations. The goal of this study was to analyze alterations in DNA methylation at the whole-genome level in liver and blood from patients with WD to investigate epigenomic alterations associated with WD diagnosis and phenotype. We used whole-genome bisulfite sequencing (WGBS) to examine distinct cohorts of WD subjects to determine whether DNA methylation could differentiate patients from healthy subjects and subjects with other liver diseases and distinguish between different WD phenotypes.

Results

WGBS analyses in liver identified 969 hypermethylated and 871 hypomethylated differentially methylated regions (DMRs) specifically identifying patients with WD, including 18 regions with genome-wide significance. WD-specific liver DMRs were associated with genes enriched for functions in folate and lipid metabolism and acute inflammatory response and could differentiate early from advanced fibrosis in WD patients. Functional annotation revealed that WD-hypermethylated liver DMRs were enriched in liver-specific enhancers, flanking active liver promoters, and binding sites of liver developmental transcription factors, including Hepatocyte Nuclear Factor 4 alpha (HNF4A), Retinoid X Receptor alpha (RXRA), Forkhead Box A1 (FOXA1), and FOXA2. DMRs associated with WD progression were also identified, including 15 with genome-wide significance. However, WD DMRs in liver were not related to large-scale changes in proportions of liver cell types. DMRs detected in blood differentiated WD patients from healthy and disease control subjects, and distinguished between patients with hepatic and neurologic WD manifestations. WD phenotype DMRs corresponded to genes enriched for functions in mental deterioration, abnormal B cell physiology, and as members of the polycomb repressive complex 1 (PRC1). 44 DMRs associated with WD phenotype tested in a small validation cohort had a predictive value of 0.9.

Conclusions

We identified a disease-mechanism relevant epigenomic signature of WD that reveals new insights into potential biomarkers and treatments for this complex monogenic disease.

Background and aims

Wilson disease (WD) is caused by mutations in the copper transporter ATP7B, with its main pathology attributed to copper-mediated oxidative damage. The limited therapeutic effect of copper chelators and the early occurrence of mitochondrial deficits, however, undermine the prevalence of this mechanism.

Methods

We characterized mitochondrial DNA copy number and mutations as well as bioenergetic deficits in blood from patients with WD and in livers of tx-j mice, a mouse model of hepatic copper accumulation. In vitro experiments with hepatocytes treated with CuSO₄ were conducted to validate in vivo studies.

Results

Here, for the first time, we characterized the bioenergetic deficits in WD as consistent with a mitochondrial DNA depletion-like syndrome. This is evidenced by enriched DNA synthesis/replication pathways in serum metabolomics and decreased mitochondrial DNA copy number in blood of WD patients as well as decreased mitochondrial DNA copy number, increased citrate synthase activity, and selective Complex IV deficit in livers of the tx-j mouse model of WD. Tx-j mice treated with the copper chelator penicillamine, methyl donor choline or both ameliorated mitochondrial DNA damage but further decreased mitochondrial DNA copy number. Experiments with copper-loaded HepG2 cells validated the concept of a direct copper-mitochondrial DNA interaction.

Conclusions

This study underlines the relevance of targeting the copper-mitochondrial DNA pool in the treatment of WD separate from the established copper-induced oxidative stress-mediated damage.

Augmented maternal care during the first postnatal week promotes life-long stress resilience and improved memory compared with the outcome of routine rearing conditions. Recent evidence suggests that this programming commences with altered synaptic connectivity of stress sensitive hypothalamic neurons. However, the epigenomic basis of the long-lived consequences is not well understood. Here, we employed whole-genome bisulfite sequencing (WGBS), RNA-sequencing (RNA-seq), and a multiplex microRNA (miRNA) assay to examine the effects of augmented maternal care on DNA cytosine methylation, gene expression, and miRNA expression. A total of 9,439 differentially methylated regions (DMRs) associated with augmented maternal care were identified in male offspring hypothalamus, as well as a modest but significant decrease in global DNA methylation. Differentially methylated and expressed genes were enriched for functions in neurotransmission, neurodevelopment, protein synthesis, and oxidative phosphorylation, as well as known stress response genes. Twenty prioritized genes were identified as highly relevant to the stress resiliency phenotype. This combined unbiased approach enabled the discovery of novel genes and gene pathways that advance our understanding of the epigenomic mechanisms underlying the effects of maternal care on the developing brain.