Kinetic modelling and meta-analysis of the B. subtilis SigA regulatory network during spore germination and outgrowth

Abstract

This study describes the meta-analysis and kinetic modelling of gene expression control by sigma factor SigA of Bacillus subtilis during germination and outgrowth based on microarray data from 14 time points. The analysis computationally models the direct interaction among SigA, SigA-controlled sigma factor genes (sigM, sigH, sigD, sigX), and their target genes. Of the >800 known genes in the SigA regulon, as extracted from databases, 311 genes were analysed, and 190 were confirmed by the kinetic model as being controlled by SigA. For the remaining genes, alternative regulators satisfying kinetic constraints were suggested. The kinetic analysis suggested another 214 genes as potential SigA targets. The modelling was able to (i) create a particular SigA-controlled gene expression network that is active under the conditions for which the expression time series was obtained, and where SigA is the dominant regulator, (ii) suggest new potential SigA target genes, and (iii) find other possible regulators of a given gene or suggest a new mechanism of its control by identifying a matching profile of unknown regulator(s). Selected predicted regulatory interactions were experimentally tested, thus validating the model.