The sequencing libraries were prepared using CEL-Seq2 (Hashimshony et al., 2016 (link)). RNA from sorted cells was extracted using TRIzol reagent (ThermoFisher, 15596018). 10 ng RNA was used for first strand cDNA synthesis using barcoded primers (the specific primers for each sample were listed in Supplementary Table S9). The second strand was synthesized by NEBNext Second Strand Synthesis Module (NEB, E6111L). The pooled dsDNA was purified with AMPure XP beads (Beckman Coulter, A63880), and subjected to in vitro transcription (IVT) using HiScribe T7 High Yield RNA Synthesis Kit (NEB, E2040S), then treated with ExoSAP-IT (Affymetrix, 78200). IVT RNA was fragmented using RNA fragmentation reagents (Ambion), and underwent another reverse transcription step using random hexamer RT primer-5’-GCC TTG GCA CCC GAG AAT TCC ANN NNN N-3’ to incorporate the second adapter. The final library was amplified with indexed primers: RP1 and RPI1 (Supplementary Table S9), and the bead purified library was quantified with 4200 TapeStation (Agilent Technologies), and paired-end sequenced on Nextseq 500 V2 (Illumina), Read 1: 15 cycles; index 1: 6 cycles; Read 2: 60 cycles.