Molecular identification of filterable bacteria and archaea in the water of acidic lakes of northern Russia

Wetland ecosystems are the natural centers of freshwater formation in northern Russia lowland landscapes. The humic acidic waters formed in bogs feed the numerous lakes of the northern regions. One milliliter of the water in these lakes contains up to 10⁴ ultrasmall microbial cells that pass through “bacterial” filters with a pore size of 0.22 μm. The vast majority of these cells do not grow on nutrient media and cannot be identified by routine cultivation-based approaches. Their identification was performed by analysis of clone libraries obtained by PCR amplification of archaeal and bacterial 16S rRNA genes from the fraction of cells collected from water filtrates of acidic lakes. Most of the obtained bacterial 16S rRNA gene sequences represented the class Betaproteobacteria and exhibited the highest homology of (94–99%) with 16S rRNA genes of representatives of the genera Herbaspirillum, Herminiimonas, Curvibacter, and Burkholderia. The archaeal 16S rRNA gene clone library comprised genes of Euryarchaeota representatives. One-third of these genes exhibited 97–99% homology to the 16S rRNA genes of taxonomically described organisms of the orders Methanobacteriales and Methanosarcinales. The rest of the cloned archaeal 16S rRNA genes were only distantly related (71–74% homology) to those in all earlier characterized archaea.     

DNA Repair Systems in Archaea: Mementos from the Last Universal Common Ancestor?

DNA repair in the Archaea is relevant to the consideration of genome maintenance and replication fidelity in the last universal common ancestor (LUCA) from two perspectives. First, these prokaryotes embody a mix of bacterial and eukaryal molecular features. Second, DNA repair proteins would have been essential in LUCA to maintain genome integrity, regardless of the environmental temperature. Yet we know very little of the basic molecular mechanisms of DNA damage and repair in the Archaea in general. Many studies on DNA repair in archaea have been conducted with hyperthermophiles because of the additional stress imposed on their macromolecules by high temperatures. In addition, of the six complete archaeal genome sequences published so far, five are thermophilic archaea. We have recently shown that the hyperthermophile Pyrococcus furiosus has an extraordinarily high capacity for repair of radiation-induced double-strand breaks and we have identified and sequenced several genes involved in DNA repair in P. furiosus. At the sequence level, only a few genes share homology with known bacterial repair genes. For instance, our phylogenetic analysis indicates that archaeal recombinases occur in two paralogous gene families, one of which is very deeply branched, and both recombinases are more closely related to the eukaryotic RAD51 and Dmc1 gene families than to the Escherichia coli recA gene. We have also identified a gene encoding a repair endo/exonuclease in the genomes of several Archaea. The archaeal sequences are highly homologous to those of the eukaryotic Rad2 family and they cluster with genes of the FEN-1 subfamily, which are known to be involved in DNA replication and repair in eukaryotes. We argue that there is a commonality of mechanisms and protein sequences, shared between prokaryotes and eukaryotes for several modes of DNA repair, reflecting diversification from a minimal set of genes thought to represent the genome of the LUCA.     

Horizontal transfer of a eukaryotic plastid-targeted protein gene to cyanobacteria

Background  

Horizontal or lateral transfer of genetic material between distantly related prokaryotes has been shown to play a major role in the evolution of bacterial and archaeal genomes, but exchange of genes between prokaryotes and eukaryotes is not as well understood. In particular, gene flow from eukaryotes to prokaryotes is rarely documented with strong support, which is unusual since prokaryotic genomes appear to readily accept foreign genes.     

Bacteria rather than Archaea dominate microbial ammonia oxidation in an agricultural soil

Agricultural ecosystems annually receive approximately 25% of the global nitrogen input, much of which is oxidized at least once by ammonia-oxidizing prokaryotes to complete the nitrogen cycle. Recent discoveries have expanded the known ammonia-oxidizing prokaryotes from the domain Bacteria to Archaea . However, in the complex soil environment it remains unclear whether ammonia oxidation is exclusively or predominantly linked to Archaea as implied by their exceptionally high abundance. Here we show that Bacteria rather than Archaea functionally dominate ammonia oxidation in an agricultural soil, despite the fact that archaeal versus bacterial amoA genes are numerically more dominant. In soil microcosms, in which ammonia oxidation was stimulated by ammonium and inhibited by acetylene, activity change was paralleled by abundance change of bacterial but not of archaeal amoA gene copy numbers. Molecular fingerprinting of amoA genes also coupled ammonia oxidation activity with bacterial but not archaeal amoA gene patterns. DNA-stable isotope probing demonstrated CO₂ assimilation by Bacteria rather than Archaea . Our results indicate that Archaea were not important for ammonia oxidation in the agricultural soil tested.     

Recent duplication and inter-locus gene conversion in major histocompatibility class II genes in a teleost,the three-spined stickleback

Using a bacterial artificial chromosome (BAC) library, we analysed a 99.5 kb genomic segment containing the major histocompatibility class II genes of a teleost, the three-spined stickleback Gasterosteus aculeatus. Experiments with G. aculeatus have provided direct evidence for balancing selection by pathogens and mate choice driving MH class II beta polymorphism. Two sets of paralogous class II alpha genes and beta genes in a tandem arrangement were identified, designated Gaac-DAA/DAB and Gaac-DBA/DBB. Expression analysis of the beta genes using single-strand conformation polymorphism revealed that both gene copies are expressed. Based on an analysis of pairwise nucleotide polymorphisms, we estimate that the duplication into two paralogous class II loci occurred only 1.2–2.4 million years ago, 1–2 orders of magnitude more recently than in other fish, bird or mammalian species. At the 3-direction of the classical MH loci, we identified another seven genes or gene fragments, two of which (small inducible cytokine, complement regulatory factor) are related to immune function in other vertebrates. None of these genes were associated with MH class II genes in zebrafish, suggesting a markedly different organisation of the MH class II region in sticklebacks, and presumably, across bony fishes. When the nucleotide substitution pattern of the novel class II beta genes was analysed together with a representative sequence sample isolated from fish in northern Germany (n=27), we found that the peptide binding region of the Gaac-DAB and Gaac-DBB loci had undergone an inter-locus gene conversion (P=0.007). In accordance, we found a 10- to 20-fold higher frequency of CpG-islands on the MH class II segment compared to other species, a feature that may be conducive for inter-locus recombination.     

Distribution and Diversity of Archaeal Ammonia Monooxygenase Genes Associated with Corals

Corals are known to harbor diverse microbial communities of Bacteria and Archaea, yet the ecological role of these microorganisms remains largely unknown. Here we report putative ammonia monooxygenase subunit A (amoA) genes of archaeal origin associated with corals. Multiple DNA samples drawn from nine coral species and four different reef locations were PCR screened for archaeal and bacterial amoA genes, and archaeal amoA gene sequences were obtained from five different species of coral collected in Bocas del Toro, Panama. The 210 coral-associated archaeal amoA sequences recovered in this study were broadly distributed phylogenetically, with most only distantly related to previously reported sequences from coastal/estuarine sediments and oceanic water columns. In contrast, the bacterial amoA gene could not be amplified from any of these samples. These results offer further evidence for the widespread presence of the archaeal amoA gene in marine ecosystems, including coral reefs.     

Gene transfer agents: phage-like elements of genetic exchange

Horizontal gene transfer is important in the evolution of bacterial and archaeal genomes. An interesting genetic exchange process is carried out by diverse phage-like gene transfer agents (GTAs) that are found in a wide range of prokaryotes. Although GTAs resemble phages, they lack the hallmark capabilities that define typical phages, and they package random pieces of the producing cell's genome. In this Review, we discuss the defining characteristics of the GTAs that have been identified to date, along with potential functions for these agents and the possible evolutionary forces that act on the genes involved in their production.     

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.     

Horizontal Transfer of Archaeal Genes into the Deinococcaceae: Detection by Molecular and Computer-Based Approaches

Members of the Deinococcaceae (e.g., Thermus, Meiothermus, Deinococcus) contain A/V-ATPases typically found in Archaea or Eukaryotes which were probably acquired by horizontal gene transfer. Two methods were used to quantify the extent to which archaeal or eukaryotic genes have been acquired by this lineage. Screening of a Meiothermus ruber library with probes made against Thermoplasma acidophilum DNA yielded a number of clones which hybridized more strongly than background. One of these contained the prolyl tRNA synthetase (RS) gene. Phylogenetic analysis shows the M. ruber and D. radiodurans prolyl RS to be more closely related to archaeal and eukaryal forms of this gene than to the typical bacterial type. Using a bioinformatics approach, putative open reading frames (ORFs) from the prerelease version of the D. radiodurans genome were screened for genes more closely related to archaeal or eukaryotic genes. Putative ORFs were searched against representative genomes from each of the three domains using automated BLAST. ORFs showing the highest matches against archaeal and eukaryotic genes were collected and ranked. Among the top-ranked hits were the A/V-ATPase catalytic and noncatalytic subunits and the prolyl RS genes. Using phylogenetic methods, ORFs were analyzed and trees assessed for evidence of horizontal gene transfer. Of the 45 genes examined, 20 showed topologies in which D. radiodurans homologues clearly group with eukaryotic or archaeal homologues, and 17 additional trees were found to show probable evidence of horizontal gene transfer. Compared to the total number of ORFs in the genome, those that can be identified as having been acquired from Archaea or Eukaryotes are relatively few (approximately 1%), suggesting that interdomain transfer is rare.     

Comparative analysis of ribosomal proteins in complete genomes: an example of reductive evolution at the domain scale

A comprehensive investigation of ribosomal genes in complete genomes from 66 different species allows us to address the distribution of r-proteins between and within the three primary domains. Thirty-four r-protein families are represented in all domains but 33 families are specific to Archaea and Eucarya, providing evidence for specialisation at an early stage of evolution between the bacterial lineage and the lineage leading to Archaea and Eukaryotes. With only one specific r-protein, the archaeal ribosome appears to be a small-scale model of the eukaryotic one in terms of protein composition. However, the mechanism of evolution of the protein component of the ribosome appears dramatically different in Archaea. In Bacteria and Eucarya, a restricted number of ribosomal genes can be lost with a bias toward losses in intracellular pathogens. In Archaea, losses implicate 15% of the ribosomal genes revealing an unexpected plasticity of the translation apparatus and the pattern of gene losses indicates a progressive elimination of ribosomal genes in the course of archaeal evolution. This first documented case of reductive evolution at the domain scale provides a new framework for discussing the shape of the universal tree of life and the selective forces directing the evolution of prokaryotes.     

Estimation of phylogenetic inconsistencies in the three domains of life

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.     

Thermotoga heats up lateral gene transfer.

The complete sequence of the bacterium Thermotoga maritima genome has revealed a large fraction of genes most closely related to those of archaeal species. This adds to the accumulating evidence that lateral gene transfer is a potent evolutionary force in prokaryotes, though questions of its magnitude remain.     

Moonlighting glutamate formiminotransferases can functionally replace 5-formyltetrahydrofolate cycloligase

5-Formyltetrahydrofolate (5-CHO-THF) is formed by a side reaction of serine hydroxymethyltransferase. Unlike other folates, it is not a one-carbon donor but a potent inhibitor of folate enzymes and must therefore be metabolized. Only 5-CHO-THF cycloligase (5-FCL) is generally considered to do this. However, comparative genomic analysis indicated (i) that certain prokaryotes lack 5-FCL, implying that they have an alternative 5-CHO-THF-metabolizing enzyme, and (ii) that the histidine breakdown enzyme glutamate formiminotransferase (FT) might moonlight in this role. A functional complementation assay for 5-CHO-THF metabolism was developed in Escherichia coli, based on deleting the gene encoding 5-FCL (ygfA). The deletion mutant accumulated 5-CHO-THF and, with glycine as sole nitrogen source, showed a growth defect; both phenotypes were complemented by bacterial or archaeal genes encoding FT. Furthermore, utilization of supplied 5-CHO-THF by Streptococcus pyogenes was shown to require expression of the native FT. Recombinant bacterial and archaeal FTs catalyzed formyl transfer from 5-CHO-THF to glutamate, with k(cat) values of 0.1-1.2 min(-1) and K(m) values for 5-CHO-THF and glutamate of 0.4-5 μM and 0.03-1 mM, respectively. Although the formyltransferase activities of these proteins were far lower than their formiminotransferase activities, the K(m) values for both substrates relative to their intracellular levels in prokaryotes are consistent with significant in vivo flux through the formyltransferase reaction. Collectively, these data indicate that FTs functionally replace 5-FCL in certain prokaryotes.     

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.     

Aggregate Size Distribution of Ammonia-Oxidizing Bacteria and Archaea at Different Landscape Positions

To quantify the spatial distribution of ammonia-oxidizing bacteria (AOB) and archaea (AOA) and to determine nitrification activity in soil aggregates along a landscape, soil samples were collected from three landscape positions (shoulder, backslope, and toeslope) at two pasture sites with contrasting climatic conditions. The abundance of AOB and AOA was estimated by quantifying their respective bacterial and archaeal amoA gene copies using real-time polymerase chain reaction. Soil organic C (SOC), total N (TN), and the potential nitrification rate (PNR) were measured in aggregate size ranges (4–1, 1–0.25, and 0.25–0.05 mm). At site 1, a decreasing trend in PNR was observed as the size of aggregates decreased. Both bacterial and archaeal amoA genes were higher in the macroaggregates (4–1 and 1–0.25 mm) than in the microaggregates (0.25–0.05 mm) along the landscape. At site 2, PNR was higher in the smallest size of aggregates. In the 0.25–0.05-mm fraction, the abundance of bacterial and archaeal amoA genes was equal to, or greater than, those found in larger aggregate sizes. The relative abundance of archaeal amoA gene and the PNR correlated with relative SOC and TN contents along the landscapes. The positive relationship between relative archaeal amoA gene abundance and PNR suggests that nitrification in the studied pastures is probably driven by ammonia-oxidizing Thaumarchaeota.     

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

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Identification of genes in the genome of the archaeon Methanosarcina mazeii that code for homologs of nuclear eukaryotic molecules involved in RNA processing

A 2.6kb fragment of chromosomal DNA from the archaeon Methanosarcina mazeii was sequenced and analyzed, and it was found to contain coding regions for three proteins that were 321, 234, and 193 amino acids (aa) in length. Homologs of the 321-aa protein were found in all archaeal genomes examined, but not in eukaryotic or bacterial genomes, with one exception in the latter. The protein with 234aa (named PrpM) was most similar to the putative protein Prp31p from Methanobacterium thermoautotrophicum, while the 193-aa protein (named FibM) was identified as an archaeal fibrillarin homolog. Prp and fibrillarin proteins are involved in RNA processing in eukaryotes, but their functions in archaea are not yet understood. The M. mazeii PrpM was also similar to three proteins from Saccharomyces cerevisiae: Prp31p, Nop56p, and Nop58p. Prp31p is a pre-mRNA processing protein, while Nop56p and Nop58p are involved in rRNA processing and interact with fibrillarin. No homologs of either protein were found in bacteria. The archaeal fibrillarin was shorter than its eukaryotic counterpart because it lacked the N-terminal glycine-arginine-rich (GAR) domain, present in most eukaryal homologs. The archaeal prp and fibrillarin gene homologs were found adjacent to each other, whereas in eukarya these genes are on separate chromosomes. Sequence signatures typical of the eukaryal molecules were identified in the M. mazeii and the other archaeal molecules studied. The close proximity of the prp and fib genes raises the possibility of a Prp-fibrillarin interaction in archaea.     

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (± standard deviation) is 197.8 ± 230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.     

Seasonal and Spatial Variability of Bacterial and Archaeal Assemblages in the Coastal Waters near Anvers Island, Antarctica

A previous report of high levels of members of the domain Archaea in Antarctic coastal waters prompted us to investigate the ecology of Antarctic planktonic prokaryotes. rRNA hybridization techniques and denaturing gradient gel electrophoresis (DGGE) analysis of the bacterial V3 region were used to study variation in Antarctic picoplankton assemblages. In Anvers Island nearshore waters during late winter to early spring, the amounts of archaeal rRNA ranged from 17.1 to 3.6% of the total picoplankton rRNA in 1996 and from 16.0 to 1.0% of the total rRNA in 1995. Offshore in the Palmer Basin, the levels of archaeal rRNA throughout the water column were higher (average, 24% of the total rRNA) during the same period in 1996. The archaeal rRNA levels in nearshore waters followed a highly seasonal pattern and markedly decreased during the austral summer at two stations. There was a significant negative correlation between archaeal rRNA levels and phytoplankton levels (as inferred from chlorophyll a concentrations) in nearshore surface waters during the early spring of 1995 and during an 8-month period in 1996 and 1997. In situ hybridization experiments revealed that 5 to 14% of DAPI (4′,6-diamidino-2-phenylindole)-stained cells were archaeal, corresponding to 0.9 × 10⁴ to 2.7 × 10⁴ archaeal cells per ml, in late winter 1996 samples. Analysis of bacterial ribosomal DNA fragments by DGGE revealed that the assemblage composition may reflect changes in water column stability, depth, or season. The data indicate that changes in Antarctic seasons are accompanied by significant shifts in the species composition of bacterioplankton assemblages and by large decreases in the relative proportion of archaeal rRNA in the nearshore water column.     

Symbiont genes in host genomes: fragments with a future?

While lateral transfer is the rule in the evolutionary history of bacterial and archaeal genes, events of transfer from prokaryotes to eukaryotes are rare. Germline-transmitted animal symbionts, such as Wolbachia pipientis, are well placed to participate in such transfers. In a recent issue of Science, Dunning Hotopp et al. identified instances of transfer of Wolbachia DNA to host genomes. It is unknown whether these transfers represent innovation in animal evolution.