The incidence of type 1 diabetes (T1D) has substantially increased over the past decade, suggesting a role for non-genetic factors such as epigenetic mechanisms in disease development. Here we present an epigenome-wide association study across 406,365 CpGs in 52 monozygotic twin pairs discordant for T1D in three immune effector cell types. We observe a substantial enrichment of differentially variable CpG positions (DVPs) in T1D twins when compared with their healthy co-twins and when compared with healthy, unrelated individuals. These T1D-associated DVPs are found to be temporally stable and enriched at gene regulatory elements. Integration with cell type-specific gene regulatory circuits highlight pathways involved in immune cell metabolism and the cell cycle, including mTOR signalling. Evidence from cord blood of newborns who progress to overt T1D suggests that the DVPs likely emerge after birth. Our findings, based on 772 methylomes, implicate epigenetic changes that could contribute to disease pathogenesis in T1D.

BAHD1 is a vertebrate protein that promotes heterochromatin formation and gene repression in association with several epigenetic regulators. However, its physiological roles remain unknown. Here, we demonstrate that ablation of the Bahd1 gene results in hypocholesterolemia, hypoglycemia and decreased body fat in mice. It also causes placental growth restriction with a drop of trophoblast glycogen cells, a reduction of fetal weight and a high neonatal mortality rate. By intersecting transcriptome data from murine Bahd1 knockout (KO) placentas at stages E16.5 and E18.5 of gestation, Bahd1-KO embryonic fibroblasts, and human cells stably expressing BAHD1, we also show that changes in BAHD1 levels alter expression of steroid/lipid metabolism genes. Biochemical analysis of the BAHD1-associated multiprotein complex identifies MIER proteins as novel partners of BAHD1 and suggests that BAHD1-MIER interaction forms a hub for histone deacetylases and methyltransferases, chromatin readers and transcription factors. We further show that overexpression of BAHD1 leads to an increase of MIER1 enrichment on the inactive X chromosome (Xi). In addition, BAHD1 and MIER1/3 repress expression of the steroid hormone receptor genes ESR1 and PGR, both playing important roles in placental development and energy metabolism. Moreover, modulation of BAHD1 expression in HEK293 cells triggers epigenetic changes at the ESR1 locus. Together, these results identify BAHD1 as a core component of a chromatin-repressive complex regulating placental morphogenesis and body fat storage and suggest that its dysfunction may contribute to several human diseases.

Phylogenetic analysis of the endo-beta-1,4-glucanase gene family of Arabidopsis and other plants revealed a clear distinction in three subfamilies (alpha, beta, and gamma). The alpha- and beta-subfamily contains proteins believed to be involved in a number of physiological roles such as elongation, ripening, and abscission. The gamma-subfamily is composed of proteins that are predicted to have a membrane-spanning domain and to be localized at the plasma membrane. Some of these proteins have been linked to cellulose biosynthesis by serving to hydrolyze a lipid-linked intermediate that acts as a primer for the elongation of beta-glucan chains during cellulose synthesis at the plasma membrane. Similar glucanases are important in cellulose biosynthesis in bacteria. Searches in the genomes of unrelated organisms that make cellulose, such as Ciona intestinalis and Dictyostelium discoideum, revealed the presence of membrane-linked endo-beta-1,4-glucanases and it is suggested that these might also have a role in cellulose synthesis.

BAH domain-containing protein 1 (BAHD1) is involved in heterochromatin formation and gene repression in human cells. BAHD1 also localizes to the inactive X chromosome (Xi), but the functional significance of this targeting is unknown. So far, research on this protein has been hampered by its low endogenous abundance and its role in epigenetic regulation remains poorly explored. In this work, we used whole-genome bisulfite sequencing (BS-seq) to compare the DNA methylation profile of HEK293 cells expressing low levels of BAHD1 (HEK-CT) to that of isogenic cells stably overexpressing BAHD1 (HEK-BAHD1). We show that increasing BAHD1 levels induces de novo DNA methylation on autosomes and a marked hypomethylation on the X chromosome (chrX). We identified 91,358 regions that have different methylation patterns in HEK-BAHD1 compared to HEK-CT cells (termed “BAHD1-DMRs”), of which 83,850 mapped on autosomes and 7508 on the X chromosome (chrX). Autosomal BAHD1-DMRs were predominantly hypermethylated and located to satellites, interspersed repeats, and intergenic regions. In contrast, BAHD1-DMRs on chrX were mainly hypomethylated and located to gene bodies and enhancers. We further found that BAHD1-DMRs display a higher-order organization by being clustered within large chromosomal domains. Half of these “BAHD1-Associated differentially methylated Domains” (BADs) overlapped with lamina-associated domains (LADs). Based on these results, we propose that BAHD1-mediated heterochromatin formation is linked to DNA methylation and may play a role in the spatial architecture of the genome.

The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.