Diversity of radA genes from cultured and uncultured archaea: comparative analysis of putative RadA proteins and their use as a phylogenetic marker

Exploring Marine Planktonic Archaea: Then and Now

In 1977, Woese and Fox leveraged molecular phylogenetic analyses of ribosomal RNAs and identified a new microbial domain of life on Earth, the Archaebacteria (now known as Archaea). At the time of their discovery, only one archaebacterial group, the strictly anaerobic methanogens, was known. But soon, other phenotypically unrelated microbial isolates were shown to belong to the Archaea, many originating from extreme habitats, including extreme halophiles, extreme thermophiles, and thermoacidophiles. Since most Archaea seemed to inhabit extreme or strictly anoxic habitats, it came as a surprise in 1992 when two new lineages of archaea were reported to be abundant in oxygen rich, temperate marine coastal waters and the deep ocean. Since that time, studies of marine planktonic archaea have revealed many more surprises, including their unexpected ubiquity, unusual symbiotic associations, unpredicted physiologies and biogeochemistry, and global abundance. In this Perspective, early work conducted on marine planktonic Archaea by my lab group and others is discussed in terms of the relevant historical context, some of the original research motivations, and surprises and discoveries encountered along the way.

Novel Genes Involved in Resistance to Both Ultraviolet Radiation and Perchlorate From the Metagenomes of Hypersaline Environments

Microorganisms that thrive in hypersaline environments on the surface of our planet are exposed to the harmful effects of ultraviolet radiation. Therefore, for their protection, they have sunscreen pigments and highly efficient DNA repair and protection systems. The present study aimed to identify new genes involved in UV radiation resistance from these microorganisms, many of which cannot be cultured in the laboratory. Thus, a functional metagenomic approach was used and for this, small-insert libraries were constructed with DNA isolated from microorganisms of high-altitude Andean hypersaline lakes in Argentina (Diamante and Ojo Seco lakes, 4,589 and 3,200 m, respectively) and from the Es Trenc solar saltern in Spain. The libraries were hosted in a UV radiation-sensitive strain of Escherichia coli (recA mutant) and they were exposed to UVB. The resistant colonies were analyzed and as a result, four clones were identified with environmental DNA fragments containing five genes that conferred resistance to UV radiation in E. coli. One gene encoded a RecA-like protein, complementing the mutation in recA that makes the E. coli host strain more sensitive to UV radiation. Two other genes from the same DNA fragment encoded a TATA-box binding protein and an unknown protein, both responsible for UV resistance. Interestingly, two other genes from different and remote environments, the Ojo Seco Andean lake and the Es Trenc saltern, encoded two hypothetical proteins that can be considered homologous based on their significant amino acid similarity (49%). All of these genes also conferred resistance to 4-nitroquinoline 1-oxide (4-NQO), a compound that mimics the effect of UV radiation on DNA, and also to perchlorate, a powerful oxidant that can induce DNA damage. Furthermore, the hypothetical protein from the Es Trenc salterns was localized as discrete foci possibly associated with damaged sites in the DNA in cells treated with 4-NQO, so it could be involved in the repair of damaged DNA. In summary, novel genes involved in resistance to UV radiation, 4-NQO and perchlorate have been identified in this work and two of them encoding hypothetical proteins that could be involved in DNA damage repair activities not previously described.

Genotypic and Lipid Analyses of Strains From the Archaeal Genus Halorubrum Reveal Insights Into Their Taxonomy, Divergence, and Population Structure

To gain a better understanding of how divergence occurs, and how taxonomy can benefit from studying natural populations, we isolated and examined 25 closely related Halorubrum strains obtained from different hypersaline communities and compared them to validly named species and other reference strains using five taxonomic study approaches: phylogenetic analysis using the 16S rRNA gene and multilocus sequencing analysis (MLSA), polar lipid profiles (PLP), average nucleotide identity (ANI) and DNA-DNA hybridization (DDH). 16S rRNA gene sequence could not differentiate the newly isolated strains from described species, while MLSA grouped strains into three major clusters. Two of those MLSA clusters distinguished candidates for new species. The third cluster with concatenated sequence identity equal to or greater than 97.5% was comprised of strains from Aran-Bidgol Lake (Iran) and solar salterns in Namibia and Spain, and two previously described species isolated from Mexico and Algeria. PLP and DDH analyses showed that Aran-Bidgol strains formed uniform populations, and that strains isolated from other geographic locations were heterogeneous and divergent, indicating that they may constitute different species. Therefore, applying only sequencing approaches and similarity cutoffs for circumscribing species may be too conservative, lumping concealed diversity into a single taxon. Further, our data support the interpretation that local populations experience unique evolutionary homogenization pressures, and once relieved of insular constraints (e.g., through migration) are free to diverge.

Role of O2 in the Growth of Rhizobium leguminosarum bv. viciae 3841 on Glucose and Succinate

Insertion sequencing (INSeq) analysis of Rhizobium leguminosarum bv. viciae 3841 (Rlv3841) grown on glucose or succinate at both 21% and 1% O₂ was used to understand how O₂ concentration alters metabolism. Two transcriptional regulators were required for growth on glucose (pRL120207 [eryD] and RL0547 [phoB]), five were required on succinate (pRL100388, RL1641, RL1642, RL3427, and RL4524 [ecfL]), and three were required on 1% O₂ (pRL110072, RL0545 [phoU], and RL4042). A novel toxin-antitoxin system was identified that could be important for generation of new plasmidless rhizobial strains. Rlv3841 appears to use the methylglyoxal pathway alongside the Entner-Doudoroff (ED) pathway and tricarboxylic acid (TCA) cycle for optimal growth on glucose. Surprisingly, the ED pathway was required for growth on succinate, suggesting that sugars made by gluconeogenesis must undergo recycling. Altered amino acid metabolism was specifically needed for growth on glucose, including RL2082 (gatB) and pRL120419 (opaA, encoding omega-amino acid:pyruvate transaminase). Growth on succinate specifically required enzymes of nucleobase synthesis, including ribose-phosphate pyrophosphokinase (RL3468 [prs]) and a cytosine deaminase (pRL90208 [codA]). Succinate growth was particularly dependent on cell surface factors, including the PrsD-PrsE type I secretion system and UDP-galactose production. Only RL2393 (glnB, encoding nitrogen regulatory protein PII) was specifically essential for growth on succinate at 1% O₂, conditions similar to those experienced by N₂-fixing bacteroids. Glutamate synthesis is constitutively activated in glnB mutants, suggesting that consumption of 2-ketoglutarate may increase flux through the TCA cycle, leading to excess reductant that cannot be reoxidized at 1% O₂ and cell death.

IMPORTANCE

Rhizobium leguminosarum, a soil bacterium that forms N₂-fixing symbioses with several agriculturally important leguminous plants (including pea, vetch, and lentil), has been widely utilized as a model to study Rhizobium-legume symbioses. Insertion sequencing (INSeq) has been used to identify factors needed for its growth on different carbon sources and O₂ levels. Identification of these factors is fundamental to a better understanding of the cell physiology and core metabolism of this bacterium, which adapts to a variety of different carbon sources and O₂ tensions during growth in soil and N₂ fixation in symbiosis with legumes.

Ancient Systems of Sodium/Potassium Homeostasis as Predecessors of Membrane Bioenergetics

Cell cytoplasm of archaea, bacteria, and eukaryotes contains substantially more potassium than sodium, and potassium cations are specifically required for many key cellular processes, including protein synthesis. This distinct ionic composition and requirements have been attributed to the emergence of the first cells in potassium-rich habitats. Different, albeit complementary, scenarios have been proposed for the primordial potassium-rich environments based on experimental data and theoretical considerations. Specifically, building on the observation that potassium prevails over sodium in the vapor of inland geothermal systems, we have argued that the first cells could emerge in the pools and puddles at the periphery of primordial anoxic geothermal fields, where the elementary composition of the condensed vapor would resemble the internal milieu of modern cells. Marine and freshwater environments generally contain more sodium than potassium. Therefore, to invade such environments, while maintaining excess of potassium over sodium in the cytoplasm, primordial cells needed means to extrude sodium ions. The foray into new, sodium-rich habitats was the likely driving force behind the evolution of diverse redox-, light-, chemically-, or osmotically-dependent sodium export pumps and the increase of membrane tightness. Here we present a scenario that details how the interplay between several, initially independent sodium pumps might have triggered the evolution of sodium-dependent membrane bioenergetics, followed by the separate emergence of the proton-dependent bioenergetics in archaea and bacteria. We also discuss the development of systems that utilize the sodium/potassium gradient across the cell membranes.

Systematic identification of gene families for use as "markers" for phylogenetic and phylogeny-driven ecological studies of bacteria and archaea and their major subgroups

With the astonishing rate that genomic and metagenomic sequence data sets are accumulating, there are many reasons to constrain the data analyses. One approach to such constrained analyses is to focus on select subsets of gene families that are particularly well suited for the tasks at hand. Such gene families have generally been referred to as “marker” genes. We are particularly interested in identifying and using such marker genes for phylogenetic and phylogeny-driven ecological studies of microbes and their communities (e.g., construction of species trees, phylogenetic based assignment of metagenomic sequence reads to taxonomic groups, phylogeny-based assessment of alpha- and beta-diversity of microbial communities from metagenomic data). We therefore refer to these as PhyEco (for phylogenetic and phylogenetic ecology) markers. The dual use of these PhyEco markers means that we needed to develop and apply a set of somewhat novel criteria for identification of the best candidates for such markers. The criteria we focused on included universality across the taxa of interest, ability to be used to produce robust phylogenetic trees that reflect as much as possible the evolution of the species from which the genes come, and low variation in copy number across taxa.

We describe here an automated protocol for identifying potential PhyEco markers from a set of complete genome sequences. The protocol combines rapid searching, clustering and phylogenetic tree building algorithms to generate protein families that meet the criteria listed above. We report here the identification of PhyEco markers for different taxonomic levels including 40 for “all bacteria and archaea”, 114 for “all bacteria (greatly expanding on the ∼30 commonly used), and 100 s to 1000 s for some of the individual phyla of bacteria. This new list of PhyEco markers should allow much more detailed automated phylogenetic and phylogenetic ecology analyses of these groups than possible previously.

Heteroduplex formation, mismatch resolution, and genetic sectoring during homologous recombination in the hyperthermophilic archaeon sulfolobus acidocaldarius

Hyperthermophilic archaea exhibit certain molecular-genetic features not seen in bacteria or eukaryotes, and their systems of homologous recombination (HR) remain largely unexplored in vivo. We transformed a Sulfolobus acidocaldariuspyrE mutant with short DNAs that contained multiple non-selected genetic markers within the pyrE gene. From 20 to 40% of the resulting colonies were found to contain two Pyr⁺ clones with distinct sets of the non-selected markers. The dual-genotype colonies could not be attributed to multiple DNAs entering the cells, or to conjugation between transformed and non-transformed cells. These colonies thus appear to represent genetic sectoring in which regions of heteroduplex DNA formed and then segregated after partial resolution of inter-strand differences. Surprisingly, sectoring was also frequent in cells transformed with single-stranded DNAs. Oligonucleotides produced more sectored transformants when electroporated as single strands than as a duplex, although all forms of donor DNA (positive-strand, negative-strand, and duplex) produced a diversity of genotypes, despite the limited number of markers. The marker patterns in the recombinants indicate that S. acidocaldarius resolves individual mismatches through un-coordinated short-patch excision followed by re-filling of the resulting gap. The conversion events that occur during transformation by single-stranded DNA do not show the strand bias necessary for a system that corrects replication errors effectively; similar events also occur in pre-formed heteroduplex electroporated into the cells. Although numerous mechanistic details remain obscure, the results demonstrate that the HR system of S. acidocaldarius can generate remarkable genetic diversity from short intervals of moderately diverged DNAs.

Recombinogenic properties of Pyrococcus furiosus strain COM1 enable rapid selection of targeted mutants

We recently reported the isolation of a mutant of Pyrococcus furiosus, COM1, that is naturally and efficiently competent for DNA uptake. While we do not know the exact nature of this mutation, the combined transformation and recombination frequencies of this strain allow marker replacement by direct selection using linear DNA. In testing the limits of its recombination efficiency, we discovered that marker replacement was possible with as few as 40 nucleotides of flanking homology to the target region. We utilized this ability to design a strategy for selection of constructed deletions using PCR products with subsequent excision, or "pop-out," of the selected marker. We used this method to construct a "markerless" deletion of the trpAB locus in the GLW101 (COM1 ΔpyrF) background to generate a strain (JFW02) that is a tight tryptophan auxotroph, providing a genetic background with two auxotrophic markers for further strain construction. The utility of trpAB as a selectable marker was demonstrated using prototrophic selection of plasmids and genomic DNA containing the wild-type trpAB alleles. A deletion of radB was also constructed that, surprisingly, had no obvious effect on either recombination or transformation, suggesting that its gene product is not involved in the COM1 phenotype. Attempts to construct a radA deletion mutation were unsuccessful, suggesting that this may be an essential gene. The ease and speed of this procedure will facilitate the construction of strains with multiple genetic changes and allow the construction of mutants with deletions of virtually any nonessential gene.

Monoclonal antibodies against the Escherichia coli DNA repair protein RadA/Sms

The RadA/Sms protein facilitates DNA repair in Escherichia coli cells damaged by UV radiation, X-rays, and chemical agents. However, the precise mechanism by which RadA/Sms aids DNA repair is unknown. Here we report the production of monoclonal antibodies (MAbs) specific for RadA/Sms for use in biochemical and physiological investigations. Histidine-tagged RadA/Sms (RadA-6xHis) was overproduced in E. coli BL21 cells transformed with the radA/sms coding region in plasmid pRSET A and purified by nickel affinity chromatography. Splenocytes from female BALB/c mice hyperimmunized with the purified protein were fused to SP2/0-Ag14 myeloma cells, and the resultant hybridomas were selected in HAT medium. MAbs were detected in hybridoma culture supernatants by indirect ELISA and Western blot analysis against purified RadA-6xHis. MAbs from four cell lines were further evaluated by Western blotting against peptide maps generated by endoproteinase Glu-C digestion of RadA-6xHis. Each of the four MAbs recognized a unique epitope on the fusion protein. Two of the MAbs (6F5 and 2A2) also detected wild-type (tagless) RadA/Sms produced from the pJS003 plasmid in E. coli K-12 cells. We anticipate that these antibodies will prove useful for the detection, isolation, and functional analysis of RadA/Sms.

A multilocus sequence analysis approach to the phylogeny and taxonomy of the Halobacteriales

Members of the order Halobacteriales are obligate extreme halophiles that belong to the domain Archaea. The classification of the Halobacteriales currently relies on a polyphasic approach, which integrates phenotypic, genotypic and chemotaxonomic characterization. However, the most utilized genetic marker for phylogeny, the 16S rRNA gene, has multiple drawbacks for use with the Halobacteriales: the species of many genera exhibit large intragenic differences between multiple ribosomal RNA operons, the gene is too conserved to discriminate reliably at the species level and it appears to be the most frequently recombined gene between closely related species. Moreover, the Halobacteriales is a rapidly expanding group due to recent successes at cultivating novel strains from a diverse set of hypersaline environments; a fast, reliable, inexpensive, portable molecular method for discriminating species is required for their investigation. Recently, multilocus sequence analysis (MLSA) has been shown to be an effective tool for strain identification and taxonomic designation, even for those taxa that experience frequent lateral gene transfer and homologous recombination. In this study, MLSA was utilized for evolutionary and taxonomic investigation of the Halobacteriales. Efficacy of the MLSA approach was tested across a hierarchical gradient using 52 halobacterial strains, representing 33 species (including names without standing in nomenclature) and 14 genera. A subset of 21 strains from the genus Haloarcula was analysed separately to test the sensitivity and relevance of the MLSA approach among closely related strains and species. The results demonstrated that MLSA differentiated individual strains, reliably grouped strains into species and species into genera and identified potential novel species and also family-like relationships. This study demonstrates that MLSA is a rapid and informative molecular method that will probably accommodate strain analysis at any taxonomic level within the Halobacteriales.

Stalking the fourth domain in metagenomic data: searching for, discovering, and interpreting novel, deep branches in marker gene phylogenetic trees

Background

Most of our knowledge about the ancient evolutionary history of organisms has been derived from data associated with specific known organisms (i.e., organisms that we can study directly such as plants, metazoans, and culturable microbes). Recently, however, a new source of data for such studies has arrived: DNA sequence data generated directly from environmental samples. Such metagenomic data has enormous potential in a variety of areas including, as we argue here, in studies of very early events in the evolution of gene families and of species.

Methodology/Principal Findings

We designed and implemented new methods for analyzing metagenomic data and used them to search the Global Ocean Sampling (GOS) Expedition data set for novel lineages in three gene families commonly used in phylogenetic studies of known and unknown organisms: small subunit rRNA and the recA and rpoB superfamilies. Though the methods available could not accurately identify very deeply branched ss-rRNAs (largely due to difficulties in making robust sequence alignments for novel rRNA fragments), our analysis revealed the existence of multiple novel branches in the recA and rpoB gene families. Analysis of available sequence data likely from the same genomes as these novel recA and rpoB homologs was then used to further characterize the possible organismal source of the novel sequences.

Conclusions/Significance

Of the novel recA and rpoB homologs identified in the metagenomic data, some likely come from uncharacterized viruses while others may represent ancient paralogs not yet seen in any cultured organism. A third possibility is that some come from novel cellular lineages that are only distantly related to any organisms for which sequence data is currently available. If there exist any major, but so-far-undiscovered, deeply branching lineages in the tree of life, we suggest that methods such as those described herein currently offer the best way to search for them.

Discontinuity and limited linkage in the homologous recombination system of a hyperthermophilic archaeon

Genetic transformation of Sulfolobus acidocaldarius by a multiply marked pyrE gene provided a high-resolution assay of homologous recombination in a hyperthermophilic archaeon. Analysis of 100 Pyr(+) transformants revealed that this recombination system could transfer each of 23 nonselected base pair substitutions to the recipient chromosome along with the selected marker. In 30% of the recombinants, donor markers were transferred as multiple blocks. In at least 40% of the recombinants, donor markers separated by 5 or 6 bp segregated from each other, whereas similar markers separated by 2 bp did not segregate. Among intermarker intervals, the frequency of recombination tract endpoints varied 40-fold, but in contrast to other recombination systems, it did not correlate with the length of the interval. The average length of donor tracts (161 bp) and the frequent generation of multiple tracts seemed generally consistent with the genetic properties observed previously in S. acidocaldarius conjugation. The efficiency with which short intervals of diverged pyrE sequence were incorporated into the genome raises questions about the threat of ectopic recombination in Sulfolobus spp. mediated by this apparently efficient yet permissive system.

The complete genome sequence of Haloferax volcanii DS2, a model archaeon

BACKGROUND

Haloferax volcanii is an easily culturable moderate halophile that grows on simple defined media, is readily transformable, and has a relatively stable genome. This, in combination with its biochemical and genetic tractability, has made Hfx. volcanii a key model organism, not only for the study of halophilicity, but also for archaeal biology in general.

METHODOLOGY/PRINCIPAL FINDINGS

We report here the sequencing and analysis of the genome of Hfx. volcanii DS2, the type strain of this species. The genome contains a main 2.848 Mb chromosome, three smaller chromosomes pHV1, 3, 4 (85, 438, 636 kb, respectively) and the pHV2 plasmid (6.4 kb).

CONCLUSIONS/SIGNIFICANCE

The completed genome sequence, presented here, provides an invaluable tool for further in vivo and in vitro studies of Hfx. volcanii.

RecA family proteins in archaea: RadA and its cousins

Recombinases of the RecA family are essential for homologous recombination and underpin genome stability, by promoting the repair of double-stranded DNA breaks and the rescue of collapsed DNA replication forks. Until now, our understanding of homologous recombination has relied on studies of bacterial and eukaryotic model organisms. Archaea provide new opportunities to study how recombination operates in a lineage distinct from bacteria and eukaryotes. In the present paper, we focus on RadA, the archaeal RecA family recombinase, and its homologues in archaea and other domains. On the basis of phylogenetic analysis, we propose that a family of archaeal proteins with a single RecA domain, which are currently annotated as KaiC, be renamed aRadC.

Discovery and characterization of a new bacterial candidate division by an anaerobic sludge digester metagenomic approach

We have constructed a large fosmid library from a mesophilic anaerobic digester and explored its 16S rDNA diversity using a high-density filter DNA–DNA hybridization procedure. We identified a group of 16S rDNA sequences forming a new bacterial lineage named WWE3 (Waste Water of Evry 3). Only one sequence from the public databases shares a sequence identity above 80% with the WWE3 group which hence cannot be affiliated to any known or candidate prokaryotic division. Despite representing a non-negligible fraction (5% of the 16S rDNA sequences) of the bacterial population of this digester, the WWE3 bacteria could not have been retrieved using the conventional 16S rDNA amplification procedure due to their unusual 16S rDNA gene sequence. WWE3 bacteria were detected by polymerase chain reaction (PCR) in various environments (anaerobic digesters, swine lagoon slurries and freshwater biofilms) using newly designed specific PCR primer sets. Fluorescence in situ hybridization (FISH) analysis of sludge samples showed that WWE3 microorganisms are oval-shaped and located deep inside sludge flocs. Detailed phylogenetic analysis showed that WWE3 bacteria form a distinct monophyletic group deeply branching apart from all known bacterial divisions. A new bacterial candidate division status is proposed for this group.

Identification homologous recombination function from haloarchaea plasmid pHH205

Homologous recombination (HR) was found to be so frequent in haloarchaea that its significance in evolution and diversity of this clade of life might have been underestimated. However, so far there has been no report on recombination function carried on plasmid. Here we report that a 4.8-kb SnaBI-PvuII digested segment from pHH205 might carry such a function. Four constructed plasmids: pUN, pUN-205, pUM and pUM-205, with pUN and pUN205 containing Nov(R) gene, pUM and pUM-205 carrying Mev(R) gene, were used to transform Haloferax volcanii DS52 (radA(-)). The results showed that only pUN-205 and pUM-205 containing the 4.8-kb SnaBI-PvuII digested segment from pHH205 were able to shift Nov(R) and Mev(R) gene into the chromosome of Haloferax volcanii DS52 through HR, whereas those in pUN and pUM could not, which indicated that the segment from pHH205 does contain a recombination function.

DinI and RecX modulate RecA?DNA structures in Escherichia coli K-12

RecA plays a central role in recombination, DNA repair and SOS induction through forming a RecA-DNA helical filament. Biochemical observations show that at low ratios to RecA, DinI and RecX stabilize and destabilize RecA-DNA filaments, respectively, and that the C-terminal 17 residues of RecA are important for RecX function. RecA-DNA filament formation was assayed in vivo using RecA-GFP foci formation in log-phase and UV-irradiated cells. In log-phase cells, dinI mutants have fewer foci than wild type and that recX mutants have more foci than wild type. A recADelta17::gfp mutant had more foci like a recX mutant. dinI recX double mutants have the same number of foci as dinI mutants alone, suggesting that dinI is epistatic to recX. After UV treatment, the dinI, recX and dinI recX mutants differed in their ability to form foci. All three mutants had fewer foci than wild type. The dinI mutant's foci persisted longer than wild-type foci. Roles of DinI and RecX after UV treatment differed from those during log-phase growth and may reflect the different DNA substrates, population of proteins or amounts during the SOS response. These experiments give new insight into the roles of these proteins.

Metagenomics: genomic analysis of microbial communities

Uncultured microorganisms comprise the majority of the planet's biological diversity. Microorganisms represent two of the three domains of life and contain vast diversity that is the product of an estimated 3.8 billion years of evolution. In many environments, as many as 99% of the microorganisms cannot be cultured by standard techniques, and the uncultured fraction includes diverse organisms that are only distantly related to the cultured ones. Therefore, culture-independent methods are essential to understand the genetic diversity, population structure, and ecological roles of the majority of microorganisms. Metagenomics, or the culture-independent genomic analysis of an assemblage of microorganisms, has potential to answer fundamental questions in microbial ecology. This review describes progress toward understanding the biology of uncultured Bacteria, Archaea, and viruses through metagenomic analyses.

Metagenomics: application of genomics to uncultured microorganisms

Metagenomics (also referred to as environmental and community genomics) is the genomic analysis of microorganisms by direct extraction and cloning of DNA from an assemblage of microorganisms. The development of metagenomics stemmed from the ineluctable evidence that as-yet-uncultured microorganisms represent the vast majority of organisms in most environments on earth. This evidence was derived from analyses of 16S rRNA gene sequences amplified directly from the environment, an approach that avoided the bias imposed by culturing and led to the discovery of vast new lineages of microbial life. Although the portrait of the microbial world was revolutionized by analysis of 16S rRNA genes, such studies yielded only a phylogenetic description of community membership, providing little insight into the genetics, physiology, and biochemistry of the members. Metagenomics provides a second tier of technical innovation that facilitates study of the physiology and ecology of environmental microorganisms. Novel genes and gene products discovered through metagenomics include the first bacteriorhodopsin of bacterial origin; novel small molecules with antimicrobial activity; and new members of families of known proteins, such as an Na(+)(Li(+))/H(+) antiporter, RecA, DNA polymerase, and antibiotic resistance determinants. Reassembly of multiple genomes has provided insight into energy and nutrient cycling within the community, genome structure, gene function, population genetics and microheterogeneity, and lateral gene transfer among members of an uncultured community. The application of metagenomic sequence information will facilitate the design of better culturing strategies to link genomic analysis with pure culture studies.

Biotechnological prospects from metagenomics

The recognition that most microorganisms in the environment cannot be cultured by standard methods stimulated the development of metagenomics, which is the genomic analysis of uncultured microorganisms. Two types of analysis have been used to obtain information from metagenomic libraries: a function-driven approach, in which metagenomic libraries are initially screened for an expressed trait, and a sequence-driven approach, in which libraries are initially screened for particular DNA sequences. New antibiotics and enzymes are among the early discoveries from metagenomics. Future refinement of methods that enrich for genes of particular function will accelerate the rate of discovery of useful molecules.

Polyphosphate kinase from activated sludge performing enhanced biological phosphorus removal

A novel polyphosphate kinase (PPK) was retrieved from an uncultivated organism in activated sludge carrying out enhanced biological phosphorus removal (EBPR). Acetate-fed laboratory-scale sequencing batch reactors were used to maintain sludge with a high phosphorus content (approximately 11% of the biomass). PCR-based clone libraries of small subunit rRNA genes and fluorescent in situ hybridization (FISH) were used to verify that the sludge was enriched in Rhodocyclus-like beta-Proteobacteria known to be associated with sludges carrying out EBPR. These organisms comprised approximately 80% of total bacteria in the sludge, as assessed by FISH. Degenerate PCR primers were designed to retrieve fragments of putative ppk genes from a pure culture of Rhodocyclus tenuis and from organisms in the sludge. Four novel ppk homologs were found in the sludge, and two of these (types I and II) shared a high degree of amino acid similarity with R. tenuis PPK (86 and 87% similarity, respectively). Dot blot analysis of total RNA extracted from sludge demonstrated that the Type I ppk mRNA was present, indicating that this gene is expressed during EBPR. Inverse PCR was used to obtain the full Type I sequence from sludge DNA, and a full-length PPK was cloned, overexpressed, and purified to near homogeneity. The purified PPK has a specific activity comparable to that of other PPKs, has a requirement for Mg(2+), and does not appear to operate in reverse. PPK activity was found mainly in the particulate fraction of lysed sludge microorganisms.

Heterogeneous yet similar introns reside in identical positions of the rRNA genes in natural isolates of the archaeon Aeropyrum pernix

Some archaeal ribosomal DNA (rDNA) introns carry homing endonuclease-like genes and are therefore assumed to propagate by "intron homing". A previous study demonstrated that three introns are located within the rRNA operon (arnSL) of Aeropyrum pernix strain K1, two of which, Ialpha and Igamma, harbor open reading frames (ORFs) encoding putative LAGLIDADG-type endonucleases. In an effort to understand further the rDNA intron distribution in natural A. pernix populations, 11 A. pernix strains were isolated from marine hydrothermal biotopes, and comparative nucleotide sequence analysis of the arnSL alleles was performed. Of the 11 isolates, eight contained multiple introns, and three patterns of intron insertion were found. Three novel introns, Idelta (62 bp in length), Ivarepsilon (122 bp) and Izeta (57 bp) were identified. They were all ORF-less, but their predicted RNA secondary structure at the exon-intron junctions was consistent with the bulge-helix-bulge motif. The insertion positions and the terminal inverted repeat sequences of Idelta and Izeta were in agreement with those of Ialpha and Igamma, respectively. This suggests that these intron variants were generated by large indels (insertions/deletions) during their evolution.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

A distinctive single-strand DNA-binding protein from the Archaeon Sulfolobus solfataricus

Single-stranded DNA binding proteins (SSBs) have been identified in all three domains of life. Here, we report the identification of a novel crenarchaeal SSB protein that is distinctly different from its euryarchaeal counterparts. Rather than comprising four DNA-binding domains and a zinc-finger motif within a single polypeptide of 645 amino acids, as for Methanococcus jannaschii, the Sulfolobus solfataricus SSB protein (SsoSSB) has a single DNA-binding domain in a polypeptide of just 148 amino acids with a eubacterial-like acidic C-terminus. SsoSSB protein was purified to homogeneity and found to form tetramers in solution, suggesting a quaternary structure analogous to that of E. coli SSB protein,despite possessing DNA-binding domains more similar to those of eukaryotic Replication Protein A (RPA). We demonstrate distributive binding of SsoSSB to ssDNA at high temperature with an apparent site size of approximately five nucleotides (nt)per monomer. Additionally, the protein is functional both in vitro and in vivo, stimulating RecA protein-mediated DNA strand-exchange and rescuing the ssb-1 lethal mutation of E. coli respectively. We discuss possible evolutionary relationships amongst the various members of the SSB/RPA family.

Molecular evolution and phylogenetic application of DMC1

The protein encoded by the single-copy nuclear gene DMC1 belongs to the recA-like group of proteins involved in meiosis. Partial nucleotide sequence, spanning exon 10 to exon 15, was used to test the applicability of the gene to phylogenetic studies in higher plants and used to assess its molecular evolution. The sequences produced from the Triticeae (Poaceae) show that most of the variation is confined to the introns. If a wider taxon sampling is used, alignment problems may be predicted. Comparisons including four complete coding sequences from GenBank reveal that the exons are more than twice as variable as rbcL, but easy to align, and hence may be valuable at higher taxonomic levels. Substitution rates are variable within the Triticeae, though local subclades show rate constancy. The relationships between exon variation and predicted protein structure are briefly discussed. In general, none of the observed nucleotide substitutions can be predicted to cause major structural or functional changes.

RadA protein from Archaeoglobus fulgidus forms rings, nucleoprotein filaments and catalyses homologous recombination

Proteins that catalyse homologous recombination have been identified in all living organisms and are essential for the repair of damaged DNA as well as for the generation of genetic diversity. In bacteria homologous recombination is performed by the RecA protein, whereas in the eukarya a related protein called Rad51 is required to catalyse recombination and repair. More recently, archaeal homologues of RecA/Rad51 (RadA) have been identified and isolated. In this work we have cloned and purified the RadA protein from the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus and characterised its in vitro activities. We show that (i) RadA protein forms ring structures in solution and binds single- but not double-stranded DNA to form nucleoprotein filaments, (ii) RadA is a single-stranded DNA-dependent ATPase at elevated temperatures, and (iii) RadA catalyses efficient D-loop formation and strand exchange at temperatures of 60-70 degrees C. Finally, we have used electron microscopy to visualise RadA-mediated joint molecules, the intermediates of homologous recombination. Intriguingly, RadA shares properties of both the bacterial RecA and eukaryotic Rad51 recombinases.

Domain analysis of an archaeal RadA protein for the strand exchange activity

Archaeal RadA, like eukaryotic Rad51 and bacterial RecA, promotes strand exchange between DNA strands with homologous sequences in vitro and is believed to participate in the homologous recombination in cells. The amino acid sequences of the archaeal RadA proteins are more similar to the eukaryotic Rad51s rather than the bacterial RecAs, and the N-terminal region containing domain I is conserved among the RadA and Rad51 proteins but is absent from RecA. To understand the structure-function relationship of RadA, we divided the RadA protein from Pyrococcus furiosus into two parts, the N-terminal one-third (RadA-n) and the residual C-terminal two-thirds (RadA-c), the latter of which contains the central core domain (domain II) of the RecA/Rad51 family proteins. RadA-c had the DNA-dependent ATPase activity and the strand exchange activity by itself, although much weaker (10%) than that of the intact RadA. These activities of RadA-c were restored to 60% of those of RadA by addition of RadA-n, indicating that the proper active structure of RadA was reconstituted in vitro. These results suggest that the basic activities of the RecA/Rad51 family proteins for homologous recombination are derived from domain II, and the N-terminal region may help to enhance the catalytic efficiencies.

Both RadA and RadB are involved in homologous recombination in Pyrococcus furiosus

RecA and Rad51 proteins are essential for homologous recombination in Bacteria and Eukarya, respectively. Homologous proteins, called RadA, have been described for Archaea. Here we present the characterization of two RecA/Rad51 family proteins, RadA and RadB, from Pyrococcus furiosus. The radA and radB genes were not induced by DNA damage resulting from exposure of the cells to gamma and UV irradiation and heat shock, suggesting that they might be constitutively expressed in this hyperthermophile. RadA had DNA-dependent ATPase, D-loop formation, and strand exchange activities. In contrast, RadB had a very weak ATPase activity that is not stimulated by DNA. This protein had a strong binding affinity for DNA, but little strand exchange activity could be detected. A direct interaction between RadA and RadB was detected by an immunoprecipitation assay. Moreover, RadB, but not RadA, coprecipitated with Hjc, a Holliday junction resolvase found in P. furiosus, in the absence of ATP. This interaction was suppressed in the presence of ATP. The Holliday junction cleavage activity of Hjc was inhibited by RadB in the absence, but not in the presence, of ATP. These results suggest that RadB has important roles in homologous recombination in Archaea and may regulate the cleavage reactions of the branch-structured DNA.

The RadA protein from a hyperthermophilic archaeon Pyrobaculum islandicum is a DNA-dependent ATPase that exhibits two disparate catalytic modes, with a transition temperature at 75 degrees C

The radA gene is an archaeal homolog of bacterial recA and eukaryotic RAD51 genes, which are critical components in homologous recombination and recombinational DNA repair. We cloned the radA gene from a hyperthermophilic archaeon, Pyrobaculum islandicum, overproduced the radA gene product in Escherichia coli and purified it to homogeneity. The purified P. islandicum RadA protein maintained its secondary structure and activities in vitro at high temperatures, up to 87 degrees C. It also showed high stability of 18.3 kcal.mol-1 (76.5 kJ.mol-1) at 25 degrees C and neutral pH. P. islandicum RadA exhibited activities typical of the family of RecA-like proteins, such as the ability to bind ssDNA, to hydrolyze ATP in a DNA-dependent manner and to catalyze DNA strand exchange. At 75 degrees C, all DNAs tested stimulated ATPase activity of the RadA. The protein exhibited a break in the Arrhenius plot of ATP hydrolysis at 75 degrees C. The cooperativity of ATP hydrolysis and ssDNA-binding ability of the protein above 75 degrees C were higher than at lower temperatures, and the activation energy of ATP hydrolysis was lower above this break point temperature. These results suggest that the ssDNA-dependent ATPase activity of P. islandicum RadA displays a temperature-dependent capacity to exist in two different catalytic modes, with 75 degrees C being the critical threshold temperature.

Efficient strand transfer by the RadA recombinase from the hyperthermophilic archaeon Desulfurococcus amylolyticus

The radA gene predicted to be responsible for homologous recombination in a hyperthermophilic archaeon, Desulfurococcus amylolyticus, was cloned, sequenced, and overexpressed in Escherichia coli cells. The deduced amino acid sequence of the gene product, RadA, was more similar to the human Rad51 protein (65% homology) than to the E. coli RecA protein (35%). A highly purified RadA protein was shown to exclusively catalyze single-stranded DNA-dependent ATP hydrolysis, which monitored presynaptic recombinational complex formation, at temperatures above 65 degrees C (catalytic rate constant of 1.2 to 2.5 min(-1) at 80 to 95 degrees C). The RadA protein alone efficiently promoted the strand exchange reaction at the range of temperatures from 80 to 90 degrees C, i.e., at temperatures approaching the melting point of DNA. It is noteworthy that both ATP hydrolysis and strand exchange are very efficient at temperatures optimal for host cell growth (90 to 92 degrees C).

Whole genome-based phylogenetic analysis of free-living microorganisms

A phylogenetic 'tree of life' has been constructed based on the observed presence and absence of families of protein-encoding genes observed in 11 complete genomes of free-living microorganisms. Past attempts to reconstruct the evolutionary relation-ships of microorganisms have been limited to sets of genes rather than complete genomes. Despite apparent rampant lateral gene transfer among microorganisms, these results indicate a single robust underlying evolutionary history for these organisms. Broadly, the tree produced is very similar to the small subunit rRNA tree although several additional phylogenetic relationships appear to be resolved, including the relationship of Archaeoglobus to the methanogens studied. This result is in contrast to notions that a robust phylogenetic reconstruction of microorganisms is impossible due to their genomes being composed of an incomprehensible amalgam of genes with complicated histories and suggests that this style of genome-wide phylogenetic analysis could become an important method for studying the ancient diversification of life on Earth. Analyses using informational and operational subsets of the genes showed that this 'tree of life' is not dependent on the phylogenetically more consistent informational genes.