Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The Archaeome's Role in Colorectal Cancer: Unveiling the DPANN Group and Investigating Archaeal Functional Signatures

BACKGROUND AND AIMS

Gut microbial imbalances are linked to colorectal cancer (CRC), but archaea's role remains underexplored. Here, using previously published metagenomic data from different populations including Austria, Germany, Italy, Japan, China, and India, we performed bioinformatic and statistical analysis to identify archaeal taxonomic and functional signatures related to CRC.

METHODS

We analyzed published fecal metagenomic data from 390 subjects, comparing the archaeomes of CRC and healthy individuals. We conducted a biostatistical analysis to investigate the relationship between Candidatus Mancarchaeum acidiphilum (DPANN superphylum) and other archaeal species associated with CRC. Using the Prokka tool, we annotated the data focusing on archaeal genes, subsequently linking them to CRC and mapping them against UniprotKB and GO databases for specific archaeal gene functions.

RESULTS

Our analysis identified enrichment of methanogenic archaea in healthy subjects, with an exception for Methanobrevibacter smithii, which correlated with CRC. Notably, CRC showed a strong association with archaeal species, particularly Natrinema sp. J7-2, Ferroglobus placidus, and Candidatus Mancarchaeum acidiphilum. Furthermore, the DPANN archaeon exhibited a significant correlation with other CRC-associated archaea (p < 0.001). Functionally, we found a marked association between MvhB-type polyferredoxin and colorectal cancer. We also highlighted the association of archaeal proteins involved in the biosynthesis of leucine and the galactose metabolism process with the healthy phenotype.

CONCLUSIONS

The archaeomes of CRC patients show identifiable alterations, including a decline in methanogens and an increase in Halobacteria species. MvhB-type polyferredoxin, linked with CRC and species like Candidatus Mancarchaeum acidiphilum, Natrinema sp. J7-2, and Ferroglobus placidus emerge as potential archaeal biomarkers. Archaeal proteins may also offer gut protection, underscoring archaea's role in CRC dynamics.

Eukaryotes inherited inositol lipids from bacteria: implications for the models of eukaryogenesis

The merger of two very different microbes, an anaerobic archaeon and an aerobic bacterium, led to the birth of eukaryotic cells. Current models hypothesize that an archaeon engulfed bacteria through external protrusions that then fused together forming the membrane organelles of eukaryotic cells, including mitochondria. Images of cultivated Lokiarchaea sustain this concept, first proposed in the inside-out model which assumes that the membrane traffic system of archaea drove the merging with bacterial cells through membrane expansions containing inositol lipids, considered to have evolved first in archaea. This assumption has been evaluated here in detail. The data indicate that inositol lipids first emerged in bacteria, not in archaea. The implications of this finding for the models of eukaryogenesis are discussed.

Quinolinate Synthase: An Example of the Roles of the Second and Outer Coordination Spheres in Enzyme Catalysis

The activation energy barrier of biochemical reactions is normally lowered by an enzyme catalyst, which directly helps the weakening of the bond(s) to be broken. In many metalloenzymes, this is a first coordination sphere effect. Besides having a direct catalytic action, enzymes can fix their reactive groups and substrates so that they are optimally positioned and also modify the water activity in the system. They can either activate substrates prior to their reaction or bind preactivated substrates, thereby drastically reducing local entropic effects. The latter type is well represented by some bisubstrate reactions, where they have been defined as "entropic traps". These can be described as "second coordination sphere" processes, but enzymes can also control the reactivity beyond this point through local conformational changes belonging to an "outer coordinate sphere" that can be modulated by substrate binding. We have chosen the [4Fe-4S] cluster-dependent enzyme quinolinate synthase to illustrate each one of these processes. In addition, this very old metalloenzyme shows low in vitro substrate binding specificity, atypical reactivity that produces dead-end products, and a unique modulation of its active site volume.

Phylogenetic Signal, Congruence, and Uncertainty across Bacteria and Archaea

Reconstruction of the Tree of Life is a central goal in biology. Although numerous novel phyla of bacteria and archaea have recently been discovered, inconsistent phylogenetic relationships are routinely reported, and many inter-phylum and inter-domain evolutionary relationships remain unclear. Here, we benchmark different marker genes often used in constructing multidomain phylogenetic trees of bacteria and archaea and present a set of marker genes that perform best for multidomain trees constructed from concatenated alignments. We use recently-developed Tree Certainty metrics to assess the confidence of our results and to obviate the complications of traditional bootstrap-based metrics. Given the vastly disparate number of genomes available for different phyla of bacteria and archaea, we also assessed the impact of taxon sampling on multidomain tree construction. Our results demonstrate that biases between the representation of different taxonomic groups can dramatically impact the topology of resulting trees. Inspection of our highest-quality tree supports the division of most bacteria into Terrabacteria and Gracilicutes, with Thermatogota and Synergistota branching earlier from these superphyla. This tree also supports the inclusion of the Patescibacteria within the Terrabacteria as a sister group to the Chloroflexota instead of as a basal-branching lineage. For the Archaea, our tree supports three monophyletic lineages (DPANN, Euryarchaeota, and TACK/Asgard), although we note the basal placement of the DPANN may still represent an artifact caused by biased sequence composition. Our findings provide a robust and standardized framework for multidomain phylogenetic reconstruction that can be used to evaluate inter-phylum relationships and assess uncertainty in conflicting topologies of the Tree of Life.

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.

An evolutionary analysis identifies a conserved pentapeptide stretch containing the two essential lysine residues for rice L-myo-inositol 1-phosphate synthase catalytic activity

A molecular evolutionary analysis of a well conserved protein helps to determine the essential amino acids in the core catalytic region. Based on the chemical properties of amino acid residues, phylogenetic analysis of a total of 172 homologous sequences of a highly conserved enzyme, L-myo-inositol 1-phosphate synthase or MIPS from evolutionarily diverse organisms was performed. This study revealed the presence of six phylogenetically conserved blocks, out of which four embrace the catalytic core of the functional protein. Further, specific amino acid modifications targeting the lysine residues, known to be important for MIPS catalysis, were performed at the catalytic site of a MIPS from monocotyledonous model plant, Oryza sativa (OsMIPS1). Following this study, OsMIPS mutants with deletion or replacement of lysine residues in the conserved blocks were made. Based on the enzyme kinetics performed on the deletion/replacement mutants, phylogenetic and structural comparison with the already established crystal structures from non-plant sources, an evolutionarily conserved peptide stretch was identified at the active pocket which contains the two most important lysine residues essential for catalytic activity.

Adaptive Evolution of Thermotoga maritima Reveals Plasticity of the ABC Transporter Network

Thermotoga maritima is a hyperthermophilic anaerobe that utilizes a vast network of ABC transporters to efficiently metabolize a variety of carbon sources to produce hydrogen. For unknown reasons, this organism does not metabolize glucose as readily as it does glucose di- and polysaccharides. The leading hypothesis implicates the thermolability of glucose at the physiological temperatures at which T. maritima lives. After a 25-day laboratory evolution, phenotypes were observed with growth rates up to 1.4 times higher than and glucose utilization rates exceeding 50% those of the wild type. Genome resequencing revealed mutations in evolved cultures related to glucose-responsive ABC transporters. The native glucose ABC transporter, GluEFK, has more abundant transcripts either as a result of gene duplication-amplification or through mutations to the operator sequence regulating this operon. Conversely, BglEFGKL, a transporter of beta-glucosides, is substantially downregulated due to a nonsense mutation to the solute binding protein or due to a deletion of the upstream promoter. Analysis of the ABC2 uptake porter families for carbohydrate and peptide transport revealed that the solute binding protein, often among the transcripts detected at the highest levels, is predominantly downregulated in the evolved cultures, while the membrane-spanning domain and nucleotide binding components are less varied. Similar trends were observed in evolved strains grown on glycerol, a substrate that is not dependent on ABC transporters. Therefore, improved growth on glucose is achieved through mutations favoring GluEFK expression over BglEFGKL, and in lieu of carbon catabolite repression, the ABC transporter network is modulated to achieve improved growth fitness.

Horizontal gene transfer and genome evolution in Methanosarcina

Background

Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae.

Results

About 5 % of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome.

The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50 % of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46 % of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons.

Conclusions

While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0393-2) contains supplementary material, which is available to authorized users.

A role for [Fe4S4] clusters in tRNA recognition--a theoretical study

Over the past several years, structural studies have led to the unexpected discovery of iron–sulfur clusters in enzymes that are involved in DNA replication/repair and protein biosynthesis. Although these clusters are generally well-studied cofactors, their significance in the new contexts often remains elusive. One fascinating example is a tryptophanyl-tRNA synthetase from the thermophilic bacterium Thermotoga maritima, TmTrpRS, that has recently been structurally characterized. It represents an unprecedented connection among a primordial iron–sulfur cofactor, RNA and protein biosynthesis. Here, a possible role of the [Fe₄S₄] cluster in tRNA anticodon-loop recognition is investigated by means of density functional theory and comparison with the structure of a human tryptophanyl-tRNA synthetase/tRNA complex. It turns out that a cluster-coordinating cysteine residue, R224, and polar main chain atoms form a characteristic structural motif for recognizing a putative 5′ cytosine or 5′ 2-thiocytosine moiety in the anticodon loop of the tRNA molecule. This motif provides not only affinity but also specificity by creating a structural and energetical penalty for the binding of other bases, such as uracil.

Reconstructed ancestral Myo-inositol-3-phosphate synthases indicate that ancestors of the Thermococcales and Thermotoga species were more thermophilic than their descendants

The bacterial genomes of Thermotoga species show evidence of significant interdomain horizontal gene transfer from the Archaea. Members of this genus acquired many genes from the Thermococcales, which grow at higher temperatures than Thermotoga species. In order to study the functional history of an interdomain horizontally acquired gene we used ancestral sequence reconstruction to examine the thermal characteristics of reconstructed ancestral proteins of the Thermotoga lineage and its archaeal donors. Several ancestral sequence reconstruction methods were used to determine the possible sequences of the ancestral Thermotoga and Archaea myo-inositol-3-phosphate synthase (MIPS). These sequences were predicted to be more thermostable than the extant proteins using an established sequence composition method. We verified these computational predictions by measuring the activities and thermostabilities of purified proteins from the Thermotoga and the Thermococcales species, and eight ancestral reconstructed proteins. We found that the ancestral proteins from both the archaeal donor and the Thermotoga most recent common ancestor recipient were more thermostable than their descendants. We show that there is a correlation between the thermostability of MIPS protein and the optimal growth temperature (OGT) of its host, which suggests that the OGT of the ancestors of these species of Archaea and the Thermotoga grew at higher OGTs than their descendants.

Inositol lipids: from an archaeal origin to phosphatidylinositol 3,5-bisphosphate faults in human disease

The last couple of decades have seen an extraordinary transformation in our knowledge and understanding of the multifarious biological roles of inositol phospholipids. Herein, I briefly consider two topics. The first is the role that recently acquired biochemical and genomic information - especially from archaeons - has played in illuminating the possible evolutionary origins of the biological employment of inositol in lipids, and some questions that these studies raise about the 'classical' biosynthetic route to phosphatidylinositol. The second is the growing recognition of the importance in eukaryotic cells of phosphatidylinositol 3,5-bisphosphate. Phosphatidylinositol 3,5-bisphosphate only entered our phosphoinositide consciousness quite recently, but it is speedily gathering a plethora of roles in diverse cellular processes and diseases thereof. These include: control of endolysosomal vesicular trafficking and of the activity of ion channels and pumps in the endolysosomal compartment; control of constitutive and stimulated protein traffic to and from plasma membrane subdomains; control of the nutrient and stress-sensing target of rapamycin complex 1 pathway (TORC1); and regulation of key genes in some central metabolic pathways.

How hyperthermophiles adapt to change their lives: DNA exchange in extreme conditions

Transfer of DNA has been shown to be involved in genome evolution. In particular with respect to the adaptation of bacterial species to high temperatures, DNA transfer between the domains of bacteria and archaea seems to have played a major role. In addition, DNA exchange between similar species likely plays a role in repair of DNA via homologous recombination, a process that is crucial under DNA damaging conditions such as high temperatures. Several mechanisms for the transfer of DNA have been described in prokaryotes, emphasizing its general importance. However, until recently, not much was known about this process in prokaryotes growing in highly thermophilic environments. This review describes the different mechanisms of DNA transfer in hyperthermophiles, and how this may contribute to the survival and adaptation of hyperthermophilic archaea and bacteria to extreme environments.

Lateral gene transfer from the dead

In phylogenetic studies, the evolution of molecular sequences is assumed to have taken place along the phylogeny traced by the ancestors of extant species. In the presence of lateral gene transfer, however, this may not be the case, because the species lineage from which a gene was transferred may have gone extinct or not have been sampled. Because it is not feasible to specify or reconstruct the complete phylogeny of all species, we must describe the evolution of genes outside the represented phylogeny by modeling the speciation dynamics that gave rise to the complete phylogeny. We demonstrate that if the number of sampled species is small compared with the total number of existing species, the overwhelming majority of gene transfers involve speciation to and evolution along extinct or unsampled lineages. We show that the evolution of genes along extinct or unsampled lineages can to good approximation be treated as those of independently evolving lineages described by a few global parameters. Using this result, we derive an algorithm to calculate the probability of a gene tree and recover the maximum-likelihood reconciliation given the phylogeny of the sampled species. Examining 473 near-universal gene families from 36 cyanobacteria, we find that nearly a third of transfer events (28%) appear to have topological signatures of evolution along extinct species, but only approximately 6% of transfers trace their ancestry to before the common ancestor of the sampled cyanobacteria. [Gene tree reconciliation; lateral gene transfer; macroevolution; phylogeny.]

COMPARATIVE MODELING OF COEVOLUTION IN COMMUNITIES OF UNICELLULAR ORGANISMS: ADAPTABILITY AND BIODIVERSITY

We propose an original program "Evolutionary constructor" that is capable of computationally efficient modeling of both population-genetic and ecological problems, combining these directions in one model of required detail level. We also present results of comparative modeling of stability, adaptability and biodiversity dynamics in populations of unicellular haploid organisms which form symbiotic ecosystems. The advantages and disadvantages of two evolutionary strategies of biota formation — a few generalists' taxa-based biota formation and biodiversity-based biota formation — are discussed.

Genome sequence of Thermotoga sp. strain RQ2, a hyperthermophilic bacterium isolated from a geothermally heated region of the seafloor near Ribeira Quente, the Azores

Thermotoga sp. strain RQ2 is probably a strain of Thermotoga maritima. Its complete genome sequence allows for an examination of the extent and consequences of gene flow within Thermotoga species and strains. Thermotoga sp. RQ2 differs from T. maritima in its genes involved in myo-inositol metabolism. Its genome also encodes an apparent fructose phosphotransferase system (PTS) sugar transporter. This operon is also found in Thermotoga naphthophila strain RKU-10 but no other Thermotogales. These are the first reported PTS transporters in the Thermotogales.

Evolution of the biosynthesis of di-myo-inositol phosphate, a marker of adaptation to hot marine environments

The synthesis of di-myo-inositol phosphate (DIP), a common compatible solute in hyperthermophiles, involves the consecutive actions of inositol-1-phosphate cytidylyltransferase (IPCT) and di-myo-inositol phosphate phosphate synthase (DIPPS). In most cases, both activities are present in a single gene product, but separate genes are also found in a few organisms. Genes for IPCT and DIPPS were found in the genomes of 33 organisms, all with thermophilic/hyperthermophilic lifestyles. Phylogeny of IPCT/DIPPS revealed an incongruent topology with 16S RNA phylogeny, thus suggesting horizontal gene transfer. The phylogenetic tree of the DIPPS domain was rooted by using phosphatidylinositol phosphate synthase sequences as out-group. The root locates at the separation of genomes with fused and split genes. We propose that the gene encoding DIPPS was recruited from the biosynthesis of phosphatidylinositol. The last DIP-synthesizing ancestor harboured separated genes for IPCT and DIPPS and this architecture was maintained in a crenarchaeal lineage, and transferred by horizontal gene transfer to hyperthermophilic marine Thermotoga species. It is plausible that the driving force for the assembly of those two genes in the early ancestor is related to the acquired advantage of DIP producers to cope with high temperature. This work corroborates the view that Archaea were the first hyperthermophilic organisms.

Comparative genomics analysis of completely sequenced microbial genomes reveals the ubiquity of N-linked glycosylation in prokaryotes

Glycosylation of proteins in prokaryotes has been known for the last few decades. Glycan structures and/or the glycosylation pathways have been experimentally characterized in only a small number of prokaryotes. Even this has become possible only during the last decade or so, primarily due to technological and methodological developments. Glycosylated proteins are diverse in their function and localization. Glycosylation has been shown to be associated with a wide range of biological phenomena. Characterization of the various types of glycans and the glycosylation machinery is critical to understand such processes. Such studies can help in the identification of novel targets for designing drugs, diagnostics, and engineering of therapeutic proteins. In view of this, the experimentally characterized pgl system of Campylobacter jejuni, responsible for N-linked glycosylation, has been used in this study to identify glycosylation loci in 865 prokaryotes whose genomes have been completely sequenced. Results from the present study show that only a small number of organisms have homologs for all the pgl enzymes and a few others have homologs for none of the pgl enzymes. Most of the organisms have homologs for only a subset of the pgl enzymes. There is no specific pattern for the presence or absence of pgl homologs vis-à-vis the 16S rRNA sequence-based phylogenetic tree. This may be due to differences in the glycan structures, high sequence divergence, horizontal gene transfer or non-orthologous gene displacement. Overall, the presence of homologs for pgl enzymes in a large number of organisms irrespective of their habitat, pathogenicity, energy generation mechanism, etc., hints towards the ubiquity of N-linked glycosylation in prokaryotes.

A structure-centric view of protein evolution, design, and adaptation

Proteins, by virtue of their central role in most biological processes, represent one of the key subjects of the study of molecular evolution. Inherent in the indispensability of proteins for living cells is the fact that a given protein can adopt a specific three-dimensional shape that is specified solely by the protein's sequence of amino acids. Over the past several decades, structural biologists have demonstrated that the array of structures that proteins may adopt is quite astounding, and this has lead to a strong interest in understanding how protein structures change and evolve over time. In this review we consider a large body of recent work that attempts to illuminate this structure-centric picture of protein evolution. Much of this work has focused on the question of how completely new protein structures (i.e., new folds or topologies) are discovered by protein sequences as they evolve. Pursuant to this question of structural innovation has been a desire to describe and understand the observation that certain types of protein structures are far more abundant than others and how this uneven distribution of proteins implicates on the process through which new shapes are discovered. We consider a number of theoretical models that have been successful at explaining this heterogeneity in protein populations and discuss the increasing amount of evidence that indicates that the process of structural evolution involves the divergence of protein sequences and structures from one another. We also consider the topic of protein designability, which concerns itself with understanding how a protein's structure influences the number of sequences that can fold successfully into that structure. Understanding and quantifying the relationship between the physical feature of a structure and its designability has been a long-standing goal of the study of protein structure and evolution, and we discuss a number of recent advances that have yielded a promising answer to this question. Finally, we review the relatively new field of protein structural phylogeny, an area of study in which information about the distribution of protein structures among different organisms is used to reconstruct the evolutionary relationships between them. Taken together, the work that we review presents an increasingly coherent picture of how these unique polymers have evolved over the course of life on Earth.

Deep phylogeny, ancestral groups and the four ages of life

Organismal phylogeny depends on cell division, stasis, mutational divergence, cell mergers (by sex or symbiogenesis), lateral gene transfer and death. The tree of life is a useful metaphor for organismal genealogical history provided we recognize that branches sometimes fuse. Hennigian cladistics emphasizes only lineage splitting, ignoring most other major phylogenetic processes. Though methodologically useful it has been conceptually confusing and harmed taxonomy, especially in mistakenly opposing ancestral (paraphyletic) taxa. The history of life involved about 10 really major innovations in cell structure. In membrane topology, there were five successive kinds of cell: (i) negibacteria, with two bounding membranes, (ii) unibacteria, with one bounding and no internal membranes, (iii) eukaryotes with endomembranes and mitochondria, (iv) plants with chloroplasts and (v) finally, chromists with plastids inside the rough endoplasmic reticulum. Membrane chemistry divides negibacteria into the more advanced Glycobacteria (e.g. Cyanobacteria and Proteobacteria) with outer membrane lipolysaccharide and primitive Eobacteria without lipopolysaccharide (deserving intenser study). It also divides unibacteria into posibacteria, ancestors of eukaryotes, and archaebacteria-the sisters (not ancestors) of eukaryotes and the youngest bacterial phylum. Anaerobic eobacteria, oxygenic cyanobacteria, desiccation-resistant posibacteria and finally neomura (eukaryotes plus archaebacteria) successively transformed Earth. Accidents and organizational constraints are as important as adaptiveness in body plan evolution.

On the need for widespread horizontal gene transfers under genome size constraint

Background

While eukaryotes primarily evolve by duplication-divergence expansion (and reduction) of their own gene repertoire with only rare horizontal gene transfers, prokaryotes appear to evolve under both gene duplications and widespread horizontal gene transfers over long evolutionary time scales. But, the evolutionary origin of this striking difference in the importance of horizontal gene transfers remains by and large a mystery.

Hypothesis

We propose that the abundance of horizontal gene transfers in free-living prokaryotes is a simple but necessary consequence of two opposite effects: i) their apparent genome size constraint compared to typical eukaryote genomes and ii) their underlying genome expansion dynamics through gene duplication-divergence evolution, as demonstrated by the presence of many tandem and block repeated genes. In principle, this combination of genome size constraint and underlying duplication expansion should lead to a coalescent-like process with extensive turnover of functional genes. This would, however, imply the unlikely, systematic reinvention of functions from discarded genes within independent phylogenetic lineages. Instead, we propose that the long-term evolutionary adaptation of free-living prokaryotes must have resulted in the emergence of efficient non-phylogenetic pathways to circumvent gene loss.

Implications

This need for widespread horizontal gene transfers due to genome size constraint implies, in particular, that prokaryotes must remain under strong selection pressure in order to maintain the long-term evolutionary adaptation of their "mutualized" gene pool, beyond the inevitable turnover of individual prokaryote species. By contrast, the absence of genome size constraint for typical eukaryotes has presumably relaxed their need for widespread horizontal gene transfers and strong selection pressure. Yet, the resulting loss of genetic functions, due to weak selection pressure and inefficient gene recovery mechanisms, must have ultimately favored the emergence of more complex life styles and ecological integration of many eukaryotes.

Reviewers

This article was reviewed by Pierre Pontarotti, Eugene V Koonin and Sergei Maslov.

An ancestral glutamine-dependent NAD(+) synthetase revealed by poor kinetic synergism

NAD(+) synthetase catalyzes the formation of NAD(+) from ATP, nicotinic acid adenine dinucleotide and ammonia. Glutamine-dependent NAD(+) synthetase obtains ammonia through the hydrolysis of glutamine to glutamate, which takes place in the glutaminase domain. The ammonia is subsequently transported to the synthetase domain through an interdomain ammonia tunnel. NAD(+) synthetase from the thermophilic bacteria Thermotoga maritima was cloned and expressed. Steady-state kinetics and stoichiometric analysis of product formation revealed an enzyme that is significantly inefficient in the synchronization of the two active sites resulting in wasteful hydrolysis of glutamine and that is not specific for glutamine over ammonia. Phylogenetic analysis of glutamine-dependent NAD(+) synthetases identifies three main groups remotely related. The T. maritima NAD(+) synthetase's group is proposed to represent the ancestral group based on the phylogenetic analysis and on the kinetic characterizations. The phylogenetic results nicely correlate also with the degree of catalytic efficiency measured for M. tuberculosis, S. cerevisiae and T. maritima NAD(+) synthetases. Furthermore, the data here reported in combination with structural data available for glutamine-dependent NAD(+) synthetase lays the foundation for further investigation on the mechanism of active site coupling in these enzymes.

A novel limb in the osmoregulatory network of Methanosarcina mazei Gö1: N(epsilon)-acetyl-beta-lysine can be substituted by glutamate and alanine

N(epsilon)-acetyl-beta-lysine is a unique compatible solute found in methanogenic archaea grown at high salinities. Deletion of the genes that encode the lysine-2,3-aminomutase (ablA) and the beta-lysine acetyltransferase (ablB) abolished the production of N(epsilon)-acetyl-beta-lysine in Methanosarcina mazei Gö1. The mutant grew well at low and intermediate salinities. Interestingly, growth at high salt (800 mM NaCl) was only slowed down but not impaired demonstrating that in M. mazei Gö1 N(epsilon)-acetyl-beta-lysine is not essential for growth at high salinities. Nuclear magnetic resonance (NMR) analysis revealed an increased glutamate pool in the mutant. In addition to alpha-glutamate, a novel solute, alanine, was produced. The intracellular alanine concentration was as high as 0.36 +/- 0.05 micromol (mg protein)-1 representing up to 18% of the total solute pool at 800 mM NaCl. The cellular alanine concentration increased with the salinity of the medium and decreased in the presence of glycine betaine in the medium, indicating that alanine is used as compatible solute by M. mazei Gö1.

Interdomain transfers of sugar transporters overcome barriers to gene expression

Horizontal gene transfer (HGT) is nature's mechanism for sharing evolved physiological traits among the members of microbial communities. The extent to which such transfers can be successful is best illustrated by the fact that Archaea-derived genes are found in many bacterial genomes, particularly those in the hyperthermophile Thermotoga maritima. The success of these intergenomic transfers depends upon the successful transcription of the newly acquired archaeal genes using a bacterial transcription machinery that does not recognize archaeal transcriptional signals. To examine how nature solves this problem, we looked to the T. maritima genome for examples of interdomain transfers. Here we lay the groundwork to examine this problem by more clearly delineating the phylogenetic history of Archaea-derived transporter genes in this genome. We find that five of these polysaccharide transporters were derived from the Archaea and one came from the Archaea after that lineage inherited it from the Bacteria. These data can be used for more detailed examinations of the recombinations that allowed these transporters to be expressed in a bacterial host. This work will guide examinations of the genome sequences from other members of the Thermotogales, which will become available.

Diversity and transcription of proteases involved in the maturation of hydrogenases in Nostoc punctiforme ATCC 29133 and Nostoc sp. strain PCC 7120

BACKGROUND

The last step in the maturation process of the large subunit of [NiFe]-hydrogenases is a proteolytic cleavage of the C-terminal by a hydrogenase specific protease. Contrary to other accessory proteins these hydrogenase proteases are believed to be specific whereby one type of hydrogenases specific protease only cleaves one type of hydrogenase. In cyanobacteria this is achieved by the gene product of either hupW or hoxW, specific for the uptake or the bidirectional hydrogenase respectively. The filamentous cyanobacteria Nostoc punctiforme ATCC 29133 and Nostoc sp strain PCC 7120 may contain a single uptake hydrogenase or both an uptake and a bidirectional hydrogenase respectively.

RESULTS

In order to examine these proteases in cyanobacteria, transcriptional analyses were performed of hupW in Nostoc punctiforme ATCC 29133 and hupW and hoxW in Nostoc sp. strain PCC 7120. These studies revealed numerous transcriptional start points together with putative binding sites for NtcA (hupW) and LexA (hoxW). In order to investigate the diversity and specificity among hydrogeanse specific proteases we constructed a phylogenetic tree which revealed several subgroups that showed a striking resemblance to the subgroups previously described for [NiFe]-hydrogenases. Additionally the proteases specificity was also addressed by amino acid sequence analysis and protein-protein docking experiments with 3D-models derived from bioinformatic studies. These studies revealed a so called "HOXBOX"; an amino acid sequence specific for protease of Hox-type which might be involved in docking with the large subunit of the hydrogenase.

CONCLUSION

Our findings suggest that the hydrogenase specific proteases are under similar regulatory control as the hydrogenases they cleave. The result from the phylogenetic study also indicates that the hydrogenase and the protease have co-evolved since ancient time and suggests that at least one major horizontal gene transfer has occurred. This co-evolution could be the result of a close interaction between the protease and the large subunit of the [NiFe]-hydrogenases, a theory supported by protein-protein docking experiments performed with 3D-models. Finally we present data that may explain the specificity seen among hydrogenase specific proteases, the so called "HOXBOX"; an amino acid sequence specific for proteases of Hox-type. This opens the door for more detailed studies of the specificity found among hydrogenase specific proteases and the structural properties behind it.

Population genomics and the bacterial species concept

In recent years, the importance of horizontal gene transfer (HGT) in bacterial evolution has been elevated to such a degree that many bacteriologists now question the very existence of bacterial species. If gene transfer is as rampant as comparative genomic studies have suggested, how could bacterial species survive such genomic fluidity? And yet, most bacteriologists recognize, and name, as species, clusters of bacterial isolates that share complex phenotypic properties. The Core Genome Hypothesis (CGH) has been proposed to explain this apparent paradox of fluid bacterial genomes associated with stable phenotypic clusters. It posits that there is a core of genes responsible for maintaining the species-specific phenotypic clusters observed throughout bacterial diversity and argues that, even in the face of substantial genomic fluidity, bacterial species can be rationally identified and named.

The archaebacterial origin of eukaryotes

The origin of the eukaryotic genetic apparatus is thought to be central to understanding the evolution of the eukaryotic cell. Disagreement about the source of the relevant genes has spawned competing hypotheses for the origins of the eukaryote nuclear lineage. The iconic rooted 3-domains tree of life shows eukaryotes and archaebacteria as separate groups that share a common ancestor to the exclusion of eubacteria. By contrast, the eocyte hypothesis has eukaryotes originating within the archaebacteria and sharing a common ancestor with a particular group called the Crenarchaeota or eocytes. Here, we have investigated the relative support for each hypothesis from analysis of 53 genes spanning the 3 domains, including essential components of the eukaryotic nucleic acid replication, transcription, and translation apparatus. As an important component of our analysis, we investigated the fit between model and data with respect to composition. Compositional heterogeneity is a pervasive problem for reconstruction of ancient relationships, which, if ignored, can produce an incorrect tree with strong support. To mitigate its effects, we used phylogenetic models that allow for changing nucleotide or amino acid compositions over the tree and data. Our analyses favor a topology that supports the eocyte hypothesis rather than archaebacterial monophyly and the 3-domains tree of life.

A unique combination of genetic systems for the synthesis of trehalose in Rubrobacter xylanophilus: properties of a rare actinobacterial TreT

Trehalose is the primary organic solute in Rubrobacter xylanophilus under all conditions tested, including those for optimal growth. We detected genes of four different pathways for trehalose synthesis in the genome of this organism, namely, the trehalose-6-phosphate synthase (Tps)/trehalose-6-phosphate phosphatase (Tpp), TreS, TreY/TreZ, and TreT pathways. Moreover, R. xylanophilus is the only known member of the phylum Actinobacteria to harbor TreT. The Tps sequence is typically bacterial, but the Tpp sequence is closely related to eukaryotic counterparts. Both the Tps/Tpp and the TreT pathways were active in vivo, while the TreS and the TreY/TreZ pathways were not active under the growth conditions tested and appear not to contribute to the levels of trehalose observed. The genes from the active pathways were functionally expressed in Escherichia coli, and Tps was found to be highly specific for GDP-glucose, a rare feature among these enzymes. The trehalose-6-phosphate formed was specifically dephosphorylated to trehalose by Tpp. The recombinant TreT synthesized trehalose from different nucleoside diphosphate-glucose donors and glucose, but the activity in R. xylanophilus cell extracts was specific for ADP-glucose. The TreT could also catalyze trehalose hydrolysis in the presence of ADP, but with a very high K(m). Here, we functionally characterize two systems for the synthesis of trehalose in R. xylanophilus, a representative of an ancient lineage of the actinobacteria, and discuss a possible scenario for the exceptional occurrence of treT in this extremophilic bacterium.

Risks from GMOs due to horizontal gene transfer

Horizontal gene transfer (HGT) is the stable transfer of genetic material from one organism to another without reproduction or human intervention. Transfer occurs by the passage of donor genetic material across cellular boundaries, followed by heritable incorporation to the genome of the recipient organism. In addition to conjugation, transformation and transduction, other diverse mechanisms of DNA and RNA uptake occur in nature. The genome of almost every organism reveals the footprint of many ancient HGT events. Most commonly, HGT involves the transmission of genes on viruses or mobile genetic elements. HGT first became an issue of public concern in the 1970s through the natural spread of antibiotic resistance genes amongst pathogenic bacteria, and more recently with commercial production of genetically modified (GM) crops. However, the frequency of HGT from plants to other eukaryotes or prokaryotes is extremely low. The frequency of HGT to viruses is potentially greater, but is restricted by stringent selection pressures. In most cases the occurrence of HGT from GM crops to other organisms is expected to be lower than background rates. Therefore, HGT from GM plants poses negligible risks to human health or the environment.

Structure-function studies of glutamate synthases: a class of self-regulated iron-sulfur flavoenzymes essential for nitrogen assimilation

Glutamate synthases play with glutamine synthetase an essential role in nitrogen assimilation processes in microorganisms, plants, and lower animals by catalyzing the net synthesis of one molecule of L-glutamate from L-glutamine and 2-oxoglutarate. They exhibit a modular architecture with a common subunit or region, which is responsible for the L-glutamine-dependent glutamate synthesis from 2-oxoglutarate. Here, a PurF- (Type II- or Ntn-) type amidotransferase domain is coupled to the synthase domain, a (beta/alpha)8 barrel containing FMN and one [3Fe-4S]0,+1 cluster, through a approximately 30 angstroms-long intramolecular tunnel for the transfer of ammonia between the sites. In bacterial and eukaryotic GltS, reducing equivalents are provided by reduced pyridine nucleotides thanks to the stable association with a second subunit or region, which acts as a FAD-dependent NAD(P)H oxidoreductase and is responsible for the formation of the two low potential [4Fe-4S]+1,+2 clusters of the enzyme. In photosynthetic cells, reduced ferredoxin is the physiological reductant. This review focus on the mechanism of cross-activation of the synthase and glutaminase reactions in response to the bound substrates and the redox state of the enzyme cofactors, as well as on recent information on the structure of the alphabeta protomer of the NADPH-dependent enzyme, which sheds light on the intramolecular electron transfer pathway between the flavin cofactors.

Structure and content of the Entamoeba histolytica genome

The intestinal parasite Entamoeba histolytica is one of the first protists for which a draft genome sequence has been published. Although the genome is still incomplete, it is unlikely that many genes are missing from the list of those already identified. In this chapter we summarise the features of the genome as they are currently understood and provide previously unpublished analyses of many of the genes.

Evolution of mal ABC transporter operons in the Thermococcales and Thermotogales

BACKGROUND

The mal genes that encode maltose transporters have undergone extensive lateral transfer among ancestors of the archaea Thermococcus litoralis and Pyrococcus furiosus. Bacterial hyperthermophiles of the order Thermotogales live among these archaea and so may have shared in these transfers. The genome sequence of Thermotoga maritima bears evidence of extensive acquisition of archaeal genes, so its ancestors clearly had the capacity to do so. We examined deep phylogenetic relationships among the mal genes of these hyperthermophiles and their close relatives to look for evidence of shared ancestry.

RESULTS

We demonstrate that the two maltose ATP binding cassette (ABC) transporter operons now found in Tc. litoralis and P. furiosus (termed mal and mdx genes, respectively) are not closely related to one another. The Tc. litoralis and P. furiosus mal genes are most closely related to bacterial mal genes while their respective mdx genes are archaeal. The genes of the two mal operons in Tt. maritima are not related to genes in either of these archaeal operons. They are highly similar to one another and belong to a phylogenetic lineage that includes mal genes from the enteric bacteria. A unique domain of the enteric MalF membrane spanning proteins found also in these Thermotogales MalF homologs supports their relatively close relationship with these enteric proteins. Analyses of genome sequence data from other Thermotogales species, Fervidobacterium nodosum, Thermosipho melanesiensis, Thermotoga petrophila, Thermotoga lettingae, and Thermotoga neapolitana, revealed a third apparent mal operon, absent from the published genome sequence of Tt. maritima strain MSB8. This third operon, mal3, is more closely related to the Thermococcales' bacteria-derived mal genes than are mal1 and mal2. F. nodosum, Ts. melanesiensis, and Tt. lettingae have only one of the mal1-mal2 paralogs. The mal2 operon from an unknown species of Thermotoga appears to have been horizontally acquired by a Thermotoga species that had only mal1.

CONCLUSION

These data demonstrate that the Tc. litoralis and P. furiosus mdx maltodextrin transporter operons arose in the Archaea while their mal maltose transporter operons arose in a bacterial lineage, but not the same lineage as the two maltose transporter operons found in the published Tt. maritima genome sequence. These Tt. maritima maltose transporters are phylogenetically and structurally similar to those found in enteric bacteria and the mal2 operon was horizontally transferred within the Thermotoga lineage. Other Thermotogales species have a third mal operon that is more closely related to the bacterial Thermococcales mal operons, but the data do not support a recent horizontal sharing of that operon between these groups.

DarkHorse: a method for genome-wide prediction of horizontal gene transfer

DarkHorse is a new approach to rapid, genome-wide identification and ranking of horizontal transfer candidate proteins.

A new approach to rapid, genome-wide identification and ranking of horizontal transfer candidate proteins is presented. The method is quantitative, reproducible, and computationally undemanding. It can be combined with genomic signature and/or phylogenetic tree-building procedures to improve accuracy and efficiency. The method is also useful for retrospective assessments of horizontal transfer prediction reliability, recognizing orthologous sequences that may have been previously overlooked or unavailable. These features are demonstrated in bacterial, archaeal, and eukaryotic examples.

Streptococcus mutans clonal variation revealed by multilocus sequence typing

Streptococcus mutans is the major pathogen of dental caries, a biofilm-dependent infectious disease, and occasionally causes infective endocarditis. S. mutans strains have been classified into four serotypes (c, e, f, and k). However, little is known about the S. mutans population, including the clonal relationships among strains of S. mutans, in relation to the particular clones that cause systemic diseases. To address this issue, we have developed a multilocus sequence typing (MLST) scheme for S. mutans. Eight housekeeping gene fragments were sequenced from each of 102 S. mutans isolates collected from the four serotypes in Japan and Finland. Between 14 and 23 alleles per locus were identified, allowing us theoretically to distinguish more than 1.2 x 10(10) sequence types. We identified 92 sequence types in these 102 isolates, indicating that S. mutans contains a diverse population. Whereas serotype c strains were widely distributed in the dendrogram, serotype e, f, and k strains were differentiated into clonal complexes. Therefore, we conclude that the ancestral strain of S. mutans was serotype c. No geographic specificity was identified. However, the distribution of the collagen-binding protein gene (cnm) and direct evidence of mother-to-child transmission were clearly evident. In conclusion, the superior discriminatory capacity of this MLST scheme for S. mutans may have important practical implications.

Effects of the nonionic surfactant tween 80 on microbial reductive dechlorination of chlorinated ethenes

Recent field studies have indicated synergistic effects of coupling microbial reductive dechlorination with physicochemical remediation (e.g., surfactant flushing) of dense nonaqueous phase liquid (DNAPL) source zones. This study explored chlorinated ethene (e.g., tetrachloroethene [PCE]) dechlorination in the presence of 50-5000 mg/L Tween 80, a nonionic surfactant employed in source zone remediation. Tween 80 did not inhibit dechlorination by four pure PCE-to-cis-1,2-dichloroethene (cis-DCE) or PCE-to-trichloroethene (TCE) dechlorinating cultures. In contrast, cis-DCE-dechlorinating Dehalococcoides isolates (strain BAV1 and strain FL2) failed to dechlorinate in the presence of Tween 80. Bio-Dechlor INOCULUM (BDI), a PCE-to-ethene dechlorinating consortium, produced cis-DCE in the presence of Tween 80, further suggesting that Tween 80 inhibits dechlorination by Dehalococcoides organisms. Quantitative real-time PCR analysis applied to BDI revealed that the number of Dehalococcoides cells decayed exponentially (R(2) = 0.85) according to the Chick-Watson disinfection model (pseudo first-order decay rate of 0.13+/-0.02 day(-1)) from an initial value of 6.6 +/-1.5 x 10(8) to 1.3+/-0.8 x 10(5) per mL of culture after 58 days of exposure to 250 mg/L Tween 80. Although Tween 80 exposure prevented ethene formation and reduced Dehalococcoides cell numbers, Dehalococcoides organisms remained viable, and dechlorination activity pist cis-DCE was recovered following the removal of Tween 80. These findings suggest that sequential Tween 80 flushing followed by microbial reductive dechlorination is a promising strategy for remediation of chlorinated ethene-impacted source zones.

Microbial biochemistry, physiology, and biotechnology of hyperthermophilic Thermotoga species

High-throughput sequencing of microbial genomes has allowed the application of functional genomics methods to species lacking well-developed genetic systems. For the model hyperthermophile Thermotoga maritima, microarrays have been used in comparative genomic hybridization studies to investigate diversity among Thermotoga species. Transcriptional data have assisted in prediction of pathways for carbohydrate utilization, iron-sulfur cluster synthesis and repair, expolysaccharide formation, and quorum sensing. Structural genomics efforts aimed at the T. maritima proteome have yielded hundreds of high-resolution datasets and predicted functions for uncharacterized proteins. The information gained from genomics studies will be particularly useful for developing new biotechnology applications for T. maritima enzymes.