Genome insights of mercury methylation among Desulfovibrio and Pseudodesulfovibrio strains

Recent advance of microbial mercury methylation in the environment

Methylmercury formation is mainly driven by microbial-mediated process. The mechanism of microbial mercury methylation has become a crucial research topic for understanding methylation in the environment. Pioneering studies of microbial mercury methylation are focusing on functional strain isolation, microbial community composition characterization, and mechanism elucidation in various environments. Therefore, the functional genes of microbial mercury methylation, global isolations of Hg methylation strains, and their methylation potential were systematically analyzed, and methylators in typical environments were extensively reviewed. The main drivers (key physicochemical factors and microbiota) of microbial mercury methylation were summarized and discussed. Though significant progress on the mechanism of the Hg microbial methylation has been explored in recent decade, it is still limited in several aspects, including (1) molecular biology techniques for identifying methylators; (2) characterization methods for mercury methylation potential; and (3) complex environmental properties (environmental factors, complex communities, etc.). Accordingly, strategies for studying the Hg microbial methylation mechanism were proposed. These strategies include the following: (1) the development of new molecular biology methods to characterize methylation potential; (2) treating the environment as a micro-ecosystem and studying them from a holistic perspective to clearly understand mercury methylation; (3) a more reasonable and sensitive inhibition test needs to be considered. KEY POINTS: • Global Hg microbial methylation is phylogenetically and functionally discussed. • The main drivers of microbial methylation are compared in various condition. • Future study of Hg microbial methylation is proposed.

Metabolically diverse microorganisms mediate methylmercury formation under nitrate-reducing conditions in a dynamic hydroelectric reservoir

Brownlee Reservoir is a mercury (Hg)-impaired hydroelectric reservoir that exhibits dynamic hydrological and geochemical conditions and is located within the Hells Canyon Complex in Idaho, USA. Methylmercury (MeHg) contamination in fish is a concern in the reservoir. While MeHg production has historically been attributed to sulfate-reducing bacteria and methanogenic archaea, microorganisms carrying the hgcA gene are taxonomically and metabolically diverse and the major biogeochemical cycles driving mercury (Hg) methylation are not well understood. In this study, Hg speciation and redox-active compounds were measured throughout Brownlee Reservoir across the stratified period in four consecutive years (2016-2019) to identify the location where and redox conditions under which MeHg is produced. Metagenomic sequencing was performed on a subset of samples to characterize the microbial community with hgcA and identify possible links between biogeochemical cycles and MeHg production. Biogeochemical profiles suggested in situ water column Hg methylation was the major source of MeHg. These profiles, combined with genome-resolved metagenomics focused on hgcA-carrying microbes, indicated that MeHg production occurs in this system under nitrate- or manganese-reducing conditions, which were previously thought to preclude Hg-methylation. Using this multidisciplinary approach, we identified the cascading effects of interannual variability in hydrology on the redox status, microbial metabolic strategies, abundance and metabolic diversity of Hg methylators, and ultimately MeHg concentrations throughout the reservoir. This work expands the known conditions conducive to producing MeHg and suggests that the Hg-methylation mitigation efforts by nitrate or manganese amendment may be unsuccessful in some locations.

Mining-impacted rice paddies select for Archaeal methylators and reveal a putative (Archaeal) regulator of mercury methylation

Methylmercury (MeHg) is a microbially produced neurotoxin derived from inorganic mercury (Hg), which accumulation in rice represents a major health concern to humans. However, the microbial control of MeHg dynamics in the environment remains elusive. Here, leveraging three rice paddy fields with distinct concentrations of Hg (Total Hg (THg): 0.21-513 mg kg-1 dry wt. soil; MeHg: 1.21-6.82 ng g-1 dry wt. soil), we resorted to metagenomics to determine the microbial determinants involved in MeHg production under contrasted contamination settings. We show that Hg methylating Archaea, along with methane-cycling genes, were enriched in severely contaminated paddy soils. Metagenome-resolved Genomes of novel putative Hg methylators belonging to Nitrospinota (UBA7883), with poorly resolved taxonomy despite high completeness, showed evidence of facultative anaerobic metabolism and adaptations to fluctuating redox potential. Furthermore, we found evidence of environmental filtering effects that influenced the phylogenies of not only hgcA genes under different THg concentrations, but also of two housekeeping genes, rpoB and glnA, highlighting the need for further experimental validation of whether THg drives the evolution of hgcAB. Finally, assessment of the genomic environment surrounding hgcAB suggests that this gene pair may be regulated by an archaeal toxin-antitoxin (TA) system, instead of the more frequently found arsR-like genes in bacterial methylators. This suggests the presence of distinct hgcAB regulation systems in bacteria and archaea. Our results support the emerging role of Archaea in MeHg cycling under mining-impacted environments and shed light on the differential control of the expression of genes involved in MeHg formation between Archaea and Bacteria.

Effect of exogenous and endogenous sulfide on the production and the export of methylmercury by sulfate-reducing bacteria

Mercury (Hg) is a global pollutant of environmental and health concern; its methylated form, methylmercury (MeHg), is a potent neurotoxin. Sulfur-containing molecules play a role in MeHg production by microorganisms. While sulfides are considered to limit Hg methylation, sulfate and cysteine were shown to favor this process. However, these two forms can be endogenously converted by microorganisms into sulfide. Here, we explore the effect of sulfide (produced by the cell or supplied exogenously) on Hg methylation. For this purpose, Pseudodesulfovibrio hydrargyri BerOc1 was cultivated in non-sulfidogenic conditions with addition of cysteine and sulfide as well as in sulfidogenic conditions. We report that Hg methylation depends on sulfide concentration in the culture and the sulfides produced by cysteine degradation or sulfate reduction could affect the Hg methylation pattern. Hg methylation was independent of hgcA expression. Interestingly, MeHg production was maximal at 0.1-0.5 mM of sulfides. Besides, a strong positive correlation between MeHg in the extracellular medium and the increase of sulfide concentrations was observed, suggesting a facilitated MeHg export with sulfide and/or higher desorption from the cell. We suggest that sulfides (exogenous or endogenous) play a key role in controlling mercury methylation and should be considered when investigating the impact of Hg in natural environments.

Roles of plant-associated microorganisms in regulating the fate of Hg in croplands: A perspective on potential pathways in maintaining sustainable agriculture

In heavy metal-contaminated croplands, plant-associated microorganisms play important roles in the adaptation of crops to heavy metals. Plant-associated microbes can interact with Hg and stimulate plant resistance to Hg toxicity, which is crucial for impeding Hg accumulation along the food chain. The roles of rhizosphere microorganisms for the improvement of plant growth and Hg resistance have drawn great research attention. However, the interactions among plant-endophyte-Hg have been neglected although they might be important for in vivo Hg detoxification. In this study, we systematically summarized 1) the roles of plant-associated microorganisms in Hg detoxification and plant growth, 2) Hg methylation and demethylation driven by plant-associated microbes, 3) the relationships between plant-associated microbes and Hg biogeochemical cycling. The possible mechanisms underlying crop-endophyte-Hg interactions were discussed, although limited studies on this aspect are available to date. The challenges and perspectives of plant-endophytes in dampening Hg phytotoxicity and controlling Hg accumulation in croplands were proposed on the basis of the present knowledge. Taken together, this work provides evidence for further understanding the interactions between soil-plant-endophyte-Hg systems and as well as new interpretations and perspectives into regulating the fate of Hg in croplands.

Orphan gene in Littorina: An unexpected role of symbionts in the host evolution

Mechanisms of reproductive isolation between closely related sympatric species are of high evolutionary significance as they may function as initial drivers of speciation and protect species integrity afterwards. Proteins involved in the establishment of reproductive barriers often evolve fast and may be key players in cessation of gene flow between the incipient species. The five Atlantic Littorina (Neritrema) species represent a notable example of recent radiation. The geographic ranges of these young species largely overlap and the mechanisms of reproductive isolation are poorly understood. In this study, we performed a detailed analysis of the reproductive protein LOSP, previously identified in Littorina. We showed that this protein is evolutionary young and taxonomically restricted to the genus Littorina. It has high sequence variation both within and between Littorina species, which is compatible with its presumable role in the reproductive isolation. The strongest differences in the LOSP structure were detected between Littorina subgenera with distinctive repetitive motifs present exclusively in the Neritrema species, but not in L. littorea. Moreover, the sequence of these repetitive structural elements demonstrates a high homology with genetic elements of bacteria, identified as components of Littorina associated microbiomes. We suggest that these elements were acquired from a symbiotic bacterial donor via horizontal genetic transfer (HGT), which is indirectly confirmed by the presence of multiple transposable elements in the LOSP flanking and intronic regions. Furthermore, we hypothesize that this HGT-driven evolutionary innovation promoted LOSP function in reproductive isolation, which might be one of the factors determining the intensive cladogenesis in the Littorina (Neritrema) lineage in contrast to the anagenesis in the L. littorea clade.

Establishment of Genome Based Criteria for Classification of the Family Desulfovibrionaceae and Proposal of Two Novel Genera, Alkalidesulfovibrio gen. nov. and Salidesulfovibrio gen. nov

Bacteria in the Desulfovibrionaceae family, which contribute to S element turnover as sulfate-reducing bacteria (SRB) and disproportionation of partially oxidized sulfoxy anions, have been extensively investigated since the importance of the sulfur cycle emerged. Novel species belonging to this taxon are frequently reported, because they exist in various environments and are easy to culture using established methods. Due to the rapid expansion of the taxon, correction and reclassification have been conducted. The development of high-throughput sequencing facilitated rapid expansion of genome sequence database. Genome-based criteria, based on these databases, proved to be potential classification standard by overcoming the limitations of 16S rRNA-based phylogeny. Although standards methods for taxogenomics are being established, the addition of a novel genus requires extensive calculations with taxa, including many species, such as Desulfovibrionaceae. Thus, the genome-based criteria for classification of Desulfovibrionaceae were established and validated in this study. The average amino-acid identity (AAI) cut-off value, 63.43 ± 0.01, was calculated to be an appropriate criterion for genus delineation of the family Desulfovibrionaceae. By applying the AAI cut-off value, 88 genomes of the Desulfovibrionaceae were divided into 27 genera, which follows the core gene phylogeny results. In this process, two novel genera (Alkalidesulfovibrio and Salidesulfovibrio) and one former invalid genus (“Psychrodesulfovibrio”) were officially proposed. Further, by applying the 95–96% average nucleotide identity (ANI) standard and the 70% digital DNA–DNA hybridization standard values for species delineation of strains that were classified as the same species, five strains have the potential to be newly classified. After verifying that the classification was appropriately performed through relative synonymous codon usage analysis, common characteristics were listed by group. In addition, by detecting metal resistance related genes via in silico analysis, it was confirmed that most strains display metal tolerance.

Bioinformatics Investigations of Universal Stress Proteins from Mercury-Methylating Desulfovibrionaceae

The presence of methylmercury in aquatic environments and marine food sources is of global concern. The chemical reaction for the addition of a methyl group to inorganic mercury occurs in diverse bacterial taxonomic groups including the Gram-negative, sulfate-reducing Desulfovibrionaceae family that inhabit extreme aquatic environments. The availability of whole-genome sequence datasets for members of the Desulfovibrionaceae presents opportunities to understand the microbial mechanisms that contribute to methylmercury production in extreme aquatic environments. We have applied bioinformatics resources and developed visual analytics resources to categorize a collection of 719 putative universal stress protein (USP) sequences predicted from 93 genomes of Desulfovibrionaceae. We have focused our bioinformatics investigations on protein sequence analytics by developing interactive visualizations to categorize Desulfovibrionaceae universal stress proteins by protein domain composition and functionally important amino acids. We identified 651 Desulfovibrionaceae universal stress protein sequences, of which 488 sequences had only one USP domain and 163 had two USP domains. The 488 single USP domain sequences were further categorized into 340 sequences with ATP-binding motif and 148 sequences without ATP-binding motif. The 163 double USP domain sequences were categorized into (1) both USP domains with ATP-binding motif (3 sequences); (2) both USP domains without ATP-binding motif (138 sequences); and (3) one USP domain with ATP-binding motif (21 sequences). We developed visual analytics resources to facilitate the investigation of these categories of datasets in the presence or absence of the mercury-methylating gene pair (hgcAB). Future research could utilize these functional categories to investigate the participation of universal stress proteins in the bacterial cellular uptake of inorganic mercury and methylmercury production, especially in anaerobic aquatic environments.

Pseudodesulfovibrio mercurii sp. nov., a mercury-methylating bacterium isolated from sediment

The sulfate-reducing, mercury-methylating strain ND132^T was isolated from the brackish anaerobic bottom sediments of Chesapeake Bay, USA. Capable of high levels of mercury (Hg) methylation, ND132^T has been widely used as a model strain to study the process and to determine the genetic basis of Hg methylation. Originally called Desulfovibrio desulfuricans ND132^T on the basis of an early partial 16S rRNA sequence, the strain has never been formally described. Phylogenetic and physiological traits place this strain within the genus Pseudodesulfovibrio, in the recently reclassified phylum Desulfobacterota (formerly Deltaproteobacteria). ND132^T is most closely related to Pseudodesulfovibrio hydrargyri BerOc1^T and Pseudodesulfovibrio indicus J2^T. Analysis of average nucleotide identity (ANI) of whole-genome sequences showed roughly 88 % ANI between P. hydrargyri BerOc1^T and ND132^T, and 84 % similarity between ND132^T and P. indicus J2^T. These cut-off scores <95 %, along with a multi-gene phylogenetic analysis of members of the family Desulfovibrionacea, and differences in physiology indicate that all three strains represent separate species. The Gram-stain-negative cells are vibrio-shaped, motile and not sporulated. ND132^T is a salt-tolerant mesophile with optimal growth in the laboratory at 32 °C, 2 % salinity, and pH 7.8. The DNA G+C content of the genomic DNA is 65.2 %. It is an incomplete oxidizer of short chain fatty acids, using lactate, pyruvate and fumarate with sulfate or sulfite as the terminal electron acceptors. ND132^T can respire fumarate using pyruvate as an electron donor. The major fatty acids are iso-C_15 : 0, anteiso-C_15 : 0, iso-C_17 : 0, iso-C_17 : 1ω9c and anteiso-C_17 : 0. We propose the classification of strain ND132^T (DSM 110689, ATCC TSD-224) as the type strain Pseudodesulfovibrio mercurii sp. nov.

Mercury Methylation Genes Identified across Diverse Anaerobic Microbial Guilds in a Eutrophic Sulfate-Enriched Lake

Mercury (Hg) methylation is a microbially mediated process that converts inorganic Hg into bioaccumulative, neurotoxic methylmercury (MeHg). The metabolic activity of methylating organisms is highly dependent on biogeochemical conditions, which subsequently influences MeHg production. However, our understanding of the ecophysiology of methylators in natural ecosystems is still limited. Here, we identified potential locations of MeHg production in the anoxic, sulfidic hypolimnion of a freshwater lake. At these sites, we used shotgun metagenomics to characterize microorganisms with the Hg-methylation gene hgcA. Putative methylators were dominated by hgcA sequences divergent from those in well-studied, confirmed methylators. Using genome-resolved metagenomics, we identified organisms with hgcA (hgcA+) within the Bacteroidetes and the recently described Kiritimatiellaeota phyla. We identified hgcA+ genomes derived from sulfate-reducing bacteria, but these accounted for only 22% of hgcA+ genome coverage. The most abundant hgcA+ genomes were from fermenters, accounting for over half of the hgcA gene coverage. Many of these organisms also mediate hydrolysis of polysaccharides, likely from cyanobacterial blooms. This work highlights the distribution of the Hg-methylation genes across microbial metabolic guilds and indicate that primary degradation of polysaccharides and fermentation may play an important but unrecognized role in MeHg production in the anoxic hypolimnion of freshwater lakes.

An Improved hgcAB Primer Set and Direct High-Throughput Sequencing Expand Hg-Methylator Diversity in Nature

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 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.

Expanded Phylogenetic Diversity and Metabolic Flexibility of Mercury-Methylating Microorganisms

Methylmercury is a potent bioaccumulating neurotoxin that is produced by specific microorganisms that methylate inorganic mercury. Methylmercury production in diverse anaerobic bacteria and archaea was recently linked to the hgcAB genes. However, the full phylogenetic and metabolic diversity of mercury-methylating microorganisms has not been fully unraveled due to the limited number of cultured experimentally verified methylators and the limitations of primer-based molecular methods. Here, we describe the phylogenetic diversity and metabolic flexibility of putative mercury-methylating microorganisms by hgcAB identification in publicly available isolate genomes and metagenome-assembled genomes (MAGs) as well as novel freshwater MAGs. We demonstrate that putative mercury methylators are much more phylogenetically diverse than previously known and that hgcAB distribution among genomes is most likely due to several independent horizontal gene transfer events. The microorganisms we identified possess diverse metabolic capabilities spanning carbon fixation, sulfate reduction, nitrogen fixation, and metal resistance pathways. We identified 111 putative mercury methylators in a set of previously published permafrost metatranscriptomes and demonstrated that different methylating taxa may contribute to hgcA expression at different depths. Overall, we provide a framework for illuminating the microbial basis of mercury methylation using genome-resolved metagenomics and metatranscriptomics to identify putative methylators based upon hgcAB presence and describe their putative functions in the environment.IMPORTANCE Accurately assessing the production of bioaccumulative neurotoxic methylmercury by characterizing the phylogenetic diversity, metabolic functions, and activity of methylators in the environment is crucial for understanding constraints on the mercury cycle. Much of our understanding of methylmercury production is based on cultured anaerobic microorganisms within the Deltaproteobacteria, Firmicutes, and Euryarchaeota. Advances in next-generation sequencing technologies have enabled large-scale cultivation-independent surveys of diverse and poorly characterized microorganisms from numerous ecosystems. We used genome-resolved metagenomics and metatranscriptomics to highlight the vast phylogenetic and metabolic diversity of putative mercury methylators and their depth-discrete activities in thawing permafrost. This work underscores the importance of using genome-resolved metagenomics to survey specific putative methylating populations of a given mercury-impacted ecosystem.