RNA sequencing of drug combinations for SL-176 and GSK-J4 in neuroblastoma cell lines

Total RNA was subjected to quality control with the Agilent TapeStation (Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s instructions. To construct libraries suitable for Illumina sequencing, the Illumina TruSeq Stranded mRNA sample preparation protocol which includes mRNA isolation, cDNA synthesis, ligation of adapters and amplification of indexed libraries was used. The yield and quality of the amplified libraries was analyzed using Qubit by ThermoFisher and the Agilent TapeStation.

The indexed cDNA libraries were normalized and combined, and the pools were sequenced on the Illumina NextSeq 550 (Illumina, San Diego, CA, USA) for a 75-cycle v2 sequencing run generating 75 bp single-end reads. Basecalling and demultiplexing was performed using CASAVA software (Illumina) with default settings generating Fastq files for further downstream mapping and analysis. Raw data was adapter and quality trimmed using Cutadapt, and aligned to the human genome build (Gene reference: Ensembl Homo_sapiens.GRCh38.104.gtf) using STAR.

Gene assignment was performed using FeatureCounts. Differential gene expression was determined using DEseq2.