from three independent donors in 24-well plates were treated with 1000 IU/mL
IFN-α, IFN-β, IFN-γ or medium alone for 10 h. RNA was
isolated by using the ReliaPrep RNA Cell Miniprep System (Promega) according to
manufacturers’ instructions. Microarrays were performed by ServiceXS Ltd.
(Leiden, Netherlands). 200 ng of total RNA were used to synthesize biotinylated
cRNA target, which was subject to QC checks, then hybridized (750 ng cRNA per
sample) to Illumina HT-12 v4 Expression Bead Chips (Illumina Inc. San Diego,
USA). Data from two Illumina HT-12 v4 Expression Bead Chips per sample were
background corrected in Illumina Beadstudio with further analysis carried out
using the Lumi and Limma Bioconductor packages in R (33 (link), 34 (link)).
Normalization of background corrected data was applied through variance
stabilizing transformation (VST) and robust spline normalization (RSN) in Lumi
(35 (link)). Testing for differential
expression between IFN-treated patient and control cells was done with linear
models and empirical Bayesian statistics, utilizing a multi-level experiment
approach through Limma. Gene lists for each comparison were generated. Data
plots were created using the ggplot2 library in R.