library_construction_protocol: RNA was extracted using a hybrid method utilizing Trizol reagent and Qiagen RNeasy kit: The tissue samples were homogenized in the presence of 1 ml Trizol and 1.5 mm stainless steel bead, using the Qiagen TissueLyser II instrument for 2 min at 25 Hz, twice. Homogenized samples were incubated for 5 min at room temperature after which 200 ml chloroform was added, mixed well, and incubated for 3 min at room temperature. Samples were centrifuged at 10,000 x g at 4°C, after which the aqueous layer was transferred to a new tube and mixed with equal volume of 70%. The samples were loaded onto Qiagen RNAeasy column, and the samples were processed according to the protocol, including the optional on column DNAse treatment.  RNA was eluted in 35 µl of RNAse free water. A Quality Control analysis was performed on the RNA samples using the Agilent Bioanalyzer to assess quality and an RNA Qubit assay was used to determine quantity.  Poly(A) RNA enrichment was conducted using NEBNext Poly(A) mRNA Magnetic Isolation Module (NEB), and the sequencing library was constructed by using the NEBNext Ultra II RNA Library Prep Kit (NEB) (P/N: E7765L, NEB, Ipswich, MA) following the manufacturer's instructions.  End repair, A-tailing, and adapter ligation was performed to generate the final cDNA library. The library quality was assessed using a Bioanalyzer and quantified using a qPCR-based method with the KAPA Library Quantification Kit (P/N: KK4873) and the QuantStudio 12K instrument. Next, prepared libraries were pooled in equimolar ratios, and the resulting pool was subjected to cluster generation using the NovaSeq 6000 System, following the manufacturer's protocols. 150 bp paired-end sequencing was performed on the NovaSeq 6000 platform targeting 50M reads per sample.  Raw sequencing data (FASTQ files) obtained from the NovaSeq 6000 was subjected to quality control analysis, including read quality assessment. Real Time Analysis Software (RTA) and NovaSeq Control Software (NCS) (1.8.0; Illumina) were used for base calling. MultiQC (v1.7; Illumina) was used for data quality assessments.