Diversity and Abundance of Sulfur-Oxidizing Bacteria in

Diversity and Abundance of Sulfur-Oxidizing Bacteria in Wadden
Sea Sediments
Revealed by DsrAB Phylogeny and dsrAB-Targeted
Real-Time PCR
Master Thesis
to obtain the academic degree Master of Science (M.Sc.)
of the International Max Planck Research School for Marine Microbiology
and the University of Bremen
Submitted by
SABINE LENK
First Reviewer: Prof. Dr. Rudolf Amann
Second Reviewer: Dr. Marc Mußmann
Submitted at: 17.03.2006
Day of defence: 31.03.2006
1
The paper at hand has been produced as a Master Thesis at the MPI for Marine
Microbiology, Department of Molecular Ecology of Prof. Dr. R. Amann during the
period of September 2005 to March 2006.
Bremen, 17.03.2006
2
Table of Contents
Table of Contents
(A) Abbreviations ............................................................................................................... 3
(B) Summary....................................................................................................................... 4
1. Introduction..................................................................................................................... 5
2. Material and Methods ................................................................................................... 14
2.1 Material ................................................................................................................... 14
2.2 Chemicals................................................................................................................ 15
2.3 Sampling ................................................................................................................. 17
2.4 Primer design .......................................................................................................... 18
2.5 Amplification of dsrAB from Wadden Sea sediments............................................ 19
2.5.1 Establishment of dsrAB clone libraries ........................................................... 19
2.5.2 Sequencing....................................................................................................... 26
2.6 Amplification of dsrAB from cultivated SOP......................................................... 28
2.7 Specificity test of dsrAB primers dsrA 240 F and dsrA 240 F* ........................... 30
2.8 Determination of dsrB copy number in Wadden Sea sediments ........................ 30
2.8.1 Primer specificity check................................................................................... 30
2.8.2 Standard Generation......................................................................................... 31
2.8.3 Quantitative Real-time PCR ............................................................................ 32
3. Results........................................................................................................................... 34
3.1 Primer optimization ................................................................................................ 34
3.2 Amplification of dsrAB from Wadden Sea sediment ............................................. 36
3.3 Phylogeny and diversity of the DsrAB from SOP .................................................. 37
3.3.1 Overall DsrAB diversity .................................................................................. 37
3.3.2 Comparative analysis of DsrAB consensus tree and 16s rRNA gene tree ...... 41
3.3.3 Comparative analysis of insertions and deletion patterns................................ 43
3.3.4 Vertical DsrAB diversity ................................................................................. 44
3.3.5 Rarefaction analysis ......................................................................................... 44
3.3.6 Comparison 16S rRNA gene identity and DsrAB identity.............................. 45
3.4 Amplification of dsrAB from cultivated SOP......................................................... 46
3.5 Primer selectivity test.............................................................................................. 46
3.6 Real-time Quantitative PCR ................................................................................... 47
4. Discussion ..................................................................................................................... 50
4.1 Primer evaluation and methodological aspects....................................................... 50
4.1.1 Evaluation of dsrAB primers for diversity studies and qRT-PCR ................. 50
4.1.2 Amplification of DsrAB gene from cultivated SOP ........................................ 51
4.2 DsrAB and 16S rRNA phylogeny of cultured SOP and evidence for LGT ........... 52
4.3 Diversity of SOP in WS sediment .......................................................................... 54
4.3.1 Overall SOP diversity ...................................................................................... 54
4.3.2 Vertical SOP diversity ..................................................................................... 56
4.4 Quantitative Real-Time PCR .................................................................................. 57
4.5 General value and limitations of DsrAB gene as molecular marker for SOP
diversity studies and qRT-PCR..................................................................................... 58
4.6 Outlook ................................................................................................................... 59
(C) Literature .................................................................................................................... 61
(D) Appendix .................................................................................................................... 66
(E) Acknowledgement ...................................................................................................... 70
(F) Statement..................................................................................................................... 71
Abbreviations
(A) Abbreviations
ARB
APS
ATP
BSA
CARD
DNA
dNTP
Dsr
dsrAB
DsrAB
et al.
F
FISH
IPTG
JSS
LB
PBS
PCR
MQ
NAD(P)
r
R
RNA
RT
SOP
S
SDS
SOP
Sox
DQR
sp.
Table
TAE
Taq
U
UK
USA
UV
V
X-Gal
Arbor
adenosine phophosulfate
adenosine triphosphate
bovine serum albuminw
Catalyzed Reporter Deposition
deoxyribonucleic acid
Desoxyribonucleosidtriphosphat
dissimilatory (bi)sulfite reductase
nucleotide sequence (gene) of the Dsr
amino acid sequence of the Dsr
et alii (und andere)
forward
Fluoreszenz in situ Hybridisierung
Isopropyl-ß-D-thiogalactosid
Janssand
Luria-Bertani
phosphate buffered saline
Polymerase Chain Reaction
MilliQ demineralized water in Millipore Quality
nicotineamid adenosine dinucleotide phosphate
ribosomal
reverse
ribonucleic acid
real-time
sulfur oxidizing prokaryotes
Svedberg
sodiumdodecyl sulfate
sulfur-oxidizing prokaryotes
sulfur oxidation
sulfide:quinone oxidoreductase
species
Table
Tris-Acetat-EDTA
Thermus aquaticus
Units
United Kingdom (Großbritannien)
United States of America (Vereinigte Staaten von Amerika)
Ultraviolett
Volt
5-Brom-4-chlor-3-indoxyl-ß-D-galactosi
3
4
Summary
(B) Summary
Microbial sulfur oxidation is an ancient metabolic process where reduced sulfur
compounds serve as electron donors for phototrophic and chemotrophic metabolism. The
role of sulfur-oxidizing prokaryotes (SOP) and the applied sulfur oxidation pathways in
marine sediments are poorly understood. Since SOP are phylogenetically diverse and are
often closely related to non-SOP, a solely 16S rRNA gene-based analysis is inadequate
for identification of SOP.
In this study the dissimilatory sulfite reductase (DsrAB), a key enzyme in both
SRP and also some SOP, was used as a phylogenetic marker gene for uncultured SOP.
For this purpose novel primer sets were developed and tested for applicability. In the
backbarrier tidal area (Janssand) of Spiekeroog Island in the German Wadden Sea sulfur
is extensively cycled between sulfate reduction and sulfur oxidation. From the respective
sediments three depth-specific clone libraries were generated. Phylogenetic analysis of 64
fully sequenced clones suggests five distinct sequence clusters of which 46 clones
affiliated with Gammaproteobacteria. In addition, one cluster of 18 clones displays a yet
unclear phylogenetic position but exhibits some similarities to Alphaproteobacteria. Most
sequences are only distantly related to known SOP indicating a large hidden diversity of
sulfur oxidizers in marine sediments. Comparative analysis of DsrAB and 16S rRNA
phylogeny of cultured SOP revealed a polyphyletic origin of the DsrAB sequences of the
alphaproteobacterial Magnetococcus sp. and Magnetospirillum spp.. These findings
strongly indicate the involvement of lateral gene transfer in the evolution of SOP.
Moreover, different treeing methods indicate a close evolutionary linkage between the
DsrAB of the bacterial sulfur oxidizing Chlorobium spp. and the archaeal sulfate reducer
Archaeoglobus fulgidus.
In addition, a novel quantitative real-time PCR assay using SYBR Green was
developed to quantify the dsrB copy number in Janssand sediments. Up to 1.85×107
copies/ml sediment were determined reaching a similar order of magnitude as previous
FISH counts of Gammaproteobacteria. The data imply that the “reverse” Dsr pathway
besides other pathways may play a major role in sulfur oxidation in Janssand sediments.
However, methodological improvements are necessary to enhance reproducibility of the
assay. Overall, the obtained results strongly support the relevance of
Gammaproteobacteria in sulfur oxidation in Janssand sediments.
Introduction
5
1. Introduction
Microbial sulfur oxidation is a key process for the biogeochemical sulfur cycle in
marine sediments and closely linked to the cycling of other elements like oxygen,
nitrogen, and carbon. The process is one of the oldest types of biological energy
conservation (Wachtershauser, 1990). Sulfur compound oxidizing prokaryotes are
phylogenetically (Figure 1) and metabolically (Table1) highly diverse (Lane et al., 1992;
Brock et al., 2006).
Figure 1: 16S rRNA tree of sulfur-oxidizing prokaryotes
Many bacteria and some archaea can carry out dissimilatory sulfur oxidation of
reduced inorganic sulfur compounds, e.g. sulfide, thiosulfate, polythionates, elemental
sulfur and sulfite (Brüser et al., 2000). The latter are employed as donors for electron
transport coupled to ATP synthesis or reduction of NAD(P)+ and other electron carriers.
Phototrophic and chemotrophic organisms exploit reduced sulfur compounds as source of
reducing equivalents (e-, H+) for assimilatory processes. In addition, chemolithotrophic
sulfur oxidizers gain energy from sulfur compound oxidation by coupling it to respiratory
electron transport.
Introduction
6
Sulfur chemolithotrophs grow primarily aerobic, i.e. using molecular oxygen as
terminal electron acceptor. However, some species (Beggiatoa sp., Thioploca sp.,
Thiobacillus denitrificans, Thiomicrospira denitrificans) oxidize H2S anaerobically
coupling it to nitrate reduction (Brock et al., 2006). SOP which exploit alternative
electron acceptors like iron or manganese have not been identified yet (Schippers and
Jorgensen, 2001, 2002). Further physiological differences can be found in the ability to
use various sulfur compounds. Consequently the biochemical sulfur oxidation pathways
are found to be variable.
Table 1: Metabolic lifestyles of sulfur oxidizers1
1
Term
Energy source
e- Donor
C source
organism
Chemolithoautotroph
chemical
H2S…
CO2
Paracoccus, Starkeya
Beggiatoa, Thiobacillus, Thiomargaritha,
Thioploca, Thiomicrospira,
Sulfolobus, Acidianus
Chemolithoheterotroph
(Mixotrophy)
chemical
H2S…
org. C
Beggiatoa, Thiothrix, Leucothrix
Photolithoautotroph
light
H2S…
CO2
Chlorobium, Chromatium, Allocromatium,
Thiocapsa Halochromatium, Rubrivivax,
Rhodobacter, Chloroflexus
Data collected from (Brüser et al., 2000; Brock et al., 2006)
The oxidation of sulfide, sulfur and polythionates is linked to thiosulfate oxidation
and can therefore easily be included into the existing thiosulfate oxidation pathways
(Brüser et al., 2000). So far several biochemical models for thiosulfate oxidation have
been proposed: the tetrathionate pathway (e.g. Thiobacillus thioparus, Acidithiobacillus
sp. and Acidiphilum sp.), the multi-enzyme-complex or SOX pathway (e.g. Paracoccus
denitrifcans, Starkeya novella) and the branched thiosulfate oxidation pathway (e.g.
Allochromatium vinosum, Thiobacillus denitrificans, Thiocapsa roseopercicina) (Brüser
et al., 2000).
The best studied archaeal sulfur oxidizer is the aerobic Acidianus ambivalens,
which has been shown to possess a sulfur oxygenase:reductase (SOR) and
thiosulfate:quinone oxidoreductase (TQO) (Sun et al., 2003; Muller et al., 2004; Urich,
2005). These enzymes catalyze the conversion of sulfur to sulfite, thiosulfate or hydrogen
sulfide and the oxidation of thiosulfate to tetrathionate, respectively (Figure 2).
Introduction
7
Figure 2: Model for sulfur oxidation in
Acidianus ambivalens. Elemental sulfur is
oxidized by a cytoplasmatic sulfur
oxygenase:reductase (SOR) yielding
thisosulfate, sulfite and hydrogen sulfide. The
latter are substrates for thiosulfate:quinone
oxidoreductase (TQO) and sulfite:acceptor
(SAOR) which are coupling the substrate
oxidation to energy conservation.
Adenylsulfate (APS) reductase and
adenylylsulfate: phosphate adenylyltransferase
(APAT), enzymes involved in the generation
of ATP from sulfite by substrate level
phosphorylation, have been identified.
Sulfide:quinone oxidoreductases (SQR), which
catalyzes the oxidation of sulfide activity has
so far not been shown, although homologues
to sqr genes can be found in the genome. CQ
Caldariella quinone (after Kletzin et al., 2004)
In the tetrathionate pathway of Acidithiobacillus and Acidiphilum spp. thiosulfate
is oxidized to the stable intermediate tetrathionate prior to cleavage by a hydrolase
(Brüser et al., 2000). For the oxidation of elemental sulfur and sulfide a different model
not involving tetrathionate has been developed (Figure 3) (Rohwerder and Sand, 2003).
Figure 3: Model for sulfur and sulfide oxidation in
Acidithiobacillus and Acidiphilum spp.. Extracellular
elemental sulfur (S8) is mobilized as persulfide sulfan
sulfur by outer membrane proteins (OMP) and
oxidized by a periplasmatic sulfur dioxygenase
(SDO). The resulting sulfite is further oxidized to
sulfate by sulfite:acceptor oxidoreductase (SOR).
Free sulfide is oxized to elemental sulfur by a
separate
sulfide:quinon
oxidoreductase.
Cyt
Cytochromes, Q Quinonens (after Rohwerder and
Sand, 2003)
Among the sulfur chemotrophs members of the genera Thiobacillus and
Paracoccus are best studied, whereas Allochromatium vinosum is the best studied
phototrophic sulfur oxidizer.
The “sulfur oxidation (SOX) pathway” or “Paracoccus sulfur oxidation (PSO)
pathway” (Kelly et al., 1997) is widespread among Alphaproteobacteria. The latter
harbors a complete thiosulfate-oxidizing multi-enzyme system (Kelly et al., 1997). The
encoding sox gene cluster (Friedrich et al., 2001) comprises 15 genes, seven of which
8
Introduction
(soxXYZABCD) encoding four essential proteins, and was first described for Paracoccus
panthotropus (Friedrich et al., 2005). Catalytic activities of the Sox proteins involve the
oxidation of sulfite, thiosulfate, sulfur and hydrogen sulfide coupled to cytochrome C
reduction without the formation of any free intermediate (Figure 4).
Figure 4: Model of thiosulfate oxidation in
Paracoccus panthotropus (Friedrich et al., 2001).
SoxAX catalyzes the oxidative and covalent transfer
of thiosulfate to the substrate carrier protein SoxYZ.
The SO3- is cleaved off by the hydrolase SoxB,
yielding sulfate and subsequently the S0 group is
oxidized to SO3- by SoxCD and again hydrolyzed by
SoxB (after Urich, 2005).
Sulfur-oxidizing members of the Beta- and Gammaproteobacteria as well as the
Chlorobiaceae (e.g Thiobacillus denitrificans, Allochromatium vinosum and Chlorobium
tepidum) possess an incomplete sox gene cluster and lack the sulfur dehydrogenase
(Friedrich et al., 2005). This points at a different pathway for oxidation of sulfur to
sulfate (Figure 5).
Figure 5: Schematic overview of key genes/enzymes
putatively associated with sulfur coumpound oxidation
in Thiobacillus denitrificans. The presence of sox and
dsr genes suggestes that different pathways enable the
oxidation of thiosulfate and sulfide to sulfate. Moreover
different enzymes may be involved in sulfur oxidation
under aerobic versus anaerobic conditions. (after Beller
et al., 2006)
Introduction
9
For phototrophs of the gamma-subclass and obligate chemotrophs of the betasubclass of Proteobacteria the branched thiosulfate oxidation pathway is established
(Figure 6). Dahl and Trüper (1994) proposed a general scheme for dissimilatory sulfur
metabolism in anaerobic phototrophic bacteria, which includes three essential steps: (1)
oxidation of sulfide or thiosulfate resulting in the formation of elemental sulfur, (2)
oxidation of sulfide or elemental sulfur to sulfite and (3) formation of sulfate as the final
product. Moreover, pioneering work on Allochromatium vinosum identified dissimilatory
sulfite reductase (DsrAB) to be involved in the anaerobic oxidation of sulfide and sulfur
(Dahl and Trüper, 1994; Pott and Dahl, 1998; Dahl et al., 2005). The respective enzyme
is well known from sulfate-reducing prokaryotes (SRP), mediating the reduction of
sulfite to sulfide (Hatchikian and Zeikus, 1983; Wolfe et al., 1994; Molitor et al., 1998).
The entire SOP dsr operon contains 15 genes, with dsrAB encoding an α2β2-structured
“reverse” operating DsrAB that catalyzes the oxidation of sulfide to sulfite (Figure 7).
Figure 6: Model of the branched thiosulfate oxidation pathway. SCS “sulfide carrier system” (after Brüser
et al., 2000)
10
Introduction
Figure 7: Two models for the interaction of Dsr proteins (Dahl et al., 2005) The reverse ctyoplasmatic
“DsrAB” oxidizes sulfide to sulfite. So far, unresolved remains how the substrate HS- is made available
from extracytoplasmatically stored sulfur and how electrons resulting from the sulfite reductase reaction are
fed into the electron transport chain. Model (1) Stored sulfur is reduced to the level of sufide ansd diffuses
as H2S or is transported in an unknown form and by an unknown mechanism over the cytoplasmatic membrane. DsrK and DsrM accept electrons from the Dsr feeding it into the photosynthetic electron flow. A
previously discovered potential NADP:acceptor oxidoreductase DsrL (not shown) could be involved in the
reductive release of sulfide from an organic perthiol carrier molecule. Model (2) DsrK is involved in the
reduction of the perthiol and obtains the required electrons from quinol via DsrM, releasing HS-. The model
can not explain the presence of the previously discovered membrane bound DsrJOP (not shown) and how
electrons from the sulfite reductase reaction are fed into the photosynthetic electron chain. (after Pott and
Dahl, 1998)
Besides the various enzymes and pathways for sulfur oxidation and its ecological
importance, so far, very little is known about the diversity and abundance of key
organisms driving sulfur oxidation in marine sediments. The polyphyly of SOP
complicates 16S rRNA based approaches. Comparative 16S rRNA sequence analysis
does not allow a conclusion on metabolic properties of the respective organisms. Hence,
only potential sulfur oxidizers can be revealed. The use of functional marker genes, i.e.
genes encoding proteins essential to a specific metabolic process, can overcome these
problems since it enables to focus on a defined physiological group. Phylogenies of key
enzymes involved in sulfur oxidation already exist for SoxB; (Petri et al., 2001),
sulfide:quinone oxidoreductases (SQR) (Theissen et al., 2003) and adenine
phosphosulfate (APS) reductase (Kuever, J., unpublished data; Rühland, C. unpublished
data).
Here, DsrAB was applied as functional marker for sulfur-oxidizing bacteria,
employing the “reverse” Dsr pathway (Schedel and Truper, 1979; Schedel et al., 1979).
The homologous protein of SRP is frequently used in environmental studies (Wagner et
al., 2005). So far an universal phylogenetic marker gene for all sulfur oxidizers is lacking.
However, the DsrAB is presumed to be the environmentally prevalent enzyme involved
Introduction
11
in sulfur compound oxidition (personal communication C. Dahl). Recently, the relevance
of DsrAB in unculturred SOP has been shown in Riftia pachyptila endosymbionts
(Markert et al., 2005), Olavius algarvensis gammaproteobacterial endosymbionts (Woyke
et al., unpublished data) and uncultured Beggiatoa sp. (Mußmann et al., unpublished
data).
Lateral gene transfer is regarded as a driving force in the evolution of metabolic
pathways among prokaryotes (Doolittle, 1999; Gogarten et al., 2002; Boucher et al.,
2003). It has frequently been reported for genes involved in sulfate reduction, eg. dsrAB
and apsA of APS reductase (Klein et al., 2001; Friedrich, 2002) and sulfur oxidation, e.g.
soxB, SQR genes and dsrAB (Petri et al., 2001; Theissen et al., 2003; Sabehi et al., 2005).
Thus detailed comparison of functional gene and 16S rRNA trees is necessary to track
putative events of horizontal transmission.
Since no fluorescence in situ hybridization probes are available that target all SOP
the determination of the in situ abundance of SOP is as yet limited to specific
phylogenetic groups such as Beggiatoa.(Mussmann et al., 2003), Thioploca (Schulz et
al., 1999), Arcobacter spp. (Llobet-Brossa et al., 1998). Moreover, available probes, e.g.
GAM660 and GSO477 target only subgroups of organisms affiliating with specific
clusters of endosymbiotic and free living hydrothermal vent associated
Gammaproteobacteria. Cultivation dependent most probable-number-counts (MPN)
revealed 0.93 to 9.30 × 105 cells/ml of chemolithoautotrophic sulfur-oxidizing bacteria in
Wadden Sea sediment. Thereof 0.93 to 9.30 × 103 cell/ml belonged to the genus
Thiomicrospira sp. (Brinkhoff et al., 1998). This numbers are probably underestimated
due to the selectivity of the cultivation approach. Thus, cultivation- as well as phylumindependent approaches to quantify SOP are essential.
Quantitative PCR (qPCR) represents a promising option. Here the copy number of
a certain marker gene is determined and thereof in situ cell number can be deduced. The
approach has already been applied for quantification of microorganisms in environmental
samples, targeting either 16S rRNA or functional marker genes (Suzuki et al., 2000;
Hermansson and Lindgren, 2001; Stubner, 2002; Kolb et al., 2003). Of course
methodological aspects like lyses efficiency in DNA extraction, PCR biases and
amplification efficiency, in addition to copy number of the respective target gene per cell
have to be considered. Consequently, qPCR is a semiquantitative method for the
assessment of bacterial abundances.
Marine sediments represent an ecosystem where extensive sulfur
transformations, including a large sulfur cycle of sulfate reduction and sulfide oxidation,
occur. The intertidal sand flats (Figure 8) at the backbarrier tidal area of Spiekeroog
Island (Janssand region, German Wadden Sea) offer an excellent sampling site for the
investigation of sedimentary SOP.
Moreover, intertidal sulfur pools (Figure 8) occur where elemental sulfur
granules accumulate in the residual water over time. This hints at a significant oxidation
of sulfide which outgases from the underlying sediment and might be catalyzed by SOP.
12
Introduction
Figure 8: Janssand sediments and sulfur pool, German Wadden Sea (M. Mußmann 2003/2005)
Previous 16S rRNA based community studies already reported sulfate reducing
Deltaproteobacteria, e.g. Desulfosarcinales and Desulfobulbaceae, to be abundant in
Wadden Sea sediments (Llobet-Brossa et al., 1998; Llobet-Brossa et al., 2002; Ishii et al.,
2004; Mussmann et al., 2005). Using quantitative CARD-FISH Ishii, Mussmann et al.
(2004) identified Gammaproteobacteria as the most abundant group in janssand
sediments. Furthermore, Arnds (2006) suggested a high abundance of potentially sulfuroxidizing Gammaproteobacteria and Roseobacter at the same site and a minor
significance of typical sulfur oxidizers such as Beggiatoa and Thiomicrospira spp..
Introduction
13
The aim of this study was to investigate the phylogeny and diversity of sulfuroxidizing bacteria employing the dissimilatory sulfite reductase (DsrAB) in Wadden Sea
sediments. For this purpose novel PCR primers specific for dsrAB of SOP were
developed. The applied primers were used to establish clone libraries from different
sediments depths and to amplify dsrAB from cultured sulfur oxidizers. Since the only
available phylogeny of DsrAB of SOP (Sabehi et al., 2005) comprises only sequences to
date, one goal was to generate a more comprehensive phylogeny by using a larger
sequence set. Moreover a phylogenetic comparison with DsrAB of SRP should be
facilitated.
In addition, a quantitative real-time PCR assay for the DsrB subunit was
developed to quantify the dsrB copy number per ml of sediment over the vertical profile.
This approach was supposed to give a first impression of the relevance of the DsrAB in
sulfur oxidation.
In detail this worked focused on the following questions:
1. What is the diversity of DsrAB of SOP in Janssand sediments?
2. What are the main groups of sulfur oxidizers?
3. Can novel groups of SOP be detected?
4. Is there a vertical stratification of certain phylotypes?
5. To what extent is the DsrAB diversity of SOP reflected on 16S rRNA level?
6. Has LGT also affected the evolutionary path of microbial sulfur oxidation?
7. How abundant are SOP employing the “reverse” Dsr pathway?
8. Is there a variation of dsrB copy number and hence cell number of SOP over the
vertical sediment profile?
9. Who are the putative key players mediating sulfur oxidation in Janssand
sediments?
Material and Methods
14
2. Material and Methods
In the following section the methodological procedure including all used material
and methods is described.
2.1 Material
Instruments:
•
Agilent 2100 Bioanalyzer
•
Autoclave Lam-201 (SANOclav, Bad Überkingen-Hausen)
•
Incubator Modell 800 (Memmert, Schwabach)
•
CYBERTECH CS1 Photodocumentation System (Cybertech, Berlin)
•
Demineralisation system MembraPure (MembraPure, Bodenheim)
•
Electrophoresis chamber (Horizon 58, GIBCO BRL, Life Technologies,
Gaithersburg, MD USA)
Kits:
•
TOPO TA Cloning Kit (Invitrogen, Karlsruhe)
•
QIAquick PCR purification kit (Fa. Quiagen, Hilden)
•
MinElute 96 UF PCR Purification Kit (Quiagen)
•
Plasmid MiniPrep96 Kit (Montage)
•
Plasmid MiniPrep96 Kit (Quiagen)
•
Power Soil DNA Kit (MO Bio Laboratories, Inc.,CA, USA)
PC-Software:
•
Agilent 2100 Bioanalyzer
•
ARB 2.4
•
BioEdit
•
Canvas 8.0.4(Microsoft)
•
ClustalW 1.8.1
•
Endnote 7
•
Excel (Microsoft)
•
Photoshop (Microsoft)
•
Sequencher 4.5.1 (GeneCodes)
•
Sequencing Analysis 5.2 (Applied Biosystems)
•
Word (Microsoft)
Material and Methods
15
•
Plates
•
OmniTray Singlewell 86 x 128 mm (Nalge Nunc International, New York, USA)
for preparation of agar plates for 96 clones
•
96-well Standard Microtiterplate (ABgene House, Surrey, UK) for long term
storage of 96 clones at 24 °C
•
Thermo-Fast® 96 Skirted mit Domed Cap Strips (ABgene House, Surrey, UK)
•
Mikroplatte 96 V (Roth, Karlsruhe)for PCR with fluorescence dye SybrGreen
•
MultiScreen®-HV (Millipore, Molsheim, Frankreich) for parallel purification of
96 Cyle Sequencing PCR products
•
Shaker KS 250 basic (IKA Labortechnik, Wilmington, USA)
•
Sequencer 3130xl Genetic Analyzer (Applied Biosystems, Foster City, USA)
•
Sequencer ABI Big Dye ABI PRISM 3100 genetic analyzer (Applied Biosystems,
Foster City, California, USA)
•
Thermocycler for PCR and Cycle Sequencing PCR
•
Mastercycler (Eppendorf, Hamburg)
•
Mastercycler gradient (Eppendorf, Hamburg)
•
Mastercycler personal (Eppendorf, Hamburg)
•
Vortexer MS1 Minishaker (IKA Labortechnik, Wilmington, USA)
•
Watherbaths
•
Typ 1086 (Memmert, Schwabach)
•
Typ 1002 (Memmert, Schwabach)
•
Centrifuges
•
Centrifuge 5810R (Eppendorf, Hamburg)
•
Centrifuge 5415R (Eppendorf, Hamburg)
•
UV-Transilluminator (Vilber Lourmat, Marne la Vallée, Frankreich)
2.2 Chemicals
•
Acetic acid (glacial) (Fluka, Buchs, Schweiz)
•
Agar-Agar (Difco, Heidelberg)
•
Agarose (Seakem, LE, Cambrex Bio Science Rockland, Inc. Rockland, ME USA)
•
Agarose, low melting (Metaphor, Southborough, USA)
Material and Methods
16
•
Ampicillin Sodiumsalt (Roth, Karlsruhe)
•
5-Brom-4-chlor-3-indolyl-ß-D-galactosid [X-Gal] (Roth, Karlsruhe)
•
Bromphenol-Blue
•
BSA solution (bovine serum albumine, Sigma, Taufkirchen)
3mg/ml in PCR water
•
Casein (Sigma-Aldrich, Steinheim)
•
3.8-Diamino-5-ethyl-6-phenylphenanthridiniumbromid
[Ethidiumbromid] (Sigma-Aldrich, Steinheim)
•
dNTP solution (Roche, Mannheim)
•
Dinatriumhydrogenphosphat Heptahydrat (Fluka, Buchs, Schweiz)
•
DNA molecular weight markers
•
96% Ethanol (Merck, Darmstadt)
•
ExoSap-IT enzymes (Bioscience, Amersham)
•
Glycerol (Fluka, Buchs, Schweiz)
•
Isopropyl-ß-D-thiogalactosid [IPTG] (Promega, Mannheim)
•
Low DNA mass ladder (Invitrogen, Breda, Netherlands)
•
Primer 50 pmol/µl (Biomerse, Ulm)
•
LiChrosolv Chromatographie Wasser (Merck, Darmstadt)
•
Na2EDTA
•
PCR Water (Sigma-Aldrich, Steinheim)
•
SephadexTM G-50 Superfine (Amersham Bioscience, Freiburg)
•
Sucrose
•
Tris (Biomol, Hamburg)
•
Taq Polymerase 5U/µl (Eppendorf, Hamburg)
•
Taq reaction buffer (10x) (Eppendorf, Hamburg)
•
Tris (Biomol, Hamburg)
•
Xylene Cyanol
•
Yeast extract (Fluka, Buchs, Schweiz)
Material and Methods
17
2.3 Sampling
An intertidal sandflat in the backbarrier tidal area of Spiekeroog Island (German
Wadden Sea, Janssand region 53°43`45``North and 7°41`80``East) served as sampling
site (Figure 9). The respective area of investigation was already described in detail in a
previous study, for comprehensive information it is referred to Arnds (2006).
North Sea
Janssand
Spiekeroog
Germany
Figure 9: Satellite image of the sampling site Janssand in the backbarrier tidal area of Spiekeroog Island
Sediment cores for the SOP diversity study were taken in Novemer 2004 and
April 2005 and DNA was extracted by M. Mußmann from three distinct depths. In
addition, an intertidal pond with precipitated sulfur was sampled in April 2005 and DNA
was extracted from the pond water. A clone library was constructed from the surface
layer 0 – 3 cm (sediment core 02/04/2005) as well as from the deeper layers 8 – 19 cm
and > 19 cm (sediment core 02/11/2004). In addition a sediment core was taken for the
quantitative survey of SOP in November 2005 frozen and cut into slices. Sediment
samples were taken in triplets from 7 distinct sediment depths, including oxic and anoxic
sediment layers and the oxic-anoxic interface as indicated by the sediment colour (Figure
10) To allow calculation of copy number per ml sediment, each sediment sample was
weighed and the mass of 0.5 ml reference sediment was determined for each depth. DNA
was extracted using the Power Soil DNA Kit (MO Bio Laboratories, Inc.,CA, USA)
according to the manufactoresrs protocol.
Material and Methods
18
LAYER
SAMPLED DEPTH [cm]
OXIC (BROWHNISH)
0.0 – 0.5
1.0 – 1.5
3.0 – 3.5
3.5 – 4.0
OXIC - ANOXIC
ANOXIC (GREYISH)
6.0 – 6.5
11.0 – 11.5
Figure 10: Sediment core for the quantitative survey of SOP taken in November 2005
2.4 Primer design
Design of appropriate primers was the precondition to facilitate the amplification
of dsrAB and dsrB from SOP. Therefore, published (www.ncbi.nlm.nih.gov/Genbank)
and unpublished dsrAB sequences (Table 2) of exclusively bacterial origin were aligned
in BioEdit using ClustalW algorithm. No archaeal dsrAB sequences were available.
Based on the sequence alignment (see appendix Figure A1 and A2) new, degenerated
primers (Table 3).were manually designed for target sites of high sequence conservation
(M. Mußmann, unpublished).
Table 2: Overview of dsrAB sequences used for primer design
Organism
Source
Allochromatium vinosum
GenBank
Alkalilimnicola ehrlichei
GenBank
Chlorobium limicola
GenBank
Chlorobium phaeobacteroides BS1
GenBank
Chlorobium phaeobacteroides DSM
GenBank
Chlorobium tepidum
GenBank
Magnetococcus sp.
GenBank
Magnetospirillum magneticum
GenBank
Magnetospirillum magnetotacticum
GenBank
Pelodictyon phaeoclathratiforme
GenBank
Prosthecochloris aestuarii
GenBank
Prosthecochloris vibrioformi
GenBank
Thiobacillus denitrificans
GenBank
Thiobacillus denitrificans
GenBank
Thiobacillus denitrificans
GenBank
Sargasso Sea uncultured bacterium
GenBank
uncultured bacterium MED1313K9
GenBank
uncultured bacterium Fosmid194A5
Dr. M. Mußmann
Olavius algarvensis gamma1 symbiont
Dr. N. Dubilier
Olavius algarvensis gamma 2 symbinot
Dr. N. Dubilier
Accession number
AAC35394/AAC35395
EAP35273/EAP35272
EAM42077/EAM42078
EAM62784/EAM62785
EAM34529/EAM34528
AAM72088/AAM72089
EAN29021/EAN29022
BAE52171/BAE52172
ZP_00053119/ZP_00053120
EAN24264/EAN24263
EAN22044/EAN22045
EAO15265/EAO15264
AAZ98438
AAZ97262
AAZ97322
EAJ64896/EAJ64897
AAY89969/AAY89968
unpublished
unpublished
unpublished
Material and Methods
19
Table 3: Overview of primers designed for dsrAB of SOP
Binding site
(Chlorobium limicola)
Primer1
Target
Sequence2 (5´
dsrA 116F
Bacteria
115 - 131
GGNCCNTGGCCBAGYTT
dsrA 240F
Bacteria
226 - 242
GGNTAYTGGAARGGYGG
dsrA 240F*
Bacteria
226 - 242
GGNTAYTGGAARGGNGG
dsrB 403F
Bacteria
385 - 401
CAYACNCARGGNTGGY
dsrB 403R
Bacteria
385 - 401
ARCCANCCYTGNGTRTG
dsrB 808R
Bacteria
781 - 797
CCDCCNACCCADATNGC
3´)
1
Reference Mußmann, unpublished; Names according to alignment positions of Chlorobium limicola.
Ambiguity base codes: R = A, G; Y = C, T; B = C, G, T; N = A, T, C, G
F – Forward primer (binds to 5´ end of the target region)
R – Reverse primer (binds to 3´ end of the target region)
2
2.5 Amplification of dsrAB from Wadden Sea sediments
2.5.1 Establishment of dsrAB clone libraries
Polymerase chain reaction
Polymerase chain reaction after (Saiki et al., 1988) was used for in vitro selective
amplification of the dsrAB of sulfide-oxidizing bacteria. This method is based on the use
of two flanking, single stranded oligonucleotid primers (~20 nt) which bind specific to
their target sequences at both ends of the amplifiable region and a thermostable Taqpolymerase. The principle of repetitive cycles, where newly synthesized strands serve as
template during the following cycles, already allows the amplification of minimal
amounts of DNA.
Material and chemicals
•
Taq reaction buffer (10x) (Eppendorf, Hamburg)
•
dNTP solution (Roche, Mannheim)
o 2.5 MM dATP
o 2.5 MM dCTP
o 2.5 MM dGTP
o 2.5 MM dTTP
•
BSA solution (bovine serum albumine, Sigma, Taufkirchen)
Material and Methods
20
o 3mg/ml in PCR water
•
Taq Polymerase 5U/µl (Eppendorf, Hamburg)
•
Primer 50 pmol/µl (Biomerse, Ulm)
•
PCR water (Sigma, Taufkirchen)
Procedure
Standard reaction of 25 µl volume for amplification of dsrAB
•
Taq buffer (10x)
2.5 µl
•
BSA (3mg/ml)
2.5 µl
•
Dntps (10mM)
2.5 µl
•
Primer F (50µM)
1.0 µl
•
Primer R (50µM)
1.0 µl
•
Taq Polymerase (5U/µl)
0.1 µl
•
Template
1 µl
•
PCR-H2O
ad 25 µl
Table 4: Program for amplification of dsrAB from environmental samples
Phase of reaction
Duration of phase [min]
Temperature [°C]
Number of cycles
Initial denaturation
4
95
1
Cycle denaturation
0.5
95
30
Annealing
0.5
55.5
30
Elongation
2
72
30
Final elongation
10
70
1
Primer test by Temperature gradient PCR
To find the optimal annealing temperature of the new, degenerated primer pair
dsrA 240F and dsrB 808R a temperature gradient PCR was performed. Positive dsrAB
PCR products from environmental DNA JSS 0-3cm and resuspended cells from
Thiobacillus denitrificans (Betaproteobacteria) served as template. PCR reactions were
prepared according to the PCR standard reaction protocol containing 0.1 µl of template
(dsrAB PCR product from environmental DNA JSS 0-3cm) and 0.3 µl of template (dsrAB
PCR product from Thiobacillus denitrificans, 1:10 dilution).
Material and Methods
21
Amplification of dsrAB from environmental DNA
The dsrAB fragment was amplified using the degenerated primer pair dsrA 240F
and dsrB 808R specific for dsrAB of sulfide-oxidizing bacteria.
Sediment as well as sulfur pond samples were amplified in five replicate reactions
to minimize stochastic PCR bias (Cavanaugh et al., 1998). Each 25 µl reaction was
prepared according to the PCR protocol (see Table ) containing 2 µl of template DNA for
the surface sediment sample JSS 0-3cm and accordingly 1µl of template for the two
deeper sediment samples JSS 8-19cm and JSS > 19cm. Amplification was carried out in
an Eppendorf Mastercycler using the conditions shown in Table 4. After completion
samples were cooled down to 12 ° C and stored at 4°C if not immediately processed
futher. The replicates of the PCR products were combined and 40 µl of each amplificate
were applied to an 1.3% agarose gel for excision of band.
Agarose gel electrophoresis
Material and chemicals
TAE buffer (50x) stock solution
•
Tris
2.0 M
•
Acetic acid (glacial)
0.95 M
•
Na2EDTA (pH 8.0)
0.5 M
The TAE buffer (1x) working solution has to be diluted 10 fold with Milli-Q.
DNA molecular weight marker
•
Low DNA mass ladder (Invitrogen, Breda, Netherlands)
Agarose gel (1.3%)
•
1.3 gram Agarose per 100ml TAE (1x) working solution
Loading buffer (6x)
•
Bromphenol-Blue (0.25%)
2.5 g/l
•
Xylene Cyanol (0.25%)
2.5 g/l
•
Sucrose (40%)
400 g/l
Ethidiumbromid staining bath
•
Ethidiumbromid
0.4 µg×ml-1 in Milli-Q
Material and Methods
22
Procedure
The gel electrophoresis was used for analysis of the fragment size. Electrophoretic
separation of nucleic acids was performed according to Sambrook et al. (2001). Samples
were prepared with > 0.2 vol.-% loading buffer and applied to an 1.3% agarose gel
covered with TAE (1x) in a horizontal electrophoresis chamber (Horizon 58, GIBCO
BRL, Life Technologies, Gaithersburg, MD USA). In addition, 2µl of DNA molecular
weight markers mixed with µl loading buffer were applied for size determination and
quantification of the fragments. For spatial separation an electrical field was applied,
according to a voltage power of 4 V/cm distance anode – cathode, on average 100 – 200
mV, for 10 to 20 min. Afterwards the agarose gel was stained for 10 to 20 min in an
Ethidiumbromid staining bath. Detection of the nucleic acids was carried out via a UVTransilluminator (Vilber Lourmat, Marne la Vallée, Frankreich). Stained agarose gels
were documented with the CYBERTECH CS1 Photodocumentation System (Cybertech,
Berlin).
Gel electrophoresis yielded gene products of different sizes including an
approximately 2kb DNA fragment for the three distinct sediment samples. The latter
band was excised and DNA was purified from the gel using the QIAquick PCR
purification kit (Fa. Quiagen, Hilden) according to the instructions of the manufacturer.
Amplification of dsrAB from the sulfur pond yielded no visible band in the agarose gel.
Cloning
Material and chemicals
TOPO-TA Cloning Kit (Invitrogen, Karlsruhe) containing
•
pCR®4-TOPO® Vector (10 ng µl-1) 3956 bp
o Triton X-100
0.1%
o Glycerol
50%
o Tris-HCl (pH 7.4)
50 mM
o EDTA
1 mM
o DTT
2 mM
o BSA
100 µg ml-1
o Phenol red
•
•
salt solution
o Sodium chloride
1.20 M
o Magnesium chloride
0.06 M
strain Escherichia coli TOP10
genotyp: lacZΔM15 hsdR lacX74 recA endA tonA
Material and Methods
•
23
SOC-medium (pH 7.3)
o Trypton
2.0%
o Yeast extract
0.5%
o Sodium chloride
10 mM
o Potassium chloride
2.5 mM
o Magnesium chloride
10 mM
o Magnesium sulfate
10 mM
o Glucose
20 mM
TOPO-TA Cloning ligation reaction
•
Salt
1µl
•
PCR product
2µl
•
sterile
3µl
•
PCR 4 – TOPO vector
1µl
Procedure
Ligation and transformation
PCR products of each depth were cloned separately using the TOPO TA Cloning
Kit (Invitrogen, Breda, Netherlands) according to the manufacturer´s protocol. In brief,
two microliters of the purified PCR product were ligated into pCR 4 – TOPO vector
incubated at room temperature for 30 min. Afterwards 2 µl of the ligation reaction were
added into one vial of chemically competent Escherichia coli cells. Incubation time on
ice was prolonged to 30 min to ensure a sufficient attachment of vector onto the cells.
Transformation was achieved via heat shock. Finally 40 and 80 µl of each transformation
reaction were spread on selective LB plates and incubated overnight at 37°C.
Selection of recombinant clones
•
LB/AMP/IPTG/x-Gal selective agar plates
•
LB medium liquid
•
LB medium liquid + 10 vol.-% glycerol
•
Ampicillin stock solution
100 µg µl-1 Ampicillin (Roth, Karlsruhe) were dissolved in MQ - water gelöst,
sterile filtered (pore size 0.2 µm) and stored at -20°C -20°C
Material and Methods
24
•
0.5 M IPTG stock solution
Isopropyl-ß-D-thiogalactosid [IPTG] (Promega, Mannheim) were dissolved in
50% ethanol, sterile filtered (pore size 0.2 µm) and stored at -20°C
•
X-Gal stock solution 100 mg ml-1 5-Chlor-4-brom-3-indol-ß-D-galactosid [X-Gal]
(Roth, Karlsruhe) were disolved in N, N-dimethyl-formamid and stored at -20°C
•
LB-Medium (pH 7.3-7.5)
o 1.0%
Casein
o 0.5%
Yeast extract
o 0.5%
Sodiumchlorid
o 1.5%
Bacto-Agar (for solid LB)
All reagents were dissolved in MQ-water and autoclaved for 45 min. Solid LB-medium
was taken out at a temperature of 80 °C to keep the solution liquid for pouring. For
plating of clones 100 µg ml-1 ampicillin, 20 µg ml-1 IPTG und 20 µg ml-1 X-Gal were
added to the LB-medium. Solutions were stored autoclaved and at 4 °C.
Cells containing plasmid inserts were selected using the blue white screening
method (Sambrook et al., 2001). After final incubation at 4°C white colonies were
transferred on new selective LB agar plates and simultaneously into 96 well-format plates
(NunclonTM, Nunc, Wiesbaden) containing 150 µl liquid LB medium. Both were
incubated overnight at 37°C, the liquid cultures on a rotary shaker at 200 rpm. For
permanent storage clones from liquid culture were transferred into each 150 µl liquid LB
containing glycerol, incubated overnight on a rotary shaker at 200 rpm and afterwards
frozen at -20°C.
Screening of clones
Material and chemicals
PCR reagents see point 2.5.1
Procedure
Screening for the dsrAB fragment was carried out by PCR using the M13 forward
and reverse vector primers (see Table ). 20 clones of each depth were randomly picked
and controlled for the correct insert size (~2 kb). Each 25 µl reaction was prepared
according to the PCR protocol (see Table ) containing 0.3 µl of primer mix and 1 µl of
liquid LB culture as template. Afterwards gel electrophoresis through a 1.3% agarose gel
was performed.
Material and Methods
25
Table 5: PCR Primers used for amplification of the dsr inserts in E. coli
Primer
M13F
M13R
dsrB 403R
Target/
Binding
site
Vector
355 - 378
Vector
199 - 222
Bacteria
403 - 419
Sequence (5´
3´)
Referenz
ACGACGTTGTAAAACGACGGCCAG
TTCACACAGGAAACAGCTATGACC
Yanisch-Perron
et al. (1983)
Mußmann,
unpublished
ARCCANCCYTGNGTRTG
F – Forward primer (binds to 5´ end of the target region)
R – Reverse primer (binds to 3´ end of the target region)
Table 6: Program for amplification of insert from pCR 4 – TOPO vector
Phase of reaction
Duration of phase [min]
Temperature [°C]
Number of cycles
Initial denaturation
4
95
1
Cycle denaturation
0.5
95
30
Annealing
30
50
30
Elongation
2
72
30
Final elongation
10
70
1
Purification of PCR products
Material and chemicals
•
•
ExoSap-IT enzymes (Bioscience, Amersham)
MinElute 96 UF PCR Purification Kit (Quiagen)
Procedure
PCR products of correct size were purified using either the MinElute 96 UF PCR
Purification Kit (Quiagen) according to he manufacturers instructions or ExoSap
purification protocol. Therefore 0.25 of ExoSap-IT enzymes and 3.25 PCR water were
combined with 1.5 µl of PCR reaction and incubated in an Eppendorf Mastercycler at
first for 30 minutes at 37°C and afterwards for 15 minutes at 80°C. 3 µl of ExoSap
reaction were used as template for sequencing.
Since neither of the PCR purification methods yielded a product from which
sequences could be obtained plasmids were extracted for sequencing of dsrAB.
Material and Methods
26
Plasmid preparation
Material and chemicals
•
Plasmid MiniPrep96 Kit (Montage)
•
Plasmid MiniPrep Kit (Quiagen)
•
2x LB medium liquid
Procedure
15 µl of liquid LB culture were transferred into 1 ml of liquid 2x LB and
incubated overnight at 37°C on a rotary shaker at 200 rpm. To obtain a compact pellet,
cells were centrifuged at 1500 g for ten minutes. LB medium was completely removed
and cells were processed further according to the Plasmid MiniPrep96 protocol (Montage)
for clones from environmental sample JSS 8-19 cm and JSS >19 cm or Plasmid MiniPrep
Kit (Quiagen) for clones from environmental sample JSS 0-3 cm.
2.5.2 Sequencing
Material and chemicals
PCR reagents see point 2.5.1
Procedure
Standard sequencing reaction of 5 µl volume for dsrAB
•
Big Dye
1.0 µl
•
Buffer
1.0 µl
•
Primer F (50µM)
0.1 µl
•
Primer R (50µM)
0.1 µl
•
Template
2.5 µl
•
PCR-H2O
ad 5 µl
Material and Methods
27
Table 7: Program for sequencing
Phase of reaction
Duration of phase
Temperature [°C]
Number of cycles
Initial denaturation
20 sec
96
1
Cycle denaturation
10 sec
5 sec
30sec
4 min
96
> 50
≤ 50
60
60
60
30 sec
21
1
Annealing
Elongation
60
Sequencing of the complete dsrAB inserts was achieved by using the M13 F and
M13R vector primers and in addition, the internal new, degenerated dsrB 403R primer.
For the latter the optimal annealing temperature was detected by performing a
temperature gradient PCR of a dsrAB insert containing clone. Vector primed sequencing
with the M13F and M13R was performed at an annealing temperature of 48 °C and 57°C
respectively for plasmids from sediment samples JSS 8-19 as well as JSS >19 and JSS
0-3 cm respectively. Sequencing with the dsrB 403R primer was performed at an
annealing temperature of 57 °C for all samples.
2.5.3. Phylogenetic analysis
The dsrAB sequences were checked against sequences in the Genbank database
using BLAST for similarity searches (Altschul et al., 1990). Partial sequences were
analyzed for sequence quality using the Sequencing Analysis software package (Applied
Biosystems). Full length sequences were generated by assembling partial sequences with
the Sequencer Version 4.5.1 software package. The dsrA and dsrB clone sequences were
aligned separately together with the respective published sequences using the ClustalX
and Bioedit software and subsequently translated into proteins. Further phylogenetic
inference of the DsrAB sequence data was performed with the ARB software package
(www.arb-home.de).
102 full length sequences from the environmental clone libraries (64 different
sequences), cultivated (14 sequences) and uncultivated (4 sequences) bacteria as well as
cultivated representatives of SRP (20 sequences (Klein et al., 2001)) and 3 partial length
sequences of uncultured SOP were entered into the dsrAB ARB database, translated into
amino acids and computational and manually aligned. In brief, first SOP DsrAB
sequences were roughly aligned against the existing SRP alignment. Afterwards the
resulting SOP DsrAB alignment was corrected on the basis of the Thiobacillus
denitrificans sequences. Finally, 20 representatives of SRP were aligned against the
refined SOP DsrAB alignement. A total number of 532 amino acid positions (alpha
subunit 318; beta subunit 214) were included in DsrAB analysis (excluding insertions and
deletions). Phylogenetic trees were calculated by performing Neighbour Joining,
Maximum Parsimony, Maximum Likelihood and Tree Puzzle analysis. For tree
reconstruction only full length sequences were used. Partial length DsrAB sequences of
Beggiatoa, Riftia pachyptila sulfur oxidizing endosymbiont and Olavius algarvensis
sulfur oxidizing gamma 1 symbiont were added to the final trees using the ARB
Parsimony Qick Add tool. To root the tree the 20 representatives of cultivated SRP (Klein
Material and Methods
28
et al., 2001) served as outgroup. In addition to the inference methods mentioned above,
we employed MrBayes phylogenetic analysis starting from the Neighbour joining tree.
Calculation parameters included 30000 generations, a print frequency of 100 and a
sample frequency of 10 resulting in a total number of 3000 trees. For final tree
construction the first 2000 trees were discarded.
2.6 Amplification of dsrAB from cultivated SOP
For completion of the phylogenetic tree dsrAB fragments should be amplified
from sulfur-oxidizing reference strains (Table 8)
Table 8: Sulfur oxidizing bacteria chosen for amplification of dsrAB
Species
Strain no/Referenz
Beggiatoa sp.
Susanne Hink
Halochromatium salexigens
DSM4395
Halothiobacillus neapolitanus
DSM581
Leucothrix mucor
Thomas Holler
Marichromatium gracile
DSM203
Magnetospirillum gryphiswaldense
PD Dr.Dirk Schüler
Rhodobacter capsulatus
DSM1710
Rubrivivax gelatinosus
DSM6859
Starkeya novella
DSM506
Thiobacillus thioparus
DSM505
Thiobaca trueperi
DSM13587
Thiocapsa sp.
DSM5653
Thiocococcus pfennigii
DSM226
Thiomicrospira psychrophila
Dr. Katrin Knittel
Thiomicrospira arctica
Dr. Katrin Knittel
Thiomicrospira kuenenii
DSM12350
Thiomicrospira frisia
DSM12351
Thiomicrospira pelophila
DSM1534
Phylum
Gamma-Proteobacteria
Gamma-Proteobacteria (PSB)
Gamma-Proteobacteria
Gamma-Proteobacteria
Gamma-Proteobacteria
Alpha-Proteobacteria
Alpha -Proteobacteria (PNSB)
Beta-Proteobacteria (PNSB)
Beta-Proteobacteria
Beta-Proteobacteria
Gamma-Proteobacteria
Gamma-Proteobacteria (PSB)
Gamma-Proteobacteria (PSB)
Gamma-Proteobacteria
~
~
~
~
Primer test by Temperature gradient PCR
Standard reaction of 25 µl volume according to point 2.5.1
•
Primer F (50µM) 16s/dsr
0.25/ 1.0 µl
•
Primer R (50µM) 16s/dsr
0.25/ 1.0 µl
Table 9: PCR Primers used for amplification of dsrAB from cultured SOP
Primer
Target
Binding site Sequence (5´
dsrA 240F*
Bacteria
240 - 256
dsrB 808R
Bacteria
808 - 824
3´)
Referenz
GGNTAYTGGAARGGNGG Mußmann,
unpublished
CCDCCNACCCADATNGC
F – Forward primer (binds to 5´ end of the target region)
R – Reverse primer (binds to 3´ end of the target region)
* modified dsrA 240 F primer containing an additional wobble at position 15
Material and Methods
29
For amplification of the dsrAB fragment from cultured SOP the new degenerated
dsrA 240F* primer was designed on the basis of an additional dsrAB sequence from
fosmid clone 194A5 (unpublished, M. Mußmann). Here, an additional degeneracy was
introduced in the previously designed primer dsrA 240F, making it fully complementary
to the respective target site. To find the optimal annealing temperature of the new dsrA
240F* and dsrB 808R primer pair an annealing temperature dependent amplification of
the dsrAB gene fragment from Magnetospirillum gryphiswaldense was performed. Cells
from pure culture were frozen and thawn and served as template for PCR using the PCR
standard reaction protocol and conditions under 2.5.1 without a final elongation.
Amplification of 16s rRNA gene and dsrAB gene fragment from reference strains
Cell pellets of the freeze-dried reference strains (obtained from DSMZ ) were
directly resuspended in 30 – 60 µl PCR water (Sigma, Taufkirchen). Living cultures were
centrifuged (Eppendorf Centrifuge 5810) at 1500 g for 10 minutes to gain a sufficient cell
pellet. Amplification of dsrAB gene fragment from the reference strains was performed
according to 2.5.1 using an annealing temperature of 57.5 and 51.3 °C respectively.
Magnetospirillum gryphiswaldense served as positive control.
To ensure the right DNA amount for potentially successful amplification of DNA
from target cells and to proof the correct template concentration, the 16s rRNA gene was
amplified simultaneously from all reference strains. Besides using the 16s rDNA primers
(Table 10), a primer specific annealing temperature (Table 11), the standard reaction
protocol (see above) and PCR conditions were analogue to dsr amplification.
Table 10: Primers used for amplification of the 16s rRNA gene
Primer
Target
Binding site Sequence (5´
3´)
8F
16s rDNA
8-24
AGAGTTTGATCMTGGC
1492 R
16s rDNA
1492-1507
TACCTTGTTACGACTT
Referenz
Muyzer et al.
1995
F – Forward primer (binds to 5´ end of the target region)
R – Reverse primer (binds to 3´ end of the target region)
Table 11: Program for amplification of dsrAB from cultured SOP
Phase of reaction
Duration of phase [min]
Temperature [°C]
Number of cycles
Initial denaturation
4
95
1
Cycle denaturation
0.5
95
35
Annealing 16s/dsr
0.5
48.0/ 57.5; 51.3
35
Elongation
2
72
35
Final elongation
10
70
1
30
Material and Methods
2.7 Specificity test of dsrAB primers dsrA 240 F and dsrA
240 F*
To test the specificity of the dsrA 240 F and dsrA 240 F* primer with regard to
the additional degeneracy PCR was performed on the Fosmid194A5 containing the dsrAB
gene fragment according to the protocol under point 2.5.1. using 0.5 µl of fosmids as
template and 35 PCR cycles. A dsrAB containing plasmid served as positive control.
2.8 Determination of dsrB copy number in Wadden Sea
sediments
To quantify the copy number of dsrAB from SOP over the vertical sediment
profile qRT -PCR based on a SybrGreen assay was employed. This method allows for
simultaneous amplification and quantification of a target template by detecting the
amplification associated fluorescence. Using PCR for quantification seems contradictory
at the first glance since this method was originally developed for massive amplification of
smallest DNA amounts, i.e. the pure qualitative detection. Nevertheless, it is based on a
exponential amplification which follows mathematical rules. Thus a quantification is
principle possible (Mülhardt, 1999).
The applied SYBR Green assay is generally based on the detection and
quantification of the respective fluorescence reporter. SYBR Green is a fluorescence dye
which intercalates with double stranded DNA. The fluorescence of SYBR Green
increases when bound to accumulating amplification product. Consequently, a higher
fluorescence signal will be detected during each subsequent PCR cycle. Since the
reaction is monitored continuously an increase in fluorescence signal can be viewed at
real-time. The parameter used for quantification, is the threshold cycle, or Ct-value,
namely, the cycle at which a significant increase in ∆Rn fluorescence/PCR product can
be observed. (Dorak). In general, homogeneous amplification of the target gene is a
precondition for reliable quantification. Since the amplification efficiency is higher for
shorter templates (Mülhardt, 1999), we amplified a 400 bp sized gene fragment of the
DrsB subunit. An absolute quantification assay, i.e. copy number was quantified by
interpolation from a standard curve, was employed for determination of the dsrAB copy
number in the environmental samples.
2.8.1 Primer specificity check
As a consequence of the unspecific binding of SybrGreen to dsDNA the measured
fluorescence reports only an increase in DNA but not necessarily an increase in specific
product. Hence specific amplification of the target gene has to be ensured. To identify
possible contaminations gel electrophoresis, sequencing and melting curve analysis were
employed.
Material and Methods
31
2.8.2 Standard Generation
For generation of an appropriate standard a dsrB clone library was established.
Therefore a dsrB gene fragment of approximately 400bp was already amplified from
environmental DNA JSS 0-3 cm using the degenerated primer pair dsrB 403F and dsrB
808R by Marc Mußmann. The PCR was performed according to the protocol under 2.5.1
using an annealing temperature of 52.9°C. Prior to further processing the stored product
was incubated under the conditions shown in Table 12 to restore sequence quality at the
end of the sequences. The PCR product was stored at -18°C.
Table 12: PCR products Cycle Program for incubation of stored
Phase of reaction
Duration of phase [min]
Temperature [°C]
Denaturation
4
95
Elongation
10
72
Final elongation
60
60
For further library construction steps listed under 2.5.1 were performed. Since
only a single band was observed the PCR product was cloned directly without excision
from the agarose gel. Vector primed sequencing of 16 clones, containing the correct size
insert, was performed using the M13 F primer according to 2.5.2.
Obtained sequences were checked against sequences in the Genbank database
(www.ncbi.nlm.nih.gov/blast/) using BLAST (Altschul et al., 1990) for similarity
searches.
Determination of copy number in standard
For generating a standard curve the respective ~ 400bp sized dsrB fragment was
amplified from a clone according to 2.5.1. The resulting PCR product of approximately
570bp was purified using the MinElute 96 UF PCR Purification Kit (Quiagen) and served
as template for serial dilutions. For the quantification of copy number in the standard the
fragment size [bp] and the DNA amount [ng/µl] was determined via the 2100 Bio
analyzer and Bio Sizing software Version A.02.12 (Agilent Technologies, USA).
Table 13: Determination of copy number in standard
Standard
Fragment size
[bp]
Concentration
[ng/µl]
Copy number
[mol of dsDNA/µl]
Copy number
per µl
sample 1
569
10.7
2,85052E-14
1,72E+10
sample 2
574
10.8
2,84895E-14
1,72E+10
Conversion of µg to pmol:
pmol of dsDNA = μg (of dsDNA) × 1515
fragment size in bp
(1)
Material and Methods
32
Conversion of pmol to absolute copy number:
mol DNA × 6, 02214 × 10 23
(2)
A two fold dilution series, comprising 11 dilution steps, was prepared starting
from a 1000fold diluted purified PCR product.
2.8.3 Quantitative Real-time PCR
Material and chemicals
•
iCycler iQTM5 Multicolour Real-Time PCR Detection systems (Bio-Rad, CA,
USA)
•
Platiunum SYBR Green qPCR Super Mix-UDG Kit (Invitrogen, Kalsruhe)
containing:
•
Platinum® SYBR® Green qPCR SuperMix-UDG
•
50 mM Magnesium Chloride (MgCl2) 1 ml 2 ×
•
ROX Reference Dye
•
20X Bovine Serum Albumin (ultrapure, non-acetylated) (1 mg/ml) 300 µl
Platiunum
•
SYBR Green qPCR Super Mix-UDG composition
•
SYBR Green I, 60 U/ml Platinum® Taq DNA polymerase, 40 mM Tris-HCl (pH
8.4), 100 mM KCl, 6 mM MgCl2, 400 µM dGTP, 400 µM dATP,
•
400 µM dCTP, 400 µM dUTP, 40 U/ml UDG, and stabilizers
Procedure
Standard reaction of 25 µl volume for Real-time Quantitative PCR of dsrAB
•
Platinum SYBR Green qPCR Super Mix-UDG
15.6 µl
•
Primer F (10µM)
1.0 µl
•
Primer R (10µM)
1.0 µl
•
Template
1 µl
•
PCR-H2O
ad 25 µl
Material and Methods
33
Table 14: Program for amplification of dsrB from standard and environmental samples
Cycle Repeats
1
2
3
4
1
1
45
81
Time Temperature
[min]
[°C]
Step
1 UDG Incubation
1 Initial denaturation
1 Cycle denaturation
2 Annealing
3 Elongation
1 Final Denaturation
2.0
2.0
0.5
0.5
0.5
0.5
50
95
95
54
72
55
PCR/Melt
DataAcquisition
Real time1
Melt curve2
1
Real time – Step at which real time data were collected
Melt curve – Step at which melt curve data were collected
2
Standards and environmental samples were analyzed in triplets (technical triplets)
to allow discrimination of technical measurement inaccuracies. To minimize pipetting
errors 3 µl of template were added to 75 µl of Mastermix (see standard reaction) and
subsequently separated into triplets using 23 µl of reaction volume. RT-qPCR was
performed according to the protocol (Table 14) The amplification associated fluorescence
was detected at each cycle during PCR. Melting curve analysis was performed after PCR
was completed over a temperature gradient of 0.5 °C temperature change starting from a
temperature of 55 °C to a final temperature of 95 °C.
Determination of dsrB copy number in Wadden Sea sediment
Data analysis and standard curve generation was performed using the Bio-RadiQ 5 software. Copy number was determined per µl of sample and manually converted
to copy number per µl of sediment. Average values of the copy numbers of the technical
triplets [copy number/µl of applied DNA] were multiplied with the dilution factor of 20
to obtain copy numbers per volume of extracted DNA [copy number/mg]. Therefore, a
DNA extraction efficiency of 100% was assumed. Some samples were not optimally
extracted and were excluded from final data set. Samples with non-optimal DNA
extraction were excluded from the final data set. Based on calculated volume of extracted
sediment [ml] values were finally converted to copy number per ml of sediment.
TM
(3)
0.5 ml
volume of extr. sediment [ml ] =
× weight of extr. sed iment [g ]
weight of 0.5 ml sed iment [g ]
copy number per ml sediment =
copy number per extr. sediment
× 1ml
volume of extr. sediment [ml ]
(4)
Data conversion and generation of the diagrams was performed using the Microsoft Excel
(Windows 2000) software.
Results
34
3. Results
3.1 Primer optimization
Prior to amplification of the dsrAB fragment from the environmental samples as
well as the cultivated SOP a temperature gradient-PCR was conducted to identify the
optimal annealing temperature of the newly developed primers. Primer pair 240 F and
dsrB 808 R were used for amplification of dsrAB from the environmental samples at an
annealing temperature of 55.5 °C (Figure 11).
environmental DNA JSS 0-3cm
Thiobacillus denitrificans
1200 bp
2000 bp
63.5
60.8
58.1
55.5
53.2
51.4
63.5
60.8
58.1
55.5
53.2
51.4
Figure 11: Temperature gradient PCR for amplification of dsrAB from Thiobacillus denitrificans and
environmental DNA JSS 0-3cm using the dsrA 240 F and dsrB 808 R primer. For primer optimisation
amplification was carried out at different annealing temperatures (T in °C).
2000 bp
NC
PC
1
2
3
4
5
Figure 12: Amplification of dsrAB from environmental DNA JSS 0-3cm. Amplification was carried out in
5 replicate reactions (number 1 – 5) at an annealing temperature of 55.5 °C. PC positive control, NC;
negative control
The dsrAB of Magnetospirillum gryphiswaldense was successfully amplified over
a wide range of annealing temperatures (Figure 13). Moreover, the 16S rRNA gene
fragment was amplified successfully from all reference strains (Figure 14 and Figure 15)
Thus the right DNA template concentration was provided. The DsrAB gene fragment
could not be amplified from the majority of the respective organisms at an annealing
temperature of 57.5°C (Figure ) A decrease in annealing temperature to 51.3° in the
following PCR also yielded no products (data not shown). Besides Magnetospirillum
gryphiswaldense successful amplification of the dsrAB fragment was achieved for
Thiobacillus denitrificans and Thiobacillus thioparus (data not shown), using the forward
primer dsrA 240F and an annealing temperature of 55.5°C.
2000 bp
51.3
52.2
53.2
54.3
55.4
56.5
57.3
58.4
NC
Figure 13: Temperature gradient PCR for amplification of dsrAB from Magnetospirillum gryphiswaldense
using the dsrA 240 F* and dsrB 808 R primer. For primer optimisation amplification was carried out at
different annealing temperatures (T in °C). The negative control was amplified at a temperature of 55.4 °C.
5
Results
35
1200 bp
2000 bp
DsrAB gene
16s rRNA gene
1200 bp
2000 bp
1
3
2
4
5
6
7
8
9
10
NC
12
Figure 14: Amplification of DsrAB gene fragment and 16s rRNA gene from cultured reference strains.
The 16s rRNA gene could be amplified from all reference strains except Rhodobacter capsulatus (3).
Amplification of the DsrAB gene fragment at an annealing temperature of 57.5 °C was only succesful for
Magnetospirillum gryphiswaldense (8). Number 1 – 10 Marichromatium gracile, Starkeya novella,
Rhodobacter capsulatus, Thiocapsa sp., Thiobacter truepii, Thiocococcus pfennigii, Thiomicrospira
pelophila, Magnetospirillum gryphiswaldense Leucothrix mucor, uncultured Beggiatoa; NC negative
control
1200 bp
2000 bp
1
2
3
4
5
6
7
8
9
10
NC
NC
Figure 15: Amplification 16s rRNA gene from cultured reference strains. The 16s rRNA gene could be
amplified from all reference strains except Thiomicrospira arctica (1). Thiomicrospira kuenenii (3) and
Starkeya novella (8). The 16s rRNA gene could be amplified from the former organisms in previous a PCR
(data not shown). Simultaneous amplification of the DsrAB gene fragment at an annealing temperature of
51.3 degree did no yield any visible band for any of the respective organisms. Number 1 – 10
Thiomicrospira arctica, Thiomicrospira psychrophila, Thiomicrospira kuenenii, Halothiobacillus
neapolitanus, Thiobacillus thioparus, Halochromatium salexigens, Rubrivivax gelatinosus, Starkeya
novella, Rhodobacter capsulatus, Magnetospirillum gryphiswaldense; NC negative control
Newly designed primers were tested for selectivity in order to detect whether the
covered sequence diversity was possibly less than the actual existing. The design of
primers on the basis of multiple sequence led to the integration of numerous
degeneracies. Due to up to 5 wobble positions per primer many variants are possible.
Hence, it was necessary to evaluate whether the primer variability was also reflected in
sequence diversity or whether certain sequences where amplified selectively.
Environmental clone dsrAB sequences were screened for the nucleotide composition at
the primer sites dsrA 240F and dsrB 808R to detect a possible selectivity of certain
primers variants. All respective nucleotides present at the wobble position of the primers
were represented in the sequences of the primer binding site. Moreover, the majority of
nucleotide combinations could be detected.
36
Results
3.2 Amplification of dsrAB from Wadden Sea sediment
A DNA fragment of ca. 2 kb size, encompassing most of the alpha and beta
subunit genes of the dissimilatory sulfite reductase, could be amplified from three distinct
sediment depths of the Wadden Sea, but not from the intertidal sulfur pond. Although
amplification of the 16s rRNA gene from the respective pond was successful for a
already low number of 20 cycles, dsrAB amplification with 35 cycles did not yield any
visible product in the electrophoresis gel. Partial and full sequences were obtained from
the respective PCR products. BLAST analysis of overall 182 dsrA partial sequences
against the GenBank database yielded 57 non DsrA sequences (Table 15), all of which
did not encode for similar protein functions as Dsr. This high fraction can be mainly
ascribed to direct sequencing without size specific prescreening of inserts. However,
some screening PCR products had the same size like the dsrA fragment and thus could
not be discriminated. Furthermore, the application of degenerated primers was reflected
in target sequence diversity, i.e. wobble position were highly divers. Hence a broad
spectrum of dsrAB sequences could be covered.
Table 15: Overview of non DsrAB BLAST hits of environmental sequences amplified with primer pair
dsrA 240 F and dsrB 808 R
Target sequence
cytochrome c5530 family protein [Methylococcus capsulatus str. Bath]
TPR repeat protein [Nitrosococcus oceani ATCC 19707]
TPR repeat [Nitrosomonas europaea ATCC 19718]
ccdB fusion protein [Cloning vector pTarg2]
Acyl-CoA dehydrogenase-like [Geobacter metallireducens GS-15]
hypothetical protein RferDRAFT_3095 [Rhodoferax ferrireducens DSM 15236]
hypothetical protein CA2559_10383 [Croceibacter atlanticus HTCC2559]
conserved hypothetical protein [Alkalilimnicola ehrlichei MLHE-1]
TonB-dependent receptor [Rhodospirillum rubrum ATCC 11170]
TonB-dependent copper receptor [Magnetococcus sp. MC-1]
LacOPZ-alpha peptide from pUC9 [Cryptosporidium hominis]
Predicted periplasmic solute-binding protein [Microbulbifer degradans 2-40]
conserved hypothetical protein [Parachlamydia sp. UWE25]
hypothetical protein AcidDRAFT_4187 [Solibacter usitatus Ellin6076]
ketol-acid reductoisomerase [Pelobacter carbinolicus DSM 2380
alpha-galactosidase [Bacillus clausii KSM-K16
sialic acid-specific 9-O-acetylesterase [Bacteroides thetaiotaomicron
GCN5-related N-acetyltransferase [Anabaena variabilis ATCC 29413]
putative integrase [Nitrosococcus oceani ATCC 19707]
O-acetylserine sulfhydrylase [Methanosarcina thermophila]
hypothetical protein PatlDRAFT_2122 [Pseudoalteromonas atlantica
urf53 [Chlorella vulgaris]
conserved hypothetical protein [Xanthomonas campestris pv. campestris str. 8004]
putative succinate dehydrogenase flavoprotein subunit [Nocardia farcinica]
probable serine/threonine protein kinase [Rhodopirellula baltica SH 1]
hypothetical protein BURPS1710b_A1869 [Burkholderia pseudomallei 1710b]
hypothetical protein SfumDRAFT_2942 [Syntrophobacter fumaroxidans MPOB]
IstB-like ATP-binding protein [Solibacter usitatus Ellin6076]
hypothetical protein BbacK_01001192 [Bartonella bacilliformis KC583]
Oxidoreductase, N-terminal [Solibacter usitatus Ellin6076]
TPR repeat [Chlorobium phaeobacteroides DSM 266]
antigen, putative [Treponema denticola ATCC 35405]
ribonuclease G [Xanthomonas axonopodis pv. citri str. 306]
Acyl-CoA dehydrogenase-like protein [Burkholderia sp. 383]
Histidine kinase A, N-terminal [Methanococcoides burtonii DSM 6242]
CoA-binding [Chloroflexus aurantiacus J-10-fl]
PREDICTED: hypothetical protein XP_896092 [Mus musculus]
hypothetical protein Npun02005532 [Nostoc punctiforme PCC 73102]
E value
3e-33
1e-86
9e-47
2e-33
e-17
1e-12
8e-36
6e-15
2e-42
6e-36
3e-12
4e-39
1e-55
7e-61
4e-125
7e-111
1e-36
4e-19
6e-79
1e-88
6e-15
7e-24
1e-24
2e-46
2e-52
0.055
0.48
2e-16
0.010
1e-40
3e-12
1e-17
4e-69
0.011
1e-23
2e-57
0.86
7e-28
Identity
107/278 (38%)
184/319 (57%)
106/223 (47%)
78/116 (67%)
52/158 (32%)
67/222 (30%)
70/108 (64%)
40/67 (59%)
121/317 (38%)
105/277 (37%)
38/40 (95%)
83/133 (62%)
110/252 (43%)
133/321 (41%)
221/269 (82%)
190/308 (61%)
105/253 (41%)
62/173 (35%)
149/226 (65%)
180/288 (62%)
55/147 (37%)
56/123 (45%)
73/199 (36%)
111/263 (42%)
21/286 (42%)
45/164 (27%)
70/259 (27%)
2/75 (56%)
18/51 (35%)
88/226 (38%)
41/120 (34%)
43/99 (43%)
135/249 (54%)
40/97 (41%)
53/107 (49%)
126/277 (45%)
20/68 (29%)
69/179 (38%)
Accession
AAU90385
ABA57571
CAD85209
CAD27778
ABB32438
EAO41905
EAP86435
EAP35670
ABC24246.1
EAN28739.1
EAL34605.1
ZP_00314664.1
CAF23477.1
EAM58117.1
ABA89756
BAD62954.1
AAO79196.1
ABA23672.1
ABA57146.1
AAG01804.1
EAO67128.1
BAA99566.1
AAY48221.1
BAD55792.1
CAD77107.1
ABA51353.1
EAO22105.1
EAM60184.1
ZP_00947034.1
EAM54819.1
EAM34521.1
AAS13216.1
AAM37615.1
ABB05792.1
EAM99616.1
EAO57496.1
XP_901185.1
ZP_00106364.2
Results
37
3.3 Phylogeny and diversity of the DsrAB from SOP
Full sequences were loaded into the ARB database and comparative sequence
analysis was performed based on a combined DsrAB dataset, comprising a maximum
number of 532 comparable amino acid positions. Overall 102 DsrAB full length
sequences from the environmental clone library (64 different sequences), cultivated (14
sequences) and uncultivated (4 sequences) bacteria as well as selected, cultured
representatives of SRP (20 sequences (Klein et al., 2001)) and additionally 3 partial
length sequences from uncultivated SOP served for reconstruction of the phylogeny of
the DsrAB.
3.3.1 Overall DsrAB diversity
DsrAB Consensus tree
For reconstruction of DsrAB phylogeny different inference methods were applied,
including Neighbor Joining, Maximum Parsimony, Maximum Likelihood and Tree
Puzzle, to achieve optimal resolution of evolutionary relationships. The obtained tree
topologies showed no major discrepancies Differences occurred in the phylogenetic
positioning of Thiobacillus denitrificans as well as the Wadden Sea clone clusters C, B
and E. Based on this results a consensus tree (Figure 16) was constructed representing the
phylogenetic relatedness as best as possible.
The DsrAB clone library contained 64 distinct sequences and was dominated by
clones which affiliated with Gammaproteobacteria (46 different sequences) constituting
72%. Eighteen of 64 clones (28%) grouped together in the monophyletic cluster E with
probable alphaproteobacterial affiliation. The majority of gammaproteobacterial
associated sequences (57%) grouped together in the monophyletic Cluster A (16 different
sequences) and B (10 different sequences) with a sequence identity between 81.9 and
99.8% and 92.3 and 99.8% respectively. Clones JSS d1_F3 and JSS d3_G8 were related
to an uncultured bacterium from the German Wadden Sea (Fosmid 194A5, Janssand)
displaying 95.5 and 94.2% identity to the fosmid sequence. Clone JSS d1_F6 clustered
with the sulfur-oxidizing gamma 1 symbiont of Olavius algarvensis (90.9% sequence
identity). Clones JSS d1_A2, JSS 2_B5 and d2_A9 JSS (96.6 – 98.4% sequence identity)
of the monophyletic cluster C were more distantly related to that sequence (84.6 – 85.4%
seqeunce identity). A fourth group of sequences formed the monophyletic Cluster D (92.5
–99.8% sequence identity) and affiliated with the sulfur oxidizing gamma 2 symbiont of
Olavius algarvensis (90.2 – 96.2% sequence identity). In addition, clone JSS d2_dB1
showed a close relationship to the gamma 2 symbiont (96.0% sequence identity). Clone
JSS d2_D7 took an intermediary position displaying 82.6% sequence identity to the
gamma 2 symbiont of Olavius algarvensis and 82.6% sequence identity to uncultured
bacterium from the Mediterranean Sea (BAC13K9). The phototrophic sulfur oxidizer
Allochromatium vinosum was found to be the closest cultured relative of the majority of
the Gammaproteobacteria-related clones. Only one sequence (JSS d2_D9) grouped
closer with Alkalilimnicola ehrlichei (73.3% sequence identity).
38
Results
The monophyletic Cluster E comprised nearly one third of all clones, including 18
distinct sequences (79.3 – 99.7% sequence identity) and hence constituted the second
largest fraction in the DsrAB clone library, accounting for 28% of all clones. Maximum
Likelihood and Neighbour Joining methods indicated an alphaproteobacterial affiliation
of this cluster. However, Maximum Parsimony yielded contradictory results. Moreover,
sequence identity between the clones and Alkalilimnicola ehrlichei was on average 1.3
higher compared to the Magnetospirillum sequences (71.2% compared to 72.5%).
Interestingly, neither of the clones affiliated with the sulfur oxidizing
Betaproteobacterium Thiobacillus denitrificans or the phototrophic chlorobi . The latter
phylum branched deepest within the SOP. Moreover, the DsrAB of two Chlorobium
phaeobacteroides strains, DSM266 and BS1, was represented within the phylogenetic
tree and branched of separately within the chlorobial cluster.
Interestingly, the closest relatives among the sulfate-reducing prokaryotes were
members of the deep-branching SRP phyla, i.e. Euryarchaeota and Firmicutes.
Archaeoglobus fulgidus affiliated closest displaying 49.4% sequence identity to
Chlorobium phaeobacteroides. Furthermore, the DsrAB sequence of the uncultered
sulfate reducing bacterium (Wadden Sea Fosmid 39f7) showed a comparable high
percentage of sequence identity (49.7%). However this does not necessarily result from
close evolutionary relatedness as reflected in the tree topology.
Figure 16: DsrAB concensus tree for the known SOP and SRP as well as environmental clone sequences
based on Neighbor Joining, Maximum Parsimony, Maximum Likelihood and Tree Puzzle phylogenetic
inference methods. Multifurcations were introduced when branching orders were not supported by all
phylogenetic inference methods used. Environmental clone sequences are colour coded according to the
depth they originated from. Phylogenetic sequence clusters are indicated by brackets.
Figure 17: DsrAB tree for the known SOP and SRP as well as environmental clone sequences as calculated
with MrBayes phylogenetic inference method. Environmental clone sequences are color coded according to
the depth they originated from. Phylogenetic sequence clusters are indicated by brackets.
(see following pages)
Results
39
40
Results
Results
41
DsrAB MrBayes tree
We applied the MrBayes inference method for phylogeny of DsrAB (Figure 17)
to check our results from the consensus tree. Overall, both trees displayed congruent tree
topologies. Slight differences occurred in the phylogenetic positioning of the sulfuroxidizing gamma 1 symbiont of Olavius algarvensis branching of between
Allochromatium vinosum (92.02% sequence identity) and the uncultured bacterium
(fosmid 194A5) from the German Wadden Sea (91.73% sequence identity). The Wadden
Sea clone Cluster B showed monophyletic affiliation with the Allochromatium vinosum
associated sequences, as before revealed by Maximum Likelihood analysis. Consistently
with results from Maximum Parsimony analyis clone d2_D7 affiliated with cluster C
(80.7 – 98.4% sequence identity), which likewise branched of differently. Furthermore,
the inferred phylogeny strongly supported the affiliation of the cluster D with the sulfur
oxidizing gamma 2 symbiont of Olavius algarvensis. Additionally, a clustering of
sequences into surface and deeper sediment layers associated subclusters became
evidently. Moreover, the topology of the MrBayes tree strongly supported the relatedness
of clone cluster E to the alphaproteobacterial Magnetospirillum spp. sequences.
Accordant to the consensus tree, SRP of the Euryarchaeota and Firmicutes affiliated
closest with the chlorobi.
3.3.2 Comparative analysis of DsrAB consensus tree and 16s
rRNA gene tree
The comparison of the tree topologies of DsrAB and 16s rRNA gene (Figure 18)
revealed one putative major LTG event within the SOP. The chlorobial affiliation of the
alphaproteobacterial Magnetococcus sp. points at horizontal transmission of the dsrAB.
Analysis of the insertion and deletion patterns yielded additional evidence.
Figure 18: Comparison of the DsrAB consensus and 16S rRNA trees. Known SOP are colour coded
according to their taxonomic classification. (see next page)
42
Results
Results
43
3.3.3 Comparative analysis of insertions and deletion patterns
As additional indications of the deduced phylogenetic relationships, particularly
with respect to LGT, insertion and deletion patterns within the DsrAB amino acid
sequences of SOP and SRP were examined (Figure 19).
Within the DsrAB sequences of the analyzed SOP Thiobacillus denitrificans
harbored three unique insertions within the alpha subunit at alignment positions 33 to 37,
68 and 252 to 267. The first insertion was also found in the DsrAB sequence of
Alkalilimnicola ehrlichei and Wadden Sea clone JSS d2_D7. One unique insertion at
alignment position 211 to 212 (DsrA) was found for Magnetococcus sp.. In addition, the
latter shared a unique deletion in the alpha subunit present in all members of the chlorobi
at alignment position 268, supporting the suggested LGT event. In contrast, the missing
amino acid was highly conserved in all other DsrAB sequences of SOP and SRP.
Interestingly, all SOP shared an insertion at alignment position 332 to 333 with
Desulforhopalus vacuolatus, Desulfobacter latus and Desulfitobacterium hafniense.
SOP associated deletions at alignment positions 211 to 226 (DsrA) and 620 to 626
(DsrB) are also present within the DsrAB sequences of Desulfotomaculum alkaliphilum,
Thermodesulfovibrio islandicus, Desulfitobacterium hafniense, Archaeoglobs fulgidus
and fosmid clone 39F2, with Desulfoarculus baarsii and Desulfobacca acetoxidans
displaying the deletion within the beta subunit as well. Overall, one deletion at alignment
positions 702 to 703 (DsrB) was found to be characteristic for the DsrAB sequences of
all analyzed SOP, in contrast, all analyzed SRP possessed a unique insertion at this
position.
D. vacuolatus
D. baarsii
D. latus
D. acetoxidans
D. variabilis
D. vulgaris
D. kuznetsovii
D. tiedjei
D. alakiliphilum
T. norvegica
T. commune
T. islandicus
D. hafniense
Fosmid clone 39F7
A. fulgidus
T. denitrificans
A. ehrlicheii
A. vinosum
Magnetococcus sp.
C. tepidum
C. limicola
M. magnetotacticum
M. magneticum
M. gryphiswaldense
30
40
....|....|....|......
PKRFPGV-----EHFHTMRAN
PNLFPGV-----EHFHTVRVN
PNLFPGV-----EHFHTMRVN
PNEFPGL-----EHFHTMRIN
PETFPGV-----AHFHTMRIN
PEKFPGV-----AHFHTVRVA
PEKFPGI-----AHFHTIRVN
PQQFPGV-----AHFHTVRVN
PEQFPNV-----AEFHTVRLN
PDLFPNV-----AHFHTMRVN
QERFPAI-----WAFHTIRVH
ADQFPGV-----AHFHTVRVN
PEDYPNL-----AAFHTVRIN
PDKWPGI-----AHFHTMRIN
GEQIPEV-----EHFHTMRIN
KDENHQPLFPEAAEFHTLRIQ
KDENGKAKFPDAAEFHTMRVM
GKVFPSSK-----EFHTVRVQ
PELFPEAK-----EFHTMRIQ
KDKFPEAA-----EFHTMRIQ
KDKIPEAA-----EFHTMRIQ
AAQFPESA-----EFHTLRIM
AAQFPESS-----EFHTLRIM
ADKHPESS-----EFHTLRIM
200
210
220
230
240
.|....|....|....|....|....|....|....|....|..
QAAVKEYIAGNYPGNGGAHSGRDWG---KFDIQKEVIDLCPTEC
QAAVAAYAGEXIPPRAAAHGARGV----KLDVQKEVXKLCPTGC
QDAVNKYVEQDPAYPANAGAHKGKDWG-PFDIEKEVINLCPTGC
QSAVKAYAGGELKPRGGGSPYA------KLDIQADVVELCPTGC
QEAVASYVGGEFAPNGGAHSGRNWG---AFDIQKEVIDLCPTRC
AEAVKAYVAGEFKPNAGAHSGRDWG---KFDIEAEVVNRCPSKC
QEAVKAYIAGDIVPRGGSHKGRDWG---KFDIQKEVINLCPSQC
QEAVKAYVGGELKPNAGAHSSRDWG---KFDIIAEVINLCPTNC
QEAVREYVNNGFD------------------IQELVIDRCPTGC
QEAVRAYVAGELVPNGGAHGFEKRP----LDIQAEVIDKCPTRC
QEAVRAYVRGELKPNAGAFSDRDWGP---FDIKKEVIDLCPTKC
QDKVKDYAKKGMN------------------IQEEIVNFCPGKC
QEAVKAYVAEGLN------------------IQAVVCDRCPTKC
KAAVTEY--AKTG----------------MDIKADVVDNCPAKC
QEAVKEYASWMDI-------------------ENEVVKLCPTGA
EEQVQAYYKAHGM----------------NDLVNDVISKCPTKA
QEEVKAFVARKGR----------------KYVVDNVISRCPTNA
QAEVKHYIADKGR----------------QYYIDNVITRCPTKA
EAEVAKFIETRGVAY--------------LTNNVITRCPTNAIV
HDEVKAWIAEKGV----------------DALVNNVINRCPTKA
RDEVTAWVEKKGM----------------DALVNQVINLCPTRA
QKEVQAKVKQLGR----------------KEFINQVIRMCPTQA
QDEVKAKVKQLGR----------------KEFINQVIRMCPTQA
QDEVKAKVAQLGR----------------KEFINQVIRMCPSQA
D. vacuolatus
D. baarsii
D. latus
D. acetoxidans
D. variabilis
D. vulgaris
D. kuznetsovii
D. tiedjei
D. alakiliphilum
T. norvegica
T. commune
T. islandicus
D. hafniense
Fosmid clone 39F7
A. fulgidus
T. denitrificans
A. ehrlicheii
A. vinosum
Magnetococcus sp.
C. tepidum
C. limicola
M. magnetotacticum
M. magneticum
M. gryphiswaldense
250
260
270
..|....|....|....|....|....|....|...
MWMEDGE-LK----------------I-DNSECTRC
MSYEGGK-LA----------------I-NNRECTRC
MKFENKE-LT----------------I-NDEHCTRC
MSFDGS-TLK----------------I-DNRNCNHC
MKYEGGK-LA----------------I-NTKECTRC
MKWDGSK-LS----------------I-DNKECVRC
MWMEDGK-LV----------------I-DNKECVRC
MCWDGSK-LE----------------I-NNAECVRC
LEWDGNE-LK----------------L-DAPECVRC
MEWDGKN-LK----------------I-WDEDCVRC
MYWDENEQKLX---------------I-NNSECNRC
ITWDGNN-LT----------------I-NNSDCLHC
LKFDAETQELS---------------V-IAEECTRC
MDWDGKE-IH----------------I-ENSYCRKC
IKWDGKE-LT----------------I-DNRECVRC
ITLVASDKFAPSETVSAANLGDGNTLCIDNKNCVRC
LSLNDDDTLD----------------V-DNGSCVRC
LSLNDDDTLD----------------V-NNRDCVRC
MNEDGKSITI------------------NNSDCVRC
IRLQDGDIDI------------------STRDCVRC
ISLKEGEMHI------------------ETRDCVRC
LGLNDDDTLD----------------V-DNKSCVRC
LGLNDDDTLD----------------V-DNKSCVRC
LALNDDDTLD----------------I-DNKSCVRC
330
340
....|....|....|....|....
FMKIEKENDYQ-ELIDLIEKIWDW
FIKMEEPYDEF---KEFVANTWDW
FVEVNADNDYK-EITEIIENIWDW
FMEMEPPYDNL---KEIAGNIMEW
FVKVEEPYDGI---KEVIESIWNW
FVAAEEPFDEI---KEVVEKIWDW
FMKAEPPYDNV---KEVIEKIWEF
FMKMEPPFDEL---KELIEKVWDW
FMKMEAPYQEV---KDLLEKIWEW
FMKMEPPYEEF---HEFVEKMWDW
FMKVEPPYDEL---KDFIYKAWDY
FMGMKPPYQEI---KDLIKNMWDW
FMKMEVEDDFQ-EFKDMIERIWEW
FYKMEPPYTEL---KDIMERIWDL
FVEVEKPYDEI---KEILEAIWDW
FMKLESEEDFE-KLNDLARNIIDF
FLPLNTEEDYD-RVVEIAEEIIDF
FKKLETEEDYE-SLVELAETIIDF
FMKMEDEDDIE-AFIDMIDSMIDY
FMKMETDEDRE-AFIELIEEIIDW
FMKMESDEDVE-TFIEQIETIIEW
FMKLETEEDFE-KLVELAHNMLDF
FMKLETEEDFA-GLVELAHNMLDF
FMKLETDEDYE-RLLEMAHNMLDF
Results
44
D. vacuolatus
D. baarsii
D. latus
D. acetoxidans
D. variabilis
D. vulgaris
D. kuznetsovii
D. tiedjei
D. alakiliphilum
T. norvegica
T. commune
T. islandicus
D. hafniense
Fosmid clone 39F7
A. fulgidus
T. denitrificans
A. ehrlicheii
A. vinosum
Magnetococcus sp.
C. tepidum
C. limicola
M. magnetotacticum
M. magneticum
M. gryphiswaldense
620
630
640
....|..|....|....|....|....|....|....
LLDTLASYGN-----SYPVGGT------GAGVTNIVH
LKKPCRAG-------NMPIGGT------GHSVTNIVH
LKNDLASRKFAGGSQKFPIGGT------GAGVTNIVH
LKAELAQR-------GIPKGGT------GHSITNIVH
LIKDLESRKFDGGSFKFPIGGT------GAGVTNIIH
LKEDLASRKFDGGSLKFPIGGT------GAGVSNIVH
LKKDLWSRKFVSGSYKFPIGGT------GASITNIVH
LIDDLKGRKFGAGSYKFPIGGT------GSGITNIVH
LIAEVKAK-------GLPVGGT------GNSITNMVH
LIDDLKSRGN-----WFPIGGT------GASVTNIVH
LIEDLKNRKHPNGSYKFSIGGT------GAGISNIVH
LLDELKAK-------KYMIGGI------GSRVSNIVH
LIAALGAK-------GLPVGGT------GASVSNIVH
LIKALKDI-------GHPVGGI------GNAISSIVH
LINEVQERV------GFPCGGTWDAVKGEYGLSNIVH
LIKALSDK-------GYPIGGT------ANSVSMIAH
LIDKLRAE-------GFPVGGT------GPSVAMISH
LIDALEGA-------GFVVGGT------QNSVAMISH
LIKALEAA-------GHPVGGT------GNSITSMAH
MIAELESL-------GFPVGGT------GMCVSAVSH
MIQELESL-------GFPVGGT------GMCVSSVSH
LVKRLEAE-------GFPVGGT------GNSVGFISH
LIERLEKD-------GFPVGGT------GNSVGFISH
LIDRLEAA-------GFPVGGT------GNSVGFISH
700
710
..|....|....|..
MCGAVHASDTAILGI
MCGAVHCSDIAILGV
MCGAVHCSDIAILGY
MCGAVHCSDIAILGV
MCGAVHCSDIAILGY
MCGAVHCSDIGVVGI
MCGAVHCSDIALLGI
MCGAVHCSDIAILGI
MCGAVHCSDIAILGI
MCGAVHCSDIAILGV
MCGAVHCSDLALVGI
MCGAVHCSDIAFVGI
MCGAAHCSDIAIVGV
MCGAVHCSDIAILGI
MCGAVHASDIAIVGI
NCGGQ--ADIAIIVQ
NCGGQ--GDIAINIQ
NCGGQ--GDIAINIQ
NCGGQ--ADIAIVVQ
NCGGQ--ADIAVVVK
NCGGQ--ADIAIVVK
NCGGQ--GDIAINVQ
NCGGQ--GDIAINVQ
NCGGQ--GDIAINVQ
Figure 19: Insertion and deletions patterns of DsrAB of selected SOP and SRP
3.3.4 Vertical DsrAB diversity
Clone libraries of dsrAB from different sediment layers were established to
elucidate vertical distribution of SOP. Overall 28 different full sequences were obtained
from the surface sample (JSS 0 – 3 cm), 21 (JSS 8 – 19cm) and 15 (JSS > 19cm) different
full sequences, respectively, from the deeper sediment layers.
Clones from all depths seemed to be randomly distributed over the both major
clone groups, the gammaproteobacterial branch as well as the Wadden Sea clone cluster
E. Nevertheless, a clustering of sequences of different depths could be observed. Within
cluster A clones JSS d2_B2 – d3_B3 (93.43 –99.8% sequence identity) of the deeper
sediment layers and clones JSS d1_G4 – d1_H5 (83.7 –99.2% sequence identity) of the
surface layer formed two distinct subclusters. A similar pattern could be found in Cluster
D were clones JSS d3_B4 – d3_F9 (89.9 – 99.7% sequence identity) and clones JSS
d1_B1 – d1_D4 (95.3 – 99.8% sequences identity) clustered together to two distinct
subgroups from the depth and surface sample, respectively. Sequences of cluster B
originated almost exclusively from the surface sediment (92.3 – 99.8% sequence
identity). As already mentioned above, clones from all sediment depths were also present
within the cluster E, in which again clones JSS d1_D1 – JSS d1_E6 (88.9 –89.0%
sequence identity) from the upper sediment layer clustered together in a distinct
subgroup.
Some clones within the clusters A, B, D and E displayed very high amino acid
sequence identity, e.g. d2_D2 and D2_G9 (99.8%), d1_C6, d1_A1 and d1_A3 (99.8%),
d1_C2 and d3_F12 (99.7%) as well as d3_E11 and d2_C6 (94.7%), which can also be
found on the nucleotide basis. These findings point at the the existence of identical clones
3.3.5 Rarefaction analysis
To estimate overall diversity of SOP in the Wadden Sea sediment a DsrAB clone
library was established from three distinct sediment depths and a total number of 112
were screened. In order to proof whether the number of retrieved sequences efficiently
covered the sequence diversity in our clone library we performed a rarefaction analysis
(Heck et al., 1975) on a combined dsrA dataset including sequences from all depths
(Figure 20). Based on sequence comparison of the dsrA fragment (~ 600 bp) 46 distinct
Results
45
sequence types were identified among 112 clone sequences. DsrAB diversity was higher
than expected as represented by the continous slope. However the rarefaction curve
slightly turned into the saturation limit. Thus we could assume coverage of the major
dsrAB sequence types by numer of anaylzed clones but further screening of the clone
library was predicted to yield additional novel dsrAB sequence types.
50
Number of sequence types
45
40
35
30
25
20
15
10
5
0
0
10
20
30
40
50
60
70
80
90
100
110
120
Number of screened clones
Figure 20: Rarefaction curve for the different dsrA sequence types. The number of sequence types is plotted versus the
number of dsrA clones. The dotted lines represent 95% confidence intervals.
3.3.6 Comparison 16S rRNA gene identity and DsrAB
identity
To evaluate the degree of sequence conservation of DrsAB and to compare
sequences identity on the basis of the functional gene with the 16S rRNA gene, identity
matrices (see appendix Table A1 and A2) were calculated for the respective nucleotide
and aminoacid sequences of representative organisms using the ARB Similarity Matrix
calculation editor. Comparative sequence analysis reavaled a high degree of sequence
conservation of the the DsrAB, with an average sequence identity of 64. 01% the
organisms analyzed. Absolute values varied between 54.2% (Thiobacillus denitrificans Chlorobium tepidum and Magnetococcus sp., respectively) and 96.41%
(Magnetospirillum magnetotacticum - Magnetospirillum magneticum). Overall, the
degree of sequence conservation of the 16s rRNA gene (82.0% average sequence
identity) exceeded the one of DsrAB. In general, aminoacid sequences were highly
conserved for species closely related on 16S rRNA gene level. For example, members of
the chlorobi, displaying a 16S rRNA gene sequence identity between 91.8% and 99. 2%
(average value 96.4%), showed comparable high values for the DsrAB, i.e. between
73.7% and 90.5% sequence identity (average value 83.9%). Similarly, the
alphaproteobacterial Magnetospirillum sp. sequences (95.8 – 99.4% 16 s rRNA gene
sequence identity, average value 97.1%) exhibited a considerable high DsrAB sequence
identity (87.5 – 96.4%, average value 90.9% ). Highest conservation of DsrAB was found
46
Results
between the closely related species Magnetospirillum magnetotacticum and
Magnetospirillum magneticum (99.4% 16s rRNA gene sequences identity) displaying
96.4% aminoacid and 92.8% nucleotide sequence identity. Although the dsrAB
nucleotide sequence in general was less conserved, displaying 41.92% average sequence
identity, a high degree of sequence conservation on the aminoacid level could also be
detected on nucleotide basis (53.4 – 81.1% sequence identity for the chlorobi, average
value 70.5%; 84.2 – 92.8% sequence identity for the Magnetospirillum sp. sequences,
average value 87.7%). Interestingly, the DsrAB of the two Chlorobium phaeobacteroides
species (95% 16 s rRNA gene sequence identity) was less conserved, displaying 74.8%
aminoacid sequence identity. However, in comparison, the nucleotide sequence identity
of 73.2% was high. Moreover the DsrAB of the Chlorobium sp. sequences (95.6%
average sequence identity 16s rRNA gene to Chlorobium phaeobacteroides) showed a
much higher sequence identity to Chlorobium phaeobacteroides DSM 266 (88.7%
average) than to Chlorobium phaeobacteroides BS 1 (74.9%). Interestingly, within the
gammaproteobacterial cluster Allochromatium vinosum displayed a lower DsrAB
sequence identity to Alkalilimnicola ehrlichei (76.1%) than to the sulfur oxidizing gamma
2 endosymbiont of Olavius algarvensis (80.6%). In contrast, the respective 16 s rRNA
gene sequence identity was higher (90.1% compared to 86.7%).
Surprisingly, the Alphaproteobacterium Magnetococcus sp. displayed a
significantly lower degree of sequence conservation on the DsrAB level (58.1% average
sequence identity) to the alfaproteobacterial sequences than to the chlorobial sequences
(65.2% average sequence identity).
3.4 Amplification of dsrAB from cultivated SOP
In addition to amplification from environmental DNA dsrAB should also be
amplified from cultivated SOP to better resolve phylogenetc relationsships. Temperature
gradient PCR yielded a dsrAB product over a wide range of annealing temperatures from
Magnetospirillum gryphiswaldense. Subsequently, primer pair 240F* and dsrB 808R
were used for amplification of dsrAB from cultivated SOP at an annealing temperature of
57.5 (Figure 14) and 51.3°C (Figure 15), respectively. Nevertheless, amplification of
dsrAB from sulfur-oxidizing reference strains was not successful in the majority of cases.
Finally, amplification of the dsrAB could only be achieved for Thiobacillus thioparus,
displaying a weak band of approximately 2000bp in the agarose gel (data not shown).
Repetition of the PCR is necessary to approve this result.
3.5 Primer selectivity test
Primer dsrAB 240F was tested for selectivity in order to detect wether the covered
sequence diversity is possibly less than the actual existing and in order to reveal a
selective amplification of dsrAB from environmental samples. Amplification of the dsrAB
from the fosmid failed with both primer combinations. Besides the positive control no
visible PCR product of dsrAB from fosmid clone 194A5 was detected in the gel for
neither of the forward primers, dsrA 240F and dsrA 240F*. To test wether primer 240F is
capable of amplifying the dsrAB fosmid target sequence and to reveal a possible
selectivity in the amplification of dsrAB from the environmental samples a repetition of
the PCR is necessary.
Results
47
3.6 Real-time Quantitative PCR
Primer specificity
In order to ensure a specific amplification of the target gene we employed gel
elctrophoresis, sequencing and melting curve analysis. Gel elctrophoresis exclusively
yielded a visible band for the respective 400bp sized dsrB fragment. Furthermore a
similarity based BLAST search (Altschul et al., 1990) of the obtained clone seqeuences
against the Genbank database (www.ncbi.nlm.nih.gov/blast) did always hit the dsrB
subunit as closest relative for all analyzed sequnces. In addition, results from melting
curve analysis (Figure 21) implied a contamination free amplification of the dsrB target
fragment. All standards as well as all samples showed an identical melting melting
temperature (Tm) of approximately 86.5°C. Neither primer dimers, nor multiple products
could be detected.
Figure 21: Melting curve for the dsrB fragment amplified from the standard (blue) and the environmental
DNA (red). The first derivative of the fluorescence plotted against temperature (-Δ(RFU)/ ΔT) is plotted
against temperature resulting in a peak at the specific melting temperature (Tm).
An approximately 600bp-sized DNA fragment amplified from Wadden Sea clone
E2 containing the respective 400bp sized dsrB fragment served as template for the
generation of an appropriate standard curve. The applied dilution series sufficiently
covered the range of the environmental samples. Triplets of the high diluted standards
showed significant variations in the treshold cycle (Figure 22).
48
Results
Figure 22: Standard curve for amplification of dsrB from the standard with primers dsrA 403F and dsrB
808R for concentration from 1.7 ×107 copies × µl-1 to 1.66 ×104 copies × µl-1
Data analysis revealed remarkable variation of the copy numbers of dsrB per ml
of extracted DNA for the technical as well as the biological triplets of the respective
depth intervals (Figure 23). Final calculated copy numbers of the dsrB fragment ranged
from 106 to 107 copies per ml sediment over the vertical profile and displayed relevant
variations between each sediment depth (Figure ) Standard deviations were considerable
high, displaying an order of magnitudes of 106 copies per ml sediment. Highest copy
numbers of 1.52 to 1.82×107 were reached in the zone of the oxic-anoxic interface in 3 to
4 cm depth and in the anoxic zone of the deepest sediment layer in 16.5 to 17cm depth
(1.85×107 copies). Close to the surface of the sediment significantly lower copy numbers
of 1.11 to 1.13×107 numbers were detected.
Results
49
copy number per µl extracted DNA
5,00E+04
5×104
1,00E+05
1×105
1,50E+05
1.5×105
2,00E+05
2×105
2,50E+05
2.5×105
3.5. - 4.0
16.5 - 17.0
depth interval [cm]
1.0 - 1.5
0,00E+00
0
Figure 23: Copy number of dsrB per µl of extraced DNA for the three replicates of a certain sediment
depth (Note, the standard deviation for the technical triplets as well as the considerable variation of copy
numbers of the biological triplets for the same sediment depth.)
copy number per ml sediment
7 1,50E+0
7 2,00E+0
7 2,50E+0
7 3,00E
0,00E+0
5,00E+06 1,00E+0
0
5×10
1×10 1.5×10
2×10
2.5×10
0
6
7
7
7
7
7
-2
0
2
depth [cm]
4
6
8
10
12
14
16
18
Figure 24: Copy number of the dsrB per ml sediment over the vertical sediment profile
Literature
50
4. Discussion
Aim of the present study was to reveal the diversity of SOP in Wadden Sea
sediments based on the establishment of a phylogeny of a functional gene essential for
sulfur oxidation, namely the DsrAB. Due to the polyphyly of SOP, this gene of the
“reverse” Dsr sulfur oxidation pathway offers an adequate tool for environmental analysis
of the SOP. Moreover, the abundance of the respective organisms should be assessed by a
newly developed dsrB targeted quantitative real- time PCR approach. The results of the
completed work will be discussed as follows:
At first methods and primer issues will be evaluated. Afterwards the DsrAB and
16S rRNA phylogeny of cultured SOP will be introduced and compared with regard to
putative events of lateral gene transfer. This is followed by a discussion of the main focus
of this study, the diversity of environmental SOP in the Janssand sediment. Subsequently
the SybrGreen based quantitative real-time PCR assay for the assessment of the dsrB
copy number is debated. Finally, a synopsis in combination with an outlook will provide
a critical view of the whole dsrAB approach and moreover propose molecular
experiments for future ecological studies.
4.1 Primer evaluation and methodological aspects
4.1.1 Evaluation of dsrAB primers for diversity studies and
qRT-PCR
In the course of the study several primers were tested to optimize amplification of
dsrAB of sulfur oxidizing bacteria. Previously, only primers were available to specifically
amplify dsrAB from SRP (Wagner et al., 1998; Klein et al., 2001) but not from SOP.
Hence, there was a clear lack to apply this functional marker molecule for diversity
studies of SOP employing the “reverse” Dsr pathway and to elucidate the evolution of
SOP and SRP. The primers dsrA 240F and dsrB 808R were chosen for the DsrAB
diversity survey since best results in amplification of environmental dsrAB were
obtained. Initially conducted test PCRs yielded a lower product with the primer
combination dsrA 116F and dsrB 808R (pers. communication M. Mußmann 2005). The
primer 240F was improved after the additional dsrAB sequence of the fosmid clone
194A5 became available. The phenomenon of numerous “non-dsrAB” amplificates is
probably associated with the use of degenerated primers and has already been reported by
Klein et al. (2001) for the dsrAB of SRP. The latter is homolog to the covered dsrAB
fragment of SOP in this study. The used primers target all published sequences of
cultured and uncultured organisms. Just a very few sequences are not fully matched.
Primer dsrB 403F was designed for quantitative real-time PCR to allow the amplification
of a rather short dsr fragment. The screening of the nucleotide composition of the
environmental clone dsrAB and dsrB sequences at the primer sites dsrA 240F, dsrB 403F
and dsrB 808R indicated no selectivity of certain primer variants. Virtually most primer
sequence variants are represented in the clone libraries indicating little bias in primer
binding. Analysis of the environmental dsrAB sequences additionally revealed that the
employed primers are potentially able to cover a wide range of sequences.
Literature
51
However, unexpectedly dsrAB of some strains such as Thiocapsa marina and
Marichromatium gracile was not amplified , although the latter organism and Thiocapsa
roseopersicina were shown to contain Dsr genes by Southern hybridization (Dahl et al.,
1999). Furthermore, dsrAB from the fosmid clone 194A5 was not targeted at all by
completely matching primers. Unfortunately the selectivity test of Primer dsrA 240F
failed. The lack of amplification of dsrAB from the fosmid 194A5 by both primers, dsrA
240F and 240F*, does not allow for a conclusion on the selectivity of primer dsrA 240F.
Primer dsrA 240F* should have amplified the respective gene fragment. The lack of
product could be due to suboptimal PCR conditions, e.g. template concentration or
annealing temperature. Overall, these results indicate that primers and PCR conditions
need to be optimized further in order to cover a maximal diversity of SOP.
4.1.2 Amplification of DsrAB gene from cultivated SOP
For comprehensive reconstruction of the DsrAB phylogeny and a detailed
resolution of the affiliations of the environmental sequences dsrAB should be amplified
from a set of reference species, encompassing several lineages of SOP. Moreover, this
refined framework should enable us to clarify whether Dsr genes, in addition to
undergoing vertical transfer, have also undergone lateral transmission.
The unsuccessful amplification from most of the reference strains might result
from non-optimal PCR conditions despite a sufficient template concentration. The chosen
annealing temperature of 57.1°C appeared optimal since dsrAB was amplified
successfully from the Magnetospirillum gryphiswaldense pure culture. Although
temperature gradient PCR indicated a broad range of annealing temperature a decrease
was not suitable for amplification. Besides PCR conditions, the overall absence of DsrAB
from the respective reference strains might be the reason for amplification failure. During
this study it became obvious, that Rhodobacter sphaeroides does not possess a DsrAB
(www.img.jgi.doe.gov/cgi-bin/pub/main.cgi?page=taxonDetail&taxon_oid=400140000).
This is also likely to be the case for the tested strain Rhodobacter capsulatus. Similarly,
the genome of Thiomicrospira crunogena does not contain Dsr genes
(http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?CMD=search&DB=genome
and
http://www.genome.jp/kegg/pathway.html).
Although the reference strains were selected on the basis of the capability to
oxidize sulfur compounds of biochemical pathways mediating the respective process
without employing DsrAB could be prevalent. For example, the examined reference
organisms Halobacillus neapolitanus, Starkeya novella, Thiobacillus thioparus and
Thiomicrospira crunogena possess the soxB gene (Petri et al., 2001), which might
involve the absence of a DsrAB and thus explain the negative PCR results. Likewise, the
whole Thiomicrospira cluster may employ a DsrAB independent sulfur oxidation
pathway. In contrast, Chlorobium limicola, Chlorobium tepidum and Pelodictyon
phaeoclathratiforme have been shown to harbour both genes (Petri et al., 2001) (Eisen et
al., 2002) and the sulfur-oxidizing gammaproteobacterium Leucothrix mucor most likely
possesses an APS reductase (Grabovich et al., 1999) (Grabovich et al., 2002) suggesting a
possible presence of DsrAB.
52
Literature
4.2 DsrAB and 16S rRNA phylogeny of cultured SOP and
evidence for LGT
All previously available DsrAB sequences of cultured and uncultured organisms
served for the built up of the phylogenetic DsrAB backbone. The addition of
environmental clone sequences stabilized the constructed evolutionary tree.
Overall, the detailed phylogenetic relationships of DsrAB among the alpha-, beta- and
gammaproteobacterial subdivions remain unresolved as indicated by the multifurcations
in the phylogenetic trees. The unique DsrAB sequence of Thiobacillus denitrificans as the
only representative of Betaproteobacteria and to the conspicuous insertion/deletion
pattern, points at a phylogenetic distinct position. Only the insertion shared with
Alkalilimnicola ehrlichei gives evidence for a gammaproteobacterial affiliation. This is
consistent with 16S rRNA gene phylogeny (Woese, 1987) which declares a concomitant
origin of these two subdivisions.
Independent of the treeing method the chlorobi along with Magnetococcus sp.
were found to be the deepest branching SOP. Assuming that the highly conserved 16S
rRNA genes have undergone no lateral transfer and thus the 16S rRNA gene based
phylogeny reflects the true organismal phylogeny (Woese, 1987) events of lateral gene
transfer can be deduced from major inconsistencies between the DsrAB and 16S rRNA
gene tree. Based on this assumption Magnetococcus sp. possess non-orthologous DsrAB
genes since its DsrAB sequence differs from the Magnetospirillum-like sequences that
group with other Proteobacteria. Consistently the putative xenologous DsrAB sequence
of Magnetococcus sp. shares a unique deletion in the alpha subunit typical for all
members of the chlorobi, which can be considered as additional signpost for the separated
phylogenetic position. Thus the phylogenetic affiliation of the non-orthologous
Magnetococcus sp. DsrAB with the presumably orthologous DsrAB of the chlorobi most
likely reflects lateral gene transfer. Analysis of the genomic GC content with regard to
variability towards the dsrAB GC content (Lawrence and Ochman, 1997) could proof this
hypothesis.
The presence of DsrAB in both sulfate-reducing and sulfur-oxidizing organisms
raises the question of its evolutionary origin and the progenitor of Dsr involved in
catalysis of sulphur-oxidation or sulfate reduction. The DsrAB harbouring ancestor of
SOP and SRP was consistently placed between the chlorobi and the deep-branching
sulfate reducers, such as Archaeoglobus fulgidus, Thermodesulfovibrio islandicus and
members of the Firmicutes. Insertion/deletion-patterns shared between all SOP and some
deep-branching SRP (here Thermodesulfovibrio islandicus, Desulfitobacterium
hafniense, Archaeoglobus fulgidus and environmental fosmid clone 39F2) give additional
evidence for the deduced phylogenetic relationship. Moreover, a synopsis of different
treeing methods suggested Archaeoglobus sp. and the chlorobi as the closest relatives of
both functional groups of prokaryotes. With regard to 16S rRNA gene phylogeny
(Woese, 1987) these findings might hint at a lateral transfer of DsrAB genes from
Archaeoglobus fulgidus to the ancestor of the chlorobi since Chlorobium tepidum is not
deeply branching within the bacterial 16S rRNA tree.
To date, the sulfate reducing archaeon Archaeoglobus fulgidus is the oldest
known organism possessing a Dsr but is supposed to have acquired its Dsr genes laterally
from a bacterial donor (Klein et al., 2001). However, the postulated bacterial origin can
not be taken for granted since the employed paralogous rooting method is controversial
Literature
53
and additional evidence from the genome of Archaeoglobus fulgidus identifying the Dsr
genes as clear xenologs is lacking. Here, a paralogous rooting approach of the DsrA and
B subunit of SOP would be worthwhile to approve the deep branching relationships of
SOP. Eisen et al. (2002) postulated lateral gene transfer of archaeal genes to the
thermophilic Chlorobium tepidum based on the findings that many genes resembled those
of the also thermophilic Archaeoglobus fulgidus. Interestingly, the sulfur oxidation genes
of Chlorobium tepidum are more similar to homologous genes in Archaeoglobus fulgidus
than to genes in the sulfur-oxidizing archaea Sulfolobus sp.. In addition, a recent analysis
of the phylogeny and alignment patterns of membrane proteins mediating Dsr dependent
sulfate reduction suggest an affiliation with both Chlorobium tepidum as well as
Archaeoglobus fulgidus (Mussmann et al., 2005). A thermophilic lifestyle of organisms
involved in LGT has already been reported from other studies (Friedrich, 2002). Further
tree topology suggests, after a putative later gene transfer event from Archaeoglobus
fulgidus to the chlorobi, a subsequent transfer of DsrAB to the proteobacterial lineages
involving vertical as well as horizontal tansmission. Indeed, Sabehi et al. (2005) recently
found an uncultured proteorhodopsin-containing bacterium (BAC) to possess a
gammaproteobacterial reverse sulfite reductase operon. In contrast, the proteorhodopsin
sequence and data from other genes supported an alphaproteobacterial affiliation. Hence,
these findings suggest an additional lateral transfer of DsrAB between the alpha- and
gammaproteobacterial subdivisions. However, such putative events within the
Proteobacteria were not detected in our study. This can probably be addressed to the
limited availability of reference sequences.
Overall, our findings yield additional evidence for the close relatedness of the
DsrAB genes involved in either sulfate reduction or sulfur oxidation and their putative
lateral transfer between these functional groups.
S0
H2S/FeS2
SO42-
Figure 25: Limited microbial sulfur cycle before 2.45 billion years (Ga). Partial sulfide oxidation and
sulfur reduction, which are commonly considered to be ancient metabolic processes were present in
addition to localized sulfate reduction. Processes that were possible after oxygen became available are
indicated by red arrows. (modified after Blank, 2004)
One key for understanding the evolutionary origin of Dsr is the early sulfur cycle
(Figure 25). Anaerobic sulfur oxidation is thought to be one of the earliest metabolic
pathways used for energy conservation by ancient prokaryotes. The assumed ancient
chloroflexi have not been shown to possess a DsrAB. In contrast, all recent anoxygenic
phototrophic SOP containing DsrAB group within evolutionary younger phyla.
Regardless of the 16S rRNA gene-based phylogeny comparative sequence
analysis of photosynthesis genes gave evidence for Chlorobium tepidum and
54
Literature
Chloroflexus aurantiacus to be each others closest relatives (Xiong et al., 2000). Thus,
Chlorobium may be more ancient than previously thought and might be a potential
candidate for an ancestral DsrAB. However, a lateral gene transfer of photosynthesis
genes was suggested for this conclusion and may weaken the argumentation (Xiong et al.,
2000; Igarashi et al., 2001). The ‘Wächtershäuser hypothesis’ postulates the formation of
pyrite and hydrogen from ferrous iron sulfide and hydrogen sulfide as ancient energy
yielding reaction for chemolithoautotroph microbial metabolism (Wachtershauser, 1990;
Brock et al., 2006). Moreover sulfur oxidation dependent reduction of ferric iron might
be a primordial, energy yielding process. However, microbial iron sulfide and pyrite
oxidation coupled to electron acceptors other than O2 or NO3-, like Mn(IV) or Fe(III)
could not been detected (Schippers and Jorgensen, 2001, 2002; Blank, 2004) so far. To
trace back the evolutionary pathway of DsrAB the reconstruction of the SOP and SRP
phylogeny should bee seen in context with geological record. Records from sulfur
isotopic fractionation indicate that large scale bacterial sulfate reduction began in marine
environment not until 2.35 billion years ago (Ga) (Cameron, 1982) followed by or nearly
concurrent with rapid oxygenation of the atmosphere about 2.2 Ga (Feng et al., 2000).
Consistent with the availability of oxygen in the biosphere was the diversification of
aerobic sulfide and pyrite oxidizing organism (Blank, 2004). Thus it seems likely that the
Dsr possessing progenitor has to be further searched among either anaerobic sulfide
oxidizers or very early sulfate reducers.
4.3 Diversity of SOP in WS sediment
4.3.1 Overall SOP diversity
One major aim of the study was to assess the diversity of SOP in the Wadden Sea
sediment based on comparative sequence analysis of a key enzyme involved in
biochemical sulfur-oxidation, namely the DsrAB.
Overall, the presented phylogenetic reconstruction of the DsrAB was based on
numerous almost full-length environmental sequences and only on a few sequences of
cultured SOP. The resulting low resolution of the DsrAB tree in comparison to the 16S
rRNA phylogeny limits the conclusiveness of our diversity analysis. For example, the
phototrophic sulfur oxidizer Allochromatium vinosum was found to be the closest
cultured relative of the vast majority of gammaproteobacterial-related sequences.
Furthermore, the putative alphaproteobacterial Wadden Sea clone cluster E affiliated with
Magnetospirillum species, which are currently the only representatives of the
Alphaproteobacteria. This phenomenon can probably be addressed to a lack of DsrAB
sequence data, rather than to close phylogenetic relatedness.
Rarefaction analysis of 112 clones suggests that the actual diversity in our clone
libraries was not completely covered. However the major groups discovered were
identical with those detected in 16S rRNA gene libraries by screening of 543 clones
(Arnds, 2006). Moreover, the analysis was based on the identification of dsrA nucleotide
sequences. Since the genetic code is degenerated unique sequence types on DNA level do
not necessarily result in distinct amino acid sequences. In addition, the high degree of
sequence identity within the clone clusters A, B, D and E suggests the existence of
identical clones which could not be discriminated due to slight sequencing errors. The
Literature
55
close relatedness of DsrAB within each Wadden Sea cluster along with the high sequence
identity suggests more limited species diversity.
Comparative amino acid sequence analysis of the DsrAB revealed a dominance of
gammaproteobacterial sequences among all three clone libraries. In addition, sequences
were identified that can not clearly be assigned to a specific proteobacterial group using
Neighbor Joining, Maximum Parsimony and Maximum Likelihood but the Bayesian
inference method supported an alphaproteobacterial affiliation. These findings are in
good accordance with previous results. 16S rRNA gene based studies (Arnds, 2006) from
the same site already reported the presence of Gammaproteobacteria which affiliated
predominantly with sulfur-oxidizing endosymbionts and some free-living SOP. Arnds
(2006) found numerous 16S rRNA sequences that clustered with the sulfur-oxidizing
gamma 1 (93 – 96% identity) and gamma 2 (96 – 98% identity) symbiont of Olavius
algarvensis. Environmental sequences clustering with the sulfur oxidizing gamma
symbionts of Olavius algarvensis were also detected in the dsrAB clone libraries.
Respectively one clone closely affiliated with the gamma 1 and gamma 2 symbiont.
Additionally, 11 sequences, forming the Wadden sea clone cluster D, grouped with the
gamma 2 symbiont. However, the detected organisms more likely represent permanent
free-living Gammaproteobacteria forming a distinct phylogenetic subgroup together with
the endosymbionts. Unfortunately, the phylogenetic position of the Wadden Sea clone
clusters A and B is of only minor resolution. Here, the phototrophic Allochromatium
vinosum represents the closest cultured relative. Since the prevalence of a
chemolithotrophic lifestyle among sedimentary sulfide oxidizers is more likely than a
phototrophic one, further analysis of pure cultures is necessary to uncover the identity of
the respective organisms.
As supported by tree topology the putative alphaproteobacterial sequences
distantly affiliated with Magnetospirillum sp.. Magnetotactic bacteria are ubiquitous in
marine sediments and their occurrence appears to be dependent on the presence of
opposing gradients of oxidized and reduced compounds, such as oxygen and reduced
sulfur species (Spring and Bazylinsky, 1999-2006). However, members of the genus
Magnetospirillum have only been isolated from freshwater habitats so far (Spring and
Bazylinsky, 1999-2006). In addition, Arnds (2006) did not detect Magnetosprillum-like
sequences. The respective alphaproteobacterial sequences instead clustered with
Roseobacter species and were also frequently detected by FISH. Indeed, members of the
Rosebacter clade harbour the ability to drive on the oxidation of inorganic sulfur
compounds (Sorokin, 1995; Buchan et al., 2005). It can be speculated whether
alphaproteobacterial sequences found in the dsrAB clone libraries belong to Rosebacterlike organisms.
The absence of Betaproteobacteria- and chlorobi-related sequences again is
consistent with findings from 16S rRNA data from the same sampling site (Arnds, 2006).
Similarly, Llobet-Brossa et al. (1998) and Ishii et. al.(2004) reported a minor contribution
of Betaproteobacteria to the bacterial community of the Wadden Sea sediments.
56
Literature
4.3.2 Vertical SOP diversity
In order to reveal potential stratifications of the SOP community over the vertical
sediment profile three clone libraries from different depths were established. Generally, it
is not possible to draw any conclusion about relative in situ abundance from clone
frequencies due to e.g. PCR bias. However, the sequence diversity of DsrAB from
Janssand sediments points at some differences in the vertical distribution of certain
phylotypes. Particularly in the Wadden Sea clone clusters A and B and less pronounced
in the cluster C some sequence subclusters mainly originate either from the oxic/suboxic
surface layers or from anoxic layers below 8 cm depths. The clustering of sequence types
obtained from the surface samples (within Wadden Sea clone clusters A, D, and E) points
at organisms which are capable to use oxygen or alternatively even nitrate as electron
acceptor in sulfur oxidation.
The affiliation of phylotypes from the deeper anoxic sediment layers with each
other (clusters A and D) may indicate sulfur-oxidizers exhibiting an anaerobic life style
such as fermentation or iron or manganese dependent sulfide oxidation. Hence, our
DsrAB based diversity study may reflect the prevalence of different lifestyles as a
consequence of adaptation to different environmental conditions and subsequent niche
differentiation. Nevertheless, the detection of DsrAB sequences in all depths does not
necessarily mean that biological sulfur oxidation proceeds throughout the sediment.
Further analysis has to prove the actual activity of the respective organisms.
The occurrence of Alpha- and Gammaproteobacteria-related sequences
throughout the vertical sediment profile is consistent with previous findings of the same
site based on FISH. Ishii et al. (2004) observed Gammaproteobacteria throughout a depth
of 40 cm. In contrast, Alphaproteobacteria were only found in the upper 10 cm. Arnds
(2006) detected numerous Alpha- and Gammaproteobateria, including potential sulfuroxidizing, up to a sampled depth of 5.5 cm.
In the intertidal sulfur pools the “reverse” Dsr pathway apparently contributes less
to sulfur oxidation, since dsrAB could not be amplified in contrast to 16S rDNA (Arnds,
2006). The minor significance of biological sulfur oxidation by Gammaproteobacteriarelated organisms in these pools is corroborated by the low abundances of
Gammaproteobacteria and lacking 16S rDNA sequences related to known SOP (Arnds,
2006). In contrast, Roseobacter accounted for 3% of all cells in these pools. Accordingly,
amplification of dsrAB should have yielded a product from the sulfur pool if the putative
alphaproteobacterial sequences of our clone libraries indeed represent Roseobacter
relatives. Overall, sulfur oxidation in the intertidal sulfur pool may proceed
predominantly chemically rather than biologically.
In conclusion, our findings are consistent with previous studies based on
comparative 16S rRNA analysis indicating that potential sulfur-oxidizing
Gammaproteobacteria are frequent in marine sediments. It can be proposed that these
organisms are mainly responsible for sulfur oxidation at Janssand site. Especially the
organisms previously targeted by the FISH probe Gam660 (Arnds, 2006) are candidates
for carriers of DsrAB. However this needs to be proven further in future research
projects.
Since it could be demonstrated that LGT is involved in the evolution of SOP the
phylogenetic inferences in the present environmental diversity study may be interpreted
cautiously.
Literature
57
4.4 Quantitative Real-Time PCR
The community ecology of sulfur oxidizers is important for understanding the
relevance of biological sulfur oxidation in the natural habitat. Due to the high
physiological and phylogenetic diversity approaches based on 16S rRNA such as FISH
does only allow the quantification of distinct organisms potentially oxidizing sulfur
compounds. In order to fill this gap a quantitative real-time PCR assay for the dsrB from
SOP was developed to get a first impression of the abundance of microorganism applying
the “reverse” Dsr pathway over the vertical sediment profile. Quantitative real-time PCR
is a semiquantitaive method since organismal abundances, in terms of total cell numbers
are not directly measured by cell counts like in microscopy based methods, e.g. FISH.
Cell abundances are rather deduced from gene copy number. This, of course, is sensitive
for biases, leading to a possible overestimation of cell number. An increase in dsrB copy
number in relation to cell number could result from growing cells, possessing more than
one copy in their genome or multiple replication of the respective gene. Indeed, a
possible replication of the Dsr genes within the genomes of certain SOP could complicate
the interpretation of the qRT-PCR data. Genomic analysis revealed that SRP (Klein et al.,
2001) and the majority of here examined SOP (www.img.jgi.doe.gov/cgibin/pub/main.cgi?page=findGenomes) harbors only a single copy of the dsrAB genes
apart from Chlorobium tepidum possesing two identical dsrAB copies (Eisen et al., 2002).
A particular position takes Thiobacillus denitrificans, harboring three putative copies of
the dsrA (Beller et al., 2006) of which probably one is functional. Hence, besides method
intrinsic biases, the deduction of cell numbers from copy number of dsrA remain critical.
Nevertheless, numerous examples have shown the value of qPCR in environmental
microbiology (Suzuki et al., 2000; Hermansson and Lindgren, 2001; Stubner, 2002; Kolb
et al., 2003).
The established qRT-PCR assay was checked for replicability and specificity.
Sequencing and further PCR based tests verified that fluorescence increase was due to
specific PCR products. However, despite the indicated high specificity of the applied
qPCR primers it should be considered that the general primer bias could result in an
underestimation of copy number. Since amplification proceeded equally for standards and
environmental DNA, as indicated by a similar slope of amplification curves, a high
efficiency of the PCR is assumed. Remarkable are the high standard deviations which
most likely result from the low amount of applied DNA. An increased amount of DNA as
PCR template could minimize pipetting errors and further decrease the threshold cycle
which will probably ensure higher detection sensitivity. For this at least the ten-fold
amount of applied sediment should be used. However, higher amounts of sediment for
DNA extraction bear the risk to introduce PCR inhibitors, e.g. salts, polysaccharides,
humic acids, which will influence the amplification efficiency.
Copy numbers of dsrB determined by qRT-PCR analysis were in the lower range
of FISH detection limit for cell counts of potential sulfur oxidizing Alpha- and
Gammaproteobacteria. Arnds (2006) reported cell numbers of 107 to 108ml-1 for the
respective organisms. Hence, the obtained maximum dsrB copy numbers are in the same
order of magnitude. Theoretically dsrB copy numbers should be as high as cell numbers
since the majority of SOP possess one copy per cell. However, the low DNA template
volumes in combination with a lyses efficiency probably below 100% and general primer
bias might be an explanation for a potential underestimation.
58
Literature
Besides overall coverage of the dsrB copy number, one aim was the assessment of
its distribution over the vertical sediment profile. The high standard deviation does not
allow a conclusion on the vertical stratification of SOP. However, highest copy numbers
were detected in the oxic – anoxic transition zone. These findings are in good match with
previous studies. Arnds (2006) also found maxima for Gammaproteobacteria (9.6 × 108
cells ml-1) and Roseobacter (2.1× 108ml-1) in the sub-oxic transition zone of sulfide and
oxygen. Maximum total cell numbers and microbial sulfate reduction rates in this zone
have already been reported by (Böttcher et al., 2000).
Thus, aerobic SOP extensively consuming available sulfide and oxygen seem to
be abundant. Moreover, the detection of dsrB sequences in anoxic sediment layers
suggests the presence of so far unidentified organisms capable to drive on sulfur
compound oxidation under anaerobic conditions. Based on comparative analysis of
DsrAB or 16S rRNA (Arnds, 2006) none of the clones affiliated with the known
denitrifying species Beggiatoa sp., Thiobacillus denitrificans or Thiomicrospira
dentrificans. SOP using manganese or possibly iron as potential electron acceptors have
not been identified yet (Schippers and Jorgensen, 2001, 2002). Unfortunately, the dsrB
gene based approach does not give any information whether the respective organisms
actively catalyze sulfur oxidation thus further analysis has to reveal this.
In summary, as a first attempt to use dsrB as a quantitative marker for qRT-PCR
the obtained results are of value. To date there is no publication attempting to quantify
sulfur oxidizers independently of their phylogeny. The data suggest that the the “reverse”
Dsr pathway may be very important for biochemical sulfur oxidation in Janssand
sediments.
4.5 General value and limitations of DsrAB gene as
molecular marker for SOP diversity studies and qRTPCR
In the context of this work the value of DsrAB as molecular marker for
environmental diversity studies is limited by the scarce number of sequences from
cultured representatives of SOP. The resulting limited overall resolution of the DsrAB
phylogeny allows only constricted information about the detailed evolutionary
relationships and the environmental diversity of SOP. Thus, additional DsrAB sequences
from known sulfur oxidizers have to be added to the phylogenetic DsrAB backbone to
better resolve the evolutionary relationships in future diversity studies. Complementation
of the phylogenetic framework of the DsrAB gene marker will allow the to link
environmental sequences to recognized lineages of SOP.
If there indeed exists a great variety of metabolic pathways facilitating biological
sulfur oxidation without relying on DsrAB, then comparative sequence analysis of the
respective enzyme will cover only part of the natural microbial diversity. Since there are
obviously several biochemical pathways of sulfur oxidation the DsrAB possessing
organisms only represent a fraction of sulfur oxidizers in a given habitat. However, the
“reverse” Dsr pathway is probably more prevalent in the environment than the sox
pathway (pers. communication C. Dahl 2005) and thus further investment is worthwhile.
Of course, alternative pathways also need to be studied to reveal which one is most
important in marine sediments.
Literature
59
Nevertheless, the problem of the applicability of a functional gene as molecular
marker also holds for other enzymes involved in dissimilatory sulfur metabolism. For
example the flavoprotein sulfide:quinone oxidoreductase (SQR), catalyzing the oxidation
of sulfide to sulfur, is widely distributed among bacteria (Theissen et al., 2003).
However, the poor sequence conservation (average amino acid sequence identity 20%)
and proposed lateral transfer of SQR genes among prokaryotes here again pose a
challenge to the use as molecular marker in phylogenetic analysis. In contrast, the level of
sequence conservation of the soxB gene, present in Alpha-, Beta- and
Gammaproteobacteria as well as Green sulfur bacteria, seems sufficient (amino acid
sequence identity 41 – 98%) but the putative involvement of LGT also holds for this
functional gene (Petri et al., 2001).
In view of the variation of biochemical sulfur oxidation pathways and the wide
distribution of different sulfur oxidizing enzymes throughout the microbial world
(Brune, 1995) it is not likely that an universal molecular marker will be available for all
sulfur oxidizing prokaryotes. Moreover the coexistence of different functional genes in
several SOP complicates a combination of approaches aimed at different molecular
markers.
4.6 Outlook
Overall, the congruence of the present results with 16S rRNA data encourages to
further use dsrAB as molecular marker in environmental diversity studies of SOP.
Therefore, the phylogenetic backbone of DsrAB has to be complemented. Primer
optimization might be necessary for the retrieval of additional DsrAB sequences data
from cultured sulfur oxidizers. In the future, additional genome sequences may be
available. Alternatively, southern blotting of the dsrAB fragment (Dahl et al., 1999)
combined with subcloning and sequencing could yield the respective sequences from
selected reference strains. Based on both enhanced primers can be designed.
Since the comparative sequences analysis of DsrAB gives no information about
the identity of uncultured environmental SOP, the direct linkage of a functional gene to a
16S rRNA phylotype is one of the main future goals. Here, application of fluorescence
associated cell sorting (FACS) (Sekar et al., 2004) of certain subgroups, e.g. hybridized
with 16S rRNA FISH probes GAM 660 or ROS 537, and subsequent amplification of
dsrAB from the sorted cell fractions is a promising approach. Moreover, screening of
dsrAB in metagenome libraries combined with screening for 16S rRNA genes and
subsequent analysis of tetranucleotide patterns (Teeling et al., 2004)could unravel the
identity of uncultured environmental sulfur oxidizers. Finally, the design of group
specific probes could elucidate the community structure of sulfur oxidizers and
furthermore allow an evaluation of the role of specific groups by quantitative FISH
(Moraru, 2006).
In future studies the expression of dsrAB may be tracked by determination of dsrB
mRNA copy number under distinct environmental conditions. Therefore, the application
of reverse transcription PCR prior to quantitative-real time PCR constitutes a promising
method (Neretin et al., 2003). In addition, immunocytochemistry using antibodies raised
against DsrAB could overcome the expression problem of DNA and mRNA targeted
assays (Pernthaler and Amann, 2004) and detect the presence of DsrAB in situ.
60
Literature
The application of the newly developed dsrB targeted qRT-PCR assay for the first time
enabled quantitative assessment of a specific physiological group of SOP. The partial
congruence with CARD-FISH and clone library data (Arnds, 2006) encourages for
further surveys. Since the relatively lower sensitivity of the SybrGreen compared to the
TaqMan assay has been reported (Stubner, 2002) the application of a probe-based
approach should be considered.
To confirm the relevance of the “reverse” Dsr pathway in Wadden Sea sediments,
the role of competing pathways such as the “sox PSO pathway” has to be explored by e.g.
soxB targeted qRT-PCR. Overall, the improvement of the current protocol is a promising
approach for the quantitative detection of SOP in sediments of the Wadden Sea and other
environments.
Literature
61
(C) Literature
Altschul, S.F., Gish, W., Miller, W., Myers, E.W., and Lipman, D.J. (1990) Basic local
alignment search tool. J. Mol. Biol. 215: 403-410.
Arnds, J. (2006) Molekularbiologische Identifizierung potentiell schwefeloxidierender
Bakterien in Oberflächensedimenten des Wattenmeeres. In Department of Biology.
Mainz: Johannes-Gutenberg Universität, p. 84.
Beller, H.R., Chain, P.S.G., Letain, T.E., Chakicherla, A., Larimer, F.W., Richardson,
P.M. et al. (2006) The genome sequence of the obligately chemolithoautotrophic,
facultatively anaerobic bacterium Thiobacillus denitfificans. Journal of Bacteriology
188: 1473-1488.
Blank, C.E. (2004) Evolutionary timing of the origins of mesophilic sulphate reduction
and oxygenic photosynthesis. a phylogenomic dating approach. Geobiology 2: 1-20.
Böttcher, M.E., Hespenheide, B., Llobet-Brossa, E., Beardsley, C., Larsen, O., Schramm,
A. et al. (2000) The biogeochemistry, stable isotope chemistry, and the microbial
community structure of a temperate intertidal mud: an integrated study. Cont. Shelf
Res. 20: 1749-1769.
Boucher, Y., Douady, C.J., Papke, R.T., Walsh, D.A., Bodreau, M.E., C.L., N. et al.
(2003) Lateral gene transfer and the origins of prokaryotic groups. Annu. Rev. Genet.
37: 283-328.
Brinkhoff, T., Santegoeds, C.M., Sahm, K., Kuever, J., and Muyzer, G. (1998) A
polyphasic approach to study the diversity and vertical distribution of sulfur-oxidizing
Thiomicrospira species in coastal sediments of the German Wadden Sea. Appl.
Environ. Microbiol. 64: 4650-4657.
Brock, T.D., Madigan, M.T., and Martinko, J.M. (2006) Biology of Microorganisms.
Upper Saddle River, N.J. [u.a.]: Prentice-Hall.
Brune, D.C. (1995) Sulfur compounds as photosynthetic electron donors. In Anoxygenic
Photosynthetic Bacteria. Bauer, C.E. (ed): Kluwer, Dordrecht, pp. 847-870.
Brüser, T., Lens, P.N.L., and Trüper, H.G. (2000) The biologicl sulfur cycle. In
Environmental technnologies to treat sulfur pollution. Lens, P.N.L. (ed). London:
IWA Publishing, pp. 47 - 76.
Buchan, A., Gonzalez, J.M., and Moran, M.A. (2005) Overview of the marine
Roseobacter lineage. Applied and Environmental Microbiology 71: 5665-5677.
Cameron, E.M. (1982) Sulfate and Sulfate Reduction in Early Precambrian Oceans.
Nature 296: 145-148.
Dahl, C., and Trüper, H.G. (1994) Enzymes of Dissimilatory Sulfide Oxidation in
Phototrophic Sulfur Bacteria. In Inorganic Microbial Sulfur Metabolism, pp. 400421.
Dahl, C., Engels, S., Pott-Sperling, A.S., Schulte, A., Sander, J., Luebbe, J. et al. (2005)
Novel Genes of the dsr Gene Cluster and Evidence for Close Interaction of Dsr
Proteins during Sulfur Oxidation in the Phototrophic Sulfur Bacterium
Allochromatium vinosum. Journal of Bacteriology 187: 1392–1404.
Dahl, C., Rakhely, G., Pott-Sperling, A.S., Fodor, B., Takacs, M., Toth, A. et al. (1999)
Genes involved in hydrogen and sulfur metabolism in phototrophic sulfur bacteria.
Fems Microbiology Letters 180: 317-324.
Doolittle, W.F. (1999) Lateral genomics. Trends in Biochemical Sciences 24: M5-M8.
62
Literature
Dorak, M.T. Real-Time PCR. In: www.dorak.info/genetics/realtime.ppt.
Eisen, J.A., Nelson, K.E., Paulsen, I.T., Heidelberg, J.F., Wu, M., Dodson, R.J. et al.
(2002) The complete genome sequence of Chlorobium tepidum TLS, a
photosynthetic, anaerobic, green-sulfur bacterium. Proceedings of the National
Academy of Sciences of the United States of America 99: 9509-9514.
Feng, L.R., Donaldson, J.A., and Holland, H.D. (2000) Alteration rinds on glacial
diamictite clasts in the Gowganda formation: Possible indicators of low atmospheric
oxygen ca. 2.3 Ga. International Geology Review 42: 684-690.
Friedrich, C.G., Rother, D., Bardischewsky, F., Quentmeier, A., and Fischer, J. (2001)
Oxidation of reduced inorganic sulfur compounds by bacteria: Emergence of a
common mechanism? Applied and Environmental Microbiology 67: 2873-2882.
Friedrich, C.G., Bardischewsky, F., Rother, D., Quentmeier, A., and Fischer, J. (2005)
Prokaryotic sulfur oxidation. Current Opinion in Microbiology 8: 253-259.
Friedrich, M.W. (2002) Phylogenetic analysis reveals multiple lateral transfers of
adenosine-5 '-phosphosulfate reductase genes among sulfate- reducing
microorganisms. Journal of Bacteriology 184: 278-289.
Gogarten, J.P., Doolittle, W.F., and Lawrence, J.G. (2002) Prokaryotic evolution in light
of gene transfer. Molecular Biology and Evolution 19: 2226-2238.
Grabovich, M.Y., Dul'tseva, N.M., and Dubinina, G.A. (2002) Carbon and sulfur
metabolism in representatives of two clusters of bacteria of the genus Leucothrix: A
comparative study. Microbiology 71: 255-261.
Grabovich, M.Y., Muntyan, M.S., Lebedeva, V.Y., Ustiyan, V.S., and Dubinina, G.A.
(1999) Lithoheterotrophic growth and electron transfer chain components of the
filamentous gliding bacterium Leucothrix mucor DSM 2157 during oxidation of
sulfur compounds. Fems Microbiology Letters 178: 155-161.
Hatchikian, E.C., and Zeikus, J.G. (1983) Characterization of a New Type of
Dissimilatory Sulfite Reductase Present in Thermodesulfobacterium-Commune.
Journal of Bacteriology 153: 1211-1220.
Heck, K.L., van Belle, G., and Simberloff, D. (1975) Explicit calculation of the
rarefaction diversity measurement and the determination of sufficient sample size.
Ecology 56: 1459-1461.
Hermansson, A., and Lindgren, P.E. (2001) Quantification of ammonia-oxidizing
bacteria in arable soil by real-time PCR. Applied and Environmental Microbiology
67: 972-976.
Igarashi, N., Harada, J., Nagashima, S., Matsuura, K., Shimada, K., and Nagashima,
K.V.P. (2001) Horizontal transfer of the photosynthesis gene cluster and operon
rearrangement in purple bacteria. Journal of Molecular Evolution 52: 333-341.
Ishii, K., Mussmann, M., MacGregor, B.J., and Amann, R. (2004) An improved
fluorescence in situ hybridization protocol for the identification of bacteria and
archaea in marine sediments. Fems Microbiology Ecology 50: 203-212.
Kelly, D.P., Shergill, J.K., Lu, W.P., and Wood, A.P. (1997) Oxidative metabolism of
inorganic sulfur compounds by bacteria. Antonie Van Leeuwenhoek International
Journal of General and Molecular Microbiology 71: 95-107.
Klein, M., Friedrich, M., Roger, A.J., Hugenholtz, P., Fishbain, S., Abicht, H. et al.
(2001) Multiple lateral transfers of dissimilatory sulfite reductase genes between
Literature
63
major lineages of sulfate-reducing prokaryotes. Journal of Bacteriology 183: 60286035.
Kletzin, A., Urich, T., Muller, F., Bandeiras, T.M., and Gomes, C.M. (2004)
Dissimilatory oxidation and reduction of elemental sulfur in thermophilic archaea.
Journal of Bioenergetics and Biomembranes 36: 77-91.
Kolb, S., Knief, C., Stubner, S., and Conrad, R. (2003) Quantitative detection of
methanotrophs in soil by novel pmoA-targeted real-time PCR assays. Applied and
Environmental Microbiology 69: 2423-2429.
Lane, D.J., Harrison, A.P.J., Stahl, D., Pace, B., Giovannoni, S.J., Olsen, G.J., and Pace,
N.R. (1992) Evolutionary relationships among sulfur- and iron-oxidizing eubacteria.
J. Bacteriol. 174: 269-278.
Lawrence, J.G., and Ochman, H. (1997) Amelioration of bacterial genomes: Rates of
change and exchange. Journal of Molecular Evolution 44: 383-397.
Llobet-Brossa, E., Rossello-Mora, R., and Amann, R. (1998) Microbial community
composition of Wadden Sea sediments as revealed by fluorescence in situ
hybridization. Appl. Environ. Microbiol. 64: 2691-2696.
Llobet-Brossa, E., Rabus, R., Böttcher, M.E., Könneke, M., Finke, N., Schramm, A. et al.
(2002) Community structure and activity of sulfate-reducing bacteria in an intertidal
surface sediment: a multi-method approach. Aquat. Microb. Ecol. 29: 211-226.
Markert, S., Arndt, C., Felbeck, H., Becher, D., Albrecht, D., Hecker, M., and Schweder,
T. (2005) Functional genome analysis of the bacterial endosymbiont from deep sea
tubeworm Riftia pachyptila. In Poster.
Molitor, M., Dahl, C., Molitor, I., Schafer, U., Speich, N., Huber, R. et al. (1998) A
dissimilatory sirohaem-sulfite-reductase-type protein from the hyperthermophilic
archaeon Pyrobaculum islandicum. Microbiology-Uk 144: 529-541.
Moraru, C.L. (2006) Development of protocols for the in situ detection of genes in
microorganisms. In IMPRS for marine Microbiology. Bremen: Bremen.
Mülhardt, C. (1999) Quantitative PCR. In Der Experimentator: Molekularbiologie.
Stuttgart, Jena, Lübeck, Ulm: Gustav Fischer Verlag, p. 224.
Muller, F.H., Bandeiras, T.M., Urich, T., Teixeira, M., Gomes, C.M., and Kletzin, A.
(2004) Coupling of the pathway of sulphur oxidation to dioxygen reduction:
characterization of a novel membrane-bound thiosulphate : quinone oxidoreductase.
Molecular Microbiology 53: 1147-1160.
Mussmann, M., Ishii, K., Rabus, R., and Amann, R. (2005) Diversity and vertical
distribution of cultured and uncultured Deltaproteobacteria in an intertidal mud flat of
the Wadden Sea. Environmental Microbiology 7: 405-418.
Mussmann, M., Schulz, H.N., Strotmann, B., Kjaer, T., Nielsen, L.P., Rossello-Mora,
R.A. et al. (2003) Phylogeny and distribution of nitrate-storing Beggiatoa spp. in
coastal marine sediments. Environ. Microbiol 5: 523-533.
Neretin, L.N., Schippers, A., Pernthaler, A., Hamann, K., Amann, R., and Jorgensen,
B.B. (2003) Quantification of dissimilatory (bi)sulphite reductase gene expression in
Desulfobacterium autotrophicum using real-time RT-PCR. Environmental
Microbiology 5: 660-671.
Pernthaler, A., and Amann, R. (2004) Simultaneous fluorescence in situ hybridization of
mRNA and rRNA in environmental bacteria. Applied and Environmental
Microbiology 70: 5426-5433.
64
Literature
Petri, R., Podgorsek, L., and Imhoff, J.F. (2001) Phylogeny and distribution of the soxB
gene among thiosulfate-oxidizing bacteria. Fems Microbiology Letters 197: 171-178.
Pott, A.S., and Dahl, C. (1998) Sirohaem sulfite reductase and other proteins encoded by
genes at the dsr locus of Chromatium vinosum are involved in the oxidation of
intracellular sulfur. Microbiology-Uk 144: 1881-1894.
Rohwerder, T., and Sand, W. (2003) The sulfane sulfur of persulfides is the actual
substrate of the sulfur-oxidizing enzymes from Acidithiobacillus and Acidiphilium
spp. Microbiology-Sgm 149: 1699-1709.
Sabehi, G., Loy, A., Jung, K.H., Partha, R., Spudich, J.L., Isaacson, T. et al. (2005) New
insights into metabolic properties of marine bacteria encoding proteorhodopsins. Plos
Biology 3: 1409-1417.
Saiki, R.K., Gelfand, D.H., Stoffel, S., Scharf, S.J., Higuchi, R., Horn, G.T. et al. (1988)
Primer-directed enzymatic amplification of DNA with a thermostable DNA
polymerase. Science 239: 487-491.
Schedel, M., and Truper, H.G. (1979) Purification of Thiobacillus-Denitrificans
Siroheme Sulfite Reductase and Investigation of Some Molecular and Catalytic
Properties. Biochimica Et Biophysica Acta 568: 454-466.
Schedel, M., Vanselow, M., and Truper, H.G. (1979) Seroheme Sulfite Reductase
Isolated from Chromatium-Vinosum - Purification and Investigation of Some of Its
Molecular and Catalytic Properties. Archives of Microbiology 121: 29-36.
Schippers, A., and Jorgensen, B.B. (2001) Oxidation of pyrite and iron sulfide by
manganese dioxide in marine sediments. Geochimica Et Cosmochimica Acta 65: 915922.
Schippers, A., and Jorgensen, B.B. (2002) Biogeochemistry of pyrite and iron sulfide
oxidation in marine sediments. Geochimica Et Cosmochimica Acta 66: 85-92.
Schulz, H.N., Strotmann, B., Gallardo, V.A., and Jørgensen, B.B. (1999) Population
study on the filamentous sulfur bacteria Thioploca spp. off the Bay of Concepcion,
Chile. Mar. Ecol. Prog. Ser. 200: 117-126.
Sekar, R., Fuchs, B.M., Amann, R., and Pernthaler, J. (2004) Flow sorting of marine
bacterioplankton after fluorescence in situ hybridization. Applied and Environmental
Microbiology 70: 6210-6219.
Sorokin, D.Y. (1995) Sulfitobacter Pontiacus Gen-Nov, Sp-Nov - a New Heterotrophic
Bacterium from the Black-Sea, Specialized on Sulfite Oxidation. Microbiology 64:
295-305.
Spring, S., and Bazylinsky, D.A. (1999-2006) The Prokaryotes online
http://141.150.157.117:8080/prokPUB/index.htm. In. New York: Springer.
Stubner, S. (2002) Enumeration of 16S rDNA of Desulfotomaculum lineage 1 in rice
field soil by real-time PCR with SybrGreen (TM) detection. Journal of
Microbiological Methods 50: 155-164.
Sun, C.W., Chen, Z.W., He, Z.G., Zhou, P.J., and Liu, S.J. (2003) Purification and
properties of the sulfur oxygenase/reductase from the acidothermophilic archaeon,
Acidianus strain S5. Extremophiles 7: 131-134.
Suzuki, M.T., Taylor, L.T., and DeLong, E.F. (2000) Quantitative analysis of smallsubunit rRNA genes in mixed microbial populations via 5 '-nuclease assays. Applied
and Environmental Microbiology 66: 4605-4614.
Literature
65
Teeling, H., Meyerdierks, A., Bauer, M., Amann, R., and Glockner, F.O. (2004)
Application of tetranucleotide frequencies for the assignment of genomic fragments.
Environmental Microbiology 6: 938-947.
Theissen, U., Hoffmeister, M., Grieshaber, M., and Martin, W. (2003) Single eubacterial
origin of eukaryotic Sulfide : Quinone oxidoreductase, a mitochondrial enzyme
conserved from the early evolution of eukaryotes during anoxic and sulfidic times.
Molecular Biology and Evolution 20: 1564-1574.
Urich, T. (2005) The sulfur oxygenase reductase from acidianus ambivalens. In Biology.
Darmstadt: Technical University, p. 151.
Wachtershauser, G. (1990) The Case for the Chemoautotrophic Origin of Life in an IronSulfur World. Origins of Life and Evolution of the Biosphere 20: 173-176.
Wagner, M., Roger, A.J., Flax, J.L., Brusseau, G.A., and Stahl, D.A. (1998) Phylogeny of
dissimilatory sulfite reductases supports an early origin of sulfate respiration. Journal
of Bacteriology 180: 2975-2982.
Wagner, M., Loy, A., Klein, M., Lee, N., Ramsing, N.B., Stahl, D.A., and Friedrich,
M.W. (2005) Functional marker genes for identification of sulfate-reducing
prokaryotes. In Environmental Microbiology, pp. 469-489.
Woese, C.R. (1987) Bacterial evolution. Microbiological Reviews 51: 221-271.
Wolfe, B.M., Lui, S.M., and Cowan, J.A. (1994) Desulfoviridin, a MultimericDissimilatory Sulfite Reductase from Desulfovibrio-Vulgaris (Hildenborough) Purification, Characterization, Kinetics and Epr Studies. European Journal of
Biochemistry 223: 79-89.
Xiong, J., Fischer, W.M., Inoue, K., Nakahara, M., and Bauer, C.E. (2000) Molecular
evidence for the early evolution of photosynthesis. Science 289: 1724-1730.
Appendix
66
(D) Appendix
Allochromatium vinosum
Fosmid 194A5
Alkalilimnicola ehrlichei
Environmental sequence IBEA_CTG_1982486
Proteorhodopsin clone MED13K9
Magnetospirillum magnetotacticum
Magnetospirillum magneticum
Magnetococcus sp. MC-1
Chlorobium limicola
Prosthecochloris aestuarii
Prosthecochloris vibrioformis
Chlorobium tepidum
Chlorobium phaeobacteroides EAM34529
Chlorobium phaeobacteroides EAM62784
Pelodictyon phaeoclathratiforme
Thiobacillus denitrificans AAZ98438
Thiobacillus denitrificans AAZ97262
Thiobacillus denitrificans AAZ97322
Olavius algarvensis Gamma2 symbiont
Olavius algarvensis Gamma1 symbiont
116f
Allochromatium vinosum
Fosmid 194A5
Alkalilimnicola ehrlichei
Environmental sequence IBEA_CTG_1982486
Proteorhodopsin clone MED13K9
Magnetospirillum magnetotacticum
Magnetospirillum magneticum
Magnetococcus sp. MC-1
Chlorobium limicola
Prosthecochloris aestuarii
Prosthecochloris vibrioformis
Chlorobium tepidum
Chlorobium phaeobacteroides EAM34529
Chlorobium phaeobacteroides EAM62784
Pelodictyon phaeoclathratiforme
Thiobacillus denitrificans AAZ98438
Thiobacillus denitrificans AAZ97262
Thiobacillus denitrificans AAZ97322
Olavius algarvensis Gamma2 symbiont
Olavius algarvensis Gamma1 symbiont
240
CTCGAGACCGGTCCCTGGCCCAGCTTCATCTCCGGCATCAAGCGTCTGCG
CTGGAAGAAGGCCCCTGGCCGAGCTTCATCAGCGGTATCAAGCGTCTGCG
CTGGAAAGCGGCCCCTGGCCCAGTTTTGTCACCGGCCTGAAGCACCTGGC
CTAGAAACTGGTCCATGGCCTAGTTTCATTAGCGGAATTAAAAAACTCAG
CTTGAGAATGGCCCATGGCCAAGCTTCATAAGCGGCATTAAAAGACTCAG
TTGGAAAGCGGTCCCTGGCCCAGCTTCGTGACCGGCTTGAAGCGTCTGGC
CTGGAAAGCGGTCCCTGGCCCAGCTTCGTGACGGGCCTGAAGCGTCTGGC
CTCGAAGAGGGACCCTGGCCCAGTTTCATCTCCGGGTTCAAAGAACTCTA
CTGGAAAGCGGCCCCTGGCCAAGCTTCGTTTCCGGCTTCAAGGAACTTGC
CTTGAAAGCGGCCCCTGGCCGAGCTTCGTATCCGGTTTCAAGGAGATGGC
CTGGAAAAAGGCCCCTGGCCAAGTTTTGTCTCTGGCTTCAAGACTCTTGC
CTTGAAAGCGGTCCGTGGCCCAGCTTCATTTCGGGCTTCAAGGCGCTTGC
CTGGAAAAAGGCCCCTGGCCGAGTTTTATTTCCGGTTTCAAAGAGCTTGC
CTTGAAAGCGGCCCCTGGCCGAGTTTTGTGACCGGTCTGAAGGACCTTGC
TTGGAAAAAGGTCCTTGGCCGAGTTTTATCAGCGGTTTCAAGGAACTTGC
CTCGAGCGCGGACCCTGGCCCAGCTTCGTCACCGGGCTCAAGCGTCTGGC
CTCGAGAAAGGCCCGTGGCCCAGCTTCGTGACCGGGCTGAAGCGTTTGGC
CTCGAGAAAGGGCCTTGGCCCAGCTTCGTGACCGGCCTCAAGCGCCTCGC
CTCGAGAACGGTCCCTGGCCAAGCTTTATTTCAGGGATCAAGCGTTTACG
TTGGAAACCGGCCCCTGGCCGAGTTTCATTTCCGGTATCAAGCGTCTGCG
GGNCCNTGGCCBAGYTT-----------------------AGACCCGCAAGGGCTACTGGAAGGGCGGTACCGTGTCGGTGTT
AGACCCGCAAGGGTTACTGGAAAGGAGGCACCGTCTCTGTCTA
AGGACCGCAAGGGTTACTGGAAGGGCGGCACCGTCAGTGTCTA
AAACCCGCAAAGGATATTGGAAGGGTGGCACTGTCAGTGTTTA
AAACACGTAAAGGGTATTGGAAAGGCGGCACAGTGAGTGTATA
AGACCCGCAAGGGCTATTGGAAGGGCGGCACGGTGGGCGTGTT
AGACCCGCAAGGGCTATTGGAAGGGCGGCACGGTGGGCGTGTT
AAACCCGCATGGGTTACTGGAAAGGTGGCGTGGTAGGTGTACA
ACACGAAAATGGGGTACTGGAAGGGCGGTCTTGTTACGGTTGA
ACACCAAAATGGGCTACTGGAAAGGCGGTCTGGTGACTGTTGA
ATACCAAGATGGGTTACTGGAAGGGTGGGCTTGTGACGGTTGA
AGACCAAGATGGGCTACTGGAAAGGCGGTCTGGTCACGGTAGA
ATACCAAGATGGGGTATTGGAAGGGCGGTCTGGTAACGGTGGA
AGACAAGGATGGGGTATTGGAAGGGCGGTCTTGTTACCGTGAA
ACACCAAGATGGGCTACTGGAAAGGCGGTCTTGTTACTGTTGA
AGACCAAGACCGGCTACTGGAAAGGCGGCACGGTCGGGGTGTT
AGACCAAGTTCGGCTACTGGAAGGGCGGCACGGTCGGCGTGTT
AGACCAAGAAAGGCTATTGGAAGGGCGGCACGGTCGGCGTGTT
AGACCCGTAAGGGTTACTGGAAAGGCGGCACCGTCAGCGTCTA
AGACCCGCAAGGGGTATTGGAAGGGTGGCACCGTTTCGGTGTA
GGNTAYTGGAARGGYGG---------------
Figure A1: Sequence alignments of primer binding sites for dsrA
Appendix
Allochromatium vinosum
Fosmid 194A5
Alkalilimnicola ehrlichei
Environmental sequence IBEA_ CTG_2027414
Environmental sequence IBEA_CTG_1982486
Environmental sequence IBEA_CTG_2018072
Proteorhodopsin clone MED13K9
Magnetospirillum magnetotacticum
Magnetospirillum magneticum
Magnetococcus sp. MC-1
Chlorobium limicola
Prosthecochloris aestuarii
Prosthecochloris vibrioformis
Chlorobium phaeobacteroides
Chlorobium tepidum
Chlorobium phaeobacteroides
Pelodictyon phaeoclathratiforme
Thiobacillus denitrificans
Olavius algarvensis Gamma2 symbiont
Olavius algarvensis Gamma1 symbiont
404_for
Allochromatium vinosum
Fosmid 194A5
Alkalilimnicola ehrlichei
Environmental sequence IBEA_ CTG_2027414
Environmental sequence IBEA_CTG_1982486
Environmental sequence IBEA_CTG_2018072
Proteorhodopsin clone MED13K9
Magnetospirillum magnetotacticum
Magnetospirillum magneticum
Magnetococcus sp. MC-1
Chlorobium limicola
Prosthecochloris aestuarii
Prosthecochloris vibrioformis
Chlorobium phaeobacteroides
Chlorobium tepidum
Chlorobium phaeobacteroides
Pelodictyon phaeoclathratiforme
Thiobacillus denitrificans
Olavius algarvensis Gamma2 symbiont
Olavius algarvensis Gamma1 symbiont
dsr_808_rev (Target sequence)
CATGATCTCGCACACCCAGGGCTGGCTGCACTGCGACATCCCG
CATGATCTCGCACACCCAGGGCTGGCTGCATTGCGACATCCCC
CATGATCAGCCACACCCAGGGCTGGCTGCACTGTGACATCCCA
TATGATTTCTCATACACAGGGCTGGTTACACTGTGATATACCA
AATGATTTCACATACTCAAGGTTGGTTACATTGTGATATACCC
TATGATTTCTCACACCCAAGGTTGGCTGCATTGTGATATTCCT
AATGATTTCACATACACAAGGATGGTTACATTGTGATATCCCA
CTTCATCAGCCATACCCAGGGCTGGCTGCACTGCGACATCCCC
CTTCATCAGCCATACCCAGGGCTGGCTGCACTGCGACATCCCC
CTCCATGGCCCACACCCAGGGCTGGTTGCACTGTGACATTCCT
TTCGGTATCGCACACGCAGGGATGGCTGCACTGCGATATTCCC
CTCGGTCTCCCATACCCAGGGTTGGCTCCATTGCGATATTCCT
GTCGGTATCGCACACCCAGGGGTGGCTCCATTGCGACATTCCC
ACCGATAGCGCATACGCAGGGATGGCTCCACTGCGATATTCCG
AGCCGTTTCGCACACCCAGGGCTGGCTGCACTGCGACATTCCG
CTCGGTGTCGCATACGCAGGGTTGGCTTCATTGCGATATTCCT
TTCTGTTTCTCATACACAAGGTTGGCTTCATTGCGATATTCCG
GATGATCGCCCACACCCAGGGCTGGCTGCACTGCGACATCCCC
GATGATATCGCACACCCAGGGCTGGCTGCACTGCGACATTCCG
CATGATCTCCCATACCCAGGGCTGGTTGCACTGCGACATCCCG
CAYACNCARGGNTGGYT---------------CGCCGA-TGCAGATCAACGACGCCGAGCACTCCAAGCTGGCCATCTGGGTCGGCGGCAACCACTCCAACGCGCGCGGCACGCCGA-TGCAGATCAACGACGCCGAACACACCAAGCTCGCCATCTGGGTGGGCGGTAACCACTCCAATGCACGTGGCACCCCCA-TGGAGATCAACGACCCGGAGCACTCCAAACTGGCGATCTGGATCGGGGGCAACCACTCCAACGCCCGTGGCGCTCCAA-TGCAAATCAATGACCCTGAGCACAGTAAACTTGCTATTTGGGTTGGTGGTAATCACTCTAATGCAAGAGGCACTCCTA-TGCAGATCAATGACCCCGAACATAGTAAATTAGCTATATGGGTTGGCGGCAACCACTCGAATGCGAGAAGCACTCCAA-TGCAAATCAATGATGCCGAACATACAAAACTTTCAATCTGGGTCGGTGGTAATCATTCAAATGCGAGAGGAACACCTA-TGCAAATTAATGATCCAGAGCATAGTAAGCTTGCAATATGGGTAGGTGGAAACCACTCAAATGCAAGGAGTACGCCCA-TGCAGATCAACGATCCCATCCATTCCAAGATCGCCATCTGGGTGGGCGGCAAGCATTCCTCTACCCGGTCCACGCCCA-TGCAGATCAACGACCCGGTCCATTCCAAGATCGCCATCTGGGTGGGCGGCAAGCACTCCTCCACCCGTTCCACTGCGA-TGGAGATCAACCATCCTGAGTACTCCAAACTGGCTATCTGGATCGGTGGTAAGAACGCCAACACCCGTTCCCCCGCCA-TGGAGATCAACCATCCCGAGCACTCCAAGTTTGCGATCTGGGTCGGCGGCAAGAACAGCAACGCCCGTTCCACTGCAA-TGGAGATCAACCATCCCGAGCATTCAAAATTCGCTATCTGGGTAGGAGGCAAGAACAGTAATGCCCGATCCACGGCAA-TGGAAATAAATCATCCCGAGCACTCCAAGTTTGCAATCTGGGTTGGCGGAAAGAACAGCAACGCCCGCTCAACTTCGA-TGGAGATCAACCATCCTGAATACTCCAGGTTTGCTATCTGGGTTGCCGGCAAGAACAGTAATGCCCGTTCAACGGCCA-TGGAGATCAACCACCCCGAGCACTCCAAGTTTGCCGTCTGGGTCGGCGGCAAAAACAGCAACGCCCGCTCGACTGCCA-TGGAGATCAACCATCCGGAACACTCAAAGTTTGCCATATGGGTAGGCGGCAAAAACAGCAATGCCCGTTCAACTTCCA-TGGAGATCAACCACCCGGAACATTCCAAATTTGCAATCTGGGTCGGCGGGAAAAACAGCAATGCCCGATCAACGCCGA-TGCAGATCAACGATCCCGAGCACTCGAAGATCGCGATCTGGGTTGGCGGCAAA----AACTCGAA-------CGCCGA-TGCAAATCAACGACCCCGAACATACCAAGCTAGCTATCTGGGTCGGCGGTAACCACTCCAACGCGCGCGGCACGCCCA-TGCAGATCAACGATGCCGAGCATTCCAAGCTGGCCATCTGGGTCGGCGGTAACCACTCCAACGCCCGCGGCA------------------------------------GCNATHTGGGTNGGHGG
Figure A2: Sequence alignments of of primer binding sites for dsrB
67
Appendix
68
Table A1: Sequence similarities1 of amino acid (lower triangle) and nucleotide sequences (upper triangle) of the DsrAB
Reference species
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
1
Beggiatoa (uncultured)
M. magnetotacticum
M. magneticum
T. denitrificans
O. algarvensis gam 2
O. algarvensis gam 1
A. vinosum
A. ehrlichei
C. phaeobacteroides DSM266
C. limicola
P. aestuarii
P. vibrioformis
P. phaeoclathratiforme
C. phaeobacteroides BS1
C. tepidum
M. spec.
M. gryphiswaldense
R. pachyptila symbiont
Sequence similarity [ % ]
1
2
3
4
25.3
24.6
23.9
68.4
92.8
32.3
67.4
96.4
33.1
56.8
64.5
64.0
77.9
66.0
65.7
63.4
68.8
69.3
67.2
79.8
67.1
68.0
65.4
74.7
69.0
69.6
66.9
50.0
56.7
57.1
55.2
45.7
56.4
56.4
54.5
52.1
57.1
57.2
55.3
50.0
58.1
57.7
54.7
48.9
57.1
57.1
54.8
52.1
55.6
56.3
55.1
54.3
55.4
55.5
54.2
39.8
57.7
58.3
54.2
65.3
87.5
88.7
63.9
85.0
63.4
63.8
55.6
5
23.9
50.0
51.1
32.2
84.6
80.6
73.3
56.2
55.3
57.8
56.3
56.6
55.6
57.5
54.5
66.4
63.6
6
23.9
43.4
45.1
36.2
48.9
91.7
74.7
60.3
60.0
61.1
59.9
60.8
58.8
59.9
64.2
69.4
50.4
7
22.5
71.7
72.9
33.2
56.8
52.9
76.1
59.1
57.8
59.4
57.8
58.7
56.5
59.4
57.6
68.6
63.0
8
26.3
31.8
33.1
35.6
32.2
35.1
32.9
58.4
57.4
58.0
58.5
58.7
57.7
57.7
56.1
70.4
68.5
9
25.9
38.6
39.3
28.1
38.9
41.4
39.5
29.5
90.5
90.0
86.5
88.9
74.8
87.6
66.3
58.5
47.8
10
28.0
40.4
41.6
30.3
39.3
42.6
41.1
30.6
81.1
90.4
86.3
86.5
73.8
87.0
64.4
57.5
45.9
11
27.6
40.3
40.9
29.6
40.5
41.9
41.7
30.4
80.1
80.9
86.5
86.5
76.0
89.8
66.8
58.5
47.8
12
28.3
39.6
40.8
31.1
39.3
43.3
41.2
33.0
55.0
56.7
56.9
85.6
75.3
85.3
65.9
57.5
46.1
13
27.3
37.7
38.7
27.3
37.7
40.1
38.9
28.5
79.3
77.6
76.8
54.2
74.5
84.3
64.1
57.6
47.0
Similarity values were calculated using the Similarity Matrix calculation editor from the ARB software package.
14
29.7
38.3
39.3
28.3
39.4
39.1
39.9
28.9
73.2
72.2
73.1
53.4
72.2
74.6
63.3
56.5
47.0
15
24.9
42.4
43.8
32.3
42.1
45.6
44.2
34.3
77.4
79.7
79.4
57.2
73.4
70.8
65.4
57.5
49.1
16
27.3
38.9
39.7
30.4
38.7
40.1
39.7
30.2
39.0
40.0
41.1
41.4
39.3
39.0
41.4
58.4
45.2
17
27.3
84.1
86.1
32.8
50.3
44.5
72.7
32.6
39.2
40.7
41.0
39.4
37.6
39.0
43.7
39.7
61.4
18
26.3
28.7
29.0
29.4
27.6
27.8
28.6
27.1
24.5
26.8
26.1
24.8
24.5
24.4
27.5
27.5
28.0
Appendix
69
Table A2: Sequence similarities1 of DsrAB amino acid (lower triangle) and 16S rRNA gene nucleotide (upper triangle) sequences.
Reference species
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
1
Beggiatoa (uncultured)
M. magnetotacticum
M. magneticum
T. denitrificans
O. algarvensis gam 2
O. algarvensis gam 1
A. vinosum
A. ehrlichei
C. phaeobacteroides DSM266
C. limicola
P. aestuarii
P. vibrioformis
P. phaeoclathratiforme
C. phaeobacteroides BS1
C. tepidum
M. spec.
M. gryphiswaldense
R. pachyptila symbiont
Sequence similarity [ % ]
1
2
3
4
82.3
82.3
82.0
68.4
99.4
81.4
67.4
96.4
80.7
56.8
64.5
64.0
77.9
66.0
65.7
63.4
68.8
69.3
67.2
79.8
67.1
68.0
65.4
74.7
69.0
69.6
66.9
50.0
56.7
57.1
55.2
45.7
56.4
56.4
54.5
52.1
57.1
57.2
55.3
50.0
58.1
57.7
54.7
48.9
57.1
57.1
54.8
52.1
55.6
56.3
55.1
54.3
55.4
55.5
54.2
39.8
57.7
58.3
54.2
65.3
87.5
88.7
63.9
85.0
63.4
63.8
55.6
5
85.6
82.2
82.0
82.6
84.6
80.6
73.3
56.2
55.3
57.8
56.3
56.6
55.6
57.5
54.5
66.4
63.6
6
87.0
83.5
83.4
82.5
88.0
91.7
74.7
60.3
60.0
61.1
59.9
60.8
58.8
59.9
64.2
69.4
50.4
7
86.3
82.7
82.5
83.1
86.7
91.9
76.1
59.1
57.8
59.4
57.8
58.7
56.5
59.4
57.6
68.6
63.0
8
87.6
84.3
84.2
82.6
88.4
91.0
90.1
58.4
57.4
58.0
58.5
58.7
57.7
57.7
56.1
70.4
68.5
9
76.5
77.2
77.0
73.2
74.5
75.7
75.7
76.7
90.5
90.0
86.5
88.9
74.8
87.6
66.3
58.5
47.8
10
76.9
76.7
76.8
73.9
76.4
76.9
76.3
78.0
96.4
90.4
86.3
86.5
73.8
87.0
64.4
57.5
45.9
11
76.1
76.8
76.8
74.2
75.3
75.9
76.1
76.6
92.8
92.3
86.5
86.5
76.0
89.8
66.8
58.5
47.8
12
76.4
77.0
77.0
73.4
75.6
76.6
76.0
77.1
93.4
94.6
93.5
85.6
75.3
85.3
65.9
57.5
46.1
13
76.3
77.8
77.8
73.2
75.2
76.5
76.7
77.7
96.3
96.1
91.8
93.4
74.5
84.3
64.1
57.6
47.0
Similarity values were calculated using the Similarity Matrix calculation editor from the ARB software package.
14
76.7
76.8
76.7
73.3
75.3
76.3
76.4
76.9
92.6
93.0
95.5
92.4
92.3
74.6
63.3
56.5
47.0
15
76.7
76.8
76.8
72.5
75.7
76.3
75.8
77.0
92.5
94.8
93.2
95.5
92.8
92.1
65.4
57.5
49.1
16
82.1
83.2
82.8
80.3
81.4
83.5
82.1
83.8
76.4
77.5
76.7
77.6
76.6
76.8
77.6
58.4
45.2
17
82.5
96.0
95.8
81.4
82.4
83.3
83.3
84.7
77.6
77.0
76.8
77.7
77.9
76.9
77.9
83.5
61.4
18
86.2
82.4
82.2
83.6
91.6
89.7
89.0
89.8
75.3
77.7
76.7
76.3
76.8
75.9
76.3
82.7
82.5
70
(E) Acknowledgement
I would like to thank my thesis supervisor Dr. Marc Mußmann for providing
such an exciting research topic and for his encouragement. His fascination for sulfur
oxidizers greatly supported this work. He always found time for helpful discussions.
Many thanks go to Prof. Dr. Rudolf Amann for setting up an excellent research
framework and for taking the part of the first reviewer.
Dr. Jens Harder is acknowledged for suggesting me for the marmic program and
thus offering me the possibility to participate in this outstanding Master degree of the
International Max Planck Research School for Marine Microbiology.
Moreover, Dr. Christiane Glöckner is very much acknowledged for caring
management of the entire marmic class.
Finally, I want to thank all members of the Molecular Ecology group for the
support and kind working atmosphere.
This work was financially supported by an educational grant of the Max Planck
Society.
Allen voran besonderen Dank meiner Familie und Falk☺ für die Unterstützung
und tatkräftiges zur Seite stehen in den gesamten Studienjahren.
71
(F) Statement
I herewith confirm that I have written this thesis unaided and that I used no other
resources than those mentioned.
Bremen 17.03.2006
Sabine Lenk