Analysis of the lambdoid prophage element e14 in the E. coli K-12 genome

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 provides insights for 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 regulatory principles is observed. These results reveal that robustness and metabolic load determine the global regulatory principles that govern gene expression processes, and these principles have broad implications for cellular function.

The DNA damage response of Escherichia coli, revisited: Differential gene expression after replication inhibition

In 1967, in this journal, Evelyn Witkin proposed the existence of a coordinated DNA damage response in Escherichia coli, which later came to be called the "SOS response." We revisited this response using the replication inhibitor azidothymidine (AZT) and RNA-Seq analysis and identified several features. We confirm the induction of classic Save our ship (SOS) loci and identify several genes, including many of the pyrimidine pathway, that have not been previously demonstrated to be DNA damage-inducible. Despite a strong dependence on LexA, these genes lack LexA boxes and their regulation by LexA is likely to be indirect via unknown factors. We show that the transcription factor "stringent starvation protein" SspA is as important as LexA in the regulation of AZT-induced genes and that the genes activated by SspA change dramatically after AZT exposure. Our experiments identify additional LexA-independent DNA damage inducible genes, including 22 small RNA genes, some of which appear to activated by SspA. Motility and chemotaxis genes are strongly down-regulated by AZT, possibly as a result of one of more of the small RNAs or other transcription factors such as AppY and GadE, whose expression is elevated by AZT. Genes controlling the iron siderophore, enterobactin, and iron homeostasis are also strongly induced, independent of LexA. We confirm that IraD antiadaptor protein is induced independent of LexA and that a second antiadaptor, IraM is likewise strongly AZT-inducible, independent of LexA, suggesting that RpoS stabilization via these antiadaptor proteins is an integral part of replication stress tolerance.

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.

SOS genes are rapidly induced while translesion synthesis polymerase activity is temporally regulated

The DNA damage inducible SOS response in bacteria serves to increase survival of the species at the cost of mutagenesis. The SOS response first initiates error-free repair followed by error-prone repair. Here, we have employed a multi-omics approach to elucidate the temporal coordination of the SOS response. Escherichia coli was grown in batch cultivation in bioreactors to ensure highly controlled conditions, and a low dose of the antibiotic ciprofloxacin was used to activate the SOS response while avoiding extensive cell death. Our results show that expression of genes involved in error-free and error-prone repair were both induced shortly after DNA damage, thus, challenging the established perception that the expression of error-prone repair genes is delayed. By combining transcriptomics and a sub-proteomics approach termed signalomics, we found that the temporal segregation of error-free and error-prone repair is primarily regulated after transcription, supporting the current literature. Furthermore, the heterology index (i.e., the binding affinity of LexA to the SOS box) was correlated to the maximum increase in gene expression and not to the time of induction of SOS genes. Finally, quantification of metabolites revealed increasing pyrimidine pools as a late feature of the SOS response. Our results elucidate how the SOS response is coordinated, showing a rapid transcriptional response and temporal regulation of mutagenesis on the protein and metabolite levels.

Word-based GWAS harnesses the rich potential of genomic data for E. coli quinolone resistance

Quinolone resistance presents a growing global health threat. We employed word-based GWAS to explore genomic data, aiming to enhance our understanding of this phenomenon. Unlike traditional variant-based GWAS analyses, this approach simultaneously captures multiple genomic factors, including single and interacting resistance mutations and genes. Analyzing a dataset of 92 genomic E. coli samples from a wastewater treatment plant in Dresden, we identified 54 DNA unitigs significantly associated with quinolone resistance. Remarkably, our analysis not only validated known mutations in gyrA and parC genes and the results of our variant-based GWAS but also revealed new (mutated) genes such as mdfA, the AcrEF-TolC multidrug efflux system, ptrB, and hisI, implicated in antibiotic resistance. Furthermore, our study identified joint mutations in 14 genes including the known gyrA gene, providing insights into potential synergistic effects contributing to quinolone resistance. These findings showcase the exceptional capabilities of word-based GWAS in unraveling the intricate genomic foundations of quinolone resistance.

A multi-scale expression and regulation knowledge base for Escherichia coli

Transcriptomic data is accumulating rapidly; thus, scalable methods for extracting knowledge from this data are critical. Here, we assembled a top-down expression and regulation knowledge base for Escherichia coli. The expression component is a 1035-sample, high-quality RNA-seq compendium consisting of data generated in our lab using a single experimental protocol. The compendium contains diverse growth conditions, including: 9 media; 39 supplements, including antibiotics; 42 heterologous proteins; and 76 gene knockouts. Using this resource, we elucidated global expression patterns. We used machine learning to extract 201 modules that account for 86% of known regulatory interactions, creating the regulatory component. With these modules, we identified two novel regulons and quantified systems-level regulatory responses. We also integrated 1675 curated, publicly-available transcriptomes into the resource. We demonstrated workflows for analyzing new data against this knowledge base via deconstruction of regulation during aerobic transition. This resource illuminates the E. coli transcriptome at scale and provides a blueprint for top-down transcriptomic analysis of non-model organisms.

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.

Molecular Signatures of the Eagle Effect Induced by the Artificial Siderophore Conjugate LP-600 in E. coli

Achieving cellular uptake is a central challenge for novel antibiotics targeting Gram-negative bacterial pathogens. One strategy is to hijack the bacterial iron transport system by siderophore-antibiotic conjugates that are actively imported into the cell. This was realized with the MECAM-ampicillin conjugate LP-600 we recently reported that was highly active against E. coli. In the present study, we investigate a paradoxical regrowth of E. coli upon treatment of LP-600 at concentrations 16-32 times above the minimum inhibitory concentration (MIC). The phenomenon, coined "Eagle-effect" in other systems, was not due to resistance formation, and it occurred for the siderophore conjugate but not for free ampicillin. To investigate the molecular imprint of the Eagle effect, a combined transcriptome and untargeted metabolome analysis was conducted. LP-600 induced the expression of genes involved in iron acquisition, SOS response, and the e14 prophage upon regrowth conditions. The Eagle effect was diminished in the presence of sulbactam, which we ascribe to a putative synergistic antibiotic action but not to β-lactamase inhibition. The study highlights the relevance of the Eagle effect for siderophore conjugates. Through the first systematic -omics investigations, it also demonstrates that the Eagle effect manifests not only in a paradoxical growth but also in unique gene expression and metabolite profiles.

A newly identified prophage-encoded gene, ymfM, causes SOS-inducible filamentation in Escherichia coli

Rod-shaped bacteria such as Escherichia coli can regulate cell division in response to stress, leading to filamentation, a process where cell growth and DNA replication continues in the absence of division, resulting in elongated cells. The classic example of stress is DNA damage which results in the activation of the SOS response. While the inhibition of cell division during SOS has traditionally been attributed to SulA in E. coli, a previous report suggests that the e14 prophage may also encode an SOS-inducible cell division inhibitor, previously named SfiC. However, the exact gene responsible for this division inhibition has remained unknown for over 35 years. A recent high-throughput over-expression screen in E. coli identified the e14 prophage gene, ymfM, as a potential cell division inhibitor. In this study, we show that the inducible expression of ymfM from a plasmid causes filamentation. We show that this expression of ymfM results in the inhibition of Z ring formation and is independent of the well characterised inhibitors of FtsZ ring assembly in E. coli, SulA, SlmA and MinC. We confirm that ymfM is the gene responsible for the SfiC phenotype as it contributes to the filamentation observed during the SOS response. This function is independent of SulA, highlighting that multiple alternative division inhibition pathways exist during the SOS response. Our data also highlight that our current understanding of cell division regulation during the SOS response is incomplete and raises many questions regarding how many inhibitors there actually are and their purpose for the survival of the organism.Importance:Filamentation is an important biological mechanism which aids in the survival, pathogenesis and antibiotic resistance of bacteria within different environments, including pathogenic bacteria such as uropathogenic Escherichia coli Here we have identified a bacteriophage-encoded cell division inhibitor which contributes to the filamentation that occurs during the SOS response. Our work highlights that there are multiple pathways that inhibit cell division during stress. Identifying and characterising these pathways is a critical step in understanding survival tactics of bacteria which become important when combating the development of bacterial resistance to antibiotics and their pathogenicity.

Mechanism of Resistance Development in E. coli against TCAT, a Trimethoprim-Based Photoswitchable Antibiotic

During the last decades, a continuous rise of multi-drug resistant pathogens has threatened antibiotic efficacy. To tackle this key challenge, novel antimicrobial therapies are needed with increased specificity for the site of infection. Photopharmacology could enable such specificity by allowing for the control of antibiotic activity with light, as exemplified by trans/cis-tetra-ortho-chloroazobenzene-trimethoprim (TCAT) conjugates. Resistance development against the on (irradiated, TCATa) and off (thermally adapted, TCATd) states of TCAT were compared to that of trimethoprim (TMP) in Escherichia coli mutant strain CS1562. Genomics and transcriptomics were used to explore the acquired resistance. Although TCAT shows TMP-like dihydrofolate reductase (DHFR) inhibition in vitro, transcriptome analyses show different responses in acquired resistance. Resistance against TCATa (on) relies on the production of exopolysaccharides and overexpression of TolC. While resistance against TCATd (off) follows a slightly different gene expression profile, both indicate hampering the entrance of the molecule into the cell. Conversely, resistance against TMP is based on alterations in cell metabolism towards a more persister-like phenotype, as well as alteration of expression levels of enzymes involved in the folate biosynthesis. This study provides a deeper understanding of the development of new therapeutic strategies and the consequences on resistance development against photopharmacological drugs.

An improved direct metamobilome approach increases the detection of larger-sized circular elements across kingdoms

Mobile genetic elements (MGEs) are instrumental in natural prokaryotic genome editing, permitting genome plasticity and allowing microbes to accumulate genetic diversity. MGEs serve as a vast communal gene pool and include DNA elements such as plasmids and bacteriophages (phages) among others. These mobile DNA elements represent a human health risk as they can introduce new traits, such as antibiotic resistance or virulence, to a bacterial strain. Sequencing libraries targeting environmental circular MGEs, referred to as metamobilomes, may broaden our current understanding of the mechanisms behind the mobility, prevalence and content of these elements. However, metamobilomics is affected by a severe bias towards small circular elements, introduced by multiple displacement amplification (MDA). MDA is typically used to overcome limiting DNA quantities after the removal of non-circular DNA during library preparations. By examining the relationship between sequencing coverage and the size of circular MGEs in paired metamobilome datasets with and without MDA, we show that larger circular elements are lost when using MDA. This study is the first to systematically demonstrate that MDA is detrimental to detecting larger-sized plasmids if small plasmids are present. It is also the first to show that MDA can be omitted when using enzyme-based DNA fragmentation and PCR in library preparation kits such as Nextera XT® from Illumina.

Bacteriophages of Shiga Toxin-Producing Escherichia coli and Their Contribution to Pathogenicity

Shiga toxins (Stx) of Shiga toxin-producing Escherichia coli (STEC) are generally encoded in the genome of lambdoid bacteriophages, which spend the most time of their life cycle integrated as prophages in specific sites of the bacterial chromosome. Upon spontaneous induction or induction by chemical or physical stimuli, the stx genes are co-transcribed together with the late phase genes of the prophages. After being assembled in the cytoplasm, and after host cell lysis, mature bacteriophage particles are released into the environment, together with Stx. As members of the group of lambdoid phages, Stx phages share many genetic features with the archetypical temperate phage Lambda, but are heterogeneous in their DNA sequences due to frequent recombination events. In addition to Stx phages, the genome of pathogenic STEC bacteria may contain numerous prophages, which are either cryptic or functional. These prophages may carry foreign genes, some of them related to virulence, besides those necessary for the phage life cycle. Since the production of one or more Stx is considered the major pathogenicity factor of STEC, we aim to highlight the new insights on the contribution of Stx phages and other STEC phages to pathogenicity.

Complete Annotated Genome Sequence of the Salmonella enterica Serovar Typhimurium LT7 Strain STK003, Historically Used in Gene Transfer Studies

The genome of Salmonella enterica serovar Typhimurium LT7 comprises a chromosome and two plasmids. One plasmid is very close to pSLT of Salmonella Typhimurium LT2; the second harbors a shufflon region. Prophage content is distinct: LT7 lacks Fels-1, while Gifsy-1 and Fels-2 show island-like divergence and likely programmed inversion, respectively.

The genome of Salmonella enterica serovar Typhimurium LT7 comprises a chromosome and two plasmids. One plasmid is very close to pSLT of Salmonella Typhimurium LT2; the second harbors a shufflon region. Prophage content is distinct: LT7 lacks Fels-1, while Gifsy-1 and Fels-2 show island-like divergence and likely programmed inversion, respectively.

Engineering Escherichia coli for methanol-dependent growth on glucose for metabolite production

Synthetic methylotrophy aims to engineer methane and methanol utilization pathways in platform hosts like Escherichia coli for industrial bioprocessing of natural gas and biogas. While recent attempts to engineer synthetic methanol auxotrophs have proved successful, these studies focused on scarce and expensive co-substrates. Here, we engineered E. coli for methanol-dependent growth on glucose, an abundant and inexpensive co-substrate, via deletion of glucose 6-phosphate isomerase (pgi), phosphogluconate dehydratase (edd), and ribose 5-phosphate isomerases (rpiAB). Since the parental strain did not exhibit methanol-dependent growth on glucose in minimal medium, we first achieved methanol-dependent growth via amino acid supplementation and used this medium to evolve the strain for methanol-dependent growth in glucose minimal medium. The evolved strain exhibited a maximum growth rate of 0.15 h^-1 in glucose minimal medium with methanol, which is comparable to that of other synthetic methanol auxotrophs. Whole genome sequencing and ¹³C-metabolic flux analysis revealed the causative mutations in the evolved strain. A mutation in the phosphotransferase system enzyme I gene (ptsI) resulted in a reduced glucose uptake rate to maintain a one-to-one molar ratio of substrate utilization. Deletion of the e14 prophage DNA region resulted in two non-synonymous mutations in the isocitrate dehydrogenase (icd) gene, which reduced TCA cycle carbon flux to maintain the internal redox state. In high cell density glucose fed-batch fermentation, methanol-dependent acetone production resulted in 22% average carbon labeling of acetone from ¹³C-methanol, which far surpasses that of the previous best (2.4%) found with methylotrophic E. coli Δpgi. This study addresses the need to identify appropriate co-substrates for engineering synthetic methanol auxotrophs and provides a basis for the next steps toward industrial one-carbon bioprocessing.

Metabolic remodeling in Escherichia coli MG1655. A prophage e14-encoded small RNA, co293, post-transcriptionally regulates transcription factors HcaR and FadR

In previous studies, we have shown the existence of metabolic remodeling in glucose-grown Escherichia coli MG1655 cells expressing the esterase Orf306 from the opd island of Sphingobium fuliginis. We now show that Orf306-dependent metabolic remodeling is due to regulation of a novel small RNA (sRNA). Endogenous propionate, produced due to the esterase/lipase activity of Orf306, repressed expression of a novel E. coli sRNA, co293. This sRNA post-transcriptionally regulates expression of the transcription factors HcaR and FadR either by inhibiting translation or by destabilizing their transcripts. Hence, repression of co293 expression elevates the levels of HcaR and FadR with consequent activation of alternative carbon catabolic pathways. HcaR activates the hca and MHP operons leading to upregulation of the phenyl propionate and hydroxy phenyl propionate (HPP) degradation pathways. Similarly, FadR stimulates the expression of the transcription factor IclR which negatively regulates the glyoxylate bypass pathway genes, aceBAK.

Experimental Evolution of Extreme Resistance to Ionizing Radiation in Escherichia coli after 50 Cycles of Selection

In previous work (D. R. Harris et al., J Bacteriol 191:5240-5252, 2009, https://doi.org/10.1128/JB.00502-09; B. T. Byrne et al., Elife 3:e01322, 2014, https://doi.org/10.7554/eLife.01322), we demonstrated that Escherichia coli could acquire substantial levels of resistance to ionizing radiation (IR) via directed evolution. Major phenotypic contributions involved adaptation of organic systems for DNA repair. We have now undertaken an extended effort to generate E. coli populations that are as resistant to IR as Deinococcus radiodurans After an initial 50 cycles of selection using high-energy electron beam IR, four replicate populations exhibit major increases in IR resistance but have not yet reached IR resistance equivalent to D. radiodurans Regular deep sequencing reveals complex evolutionary patterns with abundant clonal interference. Prominent IR resistance mechanisms involve novel adaptations to DNA repair systems and alterations in RNA polymerase. Adaptation is highly specialized to resist IR exposure, since isolates from the evolved populations exhibit highly variable patterns of resistance to other forms of DNA damage. Sequenced isolates from the populations possess between 184 and 280 mutations. IR resistance in one isolate, IR9-50-1, is derived largely from four novel mutations affecting DNA and RNA metabolism: RecD A90E, RecN K429Q, and RpoB S72N/RpoC K1172I. Additional mechanisms of IR resistance are evident.IMPORTANCE Some bacterial species exhibit astonishing resistance to ionizing radiation, with Deinococcus radiodurans being the archetype. As natural IR sources rarely exceed mGy levels, the capacity of Deinococcus to survive 5,000 Gy has been attributed to desiccation resistance. To understand the molecular basis of true extreme IR resistance, we are using experimental evolution to generate strains of Escherichia coli with IR resistance levels comparable to Deinococcus Experimental evolution has previously generated moderate radioresistance for multiple bacterial species. However, these efforts could not take advantage of modern genomic sequencing technologies. In this report, we examine four replicate bacterial populations after 50 selection cycles. Genomic sequencing allows us to follow the genesis of mutations in populations throughout selection. Novel mutations affecting genes encoding DNA repair proteins and RNA polymerase enhance radioresistance. However, more contributors are apparent.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in 'grounding' of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as 'genetic buffer zones' thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

Genome rearrangements in Escherichia coli during de novo acquisition of resistance to a single antibiotic or two antibiotics successively

Background

The ability of bacteria to acquire resistance to antibiotics relies to a large extent on their capacity for genome modification. Prokaryotic genomes are highly plastic and can utilize horizontal gene transfer, point mutations, and gene deletions or amplifications to realize genome expansion and rearrangements. The contribution of point mutations to de novo acquisition of antibiotic resistance is well-established. In this study, the internal genome rearrangement of Escherichia coli during to de novo acquisition of antibiotic resistance was investigated using whole-genome sequencing.

Results

Cells were made resistant to one of the four antibiotics and subsequently to one of the three remaining. This way the initial genetic rearrangements could be documented together with the effects of an altered genetic background on subsequent development of resistance. A DNA fragment including ampC was amplified by a factor sometimes exceeding 100 as a result of exposure to amoxicillin. Excision of prophage e14 was observed in many samples with a double exposure history, but not in cells exposed to a single antibiotic, indicating that the activation of the SOS stress response alone, normally the trigger for excision, was not sufficient to cause excision of prophage e14. Partial deletion of clpS and clpA occurred in strains exposed to enrofloxacin and tetracycline. Other deletions were observed in some strains, but not in replicates with the exact same exposure history. Various insertion sequence transpositions correlated with exposure to specific antibiotics.

Conclusions

Many of the genome rearrangements have not been reported before to occur during resistance development. The observed correlation between genome rearrangements and specific antibiotic pressure, as well as their presence in independent replicates indicates that these events do not occur randomly. Taken together, the observed genome rearrangements illustrate the plasticity of the E. coli genome when exposed to antibiotic stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5353-y) contains supplementary material, which is available to authorized users.

Genomewide phenotypic analysis of growth, cell morphogenesis, and cell cycle events in Escherichia coli

Cell size, cell growth, and cell cycle events are necessarily intertwined to achieve robust bacterial replication. Yet, a comprehensive and integrated view of these fundamental processes is lacking. Here, we describe an image‐based quantitative screen of the single‐gene knockout collection of Escherichia coli and identify many new genes involved in cell morphogenesis, population growth, nucleoid (bulk chromosome) dynamics, and cell division. Functional analyses, together with high‐dimensional classification, unveil new associations of morphological and cell cycle phenotypes with specific functions and pathways. Additionally, correlation analysis across ~4,000 genetic perturbations shows that growth rate is surprisingly not predictive of cell size. Growth rate was also uncorrelated with the relative timings of nucleoid separation and cell constriction. Rather, our analysis identifies scaling relationships between cell size and nucleoid size and between nucleoid size and the relative timings of nucleoid separation and cell division. These connections suggest that the nucleoid links cell morphogenesis to the cell cycle.

Transcriptional responses of Escherichia coli during recovery from inorganic or organic mercury exposure

Background

The protean chemical properties of mercury have long made it attractive for diverse applications, but its toxicity requires great care in its use, disposal, and recycling. Mercury occurs in multiple chemical forms, and the molecular basis for the distinct toxicity of its various forms is only partly understood. Global transcriptomics applied over time can reveal how a cell recognizes a toxicant and what cellular subsystems it marshals to repair and recover from the damage. The longitudinal effects on the transcriptome of exponential phase E. coli were compared during sub-acute exposure to mercuric chloride (HgCl₂) or to phenylmercuric acetate (PMA) using RNA-Seq.

Results

Differential gene expression revealed common and distinct responses to the mercurials throughout recovery. Cultures exhibited growth stasis immediately after each mercurial exposure but returned to normal growth more quickly after PMA exposure than after HgCl₂ exposure. Correspondingly, PMA rapidly elicited up-regulation of a large number of genes which continued for 30 min, whereas fewer genes were up-regulated early after HgCl₂ exposure only some of which overlapped with PMA up-regulated genes. By 60 min gene expression in PMA-exposed cells was almost indistinguishable from unexposed cells, but HgCl₂ exposed cells still had many differentially expressed genes. Relative expression of energy production and most metabolite uptake pathways declined with both compounds, but nearly all stress response systems were up-regulated by one or the other mercurial during recovery.

Conclusions

Sub-acute exposure influenced expression of ~45% of all genes with many distinct responses for each compound, reflecting differential biochemical damage by each mercurial and the corresponding resources available for repair. This study is the first global, high-resolution view of the transcriptional responses to any common toxicant in a prokaryotic model system from exposure to recovery of active growth. The responses provoked by these two mercurials in this model bacterium also provide insights about how higher organisms may respond to these ubiquitous metal toxicants.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4413-z) contains supplementary material, which is available to authorized users.

Interaction of Type IV Toxin/Antitoxin Systems in Cryptic Prophages of Escherichia coli K-12

Toxin/antitoxin (TA) systems are widespread in prokaryotic chromosomes and in mobile genetic elements including plasmids and prophages. The first characterized Type IV TA system CbtA/CbeA was found in cryptic prophage CP4-44 in Escherichia coli K-12. Two homologous TA loci of CbtA/CbeA also reside in cryptic prophages of E. coli K-12, YkfI/YafW in CP4-6 and YpjF/YfjZ in CP4-57. In this study, we demonstrated that YkfI and YpjF inhibited cell growth and led to the formation of "lemon-shaped" cells. Prolonged overproduction of YkfI led to the formation of "gourd-shaped" cells and immediate cell lysis. YafW and YfjZ can neutralize the toxicity of YkfI or YpjF. Furthermore, we found that YkfI and YpjF interacted with cell division protein FtsZ in E. coli, but ectopic expression in Pseudomonas and Shewanella did not cause the formation of "lemon-shaped" cells. Moreover, deletion of all of the three toxin genes together decreased resistance to oxidative stress and deletion of the antitoxin genes increased early biofilm formation. Collectively, these results demonstrated that the homologous Type IV TA systems in E. coli may target cell division protein FtsZ in E. coli and may have different physiological functions in E. coli.

Genome Sequencing Analysis of Atypical Shigella flexneri Isolated in Korea

OBJECTIVES

An atypical Shigella flexneri strain with a plural agglutination pattern [i.e., reacting not only with serum samples containing type antigen II but also with serum samples containing group antigens (3)4 and 7(8)] was selected for genome sequencing, with the aim of obtaining additional comparative information about such strains.

METHODS

The genomic DNA of atypical S. flexneri strain NCCP 15744 was sequenced using an Ion Torrent PGM sequencing machine (Life Technologies, USA). The raw sequence data were preprocessed and reference-assembled in the CLC Assembly Cell software (version 4.0.6; CLC bio, USA).

RESULTS

Ion Torrent sequencing produced 1,450,025 single reads with an average length of 144 bp, totaling ~209 Mbp. The NCCP 15744 genome is composed of one chromosome and four plasmids and contains a gtrX gene. Among the published genome sequences of S. flexneri strains, including 2457T, Sf301, and 2002017, strain NCCP 15744 showed high similarity with strain 2002017. The differences between NCCP 15744 and 2002017 are as follows: i) NCCP 15744 carries four plasmids whereas 2002017 carries five; ii) 19 genes (including CI, CII, and cro) were lost in the SHI-O genomic island of NCCP 15744 and six genes were gained as compared with strain 2002017.

CONCLUSION

Strain NCCP 15744 is genetically similar to 2002017, but these two strains have different multilocus sequence types and serotypes. The exact reason is unclear, but the 19 lost genes may be responsible for the atypical seroconversion of strain NCCP 15744.

Systematic identification and quantification of phase variation in commensal and pathogenic Escherichia coli

Bacteria have been shown to generate constant genetic variation in a process termed phase variation. We present a tool based on whole genome sequencing that allows detection and quantification of coexisting genotypes mediated by genomic inversions in bacterial cultures. We tested our method on widely used strains of Escherichia coli, and detected stable and reproducible phase variation in several invertible loci. These are shown here to be responsible for maintaining constant variation in populations grown from a single colony. Applying this tool on other bacterial strains can shed light on how pathogens adjust to hostile environments by diversifying their genomes.

A replication-inhibited unsegregated nucleoid at mid-cell blocks Z-ring formation and cell division independently of SOS and the SlmA nucleoid occlusion protein in Escherichia coli

Chromosome replication and cell division of Escherichia coli are coordinated with growth such that wild-type cells divide once and only once after each replication cycle. To investigate the nature of this coordination, the effects of inhibiting replication on Z-ring formation and cell division were tested in both synchronized and exponentially growing cells with only one replicating chromosome. When replication elongation was blocked by hydroxyurea or nalidixic acid, arrested cells contained one partially replicated, compact nucleoid located mid-cell. Cell division was strongly inhibited at or before the level of Z-ring formation. DNA cross-linking by mitomycin C delayed segregation, and the accumulation of about two chromosome equivalents at mid-cell also blocked Z-ring formation and cell division. Z-ring inhibition occurred independently of SOS, SlmA-mediated nucleoid occlusion, and MinCDE proteins and did not result from a decreased FtsZ protein concentration. We propose that the presence of a compact, incompletely replicated nucleoid or unsegregated chromosome masses at the normal mid-cell division site inhibits Z-ring formation and that the SOS system, SlmA, and MinC are not required for this inhibition.

Harnessing single cell sorting to identify cell division genes and regulators in bacteria

Cell division is an essential cellular process that requires an array of known and unknown proteins for its spatial and temporal regulation. Here we develop a novel, high-throughput screening method for the identification of bacterial cell division genes and regulators. The method combines the over-expression of a shotgun genomic expression library to perturb the cell division process with high-throughput flow cytometry sorting to screen many thousands of clones. Using this approach, we recovered clones with a filamentous morphology for the model bacterium, Escherichia coli. Genetic analysis revealed that our screen identified both known cell division genes, and genes that have not previously been identified to be involved in cell division. This novel screening strategy is applicable to a wide range of organisms, including pathogenic bacteria, where cell division genes and regulators are attractive drug targets for antibiotic development.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the 'Big Bang' theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

Viral ancestors of antiviral systems

All life must survive their corresponding viruses. Thus antiviral systems are essential in all living organisms. Remnants of virus derived information are also found in all life forms but have historically been considered mostly as junk DNA. However, such virus derived information can strongly affect host susceptibility to viruses. In this review, I evaluate the role viruses have had in the origin and evolution of host antiviral systems. From Archaea through bacteria and from simple to complex eukaryotes I trace the viral components that became essential elements of antiviral immunity. I conclude with a reexamination of the ‘Big Bang’ theory for the emergence of the adaptive immune system in vertebrates by horizontal transfer and note how viruses could have and did provide crucial and coordinated features.

Cryptic prophages help bacteria cope with adverse environments

Phages are the most abundant entity in the biosphere and outnumber bacteria by a factor of 10. Phage DNA may also constitute 20% of bacterial genomes; however, its role is ill defined. Here, we explore the impact of cryptic prophages on cell physiology by precisely deleting all nine prophage elements (166 kbp) using Escherichia coli. We find that cryptic prophages contribute significantly to resistance to sub-lethal concentrations of quinolone and β-lactam antibiotics primarily through proteins that inhibit cell division (for example, KilR of rac and DicB of Qin). Moreover, the prophages are beneficial for withstanding osmotic, oxidative and acid stresses, for increasing growth, and for influencing biofilm formation. Prophage CPS-53 proteins YfdK, YfdO and YfdS enhanced resistance to oxidative stress, prophages e14, CPS-53 and CP4-57 increased resistance to acid, and e14 and rac proteins increased early biofilm formation. Therefore, cryptic prophages provide multiple benefits to the host for surviving adverse environmental conditions.

Up to 20% of bacterial genomes are made up of cryptic prophages, but their function is relatively unknown. In this study, the authors demonstrate that prophages influence the response of the host cell to stress and provide a competitive growth advantage in the presence of antibiotics.

Single-gene deletion mutants of Escherichia coli with altered sensitivity to bicyclomycin, an inhibitor of transcription termination factor Rho

We have screened the entire KEIO collection of 3,985 single-gene knockouts in Escherichia coli for increased susceptibility or resistance to the antibiotic bicyclomycin (BCM), a potent inhibitor of the transcription termination factor Rho. We also compared the results to those of a recent study we conducted with a large set of antibiotics (A. Liu et al., Antimicrob. Agents Chemother. 54:1393-1403, 2010). We find that deletions of many different types of genes increase sensitivity to BCM. Some of these are involved in multidrug sensitivity/resistance, whereas others are specific for BCM. Mutations in a number of DNA recombination and repair genes increase BCM sensitivity, indicating that DNA damage leading to single- and double-strand breaks is a downstream effect of Rho inhibition. MDS42, which is deleted for all cryptic prophages and insertion elements (G. Posfai et al., Science 312:1044-1046, 2006), or W3102 deleted for the rac prophage-encoded kil gene, are partially resistant to BCM (C. J. Cardinale et al., Science 230:935-938, 2008). Deletion of cryptic prophages also overcomes the increased BCM sensitivity in some but not all mutants examined here. Deletion of the hns gene renders the cell more sensitive to BCM even in the Δkil or MDS42 background. This suggests that BCM activates additional modes of cell death independent of Kil and that these could provide a target to potentiate BCM killing.

Genetic basis of growth adaptation of Escherichia coli after deletion of pgi, a major metabolic gene

Bacterial survival requires adaptation to different environmental perturbations such as exposure to antibiotics, changes in temperature or oxygen levels, DNA damage, and alternative nutrient sources. During adaptation, bacteria often develop beneficial mutations that confer increased fitness in the new environment. Adaptation to the loss of a major non-essential gene product that cripples growth, however, has not been studied at the whole-genome level. We investigated the ability of Escherichia coli K-12 MG1655 to overcome the loss of phosphoglucose isomerase (pgi) by adaptively evolving ten replicates of E. coli lacking pgi for 50 days in glucose M9 minimal medium and by characterizing endpoint clones through whole-genome re-sequencing and phenotype profiling. We found that 1) the growth rates for all ten endpoint clones increased approximately 3-fold over the 50-day period; 2) two to five mutations arose during adaptation, most frequently in the NADH/NADPH transhydrogenases udhA and pntAB and in the stress-associated sigma factor rpoS; and 3) despite similar growth rates, at least three distinct endpoint phenotypes developed as defined by different rates of acetate and formate secretion. These results demonstrate that E. coli can adapt to the loss of a major metabolic gene product with only a handful of mutations and that adaptation can result in multiple, alternative phenotypes.

Roles of RubisCO and the RubisCO-like protein in 5-methylthioadenosine metabolism in the Nonsulfur purple bacterium Rhodospirillum rubrum

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) catalyzes the assimilation of atmospheric CO(2) into organic matter and is thus central to the existence of life on earth. The beginning of the 2000s was marked by the discovery of a new family of proteins, the RubisCO-like proteins (RLPs), which are structural homologs of RubisCO. RLPs are unable to catalyze CO(2) fixation. The RLPs from Chlorobaculum tepidum, Bacillus subtilis, Geobacillus kaustophilus, and Microcystis aeruginosa have been shown to participate in sulfur metabolism. Whereas the precise function of C. tepidum RLP is unknown, the B. subtilis, G. kaustophilus, and M. aeruginosa RLPs function as tautomerases/enolases in a methionine salvage pathway (MSP). Here, we show that the form II RubisCO enzyme from the nonsulfur purple bacterium Rhodospirillum rubrum is also able to function as an enolase in vivo as part of an MSP, but only under anaerobic conditions. However, unlike B. subtilis RLP, R. rubrum RLP does not catalyze the enolization of 2,3-diketo-5-methylthiopentyl-1-phosphate. Instead, under aerobic growth conditions, R. rubrum RLP employs another intermediate of the MSP, 5-methylthioribulose-1-phosphate, as a substrate, resulting in the formation of different products. To further determine the interrelationship between RubisCOs and RLPs (and the potential integration of cellular carbon and sulfur metabolism), the functional roles of both RubisCO and RLP have been examined in vivo via the use of specific knockout strains and complementation studies of R. rubrum. The presence of functional, yet separate, MSPs in R. rubrum under both aerobic (chemoheterotrophic) and anaerobic (photoheterotrophic) growth conditions has not been observed previously in any organism. Moreover, the aerobic and anaerobic sulfur salvage pathways appear to be differentially controlled, with novel and previously undescribed steps apparent for sulfur salvage in this organism.

Differential binding of Escherichia coli McrA protein to DNA sequences that contain the dinucleotide m5CpG

The Escherichia coli McrA protein, a putative C⁵-methylcytosine/C⁵-hydroxyl methylcytosine-specific nuclease, binds DNA with symmetrically methylated HpaII sequences (Cm5CGG), but its precise recognition sequence remains undefined. To determine McrA’s binding specificity, we cloned and expressed recombinant McrA with a C-terminal StrepII tag (rMcrA-S) to facilitate protein purification and affinity capture of human DNA fragments with m5C residues. Sequence analysis of a subset of these fragments and electrophoretic mobility shift assays with model methylated and unmethylated oligonucleotides suggest that N(Y > R) m5CGR is the canonical binding site for rMcrA-S. In addition to binding HpaII-methylated double-stranded DNA, rMcrA-S binds DNA containing a single, hemimethylated HpaII site; however, it does not bind if A, C, T or U is placed across from the m5C residue, but does if I is opposite the m5C. These results provide the first systematic analysis of McrA’s in vitro binding specificity.

Searching for a "hidden" prophage in a marine bacterium

Prophages are common in many bacterial genomes. Distinguishing putatively viable prophages from nonviable sequences can be a challenge, since some prophages are remnants of once-functional prophages that have been rendered inactive by mutational changes. In some cases, a putative prophage may be missed due to the lack of recognizable prophage loci. The genome of a marine roseobacter, Roseovarius nubinhibens ISM (hereinafter referred to as ISM), was recently sequenced and was reported to contain no intact prophage based on customary bioinformatic analysis. However, prophage induction experiments performed with this organism led to a different conclusion. In the laboratory, virus-like particles in the ISM culture increased more than 3 orders of magnitude following induction with mitomycin C. After careful examination of the ISM genome sequence, a putative prophage (ISM-pro1) was identified. Although this prophage contains only minimal phage-like genes, we demonstrated that this "hidden" prophage is inducible. Genomic analysis and reannotation showed that most of the ISM-pro1 open reading frames (ORFs) display the highest sequence similarity with Rhodobacterales bacterial genes and some ORFs are only distantly related to genes of other known phages or prophages. Comparative genomic analyses indicated that ISM-pro1-like prophages or prophage remnants are also present in other Rhodobacterales genomes. In addition, the lysis of ISM by this previously unrecognized prophage appeared to increase the production of gene transfer agents (GTAs). Our study suggests that a combination of in silico genomic analyses and experimental laboratory work is needed to fully understand the lysogenic features of a given bacterium.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Directed evolution of ionizing radiation resistance in Escherichia coli

We have generated extreme ionizing radiation resistance in a relatively sensitive bacterial species, Escherichia coli, by directed evolution. Four populations of Escherichia coli K-12 were derived independently from strain MG1655, with each specifically adapted to survive exposure to high doses of ionizing radiation. D(37) values for strains isolated from two of the populations approached that exhibited by Deinococcus radiodurans. Complete genomic sequencing was carried out on nine purified strains derived from these populations. Clear mutational patterns were observed that both pointed to key underlying mechanisms and guided further characterization of the strains. In these evolved populations, passive genomic protection is not in evidence. Instead, enhanced recombinational DNA repair makes a prominent but probably not exclusive contribution to genome reconstitution. Multiple genes, multiple alleles of some genes, multiple mechanisms, and multiple evolutionary pathways all play a role in the evolutionary acquisition of extreme radiation resistance. Several mutations in the recA gene and a deletion of the e14 prophage both demonstrably contribute to and partially explain the new phenotype. Mutations in additional components of the bacterial recombinational repair system and the replication restart primosome are also prominent, as are mutations in genes involved in cell division, protein turnover, and glutamate transport. At least some evolutionary pathways to extreme radiation resistance are constrained by the temporally ordered appearance of specific alleles.

The HU regulon is composed of genes responding to anaerobiosis, acid stress, high osmolarity and SOS induction

BACKGROUND

The Escherichia coli heterodimeric HU protein is a small DNA-bending protein associated with the bacterial nucleoid. It can introduce negative supercoils into closed circular DNA in the presence of topoisomerase I. Cells lacking HU grow very poorly and display many phenotypes.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed the transcription profile of every Escherichia coli gene in the absence of one or both HU subunits. This genome-wide in silico transcriptomic approach, performed in parallel with in vivo genetic experimentation, defined the HU regulon. This large regulon, which comprises 8% of the genome, is composed of four biologically relevant gene classes whose regulation responds to anaerobiosis, acid stress, high osmolarity, and SOS induction.

CONCLUSIONS/SIGNIFICANCE

The regulation a large number of genes encoding enzymes involved in energy metabolism and catabolism pathways by HU explains the highly pleiotropic phenotype of HU-deficient cells. The uniform chromosomal distribution of the many operons regulated by HU strongly suggests that the transcriptional and nucleoid architectural functions of HU constitute two aspects of a unique protein-DNA interaction mechanism.

Cell death upon epigenetic genome methylation: a novel function of methyl-specific deoxyribonucleases

BACKGROUND

Alteration in epigenetic methylation can affect gene expression and other processes. In Prokaryota, DNA methyltransferase genes frequently move between genomes and present a potential threat. A methyl-specific deoxyribonuclease, McrBC, of Escherichia coli cuts invading methylated DNAs. Here we examined whether McrBC competes with genome methylation systems through host killing by chromosome cleavage.

RESULTS

McrBC inhibited the establishment of a plasmid carrying a PvuII methyltransferase gene but lacking its recognition sites, likely through the lethal cleavage of chromosomes that became methylated. Indeed, its phage-mediated transfer caused McrBC-dependent chromosome cleavage. Its induction led to cell death accompanied by chromosome methylation, cleavage and degradation. RecA/RecBCD functions affect chromosome processing and, together with the SOS response, reduce lethality. Our evolutionary/genomic analyses of McrBC homologs revealed: a wide distribution in Prokaryota; frequent distant horizontal transfer and linkage with mobility-related genes; and diversification in the DNA binding domain. In these features, McrBCs resemble type II restriction-modification systems, which behave as selfish mobile elements, maintaining their frequency by host killing. McrBCs are frequently found linked with a methyltransferase homolog, which suggests a functional association.

CONCLUSIONS

Our experiments indicate McrBC can respond to genome methylation systems by host killing. Combined with our evolutionary/genomic analyses, they support our hypothesis that McrBCs have evolved as mobile elements competing with specific genome methylation systems through host killing. To our knowledge, this represents the first report of a defense system against epigenetic systems through cell death.

Photorhabdus luminescens genes induced upon insect infection

BACKGROUND

Photorhabdus luminescens is a Gram-negative luminescent enterobacterium and a symbiote to soil nematodes belonging to the species Heterorhabditis bacteriophora. P.luminescens is simultaneously highly pathogenic to insects. This bacterium exhibits a complex life cycle, including one symbiotic stage characterized by colonization of the upper nematode gut, and a pathogenic stage, characterized by release from the nematode into the hemocoel of insect larvae, resulting in rapid insect death caused by bacterial toxins. P. luminescens appears to sense and adapt to the novel host environment upon changing hosts, which facilitates the production of factors involved in survival within the host, host-killing, and -exploitation.

RESULTS

A differential fluorescence induction (DFI) approach was applied to identify genes that are up-regulated in the bacterium after infection of the insect host Galleria mellonella. For this purpose, a P. luminescens promoter-trap library utilizing the mCherry fluorophore as a reporter was constructed, and approximately 13,000 clones were screened for fluorescence induction in the presence of a G. mellonella larvae homogenate. Since P. luminescens has a variety of regulators that potentially sense chemical molecules, like hormones, the screen for up-regulated genes or operons was performed in vitro, excluding physicochemical signals like oxygen, temperature or osmolarity as variables. Clones (18) were obtained exhibiting at least 2.5-fold induced fluorescence and regarded as specific responders to insect homogenate. In combination with a bioinformatics approach, sequence motifs were identified in these DNA-fragments that are similar to 29 different promoters within the P. luminescens genome. By cloning each of the predicted promoters upstream of the reporter gene, induction was verified for 27 promoters in vitro, and for 24 promoters in viable G. mellonella larvae. Among the validated promoters are some known to regulate the expression of toxin genes, including tccC1 (encoding an insecticidal toxin complex), and others encoding putative toxins. A comparably high number of metabolic genes or operons were observed to be induced upon infection; among these were eutABC, hutUH, and agaZSVCD, which encode proteins involved in ethanolamine, histidine and tagatose degradation, respectively. The results reflect rearrangements in metabolism and the use of other metabolites available from the insect. Furthermore, enhanced activity of promoters controlling the expression of genes encoding enzymes linked to antibiotic production and/or resistance was observed. Antibiotic production and resistance may influence competition with other bacteria, and thus might be important for a successful infection. Lastly, several genes of unknown function were identified that may represent novel pathogenicity factors.

CONCLUSION

We show that a DFI screen is useful for identifying genes or operons induced by chemical stimuli, such as diluted insect homogenate. A bioinformatics comparison of motifs similar to known promoters is a powerful tool for identifying regulated genes or operons. We conclude that signals for the regulation of those genes or operons induced in P. luminescens upon insect infection may represent a wide variety of compounds that make up the insect host. Our results provide insight into the complex response to the host that occurs in a bacterial pathogen, particularly reflecting the potential for metabolic shifts and other specific changes associated with virulence.

Identification of a novel prophage-like gene cluster actively expressed in both virulent and avirulent strains of Leptospira interrogans serovar Lai

DNA microarray analysis was used to compare the differential gene expression profiles between Leptospira interrogans serovar Lai type strain 56601 and its corresponding attenuated strain IPAV. A 22-kb genomic island covering a cluster of 34 genes (i.e., genes LA0186 to LA0219) was actively expressed in both strains but concomitantly upregulated in strain 56601 in contrast to that of IPAV. Reverse transcription-PCR assays proved that the gene cluster comprised five transcripts. Gene annotation of this cluster revealed characteristics of a putative prophage-like remnant with at least 8 of 34 sequences encoding prophage-like proteins, of which the LA0195 protein is probably a putative prophage CI-like regulator. The transcription initiation activities of putative promoter-regulatory sequences of transcripts I, II, and III, all proximal to the LA0195 gene, were further analyzed in the Escherichia coli promoter probe vector pKK232-8 by assaying the reporter chloramphenicol acetyltransferase (CAT) activities. The strong promoter activities of both transcripts I and II indicated by the E. coli CAT assay were well correlated with the in vitro sequence-specific binding of the recombinant LA0195 protein to the corresponding promoter probes detected by the electrophoresis mobility shift assay. On the other hand, the promoter activity of transcript III was very low in E. coli and failed to show active binding to the LA0195 protein in vitro. These results suggested that the LA0195 protein is likely involved in the transcription of transcripts I and II. However, the identical complete DNA sequences of this prophage remnant from these two strains strongly suggests that possible regulatory factors or signal transduction systems residing outside of this region within the genome may be responsible for the differential expression profiling in these two strains.

Identification of prophages in bacterial genomes by dinucleotide relative abundance difference

Background

Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful.

Methodology

Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method.

Conclusions

The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences.

Contribution of the SitABCD, MntH, and FeoB metal transporters to the virulence of avian pathogenic Escherichia coli O78 strain chi7122

The roles of SitABCD, MntH, and FeoB metal transporters in the virulence of avian pathogenic Escherichia coli (APEC) O78 strain chi7122 were assessed using isogenic mutants in chicken infection models. In a single-strain infection model, compared to chi7122, the Deltasit strain demonstrated reduced colonization of the lungs, liver, and spleen. Complementation of the Deltasit strain restored virulence. In a coinfection model, compared to the virulent APEC strain, the Deltasit strain demonstrated mean 50-fold, 126-fold, and 25-fold decreases in colonization of the lungs, liver, and spleen, respectively. A DeltamntH Deltasit strain was further attenuated, demonstrating reduced persistence in blood and mean 1,400-fold, 954-fold, and 83-fold reduced colonization in the lungs, liver, and spleen, respectively. In coinfections, the DeltafeoB Deltasit strain demonstrated reduced persistence in blood but increased colonization of the liver. The DeltamntH, DeltafeoB, and DeltafeoB DeltamntH strains were as virulent as the wild type in either of the infection models. Strains were also tested for sensitivity to oxidative stress-generating agents. The DeltamntH Deltasit strain was the most sensitive strain and was significantly more sensitive than the other strains to hydrogen peroxide, plumbagin, and paraquat. sit sequences were highly associated with APEC and human extraintestinal pathogenic E. coli compared to commensal isolates and diarrheagenic E. coli. Comparative genomic analyses also demonstrated that sit sequences are carried on conjugative plasmids or associated with phage elements and were likely acquired by distinct genetic events among pathogenic E. coli and Shigella sp. strains. Overall, the results demonstrate that SitABCD contributes to virulence and, together with MntH, to increased resistance to oxidative stress.