A possible role of Drosophila CTCF in mitotic bookmarking and maintaining chromatin domains during the cell cycle
© Shen et al. 2015
Received: 20 January 2015
Accepted: 20 May 2015
Published: 27 May 2015
The CCCTC-binding factor (CTCF) is a highly conserved insulator protein that plays various roles in many cellular processes. CTCF is one of the main architecture proteins in higher eukaryotes, and in combination with other architecture proteins and regulators, also shapes the three-dimensional organization of a genome. Experiments show CTCF partially remains associated with chromatin during mitosis. However, the role of CTCF in the maintenance and propagation of genome architectures throughout the cell cycle remains elusive.
We performed a comprehensive bioinformatics analysis on public datasets of Drosophila CTCF (dCTCF). We characterized dCTCF-binding sites according to their occupancy status during the cell cycle, and identified three classes: interphase-mitosis-common (IM), interphase-only (IO) and mitosis-only (MO) sites. Integrated function analysis showed dCTCF-binding sites of different classes might be involved in different biological processes, and IM sites were more conserved and more intensely bound. dCTCF-binding sites of the same class preferentially localized closer to each other, and were highly enriched at chromatin syntenic and topologically associating domains boundaries.
Our results revealed different functions of dCTCF during the cell cycle and suggested that dCTCF might contribute to the establishment of the three-dimensional architecture of the Drosophila genome by maintaining local chromatin compartments throughout the whole cell cycle.
Referred to CTCF, the CCCTC-binding factor is a zinc finger protein highly conserved from Drosophila to human, and is the only known insulator protein in vertebrates [1, 2]. Initially discovered as a repressor of the chicken c-myc gene , CTCF is reported to be involved in many cellular processes, including transcription activation and repression, chromosome insulation, X-chromosome inactivation, DNA replication, and nucleosome positioning [4–8].
The many functions of CTCF can now be viewed in the context of genome-wide analyses. Researchers have identified hundreds of thousands of CTCF-binding sites across the genomes in different tissues of different species [9–12]. The widespread distribution of CTCF can be attributed to the interactions between the zinc finger domains of this protein and specific DNA sequences; CTCF can bind divergent sequences by using different combinations of its 11 zinc fingers . The canonical CTCF binding motif is 20 bp ; however, using new technique and large-scale data, researchers identified a 33/34-mer two-part CTCF motif in mammals [11, 14]. CTCF-binding sites with larger motifs usually show stronger ChIP signal enrichment and are more conserved [11, 13, 14].
With the development of in vivo imaging techniques and molecular methods based on proximity ligation (3C, 4C, Hi-C, etc.) [15–18], emerging evidence suggests that genomes are dynamically organized at multiple structural levels, and that the hierarchical three-dimensional structure of chromatin is remarkably important for cellular function . It is possible that CTCF, through using different combinations of zinc fingers, interacting with different protein partners, and the last but not the least, employing various post-translational modifications, could mediate extensive intra- and inter-chromatin interactions [7, 8]. Furthermore, research strongly suggests that CTCF clusters with other architecture proteins, and that CTCF-binding sites are enriched at topologically associating domain boundaries in mammalian and Drosophila genomes [20–22]. Thus, it is likely that CTCF plays a conserved role in chromatin domain organization. In addition, CTCF may be the main component of the heritable epigenetic system, regulating the interplay between DNA methylation, nuclear architecture, and lineage-specific gene expression.
Recently, there has been a growing interest in how the transcription program is re-established during mitosis. Several transcription factors, including CTCF, have been documented to remain bound to mitotic chromatin [23–29], and CTCF is also reported to function during the entire cell cycle .
However, changes of CTCF-binding sites during the cell cycle, and functions of this protein at different cell cycle phases, remain largely uncharacterized. Could CTCF act as mitotic bookmarkers that help establish, maintain and propagate the genomic topological organization during the cell cycle is also unknown. In this report, we analysed dCTCF binding site in Drosophila genome using public available datasets, and identified sites that are bound in interphase and mitosis, only during mitosis and only during interphase. Further, we found differences in conservation, binding motives, and GO enrichments among these three classes of dCTCF-binding sites. In addition, we observed that dCTCF-binding sites of the same class preferentially localized closer to each other, and were highly enriched at chromatin syntenic and topologically associating domains boundaries. Thus, dCTCF might contribute to the three-dimensional architecture of the Drosophila genome by maintaining local chromatin compartments throughout the whole cell cycle.
Cell cycle phase-specific dCTCF binding sites
Our de novo motif analyses revealed that dCTCF-binding sites were strongly enriched for the conventional dCTCF motif (Fig. 1b, module A) and also another slightly enriched motif (Fig. 1b, module B), which was previously reported in human and mouse . We found this two-part motif was present in 60 % dCTCF IM sites and 51 % IO sites, but only 29 % MO sites. It was interesting that MO sites demonstrated low enrichment for this motif, but rather contained “GCW” repeats. By examining the GC content of all dCTCF-binding sites (Additional file 2: Figure S1B), we found dCTCF preferentially bound to GC-rich regions during mitosis. These data indicate dCTCF bound to DNA through different mechanism during interphase and mitosis phase. In addition, we identified a consensus Su(Hw) motif within CP190 IO sites, and a dCTCF motif within CP190 IM sites (Additional file 2: Figure S1C).
dCTCF IM sites were tightly bound and highly conserved
During both interphase and mitosis, the average ChIP signal of dCTCF IM sites was significantly higher than that of IO or MO sites (Fig. 1d), which showed that IM sites were more tightly bound by dCTCF during the whole cell cycle. Meanwhile, motif analysis discovered that most dCTCF IM sites contained a two-module motif, indicating the relative high binding strength.
We performed the conservation analysis of dCTCF among different Drosophila species, and observed that dCTCF IM sites were significantly more conserved than IO sites (p < 0.0001, Fig. 1c). As no available data existed regarding MO sites in Drosophila species aside from D. melanogaster, it was hard to confirm whether MO sites were conserved between among different fly species.
Different dCTCF-binding sites play different roles during the cell cycle
To gain insight into the functional differences of IM, IO, and MO sites, we performed GO analysis for genes associated with dCTCF-binding sites using DAVID  (Fig. 2b, Additional file 3: Table S2). Genes associated with IM sites showed a broad range of ontologies, indicating that IM sites were involved in many biological processes during the whole cell cycle. On the contrary, IO and MO showed distinct enriched ontologies. Particularly, IO sites were enriched at genes involved in neuron differentiation and development, which is possibly the result of the original ChIP data being collected from embryonic D. melanogaster Kc167 cells, a cell type in which dCTCF may be involved in the regulation of this major cellular process. MO sites were enriched at genes involved in microtubule cytoskeleton organization and cell cycle, indicating that these dCTCF-binding sites might play a key role in mitotic progression.
We subsequently evaluated co-occupancy status of dCTCF with other architecture proteins and epigenetic markers (Additional file 4: Figure S2). It was found dCTCF, BEAF-32, and CP190 co-localized mainly in gene TSS regions, and usually had high H3K4me3 and H3K27ac signals, which were also high around dCTCF MO sites distributed in coding exons during interphase. Although the intensity of these signals during mitosis was unknown, these results implied that dCTCF MO sites could be involved in specific regulatory functions during mitosis. Interestingly, BEAF-32 and CP190 were barely detectable near dCTCF MO sites. As dCTCF demonstrated more MO sites compared to BEAF-32 and CP190, it is likely that dCTCF functions differently from other insulators during mitosis.
dCTCF-binding sites of the same class tend to maintain local compartments
When we compared with syntenic blocks, we further found that dCTCF-binding sites were also in fact enriched at syntenic boundaries (Fig. 4b, c). These results showed that dCTCF-binding sites were enriched at not only higher-order spatial physical and functional domain boundaries, but also evolutionary conserved linear domain boundaries. With similar method to determine specific TAD boundaries afore mentioned, we finally checked distribution of IM/IO/MO syntenic block boundaries and random ones. Once again, we found that syntenic block boundaries belonging to the same group tend to be closer than random sets (Fig. 4e, f). With these findings, we hypothesized that dCTCF could contribute to the formation of boundaries at various chromatin domain levels, and that dCTCF-binding sites of the same class may maintain local chromatin compartments.
Here, we investigated the genome-wide distributions of three types of dCTCF-binding sites (i.e. IM, IO, and MO) in D. melanogaster. A large proportion of binding sites remained bound by dCTCF during mitosis. dCTCF IM sites usually contained two-module motifs, which was consistent with their stronger ChIP enrichment and higher conservation, implying that dCTCF IM sites might play important roles during the cell cycle and need to be retained during mitosis. Recent work in human cells revealed that, contrary to classical insulator function, most CTCF-binding sites mediate promoter-enhancer communication . Coincidently, we found a larger proportion of IM than IO sites that located within TSS. It should be interesting to further examine whether TSS, and hence promoter-associated dCTCF IM sites actually promote rapid activation of gene expression during M/G1 transition.
The proportion dCTCF MO binding sites was much larger than that of recently reported mitotic bookmarking factors, for examples, MO sites for Myc, Polycomb Group (PcG) proteins, GATA1, FoxA1, and RBPJ accounted for only 0 % , 0 % , 9.9 % , 13.7 % , and 10.6 % , respectively; and this proportion is also higher than other architecture proteins in D. melanogaster, such as BEAF-32 and CP190, which were only 1 % and 4 %, respectively (Additional file 1: Table S1, Additional file 2: Figure S1A). Besides, these novel dCTCF MO sites demonstrated particular binding motives, agreeing with previous reports on GATA1 and FoxaA1, which also occupy mitotic-specific binding sites with different sequence features (compared to interphase binding sites) [27, 28]. As chromatin is condensed during mitosis, it is possible that the recognition pattern of certain transcription factors might change during this particular phase of the cell cycle. Such a substantial proportion of dCTCF MO sites and special binding motif promoted us to suspect particular roles these site might play in mitosis. Consequently, GO analysis implied that dCTCF MO sites were enriched at genes involved in microtubule cytoskeleton organization and cell cycle, indicating that these sites mediated dCTCF functions in mitotic progression, though the underlining molecular mechanism remains unclear. We observed that dCTCF preferentially bound to exon regions during mitosis, possibly due to reduction of open chromatin regions during mitosis. CTCF has been shown to bind to exons and plays a role in alternative splicing , whether dCTCF functions in a similar manner during mitosis in Drosophila remains to be an interesting open question.
Previous studies have suggested that dCTCF and Su(Hw) do not interact with one another directly, but rather co-localize with CP190 individually at distinct loci . Our motif analysis identified a consensus Su(Hw) motif within CP190 IO sites, and a dCTCF motif within CP190 IM sites (Additional file 2: Figure S1C), as Su(Hw) was not typically present on chromosomes during mitosis , we postulated that CP190 recruited by Su(Hw) to chromosomes may be erased during mitosis while CP190 recruited by dCTCF may be maintained. Thus providing an extended view of these architecture proteins.
Several investigations have observed that a portion of the FoxA1, RBPJ, and H2A.Z binding peaks obtained by ChIP-seq appeared to shift in their genomic localization when mitotic and asynchronous interphase cells were compared. It was previously proposed that such peak shifting facilitated the rebinding of specific proteins to real target sites during the exit from M phase, which was accomplished by the partial maintenance or storage of mitotic bookmarking DNA binding proteins on mitotic chromatin. However, we did not observe a clear peak shifting phenomenon for dCTCF, which implied that the binding of dCTCF to MO sites might have functions other than the storage of this protein on mitotic chromatin. We suspected that dCTCF was maintained on mitotic chromatin by a different yet unknown mechanism. Further experiments, such as comparison of ChIP-seq results from different cell types of different cell cycle phases, will provide more clues about the binding mechanism and functions of CTCF during the cell cycle in different organisms.
Previous studies have demonstrated that partial dCTCF-binding sites serve as boundary elements. Our analysis confirmed this result, and found that dCTCF-binding sites were not only enriched at topological and functional domain boundaries, but also at evolutionarily conserved domain boundaries (i.e., syntenic blocks), which extended the definition of the “genomic domain”. We also found that cell cycle phase specific dCTCF-binding sites of the same class were preferentially bound closer to one another, and a similar result was observed for domain boundaries containing dCTCF-binding sites. Although MO boundaries of TAD did not show the same phenomenon, partially because of the lack of TAD data pertaining to mitosis, and recent Hi-C experiments suggested TADs were mostly lost during mitosis  in human cells. However, these results did not conclusively exclude the possibility of the existence of potential TAD boundaries. In vivo live imaging of potential domain boundaries during the cell cycle might help to rectify this issue. We proposed a model of dCTCF loci during different phases of the cell cycle (Fig. 3d), which indicated that cell cycle phase specific dCTCF-binding sites help to establish cell cycle phase specific chromatin domains. Further experiments such as cell cycle specific dCTCF knock down, or the degradation of this protein, might provide direct biochemical evidence for this model.
We proposed that dCTCF, by collaboration with other architecture proteins and sequence-specific transcription factors, may partially maintain local chromatin domain structures; such structures could be propagated through cell division, and facilitate the rapid re-establishment of chromatin structures and the re-activation of some genes during the M/G1 phase transition.
We determine that dCTCF remains bound to chromatin throughout the cell cycle, and demonstrate several cell cycle phase specific sites with different sequence features. Each class of dCTCF-binding sites (i.e., IM, IO, and MO) plays different roles at different phases of the cell cycle. Specifically, IM sites appear to be involved in many biological processes during the entire cell cycle, whereas MO sites may play a key role only during mitosis. IM sites are also unique among dCTCF-binding sites, as they are highly conserved between Drosophila species, and are intensely bound by dCTCF. In addition to providing insight into the roles of dCTCF during the cell cycle, our results also indicate that dCTCF likely contributes to the formation of various chromatin domain boundaries, as cell cycle phase specific dCTCF-binding sites are closer to one another and appear to maintain cell cycle phase specific local compartments.
The Drosophila melanogaster dCTCF, BEAF-32 and CP190 ChIP-seq data sets were downloaded from Yang et al. , (deposited in NCBI’s Gene Expression Omnibus (GEO) under accession numbers GSE32584). For cell synchronization, Yang et al. treated the cells with hydroxyurea and nocodazole. The enrichment of the mitotic and interphase cell populations showed 97–99 % purity. Histone modification data were from Kellner et al. , (GSE36374). dCTCF ChIP-seq data from other Drosophila species, used for the conservation analysis, were from Ni et al. , (GSE24449). Chromatin topologically associating domain boundaries were identified in Sexton et al. , and regions of conserved synteny across 12 Drosophila species were identified in Bhutkar et al. . Syntenic blocks conserved between human and D. melanogaster were downloaded from Sinha et al. . We only used genomic features data (e.g., genes and binding sites) available for chromosomes 2 L, 2R, 3 L, 3R, 4, and X.
For each insulator, we did de novo motif searching by MEME  using the top 1000 ChIP-seq-enriched binding site DNA sequences. And we also randomly chose 500 sequences from IM/IO/MO sites separately as candidate sets. To calculate the occurrence (%) of the two-part motif in each class of dCTCF-binding sites, we used FIMO tools from the MEME software suite. We downloaded motif matrix from JASPAR  for comparison.
We mapped all the non-D. melanogaster species binding sites into the D. melanogaster genome with liftOver  using default parameters, except for a match of 0.5. The number of D. melanogaster-binding sites overlapping with each of the non-D. melanogaster liftOver binding sites was counted.
Calculation of distances between specific dCTCF-binding sites or domain boundaries
We calculated the shortest distance between two specific sites for each dCTCF-binding sites, that is, we first chose one dCTCF-binding site as the bait, then calculated the distance between the bait and the nearest other site. We used the control sets by randomly shuffling the location of all sites along chromosomes. We termed any domain boundary containing dCTCF IM/IO/MO sites as “IM/IO/MO boundary”, and calculated the minimum domain numbers between specific boundaries. We permuted the dCTCF-binding sites at boundaries as control sets. The distribution of shortest distances, or minimum domain numbers of control sets, were computed 1,000 times followed by normalization.
This research was supported by National Basic Research Program of China [2013CB966802], National Natural Science Foundation of China [31370762, 31030026, 31071119, 81372218] and National Key Technology R&D Program of China [2012BAI01B07].
- Ohlsson R, Renkawitz R, Lobanenkov V. CTCF is a uniquely versatile transcription regulator linked to epigenetics and disease. Trends Genet. 2001;17:520–7.PubMedView ArticleGoogle Scholar
- Moon H, Filippova G, Loukinov D, Pugacheva E, Chen Q, Smith ST, et al. CTCF is conserved from Drosophila to humans and confers enhancer blocking of the Fab-8 insulator. EMBO Rep. 2005;6:165–70.PubMed CentralPubMedView ArticleGoogle Scholar
- Lobanenkov V, Nicolas R, Adler V, Paterson H, Klenova E, Polotskaja A, et al. A novel sequence-specific DNA binding protein which interacts with three regularly spaced direct repeats of the CCCTC-motif in the 5'-flanking sequence of the chicken c-myc gene. Oncogene. 1990;5:1743–53.PubMedGoogle Scholar
- Bell AC, West AG, Felsenfeld G. The protein CTCF is required for the enhancer blocking activity of vertebrate insulators. Cell. 1999;98:387–96.PubMedView ArticleGoogle Scholar
- Xu N, Donohoe ME, Silva SS, Lee JT. Evidence that homologous X-chromosome pairing requires transcription and Ctcf protein. Nat Genet. 2007;39:1390–6.PubMedView ArticleGoogle Scholar
- Fu Y, Sinha M, Peterson CL, Weng Z. The insulator binding protein CTCF positions 20 nucleosomes around its binding sites across the human genome. PLoS Genet. 2008;4:e1000138.PubMed CentralPubMedView ArticleGoogle Scholar
- Phillips JE, Corces VG. CTCF: master weaver of the genome. Cell. 2009;137:1194–211.PubMed CentralPubMedView ArticleGoogle Scholar
- Ong CT, Corces VG. CTCF: an architectural protein bridging genome topology and function. Nat Rev Genet. 2014;15:234–46.PubMedView ArticleGoogle Scholar
- Jothi R, Cuddapah S, Barski A, Cui K, Zhao K. Genome-wide identification of in vivo protein-DNA binding sites from ChIP-Seq data. Nucleic Acids Res. 2008;36:5221–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin D, Pantoja C, Fernández Miñán A, Valdes-Quezada C, Moltó E, Matesanz F, et al. Genome-wide CTCF distribution in vertebrates defines equivalent sites that aid the identification of disease-associated genes. Nat Struct Mol Biol. 2011;18:708–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmidt D, Schwalie PC, Wilson MD, Ballester B, Gonçalves Â, Kutter C, et al. Waves of retrotransposon expansion remodel genome organization and CTCF binding in multiple mammalian lineages. Cell. 2012;148:335–48.PubMed CentralPubMedView ArticleGoogle Scholar
- Ni X, Zhang YE, Nègre N, Chen S, Long M, White KP. Adaptive evolution and the birth of CTCF binding sites in the Drosophila genome. PLoS Biol. 2012;10:e1001420.PubMed CentralPubMedView ArticleGoogle Scholar
- Nakahashi H, Kwon KK, Resch W, Vian L, Dose M, Stavreva D, et al. A genome-wide map of CTCF multivalency redefines the CTCF code. Cell Rep. 2013;3:1678–89.PubMed CentralPubMedView ArticleGoogle Scholar
- Rhee HS, Pugh BF. Comprehensive genome-wide protein-DNA interactions detected at single-nucleotide resolution. Cell. 2011;147:1408–19.PubMed CentralPubMedView ArticleGoogle Scholar
- Hübner MR, Spector DL. Chromatin dynamics. Annu Rev Biophys. 2010;39:471–89.PubMed CentralPubMedView ArticleGoogle Scholar
- Dekker J, Rippe K, Dekker M, Kleckner N. Capturing chromosome conformation. Science. 2002;295:1306–11.PubMedView ArticleGoogle Scholar
- Zhao Z, Tavoosidana G, Sjölinder M, Göndör A, Mariano P, Wang S, et al. Circular chromosome conformation capture (4C) uncovers extensive networks of epigenetically regulated intra- and interchromosomal interactions. Nat Genet. 2006;38:1341–7.PubMedView ArticleGoogle Scholar
- Lieberman-Aiden E, Van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326:289–93.PubMed CentralPubMedView ArticleGoogle Scholar
- Göndör A, Ohlsson R. Chromosome crosstalk in three dimensions. Nature. 2009;461:212–7.PubMedView ArticleGoogle Scholar
- Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Sexton T, Yaffe E, Kenigsberg E, Bantignies F, Leblanc B, Hoichman M, et al. Three-dimensional folding and functional organization principles of the Drosophila genome. Cell. 2012;148:458–72.PubMedView ArticleGoogle Scholar
- Van Bortle K, Nichols MH, Li L, Ong C-T, Takenaka N, Qin Z, et al. Insulator function and topological domain border strength scale with architectural protein occupancy. Genome Biol. 2014;15:R82.PubMed CentralPubMedView ArticleGoogle Scholar
- Burke LJ, Zhang R, Bartkuhn M, Tiwari VK, Tavoosidana G, Kurukuti S, et al. CTCF binding and higher order chromatin structure of the H19 locus are maintained in mitotic chromatin. EMBO J. 2005;24:3291–300.PubMed CentralPubMedView ArticleGoogle Scholar
- Nekrasov M, Amrichova J, Parker BJ, Soboleva TA, Jack C, Williams R, et al. Histone H2A.Z inheritance during the cell cycle and its impact on promoter organization and dynamics. Nat Struct Mol Biol. 2012;19:1076–83.PubMedView ArticleGoogle Scholar
- Yang J, Sung E, Donlin-Asp PG, Corces VG. A subset of Drosophila Myc sites remain associated with mitotic chromosomes colocalized with insulator proteins. Nat Commun. 2013;4:1464.PubMed CentralPubMedView ArticleGoogle Scholar
- Follmer NE, Wani AH, Francis NJ. A polycomb group protein is retained at specific sites on chromatin in mitosis. PLoS Genet. 2012;8:e1003135.PubMed CentralPubMedView ArticleGoogle 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 CentralPubMedView ArticleGoogle Scholar
- Caravaca JM, Donahue G, Becker JS, He X, Vinson C, Zaret KS. Bookmarking by specific and nonspecific binding of FoxA1 pioneer factor to mitotic chromosomes. Genes Dev. 2013;27:251–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Lake RJ, Tsai PF, Choi I, Won KJ, Fan HY. RBPJ, the major transcriptional effector of Notch signaling, remains associated with chromatin throughout mitosis, suggesting a role in mitotic bookmarking. PLoS Genet. 2014;10:e1004204.PubMed CentralPubMedView ArticleGoogle Scholar
- Heath H, de Almeida CR, Sleutels F, Dingjan G, van de Nobelen S, Jonkers I, et al. CTCF regulates cell cycle progression of alphabeta T cells in the thymus. EMBO J. 2008;27:2839–50.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- Sanyal A, Lajoie BR, Jain G, Dekker J. The long-range interaction landscape of gene promoters. Nature. 2012;489:109–13.PubMed CentralPubMedView ArticleGoogle Scholar
- Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, et al. CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011;479:74–9.PubMedView ArticleGoogle Scholar
- Bushey AM, Ramos E, Corces VG. Three subclasses of a Drosophila insulator show distinct and cell type-specific genomic distributions. Genes Dev. 2009;23:1338–50.PubMed CentralPubMedView ArticleGoogle Scholar
- Naumova N, Imakaev M, Fudenberg G, Zhan Y, Lajoie BR, Mirny LA, et al. Organization of the mitotic chromosome. Science. 2013;342:948–53.PubMed CentralPubMedView ArticleGoogle Scholar
- Kellner WA, Ramos E, Van Bortle K, Takenaka N, Corces VG. Genome-wide phosphoacetylation of histone H3 at Drosophila enhancers and promoters. Genome Res. 2012;22:1081–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhutkar A, Schaeffer SW, Russo SM, Xu M, Smith TF, Gelbart WM. Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes. Genetics. 2008;179:1657–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Sinha AU, Meller J. Cinteny: flexible analysis and visualization of synteny and genome rearrangements in multiple organisms. BMC Bioinformatics. 2007;8:82.PubMed CentralPubMedView ArticleGoogle Scholar
- Bailey TL, Elkan C. Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994;2:28–36.PubMedGoogle Scholar
- Mathelier A, Zhao X, Zhang AW, Parcy F, Worsley-Hunt R, Arenillas DJ, et al. JASPAR 2014: an extensively expanded and updated open-access database of transcription factor binding profiles. Nucleic Acids Res. 2014;42:D142–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Karolchik D, Barber GP, Casper J, Clawson H, Cline MS, Diekhans M, et al. The UCSC Genome Browser database: 2014 update. Nucleic Acids Res. 2014;42:D764–70.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.