Figures and data
Chromatin accessibility dynamic during GABAeric interneuron differentiation
A. Schematic illustration of GABAergic interneuron differentiation protocol. Cells from iPSC differentiation were collected for ATAC-Seq. Interneuron differentiation was grouped into four stages: iPSC (Day 0), neural progenitor cells; NPCs (Day 19), immature neurons, imN (Day 35), GABAergic interneurons, GABA (Day 65).
B. Heatmap of four Clusters identified by time series analysis using temporal changes in chromatin accessibility compared to D0 (Day 0) and other points. The signal strengths were denoted by color intensities.
C. Line plot of four clusters showing the average changes in chromatin accessibility.
D. Barplot for genomic distribution of differential chromatin accessible regions for each cluster.
E. Heatmap for top 10 TFs enrichment in each cluster. The significance was denoted by color intensities. O-S-E-N = OCT4-SOX2-TCF-NANOG.
F. Bubble plot of KEGG pathway enrichment for each cluster. p-value and enrichment were indicated.
G. Genome browser view showing representative differential chromatin accessible regions at the indicated gene loci.
Chromatin accessibility dynamic during GABAergic differentiation in Dravet Syndrome patients
A. Schematic Illustration of the design for GABAergic differentiation in Dravet syndrome patients. Cells from iPSC differentiation were collected for ATAC-Seq. GABAergic differentiation in Dravet Syndrome patients was grouped into four stages: iPSC (Day 0), NPCs (Day 19), imN (Day 35), GABA(Day 65).
B. Heatmap of four clusters identified by time series analysis using temporal changes in chromatin accessibility compared D0 (Day 0) and other points. The signal strengths were denoted by color intensities
C. Line plot of four clusters showing the average changes in chromatin accessibility.
D. Barplot for genomic distribution of differential chromatin accessible regions for each cluster.
E. Heatmap for top 10 TFs enrichment in each cluster. The significance was denoted by color intensities. O-S-E-N = OCT4-SOX2-TCF-NANOG.
F. Bubble plot of KEGG pathway enrichment for each cluster. p-value and enrichment were indicated.
G. Genome browser view showing representative differential chromatin accessible regions at the indicated gene loci.
Common chromatin accessibility features between control and Dravet Syndrome patients during GABAergic differentiation
A. Heatmap for common chromatin accessible regions detected between control and Dravet Syndrome (DS) group during differentiation. Four clusters were obtained from time series analysis. Representative genes are labeled on the right side of heatmap.
B. Barplot for genomic distribution of chromatin accessible regions for each cluster.
C. Bubble plot of KEGG pathway enrichment for each cluster. p-value and enrichment were indicated.
D. Genome browser view showing representative chromatin accessible regions at the indicated gene loci.
E. Transcription factors (TFs) deviation Z scores heatmap for unique TF enriched in each cluster. Representative TFs are labeled on the right side of heatmap.
F. Representative TFs enrichment (deviation Z score) dynamics at ATAC-Seq peaks in control and Dravet Syndrome (DS) group during differentiation.
Unique chromatin accessibility features between control and Dravet Syndrome patients during GABAergic differentiation
A. Principal component analysis (PCA) plot of all ATAC-Seq from control and Dravet Syndrome patients during GABAergic differentiation.
B. Heatmap showing the differential chromatin accessible peaks for comparison at each time point between the control and Dravet Syndrome (DS) patients. Representative genes are labeled on the right side of heatmap.
C. Barplot for genomic features of unique chromatin accessible regions at each time point for the control and Dravet Syndrome (DS) patients.
D. Bubble plot of KEGG pathway enrichment at GABAeric interneuron (Day 65) for the control and Dravet Syndrome (DS) patients. p-value and enrichment score were indicated.
E. Heatmap showing transcription factors (TFs) enrichment at differential chromatin accessible peaks for comparison at each time point between the control and Dravet Syndrome (DS) patients. Representative TFs are labeled on the right side of heatmap.
F. Representative TFs enrichment (deviation Z score) dynamics at ATAC-Seq peaks in control and Dravet Syndrome (DS) group during differentiation.
Chromatin accessibility response in valproic acid treatment of GABAergic interneuron
A. Schematic illustration showing measurement of chromatin accessibility response in valproic acid (VPA) treatment of GABAergic interneurons.
B. Principal component analysis (PCA) plot of ATAC-Seq with and without VPA treatment of GABAergic interneuron in control and Dravet Syndrome patients.
C. Heatmap of differential chromatin accessible peaks from the VPA responsible group, consisting of two control samples and one Dravet Syndrome (DS) patient sample.
D. Barplot for genomic distribution of VPA responsible accessible regions from differential conditions.
E. Bubble plot of KEGG pathway enrichment for each condition. p-value and enrichment were indicated.
F. TF deviation Z score heatmap for top 50 transcription factors (TFs) enriched at responsible chromatin regions from differential conditions.
G. Genome browser view showing representative VPA responsible chromatin accessible regions at the indicated gene loci.
Quality control of ATAC-Seq samples during GABAeric interneuron differentiation.
A. Representative pictures of immunofluorescent staining of Ctl-iPSC derived GABAergic neurons at D65 (Ctl-GABA). The Ctl-GABA cells form networks and stain positive for the neuronal markers MAP2 and double cortin (DCX), the GABAergic marker γ-aminobutyric acid (GABA) and the interneuronal markers somatostatin (SST), parvalbumin (PV) and glutamate decarboxylase 1 (GAD1). Size bar 100 μm.
B. Pearson correlation matrix from technical duplicates of ATAC-Seq at different time points.
C. Quantification of fraction of reads in peaks (FRiP) of ATAC-Seq in each condition. The dotted line indicates the cutoff used.
D. Transcription start site (TSS) enrichment score of ATAC-Seq in each sample.
E. Genomic distribution of chromatin accessible peaks for each sample.
F. Principal component analysis of ATAC-Seq from each sample.
Chromatin accessibility clusters identification during GABAeric interneuron differentiation.
A. Volcano plots show differential accessible chromatin peaks between D0 and different time points. D0 (n = 6046) vs. D19 (n = 3866); D0 (n = 8880) vs. D35 (n = 7865); D0 (n = 7870) vs. D65 (n = 7776). The dotted lines indicate the cutoff used to identify differential ATAC-Seq peaks.
B. Variance with different numbers of clusters derived from chromatin accessibility during GABAeric interneuron differentiation
C. The scaled average chromatin accessibility during GABAeric interneuron differentiation from each cluster.
D. The tendency of chromatin accessibility during GABAeric interneuron differentiation from each cluster.
Quality control of ATAC-Seq samples during GABAeric interneuron differentiation in Dravet Syndrome patients.
A. Representative pictures of immunofluorescent staining of DS-iPSC derived GABAergic neurons at D65 (DS-GABA). The DS-GABA cells form networks and stain positive for the neuronal markers MAP2 and double cortin (DCX), the GABAergic marker γ-aminobutyric acid (GABA) and the interneuronal markers somatostatin (SST), parvalbumin (PV) and glutamate decarboxylase 1 (GAD1). Size bar 100 μm.
B. Pearson correlation matrix from technical duplicates of ATAC-Seq from all samples.
C. Quantification of fraction of reads in peaks (FRiP) of ATAC-Seq in each condition. The dotted line indicates the cutoff was used.
D. Transcription start site (TSS) enrichment score of ATAC-Seq in each sample.
E. Genomic distribution of chromatin accessible regions for each sample.
F. Principal component analysis of ATAC-Seq from each sample.
Chromatin accessibility clusters identification during GABAeric interneuron differentiation in Dravet Syndrome patients.
A. Volcano plots show differential accessible chromatin peaks between D0 and different time points. D0 (n = 7798) vs. D19 (n = 5183); D0 (n = 10154) vs. D35 (n = 7239); D0 (n = 6909) vs. D65 (n = 5667). The dotted lines indicate the cutoff used to identify differential ATAC-Seq peaks.
B. Variance with different numbers of clusters derived from chromatin accessibility during GABAeric interneuron differentiation in Dravet Syndrome patients.
C. The scaled average chromatin accessibility during GABAeric interneuron differentiation from each cluster in Dravet Syndrome patients.
D. The tendency of chromatin accessibility during GABAeric interneuron differentiation from each cluster in Dravet Syndrome patients.
Identification of common chromatin accessibility between control and Dravet Syndrome patients during GABAeric interneuron differentiation.
A. Venn diagrams showing the shared and distinct accessible chromatin peaks between control and patient samples.
B. Variance with different numbers of clusters derived from common accessible chromatin regions between control and Dravet Syndrome patients.
C. The tendency of chromatin accessibility from common accessible chromatin regions between control and Dravet Syndrome patients during GABAeric interneuron differentiation.
Identification of unique accessible chromatin region between control and Dravet Syndrome patients during GABAeric interneuron differentiation.
A. Volcano plot of differential accessible chromatin peaks between control and Dravet Syndrome (DS) patients at different stage of differentiation (NPC (Day 19), imN (Day 35) and GABAeric interneuron (Day 65)). Dotted lines indicate cutoff from fold change and false discovery rate.
B. KEGG signal pathway enrichments from unique accessible chromatin region in control and Dravet Syndrome patients at NPC (Day 19), imN (Day 35).
Valproic acid treatment reshapes chromatin accessibility of GABAeric interneuron.
A. Pearson correlation matrix from technical duplicates of ATAC-Seq of VPA treated control and Dravet Syndrome patient group.
B. Distance matrix of ATAC-Seq from VPA treated control and Dravet Syndrome patient group.
C. Quantification of fraction of reads in peaks (FRiP) of ATAC-Seq in each condition.
D. Transcription start site (TSS) enrichment score of ATAC-Seq in each sample.
E. Heatmap of differential chromatin accessible regions from VPA non-responsible group, consisting of one control sample and one Dravet Syndrome (DS) patient sample.
F. Genomic distribution of chromatin accessible regions for each individual sample.
