Materials and methods
Plasmids
pHes5, a 764bp fragment of the mouse Hes5-promoter including the 5’UTR (Ohtsuka et al.,
2006), pTa1, a 1097bp fragment of the mouse Tuba1a promoter (Gal et al., 2006; Gloster et
al., 1994; Mizutani et al., 2007), pNeuroD1-Cre and pCAG-Venus were described before
(van den Ameele et al., 2022).
Dam-fusion proteins: The H3K4me3 and H3K27me3 chromatin binding proteins (TAF3
PHD amino acids 856-929 (Kungulovski et al., 2016), the DNMT3A PWWP domain and
CBX7 Chromo Domain amino acids 7-61 (Kungulovski et al., 2014)) were amplified by RT-
PCR from cDNA of human ESC-derived NSCs. TAF3 PHD was cloned C-terminal, and
CBX7 CD N-terminal of Dam. The single chain antibodies against H3K9ac (19e5; gift of Y.
Sato (Sato et al., 2013)), H3K9me3 (JH9-6 from clone 309M3-A; gift of S. Koide (Hattori et
al., 2013)) and H4K20me1 (15f11; gift of Y. Sato (Sato et al., 2016)) were amplified by PCR
from plasmids and cloned N-terminal of Dam. A nuclear localization signal (2xSV40-NLS)
was inserted between Dam and all mintbodies or CBX7 CD and Dam. When this NLS was
omitted from 19e5-Dam or placed N-terminal of 19e5, methylation efficiency in Drosophila
NSCs was reduced (data not shown).
Intron-dam constructs: i1Dam and i2Dam were described previously (Wade et al., 2021;van
den Ameele et al., 2022). I1Dam is intron 3 of mouse IghE (Lacy-Hulbert et al., 2001),
placed between the 3rd and 4th helix (BamHI site) of the DNA-binding domain of the Dam
methylase (Horton et al., 2006). I2Dam is a modified version of the Promega chimeric intron
sequence (Bothwell et al., 1981; Promega, 2009; Senapathy et al., 1990). Modifications to the
intron and the exon sequence at the splice junctions of i2Dam were made to optimize in silico
predicted splicing efficiency, remove a weak predicted acceptor site and insert an XmaI site.
Splice site predictions were performed with NNSPLICE0.9 (Reese et al., 1997) and
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Spliceport (Dogan et al., 2007). Both constructs were cloned into pCAG-IRES-GFP (pCIG,
gift from P. Vanderhaeghen) to obtain pCAG-mcherry-i1/2Dam. Introducing an intron into
the Dam CDS has the additional benefit that it greatly improves cloning efficiency, by
preventing toxicity in bacteria.
FloxDam construct: floxDam was described previously (van den Ameele et al., 2022) and
contains Lox71 and Lox61 sites (Albert et al., 1995) respectively upstream of mCherry and
within intron2. The two Lox-sites and intervening sequence are inverted to obtain pCAG-
flox2Dam. Cre recombinase activity will result in unidirectional reversion of this cassette,
thus reconstituting the i2Dam construct. Placing the Lox-site within the intron has the
advantage that splicing will remove the remaining lox-site from the Dam ORF, which
otherwise interferes with methylation efficiency (data not shown).
Xenopus plasmids: Dam and Dam-TAF3phd were cloned into the EcoRI-NotI sites of pCS2-
HA+. mRNA was synthetized in vitro using the MEGAscript SP6 Kit (Ambion, AM1330M)
following the manufacturer’s instructions.
Drosophila plasmids: pUAST-attB-mCherry-i2Dam was cloned by replacing Dam in
pUAST-attB-LT3-Ndam (Southall et al., 2013) with i2Dam from pHes5-mCherry-i2Dam. All
Dam-fusion proteins in Drosophila were fused to Dam, apart from the CBX7 Chromo
Domain, which was fused to i2Dam.
Dam-positive bacteria were 10-beta Competent E. coli (NEB C3019H). Dam-negative
bacteria were dam-/dcm- competent E. coli (NEB C2925H). Plasmids for IUE and transient
transfection DamID were prepared from Dam-negative bacteria with Endofree Plasmid Maxi
kit (Qiagen 12362).
Mice and in utero electroporation (IUE)
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
All mouse husbandry and experiments were carried out in a Home Office-designated facility,
according to the UK Home Office guidelines upon approval by the local ethics committee
(Project Licence PPL70/8727). Experiments were done in wild-type MF1 mice. Timed
natural matings were used, where noon of the day of plug-identification was E0.5. IUE was
performed as previously described (Dimidschstein et al., 2013; Saito and Nakatsuji, 2001) at
E13.5 with 50ms, 40V unipolar pulses (BTX ECM830) using CUY650P5 electrodes
(Sonidel). DamID plasmids were injected at 1µg/µl together with pCAG-Venus at 0.25µg/µl.
Embryos were harvested after 24 hours (E14.5) for TaDa with pHes5 and pTa1 or after 72
hours (E16.5) for TaDa with pND1-Cre.
Xenopus husbandry and mRNA injection
Mature Xenopus laevis males and females were obtained from Nasco. Work with Xenopus
laevis is covered under the Home Office Project Licence PPL70/8591. Frog husbandry and
all experiments were performed according to the relevant regulatory standards, essentially as
previously described (Hörmanseder et al., 2017). For mRNA injection, eggs were in vitro
fertilized, dejellied using 2% Cystein solution in 0.1xMMR, pH 7.8, washed 3 times with
0.1x MMR and transferred into 0.5x MMR for injections. 2-3 pg of mRNA was used per
injection. Embryos were cultured at 23°C and collected at gastrula stage 11 or neurula stage
21. Genomic DNA extraction for DamID library preparation was performed using a Qiagen
DNAeasy blood and tissue kit, as per manufacturer instructions.
Fly husbandry and transgenesis
UAS-mCherry-i2Dam was injected into y,sc,v,nos-phiC31;attP40;+ (Bl#25709) and
y,sc,v,nos-phiC31;+;attP2 (Bl#25710) embryos. All other Drosophila constructs were only
injected into y,sc,v,nos-phiC31;attP40;+ (Bl#25709) embryos. Male transgenic TaDa flies
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
were crossed to w1118;Worniu-GAL4 (from (Albertson, 2004));Tub-Gal80ts (active in NSCs)
virgins. Flies were reared in cages at 25°C, embryos were collected on food plates for 3 hr
and transferred to 29°C for 24 hours before dissection at third instar.
DamID-seq
For Drosophila DamID, between 20 and 50 larval brains were dissected. For DamID on
mouse cortex after IUE, embryos were cooled on ice, and the electroporated cortex was
identified with a fluorescent binocular microscope (Fig. 2B). Meninges were removed and
the electroporated region microdissected. For Xenopus DamID, 5 injected embryos were
pooled for genomic DNA extraction, and one fifth of the genomic DNA was processed for
DamID. All samples were processed for DamID as described previously (Marshall et al.,
2016). DamID fragments were prepared for Illumina sequencing according to a modified
TruSeq protocol (Marshall et al., 2016). All sequencing was performed as single end 50 bp
reads generated by the Gurdon Institute NGS Core using an Illumina HiSeq 1500.
DamID-seq data processing
Analysis of fastq-files from DamID-experiments was performed with the damidseq pipeline
script (Marshall and Brand, 2015) that maps reads to an indexed bowtie2 genome, bins into
GATC-fragments according to GATC-sites and normalises reads against a Dam-only control.
Preprocessed fastq-files from Drosophila, Xenopus and mouse were mapped respectively to
the mm10 (GRCm38.p6), dm6 (BDGP6) or Xla.v9.1 genome assemblies. Reads were binned
into fragments delineated by 5’-GATC-3’ motifs (GATC-bins). Individual replicates (see Fig.
S1 for n-numbers) for the Dam-fusion constructs were normalized against separate Dam-only
replicates with a modified version of the damidseq_pipeline (Marshall and Brand, 2015)
(RPM normalization, 300 bp binsize) and all resulting binding profiles for one Dam-fusion
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
construct were quantile normalized to each other (Marshall and Brand, 2015). The resulting
logarithmic profiles in bedgraph format were averaged for all GATC-bins across the genome
and subsequently backtransformed (“unlog”). Files were converted to the bigwig file format
with bedGraphToBigWig (v4) for visualization with the Integrative Genomics Viewer (IGV)
(v2.4.19).
Peak calling and analysis
Macs2 (v2.1.2) (Zhang et al., 2008) was used to call broad peaks for every dam-fusion/dam-
only pair on the set of *.bam-files generated by the damidseq_pipeline, using Dam-only as
control. Peaks were filtered stringently for FDR<10-2 (mouse) or FDR<10-5 (Drosophila and
Xenopus) and were only considered if present in all pairwise comparisons for a particular
experimental condition. Peaks were merged when they were within 1kb from each other with
the merge-function from bedtools (v2.26.0). A similar approach was used for the published
ChIP-seq datasets, without control. Enrichment of publicly available histone modification
datasets from ENCODE or modENCODE for mouse peak coordinates was determined with i-
cisTarget (Imrichová et al., 2015). Overlapping features were detected by uploading peak
coordinates to the i-cisTarget online platform (Gene annotation: RefSeq r70; Database: v5.0)
after conversion of mouse mm10 peak coordinates to mm9 with the online UCSC liftOver
tool. Normalized Enrichment Scores were plotted as bar graphs in Microsoft Excel.
Associating peaks with the defined genomic features was performed using ChIPseeker
(v3.10), with plotAnnoBar function and promoter defined as 2kb upstream of the
transcription start site.
Genome wide correlation
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Correlation of signal intensity between TaDa conditions and with ChIP-seq datasets from
mouse E14.5 forebrain (ENCODE) across the genome was done with the
multiBigwigSummary and plotCorrelation functions from the deepTools suite (v3.1.3).
Default bin size of 10kb was used, which allows smoothing of the signal and accounts for
variability in peak centre location between TaDa and ChIP-seq. For genome-wide correlation
between replicates a bin size of 1kb (-bs 1000) was used. Spearman’s correlation coefficient
was calculated to account for large differences in signal intensity between TaDa and ChIP-
seq.
Gene coordinates and gene expression data
Drosophila NSC gene expression came from (Yang et al., 2016). For genes expressed in
NSCs, all genes with TPM>10 in all three replicates from sorted antennal lobe NSCs were
included; for genes not expressed in NSCs, all genes with 0≤TPM<1 in all three replicates
were included. Mouse RGC gene expression came from (Florio et al., 2015). The top or
bottom 1500 genes from the aRG bulk RNA seq dataset (average of 4 replicates) were used
for further analysis. Mouse gene coordinates of protein-coding transcripts defined by ensembl
were downloaded from biomaRt (v2.38.0). Drosophila gene coordinates of “canonical genes”
were obtained from the UCSC Table Browser (Karolchik et al., 2004).
Data visualization
Genome browser views were generated using the Integrative Genomics Viewer IGV
(v2.4.19) (Robinson et al., 2011) with the midline for TaDa ratio tracks set at 1 and for ChIP-
seq set at 0. Peak or gene coordinates were saved in .bed format and supplied as features to
Seqplots (v1.12.1) (Stempor and Ahringer, 2016); quantile normalized, averaged and
backtransformed TaDa profiles were provided in bigwig format. plotAverage and
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
plotHeatmap were used to visualize and average the binding intensities across all supplied
coordinates with bin-size of 50bp. Clustering of histone mark signal over gene bodies
(including 5kb up- and downstream) was done with the built-in clustering tool of SeqPlots
(v1.12.1) (Stempor and Ahringer, 2016), using k-means clustering into 10 separate clusters.
Figures were assembled in Adobe Illustrator.
Previously published datasets
In situ hybridization images were from the Eurexpress database (Diez-Roux et al., 2011).
Histone mark ChIP-seq datasets from mouse E14.5 forebrain (ENCFF002AME for
H3K4me3; ENCFF002ANR for H3K9ac; ENCFF002DQB for H3K9me3; ENCFF002AFY
for H3K27me3) was downloaded respectively from the ENCODE (Bernstein et al., 2012)
portal as fastq-files and remapped to the mm10 genome assemblies using the just_align
function of the damidseq pipeline (Marshall and Brand, 2015). Drosophila NSC RNA PolII
binding data came from GSE77860 (Marshall et al., 2017). Mouse RGC gene expression
came from GSE65000 (Florio et al., 2015).
Data availability
All sequencing data is available from GEO.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
References
Albert, H., Dale, E. C. C., Lee, E. and Ow, D. W. W. (1995). Site‐specific integration of
DNA into wild‐type and mutant lox sites placed in the plant genome. Plant J.
Albertson, R. (2004). Scribble protein domain mapping reveals a multistep localization
mechanism and domains necessary for establishing cortical polarity. J. Cell Sci. 117,
6061–6070.
Bannister, A. J. and Kouzarides, T. (2011). Regulation of chromatin by histone
modifications. Cell Res. 21(3), 381-395.
Bernstein, B. E., Birney, E., Dunham, I., Green, E. D., Gunter, C., Snyder, M. and
ENCODE Project Consortium (2012). An integrated encyclopedia of DNA elements
in the human genome. Nature 489, 57–74.
Binek, A., Rojo, D., Godzien, J., Rupérez, F. J., Nuñez, V., Jorge, I., Ricote, M.,
Vázquez, J. and Barbas, C. (2018). Flow cytometry has a significant impact on the
cellular metabolome. J. Proteome Res. acs.jproteome.8b00472.
Bothwell, A. L. M. L. M., Paskind, M., Reth, M., Imanishi-Kari, T., Rajewsky, K. and
Baltimore, D. (1981). Heavy chain variable region contribution to the NPb family of
antibodies: somatic mutation evident in a γ2a variable region. Cell. 24(3), 625-637.
Brand, A. H. and Perrimon, N. (1993). Targeted gene expression as a means of altering cell
fates and generating dominant phenotypes. Development 118, 401–15.
Celniker, S. E., Dillon, L. A. L., Gerstein, M. B., Gunsalus, K. C., Henikoff, S., Karpen,
G. H., Kellis, M., Lai, E. C., Lieb, J. D., MacAlpine, D. M., et al. (2009). Unlocking
the secrets of the genome. Nature. 459, 927-930.
Cheetham, S. W., Gruhn, W. H., van den Ameele, J., Krautz, R., Southall, T. D.,
Kobayashi, T., Surani, M. A. and Brand, A. H. (2018). Targeted DamID reveals
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
differential binding of mammalian pluripotency factors. Development 145, dev.170209.
Chirichella, M., Lisi, S., Fantini, M., Goracci, M., Calvello, M., Brandi, R., Arisi, I.,
D’Onofrio, M., Di Primio, C. and Cattaneo, A. (2017). Post-translational selective
intracellular silencing of acetylated proteins with de novo selected intrabodies. Nat.
Methods. 14(3), 279-282.
Diez-Roux, G., Banfi, S., Sultan, M., Geffers, L., Anand, S., Rozado, D., Magen, A.,
Canidio, E., Pagani, M., Peluso, I., et al. (2011). A high-resolution anatomical atlas of
the transcriptome in the mouse embryo. PLoS Biol. journal.pbio.1000582.
Dimidschstein, J., Passante, L., Dufour, A., van den Ameele, J., Tiberi, L., Hrechdakian,
T., Adams, R., Klein, R., Lie, D. C., Jossin, Y., et al. (2013). Ephrin-B1 controls the
columnar distribution of cortical pyramidal neurons by restricting their tangential
migration. Neuron 79, 1123–1135.
Dogan, R. I., Getoor, L., Wilbur, W. J. and Mount, S. M. (2007). SplicePort-An
interactive splice-site analysis tool. Nucleic Acids Res. 35, W285-91.
Florio, M., Albert, M., Taverna, E., Namba, T., Brandl, H., Lewitus, E., Haffner, C.,
Sykes, A., Wong, F. K., Peters, J., et al. (2015). Human-specific gene ARHGAP11B
promotes basal progenitor amplification and neocortex expansion. Science. 347(6229),
1465-70.
Gal, J. S., Morozov, Y. M., Ayoub, A. E., Chatterjee, M., Rakic, P. and Haydar, T. F.
(2006). Molecular and morphological heterogeneity of neural precursors in the mouse
neocortical proliferative zones. J Neurosci 26, 1045–1056.
Gloster, a, Wu, W., Speelman, A., Weiss, S., Causing, C., Pozniak, C., Reynolds, B.,
Chang, E., Toma, J. G. and Miller, F. D. (1994). The T alpha 1 alpha-tubulin
promoter specifies gene expression as a function of neuronal growth and regeneration in
transgenic mice. J. Neurosci. 14(12), 7319-30.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Gurdon, J. B. (1960). The developmental capacity of nuclei taken from differentiating
endoderm cells of Xenopus laevis. J. Embryol. Exp. Morphol. 8 (4), 505-526.
Hand, R., Bortone, D., Mattar, P., Nguyen, L., Heng, J. I., Guerrier, S., Boutt, E.,
Peters, E., Barnes, A. P., Parras, C., et al. (2005). Phosphorylation of Neurogenin2
specifies the migration properties and the dendritic morphology of pyramidal neurons in
the neocortex. Neuron 48, 45–62.
Handley, A., Schauer, T., Ladurner, A. G. and Margulies, C. E. (2015). Designing Cell-
Type-Specific Genome-wide Experiments. Mol. Cell 58, 621–631.
Harmansa, S. and Affolter, M. (2018). Protein binders and their applications in
developmental biology. Dev. 26, 145(2).
Hörmanseder, E., Simeone, A., Allen, G. E., Bradshaw, C. R., Figlmüller, M., Gurdon,
J. and Jullien, J. (2017). H3K4 Methylation-Dependent Memory of Somatic Cell
Identity Inhibits Reprogramming and Development of Nuclear Transfer Embryos. Cell
Stem Cell. 21(1), 135-143.
Horton, J. R., Liebert, K., Bekes, M., Jeltsch, A. and Cheng, X. (2006). Structure and
Substrate Recognition of the Escherichia coli DNA Adenine Methyltransferase. J. Mol.
Biol. 358, 559–570.
Imrichová, H., Hulselmans, G., Atak, Z. K., Potier, D. and Aerts, S. (2015). I-cisTarget
2015 update: Generalized cis-regulatory enrichment analysis in human, mouse and fly.
Nucleic Acids Res. 43, W57–W64.
Johnson, D. S., Mortazavi, A., Myers, R. M. and Wold, B. (2007). Genome-wide mapping
of in vivo protein-DNA interactions. Science. 316(5830), 1497-1502.
Karolchik, D., Hinrichs, A. S., Furey, T. S., Roskin, K. M., Sugnet, C. W., Haussler, D.
and Kent, W. J. (2004). The UCSC Table Browser data retrieval tool. Nucleic Acids
Res. 32, 493D – 496.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Kungulovski, G., Kycia, I., Tamas, R., Jurkowska, R. Z., Kudithipudi, S., Henry, C.,
Reinhardt, R., Labhart, P. and Jeltsch, A. (2014). Application of histone
modification-specific interaction domains as an alternative to antibodies. Genome Res.
24(11), 1842-1853.
Kungulovski, G., Mauser, R., Reinhardt, R. and Jeltsch, A. (2016). Application of
recombinant TAF3 PHD domain instead of anti-H3K4me3 antibody. Epigenetics and
Chromatin. 9, 11.
Lacy-Hulbert, A., Thomas, R., Li, X. P. P., Lilley, C. E. E., Coffin, R. S. S. and Roes, J.
(2001). Interruption of coding sequences by heterologous introns can enhance the
functional expression of recombinant genes. Gene Ther. 8(8), 649-653.
Landt, S. G., Marinov, G. K., Kundaje, A., Kheradpour, P., Pauli, F., Batzoglou, S.,
Bernstein, B. E., Bickel, P., Brown, J. B., Cayting, P., et al. (2012). ChIP-seq
guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 22,
1813-1831.
Marshall, O. J. and Brand, A. H. (2015). Damidseq-pipeline: An automated pipeline for
processing DamID sequencing datasets. Bioinformatics 31, 3371–3373.
Marshall, O. J. and Brand, A. H. (2017). Chromatin state changes during neural
development revealed by in vivo cell-type specific profiling. Nat. Commun. 8, 2271.
Marshall, O. J., Southall, T. D., Cheetham, S. W. and Brand, A. H. (2016). Cell-type-
specific profiling of protein-DNA interactions without cell isolation using targeted
DamID with next-generation sequencing. Nat. Protoc. 11, 1586–1598.
McGuire, S. E., Le, P. T., Osborn, A. J., Matsumoto, K. and Davis, R. L. (2003).
Spatiotemporal Rescue of Memory Dysfunction in Drosophila. Science. 302(5651),
1765-8.
Mizutani, K., Yoon, K., Dang, L., Tokunaga, A. and Gaiano, N. (2007). Differential
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Notch signalling distinguishes neural stem cells from intermediate progenitors. Nature
449, 351–355.
Ohtsuka, T., Imayoshi, I., Shimojo, H., Nishi, E., Kageyama, R. and McConnell, S. K.
(2006). Visualization of embryonic neural stem cells using Hes promoters in transgenic
mice. Mol Cell Neurosci 31, 109–122.
Potier, D., Davie, K., Hulselmans, G., NavalSanchez, M., Haagen, L., Huynh-Thu, V. A.,
Koldere, D., Celik, A., Geurts, P., Christiaens, V., et al. (2014). Mapping Gene
Regulatory Networks in Drosophila Eye Development by Large-Scale Transcriptome
Perturbations and Motif Inference. Cell Rep. 9, 2290–2303.
Promega (2009). pCI and pSI Mammalian Expression Vectors Technical Bulletin.
Rang, F. J., de Luca, K. L., de Vries, S. S., Valdes-Quezada, C., Boele, E., Nguyen, P. D.,
Guerreiro, I., Sato, Y., Kimura, H., Bakkers, J. & Kind, J. (2022). Single-cell
profiling of transcriptome and histone modifications with EpiDamID. Molecular Cell.
82(10), 1956-1970.
Reese, M. G., Eeckman, F. H., Kulp, D. and Haussler, D. (1997). Improved Splice Site
Detection in Genie. J. Comput. Biol. 4, 311–323.
Robinson, J. T., Thorvaldsdóttir, H., Winckler, W., Guttman, M., Lander, E. S., Getz,
G. and Mesirov, J. P. (2011). Integrative genomics viewer. Nat. Biotechnol. 29, 24–6.
Saito, T. and Nakatsuji, N. (2001). Efficient gene transfer into the embryonic mouse brain
using in vivo electroporation. Dev. Biol. 240, 237–246.
Sato, Y., Mukai, M., Ueda, J., Muraki, M., Stasevich, T. J., Horikoshi, N., Kujirai, T.,
Kita, H., Kimura, T., Hira, S., et al. (2013). Genetically encoded system to track
histone modification in vivo. Sci. Rep. 3.
Sato, Y., Kujirai, T., Arai, R., Asakawa, H., Ohtsuki, C., Horikoshi, N., Yamagata, K.,
Ueda, J., Nagase, T., Haraguchi, T., et al. (2016). A Genetically Encoded Probe for
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Live-Cell Imaging of H4K20 Monomethylation. J. Mol. Biol. 428, 3885–3902.
Senapathy, P., Shapiro, M. B. B. and Harris, N. L. L. (1990). Splice junctions, branch
point sites, and exons: Sequence statistics, identification, and applications to genome
project. Methods Enzymol.
Skene, P. J. and Henikoff, S. (2017). An efficient targeted nuclease strategy for high-
resolution mapping of DNA binding sites. Elife.
Southall, T. D., Gold, K. S., Egger, B., Davidson, C. M., Caygill, E. E., Marshall, O. J.
and Brand, A. H. (2013). Cell-type-specific profiling of gene expression and chromatin
binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev.
Cell 26, 101–12.
Stempor, P. and Ahringer, J. (2016). SeqPlots - Interactive software for exploratory data
analyses, pattern discovery and visualization in genomics. Wellcome Open Res. 1, 14.
Tang, J. L. Y., Hakes, A. E., Krautz, R., Suzuki, T., Contreras, E. G., Fox, P. M. and
Brand, A. H. (2022). NanoDam identifies Homeobrain (ARX) and Scarecrow
(NKX2.1) as conserved temporal factors in the Drosophila central brain and visual
system. Dev. Cell. 57(9), 1193-1207.
van den Ameele, J., Cheetham, S. W., Krautz, R., Donovan, A. P. A., Llora-Batlle, O.,
Yakob, R. and Brand, A. H. (2022). Reduced chromatin accessibility correlates with
resistance to Notch activation. Nat. Commun. 13(1), 2210.
van den Brink, S. C., Sage, F., Vértesy, Á., Spanjaard, B., Peterson-Maduro, J., Baron,
C. S., Robin, C. and van Oudenaarden, A. (2017). Single-cell sequencing reveals
dissociation-induced gene expression in tissue subpopulations. Nat. Methods 14, 935–
936.
van Steensel, B. and Henikoff, S. (2000). Identification of in vivo DNA targets of chromatin
proteins using tethered dam methyltransferase. Nat. Biotechnol. 18, 424–8.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Villaseñor, R., Pfaendler, R., Ambrosi, C., Butz, S., Giuliani, S., Bryan, E., Sheahan, T.
W., Gable, A. L., Schmolka, N., Manzo, M., et al. (2020). ChromID identifies the
protein interactome at chromatin marks. Nat. Biotechnol. 38, 728-736.
Vogel, M. J., Peric-Hupkes, D. and van Steensel, B. (2007). Detection of in vivo protein–
DNA interactions using DamID in mammalian cells. Nat. Protoc. 2, 1467–1478.
Wade, A. A., van den Ameele, J., Cheetham, S. W., Yakob, R., Brand, A. H. and Nord,
A. S. (2021). In vivo targeted DamID identifies CHD8 genomic targets in fetal mouse
brain. iScience. 24(11), 103234.
Yang, C.-P., Fu, C.-C., Sugino, K., Liu, Z., Ren, Q., Liu, L.-Y., Yao, X., Lee, L. P. and
Lee, T. (2016). Transcriptomes of lineage-specific Drosophila neuroblasts profiled by
genetic targeting and robotic sorting. Development 143, 411–421.
Zhang, Y., Liu, T., Meyer, C. A., Eeckhoute, J., Johnson, D. S., Bernstein, B. E.,
Nussbaum, C., Myers, R. M., Brown, M., Li, W., et al. (2008). Model-based analysis
of ChIP-Seq (MACS). Genome Biol. 9, R137.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Figure legends
Figure 1. Profiling histone marks with TaDa. (A) Schematic showing methods used for
Targeted DamID profiling of histone modifications. (B) Illustration demonstrating fusion of
Dam to chromatin binding domains for H3K4me3 or H3K27me3 and to mintbodies (scFvs)
against H3K9ac or H4K20me1. (C) Histone marks near the Mira locus (shaded) in
Drosophila third instar NSCs (D) Histone marks near the ANT-C locus (shaded) in
Drosophila third instar NSCs. (E) Average (±s.e.m.) signal intensity and (F) intensity of
TaDa signal across genes (transcription start/end site, TSS/TES, ±5kb) expressed in third
instar larval NSCs. (G) TaDa profiles for PolII, H3K27me3 and polycomb (Pc) around the
Antennapedia complex.
Figure 2. TaDa in mouse RGCs. A-C. Schematic of in utero TaDa in mouse RGCs. (D)
TaDa profiles at FoxG1 in mouse RGCs. (E) TaDa profiles at the HoxB locus in mouse
RGCs. (F) Schematic of TaDa during mouse cortical neurogenesis. (G) H3K4me3 TaDa
profiles at the Nestin locus throughout neurogenesis. (H-I) H3K4me3 TaDa profiles at the
NeuroG2 locus throughout neurogenesis and the mRNA expression pattern of NeuroG2 at
E14.5. (J-K) H3K4me3 TaDa profiles at the Dcx locus throughout neurogenesis and the
mRNA expression pattern of Dcx at E14.5.
Figure 3. Cell-type specific histone mark TaDa in Xenopus embryos. (A) Schematic of
H3K4me3 TaDa during Xenopus embryogenesis after injection of in vitro transcribed mRNA
in 1-cell stage embryos. (B-D) H3K4me3 TaDa profiles at the indicated developmental stages
on the Meox2 and Sox17 loci.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Supplementary figure legends
Supplementary Fig. 1. TaDa of histone modifications in Drosophila third instar larval
NSCs. (A) Pearson correlation coefficients between all TaDa replicates from Drosophila third
instar larval NSCs and Dam-only conditions. (B) Genomic feature distribution of TaDa
profiles in Drosophila third instar larval NSCs. (C) TaDa profiles in Drosophila third instar
larval NSCs at the asense locus (shaded). (D) TaDa profiles in Drosophila third instar larval
NSCs at the vnd locus (shaded). (E-F) Pearson correlation coefficients between TaDa using
chromatin readers and previously profiled chromatin binding proteins (Marshall & Brand.,
2017).
Supplementary figure 2. TaDa of histone modifications in mouse RGCs. (A-C)
Schematic overview of conventional TaDa constructs and those compatible with transient
transfection, by introducing an intron or an inversion flanked by lox-sites for recombination.
(D) Pearson correlation coefficients between all mouse RGC TaDa replicates and the Dam-
only conditions. (E) Correlations of TaDa for histone modifications in mouse RGCs to ChIP-
seq data from the ENCODE consortium database. (F) TaDa (top) and ChIP-seq (ENCODE)
profiles for histone modifications near the Nrarp genomic locus (shaded). (G) Genomic
feature distribution of histone mark TaDa profiles in mouse RGCs. (H) Average signal
intensity (±s.e.m.) of TaDa signal at genes (TSS to TES ±5kb) expressed in mouse RGCs.
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Figure 1
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Figure 2
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
Figure 3
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint
.CC-BY-NC-ND 4.0 International licensemade available under a
(which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is
The copyright holder for this preprintthis version posted April 11, 2024. ; https://doi.org/10.1101/2024.04.11.589050doi: bioRxiv preprint