Category: Single Cell Sequencing
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 7: Genome-Wide Differential TF Activity with chromVAR and limma
Meta description: Learn genome-wide differential transcription factor activity analysis for scATAC-seq: donor pseudobulk, chromVAR motif deviations and limma moderated t-tests, with R code. Stop hand-picking motifs. Score all 633 JASPAR2020 transcription factor motifs in CD14 monocytes with chromVAR, test severe versus healthy at the level of donors with limma, and learn how to read the…
//
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 6: Footprinting Analysis
Go beyond “which motifs are enriched” and ask “is the transcription factor actually sitting there?” — run scATAC-seq footprinting two ways, with Signac Footprint() on fragments and TOBIAS on cell-type pseudobulk BAMs, and learn how far each answer can be trusted. Part 5 ended with a list of transcription factor families. In CD14 monocytes, the…
//
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 5: Differential Accessibility Analysis
Compare chromatin accessibility between disease and health inside a single cell type — with a single-cell test, a pseudobulk test, motif enrichment, and pathway analysis — and learn why the two tests disagree. Part 4 ended with 30,571 cells carrying biological names: CD14 monocytes, CD8 T effectors, B cells, NK cells. Naming cells is not…
//
-
How to Analyze Single-Cell ATAC-seq Data: A Complete Beginner’s Guide Part 4: Cell Type Identification
You have 30,571 cells sorted into 17 clusters. Every cluster is a real group of cells that share an open-chromatin landscape. And not one of them has a name. That is where Part 3 left off, and it is where every scATAC-seq analysis stalls if you stop there. “Cluster 4 is more accessible than cluster…
//
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 3: Integration and Clustering
Build a consensus peak set, normalize sparse chromatin with TF-IDF, correct batch effects with Harmony and reciprocal LSI, and find biologically meaningful clusters Prerequisites: This tutorial assumes you have completed Part 1: From FASTQ to Peaks and Part 2: Thorough Quality Control with Signac. You should have the four <sample>_qc.rds objects Part 2 produced, the…
//
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 2: Thorough Quality Control with Signac
Separate real nuclei from empty droplets, dead cells, doublets, and artifacts using TSS enrichment, nucleosome signal, FRiP, blacklist ratio, and AMULET. In Part 1 you took raw FASTQ files from a real COVID-19 case/control cohort, ran cellranger-atac count, and ended up with a fragments file, a set of peaks, and a peak-by-barcode matrix. Cell Ranger…
//
-
How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 1: From FASTQ to Peaks
Learn what chromatin accessibility tells you, then take a real 10x Genomics scATAC-seq dataset from raw FASTQ files all the way to a peak-by-cell matrix with Cell Ranger ATAC and MACS3. Prerequisites: Comfort on the Linux command line (navigating folders, running commands) and roughly 100 GB of free disk space. No prior epigenomics or programming…
//
-
How to Choose Your scRNA-seq QC Tools (Part2-2): SoupX vs DecontX and DoubletFinder vs scDblFinder
A decision-focused quick tip for picking the right ambient RNA correction and doublet detection tools — run side by side on real 10x PBMC data Introduction: Same Job, Different Tools In Part 2 of this series, we built a full quality control (QC) workflow that corrected ambient RNA contamination and removed doublets before any downstream…
//
-
How to Analyze Single-Cell RNA-seq Data — Complete Beginner’s Guide Part 17: Infer Signaling Pathway Activity with decoupleR and PROGENy
You already know the decoupleR engine from Part 16. Swap the network, swap one method, and you go from “which transcription factors are active” to “which signaling pathways are active” — in about half the code. In Part 16 you inferred per-cell transcription factor activity with decoupleR and CollecTRI, and the analysis landed on a…
//
-
How to Analyze Single-Cell RNA-seq Data — Complete Beginner’s Guide Part 16: Build Gene Regulatory Networks with decoupleR and CollecTRI
Infer per-cell transcription factor activity from a curated, signed regulatory network — no de novo network learning, no motif databases, just one fast linear model. In Part 15 you learned to redraw the entire scRNA-seq series with cleaner code. Now we return to biology with a new question: which transcription factors are driving each cell…
//
Search
Categories
- ATAC-seq (7)
- bulk RNA-seq (27)
- chromatin accessibility (14)
- Database (4)
- Epigenetics (14)
- Genomics (12)
- HPC (6)
- Insights (1)
- Metagenomics (1)
- RNA-seq (20)
- Scientific Programming (6)
- Single Cell Sequencing (27)
- Transcriptomics (28)
Recent Posts
- Bioinformatics Isn’t Dead. The Pipeline-Runner Is.
- How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 7: Genome-Wide Differential TF Activity with chromVAR and limma
- How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 6: Footprinting Analysis
- How to Analyze Single-Cell ATAC-seq Data — A Complete Beginner’s Guide Part 5: Differential Accessibility Analysis
Tags
Alternative Splicing Analysis ATAC-seq BAM ChIP-seq chromatin accessibility DESeq2 Differential Expression edgeR FASTQ GATK Mutect2 gene expression GSE282769 heatmap HOMER HPC Isoform limma MACS2 miRNA miRNA-seq MSigDB Normalization peak calling Pixi bioinformatics environment RNA-seq single-cell epigenomics beginners guide SLURM somatic mutations Transcript whole genome sequencing