SRX23765823 GSM8115270_r1 SRP492046 PRJNA1081392 SRS20591953 GSM8115270 GSM8115270 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765824 GSM8115271_r1 SRP492046 PRJNA1081392 SRS20593299 GSM8115271 GSM8115271 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765825 GSM8115272_r1 SRP492046 PRJNA1081392 SRS20591954 GSM8115272 GSM8115272 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765826 GSM8115273_r1 SRP492046 PRJNA1081392 SRS20591956 GSM8115273 GSM8115273 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765827 GSM8115274_r1 SRP492046 PRJNA1081392 SRS20591955 GSM8115274 GSM8115274 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765828 GSM8115275_r1 SRP492046 PRJNA1081392 SRS20592753 GSM8115275 GSM8115275 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765829 GSM8115276_r1 SRP492046 PRJNA1081392 SRS20591957 GSM8115276 GSM8115276 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765830 GSM8115277_r1 SRP492046 PRJNA1081392 SRS20591958 GSM8115277 GSM8115277 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765831 GSM8115278_r1 SRP492046 PRJNA1081392 SRS20591959 GSM8115278 GSM8115278 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765832 GSM8115279_r1 SRP492046 PRJNA1081392 SRS20591960 GSM8115279 GSM8115279 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765833 GSM8115280_r1 SRP492046 PRJNA1081392 SRS20593681 GSM8115280 GSM8115280 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765834 GSM8115281_r1 SRP492046 PRJNA1081392 SRS20591961 GSM8115281 GSM8115281 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765835 GSM8115282_r1 SRP492046 PRJNA1081392 SRS20591962 GSM8115282 GSM8115282 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765836 GSM8115283_r1 SRP492046 PRJNA1081392 SRS20591963 GSM8115283 GSM8115283 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765837 GSM8115284_r1 SRP492046 PRJNA1081392 SRS20591964 GSM8115284 GSM8115284 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000 SRX23765838 GSM8115285_r1 SRP492046 PRJNA1081392 SRS20591965 GSM8115285 GSM8115285 RNA-Seq TRANSCRIPTOMIC cDNA Sampling was carried out 6 and 24 hpt and for each replicate 20 leaves were pooled (4 leaves/plant) and they were immediately frozen in liquid nitrogen. The material was finely grounded and homogenized with a Tissue Lyser II. RNA was extracted from 100 mg of the grounded leaf material coming from each replicate using the PureLink™ Plant RNA Reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer's instructions. The remaining DNA was digested with the TURBO DNA-freeTM Kit (Invitrogen Life Technologies, Carlsbad, CA, USA). RNA “TruSeq Stranded mRNA Sample Prep kit” (Illumina, San Diego, CA) was used for library preparation following the manufacturer's instructions, starting with 1-2 µg of good quality RNA (RIN > 7) as input. The RNA was fragmented 3 min at 94°C and every purification step was performed using 0.81X Agencourt AMPure XP beads. Both RNA samples and final libraries were quantified using the Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA) and quality tested by Agilent 2100 Bioanalyzer RNA Nano assay (Agilent technologies, Santa Clara, CA, USA). Libraries were then processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions and sequenced on paired-end (2x150 bp, 30M reads per sample) at the multiplexing level requested on NovaSeq6000 (Illumina, San Diego, CA, USA). Illumina NovaSeq 6000