Ultra-Deep Bisulfite Sequencing to Detect Specific DNA Methylation Patterns of Minor Cell Types in Heterogeneous Cell Populations: An Example of the Pituitary Tissue

Epigenome editing revealed the role of DNA methylation of T-DMR/CpG island shore on Runx2 transcription

RUNX2 is a transcription factor crucial for bone formation. Mutant mice with varying levels of Runx2 expression display dosage-dependent skeletal abnormalities, underscoring the importance of Runx2 dosage control in skeletal formation. RUNX2 activity is regulated by several molecular mechanisms, including epigenetic modification such as DNA methylation. In this study, we investigated whether targeted repressive epigenome editing including hypermethylation to the Runx2-DMR/CpG island shore could influence Runx2 expression using Cas9-based epigenome-editing tools. Through the transient introduction of CRISPRoff-v2.1 and gRNAs targeting Runx2-DMR into MC3T3-E1 cells, we successfully induced hypermethylation of the region and concurrently reduced Runx2 expression during osteoblast differentiation. Although the epigenome editing of Runx2-DMR did not impact the expression of RUNX2 downstream target genes, these results indicate a causal relationship between the epigenetic status of the Runx2-DMR and Runx2 transcription. Additionally, we observed that hypermethylation of the Runx2-DMR persisted for at least 24 days under growth conditions but decreased during osteogenic differentiation, highlighting an endogenous DNA demethylation activity targeting the Runx2-DMR during the differentiation process. In summary, our study underscore the usefulness of the epigenome editing technology to evaluate the function of endogenous genetic elements and revealed that the Runx2-DMR methylation is actively regulated during osteoblast differentiation, subsequently could influence Runx2 expression.

Identification of cell type-specific methylation signals in bulk whole genome bisulfite sequencing data

BACKGROUND

The traditional approach to studying the epigenetic mechanism CpG methylation in tissue samples is to identify regions of concordant differential methylation spanning multiple CpG sites (differentially methylated regions). Variation limited to single or small numbers of CpGs has been assumed to reflect stochastic processes. To test this, we developed software, Cluster-Based analysis of CpG methylation (CluBCpG), and explored variation in read-level CpG methylation patterns in whole genome bisulfite sequencing data.

RESULTS

Analysis of both human and mouse whole genome bisulfite sequencing datasets reveals read-level signatures associated with cell type and cell type-specific biological processes. These signatures, which are mostly orthogonal to classical differentially methylated regions, are enriched at cell type-specific enhancers and allow estimation of proportional cell composition in synthetic mixtures and improved prediction of gene expression. In tandem, we developed a machine learning algorithm, Precise Read-Level Imputation of Methylation (PReLIM), to increase coverage of existing whole genome bisulfite sequencing datasets by imputing CpG methylation states on individual sequencing reads. PReLIM both improves CluBCpG coverage and performance and enables identification of novel differentially methylated regions, which we independently validate.

CONCLUSIONS

Our data indicate that, rather than stochastic variation, read-level CpG methylation patterns in tissue whole genome bisulfite sequencing libraries reflect cell type. Accordingly, these new computational tools should lead to an improved understanding of epigenetic regulation by DNA methylation.

DNA methylation ambiguity in the Fibrillin-1 (FBN1) CpG island shore possibly involved in Marfan syndrome

Fibrillin-1 (FBN1) is responsible for haploinsufficient and autosomal dominant Marfan syndrome. Even in the same Marfan pedigree, penetrance and expressivity in heterozygous individuals can differ and result in variable disease onset and severity. Thus, other factors in addition to mutations in FBN1 are likely to contribute to the disease. In this study, we examined the regulation of FBN1 in porcine Marfan syndrome model, focusing on DNA methylation patterns distinguishable as wild-type (WT) and FBN1 null (KO) alleles in heterozygous cells. Most importantly, the ratio of the transcriptionally active hypomethylated WT allele was altered during cellular passage and highly correlated with FBN1 mRNA level compared with that in the KO allele. Transcribed FBN1 RNA from the KO allele was abolished after splicing coupled with translational initiation, suggesting that the functional FBN1 mRNA levels were affected by DNA methylation of the WT allele.

DNA methylation in AgRP neurons regulates voluntary exercise behavior in mice

DNA methylation regulates cell type-specific gene expression. Here, in a transgenic mouse model, we show that deletion of the gene encoding DNA methyltransferase Dnmt3a in hypothalamic AgRP neurons causes a sedentary phenotype characterized by reduced voluntary exercise and increased adiposity. Whole-genome bisulfite sequencing (WGBS) and transcriptional profiling in neuronal nuclei from the arcuate nucleus of the hypothalamus (ARH) reveal differentially methylated genomic regions and reduced expression of AgRP neuron-associated genes in knockout mice. We use read-level analysis of WGBS data to infer putative ARH neural cell types affected by the knockout, and to localize promoter hypomethylation and increased expression of the growth factor Bmp7 to AgRP neurons, suggesting a role for aberrant TGF-β signaling in the development of this phenotype. Together, these data demonstrate that DNA methylation in AgRP neurons is required for their normal epigenetic development and neuron-specific gene expression profiles, and regulates voluntary exercise behavior.

AgRP neurons in the hypothalamic arcuate nucleus (ARH) are involved in regulating hunger and energy balance. Here the authors show that knockout of the DNA methyltransferase Dnmt3a in AgRP neurons of the ARH leads to a reduction in voluntary exercise along with numerous epigenetic and gene expression changes in ARH neurons.

A method for identifying allele-specific hydroxymethylation

We previously identified sequence-dependent allele-specific methylation (sd-ASM) in adult human peripheral blood leukocytes, in which ASM occurs in cis depending on adjacent polymorphic sequences. A number of groups have identified sd-ASM sites in the human and mouse genomes, illustrating the prevalence of sd-ASM in mammalian genomes. In addition, sd-ASM can lead to sequence-dependent allele-specific expression of neighbouring genes. Imprinted genes also often exhibit parent-of-origin-dependent allele-specific methylation (pd-ASM), which causes parent-of-origin-dependent allele-specific expression. However, whether most of the already known sd-ASM and pd-ASM sites are methylated or hydroxymethylated remains unclear due to technical restrictions. Accordingly, a novel method that enables examination of allelic methylation and hydroxymethylation status and also overcomes the drawbacks of conventional methods is needed. Such a method could also be used to elucidate the mechanisms underlying polymorphism-associated inter-individual differences in disease susceptibility and the mechanism of genomic imprinting. Here, we developed a simple method to determine allelic hydroxymethylation status and identified novel sequence- and parent-of-origin-dependent allele-specific hydroxymethylation sites. Correlation analyses of TF binding sequences and methylation or hydroxymethylation between three mouse strains revealed the involvement of Pax5 in strain-specific methylation and hydroxymethylation in exon 7 of Pdgfrb.

Establishment of DNA methylation patterns of the Fibrillin1 (FBN1) gene in porcine embryos and tissues

DNA methylation in transcriptional regulatory regions is crucial for gene expression. The DNA methylation status of the edges of CpG islands, called CpG island shore, is involved in tissue/cell-type-specific gene expression. Haploinsufficiency diseases are caused by inheritance of one mutated null allele and are classified as autosomal dominant. However, in the same pedigree, phenotypic variances are observed despite the inheritance of the identical mutated null allele, including Fibrillin1 (FBN1), which is responsible for development of the haploinsufficient Marfan disease. In this study, we examined the relationship between gene expression and DNA methylation patterns of the FBN1 CpG island shore focusing on transcriptionally active hypomethylated alleles (Hypo-alleles). No difference in the DNA methylation level of FBN1 CpG island shore was observed in porcine fetal fibroblast (PFF) and the liver, whereas FBN1 expression was higher in PFF than in the liver. However, Hypo-allele ratio of the FBN1 CpG island shore in PFF was higher than that in the liver, indicating that Hypo-allele ratio of the FBN1 CpG island shore likely correlated with FBN1 expression level. In addition, oocyte-derived DNA hypermethylation in preimplantation embryos was erased until the blastocyst stage, and re-methylation of the FBN1 CpG island shore was observed with prolonged in vitro culture of blastocysts. These results suggest that the establishment of the DNA methylation pattern within the FBN1 CpG island shore occurs after the blastocyst stage, likely during organogenesis. In conclusion, Hypo-allele ratios of the FBN1 CpG island shore correlated with FBN1 expression levels in porcine tissues.