A sample was collected from the whole white mat, which was directly attached to the fountain substrate, in the upper basin of the Alpena library fountain site (45 3.747′N, 83 25.872′W) on July 24th, 2012, and preserved in RNAlater (Ambion) as instructed by the manufacturer. Sulfide concentrations were determined by the methylene blue method (Cline, 1969) using Procedure 1 described by Reese et al. (2011). Dissolved oxygen was measured by a YSI Pro Dissolved Oxygen probe (Digital Professional Series). At the Isolated Sinkhole site (45 10.711′N, 83 09.173′W), a sample of white mat was skimmed off of sediment from a box core collected on September 6th 2008 on board the R/V Laurentian using a 60 ml syringe with attached flexible Tygon tubing, and immediately frozen at -20°C. After thawing, the white mat was manually separated from surrounding sediment and washed gently with a filter-sterilized NaCl solution to remove visible sediment particles. We note that even such gentle washing could reduce the diversity detected, and since it was only applied to the Isolated Sinkhole sample (where attached sediments were present), comparisons between the two samples could theoretically have been affected. However, it is unlikely to influence the results of our main focus, the large sulfur-oxidizing bacteria. DNA from both samples was extracted using the FastDNA Spin Kit for Soil (MP Biomedicals). Both samples were then sequenced in September of 2012 by standard methods using the Illumina HiSeq with paired-end 100 base-length reads at the University of Michigan DNA Sequencing Core. Sequence reads were deposited to the NCBI Sequence Read Archive (see Supplementary Table S1 for accession numbers).

16S rRNA gene sequences were assembled from the metagenomic datasets using EMIRGE (Miller et al., 2011). EMIRGE assemblies were performed with 120 iterations against the Silva small subunit reference database v. 111 (Pruesse et al., 2007) using quality filtered and trimmed reads as described below. Nine of the resulting 22 EMIRGE sequences had top BLAST hits within the Beggiatoaceae. These sequences were screened for chimeric sequences using Pintail (Ashelford et al., 2005) in mothur v. 1.33 (Schloss et al., 2009), and no chimeras were detected. Initial difficulties in resolving the phylogeny of one the EMIRGE sequences suggested the possible presence of an intron. The EMIRGE sequences were aligned with 16S rRNA gene sequences from taxa within the Beggiatoaceae that contain introns, as well as taxa whose 16S rRNA gene sequences do not contain introns, using Geneious (9.1.3) (Kearse et al., 2012). This alignment identified homologs to introns observed by Salman et al. (2012). The introns were manually removed and the 16S rRNA genes were aligned using the SINA aligner (Pruesse et al., 2012). The alignment was manually corrected in ARB (Ludwig et al., 2004) before phylogenetic analyses. Alignments were end-clipped so that all sequences were of equal length and positions with >50% gaps were removed (final alignment 28- 1271 E. coli base positions). Tree building was performed within Geneious. A neighbor joining phylogenetic tree with Jukes-Cantor distance correction (Jukes and Cantor, 1969; Saitou and Nei, 1987) and a heuristic search maximum parsimony tree (Fitch, 1971) were constructed with PAUP 4.0a148 (Swofford, 2003). Both trees were bootstrapped (Felsenstein, 1985) with 1000 replicates and consensus trees were generated. A consensus Bayesian interference phylogenic tree was constructed with MrBayes 3.2.6 with 1,000,000 generational runs (Huelsenbeck and Ronquist, 2001) using six substitution categories (aka. the general time reversible (GTR) model), gamma distribution and an invariant site nucleotide substitution rate (G++I). The trees were saved every 100 generations and posterior probabilities calculated after discarding the first 20% of trees. For maximum likelihood analyses, the GTR nucleotide substitution model with G+I rate variation was selected with JModelTest v.2.1.10 (Darriba et al., 2012). The Maximum-Likelihood (ML) analysis was computed in RAxML v.7.2.8 (Stamatakis, 2006) with 1000 bootstrap replicates. The ML tree was selected for the phylogenic distance tree. Branches were collapsed when ML bootstrap values fell below 50%.

rpoB phylogenies were calculated using full-length or nearly full-length rpoB gene sequences from the metagenomes, and from freshwater and marine Beggiatoaceae available from isolate genomes and metagenomes at the Joint Genome Institute Integrated Microbial Genomes and Metagenomics (IMG/M) website (Markowitz et al., 2008). Nucleotide sequences were translated in ARB (Ludwig et al., 2004). Protein sequences were aligned using M-Coffee (Moretti et al., 2007), re-imported into ARB, and the corresponding nucleotide sequences re-aligned according to the amino acid alignments. Additional sequences were aligned with the ARB integrated aligner. Maximum likelihood analyses were computed using the rapid bootstrap algorithm (1000 replicates) in RAxML v.8.0.24 (Stamatakis, 2006). All positions with missing information and gaps were masked from the alignments prior to analysis, for final lengths of 989 and 2967 columns for the protein and nucleotide alignments, respectively. Evolutionary models were selected using the corrected Akaike information criterion (AICc) with ProtTest v.2.4 (Abascal et al., 2005) and JModelTest v2.1.5 (Darriba et al., 2012). Maximum likelihood analyses of protein sequences were calculated with the LG amino acid replacement matrix (Le and Gascuel, 2008) and gamma-distributed rates, with proportions of invariant sites and the alpha parameter estimated from the data. Maximum likelihood analyses of nucleotide sequences were calculated with the general time reversible model, gamma-distributed rates, and proportions of invariant sites and the alpha parameter estimated from the data.

Sequence reads were dereplicated (100% identity over 100% lengths), then trimmed with the adaptive read trimmer, Sickle (Joshi and Fass, 2011). De novo assembly of both metagenomic data sets was done using IDBA UD with a minimum kmer of 52, maximum kmer of 92, and a step size of 8 (Peng et al., 2012). Emergent self-organizing maps of tetranucleotide frequencies were then used to sort scaffolds of length greater than 4,000 bp into genomic bins as described previously (Dick et al., 2009). The Joint Genome Institute Integrated Microbial Genomes and Metagenomics (IMG/M) system was used to annotate genes on assembled scaffolds in the bins of interest (Markowitz et al., 2008). Annotations of genes of interest in these bins were confirmed manually using BLASTp to the National Center for Biotechnology Information non-redundant database with cutoffs of 150 for bitscore and 30% for percent identity (Pruitt et al., 2005), and results were evaluated manually. Scaffold sequences and gene annotations were deposited to NCBI (BioProject PRJNA340051 and PRJNA340052) and Integrated Microbial Genomes (Genome ID 3300003868 and 3300003874), respectively (see Supplementary Table S1 for accession numbers of individual genome bins).

Bin completeness was estimated by searching each bin for 36 conserved phylogenetic markers (Ciccarelli et al., 2006). The percentage completeness reported is the number of these conserved phylogenetic markers present out of the expected 36. The taxonomic identification of each bin was determined through a combination of methods including 16S identity, IMG/M, and PhyloSift (Darling et al., 2014). We first conducted a BLASTn search of scaffolds against the Silva Small Subunit RNA database release 119 using cutoffs of 700bp for alignment length and 80% for percent identity (Quast et al., 2013). Abundance of unbinned 16S rRNA genes identified through this BLASTn was then determined through a BLASTn search of all reads against 16S rRNA gene containing portions of scaffolds (Quast et al., 2013). This was then compared to the bin abundances of corresponding taxonomy in order to match unbinned 16S rRNA genes to bins. All bins were further identified using conserved phylogenetic markers as annotated by IMG/M (Ciccarelli et al., 2006). The most resolved taxonomy with ≥50% of gene sequences matching to it was reported (except in the case of Bin A7, in which Rhodoferax ferrireducens was reported due to its co-binning with that reference genome). These bin identities were confirmed and resolved using PhyloSift database of 37 universal, conserved gene families (Darling et al., 2014). The percent identity reported for Phylosift analyses refers to the percent of scaffolds that were identified as that taxonomy in the PhyloSift database. Bin abundance was found by mapping all reads to scaffolds of each bin using BWA with default parameters in order to determine coverage (Li and Durbin, 2009). The coverage of each scaffold was normalized for scaffold length and averaged to get mean genome coverage for the bin, which was used as a measure of abundance. In addition, scaffolds from Bin A8 were mapped to a published Thiothrix lacustris genome (NCBI accession NZ_JHYQ01000054.1) (Chernousova et al., 2009) using the same method (Li and Durbin, 2009).

The metagenomes of both mat communities contained abundant members of known groups of LSB. The Isolated Sinkhole metagenome was dominated by two LSB genomic bins. One bin (IS2) represents a near-complete genome that was classified as Beggiatoa by both PhyloSift and 16S rRNA gene similarity (Table 1). The second (IS1) was classified by Phylosift as family Thiotrichaceae, but a fragment (1407 bp) of the 16S rRNA gene sequence had a closest match (89% sequence identity) to Thiomargarita namibiensis (Salman et al., 2011), which belongs to one of 8 currently recognized genera within the Beggiatoaceae. The 16S rRNA gene sequences associated with these two bins shared 84% sequence identity. The PhyloSift classification uses NCBI taxonomy, based on the system of Garrity et al. (2005), which proposed the Thiotrichaceae as a new polyphyletic family that encompasses the three traditional families of LSB, Leucotrichaceae, Achromatiaceae, and Beggiatoaceae. However, more recent work supports the traditional three family system (Salman et al., 2011, 2013) and thus would consider both bins IS1 and IS2 as members of the Beggiatoaceae.

One of the abundant genomic bins in the Alpena fountain (Bin A8) was closely related to T. lacustris (Table 1), an LSB genus known to form long filaments by attaching to surfaces in habitats with sulfidic flowing water (Larkin and Strohl, 1983), such as the Alpena fountain (Figure 1). This identification was supported by PhyloSift classification, by the co-binning of these scaffolds with the T. lacustris reference genome (Figure 2B), and by mapping of the reads to the T. lacustris genome at 157x average coverage. An unbinned 16S gene was found to have a 99% sequence similarity to T. lacustris (Chernousova et al., 2009), indicating it is the same species (Yarza et al., 2014).

In addition to LSB, the Isolated Sinkhole and Alpena fountain also had putative sulfate-reducing Deltaproteobacteria. Bins assigned to Desulfobacteraceae and Desulfobulbaceae were at moderate (Bin IS3; 13%) to low (Bins A2 and A6; 2% and 1%) relative abundances (Table 1). Members of these two families are strict anaerobes that reduce sulfate to sulfide (Kuever, 2014a,b). Other members of the microbial mat communities were different between the two sites. The Alpena fountain community is dominated by a well-defined genomic bin (A3; Figure 2) whose taxonomic identity could not be resolved by PhyloSift or IMG/M taxonomic profiling. Twenty-six percentage of its gene sequences matched to Archaea and 19% matched to Bacteria. Of the archaeal matches, the most common hits (24%) were to Euryarchaeota. An unbinned 16S gene with 97% sequence similarity to an uncultured SM1 euryarchaeon from cold sulfur springs that is highly divergent from cultivated archaea (Rudolph et al., 2001) was found to dominate the 16S rRNA gene sequence reads (Supplementary Figure S1) of the fountain, paralleling the dominance of genomic bin A3 of the fountain metagenome (Table 1). Members of this lineage, which was recently proposed to be a new euryarchaeal order, Candidatus Altiarchaeles (Probst et al., 2014b), form a ‘string of pearls’ morphology and are closely associated with Thiothrix filaments (Moissl et al., 2003) by attachment with nano-grappling structures (Moissl et al., 2005; Perras et al., 2014). Given the simplicity of the Alpena fountain community and the parallel dominance of the 16S gene and whole-genome scaffolds, this genomic bin can be confidently assigned as SM1 euryarchaeon.

Also present at high abundance in the Alpena fountain microbial mat was a genomic bin related to R. ferrireducens (Bin A7), a facultatively anaerobic betaproteobacterium that oxidizes organic carbon with the reduction of Fe(III) (Finneran et al., 2003). At lower abundance were two unknown groups of bacteria that PhyloSift analyses suggest are Verrucomicrobia (Bin1, 50% identity) and Elusimicrobia (Bin5, 87% identity) (Table 1). Also present at low abundance was a genomic bin for candidate division OD1 (Bin A4) (Harris et al., 2004), which has been found mainly in sulfur-rich, anoxic environments (Peura et al., 2012) and was recently proposed to be named Parcubacteria (Rinke et al., 2013).

To better identify the LSB, phylogenetic analysis was conducted on 16S rRNA and rpoB gene sequences assembled from shotgun metagenomic data. Initial alignments revealed the presence of two different large introns in two of the 16S rRNA genes (Figure 3). These introns have high sequence similarity to those previously identified in the marine LSB T. namibiensis (Salman et al., 2012). Such self-splicing introns are widespread among marine LSB, likely hinder PCR amplification of 16S rRNA genes, and may be horizontally transferred (Salman et al., 2012; MacGregor et al., 2013a; Flood et al., 2016). The introns found here are the first to be identified in non-marine LSB, and given their high sequence similarity to T. namibiensis, they are likely the result of horizontal gene transfer. After removal of the introns the remainder of the 16S rRNA gene sequence was most similar to other sequences from Isolated Sinkhole that may represent a novel group in the Beggiatoaceae (see below).

Phylogenetic analysis showed that the LSB sequences from Alpena Fountain and Isolated Sinkhole fell within two different clades (Figure 4). As expected, the Alpena fountain sequences clustered with Thiothrix sequences. The Isolated Sinkhole sequences formed novel clades within the Beggiatoaceae. Several branches were unresolved, consistent with previous phylogenetic analyses of these groups (Salman et al., 2013), but results suggest that the Isolated Sinkhole sequences represent new genera that are intermediate between the freshwater and marine Beggiatoaceae (Figure 4). Similar results were obtained from phylogenetic analysis of rpoB gene and protein sequences (Supplementary Figure S2). This phylogenetic placement might reflect the salinity of the groundwater fluids at the Alpena fountain and Isolated Sinkhole, which is ∼5% of seawater (Biddanda et al., 2012). However, the absence of such phylogenetically distinct LSB in other environments with comparable salinity, such as the Frasassi and Acquasanta cave systems (Macalady et al., 2008; Jones et al., 2010; Hamilton et al., 2015), suggests that other environmental factors may select for such organisms. Hereafter we refer to these organisms as Isolated Sinkhole Beggiatoaceae (ISB).

Genes for dissimilatory sulfur metabolism were prevalent in the metagenomes of microbial mats at both sites. LSB are characterized by their capability to oxidize reduced sulfur compounds for chemolithoautotrophic or mixotrophic growth (Teske and Nelson, 2006). As expected, sulfur oxidation pathways were encoded in the genome bins of Thiothrix (A8) and ISB (IS1 and IS2). The Thiothrix bin had genes for flavocytochrome c sulfide dehydrogenase (fcc), sulfide quinone oxidoreductase (sqr), the sox system (soxABXYZ) and the reverse dissimilatory sulfite reductase pathway (dsr) (Supplementary Table S2 and Figure 5), consistent with the ability of T. lacustris to use either sulfide, thiosulfate, or intracellular elemental sulfur globules as electron donors during lithotrophic growth (Chernousova et al., 2009). Bin A8 also contains genes for polysulfide reductase (psr), respiratory tetrathionate reductase and thiosulfate reductase (phs), which disproportionates thiosulfate into sulfide and sulfite during the respiration of thiosulfate (Akegawa et al., 1985; Stoffels et al., 2012). Thiosulfate reductase genes are also present in Beggiatoa genomes (Mußmann et al., 2007), and thiosulfate disproportionation is a key process in some environments with active sulfur cycling (Jørgensen and Bak, 1991). The presence of genes for both oxidation (sox) and respiration/disproportionation (phs) suggests that it may be a key intermediate in sulfur cycling in the Thiothrix mats studied here.

Only one of the two ISB bins from the Isolated Sinkhole (IS2) was near complete (Table 1). This genomic bin also contained genes for oxidation of sulfide (sqr and fcc), thiosulfate (soxABXY) and the reverse Dsr pathway (dsr, and also sat and apr) (Supplementary Table S2 and Figure 5), consistent with previously sequenced genomes of the Beggiatoaceae (Mußmann et al., 2007; MacGregor et al., 2013b; Flood et al., 2016; Winkel et al., 2016). It also had genes annotated as tetrathionate reductase. The availability of multiple pathways for oxidation of reduced sulfur species in the LSB at both the Alpena fountain and Isolated Sinkhole is consistent with a flexible lifestyle suited to a niche of varying availability of electron donors and acceptors (Klatt and Polerecky, 2015; Winkel et al., 2016).

In addition to lithotrophic growth on reduced sulfur compounds, genes for Ni-Fe hydrogenase in both the Thiothrix and ISB bins suggest that these organisms can also use hydrogen as an electron donor. These genes are commonly present in a wide range of organisms and habitats (Greening et al., 2016), including sulfur-oxidizing bacteria (Anantharaman et al., 2013). Members of the Beggiatoaceae are known to use H₂ as a supplemental energy source (Schmidt et al., 1987; Kreutzmann and Schulz-Vogt, 2016; Winkel et al., 2016). H₂ is produced by fermentation and may be present in terrestrial sulfidic springs (Aragno, 1992).

Both the Thiothrix and ISB genomic bins contained evidence of potential for aerobic respiration and nitrate reduction. Two different types of high affinity terminal oxidases for respiration at low-O₂ concentrations were identified: cytochrome cbb₃ oxidase and cytochrome bd oxidase (Morris and Schmidt, 2013). The Thiothrix bin (A8) contained both of these terminal oxidases, whereas the ISB bin contained only the cbb₃-type oxidase (Supplementary Table S2 and Figure 5). The presence of these high-affinity terminal oxidases suggest that LSB at each of these sites conduct aerobic respiration at low concentrations of O₂, which is likely introduced through diffusion and mixing with the atmosphere (Alpena fountain) or with overlying lake water (Isolated Sinkhole). These results are consistent with the microoxic lifestyle of many Beggiatoa spp. (Nelson et al., 1986b). To evaluate the possibility that these LSB are also capable of denitrification, we also searched for genes for reduction of nitrate, nitrite, and nitrous and nitric oxides (Supplementary Table S3; Philippot, 2002). Nar genes for dissimilatory nitrate reduction were found in the LSB genomic bins of both sites (ISB in the Isolated Sinkhole and T. lacustris in the fountain), as were norD and norQ genes for nitric oxide reduction. However, a homolog of the norB gene was only found in the Thiothrix bin (A8), and genes for dissimilatory reduction of nitrite or nitrous oxide were not found. The pathway for dissimilatory reduction of nitrate to ammonia (nirS and cytochrome c nitrite reductase), which is present in Thiomargarita along with genes for classical denitrification to N₂ (Flood et al., 2016), was also not present in the LSB genomes analyzed here. Thus, we find genomic evidence for only partial denitrification, including reduction of nitrate and potentially nitric oxide. The resulting accumulation of nitrite could necessitate the detoxifying nitrite-reducing enzymes in the sulfate-reducing bacteria (see below). In culture, the type strain for T. lacustris is capable of nitrate reduction (Chernousova et al., 2009) but anaerobic growth on nitrate only occurs in the presence of thiosulfate, not organic substrates (Trubitsyn et al., 2013).

Large sulfur bacteria including Thiothrix and Beggiatoa species have been shown to grow by autotrophy, heterotrophy, and mixotrophy (Nelson and Jannasch, 1983; Nelson et al., 1986a; Schmidt et al., 1987; Chernousova et al., 2009). In low oxygen, sulfidic environments, Beggiatoa prefer lithoautotrophic sulfide oxidation (Schmidt et al., 1987). We found genetic evidence of multiple carbon fixation pathways (Berg et al., 2010) in the LSB genomic bins recovered here (Supplementary Table S2 and Figure 5). The Thiothrix bin had genetic markers of the Calvin cycle (rbc genes) and reductive acetyl-CoA pathway (acyl-CoA synthetase). The reductive acetyl-CoA pathway contains O₂-sensitive enzymes and is only used by anaerobic or microaerophilic microbes (Hügler and Sievert, 2011). The IS2 ISB bin contained key genes for the 3-hydroxypropionate bicycle (malonyl-CoA reductase, propionyl-CoA synthase and acetyl-CoA carboxylase) and the Calvin cycle (rbc). We also found numerous genes annotated as organic carbon transporters in these bins, including transporters for acetate and dicarboxylates. Also present were genes for acetate kinase and acetyl CoA synthetase, which could channel acetate into the general metabolism and are found in Beggiatoa genomes (Mußmann et al., 2007; MacGregor et al., 2013b; Flood et al., 2016; Winkel et al., 2016). Genes for oxidation of glycolate, a common photosynthetic exudate, were also found, suggesting an interaction with cyanobacteria (Mußmann et al., 2007). Although no cyanobacteria were detected at Isolated Sinkhole, close relatives of these LSB are tightly associated with cyanobacteria at the Middle Island Sinkhole (Voorhies et al., 2012). Alternatively, glycolate oxidases have been implicated in the consumption of glycolate produced through the oxygenase activity of RubisCO in LSB (Winkel et al., 2016). Many glycolate-related genes were also found in the Thiothrix genomic bin (Supplementary Table S2). Overall, these results are consistent with the LSB being capable of both autotrophy and heterotrophy, and perhaps mixotrophy.

Abundant members of both communities were found to have pathway for dissimilatory sulfate reduction via the sat, apr, and dsr genes (Supplementary Table S2). This suggests that in addition to sulfide present in the groundwater as it emerges from the subsurface, sulfide for the LSB could be sourced from the bacterial reduction of sulfate and/or intermediate sulfur species within the mat (Figure 5). The genomic bins that contained genes for sulfate reduction were from Deltaproteobacteria groups well known to reduce sulfate. In the Alpena fountain, the genome bin (A2) identified as Desulfobulbaceae was the only genomic bin with the entire pathway for reduction of sulfate to sulfide. In addition to producing sulfide via the Dsr pathway, genes for tetrathionate reductase and thiosulfate reductase suggest that sulfide can also be produced via alternative routes involving the reduction and disproportionation of sulfur intermediates (Figure 5). It also contained genes for the use of O₂ as an electron acceptor, H₂ as an electron donor, a wide variety of organic carbon substrates as a carbon source (Supplementary Table S2). This bin also had genes potentially encoding the reductive acetyl-CoA pathway (Wood-Ljungdahl) for carbon fixation (Supplementary Table S2). In addition to using this pathway for carbon fixation, it can also be run in reverse and coupled to sulfate reduction to generate energy (Ragsdale and Pierce, 2008). A partial bin of Desulfobacteraceae in the Alpena fountain (A6) also had the genes for sulfite and O₂ reduction, H₂ oxidation, carbon fixation, and heterotrophy, but it did not have the sat or apr genes for sulfate reduction. However, a near-complete genomic bin of Desulfobacteraceae in the Isolated Sinkhole (IS3) contained the complete pathway along with genes for disproportionation of thiosulfate and the same complement of genes for H₂ oxidation, O₂ reduction, carbon fixation, and heterotrophy (Figure 5). This bin also contained a 16S rRNA gene sequence that was 97% similar to a sequence from an uncultured deltaproteobacterium of the Frasassi Cave biofilms. This sequence clusters phylogenetically with Desulfonema and appears to be filamentous and to have close physical associations with Beggiatoa and Thiothrix (Macalady et al., 2006). Desulfobacteraceae and Desulfobulbaceae have also been observed at the Middle Island Sinkhole, in mat consortia of cyanobacteria and LSB (Nold et al., 2010; Kinsman-Costello et al., 2017). Both Deltaproteobacteria bins from the Alpena fountain (A2 and A6) had genes for cytochrome c nitrite reductase, which reduces nitrite to ammonium and is used by sulfate-reducing bacteria for detoxification of nitrite, often in association with sulfide-oxidizing, nitrate-reducing bacteria (Haveman et al., 2004). Bin IS3 also had the complete pathway for nitrate reduction to ammonium via the reverse hydroxylamine:ubiquinone reductase module pathway (Supplementary Table S2; Hanson et al., 2013).

The genome-predicted phenotypes of H₂-based autotrophy, organoheterotrophy, respiration of O₂ have all been observed in cultured sulfate-reducing bacteria including Desulfobulbaceae and Desulfobacteraceae (Kuever, 2014a,b). Ni-Fe hydrogenase oxidizes H₂, which is an important electron donor for anaerobic sulfate-reducing bacteria such as the deltaproteobacterium Desulfovibrio (Eidsness et al., 1989). Genes for cytochrome oxidases may indicate that the sulfate-reducing bacteria can use oxygen as a terminal electron acceptor for carbon metabolism (Canfield and Des Marais, 1991). Generally this process is linked to ATP production, but not growth (Dilling and Cypionka, 1990; Dannenberg et al., 1992; Sass et al., 2002), although weak growth has been observed (Marschall et al., 1993; Sigalevich and Cohen, 2000).

Whereas the Isolated Sinkhole metagenome was dominated by just the two LSB bins and one sulfate-reducing bacterium bin, the Alpena fountain featured additional near-complete genomes from other community members. The genomic bin from the SM1 euryarchaeon accounted for 71% of all mapped metagenomic reads from the fountain (Table 1). This lineage is widespread is sulfide springs (Rudolph et al., 2004) but its physiology and potential metabolism have been enigmatic (Probst et al., 2013). Based on close physical associations, syntrophic relationships with Thiothrix (Moissl et al., 2002) and sulfate reducing bacteria (Probst et al., 2013, 2014a) have been proposed. Recent genomic data derived from SM1 populations in flowing, oxygen-depleted groundwater in Germany and stratified suboxic/anoxic groundwater in the USA indicate that SM1 can fix carbon via a novel reductive acetyl-CoA (Wood-Ljungdahl) pathway, but no clear genetic evidence of electron donors or acceptors for energy metabolism was uncovered (Probst et al., 2014a). The Alpena fountain SM1 population also has key genes for the reductive acetyl-CoA pathway for carbon fixation, including acetyl-coenzyme A synthetase and CO dehydrogenase/acetyl-CoA synthease (Supplementary Table S2). We also found genes for a Ni,Fe-hydrogenase, suggesting an electron donor for lithotrophy, and genes putatively annotated as transporters for amino acids, polysaccharides, and C4-dicarboxylates, suggesting potential for heterotrophy as well (Figure 5 and Supplementary Table S2). Experimental studies are required to confirm these genome-generated hypotheses.

The next most abundant bin in the Alpena fountain was A7, which accounted for 14% of mapped metagenomic reads and was identified as a Rhodoferax species based on PhyloSift, 16S rRNA gene sequence, and co-binning with the R. ferrireducens genome (Figure 2). Consistent with the known metabolic versatility of R. ferrireducens (Risso et al., 2009; Yao et al., 2013), we found genes for transporters of organic compounds, sulfide and thiosulfate oxidation, aerobic respiration, dissimilatory nitrate reduction, and autotrophy (Supplementary Table S2). Other genomic bins (A1, A4, A5) were present at low relative abundance and the few genes of functional significance noted were related to uptake of organic carbon compounds (Supplementary Table S2).