Whole-genome comparison clarifies close phylogenetic relationships between the phyla Dictyoglomi and Thermotogae

Identification of a deep-branching thermophilic clade sheds light on early bacterial evolution

It has been proposed that early bacteria, or even the last universal common ancestor of all cells, were thermophilic. However, research on the origin and evolution of thermophily is hampered by the difficulties associated with the isolation of deep-branching thermophilic microorganisms in pure culture. Here, we isolate a deep-branching thermophilic bacterium from a deep-sea hydrothermal vent, using a two-step cultivation strategy ("Subtraction-Suboptimal", StS) designed to isolate rare organisms. The bacterium, which we name Zhurongbacter thermophilus 3DAC, is a sulfur-reducing heterotroph that is phylogenetically related to Coprothermobacterota and other thermophilic bacterial groups, forming a clade that seems to represent a major, early-diverging bacterial lineage. The ancestor of this clade might be a thermophilic, strictly anaerobic, motile, hydrogen-dependent, and mixotrophic bacterium. Thus, our study provides insights into the early evolution of thermophilic bacteria.

Disentangling the lipid divide: Identification of key enzymes for the biosynthesis of membrane-spanning and ether lipids in Bacteria

Bacterial membranes are composed of fatty acids (FAs) ester-linked to glycerol-3-phosphate, while archaea have membranes made of isoprenoid chains ether-linked to glycerol-1-phosphate. Many archaeal species organize their membrane as a monolayer of membrane-spanning lipids (MSLs). Exceptions to this "lipid divide" are the production by some bacterial species of (ether-bound) MSLs, formed by tail-to-tail condensation of FAs resulting in the formation of (iso) diabolic acids (DAs), which are the likely precursors of paleoclimatological relevant branched glycerol dialkyl glycerol tetraether molecules. However, the enzymes responsible for their production are unknown. Here, we report the discovery of bacterial enzymes responsible for the condensation reaction of FAs and for ether bond formation and confirm that the building blocks of iso-DA are branched iso-FAs. Phylogenomic analyses of the key biosynthetic genes reveal a much wider diversity of potential MSL (ether)-producing bacteria than previously thought, with importantt implications for our understanding of the evolution of lipid membranes.

Metaproteomics reveals enzymatic strategies deployed by anaerobic microbiomes to maintain lignocellulose deconstruction at high solids

Economically viable production of cellulosic biofuels requires operation at high solids loadings—on the order of 15 wt%. To this end we characterize Nature’s ability to deconstruct and utilize mid-season switchgrass at increasing solid loadings using an anaerobic methanogenic microbiome. This community exhibits undiminished fractional carbohydrate solubilization at loadings ranging from 30 g/L to 150 g/L. Metaproteomic interrogation reveals marked increases in the abundance of specific carbohydrate-active enzyme classes. Significant enrichment of auxiliary activity family 6 enzymes at higher solids suggests a role for Fenton chemistry. Stress-response proteins accompanying these reactions are similarly upregulated at higher solids, as are β-glucosidases, xylosidases, carbohydrate-debranching, and pectin-acting enzymes—all of which indicate that removal of deconstruction inhibitors is important for observed undiminished solubilization. Our work provides insights into the mechanisms by which natural microbiomes effectively deconstruct and utilize lignocellulose at high solids loadings, informing the future development of defined cultures for efficient bioconversion.

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Multidomain ribosomal protein trees and the planctobacterial origin of neomura (eukaryotes, archaebacteria)

Palaeontologically, eubacteria are > 3× older than neomura (eukaryotes, archaebacteria). Cell biology contrasts ancestral eubacterial murein peptidoglycan walls and derived neomuran N-linked glycoprotein coats/walls. Misinterpreting long stems connecting clade neomura to eubacteria on ribosomal sequence trees (plus misinterpreted protein paralogue trees) obscured this historical pattern. Universal multiprotein ribosomal protein (RP) trees, more accurate than rRNA trees, are taxonomically undersampled. To reduce contradictions with genically richer eukaryote trees and improve eubacterial phylogeny, we constructed site-heterogeneous and maximum-likelihood universal three-domain, two-domain, and single-domain trees for 143 eukaryotes (branching now congruent with 187-protein trees), 60 archaebacteria, and 151 taxonomically representative eubacteria, using 51 and 26 RPs. Site-heterogeneous trees greatly improve eubacterial phylogeny and higher classification, e.g. showing gracilicute monophyly, that many 'rDNA-phyla' belong in Proteobacteria, and reveal robust new phyla Synthermota and Aquithermota. Monoderm Posibacteria and Mollicutes (two separate wall losses) are both polyphyletic: multiple outer membrane losses in Endobacteria occurred separately from Actinobacteria; neither phylum is related to Chloroflexi, the most divergent prokaryotes, which originated photosynthesis (new model proposed). RP trees support an eozoan root for eukaryotes and are consistent with archaebacteria being their sisters and rooted between Filarchaeota (=Proteoarchaeota, including 'Asgardia') and Euryarchaeota sensu-lato (including ultrasimplified 'DPANN' whose long branches often distort trees). Two-domain trees group eukaryotes within Planctobacteria, and archaebacteria with Planctobacteria/Sphingobacteria. Integrated molecular/palaeontological evidence favours negibacterial ancestors for neomura and all life. Unique presence of key pre-neomuran characters favours Planctobacteria only as ancestral to neomura, which apparently arose by coevolutionary repercussions (explained here in detail, including RP replacement) of simultaneous outer membrane and murein loss. Planctobacterial C-1 methanotrophic enzymes are likely ancestral to archaebacterial methanogenesis and β-propeller-α-solenoid proteins to eukaryotic vesicle coats, nuclear-pore-complexes, and intraciliary transport. Planctobacterial chaperone-independent 4/5-protofilament microtubules and MamK actin-ancestors prepared for eukaryote intracellular motility, mitosis, cytokinesis, and phagocytosis. We refute numerous wrong ideas about the universal tree.

Proposal for a new classification of a deep branching bacterial phylogenetic lineage: transfer of Coprothermobacter proteolyticus and Coprothermobacter platensis to Coprothermobacteraceae fam. nov., within Coprothermobacterales ord. nov., Coprothermobacteria classis nov. and Coprothermobacterota phyl. nov. and emended description of the family Thermodesulfobiaceae

The genus Coprothermobacter (initially named Thermobacteroides) is currently placed within the phylum Firmicutes. Early 16S rRNA gene based phylogenetic studies pointed out the great differences between Coprothermobacter and other members of the Firmicutes, revealing that it constitutes a new deep branching lineage. Over the years, several studies based on 16S rRNA gene and whole genome sequences have indicated that Coprothermobacter is very distant phylogenetically to all other bacteria, supporting its placement in a distinct deeply rooted novel phylum. In view of this, we propose its allocation to the new family Coprothermobacteraceae within the novel order Coprothermobacterales, the new class Coprothermobacteria, and the new phylum Coprothermobacterota, and an emended description of the family Thermodesulfobiaceae.

Draft Genome Sequence of a Dictyoglomus sp. from an Enrichment Culture of a New Zealand Geothermal Spring

A draft genome of a novel Dictyoglomus sp., NZ13-RE01, was obtained from a New Zealand hot spring enrichment culture. The 1,927,012-bp genome is similar in both size and G+C content to other Dictyoglomus spp. Like its relatives, Dictyoglomus sp. NZ13-RE01 encodes many genes involved in complex carbohydrate metabolism.

Enhanced features of Dictyoglomus turgidum Cellulase A engineered with carbohydrate binding module 11 from Clostridium thermocellum

Lignocellulosic biomass (LCB) is a low-cost and abundant source of fermentable sugars. Enzymatic hydrolysis is one of the main ways to obtain sugars from biomass, but most of the polysaccharide-degrading enzymes are poorly efficient on LCB and cellulases with higher performances are required. In this study, we designed a chimeric protein by adding the carbohydrate binding module (CBM) of the cellulosomal enzyme CtLic26A-Cel5E (endoglucanase H or CelH) from Clostridium (Ruminiclostridium) thermocellum to the C-terminus of Dtur CelA, an interesting hyperthermostable endoglucanase from Dictyoglomus turgidum. The activity and binding rate of both native and chimeric enzyme were evaluated on soluble and insoluble polysaccharides. The addition of a CBM resulted in a cellulase with enhanced stability at extreme pHs, higher affinity and activity on insoluble cellulose.

A review of bioinformatics platforms for comparative genomics. Recent developments of the EDGAR 2.0 platform and its utility for taxonomic and phylogenetic studies

The rapid development of next generation sequencing technology has greatly increased the amount of available microbial genomes. As a result of this development, there is a rising demand for fast and automated approaches in analyzing these genomes in a comparative way. Whole genome sequencing also bears a huge potential for obtaining a higher resolution in phylogenetic and taxonomic classification. During the last decade, several software tools and platforms have been developed in the field of comparative genomics. In this manuscript, we review the most commonly used platforms and approaches for ortholog group analyses with a focus on their potential for phylogenetic and taxonomic research. Furthermore, we describe the latest improvements of the EDGAR platform for comparative genome analyses and present recent examples of its application for the phylogenomic analysis of different taxa. Finally, we illustrate the role of the EDGAR platform as part of the BiGi Center for Microbial Bioinformatics within the German network on Bioinformatics Infrastructure (de.NBI).

Evolutionary Relationships and Taxa-Specific Conserved Signature Indels Among Cellulases of Archaea, Bacteria, and Eukarya

The cellulases from different cellulolytic organisms have evolutionary relationships, which range from single-celled prokaryotes to the complex eukaryotes of the living world. This in silico analysis revealed the presence of a conserved cellulase domain along with evolutionary relationships among cellulases from several species of Archaea, Bacteria, and Eukarya. The amino acid sequences of cellulases from Archaea and Bacteria showed closer identity with their domain or phylum members that provided insights into convergent and divergent evolution of cellulases from other enzymes with different substrate specificities. Evolutionary relatedness was also observed in phylogenetic trees among a number of cellulase sequences of diverse taxa. In cellulases, propensity for alanine, glycine, leucine, serine, and threonine was high, but low for cysteine, histidine, and methionine. Catalytic aspartic acid had a higher propensity than glutamic acid, and both were involved in regular expression patterns. Characteristic group and multigroup-specific conserved signature indels located in the catalytic domains of cellulases were observed that further clarified evolutionary relationships. These indels can be distinctive molecular tools for understanding phylogeny and identification of unknown cellulolytic species of common evolutionary descent in different environments.

The Complete Genome Sequence of Hyperthermophile Dictyoglomus turgidum DSM 6724™ Reveals a Specialized Carbohydrate Fermentor

Here we report the complete genome sequence of the chemoorganotrophic, extremely thermophilic bacterium, Dictyoglomus turgidum, which is a Gram negative, strictly anaerobic bacterium. D. turgidum and D. thermophilum together form the Dictyoglomi phylum. The two Dictyoglomus genomes are highly syntenic, and both are distantly related to Caldicellulosiruptor spp. D. turgidum is able to grow on a wide variety of polysaccharide substrates due to significant genomic commitment to glycosyl hydrolases, 16 of which were cloned and expressed in our study. The GH5, GH10, and GH42 enzymes characterized in this study suggest that D. turgidum can utilize most plant-based polysaccharides except crystalline cellulose. The DNA polymerase I enzyme was also expressed and characterized. The pure enzyme showed improved amplification of long PCR targets compared to Taq polymerase. The genome contains a full complement of DNA modifying enzymes, and an unusually high copy number (4) of a new, ancestral family of polB type nucleotidyltransferases designated as MNT (minimal nucleotidyltransferases). Considering its optimal growth at 72°C, D. turgidum has an anomalously low G+C content of 39.9% that may account for the presence of reverse gyrase, usually associated with hyperthermophiles.

Phylogeny Reconstruction with Alignment-Free Method That Corrects for Horizontal Gene Transfer

Advances in sequencing have generated a large number of complete genomes. Traditionally, phylogenetic analysis relies on alignments of orthologs, but defining orthologs and separating them from paralogs is a complex task that may not always be suited to the large datasets of the future. An alternative to traditional, alignment-based approaches are whole-genome, alignment-free methods. These methods are scalable and require minimal manual intervention. We developed SlopeTree, a new alignment-free method that estimates evolutionary distances by measuring the decay of exact substring matches as a function of match length. SlopeTree corrects for horizontal gene transfer, for composition variation and low complexity sequences, and for branch-length nonlinearity caused by multiple mutations at the same site. We tested SlopeTree on 495 bacteria, 73 archaea, and 72 strains of Escherichia coli and Shigella. We compared our trees to the NCBI taxonomy, to trees based on concatenated alignments, and to trees produced by other alignment-free methods. The results were consistent with current knowledge about prokaryotic evolution. We assessed differences in tree topology over different methods and settings and found that the majority of bacteria and archaea have a core set of proteins that evolves by descent. In trees built from complete genomes rather than sets of core genes, we observed some grouping by phenotype rather than phylogeny, for instance with a cluster of sulfur-reducing thermophilic bacteria coming together irrespective of their phyla. The source-code for SlopeTree is available at: http://prodata.swmed.edu/download/pub/slopetree_v1/slopetree.tar.gz.

Due to their lack of distinct morphological features, bacteria and archaea were extremely difficult to classify until technology was developed to obtain their DNA sequences; these sequences could then be compared to estimate evolutionary relationships. Now, due to technological advances, there is a flood of available sequences from a wide variety of organisms. These advances have spurred the development of algorithms which can estimate evolutionary relationships using whole genomes, in contrast to the more traditional methods which used single genes earlier and now typically use groups of conserved genes. However, there are many challenges when attempting to infer evolutionary relationships, in particular horizontal gene transfer, where DNA is transferred from one organism to another, resulting in an organism’s genome containing DNA that does not reflect its evolution by descent. We developed a new whole-genome method for estimating evolutionary distances which identifies and corrects for horizontal transfer. We found that for SlopeTree and all other whole-genome methods we applied, horizontal transfer causes some evolutionary distances to be grossly underestimated, and that our correction corrects for this.

Exploring membrane respiratory chains

Acquisition of energy is central to life. In addition to the synthesis of ATP, organisms need energy for the establishment and maintenance of a transmembrane difference in electrochemical potential, in order to import and export metabolites or to their motility. The membrane potential is established by a variety of membrane bound respiratory complexes. In this work we explored the diversity of membrane respiratory chains and the presence of the different enzyme complexes in the several phyla of life. We performed taxonomic profiles of the several membrane bound respiratory proteins and complexes evaluating the presence of their respective coding genes in all species deposited in KEGG database. We evaluated 26 quinone reductases, 5 quinol:electron carriers oxidoreductases and 18 terminal electron acceptor reductases. We further included in the analyses enzymes performing redox or decarboxylation driven ion translocation, ATP synthase and transhydrogenase and we also investigated the electron carriers that perform functional connection between the membrane complexes, quinones or soluble proteins. Our results bring a novel, broad and integrated perspective of membrane bound respiratory complexes and thus of the several energetic metabolisms of living systems. This article is part of a Special Issue entitled 'EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2-6, 2016', edited by Prof. Paolo Bernardi.

Evolutionary relationships of completely sequenced Clostridia species and close relatives

The class Clostridia in the phylum Firmicutes includes a very heterogeneous assemblage of bacteria. Their evolutionary relationships are not well established; revisions of their phylogenetic placements based on comparative studies of 16S rRNA gene sequences are in progress as genome sequence information accumulates. In this work, phylogenetic trees were reconstructed based on 21 concatenated ribosomal protein sequences using Bayesian and maximum-likelihood methods. Both trees consistently indicate that the Halanaerobiales is a deeply branching order among the class Clostridia. The rest of the Clostridia species are grouped into 10 monophyletic clusters, most of which are comprised of two or three orders and families according to the current Clostridial taxonomy. The maximum-likelihood tree placed Coprothermobacter proteolyticus and Thermodesulfobium narugense in the class Clostridia in accordance with the current taxonomy, in which these two bacteria are assigned to the family Thermodesulfobiaceae. However, the Bayesian tree placed these two bacteria at the boundary between the Firmicutes and Actinobacteria. A gene arrangement that is present uniquely in the Firmicutes species was identified. Both Coprothermobacter proteolyticus and Thermodesulfobium narugense do not have this arrangement characteristic of the Firmicutes. On the basis of the Bayesian tree and gene arrangement comparison, it is suggested that Coprothermobacter proteolyticus and Thermodesulfobium narugense should be placed outside the phylum Firmicutes.

Thermophilic anaerobic digestion of thermal pretreated sludge: role of microbial community structure and correlation with process performances

Thermal hydrolysis pretreatment coupled with Thermophilic Anaerobic Digestion (TAD) for Waste Activated Sludge (WAS) treatment is a promising combination to improve biodegradation kinetics during stabilization. However, to date there is a limited knowledge of the anaerobic biomass composition and its impact on TAD process performances. In this study, the structure and dynamics of the microbial communities selected in two semi-continuous anaerobic digesters, fed with untreated and thermal pretreated sludge, were investigated. The systems were operated for 250 days at different organic loading rate. 16S rRNA gene clonal analysis and Fluorescence In Situ Hybridization (FISH) analyses allowed us to identify the majority of bacterial and archaeal populations. Proteolytic Coprothermobacter spp. and hydrogenotrophic Methanothermobacter spp. living in strict syntrophic association were found to dominate in TAD process. The establishment of a syntrophic proteolytic pathway was favoured by the high temperature of the process and enhanced by the thermal pretreatment of the feeding sludge. Proteolytic activity, alone or with thermal pretreatment, occurred during TAD as proven by increasing concentration of soluble ammonia and soluble COD (sCOD) during the process. However, the availability of a readily biodegradable substrate due to pretreatment allowed to significant sCOD removals (more than 55%) corresponding to higher biogas production in the reactor fed with thermal pretreated sludge. Microbial population dynamics analysed by FISH showed that Coprothermobacter and Methanothermobacter immediately established a stable syntrophic association in the reactor fed with pretreated sludge in line with the overall improved TAD performances observed under these conditions.

Mesoaciditoga lauensis gen. nov., sp. nov., a moderately thermoacidophilic member of the order Thermotogales from a deep-sea hydrothermal vent

A novel moderately thermophilic, heterotrophic bacterium was isolated from a deep-sea hydrothermal vent deposit from the Mariner field along the Eastern Lau Spreading Center of the south-western Pacific Ocean. Cells were short motile rods (about 0.4×0.8 µm) that occurred singly or in pairs and were surrounded by a sheath-like membrane or 'toga'. The cells grew between 45 and 65 °C (optimum 57-60 °C) and at pH 4.1-6.0 (optimum pH 5.5-5.7) and grew optimally at 3 % (w/v) NaCl. The isolate grew on a range of carbon and proteinaceous substrates and reduced sulfur. The G+C content of the DNA was about 45 mol%. Phylogenetic analysis of the 16S rRNA gene sequence placed the new isolate as a deeply diverging lineage within the order Thermotogales. Based on the physiological, morphological and phylogenetic data, the isolate represents a novel species of a new genus with the proposed name Mesoaciditoga lauensis gen. nov., sp. nov. The type strain of Mesoaciditoga lauensis is cd-1655R(T) ( = DSM 25116(T) = OCM 1212(T)).

Phylogeny of bacterial and archaeal genomes using conserved genes: supertrees and supermatrices

Over 3000 microbial (bacterial and archaeal) genomes have been made publically available to date, providing an unprecedented opportunity to examine evolutionary genomic trends and offering valuable reference data for a variety of other studies such as metagenomics. The utility of these genome sequences is greatly enhanced when we have an understanding of how they are phylogenetically related to each other. Therefore, we here describe our efforts to reconstruct the phylogeny of all available bacterial and archaeal genomes. We identified 24, single-copy, ubiquitous genes suitable for this phylogenetic analysis. We used two approaches to combine the data for the 24 genes. First, we concatenated alignments of all genes into a single alignment from which a Maximum Likelihood (ML) tree was inferred using RAxML. Second, we used a relatively new approach to combining gene data, Bayesian Concordance Analysis (BCA), as implemented in the BUCKy software, in which the results of 24 single-gene phylogenetic analyses are used to generate a "primary concordance" tree. A comparison of the concatenated ML tree and the primary concordance (BUCKy) tree reveals that the two approaches give similar results, relative to a phylogenetic tree inferred from the 16S rRNA gene. After comparing the results and the methods used, we conclude that the current best approach for generating a single phylogenetic tree, suitable for use as a reference phylogeny for comparative analyses, is to perform a maximum likelihood analysis of a concatenated alignment of conserved, single-copy genes.

PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes

New microbial genomes are constantly being sequenced, and it is crucial to accurately determine their taxonomic identities and evolutionary relationships. Here we report PhyloPhlAn, a new method to assign microbial phylogeny and putative taxonomy using >400 proteins optimized from among 3,737 genomes. This method measures the sequence diversity of all clades, classifies genomes from deep-branching candidate divisions through closely related subspecies and improves consistency between phylogenetic and taxonomic groupings. PhyloPhlAn improved taxonomic accuracy for existing and newly sequenced genomes, detecting 157 erroneous labels, correcting 46 and placing or refining 130 new genomes. We provide examples of accurate classifications from subspecies (Sulfolobus spp.) to phyla, and of preliminary rooting of deep-branching candidate divisions, including consistent statistical support for Caldiserica (formerly candidate division OP5). PhyloPhlAn will thus be useful for both phylogenetic assessment and taxonomic quality control of newly sequenced genomes. The final phylogenies, conserved protein sequences and open-source implementation are available online.

Genomic determinants of sporulation in Bacilli and Clostridia: towards the minimal set of sporulation-specific genes

Three classes of low-G+C Gram-positive bacteria (Firmicutes), Bacilli, Clostridia and Negativicutes, include numerous members that are capable of producing heat-resistant endospores. Spore-forming firmicutes include many environmentally important organisms, such as insect pathogens and cellulose-degrading industrial strains, as well as human pathogens responsible for such diseases as anthrax, botulism, gas gangrene and tetanus. In the best-studied model organism Bacillus subtilis, sporulation involves over 500 genes, many of which are conserved among other bacilli and clostridia. This work aimed to define the genomic requirements for sporulation through an analysis of the presence of sporulation genes in various firmicutes, including those with smaller genomes than B. subtilis. Cultivable spore-formers were found to have genomes larger than 2300 kb and encompass over 2150 protein-coding genes of which 60 are orthologues of genes that are apparently essential for sporulation in B. subtilis. Clostridial spore-formers lack, among others, spoIIB, sda, spoVID and safA genes and have non-orthologous displacements of spoIIQ and spoIVFA, suggesting substantial differences between bacilli and clostridia in the engulfment and spore coat formation steps. Many B. subtilis sporulation genes, particularly those encoding small acid-soluble spore proteins and spore coat proteins, were found only in the family Bacillaceae, or even in a subset of Bacillus spp. Phylogenetic profiles of sporulation genes, compiled in this work, confirm the presence of a common sporulation gene core, but also illuminate the diversity of the sporulation processes within various lineages. These profiles should help further experimental studies of uncharacterized widespread sporulation genes, which would ultimately allow delineation of the minimal set(s) of sporulation-specific genes in Bacilli and Clostridia.

Phylogenomics of prokaryotic ribosomal proteins

Archaeal and bacterial ribosomes contain more than 50 proteins, including 34 that are universally conserved in the three domains of cellular life (bacteria, archaea, and eukaryotes). Despite the high sequence conservation, annotation of ribosomal (r-) protein genes is often difficult because of their short lengths and biased sequence composition. We developed an automated computational pipeline for identification of r-protein genes and applied it to 995 completely sequenced bacterial and 87 archaeal genomes available in the RefSeq database. The pipeline employs curated seed alignments of r-proteins to run position-specific scoring matrix (PSSM)-based BLAST searches against six-frame genome translations, mitigating possible gene annotation errors. As a result of this analysis, we performed a census of prokaryotic r-protein complements, enumerated missing and paralogous r-proteins, and analyzed the distributions of ribosomal protein genes among chromosomal partitions. Phyletic patterns of bacterial and archaeal r-protein genes were mapped to phylogenetic trees reconstructed from concatenated alignments of r-proteins to reveal the history of likely multiple independent gains and losses. These alignments, available for download, can be used as search profiles to improve genome annotation of r-proteins and for further comparative genomics studies.

Evolution of lysine biosynthesis in the phylum deinococcus-thermus

Thermus thermophilus biosynthesizes lysine through the α-aminoadipate (AAA) pathway: this observation was the first discovery of lysine biosynthesis through the AAA pathway in archaea and bacteria. Genes homologous to the T. thermophilus lysine biosynthetic genes are widely distributed in bacteria of the Deinococcus-Thermus phylum. Our phylogenetic analyses strongly suggest that a common ancestor of the Deinococcus-Thermus phylum had the ancestral genes for bacterial lysine biosynthesis through the AAA pathway. In addition, our findings suggest that the ancestor lacked genes for lysine biosynthesis through the diaminopimelate (DAP) pathway. Interestingly, Deinococcus proteolyticus does not have the genes for lysine biosynthesis through the AAA pathway but does have the genes for lysine biosynthesis through the DAP pathway. Phylogenetic analyses of D. proteolyticus lysine biosynthetic genes showed that the key gene cluster for the DAP pathway was transferred horizontally from a phylogenetically distant organism.

Genome Signature Difference between Deinococcus radiodurans and Thermus thermophilus

The extremely radioresistant bacteria of the genus Deinococcus and the extremely thermophilic bacteria of the genus Thermus belong to a common taxonomic group. Considering the distinct living environments of Deinococcus and Thermus, different genes would have been acquired through horizontal gene transfer after their divergence from a common ancestor. Their guanine-cytosine (GC) contents are similar; however, we hypothesized that their genomic signatures would be different. Our findings indicated that the genomes of Deinococcus radiodurans and Thermus thermophilus have different tetranucleotide frequencies. This analysis showed that the genome signature of D. radiodurans is most similar to that of Pseudomonas aeruginosa, whereas the genome signature of T. thermophilus is most similar to that of Thermanaerovibrio acidaminovorans. This difference in genome signatures may be related to the different evolutionary backgrounds of the 2 genera after their divergence from a common ancestor.