Updated clusters of orthologous genes for Archaea: a complex ancestor of the Archaea and the byways of horizontal gene transfer

A Novel Transcriptional Regulator Related to Thiamine Phosphate Synthase Controls Thiamine Metabolism Genes in Archaea

Thiamine (vitamin B₁) is a precursor of thiamine pyrophosphate (TPP), an essential coenzyme in the central metabolism of all living organisms. Bacterial thiamine biosynthesis and salvage genes are controlled at the RNA level by TPP-responsive riboswitches. In Archaea, TPP riboswitches are restricted to the Thermoplasmatales order. Mechanisms of transcriptional control of thiamine genes in other archaeal lineages remain unknown. Using the comparative genomics approach, we identified a novel family of transcriptional regulators (named ThiR) controlling thiamine biosynthesis and transport genes in diverse lineages in the Crenarchaeota phylum as well as in the Halobacteria and Thermococci classes of the Euryarchaeota ThiR regulators are composed of an N-terminal DNA-binding domain and a C-terminal ligand-binding domain, which is similar to the archaeal thiamine phosphate synthase ThiN. By using comparative genomics, we predicted ThiR-binding DNA motifs and reconstructed ThiR regulons in 67 genomes representing all above-mentioned lineages. The predicted ThiR-binding motifs are characterized by palindromic symmetry with several distinct lineage-specific consensus sequences. In addition to thiamine biosynthesis genes, the reconstructed ThiR regulons include various transporters for thiamine and its precursors. Bioinformatics predictions were experimentally validated by in vitro DNA-binding assays with the recombinant ThiR protein from the hyperthermophilic archaeon Metallosphaera yellowstonensis MK1. Thiamine phosphate and, to some extent, TPP and hydroxyethylthiazole phosphate were required for the binding of ThiR to its DNA targets, suggesting that ThiR is derepressed by limitation of thiamine phosphates. The thiamine phosphate-binding residues previously identified in ThiN are highly conserved in ThiR regulators, suggesting a conserved mechanism for effector recognition.

IMPORTANCE

Thiamine pyrophosphate is a cofactor for many essential enzymes for glucose and energy metabolism. Thiamine or vitamin B₁ biosynthesis and its transcriptional regulation in Archaea are poorly understood. We applied the comparative genomics approach to identify a novel family of regulators for the transcriptional control of thiamine metabolism genes in Archaea and reconstructed the respective regulons. The predicted ThiR regulons in archaeal genomes control the majority of thiamine biosynthesis genes. The reconstructed regulon content suggests that numerous uptake transporters for thiamine and/or its precursors are encoded in archaeal genomes. The ThiR regulon was experimentally validated by DNA-binding assays with Metallosphaera spp. These discoveries contribute to our understanding of metabolic and regulatory networks involved in vitamin homeostasis in diverse lineages of Archaea.

Biology of archaea from a novel family Cuniculiplasmataceae (Thermoplasmata) ubiquitous in hyperacidic environments

The order Thermoplasmatales (Euryarchaeota) is represented by the most acidophilic organisms known so far that are poorly amenable to cultivation. Earlier culture-independent studies in Iron Mountain (California) pointed at an abundant archaeal group, dubbed ‘G-plasma’. We examined the genomes and physiology of two cultured representatives of a Family Cuniculiplasmataceae, recently isolated from acidic (pH 1–1.5) sites in Spain and UK that are 16S rRNA gene sequence-identical with ‘G-plasma’. Organisms had largest genomes among Thermoplasmatales (1.87–1.94 Mbp), that shared 98.7–98.8% average nucleotide identities between themselves and ‘G-plasma’ and exhibited a high genome conservation even within their genomic islands, despite their remote geographical localisations. Facultatively anaerobic heterotrophs, they possess an ancestral form of A-type terminal oxygen reductase from a distinct parental clade. The lack of complete pathways for biosynthesis of histidine, valine, leucine, isoleucine, lysine and proline pre-determines the reliance on external sources of amino acids and hence the lifestyle of these organisms as scavengers of proteinaceous compounds from surrounding microbial community members. In contrast to earlier metagenomics-based assumptions, isolates were S-layer-deficient, non-motile, non-methylotrophic and devoid of iron-oxidation despite the abundance of methylotrophy substrates and ferrous iron in situ, which underlines the essentiality of experimental validation of bioinformatic predictions.

Inevitability of Genetic Parasites

Almost all cellular life forms are hosts to diverse genetic parasites with various levels of autonomy including plasmids, transposons and viruses. Theoretical modeling of the evolution of primordial replicators indicates that parasites (cheaters) necessarily evolve in such systems and can be kept at bay primarily via compartmentalization. Given the (near) ubiquity, abundance and diversity of genetic parasites, the question becomes pertinent: are such parasites intrinsic to life? At least in prokaryotes, the persistence of parasites is linked to the rate of horizontal gene transfer (HGT). We mathematically derive the threshold value of the minimal transfer rate required for selfish element persistence, depending on the element duplication and loss rates as well as the cost to the host. Estimation of the characteristic gene duplication, loss and transfer rates for transposons, plasmids and virus-related elements in multiple groups of diverse bacteria and archaea indicates that most of these rates are compatible with the long term persistence of parasites. Notably, a small but non-zero rate of HGT is also required for the persistence of non-parasitic genes. We hypothesize that cells cannot tune their horizontal transfer rates to be below the threshold required for parasite persistence without experiencing highly detrimental side-effects. As a lower boundary to the minimum DNA transfer rate that a cell can withstand, we consider the process of genome degradation and mutational meltdown of populations through Muller's ratchet. A numerical assessment of this hypothesis suggests that microbial populations cannot purge parasites while escaping Muller's ratchet. Thus, genetic parasites appear to be virtually inevitable in cellular organisms.

Genomic composition and dynamics among Methanomicrobiales predict adaptation to contrasting environments

Members of the order Methanomicrobiales are abundant, and sometimes dominant, hydrogenotrophic (H₂-CO₂ utilizing) methanoarchaea in a broad range of anoxic habitats. Despite their key roles in greenhouse gas emissions and waste conversion to methane, little is known about the physiological and genomic bases for their widespread distribution and abundance. In this study, we compared the genomes of nine diverse Methanomicrobiales strains, examined their pangenomes, reconstructed gene flow and identified genes putatively mediating their success across different habitats. Most strains slowly increased gene content whereas one, Methanocorpusculum labreanum, evidenced genome downsizing. Peat-dwelling Methanomicrobiales showed adaptations centered on improved transport of scarce inorganic nutrients and likely use H⁺ rather than Na⁺ transmembrane chemiosmotic gradients during energy conservation. In contrast, other Methanomicrobiales show the potential to concurrently use Na⁺ and H⁺ chemiosmotic gradients. Analyses also revealed that the Methanomicrobiales lack a canonical electron bifurcation system (MvhABGD) known to produce low potential electrons in other orders of hydrogenotrophic methanogens. Additional putative differences in anabolic metabolism suggest that the dynamics of interspecies electron transfer from Methanomicrobiales syntrophic partners can also differ considerably. Altogether, these findings suggest profound differences in electron trafficking in the Methanomicrobiales compared with other hydrogenotrophs, and warrant further functional evaluations.

Horizontal gene transfer: essentiality and evolvability in prokaryotes, and roles in evolutionary transitions

The wide spread of gene exchange and loss in the prokaryotic world has prompted the concept of ‘lateral genomics’ to the point of an outright denial of the relevance of phylogenetic trees for evolution. However, the pronounced coherence congruence of the topologies of numerous gene trees, particularly those for (nearly) universal genes, translates into the notion of a statistical tree of life (STOL), which reflects a central trend of vertical evolution. The STOL can be employed as a framework for reconstruction of the evolutionary processes in the prokaryotic world. Quantitatively, however, horizontal gene transfer (HGT) dominates microbial evolution, with the rate of gene gain and loss being comparable to the rate of point mutations and much greater than the duplication rate. Theoretical models of evolution suggest that HGT is essential for the survival of microbial populations that otherwise deteriorate due to the Muller’s ratchet effect. Apparently, at least some bacteria and archaea evolved dedicated vehicles for gene transfer that evolved from selfish elements such as plasmids and viruses. Recent phylogenomic analyses suggest that episodes of massive HGT were pivotal for the emergence of major groups of organisms such as multiple archaeal phyla as well as eukaryotes. Similar analyses appear to indicate that, in addition to donating hundreds of genes to the emerging eukaryotic lineage, mitochondrial endosymbiosis severely curtailed HGT. These results shed new light on the routes of evolutionary transitions, but caution is due given the inherent uncertainty of deep phylogenies.

Genomics-informed isolation and characterization of a symbiotic Nanoarchaeota system from a terrestrial geothermal environment

Biological features can be inferred, based on genomic data, for many microbial lineages that remain uncultured. However, cultivation is important for characterizing an organism's physiology and testing its genome-encoded potential. Here we use single-cell genomics to infer cultivation conditions for the isolation of an ectosymbiotic Nanoarchaeota (‘Nanopusillus acidilobi') and its host (Acidilobus, a crenarchaeote) from a terrestrial geothermal environment. The cells of ‘Nanopusillus' are among the smallest known cellular organisms (100–300 nm). They appear to have a complete genetic information processing machinery, but lack almost all primary biosynthetic functions as well as respiration and ATP synthesis. Genomic and proteomic comparison with its distant relative, the marine Nanoarchaeum equitans illustrate an ancient, common evolutionary history of adaptation of the Nanoarchaeota to ectosymbiosis, so far unique among the Archaea.

Many microbial lineages have not yet been cultured, which hampers our understanding of their physiology. Here, Wurch et al. use single-cell genomics to infer cultivation conditions for the isolation of a tiny ectosymbiotic nanoarchaeon and its crenarchaeota host from a geothermal spring.

Origin of eukaryotes from within archaea, archaeal eukaryome and bursts of gene gain: eukaryogenesis just made easier?

The origin of eukaryotes is a fundamental, forbidding evolutionary puzzle. Comparative genomic analysis clearly shows that the last eukaryotic common ancestor (LECA) possessed most of the signature complex features of modern eukaryotic cells, in particular the mitochondria, the endomembrane system including the nucleus, an advanced cytoskeleton and the ubiquitin network. Numerous duplications of ancestral genes, e.g. DNA polymerases, RNA polymerases and proteasome subunits, also can be traced back to the LECA. Thus, the LECA was not a primitive organism and its emergence must have resulted from extensive evolution towards cellular complexity. However, the scenario of eukaryogenesis, and in particular the relationship between endosymbiosis and the origin of eukaryotes, is far from being clear. Four recent developments provide new clues to the likely routes of eukaryogenesis. First, evolutionary reconstructions suggest complex ancestors for most of the major groups of archaea, with the subsequent evolution dominated by gene loss. Second, homologues of signature eukaryotic proteins, such as actin and tubulin that form the core of the cytoskeleton or the ubiquitin system, have been detected in diverse archaea. The discovery of this ‘dispersed eukaryome’ implies that the archaeal ancestor of eukaryotes was a complex cell that might have been capable of a primitive form of phagocytosis and thus conducive to endosymbiont capture. Third, phylogenomic analyses converge on the origin of most eukaryotic genes of archaeal descent from within the archaeal evolutionary tree, specifically, the TACK superphylum. Fourth, evidence has been presented that the origin of the major archaeal phyla involved massive acquisition of bacterial genes. Taken together, these findings make the symbiogenetic scenario for the origin of eukaryotes considerably more plausible and the origin of the organizational complexity of eukaryotic cells more readily explainable than they appeared until recently.

Eukaryotes first: how could that be?

In the half century since the formulation of the prokaryote : eukaryote dichotomy, many authors have proposed that the former evolved from something resembling the latter, in defiance of common (and possibly common sense) views. In such 'eukaryotes first' (EF) scenarios, the last universal common ancestor is imagined to have possessed significantly many of the complex characteristics of contemporary eukaryotes, as relics of an earlier 'progenotic' period or RNA world. Bacteria and Archaea thus must have lost these complex features secondarily, through 'streamlining'. If the canonical three-domain tree in which Archaea and Eukarya are sisters is accepted, EF entails that Bacteria and Archaea are convergently prokaryotic. We ask what this means and how it might be tested.

Losing Complexity: The Role of Simplification in Macroevolution

Macroevolutionary patterns can be produced by combinations of diverse and even oppositional dynamics. A growing body of data indicates that secondary simplifications of molecular and cellular structures are common. Some major diversifications in eukaryotes have occurred because of loss and minimalisation; numerous episodes in prokaryote evolution have likewise been driven by the reduction of structure. After examining a range of examples of secondary simplification and its consequences across the tree of life, we address how macroevolutionary explanations might incorporate simplification as well as complexification, and adaptive as well as nonadaptive dynamics.

Diversity and Evolution of Type IV pili Systems in Archaea

Many surface structures in archaea including various types of pili and the archaellum (archaeal flagellum) are homologous to bacterial type IV pili systems (T4P). The T4P consist of multiple proteins, often with poorly conserved sequences, complicating their identification in sequenced genomes. Here we report a comprehensive census of T4P encoded in archaeal genomes using sensitive methods for protein sequence comparison. This analysis confidently identifies as T4P components about 5000 archaeal gene products, 56% of which are currently annotated as hypothetical in public databases. Combining results of this analysis with a comprehensive comparison of genomic neighborhoods of the T4P, we present models of organization of 10 most abundant variants of archaeal T4P. In addition to the differentiation between major and minor pilins, these models include extra components, such as S-layer proteins, adhesins and other membrane and intracellular proteins. For most of these systems, dedicated major pilin families are identified including numerous stand alone major pilin genes of the PilA family. Evidence is presented that secretion ATPases of the T4P and cognate TadC proteins can interact with different pilin sets. Modular evolution of T4P results in combinatorial variability of these systems. Potential regulatory or modulating proteins for the T4P are identified including KaiC family ATPases, vWA domain-containing proteins and the associated MoxR/GvpN ATPase, TFIIB homologs and multiple unrelated transcription regulators some of which are associated specific T4P. Phylogenomic analysis suggests that at least one T4P system was present in the last common ancestor of the extant archaea. Multiple cases of horizontal transfer and lineage-specific duplication of T4P loci were detected. Generally, the T4P of the archaeal TACK superphylum are more diverse and evolve notably faster than those of euryarchaea. The abundance and enormous diversity of T4P in hyperthermophilic archaea present a major enigma. Apparently, fundamental aspects of the biology of hyperthermophiles remain to be elucidated.

Regulation of genetic flux between bacteria by restriction-modification systems

Restriction-modification (R-M) systems are often regarded as bacteria's innate immune systems, protecting cells from infection by mobile genetic elements (MGEs). Their diversification has been recently associated with the emergence of particularly virulent lineages. However, we have previously found more R-M systems in genomes carrying more MGEs. Furthermore, it has been suggested that R-M systems might favor genetic transfer by producing recombinogenic double-stranded DNA ends. To test whether R-M systems favor or disfavor genetic exchanges, we analyzed their frequency with respect to the inferred events of homologous recombination and horizontal gene transfer within 79 bacterial species. Genetic exchanges were more frequent in bacteria with larger genomes and in those encoding more R-M systems. We created a recognition target motif predictor for Type II R-M systems that identifies genomes encoding systems with similar restriction sites. We found more genetic exchanges between these genomes, independently of their evolutionary distance. Our results reconcile previous studies by showing that R-M systems are more abundant in promiscuous species, wherein they establish preferential paths of genetic exchange within and between lineages with cognate R-M systems. Because the repertoire and/or specificity of R-M systems in bacterial lineages vary quickly, the preferential fluxes of genetic transfer within species are expected to constantly change, producing time-dependent networks of gene transfer.

Adaptive engineering of a hyperthermophilic archaeon on CO and discovering the underlying mechanism by multi-omics analysis

The hyperthermophilic archaeon Thermococcus onnurineus NA1 can grow and produce H₂ on carbon monoxide (CO) and its H₂ production rates have been improved through metabolic engineering. In this study, we applied adaptive evolution to enhance H₂ productivity. After over 150 serial transfers onto CO medium, cell density, CO consumption rate and H₂ production rate increased. The underlying mechanism for those physiological changes could be explained by using multi-omics approaches including genomic, transcriptomic and epigenomic analyses. A putative transcriptional regulator was newly identified to regulate the expression levels of genes related to CO oxidation. Transcriptome analysis revealed significant changes in the transcript levels of genes belonging to the categories of transcription, translation and energy metabolism. Our study presents the first genome-scale methylation pattern of hyperthermophilic archaea. Adaptive evolution led to highly enhanced H₂ productivity at high CO flow rates using synthesis gas produced from coal gasification.

The changing view of eukaryogenesis – fossils, cells, lineages and how they all come together

Eukaryogenesis – the emergence of eukaryotic cells – represents a pivotal evolutionary event. With a fundamentally more complex cellular plan compared to prokaryotes, eukaryotes are major contributors to most aspects of life on Earth. For decades, we have understood that eukaryotic origins lie within both the Archaea domain and α-Proteobacteria. However, it is much less clear when, and from which precise ancestors, eukaryotes originated, or the order of emergence of distinctive eukaryotic cellular features. Many competing models for eukaryogenesis have been proposed, but until recently, the absence of discriminatory data meant that a consensus was elusive. Recent advances in paleogeology, phylogenetics, cell biology and microbial diversity, particularly the discovery of the ‘Candidatus Lokiarcheaota’ phylum, are now providing new insights into these aspects of eukaryogenesis. The new data have allowed finessing the time frame during which the events of eukaryogenesis occurred, a more precise identification of the contributing lineages and their likely biological features. The new data have allowed finessing of the time frame during which the events of eukaryogenesis occurred, a more precise identification of the contributing lineages and clarification of their probable biological features.

Early Microbial Evolution: The Age of Anaerobes

In this article, the term "early microbial evolution" refers to the phase of biological history from the emergence of life to the diversification of the first microbial lineages. In the modern era (since we knew about archaea), three debates have emerged on the subject that deserve discussion: (1) thermophilic origins versus mesophilic origins, (2) autotrophic origins versus heterotrophic origins, and (3) how do eukaryotes figure into early evolution. Here, we revisit those debates from the standpoint of newer data. We also consider the perhaps more pressing issue that molecular phylogenies need to recover anaerobic lineages at the base of prokaryotic trees, because O2 is a product of biological evolution; hence, the first microbes had to be anaerobes. If molecular phylogenies do not recover anaerobes basal, something is wrong. Among the anaerobes, hydrogen-dependent autotrophs--acetogens and methanogens--look like good candidates for the ancestral state of physiology in the bacteria and archaea, respectively. New trees tend to indicate that eukaryote cytosolic ribosomes branch within their archaeal homologs, not as sisters to them and, furthermore tend to root archaea within the methanogens. These are major changes in the tree of life, and open up new avenues of thought. Geochemical methane synthesis occurs as a spontaneous, abiotic exergonic reaction at hydrothermal vents. The overall similarity between that reaction and biological methanogenesis fits well with the concept of a methanogenic root for archaea and an autotrophic origin of microbial physiology.

Comparative genomic analysis of evolutionarily conserved but functionally uncharacterized membrane proteins in archaea: Prediction of novel components of secretion, membrane remodeling and glycosylation systems

A systematic comparative genomic analysis of all archaeal membrane proteins that have been projected to the last archaeal common ancestor gene set led to the identification of several novel components of predicted secretion, membrane remodeling, and protein glycosylation systems. Among other findings, most crenarchaea have been shown to encode highly diverged orthologs of the membrane insertase YidC, which is nearly universal in bacteria, eukaryotes, and euryarchaea. We also identified a vast family of archaeal proteins, including the C-terminal domain of N-glycosylation protein AglD, as membrane flippases homologous to the flippase domain of bacterial multipeptide resistance factor MprF, a bifunctional lysylphosphatidylglycerol synthase and flippase. Additionally, several proteins were predicted to function as membrane transporters. The results of this work, combined with our previous analyses, reveal an unexpected diversity of putative archaeal membrane-associated functional systems that remain to be functionally characterized. A more general conclusion from this work is that the currently available collection of archaeal (and bacterial) genomes could be sufficient to identify (almost) all widespread functional modules and develop experimentally testable predictions of their functions.

Nanoarchaeota, Their Sulfolobales Host, and Nanoarchaeota Virus Distribution across Yellowstone National Park Hot Springs

Nanoarchaeota are obligate symbionts with reduced genomes first described from marine thermal vent environments. Here, both community metagenomics and single-cell analysis revealed the presence of Nanoarchaeota in high-temperature (∼90°C), acidic (pH ≈ 2.5 to 3.0) hot springs in Yellowstone National Park (YNP) (United States). Single-cell genome analysis of two cells resulted in two nearly identical genomes, with an estimated full length of 650 kbp. Genome comparison showed that these two cells are more closely related to the recently proposed Nanobsidianus stetteri from a more neutral YNP hot spring than to the marine Nanoarchaeum equitans. Single-cell and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) analysis of environmental hot spring samples identified the host of the YNP Nanoarchaeota as a Sulfolobales species known to inhabit the hot springs. Furthermore, we demonstrate that Nanoarchaeota are widespread in acidic to near neutral hot springs in YNP. An integrated viral sequence was also found within one Nanoarchaeota single-cell genome and further analysis of the purified viral fraction from environmental samples indicates that this is likely a virus replicating within the YNP Nanoarchaeota.

A novel archaeal DNA repair factor that acts with the UvrABC system to repair mitomycin C-induced DNA damage in a PCNA-dependent manner

The sliding clamp proliferating cell nuclear antigen (PCNA) plays a vital role in a number of DNA repair pathways in eukaryotes and archaea by acting as a stable platform onto which other essential protein factors assemble. Many of these proteins interact with PCNA via a short peptide sequence known as a PIP (PCNA interacting protein) motif. Here we describe the identification and functional analysis of a novel PCNA interacting protein NreA that is conserved in the archaea and that has a PIP motif at its C-terminus. Using the genetically tractable euryarchaeon Haloferax volcanii as a model system, we show that the NreA protein is not required for cell viability but that loss of NreA (or replacement of the wild-type protein with a truncated version lacking the C-terminal PIP motif) results in an increased sensitivity to the DNA damaging agent mitomycin C (MMC) that correlates with delayed repair of MMC-induced chromosomal DNA damage monitored by pulsed-field gel electrophoresis. Genetic epistasis analysis in Hfx. volcanii suggests that NreA works together with the UvrABC proteins in repairing DNA damage resulting from exposure to MMC. The wide distribution of NreA family members implies an important role for the protein in DNA damage repair in all archaeal lineages.

Archaeal ancestors of eukaryotes: not so elusive any more

The origin of eukaryotes is one of the hardest problems in evolutionary biology and sometimes raises the ominous specter of irreducible complexity. Reconstruction of the gene repertoire of the last eukaryotic common ancestor (LECA) has revealed a highly complex organism with a variety of advanced features but no detectable evolutionary intermediates to explain their origin. Recently, however, genome analysis of diverse archaea led to the discovery of apparent ancestral versions of several signature eukaryotic systems, such as the actin cytoskeleton and the ubiquitin network, that are scattered among archaea. These findings inspired the hypothesis that the archaeal ancestor of eukaryotes was an unusually complex form with an elaborate intracellular organization. The latest striking discovery made by deep metagenomic sequencing vindicates this hypothesis by showing that in phylogenetic trees eukaryotes fall within a newly identified archaeal group, the Lokiarchaeota, which combine several eukaryotic signatures previously identified in different archaea. The discovery of complex archaea that are the closest living relatives of eukaryotes is most compatible with the symbiogenetic scenario for eukaryogenesis.

Growth-Phase-Specific Modulation of Cell Morphology and Gene Expression by an Archaeal Histone Protein

In all three domains of life, organisms use nonspecific DNA-binding proteins to compact and organize the genome as well as to regulate transcription on a global scale. Histone is the primary eukaryotic nucleoprotein, and its evolutionary roots can be traced to the archaea. However, not all archaea use this protein as the primary DNA-packaging component, raising questions regarding the role of histones in archaeal chromatin function. Here, quantitative phenotyping, transcriptomic, and proteomic assays were performed on deletion and overexpression mutants of the sole histone protein of the hypersaline-adapted haloarchaeal model organism Halobacterium salinarum. This protein is highly conserved among all sequenced haloarchaeal species and maintains hallmark residues required for eukaryotic histone functions. Surprisingly, despite this conservation at the sequence level, unlike in other archaea or eukaryotes, H. salinarum histone is required to regulate cell shape but is not necessary for survival. Genome-wide expression changes in histone deletion strains were global, significant but subtle in terms of fold change, bidirectional, and growth phase dependent. Mass spectrometric proteomic identification of proteins from chromatin enrichments yielded levels of histone and putative nucleoid-associated proteins similar to those of transcription factors, consistent with an open and transcriptionally active genome. Taken together, these data suggest that histone in H. salinarum plays a minor role in DNA compaction but important roles in growth-phase-dependent gene expression and regulation of cell shape. Histone function in haloarchaea more closely resembles a regulator of gene expression than a chromatin-organizing protein like canonical eukaryotic histone.

Histones comprise the major protein component of eukaryotic chromatin and are required for both genome packaging and global regulation of expression. The current paradigm maintains that archaea whose genes encode histone also use these proteins to package DNA. In contrast, here we demonstrate that the sole histone encoded in the genome of the salt-adapted archaeon Halobacterium salinarum is both unessential and unlikely to be involved in DNA compaction despite conservation of residues important for eukaryotic histones. Rather, H. salinarum histone is required for global regulation of gene expression and cell shape. These data are consistent with the hypothesis that H. salinarum histone, strongly conserved across all other known salt-adapted archaea, serves a novel role in gene regulation and cell shape maintenance. Given that archaea possess the ancestral form of eukaryotic histone, this study has important implications for understanding the evolution of histone function.

Assessment of the Carbon Monoxide Metabolism of the Hyperthermophilic Sulfate-Reducing Archaeon Archaeoglobus fulgidus VC-16 by Comparative Transcriptome Analyses

The hyperthermophilic, sulfate-reducing archaeon, Archaeoglobus fulgidus, utilizes CO as an energy source and it is resistant to the toxic effects of high CO concentrations. Herein, transcription profiles were obtained from A. fulgidus during growth with CO and sulfate or thiosulfate, or without an electron acceptor. This provided a basis for a model of the CO metabolism of A. fulgidus. The model suggests proton translocation by “Mitchell-type” loops facilitated by Fqo catalyzing a Fd_red:menaquinone oxidoreductase reaction, as the major mode of energy conservation, rather than formate or H₂ cycling during respiratory growth. The bifunctional CODH (cdhAB-2) is predicted to play an ubiquitous role in the metabolism of CO, and a novel nitrate reductase-associated respiratory complex was induced specifically in the presence of sulfate. A potential role of this complex in relation to Fd_red and APS reduction is discussed. Multiple membrane-bound heterodisulfide reductase (DsrMK) could promote both energy-conserving and non-energy-conserving menaquinol oxidation. Finally, the FqoF subunit may catalyze a Fd_red:F₄₂₀ oxidoreductase reaction. In the absence of electron acceptor, downregulation of F₄₂₀H₂ dependent steps of the acetyl-CoA pathway is linked to transient formate generation. Overall, carboxidotrophic growth seems as an intrinsic capacity of A. fulgidus with little need for novel resistance or respiratory complexes.

The serine/threonine kinase 33 is present and expressed in palaeognath birds but has become a unitary pseudogene in neognaths about 100 million years ago

Background

Serine/threonine kinase 33 (STK33) has been shown to be conserved across all major vertebrate classes including reptiles, mammals, amphibians and fish, suggesting its importance within vertebrates. It has been shown to phosphorylate vimentin and might play a role in spermatogenesis and organ ontogenesis. In this study we analyzed the genomic locus and expression of stk33 in the class Aves, using a combination of large scale next generation sequencing data analysis and traditional PCR.

Results

Within the subclass Palaeognathae we analyzed the white-throated tinamou (Tinamus guttatus), the African ostrich (Struthio camelus) and the emu (Dromaius novaehollandiae). For the African ostrich we were able to generate a 62,778 bp long genomic contig and an mRNA sequence that encodes a protein showing highly significant similarity to STK33 proteins from other vertebrates. The emu has been shown to encode and transcribe a functional STK33 as well. For the white-throated tinamou we were able to identify 13 exons by sequence comparison encoding a protein similar to STK33 as well.

In contrast, in all 28 neognath birds analyzed, we could not find evidence for the existence of a functional copy of stk33 or its expression. In the genomes of these 28 bird species, we found only remnants of the stk33 locus carrying several large genomic deletions, leading to the loss of multiple exons. The remaining exons have acquired various indels and premature stop codons.

Conclusions

We were able to elucidate and describe the genomic structure and the transcription of a functional stk33 gene within the subclass Palaeognathae, but we could only find degenerate remnants of stk33 in all neognath birds analyzed. This led us to the conclusion that stk33 became a unitary pseudogene in the evolutionary history of the class Aves at the paleognath-neognath branch point during the late cretaceous period about 100 million years ago. We hypothesize that the pseudogenization of stk33 might have become fixed in neognaths due to either genetic redundancy or a non-orthologous gene displacement and present potential candidate genes for such an incident.

Electronic supplementary material

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

Application of Subspace Clustering in DNA Sequence Analysis

Identification and clustering of orthologous genes plays an important role in developing evolutionary models such as validating convergent and divergent phylogeny and predicting functional proteins in newly sequenced species of unverified nucleotide protein mappings. Here, we introduce an application of subspace clustering as applied to orthologous gene sequences and discuss the initial results. The working hypothesis is based upon the concept that genetic changes between nucleotide sequences coding for proteins among selected species and groups may lie within a union of subspaces for clusters of the orthologous groups. Estimates for the subspace dimensions were computed for a small population sample. A series of experiments was performed to cluster randomly selected sequences. The experimental design allows for both false positives and false negatives, and estimates for the statistical significance are provided. The clustering results are consistent with the main hypothesis. A simple random mutation binary tree model is used to simulate speciation events that show the interdependence of the subspace rank versus time and mutation rates. The simple mutation model is found to be largely consistent with the observed subspace clustering singular value results. Our study indicates that the subspace clustering method may be applied in orthology analysis.

Comparative genomics reveals conserved positioning of essential genomic clusters in highly rearranged Thermococcales chromosomes

The genomes of the 21 completely sequenced Thermococcales display a characteristic high level of rearrangements. As a result, the prediction of their origin and termination of replication on the sole basis of chromosomal DNA composition or skew is inoperative. Using a different approach based on biologically relevant sequences, we were able to determine oriC position in all 21 genomes. The position of dif, the site where chromosome dimers are resolved before DNA segregation could be predicted in 19 genomes. Computation of the core genome uncovered a number of essential gene clusters with a remarkably stable chromosomal position across species, in sharp contrast with the scrambled nature of their genomes. The active chromosomal reorganization of numerous genes acquired by horizontal transfer, mainly from mobile elements, could explain this phenomenon.

•

Thermococcales chromosomal landmarks were uncovered using biologically relevant sequences.

•

Core genomes procedures predict integration of mobile elements on Thermococcales chromosomes.

•

Thermococcales genomes are highly rearranged but core clusters positions remain invariable.

•

Thermococcales core genes are more expressed and predominantly encoded on the leading strand.

Uncultivated thermophiles: current status and spotlight on 'Aigarchaeota'

Meta-analysis of cultivation-independent sequence data shows that geothermal systems host an abundance of novel organisms, representing a vast unexplored phylogenetic and functional diversity among yet-uncultivated thermophiles. A number of thermophiles have recently been interrogated using metagenomic and/or single-cell genomic approaches, including members of taxonomic groups that inhabit both thermal and non-thermal environments, such as 'Acetothermia' (OP1) and 'Atribacteria' (OP9/JS1), as well as the exclusively thermophilic lineages 'Korarchaeota', 'Calescamantes' (EM19), 'Fervidibacteria' (OctSpA1-106), and 'Aigarchaeota' (HWCG-I). The 'Aigarchaeota', a sister lineage to the Thaumarchaeota, likely includes both hyperthermophiles and moderate thermophiles. They inhabit terrestrial, marine, and subsurface thermal environments and comprise at least nine genus-level lineages, several of which are globally distributed.

A primase subunit essential for efficient primer synthesis by an archaeal eukaryotic-type primase

Archaea encode a eukaryotic-type primase comprising a catalytic subunit (PriS) and a noncatalytic subunit (PriL). Here we report the identification of a primase noncatalytic subunit, denoted PriX, from the hyperthermophilic archaeon Sulfolobus solfataricus. Like PriL, PriX is essential for the survival of the organism. The crystallographic analysis complemented by sensitive sequence comparisons shows that PriX is a diverged homologue of the C-terminal domain of PriL but lacks the iron-sulfur cluster. Phylogenomic analysis provides clues on the origin and evolution of PriX. PriX, PriL and PriS form a stable heterotrimer (PriSLX). Both PriSX and PriSLX show far greater affinity for nucleotide substrates and are substantially more active in primer synthesis than the PriSL heterodimer. In addition, PriL, but not PriX, facilitates primer extension by PriS. We propose that the catalytic activity of PriS is modulated through concerted interactions with the two noncatalytic subunits in primer synthesis.

Bacterial gene import and mesophilic adaptation in archaea

It is widely believed that the archaeal ancestor was hyperthermophilic, but during archaeal evolution, several lineages - including haloarchaea and their sister methanogens, the Thaumarchaeota, and the uncultured Marine Group II and Marine Group III Euryarchaeota (MGII/III) - independently adapted to lower temperatures. Recent phylogenomic studies suggest that the ancestors of these lineages were recipients of massive horizontal gene transfer from bacteria. Many of the acquired genes, which are often involved in metabolism and cell envelope biogenesis, were convergently acquired by distant mesophilic archaea. In this Opinion article, we explore the intriguing hypothesis that the import of these bacterial genes was crucial for the adaptation of archaea to mesophilic lifestyles.

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 manual curation strategy to improve genome annotation: application to a set of haloarchael genomes

Genome annotation errors are a persistent problem that impede research in the biosciences. A manual curation effort is described that attempts to produce high-quality genome annotations for a set of haloarchaeal genomes (Halobacterium salinarum and Hbt. hubeiense, Haloferax volcanii and Hfx. mediterranei, Natronomonas pharaonis and Nmn. moolapensis, Haloquadratum walsbyi strains HBSQ001 and C23, Natrialba magadii, Haloarcula marismortui and Har. hispanica, and Halohasta litchfieldiae). Genomes are checked for missing genes, start codon misassignments, and disrupted genes. Assignments of a specific function are preferably based on experimentally characterized homologs (Gold Standard Proteins). To avoid overannotation, which is a major source of database errors, we restrict annotation to only general function assignments when support for a specific substrate assignment is insufficient. This strategy results in annotations that are resistant to the plethora of errors that compromise public databases. Annotation consistency is rigorously validated for ortholog pairs from the genomes surveyed. The annotation is regularly crosschecked against the UniProt database to further improve annotations and increase the level of standardization. Enhanced genome annotations are submitted to public databases (EMBL/GenBank, UniProt), to the benefit of the scientific community. The enhanced annotations are also publically available via HaloLex.

Giving microbes their due – animal life in a microbially dominant world

The new technology of next-generation sequencing is changing our perceptions of the form and function of the biological world. The emerging data reveal an array of microbes that is more vast and more central to all biological processes than previously appreciated. Further, evidence is accumulating that the alliances of microbes with one another and with constituents of the macrobiological world are critical for the health of the biosphere. This contribution summarizes the basic arguments as to why, when considering the biochemical adaptations of animals, we should integrate the roles of their microbial partners.

Lengths of Orthologous Prokaryotic Proteins Are Affected by Evolutionary Factors

Proteins of the same functional family (for example, kinases) may have significantly different lengths. It is an open question whether such variation in length is random or it appears as a response to some unknown evolutionary driving factors. The main purpose of this paper is to demonstrate existence of factors affecting prokaryotic gene lengths. We believe that the ranking of genomes according to lengths of their genes, followed by the calculation of coefficients of association between genome rank and genome property, is a reasonable approach in revealing such evolutionary driving factors. As we demonstrated earlier, our chosen approach, Bubble-sort, combines stability, accuracy, and computational efficiency as compared to other ranking methods. Application of Bubble Sort to the set of 1390 prokaryotic genomes confirmed that genes of Archaeal species are generally shorter than Bacterial ones. We observed that gene lengths are affected by various factors: within each domain, different phyla have preferences for short or long genes; thermophiles tend to have shorter genes than the soil-dwellers; halophiles tend to have longer genes. We also found that species with overrepresentation of cytosines and guanines in the third position of the codon (GC3 content) tend to have longer genes than species with low GC3 content.

Systems biology approaches to defining transcription regulatory networks in halophilic archaea

To survive complex and changing environmental conditions, microorganisms use gene regulatory networks (GRNs) composed of interacting regulatory transcription factors (TFs) to control the timing and magnitude of gene expression. Genome-wide datasets; such as transcriptomics and protein-DNA interactions; and experiments such as high throughput growth curves; facilitate the construction of GRNs and provide insight into TF interactions occurring under stress. Systems biology approaches integrate these datasets into models of GRN architecture as well as statistical and/or dynamical models to understand the function of networks occurring in cells. Previously, these types of studies have focused on traditional model organisms (e.g. Escherichia coli, yeast). However, recent advances in archaeal genetics and other tools have enabled a systems approach to understanding GRNs in these relatively less studied archaeal model organisms. In this report, we outline a systems biology workflow for generating and integrating data focusing on the TF regulator. We discuss experimental design, outline the process of data collection, and provide the tools required to produce high confidence regulons for the TFs of interest. We provide a case study as an example of this workflow, describing the construction of a GRN centered on multi-TF coordinate control of gene expression governing the oxidative stress response in the hypersaline-adapted archaeon Halobacterium salinarum.

The two-domain tree of life is linked to a new root for the Archaea

One of the most fundamental questions in evolutionary biology is the origin of the lineage leading to eukaryotes. Recent phylogenomic analyses have indicated an emergence of eukaryotes from within the radiation of modern Archaea and specifically from a group comprising Thaumarchaeota/"Aigarchaeota" (candidate phylum)/Crenarchaeota/Korarchaeota (TACK). Despite their major implications, these studies were all based on the reconstruction of universal trees and left the exact placement of eukaryotes with respect to the TACK lineage unclear. Here we have applied an original two-step approach that involves the separate analysis of markers shared between Archaea and eukaryotes and between Archaea and Bacteria. This strategy allowed us to use a larger number of markers and greater taxonomic coverage, obtain high-quality alignments, and alleviate tree reconstruction artifacts potentially introduced when analyzing the three domains simultaneously. Our results robustly indicate a sister relationship of eukaryotes with the TACK superphylum that is strongly associated with a distinct root of the Archaea that lies within the Euryarchaeota, challenging the traditional topology of the archaeal tree. Therefore, if we are to embrace an archaeal origin for eukaryotes, our view of the evolution of the third domain of life will have to be profoundly reconsidered, as will many areas of investigation aimed at inferring ancestral characteristics of early life and Earth.

Complex archaea that bridge the gap between prokaryotes and eukaryotes

The origin of the eukaryotic cell remains one of the most contentious puzzles in modern biology. Recent studies have provided support for the emergence of the eukaryotic host cell from within the archaeal domain of life, but the identity and nature of the putative archaeal ancestor remain a subject of debate. Here we describe the discovery of 'Lokiarchaeota', a novel candidate archaeal phylum, which forms a monophyletic group with eukaryotes in phylogenomic analyses, and whose genomes encode an expanded repertoire of eukaryotic signature proteins that are suggestive of sophisticated membrane remodelling capabilities. Our results provide strong support for hypotheses in which the eukaryotic host evolved from a bona fide archaeon, and demonstrate that many components that underpin eukaryote-specific features were already present in that ancestor. This provided the host with a rich genomic 'starter-kit' to support the increase in the cellular and genomic complexity that is characteristic of eukaryotes.

The Turbulent Network Dynamics of Microbial Evolution and the Statistical Tree of Life

The wide spread and high rate of gene exchange and loss in the prokaryotic world translate into “network genomics”. The rates of gene gain and loss are comparable with the rate of point mutations but are substantially greater than the duplication rate. Thus, evolution of prokaryotes is primarily shaped by gene gain and loss. These processes are essential to prevent mutational meltdown of microbial populations by stopping Muller’s ratchet and appear to trigger emergence of major novel clades by opening up new ecological niches. At least some bacteria and archaea seem to have evolved dedicated devices for gene transfer. Despite the dominance of gene gain and loss, evolution of genes is intrinsically tree-like. The significant coherence between the topologies of numerous gene trees, particularly those for (nearly) universal genes, is compatible with the concept of a statistical tree of life, which forms the framework for reconstruction of the evolutionary processes in the prokaryotic world.

Archaeal Clusters of Orthologous Genes (arCOGs): An Update and Application for Analysis of Shared Features between Thermococcales, Methanococcales, and Methanobacteriales

With the continuously accelerating genome sequencing from diverse groups of archaea and bacteria, accurate identification of gene orthology and availability of readily expandable clusters of orthologous genes are essential for the functional annotation of new genomes. We report an update of the collection of archaeal Clusters of Orthologous Genes (arCOGs) to cover, on average, 91% of the protein-coding genes in 168 archaeal genomes. The new arCOGs were constructed using refined algorithms for orthology identification combined with extensive manual curation, including incorporation of the results of several completed and ongoing research projects in archaeal genomics. A new level of classification is introduced, superclusters that unit two or more arCOGs and more completely reflect gene family evolution than individual, disconnected arCOGs. Assessment of the current archaeal genome annotation in public databases indicates that consistent use of arCOGs can significantly improve the annotation quality. In addition to their utility for genome annotation, arCOGs also are a platform for phylogenomic analysis. We explore this aspect of arCOGs by performing a phylogenomic study of the Thermococci that are traditionally viewed as the basal branch of the Euryarchaeota. The results of phylogenomic analysis that involved both comparison of multiple phylogenetic trees and a search for putative derived shared characters by using phyletic patterns extracted from the arCOGs reveal a likely evolutionary relationship between the Thermococci, Methanococci, and Methanobacteria. The arCOGs are expected to be instrumental for a comprehensive phylogenomic study of the archaea.

No evidence of inhibition of horizontal gene transfer by CRISPR-Cas on evolutionary timescales

The CRISPR (clustered, regularly, interspaced, short, palindromic repeats)–Cas (CRISPR-associated genes) systems of archaea and bacteria provide adaptive immunity against viruses and other selfish elements and are believed to curtail horizontal gene transfer (HGT). Limiting acquisition of new genetic material could be one of the sources of the fitness cost of CRISPR–Cas maintenance and one of the causes of the patchy distribution of CRISPR–Cas among bacteria, and across environments. We sought to test the hypothesis that the activity of CRISPR–Cas in microbes is negatively correlated with the extent of recent HGT. Using three independent measures of HGT, we found no significant dependence between the length of CRISPR arrays, which reflects the activity of the immune system, and the estimated number of recent HGT events. In contrast, we observed a significant negative dependence between the estimated extent of HGT and growth temperature of microbes, which could be explained by the lower genetic diversity in hotter environments. We hypothesize that the relevant events in the evolution of resistance to mobile elements and proclivity for HGT, to which CRISPR–Cas systems seem to substantially contribute, occur on the population scale rather than on the timescale of species evolution.

Extending the conserved phylogenetic core of archaea disentangles the evolution of the third domain of life

Initial studies of the archaeal phylogeny relied mainly on the analysis of the RNA component of the small subunit of the ribosome (SSU rRNA). The resulting phylogenies have provided interesting but partial information on the evolutionary history of the third domain of life because SSU rRNA sequences do not contain enough phylogenetic signal to resolve all nodes of the archaeal tree. Thus, many relationships, and especially the most ancient ones, remained elusive. Moreover, SSU rRNA phylogenies can be heavily biased by tree reconstruction artifacts. The sequencing of complete genomes allows using a variety of protein markers as an alternative to SSU rRNA. Taking advantage of the recent burst of archaeal complete genome sequences, we have carried out an in-depth phylogenomic analysis of this domain. We have identified 200 new protein families that, in addition to the ribosomal proteins and the subunits of the RNA polymerase, form a conserved phylogenetic core of archaeal genes. The accurate analysis of these markers combined with desaturation approaches shed new light on the evolutionary history of Archaea and reveals that several relationships recovered in recent analyses are likely the consequence of tree reconstruction artifacts. Among others, we resolve a number of important relationships, such as those among methanogens Class I, and we propose the definition of two new superclasses within the Euryarchaeota: Methanomada and Diaforarchaea.

A regulatory hierarchy controls the dynamic transcriptional response to extreme oxidative stress in archaea

Networks of interacting transcription factors are central to the regulation of cellular responses to abiotic stress. Although the architecture of many such networks has been mapped, their dynamic function remains unclear. Here we address this challenge in archaea, microorganisms possessing transcription factors that resemble those of both eukaryotes and bacteria. Using genome-wide DNA binding location analysis integrated with gene expression and cell physiological data, we demonstrate that a bacterial-type transcription factor (TF), called RosR, and five TFIIB proteins, homologs of eukaryotic TFs, combinatorially regulate over 100 target genes important for the response to extremely high levels of peroxide. These genes include 20 other transcription factors and oxidative damage repair genes. RosR promoter occupancy is surprisingly dynamic, with the pattern of target gene expression during the transition from rapid growth to stress correlating strongly with the pattern of dynamic binding. We conclude that a hierarchical regulatory network orchestrated by TFs of hybrid lineage enables dynamic response and survival under extreme stress in archaea. This raises questions regarding the evolutionary trajectory of gene networks in response to stress.

Complex circuits of genes rather than a single gene underlie many important processes such as disease, development, and cellular damage repair. Although the wiring of many of these circuits has been mapped, how circuits operate in real time to carry out their functions is poorly understood. Here we address these questions by investigating the function of a gene circuit that responds to reactive oxygen species damage in archaea, microorganisms that represent the third domain of life. Members of this domain of life are excellent models for investigating the function and evolution of gene circuits. Components of archaeal regulatory machinery driving gene circuits resemble those of both bacteria and eukaryotes. Here we demonstrate that regulatory proteins of hybrid ancestry collaborate to control the expression of over 100 genes whose products repair cellular damage. Among these are other regulatory proteins, setting up a stepwise hierarchical circuit that controls damage repair. Regulation is dynamic, with gene targets showing immediate response to damage and restoring normal cellular functions soon thereafter. This study demonstrates how strong environmental forces such as stress may have shaped the wiring and dynamic function of gene circuits, raising important questions regarding how circuits originated over evolutionary time.

Identification of horizontally transferred genes in the genus Colletotrichum reveals a steady tempo of bacterial to fungal gene transfer

BACKGROUND

Horizontal gene transfer (HGT) is the stable transmission of genetic material between organisms by means other than vertical inheritance. HGT has an important role in the evolution of prokaryotes but is relatively rare in eukaryotes. HGT has been shown to contribute to virulence in eukaryotic pathogens. We studied the importance of HGT in plant pathogenic fungi by identifying horizontally transferred genes in the genomes of three members of the genus Colletotrichum.

RESULTS

We identified eleven HGT events from bacteria into members of the genus Colletotrichum or their ancestors. The HGT events include genes involved in amino acid, lipid and sugar metabolism as well as lytic enzymes. Additionally, the putative minimal dates of transference were calculated using a time calibrated phylogenetic tree. This analysis reveals a constant flux of genes from bacteria to fungi throughout the evolution of subphylum Pezizomycotina.

CONCLUSIONS

Genes that are typically transferred by HGT are those that are constantly subject to gene duplication and gene loss. The functions of some of these genes suggest roles in niche adaptation and virulence. We found no evidence of a burst of HGT events coinciding with major geological events. In contrast, HGT appears to be a constant, albeit rare phenomenon in the Pezizomycotina, occurring at a steady rate during their evolution.

Rooting the domain archaea by phylogenomic analysis supports the foundation of the new kingdom Proteoarchaeota

The first 16S rRNA-based phylogenies of the Archaea showed a deep division between two groups, the kingdoms Euryarchaeota and Crenarchaeota. This bipartite classification has been challenged by the recent discovery of new deeply branching lineages (e.g., Thaumarchaeota, Aigarchaeota, Nanoarchaeota, Korarchaeota, Parvarchaeota, Aenigmarchaeota, Diapherotrites, and Nanohaloarchaeota) which have also been given the same taxonomic status of kingdoms. However, the phylogenetic position of some of these lineages is controversial. In addition, phylogenetic analyses of the Archaea have often been carried out without outgroup sequences, making it difficult to determine if these taxa actually define lineages at the same level as the Euryarchaeota and Crenarchaeota. We have addressed the question of the position of the root of the Archaea by reconstructing rooted archaeal phylogenetic trees using bacterial sequences as outgroup. These trees were based on commonly used conserved protein markers (32 ribosomal proteins) as well as on 38 new markers identified through phylogenomic analysis. We thus gathered a total of 70 conserved markers that we analyzed as a concatenated data set. In contrast with previous analyses, our trees consistently placed the root of the archaeal tree between the Euryarchaeota (including the Nanoarchaeota and other fast-evolving lineages) and the rest of archaeal species, which we propose to class within the new kingdom Proteoarchaeota. This implies the relegation of several groups previously classified as kingdoms (e.g., Crenarchaeota, Thaumarchaeota, Aigarchaeota, and Korarchaeota) to a lower taxonomic rank. In addition to taxonomic implications, this profound reorganization of the archaeal phylogeny has also consequences on our appraisal of the nature of the last archaeal ancestor, which most likely was a complex organism with a gene-rich genome.

Potential impact on kidney infection: a whole-genome analysis of Leptospira santarosai serovar Shermani

Leptospira santarosai serovar Shermani is the most frequently encountered serovar, and it causes leptospirosis and tubulointerstitial nephritis in Taiwan. This study aims to complete the genome sequence of L. santarosai serovar Shermani and analyze the transcriptional responses of L. santarosai serovar Shermani to renal tubular cells. To assemble this highly repetitive genome, we combined reads that were generated from four next-generation sequencing platforms by using hybrid assembly approaches to finish two-chromosome contiguous sequences without gaps by validating the data with optical restriction maps and Sanger sequencing. Whole-genome comparison studies revealed a 28-kb region containing genes that encode transposases and hypothetical proteins in L. santarosai serovar Shermani, but this region is absent in other pathogenic Leptospira spp. We found that lipoprotein gene expression in both L. santarosai serovar Shermani and L. interrogans serovar Copenhageni were upregulated upon interaction with renal tubular cells, and LSS19962, a L. santarosai serovar Shermani-specific gene within a 28-kb region that encodes hypothetical proteins, was upregulated in L. santarosai serovar Shermani-infected renal tubular cells. Lipoprotein expression during leptospiral infection might facilitate the interactions of leptospires within kidneys. The availability of the whole-genome sequence of L. santarosai serovar Shermani would make it the first completed sequence of this species, and its comparison with that of other Leptospira spp. may provide invaluable information for further studies in leptospiral pathogenesis.

Expanded microbial genome coverage and improved protein family annotation in the COG database

Microbial genome sequencing projects produce numerous sequences of deduced proteins, only a small fraction of which have been or will ever be studied experimentally. This leaves sequence analysis as the only feasible way to annotate these proteins and assign to them tentative functions. The Clusters of Orthologous Groups of proteins (COGs) database (http://www.ncbi.nlm.nih.gov/COG/), first created in 1997, has been a popular tool for functional annotation. Its success was largely based on (i) its reliance on complete microbial genomes, which allowed reliable assignment of orthologs and paralogs for most genes; (ii) orthology-based approach, which used the function(s) of the characterized member(s) of the protein family (COG) to assign function(s) to the entire set of carefully identified orthologs and describe the range of potential functions when there were more than one; and (iii) careful manual curation of the annotation of the COGs, aimed at detailed prediction of the biological function(s) for each COG while avoiding annotation errors and overprediction. Here we present an update of the COGs, the first since 2003, and a comprehensive revision of the COG annotations and expansion of the genome coverage to include representative complete genomes from all bacterial and archaeal lineages down to the genus level. This re-analysis of the COGs shows that the original COG assignments had an error rate below 0.5% and allows an assessment of the progress in functional genomics in the past 12 years. During this time, functions of many previously uncharacterized COGs have been elucidated and tentative functional assignments of many COGs have been validated, either by targeted experiments or through the use of high-throughput methods. A particularly important development is the assignment of functions to several widespread, conserved proteins many of which turned out to participate in translation, in particular rRNA maturation and tRNA modification. The new version of the COGs is expected to become an important tool for microbial genomics.

Origins of major archaeal clades correspond to gene acquisitions from bacteria

The mechanisms that underlie the origin of major prokaryotic groups are poorly understood. In principle, the origin of both species and higher taxa among prokaryotes should entail similar mechanisms — ecological interactions with the environment paired with natural genetic variation involving lineage-specific gene innovations and lineage-specific gene acquisitions^1,2,3,4. To investigate the origin of higher taxa in archaea, we have determined gene distributions and gene phylogenies for the 267,568 protein coding genes of 134 sequenced archaeal genomes in the context of their homologs from 1,847 reference bacterial genomes. Archaea-specific gene families define 13 traditionally recognized archaeal higher taxa in our sample. Here we report that the origins of these 13 groups unexpectedly correspond to 2,264 group-specific gene acquisitions from bacteria. Interdomain gene transfer is highly asymmetric, transfers from bacteria to archaea are more than 5-fold more frequent than vice versa. Gene transfers identified at major evolutionary transitions among prokaryotes specifically implicate gene acquisitions for metabolic functions from bacteria as key innovations in the origin of higher archaeal taxa.

The neomuran revolution and phagotrophic origin of eukaryotes and cilia in the light of intracellular coevolution and a revised tree of life

Three kinds of cells exist with increasingly complex membrane-protein targeting: Unibacteria (Archaebacteria, Posibacteria) with one cytoplasmic membrane (CM); Negibacteria with a two-membrane envelope (inner CM; outer membrane [OM]); eukaryotes with a plasma membrane and topologically distinct endomembranes and peroxisomes. I combine evidence from multigene trees, palaeontology, and cell biology to show that eukaryotes and archaebacteria are sisters, forming the clade neomura that evolved ~1.2 Gy ago from a posibacterium, whose DNA segregation and cell division were destabilized by murein wall loss and rescued by the evolving novel neomuran endoskeleton, histones, cytokinesis, and glycoproteins. Phagotrophy then induced coevolving serial major changes making eukaryote cells, culminating in two dissimilar cilia via a novel gliding-fishing-swimming scenario. I transfer Chloroflexi to Posibacteria, root the universal tree between them and Heliobacteria, and argue that Negibacteria are a clade whose OM, evolving in a green posibacterium, was never lost.

Genomes in turmoil: quantification of genome dynamics in prokaryote supergenomes

BACKGROUND

Genomes of bacteria and archaea (collectively, prokaryotes) appear to exist in incessant flux, expanding via horizontal gene transfer and gene duplication, and contracting via gene loss. However, the actual rates of genome dynamics and relative contributions of different types of event across the diversity of prokaryotes are largely unknown, as are the sizes of microbial supergenomes, i.e. pools of genes that are accessible to the given microbial species.

RESULTS

We performed a comprehensive analysis of the genome dynamics in 35 groups (34 bacterial and one archaeal) of closely related microbial genomes using a phylogenetic birth-and-death maximum likelihood model to quantify the rates of gene family gain and loss, as well as expansion and reduction. The results show that loss of gene families dominates the evolution of prokaryotes, occurring at approximately three times the rate of gain. The rates of gene family expansion and reduction are typically seven and twenty times less than the gain and loss rates, respectively. Thus, the prevailing mode of evolution in bacteria and archaea is genome contraction, which is partially compensated by the gain of new gene families via horizontal gene transfer. However, the rates of gene family gain, loss, expansion and reduction vary within wide ranges, with the most stable genomes showing rates about 25 times lower than the most dynamic genomes. For many groups, the supergenome estimated from the fraction of repetitive gene family gains includes about tenfold more gene families than the typical genome in the group although some groups appear to have vast, 'open' supergenomes.

CONCLUSIONS

Reconstruction of evolution for groups of closely related bacteria and archaea reveals an extremely rapid and highly variable flux of genes in evolving microbial genomes, demonstrates that extensive gene loss and horizontal gene transfer leading to innovation are the two dominant evolutionary processes, and yields robust estimates of the supergenome size.

Global phylogenomic analysis disentangles the complex evolutionary history of DNA replication in archaea

The archaeal machinery responsible for DNA replication is largely homologous to that of eukaryotes and is clearly distinct from its bacterial counterpart. Moreover, it shows high diversity in the various archaeal lineages, including different sets of components, heterogeneous taxonomic distribution, and a large number of additional copies that are sometimes highly divergent. This has made the evolutionary history of this cellular system particularly challenging to dissect. Here, we have carried out an exhaustive identification of homologs of all major replication components in over 140 complete archaeal genomes. Phylogenomic analysis allowed assigning them to either a conserved and probably essential core of replication components that were mainly vertically inherited, or to a variable and highly divergent shell of extra copies that have likely arisen from integrative elements. This suggests that replication proteins are frequently exchanged between extrachromosomal elements and cellular genomes. Our study allowed clarifying the history that shaped this key cellular process (ancestral components, horizontal gene transfers, and gene losses), providing important evolutionary and functional information. Finally, our precise identification of core components permitted to show that the phylogenetic signal carried by DNA replication is highly consistent with that harbored by two other key informational machineries (translation and transcription), strengthening the existence of a robust organismal tree for the Archaea.

Dark matter in archaeal genomes: a rich source of novel mobile elements, defense systems and secretory complexes

Microbial genomes encompass a sizable fraction of poorly characterized, narrowly spread fast-evolving genes. Using sensitive methods for sequences comparison and protein structure prediction, we performed a detailed comparative analysis of clusters of such genes, which we denote "dark matter islands", in archaeal genomes. The dark matter islands comprise up to 20% of archaeal genomes and show remarkable heterogeneity and diversity. Nevertheless, three classes of entities are common in these genomic loci: (a) integrated viral genomes and other mobile elements; (b) defense systems, and (c) secretory and other membrane-associated systems. The dark matter islands in the genome of thermophiles and mesophiles show similar general trends of gene content, but thermophiles are substantially enriched in predicted membrane proteins whereas mesophiles have a greater proportion of recognizable mobile elements. Based on this analysis, we predict the existence of several novel groups of viruses and mobile elements, previously unnoticed variants of CRISPR-Cas immune systems, and new secretory systems that might be involved in stress response, intermicrobial conflicts and biogenesis of novel, uncharacterized membrane structures.

Evolution of replicative DNA polymerases in archaea and their contributions to the eukaryotic replication machinery

The elaborate eukaryotic DNA replication machinery evolved from the archaeal ancestors that themselves show considerable complexity. Here we discuss the comparative genomic and phylogenetic analysis of the core replication enzymes, the DNA polymerases, in archaea and their relationships with the eukaryotic polymerases. In archaea, there are three groups of family B DNA polymerases, historically known as PolB1, PolB2 and PolB3. All three groups appear to descend from the last common ancestors of the extant archaea but their subsequent evolutionary trajectories seem to have been widely different. Although PolB3 is present in all archaea, with the exception of Thaumarchaeota, and appears to be directly involved in lagging strand replication, the evolution of this gene does not follow the archaeal phylogeny, conceivably due to multiple horizontal transfers and/or dramatic differences in evolutionary rates. In contrast, PolB1 is missing in Euryarchaeota but otherwise seems to have evolved vertically. The third archaeal group of family B polymerases, PolB2, includes primarily proteins in which the catalytic centers of the polymerase and exonuclease domains are disrupted and accordingly the enzymes appear to be inactivated. The members of the PolB2 group are scattered across archaea and might be involved in repair or regulation of replication along with inactivated members of the RadA family ATPases and an additional, uncharacterized protein that are encoded within the same predicted operon. In addition to the family B polymerases, all archaea, with the exception of the Crenarchaeota, encode enzymes of a distinct family D the origin of which is unclear. We examine multiple considerations that appear compatible with the possibility that family D polymerases are highly derived homologs of family B. The eukaryotic DNA polymerases show a highly complex relationship with their archaeal ancestors including contributions of proteins and domains from both the family B and the family D archaeal polymerases.