library_construction_protocol: A small piece of skin tissue was flash frozen (for RNA isolation from the whole skin tissue) and rest of the skin tissue was enzymatically digested using collagenase-trypsin method. After preparing the single cell suspension, dead cells were stained with Propidium Iodide (PI). Live GFP+ keratinocytes were FACS-purified and directly collected in TRIzol LS Total RNAs were isolated following the manufacturer's protocol with the following additions: after the first phase separation, additional chloroform extraction step of the aqueous layer in Phase-lock Gel heavy tubes (Quanta Biosciences); addition of 1ul Glycoblue (ThermoFisher) immediately prior to isopropanol precipitation; two washes of the RNA pellet with 75% ethanol. Quality of each sample was confirmed by measuring RNA concentration using Qubit RNA HS kit (ThermoFisher) and by determining RNA integrity with Fragment Analyzer (Agilent). Ribosomal RNA was subtracted by hybridization from total RNA samples using the NEBNext rRNA Depletion Kit (New England Biolabs). TruSeq-barcoded RNAseq libraries were generated with the NEBNext Ultra II RNA Library Prep Kit (New England Biolabs) and libraries were  sequenced on a NovaSeq with 2x150bp paired-end settings.