The homeobox transcription factor DUXBL controls exit from totipotency.

In mice, exit from the totipotent two-cell (2C) stage embryo requires silencing of the 2C-associated transcriptional program. However, the molecular mechanisms involved in this process remain poorly understood. Here we demonstrate that the 2C-specific transcription factor double homeobox protein (DUX) mediates an essential negative feedback loop by inducing the expression of DUXBL to promote this silencing. We show that DUXBL gains accessibility to DUX-bound regions specifically upon DUX expression. Furthermore, we determine that DUXBL interacts with TRIM24 and TRIM33, members of the TRIM superfamily involved in gene silencing, and colocalizes with them in nuclear foci upon DUX expression. Importantly, DUXBL overexpression impairs 2C-associated transcription, whereas Duxbl inactivation in mouse embryonic stem cells increases DUX-dependent induction of the 2C-transcriptional program. Consequently, DUXBL deficiency in embryos results in sustained expression of 2C-associated transcripts leading to early developmental arrest. Our study identifies DUXBL as an essential regulator of totipotency exit enabling the first divergence of cell fates.

Decoding complexity in biomolecular recognition of DNA i-motifs with microarrays.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

Direct ionic stress sensing and mitigation by the transcription factor NFAT5.

Homeostatic control of intracellular ionic strength is essential for protein, organelle and genome function, yet mechanisms that sense and enable adaptation to ionic stress remain poorly understood in animals. We find that the transcription factor NFAT5 directly senses solution ionic strength using a C-terminal intrinsically disordered region. Both in intact cells and in a purified system, NFAT5 forms dynamic, reversible biomolecular condensates in response to increasing ionic strength. This self-associative property, conserved from insects to mammals, allows NFAT5 to accumulate in the nucleus and activate genes that restore cellular ion content. Mutations that reduce condensation or those that promote aggregation both reduce NFAT5 activity, highlighting the importance of optimally tuned associative interactions. Remarkably, human NFAT5 alone is sufficient to reconstitute a mammalian transcriptional response to ionic or hypertonic stress in yeast. Thus NFAT5 is both the sensor and effector of a cell-autonomous ionic stress response pathway in animal cells.

A novel GATA2 distal enhancer mutation results in MonoMAC syndrome in 2 second cousins.

Mutations in the transcription factor GATA2 can cause MonoMAC syndrome, a GATA2 deficiency disease characterized by several findings, including disseminated nontuberculous mycobacterial infections, severe deficiencies of monocytes, natural killer cells, and B lymphocytes, and myelodysplastic syndrome. GATA2 mutations are found in âˆ¼90% of patients with a GATA2 deficiency phenotype and are largely missense mutations in the conserved second zinc-finger domain. Mutations in an intron 5 regulatory enhancer element are also well described in GATA2 deficiency. Here, we present a multigeneration kindred with the clinical features of GATA2 deficiency but lacking an apparent GATA2 mutation. Whole genome sequencing revealed a unique adenine-to-thymine variant in the GATA2 -110 enhancer 116,855 bp upstream of the GATA2 ATG start site. The mutation creates a new E-box consensus in position with an existing GATA-box to generate a new hematopoietic regulatory composite element. The mutation segregates with the disease in several generations of the family. Cell type-specific allelic imbalance of GATA2 expression was observed in the bone marrow of a patient with higher expression from the mutant-linked allele. Allele-specific overexpression of GATA2 was observed in CRISPR/Cas9-modified HL-60 cells and in luciferase assays with the enhancer mutation. This study demonstrates overexpression of GATA2 resulting from a single nucleotide change in an upstream enhancer element in patients with MonoMAC syndrome. Patients in this study were enrolled in the National Institute of Allergy and Infectious Diseases clinical trial and the National Cancer Institute clinical trial (both trials were registered at www.clinicaltrials.gov as #NCT01905826 and #NCT01861106, respectively).

DECODING COMPLEXITY IN BIOMOLECULAR RECOGNITION OF DNA I-MOTIFS.

DNA i-motifs (iMs) are non-canonical C-rich secondary structures implicated in numerous cellular processes. Though iMs exist throughout the genome, our understanding of iM recognition by proteins or small molecules is limited to a few examples. We designed a DNA microarray containing 10,976 genomic iM sequences to examine the binding profiles of four iM-binding proteins, mitoxantrone, and the iMab antibody. iMab microarray screens demonstrated that pH 6.5, 5% BSA buffer was optimal, and fluorescence was correlated with iM C-tract length. hnRNP K broadly recognizes diverse iM sequences, favoring 3-5 cytosine repeats flanked by thymine-rich loops of 1-3 nucleotides. Array binding mirrored public ChIP-Seq datasets, in which 35% of well-bound array iMs are enriched in hnRNP K peaks. In contrast, other reported iM-binding proteins had weaker binding or preferred G-quadruplex (G4) sequences instead. Mitoxantrone broadly binds both shorter iMs and G4s, consistent with an intercalation mechanism. These results suggest that hnRNP K may play a role in iM-mediated regulation of gene expression in vivo, whereas hnRNP A1 and ASF/SF2 are possibly more selective in their binding preferences. This powerful approach represents the most comprehensive investigation of how biomolecules selectively recognize genomic iMs to date.

A Knock-In Mouse Model of Thymoma With the GTF2I L424H Mutation.

INTRODUCTION: The pathogenesis of thymic epithelial tumors remains largely unknown. We previously identified GTF2I L424H as the most frequently recurrent mutation in thymic epithelial tumors. Nevertheless, the precise role of this mutation in tumorigenesis of thymic epithelial cells is unclear. METHODS: To investigate the role of GTF2I L424H mutation in thymic epithelial cells in vivo, we generated and characterized a mouse model in which the Gtf2i L424H mutation was conditionally knocked-in in the Foxn1+ thymic epithelial cells. Digital spatial profiling was performed on thymomas and normal thymic tissues with GeoMx-mouse whole transcriptome atlas. Immunohistochemistry staining was performed using both mouse tissues and human thymic epithelial tumors. RESULTS: We observed that the Gtf2i mutation impairs development of the thymic medulla and maturation of medullary thymic epithelial cells in young mice and causes tumor formation in the thymus of aged mice. Cell cycle-related pathways, such as E2F targets and MYC targets, are enriched in the tumor epithelial cells. Results of gene set variation assay analysis revealed that gene signatures of cortical thymic epithelial cells and thymic epithelial progenitor cells are also enriched in the thymomas of the knock-in mice, which mirrors the human counterparts in The Cancer Genome Atlas database. Immunohistochemistry results revealed similar expression pattern of epithelial cell markers between mouse and human thymomas. CONCLUSIONS: We have developed and characterized a novel thymoma mouse model. This study improves knowledge of the molecular drivers in thymic epithelial cells and provides a tool for further study of the biology of thymic epithelial tumors and for development of novel therapies.

Spatial meta-transcriptomics reveal associations of intratumor bacteria burden with lung cancer cells showing a distinct oncogenic signature.

BACKGROUND: The lung intratumor microbiome influences lung cancer tumorigenesis and treatment responses, but detailed data on the extent, location, and effects of microbes within lung tumors are missing, information needed for improved prognosis and treatment. METHODS: To address this gap, we developed a novel spatial meta-transcriptomic method simultaneously detecting the expression level of 1,811 host genes and 3 microbe targets (bacteria, fungi, and cytomegalovirus). After rigorous validation, we analyzed the spatial meta-transcriptomic profiles of tumor cells, T cells, macrophages, other immune cells, and stroma in surgically resected tumor samples from 12 patients with early-stage lung cancer. RESULTS: Bacterial burden was significantly higher in tumor cells compared with T cells, macrophages, other immune cells, and stroma. This burden increased from tumor-adjacent normal lung and tertiary lymphoid structures to tumor cells to the airways, suggesting that lung intratumor bacteria derive from the latter route of entry. Expression of oncogenic ß-catenin was strongly correlated with bacterial burden, as were tumor histological subtypes and environmental factors. CONCLUSIONS: Intratumor bacteria were enriched with tumor cells and associated with multiple oncogenic pathways, supporting a rationale for reducing the local intratumor microbiome in lung cancer for patient benefit. TRIAL REGISTRATION NUMBER: NCT00242723, NCT02146170.

Altering the Double-Stranded DNA Specificity of the bZIP Domain of Zta with Site-Directed Mutagenesis at N182.

Zta, the Epstein-Barr virus bZIP transcription factor (TF), binds both unmethylated and methylated double-stranded DNA (dsDNA) in a sequence-specific manner. We studied the contribution of a conserved asparagine (N182) to sequence-specific dsDNA binding to four types of dsDNA: (i) dsDNA with cytosine in both strands ((DNA(C|C)), (ii, iii) dsDNA with 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and cytosine in the second strand ((DNA(5mC|C) and DNA(5hmC|C)), and (iv) dsDNA with methylated cytosine in both strands in all CG dinucleotides ((DNA(5mCG)). We replaced asparagine with five similarly sized amino acids (glutamine (Q), serine (S), threonine (T), isoleucine (I), or valine (V)) and used protein binding microarrays to evaluate sequence-specific dsDNA binding. Zta preferentially binds the pseudo-palindrome TRE (AP1) motif (T-4G-3A-2G/C 0T2C3A4 ). Zta (N182Q) changes binding to A3 in only one half-site. Zta(N182S) changes binding to G3 in one or both halves of the motif. Zta(N182S) and Zta(N182Q) have 34- and 17-fold weaker median dsDNA binding, respectively. Zta(N182V) and Zta(N182I) have increased binding to dsDNA(5mC|C). Molecular dynamics simulations rationalize some of these results, identifying hydrogen bonds between glutamine and A3 , but do not reveal why serine preferentially binds G3 , suggesting that entropic interactions may mediate this new binding specificity.

ASXL1 and STAG2 are common mutations in GATA2 deficiency patients with bone marrow disease and myelodysplastic syndrome.

Patients with GATA2 deficiencyharbor de novo or inherited germline mutations in the GATA2 transcription factor gene, predisposing them to myeloid malignancies. There is considerable variation in disease progression, even among family members with the same mutation in GATA2. We investigated somatic mutations in 106 patients with GATA2 deficiency to identify acquired mutations that are associated with myeloid malignancies. Myelodysplastic syndrome (MDS) was the most common diagnosis (∼44%), followed by GATA2 bone marrow immunodeficiency disorder (G2BMID; ∼37%). Thirteen percent of the cohort had GATA2 mutations but displayed no disease manifestations. There were no correlations between age or sex with disease progression or survival. Cytogenetic analyses showed a high incidence of abnormalities (∼43%), notably trisomy 8 (∼23%) and monosomy 7 (∼12%), but the changes did not correlate with lower survival. Somatic mutations in ASXL1 and STAG2 were detected in ∼25% of patients, although the mutations were rarely concomitant. Mutations in DNMT3A were found in ∼10% of patients. These somatic mutations were found similarly in G2BMID and MDS, suggesting clonal hematopoiesis in early stages of disease, before the onset of MDS. ASXL1 mutations conferred a lower survival probability and were more prevalent in female patients. STAG2 mutations also conferred a lower survival probability, but did not show a statistically significant sex bias. There was a conspicuous absence of many commonly mutated genes associated with myeloid malignancies, including TET2, IDH1/2, and the splicing factor genes. Notably, somatic mutations in chromatin-related genes and cohesin genes characterized disease progression in GATA2 deficiency.

The ETS transcription factor ERF controls the exit from the naïve pluripotent state in a MAPK-dependent manner.

The naïve epiblast transitions to a pluripotent primed state during embryo implantation. Despite the relevance of the FGF pathway during this period, little is known about the downstream effectors regulating this signaling. Here, we examined the molecular mechanisms coordinating the naïve to primed transition by using inducible ESC to genetically eliminate all RAS proteins. We show that differentiated RASKO ESC remain trapped in an intermediate state of pluripotency with naïve-associated features. Elimination of the transcription factor ERF overcomes the developmental blockage of RAS-deficient cells by naïve enhancer decommissioning. Mechanistically, ERF regulates NANOG expression and ensures naïve pluripotency by strengthening naïve transcription factor binding at ESC enhancers. Moreover, ERF negatively regulates the expression of the methyltransferase DNMT3B, which participates in the extinction of the naïve transcriptional program. Collectively, we demonstrated an essential role for ERF controlling the exit from naïve pluripotency in a MAPK-dependent manner during the progression to primed pluripotency.

CTCF is a barrier for 2C-like reprogramming.

Totipotent cells have the ability to generate embryonic and extra-embryonic tissues. Interestingly, a rare population of cells with totipotent-like potential, known as 2 cell (2C)-like cells, has been identified within ESC cultures. They arise from ESC and display similar features to those found in the 2C embryo. However, the molecular determinants of 2C-like conversion have not been completely elucidated. Here, we show that the CCCTC-binding factor (CTCF) is a barrier for 2C-like reprogramming. Indeed, forced conversion to a 2C-like state by the transcription factor DUX is associated with DNA damage at a subset of CTCF binding sites. Depletion of CTCF in ESC efficiently promotes spontaneous and asynchronous conversion to a 2C-like state and is reversible upon restoration of CTCF levels. This phenotypic reprogramming is specific to pluripotent cells as neural progenitor cells do not show 2C-like conversion upon CTCF-depletion. Furthermore, we show that transcriptional activation of the ZSCAN4 cluster is necessary for successful 2C-like reprogramming. In summary, we reveal an unexpected relationship between CTCF and 2C-like reprogramming.

REL Domain of NFATc2 Binding to Five Types of DNA Using Protein Binding Microarrays.

NFATc2 is a DNA binding protein in the Rel family transcription factors, which binds a CGGAA motif better when both cytosines in the CG dinucleotide are methylated. Using protein binding microarrays (PBMs), we examined the DNA binding of NFATc2 to three additional types of DNA: single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) with either 5-methylcytosine (5mC, M) or 5-hydroxymethylcytosine (5hmC, H) in one strand and a cytosine in the second strand. ATTTCCAC, the complement of the core GGAA motif, is better bound as ssDNA compared to dsDNA. dsDNA containing the 5-mer CGGAA with either 5mC or 5hmC in one DNA strand is bound stronger than CGGAA. In contrast, the reverse complement TTCCG is bound weaker when it contains 5mC. Analysis of the available NFATc2:dsDNA complexes rationalizes these PBM data.

bZIP Dimers CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun Prefer to Bind to Some Double-Stranded DNA Sequences Containing 5-Formylcytosine and 5-Carboxylcytosine.

In mammalian cells, 5-methylcytosine (5mC) occurs in genomic double-stranded DNA (dsDNA) and is enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), then to 5-formylcytosine (5fC), and finally to 5-carboxylcytosine (5caC). These cytosine modifications are enriched in regulatory regions of the genome. The effect of these oxidative products on five bZIP dimers (CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun) binding to five types of dsDNA was measured using protein binding microarrays. The five dsDNAs contain either cytosine in both DNA strands or cytosine in one strand and either 5mC, 5hmC, 5fC, or 5caC in the second strand. Some sequences containing the CEBP half-site GCAA are bound more strongly by all five bZIP domains when dsDNA contains 5mC, 5hmC, or 5fC. dsDNA containing 5caC in some TRE (AP-1)-like sequences, e.g., TGACTAA, is better bound by Zta, ATF3|cJun, and cFos|cJun.

Custom G4 Microarrays Reveal Selective G-Quadruplex Recognition of Small Molecule BMVC: A Large-Scale Assessment of Ligand Binding Selectivity.

G-quadruplexes (G4) are considered new drug targets for human diseases such as cancer. More than 10,000 G4s have been discovered in human chromatin, posing challenges for assessing the selectivity of a G4-interactive ligand. 3,6-bis(1-Methyl-4-vinylpyridinium) carbazole diiodide (BMVC) is the first fluorescent small molecule for G4 detection in vivo. Our previous structural study shows that BMVC binds to the MYC promoter G4 (MycG4) with high specificity. Here, we utilize high-throughput, large-scale custom DNA G4 microarrays to analyze the G4-binding selectivity of BMVC. BMVC preferentially binds to the parallel MycG4 and selectively recognizes flanking sequences of parallel G4s, especially the 3'-flanking thymine. Importantly, the microarray results are confirmed by orthogonal NMR and fluorescence binding analyses. Our study demonstrates the potential of custom G4 microarrays as a platform to broadly and unbiasedly assess the binding selectivity of G4-interactive ligands, and to help understand the properties that govern molecular recognition.

Custom DNA Microarrays Reveal Diverse Binding Preferences of Proteins and Small Molecules to Thousands of G-Quadruplexes.

Single-stranded DNA (ssDNA) containing four guanine repeats can form G-quadruplex (G4) structures. While cellular proteins and small molecules can bind G4s, it has been difficult to broadly assess their DNA-binding specificity. Here, we use custom DNA microarrays to examine the binding specificities of proteins, small molecules, and antibodies across ∼15,000 potential G4 structures. Molecules used include fluorescently labeled pyridostatin (Cy5-PDS, a small molecule), BG4 (Cy5-BG4, a G4-specific antibody), and eight proteins (GST-tagged nucleolin, IGF2, CNBP, FANCJ, PIF1, BLM, DHX36, and WRN). Cy5-PDS and Cy5-BG4 selectively bind sequences known to form G4s, confirming their formation on the microarrays. Cy5-PDS binding decreased when G4 formation was inhibited using lithium or when ssDNA features on the microarray were made double-stranded. Similar conditions inhibited the binding of all other molecules except for CNBP and PIF1. We report that proteins have different G4-binding preferences suggesting unique cellular functions. Finally, competition experiments are used to assess the binding specificity of an unlabeled small molecule, revealing the structural features in the G4 required to achieve selectivity. These data demonstrate that the microarray platform can be used to assess the binding preferences of molecules to G4s on a broad scale, helping to understand the properties that govern molecular recognition.

The bZIP mutant CEBPB (V285A) has sequence specific DNA binding propensities similar to CREB1.

The bZIP homodimers CEBPB and CREB1 bind DNA containing methylated cytosines differently. CREB1 binds stronger to the C/EBP half-site GCAA when the cytosine is methylated. For CEBPB, methylation of the same cytosine does not affect DNA binding. The X-ray structure of CREB1 binding the half site GTCA identifies an alanine in the DNA binding region interacting with the methyl group of T, structurally analogous to the methyl group of methylated C. This alanine is replaced with a valine in CEBPB. To explore the contribution of this amino acid to binding with methylated cytosine of the GCAA half-site, we made the reciprocal mutants CEBPB(V285A) and CREB1(A297V) and used protein binding microarrays (PBM) to examine binding to four types of double-stranded DNA (dsDNA): 1) DNA with cytosine in both strands (DNA(C|C)), 2) DNA with 5-methylcytosine (M) in one strand and cytosine in the second strand (DNA(M|C)), 3) DNA with 5-hydroxymethylcytosine (H) in one strand and cytosine in the second strand (DNA(H|C)), and 4) DNA with both cytosines in all CG dinucleotides containing 5-methylcytosine (DNA(5mCG)). When binding to DNA(C|C), CEBPB (V285A) preferentially binds the CRE consensus motif (TGACGTCA), similar to CREB1. The reciprocal mutant, CREB1(A297V) binds DNA with some similarity to CEBPB, with strongest binding to the methylated PAR site 8-mer TTACGTAA. These data demonstrate that V285 residue inhibits CEBPB binding to methylated cytosine of the GCAA half-site.

GABPα and CREB1 Binding to Double Nucleotide Polymorphisms of Their Consensus Motifs and Cooperative Binding to the Composite ETS ⇔ CRE Motif (ACCGGAAGTGACGTCA).

Previously, cooperative binding of the bZIP domain of CREB1 and the ETS domain of GABPα was observed for the composite DNA ETS ⇔ CRE motif (A 0 C 1 C 2 G 3 G 4 A 5 A 6 G 7 T 8 G 9 A 10 C 11 G 12 T 13 C 14 A 15 ). Single nucleotide polymorphisms (SNPs) at the beginning and end of the ETS motif (ACCGGAAGT) increased cooperative binding. Here, we use an Agilent microarray of 60-mers containing all double nucleotide polymorphisms (DNPs) of the ETS ⇔ CRE motif to explore GABPα and CREB1 binding to their individual motifs and their cooperative binding. For GABPα, all DNPs were bound as if each SNP acted independently. In contrast, CREB1 binding to some DNPs was stronger or weaker than expected, depending on the locations of each SNP. CREB1 binding to DNPs where both SNPs were in the same half site, T 8 G 9 A 10 or T 13 C 14 A 15 , was greater than expected, indicating that an additional SNP cannot destroy binding as much as expected, suggesting that an individual SNP is enough to abolish sequence-specific DNA binding of a single bZIP monomer. If a DNP contains SNPs in each half site, binding is weaker than expected. Similar results were observed for additional ETS and bZIP family members. Cooperative binding between GABPα and CREB1 to the ETS ⇔ CRE motif was weaker than expected except for DNPs containing A 7 and SNPs at the beginning of the ETS motif.

Replacing C189 in the bZIP domain of Zta with S, T, V, or A changes DNA binding specificity to four types of double-stranded DNA.

Zta is a bZIP transcription factor (TF) in the Epstein-Barr virus that binds unmethylated and methylated DNA sequences. Substitution of cysteine 189 of Zta to serine (Zta(C189S)) results in a virus that is unable to execute the lytic cycle, which was attributed to a change in binding to methylated DNA sequences. To learn more about the role of this position in defining sequence-specific DNA binding, we mutated cysteine 189 to four other amino acids, producing Zta(C189S), Zta(C189T), Zta(C189A), and Zta(C189V) mutants. Zta and mutants were used in protein binding microarray (PBM) experiments to evaluate sequence-specific DNA binding to four types of double-stranded DNA (dsDNA): 1) with cytosine in both strands (DNA(C|C)), 2) with 5-methylcytosine (5mC) in one strand and cytosine in the second strand (DNA(5mC|C)), 3) with 5-hydroxymethylcytosine (5hmC) in one strand and cytosine in the second strand (DNA(5hmC|C)), and 4) with both cytosines in all CG dinucleotides containing 5mC (DNA(5mCG)). Zta(C189S) and Zta(C189T) bound the TRE (AP-1) motif (TGAG/CTCA) more strongly than wild-type Zta, while binding to other sequences, including the C/EBP half site GCAA was reduced. Binding of Zta(C189S) and Zta(C189T) to DNA containing modified cytosines (DNA(5mC|C), DNA(5hmC|C), and DNA(5mCG)) was reduced compared to Zta. Zta(C189A) and Zta(C189V) had higher non-specific binding to all four types of DNA. Our data suggests that position C189 in Zta impacts sequence-specific binding to DNA containing modified and unmodified cytosine.

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.

5-Methylcytosine (5mC) and 5-Hydroxymethylcytosine (5hmC) Enhance the DNA Binding of CREB1 to the C/EBP Half-Site Tetranucleotide GCAA.

In human and mouse stem cells and brain, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) can occur outside of CG dinucleotides. Using protein binding microarrays (PBMs) containing 60-mer DNA probes, we evaluated the effect of 5mC and 5hmC on one DNA strand on the double-stranded DNA binding of the mouse B-ZIP transcription factors (TFs) CREB1, ATF1, and JUND. 5mC inhibited binding of CREB1 to the canonical CRE half-site |GTCA but enhanced binding to the C/EBP half-site |GCAA. 5hmC inhibited binding of CREB1 to all 8-mers except TGAT|GCAA, where binding is enhanced. We observed similar DNA binding patterns with ATF1, a closely related B-ZIP domain. In contrast, both 5mC and 5hmC inhibited binding of JUND. These results identify new DNA sequences that are well-bound by CREB1 and ATF1 only when they contain 5mC or 5hmC. Analysis of two X-ray structures examines the consequences of 5mC and 5hmC on DNA binding by CREB and FOS|JUN.