The Essential Genome of Escherichia coli K-12

The intrinsic macrolide resistome of Escherichia coli

Intrinsic resistance to macrolides in Gram-negative bacteria is primarily attributed to the low permeability of the outer membrane, though the underlying genetic and molecular mechanisms remain to be fully elucidated. Here, we used transposon directed insertion-site sequencing (TraDIS) to identify chromosomal non-essential genes involved in Escherichia coli intrinsic resistance to a macrolide antibiotic, tilmicosin. We constructed two highly saturated transposon mutant libraries of >290,000 and >390,000 unique Tn5 insertions in a clinical enterotoxigenic strain (ETEC5621) and in a laboratory strain (K-12 MG1655), respectively. TraDIS analysis identified genes required for growth of ETEC5621 and MG1655 under 1/8 MIC (n = 15 and 16, respectively) and 1/4 MIC (n = 38 and 32, respectively) of tilmicosin. For both strains, 23 genes related to lipopolysaccharide biosynthesis, outer membrane assembly, the Tol-Pal system, efflux pump, and peptidoglycan metabolism were enriched in the presence of the antibiotic. Individual deletion of genes (n = 10) in the wild-type strains led to a 64- to 2-fold reduction in MICs of tilmicosin, erythromycin, and azithromycin, validating the results of the TraDIS analysis. Notably, deletion of surA or waaG, which impairs the outer membrane, led to the most significant decreases in MICs of all three macrolides in ETEC5621. Our findings contribute to a genome-wide understanding of intrinsic macrolide resistance in E. coli, shedding new light on the potential role of the peptidoglycan layer. They also provide an in vitro proof of concept that E. coli can be sensitized to macrolides by targeting proteins maintaining the outer membrane such as SurA and WaaG.

Synthetic genomes unveil the effects of synonymous recoding

Engineering the genetic code of an organism provides the basis for (i) making any organism safely resistant to natural viruses and (ii) preventing genetic information flow into and out of genetically modified organisms while (iii) allowing the biosynthesis of genetically encoded unnatural polymers 1-4 . Achieving these three goals requires the reassignment of multiple of the 64 codons nature uses to encode proteins. However, synonymous codon replacement-recoding-is frequently lethal, and how recoding impacts fitness remains poorly explored. Here, we explore these effects using whole-genome synthesis, multiplexed directed evolution, and genome-transcriptome-translatome-proteome co-profiling on multiple recoded genomes. Using this information, we assemble a synthetic Escherichia coli genome in seven sections using only 57 codons to encode proteins. By discovering the rules responsible for the lethality of synonymous recoding and developing a data-driven multi-omics-based genome construction workflow that troubleshoots synthetic genomes, we overcome the lethal effects of 62,007 synonymous codon swaps and 11,108 additional genomic edits. We show that synonymous recoding induces transcriptional noise including new antisense RNAs, leading to drastic transcriptome and proteome perturbation. As the elimination of select codons from an organism's genetic code results in the widespread appearance of cryptic promoters, we show that synonymous codon choice may naturally evolve to minimize transcriptional noise. Our work provides the first genome-scale description of how synonymous codon changes influence organismal fitness and paves the way for the construction of functional genomes that provide genetic firewalls from natural ecosystems and safely produce biopolymers, drugs, and enzymes with an expanded chemistry.

Bacterial genome reduction for optimal chassis of synthetic biology: a review

Bacteria with streamlined genomes, that harbor full functional genes for essential metabolic networks, are able to synthesize the desired products more effectively and thus have advantages as production platforms in industrial applications. To obtain streamlined chassis genomes, a large amount of effort has been made to reduce existing bacterial genomes. This work falls into two categories: rational and random reduction. The identification of essential gene sets and the emergence of various genome-deletion techniques have greatly promoted genome reduction in many bacteria over the past few decades. Some of the constructed genomes possessed desirable properties for industrial applications, such as: increased genome stability, transformation capacity, cell growth, and biomaterial productivity. The decreased growth and perturbations in physiological phenotype of some genome-reduced strains may limit their applications as optimized cell factories. This review presents an assessment of the advancements made to date in bacterial genome reduction to construct optimal chassis for synthetic biology, including: the identification of essential gene sets, the genome-deletion techniques, the properties and industrial applications of artificially streamlined genomes, the obstacles encountered in constructing reduced genomes, and the future perspectives.

Noise robustness and metabolic load determine the principles of central dogma regulation

The processes of gene expression are inherently stochastic, even for essential genes required for growth. How does the cell maximize fitness in light of noise? To answer this question, we build a mathematical model to explore the trade-off between metabolic load and growth robustness. The model predicts novel principles of central dogma regulation: Optimal protein expression levels for many genes are in vast overabundance. Essential genes are transcribed above a lower limit of one message per cell cycle. Gene expression is achieved by load balancing between transcription and translation. We present evidence that each of these novel regulatory principles is observed. These results reveal that robustness and metabolic load determine the global regulatory principles that govern central dogma processes, and these principles have broad implications for cellular function.

Inducible transposon mutagenesis for genome-scale forward genetics

Transposon insertion sequencing (Tn-seq) is a powerful method for genome-scale functional genetics in bacteria. However, its effectiveness is often limited by a lack of mutant diversity, caused by either inefficient transposon delivery or stochastic loss of mutants due to population bottlenecks. Here, we introduce "InducTn-seq", which leverages inducible mutagenesis for temporal control of transposition. InducTn-seq generates millions of transposon mutants from a single colony, enabling the sensitive detection of subtle fitness defects and transforming binary classifications of gene essentiality into a quantitative fitness measurement across both essential and non-essential genes. Using a mouse model of infectious colitis, we show that InducTn-seq bypasses a highly restrictive host bottleneck to generate a diverse transposon mutant population from the few cells that initiate infection, revealing the role of oxygen-related metabolic plasticity in pathogenesis. Overall, InducTn-seq overcomes the limitations of traditional Tn-seq, unlocking new possibilities for genome-scale forward genetic screens in bacteria.

RpoS activates formation of Salmonella Typhi biofilms and drives persistence in the gall bladder

The development of strategies for targeting the asymptomatic carriage of Salmonella Typhi in chronic typhoid patients has suffered owing to our basic lack of understanding of the molecular mechanisms that enable the formation of S. Typhi biofilms. Traditionally, studies have relied on cholesterol-attached biofilms formed by a closely related serovar, Typhimurium, to mimic multicellular Typhi communities formed on human gallstones. In long-term infections, S. Typhi adopts the biofilm lifestyle to persist in vivo and survive in the carrier state, ultimately leading to the spread of infections via the fecal-oral route of transmission. In the present work, we studied S. Typhi biofilms directly, applied targeted as well as genome-wide genetic approaches to uncover unique biofilm components that do not conform to the CsgD-dependent pathway established in S. Typhimurium. We undertook a genome-wide Tn5 mutation screen in a highly successful parental lineage of S. Typhi, strain H58, in gallstone-mimicking conditions. We generated New Generation Sequencing libraries based on the ClickSeq technology to identify the key regulators, IraP and RpoS, and the matrix components as Sth fimbriae, Vi capsule and lipopolysaccharide. We discovered that the starvation sigma factor, RpoS, was required for the transcriptional activation of matrix-encoding genes in vitro, and for S. Typhi colonization in persistent infections in vivo, using a heterologous fish larval model. An rpoS null mutant failed to colonize the gall bladder in chronic zebrafish infections. Overall, our work uncovered a novel RpoS-driven paradigm for the formation of cholesterol-attached Typhi biofilms, and emphasized the role(s) of stress signaling pathways for adaptation in chronic infections. Our identification of the biofilm regulators in S. Typhi paves the way for the development of drugs against typhoid carriage, which will ultimately control the increased incidence of gall bladder cancer in typhoid carriers.

Genetic Analysis of the ts-Lethal Mutant Δpa0665/pTS-pa0665 Reveals Its Role in Cell Morphology and Oxidative Phosphorylation in Pseudomonas aeruginosa

Pa0665 in Pseudomonas aeruginosa shares homologous sequences with that of the essential A-type iron-sulfur (Fe-S) cluster insertion protein ErpA in Escherichia coli. However, its essentiality in P. aeruginosa and its complementation with E. coli erpA has not been experimentally examined. To fulfill this task, we constructed plasmid-based ts-mutant Δpa0665/pTS-pa0665 using a three-step protocol. The mutant displayed growth defects at 42 °C, which were complemented by expressing ec.erpA. Microscopic observations indicated a petite cell phenotype for Δpa0665/pTS-pa0665 at 42 °C, correlated with the downregulation of the oprG gene. RNA sequencing revealed significant transcriptional changes in genes associated with the oxidative phosphorylation (OXPHOS) system, aligning with reduced ATP levels in Δpa0665/pTS-pa0665 under 42 °C. Additionally, the ts-mutant showed heightened sensitivity to H2O2 at 42 °C. Overall, our study demonstrates the essential role of pa0665 for OXPHOS function and is complemented by ec.erpA. We propose that the plasmid-based ts-allele is useful for genetic analysis of essential genes of interest in P. aeruginosa.

Two synonymous single-nucleotide polymorphisms promoting fluoroquinolone resistance of Escherichia coli in the environment

Single-nucleotide polymorphism (SNP) is one of the core mechanisms that respond to antibiotic resistance of Escherichia coli (E. coli), which is a major issue in environmental pollution. A specific type of SNPs, synonymous SNPs, have been generally considered as the "silent" SNPs since they do not change the encoded amino acid. However, the impact of synonymous SNPs on mRNA splicing, nucleo-cytoplasmic export, stability, and translation was gradually discovered in the last decades. Figuring out the mechanism of synonymous SNPs in regulating antibiotic resistance is critical to improve antimicrobial therapy strategies in clinics and biological treatment strategies of antibiotic-resistant E. coli-polluted materials. With our newly designed antibiotic resistant SNPs prediction algorithm, Multilocus Sequence Type based Identification for Phenotype-single nucleotide polymorphism Analysis (MIPHA), and in vivo validation, we identified 2 important synonymous SNPs 522 G>A and 972 C>T, located at hisD gene, which was previously predicted as a fluoroquinolone resistance-related gene without a detailed mechanism in the E. coli samples with environmental backgrounds. We first discovered that hisD causes gyrA mutation via the upregulation of sbmC and its downstream gene umuD. Moreover, those 2 synonymous SNPs of hisD cause its own translational slowdown and further reduce the expression levels of sbmC and its downstream gene umuD, making the fluoroquinolone resistance determining region of gyrA remains unmutated, ultimately causing the bacteria to lose their ability to resist drugs. This study provided valuable insight into the role of synonymous SNPs in mediating antibiotic resistance of bacteria and a new perspective for the treatment of environmental pollution caused by drug-resistant bacteria.

Long range segmentation of prokaryotic genomes by gene age and functionality

Bacterial and archaeal genomes encompass numerous operons that typically consist of two to five genes. On larger scales, however, gene order is poorly conserved through the evolution of prokaryotes. Nevertheless, non-random localization of different classes of genes on prokaryotic chromosomes could reflect important functional and evolutionary constraints. We explored the patterns of genomic localization of evolutionarily conserved (ancient) and variable (young) genes across the diversity of bacteria and archaea. Nearly all bacterial and archaeal chromosomes were found to encompass large segments of 100-300 kilobases that were significantly enriched in either ancient or young genes. Similar clustering of genes with lethal knockout phenotype (essential genes) was observed as well. Mathematical modeling of genome evolution suggests that this long-range gene clustering in prokaryotic chromosomes reflects perpetual genome rearrangement driven by a combination of selective and neutral processes rather than evolutionary conservation.

Expanding the flexibility of base editing for high-throughput genetic screens in bacteria

Genome-wide screens have become powerful tools for elucidating genotype-to-phenotype relationships in bacteria. Of the varying techniques to achieve knockout and knockdown, CRISPR base editors are emerging as promising options. However, the limited number of available, efficient target sites hampers their use for high-throughput screening. Here, we make multiple advances to enable flexible base editing as part of high-throughput genetic screening in bacteria. We first co-opt the Streptococcus canis Cas9 that exhibits more flexible protospacer-adjacent motif recognition than the traditional Streptococcus pyogenes Cas9. We then expand beyond introducing premature stop codons by mutating start codons. Next, we derive guide design rules by applying machine learning to an essentiality screen conducted in Escherichia coli. Finally, we rescue poorly edited sites by combining base editing with Cas9-induced cleavage of unedited cells, thereby enriching for intended edits. The efficiency of this dual system was validated through a conditional essentiality screen based on growth in minimal media. Overall, expanding the scope of genome-wide knockout screens with base editors could further facilitate the investigation of new gene functions and interactions in bacteria.

Technical considerations for cost-effective transposon directed insertion-site sequencing (TraDIS)

Transposon directed insertion-site sequencing (TraDIS), a variant of transposon insertion sequencing commonly known as Tn-Seq, is a high-throughput assay that defines essential bacterial genes across diverse growth conditions. However, the variability between laboratory environments often requires laborious, time-consuming modifications to its protocol. In this technical study, we aimed to refine the protocol by identifying key parameters that can impact the complexity of mutant libraries. Firstly, we discovered that adjusting electroporation parameters including transposome concentration, transposome assembly conditions, and cell densities can significantly improve the recovery of viable mutants for different Escherichia coli strains. Secondly, we found that post-electroporation conditions, such as recovery time and the use of different mediums for selecting mutants may also impact the complexity of viable mutants in the library. Finally, we developed a simplified sequencing library preparation workflow based on a Nextera-TruSeq hybrid design where ~ 80% of sequenced reads correspond to transposon-DNA junctions. The technical improvements presented in our study aim to streamline TraDIS protocols, making this powerful technique more accessible for a wider scientific audience.

A Targeted Genome-scale Overexpression Platform for Proteobacteria

Targeted, genome-scale gene perturbation screens using Clustered Regularly Interspaced Short Palindromic Repeats interference (CRISPRi) and activation (CRISPRa) have revolutionized eukaryotic genetics, advancing medical, industrial, and basic research. Although CRISPRi knockdowns have been broadly applied in bacteria, options for genome-scale overexpression face key limitations. Here, we develop a facile approach for genome-scale gene overexpression in bacteria we call, "CRISPRtOE" (CRISPR transposition and OverExpression). We create a platform for comprehensive gene targeting using CRISPR-associated transposition (CAST) and show that transposition occurs at a higher frequency in non-transcribed DNA. We then demonstrate that CRISPRtOE can upregulate gene expression in Proteobacteria with medical and industrial relevance by integrating synthetic promoters of varying strength upstream of target genes. Finally, we employ CRISPRtOE screening at the genome-scale in Escherichia coli, recovering known antibiotic targets and genes with unexplored roles in antibiotic function. We envision that CRISPRtOE will be a valuable overexpression tool for antibiotic mode of action, industrial strain optimization, and gene function discovery in bacteria.

Extraribosomal Functions of Bacterial Ribosomal Proteins-An Update, 2023

Ribosomal proteins (r-proteins) are abundant, highly conserved, and multifaceted cellular proteins in all domains of life. Most r-proteins have RNA-binding properties and can form protein-protein contacts. Bacterial r-proteins govern the co-transcriptional rRNA folding during ribosome assembly and participate in the formation of the ribosome functional sites, such as the mRNA-binding site, tRNA-binding sites, the peptidyl transferase center, and the protein exit tunnel. In addition to their primary role in a cell as integral components of the protein synthesis machinery, many r-proteins can function beyond the ribosome (the phenomenon known as moonlighting), acting either as individual regulatory proteins or in complexes with various cellular components. The extraribosomal activities of r-proteins have been studied over the decades. In the past decade, our understanding of r-protein functions has advanced significantly due to intensive studies on ribosomes and gene expression mechanisms not only in model bacteria like Escherichia coli or Bacillus subtilis but also in little-explored bacterial species from various phyla. The aim of this review is to update information on the multiple functions of r-proteins in bacteria.

Machine learning identifies key metabolic reactions in bacterial growth on different carbon sources

Carbon source-dependent control of bacterial growth is fundamental to bacterial physiology and survival. However, pinpointing the metabolic steps important for cell growth is challenging due to the complexity of cellular networks. Here, the elastic net model and multilayer perception model that integrated genome-wide gene-deletion data and simulated flux distributions were constructed to identify metabolic reactions beneficial or detrimental to Escherichia coli grown on 30 different carbon sources. Both models outperformed traditional in silico methods by identifying not just essential reactions but also nonessential ones that promote growth. They successfully predicted metabolic reactions beneficial to cell growth, with high convergence between the models. The models revealed that biosynthetic pathways generally promote growth across various carbon sources, whereas the impact of energy-generating pathways varies with the carbon source. Intriguing predictions were experimentally validated for findings beyond experimental training data and the impact of various carbon sources on the glyoxylate shunt, pyruvate dehydrogenase reaction, and redundant purine biosynthesis reactions. These highlight the practical significance and predictive power of the models for understanding and engineering microbial metabolism.

Reduction-to-synthesis: the dominant approach to genome-scale synthetic biology

Advances in systems and synthetic biology have propelled the construction of reduced bacterial genomes. Genome reduction was initially focused on exploring properties of minimal genomes, but more recently it has been deployed as an engineering strategy to enhance strain performance. This review provides the latest updates on reduced genomes, focusing on dual-track approaches of top-down reduction and bottom-up synthesis for their construction. Using cases from studies that are based on established industrial workhorse strains, we discuss the construction of a series of synthetic phenotypes that are candidates for biotechnological applications. Finally, we address the possible uses of reduced genomes for biotechnological applications and the needed future research directions that may ultimately lead to the total synthesis of rationally designed genomes.

Genome-wide identification of fitness-genes in aminoglycoside-resistant Escherichia coli during antibiotic stress

Resistance against aminoglycosides is widespread in bacteria. This study aimed to identify genes that are important for growth of E. coli during aminoglycoside exposure, since such genes may be targeted to re-sensitize resistant E. coli to treatment. We constructed three transposon mutant libraries each containing > 230.000 mutants in E. coli MG1655 strains harboring streptomycin (aph(3″)-Ib/aph(6)-Id), gentamicin (aac(3)-IV), or neomycin (aph(3″)-Ia) resistance gene(s). Transposon Directed Insertion-site Sequencing (TraDIS), a combination of transposon mutagenesis and high-throughput sequencing, identified 56 genes which were deemed important for growth during streptomycin, 39 during gentamicin and 32 during neomycin exposure. Most of these fitness-genes were membrane-located (n = 55) and involved in either cell division, ATP-synthesis or stress response in the streptomycin and gentamicin exposed libraries, and enterobacterial common antigen biosynthesis or magnesium sensing/transport in the neomycin exposed library. For validation, eight selected fitness-genes/gene-clusters were deleted (minCDE, hflCK, clsA and cpxR associated with streptomycin and gentamicin resistance, and phoPQ, wecA, lpp and pal associated with neomycin resistance), and all mutants were shown to be growth attenuated upon exposure to the corresponding antibiotics. In summary, we identified genes that are advantageous in aminoglycoside-resistant E. coli during antibiotic stress. In addition, we increased the understanding of how aminoglycoside-resistant E. coli respond to antibiotic exposure.

Changing fitness effects of mutations through long-term bacterial evolution

The distribution of fitness effects of new mutations shapes evolution, but it is challenging to observe how it changes as organisms adapt. Using Escherichia coli lineages spanning 50,000 generations of evolution, we quantify the fitness effects of insertion mutations in every gene. Macroscopically, the fraction of deleterious mutations changed little over time whereas the beneficial tail declined sharply, approaching an exponential distribution. Microscopically, changes in individual gene essentiality and deleterious effects often occurred in parallel; altered essentiality is only partly explained by structural variation. The identity and effect sizes of beneficial mutations changed rapidly over time, but many targets of selection remained predictable because of the importance of loss-of-function mutations. Taken together, these results reveal the dynamic-but statistically predictable-nature of mutational fitness effects.

Randomly barcoded transposon mutant libraries for gut commensals I: Strategies for efficient library construction

Randomly barcoded transposon mutant libraries are powerful tools for studying gene function and organization, assessing gene essentiality and pathways, discovering potential therapeutic targets, and understanding the physiology of gut bacteria and their interactions with the host. However, construction of high-quality libraries with uniform representation can be challenging. In this review, we survey various strategies for barcoded library construction, including transposition systems, methods of transposon delivery, optimal library size, and transconjugant selection schemes. We discuss the advantages and limitations of each approach, as well as factors to consider when selecting a strategy. In addition, we highlight experimental and computational advances in arraying condensed libraries from mutant pools. We focus on examples of successful library construction in gut bacteria and their application to gene function studies and drug discovery. Given the need for understanding gene function and organization in gut bacteria, we provide a comprehensive guide for researchers to construct randomly barcoded transposon mutant libraries.

Mobile CRISPRi moves through the complexity of bacterial genetics

In this issue of Cell Reports Methods, Rachwalski et al. describe a high-throughput method to screen genetic interactions in bacteria using a conjugative CRISPR interference (CRISPRi) plasmid. The method enables systematic studies of gene essentiality in diverse genomic and environmental contexts and is applicable to Escherichia coli as well as other bacteria.

A mobile CRISPRi collection enables genetic interaction studies for the essential genes of Escherichia coli

Advances in gene editing, in particular CRISPR interference (CRISPRi), have enabled depletion of essential cellular machinery to study the downstream effects on bacterial physiology. Here, we describe the construction of an ordered E. coli CRISPRi collection, designed to knock down the expression of 356 essential genes with the induction of a catalytically inactive Cas9, harbored on the conjugative plasmid pFD152. This mobile CRISPRi library can be conjugated into other ordered genetic libraries to assess combined effects of essential gene knockdowns with non-essential gene deletions. As proof of concept, we probed cell envelope synthesis with two complementary crosses: (1) an Lpp deletion into every CRISPRi knockdown strain and (2) the lolA knockdown plasmid into the Keio collection. These experiments revealed a number of notable genetic interactions for the essential phenotype probed and, in particular, showed suppressing interactions for the loci in question.

Structural assembly of the bacterial essential interactome

The study of protein interactions in living organisms is fundamental for understanding biological processes and central metabolic pathways. Yet, our knowledge of the bacterial interactome remains limited. Here, we combined gene deletion mutant analysis with deep-learning protein folding using AlphaFold2 to predict the core bacterial essential interactome. We predicted and modeled 1402 interactions between essential proteins in bacteria and generated 146 high-accuracy models. Our analysis reveals previously unknown details about the assembly mechanisms of these complexes, highlighting the importance of specific structural features in their stability and function. Our work provides a framework for predicting the essential interactomes of bacteria and highlight the potential of deep-learning algorithms in advancing our understanding of the complex biology of living organisms. Also, the results presented here offer a promising approach to identify novel antibiotic targets.

Investigating Pseudomonas aeruginosa Gene Function During Pathogenesis Using Mobile-CRISPRi

CRISPR interference (CRISPRi) is a robust gene silencing technique that is ideal for targeting essential and conditionally essential (CE) genes. CRISPRi is especially valuable for investigating gene function in pathogens such as P. aeruginosa where essential and CE genes underlie clinically important phenotypes such as antibiotic susceptibility and virulence. To facilitate the use of CRISPRi in diverse bacteria-including P. aeruginosa-we developed a suite of modular, mobilizable, and integrating vectors we call, "Mobile-CRISPRi." We further optimized Mobile-CRISPRi for use in P. aeruginosa mouse models of acute lung infection by expressing the CRISPRi machinery at low levels constitutively, enabling partial knockdown of essential and CE genes without the need for an exogenous inducer. Here, we describe protocols for creating Mobile-CRISPRi knockdown strains and testing their phenotypes in a mouse pneumonia model of P. aeruginosa infection. In addition, we provide comprehensive guide RNA designs to target genes in common laboratory strains of P. aeruginosa and other Pseudomonas species.

Engineering of a Substrate Affinity Reduced S-Adenosyl-methionine Synthetase as a Novel Biosensor for Growth-Coupling Selection of L-Methionine Overproducers

Biosensors are powerful tools for monitoring specific metabolites or controlling metabolic flux towards the products in a single cell, which play important roles in microbial cell factory construction. Despite their potential role in metabolic flux monitoring, the development of biosensors for small molecules is still limited. Reported biosensors often exhibit bottlenecks of poor specificity and a narrow dynamic range. Moreover, fine-tuning the substrate binding affinity of a crucial enzyme can decrease its catalytic activity, which ultimately results in the repression of the corresponding essential metabolite biosynthesis and impairs cell growth. However, increasing intracellular substrate concentration can elevate the availability of the essential metabolite and may lead to restore cellular growth. Herein, a new strategy was proposed for constructing whole-cell biosensors based on enzyme encoded by essential gene that offer inherent specificity and universality. Specifically, S-adenosyl-methionine synthetase (MetK) in E. coli was chosen as the crucial enzyme, and a series of MetK variants were identified that were sensitive to L-methionine concentration. This occurrence enabled the engineered cell to sense L-methionine and exhibit L-methionine dose-dependent cell growth. To improve the biosensor's dynamic range, an S-adenosyl-methionine catabolic enzyme was overexpressed to reduce the intracellular availability of S-adenosyl-methionine. The resulting whole-cell biosensor effectively coupled the intracellular concentration of L-methionine with growth and was successfully applied to select strains with enhanced L-methionine biosynthesis from random mutagenesis libraries. Overall, our study presents a universal strategy for designing and constructing growth-coupled biosensors based on crucial enzyme, which can be applied to select strains overproducing high value-added metabolites in cellular metabolism.

Identification of recG genetic interactions in Escherichia coli by transposon sequencing

IMPORTANCE

DNA damage and subsequent DNA repair processes are mutagenic in nature and an important driver of evolution in prokaryotes, including antibiotic resistance development. Genetic screening approaches, such as transposon sequencing (Tn-seq), have provided important new insights into gene function and genetic relationships. Here, we employed Tn-seq to gain insight into the function of the recG gene, which renders Escherichia coli cells moderately sensitive to a variety of DNA-damaging agents when they are absent. The reported recG genetic interactions can be used in combination with future screens to aid in a more complete reconstruction of DNA repair pathways in bacteria.

Mechanism of target site selection by type V-K CRISPR-associated transposases

CRISPR-associated transposases (CASTs) repurpose nuclease-deficient CRISPR effectors to catalyze RNA-guided transposition of large genetic payloads. Type V-K CASTs offer potential technology advantages but lack accuracy, and the molecular basis for this drawback has remained elusive. Here, we reveal that type V-K CASTs maintain an RNA-independent, "untargeted" transposition pathway alongside RNA-dependent integration, driven by the local availability of TnsC filaments. Using cryo-electron microscopy, single-molecule experiments, and high-throughput sequencing, we found that a minimal, CRISPR-less transpososome preferentially directs untargeted integration at AT-rich sites, with additional local specificity imparted by TnsB. By exploiting this knowledge, we suppressed untargeted transposition and increased type V-K CAST specificity up to 98.1% in cells without compromising on-target integration efficiency. These findings will inform further engineering of CAST systems for accurate, kilobase-scale genome engineering applications.

Essential gene complement of Planctopirus limnophila from the bacterial phylum Planctomycetes

Planctopirus limnophila belongs to the bacterial phylum Planctomycetes, a relatively understudied lineage with remarkable cell biology features. Here, we report a genome-wide analysis of essential gene content in P. limnophila. We show that certain genes involved in peptidoglycan synthesis or cell division, which are essential in most other studied bacteria, are not essential for growth under laboratory conditions in this species. We identify essential genes likely involved in lipopolysaccharide biosynthesis, consistent with the view of Planctomycetes as diderm bacteria, and highlight other essential genes of unknown functions. Furthermore, we explore potential stages of evolution of the essential gene repertoire in Planctomycetes and the related phyla Verrucomicrobia and Chlamydiae. Our results provide insights into the divergent molecular and cellular biology of Planctomycetes.

Enhancing isoprenol production by systematically tuning metabolic pathways using CRISPR interference in E. coli

Regulation of metabolic gene expression is crucial for maximizing bioproduction titers. Recent engineering tools including CRISPR/Cas9, CRISPR interference (CRISPRi), and CRISPR activation (CRISPRa) have enabled effective knock-out, knock-down, and overexpression of endogenous pathway genes, respectively, for advanced strain engineering. CRISPRi in particular has emerged as a powerful tool for gene repression through the use of a deactivated Cas9 (dCas9) protein and target guide RNA (gRNA). By constructing gRNA arrays, CRISPRi has the capacity for multiplexed gene downregulation across multiple orthogonal pathways for enhanced bioproduction titers. In this study, we harnessed CRISPRi to downregulate 32 essential and non-essential genes in E. coli strains heterologously expressing either the original mevalonate pathway or isopentenyl diphosphate (IPP) bypass pathway for isoprenol biosynthesis. Isoprenol remains a candidate bioproduct both as a drop-in blend additive and as a precursor for the high-performance sustainable aviation fuel, 1,4-dimethylcyclooctane (DMCO). Of the 32 gRNAs targeting genes associated with isoprenol biosynthesis, a subset was found to vastly improve product titers. Construction of a multiplexed gRNA library based on single guide RNA (sgRNA) performance enabled simultaneous gene repression, yielding a 3 to 4.5-fold increase in isoprenol titer (1.82 ± 0.19 g/L) on M9-MOPS minimal medium. We then scaled the best performing CRISPRi strain to 2-L fed-batch cultivation and demonstrated translatable titer improvements, ultimately obtaining 12.4 ± 1.3 g/L isoprenol. Our strategy further establishes CRISPRi as a powerful tool for tuning metabolic flux in production hosts and that titer improvements are readily scalable with potential for applications in industrial bioproduction.

Microbial Primer: Transposon directed insertion site sequencing (TraDIS): A high throughput method for linking genotype to phenotype

Genetic screens are a key tool for linking phenotype and genotype. Transposon mutagenesis was one of the first genetic methodologies to associate genetic loci with phenotypes. The advent of next-generation sequencing transformed the use of this technique allowing rapid interrogation of whole genomes for genes that correlate with phenotype. One method is transposon directed insertion-site sequencing (TraDIS). Here we describe the method, recent developments in technology, and the advantages and disadvantages of this method compared to other genetic screening tools.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Construction and application of the conditionally essential gene knockdown library in Klebsiella pneumoniae to screen potential antimicrobial targets and virulence genes via Mobile-CRISPRi-seq

Klebsiella pneumoniae is a ubiquitous human pathogen, and its clinical treatment faces two major challenges: multidrug resistance and the pathogenesis of hypervirulent K. pneumoniae. The discovery and study of conditionally essential (CE) genes that can function as potential antimicrobial targets has always been a research concern due to their restriction in the development of novel antibiotics. However, the lack of essential functional genomic data has hampered the study of the mechanisms of essential genes related to antimicrobial susceptibility. In this study, we developed a pooled CE genes mobile clustered regularly interspaced short palindromic repeat (CRISPR) interference screening method (Mobile-CRISPRi-seq) for K. pneumoniae to identify genes that play critical roles in antimicrobial fitness in vitro and host immunity in vivo. Targeting 870 predicted CE genes in K. pneumoniae, Mobile-CRISPRi-seq uncovered the depletion of tetrahydrofolate synthesis pathway genes folB and folP under trimethoprim pressure. Our screening also identified genes waaE and fldA related to polymyxin and β-lactam susceptibility by applying a screening strategy based on Mobile-CRISPRi-seq and comparative genomics. Furthermore, using a mouse infection model and Mobile-CRISPRi-seq, multiple virulence genes were identified, and among these genes, pal, yciS, and ribB were demonstrated to contribute to the pathogenesis of K. pneumoniae. This study provides a simple, rapid, and effective platform for screening potential antimicrobial targets and virulence genes in K. pneumoniae, and this broadly applicable system can be expanded for high-throughput functional gene study in multiple pathogenic bacteria, especially in gram-negative bacteria. IMPORTANCE The discovery and investigation of conditionally essential (CE) genes that can function as potential antimicrobial targets has always been a research concern because of the restriction of antimicrobial targets in the development of novel antibiotics. In this study, we developed a pooled CE gene-wide mobile clustered regularly interspaced short palindromic repeat (CRISPR) interference sequencing (Mobile-CRISPRi-seq) strategy in Klebsiella pneumoniae to identify genes that play critical roles in the fitness of antimicrobials in vitro and host immunity in vivo. The data suggest a robust tool to screen for loss-of-function phenotypes in a pooled gene knockdown library in K. pneumoniae, and Mobile-CRISPRi-seq may be expanded to multiple bacteria for screening and identification of genes with crucial roles in the fitness of antimicrobials and hosts.

Strategies to identify and edit improvements in synthetic genome segments episomally

Genome engineering projects often utilize bacterial artificial chromosomes (BACs) to carry multi-kilobase DNA segments at low copy number. However, all stages of whole-genome engineering have the potential to impose mutations on the synthetic genome that can reduce or eliminate the fitness of the final strain. Here, we describe improvements to a multiplex automated genome engineering (MAGE) protocol to improve recombineering frequency and multiplexability. This protocol was applied to recoding an Escherichia coli strain to replace seven codons with synonymous alternatives genome wide. Ten 44 402-47 179 bp de novo synthesized DNA segments contained in a BAC from the recoded strain were unable to complement deletion of the corresponding 33-61 wild-type genes using a single antibiotic resistance marker. Next-generation sequencing (NGS) was used to identify 1-7 non-recoding mutations in essential genes per segment, and MAGE in turn proved a useful strategy to repair these mutations on the recoded segment contained in the BAC when both the recoded and wild-type copies of the mutated genes had to exist by necessity during the repair process. Finally, two web-based tools were used to predict the impact of a subset of non-recoding missense mutations on strain fitness using protein structure and function calls.

Mapping temperature-sensitive mutations at a genome scale to engineer growth switches in Escherichia coli

Temperature-sensitive (TS) mutants are a unique tool to perturb and engineer cellular systems. Here, we constructed a CRISPR library with 15,120 Escherichia coli mutants, each with a single amino acid change in one of 346 essential proteins. 1,269 of these mutants showed temperature-sensitive growth in a time-resolved competition assay. We reconstructed 94 TS mutants and measured their metabolism under growth arrest at 42°C using metabolomics. Metabolome changes were strong and mutant-specific, showing that metabolism of nongrowing E. coli is perturbation-dependent. For example, 24 TS mutants of metabolic enzymes overproduced the direct substrate metabolite due to a bottleneck in their associated pathway. A strain with TS homoserine kinase (ThrBF267D ) produced homoserine for 24 h, and production was tunable by temperature. Finally, we used a TS subunit of DNA polymerase III (DnaXL289Q ) to decouple growth from arginine overproduction in engineered E. coli. These results provide a strategy to identify TS mutants en masse and demonstrate their large potential to produce bacterial metabolites with nongrowing cells.

Advancements in synthetic biology-based bacterial cancer therapy: A modular design approach

Synthetic biology aims to program living bacteria cells with artificial genetic circuits for user-defined functions, transforming them into powerful tools with numerous applications in various fields, including oncology. Cancer treatments have serious side effects on patients due to the systemic action of the drugs involved. To address this, new systems that provide localized antitumoral action while minimizing damage to healthy tissues are required. Bacteria, often considered pathogenic agents, have been used as cancer treatments since the early 20th century. Advances in genetic engineering, synthetic biology, microbiology, and oncology have improved bacterial therapies, making them safer and more effective. Here we propose six modules for a successful synthetic biology-based bacterial cancer therapy, the modules include Payload, Release, Tumor-targeting, Biocontainment, Memory, and Genetic Circuit Stability Module. These will ensure antitumor activity, safety for the environment and patient, prevent bacterial colonization, maintain cell stability, and prevent loss or defunctionalization of the genetic circuit.

Genome-wide promoter responses to CRISPR perturbations of regulators reveal regulatory networks in Escherichia coli

Elucidating genome-scale regulatory networks requires a comprehensive collection of gene expression profiles, yet measuring gene expression responses for every transcription factor (TF)-gene pair in living prokaryotic cells remains challenging. Here, we develop pooled promoter responses to TF perturbation sequencing (PPTP-seq) via CRISPR interference to address this challenge. Using PPTP-seq, we systematically measure the activity of 1372 Escherichia coli promoters under single knockdown of 183 TF genes, illustrating more than 200,000 possible TF-gene responses in one experiment. We perform PPTP-seq for E. coli growing in three different media. The PPTP-seq data reveal robust steady-state promoter activities under most single TF knockdown conditions. PPTP-seq also enables identifications of, to the best of our knowledge, previously unknown TF autoregulatory responses and complex transcriptional control on one-carbon metabolism. We further find context-dependent promoter regulation by multiple TFs whose relative binding strengths determined promoter activities. Additionally, PPTP-seq reveals different promoter responses in different growth media, suggesting condition-specific gene regulation. Overall, PPTP-seq provides a powerful method to examine genome-wide transcriptional regulatory networks and can be potentially expanded to reveal gene expression responses to other genetic elements.

Ribosome inactivation by Escherichia coli GTPase RsgA inhibits T4 phage

Introduction

Bacteria must combat phages, and myriad bacterial anti-phage systems have been discovered that reduce host metabolism, for example, by depleting energetic compounds like ATP and NAD+. Hence, these systems indirectly inhibit protein production. Surprisingly, direct reduction of ribosome activity has not been demonstrated to thwart phage.

Methods

Here, by producing each of the 4,287 Escherichia coli proteins and selecting for anti-phage activity that leads to enhanced growth, we investigated the role of host proteins in phage inhibition.

Results and discussion

We identified that E. coli GTPase RsgA inhibits lytic phage T4 by inactivating ribosomes.

Cell Free Bacteriophage Synthesis from Engineered Strains Improves Yield

Phage therapy to treat life-threatening drug-resistant infections has been hampered by technical challenges in phage production. Cell-free bacteriophage synthesis (CFBS) can overcome the limitations of standard phage production methods by manufacturing phage virions in vitro. CFBS mimics intracellular phage assembly using transcription/translation machinery (TXTL) harvested from bacterial lysates and combined with reagents to synthesize proteins encoded by a phage genomic DNA template. These systems may enable rapid phage production and engineering to accelerate phages from bench-to-bedside. TXTL harvested from wild type or commonly used bacterial strains was not optimized for bacteriophage production. Here, we demonstrate that TXTL from genetically modified E. coli BL21 can be used to enhance phage T7 yields in vitro by CFBS. Expression of 18 E. coli BL21 genes was manipulated by inducible CRISPR interference (CRISPRi) mediated by nuclease deficient Cas12a from F. novicida (dFnCas12a) to identify genes implicated in T7 propagation as positive or negative effectors. Genes shown to have a significant effect were overexpressed (positive effectors) or repressed (negative effectors) to modify the genetic background of TXTL harvested for CFBS. Phage T7 CFBS yields were improved by up to 10-fold in vitro through overexpression of translation initiation factor IF-3 (infC) and small RNAs OxyS and CyaR and by repression of RecC subunit exonuclease RecBCD. Continued improvement of CFBS will mitigate phage manufacturing bottlenecks and lower hurdles to widespread adoption of phage therapy.

Mechanism of target site selection by type V-K CRISPR-associated transposases

Unlike canonical CRISPR-Cas systems that rely on RNA-guided nucleases for target cleavage, CRISPR-associated transposases (CASTs) repurpose nuclease-deficient CRISPR effectors to facilitate RNA-guided transposition of large genetic payloads. Type V-K CASTs offer several potential upsides for genome engineering, due to their compact size, easy programmability, and unidirectional integration. However, these systems are substantially less accurate than type I-F CASTs, and the molecular basis for this difference has remained elusive. Here we reveal that type V-K CASTs undergo two distinct mobilization pathways with remarkably different specificities: RNA-dependent and RNA-independent transposition. Whereas RNA-dependent transposition relies on Cas12k for accurate target selection, RNA-independent integration events are untargeted and primarily driven by the local availability of TnsC filaments. The cryo-EM structure of the untargeted complex reveals a TnsB-TnsC-TniQ transpososome that encompasses two turns of a TnsC filament and otherwise resembles major architectural aspects of the Cas12k-containing transpososome. Using single-molecule experiments and genome-wide meta-analyses, we found that AT-rich sites are preferred substrates for untargeted transposition and that the TnsB transposase also imparts local specificity, which collectively determine the precise insertion site. Knowledge of these motifs allowed us to direct untargeted transposition events to specific hotspot regions of a plasmid. Finally, by exploiting TnsB's preference for on-target integration and modulating the availability of TnsC, we suppressed RNA-independent transposition events and increased type V-K CAST specificity up to 98.1%, without compromising the efficiency of on-target integration. Collectively, our results reveal the importance of dissecting target site selection mechanisms and highlight new opportunities to leverage CAST systems for accurate, kilobase-scale genome engineering applications.

The one-message-per-cell-cycle rule: A conserved minimum transcription level for essential genes

The inherent stochasticity of cellular processes leads to significant cell-to-cell variation in protein abundance. Although this noise has already been characterized and modeled, its broader implications and significance remain unclear. In this paper, we revisit the noise model and identify the number of messages transcribed per cell cycle as the critical determinant of noise. In yeast, we demonstrate that this quantity predicts the non-canonical scaling of noise with protein abundance, as well as quantitatively predicting its magnitude. We then hypothesize that growth robustness requires an upper ceiling on noise for the expression of essential genes, corresponding to a lower floor on the transcription level. We show that just such a floor exists: a minimum transcription level of one message per cell cycle is conserved between three model organisms: Escherichia coli, yeast, and human. Furthermore, all three organisms transcribe the same number of messages per gene, per cell cycle. This common transcriptional program reveals that robustness to noise plays a central role in determining the expression level of a large fraction of essential genes, and that this fundamental optimal strategy is conserved from E. coli to human cells.

The one-message-per-cell-cycle rule: A conserved minimum transcription level for essential genes

The inherent stochasticity of cellular processes leads to significant cell-to-cell variation in protein abundance. Although this noise has already been characterized and modeled, its broader implications and significance remain unclear. In this paper, we revisit the noise model and identify the number of messages transcribed per cell cycle as the critical determinant of noise. In yeast, we demonstrate that this quantity predicts the non-canonical scaling of noise with protein abundance, as well as quantitatively predicting its magnitude. We then hypothesize that growth robustness requires an upper ceiling on noise for the expression of essential genes, corresponding to a lower floor on the transcription level. We show that just such a floor exists: a minimum transcription level of one message per cell cycle is conserved between three model organisms: Escherichia coli, yeast, and human. Furthermore, all three organisms transcribe the same number of messages per gene, per cell cycle. This common transcriptional program reveals that robustness to noise plays a central role in determining the expression level of a large fraction of essential genes, and that this fundamental optimal strategy is conserved from E. coli to human cells.

CRISPR-Associated Transposase for Targeted Mutagenesis in Diverse Proteobacteria

Genome editing tools, through the disruption of an organism's native genetic material or the introduction of non-native DNA, facilitate functional investigations to link genotypes to phenotypes. Transposons have been instrumental genetic tools in microbiology, enabling genome-wide, randomized disruption of genes and insertions of new genetic elements. Due to this randomness, identifying and isolating particular transposon mutants (i.e., those with modifications at a genetic locus of interest) can be laborious, often requiring one to sift through hundreds or thousands of mutants. Programmable, site-specific targeting of transposons became possible with recently described CRISPR-associated transposase (CASTs) systems, allowing the streamlined recovery of desired mutants in a single step. Like other CRISPR-derived systems, CASTs can be programmed by guide-RNA that is transcribed from short DNA sequence(s). Here, we describe a CAST system and demonstrate its function in bacteria from three classes of Proteobacteria. A dual plasmid strategy is demonstrated: (i) CAST genes are expressed from a broad-host-range replicative plasmid and (ii) guide-RNA and transposon are encoded on a high-copy, suicidal pUC plasmid. Using our CAST system, single-gene disruptions were performed with on-target efficiencies approaching 100% in Beta- and Gammaproteobacteria (Burkholderia thailandensis and Pseudomonas putida, respectively). We also report a peak efficiency of 45% in the Alphaproteobacterium Agrobacterium fabrum. In B. thailandensis, we performed simultaneous co-integration of transposons at two different target sites, demonstrating CAST's utility in multilocus strategies. The CAST system is also capable of high-efficiency large transposon insertion totaling over 11 kbp in all three bacteria tested. Lastly, the dual plasmid system allowed for iterative transposon mutagenesis in all three bacteria without loss of efficiency. Given these iterative capabilities and large payload capacity, this system will be helpful for genome engineering experiments across several fields of research.

Analysis of the Plasmid-Based ts-Mutant ΔfabA/pTS-fabA Reveals Its Lethality under Aerobic Growth Conditions That Is Suppressed by Mild Overexpression of desA at a Restrictive Temperature in Pseudomonas aeruginosa

It is uncertain whether PA1610|fabA is essential or dispensable for growth on LB-agar plates under aerobic conditions in Pseudomonas aeruginosa PAO1. To examine its essentiality, we disrupted fabA in the presence of a native promoter-controlled complementary copy on ts-plasmid. In this analysis, we showed that the plasmid-based ts-mutant ΔfabA/pTS-fabA failed to grow at a restrictive temperature, consistent with the observation by Hoang and Schweizer (T. T. Hoang, H. P. Schweizer, J Bacteriol 179:5326-5332, 1997, https://doi.org/10.1128/jb.179.17.5326-5332.1997), and expanded on this by showing that ΔfabA exhibited curved cell morphology. On the other hand, strong induction of fabA-OE or PA3645|fabZ-OE impeded the growth of cells displaying oval morphology. Suppressor analysis revealed a mutant sup gene that suppressed a growth defect but not cell morphology of ΔfabA. Genome resequencing and transcriptomic profiling of sup identified PA0286|desA, whose promoter carried a single-nucleotide polymorphism (SNP), and transcription was significantly upregulated (level increase of >2-fold, P < 0.05). By integration of the SNP-bearing promoter-controlled desA gene into the chromosome of ΔfabA/pTS-fabA, we showed that the SNP is sufficient for ΔfabA to phenocopy the sup mutant. Furthermore, mild induction of the araC-PBAD-controlled desA gene but not desB rescued ΔfabA. These results validated that mild overexpression of desA fully suppressed the lethality but not the curved cell morphology of ΔfabA. Similarly, Zhu et al. (Zhu K, Choi K-H, Schweizer HP, Rock CO, Zhang Y-M, Mol Microbiol 60:260-273, 2006, https://doi.org/10.1111/j.1365-2958.2006.05088.x) showed that multicopy desA partially alleviated the slow growth phenotype of ΔfabA, the difference in which was that ΔfabA was viable. Taken together, our results demonstrate that fabA is essential for aerobic growth. We propose that the plasmid-based ts-allele is useful for exploring the genetic suppression interaction of essential genes of interest in P. aeruginosa. IMPORTANCE Pseudomonas aeruginosa is an opportunistic pathogen whose multidrug resistance demands new drug development. Fatty acids are essential for viability, and essential genes are ideal drug targets. However, the growth defect of essential gene mutants can be suppressed. Suppressors tend to be accumulated during the construction of essential gene deletion mutants, hampering the genetic analysis. To circumvent this issue, we constructed a deletion allele of fabA in the presence of a native promoter-controlled complementary copy in the ts-plasmid. In this analysis, we showed that ΔfabA/pTS-fabA failed to grow at a restrictive temperature, supporting its essentiality. Suppressor analysis revealed desA, whose promoter carried a SNP and whose transcription was upregulated. We validated that both the SNP-bearing promoter-controlled and regulable PBAD promoter-controlled desA suppressed the lethality of ΔfabA. Together, our results demonstrate that fabA is essential for aerobic growth. We propose that plasmid-based ts-alleles are suitable for genetic analysis of essential genes of interest.

Mutations identified in engineered Escherichia coli with a reduced genome

Characterizing genes that regulate cell growth and survival in model organisms is important for understanding higher organisms. Construction of strains harboring large deletions in the genome can provide insights into the genetic basis of cell growth compared with only studying wild-type strains. We have constructed a series of genome-reduced strains with deletions spanning approximately 38.9% of the E. coli chromosome. Strains were constructed by combining large deletions in chromosomal regions encoding nonessential gene groups. We also isolated strains Δ33b and Δ37c, whose growth was partially restored by adaptive laboratory evolution (ALE). Genome sequencing of nine strains, including those selected following ALE, identified the presence of several Single Nucleotide Variants (SNVs), insertions, deletions, and inversions. In addition to multiple SNVs, two insertions were identified in ALE strain Δ33b. The first was an insertion at the promoter region of pntA, which increased cognate gene expression. The second was an insertion sequence (IS) present in sibE, encoding the antitoxin in a toxin-antitoxin system, which decreased expression of sibE. 5 strains of Δ37c independently isolated following ALE harboring multiple SNVs and genetic rearrangements. Interestingly, a SNV was identified in the promoter region of hcaT in all five strains, which increased hcaT expression and, we predict, rescued the attenuated Δ37b growth. Experiments using defined deletion mutants suggested that hcaT encodes a 3-phenylpropionate transporter protein and is involved in survival during stationary phase under oxidative stress. This study is the first to document accumulation of mutations during construction of genome-reduced strains. Furthermore, isolation and analysis of strains derived from ALE in which the growth defect mediated by large chromosomal deletions was rescued identified novel genes involved in cell survival.

Inference of transcriptome signatures of Escherichia coli in long-term stationary phase

"Non-growing" is a dominant life form of microorganisms in nature, where available nutrients and resources are limited. In laboratory culture systems, Escherichia coli can survive for years under starvation, denoted as long-term stationary phase, where a small fraction of cells manages to survive by recycling resources released from nonviable cells. Although the physiology by which viable cells in long-term stationary phase adapt to prolonged starvation is of great interest, their genome-wide response has not been fully understood. In this study, we analyzed transcriptional profiles of cells exposed to the supernatant of 30-day long-term stationary phase culture and found that their transcriptome profiles displayed several similar responses to those of cells in the 16-h short-term stationary phase. Nevertheless, our results revealed that cells in long-term stationary phase supernatant exhibit higher expressions of stress-response genes such as phage shock proteins (psp), and lower expressions of growth-related genes such as ribosomal proteins than those in the short-term stationary phase. We confirmed that the mutant lacking the psp operon showed lower survival and growth rate in the long-term stationary phase culture. This study identified transcriptional responses for stress-resistant physiology in the long-term stationary phase environment.

Application of TraDIS to define the core essential genome of Campylobacter jejuni and Campylobacter coli

Campylobacter species are the major cause of bacterial gastroenteritis. As there is no effective vaccine, combined with the rapid increase in antimicrobial resistant strains, there is a need to identify new targets for intervention. Essential genes are those that are necessary for growth and/or survival, making these attractive targets. In this study, comprehensive transposon mutant libraries were created in six C. jejuni strains, four C. coli strains and one C. lari and C. hyointestinalis strain, allowing for those genes that cannot tolerate a transposon insertion being called as essential. Comparison of essential gene lists using core genome analysis can highlight those genes which are common across multiple strains and/or species. Comparison of C. jejuni and C. coli, the two species that cause the most disease, identified 316 essential genes. Genes of interest highlighted members of the purine pathway being essential for C. jejuni whilst also finding that a functional potassium uptake system is essential. Protein-protein interaction networks using these essential gene lists also highlighted proteins in the purine pathway being major 'hub' proteins which have a large number of interactors across the network. When adding in two more species (C. lari and C. hyointestinalis) the essential gene list reduces to 261. Within these 261 essential genes, there are many genes that have been found to be essential in other bacteria. These include htrB and PEB4, which have previously been found as core virulence genes across Campylobacter species in other studies. There were 21 genes which have no known function with eight of these being associated with the membrane. These surface-associated essential genes may provide attractive targets. The essential gene lists presented will help to prioritise targets for the development of novel therapeutic and preventative interventions.

Whole-genome sequencing analysis of two heat-evolved Escherichia coli strains

BACKGROUND

High temperatures cause a suite of problems for cells, including protein unfolding and aggregation; increased membrane fluidity; and changes in DNA supercoiling, RNA stability, transcription and translation. Consequently, enhanced thermotolerance can evolve through an unknown number of genetic mechanisms even in the simple model bacterium Escherichia coli. To date, each E. coli study exploring this question resulted in a different set of mutations. To understand the changes that can arise when an organism evolves to grow at higher temperatures, we sequenced and analyzed two previously described E. coli strains, BM28 and BM28 ΔlysU, that have been laboratory adapted to the highest E. coli growth temperature reported to date.

RESULTS

We found three large deletions in the BM28 and BM28 ΔlysU strains of 123, 15 and 8.5 kb in length and an expansion of IS10 elements. We found that BM28 and BM28 ΔlysU have considerably different genomes, suggesting that the BM28 culture that gave rise to BM28 and BM28 ΔlysU was a mixed population of genetically different cells. Consistent with published findings of high GroESL expression in BM28, we found that BM28 inexplicitly carries the groESL bearing plasmid pOF39 that was maintained simply by high-temperature selection pressure. We identified over 200 smaller insertions, deletions, single nucleotide polymorphisms and other mutations, including changes in master regulators such as the RNA polymerase and the transcriptional termination factor Rho. Importantly, this genome analysis demonstrates that the commonly cited findings that LysU plays a crucial role in thermotolerance and that GroESL hyper-expression is brought about by chromosomal mutations are based on a previous misinterpretation of the genotype of BM28.

CONCLUSIONS

This whole-genome sequencing study describes genetically distinct mechanisms of thermotolerance evolution from those found in other heat-evolved E. coli strains. Studying adaptive laboratory evolution to heat in simple model organisms is important in the context of climate change. It is important to better understand genetic mechanisms of enhancing thermotolerance in bacteria and other organisms, both in terms of optimizing laboratory evolution methods for various organisms and in terms of potential genetic engineering of organisms most at risk or most important to our societies and ecosystems.

A swapped genetic code prevents viral infections and gene transfer

Engineering the genetic code of an organism has been proposed to provide a firewall from natural ecosystems by preventing viral infections and gene transfer1-6. However, numerous viruses and mobile genetic elements encode parts of the translational apparatus7-9, potentially rendering a genetic-code-based firewall ineffective. Here we show that such mobile transfer RNAs (tRNAs) enable gene transfer and allow viral replication in Escherichia coli despite the genome-wide removal of 3 of the 64 codons and the previously essential cognate tRNA and release factor genes. We then establish a genetic firewall by discovering viral tRNAs that provide exceptionally efficient codon reassignment allowing us to develop cells bearing an amino acid-swapped genetic code that reassigns two of the six serine codons to leucine during translation. This amino acid-swapped genetic code renders cells resistant to viral infections by mistranslating viral proteomes and prevents the escape of synthetic genetic information by engineered reliance on serine codons to produce leucine-requiring proteins. As these cells may have a selective advantage over wild organisms due to virus resistance, we also repurpose a third codon to biocontain this virus-resistant host through dependence on an amino acid not found in nature10. Our results may provide the basis for a general strategy to make any organism safely resistant to all natural viruses and prevent genetic information flow into and out of genetically modified organisms.

About the dark corners in the gene function space of Escherichia coli remaining without illumination by scientific literature

BACKGROUND

Although Escherichia coli (E. coli) is the most studied prokaryote organism in the history of life sciences, many molecular mechanisms and gene functions encoded in its genome remain to be discovered. This work aims at quantifying the illumination of the E. coli gene function space by the scientific literature and how close we are towards the goal of a complete list of E. coli gene functions.

RESULTS

The scientific literature about E. coli protein-coding genes has been mapped onto the genome via the mentioning of names for genomic regions in scientific articles both for the case of the strain K-12 MG1655 as well as for the 95%-threshold softcore genome of 1324 E. coli strains with known complete genome. The article match was quantified with the ratio of a given gene name's occurrence to the mentioning of any gene names in the paper. The various genome regions have an extremely uneven literature coverage. A group of elite genes with ≥ 100 full publication equivalents (FPEs, FPE = 1 is an idealized publication devoted to just a single gene) attracts the lion share of the papers. For K-12, ~ 65% of the literature covers just 342 elite genes; for the softcore genome, ~ 68% of the FPEs is about only 342 elite gene families (GFs). We also find that most genes/GFs have at least one mentioning in a dedicated scientific article (with the exception of at least 137 protein-coding transcripts for K-12 and 26 GFs from the softcore genome). Whereas the literature growth rates were highest for uncharacterized or understudied genes until 2005-2010 compared with other groups of genes, they became negative thereafter. At the same time, literature for anyhow well-studied genes started to grow explosively with threshold T10 (≥ 10 FPEs). Typically, a body of ~ 20 actual articles generated over ~ 15 years of research effort was necessary to reach T10. Lineage-specific co-occurrence analysis of genes belonging to the accessory genome of E. coli together with genomic co-localization and sequence-analytic exploration hints previously completely uncharacterized genes yahV and yddL being associated with osmotic stress response/motility mechanisms.

CONCLUSION

If the numbers of scientific articles about uncharacterized and understudied genes remain at least at present levels, full gene function lists for the strain K-12 MG1655 and the E. coli softcore genome are in reach within the next 25-30 years. Once the literature body for a gene crosses 10 FPEs, most of the critical fundamental research risk appears overcome and steady incremental research becomes possible.

Revealing Causes for False-Positive and False-Negative Calling of Gene Essentiality in Escherichia coli Using Transposon Insertion Sequencing

The massive sequencing of transposon insertion mutant libraries (Tn-Seq) represents a commonly used method to determine essential genes in bacteria. Using a hypersaturated transposon mutant library consisting of 400,096 unique Tn insertions, 523 genes were classified as essential in Escherichia coli K-12 MG1655. This provided a useful genome-wide gene essentiality landscape for rapidly identifying 233 of 301 essential genes previously validated by a knockout study. However, there was a discrepancy in essential gene sets determined by conventional gene deletion methods and Tn-Seq, although different Tn-Seq studies reported different extents of discrepancy. We have elucidated two causes of this discrepancy. First, 68 essential genes not detected by Tn-Seq contain nonessential subgenic domains that are tolerant to transposon insertion, which leads to the false assignment of an essential gene as a nonessential or dispensable gene. These genes exhibited a high level of transposon insertion in their subgenic nonessential domains. In contrast, 290 genes were additionally categorized as essential by Tn-Seq, although their knockout mutants were available. The comparative analysis of Tn-Seq and high-resolution footprinting of nucleoid-associated proteins (NAPs) revealed that a protein-DNA interaction hinders transposon insertion. We identified 213 false-positive genes caused by NAP-genome interactions. These two limitations have to be considered when addressing essential bacterial genes using Tn-Seq. Furthermore, a comparative analysis of high-resolution Tn-Seq with other data sets is required for a more accurate determination of essential genes in bacteria. IMPORTANCE Transposon mutagenesis is an efficient way to explore gene essentiality of a bacterial genome. However, there was a discrepancy between the essential gene set determined by transposon mutagenesis and that determined using single-gene knockout strains. In this study, we generated a hypersaturated Escherichia coli transposon mutant library comprising approximately 400,000 different mutants. Determination of transposon insertion sites using next-generation sequencing provided a high-resolution essentiality landscape of the E. coli genome. We identified false negatives of essential gene discovery due to the permissive insertion of transposons in the C-terminal region. Comparisons between the transposon insertion landscape with binding profiles of DNA-binding proteins revealed interference of nucleoid-associated proteins to transposon insertion, generating false positives of essential gene discovery. Consideration of these findings is required to avoid the misinterpretation of transposon mutagenesis results.

Role of Stress-Induced Proteins RpoS and YicC in the Persistence of Salmonella enterica subsp. enterica Serotype Typhimurium in Tomato Plants

Understanding the functional role of bacterial genes in the persistence of Salmonella in plant organs can facilitate the development of agricultural practices to mitigate food safety risks associated with the consumption of fresh produce contaminated with Salmonella spp. Our study showed that Salmonella enterica subsp. enterica serotype Typhimurium (strain MDD14) persisted less in inoculated tomato plants than other Salmonella Typhimurium strains tested (JSG210, JSG626, JSG634, JSG637, JSG3444, and EV030415; P < 0.01). In-vitro assays performed in limited-nutrient conditions (growth rate, biofilm production, and motility) were inconclusive in explaining the in-planta phenotype observed with MDD14. Whole-genome sequencing combined with non-synonymous single nucleotide variations analysis was performed to identify genomic differences between MDD14 and the other Salmonella Typhimurium strains. The genome of MDD14 contained a truncated version (123 bp N-terminal) of yicC and a mutated version of rpoS (two non-synonymous substitutions, i.e., G66E and R82C), which are two stress-induced proteins involved in iron acquisition, environmental sensing, and cell envelope integrity. The rpoS and yicC genes were deleted in Salmonella Typhimurium JSG210 with the Lambda Red recombining system. Both mutants had limited persistence in tomato plant organs, similar to that of MDD14. In conclusion, we demonstrated that YicC and RpoS are involved in the persistence of Salmonella in tomato plants in greenhouse conditions and, thus, could represent potential targets to mitigate persistence of Salmonella spp. in planta. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Dynamics of high hydrostatic pressure resistance development in RpoS-deficient Escherichia coli

High hydrostatic pressure (HHP) treatment is one of the most widely accepted non-thermal food processing methods, but HHP-resistance development in pathogenic or spoilage bacteria might compromise the safety and stability of HHP-treated foods. Charting the possible routes and mechanisms of HHP resistance development in foodborne bacteria is therefore essential to anticipate or prevent the appearance of resistant variants. While upregulation of the RpoS-governed general stress response is a well-established route for increased HHP resistance in Escherichia coli, previous work revealed that mutations causing attenuated cAMP/CRP activity or aggregation-prone TnaA variants can evolve to overcome the HHP-hypersensitivity of an E. coli ΔrpoS mutant. In this study, further directed evolution and genetic analysis approaches allowed us to demonstrate that both kinds of mutants tend to co-emerge and compete with each other in E. coli ΔrpoS populations evolving towards HHP resistance, because of the higher HHP resistance of cAMP/CRP mutants and the faster growth rate of the TnaA mutants. Moreover, closer scrutiny of evolving populations revealed RpoS, cAMP/CRP and TnaA independent routes of HHP resistance development, based on downregulation of YegW or RppH activity.