Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Geomicrobiology | Spotlight

Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs

Caitlin M. Gionfriddo, Matthew B. Stott, Jean F. Power, Jacob M. Ogorek, David P. Krabbenhoft, Ryan Wick, Kathryn Holt, Lin-Xing Chen, Brian C. Thomas, Jillian F. Banfield, John W. Moreau
Isaac Cann, Editor
Caitlin M. Gionfriddo
aSchool of Earth Sciences, The University of Melbourne, Parkville, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Caitlin M. Gionfriddo
Matthew B. Stott
bGNS Science, Wairakei Research Centre, Taupo, New Zealand
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jean F. Power
bGNS Science, Wairakei Research Centre, Taupo, New Zealand
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jacob M. Ogorek
cWisconsin Water Science Center, U.S. Geological Survey, Middleton, Wisconsin, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David P. Krabbenhoft
cWisconsin Water Science Center, U.S. Geological Survey, Middleton, Wisconsin, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ryan Wick
dDepartment of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ryan Wick
Kathryn Holt
dDepartment of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kathryn Holt
Lin-Xing Chen
eDepartment of Earth and Planetary Science, UC Berkeley, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brian C. Thomas
eDepartment of Earth and Planetary Science, UC Berkeley, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jillian F. Banfield
aSchool of Earth Sciences, The University of Melbourne, Parkville, Victoria, Australia
eDepartment of Earth and Planetary Science, UC Berkeley, Berkeley, California, USA
fDepartment of Environmental Science, Policy, and Management, UC Berkeley, Berkeley, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John W. Moreau
aSchool of Earth Sciences, The University of Melbourne, Parkville, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John W. Moreau
Isaac Cann
University of Illinois at Urbana—Champaign
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.00176-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Geothermal systems emit substantial amounts of aqueous, gaseous, and methylated mercury, but little is known about microbial influences on mercury speciation. Here, we report results from genome-resolved metagenomics and mercury speciation analysis of acidic warm springs in the Ngawha Geothermal Field (<55°C, pH <4.5), Northland Region, Aotearoa New Zealand. Our aim was to identify the microorganisms genetically equipped for mercury methylation, demethylation, or Hg(II) reduction to volatile Hg(0) in these springs. Dissolved total and methylated mercury concentrations in two adjacent springs with different mercury speciation ranked among the highest reported from natural sources (250 to 16,000 ng liter−1 and 0.5 to 13.9 ng liter−1, respectively). Total solid mercury concentrations in spring sediments ranged from 1,274 to 7,000 μg g−1. In the context of such ultrahigh mercury levels, the geothermal microbiome was unexpectedly diverse and dominated by acidophilic and mesophilic sulfur- and iron-cycling bacteria, mercury- and arsenic-resistant bacteria, and thermophilic and acidophilic archaea. By integrating microbiome structure and metagenomic potential with geochemical constraints, we constructed a conceptual model for biogeochemical mercury cycling in geothermal springs. The model includes abiotic and biotic controls on mercury speciation and illustrates how geothermal mercury cycling may couple to microbial community dynamics and sulfur and iron biogeochemistry.

IMPORTANCE Little is currently known about biogeochemical mercury cycling in geothermal systems. The manuscript presents a new conceptual model, supported by genome-resolved metagenomic analysis and detailed geochemical measurements. The model illustrates environmental factors that influence mercury cycling in acidic springs, including transitions between solid (mineral) and aqueous phases of mercury, as well as the interconnections among mercury, sulfur, and iron cycles. This work provides a framework for studying natural geothermal mercury emissions globally. Specifically, our findings have implications for mercury speciation in wastewaters from geothermal power plants and the potential environmental impacts of microbially and abiotically formed mercury species, particularly where they are mobilized in spring waters that mix with surface or groundwaters. Furthermore, in the context of thermophilic origins for microbial mercury volatilization, this report yields new insights into how such processes may have evolved alongside microbial mercury methylation/demethylation and the environmental constraints imposed by the geochemistry and mineralogy of geothermal systems.

INTRODUCTION

Geothermal springs and fumaroles emit substantial amounts of aqueous and gaseous mercury (Hg) (1). Aqueous Hg concentrations in these systems often exceed 100 ng liter−1, and total Hg levels can approach 25 μg liter−1 (2–4). Despite these ultrahigh mercury levels, few studies have examined biotic and abiotic mechanisms for Hg transformations or Hg speciation in geothermal springs. Specifically, the potential for native thermophiles to mediate mercury transformations, i.e., reduction of Hg(II) to Hg(0) or methylation/demethylation of mercury (to CH3Hg+ [MeHg] or to Hg[II], respectively) (5–7), remains poorly understood.

Hg species have a high binding affinity to thiols and lipids, causing damage to proteins, enzymes, and nucleic acids. They can inhibit microorganisms at submicromolar concentrations (8). Microorganisms living in environments with elevated Hg (>100 ng liter−1) commonly possess genes encoding Hg resistance (9, 10). The mer operon is used by many bacteria and archaea to detoxify Hg(II), by converting it to volatile Hg(0) (7, 10, 11). Interestingly, mer has a phylogenetic origin in the thermophiles (12, 13). Additionally, microbes carrying the organomercurial lyase-encoding merB gene as part of the mer operon are able to detoxify organic Hg compounds, in tandem with mercuric reductase (MerA), to produce methane (CH4) and Hg(0) (9, 14). Finally, some anaerobic bacteria are suspected to oxidatively demethylate Hg species independent of the mer pathway, producing carbon dioxide (CO2) and Hg(II) (15), although a biochemical pathway has not been identified for mer-independent demethylation.

The efficiency of microbial methylation of Hg(II) to MeHg is still largely unknown in geothermal spring ecosystems, particularly under acidic conditions (<55°C, pH <4.5) (4, 5, 16). However, the hgcAB genes required by microorganisms to methylate Hg (17) have been reported from many environments, including wetland sediments, rice paddy soils, thawing permafrost, hypersaline and hypersulfidic waters, soda lakes, and geothermal systems (16, 18–20). These environments typically host abundant sulfate- and iron-reducing bacteria (Deltaproteobacteria), as well as methanogenic and acetogenic Methanomicrobia (Euryarchaeota), Chloroflexi, and Firmicutes, all of which contain species capable of Hg methylation (18, 21–23). Metagenomic analyses have also identified hgcAB genes in Chrysiogenetes, “Candidatus Atribacteria” (candidate phylum OP9), and candidate phylum ACD79 (16). Although Hg methylation is most often associated with sulfate reduction, the existence of environmental triggers or controls on Hg methylation remains poorly understood as does the evolution and phylogenetic distribution of hgcAB genes.

By comparison, in nongeothermal environments, elevated Hg concentrations can inhibit microbial Hg methylation (24) and lead to conditions favoring MeHg demethylation (25) or methylation-demethylation cycles (26–28). However, in-depth understanding of Hg methylation in geothermal springs, similarly to acid mine drainage (AMD) (29), has not been established. Conversely, acidophiles in AMD systems have been well studied with respect to potential for metal and sulfur cycling (30), but with less focus on Hg speciation and methylation. Here, we combined metagenomic and geochemical speciation analyses to understand Hg transformations in the context of biogeochemical cycling (e.g., S and Fe) in an acidic warm spring microbiome. To understand physicochemical constraints on microbial Hg transformations, we studied Hg speciation across a gradient of environmental factors that can influence microbiome composition and activity. We compare our findings to those of previous geothermal spring studies to refine the conceptual model for geothermal Hg cycling.

RESULTS

Water and sediment chemistry.Chemical and physical data for each spring are shown in Tables 1 and 2. Chloride concentrations exceeded 400 mg liter−1 in all springs except Cub Bath (TS2) at 29.8 mg liter−1. Across the Ngawha Geothermal Field (NGF), sulfate concentrations varied from 9.3 to 1,200 mg liter−1, with the highest concentrations found in springs of pH <4 (Table 1). In Tiger Springs (TS), Cub Bath had the lowest sulfate concentration (319 mg liter−1), while those at the other TS sites were >1,000 mg liter−1. Sulfide concentrations were nearly 3.5× higher in Cub Bath (6.45 mg liter−1) than in Tiger Bath (TS1) (1.82 mg liter−1). The total iron concentration in Tiger Bath (12.7 mg liter−1) in October 2011 was twice that of Cub Bath (6.06 mg liter−1).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Anion and cation measurements in filtered water from hot springs in the Ngawha Geothermal Field

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Mercury analyses of NGF hot springs

Total Hg (HgT) and total MeHg (MeHgT) measurements of filtered waters from NGF hot springs sampled in April 2011 are shown in Table 2. In comparison to background levels for freshwater systems, HgT was relatively high at ∼250 to 16,000 ng liter−1 compared to 0.4 to 74 ng liter−1 for nongeothermal lakes and 1 to 7 ng liter−1 for rivers and streams (31). Nonthermal waters in the Ngawha region are reported to contain 300 to 500 ng liter−1 Hg (32). Previously reported values for dissolved Hg in NGF thermal waters range from 1,000 to 350,000 ng liter−1 (32). In this study, the highest levels of HgT were found in the Tiger Springs region (all >1,000 ng liter−1) and from Cinnabar Bath (GN1) in the Ginn Ngawha Spa (1,460 ng liter−1) (see Fig. S3 in the supplemental material). Methylmercury levels ranged from 0.5 to 14.0 ng liter−1; however, elevated concentrations of MeHgT did not correlate with higher HgT concentrations. Concentrations of HgT and MeHgT in filtered water samples from NGF were compared to reference values (Table 2). MeHg levels in filtered water samples from Yellowstone National Park (YNP) were substantially lower than those observed at NGF, with most YNP springs having concentrations below the detection limit (0.013 ng liter−1). Those with detectable MeHg were still relatively low at 0.026 to 0.080 ng liter−1 (5). Only three sites sampled in the NGF exhibited a significant proportion of HgT as MeHg (∼1% [vol/vol]): Cub Bath (TS2), Kotanitanga Bath (NS1B), and a drainage pool adjacent to Favourite Bath (NS3B). Methylmercury levels were nearly an order of magnitude greater in Cub Bath (TS2) than in Tiger Bath (TS1), even though the total mercury levels were roughly the same (∼1,000 ng liter−1).

Sediments collected from the Tiger Springs area in October 2011 were analyzed for HgT (Table 2). Total Hg concentrations were highest in the hole adjacent to Tiger Bath (TS3) at ∼7,000 μg g−1. Tiger Bath (TS1) sediments had approximately one-half the solid HgT of TS3, at 3,467 μg g−1, and Cub Bath sediments had approximately one-third the solid HgT of Tiger Bath (1,274 μg g−1). Gaseous Hg(0) emissions were also recorded from the baths (Table 2). The Hg(0) concentration measured above Tiger Bath (25.0 ng liter−1) was greater than that of Cub Bath (4.34 ng liter−1) but also varied across the bath (4.38 to 25.0 ng liter−1). Furthermore, these values are an order of magnitude lower than previously reported concentrations of fumarolic Hg(0) in the Tiger Springs area at 13.5 to 276 μg liter−1 and 710 μg liter−1 from Tiger Bath specifically (32).

Microbial diversity from metagenomic data sets.Phylogenetic analysis of ribosomal marker protein S3 from assembled Tiger and Cub Bath metagenomes identified members of Deltaproteobacteria, Gammaproteobacteria, Thermotogae, “Candidatus Micrarchaeota,” “Candidatus Parvarchaeota,” Thermoplasmata, and other Euryarchaeota (see Fig. S4 and S5). Furthermore, phylogenetic analysis revealed a greater breadth in the phyla of genomes resolved from Cub Bath compared to that from Tiger Bath, with Verrucomicrobia, Acidobacteria, Firmicutes, Planctomycetes, Alphaproteobacteria, and Betaproteobacteria identified in Cub Bath only (Fig. S4). Ribosomal proteins from Thermoprotei were identified only in assembled Tiger Bath metagenomic data (see Fig. S5). Genome binning resulted in acidophilic and thermophilic bacteria and archaea dominating the Tiger and Cub Bath metagenomic data sets, with the average coverage of scaffolded contigs within genome bins ranging from 7.1 to 789× (Fig. 1; Table S2). The highest coverage bins in each bath were related to Acidithiobacillus, Thermotogales, and Thermoplasma (Fig. 1).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Rank abundance by scaffold coverage of ribosomal protein S3 within binned and unbinned genomes from Tiger (NW1) and Cub (NW2) Bath metagenomes. Analyses and annotations were performed in ggKbase (https://ggkbase.berkeley.edu/). Genomic bin phylogeny used ribosomal S3 proteins from genomic bins, while scaffold phylogeny to the lowest common ancestor is given for unbinned ribosomal proteins. When multiple ribosomal protein S3s had the same taxonomic classification, the average coverage is shown, with error bars representing standard deviations. Values are ranked by NW2 coverage.

Genome bins were comprised of aerobes as well as obligate and facultative anaerobes capable of sulfur oxidation (Acidithiobacillus spp. and Thiomonas), sulfur reduction (Desulfurella acetivorans and Granulicella mallensis), iron oxidation (Ferroplasma and Acidithiobacillus ferrivorans), iron reduction (Acidobacterium capsulatum), methane oxidization (Methylacidiphilum), and acetate oxidation (Desulfurella acetivorans) (Fig. 2). To elucidate important biogeochemical links to Hg cycles mediated by these microbial phylotypes, metagenomes were searched for genes encoding Hg, sulfur, sulfate, and methane cycling (Fig. 2). We note here that genes involved in methanogenesis (notably mcrA) and methane oxidation (pmoA) were nearly absent from metagenomes (Fig. 2).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Heat map showing functional proteins from various biogeochemical pathways associated with Hg within genomic bins and unbinned (NGAWHA_1_UNK and NGAWHA_2_UNK) scaffolds from Tiger (NW1) and Cub (NW2) Bath metagenomes. Intensity of color refers to number of genes from each bin that encode the enzyme, operon, or pathway involved in mercury, sulfur, or methane cycling; values are provided in cells for reference. Genomic bins were ordered by the consensus coverage for all scaffolds within the bin (highest to lowest). Analyses and annotations were performed in ggKbase (https://ggkbase.berkeley.edu/).

Mer operon.Assembled metagenomes were screened for genes belonging to the mer operon that encode mercuric reductase (merA), organomercurial lyase (merB), a periplasmic protein (merP), and inner membrane proteins involved in Hg(II) transport (merT, merC, merE, merF, and merG), as well as one or more regulatory proteins (merR and merD) (13). At Ngawha, scaffolded mer genes were often encoded in the higher coverage genome bins of each bath, Acidithiobacillus spp., Thiomonas, and Thermotogae (see Table S5; Fig. 2). The coverage of mer scaffolds ranged from 3.9 to 836×, with scaffold lengths of 1,021 to 101,352 bp, indicating that a significant number of reads mapped to each sequence. A high fraction of reads (3.24E−04 to 1.94E−04) (see Table S1) from each metagenome were predicted to encode mercuric reductase (MerA) using the hidden Markov model (HMM). BLASTP analysis of HMM search outputs indicated that the HMM was insufficient for filtering sequences that are predicted to encode MerA paralogs dihydrolipoamide dehydrogenase and pyridine nucleotide-disulfide oxidoreductase. Therefore, the HMM likely overestimates MerA abundance encoded by the raw reads. MerA homologues identified in the assembled Tiger and Cub Bath metagenomes using ggKbase annotations were related to both Archaea (Euryarchaeota) and Bacteria (Proteobacteria and Bacteroidetes) (Fig. 3). Most archaeal MerA were related to Thermoplasma; however, four MerA homologues branched deeply from known archaeal homologues and were related (<43%, BLASTP sequence alignment) to MerA from “Candidatus Methanoperedens nitroreducens” (NW1 scaffolds 675 and 247, and NW2 scaffolds 1319 and 13696) (Fig. 3). Most MerA homologues from Cub Bath were related to Acidithiobacillus spp. (Fig. 3). MerA homologues identified in the metagenomes were from operons also encoding MerP, MerT, and MerR (Table S5). Several scaffolds related to Euryarchaeota (Thermoplasma) encoded just for MerA and MerP (NW1_scaffold_26, NW1_scaffold_675, and NW2_scaffold_876) (Table S5), consistent with Thermoplasma mer genes sequenced from acid mine drainage (AMD) (33).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Maximum likelihood tree showing MerA phylogeny of 31 sequences pulled from assembled Tiger (NW1) and Cub (NW2) Bath metagenomes using ggKbase. Genomic bins or scaffold ID (when MerA was unbinned) are given in bold font. Included in analysis are 113 MerA homologues, including 46 sequences from Yellowstone National Park (YNP) metagenomes (82). Trees were constructed using Le Gascuel amino-acid substitution model with gamma distribution in MEGA6 (77). Data were bootstrapped with 100 replications. The initial neighbor-joining tree was constructed with pairwise distances estimated using a JTT model. Positions with less than 90% site coverage (e.g., alignment gaps, missing data, or ambiguous bases) were excluded. A total of 419 positions were used in the final data set.

Two similar but distinct merB genes were identified in assembled Tiger and Cub Bath metagenomes, related to Acidithiobacillus caldus and Thioalkalivibrio spp., respectively (see Fig. S6). The merB gene from Tiger Bath (NW1_scaffold_113) was linked to a genome bin identified as Acidithiobacillus caldus, with an average scaffold coverage of 750× (Table S5). The same scaffold (NW1_scaffold_113) contained two merR genes that were convergent and divergent, respectively, to mer genes that encode MerT, MerP, and MerA. The merB gene from Cub Bath (NW2_scaffold_17399) was most closely related to merB from Thioalkavibrio spp.; however, the scaffold itself was unbinned, with an overall taxonomic identification of Acidithiobacillus ferrivorans (Table S5). The scaffold also contained merR and merA genes related to Acidithiobacillus ferrivorans (WP_035195121). The translated organomercurial lyases (MerB) from Tiger and Cub Baths aligned to conserved cysteines (34), indicative of their true functionality (see Fig. S7).

Biological Hg methylation.Both Tiger and Cub Bath metagenomic read sets were searched for sequences sharing homology to Hg methylation genes (hgcA and hgcB) (17). A small fraction of reads (2.16E−07) were identified as fragments of hgcA sequences in the Cub Bath metagenome (Table S1). The nucleotide reads were each 100 bp in length and, when translated, aligned to several regions of HgcA from known methylators, including the highly conserved G(I/V)NVWC region of the HgcA protein (17) (see Fig. S8). Several of the reads aligned to one another, and a composite amino acid sequence (59 amino acids [aa] in length) is shown in Fig. S8. Phylogenetic analyses of the translated reads revealed two distinct hgcA-like genes in the Cub Bath metagenome (Fig. 4). One of the reads closely matched (BLASTP search) to pterin-binding regions of HgcA-like proteins from marine bacteria Streptomyces sp. CNQ-509 and Nitrospina spp. (E value, 8E−05; 93% sequence coverage, 55%) and to fused HgcAB proteins primarily found in thermophilic archaeal and bacterial genomes, Thermococcus sp. strain EP1, Kosmotoga pacifica, Pyrococcus furiosus, and Methanococcoides methylutens (E values, 3E−06 to 8E−08; 100% sequence coverage, 65% to 75% sequence identification [ID]). Whether these microbes with genomes encoding an HgcAB fused protein are capable of Hg methylation is unknown (16). When tested for Hg methylation capability, both Pyrococcus furiosus and Methanococcoides methylutens were unable to produce MeHg at levels higher than controls (16, 35).

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Maximum likelihood tree showing HgcA phylogeny of reads obtained from the Cub Bath (NW2) metagenome using HMM search. Included in analysis are 183 sequences, including 56 HgcA homologues from YNP metagenomes, pulled from JGI (see Table S3 in the supplemental material). Amino acid sequences are compared to HgcA homologues from known and predicted methylators. Also included are HgcAB fused proteins from hyperthermophilic bacteria and archaea and carbon monoxide dehydrogenase/acetyl coenzyme A (acetyl-CoA) synthase subunit gamma (HgcA paralogs) from nonmethylators. Trees were inferred from the Le Gascuel amino acid substitution model with gamma distribution in MEGA6 (77). Data were bootstrapped with 100 replications. The initial neighbor-joining tree was constructed with pairwise distances estimated using a JTT model. Positions with less than 93% site coverage (e.g., alignment gaps, missing data, or ambiguous bases) were excluded. A total of 55 positions were used in the final data set. Groups (I to IX) designate distinct subtrees that contain HgcA from hot spring metagenomes. Representative reference sequences are labeled within each subtree, while phylogeny of all branches is indicated by color.

The second set of hgcA reads from Cub Bath (reads 2 to 7) (Fig. 4), including the composite sequence, were aligned using BLASTP to the pterin-binding region of HgcA proteins from known and predicted Hg methylators Desulfosporosinus youngiae, Clostridium cellobioparum, Desulfosporosinus sp. strain Tol-M, and Desulfosporosinus sp. BRH-c37 (E values, 4E−21 to 7E−22; 100% query cover; 71 to 75% ID). Phylogenetic analyses (Fig. 4) of the composite sequence revealed homology to HgcA from Nitrospirae bacterium HCH-1 (GenBank ID LNQR00000000.1) and to HgcA from hot spring metagenomes. The Cub Bath hgcA sequences could not be binned; therefore, additional taxonomic information about hgcA carrying genomes was not obtained. Taxonomic assignments of Hg methylators from HgcA phylogeny alone are not conclusive. The hgcA sequences pulled from assembled YNP and British Columbia hot spring metagenomes from JGI-IMG databases span a wide range of diverse phyla and are predicted to encode HgcAs that relate closely to those from Deltaproteobacteria, Firmicutes, Nitrospirae, Planctomycetes, and Euryarchaeota (Methanomicrobia and Thermoplasmata) (Fig. 4). BLASTP searches of translated amino acid sequences revealed common closest matches (<83% ID) to Clostridium straminisolvens, Deltaproteobacterium sp. NaphS2, Nitrospirae bacteria SG8-3 and HCH-1, Phycisphaerae bacterium SG8-4, Smithella spp., Methanocella arvoryzae, and Thermoplasmatales archaeon DG-70-1. Of these closest relatives, only Thermoplasmatales archaeon DG-70-1 is from a thermophilic phylum (Thermoplasmata); however, this strain was isolated from an anaerobic, moderately halophilic, and mesophilic aquatic environment (36) rather than a geothermal spring.

A high fraction of reads from each metagenome aligned to hgcB sequences from known and predicted methylators (2.08E−06 to 3.02E−06) (Table S1). However, BLAST analysis of the translated hgcB-like genes resulted in closest matches to ferredoxin-encoding gene sequences in both Tiger and Cub Bath. To differentiate between other ferredoxin proteins and HgcB homologues, the translated hgcB-like genes were searched for the conserved C(M/I)ECGAC site in HgcB (17). A complete hgcB gene was identified in the assembled Cub Bath metagenome (NW2_scaffold_10600_2) (see supplemental material). A translated BLASTP search of the sequence revealed the closest matches (E values, 2E−26 to 3E−29; 91% to 98% coverage; 53% to 57% identity) to HgcB were from predicted Hg methylators Dehalobacter spp., Clostridium straminisolvens, Nitrospirae bacterium HCH-1, and Syntrophobotulus glycolicus. Other features on the assembled contig that contains the hgcB gene (NW2_idba_contig_10600) include a partial hgcA gene upstream from hgcB, and downstream genes encode a DsrE/DsrF-like family protein (involved in intracellular sulfur reduction) related to Planctomycetes, a copper-binding protein related to Aquificae, a putative transcriptional regulator related to ArsR repressor, and an arsenical pump membrane protein (ArsB), both related to Firmicutes (see Fig. S9).

Biological sulfur species cycling in Tiger and Cub Baths.The availability and speciation of sulfur compounds within geothermal environments constrains the strategies employed by microorganisms to deal with Hg in these ecosystems. For example, Hg methylation is impeded in ecosystems with elevated aqueous sulfide concentrations (37). Conversely, the availability of soluble Hg to microorganisms is limited by the solubility and/or microbially mediated dissolution of cinnabar/metacinnabar (38). To elucidate these interactions, the Cub and Tiger Bath metagenomes were searched for genes related to dissimilatory sulfite/sulfate reduction using dsrAB genes that encode subunits A and B of the dissimilatory bi(sulfite) reductase enzyme (DsrAB) (39) and the sox pathway genes [soxR, soxXA, soxYZ, soxB, sox(CD)2, soxG] that encode enzymes used in hydrogen sulfide and elemental sulfur oxidation (40).

Complete or near complete dsrAB sequences were recovered from assembled Tiger and Cub Bath metagenomes using ggKbase annotations (see Table S6). The genes were often present on scaffolds containing other genes encoding proteins involved in sulfite/sulfate reduction to sulfide. Both dsrC and dsrD genes were present on scaffolds with dsrAB in NW1_Thaumarchaeota_unknown_1, NW2_Deltaproteobacteria_41_12, and NW2_Acidobacterium_capsulatum_related. Phylogenetic analyses of the translated genes indicated six distinct groups of DsrAB sequences in Tiger and Cub Baths (Fig. 5; Table S6), groups I to VI. Three distinct groups of reductive bacterial type DsrAB groups (I to III) were identified, most closely related to known thermophilic sulfur-reducing bacteria such as Desulfurella acetivorans as well as uncultured Gemmatimonas sp. strain Sg8-17. Groups IV and V contained sequences related to reductive archaeal DsrAB in genomic bins related to Thermoplasmata, Thermoplasma acidophilum, and Aciduliprofundum. Group VI DsrABs contained sequences distantly related to “Candidatus Rokubacteria” CSP1-6, Caldivirga maquilingensis, and Thermodesulfobacterium spp. (WP_051754629). Group V DsrABs were likely uncultured reductive archaeal type DsrABs related to Vulcanisaeta; unbinned sequences were from scaffolds with a majority of sequences related to Archaea. Group VI DsrABs were most closely related (>50%) to DsrAB from the sulfate-reducer “Candidatus Rokubacteria” CSP1-6 (41); the sequence was distinct from DsrA of group IV and is from genomic bin NW2_Planctomycetes_54_12 (Table S6). The DsrAB phylogenies as represented in Fig. 5 agree with the phylogenetic annotations of the scaffolds from which the genes were obtained (Table S6). Importantly, while there is evidence for bacterial and archaeal sulfate and sulfur reduction within Tiger and Cub Baths, none of the DsrAB groups detected (groups I to VI) (Fig. 5) contain known Hg methylators from the primary sulfate-reducing phyla, Deltaproteobacteria and Firmicutes (18), and there were no hgcAB genes present in sulfate-reducing genomic bins from Tiger or Cub Bath.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Phylogenetic analysis of DsrAB by maximum likelihood method. The 15 DsrAB homologues in this study were pulled from IDBA-UD assembled Tiger and Cub metagenomes using ggKbase. They were compared to 218 reference DsrAB sequences pulled from the Dome database, including 15 representative DsrAB sequences from thermophilic environments (39). There were six distinct groups of DsrAB homologues found at Ngawha. The tree was constructed using MEGA6 (77) with the Le Gascuel 2008 model with gamma distribution; pairwise distances were estimated using a JTT model. All positions with less than 95% site coverage were eliminated, including alignment gaps, missing data, and ambiguous bases. There was a total of 492 positions in the final data set, with 100 bootstrap replicates.

Genes encoding sulfur oxidation (sox) in acidophilic Betaproteobacteria and Gammaproteobacteria were detected in both Tiger and Cub Bath metagenomes and were related to Acidithiobacillus spp. and Thiomonas spp. (Fig. 2; Table S7). Furthermore, sequences from archaeal and bacterial sox pathways, including from acidophiles Acidiphilium and Acidocella, were present in the unbinned metagenomic data (Table S7).

DISCUSSION

Geothermal systems provide an environment in which relationships between the chemical and physical processes controlling Hg speciation and bioavailability and microbial Hg transformations (4, 7) remain poorly understood. Temperature and pH constitute major drivers of microbial diversity in geothermal springs, with pH contributing to a greater extent (42, 43). Indeed, previous studies show that acidic geothermal spring communities appear quite distinct from those of neutral and alkali springs, irrespective of temperature (42, 43). In our study, despite their mutual close proximity and broadly similar physicochemical properties and dissolved total Hg concentrations, the Tiger and Cub Baths of the NGF hosted very distinct microbiomes (see Fig. S10 in the supplemental material). The greater diversity in genomic bins representative of the Cub Bath microbiome than of Tiger Bath (Fig. 1) may have promoted Fe and S redox cycling to a greater extent, which in turn could facilitate the dissolution of metal sulfides such as pyrite or cinnabar. Subsequently, this dissolution could have increased the bioavailability of Hg(II) for hgcAB+ equipped microorganisms by oxidizing reduced S and increasing dissolved Hg(II) (Fig. 6).

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Conceptual model of biogeochemical cycling of mercury (Hg), sulfur (S), and iron (Fe) in Hg-enriched, sulfidic, low pH mesothermal springs. Gaseous elemental mercury [Hg(0)] [as well as Hg(II)] from deep geological sources enters the surface waters of the springs where it becomes oxidized to Hg(II) (enhanced by chloride [Cl−]) (32) and then complexes with sulfides (S2-) to produce cinnabar (red rhombohedral symbols; HgS(s)). Here, S2− is indicative of all reduced sulfide species. Round icons represent microbially mediated reactions, white are primarily aerobic-associated mechanisms, and black are primarily anaerobic. Sulfur-oxidizing bacteria (SOB) equipped with the Sox pathway along with Fe-oxidizing bacteria and archaea (FeOBA) are able to enhance dissolution of metal sulfides, such as pyrite (silver rhombohedral symbols; FeS2) and HgS(s). Sulfate-reducing bacteria (SRB) and Fe-reducing bacteria (FeRB) further mediate the redox chemistry of S, Fe, and Hg. As Hg(II) becomes bioavailable to microbes, it can be reduced to Hg(0) by microbes equipped with mercuric reductase (MerA) or methylated to MeHg by HgcAB-equipped microbes. MeHg can be demethylated by MerB-equipped microbes to Hg(II) and CH, and then reduced to Hg(0) by MerA. At the surface of the springs, photoreduction can also contribute to Hg(II) reduction to Hg(0), as well as the degradation of MeHg. Photolytic oxidation may also counter Hg(0) volatilization from surface waters, keeping Hg(II) in spring water to be transformed by microbes or partitioned to sulfide minerals. Advective mixing of spring waters ensures that Hg species travel across the redox boundaries that likely partition the Mer-equipped microbes to oxic surface waters from the HgcAB-equipped microbes that likely occupy the anaerobic sediment/water boundary.

NGF genomes featured similar metabolic capabilities to those recovered from AMD. The oxidation of reduced sulfur species by both aerobic and anaerobic chemolithotrophic microorganisms produces electrons utilized in respiration and CO2 assimilation (44), critical processes for living in AMD and/or geothermal systems (44, 45). Dominant members of the Tiger Springs’ communities, Acidithiobacillus spp. and Thiomonas spp., can utilize various metal- and sulfur-oxidizing enzymes, pathways, electron transport mechanisms, and substrates (44) to sustain activity. Acidithiobacillus ferrivorans is an obligate chemolithoautotroph and facultative anaerobe that oxidizes Fe(II); some strains are also able to utilize sulfur, thiosulfate, tetrathionate, and pyrite (46). Acidithiobacillus ferrooxidans can utilize metal sulfides to support growth (44), and Acidithiobacillus caldus is also capable of oxidizing reduced inorganic sulfur species, producing sulfate via the sox pathway (47). A number of other NGF genome bins, including several associated with Thermoplasma and Thiomonas spp., were equipped to respire using either sulfur or organic carbon (30), and Thiomonas can also oxidize arsenite [As(III)] to arsenate [As(V)] (48). The Thiomonas-like genome bin from Cub Bath showed evidence for the presence of both sox and arsenite oxidase pathways (see Fig. S2).

The microbiomes of Tiger Springs were also similar to those found in AMD with respect to stress resistance/response mechanisms for acid and heavy metals (e.g., see references 30 and 49), featuring mechanisms for pH homeostasis, for example. The highest genome coverage in each bath was associated with Acidithiobacillus caldus, a microorganism equipped with multiple heavy metal resistance pathways (ars, mer, czc, and tellurite resistance). The mer genes in Tiger Spring metagenomes were predominately found in aerobic bacterial and archaeal genome bins (Fig. 2; Table S5), particularly the mesophilic acidophiles Thermoplasma and Acidithiobacillus (Fig. 3). This finding was in sharp contrast to the diversity of merA genes in YNP metagenomic data sets, which were principally from archaeal taxa (e.g., Sulfolobales, Acidilobales, and DPANN). Based on the observations of Geesey et al. (4) that archaea dominated acid springs with high Hg content, and given that MerA homologues are often encoded in the genomes of acidophilic archaea (4, 12, 13), we expected that archaeal MerA would dominate in Tiger and Cub Baths.

While Tiger and Cub Baths are considered acidic (pH <4), they are substantially lower in temperature than springs studied in YNP and the western United States; these different observations may be due to a number of factors. First, Geesey et al. (4) noted that the number of bacterial MerA homologues detected in acid springs increases with decreasing temperature (from >73°C to <55°C). Thus, the lower temperatures at NGF may explain the presence of a larger number of bacterial MerA homologues. Second, bacterial taxa are also known to be rare in higher temperature (>65 °C) acidic (<pH 4) geothermal ecosystems (50, 51). Thus, the minimal distribution of bacteria with mer genes was most likely a function of low thermophilic acidophile diversity rather than the absence of a taxonomic capacity to transform Hg.

A notable feature of Cub Bath was the higher percentage of HgT present as MeHgT than in the Tiger Bath. We speculate that this finding reflects a greater degree of microbially mediated turnover of aqueous Hg(II) to MeHg in the former spring. However, the difference in MeHg levels between the baths could also be indicative of higher demethylation rates in Tiger Bath than in Cub Bath. While merB genes were identified in both baths, the sequencing coverage of merB scaffolds was higher in Tiger Bath (836×) than in Cub Bath (18×) (Table S5). While abiotic methylation in geothermal waters is not yet well understood (52), such processes are not likely to account for the higher MeHg levels in Cub Bath relative to those in Tiger Bath. Furthermore, NGF MeHg values were 1 to 2 orders of magnitude greater than values recorded at some YNP hot springs (5, 6). An analysis of NGF metagenomes found that hgcA and hgcB genes were only detected in Cub Bath and not Tiger Bath.

Cub Bath Hg methylation genes are most likely bacterial, belonging to a sulfate-reducing clade closely related to Firmicutes and Deltaproteobacteria, possibly from the Nitrospirae (Fig. S9). Microbial Hg methylation via active sulfate-reducing microorganisms is consistent with the greater observed amount of MeHg (∼1.3% of dissolved total Hg) alongside higher concentrations of reduced sulfur (3.5×) in Cub Bath. Of particular note was an hgcAB+ bacterium in Cub Bath metagenome (NW2_unbinned) that may represent a novel acidophilic, possibly sulfate-reducing, Hg methylator equipped with heavy metal resistance. This genome harbored hgcA and hgcB on a scaffold with a putative copper chaperone (HMA/CopZ), arsenical resistance operon repressor (ArsR), and arsenical pump membrane protein (ArsB) (Fig. S9). There are few reports of arsenate reduction by sulfate-reducing bacteria (53), although putative ars genes have been found in several species of Desulfovibrio, Desulfosporosinus, Desulfomicrobium, and Desulfotalea. Genomes of known and predicted Hg-methylating bacteria and archaea, as well as closely related non-Hg-methylating bacteria and archaea, were searched for homologous proteins to those encoded on NW2_scaffold_10600 (see Table S8). Notably, ArsR family transcriptional regulators are encoded directly upstream of the genes for HgcA in Desulfovibrio desulfuricans ND132 (DND132_1054) and Desulfomicrobium baculatum (Dbac_0377). While homologous ArsR proteins are common in genomes containing HgcA (Table S8), no genomes encode all of these proteins, although several genomes of known and predicted methylators encode copper chaperones and arsenic resistance proteins, providing a measure of confidence in genome assembly results. Notably, three homologous proteins are encoded in the genomes of Desulfosporosinus spp., often found in sulfate-rich, heavy-metal-contaminated low-pH environments (54), and Desulfosporosinus acidiphilus, with an optimum growth pH of 3.6 to 5.5, was the first acidophile observed to methylate Hg (18). Therefore, we infer that the Cub Bath bacterium is likely a similar taxon capable of both Hg methylation and As(V) reduction.

One of the strongest influences on Hg speciation and bioavailability (for methylation or volatilization) in acidic sulfidic hot springs is the formation and dissolution of cinnabar (HgS(s)), which in turn is impacted by the activity of sulfur- and iron-oxidizing microorganisms (5, 32, 37, 55, 56). A previous study (32) identified HgS(s) as the most abundant and widespread Hg-bearing mineral in the Ngawha region, an observation we confirmed by X-ray diffraction (XRD) analysis from a topsoil sample in the Tiger Springs area (Fig. S2). The observed physicochemical conditions in Tiger and Cub Baths were on the cusp of HgS(s) formation/dissolution, at pH 2 to 3 and pH >4 (57). Thus, the bioavailability of Hg(II) in the acid warm springs of the NGF reflected Hg solubility in the context of sulfur speciation and acidic conditions (37, 38). In aerobic sediments, sulfur- and iron-oxidizing bacteria and archaea may enhance the dissolution of cinnabar and increase Hg(II) bioavailability. Like sulfide, chloride can also affect Hg bioavailability to methylating microorganisms; methylation rates have been shown to correlate inversely with increasing chloride (Cl−) concentration (58, 59). Chloride concentrations in Tiger Bath (470 mg liter−1) were higher than those of typical freshwaters (47.9 mg liter−1), and in combination with low pH, could have decreased net Hg methylation rates without impacting the viability of methylating microorganisms (58). Other factors that may influence Hg bioavailability include dissolved organic material (DOM), thiol-DOM interactions (60–62), and turbidity, which likely limits photolytic Hg transformations.

Together, these findings can account for the nearly 10× amount of filtered MeHg concentrations in Cub Bath, from which hgcAB sequences were recovered, compared to that in Tiger Bath, despite their near identical HgT concentrations in filtered water samples and the nearly 3× higher solid Hg content of Tiger Bath. In contrast, no hgcA reads were detected in the Tiger Bath metagenome at the sequencing depth of our study. Thus, bioavailable Hg may be getting enzymatically reduced or recomplexed by Hg-binding ligands (e.g., Cl−, DOM, and S2−). Indeed, the low flux of Hg(0) from the surface of the bath suggests that most of the microbially reduced Hg(II) remains dissolved and potentially is continuously cycled between Hg(0) and Hg(II) redox states. Advective mixing of anoxic and oxic waters would also promote Hg cycling and exchange between aerobes and anaerobes and, likewise, mer- and hgcAB-carrying microbes, in the acid warm springs (Fig. 6).

While acidophilic microbial mats in YNP have been shown to accumulate MeHg and may actively methylate Hg (5, 6), the relative difference in pHs between the two NGF springs is minimal and therefore unlikely to account for the difference in recovery of hgcAB genes. An alternative explanation may be found in the Hg(T) concentration in the sediments of both springs. The total solid Hg concentration in Tiger Bath was nearly three times higher than that of Cub Bath (3,467 μg g−1 and 1,274 μg g−1, respectively). Similarly, a small nascent hot spring located adjacent to Tiger Bath (TS3) had extremely high concentrations of both total dissolved and solid Hg (16,700 ng liter−1 and 7,000 μg g−1, respectively). Taken together, these observations suggest that the area immediately adjacent to Tiger Bath is a highly localized Hg “hot spot,” with elevated Hg levels selecting for microbial genomes encoding Hg resistance. Hg methylation capability apparently does not confer Hg(II) resistance (63), and with the exception of two Geobacter spp., mer operon genes appear to be absent in genomes of known hgcAB+ microorganisms (14, 27). These observations are further corroborated by spike-in experiments in river water sediments that showed microbial Hg methylation was inhibited by Hg concentrations as low as 15.3 μg g−1 (24). Another plausible control on Hg methylation in the two baths may have been temperature, which was recorded as nearly 10°C different between baths during each sampling season (Table 2). Similar differences in temperature were recorded in the sediments (Tiger Bath at 51.1 to 68.6°C and Cub Bath at 43.5 to 60.3°C). Principal-component analysis of covariance between geochemical parameters measured in all springs across the NGF (Tables 1 and 2) indicates an inverse relationship between MeHg and spring temperature (see Fig. S11). Previous work shows that known hgcAB+ methylators and predicted methylators are almost exclusively mesophiles with growth optima at ∼30°C (18). Known exceptions include the thermophile and known Hg-methylator Desulfacinum hydrothermale (64), psychrotolerant and psychrophilic predicted methylators Geopsychrobacter electrodiphilus and Methanolobus psychrophilus R15, and a group of hyperthermophiles possessing a fused hgcAB-like gene of unknown methylation functionality (16).

Conceptual model of Hg speciation in warm acidic hot springs.In the NGF acid warm springs, Hg biogeochemical cycling reflects the detected or inferred microbially mediated Hg transformations as well as constraints imposed by metal sulfide solubility and vigorous microbially mediated reduced S and metal oxidation reactions. The main abiotic and biotic controls on Hg cycling in acidic mesothermal springs are depicted in the conceptual model shown in Fig. 6. The bioavailability of Hg to Mer- and HgcAB-equipped microorganisms is controlled by cinnabar precipitation and dissolution, which in turn is influenced by pH and the presence of sulfur- and iron-oxidizing and sulfate- and iron-reducing bacteria and archaea. Bioavailable Hg(II) is methylated to MeHg by microbes, putatively sulfate-reducing bacteria, equipped with HgcAB. This process is limited by demethylation of MeHg (via MerB) and reduction (via MerA) of Hg(II) to Hg(0) by facultative anaerobic and aerobic iron-cycling and sulfur-oxidizing microorganisms. The volatile Hg(0) may evolve from spring waters or be photooxidized and recycled to Hg(II). A significant sink for Hg(II) within the springs involves formation of solid HgS (metacinnabar and cinnabar) and, potentially, the adsorption of MeHg to sediments or particulate organic matter (which was not assessed here). However, under the acidic conditions, most MeHg and Hg(II) should remain dissolved and could be continuously cycled between methylating and demethylating microorganisms. Importantly, temperature (>50°C) and elevated Hg(T) concentrations will restrict microbial methylation of bioavailable Hg(II). However, as geothermal inputs mix with cooler water, the microbial Hg methylation potential increases. Therefore, surface waters and groundwaters that receive geothermal inputs, such as catchment waterways downgradient from NGF springs and discharges from hydrothermal power plants, may be important environmental pathways for MeHg mobilization and bioaccumulation.

MATERIALS AND METHODS

Site description.The Ngawha Geothermal Field (NGF) constituted the field site for investigating Hg cycling in a low pH (<4.5), elevated Hg (>100 ng liter−1), and sulfide-rich (>0.1 mg liter−1) environment (see Fig. S1 in the supplemental material). Mercury ore deposits in the NGF occur as cinnabar, metacinnabar, and native Hg(0) in association with active hot springs, fumaroles, and mud pools (65). Elemental Hg [mainly gaseous Hg(0)] travels from deep geological sources to the surface, either in hydrothermal fluids or geothermal gases, where it reacts with oxygen in the presence of chloride to form Hg(II); Hg(II) in turn reacts with dissolved sulfide (biogenic) to precipitate cinnabar (32). Roughly 33,000 kg of cinnabar ore was mined from the Tiger Springs area at Ngawha during the first half of the 20th century (65). Tiger Springs and other areas within the NGF still host an active geothermal system that releases approximately 530 kg of HgT annually, ∼44% of which is thought to be emitted to the atmosphere (32). The remaining Hg resides in the local surficial waters and sediments (32). Cinnabar precipitation was confirmed by X-ray diffraction analysis of Tiger Springs sediments and nearby topsoil (Fig. S2).

Sampling techniques.Several springs from three areas, Tiger Springs (TS), Ginn Ngawha Spa (GN), and Ngawha Springs Baths (NS), were sampled in April and October of 2011. Many of the springs were edged with boards for use as soaking baths. Water samples for Hg analyses were filtered through 0.45-μm membrane syringe filters, preserved with 1% (vol/vol) reagent-grade HCl, and stored in acid-washed high-density polyethylene (HDPE) bottles in the dark. Filtered water samples (0.45 μm) for anion analysis were stored in sterile 50-ml plastic Falcon tubes. All filtered samples were stored at 4°C. Redox potential (millivolts) and pH measurements were taken at each sampling site using an Orion Model 250A portable meter with a glass pH electrode; spring water temperature measurements were also taken at the time of sample collection. Sediment samples were collected from the floor or wall of each spring, as well as from bulk water samples, in sterile 50-ml Falcon tubes. The use of sediments collected on the bottom of the boarded spring bath and suspended in water allowed us to capture a mix of aerobic and anaerobic members of the microbial community as well as those that live across a range of temperatures, from the cooler waters (<45°C) to the hotter sediments (>55°C). Samples were stored on ice for ∼4 days until they could be transferred to laboratory storage at −80°C (i.e., transported from New Zealand to Melbourne, Australia). Sediment samples were used for HgT analyses and for whole-community DNA extractions.

Hg and MeHg analyses.Total Hg and MeHg concentrations of filtered waters and freeze-dried sediments collected in April and October 2011 were measured at the Wisconsin Water Science Center (WWSC, U.S. Geological Survey, Middleton, Wisconsin) on a Perkin-Elmer Elan 9000 quadrupole inductively coupled plasma mass spectrometer (ICP-MS) and a Brooks Rand atomic fluorescence spectrophotometer model III, respectively. Filtered water samples were analyzed within 6 months of sampling at the Wisconsin Mercury Research Lab (WMRL) of the Wisconsin Water Science Center (USGS, Middleton, Wisconsin). Filtered water samples analyzed for HgT species were treated with a BrCl solution to ensure all Hg species in the sample were oxidized to Hg(II). Prior to analysis, SnCl2 solution was added to the vials, which reduced the Hg(II) species to volatile Hg(0). The samples were then ethylated, purged with argon gas, and analyzed by gas chromatography (GC) (using Brooks Rand Autosampler and Total-Hg Purge and Trap system) in tandem with atomic fluorescence spectrometry. Methylmercury analysis of filtered waters was determined by distillation, gas chromatography separation, and speciated isotope dilution mass spectrometry using ICP-MS, according to USGS method 01–445 and WMRL standard protocols (66). Four blanks and two duplicate spikes were included in each run for quality assurance. Method detection limits for total and methylated mercury were 0.007 ng and 0.03 to 1.2 ng liter−1, respectively (depending on the dilution factor required for each sample).

Sediment samples (∼5 g [wet weight]) were freeze-dried overnight on a Heto-Drywinner vacuum system before shipment in sterile glass vials to the WMRL. For solid HgT analysis, freeze-dried samples were digested with 3:1 HCl/HNO3 overnight in a Teflon vessel. The digested sample was then oxidized with BrCl solution. Total Hg analysis was then performed by using the same procedure as for filtered water. The detection limit for the solid HgT analysis was 0.2 ng. Field blanks for Hg analyses were prepared using ultrapure reaction-grade water spiked with 1% (vol/vol) ultrapure HCl stored in each type of sampling material until analysis. The HgT and MeHgT values for field blanks were 0.32 ng liter−1 and 0.57 ng liter−1, respectively.

Hg(0) analysis.Vapor samples were collected in October 2011 at a height of 5 to 10 cm over Tiger and Cub Baths (Tiger Springs area) with SKC Anasorb sorbent tubes (model C300). These sorbent tubes are typically used to measure passive exposure to mercury in industrial settings (67, 68) as a flow rate of 2 liters min−1. The volume of air sampled was regulated using an SKC sample air pump (PCXR4). Gas samples were passed through a soda lime trap before collection on the hopcalite sorbent to trap excess condensation and neutralize acid. Once used for sampling, sorbent tubes were sealed with Teflon tape and sent to ChemCentre (WA, Australia) to be analyzed per NIOSH method 6009. The sorbent material was dissolved and oxidized in 1:1 HNO3/HCl. This solution was then diluted with distilled (DI) water, and immediately before analysis on a cold vapor atomic absorption spectrometer (CVAAS), 10% SnCl2 was added to reduce all Hg(II) to Hg(0). This Hg(0) was then purged into the CVAAS analyzer (detection limit = 0.01 μg) in an argon gas stream. Sampling of Tiger and Cub Baths was performed under similar conditions. First, a blank Anasorb tube was exposed to ambient levels of Hg (unsealed, with no air pumped through). An additional blank, that remained sealed, was also included in analyses. Blank samples were below the method detection limit (<0.01 μg). Samples taken from the same site, with various pump times (20 to 30 min) did not yield similar adsorbed Hg concentrations. A positive correlation between total Hg adsorbed versus volume of air pumped through the adsorbent trap (r2 = 0.976, n = 4), indicated that saturation of the adsorbent was not reached.

Common ion analysis.Common anion concentrations (F−, Cl−, Br−, NO3−, and SO42−) were measured in unacidified filtered water samples collected in April 2011 using a Dionex DX-120 ion chromatograph with an IonPAC As14 column (4 mm by 250 mm) in the Department of Chemistry at the University of Melbourne. Samples were stored in the dark at 4°C for 4 months between sample collection and analysis. The instrument was set to the following conditions: flow rate, 1.4 ml min−1; eluent, 4.8 mM sodium carbonate and 0.6 mM sodium bicarbonate; injection volume, 25 μl. Chromatographs were viewed on PeakNet Run System I software. The method detection limit for ion chromatography (IC) analyses was 0.0005 mg liter−1. Concentrations were determined by a five-point calibration curve generated from Dionex Combined Seven Anion Standard I. The 5× dilution standard was run in duplicates to ensure precision of the method, and both sample and method blanks were produced using Milli-Q water at the time of sampling and during analysis, respectively. Determination of sulfide in water samples was performed on filtered water samples prepared in the field for analysis in October 2011 using the methylene blue method (69) and a Hach DR 2800 spectrophotometer (method 8131). Sulfide measurements were performed in triplicates for each sample. Total Fe was measured in filtered acidified water samples collected in October 2011 using the Ferrozine method (70) with prepared reagents coordinating with those for Hach method 8147. The method detection limit for sulfide measurements was 0.02 mg liter−1 and for Fe measurements was 0.5 mg liter−1. The samples had been stored in the dark at 4°C prior to analysis. Sample blanks were prepared at time of sample collection using Milli-Q water.

DNA extraction.DNA was extracted from sediments collected from Tiger (TS1) and Cub (TS2) Baths and surrounding sites (TS3 and TS4) in October 2011 using the MO-BIO PowerMax DNA isolation kit, using an alteration to the manufacturer’s protocol. Approximately 5 g of sediment was used for extraction. Sediments were first treated with 5 ml of 10 mM Tris-Cl (solution C6), the samples were then centrifuged at 6,000 × g for 5 min, and the supernatant was decanted. Then, bead solution C1 was added, and the alternative lysis method in the protocol was followed, which replaces the 10-min bead-beating step with incubation in a water bath at 65°C for 30 min, followed by vortexing the sample for 1 min (71). This alternative lysis method was used to reduce shearing of the DNA. Duplicate DNA extractions were performed to ensure a sufficient mass of DNA was extracted for each sample. Replicate DNA extractions were concentrated onto the same spin filter and cleaned using a QIAquick PCR purification kit (Qiagen, California).

Metagenomic analysis.Approximately 113 to 216 ng of genomic DNA extracted from samples collected in October 2011 (TS1A and TS2A) was barcoded by sample and prepared for sequencing with the NexteraXT kit, and ten samples were run on a single lane of an Illumina HiSeq 2500 with 2 × 100-bp paired-end sequencing (Australian Genome Research Facility, Melbourne, AU). This produced ∼9 Gbp of sequence data for two Ngawha metagenomes, Tiger Bath (NW1) and Cub Bath (NW2), and 43 Gbp for all ten samples (Table S1). Sequences were binned by bar code, quality filtered, and trimmed to remove Illumina adapters using Trimmomatic (72). Metagenomes were examined after assembling short reads into longer contigs with IDBA-UD (73). Assembled contigs were uploaded and binned for analysis using ggKbase and have been made publicly available (http://ggkbase.berkeley.edu/).

Detection of Hg cycling genes in metagenomes.The metagenomic read sets were screened directly for sequences sharing homology with hgcA and hgcB using a hidden Markov model (HMM) method described previously (74). Metagenomic read sets were also screened directly for sequence homology with merA using an HMM search built previously (74). Assembled metagenomic sequences were searched for mer operon, hgcA, and hgcB genes using ggKbase. Genes were annotated in ggKbase by BLAST searches against the NCBI database (75). Multiple merA sequences were extracted from ggKbase, translated to amino acids using the bacterial translation table, and aligned using ClustalW (76) in MEGA6 (77) with MerA reference sequences from confirmed mer operon-possessing microorganisms (13).

Genomic binning and phylogenetic analyses.Genomic binning of metagenomic data was performed using ggKbase, based on scaffold coverage, GC content of scaffold sequences, and common taxonomy of the contigs. Bins were further refined with emergent self-organizing mapping (ESOM) (78). Scaffold coverage was calculated using Bowtie2 to map reads to the assembled sequences (79). Genome completeness was estimated from the presence of bacterial single-copy genes (51 in total) or archaeal single-copy genes (38 in total) (Table S2). Several bins, namely, NW1_Thermoplasmata_unknown_1, NW2_Desulfurella_acetivorans_33_49, and NW2_Rhodospirillales_68_7, are likely multigenome bins that could not be refined to single genomes from coverage and average GC splits. Information regarding genomic bins for each metagenome, as well as unbinned scaffolds, is provided in Table S2. Taxonomy was assigned to the metagenome-assembled genomes based on consensus classification of contigs. There are 34 ribosomal proteins considered universal among bacteria, archaea, and eukaryotes and which can be used to infer phylogeny (80). In this study, we used ribosomal protein S3, a single-copy gene present in every genomic bin from this study, to compare phylogenies and coverage across each metagenome-assembled genome. Ribosomal protein sequences were pulled from ggKbase and then aligned to translated ribosomal protein S3 sequences from the NCBI nonredundant protein database.

Hg cycling genes from publicly available hot spring metagenomes.To compare Hg cycling genes from NGF to those of other geothermal settings, assembled metagenomes from Yellowstone National Park (YNP) were searched for hgcA (Table S3) and merA genes (Table S4). Included in the phylogenetic analysis were HgcA protein sequences extracted from YNP metagenomic data sets (https://img.jgi.doe.gov) (Table S3). Sequences were obtained by querying all assembled hot spring metagenomic data sets for sequences sharing homology to carbon monoxide dehydrogenase (pfam03599). HgcA proteins share homology with carbon monoxide dehydrogenases and are often misannotated as such in microbial genomes (17). HgcA sequences were differentiated from carbon monoxide dehydrogenases using the HMM described above, with an inclusion value cutoff of 1E−7, and only sequences that included the conserved cap-helix region of HgcA were kept for analyses. Of a total of 234 hot spring metagenomes (of which 216 were from YNP springs) searched (as of July 2016), there were 4,520 matches to sequences annotated as carbon monoxide dehydrogenases (pfam03599) (2,148 in YNP metagenomes). Of these sequences, 102 were identified as HgcA sequences from 18 different metagenomes (Table S3). Included in the phylogenetic analyses were 65 HgcA sequences found in 15 YNP metagenomes. The hgcA genes were found in YNP metagenomes sampled from four sites: Mushroom Springs (Gp0111644–45 and Gp0057794), Octopus Springs (Gp0057360, Gp0057796, Gp00111632, Gp0111634, Gp0111638, Gp0111642, and Gp0111646), Obsidian Pool (Gp0056876), and Fairy Spring (Gp0051404). Environmental metadata were only provided for a small number of these metagenomes. Sequences from Fairy Geyser Spring (Gp0051404) are from an alkaline (pH 9 to 9.2) and mesophilic (33.3 to 36°C) phototroph-dominated mat, while sequences from Mushroom Spring (Gp0111644–45 and Gp0057794) were from the undermat layer (∼3 to 5 mm) of an anoxygenic and phototrophic microbial mat in an effluent channel of the alkaline bath. The water above the mat had a recorded temperature of 60°C (81). A separate metagenomic study has reported temperature and pH measurements for Mushroom Spring (60°C, pH 8.2), Obsidian Pool (56°C, pH 5.7), and Octopus Spring (80 to 82°C, pH 7.9) (82). The additional 37 hgcA sequences were from assembled metagenomes sequenced from Dewar Creek Spring (77°C, pH 8.0) and Larsen North Spring, in British Columbia, Canada (83).

Data availability.The unassembled raw metagenomic sequencing data from this study are publicly available from the NCBI Sequence Read Archive (SRA). The assembled (IDBA-UD) metagenomic contigs from this study are publicly available from the NCBI Whole Genome Shotgun (WGS) project. The sequencing data can be found under the BioProject accession number PRJNA622280, with the BioSample accession numbers SAMN14506364 for the Tiger Bath and SAMN14506365 for Cub Bath metagenomes. This WGS project has been deposited at DDBJ/ENA/GenBank under the accession numbers JABEBW000000000 and JABEBX000000000.

ACKNOWLEDGMENTS

Samples were collected with permission from landowners, mana whenua, and kaitiaki for a range geothermal features in the Ngawha Geothermal Field. We thank these groups for their support. We also thank Karen Houghton at GNS and Mike Tate, John DeWild, and Charlie Thompson at the USGS Mercury Lab for their analytical and logistical support.

J.W.M. and C.M.G. were funded in part by an Early Career Researcher Award to J.W.M. from The University of Melbourne. J.F.P. and M.B.S. were funded by the geothermal science funding from GNS Science.

FOOTNOTES

    • Received 3 February 2020.
    • Accepted 7 May 2020.
    • Accepted manuscript posted online 15 May 2020.
  • Supplemental material is available online only.

  • Copyright © 2020 American Society for Microbiology.

All Rights Reserved.

REFERENCES

  1. 1.↵
    1. Gustin MS,
    2. Lindberg SE,
    3. Weisberg PJ
    . 2008. An update on the natural sources and sinks of atmospheric mercury. Appl Geochem 23:482–493. doi:10.1016/j.apgeochem.2007.12.010.
    OpenUrlCrossRefWeb of Science
  2. 2.↵
    1. Ball JW,
    2. McCleskey RB,
    3. Nordstrom DK,
    4. Holloway JM
    . 2008. Water-chemistry data for selected springs, geysers, and streams in Yellowstone National Park, Wyoming, 2003-2005. US Geological Survey, Reston, VA.
  3. 3.↵
    1. Mariner R,
    2. Presser T,
    3. Evans WC
    . 1977. Hot springs of the central Sierra Nevada, California. US Geological Survey, Reston, VA.
  4. 4.↵
    1. Geesey GG,
    2. Barkay T,
    3. King S
    . 2016. Microbes in mercury-enriched geothermal springs in western North America. Sci Total Environ 569-570:321–331. doi:10.1016/j.scitotenv.2016.06.080.
    OpenUrlCrossRef
  5. 5.↵
    1. King SA,
    2. Behnke S,
    3. Slack K,
    4. Krabbenhoft DP,
    5. Nordstrom DK,
    6. Burr MD,
    7. Striegl RG
    . 2006. Mercury in water and biomass of microbial communities in hot springs of Yellowstone National Park, USA. Appl Geochem 21:1868–1879. doi:10.1016/j.apgeochem.2006.08.004.
    OpenUrlCrossRef
  6. 6.↵
    1. Boyd ES,
    2. King S,
    3. Tomberlin JK,
    4. Nordstrom DK,
    5. Krabbenhoft DP,
    6. Barkay T,
    7. Geesey GG
    . 2009. Methylmercury enters an aquatic food web through acidophilic microbial mats in Yellowstone National Park, Wyoming. Environ Microbiol 11:950–959. doi:10.1111/j.1462-2920.2008.01820.x.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Wang Y,
    2. Boyd E,
    3. Crane S,
    4. Lu-Irving P,
    5. Krabbenhoft D,
    6. King S,
    7. Dighton J,
    8. Geesey G,
    9. Barkay T
    . 2011. Environmental conditions constrain the distribution and diversity of archaeal merA in Yellowstone National Park, Wyoming, U.S.A. Microb Ecol 62:739–752. doi:10.1007/s00248-011-9890-z.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Robinson JB,
    2. Tuovinen OH
    . 1984. Mechanisms of microbial resistance and detoxification of mercury and organomercury compounds: physiological, biochemical, and genetic analyses. Microbiol Rev 48:95–124. doi:10.1128/MMBR.48.2.95-124.1984.
    OpenUrlFREE Full Text
  9. 9.↵
    1. Schaefer JK,
    2. Yagi J,
    3. Reinfelder JR,
    4. Cardona T,
    5. Ellickson KM,
    6. Tel-Or S,
    7. Barkay T
    . 2004. Role of the bacterial organomercury lyase (MerB) in controlling methylmercury accumulation in mercury-contaminated natural waters. Environ Sci Technol 38:4304–4311. doi:10.1021/es049895w.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Simbahan J,
    2. Kurth E,
    3. Schelert J,
    4. Dillman A,
    5. Moriyama E,
    6. Jovanovich S,
    7. Blum P
    . 2005. Community analysis of a mercury hot spring supports occurrence of domain-specific forms of mercuric reductase. Appl Environ Microbiol 71:8836–8845. doi:10.1128/AEM.71.12.8836-8845.2005.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Chatziefthimiou AD,
    2. Crespo-Medina M,
    3. Wang Y,
    4. Vetriani C,
    5. Barkay T
    . 2007. The isolation and initial characterization of mercury resistant chemolithotrophic thermophilic bacteria from mercury rich geothermal springs. Extremophiles 11:469–479. doi:10.1007/s00792-007-0065-2.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Barkay T,
    2. Kritee K,
    3. Boyd E,
    4. Geesey G
    . 2010. A thermophilic bacterial origin and subsequent constraints by redox, light and salinity on the evolution of the microbial mercuric reductase. Environ Microbiol 12:2904–2917. doi:10.1111/j.1462-2920.2010.02260.x.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    1. Boyd E,
    2. Barkay T
    . 2012. The mercury resistance operon: from an origin in a geothermal environment to an efficient detoxification machine. Front Microbiol 3:349. doi:10.3389/fmicb.2012.00349.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Osborn AM,
    2. Bruce KD,
    3. Strike P,
    4. Ritchie DA
    . 1997. Distribution, diversity and evolution of the bacterial mercury resistance (mer) operon. FEMS Microbiol Rev 19:239–262. doi:10.1111/j.1574-6976.1997.tb00300.x.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    1. Oremland RS,
    2. Culbertson CW,
    3. Winfrey MR
    . 1991. Methyl mercury decomposition in sediments and bacterial cultures: involvement of methanogens and sulfate reducers in oxidative demethylation. Appl Environ Microbiol 57:130–137. doi:10.1128/AEM.57.1.130-137.1991.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Podar M,
    2. Gilmour CC,
    3. Brandt CC,
    4. Soren A,
    5. Brown SD,
    6. Crable BR,
    7. Palumbo AV,
    8. Somenahally AC,
    9. Elias DA
    . 2015. Global prevalence and distribution of genes and microorganisms involved in mercury methylation. Sci Adv 1:e1500675. doi:10.1126/sciadv.1500675.
    OpenUrlFREE Full Text
  17. 17.↵
    1. Parks JM,
    2. Johs A,
    3. Podar M,
    4. Bridou R,
    5. Hurt RA,
    6. Smith SD,
    7. Tomanicek SJ,
    8. Qian Y,
    9. Brown SD,
    10. Brandt CC,
    11. Palumbo AV,
    12. Smith JC,
    13. Wall JD,
    14. Elias DA,
    15. Liang L
    . 2013. The genetic basis for bacterial mercury methylation. Science 339:1332–1335. doi:10.1126/science.1230667.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Gilmour CC,
    2. Podar M,
    3. Bullock AL,
    4. Graham AM,
    5. Brown SD,
    6. Somenahally AC,
    7. Johs A,
    8. Hurt RA, Jr,
    9. Bailey KL,
    10. Elias DA
    . 2013. Mercury methylation by novel microorganisms from new environments. Environ Sci Technol 47:11810–11820. doi:10.1021/es403075t.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Liu Y-R,
    2. Yu R-Q,
    3. Zheng Y-M,
    4. He J-Z
    . 2014. Analysis of the microbial community structure by monitoring an Hg methylation gene (hgcA) in paddy soils along an Hg gradient. Appl Environ Microbiol 80:2874–2879. doi:10.1128/AEM.04225-13.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Schaefer J,
    2. Kronberg RM,
    3. Morel F,
    4. Skyllberg U
    . 2014. Detection of a key Hg methylation gene, hgcA, in wetland soils. Environ Microbiol Rep 6:441–447. doi:10.1111/1758-2229.12136.
    OpenUrlCrossRef
  21. 21.↵
    1. Compeau G,
    2. Bartha R
    . 1985. Sulfate-reducing bacteria: principal methylators of mercury in anoxic estuarine sediment. Appl Environ Microbiol 50:498–502. doi:10.1128/AEM.50.2.498-502.1985.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Fleming EJ,
    2. Mack EE,
    3. Green PG,
    4. Nelson DC
    . 2006. Mercury methylation from unexpected sources: molybdate-inhibited freshwater sediments and an iron-reducing bacterium. Appl Environ Microbiol 72:457–464. doi:10.1128/AEM.72.1.457-464.2006.
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    1. Kerin EJ,
    2. Gilmour CC,
    3. Roden E,
    4. Suzuki M,
    5. Coates J,
    6. Mason R
    . 2006. Mercury methylation by dissimilatory iron-reducing bacteria. Appl Environ Microbiol 72:7919–7921. doi:10.1128/AEM.01602-06.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Chen Y,
    2. Bonzongo J,
    3. Miller G
    . 1996. Levels of methylmercury and controlling factors in surface sediments of the Carson River system, Nevada. Environ Pollut 92:281–287. doi:10.1016/0269-7491(95)00112-3.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Oremland RS,
    2. Miller LG,
    3. Dowdle P,
    4. Connell T,
    5. Barkay T
    . 1995. Methylmercury oxidative degradation potentials in contaminated and pristine sediments of the Carson River, Nevada. Appl Environ Microbiol 61:2745–2753. doi:10.1128/AEM.61.7.2745-2753.1995.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Pak KR,
    2. Bartha R
    . 1998. Mercury methylation and demethylation in anoxic lake sediments and by strictly anaerobic bacteria. Appl Environ Microbiol 64:1013–1017. doi:10.1128/AEM.64.3.1013-1017.1998.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Lu X,
    2. Liu Y,
    3. Johs A,
    4. Zhao L,
    5. Wang T,
    6. Yang Z,
    7. Lin H,
    8. Elias DA,
    9. Pierce EM,
    10. Liang L,
    11. Barkay T,
    12. Gu B
    . 2016. Anaerobic mercury methylation and demethylation by Geobacter bemidjiensis Bem. Environ Sci Technol 50:4366–4373. doi:10.1021/acs.est.6b00401.
    OpenUrlCrossRef
  28. 28.↵
    1. Avramescu M-L,
    2. Yumvihoze E,
    3. Hintelmann H,
    4. Ridal J,
    5. Fortin D,
    6. Lean DR
    . 2011. Biogeochemical factors influencing net mercury methylation in contaminated freshwater sediments from the St. Lawrence River in Cornwall, Ontario, Canada. Sci Total Environ 409:968–978. doi:10.1016/j.scitotenv.2010.11.016.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Rytuba JJ
    . 2003. Mercury from mineral deposits and potential environmental impact. Environ Geol 43:326–338. doi:10.1007/s00254-002-0629-5.
    OpenUrlCrossRef
  30. 30.↵
    1. Baker BJ,
    2. Banfield JF
    . 2003. Microbial communities in acid mine drainage. FEMS Microbiol Ecol 44:139–152. doi:10.1016/S0168-6496(03)00028-X.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    EPA. 2001. Water quality criterion for the protection of human health: methylmercury: final. Office of Science and Technology, Office of Water, U.S. Environmental Protection Agency, Washington, DC.
  32. 32.↵
    1. Davey HA,
    2. van Moort JC
    . 1986. Current mercury deposition at Ngawha Springs, New Zealand. Appl Geochem 1:75–93. doi:10.1016/0883-2927(86)90039-9.
    OpenUrlCrossRef
  33. 33.↵
    1. Yelton AP,
    2. Comolli LR,
    3. Justice NB,
    4. Castelle C,
    5. Denef VJ,
    6. Thomas BC,
    7. Banfield JF
    . 2013. Comparative genomics in acid mine drainage biofilm communities reveals metabolic and structural differentiation of co-occurring archaea. BMC Genomics 14:485. doi:10.1186/1471-2164-14-485.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Pitts KE,
    2. Summers AO
    . 2002. The roles of thiols in the bacterial organomercurial lyase (MerB). Biochemistry 41:10287–10296. doi:10.1021/bi0259148.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Gilmour CC,
    2. Bullock AL,
    3. McBurney A,
    4. Podar M,
    5. Elias DA
    . 2018. Robust mercury methylation across diverse methanogenic Archaea. mBio 9:e02403-17. doi:10.1128/mBio.02403-17.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Baker BJ,
    2. Lazar CS,
    3. Teske AP,
    4. Dick GJ
    . 2015. Genomic resolution of linkages in carbon, nitrogen, and sulfur cycling among widespread estuary sediment bacteria. Microbiome 3:14. doi:10.1186/s40168-015-0077-6.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Gilmour CC,
    2. Henry EA
    . 1991. Mercury methylation in aquatic systems affected by acid deposition. Environ Pollut 71:131–169. doi:10.1016/0269-7491(91)90031-Q.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Jew A,
    2. Behrens S,
    3. Rytuba J,
    4. Kappler A,
    5. Spormann A,
    6. Brown G, Jr
    . 2014. Microbially enhanced dissolution of HgS in an acid mine drainage system in the California Coast Range. Geobiology 12:20–33. doi:10.1111/gbi.12066.
    OpenUrlCrossRef
  39. 39.↵
    1. Müller AL,
    2. Kjeldsen KU,
    3. Rattei T,
    4. Pester M,
    5. Loy A
    . 2015. Phylogenetic and environmental diversity of DsrAB-type dissimilatory (bi)sulfite reductases. ISME J 9:1152–1165. doi:10.1038/ismej.2014.208.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Friedrich CG,
    2. Bardischewsky F,
    3. Rother D,
    4. Quentmeier A,
    5. Fischer J
    . 2005. Prokaryotic sulfur oxidation. Curr Opin Microbiol 8:253–259. doi:10.1016/j.mib.2005.04.005.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Hug LA,
    2. Baker BJ,
    3. Anantharaman K,
    4. Brown CT,
    5. Probst AJ,
    6. Castelle CJ,
    7. Butterfield CN,
    8. Hernsdorf AW,
    9. Amano Y,
    10. Ise K,
    11. Suzuki Y,
    12. Dudek N,
    13. Relman DA,
    14. Finstad KM,
    15. Amundson R,
    16. Thomas BC,
    17. Banfield JF
    . 2016. A new view of the tree of life. Nat Microbiol 1:16048. doi:10.1038/nmicrobiol.2016.48.
    OpenUrlCrossRef
  42. 42.↵
    1. Power JF,
    2. Carere CR,
    3. Lee CK,
    4. Wakerley GLJ,
    5. Evans DW,
    6. Button M,
    7. White D,
    8. Climo MD,
    9. Hinze AM,
    10. Morgan XC,
    11. McDonald IR,
    12. Cary SC,
    13. Stott MB
    . 2018. Microbial biogeography of 925 geothermal springs in New Zealand. Nat Commun 9:2876. doi:10.1038/s41467-018-05020-y.
    OpenUrlCrossRef
  43. 43.↵
    1. Sharp CE,
    2. Brady AL,
    3. Sharp GH,
    4. Grasby SE,
    5. Stott MB,
    6. Dunfield PF
    . 2014. Humboldt’s spa: microbial diversity is controlled by temperature in geothermal environments. ISME J 8:1166–1174. doi:10.1038/ismej.2013.237.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Ghosh W,
    2. Dam B
    . 2009. Biochemistry and molecular biology of lithotrophic sulfur oxidation by taxonomically and ecologically diverse bacteria and archaea. FEMS Microbiol Rev 33:999–1043. doi:10.1111/j.1574-6976.2009.00187.x.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Huber JA,
    2. Butterfield DA,
    3. Baross JA
    . 2003. Bacterial diversity in a subseafloor habitat following a deep-sea volcanic eruption. FEMS Microbiol Ecol 43:393–409. doi:10.1111/j.1574-6941.2003.tb01080.x.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    1. Barahona S,
    2. Dorador C,
    3. Zhang R,
    4. Aguilar P,
    5. Sand W,
    6. Vera M,
    7. Remonsellez F
    . 2014. Isolation and characterization of a novel Acidithiobacillus ferrivorans strain from the Chilean Altiplano: attachment and biofilm formation on pyrite at low temperature. Res Microbiol 165:782–793. doi:10.1016/j.resmic.2014.07.015.
    OpenUrlCrossRef
  47. 47.↵
    1. Chen L,
    2. Ren Y,
    3. Lin J,
    4. Liu X,
    5. Pang X,
    6. Lin J
    . 2012. Acidithiobacillus caldus sulfur oxidation model based on transcriptome analysis between the wild type and sulfur oxygenase reductase defective mutant. PLoS One 7:e39470. doi:10.1371/journal.pone.0039470.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Duquesne K,
    2. Lieutaud A,
    3. Ratouchniak J,
    4. Muller D,
    5. Lett MC,
    6. Bonnefoy V
    . 2008. Arsenite oxidation by a chemoautotrophic moderately acidophilic Thiomonas sp.: from the strain isolation to the gene study. Environ Microbiol 10:228–237. doi:10.1111/j.1462-2920.2007.01447.x.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Chen L-X,
    2. Méndez-García C,
    3. Dombrowski N,
    4. Servín-Garcidueñas LE,
    5. Eloe-Fadrosh EA,
    6. Fang B-Z,
    7. Luo Z-H,
    8. Tan S,
    9. Zhi X-Y,
    10. Hua Z-S,
    11. Martinez-Romero E,
    12. Woyke T,
    13. Huang L-N,
    14. Sánchez J,
    15. Peláez AI,
    16. Ferrer M,
    17. Baker BJ,
    18. Shu W-S
    . 2018. Metabolic versatility of small archaea Micrarchaeota and Parvarchaeota. ISME J 12:756–775. doi:10.1038/s41396-017-0002-z.
    OpenUrlCrossRef
  50. 50.↵
    1. Colman DR,
    2. Poudel S,
    3. Hamilton TL,
    4. Havig JR,
    5. Selensky MJ,
    6. Shock EL,
    7. Boyd ES
    . 2018. Geobiological feedbacks and the evolution of thermoacidophiles. ISME J 12:225–236. doi:10.1038/ismej.2017.162.
    OpenUrlCrossRef
  51. 51.↵
    1. Ward L,
    2. Taylor MW,
    3. Power JF,
    4. Scott BJ,
    5. McDonald IR,
    6. Stott MB
    . 2017. Microbial community dynamics in Inferno Crater Lake, a thermally fluctuating geothermal spring. ISME J 11:1158–1167. doi:10.1038/ismej.2016.193.
    OpenUrlCrossRef
  52. 52.↵
    1. Nagase H,
    2. Ose Y,
    3. Sato T,
    4. Ishikawa T
    . 1984. Mercury methylation by compounds in humic material. Sci Total Environ 32:147–156. doi:10.1016/0048-9697(84)90127-X.
    OpenUrlCrossRef
  53. 53.↵
    1. Barton LL,
    2. Tomei-Torres FA,
    3. Xu H,
    4. Zocco T
    . 2015. Metabolism of metals and metalloids by the sulfate-reducing bacteria, p 57–83. In Saffarini D (ed), Bacteria-metal interactions. Springer, Cham, Switzerland.
  54. 54.↵
    1. Pester M,
    2. Knorr K-H,
    3. Friedrich MW,
    4. Wagner M,
    5. Loy A
    . 2012. Sulfate-reducing microorganisms in wetlands–fameless actors in carbon cycling and climate change. Front Microbiol 3:72. doi:10.3389/fmicb.2012.00072.
    OpenUrlCrossRefPubMed
  55. 55.↵
    1. Baldi F,
    2. Olson GJ
    . 1987. Effects of cinnabar on pyrite oxidation by Thiobacillus ferrooxidans and cinnabar mobilization by a mercury-resistant strain. Appl Environ Microbiol 53:772–776. doi:10.1128/AEM.53.4.772-776.1987.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Benoit JM,
    2. Gilmour CC,
    3. Mason RP,
    4. Heyes A
    . 1999. Sulfide controls on mercury speciation and bioavailability to methylating bacteria in sediment. Environ Sci Technol 33:951–957. doi:10.1021/es9808200.
    OpenUrlCrossRefWeb of Science
  57. 57.↵
    1. Barnett MO,
    2. Harris LA,
    3. Turner RR,
    4. Stevenson RJ,
    5. Henson TJ,
    6. Melton RC,
    7. Hoffman DP
    . 1997. Formation of mercuric sulfide in soil. Environ Sci Technol 31:3037–3043. doi:10.1021/es960389j.
    OpenUrlCrossRefWeb of Science
  58. 58.↵
    1. Barkay T,
    2. Gillman M,
    3. Turner RR
    . 1997. Effects of dissolved organic carbon and salinity on bioavailability of mercury. Appl Environ Microbiol 63:4267–4271. doi:10.1128/AEM.63.11.4267-4271.1997.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Boyd ES,
    2. Yu R-Q,
    3. Barkay T,
    4. Hamilton TL,
    5. Baxter BK,
    6. Naftz DL,
    7. Marvin-DiPasquale M
    . 2017. Effect of salinity on mercury methylating benthic microbes and their activities in Great Salt Lake, Utah. Sci Total Environ 581–582:495–506. doi:10.1016/j.scitotenv.2016.12.157.
    OpenUrlCrossRef
  60. 60.↵
    1. Moreau JW,
    2. Gionfriddo CM,
    3. Krabbenhoft DP,
    4. Ogorek JM,
    5. DeWild JF,
    6. Aiken GR,
    7. Roden EE
    . 2015. The effect of natural organic matter on mercury methylation by Desulfobulbus propionicus 1pr3. Front Microbiol 6:1389. doi:10.3389/fmicb.2015.01389.
    OpenUrlCrossRef
  61. 61.↵
    1. Graham AM,
    2. Aiken GR,
    3. Gilmour CC
    . 2012. Dissolved organic matter enhances microbial mercury methylation under sulfidic conditions. Environ Sci Technol 46:2715–2723. doi:10.1021/es203658f.
    OpenUrlCrossRefPubMed
  62. 62.↵
    1. Ndu U,
    2. Christensen GA,
    3. Rivera NA,
    4. Gionfriddo CM,
    5. Deshusses MA,
    6. Elias DA,
    7. Hsu-Kim H
    . 2018. Quantification of mercury bioavailability for methylation using diffusive gradient in thin-film samplers. Environ Sci Technol 52:8521–8529. doi:10.1021/acs.est.8b00647.
    OpenUrlCrossRef
  63. 63.↵
    1. Gilmour CC,
    2. Elias DA,
    3. Kucken AM,
    4. Brown SD,
    5. Palumbo AV,
    6. Schadt CW,
    7. Wall JD
    . 2011. Sulfate-reducing bacterium Desulfovibrio desulfuricans ND132 as a model for understanding bacterial mercury methylation. Appl Environ Microbiol 77:3938–3951. doi:10.1128/AEM.02993-10.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Malcolm EG,
    2. Schaefer JK,
    3. Ekstrom EB,
    4. Tuit CB,
    5. Jayakumar A,
    6. Park H,
    7. Ward BB,
    8. Morel FM
    . 2010. Mercury methylation in oxygen deficient zones of the oceans: no evidence for the predominance of anaerobes. Mar Chem122:11–19. doi:10.1016/j.marchem.2010.08.004.
    OpenUrlCrossRef
  65. 65.↵
    1. Christie T,
    2. Brathwaite B
    . 2002. Mineral commodity report 8: mercury, p 32. Institute of Geological and Nuclear Sciences Ltd., Lower Hutt, New Zealand.
  66. 66.↵
    1. Olson ML,
    2. DeWild JF
    . 1999. Techniques for the collection and species-specific analysis of low levels of mercury in water, sediment, and biota, p 191–200. In Morganwalp DW, Buxton HT (ed), U.S. Geological Survey Toxic Substances Hydrology Program—Proceedings of the Technical Meeting, Charleston, South Carolina, March 8–12, 1999—Volume 2 of 3—Contamination of Hydrologic Systems and Related Ecosystems: U.S. Geological Survey Water-Resources Investigations Report 99–4018. US Geological Survey, Reston, VA.
  67. 67.↵
    1. Takaya M,
    2. Joeng JY,
    3. Ishihara N,
    4. Serita F,
    5. Kohyama N
    . 2006. Field evaluation of mercury vapor analytical methods: comparison of the “double amalgam method” and ISO 17733. Ind Health 44:287–290. doi:10.2486/indhealth.44.287.
    OpenUrlCrossRefPubMed
  68. 68.↵
    1. Gionfriddo C,
    2. Ogorek J,
    3. Butcher M,
    4. Krabbenhoft D,
    5. Moreau J
    . 2015. Mercury distribution and mobility at the abandoned Puhipuhi mercury mine, Northland, New Zealand. N Z J Geol Geophys 58:78–87. doi:10.1080/00288306.2014.979840.
    OpenUrlCrossRef
  69. 69.↵
    1. Cline JD
    . 1969. Spectrophotometric determination of hydrogen sulfide in natural waters 1. Limnol Oceanogr 14:454–458. doi:10.4319/lo.1969.14.3.0454.
    OpenUrlCrossRef
  70. 70.↵
    1. Stookey LL
    . 1970. Ferrozine—a new spectrophotometric reagent for iron. Anal Chem 42:779–781. doi:10.1021/ac60289a016.
    OpenUrlCrossRef
  71. 71.↵
    1. Zhou J,
    2. Bruns MA,
    3. Tiedje JM
    . 1996. DNA recovery from soils of diverse composition. Appl Environ Microbiol 62:316–322. doi:10.1128/AEM.62.2.316-322.1996.
    OpenUrlAbstract/FREE Full Text
  72. 72.↵
    1. Bolger AM,
    2. Lohse M,
    3. Usadel B
    . 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. doi:10.1093/bioinformatics/btu170.
    OpenUrlCrossRefPubMedWeb of Science
  73. 73.↵
    1. Peng Y,
    2. Leung HC,
    3. Yiu S-M,
    4. Chin FY
    . 2012. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 28:1420–1428. doi:10.1093/bioinformatics/bts174.
    OpenUrlCrossRefPubMedWeb of Science
  74. 74.↵
    1. Gionfriddo CM,
    2. Tate MT,
    3. Wick RR,
    4. Schultz MB,
    5. Zemla A,
    6. Thelen MP,
    7. Schofield R,
    8. Krabbenhoft DP,
    9. Holt KE,
    10. Moreau JW
    . 2016. Microbial mercury methylation in Antarctic sea ice. Nat Microbiol 1:16127. doi:10.1038/nmicrobiol.2016.127.
    OpenUrlCrossRef
  75. 75.↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    . 1990. Basic local alignment search tool. J Mol Biol 215:403–410. doi:10.1016/S0022-2836(05)80360-2.
    OpenUrlCrossRefPubMedWeb of Science
  76. 76.↵
    1. Larkin MA,
    2. Blackshields G,
    3. Brown NP,
    4. Chenna R,
    5. McGettigan PA,
    6. McWilliam H,
    7. Valentin F,
    8. Wallace IM,
    9. Wilm A,
    10. Lopez R,
    11. Thompson JD,
    12. Gibson TJ,
    13. Higgins DG
    . 2007. ClustalW and ClustalX version 2. Bioinformatics 23:2947–2948. doi:10.1093/bioinformatics/btm404.
    OpenUrlCrossRefPubMedWeb of Science
  77. 77.↵
    1. Tamura K,
    2. Stecher G,
    3. Peterson D,
    4. Filipski A,
    5. Kumar S
    . 2013. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30:2725–2729. doi:10.1093/molbev/mst197.
    OpenUrlCrossRefPubMedWeb of Science
  78. 78.↵
    1. Ultsch A,
    2. Mörchen F
    . 2005. ESOM-Maps: tools for clustering, visualization, and classification with Emergent SOM. University Marburg, Marburg, Germany.
  79. 79.↵
    1. Langmead B,
    2. Salzberg SL
    . 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357–359. doi:10.1038/nmeth.1923.
    OpenUrlCrossRefPubMedWeb of Science
  80. 80.↵
    1. Yutin N,
    2. Puigbò P,
    3. Koonin EV,
    4. Wolf YI
    . 2012. Phylogenomics of prokaryotic ribosomal proteins. PLoS One 7:e36972. doi:10.1371/journal.pone.0036972.
    OpenUrlCrossRefPubMed
  81. 81.↵
    1. Thiel V,
    2. Wood JM,
    3. Olsen MT,
    4. Tank M,
    5. Klatt CG,
    6. Ward DM,
    7. Bryant DA
    . 2016. The dark side of the mushroom spring microbial mat: life in the shadow of chlorophototrophs. I. Microbial diversity based on 16S rRNA gene amplicons and metagenomic sequencing. Front Microbiol 7:919. doi:10.3389/fmicb.2016.00919.
    OpenUrlCrossRef
  82. 82.↵
    1. Inskeep WP,
    2. Jay ZJ,
    3. Tringe SG,
    4. Herrgard M,
    5. Rusch DB
    , YNP Metagenome Project Steering Committee and Working Group Members. 2013. The YNP metagenome project: environmental parameters responsible for microbial distribution in the Yellowstone geothermal ecosystem. Front Microbiol 4:67. doi:10.3389/fmicb.2013.00067.
    OpenUrlCrossRefPubMed
  83. 83.↵
    1. Eloe-Fadrosh EA,
    2. Paez-Espino D,
    3. Jarett J,
    4. Dunfield PF,
    5. Hedlund BP,
    6. Dekas AE,
    7. Grasby SE,
    8. Brady AL,
    9. Dong H,
    10. Briggs BR,
    11. Li W-J,
    12. Goudeau D,
    13. Malmstrom R,
    14. Pati A,
    15. Pett-Ridge J,
    16. Rubin EM,
    17. Woyke T,
    18. Kyrpides NC,
    19. Ivanova NN
    . 2016. Global metagenomic survey reveals a new bacterial candidate phylum in geothermal springs. Nat Commun 7:10476. doi:10.1038/ncomms10476.
    OpenUrlCrossRefPubMed
  84. 84.↵
    1. Hoggins FE,
    2. Brooks RR
    . 1973. Natural dispersion of mercury from Puhipuhi, northland, New Zealand. N Z J Mar Freshwater Res 7:125–132. doi:10.1080/00288330.1973.9515459.
    OpenUrlCrossRef
  85. 85.↵
    1. Krabbenhoft D,
    2. Engstrom D,
    3. Gilmour C,
    4. Harris R,
    5. Hurley J,
    6. Mason R
    . 2007. Monitoring and evaluating trends in sediment and water indicators, p 47–87. In Murray MW, Saltman T, Harris R, Mason R, Reash R, Krabbenhoft DP (ed), Ecosystem responses to mercury contamination: indicators of change. CRC Press, Boca Raton, FL.
PreviousNext
Back to top
Download PDF
Citation Tools
Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs
Caitlin M. Gionfriddo, Matthew B. Stott, Jean F. Power, Jacob M. Ogorek, David P. Krabbenhoft, Ryan Wick, Kathryn Holt, Lin-Xing Chen, Brian C. Thomas, Jillian F. Banfield, John W. Moreau
Applied and Environmental Microbiology Jul 2020, 86 (15) e00176-20; DOI: 10.1128/AEM.00176-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs
Caitlin M. Gionfriddo, Matthew B. Stott, Jean F. Power, Jacob M. Ogorek, David P. Krabbenhoft, Ryan Wick, Kathryn Holt, Lin-Xing Chen, Brian C. Thomas, Jillian F. Banfield, John W. Moreau
Applied and Environmental Microbiology Jul 2020, 86 (15) e00176-20; DOI: 10.1128/AEM.00176-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

biogeochemistry
geothermal
hgcAB
MerA
mercuric ion detoxification
mercury
metagenomics
methylmercury

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336