1,459
Views
0
CrossRef citations to date
0
Altmetric
Research Article

Integrative CUT&Tag-RNA-Seq Analysis of Histone Variant MacroH2A1-Dependent Orchestration of Human Induced Pluripotent Stem Cell Reprogramming

, , ORCID Icon, ORCID Icon, , , , ORCID Icon & ORCID Icon show all
Pages 863-877 | Received 24 Jul 2023, Accepted 29 Sep 2023, Published online: 17 Oct 2023

Figures & data

Figure 1. CUT&Tag and RNA-Seq integrated data analysis workflow.
Figure 1. CUT&Tag and RNA-Seq integrated data analysis workflow.
Figure 2. MacroH2A1.1 and macroH2A1.2 occupancy profiles of differentially expressed genes and their regulatory regions.

(A) Volcano plots of the significant differentially expressed genes upon histone binding for macroH2A1.1 (left) and macroH2A1.2 (right). (B) Distribution of peaks in differentially expressed genes relative to the (left) gene structure and (right) distance from the nearest TSS.

Figure 2. MacroH2A1.1 and macroH2A1.2 occupancy profiles of differentially expressed genes and their regulatory regions. (A) Volcano plots of the significant differentially expressed genes upon histone binding for macroH2A1.1 (left) and macroH2A1.2 (right). (B) Distribution of peaks in differentially expressed genes relative to the (left) gene structure and (right) distance from the nearest TSS.
Figure 3. Stacked bar plot of normalized peak frequency of macroH2A1.1 and macroH2A1.2 for each human chromosome.

Bars are split into gene regions.

Figure 3. Stacked bar plot of normalized peak frequency of macroH2A1.1 and macroH2A1.2 for each human chromosome.Bars are split into gene regions.
Figure 4. MacroH2A1.1 and macroH2A1.2 genome binding to differentially expressed genes involved in functions related to the three embryonic germ layers.

(A) Circos plot representing differentially expressed genes upon histone binding (UHB-DEGs) for macroH2A1.1 and macroH2A1.2, and several features related to their expression, gene ontology (GO) terms and CUT<&Tag peaks. From the outer to the inner tracks of the Circos plots: chromosomes, CUT<&Tag peaks, heatmap of fold-change expression values (red and blue mean upregulated and downregulated genes in macroH2A1.1 or macroH2A1.2 overexpressing cells), UHB-DEGs whose biological processes relate to mesoderm (light blue links), endoderm (light purple) and ectoderm (light green). Gray links were genes associated with more than one germ layer. (B) 100% stacked-bars of the frequencies of GO terms grouped by germ layers. (C) Venn diagram of UHB-DEGs associated with the three germ layers.

Figure 4. MacroH2A1.1 and macroH2A1.2 genome binding to differentially expressed genes involved in functions related to the three embryonic germ layers. (A) Circos plot representing differentially expressed genes upon histone binding (UHB-DEGs) for macroH2A1.1 and macroH2A1.2, and several features related to their expression, gene ontology (GO) terms and CUT<&Tag peaks. From the outer to the inner tracks of the Circos plots: chromosomes, CUT<&Tag peaks, heatmap of fold-change expression values (red and blue mean upregulated and downregulated genes in macroH2A1.1 or macroH2A1.2 overexpressing cells), UHB-DEGs whose biological processes relate to mesoderm (light blue links), endoderm (light purple) and ectoderm (light green). Gray links were genes associated with more than one germ layer. (B) 100% stacked-bars of the frequencies of GO terms grouped by germ layers. (C) Venn diagram of UHB-DEGs associated with the three germ layers.
Supplemental material

Supplemental Information

Download Zip (4 MB)

Data sharing statement

The RNA-seq data [Citation22] was deposited into GEO with the dataset identifier GSE164396. Genomics sequences obtained in this study by CUT&Tag has been deposited into GEO with the dataset identifier GSE214013.