Predicting the location and strength of promoters from genomic sequence requires accurate sequenced-based promoter models. We present the first model of a full-length bacterial promoter, encompassing both upstream sequences (UP-elements) and core promoter modules, based on a set of 60 promoters dependent on σ(E), an alternative ECF-type σ factor. UP-element contribution, best described by the length and frequency of A- and T-tracts, in combination with a PWM-based core promoter model, accurately predicted promoter strength both in vivo and in vitro. This model also distinguished active from weak/inactive promoters. Systematic examination of promoter strength as a function of RNA polymerase (RNAP) concentration revealed that UP-element contribution varied with RNAP availability and that the σ(E) regulon is comprised of two promoter types, one of which is active only at high concentrations of RNAP. Distinct promoter types may be a general mechanism for increasing the regulatory capacity of the ECF group of alternative σ's. Our findings provide important insights into the sequence requirements for the strength and function of full-length promoters and establish guidelines for promoter prediction and for forward engineering promoters of specific strengths.

Bacteria often cope with environmental stress by inducing alternative sigma (sigma) factors, which direct RNA polymerase to specific promoters, thereby inducing a set of genes called a regulon to combat the stress. To understand the conserved and organism-specific functions of each sigma, it is necessary to be able to predict their promoters, so that their regulons can be followed across species. However, the variability of promoter sequences and motif spacing makes their prediction difficult. We developed and validated an accurate promoter prediction model for Escherichia coli sigmaE, which enabled us to predict a total of 89 unique sigmaE-controlled transcription units in E. coli K-12 and eight related genomes. SigmaE controls the envelope stress response in E. coli K-12. The portion of the regulon conserved across genomes is functionally coherent, ensuring the synthesis, assembly, and homeostasis of lipopolysaccharide and outer membrane porins, the key constituents of the outer membrane of Gram-negative bacteria. The larger variable portion is predicted to perform pathogenesis-associated functions, suggesting that sigmaE provides organism-specific functions necessary for optimal host interaction. The success of our promoter prediction model for sigmaE suggests that it will be applicable for the prediction of promoter elements for many alternative sigma factors.

A systems-level understanding of Gram-positive bacteria is important from both an environmental and health perspective and is most easily obtained when high-quality, validated genomic resources are available. To this end, we constructed two ordered, barcoded, erythromycin-resistance- and kanamycin-resistance-marked single-gene deletion libraries of the Gram-positive model organism, Bacillus subtilis. The libraries comprise 3,968 and 3,970 genes, respectively, and overlap in all but four genes. Using these libraries, we update the set of essential genes known for this organism, provide a comprehensive compendium of B. subtilis auxotrophic genes, and identify genes required for utilizing specific carbon and nitrogen sources, as well as those required for growth at low temperature. We report the identification of enzymes catalyzing several missing steps in amino acid biosynthesis. Finally, we describe a suite of high-throughput phenotyping methodologies and apply them to provide a genome-wide analysis of competence and sporulation. Altogether, we provide versatile resources for studying gene function and pathway and network architecture in Gram-positive bacteria.

DNA topoisomerases play essential roles in chromosome organization and replication. Most bacteria possess multiple topoisomerases which have specialized functions in the control of DNA supercoiling or in DNA catenation/decatenation during recombination and chromosome segregation. DNA topoisomerase I is required for the relaxation of negatively supercoiled DNA behind the transcribing RNA polymerase. Conflicting results have been reported on the essentiality of the topA gene encoding topoisomerase I in the model bacterium Bacillus subtilis. In this work, we have studied the requirement for topoisomerase I in B. subtilis. All stable topA mutants carried different chromosomal amplifications of the genomic region encompassing the parEC operon encoding topoisomerase IV. Using a fluorescent amplification reporter system we observed that each individual topA mutant had acquired such an amplification. Eventually, the amplifications were replaced by a point mutation in the parEC promoter region which resulted in a fivefold increase of parEC expression. In this strain both type I topoisomerases, encoded by topA and topB, were dispensable. Our results demonstrate that topoisomerase IV at increased expression is necessary and sufficient to take over the function of type 1A topoisomerases.

Extracytoplasmic function σ factors (ECFs) represent one of the major bacterial signal transduction mechanisms in terms of abundance, diversity and importance, particularly in mediating stress responses. Here, we performed a comprehensive phylogenetic analysis of this protein family by scrutinizing all proteins in the NCBI database. As a result, we identified an average of ∼10 ECFs per bacterial genome and 157 phylogenetic ECF groups that feature a conserved genetic neighborhood and a similar regulation mechanism. Our analysis expands previous classification efforts ∼50-fold, enriches many original ECF groups with previously unclassified proteins and identifies 22 entirely new ECF groups. The ECF groups are hierarchically related to each other and are further composed of subgroups with closely related sequences. This two-tiered classification allows for the accurate prediction of common promoter motifs and the inference of putative regulatory mechanisms across subgroups composing an ECF group. This comprehensive, high-resolution description of the phylogenetic distribution of the ECF family, together with the massive expansion of classified ECF sequences and an openly accessible data repository called 'ECF Hub' (https://www.computational.bio.uni-giessen.de/ecfhub), will serve as a powerful hypothesis-generator to guide future research in the field.

Essential genes are the hubs of cellular networks, but lack of high-throughput methods for titrating gene expression has limited our understanding of the fitness landscapes against which their expression levels are optimized. We developed a modified CRISPRi system leveraging the predictable reduction in efficacy of imperfectly matched sgRNAs to generate defined levels of CRISPRi activity and demonstrated its broad applicability. Using libraries of mismatched sgRNAs predicted to span the full range of knockdown levels, we characterized the expression-fitness relationships of most essential genes in Escherichia coli and Bacillus subtilis. We find that these relationships vary widely from linear to bimodal but are similar within pathways. Notably, despite ∼2 billion years of evolutionary separation between E. coli and B. subtilis, most essential homologs have similar expression-fitness relationships with rare but informative differences. Thus, the expression levels of essential genes may reflect homeostatic or evolutionary constraints shared between the two organisms.

Kasugamycin (KSG) is an aminoglycoside antibiotic widely used in agriculture and exhibits considerable medical potential. Previous studies suggested that KSG interferes with translation by blocking binding of canonical messenger RNA (mRNA) and initiator transfer tRNA (tRNA) to the small ribosomal subunit, thereby preventing initiation of protein synthesis. Here, by using genome-wide approaches, we show that KSG can interfere with translation even after the formation of the 70S initiation complex on mRNA, as the extent of KSG-mediated translation inhibition correlates with increased occupancy of start codons by 70S ribosomes. Even at saturating concentrations, KSG does not completely abolish translation, allowing for continuing expression of some Escherichia coli proteins. Differential action of KSG significantly depends on the nature of the mRNA residue immediately preceding the start codon, with guanine in this position being the most conducive to inhibition by the drug. In addition, the activity of KSG is attenuated by translational coupling as genes whose start codons overlap with the coding regions or the stop codons of the upstream cistrons tend to be less susceptible to drug-mediated inhibition. Altogether, our findings reveal KSG as an example of a small ribosomal subunit-targeting antibiotic with a well-pronounced context specificity of action.

Ceragenins are a family of synthetic amphipathic molecules designed to mimic the properties of naturally occurring cationic antimicrobial peptides (CAMPs). Although ceragenins have potent antimicrobial activity, whether their mode of action is similar to that of CAMPs has remained elusive. Here, we reported the results of a comparative study of the bacterial responses to two well-studied CAMPs, LL37 and colistin, and two ceragenins with related structures, CSA13 and CSA131. Using transcriptomic and proteomic analyses, we found that Escherichia coli responded similarly to both CAMPs and ceragenins by inducing a Cpx envelope stress response. However, whereas E. coli exposed to CAMPs increased expression of genes involved in colanic acid biosynthesis, bacteria exposed to ceragenins specifically modulated functions related to phosphate transport, indicating distinct mechanisms of action between these two classes of molecules. Although traditional genetic approaches failed to identify genes that confer high-level resistance to ceragenins, using a Clustered Regularly Interspaced Short Palindromic Repeats interference (CRISPRi) approach we identified E. coli essential genes that when knocked down modify sensitivity to these molecules. Comparison of the essential gene-antibiotic interactions for each of the CAMPs and ceragenins identified both overlapping and distinct dependencies for their antimicrobial activities. Overall, this study indicated that, while some bacterial responses to ceragenins overlap those induced by naturally occurring CAMPs, these synthetic molecules target the bacterial envelope using a distinctive mode of action. IMPORTANCE The development of novel antibiotics is essential because the current arsenal of antimicrobials will soon be ineffective due to the widespread occurrence of antibiotic resistance. The development of naturally occurring cationic antimicrobial peptides (CAMPs) for therapeutics to combat antibiotic resistance has been hampered by high production costs and protease sensitivity, among other factors. The ceragenins are a family of synthetic CAMP mimics that kill a broad spectrum of bacterial species but are less expensive to produce, resistant to proteolytic degradation, and seemingly resistant to the development of high-level resistance. Determining how ceragenins function may identify new essential biological pathways of bacteria that are less prone to the development of resistance and will further our understanding of the design principles for maximizing the effects of synthetic CAMPs.

All cells must adapt to rapidly changing conditions. The heat shock response (HSR) is an intracellular signaling pathway that maintains proteostasis (protein folding homeostasis), a process critical for survival in all organisms exposed to heat stress or other conditions that alter the folding of the proteome. Yet despite decades of study, the circuitry described for responding to altered protein status in the best-studied bacterium, E. coli, does not faithfully recapitulate the range of cellular responses in response to this stress. Here, we report the discovery of the missing link. Surprisingly, we found that σ(32), the central transcription factor driving the HSR, must be localized to the membrane rather than dispersed in the cytoplasm as previously assumed. Genetic analyses indicate that σ(32) localization results from a protein targeting reaction facilitated by the signal recognition particle (SRP) and its receptor (SR), which together comprise a conserved protein targeting machine and mediate the cotranslational targeting of inner membrane proteins to the membrane. SRP interacts with σ(32) directly and transports it to the inner membrane. Our results show that σ(32) must be membrane-associated to be properly regulated in response to the protein folding status in the cell, explaining how the HSR integrates information from both the cytoplasm and bacterial cell membrane.

The Gram-positive bacterium Bacillus subtilis uses serine not only as a building block for proteins but also as an important precursor in many anabolic reactions. Moreover, a lack of serine results in the initiation of biofilm formation. However, excess serine inhibits the growth of B. subtilis. To unravel the underlying mechanisms, we isolated suppressor mutants that can tolerate toxic serine concentrations by three targeted and non-targeted genome-wide screens. All screens as well as genetic complementation in Escherichia coli identified the so far uncharacterized permease YbeC as the major serine transporter of B. subtilis. In addition to YbeC, the threonine transporters BcaP and YbxG make minor contributions to serine uptake. A strain lacking these three transporters was able to tolerate 100 mM serine whereas the wild type strain was already inhibited by 1 mM of the amino acid. The screen for serine-resistant mutants also identified mutations that result in increased serine degradation and in increased expression of threonine biosynthetic enzymes suggesting that serine toxicity results from interference with threonine biosynthesis.

CRISPR interference is an increasingly popular method for perturbing gene expression. Guided by single-guide RNAs (sgRNAs), nuclease-deficient Cas9 proteins bind to specific DNA sequences and hinder transcription. Specificity is achieved through complementarity of the sgRNAs to the DNA. Changing complementarity by introducing single-nucleotide mismatches can be exploited to tune knockdown. Here, we present a computational pipeline to identify sgRNAs targeting specific genes in a bacterial genome, filter them, and titrate their activity by introducing mismatches. For complete details on the use and execution of this protocol, please refer to Hawkins et al. (2020).

To maintain cellular structure and integrity during division, Gram-negative bacteria must carefully coordinate constriction of a tripartite cell envelope of inner membrane, peptidoglycan (PG), and outer membrane (OM). It has remained enigmatic how this is accomplished. Here, we show that envelope machines facilitating septal PG synthesis (PBP1B-LpoB complex) and OM constriction (Tol system) are physically and functionally coordinated via YbgF, renamed CpoB (Coordinator of PG synthesis and OM constriction, associated with PBP1B). CpoB localizes to the septum concurrent with PBP1B-LpoB and Tol at the onset of constriction, interacts with both complexes, and regulates PBP1B activity in response to Tol energy state. This coordination links PG synthesis with OM invagination and imparts a unique mode of bifunctional PG synthase regulation by selectively modulating PBP1B cross-linking activity. Coordination of the PBP1B and Tol machines by CpoB contributes to effective PBP1B function in vivo and maintenance of cell envelope integrity during division.

The sigma factors are the key regulators of bacterial transcription. ECF (extracytoplasmic function) sigma's are the largest and most divergent group of sigma(70) family members. ECF sigma's are normally sequestered in an inactive complex by their specific anti-sigma factor, which often spans the inner membrane. Here, we determined the 2 A resolution crystal structure of the Escherichia coli ECF sigma factor sigma(E) in an inhibitory complex with the cytoplasmic domain of its anti-sigma, RseA. Despite extensive sequence variability, the two major domains of sigma(E) are virtually identical in structure to the corresponding domains of other sigma(70) family members. In combination with a model of the sigma(E) holoenzyme and biochemical data, the structure reveals that RseA functions by sterically occluding the two primary binding determinants on sigma(E) for core RNA polymerase.

Temperature influences the structural and functional properties of cellular components, necessitating stress responses to restore homeostasis following temperature shift. Whereas the circuitry controlling the heat shock response is well understood, that controlling the E. coli cold shock adaptation program is not. We found that during the growth arrest phase (acclimation) that follows shift to low temperature, protein synthesis increases, and open reading frame (ORF)-wide mRNA secondary structure decreases. To identify the regulatory system controlling this process, we screened for players required for increased translation. We identified a two-member mRNA surveillance system that enables recovery of translation during acclimation: RNase R assures appropriate mRNA degradation and the Csps dynamically adjust mRNA secondary structure to globally modulate protein expression level. An autoregulatory switch in which Csps tune their own expression to cellular demand enables dynamic control of global translation. The universality of Csps in bacteria suggests broad utilization of this control mechanism.

Bacterial mRNAs are organized into operons consisting of discrete open reading frames (ORFs) in a single polycistronic mRNA. Individual ORFs on the mRNA are differentially translated, with rates varying as much as 100-fold. The signals controlling differential translation are poorly understood. Our genome-wide mRNA secondary structure analysis indicated that operonic mRNAs are comprised of ORF-wide units of secondary structure that vary across ORF boundaries such that adjacent ORFs on the same mRNA molecule are structurally distinct. ORF translation rate is strongly correlated with its mRNA structure in vivo, and correlation persists, albeit in a reduced form, with its structure when translation is inhibited and with that of in vitro refolded mRNA. These data suggest that intrinsic ORF mRNA structure encodes a rough blueprint for translation efficiency. This structure is then amplified by translation, in a self-reinforcing loop, to provide the structure that ultimately specifies the translation of each ORF.

The explosion of sequence information in bacteria makes developing high-throughput, cost-effective approaches to matching genes with phenotypes imperative. Using E. coli as proof of principle, we show that combining large-scale chemical genomics with quantitative fitness measurements provides a high-quality data set rich in discovery. Probing growth profiles of a mutant library in hundreds of conditions in parallel yielded > 10,000 phenotypes that allowed us to study gene essentiality, discover leads for gene function and drug action, and understand higher-order organization of the bacterial chromosome. We highlight new information derived from the study, including insights into a gene involved in multiple antibiotic resistance and the synergy between a broadly used combinatory antibiotic therapy, trimethoprim and sulfonamides. This data set, publicly available at http://ecoliwiki.net/tools/chemgen/, is a valuable resource for both the microbiological and bioinformatic communities, as it provides high-confidence associations between hundreds of annotated and uncharacterized genes as well as inferences about the mode of action of several poorly understood drugs.

Cells react to their environment through gene regulatory networks. Network integrity requires minimization of undesired crosstalk between their biomolecules. Similar constraints also limit the use of regulators when building synthetic circuits for engineering applications. Here, we mapped the promoter specificities of extracytoplasmic function (ECF) σs as well as the specificity of their interaction with anti-σs. DNA synthesis was used to build 86 ECF σs (two from every subgroup), their promoters, and 62 anti-σs identified from the genomes of diverse bacteria. A subset of 20 σs and promoters were found to be highly orthogonal to each other. This set can be increased by combining the -35 and -10 binding domains from different subgroups to build chimeras that target sequences unrepresented in any subgroup. The orthogonal σs, anti-σs, and promoters were used to build synthetic genetic switches in Escherichia coli. This represents a genome-scale resource of the properties of ECF σs and a resource for synthetic biology, where this set of well-characterized regulatory parts will enable the construction of sophisticated gene expression programs.

Functional genomics efforts face tradeoffs between number of perturbations examined and complexity of phenotypes measured. We bridge this gap with Perturb-seq, which combines droplet-based single-cell RNA-seq with a strategy for barcoding CRISPR-mediated perturbations, allowing many perturbations to be profiled in pooled format. We applied Perturb-seq to dissect the mammalian unfolded protein response (UPR) using single and combinatorial CRISPR perturbations. Two genome-scale CRISPR interference (CRISPRi) screens identified genes whose repression perturbs ER homeostasis. Subjecting ∼100 hits to Perturb-seq enabled high-precision functional clustering of genes. Single-cell analyses decoupled the three UPR branches, revealed bifurcated UPR branch activation among cells subject to the same perturbation, and uncovered differential activation of the branches across hits, including an isolated feedback loop between the translocon and IRE1α. These studies provide insight into how the three sensors of ER homeostasis monitor distinct types of stress and highlight the ability of Perturb-seq to dissect complex cellular responses.

Fighting antibiotic resistance requires a deeper understanding of the genetic factors that determine the antibiotic susceptibility of bacteria. Here we describe a chemical-genomic screen in Escherichia coli K-12 that was designed to discover new aspects of antibiotic resistance by focusing on a set of 26 antibiotics and other stresses with poorly characterized mode-of-action and determinants of resistance. We show that the screen identifies new resistance determinants for these antibiotics including a common signature from two antimicrobials, kasugamycin and blasticidin S, used to treat crop diseases like rice blast and fire blight. Following this signature, we further investigated the mechanistic basis for susceptibility to kasugamycin and blasticidin S in E. coli using both genetic and biochemical approaches. We provide evidence that these compounds hijack an overlapping set of peptide ABC-importers to enter the bacterial cell. Loss of uptake may be an underappreciated mechanism for the development of kasugamycin resistance in bacterial plant pathogens.

In enteric bacteria, the transcription factor σ(E) maintains membrane homeostasis by inducing synthesis of proteins involved in membrane repair and two small regulatory RNAs (sRNAs) that down-regulate synthesis of abundant membrane porins. Here, we describe the discovery of a third σ(E)-dependent sRNA, MicL (mRNA-interfering complementary RNA regulator of Lpp), transcribed from a promoter located within the coding sequence of the cutC gene. MicL is synthesized as a 308-nucleotide (nt) primary transcript that is processed to an 80-nt form. Both forms possess features typical of Hfq-binding sRNAs but surprisingly target only a single mRNA, which encodes the outer membrane lipoprotein Lpp, the most abundant protein of the cell. We show that the copper sensitivity phenotype previously ascribed to inactivation of the cutC gene is actually derived from the loss of MicL and elevated Lpp levels. This observation raises the possibility that other phenotypes currently attributed to protein defects are due to deficiencies in unappreciated regulatory RNAs. We also report that σ(E) activity is sensitive to Lpp abundance and that MicL and Lpp comprise a new σ(E) regulatory loop that opposes membrane stress. Together MicA, RybB, and MicL allow σ(E) to repress the synthesis of all abundant outer membrane proteins in response to stress.