library_construction_protocol: The skin tissue was homogenized in TRIzol™ using Precellys® beads.  Total RNA was extracted with the PureLink™ RNA Mini Kit, including optional DNase treatment. RNA quality and concentration were assessed using LabChip® GXII system, with an average RQS of 7.1 ± 0.9. The isolated total RNA was further purified with Ampure XP beads (Beckman Coulter, Brea, USA) and quality controlled on a Bioanalyzer (Agilent, Santa Clara, USA) RNA nano chip. 100 ng of the purified total RNA was converted into RNA-seq libraries using the CorALL kit (Lexogen, Vienna, Austria) according to the manufacturer's instructions.