Functional and evolutionary inference in gene networks: Does topology matter?

The relationship between the topology of a biological network and its functional or evolutionary properties has attracted much recent interest. It has been suggested that most, if not all, biological networks are 'scale free.' That is, their connections follow power-law distributions, such that there are very few nodes with very many connections and vice versa. The number of target genes of known transcriptional regulators in the yeast, Saccharomyces cerevisiae, appears to follow such a distribution, as do other networks, such as the yeast network of protein-protein interactions. These findings have inspired attempts to draw biological inferences from general properties associated with scale-free network topology. One often cited general property is that, when compromised, highly connected nodes will tend to have a larger effect on network function than sparsely connected nodes. For example, more highly connected proteins are more likely to be lethal when knocked out. However, the correlation between lethality and connectivity is relatively weak, and some highly connected proteins can be removed without noticeable phenotypic effect. Similarly, network topology only weakly predicts the response of gene expression to environmental perturbations. Evolutionary simulations of gene-regulatory networks, presented here, suggest that such weak or non-existent correlations are to be expected, and are likely not due to inadequacy of experimental data. We argue that 'top-down' inferences of biological properties based on simple measures of network topology are of limited utility, and we present simulation results suggesting that much more detailed information about a gene's location in a regulatory network, as well as dynamic gene-expression data, are needed to make more meaningful functional and evolutionary predictions. Specifically, we find in our simulations that: (1) the relationship between a gene's connectivity and its fitness effect upon knockout depends on its equilibrium expression level; (2) correlation between connectivity and genetic variation is virtually non-existent, yet upon independent evolution of networks with identical topologies, some nodes exhibit consistently low or high polymorphism; and (3) certain genes show low polymorphism yet high divergence among independent evolutionary runs. This latter pattern is generally taken as a signature of positive selection, but in our simulations its cause is often neutral coevolution of regulatory inputs to the same gene.