Detecting cellular reprogramming determinants by differential stability analysis of gene regulatory networks

Cellular differentiation and reprogramming are processes that are carefully orchestrated by the activation and repression of specific sets of genes. An increasing amount of experimental results show that despite the large number of genes participating in transcriptional programs of cellular phenotypes, only few key genes, which are coined here as reprogramming determinants, are required to be directly perturbed in order to induce cellular reprogramming. However, identification of reprogramming determinants still remains a combinatorial problem, and the state-of-art methods addressing this issue rests on exhaustive experimentation or prior knowledge to narrow down the list of candidates.

Here we present a computational method, without any preliminary selection of candidate genes, to identify reduced subsets of genes, which when perturbed can induce transitions between cellular phenotypes. The method relies on the expression profiles of two stable cellular phenotypes along with a topological analysis stability elements in the gene regulatory network that are necessary to cause this multi-stability. Since stable cellular phenotypes can be considered as attractors of gene regulatory networks, cell fate and cellular reprogramming involves transition between these attractors, and therefore current method searches for combinations of genes that are able to destabilize a specific initial attractor and stabilize the final one in response to the appropriate perturbations.

The method presented here represents a useful framework to assist researchers in the field of cellular reprogramming to design experimental strategies with potential applications in the regenerative medicine and disease modelling.