Networks of placental DNA methylation correlate with maternal serum PCB concentrations and child neurodevelopment.

BACKGROUND: Gestational exposure to polychlorinated biphenyls (PCBs) has been associated with elevated risk for neurodevelopmental disorders. Placental epigenetics may serve as a potential mechanism of risk or marker of altered placental function. Prior studies have associated differential placental DNA methylation with maternal PCB exposure or with increased risk of autism spectrum disorder (ASD). However, sequencing-based placental methylomes have not previously been tested for simultaneous associations with maternal PCB levels and child neurodevelopmental outcomes. OBJECTIVES: We aimed to identify placental DNA methylation patterns associated with maternal PCB levels and child neurodevelopmental outcomes in the high-risk ASD MARBLES cohort. METHODS: We measured 209 PCB congeners in 104 maternal serum samples collected at delivery. We identified networks of DNA methylation from 147 placenta samples using the Comethyl R package, which performs weighted gene correlation network analysis for whole genome bisulfite sequencing data. We tested placental DNA methylation modules for association with maternal serum PCB levels, child neurodevelopment, and other participant traits. RESULTS: PCBs 153 + 168, 170, 180 + 193, and 187 were detected in over 50% of maternal serum samples and were highly correlated with one another. Consistent with previous findings, maternal age was the strongest predictor of serum PCB levels, alongside year of sample collection, pre-pregnancy BMI, and polyunsaturated fatty acid levels. Twenty seven modules of placental DNA methylation were identified, including five which significantly correlated with one or more PCBs, and four which correlated with child neurodevelopment. Two modules associated with maternal PCB levels as well as child neurodevelopment, and mapped to CSMD1 and AUTS2, genes previously implicated in ASD and identified as differentially methylated regions in mouse brain and placenta following gestational PCB exposure. CONCLUSIONS: Placental DNA co-methylation modules were associated with maternal PCBs and child neurodevelopment. Methylation of CSMD1 and AUTS2 could be markers of altered placental function and/or ASD risk following maternal PCB exposure.

Placental methylome reveals a 22q13.33 brain regulatory gene locus associated with autism.

BACKGROUND: Autism spectrum disorder (ASD) involves complex genetics interacting with the perinatal environment, complicating the discovery of common genetic risk. The epigenetic layer of DNA methylation shows dynamic developmental changes and molecular memory of in utero experiences, particularly in placenta, a fetal tissue discarded at birth. However, current array-based methods to identify novel ASD risk genes lack coverage of the most structurally and epigenetically variable regions of the human genome. RESULTS: We use whole genome bisulfite sequencing in placenta samples from prospective ASD studies to discover a previously uncharacterized ASD risk gene, LOC105373085, renamed NHIP. Out of 134 differentially methylated regions associated with ASD in placental samples, a cluster at 22q13.33 corresponds to a 118-kb hypomethylated block that replicates in two additional cohorts. Within this locus, NHIP is functionally characterized as a nuclear peptide-encoding transcript with high expression in brain, and increased expression following neuronal differentiation or hypoxia, but decreased expression in ASD placenta and brain. NHIP overexpression increases cellular proliferation and alters expression of genes regulating synapses and neurogenesis, overlapping significantly with known ASD risk genes and NHIP-associated genes in ASD brain. A common structural variant disrupting the proximity of NHIP to a fetal brain enhancer is associated with NHIP expression and methylation levels and ASD risk, demonstrating a common genetic influence. CONCLUSIONS: Together, these results identify and initially characterize a novel environmentally responsive ASD risk gene relevant to brain development in a hitherto under-characterized region of the human genome.

Comethyl: a network-based methylome approach to investigate the multivariate nature of health and disease.

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.

Low-pass whole genome bisulfite sequencing of neonatal dried blood spots identifies a role for RUNX1 in Down syndrome DNA methylation profiles.

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.

Expression Changes in Epigenetic Gene Pathways Associated With One-Carbon Nutritional Metabolites in Maternal Blood From Pregnancies Resulting in Autism and Non-Typical Neurodevelopment.

The prenatal period is a critical window for the development of autism spectrum disorder (ASD). The relationship between prenatal nutrients and gestational gene expression in mothers of children later diagnosed with ASD or non-typical development (Non-TD) is poorly understood. Maternal blood collected prospectively during pregnancy provides insights into the effects of nutrition, particularly one-carbon metabolites, on gene pathways and neurodevelopment. Genome-wide transcriptomes were measured with microarrays in 300 maternal blood samples in Markers of Autism Risk in Babies-Learning Early Signs. Sixteen different one-carbon metabolites, including folic acid, betaine, 5'-methyltretrahydrofolate (5-MeTHF), and dimethylglycine (DMG) were measured. Differential expression analysis and weighted gene correlation network analysis (WGCNA) were used to compare gene expression between children later diagnosed as typical development (TD), Non-TD and ASD, and to one-carbon metabolites. Using differential gene expression analysis, six transcripts (TGR-AS1, SQSTM1, HLA-C, and RFESD) were associated with child outcomes (ASD, Non-TD, and TD) with genome-wide significance. Genes nominally differentially expressed between ASD and TD significantly overlapped with seven high confidence ASD genes. WGCNA identified co-expressed gene modules significantly correlated with 5-MeTHF, folic acid, DMG, and betaine. A module enriched in DNA methylation functions showed a suggestive protective association with folic acid/5-MeTHF concentrations and ASD risk. Maternal plasma betaine and DMG concentrations were associated with a block of co-expressed genes enriched for adaptive immune, histone modification, and RNA processing functions. These results suggest that the prenatal maternal blood transcriptome is a sensitive indicator of gestational one-carbon metabolite status and changes relevant to children's later neurodevelopmental outcomes. LAY SUMMARY: Pregnancy is a time when maternal nutrition could interact with genetic risk for autism spectrum disorder. Blood samples collected during pregnancy from mothers who had a prior child with autism were examined for gene expression and nutrient metabolites, then compared to the diagnosis of the child at age three. Expression differences in gene pathways related to the immune system and gene regulation were observed for pregnancies of children with autism and non-typical neurodevelopment and were associated with maternal nutrients.

Cord blood DNA methylome in newborns later diagnosed with autism spectrum disorder reflects early dysregulation of neurodevelopmental and X-linked genes.

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.

mtDNA depletion-like syndrome in Wilson 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 CuSO4 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.

A meta-analysis of two high-risk prospective cohort studies reveals autism-specific transcriptional changes to chromatin, autoimmune, and environmental response genes in umbilical cord blood.

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 (log2(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.

Placental DNA methylation levels at CYP2E1 and IRS2 are associated with child outcome in a prospective autism study.

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.

Whole genome bisulfite sequencing of Down syndrome brain reveals regional DNA hypermethylation and novel disorder insights.

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.

Epigenomic signatures in liver and blood of Wilson disease patients include hypermethylation of liver-specific enhancers.

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.

Epigenetic changes of the thioredoxin system in the tx-j mouse model and in patients with Wilson disease.

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.

Snord116-dependent diurnal rhythm of DNA methylation in mouse cortex.

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.

Experience-dependent neuroplasticity of the developing hypothalamus: integrative epigenomic approaches.

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.

Design, synthesis and identification of a new class of triarylmethyl amine compounds as inhibitors of apolipoprotein E production.

We have identified a new class of triarylmethyl amine compounds that can inhibit apolipoprotein E (apoE) production. ApoE is a cholesterol- and lipid-carrier protein implicated in aging, atherosclerosis, Alzheimer's Disease (AD), and other neurological and lipid-related disorders. Attenuation of apoE production is generally considered to be of therapeutic value. A majority of the apoE in the brain is produced by astrocytes. Here, we describe the design, synthesis, and biological screening of a small library of compounds that led to the identification of four triarylmethyl amines as potent inhibitors of apoE production in CCF-STTG1 astrocytoma cells.