library_construction_protocol: Single cell suspensions were prepared in calcium and magnesium-free 1X PBS containing 0.04% weight/volume BSA non-acetylated(400ug/ml), were loaded on a Next GEM Chip A (PN-1000009) together with Chromium Single Cell 3' Gel Bead Kit v2 (PN 120235) and single cells were captured on a 10X Genomics Chromium controller with a recovery target of 5000 cells per sample.  cDNAs were generated following the 10X Genomics protocol. cDNAs were generated following the 10X Genomics protocol. cDNAs were size selected using Ampure XP beads from Beckman Coulter (A63881), and their quality was checked with a Bioanalyzer High Sensitivity DNA Kit from Agilent (5067-4626). One quarter of the total cDNA was used to generate libraries using Chromium Single Cell 3' Library Kit v2 (PN-120234) and barcoded using the Chromium Multiplex Kit, 96 rxns (PN-120262) following the 10X protocol. Libraries were size selected using Ampure XP beads and their size distribution were assessed with the Bioanalyzer High Sensitivity DNA Kit (average size of 470 bp). Libraries were sequenced on Illumina Novaseq S Prime flowcell with a depth of 57,000-94,000 PE28-91 reads per cell. scRNAseq, paired-end sequencing on NovaSeq 6000