design: For one sample, single cells of two mice were combined. Single cell library preparation and sequencing were conducted by Singulomics (Bronx, USA). Cryopreserved viable single cell suspensions were thawed, washed, resuspended in cell culture media with 0.04% bovine serum albumin and counted. Viable cell suspensions were then loaded into the Chromium Controller (10x Genomics) to generate gel beads-in-emulsion (GEM), with each GEM containing a single cell as well as barcoded oligonucleotides. We targeted 10,000 cells to be captured per sample. Next, the GEMs were placed in the SimpliAmp 96-well Thermal Cycler (Thermo Fisher Scientific) and reverse transcription was performed in each GEM (GEM-RT). After the reaction, the complementary cDNA was amplified and cleaned using Silane DynaBeads (10x Genomics) and the SPRI select Reagent kit (Beckman Coulter). Amplified full-length cDNAs from poly-adenylated mRNA were then used to generate a 3 Gene Expression library (Chromium Next GEM 3 Single Cell Reagent kits v3.1, dual index) following the manufacturer's instructions (10x Genomics). Amplified cDNAs and the libraries were measured using Qubit dsDNA HS assay (Thermo Fisher Scientific) and quality was assessed using BioAnalyzer (Agilent Technologies). The libraries were sequenced with ~200 million PE150 reads per sample on Illumina NovaSeq 6000, and fastq files of two samples were generated using Illumina bcl2fastq. Fastq files were subsequently processed using 10x Genomics Cell Ranger analytical pipeline (v4.0.0) and mouse mm10 reference.