Defining hierarchical protein interaction networks from spectral analysis of bacterial proteomes

Cellular phenotypes emerge from a hierarchy of molecular interactions: proteins interact to form complexes, pathways, and phenotypes. We show that hierarchical networks of protein interactions can be extracted from the statistical pattern of proteome variation as measured across thousands of bacteria and that these hierarchies reflect the emergence of complex bacterial phenotypes. We describe the mathematics underlying our statistical approach and validate our results through gene-set enrichment analysis and comparison to existing experimentally-derived hierarchical databases. We demonstrate the biological utility of our unbiased hierarchical models by creating a model of motility in Pseudomonas aeruginosa and using it to discover a previously unappreciated genetic effector of twitch-based motility. Overall, our approach, SCALES (Spectral Correlation Analysis of Layered Evolutionary Signals), predicts hierarchies of protein interaction networks describing emergent biological function using only the statistical pattern of bacterial proteome variation.