Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.Sequencing and comparative analysis of multiple diverse genomes is revolutionizing contemporary biology by providing a framework for interpreting and predicting the physiologic properties of an organism. A variety of emerging postgenomic techniques such as genome-wide expression profiling and monitoring of macromolecular complex formation can reveal the detailed molecular compositions of cells. New computational approaches to exploring the inherent organization of cellular networks, the mode and dynamics of interactions among cellular constituents, are in early stages of development (14,22,23). These techniques allow us to begin unraveling a major paradigm of cellular biology: how biological properties arise from the large number of components making up an individual cell.

Recent advances enabled by the Hi-C technique have unraveled many principles of chromosomal folding that were subsequently linked to disease and gene regulation. In particular, Hi-C revealed that chromosomes of animals are organized into topologically associating domains (TADs), evolutionary conserved compact chromatin domains that influence gene expression. Mechanisms that underlie partitioning of the genome into TADs remain poorly understood. To explore principles of TAD folding in Drosophila melanogaster, we performed Hi-C and poly(A) + RNA-seq in four cell lines of various origins (S2, Kc167, DmBG3-c2, and OSC). Contrary to previous studies, we find that regions between TADs (i.e., the inter-TADs and TAD boundaries) in Drosophila are only weakly enriched with the insulator protein dCTCF, while another insulator protein Su(Hw) is preferentially present within TADs. However, Drosophila inter-TADs harbor active chromatin and constitutively transcribed (housekeeping) genes. Accordingly, we find that binding of insulator proteins dCTCF and Su(Hw) predicts TAD boundaries much worse than active chromatin marks do. Interestingly, inter-TADs correspond to decompacted inter-bands of polytene chromosomes, whereas TADs mostly correspond to densely packed bands. Collectively, our results suggest that TADs are condensed chromatin domains depleted in active chromatin marks, separated by regions of active chromatin. We propose the mechanism of TAD self-assembly based on the ability of nucleosomes from inactive chromatin to aggregate, and lack of this ability in acetylated nucleosomal arrays. Finally, we test this hypothesis by polymer simulations and find that TAD partitioning may be explained by different modes of inter-nucleosomal interactions for active and inactive chromatin.

Using comparative analysis of genes, operons, and regulatory elements, we describe the cobalamin (vitamin B 12 ) biosynthetic pathway in available prokaryotic genomes. Here we found a highly conserved RNA secondary structure, the regulatory B12 element, which is widely distributed in the upstream regions of cobalamin biosynthetic/transport genes in eubacteria. In addition, the binding signal (CBL-box) for a hypothetical B 12 regulator was identified in some archaea. A search for B12 elements and CBL-boxes and positional analysis identified a large number of new candidate B 12 -regulated genes in various prokaryotes. Among newly assigned functions associated with the cobalamin biosynthesis, there are several new types of cobalt transporters, ChlI and ChlD subunits of the CobN-dependent cobaltochelatase complex, cobalt reductase BluB, adenosyltransferase PduO, several new proteins linked to the lower ligand assembly pathway, L-threonine kinase PduX, and a large number of other hypothetical proteins. Most missing genes detected within the cobalamin biosynthetic pathways of various bacteria were identified as nonorthologous substitutes. The variable parts of the cobalamin metabolism appear to be the cobalt transport and insertion, the CobG/CbiG-and CobF/CbiD-catalyzed reactions, and the lower ligand synthesis pathway. The most interesting result of analysis of B12 elements is that B 12 -independent isozymes of the methionine synthase and ribonucleotide reductase are regulated by B12 elements in bacteria that have both B 12 -dependent and B 12 -independent isozymes. Moreover, B 12 regulons of various bacteria are thought to include enzymes from known B 12 -dependent or alternative pathways. Cobalamin (CBL),1 along with chlorophyll, heme, siroheme, and coenzyme F 430 , constitute a class of the most structurally complex cofactors synthesized by bacteria. The distinctive feature of these cofactors is their tetrapyrrole-derived framework with a centrally chelated metal ion (cobalt, magnesium, iron, or nickel). Methylcobalamin and Ado-CBL, two derivatives of vitamin B 12 (cyanocobalamin) with different upper axial ligands, are essential cofactors for several important enzymes that catalyze a variety of transmethylation and rearrangement reactions. Among the most prominent vitamin B 12 -dependent enzymes in bacteria and archaea are the methionine synthase isozyme MetH from enteric bacteria; the ribonucleotide reductase isozyme NrdJ from deeply rooted eubacteria and archaea; diol dehydratases and ethanolamine ammonia lyase from enteric bacteria involved in anaerobic glycerol, 1,2-propanediol, and ethanolamine fermentation; glutamate and methylmalonyl-CoA mutases from clostridia and streptomycetes; and various CBL-dependent methyltransferases from methane-producing archaea (1-5).Most prokaryotic organisms as well as animals (including humans) and protists have enzymes that require CBL as cofactor, whereas plants and fungi are thought not to use it. Among the CBL-utilizing organisms, only some bacterial and archaeal species ...

The specific and tightly controlled transport of numerous nutrients and metabolites across cellular membranes is crucial to all forms of life. However, many of the transporter proteins involved have yet to be identified, including the vitamin transporters in various human pathogens, whose growth depends strictly on vitamin uptake. Comparative analysis of the ever-growing collection of microbial genomes coupled with experimental validation enables the discovery of such transporters. Here, we used this approach to discover an abundant class of vitamin transporters in prokaryotes with an unprecedented architecture. These transporters have energy-coupling modules comprised of a conserved transmembrane protein and two nucleotide binding proteins similar to those of ATP binding cassette (ABC) transporters, but unlike ABC transporters, they use small integral membrane proteins to capture specific substrates. We identified 21 families of these substrate capture proteins, each with a different specificity predicted by genome context analyses. Roughly half of the substrate capture proteins (335 cases) have a dedicated energizing module, but in 459 cases distributed among almost 100 gram-positive bacteria, including numerous human pathogens, different and unrelated substrate capture proteins share the same energy-coupling module. The shared use of energy-coupling modules was experimentally confirmed for folate, thiamine, and riboflavin transporters. We propose the name energycoupling factor transporters for the new class of membrane transporters.Transport proteins residing in the cytoplasmic membrane allow the selective uptake and efflux of solutes and are essential for cellular growth and metabolism (20). Reflecting the importance of transporters, between 3% and 16% of the genes in prokaryote genomes are predicted to encode transporter proteins (26). These transporters form numerous families that are diverse in structure, energy-coupling mechanisms, and substrate specificities (25). As only a small fraction of predicted transporter proteins have known substrates, the functional prediction and annotation of the specificities of transporter proteins in the rapidly growing number of sequenced genomes represent a substantial challenge (25, 36). For example, the uptake of many cofactors and their precursors is essential for the growth of various pathogenic bacteria whose genomes are sequenced, but the transport proteins involved have not yet been identified. The use of computational comparative genomic techniques including gene colocalization, cooccurrence, and coregulation analyses combined with experimental assays is a powerful approach to identify novel transporters and to uncover their cellular role (for a recent review, see reference 11).The starting point for the present analysis was our recent discovery of multicomponent transport systems for the vitamin biotin (BioYNM) and the transition metals nickel (NikMNQO) and cobalt (CbiMNQO) (14,30). These transporters all have substrate-specific components (S components), which...

Vitamin B 1 in its active form thiamin pyrophosphate is an essential coenzyme that is synthesized by coupling of pyrimidine (hydroxymethylpyrimidine; HMP) and thiazole (hydroxyethylthiazole) moieties in bacteria. Using comparative analysis of genes, operons, and regulatory elements, we describe the thiamin biosynthetic pathway in available bacterial genomes. The previously detected thiamin-regulatory element, thi box (Miranda-Rios, J., Navarro, M., and Soberon, M. (2001) Proc. Natl. Acad. Sci. U. S. A. 98, 9736 -9741), was extended, resulting in a new, highly conserved RNA secondary structure, the THI element, which is widely distributed in eubacteria and also occurs in some archaea. Search for THI elements and analysis of operon structures identified a large number of new candidate thiamin-regulated genes, mostly transporters, in various prokaryotic organisms. In particular, we assign the thiamin transporter function to yuaJ in the Bacillus/Clostridium group and the HMP transporter function to an ABC transporter thiXYZ in some proteobacteria and firmicutes. By analogy to the model of regulation of the riboflavin biosynthesis, we suggest thiamin-mediated regulation based on formation of alternative RNA structures involving the THI element. Either transcriptional or translational attenuation mechanism may operate in different taxonomic groups, dependent on the existence of putative hairpins that either act as transcriptional terminators or sequester translation initiation sites. Based on analysis of co-occurrence of the thiamin biosynthetic genes in complete genomes, we predict that eubacteria, archaea, and eukaryota have different pathways for the HMP and hydroxyethylthiazole biosynthesis.

Zinc is an important component of many proteins, but in large concentrations it is poisonous to the cell. Thus its transport is regulated by zinc repressors ZUR of proteobacteria and Gram-positive bacteria from the Bacillus group and AdcR of bacteria from the Streptococcus group. Comparative computational analysis allowed us to identify binding signals of ZUR repressors GAAATGTTATANTATAA-CATTTC for ␥-proteobacteria, GTAATGTAATAACATTAC for the Agrobacterium group, GATATGTTATAACATATC for the Rhododoccus group, TAAATCGTAATNATTACGATTTA for Gram-positive bacteria, and TTAACYRGTTAA of the streptococcal AdcR repressor. In addition to known transporters and their paralogs, zinc regulons were predicted to contain a candidate component of the ATP binding cassette, zinT (b1995 in Escherichia coli and yrpE in Bacillus subtilis). Candidate AdcR-binding sites were identified upstream of genes encoding pneumococcal histidine triad (PHT) proteins from a number of pathogenic streptococci. Protein functional analysis of this family suggests that PHT proteins are involved in the invasion process. Finally, repression by zinc was predicted for genes encoding a variety of paralogs of ribosomal proteins. The original copies of all these proteins contain zinc-ribbon motifs and thus likely bind zinc, whereas these motifs are destroyed in zinc-regulated paralogs. We suggest that the induction of these paralogs in conditions of zinc starvation leads to their incorporation in a fraction of ribosomes instead of the original ribosomal proteins; the latter are then degraded with subsequent release of some zinc for the utilization by other proteins. Thus we predict a mechanism for maintaining zinc availability for essential enzymes.

Bacterial response to nitric oxide (NO) is of major importance since NO is an obligatory intermediate of the nitrogen cycle. Transcriptional regulation of the dissimilatory nitric oxides metabolism in bacteria is diverse and involves FNR-like transcription factors HcpR, DNR, and NnrR; two-component systems NarXL and NarQP; NO-responsive activator NorR; and nitrite-sensitive repressor NsrR. Using comparative genomics approaches, we predict DNA-binding motifs for these transcriptional factors and describe corresponding regulons in available bacterial genomes. Within the FNR family of regulators, we observed a correlation of two specificity-determining amino acids and contacting bases in corresponding DNA recognition motif. Highly conserved regulon HcpR for the hybrid cluster protein and some other redox enzymes is present in diverse anaerobic bacteria, including Clostridia, Thermotogales, and delta-proteobacteria. NnrR and DNR control denitrification in alpha- and beta-proteobacteria, respectively. Sigma-54-dependent NorR regulon found in some gamma- and beta-proteobacteria contains various enzymes involved in the NO detoxification. Repressor NsrR, which was previously known to control only nitrite reductase operon in Nitrosomonas spp., appears to be the master regulator of the nitric oxides' metabolism, not only in most gamma- and beta-proteobacteria (including well-studied species such as Escherichia coli), but also in Gram-positive Bacillus and Streptomyces species. Positional analysis and comparison of regulatory regions of NO detoxification genes allows us to propose the candidate NsrR-binding motif. The most conserved member of the predicted NsrR regulon is the NO-detoxifying flavohemoglobin Hmp. In enterobacteria, the regulon also includes two nitrite-responsive loci, nipAB (hcp-hcr) and nipC (dnrN), thus confirming the identity of the effector, i.e. nitrite. The proposed NsrR regulons in Neisseria and some other species are extended to include denitrification genes. As the result, we demonstrate considerable interconnection between various nitrogen-oxides-responsive regulatory systems for the denitrification and NO detoxification genes and evolutionary plasticity of this transcriptional network.

Clostridium difficile, a Gram positive, anaerobic, spore-forming bacterium is an emergent pathogen and the most common cause of nosocomial diarrhea. Although transmission of C. difficile is mediated by contamination of the gut by spores, the regulatory cascade controlling spore formation remains poorly characterized. During Bacillus subtilis sporulation, a cascade of four sigma factors, σF and σG in the forespore and σE and σK in the mother cell governs compartment-specific gene expression. In this work, we combined genome wide transcriptional analyses and promoter mapping to define the C. difficile σF, σE, σG and σK regulons. We identified about 225 genes under the control of these sigma factors: 25 in the σF regulon, 97 σE-dependent genes, 50 σG-governed genes and 56 genes under σK control. A significant fraction of genes in each regulon is of unknown function but new candidates for spore coat proteins could be proposed as being synthesized under σE or σK control and detected in a previously published spore proteome. SpoIIID of C. difficile also plays a pivotal role in the mother cell line of expression repressing the transcription of many members of the σE regulon and activating sigK expression. Global analysis of developmental gene expression under the control of these sigma factors revealed deviations from the B. subtilis model regarding the communication between mother cell and forespore in C. difficile. We showed that the expression of the σE regulon in the mother cell was not strictly under the control of σF despite the fact that the forespore product SpoIIR was required for the processing of pro-σE. In addition, the σK regulon was not controlled by σG in C. difficile in agreement with the lack of pro-σK processing. This work is one key step to obtain new insights about the diversity and evolution of the sporulation process among Firmicutes.