Low-level repressive histone marks fine-tune gene transcription in neural stem cells

  1. Arjun Rajan
  2. Lucas Anhezini
  3. Noemi Rives-Quinto
  4. Jay Y Chhabra
  5. Megan C Neville
  6. Elizabeth D Larson
  7. Stephen F Goodwin
  8. Melissa M Harrison
  9. Cheng-Yu Lee  Is a corresponding author
  1. Life Sciences Institute, University of Michigan-Ann Arbor, United States
  2. Centre for Neural Circuits and Behaviour, University of Oxford, United Kingdom
  3. Department of Biomolecular Chemistry, University of Wisconsin-Madison, United States
  4. Department of Cell and Developmental Biology, University of Michigan Medical School, United States
  5. Division of Genetic Medicine, Department of Internal Medicine, University of Michigan Medical School, United States
  6. Rogel Cancer Center, University of Michigan Medical School, United States
7 figures, 1 table and 3 additional files

Figures

Figure 1 with 1 supplement
A single-cell gene expression atlas of type II neuroblast lineages.

(A) Summary illustration of gene and Gal4 driver expression patterns in the type II neuroblast lineage. The type II NB Gal4 driver: Wor-Gal4,Ase-Gal80. imm INP driver: R9D11-Gal4. (B) Harmonization of the scRNA-seq dataset from the entire type II neuroblast (NB) lineage generated in this study (blue) and our previously published scRNA-seq dataset which were limited to INPs and their progeny (orange). The genotype of larval brains used for scRNA-seq in this study: UAS-dcr2; Wor-Gal4, Ase-Gal80; UAS-RFP::stinger. (C) UMAPs of known cell-type-specific marker genes. Color intensity indicated scaled (log1p) gene expression value. (D) Pseudotime analysis starting from cells enriched for dpn, pnt, and RFP transcripts. (E) Left: Leiden clustering of the scRNA-seq atlas. Right: Representative UMAPs of dynamically expressed transcription factors from clusters 14 (NBs and immature INPs) and 1 (INPs). Color intensity indicated scaled (log1p) gene expression value. (F) Annotated gene expression atlas of a wild-type type II neuroblast lineage.

Figure 1—figure supplement 1
Quality control data for the scRNA-seq atlas.

(A) The expression pattern of the type II NB Gal4 driver (Wor-Gal4,Ase-Gal80) used in fluorescently labeling and sorting cell types in the type II neuroblast lineage. RFP is detectable in the nuclei of all type II neuroblasts and their progeny. (B) Violin plots showing (Top) number of genes or (Bottom) mitochondrial UMI percentage for Leiden clusters shown in 1 F. Clusters with low number of genes and higher mitochondrial UMI percentage suggest low-quality or dying cells.

Figure 2 with 1 supplement
Fruc functions through transcriptional repression to regulate stemness gene expression.

(A) Top: fru and dpn mRNAs are highly enriched in neuroblasts in cluster 14 of the scRNA-seq dataset, but only dpn mRNAs are detected in INPs in cluster 1. Bottom: Domains in Fru protein isoforms. ZF: zinc-finger DNA-binding domain. (B) Fru protein is detected in the neuroblast but not in INPs in a GFP-marked type II neuroblast lineage clone. The genotype used in this experiment is Elav-Gal4,UAS-mCD8::GFP,hs-flp; FRT82B,Tub-Gal80/FRT82B. (C) Endogenously expressed FruC tagged by a Myc epitope (fruC::Myc) is detected in type I and II neuroblasts but not in their differentiating progeny. (D–E) Immature INPs in GFP-marked wild-type type II neuroblast clones never show detectable Dpn expression (100%). 1–2 Ase- immature INP per fru-null (frusat15) type II neuroblast clone show detectable Dpn expression (100%). The genotype in this experiment is Elav-Gal4,UAS-mCD8::GFP,hs-flp; FRT82B,Tub-Gal80/FRT82B,fruSat15. (F–I) The newborn immature INP marked by cortical Miranda staining show undetectable Dpn and E(spl)mγ::GFP expression in the majority of wild-type clones (85.7%). Reducing fruC function leads to ectopic Dpn and E(spl)mγ::GFP expression in the newborn immature INP in most type II neuroblast lineages (94.4%). NB-Gal4: Wor-Gal4. (J–L) Type II neuroblasts overexpressing Fruc display characteristics of differentiation including a reduced cell diameter and aberrant Ase expression. Quantification of the cell diameter is shown in L. (type II NB>: Wor-Gal4, Ase-Gal80). (M–P) Overexpressing full-length Fruc or a constitutive transcriptional repressor form of FruC (Fruc,zf::ERD) is sufficient to partially restore differentiation in brat-null type II neuroblast clones (brat11/Df(2 L)Exel8040,hs-flp; Act5C-Gal4>FRT > FRT>UAS-GFP/UAS-fruC or ERD::fruC,zf). The percentage of GMCs per clone is shown in P. Yellow dashed line encircles a type II neuroblast lineage. White dotted line separates optic lobe from brain. white arrow: type II neuroblast; white arrowhead: Ase- immature INP; yellow arrow: Ase+ immature INP; yellow arrowhead: INP; magenta arrow: type I neuroblast; magenta arrowhead: GMC. Scale bars: 10 μm. p-values: ***<0.0005, and ****<0.00005.

Figure 2—figure supplement 1
FruA and FruB are ubiquitously expressed in neuroblast lineages.

(A) Endogenously expressed FruA::Myc or (B) FruB::Myc are ubiquitously detected in type I & II neuroblasts or their differentiating progeny. White dotted line separates the optic lobe from the brain. White dashed line separates brain from the optic lobe. white arrow: type II neuroblast; white arrowhead: Ase- immature INP; yellow arrow: Ase+ immature INP; yellow arrowhead: INP; magenta arrow: type I neuroblast; magenta arrowhead: GMC.

Figure 3 with 1 supplement
Fruc preferentially binds regulatory elements of genes uniquely expressed in neuroblasts.

(A) Genomic binding distribution of FruC-bound peaks (total # of peaks shown in parentheses) from CUT&RUN or random (set of Fruc peaks shuffled to randomly determined places in the genome) in type II neuroblast-enriched chromatin from brat-null brains (brat11/Df(2L)Exel8040). Fruc preferentially binds promoters. (B) Heatmap is centered on promoters or regulatory regions showing accessible chromatin as defined by ATAC-seq with 500 bp flanking regions and ordered by signal intensity of Fruc::Myc binding. (C) Representative z score-normalized genome browser tracks showing regions with accessible chromatin (ATAC-seq) and bound by Fruc::Myc (detected by the Myc antibody or FruCOM antibody), Su(H), Trl, or IgG at Notch, dpn, and E(spl)mγ loci. (D) Genes in cluster 14 from the scRNA-seq dataset were separated into neuroblast-enriched genes (right), immature INP-enriched genes (left), and invariant genes. The middle circle is the set of genes bound by Fruc or Su(H) (shown in parentheses). UMAPs show gene enrichment score for Fruc-bound genes uniquely expressed in neuroblasts (right), uniquely expressed in immature INPs (left), and ubiquitously expressed (invariant genes) (middle). (E) Percentage of genes defined in D bound by either Fruc, Su(H), random Fruc peaks, or random Su(H) peaks (set of Su(H) peaks shuffled to randomly determined places in the genome). (F–G) Genomic binding distribution of identified Su(H)-bound peaks (total # of peaks shown in parentheses) from CUT&RUN in type II neuroblast-enriched chromatin. Heatmap is centered on promoters or regulatory regions and ordered by signal intensity of Su(H) binding. (H–I) Genomic binding distribution of identified Trl-bound peaks (total # of peaks shown in parentheses) from CUT&RUN in type II neuroblast-enriched chromatin. Heatmap is centered on promoters or regulatory regions and ordered by signal intensity of Trl binding.

Figure 3—figure supplement 1
Quality control for determining FruC-, Su(H)- and Trl-binding using CUT&RUN.

(A) FruC::Myc is detected in all type II neuroblasts (marked by Dpn expression) in brat-null (brat11/Df(2L)Exel8040) brains homozygous for the fruC::Myc knock-in allele. (B) Genome-wide occupancy of FruC::Myc in 10 kb regions in type II neuroblasts determined by the Myc antibody and Frucom antibody are highly correlated. (C) Representative z score-normalized genome browser tracks showing chromatin accessibility (ATAC-seq) and regions bound by Fruc::Myc, FruCOM, Su(H), Trl, or IgG at tll and klu loci. (E) Top motifs identified by iCisTarget from 200 bp regions centered on fruc peak summits. Trl was identified as a top motif using all peaks or only regulator peaks, with NES scores shown in parentheses. M1BP was identified as a top motif using promoter peaks, with NES scores shown in parentheses. The percent of all, promoter, regulatory, or random peaks containing at least one motif was calculated by finding the corresponding motif distribution in the genome with HOMER. (D) Su(H) or (F) Trl is detected in all type II neuroblasts (marked by Dpn expression) in brat-null brains homozygous for the fruC::Myc knock-in allele. Scale bars: 10 μm.

Figure 4 with 1 supplement
Reduced fru function increases Notch pathway gene transcription.

(A) Loss of fruC function (fruΔC/Aj96u3) enhances the supernumerary neuroblast phenotype in brathypo (bratDG19310/11) brains heterozygous for a fru deletion (fruAj96u3/+), but loss of fruA (fruΔA/Aj96u3) or fruB (fruΔB/Aj96u3) function does not. (B–D) Overexpressing two copies of dominant-negative mam transgenes in immature INPs partially suppresses the supernumerary neuroblast phenotype in numbhypo (numbex112/15) brains. The genotype used in this experiment is numbEx112/R9D11-Gal4,numb15 or numbEx112/R9D11-Gal4,numb15; UAS-mamDN/UAS-mamDN. (E) Loss of fruC function (fruΔC/Aj96u3) enhances the supernumerary neuroblast phenotype in numbhypo brains heterozygous for a fru deletion (fruAj96u3/+), but loss of fruA (fruΔA/Aj96u3) or fruB (fruΔB/Aj96u3) function does not. (F–N) sm-FISH using intron probes confirms increased Notch and Notch target gene (dpn and klu) transcription in fruC-mutant neuroblasts comparing with control neuroblasts. CycE nascent transcripts serve as a control because CycE is not a Notch target gene, and CycE transcription is unaffected by fruC knockdown. The genotype used in this experiments is Wor-Gal4 or Wor-Gal4/UAS-fruCRNAi (O) Quantification of Notch, dpn, klu and CycE nascent transcript foci in control versus fruC-mutant neuroblasts. sm-FISH signals were counted in 8 dorsal-most neuroblasts (>6 μm in diameter) per brain lobe. CycE: 1±0.97 (n=128 neuroblasts) in control; 1.16±0.91 (n=96 neuroblasts) in fruCRNAi. dpn: 2.33±1.65 (n=88 neuroblasts) in control; 6.52±3.55 (n=104 neuroblasts) in frucRNAi. klu: 1.12±1.01 (n=80 neuroblasts) in control; 1.79±1.43 (n=112 neuroblasts) in frucRNAi. Notch: 1.06±0.88 (n=112 neuroblasts) in control; 1.77±1.07 (n=104 neuroblasts) in frucRNAi. (P) Knocking down E(spl)mγ function by RNAi in immature INPs suppresses increased supernumerary formation in numbhypo brains heterozygous for fru while having no effect on numbhypo brains alone. The genotypes used in this experiment are numbEx112/R9D11-Gal4,numb15; UAS-E(spl)mγRNAi/+ or numbEx112/R9D11-Gal4,numb15; UAS-E(spl)mγRNAi/fruAj96u3/+. White dashed line separates brain from the optic lobe. Yellow dashed line encircles a neuroblast. white arrow: type II neuroblast. Scale bars: 10 μm. p-value: NS: non-significant, *<0.05, **<0.005, ***<0.0005, and ****<0.00005.

Figure 4—figure supplement 1
Loss of fruC function enhances supernumerary neuroblast formation in brathypo and numbhypo brains.

(A–E) Representative images of brathypo (brat bratDG19310/11) brains alone, heterozygous for a fru deletion (fruAj96u3/+) or homozygous mutant for a specific fru isoform (fruΔA/Aj96u3, fruΔB/Aj96u3 or fruΔC/Aj96u3). (F) The numbEx allele was generated by imprecisely excising the numbNP2301 transposable element inserted in the 5’-regulatory region of numb. Numb protein remains detectable in numbEx homozygous neuroblasts, but is undetectable in numb-null neuroblasts. (G) Quantification of total type II neuroblasts per brain lobe in numbhypo larvae alone or heterozygous for a fru deletion that overexpress two copies of UAS-mamDN transgenes in immature INPs. The genotype used in this experiment is numbEx112/R9D11-Gal4,numb15, numbEx112/R9D11-Gal4,numb15; UAS-mamDN/UAS-mamDN, numbEx112/R9D11-Gal4,numb15; fruAj96u3/+ or numbEx112/R9D11-Gal4,numb15; fruAj96u3/UAS-mamDN. (H–L) Representative images of numbhypo (numbex112/15) brains alone, heterozygous for a fru deletion (fruAj96u3/+) or homozygous mutant for a specific fru isoform (fruΔA/Aj96u3, fruΔB/Aj96u3 or fruΔC/Aj96u3). (M) The CycE, dpn, klu, and N isoform used for generating intron probes used in sm-FISH experiment shown in Figure 4F–O. The intron used for probe generation is highlighted. Yellow dashed line separates brain from the optic lobe. p-value: ***<0.0005, and ****<0.00005.

Figure 5 with 1 supplement
Low levels of H3K27me3 are enriched in FruC-bound regions.

(A) Representative z score-normalized genome browser tracks showing FruC-binding and the enrichment of H3K4me3, H3K27ac, and H3K27me3 at the Notch locus in type II neuroblasts in control (brat11/Df(2L)Exel8040) or fruC-null (brat11/Df(2L)Exel8040; fruΔC/Aj96u3) brains. Left: Zoomed-out images showing nearest heterochromatin domains. Right: Zoomed-in images showing enrichment of histone marks in FruC-bound regions. (B) Representative z score-normalized genome browser track showing FruC-binding and the H3K27me3 throughout the chromosome arm 2 R in type II neuroblasts in control or fruC-null brains. Fru peaks are shown along with Pc domain regions called on data using similar strategy as previously used to call canonical Pc domains (Brown et al., 2018). (C) Volcano plot showing fold-change of H3K27me3 signal in overall genomic regions in fruC-null brains versus control brains. (D) Volcano plot showing fold-change of H3K27me3 signal in regions bound by Fruc in fruC-null brains versus control brains. (E) Heatmaps are centered on Fruc summits with 2 kb flanking regions in genes uniquely transcribed in type II neuroblasts in control or fruC-null brains and ordered by signal intensity of H3K27me3 enrichment calculated from TMM-normalized tracks. (F) Left: Density plots showing proportion of all 500 bp regions in the genome not bound by Fruc covered by different amounts of H3K27me3 reads in fruC-null brains vs control brains. Right: Density plots showing proportion of all 500 bp Fruc bound regions covered by different amounts of H3K27me3 reads in fruC-null brains vs control brains. (G) Dotplot representing coverage of each Fruc peak not in Pc domains vs coverage of each Pc domain. The horizontal line in the volcano plot represents -log10(0.05)=1.301. All genes above this line have a FDR <0.05.

Figure 5—figure supplement 1
Levels of active histone marks are unchanged FruC-bound regions.

(A–B) Loss of fruC function does not further exacerbate the supernumerary neuroblast phenotype in brat-null brains. Scale bars: 10 μm. (C) Volcano plot showing fold-change of H3K27ac signal in overall genomic regions in fruC-null brains versus control brains. (D) Volcano plot showing fold-change of H3K27ac signal in regions bound by Fruc in fruC-null brains versus control brains. (E) Representative z score-normalized genome browser tracks showing FruC-binding and the enrichment of H3K4me3, H3K27ac, and H3K27me3 in dpn and E(spl)mγ loci in type II neuroblasts in control (brat11/Df(2L)Exel8040) or fruC-null (brat11/Df(2L)Ex040; fruΔC/Aj96u3) brains. Left: Zoomed-out images showing nearest heterochromatin domains. Right: Zoomed-in images showing enrichment of histone marks in FruC-bound regions. (F) Volcano plot showing fold-change of H3K4me3 signal in overall genomic regions in fruC-null brains versus control brains. (G) Volcano plot showing fold-change of H3K4me3 signal in regions bound by Fruc in fruC-null brains versus control brains. (H) Heatmaps are centered on Fruc summits with 2 kb flanking regions in genes transcribed throughout the type II neuroblast lineage in control or fruC-null brains and ordered by signal intensity of H3K27me3 enrichment calculated from TMM-normalized tracks. The horizontal line in the volcano plot represents -log10(0.05)=1.301. All genes above this line have a FDR <0.05.

PRC2 subunits bind a high percentage of neuroblast-specific genes.

(A) Representative z score-normalized genome browser tracks showing regions bound by FruC, Su(z)12, and Caf-1 in type II neuroblasts in control (brat11/Df(2L)Exel8040) brains. Left: Zoomed-out images of the loci. Right: Zoomed-in images showing enrichment of PRC2 subunits in FruC-bound regions. (B–C) Heatmaps are centered on Fruc summits with 2 kb flanking regions in genes uniquely transcribed in type II neuroblasts in control brains and ordered by average signal intensity of Caf-1 and Su(z)12. Heatmap intensity is calculated from z score-normalized tracks. (D–F) The heterozygosity of E(z) (E(z)731/+) or Caf-1 (Caf-1short/+) enhances the supernumerary neuroblast phenotype in numbhypo (numbex112/15) brains. (G) Quantification of total type II neuroblasts per brain lobe in numbhypo (numbex112/15) brain alone or heterozygous for E(z) (E(z)731/+) or Caf-1 (Caf-1short/+).Scale bars: 10 μm.

FruC likely functions together with PRC2 to dampen the expression of stemness genes by promoting low levels of H3K27me3 at their cis-regulatory elements.

Loss of fruC functions leads to reduced repressive histone marks and increased stemness gene expression in neuroblasts.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (D. melanogaster)P{GawB}numbNP2301Kyoto Drosophila Stock Center104153
Genetic reagent (D. melanogaster)numbex115This studysee "Fly genetics and transgenes" in Material and Methods
Genetic reagent (D. melanogaster)fruA::Mycvon Philipsborn et al., 2014
Genetic reagent (D. melanogaster)fruB::Mycvon Philipsborn et al., 2014
Genetic reagent (D. melanogaster)fruC::Mycvon Philipsborn et al., 2014
Genetic reagent (D. melanogaster)brat11Arama et al., 2000stock #: 97265
Genetic reagent (D. melanogaster)bratDG19310Xiao et al., 2012
Genetic reagent (D. melanogaster)Df(2 L)Exel8040Bloomington Drosophila Stock Centerstock #: 7847
Genetic reagent (D. melanogaster)numb15Berdnik et al., 2002
Genetic reagent (D. melanogaster)frusat15Ito et al., 1996
Genetic reagent (D. melanogaster)fruAJ96u3Song et al., 2002
Genetic reagent (D. melanogaster)fruΔANeville et al., 2014
Genetic reagent (D. melanogaster)fruΔBNeville et al., 2014
Genetic reagent (D. melanogaster)fruΔCBilleter et al., 2006
Genetic reagent (D. melanogaster)E(z)731Anderson et al., 2011
Genetic reagent (D. melanogaster)Caf-1shortAnderson et al., 2011
Genetic reagent (D. melanogaster)UAS-fruC::MycThis studysee "Fly genetics and transgenes" in Material and Methods
Genetic reagent (D. melanogaster)UAS-fruC::ERD::MycThis studysee "Fly genetics and transgenes" in Material and Methods
Genetic reagent (D. melanogaster)UAS-dcr2; Wor-Gal4, Ase-Gal80; UAS-Stinger::RFPReichardt et al., 2018
Genetic reagent (D. melanogaster)P{GawB}elavC155, P{UAS-mCD8::GFP.L}Ptp4ELL4, P{hsFLP}1Bloomington Drosophila Stock Centerstock #: 5146
Genetic reagent (D. melanogaster)P{GAL4-Act5C(FRT.CD2).P}SBloomington Drosophila Stock Centerstock #: 4780
Genetic reagent (D. melanogaster)y[1] w[*]; P{ry[+t7.2]=neoFRT}82B P{w[+mC]=tubP-GAL80}LL3Bloomington Drosophila Stock Centerstock #: 5135
Genetic reagent (D. melanogaster)w[1118]; P{y[+t7.7] w[+mC]=GMR9D11-GAL4}attP2Bloomington Drosophila Stock Centerstock #: 40731
Genetic reagent (D. melanogaster)UAS-E(spl)mγRNAiBloomington Drosophila Stock Centerstock #: 25978
Antibodyanti-GFP
(chicken polyclonal)
Aves LabsCat# GFP-1010
RRID:AB_2307313
IF (1:2000)
Antibodyanti-Ase (rabbit polyclonal)Weng et al., 2010IF (1:400)
Antibodyanti-FruCOM (rabbit polyclonal)Ito et al., 2012IF (1:500)
CUT-&-RUN (1 μl)
Antibodyanti-Trl (rabbit polyclonal)Judd et al., 2021IF (1:500)
CUT-&-RUN (1 μl)
Antibodyanti-cMyc
(mouse polyclonal)
SigmaSKU: M4439IF (1:200)
Antibodyanti-Su(H) (mouse polyclonal)Santa CruzSKU: 398453IF (1:100),
Antibodyanti-Pros (mouse monoclonal)Lee et al., 2006aIF (1:500)
Antibodyanti-dpn (rat monoclonal, ascite)Weng et al., 2010IF (1:1000)
Antibodyanti-Mira (rat monoclonal)Lee et al., 2006aIF (1:100)
Antibodyanti-cMyc (goat polyclonal)AbcamCat#: ab9132CUT-&-RUN (1 μl)
Antibodyanti-Caf-1 (rabbit polyclonal)Tyler et al., 1996CUT-&-RUN (1 μl)
Antibodyanti-Su(z)12 (rabbit polyclonal)Loubière et al., 2016CUT-&-RUN (1 μl)
Antibodyanti-IgG (rabbit polyclonal)EpiCypherCUT-&-RUN (1 μl)
Antibodyanti-H3K9me3 (rabbit polyclonal)AbcamCat#: ab8898CUT-&-RUN (1 μl)
Antibodyanti-H3K4me3 (rabbit polyclonal)Active MotifCat#: 39159CUT-&-RUN (1 μl)
Antibodyanti-H3K27me3 (rabbit polyclonal)Sigma AldrichCat#: 07–449CUT-&-RUN (1 μl)
Antibodyanti-H3K27ac (rabbit polyclonal)Active MotifCat#: 39136CUT-&-RUN (1 μl)
Chemical compoundRhodamine PhalloidinInvitrogenCat#: R415IF (1:100)
Chemical compoundAlexa Fluor Plus 405 PhalloidinInvitrogenCat#: A30104
Chemical compoundDRAQ5AbcamCat#: ab108410
Chemical compoundPapainMillipore SigmaCat#: P4762-25MG
Chemical compoundCollagenase type IMillipore SigmaCat#: SCR103
Chemical compoundE-64Millipore SigmaCat#: E3132-1MG
Chemical compoundFetal Bovine SerumMillipore SigmaCat#: F0926-50ML
Chemical compoundSchneider’s MediaMillipore SigmaCat#: S0146-500ML
Chemical compoundProLong Gold Antifade MountantInvitrogenCat#: P36930
Chemical compoundProLong Gold Antifade Mountant with DNA Stain DAPIInvitrogenCat#: P36935
Chemical compoundAgencourt AMPure XP - 5 mLBeckman CoulterCat#: A63880
Chemical compoundSodium butyrateSigma AldrichCat#: B5887
Commercial assay or kit10x chromium v3 single-cell gene expression kit10x GenomicsCat# 1000154
Commercial assay or kitB1-labeled HCR RNA-FISH intron probe (Notch)Molecular Instruments, Incdm6, chrX: 3159233–3163479
Commercial assay or kitB1-labeled HCR RNA-FISH intron probe (dpn)Molecular Instruments, Incdm6, chr2R: 8230502–8231599 & 8231694–8231923
Commercial assay or kitB1-labeled HCR RNA-FISH intron probe (klu)Molecular Instruments, Incdm6, chr3L: 10985072–10996227
Commercial assay or kitB1-labeled HCR RNA-FISH intron probe (CycE)Molecular Instruments, Incdm6, chr2L: 15731785–15737628
Commercial assay or kitHCR RNA-FISH amplifiersMolecular Instruments, Inc
Commercial assay or kitCUTANA ChIC/CUT&RUN KitEpicypherSKU: 14–1048
Commercial assay or kitNEBNext Ultra II DNA Library Prep Kit for IlluminaNew England BiolabsCat#: E7645L
Commercial assay or kitNEBNext Multiplex Oligos for Illumina (96 Unique Dual Index Primer Pairs)New England BiolabsCat#: E6440L
Software, algorithmScript FilesThis StudyGITHUB
Software, algorithmCell Ranger (6.0.1)10x GenomicsRRID:SCR_017344
Software, algorithmFiji/ImageJSchindelin et al., 2012RRID:SCR_002285 https://fiji.sc/
Software, algorithmscanpy scRNA-seq analysis softwareWolf et al., 2018RRID:SCR_018139
Software, algorithmHarmonyKorsunsky et al., 2019
Software, algorithmBioinfokitRenesh Bedre, 2020, March 5
Software, algorithmMatplotlibhttps://zenodo.org/badge/
https://doi.org/10.5281/zenodo.592536.svg
RRID:SCR_008624
Software, algorithmSeabornhttps://doi.org/10.21105/joss.03021RRID:SCR_018132
Software, algorithmFastQCWingett and Andrews, 2018RRID:SCR_014583
Software, algorithmcutadaptMartin, 2011RRID: SCR_011841
Software, algorithmSamtoolsLi et al., 2009RRID: SCR_002105
Software, algorithmMACS2Zhang et al., 2008RRID: SCR_013291
Software, algorithmCUT-RUNTools-2.0Yu et al., 2021
Software, algorithmbedtoolsQuinlan and Hall, 2010RRID: SCR_008848
Software, algorithmHOMERHeinz et al., 2010RRID: SCR_010881
Software, algorithmdeeptoolsRamírez et al., 2016RRID: SCR_016366
Software, algorithmIntegrative Genomics ViewerRobinson et al., 2011RRID: SCR_011793
Software, algorithmGoPeaksYashar et al., 2022
Software, algorithmSICERZang et al., 2009RRID:SCR_010843
Software, algorithmfeatureCountsLiao et al., 2014RRID: SCR_012919
Software, algorithmEdgeRRobinson et al., 2010RRID: SCR_012802
Software, algorithmdiffRepsShen et al., 2013RRID:SCR_010873
Software, algorithmiCisTargetImrichová et al., 2015
Software, algorithmXSTREMEGrant and Bailey, 2021RRID:SCR_001783
Software, algorithmLAS AFLeica MicrosystemsRRID:SCR_013673

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Arjun Rajan
  2. Lucas Anhezini
  3. Noemi Rives-Quinto
  4. Jay Y Chhabra
  5. Megan C Neville
  6. Elizabeth D Larson
  7. Stephen F Goodwin
  8. Melissa M Harrison
  9. Cheng-Yu Lee
(2023)
Low-level repressive histone marks fine-tune gene transcription in neural stem cells
eLife 12:e86127.
https://doi.org/10.7554/eLife.86127