Characterization of Whole Blood Gene Expression Profiles as a Sequel to Globin mRNA Reduction in Patients with Sickle Cell Disease

Global transcriptome analysis of whole blood RNA using microarrays has been proven to be challenging due to the high abundance of globin transcripts that constitute 70% of whole blood mRNA. This is a particular problem in patients with sickle cell disease, secondary to the high abundance of globin-expressing nucleated red blood cells and reticulocytes in the circulation. In order to accurately measure the steady state blood transcriptome in sickle cell patients we evaluated the efficacy of reducing globin transcripts in PAXgene stabilized RNA for genome-wide transcriptome analyses using microarrays. We demonstrate here by both microarrays and Q-PCR that the globin mRNA depletion method resulted in 5565 fold reduction in globin transcripts in whole blood collected from healthy volunteers and sickle cell disease patients. This led to an improvement in microarray data quality by reducing data variability, with increased detection rate of expressed genes and improved overlap with the expression signatures of isolated peripheral blood mononuclear (PBMC) preparations. Analysis of differences between the whole blood transcriptome and PBMC transcriptome revealed important erythrocyte genes that participate in sickle cell pathogenesis and compensation. The combination of globin mRNA reduction after whole-blood RNA stabilization represents a robust clinical research methodology for the discovery of biomarkers for hematologic diseases.