Differential gene expression analysis of sickle cell anemia in steady and crisis state

Sickle cell anemia is one of the most prevalent genetic diseases worldwide, showing great clinical heterogeneity. This study compared the gene expression patterns between sickle cell anemia pediatric patients in steady state and in crisis state, as compared to age-paired, healthy individuals. RNA sequencing was performed from these groups of patients/controls using Illumina HiSeq 2500 equipment. The resulting differentially expressed genes were loaded into QIAGEN's ingenuity pathway analysis. The results showed that EIF2 pathway and NRF2-mediated oxidative stress-response pathways were more highly activated both in steady state and in crisis patients, as compared to healthy individuals. In addition, we found increased activation of eIF4 and p70S6K signaling pathways in crisis state compared to healthy individuals. The transcription factor GATA-1 was found exclusively in steady state while SPI was found exclusively in crisis state. IL6 and VEGFA were found only in crisis state, while IL-1B was found exclusively in steady state. The regulator effects analysis revealed IgG1 as an upstream regulator in steady state compared to healthy individuals, resulting in invasion of prostate cancer cell lines as the disease/function outcome. For crisis-state patients versus healthy individuals, two networks of regulator effects revealed STAT1, CD40LG, TGM2, IRF7, IRF4, and IRF1 acting as upstream regulators, resulting in disease/function outcomes, including engulfment of cells and aggregation of blood cells and inflammation of joints. Our results indicated genes and pathways that can provide clues on the molecular events involved in the severity of sickle cell disease.