Analysis of DNA methylation acquisition at the imprinted Dlk1 locus reveals asymmetry at CpG dyads

Maternal RNA transcription in Dlk1-Dio3 domain is critical for proper development of the mouse placental vasculature

The placenta is a unique organ for ensuring normal embryonic growth in the uterine. Here, we found that maternal RNA transcription in Dlk1-Dio3 imprinted domain is essential for placentation. PolyA signals were inserted into Gtl2 to establish a mouse model to prevent the expression of maternal RNAs in the domain. The maternal allele knock-in (MKI) and homozygous (HOMO) placentas showed an expanded junctional zone, reduced labyrinth and poor vasculature impacting both fetal and maternal blood spaces. The MKI and HOMO models displayed dysregulated gene expression in the Dlk1-Dio3 domain. In situ hybridization detected Dlk1, Gtl2, Rtl1, miR-127 and Rian dysregulated in the labyrinth vasculature. MKI and HOMO induced Dlk1 to lose imprinting, and DNA methylation changes of IG-DMR and Gtl2-DMR, leading to abnormal gene expression, while the above changes didn't occur in paternal allele knock-in placentas. These findings demonstrate that maternal RNAs in the Dlk1-Dio3 domain are involved in placental vasculature, regulating gene expression, imprinting status and DNA methylation.

Epigenetic control and genomic imprinting dynamics of the Dlk1-Dio3 domain

Genomic imprinting is an epigenetic process whereby genes are monoallelically expressed in a parent-of-origin-specific manner. Imprinted genes are frequently found clustered in the genome, likely illustrating their need for both shared regulatory control and functional inter-dependence. The Dlk1-Dio3 domain is one of the largest imprinted clusters. Genes in this region are involved in development, behavior, and postnatal metabolism: failure to correctly regulate the domain leads to Kagami-Ogata or Temple syndromes in humans. The region contains many of the hallmarks of other imprinted domains, such as long non-coding RNAs and parental origin-specific CTCF binding. Recent studies have shown that the Dlk1-Dio3 domain is exquisitely regulated via a bipartite imprinting control region (ICR) which functions differently on the two parental chromosomes to establish monoallelic expression. Furthermore, the Dlk1 gene displays a selective absence of imprinting in the neurogenic niche, illustrating the need for precise dosage modulation of this domain in different tissues. Here, we discuss the following: how differential epigenetic marks laid down in the gametes cause a cascade of events that leads to imprinting in the region, how this mechanism is selectively switched off in the neurogenic niche, and why studying this imprinted region has added a layer of sophistication to how we think about the hierarchical epigenetic control of genome function.

Methylation is maintained specifically at imprinting control regions but not other DMRs associated with imprinted genes in mice bearing a mutation in the Dnmt1 intrinsically disordered domain

Differential methylation of imprinting control regions in mammals is essential for distinguishing the parental alleles from each other and regulating their expression accordingly. To ensure parent of origin-specific expression of imprinted genes and thereby normal developmental progression, the differentially methylated states that are inherited at fertilization must be stably maintained by DNA methyltransferase 1 throughout subsequent somatic cell division. Further epigenetic modifications, such as the acquisition of secondary regions of differential methylation, are dependent on the methylation status of imprinting control regions and are important for achieving the monoallelic expression of imprinted genes, but little is known about how imprinting control regions direct the acquisition and maintenance of methylation at these secondary sites. Recent analysis has identified mutations that reduce DNA methyltransferase 1 fidelity at some genomic sequences but not at others, suggesting that it may function differently at different loci. We examined the impact of the mutant DNA methyltransferase 1 P allele on methylation at imprinting control regions as well as at secondary differentially methylated regions and non-imprinted sequences. We found that while the P allele results in a major reduction in DNA methylation levels across the mouse genome, methylation is specifically maintained at imprinting control regions but not at their corresponding secondary DMRs. This result suggests that DNA methyltransferase 1 may work differently at imprinting control regions or that there is an alternate mechanism for maintaining methylation at these critical regulatory regions and that maintenance of methylation at secondary DMRs is not solely dependent on the methylation status of the ICR.

Tissue-specific Grb10/Ddc insulator drives allelic architecture for cardiac development

Allele-specific expression of imprinted gene clusters is governed by gametic DNA methylation at master regulators called imprinting control regions (ICRs). Non-gametic or secondary differentially methylated regions (DMRs) at promoters and exonic regions reinforce monoallelic expression but do not control an entire cluster. Here, we unveil an unconventional secondary DMR that is indispensable for tissue-specific imprinting of two previously unlinked genes, Grb10 and Ddc. Using polymorphic mice, we mapped an intronic secondary DMR at Grb10 with paternal-specific CTCF binding (CBR2.3) that forms contacts with Ddc. Deletion of paternal CBR2.3 removed a critical insulator, resulting in substantial shifting of chromatin looping and ectopic enhancer-promoter contacts. Destabilized gene architecture precipitated abnormal Grb10-Ddc expression with developmental consequences in the heart and muscle. Thus, we redefine the Grb10-Ddc imprinting domain by uncovering an unconventional intronic secondary DMR that functions as an insulator to instruct the tissue-specific, monoallelic expression of multiple genes-a feature previously ICR exclusive.

Deletion of Meg8-DMR Enhances Migration and Invasion of MLTC-1 Depending on the CTCF Binding Sites

The Dlk1-Dio3 imprinted domain on mouse chromosome 12 contains three well-characterized paternally methylated differentially methylated regions (DMRs): IG-DMR, Gtl2-DMR, and Dlk1-DMR. These DMRs control the expression of many genes involved in embryonic development, inherited diseases, and human cancer in this domain. The first maternal methylation DMR discovered in this domain was the Meg8-DMR, the targets and biological function of which are still unknown. Here, using an enhancer-blocking assay, we first dissected the functional parts of the Meg8-DMR and showed that its insulator activity is dependent on the CCCTC-binding factor (CTCF) in MLTC-1. Results from RNA-seq showed that the deletion of the Meg8-DMR and its compartment CTCF binding sites, but not GGCG repeats, lead to the downregulation of numerous genes on chromosome 12, in particular the drastically reduced expression of Dlk1 and Rtl1 in the Dlk1-Dio3 domain, while differentially expressed genes are enriched in the MAPK pathway. In vitro assays revealed that the deletion of the Meg8-DMR and CTCF binding sites enhances cell migration and invasion by decreasing Dlk1 and activating the Notch1-Rhoc-MAPK/ERK pathway. These findings enhance research into gene regulation in the Dlk1-Dio3 domain by indicating that the Meg8-DMR functions as a long-range regulatory element which is dependent on CTCF binding sites and affects multiple genes in this domain.

Effects of paternal exposure to cigarette smoke on sperm DNA methylation and long-term metabolic syndrome in offspring

BACKGROUND

Although paternal exposure to cigarette smoke may contribute to obesity and metabolic syndrome in offspring, the underlying mechanisms remain uncertain.

METHODS

In the present study, we analyzed the sperm DNA-methylation profiles in tobacco-smoking normozoospermic (SN) men, non-tobacco-smoking normozoospermic (N) men, and non-smoking oligoasthenozoospermic (OA) men. Using a mouse model, we also analyzed global methylation and differentially methylated regions (DMRs) of the DLK1 gene in paternal spermatozoa and the livers of progeny. In addition, we quantified DLK1 expression, executed an intra-peritoneal glucose tolerance test (IPGTT), measured serum metabolites, and analyzed liver lipid accumulation in the F1 offspring.

RESULTS

Global sperm DNA-methylation levels were significantly elevated (p  < 0.05) in the SN group, and the methylation patterns were different among N, SN, and OA groups. Importantly, the methylation level of the DLK1 locus (cg11193865) was significantly elevated in the SN group compared to both N and OA groups (p  < 0.001). In the mouse model, the group exposed to cigarette smoke extract (CSE) exhibited a significantly higher global methylation DNA level in spermatozoa (p  < 0.001) and on the DMR sites of Dlk1 in 10-week-old male offspring (p  < 0.05), with a significant increase in Dlk1 expression in their livers (p  < 0.001). In addition, IPGTT and LDL levels were significantly altered (p  < 0.001), with elevated liver fat accumulation (p  < 0.05) in F1 offspring.

CONCLUSION

Paternal exposure to cigarette smoke led to increased global methylation of sperm DNA and alterations to the DMR of the DLK1 gene in the F1 generation, which may be inherited parentally and may perturb long-term metabolic function.

TRIM28 Regulates Dlk1 Expression in Adipogenesis

The tripartite motif-containing protein 28 (TRIM28) is a transcription corepressor, interacting with histone deacetylase and methyltransferase complexes. TRIM28 is a crucial regulator in development and differentiation. We would like to investigate its function and regulation in adipogenesis. Knockdown of Trim28 by transducing lentivirus-carrying shRNAs impairs the differentiation of 3T3-L1 preadipocytes, demonstrated by morphological observation and gene expression analysis. To understand the molecular mechanism of Trim28-mediated adipogenesis, the RNA-seq was performed to find out the possible Trim28-regulated genes. Dlk1 (delta-like homolog 1) was increased in Trim28 knockdown 3T3-L1 cells both untreated and induced to differentiation. Dlk1 is an imprinted gene and known as an inhibitor of adipogenesis. Further knockdown of Dlk1 in Trim28 knockdown 3T3-L1 would rescue cell differentiation. The epigenetic analysis showed that DNA methylation of Dlk1 promoter and differentially methylated regions (DMRs) was not altered significantly in Trim28 knockdown cells. However, compared to control cells, the histone methylation on the Dlk1 promoter was increased at H3K4 and decreased at H3K27 in Trim28 knockdown cells. Finally, we found Trim28 might be recruited by transcription factor E2f1 to regulate Dlk1 expression. The results imply Trim28-Dlk1 axis is critical for adipogenesis.

Loss of TSC complex enhances gluconeogenesis via upregulation of Dlk1-Dio3 locus miRNAs

Loss of the tumor suppressor tuberous sclerosis complex 1 (Tsc1) in the liver promotes gluconeogenesis and glucose intolerance. We asked whether this could be attributed to aberrant expression of small RNAs. We performed small-RNA sequencing on liver of Tsc1-knockout mice, and found that miRNAs of the delta-like homolog 1 (Dlk1)-deiodinase iodothyronine type III (Dio3) locus are up-regulated in an mTORC1-dependent manner. Sustained mTORC1 signaling during development prevented CpG methylation and silencing of the Dlk1-Dio3 locus, thereby increasing miRNA transcription. Deletion of miRNAs encoded by the Dlk1-Dio3 locus reduced gluconeogenesis, glucose intolerance, and fasting blood glucose levels. Thus, miRNAs contribute to the metabolic effects observed upon loss of TSC1 and hyperactivation of mTORC1 in the liver. Furthermore, we show that miRNA is a downstream effector of hyperactive mTORC1 signaling.

Hemimethylation of CpG dyads is characteristic of secondary DMRs associated with imprinted loci and correlates with 5-hydroxymethylcytosine at paternally methylated sequences

BACKGROUND

In mammals, the regulation of imprinted genes is controlled by differential methylation at imprinting control regions which acquire parent of origin-specific methylation patterns during gametogenesis and retain differences in allelic methylation status throughout fertilization and subsequent somatic cell divisions. In addition, many imprinted genes acquire differential methylation during post-implantation development; these secondary differentially methylated regions appear necessary to maintain the imprinted expression state of individual genes. Despite the requirement for both types of differentially methylated sequence elements to achieve proper expression across imprinting clusters, methylation patterns are more labile at secondary differentially methylated regions. To understand the nature of this variability, we analyzed CpG dyad methylation patterns at both paternally and maternally methylated imprinted loci within multiple imprinting clusters.

RESULTS

We determined that both paternally and maternally methylated secondary differentially methylated regions associated with imprinted genes display high levels of hemimethylation, 29-49%, in comparison to imprinting control regions which exhibited 8-12% hemimethylation. To explore how hemimethylation could arise, we assessed the differentially methylated regions for the presence of 5-hydroxymethylcytosine which could cause methylation to be lost via either passive and/or active demethylation mechanisms. We found enrichment of 5-hydroxymethylcytosine at paternally methylated secondary differentially methylated regions, but not at the maternally methylated sites we analyzed in this study.

CONCLUSIONS

We found high levels of hemimethylation to be a generalizable characteristic of secondary differentially methylated regions associated with imprinted genes. We propose that 5-hydroxymethylcytosine enrichment may be responsible for the variability in methylation status at paternally methylated secondary differentially methylated regions associated with imprinted genes. We further suggest that the high incidence of hemimethylation at secondary differentially methylated regions must be counteracted by continuous methylation acquisition at these loci.

Transit amplifying cells coordinate mouse incisor mesenchymal stem cell activation

Stem cells (SCs) receive inductive cues from the surrounding microenvironment and cells. Limited molecular evidence has connected tissue-specific mesenchymal stem cells (MSCs) with mesenchymal transit amplifying cells (MTACs). Using mouse incisor as the model, we discover a population of MSCs neibouring to the MTACs and epithelial SCs. With Notch signaling as the key regulator, we disclose molecular proof and lineage tracing evidence showing the distinct MSCs contribute to incisor MTACs and the other mesenchymal cell lineages. MTACs can feedback and regulate the homeostasis and activation of CL-MSCs through Delta-like 1 homolog (Dlk1), which balances MSCs-MTACs number and the lineage differentiation. Dlk1's function on SCs priming and self-renewal depends on its biological forms and its gene expression is under dynamic epigenetic control. Our findings can be validated in clinical samples and applied to accelerate tooth wound healing, providing an intriguing insight of how to direct SCs towards tissue regeneration.

Overexpression of microRNAs from the Gtl2-Rian locus contributes to postnatal death in mice

The Dlk1-Dio3 imprinted domain functions in embryonic development but the roles of noncoding RNAs expressed from this domain remain unclear. We addressed this question by generating transgenic (TG) mice harbouring a BAC carrying IG-DMR (intergenic-differentially methylated region), Gtl2-DMR, Gtl2, Rtl1/Rtl1as, and part of Rian. High postnatal lethality (>85%) of the BAC-TG pups was observed in the maternally transmitted individuals (MAT-TG), but not following paternal transmission (PAT-TG). The DNA methylation status of IG-DMR and Gtl2-DMR in the BAC-allele was paternally imprinted similar to the genomic allele. The mRNA-Seq and miRNA-Seq analysis revealed marked expression changes in the MAT-TG, with 1,500 upregulated and 2,131 downregulated genes. The long noncoding RNAs and 12 miRNAs containing the BAC locus were markedly enhanced in the MAT-TG. We identified the 24 target genes of the overexpressed miRNAs and confirmed the downregulation in the MAT-TG. Notably, overexpression of mir770, mir493, and mir665 from Gtl2 in the MAT-TG embryos led to decreased expression of the 3 target genes, Col5a1, Pcgf2, and Clip2. Our results suggest that decreased expression of the 3 target genes concomitant with overexpression of the miRNAs within Gtl2 may be involved in the postnatal death in the MAT-TG. Because this imprinted domain is well conserved between mice and humans, the results of genetic and molecular analysis in mice hold important implications for related human disorders such as Temple syndrome.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Detecting rare asymmetrically methylated cytosines and decoding methylation patterns in the honeybee genome

Context-dependent gene expression in eukaryotes is controlled by several mechanisms including cytosine methylation that primarily occurs in the CG dinucleotides (CpGs). However, less frequent non-CpG asymmetric methylation has been found in various cell types, such as mammalian neurons, and recent results suggest that these sites can repress transcription independently of CpG contexts. In addition, an emerging view is that CpG hemimethylation may arise not only from deregulation of cellular processes but also be a standard feature of the methylome. Here, we have applied a novel approach to examine whether asymmetric CpG methylation is present in a sparsely methylated genome of the honeybee, a social insect with a high level of epigenetically driven phenotypic plasticity. By combining strand-specific ultra-deep amplicon sequencing of illustrator genes with whole-genome methylomics and bioinformatics, we show that rare asymmetrically methylated CpGs can be unambiguously detected in the honeybee genome. Additionally, we confirm differential methylation between two phenotypically and reproductively distinct castes, queens and workers, and offer new insight into the heterogeneity of brain methylation patterns. In particular, we challenge the assumption that symmetrical methylation levels reflect symmetry in the underlying methylation patterns and conclude that hemimethylation may occur more frequently than indicated by methylation levels. Finally, we question the validity of a prior study in which most of cytosine methylation in this species was reported to be asymmetric.

Asymmetric DNA methylation of CpG dyads is a feature of secondary DMRs associated with the Dlk1/Gtl2 imprinting cluster in mouse

Background

Differential DNA methylation plays a critical role in the regulation of imprinted genes. The differentially methylated state of the imprinting control region is inherited via the gametes at fertilization, and is stably maintained in somatic cells throughout development, influencing the expression of genes across the imprinting cluster. In contrast, DNA methylation patterns are more labile at secondary differentially methylated regions which are established at imprinted loci during post-implantation development. To investigate the nature of these more variably methylated secondary differentially methylated regions, we adopted a hairpin linker bisulfite mutagenesis approach to examine CpG dyad methylation at differentially methylated regions associated with the murine Dlk1/Gtl2 imprinting cluster on both complementary strands.

Results

We observed homomethylation at greater than 90% of the methylated CpG dyads at the IG-DMR, which serves as the imprinting control element. In contrast, homomethylation was only observed at 67–78% of the methylated CpG dyads at the secondary differentially methylated regions; the remaining 22–33% of methylated CpG dyads exhibited hemimethylation.

Conclusions

We propose that this high degree of hemimethylation could explain the variability in DNA methylation patterns at secondary differentially methylated regions associated with imprinted loci. We further suggest that the presence of 5-hydroxymethylation at secondary differentially methylated regions may result in hemimethylation and methylation variability as a result of passive and/or active demethylation mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0138-0) contains supplementary material, which is available to authorized users.