DISSERTATION
Titel der Dissertation
Community structure and distribution of functional microbial groups within two complex environments: microorganisms associated with marine sponges and potential sulfur-compounds-reducing microorganisms in peatlands.
angestrebter akademischer Grad
Doktorin der Naturwissenschaften (Dr. rer.nat.)
Verfasserin / Verfasser: Doris Steger Matrikel-Nummer: 9803561 Dissertationsgebiet (lt. Stu- 444 Ökologie dienblatt): Betreuerin / Betreuer: Univ.-Prof. Dr. Michael Wagner
Wien, am 18 März 2010
Formular Nr.: A.04
Contents
Outline .. ....1
Part I Sulfate reducers in peatlands
Introduction 5
Microorganisms with Novel Dissimilatory (Bi)Sulfite Reductase Genes are Widespread in Geographically Distant, Low Sulfate Peatlands .21
Part II Hybridization regimes in DNA microarray experiments
Introduction .. ...79
Spatial Gradients of Target Nucleic Acids Govern DNA Microarray Hybridization Assays ..87
Part III Microorganisms in marine Sponges
Introduction 113
Sponge Associated Microorganisms: Evolution, Ecology, and Biotechnological Potential .125
Diversity and Mode of Transmission of Ammoniaoxidizing Archaea in Marine Sponges .195
Part IV Summary .. ..205
Appendix ..211
Outline
Part I of this thesis describes the analysis of dsrAB carrying microorganisms in terrestrial wetlands. Newly developed tools such as a functional gene microarray, quantitative PCR assays and dsrB based denaturing gradient gel electrophoresis as well as clone libraries were applied to unravel the widespread occurrence of deep branching dsrAB lineages, which outnumbered characterized dsrAB lineages in the Schlöppnerbrunnen fen system.
Part II validates the signal variability in DNA microarray hybridizations. A simplified microarray format and numerical simulation tests were applied to investigate systematic variations in detected signal intensities and to explore underlying mechanisms of observed spatial gradients.
Part III gives a comprehensive overview on the community structure of sponge associated microorganisms. A meta analysis of all public available sponge derived 16S rRNA gene sequences as well as of new sequences retrieved from so far unstudied marine sponges was conducted. Moreover, the mode of transmission of Archaea probably involved in ammonia oxidation was analyzed by 16S rRNA and ammonia monooxygenase gene based clone libraries. The presence of detected putative ammonia oxidizing archaea (AOA) in adult sponges and corresponding sponge larvae was confirmed by fluorescence in situ hybridization.
Part IV includes a summary of the findings obtained in this thesis
The Appendix includes a list of the author’s publications, oral and poster presentations, the acknowledgment, and the curriculum vitae.
1
Part I
Sulfate reducers in peatlands
Sulfate reducers in peatlands
Introduction
1. Peatlands and their role in the global carbon cycle
Wetlands are a common part of the global landscape, with mainly three distinguishing characteristics to other landform units: (i) the permanent or periodic presence of a water table at, near or above the land surface, which determines (ii) a unique soil development that is frequently characterized by low oxygen content and (iii) specialized plant communities that are adapted to these environments (Cowardin et al. 1995). Peatlands, also called organic wetlands, build up a substantial fraction of these systems and are probably the most widespread type of wetlands by occupying globally an area of about 4*106 km2 (Joosten et al. 2002). They represent ecosystems where net primary production exceeded the rate of decomposition and as a result predominately organic, carbon rich matter, so-called peat, accumulated over time. Their persistence depends on a constant long-term water input, and the origin of this water influences the form and function of peatlands. Based on their main water source, peatlands can be divided in predominantly groundwater-fed fens (minerotrophic peatlands) and in rainwater-fed bogs (ombrotrophic peatlands) (Bridgham et al. 1996). The source of water considerably influences the geochemistry of the peatland. Minerotrophic fens receive nutrient-input from incoming mineralized groundwater and thus higher concentrations of cations and anions are present compared to exclusively precipitation-fed bogs. Bogs are oligotrophic, as isolation from surrounding mineral soil results in low nutrient and mineral content, which is often associated with low pH values (< 5) (Bridgham et al. 1996). Dominating vegetation forms are Sphagnum mosses, sedges and ericaceous shrubs, which are adapted to these conditions (Vanbreemen 1995).
Peatlands store about 15–30% of the world’s soil carbon (Turunen et al. 2002; Davidson et al. 2006). Thus, they are massive deposits of carbon, acting as a net carbon sink and contributing to global cooling in the past 8,000 to 11,000 years (Frolking et al. 2007). Undisturbed peatlands are likely to continue functioning as net carbon sinks at the decadal scale, despite the large interannual variability of individual peatlands (Moore et al. 1998). However, environmental and anthropogenic factors, such as climate change and land-use, affect the carbon balance of these ecosystems and can lead to substantial carbon loss (Dise 2009). Primary carbon inputs in
5 Part I peatlands derive from plant assimilated carbon dioxide, whereby energy rich organic compounds are deposited as litter or released below ground by plant roots. Carbon loss is dominated by the emission of the gases carbon dioxide and methane (Keller et al. 2007).
complex polymers
exoenzymes denitrifiers - NO3 → N2 monomers manganese reducers Mn4+ → Mn2+
iron reducers
redox potential primary fermenters Fe3+ → Fe2+ Acidiphilium-, Geobacter-, Clostridiaceae, Enterobacteriaceae, and Geothrix-related spp. Actinomycetales,
sulfate reducers 2+ 2- SO4 → S Desulfomonile-, Desulforhopalus- and secondary fermenters acetogens Desulfosporosinus–related spp. Syntrophobacter-related spp. Clostridiaceae, Acidaminococcaceae methanogens
CO2 → CH4 Methanobacterium-, Methanosaeta-, Methanosarcina- and Methanospirillum-related spp.
Figure 1: Carbon mineralization processes in peatlands. Extracellular enzymes hydrolyze complex organic matter to monomers. Organic substrates, methanol or carbon dioxide are oxidized to acetate by acetogens or degraded by primary fermenters to organic acids, alcohols and amino acids. Resulting products are available either for secondary fermentation (syntrophic microorganisms) or for the microbial mediated sequential reduction of nitrate, manganese, iron (III), sulfate and finally carbon dioxide. Microorganisms identified in the Schlöppnerbrunnen fen are displayed.
Soil microorganisms that use either oxygen or other electron acceptors as energy sources to decompose organic matter to carbon dioxide, methane and dissolved organic carbon play a major role in carbon effluxes from peatlands. As oxygen supply in peatlands is commonly restricted to small volumes, mainly anaerobic microbial mineralization using alternative electron acceptors takes place (Figure 1). Crucial factors influencing these decomposition processes are on one hand the availability of electron acceptors and on the other hand nutrient content and degradability of organic matter. In this context, the activity and potential environment-mediated inhibition of extracellular enzymes play a key role for substrate supply and microbial growth, as they hydrolyze complex organic matter to low molecular weight compounds (Makoi et al. 2008). Subsequently, acetogens are capable to oxidize a range of organic substrates (e.g. glucose),
6 Sulfate reducers in peatlands methanol or carbon dioxide with hydrogen as electron donor to produce acetate. Furthermore, primary fermenters degrade sugars to fermentation products, such as formate, lactate, butyrate and propionate. Resulting end products of these transformations are then available either for secondary fermentation catalyzed by syntrophic microorganisms or for several microbially mediated anaerobic reduction processes. Generally, a sequential reduction of nitrate, manganese, iron (III), sulfate and finally carbon dioxide is assumed, whereby respective respiration processes are regulated by differences in energy yields of the reactions (Achtnich et al. 1995). However, several studies have shown that respective processes can occur spatially and temporally at the same time in these heterogeneous and highly structured systems (Wieder et al. 1990; Yavitt et al. 1990; Koretsky et al. 200 ). Processes using inorganic electron acceptors often do not fully explain total anaerobic mineralized carbon amounts in peatlands, and thus the presence of other organic electron acceptors, such as humic acids, and fermentation processes coupled to methanogenesis are supposed to contribute substantially to the decomposition of organic compounds (Lovley et al. 1996; Vile et al. 2003; Keller et al. 200 ; Wüst et al. 2009).
2. Biogeochemical cycling of sulfur compounds in peatlands
Microbial mediated transformations of sulfur compounds are closely linked to the carbon cycle as reducing processes are associated with organic matter utilization in anoxic habitats such as peatlands. Peatlands have water saturated anoxic and overlying oxic zones with fluctuating water tables, and thus dynamic sulfur cycling is expected (Figure 2).
7 Part I
2 O SO4 2 H2S DMS
2 atmospheric SO4 deposition
su rf ac e inf low O2
2 2 oxic R―SO4 SO4 H2S 2- SO4 assimilation Absorption on DOM minerals CBS 2 2 S2O3 SO4 g rou ndw ate r anoxic flow 2 2 SO4 H2S SO4
Fe3+ RIS
FeS2 ,FeS DMS
2 Figure 2: Biogeochemical sulfur cycle in peatlands. Sulfate (SO4 ) can be assimilated, 2 2 adsorbed on minerals or forms organic sulfate esters (R SO4 ). Microbial dissimilatory SO4 reduction results in hydrogen sulfide (H2S) production. In oxic zones, H2S gets reoxidized chemically to sulfate. H2S can also react to organic or carbon bonded sulfur (CBS) or inorganic reduced sulfur (RIS). RIS and CBS pools can be re-oxidized either aerobically or anaerobically. Chemical oxidation of H2S by dissolved organic matter (DOM) to thiosulfate, which subsequently can be utilized by microorganisms for thiosulfate reduction or disproportionation 2- into H2S and SO4 are proposed. Microbial mediated dimethyl sulfide (DMS) production in peatlands is suggested, which would be available for microbial reduction. Exclusively microbial driven transformations are indicated with dotted lines.
2- In bogs, sulfur is mainly introduced as sulfate (SO4 ) by atmospheric deposition, while in fens surface-water and groundwater inflows can contribute considerably to the sulfate pools. Sulfate can be assimilated by plants and microorganisms, adsorbed on soil minerals or forms organic 2- sulfate esters (R-SO4 ), which may hydrolyze under sulfate-poor conditions (Mandernack et al. 2000; Clark et al. 2005). In the anoxic zone, sulfate-reducing microorganisms (SRM) use sulfate to gain energy via dissimilatory sulfate reduction resulting in the production of sulfide. Sulfide reacts with hydrogen to hydrogen sulfide (H2S), which diffuses into the atmosphere or gets re- oxidized chemically to sulfate in oxic zones (Clark et al. 2005). H2S can also react with organic matter to form organic or carbon bonded sulfur (CBS) (Wieder et al. 1988) or with inorganic compounds, such as iron, contributing to the inorganic reduced sulfur (RIS) pool (Chapman 2001). RIS and CBS pools tend to be unstable in peatlands and can be re-oxidized either aerobically with oxygen if the water table drops (Mandernack et al. 2000), or anaerobically,
8 Sulfate reducers in peatlands using iron(III) or other compounds as anaerobic electron acceptor (Wieder et al. 1988; Blodau et al. 200 ). An alternative pathway could be chemical oxidation of H2S by dissolved organic matter to thiosulfate (Heitmann et al. 2006; Heitmann et al. 200 ), for example by quinone moieties of humic acids (Blodau et al. 200 ). Subsequently, thiosulfate can be utilized by microorganisms for thiosulfate reduction to H2S or oxidation to sulfate, respectively, or might be disproportionated into H2S and sulfate and may thus maintain the activity of present SRM (Jorgensen 1990). Finally, microbial mediated dimethyl sulfide production in peatlands is suggested (Kiene et al. 1995). Nevertheless, emissions of dimethyl sulfide from peatlands may be insignificant as it seems to be degraded rapidly by SRM and methanogens in freshwater sediments (Lomans et al. 1999).
3. The guild of sulfate-reducing microorganisms and their diversity in wetland systems
SRM represent a phylogenetically and metabolically diverse functional guild of microorganisms. Based on phylogenetic analysis of their 16S rRNA genes, at present, recognized SRM are affiliated to seven different phylogenetic phyla, two within the Archaea (i.e. the euryarchaeotal genus Archaeoglobus (Stetter 19 ) and the crenarchaeotal genera Caldivirga (Itoh et al. 1999) and Thermocladium (Itoh et al. 199 )), and five within the Bacteria. The major fraction of sulfate reducing bacteria belongs to the class Deltaproteobacteria, followed by the Firmicutes harboring the sulfate reducing genera Desulfotomaculum (Skerman et al. 19 0), Desulfosporosinus (Stackebrandt et al. 199 ), and Desulfosporomusa (Sass et al. 2004). Additionally, thermophilic SRM affiliated to Nitrospira (Thermodesulfovibrio spp. (Henry et al. 1994)), Thermodesulfobacteria (i.e. the genera Thermodesulfobacterium (Zeikus et al. 19 3) and Thermodesulfatator (Moussard et al. 2004)) and the family Thermodesulfobiaceae (Thermodesulfobium narugense (Mori et al. 2003)) were isolated from various environments.
SRM are characterized by their capability to gain energy for growth and metabolism by coupling the oxidation of organic compounds or hydrogen to the reduction of sulfate to sulfide in a process termed dissimilatory sulfate reduction (Rabus et al. 2000). The multi step process of dissimilatory sulfate reduction takes place in the cytoplasm or in association with the inner side of the cytoplasmic membrane and involves the transfer of eight electrons. Before being
9 Part I reduced, sulfate has to be activated with adenosine 5’ phosphate (ATP) by an ATP sulfurylase resulting in adenosine 5’ phosphosulfate (APS) and pyrophosphate. Subsequently, APS is converted to (bi ) sulfite by an APS reductase, thereby two electrons are transferred to the sulfur ion and adenosine 5’ monophosphate is released. The reduction of sulfite to sulfide involves the transfer of 6 electrons and is catalyzed by a dissimilatory (bi)sulfite reductase (DSR). The electron transport reactions lead to the formation of a proton motive force that drives ATP synthesis by an ATPase (Rabus et al. 2000). Many SRM are capable to grow also by using numerous other electron acceptors, such as sulfite, thiosulfate, dimethylsulfoxide, nitrate or elemental sulfur (Dalsgaard et al. 1994; Moura et al. 199 ; Sass et al. 2009), which are reduced to sulfide as well. Based on their capability to oxidize low molecular mass organic compounds either completely or incompletely, two main physiological groups are roughly distinguished: SRM involved in incomplete oxidation of organic compounds resulting in acetate as an end product, or those capable of complete oxidation resulting in carbon dioxide (Muyzer et al. 200 ). Thereby, dissimilatory sulfate reduction leads to a relatively low energy gain, as for example the free energy change (∆G) of complete oxidation of acetate or lactate is 4 or 128 kJ/mol, respectively (Rabus et al. 2000). In the absence of sulfate, several SRM can also switch to the fermentation of organic substrates, or shift to a syntrophic lifestyle, such as Syntrophobacter species (de Bok et al. 2002).
Several studies used 16S rRNA gene based approaches to investigate SRM diversity and distribution in complex systems. For instance, denaturing gradient gel electrophoresis targeting the 16S rRNA genes were used to detect selected groups of SRM in marine water columns and cyanobacterial mats (Teske et al. 1996; Teske et al. 1998), in industrial bioreactors (Dar et al. 2005) or in industrial paper mill (Maukonen et al. 2006). Furthermore, several 16S rRNA gene targeting probes for fluorescence in situ hybridization for the specific detection of different phylogenetic groups of sulfate reducing bacteria are available (Lücker et al. 200 ). As high throughput method, a DNA microarray targeting all 16S rRNA genes of recognized sulfate reducers was used to characterize the diversity of sulfate reducing communities in tooth pockets and cyanobacterial mats, as well as in peatlands (Loy et al. 2002; Loy et al. 2004). To date, the application of this so called SRM phylochip was the only study used to investigate specifically the 16S rRNA gene based composition of SRM in wetland systems.
The polyphyletic origin of SRM restricts diversity studies by comparative analysis of their 16S rRNA genes as no single pair of oligonucleotides allows specific detection of all SRM using the
10 Sulfate reducers in peatlands
16S rRNA genes as targets. As alternative, phylogenetic marker targeting functional genes encoding for enzymes involved in the metabolic pathway of sulfate reduction can be used to identify SRM. For instance, the genes encoding for the APS reductase (apsA) as well as for the DSR (dsrAB) were used in several studies to unravel SRM diversity (Hipp et al. 199 ; Friedrich 2002; Wagner et al. 2005). The presence of both enzymes is not restricted exclusively to SRM, but they share these enzymes with members of the guild of sulfur oxidizing microorganisms (SOM). In contrast to the APS reductase, where apsA of SRM and SOM exhibit a patchy affiliation in comparative phylogenetic analyses (Meyer et al. 200 ), dsrAB of SRM and SOM are clearly separated in two distinct lineages, allowing a direct assignment of new environmental gene sequences to one of these guilds (Loy et al. 2009). DsrAB genes are present in all known SRM, however, some other DSR carrying microorganisms like Sporotomaculum (Qiu et al. 2003) and Pelotomaculum (Imachi et al. 200 ), whose dsrAB genes are phylogenetically not clearly separated from those of SRM, are not able to reduce sulfate but live only in syntrophic communities. Beside sulfate, dsrAB carrying microorganisms can use several alternative inorganic and organic sulfur electron acceptors for energy conservation (Rabus et al. 2000). In the absence of a suitable electron acceptor, they are also able to use fermentation to gain energy. Accordingly, dsrAB based surveys identify not exclusively SRM, but include syntrophic microorganisms as well as sulfite, thiosulfate and/or organosulfonates reducers.
Numerous dsrAB based surveys were applied to explore SRM community composition and distribution. For instance, they are well studied in marine ecosystems (Dhillon et al. 2003; Nakagawa et al. 2006; Leloup et al. 200 ; Leloup et al. 2009), where high sulfate concentrations of up to 2 mM favor their growth and as a result dissimilatory sulfate reduction was shown to be one of the main mineralization processes in marine sediments (Jørgensen 1982). Their widespread occurrence was also demonstrated in brackish systems such as estuarine sediments (Joulian et al. 2001; Jiang et al. 2009), in freshwater habitats such as lake sediments (Karr et al. 2005; Vladar et al. 2008) and in wastewater treatment plants (Schramm et al. 1999; Dar et al. 200 ). In terrestrial ecosystems, research focuses mainly on specialized waterlogged soils like rice paddies (Liu et al. 2009), whereas only few studies investigated the SRM community composition in natural freshwater wetlands such as peatlands. For instance, sulfate reducing assemblages in the Everglades were investigated using dsrAB based clone libraries and terminal restriction fragment length polymorphism (T RFLP) analysis targeting the DSR subunit A (dsrA) (Castro et al. 2002; Castro et al. 2005). Sulfate reduction rates revealed higher
11 Part I activities of SRM in eutrophic zones, which were influenced by sulfate input from agricultural runoff, than in more pristine soils (Castro et al. 2002). There, a highly diverse dsrAB composition was detected with the majority of DSR sequences affiliated to different gram- positive Desulfotomaculum spp. and to sequences of yet uncultured microorganisms. In the eutrophic regions, most of the Desulfotomaculum associated sequences were related to microorganisms that are able to completely oxidize their organic substrates, whereas in the pristine regions all Desulfotomaculum-like sequences were exclusively related to incomplete oxidizers (Castro et al. 2002). Additionally, indications were given that sulfate concentration is not the key driver for differences between SRM assemblages, but a selection by the type and amount of electron donor was suggested (Castro et al. 2005). In a more recent approach, spatial relationships between dsrAB diversity and concentrations of metals and sulfur were studied in peatland soils in western New York uncovering a close dependency of dsrAB based community composition and redox stratification. DsrAB genes were mainly detected in deeper soil layers, where sulfur accumulation was evident. In the top peat layers, where manganese accumulation indicated oxic conditions, no dsrAB genes were detected (Martinez et al. 2007).
4. The Schlöppnerbrunnen fen as a model system
The main study site of the first part of this thesis, the minerotrophic fen Schlöppnerbrunnen II, lies in the north-eastern part of Bavaria within the forested Lehstenbach catchment (Fichtelgebirge, Germany). The total sulfur concentration in the fen soil as well as in the surrounding forest and mineral soils is in the upper range of the values reported for forest soils in Europe (~3 – 6.5 g*kg-1) (Prietzel et al. 2009). Most likely these high concentrations are caused from intensive deposition of atmospheric sulfur originating from combustion of soft coal in Eastern Europe. Although there has been a substantial reduction in air sulfur emissions since the late 1980s (Klemm et al. 1999), desorbed sulfate is still transported via groundwater flow from upland soils into the groundwater-fed wetland. Most of the introduced sulfate is retained by accumulation of inorganic and organic reduced sulfur species (RIS and CBS) in the soil (Prietzel et al. 2007), whereby estimates of RIS amounts range from <10% (Paul et al. 2006) to ~35% of total sulfur (Prietzel et al. 2009). Intermediate forms, such as sulfones, are abundant and represent 32% - 41% of the sulfur pool and oxidized sulfur species account for 5% - 23% of the total sulfur pool, respectively (Prietzel et al. 2009). Additionally, sulfate concentrations of soil
12 Sulfate reducers in peatlands solutions demonstrated unexpected stable concentrations over time, which points to a well buffered system regarding sulfate dynamics (Alewell et al. 2006).
Biochemical analysis of several studies led to the conclusion that Schlöppnerbrunnen II is dominated by alternating oxidation–reduction cycles (Alewell et al. 2006; Paul et al. 2006; Prietzel et al. 2009), possibly caused by oxygen diffusion from roots of vascular plants, which are common in the fen (Paul et al. 2006). It has been shown that the habitat is characterized by the co existence of redox processes, such as nitrate, sulfate, iron, manganese reduction and methanogenesis at a small spatial scale of few cm (Alewell et al. 2006; Paul et al. 2006; Alewell et al. 2008). The presence of high dissolved organic carbon concentrations could favor the simultaneous presence of different reduction processes, as electron donors are possibly not reaction limiting (Alewell et al. 2008). Additionally, in the presence of oxygen, organic and inorganic sulfur species will be re oxidized (Paul et al. 2006) resulting in mobile sulfate formation. Sulfate can to some extent be adsorpted to iron oxide, which is present in high concentrations in the system (Abdelouas et al. 2000; Alewell et al. 2006; Alewell et al. 2008) or is used as electron acceptor by SRM (Loy et al. 2003; Loy et al. 2004; Schmalenberger et al. 200 ). Part of the sulfate is most likely lost from the fen with the drainage water. Consequently, although it was shown that the Schlöppnerbrunnen II fen can act temporarily as a hot spot for dissimilatory sulfate reduction (Alewell et al. 2001), long-term maintenance of sulfur species within the system is not expected (Paul et al. 2006; Prietzel et al. 2009).
The microbial diversity in the Schlöppnerbrunnen fen systems was examined in several studies. The application of the so-called SRM phylochip was the only study used to investigate specifically the 16S rRNA gene based composition of SRM in the Schlöppnerbrunnen fen system, whereby only Syntrophobacter- and Desulfomonile- like bacteria could be identified. Three further studies investigated the Schlöppnerbrunnen fen system using general 16S rRNA gene targeting oligonucleotides. A two-dimensional single strand conformation polymorphism (SSCP) approach revealed significant differences between the microbial communities of oxic and anoxic zones revealing sequences mainly affiliated to Proteobacteria and Acidobacteria as dominant bacterial groups (Schmalenberger et al. 2008). The application of 16S rRNA based stable isotope probing identified, in addition to the aforementioned groups, sequences related to Actinobacteria and Clostridia as active anaerobic consumers of monosaccharides like xylose and glucose (Hamberger et al. 2008). Furthermore, uncultured Archaea related to Methanosarcinaceae, Methanobacteriaceae and Crenarchaeota were identified in the heavy
13 Part I fraction of the anoxic slurries. Additional anoxic microcosm studies retrieved further formate consuming and obligate acetoclastic methanogens, related to the families Methanomicrobiaceae and Methanosaetaceae, as potential trophic partners for moderately acid tolerant fermenters (Wüst et al. 2009). Furthermore, dsrAB based clone libraries and T RFLP analyses were engaged to investigate diversity patterns of SRM in the Schlöppnerbrunnen fen. DsrAB based analyses revealed a highly diverse community structure with many novel dsrAB sequences unrelated to known SRM (Loy et al. 2004; Schmalenberger et al. 200 ). T RFLP analysis distinguished a high diverse community pattern over a depth gradient of 0 cm to 50 cm clustering in three subgroups, namely the upper 10 cm, 15 – 25 cm and 30 – 50 cm depth zones. The environment in the upper fen soil is controlled by the rhizosphere, which permits cation exchange, O2 leakage and release of organic acids. Additionally, due to fluctuations in the water table, the upper layers are exposed to oxygenation and drying, which enhance chemical or microbial renewal of reducible electron acceptors. Thus, detected SRM have to be
O2 tolerant and – as sulfate concentrations were very low – a syntrophic lifestyle of the resident SRM was suggested (Schmalenberger et al. 200 ).
It was the aim of the first part of this thesis to contribute to a better understanding of the community structure and distribution of dsrAB containing microorganisms, which are possibly involved in the biochemical sulfur cycling in terrestrial wetlands. The description of microbial diversity is the first essential step for the exploration of the composition, interactions, and dynamics of microbial communities within their natural environment. Potent and high throughput techniques are needed to compile habitat specific microbial inventories and thus, novel tools for rapid community screening, such as a dsrB based DGGE assay and a newly developed functional gene array, were applied.
14 Sulfate reducers in peatlands
5. References
Abdelouas, A., W. Lutze, et al. (2000). "Biological reduction of uranium in groundwater and subsurface soil." Sci. Total. Environ. 250(1 3): 21 35.
Achtnich, C., F. Bak, et al. (1995). "Competition for Electron-Donors among Nitrate Reducers, Ferric Iron Reducers, Sulfate Reducers, and Methanogens in Anoxic Paddy Soil." Biology and Fertility of Soils 19(1): 65-72.
Alewell, C. and M. Novak (2001). "Spotting zones of dissimilatory sulfate reduction in a forested catchment: the 34S-35S approach." Environ. Pollut. 112(3): 369-377.
Alewell, C., S. Paul, et al. (2006). "Characterizing the Redox Status in Three Different Forested Wetlands with Geochemical Data." Environ Sci Technol 40(24): 7609-7615.
Alewell, C., S. Paul, et al. (2008). "Co-regulation of redox processes in freshwater wetlands as a function of organic matter availability?" Sci Total Environ 404(2-3): 335-42.
Blodau, C., B. Mayer, et al. (2007). "Support for an anaerobic sulfur cycle in two Canadian peatland soils." J. Geophys. Res. 112: G02004.
Bridgham, S. D., J. Pastor, et al. (1996). "Multiple limiting gradients in peatlands: A call for a new paradigm." Wetlands 16(1): 45-65.
Castro, H., S. Newman, et al. (2005). "Distribution and stability of sulfate-reducing prokaryotic and hydrogenotrophic methanogenic assemblages in nutrient-impacted regions of the Florida Everglades." Appl Environ Microbiol 71(5): 2695-704.
Castro, H., K. R. Reddy, et al. (2002). "Composition and function of sulfate-reducing prokaryotes in eutrophic and pristine areas of the Florida Everglades." Appl. Environ. Microbiol. 68(12): 6129-6137.
Chapman, S. J. (2001). "Sulphur forms in open and afforested areas of two scottish peatlands." Water Air and Soil Pollution 128(1-2): 23-39.
Clark, J. M., P. J. Chapman, et al. (2005). "Influence of drought-induced acidification on the mobility of dissolved organic carbon in peat soils." Global Change Biology 11(5): 791-809.
Cowardin, L. M. and F. C. Golet (1995). "Us-Fish-and-Wildlife-Service 1979 Wetland Classification - a Review." Vegetatio 118(1-2): 139-152.
Dalsgaard, T. and F. Bak (1994). "Nitrate Reduction in a Sulfate-Reducing Bacterium, Desulfovibrio- Desulfuricans, Isolated from Rice Paddy Soil - Sulfide Inhibition, Kinetics, and Regulation." Applied and Environmental Microbiology 60(1): 291-297.
15 Part I
Dar, S. A., J. G. Kuenen, et al. (2005). "Nested PCR-denaturing gradient gel electrophoresis approach to determine the diversity of sulfate-reducing bacteria in complex microbial communities." Applied and Environmental Microbiology 71(5): 2325-2330.
Dar, S. A., L. Yao, et al. (2007). "Analysis of diversity and activity of sulfate-reducing bacterial communities in sulfidogenic bioreactors using 16S rRNA and dsrB genes as molecular markers." Appl Environ Microbiol 73(2): 594-604.
Davidson, E. A. and I. A. Janssens (2006). "Temperature sensitivity of soil carbon decomposition and feedbacks to climate change." Nature 440(7081): 165-73. de Bok, F. A., M. L. Luijten, et al. (2002). "Biochemical evidence for formate transfer in syntrophic propionate-oxidizing cocultures of Syntrophobacter fumaroxidans and Methanospirillum hungatei." Appl Environ Microbiol 68(9): 4247-52.
Dhillon, A., A. Teske, et al. (2003). "Molecular characterization of sulfate-reducing bacteria in the Guaymas Basin." Appl. Environ. Microbiol. 69(5): 2765-2772.
Dise, N. B. (2009). "Peatland response to global change." Science 326(5954): 810-1.
Friedrich, M. W. (2002). "Phylogenetic analysis reveals multiple lateral transfers of adenosine-5'- phosphosulfate reductase genes among sulfate-reducing microorganisms." J. Bacteriol. 184(1): 278-289.
Frolking, S. and N. T. Roulet (2007). "Holocene radiative forcing impact of northern peatland carbon accumulation and methane emissions." Glob. Change Biol. 13: 1079-1088.
Hamberger, A., M. A. Horn, et al. (2008). "Anaerobic consumers of monosaccharides in a moderately acidic fen." Appl Environ Microbiol 74(10): 3112-20.
Heitmann, T. and C. Blodau (2006). "Oxidation and incorporation of hydrogen sulfide by dissolved organic matter." Chemical Geology 235(1-2): 12-20.
Heitmann, T., T. Goldhammer, et al. (2007). "Electron transfer of dissolved organic matter and its potential significance for anaerobic respiration in a northern bog." Global Change Biology 13(8): 1771- 1785.
Henry, E. A., R. Devereux, et al. (1994). "Characterization of a new thermophilic sulfate-reducing bacterium Thermodesulfovibrio yellowstonii, gen. nov. and sp. nov.: its phylogenetic relationship to Thermodesulfobacterium commune and their origins deep within the bacterial domain." Arch. Microbiol. 161(1): 62-69.
Hipp, W. M., A. S. Pott, et al. (1997). "Towards the phylogeny of APS reductases and sirohaem sulfite reductases in sulfate-reducing and sulfur-oxidizing prokaryotes." Microbiology 143(Pt 9): 2891-2902.
Imachi, H., S. Sakai, et al. (2007). "Pelotomaculum propionicicum sp. nov., an anaerobic, mesophilic, obligately syntrophic, propionate-oxidizing bacterium." Int J Syst Evol Microbiol 57(Pt 7): 1487-92.
Itoh, T., K.-I. Suzuki, et al. (1999). "Caldivirga maquilingensis gen. nov., sp. nov., a new genus of rod- shaped crenarchaeote isolated from a hot spring in the Philippines." Int. J. Syst. Bacteriol. 49: 1157-1163.
16 Sulfate reducers in peatlands
Itoh, T., K. Suzuki, et al. (1998). "Thermocladium modestius gen. nov., sp. nov., a new genus of rod shaped, extremely thermophilic crenarchaeote." Int J Syst Bacteriol 48 Pt 3: 8 9 7.
Jiang, L., Y. Zheng, et al. (2009). "Vertical distribution and diversity of sulfate reducing prokaryotes in the Pearl River estuarine sediments, Southern China." FEMS Microbiol Ecol 70(2): 93-106.
Joosten, H. and D. Clarke (2002). Wise Use of Mires and Peatlands-Background and Principles Including a Framework for Decision-Making., Int. Mire Conserv. Group and Int. Peat Soc.
Jorgensen, B. B. (1990). "A Thiosulfate Shunt in the Sulfur Cycle of Marine-Sediments." Science 249(4965): 152-154.
Jørgensen, B. B. (1982). "Mineralization of organic matter in the sea-bed - the role of sulphate reduction." Nature 296: 643–645.
Joulian, C., N. B. Ramsing, et al. (2001). "Congruent phylogenies of most common small-subunit rRNA and dissimilatory sulfite reductase gene sequences retrieved from estuarine sediments." Appl. Environ. Microbiol. 67(7): 3314-3318.
Karr, E. A., W. M. Sattley, et al. (2005). "Diversity and distribution of sulfate-reducing bacteria in permanently frozen Lake Fryxell, McMurdo Dry Valleys, Antarctica." Appl Environ Microbiol 71(10): 6353- 9.
Keller, J. K. and S. D. Bridgham (2007). "Pathways of anaerobic carbon cycling across an ombrotrophic- minerotrophic peatland gradient." Limnol. Oceanogr. 52(1): 96 - 107.
Kiene, R. P. and M. E. Hines (1995). "Microbial Formation of Dimethyl Sulfide in Anoxic Sphagnum Peat." Applied and Environmental Microbiology 61(7): 2720-2726.
Klemm, O. and H. Lange (1999). "Trends of air pollution in the Fichtelgebirge Mountains, Bavaria." Environ Sci Pollut Res Int 6(4): 193-9.
Koretsky, C. M., M. Haveman, et al. (2007). "Spatial variation of redox and trace metal geochemistry in a minerotrophic fen." Biogeochemistry 86: 33-62.
Leloup, J., H. Fossing, et al. (2009). "Sulfate-reducing bacteria in marine sediment (Aarhus Bay, Denmark): abundance and diversity related to geochemical zonation." Environmental Microbiology 11(5): 1278-1291.
Leloup, J., A. Loy, et al. (2007). "Diversity and abundance of sulfate-reducing microorganisms in the sulfate and methane zones of a marine sediment, Black Sea." Environ Microbiol 9(1): 131-42.
Liu, X. Z., L. M. Zhang, et al. (2009). "Abundance and community structure of sulfate reducing prokaryotes in a paddy soil of southern China under different fertilization regimes." Soil Biology & Biochemistry 41(4): 687-694.
Lomans, B. P., H. J. Op den Camp, et al. (1999). "Role of methanogens and other bacteria in degradation of dimethyl sulfide and methanethiol in anoxic freshwater sediments." Appl Environ Microbiol 65(5): 2116- 21.
17 Part I
Lovley, D. R., J. D. Coates, et al. (1996). "Humic substances as electron acceptors for microbial respiration." Nature 382(6590): 445-448.
Loy, A., S. Duller, et al. (2009). "Reverse dissimilatory sulfite reductase as phylogenetic marker for a subgroup of sulfur-oxidizing prokaryotes." Environ Microbiol 11(2): 289-99.
Loy, A., K. Küsel, et al. (2003). "Diversity of sulfate-reducing prokaryotes in an acidic fen: Comparison of 16S rRNA gene-based oligonucleotide microarray hybridization with dissimilatory (bi)sulfite reductase gene surveys." Appl. Environ. Microbiol. submitted.
Loy, A., K. Küsel, et al. (2004). "Microarray and functional gene analyses of sulfate-reducing prokaryotes in low-sulfate, acidic fens reveal cooccurrence of recognized genera and novel lineages." Appl Environ Microbiol 70(12): 6998-7009.
Loy, A., A. Lehner, et al. (2002). "Oligonucleotide microarray for 16S rRNA gene-based detection of all recognized lineages of sulfate-reducing prokaryotes in the environment." Appl. Environ. Microbiol. 68(10): 5064-5081.
Lücker, S., D. Steger, et al. (2007). "Improved 16S rRNA-targeted probe set for analysis of sulfate- reducing bacteria by fluorescence in situ hybridization." J Microbiol Methods 69(3): 523-8.
Makoi, J. H. J. R. and P. A. Ndakidemi (2008). "Selected soil enzymes: Examples of their potential roles in the ecosystem." African Journal of Biotechnology 7(3): 181-191.
Mandernack, K. W., L. Lynch, et al. (2000). "Sulfur cycling in wetland peat of the New Jersey Pinelands and its effect on stream water chemistry." Geochimica Et Cosmochimica Acta 64(23): 3949-3964.
Martinez, C. E., C. Yanez, et al. (2007). "Biogeochemistry of metalliferous peats: Sulfur speciation and depth distributions of dsrAB genes and Cd, Fe, Mn, S, and Zn in soil cores." Environmental Science & Technology 41(15): 5323-5329.
Maukonen, J., M. Saarela, et al. (2006). "Desulfovibrionales-related bacteria in a paper mill environment as detected with molecular techniques and culture." Journal of Industrial Microbiology & Biotechnology 33(1): 45-54.
Meyer, B. and J. Kuever (2007). "Phylogeny of the alpha and beta subunits of the dissimilatory adenosine-5'-phosphosulfate (APS) reductase from sulfate-reducing prokaryotes - origin and evolution of the dissimilatory sulfate-reduction pathway." Microbiology-Sgm 153: 2026-2044.
Moore, T. R., N. T. Roulet, et al. (1998). "Uncertainty in predicting the effect of climatic change on the carbon cycling of Canadian peatlands." Climatic Change 40(2): 229-245.
Mori, K., H. Kim, et al. (2003). "A novel lineage of sulfate-reducing microorganisms: Thermodesulfobiaceae fam. nov., Thermodesulfobium narugense, gen. nov., sp. nov., a new thermophilic isolate from a hot spring." Extremophiles 7(4): 283-90.
Moura, I., S. Bursakov, et al. (1997). "Nitrate and nitrite utilization in sulfate-reducing bacteria." Anaerobe 3: 279-290.
18 Sulfate reducers in peatlands
Moussard, H., S. L'Haridon, et al. (2004). "Thermodesulfatator indicus gen. nov., sp. nov., a novel thermophilic chemolithoautotrophic sulfate reducing bacterium isolated from the Central Indian Ridge." Int J Syst Evol Microbiol 54(Pt 1): 227-33.
Muyzer, G. and A. J. Stams (2008). "The ecology and biotechnology of sulphate-reducing bacteria." Nat Rev Microbiol 6(6): 441-54.
Nakagawa, S., F. Inagaki, et al. (2006). "Microbial community in black rust exposed to hot ridge flank crustal fluids." Applied and Environmental Microbiology 72(10): 6789-6799.
Paul, S., K. Kusel, et al. (2006). "Reduction processes in forest wetlands: Tracking down heterogeneity of source/sink functions with a combination of methods." Soil Biology & Biochemistry 38(5): 1028-1039.
Prietzel, J., J. Thieme, et al. (2007). "Sulfur K-edge XANES spectroscopy reveals differences in sulfur speciation of bulk soils, humic acid, fulvic acid, and particle size separates." Soil Biology & Biochemistry 39(4): 877-890.
Prietzel, J., J. Thieme, et al. (2009). "Sulfur speciation in well-aerated and wetland soils in a forested catchment assessed by sulfur K-edge X-ray absorption near-edge spectroscopy (XANES)." Journal of Plant Nutrition and Soil Science-Zeitschrift Fur Pflanzenernahrung Und Bodenkunde 172(3): 393-403.
Qiu, Y. L., Y. Sekiguchi, et al. (2003). "Sporotomaculum syntrophicum sp. nov., a novel anaerobic, syntrophic benzoate-degrading bacterium isolated from methanogenic sludge treating wastewater from terephthalate manufacturing." Arch Microbiol 179(4): 242-9.
Rabus, R., T. Hansen, et al. (2000). Dissimilatory sulfate- and sulfur-reducing prokaryotes. The Prokaryotes: an evolving electronic resource for the microbiological community. M. Dworkin, S. Falkow, E. Rosenberg, K.-H. Schleifer and E. Stackebrandt. New York, Springer-Verlag.
Sass, H., J. Overmann, et al. (2004). "Desulfosporomusa polytropa gen. nov., sp. nov., a novel sulfate- reducing bacterium from sediments of an oligotrophic lake." Arch Microbiol 182(2-3): 204-11.
Sass, H., S. Ramamoorthy, et al. (2009). "Desulfovibrio idahonensis sp. nov., sulfate-reducing bacteria isolated from a metal(loid)-contaminated freshwater sediment." Int J Syst Evol Microbiol 59(Pt 9): 2208- 14.
Schmalenberger, A., H. L. Drake, et al. (2007). "High unique diversity of sulfate-reducing prokaryotes characterized in a depth gradient in an acidic fen." Environ Microbiol 9(5): 1317-28.
Schmalenberger, A., C. C. Tebbe, et al. (2008). "Two-dimensional single strand conformation polymorphism (SSCP) of 16S rRNA gene fragments reveals highly dissimilar bacterial communities in an acidic fen." European Journal of Soil Biology 44(5-6): 495-500.
Schramm, A., C. M. Santegoeds, et al. (1999). "On the occurrence of anoxic microniches, denitrification, and sulfate reduction in aerated activated sludge." Appl. Environ. Microbiol. 65(9): 4189-4196.
Skerman, V. B. D., V. McGowan, et al. (1980). "Approved Lists of Bacterial Names." Int J Syst Bacteriol 30: 225-420.
19 Part I
Stackebrandt, E., C. Sproer, et al. (1997). "Phylogenetic analysis of the genus Desulfotomaculum: evidence for the misclassification of Desulfotomaculum guttoideum and description of Desulfotomaculum orientis as Desulfosporosinus orientis gen. nov., comb. nov." Int. J. Syst. Bacteriol. 47(4): 1134-1139.
Stetter, K. O. (1988). "Archaeoglobus fulgidus gen. nov., sp. nov. a new taxon of extremely thermophilic Archaebacteria." Syst Appl. Microbiol. 10: 172-173.
Teske, A., N. B. Ramsing, et al. (1998). "Sulfate-reducing bacteria and their activities in cyanobacterial mats of Solar Lake (Sinai, Egypt)." Appl. Environ. Microbiol. 64(8): 2943-2951.
Teske, A., C. Wawer, et al. (1996). "Distribution of sulfate-reducing bacteria in a stratified fjord (Mariager Fjord, Denmark) as evaluated by most-probable-number counts and denaturing gradient gel electrophoresis of PCR-amplified ribosomal DNA fragments." Appl. Environ. Microbiol. 62(4): 1405-1415.
Turunen, J., E. Tomppo, et al. (2002). "Estimating carbon accumulation rates of undrained mires in Finland–application to boreal and subarctic regions." Holocene 12: 69-80.
Vanbreemen, N. (1995). "How Sphagnum Bogs down Other Plants." Trends in Ecology & Evolution 10(7): 270-275.
Vile, M. A., S. D. Bridgham, et al. (2003). "Response of anaerobic carbon mineralization rates to sulfate amendments in a boreal peatland." Ecological Applications 13(3): 720-734.
Vladar, P., A. Rusznyak, et al. (2008). "Diversity of sulfate-reducing bacteria inhabiting the rhizosphere of Phragmites australis in Lake Velencei (Hungary) revealed by a combined cultivation-based and molecular approach." Microb Ecol 56(1): 64-75.
Wagner, M., A. Loy, et al. (2005). "Functional marker genes for identification of sulfate-reducing prokaryotes." Methods Enzymol 397: 469-89.
Wieder, R. K. and G. E. Lang (1988). "Cycling of Inorganic and Organic Sulfur in Peat from Big Run Bog, West-Virginia." Biogeochemistry 5(2): 221-242.
Wieder, R. K., J. B. Yavitt, et al. (1990). "Methane production and sulfate reduction in two Appalachian peatlands." Biogeochemistry 10: 81-104.
Wüst, P. K., M. A. Horn, et al. (2009). "Trophic links between fermenters and methanogens in a moderately acidic fen soil." Environ Microbiol.
Yavitt, J. B. and G. E. Lang (1990). "Methane production in contrasting wetland sites: Response to organic-chemical components of peat and to sulfate reduction." Geomicrobiology 8(1): 27-46.
Zeikus, J. G., M. A. Dawson, et al. (1983). "Microbial ecology of volcanic sulphidogenesis: isolation and characterization of Thermodesulfobacterium commune gen. nov. and sp. nov." J. Gen. Microbiol. 129: 1159-1169.
20 Microorganisms with Novel Dissimilatory (Bi)Sulfite Reductase Genes are Widespread in Geographically Distant, Low-Sulfate Peatlands
Doris Steger1, Cecilia Wentrup1#, Christina Braunegger1, Pinsurang Deevong1,2, Manuel Hofer1, Michael Pester1, Michael Wagner1 and Alexander Loy1* 1Department of Microbial Ecology, University of Vienna, Althanstrasse 14, A-1090 Wien, Austria 2Department of Microbiology, Kasetsart University, 50 Phaholyothin Road, Jatujak, Bangkok 10900, Thailand #Current address: Department of Molecular Ecology, Max-Planck Institute for Marine Microbiology, Celsiusstr. 1, D-28359 Bremen, Germany
Draft not corrected by all co authors
Concept by D.S. and A.L. Microarray analyses performed by D.S. and M.H., DGGE analysis and dsrAB gene library construction carried out by D.S. and C.B. C.W. constructed 16S rRNA gene libraries, P.D. and C.W. performed the qPCR assays. Phylogenetic analyses done by D.S. D.S. and A.L. wrote the manuscript.
Part I
Abstract
Schlöppnerbrunnen peatlands of the Lehstenbach catchment (Germany) house so far unidentified microorganisms with phylogenetically novel variants of the dissimilatory (bi)sulfite reductase genes dsrAB. These genes are generally characteristic for microorganisms that reduce sulfate, sulfite or some organosulfonates for energy conservation, but can also be present in anaerobic syntrophs. The abundance, community dynamics, and biogeography of dsrAB carrying microorganisms in peatlands were unknown. Soils from different depths were sampled at three time points over a six year period to analyze the diversity and distribution of dsrAB containing microorganisms at the Schlöppnerbrunnen II fen site by newly developed functional gene microarray and quantitative PCR assays. Members of novel, uncultivated dsrAB lineages (approximately representing species level groups) (i) dominated a temporally stable but spatially diverse dsrAB community and (ii) represented ‘core’ members (relative abundances of 1 2%; dsrA copy number versus the total number of bacterial and archaeal 16S rRNA genes) of the autochthonous microbial community in this fen. A bacterial and archaeal 16S rRNA gene inventory of a fen soil sample with relatively high abundance of two novel dsrAB lineages was dominated by clone sequences of the phylum Acidobacteria. Members of acidobacterial subdivisions 1, 2, and 3 were present in the Schlöppnerbrunnen fen at relatively high 16S rRNA gene abundances of 5.6 26.3%, 6 15.5%, and 3.6 13.1%, respectively, with low depth dependent variation, suggesting an important biogeochemical function in oxic and suboxic peat layers. Correlation analysis of quantitative PCR data provided no indications that members of the three acidobacterial subdivisions were carriers of the two novel dsrAB variants. Finally, denaturing gradient gel electrophoresis (DGGE) and clone library based comparison of the dsrAB diversity in soils from a wet meadow, three bogs, and five fens of various geographic locations (distance ~1 400 km), identified one Syntrophobacter-related and nine uncultivated dsrAB lineages to be widespread in different peatlands. Signatures of biogeography in dsrB DGGE data were not correlated with geographic distance, but could largely be explained by soil pH and wetland type; implying that distribution of dsrAB carrying microorganisms in wetlands on the scale of a few hundred kilometers is not limited by dispersal but determined by contemporary environmental conditions.
22 Sulfate reducers in peatlands
Introduction
Peatlands contain 15% to 30% of global soil carbon (24, 114) and represent a net carbon sink that contributed to global cooling in the past 8,000 to 11,000 years (35). While peatlands are generally rather resilient to external perturbation, long term global changes such as warming, decreased precipitation, and increased atmospheric deposition of reactive nitrogen and sulfur compounds might induce transformation of peatlands to new ecosystem types, accompanied by unforeseeable changes in the carbon balance (30). Carbon loss from peatlands is mainly mediated through anaerobic microbial decomposition of organic matter to the greenhouse gases carbon dioxide and methane (55). Ten to twenty percent of globally emitted methane derives from peatlands (15), (119, 121). Primary and secondary fermentation and methanogenesis are considered to be the main carbon degradation processes because of the absence or limited availability of alternative electron acceptors. However, other microbial processes such as denitrification and dissimilatory iron and sulfate reduction can occur together with methanogenesis in the same peat soil fraction and contribute considerably to anaerobic carbon mineralization (63). Fluctuations in environmental conditions on short and long term scales govern trophic interdependencies among the involved microorganisms. Transitions between synergistic (e.g. syntrophic interspecies transfer of hydrogen/formate) and antagonstic (e.g. competition for the same substrates) microbial interactions determine the extent of carbon flow diversion away from methanogenesis. A prime example is the suppression of microorganisms catalyzing methanogenic carbon degradation by sulfate reducing microorganisms (SRM) that are energetically favored in the competition for substrates such as acetate, alcohols, and hydrogen (36, 115, 116). While sulfate concentrations are generally low in peatlands (10 to 300 M), ongoing sulfate reduction at rates (2.5 to 180 nmol cm 3 d 1) that are comparable to rates in sulfate rich environments is fuelled by an anoxic recycling of sulfur compounds via the 'thiosulfate shunt' (10). Alternative replenishment of the sulfate pool by reoxidation of reduced sulfur species in the presence of oxygen is dependent on vegetation type and alternating periods with precipitation and drought (31, 92, 124). In addition, increasing global atmospheric sulfur pollution and acid precipitation contributes to terrestrial sulfate pools and is predicted to repress methane emissions from peatlands by up to 15% within the thirst third of this century (37). Given the significance of dissimilatory sulfate reduction, it is surprising that most information about the identity of microorganisms that catalyze this process in peatlands derives from studies of a single model fen system (Schlöppnerbrunnen) located in the forested Lehstenbach
23 Part I catchment (Bavaria, Germany). Different redox processes such as fermentation (43), methanogenesis (44), Fe(III) reduction (100), and sulfate reduction (1, 75) are ongoing on the study site (3). Atmospheric deposition of sulfur originating from combustion of soft coal in Eastern Europe until the 1990s led to the accumulation of sulfur species in the soils of this catchment. Although air pollution decreased in recent years (61), previously stored soil sulfate can desorb and is transported from upland soil into the groundwater fed fen where it drives dissimilatory sulfate reduction (Alewell 2000). DNA stable isotope probing using in situ concentrations of typical 13C labeled degradation intermediates (mixture of lactate, acetate, formate, and propionate) has shown that a low abundance Desulfosporosinus species, representing on average only 0.006% of the total bacterial and archaeal 16S rRNA genes, could be responsible for a major part of sulfate reduction in the studied fen (Pester et al. submitted). In addition, other microorganisms that are potentially involved in sulfate reduction were previously detected in this fen using 16S rRNA gene and dsrAB based diversity analyses. Few dsrAB sequences were affiliated with the SRM genera Desulfomonile and Syntrophobacter, while most other sequences could derive from new taxa, as these novel dsrAB lineages have no close cultivated relative (75, 106) Pester et al. submitted). Microorganisms that respire sulfite or sulfate anaerobically depend on the dsrAB encoded key enzyme dissimilatory (bi)sulfite reductase for energy conservation and thus these genes have been widely used as markers for PCR based molecular diversity studies of this guild (29, 58, 67, 95, 117). However, also some organisms that are phylogenetically related to SRM, but have seemingly lost the ability for sulfite/sulfate reduction, can harbor dsrAB. The dsrAB sequences of these organosulfonate reducers (8, 19, 20) or syntrophs (48, 98) can be amplified by the commonly used DSR1F DSR4R PCR primer mix (74). Consistent with the notion that dsrAB can derive from non SRM, DNA stable isotope probing of dsrAB in the model peatland did not indicate that the novel lineages are linked to lactate , acetate , formate or propionate dependent sulfate reduction (Pester et al. submitted). Besides their unknown ecophysiological function in the fen, additional questions regarding the biology of these enigmatic dsrAB containing microorganisms remain unanswered: what is their taxonomic affiliation and are they endemic only to this particular fen site or more widely distributed in different types of wetlands? Based on the repeated detection of the same or similar operational taxonomic units (OTUs) of dsrAB (75, 106) we hypothesize that some of these yet unidentified dsrAB containing microorganisms are part of the autochthonous microbial community in the Schlöppnerbrunnen model peatland. Using newly developed functional gene array (FGA) and quantitative real time PCR (qPCR) assays, we analyzed the distribution and abundance of individual dsrAB lineages
24 Sulfate reducers in peatlands in different depths of the peatland soil in samples collected at three time points over a period of six years. For linking of dsrAB and phylogenetic sequence information obtained from clone library analyses of bacterial and archaeal 16S rRNA genes, we tested if abundances of selected dsrAB OTUs correlated with the abundances of the dominant 16S rRNA gene sequence clusters in the fen. We further investigated the diversity of dsrAB in nine geographically separated wetlands, representing three different wetland types, for patterns in biogeography and endemism by using dsrB denaturing gradient gel electrophoresis (DGGE) and dsrAB clone library analyses.
Materials and Methods
Description of sites and sampling of soil.
The geographic distribution and major characteristics of the nine wetland sites analyzed in this study are summarized in Table S1 and Figure S1. The main sampling site, the acidic lowland fen Schlöppnerbrunnen II, is located in the Lehstenbach catchment in the Fichtelgebirge mountains (Bavaria, Germany) (see (2, 64, 75) for a detailed site description). Samples collected in 2001 (July 24th), were retrieved from a single soil core, which was divided in four depth sections (0 7.5, 7.5 15, 15 22.5, and 22.5 30 cm) (75). In 2004 (September 21st) and in 2007 (May 16th), soil cores (diameter ~8 cm) were collected from three random locations within the site (approximate distance between locations was 15 m) and were subsequently divided in four depth sections (0 5, 5 10, 10 15, and 15 20 cm). Two additional depth sections (20 30 and 30 40 cm) were collected from two locations in 2007. Besides Schlöppnerbrunnen II, triplicate soil cores (two depth sections ~0 20 and ~20 40 cm) were taken from random locations at eight additional wetland sites in Italy and Austria (Table S1). Immediately at the sampling site, samples for molecular analyses were homogenized, frozen on dry ice for transportation, and stored at 80 °C upon arrival in the laboratory.
DNA extraction.
For microarray, DGGE and qPCR analysis, DNA was extracted from approximately 250 mg (wet weight) of each soil sample using the Power SoilTM DNA Kit (MoBio Laboratories, Solana Beach, CA). For the preparation of a 16S rRNA gene library from Schlöppnerbrunnen II fen soil
25 Part I
(sampled in 2004, 10 15 cm soil section), genomic DNA was additionally extracted from triplicate soil samples using the protocol of Zhou et al. (127) and a slightly modified version of the protocol of Griffiths et al. (40). In contrast to the original protocol, precipitation of nucleic acids was performed with 0.1 volume of sodium acetate (pH 5.2) and 0.6 volume of isopropanol for 2 h at room temperature.
Fluorescence labeling of target genes for microarray hybridization.
An approximately 1.9 kb fragment of dsrAB was PCR amplified from 5 ng reference clone DNA or 50 ng environmental DNA using the degenerate primers DSR1Fmix (equimolar mixture of 10 M of each DSR1F variant) and DSR4Rmix (equimolar mixture of 10 M of each DSR4R variant) (Table S2). 16S rRNA and nucleotide transport protein (ntt) reference genes, targeted by control probes, were amplified by primer pairs listed in Table S3 (50 M of each primer). All forward primers contained a T3 promoter site sequence (5’ AATTAACCCTCACTAAAGGG 3’) at the 5′ end to enable T3 RNA polymerase-based reverse transcription labeling of the PCR products. PCR mixtures (50 l) were prepared by using 2 mM of each deoxynucleoside triphosphate, 10x
Taq buffer, 2 mM MgCl2 and 2 U of Taq DNA polymerase (Fermentas Inc., Hanover, MD, USA). For amplification of environmental samples, 1 l of bovine serum albumine (20 g/ l, New England BioLabs Inc., Beverly, MA, USA) was additionally added to each reaction to enhance PCR efficiency. Reference and control genes were amplified using an initial denaturation step at 95°C for 3 min, followed by 30 cycles of denaturation at 95°C for 30 s, annealing at 48°C (dsrAB), 52°C (16S rRNA gene) or 55°C (ntt) for 30 s, and elongation at 72°C for 1 min 10 sec. The cycling was completed by a final elongation step at 72°C for 3 min. PCR amplification of dsrAB from environmental DNA extracts was carried out by two successive “hot start” PCRs to minimize unspecific amplification products (23). The first PCR was performed in touch down mode by using an initial denaturation step at 95°C for 3 min, followed by 10 cycles of denaturation at 95°C for 30 s, annealing at 58 48°C for 30 s (after each cycle the temperature was reduced by 1°C), and elongation at 72°C for 1 min 10 s. Additional 10 PCR cycles at a constant annealing temperature of 48°C were performed prior to the final elongation step at 72°C for 3 min. PCR products were examined by 1% agarose gel electrophoresis for presence and sizes of amplicons. PCR products were purified from the gel using the Montage™ DNA Gel Extraction kit (Millipore, Bedford, MA, USA) according to the manufacturer’s instructions. One microliter of purified PCR product was reamplified in a second PCR applying an initial denaturation step at 95°C for 3 min, followed by 20 cycles of denaturation at 95°C for 30 s,
26 Sulfate reducers in peatlands annealing at 48°C for 30 s, and elongation at 72°C for 1 min 10 sec. The cycling was completed by a final elongation step at 72°C for 3 min. PCR products were purified using the QIAquick PCR purification kit (QIAgen, Hilden, Germany). Target labeling was achieved by in vitro transcription according to the following protocol. 500 ng of purified PCR product, 4 l 5x T3 RNA polymerase buffer (Fermentas), 2 l dithiothreitol (100 mM), 0.5 l RNAsin (40 U/ l) (Promega GmbH, Mannheim, Germany), 1 l each of ATP, CTP, and GTP (each 10 mM), 0.5 l UTP (10 mM), 2 l T3 RNA polymerase (20 U/ l) (Fermentas) and 0.75 l Cy3 UTP (5 mM) were adjusted with RNase free water to a total volume of 20 l and incubated at 37°C for 4 hours. The DNA template was subsequently degraded by adding 2 U DNase I (Fermentas) and incubation at 37°C for 15 min. Enzymatic digestion was stopped with 2 l of EDTA (25 mM, Fermentas). Following adjustment of the volume to 100 l with TE buffer (10 mM Tris HCl, pH 7.5; 1 mM EDTA), RNA was precipitated with 10 l of 5M NaCl and 300 l of ethanol. RNA was washed with 500 l of ice cold 70% ethanol and resuspended in 50 l TE buffer. RNA concentration and the amount of incorporated dye were measured with a ND 1000 spectrophotometer (Nanodrop, Inc.). Labeled RNA was fragmented by incubating with 10 mM ZnCl2 and 20 mM Tris/HCl (pH 7.4) at 60°C for 30 min. Fragmentation was stopped by the addition of 10 mM EDTA pH 8.0. Labeled RNA was divided in several aliquots that were stored at 20°C.
Probe design and manufacturing of microarrays.
Probes targeting dsrA or dsrB were designed using the probe tools of the ARB software package (78) and a dsrAB ARB database containing approximately 500 sequences (>1500 nucleotides) from pure cultures (74, 129) and environmental studies. Based on a distance matrix tree of all previously published dsrAB clone sequences retrieved from Schlöppnerbrunnen I+II system (75, 106), 146 probes with an average length of 30 bases were designed targeting different clone groups from the study site (Figure S2). Thermodynamic properties of the probes were calculated using the two state hybridization server of DINAMelt (settings: linear DNA, 37°C, 1 M Na+, 0 M Mg2+, 0.1 mM strand concentration) (80, 81). Probes lengths were adjusted to obtain similar theoretical free energy (deltaG) values. Final probes had a length of 27 to 34 nucleotides with an average deltaG of 41.0 ± 0.7 kcal/mol. The in silico specificity and number of weighted mismatches to non target dsrAB sequences of each probe were determined using probeCheck (73). Oligonucleotides for microarray spotting were synthesized by Microsynth (Balgach, Switzerland). The 5’ end of each oligonucleotide probe
27 Part I carried a T spacer consisting of 30 dTTP molecules and the 5’ terminal dTTP was aminated to allow covalent coupling of the oligonucleotides to aldehyde group coated VSS 25 glass slides (CEL Associates, Houston, Tex.). Each probe was adjusted to a concentration of 50 pmol/ l in 50% dimethyl sulfoxide and printed on the slide by using a BioRobotics MicroGrid spotter (Genomics Solutions, Ann Arbor, MI, USA) under constant temperature (20°C) and humidity (min. 50%) conditions. Each microarray contained three and six replicate spots of each dsrA/dsrB targeted and control probe, respectively. Freshly spotted DNA microarrays were incubated overnight at room temperature in a humid chamber. Slides were further treated with sodium borohydride as described previously (76).
Microarray hybridization.
For clones 80 ng and for environmental samples 500 ng of labeled dsrAB RNA fragments and defined amounts of labeled RNA targeting the control probes (Table S3) were mixed with 120 l of hybridization buffer [0.1% SDS, 0.1% 100x Denhardt’s reagent (Invitrogen), 6x SSC and 15% formamide] and incubated at 65°C for 5 15 min. The ThermoTWISTER incubator (QUANTIFOIL Instruments GmbH, Jena, Germany) and the HybriWell Sealing System HBW2222 FL (Grace BioLabs) were used for hybridization. Microarrays were hybridized for 17h at 55°C under continuous shaking at 400 rpm. Following hybridization, slides were washed by shaking at room temperature for 5 min in 2x SSC/0.1% SDS, for 5 min in 0.1x SSC, and finally for 20 s in ice cold double distilled water. Slides were dried by centrifugation (3 min, 300 g), stored at room temperature in the dark, and scanned the same day.
Scanning and data analysis.
Fluorescence images were recorded by scanning the slides at three lines to average and at 10 m resolution with a GenePix Personal 4100A array scanner (Axon Instruments, Molecular Devices Corporation, Sunnyvale, CA, USA). Gain of the photo multiplier tube was adjusted to record images with signal intensities just below the saturation level. Scanned images were saved as multilayer tiff images and analyzed with the GenePix Pro 6.0 software (Axon Instruments). Low quality hybridizations were repeated.
28 Sulfate reducers in peatlands
Microarray hybridization results were first corrected for differences in the local background according to:
−1 SBRPi = S Pi × BPi where SBRPi is the signal to noise ratio of the probe spot Pi, SPi is the median pixel intensity of the specific probe spot and BPi is the median pixel intensity of the local background area around probe Pi. To account for variations between different hybridizations, the mean signal to noise ratio x SBRCi of ten internal control oligonucleotides (dsrCONT1 to 4, dsrCONT7 to 12) and differences in labeling efficiency were used to normalize the SBRPi using the following formula:
nSBR = SBR × x SBR −1 ×[dye]−1 ×[template]×100 Pi ( Pi Ci ) where nSBRPi is the normalized signal to noise ratio of the probe spot Pi, [dye] is the concentration of incorporated Cy3 dye molecules in pmol/ l and [template] is the concentration of labeled RNA in ng/ l. Heatmaps of the microarray results were generated using the software visualization tool JColorGrid (53). Non metric multidimensional scaling ordinations based on Bray Curtis similarities between microarray results (based on the mean nSNRs of all probes in the respective probe target groups, Figure 1A) were performed using the software PRIMER 5 (http://www.primer e.com/) (18).
Quantitative real time PCR.
QPCR assays were performed using an iCycler IQ Thermocycler (Biorad Laboratories GmbH, München) and the Platinum® SYBR Green I qPCR Super Mix UDG (Invitrogen Corporation, Carlsbad, CA, USA) according to the manufacturer’s instructions. 16S rRNA gene and dsrA targeted primers were designed using the “Probe_Design” and “Probe_Match” tools of the ARB program package. Potential formation of primer dimers and theoretical melting temperatures were determined using the open source program Primer3 (http://primer3.sourceforge.net) (103) and the Oligonucleotide Properties Calculator OligoCalc (http://www.basic.northwestern.edu/biotools/oligocalc.html) (56), respectively. Primers, assay performances, and cycling conditions are given in Table 1. Thermal cycling was initiated by a denaturation step at 94°C for 3 min, followed by 40 45 cycles of denaturation at 94°C for 40 s, annealing at the specific temperature for 40 s, and elongation at 72°C for 40 45 s. PCR products were used as standards for assay calibration and were amplified from the cloned sequences given in Table 1.
29 7 6 5 6 6 6 7 6 a (target (target 10 10 10 10 10 10 10 10 2 2 2 2 2 2 2 2 Dynamic range range Dynamic genes/reaction) 10 10 10 10 10 10 10 10 R² 1.00 Linearity Linearity qPCR performance qPCR (%) Efficiency Efficiency ± 3.083.2 0.98 ± 0.681.5 0.99 ± 2.386.1 0.99 89.8 ± ± 0.489.8 0.99 ± 0.683.1 0.99 73.6 ± ± 1.473.6 0.98 ± 0.177.5 0.99 90.0 ± ± 2.8 90.0 P1K30f P1K30f P5K23f AY167468 AY167480 or clone names) clone or P4K21f, P2K11f P4K21f, X70905, X70906 X70905, X70905, X70906 X70905, AY167478, AY167465 AY167478, AY167467, AY167479, AY167479, AY167467, AY167466, AY167473, Templates for standard standard for Templates synthesis (accession n° n° (accession synthesis P4K3f, P5K18f, P5K16f, P5K18f, P4K3f, P4K23f P4K4f P7K24f et al., al., et submitted Reference Thisstudy Thisstudy Thisstudy Thisstudy Thisstudy Pester Pester 190 330 330 360 160 120 (bp) Approx. length length Approx. of PCR product PCRproduct of 750 250 250 500 [nM] 1000 conc. conc. Primer Primer 69 64 52 68.5 100062.5 360 1000 Thisstudy 100 Thisstudy temp (°C) temp Annealing Annealing and and and related related and OTU6 64 250 OTU2 OTU1 64 64 250 group 1 group 2 group 3 group
bacteria Archaea ntrophobacter ntrophobacter Target group Target Bacteria Bacteria y Acidobacteria Acidobacteria Acidobacteria Acidobacteria spp. S gene dsrA dsrA dsrA Target Target 16S rRNA 16S rRNA 16S rRNA 16S rRNA 16S 16S rRNA 16S content (%) GC(%) length Primer Primer Sequence (5' 3') (5' Sequence imers used for imers quantitativeused real time PCR. r P
Name Table 1. DsrA243FbDsrA561Rb CACGAC CACCGA TCA GCT DsrA243Fa GTA CTTGGT CAS TTC GGC C DsrA561Ra CACCACCACCGA TGA ACT DsrA216F 18 ATAGGC CTT GGC YGC BAC C 19DsrA561Re CUCGGGAAC 61 AGA CAUCGU U SYBAC836F 18 58 ATAGGC TTT GACCGC TGC 19 C GGGTAC TCA TTC SYBAC986RCTGCTG TG 68 58 19 56 CCGGGG ATG TCA AGC CCA 1389F 19 20 53 1492R 58 55 18 TG TACACA CCG CCC GTTAC GGY CTT GTTACG ACT T 67 16 16 44 50 63 Acid303FaAcid303Fb GCGCACGGM CACACT GGA Acid657R GCGCGC GGC CACACT GGAAcid702Fa ATTCCACKC ACCTCT CCC AY Acid702Fb 18 AGATAT CTG CAGCC GAA CAY Acid805R 18 67 72 AGATAT CCG CAG GAA CAT CC 19Acid306F 20 78 CTG ATSGTT Acid493RTAGGGC TAG 50 60 20 45 50 CACGGC CACACT GGC AC AGT TAGCCG CAG CTK CTT CT 50 18 20 18 50 50 55 71
30 Sulfate reducers in peatlands
The specific annealing temperature for each primer pair was determined by gradient PCR using perfectly matching target genes and a selection of non target clones with mismatches in the primer target site as templates. For each PCR product, a single band of the expected size was observed by agarose gel electrophoresis. Five nanogram of soil DNA was used as template per PCR. The specificity of each qPCR assay was confirmed by melting curve analyses of the sample derived and respective reference clone derived PCR products. Additionally, environmental PCR products were cloned and the identities of a few selected clones were verified by sequencing. Absence of PCR inhibitory substances was confirmed by qPCR analyses (using primers 1389F and 1492R) of dilution series of five selected soil DNA extracts, showing similar efficiencies and correlation coefficients as the standards.
Denaturing gradient gel electrophoresis.
Triplicate soil cores from nine different wetlands (including Schlöppnerbrunnen fen samples taken in 2004) were analyzed using newly developed forward primers for dsrB based DGGE. DNA extracts obtained from the different depth fractions of one soil core were pooled prior to PCR. Preparation of GC clamp tagged dsrB amplicons was carried out using a nested PCR approach. First, dsrAB fragments were amplified from 20 ng soil DNA in technical duplicates using the degenerate primers DSR1Fmix and DSR4Rmix (Table S2) for hot start, touchdown PCR as described above, except that only 20 cycles (10 cycles with annealing at 58 48°C, followed by 10 cycles at 48°C) were performed. The dsrAB amplicons of duplicate PCRs were mixed, purified by gel electrophoresis using the Montage™ DNA Gel Extraction kit (Millipore), and used as template for the dsrB targeted PCR. Five different PCR reactions per template were performed using the degenerate primers DSR4Rmix (10 M each variant) and the GC clamp carrying DSR1728Fmixes A to E (Table S5). Reaction mixes (total volume 100 l) were prepared with 2 mM of each deoxynucleoside triphosphate, 10x Taq buffer, 2 mM MgCl2 and 2 U of Taq DNA polymerase (Fermentas), 2 l of bovine serum albumine (20 mg/ml, New England BioLabs) and 2 l of template. Amplification was started by an initial denaturation step at 95°C for 3 min, followed by 20 cycles of denaturation at 95°C for 30 s, annealing at 55°C for 30 s, and elongation at 72°C for 20 sec. Cycling was completed by a final elongation step at 72°C for 3 min. Successfully amplified PCR products of each sample were pooled and concentrated to ~ 50 l using a vacuum centrifuge. Additionally, three dsrAB plasmids were subjected to the described procedure and used as an internal standard to allow comparison between DGGE gels. DGGE was performed as described earlier (86). In brief, denaturing
31 Part I gradient of 30 to 80% denaturants (100% denaturants mixture consist of 7 M urea and 40% formamide) was used in an 8% (w/v) polyacrylamide gel. DGGE was performed in 1x Tris acetate EDTA buffer at 60°C and at a voltage of 150 V for 6 h. Following electrophoresis, the gels were incubated for 60 min in a SYBR® Green I solution. For identification, individual bands were excised, reamplified, purified using the QIAquick PCR purification kit (Qiagen), and sequenced directly. Principal component analysis (PCA) (Canoco for Windows 4.5 (112)) was performed using a presence/absence matrix of DGGE bands and environmental variables, i.e. pH, sulfate and nitrate). Schlöppnerbrunnen II was excluded from PCA as corresponding environmental parameters were not determined for these samples. In addition, the DGGE band presence/absence matrix, environmental parameters and geographical data were used for calculation of separate Bray Curtis similarity matrices. The similarity matrices were tested pairwise for similar patterns between wetlands using the RELATE routine of the Primer 5 program (18) (999 permutation, Spearman’s rho).
Clone library construction and comparative sequence analysis.
Bacterial and archaeal 16S rRNA gene clone libraries were constructed from Schlöppnerbrunnen II fen soil (sampled in 2004, 10 15 cm soil section) DNA extracts using the primer pairs 616V 1492R (77) and Arch21F 1492R (27, 76), respectively. Prior to PCR, DNA extracts obtained from the three soil cores by using three different isolation methods (n=9) were mixed in equal parts. Cycling started with an initial denaturation step at 95°C for 5 min, followed by 25 cycles (Bacteria) or 30 cycles (Archaea) of denaturation at 95°C for 40 s, annealing at 52°C (Bacteria) or 56°C (Archaea) for 40 s, and elongation at 72°C for 1 min 30 sec. PCR was completed by a final elongation step at 72°C for 10 min. After TOPO TA cloning (Invitrogen) and sequencing, sequences were checked for chimeras with the software tools Bellerophon (46) and Pintail (http://www.bioinformatics toolkit.org/Web Pintail/) (5) and by comparing the phylogenetic placement of different parts of each 16S rRNA gene sequence (i.e. base positions 26 747 and 747 1471, numbering according to E. coli). A dsrAB clone library was construction from pooled DNA extracts of all Rasner Möser fen soil subsamples using the degenerate primers DSR1Fmix and DSR4Rmix (Table S2). Touchdown PCR was performed as described above, with the exception that 25 cycles (10 cycles with annealing at 58 48°C, followed by 15 cycles at 48°C). The cycling was completed by a final elongation step at 72°C for 10 min. The presence and sizes of the amplification products were
32 Sulfate reducers in peatlands determined by agarose (1%) gel electrophoresis. Prior to cloning using the TOPO XL cloning kit (Invitrogen), PCR products were purified from the gel using the Montage™ DNA Gel Extraction kit (Millipore). Operational taxonomic units definitions and Chao1 (17) estimates of richness were carried out using the software package DOTUR (105). Coverage was calculated according to the method described by Good (39). Phylogenetic analyses were performed by using distance matrix, maximum parsimony and maximum likelihood methods implemented in the ARB program package (78). 16S rRNA sequences were analyzed using the SILVA 16S rRNA database SSU Ref version 98 (97) and the SILVA Web Aligner SINA was applied to identify the two next related sequences of each fen soil clone. The bacterial 16S rRNA tree was inferred using PHYML and a 50% conservation filter for Bacteria (1,248 valid alignment positions). The iTOL tool was used for tree visualization (69). The archaeal 16S rRNA tree was calculated using TreePuzzle (111) and a 50% conservation filter for Archaea (1,229 valid alignment positions). Deduced DsrAB and DsrB amino acid sequences were grouped in OTUs based on an identity threshold of 90% (approximately defining species level groups, (58, 75). Phylogenetic analyses were performed using an indel filter (543 valid alignment positions after exclusion of variable regions with insertions and deletions). A consensus tree was drawn based on PHYLIP ProML (JTT), protein parsimony, and distance matrix (FITCH, JTT, global rearrangements, and randomized input order of sequences) analyses of sequences >500 amino acids. Sequences <500 amino acids were individually added to the trees without changing the overall topology by using of the ARB "parsimony interactive" method.
Sequence accession numbers. The sequences determined in this study have been deposited in the GenBank database under accession numbers GU127834 to GU127875 (archaeal 16S rRNA gene clones), GU127668 to GU127833 (bacterial 16S rRNA gene clones), GU127936 to GU127971 (dsrAB clones), GU127876 to GU127935 (dsrB clones).
33 Part I
Results
Temporal and spatial dynamics of dsrAB containing microorganisms in the Schlöppnerbrunnen II fen. A habitat specific FGA, consisting of 146 oligonucleotide probes, was developed and optimized for diversity analysis of dsrAB containing microorganisms in the Schlöppnerbrunnen fen system (Table S4, Figure S2) (75, 106). In a first step, specific hybridization conditions were established by performing a series of hybridizations with clone dsrSBI 3 (AY167467) at increasing formamide concentrations (0 to 60%, at 5% increments) in the hybridization buffer. A formamide concentration of 15% was selected for all subsequent hybridizations as the best compromise between sensitivity (i.e. high signal intensity of perfectly matched probes) and specificity (i.e. high signal intensity ratios between perfectly matched and mismatched probes). In the second step, the specificity of the dsrAB FGA was further evaluated by individual hybridization with 36 dsrAB reference clones. This set of reference clones contained at least one perfectly matched target for 120 probes. 26 probes could not be fully tested, because perfectly matching reference clones were not available. 98 of the 120 tested probes gave a positive signal only with the perfectly matching target clone. False positive hybridizations only occurred between probes and non target sequences with less than two weighted mismatches. Two probes (dsrB40 and dsrA127) were false negative and were thus excluded from further analysis. In summary, only 0.57% and 0.04% of the 5256 individual probe target hybridizations were false positive and false negative, respectively (Table S6). The nSNR values of perfectly matched probe target duplexes of the final probe set ranged from 0.04 to 2.2 (factor 55), showing that the signal intensities of individual probes varied considerably. In the third step, the sensitivity of the dsrAB FGA was assessed by a series of hybridizations containing different concentrations of the clone dsrSBI 36 (AY167469). One copy, 10, 100, 1000, 10000, and 100000 copies of the plasmid were added to 60 ng (corresponding to approximately 1x107 microbial genomes; assuming an average genome size of 5 Mbp) aliquots of fen soil DNA (sampled in 2007, depth 15 20 cm) prior to PCR amplification and further target preparation. A minimum of 1000 plasmids was required for positive detection of clone dsrSBI 36 by dsrAB FGA hybridization. Given that most dsrAB containing microorganisms only have a single copy of dsrAB in their genome, this corresponds to a detection limit of 0.01% of the total microbial community for this target
34 Sulfate reducers in peatlands organism (nSBR values of probes targeting clone dsrSBI 36 ranged from 0.17 0.41 in the specificity test). Additionally, the spot to spot variability (from triplicate spots of the same slide) and the slide to slide variability (from replicate spots of different slides) were consulted to evaluate the robustness of the microarray analysis. The spot to spot variability of five randomly chosen arrays, given as mean coefficient of variation of SBR values, was 3.8 ± 0.7% between triplicate spots on a given slide and 7.4 ± 3.7% between spots on three separate slides hybridized with the same fluorescently labeled target RNA. Technical variability of environmental microarray analyses, as determined by hybridizations with labeled RNA prepared from the same environmental sample (year 2001, depth 0 10 cm) but from separate DNA extractions and labeling reactions, showed a mean coefficient of variation of 15.3% based on nSBR values. This variability includes the composite methodological biases introduced by DNA extraction, PCR, labeling, and hybridization.
The newly developed Schlöppnerbrunnen fen specific dsrAB FGA was applied to analyze the spatial and temporal distribution of dsrAB at the site Schlöppnerbrunnen II. Microarray analyses of soil samples obtained from different depths (0 to 40 cm) were performed in technical triplicates (i.e. analysis of three separate DNA extracts obtained from the same soil core) for the year 2001 and in biological triplicates (i. e. analysis of three soil cores) for the years 2004 and 2007. Microarray procedures for target preparation, hybridization, and data analysis were kept identical for all soil samples to allow comparison between different hybridizations. Soil depth profiles of dsrAB FGA hybridization patterns that were obtained for the three years were generally similar (Figure 1A) average Bray Curtis similarity of 59±18% among all samples). However, in deeper soil layers the dsrAB diversity was considerably greater than in the top soil layer, where only a few probe target groups were detected by the FGA. In accordance with this observation, multidimensional scaling (MDS) analysis of Bray Curtis similarities between microarray hybridization patterns clearly separated the top soil layers from all other samples (Figure 1B). Members of the novel dsrAB OTUs 1 and 4 (75), which are only distantly related to Desulfobacca acetoxidans (DsrAB identities of 69 77%), were present in most deeper soil fractions, as indicated by positive probe signals for probe target groups 1 and 2 9, respectively.
35 Part I
A 8 7 9 1 2 3 4 5 6 11 14 15 17 18 20 21 24 28 29 31 34 35 40 43 46 49 52 55 56 57 10 13 16 19 22 23 25 26 27 30 32 33 36 37 38 39 41 42 44 45 47 48 50 51 53 54 58 2001 12 0 – 10 cm 10 – 17.5 cm > 0.39 17.5 – 25 cm 0.36 0.39 25 – 30 cm 0.33 0.36 0.30 0.33 2004 0.27 0.30 0 – 5 cm 0.24 0.27 5 – 10 cm 0.21 0.24 10 – 15 cm 0.18 0.21 15 – 20 cm 0.12 0.18 0.12 0.15 2007 0.09 0.12 0 – 5 cm 0.06 0.09 5 – 10 cm 0.03 0.06 10 – 15 cm < 0.03 15 – 20 cm 20 – 30 cm 30 – 40 cm OTU6 OTU3 OTU7 OTU2 OTU1 OTU8 OTU9 OTU5 OTU4 OTU11
B
Stress: 0.04 20 30 cm 30 40 cm>83% 15 20 cm >55%
10 15 cm 25 30 cm 10 17.5 cm
17.5 25 cm
5 10 cm 5 10 cm >51% 10 15 cm >71% 15 20 cm 0 5 cm 0 10 cm
0 5 cm
20012004 2007
Figure 1. Microarray based dsrAB diversity analysis of Schlöppnerbrunnen II fen soils sampled from different depths in the years 2001, 2004, and 2007. (A) Results of microarray hybridizations are displayed as mean nSBR of all probes within a probe target group and are averaged between triplicate hybridizations (results of individual replicate analyses are presented in Figure S3). The color code translates into different mean nSBR values. Probe target groups are arranged phylogenetically according to their position in schematic dsrAB neighbor joining tree. Affiliation of probe target groups to previously determined dsrAB OTUs (75) is indicated. (B) Multidimensional scaling plot based on Bray Curtis similarities between microarray hybridization patterns shown in panel A. Minimal Bray Curtis similarities of all samples in a given group of samples is indicated in color.
36 Sulfate reducers in peatlands
The following dsrAB groups with no close, cultivated relatives were also detected frequently in some of the deeper soil layers: Members of OTU 2 (probe target group 49) and the related probe target groups 45, 46, and 47, probe target groups 41, 42, and 43, members of OTU 3 (probe target groups 31 33) and OTU 11 (probe target groups 13 15). Microorganisms related to the genera Desulfomonile (probe target group 28) and Syntrophobacter (OTU 6; probe target groups 21, 23, 25) were only occasionally detected by dsrAB FGA analysis.
For confirmation of microarray results, dsrA targeted qPCR assays for OTU 1 (probe target group 1), OTU 2 (probe target group 49), and OTU 6 (probe target groups 21, 23 and 25) were developed and applied to determine the abundance of these OTUs in the Schlöppnerbrunnen II soil samples. The average total number of bacterial and archaeal 16S rRNA genes for all samples from different soil depths and years was 3.6±1.6 x 107 copies per gram of wet soil, with total numbers being slightly higher in upper soil layers (Figure 2).
Copy number ratios Relative abundance (%) OTU1 OTU2 OTU1 OTU6 OTU6 OTU2 OTU 1 OTU2 OTU 6 16.7 2.0 8.6 2001 0 – 10 cm 0.05 0.01 0.00 116.4 25.4 4.6 10 – 17.5 cm 1.05 0.23 0.01 87.0 93.8 0.9 17.5 – 25 cm 1.50 1.61 0.02 25 – 30 cm 8.0 5.7 1.4 0.41 0.29 0.05
30.5 1.9 16.3 2004 0 – 5 cm 0.09 0.01 0.00 6.9 2.5 2.7 5 – 10 cm 0.09 0.03 0.01 46.6 34.6 1.3 10 – 15 cm 1.27 0.94 0.03 5.8 14.0 0.4 15 – 20 cm 0.48 1.14 0.08 18.5 12.4 1.5 2007 0 – 5 cm 0.05 0.04 0.00 16.1 31.7 0.5 5 – 10 cm 0.09 0.17 0.01 8.0 42.0 0.2 10 – 15 cm 0.26 1.35 0.03 2.2 48.1 0.0 15 – 20 cm 0.08 1.70 0.04 2.5 26.6 0.1 20 – 30 cm 0.11 1.22 0.05 30 – 40 cm 1.8 17.4 0.1 0.07 0.69 0.04
1 10 100 1000 104 105 106 107 108 109 Gene copy number per gram of wet soil
Figure 2. Absolute and relative abundances of three dsrAB-OTUs in different depths of the Schlöppnerbrunnen II fen site in the years 2001, 2004, and 2007, as determined by qPCR. Error bars are the standard deviations of the mean for the three replicates. Black circles indicate total number of bacterial and archaeal 16S rRNA genes. Copy numbers of dsrA of OTUs 1, 2, and 6 are displayed as white, light grey, and dark grey bars, respectively. Additionally, the relative abundance of each dsrAB OTU (given in % of the total number of bacterial and archaeal 16S rRNA genes) and copy number ratios between individual dsrAB OTUs are shown for all years and depths.
37 Part I
In contrast, dsrA copy numbers of all three target OTUs were lowest in the top soil layers and generally increased with soil depth, although OTU 1 and OTU 2 were represented by slightly lower dsrA copy numbers in the deepest soil layers. In all years and soil depths, dsrA of the novel OTUs 1 and 2 outnumbered dsrA of the Syntrophobacter related OTU 6; with ratios ranging from 1.8 to 116.4 for OTU 1 and from 1.9 to 93.8 for OTU 2 (Figure 2). The mean dsrA copy number per g fen soil was 8.9 ± 0.6 x 104, for OTU 1, 2.1 ± 0.07 x 105 for OTU 2, and 7.9 ± 0.6 x 103, for the OTU 6. The mean relative abundances (percent ratio of dsrA copy numbers to total bacterial and archaeal 16S rRNA gene copy numbers averaged over soil depths) of OTUs 1, 2, and 6 were 0.75%, 0.53%, and 0.02% in the year 2001, 0.48%, 0.53%, and 0.03% in the year 2004, and 0.11%, 0.86%, and 0.03% in 2007, respectively. For Syntrophobacter related bacteria (represented by OTU 6), an additional 16S rRNA gene targeted qPCR assay was developed based on previously published 16S rRNA gene clone sequences from the Schlöppnerbrunnen fen sites (75). As expected, 16S rRNA gene and dsrA based qPCR data of the Syntrophobacter related target group showed a high correlation (Pearson correlation coefficient of 0.854, p<0.01), confirming that dsrAB of OTU 6 indeed derived from Syntrophobacter related bacteria in the Schlöppnerbrunnen fen soil samples.
Comparisons of dsrA copy numbers to the nSBR of the corresponding probe target groups of the dsrAB FGA demonstrated significant correlation (p<0.01) between qPCR and microarray data; with Pearson coefficients of 0.811 for OTU 1 (with probe target group 1), 0.756 for OTU 2 (with probe target group 49), and 0.734 (with probe target group 21) and 0.863 (with probe target group 23) for OTU 6.
Bacterial and archaeal 16S rRNA gene inventories of the Schlöppnerbrunnen II fen. Soil sampled in the year 2004 from 10 to 15 cm depth contained members of the novel OTUs 1 and 2 at a relatively high abundance and thus was chosen as a representative sample for phylogenetic identification of the abundant Schlöppnerbrunnen II fen soil microorganism. A bacterial and an archaeal 16S rRNA gene clone library were constructed and phylogenetically analyzed. After removing chimeras from the bacterial (n=9) and the archaeal (n=1) dataset, 16S rRNA gene sequences with identities ≥99% (approximately defining species level phylotypes; (109)) were grouped in OTUs/phylotypes.
38 Sulfate reducers in peatlands
Verrucomicrobia Elusimicrobia
6 Nitrospirae
1
8
4
VS
N
M
A D
M16
AB2 D A A
A Q Y2110 M
S15
M1 EU 4 162457 acidic spha 2 TM6 40 e
4 5 ent clone M
1
0 5
624 2 644 9 lon
7 9 9 8 73 X9 ac 2 P 2 r c Planctomycetes 2 ei1 6 23 1 93 tu V G 4 te id S 3 hr a a J DQ404823 contam C 2 L 22 92 nd ic avaricum EU335362 soil aggregate clone 1 AF52399 cid deg agm D w s m ne b r ra e d 10 0 Q o 66 ice p m 01CT ic na 7 0 h D m cl e 18 inutum P un s ite 1 radin l a i riu n 9 p o 32 7 e De ph g r o 62 B1 a ater pond sedim 935 a s r 17 dd n s clo d 0 G8 agn gnum peat bog clo H 5 fore u m m B AY271 (1) lta
hiz (1) y AB23877 ow g e (1) e y m pe s bacte nt FW0 P ite g (1) so e J up c e 0 (1) biu m ferriphilu dd um osp to fe n E (1) s nrichmen
lo a f m clon la e s i 325 05 ro iment clone F il iso en o (4) acC lo AF0 (1) s n s ne L p d 5162 254 a e M nd s ted wetla p il fe inated sediment clone 655 fen soil clone P7K5f face soil clone S J147 tiv h f gne e t b (3) 6 Sarobetsu Mire peat soil clone HSM SS 013 c w c soil pirillu a (2) s e n e 4 dg oil c e B HS late imic f 7 (1) fe lo oil c e BE A a re ra og d 7 s s t t D B n 23 grassla527 t b r n n s slu 6 EU043589 dune field clone GASP 0KA 526 F06 oil ioa 2K s lon ifer sed ne k eam clone s M te lu y O fe s e clo (2) D kerm lo clo K28 0 lo o t clo d EF018763 aspen rhizosphere clone Amb 1 Alte o A ed n s E il c 7 Q a T c n ce p g e P ina t il c rad il c n it n lon 5 P (1) e P3 ri 4 0 n ad ce clo (1) ne B1 o clo 8 m utu c e 5 d clo 5 f cC u o 7 so a n 0 ed aqu lo op B p lon R 1 0 (1) il ta 4 t s n e e V ne B P5 4K13f fen W clo 6 4 ro 20 ndidatu 1 n ne 1 a s sp. 829 Leptos n urfa 8 nd s s cl P one P Ba B 76318 suboxic freshw 3 55 e K10 activa coccus a e 0 l s inat 6 ia mu M fen soil clone P8K on BacB 007 ne F K 89 9 c S 038 co 5K (1) s B 7 b pe C12 56 2 5193 2 (2) f P4K 2 3 l Ca SP 5 u e ores K 23 i P oil clo o P0 9 soil 1 AJ231195 Isosphaera pallid str. e B 4 lone 0f 132 J58 38 il c f o obic 8 f F3 7 fen soil clone P7K18f 33 e DQ6 3 C AY622271 subsur s mir fen s P7 Q W4 A Q1 cin 1 f fen soil clone P2K10f 2 s er A il c K t lo 1 fen Q1 47 u 2 D n lon a 7 (1) group 3 V so K 1 EU 76 s (1) AKIW863 ne gar ne P iphila 5 D c t eSm fe 2 so P oil cl 17 f 40 oil contam 1f (1) il X71838 3 2 cl FC 7 fen soil clone P8K21f be olyticu 951 oil o DQ984569 soilfen clone soil clone IFD 124 P9K2 f DQ138961 ricen straw soil clonepaddy soils clone JG29539 one o 7 DQ A o n PT43 K 1 fe 4810 778 E AJ390478 rice roots clone PR 3 n 6 n e 1 DQ088e Sar U135 e F 335 il cl X62911 Planctom F f f o (2) P s E U 28 (1) A 17f s 286 g fen soil clone P1K3 6K C 2 E 7 (1) EU266 group 15 4 62 il clone P1K 13 fen o E ra 7 2 120 F0 ss pra f CP000473 Solibacter usitatus Ellin6076 W 3276 16S 1224 AB (1) fen s e F 7 fen soil a ment clone 35 19 n 1A 8 cy AY963300 forestirie s soil clone A fen soil clone P4K21f lo 1 31 yces lim 9f AY326547 Amazon soil clone 55 (1) c e 36 ano nd on G (1) oil clo EU335400 soil clone BacC u 028 cl S J bacte clone P2K7f DQ129640 urban aerosol clone etla il 13 fe nophilus fen soil clone P7K19f w so 11 F/G EF rial m n ne 1f P7K23f e C 494 so FFC R fen soil clone P2K11f ed c E e K il 7 st fa e lon 37 at c c DQ833485 sedi re ur f (1) lon c 5 2 lo H284 fo s 24 ile Riv lone ne AF465658 thermal8 soilu bcloneK YNPRH5il c p B 1 er P 8 s 5 so ste 045 P gra HAV 1 EU33537339 soil clone2 BacBP re096 (1) a su 050 n (1) K2 1 fen soil2 clone87 ne tu 23f g w Ns SM DQ303 itic f F5 6 lo as 5K U H la fen s Omat6H04 A 90 il c p P inin e D ne 336 Tu U nd Q o 03 ne m n o 35000 sc oil c D n s 3 m lo A ap 8 fe 83 niu t soil cl Y9 ber m B e lone 0 ra t soil c 13 rady clo P9 Q s (2) e pea 524 ag rhiz ne D n soil clo57 u re 0 fo natum clone N K2 fe 5 fo u mir 18 res ob R 3 2956 26 ssu t so ium elk.1.125 f AJ 13 bets (2) UN AY39 (1) il cl aro 26f D 543 fen soilo clon T anii AY9 6 S 6K one 2 e DU M11 36 72 P il cl DQ45 pa fen soil62 clonelone P4K4ft so AJ stu ne P8K22fNs 2 il c s 519 15 re s su3 AB so fore 652 07 fo oil c 48 en group 2 uraniu res lone f 13386 (2) m m (1) fen soilt soil clone clone P6K14f EFB1 Y9 K24f ill tailin 113 A ne P7 5f (2) g s AY913248oil clo forestP 4soilK clone DUNssu029 oil c AC6 n s one lone G 8 fe il cl 0 (1) fe uB en so acB 08 EF4 n soi H2 A f e B SP1 92927 l clon D 2 oil clon 8 WM D soil c e P4K 3 DQ984565 soil cloneK12 IFDf (2) 86 ne B0 Q2284 lone 27f 335372 se P7 oil clo 02 sub JH W EU il clon adow s glacial c HS127 n so me C65 EU3353 lone JE fe 7 tundra e FA 13 soil clo G.d8 45070 oil clon ne BacB DQ 4 forest s 16 (1) fen soil c 082 Q45150 e HSM SS 0 lone P2K19 D clone BacB 017t soil clon (1) fen s f EU335295 soil u Mire pea oil clone P4 9 Sarobets AF5239 K29f AB23877 50 forested wetla f (13 AY9 nd clone FW93 fen soil clone P2K9 t bog clone B76 13376 forest soil clo sphagnum pea ne DUNssu169 AM162428 acidic (2) fen soil clone P4K25f fen soil clone P5K27f (1) AY234551 soil isolate Ellin5134 DQ451505 forest soil clone FAC66 lev 16S 1212 AB262724 Sarobetsu m F019854 aspen rhizosphere clone E ire peat soil clone HSM050P B 13 E brifaciens fen soil clone P4K23f D86512 Acidisphaera ru (1) D2 6171 Acidobacte e RB215 fen soil clone P9K25frium capsulatum DQ386264 Siderooxidanssp lithoautotrophicushere clon ragmites rhizo e P2K2f EF5158 (8) 0323 Ph fen soil clon fen soil 77clone gras P5K16fslan AB24 ortensis 10% d soil clo (5) 11f ne FCPT602 omonas he P3K EU150198 spr echlor il clon 6 (10 AY277621 D fen so SS 9 AY96 uce fir fo clone fe 3461 rest clo (4) ctor 5K33f n soi forest soil c ne NiA biorea ne P fe l clone lo SF 9 eded il clo 6059 n so P5K ne BS e se fen so e 65 fen s il clo 18f 23 sludg lon oil c ne P5K3 (1) 5900 (1) ent c lo 1f AY94 edim 2K1f EU3 ne P (2) ted s ne P fen s353 3K1f ina I 27 oil clo67 s (6) tam il clo X AB238777 oil c con HAB ne P4Klon 4614 fen so AB262723 Sarobetsu miree Bac peat soil clone HSM050P B 1 DQ40 (1) ne C 13f 3f B 0 e clo P2K fe Sarobetsu (6)Mire peat55 soil clone HSM SS 0 rin ne n soil 8 ma cidia frumentensislo fen c 4091 il c WD260nes soil clolon AF526927 Mars OdysseyJ2 encapsulation facilityn so clone T5 3ge 17 EU335320 so e P4 A fe li 2 K lca 4 5f AY ne P1 2 (2) d soil cloneas a D1 3K fen 39 f on ne P 1 5 K30f(3) clo ne 4b s 450 p il clone BacB 043 dom oil EU oil c eu s il clo e K fen 33 lon ast (2) dy so clon 65596317f s 531 e P ure 14 AJ292673201 pollute Ps pad fen K AY oil c 9 so s CQ778955 Candidatus Endoecteinas40 ne P4 f 963 3K oi U2 ded t bog olinii72 en lone il clo 2f l c E o (1) f E s 394 lone 2 flo pea ent clolone A 1 0 U o P5 (1) 1 il c J 2 AY9633416 il forest fosoil clone AH48ne E 783 num im o S K 8 cl K1 BacB 0 B113 g ed s e 2 fen soil clone0 P3K14fon re J61 n P 46 e st soil 1f 1 A ed s fen ne f EU150219 spruce1 fir forestP clon 7 spha at clo fen soil clone P1K14ffo 4 (1) 35 5 (1) res K1f clo rtium fen soil clone P3K4 t s ne F5248 ntamin fen soil clone P4K28f (4) AS A nso soil clo oil 4 co o fen AY963454 clo 26 83 X70905g c Syntrophobacter w fen soil clone P5K3f n in DQ451448 forest soil clone FA e S rm (2) 9 DQ404 fo EU680451 forest soil clone S5 169 1 DQ451447 forest soil clone (1) 3 46 ans fen soil clone P5K6 7f fe forest soil clone BS16 1 7 n soil c tr y clone Er LLAYS 99 F E B D F C (1) ithella propionica 1K 2 f (7) T L42613 Desulforhopalus vacuolatus f fen Q 51 2 e EU clone mbI eb P 489 4 6 (1) 20 A 5 e NiA SF 31 ent clone 655973n on S fe soil clo 2 lo Acidobacteria 950 AF200962 arctic ocean isolate LSv22 l B2 150 071 03 (1) m c tepidum 55 EU234214 river clone C3 n group1 fen soil clone P1K24f 5K n J00 il CP acilis A e P A r U 1 o P 04 fen soil clo Y so 5 gr (1) AC 20 0 2 AJ831750 chemocline clone 172 oil clos g 4f X92709 peat bog75 clone21 TM232 8 1 tu as 3 (3) e F r F F f D il c ne AB K8 AF126282 Sm d n 8 2 3 l te fe Q 647 ne n sl n s S D 41 5 L spru EU542537 soil slurr K8 f fe lo dra 9K2 is A a f fe fe n 45 ne AS 9 D Q 26 P5K22 fie E fen soil clone P7K25f n 6 A 2 1 AB AX814126 D inated sedi aliangium lone 1B P J n s 91 be n il clo bac 17 a (1) m f n s A d s Q129128 forest soil clone 6 4 c lo 0 8 U 53 e P me 6 P5 M943 1513 fore dy oil clo o c 02 71 s oil 2 5 c d h 6 tob C 11 c e 2 D o sto ne 14f 0 72 e (1) oil dens 5K S 07 o 9 y ne oil clon 6 id 6 fir fo FAC 9 il 67 a i pa s P V U clo l 87 e B nze nd s on lo 7 2 ne P9K22fim K a d l c f c fen s clo MO 7K e C 24 tu ow s n c e M 16 S 1 l la t soil c c e RB 6 l c 0 ana 5 o HS 6 o ne P illu (1) 8 est l ne i 1f s n ri il c n 4 c n 1 9 AB426191 lotus field soil ss e n aro ne c rest clo o 9 r e 1 oi o DK 7 u soil 0 4 l lo n e s clo o for b o AB062751 H (1) 8 ne r st soil clo F fen soil clone P6K12f c e P5K n ore s s c b (1) v EU134562 prairie soil clone FFCH3849 n A c Ca a il 31 6 4 l P 1 b C n or dra il clo iu cter c i og niu ontam itu 7 e a c 3 gra 3 f fe ehalococcoides ethenogen 7 K e PP 90 re a t clo o m lo 89 fe n rb 49 ldilin e t DQ404843 contam 7 2 K2 13f lin 5 so s s n o t b m ne 04 1 DQ (1) 3 l al fo m m f 4 he on u e NiA U en c xib n oil clo u 03 10f (1) er p 7f m F0 E 16 45 633 ea s e s 86 (1) le fe s ina a 5 ted (2) e r J adow s (2) p d n o ke r i F dim c inin a (2) re 2 DQ A e e FA o SF E e ted
izo ae W m fen soil clone P4K7 te ina x AY9 la st s (2) p ous se
h fen soil clone P (1) u a i d g w da m M (1) r n ea g 2 046 F d se in e r t s a
o S oil t C7 ns roundwa 4 c oil c ag (1) m t p P1 on (2) n as Deltaproteobacteria CW a ite 078 s hi
o c 4 p p oi d m c te r nta lo
la lone m g im l i AB sph o clo n T SM03 H07 te ne p a c
a ent AY918928 Prolixibacter bellariivorans ile
hr bo ter clo idic F r ne Beta/Gammaproteobacteria cid P on hw s a 8 0 b c
ca
r o 6 9 P AF170103 Chloroherpeton thalassiu 6 H es lo o il c 3 0 ac WM S fr ne ne ZZ 28 8
dr es M 04 4 lo DQ404733 contaminated sediment clone 655073 ic Alphaproteobacteria 40 L 2 0 x droca ne
hy SP 62 5
V B F5 8 0P 14C JG
hy E AB2 ubo 1 98 B a
AM1 7 13
c 3 3 4 s
9
5 0a 1
2
3 1 42
216 KF
2 DQ Firmicutes 3
AY
Q676 2
D Bacteroidetes Actinobacteria Chlorobi Chloroflexi
Figure 3. 16S rRNA gene based maximum likelihood (PhyML) tree showing bacterial sequences retrieved from Schlöppnerbrunnen II fen soil. The tree includes for each OTU a representative clone (shown in bold), sequences from the next closest relatives, and representative cultivated microorganisms in the ARB SILVA database. For each OTU the number of retrieved clones is depicted in brackets. For the phylum Acidobacteria, the grouping of subdivisions 1, 2, 3, and 15 (nomenclature according (7)) is shown. Scale bar indicates 10% estimated sequence divergence.
39 Part I
The 166 bacterial and 42 archaeal clones formed 85 and 18 OTUs and covered 69% and 76% of the expected diversity, respectively. Chao1 richness analyses estimated a minimal number of 160 bacterial and 25 archaeal OTUs. The bacterial OTUs affiliated with members of the eleven described phyla Actinobacteria, Acidobacteria, Bacteroidetes, Chloroflexi, Chlorobi, Elusimicrobia, Firmicutes, Nitrospirae, Planctomycetes, Proteobacteria (including alpha , beta , gamma and delta classes), and Verrucomicrobia, and the uncultured candidate phylum TM6 (Figure 3). Most sequences (53.6%), corresponding to 33 OTUs, derived from the phylum Acidobacteria. Most OTUs (n=77) showed low similarity (<90%) to sequences of cultured members of the respective phyla. Only 8 OTUs, represented by the following clones, had 16S rRNA sequence identities of >95% (approximately defining genus level phylotypes): P2K20F (99% to Desulforhopalus vacuolatus, L42613, Deltaproteobacteria), P9K23f (97.5% to Bradyrhizobium elkanii, U35000, Alphaproteobacteria), P5K6f (96.1% to Smithella propionica, AF126282, Deltaproteobacteria), P2K2f (95.78% to Siderooxidans lithoautotrophicus, DQ386264, Betaproteobacteria), P6K11f (95.4% to Desulfosporosinus orientis, Y11570, Firmicutes), P1K21f (95.2% to Solibacter usitatus Ellin6076, AY673303, Acidobacteria), P4K25f (95.1% to isolate Ellin5134, AY234551, Alphaproteobacteria), and P3K5f (95% to Syntrophobacter wolinii, X90605/06, Deltaproteobacteria). The next relatives to the recovered 16S rRNA sequences of this study originated from a broad range of habitats. However, the majority of sequences (68%, n=87, corresponding to 64 OTUs) were closest related to sequences from soil, with several clones (n=29, corresponding to 20 OTUs) being related to sequences obtained from wetlands or water saturated soil habitats e.g. forested wetland impacted by reject coal (14), rice paddy field (89), peat bogs (26, 102), alpine tundra and wet meadows (21). Furthermore, 46 sequences (28%) had 16S rRNA sequence identities of >95% to sequences retrieved from heavy DNA fractions after stable isotope probing of Schlöppnerbrunnen II fen soil mesocosms that were incubated with and without sulfate (ZITAT MICHA). Most sequences were in common to sequences obtained from both incubations with and without sulfate and affiliated to Acidobacteria (n=36, 95.1% to 99.7%), Planctomycetes (n=1, 96.4%), Actinobacteria (n=1, 95.1 to 95.5%) and Alphaproteobacteria (95.4%). However, similar Nitrospirae sequences (n=6, 97 to 99.4%) and Firmicutes sequences (n=1, 96.6 to 97.6) were detected only in incubations with sulfate. With the exception of acidobacterial clones P5K31f (98.2% to AM270582) and P1K30f (99% to AM773947) and the Chlorobi clone P5K17f (96% to AM773938), most sequences from our study represent genus level phylotypes that have not been detected by three other 16S rRNA based studies of the Schlöppnerbrunnen fen system (43, 75, 125).
40 Sulfate reducers in peatlands
72 DQ278152, soil clone RotF 135iia 80 DQ278155, soil clone OdenF 150iiia 100 DQ278153, soil clone OdenF 150iia Schlöppnerbrunnen clone A_P1K6f (2) 76 100 77 Schlöppnerbrunnen clone A_P2K26f (1) Schlöppnerbrunnen clone A_P2K16f (3) 100 Schlöppnerbrunnen clone A_P2K32f (10) 99 97 AB243803, rice paddy field clone NRP L AB243793, rice paddy field clone NRP B Schlöppnerbrunnen clone A_P1K8f (2) 88 100 100 DQ901270, acid copper mine clone K1G 43 AB262710, Sarobetsu mire clone HSM050P A 9 DQ223194, subsurface water clone EV221H2111601H196 95 90 DQ223192, subsurface water clone EV221H2111601H168 94 100 Schlöppnerbrunnen clone A_P2K3f (2) Schlöppnerbrunnen clone A_P2K18f (1) AB262708, Sarobetsu mire clone HSM050P A 7 AB243795, rice paddy field clone NRP D Schlöppnerbrunnen clone A_P2K15f (1) 91 EU239960, Candidatus Nitrosocaldus yellowstonii EU281336, Candidatus Nitrososphaera gargensis 85 AB294263, groundwater clone SWA07 Schlöppnerbrunnen clone A_P2K11f (9) AB050208, gold mine clone SAGMA D EU309868, rhizosphere clone LAR_Cren_60 AJ576209, leachate landfill clone GZK8 97 EU280219, salt marsh sediment clone HQSAT_10F8 96 Schlöppnerbrunnen clone A_P2K21f (1) 98 Schlöppnerbrunnen clone A_P2K14f (2) Schlöppnerbrunnen clone A_P2K30f (1) Schlöppnerbrunnen clone A_P3K4f (1) Schlöppnerbrunnen clone A_P2K20f (1) U59986, sediment clone pGrfC26 AJ548934, oligotrophic fen cloneFenO 16S EU155995, rich minerotrophic fen clone MH1492_3B 76 EU155994, rich minerotrophic fen clone MH1492_12C Schlöppnerbrunnen clone A_P3K3f (1) EU155991, rich minerotrophic fen clone MH1492_4A EU155993, rich minerotrophic fen clone MH1492_B9G 100 DQ060322, Ignisphaera sp. Tok10A.S1 Crenarchaeota X14835, Thermofilum pendens 100 EU155939, rich minerotrophic fen clone MHLsu47_B8A 100 Schlöppnerbrunnen clone A_P2K23f (1) 100 EU155938, rich minerotrophic fen clone MHLsu47_B11B 100 CP001338, Candidatus Methanosphaerula palustris E19c CP000254, Methanospirillum hungatei JF1 98 100 EU155965, rich minerotrophic fen clone MH1492_B8C 100 Schlöppnerbrunnen clone A_P1K22f (1) 99 EU155964, rich minerotrophic fen clone MH1492_B6H 100 X16932, Methanosaeta concilii 100 Schlöppnerbrunnen clone A_P2K10f (2) EU155955, rich minerotrophic fen clone MH1492_B9E EU155912, rich minerotrophic fen clone MHLsu47_22D 10% Euryarchaeota
Figure 4. 16S rRNA gene based maximum likelihood (TreePuzzle) tree showing archaeal sequences retrieved from Schlöppnerbrunnen II fen soil. The tree includes for each OTU a representative clone (shown in bold), sequences from the next closest relatives, and representative cultivated microorganisms in the ARB SILVA database. For each OTU the number of retrieved clones is depicted in brackets. TreePuzzle support values over 70 are given at the respective nodes. Scale bar indicates 10% estimated sequence divergence.
Most archaeal OTUs (n=15) were affiliated with the phylum Crenarchaeota, while only 3 OTUs branched within the Euryarchaeota (Figure 4). While none of the crenarchaeal OTUs from the Schlöppnerbrunnen II fen was closely related to a cultivated microorganism (16S rRNA sequence identities <80%), most OTUs grouped with other soil and wetland derived sequences. Clone A_P2K11f, representing an OTU consisting of nine sequences, was affiliated
41 Part I with Nitrososphaera gargensis, an archaeum recently proposed to belong to a novel archaeal phylum, the Thaumarchaeota (13) (Spang et al. submitted). The euryarchaeal OTUs were closely related (sequence identities of 98 99%) to clones from a rich minerotrophic, pH neutral fen located in central New York State (16). The closest described relatives of clones A_P2K23f and A_P2K10f were Candidatus Methanosphaerula palustris and Methanosaeta concilii (each with a sequence similarity of 94%), respectively. In contrast to the bacterial sequences, a considerable proportion of archaeal sequences (43%) had highest similarities to sequences from a previously published 16S rRNA gene inventory of the Schlöppnerbrunnen fen system (125).
Abundance of Acidobacteria in the Schlöppnerbrunnen II fen. Sequences of the phylum Acidobacteria predominated in the bacterial 16S rRNA gene library. Development and application of specific 16S rRNA gene targeted qPCR assays for the three main acidobacterial subdivisions/groups 1, 2, and 3 (Fig. 3, Table 1) confirmed that each of the three groups constituted a major part of the Schlöppnerbrunnen II fen soil microbial community in all soil depths and years (Figure 5).
Total bacterial and archaeal 16S rRNA gene copy number per gram of wet soil 1 10 100 1000 104 105 106 107 108 2004 0 5 cm 5 10 cm 10 15 cm
15 20 cm 2007 0 5 cm
5 10 cm 10 15 cm 15 20 cm 20 30 cm
30 40 cm
0% 10% 20% 30% 40% Relative 16S rRNA gene abundance of acidobacterial groups Figure 5. Relative abundance of acidobacterial subdivisions/groups in different depths of the Schlöppnerbrunnen II fen site in the years 2004 and 2007, as determined by of 16S rRNA gene targeted qPCR. Black circles indicate total numbers of bacterial and archaeal 16S rRNA genes. Relative abundances (given in % of the total number of bacterial and archaeal 16S rRNA genes) of acidobacterial groups 1, 2, and 3 are indicated in black, grey, and white bars, respectively. Error bars are the standard deviations of the mean for the three biological replicates.
42 Sulfate reducers in peatlands
Members of group 1 were most abundant, followed by members of group 2 and 3: Mean copy numbers per gram wet soil were 4.3 ± 1.3 x 106, 3.9 ± 3.6 x 106, and 2.4 ± 0.4 x 106 in the year 2004 and 9 ± 2.2 x 106, 5.6 ± 5 x 106, and 4.9 ± 0.9 x 106 in 2007 for group 1, group 2, and group 3, respectively. The relative abundance was 10.3% ± 5.1, 8.5% ± 1.9, and 5.0% ± 1 in the year 2004 and 20.0% ± 3.5, 12.1% ± 6.1, and 3.5% ± 1.7 in 2007 for the groups 1, 2, and 3, respectively. In contrast to Syntrophobacter related microorganisms (see above), 16S rRNA gene abundance of any of the three acidobacterial groups was not correlated with the dsrA abundance of any of the novel OTUs 1 and 2 (p=0.123 0.875, Pearson coefficients 0.058 0.521).
Diversity and biogeography of dsrAB in geographically separated wetlands. In order to evaluate if the dsrAB containing microorganisms in the Schlöppnerbrunnen fens are endemic or show a wider geographic distribution, a modified DGGE assay was applied for fingerprinting of dsrB in nine wetlands located in Austria, Italy and Germany (Table S1, Figure S1). In comparison to the original DGGE forward primer DSRp2060F (38), which covers only 29% of all sequences with the primer binding site in the dsrAB database, the newly developed DSR1728F primer mix has a significantly improved in silico coverage of 95% (Table S5). This modified dsrB DGGE analysis yielded 11 to 40 bands for the different wetland soil samples. The lowest and highest numbers of DGGE bands were detected in the Berndorf wet meadow soil and the Roßbrand II fen soil, respectively. With exception of the Roßbrand peats I and II, located in close distance (~1 km) to each other, the dsrB DGGE profiles between replicate soil cores were more similar to each other than dsrB DGGE profiles between different wetland sites (Figure 6). Furthermore, all precipitation fed bogs showed greater similarity in dsrB DGGE banding pattern to each other than to other wetland sites. PCA resulted in four principal components (PC), which explained 74.1% of the DGGE based community composition data, whereby the first two ordination axes (Figure 6) explained 48.8% of the variation. Correlation analysis between principal components and environmental factors revealed that PC1 was strongly negatively correlated with pH (correlation coefficient of 0.9963), whereas PC2 was moderately negatively correlated with sulfate (correlation coefficient of 0.7286) and nitrate (correlation coefficient of 0.6834). Furthermore, the different peatland sites showed no major relationship based on their geographical position. In contrast, the measured environmental factors (pH, concentrations of sulfate and nitrate) explained significantly the differences in the dsrB DGGE profiles of the investigated wetlands (rho=0.568, p<0.001).
43 Part I
1.0 ABC Schallhof (F)
Rasner Möser (F) 0.5 ABC
Krähmoos (B)
PC 2 PC pH AC B 0.0 Roßbrand II (F) Schremser Hochmoor (B) BGroße A A Heide (B) AB C BC C A B C Roßbrand II (F) Roßbrand I (F) 0.5 ABC Berndorf (WM) Nitrate Sulfate