A workflow for characterizing stage-specific transcription factor expression patterns during control/endogenous differentiation using bulk time-series RNA-seq data.
The full analysis is available as a GitHub Pages site: https://zohebkhan1.github.io/temporal-transcription-factor-identification/
This workflow uses publicly available data from:
Martinez JL et al. Transcriptional consequences of trisomy 21 on neural induction. Front Cell Neurosci. 2024;18:1341141. doi:10.3389/fncel.2024.1341141
Stage-specific transcription factor expression modules identified via hierarchical clustering, sorted by peak expression timing.
Mean expression trends with confidence intervals for each TF cluster across differentiation timepoints.
Overlay of all cluster expression trajectories showing the temporal wave pattern.
Comparison of TF cluster expression between control (D21) and experimental (T21) samples.
The following R packages are required:
- DESeq2
- AnnotationDbi
- org.Hs.eg.db
- edgeR
- dplyr
- tidyr
- tibble
- ggplot2
- pheatmap
- RColorBrewer
- viridis
- patchwork
Install Bioconductor packages:
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install(c("DESeq2", "AnnotationDbi", "org.Hs.eg.db", "edgeR"))Install CRAN packages:
install.packages(c("dplyr", "tidyr", "tibble", "ggplot2", "pheatmap", "RColorBrewer", "viridis", "patchwork"))


