Abstract
Members of the Bacillus cereus group are ubiquitously present in the environment and can adapt to a wide range of environmental fluctuations. In bacteria, these adaptive responses are generally mediated by two-component signal transduction systems (TCSs), which consist of a histidine kinase (HK) and its cognate response regulator (RR). With the use of in silico techniques, a complete set of HKs and RRs was recovered from eight completely sequenced B. cereus group genomes. By applying a bidirectional best-hits method combined with gene neighbourhood analysis, a footprint of these proteins was made. Around 40 HK-RR gene pairs were detected in each member of the B. cereus group. In addition, each member contained many HK and RR genes not encoded in pairs ("orphans"). Classification of HKs and RRs based on their enzymic domains together with the analysis of two neighbour-joining trees of these domains revealed putative interaction partners for most of the "orphans". Putative biological functions, including involvement in virulence and host-microbe interactions, were predicted for the B. cereus group HKs and RRs by comparing them with those of B. subtilis and other micro-organisms. Remarkably, B. anthracis appeared to lack specific HKs and RRs and was found to contain many truncated, putatively non-functional, HK and RR genes. It is hypothesized that specialization of B. anthracis as a pathogen could have reduced the range of environmental stimuli to which it is exposed. This may have rendered some of its TCSs obsolete, ultimately resulting in the deletion of some HK and RR genes.
