Tissue-specific epigenetics in gene neighborhoods: myogenic transcription factor genes

Myogenic regulatory factor (MRF) genes, MYOD1, MYOG, MYF6 and MYF5, are critical for the skeletal muscle lineage. Here, we used various epigenome profiles from human myoblasts (Mb), myotubes (Mt), muscle and diverse non-muscle samples to elucidate the involvement of multigene neighborhoods in the regulation of MRF genes. We found more far-distal enhancer chromatin associated with MRF genes in Mb and Mt than previously reported from studies in mice. For the MYF5/MYF6 gene-pair, regions of Mb-associated enhancer chromatin were located throughout the adjacent 236-kb PTPRQ gene even though Mb expressed negligible amounts of PTPRQ mRNA. Some enhancer chromatin regions inside PTPRQ in Mb were also seen in PTPRQ mRNA-expressing non-myogenic cells. This suggests dual-purpose PTPRQ enhancers that upregulate expression of PTPRQ in non-myogenic cells and MYF5/MYF6 in myogenic cells. In contrast, the myogenic enhancer chromatin regions distal to MYOD1 were intergenic and up to 19 kb long. Two of them contain small, known MYOD1 enhancers, and one displayed an unusually high level of 5-hydroxymethylcytosine in a quantitative DNA hydroxymethylation assay. Unexpectedly, three regions of MYOD1-distal enhancer chromatin in Mb and Mt overlapped enhancer chromatin in umbilical vein endothelial cells, which might upregulate a distant gene (PIK3C2A). Lastly, genes surrounding MYOG were preferentially transcribed in Mt, like MYOG itself, and exhibited nearby myogenic enhancer chromatin. These neighboring chromatin regions may be enhancers acting in concert to regulate myogenic expression of multiple adjacent genes. Our findings reveal the very different and complex organization of gene neighborhoods containing closely related transcription factor genes.

Small RNA-mediated DNA (cytosine-5) methyltransferase 1 inhibition leads to aberrant DNA methylation

Mammalian cells contain copious amounts of RNA including both coding and noncoding RNA (ncRNA). Generally the ncRNAs function to regulate gene expression at the transcriptional and post-transcriptional level. Among ncRNA, the long ncRNA and small ncRNA can affect histone modification, DNA methylation targeting and gene silencing. Here we show that endogenous DNA methyltransferase 1 (DNMT1) co-purifies with inhibitory ncRNAs. MicroRNAs (miRNAs) bind directly to DNMT1 with high affinity. The binding of miRNAs, such as miR-155-5p, leads to inhibition of DNMT1 enzyme activity. Exogenous miR-155-5p in cells induces aberrant DNA methylation of the genome, resulting in hypomethylation of low to moderately methylated regions. And small shift of hypermethylation of previously hypomethylated region was also observed. Furthermore, hypomethylation led to activation of genes. Based on these observations, overexpression of miR-155-5p resulted in aberrant DNA methylation by inhibiting DNMT1 activity, resulting in altered gene expression.

Notch signaling genes: myogenic DNA hypomethylation and 5-hydroxymethylcytosine

Notch intercellular signaling is critical for diverse developmental pathways and for homeostasis in various types of stem cells and progenitor cells. Because Notch gene products need to be precisely regulated spatially and temporally, epigenetics is likely to help control expression of Notch signaling genes. Reduced representation bisulfite sequencing (RRBS) indicated significant hypomethylation in myoblasts, myotubes, and skeletal muscle vs. many nonmuscle samples at intragenic or intergenic regions of the following Notch receptor or ligand genes: NOTCH1, NOTCH2, JAG2, and DLL1. An enzymatic assay of sites in or near these genes revealed unusually high enrichment of 5-hydroxymethylcytosine (up to 81%) in skeletal muscle, heart, and cerebellum. Epigenetics studies and gene expression profiles suggest that hypomethylation and/or hydroxymethylation help control expression of these genes in heart, brain, myoblasts, myotubes, and within skeletal muscle myofibers. Such regulation could promote cell renewal, cell maintenance, homeostasis, and a poised state for repair of tissue damage.

Methyllysine reader plant homeodomain (PHD) finger protein 20-like 1 (PHF20L1) antagonizes DNA (cytosine-5) methyltransferase 1 (DNMT1) proteasomal degradation

Inheritance of DNA cytosine methylation pattern during successive cell division is mediated by maintenance DNA (cytosine-5) methyltransferase 1 (DNMT1). Lysine 142 of DNMT1 is methylated by the SET domain containing lysine methyltransferase 7 (SET7), leading to its degradation by proteasome. Here we show that PHD finger protein 20-like 1 (PHF20L1) regulates DNMT1 turnover in mammalian cells. Malignant brain tumor (MBT) domain of PHF20L1 binds to monomethylated lysine 142 on DNMT1 (DNMT1K142me1) and colocalizes at the perinucleolar space in a SET7-dependent manner. PHF20L1 knockdown by siRNA resulted in decreased amounts of DNMT1 on chromatin. Ubiquitination of DNMT1K142me1 was abolished by overexpression of PHF20L1, suggesting that its binding may block proteasomal degradation of DNMT1K142me1. Conversely, siRNA-mediated knockdown of PHF20L1 or incubation of a small molecule MBT domain binding inhibitor in cultured cells accelerated the proteasomal degradation of DNMT1. These results demonstrate that the MBT domain of PHF20L1 reads and controls enzyme levels of methylated DNMT1 in cells, thus representing a novel antagonist of DNMT1 degradation.

DNA methylation and differentiation: HOX genes in muscle cells

BACKGROUND

Tight regulation of homeobox genes is essential for vertebrate development. In a study of genome-wide differential methylation, we recently found that homeobox genes, including those in the HOX gene clusters, were highly overrepresented among the genes with hypermethylation in the skeletal muscle lineage. Methylation was analyzed by reduced representation bisulfite sequencing (RRBS) of postnatal myoblasts, myotubes and adult skeletal muscle tissue and 30 types of non-muscle-cell cultures or tissues.

RESULTS

In this study, we found that myogenic hypermethylation was present in specific subregions of all four HOX gene clusters and was associated with various chromatin epigenetic features. Although the 3' half of the HOXD cluster was silenced and enriched in polycomb repression-associated H3 lysine 27 trimethylation in most examined cell types, including myoblasts and myotubes, myogenic samples were unusual in also displaying much DNA methylation in this region. In contrast, both HOXA and HOXC clusters displayed myogenic hypermethylation bordering a central region containing many genes preferentially expressed in myogenic progenitor cells and consisting largely of chromatin with modifications typical of promoters and enhancers in these cells. A particularly interesting example of myogenic hypermethylation was HOTAIR, a HOXC noncoding RNA gene, which can silence HOXD genes in trans via recruitment of polycomb proteins. In myogenic progenitor cells, the preferential expression of HOTAIR was associated with hypermethylation immediately downstream of the gene. Other HOX gene regions also displayed myogenic DNA hypermethylation despite being moderately expressed in myogenic cells. Analysis of representative myogenic hypermethylated sites for 5-hydroxymethylcytosine revealed little or none of this base, except for an intragenic site in HOXB5 which was specifically enriched in this base in skeletal muscle tissue, whereas myoblasts had predominantly 5-methylcytosine at the same CpG site.

CONCLUSIONS

Our results suggest that myogenic hypermethylation of HOX genes helps fine-tune HOX sense and antisense gene expression through effects on 5' promoters, intragenic and intergenic enhancers and internal promoters. Myogenic hypermethylation might also affect the relative abundance of different RNA isoforms, facilitate transcription termination, help stop the spread of activation-associated chromatin domains and stabilize repressive chromatin structures.

Early de novo DNA methylation and prolonged demethylation in the muscle lineage

Myogenic cell cultures derived from muscle biopsies are excellent models for human cell differentiation. We report the first comprehensive analysis of myogenesis-specific DNA hyper- and hypo-methylation throughout the genome for human muscle progenitor cells (both myoblasts and myotubes) and skeletal muscle tissue vs. 30 non-muscle samples using reduced representation bisulfite sequencing. We also focused on four genes with extensive hyper- or hypo-methylation in the muscle lineage (PAX3, TBX1, MYH7B/MIR499 and OBSCN) to compare DNA methylation, DNaseI hypersensitivity, histone modification, and CTCF binding profiles. We found that myogenic hypermethylation was strongly associated with homeobox or T-box genes and muscle hypomethylation with contractile fiber genes. Nonetheless, there was no simple relationship between differential gene expression and myogenic differential methylation, rather only for subsets of these genes, such as contractile fiber genes. Skeletal muscle retained ~30% of the hypomethylated sites but only ~3% of hypermethylated sites seen in myogenic progenitor cells. By enzymatic assays, skeletal muscle was 2-fold enriched globally in genomic 5-hydroxymethylcytosine (5-hmC) vs. myoblasts or myotubes and was the only sample type enriched in 5-hmC at tested myogenic hypermethylated sites in PAX3/CCDC140 andTBX1. TET1 and TET2 RNAs, which are involved in generation of 5-hmC and DNA demethylation, were strongly upregulated in myoblasts and myotubes. Our findings implicate de novo methylation predominantly before the myoblast stage and demethylation before and after the myotube stage in control of transcription and co-transcriptional RNA processing. They also suggest that, in muscle, TET1 or TET2 are involved in active demethylation and in formation of stable 5-hmC residues.

Dual binding of chromomethylase domains to H3K9me2-containing nucleosomes directs DNA methylation in plants

DNA methylation and histone modification exert epigenetic control over gene expression. CHG methylation by CHROMOMETHYLASE3 (CMT3) depends on histone H3K9 dimethylation (H3K9me2), but the mechanism underlying this relationship is poorly understood. Here, we report multiple lines of evidence that CMT3 interacts with H3K9me2-containing nucleosomes. CMT3 genome locations nearly perfectly correlated with H3K9me2, and CMT3 stably associated with H3K9me2-containing nucleosomes. Crystal structures of maize CMT3 homolog ZMET2, in complex with H3K9me2 peptides, showed that ZMET2 binds H3K9me2 via both bromo adjacent homology (BAH) and chromo domains. The structures reveal an aromatic cage within both BAH and chromo domains as interaction interfaces that capture H3K9me2. Mutations that abolish either interaction disrupt CMT3 binding to nucleosomes and show a complete loss of CMT3 activity in vivo. Our study establishes dual recognition of H3K9me2 marks by BAH and chromo domains and reveals a distinct mechanism of interplay between DNA methylation and histone modification.

Abstract LB-381: Biochemical characterization of recombinant β-glucosyltransferase and analysis of global 5-hydroxymethylcytosine in unique genomes

5-hydroxymethylcytosine (5-hmC) is an enzymatic oxidative product of 5-methylcytosine (5-mC). Ten Eleven Translocation (TET) family of enzymes catalyze the conversion of 5-mC to 5-hmC. Phage encoded glucosyltransferases are known to glucosylate 5-hmC, which can be utilized to detect and analyze the 5-hmC an epigenetic mark in the mammalian epigenome. Here we have performed a detailed biochemical characterization and steady-state kinetic parameter analysis of T4 phage β-glucosyltransferase (β-GT). Recombinant β-GT glucosylates 5-hmC DNA in a non-processive manner and binding to either 5-hmC DNA or uridine diphosphoglucose (UDP-glucose) substrates is random, with both binary complexes being catalytically competent. Product inhibition studies with β-GT demonstrated that UDP is a competitive inhibitor to UDP-glucose and a mixed inhibitor to 5-hmC DNA. Similarly, the glucosylated-5-hmC (5-ghmC) DNA is a competitive inhibitor to 5-hmC DNA and mixed inhibitor to UDP-glucose. 5-hmC DNA binds ∼10 fold stronger to the β-GT enzyme when compared to its glucosylated product. The numbers of 5-hmC on target sequences influenced the turnover numbers for recombinant β-GT. Furthermore, we have utilized recombinant β-GT to estimate global 5-hmC content in a variety of genomic DNAs. Most of the genomic DNAs derived from vertebrate tissue/cell-lines contained 5-hmC. DNA from mouse, human and bovine brains displayed 0.5 to 0.9% of the total nucleotides as 5-hmC, which was higher compared to other tissues. A comparison between cancer and healthy tissue genomes suggested a lower percentage of 5-hmC in cancer, which may reflect the global hypomethylation of 5-mC observed during oncogenesis.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr LB-381. doi:1538-7445.AM2012-LB-381

Biochemical characterization of recombinant β-glucosyltransferase and analysis of global 5-hydroxymethylcytosine in unique genomes

5-Hydroxymethylcytosine (5-hmC) is an enzymatic oxidative product of 5-methylcytosine (5-mC). The Ten Eleven Translocation (TET) family of enzymes catalyze the conversion of 5-mC to 5-hmC. Phage-encoded glucosyltransferases are known to glucosylate 5-hmC, which can be utilized to detect and analyze the 5-hmC as an epigenetic mark in the mammalian epigenome. Here we have performed a detailed biochemical characterization and steady-state kinetic parameter analysis of T4 phage β-glucosyltransferase (β-GT). Recombinant β-GT glucosylates 5-hmC DNA in a nonprocessive manner, and binding to either 5-hmC DNA or uridine diphosphoglucose (UDP-glucose) substrates is random, with both binary complexes being catalytically competent. Product inhibition studies with β-GT demonstrated that UDP is a competitive inhibitor with respect to UDP-glucose and a mixed inhibitor with respect to 5-hmC DNA. Similarly, the glucosylated-5-hmC (5-ghmC) DNA is a competitive inhibitor with respect to 5-hmC DNA and mixed inhibitor with respect to UDP-glucose. 5-hmC DNA binds ∼10 fold stronger to the β-GT enzyme when compared to its glucosylated product. The numbers of 5-hmC on target sequences influenced the turnover numbers for recombinant β-GT. Furthermore, we have utilized recombinant β-GT to estimate global 5-hmC content in a variety of genomic DNAs. Most of the genomic DNAs derived from vertebrate tissue and cell lines contained 5-hmC. DNA from mouse, human, and bovine brains displayed 0.5–0.9% of the total nucleotides as 5-hmC, which was higher compared to the levels found in other tissues. A comparison between cancer and healthy tissue genomes suggested a lower percentage of 5-hmC in cancer, which may reflect the global hypomethylation of 5-mC observed during oncogenesis.

The E-box binding factors Max/Mnt, MITF, and USF1 act coordinately with FoxO to regulate expression of proapoptotic and cell cycle control genes by phosphatidylinositol 3-kinase/Akt/glycogen synthase kinase 3 signaling

Phosphatidylinositol (PI) 3-kinase/Akt signaling plays a critical role in cell proliferation and survival, partly by regulation of FoxO transcription factors. Previous work using global expression profiling indicated that inhibition of PI 3-kinase in proliferating cells led to induction of genes that promote cell cycle arrest and apoptosis. The upstream regulatory regions of these genes had binding sites not only for FoxO, but also for Myc/Max transcription factors. In the present study, we have addressed the role of Myc family members and related E-box-binding proteins in the regulation of these genes. Chromatin immunoprecipitations and RNA interference indicated that transcription was repressed by Max-Mnt-Sin3a-histone deacetylase complexes in proliferating cells. Inhibition of PI 3-kinase led to a loss of Max/Mnt binding and transcriptional induction by MITF and USF1, as well as FoxO. Both MITF and USF1 were activated by glycogen synthase kinase (GSK) 3, with GSK3 phosphorylation sites on USF1 identified as the previously described activating site threonine 153 as well as serine 186. siRNA against MITF as well as against FoxO3a protected cells from apoptosis following PI 3-kinase inhibition. These results define a novel E-box-regulated network that functions coordinately with FoxO to regulate transcription of apoptotic and cell cycle regulatory genes downstream of PI 3-kinase/Akt/GSK3 signaling.

Phosphatidylinositol 3-kinase signaling in proliferating cells maintains an anti-apoptotic transcriptional program mediated by inhibition of FOXO and non-canonical activation of NFkappaB transcription factors

BACKGROUND

Phosphatidylinositol (PI) 3-kinase is activated by a variety of growth factor receptors and the PI 3-kinase/Akt signaling pathway is a key regulator of cell proliferation and survival. The downstream targets of PI 3-kinase/Akt signaling include direct regulators of cell cycle progression and apoptosis as well as a number of transcription factors. Growth factor stimulation of quiescent cells leads to robust activation of PI 3-kinase, induction of immediate-early genes, and re-entry into the cell cycle. A lower level of PI 3-kinase signaling is also required for the proliferation and survival of cells maintained in the presence of growth factors, but the gene expression program controlled by PI 3-kinase signaling in proliferating cells has not been elucidated.

RESULTS

We used microarray analyses to characterize the changes in gene expression resulting from inhibition of PI 3-kinase in proliferating cells. The genes regulated by inhibition of PI 3-kinase in proliferating cells were distinct from genes induced by growth factor stimulation of quiescent cells and highly enriched in genes that regulate programmed cell death. Computational analyses followed by chromatin immunoprecipitations demonstrated FOXO binding to both previously known and novel sites in promoter regions of approximately one-third of the up-regulated genes, consistent with activation of FOXO1 and FOXO3a in response to inhibition of PI 3-kinase. NFkappaB binding sites were similarly identified in promoter regions of over one-third of the down-regulated genes. RelB was constitutively bound to promoter regions in cells maintained in serum, however binding decreased following PI 3-kinase inhibition, indicating that PI 3-kinase signaling activates NFkappaB via the non-canonical pathway in proliferating cells. Approximately 70% of the genes targeted by FOXO and NFkappaB regulate cell proliferation and apoptosis, including several regulators of apoptosis that were not previously known to be targeted by these transcription factors.

CONCLUSION

PI 3-kinase signaling in proliferating cells regulates a novel transcriptional program that is highly enriched in genes that regulate apoptosis. At least one-third of these genes are regulated either by FOXO transcription factors, which are activated following PI 3-kinase inhibition, or by RelB, which is activated by PI 3-kinase via the non-canonical pathway in proliferating cells.

E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints

The E2F transcription factor family is known to play a key role in the timely expression of genes required for cell cycle progression and proliferation, but only a few E2F target genes have been identified. We explored the possibility that E2F regulators play a broader role by identifying additional genes bound by E2F in living human cells. A protocol was developed to identify genomic binding sites for DNA-binding factors in mammalian cells that combines immunoprecipitation of cross-linked protein-DNA complexes with DNA microarray analysis. Among approximately 1200 genes expressed during cell cycle entry, we found that the promoters of 127 were bound by the E2F4 transcription factor in primary fibroblasts. A subset of these targets was also bound by E2F1. Most previously identified target genes known to have roles in DNA replication and cell cycle control and represented on the microarray were confirmed by this analysis. We also identified a remarkable cadre of genes with no previous connection to E2F regulation, including genes that encode components of the DNA damage checkpoint and repair pathways, as well as factors involved in chromatin assembly/condensation, chromosome segregation, and the mitotic spindle checkpoint. Our data indicate that E2F directly links cell cycle progression with the coordinate regulation of genes essential for both the synthesis of DNA as well as its surveillance.