A metagenomic-based survey of microbial (de)halogenation potential in a German forest soil

In soils halogens (fluorine, chlorine, bromine, iodine) are cycled through the transformation of inorganic halides into organohalogen compounds and vice versa. There is evidence that these reactions are microbially driven but the key enzymes and groups of microorganisms involved are largely unknown. Our aim was to uncover the diversity, abundance and distribution of genes encoding for halogenating and dehalogenating enzymes in a German forest soil by shotgun metagenomic sequencing. Metagenomic libraries of three soil horizons revealed the presence of genera known to be involved in halogenation and dehalogenation processes such as Bradyrhizobium or Pseudomonas. We detected a so far unknown diversity of genes encoding for (de)halogenating enzymes in the soil metagenome including specific and unspecific halogenases as well as metabolic and cometabolic dehalogenases. Genes for non-heme, no-metal chloroperoxidases and haloalkane dehalogenases were the most abundant halogenase and dehalogenase genes, respectively. The high diversity and abundance of (de)halogenating enzymes suggests a strong microbial contribution to natural halogen cycling. This was also confirmed in microcosm experiments in which we quantified the biotic formation of chloroform and bromoform. Knowledge on microorganisms and genes that catalyze (de)halogenation reactions is critical because they are highly relevant to industrial biotechnologies and bioremediation applications.

by enzymes. A major group of halogenating enzymes are the haloperoxidases which unspecifically halogenate organic matter using hydrogen peroxide and a halogen ion (Cl − , Br − , I − ) as substrate [34][35][36][37] . Based on their cofactors they can be classified into heme-dependent haloperoxidases 38 and vanadium-dependent haloperoxidases 39 . Perhydrolases, or non-heme, no-metal haloperoxidases also require hydrogen peroxide and catalyze unspecific halogenation reactions but do not contain any metal cofactors 36 . Beside the haloperoxidases also halogenases with specific and regioselective halogenation reaction mechanisms exist. Flavin-dependent halogenases are involved in bacterial secondary metabolism, e.g. antibiotic syntheses 40 . Another class of specific halogenases are the alpha-ketoglutarate-dependent halogenases 41 . One known halogenase, a bacterial fluorinase, is able to fluorinate S-adenosyl-L-methionine via a nucleophilic mechanism 42 . Furthermore, methyltransferases of plants, fungi and algae 43 are known to form halomethanes. Since organohalogen compounds are prominent environmental pollutants, their biotic degradation has been studied intensely in the past decades and a variety of different dehalogenation pathways including hydrolytic dehalogenation, dehydrohalogenation, thiolytic dehalogenation and intramolecular substitution have been described 36,44 . Dehalogenation of halomethanes by methyltransfer was described for bacterial methyltransferases 45,46 . Microorganisms can use organohalogens either as carbon source (metabolic degradation) 31 or they are co-metabolically degraded during the degradation of primary substrates such as methane 47 . Metabolic and cometabolic degradation of organohalogens are possible under oxic and anoxic conditions 31,48 . Organohalogens can even be used as terminal electron acceptor in a metabolism called organohalide respiration 49 . Numerous pathways and enzymes involved in biotic halogenation and dehalogenation reactions have been identified. But so far little is known about the natural diversity and abundance of the different groups of halogenating and dehalogenating enzymes. It is further not well understood which genes and microorganisms are the main contributors to biotic halogen cycling 27,28 . Natural halogenation in soils is widespread and not only restricted to forest soils. It also occurs in grasslands and agricultural soils and the microbial chlorination and dechlorination of soil organic matter seems to be an ubiquitous phenomenon 50 . Knowledge on the microbial potential for halogenation and dehalogenation reactions in soils is important, since soils act as important sources of volatile organohalogens (e.g. CHCl 3 51 ), as well as sinks for natural and anthropogenic organohalogen compounds 32 . Here we combined geochemical analyses with microcosm experiments and shotgun metagenomics to unravel the natural diversity and relative abundance of genes encoding for halogenating and dehalogenating enzymes in a forest soil.
At the sampling site two replicate soil profiles were sampled within a distance of 2 m from each other. Bulk soil samples for each profile were collected from the three distinguishable horizons of the top 40 cm, homogenized and stored at − 80 °C for genetic analysis. For biogeochemical analysis bulk samples of the two soil profiles were mixed, homogenized and stored at 4 °C. Samples were taken in October 2013.
Geochemical analysis. For water content determination, fresh soil samples were weighed and subsequently dried at 105 °C until weight stability. pH was measured in a suspension of 10 g air dried soil in 25 mL of a 0.01 M CaCl 2 -solution. For determination of leachable chloride and leachable organic carbon, 10 g of soil were mixed with 100 mL deionized water and shaken at 150 rpm for 24 h on a rotary shaker. Samples were centrifuged for 5 minutes at 4000 × g and filtered through a 0.45 μ m pore size cellulose ester filter (Millex HA filter, EMD Millipore Corporation, USA). Dissolved organic carbon was measured with a High TOC Elementar system (Elementar Analysensysteme GmbH, Hanau, Germany) and chloride was determined by ion chromatography (Dionex DX 120, Thermo Scientific, Sunnyvale, CA, USA). For total organic carbon analysis soil samples were dried at 40 °C and sieved (2 mm mesh) to exclude large roots and stones. The organic carbon content was determined by heat combustion (1150 °C) and thermal conductivity analysis on a CNS element analyzer (Elementar Vario EL III, Elementar Analysensysteme GmbH, Hanau, Germany). Adsorbable organic halogen (AOX) content in the soil samples was determined according to the standard protocol (DIN EN ISO 9562) for soil leachates (DIN EN 12457-4) at the Laboratory for Environmental and Product Analytics (DEKRA GmbH, Halle, Germany).

Detection of volatile organohalogen compounds.
Microcosm experiments to quantify formation of volatile organohalogen compounds (VOX) in the soil horizons via GC-MS were set up in triplicates per soil horizon as follows: 3.5 g of native soil were incubated with 8.5 mL of sterile deionized water and incubated for 1 h at 30 °C in the dark prior to VOX quantification. Details on incubation conditions and GC-MS measurements have been published previously 15 . DNA extraction. Three different methods were applied to extract genomic DNA from the replicate soil samples. 10 g of soil were used for extraction with the PowerMax ® Soil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, USA). Furthermore we applied a microwave-based extraction method 54 with the following modifications: all steps were up-scaled for the extraction of DNA from 6 g of soil and DNA was precipitated by mixing the supernatant from the chloroform-isoamylalcohol-extraction with an equal amount of isopropanol followed by a 1 h incubation step at room temperature. The third DNA extraction protocol was based on a sodium-dodecyl-sulfate method combined with freeze-thawing, protein digestion and chloroform-isoamylalcohol extraction 55 . Since the DNA extracts of the latter two methods were still of brownish color, DNA was further purified by agarose gel electrophoresis using 0.7% agarose gels. High molecular weight Quality processing, sequence alignment, taxonomic and functional analysis. Quality processing was performed using the Metagenomic RAST server 56 . Quality processing included trimming of low quality bases with the SolexaQA software package 57 and a phred score of 30 as the lowest cutoff for a high quality base. Subsequently artificial duplicate reads produced by sequencing artifacts 58 were removed with a k-mer based approach. For annotation unassembled reads were aligned against the non-redundant NCBI Reference Sequence (RefSeq) protein database using the program DIAMOND 59 with a minimum percentage identity cutoff of 70% for protein sequences and an e-value cutoff of 1 × 10 −10 . The top 50 hits matching the cutoff criteria for each read were retained for further analysis. Phylogenetic analysis was performed in MEGAN 5 60 using the Lowest Common Ancestor (LCA) algorithm only considering hits within the top 1% of the best bit score and a minimum bit score of 50. The LCA algorithm assigns species-specific sequences to specific taxa. Sequences that are conserved among different species (e.g. as consequence of horizontal gene transfer) will only be assigned to taxa of higher rank 60 . Nonetheless, it is very difficult to directly prove that a given (de)halogenase gene appears in a specific microbial taxon. Whenever we mention a specific species name in the results and discussion we refer to bacteria, archaea, or eukarya that contain a (de)halogenase gene closely related to the (de)halogenase gene of the respective species. Functional analysis using MEGAN 5 was based on the KEGG database and classification 61 . Each of the top 50 RefSeq hits for a read was mapped to a KEGG orthology (KO) group by identifying the best hit for a reference sequence for which a KO assignment is known. For the final assignment of a read to a KO group the KO assignment with the highest bitscore (best hit) of the assignments for the top 50 hits per read was selected. Reads related to genes of halogenating and dehalogenating enzymes were identified by analyzing reads assigned to KO groups for halogenating and dehalogenating enzymes. Since KO groups do not cover all halogenating and dehalogenating enzymes, we additionally aligned all reads with no hits to KEGG against specific databases for halogenating and dehalogenating enzymes using DIAMOND and the same cutoffs as for the RefSeq-annotation. Specific databases were created by searching the protein databases UniProt 62 and Peroxibase 63 for halogenating and dehalogenating enzymes. The specific databases include only enzymes of organisms for which halogenation or dehalogenation activity had been experimentally proven and published. "Putative enzymes" were not considered. KEGG hits and specific database hits were combined for relative abundance calculation. Abundances of functional genes were normalized to the total number of reads in the corresponding library and expressed as hits per million metagenomic reads. Abundance calculations for taxonomic groups were expressed relative to the number of all reads with a taxonomic assignment in the metagenomic library. Sequencing reads of the 12 metagenomic libraries are publically available via the MG-RAST metagenomic analysis server under project ID number 11442.

Statistical analysis.
Statistical comparison of the abundance of functional features between the soil horizons was performed using STAMP 64 applying Analysis of Variance (ANOVA) as statistical test combined with the Tukey-Kramer method as post-hoc test. If the p-value for the 95% confidence interval was below 0.05, differences were considered significant. In statistical analyses each soil horizon included the data for the forward and reverse reads of the duplicate metagenome libraries (n = 4). To visualize differences in gene abundance between the soil horizons row z-sores were calculated in R 65 . Row z-scores represent the numbers of standard deviations a value differs from the mean.

Results and Discussion
Geochemical potential for natural halogenation and dehalogenation reactions. Total organic carbon and water-leachable and therefore potentially bioavailable carbon were highest in the Of-horizon with 301 g/kg dry soil and 619.4 mg/kg dry soil, respectively (Table 1). Water-leachable AOX was highest in the Of-horizon with 0.48 mg/kg dry soil and decreased with soil depth. The performed AOX measurements only provide information on the water-leachable AOX-compounds. However, it is important to note that also the non-soluble fraction of the soil matrix contains halogenated organic compounds.
Soluble AOX gradients correlated with organic carbon and chloride gradients in the Schoenbuch soil. Especially the Of-horizon in the Schoenbuch forest was characterized by a high content of weathering plant material. Transformation of inorganic chloride during humification of plant material leads to the rapid formation of stable and less volatile aromatic organohalogen compounds 66 . Our results support previous findings in the way that the presence of both organic carbon and halide ions stimulate natural halogenation and dehalogenation reactions in soil and that elevated organic matter contents accelerates chlorination rates 67 .
Formation of volatile organohalogens in soil microcosm experiments. Besides AOX we followed the natural formation of volatile organohalogen compounds (VOX) in soil from the Schoenbuch forest. We observed the formation of chloroform (CHCl 3 ) and bromoform (CHBr 3 ) in soil microcosms after 1 h of incubation (Fig. 2). Highest VOX concentrations were observed for the Of-horizon with 2.8 ± 0.2 and 3.4 ± 0.3 μ g/kg dry soil for chloroform and bromoform, respectively.
Soils are a known natural source of chloroform 22,68 . Furthermore it was demonstrated that in presence of inorganic bromide the formation of bromoform in soil is detectable 69 . In all three soil horizons of the Schoenbuch forest bromide could not be detected by ion chromatography. However, the formation of bromoform was observed in the OF-and Ah-horizon suggesting the presence of sufficient amounts of bromide for microbial bromoform formation. Especially soils with high organic carbon content due to decaying plant material and a rich humic layer were prominent sources for chloroform 68 . This was confirmed for the Schoenbuch forest soil, from which the emissions of trihalomethanes were highest in the organic rich Of-horizon. A recent study on chloroform formation from humic substances in soils using stable isotope analysis suggested microbial formation via extracellular chloroperoxidases as potential source of VOX formation 70 . The predominance of microbial chlorination over abiotic chlorination reactions in forest soils was demonstrated by a clear temperature sensitivity of the observed chlorination reactions 29 and the significantly lower chlorination of organic matter in autoclaved and/or gamma sterilized soils 30 . Also microcosm studies on microbial dehalogenation revealed that both, anaerobic dehalogenation 32 and aerobic dehalogenation 33 of organohalogen compounds by microorganisms prevailed over abiotic reactions in the investigated soils. Both microbial halogenation and dehalogenation reactions in soils contribute to the natural halogen cycling, but so far the diversity and abundance of the involved microorganisms and enzymes have not been studied in great detail 27,28 . Since soluble and volatile organohalogen compounds were detectable Horizon Water content (%) pH TOC (g/kg) a Leachable OC (mg/kg) a Leachable AOX (mg/kg) a Leachable Cl − (mg/kg) a General information on the Schoenbuch metagenome. Metagenomic sequencing of two replicate samples per soil horizon resulted in a total of 38.8 million reads with a read length of 300 bp. After quality processing a total of 36.2 million high quality reads were used for taxonomic and functional analysis. Detailed sequencing statistics for the metagenome libraries of the duplicate soil samples are given in Table S1. Taxonomic classification was possible for 20.4% of the metagenomic reads, whereas functional annotation was possible for 8.5% of the reads. Since our study focuses on the microbial halogen cycle only sequences related to Bacteria, Archaea or Fungi were considered in our analysis. Of the reads that could be taxonomically assigned 99.5% were related to Bacteria, whereas 0.1% and 0.4% were related to Archaea and Fungi, respectively. The higher relative abundance of Bacteria over Archaea was confirmed by quantifying 16S rRNA gene copy numbers of both domains by qPCR (results shown in Table S2). 16S rRNA gene copy numbers in the three soil horizons were approximately three orders of magnitude higher for bacterial 16S rRNA genes compared to archaeal 16S rRNA genes. However, 16S rRNA gene copy numbers were not corrected for ribosomal rRNA gene operon numbers. Strong predominance of bacterial over archaeal reads in soil metagenomic libraries has also been demonstrated in a cross-metagenomic survey of 16 different soil samples 71 and metagenomics analyses of permafrost soils 72 . Bacterial reads in the Schoenbuch    (Fig. 3). Further, reads affiliated to Bacteroidetes, Actinobacteria and Verrucomicrobia constituted considerable fractions of all bacterial reads.
The dominant bacterial phyla in the Schoenbuch forest soil are typical members of soil microbial communities and represented the majority of the bacterial reads in metagenomes of e.g. desert and forest soils 71 , tallgrass prairie soils 73 and a boreal forest soil 74 . Functional metagenomic reads were mainly associated with the KEGG subsystem metabolism (43.9-45.2%) or could not directly be grouped within one of the KEGG subsystems (29.2-29.9%) ( Figure S1).

Identification of microorganisms and enzymes possibly involved in natural halogen cycling in
Schoenbuch forest soil. We screened the metagenome for microorganisms that are known to possess genes encoding for enzymes that perform halogenation or dehalogenation reactions or for which halogenation and dehalogenation reactions have been confirmed by experimental approaches (Table 2). Relative abundances of these taxa were calculated on the genus rank, since taxonomic classification at the species or strain level is not reliable for short metagenomic reads.
Bradyrhizobium and Burkholderia were the most abundant genera possessing genes for both, halogenating and dehalogenating enzymes. With the exception of eleven fungal genera and one archaeal genus all other genera belonged to the Bacteria indicating that halogen cycling might be mainly bacteria driven in the investigated forest soil. Most taxa in Table 2 are facultative aerobic microorganisms suggesting the prevalence of aerobic halogenation and dehalogenation pathways. Anaerobic bacteria known for reductive dehalogenation such as Dehalococcoides or Dehalobacter were less abundant, probably because the top 40 cm of the Schoenbuch forest soil were mainly oxic. Nonetheless, anoxic microsites in water filled micropores could sustain growth and activity of reductively dehalogenating microorganisms even in primarily oxic soil horizons.
In order to assess the genetic potential for microbial halogenation and dehalogenation reactions in the Schoenbuch forest soil, we tried to identify reads that encode for halogenating and dehalogenating enzymes. Their relative abundances in the metagenomic libraries of the duplicate samples of each soil horizon are displayed in Fig. 4. The applied metagenomic approach revealed a high genetic diversity for halogenating and dehalogenating enzymes covering a variety of different halogenation and dehalogenation mechanisms. Most retrieved halogenase genes encoded for enzymes with oxidative halogenation mechanisms. Also Vaillancourt et al. described that oxidative halogenation pathways predominate in many ecosystems 37 . Furthermore experiments on the chlorination of organic matter in forest soils suggested oxygen-dependent enzymes driving the biotic chlorination in soils 29 . For dehalogenating enzymes a variety of oxidative and reductive dehalogenation reactions are known. The majority of the dehalogenase genes we found in the Schoenbuch soil metagenome were related to hydrolytic or oxidative dehalogenases 31 . The only reductive dehalogenase genes we identified were related to a pceA gene encoding for a reductive dehalogenase that catalyses the dechlorination of perchloroethene and trichloroethene 75 . The relative abundances of the most abundant halogenase and dehalogenase genes were in the same order of magnitude as functional genes involved in microbial nitrogen cycling (nosZ, nif-genes) or housekeeping genes such as e.g. DNA or RNA polymerases. The fact that halogenase and dehalogenase genes occurred at relative abundances similar to essential soil microbial community functions emphasizes the importance of these enzymes for (de) halogenation reactions in forest soils and suggests a major role of bacteria in the cycling of halogens in soils.
The heatmap in Fig. 4 shows that the relative abundances of halogenase and dehalogenase genes in the two replicate metagenomic libraries were in the same order of magnitude and follow the same trends with soil depth. Therefore, we combined the data of both libraries for further analysis. The variance between the replicate libraries is then reflected by the given standard deviation.
We verified the relative abundances of selected halogenase, dehalogenase and reference genes involved in nitrogen cycling in the different soils horizons of the metagenome data set by qPCR. For the four selected genes  Table 2. Mean abundance of taxa known to possess enzymes for biotic halogenation or dehalogenation reactions. Relative abundances are given on the genus level. Dehalo: Genus comprises species with genetic dehalogenation potential. Halo: Genus comprises species with genetic halogenation potential. Some genera contain species that possess the genetic potential to perform both halogenation and dehalogenation reactions.
Of major interest with respect to halogenating enzymes is the proportion of genes encoding for either specific or unspecific halogenases. Genes for unspecific halogenases represented 86.7-93.5% of the total halogenase reads whereas genes for specific halogenases represented 6.5-13.3% (Fig. 5A). Unspecific halogenases increased significantly with sediment depth although differences between the Ah-and IIP-horizon were not significant. All unspecific halogenases were haloperoxidases. The higher proportion of unspecific halogenases in the deeper soil horizons might be related to their ability of reducing hydrogen peroxide to water. The rhizosphere at the The relative abundance of genes of the nitrogen cycle and of selected housekeeping genes is given as reference. Functional assignments are based on 70% amino acid sequence identity and an e-value of 1 × 10 −10 . The color code represents the row z-score, the number of standard deviations a value differs from the mean. Numeric values within the heatmap represent the relative abundance in hits per million metagenomic reads. In samples with a relative abundance of 0.0 no reads for the corresponding enzyme were found. sampling site was located at the intersection of the Ah-and IIP-horizon. Also the high abundance of nitrogen fixation genes (nif-genes) locates the rhizosphere near the IIP-horizon (Fig. 4). Haloperoxidases could be used by microorganisms as defence against oxidative stress induced by reactive oxygen species released by plants to antagonize pathogens and rhizosphere infections 76 . Specific halogenases such as flavin-dependent halogenases are involved in secondary metabolism, e.g. antibiotic synthesis 40 . 16S rRNA gene copy numbers and organic carbon content were highest in the Of-horizon suggesting that microbial competition and the necessity for production of antimicrobial agents might be high in this soil layer. This might also be a potential explanation for the high proportion of specific halogenase genes in the Of-horizon. The proportion of genes for metabolic dehalogenases significantly increased with soil depth. Metabolic dehalogenases constituted the major fraction of all dehalogenase assigned reads (59.5-71.7%), while genes encoding for cometabolic dehalogenases were less abundant (28.9-40.5%) (Fig. 5B). Many metabolic and cometabolic dehalogenases have a broad substrate specificity, e.g. the methane monooxygenases or haloalkane dehalogenases 31 .
Microorganisms using organohalogens as carbon source or electron acceptor are therefore not necessarily restricted by the availability of their primary substrate for their dehalogenating enzymes. This makes correlations between the abundance of dehalogenase genes and specific organohalogens formed during microbial halogenation reactions difficult. The relative proportion of cometabolic dehalogenases was highest in the Of-horizon. These enzymes utilize non-halogenated organic compounds as substrate and since the Of-horizon had the highest content of organic carbon, this might explain that organisms possessing monooxygenases or dioxygenases are abundant in this horizon, where they can utilize the available aromatic compounds, e.g. phenolic breakdown products of lignin degradation 77 . The ratio of halogenase to dehalogenase genes (Fig. 5C) revealed a higher abundance of dehalogenase genes in all soil horizons, whereas the ratio was closest to 1 in the Of-horizon (0.71). The observed ratio only displays the genetic potential for enzymatic halogenation or dehalogenation. Since gene expression and protein synthesis are dependent on many factors and differ strongly between different genes the relative abundance of functional genes in metagenomic datasets is no indicator of the importance of a certain function or activity in a given sample. However, since we quantified the net release of chloroform in all laboratory soil microcosms, chloroform formation must have been higher than chloroform degradation in all soil horizons.
For each halogenase and dehalogenase subgroup as classified in Fig. 4 we further investigated the distribution of the most abundant subgroup within the soil profiles. The most abundant unspecific dehalogenases were the non-heme, no metal chloroperoxidases (Fig. 6A) with 93.6-99.0 hits per million reads. No significant differences in abundance between the three soil horizons were detected. In general chloroperoxidases oxidize halides in the presence of hydrogen peroxide to the corresponding hypohalous acid, responsible for the unspecific halogenation of electron-rich organic matter 78 . The non-heme, no metal chloroperoxidases also use hydrogen peroxide for their halogenation mechanism 79 . These enzymes are also referred to as perhydrolases and show like the heme-or Ratio of halogenase and dehalogenase gene abundance in the soil horizons. A ratio of 1 represents an equal abundance and a ratio below 1 a higher abundance of dehalogenase genes. Horizons were compared using ANOVA and statistical significant differences are marked by asterisks (* p < 0.05, * * p < 0.01, * * * p < 0.001). x indicates no significant differences (p > 0.05). Only comparisons for adjacent soil horizons are shown.
vanadium containing chloroperoxidases no substrate specificity 36 . Therefore, a role in synthesis of specific compounds can be excluded. It was hypothesized that chloroperoxidases might be involved in microbial antagonism through the production of reactive chlorine species as antimicrobial agents 80 . As mentioned above the ability to reduce hydrogen peroxide to water suggests a role in oxidative stress response by microorganisms associated with the rhizosphere of plants. Many plant-associated organisms, e.g. Sinorhizobium meliloti 81 , possess non-heme, no metal chloroperoxidases 76 . Due to their reaction mechanism it is likely that chloroperoxidases are involved in the formation of chloroform in soils 51,70 . The high abundance of genes for chloroperoxidases in the Schoenbuch forest soil might be a possible explanation for the observed formation of chloroform in our microcosm experiments. Since chloroperoxidases can also use bromide 37 , these enyzmes might also play a role in the formation of bromoform.
Halogenase genes such as the flavin-dependent halogenase prnA were the most abundant specific halogenases in the dataset. The Of-horizon revealed a significantly higher abundance of prnA genes (13.6 hits per million metagenomic reads) compared to the Ah and IIP-horizon (Fig. 6B). The prnA gene encodes for a tryptophan-7 halogenase which is together with a second flavin-dependent halogenase (prnC) involved in the biosynthesis of the antifungal antibiotic pyrrolnitrin 40 . As discussed above the role of the PrnA enzyme in antibiotic synthesis might explain the higher abundance of the prnA gene in the Of-horizon. The carbon content, bacterial and archaeal cell numbers (as approximated by 16S rRNA gene copy numbers), and the abundance of fungal metagenomic reads were highest in the Of-horizon. Therefore this soil horizon might constitute a microbial habitat in which the genetic potential for the production of antimicrobials could provide a competitive advantage which might explain the increased abundance of genes involved in synthesis of an antifungal antibiotic such as pyrrolnitrin.
The abundance of haloalkane dehalogenase genes was highest in the IIP-horizon (98.0 hits per million reads) and decreased significantly in the Ah-and Of-horizon (p < 0.001) (Fig. 6C). Haloalkane dehalogenases are hydrolytic dehalogenases and have a broad substrate spectrum including various chlorinated and brominated aliphatic compounds 82 . The IIP-horizon was characterized by the lowest soluble AOX content which might be indicative of active dehalogenation mechanisms, also suggested by the high genetic potential for metabolic dehalogenation reactions in the IIP-horizon. The lower emissions of chloroform and the absence of bromoform formation in soil from the IIP-horizon could be due to an elevated activity of dehalogenating enzymes involved in the degradation of chloroform and bromoform.
Methane monooxygenase genes represented the most prominent genes among the cometabolic dehalogenases (Fig. 6D). Their abundance was highest in the IIP-horizon (15.8 hits per million reads). Comparison between directly adjacent horizons revealed no significant differences in abundance but the difference in abundance between the Oh-and IIP-horizon was significant (p < 0.01). Methanotrophs primarily use methane monooxygenases to catalyze the oxidation of methane to methanol but the enzyme can also oxidize a wide range of alkanes Figure 6. Abundance in hits per million metagenomic reads of non-heme, no metal chloroperoxidase genes (A), flavin-dependent halogenase genes (prnA) (B), haloalkane dehalogenase genes (C) and methane monooxygenase genes (D). The four enzymes are the most abundant representatives of unspecific and specific halogenases and metabolic and cometabolic dehalogenases, respectively. Horizons were compared using ANOVA and statistical significant differences are marked by asterisks (* p < 0.05, * * p < 0.01, * * * p < 0.001). x indicates no significant differences (p > 0.05). Only comparisons for adjacent soil horizons are displayed. and alkenes 83 . Methane monooxygenases are also known for cometabolic oxidation of halogenated alkenes 84 and alkanes such as chloroform 85 . The high abundance of methane monooxygenase genes in the IIP-horizon suggests the occurrence of high numbers of methanotrophic bacteria in this soil horizon. Following rainfalls parts of the IIP-horizon might quickly turn anoxic which could promote methanogenesis and the activity of methanotrophs in oxic zones since activities of methanotrophic bacteria and methanogenic archaea in soil are known to be correlated 86 .

Summary and Outlook
The microbial mechanisms driving the halogen cycle in soils are mainly unknown. Therefore knowledge on formation and degradation processes is important to evaluate the role of soils as sinks or sources of organohalogen compounds 87 . The metagenomic survey conducted in this study revealed a tremendous diversity and high abundance of genes encoding for halogenating and dehalogenating enzymes in the investigated soil. Although we did not analyze gene expression or enzyme activity we could show that the studied forest soil harbours the genetic potential for specific and unspecific halogenation as well as metabolic and cometabolic dehalogenation activities with a clear predominance of oxidative bacterial reaction pathways. However, the relative contribution of the different enzymatic groups to the overall cycling of organic and inorganic halogens in the Schoenbuch requires further study. Here we demonstrated that metagenomics allows for the identification of the diversity and relative abundance of enzymatic halogenating and dehalogenating reaction mechanisms in soils that build the basis for further investigation of microbial halogen cycling. Since halogenating and dehalogenating enzymes use different reaction mechanisms for the (de)halogenation of organic matter the contribution of individual enzymatic mechanisms to overall halogen cycling should be further elucidated by stable chlorine isotope fractionation as recently demonstrated 70,88 . The combination of omics approaches, laboratory microcosm experiments, and stable isotope analysis constitutes a powerful set of tools to further investigate the microbial contribution to natural halogenation and dehalogenation reactions in soils.