A range of tools are available to model, simulate and analyze gene regulatory networks (GRNs). However, these tools provide limited ability to control network topology, system dynamics, design of experiments, data properties, or noise characteristics. Independent control of these properties is the key to drawing conclusions on which inference method to use and what result to expect from it, as well as obtaining desired approximations of real biological systems. To draw conclusions on the relation between a network or data property and the performance of an inference method in simulations, system approximations with varying properties are needed. We present a Matlab package \gs for generation and analysis of networks and data in a dynamical systems framework with focus on the ability to vary properties. It supplies not only essential components that have been missing, but also wrappers to existing tools in common use. In particular, it contains tools for controlling and analyzing network topology (random, small-world, scale-free), stability of linear time-invariant systems, signal to noise ratio (SNR), and Interampatteness. It also contains methods for design of perturbation experiments, bootstrapping, analysis of linear dependence, sample selection, scaling of the SNR, and performance evaluation. GeneSPIDER offers control of network and data properties in simulations, together with tools to analyze these properties and draw conclusions on the quality of inferred GRNs. It can be fetched freely from the online =git= repository https://bitbucket.org/sonnhammergrni/genespider.
To understand how the components of a complex system like the biological cell interact and regulate each other, we need to collect data for how the components respond to system perturbations. Such data can then be used to solve the inverse problem of inferring a network that describes how the pieces influence each other. The work in this thesis deals with modelling the cell regulatory system, often represented as a network, with tools and concepts derived from systems biology. The first investigation focuses on network sparsity and algorithmic biases introduced by penalised network inference procedures. Many contemporary network inference methods rely on a sparsity parameter such as the L1 penalty term used in the LASSO. However, a poor choice of the sparsity parameter can give highly incorrect network estimates. In order to avoid such poor choices, we devised a method to optimise the sparsity parameter, which maximises the accuracy of the inferred network. We showed that it is effective on in silico data sets with a reasonable level of informativeness and demonstrated that accurate prediction of network sparsity is key to elucidate the correct network parameters. The second investigation focuses on how knowledge from association networks can be transferred to regulatory network inference procedures. It is common that the quality of expression data is inadequate for reliable gene regulatory network inference. Therefore, we constructed an algorithm to incorporate prior knowledge and demonstrated that it increases the accuracy of network inference when the quality of the data is low. The third investigation aimed to understand the influence of system and data properties on network inference accuracy. L1 regularisation methods commonly produce poor network estimates when the data used for inference is ill-conditioned, even when the signal to noise ratio is so high that all links in the network can be proven to exist for the given significance. In this study we elucidated some general principles for under what conditions we expect strongly degraded accuracy. Moreover, it allowed us to estimate expected accuracy from conditions of simulated data, which was used to predict the performance of inference algorithms on biological data. Finally, we built a software package GeneSPIDER for solving problems encountered during previous investigations. The software package supports highly controllable network and data generation as well as data analysis and exploration in the context of network inference.
At the time of the doctoral defense, the following paper was unpublished and had a status as follows: Paper 4: Manuscript.