Elucidation of regulatory modes for five two-component systems in Escherichia coli reveals novel relationships

Abstract

AbstractEscherichia coli uses two-component systems (TCSs) to respond to environmental signals. TCSs affect gene expression and are parts of E. coli’s global transcriptional regulatory network (TRN). Here, we identified the regulons of five TCSs in E. coli MG1655: BaeSR and CpxAR, which were stimulated by ethanol stress; KdpDE and PhoRB, induced by limiting potassium and phosphate, respectively; and ZraSR, stimulated by zinc. We analyzed RNA-seq data using independent component analysis (ICA). ChIP-exo data was used to validate condition-specific target gene binding sites. Based on this data we (1) identify the target genes for each TCS; (2) show how the target genes are transcribed in response to stimulus; and (3) reveal novel relationships between TCSs, which indicate non-cognate inducers for various response regulators, such as BaeR to iron starvation, CpxR to phosphate limitation, and PhoB and ZraR to cell envelope stress. Our understanding of the TRN in E. coli is thus notably expanded.ImportanceE. coli is a common commensal microbe found in human gut microenvironment; however, some strains cause diseases like diarrhea, urinary tract infections and meningitis. E. coli’s two-component system (TCS) modulates target gene expression, specially related to virulence, pathogenesis and anti-microbial peptides, in response to environmental stimuli. Thus, it is of utmost importance to understand the transcriptional regulation of the TCSs to infer its environmental adaptation and disease pathogenicity. Utilizing a combinatorial approach integrating RNAseq, independent component analysis, ChIP-exo and data mining, we show that TCSs have five different modes of transcriptional regulation. Our data further highlights non-cognate inducers of TCSs emphasizing cross-regulatory nature of TCSs in E. coli and suggests that TCSs may have a role beyond their cognate functionalities. In summary, these results when further incorporated with genome scale metabolic models can lead to understanding of metabolic capabilities of bacteria and correctly predict complex phenotype under diverse conditions.