Ximiao He

CG methylated microarrays identify a novel methylated sequence bound by the CEBPB|ATF4 heterodimer that is active in vivo

To evaluate the effect of CG methylation on DNA binding of sequence-specific B-ZIP transcription factors (TFs) in a high-throughput manner, we enzymatically methylated the cytosine in the CG dinucleotide on protein binding microarrays. Two Agilent DNA array designs were used. One contained 40,000 features using de Bruijn sequences where each 8-mer occurs 32 times in various positions in the DNA sequence. The second contained 180,000 features with each CG containing 8-mer occurring three times. The first design was better for identification of binding motifs, while the second was better for quantification. Using this novel technology, we show that CG methylation enhanced binding for CEBPA and CEBPB and inhibited binding for CREB, ATF4, JUN, JUND, CEBPD, and CEBPG. The CEBPB|ATF4 heterodimer bound a novel motif CGAT|GCAA 10-fold better when methylated. The electrophoretic mobility shift assay (EMSA) confirmed these results. CEBPB ChIP-seq data using primary female mouse dermal fibroblasts with 50× methylome coverage for each strand indicate that the methylated sequences well-bound on the arrays are also bound in vivo. CEBPB bound 39% of the methylated canonical 10-mers ATTGC|GCAAT in the mouse genome. After ATF4 protein induction by thapsigargin which results in ER stress, CEBPB binds methylated CGAT|GCAA in vivo, recapitulating what was observed on the arrays. This methodology can be used to identify new methylated DNA sequences preferentially bound by TFs, which may be functional in vivo.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding

BackgroundChromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS.ResultsWe used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGAC/GTCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---T/CGT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action.ConclusionOur analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.

Methylated Cytosines Mutate to Transcription Factor Binding Sites that Drive Tetrapod Evolution.

Abstract In mammals, the cytosine in CG dinucleotides is typically methylated producing 5-methylcytosine (5mC), a chemically less stable form of cytosine that can spontaneously deaminate to thymidine resulting in a T•G mismatched base pair. Unlike other eukaryotes that efficiently repair this mismatched base pair back to C•G, in mammals, 5mCG deamination is mutagenic, sometimes producing TG dinucleotides, explaining the depletion of CG dinucleotides in mammalian genomes. It was suggested that new TG dinucleotides generate genetic diversity that may be critical for evolutionary change. We tested this conjecture by examining the DNA sequence properties of regulatory sequences identified by DNase I hypersensitive sites (DHSs) in human and mouse genomes. We hypothesized that the new TG dinucleotides generate transcription factor binding sites (TFBS) that become tissue-specific DHSs (TS-DHSs). We find that 8-mers containing the CG dinucleotide are enriched in DHSs in both species. However, 8-mers containing a TG and no CG dinucleotide are preferentially enriched in TS-DHSs when compared with 8-mers with neither a TG nor a CG dinucleotide. The most enriched 8-mer with a TG and no CG dinucleotide in tissue-specific regulatory regions in both genomes is the AP-1 motif (TGAC/GTCAN), and we find evidence that TG dinucleotides in the AP-1 motif arose from CG dinucleotides. Additional TS-DHS-enriched TFBS containing the TG/CA dinucleotide are the E-Box motif (GCAGCTGC), the NF-1 motif (GGCA—TGCC), and the GR (glucocorticoid receptor) motif (G-ACA—TGT-C). Our results support the suggestion that cytosine methylation is mutagenic in tetrapods producing TG dinucleotides that create TFBS that drive evolution.

High-resolution genome-wide DNA methylation maps of mouse primary female dermal fibroblasts and keratinocytes

BackgroundGenome-wide DNA methylation at a single nucleotide resolution in different primary cells of the mammalian genome helps to determine the characteristics and functions of tissue-specific hypomethylated regions (TS-HMRs). We determined genome-wide cytosine methylation maps at 91X and 36X coverage of newborn female mouse primary dermal fibroblasts and keratinocytes and compared with mRNA-seq gene expression data.ResultsThese high coverage methylation maps were used to identify HMRs in both cell types. A total of 2.91% of the genome are in keratinocyte HMRs, and 2.15% of the genome are in fibroblast HMRs with 1.75% being common. Half of the TS-HMRs are extensions of common HMRs, and the remaining are unique TS-HMRs. Four levels of CG methylation are observed: 1) total unmethylation for CG dinucleotides in HMRs in CGIs that are active in all tissues; 2) 10% to 40% methylation for TS-HMRs; 3) 60% methylation for TS-HMRs in cells types where they are not in HMRs; and 4) 70% methylation for the nonfunctioning part of the genome. SINE elements are depleted inside the TS-HMRs, while highly enriched in the surrounding regions. Hypomethylation at the last exon shows gene repression, while demethylation toward the gene body positively correlates with gene expression. The overlapping HMRs have a more complex relationship with gene expression. The common HMRs and TS-HMRs are each enriched for distinct Transcription Factor Binding Sites (TFBS). C/EBPβ binds to methylated regions outside of HMRs while CTCF prefers to bind in HMRs, highlighting these two parts of the genome and their potential interactions.ConclusionsKeratinocytes and fibroblasts are of epithelial and mesenchymal origin. High-resolution methylation maps in these two cell types can be used as reference methylomes for analyzing epigenetic mechanisms in several diseases including cancer.Please see related article at the following link: http://www.epigeneticsandchromatin.com/content/7/1/34

Overlapping ETS and CRE Motifs ((G/C)CGGAAGTGACGTCA) preferentially bound by GABPα and CREB proteins.

Previously, we identified 8-bps long DNA sequences (8-mers) that localize in human proximal promoters and grouped them into known transcription factor binding sites (TFBS). We now examine split 8-mers consisting of two 4-mers separated by 1-bp to 30-bps (X4-N1-30-X4) to identify pairs of TFBS that localize in proximal promoters at a precise distance. These include two overlapping TFBS: the ETS⇔ETS motif (C/GCCGGAAGCGGAA) and the ETS⇔CRE motif (C/GCGGAAGTGACGTCAC). The nucleotides in bold are part of both TFBS. Molecular modeling shows that the ETS⇔CRE motif can be bound simultaneously by both the ETS and the B-ZIP domains without protein-protein clashes. The electrophoretic mobility shift assay (EMSA) shows that the ETS protein GABPα and the B-ZIP protein CREB preferentially bind to the ETS⇔CRE motif only when the two TFBS overlap precisely. In contrast, the ETS domain of ETV5 and CREB interfere with each other for binding the ETS⇔CRE. The 11-mer (CGGAAGTGACG), the conserved part of the ETS⇔CRE motif, occurs 226 times in the human genome and 83% are in known regulatory regions. In vivo GABPα and CREB ChIP-seq peaks identified the ETS⇔CRE as the most enriched motif occurring in promoters of genes involved in mRNA processing, cellular catabolic processes, and stress response, suggesting that a specific class of genes is regulated by this composite motif.

The Epstein-Barr Virus B-ZIP Protein Zta Recognizes Specific DNA Sequences Containing 5-Methylcytosine and 5-Hydroxymethylcytosine

The Epstein-Barr virus (EBV) B-ZIP transcription factor Zta binds to many DNA sequences containing methylated CG dinucleotides. Using protein binding microarrays (PBMs), we analyzed the sequence specific DNA binding of Zta to four kinds of double-stranded DNA (dsDNA): (1) DNA containing cytosine in both strands, (2) DNA with 5-methylcytosine (5mC) in one strand and cytosine in the second strand, (3) DNA with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand, and (4) DNA in which both cytosines in all CG dinucleotides contain 5mC. We compared these data to PBM data for three additional B-ZIP proteins (CREB1 and CEBPB homodimers and cJun|cFos heterodimers). With cytosine, Zta binds the TRE motif TGAC/GTCA as previously reported. With CG dinucleotides containing 5mC on both strands, many TRE motif variants containing a methylated CG dinucleotide at two positions in the motif, such as MGAGTCA and TGAGMGA (where M = 5mC), were preferentially bound. 5mC inhibits binding of Zta to both TRE motif half-sites GTCA and CTCA. Like the CREB1 homodimer, the Zta homodimer and the cJun|cFos heterodimer more strongly bind the C/EBP half-site tetranucleotide GCAA when it contains 5mC. Zta also binds dsDNA sequences containing 5hmC in one strand, although the effect is less dramatic than that observed for 5mC. Our results identify new DNA sequences that are well-bound by the viral B-ZIP protein Zta only when they contain 5mC or 5hmC, uncovering the potential for discovery of new viral and host regulatory programs controlled by EBV.

GABPα Binding to Overlapping ETS and CRE DNA Motifs Is Enhanced by CREB1: Custom DNA Microarrays

To achieve proper spatiotemporal control of gene expression, transcription factors cooperatively assemble onto specific DNA sequences. The ETS domain protein monomer of GABPα and the B-ZIP domain protein dimer of CREB1 cooperatively bind DNA only when the ETS (C/GCGGAAGT) and CRE (GTGACGTCAC) motifs overlap precisely, producing the ETS↔CRE motif (C/GCGGAAGTGACGTCAC). We designed a Protein Binding Microarray (PBM) with 60-bp DNAs containing four identical sectors, each with 177,440 features that explore the cooperative interactions between GABPα and CREB1 upon binding the ETS↔CRE motif. The DNA sequences include all 15-mers of the form C/GCGGA—–CG—, the ETS↔CRE motif, and all single nucleotide polymorphisms (SNPs), and occurrences in the human and mouse genomes. CREB1 enhanced GABPα binding to the canonical ETS↔CRE motif CCGGAAGT two-fold, and up to 23-fold for several SNPs at the beginning and end of the ETS motif, which is suggestive of two separate and distinct allosteric mechanisms of cooperative binding. We show that the ETS-CRE array data can be used to identify regions likely cooperatively bound by GABPα and CREB1 in vivo, and demonstrate their ability to identify human genetic variants that might inhibit cooperative binding.

Nucleosomes are enriched at the boundaries of hypomethylated regions (HMRs) in mouse dermal fibroblasts and keratinocytes

BackgroundThe interplay between epigenetic modifications and chromatin structure are integral to our understanding of genome function. Methylation of cytosine (5mC) at CG dinucleotides, traditionally associated with transcriptional repression, is the most highly studied chemical modification of DNA, occurring at over 70% of all CG dinucleotides in the genome. Hypomethylated regions (HMRs) often occur in CG islands (CGIs), however, they also occur outside of CGIs and function as cell-type specific enhancers. During the process of differentiation, reorganization of chromatin and nucleosome arrangement at regulatory regions is thought to occur in order for the establishment of cell-type specific transcriptional programs. However, the specifics regarding the organization of nucleosomes at HMRs and the potential mechanisms regulating nucleosome occupancy in these regions are unknown. Here, we have investigated nucleosome organization around hypomethylated regions (HMRs) identified in two mouse primary cells.ResultsMicroccocal nuclease (MNase) digested mononucleosomes from primary cultures of new-born female mouse dermal fibroblasts and keratinocytes were mapped and compared to the HMRs obtained from single base-pair resolution methylomes. In both cell types, we find that nucleosomes are enriched at HMR boundaries. In contrast to the nucleosomes found at boundaries of HMRs in CGIs, HMRs outside of CGIs are calculated to be preferentially bound by nucleosomes, with phased nucleosomes propagating into the methylated region. Nucleosomes are enriched at the tissue-specific HMRs (TS-HMR) boundaries in both cell types suggesting that nucleosome organization surrounding HMR boundaries is independent of methylation status. In addition, we find potential transcription factor (TF) binding sites (E-box motifs) enriched in non-CGI TS-HMR boundaries.ConclusionsOur results show that intrinsic nucleosome occupancy score (INOS) positively correlate with the nucleosome organization surrounding non-CGI TS-HMRs, suggesting that DNA sequence plays a role in the establishment of HMRs in the genome. Since nucleosomes impact all processes involving the genome, our results provide a link between epigenetic modifications, chromatin structure, and regulatory function.