DamID profiling of dynamic Polycomb-binding sites in Drosophila imaginal disc development and tumorigenesis
© The Author(s) 2018
Received: 20 January 2018
Accepted: 21 May 2018
Published: 5 June 2018
Tracking dynamic protein–chromatin interactions in vivo is key to unravel transcriptional and epigenetic transitions in development and disease. However, limited availability and heterogeneous tissue composition of in vivo source material impose challenges on many experimental approaches.
Here we adapt cell-type-specific DamID-seq profiling for use in Drosophila imaginal discs and make FLP/FRT-based induction accessible to GAL driver-mediated targeting of specific cell lineages. In a proof-of-principle approach, we utilize ubiquitous DamID expression to describe dynamic transitions of Polycomb-binding sites during wing imaginal disc development and in a scrib tumorigenesis model. We identify Atf3 and Ets21C as novel Polycomb target genes involved in scrib tumorigenesis and suggest that target gene regulation by Atf3 and AP-1 transcription factors, as well as modulation of insulator function, plays crucial roles in dynamic Polycomb-binding at target sites. We establish these findings by DamID-seq analysis of wing imaginal disc samples derived from 10 larvae.
Our study opens avenues for robust profiling of small cell population in imaginal discs in vivo and provides insights into epigenetic changes underlying transcriptional responses to tumorigenic transformation.
Understanding the in vivo dynamics of DNA binding by chromatin regulatory proteins is key to elucidate the molecular basis of cell behaviours ranging from differentiation to adaptation and plasticity. The model system Drosophila has contributed tremendously to our understanding of chromatin dynamics during developmental transitions, stem cell differentiation and also tumorigenesis. Yet, like other in vivo model systems, the small size and the heterogeneous fate composition of Drosophila tissues still pose challenges to the detailed tracking of DNA binding sites in different cell populations and lineages in vivo.
Several experimental approaches to overcome these challenges have been developed. For example, chromatin immunoprecipitation (ChIP) protocols use fluorescence-activated cell sorting (FACS) or immunoprecipitation (IP)-based methods to isolate Drosophila cell populations from tissues [1–4]. These approaches, however, still require a significant amount of input material for reproducible results, which has prevented these methods from being used in contexts where small source tissues, such as imaginal discs, are routinely isolated by manual dissection. Alternatively, recent publications establish cell-type-specific DamID profiling in Drosophila brains [5–8]. PCR-amplified tracking of adenine methylation (m6A) conferred by DamID to GATC sequence motifs and the absence of IP steps significantly reduces the input material required for DamID . Moreover, m6A is only generated in cell types expressing DamID constructs; therefore, DamID protocols do not necessitate to physically isolate cell populations from complex tissues [5, 6]. Thus, DamID is a very attractive technology to profile small and even rare cell populations in vivo.
We wanted to adapt the inducible FRT/FLP-out DamID system described for Drosophila brains  to cell-type-specific profiling in imaginal discs. These small tissues have a rich history as model to study developmental patterning, tumorigenesis and regeneration  but are mostly accessed by manual dissection for experimental analysis. We wanted to establish versatility of targeting DamID expression to specific cell types by enabling the use of GAL4 driver lines available in these tissues. While the TaDa-DamID system [6, 8] also utilizes cell-type-specific targeting by GAL4 drivers, TaDa depends on acute expression patterns of a chosen GAL4 driver at the time of analysis. In contrast, we aimed to target DamID to specific cell lineages enabling tracking of DNA binding sites in parental and descendant populations—independent of whether the GAL4 driver used was still active in descendant cells. Furthermore, while the FRT/FLP-out DamID has been suggested to be compatible with GAL4-dependent targeting , its cell-type specificity and experimental feasibility have not yet been tested. Finally, we sought to establish a proof of principle that a limiting amount of manually dissected imaginal disc material is sufficient to sensitively detect changes in DNA binding activity in development and disease.
More specifically, we asked whether DamID may be suitable to track the epigenetic regulator Polycomb (Pc) in wing imaginal discs (WIDs) during different developmental stages and tumorigenic transformation. Polycomb is the founding member of the Polycomb group (PcG) family of proteins who form different complexes, such as the Polycomb Repressive Complexes 1 and 2 (PRC1 and PRC2). PcG proteins epigenetically silence genes required for fate specification, cell cycle progression and tissue growth by modulating multiple histone modifications [11–15]. Previous studies demonstrated that PcG protein binding sites change dynamically throughout early embryonic development and suggested that a number of Pc target genes, like JAK/STAT cytokines of the unpaired (upd) family, may be silenced by Pc to suppress tumorigenesis [16–20]. In fact, a significant overlap between PcG target genes and genes upregulated in neoplastic WIDs mutant for the epithelial polarity regulator scribbled (scrib) has been described . However, direct experimental evidence for dynamic Pc-binding at co-regulated candidate genes is still outstanding.
We report here the co-regulation of multiple oncogenic genes by dynamic Pc-binding, while also identifying at least two novel Pc target genes involved in scrib tumorigenesis. We furthermore suggest that gene regulation by Atf3 and AP1 transcription factors as well as modulation of insulator function plays crucial roles in dynamic PcG behaviour. We establish these findings by DamID-seq analysis of wing imaginal discs samples derived from as little as 10 larvae. We furthermore describe a versatile GAL4-driven cell lineage-specific DamID system that can be used for DamID-seq profiling in many Drosophila tissue.
Establishment of versatile GAL4-dependent control of cell lineage-specific DamID
We wanted to optimize this inducible DamID system for flexible cell-type-specific targeting by the rich repertoire of GAL4 driver lines available. We thus screened a number of UAS-FLP constructs from different sources for their ability to mediate efficient removal of the FRT-flanked transcriptional terminator cassette. Moreover, we specifically searched for a UAS-FLP line that did not show leaky expression in the absence of a GAL4 driver to prevent unspecific removal of the terminator cassette. Indeed, combining a UAS-FLP(JD2) transgene  with the inducible DamID system caused GAL4-independed removal of the terminator cassette (Additional file 1: Fig. S1B′). In contrast, a UAS-FLP(EXEL) transgene  did not induce removal of the terminator cassette in WIDs in the absence of a GAL4 driver (Additional file 1: Fig. S1B″). Only combining a DamID;UAS-FLP(EXEL) line with a rotund(rn)GAL4 driver caused partial removal of the terminator cassette in WIDs, consistent with the restricted expression of rnGAL4 in the central domain of the disc (Additional file 1: Fig. S1B″). This region was visualized using the G-trace system (Additional file 1: Fig. S1C) , which maps cell lineage history and real-time expression of a GAL4 driver of choice. To prove that Dam and Dam-Pc fusion proteins are really expressed in a cell-type-specific and GAL4/UAS-FLP(EXEL)-dependent manner, we sought to visualize expression of the Myc-tag encoded by both constructs [5, 7]. To this end, we induced removal of the terminator cassette by crossing a stable DamID;UAS-FLP(EXEL) line to a patched(ptc) GAL4 driver. ptcGAL4 is active in a row of cells anterior to the anterior–posterior compartment boundary in WIDs (Additional file 1: Fig. S1C′). However, most of the anterior compartment derives from cells that had expressed ptc earlier during development (Additional file 1: Fig. S1C′). Thus, the early removal of the terminator cassette during development under the control of ptcGAL4 is expected to cause expression of Myc-tagged Dam and Dam-Pc proteins in all cells of the anterior WID compartment. Notably, Dam and Dam-Pc proteins expressed under the control of the heat-shock promoter are present at undetectable levels if flies were kept at 21 °C. However, if boosted by a heat shock (see Experimental procedures), high expression of the Myc-tag could be detected specifically in the anterior compartment, if FLP expression was induced by ptcGAL4 (Fig. 1c, d). Importantly, Myc-tag expression was completely absent in the posterior compartment. Similarly, when DamID was induced using the posterior compartment driver engrailed(en)GAL4, boosted expression of the Myc-tag was exclusively detected in the posterior compartment (data not shown). These results indicate that UAS-FLP(EXEL) allows for the specific and flexible induction of cell-type-specific DamID in WIDs under the versatile control of cell-type-specific GAL4 drivers.
High expression levels of Dam are known to interfere with DamID specificity  and viability  (Additional file 1: Fig.S1D, E). Therefore, to understand whether expression of Dam by the low basal activity of the Hsp70 promoter at 21 °C is suitable for DamID profiling by maintaining wing disc cell viability, we monitored the occurrence of mitosis and apoptosis by immunodetection of phospho-H3S10 (pH3) and the activated effector caspase Dcp-1, respectively. No differences in mitotic or apoptotic activity between the anterior and posterior compartment could be observed when larvae were maintained at 21 °C and the terminator cassette was removed under the control of ptcGAL4/UAS-FLP(EXEL) (Fig. 1e–f″). Furthermore, immunodetection of developmental regulators such as Ptc itself (Fig. 1g–g″) or wingless (Wg) (Fig. 1h–h″) revealed appropriate patterning activity, and adult wings arising from these discs displayed only subtle alterations, such as extra vein tissues (Fig. 1i–i″). Combined these results suggest that inducible DamID profiling does not interfere with WID viability and developmental progression and thus presents an excellent option for cell-type-specific mapping of DNA binding sites in WIDs in vivo.
DamID and ChIP profiles of Polycomb-binding sites correlate
We first induced ubiquitous expression of Dam and Dam-Pc in whole larvae using a FLP under the control of a heat-shock promoter (hsflp). We isolated and amplified methylated genomic DNA from WIDs of 10 WT or scrib 1 third-instar larvae expressing either Dam alone or a Dam-Pc fusion protein (Fig. 2b) and generated NGS libraries using protocols devoid of additional PCR amplification steps to avoid PCR biases (see Experimental procedures).
The PCR-free NGS library preparation from 20 WIDs generated sequencing profiles with relatively low correlation coefficients across replicates (Additional file 1: Fig. S1F), likely due to high noise in profiles. However, assessment of multiple reproducibility parameters, such as correlation coefficients (Additional file 1: Fig. S1F), hierarchical clustering approaches using 94 DamID-Seq profiles (Additional file 1: Fig. S3A) and autocorrelation of neighbouring GATC sites at Lag 2 (Additional file 1: Fig. S3B) , revealed that technical replicates within genotypes are always more similar to each other than replicates across genotypes. Thus, PCR-free DamID-seq libraries can reproducibly reveal DNA binding profiles for small in vivo tissue samples.
While a subset of PcG target genes was previously reported to be upregulated in scrib 1 WIDs , we found that total levels of H3K27 modifications were comparable between WT and scrib 1 WIDs (Additional file 1: Fig. S4A). Our DamID-seq profiles confirmed that Pc-binding at individual sites (as defined by any genomic sequences flanked by Dam-targeted GATC motifs, also referred to as GATC fragments hereafter) was not globally altered in scrib 1 (Fig. 2c). Indeed, when the genome-wide distribution of Pc-binding intensities at these sites was compared, the correlation between WT and scrib 1 discs (Pearson’s correlation, r = 0.47, Fig. 2d) was only slightly lower than for biological replicates (Pearson’s correlation r = 0.51, Additional file 1: Fig.S1F). Importantly, broad binding of Pc to the Bithorax complex (BX-C) observed in Pc-DamID profiles could also be detected in Pc ChIP profiles from S2 cells, DmBG3 cells and whole embryo (Fig. 2e) [27, 28]. The Pearson’s correlation coefficients calculated for a comparison of the genome-wide Pc-binding intensities at individual GATC fragments in our Pc-DamID-seq and the corresponding GATC fragments in individual Pc ChIP-chip profiles ranged from 0.25 to 0.4 (Fig. 2e′). This finding is in agreement with previous comparisons of the two techniques [29–31] (for example Pearson’s correlation r = 0.37 in ). Our analysis thus indicates that DamID-seq is a suitable method to reveal DNA binding profiles of Polycomb in WID in vivo.
Polycomb-binding is altered only at a subset of target sites in scrib 1 wing discs
To understand whether alterations in Pc-binding at specific target genes may contribute to tumour phenotypes in scrib 1 disc, we performed a three-state hidden Markov model (HMM) analysis of Pc-binding at individual GATC fragments to define ‘depleted’, ‘intermediate’ and ‘enriched’ Pc-binding states and analysed transitions between these states when comparing scrib 1 to WT WIDs (see Experimental procedures, Additional file 1: Fig. S4B, Additional files 2: SF2 and 3: SF3) [26, 32–35]. As expected, we obtained three possible clusters that described the changes between the two profiles, namely (1) ‘no change’, which defined GATC fragments that did not vary in their Pc-binding classification between WT and scrib 1 profiles, irrespective of whether these sites were bound by Pc in WT and scrib 1 WIDs or not; (2) ‘loss’ defined GATC fragments, which were bound by Pc in WT but not in scrib 1 discs; and (3) ‘gain’ defined GATC fragments, which were not bound by Pc in WT but in scrib 1 WID samples. This analysis revealed that about 11% of ‘intermediate’ and ‘enriched’ Pc-binding states present in WT were lost in scrib 1 WIDs and about 18% of scrib 1 ‘intermediate’ and ‘enriched’ Pc-binding states were arising de novo (Fig. 2f). This suggests that Pc-binding dynamics are altered in a loci-specific manner in scrib 1 discs.
Polycomb-binding at differentiation and tumour-associated targets is altered in scrib 1 discs
Our approach indicated the presence of multiple genes associated with transcriptional upregulation upon loss of Pc-binding (group III) and with transcriptional repression upon gain of Pc-binding (group II) in scrib 1 (Fig. 3b), which is consistent with the described function of Pc in gene silencing [11–15]. We thus focused our subsequent analysis on these genes.
Surprisingly, group II included genes implicated in axon guidance, for example dsx, Lrt, caps, PlexB, pdm3, Toll-7 and Fas3 (Fig. 3b), possibly reflecting a failure to develop wing and thorax sensory neurons. While all group II genes gained Pc-binding for at least one taGATCf in scrib 1 discs, we wanted to provide additional evidence for a role of PcG in regulating their expression. An analysis of transcript levels in WIDs mutant for the PRC1 components Psc/Su(z)2  revealed that specifically dsx, Toll-7 and the neuronal Notch target pnt were upregulated upon loss of repressive PcG complex function (Fig. 3c, Additional file 1: Fig. S4E). This suggests that at least a subset of group II genes are bona fide Pc target genes.
Despite these different behaviours with respect to H3K27me3 modifications, Ets21C, Atf3, Ilp8, SOCS36E, upd3 and chinmo are all upregulated upon loss of repressive PRC1 complex function in Psc/Su(z)2 mutant WIDs, demonstrating a role for Pc in silencing these tissue-stress-responsive genes in wild-type WIDs (Additional file 1: Fig. S3D). Thus, we identify at least three tumour-associated genes as novel bona fide Pc target genes and imply that the tumour-suppressive function of PcG proteins  integrates with regulation by the two important tumour-promoting pathways JNK and JAK/STAT.
Modulation of Polycomb-binding and target gene expression is associated with enrichment of specific regulatory elements
A question we wanted to address is how epigenetic mechanisms may intersect with changes in signalling environment of cells, and more specifically, how Pc-binding may be affected by cross-talk with transcription factors that act as effectors of signalling cascades activated during tumorigenesis. Thus, to advance our insight into how gain or loss of Pc-binding in scrib 1 WIDs may be regulated, we analysed GATC fragments classified by the three-state HMM analysis to be ‘enriched’ in Pc-binding, for predicted transcription factor binding motifs or modENCODE-identified chromatin domains [27, 46] using i-cisTarget  (see Experimental procedures). In parallel, we performed an i-cisTarget on GATC fragments classified as gain or loss of ‘enriched’ Pc-binding states in scrib 1 WIDs (Fig. 4a). As expected, Pc-bound GATC fragments in WT were enriched for PRC1 and PRC2 binding, as well as H3K27me3 and H3K9me3 modifications (Fig. 4a′, a″). In contrast, regulatory regions exhibiting dynamic Pc-binding transitions in scrib 1 displayed high NES scores for RNA-mediated silencing machineries (Piwi, Ago2), transcriptional activation by histone acetylation (Nejire/CBP) or recruitment of RNAPol II (Fig. 4a′), all of which may cooperate with CTCF (Fig. 4a′) in insulator-dependent transcriptional regulation and spatial organization of chromatin [48–52]. Interestingly, histone modifications previously observed to occur at genes that are expressed, but importantly, at intermediate levels , were also detected at dynamic Pc-binding sites (Fig. 4a″). This suggests that Pc target genes, which experience altered Pc-binding in scrib 1 , may be subject to transcriptional modulation rather than absolute repression by Pc.
Next, we wondered whether tumour-associated transcripts upregulated upon loss of Pc-binding in scrib 1 (group III, Fig. 3b) were characterized by a specific signature of regulatory elements. We thus repeated an i-cisTarget analysis for the presumptive regulatory region spanning the transcriptional start site (TSS)(− 2.5 kb ~ + 1 kb) of genes belonging to group III (Fig. 4b). Strikingly, AP-1 (Jra/Kay), Atf3, Cnc and Lola-binding motifs enriched in group III loci (Fig. 4b′, Additional file 5: Table S2) and align with the stress-dependent activation of chinmo, Atf3, Ets21C, Ilp8, upd3 and SOCS36E associated with high JNK and JAK/STAT activity during wound healing, regeneration and tumorigenesis [38, 44, 54–57].
We repeated an i-cisTarget analysis for group II genes, whose transcripts were downregulated upon gain of Pc-binding in scrib 1 (Fig. 4b) to ask how Polycomb may be recruited to these sites. In agreement with the observation that group II genes were enriched for axon guidance targets, we found that transcription factors specifically expressed in neurons, such as Jumu and CG12299, were enriched in regulatory regions of group II (Fig. 4b′, Additional file 5: Table S2). Importantly, however, wing patterning regulators, such as the transcription factor Rn and the Dpp/TGF-β signalling effectors Med and Mad, were also enriched, confirming that wing differentiation is affected in a Polycomb-dependent manner in scrib 1 WID (Fig. 4b′) . These data, however, may indicate that transcriptional downregulation of genetic circuits involved in neuronal and wing disc patterning promotes binding of Pc to these target genes.
Based on our finding that GATC fragments gaining Pc-binding in scrib 1 were enriched for CTCF (Fig. 4a′), we asked whether insulator elements locate to group II genes. Strikingly, 71% of group II genes contained Flybase-mapped class I and II insulator elements within their gene body. In contrast, insulator features mapped to only 19% of group III genes. This suggests that insulator-dependent modulation of Pc function or Pc-dependent modulation of insulator function may have important consequences for Pc-targeted gene expression in scrib 1 .
Polycomb-binding transitions fail in scrib 1 imaginal discs development
A subsequent i-cisTarget analysis of young WID profiles revealed that Pc-‘enriched’ GATC fragments in young WIDs displayed PRC1 and PRC2-binding, confirming that they are canonical Pc target sites (Fig. 5d). GATC fragments that specifically lost ‘enriched’ Pc-binding in late development scored high for binding sites of the wing differentiation regulators nubbin (Nub) and scalloped (Sd) (Fig. 5d), reflecting the expansion of the central wing domain. GATC fragments that gained Pc-binding in late development were enriched in binding sites for Atf3 and Adf1 (Fig. 5d). Adf1 was recently identified to be critical for recruitment and tethering of Pc to target sites . The enrichment of Atf3 motifs may suggest that Atf3 target genes are increasingly silenced as wing discs development progresses, which has indeed been observed for Atf3 expression . This may also have important implications for the reduction in regenerative capacity previously attributed to Pc silencing of critical tissue-stress-responsive enhancers in late WIDs .
However, GATC fragments with dynamic Pc transitions during development were also enriched for CTCF and Su(Hw) insulator components, as well as for the histone demethylase Lsd1. Combined, these invoke earlier observations of insulator signatures at dynamic Pc-targeted sites (Fig. 4a′) and imply that Pc-binding dynamics at insulator elements, which are critical for organization of chromatin in the nucleus [48–52], are crucial to Pc function during differentiation. Intriguingly, a detailed analysis of our DamID profiles revealed that the Pc-bound GATC fragment sizes recovered from earlier developmental stages were larger than those recovered from late imaginal discs (Additional file 1: Fig. S5). Moreover, in scrib 1 datasets, GATC fragment sizes occupied an intermediate distribution (Additional file 1: Fig. S5). The size range differences cannot be recapitulated by Dam profiles alone (data not shown). It may suggest that Pc-binding to genome regions characterized by different GATC motif frequencies is developmentally regulated and may reflect different distributions at promoters, introns or intergenic regions. However, it may also suggest a link between changes to Pc-binding and chromatin accessibility, where chromatin compaction during development may reduce the likelihood of distant GATC motifs to be methylated by Pc-Dam fusion proteins.
As a consequence of the limited availability and accessibility of sample material, in vivo ChIP protocols are technically challenging . Here, we report that DamID sensitively and reproducibly detects Pc-binding differences in wing imaginal discs with input samples derived from just 10 larvae. We propose that the lower limit necessary for good quality DamID profiles of imaginal discs is even less. For example, we specifically omitted PCR amplifications during preparation of NGS libraries to avoid oversampling of PCR biases. Consequently, we eliminated an opportunity to amplify weak signals to detectable levels. Indeed, published DamID-seq protocols report PCR amplification of NGS libraries without adverse effects [5, 8].
By targeting an ectopic signature to specific cells, FRT/FLP-out DamID circumvents the challenges of in vivo ChIP approaches that require the researcher to purify cell-type-specific nuclei from complex tissues. For this purpose, previously described cell-type-specific DamID systems rely either on the real-time expression patterns of GAL4 drivers (TaDa) or on cell-type-specific promoters that directly drive the expression of a FLP to achieve cell-type specificity [5–8]. In contrast, we describe a cell lineage-specific DamID system by utilizing a specific UAS-FLP(EXEL) that can be combined with any GAL4 driver for maximum flexibility to permanently target DamID to different cell types and their descendants. Genetic strategies based on individual GAL4 drivers can be optimized and validated by G-trace analysis to reveal temporal and spatial patterns of the GAL4-targeted lineage. Combined, the approach reported here opens the opportunity to track transitions of DNA binding sites in parent and daughter cell populations of a cell lineage over time.
Here we demonstrate that DamID sensitively detects significant changes in Pc-binding between three different source samples. While Pc silencing is not globally altered in a scrib 1 mutant background, the transcriptional changes that correlated with altered Pc-binding at specific loci allowed us to identify three novel Pc target genes (Atf3, Ets21C, Ilp8), which are implicated in tissues stress responses and tumour growth in many proliferating tissues [36–38, 43, 44]. We find that Atf3, AP-1 (Jra/Kay) and Lola-binding sites are enriched at these genes that are activated in scrib 1 mutant discs, suggesting that these transcriptional regulators [38, 44, 54, 55] may oppose Pc silencing to activate a PcG target gene network in tissue repair and tumorigenic transformation. Curiously, transcript levels of core PcG components are downregulated by stress-induced JNK signalling  and two core PRC1 transcripts are mildly reduced in scrib 1 WIDs . This downregulation of PcG may sensitize Pc target genes, such as Atf3, Ets21C, Ilp8, upd3, SOCS36E and chinmo, for activation in stress-induced or tumorigenic contexts.
Our findings furthermore imply the high correlation between scrib 1 and younger WID profiles indicates that a failure of scrib 1 WID to undergo Pc-dependent fate differentiation contribute to scrib 1 phenotypes. Our analysis furthermore implies that such developmental transitions mediated by Pc may be associated with insulator dynamics that could mediate global changes to accessibility of Pc-regulated chromatin domains. How insulator dynamics may regulate dynamic Pc-binding during development needs to be clarified in future studies. Similarly, while our analysis focused on Pc dynamics in different tissue states, a recent study highlights large scale remodelling of HP1-dependent chromatin and of silent ‘black’ chromatin states in developmental transitions of neuron, which are also likely to play a role in imaginal disc development and tumorigenesis .
All stocks and experimental crosses were maintained on standard fly food at 18 °C or 25 °C unless otherwise specified. The following transgenes and fly lines were used in this study:
ptcGAL4 and ptcGAL4, tubGAL80 ts /CyO;
rnGAL4; and rn[GAL4-DeltaS], tubGAL80 ts /TM6c
en-GAL4, UAS-GFP; tub-GAL80 ts
scrib 1 ;
hsflp 122 ;
Organismal induction of DamID constructs
Development of embryos was synchronized by an 8-h egg collection on standard fly food at 21 °C. FLP expression, which was controlled by a heat-shock promoter (hsflp), was induced by a 1-h temperature shift to 37 °C in a water bath. To analyse DamID profiles of young WIDs, a heat shock was performed at 3 days after egg lay (AEL). To analyse DamID profiles of late WIDs, a heat shock was performed at 5 days AEL. To account for the developmental delay characteristic of scrib 1 homozygous animals, scrib 1 larvae were heat-shocked at 6 days AEL. Afterwards, larvae were kept at 21 °C to maintain a low basal activity of the Hsp70 heat-shock promoter driving expression of Dam and Dam-Pc transcripts. Wing imaginal discs were dissected 48 h after induction of hsflp. Genomic excision of the STOP cassette from DamID constructs as a result of FLP activity was tested with regular PCR protocols on gDNA extracted from WIDs (see below) using the primers hhsp-int (actgcaactactgaaatctgc) and Dam-r (cgctattgatatcggcaagg).
Tissue dissection and genomic DNA extraction
Ten Drosophila larvae were dissected in cold Shields and Sang M3 medium, and WIDs were collected in 1.5-ml tubes on ice. Discs were resuspended in a total volume of 400 µl lysis buffer (10 mM Tris–HCl pH 8.0; 10 mM EDTA pH 8.0; 100 mM NaCl; 0.5% SDS) with proteinase K (20 mg/ml, NEB) and incubated for 4 h at 55 °C. Phenol–chloroform purification and RNase A (QIAGEN) digestion were followed up by a standard ethanol precipitation to obtain pure DNA. Each sample was subsequently run on 1% agarose gel to confirm DNA integrity and to estimate DNA concentrations. DNA from control and experimental samples was isolated at the same time and processed in parallel.
DamID sample processing, PCR and NGS library preparation
Isolation of genomic DNA (gDNA) from WIDs is described above. For each condition and stage, two independent biological samples were processed and analysed as described in  with minor changes. Briefly, after gDNA extraction, 600 µg of gDNA was digested with DpnI restriction enzyme (10 U, New England Biolabs) with CutSmart buffer (New England Biolabs) in a total volume of 10 µl at 37 °C for 6 h. DpnI digestion was terminated with heat inactivation at 80 °C for 20 min. Digested fragments were ligated to 12.5 pmol DamID adapters with T4 ligase (Roche) with T4 ligase buffer in a total volume of 20 µl for 16 h at 16 °C. Ligated gDNA fragments were subsequently digested with DpnII (10 U, New England Biolabs) in DpnII buffer (New England Biolabs) in a total volume of 50 µl for 1 h at 37 °C. Ten microlitres of DpnII digested products was amplified by PCR using MyTaq Red Mix (Bioline) with 50 µM Adr-PCR primers in a total volume of 50 µl. PCR program: 10 min at 68 °C; 1 min at 94 °C, 5 min at 65 °C, 15 min at 68 °C; 1 min at 94 °C, 1 min at 65 °C, 10 min at 68 °C—repeated 3X; 1 min at 94 °C, 1 min at 65 °C, 2 min at 68 °C—repeated (17X). Twelve microlitres of PCR products was run on 1.5% agarose gel to examine the expected DNA smear. Primers and adaptors sequences are described in . PCR products were purified using QIAquick PCR purification kit (QIAGEN) according to manufactures protocol. Samples were eluted in 50 µl of nuclease-free water. After purification, DNA concentration was determined with Qubit Fluorometric Quantitation (ThermoFisher) and adjusted to 20 ng/µl for all samples prior to libraries preparation for NGS. One microgram of DNA was transferred to a microTUBE AFA Fiber Screw-Cap 6 × 16 mm (Covaris) and sheared to an average size of around 350 bp, using a Covaris M220 focused-ultrasonicator with the following settings: duty factor = 20%, peak incident power = 50 W, cycles per burst = 200, time = 55 s, temperature = 6 °C. Illumina TruSeq PCR-free LT library preparation kit (Illumina) was used to obtain DamID-seq library according to manufactures protocol. Next-generation sequencing was run on Illumina GenomeAnalyzer IIx cBot machine. fastq file analysis was performed according to methods described in .
Bioinformatic tools—general information
Bioinformatic analysis was performed using R (v. 3.4.0) (https://www.r-project.org/) and bedtools (v. 2.26.0) software (http://bedtools.readthedocs.io/en/latest/#). Analysis for enriched regulatory elements was performed using i-cisTarget (https://gbiomed.kuleuven.be/apps/lcb/i-cisTarget/index.php) .
Identification and characterization of Pc-bound target sites
DamID-seq fastq files were processed as described previously  with the following two modifications. The mapping of reads onto GATC fragments by the software ‘HTSeq-count’ was performed with a higher stringency criterion (by using the ‘intersection_strict’ instead of ‘union’ overlap resolution mode). GATC fragments showing highly discordant values between replicates were excluded from the analysis as described .
Pc-binding sites (‘bound’ targets) were identified based on Dam-normalized log2-transformed DamID-seq profiles by fitting a three-state hidden Markov model (HMM) to define ‘enriched’, ‘intermediate’ and ‘depleted’ Pc-binding states for each GATC fragment, as described previously (Additional files 2: SF2, 3: SF3, 6: SF4) [26, 32–35]. We chose a three-state model to avoid random assignment of intermediate binding to either ‘enriched’ or ‘depleted’ states . Thus, while ‘intermediate’ states could arise for any biological, genetic or technical reasons, we could distinguish them in our analysis.
As the lengths of the genomic GATC fragments (bins) are not of equal size, we used the BioHMM algorithm, a heterogeneous HMM, which takes into account the distance between adjacent bins . This algorithm was previously implemented in the Bioconductor package snapCGH [67, 68]. We adapted the BioHMM algorithm for identification of three Pc-binding states (‘enriched’, ‘intermediate’ and ‘depleted’). The R code of adapted BioHMM algorithm is provided as Additional file 7: SF1. The three-state HMM analysis outputs of each GATC fragment were compared between ‘WT’ and ‘scrib’ datasets, as well as between ‘early’ and ‘late’ development in WT, to assess the dynamics of Polycomb-binding between two samples. To maintain the directionality of differences, the result of this comparison was reported as either ‘gain’, ‘loss’ or ‘no change’ for each GATC fragment between ‘enriched’, ‘intermediate’ and ‘depleted’ HMM states.
RNA-seq and ChIP-chip data analysis
RNA-seq datasets were obtained from . Genes were selected for further analysis according to the statistical significance (adjusted p val< 0.05) and subsequently divided in upregulated and downregulated expression according to the change in transcript levels. Differential gene expression was provided as log2 of the fold change between WT, scrib 1 and Psc/Su(z)2 XL26 datasets. ChIP-chip datasets were downloaded from the modENCODE repository (http://www.modencode.org/): Pc in S2 cells (ID 3791), Pc in DmBG3 cells (ID 325), Pc in embryo (ID 3957). Sequence overlap of microarray probe sequences in ChIP-chip datasets and Dam-normalized GATC fragments in DamID-Seq datasets was analysed using bedtool intersect function. Pearson’s correlation between DamID-seq and ChIP-chip data was calculated by correlating the intensity of Dam-normalized Pc-binding at each GATC fragment in either WT or scrib 1 datasets to the intensity of Pc-binding at the corresponding microarray probe for the respective Pc ChIP-chip analysis from S2 cells, DmBG3 cells or embryo.
Transcription-associated GATC fragments (taGATCf)
Regulatory regions associated with genes differentially expressed in scrib 1 were defined as genomic regions spanning 2.5 kb upstream to 1.5 kb downstream of the transcriptional start sites (TSS) of the selected genes. Briefly, the coordinates of the regulatory regions were calculated from the TSS coordinates and the strand on which the TSS mapped on. This information was acquired from Flybase (Batch Download, http://flybase.org) (genome annotation dm6) using the FB.ID of all differentially expressed genes. Subsequently, genome coordinates of GATC fragments were converted into the appropriate genome annotation (dm3 → dm6, LiftOver tool—UCSC, https://genome.ucsc.edu/cgi-bin/hgLiftOver) and mapped to the regulatory regions using the intersect function in bedtools2 (no limitations were considered on the amount of overlap between the two coordinates’ sets). Only GATC fragments that overlapped with selected regulatory regions were defined as transcription-associated GATC fragments (taGATCf) and used for the comparative analysis of DamID-seq and RNA-seq data in wild-type and scrib 1 WIDs (Fig. 3a). Subsequently, regulatory regions mapping to upregulated or downregulated genes were further subdivided according to transitions in Pc-binding at each of their associated taGATCf (‘gain’, ‘loss’ or ‘no change’ for each GATC fragment between ‘enriched’, ‘intermediate’ and ‘depleted’ HMM states). The entire regulatory region was subsequently classified as gain in Pc-binding, if one or more taGATCf within this region ‘gained’ Pc-binding and other taGATC fragments displayed ‘no change’. Conversely, a regulatory region was classified as loss in Pc-binding, if one or more taGATCf within this region ‘lost’ Pc-binding and other taGATC fragments displayed ‘no change’. Finally, regulatory regions which contained a mix of taGATCf with both gain and loss HHM states were classified as mixed (m1-low mRNA levels, m2-high mRNA levels, Additional file 1: Fig. S3.D) and not considered in subsequent analysis.
As a result, the described method subdivides regulatory regions into the following four groups: group I (RNA—up regulated; Pc-binding—gain); group II (RNA—down regulated; Pc-binding—gain); group III (RNA—up regulated; Pc-binding—loss); and group IV (RNA—down regulated; Pc-binding—loss).
Analysis for enriched regulatory elements using i-cisTarget
We performed our i-cisTarget analysis adhering to an enrichment score threshold = 2 and rank threshold = 10,000. We defined significantly enriched motifs by setting the normalized enrichment scores (NES) > 3. For factors with multiple enriched motifs, we selected only the one with the highest NES. The following features (Databases 3.0 of i-cisTarget) were selected during the analysis: PWMs, TF binding sites, non-TF binding sites, histone modifications. These parameters were common to all icis-Target analysis.
icis-Target analysis on GATC fragments with assigned HMM transitions:
Figure 4a–a″: this analysis was performed on pools of GATC fragments with the following defined HMM transition states: gain or loss of Pc-binding in scrib 1 by transitioning in and out of ‘enriched’ HMM states, and no change in Pc-binding in scrib 1 by staying ‘enriched’ (excluding ‘depleted’ and ‘intermediate’ HMM states).
Figure 5d: one analysis was performed on GATC fragments that were defined as ‘enriched’ in WT Early profiles after the three-state HMM analysis (Fig. 5d). Another analysis was performed on pools of GATC fragments with gain or loss of Pc-binding in ‘WT Late’ discs by in and out of ‘enriched’ HMM states.
icis-Target analysis on taGATCf fragments mapping to the presumptive regulatory region of Pc-targeted genes:
Sequences of all regulatory regions established for the analysis of taGATCf were first converted to a genome annotation suitable for icis-Target analysis (dm6 → dm3) and then subdivided into their respective group (group I, group II, group III and group IV). The icis-Target analysis was performed on groups II and III independently (Fig. 4c–c″).
To detect the Myc-tagged Dam proteins, expression of Dam and Dam-Pc constructs was boosted by a heat shock for 1 h at 37 °C 6 h prior to dissection to strongly induce the Hsp70 promoter. This heat shock induces abnormally high Dam and Dam-Pc expression levels that can be detected by immunohistochemistry but are unsuitable for genomic DamID profiling and reduce cell viability. Larvae were dissected and cuticles were fixed for 15 min at room temperature in 4% paraformaldehyde (PFA). Washing steps were performed in 0.1% Triton X-100/PBS (PBT). The following antibodies were incubated overnight at 4 °C: rabbit α-Dcp-1 (1:500, Cell Signalling), mouse α-H3S10p (1:2000, Abcam), mouse α-Myc (1:50, DSHB). Secondary antibodies (Molecular Probes), DAPI and phalloidin-TRITC (Sigma) were incubated at room temperature for 2 h. Experimental and control samples were processed together and imaged on the same microscope (Leica TCS SP-5).
Adult wing imaging
Adult flies were collected 12 h after eclosion and stored in 2-propanol. Wings were dissected and mounted in Euparal (Sigma) on regular slides for microscopy. Imaging was done using a stereoscopic zoom microscope (Nikon, SMZ745).
MLF, GG and AKC designed the experiments. MLF, GG, AC, LB, SK and HB performed the experiments. MLF, GG, AP, BS and AKC analysed the data. MLF and AKC wrote the paper. All authors read and approved the final manuscript.
We thank the reviewers for critical comments on the manuscript. We thank B. Bunker, D. Bilder, T. Straub, A. Ivankin for technical help and advice with bioinformatics and data analysis. We thank the Bloomington Stock Centre and DSHB for providing fly stocks and antibodies. We thank the IMPRS-LS and LSM Munich graduate schools for supporting our students.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during the current study will be made available in the GEO repository. They are submitted for manuscript review as additional files (Additional file 8: DamID_HMM_Dm3.txt, Additional file 9: DamID_Raw_Counts_Dm3.txt) mapped on Dm3 Drosophila genome annotation.
Consent for publication
Ethics approval and consent to participate
Funding for this work was provided by the DFG (CL490-1 to AKC) and Russian Fundamental Scientific Research Program (0310-2018-0009 to AP) and the Russian Science Foundation (16-14-10288 to AVP).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Deal RB, Henikoff S. The INTACT method for cell type-specific gene expression and chromatin profiling in Arabidopsis thaliana. Nat Protoc. 2011;6(1):56–68.View ArticlePubMedGoogle Scholar
- Steiner FA, Talbert PB, Kasinathan S, Deal RB, Henikoff S. Cell-type-specific nuclei purification from whole animals for genome-wide expression and chromatin profiling. Genome Res. 2012;22(4):766–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonn S, Zinzen RP, Perez-Gonzalez A, Riddell A, Gavin AC, Furlong EE. Cell type-specific chromatin immunoprecipitation from multicellular complex samples using BiTS-ChIP. Nat Protoc. 2012;7(5):978–94.View ArticlePubMedGoogle Scholar
- Schauer T, Schwalie PC, Handley A, Margulies CE, Flicek P, Ladurner AG. CAST-ChIP maps cell-type-specific chromatin states in the Drosophila central nervous system. Cell reports. 2013;5(1):271–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Pindyurin AV, Pagie L, Kozhevnikova EN, van Arensbergen J, van Steensel B. Inducible DamID systems for genomic mapping of chromatin proteins in Drosophila. Nucleic Acids Res. 2016;44(12):5646–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Southall TD, Gold KS, Egger B, Davidson CM, Caygill EE, Marshall OJ, Brand AH. Cell-type-specific profiling of gene expression and chromatin binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev Cell. 2013;26(1):101–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Pindyurin AV. Genome-wide cell type-specific mapping of in vivo chromatin protein binding using an FLP-Inducible DamID system in Drosophila. Methods Mol Biol. 2017;1654:99–124.View ArticlePubMedGoogle Scholar
- Marshall OJ, Southall TD, Cheetham SW, Brand AH. Cell-type-specific profiling of protein-DNA interactions without cell isolation using targeted DamID with next-generation sequencing. Nat Protoc. 2016;11(9):1586–98.View ArticlePubMedGoogle Scholar
- Kind J, Pagie L, Ortabozkoyun H, Boyle S, de Vries SS, Janssen H, Amendola M, Nolen LD, Bickmore WA, van Steensel B. Single-cell dynamics of genome-nuclear lamina interactions. Cell. 2013;153(1):178–92.View ArticlePubMedGoogle Scholar
- Beira JV, Paro R. The legacy of Drosophila imaginal discs. Chromosoma. 2016;125(4):573–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Steffen PA, Ringrose L. What are memories made of? How Polycomb and Trithorax proteins mediate epigenetic memory. Nat Rev Mol Cell Biol. 2014;15(5):340–56.View ArticlePubMedGoogle Scholar
- Chittock EC, Latwiel S, Miller TC, Muller CW. Molecular architecture of polycomb repressive complexes. Biochem Soc Trans. 2017;45(1):193–205.View ArticlePubMedPubMed CentralGoogle Scholar
- Muller J, Verrijzer P. Biochemical mechanisms of gene regulation by polycomb group protein complexes. Curr Opin Genet Dev. 2009;19(2):150–8.View ArticlePubMedGoogle Scholar
- Geisler SJ, Paro R. Trithorax and Polycomb group-dependent regulation: a tale of opposing activities. Development. 2015;142(17):2876–87.View ArticlePubMedGoogle Scholar
- Sparmann A, van Lohuizen M. Polycomb silencers control cell fate, development and cancer. Nat Rev Cancer. 2006;6(11):846–56.View ArticlePubMedGoogle Scholar
- Classen AK, Bunker BD, Harvey KF, Vaccari T, Bilder D. A tumor suppressor activity of Drosophila Polycomb genes mediated by JAK-STAT signaling. Nat Genet. 2009;41(10):1150–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Bunker BD, Nellimoottil TT, Boileau RM. Classen AK. Bilder D: The transcriptional response to tumorigenic polarity loss in Drosophila. eLife; 2015. p. 4.Google Scholar
- Schwartz YB, Kahn TG, Nix DA, Li XY, Bourgon R, Biggin M, Pirrotta V. Genome-wide analysis of Polycomb targets in Drosophila melanogaster. Nat Genet. 2006;38(6):700–5.View ArticlePubMedGoogle Scholar
- Ringrose L. Polycomb comes of age: genome-wide profiling of target sites. Curr Opin Cell Biol. 2007;19(3):290–7.View ArticlePubMedGoogle Scholar
- Schwartz YB, Pirrotta V. Polycomb silencing mechanisms and the management of genomic programmes. Nat Rev Genet. 2007;8(1):9–22.View ArticlePubMedGoogle Scholar
- Duffy JB, Harrison DA, Perrimon N. Identifying loci required for follicular patterning using directed mosaics. Development. 1998;125(12):2263–71.PubMedGoogle Scholar
- Thibault ST, Singer MA, Miyazaki WY, Milash B, Dompe NA, Singh CM, Buchholz R, Demsky M, Fawcett R, Francis-Lang HL, Ryner L, Cheung LM, Chong A, Erickson C, Fisher WW, Greer K, Hartouni SR, Howie E, Jakkula L, Joo D, Killpack K, Laufer A, Mazzotta J, Smith RD, Stevens LM, Stuber C, Tan LR, Ventura R, Woo A, Zakrajsek I, Zhao L, Chen F, Swimmer C, Kopczynski C, Duyk G, Winberg ML, Margolis J. A complementary transposon tool kit for Drosophila melanogaster using P and piggyBac. Nat. Genet. 2004;36(3):283–7.View ArticlePubMedGoogle Scholar
- Evans CJ, Olson JM, Ngo KT, Kim E, Lee NE, Kuoy E, Patananan AN, Sitz D, Tran P, Do MT, et al. G-TRACE: rapid Gal4-based cell lineage analysis in Drosophila. Nat Methods. 2009;6(8):603–5.View ArticlePubMedPubMed CentralGoogle Scholar
- van Steensel B, Delrow J, Henikoff S. Chromatin profiling using targeted DNA adenine methyltransferase. Nat Genet. 2001;27(3):304–8.View ArticlePubMedGoogle Scholar
- Hariharan IK, Bilder D. Regulation of imaginal disc growth by tumor-suppressor genes in Drosophila. Annu Rev Genet. 2006;40:335–61.View ArticlePubMedGoogle Scholar
- van Bemmel JG, Filion GJ, Rosado A, Talhout W, de Haas M, van Welsem T, van Leeuwen F, van Steensel B. A network model of the molecular organization of chromatin in Drosophila. Mol Cell. 2013;49(4):759–71.View ArticlePubMedGoogle Scholar
- Negre N, Brown CD, Ma L, Bristow CA, Miller SW, Wagner U, Kheradpour P, Eaton ML, Loriaux P, Sealfon R, et al. A cis-regulatory map of the Drosophila genome. Nature. 2011;471(7339):527–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Kharchenko PV, Alekseyenko AA, Schwartz YB, Minoda A, Riddle NC, Ernst J, Sabo PJ, Larschan E, Gorchakov AA, Gu T, et al. Comprehensive analysis of the chromatin landscape in Drosophila melanogaster. Nature. 2011;471(7339):480–5.View ArticlePubMedGoogle Scholar
- Shimbo T, Du Y, Grimm SA, Dhasarathy A, Mav D, Shah RR, Shi H, Wade PA. MBD3 localizes at promoters, gene bodies and enhancers of active genes. PLoS Genet. 2013;9(12):e1004028.View ArticlePubMedPubMed CentralGoogle Scholar
- Negre N, Hennetin J, Sun LV, Lavrov S, Bellis M, White KP, Cavalli G. Chromosomal distribution of PcG proteins during Drosophila development. PLoS Biol. 2006;4(6):e170.View ArticlePubMedPubMed CentralGoogle Scholar
- Tolhuis B, de Wit E, Muijrers I, Teunissen H, Talhout W, van Steensel B, van Lohuizen M. Genome-wide profiling of PRC1 and PRC2 Polycomb chromatin binding in Drosophila melanogaster. Nat Genet. 2006;38(6):694–9.View ArticlePubMedGoogle Scholar
- Filion GJ, van Bemmel JG, Braunschweig U, Talhout W, Kind J, Ward LD, Brugman W, de Castro IJ, Kerkhoven RM, Bussemaker HJ, et al. Systematic protein location mapping reveals five principal chromatin types in Drosophila cells. Cell. 2010;143(2):212–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Meuleman W, Peric-Hupkes D, Kind J, Beaudry JB, Pagie L, Kellis M, Reinders M, Wessels L, van Steensel B. Constitutive nuclear lamina-genome interactions are highly conserved and associated with A/T-rich sequence. Genome Res. 2013;23(2):270–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Babu MM. Computational approaches to study transcriptional regulation. Biochem Soc Trans. 2008;36(Pt 4):758–65.View ArticlePubMedGoogle Scholar
- Ali M, Rincon-Arano H, Zhao W, Rothbart SB, Tong Q, Parkhurst SM, Strahl BD, Deng LW, Groudine M, Kutateladze TG. Molecular basis for chromatin binding and regulation of MLL5. Proc Natl Acad Sci USA. 2013;110(28):11296–301.View ArticlePubMedPubMed CentralGoogle Scholar
- Toggweiler J, Willecke M, Basler K. The transcription factor Ets21C drives tumor growth by cooperating with AP-1. Sci Rep. 2016;6:34725.View ArticlePubMedPubMed CentralGoogle Scholar
- Kulshammer E, Mundorf J, Kilinc M, Frommolt P, Wagle P, Uhlirova M. Interplay among Drosophila transcription factors Ets21c, Fos and Ftz-F1 drives JNK-mediated tumor malignancy. Dis Model Mech. 2015;8(10):1279–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Atkins M, Potier D, Romanelli L, Jacobs J, Mach J, Hamaratoglu F, Aerts S, Halder G. An ectopic network of transcription factors regulated by hippo signaling drives growth and invasion of a malignant tumor model. Curr Biol. 2016;26(16):2101–13.View ArticlePubMedGoogle Scholar
- Colombani J, Andersen DS, Leopold P. Secreted peptide Dilp8 coordinates Drosophila tissue growth with developmental timing. Science. 2012;336(6081):582–5.View ArticlePubMedGoogle Scholar
- Garelli A, Gontijo AM, Miguela V, Caparros E, Dominguez M. Imaginal discs secrete insulin-like peptide 8 to mediate plasticity of growth and maturation. Science. 2012;336(6081):579–82.View ArticlePubMedGoogle Scholar
- Amoyel M, Anderson AM, Bach EA. JAK/STAT pathway dysregulation in tumors: a Drosophila perspective. Semin Cell Dev Biol. 2014;28:96–103.View ArticlePubMedGoogle Scholar
- Loubiere V, Delest A, Thomas A, Bonev B, Schuettengruber B, Sati S, Martinez AM, Cavalli G. Coordinate redeployment of PRC1 proteins suppresses tumor formation during Drosophila development. Nat Genet. 2016;48(11):1436–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Flaherty MS, Salis P, Evans CJ, Ekas LA, Marouf A, Zavadil J, Banerjee U, Bach EA. chinmo is a functional effector of the JAK/STAT pathway that regulates eye development, tumor formation, and stem cell self-renewal in Drosophila. Dev Cell. 2010;18(4):556–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Davie K, Jacobs J, Atkins M, Potier D, Christiaens V, Halder G, Aerts S. Discovery of transcription factors and regulatory regions driving in vivo tumor development by ATAC-seq and FAIRE-seq open chromatin profiling. PLoS Genet. 2015;11(2):e1004994.View ArticlePubMedPubMed CentralGoogle Scholar
- Santabarbara-Ruiz P, Lopez-Santillan M, Martinez-Rodriguez I, Binagui-Casas A, Perez L, Milan M, Corominas M, Serras F. ROS-induced JNK and p38 signaling is required for unpaired cytokine activation during drosophila regeneration. PLoS Genet. 2015;11(10):e1005595.View ArticlePubMedPubMed CentralGoogle Scholar
- Roy S, Ernst J, Kharchenko PV, Kheradpour P, Negre N, Eaton ML, Landolin JM, Bristow CA, Ma L, Lin MF, et al. Identification of functional elements and regulatory circuits by Drosophila modENCODE. Science. 2010;330(6012):1787–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Verfaillie A, Imrichova H, Janky R, Aerts S. iRegulon and i-cisTarget: reconstructing regulatory networks using motif and track enrichment. Curr Protoc Bioinformatics. 2015;52(2):1611–39.Google Scholar
- Moshkovich N, Nisha P, Boyle PJ, Thompson BA, Dale RK, Lei EP. RNAi-independent role for Argonaute2 in CTCF/CP190 chromatin insulator function. Genes Dev. 2011;25(16):1686–701.View ArticlePubMedPubMed CentralGoogle Scholar
- Philip P, Boija A, Vaid R, Churcher AM, Meyers DJ, Cole PA, Mannervik M, Stenberg P. CBP binding outside of promoters and enhancers in Drosophila melanogaster. Epigenetics Chromatin. 2015;8:48.View ArticlePubMedPubMed CentralGoogle Scholar
- Cubenas-Potts C, Corces VG. Architectural proteins, transcription, and the three-dimensional organization of the genome. FEBS Lett. 2015;589(20 Pt A):2923–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Li HB, Muller M, Bahechar IA, Kyrchanova O, Ohno K, Georgiev P, Pirrotta V. Insulators, not Polycomb response elements, are required for long-range interactions between Polycomb targets in Drosophila melanogaster. Mol Cell Biol. 2011;31(4):616–25.View ArticlePubMedGoogle Scholar
- Li HB, Ohno K, Gui H, Pirrotta V. Insulators target active genes to transcription factories and polycomb-repressed genes to polycomb bodies. PLoS Genet. 2013;9(4):e1003436.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosenfeld JA, Wang Z, Schones DE, Zhao K, DeSalle R, Zhang MQ. Determination of enriched histone modifications in non-genic portions of the human genome. BMC Genom. 2009;10:143.View ArticleGoogle Scholar
- Ferres-Marco D, Gutierrez-Garcia I, Vallejo DM, Bolivar J, Gutierrez-Avino FJ, Dominguez M. Epigenetic silencers and Notch collaborate to promote malignant tumours by Rb silencing. Nature. 2006;439(7075):430–6.View ArticlePubMedGoogle Scholar
- Blanco E, Ruiz-Romero M, Beltran S, Bosch M, Punset A, Serras F, Corominas M. Gene expression following induction of regeneration in Drosophila wing imaginal discs. Expression profile of regenerating wing discs. BMC Dev Biol. 2010;10:94.View ArticlePubMedPubMed CentralGoogle Scholar
- Doggett K, Turkel N, Willoughby LF, Ellul J, Murray MJ, Richardson HE, Brumby AM. BTB-zinc finger oncogenes are required for ras and notch-driven tumorigenesis in Drosophila. PLoS ONE. 2015;10(7):e0132987.View ArticlePubMedPubMed CentralGoogle Scholar
- Brock AR, Seto M, Smith-Bolton RK. Cap-n-collar promotes tissue regeneration by regulating ROS and JNK signaling in the Drosophila melanogaster wing imaginal Disc. Genetics. 2017;206(3):1505–20.View ArticlePubMedGoogle Scholar
- Orsi GA, Kasinathan S, Hughes KT, Saminadin-Peter S, Henikoff S, Ahmad K. High-resolution mapping defines the cooperative architecture of Polycomb response elements. Genome Res. 2014;24(5):809–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Sekyrova P, Bohmann D, Jindra M, Uhlirova M. Interaction between Drosophila bZIP proteins Atf3 and Jun prevents replacement of epithelial cells during metamorphosis. Development. 2010;137(1):141–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Harris RE, Setiawan L, Saul J, Hariharan IK. Localized epigenetic silencing of a damage-activated WNT enhancer limits regeneration in mature Drosophila imaginal discs. Elife. 2016;. https://doi.org/10.7554/eLife.11588.Google Scholar
- Papp B, Muller J. Histone trimethylation and the maintenance of transcriptional ON and OFF states by trxG and PcG proteins. Genes Dev. 2006;20(15):2041–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee N, Maurange C, Ringrose L, Paro R. Suppression of Polycomb group proteins by JNK signalling induces transdetermination in Drosophila imaginal discs. Nature. 2005;438(7065):234–7.View ArticlePubMedGoogle Scholar
- Marshall OJ, Brand AH. Chromatin state changes during neural development revealed by in vivo cell-type specific profiling. Nat Commun. 2017;8(1):2271.View ArticlePubMedPubMed CentralGoogle Scholar
- Vogel MJ, Peric-Hupkes D, van Steensel B. Detection of in vivo protein-DNA interactions using DamID in mammalian cells. Nat Protoc. 2007;2(6):1467–78.View ArticlePubMedGoogle Scholar
- Maksimov DA, Laktionov PP, Belyakin SN. Data analysis algorithm for DamID-seq profiling of chromatin proteins in Drosophila melanogaster. Chromosome Res. 2016;24(4):481–94.View ArticlePubMedGoogle Scholar
- Marioni JC, Thorne NP, Tavare S. BioHMM: a heterogeneous hidden Markov model for segmenting array CGH data. Bioinformatics. 2006;22(9):1144–6.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith ML, Marioni JC, Hardcastle TJ, Thorne NP. snapCGH: Segmentation. Bioconductor: Normalization and processing of aCGH data users’ guide; 2006.Google Scholar