Large-scale identification of essential Salmonella genes by trapping lethal insertions

Essential Genes Discovery in Microorganisms by Transposon-Directed Sequencing (Tn-Seq): Experimental Approaches, Major Goals, and Future Perspectives

Essential genes are crucial for microbial viability, playing key roles in both the primary and secondary metabolism. Since mutations in these genes can threaten organism viability, identifying them is challenging. Conditionally essential genes are required only under specific conditions and are important for functions such as virulence, immunity, stress survival, and antibiotic resistance. Transposon-directed sequencing (Tn-Seq) has emerged as a powerful method for identifying both essential and conditionally essential genes. In this review, we explored Tn-Seq workflows, focusing on eubacterial species and some yeast species. A comparison of 14 eubacteria species revealed 133 conserved essential genes, including those involved in cell division (e.g., ftsA, ftsZ), DNA replication (e.g., dnaA, dnaE), ribosomal function, cell wall synthesis (e.g., murB, murC), and amino acid synthesis (e.g., alaS, argS). Many other essential genes lack clear orthologues across different microorganisms, making them specific to each organism studied. Conditionally essential genes were identified in 18 bacterial species grown under various conditions, but their conservation was low, reflecting dependence on specific environments and microorganisms. Advances in Tn-Seq are expected to reveal more essential genes in the near future, deepening our understanding of microbial biology and enhancing our ability to manipulate microbial growth, as well as both the primary and secondary metabolism.

Dual transcriptome based reconstruction of Salmonella-human integrated metabolic network to screen potential drug targets

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a highly adaptive pathogenic bacteria with a serious public health concern due to its increasing resistance to antibiotics. Therefore, identification of novel drug targets for S. Typhimurium is crucial. Here, we first created a pathogen-host integrated genome-scale metabolic network by combining the metabolic models of human and S. Typhimurium, which we further tailored to the pathogenic state by the integration of dual transcriptome data. The integrated metabolic model enabled simultaneous investigation of metabolic alterations in human cells and S. Typhimurium during infection. Then, we used the tailored pathogen-host integrated genome-scale metabolic network to predict essential genes in the pathogen, which are candidate novel drug targets to inhibit infection. Drug target prioritization procedure was applied to these targets, and pabB was chosen as a putative drug target. It has an essential role in 4-aminobenzoic acid (PABA) synthesis, which is an essential biomolecule for many pathogens. A structure based virtual screening was applied through docking simulations to predict candidate compounds that eliminate S. Typhimurium infection by inhibiting pabB. To our knowledge, this is the first comprehensive study for predicting drug targets and drug like molecules by using pathogen-host integrated genome-scale models, dual RNA-seq data and structure-based virtual screening protocols. This framework will be useful in proposing novel drug targets and drugs for antibiotic-resistant pathogens.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Genetic Determinants in Salmonella enterica Serotype Typhimurium Required for Overcoming In Vitro Stressors in the Mimicking Host Environment

Salmonella enterica serotype Typhimurium, a nontyphoidal Salmonella (NTS), results in a range of enteric diseases, representing a major disease burden worldwide. There is still a significant portion of Salmonella genes whose mechanistic basis to overcome host innate defense mechanisms largely remains unknown. Here, we have applied transposon insertion sequencing (Tn-seq) method to unveil the genetic factors required for the growth or survival of S. Typhimurium under various host stressors simulated in vitro. A highly saturating Tn5 library of S. Typhimurium 14028s was subjected to selection during growth in the presence of short-chain fatty acid (100 mM propionate), osmotic stress (3% NaCl), or oxidative stress (1 mM H2O2) or survival in extreme acidic pH (30 min in pH 3) or starvation (12 days in 1× phosphate-buffered saline [PBS]). We have identified a total of 339 conditionally essential genes (CEGs) required to overcome at least one of these conditions mimicking host insults. Interestingly, all eight genes encoding FoF1-ATP synthase subunit proteins were required for fitness in all five stresses. Intriguingly, a total of 88 genes in Salmonella pathogenicity islands (SPI), including SPI-1, SPI-2, SPI-3, SPI-5, SPI-6, and SPI-11, are also required for fitness under the in vitro conditions. Additionally, by comparative analysis of the genes identified in this study and the genes previously shown to be required for in vivo fitness, we identified novel genes (marBCT, envF, barA, hscA, rfaQ, rfbI, and the genes encoding putative proteins STM14_1138, STM14_3334, STM14_4825, and STM_5184) that have compelling potential for the development of vaccines and antibacterial drugs to curb Salmonella infection. IMPORTANCE Salmonella enterica serotype Typhimurium is a major human bacterial pathogen that enters the food chain through meat animals asymptomatically carrying this pathogen. Despite the rich genome sequence data, a significant portion of Salmonella genes remain to be characterized for their potential contributions to virulence. In this study, we used transposon insertion sequencing (Tn-seq) to elucidate the genetic factors required for growth or survival under various host stressors, including short-chain fatty acids, osmotic stress, oxidative stress, extreme acid, and starvation. Among the total of 339 conditionally essential genes (CEGs) that are required under at least one of these five stress conditions were 221 previously known virulence genes required for in vivo fitness during infection in at least one of four animal species, including mice, chickens, pigs, and cattle. This comprehensive map of virulence phenotype-genotype in S. Typhimurium provides a roadmap for further interrogation of the biological functions encoded by the genome of this important human pathogen to survive in hostile host environments.

The enigmatic biology of rickettsiae: recent advances, open questions and outlook

Rickettsiae are obligate intracellular bacteria that can cause life-threatening illnesses and are among the oldest known vector-borne pathogens. Members of this genus are extraordinarily diverse and exhibit a broad host range. To establish intracellular infection, Rickettsia species undergo complex, multistep life cycles that are encoded by heavily streamlined genomes. As a result of reductive genome evolution, rickettsiae are exquisitely tailored to their host cell environment but cannot survive extracellularly. This host-cell dependence makes for a compelling system to uncover novel host-pathogen biology, but it has also hindered experimental progress. Consequently, the molecular details of rickettsial biology and pathogenesis remain poorly understood. With recent advances in molecular biology and genetics, the field is poised to start unraveling the molecular mechanisms of these host-pathogen interactions. Here, we review recent discoveries that have shed light on key aspects of rickettsial biology. These studies have revealed that rickettsiae subvert host cells using mechanisms that are distinct from other better-studied pathogens, underscoring the great potential of the Rickettsia genus for revealing novel biology. We also highlight several open questions as promising areas for future study and discuss the path toward solving the fundamental mysteries of this neglected and emerging human pathogen.

DEG 15, an update of the Database of Essential Genes that includes built-in analysis tools

Essential genes refer to genes that are required by an organism to survive under specific conditions. Studies of the minimal-gene-set for bacteria have elucidated fundamental cellular processes that sustain life. The past five years have seen a significant progress in identifying human essential genes, primarily due to the successful use of CRISPR/Cas9 in various types of human cells. DEG 15, a new release of the Database of Essential Genes (www.essentialgene.org), has provided major advancements, compared to DEG 10. Specifically, the number of eukaryotic essential genes has increased by more than fourfold, and that of prokaryotic ones has more than doubled. Of note, the human essential-gene number has increased by more than tenfold. Moreover, we have developed built-in analysis modules by which users can perform various analyses, such as essential-gene distributions between bacterial leading and lagging strands, sub-cellular localization distribution, enrichment analysis of gene ontology and KEGG pathways, and generation of Venn diagrams to compare and contrast gene sets between experiments. Additionally, the database offers customizable BLAST tools for performing species- and experiment-specific BLAST searches. Therefore, DEG comprehensively harbors updated human-curated essential-gene records among prokaryotes and eukaryotes with built-in tools to enhance essential-gene analysis.

Environmental flexibility does not explain metabolic robustness

Cells show remarkable resilience against genetic and environmental perturbations. However, its evolutionary origin remains obscure. In order to leverage methods of systems biology for examining cellular robustness, a computationally accessible way of quantification is needed. Here, we present an unbiased metric of structural robustness in genome-scale metabolic models based on concepts prevalent in reliability engineering and fault analysis. The probability of failure (PoF) is defined as the (weighted) portion of all possible combinations of loss-of-function mutations that disable network functionality. It can be exactly determined if all essential reactions, synthetic lethal pairs of reactions, synthetic lethal triplets of reactions etc. are known. In theory, these minimal cut sets (MCSs) can be calculated for any network, but for large models the problem remains computationally intractable. Herein, we demonstrate that even at the genome scale only the lowest-cardinality MCSs are required to efficiently approximate the PoF with reasonable accuracy. Building on an improved theoretical understanding, we analysed the robustness of 489 E. coli, Shigella, Salmonella, and fungal genome-scale metabolic models (GSMMs). In contrast to the popular "congruence theory", which explains the origin of genetic robustness as a byproduct of selection for environmental flexibility, we found no correlation between network robustness and the diversity of growth-supporting environments. On the contrary, our analysis indicates that amino acid synthesis rather than carbon metabolism dominates metabolic robustness.

Symmetric activity of DNA polymerases at and recruitment of exonuclease ExoR and of PolA to the Bacillus subtilis replication forks

DNA replication forks are intrinsically asymmetric and may arrest during the cell cycle upon encountering modifications in the DNA. We have studied real time dynamics of three DNA polymerases and an exonuclease at a single molecule level in the bacterium Bacillus subtilis. PolC and DnaE work in a symmetric manner and show similar dwell times. After addition of DNA damage, their static fractions and dwell times decreased, in agreement with increased re-establishment of replication forks. Only a minor fraction of replication forks showed a loss of active polymerases, indicating relatively robust activity during DNA repair. Conversely, PolA, homolog of polymerase I and exonuclease ExoR were rarely present at forks during unperturbed replication but were recruited to replications forks after induction of DNA damage. Protein dynamics of PolA or ExoR were altered in the absence of each other during exponential growth and during DNA repair, indicating overlapping functions. Purified ExoR displayed exonuclease activity and preferentially bound to DNA having 5′ overhangs in vitro. Our analyses support the idea that two replicative DNA polymerases work together at the lagging strand whilst only PolC acts at the leading strand, and that PolA and ExoR perform inducible functions at replication forks during DNA repair.

The fitness landscape of the African Salmonella Typhimurium ST313 strain D23580 reveals unique properties of the pBT1 plasmid

We have used a transposon insertion sequencing (TIS) approach to establish the fitness landscape of the African Salmonella enterica serovar Typhimurium ST313 strain D23580, to complement our previous comparative genomic and functional transcriptomic studies. We used a genome-wide transposon library with insertions every 10 nucleotides to identify genes required for survival and growth in vitro and during infection of murine macrophages. The analysis revealed genomic regions important for fitness under two in vitro growth conditions. Overall, 724 coding genes were required for optimal growth in LB medium, and 851 coding genes were required for growth in SPI-2-inducing minimal medium. These findings were consistent with the essentiality analyses of other S. Typhimurium ST19 and S. Typhi strains. The global mutagenesis approach also identified 60 sRNAs and 413 intergenic regions required for growth in at least one in vitro growth condition. By infecting murine macrophages with the transposon library, we identified 68 genes that were required for intra-macrophage replication but did not impact fitness in vitro. None of these genes were unique to S. Typhimurium D23580, consistent with a high conservation of gene function between S. Typhimurium ST313 and ST19 and suggesting that novel virulence factors are not involved in the interaction of strain D23580 with murine macrophages. We discovered that transposon insertions rarely occurred in many pBT1 plasmid-encoded genes (36), compared with genes carried by the pSLT-BT virulence plasmid and other bacterial plasmids. The key essential protein encoded by pBT1 is a cysteinyl-tRNA synthetase, and our enzymological analysis revealed that the plasmid-encoded CysRSpBT1 had a lower ability to charge tRNA than the chromosomally-encoded CysRSchr enzyme. The presence of aminoacyl-tRNA synthetases in plasmids from a range of Gram-negative and Gram-positive bacteria suggests that plasmid-encoded essential genes are more common than had been appreciated.

ePath: an online database towards comprehensive essential gene annotation for prokaryotes

Experimental techniques for identification of essential genes (EGs) in prokaryotes are usually expensive, time-consuming and sometimes unrealistic. Emerging in silico methods provide alternative methods for EG prediction, but often possess limitations including heavy computational requirements and lack of biological explanation. Here we propose a new computational algorithm for EG prediction in prokaryotes with an online database (ePath) for quick access to the EG prediction results of over 4,000 prokaryotes ( https://www.pubapps.vcu.edu/epath/ ). In ePath, gene essentiality is linked to biological functions annotated by KEGG Ortholog (KO). Two new scoring systems, namely, E_score and P_score, are proposed for each KO as the EG evaluation criteria. E_score represents appearance and essentiality of a given KO in existing experimental results of gene essentiality, while P_score denotes gene essentiality based on the principle that a gene is essential if it plays a role in genetic information processing, cell envelope maintenance or energy production. The new EG prediction algorithm shows prediction accuracy ranging from 75% to 91% based on validation from five new experimental studies on EG identification. Our overall goal with ePath is to provide a comprehensive and reliable reference for gene essentiality annotation, facilitating the study of those prokaryotes without experimentally derived gene essentiality information.

YwqL (EndoV), ExoA and PolA act in a novel alternative excision pathway to repair deaminated DNA bases in Bacillus subtilis

DNA deamination generates base transitions and apurinic/apyrimidinic (AP)-sites which are potentially genotoxic and cytotoxic. In Bacillus subtilis uracil can be removed from DNA by the uracil DNA-glycosylase through the base excision repair pathway. Genetic evidence suggests that B. subtilis YwqL, a homolog of Endonuclease-V (EndoV), acts on a wider spectrum of deaminated bases but the factors that complete this pathway have remained elusive. Here, we report that a purified His₆-YwqL (hereafter BsEndoV) protein had in vitro endonuclease activity against double-stranded DNAs containing a single uracil (U), hypoxanthine (Hx), xanthine (X) or an AP site. Interestingly, while BsEndoV catalyzed a single strand break at the second phosphodiester bond towards the 3'-end of the U and AP lesions, there was an additional cleavage of the phosphodiester bond preceding the Hx and X lesions. Remarkably, the repair event initiated by BsEndoV on Hx and X, was completed by a recombinant B. subtilis His₆-DNA polymerase A (BsPolA), but not on BsEndoV-processed U and AP lesions. For the latter lesions a second excision event performed by a recombinant B. subtilis His₆-ExoA (BsExoA) was necessary before completion of their repair by BsPolA. These results suggest the existence of a novel alternative excision repair pathway in B. subtilis that counteracts the genotoxic effects of base deamination. The presence of this novel pathway in vivo in B. subtilis was also supported by analysis of effects of single or multiple deletions of exoA, endoV and polA on spontaneous mutations in growing cells, and the sensitivity of growing wild-type and mutant cells to a DNA deaminating agent.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

Nitric Oxide Disrupts Zinc Homeostasis in Salmonella enterica Serovar Typhimurium

Nitric oxide (NO·) produced by mammalian cells exerts antimicrobial actions that result primarily from the modification of protein thiols (S-nitrosylation) and metal centers. A comprehensive approach was used to identify novel targets of NO· in Salmonella enterica serovar Typhimurium (S. Typhimurium). Newly identified targets include zinc metalloproteins required for DNA replication and repair (DnaG, PriA, and TopA), protein synthesis (AlaS and RpmE), and various metabolic activities (ClpX, GloB, MetE, PepA, and QueC). The cytotoxic actions of free zinc are mitigated by the ZntA and ZitB zinc efflux transporters, which are required for S. Typhimurium resistance to zinc overload and nitrosative stress in vitro Zinc efflux also ameliorates NO·-dependent zinc mobilization following internalization by activated macrophages and is required for virulence in NO·-producing mice, demonstrating that host-derived NO· causes zinc stress in intracellular bacteria.IMPORTANCE Nitric oxide (NO·) is produced by macrophages in response to inflammatory stimuli and restricts the growth of intracellular bacteria. Mechanisms of NO·-dependent antimicrobial actions are incompletely understood. Here, we show that zinc metalloproteins are important targets of NO· in Salmonella, including the DNA replication proteins DnaG and PriA, which were hypothesized to be NO· targets in earlier studies. Like iron, zinc is a cofactor for several essential proteins but is toxic at elevated concentrations. This study demonstrates that NO· mobilizes free zinc in Salmonella and that specific efflux transporters ameliorate the cytotoxic effects of free zinc during infection.

Emerging and evolving concepts in gene essentiality

Gene essentiality is a founding concept of genetics with important implications in both fundamental and applied research. Multiple screens have been performed over the years in bacteria, yeasts, animals and more recently in human cells to identify essential genes. A mounting body of evidence suggests that gene essentiality, rather than being a static and binary property, is both context dependent and evolvable in all kingdoms of life. This concept of a non-absolute nature of gene essentiality changes our fundamental understanding of essential biological processes and could directly affect future treatment strategies for cancer and infectious diseases.

The Suf Iron-Sulfur Cluster Biosynthetic System Is Essential in Staphylococcus aureus, and Decreased Suf Function Results in Global Metabolic Defects and Reduced Survival in Human Neutrophils

Staphylococcus aureus remains a causative agent for morbidity and mortality worldwide. This is in part a result of antimicrobial resistance, highlighting the need to uncover novel antibiotic targets and to discover new therapeutic agents. In the present study, we explored the possibility that iron-sulfur (Fe-S) cluster synthesis is a viable antimicrobial target. RNA interference studies established that Suf (sulfur mobilization)-dependent Fe-S cluster synthesis is essential in S. aureus We found that sufCDSUB were cotranscribed and that suf transcription was positively influenced by sigma factor B. We characterized an S. aureus strain that contained a transposon inserted in the intergenic space between sufC and sufD (sufD*), resulting in decreased transcription of sufSUB Consistent with the transcriptional data, the sufD* strain had multiple phenotypes associated with impaired Fe-S protein maturation. They included decreased activities of Fe-S cluster-dependent enzymes, decreased growth in media lacking metabolites that require Fe-S proteins for synthesis, and decreased flux through the tricarboxylic acid (TCA) cycle. Decreased Fe-S cluster synthesis resulted in sensitivity to reactive oxygen and reactive nitrogen species, as well as increased DNA damage and impaired DNA repair. The sufD* strain also exhibited perturbed intracellular nonchelated Fe pools. Importantly, the sufD* strain did not exhibit altered exoprotein production or altered biofilm formation, but it was attenuated for survival upon challenge by human polymorphonuclear leukocytes. The results presented are consistent with the hypothesis that Fe-S cluster synthesis is a viable target for antimicrobial development.

Construction of a minimal genome as a chassis for synthetic biology

Microbial diversity and complexity pose challenges in understanding the voluminous genetic information produced from whole-genome sequences, bioinformatics and high-throughput ‘-omics’ research. These challenges can be overcome by a core blueprint of a genome drawn with a minimal gene set, which is essential for life. Systems biology and large-scale gene inactivation studies have estimated the number of essential genes to be ∼300–500 in many microbial genomes. On the basis of the essential gene set information, minimal-genome strains have been generated using sophisticated genome engineering techniques, such as genome reduction and chemical genome synthesis. Current size-reduced genomes are not perfect minimal genomes, but chemically synthesized genomes have just been constructed. Some minimal genomes provide various desirable functions for bioindustry, such as improved genome stability, increased transformation efficacy and improved production of biomaterials. The minimal genome as a chassis genome for synthetic biology can be used to construct custom-designed genomes for various practical and industrial applications.

Bacteriophages and phage-inspired nanocarriers for targeted delivery of therapeutic cargos

The main goal of drug delivery systems is to target therapeutic cargoes to desired cells and to ensure their efficient uptake. Recently a number of studies have focused on designing bio-inspired nanocarriers, such as bacteriophages, and synthetic carriers based on the bacteriophage structure. Bacteriophages are viruses that specifically recognize their bacterial hosts. They can replicate only inside their host cell and can act as natural gene carriers. Each type of phage has a particular shape, a different capacity for loading cargo, a specific production time, and their own mechanisms of supramolecular assembly, that have enabled them to act as tunable carriers. New phage-based technologies have led to the construction of different peptide libraries, and recognition abilities provided by novel targeting ligands. Phage hybridization with non-organic compounds introduces new properties to phages and could be a suitable strategy for construction of bio-inorganic carriers. In this review we try to cover the major phage species that have been used in drug and gene delivery systems, and the biological application of phages as novel targeting ligands and targeted therapeutics.

Analysis of Mycobacterium avium subsp. paratuberculosis mutant libraries reveals loci-dependent transposition biases and strategies for novel mutant discovery

Mycobacterium avium subsp. paratuberculosis (MAP), the aetiological agent of Johne's disease, is one of the most important bacterial pathogens in ruminants. A thorough understanding of MAP pathogenesis is needed to develop new vaccines and diagnostic tests. The generation of comprehensive random transposon mutant libraries is a fundamental genetic technology to determine the role of genes in physiology and pathogenesis. In this study, whole MAP genome analysis compared the insertion sites for the mycobacterial transposon Tn5367 derived from the Mycobacterium smegmatis insertion sequence IS1096 and the mariner transposon MycoMarT7 carrying the Himar1 transposase. We determined that only MycoMarT7 provides a random representation of insertions in 99 % of all MAP genes. Analysis of the MAP K-10 genome indicated that 710 of all ORFs do not possess IS1096 recognition sites, while only 37 do not have the recognition site for MycoMarT7. Thus, a significant number of MAP genes remain underrepresented in insertion libraries from IS1096-derived transposons. Analysis of MycoMarT7 and Tn5367 mutants showed that Tn5367 has a predilection to insert within intergenic regions, suggesting that MycoMarT7 is the more adequate for generating a comprehensive library. However, we uncovered the novel finding that both transposons have loci-dependent biases, with Tn5367 being the most skewed. These loci-dependent transposition biases led to an underestimation of the number of independent mutants required to generate a comprehensive mutant library, leading to an overestimation of essential genes. Herein, we also demonstrated a useful platform for gene discovery and analysis by isolating three novel mutants for each transposon.

Homology modeling, functional annotation and comparative genomics of outer membrane protein H of Pasteurella multocida

Pasteurella multocida is an important pathogen of animals and humans. Outer Membrane Protein (Omp) H is a major conserved protein in the envelope of P. multocida and has been commonly targeted as a protective antigen. However, not much is known about its structure and function due to the difficulties that are typically associated with obtaining sufficient amounts of purified prokaryotic transmembrane proteins. The present work is aimed at studying the OmpH using an in silico approach and consolidate the findings in light of existing experimental evidences. Our study describes the first 3D model of the P. multocida OmpH obtained through a combination of several in silico modeling approaches. From our results, OmpH of P. multocida could be classified as a homotrimeric, 16 stranded, β-barrel porin involved in the non-specific transport of small, hydrophilic molecules, serving essential osmoregulatory function. Moreover, very small homologous sequences could be identified in the host proteome, strengthening the probability of a successful OmpH-based vaccine against the pathogen with remote chances of cross-reaction to host proteins.

Systems biology perspectives on minimal and simpler cells

The concept of the minimal cell has fascinated scientists for a long time, from both fundamental and applied points of view. This broad concept encompasses extreme reductions of genomes, the last universal common ancestor (LUCA), the creation of semiartificial cells, and the design of protocells and chassis cells. Here we review these different areas of research and identify common and complementary aspects of each one. We focus on systems biology, a discipline that is greatly facilitating the classical top-down and bottom-up approaches toward minimal cells. In addition, we also review the so-called middle-out approach and its contributions to the field with mathematical and computational models. Owing to the advances in genomics technologies, much of the work in this area has been centered on minimal genomes, or rather minimal gene sets, required to sustain life. Nevertheless, a fundamental expansion has been taking place in the last few years wherein the minimal gene set is viewed as a backbone of a more complex system. Complementing genomics, progress is being made in understanding the system-wide properties at the levels of the transcriptome, proteome, and metabolome. Network modeling approaches are enabling the integration of these different omics data sets toward an understanding of the complex molecular pathways connecting genotype to phenotype. We review key concepts central to the mapping and modeling of this complexity, which is at the heart of research on minimal cells. Finally, we discuss the distinction between minimizing the number of cellular components and minimizing cellular complexity, toward an improved understanding and utilization of minimal and simpler cells.

Gas chromatography-mass spectrometry-based metabolite profiling of Salmonella enterica serovar Typhimurium differentiates between biofilm and planktonic phenotypes

The aim of this study was to utilize gas chromatography coupled with mass spectrometry (GC-MS) to compare and identify patterns of biochemical change between Salmonella cells grown in planktonic and biofilm phases and Salmonella biofilms of different ages. Our results showed a clear separation between planktonic and biofilm modes of growth. The majority of metabolites contributing to variance between planktonic and biofilm supernatants were identified as amino acids, including alanine, glutamic acid, glycine, and ornithine. Metabolites contributing to variance in intracellular profiles were identified as succinic acid, putrescine, pyroglutamic acid, and N-acetylglutamic acid. Principal-component analysis revealed no significant differences between the various ages of intracellular profiles, which would otherwise allow differentiation of biofilm cells on the basis of age. A shifting pattern across the score plot was illustrated when analyzing extracellular metabolites sampled from different days of biofilm growth, and amino acids were again identified as the metabolites contributing most to variance. An understanding of biofilm-specific metabolic responses to perturbations, especially antibiotics, can lead to the identification of novel drug targets and potential therapies for combating biofilm-associated diseases. We concluded that under the conditions of this study, GC-MS can be successfully applied as a high-throughput technique for "bottom-up" metabolomic biofilm research.

Gene essentiality analysis based on DEG 10, an updated database of essential genes

The database of essential genes (DEG, available at http://www.essentialgene.org), constructed in 2003, has been timely updated to harbor essential-gene records of bacteria, archaea, and eukaryotes. DEG 10, the current release, includes not only essential protein-coding genes determined by genome-wide gene essentiality screens but also essential noncoding RNAs, promoters, regulatory sequences, and replication origins. Therefore, DEG 10 includes essential genomic elements under different conditions in three domains of life, with customizable BLAST tools. Based on the analysis of DEG 10, we show that the percentage of essential genes in bacterial genomes exhibits an exponential decay with increasing genome sizes. The functions, ATP binding (GO:0005524), GTP binding (GO:0005525), and DNA-directed RNA polymerase activity (GO:0003899), are likely required for organisms across life domains.

Microarray transposon tracking for the mapping of conditionally essential genes in Campylobacter jejuni

Although whole genome approaches to the study of bacteria have grown substantially in the past decade, there is still a need for quick and easy methods for the determination of which genes are essential for the growth of these bacteria under specific growth conditions. There are numerous methods to accomplish this depending on the resources and equipment available, each with their own advantages and disadvantages. Here we describe a method we successfully employed to map the essential genes of Campylobacter jejuni using a microarray transposon tracking approach where we constructed a saturated transposon mutant library in the C. jejuni strain NCTC11168 and used a genomic microarray approach to identify genes lacking transposon insertions under standard laboratory growth conditions. With a fully saturated library, the absence of transposon insertions can be used as an indicator of a gene essential for the survival and growth for the conditions used for the mutant library.

Energetics of pathogenic bacteria and opportunities for drug development

The emergence and spread of drug-resistant pathogens and our inability to develop new antimicrobials to overcome resistance has inspired scientists to consider new targets for drug development. Cellular bioenergetics is an area showing promise for the development of new antimicrobials, particularly in the discovery of new anti-tuberculosis drugs where several new compounds have entered clinical trials. In this review, we have examined the bioenergetics of various bacterial pathogens, highlighting the versatility of electron donor and acceptor utilisation and the modularity of electron transport chain components in bacteria. In addition to re-examining classical concepts, we explore new literature that reveals the intricacies of pathogen energetics, for example, how Salmonella enterica and Campylobacter jejuni exploit host and microbiota to derive powerful electron donors and sinks; the strategies Mycobacterium tuberculosis and Pseudomonas aeruginosa use to persist in lung tissues; and the importance of sodium energetics and electron bifurcation in the chemiosmotic anaerobe Fusobacterium nucleatum. A combination of physiological, biochemical, and pharmacological data suggests that, in addition to the clinically-approved target F1Fo-ATP synthase, NADH dehydrogenase type II, succinate dehydrogenase, hydrogenase, cytochrome bd oxidase, and menaquinone biosynthesis pathways are particularly promising next-generation drug targets. The realisation of cellular energetics as a rich target space for the development of new antimicrobials will be dependent upon gaining increased understanding of the energetic processes utilised by pathogens in host environments and the ability to design bacterial-specific inhibitors of these processes.

Stress response of Salmonella enterica serovar typhimurium to acidified nitrite

The antimicrobial action of the curing agent sodium nitrite (NaNO2), which is added as a preservative to raw meat products, depends on its conversion to nitric oxide and other reactive nitrogen species under acidic conditions. In this study, we used RNA sequencing to analyze the acidified-NaNO2 shock and adaptive responses of Salmonella enterica serovar Typhimurium, a frequent contaminant in raw meat, considering parameters relevant for the production of raw-cured sausages. Upon a 10-min exposure to 150 mg/liter NaNO2 in LB (pH 5.5) acidified with lactic acid, genes involved in nitrosative-stress protection, together with several other stress-related genes, were induced. In contrast, genes involved in translation, transcription, replication, and motility were downregulated. The induction of stress tolerance and the reduction of cell proliferation obviously promote survival under harsh acidified-NaNO2 stress. The subsequent adaptive response was characterized by upregulation of NsrR-regulated genes and iron uptake systems and by downregulation of genes involved in anaerobic respiratory pathways. Strikingly, amino acid decarboxylase systems, which contribute to acid tolerance, displayed increased transcript levels in response to acidified NaNO2. The induction of systems known to be involved in acid resistance indicates a nitrite-mediated increase in the level of acid stress. Deletion of cadA, which encodes lysine decarboxylase, resulted in increased sensitivity to acidified NaNO2. Intracellular pH measurements using a pH-sensitive green fluorescent protein (GFP) variant showed that the cytoplasmic pH of S. Typhimurium in LB medium (pH 5.5) is decreased upon the addition of NaNO2. This study provides the first evidence that intracellular acidification is an additional antibacterial mode of action of acidified NaNO2.

Species-selective killing of bacteria by antimicrobial peptide-PNAs

Broad-spectrum antimicrobials kill indiscriminately, a property that can lead to negative clinical consequences and an increase in the incidence of resistance. Species-specific antimicrobials that could selectively kill pathogenic bacteria without targeting other species in the microbiome could limit these problems. The pathogen genome presents an excellent target for the development of such antimicrobials. In this study we report the design and evaluation of species-selective peptide nucleic acid (PNA) antibacterials. Selective growth inhibition of B. subtilis, E. coli, K. pnuemoniae and S. enterica serovar Typhimurium in axenic or mixed culture could be achieved with PNAs that exploit species differences in the translation initiation region of essential genes. An S. Typhimurium-specific PNA targeting ftsZ resulted in elongated cells that were not observed in E. coli, providing phenotypic evidence of the selectivity of PNA-based antimicrobials. Analysis of the genomes of E. coli and S. Typhimurium gave a conservative estimate of >150 PNA targets that could potentially discriminate between these two closely related species. This work provides a basis for the development of a new class of antimicrobial with a tuneable spectrum of activity.

Non-essential genes form the hubs of genome scale protein function and environmental gene expression networks in Salmonella enterica serovar Typhimurium

BACKGROUND

Salmonella Typhimurium is an important pathogen of human and animals. It shows a broad growth range and survives in harsh conditions. The aim of this study was to analyze transcriptional responses to a number of growth and stress conditions as well as the relationship of metabolic pathways and/or cell functions at the genome-scale-level by network analysis, and further to explore whether highly connected genes (hubs) in these networks were essential for growth, stress adaptation and virulence.

RESULTS

De novo generated as well as published transcriptional data for 425 selected genes under a number of growth and stress conditions were used to construct a bipartite network connecting culture conditions and significantly regulated genes (transcriptional network). Also, a genome scale network was constructed for strain LT2. The latter connected genes with metabolic pathways and cellular functions. Both networks were shown to belong to the family of scale-free networks characterized by the presence of highly connected nodes or hubs which are genes whose transcription is regulated when responding to many of the assayed culture conditions or genes encoding products involved in a high number of metabolic pathways and cell functions.The five genes with most connections in the transcriptional network (wraB, ygaU, uspA, cbpA and osmC) and in the genome scale network (ychN, siiF (STM4262), yajD, ybeB and dcoC) were selected for mutations, however mutagenesis of ygaU and ybeB proved unsuccessful. No difference between mutants and the wild type strain was observed during growth at unfavorable temperatures, pH values, NaCl concentrations and in the presence of H2O2. Eight mutants were evaluated for virulence in C57/BL6 mice and none differed from the wild type strain. Notably, however, deviations of phenotypes with respect to the wild type were observed when combinations of these genes were deleted.

CONCLUSION

Network analysis revealed the presence of hubs in both transcriptional and functional networks of S. Typhimurium. Hubs theoretically confer higher resistance to random mutation but a greater susceptibility to directed attacks, however, we found that genes that formed hubs were dispensable for growth, stress adaptation and virulence, suggesting that evolution favors non-essential genes as main connectors in cellular networks.

DEG 10, an update of the database of essential genes that includes both protein-coding genes and noncoding genomic elements

The combination of high-density transposon-mediated mutagenesis and high-throughput sequencing has led to significant advancements in research on essential genes, resulting in a dramatic increase in the number of identified prokaryotic essential genes under diverse conditions and a revised essential-gene concept that includes all essential genomic elements, rather than focusing on protein-coding genes only. DEG 10, a new release of the Database of Essential Genes (available at http://www.essentialgene.org), has been developed to accommodate these quantitative and qualitative advancements. In addition to increasing the number of bacterial and archaeal essential genes determined by genome-wide gene essentiality screens, DEG 10 also harbors essential noncoding RNAs, promoters, regulatory sequences and replication origins. These essential genomic elements are determined not only in vitro, but also in vivo, under diverse conditions including those for survival, pathogenesis and antibiotic resistance. We have developed customizable BLAST tools that allow users to perform species- and experiment-specific BLAST searches for a single gene, a list of genes, annotated or unannotated genomes. Therefore, DEG 10 includes essential genomic elements under different conditions in three domains of life, with customizable BLAST tools.

Comprehensive assignment of roles for Salmonella typhimurium genes in intestinal colonization of food-producing animals

Chickens, pigs, and cattle are key reservoirs of Salmonella enterica, a foodborne pathogen of worldwide importance. Though a decade has elapsed since publication of the first Salmonella genome, thousands of genes remain of hypothetical or unknown function, and the basis of colonization of reservoir hosts is ill-defined. Moreover, previous surveys of the role of Salmonella genes in vivo have focused on systemic virulence in murine typhoid models, and the genetic basis of intestinal persistence and thus zoonotic transmission have received little study. We therefore screened pools of random insertion mutants of S. enterica serovar Typhimurium in chickens, pigs, and cattle by transposon-directed insertion-site sequencing (TraDIS). The identity and relative fitness in each host of 7,702 mutants was simultaneously assigned by massively parallel sequencing of transposon-flanking regions. Phenotypes were assigned to 2,715 different genes, providing a phenotype–genotype map of unprecedented resolution. The data are self-consistent in that multiple independent mutations in a given gene or pathway were observed to exert a similar fitness cost. Phenotypes were further validated by screening defined null mutants in chickens. Our data indicate that a core set of genes is required for infection of all three host species, and smaller sets of genes may mediate persistence in specific hosts. By assigning roles to thousands of Salmonella genes in key reservoir hosts, our data facilitate systems approaches to understand pathogenesis and the rational design of novel cross-protective vaccines and inhibitors. Moreover, by simultaneously assigning the genotype and phenotype of over 90% of mutants screened in complex pools, our data establish TraDIS as a powerful tool to apply rich functional annotation to microbial genomes with minimal animal use.

Salmonella Typhimurium is a major cause of human diarrhoeal infections, usually acquired from chickens, pigs, cattle, or their products. To understand the basis of persistence and pathogenesis in these reservoir hosts, and to inform the design of novel vaccines and treatments, we generated a library of 7,702 S. Typhimurium mutants, each bearing an insertion at a random position in the genome. Using DNA sequencing, we identified the disrupted gene in each mutant and determined its relative abundance in a laboratory culture and after experimental infection of mice, chickens, pigs, and cattle. The method allowed large numbers of mutants to be investigated simultaneously, drastically reducing the number of animals required to perform a comprehensive screen. We identified mutants that grow in culture but do not survive in one or more of the animals. The genes disrupted in these mutants are inferred to be important for the infection process. Most of these genes were required in all three food-producing animals, but smaller subsets of genes may mediate persistence in a specific host species. The data provide the most comprehensive map of virulence-associated genes for any bacterial pathogen in natural hosts and are highly relevant for the design of control strategies.

Underlying functional genomics of fat deposition in adipose tissue

The objective of this study was to gain insight into the underlying mechanisms of fat deposition. Two sheep breeds with large fat-tail (Lori-Bakhtiari) and with thin-tail (Zel) were used as models. To determine important and key candidate lipid metabolism related genes, comparative genomic approaches were employed. Gene expression profiles of adipose tissues were analyzed in human, pig, and cattle by express sequence tag (EST) analysis. EST analysis determined 65, 102 and 125 transcripts in human, pig and cattle respectively with at least 10 fold over-expression in the adipose tissue. Based on our comparative functional genomic analysis, seven genes were more abundant and common in investigated mammalian adipose tissues promising a conserved novel gene network in mammalian lipid metabolism. The candidate genes including fatty acid binding protein 4 (FABP4), fatty acid synthase (FASN), Stearoyl-CoA desaturase (SCD) and Lipoprotein lipase (LPL) were selected for further gene expression investigation within two sheep breeds. The real time PCR results showed that among the genes tested, FABP4 was expressed at higher levels than the others. The expression of FABP4 was significantly higher in the fat-tail of Lori-Bakhtiari than in the fat-tail and visceral adipose tissues of Zel (P<0.05). The findings suggest that the FABP4 gene expression in the fat-tail is an important index of fat deposition.

A comparison of dense transposon insertion libraries in the Salmonella serovars Typhi and Typhimurium

Salmonella Typhi and Typhimurium diverged only ∼50 000 years ago, yet have very different host ranges and pathogenicity. Despite the availability of multiple whole-genome sequences, the genetic differences that have driven these changes in phenotype are only beginning to be understood. In this study, we use transposon-directed insertion-site sequencing to probe differences in gene requirements for competitive growth in rich media between these two closely related serovars. We identify a conserved core of 281 genes that are required for growth in both serovars, 228 of which are essential in Escherichia coli. We are able to identify active prophage elements through the requirement for their repressors. We also find distinct differences in requirements for genes involved in cell surface structure biogenesis and iron utilization. Finally, we demonstrate that transposon-directed insertion-site sequencing is not only applicable to the protein-coding content of the cell but also has sufficient resolution to generate hypotheses regarding the functions of non-coding RNAs (ncRNAs) as well. We are able to assign probable functions to a number of cis-regulatory ncRNA elements, as well as to infer likely differences in trans-acting ncRNA regulatory networks.

Expression divergence between Escherichia coli and Salmonella enterica serovar Typhimurium reflects their lifestyles

Escherichia coli K12 is a commensal bacteria and one of the best-studied model organisms. Salmonella enterica serovar Typhimurium, on the other hand, is a facultative intracellular pathogen. These two prokaryotic species can be considered related phylogenetically, and they share a large amount of their genetic material, which is commonly termed the "core genome." Despite their shared core genome, both species display very different lifestyles, and it is unclear to what extent the core genome, apart from the species-specific genes, plays a role in this lifestyle divergence. In this study, we focus on the differences in expression domains for the orthologous genes in E. coli and S. Typhimurium. The iterative comparison of coexpression methodology was used on large expression compendia of both species to uncover the conservation and divergence of gene expression. We found that gene expression conservation occurs mostly independently from amino acid similarity. According to our estimates, at least more than one quarter of the orthologous genes has a different expression domain in E. coli than in S. Typhimurium. Genes involved with key cellular processes are most likely to have conserved their expression domains, whereas genes showing diverged expression are associated with metabolic processes that, although present in both species, are regulated differently. The expression domains of the shared "core" genome of E. coli and S. Typhimurium, consisting of highly conserved orthologs, have been tuned to help accommodate the differences in lifestyle and the pathogenic potential of Salmonella.

From essential to persistent genes: a functional approach to constructing synthetic life

A central undertaking in synthetic biology (SB) is the quest for the 'minimal genome'. However, 'minimal sets' of essential genes are strongly context-dependent and, in all prokaryotic genomes sequenced to date, not a single protein-coding gene is entirely conserved. Furthermore, a lack of consensus in the field as to what attributes make a gene truly essential adds another aspect of variation. Thus, a universal minimal genome remains elusive. Here, as an alternative to defining a minimal genome, we propose that the concept of gene persistence can be used to classify genes needed for robust long-term survival. Persistent genes, although not ubiquitous, are conserved in a majority of genomes, tend to be expressed at high levels, and are frequently located on the leading DNA strand. These criteria impose constraints on genome organization, and these are important considerations for engineering cells and for creating cellular life-like forms in SB.

Global transcriptional analysis of dehydrated Salmonella enterica serovar Typhimurium

Despite the scientific and industrial importance of desiccation tolerance in Salmonella, knowledge regarding its genetic basis is still scarce. In the present study, we performed a transcriptomic analysis of dehydrated and water-suspended Salmonella enterica serovar Typhimurium using microarrays. Dehydration induced expression of 90 genes and downregulated that of 7 genes. Ribosomal structural genes represented the most abundant functional group with a relatively higher transcription during dehydration. Other main induced functional groups included genes involved in amino acid metabolism, energy production, ion transport, transcription, and stress response. The highest induction was observed in the kdpFABC operon, encoding a potassium transport channel. Knockout mutations were generated in nine upregulated genes. Five mutants displayed lower tolerance to desiccation, implying the involvement of the corresponding genes in the adaptation of Salmonella to desiccation. These included genes encoding the isocitrate-lyase AceA, the lipid A biosynthesis palmitoleoyl-acyltransferase Ddg, the modular iron-sulfur cluster scaffolding protein NifU, the global regulator Fnr, and the alternative sigma factor RpoE. Notably, these proteins were previously implicated in the response of Salmonella to oxidative stress, heat shock, and cold shock. A strain with a mutation in the structural gene kdpA had a tolerance to dehydration comparable to that of the parent strain, implying that potassium transport through this system is dispensable for early adaptation to the dry environment. Nevertheless, this mutant was significantly impaired in long-term persistence during cold storage. Our findings indicate the involvement of a relatively small fraction of the Salmonella genome in transcriptional adjustment from water to dehydration, with a high prevalence of genes belonging to the protein biosynthesis machinery.

Identification of essential genes of the periodontal pathogen Porphyromonas gingivalis

Background

Porphyromonas gingivalis is a Gram-negative anaerobic bacterium associated with periodontal disease onset and progression. Genetic tools for the manipulation of bacterial genomes allow for in-depth mechanistic studies of metabolism, physiology, interspecies and host-pathogen interactions. Analysis of the essential genes, protein-coding sequences necessary for survival of P. gingivalis by transposon mutagenesis has not previously been attempted due to the limitations of available transposon systems for the organism. We adapted a Mariner transposon system for mutagenesis of P. gingivalis and created an insertion mutant library. By analyzing the location of insertions using massively-parallel sequencing technology we used this mutant library to define genes essential for P. gingivalis survival under in vitro conditions.

Results

In mutagenesis experiments we identified 463 genes in P. gingivalis strain ATCC 33277 that are putatively essential for viability in vitro. Comparing the 463 P. gingivalis essential genes with previous essential gene studies, 364 of the 463 are homologues to essential genes in other species; 339 are shared with more than one other species. Twenty-five genes are known to be essential in P. gingivalis and B. thetaiotaomicron only. Significant enrichment of essential genes within Cluster of Orthologous Groups ‘D’ (cell division), ‘I’ (lipid transport and metabolism) and ‘J’ (translation/ribosome) were identified. Previously, the P. gingivalis core genome was shown to encode 1,476 proteins out of a possible 1,909; 434 of 463 essential genes are contained within the core genome. Thus, for the species P. gingivalis twenty-two, seventy-seven and twenty-three percent of the genome respectively are devoted to essential, core and accessory functions.

Conclusions

A Mariner transposon system can be adapted to create mutant libraries in P. gingivalis amenable to analysis by next-generation sequencing technologies. In silico analysis of genes essential for in vitro growth demonstrates that although the majority are homologous across bacterial species as a whole, species and strain-specific subsets are apparent. Understanding the putative essential genes of P. gingivalis will provide insights into metabolic pathways and niche adaptations as well as clinical therapeutic strategies.

High-throughput comparison of gene fitness among related bacteria

Background

The contribution of a gene to the fitness of a bacterium can be assayed by whether and to what degree the bacterium tolerates transposon insertions in that gene. We use this fact to compare the fitness of syntenic homologous genes among related Salmonella strains and thereby reveal differences not apparent at the gene sequence level.

Results

A transposon Tn5 derivative was used to construct mutants in Salmonella Typhimurium ATCC14028 (STM1) and Salmonella Typhi Ty2 (STY1), which were then grown in rich media. The locations of 234,152 and 53,556 integration sites, respectively, were mapped by sequencing. These data were compared to similar data available for a different Ty2 isolate (STY2) and essential genes identified in E. coli K-12 (ECO). Of 277 genes considered essential in ECO, all had syntenic homologs in STM1, STY1, and STY2, and all but nine genes were either devoid of transposon insertions or had very few. For three of these nine genes, part of the annotated gene lacked transposon integrations (yejM, ftsN and murB). At least one of the other six genes, trpS, had a potentially functionally redundant gene encoded elsewhere in Salmonella but not in ECO. An additional 165 genes were almost entirely devoid of transposon integrations in all three Salmonella strains examined, including many genes associated with protein and DNA synthesis. Four of these genes (STM14_1498, STM14_2872, STM14_3360, and STM14_5442) are not found in E. coli. Notable differences in the extent of gene selection were also observed among the three different Salmonella isolates. Mutations in hns, for example, were selected against in STM1 but not in the two STY strains, which have a defect in rpoS rendering hns nonessential.

Conclusions

Comparisons among transposon integration profiles from different members of a species and among related species, all grown in similar conditions, identify differences in gene contributions to fitness among syntenic homologs. Further differences in fitness profiles among shared genes can be expected in other selective environments, with potential relevance for comparative systems biology.

Superessential reactions in metabolic networks

The metabolic genotype of an organism can change through loss and acquisition of enzyme-coding genes, while preserving its ability to survive and synthesize biomass in specific environments. This evolutionary plasticity allows pathogens to evolve resistance to antimetabolic drugs by acquiring new metabolic pathways that bypass an enzyme blocked by a drug. We here study quantitatively the extent to which individual metabolic reactions and enzymes can be bypassed. To this end, we use a recently developed computational approach to create large metabolic network ensembles that can synthesize all biomass components in a given environment but contain an otherwise random set of known biochemical reactions. Using this approach, we identify a small connected core of 124 reactions that are absolutely superessential (that is, required in all metabolic networks). Many of these reactions have been experimentally confirmed as essential in different organisms. We also report a superessentiality index for thousands of reactions. This index indicates how easily a reaction can be bypassed. We find that it correlates with the number of sequenced genomes that encode an enzyme for the reaction. Superessentiality can help choose an enzyme as a potential drug target, especially because the index is not highly sensitive to the chemical environment that a pathogen requires. Our work also shows how analyses of large network ensembles can help understand the evolution of complex and robust metabolic networks.

Genetic analysis of 15 protein folding factors and proteases of the Escherichia coli cell envelope

Each cell hosts thousands of proteins that vary greatly in abundance, structure, and chemical properties. To ensure that all proteins are biologically active and properly localized, efficient quality control systems have evolved. While the structure, function, and regulation of some individual protein folding factors and proteases were resolved up to atomic resolution, others remain poorly characterized. In addition, little is known about which factors are required for viability under specific stress conditions. We therefore determined the physiological implications of 15 factors of the E. coli cell envelope by an integrated genetic approach comprising phenotypic analyses. Our data indicate that surA and tsp null mutations are a lethal combination in rich medium, that surA dsbA and surA dsbC double mutants are temperature sensitive, and that surA ptrA, surA yfgC, dsbA fkpA, degP tsp, degP ppiD, tsp ppiD, and degP dsbA double mutants are temperature sensitive in rich medium containing 0.5 M NaCl, while degP dsbA, degP yfgC, tsp ydgD, and degP tsp double mutants do not grow in the presence of SDS/EDTA. Furthermore, we show that in degP dsbA, degP tsp, and degP yfgC double mutants a subpopulation of LamB exists as unfolded monomers. In addition, dsbA null mutants expressed lower levels of the outer membrane proteins LptD, LamB, FhuA, and OmpW while FhuA levels were reduced in surA single and degP ppiD double mutants. Lower FhuA levels in degP ppiD strains depend on Tsp, since in a tsp degP ppiD triple mutant FhuA levels are restored.

[Current status of theoretical studies on essential genes in microbes]

Essential genes are indispensable for the survival of an organism in optimal conditions. Recently, study on essential gene is becoming a hot topic of microbiology, genomics, and bioinformatics. This paper described the experiments that determined essential genes in some microbes and the theoretical researches on essential genes were reviewed. The major content contained comparison of essential genes and non-essential genes based on information on evolutionary conservation and sequence composition, and in silico prediction of essential genes, and analysis of the chromosomal distributions of essential genes. Finally, related progresses were concluded and the open problems were pointed out.

OGEE: an online gene essentiality database

OGEE is an Online GEne Essentiality database. Its main purpose is to enhance our understanding of the essentiality of genes. This is achieved by collecting not only experimentally tested essential and non-essential genes, but also associated gene features such as expression profiles, duplication status, conservation across species, evolutionary origins and involvement in embryonic development. We focus on large-scale experiments and complement our data with text-mining results. Genes are organized into data sets according to their sources. Genes with variable essentiality status across data sets are tagged as conditionally essential, highlighting the complex interplay between gene functions and environments. Linked tools allow the user to compare gene essentiality among different gene groups, or compare features of essential genes to non-essential genes, and visualize the results. OGEE is freely available at http://ogeedb.embl.de.

Genome-wide essential gene identification in Streptococcus sanguinis

A clear perception of gene essentiality in bacterial pathogens is pivotal for identifying drug targets to combat emergence of new pathogens and antibiotic-resistant bacteria, for synthetic biology, and for understanding the origins of life. We have constructed a comprehensive set of deletion mutants and systematically identified a clearly defined set of essential genes for Streptococcus sanguinis. Our results were confirmed by growing S. sanguinis in minimal medium and by double-knockout of paralogous or isozyme genes. Careful examination revealed that these essential genes were associated with only three basic categories of biological functions: maintenance of the cell envelope, energy production, and processing of genetic information. Our finding was subsequently validated in two other pathogenic streptococcal species, Streptococcus pneumoniae and Streptococcus mutans and in two other gram-positive pathogens, Bacillus subtilis and Staphylococcus aureus. Our analysis has thus led to a simplified model that permits reliable prediction of gene essentiality.

Enzymes are enriched in bacterial essential genes

Essential genes, those indispensable for the survival of an organism, play a key role in the emerging field, synthetic biology. Characterization of functions encoded by essential genes not only has important practical implications, such as in identifying antibiotic drug targets, but can also enhance our understanding of basic biology, such as functions needed to support cellular life. Enzymes are critical for almost all cellular activities. However, essential genes have not been systematically examined from the aspect of enzymes and the chemical reactions that they catalyze. Here, by comprehensively analyzing essential genes in 14 bacterial genomes in which large-scale gene essentiality screens have been performed, we found that enzymes are enriched in essential genes. Essential enzymes have overrepresented ligases (especially those forming carbon-oxygen bonds and carbon-nitrogen bonds), nucleotidyltransferases and phosphotransferases, while have underrepresented oxidoreductases. Furthermore, essential enzymes tend to associate with more gene ontology domains. These results, from the aspect of chemical reactions, provide further insights into the understanding of functions needed to support natural cellular life, as well as synthetic cells, and provide additional parameters that can be integrated into gene essentiality prediction algorithms.

Human body temperature and new approaches to constructing temperature-sensitive bacterial vaccines

Many of the live human and animal vaccines that are currently in use are attenuated by virtue of their temperature-sensitive (TS) replication. These vaccines are able to function because they can take advantage of sites in mammalian bodies that are cooler than the core temperature, where TS vaccines fail to replicate. In this article, we discuss the distribution of temperature in the human body, and relate how the temperature differential can be exploited for designing and using TS vaccines. We also examine how one of the coolest organs of the body, the skin, contains antigen-processing cells that can be targeted to provoke the desired immune response from a TS vaccine. We describe traditional approaches to making TS vaccines, and highlight new information and technologies that are being used to create a new generation of engineered TS vaccines. We pay particular attention to the recently described technology of substituting essential genes from Arctic bacteria for their homologues in mammalian pathogens as a way of creating TS vaccines.

Gene dispensability

Genome-wide mutagenesis studies indicate that up to about 90% of genes in bacteria and 80% in eukaryotes can be inactivated individually leaving an organism viable, often seemingly unaffected. Several strategies are used to learn what these apparently dispensable genes contribute to fitness. Assays of growth under hundreds of physical and chemical stresses are among the most effective experimental approaches. Comparative studies of genomic DNA sequences continue to be valuable in discriminating between the core bacterial genome and the more variable niche-specific genes. The concept of the core genome appears currently unfeasible for eukaryotes but progress has been made in understanding why they contain numerous gene duplicates.

Identification of essential genes in C. jejuni genome highlights hyper-variable plasticity regions

A microarray transposon-based tracking approach was used to identify Campylobacter jejuni genes which are required for cell growth at 37°C, under a microaerophilic atmosphere and on a rich Mueller-Hinton medium. A transposon-based mutant library, comprised of 7,201 individual mutants was constructed, representing 4.48× coverage of the genome. An analysis of genes lacking a transposon insertion revealed 195 essential gene candidates. The function of these genes represent many of the expected core functions of the cell, such as energy metabolism, macromolecule and cofactor biosynthesis, cell structural proteins as well as basic cell processes. Forty-nine hypothetical proteins were also identified, further underlining the importance of currently unknown proteins and pathways within C. jejuni. Unlike other bacteria, the essential genes were not uniformly distributed along the chromosome with three main regions lacking essential genes. These particular regions corresponded to known hyper-variable plasticity regions of C. jejuni genome indicating, as expected, that these regions are dispensable in any given C. jejuni strain. Overall, this work identified dispensable and essential genes in C. jejuni that will ultimately lead to a better understanding of Campylobacter physiology.

Temperature-sensitive bacterial pathogens generated by the substitution of essential genes from cold-loving bacteria: potential use as live vaccines

Temperature-sensitive (TS) viruses have been used for decades as vaccines capable of limited replication in their hosts. Although attenuated bacteria, such as the Bacille Calmette-Guérin anti-tuberculosis vaccine, have been used for almost a century, it is only recently that there has been progress in using TS bacterial strains as live vaccines. Decades of work on essential bacterial genes and the recent explosion in the number of available bacterial genomic sequences set the groundwork for the identification of essential genes from diverse bacteria. This knowledge has allowed for the substitution of essential genes from cold-loving bacteria into the chromosomes of pathogenic bacteria. Many of these gene substitutions generated TS pathogenic bacterial strains, and some were demonstrated to provide protective immunity in mice. This work opens the possibility of engineering many pathogenic bacteria to create TS strains that can be used as vaccines.

Building the blueprint of life

With recent breakthroughs in experimental microbiology making it possible to synthesize and implant an entire genome to create a living cell, the challenge of constructing a working blueprint for the first truly minimal synthetic organism is more important than ever. Here we review the significant progress made in the design and creation of a minimal organism. We discuss how comparative genomes, gene essentiality data, naturally small genomes, and metabolic modeling are all being applied to produce a catalogue of the biological functions essential for life. We compare the minimal gene sets from three published sources with functions identified in 13 existing gene essentiality datasets. We examine how genome-scale metabolic models have been applied to design a minimal metabolism for growth in simple and complex media. Additionally, we survey the progress of efforts to construct a minimal organism, either through implementation of combinatorial deletions in Bacillus subtilis and Escherichia coli or through the synthesis and implantation of synthetic genomes.

Deciphering the intracellular metabolism of Listeria monocytogenes by mutant screening and modelling

Background

The human pathogen Listeria monocytogenes resides and proliferates within the cytoplasm of epithelial cells. While the virulence factors essentially contributing to this step of the infection cycle are well characterized, the set of listerial genes contributing to intracellular replication remains to be defined on a genome-wide level.

Results

A comprehensive library of L. monocytogenes strain EGD knockout mutants was constructed upon insertion-duplication mutagenesis, and 1491 mutants were tested for their phenotypes in rich medium and in a Caco-2 cell culture assay. Following sequencing of the plasmid insertion site, 141 different genes required for invasion of and replication in Caco-2 cells were identified. Ten in-frame deletion mutants were constructed that confirmed the data. The genes with known functions are mainly involved in cellular processes including transport, in the intermediary metabolism of sugars, nucleotides and lipids, and in information pathways such as regulatory functions. No function could be ascribed to 18 genes, and a counterpart of eight genes is missing in the apathogenic species L. innocua. Mice infection studies revealed the in vivo requirement of IspE (Lmo0190) involved in mevalonate synthesis, and of the novel ABC transporter Lmo0135-0137 associated with cysteine transport. Based on the data of this genome-scale screening, an extreme pathway and elementary mode analysis was applied that demonstrates the critical role of glycerol and purine metabolism, of fucose utilization, and of the synthesis of glutathione, aspartate semialdehyde, serine and branched chain amino acids during intracellular replication of L. monocytogenes.

Conclusion

The combination of a genetic screening and a modelling approach revealed that a series of transporters help L. monocytogenes to overcome a putative lack of nutrients within cells, and that a high metabolic flexibility contributes to the intracellular replication of this pathogen.

Identifying essential genes in bacterial metabolic networks with machine learning methods

Background

Identifying essential genes in bacteria supports to identify potential drug targets and an understanding of minimal requirements for a synthetic cell. However, experimentally assaying the essentiality of their coding genes is resource intensive and not feasible for all bacterial organisms, in particular if they are infective.

Results

We developed a machine learning technique to identify essential genes using the experimental data of genome-wide knock-out screens from one bacterial organism to infer essential genes of another related bacterial organism. We used a broad variety of topological features, sequence characteristics and co-expression properties potentially associated with essentiality, such as flux deviations, centrality, codon frequencies of the sequences, co-regulation and phyletic retention. An organism-wise cross-validation on bacterial species yielded reliable results with good accuracies (area under the receiver-operator-curve of 75% - 81%). Finally, it was applied to drug target predictions for Salmonella typhimurium. We compared our predictions to the viability of experimental knock-outs of S. typhimurium and identified 35 enzymes, which are highly relevant to be considered as potential drug targets. Specifically, we detected promising drug targets in the non-mevalonate pathway.

Conclusions

Using elaborated features characterizing network topology, sequence information and microarray data enables to predict essential genes from a bacterial reference organism to a related query organism without any knowledge about the essentiality of genes of the query organism. In general, such a method is beneficial for inferring drug targets when experimental data about genome-wide knockout screens is not available for the investigated organism.