Interim report on genomics of Escherichia coli

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

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.

Genome variations associated with viral susceptibility and calcification in Emiliania huxleyi

Emiliania huxleyi, a key player in the global carbon cycle is one of the best studied coccolithophores with respect to biogeochemical cycles, climatology, and host-virus interactions. Strains of E. huxleyi show phenotypic plasticity regarding growth behaviour, light-response, calcification, acidification, and virus susceptibility. This phenomenon is likely a consequence of genomic differences, or transcriptomic responses, to environmental conditions or threats such as viral infections. We used an E. huxleyi genome microarray based on the sequenced strain CCMP1516 (reference strain) to perform comparative genomic hybridizations (CGH) of 16 E. huxleyi strains of different geographic origin. We investigated the genomic diversity and plasticity and focused on the identification of genes related to virus susceptibility and coccolith production (calcification). Among the tested 31940 gene models a core genome of 14628 genes was identified by hybridization among 16 E. huxleyi strains. 224 probes were characterized as specific for the reference strain CCMP1516. Compared to the sequenced E. huxleyi strain CCMP1516 variation in gene content of up to 30 percent among strains was observed. Comparison of core and non-core transcripts sets in terms of annotated functions reveals a broad, almost equal functional coverage over all KOG-categories of both transcript sets within the whole annotated genome. Within the variable (non-core) genome we identified genes associated with virus susceptibility and calcification. Genes associated with virus susceptibility include a Bax inhibitor-1 protein, three LRR receptor-like protein kinases, and mitogen-activated protein kinase. Our list of transcripts associated with coccolith production will stimulate further research, e.g. by genetic manipulation. In particular, the V-type proton ATPase 16 kDa proteolipid subunit is proposed to be a plausible target gene for further calcification studies.

Two rep genes in small cryptic plasmid pKST21 of Escherichia coli

The complete nucleotide sequence of a small cryptic plasmid pKST21 from Escherichia coli was determined. This plasmid is 1,460 bp long with an overall GC content of 51 %. Based on sequence analysis, the presence of two segments with different average GC density was observed. The segment with higher GC content revealed 98-90 % similarity to several small plasmids of E. coli and to pCR1 from Gram-positive Corynebacterium renale. Plasmid pKST21 possesses two conversely oriented open reading frames encoding proteins with a high degree of amino acid identity to Rep proteins involved in replication. ORF1 encodes replication protein similar to RepA protein of Bartonella tribocorum or Bacillus cereus plasmids or to the putative plasmid Rep protein from ecologically close Selenomonas ruminantium. ORF2 similarly encodes a replication protein, which shares 97 % homology with Rep protein from C. renale. Genetic diversity observed in plasmid pKST21 indicates a mosaic structure of the plasmid with different segments acquired from different sources. Deletion analysis showed that both fragments carrying the repA and repB genes are necessary for the replication of pKST21 in E. coli. The presence of plasmid with the same gene composition was revealed in 14 % of tested E. coli isolates from the rumen of sheep. All these strains produced identical ERIC-PCR profiles indicating isogenic origin of the strain and lack of horizontal gene transfer of pKST21 plasmid.

Functional genomics via metabolic footprinting: monitoring metabolite secretion by Escherichia coli tryptophan metabolism mutants using FT-IR and direct injection electrospray mass spectrometry

We sought to test the hypothesis that mutant bacterial strains could be discriminated from each other on the basis of the metabolites they secrete into the medium (their ‘metabolic footprint’), using two methods of ‘global’ metabolite analysis (FT–IR and direct injection electrospray mass spectrometry). The biological system used was based on a published study of Escherichia coli tryptophan mutants that had been analysed and discriminated by Yanofsky and colleagues using transcriptome analysis. Wild-type strains supplemented with tryptophan or analogues could be discriminated from controls using FT–IR of 24 h broths, as could each of the mutant strains in both minimal and supplemented media. Direct injection electrospray mass spectrometry with unit mass resolution could also be used to discriminate the strains from each other, and had the advantage that the discrimination required the use of just two or three masses in each case. These were determined via a genetic algorithm. Both methods are rapid, reagentless, reproducible and cheap, and might beneficially be extended to the analysis of gene knockout libraries.

Bacterial structural genomics initiative: overview of methods and technologies applied to the process of structure determination

The focus over the last several years on increasing the number of three-dimensional structures of macromolecules by implementation of high throughput methodology has led to the establishment of dedicated structural genomics programs around the world. These worldwide efforts have in turn led to development of novel, parallelized approaches to cloning, expression, purification, and crystallization of proteins. This chapter describes in some detail the approaches and protocols that have been implemented in the Bacterial Structural Genomics Initiative.

Introduction and Perspectives

Classically, metabolism was investigated by studying molecular characteristics of enzymes and their regulators in isolation. This reductionistic approach successfully established mechanistic relationships with the immediate interacting neighbors and allowed reconstruction of network structures. Severely underdeveloped was the ability to make precise predictions about the integrated operation of pathways and networks that emerged from the typically nonlinear and complex interactions of proteins and metabolites. The burden of metabolic engineering is a consequence of this fact-one cannot yet predict with any certainty precisely what needs to be engineered to produce more complex phenotypes. What was and still is missing are concepts, methods, and algorithms to integrate data and information into a quantitatively coherent whole, as well as theoretical concepts to reliably predict the consequence of environmental stimuli or genetic interventions. This introduction and perspective to Domain 3, Metabolism and Metabolic Fluxes, starts with a brief overview of the panoply of global measurement technologies that herald the dawning of systems biology and whose impact on metabolic research is apparent throughout the Domain 3. In the middle section, applications to Escherichia coli are used to illustrate general concepts and successes of computational methods that approach metabolism as a network of interacting elements, and thus have potential to fill the gap in quantitative data and information integration. The final section highlights prospective focus areas for future metabolic research, including functional genomics, eludication of evolutionary principles, and the integration of metabolism with regulatory networks.

A novel data mining method to identify assay-specific signatures in functional genomic studies

Background:

The highly dimensional data produced by functional genomic (FG) studies makes it difficult to visualize relationships between gene products and experimental conditions (i.e., assays). Although dimensionality reduction methods such as principal component analysis (PCA) have been very useful, their application to identify assay-specific signatures has been limited by the lack of appropriate methodologies. This article proposes a new and powerful PCA-based method for the identification of assay-specific gene signatures in FG studies.

Results:

The proposed method (PM) is unique for several reasons. First, it is the only one, to our knowledge, that uses gene contribution, a product of the loading and expression level, to obtain assay signatures. The PM develops and exploits two types of assay-specific contribution plots, which are new to the application of PCA in the FG area. The first type plots the assay-specific gene contribution against the given order of the genes and reveals variations in distribution between assay-specific gene signatures as well as outliers within assay groups indicating the degree of importance of the most dominant genes. The second type plots the contribution of each gene in ascending or descending order against a constantly increasing index. This type of plots reveals assay-specific gene signatures defined by the inflection points in the curve. In addition, sharp regions within the signature define the genes that contribute the most to the signature. We proposed and used the curvature as an appropriate metric to characterize these sharp regions, thus identifying the subset of genes contributing the most to the signature. Finally, the PM uses the full dataset to determine the final gene signature, thus eliminating the chance of gene exclusion by poor screening in earlier steps. The strengths of the PM are demonstrated using a simulation study, and two studies of real DNA microarray data – a study of classification of human tissue samples and a study of E. coli cultures with different medium formulations.

Conclusion

We have developed a PCA-based method that effectively identifies assay-specific signatures in ranked groups of genes from the full data set in a more efficient and simplistic procedure than current approaches. Although this work demonstrates the ability of the PM to identify assay-specific signatures in DNA microarray experiments, this approach could be useful in areas such as proteomics and metabolomics.

New multiple-deletion method for the Corynebacterium glutamicum genome, using a mutant lox sequence

Due to the difficulty of multiple deletions using the Cre/loxP system, a simple, markerless multiple-deletion method based on a Cre/mutant lox system combining a right-element (RE) mutant lox site with a left-element (LE) mutant lox site was employed for large-scale genome rearrangements in Corynebacterium glutamicum. Eight distinct genomic regions that had been identified previously by comparative analysis of C. glutamicum R and C. glutamicum 13032 genomes were targeted for deletion. By homologous recombination, LE and RE mutant lox sites were integrated at each end of a target region. Highly efficient and accurate deletions between the two chromosomal mutant lox sites in the presence of Cre recombinase were realized. A deletion mutant lacking 190 kb of chromosomal regions, encoding a total of 188 open reading frames (ORFs), was obtained. These deletions represent the largest genomic excisions in C. glutamicum reported to date. Despite the loss of numerous predicted ORFs, the mutant exhibited normal growth under standard laboratory conditions. The Cre/loxP system using a pair of mutant lox sites provides a new, efficient genome rearrangement technique for C. glutamicum. It should facilitate the understanding of genome functions of microorganisms.

Multiple large segment deletion method for Corynebacterium glutamicum

A precise and scarless genome excision method, employing the Cre/loxP system in concert with double-strand break (DSB)-stimulated intramolecular recombination was developed. The DSBs were mediated by the restriction endonuclease, I-SceI. It permitted multiple deletions of independent 14-, 43-, and 10-kb-long genomic regions on the Corynebacterium glutamicum genome. Accuracy of deletion was confirmed by the loss of marker genes, PCR, and sequencing of new genome joints. Eleven, 58, and 4 genes were predicted on the 14-, 43-, and 10-kb deleted regions, respectively. Although the resultant mutant lost a total of 67 kb encoding 73 genes, it still exhibited normal growth under standard laboratory conditions. Such a large segment deletion method in which multiple, successive deletions are possible is useful for genome engineering.

Large-scale engineering of the Corynebacterium glutamicum genome

The engineering of Corynebacterium glutamicum is important for enhanced production of biochemicals. To construct an improved C. glutamicum genome, we developed a precise genome excision method based on the Cre/loxP recombination system and successfully deleted 11 distinct genomic regions identified by comparative analysis of C. glutamicum genomes. Despite the loss of several predicted open reading frames, the mutant cells exhibited normal growth under standard laboratory conditions. With a total of 250 kb (7.5% of the genome), the 11 genomic regions were loaded with cryptic prophages, transposons, and genes of unknown function which were dispensable for cell growth, indicating recent horizontal acquisitions to the genome. This provides an interesting background for functional genomic studies and can be used in the improvement of cell traits.

Ribosomal RNA probes and microarrays: their potential use in assessing microbial biodiversity

The awareness that global biological diversity is affected by numerous, mostly human-made threats has made biodiversity assessment an important scientific issue for decades. Biodiversity includes different levels of complexity, such as community diversity, habitat diversity, genetic diversity, and species diversity. The application of molecular methods to answer ecological questions permits issues of biodiversity to be addressed at all levels. Microorganisms dominate global biological diversity in terms of their species numbers. However, their small size and limited morphological features make it challenging to obtain a comprehensive view of their biodiversity. The application of ribosomal RNA (rRNA) probes contributes significantly to the assessment of biodiversity at the molecular level. DNA microarrays offer a great potential to facilitate the application of molecular probes and other DNA analytical methods to answer ecological and biodiversity questions. We provide an introduction into the application of rRNA probes and DNA microarrays for the assessment of microbial biodiversity, as well as protocols for the implementation of DNA microarrays.

EchoBASE: an integrated post-genomic database for Escherichia coli

EchoBASE (http://www.ecoli-york.org) is a relational database designed to contain and manipulate information from post-genomic experiments using the model bacterium Escherichia coli K-12. Its aim is to collate information from a wide range of sources to provide clues to the functions of the approximately 1500 gene products that have no confirmed cellular function. The database is built on an enhanced annotation of the updated genome sequence of strain MG1655 and the association of experimental data with the E.coli genes and their products. Experiments that can be held within EchoBASE include proteomics studies, microarray data, protein–protein interaction data, structural data and bioinformatics studies. EchoBASE also contains annotated information on ‘orphan’ enzyme activities from this microbe to aid characterization of the proteins that catalyse these elusive biochemical reactions.

Specialized microbial databases for inductive exploration of microbial genome sequences

Background

The enormous amount of genome sequence data asks for user-oriented databases to manage sequences and annotations. Queries must include search tools permitting function identification through exploration of related objects.

Methods

The GenoList package for collecting and mining microbial genome databases has been rewritten using MySQL as the database management system. Functions that were not available in MySQL, such as nested subquery, have been implemented.

Results

Inductive reasoning in the study of genomes starts from "islands of knowledge", centered around genes with some known background. With this concept of "neighborhood" in mind, a modified version of the GenoList structure has been used for organizing sequence data from prokaryotic genomes of particular interest in China. GenoChore , a set of 17 specialized end-user-oriented microbial databases (including one instance of Microsporidia, Encephalitozoon cuniculi, a member of Eukarya) has been made publicly available. These databases allow the user to browse genome sequence and annotation data using standard queries. In addition they provide a weekly update of searches against the world-wide protein sequences data libraries, allowing one to monitor annotation updates on genes of interest. Finally, they allow users to search for patterns in DNA or protein sequences, taking into account a clustering of genes into formal operons, as well as providing extra facilities to query sequences using predefined sequence patterns.

Conclusion

This growing set of specialized microbial databases organize data created by the first Chinese bacterial genome programs (ThermaList, Thermoanaerobacter tencongensis, LeptoList, with two different genomes of Leptospira interrogans and SepiList, Staphylococcus epidermidis) associated to related organisms for comparison.

Escherichia coli and Salmonella 2000: the View From Here

In 1995, an editorial in Science (267:1575) commented that predictions made some 25 years previously regarding "Biology and the Future of Man" were largely fulfilled but that "the most revolutionary and unexpected findings were not predicted." We would be glad to do as well! As we stated at the beginning, our work as editors of the Escherichia coli and Salmonella book did not endow us with special powers of prophecy but it does permit us to express our excitement for the future. In our opinion, E. coli and S. enterica will continue to play a central role in biological research. This is not because they are intrinsically more interesting than any other bacteria, as we believe that all bacteria are equally interesting. However, knowledge builds on knowledge, and it is here that these two species continue to have a large edge not only over other microorganisms but also, for some time to come, over all other forms of life. It is interesting in this connection that biotechnology, having made detours through other microorganisms, always seems to return to E. coli.

Deciphering the function of an ORF: Salmonella enterica DeoM protein is a new mutarotase specific for deoxyribose

We identified in Salmonella enterica serovar Typhi a cluster of four genes encoding a deoxyribokinase (DeoK), a putative permease (DeoP), a repressor (DeoQ), and an open reading frame encoding a 337 amino acid residues protein of unknown function. We show that the latter protein, called DeoM, is a hexamer whose synthesis is increased by a factor over 5 after induction with deoxyribose. The CD spectrum of the purified recombinant protein indicated a dominant contribution of betatype secondary structure and a small content of alpha-helix. Temperature and guanidinium hydrochloride induced denaturation of DeoM indicated that the hexamer dissociation and monomer unfolding are coupled processes. DeoM exhibits 12.5% and 15% sequence identity with galactose mutarotase from Lactococcus lactis and respectively Escherichia coli, which suggested that these three proteins share similar functions. Polarimetric experiments demonstrated that DeoM is a mutarotase with high specificity for deoxyribose. Site-directed mutagenesis of His183 in DeoM, corresponding to a catalytically active residue in GalM, yielded an almost inactive deoxyribose mutarotase. DeoM was crystallized and diffraction data collected for two crystal systems, confirmed its hexameric state. The possible role of the protein and of the entire gene cluster is discussed in connection with the energy metabolism of S. enterica under particular growth conditions.

Model systems in drug discovery: chemical genetics meets genomics

Animal model systems are an intricate part of the discovery and development of new medicines. The sequencing of not only the human genome but also those of the various pathogenic bacteria, the nematode Caenorhabditis elegans, the fruitfly Drosophila, and the mouse has enabled the discovery of new drug targets to push forward at an unprecedented pace. The knowledge and tools in these "model" systems are allowing researchers to carry out experiments more efficiently and are uncovering previously hidden biological connections. While the history of bacteria, yeast, and mice in drug discovery are long, their roles are ever evolving. In contrast, the history of Drosophila and C. elegans at pharmaceutical companies is short. We will briefly review the historic role of each model organism in drug discovery and then update the readers as to the abilities and liabilities of each model within the context of drug development.

Evolution of the protein repertoire

Most proteins have been formed by gene duplication, recombination, and divergence. Proteins of known structure can be matched to about 50% of genome sequences, and these data provide a quantitative description and can suggest hypotheses about the origins of these processes.

Regulation of the isofunctional genes ubiD and ubiX of the ubiquinone biosynthetic pathway of Escherichia coli

The expressions of the isofunctional genes ubiD and ubiX of the ubiquinone biosynthetic pathway of Escherichia coli were compared under a variety of growth conditions and in several genetic backgrounds. LacZ operon fusions were constructed and were inserted in single copies into strain MC4100 and into its fnr, arcA or hemA carrying derivatives. During aerobic growth the expressions of both ubiD and ubiX depended on the carbon source: succinate>glycerol>glucose. Mutations in fnr, arcA or hemA increased the expressions of both genes. During anaerobic growth in LB medium glucose strongly inhibited the expression of ubiD but not of ubiX.

A microarray-based antibiotic screen identifies a regulatory role for supercoiling in the osmotic stress response of Escherichia coli

Changes in DNA supercoiling are induced by a wide range of environmental stresses in Escherichia coli, but the physiological significance of these responses remains unclear. We now demonstrate that an increase in negative supercoiling is necessary for transcriptional activation of a large subset of osmotic stress-response genes. Using a microarray-based approach, we have characterized supercoiling-dependent gene transcription by expression profiling under conditions of high salt, in conjunction with the microbial antibiotics novobiocin, pefloxacin, and chloramphenicol. Algorithmic clustering and statistical measures for gauging cellular function show that this subset is enriched for genes critical in osmoprotectant transport/synthesis and rpoS-driven stationary phase adaptation. Transcription factor binding site analysis also supports regulation by the global stress sigma factor rpoS. In addition, these studies implicate 60 uncharacterized genes in the osmotic stress regulon, and offer evidence for a broader role for supercoiling in the control of stress-induced transcription.

Induction of SOS response, cellular efflux and oxidative stress response genes by chlorambucil in DNA repair-deficient Escherichia coli cells (ada, ogt and mutS)

Chlorambucil (CLB) is a bifunctional alkylating drug widely used as an anticancer agent and as an immunosuppressant. It is known to be mutagenic, teratogenic and carcinogenic. The cellular actions of CLB have remained poorly investigated. It is very likely that DNA damage and its repair are the key elements determining the destiny of CLB-exposed cells. We investigated the role of two specific DNA repair pathways involved in CLB-induced mutagenicity and gene expression changes by using Escherichia coli strains lacking either (i) two DNA methyltransferase functions (O(6)-methylguanine-DNA methyltransferase I (ada) and II (ogt)), or (ii) mismatch repair (MutS (mutS)). Mutagenicity was determined as the development of ciproxin and rifampicin resistance and the gene expression changes were assessed using expression profiling of all E. coli 4290 open reading frames (ORFs) by cDNA array. Chlorambucil-induced mutants in mutS cells, implying the importance of mismatch repair in preventing CLB-induced mutations. It also induced mutants in the ada, ogt strain, but to a lesser extent than in the wild-type strain. The simultaneous upregulation of several genes of the SOS response, cellular efflux and oxidative stress response, was demonstrated in both of the DNA repair-deficient strains but not in the wild-type cells. These and our previous results show that single-gene knock-out cells use specific gene regulation strategies to avoid mutations and cell death induced by agents such as chlorambucil.

ASAP, a systematic annotation package for community analysis of genomes

ASAP (a systematic annotation package for community analysis of genomes) is a relational database and web interface developed to store, update and distribute genome sequence data and functional characterization (https://asap.ahabs.wisc.edu/annotation/php/ASAP1.htm). ASAP facilitates ongoing community annotation of genomes and tracking of information as genome projects move from preliminary data collection through post-sequencing functional analysis. The ASAP database includes multiple genome sequences at various stages of analysis, corresponding experimental data and access to collections of related genome resources. ASAP supports three levels of users: public viewers, annotators and curators. Public viewers can currently browse updated annotation information for Escherichia coli K-12 strain MG1655, genome-wide transcript profiles from more than 50 microarray experiments and an extensive collection of mutant strains and associated phenotypic data. Annotators worldwide are currently using ASAP to participate in a community annotation project for the Erwinia chrysanthemi strain 3937 genome. Curation of the E. chrysanthemi genome annotation as well as those of additional published enterobacterial genomes is underway and will be publicly accessible in the near future.

Structural and nucleotide-binding properties of YajQ and YnaF, two Escherichia coli proteins of unknown function

Structural genomics is a new approach in functional assignment of proteins identified via whole-genome sequencing programs. Its rationale is that nonhomologous proteins performing similar or related biological functions might have similar tertiary structure. We used dye pseudoaffinity chromatography, two-dimensional gel electrophoresis, and mass spectrometry to identify two novel Escherichia coli nucleotide-binding proteins, YnaF and YajQ. YnaF exhibited significant sequence identity with MJ0577, an ATP-binding protein from a hyperthermophile (Methanococcus jannaschii), and with UspA, a protein from Haemophilus influenzae that belongs to the Universal Stress Protein family. YnaF conserves the ATP-binding site and the dimeric structure observed in the crystal of MJ0577. The protein YajQ, present in many bacterial genomes, is missing in eukaryotes. In the absence of significant similarities of YajQ to any solved structure, we determined its structural and ligand-binding properties by NMR and isothermal titration calorimetry. We demonstrate that YajQ is composed of two domains, each centered on a beta-sheet, that are connected by two helical segments. NMR studies, corroborated with local sequence conservation among YajQ homologs in various bacteria, indicate that one of the beta-sheets is mostly involved in biological activity.

Specificity of chloroplast-localized peptide deformylases as determined with peptide analogs of chloroplast-translated proteins

Peptide deformylase (DEF; EC 3.5.1.88) removes the N-formyl group from nascent polypeptides. Two nuclear-encoded DEFs in Arabidopsis thaliana (At) are localized to chloroplasts, and thus, the N-termini of chloroplast-translated proteins may be a consequence of AtDEFs' substrate specificity. Using peptide analogs of select chloroplast-translated proteins, AtDEF1 activity was as much as 100-fold lower than AtDEF2 activity and showed little variance with peptide sequence. However, AtDEF2 activity was significantly influenced by peptide sequence, with the most efficiently processed substrate mimicking the N-terminus of the nascent D1 polypeptide, a core protein of photosystem II. Though AtDEF2's specificity was predictive of N-formyl retention for some chloroplast proteins, exceptions suggests that additional factors in vivo aid in determining the retention of an N-formyl group.

Physiological genomics of Escherichia coli protein families

The well-researched Escherichia coli genome offers the opportunity to explore the value of using protein families within a single organism to enrich functional annotation procedures and to study mechanisms of protein evolution. Having identified multimodular proteins resulting from gene fusion, and treated each module as a separate protein, nonoverlapping sequence-similar families in E. coli could be assembled. Of 3,902 proteins of length 100 residues or more, 2,415 clustered into 609 protein families. The relatedness of function among members of each family was dissected in detail. Data on paralogous protein families provides valuable information in attributing putative function to unknown genes, supplementing existing function annotation. Enzymes, transporters, and regulators represent the three major types of proteins in E. coli. They are shown to have distinctive patterns in gene duplication and divergence and gene fusion, suggesting that details of protein evolution have been different for genes in these categories. Data for the complete list of paralogous protein families and updated functional annotation for E. coli K-12 are accessible in GenProtEC (http://genprotec.mbl.edu).

Engineering a reduced Escherichia coli genome

Our goal is to construct an improved Escherichia coli to serve both as a better model organism and as a more useful technological tool for genome science. We developed techniques for precise genomic surgery and applied them to deleting the largest K-islands of E. coli, identified by comparative genomics as recent horizontal acquisitions to the genome. They are loaded with cryptic prophages, transposons, damaged genes, and genes of unknown function. Our method leaves no scars or markers behind and can be applied sequentially. Twelve K-islands were successfully deleted, resulting in an 8.1% reduced genome size, a 9.3% reduction of gene count, and elimination of 24 of the 44 transposable elements of E. coli. These are particularly detrimental because they can mutagenize the genome or transpose into clones being propagated for sequencing, as happened in 18 places of the draft human genome sequence. We found no change in the growth rate on minimal medium, confirming the nonessential nature of these islands. This demonstration of feasibility opens the way for constructing a maximally reduced strain, which will provide a clean background for functional genomics studies, a more efficient background for use in biotechnology applications, and a unique tool for studies of genome stability and evolution.

Chlorambucil-induced high mutation rate and suicidal gene downregulation in a base excision repair-deficient Escherichia coli strain

Chlorambucil (CLB; N,N-bis(2-chloroethyl)-p-aminophenylbutyric acid) is a bifunctional alkylating agent widely used as an anticancer drug and also as an immunosuppressant. Its chemical structure and clinical experience indicate that CLB is mutagenic and carcinogenic. We have investigated the ability of CLB to induce mutations and gene expression changes in the wild-type (WT) Escherichia coli strain AB1157 and in the base excision repair-deficient (alkA1, tag-1) E. coli strain MV1932 using a rifampicin (rif) forward mutation system and a cDNA array method. The results showed that CLB is a potent mutagen in MV1932 cells compared with the E. coli WT strain AB1157, emphasizing the role of 3-methyladenine DNA glycosylases I and II in protecting the cells from CLB-induced DNA damage and subsequent mutations. Global gene expression profiling revealed that nine genes in WT E. coli and 100 genes in MV1932, of a total of 4290 genes, responded at least 2.5-fold to CLB. Interestingly, all of these MV1932 genes were downregulated, while 22% were upregulated in WT cells. The downregulated genes in MV1932 represented most (19/23) functional categories, and unexpectedly, many of them code for proteins responsible for genomic integrity. These include: (i) RecF (SOS-response, adaptive mutation), (ii) RecC (resistance to cross-linking agents), (iii) HepA (DNA repair, a possible substitute of RecBCD), (iv) Ssb (DNA recombination repair, controls RecBCD), and (v) SbcC (genetic recombination). Our results strongly suggest that in addition to the DNA damage itself, the downregulation of central protecting genes is responsible for the decreased cell survival (demonstrated in a previous work) and the increased mutation rate (this work) of DNA repair-deficient cells, when exposed to CLB.

Metabolic flux responses to pyruvate kinase knockout in Escherichia coli

The intracellular carbon flux distribution in wild-type and pyruvate kinase-deficient Escherichia coli was estimated using biosynthetically directed fractional 13C labeling experiments with [U-13C6]glucose in glucose- or ammonia-limited chemostats, two-dimensional nuclear magnetic resonance (NMR) spectroscopy of cellular amino acids, and a comprehensive isotopomer model. The general response to disruption of both pyruvate kinase isoenzymes in E. coli was a local flux rerouting via the combined reactions of phosphoenolpyruvate (PEP) carboxylase and malic enzyme. Responses in the pentose phosphate pathway and the tricarboxylic acid cycle were strongly dependent on the environmental conditions. In addition, high futile cycling activity via the gluconeogenic PEP carboxykinase was identified at a low dilution rate in glucose-limited chemostat culture of pyruvate kinase-deficient E. coli, with a turnover that is comparable to the specific glucose uptake rate. Furthermore, flux analysis in mutant cultures indicates that glucose uptake in E. coli is not catalyzed exclusively by the phosphotransferase system in glucose-limited cultures at a low dilution rate. Reliability of the flux estimates thus obtained was verified by statistical error analysis and by comparison to intracellular carbon flux ratios that were independently calculated from the same NMR data by metabolic flux ratio analysis.

Genomic functional annotation using co-evolution profiles of gene clusters

BACKGROUND

The current speed of sequencing already exceeds the capability of annotation, creating a potential bottleneck. A large proportion of the genes in microbial genomes remains uncharacterized. Here we propose a new method for functional annotation using the conservation patterns of gene clusters. If several gene clusters show the same coevolution pattern across different genomes it is reasonable to infer they are functionally related. The gene cluster phylogenetic profile integrates chromosomal proximity information and phylogenetic profile information and allows us to infer functional dependences between the gene clusters even at great distance on the chromosome.

RESULTS

As a proof of concept, we applied our method to the genome of Escherichia coli K12 strain. Our method establishes functional relationships among 176 gene clusters, comprising 738 E. coli genes. The accuracy of pair phylogenetic profiles was compared with the single-gene phylogenetic profile and was shown to be higher. As a result, we are able to suggest functional roles for several previously unknown genes or unknown genomic regions in E. coli. We also examined the robustness of coevolution signals across a larger set of genomes and suggest a possible upper limit of accuracy for the phylogenetic profile methods.

CONCLUSIONS

The higher-order phylogenetic profiles, such as the gene-pair phylogenetic profiles, can detect functional dependences that are missed by using conventional single-gene phylogenetic profile or the chromosomal proximity method only. We show that the gene-pair phylogenetic profile is more accurate than the single-gene phylogenetic profiles.

Microarrays in ecology and evolution: a preview

Microarray technology provides a new tool with which molecular ecologists and evolutionary biologists can survey genome-wide patterns of gene expression within and among species. New analytical approaches based on analysis of variance will allow quantification of the contributions of among individual variation, genotype, sex, microenvironment, population structure, and geography to variation in gene expression. Applications of this methodology are reviewed in relation to studies of mechanisms of adaptation and divergence; delineation of developmental and physiological pathways and networks; characterization of quantitative genetic parameters at the level of transcription ('quantitative genomics'); molecular dissection of parasitism and symbiosis; and studies of the diversification of gene content. Establishment of microarray resources is neither prohibitively expensive nor technologically demanding, and a commitment to development of gene expression profiling methods for nonmodel organisms could have a tremendous impact on molecular and genetic research at the interface of organismal and population biology.

Biodegradation of aromatic compounds by Escherichia coli

Although Escherichia coli has long been recognized as the best-understood living organism, little was known about its abilities to use aromatic compounds as sole carbon and energy sources. This review gives an extensive overview of the current knowledge of the catabolism of aromatic compounds by E. coli. After giving a general overview of the aromatic compounds that E. coli strains encounter and mineralize in the different habitats that they colonize, we provide an up-to-date status report on the genes and proteins involved in the catabolism of such compounds, namely, several aromatic acids (phenylacetic acid, 3- and 4-hydroxyphenylacetic acid, phenylpropionic acid, 3-hydroxyphenylpropionic acid, and 3-hydroxycinnamic acid) and amines (phenylethylamine, tyramine, and dopamine). Other enzymatic activities acting on aromatic compounds in E. coli are also reviewed and evaluated. The review also reflects the present impact of genomic research and how the analysis of the whole E. coli genome reveals novel aromatic catabolic functions. Moreover, evolutionary considerations derived from sequence comparisons between the aromatic catabolic clusters of E. coli and homologous clusters from an increasing number of bacteria are also discussed. The recent progress in the understanding of the fundamentals that govern the degradation of aromatic compounds in E. coli makes this bacterium a very useful model system to decipher biochemical, genetic, evolutionary, and ecological aspects of the catabolism of such compounds. In the last part of the review, we discuss strategies and concepts to metabolically engineer E. coli to suit specific needs for biodegradation and biotransformation of aromatics and we provide several examples based on selected studies. Finally, conclusions derived from this review may serve as a lead for future research and applications.

LuxArray, a high-density, genomewide transcription analysis of Escherichia coli using bioluminescent reporter strains

A sequenced collection of plasmid-borne random fusions of Escherichia coli DNA to a Photorhabdus luminescens luxCDABE reporter was used as a starting point to select a set of 689 nonredundant functional gene fusions. This group, called LuxArray 1.0, represented 27% of the predicted transcriptional units in E. coli. High-density printing of the LuxArray 1.0 reporter strains to membranes on agar plates was used for simultaneous reporter gene assays of gene expression. The cellular response to nalidixic acid perturbation was analyzed using this format. As expected, fusions to promoters of LexA-controlled SOS-responsive genes dinG, dinB, uvrA, and ydjM were found to be upregulated in the presence of nalidixic acid. In addition, six fusions to genes not previously known to be induced by nalidixic acid were also reproducibly upregulated. The responses of two of these, fusions to oraA and yigN, were induced in a LexA-dependent manner by both nalidixic acid and mitomycin C, identifying these as members of the LexA regulon. The responses of the other four were neither induced by mitomycin C nor dependent on lexA function. Thus, the promoters of ycgH, intG, rihC, and a putative operon consisting of lpxA, lpxB, rnhB, and dnaE were not generally DNA damage responsive and represent a more specific response to nalidixic acid. These results demonstrate that cellular arrays of reporter gene fusions are an important alternative to DNA arrays for genomewide transcriptional analyses.

MultiFun, a multifunctional classification scheme for Escherichia coli K-12 gene products

An enriched classification system for cellular functions of gene products of Escherichia coli K-12 was developed based on the initial classification by Riley. In the new classification scheme, MultiFun, cellular functions are divided into 10 major categories: Metabolism, Information Transfer, Regulation, Transport, Cell Processes, Cell Structure, Location, Extra-chromosomal Origin, DNA Site, and Cryptic Gene. These major categories are further sub-divided into a hierarchical scheme. Two thousand nine hundred twenty-two gene products of E. coli K-12 were assigned to one or more functions depending on the role they play in the cell. Functional assignments were made to 66% of E. coli gene products, ranging from 1 to 16 assignments per gene product. The expansion of cellular function categories and the assignment to more than one category (multifunction) provides a more complete description of the gene products and their roles and hence better reflects the functional complexity of organisms. We believe this classification system will be useful in the field of genome analysis, both for annotation purposes and for comparative studies. The functional classification scheme and the cellular function assignments made to E. coli gene products can be accessed from the web at the databases GenProtEC (http://genprotec.mbl.edu) and EcoCyc (http://www.ecocyc.org).

Divergence of function in sequence-related groups of Escherichia coli proteins

The most prominent mechanism of molecular evolution is believed to have been duplication and divergence of genes. Proteins that belong to sequence-related groups in any one organism are candidates to have emerged from such a process and to share a common ancestor. Groups of proteins in Escherichia coli having sequence similarity are mostly composed of proteins with closely related function, but some groups comprise proteins with unrelated functions. In order to understand how function can change while sequences remain similar, we have examined some of these groups in detail. The enzymes analyzed in this work include representatives of amidotransferases, phosphotransferases, decarboxylases, and others. Most sequence-related groups contain enzymes that are in the same classes of Enzyme Commission (EC) numbers. We have concentrated on groups that are heterogeneous in that respect, and also on groups containing more than one enzyme of any pathway. We find that although the EC number may differ, the reaction chemistry of these sequence-related proteins is the same or very similar. Some of these families illustrate how diversification has taken place in evolution, using common features of either reaction chemistry or ligand specificity, or both, to create catalysts for different kinds of biochemical reactions. This information has relevance to the area of functional genomics in which the activities of gene products of unknown reading frames are attributed by analogy to the functions of sequence-related proteins of known function.

Escherichia coli and Salmonella 2000: the view from here

Five years after the publication of the second edition of the reference book Escherichia coli and Salmonella: Cellular and Molecular Biology, and on the eve of launching a successor venture, the editors and colleagues examine where we stand in our quest for an understanding of these organisms. The main areas selected for this brief inquiry are genomics, evolution, molecular multifunctionality, functional backups, regulation of gene expression, cell biology, sensing of the environment, and ecology.

A functional update of the Escherichia coli K-12 genome

BACKGROUND

Since the genome of Escherichia coli K-12 was initially annotated in 1997, additional functional information based on biological characterization and functions of sequence-similar proteins has become available. On the basis of this new information, an updated version of the annotated chromosome has been generated.

RESULTS

The E. coli K-12 chromosome is currently represented by 4,401 genes encoding 116 RNAs and 4,285 proteins. The boundaries of the genes identified in the GenBank Accession U00096 were used. Some protein-coding sequences are compound and encode multimodular proteins. The coding sequences (CDSs) are represented by modules (protein elements of at least 100 amino acids with biological activity and independent evolutionary history). There are 4,616 identified modules in the 4,285 proteins. Of these, 48.9% have been characterized, 29.5% have an imputed function, 2.1% have a phenotype and 19.5% have no function assignment. Only 7% of the modules appear unique to E. coli, and this number is expected to be reduced as more genome data becomes available. The imputed functions were assigned on the basis of manual evaluation of functions predicted by BLAST and DARWIN analyses and by the MAGPIE genome annotation system.

CONCLUSIONS

Much knowledge has been gained about functions encoded by the E. coli K-12 genome since the 1997 annotation was published. The data presented here should be useful for analysis of E. coli gene products as well as gene products encoded by other genomes.