Search

Summary Xenorhabdus bovienii (SS-2004) bacteria reside in the intestine of the infective-juvenile (IJ) stage of the entomopathogenic nematode, Steinernema jollieti. The recent sequencing of the X. bovienii genome facilitates its use as a model to understand host – symbiont interactions. To provide a biological foundation for such studies, we characterized X. bovienii in vitro and host interaction phenotypes. Within the nematode host X. bovienii was contained within a membrane bound envelope that also enclosed the nematode-derived intravesicular structure. Steinernema jollieti nematodes cultivated on mixed lawns of X. bovienii expressing green or DsRed fluorescent proteins were predominantly colonized by one or the other strain, suggesting the colonizing population is founded by a few cells. Xenorhabdus bovienii exhibits phenotypic variation between orange-pigmented primary form and cream-pigmented secondary form. Each form can colonize IJ nematodes when cultured in vitro on agar. However, IJs did not develop or emerge from Galleria mellonella insects infected with secondary form. Unlike primary-form infected insects that were soft and flexible, secondary-form infected insects retained a rigid exoskeleton structure. Xenorhabdus bovienii primary and secondary form isolates are virulent towards Manduca sexta and several other insects. However, primary form stocks present attenuated virulence, suggesting that X. bovienii, like Xenorhabdus nematophila may undergo virulence modulation.

Threonine (Thr) is one of a few limiting essential amino acids (EAAs) in the animal feed industry, and its level in feed rations can impact production of important meat sources, such as swine and poultry. Threonine as well as EAAs lysine (Lys) and methionine (Met) are all synthesized via the aspartate family pathway. Here, we report a successful strategy to produce high free threonine soybean seed via identification of a feedback-resistant aspartate kinase (AK) enzyme that can be over-expressed in developing soybean seed. Towards this goal, we have purified and biochemically characterized AK from the enteric bacterium Xenorhabdus bovienii (Xb). Site-directed mutagenesis of XbAK identified two key regulatory residues Glu-257 and Thr-359 involved in lysine inhibition. Three feedback-resistant alleles, XbAK_T359I, XbAK_E257K and XbAK_E257K/T359I, have been generated. This study is the first to kinetically characterize the XbAK enzyme and provide biochemical and transgenic evidence that Glu-257 near the catalytic site is a critical residue for the allosteric regulation of AK. Furthermore, seed-specific expression of the feedback-resistant XbAK_T359I or XbAK_E257K allele results in increases of free Thr levels of up to 100-fold in R(1) soybean seed when compared to wild-type. Expression of feedback-sensitive wild-type AK did not substantially impact seed Thr content. In addition to high Thr, transgenic seed also showed substantial increases in other major free amino acid (FAA) levels, resulting in an up to 3.5-fold increase in the total FAA content. The transgenic seed was normal in appearance and germinated well under greenhouse conditions.

A postulate of systems biology is that, within the sum total of its genomic data, a network map of all functional interactions exists for each organism ([5][1]), and this “interaction map” can be assembled from that data agglomeration through the application of integration algorithms ([9][2], [

The symbiotic pathogenic bacterium Xenorhabdus nematophila produces two distinct intracellular inclusion bodies. The pixA gene, which encodes the 185-residue methionine-rich PixA inclusion body protein, was analyzed in the present study. The pixA gene was optimally expressed under stationary-phase conditions but its expression did not require RpoS. Analysis of a pixA mutant strain showed that PixA was not required for virulence towards the insect host or for colonization of or survival within the nematode host, and was not essential for nematode reproduction. The pixA gene was not present in the genome of Xenorhabdus bovienii , which also produces proteinaceous inclusions, indicating that PixA is specifically produced in X. nematophila .

Currently, Myxococcus xanthus is the only species of myxobacteria to mature into the postgenomic phase, and therefore, this chapter focuses on this organism. The chapter begins with a discussion with genomic screens using transposable elements, continues with the paralogous analysis of gene families, and concludes with whole-genome yeast two-hybrid (Y2H) assays. Availability of the M. xanthus genome sequence makes the paralogous characterization of gene families possible. Both of the findings presented in the chapter are significant because they represent how the application of a genome sequence can aid in the identification of genes involved in specific processes such as development and vegetative growth. All of the postgenomics methods presented in the chapter focus on using the genome sequence to facilitate the characterization of genes using experimental techniques. One of the ultimate goals of myxobacterial genomics is to determine the complete set of genetic networks that contribute to the complex life cycle of the myxobacteria. Making sense of these data will require the concerted effort of the whole community of those involved in myxobacterial research to pool resources and share data. Exploratory work using microarrays is currently under way in an effort to identify the temporal expression of genes under a host of different environmental conditions. This work has thus far yielded insight into the number of genes associated with transcriptional activators (TAs) under vegetative conditions.

While Barry's career has focused on the analysis of genomes, that mantra still rings in his head, and he has found that what is true for genetics is also true for genomics. The link between genotype and phenotype provided the framework for most genetics research prior to the genomics era. Biology and most sciences develop in a standard technology arc: description, simple testing, complex testing, and finally full-scale implementation. The first papers in the genomics field described the basic attributes of the technology and genome architecture. Microbiology is currently undergoing a revolution brought on by large-scale sequencing technology and the reality of dauntingly large data sets. Just 4 years after the Haemophilus genome was sequenced, forward-thinking biologists like Lee Hartwell began proposing a new paradigm, initially called “modular biology,” for genomics research. These ideas and frameworks have since launched the fields of systems biology, synthetic biology, and functional genomics. Much like the early genome papers in the mid-1990s, the earliest metagenomics papers also reported on our ability to sequence and make sense of environmental metagenomes. As a result, comparative metagenomics approaches may be the most fruitful approaches for analyzing such data sets. The authors also believe that only by studying the full continuum of metabolic and ecological interactions in a niche, and their reflection in the total genetic complement of all organisms of that niche (what we call the “nicheome”), can we understand the forces shaping microbial evolution.

Agrobacterium tumefaciens is a plant pathogen capable of transferring a defined segment of DNA to a host plant, generating a gall tumor. Replacing the transferred tumor-inducing genes with exogenous DNA allows the introduction of any desired gene into the plant. Thus, A. tumefaciens has been critical for the development of modern plant genetics and agricultural biotechnology. Here we describe the genome of A. tumefaciens strain C58, which has an unusual structure consisting of one circular and one linear chromosome. We discuss genome architecture and evolution and additional genes potentially involved in virulence and metabolic parasitism of host plants.

It is generally agreed that cytochrome c biogenesis requires that the apocytochrome and heme be transported separately to their site of function and assembly. In bacteria, this is outside the cytoplasmic membrane, whereby the apocytochromes c use sec-dependent signals for their translocation. Two different hypotheses have recently emerged as to how heme is exported: one involves an helABCD-encoded ATP binding cassette (ABC) transporter complex and the second does not. The second hypothesis concludes that an (HelAB)2 heterodimeric ABC transporter does not transport heme but some other substrate for cytochrome c biogenesis. The evidence supporting each of these two hypotheses and the role of this ABC transporter is discussed.

Oxidation-reduction titrations for the active-site disulfide/dithiol couples of the helX- and ccl2-encoded proteins involved in cytochrome c biogenesis in the purple non-sulfur bacterium Rhodobacter capsulatus have been carried out. The R. capsulatus HelX and Ccl2 proteins are predicted to function as part of a dithiol/disulfide cascade that reduces a disulfide on the apocytochromes c so that two cysteine thiols are available to form thioether linkages between the heme prosthetic group and the protein. Oxidation-reduction midpoint potential (E(m)) values, at pH 7.0, of -300 +/- 10 and -210 +/- 10 mV were measured for the HelX and Ccl2 (a soluble, truncated form of Ccl2) R. capsulatus proteins, respectively. Titrations of the disulfide/dithiol couple of a peptide designed to serve as a model for R. capsulatus apocytochrome c(2) have also been carried out, and an E(m) value of -170 +/- 10 mV was measured for the model peptide at pH 7.0. E(m) versus pH plots for HelX, Ccl2, and the apocytochrome c(2) model peptide were all linear over the pH range from 5.0 to 8.0, with the -59 mV/pH unit slope expected for a reaction in which two protons are taken up for each disulfide that is reduced. These results provide thermodynamic support for the proposal that HelX reduces Ccl2 and that reduced Ccl2, in turn, serves as the reductant for the production of the two thiols of the CysXxxYyyCysHis heme-binding motif of the apocytochromes.

Heme proteins play pivotal roles in a wealth of biological processes. Despite this, the molecular mechanisms by which heme traverses bilayer membranes for use in biosynthetic reactions are unknown. The biosynthesis of c -type cytochromes requires that heme is transported to the bacterial periplasm or mitochondrial intermembrane space where it is covalently ligated to two reduced cysteinyl residues of the apocytochrome. Results herein suggest that a family of integral membrane proteins in prokaryotes, protozoans, and plants act as transmembrane heme delivery systems for the biogenesis of c -type cytochromes. The complete topology of a representative from each of the three subfamilies was experimentally determined. Key histidinyl residues and a conserved tryptophan-rich region (designated the WWD domain) are positioned at the site of cytochrome c assembly for all three subfamilies. These histidinyl residues were shown to be essential for function in one of the subfamilies, an ABC transporter encoded by helABCD . We believe that a directed heme delivery pathway is vital for the synthesis of cytochromes c , whereby heme iron is protected from oxidation via ligation to histidinyl residues within the delivery proteins.

The photosynthetic bacterium Rhodobacter capsulatus synthesizes c-type cytochromes under a variety of growth conditions. For example, under aerobic growth, c-type cytochromes are synthesized as part of an electron transport pathway, using oxygen as the terminal electron acceptor. Anaerobically in the light, R. capsulatus requires cytochrome bc1 and other c-type cytochromes for the photosynthetic electron transport pathway. It is shown here that the ccl1 and ccl2 genes of R. capsulatus are required for the synthesis of all c-type cytochromes, including the cytochrome c' protein of unknown function but of structural similarity to cytochrome b562. Polar and nonpolar mutations constructed in each gene demonstrated that the ccl12 genes form an operon. Expression of the ccl12 genes was examined by using lacZ and phoA fusions as translational reporters. Primer extension analysis was used to determine transcriptional control and the start site of the ccl12 promoter. Finally, antiserum to the Ccl2 protein was used to quantitate levels of Ccl2 under six different growth conditions. The Ccl2 protein is present at 20-fold-higher levels under conditions where oxygen is present. In contrast, other cytochromes c biogenesis proteins, HelA and HelX, previously shown to be part of an helABCDX operon, are at relatively similar levels under these six growth conditions. This discovery is discussed in terms of the physiology and evolution of cytochromes c biogenesis, with particular attention to oxidative environments.

Background Western corn rootworm (WCR) is one of the most significant insect pests of maize in North America. WCR has dramatically increased its range in the last century, invading key maize production areas in the US and abroad. In addition, this species has a history of evolving traits that allow it to escape various control options. Improved genetic and genomic resources are crucial tools for understanding population history and the genetic basis of trait evolution. Here we produce and analyze a transcriptome assembly for WCR. We also perform whole genome population resequencing, and combine these resources to better understand the evolutionary history of WCR. Results The WCR transcriptome assembly presented here contains approximately 16,000 unigenes, many of which have high similarity to other insect species. Among these unigenes we found several gene families that have been implicated in insecticide resistance in other species. We also identified over 500,000 unigene based SNPs among 26 WCR populations. We used these SNPs to scan for outliers among the candidate genes, and to understand how population processes have shaped genetic variation in this species. Conclusions This study highlights the utility of transcriptomic and genomic resources as foundational tools for dealing with highly adaptive pest species. Using these tools we identified candidate gene families for insecticide resistance and reveal aspects of WCR population history in light of the species’ recent range expansion. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-195) contains supplementary material, which is available to authorized users.

The cydDC operon of Escherichia coli encodes an ATP-dependent transporter of unknown function that is required for cytochrome bd synthesis. Strains containing defects in either the cydD or cydC gene also demonstrate hypersensitivity to growth at high temperatures and the inability to exit the stationary phase at 37 degrees C. We wished to determine what is responsible for these hypersensitive phenotypes and whether they are due to a lack of the CydDC proteins or a defect of the cytochrome bd encoded by the cydAB genes. Using both K-12- and B-type strains of E. coli, we have compared the phenotypes of isogenic cydAB mutants and cydC mutants. In both K-12- and B-type backgrounds, the hypersensitive phenotypes are due to defects of cytochrome bd activity and not defects of the cydDC genes. We also found that the temperature-sensitive growth phenotypes can be suppressed by exogenous reducing agents, such as glutathione and cysteine. Strikingly, even the enzymes catalase and superoxide dismutase, when added exogenously, can correct the temperature-sensitive and stationary phase arrest phenotypes. We propose that the temperature sensitive growth phenotypes are due to a buildup of diffusible oxygen radicals brought on by the absence of cytochrome bd.

Klebsiella pneumoniae can use nitrate and nitrite as sole nitrogen sources through the nitrate assimilation pathway. We previously identified structural genes for assimilatory nitrate and nitrite reductases, nasA and nasB, respectively. We report here our further identification of four genes, nasFEDC, upstream of the nasBA genes. The nasFEDCBA genes probably form an operon. Mutational and complementation analyses indicated that both the nasC and nasA genes are required for nitrate assimilation. The predicted NASC protein is homologous to a variety of NADH-dependent oxidoreductases. Thus, the NASC protein probably mediates electron transfer from NADH to the NASA protein, which contains the active site for nitrate reduction. The deduced NASF, NASE, and NASD proteins are homologous to the NRTA, NRTB, and NRTD proteins, respectively, that are involved in nitrate uptake in Synechococcus sp. (T. Omata, X. Andriesse, and A. Hirano, Mol. Gen. Genet. 236:193-202, 1993). Mutational and complementation studies indicated that the nasD gene is required for nitrate but not nitrite assimilation. By analogy with the Synechococcus nrt genes, we propose that the nasFED genes are involved in nitrate transport in K. pneumoniae.

Klebsiella pneumoniae can use nitrate and nitrite as sole nitrogen sources during aerobic growth. Assimilatory nitrate and nitrite reductases convert nitrate through nitrite to ammonium. We report here the molecular cloning of the nasA and nasB genes, which encode assimilatory nitrate and nitrite reductase, respectively. These genes are tightly linked and probably form a nasBA operon. In vivo protein expression and DNA sequence analysis revealed that the nasA and nasB genes encode 92- and 104-kDa proteins, respectively. The NASA polypeptide is homologous to other prokaryotic molybdoenzymes, and the NASB polypeptide is homologous to eukaryotic and prokaryotic NADH-nitrite reductases. The narL gene product positively regulates expression of the structural genes for respiratory nitrate reductase, narGHJI. Surprisingly, we found that the nasBA operon is tightly linked to the narL-narGHJI region in K. pneumoniae, even though the nitrate assimilatory and respiratory enzymes serve different physiological functions.

Siroheme, a cofactor of both sulfite and nitrite reductase in Salmonella typhimurium, requires the cysG gene for its synthesis. Three steps are required to synthesize siroheme from uroporphyrinogen III, the last common intermediate in the heme and siroheme pathways. All previously characterized cysG mutants were shown to be defective for the synthesis of cobalamin (B12), which shares a common precursor with siroheme. Since few cysG auxotrophs had been previously analyzed and since there is no evidence of siroheme mutants outside of the cysG region, we sought to expand the analysis of the region by isolating more mutations and studying the transcriptional regulation of the cysG gene using lacZ fusions. We isolated and analyzed 66 cysG auxotrophs. All were defective for both siroheme and cobalamin synthesis. Five exceptional mutants were partially defective for the synthesis of both and appear to be leaky. Complementation tests with tandem duplications suggest that the mutations causing the Cys auxotrophy affect only one cistron. The cysG gene is transcribed in a clockwise direction; this was demonstrated by a method that permits determining the orientation of two genes of unknown orientation provided their relative map order is known. The cysG gene was not part of the cysteine regulon, but had a substantial basal level of expression which was induced fivefold when cells were grown anaerobically on nitrite. Finally, we used Mud-generated duplications to genetically determine the organization of the cysG and nirB genes.

Department of Plant Pathology and Microbiology, University of Arizona, Tucson, AZ, USA; Department of Biochemistry, Virginia Tech, Blacksburg, VA, USA; Hiram College, Hiram, OH, USA; Monsanto Company, St. Louis, MO, USA; Virginia Bioinformatics Institute, Virginia Tech, Blacksburg, VA, USA; Department of Applied Biological Sciences, Arizona State University, Tempe, AZ, USA; Department of Biology, Seattle Pacific University, Seattle, WA, USA

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced δ-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell–cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Accurate determination of functional interactions among proteins at the genome level remains a challenge for genomic research. Here we introduce a genome-scale approach to functional protein annotation--phylogenomic mapping--that requires only sequence data, can be applied equally well to both finished and unfinished genomes, and can be extended beyond single genomes to annotate multiple genomes simultaneously. We have developed and applied it to more than 200 sequenced bacterial genomes. Proteins with similar evolutionary histories were grouped together, placed on a three dimensional map and visualized as a topographical landscape. The resulting phylogenomic maps display thousands of proteins clustered in mountains on the basis of coinheritance, a strong indicator of shared function. In addition to systematic computational validation, we have experimentally confirmed the ability of phylogenomic maps to predict both mutant phenotype and gene function in the delta proteobacterium Myxococcus xanthus.

A search of the M1genome sequence, which includes 97% of theMyxococcus xanthusgenes, identified 53 sequence homologs of σ54-dependent enhancer binding proteins (EBPs). A DNA microarray was constructed from the M1genome that includes those homologs and 318 otherM. xanthusgenes for comparison. To screen the developmental program with this array, an RNA extract from growing cells was compared with one prepared from developing cells at 12 h. Previous reporter studies had shown thatM. xanthushas initiated development and has begun to express many developmentally regulated genes by 12 h. The comparison revealed substantial increases in the expression levels of 11 transcription factors that may respond to environmental stimuli. Six of the 53 EBP homologs were expressed at significantly higher levels at 12 h of development than during growth. Three were previously unknown genes, and they were inactivated to look for effects on fruiting body development. One knockout mutant produced fruiting bodies of abnormal shape that depended on the composition of the medium.

An analysis of thirty-three genomes of selected bacteria for the presence of specific respiratory pathways and cytochrome c biogenesis systems has led to observations on respiration and biogenesis. A table summarizing these results is presented. The data suggested that Bordetella pertussis would be an excellent genetic model to study the System II cytochrome c biogenesis pathway. These observations are discussed and the results of genetic studies on System II biogenesis in B. pertussis are presented as a case for the power of comparative genomics. System II is present in organisms as diverse as Helicobacter, Neisseria, Porphyromonas, mycobacteria, cyanobacteria, and plants (chloroplasts), indicating this pathway's prominence and that horizontal transfer of system II (and/or System I) must have occurred on multiple occasions.

The past 10 years have heralded remarkable progress in the understanding of the biogenesis of c-type cytochromes. The hallmark of c-type cytochrome synthesis is the covalent ligation of haem vinyl groups to two cysteinyl residues of the apocytochrome (at a Cys-Xxx-Yyy-Cys-His signature motif). From genetic, genomic and biochemical studies, it is clear that three distinct systems have evolved in nature to assemble this ancient protein. In this review, common principles of assembly for all systems and the molecular mechanisms predicted for each system are summarized. Prokaryotes, plant mitochondria and chloroplasts use either system I or II, which are each predicted to use dedicated mechanisms for haem delivery, apocytochrome ushering and thioreduction. Accessory proteins of systems I and II co-ordinate the positioning of these two substrates at the membrane surface for covalent ligation. The third system has evolved specifically in mitochondria of fungi, invertebrates and vertebrates. For system III, a pivotal role is played by an enzyme called cytochrome c haem lyase (CCHL) in the mitochondrial intermembrane space.

Strains of Escherichia coli containing mutations in the cydDC genes are defective for synthesis of the heme proteins cytochrome bd and c-type cytochromes. The cydDC genes encode a putative heterodimeric ATP-binding cassette transporter that has been proposed to act as an exporter of heme to the periplasm. To more fully understand the role of this transporter (and other factors) in heme protein biosynthesis, we developed plasmids that produce various heme proteins (e.g., cytochrome b5, cytochrome b562, and hemoglobin) in the periplasm of E. coli. By using these reporters, it was shown that the steady-state levels of polypeptides of heme proteins known to be stable without heme (e.g., cytochrome b5 and hemoglobin apoprotein) are significantly reduced in a cydC mutant. Exogenous addition of hemin to the cydC mutant still resulted in < 10% of wild-type steady-state levels of apohemoglobin in the periplasm. Since the results of heme reporter studies are not consistent with lower heme availability (i.e., heme export) in a cydC mutant, we analyzed other properties of the periplasm in cydC mutants and compared them with those of the periplasm in cydAB (encoding cytochrome bd) mutants and wild-type cells. Our results led us to favor a hypothesis whereby cydDC mutants are defective in the reduction environment within the periplasmic space. Such an imbalance could lead to defects in the synthesis of heme-liganded proteins. The heme reporters were also used to analyze strains of E. coli with a defect in genes encoding homologs of a different ABC transporter (helABC). The helABC genes have previously been shown to be required for the assembly of c-type cytochromes in Rhodobacter capsulatus (R. G. Kranz, J. Bacteriol. 171:456-464, 1989; D. L. Beckman, D. R. Trawick, and R. G. Kranz, Genes Dev. 6:268-283, 1992). This locus was shown to be essential in E. coli for endogenous cytochrome c biogenesis but not cytochrome b562 synthesis. Consistent with these and previous results, it is proposed that the HelABC transporter is specifically involved in heme export for ligation (hel). This class of periplasmic cytochromes is proposed to require heme liganding before undergoing correct folding.

Klebsiella pneumoniae can use nitrate and nitrite as sole nitrogen sources through the nitrate assimilatory pathway. The structural genes for assimilatory nitrate and nitrite reductases together with genes necessary for nitrate transport form an operon, nasFEDCBA. Expression of the nasF operon is regulated both by general nitrogen control and also by nitrate or nitrite induction. We have identified a gene, nasR, that is necessary for nitrate and nitrite induction. The nasR gene, located immediately upstream of the nasFEDCBA operon, encodes a 44-kDa protein. The NasR protein shares carboxyl-terminal sequence similarity with the AmiR protein of Pseudomonas aeruginosa, the positive regulator of amiE (aliphatic amidase) gene expression. In addition, we present evidence that the nasF operon is not autogenously regulated.

Background Lateral gene transfer (LGT) is thought to promote speciation in bacteria, though well-defined examples have not been put forward. Methodology/principle findings We examined the evolutionary history of the genes essential for a trait that defines a phylogenetic order, namely fruiting body development of the Myxococcales. Seventy-eight genes that are essential for Myxococcus xanthus development were examined for LGT. About 73% of the genes exhibit a phylogeny similar to that of the 16S rDNA gene and a codon bias consistent with other M. xanthus genes suggesting vertical transmission. About 22% have an altered codon bias and/or phylogeny suggestive of LGT. The remaining 5% are unique. Genes encoding signal production and sensory transduction were more likely to be transmitted vertically with clear examples of duplication and divergence into multigene families. Genes encoding metabolic enzymes were frequently acquired by LGT. Myxobacteria exhibit aerobic respiration unlike most of the delta Proteobacteria. M. xanthus contains a unique electron transport pathway shaped by LGT of genes for succinate dehydrogenase and three cytochrome oxidase complexes. Conclusions/significance Fruiting body development depends on genes acquired by LGT, particularly those involved in polysaccharide production. We suggest that aerobic growth fostered innovation necessary for development by allowing myxobacteria access to a different gene pool from anaerobic members of the delta Proteobacteria. Habitat destruction and loss of species diversity could restrict the evolution of new bacterial groups by limiting the size of the prospective gene pool.

Automated DNA sequencing technology is so rapid that analysis has become the rate-limiting step. Hundreds of prokaryotic genome sequences are publicly available, with new genomes uploaded at the rate of approximately 20 per month. As a result, this growing body of genome sequences will include microorganisms not previously identified, isolated, or observed. We hypothesize that evolutionary pressure exerted by an ecological niche selects for a similar genetic repertoire in those prokaryotes that occupy the same niche, and that this is due to both vertical and horizontal transmission. To test this, we have developed a novel method to classify prokaryotes, by calculating their Pfam protein domain distributions and clustering them with all other sequenced prokaryotic species. Clusters of organisms are visualized in two dimensions as ‘mountains’ on a topological map. When compared to a phylogenetic map constructed using 16S rRNA, this map more accurately clusters prokaryotes according to functional and environmental attributes. We demonstrate the ability of this map, which we term a “niche map”, to cluster according to ecological niche both quantitatively and qualitatively, and propose that this method be used to associate uncharacterized prokaryotes with their ecological niche as a means of predicting their functional role directly from their genome sequence.