Genome-wide identification of transcription start sites, promoters and transcription factor binding sites in E. coli

Lessons from RatA: Why the Basics in Molecular Biology Are Still Crucial!

Since the first bacterial genomes were sequenced and annotated over 25 years ago, sequencing technologies have rapidly advanced in both speed and cost efficiency. To date, over two million annotated bacterial genomes have been deposited in the National Center for Biotechnology Information (NCBI) database. Yet, there are many genes with unknown functions and, furthermore, conflicting results have been published for many investigated genes. One example is the ratA (or pasT) gene from Escherichia coli (E. coli) K-12 strains. Initially identified as a ribosome-targeting toxin, later studies described RatA as the bacterial homolog of the mitochondrial Coq10 protein and, therefore, beneficial for E. coli cells during aerobic growth. This study shows that these conflicting results originated from a mis-annotation of the start codon in the genomic sequence. Overexpression of the ratA gene as currently annotated leads to the synthesis of two RatA protein variants, a toxic and a non-toxic one. This study further identifies the endogenous ratA promoter and shows that only the shorter, non-toxic variant of RatA is synthesized during different growth phases specifically under aerobic conditions. Our findings thereby not only solidify the role of RatA in E. coli, but also demonstrate the importance of first validating the basics of molecular biology when investigating a previously poorly described gene, even in times of advanced high-throughput techniques.

The emergence and evolution of gene expression in genome regions replete with regulatory motifs

Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called -10 and -35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 'promoter islands', DNA sequences enriched with -10 and -35 boxes. We mutagenize these starting 'parent' sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new -10 and -35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all -10 and -35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new -10 and -35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that -10 and -35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.

Deciphering transcript architectural complexity in bacteria and archaea

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Synonymous codon substitutions modulate transcription and translation of a divergent upstream gene by modulating antisense RNA production

Synonymous codons were originally viewed as interchangeable, with no phenotypic consequences. However, substantial evidence has now demonstrated that synonymous substitutions can perturb a variety of gene expression and protein homeostasis mechanisms, including translational efficiency, translational fidelity, and cotranslational folding of the encoded protein. To date, most studies of synonymous codon-derived perturbations have focused on effects within a single gene. Here, we show that synonymous codon substitutions made far within the coding sequence of Escherichia coli plasmid-encoded chloramphenicol acetyltransferase (cat) can significantly increase expression of the divergent upstream tetracycline resistance gene, tetR. In four out of nine synonymously recoded cat sequences tested, expression of the upstream tetR gene was significantly elevated due to transcription of a long antisense RNA (asRNA) originating from a transcription start site within cat. Surprisingly, transcription of this asRNA readily bypassed the native tet transcriptional repression mechanism. Even more surprisingly, accumulation of the TetR protein correlated with the level of asRNA, rather than total tetR RNA. These effects of synonymous codon substitutions on transcription and translation of a neighboring gene suggest that synonymous codon usage in bacteria may be under selection to both preserve the amino acid sequence of the encoded gene and avoid DNA sequence elements that can significantly perturb expression of neighboring genes. Avoiding such sequences may be especially important in plasmids and prokaryotic genomes, where genes and regulatory elements are often densely packed. Similar considerations may apply to the design of genetic circuits for synthetic biology applications.

Escherichia coli non-coding regulatory regions are highly conserved

Microbial genome sequences are rapidly accumulating, enabling large-scale studies of sequence variation. Existing studies primarily focus on coding regions to study amino acid substitution patterns in proteins. However, non-coding regulatory regions also play a distinct role in determining physiologic responses. To investigate intergenic sequence variation on a large-scale, we identified non-coding regulatory region alleles across 2350 Escherichia coli strains. This 'alleleome' consists of 117 781 unique alleles for 1169 reference regulatory regions (transcribing 1975 genes) at single base-pair resolution. We find that 64% of nucleotide positions are invariant, and variant positions vary in a median of just 0.6% of strains. Additionally, non-coding alleles are sufficient to recover E. coli phylogroups. We find that core promoter elements and transcription factor binding sites are significantly conserved, especially those located upstream of essential or highly-expressed genes. However, variability in conservation of transcription factor binding sites is significant both within and across regulons. Finally, we contrast mutations acquired during adaptive laboratory evolution with wild-type variation, finding that the former preferentially alter positions that the latter conserves. Overall, this analysis elucidates the wealth of information found in E. coli non-coding sequence variation and expands pangenomic studies to non-coding regulatory regions at single-nucleotide resolution.

Construction and characterization of a promoter library with varying strengths to enhance acetoin production from xylose in Serratia marcescens

Serratia marcescens is utilized as a significant enterobacteria in the production of various high-value secondary metabolites. Acetoin serves as a crucial foundational compound of development and finds application in a broad range of fields. Furthermore, S. marcescens HBQA-7 is capable of utilizing xylose as its exclusive carbon source for acetoin production. The objective of this study was to utilize a constitutive promoter screening strategy to enhance both xylose utilization and acetoin production in S. marcescens HBQA-7. By utilizing RNA-seq, we identified the endogenous constitutive promoter P6 that is the most robust, which facilitated the overexpression of the sugar transporter protein GlfL445I, α-acetyl lactate synthase, and α-acetyl lactate decarboxylase, respectively. The resultant recombinant strains exhibited enhanced xylose utilization rates and acetoin yields. Subsequently, a recombinant plasmid, denoted as pBBR1MCS-P6-glfL445IalsSalsD, was constructed, simultaneously expressing the aforementioned three genes. The resulting recombinant strain, designated as S3, demonstrated a 1.89-fold boost in xylose consumption rate compared with the original strain during shake flask fermentation. resulting in the accumulation of 7.14 g/L acetoin in the final fermentation medium. Subsequently, in a 5 L fermenter setup, the acetoin yield reached 48.75 g/L, corresponding to a xylose-to-acetoin conversion yield of 0.375 g/g.

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

New layers of regulation of the general stress response sigma factor RpoS

The general stress response (GSR) sigma factor RpoS from Escherichia coli has emerged as one of the key paradigms for study of how numerous signal inputs are accepted at multiple levels into a single pathway for regulation of gene expression output. While many studies have elucidated the key pathways controlling the production and activity of this sigma factor, recent discoveries have uncovered still more regulatory mechanisms which feed into the network. Moreover, while the regulon of this sigma factor comprises a large proportion of the E. coli genome, the downstream expression levels of all the RpoS target genes are not identically affected by RpoS upregulation but respond heterogeneously, both within and between cells. This minireview highlights the most recent developments in our understanding of RpoS regulation and expression, in particular those which influence the regulatory network at different levels from previously well-studied pathways.

Balanced cell division is secured by two different regulatory sites in OxyS RNA

The hydrogen peroxide-induced small RNA OxyS has been proposed to originate from the 3' UTR of a peroxide mRNA. Unexpectedly, phylogenetic OxyS targetome predictions indicate that most OxyS targets belong to the category of "cell cycle," including cell division and cell elongation. Previously, we reported that Escherichia coli OxyS inhibits cell division by repressing expression of the essential transcription termination factor nusG, thereby leading to the expression of the KilR protein, which interferes with the function of the major cell division protein, FtsZ. By interfering with cell division, OxyS brings about cell-cycle arrest, thus allowing DNA damage repair. Cell division and cell elongation are opposing functions to the extent that inhibition of cell division requires a parallel inhibition of cell elongation for the cells to survive. In this study, we report that in addition to cell division, OxyS inhibits mepS, which encodes an essential peptidoglycan endopeptidase that is responsible for cell elongation. Our study indicates that cell-cycle arrest and balancing between cell division and cell elongation are important and conserved functions of the oxidative stress-induced sRNA OxyS.

Promoter capture drives the emergence of proto-genes in Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli Long-Term Evolution Experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time-span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, thereby serving as raw material for new gene emergence. Most proto-genes result either from insertion element activity or chromosomal translocations that fused pre-existing regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, persist stably, and can serve as potential substrates for new gene formation.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Antisense transcription and its roles in adaption to environmental stress in E. coli

It has been reported that a highly varying proportion (1% ∼ 93%) of genes in various prokaryotes have antisense RNA (asRNA) transcription. However, the extent of the pervasiveness of asRNA transcription in the well-studied E. coli K12 strain has thus far been an issue of debate. Furthermore, very little is known about the expression patterns and functions of asRNAs under various conditions. To fill these gaps, we determined the transcriptomes and proteomes of E. coli K12 at multiple time points in five culture conditions using strand-specific RNA-seq, differential RNA-seq, and quantitative mass spectrometry methods. To reduce artifacts of possible transcriptional noise, we identified asRNA using stringent criteria with biological replicate verification and transcription start sites (TSSs) information included. We identified a total of 660 asRNAs, which were generally short and largely condition-dependently transcribed. We found that the proportions of the genes which had asRNA transcription highly depended on the culture conditions and time points. We classified the transcriptional activities of the genes in six transcriptional modes according to their relative levels of asRNA to mRNA. Many genes changed their transcriptional modes at different time points of the culture conditions, and such transitions can be described in a well-defined manner. Intriguingly, the protein levels and mRNA levels of genes in the sense-only/sense-dominant mode were moderately correlated, but the same was not true for genes in the balanced/antisense-dominant mode, in which asRNAs were at a comparable or higher level to mRNAs. These observations were further validated by western blot on candidate genes, where an increase in asRNA transcription diminished gene expression in one case and enhanced it in another. These results suggest that asRNAs may directly or indirectly regulate translation by forming duplexes with cognate mRNAs. Thus, asRNAs may play an important role in the bacterium's responses to environmental changes during growth and adaption to different environments.

IMPORTANCE

The cis -antisense RNA (asRNA) is a type of understudied RNA molecules in prokaryotes, which is believed to be important in regulating gene expression. Our current understanding of asRNA is constrained by inconsistent reports about its identification and properties. These discrepancies are partially caused by a lack of sufficient samples, biological replicates, and culture conditions. This study aimed to overcome these disadvantages and identified 660 putative asRNAs using integrated information from strand-specific RNA-seq, differential RNA-seq, and mass spectrometry methods. In addition, we explored the relative expression between asRNAs and sense RNAs and investigated asRNA regulated transcriptional activity changes over different culture conditions and time points. Our work strongly suggests that asRNAs may play a crucial role in bacterium's responses to environmental changes during growth and adaption to different environments.

CvkR is a MerR-type transcriptional repressor of class 2 type V-K CRISPR-associated transposase systems

Certain CRISPR-Cas elements integrate into Tn7-like transposons, forming CRISPR-associated transposon (CAST) systems. How the activity of these systems is controlled in situ has remained largely unknown. Here we characterize the MerR-type transcriptional regulator Alr3614 that is encoded by one of the CAST (AnCAST) system genes in the genome of cyanobacterium Anabaena sp. PCC 7120. We identify a number of Alr3614 homologs across cyanobacteria and suggest naming these regulators CvkR for Cas V-K repressors. Alr3614/CvkR is translated from leaderless mRNA and represses the AnCAST core modules cas12k and tnsB directly, and indirectly the abundance of the tracr-CRISPR RNA. We identify a widely conserved CvkR binding motif 5'-AnnACATnATGTnnT-3'. Crystal structure of CvkR at 1.6 Å resolution reveals that it comprises distinct dimerization and potential effector-binding domains and that it assembles into a homodimer, representing a discrete structural subfamily of MerR regulators. CvkR repressors are at the core of a widely conserved regulatory mechanism that controls type V-K CAST systems.

An RpoN-dependent PEP-CTERM gene is involved in floc formation of an Aquincola tertiaricarbonis strain

BACKGROUND

The floc is a characteristic of microbial aggregate growth, displaying cloudy suspensions in water. Floc formation has been demonstrated in a series of bacteria and the floc-forming bacteria play a crucial role in activated sludge (AS) process widely used for municipal sewage and industrial wastewater treatment over a century. It has been demonstrated that some exopolysaccharide biosynthesis genes and the sigma factor (sigma54 or rpoN) were required for floc forming in some bacteria. However, the mechanism underlying the floc formation stills need to be elucidated.

RESULTS

In this study, we demonstrate that a TPR (Tetratricopeptide repeats) protein-encoding gene prsT is required for floc formation of Aquincola tertiaricarbonis RN12 and an upstream PEP-CTERM gene (designated pepA), regulated by RpoN1, is involved in its floc formation but not swarming motility and biofilm formation. Overexpression of PepA could rescue the floc-forming phenotype of the rpoN1 mutant by decreasing the released soluble exopolysaccharides and increasing the bound polymers.

CONCLUSION

Our results indicate that the wide-spread PEP-CTERM proteins play an important role in the self-flocculation of bacterial cells and may be a component of extracellular polymeric substances required for floc-formation.

Systems engineering of Escherichia coli for high-level shikimate production

To further increase the production efficiency of microbial shikimate, a valuable compound widely used in the pharmaceutical and chemical industries, ten key target genes contributing to shikimate production were identified by exploiting the enzyme constraint model ec_iML1515, and subsequently used for promoting metabolic flux towards shikimate biosynthesis in the tryptophan-overproducing strain Escherichia coli TRP0. The engineered E. coli SA05 produced 78.4 g/L shikimate via fed-batch fermentation. Deletion of quinate dehydrogenase and introduction of the hydroaromatic equilibration-alleviating shikimate dehydrogenase mutant AroE^T61W/L241I reduced the contents of byproducts quinate (7.5 g/L) and 3-dehydroshikimic acid (21.4 g/L) by 89.1% and 52.1%, respectively. Furthermore, a high concentration shikimate responsive promoter P_rpoS was recruited to dynamically regulate the expression of the tolerance target ProV to enhance shikimate productivity by 23.2% (to 2 g/L/h). Finally, the shikimate titer was increased to 126.4 g/L, with a yield of 0.50 g/g glucose and productivity of 2.63 g/L/h, using a 30-L fermenter and the engineered strain E. coli SA09. This is, to the best of our knowledge, the highest reported shikimate titer and productivity in E. coli.

Zur and zinc increase expression of E. coli ribosomal protein L31 through RNA-mediated repression of the repressor L31p

Bacteria can adapt in response to numerous stress conditions. One such stress condition is zinc depletion. The zinc-sensing transcription factor Zur regulates the way numerous bacterial species respond to severe changes in zinc availability. Under zinc sufficient conditions, Zn-loaded Zur (Zn2-Zur) is well-known to repress transcription of genes encoding zinc uptake transporters and paralogues of a few ribosomal proteins. Here, we report the discovery and mechanistic basis for the ability of Zur to up-regulate expression of the ribosomal protein L31 in response to zinc in E. coli. Through genetic mutations and reporter gene assays, we find that Zur achieves the up-regulation of L31 through a double repression cascade by which Zur first represses the transcription of L31p, a zinc-lacking paralogue of L31, which in turn represses the translation of L31. Mutational analyses show that translational repression by L31p requires an RNA hairpin structure within the l31 mRNA and involves the N-terminus of the L31p protein. This work uncovers a new genetic network that allows bacteria to respond to host-induced nutrient limiting conditions through a sophisticated ribosomal protein switching mechanism.

Regulation of constant cell elongation and Sfm pili synthesis in Escherichia coli via two active forms of FimZ orphan response regulator

Escherichia coli (E. coli) has multiple copies of the chaperone-usher (CU) pili operon in five fimbria groups: CU pili, curli, type IV pili, type III secretion pili, and type IV secretion pili. Commensal E. coli K-12 contains 12 CU pili operons. Among these operons, Sfm is expressed by the sfmACDHF operon. Transcriptome analyses, reporter assays, and chromatin immunoprecipitation PCR analyses reported that FimZ directly binds to and activates the sfmA promoter, transcribing sfmACDHF. In addition, FimZ regularly induces constant cell elongation in E. coli, which is required for F-type ATPase function. The bacterial two-hybrid system showed a specific interaction between FimZ and the α subunit of the cytoplasmic F1 domain of F-type ATPase. Studies performed using mutated FimZs have revealed two active forms, I and II. Active form I is required for constant cell elongation involving amino acid residues K106 and D109. Active form II additionally required D56, a putative phosphorylation site, to activate the sfmA promoter. The chromosomal fimZ was hardly expressed in parent strain but functioned in phoB and phoP double-gene knockout strains. These insights may help to understand bacterial invasion restricted host environments by the sfm γ-type pili.

The regulon of Brucella abortus two-component system BvrR/BvrS reveals the coordination of metabolic pathways required for intracellular life

Brucella abortus is a facultative intracellular pathogen causing a severe zoonotic disease worldwide. The two-component regulatory system (TCS) BvrR/BvrS of B. abortus is conserved in members of the Alphaproteobacteria class. It is related to the expression of genes required for host interaction and intracellular survival. Here we report that bvrR and bvrS are part of an operon composed of 16 genes encoding functions related to nitrogen metabolism, DNA repair and recombination, cell cycle arrest, and stress response. Synteny of this genomic region within close Alphaproteobacteria members suggests a conserved role in coordinating the expression of carbon and nitrogen metabolic pathways. In addition, we performed a ChIP-Seq analysis after exposure of bacteria to conditions that mimic the intracellular environment. Genes encoding enzymes at metabolic crossroads of the pentose phosphate shunt, gluconeogenesis, cell envelope homeostasis, nucleotide synthesis, cell division, and virulence are BvrR/BvrS direct targets. A 14 bp DNA BvrR binding motif was found and investigated in selected gene targets such as virB1, bvrR, pckA, omp25, and tamA. Understanding gene expression regulation is essential to elucidate how Brucella orchestrates a physiological response leading to a furtive pathogenic strategy.

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

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

Genome-wide transcription start site mapping in the facultative intracellular pathogen Brucella melitensis by Capping-seq

Brucella melitensis (B. melitensis) is an important facultative intracellular bacterium that causes global zoonotic diseases. Continuous intracellular survival and replication are the main obstruction responsible for the accessibility of prevention and treatment of brucellosis. Bacteria respond to complex environment by regulating gene expression. Many regulatory factors function at loci where RNA polymerase initiates messenger RNA synthesis. However, limited gene annotation is a current obstacle for the research on expression regulation in bacteria. To improve annotation and explore potential functional sites, we proposed a novel genome-wide method called Capping-seq for transcription start site (TSS) mapping in B. melitensis. This technique combines capture of capped primary transcripts with Single Molecule Real-Time (SMRT) sequencing technology. We identified 2,369 TSSs at single nucleotide resolution by Capping-seq. TSSs analysis of Brucella transcripts showed a preference of purine on the TSS positions. Our results revealed that -35 and -10 elements of promoter contained consensus sequences of TTGNNN and TATNNN, respectively. The 5' ends analysis showed that 57% genes are associated with more than one TSS and 47% genes contain long leader regions, suggested potential complex regulation at the 5' ends of genes in B. melitensis. Moreover, we identified 52 leaderless genes that are mainly involved in the metabolic processes. Overall, Capping-seq technology provides a unique solution for TSS determination in prokaryotes. Our findings develop a systematic insight into the primary transcriptome characterization of B. melitensis. This study represents a critical basis for investigating gene regulation and pathogenesis of Brucella.

In Silico Study of Mercury Resistance Genes Extracted from Pseudomonas spp. Involved in Bioremediation: Understanding the Promoter Regions and Regulatory Elements

Microbial genes and their product were diverse and beneficial for heavy metal bioremediation from the contaminated sites. Screening of genes and gene products plays a significant role in the detoxification of pollutants. Understanding of the promoter region and its regulatory elements is a vital implication of microbial genes. To the best of our knowledge, there is no in silico study reported so far on mer gene families used for heavy metal bioremediation. The motif distribution was observed densely upstream of the TSSs (transcription start sites) between +1 and -350 bp and sparsely distributed above -350 bp, according to the current study. MEME identified the best common candidate motifs of TFs (transcription factors) binding with the lowest e value (7.2 $e$ -033) and is the most statistically significant candidate motif. The EXPREG output of the 11 TFs with varying degrees of function such as activation, repression, transcription, and dual purposes was thoroughly examined. Data revealed that transcriptional gene regulation in terms of activation and repression was observed at 36.4% and 54.56%, respectively. This shows that most TFs are involved in transcription gene repression rather than activation. Likewise, EXPREG output revealed that transcriptional conformational modes, such as monomers, dimers, tetramers, and other factors, were also analyzed. The data indicated that most of the transcriptional conformation mode was dual, which accounts for 96%. CpG island analysis using online and offline tools revealed that the gene body had fewer CpG islands compared to the promoter regions. Understanding the common candidate motifs, transcriptional factors, and regulatory elements of the mer operon gene cluster using a machine learning approach could help us better understand gene expression patterns in heavy metal bioremediation.

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

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

Genome-wide promoter assembly in E. coli measured at single-base resolution

When detected at single-base-pair resolution, the genome-wide location, occupancy level, and structural organization of DNA-binding proteins provide mechanistic insights into genome regulation. Here we use ChIP-exo to provide a near-base-pair resolution view of the epigenomic organization of the Escherichia coli transcription machinery and nucleoid structural proteins at the time when cells are growing exponentially and upon rapid reprogramming (acute heat shock). We examined the site specificity of three sigma factors (RpoD/σ70, RpoH/σ32, and RpoN/σ54), RNA polymerase (RNAP or RpoA, -B, -C), and two nucleoid proteins (Fis and IHF). We suggest that DNA shape at the flanks of cognate motifs helps drive site specificity. We find that although RNAP and sigma factors occupy active cognate promoters, RpoH and RpoN can occupy quiescent promoters without the presence of RNAP. Thus, promoter-bound sigma factors can be triggered to recruit RNAP by a mechanism that is distinct from an obligatory cycle of free sigma binding RNAP followed by promoter binding. These findings add new dimensions to how sigma factors achieve promoter specificity through DNA sequence and shape, and further define mechanistic steps in regulated genome-wide assembly of RNAP at promoters in E. coli.

The Escherichia coli serS gene promoter region overlaps with the rarA gene

Deletion of the entire gene encoding the RarA protein of Escherichia coli results in a growth defect and additional deficiencies that were initially ascribed to a lack of RarA function. Further work revealed that most of the effects reflected the presence of sequences in the rarA gene that affect expression of the downstream gene, serS. The serS gene encodes the seryl aminoacyl-tRNA synthetase. Decreases in the expression of serS can trigger the stringent response. The sequences that affect serS expression are located in the last 15 nucleotides of the rarA gene.

Mutants lacking global regulators, fis and arcA, in Escherichia coli enhanced growth fitness under acetate metabolism by pathway reprogramming

Global regulatory transcription factors play a significant role in controlling microbial metabolism under genetic and environmental perturbations. A system-level effect of carbon sources such as acetate on microbial metabolism under disrupted global regulators has not been well established. Acetate is one of the major substrates available in various nutrient niches such as the mammalian gut and a keto diet. A substantial amount of acetate gets secreted in aerobic metabolism. Therefore, investigating the study on acetate metabolism is highly significant. It is known that the global regulators fis and arcA regulate acetate uptake genes in E. coli under glucose conditions. This study deciphered the growth and flux distribution of E. coli transcription regulatory knockouts Δfis, ΔarcA and double deletion mutant, ΔarcAΔfis under acetate using ¹³C-metabolic flux analysis (MFA), which has not been investigated before. We observed that the mutants exhibited an expeditious growth rate (~ 1.2-1.6-fold) with a proportionate increase in acetate uptake rates compared to the wild type. ¹³C-MFA displayed the distinct metabolic reprogramming of intracellular fluxes via the TCA cycle, anaplerotic pathway and gluconeogenesis, which conferred an advantage of a faster growth rate with better carbon usage in all the mutants. This resulted in higher metabolic fluxes through the TCA cycle (~ 18-90%), lower gluconeogenesis (~ 15-35%) and higher CO₂ and ATP production with the proportional increase in growth rate. The study reveals a novel insight by stating the sub-optimality of the wild-type strain grown under acetate substrate aerobically. These mutant strains efficiently oxidize acetate, thus acting as potential candidates for the biosynthesis of isoprenoids, biofuels, vitamins and various pharmaceutical products.Key Points• Mutants exhibited a better balance between energy and precursor synthesis than WT.• Leveraged in the unravelling of regulatory control under various nutrient shifts.• Metabolic readjustment resulted in optimal biomass requirement and faster growth.

In silico analysis of promoter regions to identify regulatory elements in TetR family transcriptional regulatory genes of Mycobacterium colombiense CECT 3035

Background

Mycobacterium colombiense is an acid-fast, non-motile, rod-shaped mycobacterium confirmed to cause respiratory disease and disseminated infection in immune-compromised patients, and lymphadenopathy in immune-competent children. It has virulence mechanisms that allow them to adapt, survive, replicate, and produce diseases in the host. To tackle the diseases caused by M. colombiense, understanding of the regulation mechanisms of its genes is important. This paper, therefore, analyzes transcription start sites, promoter regions, motifs, transcription factors, and CpG islands in TetR family transcriptional regulatory (TFTR) genes of M. colombiense CECT 3035 using neural network promoter prediction, MEME, TOMTOM algorithms, and evolutionary analysis with the help of MEGA-X.

Results

The analysis of 22 protein coding TFTR genes of M. colombiense CECT 3035 showed that 86.36% and 13.64% of the gene sequences had one and two TSSs, respectively. Using MEME, we identified five motifs (MTF1, MTF2, MTF3, MTF4, and MTF5) and MTF1 was revealed as the common promoter motif for 100% TFTR genes of M. colombiense CECT 3035 which may serve as binding site for transcription factors that shared a minimum homology of 95.45%. MTF1 was compared to the registered prokaryotic motifs and found to match with 15 of them. MTF1 serves as the binding site mainly for AraC, LexA, and Bacterial histone-like protein families. Other protein families such as MATP, RR, σ-70 factor, TetR, LytTR, LuxR, and NAP also appear to be the binding candidates for MTF1. These families are known to have functions in virulence mechanisms, metabolism, quorum sensing, cell division, and antibiotic resistance. Furthermore, it was found that TFTR genes of M. colombiense CECT 3035 have many CpG islands with several fragments in their CpG islands. Molecular evolutionary genetic analysis showed close relationship among the genes.

Conclusion

We believe these findings will provide a better understanding of the regulation of TFTR genes in M. colombiense CECT 3035 involved in vital processes such as cell division, pathogenesis, and drug resistance and are likely to provide insights for drug development important to tackle the diseases caused by this mycobacterium. We believe this is the first report of in silico analyses of the transcriptional regulation of M. colombiense TFTR genes.

Regulation of Leaderless mRNA Translation in Bacteria

In bacteria, the translation of genetic information can begin through at least three different mechanisms: canonical or Shine-Dalgarno-led initiation, readthrough or 70S scanning initiation, or leaderless initiation. Here, we discuss the main features and regulation of the last, which is characterized mainly by the ability of 70S ribosomal particles to bind to AUG located at or near the 5′ end of mRNAs to initiate translation. These leaderless mRNAs (lmRNAs) are rare in enterobacteria, such as Escherichia coli, but are common in other bacteria, such as Mycobacterium tuberculosis and Deinococcus deserti, where they may represent more than 20% and even up to 60% of the genes. Given that lmRNAs are devoid of a 5′ untranslated region and the Shine-Dalgarno sequence located within it, the mechanism of translation regulation must depend on molecular strategies that are different from what has been observed in the Shine-Dalgarno-led translation. Diverse regulatory mechanisms have been proposed, including the processing of ribosomal RNA and changes in the abundance of translation factors, but all of them produce global changes in the initiation of lmRNA translation. Thus, further research will be required to understand how the initiation of the translation of particular lmRNA genes is regulated.

Tailoring the evolution of BL21(DE3) uncovers a key role for RNA stability in gene expression toxicity

Gene expression toxicity is an important biological phenomenon and a major bottleneck in biotechnology. Escherichia coli BL21(DE3) is the most popular choice for recombinant protein production, and various derivatives have been evolved or engineered to facilitate improved yield and tolerance to toxic genes. However, previous efforts to evolve BL21, such as the Walker strains C41 and C43, resulted only in decreased expression strength of the T7 system. This reveals little about the mechanisms at play and constitutes only marginal progress towards a generally higher producing cell factory. Here, we restrict the solution space for BL21(DE3) to evolve tolerance and isolate a mutant strain Evo21(DE3) with a truncation in the essential RNase E. This suggests that RNA stability plays a central role in gene expression toxicity. The evolved rne truncation is similar to a mutation previously engineered into the commercially available BL21Star(DE3), which challenges the existing assumption that this strain is unsuitable for expressing toxic proteins. We isolated another dominant mutation in a presumed substrate binding site of RNase E that improves protein production further when provided as an auxiliary plasmid. This makes it easy to improve other BL21 variants and points to RNases as prime targets for cell factory optimisation.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Widespread divergent transcription from bacterial and archaeal promoters is a consequence of DNA-sequence symmetry

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Regulation of Expression of the TIR-Containing Protein C Gene of the Uropathogenic Escherichia coli Strain CFT073

The uropathogenic Escherichia coli strain CFT073 causes kidney abscesses in mice Toll/interleukin-1 receptor domain-containing protein C (TcpC) dependently and the corresponding gene is present in around 40% of E. coli isolates of pyelonephritis patients. It impairs the Toll-like receptor (TLR) signaling chain and the NACHT leucin-rich repeat PYD protein 3 inflammasome (NLRP3) by binding to TLR4 and myeloid differentiation factor 88 as well as to NLRP3 and caspase-1, respectively. Overexpression of the tcpC gene stopped replication of CFT073. Overexpression of several tcpC-truncation constructs revealed a transmembrane region, while its TIR domain induced filamentous bacteria. Based on these observations, we hypothesized that tcpC expression is presumably tightly controlled. We tested two putative promoters designated P1 and P2 located at 5′ of the gene c2397 and 5′ of the tcpC gene (c2398), respectively, which may form an operon. High pH and increasing glucose concentrations stimulated a P2 reporter construct that was considerably stronger than a P1 reporter construct, while increasing FeSO₄ concentrations suppressed their activity. Human urine activated P2, demonstrating that tcpC might be induced in the urinary tract of infected patients. We conclude that P2, consisting of a 240 bp region 5′ of the tcpC gene, represents the major regulator of tcpC expression.

Identification of Gradient Promoters of Gluconobacter oxydans and Their Applications in the Biosynthesis of 2-Keto-L-Gulonic Acid

The acetic acid bacterium Gluconobacter oxydans is known for its unique incomplete oxidation and therefore widely applied in the industrial production of many compounds, e.g., 2-keto-L-gulonic acid (2-KLG), the direct precursor of vitamin C. However, few molecular tools are available for metabolically engineering G. oxydans, which greatly limit the strain development. Promoters are one of vital components to control and regulate gene expression at the transcriptional level for boosting production. In this study, the low activity of SDH was found to hamper the high yield of 2-KLG, and enhancing the expression of SDH was achieved by screening the suitable promoters based on RNA sequencing data. We obtained 97 promoters from G. oxydans’s genome, including two strong shuttle promoters and six strongest promoters. Among these promoters, P₃₀₂₂ and P₀₉₄₃ revealed strong activities in both Escherichia coli and G. oxydans, and the activity of the strongest promoter (P₂₇₀₃) was about threefold that of the other reported strong promoters of G. oxydans. These promoters were used to overexpress SDH in G. oxydans WSH-003. The titer of 2-KLG reached 3.7 g/L when SDH was under the control of strong promoters P₂₀₅₇ and P₂₇₀₃. This study obtained a series of gradient promoters, including two strong shuttle promoters, and expanded the toolbox of available promoters for the application in metabolic engineering of G. oxydans for high-value products.

Knowns and Unknowns of Vitamin B6 Metabolism in Escherichia coli

Vitamin B₆ is an ensemble of six interconvertible vitamers: pyridoxine (PN), pyridoxamine (PM), pyridoxal (PL), and their 5'-phosphate derivatives, PNP, PMP, and PLP. Pyridoxal 5'-phosphate is a coenzyme in a variety of enzyme reactions concerning transformations of amino and amino acid compounds. This review summarizes all known and putative PLP-binding proteins found in the Escherichia coli MG1655 proteome. PLP can have toxic effects since it contains a very reactive aldehyde group at its 4' position that easily forms aldimines with primary and secondary amines and reacts with thiols. Most PLP is bound either to the enzymes that use it as a cofactor or to PLP carrier proteins, protected from the cellular environment but at the same time readily transferable to PLP-dependent apoenzymes. E. coli and its relatives synthesize PLP through the seven-step deoxyxylulose-5-phosphate (DXP)-dependent pathway. Other bacteria synthesize PLP in a single step, through a so-called DXP-independent pathway. Although the DXP-dependent pathway was the first to be revealed, the discovery of the widespread DXP-independent pathway determined a decline of interest in E. coli vitamin B₆ metabolism. In E. coli, as in most organisms, PLP can also be obtained from PL, PN, and PM, imported from the environment or recycled from protein turnover, via a salvage pathway. Our review deals with all aspects of vitamin B₆ metabolism in E. coli, from transcriptional to posttranslational regulation. A critical interpretation of results is presented, in particular, concerning the most obscure aspects of PLP homeostasis and delivery to PLP-dependent enzymes.

Rewiring the specificity of extracytoplasmic function sigma factors

Bacterial genomes are being sequenced at an exponentially increasing rate, but our inability to decipher their transcriptional wiring limits our ability to derive new biology from these sequences. De novo determination of regulatory interactions requires accurate prediction of regulators' DNA binding and precise determination of biologically significant binding sites. Here we address these challenges by solving the DNA-specificity code of extracytoplasmic function sigma factors (ECF σs), a major family of bacterial regulators, and determining their putative regulons. We generated an aligned collection of ECF σs and their promoters by leveraging the autoregulatory nature of ECF σs as a means of promoter discovery and analyzed it to identify and characterize the conserved amino acid-nucleotide interactions that determine promoter specificity. This enabled de novo prediction of ECF σ specificity, which we combined with a statistically rigorous phylogenetic footprinting pipeline based on precomputed orthologs to predict the direct targets of ∼67% of ECF σs. This global survey indicated that some ECF σs are conserved global regulators controlling many genes throughout the genome, which are important under many conditions, while others are local regulators, controlling a few closely linked genes in response to specific stimuli in select species. This analysis reveals important organizing principles of bacterial gene regulation and presents a conceptual and computational framework for deciphering gene regulatory networks.

The small regulatory RNA Lpr10 regulates the expression of RpoS in Legionella pneumophila

Legionella pneumophila (Lp) is a waterborne bacterium able to infect human alveolar macrophages, causing Legionnaires' disease. Lp can survive for several months in water, while searching for host cells to grow in, such as ciliates and amoeba. In Lp, the sigma factor RpoS is essential for survival in water. A previous transcriptomic study showed that RpoS positively regulates the small regulatory RNA Lpr10. In the present study, deletion of lpr10 results in an increased survival of Lp in water. Microarray analysis and RT-qPCR revealed that Lpr10 negatively regulates the expression of RpoS in the postexponential phase. Electrophoretic mobility shift assay and in-line probing showed that Lpr10 binds to a region upstream of the previously identified transcription start sites (TSS) of rpoS. A third putative transcription start site was identified by primer extension analysis, upstream of the Lpr10 binding site. In addition, nlpD TSS produces a polycistronic mRNA including the downstream gene rpoS, indicating a fourth TSS for rpoS. Our results suggest that the transcripts from the third and fourth TSS are negatively regulated by the Lpr10 sRNA. Therefore, we propose that Lpr10 is involved in a negative regulatory feedback loop to maintain expression of RpoS to an optimal level.

The BB0345 Hypothetical Protein of Borrelia burgdorferi Is Essential for Mammalian Infection

During the natural enzootic life cycle of Borrelia burgdorferi (also known as Borreliella burgdorferi), the bacteria must sense conditions within the vertebrate and arthropod and appropriately regulate expression of genes necessary to persist within these distinct environments. bb0345 of B. burgdorferi encodes a hypothetical protein of unknown function that is predicted to contain an N-terminal helix-turn-helix (HTH) domain. Because HTH domains can mediate protein-DNA interactions, we hypothesized that BB0345 might represent a previously unidentified borrelial transcriptional regulator with the ability to regulate events critical for the B. burgdorferi enzootic cycle. To study the role of BB0345 within mammals, we generated a bb0345 mutant and assessed its virulence potential in immunocompetent mice. The bb0345 mutant was able to initiate localized infection and disseminate to distal tissues but was cleared from all sites by 14 days postinfection. In vitro growth curve analyses revealed that the bb0345 mutant grew similar to wild-type bacteria in standard Barbour-Stoenner-Kelley II (BSK-II) medium; however, the mutant was not able to grow in dilute BSK-II medium or dialysis membrane chambers (DMCs) implanted in rats. Proteinase K accessibility assays and whole-cell partitioning indicated that BB0345 was intracellular and partially membrane associated. Comparison of protein production profiles between the wild-type parent and the bb0345 mutant revealed no major differences, suggesting BB0345 may not be a global transcriptional regulator. Taken together, these data show that BB0345 is essential for B. burgdorferi survival in the mammalian host, potentially by aiding the spirochete with a physiological function that is required by the bacterium during infection.

Introducing differential RNA-seq mapping to track the early infection phase for Pseudomonas phage ɸKZ

As part of the ongoing renaissance of phage biology, more phage genomes are becoming available through DNA sequencing. However, our understanding of the transcriptome architecture that allows these genomes to be expressed during host infection is generally poor. Transcription start sites (TSSs) and operons have been mapped for very few phages, and an annotated global RNA map of a phage – alone or together with its infected host – is not available at all. Here, we applied differential RNA-seq (dRNA-seq) to study the early, host takeover phase of infection by assessing the transcriptome structure of Pseudomonas aeruginosa jumbo phage ɸKZ, a model phage for viral genetics and structural research. This map substantially expands the number of early expressed viral genes, defining TSSs that are active ten minutes after ɸKZ infection. Simultaneously, we record gene expression changes in the host transcriptome during this critical metabolism conversion. In addition to previously reported upregulation of genes associated with amino acid metabolism, we observe strong activation of genes with functions in biofilm formation (cdrAB) and iron storage (bfrB), as well as an activation of the antitoxin ParD. Conversely, ɸKZ infection rapidly down-regulates complexes IV and V of oxidative phosphorylation (atpCDGHF and cyoABCDE). Taken together, our data provide new insights into the transcriptional organization and infection process of the giant bacteriophage ɸKZ and adds a framework for the genome-wide transcriptomic analysis of phage–host interactions.

Studies on the Regulation of (p)ppGpp Metabolism and Its Perturbation Through the Over-Expression of Nudix Hydrolases in Escherichia coli

Stringent response mediated by modified guanosine nucleotides is conserved across bacteria and is regulated through the Rel/Spo functions. In Escherichia coli, RelA and SpoT proteins synthesize the modified nucleotides ppGpp and pppGpp, together referred to as (p)ppGpp. SpoT is also the primary (p)ppGpp hydrolase. In this study, using hypomorphic relA alleles, we provide experimental evidence for SpoT-mediated negative regulation of the amplification of RelA-dependent stringent response. We investigated the kinetics of ppGpp degradation in cells recovering from stringent response in the complete absence of SpoT function. We found that, although greatly diminished, there was slow ppGpp degradation and growth resumption after a lag period, concomitant with decrease in ppGpp pool. We present evidence for reduction in the ppGpp degradation rate following an increase in pppGpp pool, during recovery from stringent response. From a genetic screen, the nudix hydrolases MutT and NudG were identified as over-expression suppressors of the growth defect of ΔspoT and ΔspoT ΔgppA strains. The effect of over-expression of these hydrolases on the stringent response to amino acid starvation and basal (p)ppGpp pool was studied. Over-expression of each hydrolase reduced the strength of the stringent response to amino acid starvation, and additionally, perturbed the ratio of ppGpp to pppGpp in strains with reduced SpoT hydrolase activity. In these strains that do not accumulate pppGpp during amino acid starvation, the expression of NudG or MutT supported pppGpp accumulation. This lends support to the idea that a reduction in the SpoT hydrolase activity is sufficient to cause the loss of pppGpp accumulation and therefore the phenomenon is independent of hydrolases that target pppGpp, such as GppA.

Deciphering the regulatory genome of Escherichia coli, one hundred promoters at a time

Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most well-studied of organisms, the bacterium Escherichia coli, for ≈65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, Reg-Seq, that links massively parallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than a E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

Deep sequencing reveals new roles for MuB in transposition immunity and target-capture, and redefines the insular Ter region of E. coli

Background

The target capture protein MuB is responsible for the high efficiency of phage Mu transposition within the E. coli genome. However, some targets are off-limits, such as regions immediately outside the Mu ends (cis-immunity) as well as the entire ~ 37 kb genome of Mu (Mu genome immunity). Paradoxically, MuB is responsible for cis-immunity and is also implicated in Mu genome immunity, but via different mechanisms. This study was undertaken to dissect the role of MuB in target choice in vivo.

Results

We tracked Mu transposition from six different starting locations on the E. coli genome, in the presence and absence of MuB. The data reveal that Mu’s ability to sample the entire genome during a single hop in a clonal population is independent of MuB, and that MuB is responsible for cis-immunity, plays a minor role in Mu genome immunity, and facilitates insertions into transcriptionally active regions. Unexpectedly, transposition patterns in the absence of MuB have helped extend the boundaries of the insular Ter segment of the E. coli genome.

Conclusions

The results in this study demonstrate unambiguously the operation of two distinct mechanisms of Mu target immunity, only one of which is wholly dependent on MuB. The study also reveals several interesting and hitherto unknown aspects of Mu target choice in vivo, particularly the role of MuB in facilitating the capture of promoter and translation start site targets, likely by displacing macromolecular complexes engaged in gene expression. So also, MuB facilitates transposition into the restricted Ter region of the genome.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

Effective CRISPRa-mediated control of gene expression in bacteria must overcome strict target site requirements

In bacterial systems, CRISPR-Cas transcriptional activation (CRISPRa) has the potential to dramatically expand our ability to regulate gene expression, but we lack predictive rules for designing effective gRNA target sites. Here, we identify multiple features of bacterial promoters that impose stringent requirements on CRISPRa target sites. Notably, we observe narrow, 2-4 base windows of effective sites with a periodicity corresponding to one helical turn of DNA, spanning ~40 bases and centered ~80 bases upstream of the TSS. However, we also identify two features suggesting the potential for broad scope: CRISPRa is effective at a broad range of σ70-family promoters, and an expanded PAM dCas9 allows the activation of promoters that cannot be activated by S. pyogenes dCas9. These results provide a roadmap for future engineering efforts to further expand and generalize the scope of bacterial CRISPRa.

Transcription closed and open complex formation coordinate expression of genes with a shared promoter region

Many genes are spaced closely, allowing coordination without explicit control through shared regulatory elements and molecular interactions. We study the dynamics of a stochastic model of a gene-pair in a head-to-head configuration, sharing promoter elements, which accounts for the rate-limiting steps in transcription initiation. We find that only in specific regions of the parameter space of the rate-limiting steps is orderly coexpression exhibited, suggesting that successful cooperation between closely spaced genes requires the coevolution of compatible rate-limiting step configuration. The model predictions are validated using in vivo single-cell, single-RNA measurements of the dynamics of pairs of genes sharing promoter elements. Our results suggest that, in E. coli, the kinetics of the rate-limiting steps in active transcription can play a central role in shaping the dynamics of gene-pairs sharing promoter elements.

Prediction of sgRNA on-target activity in bacteria by deep learning

BACKGROUND

One of the main challenges for the CRISPR-Cas9 system is selecting optimal single-guide RNAs (sgRNAs). Recently, deep learning has enhanced sgRNA prediction in eukaryotes. However, the prokaryotic chromatin structure is different from eukaryotes, so models trained on eukaryotes may not apply to prokaryotes.

RESULTS

We designed and implemented a convolutional neural network to predict sgRNA activity in Escherichia coli. The network was trained and tested on the recently-released sgRNA activity dataset. Our convolutional neural network achieved excellent performance, yielding average Spearman correlation coefficients of 0.5817, 0.7105, and 0.3602, respectively for Cas9, eSpCas9 and Cas9 with a recA coding region deletion. We confirmed that the sgRNA prediction models trained on prokaryotes do not apply to eukaryotes and vice versa. We adopted perturbation-based approaches to analyze distinct biological patterns between prokaryotic and eukaryotic editing. Then, we improved the predictive performance of the prokaryotic Cas9 system by transfer learning. Finally, we determined that potential off-target scores accumulated on a genome-wide scale affect on-target activity, which could slightly improve on-target predictive performance.

CONCLUSIONS

We developed convolutional neural networks to predict sgRNA activity for wild type and mutant Cas9 in prokaryotes. Our results show that the prediction accuracy of our method is improved over state-of-the-art models.