Abstract
Pseudomonas aeruginosa possesses one of the most complex bacterial regulatory networks, which largely contributes to its success as a pathogen. However, most of its transcription factors (TFs) are still uncharacterized and the potential intra-species variability in regulatory networks has been mostly ignored so far. Here, we used DAP-seq to map the genome-wide binding sites of all 55 DNA-binding two-component systems (TCSs) response regulators (RRs) across the three major P. aeruginosa lineages. The resulting networks encompass about 40% of all genes in each strain and contain numerous new regulatory interactions across most major physiological processes. Strikingly, about half of the detected targets are specific to only one or two strains, revealing a previously unknown large functional diversity of TFs within a single species. Three main mechanisms were found to drive this diversity, including differences in accessory genome content, as exemplified by the strain-specific plasmid in IHMA87 outlier strain which harbors numerous binding sites of conserved chromosomally-encoded RRs. Additionally, most RRs display potential auto-regulation or RR-RR cross-regulation, bringing to light the vast complexity of this network. Overall, we provide the first complete delineation of the TCSs regulatory network in P. aeruginosa that will represent an important resource for future studies on this pathogen.