mRNA was enriched from 1 µg total RNA using NEBNext® Poly(A) mRNA Magnetic Isolation Module (New England Biolabs E7490L). The library was prepared from the mRNA enriched fraction using NEBNext® Ultra™ II Directional RNA Library Prep (New England Biolabs E7765L) according to the manufacturer’s instructions. The library was sequenced on the Illumina HiSeq 2500 platform, generating 150 bp end reads. The reads were trimmed for adapters using Trimmomatic v0.3661 (link) and aligned to GRCh38 genome assembly using STAR aligner v2.662 (link) and gene-wise read coverage was quantified using Feature Counts from the Subread package v1.6363 (link). All downstream analysis was performed using R 3.5.164 . Genes with TPM below 1 in every sample were discarded. DESeq265 (link) was used to identify differentially expressed genes. Gene ontology pathway analysis were done using DAVID v6.866 (link),67 (link), Enrichr68 (link),69 (link), and Panther v15.070 .
Free full text: Click here