Comprehensive mapping of the Escherichia coli flagellar regulatory network

Motility-activating mutations upstream of flhDC reduce acid shock survival of Escherichia coli

Many neutralophilic bacterial species try to evade acid stress with an escape strategy, which is reflected in the increased expression of genes coding for flagellar components. Extremely acid-tolerant bacteria, such as Escherichia coli, survive the strong acid stress, e.g., in the stomach of vertebrates. Recently, we were able to show that the induction of motility genes in E. coli is strictly dependent on the degree of acid stress, i.e., they are induced under mild acid stress but not under severe acid stress. However, it was not known to what extent fine-tuned expression of motility genes is related to fitness and the ability to survive periods of acid shock. In this study, we demonstrate that the expression of FlhDC, the master regulator of flagellation, is inversely correlated with the acid shock survival of E. coli. We encountered this phenomenon when analyzing mutants from the Keio collection, in which the expression of flhDC was altered by an insertion sequence element. These results suggest a fitness trade-off between acid tolerance and motility.IMPORTANCEEscherichia coli is extremely acid-resistant, which is crucial for survival in the gastrointestinal tract of vertebrates. Recently, we systematically studied the response of E. coli to mild and severe acidic conditions using Ribo-Seq and RNA-Seq. We found that motility genes are induced at pH 5.8 but not at pH 4.4, indicating stress-dependent synthesis of flagellar components. In this study, we demonstrate that motility-activating mutations upstream of flhDC, encoding the master regulator of flagella genes, reduce the ability of E. coli to survive periods of acid shock. Furthermore, we show an inverse correlation between motility and acid survival using a chromosomal isopropyl β-D-thio-galactopyranoside (IPTG)-inducible flhDC promoter and by sampling differentially motile subpopulations from swim agar plates. These results reveal a previously undiscovered trade-off between motility and acid tolerance and suggest a differentiation of E. coli into motile and acid-tolerant subpopulations, driven by the integration of insertion sequence elements.

TusDCB, a sulfur transferase complex involved in tRNA modification, contributes to UPEC pathogenicity

tRNA modifications play a crucial role in ensuring accurate codon recognition and optimizing translation levels. While the significance of these modifications in eukaryotic cells for maintaining cellular homeostasis and physiological functions is well-established, their physiological roles in bacterial cells, particularly in pathogenesis, remain relatively unexplored. The TusDCB protein complex, conserved in γ-proteobacteria like Escherichia coli, is involved in sulfur modification of specific tRNAs. This study focused on the role of TusDCB in the virulence of uropathogenic E. coli (UPEC), a bacterium causing urinary tract infections. The findings indicate that TusDCB is essential for optimal production of UPEC's virulence factors, including type 1 fimbriae and flagellum, impacting the bacterium's ability to aggregate in bladder epithelial cells. Deletion of tusDCB resulted in decreased virulence against urinary tract infection mice. Moreover, mutant TusDCB lacking sulfur transfer activity and tusE- and mnmA mutants revealed the indispensability of TusDCB's sulfur transfer activity for UPEC pathogenicity. The study extends its relevance to highly pathogenic, multidrug-resistant strains, where tusDCB deletion reduced virulence-associated bacterial aggregation. These insights not only deepen our understanding of the interplay between tRNA sulfur modification and bacterial pathogenesis but also highlight TusDCB as a potential therapeutic target against UPEC strains resistant to conventional antimicrobial agents.

FlhE functions as a chaperone to prevent formation of periplasmic flagella in Gram-negative bacteria

The bacterial flagellum is an organelle utilized by many Gram-negative bacteria to facilitate motility. The flagellum is composed of a several µm long, extracellular filament that is connected to a cytoplasmic rotor-stator complex via a periplasmic rod. Composed of ∼20 structural proteins, ranging from a few subunits to several thousand building blocks, the flagellum is a paradigm of a complex macromolecular structure that utilizes a highly regulated assembly process. This process is governed by multiple checkpoints that ensure an ordered gene expression pattern coupled to the assembly of the various flagellar building blocks in order to produce a functional flagellum. Using epifluorescence, super-resolution STED and transmission electron microscopy, we discovered that in Salmonella , the absence of one periplasmic protein, FlhE, prevents proper flagellar morphogenesis and results in the formation of periplasmic flagella. The periplasmic flagella disrupt cell wall synthesis, leading to a loss of the standard cell morphology resulting in cell lysis. We propose a model where FlhE functions as a periplasmic chaperone to control assembly of the periplasmic rod to prevent formation of periplasmic flagella. Our results highlight that bacteria evolved sophisticated regulatory mechanisms to control proper flagellar assembly and minor deviations from this highly regulated process can cause dramatic physiological consequences.

Flagella-mediated adhesion of Escherichia coli O157:H7 to surface of stainless steel, glass and fresh produces during sublethal injury and recovery

E. coli O157:H7 can be induced into sublethally injured (SI) state by lactic acid (LA) and regain activity in nutrient environments. This research clarified the role of flagella-related genes (fliD, fliS, cheA and motA) in adhesion of E. coli O157:H7 onto stainless steel, glass, lettuce, spinach, red cabbage and cucumber during LA-induced SI and recovery by plate counting. Results of adhesion showed improper flagellar rotation caused by the deletion of motA resulting in the decreased adhesion. Motility of wildtype determined by diameter of motility halo decreased in SI state and repaired with recovery time increasing, lagging behind changes in expression of flagella-related genes. Flagellar function-impaired strains all exhibited non-motile property. Thus, we speculated that flagella-mediated motility is critical in early stage of adhesion. We also found the effects of Fe2+, Ca2+ and Mn2+ on adhesion or motility of wildtype was independent of bacterial states. However, the addition of Ca2+ and Mn2+ did not affect motility of flagellar function-impaired strains as they did on wildtype. This research provides new insights to understand the role of flagella and cations in bacterial adhesion, which will aid in development of anti-adhesion agents to reduce bio-contamination in food processing.

Bacterial RNase III: Targets and physiology

Bacteria can rapidly adapt to changes in their environment thanks to the innate flexibility of their genetic expression. The high turnover rate of RNAs, in particular messenger and regulatory RNAs, provides an important contribution to this dynamic adjustment. Recycling of RNAs is ensured by ribonucleases, among which RNase III is the focus of this review. RNase III enzymes are highly conserved from prokaryotes to eukaryotes and have the specific ability to cleave double-stranded RNAs. The role of RNase III in bacterial physiology has remained poorly explored for a long time. However, transcriptomic approaches recently uncovered a large impact of RNase III in gene expression in a wide range of bacteria, generating renewed interest in the physiological role of RNase III. In this review, we first describe the RNase III targets identified from global approaches in 8 bacterial species within 4 Phyla. We then present the conserved and unique functions of bacterial RNase III focusing on growth, resistance to stress, biofilm formation, motility and virulence. Altogether, this review highlights the underestimated impact of RNase III in bacterial adaptation.

A comprehensive list of genes required for the efficient conjugation of plasmid Rts1 was determined by systematic deletion analysis

While conjugation-related genes have been identified in many plasmids by genome sequencing, functional analyses have not yet been performed in most cases, and a full set of conjugation genes has been identified for only a few plasmids. Rts1, a prototype IncT plasmid, is a conjugative plasmid that was originally isolated from Proteus vulgaris. Here, we conducted a systematic deletion analysis of Rts1 to fully understand its conjugation system. Through this analysis along with complementation assays, we identified 32 genes that are required for the efficient conjugation of Rts1 from Escherichia coli to E. coli. In addition, the functions of the 28 genes were determined or predicted; 21 were involved in mating-pair formation, three were involved in DNA transfer and replication, including a relaxase gene belonging to the MOBH12 family, one was involved in coupling, and three were involved in transcriptional regulation. Among the functionally well-analysed conjugation systems, most of the 28 genes showed the highest similarity to those of the SXT element, which is an integrative conjugative element of Vibrio cholerae. The Rts1 conjugation gene set included all 23 genes required for the SXT system. Two groups of plasmids with conjugation systems nearly identical or very similar to that of Rts1 were also identified.

Comparative genomic analysis of uropathogenic Escherichia coli strains from women with recurrent urinary tract infection

Introduction

Recurrent urinary tract infections (RUTIs) caused by uropathogenic Escherichia coli are costly public health problems impacting patients' quality of life.

Aim

In this work, a comparative genomics analysis of three clinical RUTI strains isolated from bladder biopsy specimens was performed.

Materials and methods

One hundred seventy-two whole genomes of urinary tract E. coli strains were selected from the NCBI database. The search for virulence factors, fitness genes, regions of interest, and genetic elements associated with resistance was manually carried out. The phenotypic characterization of antibiotic resistance, haemolysis, motility, and biofilm formation was performed. Moreover, adherence and invasion assays with human bladder HTB-5 cells, and transmission electron microscopy (TEM) were performed.

Results

The UTI-1_774U and UTI-3_455U/ST1193 strains were associated with the extraintestinal pathotypes, and the UTI-2_245U/ST295 strain was associated with the intestinal pathotype, according to a phylogenetic analysis of 172 E. coli urinary strains. The three RUTI strains were of clinical, epidemiological, and zoonotic relevance. Several resistance genes were found within the plasmids of these strains, and a multidrug resistance phenotype was revealed. Other virulence genes associated with CFT073 were not identified in the three RUTI strains (genes for type 1 and P fimbriae, haemolysin hlyA, and sat toxin). Quantitative adherence analysis showed that UTI-1_774U was significantly (p < 0.0001) more adherent to human bladder HTB-5 cells. Quantitative invasion analysis showed that UTI-2_245U was significantly more invasive than the control strains. No haemolysis or biofilm activity was detected in the three RUTI strains. The TEM micrographs showed the presence of short and thin fimbriae only in the UTI-2_245U strain.

Conclusion

The high variability and genetic diversity of the RUTI strains indicate that are a mosaic of virulence, resistance, and fitness genes that could promote recurrence in susceptible patients.

σE of Streptomyces coelicolor can function both as a direct activator or repressor of transcription

σ factors are considered as positive regulators of gene expression. Here we reveal the opposite, inhibitory role of these proteins. We used a combination of molecular biology methods and computational modeling to analyze the regulatory activity of the extracytoplasmic σE factor from Streptomyces coelicolor. The direct activator/repressor function of σE was then explored by experimental analysis of selected promoter regions in vivo. Additionally, the σE interactome was defined. Taken together, the results characterize σE, its regulation, regulon, and suggest its direct inhibitory function (as a repressor) in gene expression, a phenomenon that may be common also to other σ factors and organisms.

In vivo RNA interactome profiling reveals 3'UTR-processed small RNA targeting a central regulatory hub

Small noncoding RNAs (sRNAs) are crucial regulators of gene expression in bacteria. Acting in concert with major RNA chaperones such as Hfq or ProQ, sRNAs base-pair with multiple target mRNAs and form large RNA-RNA interaction networks. To systematically investigate the RNA-RNA interactome in living cells, we have developed a streamlined in vivo approach iRIL-seq (intracellular RIL-seq). This generic approach is highly robust, illustrating the dynamic sRNA interactomes in Salmonella enterica across multiple stages of growth. We have identified the OmpD porin mRNA as a central regulatory hub that is targeted by a dozen sRNAs, including FadZ cleaved from the conserved 3'UTR of fadBA mRNA. Both ompD and FadZ are activated by CRP, constituting a type I incoherent feed-forward loop in the fatty acid metabolism pathway. Altogether, we have established an approach to profile RNA-RNA interactomes in live cells, highlighting the complexity of RNA regulatory hubs and RNA networks.

Hierarchy and interconnected networks in the WhiB7 mediated transcriptional response to antibiotic stress in Mycobacterium abscessus

Mycobacterium abscessus is intrinsically resistant to antibiotics effective against other pathogenic mycobacteria largely due to the drug-induced expression of genes that confer resistance. WhiB7 is a major hub controlling the induction of resistance to ribosome-targeting antibiotics. It activates the expression of >100 genes, 7 of which are known determinants of drug resistance; the function of most genes within the regulon is however unknown, but some conceivably encode additional mechanisms of resistance. Furthermore, the hierarchy of gene expression within the regulon, if any, is poorly understood. In the present work we have identified 56 WhiB7 binding sites using chromatin immunoprecipitation sequencing (CHIP-Seq) which accounts for the WhiB7-dependent upregulation of 72 genes, and find that M. abscessus WhiB7 functions exclusively as a transcriptional activator at promoters recognized by σA/σB. We have investigated the role of 18 WhiB7 regulated genes in drug resistance. Our results suggest that while some genes within the regulon (eg. erm41, hflX, eis2 and the ABCFs) play a major role in resistance, others make smaller contributions (eg. MAB_4324c and MAB_1409c) and the observed hypersensitivity ΔMabwhiB7 is a cumulative effect of these individual contributions. Moreover, our CHIP-Seq data implicate additional roles of WhiB7 induced genes beyond antibiotic resistance. Finally, we identify a σH-dependent network in aminoglycoside and tigecycline resistance which is induced upon drug exposure and is further activated by WhiB7 demonstrating the existence of a crosstalk between components of the WhiB7-dependent and -independent circuits.

Cross-talk between QseBC and PmrAB two-component systems is crucial for regulation of motility and colistin resistance in Enteropathogenic Escherichia coli

The quorum sensing two-component system (TCS) QseBC has been linked to virulence, motility and metabolism regulation in multiple Gram-negative pathogens, including Enterohaemorrhagic Escherichia coli (EHEC), Uropathogenic E. coli (UPEC) and Salmonella enterica. In EHEC, the sensor histidine kinase (HK) QseC detects the quorum sensing signalling molecule AI-3 and also acts as an adrenergic sensor binding host epinephrine and norepinephrine. Downstream changes in gene expression are mediated by phosphorylation of its cognate response regulator (RR) QseB, and 'cross-talks' with non-cognate regulators KdpE and QseF to activate motility and virulence. In UPEC, cross-talk between QseBC and TCS PmrAB is crucial in the regulation and phosphorylation of QseB RR that acts as a repressor of multiple pathways, including motility. Here, we investigated QseBC regulation of motility in the atypical Enteropathogenic E. coli (EPEC) strain O125ac:H6, causative agent of persistent diarrhoea in children, and its possible cross-talk with the KdpDE and PmrAB TCS. We showed that in EPEC QseB acts as a repressor of genes involved in motility, virulence and stress response, and in absence of QseC HK, QseB is likely activated by the non-cognate PmrB HK, similarly to UPEC. We show that in absence of QseC, phosphorylated QseB activates its own expression, and is responsible for the low motility phenotypes seen in a QseC deletion mutant. Furthermore, we showed that KdpD HK regulates motility in an independent manner to QseBC and through a third unidentified party different to its own response regulator KdpE. We showed that PmrAB has a role in iron adaptation independent to QseBC. Finally, we showed that QseB is the responsible for activation of colistin and polymyxin B resistance genes while PmrA RR acts by preventing QseB activation of these resistance genes.

Promoter selectivity of the RhlR quorum-sensing transcription factor receptor in Pseudomonas aeruginosa is coordinated by distinct and overlapping dependencies on C4-homoserine lactone and PqsE

Quorum sensing is a mechanism of bacterial cell-cell communication that relies on the production and detection of small molecule autoinducers, which facilitate the synchronous expression of genes involved in group behaviors, such as virulence factor production and biofilm formation. The Pseudomonas aeruginosa quorum sensing network consists of multiple interconnected transcriptional regulators, with the transcription factor, RhlR, acting as one of the main drivers of quorum sensing behaviors. RhlR is a LuxR-type transcription factor that regulates its target genes when bound to its cognate autoinducer, C4-homoserine lactone, which is synthesized by RhlI. RhlR function is also regulated by the metallo-β-hydrolase enzyme, PqsE. We recently showed that PqsE binds RhlR to alter its affinity for promoter DNA, a new mechanism of quorum-sensing receptor activation. Here, we perform ChIP-seq analyses of RhlR to map the binding of RhlR across the P. aeruginosa genome, and to determine the impact of C4-homoserine lactone and PqsE on RhlR binding to different sites across the P. aeruginosa genome. We identify 40 RhlR binding sites, all but three of which are associated with genes known to be regulated by RhlR. C4-homoserine lactone is required for maximal binding of RhlR to many of its DNA sites. Moreover, C4-homoserine lactone is required for maximal RhlR-dependent transcription activation from all sites, regardless of whether it impacts RhlR binding to DNA. PqsE is required for maximal binding of RhlR to many DNA sites, with similar effects on RhlR-dependent transcription activation from those sites. However, the effects of PqsE on RhlR specificity are distinct from those of C4-homoserine lactone, and PqsE is sufficient for RhlR binding to some DNA sites in the absence of C4-homoserine lactone. Together, C4-homoserine lactone and PqsE are required for RhlR binding at the large majority of its DNA sites. Thus, our work reveals three distinct modes of activation by RhlR: i) when RhlR is unbound by autoinducer but bound by PqsE, ii) when RhlR is bound by autoinducer but not bound by PqsE, and iii) when RhlR is bound by both autoinducer and PqsE, establishing a stepwise mechanism for the progression of the RhlR-RhlI-PqsE quorum sensing pathway in P. aeruginosa.

σ28-dependent small RNA regulation of flagella biosynthesis

Flagella are important for bacterial motility as well as for pathogenesis. Synthesis of these structures is energy intensive and, while extensive transcriptional regulation has been described, little is known about the posttranscriptional regulation. Small RNAs (sRNAs) are widespread posttranscriptional regulators, most base pairing with mRNAs to affect their stability and/or translation. Here, we describe four UTR-derived sRNAs (UhpU, MotR, FliX and FlgO) whose expression is controlled by the flagella sigma factor σ28 (fliA) in Escherichia coli. Interestingly, the four sRNAs have varied effects on flagellin protein levels, flagella number and cell motility. UhpU, corresponding to the 3´ UTR of a metabolic gene, likely has hundreds of targets including a transcriptional regulator at the top flagella regulatory cascade connecting metabolism and flagella synthesis. Unlike most sRNAs, MotR and FliX base pair within the coding sequences of target mRNAs and act on ribosomal protein mRNAs connecting ribosome production and flagella synthesis. The study shows how sRNA-mediated regulation can overlay a complex network enabling nuanced control of flagella synthesis.

Mechanism of Escherichia coli Lethality Caused by Overexpression of flhDC, the Flagellar Master Regulator Genes, as Revealed by Transcriptome Analysis

The flhDC operon of Escherichia coli encodes a transcription factor that initiates flagella synthesis, elevates flagella construction and enhances cell motility, which all are energetically costly and highly regulated processes. In this study, we found that overexpression of flhDC genes from a strong regulatable pN15E6 plasmid could inhibit the growth of E. coli host cells and even eventually cause death. We used transcriptome analysis to investigate the mechanism of flhDC overexpression lethal to host bacteria. The results showed that a total of 568 differentially expressed genes (DEGs), including 378 up-regulated genes and 190 down-regulated genes were detected when the flhDC genes were over-expressed. Functional enrichment analysis results showed that the DEGs are related to a series of crucial biomolecular processes, including flagella synthesis, oxidative phosphorylation and pentose phosphate pathways, etc. We then examined, using RT-qPCR, the expression of key genes of the oxidative phosphorylation pathway at different time points after induction. Results showed that their expression increased in the early stage and decreased afterward, which was suggested to be the result of feedback on the overproduction of ROS, a strong side effect product of the elevated oxidative phosphorylation process. To further verify the level of ROS output, flhDC over-expressed bacteria cells were stained with DCHF-DA and a fluorescence signal was detected using flow cytometry. Results showed that the level of ROS output was higher in cells with over-expressed flhDC than in normal controls. Besides, we found upregulation of other genes (recN and zwf) that respond to ROS damage. This leads to the conclusion that the bacterial death led by the overexpression of flhDC genes is caused by damage from ROS overproduction, which leaked from the oxidative phosphorylation pathway.

SwrA extends DegU over an UP element to activate flagellar gene expression in Bacillus subtilis

SwrA activates flagellar gene expression in Bacillus subtilis to increase the frequency of motile cells in liquid and elevate flagellar density to enable swarming over solid surfaces. Here we use ChIP-seq to show that SwrA interacts with many sites on the chromosome in a manner that depends on the response regulator DegU. We identify a DegU-specific inverted repeat DNA sequence and show that SwrA synergizes with phosphorylation to increase DegU DNA binding affinity. We further show that SwrA increases the size of the DegU footprint expanding the region bound by DegU towards the promoter. The location of the DegU inverted repeat was critical and moving the binding site closer to the promoter impaired transcription more that could be explained by deactivation. We conclude that SwrA/DegU forms a heteromeric complex that enables both remote binding and interaction between the activator and RNA polymerase in the context of an interceding UP element. We speculate that multimeric activators that resolve cis-element spatial conflicts are common in bacteria and likely act on flagellar biosynthesis loci and other long operons of other multi-subunit complexes.

IMPORTANCE

In Bacteria, the sigma subunit of RNA polymerase recognizes specific DNA sequences called promoters that determine where gene transcription begins. Some promoters also have sequences immediately upstream called an UP element that is bound by the alpha subunit of RNA polymerase and is often necessary for transcription. Finally, promoters may be activated by transcription factors that bind DNA specific sequences and help recruit RNA polymerase to weak promoter elements. Here we show that the promoter for the 32 gene long flagellar operon in Bacillus subtilis requires an UP element and is activated by a heteromeric transcription factor of DegU and SwrA. Our evidence suggests that SwrA oligomerizes DegU over the DNA to allow RNA polymerase to interact with DegU and the UP element simultaneously. Heteromeric activator complexes are known but poorly-understood in bacteria and we speculate they may be needed to resolve spatial conflicts in the DNA sequence.

The PapB/FocB family protein TosR acts as a positive regulator of flagellar expression and is required for optimal virulence of uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) is a major causative agent of urinary tract infections. The bacteria internalize into the uroepithelial cells, where aggregate and form microcolonies. UPEC fimbriae and flagella are important for the formation of microcolonies in uroepithelial cells. PapB/FocB family proteins are small DNA-binding transcriptional regulators consisting of approximately 100 amino acids that have been reported to regulate the expression of various fimbriae, including P, F1C, and type 1 fimbriae, and adhesins. In this study, we show that TosR, a member of the PapB/FocB family is the activator of flagellar expression. The tosR mutant had similar expression levels of type 1, P and F1C fimbriae as the parent strain, but flagellar production was markedly lower than in the parent strain. Flagellin is a major component of flagella. The gene encoding flagellin, fliC, is transcriptionally activated by the sigma factor FliA. The fliA expression is induced by the flagellar master regulator FlhDC. The flhD and flhC genes form an operon. The promoter activity of fliC, fliA and flhD in the tosR mutant was significantly lower than in the parent strain. The purified recombinant TosR does not bind to fliC and fliA but to the upstream region of the flhD gene. TosR is known to bind to an AT-rich DNA sequence consisting of 29 nucleotides. The characteristic AT-rich sequence exists 550-578 bases upstream of the flhD gene. The DNA fragment lacking this sequence did not bind TosR. Furthermore, loss of the tosR gene reduced motility and the aggregation ability of UPEC in urothelial cells. These results indicate that TosR is a transcriptional activator that increases expression of the flhDC operon genes, contributing to flagellar expression and optimal virulence.

Genome-Wide Mapping of the Escherichia coli PhoB Regulon Reveals Many Transcriptionally Inert, Intragenic Binding Sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of, genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert and hence are tolerated as genomic "noise." IMPORTANCE Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of nonregulatory, intragenic binding sites for transcription factors and that these binding sites are not under selective pressure.

Online Omics Platform Expedites Industrial Application of Halomonas bluephagenesis TD1.0

Multi-omic data mining has the potential to revolutionize synthetic biology especially in non-model organisms that have not been extensively studied. However, tangible engineering direction from computational analysis remains elusive due to the interpretability of large datasets and the difficulty in analysis for non-experts. New omics data are generated faster than our ability to use and analyse results effectively, resulting in strain development that proceeds through classic methods of trial-and-error without insight into complex cell dynamics. Here we introduce a user-friendly, interactive website hosting multi-omics data. Importantly, this new platform allows non-experts to explore questions in an industrially important chassis whose cellular dynamics are still largely unknown. The web platform contains a complete KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis derived from principal components analysis, an interactive bio-cluster heatmap analysis of genes, and the Halomonas TD1.0 genome-scale metabolic (GEM) model. As a case study of the effectiveness of this platform, we applied unsupervised machine learning to determine key differences between Halomonas bluephagenesis TD1.0 cultivated under varied conditions. Specifically, cell motility and flagella apparatus are identified to drive energy expenditure usage at different osmolarities, and predictions were verified experimentally using microscopy and fluorescence labelled flagella staining. As more omics projects are completed, this landing page will facilitate exploration and targeted engineering efforts of the robust, industrial chassis H bluephagenesis for researchers without extensive bioinformatics background.

Hierarchy and networks in the transcriptional response of Mycobacterium abscessus to antibiotics

Mycobacterium abscessus causes acute and chronic pulmonary infection in patients with chronic lung damage. It is intrinsically resistance to antibiotics effective against other pathogenic mycobacteria largely due to the drug-induced expression of genes that confer resistance. Induction of genes upon exposure to ribosome targeting antibiotics proceeds via WhiB7-dependent and -independent pathways. WhiB7 controls the expression of >100 genes, a few of which are known determinants of drug resistance. The function of the vast majority of genes within the regulon is unknown, but some conceivably encode additional mechanisms of resistance. Furthermore, the hierarchy of gene expression within the regulon, if any, is poorly understood. In the present work we have identified 56 WhiB7 binding sites using chromatin immunoprecipitation sequencing (CHIP-Seq) which accounts for the WhiB7-dependent upregulation of 70 genes, and find that M. abscessus WhiB7 functions exclusively as a transcriptional activator at promoters recognized by σ A /σ B We have investigated the role of 18 WhiB7 regulated genes in drug resistance and demonstrated the role of MAB_1409c and MAB_4324c in aminoglycoside resistance. Further, we identify a σ H -dependent pathway in aminoglycoside and tigecycline resistance which is induced upon drug exposure and is further activated by WhiB7 demonstrating the existence of a crosstalk between components of the WhiB7-dependent and -independent circuits.

Abstract Importance

The induction of multiple genes that confer resistance to structurally diverse ribosome-targeting antibiotics is funneled through the induction of a single transcriptional activator, WhiB7, by antibiotic-stalled ribosomes. This poses a severe restriction in M. abscessus therapy as treatment with one ribosome-targeting antibiotic confers resistance to all other ribosome-targeting antibiotics. Here we uncover the intricacies of the WhiB7 regulatory circuit, identify three previously unknown determinants of aminoglycoside resistance and unveil a communication between WhiB7 dependent and independent components. This not only expands our understanding of the antibiotic resistance potential of M. abscessus but can also inform the development of much needed therapeutic options.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

The E. coli pathobiont LF82 encodes a unique variant of σ70 that results in specific gene expression changes and altered phenotypes

LF82, an adherent invasive Escherichia coli pathobiont, is associated with ileal Crohn's disease, an inflammatory bowel disease of unknown etiology. Although LF82 contains no virulence genes, it carries several genetic differences, including single nucleotide polymorphisms (SNPs), that distinguish it from nonpathogenic E. coli. We have identified and investigated an extremely rare SNP that is within the highly conserved rpoD gene, encoding σ70, the primary sigma factor for RNA polymerase. We demonstrate that this single residue change (D445V) results in specific transcriptome and phenotypic changes that are consistent with multiple phenotypes observed in LF82, including increased antibiotic resistance and biofilm formation, modulation of motility, and increased capacity for methionine biosynthesis. Our work demonstrates that a single residue change within the bacterial primary sigma factor can lead to multiple alterations in gene expression and phenotypic changes, suggesting an underrecognized mechanism by which pathobionts and other strain variants with new phenotypes can emerge.

Genome-wide mapping of the Escherichia coli PhoB regulon reveals many transcriptionally inert, intragenic binding sites

Genome-scale analyses have revealed many transcription factor binding sites within, rather than upstream of genes, raising questions as to the function of these binding sites. Here, we use complementary approaches to map the regulon of the Escherichia coli transcription factor PhoB, a response regulator that controls transcription of genes involved in phosphate homeostasis. Strikingly, the majority of PhoB binding sites are located within genes, but these intragenic sites are not associated with detectable transcription regulation and are not evolutionarily conserved. Many intragenic PhoB sites are located in regions bound by H-NS, likely due to shared sequence preferences of PhoB and H-NS. However, these PhoB binding sites are not associated with transcription regulation even in the absence of H-NS. We propose that for many transcription factors, including PhoB, binding sites not associated with promoter sequences are transcriptionally inert, and hence are tolerated as genomic "noise".

IMPORTANCE

Recent studies have revealed large numbers of transcription factor binding sites within the genes of bacteria. The function, if any, of the vast majority of these binding sites has not been investigated. Here, we map the binding of the transcription factor PhoB across the Escherichia coli genome, revealing that the majority of PhoB binding sites are within genes. We show that PhoB binding sites within genes are not associated with regulation of the overlapping genes. Indeed, our data suggest that bacteria tolerate the presence of large numbers of non-regulatory, intragenic binding sites for transcription factors, and that these binding sites are not under selective pressure.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, E_max, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with E_max but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and E_max values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

The stationary phase-specific sRNA FimR2 is a multifunctional regulator of bacterial motility, biofilm formation and virulence

Bacterial pathogens employ a plethora of virulence factors for host invasion, and their use is tightly regulated to maximize infection efficiency and manage resources in a nutrient-limited environment. Here we show that during Escherichia coli stationary phase the 3' UTR-derived small non-coding RNA FimR2 regulates fimbrial and flagellar biosynthesis at the post-transcriptional level, leading to biofilm formation as the dominant mode of survival under conditions of nutrient depletion. FimR2 interacts with the translational regulator CsrA, antagonizing its functions and firmly tightening control over motility and biofilm formation. Generated through RNase E cleavage, FimR2 regulates stationary phase biology by fine-tuning target mRNA levels independently of the chaperones Hfq and ProQ. The Salmonella enterica orthologue of FimR2 induces effector protein secretion by the type III secretion system and stimulates infection, thus linking the sRNA to virulence. This work reveals the importance of bacterial sRNAs in modulating various aspects of bacterial physiology including stationary phase and virulence.

Exopolysaccharides amylovoran and levan contribute to sliding motility in the fire blight pathogen Erwinia amylovora

Erwinia amylovora, the causative agent of fire blight, uses flagella-based motilities to translocate to host plant natural openings; however, little is known about how this bacterium migrates systemically in the apoplast. Here, we reveal a novel surface motility mechanism, defined as sliding, in E. amylovora. Deletion of flagella assembly genes did not affect this movement, whereas deletion of biosynthesis genes for the exopolysaccharides (EPSs) amylovoran and levan resulted in non-sliding phenotypes. Since EPS production generates osmotic pressure that potentially powers sliding, we validated this mechanism by demonstrating that water potential positively contributes to sliding. In addition, no sliding was observed when the water potential of the surface was lower than -0.5 MPa. Sliding is a passive motility mechanism. We further show that the force of gravity plays a critical role in directing E. amylovora sliding on unconfined surfaces but has a negligible effect when cells are sliding in confined microcapillaries, in which EPS-dependent osmotic pressure acts as the main force. Although amylovoran and levan are both required for sliding, we demonstrate that they exhibit different roles in bacterial communities. In summary, our study provides fundamental knowledge for a better understanding of mechanisms that drive bacterial sliding motility.

Coordination of gene expression with cell size enables Escherichia coli to efficiently maintain motility across conditions

To swim and navigate, motile bacteria synthesize a complex motility machinery involving flagella, motors, and a sensory system. A myriad of studies has elucidated the molecular processes involved, but less is known about the coordination of motility expression with cellular physiology: In Escherichia coli, motility genes are strongly up-regulated in nutrient-poor conditions compared to nutrient-replete conditions; yet a quantitative link to cellular motility has not been developed. Here, we systematically investigated gene expression, swimming behavior, cell growth, and available proteomics data across a broad spectrum of exponential growth conditions. Our results suggest that cells up-regulate the expression of motility genes at slow growth to compensate for reduction in cell size, such that the number of flagella per cell is maintained across conditions. The observed four or five flagella per cell is the minimum number needed to keep the majority of cells motile. This simple regulatory objective allows E. coli cells to remain motile across a broad range of growth conditions, while keeping the biosynthetic and energetic demands to establish and drive the motility machinery at the minimum needed. Given the strong reduction in flagella synthesis resulting from cell size increases at fast growth, our findings also provide a different physiological perspective on bacterial cell size control: A larger cell size at fast growth is an efficient strategy to increase the allocation of cellular resources to the synthesis of those proteins required for biomass synthesis and growth, while maintaining processes such as motility that are only needed on a per-cell basis.

Dual Regulatory Role Exerted by Cyclic Dimeric GMP To Control FsnR-Mediated Bacterial Swimming

Bacterial motility has great medical and ecological significance because of its essential role in bacterial survival and pathogenesis. Cyclic dimeric GMP (c-di-GMP), a second messenger in bacteria, is the predominant regulator of flagellar synthesis and motility and possesses turnover mechanisms that have been thoroughly investigated. Therefore, much attention has been focused on identifying the upstream stimulatory signals and downstream modules that respond to altered c-di-GMP levels. Here, we systematically analyzed c-di-GMP cyclases and phosphodiesterases in Stenotrophomonas maltophilia to screen for motility regulators. Of these enzymes, we identified and characterized a new phosphodiesterase named SisP, which was found to facilitate bacterial swimming upon stimulation with ferrous iron. SisP-mediated degradation of c-di-GMP leads to FsnR-dependent transcription of flagellar genes. Remarkably, c-di-GMP controls FsnR via two independent mechanisms: by direct binding and indirectly by modulating its phosphorylation state. In this study, we deciphered a novel "one stone, two birds" regulatory strategy of c-di-GMP and uncovered the signal that stimulates c-di-GMP hydrolysis. Facilitation of bacterial swimming motility by ferrous iron might contribute to the higher risk of bacterial infection in acutely ill patients. IMPORTANCE Stenotrophomonas maltophilia has become a great threat to human health because of the high mortality of infected patients. Swimming motility plays a crucial role in regulating bacterial virulence and adaptation. However, limited progress has been made in cyclic dimeric GMP (c-di-GMP) controlling swimming motility of S. maltophilia. Here, we characterized c-di-GMP turnover enzymes encoded by S. maltophilia and dissected the regulatory details of a phosphodiesterase named SisP. We demonstrated that SisP degrades c-di-GMP to fully activate FsnR through directly releasing FsnR from the FsnR-c-di-GMP complex and indirectly increasing its phosphorylation level. This finding uncovered a quantitative, rather than an on-off, regulatory manner employed by c-di-GMP to regulate activities of its effectors. Identification of the specific activation of SisP by ferrous iron proposes SisP as a putative drug-target for controlling bacterial infection and ferrous iron at the wounds or cuts as a putative factor contributing to the higher risk of bacterial infection.

Regulatory Interplay between RNase III and Antisense RNAs in E. coli: the Case of AsflhD and FlhD, Component of the Master Regulator of Motility

In order to respond to ever-changing environmental cues, bacteria display resilient regulatory mechanisms controlling gene expression. At the post-transcriptional level, this is achieved by a combination of RNA-binding proteins, such as ribonucleases (RNases), and regulatory RNAs, including antisense RNAs (asRNAs). Bound to their complementary mRNA, asRNAs are primary targets for the double-strand-specific endoribonuclease, RNase III. Taking advantage of our own and previously published transcriptomic data sets obtained in strains inactivated for RNase III, we selected several candidate asRNAs and confirmed the existence of RNase III-sensitive asRNAs for crp, ompR, phoP, and flhD genes, all encoding global regulators of gene expression in Escherichia coli. Using FlhD, a component of the master regulator of motility (FlhD₄C₂), as our model, we demonstrate that the asRNA AsflhD, transcribed from the coding sequence of flhD, is involved in the fine-tuning of flhD expression and thus participates in the control of motility.

Clinical Escherichia coli: From Biofilm Formation to New Antibiofilm Strategies

Escherichia coli is one of the species most frequently involved in biofilm-related diseases, being especially important in urinary tract infections, causing relapses or chronic infections. Compared to their planktonic analogues, biofilms confer to the bacteria the capacity to be up to 1000-fold more resistant to antibiotics and to evade the action of the host’s immune system. For this reason, biofilm-related infections are very difficult to treat. To develop new strategies against biofilms, it is important to know the mechanisms involved in their formation. In this review, the different steps of biofilm formation in E. coli, the mechanisms of tolerance to antimicrobials and new compounds and strategies to combat biofilms are discussed.

RegulonDB 11.0: Comprehensive high-throughput datasets on transcriptional regulation in Escherichia coli K-12

Genomics has set the basis for a variety of methodologies that produce high-throughput datasets identifying the different players that define gene regulation, particularly regulation of transcription initiation and operon organization. These datasets are available in public repositories, such as the Gene Expression Omnibus, or ArrayExpress. However, accessing and navigating such a wealth of data is not straightforward. No resource currently exists that offers all available high and low-throughput data on transcriptional regulation in Escherichia coli K-12 to easily use both as whole datasets, or as individual interactions and regulatory elements. RegulonDB (https://regulondb.ccg.unam.mx) began gathering high-throughput dataset collections in 2009, starting with transcription start sites, then adding ChIP-seq and gSELEX in 2012, with up to 99 different experimental high-throughput datasets available in 2019. In this paper we present a radical upgrade to more than 2000 high-throughput datasets, processed to facilitate their comparison, introducing up-to-date collections of transcription termination sites, transcription units, as well as transcription factor binding interactions derived from ChIP-seq, ChIP-exo, gSELEX and DAP-seq experiments, besides expression profiles derived from RNA-seq experiments. For ChIP-seq experiments we offer both the data as presented by the authors, as well as data uniformly processed in-house, enhancing their comparability, as well as the traceability of the methods and reproducibility of the results. Furthermore, we have expanded the tools available for browsing and visualization across and within datasets. We include comparisons against previously existing knowledge in RegulonDB from classic experiments, a nucleotide-resolution genome viewer, and an interface that enables users to browse datasets by querying their metadata. A particular effort was made to automatically extract detailed experimental growth conditions by implementing an assisted curation strategy applying Natural language processing and machine learning. We provide summaries with the total number of interactions found in each experiment, as well as tools to identify common results among different experiments. This is a long-awaited resource to make use of such wealth of knowledge and advance our understanding of the biology of the model bacterium E. coli K-12.

The Rcs System Contributes to the Motility Defects of the Twin-Arginine Translocation System Mutant of Extraintestinal Pathogenic Escherichia coli

Flagellum-mediated bacterial motility is important for bacteria to take up nutrients, adapt to environmental changes, and establish infection. The twin-arginine translocation system (Tat) is an important protein export system, playing a critical role in bacterial physiology and pathogenesis. It has been observed for a long time that the Tat system is critical for bacterial motility. However, the underlying mechanism remains unrevealed. In this study, a comparative transcriptomics analysis was performed with extraintestinal pathogenic Escherichia coli (ExPEC), which identified a considerable number of genes differentially expressed when the Tat system was disrupted. Among them, a large proportion of flagellar biosynthesis genes showed downregulation, indicating that transcription regulation plays an important role in mediating the motility defects. We further identified three Tat substrate proteins, MdoD, AmiA, and AmiC, that were responsible for the nonmotile phenotype. The Rcs system was deleted in the Δtat, the ΔmdoD, and the ΔamiAΔamiC strains, which restored the motility of ΔmdoD and partially restored the motility of Δtat and ΔamiAΔamiC. The flagella were also observed in all of the ΔtatΔrcsDB, ΔmdoDΔrcsDB, and ΔamiAΔamiCΔrcsDB strains, but not in the Δtat, ΔmdoD, and ΔamiAΔamiC strains, by using transmission electron microscopy. Quantitative reverse transcription-PCR data revealed that the regulons of the Rcs system displayed differential expression in the tat mutant, indicating that the Rcs signaling was activated. Our results suggest that the Rcs system plays an important role in mediating the motility defects of the tat mutant of ExPEC. IMPORTANCE The Tat system is an important protein export system critical for bacterial physiology and pathogenesis. It has been observed for a long time that the Tat system is critical for bacterial motility. However, the underlying mechanism remains unrevealed. In this study, we combine transcriptomics analysis and bacterial genetics, which reveal that transcription regulation plays an important role in mediating the motility defects of the tat mutant of extraintestinal pathogenic Escherichia coli. The Tat substrate proteins responsible for the motility defects are identified. We further show that the Rcs system contributes to the motility suppression. We for the first time reveal the link between the Tat system and bacterial motility, which is important for understanding the physiological functions of the Tat system.

Multiple in vivo roles for the C-terminal domain of the RNA chaperone Hfq

Hfq, a bacterial RNA chaperone, stabilizes small regulatory RNAs (sRNAs) and facilitates sRNA base-pairing with target mRNAs. Hfq has a conserved N-terminal domain and a poorly conserved disordered C-terminal domain (CTD). In a transcriptome-wide examination of the effects of a chromosomal CTD deletion (Hfq1-65), the Escherichia coli mutant was most defective for the accumulation of sRNAs that bind the proximal and distal faces of Hfq (Class II sRNAs), but other sRNAs also were affected. There were only modest effects on the levels of mRNAs, suggesting little disruption of sRNA-dependent regulation. However, cells expressing Hfq lacking the CTD in combination with a weak distal face mutation were defective for the function of the Class II sRNA ChiX and repression of mutS, both dependent upon distal face RNA binding. Loss of the region between amino acids 66-72 was critical for this defect. The CTD region beyond amino acid 72 was not necessary for distal face-dependent regulation, but was needed for functions associated with the Hfq rim, seen most clearly in combination with a rim mutant. Our results suggest that the C-terminus collaborates in various ways with different binding faces of Hfq, leading to distinct outcomes for individual sRNAs.

Does Silver in Different Forms Affect Bacterial Susceptibility and Resistance? A Mechanistic Perspective

The exceptional increase in antibiotic resistance in past decades motivated the scientific community to use silver as a potential antibacterial agent. However, due to its unknown antibacterial mechanism and the pattern of bacterial resistance to silver species, it has not been revolutionized in the health sector. This study deciphers mechanistic aspects of silver species, i.e., ions and lysozyme-coated silver nanoparticles (L-Ag NPs), against E. coli K12 through RNA sequencing analysis. The obtained results support the reservoir nature of nanoparticles for the controlled release of silver ions into bacteria. This study differentiates between the antibacterial mechanism of silver species by discussing the pathway of their entry in bacteria, sequence of events inside cells, and response of bacteria to overcome silver stress. Controlled release of ions from L-Ag NPs not only reduces bacterial growth but also reduces the likelihood of resistance development. Conversely, direct exposure of silver ions, leads to rapid activation of the bacterial defense system leading to development of resistance against silver ions, like the well-known antibiotic resistance problem. These findings provide valuable insight on the mechanism of silver resistance and antibacterial strategies deployed by E. coli K12, which could be a potential target for the generation of aim-based and effective nanoantibiotics.

Disentangling direct from indirect relationships in association networks

Networks are vital tools for understanding and modeling interactions in complex systems in science and engineering, and direct and indirect interactions are pervasive in all types of networks. However, quantitatively disentangling direct and indirect relationships in networks remains a formidable task. Here, we present a framework, called iDIRECT (Inference of Direct and Indirect Relationships with Effective Copula-based Transitivity), for quantitatively inferring direct dependencies in association networks. Using copula-based transitivity, iDIRECT eliminates/ameliorates several challenging mathematical problems, including ill-conditioning, self-looping, and interaction strength overflow. With simulation data as benchmark examples, iDIRECT showed high prediction accuracies. Application of iDIRECT to reconstruct gene regulatory networks in Escherichia coli also revealed considerably higher prediction power than the best-performing approaches in the DREAM5 (Dialogue on Reverse Engineering Assessment and Methods project, #5) Network Inference Challenge. In addition, applying iDIRECT to highly diverse grassland soil microbial communities in response to climate warming showed that the iDIRECT-processed networks were significantly different from the original networks, with considerably fewer nodes, links, and connectivity, but higher relative modularity. Further analysis revealed that the iDIRECT-processed network was more complex under warming than the control and more robust to both random and target species removal (P < 0.001). As a general approach, iDIRECT has great advantages for network inference, and it should be widely applicable to infer direct relationships in association networks across diverse disciplines in science and engineering.

Temperature-Dependent Influence of FliA Overexpression on PHL628 E. coli Biofilm Growth and Composition

Biofilm growth and survival pose a problem in both medical and industrial fields. Bacteria in biofilms are more tolerant to antibiotic treatment due to the inability of antibiotics to permeate to the bottom layers of cells in a biofilm and the creation of altered microenvironments of bacteria deep within the biofilm. Despite the abundance of information we have about E. coli biofilm growth and maturation, we are still learning how manipulating different signaling pathways influences the formation and fitness of biofilm. Understanding the impact of signaling pathways on biofilm formation may narrow the search for novel small molecule inhibitors or activators that affect biofilm production and stability. Here, we study the influence of the minor sigma transcription factor FliA (RpoF, sigma-28), which controls late-stage flagellar assembly and chemotaxis, on biofilm production and composition at various temperatures in the E. coli strain PHL628, which abundantly produces the extracellular structural protein curli. We examined FliA's influence on external cellular structures like curli and flagella and the biomolecular composition of the biofilm's extracellular polymeric substance (EPS) using biochemical assays, immunoblotting, and confocal laser scanning microscopy (CLSM). At 37°C, FliA overexpression results in the dramatic growth of biofilm in polystyrene plates and more modest yet significant biofilm growth on silica slides. We observed no significant differences in curli concentration and carbohydrate concentration in the EPS with FliA overexpression. Still, we did see significant changes in the abundance of EPS protein using CLSM at higher growth temperatures. We also noticed increased flagellin concentration, a major structural protein in flagella, occurred with FliA overexpression, specifically in planktonic cultures. These experiments have aided in narrowing our focus to FliA's role in changing the protein composition of the EPS, which we will examine in future endeavors.

Bacterial RNA chaperones and chaperone-like riboregulators: behind the scenes of RNA-mediated regulation of cellular metabolism

In all domains of life, RNA chaperones safeguard and guide the fate of the cellular RNA pool. RNA chaperones comprise structurally diverse proteins that ensure proper folding, stability, and ribonuclease resistance of RNA, and they support regulatory activities mediated by RNA. RNA chaperones constitute a topologically diverse group of proteins that often present an unstructured region and bind RNA with limited nucleotide sequence preferences. In bacteria, three main proteins - Hfq, ProQ, and CsrA - have been shown to regulate numerous complex processes, including bacterial growth, stress response and virulence. Hfq and ProQ have well-studied activities as global chaperones with pleiotropic impact, while CsrA has a chaperone-like role with more defined riboregulatory function. Here, we describe relevant novel insights into their common features, including RNA binding properties, unstructured domains, and interplay with other proteins important to RNA metabolism.

Gain of Spontaneous clpX Mutations Boosting Motility via Adaption to Environments in Escherichia coli

Motility is finely regulated and is crucial to bacterial processes including colonization and biofilm formation. There is a trade-off between motility and growth in bacteria with molecular mechanisms not fully understood. Hypermotile Escherichia coli could be isolated by evolving non-motile cells on soft agar plates. Most of the isolates carried mutations located upstream of the flhDC promoter region, which upregulate the transcriptional expression of the master regulator of the flagellum biosynthesis, FlhDC. Here, we identified that spontaneous mutations in clpX boosted the motility of E. coli largely, inducing several folds of changes in swimming speed. Among the mutations identified, we further elucidated the molecular mechanism underlying the ClpX^V78F mutation on the regulation of E. coli motility. We found that the V78F mutation affected ATP binding to ClpX, resulting in the inability of the mutated ClpXP protease to degrade FlhD as indicated by both structure modeling and in vitro protein degradation assays. Moreover, our proteomic data indicated that the ClpX^V78F mutation elevated the stability of known ClpXP targets to various degrees with FlhD as one of the most affected. In addition, the specific tag at the C-terminus of FlhD being recognized for ClpXP degradation was identified. Finally, our transcriptome data characterized that the enhanced expression of the motility genes in the ClpX^V78F mutations was intrinsically accompanied by the reduced expression of stress resistance genes relating to the reduced fitness of the hypermotile strains. A similar pattern was observed for previously isolated hypermotile E. coli strains showing high expression of flhDC at the transcriptional level. Hence, clpX appears to be a hot locus comparable to the upstream of the flhDC promoter region evolved to boost bacterial motility, and our finding provides insight into the reduced fitness of the hypermotile bacteria.

Efficient markerless integration of genes in the chromosome of probiotic E. coli Nissle 1917 by bacterial conjugation

The probiotic strain Escherichia coli Nissle 1917 (EcN) is a common bacterial chassis in synthetic biology developments for therapeutic applications given its long track record of safe administration in humans. Chromosomal integration of the genes of interest (GOIs) in the engineered bacterium offers significant advantages in genetic stability and to control gene dose, but common methodologies relying on the transformation of EcN are inefficient. In this work, we implement in EcN the use of bacterial conjugation in combination with markerless genome engineering to efficiently insert multiple GOIs at different loci of EcN chromosome, leaving no antibiotic resistance genes, vector sequences or scars in the modified bacterium. The resolution of cointegrants that leads to markerless insertion of the GOIs requires expression of I-SceI endonuclease and its efficiency is enhanced by λ Red proteins. We show the potential of this strategy by integrating different genes encoding fluorescent and bioluminescent reporters (i.e. GFP, mKate2, luxCDABE) both individually and sequentially. We also demonstrate its application for gene deletions in EcN (ΔflhDC) and to replace the endogenous regulation of chromosomal locus (i.e. flhDC) by heterologous regulatory elements (e.g. tetR-Ptet) in order to have an ectopic control of gene expression in EcN with an external inducer to alter bacterial behaviour (e.g. flagellar motility). Whole-genome sequencing confirmed the introduction of the designed modifications without off-target alterations in the genome. This straightforward approach accelerates the generation of multiple modifications in EcN chromosome for the generation of living bacterial therapeutics.

Multiple Copies of flhDC in Paraburkholderia unamae Regulate Flagellar Gene Expression, Motility, and Biofilm Formation

FlhDC is a heterohexameric complex that acts as a master regulator of flagellar biosynthesis genes in numerous bacteria. Previous studies have identified a single flhDC operon encoding this complex. However, we found that two flhDC loci are present throughout Paraburkholderia, and two additional flhC copies are also present in Paraburkholderia unamae. Systematic deletion analysis in P. unamae of the different flhDC copies showed that one of the operons, flhDC1, plays the predominant role, with deletion of its genes resulting in a severe inhibition of motility and biofilm formation. Expression analysis using promoter-lacZ fusions and real-time quantitative PCR support the primary role of flhDC1 in flagellar gene regulation, with flhDC2 a secondary contributor. Phylogenetic analysis shows the presence of the flhDC1 and flhDC2 operons throughout Paraburkholderia. In contrast, Burkholderia and other bacteria only carry the copy syntenous with flhDC2. The variations in impact each copy of flhDC has on downstream processes indicate that regulation of FlhDC in P. unamae, and likely other Paraburkholderia species, is regulated at least in part by the presence of multiple copies of these genes. IMPORTANCE Motility is important in the colonization of plant roots by beneficial and pathogenic bacteria, with flagella playing essential roles in host cell adhesion, entrance, and biofilm formation. Flagellar biosynthesis is energetically expensive. Its complex regulation by the FlhDC master regulator is well studied in peritrichous flagella expressing enterics. We report the unique presence throughout Paraburkholderia of multiple copies of flhDC. In P. unamae, the flhDC1 copy showed higher expression and a greater effect on swim motility, flagellar development, and regulation of downstream genes, than the flhDC2 copy that is syntenous to flhDC in Escherichia coli and pathogenic Burkholderia spp. The flhDC genes have evolved differently in these plant-growth-promoting bacteria, giving an additional layer of complexity in gene regulation by FlhDC.

OsnR is an autoregulatory negative transcription factor controlling redox-dependent stress responses in Corynebacterium glutamicum

Background

Corynebacterium glutamicum is used in the industrial production of amino acids and nucleotides. During the course of fermentation, C. glutamicum cells face various stresses and employ multiple regulatory genes to cope with the oxidative stress. The osnR gene plays a negative regulatory role in redox-dependent oxidative-stress responses, but the underlying mechanism is not known yet.

Results

Overexpression of the osnR gene in C. glutamicum affected the expression of genes involved in the mycothiol metabolism. ChIP-seq analysis revealed that OsnR binds to the promoter region of multiple genes, including osnR and cg0026, which seems to function in the membrane-associated redox metabolism. Studies on the role of the osnR gene involving in vitro assays employing purified OsnR proteins and in vivo physiological analyses have identified that OsnR inhibits the transcription of its own gene. Further, oxidant diamide stimulates OsnR-binding to the promoter region of the osnR gene. The genes affected by the overexpression of osnR have been found to be under the control of σ^H. In the osnR-overexpressing strain, the transcription of sigH is significantly decreased and the stimulation of sigH transcription by external stress is lost, suggesting that osnR and sigH form an intimate regulatory network.

Conclusions

Our study suggests that OsnR not only functions as a transcriptional repressor of its own gene and of those involved in redox-dependent stress responses but also participates in the global transcriptional regulation by controlling the transcription of other master regulators, such as sigH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01693-1.

Saturation mutagenesis charts the functional landscape of Salmonella ProQ and reveals a gene regulatory function of its C-terminal domain

The global RNA-binding protein ProQ has emerged as a central player in post-transcriptional regulatory networks in bacteria. While the N-terminal domain (NTD) of ProQ harbors the major RNA-binding activity, the role of the ProQ C-terminal domain (CTD) has remained unclear. Here, we have applied saturation mutagenesis coupled to phenotypic sorting and long-read sequencing to chart the regulatory capacity of Salmonella ProQ. Parallel monitoring of thousands of ProQ mutants allowed mapping of critical residues in both the NTD and the CTD, while the linker separating these domains was tolerant to mutations. Single amino acid substitutions in the NTD associated with abolished regulatory capacity strongly align with RNA-binding deficiency. An observed cellular instability of ProQ associated with mutations in the NTD suggests that interaction with RNA protects ProQ from degradation. Mutation of conserved CTD residues led to overstabilization of RNA targets and rendered ProQ inert in regulation, without affecting protein stability in vivo. Furthermore, ProQ lacking the CTD, although binding competent, failed to protect an mRNA target from degradation. Together, our data provide a comprehensive overview of residues important for ProQ-dependent regulation and reveal an essential role for the enigmatic ProQ CTD in gene regulation.

Dissecting psa Locus Regulation in Yersinia pestis

The pH 6 antigen (PsaA) of Yersinia pestis is a virulence factor that is expressed in response to high temperature (37°C) and low pH (6.0). Previous studies have implicated the PsaE and PsaF regulators in the temperature- and pH-dependent regulation of psaA. Here, we show that PsaE levels are themselves controlled by pH and temperature, explaining the regulation of psaA. We identify hundreds of binding sites for PsaE across the Y. pestis genome, with the majority of binding sites located in intergenic regions bound by the nucleoid-associated protein H-NS. However, we detect direct regulation of only two transcripts by PsaE, likely due to displacement of H-NS from the corresponding promoter regions; our data suggest that most PsaE binding sites are nonregulatory or that they require additional environmental cues. We also identify the precise binding sites for PsaE that are required for temperature- and pH-dependent regulation of psaA and psaE. Thus, our data reveal the critical role that PsaE plays in the regulation of psaA and suggest that PsaE may have many additional regulatory targets. IMPORTANCE Y. pestis, the etiologic agent of plague, has been responsible for high mortality in several epidemics throughout human history. The plague bacillus has been used as a biological weapon during human history and is currently one of the most likely biological threats. PsaA and PsaE appear to play important roles during Y. pestis infection. Understanding their regulation by environmental cues would facilitate a solution to impede Y. pestis infection.

Crucial role of Salmonella genomic island 1 master activator in the parasitism of IncC plasmids

IncC conjugative plasmids and the multiple variants of Salmonella Genomic Island 1 (SGI1) are two functionally interacting families of mobile genetic elements commonly associated with multidrug resistance in the Gammaproteobacteria. SGI1 and its siblings are specifically mobilised in trans by IncC conjugative plasmids. Conjugative transfer of IncC plasmids is activated by the plasmid-encoded master activator AcaCD. SGI1 carries five AcaCD-responsive promoters that drive the expression of genes involved in its excision, replication, and mobilisation. SGI1 encodes an AcaCD homologue, the transcriptional activator complex SgaCD (also known as FlhDCSGI1) that seems to recognise and activate the same SGI1 promoters. Here, we investigated the relevance of SgaCD in SGI1's lifecycle. Mating assays revealed the requirement for SgaCD and its IncC-encoded counterpart AcaCD in the mobilisation of SGI1. An integrative approach combining ChIP-exo, Cappable-seq, and RNA-seq confirmed that SgaCD activates each of the 18 AcaCD-responsive promoters driving the expression of the plasmid transfer functions. A comprehensive analysis of the activity of the complete set of AcaCD-responsive promoters of SGI1 and the helper IncC plasmid was performed through reporter assays. qPCR and flow cytometry assays revealed that SgaCD is essential to elicit the excision and replication of SGI1 and destabilise the helper IncC plasmid.

Diversity in the swimming motility and flagellar regulon structure of uropathogenic Morganella morganii strains

In current times, the opportunistic pathogen Morganella morganii is increasingly becoming a cause of urinary tract infections. The condition has been further complicated by the multiple drug resistance of most isolates. Swimming motility plays an important role in the development of urinary tract infections, allowing bacteria to colonize the upper urinary tract. We determined the differences between the growth, swimming motility, and biofilm formation of two M. morganii strains MM 1 and MM 190 isolated from the urine of patients who had community-acquired urinary tract infections. MM 190 showed a lower growth rate but better-formed biofilms in comparison to MM 1. In addition, MM 190 possessed autoaggregation abilities. It was found that a high temperature (37 °C) inhibits the flagellation of strains and makes MM 190 less motile. At the same time, the MM 1 strain maintained its rate of motility at this temperature. We demonstrated that urea at a concentration of 1.5% suppresses the growth and swimming motility of both strains. Genome analysis showed that MM 1 has a 17.7-kb-long insertion in flagellar regulon between fliE and glycosyltransferase genes, which was not identified in corresponding loci of MM 190 and 9 other M. morganii strains with whole genomes. Both strains carry two genes encoding flagellin, which may indicate flagellar antigen phase variation. However, the fliC2 genes have only 91% identity to each other and exhibit some variability in the regulatory region. We assume that all these differences influence the swimming motility of the strains.

Two homologous Salmonella serogroup C1-specific genes are required for flagellar motility and cell invasion

Background

Salmonella is a major bacterial pathogen associated with a large number of outbreaks of foodborne diseases. Many highly virulent serovars that cause human illness belong to Salmonella serogroup C1, and Salmonella ser. Choleraesuis is a prominent cause of invasive infections in Asia. Comparative genomic analysis in our previous study showed that two homologous genes, SC0368 and SC0595 in Salmonella ser. Choleraesuis were unique to serogroup C1. In this study, two single-deletion mutants (Δ0368 and Δ0595) and one double-deletion mutant (Δ0368Δ0595) were constructed based on the genome. All these mutants and the wild-type strain were subjected to RNA-Seq analysis to reveal functional relationships of the two serogroup C1-specific genes.

Results

Data from RNA-Seq indicated that deletion of SC0368 resulted in defects in motility through repression of σ²⁸ in flagellar regulation Class 3. Consistent with RNA-Seq data, results from transmission electron microcopy (TEM) showed that flagella were not present in △0368 and △0368△0595 mutants resulting in both swimming and swarming defects. Interestingly, the growth rates of two non-motile mutants △0368 and △0368△0595 were significantly greater than the wild-type, which may be associated with up-regulation of genes encoding cytochromes, enhancing bacterial proliferation. Moreover, the △0595 mutant was significantly more invasive in Caco-2 cells as shown by bacterial enumeration assays, and the expression of lipopolysaccharide (LPS) core synthesis-related genes (rfaB, rfaI, rfaQ, rfaY, rfaK, rfaZ) was down-regulated only in the △0368△0595 mutant. In addition, this study also speculated that these two genes might be contributing to serotype conversion for Salmonella C1 serogroup based on their apparent roles in biosynthesis of LPS and the flagella.

Conclusion

A combination of biological and transcriptomic (RNA-Seq) analyses has shown that the SC0368 and SC0595 genes are involved in biosynthesis of flagella and complete LPS, as well as in bacterial growth and virulence. Such information will aid to revealing the role of these specific genes in bacterial physiology and evolution within the serogroup C1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07759-z.

Genomic Analysis of the Endophytic Stenotrophomonas Strain 169 Reveals Features Related to Plant-Growth Promotion and Stress Tolerance

Plant-associated Stenotrophomonas isolates have great potential for plant growth promotion, especially under stress conditions, due to their ability to promote tolerance to abiotic stresses such as salinity or drought. The endophytic strain Stenotrophomonas sp. 169, isolated from a field-grown poplar, increased the growth of inoculated in vitro plants, with a particular effect on root development, and was able to stimulate the rooting of poplar cuttings in the greenhouse. The strain produced high amounts of the plant growth-stimulating hormone auxin under in vitro conditions. The comparison of the 16S rRNA gene sequences and the phylogenetic analysis of the core genomes showed a close relationship to Stenotrophomonas chelatiphaga and a clear separation from Stenotrophomonas maltophilia. Whole genome sequence analysis revealed functional genes potentially associated with attachment and plant colonization, growth promotion, and stress protection. In detail, an extensive set of genes for twitching motility, chemotaxis, flagella biosynthesis, and the ability to form biofilms, which are connected with host plant colonization, could be identified in the genome of strain 169. The production of indole-3-acetic acid and the presence of genes for auxin biosynthesis pathways and the spermidine pathway could explain the ability to promote plant growth. Furthermore, the genome contained genes encoding for features related to the production of different osmoprotective molecules and enzymes mediating the regulation of stress tolerance and the ability of bacteria to quickly adapt to changing environments. Overall, the results of physiological tests and genome analysis demonstrated the capability of endophytic strain 169 to promote plant growth. In contrast to related species, strain 169 can be considered non-pathogenic and suitable for biotechnology applications.

Roles of OmpX, an Outer Membrane Protein, on Virulence and Flagellar Expression in Uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) is a major pathogen that causes urinary tract infection (UTI). This bacterium adheres to and internalizes within urinary tract cells, where it aggregates and subsequently forms biofilm-like multicellular colonies that protect UPEC from antimicrobial agents and the host's immune system. Here, we show that OmpX, an outer membrane protein, plays a role in the pathogenesis of UPEC in renal cells. Deletion of ompX decreased bacterial internalization and aggregation within kidney epithelial cells and also impaired the colonization of mouse urinary tracts, but the ompX mutant still adhered to the epithelial cells at a level similar to that of the parent strain. FlhD, the master regulator of flagellum-related genes, had a low expression level in the ompX mutant compared to the parent strain, and the ompX mutant exhibited defective motility due to lower flagellar production than the parent strain. The fliC mutant, which lacks flagella, exhibited lower levels of bacterial internalization and aggregation than the parent strain. Additional deletion of ompX in the fliC mutant did not further decrease bacterial internalization. These combined results suggest that OmpX contributes to flagellar production in UPEC and then sustains UPEC virulence associated with bacterial internalization and aggregation within urinary tract cells and colonization in the urinary tract.