Overview of ChIP-seq results.

(A). Density of all TFs (green) and ChIPed TFs (orange) in this study throughout the P. aeruginosa genome. (B). Annotation heatmap of all peak distribution with 6 locations: Upstream, OverlapStart, Inside, OverlapEnd, Downstream, and IncludeFeature. (C). Peak distance to the translational start site (TSS) of each DBD family. (D). Treemap of the 172 TFs peak numbers based on DBD family. Each box’s size represents the family’s size (number of peaks), and the explained variance of each DBD type means the color shades of each box. DBD families of ChIPed TFs are classified into 20 different categories: LuxR, LysR, Two DBDs (Two DNA binding domains), AraC, TetR, ArsR, CRP, OmpR, GntR, MarR, AsnC, Cro/CI, TyrR, Rrf2, MerR, IclR, Fis, RpiR, DeoR, and undetectable. (E). The dot plot shows the top 10 gene ontology (GO) terms from the PseudoCAP annotation of ChIPed TFs. The size of the dots indicates the significance of each functional category, quantified by −log10 (p. adjust).

Hierarchical networks and network motifs.

(A). Overview of hierarchy regulatory network after removing TFs with degree O+I <= 10. Nodes indicate TFs, and the color represents the hierarchical level. The top level was highlighted in red, the middle in yellow, and the bottom highlighted in blue. The edges with arrows indicate the regulatory direction. Gray means downward-pointing, and red means upward-pointing. (B). The auto-regulator motif of 9 TFs. (C). Three transcription-factor motifs occurrences with five basic triangular motifs and eight toggle switch motifs. The circle presents TF, and the arrow indicates the regulatory direction. (D). Alluvial diagram reveals basic triangle motif 3 (n = 1956) depending on DBD types. The color of splines is highlighted in different DBD families, and the name of DBD families are labeled.

Core co-association regulatory networks.

(A). Core clusters of significant co-binding patterns of TFs. Each TF is highlighted in different colors based on DBD types. The co-association score by pair of TFs was calculated by Jaccard statistics, which measures the ratio of the number of base pairs in overlapped binding peaks on both TFs to the number of base pairs in their union. (B). The histogram’s overlapped target genes of TFs in cluster 3 represent the number of target peaks in the individual/overlapped set. (C). The network of co-regulation of TFs in cluster 3 with co-bound targets of more than 4. (D-E). Genome browser view of TFs in cluster 3 binding intensities at the PA2504 and pqsH locus.

Newly identified virulence-related master regulators.

(A). Overview of all identified master regulators related to 6 pathways, including QS, motility, biofilm, siderophore, pyocyanin, and ROS. Each circle represents one TF, with the height indicating the significance, quantified by -log10 (p.adjust), and the size indicates the number of targets associated with the virulence pathway. (B). Intersection of master regulators in 6 virulent pathways. The bar chart on the top shows the number of intersections, while the matrix below indicates which TFs are involved in specific biological processes such as siderophore production, pyocyanin production, biofilm formation, ROS response, QS, and motility.

Overview of transcriptional network of TFs in P. aeruginosa.

(A). The interaction network of virulence-related master regulators in P. aeruginosa. The ten virulence pathways are highlighted in different colors. The size of the circle represents the degree, and the width of the edges indicates the overlapped number of TFs between two pathways. (B). The target annotated using the COG database shows four orthology classification profiling from different DBD types TFs. The size of rectangle indicates the number of targets.

Conservation and variability of TFs in PAO1.

(A). The pie chart shows the proportions for all genes. (B). The pie chart of proportion for all TFs identified from PAO1. (C). The bar plot of proportion for non-core TFs. Core genes indicate genes existed in 99% ∼ 100% strains; Soft core genes indicate genes existed in 95% ∼ 99% strains; Shell genes indicate genes existed in 15% ∼ 95% strains; Cloud genes indicate genes existed in 0% ∼ 15% strains. (D). The conservation and evolutionary relationship of all 373 TFs in PAO1 among bacteria, archaebacteria, fungi, plants, and animals. The conservation value was normalized after blastp alignment82. The phylogenetic trees were constructed using MEGA1183 and plotted via R package ggtree84.