Significant Association between Sulfate-Reducing Bacteria and Uranium-Reducing Microbial Communities as Revealed by a Combined Massively Parallel Sequencing-Indicator Species Approach

Massively parallel sequencing has provided a more affordable and high-throughput method to study microbial communities, although it has mostly been used in an exploratory fashion. We combined pyrosequencing with a strict indicator species statistical analysis to test if bacteria specifically responded to ethanol injection that successfully promoted dissimilatory uranium(VI) reduction in the subsurface of a uranium contamination plume at the Oak Ridge Field Research Center in Tennessee. Remediation was achieved with a hydraulic flow control consisting of an inner loop, where ethanol was injected, and an outer loop for flow-field protection. This strategy reduced uranium concentrations in groundwater to levels below 0.126 μM and created geochemical gradients in electron donors from the inner-loop injection well toward the outer loop and downgradient flow path. Our analysis with 15 sediment samples from the entire test area found significant indicator species that showed a high degree of adaptation to the three different hydrochemical-created conditions. Castellaniella and Rhodanobacter characterized areas with low pH, heavy metals, and low bioactivity, while sulfate-, Fe(III)-, and U(VI)-reducing bacteria (Desulfovibrio, Anaeromyxobacter, and Desulfosporosinus) were indicators of areas where U(VI) reduction occurred. The abundance of these bacteria, as well as the Fe(III) and U(VI) reducer Geobacter, correlated with the hydraulic connectivity to the substrate injection site, suggesting that the selected populations were a direct response to electron donor addition by the groundwater flow path. A false-discovery-rate approach was implemented to discard false-positive results by chance, given the large amount of data compared.Massively parallel sequencing has increased our ability to study microbial communities to a greater depth and at decreased sequencing costs to an extent that replication and gradient interrogation are now reasonably attainable. However, this massive throughput has mostly been used in exploratory studies, given the challenges to analysis of the big data sets generated and the relative novelty of the technique. To date, no report of a study that has used this method to describe the microbial community over a large area influenced by complicated hydrogeochemical factors during bioremediation has been published. Here, we used pyrosequencing technology complemented with a hypothesis-based approach to identify bacteria associated with biostimulation of U(VI) reduction at Area 3 of the U.S. Department of Energy''s (DOE''s) Oak Ridge Field Research Center (FRC) at Oak Ridge, TN.The Oak Ridge FRC is one of the most-studied sites for uranium bioremediation (2, 8, 19-22, 27, 37, 45-48). Previously used as a uranium enrichment plant, the site remains contaminated with depleted uranium, nitrate, and acidity. To deal with uranium contamination, dissimilatory metal reduction has been studied as an alternative that reduces risk by converting toxic soluble metals and radionuclides to insoluble, less toxic forms (2, 3, 16, 21, 26, 45). For example, some microbes can use metals such as Cr(VI), Se(VI), and the radionuclides U(VI) and Tc(VII) as final electron acceptors, producing a reduced insoluble species, thus blocking dispersal and reducing bioavailability.The ability to reduce U(VI) to U(IV) has been found in several unrelated phylogenetic groups, i.e., Delta-, Beta-, and Gammaproteobacteria, Firmicutes, Deinococci, and Actinobacteria, among others (42). Most previous studies have focused on the Fe(III)-reducing bacteria (FRB), especially Geobacter, and the sulfate-reducing bacteria (SRB), especially Desulfovibrio. Uranium(VI) reduction for bioremediation purposes has been tested and confirmed in laboratory-scale experiments using serum bottles (13, 18, 48), microcosms (23, 32), sediment columns (14, 43), and in situ field studies (3, 21, 41, 45), with the last one demonstrating the feasibility of U(VI) remediation and the correlation of U(VI) reduction with FRB (3, 6, 18, 31, 41) or SRB (40), or both (8, 19, 49).During field studies at Area 3 of the Oak Ridge site, a hydraulic control system together with ethanol injection successfully promoted U(VI) reduction from 5 μM to levels below U.S. Environmental Protection Agency (EPA) maximum contaminant levels (MCLs) for drinking water (0.126 μM) over a 2-year period (46). Reduction of U(VI) to U(IV) was confirmed by X-ray absorption near edge structure (XANES) (22, 46). Previous microbial surveys of sediments and groundwater from Area 3 wells by the use of 16S rRNA gene clone libraries detected genera known to harbor U(VI)-reducing members, such as Geobacter, Desulfovibrio, Anaeromyxobacter, Desulfosporosinus, and Acidovorax, after U(VI) reduction was established (8, 19). In one study, microbial counts from sediments were correlated with the hydraulic path, suggesting differences in organic carbon availability throughout Area 3 (8). The study that tracked the groundwater microbial communities of four locations of Area 3 over a 1.5-year period during ethanol stimulation found that nitrate, uranium, sulfide, and ethanol were correlated with particular bacterial populations and that the engineering control of dissolved oxygen and delivered nutrients was also significant in explaining the microbial community variability (19). However, the analysis of communities has been focused on limited wells and the community of the entire test area has not been characterized.On the basis of the previous results, we further hypothesized that the hydrological control strategy employed for the remediation of the site constrained the geochemistry of the site by controlling the distribution of organic carbon substrates and other nutrients and that this in turn selected a characteristic microbial community that was distinguishable from its surrounding community. We used massively parallel sequencing of 16S rRNA genes from sediments of 15 wells to characterize the microbial communities along hydrological gradients from the microbiologically active and hydraulically protected inner-loop zone to less active and still contaminated areas outside the treatment area and downgradient. Our sediment-sampling strategy allows a more precise spatial characterization than the use of groundwater samples, where filtering large volumes of water is often required, and also because samples of the attached communities can differ from the planktonic ones, as expected in oligotrophic aquifers (15), such as this site. The deeper sequencing allowed a more extensive survey of the communities, higher confidence in the detection of less dominant but significant members, and a more statistically robust indicator species assessment. We were able to detect groups significantly associated with U(VI) reduction and to explain differences in community structure with hydrogeochemical conditions.