Haplotype-dependent allele-specific methylation (hap-ASM) can impact disease susceptibility, but maps of this phenomenon using stringent criteria in disease-relevant tissues remain sparse. Here we apply array-based and Methyl-Seq approaches to multiple human tissues and cell types, including brain, purified neurons and glia, T lymphocytes, and placenta, and identify 795 hap-ASM differentially methylated regions (DMRs) and 3,082 strong methylation quantitative trait loci (mQTLs), most not previously reported. More than half of these DMRs have cell type-restricted ASM, and among them are 188 hap-ASM DMRs and 933 mQTLs located near GWAS signals for immune and neurological disorders. Targeted bis-seq confirmed hap-ASM in 12/13 loci tested, including CCDC155, CD69, FRMD1, IRF1, KBTBD11, and S100A(∗)-ILF2, associated with immune phenotypes, MYT1L, PTPRN2, CMTM8 and CELF2, associated with neurological disorders, NGFR and HLA-DRB6, associated with both immunological and brain disorders, and ZFP57, a trans-acting regulator of genomic imprinting. Polymorphic CTCF and transcription factor (TF) binding sites were over-represented among hap-ASM DMRs and mQTLs, and analysis of the human data, supplemented by cross-species comparisons to macaques, indicated that CTCF and TF binding likelihood predicts the strength and direction of the allelic methylation asymmetry. These results show that hap-ASM is highly tissue specific; an important trans-acting regulator of genomic imprinting is regulated by this phenomenon; and variation in CTCF and TF binding sites is an underlying mechanism, and maps of hap-ASM and mQTLs reveal regulatory sequences underlying supra- and sub-threshold GWAS peaks in immunological and neurological disorders.

Endometriosis is characterized by growth of endometrial-like tissue outside the uterine cavity. Since its pathogenesis may involve epigenetic changes, we used Illumina 450K Methylation Beadchips to profile CpG methylation in endometriosis stromal cells compared to stromal cells from normal endometrium. We validated and extended the Beadchip data using bisulfite sequencing (bis-seq), and analyzed differential methylation (DM) at the CpG-level and by an element-level classification for groups of CpGs in chromatin domains. Genes found to have DM included examples encoding transporters (SLC22A23), signaling components (BDNF, DAPK1, ROR1, and WNT5A) and transcription factors (GATA family, HAND2, HOXA cluster, NR5A1, OSR2, TBX3). Intriguingly, among the TF genes with DM we also found JAZF1, a proto-oncogene affected by chromosomal translocations in endometrial stromal tumors. Using RNA-Seq we identified a subset of the DM genes showing differential expression (DE), with the likelihood of DE increasing with the extent of the DM and its location in enhancer elements. Supporting functional relevance, treatment of stromal cells with the hypomethylating drug 5aza-dC led to activation of DAPK1 and SLC22A23 and repression of HAND2, JAZF1, OSR2, and ROR1 mRNA expression. We found that global 5hmC is decreased in endometriotic versus normal epithelial but not stroma cells, and for JAZF1 and BDNF examined by oxidative bis-seq, found that when 5hmC is detected, patterns of 5hmC paralleled those of 5mC. Together with prior studies, these results define a consistent epigenetic signature in endometriosis stromal cells and nominate specific transcriptional and signaling pathways as therapeutic targets.

No abstract