Abstract
BACKGROUND: Chronic lymphocytic leukemia (CLL) is an incurable malignancy of mature B-lymphocytes, characterized as being a heterogeneous disease with variable clinical manifestation and survival. Mutational statuses of rearranged immunoglobulin heavy chain variable (IGVH) genes has been consider one of the most important prognostic factors in CLL, but despite of its proven value to predict the course of the disease, the regulatory programs and biological mechanisms responsible for the differences in clinical behavior are poorly understood. METHODS: In this study, (i) we performed differential gene expression analysis between the IGVH statuses using multiple and independent CLL cohorts in microarrays platforms, based on this information, (ii) we constructed a simplified protein-protein interaction (PPI) network and (iii) investigated its structure and critical genes. This provided the basis to (iv) develop a Boolean model, (v) infer biological regulatory mechanism and (vi) performed perturbation simulations in order to analyze the network in dynamic state. RESULTS: The result of topological analysis and the Boolean model showed that the transcriptional relationships of IGVH mutational status were determined by specific regulatory proteins (PTEN, FOS, EGR1, TNF, TGFBR3, IFGR2 and LPL). The dynamics of the network was controlled by attractors whose genes were involved in multiple and diverse signaling pathways, which may suggest a variety of mechanisms related with progression occurring over time in the disease. The overexpression of FOS and TNF fixed the fate of the system as they can activate important genes implicated in the regulation of process of adhesion, apoptosis, immune response, cell proliferation and other signaling pathways related with cancer. CONCLUSION: The differences in prognosis prediction of the IGVH mutational status are related with several regulatory hubs that determine the dynamic of the system.