Cohesin-dependence of neuronal gene expression relates to chromatin loop length

  1. Lesly Calderon
  2. Felix D Weiss
  3. Jonathan A Beagan
  4. Marta S Oliveira
  5. Radina Georgieva
  6. Yi-Fang Wang
  7. Thomas S Carroll
  8. Gopuraja Dharmalingam
  9. Wanfeng Gong
  10. Kyoko Tossell
  11. Vincenzo de Paola
  12. Chad Whilding
  13. Mark A Ungless
  14. Amanda G Fisher
  15. Jennifer E Phillips-Cremins  Is a corresponding author
  16. Matthias Merkenschlager  Is a corresponding author
  1. MRC London Institute of Medical Sciences, Imperial College London, United Kingdom
  2. Institute of Clinical Sciences, Faculty of Medicine, Imperial College, United Kingdom
  3. Department of Bioengineering, University of Pennsylvania, United States
  4. Epigenetics Program, Perelman School of Medicine, University of Pennsylvania, United States
  5. Department of Genetics, Perelman School of Medicine, University of Pennsylvania, United States
9 figures, 1 table and 7 additional files

Figures

Figure 1 with 2 supplements
Conditional cohesin deletion in post-mitotic neurons.

(a) E17.5–E18.5 cortices were dissociated and plated on poly-D-lysine. After 10 days, cultures were stained for pan neuronal (MAP2), astrocyte (GFAP), and microglia (IBA1) markers, and cell type …

Figure 1—source data 1

Figure 1: Conditional cohesin deletion in post-mitotic neurons.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig1-data1-v2.xlsx
Figure 1—figure supplement 1
Rad21 NexCre RiboTag validation.

(a) RAD21 and LAMIN B protein expression in Rad21+/+ NexCre and Rad21lox/lox NexCre cortical explant cultures was analyzed by fluorescent immunoblots. One representative blot of six is shown here, …

Figure 1—figure supplement 2
Restoration of cohesin rescues chromatin loops.

(a) Transient RAD21 depletion by Dox-inducible TEV expression and recovery after Dox washout (Weiss et al., 2021). RAD21-TEV western blots (left) and heat maps of chromatin contacts at the …

Figure 2 with 1 supplement
Loss of cohesin from immature post-mitotic neurons perturbs neuronal gene expression.

(a) Volcano plot representing log2 fold-change (FC) versus significance (-log10 of adjusted p values) of downregulated genes (1028) and upregulated genes (572) in RiboTag RNA-seq of Rad21lox/lox NexC…

Figure 2—source data 1

Figure 2: Loss of cohesin from immature post-mitotic neurons perturbs neuronal gene expression.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig2-data1-v2.xlsx
Figure 2—figure supplement 1
Gene expression in Rad21lox/lox NexCre neurons.

(a) Examples of deregulated genes in Rad21lox/lox NexCre neurons. Genes associated with autism spectrum disorders are highlighted in red. (c) Overlap between human genes associated with autism …

Figure 3 with 1 supplement
Cohesin contributes to the maturation of post-mitotic neurons.

(a) Top: Schema of cortical layers (Greig et al., 2013) showing subplate (SP), layer 6 (VI), layer 5 (V), the cortical plate (CP), and the marginal zone (MZ). Middle: Immunofluorescence analysis of …

Figure 3—source data 1

Figure 3: Cohesin contributes to the maturation of post-mitotic neurons.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig3-data1-v2.xlsx
Figure 3—figure supplement 1
Impact of cohesin loss in immature post-mitotic neurons in vivo.

(a) Expected Mendelian ratios and observed percentages of live Rad21+/+ NexCre, Rad21lox/+ NexCre, Rad21lox/lox NexCre mice at the indicated developmental stages, n=217. (b) Immunofluorescence …

Acute cohesin depletion deregulates ARG expression.

(a) Western blot documenting acute RAD21 depletion by 4-OHT-inducible RAD-TEV cleavage (left). GSEA of the gene set downregulated (DEseq2, adj. p<0.05) in RAD21-TEV neurons in Rad21lox/lox NexCre

Figure 5 with 1 supplement
Activity-regulated neuronal gene (ARG) classes differ in their reliance on cohesin.

(a) Pie charts show the expression of constitutive genes (top), immediate early genes (IEGs) (center), and late response genes (LRGs) (bottom) in Rad21 NexCre neurons under four different …

Figure 5—source data 1

Figure 5: Activity-regulated neuronal gene (ARG) classes differ in their reliance on cohesin.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig5-data1-v2.xlsx
Figure 5—figure supplement 1
Gene expression and genotype interaction analysis.

(a) The expression of constitutive and previously defined activity-regulated genes (Kim et al., 2010) was determined by RNA-seq in explant cultures of Rad21lox/lox NexCre neurons at baseline, in the …

Figure 6 with 1 supplement
The genomic distance traversed by chromatin contacts formed by neuronal genes predicts whether or not cohesin is required for their full expression.

(a) The span of Hi-C loops (left), Hi-C loops with CTCF bound to at least one of the loop anchors (middle) and Hi-C loops between promoters and inducible enhancers (right) for Immediate early genes …

Figure 6—source data 1

Figure 6: The genomic distance traversed by chromatin contacts formed by neuronal genes predicts whether or not cohesin is required for their full expression.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig6-data1-v2.xlsx
Figure 6—figure supplement 1
Promoter binding of CTCF or cohesin does not distinguish cohesin-dependent from cohesin-independent genes.

(a) Presence of CTCF (Bonev et al., 2017), RAD21 (Fujita et al., 2017), or cohesin-non-CTCF binding at gene promoters (TSS ± 2 kb) at immediate early genes (IEGs), late response genes (LRGs) that …

Figure 7 with 3 supplements
Fos enhancer-promoter contacts are robustly induced in cohesin-deficient neurons.

(a) Expression of the immediate early genes (IEG) Fos at baseline, after TTX/D-AP5 (TTX), and KCl-stimulation (left, mean log2-transformed counts from 3 biological replicates, * adj. p<0.05). (b) …

Figure 7—source data 1

Figure 7: Fos enhancer-promoter contacts are robustly induced in cohesin-deficient neurons.

https://cdn.elifesciences.org/articles/76539/elife-76539-fig7-data1-v2.xlsx
Figure 7—figure supplement 1
Inducible gene expression in cohesin-deficient neurons.

(a) Examples of Activity-regulated neuronal gene (ARG) expression at baseline, after TTX/D-AP5 (TTX), and KCl-stimulation. Mean log2-transformed counts from three biological replicates (* adj. p<0.05…

Figure 7—figure supplement 2
Contacts between the Arc promoter and an inducible enhancer in wild -type and cohesin-deficient neurons.

(a) Expression of Arc mRNA at baseline, after TTX/D-AP5 (TTX), and KCl-stimulation (top, mean log2-transformed counts from three biological RNA-seq replicates, * adj. p<0.05) and at the indicated …

Figure 7—figure supplement 3
Replicate 5 C experiments.

(a) Top: Interaction frequency zoom-in heatmaps of 250 kb region surrounding the Fos gene obtained by 5C analysis of chr12 86201802–87697802. Dashed lines and arrow heads mark major CTCF binding …

Figure 8 with 1 supplement
Bdnf enhancer-promoter contacts are weakened in the absence of cohesin.

(a) Total Bdnf transcripts at baseline, after TTX/D-AP5 (TTX), and KCl-stimulation (left, mean log2-transformed counts from three biological replicates, * adj. p<0.05). (b) Bdnf promoter IV …

Figure 8—figure supplement 1
Quantification of observed and distance-corrected loop strength.

(a) Observed values loop strength, as presented throughout the paper. (b) Observed values loop strength corrected for the distance-dependent background signal and shown on the same y-axis scale for …

Author response image 1
No scaling between gene length and cohesin-dependence in Rad21 NexCre neurons.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (Mus musculus)Rad21lox
Rad21tm1.1Mmk
DOI: 10.1038/nature10312MGI:5293824
Genetic reagent (Mus musculus)NexCre
Neurod6tm1(cre)Kan
DOI: 10.1002/dvg.20256MGI:2668659
Genetic reagent (Mus musculus)Rpl22(HA)lox (RiboTag)
Rpl22tm1.1Psam
DOI: 10.1073/pnas.0907143106MGI:4355967
Genetic reagent (Mus musculus)Rad21tev
Rad21tm1.1Kktk
DOI: 10.1101/gad.605910MGI:4840469
antibodyanti-RAD21 (rabbit polyclonal)AbcamCat #. ab154769WB: (dilution 1:1000)
IF: (dilution 1:500)
antibodyanti-LAMIN B (goat polyclonal)Santa Cruz BiotechnologyCat #. sc-6216WB: (dilution 1:10,000)
antibodyanti-rabbit IgG (H + L) Alexa Fluor 680 (goat polyclonal)ThermoFisher ScientificCat #. A-21109WB: (dilution 1:10,000)
antibodyanti-goat IgG (H + L) Alexa Fluor 680 (donkey polyclonal)ThermoFisher ScientificCat #. A-21084WB: (dilution 1:10,000)
antibodyanti-HA (rabbit polyclonal)SigmaCat #. H6908polysome immunoprecipitation
antibodyanti-GFAP (rabbit polyclonal)WakoCat #. Z0334IF: (dilution 1:500)
antibodyanti-MAP2 (chicken polyclonal)AbcamCat #. ab611203IF: (dilution 1:5000)
antibodyanti-GAD67 (mouse monoclonal)MilliporeCat #. MAB5406IF: (dilution 1:500)
antibodyanti-HA (mouse monoclonal)CovanceCat #. MMS-101RIF: (dilution 1:1000)
antibodyIBA1 (rabbit polyclonal)WakoCat #. 019–19741IF: (dilution 1:250)
antibodyanti-TUBB3 (Tuj1, mouse monoclonal)BiolegendCat #. 801,202IF: (dilution 1:500)
antibodyanti-gamma-H2AX (rabbit polyclonal)Bethyl LaboratoriesCat #. A300-081AIF: (dilution 1:3000)
antibodyanti-Cleaved Caspase-3 (Asp175) (rabbit polyclonal)Cell signallingCat #. 9,661IF: (dilution 1:400)
antibodyanti-TBR1 (rabbit polyclonal)AbcamCat #. ab31940IF: (dilution 1:1000)
antibodyanti-CTIP2 (25B6, rat monoclonal)AbcamCat #. ab18465IF: (dilution 1:500)
antibodyanti-CUX-1 (rabbit polyclonal)Santa Cruz BiotechnologyCat #. sc-13024IF: (dilution 1:400)
antibodyanti–Phospho-Histone H3 S10 Alexa Fluor 647 conjugate (rabbit polyclonal)Cell signallingCat #. 9,716IF: (dilution 1:50)
antibodyanti-rabbit IgG (H + L) Alexa Fluor 647 (goat polyclonal)ThermoFisher ScientificCat #. A-21244IF: (dilution 1:500)
antibodyanti-Rabbit IgG (H + L) Alexa Fluor 568 (goat polyclonal)ThermoFisher ScientificCat #. A-11011IF: (dilution 1:500)
antibodygoat anti-mouse IgG (H + L) Alexa Fluor 488 (goat polyclonal)ThermoFisher ScientificA-11001IF: (dilution 1:500)
antibodyanti-chicken IgY (H + L) Alexa Fluor 568 (goat polyclonal)Abcamab175711IF: (dilution 1:500)
Software, algorithmImageJ software(http://imagej.nih.gov/ij/)
Software, algorithmGraphPad Prism software(https://graphpad.com)
Software, algorithmFilamentTracer,
Imaris software, Bitplane AG
https://imaris.oxinst.com
Software, algorithmGSEA Desktop v3.0https://www.gsea-msigdb.org/gsea/index.jsp
Software, algorithmLeica Application Suite X (LAS X, v2.7) softwarehttps://www.leica-microsystems.com/products/microscope-software/p/leica-las-x-ls/
Software, algorithmCellProfiler v2.2https://cellprofiler.org
Software, algorithmImage Studio Software (v5.2)Li-cor Image Studio https://www.licor.com/bio/image-studio/

Additional files

Supplementary file 1

Gene expression analysis NexCre Ribotag RNA-seq Rad21 lox/lox minus wild-type.

https://cdn.elifesciences.org/articles/76539/elife-76539-supp1-v2.xlsx
Supplementary file 2

Gene expression analysis total RNA-seq NexCre Rad21 lox/lox minus wild-type baseline ('at rest').

https://cdn.elifesciences.org/articles/76539/elife-76539-supp2-v2.xlsx
Supplementary file 3

Gene ontology analysis of NexCre Ribotag RNA-seq Rad21 lox/lox minus wild-type.

https://cdn.elifesciences.org/articles/76539/elife-76539-supp3-v2.xlsx
Supplementary file 4

Gene expression analysis RAD21-TEV 24 h cleaved vs control at baseline.

https://cdn.elifesciences.org/articles/76539/elife-76539-supp4-v2.xlsx
Supplementary file 5

Gene expression analysis total RNA-seq NexCre Rad21 lox/lox minus wild-type TTX, KCl 1 h, KCl 6 h.

https://cdn.elifesciences.org/articles/76539/elife-76539-supp5-v2.xlsx
Supplementary file 6

Gene expression analysis RAD21-TEV 24 h cleaved vs control BDNF 30 min and BDNF 120 min.

https://cdn.elifesciences.org/articles/76539/elife-76539-supp6-v2.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/76539/elife-76539-transrepform1-v2.docx

Download links