Dynamic binding of transcription factors (TFs) to regulatory elements controls transcriptional states throughout organism development. Epigenetics modifications, such as DNA methylation mostly within cytosine-guanine dinucleotides (CpGs), have the potential to modulate TF binding to DNA. Although DNA methylation has long been thought to repress TF binding, a more recent model proposes that TF binding can also inhibit DNA methylation. Here, we review the possible scenarios by which DNA methylation and TF binding affect each other. Further in vivo experiments will be required to generalize these models.
Multicellular organisms establish and maintain different transcriptional states in diverse cell types through dynamic and specific regulation of gene expression. This regulation is mediated by transcription factors (TFs) binding to regulatory elements . TFs are often defined as any protein able to bind DNA and affect gene expression. They are composed of a DNA-binding domain, a trans-activating domain, which binds other proteins and mediates activator or repressor functions and an optional signal-sensing domain, which regulates their activity (Figure 1A). TFs bind to DNA, through the recognition of short specific DNA sequence motifs (Figure 1B). In the past years, much effort has been invested into the identification of TF motifs, using in vitro experiments such as electrophoretic mobility shift assay (EMSA) , high-throughput systematic evolution of ligands by exponential enrichment (HT-SELEX)  and protein binding microarrays (PBM) . As a result, thousands of TF consensus sequence motifs (represented as sequence logos; Figure 1B) have been stored in databases such as JASPAR  and have been used as valuable input for further functional studies. However, while hundreds of thousands of occurrences of specific TF motifs are found within genomes, only a few thousands are bound in vivo in a context-specific manner, as assessed by chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq). Therefore, the specificity of TF binding must be controlled by other means. Combinatorial binding of TFs is a hallmark of regulatory regions, where TF binding relies on the presence of its motif, as well as motifs of other co-expressed TFs  (Figure 1C). More recently, the 3D structure of DNA or ‘DNA shape’, characterized by physical features of DNA major or minor grooves, has been thought to affect TF binding  (Figure 1D).
TF binding to DNA
The physical access of TFs to DNA can further be modulated by epigenetic regulation such as nucleosome positioning, histone modifications or DNA methylation  (Figure 1E). Early on, DNA methylation of cytosines, mostly within cytosine-guanine dinucleotides (CpGs) (Box 1), emerged as a mechanism to modulate TF affinity. DNA methylation was described as a mark repressing transcription, where the presence of DNA methylation at CpG-rich gene promoters, called CpG islands, would block TF binding leading to gene silencing . With the emergence of high-throughput sequencing and the development of new methods to better profile TF binding and DNA methylation, our views on the sensitivity of TFs to DNA methylation have started to change. DNA methylation was found to be abundant throughout the genome except at active regulatory regions bound by TFs, such as gene promoters and distal enhancers [19–21]. In parallel, pioneer TFs have been described as able to bind non-accessible chromatin and recruit remodeling complexes, which can displace or remove nucleosomes, thus facilitating the context-specific binding of TFs . Together with the observation of a strong anti-correlation between TF binding and DNA methylation patterns genome wide, a new actively discussed hypothesis has emerged where DNA methylation could be a consequence rather than a cause of TF binding and transcriptional activity .
DNA methylation is a chemical modification conserved across evolution from archeas to higher vertebrates such as human, mouse, zebrafish, plants or sea squirts although sporadic or even absent from invertebrate species such as fruit flies, worms and yeast . DNA methylation affects cytosine residues, where a methyl group is added to the fifth carbon atom, thus forming a 5-methylcytosine (5mC). While in most vertebrates, DNA methylation occurs mainly in a CpG context, in other species such as plants, it also occurs in a CHG or CHH context .
DNA methylation is mediated by the DNA methyltransferases (DNMTs): DNMT3A, DNMT3B, and their cofactor DNMT3L, establish 5mC de novo, while DNMT1, with its co-cofactor UHRF1, maintains DNA methylation patterns during replication [10,11]. 5mC is a reversible modification. It can be removed passively during successive rounds of replication through an intermediate state of hemimethylation . Recently, Ten-Eleven Translocation (TET) enzymes have been identified as drivers of active demethylation through an intermediate state of hydroxymethylation (and further formyl- and carboxylmethylation) [13–15].
While the chemical evidence of cytosine methylation was initially discovered in 1948 , the first functional evidence that it could regulate physiological processes emerged in 1975 . DNA methylation is now known to be implicated in regulating gene expression in many biological processes such as development, cellular differentiation and cancer .
In this review, we summarize the progress and evolution of our understanding of the sensitivity of TFs to DNA methylation, which appears to be factor-specific and condition-dependent. We investigate the determinants of TF sensitivity to DNA methylation, the impact of TF binding on DNA methylation patterns and discuss how these principles can be generalized.
Impact of DNA methylation on TF binding
DNA methylation represses TF binding
The sensitivity of TF to DNA methylation started being investigated in vitro in the late 80s and already different TFs appeared to have different sensitivities . Several TFs were identified as sensitive to DNA methylation by EMSA (Figure 2A). Methylation of a CpG (mCpG) central to the MLTF (also called USF) motif prevented binding of the TFs and inhibited the expression of the adenovirus major late promoter, whereas methylation of a CpG 6 base pairs away had no effect . Methylation of the CREB binding sites also resulted in a loss of TF binding and transcriptional activity . Other TFs, such as AP-2, MYC, E2F, NF-kB and ETS were repressed from binding by mCpG within their binding sites [24,27].
Possible scenarios by which DNA methylation could impact TF binding
The zinc finger CTCF is one of the most studied TFs. It was first identified in chickens  and shortly after in humans, as a negative regulator of the c-myc gene  and an activator of the amyloid β protein precursor (AmBP) . It was demonstrated to act as an insulator on the chicken β-globin locus, blocking enhancers from activating distal genes . By dimerizing and tethering loop formation, it acts on the genome architecture and regulates gene expression [32,33].
By investigating imprinted genes, CTCF was found to be sensitive to DNA methylation at the Igf2-H19 locus in mouse. When bound to the unmethylated maternal allele, CTCF acts as an insulator, restricting the action of the downstream enhancer to the H19 gene. Whereas, on the methylated paternal allele, CTCF cannot bind and the enhancer activates the Igf2 gene . It was later shown by point mutations of the four CTCF sites in vivo that impaired CTCF binding led to an increase in local DNA methylation in newborn mice although CTCF binding was not required for establishing an unmethylated state during oogenesis . Additionally, methylation of two CTCF motifs at the myotonic dystrophy type 1 (DM1) locus was also shown to abolish CTCF binding thus altering the expression of the DMPK and SIX5 genes .
Ten years later, the development of high-throughput sequencing approaches to probe DNA methylation patterns (whole-genome bisulfite sequencing) and TF binding (ChIP-seq), enabled to better study CTCF binding genome wide and challenged traditional views. In mouse embryonic stem (ES) cells, CTCF binding sites were found to be mainly located in regions with no or low levels of DNA methylation . However, CTCF did not bind additional sites in absence of DNA methylation (DNMTs triple knockout), except at known imprinted loci, suggesting that CTCF binding was not repressed by DNA methylation genome wide in vivo. Validations using stable insertions of a reporter construct containing the CTCF motif showed that CTCF could bind methylated DNA and led to local demethylation whereas mutated CTCF motifs were not bound and remained methylated . Therefore, in contrary to what was expected, CTCF appeared to be insensitive to DNA methylation genome wide in vivo , which was further confirmed in HCT116 cells . CTCF profiling across 19 different cell types showed that while CTCF binds distinct sites in different cell types, 41% of its variable binding sites are linked to DNA methylation . When looking at CTCF motif occurrences within the genome, 25% do contain CpGs  whereas 45% of those located within CTCF ChIP-seq binding sites do. This highlights the fact that only a fraction of putative CTCF binding sites are potentially affected by DNA methylation. Recent reports investigated the contribution of specific CpGs to binding methylated DNA using structural information and a binding affinity assay (Methyl-Spec-seq) and could identify that it is the methylated cytosine at position 5 in the JASPAR motif that specifically inhibit CTCF binding [39,40].
In summary, although CTCF was originally described as a methylation-sensitive factor at the imprinted Igf2-H19 locus, only a limited set of in vivo CTCF binding sites, presumably harboring CpGs in their motifs, will be sensitive to DNA methylation. Its flexible binding sites might help regulate complex developmental and cellular processes. CTCF thus perfectly illustrates the fact that TF sensitivity to DNA methylation, despite its early discovery, remains an open question. CTCF motif JASPAR ID: MA0139.1.
The development of high-throughput technologies later enabled to test the sensitivity of many more TFs in vitro. A quantitative mass-spectrometry approach identified ZBTB2, JUND, CREB1, ATF7 as preferentially bound to unmethylated cytosines over methylated ones . An approach using methylated binding microarrays found that DNA methylation inhibited binding of the basic leucine zipper (BZIP) TFs CREB, ATF4, JUN, JUND, CEBPD and CEBPG . TFs identified as being repressed from binding by DNA methylation were in good agreement in both studies and earlier ones.
NRF1 was the first TF shown to be sensitive to DNA methylation genome-wide in vivo . Sensitive TFs were identified in mouse embryonic stem (ES) cells by profiling open chromatin regions in presence and absence of DNA methylation (using DNMTs triple knockout (TKO) cells) as they were expected to bind new accessible sites in absence of DNA methylation. A motif analysis identified NRF1 as well as MYC/USF/CREB and GABPA/ETS as candidate methylation-sensitive TFs. NRF1 was then validated in vivo by ChIP-seq, where it could bind many new sites in absence of DNA methylation. Further, validations using stable insertions of a reporter construct containing the NRF1 motif showed that it could only bind its unmethylated motif and was not able to bind when methylated .
More recent in vitro approaches using methylation-sensitive SELEX expanded the catalog of TFs known to be repressed from binding by DNA methylation. The first study investigated 519 TFs using methyl-SELEX and bisulfite-SELEX and found that 23% (117 out of 519) were inhibited by mCpG (called methyl-minus) . Their global analysis could identify TF families that tend to be inhibited by mCpG such as basic helix–loop–helix (BHLH), BZIP and ETS TFs. The majority (96 TFs, 82%) had CpGs in their original motif and validated known sensitive TFs such as MYC, USF, CREB, ATF, AP (JUN, FOS), E2F, ETS although NRF1 was not mentioned in their study. However, profiling of MYC binding in vivo in cells lacking DNA methylation (ChIP-seq in DNMTs TKO) showed that DNA methylation had only a minimal effect on its binding sites .
Another study probed ATF4 sensitivity using methyl-SELEX . Since the motif does not have a prominent CpG, they found that motifs with no CpG showed no preferential binding for methylated or unmethylated DNA, motifs with CpGs in the center were not bound when methylated and motifs with CpGs in the flank were bound when methylated . This confirms that not only an mCpG but also its position within the motif is critical for inhibiting TF binding.
In parallel, an in vitro study in plants identified 234 TFs (72% out of 327 tested) to be inhibited from binding by mCpGs  although plants have a different TF repertoire.
Most studies agree that inhibition by an mCpG might be due to steric hindrance of TF binding [26,44]. More recently, DNA shape has been described as an additional feature of TF binding  and roll and propeller twist DNA shapes have been found to be strongly affected by mCpGs . Additionally, nucleosome positioning and histone modifications are linked to DNA methylation and might also affect TF binding [48,49].
These different studies identified several TFs to be inhibited from binding by mCpGs and therefore sensitive to DNA methylation, suggesting a widespread mechanism (Figure 2A). Of note, most of these observations came from in vitro studies, and few in vivo studies confirm their sensitivity. Therefore, the functional impact of this sensitivity on a genome-wide scale remains to be further investigated.
DNA methylation promotes TF binding
In parallel to the discovery of TFs that are inhibited by mCpGs, proteins that recognize mCpGs specifically through a methyl-binding domain (MBD) were identified: MeCp2, MBD1, MBD2 and MBD4 [50–52]. However, this recognition can be considered as independent of the underlying DNA sequence, as opposed to TFs that recognize specific DNA sequence motifs .
Later, both individual and high-throughput in vitro studies have identified sequence-specific TFs that bind mCpGs (Figure 2B). They are also described as sensitive to DNA methylation since they require an mCpG to bind as opposed to the sensitive TFs that are inhibited by mCpGs. A quantitative mass-spectrometry approach identified 19 proteins as binding preferentially to mCpGs over unmethylated ones (MeCP2, MBD1, MBD4, UHRF1, RFX1/5, ZFHX3,KLF2/4/5) although it does not account for sequence specificity . A competition assay on methylated protein microarrays identified 41 TFs and 6 cofactors (3% of 1321 TFs and 210 cofactors tested) to preferentially bind motifs with mCpGs although in presence of a ten-fold excess of unmethylated sequences . However, the majority recognized several distinct motifs and only 22 recognized fewer than three different motifs. Eight out of eleven were validated by EMSA (including the TFs ARNT2, DIDO1, MEF2A and HOXA9) implying a 27% false-positive rate . NRF1 was found to bind both states, but more strongly to methylated than unmethylated motifs , although it was later shown to be inhibited from binding by mCpG in vivo . A study in plants identified 14 TFs (4.3% out of 327 tested), which preferentially bind methylated motifs . More recently, an approach using methyl-SELEX and bisulfite-SELEX found that 34% (175 out of 519) of the tested TFs could bind mCpGs (called methyl-plus) such as KAISO/ZBTB33, CEBPB/E/G, KLF, OCT4, HOX, PAX or SP1 . However, only 49% of these (85 out of 175) had a CpG in their canonical motif whereas the others recognized a weaker methylated site . Another approach called Methyl-Spec-seq, measuring the effect of mCpGs on TF binding affinity at every position within a binding site, could quantify the specific positions that affected ZFP57, CTCF, BATF1, GLI1 and HOXB13 binding, including hemimethylation of one of the two strands .
The zinc-finger Cys2His2-like (C2H2) TF KAISO (also called ZBTB33) was first described to bind mCpGs by EMSA in vitro  as well as by a structural report that showed the molecular basis for KAISO bi-modal recognition of both unmethylated and methylated CpGs . Re-analysis of KAISO binding sites and DNA methylation patterns in vivo suggests that KAISO does not bind methylated DNA but rather highly active promoters marked with high levels of acetylated histones , although this interpretation does not take DNA methylation dynamics into account.
Other zinc-finger proteins were identified as readers of methylated DNA , such as KAISO-like ZBTB4 and ZBTB38 in transient transfections in mice . ZFP57 is a well-known example of a TF that binds methylated DNA at imprinted regions in the mouse genome in vivo [59–61] and its mCpG binding preference was shown to be asymmetric .
Another zinc-finger family protein KLF4 was identified as binding mCpGs in a proteomics-based approach and by DNA pull-down . Re-analysis of KLF4 binding sites in mouse ES cells in vivo identified 18.5% as methylated . It was also found by a methylated microarray approach although it could bind both methylated and unmethylated sites, but displaying different sequence preferences . Re-analysis of KLF4 binding sites in human ES cells in vivo identified that of the KLF4 binding sites having a CpG, 48% were methylated , but represented only 3% of all binding sites . Probing the methylation levels of KLF4 binding sites at four different loci in vivo by ChIP bisulfite sequencing found that it could bind two unmethylated sites (TACpGCC) and two methylated sites (CCmCpGCC) .
Several members of the BZIP CEBP TF family were found to bind mCpGs. CEBPA was shown to bind mCpG within the CRE motif by EMSA in vitro . An approach using methylated binding microarrays found that mCpGs promoted binding of CEBPA and CEBPB although CEBPD and CEBPG, which bind similar motifs, were inhibited . Profiling of CEBPB by ChIP-seq in vivo identified only 11% of its methylated motifs as bound, in contrast with 54% of its unmethylated motifs, located in open-chromatin regions . However, TFs are known not to bind all their motif occurrences in a certain context. Further, a similar re-analysis identified 25% of the CEBPB binding sites as methylated , which is surprisingly high since most TF binding sites are located in unmethylated open chromatin regions. More recently, a methyl-SELEX approach identified CEBPB only as weakly binding to mCpGs (called methyl-plus) as well as CEBPE and CEBPG . A different methyl-SELEX approach found that CEBPB could bind both methylated and unmethylated sequences suggesting that CEBPB could tolerate DNA methylation . The methyl-SELEX approach identified several other TFs that could bind to mCpGs . OCT4 (also called POU5F1) was classified as a methyl-plus TFs although it does not have a CpG in its canonical motif. They further tested its sensitivity in vivo by profiling OCT4 binding by ChIP-seq in WT and DNMTs TKO mouse ES cells and could identify a few sites that lost OCT4 binding in absence of DNA methylation, suggesting that OCT4 requires DNA methylation at these sites in vivo .
Additionally, several HOX TFs were also classified as methyl-plus TFs, with some containing CpGs in their motifs. They showed that HOXC11 could specifically drive luciferase activity of an exogenously inserted construct containing its motif only when methylated . HOXA9 was also previously found to bind mCpGs by EMSA . Structural reports further showed the recognition of HOX TFs to mCpGs, such as HOXB13  and the PBX-HOXA9 complex . However, in the case of HOXB13, only the mCpG on the top strand contributes to binding whereas the other strand does not .
Structural reports proposed mechanisms for TF binding to mCpGs. Studies on the HOX TFs suggest that a mCpG in their motif could mimic a thymidine base, which could explain different sensitivities among HOX paralogs, and could be generalized to other TFs [44,45]. Other structural reports suggest that the binding of several TFs to mCpGs such as KAISO, ZFP57 and KLF4 depends on an arginine preceding the first zinc-binding histidine (called the arginine-histidine (RH) motif) , although the presence of an RH motif in zinc-finger proteins may not be a good predictor of mCpG binding [65,66].
These different studies identified many TFs as able to bind mCpGs and therefore were sensitive to DNA methylation (Figure 2B), suggesting a widespread mechanism. However, most results report in vitro affinities and some are contradictory. Recent studies have compiled the methylation status of TF binding sites from in vivo datasets [67,68] although those analyses are static and only correlative. Therefore, the functionality of TFs binding to mCpGs remains to be further investigated experimentally.
DNA methylation does not affect TF binding
Alongside the discovery of TFs sensitive to DNA methylation, others appeared not to be affected and are therefore called insensitive to DNA methylation (Figure 2C). In 1988, SP1 was the first TF to be described as insensitive to mCpGs located both at the center and at the periphery of the SP1 motif by EMSA in vitro [69,70]. However, a later study found that an mCpG affected SP1 binding in vitro and that the aberrant methylation of the retinoblastoma gene promoter in cancer was suggested to prevent SP1 binding in vivo . The YY1 TF was then identified as insensitive to DNA methylation at the Surf genes promoter whereas ETS TFs were blocked by mCpGs .
More recently high-throughput in vitro approaches identified more TFs as insensitive to DNA methylation. A study in plants identified 79 insensitive TFs (24% out of 327 tested) . An approach using methyl-SELEX and bisulfite-SELEX found that 40% (202 out of 519) of the tested TFs were not affected and the majority (84%, 169 out of 202) did not have CpGs in their motifs .
The first evidence for TF insensitivity in vivo came in 2011, when upon removal of DNA methylation (using DNMTs TKO cells), CTCF binding sites were globally unaltered, suggesting that DNA methylation was not preventing CTCF binding in WT mouse ES cells . This was surprising knowing that CTCF was a well-known example of TF sensitivity to DNA methylation (Box 2) and that a similar approach identified NRF1 as sensitive . CTCF as well as REST were then validated using stable insertions of methylated reporters containing their motifs where they could bind methylated regions and lead to local demethylation . In fact, relatively few new regions of open chromatin bound by TFs were identified in absence of DNA methylation in mouse ES cells suggesting that most TFs did not seem to be affected by DNA methylation in vivo in this cell type .
So far, the main feature predicting TF sensitivity appears to be the presence and position of CpGs within each TF motif, where only TFs that recognize a specific CpG within a motif could be affected when this CpG is methylated. Even though some TFs can recognize alternative motifs containing mCpGs (e.g. OCT4), those binding sites are not prevalent in vivo. When analyzing the non-redundant JASPAR CORE 2018 vertebrate motifs , and assuming that the 579 motifs reflect well the 1600 human TFs , we found that 70% of the TF motifs do not contain prominent CpGs and are therefore more likely not to be affected by DNA methylation. Only 30% of the motifs contain CpGs and might therefore be directly sensitive to DNA methylation including 3% with two CpGs such as NRF1 (Figure 3).
Proportion of TFs having CpGs in their motifs
Impact of TF binding on DNA methylation
TF binding promotes DNA methylation
Although the DNMT enzymes have long been identified as writers of DNA methylation, the players involved in the recruitment of DNMT3A and DNMT3B to establish DNA methylation de novo remains elusive. Some TFs have been shown to recruit DNMTs and promote DNA methylation at specific sites (Figure 4A). For example, this was shown for the leukemia-promoting PML-RAR fusion protein at the RARβ2 promoter , MYC through association with MIZ1 at the p21Cip1 gene , PU.1 in reporter constructs , E2F6 at germline-specific genes , ZFP354B at the FAS promoter in NIH 3T3 cells  and ZNF304 at the INK4-ARF promoter in human ES cells . However, those are few examples by distinct TFs at specific sites, some of which could result from indirect effects. Considering that in most cell types, most of the genome is methylated [19–21], the establishment of DNA methylation de novo could also be unspecific to the DNA sequence and rely on the chromatin context and other regulatory mechanisms or be a default state .
Possible scenarios by which TF binding could impact DNA methylation
TF binding maintains DNA methylation states
Once DNA methylation states are established, mCpG readers, such as the MBD family members, are thought to induce transcriptional repression through chromatin remodeling and maintain DNA methylation states  (Figure 4B). However, MBD binding was shown to correlate with CpG content  and might therefore be limited to regions with high mCpG densities, often located at CpG-island gene promoters . Accumulation of MBD binding to mCpG is also thought to protect from TF binding and DNA demethylation by steric interference but remains to be shown. For example, the inhibition of NRF1 binding genome-wide by DNA methylation in vivo genome-wide is independent of mCpG and MeCP2 density .
Sequence-specific TFs might also be involved in maintaining DNA methylation states. For example, ZFP57, which binds to mCpGs, was found to be necessary for the maintenance of DNA methylation and the histone 3 lysine 9 tri-methylation (H3K9me3) histone mark at these sites . However, more examples of such a mechanism remain to be found and TFs able to bind mCpGs have recently been suggested to induce demethylation.
TF binding triggers demethylation
Recently, a model where TFs would instruct DNA methylation patterns has emerged (Figure 4C). TFs that are insensitive or bind to mCpGs were shown to induce active demethylation through ten-eleven translocations (TETs) recruitment [13,14] (Box 1). This was first shown for CTCF and REST using stable insertions of methylated reporters in mouse ES cells where they could bind and induce local demethylation . It was later demonstrated for PU.1 by motif analyses in differentially methylated regions , for RUNX1, RUNX3, GATA2, CEBPB, MAFB, NR4A2, MYOD1, CEBPA and TBX5 by methylation array , for CEBPA, KLF4 and TFCP2L1 by profiling DNA hydroxymethylation  and for EGR1 by ChIP-seq .
The ability of TFs to mediate DNA demethylation was suggested for TFs insensitive to DNA methylation, such as REST, to enable binding of TFs otherwise repressed by DNA methylation, such as NRF1, and speculated to be a useful feature of pioneer TFs . Pioneer TFs, such as GATA4, KLF4 or OCT4, were also proposed to bind mCpGs and then thought to induce DNA demethylation [44,62]. This would explain why TFs that bind mCpGs appear to have unmethylated binding sites in genomic studies. However, most of the TFs described as pioneers or reprogramming do not contain CpGs in their canonical motifs (e.g. OCT4, KLF4, FOXA1, GATA3, CEBPA; Figure 3) and might therefore be insensitive to DNA methylation and able to induce local demethylation of the surrounding region.
TF binding maintains unmethylated states
Once a region is unmethylated, binding of TFs was shown to protect from DNA methylation (Figure 4D). CxxC zinc-finger proteins, such as CXXC1/CFP1, FBXL19 or TET1, were identified to bind clustered unmethylated CpGs within CpG islands, recruit chromatin remodelers and help maintaining unmethylated regions [84–86].
Sequence-specific TFs have also been shown to protect unmethylated regions from de novo methylation. Early on, SP1 was shown to bind unmethylated CpG islands and mutations of its motif resulted in increased DNA methylation at the Aprt gene promoter [87,88] and later at the Gtf2a1l promoter . CTCF was also shown to maintain unmethylated regions at the Igf2-H19 locus  (Box 2).
The maintenance of unmethylated regions by TFs was mainly shown at CpG islands that mark most housekeeping genes in vertebrate genomes . In fact, CpG islands are thought to have emerged during evolution due to high rates of methylated cytosine deamination into thymine (5mC to T), which affected most of the methylated genome except at constitutively active unmethylated promoters .
In CpG-poor regions, reintroduction of de novo DNA methylation could outcompete NRF1 binding suggesting that binding of TFs sensitive to DNA methylation cannot maintain unmethylated states . However, binding of insensitive TFs should both be able to induce demethylation and maintain unmethylated DNA at active regulatory regions (e.g. CTCF, REST).
DNA methylation patterns are highly dynamic during development, cellular differentiation and disease states [92–95]. Since they strongly anti-correlate with chromatin accessibility and TF binding, regions with no or low DNA methylation levels can be good predictors of TF binding sites [20,96]. However, this mere anti-correlation does not inform us about which of the DNA methylation or TF binding comes first.
The established model proposes that DNA methylation patterns instruct TF binding and act as a transcriptional repressor, although the functions of DNA methylation still remain to be further investigated . Recent studies support this model and several TFs were identified as sensitive to DNA methylation, such as NRF1 in vivo genome wide  (Figure 2A). However, although DNA methylation can indeed inhibit TF binding, the extent of this repression in vivo seems limited in mouse ES cells, and remains to be explored in differentiated cell types.
More recently, a new model has emerged where TF binding instructs DNA methylation patterns. Recent studies support this model where TFs, that are insensitive to DNA methylation (e.g. REST or CTCF ) or recognize mCpGs (e.g. OCT4 ), bind methylated regions, induce demethylation through the recruitment of TET enzymes and maintain regions unmethylated by protecting them from de novo methylation (Figure 4C,D). This ability for TFs to induce demethylation in turn enables binding of sensitive TFs to these sites (e.g. NRF1 ), which could explain why DNA methylation seems to have a minimal effect on the binding of some sensitive TFs (e.g. MYC ). This feature is suggested to be useful for pioneer TFs, to bind closed methylated DNA inducing chromatin remodeling, leading to dynamic patterns of TF binding and gene expression.
Although the sensitivity of TFs to DNA methylation was mainly assessed in vitro, it appears to be dependent on the presence of mCpGs at specific positions within TF motifs (Figure 3). Therefore, TF sensitivity is factor- and condition-specific, where the same TF could be sensitive at sites with mCpGs and insensitive when CpGs are unmethylated (e.g. CTCF  (Box 2), OCT4 ).
In order to generalize these models, further in vivo experiments will be required. However, the dynamics of DNA methylation changes driven by TFs cannot be assessed in static conditions where a TF binding site will always appear as unmethylated. Therefore, perturbation or kinetic experiments, facilitated by genome editing tools , will be essential to better understand the interplay between TF binding and DNA methylation.
DNA methylation and TF binding patterns anti-correlate.
Established model: DNA methylation represses TF binding.
Emerging model: TF binding to methylated regions induces demethylation.
TF sensitivity to DNA methylation depends on the position of methylated CpGs within their binding motifs.
The authors thank M. Weber, D. Shlyueva and R. Grand for feedback on the review. We apologize to all colleagues whose work could not be cited due to space limitations.
The authors declare that there are no competing interests associated with the manuscript.
This work was supported by the Systems Biology Cancer Plan grant from the ITMO Cancer AVIESAN (French National Alliance for Life Sciences and Health) [grant number 18CB008-00]; and the Initiative of Excellence IDEX-Unistra from the French national programme ‘Investment for the future’ [grant number ANR-10-IDEX-0002-02].
E. Héberlé and A.F. Bardet wrote the manuscript and edited the figures.
activating transcription factor 4
basic leucine zipper
CCAAT enhancer binding protein β, η, γ
chromatin immunoprecipitation followed by high throughput sequencing
cyclic adenosine 3′,5′-monophosphate response element-binding protein
E2 transcription factor
electrophoretic mobility shift assay
insulin-like growth factor 1
- JUND (AP-1)
JunD proto-oncogene, AP-1 transcription factor subunit
methyl-binding domain 1/2/4
methylation of CpG
methyl-CpG binding protein 2
nuclear respiratory factor 1
octamer binding protein 4/POU class 5 homeobox 1
spleen focus integration oncogene 1
retinoic acid receptor
RE1 silencing transcription factor
regulatory factor X1/5
specificity protein 1
DNMT triple knockout
zinc finger homeobox
zinc finger protein 57 homolog