Genomic analysis of the uncultivated marine crenarchaeote Cenarchaeum symbiosum

Assembling the marine metagenome, one cell at a time

The difficulty associated with the cultivation of most microorganisms and the complexity of natural microbial assemblages, such as marine plankton or human microbiome, hinder genome reconstruction of representative taxa using cultivation or metagenomic approaches. Here we used an alternative, single cell sequencing approach to obtain high-quality genome assemblies of two uncultured, numerically significant marine microorganisms. We employed fluorescence-activated cell sorting and multiple displacement amplification to obtain hundreds of micrograms of genomic DNA from individual, uncultured cells of two marine flavobacteria from the Gulf of Maine that were phylogenetically distant from existing cultured strains. Shotgun sequencing and genome finishing yielded 1.9 Mbp in 17 contigs and 1.5 Mbp in 21 contigs for the two flavobacteria, with estimated genome recoveries of about 91% and 78%, respectively. Only 0.24% of the assembling sequences were contaminants and were removed from further analysis using rigorous quality control. In contrast to all cultured strains of marine flavobacteria, the two single cell genomes were excellent Global Ocean Sampling (GOS) metagenome fragment recruiters, demonstrating their numerical significance in the ocean. The geographic distribution of GOS recruits along the Northwest Atlantic coast coincided with ocean surface currents. Metabolic reconstruction indicated diverse potential energy sources, including biopolymer degradation, proteorhodopsin photometabolism, and hydrogen oxidation. Compared to cultured relatives, the two uncultured flavobacteria have small genome sizes, few non-coding nucleotides, and few paralogous genes, suggesting adaptations to narrow ecological niches. These features may have contributed to the abundance of the two taxa in specific regions of the ocean, and may have hindered their cultivation. We demonstrate the power of single cell DNA sequencing to generate reference genomes of uncultured taxa from a complex microbial community of marine bacterioplankton. A combination of single cell genomics and metagenomics enabled us to analyze the genome content, metabolic adaptations, and biogeography of these taxa.

Environmental factors shaping the ecological niches of ammonia-oxidizing archaea

For more than 100 years it was believed that bacteria were the only group responsible for the oxidation of ammonia. However, recently, a new strain of archaea bearing a putative ammonia monooxygenase subunit A (amoA) gene and able to oxidize ammonia was isolated from a marine aquarium tank. Ammonia-oxidizing archaea (AOA) were subsequently discovered in many ecosystems of varied characteristics and even found as the predominant causal organisms in some environments. Here, we summarize the current knowledge on the environmental conditions related to the presence of AOA and discuss the possible site-related properties. Considering these data, we deduct the possible niches of AOA based on pH, sulfide and phosphate levels. It is proposed that the AOA might be important actors within the nitrogen cycle in low-nutrient, low-pH, and sulfide-containing environments.

Unique archaeal assemblages in the Arctic Ocean unveiled by massively parallel tag sequencing

The Arctic Ocean plays a critical role in controlling nutrient budgets between the Pacific and Atlantic Ocean. Archaea are key players in the nitrogen cycle and in cycling nutrients, but their community composition has been little studied in the Arctic Ocean. Here, we characterize archaeal assemblages from surface and deep Arctic water masses using massively parallel tag sequencing of the V6 region of the 16S rRNA gene. This approach gave a very high coverage of the natural communities, allowing a precise description of archaeal assemblages. This first taxonomic description of archaeal communities by tag sequencing reported so far shows that it is possible to assign an identity below phylum level to most (95%) of the archaeal V6 tags, and shows that tag sequencing is a powerful tool for resolving the diversity and distribution of specific microbes in the environment. Marine group I Crenarchaeota was overall the most abundant group in the Arctic Ocean and comprised between 27% and 63% of all tags. Group III Euryarchaeota were more abundant in deep-water masses and represented the largest archaeal group in the deep Atlantic layer of the central Arctic Ocean. Coastal surface waters, in turn, harbored more group II Euryarchaeota. Moreover, group II sequences that dominated surface waters were different from the group II sequences detected in deep waters, suggesting functional differences in closely related groups. Our results unveiled for the first time an archaeal community dominated by group III Euryarchaeota and show biogeographical traits for marine Arctic Archaea.

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(2) assimilation by Bacteria rather than Archaea. Our results indicate that Archaea were not important for ammonia oxidation in the agricultural soil tested.

A bioinformatician's guide to metagenomics

As random shotgun metagenomic projects proliferate and become the dominant source of publicly available sequence data, procedures for the best practices in their execution and analysis become increasingly important. Based on our experience at the Joint Genome Institute, we describe the chain of decisions accompanying a metagenomic project from the viewpoint of the bioinformatic analysis step by step. We guide the reader through a standard workflow for a metagenomic project beginning with presequencing considerations such as community composition and sequence data type that will greatly influence downstream analyses. We proceed with recommendations for sampling and data generation including sample and metadata collection, community profiling, construction of shotgun libraries, and sequencing strategies. We then discuss the application of generic sequence processing steps (read preprocessing, assembly, and gene prediction and annotation) to metagenomic data sets in contrast to genome projects. Different types of data analyses particular to metagenomes are then presented, including binning, dominant population analysis, and gene-centric analysis. Finally, data management issues are presented and discussed. We hope that this review will assist bioinformaticians and biologists in making better-informed decisions on their journey during a metagenomic project.

Tri-split tRNA is a transfer RNA made from 3 transcripts that provides insight into the evolution of fragmented tRNAs in archaea

Transfer RNA (tRNA) is essential for decoding the genome sequence into proteins. In Archaea, previous studies have revealed unique multiple intron-containing tRNAs and tRNAs that are encoded on 2 separate genes, so-called split tRNAs. Here, we discovered 10 fragmented tRNA genes in the complete genome of the hyperthermoacidophilic Archaeon Caldivirga maquilingensis that are individually transcribed and further trans-spliced to generate all of the missing tRNAs encoding glycine, alanine, and glutamate. Notably, the 3 mature tRNA(Gly)'s with synonymous codons are created from 1 constitutive 3' half transcript and 4 alternatively switching transcripts, representing tRNA made from a total of 3 transcripts named a "tri-split tRNA." Expression and nucleotide sequences of 10 split tRNA genes and their joined tRNA products were experimentally verified. The intervening sequences of split tRNA have high identity to tRNA intron sequences located at the same positions in intron-containing tRNAs in related Thermoproteales species. This suggests that an evolutionary relationship between intron-containing and split tRNAs exists. Our findings demonstrate the first example of split tRNA genes in a free-living organism and a unique tri-split tRNA gene that provides further insight into the evolution of fragmented tRNAs.

Vertical structure of archaeal communities and the distribution of ammonia monooxygenase A gene variants in two meromictic High Arctic lakes

The distribution of archaeal amoA and 16S rRNA genes was evaluated in two marine-derived, meromictic lakes in the Canadian High Arctic: Lake A and Lake C1 on the northern coast of Ellesmere Island. The amoA gene was recorded in both lakes, with highest copy numbers in the oxycline. Sequence analysis showed that amoA from the two lakes shared 94% similarity, indicating at least two phylogenetically distinct clusters. Clone libraries of archaeal 16S rRNA genes from Lake A revealed strong vertical differences in archaeal community diversity and composition down the water column. The oxic layer was dominated by one group of Euryarchaeota affiliated to the Lake Dagow Sediment (LDS) cluster. This group was absent from the oxycline, which had an extremely low archaeal diversity of two phylotypes. Both belonged to the Crenarchaeota Marine Group I (MGI), the marine group that has been linked to archaeal amoA; however, there was a low ratio of amoA to MGI copy numbers, suggesting that many MGI Archaea did not carry the amoA gene. The anoxic zone contained representatives of the RC-V (Rice Cluster-V) and LDS clusters of Euryarchaeota. These results show the strong vertical differentiation of archaeal communities in polar meromictic lakes, and they suggest archaeal nitrification within the oxycline of these highly stratified waters.

A DNA topoisomerase IB in Thaumarchaeota testifies for the presence of this enzyme in the last common ancestor of Archaea and Eucarya

DNA topoisomerase IB (TopoIB) was thought for a long time to be a eukaryotic specific enzyme. A shorter version was then found in viruses and later on in several bacteria, but not in archaea. Here, we show that a eukaryotic-like TopoIB is present in the recently sequenced genomes of two archaea of the newly proposed phylum Thaumarchaeota. Phylogenetic analyses suggest that a TopoIB was present in the last common ancestor of Archaea and Eucarya. This finding indicates that the last common ancestor of Archaea and Eucarya may have harboured a DNA genome.

Environmental detection of octahaem cytochrome c hydroxylamine/hydrazine oxidoreductase genes of aerobic and anaerobic ammonium-oxidizing bacteria

Bacterial aerobic ammonium oxidation and anaerobic ammonium oxidation (anammox) are important processes in the global nitrogen cycle. Key enzymes in both processes are the octahaem cytochrome c (OCC) proteins, hydroxylamine oxidoreductase (HAO) of aerobic ammonium-oxidizing bacteria (AOB), which catalyses the oxidation of hydroxylamine to nitrite, and hydrazine oxidoreductase (HZO) of anammox bacteria, which converts hydrazine to N(2). While the genomes of AOB encode up to three nearly identical copies of hao operons, genome analysis of Candidatus'Kuenenia stuttgartiensis' showed eight highly divergent octahaem protein coding regions as possible candidates for the HZO. Based on their phylogenetic relationship and biochemical characteristics, the sequences of these eight gene products grouped in three clusters. Degenerate primers were designed on the basis of available gene sequences with the aim to detect hao and hzo genes in various ecosystems. The hao primer pairs amplified gene fragments from 738 to 1172 bp and the hzo primer pairs amplified gene fragments from 289 to 876 bp in length, when tested on genomic DNA isolated from a variety of AOB and anammox bacteria. A selection of these primer pairs was also used successfully to amplify and analyse the hao and hzo genes in community DNA isolated from different ecosystems harbouring both AOB and anammox bacteria. We propose that OCC protein-encoding genes are suitable targets for molecular ecological studies on both aerobic and anaerobic ammonium-oxidizing bacteria.

The archaebacterial origin of eukaryotes

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

Major gradients in putatively nitrifying and non-nitrifying Archaea in the deep North Atlantic

Aerobic nitrification of ammonia to nitrite and nitrate is a key process in the oceanic nitrogen cycling mediated by prokaryotes. Apart from Bacteria belonging to the beta- and gamma-Proteobacteria involved in the first nitrification step, Crenarchaeota have recently been recognized as main drivers of the oxidation of ammonia to nitrite in soil as well as in the ocean, as indicated by the dominance of archaeal ammonia monooxygenase (amoA) genes over bacterial amoA. Evidence is accumulating that archaeal amoA genes are common in a wide range of marine systems. Essentially, all these reports focused on surface and mesopelagic (200-1,000 m depth) waters, where ammonia concentrations are higher than in waters below 1,000 m depth. However, Crenarchaeota are also abundant in the water column below 1,000 m, where ammonia concentrations are extremely low. Here we show that, throughout the North Atlantic Ocean, the abundance of archaeal amoA genes decreases markedly from subsurface waters to 4,000 m depth, and from subpolar to equatorial deep waters, leading to pronounced vertical and latitudinal gradients in the ratio of archaeal amoA to crenarchaeal 16S ribosomal RNA (rRNA) genes. The lack of significant copy numbers of amoA genes and the very low fixation rates of dark carbon dioxide in the bathypelagic North Atlantic suggest that most bathypelagic Crenarchaeota are not autotrophic ammonia oxidizers: most likely, they utilize organic matter and hence live heterotrophically.

The dynamic genetic repertoire of microbial communities

Community genomic data have revealed multiple levels of variation between and within microbial consortia. This variation includes large-scale differences in gene content between ecosystems as well as within-population sequence heterogeneity. In the present review, we focus specifically on how fine-scale variation within microbial and viral populations is apparent from community genomic data. A major unresolved question is how much of the observed variation is due to neutral vs. adaptive processes. Limited experimental data hint that some of this fine-scale variation may be in part functionally relevant, whereas sequence-based and modeling analyses suggest that much of it may be neutral. While methods for interpreting population genomic data are still in their infancy, we discuss current interpretations of existing datasets in the light of evolutionary processes and models. Finally, we highlight the importance of virus–host dynamics in generating and shaping within-population diversity.

Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea

Ammonia-oxidizing archaea (AOA) have recently been found to be potentially important in nitrogen cycling in a variety of environments, such as terrestrial soils, wastewater treatment reactors, marine waters and sediments, and especially in estuaries, where high input of anthropogenic nitrogen is often experienced. The sedimentary AOA diversity, community structure and spatial distribution in the Changjiang Estuary and the adjacent East China Sea were studied. Multivariate statistical analysis indicated that the archaeal amoA genotype communities could be clustered according to sampling transects, and the station located in an estuarine mixing zone harboured a distinct AOA community. The distribution of AOA communities correlated significantly with the gradients of surface-water salinity and sediment sorting coefficient. The spatial distribution of putative soil-related AOA in certain sampling stations indicated a strong impact of the Changjiang freshwater discharge on the marine benthic microbial ecosystem. Besides freshwater, nutrients, organic matter and suspended particles, the Changjiang Diluted Water might also contribute to the transport of terrestrial archaea into the seawater and sediments along its flow path.

Ammonia-oxidizing Crenarchaeota and nitrification inside the tissue of a colonial ascidian

Marine Crenarchaeota represent an abundant component of the oceanic microbiota that play an important role in the global nitrogen cycle. Here we report the association of the colonial ascidian Cystodytes dellechiajei with putative ammonia-oxidizing Crenarchaeota that could actively be involved in nitrification inside the animal tissue. As shown by 16S rRNA gene analysis, the ascidian-associated Crenarchaeota were phylogenetically related to Nitrosopumilus maritimus, the first marine archaeon isolated in pure culture that grows chemolithoautotrophically oxidizing ammonia to nitrite aerobically. Catalysed reporter deposition (CARD)-FISH revealed that the Crenarchaeota were specifically located inside the tunic tissue of the colony, where moreover the expression of amoA gene was detected. The amoA gene encodes the alpha-subunit of ammonia monooxygenase, which is involved in the first step of nitrification, the oxidation of ammonia to nitrite. Sequencing of amoA gene showed that they were phylogenetically related to amoA genes of N. maritimus and other putative ammonia-oxidizing marine Crenarchaeota. In order to track the suspected nitrification activity inside the ascidian colony under in vivo conditions, microsensor profiles were measured through the tunic tissue. Net NO(x) production was detected in the tunic layer 1200-1800 microm with rates of 58-90 nmol cm(-3) h(-1). Oxygen and pH microsensor profiles showed that the layer of net NO(x) production coincided with O(2) concentrations of 103-116 microM and pH value of 5.2. Together, molecular and microsensor data indicate that Crenarchaeota could oxidize ammonia to nitrite aerobically, and thus be involved in nitrification inside the ascidian tissue.

Crenarchaeota and their role in the nitrogen cycle in a subsurface radioactive thermal spring in the Austrian Central Alps

Previous results from a 16S rRNA gene library analysis showed high diversity within the prokaryotic community of a subterranean radioactive thermal spring, the "Franz-Josef-Quelle" (FJQ) in Bad Gastein, Austria, as well as evidence for ammonia oxidation by crenarchaeota. This study reports further characterization of the community by denaturing gradient gel electrophoresis (DGGE) analysis, fluorescence in situ hybridization (FISH), and semiquantitative nitrification measurements. DGGE bands from three types of samples (filtered water, biofilms on glass slides, and naturally grown biofilms), including samples collected at two distinct times (January 2005 and July 2006), were analyzed. The archaeal community consisted mainly of Crenarchaeota of the soil-subsurface-freshwater group (group 1.1b) and showed a higher diversity than in the previous 16S rRNA gene library analysis, as was also found for crenarchaeal amoA genes. No bacterial amoA genes were detected. FISH analysis of biofilms indicated the presence of archaeal cells with an abundance of 5.3% (+/-4.5%) in the total 4',6-diamidino-2-phenylindole (DAPI)-stained community. Microcosm experiments of several weeks in duration showed a decline of ammonium that correlated with an increase of nitrite, the presence of crenarchaeal amoA genes, and the absence of bacterial amoA genes. The data suggested that only ammonia-oxidizing archaea (AOA) perform the first step of nitrification in this 45 degrees C environment. The crenarchaeal amoA gene sequences grouped within a novel cluster of amoA sequences from the database, originating from geothermally influenced environments, for which we propose the designation "thermal spring" cluster and which may be older than most AOA from soils on earth.

Phylogenetic diversity and spatial distribution of the microbial community associated with the Caribbean deep-water sponge Polymastia cf. corticata by 16S rRNA, aprA, and amoA gene analysis

Denaturing gradient gel electrophoresis (DGGE)-based analyses of 16S rRNA, aprA, and amoA genes demonstrated that a phylogenetically diverse and complex microbial community was associated with the Caribbean deep-water sponge Polymastia cf. corticata Ridley and Dendy, 1887. From the 38 archaeal and bacterial 16S rRNA phylotypes identified, 53% branched into the sponge-specific, monophyletic sequence clusters determined by previous studies (considering predominantly shallow-water sponge species), whereas 26% appeared to be P. cf. corticata specifically associated microorganisms ("specialists"); 21% of the phylotypes were confirmed to represent seawater- and sediment-derived proteobacterial species ("contaminants") acquired by filtration processes from the host environment. Consistently, the aprA and amoA gene-based analyses indicated the presence of environmentally derived sulfur- and ammonia-oxidizers besides putative sponge-specific sulfur-oxidizing Gammaproteobacteria and Alphaproteobacteria and a sulfate-reducing archaeon. A sponge-specific, endosymbiotic sulfur cycle as described for marine oligochaetes is proposed to be also present in P. cf. corticata. Overall, the results of this work support the recent studies that demonstrated the sponge species specificity of the associated microbial community while the biogeography of the host collection site has only a minor influence on the composition. In P. cf. corticata, the specificity of the sponge-microbe associations is even extended to the spatial distribution of the microorganisms within the sponge body; distinct bacterial populations were associated with the different tissue sections, papillae, outer and inner cortex, and choanosome. The local distribution of a phylotype within P. cf. corticata correlated with its (1) phylogenetic affiliation, (2) classification as sponge-specific or nonspecifically associated microorganism, and (3) potential ecological role in the host sponge.

Genomic patterns of recombination, clonal divergence and environment in marine microbial populations

Microorganisms represent the largest reservoir of biodiversity on Earth, both in numbers and total genetic diversity, but it remains unclear whether this biodiversity is organized in discrete units that correspond to ecologically coherent species. To further explore this question, we examined patterns of genomic diversity in sympatric microbial populations. Analyses of a total of approximately 200 Mb of microbial community genomic DNA sequence recovered from 4000 m depth in the Pacific Ocean revealed discrete sequence-defined populations of Bacteria and Archaea, with intrapopulation genomic sequence divergence ranging from approximately 1% to approximately 6%. The populations appeared to be maintained, at least in part, by intrapopulation genetic exchange (homologous recombination), although the frequency of recombination was estimated to be about three times lower than that observed previously in thermoacidophilic archaeal biofilm populations. Furthermore, the genotypes of a given population were clearly distinguishable from their closest co-occurring relatives based on their relative abundance in situ. The genetic distinctiveness and the matching sympatric abundances imply that these genotypes share similar ecophysiological properties, and therefore may represent fundamental units of microbial diversity in the deep sea. Comparisons to surface-dwelling relatives of the Sargasso Sea revealed that distinct sequence-based clusters were not always detectable, presumably due to environmental variations, further underscoring the important relationship between environmental contexts and genetic mechanisms, which together shape and sustain microbial population structure.

Genomic sequencing of single microbial cells from environmental samples

Recently developed techniques allow genomic DNA sequencing from single microbial cells [Lasken RS: Single-cell genomic sequencing using multiple displacement amplification. Curr Opin Microbiol 2007, 10:510-516]. Here, we focus on research strategies for putting these methods into practice in the laboratory setting. An immediate consequence of single-cell sequencing is that it provides an alternative to culturing organisms as a prerequisite for genomic sequencing. The microgram amounts of DNA required as template are amplified from a single bacterium by a method called multiple displacement amplification (MDA) avoiding the need to grow cells. The ability to sequence DNA from individual cells will likely have an immense impact on microbiology considering the vast numbers of novel organisms, which have been inaccessible unless culture-independent methods could be used. However, special approaches have been necessary to work with amplified DNA. MDA may not recover the entire genome from the single copy present in most bacteria. Also, some sequence rearrangements can occur during the DNA amplification reaction. Over the past two years many research groups have begun to use MDA, and some practical approaches to single-cell sequencing have been developed. We review the consensus that is emerging on optimum methods, reliability of amplified template, and the proper interpretation of 'composite' genomes which result from the necessity of combining data from several single-cell MDA reactions in order to complete the assembly. Preferred laboratory methods are considered on the basis of experience at several large sequencing centers where >70% of genomes are now often recovered from single cells. Methods are reviewed for preparation of bacterial fractions from environmental samples, single-cell isolation, DNA amplification by MDA, and DNA sequencing.

A korarchaeal genome reveals insights into the evolution of the Archaea

The candidate division Korarchaeota comprises a group of uncultivated microorganisms that, by their small subunit rRNA phylogeny, may have diverged early from the major archaeal phyla Crenarchaeota and Euryarchaeota. Here, we report the initial characterization of a member of the Korarchaeota with the proposed name, "Candidatus Korarchaeum cryptofilum," which exhibits an ultrathin filamentous morphology. To investigate possible ancestral relationships between deep-branching Korarchaeota and other phyla, we used whole-genome shotgun sequencing to construct a complete composite korarchaeal genome from enriched cells. The genome was assembled into a single contig 1.59 Mb in length with a G + C content of 49%. Of the 1,617 predicted protein-coding genes, 1,382 (85%) could be assigned to a revised set of archaeal Clusters of Orthologous Groups (COGs). The predicted gene functions suggest that the organism relies on a simple mode of peptide fermentation for carbon and energy and lacks the ability to synthesize de novo purines, CoA, and several other cofactors. Phylogenetic analyses based on conserved single genes and concatenated protein sequences positioned the korarchaeote as a deep archaeal lineage with an apparent affinity to the Crenarchaeota. However, the predicted gene content revealed that several conserved cellular systems, such as cell division, DNA replication, and tRNA maturation, resemble the counterparts in the Euryarchaeota. In light of the known composition of archaeal genomes, the Korarchaeota might have retained a set of cellular features that represents the ancestral archaeal form.

Phylogenetic analyses of ribosomal DNA-containing bacterioplankton genome fragments from a 4000 m vertical profile in the North Pacific Subtropical Gyre

High-throughput identification of rRNA gene-containing clones in large insert metagenomic libraries is difficult, because of the high background of host ribosomal RNA (rRNA) and rRNA genes. To address this challenge, a membrane hybridization method was developed to identify all bacterial small subunit rRNA-containing fosmid clones of microbial community DNA from seven different depths in the North Pacific Subtropical Gyre. Out of 101,376 clones screened, 751 rDNA-containing clones were identified that grouped in approximately 60 different clades. Several rare sequences only remotely related to known groups were detected, including a Wolbachia-related sequence containing a putative intron or intervening sequence, as well as seven sequences from Order Myxococcales not previously detected in pelagic habitats. Stratified, depth-specific population structure was evident within both cultured and uncultured lineages. Conversely, some eurybathyal members of the genera Alcanivorax and Rhizobium shared identical small subunit ribosomal DNA sequences that were distributed from surface waters to the 4000 m depth. Comparison with similar analyses in Monterey Bay microbial communities revealed previously recognized, as well as some distinctive, depth-stratified partitioning that distinguished coastal from open ocean bacterioplankton populations. While some bias was evident in fosmid clone recovery in a few particular lineages, the overall phylogenetic group recovery and distributions were consistent with previous studies, as well as with direct shotgun sequence data from the same source DNA.

Microbial diversity and the genetic nature of microbial species

The earth contains a huge number of largely uncharacterized Bacteria and Archaea. Microbiologists are struggling to summarize their genetic diversity and classify them, which has resulted in heated debates on methods for defining species, mechanisms that lead to speciation and whether microbial species even exist. This Review proposes that decisions on the existence of species and methods to define them should be guided by a method-free species concept that is based on cohesive evolutionary forces. It summarizes current approaches to defining species and the problems of these approaches, and presents selected examples of the population genetic patterns at and below the species level.

Nitrification in terrestrial hot springs of Iceland and Kamchatka

Archaea have been detected recently as a major and often dominant component of the microbial communities performing ammonia oxidation in terrestrial and marine environments. In a molecular survey of archaeal ammonia monooxygenase (AMO) genes in terrestrial hot springs of Iceland and Kamchatka, the amoA gene encoding the alpha-subunit of AMO was detected in a total of 14 hot springs out of the 22 investigated. Most of these amoA-positive hot springs had temperatures between 82 and 97 degrees C and pH range between 2.5 and 7. In phylogenetic analyses, these amoA genes formed three independent lineages within the known sequence clusters of marine or soil origin. Furthermore, in situ gross nitrification rates in Icelandic hot springs were estimated by the pool dilution technique directly on site. At temperatures above 80 degrees C, between 56 and 159 mumol NO(3)(-) L(-1) mud per day was produced. Furthermore, addition of ammonium to the hot spring samples before incubation yielded a more than twofold higher potential nitrification rate, indicating that the process was limited by ammonia supply. Our data provide evidence for an active role of archaea in nitrification of hot springs in a wide range of pH values and at a high temperature.

Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota

The archaeal domain is currently divided into two major phyla, the Euryarchaeota and Crenarchaeota. During the past few years, diverse groups of uncultivated mesophilic archaea have been discovered and affiliated with the Crenarchaeota. It was recently recognized that these archaea have a major role in geochemical cycles. Based on the first genome sequence of a crenarchaeote, Cenarchaeum symbiosum, we show that these mesophilic archaea are different from hyperthermophilic Crenarchaeota and branch deeper than was previously assumed. Our results indicate that C. symbiosum and its relatives are not Crenarchaeota, but should be considered as a third archaeal phylum, which we propose to name Thaumarchaeota (from the Greek 'thaumas', meaning wonder).

A moderately thermophilic ammonia-oxidizing crenarchaeote from a hot spring

The recent discovery of ammonia-oxidizing archaea (AOA) dramatically changed our perception of the diversity and evolutionary history of microbes involved in nitrification. In this study, a moderately thermophilic (46 degrees C) ammonia-oxidizing enrichment culture, which had been seeded with biomass from a hot spring, was screened for ammonia oxidizers. Although gene sequences for crenarchaeotal 16S rRNA and two subunits of the ammonia monooxygenase (amoA and amoB) were detected via PCR, no hints for known ammonia-oxidizing bacteria were obtained. Comparative sequence analyses of these gene fragments demonstrated the presence of a single operational taxonomic unit and thus enabled the assignment of the amoA and amoB sequences to the respective 16S rRNA phylotype, which belongs to the widely distributed group I.1b (soil group) of the Crenarchaeota. Catalyzed reporter deposition (CARD)-FISH combined with microautoradiography (MAR) demonstrated metabolic activity of this archaeon in the presence of ammonium. This finding was corroborated by the detection of amoA gene transcripts in the enrichment. CARD-FISH/MAR showed that the moderately thermophilic AOA is highly active at 0.14 and 0.79 mM ammonium and is partially inhibited by a concentration of 3.08 mM. The enriched AOA, which is provisionally classified as "Candidatus Nitrososphaera gargensis," is the first described thermophilic ammonia oxidizer and the first member of the crenarchaeotal group I.1b for which ammonium oxidation has been verified on a cellular level. Its preference for thermophilic conditions reinvigorates the debate on the thermophilic ancestry of AOA.

Genome sequence of Thermofilum pendens reveals an exceptional loss of biosynthetic pathways without genome reduction

We report the complete genome of Thermofilum pendens, a deeply branching, hyperthermophilic member of the order Thermoproteales in the archaeal kingdom Crenarchaeota. T. pendens is a sulfur-dependent, anaerobic heterotroph isolated from a solfatara in Iceland. It is an extracellular commensal, requiring an extract of Thermoproteus tenax for growth, and the genome sequence reveals that biosynthetic pathways for purines, most amino acids, and most cofactors are absent. In fact, T. pendens has fewer biosynthetic enzymes than obligate intracellular parasites, although it does not display other features that are common among obligate parasites and thus does not appear to be in the process of becoming a parasite. It appears that T. pendens has adapted to life in an environment rich in nutrients. T. pendens was known previously to utilize peptides as an energy source, but the genome revealed a substantial ability to grow on carbohydrates. T. pendens is the first crenarchaeote and only the second archaeon found to have a transporter of the phosphotransferase system. In addition to fermentation, T. pendens may obtain energy from sulfur reduction with hydrogen and formate as electron donors. It may also be capable of sulfur-independent growth on formate with formate hydrogen lyase. Additional novel features are the presence of a monomethylamine:corrinoid methyltransferase, the first time that this enzyme has been found outside the Methanosarcinales, and the presence of a presenilin-related protein. The predicted highly expressed proteins do not include proteins encoded by housekeeping genes and instead include ABC transporters for carbohydrates and peptides and clustered regularly interspaced short palindromic repeat-associated proteins.

"Candidatus Cloacamonas acidaminovorans": genome sequence reconstruction provides a first glimpse of a new bacterial division

Many microorganisms live in anaerobic environments. Most of these microorganisms have not yet been cultivated. Here, we present, from a metagenomic analysis of an anaerobic digester of a municipal wastewater treatment plant, a reconstruction of the complete genome of a bacterium belonging to the WWE1 candidate division. In silico proteome analysis indicated that this bacterium might derive most of its carbon and energy from the fermentation of amino acids, and hence, it was provisionally classified as "Candidatus Cloacamonas acidaminovorans." "Candidatus Cloacamonas acidaminovorans" is probably a syntrophic bacterium that is present in many anaerobic digesters. This report highlights how environmental sequence data might provide genomic and functional information about a new bacterial clade whose members are involved in anaerobic digestion.

New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation

Microbial activities drive the global nitrogen cycle, and in the past few years, our understanding of nitrogen cycling processes and the micro-organisms that mediate them has changed dramatically. During this time, the processes of anaerobic ammonium oxidation (anammox), and ammonia oxidation within the domain Archaea, have been recognized as two new links in the global nitrogen cycle. All available evidence indicates that these processes and organisms are critically important in the environment, and particularly in the ocean. Here we review what is currently known about the microbial ecology of anaerobic and archaeal ammonia oxidation, highlight relevant unknowns and discuss the implications of these discoveries for the global nitrogen and carbon cycles.

Molecular and biogeochemical evidence for ammonia oxidation by marine Crenarchaeota in the Gulf of California

Nitrification plays an important role in marine biogeochemistry, yet efforts to link this process to the microorganisms that mediate it are surprisingly limited. In particular, ammonia oxidation is the first and rate-limiting step of nitrification, yet ammonia oxidation rates and the abundance of ammonia-oxidizing bacteria (AOB) have rarely been measured in tandem. Ammonia oxidation rates have not been directly quantified in conjunction with ammonia-oxidizing archaea (AOA), although mounting evidence indicates that marine Crenarchaeota are capable of ammonia oxidation, and they are among the most abundant microbial groups in the ocean. Here, we have directly quantified ammonia oxidation rates by 15N labeling, and AOA and AOB abundances by quantitative PCR analysis of ammonia monooxygenase subunit A (amoA) genes, in the Gulf of California. Based on markedly different archaeal amoA sequence types in the upper water column (60 m) and oxygen minimum zone (OMZ; 450 m), novel amoA PCR primers were designed to specifically target and quantify 'shallow' (group A) and 'deep' (group B) clades. These primers recovered extensive variability with depth. Within the OMZ, AOA were most abundant where nitrification may be coupled to denitrification. In the upper water column, group A tracked variations in nitrogen biogeochemistry with depth and between basins, whereas AOB were present in relatively low numbers or undetectable. Overall, 15NH4+ oxidation rates were remarkably well correlated with AOA group A amoA gene copies (r2=0.90, P<0.001), and with 16S rRNA gene copies from marine Crenarchaeota (r2=0.85, P<0.005). These findings represent compelling evidence for an archaeal role in oceanic nitrification.

High intraspecific recombination rate in a native population of Candidatus pelagibacter ubique (SAR11)

Recombination is an important process in microbial evolution. Rates of recombination with extracellular DNA matter because models of microbial population structure are profoundly influenced by the degree to which recombination is occurring within the population. Low rates of recombination may be sufficient to ensure the lateral propagation of genes that have a high selective advantage without disrupting the clonal pattern of inheritance for other genes. High rates of recombination potentially can obscure clonal patterns, leading to linkage equilibrium, and give microbial populations a population genetic structure more akin to sexually interbreeding eukaryotic populations. We examined eight loci from nine strains of candidatus Pelagibacter ubique (SAR11), isolated from a single 2L niskin sample of natural seawater, for evidence of genetic recombination between strains. The Shimodaira-Hasegawa test revealed significant phylogenetic incongruence in seven of the genes, indicating that frequent recombination obscures phylogenetic signals from the linear inheritance of genes in this population. Statistical evidence for intragenic recombination was found for six loci. An informative sites matrix showed extensive evidence for a widespread breakdown of linkage disequilibrium. Although the mechanisms of genetic transfer in native SAR11 populations are unknown, we measured recombination rates, rho, that are much higher than point mutation rates, theta, as a source of genetic diversity in this clade. The eukaryotic model of species sharing a common pool of alleles is more apt for this SAR11 population than a strictly clonal model of inheritance in which allelic diversity is controlled by periodic selection.

Diversity and mode of transmission of ammonia-oxidizing archaea in marine sponges

The model marine crenarchaeote 'Cenarchaeum symbiosum' is until now the only ammonia-oxidizing archaeon known from a marine sponge. Here, phylogenetic analyses based on the 16S rRNA and ammonia monooxygenase subunit A (amoA) genes revealed the presence of putative ammonia-oxidizing archaea (AOA) in a diverse range of sponges from the western Pacific, Caribbean and Mediterranean. amoA diversity was limited even between different oceans, with many of the obtained sequences (75.9%; n(total) = 83) forming a monophyletic, apparently sponge- (and coral-) specific lineage, analogous to those previously inferred from comparative 16S rRNA gene studies of sponge-associated microbes. The presence of AOA in sponge larvae, as detected by 16S rRNA and amoA PCR assays as well as by fluorescence in situ hybridization, suggests they are vertically transmitted and thus might be of importance for ammonia detoxification within the sponge.

Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea

Background

An evolutionary classification of genes from sequenced genomes that distinguishes between orthologs and paralogs is indispensable for genome annotation and evolutionary reconstruction. Shortly after multiple genome sequences of bacteria, archaea, and unicellular eukaryotes became available, an attempt on such a classification was implemented in Clusters of Orthologous Groups of proteins (COGs). Rapid accumulation of genome sequences creates opportunities for refining COGs but also represents a challenge because of error amplification. One of the practical strategies involves construction of refined COGs for phylogenetically compact subsets of genomes.

Results

New Archaeal Clusters of Orthologous Genes (arCOGs) were constructed for 41 archaeal genomes (13 Crenarchaeota, 27 Euryarchaeota and one Nanoarchaeon) using an improved procedure that employs a similarity tree between smaller, group-specific clusters, semi-automatically partitions orthology domains in multidomain proteins, and uses profile searches for identification of remote orthologs. The annotation of arCOGs is a consensus between three assignments based on the COGs, the CDD database, and the annotations of homologs in the NR database. The 7538 arCOGs, on average, cover ~88% of the genes in a genome compared to a ~76% coverage in COGs. The finer granularity of ortholog identification in the arCOGs is apparent from the fact that 4538 arCOGs correspond to 2362 COGs; ~40% of the arCOGs are new. The archaeal gene core (protein-coding genes found in all 41 genome) consists of 166 arCOGs. The arCOGs were used to reconstruct gene loss and gene gain events during archaeal evolution and gene sets of ancestral forms. The Last Archaeal Common Ancestor (LACA) is conservatively estimated to possess 996 genes compared to 1245 and 1335 genes for the last common ancestors of Crenarchaeota and Euryarchaeota, respectively. It is inferred that LACA was a chemoautotrophic hyperthermophile that, in addition to the core archaeal functions, encoded more idiosyncratic systems, e.g., the CASS systems of antivirus defense and some toxin-antitoxin systems.

Conclusion

The arCOGs provide a convenient, flexible framework for functional annotation of archaeal genomes, comparative genomics and evolutionary reconstructions. Genomic reconstructions suggest that the last common ancestor of archaea might have been (nearly) as advanced as the modern archaeal hyperthermophiles. ArCOGs and related information are available at: .

Reviewers

This article was reviewed by Peer Bork, Patrick Forterre, and Purificacion Lopez-Garcia.

Uncultured archaea in deep marine subsurface sediments: have we caught them all?

Deep marine subsurface sediments represent a novel archaeal biosphere with unknown physiology; the sedimentary subsurface harbors numerous novel phylogenetic lineages of archaea that are at present uncultured. Archaeal 16S rRNA analyses of deep subsurface sediments demonstrate their global occurrence and wide habitat range, including deep subsurface sediments, methane seeps and organic-rich coastal sediments. These subsurface archaeal lineages were discovered by PCR of extracted environmental DNA; their detection ultimately depends on the specificity of the archaeal PCR 16S rRNA primers. Surprisingly high mismatch frequencies for some archaeal PCR primers result in amplification bias against the corresponding archaeal lineages; this review presents some examples. Obviously, most archaeal 16S rRNA PCR primers were developed either before the discovery of these deep subsurface archaeal lineages, or without taking their sequence variants into account. PCR surveys with multiple primer combinations, revision and updates of primers whenever possible, and increasing use of PCR-independent methods in molecular microbial ecology will contribute to a more comprehensive view of subsurface archaeal communities.

Natural variation in SAR11 marine bacterioplankton genomes inferred from metagenomic data

BACKGROUND

One objective of metagenomics is to reconstruct information about specific uncultured organisms from fragmentary environmental DNA sequences. We used the genome of an isolate of the marine alphaproteobacterium SAR11 ('Candidatus Pelagibacter ubique'; strain HTCC1062), obtained from the cold, productive Oregon coast, as a query sequence to study variation in SAR11 metagenome sequence data from the Sargasso Sea, a warm, oligotrophic ocean gyre.

RESULTS

The average amino acid identity of SAR11 genes encoded by the metagenomic data to the query genome was only 71%, indicating significant evolutionary divergence between the coastal isolates and Sargasso Sea populations. However, an analysis of gene neighbors indicated that SAR11 genes in the Sargasso Sea metagenomic data match the gene order of the HTCC1062 genome in 96% of cases (> 85,000 observations), and that rearrangements are most frequent at predicted operon boundaries. There were no conserved examples of genes with known functions being found in the coastal isolates, but not the Sargasso Sea metagenomic data, or vice versa, suggesting that core regions of these diverse SAR11 genomes are relatively conserved in gene content. However, four hypervariable regions were observed, which may encode properties associated with variation in SAR11 ecotypes. The largest of these, HVR2, is a 48 kb region flanked by the sole 5S and 23S genes in the HTCC1062 genome, and mainly encodes genes that determine cell surface properties. A comparison of two closely related 'Candidatus Pelagibacter' genomes (HTCC1062 and HTCC1002) revealed a number of "gene indels" in core regions. Most of these were found to be polymorphic in the metagenomic data and showed evidence of purifying selection, suggesting that the same "polymorphic gene indels" are maintained in physically isolated SAR11 populations.

CONCLUSION

These findings suggest that natural selection has conserved many core features of SAR11 genomes across broad oceanic scales, but significant variation was found associated with four hypervariable genome regions. The data also led to the hypothesis that some gene insertions and deletions might be polymorphisms, similar to allelic polymorphisms.

Prokaryotic taxonomy and phylogeny in the genomic era: advancements and challenges ahead

Advancing prokaryotic taxonomy constitutes a contemporary academic challenge as well as practical necessity. Genome sequencing has greatly facilitated the evaluation of the current taxonomic system and the development of simpler, more portable and accurate, sequence-based alternatives to substitute for the traditional cumbersome methods. Studies based on the former genome-enabled methods reveal that existing taxonomic designations, including the species level, correspond frequently to a continuum of genetic diversity as opposed to natural groupings (e.g. biological species). Improving such artificial and often ambiguous taxonomic designations, however, will require larger genomic datasets and more carefully designed sampling of natural populations. Only then can the promise for a superior genome-based taxonomy materialize.

Primary producing prokaryotic communities of brine, interface and seawater above the halocline of deep anoxic lake L'Atalante, Eastern Mediterranean Sea

Meso- and bathypelagic ecosystems represent the most common marine ecological niche on Earth and contain complex communities of microorganisms that are for the most part ecophysiologically poorly characterized. Gradients of physico-chemical factors (for example, depth-related gradients of light, temperature, salinity, nutrients and pressure) constitute major forces shaping ecosystems at activity 'hot spots' on the ocean floor, such as hydrothermal vents, cold seepages and mud volcanoes and hypersaline lakes, though the relationships between community composition, activities and environmental parameters remain largely elusive. We report here results of a detailed study of primary producing microbial communities in the deep Eastern Mediterranean Sea. The brine column of the deep anoxic hypersaline brine lake, L'Atalante, the overlying water column and the brine-seawater interface, were characterized physico- and geochemically, and microbiologically, in terms of their microbial community compositions, functional gene distributions and [(14)C]bicarbonate assimilation activities. The depth distribution of genes encoding the crenarchaeal ammonia monooxygenase alpha subunit (amoA), and the bacterial ribulose-1,5-biphosphate carboxylase/oxygenase large subunit (RuBisCO), was found to coincide with two different types of chemoautotrophy. Meso- and bathypelagic microbial communities were enriched in ammonia-oxidizing Crenarchaeota, whereas the autotrophic community at the oxic/anoxic interface of L'Atalante lake was dominated by Epsilonproteobacteria and sulfur-oxidizing Gammaproteobacteria. These autotrophic microbes are thus the basis of the food webs populating these deep-sea ecosystems.

Divide and conquer: enriching environmental sequencing data

BACKGROUND

In environmental sequencing projects, a mix of DNA from a whole microbial community is fragmented and sequenced, with one of the possible goals being to reconstruct partial or complete genomes of members of the community. In communities with high diversity of species, a significant proportion of the sequences do not overlap any other fragment in the sample. This problem will arise not only in situations with a relatively even distribution of many species, but also when the community in a particular environment is routinely dominated by the same few species. In the former case, no genomes may be assembled at all, while in the latter case a few dominant species in an environment will always be sequenced at high coverage to the detriment of coverage of the greater number of sparse species.

METHODS AND RESULTS

Here we show that, with the same global sequencing effort, separating the species into two or more sub-communities prior to sequencing can yield a much higher proportion of sequences that can be assembled. We first use the Lander-Waterman model to show that, if the expected percentage of singleton sequences is higher than 25%, then, under the uniform distribution hypothesis, splitting the community is always a wise choice. We then construct simulated microbial communities to show that the results hold for highly non-uniform distributions. We also show that, for the distributions considered in the experiments, it is possible to estimate quite accurately the relative diversity of the two sub-communities.

CONCLUSION

Given the fact that several methods exist to split microbial communities based on physical properties such as size, density, surface biochemistry, or optical properties, we strongly suggest that groups involved in environmental sequencing, and expecting high diversity, consider splitting their communities in order to maximize the information content of their sequencing effort.

Sponge-associated microorganisms: evolution, ecology, and biotechnological potential

Marine sponges often contain diverse and abundant microbial communities, including bacteria, archaea, microalgae, and fungi. In some cases, these microbial associates comprise as much as 40% of the sponge volume and can contribute significantly to host metabolism (e.g., via photosynthesis or nitrogen fixation). We review in detail the diversity of microbes associated with sponges, including extensive 16S rRNA-based phylogenetic analyses which support the previously suggested existence of a sponge-specific microbiota. These analyses provide a suitable vantage point from which to consider the potential evolutionary and ecological ramifications of these widespread, sponge-specific microorganisms. Subsequently, we examine the ecology of sponge-microbe associations, including the establishment and maintenance of these sometimes intimate partnerships, the varied nature of the interactions (ranging from mutualism to host-pathogen relationships), and the broad-scale patterns of symbiont distribution. The ecological and evolutionary importance of sponge-microbe associations is mirrored by their enormous biotechnological potential: marine sponges are among the animal kingdom's most prolific producers of bioactive metabolites, and in at least some cases, the compounds are of microbial rather than sponge origin. We review the status of this important field, outlining the various approaches (e.g., cultivation, cell separation, and metagenomics) which have been employed to access the chemical wealth of sponge-microbe associations.

Genetic exchange across a species boundary in the archaeal genus ferroplasma

Speciation as the result of barriers to genetic exchange is the foundation for the general biological species concept. However, the relevance of genetic exchange for defining microbial species is uncertain. In fact, the extent to which microbial populations comprise discrete clusters of evolutionarily related organisms is generally unclear. Metagenomic data from an acidophilic microbial community enabled a genomewide, comprehensive investigation of variation in individuals from two coexisting natural archaeal populations. Individuals are clustered into species-like groups in which cohesion appears to be maintained by homologous recombination. We quantified the dependence of recombination frequency on sequence similarity genomewide and found a decline in recombination with increasing evolutionary distance. Both inter- and intralineage recombination frequencies have a log-linear dependence on sequence divergence. In the declining phase of interspecies genetic exchange, recombination events cluster near the origin of replication and are localized by tRNAs and short regions of unusually high sequence similarity. The breakdown of genetic exchange with increasing sequence divergence could contribute to, or explain, the establishment and preservation of the observed population clusters in a manner consistent with the biological species concept.

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.

An environmental genomics perspective on the diversity and function of marine sponge-associated microbiota

Numerous marine demosponges serve as habitats to phylogenetically complex microbial communities, with population densities exceeding those of seawater by several orders of magnitude. 16S rRNA gene-based studies have enabled a detailed phylogenetic description of the microorganisms associated with sponges, whereas environmental genomics analyses are beginning to reveal insights into their metabolic and physiological properties. Additionally, metagenomic approaches provide access to functional genes and gene clusters, thereby paving the way for the heterologous expression of novel bioactive substances from microbial symbionts of marine invertebrates.