Using pyrosequencing to shed light on deep mine microbial ecology

  • Robert A Edwards1, 2, 3, 4Email author,

    Affiliated with

    • Beltran Rodriguez-Brito1, 3,

      Affiliated with

      • Linda Wegley1,

        Affiliated with

        • Matthew Haynes1,

          Affiliated with

          • Mya Breitbart1,

            Affiliated with

            • Dean M Peterson5,

              Affiliated with

              • Martin O Saar6,

                Affiliated with

                • Scott Alexander6,

                  Affiliated with

                  • E Calvin AlexanderJr6 and

                    Affiliated with

                    • Forest Rohwer1, 2

                      Affiliated with

                      BMC Genomics20067:57

                      DOI: 10.1186/1471-2164-7-57

                      Received: 04 November 2005

                      Accepted: 20 March 2006

                      Published: 20 March 2006



                      Contrasting biological, chemical and hydrogeological analyses highlights the fundamental processes that shape different environments. Generating and interpreting the biological sequence data was a costly and time-consuming process in defining an environment. Here we have used pyrosequencing, a rapid and relatively inexpensive sequencing technology, to generate environmental genome sequences from two sites in the Soudan Mine, Minnesota, USA. These sites were adjacent to each other, but differed significantly in chemistry and hydrogeology.


                      Comparisons of the microbes and the subsystems identified in the two samples highlighted important differences in metabolic potential in each environment. The microbes were performing distinct biochemistry on the available substrates, and subsystems such as carbon utilization, iron acquisition mechanisms, nitrogen assimilation, and respiratory pathways separated the two communities. Although the correlation between much of the microbial metabolism occurring and the geochemical conditions from which the samples were isolated could be explained, the reason for the presence of many pathways in these environments remains to be determined. Despite being physically close, these two communities were markedly different from each other. In addition, the communities were also completely different from other microbial communities sequenced to date.


                      We anticipate that pyrosequencing will be widely used to sequence environmental samples because of the speed, cost, and technical advantages. Furthermore, subsystem comparisons rapidly identify the important metabolisms employed by the microbes in different environments.


                      Banded iron formations started appearing ~3,700 million years ago when localized sea floor cyanobacterial photosynthesis raised oxygen concentrations high enough that dissolved iron precipitated. That iron powered the industrial revolution. The Soudan Iron Mine in Minnesota, USA was active from 1884 to 1962, and during this period 17.9 million tons of iron ore, primarily hematite, were removed. Nowadays the mine is used as a state park and as a facility for high-energy physics experiments.

                      Metagenomics is a term used to describe "the functional and sequence-based analysis of the collective microbial genomes contained in an environmental sample"[1, 2]. Random shotgun sequencing of DNA from natural communities has been used to characterize seawater, sediment, and fecal viral communities [25], as well as the microbial communities in soil, whale falls, seawater and the Iron Mountain acid mine drainage (AMD) [610]. Comparative metagenomics was introduced recently[6], identifying those sets of genes that distinguish environmental samples. For example, samples from the surface of the ocean contain significantly more photosynthetic genes than soil or other samples[6, 8, 10]. We have used comparative metagenomics to characterize the metabolic potential of different environments, and identify those genes, pathways, and subsystems that are more common in any particular environment [11].

                      Most current sequencing is a modification of the classical Sanger method, where extending DNA fragments are stopped by the random incorporation of a fluorescently labeled ddNTP. The different-sized fragments are then separated using capillary gel electrophoresis and detected with a LASER. Pyrosequencing is a fundamentally different methodology because only one dNTP is added into the reaction at a time [1214]. If there is a complementary base, then the DNA polymerase catalyzes the reaction and releases pyrophosphate. ATP sulfurylase uses the pyrophosphate to produce ATP in the presence of adenosine 5' phosphosulfate (APS). A Charge-Coupled Device (CCD) measures the light produced when the ATP is used by luciferase to convert luciferin to oxyluciferin. 454 Life Sciences has scaled this process up to be massively parallel, determining the composition of more than 300,000 sequences at once, for approximately the same price as 96 to 192 sequencing reactions performed using traditional chemistries[12]. In addition to the massive parallelization, the 454 technology does not require cloning of the environmental samples, thus eliminating many of the problems that are associated with this step of metagenomics[2].

                      This report describes the first application of pyrosequencing to environmental samples. From this sequence data, we identify the 16S rDNA sequences present in the sample, and apply new annotation methods to this data using the SEED database[15]. This paper also describes a comprehensive statistical treatment of the genes identified in each sample using a completely novel methodology that exploits the differences between metagenome sequences. We demonstrate that completely unique microbial communities inhabit proximate environments joined by a common watercourse, and that using metagenomics we can identify the unique metabolic potentials prevalent in each environment such as their mechanisms of iron acquisition and respiration. The integration of pyrosequencing, subsystems analysis, comparative metagenomics, statistics, hydrogeology, and chemistry provides a comprehensive systems analysis of the Soudan Mine.

                      Results and discussion

                      Description of the environmental samples

                      The two environments sampled within the Soudan Mine are shown in Figure 1. Groundwater is not very abundant in the banded iron formations of the mine at our sampling depth of 714 m below the surface. However, small amounts of water emerge steadily from exploration boreholes that extend downward from the deepest level of the mine. The two sets of samples were collected from water trickling out of two such boreholes that are separated by about 100 meters (Figure 1). This water is a calcium, sodium, chloride solution about twice as salty as seawater. It is anoxic, with up to 150 ppm of dissolved ferrous iron and variable enrichments of several trace elements. At both locations, the water emerging from the boreholes produces cm-scale "Black" environments that appear to extend down into the borehole. The water flowing away from each borehole, on the floor of the mine tunnel, is exposed to the oxygenated mine atmosphere, and transitions to a sequence of "Red" environments within a few cm of the orifices. The oxidized environments are continuously fed by anoxic water flowing from the boreholes. The water in the borehole, which yielded the Black sample, as well as a number of similar sites found throughout the mine, has a pH of 6.70 and redox potential of -142 mV. Some of the Black areas are associated with bubbling of gas. The Black sediment contained 5.8 × 105 microbes per ml. X-ray diffraction analyses of the minerals in this area show that chlorite-serpentine [(Mg,Al)6(Si,Al)4O10(OH)8], clinochlore, ferroan [(Mg,Fe)6(Si,Al)4O10(OH)8], quartz, and silinaite [LiNaSiO5·HCl] are present in the Black sediments. Water slowly flows from the borehole into the stream running down the main mine tunnel. As the water comes in contact with oxygen in the passageway, the pH rapidly decreases to 4.37 and redox potential increases to -8 mV. The Red sample contained 1.2 × 106 microbes per ml, and these sediments include goethite [FeO(OH)], followed by szaibelyite [(Mg,Mn)BO2(OH)], and sussexite [(Mn,Mg)BO2(OH)].
                      Figure 1

                      Sampling from the Soudan Mine. The Soudan Mine is an Algoma-type Iron Formation rich in hematite. Panel A shows a cross-section of the mine looking East-North-East at 78.5°. Panel B depicts a three dimensional view of the mine, including the cross-section shown in Panel A, and with the sampling sites shown for the "Red" and "Black" samples. Panel C shows the overall flow of water in the mine at level 27, located 714 meters below the surface (Panel D). Panels E and F show a close up of the two sampling sites.

                      The first two pyrosequences of environmental samples

                      DNA was purified from the two samples, amplified using the GenomiPhi procedure (GE Healthcare, Piscataway, NJ), and then sequenced by 454 Life Sciences. A summary of the sequence characteristics determined using the pyrosequencing technique is shown in Table 1. The raw sequence reads and quality scores [see Additional files 6 and 7] are provided in compressed format.
                      Table 1

                      Summary of pyrosequence data from the Soudan Mine


                      Red Sample

                      Black Sample

                      Number of Sequences



                      Total Length of Sequences

                      35,439,683 bp

                      38,502,057 bp

                      Average Length of Sequences

                      106.0 bp

                      99.1 bp

                      Average Quality Score1






                      1 The quality score of each base was provided by 454 Life Sciences, and is analogous to the Phred score of Sanger Sequencing methods [26]. The value cited here is the average of the mean quality score per sequence.

                      2The skew was calculated by comparing the dinucleotide frequencies within each library. A similar analysis performed on Bacterial and Eukaryotic genomes sampled at random yielded entropies of 2.63 and 2.66 respectively.

                      The two samples produced more than 70 Mbp of sequence data from over 700,000 sequences, and there was no significant skew in the sequence data (as measured by dinucleotide frequency) when the data generated by pyrosequencing was compared to complete genome sequences.

                      16S rDNA analysis of the samples

                      The two sequence libraries were compared to the 16S rDNA database from the Ribosomal Database Project[16]. As shown in Figure 2, the Black sample was dominated by Actinomycetales such as Brevibacterium and Corynebacterium that volatilize sulfur via an organic intermediate and can also break down complex heterocyclic and polycyclic ring structures[17, 18]. In contrast, members of the Chromatiales, including the genera Chromatiaceae, Thiobacillus, and Halothiobacillus, dominate the Red sample. These chemoautotrophic Bacteria often use the Calvin-Benson-Bassham cycle to fix CO2 through the oxidation of iron or sulfur, and consequently they would be expected to be present in samples from an iron-rich deposit. These two communities are fundamentally different both from each other and from the community identified in the Iron Mountain metagenome[7]. The community in the Red sample has a much higher species richness than the Black sample, and the differences between the Soudan and Iron mines reflect the iron composition (hematite versus pyrite), temperature, and pH of the various environments[7].
                      Figure 2

                      Composition of the 16S rDNA sequences from the two samples and comparison of 16S sequences from the 454 libraries and a traditional clone library. The percentage of all sequences from each library in each of the orders is shown for the 454-sequenced Black sample (solid black bars; n = 24), the 454 sequenced red sample (solid red bars; n = 76), and the PCR amplified clone library (hatched red bars; n = 91).

                      A16S clone library was created from the Red sample to validate the 454 sequencing approach. Ninety-six clones were sequenced using traditional techniques, and compared to the 16S rDNA database from the Ribosomal Database Project [16]. The congruity between the 16S genes sequenced in the 454 library and the 16S sequences from the clone library, as shown in Fig. 2, is quite remarkable.

                      We also used the 16S sequences to evaluate the randomness of the library. An analysis of 160 bacterial genome sequences in the SEED database [15, 19] with annotated 16S genes showed that about 1 in 105 bases is from a 16S gene. Based on this estimate, as a rule of thumb the Soudan samples are expected to contain approximately 3,000 bases of 16S sequence in total, or approximately 30 sequences. Twenty four sequences were found to have significant similarity (with an E value less than 1 × 10-5 and a match of 50 bp or more) to 16S rDNA from the Black sample and seventy six sequences were found to have significant similarity to 16S rDNA from the Red sample.

                      Metabolic potential from the metagenome library

                      Sequences from both libraries were compared to the SEED database, a curated database of microbial genomes [15, 19]. The annotations using the SEED interface primarily occur through the development of subsystems, a technique pioneered by the Fellowship for Interpretation of Genomes[15, 20]. Subsystems are groups of genes that function together, such as the genes whose products are involved in a metabolic pathway, or the group of genes whose products make a cellular structure. A summary of the subsystem hits are shown in Figure 3, and all matches to subsystems are provided as supplemental data [see Additional files 1, 2 and 3]. These subsystems show that the pyrosequencing generates sequences that represent a large swathe of central metabolism in each of the environments. Common metabolic potential that is expected to be present in sulfur-utilizing chemoautotrophs is represented in the mine libraries, including the Calvin-Benson cycle, inorganic sulfur assimilation, amino acid biosynthetic genes, and so on. The comparison of the subsystem similarities suggested the simple hypothesis that groups of genes (or subsystems) important to a particular environment will be enriched in that environment. To distinguish between ecologically important differences and differences caused by sampling error, a method was devised to identify those subsystems that are statistically significantly overrepresented in one sample when compared to another [11].
                      Figure 3

                      Subsystems in the Red and Black Samples. The occurrence of classes of subsystems is shown as a percent of all subsystems in each sample for the Red and Black samples. Notes and abbreviations: The subsystem class "Glu, Asp" also contains Gln and Asn. The subsystem class "Lys, Thr" also contains Met and Cys. CHO: Carbohydrates; sacch: saccharides; Extracell. Poly: Extracellular polysaccharides; Myco: Mycobacterial cell wall; Gm: Gram stain positive (+) or negative (-); Clust: clusters; RFN: Riboflavin; T: Transporters; Mot: Motility; N: Nitrogen; Resp: Respiration; e-: electron; S: Sulfur.

                      Subsystems enriched in the Black or Red samples

                      Table 2 shows subsystems that were determined to be statistically more common, with 95% confidence, in either the Red or Black samples from the Soudan mine. The subsystems that are overrepresented in a metagenome can yield significant insights into the microbial ecology of the environment. A few specific examples are detailed below.
                      Table 2

                      Subsystems statistically more likely to be present in either the Red or Black samples. These subsystems are more frequently found among sequences from either the Red or Black samples with a sample size of 5,000 proteins, 20,000 repeated samples, and P < 0.05.

                      Red Sample (Oxidized, pH4.37, E h -8)

                      Black Sample (Reduced, pH 6.70, E h -142)

                      Amino Acids and Derivatives

                      Arginine biosynthesis

                      Urea decomposition

                      Tryptophan synthesis

                      Chorismate synthesis

                      Asp-Glu-tRNA(Asn-Gln) transamidation

                      Branched-chain amino acid biosynthesis

                      Histidine biosynthesis

                      Isoleucine degradation


                      Leucine biosynthesis


                      Leucine degradation and HMG-CoA metabolism


                      Valine degradation


                      Methionine salvage



                      Glyoxylate synthesis

                      Cell Division and Cell Cycle



                      Cell Wall and Capsule


                      N-linked glycosylation in Bacteria


                      Teichoic acid biosynthesis

                      Cofactors, Vitamins, Prosthetic Groups, Pigments

                      Folate biosynthesis

                      Coenzyme A biosynthesis in pathogens

                      Methylglyoxal metabolism

                      Molybdopterin biosynthesis

                      Pyruvate metabolism I: anaplerotic rx, PEP


                      Ubiquinone biosynthesis

                      Polyisoprenoid biosynthesis

                      Ubiquinone menaquinone-cytochrome c reductase

                      NAD and NADP cofactor biosynthesis global

                      Riboflavin metabolism

                      Coenzyme PQQ synthesis


                      Pyrroloquinoline quinone biosynthesis


                      Siderophore enterobactin biosynthesis


                      Siderophore enterobactin biosynthesis and ferric enterobactin transport


                      Thiamin biosynthesis

                      DNA metabolism

                      DNA repair, bacterial


                      Fatty Acids and Lipids

                      Fatty acid metabolism

                      Glycerolipid and glycerphospholipid metabolism

                      Fatty acid oxidation pathway


                      Membrane Transport

                      ABC transporter maltose

                      ABC transporter ferrichrome


                      ABC transporter heme


                      CbiQO-type ABC transporter systems


                      Sodium hydrogen antiporter

                      Metabolism of aromatic compounds

                      Phenylacetate pathway of aromatic compound degradation

                      Homogentisate pathway of aromatic compound degradation

                      Motility and Chemotaxis

                      Bacterial chemotaxis




                      Nitrogen Metabolism



                      Nucleosides and Nucleotides

                      De novo purine biosynthesis


                      Ribonucleotide reduction


                      Protein Metabolism

                      Ribosome LSU bacterial

                      Phenylpropionate degradation

                      Ribosome SSU bacterial


                      Translation factors bacterial


                      Universal GTPases


                      Protein degradation



                      F0F1-type ATP synthase

                      NiFe hydrogenase maturation

                      Terminal cytochrome C oxidases




                      Membrane-bound Ni, Fe-hydrogenase


                      Na(+)-translocating NADH-quinone oxidoreductase and rnf-like group of electron transport complexes


                      Respiratory complex I


                      Respiratory dehydrogenases 1


                      RNA metabolism

                      Polyadenylation bacterial


                      RNA polymerase bacterial


                      tRNA aminoacylation


                      Stress response

                      Glutathione redox metabolism


                      ppGpp biosynthesis


                      Sulfur Metabolism

                      Sulfate assimilation



                      Resistance to fluoroquinolones

                      Several subsystems involved in iron uptake and utilization such as siderophores and ABC transporters for ferrichrome and heme are more common in the Black sample. The overall concentration of iron at the two sites was similar (Table 3; Figure 5). However, the iron in the Black sample is present as either Fe2+ dissolved in the water or as ferroan [(Mg,Fe)6(Si,Al)4O10(OH)8]. In either case, the ferrous iron can not be assimilated biologically, and the microbes are forced to scavenge for the limited ferric iron (Fe3+) available. In contrast, in the Red sample, goethite [FeO(OH)] is present and ferric iron is more readily available for biological utilization. The Black sample is enriched for amino acid degradation pathways and microbes may be assimilating nitrogen or carbon through these pathways. It is not currently apparent from where free amino acids would be supplied.
                      Table 3

                      Water chemistry from Soudan Mine. No significant differences were found for Ca, Mg, Na, K, Li, Al, Mn, Sr, Ba, Si, Cr, Co, Ni, Cu, Zn, As, Se, Rb, Cd, Cs, Pb, total alkalitity, lactate, acetate, formate, chlorate, oxalate, and trace elements.




                      Temp (°C)






                      redox (mV)



                      Fe (ppm)



                      Total N (ppm)



                         • NH4



                         • NO3



                         • NO2



                      SO4 (ppm)



                      PO4 (ppm)



                      B (ppm)



                      Mo (ppm)



                      W (ppm)



                      Tl (ppm)



                      U (ppm)



                      Figure 5

                      Cations and Anions found in the Soudan Mine. The pie chart shows the abundance of cations and anions found in the mine. The numbers in parentheses are the concentrations (in ppm) of each ion in the "Black" and "Red" samples respectively. The minor ions are shown expanded in the rightmost pie.

                      The respiratory complexes and cytochrome-C oxidases are more commonly found in the sample from the oxidized environment (the Red sample; Table 2). Respiration proceeds via multiple electron transfer steps (Figure 6). In an aerobic environment, electrons are passed from hydrogenases to quinones (e.g., ubiquinone, quinone, menaquinone, and plastoquinone) and then to cytochromes resulting in the conversion of oxygen to water. In anaerobic environments the electrons are shuffled through nitrite and nitrate reductases, reducing NO3 first to NO2 and then to N2 gas. The Black sample is enriched for these denitrification genes suggesting that the latter pathway predominates while the Red sample is enriched for components of the aerobic respiratory pathway. Moreover, the Black sample had a lower concentration of free nitrate than the Red sample, presumably because nitrate is being used as an electron acceptor during respiration (although nitrite was below the level of detection in both samples; Table 3).
                      Figure 6

                      Respiration in aerobic and anaerobic environments. Among other potential pathways in the Soudan mine, electrons are transferred from hydrogenases to either cytochromes and then to oxygen to produce water in an oxidative environment, or via nitrate and nitrite reductases (denitrification) in anaerobic environments. Genes encoding the hydrogenases, respiratory complexes, and terminal cytochromes of the aerobic sample were significantly more abundant in the Red (oxidized) sample, while genes encoding the hydrogenases and denitrification genes were more abundant in the Black (reduced) sample. After Vassieva, O. [25]

                      This analysis demonstrates that by combining pyrosequencing, subsystems analysis, and comparative metagenomics the microbiology of different environments can be correlated with the chemistry and hydrogeology of those environments to identify significant ecological differences between them.

                      Comparisons between Soudan and Iron Mountain communities

                      A previous study used Sanger sequencing to determine the metagenome of the Iron Mountain community[7]. The environmental differences (such as the difference in temperature) account for the predominant differences between the microbial communities. The organismal differences are reflected in the individual biochemistries of the samples [see Additional files 4 and 5]. For example, the AMD metagenome contains significantly more occurrences of Archaea-specific subsystems such as those involved in protein biosynthesis than the Soudan samples. The AMD sample has a preference for CO2 fixation and simple carbohydrate metabolism when compared to either of the Soudan samples. There are also many currently unexplained differences between subsystems found in these environments that must relate the biology of the organisms to the chemistry of the environment.

                      Comparisons between Soudan and other metagenome sequences

                      The SEED database used for these studies contained 351 subsystems. The vast majority (83%) of subsystems were present in one or more of the sequenced metagenomes, and over half (52%) of the subsystems are present in every metagenome. A comparison of the subsystem classification reveals trends between the metagenomes (Figure 4). For example, oxygenic photosynthesis is prevalent in samples that are naturally illuminated such as the Sargasso Sea[10]. This analysis also suggested that phosphorous metabolism is more prevalent in oceanic surfaces rather than terrestrial environments. Comparisons of the Minnesota Farm metagenome[6] with the Soudan Mine metagenomes, also from Minnesota, showed important differences in the production and consumption of secondary metabolites, membrane transport, and fatty acid metabolism. The complete lists of statistically significantly different subsystems between both Red or Black samples and each of the previously published metagenomes are supplied as supplemental material [see Additional files 4 and 5].
                      Figure 4

                      Subsystems present in different metagenome sequences. The subsystems present in the Soudan samples, the Iron Mountain AMD sample, the Minnesota Farm and the Sargasso Sea are shown grouped by family. The red x corresponds to very low abundance or complete absence of that family of subsystems. The size of the circle represents the proportion of sequences seen within that family of subsystems.


                      This is the first metagenome analysis performed using pyrosequencing, which is approximately 10 to 30 times cheaper than current Sanger sequencing. Pyrosequencing also eliminates the need for cloning, thus removing the potential for both aberrant recombinants in the surrogate host and for cloning-related artifacts such as counterselection against potentially toxic genes such as those found on phages[2]. The main concerns with current pyrosequencing technology are the short length of sequence fragments (average of 105 bp in this study), and the requirement to use whole genome amplification to generate sufficient DNA for sequencing from environmental libraries The former may make it difficult to accurately assemble genomes in the absence of a scaffold, while the later may bias these analyses. Our preliminary unpublished data suggests that the whole genome amplification bias is minimal, and is preferentially towards the ends of linear pieces of DNA (Haynes, Rayhawk, Edwards, Rohwer; unpublished). Since these biases are applied equally to both libraries, they will be negated during the comparative study to highlight differences between metagenomes. Nonetheless, the short fragments are sufficient to determine statistically significant differences between metagenomes that reflect the most likely biology occurring in each environment. The low cost, high yield of pyrosequencing combined with statistical analyses on the abundance of subsystems in the samples allows the rapid identification of key processes driving the metabolism of different environments.

                      The systems approach of integrating biology, chemistry, and geology has yielded significant insights into the metabolism of two environments in the Soudan Mine. The oxidized sample is using aerobic respiratory pathways while the reduced sample is using anaerobic pathways. Nitrogen assimilation, iron acquisition, and sulfur metabolism are all differentiated between these two samples from close proximity within the same mine. However, many more significant differences between the samples remain unexplained by our current knowledge of bacterial physiology and metabolism. Explaining these differences will be a grand challenge for the future. By combining pyrosequencing, subsystems analysis, comparative metagenomics, and statistics, Occam has used his razor on metagenomics.


                      Sample collections, microbial enumeration, and DNA extraction

                      Samples were collected from several sites in the Soudan Mine. This analysis concerns the sample collection at two sites on Level 27 (714 m below the surface; Figure 1). Water and sediments were sampled from the two locations shown in Figure 1 giving the "Black" (reduced) sample and "Red" (oxidized) sample. Microbes were concentrated from these samples by filtration with 0.22 μm Sterivex units. Microbial counts were enumerated by staining the samples with SYBR-Gold (Invitrogen, Carlsbad, CA) and visualization with an epifluorescent microscope [21]. DNA was extracted from the microbial sample using either the Ultra Clean Soil DNA Kit or Power Soil Kit (MolBio, Boulder, CO). The DNA was amplified with GenomiPhi (GE Healthcare, Piscataway, NJ) in an Eppendorf thermal cycler (Eppendorf, Westbury, NY) using multiple reactions containing 50–100 ng of the isolated DNA as template and the manufacturer's recommended protocols. After amplification, the resulting DNA was purified with silica columns (Qiagen, Valencia, CA) and concentrated by ethanol precipitation. The DNA was resuspended in water to a final concentration of 0.3 mg/ml. Approximately 10 μg of each sample was sequenced using the pyrosequencing technology (454 Life Sciences, Branford, CT).

                      Bacterial-specific primers 27F (5'-AGAGTTTGATCMTGGCTCAG-3') and the universal 1492R primer (5'-TACGGYTACCTTGTTACGACTT-3') [22] were used to amplify the 16S rDNA genes. PCR products were cloned into the pCR®4-TOPO® vector as recommended by the manufacturer (Invitrogen, Carlsbad, CA).

                      Water and mineral analyses

                      Water samples were collected by filtering the water through 0.2 μm filters into clean bottles. Field measurements of pH, Eh, temperature and conductivity were conducted in situ. The sediment samples were collected as slurries with a pipette resting on and in the sediments. Those slurries were transferred to clean centrifuge tubes, allowed to settle by gravity and then the fluid was decanted.

                      Major anions in the water were determined by GC (Dionex IGS-2000, Sunnyvale, CA) and major and trace elements by ICP/MS (Thermo Electron PQ ExCell, Franklin, MA). The mineral identifications are based on XRD (Bruker-AX D500 X-ray Diffractometer, Germany) measurements. The X-ray peaks were relatively small. Much of the sediment was apparently not well crystallized.

                      Sequence analysis

                      The unassembled sequences provided by 454 were compared to the SEED database using the BLASTX algorithm on the Teragrid cluster at Argonne National Laboratories[15, 23]. All BLAST searches were performed using an expect value cutoff of 1 × 10-5. At this cutoff approximately 3 of the observed hits would be expected to occur at random[23].

                      The BLASTN algorithm was used to identify 16S genes from release 9 of the RDP database [16, 24]. These BLAST searches were also performed using an expect value cutoff of 1 × 10-5 and a minimum sequence match length of 50 nt.

                      Statistical analyses of metagenome datasets

                      The statistical analysis of subsystems present in each sample was performed essentially as described elsewhere [11]. The presence or absence of subsystems between two data sets was determined using 20,000 replicates of samples of 5,000 subsystems each. The 95% confidence interval for the median was constructed using the 0.025 and 0.975 percentiles.



                      The authors are grateful to Bill Miller, Director of the Soudan Facility, for arranging the sampling trips, ploughing through paperwork, and bringing these fascinating microbial communities to our attention. Thanks to Jim Essig, for guiding us to Level 10, arranging our sampling protocols, and generally helping out. Tony Zavodnick, Paul Paulisich, and Jack Zorman provided invaluable guidance and explanations of mining at the Soudan Mine site and the whole Soudan mine and facility crews for making our sampling trip so enjoyable. In addition, we thank Robert Olson (Argonne National Labs) for assistance with the computational analysis. This work was supported by a grant NSF DEB-BE 04-21955 from the NSF Biocomplexity program (to FR).

                      Authors’ Affiliations

                      Department of Biology, San Diego State University
                      Center for Microbial Sciences, San Diego State University
                      Computational Science Research Center, San Diego State University
                      Fellowship for Interpretation of Genomes
                      Natural Resources Research Institute,Department of Geological Sciences, University of Minnesota
                      Department of Geology and Geophysics, University of Minnesota


                      1. Riesenfeld CS, Schloss PD, Handelsman J: Metagenomics: genomic analysis of microbial communities. Annu Rev Genet 2004, 38:525–552.View ArticlePubMed
                      2. Edwards RA, Rohwer F: Viral metagenomics. Nat Rev Microbiol 2005,3(6):504–510.View ArticlePubMed
                      3. Breitbart M, Felts B, Kelley S, Mahaffy JM, Nulton J, Salamon P, Rohwer F: Diversity and population structure of a near-shore marine-sediment viral community. Proc R Soc Lond B Biol Sci 2004,271(1539):565–574.View Article
                      4. Breitbart M, Hewson I, Felts B, Mahaffy JM, Nulton J, Salamon P, Rohwer F: Metagenomic analyses of an uncultured viral community from human feces. J Bacteriol 2003,185(20):6220–6223.View ArticlePubMed
                      5. Breitbart M, Salamon P, Andresen B, Mahaffy JM, Segall AM, Mead D, Azam F, Rohwer F: Genomic analysis of uncultured marine viral communities. Proc Natl Acad Sci U S A 2002,99(22):14250–14255.View ArticlePubMed
                      6. Tringe SG, von Mering C, Kobayashi A, Salamov AA, Chen K, Chang HW, Podar M, Short JM, Mathur EJ, Detter JC, Bork P, Hugenholtz P, Rubin EM: Comparative metagenomics of microbial communities. Science 2005,308(5721):554–557.View ArticlePubMed
                      7. Tyson GW, Chapman J, Hugenholtz P, Allen EE, Ram RJ, Richardson PM, Solovyev VV, Rubin EM, Rokhsar DS, Banfield JF: Community structure and metabolism through reconstruction of microbial genomes from the environment. Nature 2004,428(6978):37–43.View ArticlePubMed
                      8. DeLong EF: Microbial community genomics in the ocean. Nat Rev Microbiol 2005,3(6):459–469.View ArticlePubMed
                      9. Cann AJ, Fandrich SE, Heaphy S: Analysis of the virus population present in equine faeces indicates the presence of hundreds of uncharacterized virus genomes. Virus Genes 2005,30(2):151–156.View ArticlePubMed
                      10. Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, Wu D, Paulsen I, Nelson KE, Nelson W, Fouts DE, Levy S, Knap AH, Lomas MW, Nealson K, White O, Peterson J, Hoffman J, Parsons R, Baden-Tillson H, Pfannkoch C, Rogers YH, Smith HO: Environmental genome shotgun sequencing of the Sargasso Sea. Science 2004,304(5667):66–74.View ArticlePubMed
                      11. Rodriguez-Brito B, Rohwer F, Edwards RA: An application of statistics to comparative metagenomics. BMC Bioinformatics 2006, in press.
                      12. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, McDade KE, McKenna MP, Myers EW, Nickerson E, Nobile JR, Plant R, Puc BP, Ronan MT, Roth GT, Sarkis GJ, Simons JF, Simpson JW, Srinivasan M, Tartaro KR, Tomasz A, Vogt KA, Volkmer GA, Wang SH, Wang Y, Weiner MP, Yu P, Begley RF, Rothberg JM: Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005.
                      13. Ronaghi M, Karamohamed S, Pettersson B, Uhlen M, Nyren P: Real-time DNA sequencing using detection of pyrophosphate release. Anal Biochem 1996,242(1):84–89.View ArticlePubMed
                      14. Ronaghi M, Uhlen M, Nyren P: A sequencing method based on real-time pyrophosphate. Science 1998,281(5375):363, 365.View ArticlePubMed
                      15. Overbeek R, Begley T, Butler R, Choudhuri J, Chuang H, Cohoon M, de Crécy-Lagard V, Diaz N, Disz T, Edwards R, Fonstein M, Frank E, Gerdes S, Glass E, Goesmann A, Hanson A, Iwata-Reuyl D, Jensen R, Jamshidi N, Krause L, Kubal M, Larsen N, Linke B, McHardy A, Meyer F, Neuweger H, Olsen G, Olson R, Osterman A, Portnoy V, Pusch G, Rodionov D, Rückert C, Steiner J, Stevens R, Thiele I, Vassieva O, Ye Y, Zagnitko O, Vonstein V: The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes . Nucleic Acids Res 2005.
                      16. Cole JR, Chai B, Farris RJ, Wang Q, Kulam SA, McGarrell DM, Garrity GM, Tiedje JM: The ribosomal database project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res 2005,33(Database issue):D294–6.View ArticlePubMed
                      17. Klubek B, Schmidt C, Burnham H: Characterization of soil bacteria that desulphurize organic sulphur compounds .1. Classification and growth studies. Microbios 1996,88(357):223–236.
                      18. Turtura GC, Perfetto A, Lorenzelli P: Microbiological investigation on black crusts from open-air stone monuments of Bologna (Italy). Microbiologica 2000,23(2):207–228.PubMed
                      19. The SEED[http://​theseed.​uchicago.​edu/​FIG/​index.​cgi]
                      20. Overbeek R, Disz T, Stevens R: The SEED: A peer-to-peer environment for genome annotation. Commun ACM 2004,47(11):46–51.View Article
                      21. Noble RT, Fuhrman JA: Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria. Aquat Microb Ecol 1998,14(2):113–118.View Article
                      22. Amann RI, Ludwig W, Schleifer KH: Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev 1995,59(1):143–169.PubMed
                      23. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403–410.PubMed
                      24. Ribosomal database project - II[http://​rdp.​cme.​msu.​edu/​]
                      25. FIG Subsystem Forum[http://​www.​subsys.​info]
                      26. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res 1998,8(3):186–194.PubMed


                      © Edwards et al. 2006

                      This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.