10X Single-cell RNA-seq were performed on four mouse samples, two KOs (only one gene was knocked out) and two WTs. KO1 and KO2 have the same weeks of age with WT1 and WT2, respectively. The purpose of this study is to investigate the difference between KO and WT, as well as the impact of Knocking-out on the expression of related genes.
The sequencing were performed via the workflow "Multiple Samples, Multiple GEM Well, One Flowcell". I ran Cellranger count four times separately for the four samples, then I don't know how to do for the next step.
Option 1. Should I use Cellranger aggr to aggregate KO1 and KO2 to a single KO output, and then WT1 and WT2 to a single WT output?
Option 2. Or should I directly use Cellranger aggr to aggregate the four samples to a single output, containing all KOs and WTs?
Option 3. Or should I use Cellranger aggr to aggregate KO1 and WT1, as well as KO2 and WT2, because KO1 and WT1 were sequenced in a same flowcell?
After Cellranger aggr (no matter option 1, 2 or 3), how should I perform the next differential analysis? Could Seurat do that?
Thanks a million for your attention,