Monitoring the prolonged Tnf stimulation in space and time with topological-functional networks

Genes in linear proximity often share regulatory inputs, expression and evolutionary patterns, even in complex eukaryote genomes with extensive intergenic sequences. Gene regulation, on the other hand, is effected through the co-ordinated activation (or suppression) of genes participating in common biological pathways, which are often transcribed from distant loci. Existing approaches for the study of gene expression focus on the functional aspect, taking positional constraints into account only marginally. In this work we propose a novel concept for the study of gene expression, through the combination of topological and functional information into bipartite networks. Starting from genome-wide expression profiles, we define extended chromosomal regions with consistent patterns of differential gene expression and then associate these domains with enriched functional pathways. By analyzing the resulting networks in terms of size, connectivity and modularity we can draw conclusions on the way genome organization may underlie the gene regulation program. Implementation of this approach in a detailed RNASeq profiling of sustained Tnf stimulation of mouse synovial fibroblasts, allowed us to identify unexpected regulatory changes taking place in the cells after 24 h of stimulation. Bipartite network analysis suggests that the cytokine response set by Tnf, progresses through two distinct transitions. An early generalization of the inflammatory response, that is followed by a late shutdown of immune-related functions and the redistribution of expression to developmental and cell adhesion pathways and distinct chromosomal regions. We show that the incorporation of topological information may provide additional insights in the complex propagation of Tnf activation.