Microbial Diversity and Biochemical Potential Encoded by Thermal Spring Metagenomes Derived from the Kamchatka Peninsula

Volcanic regions contain a variety of environments suitable for extremophiles. This study was focused on assessing and exploiting the prokaryotic diversity of two microbial communities derived from different Kamchatkian thermal springs by metagenomic approaches. Samples were taken from a thermoacidophilic spring near the Mutnovsky Volcano and from a thermophilic spring in the Uzon Caldera. Environmental DNA for metagenomic analysis was isolated from collected sediment samples by direct cell lysis. The prokaryotic community composition was examined by analysis of archaeal and bacterial 16S rRNA genes. A total number of 1235 16S rRNA gene sequences were obtained and used for taxonomic classification. Most abundant in the samples were members of Thaumarchaeota, Thermotogae, and Proteobacteria. The Mutnovsky hot spring was dominated by the Terrestrial Hot Spring Group, Kosmotoga, and Acidithiobacillus. The Uzon Caldera was dominated by uncultured members of the Miscellaneous Crenarchaeotic Group and Enterobacteriaceae. The remaining 16S rRNA gene sequences belonged to the Aquificae, Dictyoglomi, Euryarchaeota, Korarchaeota, Thermodesulfobacteria, Firmicutes, and some potential new phyla. In addition, the recovered DNA was used for generation of metagenomic libraries, which were subsequently mined for genes encoding lipolytic and proteolytic enzymes. Three novel genes conferring lipolytic and one gene conferring proteolytic activity were identified.


Introduction
Sites of volcanic activity can be found all over the world and even under the sea. Volcanic regions provide a variety of different environments for extremophilic archaeal and bacterial microorganisms. Well-known examples of such extreme environments are terrestrial surface hot springs. With respect to geographical, physical, environmental, and chemical characteristics, hot springs are unique sites for extremophilic microorganisms [1][2][3]. Extremophiles inhabiting hot springs are considered to be the closest living descendants of the earliest life forms on Earth [4,5]. Therefore, these springs provide insights into the origin and evolution of life. In addition, thermophiles and hyperthermophiles produce a variety of hydrolytic enzymes such as lipases, glycosidases, peptidases and other biomolecules, which are of industrial interest [6][7][8]. For example, Hotta et al. [9] found an extremely stable carboxylesterase in the hyperthermophilic archaeon Pyrobaculum calidifontis VA1, and Arpigny et al. [10] identified a novel heat-stable lipolytic enzyme in Sulfolobus acidocaldarius DSM 639.
Especially in extreme environments, most microorganisms are reluctant to cultivation-based approaches [11,12]. Therefore, culture-independent metagenomic strategies are promising approaches to assess the phylogenetic composition and functional potential of microbial communities living in extreme environments [7,13,14]. For example, Simon et al. studied the prokaryotic community in glacier ice and found a 2 Archaea highly diverse bacterial community [15]. In 1998, Hugenholtz et al. [1] investigated the bacterial diversity in the Obsidian Pool in Yellowstone National Park and identified several new bacterial candidate divisions. The same pool and two others were studied later by Meyer-Dombard et al. [16]. They encountered diverse bacterial and archaeal communities in all three hot springs.
In the present study, we investigated the phylogenetic composition and metabolic potential of two microbial communities derived from two extreme sites of the Kamchatka peninsula, which is located in the Far East of Russia. The Kamchatka peninsula comprises an area of approximately 472,300 km 2 and is described as the land of fire by its first explorers due to the high density of volcanoes and associated volcanic phenomena. For example, the largest active volcano of the northern hemisphere, the Klyuchevskaya Sopka, is located on the Kamchatka Peninsula. Sediment samples analyzed in this study were taken from two hot springs providing a thermoacidophilic (70 ∘ C, pH 3.5-4) or a thermophilic (81 ∘ C, pH 7.2-7.4) environment. The composition of the prokaryotic communities of the two Kamchatkian hot springs was assessed by 16S rRNA gene analysis. In addition, metagenomic libraries were generated and screened for novel biocatalysts.

Sampling and DNA Extraction.
Two sediment samples were taken from the hot springs located on the Kamchatka peninsula in summer 2001. The first sample was collected from a thermoacidophilic spring (70 ∘ C, pH 3.5-4) at the Mutnovsky volcano (52.453 N, 158.195 E). The second sample was taken from a thermophilic spring (81 ∘ C, pH 7.2-7.4) in the Uzon Caldera (54.5 N, 159.967 E). The chemical analysis of both sediment samples is shown in Table 1.
DNA was extracted as described by Zhou et al., 1996 [17]. The concentration of the recovered DNA was quantified using a NanoDrop ND-1000 spectrophotometer (PEQLAB, Erlangen, Germany).

Amplification of 16S rRNA Genes and Generation Clone
Libraries. To assess the prokaryotic community structure, archaeal and bacterial 16S rRNA genes were amplified by PCR and analyzed. The PCR reaction mixture (50 L) contained 2.5 L of 10-fold Mg-free Taq polymerase buffer, 200 M of each of the four deoxynucleoside triphosphates, 1.75 mM MgCl 2 , 0.4 M of each primer, 1 U of Taq DNA polymerase (Fermentas, St. Leon-Rot, Germany), and approximately 25 ng of recovered DNA as a template. Prokaryotic 16S rRNA genes were amplified with the following set of primers: 8F 5 -AGAGTTTGATCMTGGC-3 [18] and 1114R 5 -GGG-TTGCGCTCGTTRC-3 [19], A800F 5 -GTAGTCCYGGCY-GTAAAC-3 [20] and A1530R 5 -GGAGGTGATCCAGCC-G-3 [21], and Arch8F 5 -TCCGGTTGATCCTGCCGG-3 [15] and Arch958R 5 -YCCGGCGTTGAMTCCAATT-3 [22]. The following thermal cycling scheme was used: initial denaturation at 94 ∘ C for 2 min, 25 cycles of denaturation at 94 ∘ C for 1.5 min, annealing at 56 ∘ C (8F and 1114R), 51 ∘ C (A800F and A1530R), or 55 ∘ C (Arch8F and 958R), followed by extension at 72 ∘ C (1 min for 1 kb). The final extension was carried out at 72 ∘ C for 10 min. Negative controls were performed by using the reaction mixture without template. The obtained PCR products were purified using the Wizard SV Gel and PCR clean up system (Promega, Madison, USA) and subsequently cloned into pCR2.1-TOPO as recommended by the manufacturer (Invitrogen, Carlsbad, USA). The resulting recombinant plasmids were used to transform Escherichia coli TOP 10 cells. A total of 1271 insert-carrying plasmids were isolated from randomly selected E. coli clones. The insert sequences were determined by the Göttingen Genomics Laboratory (Göttingen, Germany).

Analysis of 16S rRNA Genes.
To assess the prokaryotic community structure, the retrieved 16S rRNA gene sequences were analyzed using QIIME [23]. The obtained 16S gene sequences were edited using gap4 [24] and initially checked for the presence of chimeric sequences using Mallard [25], Bellerophon [26], and Chimera Check [27]. Remaining sequences were clustered employing the UCLUST algorithm [28] and the following QIIME scripts: pick otus.py and pick rep set.py. The sequences were clustered in operational taxonomic units (OTUs) at 1, 3, and 20% genetic dissimilarity. The phylogenetic composition of the prokaryotic communities in both samples was determined using the QIIME assign taxonomy.py script. A BLAST alignment [29] against the most recent SILVA ARB database [30] was performed. Sequences were classified with respect to the taxonomy of their best hit in the ARB database. Finally, OTU tables were generated. Rarefaction curves, Shannon indices [31], and Chao1 indices [32] were calculated employing QIIME. In addition, the maximal number of OTUs ( max ) was estimated for each sample using the Michaelis-Menten-fit alpha diversity metrics included in the QIIME software package.
One sequence per OTU (1% genetic distance) was further used for the construction of phylogenetic trees. Sequences were imported into the most recent SSU Ref SILVA database of the ARB program package [33]. Multiple sequence alignments were checked manually and improved by employing the ARB editor tool. Phylogenetic trees were created by employing the maximum parsimony algorithm implemented in ARB. The robustness of obtained tree topologies was evaluated by bootstrap analysis with 100 resemblings.

Construction of Small-Insert Metagenomic Libraries.
To exploit the biochemical potential, metagenomic small-insert libraries were generated. Due to the low DNA recovery, starting material for the generation of these libraries was obtained by multiple displacement amplification employing the GenomiPhi V2 DNA Amplification Kit (GE Healthcare, Munich, Germany). To improve cloning efficiency, hyperbranched structures were resolved and the DNA was inserted into pCR-XL-TOPO (Invitrogen) as described by Simon et al. [34]. In this way, two metagenomic libraries were generated.

Screening for Hydrolytic Activity and Identification of
Corresponding Genes. The constructed metagenomic libraries were screened for genes conferring lipolytic or proteolytic activity using a function-driven approach. The constructed libraries were used to transform E. coli DH5 cells. Recombinant cells were plated on LB agar plates containing either 1% tributyrin (lipolytic activity) [35] or 2% skim milk (proteolytic activity) [36]. Plates were incubated at 37 ∘ C for up to two weeks. Hydrolytic activity is indicated by halo formation. To determine the substrate specificity of the protease-producing clones, skim milk was replaced by 0.3% (w/v) azocasein, azoalbumin, or elastin-Congo red. Insert sequences of recombinant plasmids derived from positive clones were determined by the Göttingen Genomics Laboratory. The retrieved insert sequences were edited employing gap4 [24], and putative ORFs were annotated using Artemis (version 11.0) [37].

Cloning of Genes Conferring Lipolytic or Proteolytic Activity into Expression Vectors and Purification of the
Corresponding Gene Products. For enzyme production of lipolytic proteins, identified genes were cloned into pET101-TOPO according to the Champion pET101 directional TOPO Expression Kit (Invitrogen). In this way, sequences encoding a His 6 tag and a V5 epitope provided by the vector were added to the 3 end of the coding regions. Alternative start codons were replaced by ATG. Escherichia coli BL21 Star (DE3) (Invitrogen) was used as a host for enzyme production. The production was performed as recommended by the manufacturer. Subsequently, recombinant cells were harvested by centrifugation, washed with Tris buffer (30 mM, pH 8.0), and resuspended in 50 mmol L −1 of sodium phosphate buffer containing 0,3 mol L −1 NaCl (pH 8.0). The cells were disrupted employing a French Press (1.38 × 10 8 Pa). Subsequently, the extract was cleared by centrifugation at 14,000 g and 4 ∘ C for 45 min. The supernatant was used as a source for soluble proteins. Recombinant proteins were purified using the Protino Ni-TED 2000 packed columns (Macherey-Nagel, Düren, Germany) as recommended by the manufacturer. Protein preparations were dialyzed using sodium phosphate buffer (50 mM, pH 7.5) to remove residual imidazole. Protein concentration in the purified sample was determined with Roti-Quant (Carl Roth, Karlsruhe, Germany) as suggested by the manufacturer. The purity of the protein preparations was analyzed by SDS-polyacrylamide gel electrophoresis according to Laemmli [38].
For protein production of the putative protease, the corresponding gene (pepBW1) was cloned into pBAD/Myc-His (Invitrogen) using a modified fusion method [39]. The gene was amplified in two PCR reactions with the following two sets of primer pairs containing synthetic sites (underlined) for cloning into the vector: pair 1, 5 -CATGGTGTTCAAT-AAATATGTCTT-3 and 5 -CTAGGGTAGACTTAACGC-3 and pair 2, 5 -GTGTTCAATAAATATGTCTTATT-3 and 5 -CGCTAGGGTAGACTTAACGC-3 . After mixing, denaturation, and hybridization, four different hybridization products are formed of which one contained overhangs complementary to vector digested with NcoI (Fermentas) and Bsp119I (Fermentas). The preligation mixture (10 L) contained 1 L O-buffer (Fermentas) and approximately 50 ng of each PCR product. To create the appropriate overhangs for cloning, the following thermal cycling scheme was used: denaturation at 95 ∘ C for 3 min, 4 cycles of reannealing at 65 ∘ C for 2 min, and reannealing at 25 ∘ C for 15 min. The ligation reaction mixture (20 L) contained the preligation mixture, 1 L of O-Buffer (Fermentas), 1 L ATP (10 mM), 1 U of T4 DNA ligase, and 10 ng of pBAD/Myc-His digested with NcoI and Bsp119I. The reaction was incubated at 16 ∘ C overnight and inactivated by heating at 65 ∘ C for 10 min. The resulting recombinant plasmids were then used to transform E. coli top 10 cells.

Characterization of Lipolytic
Activity. Determination of enzyme specificity against different triacylglycerides was determined by growing E. coli BL21 (DE3) harboring the recombinant plasmids on LB agar plates containing triacylglycerides with different chain length (C4 to C18). Plates were supplemented with IPTG to a final concentration of 0.1 mM to induce gene expression.

Archaea
For quantitative analysis, p-nitrophenyl esters with various chain lengths were used as descripted by Rashamuse et al. [40]. Routine esterase activity assays were performed by measuring the release of p-nitrophenol from a p-nitrophenyl (p-NP) ester at 410 nm using a Cary 100 UV-Vis spectrophotometer (Varian, Palo Alto, USA) with a Peltier temperature controller. Unless otherwise described, enzyme activity was measured at 50 ∘ C in Tris-HCl (50 mmol L −1 ; pH 7.5) with p-NP caprylate (1 mmol L −1 ; dissolved in 2-propanol) as a substrate. p-NP caprylate was used as a substrate in the standard assay because of its stability at high temperatures and alkaline pH values. Enzyme activity of EstBW2 was determined with p-NP butyrate as the enzyme showed no activity with the other substrates tested. Therefore, temperature and pH dependence of EstBW2 were not measured above 75 ∘ C and pH 8. All measurements were performed in triplicate.
To determine substrate specificity, enzyme activity was measured at standard assay conditions employing the following p-NP esters of various chain lengths: The temperature dependence of enzyme activity was determined between 20 and 95 ∘ C under standard assay conditions. To compensate temperature effects on pH values, buffers were preheated to set-point temperature and adjusted using Tris buffer (50 mmol L −1 ). Thermostability was measured by incubating the enzyme at different temperatures over various time periods. Enzyme activity was subsequently measured under standard assay conditions. Optimal pH values for enzyme activity were measured under standard assay conditions employing different overlapping buffer solutions (50 mM): sodium acetate buffer (pH 4 and 5), sodium phosphate buffer (pH 5, 6, and 7), Tris-HCl (pH 7, 8, and 9), CHES (pH 9 and 10), and CAPS (pH 10 and 11).
The effect of different detergents on enzyme activity was determined by under standard assay conditions in the presence of 1 mM AgNO 3 , 1 mM CaCl 2 , 1 mM CoCl 2 , 1 mM cetyltrimethylammonium bromide (CTAB), 1 mM CuCl 2 , 1 mM ethylenediaminetetraacetic acid (EDTA), 1 mM FeCl 3 , 1 mM KCl, 1 mM MgCl 2 , 1 mM MnCl 2 , 1 mM NaCl, 1 mM NiSO 4 , 1 mM Sodium Dodecyl Sulfate (SDS), 1 mM ZnCl 2 , 0.01% (v/v) Tween 80, or 0.01% (v/v) 2-Mercaptoethanol. The serine dependence of the recovered lipolytic enzymes was validated under standard assay conditions by incubation in the presence of 1 mM phenylmethylsulfonyl fluoride (PMSF). In addition, we analyzed the significance of the determined effects on enzyme activity. As the enzyme tests were performed in triplicate, we assumed that all measured enzyme activities were normally distributed. The variance homogeneity was tested employing the F-test, and the significance of the detergent effects was subsequently tested either with the Student's t-test (homogenous variances) or the Welch's t test (heterogeneous variances). All statistical analyses were performed in R [41].

Sampling and Chemical Properties of the Investigated
Sediments. Sediment samples were collected from two Kamchatkian hot springs. The springs were located near the Mutnovsky volcano (Mutnovsky sample) and in the Uzon Caldera (Uzon sample), which represent a thermoacidophilic (70 ∘ C, pH 3.5-4) and a thermophilic (81 ∘ C, pH 7.2-7.4) environment, respectively. Both investigated sediments were chemically distinct from each other ( Table 1). The Mutnovsky sample contained higher Al, Ca, Co, Cu, Fe, Pb, and Zn concentrations than the Uzon sample. For As, B, Ba, K, Mn, and Na concentrations, the opposite was recorded. The concentrations of Cr, Mg, Ni, Sr, P, Ti, and V and the total organic carbon contents were almost identical in both samples.

Isolation of Metagenomic DNA and Construction of Metagenomic Libraries.
To assess the prokaryotic diversity and metabolic potential by metagenomic approaches, environmental DNA was extracted from both samples. Approximately, 2.7 g DNA per 10 g sediment was recovered from both samples. After removal of remaining salts, archaeal and bacterial 16S rRNA genes were amplified from the purified DNA. The resulting PCR products were used for the generation of 16S rRNA gene libraries. A total of 1271 clones were sequenced from these libraries. After quality filtering and removal of potential chimeric sequences, 1235 high-quality 16S rRNA gene sequences were obtained (536 for the Mutnovsky sample, 699 for the Uzon sample). The DNA from both samples was also used to construct metagenomic libraries. The Mutnovsky library comprised approximately 479,000 plasmids with an average insert size of 5.3 kb. The percentage of insert-carrying plasmids was 74%. The Uzon library consisted of approximately 117,000 plasmids with an average insert size of 4 kb. The percentage of insert-carrying plasmids was 85%. In summary, the generated small-insert metagenomic libraries harbored approximately 2.27 Gbp of cloned environmental DNA.

Archaeal Community Structures.
We were able to assign 265 16S rRNA gene sequences of both samples to the domain Archaea. The classified sequences were affiliated to four different archaeal phyla (Figure 1). The Thaumarchaeota was the most abundant archaeal phylum in both samples (57% and 68% of all sequences, resp.). Most of the sequences were affiliated to Miscellaneous Crenarchaeotic Group nowadays belonging to the recently proposed Thaumarchaeota (37%). Another abundant thaumarchaeotic group was the Terrestrial Hot Spring Group (24.7%). The majority of the remaining sequences of the Mutnovsky sample were affiliated to uncultured members of the Euryarchaeota (34.7%), which were only detected in this sample. The majority of the remaining Uzon sequences belonged to the Crenarchaeota (32.1%). Sequences were affiliated to known genera such as Sulfophobococcus (12.1%), Thermofilum (6.8%), Ignisphaera (5.3%), and to Desulfurococcus kamchatkensis (4.2%). The archaeal phylum Korarchaeota was only identified in the Uzon sample (1 sequence).

Diversity and Species Richness of Archaeal Communities.
To determine the archaeal diversity and richness, rarefaction analyses were performed with QIIME. The observed OTU numbers in the Mutnovsky sample and the Uzon sample were 33 and 13 (1% genetic distance), 25 and 11 (3% genetic distance), and 7 and 5 (20% genetic distance), respectively ( Table 2). The maximal expectable number of clusters for both samples was determined based on the Michaelis-Menten fit metrics. On average, more than 90% of the entire archaeal community was covered by the surveying effort. Shannon indices of the Mutnovsky and Uzon sample were 1.83 and 2.96 (1% genetic distance), 2.70 and 1.83 (3% genetic distance), and 0.86 and 1.87 (20% genetic distance), respectively. These indicated low archaeal diversity in the investigated samples.
Comparison of the rarefaction analyses with the number of OTUs determined by Chao1 richness estimator revealed that, at 1 and 3% genetic distance, the rarefaction curves were almost saturated ( Figure 4). Thus, the majority of the estimated richness was recovered by the surveying effort ( Table 2).

Diversity and Species Richness of Bacterial Communities.
The observed OTU numbers in both hot springs were 17 and 50 (1% genetic distance), 12 and 42 (3% genetic distance), and 10 and 11 (20% genetic distance) in the Mutnovsky sample and the Uzon sample, respectively ( Table 2). Analysis of the maximal expectable number of clusters indicated that more than 94% of the entire bacterial community was recovered by the surveying effort. Correspondingly, comparison of the rarefaction analyses with the number of OTUs determined by Chao1 richness estimator revealed that at 1%, 3%, and 20% genetic distance, the rarefaction curves were almost saturated (Table 2, Figure 4).

Screening of Metagenomic
Libraries. The two generated metagenomic libraries were employed in a function-based screening to identify novel lipolytic and proteolytic enzymes. Three novel genes encoding lipolytic enzymes (plpBW1, estBW1, and estBW2) and one gene encoding a proteolytic enzyme (pepBW1) were identified during the screening of the metagenomic library derived from the Uzon sample. No hydrolytic enzymes were identified within the metagenomic library derived from the Mutnovsky sample.
The closest relatives of all identified protein sequences originated from known thermophiles. They were similar to uncharacterized putative gene products derived from Desulfurococcus kamchatkensis (PepBW1), Sulfurihydrogenibium azorense (PlpBW1 and EstBW2), and Thermobaculum terrenum (EstBW1) (  PlpBW1 was affiliated to the patatin-like proteins (PLPs). Four conserved domains are described for this enzyme type [42], which could all be identified within the amino acid sequence (data not shown). Interestingly, PLPs do not possess a catalytic triad. The lipolytic activity is conferred by a catalytic dyad formed by a serine residue and an aspartate residue [42]. EstBW1 and EstBW2 were affiliated to family V of lipolytic enzymes according to the classification system of Arpigny and Jaeger [43]. PepBW1 was classified employing the MEROPS database [44]. It was affiliated to the subtilisin family (family S8). Interestingly, the pepBW1 gene sequence was almost identical to that of a putative gene encoding a serine peptidase of Desulfurococcus kamchatkiensis (Table 3); Desulfurococcus kamchatkiensis belongs to the Crenarchaeota and was also isolated from a thermal spring within the Uzon Caldera [45]. In addition, the 16S rRNA gene sequence of this species was found in our 16S analysis of the Uzon sample.

Characterization of Recombinant Enzymes.
To characterize all recombinant proteins, the genes conferring lipolytic and proteolytic activity were cloned into expression vectors.
The recombinant E. coli strain containing PepBW1 was tested towards different proteins and showed proteolytic activity 8 Archaea Table 2: Prokaryotic diversity and richness at 1, 3, and 20% genetic distance. Numbers of observed OTUs as well as Shannon and Chao1 indices were calculated with QIIME [16]. The maximal number of OTUs ( max ) was calculated using the Michaelis-Menten-fit diversity metrics implemented in the QIIME package. Coverage was determined by dividing the observed number of OTUs with max .  with skim milk and elastin-Congo red but not with azoalbumin or azocasein. The activities of the recombinant lipolytic proteins were tested towards different triacylglycerides. All proteins showed activity with tributyrin as substrate. In addition, PlpBW1 showed activity with long-chain triacylglycerides, up to trimyristin (C14). Hydrolysis of different p-nitrophenyl esters was used to further analyze the substrate specificity (Figure 5(a)). PlpBW1 and EstBW1 showed highest activity with p-NP acetate and p-NP butyrate, respectively. Both enzymes exhibit activity towards all tested p-NP esters, except p-NP palmitate. The activity decreased with increasing chain length. In contrast, EstBW2 showed only activity towards p-NP butyrate. Specific activities under standard assay conditions using the optimal substrate were 2.6 ± 0.3 U/mg (PlpBW1), 2.33 ± 0.32 U/mg (EstBW1), and 1.89 ± 0.21 U/mg (EstBW2). Based on the results, all three lipolytic enzymes are most likely carboxylesterases and not lipases.
All lipolytic enzymes were active over a wide temperature range. PlpBW1, EstBW1, and EstBW2 retained a minimum of 50% activity from 60 to 90 ∘ C, 65 to 95 ∘ C, and 40 to 75 ∘ C, respectively ( Figure 5(b)). Maximal activities were recorded for PlpBW1 at 85 ∘ C, for EstBW1 at 90 ∘ C, and for EstBW2 at 65 ∘ C. We further determined the stability of the three lipolytic enzymes with respect to different temperatures. The half-lives of PlpBW1 were 45 min at 70 ∘ C, 15 min at 80 ∘ C, and 5 min at 90 ∘ C. EstBW1 exhibited half-lives of 5 h at 70 ∘ C, 2.5 h at 80 ∘ C, and a remarkable half-life of 15 min at 90 ∘ C. EstBW2 was less stable at 90 ∘ C (7 min half-live), but the activity was almost unaffected by 5 h of incubation at 70 ∘ C and 80 ∘ C (data not shown).
The pH effect on enzyme activity was measured at pH values ranging from 4 to 11 ( Figure 5(c)). All enzymes exhibited high activity at neutral or alkaline pH values. Maximal activities were determined at pH 10 (PlpBW1) and pH 7 (EstBW1 and EstBW2).
Addition of EDTA, KCl, or NaCl to the reaction mixture had no significant effect on enzyme activity ( > 0.05), whereas all other tested detergents exhibited an effect on the activity of at least one of the recovered enzymes (Table 4). CTAB, Tween 80, and ZnCl 2 impacted the activities of all three enzymes significantly ( < 0.05). PlpBW1 showed a more than 2.5-fold higher activity in presence of CTAB, whereas EstBW1 and EstBW2 displayed a loss in activity. Tween 80 increased enzyme activity of PlpBW1 and EstBW2 but not that of EstBW1. The addition of ZnCl 2 decreased the activity of all three recombinant enzymes. Lipolytic activity of EstBW1 and EstBW2 was completely inhibited by the phenylmethylsulfonyl fluoride, indicating the presence of a serine residue at the active site of both enzymes. Interestingly, the activity of PlpBW1 was not affected by PMSF.

Prokaryotic Community Composition in the Kamchatkian
Springs. The number of metagenomic studies has been rapidly increased over the past years. Metagenomics has been Curves were calculated at 1%, 3%, and 20% genetic distance level employing QIIME [23].
employed to assess and exploit the biodiversity of many habitats including environments of extremophiles [1,15,16,46,47]. In this study, we investigated the prokaryotic diversity of two hot springs located on the Kamchatka peninsula. We found different bacterial and archaeal communities at both sites, which were dominated by Proteobacteria, Thermotogae, and Thaumarchaeota. Jackson et al. (2001) studied a mat derived from the Norris Geyser Basin, an acidic thermal spring in the Yellowstone National Park [46]. They found community pattern comparable to that in the Mutnovsky sample with one difference. They were not able to identify Thaumarchaeota, which is not surprising as this phylum was first proposed in 2008 [48]. Members of this phylum are not restricted to thermophilic habitats as they were originally described as mesophilic Crenarchaeota [48][49][50]. A study by Meyer-Dombard et al. (2005) investigated the prokaryotic community in three thermal springs in the Yellowstone National Park (the Silvan Spring, the Bison Pool, and the Obsidian Pool) [16]. Whereas the other pools have a rather neutral milieu, the Silvan Spring has a low pH of 5. However, the prokaryotic community structure of this acidic spring was different to that found in the acidic Mutnovsky spring sample. Meyer-Dombard et al. identified the Crenarchaeota as the most abundant archaeal group, whereas Thaumarchaeota were the most abundant group in our study. A more recent study on prokaryotic community composition of hot springs on the Tibetan Plateau also found Thaumarchaeota as the dominant archaeal group [48].
Analysis of the Uzon sample revealed a more diverse prokaryotic community than in the Mutnovsky sample. Only two OTUs at 1% genetic distance were shared, whereas all the other OTUs were unique for each sample (Figures 2  and 3). The observed differences in community composition between the two sampling sites might be due to the different Relative activity (%) 25  temperatures and pH values at the sites. However, Huang et al. (2011) found no statistical correlation between temperature and diversity [48]. Despite the geographical separation, the Obsidian Pool and the Uzon Caldera hot spring share a very similar community structure, as almost the same dominant archaeal and bacterial groups were identified [16]. In addition, some rare phyla were present in both samples, that is, the Korarchaeota. This phylum is a relatively new phylum first described by Barns et al. in 1996 [51]. Another rare bacterial group found in both samples was the Candidate division OP9. These results confirm also the presumption proposed for other extreme environments that similar environmental conditions result in similar microbial communities [15].

Hydrolytic Enzymes.
In the present study, we were able to identify three novel lipolytic enzymes and one proteolytic enzyme. The determined optimal temperatures and pH values reflect the environmental conditions of the samples used for DNA isolation, indicating that the environment shapes the characteristics of the enzymes. Correspondingly, the characterized lipolytic enzymes (PlpBW1, EstBW1, and EstBW2) showed features similar to those of other metagenomederived esterases, which were identified in thermophilic sites. Rhee et al. (2005) identified a thermophilic esterase in metagenomic libraries generated from hot spring and mud hole DNA [52]. The enzyme was active from 30 to 95 ∘ C and exhibited an optimal pH value of approximately 6.0.  derived from a Thailand hot spring and identified two novel lipolytic enzymes, of which one was also characterized as a patatin-like protein [53]. To our knowledge, PlpBW1 is the second reported patatin-like protein derived from a hot spring metagenomic library up to now [53]. In contrast to most other lipolytic enzymes containing a serine residue in the active site, PlpBW1 is not inactivated by the inhibitor PMSF [51,54]. The effect of Zn 2+ ions recorded for all recombinant enzymes investigated in this study was also mentioned for esterases studied by Chu et al. [54]. They also recorded a decrease of activity in presence of Zn 2+ ions. The activity of the recovered lipolytic enzymes was positively or negatively influenced by addition of CTAB or Tween 80 (Table 4). It has been shown that these detergents can either promote or decrease activity of lipolytic enzymes by formation of micellar aggregates and monomers which then interact with hydrophobic parts of the enzymes [55].
In addition to the esterases, metagenomic libraries were mined for proteolytic activity. The identified serine peptidase, PepBW1, is the first metagenome-derived peptidase from a thermophilic environment. As PepBW1 is derived from an archaeal organism, it illustrates that screening in a heterologous host can be successful, even if the target gene originates from a different domain of life [56].

Ecology of Hot Springs.
As most studies on ecology of hot springs are targeting the prokaryotic diversity, for example, via 16S rRNA gene or other marker gene analyses, little is known on the global relevance of these extremophilic communities. Burgess et al. (2011) studied two thermal pools in the Uzon Caldera by 16S rRNA gene analysis and related some community members to different archaeal and bacterial groups, which might play a role in cycling of C, N, and S [57].
The vital role of Archaea in N 2 fixation and denitrification is well established [58]. The first step of nitrification, ammonium oxidation, was originally thought to be restricted to some Proteobacteria [59]. However, recent metagenomic studies provided evidence that Archaea are capable to oxidize ammonium to nitrate [58,59]. Until recently, methanogenic Euryarchaeota were thought to be the only archaeal group of global relevance for element cycling. This presumption changed with the discovery of ammonia-oxidizing archaea [60], which are affiliated to the recently proposed phylum Thaumarchaeota. Members of this phylum contribute significantly to the global N cycle, as their high abundance and extremely low substrate threshold provides compelling evidence for a dominant role as ammonia oxidizers in open oceans [60]. In our study, we identified diverse thaumarchael groups in both investigated sediment samples. Thus, hot springs may also play a major role in the global N cycle.