Message RNA was purified from PBMCs, labeled and hybridized to Affymetrix Exon 1.0 ST Array for 35 patients of Howard cohort (Zhang et al. 2014a (link)), and to Affymetrix Gene 2.0 ST Array for 134 SCD patients selected from UIC cohort. Probe sequences were aligned to human genome assembly GRCh37 to select for probes with unique perfect alignment. Probes that interrogate multiple gene transcripts and that contain SNPs with ≥1% minor allele frequency in dbSNP dataset were removed. Probe intensities were log2 transformed, background corrected and quantile normalized. Probe intensity was subtracted by the corresponding probe mean across samples. Gene expression level was summarized as mean intensity across probes within gene, using Gencode version 19. Batch effects of RNA labeling and array hybridization were adjusted using an empirical Bayes method (Johnson et al. 2007 (link)).