Key Factors Controlling Microbial Distribution on a DNAPL Source Area

doi:10.21203/rs.3.rs-400510/v1

Download PDF

Research Article

Key Factors Controlling Microbial Distribution on a DNAPL Source Area

https://doi.org/10.21203/rs.3.rs-400510/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 05 Aug, 2021

Read the published version in Environmental Science and Pollution Research →

You are reading this latest preprint version

The chlorinated solvents are among the common groundwater contaminants that show high complexity in their distribution in the subsoil. Microorganisms play a vital role in the natural attenuation of chlorinated solvents. Up to now, how the in situ soil microbial community responds to chlorinated solvents contamination remains unclear. In this study, the microbial community distribution within two boreholes located in the source area of perchloroethene (PCE) were investigated via terminal restriction fragment length polymorphism (T-RFLP) and clone library analysis. Microbial data were related to the lithological and geochemical data and the concentration and isotopic composition of chloroethenes to stablish the key factors controlling the distribution of the microbial communities. The results indicated that Proteobacteria, Actinobacteria and Firmicutes were the most abundant phylums in the sediment. The statistical correlation with the environmental data proved that fine granulometry, oxygen tolerance, terminal electron acceptor processes and toxicity control microbial structure. This study improves our understanding of how the microbial community in the subsoil responds to high concentration of chlorinated solvents.

Environmental Engineering

Environmental Policy

T-RFLP

toxicity

DNAPL

microbial heterogeneity

perchloroethene

The subsoil is a heterogeneous medium with many fluctuating parameters that affect the growth and survival of microorganisms (Paul and Clark 1989). Geological factors, like pore size and the interconnectivity of the sediments (Puigserver et al. 2020; Mahmoudi et al. 2020), and the biogeochemical composition of the sediments, like the content and composition of the organic matter and the presence of metals (Van Horn et al. 2013), control the variation in the microbial communities across the subsurface environments. In saturated zones, the velocity, hydrochemical parameters and temperature of the groundwater, as well as the composition of the planktonic microbial communities, also affect structure of the microbial communities along the sediment (Velasco Ayuso et al. 2009; Guo et al. 2019). In pollutant episodes, the distribution of contaminants and their daughter products is controlled by geological, hydrogeological and biogeochemical parameters (Guilbeault et al. 2005; Hartog et al. 2010; Puigserver et al. 2013), and this, in turn, conditions the microbial communities (Griebler and Lueders, 2009; Schmidt et al., 2017). Specifically, the type and distribution of contaminants, and their toxicity and biodegradability are key factors that explain changes in the structure of the microbial communities (Rossi et al. 2012; Puigserver et al. 2013, 2020).

Chlorinated solvents, which are types of DNAPL (dense non-aqueous phase liquids), are among the common groundwater contaminants that show high complexity in their distribution in the subsoil and affect the composition of the microbial communities (Puigserver et al. 2020). Chlorinated contaminants are accidently released into the environment at many industrial and urban sites worldwide and are highly toxic (He et al. 2015). They migrate as a free phase through the porosity, are heterogeneously distributed as pools and a residual phase between pores, and are sorbed in organic matter and fine materials due to molecular diffusion (Parker et al. 2003; Chapman and Parker 2005). Together, these migrated compounds form the source zone (Parker et al. 2003). The compounds dissolved in groundwater create large contamination plumes and can be volatilised, remaining within the gaseous matrix or dissolving again in the water (Mackay et al. 2006). The morphology of the contamination source areas is what conditions the plume of volatilised and dissolved contaminants (Pankow and Cherry 1996). In remediation strategies for chlorinated solvents, the source zone is treated in the first stage, and the plume is treated afterwards (or at the same time).Chemical and physical strategies have mostly been applied as source zone remediation strategies, rather than the biological strategies that are mainly used in the plume (Stroo 2010; Stroo et al. 2012). In recent years, however, studies have highlighted the potential for applying bioremediation strategies in source zones as well (Herrero et al., 2019; Sung and Ritalahti, 2003; Yang and McCarty, 2000, 2002).

Under anoxic conditions, chlorinated solvents such as chloroethenes are mainly degrade via reductive dehalogenation, which involves the sequential reduction of these compounds (Smidt and de Vos 2004). This process requires increasingly reductive conditions (Wiedemeier et al. 1999; Bradley 2003), and methanogenic conditions for their complete reductive dehalogenation to an inert compound (Hata et al. 2004). Organohalide-respiring bacteria (OHRB) are mainly responsible for the reductive dehalogenation (Adrian and Löffler 2016). OHRB are usually a small percentage of the total bacterial community, compared to the fermenting and sulphate-reducing populations, that are in much greater abundance (Ndon et al. 2000; He et al. 2002; Men et al. 2012).

The complex interactions between dehalogenating microorganisms and the structure of the microbial community have been of interest during the last two decades (e.g. Atashgahi et al., 2017; Balaban et al., 2019; Dojka et al., 1998; Fennell et al., 2001; Flynn et al., 2000; Freeborn et al., 2005; Hendrickson et al., 2002; Hohnstock-Ashe et al., 2001; Rossi et al., 2012, 2009). In addition, a better understanding of the relationship between the structure of the microbial community and the dehalogenators will lead to the development and optimisation of bioremediation strategies. Of particular interest are the contact areas between two different geological units, as shifts in microbial communities and biogeochemical processes are expected in these areas (McMahon and Chapelle 2008; Puigserver et al. 2013, 2016; Griebler and Avramov 2015)

A perchloroethene (PCE) source area in a site of alluvial fans whit a highly heterogeneous geological, biogeochemical and contaminant distribution was chosen in order to determine the main factors that condition the structure of the microbial communities in the source zone of chlorinated solvents. The goal was to identify the factors that affect the structure of the microbial communities, enabling a more detailed definition of the conceptual model of a contamination episode. The identification of these factors has the potential to improve the efficiency of bioremediation strategies. The studied factors are related to the granulometry, biogeochemical processes, and distribution of chloroethenes in the sediments, and to the hydrochemistry of the aquifer. The specific objectives of the research were: 1) to characterise the microbial distribution of two boreholes in the source area, particularly sampling the contact areas, where microbial shifts are expected; and 2) to assess the main factors that affect the composition of the microbial communities.

2.1. Site description, core sampling and conservation protocol

The field site is in an industrial area in Vilafant (Alt Empordà, NE Spain), approximately 150 km to the north of Barcelona. The aquifer consists of Pliocene prograding alluvial fan deposits, and a PCE-DNAPL source was detected by the Catalan Water Agency (ACA) in the transition zone to a basal aquitard.

The drilling method and the general sampling and conservation protocol are described by Puigserver et al., (2016). The core sampling was performed take into consideration the lithological and textural changes and following the criteria indicated by Guilbeault et al., (2005) and Puigserver et al. (2013). A total of 60 samples were taken from the F1UB borehole and 115 samples from the F2UB borehole. These samples were taken from the central part of the borehole with sterile tools; they were placed inside a sterile container with distilled water, and were immediately frozen to below − 20°C.

A total of 29 samples, 15 from the F1UB borehole and 14 from the F2UB borehole, were selected for molecular analysis. The selection criteria were based on the detailed geological characterisation, the concentration of chloroethenes in the porewater, and the concentration of organic carbon, iron and manganese in the sediment.

The groundwater sampling of the two multilevel wells located in the F1UB and F2UB boreholes, and the hydrochemical analyses, are described by Herrero et al. (2021).

2.2. Environmental data analysis and treatment

The environmental data used in the correlation with the microbial data are the particle size of the sediment, the concentration of organic carbon, Fe, Mn and chloroethenes, and the isotopic composition of the PCE. The content of organic carbon (Corg), Fe and Mn sorbed in the fine fraction of sediments, and the chloroethene concentration analysis and calculations in the porewater are described by Puigserver et al. (2016). The analysis of the isotopic composition of the PCE in the porewater is described by Herrero et al. (2021b).

To assess the presence or absence of the process of reductive dehalogenation, a qualitative variable was developed. It was determined that if a daughter product of PCE and/or isotopically enriched PCE in the porewater was present, reductive dehalogenation was occurring. To assess the toxicity, a new variable of the sum of all chloroethenes was used.

2.3. Molecular analysis

The analyses were performed at the Helmholtz Centre for Environmental Research-UFZ (Leipzig, Germany). Genomic DNA was extracted from 1.1 g of sediment with the NucleoSpin® Soil of Macherey & Nagel, following the manufacturer’s protocol, to perform terminal restriction fragment length polymorphism (T-RFLP) and clone library analysis.

Polymerase chain reaction (PCR) was performed according to Puigserver et al. (2016). The PCR product was purified using the purification Kit Wizard® for Genomic DNA (Promega). A total of 50 ng of purified DNA was restricted with three different restriction enzymes (HaeIII, HhaI and MspI, Thermo Scientific). The dry DNA was dissolved with Hi-DiTM Formamid (Applied Biosystems) with the standard GeneScan™ 500 ROX™, and was analysed with the ABI 3100 Genetic Analyser (Applied Biosystems) and Genemapper 3.7 Software (Applied Biosystems). Duplicates of each sample were analysed, and consequently, six results were obtained for each sample.

2.4. Data treatment

To validate the results, all the duplicates for each restriction enzyme were checked and the test was repeated if the results were not conclusive. All restriction fragments (RFs) smaller than 50 bp and having a proportion smaller than 1% of the total area were eliminated.

The microbial richness of each sample was determined by the number of valid RFs for each restriction enzyme. The degree of development was determined by the total area of the valid RFs. Quantification of the population density using the T-RFLP technique allows a comparison of the degree of development for each of the samples. Several authors (Bruce 1997; Liu et al. 1997) recommend treating this measure as a semiquantitative value, as T-RFLP analysis is subject to all the biases inherent in any PCR approach. Following these directions, the average of the areas of the different restriction enzymes was averaged again and transformed on a scale of 1 to 10.

The results obtained using T-RFLP allowed the degree of similarity to be established from a cluster analysis. Ward’s algorithm was chosen, as it considers the peaks and the percentage area of each (Murtagh and Legendre 2011). This method is based on integrating the different individuals (in this case, microbial communities) into clusters producing the minimum difference, in terms of the percentage area and the number of peaks. The method raises all possible fusions at each stage and selects the one that maximises homogeneity (or minimises heterogeneity): it calculates the centroids of the groups resulting from the possible fusions, calculates the distance to the centroid of all group observations, and chooses the solution with the smaller total quadratic sum.

The correlation between the most abundant and common RFs for the different restriction enzymes was determined by statistical treatment. Next, the phylogenetic assignment tool (PAT, Kent et al. (2003)), was used to assign potential taxa to each set of three RF signals. The PAT database was amplified with the data from the In silico project (Bikandi et al. 2004). The reported data of RFs are the means of the results for the three restriction enzymes.

2.5 Conceptual model

The PCE source zone has been characterised by Puigserver et al. (2016) as having five hydrostratigraphic units (the unsaturated zone (UZ), the discontinuous confining aquitard (UDTA), the upper aquifer (UPA), the transition zone to the basal aquitard (TZBA) and the basal aquitard (BA)). In addition, the distribution of chloroethenes, Fe, Mn, Corg and the richness of microbial communities have been characterised (Puigserver et al. 2016). Herrero et al. (2021b) present a new compound-specific isotope analysis (CSIA) method for chlorinated solvents in the porewater, applied to the three saturated hydrostratigraphic units (UPA, TZBA and BA).

The UZ is composed of coarse, medium and fine gravel and sand, with a silty-clayey matrix exceeding 50%. A heterogeneous distribution of PCE was detected in this unit, with a pronounced maximum of 10,385 µg/L in the porewater in the F1UB borehole, and evidence of dehalogenation above this. The UDTA consists mainly of clays, with a network of subvertical microfractures. Increased PCE was detected in the porewater in comparison with the upper unit (UZ) and the lower unit (UPA), with concentrations around 300 µg/L. The UPA is composed of gravel and coarse sand with about 15% of the levels having a silty-clayey matrix. PCE concentrations in this zone are low, with the exception of one silty-clayey matrix level which shows a concentration of 1150 µg/L in the porewater in the F1UB borehole, and there is evidence of dehalogenation (enriched δ¹³C_PCE in the porewater). The groundwater presents oxic conditions, with dissolved oxygen concentrations of around 10 mg/L. There is also evidence of denitrification, Mn reduction and reductive dehalogenation in the upper part of the UPA (Herrero et al. 2021a). The TZBA is made up of gravel and coarse sand alternating with numerous layers of medium to fine sand and silt on a centimetre to decimetre scale of limited horizontal extension. The F1UB borehole has more levels with a clayey-silty matrix (about 90%) than the F2UB borehole, which has about 30%. A residual pool of PCE was found in the TZBA contact with the BA, with maximum concentrations in the porewater of 18,175 and 6,409 µg/L of PCE in F1UB and F2UB, respectively. Reductive dehalogenation was found to be active above this maximum PCE. The groundwater presents more reductive redox conditions as it becomes deeper, and denitrification, Mn and Fe reduction, sulphate reduction and reductive dehalogenation of PCE and trichloroethene (TCE) were detected. No anaerobic processes were detected in the upper zone of the TZBA, and the conditions remain oxic, with a concentration of dissolved oxygen around 8 mg/L (Herrero et al. 2021a). The BA is composed of fine sands and laminar silts that are microfractured (with subvertical fractures). The distribution of PCE and TCE within the BA is ruled by the presence of vertical microfractures and stratification planes in the very fine sands with a silty-clayey matrix (Herrero et al. 2021b).

Microbial richness and the semiquantitative measure of degree of development indicate hot spots of biogeochemical activity, as well as areas with the inhibition and/or specialisation of some microbial populations. Figure 1 shows the variation of richness and degree of development in depth and related to the hydrostratigraphic units and the distribution of DNAPL. There is a clear correlation between richness and degree of development in the F1UB borehole, while this correlation is more diffuse in F2UB. The richness and degree of development show relative maximums in the joint points of the BA (12BA-F1 and 14BA-F1) and in the TZBA above the residual pool (9TZBA-F1 and 6TZBA-F2), while the richness shows relative maximums unrelated to the degree of development in the contact areas between the two different hydrostratigraphic units (1UZ-F2 and 3UPA-F2). Generally, the decrease of richness and degree of development is related to the presence of a high amount of PCE (2UZ-F1, 10TZBA-F1, 11TZBA-F1, 7TZBA-F2, 8BA-F2 and 9BA-F2) and with the hydrostratigraphic units of finer material (5UDTA-F1, 13BA-F1, 15BA-F1, 2UDTA-F2 and 11BA-F2).

3.1 Comparison of microbial community structure

The T-RFLP results for the three restriction enzymes of the boreholes F1UB (Figure SI 1) and F2UB (Figure SI 2) were highly complex. Within the 29 sediment samples, 40, 70 and 37 different RFs higher than 1% of the total area were detected for MspI, HhaI and HaeIII, respectively (Figures SI 1 and 2). For the sediment and groundwater samples, 47, 77 and 45 different RFs for MspI, HhaI and HaeIII, respectively (Figures SI 1, 2 and 3) were detected, with at least eight new RFs in the groundwater. The statistical of cluster analysis treatment (Figures SI 4, 5 and 6) showed no conclusive results, as the samples were grouped in different clusters for each enzyme. Specifically, the microbial communities with higher richness were located in different clusters for each restriction enzyme. The distribution of the samples within each cluster proved that there were some similarities between communities, but none could be explained by its location or a single set of environmental variables. The environmental data tested included the total amount of chloroethene, PCE, TCE, cis dichloroethane (cisDCE), Corg, Mn and Fe, the depth and predominant lithology. There was no dominant variable controlling the composition and distribution of the microbial communities, but there was a group of variables that differed in importance depending on the location. Microbial communities in the presence of contaminants developed in a more complex way given the increased heterogeneity of the medium. Communities developed and had different metabolism depending on: 1) the characteristics of the surface to which they were attached; 2) the balance of nutrients and contaminants between the solid, liquid and gaseous phases; and 3) the concentration of the pollutants. In these cases, we are talking about heterogeneity on a centimetre scale.

A comparative analysis of the three restriction enzymes of each sample allowed the fingerprints of the 10 most abundant populations to be determined (Table 1). These ten populations were selected as they are the dominant population in at least one of the characterized communities. Each set of three RFs was analysed using the philogenetic assignment tool (PAT, Kent et al., 2003) and double checked with the results of the clone library (Puigserver et al. 2016).

RF1 was identified by PAT and clone library as Propionibacterium acnes, an anaerobic microorganism that produces propionic acid by fermentation (Green 1992), and is also related to the reductive dehalogenation of PCE and TCE (Chang et al. 2011; Moreno et al. 2011). This RF was found almost ubiquitously in the whole study area, in sediments and groundwater, and especially in the contact areas of the different hydrostratigraphic units and the upper and lower levels of the pool of PCE (Fig. 2A and B). RF2 by PAT was identified as Acidithiobacillus ferrooxidans, a facultative aerobic organism capable of reducing Fe³⁺ (Ohmura et al. 2002). This RF was distributed heterogeneously along the two boreholes, although it was related to the most oxidant conditions of the upper part of the aquifer and of the unsaturated zone of F1UB (Fig. 2A) and to the groundwater of the upper and lower part of the aquifer (Fig. 2B). RF3 was identified by PAT as Streptomyces sp. or Arthrobacter sp.. This RF was distributed homogeneously in the UDTA and UPA in the F2UB borehole (Fig. 2A). RF4 was identified by PAT as Streptococcus sp. and was positively identified by the clone library as Aerococcus viridans. The Aerococcus genera is microaerophilic (Vela et al. 2007) and autochthonous to groundwater (Cruz-Perez et al. 1996). There was a high proportion of RF4 in the UDTA, the TZBA and the BA, while it was practically absent from the UZ and in the UPA, and it increased in the TZBA (Fig. 2A). Also, RF4 was found in the groundwater at the centre of the aquifer (Fig. 2.B). RF5 was identified by PAT as Staphylococcus sp. and was positively identified by the clone library as Aeribacillus pallidus. RF5 was found in the UZ, the UDTA and the UPA of F1UB, and in the BA of F2UB (Fig. 2A), and in the groundwater of the TZBA of F2UB (Fig. 2B). RF6 was identified by PAT as Microbacterium sp. or Terrabacter sp. and was found at the base of the TZBA and the BA of F1UB (Fig. 2A) and in the groundwater of the central part of the aquifer in F1UB (Fig. 2B). RF7 was identified by PAT and the clone library as Acinetobacter junii, an aerobic bacterium that is found ubiquitously in the soil and water and is able to degrade a wide variety of organic compounds (Towner 1992). RF7 was found mainly in the UZ of F1UB, and in the interphase of the UDTA and the UPA and the upper part of the BA of F2UB (Fig. 2A). RF8 was identified by PAT as Haemophilus sp.. RF8 was detected at the interphase between the UDTA and the UPA and in the upper part of the BA of the F2UB borehole and in F1UB at the bottom part of the BA (Fig. 2A). RF9 was not identified by PAT and was located mainly in the UDTA of F1UB and the BA of both boreholes (Fig. 2A). RF10 was identified by PAT and the clone library as Variovorax paradoxus, an aerobic bacterium related to oxidative dehalogenation (Futamata et al. 2005; Humphries et al. 2005). RF10 was located above the peak of PCE of the UZ and the TZBA of F1UB and in the UPA and the BA of F2UB (Fig. 2A).

The RFs (Fig. 2) were more easily connected to the environmental data, as it was possible to identify the main factors determining the distribution of the microbial populations, and therefore, the structure of the microbial communities. These factors were grouped into four groups: geological factors (majority granulometry, percentage of fines), hydrogeological factors (capacity to be transported in an aqueous medium), terminal electron-accepting processes (TEAP, e.g. Corg, Mn, Fe, metabolism of the identified populations), and conditioning factors due to the presence of contamination (concentration of PCE and evidence of its degradation).

Table 1

Most abundant microbial populations, quantified by the restriction fragments (RFs) and identified by the phylogenetic assignment tool (PAT). +: identified by clone library and PAT. *identical fingerprint to bacteria in anaerobic fermentation reactor (GU454879.1.1495), microbial biofilm (DQ499314.1.1492), and groundwater contaminated with nitric acid bearing uranium waste (AY662046.1.1527), among others.
RF	HaeIII	HhaI	MspI	Bacteria	Phylum
1	62.5	675	165	Propionibacterium acnes (+)	Actinobacteria
2	251	205	485	Acidithiobacillus ferrooxidans	γ-Proteobacteria
3	226	468	160	Streptomyces, Arthrobacter	Actinobacteria
4	308	585	560	Streptococcus, Aerococcus viridans (+)	Firmicutes (Bacilli)
5	308	236	153	Aeribacillus pallidus (+), Staphylococcus sp.	Firmicutes (Bacilli)
6	230	143	279	Microbacterium sp., Terrabacter sp.	Actinobacteria
7	253	207	491	Acinetobacter junii (+)	γ-Proteobacteria
8	204	363	491	Haemophilus sp.	γ-Proteobacteria
9	196	204	140	Uncultered bacterium*.
10	217	62	485	Variovorax paradoxus (+)	β- Proteobacteria

3.2 Geological factors

There is a bivariate correlation between the distributions of fine materials, from fine sand to clay (diameter less than 0.25 mm), and certain microorganisms (Fig. 3A, B and C). The RF4, RF6 and RF9 populations were mainly found in the hydrostratigraphic units with more fine materials (UDTA and BA), and in the UPA and TZBA levels with more fine materials (Fig. 2). In fact, RF4 and RF9 were not detected in samples with less than 40% of fine materials (Fig. 3A and C), and RF6 was only found in the levels with a minimum of 80% of fine materials (Fig. 3B). Other microorganisms, such as RF1, RF2 and RF5, do not show a dependence on sediment granulometry and were distributed throughout the different hydrostratigraphic units (Fig. 2).

The absence or very low proportion of specific microbial populations (RF4, RF6 and RF9) in the UZ, the UPA and the coarser levels of the TZBA, compared to a higher proportion in the finer particle size levels such as the UDTA and the BA, can be explained in several ways. On the one hand, the finer levels are those with a higher proportion of organic carbon, and the populations may therefore be adapted to its degradation. Another explanation could be the non-dependence of the nutrient supply (bioavailability) on saturated sections that are hydraulically more conductive. This would mean that there are populations more capable of taking advantage of the contributions of groundwater. Another explanation may be that these populations are not adapted to changes in the physical parameters (such as temperature) or the hydrochemical parameters (such as dissolved oxygen, dissolved organic matter, redox potential, phosphates, nitrates, among others) of the groundwater.

On the other hand, all populations that were identified in the saturated zones (UPA and TZBA) were found in the BA. It can therefore be concluded that pore size is not a limiting factor in the distribution of the majority populations. This differs from other studies (e.g. Puigserver et al., 2020) that found that, pore size limits the colonisation of some bacteria at the finest levels.

3.3 Hydrogeological factors

The ability of bacteria to colonise sediment through the flow of groundwater is another factor that explains the distribution of microbial populations. As can be seen in Figs. 2A and B, RF3, RF7, RF8 and RF10 were only found in sediments and not in groundwater. In the case of RF7 and RF8 this can be explained by their low presence in sediments in the UPA and TZBA, and RF3 and RF10 may not be able to survive in planktonic form or to form floccules.

The presence of RF1, RF2, and RF5 in the groundwater seems to demonstrate that these microorganisms can colonise other areas of the aquifer, either as floccules, planktonic cells, or attached to clays or silts (Griebler and Lueders 2009). These populations are also related to active biogeochemical processes (denitrification and reduction of Mn and Fe), as they were found in the upper and lower part of the aquifer. These two zones have been defined as ecotones by Herrero et al. (2021).

The presence of RF4, RF6 and RF9 in the groundwater may be related to whether these bacteria are attached to clays or silts in suspension in the groundwater. On the one hand, these populations are related to the fine materials (previous section) and, on the other hand, they were mostly detected in the centre of the aquifer, where no biogeochemical process was detected at a hydrochemical level (Herrero et al. 2021a). In relation to this, Zhao et al. (2012) show that Streptococcus (RF4) is able to adhere and travel in clay size particles.

3.4 Oxygen Tolerance

The tolerance of microorganisms to fluctuating oxygen levels is a limiting factor. The aquifer (UPA and TZBA) has dissolved oxygen concentrations that vary in depth and time (from 12.30 to 0.12 mg/L). Although the medium is generally oxic, there are micro-niches with gradations of oxygen concentration and redox conditions on a millimetre scale, which allow anaerobic microorganisms to metabolise (Rivett et al. 2008; Perović et al. 2017). These gradations are more important when there is more geological heterogeneity, as is the case in the TZBA compared to the UPA (Puigserver et al. 2016). In fact, denitrification and the reduction of Mn were detected in the upper part of the UPA, and denitrification, the reduction of Mn and Fe and sulphate-reduction in the lower part of the TZBA (Herrero et al. 2021a).

Under these conditions, the widely distributed populations of RF1, RF2, RF3 and RF4 in the boreholes are identified as facultative microorganisms. Propionibacterium acnes (RF1 Table 1) is mostly considered to be an anaerobic bacterium, although some strains have been identified as facultative or microaerophilic (Stackebrandt et al., 2006). Acidithiobacillus ferrooxidans, (RF2 Table 1) is a facultative aerobic organism that, in the absence of oxygen, is able to use Fe³⁺ as a final electron acceptor (Ohmura et al. 2002). RF3, identified as an aerobic bacterium of the genera Streptomyces and/or Arthrobacter, as Streptomyces sp., is capable of growing under microaerobic conditions and surviving under anaerobic conditions (Van Keulen et al. 2007), and Arthrobacter sp. can grow under anaerobic conditions using fermentation and nitrate ammonification (Eschbach et al. 2003). RF4 was identified as Streptococcus sp. and/or Aerococcus sp. Streptococcus is a facultative organism (Hardie and Whiley 2006), probably derived from agricultural fertilisers that have adapted to the environment (Zhao et al., 2012) and Aerococcus sp. is an aerobic facultative organism (Das and Kazy 2014).

In the oxygen and redox conditions detected, it is possible that biofilms were present, given the capacity of Propionibacterium sp. (Tyner and Patel 2016), Streptomyces sp. (Liermann et al. 2000), Terrabacter sp. (Piazza et al. 2019) and RF10 (DQ499314.1.1492) among others, to produce them. The formation of biofilm would allow a gradient of redox potential and oxygen, which would allow anaerobic microorganisms to have an active metabolism (Davey and O’Toole, 2000).

3.5 Anaerobic TEAP– reduction of Fe and reductive dehalogenation

The ability of microorganisms to reduce and/or oxidise Mn and Fe is another factor that determines the distribution of microbial populations. The complexity of the processes of the reduction and oxidation of Mn and Fe and the formation of new minerals has not allowed any statistical correlation to be found between any RF and the total Mn and Fe content in the sediment. However, the identification of several populations capable of reducing and/or oxidising these metals is well known. Acidithiobacillus ferrooxidans (RF2) oxidises Fe²⁺ under aerobic conditions, and under anaerobic conditions it is capable of reducing Fe³⁺ (Ohmura et al. 2002). Terrabacter sp. (RF6) is related to the ability to oxidise Mn, and to microbial communities that oxidise Fe (Piazza et al. 2019). Staphylococcus sp. (RF5) and Arthrobacter sp. (RF3) have the capacity to reduce Fe³⁺ (Paul et al. 2015).

The reductive dehalogenation of chloroethenes occurs in environments in which there are anaerobic TEAPs (Nijenhuis and Kuntze 2016). The presence or absence of reductive dehalogenation processes can be identified from an increase in metabolic rates (e.g., an increase of TCE with respect to PCE or an increase of cisDCE with respect to TCE (Puigserver et al., 2016)) and the presence of isotopically enriched PCE (Herrero et al. 2021b). The bivariant correlation of RF2 and RF10 with the process of reductive dehalogenation (Fig. 3H and I) does not imply that these populations can develop such a process. RF2 (Acidithiobacillus ferrooxidans) is related to Fe³⁺ reduction (Ohmura et al. 2002) and RF10 to an unidentified bacterium found in an anaerobic bioreactor (GU454879. 1.1495). Consequently, it is assumed that this relationship is due to the more anoxic conditions in which these populations are found.

3.6 Factors arising from the presence of contamination

Toxicity, evaluated via the sum of chloroethenes (CE) in the porewater, was evident for the 10TZBA-F1 (18.900 µmol CE/L) and 2UZ-F1 (10.500 µmol CE/L) samples, was lower in the 12BA-F1 (4.760 µmol CE/L) sample and was not detected in the other samples, where the concentration was lower than 2.500 µmol CE/L. Toxicity is one of the variables that decreases microbial diversity and the degree of development (10TZBA-F1 and 2UZ-F1, Fig. 3F and G, had lower values than the adjacent microbial communities). The same effect was detected in the most abundant populations of the site, RF1 and RF2 (Fig. 3H and I).

On the other hand, a relative increase in RF3 was detected in 10TZBA-F1 and in 7TZBA-F2, with the maximums of PCE in the TZBA, and of RF5 in 2UZ-F1, and a maximum of PCE in the UZ (Sect. 2.5 and Fig. 1). This increase is attributed to an absence of the toxicity effect in RF3 and RF5 and to inhibition by toxicity in the other populations.

The most abundant phylums in the subsoil were Proteobacteria, Actinobacteria and Firmicutes. The distribution of microbial communities in the sediment in the source zone of a chlorinated solvent contamination is highly complex. This distribution can be explained by a group of environmental variables that differ in importance depending on their location, given the high degree of geological and biogeochemical heterogeneity and the complex distribution of the contaminants. Communities develop differently depending on the characteristics of the surface to which they are attached, the biogeochemical conditions of the environment and the toxicity of the pollutants.

The percentage of fine materials, the capacity of the microorganisms to be transported in an aqueous environment, tolerance to changes in the concentration of dissolved oxygen, capacity to perform TEAP and toxicity are the factors that are identified as affecting the majority of the populations in this study.

The complexity of the structure of microbial communities in the sediment and the differences from the microbial communities in the groundwater point to the importance of studying these microbial communities when selecting of a bioremediation strategy and predicting the response of the microbial communities.

Ethical Approval: not applicable

-Consent to Participate: not applicable

-Consent to Publish: not applicable

-Authors Contributions: All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Jofre Herrero, Diana Puigserver and José Maria Carmona, except for the molecular data, that the data collection and analysis were performed by Jofre Herrero, Ivonne Nijenhuis and Kevin Kuntze. The first draft of the manuscript was written by Jofre Herrero and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Funding: This study was funded by projects CTM, 2005-07824 and CGL, 2008-02164/BTE (Spanish Ministry of Education).

-Competing Interests: the authors declare that they have no conflict of interest.

-Availability of data and materials: All data generated or analysed during this study are included in this published article [and its supplementary information files], except for the geochemical data, that could be found in Puigserver et al., 2016, and hydrochemical data, that could be found in Herrero et al., 2021. The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Adrian L, Löffler FE (eds) (2016) Organohalide-Respiring Bacteria. Springer Berlin Heidelberg, Berlin, Heidelberg
Atashgahi S, Lu Y, Zheng Y, et al (2017) Geochemical and microbial community determinants of reductive dechlorination at a site biostimulated with glycerol. Environ Microbiol 19:968–981. doi: 10.1111/1462-2920.13531
Balaban N, Yankelzon I, Adar E, et al (2019) The Spatial Distribution of the Microbial Community in a Contaminated Aquitard below an Industrial Zone. Water 11:2128. doi: 10.3390/w11102128
Bikandi J, San Millán R, Rementeria A, Garaizar J (2004) In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction. Bioinformatics 20:798–9. doi: 10.1093/bioinformatics/btg491
Bradley PM (2003) History and Ecology of Chloroethene Biodegradation: A Review. Bioremediat J 7:81–109. doi: 10.1080/713607980
Bruce KD (1997) Analysis of mer gene subclasses within bacterial communities in soils and sediments resolved by fluorescent-PCR-restriction fragment length polymorphism profiling. Appl Environ Microbiol 63:4914–4919
Chang Y-C, Ikeutsu K, Toyama T, et al (2011) Isolation and characterization of tetrachloroethylene-and cis-1, 2-dichloroethylene-dechlorinating propionibacteria. J Ind Microbiol Biotechnol 38:1667–1677
Chapman SW, Parker BL (2005) Plume persistence due to aquitard back diffusion following dense nonaqueous phase liquid source removal or isolation. Water Resour Res 41:. doi: 10.1029/2005WR004224
Cruz-Perez P, Stetzenbach LD, Alvarez AJ (1996) PCR detection of groundwater bacteria associated with colloidal transport
Das R, Kazy SK (2014) Microbial diversity, community composition and metabolic potential in hydrocarbon contaminated oily sludge: Prospects for in situ bioremediation. Environ Sci Pollut Res 21:7369–7389. doi: 10.1007/s11356-014-2640-2
Davey ME, O’toole GA (2000) Microbial biofilms: from ecology to molecular genetics. Microbiol Mol Biol Rev 64:847–867
Dojka MA, Hugenholtz P, Haack SK, Pace NR (1998) Microbial diversity in a hydrocarbon-and chlorinated-solvent-contaminated aquifer undergoing intrinsic bioremediation. Appl Environ Microbiol 64:3869–3877
Eschbach M, MÃ¶bitz H, Rompf A, Jahn D (2003) Members of the genus Arthrobacter grow anaerobically using nitrate ammonification and fermentative processes: anaerobic adaptation of aerobic bacteria abundant in soil. FEMS Microbiol Lett 223:227–230. doi: 10.1016/S0378-1097(03)00383-5
Fennell DE, Carroll AB, Gossett JM, Zinder SH (2001) Assessment of indigenous reductive dechlorinating potential at a TCE-contaminated site using microcosms, polymerase chain reaction analysis, and site data. Environ Sci Technol 35:1830–1839
Flynn SJ, Löffler FE, Tiedje JM (2000) Microbial community changes associated with a shift from reductive dechlorination of PCE to reductive dechlorination of cis-DCE and VC. Environ Sci Technol 34:1056–1061
Freeborn RA, West KA, Bhupathiraju VK, et al (2005) Phylogenetic analysis of TCE-dechlorinating consortia enriched on a variety of electron donors. Environ Sci Technol 39:8358–8368
Futamata H, Nagano Y, Watanabe K, Hiraishi A (2005) Unique kinetic properties of phenol-degrading Variovorax strains responsible for efficient trichloroethylene degradation in a chemostat enrichment culture. Appl Environ Microbiol 71 (2):904–911. doi: 10.1128/AEM.71.2.904
Green PN (1992) The Genus Methylobacterium. The Prokaryotes 2342-2349. Edited by A. Balows, H. G. Trüper, M. D
Griebler C, Avramov M (2015) Groundwater ecosystem services: a review. Freshw Sci 34:355–367. doi: 10.1086/679903
Griebler C, Lueders T (2009) Microbial biodiversity in groundwater ecosystems. Freshw. Biol.
Guilbeault M a, Parker BL, Cherry J a (2005) Mass and flux distributions from DNAPL zones in sandy aquifers. Ground Water 43:70–86. doi: 10.1111/j.1745-6584.2005.tb02287.x
Guo L, Wang G, Sheng Y, et al (2019) Groundwater microbial communities and their connection to hydrochemical environment in Golmud, Northwest China. Sci Total Environ 695:133848. doi: 10.1016/j.scitotenv.2019.133848
Hardie JM, Whiley RA (2006) The Genus Streptococcus-Oral. In: Dworkin M, Falkow S, Rosenberg E, et al. (eds) The Prokaryotes volume 4, 3rd edn. Springer-Verlag Inc., Singapore, pp 1421-1449. Edited by A. Balows, H. G. Trüper, M. D
Hartog N, Cho J, Parker BL, Annable MD (2010) Characterization of a heterogeneous DNAPL source zone in the Borden aquifer using partitioning and interfacial tracers: Residual morphologies and background sorption. J Contam Hydrol 115:79–89. doi: 10.1016/j.jconhyd.2010.04.004
Hata J, Miyata N, Kim E-S, et al (2004) Anaerobic degradation of< i> cis</i>-1, 2-dichloroethylene and vinyl chloride by< i> Clostridium</i> sp. strain DC1 isolated from landfill leachate sediment. J Biosci Bioeng 97:196–201
He J, Sung Y, Dollhopf ME, et al (2002) Acetate versus hydrogen as direct electron donors to stimulate the microbial reductive dechlorination process at chloroethene-contaminated sites. Environ Sci Technol 36:3945–3952. doi: 10.1021/es025528d
He YT, Wilson JT, Su C, Wilkin RT (2015) Review of Abiotic Degradation of Chlorinated Solvents by Reactive Iron Minerals in Aquifers. Groundw Monit Remediat 35:57–75. doi: 10.1111/gwmr.12111
Hendrickson ER, Payne JA, Young RM, et al (2002) Molecular analysis of Dehalococcoides 16S ribosomal DNA from chloroethene-contaminated sites throughout North America and Europe. Appl Environ Microbiol 68:485–495
Herrero J, Puigserver D, Nijenhuis I, et al (2019) Combined use of ISCR and biostimulation techniques in incomplete processes of reductive dehalogenation of chlorinated solvents. Sci Total Environ 648:819–829. doi: 10.1016/J.SCITOTENV.2018.08.184
Herrero J, Puigserver D, Nijenhuis I, et al (2021a) The role of ecotones in the dehalogenation of chloroethenes in alluvial fan aquifers. Environ Sci Pollut Res. doi: http://dx.doi.org/10.1007/s11356-021-12538-0
Herrero J, Puigserver D, Parker BL, Carmona JM (2021b) A new Method for determining Compound Specific Carbon Isotope of Chlorinated Solvents in Porewater. Groundw Monit Remediat. doi: 10.1111/gwmr.12435
Hohnstock-Ashe AM, Plummer SM, Yager RM, et al (2001) Further biogeochemical characterization of a trichloroethene-contaminated fractured dolomite aquifer: electron source and microbial communities involved in reductive dechlorination. Environ Sci Technol 35:4449–4456
Humphries JA, Ashe AMH, Smiley JA, Johnston CG (2005) Microbial community structure and trichloroethylene degradation in groundwater. J Microbiol 51:433–439. doi: 10.1139/W05-025
Kent AD, Smith DJ, Benson BJ, Triplett EW (2003) Web-based phylogenetic assignment tool for analysis of terminal restriction fragment length polymorphism profiles of microbial communities. Appl Environ Microbiol 69:6768–6776
Liermann LJ, Barnes AS, Kalinowski BE, et al (2000) Microenvironments of pH in biofilms grown on dissolving silicate surfaces. Chem Geol 171:1–16. doi: 10.1016/S0009-2541(00)00202-3
Liu W-T, Marsh TL, Cheng H, Forney LJ (1997) Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl Environ Microbiol 63:4516–4522
Mackay D, Shiu WY, Ma KC, Lee SC (2006) Handbook of physical-chemical properties and environmental fate for organic chemicals: Volume I: Introduction and hydrocarbons. CRC/Taylor & Francis
Mahmoudi D, Rezaei M, Ashjari J, et al (2020) Impacts of stratigraphic heterogeneity and release pathway on the transport of bacterial cells in porous media. Sci Total Environ 729:138804. doi: 10.1016/j.scitotenv.2020.138804
McMahon PB, Chapelle FH (2008) Redox processes and water quality of selected principal aquifer systems. Ground Water 46:259–71. doi: 10.1111/j.1745-6584.2007.00385.x
Men Y, VerBerkmoes NC, Shah NW, et al (2012) Sustainable syntrophic growth of Dehalococcoides ethenogenes strain 195 with Desulfovibrio vulgaris Hildenborough and Methanobacterium congolense: global transcriptomic and proteomic analyses. ISME J 6:410–421
Moreno B, Nogales R, Macci C, et al (2011) Microbial eco-physiological profiles to estimate the biological restoration of a trichloroethylene-contaminated soil. Ecol Indic 11:1563–1571
Murtagh F, Legendre P (2011) Ward’s Hierarchical Clustering Method: Clustering Criterion and Agglomerative Algorithm. arXiv:11116285
Ndon UJ, Randall AA, Khouri TZ (2000) Reductive dechlorination of tetrachloroethylene by soil sulfate-reducing microbes under various electron donor conditions. Environ Monit Assess 60:329–336. doi: 10.1023/A:1006186220604
Nijenhuis I, Kuntze K (2016) Anaerobic microbial dehalogenation of organohalides — state of the art and remediation strategies. Curr Opin Biotechnol 38:33–38. doi: 10.1016/J.COPBIO.2015.11.009
Ohmura N, Sasaki K, Matsumoto N, Saiki H (2002) Anaerobic Respiration Using Fe3+ S0, and H2 in the Chemolithoautotrophic Bacterium Acidithiobacillus ferrooxidans. J Bacteriol 184:2081–2087. doi: 10.1128/JB.184.8.2081
Pankow JF, Cherry JA (1996) Dense chlorinated solvents and other DNAPLs in groundwater: History, behavior, and remediation. Waterloo Press
Parker BL, Cherry JA, Chapman SW, Guilbeault MA (2003) Review and analysis of chlorinated solvent dense nonaqueous phase liquid distributions in five sandy aquifers. Vadose Zo J 2:116–137. doi: 10.2113/2.2.116
Paul D, Kazy SK, Banerjee T Das, et al (2015) Arsenic biotransformation and release by bacteria indigenous to arsenic contaminated groundwater. Bioresour Technol 188:14–23. doi: 10.1016/j.biortech.2015.02.039
Paul EA, Clark FE (1989) Soil Microbiology and Biochemistry, Academic. San Diego, Calif
Perović M, Obradović V, Kovačević S, et al (2017) Indicators of Groundwater Potential for Nitrate Transformation in a Reductive Environment. Water Environ Res 89:4–16. doi: 10.2175/106143016x14733681696121
Piazza A, Casalini LC, Pacini VA, et al (2019) Environmental bacteria involved in manganese(II) oxidation and removal from groundwater. Front Microbiol 10:1–13. doi: 10.3389/fmicb.2019.00119
Puigserver D, Carmona JM, Cortés A, et al (2013) Subsoil heterogeneities controlling porewater contaminant mass and microbial diversity at a site with a complex pollution history. J Contam Hydrol 144:1–19. doi: https://doi.org/10.1016/j.jconhyd.2012.10.009
Puigserver D, Herrero J, Parker BL, Carmona JM (2020) Natural attenuation of pools and plumes of carbon tetrachloride and chloroform in the transition zone to bottom aquitards and the microorganisms involved in their degradation. Sci Total Environ 712:135679. doi: 10.1016/j.scitotenv.2019.135679
Puigserver D, Herrero J, Torres M, et al (2016) Reductive dechlorination in recalcitrant sources of chloroethenes in the transition zone between aquifers and aquitards. Environ Sci Pollut Res 23:18724–18741. doi: 10.1007/s11356-016-7068-4
Rivett MO, Buss SR, Morgan P, et al (2008) Nitrate attenuation in groundwater: A review of biogeochemical controlling processes. Water Res. 42:4215–4232
Rossi P, Gillet F, Rohrbach E, et al (2009) Statistical assessment of variability of terminal restriction fragment length polymorphism analysis applied to complex microbial communities. Appl Environ Microbiol 75:7268–7270. doi: 10.1128/AEM.00135-09
Rossi P, Shani N, Kohler F, et al (2012) Ecology and Biogeography of Bacterial Communities Associated with Chloroethene-Contaminated Aquifers. Front Microbiol 3:260. doi: 10.3389/fmicb.2012.00260
Smidt H, de Vos WM (2004) Anaerobic microbial dehalogenation. AnnuRevMicrobiol 58:43–73. doi: 10.1146/annurev.micro.58.030603.123600
Stroo HF (2010) Bioremediation of chlorinated solvent plumes. In: In situ Remediation of Chlorinated Solvent Plumes. pp 309–324
Stroo HF, Leeson A, Marqusee JA, et al (2012) Chlorinated ethene source remediation: lessons learned. Environ Sci Technol 46:6438–47. doi: 10.1021/es204714w
Sung Y, Ritalahti K (2003) Characterization of two tetrachloroethene-reducing, acetate-oxidizing anaerobic bacteria and their description as Desulfuromonas michiganensis sp. nov. Appl … 69:2964–2974. doi: 10.1128/AEM.69.5.2964
Towner KJ (1992) The Genus Acinetobacter. The Prokaryotes 3137-3143. Edited by A. Balows, H. G. Trüper, M. D
Tyner H, Patel R (2016) Propionibacterium acnes biofilm - A sanctuary for Staphylococcus aureus? Anaerobe 40:63–67. doi: 10.1016/j.anaerobe.2016.05.014
Van Horn DJ, Van Horn ML, Barrett JE, et al (2013) Factors Controlling Soil Microbial Biomass and Bacterial Diversity and Community Composition in a Cold Desert Ecosystem: Role of Geographic Scale. PLoS One 8:e66103. doi: 10.1371/journal.pone.0066103
Van Keulen G, Alderson J, White J, Sawers RG (2007) The obligate aerobic actinomycete Streptomyces coelicolor A3(2) survives extended periods of anaerobic stress. Environ Microbiol 9:3143–3149. doi: 10.1111/j.1462-2920.2007.01433.x
Vela AI, García N, Latre M V, et al (2007) Aerococcus suis sp. nov., isolated from clinical specimens from swine. Int J Syst Evol Microbiol 57:1291–4. doi: 10.1099/ijs.0.64537-0
Velasco Ayuso S, Acebes P, López-Archilla AI, et al (2009) Environmental factors controlling the spatiotemporal distribution of microbial communities in a coastal, sandy aquifer system (Doñana, southwest Spain). Hydrogeol J 17:767–780. doi: 10.1007/s10040-008-0397-8
Wiedemeier T, Rifai H, Newell C, Wilson J (1999) Natural attenuation of fuels and chlorinated solvents in the subsurface. Wiley, New York
Yang Y, McCarty PL (2000) Biologically enhanced dissolution of tetrachloroethene DNAPL. Environ Sci Technol 34:2979–2984. doi: 10.1021/es991410u
Yang Y, McCarty PL (2002) Comparison between donor substrates for biologically enhanced tetrachloroethene DNAPL dissolution. Environ Sci Technol 36:3400–3404. doi: 10.1021/es011408e
Zhao WQ, Liu X, Huang QY, et al (2012) Sorption of Streptococcus suis on various soil particles from an Alfisol and effects on pathogen metabolic activity. Eur J Soil Sci 63:558–564

Supplementaryinformation.docx

Download PDF

Journal Publication

published 05 Aug, 2021

Read the published version in Environmental Science and Pollution Research →

Reviews received at journal
14 Apr, 2021
Reviewers invited by journal
13 Apr, 2021
Editor assigned by journal
05 Apr, 2021
First submitted to journal
01 Apr, 2021

You are reading this latest preprint version

Key Factors Controlling Microbial Distribution on a DNAPL Source Area

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Material And Methods

2.1. Site description, core sampling and conservation protocol

2.2. Environmental data analysis and treatment

2.3. Molecular analysis

2.4. Data treatment

2.5 Conceptual model

3. Results And Discussion

3.1 Comparison of microbial community structure

3.2 Geological factors

3.3 Hydrogeological factors

3.4 Oxygen Tolerance

3.5 Anaerobic TEAP– reduction of Fe and reductive dehalogenation

3.6 Factors arising from the presence of contamination

4. Conclusions

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1