GATA3-dependent cellular reprogramming requires activation-domain dependent recruitment of a chromatin remodeler
© Takaku et al. 2016
Received: 30 October 2015
Accepted: 11 February 2016
Published: 27 February 2016
Transcription factor-dependent cellular reprogramming is integral to normal development and is central to production of induced pluripotent stem cells. This process typically requires pioneer transcription factors (TFs) to induce de novo formation of enhancers at previously closed chromatin. Mechanistic information on this process is currently sparse.
Here we explore the mechanistic basis by which GATA3 functions as a pioneer TF in a cellular reprogramming event relevant to breast cancer, the mesenchymal to epithelial transition (MET). In some instances, GATA3 binds previously inaccessible chromatin, characterized by stable, positioned nucleosomes where it induces nucleosome eviction, alters local histone modifications, and remodels local chromatin architecture. At other loci, GATA3 binding induces nucleosome sliding without concomitant generation of accessible chromatin. Deletion of the transactivation domain retains the chromatin binding ability of GATA3 but cripples chromatin reprogramming ability, resulting in failure to induce MET.
These data provide mechanistic insights into GATA3-mediated chromatin reprogramming during MET, and suggest unexpected complexity to TF pioneering. Successful reprogramming requires stable binding to a nucleosomal site; activation domain-dependent recruitment of co-factors including BRG1, the ATPase subunit of the SWI/SNF chromatin remodeling complex; and appropriate genomic context. The resulting model provides a new conceptual framework for de novo enhancer establishment by a pioneer TF.
Transcription factors (TFs) are utilized throughout life as a means of regulating transcriptional output. The human genome encodes somewhere between 1,500 and 2,000 TFs , making them one of the more abundant protein classes. These proteins recognize DNA motifs in and around genes and alter RNA synthesis through a variety of mechanisms [2, 3]. Classification based on structures of their DNA-binding domains provides cardinal insights into their modes of DNA recognition. While structural studies have advanced our understanding of how TFs interact with DNA, genomic studies of their localization indicate that the vast majority of potential TF binding sites are not occupied in vivo [4–6].
In eukaryotes, genomic DNA is packaged into nucleosomes, the fundamental unit of chromatin structure. Nucleosomal architecture presents a barrier to successful protein-DNA interactions for many TFs [7–9]. Pioneer TFs represent an exception to this general rule; these proteins have the ability to interact productively with DNA in a nucleosomal context. Pioneer factors have at least three fundamental properties as recently described by Iwafuchi-Doi and Zaret : (1) they can interact with their cognate recognition sequences in nuclease-resistant chromatin prior to activation of transcription; (2) they increase local chromatin accessibility as a prelude to productive binding of other factors; and (3) they have integral roles in the establishment of cell lineage. FOXA1 represents a prototypical example of a pioneer TF; FOXA1 can bind nucleosomal DNA in vitro, and target closed chromatin regions in vivo which subsequently become accessible and permissive for productive binding of other factors [5, 11–13]. FOXA1 is required for multiple developmental lineage commitment choices including those leading to development of liver, lung, and mammary gland [13–17]. Other TFs involved in reprogramming of somatic cells to pluripotency such as Oct4, Sox2, and Klf4, also function as pioneer factors utilizing unique adaptations of their DNA recognition surfaces to accommodate nucleosome structure, which likely represents a common genome interrogation mechanism of the pioneer class [18, 19]. FOXA1 has structural similarity to linker histone [20, 21], therefore it is thought to compete with linker histone . However validated pioneer factors exhibit large structural diversity in sequence-recognition domains, mechanism(s) underlying chromatin opening may vary.
DNA-binding proteins that have zinc finger domains compose the most abundant protein family in eukaryotic genomes . The vertebrate GATA family proteins, which possess two zinc-finger domains, play critical roles across development and differentiation, contributing to lineage commitment. They also possess several hallmark features, as a family, of pioneer TFs: several family members can bind to nucleosomal DNA and induce accessibility, both in vivo and in vitro [11, 23–29]. Importantly the DNA-binding domains of GATA family proteins are structurally distinct from FOXA1, implying that mechanism(s) to create open chromatin may differ. While elegant evidence supports an intrinsic capacity of FOXA1 to remodel chromatin without the assistance of co-factors , it remains unclear how and whether GATA proteins induce similar structural alterations.
In this study, we focus on the human GATA3 protein, which is the only human GATA family member whose DNA-binding domain has been crystallized . GATA3 is a key regulator of multiple developmental pathways including mammary epithelial cell differentiation, T lymphocyte development, and trophoblast development [31–34]. GATA3 is implicated in breast cancer progression, and recently has been identified as one of the most frequently mutated genes in breast cancer . A considerable number of clinical studies have demonstrated that reduced expression level of GATA3 correlates with poor prognosis . Consistent with the in vivo data, GATA3 exerts metastasis suppressive functions when expressed exogenously in a human basal-like breast cancer cell line by inducing the mesenchymal-to-epithelial transition (MET) [37–39]. Although this transition is accompanied by gene expression reprogramming including microRNA induction , the molecular mechanisms underlying this cell reprogramming process at the chromatin level are incompletely understood. In other contexts, GATA3 is known to reprogram fibroblasts to induced pluripotent stem cells (iPSCs)  as well as induced trophoblast stem cells (iTSCs) [41, 42] in collaboration with other TFs, implying a potential role for GATA3 in chromatin reprogramming.
In order to investigate the impact(s) of GATA3 on chromatin, we utilized the MET model in the MDA-MB-231 breast cancer cell line. Consistent with previous results, GATA3 expression induced mesenchymal-to-epithelial transition at the cellular and molecular level and mediated reprogramming of chromatin at genes integral to MET. Our studies provide experimental evidence that GATA3 can direct alterations in local chromatin during the MET reprogramming process, including pioneer factor activity, that differ depending on local context and that these alterations proceed through an intermediate consisting of nucleosomal DNA bound by a pioneer factor.
GATA3 binding sites acquire open chromatin structure and enhancer-like histone modifications
Pioneer factors are characterized by their ability to bind to previously ‘closed’ chromatin and elicit a chromatin state change. To evaluate the impact of GATA3 on chromatin structure, we measured chromatin accessibility in MDA-MB-231 cells with and without GATA3 expression using ATAC-seq . At the individual locus level, many GATA3-binding sites including MET associated gene loci exhibited increased accessibility in GATA3-expressing cells (hereafter MDA231-GATA3 cells) (Fig 1b; Additional file 1: Figure S1F). Across the genome, TSS-distal peaks (greater than 1 kb from nearest TSS) exhibited a significant difference (P values < e-10) between control and MDA231-GATA3 cells (Fig. 1c). Randomly selected GATA motif containing regions showed similar levels of transposase sensitivity between control and GATA3-expressing cells, indicating that overall transposase integration efficiencies were similar (Additional file 1: Figure S1G). To dissect the function of chromatin domains bound by GATA3, we investigated histone modifications typically associated with active chromatin at enhancers, such as H3K4 monomethylation (H3K4me1) and H3K27 acetylation (H3K27ac). Similar to the ATAC-seq results, H3K4me1 and H3K27ac levels in MDA-231-GATA3 cells were drastically increased at TSS distal GATA3 peaks (Fig. 1d, e).
Chromatin structure based classification reveals a pioneer class of binding sites
To elucidate the correlation between gene expression and GATA3 binding, we assigned GATA3 peaks to target genes using a simple distance metric, and measured the impact of GATA3 on gene expression using RNA-seq. GATA3-binding events in category G1 were enriched for increased gene expression, and the ratio of the genes categorized as upregulated was higher in G1 than that in G2 and G3 (Additional file 1: Figure S3A-D). Interestingly, target genes in G1 contain transcripts indicative of reprogramming of MDA-MB-231 to an epithelial identity, such as KRT8, KRT18, CDH1, and CLDN3 (Additional file 2: Table S1). When analyzed over time following exogenous expression of GATA3 using lentivirus, a subset of those genes were upregulated within 48-96 h after the induction of GATA3 (Additional file 1: Figure S3E). Functional pathway analysis of genes found in G1 featured enrichment in genes integral to the MET and to maintenance of epithelial identity (Additional file 1: Figure S3F). These functional categories were not enriched in G3 (Additional file 1: Figure S3G, Additional file 2: Table S1), suggesting that GATA3 binding and chromatin opening at G1 sites is integral to the MET, consistent with function as a pioneer.
GATA3 remodels local nucleosome structure
Mechanistically, it was unclear why GATA3 binding at previously inaccessible loci, as detected by ChIP-seq, should sometimes, but not always, result in formation of accessible chromatin. Therefore, we sought to understand why GATA3 could open chromatin structure only at G1 sites but not at G3 sites. To address this question, we compared the enrichment of GATA3 between these groups. The open chromatin loci within categories G1 and G2 had similar enrichment of GATA3; the transposase-inaccessible loci bound by GATA3 in category G3 clearly displayed lower enrichment (Fig. 2e; Additional file 1: Figure S2J), suggesting a fundamental difference in the interaction of GATA3 with its cognate recognition element at these two classes of sites. These observed differences in ChIP-seq read level could result from differences in inherent affinity driven by local DNA sequence. Alternatively, they could result from interaction of GATA3 with nucleosomal DNA (G3) fundamentally differing from interaction with accessible chromatin (G1 and G2). Finally, the differences in enrichment could result from differential association with co-factors, such as the Friend of GATA proteins [45, 46] or other partners. Investigating the genes differentially expressed in our system (Additional file 1: Figure S1C), 108 known TFs were upregulated in MDA-231-GATA3 cells, and two TF motifs were co-enriched in category G1, indicating the potential for cooperative action at a subset of sites (Additional file 3: Table S2).
Qualitatively, consensus motifs derived from the loci in G1 did not differ substantially from those in G3 (Additional file 1: Table S3), suggesting differences in binding do not result from deviation of recognition elements from consensus. Accordingly, we hypothesized that the observed difference in enrichment might result from either difference in co-factor association or from difference in substrate quality when comparing nucleosomal DNA to disrupted chromatin. Either of these two possible outcomes should result in an alteration in apparent binding affinity at categories G1 and G2 when compared to category G3. To test this hypothesis, we asked whether increasing thermodynamic entropy would differentially impact category G3 as opposed to categories G1 and G2. Accordingly, we briefly incubated nuclei at 37o C or on ice in physiological salt prior to cross linking and performed ChIP-seq. We observed that binding events in category G3, and only G3, were significantly diminished by increasing temperature (Fig. 2f, g, Additional file 1: Figure S2K, S2L). These data suggest that the on and off rate constants at category G3 differ from those at categories G1 and G2. This scenario is consistent with a model wherein category G3 represents GATA3 bound at its recognition site in the context of canonical nucleosomal architecture whereas categories G1 and G2 fundamentally differ in the nature of the substrate. Utilizing our previous ChIP-seq data  and ENCODE DNase-seq data , we confirmed the presence of a category of G3-like GATA3 bound sites in luminal breast cancer cells (T47D and MCF7 cells) (Additional file 1: Figure S2M, S2N). Of the GATA3 peaks, 51 % (in the case of MCF7) and 37 % (in the case of T47D) did not overlap with DNase-seq peaks in these cells, suggesting the existence of loci that are bound by GATA3 but not accessible, like category G3 in the MDA-231 reprogramming system, is unlikely to be an artefact of expression out of normal context.
The cellular data above suggested that GATA3 may bind to nucleosomes within regions of closed chromatin. To test this hypothesis, we first generated a model structure of a GATA3-bound nucleosome by utilizing crystal structure information [30, 47] for the nucleosome and for the GATA3 DNA binding domain (DBD) (Fig. 3d). In the resulting model, the GATA3-DNA complex was aligned by the ‘GATA’ sequence found within the 601-nucleosome positioning sequence . An appreciable structural barrier to binding was not observed in this model. Therefore, we reconstituted nucleosomes on 601-based DNAs with or without the consensus GATA motif using recombinant human histones (Fig. 3d, e). We also purified the DNA binding domain of GATA3 from engineered bacteria, and performed gel shift analyses (Fig. 3f, g). In the presence of its binding motif, GATA3 displayed discrete shifted bands, suggesting specific binding to the nucleosomal target (Fig. 3g, lanes 1 to 3). However, in the absence of the motif, the binding affinity was significantly reduced and appeared to be non-specific (Fig. 3g, lanes 5 to 7). In both cases, the use of excess amounts of protein resulted in non-specific nucleosome binding (Fig. 3g, lanes 4 and 8). We thus concluded that GATA3 can interact productively with its target recognition sequence within the context of a nucleosome, consistent with function as a pioneer.
Transactivation domain is required for GATA3-mediated chromatin remodeling
GATA3 interacts with BRG1, a chromatin remodeling factor in solution and on chromatin
The inability of mutant GATA3 to remodel chromatin suggested that either the transactivation domain has chromatin remodeling activity similar to the carboxyl-terminal domain of FOXA1 , or that this region functions as a platform to recruit co-factors. Enrichment analysis of genes assigned to the G1 group identified the BRG1/SMARCA4 chromatin remodeling factor as a candidate co-factor (Additional file 1: Table S4). BRG1 is also known to function with the other GATA family proteins such as GATA1 and GATA4 [23, 50, 51] and is characterized by its ability to evict nucleosomes [52–54]. Therefore we asked whether BRG1 and GATA3 had biochemical interactions. We prepared conventional nuclear extracts from MDA-MB-231 cells expressing GATA3. Co-immunoprecipitation analysis indicated solution interactions between GATA3 and BRG1 in the MDA-MB-231 context (Fig. 6a). This interaction was further confirmed by co-immunoprecipitation of the endogenous proteins from T47D cells (Fig. 6b). Accordingly, we explored BRG1 accumulation at GATA3 binding sites in our MDA-MB-231 system using ChIP-seq for endogenous BRG1. We observed consistent enrichment of BRG1 at GATA3 binding sites in cells expressing wild-type, but not mutant, GATA3 at loci where GATA3 functions as a pioneer and evicts nucleosomes (Fig. 6c-e; Additional file 1: Figure S6E, S6F). These data demonstrate that BRG1 recruitment correlates with GATA3-dependent nucleosome eviction in a manner dependent on local context and on the GATA3 amino terminal activation domain.
Information content embedded in eukaryotic chromosomes modulates transcriptional output, influencing the transcriptional patterns integral to biological processes. A critical unit of this information is the TF-DNA complex, which can modulate the gene expression network at multiple levels. TFs receive instruction from local DNA sequence, reading the DNA backbone and the functional groups protruding into the major and minor grooves to recognize specific binding sites in genomic DNA. In eukaryotic nuclei, this paradigm is complicated by the packaging of DNA into chromatin which regulates access to the chemical moieties on DNA that drive productive TF-DNA recognition events. Integral to the elaboration of cell type and specialized function is a requirement for certain TFs to penetrate previously inaccessible chromatin leading to new patterns of transcriptional output.
In this report, we provide detailed analysis of GATA3-mediated chromatin reprogramming in the context of MET in a system relevant to breast cancer; we observed several key outcomes. At approximately 25 % of bound loci (category G1), GATA3 has the capacity to pioneer new binding sites in previously inaccessible chromatin that features high nucleosome occupancy over the GATA motif prior to GATA3 expression. At approximately 45 % of bound loci (category G2), GATA3 binds to a motif located within accessible chromatin where it appears to modestly increase the width of the accessible region. Surprisingly, we observed another 25 % of bound loci (category G3) where GATA3 occupancy is clearly detected by ChIP-seq, but where the local chromatin architecture remains refractory to transposition following GATA3 binding. Biochemical experiments suggested that protein-DNA interactions at the latter category of bound loci – where GATA3 binds to inaccessible chromatin and fails to remodel local architecture to a transposition-permissible state – were fundamentally different than the other two major categories. GATA3 could be functioning to alter chromatin and transcriptional status at both category G1 and G2, however category G2 sites were already open and did not exhibit drastic differences in chromatin architecture following expression of GATA3. Because of this limitation, we focused on the functional properties of binding events in category G1 and G3. We interpret the data to support a model for GATA3 interaction with chromatin. We envision that the protein scans the chromatin fiber for consensus binding motifs where the chemical moieties involved in productive interaction are positioned relative to the nucleosome surface in a manner providing access for interrogation by protein. Recognition of a productive binding site within a nucleosomal context engenders formation of a GATA3/nucleosome complex, a species we are able to recapitulate in vitro (Fig. 3d-g). At a subset of such sites, GATA3 is able to initiate nucleosome eviction over the center of the binding site with subsequent generation of flanking arrays of positioned nucleosomes bearing the hallmarks of an active enhancer (Fig. 2). Such sites are characterized by the GATA3-dependent generation of chromatin accessible to the structural probes sampled here, and by acquisition of active enhancer histone marks. At other sites, GATA3 is able to productively bind as measured by ChIP-Seq. In contrast to category G1 loci, these category G3 sites are characterized both before and after GATA3 binding by a nucleosome array which slightly slide relative to DNA sequence following expression of GATA3. These sites exhibit no evidence of nucleosome eviction over the binding site and do not generate accessible chromatin. We propose that association of co-factors to a GATA3/nucleosome complex differs in a context-dependent manner and that the nature of specific co-factors fundamentally alters outcome. At pioneering loci, the action of co-factors including the SWI/SNF complex alters the protein/substrate interaction, leading to nucleosome eviction and generation of accessible chromatin (Fig. 2e, f). At category G3, other co-factors lead to nucleosome sliding, repositioning an array of nucleosomes relative to DNA sequence without generation of accessibility. In both instances, alterations in local chromatin are dependent on the presence of the activation domain.
The model proposed here envisions two distinct steps to pioneering activity by GATA3 – binding followed by local remodeling – and differs from that explaining the pioneering activity of FOXA1. Elegant biochemical and cellular data suggest that FOXA1 has the intrinsic capacity to elicit local perturbation of chromatin structure without the involvement of additional co-factors . Pioneering activity by GATA3, on the other hand, correlates with the capacity to recruit the mammalian SWI/SNF chromatin remodeling complex and requires the presence of the transcriptional activation domain (Figs. 5 and 6). While local remodeling by FOXA1 appears inextricably linked to productive binding, GATA3 appears to require co-factors for this function. A critical question arises as to why GATA3 binding leads to distinct outcomes at different loci despite the observation of a biochemical interaction between GATA3 and BRG1. It is possible that this interaction alone is not sufficient to drive the initiation of nucleosome eviction at all loci; it is also conceivable that local architectural constraints influence interaction with other co-factors (Additional file 3: Table S2). Whatever the proximate cause, the differences in local nucleosome architecture around GATA3 binding sites prior to GATA3 expression at loci that differ in outcome (nucleosome sliding versus eviction) (Fig. 3b, c) suggest architectural constraints may distinguish the two types of events. Architectural influences of DNA, and by extension chromatin, on outcomes of TF mediated events have been elegantly documented in the case of the glucocorticoid receptor .
GATA3 has been identified as one of the most frequently mutated genes in breast cancer [35, 56]. GATA3 mutations accumulate in a striking pattern at the C-terminal region of the protein including the second zinc finger [56, 57]. Mutations in the DNA-binding domain have been shown to affect chromatin binding activity and motif recognition  and are likely to lead to alterations in genomic distribution. The outcomes of mutations at the far carboxyl terminus of GATA3, an abundant class of cancer-specific mutants, are currently understudied. Because most cancer-derived mutants have an intact N-terminal domain, their abnormal distribution may induce inappropriate pioneering activity leading to reprogrammed gene expression. It is also conceivable that the cancer-specific mutations in GATA3 alter co-factor recruitment and interaction (either positively or negatively) leading to novel or impaired pioneering functions. In luminal breast cancer cells, GATA3 is known to functionally interact with the TFs FOXA1 and estrogen receptor-α [28, 58], leading to the potential for GATA3 mutations to alter the transcriptional network of these factors. Given the capacity of GATA3 to catalyze cellular reprogramming relevant to cancer biology, characterization of the impact of cancer-specific mutations in this pioneering TF represents a field or enquiry with major clinical ramifications.
Our data provide new experimental evidence that GATA3 functions as a pioneer factor during cellular reprogramming using the MET of breast cancer cells as a model. Dissection of chromatin features during this process reveals key mechanistic features required to generate active chromatin structure: binding to a nucleosomal locus in the appropriate context and transactivation domain dependent co-factor recruitment. These findings provide novel insights into the mechanism of de novo enhancer establishment by a pioneer TF.
Generation of stable cell line
Lentiviruses were produced by transient transfection of 293T cells with pHAGE, psPAX2, and pMD2.G plasmids. psPAX2 and pMD2.G were gifts from Didier Trono (Addgene plasmid #12260, #12259). The virus-containing medium and polybrene (at a final concentration of 8 μg/mL) were added to MDA-MB-231 cells. Infected cells were selected with 0.5 μg/mL puromycin. Only early passage cells (passages 3 to 7) were used in all experiments. Cells expressing either GFP or the tag alone (including the TEV protease recognition site) were used as a control cell line.
Nucleosome binding assay
The reconstituted nucleosomes and GATA3 DBD were incubated at room temperature for 30 min in the reaction buffer containing 30 mM Tris-HCl (pH 7.5), 120 mM NaCl, 25 mM KCl, 4 % glycerol, 1.5 % Ficoll 400, 0.4 mM zinc sulphate, and 0.1 mM EDTA. The complex and free nucleosomes were separated on 6 % native-polyacrylamide gel, and visualized by ethidium bromide staining.
ATAC-seq sample preparation and data analysis
The sequence library for ATAC-seq was prepared as previously described  with several modifications (see Additional file 1: supplementary methods). Reads were mapped to hg19 genome using Bowtie 0.12.8  with the same parameter used in . Only non-duplicate reads were used for the subsequent analysis. The offset parameters suggested by Buenrostro et al. were applied to each mapped read . Reads from three biological replicates were then merged for use in all subsequent analysis.
Native mononucleosomes were prepared as previously described  with several modifications (see Additional file 1: supplementary methods). Briefly, cells were treated with a hypotonic buffer for 15 min on ice. The isolated chromatin was digested with 4 U of MNase at 37 °C for 4 min. Digested chromatin was homogenized by passing through a 27G ½ needle (8 strokes), and the debris was removed by centrifugation. DNA was purified using a spin column. Mononucleosomes were separated by agarose gel electrophoresis, and purified by a QIAGEN gel extraction kit. Sequencing libraries were prepared by NEXTflex Rapid DNA-Seq kit. Reads were mapped to hg19 genome, and duplicate reads were removed using MarkDuplicates.jar from picard-tools-1.107 package (http://broadinstitute.github.io/picard/). All paired-end reads were converted to a single fragment for the metagene plot analysis.
Chromatin immunoprecipitation (ChIP) for histones and BRG1
ChIP experiments were performed as previously described  with the following modifications. Cells were fixed with 1 % formaldehyde either at room temperature for 10 min (histone ChIP) or at 37 °C for 20 min (BRG1 ChIP), and quenched with glycine. For BRG1 ChIP, fixed cells were further treated with a hypotonic buffer. Cell pellets were resuspended with lysis buffer A (20 mM Tris-HCl pH 8.0, 2 mM EDTA, 0.5 mM EGTA, 0.5 mM PMSF, 5 mM sodium butyrate, 1 % SDS and protease inhibitor cocktail) for histone ChIP or lysis buffer B (same as buffer A, except 0.1 % SDS) for BRG1 ChIP. Chromatin was fragmented by sonication with Bioruptor (histone ChIP), or Covaris S220 (BRG1 ChIP), then diluted to adjust the SDS concentration to 0.1 % or 0.05 %. Immunoprecipitation was performed with indicated antibodies (Additional file 1: Table S5). Dynabeads Protein A and G were used to capture antibodies, and beads were washed with 1 mL of low salt (20 mM Tris-HCl pH 8.0, 150 mM NaCl, 2 mM EDTA, 1 % Triton X-100, 0.1 % SDS), high salt (same as low salt buffer, except 500 mM NaCl), and LiCl buffer (Tris-HCl pH 8.0, 250 mM LiCl, 2 mM EDTA, 1 % NP-40, 1 % (wt/vol) sodium deoxycholate). Eluted DNA was treated with RNase at 65 °C, followed by the incubation with proteinase K. DNA was purified by AMPure XP. Histone and BRG1 ChIP experiments were conducted with two and three biological replicates, respectively (Additional file 1: Table S5).
ChIP for GATA3
Cells were treated with CSK buffer 2 (same as described above, except containing 1 mM EGTA, 2 mM EDTA, 1 mM PMSF, and protease inhibitor cocktail) for 5 min on ice or at 37 °C. After centrifugation at 100 g for 5 min at 4 °C, nuclear pellets were washed with CSK buffer 3 (same as CSK buffer1, except that Triton X-100 was omitted). After centrifugation, nuclei were resuspended with CSK buffer 3, then an equal volume of 2 % formaldehyde containing CSK buffer 3 was added to the solution, followed by the incubation at 4 °C for 15 min. The fixation was quenched by addition of glycine, and the fixed nuclei were washed twice with CSK buffer 3. The nuclei were lysed with 0.5 % SDS containing lysis buffer A, and chromatin was fragmented by sonication with Bioruptor. The fragmented chromatin was diluted with nine-fold volume of IP buffer (20 mM Tris-HCl pH 8.0, 150 mM NaCl, 2 mM EDTA, 10 % glycerol, 0.5 % Triton X-100, 0.5 mM PMSF, 5 mM sodium butyrate, and protease inhibitor cocktail). Immunoprecipitation was performed as described above with two biological replicates (Additional file 1: Table S5). Homemade GATA3 antibody (prepared in Wade laboratory) was used for ChIP experiments , and Ty1 antibody was prepared in Stunnenberg laboratory .
Reads were filtered based on a mean base quality score >20, and mapped to hg19 genome using Bowtie 0.12.8 . In order to minimize PCR amplification bias, duplicate reads were removed using MarkDuplicates.jar from picard-tools-1.107 package. Reads from the biological replicates were merged for use in all subsequent analysis. All paired-end reads were converted to a single fragment for the metagene plot analysis and visualization on the UCSC Genome Browser. After normalization to the indicated counts (Additional file 1: Table S6), bigWig files were generated to visualize the genomic coverage.
Peak call and motif analysis
Peak calling and motif analyses were performed using HOMER v4.1 with default parameters . Motif length 8, 10, and 12 bases were used for de novo motif analysis. Based on the distance from the midpoint of a peak to the nearest RefSeq transcription start site, we defined all peaks as either TSS-proximal (distance: less than or equal to 1 kb) or distal (distance: more than 1 kb) peaks.
GATA3 peak classification based on ATAC-seq
EdgeR  was used to divide the GATA3 peaks into three groups (signal stronger: G1, signal conserved, and signal reduced: G3) based on the signal changes of ATAC-seq before and after GATA3-induced epithelial transition. After normalization to 30 million total reads, ATAC-seq read counts in the 200 bp window centered on GATA3 peaks were collected. The following parameters were used for the peak classification: FDR <0.1 and fold change >2. The signal conserved group was further subclassified into signal positive (G2) and negative (G3) groups based on a threshold of 20 counts in the +/- 1 kb region flanking the peak midpoint. In the heatmaps (Figs. 2a, 5a, and 6d, Additional file 1: Figure S5A), the peaks are ordered by decreasing weighted ATAC signal within each of the four peak classes.
Box plots were generated by with bwplot (lattice R package; for randomly selected GATA motif containing genomic regions) or Prism 6 (GraphPad Software). Mann-Whitney test was used to compare data groups. Asterisk indicates P value <0.0001.
Availability of supporting data
The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE72141 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE72141).
Ethical approval was not needed for this study.
We gratefully acknowledge Jason Williams, Jeff Tucker, Christine Malone, Robert Petrovich, and Jason Malphurs for outstanding technical assistance during the course of this work, and Guang Hu for expression vectors. This manuscript was substantially improved by critical comments from Paula Vertino, Guang Hu, and Raja Jothi. This work was supported, in part, by the Intramural Research Program of the National Institute of Environmental Health Sciences, NIH (ES101965 to PAW; ES071006 to TKA; and ES043010 to LP), by MEXT KAKENHI Grant Number 25116002 (to HK) and by a Fellowship Award from the Japan Society for the Promotion of Science (to MT).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Yusuf D, Butland SL, Swanson MI, Bolotin E, Ticoll A, Cheung WA, et al. The transcription factor encyclopedia. Genome Biol. 2012;13:R24.PubMed CentralView ArticlePubMedGoogle Scholar
- The ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.PubMed CentralView ArticleGoogle Scholar
- Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10:252–63.View ArticlePubMedGoogle Scholar
- Geertz M, Maerkl SJ. Experimental strategies for studying transcription factor-DNA binding specificities. Brief Funct Genomics. 2010;9:362–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Zaret KS, Carroll JS. Pioneer transcription factors: establishing competence for gene expression. Genes Dev. 2011;25:2227–41.PubMed CentralView ArticlePubMedGoogle Scholar
- Shlyueva D, Stampfel G, Stark A. Transcriptional enhancers: from properties to genome-wide predictions. Nat Rev Genet. 2014;15:272–86.View ArticlePubMedGoogle Scholar
- Adams CC, Workman JL. Binding of disparate transcriptional activators to nucleosomal DNA is inherently cooperative. Mol Cell Biol. 1995;15:1405–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Cirillo LA, Zaret KS. An early developmental transcription factor complex that is more stable on nucleosome core particles than on free DNA. Mol Cell. 1999;4:961–9.View ArticlePubMedGoogle Scholar
- Kornberg RD, Lorch Y. Irresistible force meets immovable object: transcription and the nucleosome. Cell. 1991;67:833–6.View ArticlePubMedGoogle Scholar
- Iwafuchi-Doi M, Zaret KS. Pioneer transcription factors in cell reprogramming. Genes Dev. 2014;28:2679–92.PubMed CentralView ArticlePubMedGoogle Scholar
- Cirillo LA, Lin FR, Cuesta I, Friedman D, Jarnik M, Zaret KS. Opening of compacted chromatin by early developmental transcription factors HNF3 (FoxA) and GATA-4. Mol Cell. 2002;9:279–89.View ArticlePubMedGoogle Scholar
- Cirillo LA, McPherson CE, Bossard P, Stevens K, Cherian S, Shim EY, et al. Binding of the winged-helix transcription factor HNF3 to a linker histone site on the nucleosome. EMBO J. 1998;17:244–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Gualdi R, Bossard P, Zheng M, Hamada Y, Coleman JR, Zaret KS. Hepatic specification of the gut endoderm in vitro: cell signaling and transcriptional control. Genes Dev. 1996;10:1670–82.View ArticlePubMedGoogle Scholar
- Bernardo GM, Lozada KL, Miedler JD, Harburg G, Hewitt SC, Mosley JD, et al. FOXA1 is an essential determinant of ERalpha expression and mammary ductal morphogenesis. Development. 2010;137:2045–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Besnard V, Wert SE, Kaestner KH, Whitsett JA. Stage-specific regulation of respiratory epithelial cell differentiation by Foxa1. Am J Physiol Lung Cell Mol Physiol. 2005;289:L750–759.View ArticlePubMedGoogle Scholar
- Fakhouri TH, Stevenson J, Chisholm AD, Mango SE. Dynamic chromatin organization during foregut development mediated by the organ selector gene PHA-4/FoxA. PLoS Genet. 2010;6.Google Scholar
- Lee CS, Friedman JR, Fulmer JT, Kaestner KH. The initiation of liver development is dependent on Foxa transcription factors. Nature. 2005;435:944–7.View ArticlePubMedGoogle Scholar
- Soufi A, Donahue G, Zaret KS. Facilitators and impediments of the pluripotency reprogramming factors’ initial engagement with the genome. Cell. 2012;151:994–1004.PubMed CentralView ArticlePubMedGoogle Scholar
- Soufi A, Garcia MF, Jaroszewicz A, Osman N, Pellegrini M, Zaret KS. Pioneer transcription factors target partial DNA motifs on nucleosomes to initiate reprogramming. Cell. 2015;161:555–68.View ArticlePubMedGoogle Scholar
- Clark KL, Halay ED, Lai E, Burley SK. Co-crystal structure of the HNF-3/fork head DNA-recognition motif resembles histone H5. Nature. 1993;364:412–20.View ArticlePubMedGoogle Scholar
- Ramakrishnan V, Finch JT, Graziano V, Lee PL, Sweet RM. Crystal structure of globular domain of histone H5 and its implications for nucleosome binding. Nature. 1993;362:219–23.View ArticlePubMedGoogle Scholar
- Ravasi T, Huber T, Zavolan M, Forrest A, Gaasterland T, Grimmond S, et al. Systematic characterization of the zinc-finger-containing proteins in the mouse transcriptome. Genome Res. 2003;13:1430–42.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu G, Schones DE, Cui K, Ybarra R, Northrup D, Tang Q, et al. Regulation of nucleosome landscape and transcription factor targeting at tissue-specific enhancers by BRG1. Genome Res. 2011;21:1650–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Kadauke S, Udugama MI, Pawlicki JM, Achtman JC, Jain DP, Cheng Y, et al. Tissue-specific mitotic bookmarking by hematopoietic transcription factor GATA1. Cell. 2012;150:725–37.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee HJ, Takemoto N, Kurata H, Kamogawa Y, Miyatake S, O’Garra A, et al. GATA-3 induces T helper cell type 2 (Th2) cytokine expression and chromatin remodeling in committed Th1 cells. J Exp Med. 2000;192:105–15.PubMed CentralView ArticlePubMedGoogle Scholar
- Ouyang W, Lohning M, Gao Z, Assenmacher M, Ranganath S, Radbruch A, et al. Stat6-independent GATA-3 autoactivation directs IL-4-independent Th2 development and commitment. Immunity. 2000;12:27–37.View ArticlePubMedGoogle Scholar
- Takemoto N, Kamogawa Y, Jun Lee H, Kurata H, Arai KI, O’Garra A, et al. Cutting edge: chromatin remodeling at the IL-4/IL-13 intergenic regulatory region for Th2-specific cytokine gene cluster. J Immunol. 2000;165:6687–91.View ArticlePubMedGoogle Scholar
- Theodorou V, Stark R, Menon S, Carroll JS. GATA3 acts upstream of FOXA1 in mediating ESR1 binding by shaping enhancer accessibility. Genome Res. 2013;23:12–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei G, Abraham BJ, Yagi R, Jothi R, Cui K, Sharma S, et al. Genome-wide analyses of transcription factor GATA3-mediated gene regulation in distinct T cell types. Immunity. 2011;35:299–311.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Bates DL, Dey R, Chen PH, Machado AC, Laird-Offringa IA, et al. DNA binding by GATA transcription factor suggests mechanisms of DNA looping and long-range gene regulation. Cell Rep. 2012;2:1197–206.PubMed CentralView ArticlePubMedGoogle Scholar
- Ng YK, George KM, Engel JD, Linzer DI. GATA factor activity is required for the trophoblast-specific transcriptional regulation of the mouse placental lactogen I gene. Development. 1994;120:3257–66.PubMedGoogle Scholar
- Ho IC, Vorhees P, Marin N, Oakley BK, Tsai SF, Orkin SH, et al. Human GATA-3: a lineage-restricted transcription factor that regulates the expression of the T cell receptor alpha gene. EMBO J. 1991;10:1187–92.PubMed CentralPubMedGoogle Scholar
- Ko LJ, Yamamoto M, Leonard MW, George KM, Ting P, Engel JD. Murine and human T-lymphocyte GATA-3 factors mediate transcription through a cis-regulatory element within the human T-cell receptor delta gene enhancer. Mol Cell Biol. 1991;11:2778–84.PubMed CentralView ArticlePubMedGoogle Scholar
- Kouros-Mehr H, Slorach EM, Sternlicht MD, Werb Z. GATA-3 maintains the differentiation of the luminal cell fate in the mammary gland. Cell. 2006;127:1041–55.PubMed CentralView ArticlePubMedGoogle Scholar
- The Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.PubMed CentralView ArticleGoogle Scholar
- Chou J, Provot S, Werb Z. GATA3 in development and cancer differentiation: cells GATA have it! J Cell Physiol. 2010;222:42–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Chou J, Lin JH, Brenot A, Kim JW, Provot S, Werb Z. GATA3 suppresses metastasis and modulates the tumour microenvironment by regulating microRNA-29b expression. Nat Cell Biol. 2013;15:201–13.PubMed CentralView ArticlePubMedGoogle Scholar
- Chu IM, Michalowski AM, Hoenerhoff M, Szauter KM, Luger D, Sato M, et al. GATA3 inhibits lysyl oxidase-mediated metastases of human basal triple-negative breast cancer cells. Oncogene. 2012;31:2017–27.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan W, Cao QJ, Arenas RB, Bentley B, Shao R. GATA3 inhibits breast cancer metastasis through the reversal of epithelial-mesenchymal transition. J Biol Chem. 2010;285:14042–51.PubMed CentralView ArticlePubMedGoogle Scholar
- Montserrat N, Nivet E, Sancho-Martinez I, Hishida T, Kumar S, Miquel L, et al. Reprogramming of human fibroblasts to pluripotency with lineage specifiers. Cell Stem Cell. 2013;13:341–50.View ArticlePubMedGoogle Scholar
- Benchetrit H, Herman S, van Wietmarschen N, Wu T, Makedonski K, Maoz N, Yom Tov N, Stave D, Lasry R, Zayat V, et al. Extensive Nuclear Reprogramming Underlies Lineage Conversion into Functional Trophoblast Stem-like Cells. Cell Stem Cell. 2015;17:543–56.View ArticlePubMedGoogle Scholar
- Kubaczka C, Senner CE, Cierlitza M, Arauzo-Bravo MJ, Kuckenberg P, Peitz M, Hemberger M, Schorle H. Direct Induction of Trophoblast Stem Cells from Murine Fibroblasts. Cell Stem Cell. 2015;17:557–68.View ArticlePubMedGoogle Scholar
- Adomas AB, Grimm SA, Malone C, Takaku M, Sims JK, Wade PA. Breast tumor specific mutation in GATA3 affects physiological mechanisms regulating transcription factor turnover. BMC Cancer. 2014;14:278.PubMed CentralView ArticlePubMedGoogle Scholar
- Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10:1213–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsang AP, Visvader JE, Turner CA, Fujiwara Y, Yu C, Weiss MJ, et al. FOG, a multitype zinc finger protein, acts as a cofactor for transcription factor GATA-1 in erythroid and megakaryocytic differentiation. Cell. 1997;90:109–19.View ArticlePubMedGoogle Scholar
- Zhou M, Ouyang W, Gong Q, Katz SG, White JM, Orkin SH, et al. Friend of GATA-1 represses GATA-3-dependent activity in CD4+ T cells. J Exp Med. 2001;194:1461–71.PubMed CentralView ArticlePubMedGoogle Scholar
- Vasudevan D, Chua EY, Davey CA. Crystal structures of nucleosome core particles containing the ‘601’ strong positioning sequence. J Mol Biol. 2010;403:1–10.View ArticlePubMedGoogle Scholar
- Lowary PT, Widom J. New DNA sequence rules for high affinity binding to histone octamer and sequence-directed nucleosome positioning. J Mol Biol. 1998;276:19–42.View ArticlePubMedGoogle Scholar
- Yang Z, Gu L, Romeo PH, Bories D, Motohashi H, Yamamoto M, et al. Human GATA-3 trans-activation, DNA-binding, and nuclear localization activities are organized into distinct structural domains. Mol Cell Biol. 1994;14:2201–12.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim SI, Bresnick EH, Bultman SJ. BRG1 directly regulates nucleosome structure and chromatin looping of the alpha globin locus to activate transcription. Nucleic Acids Res. 2009;37:6019–27.PubMed CentralView ArticlePubMedGoogle Scholar
- Lickert H, Takeuchi JK, Von Both I, Walls JR, McAuliffe F, Adamson SL, et al. Baf60c is essential for function of BAF chromatin remodelling complexes in heart development. Nature. 2004;432:107–12.View ArticlePubMedGoogle Scholar
- Gutierrez JL, Chandy M, Carrozza MJ, Workman JL. Activation domains drive nucleosome eviction by SWI/SNF. EMBO J. 2007;26:730–40.PubMed CentralView ArticlePubMedGoogle Scholar
- Lorch Y, Zhang M, Kornberg RD. Histone octamer transfer by a chromatin-remodeling complex. Cell. 1999;96:389–92.View ArticlePubMedGoogle Scholar
- Phelan ML, Schnitzler GR, Kingston RE. Octamer transfer and creation of stably remodeled nucleosomes by human SWI-SNF and its isolated ATPases. Mol Cell Biol. 2000;20:6380–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Watson LC, Kuchenbecker KM, Schiller BJ, Gross JD, Pufall MA, Yamamoto KR. The glucocorticoid receptor dimer interface allosterically transmits sequence-specific DNA signals. Nat Struct Mol Biol. 2013;20:876–83.PubMed CentralView ArticlePubMedGoogle Scholar
- Usary J, Llaca V, Karaca G, Presswala S, Karaca M, He X, et al. Mutation of GATA3 in human breast tumors. Oncogene. 2004;23:7669–78.View ArticlePubMedGoogle Scholar
- Takaku M, Grimm SA, Wade PA. GATA3 in Breast Cancer: Tumor Suppressor or Oncogene? Gene Expr. 2015;16:163–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong SL, Li G, Loh SL, Sung WK, Liu ET. Cellular reprogramming by the conjoint action of ERalpha, FOXA1, and GATA3 to a ligand-inducible growth state. Mol Syst Biol. 2011;7:526.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.PubMed CentralView ArticlePubMedGoogle Scholar
- Cui K, Zhao K. Genome-wide approaches to determining nucleosome occupancy in metazoans using MNase-Seq. Methods Mol Biol. 2012;833:413–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Shimbo T, Du Y, Grimm SA, Dhasarathy A, Mav D, Shah RR, et al. MBD3 localizes at promoters, gene bodies and enhancers of active genes. PLoS Genet. 2013;9:e1004028.PubMed CentralView ArticlePubMedGoogle Scholar
- Menafra R, Brinkman AB, Matarese F, Franci G, Bartels SJ, Nguyen L, et al. Genome-wide binding of MBD2 reveals strong preference for highly methylated loci. PLoS One. 2014;9:e99603.PubMed CentralView ArticlePubMedGoogle Scholar
- Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.PubMed CentralView ArticlePubMedGoogle Scholar