The gene pair hgcAB is essential for microbial mercury methylation. Our understanding of its abundance and diversity in nature is rapidly evolving. In this study we developed a new broad-range primer set for hgcAB , plus an expanded hgcAB reference library, and used these to characterize Hg-methylating communities from diverse environments. We applied this new Hg-methylator database to assign taxonomy to hgcA sequences from clone, amplicon, and metagenomic datasets. We evaluated potential biases introduced in primer design, sequence length, and classification, and suggest best practices for studying Hg-methylator diversity. Our study confirms the emerging picture of an expanded diversity of HgcAB-encoding microbes in many types of ecosystems, with abundant putative mercury methylators Nitrospirae and Chloroflexi in several new environments including salt marsh and peat soils. Other common microbes encoding HgcAB included Phycisphaerae , Aminicenantes , Spirochaetes , and Elusimicrobia. Combined with high-throughput amplicon specific sequencing, the new primer set also indentified novel hgcAB sequences similar to Lentisphaerae , Bacteroidetes , Atribacteria , and candidate phyla WOR-3 and KSB1 bacteria. Gene abundance data also corroborate the important role of two “classic” groups of methylators ( Deltaproteobacteria and Methanomicrobia ) in many environments, but generally show a scarcity of hgcAB + Firmicutes . The new primer set was developed to specifically target hgcAB sequences found in nature, reducing degeneracy and providing increased sensitivity while maintaining broad diversity capture. We evaluated mock communities to confirm primer improvements, including culture spikes to environmental samples with variable DNA extraction and PCR amplification efficiencies. For select sites, this new workflow was combined with direct high-throughput hgcAB sequencing. The hgcAB diversity generated by direct amplicon sequencing confirmed the potential for novel Hg-methylators previously identified using metagenomic screens. A new phylogenetic analysis using sequences from freshwater, saline, and terrestrial environments showed Deltaproteobacteria HgcA sequences generally clustered among themselves, while metagenome-resolved HgcA sequences in other phyla tended to cluster by environment, suggesting horizontal gene transfer into many clades. HgcA from marine metagenomes often formed distinct subtrees from those sequenced from freshwater ecosystems. Overall the majority of HgcA sequences branch from a cluster of HgcAB ...

Metagenomic studies on geothermal environments have been central in recent discoveries on the diversity of archaeal methane and alkane metabolism. Here, we investigated methanogenic populations inhabiting terrestrial geothermal features in Yellowstone National Park (YNP) by combining amplicon sequencing with metagenomics and mesocosm experiments. Detection of methyl-coenzyme M reductase subunit A (mcrA) gene amplicons demonstrated a wide diversity of Mcr-encoding archaea inhabit geothermal features with differing physicochemical regimes across YNP. From three selected hot springs we recovered twelve Mcr-encoding metagenome assembled genomes (MAGs) affiliated with lineages of cultured methanogens as well as Candidatus (Ca.) Methanomethylicia, Ca. Hadesarchaeia, and Archaeoglobi. These MAGs encoded the potential for hydrogenotrophic, aceticlastic, hydrogen-dependent methylotrophic methanogenesis, or anaerobic short-chain alkane oxidation. While Mcr-encoding archaea represent minor fractions of the microbial community of hot springs, mesocosm experiments with methanogenic precursors resulted in the stimulation of methanogenic activity and the enrichment of lineages affiliated with Methanosaeta and Methanothermobacter as well as with uncultured Mcr-encoding archaea including Ca. Korarchaeia, Ca. Nezhaarchaeia, and Archaeoglobi. We revealed that diverse Mcr-encoding archaea with the metabolic potential to produce methane from different precursors persist in the geothermal environments of YNP and can be enriched under methanogenic conditions. This study highlights the importance of combining environmental metagenomics with laboratory-based experiments to expand our understanding of uncultured Mcr-encoding archaea and their potential impact on microbial carbon transformations in geothermal environments and beyond.

Metagenomic studies on geothermal environments have been central in recent discoveries on the diversity of archaeal methane and alkane metabolism. Here, we investigated the methanogenic populations inhabiting terrestrial geothermal features in Yellowstone National Park (YNP) by combining amplicon sequencing with metagenomics and mesocosm experiments. Detection of gene amplicons of methyl-coenzyme M reductase subunit A (mcrA) indicated a wide diversity of Mcr-encoding archaea across geothermal features with differing physicochemical regimes. From three selected hot springs we recovered twelve Mcr-encoding metagenome assembled genomes (MAGs) affiliated with lineages of cultured methanogens as well as Candidatus (Ca.) Methanomethylicia, Ca. Hadesarchaeia, and Archaeoglobi. These MAGs encoded the potential for hydrogenotrophic, aceticlastic, or hydrogen-dependent methylotrophic methanogenesis, or anaerobic short-chain alkane oxidation. While Mcr-encoding archaea represented a minor fraction of the microbial community of hot springs, mesocosm experiments with methanogenic precursors resulted in stimulation of methanogenic activity and the enrichment of lineages affiliated with Methanosaeta and Methanothermobacter as well as with uncultured Mcr-encoding archaea including Ca. Korarchaeia, Ca. Nezhaarchaeia, and Archaeoglobi. Altogether, we revealed that diverse Mcr-encoding populations with the metabolic potential to produce methane from different precursors persist in the geothermal environments of YNP. This study highlights the importance of combining environmental metagenomics with laboratory-based experiments to expand our understanding of uncultured Mcr-encoding archaea and their potential impact on microbial carbon transformations in geothermal environments and beyond.

The gene pair hgcAB is essential for microbial mercury methylation. Our understanding of its abundance and diversity in nature is rapidly evolving. In this study we developed a new broadrange primer set for hgcAB, plus an expanded hgcAB reference library, and used these to characterize Hg-methylating communities from diverse environments. We applied this new Hgmethylator database to assign taxonomy to hgcA sequences from clone, amplicon, and metagenomic datasets. We evaluated potential biases introduced in primer design, sequence length, and classification, and suggest best practices for studying Hg-methylator diversity. Our study confirms the emerging picture of an expanded diversity of HgcAB-encoding microbes in many types of ecosystems, with abundant putative mercury methylators Nitrospirae and Chloroflexi in several new environments including salt marsh and peat soils. Other common microbes encoding HgcAB included Phycisphaerae, Aminicenantes, Spirochaetes, and Elusimicrobia. Gene abundance data also corroborate the important role of two "classic" groups of methylators (Deltaproteobacteria and Methanomicrobia) in many environments, but generally show a scarcity of hgcAB+ Firmicutes. The new primer set was developed to specifically target hgcAB sequences found in nature, reducing degeneracy and providing increased sensitivity while maintaining broad diversity capture. We evaluated mock communities to confirm primer improvements, including culture spikes to environmental samples with variable DNA extraction and PCR amplification efficiencies. For select sites, this new workflow was combined with direct high-throughput hgcAB sequencing. The hgcAB diversity generated by direct amplicon sequencing confirmed the potential for novel Hg-methylators previously identified using metagenomic screens. A new phylogenetic analysis using sequences from freshwater, saline, and terrestrial environments showed Deltaproteobacteria HgcA sequences generally clustered among themselves, while metagenome-resolved HgcA sequences in other phyla tended to cluster by environment, suggesting horizontal gene transfer into many clades. HgcA from marine metagenomes often formed distinct subtrees from those sequenced from freshwater ecosystems. Overall the majority of HgcA sequences branch from a cluster of HgcAB fused proteins related to Thermococci, Atribacteria (candidate division OP9), Aminicenantes (OP8), and Chloroflexi. The improved primer set and library, combined with direct amplicon sequencing, provide a significantly improved assessment of the abundance and diversity of hgcAB+ microbes in nature. Contribution to the Field StatementThe gene pair hgcAB is essential for microbial production of the neurotoxin methylmercury. In recent years these genes have been used as biomarkers to determine the potential of a microbiome to generate methylmercury via PCR amplification using degenerate primers from several research groups. However, improved techniques for capturing hgcAB diversity are necessary for identifying the major environmental p...

Geothermal environments, including terrestrial hot springs and deep-sea hydrothermal sediments, often contain many poorly understood lineages of archaea. Here, we recovered ten metagenome-assembled genomes (MAGs) from geothermal sediments and propose that they constitute a new archaeal class within the TACK superphylum, “Candidatus Culexarchaeia”, named after the Culex Basin in Yellowstone National Park. Culexarchaeia harbor distinct sets of proteins involved in key cellular processes that are either phylogenetically divergent or are absent from other closely related TACK lineages, with a particular divergence in cell division and cytoskeletal proteins. Metabolic reconstruction revealed that Culexarchaeia have the capacity to metabolize a wide variety of organic and inorganic substrates. Notably, Culexarchaeia encode a unique modular, membrane associated, and energy conserving [NiFe]-hydrogenase complex that potentially interacts with heterodisulfide reductase (Hdr) subunits. Comparison of this [NiFe]-hydrogenase complex with similar complexes from other archaea suggests that interactions between membrane associated [NiFe]-hydrogenases and Hdr may be more widespread than previously appreciated in both methanogenic and non-methanogenic lifestyles. The analysis of Culexarchaeia further expands our understanding of the phylogenetic and functional diversity of lineages within the TACK superphylum and the ecology, physiology, and evolution of these organisms in extreme environments.