AS3MT-mediated tolerance to arsenic evolved by multiple independent horizontal gene transfers from bacteria to eukaryotes

Lipophilic arsenic compounds in the cultured green alga Chlamydomonas reinhardtii

In this study, arsenic (As) speciation was investigated in the freshwater alga Chlamydomonas reinhardtii treated with 20 μg/L arsenate using fractionation as well as ICP-MS/ESI-MS analyses and was compared with the known As metabolite profile of wild-grown Saccharina latissima. While the total As accumulation in C. reinhardtii was about 85% lower than in S. latissima, the relative percentage of arsenolipids was significantly higher in C. reinhardtii (57.0% vs. 5.01%). As-containing hydrocarbons and phospholipids dominated the hydrophobic As profile in S. latissima, but no As-containing hydrocarbons were detectable in C. reinhardtii. Instead for the first time, an arsenoriboside-containing phytol (AsSugPhytol) was found to dominate the hydrophobic arsenicals of C. reinhardtii. Interestingly, this compound and its relatives had so far been only found in green marine microalgae, open sea plankton (mixed assemblage), and sediments but not in brown or red macroalgae. This compound family might therefore relate to differences in the arsenic metabolism between the algae phyla.

Genome-wide signatures of adaptation to extreme environments in red algae

The high temperature, acidity, and heavy metal-rich environments associated with hot springs have a major impact on biological processes in resident cells. One group of photosynthetic eukaryotes, the Cyanidiophyceae (Rhodophyta), has successfully thrived in hot springs and associated sites worldwide for more than 1 billion years. Here, we analyze chromosome-level assemblies from three representative Cyanidiophyceae species to study environmental adaptation at the genomic level. We find that subtelomeric gene duplication of functional genes and loss of canonical eukaryotic traits played a major role in environmental adaptation, in addition to horizontal gene transfer events. Shared responses to environmental stress exist in Cyanidiales and Galdieriales, however, most of the adaptive genes (e.g., for arsenic detoxification) evolved independently in these lineages. Our results underline the power of local selection to shape eukaryotic genomes that may face vastly different stresses in adjacent, extreme microhabitats.

The evolution of plant proton pump regulation via the R domain may have facilitated plant terrestrialization

Plasma membrane (PM) H⁺-ATPases are the electrogenic proton pumps that export H⁺ from plant and fungal cells to acidify the surroundings and generate a membrane potential. Plant PM H⁺-ATPases are equipped with a C‑terminal autoinhibitory regulatory (R) domain of about 100 amino acid residues, which could not be identified in the PM H⁺-ATPases of green algae but appeared fully developed in immediate streptophyte algal predecessors of land plants. To explore the physiological significance of this domain, we created in vivo C-terminal truncations of autoinhibited PM H⁺‑ATPase2 (AHA2), one of the two major isoforms in the land plant Arabidopsis thaliana. As more residues were deleted, the mutant plants became progressively more efficient in proton extrusion, concomitant with increased expansion growth and nutrient uptake. However, as the hyperactivated AHA2 also contributed to stomatal pore opening, which provides an exit pathway for water and an entrance pathway for pests, the mutant plants were more susceptible to biotic and abiotic stresses, pathogen invasion and water loss, respectively. Taken together, our results demonstrate that pump regulation through the R domain is crucial for land plant fitness and by controlling growth and nutrient uptake might have been necessary already for the successful water-to-land transition of plants.

Localized intensification of arsenic methylation within landfill leachate-saturated zone

Leachate-saturated zone (LSZ) of landfills is a complicated biogeochemical hotspot due to the continuous input of electron donors and acceptors from the top refuse layer with leachate migration. In this study, the methylation behavior of the arsenic (As) was investigated. The results indicate that As-methylation processes are influenced by temperature fields in LSZ. The dimethylarsinic acid biotransformation capability can be enhanced with an increase in temperature. Microbial diversity, quantification of functional gene (arsM), and co-occurrence network analysis further characterized the drivers of As methylation in LSZ. As-biogeochemical cycle pathways, as well as As-functional gene distribution among different temperature fields, were modeled on the basis of KEGG annotation. Binning analysis was further employed to assemble As-methylated metagenomes, enabling the identification of novel species for As methylation in landfills. Then, 87 high-quality draft metagenome-assembled genomes (MAGs) were reconstructed from LSZ refuse samples; nearly 15 % (13 of 87) belonged to putative As-methylates functional MAGs. Combined with the model of the As-biogeochemical cycle, nine putative functional species could complete methylation processes alone. The findings of this study highlighted the temperature influence on the As-methylation behavior in LSZ and could facilitate the management of As contamination in landfills.

Patterns of microbial arsenic detoxification genes in low-arsenic continental paddy soils

Microbes mediate the arsenic detoxification in paddy soils, determining the fate of arsenic in soils and its availability to rice plants, yet little is known about the structures and abundances of functional genes as well as the driving forces in low-arsenic paddy fields. To depict the arsenic detoxification functional gene patterns, 429 soil samples were collected from 39 paddy fields across four climatic zones in China, with the arsenic contents ranged from 9.76 to 19.74 mg kg^-1. GeoChip, a microarray-based metagenomic technique, was used to analyze the functional genes involved in arsenic detoxification. A total of three arsenic detoxification gene families were detected, aoxB, arxA (arsenite oxidase), and arsM (methyltransferase). Both the diversity and abundance of functional genes varied significantly among sampling sites (p < 0.05) and decreased along the arsenic gradient. Arsenic detoxification genes were carried by bacteria, archaea, and eukaryotes. Redundancy analysis showed that soil samples were grouped according to both climatic zones they located in and arsenic gradients at the continental scale. Soil pH, average annual temperature (AAT), arsenic, annual average precipitation (AAP), and CEC were the most important factors in shaping the functional structure. Structural equation modeling showed that AAT (r = 0.21), pH (r = -0.20), and arsenic contents (r = -0.11) directly affected the arsenic detoxification gene abundances. These findings provide an overall picture of microbial communities involved in arsenic detoxification in paddy soils and reveal the importance of climatic factors in shaping functional genes across a large spatial scale.

Synonymous single nucleotide polymorphism in arsenic (+3) methyltransferase of the Western mosquitofish (Gambusia affinis) and its gene expression among field populations

Naturally occurring arsenic is toxic at extremely low concentrations, yet some species persist even in high arsenic environments. We wanted to test if these species show evidence of evolution associated with arsenic exposure. To do this, we compared allelic variation across 872 coding nucleotides of arsenic (+3) methyltransferase (as3mt) and whole fish as3mt gene expression from three field populations of Gambusia affinis, from water sources containing low (1.9 ppb), medium-low (3.3 ppb), and high (15.7 ppb) levels of arsenic. The high arsenic site exceeds the US EPA's Maximum Contamination Level for drinking water. Medium-low and high populations exhibited homozygosity, and no sequence variation across all animals sampled. Eleven of 24 fish examined (45.8%) in the low arsenic population harbored synonymous single nucleotide polymorphisms (SNPs) in exons 4 and/or 10. SNP presence in the low arsenic population was not associated with differences in as3mt transcript levels compared to fish from the medium-low site, where SNPs were noted; however, as3mt expression in fish from the high arsenic concentration site was significantly lower than the other two sites. Low sequence variation in fish populations from sites with medium-low and high arsenic concentrations suggests greater selective pressure on this allele, while higher variation in the low population suggests a relaxed selection. Our results suggest gene regulation associated with arsenic detoxification may play a more crucial role in influencing responses to arsenic than polymorphic gene sequence. Understanding microevolutionary processes to various contaminants require the evaluation of multiple populations across a wide range of pollution exposures.

Arsenic methylation - Lessons from three decades of research

Between 1990 and 2020, our understanding of the significance of arsenic biomethylation changed in remarkable ways. At the beginning of this period, the conversion of inorganic arsenic into mono- and di-methylated metabolites was viewed primarily as a process that altered the kinetic behavior of arsenic. By increasing the rate of clearance of arsenic, the formation of methylated metabolites reduced exposure to this toxin; that is, methylation was detoxification. By 2020, it was clear that at least some of the toxic effects associated with As exposure depended on formation of methylated metabolites containing trivalent arsenic. Because the trivalent oxidation state of arsenic is associated with increased potency as a cytotoxin and clastogen, these findings were consistent with methylation-related changes in the dynamic behavior of arsenic. That is, methylation was activation. Our current understanding of the role of methylation as a modifier of kinetic and dynamic behaviors of arsenic is the product of research at molecular, cellular, organismic, and population levels. This information provides a basis for refining our estimates of risk associated with long term exposure to inorganic arsenic in environmental media, food, and water. This report summarizes the growth of our knowledge of enzymatically catalyzed methylation of arsenic over this period and considers the prospects for new discoveries.

Origins, fate, and actions of methylated trivalent metabolites of inorganic arsenic: progress and prospects

The toxic metalloid inorganic arsenic (iAs) is widely distributed in the environment. Chronic exposure to iAs from environmental sources has been linked to a variety of human diseases. Methylation of iAs is the primary pathway for metabolism of iAs. In humans, methylation of iAs is catalyzed by arsenic (+ 3 oxidation state) methyltransferase (AS3MT). Conversion of iAs to mono- and di-methylated species (MAs and DMAs) detoxifies iAs by increasing the rate of whole body clearance of arsenic. Interindividual differences in iAs metabolism play key roles in pathogenesis of and susceptibility to a range of disease outcomes associated with iAs exposure. These adverse health effects are in part associated with the production of methylated trivalent arsenic species, methylarsonous acid (MAs^III) and dimethylarsinous acid (DMAs^III), during AS3MT-catalyzed methylation of iAs. The formation of these metabolites activates iAs to unique forms that cause disease initiation and progression. Taken together, the current evidence suggests that methylation of iAs is a pathway for detoxification and for activation of the metalloid. Beyond this general understanding of the consequences of iAs methylation, many questions remain unanswered. Our knowledge of metabolic targets for MAs^III and DMAs^III in human cells and mechanisms for interactions between these arsenicals and targets is incomplete. Development of novel analytical methods for quantitation of MAs^III and DMAs^III in biological samples promises to address some of these gaps. Here, we summarize current knowledge of the enzymatic basis of MAs^III and DMAs^III formation, the toxic actions of these metabolites, and methods available for their detection and quantification in biomatrices. Major knowledge gaps and future research directions are also discussed.

Predicted AS3MT Proteins Methylate Arsenic and Support Two Major Phylogenetic AS3MT Groups

Inorganic arsenic is one of the most toxic and carcinogenic substances in the environment, but many organisms, including humans, methylate inorganic arsenic to mono-, di-, and trimethylated arsenic metabolites, which the organism can excrete. In humans and other eukaryotic organisms, the arsenite methyltransferase (AS3MT) protein methylates arsenite. AS3MT sequences from eukaryotic organisms group phylogenetically with predicted eubacterial AS3MT sequences, which has led to the suggestion that AS3MT was acquired from eubacteria by multiple events of horizontal gene transfer. In this study, we evaluated whether 55 (out of which 47 were predicted based on protein sequence similarity) sequences encoding putative AS3MT orthologues in 47 species from different kingdoms can indeed methylate arsenic. Fifty-three of the proteins showed arsenic methylating capacity. For example, the predicted AS3MT of the human gut bacterium Faecalibacterium prausnitzii methylated arsenic efficiently. We performed a kinetic analysis of 14 AS3MT proteins representing two phylogenetically distinct clades (Group 1 and 2) that each contain both eubacterial and eukaryotic sequences. We found that animal and bacterial AS3MTs in Group 1 rarely produce trimethylated arsenic, whereas Hydra vulgaris and the bacterium Rhodopseudomonas palustris in Group 2 produce trimethylated arsenic metabolites. These findings suggest that animals during evolution have acquired different arsenic methylating phenotypes from different bacteria. Further, it shows that humans carry two bacterial systems for arsenic methylation: one bacterium-derived AS3MT from Group 1 incorporated in the human genome and one from Group 2 in F. prausnitzii present in the gut microbiome.

Arsenic Methyltransferase and Methylation of Inorganic Arsenic

Arsenic occurs naturally in the environment, and exists predominantly as inorganic arsenite (As (III) and arsenate As (V)). Arsenic contamination of drinking water has long been recognized as a major global health concern. Arsenic exposure causes changes in skin color and lesions, and more severe health conditions such as black foot disease as well as various cancers originating in the lungs, skin, and bladder. In order to efficiently metabolize and excrete arsenic, it is methylated to monomethylarsonic and dimethylarsinic acid. One single enzyme, arsenic methyltransferase (AS3MT) is responsible for generating both metabolites. AS3MT has been purified from several mammalian and nonmammalian species, and its mRNA sequences were determined from amino acid sequences. With the advent of genome technology, mRNA sequences of AS3MT have been predicted from many species throughout the animal kingdom. Horizontal gene transfer had been postulated for this gene through phylogenetic studies, which suggests the importance of this gene in appropriately handling arsenic exposures in various organisms. An altered ability to methylate arsenic is dependent on specific single nucleotide polymorphisms (SNPs) in AS3MT. Reduced AS3MT activity resulting in poor metabolism of iAs has been shown to reduce expression of the tumor suppressor gene, p16, which is a potential pathway in arsenic carcinogenesis. Arsenic is also known to induce oxidative stress in cells. However, the presence of antioxidant response elements (AREs) in the promoter sequences of AS3MT in several species does not correlate with the ability to methylate arsenic. ARE elements are known to bind NRF2 and induce antioxidant enzymes to combat oxidative stress. NRF2 may be partly responsible for the biotransformation of iAs and the generation of methylated arsenic species via AS3MT. In this article, arsenic metabolism, excretion, and toxicity, a discussion of the AS3MT gene and its evolutionary history, and DNA methylation resulting from arsenic exposure have been reviewed.

Transcriptional and translational S-box riboswitches differ in ligand-binding properties

There are a number of riboswitches that utilize the same ligand-binding domain to regulate transcription or translation. S-box (SAM-I) riboswitches, including the riboswitch present in the Bacillus subtilis metI gene, which encodes cystathionine γ-synthase, regulate the expression of genes involved in methionine metabolism in response to SAM, primarily at the level of transcriptional attenuation. A rarer class of S-box riboswitches is predicted to regulate translation initiation. Here we identified and characterized a translational S-box riboswitch in the metI gene from Desulfurispirillum indicum The regulatory mechanisms of riboswitches are influenced by the kinetics of ligand interaction. The half-life of the translational D. indicum metI RNA-SAM complex is significantly shorter than that of the transcriptional B. subtilis metI RNA. This finding suggests that, unlike the transcriptional RNA, the translational metI riboswitch can make multiple reversible regulatory decisions. Comparison of both RNAs revealed that the second internal loop of helix P3 in the transcriptional RNA usually contains an A residue, whereas the translational RNA contains a C residue that is conserved in other S-box RNAs that are predicted to regulate translation. Mutational analysis indicated that the presence of an A or C residue correlates with RNA-SAM complex stability. Biochemical analyses indicate that the internal loop sequence critically determines the stability of the RNA-SAM complex by influencing the flexibility of residues involved in SAM binding and thereby affects the molecular mechanism of riboswitch function.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

Recent population genomic insights into the genetic basis of arsenic tolerance in humans: the difficulties of identifying positively selected loci in strongly bottlenecked populations

Recent advances in genomics have enabled researchers to shed light on the evolutionary processes driving human adaptation, by revealing the genetic architectures underlying traits ranging from lactase persistence, to skin pigmentation, to hypoxic response, to arsenic tolerance. Complicating the identification of targets of positive selection in modern human populations is their complex demographic history, characterized by population bottlenecks and expansions, population structure, migration, and admixture. In particular, founder effects and recent strong population size reductions, such as those experienced by the indigenous peoples of the Americas, have severe impacts on genetic variation that can lead to the accumulation of large allele frequency differences between populations due to genetic drift rather than natural selection. While distinguishing the effects of demographic history from selection remains challenging, neglecting neutral processes can lead to the incorrect identification of candidate loci. We here review the recent population genomic insights into the genetic basis of arsenic tolerance in Andean populations, and utilize this example to highlight both the difficulties pertaining to the identification of local adaptations in strongly bottlenecked populations, as well as the importance of controlling for demographic history in selection scans.

Functional evolution of nodulin 26-like intrinsic proteins: from bacterial arsenic detoxification to plant nutrient transport

Nodulin 26-like intrinsic proteins (NIPs) play essential roles in transporting the nutrients silicon and boron in seed plants, but the evolutionary origin of this transport function and the co-permeability to toxic arsenic remains enigmatic. Horizontal gene transfer of a yet uncharacterised bacterial AqpN-aquaporin group was the starting-point for plant NIP evolution. We combined intense sequence, phylogenetic and genetic context analyses and a mutational approach with various transport assays in oocytes and plants to resolve the transorganismal and functional evolution of bacterial and algal and terrestrial plant NIPs and to reveal their molecular transport specificity features. We discovered that aqpN genes are prevalently located in arsenic resistance operons of various prokaryotic phyla. We provided genetic and functional evidence that these proteins contribute to the arsenic detoxification machinery. We identified NIPs with the ancestral bacterial AqpN selectivity filter composition in algae, liverworts, moss, hornworts and ferns and demonstrated that these archetype plant NIPs and their prokaryotic progenitors are almost impermeable to water and silicon but transport arsenic and boron. With a mutational approach, we demonstrated that during evolution, ancestral NIP selectivity shifted to allow subfunctionalisations. Together, our data provided evidence that evolution converted bacterial arsenic efflux channels into essential seed plant nutrient transporters.

Roles of plasma membrane proton ATPases AHA2 and AHA7 in normal growth of roots and root hairs in Arabidopsis thaliana

Plasma membrane H⁺ -ATPase pumps build up the electrochemical H⁺ gradients that energize most other transport processes into and out of plant cells through channel proteins and secondary active carriers. In Arabidopsis thaliana, the AUTOINHIBITED PLASMA MEMBRANE H⁺ -ATPases AHA1, AHA2 and AHA7 are predominant in root epidermal cells. In contrast to other H⁺ -ATPases, we find that AHA7 is autoinhibited by a sequence present in the extracellular loop between transmembrane segments 7 and 8. Autoinhibition of pump activity was regulated by extracellular pH, suggesting negative feedback regulation of AHA7 during establishment of an H⁺ gradient. Due to genetic redundancy, it has proven difficult to test the role of AHA2 and AHA7, and mutant phenotypes have previously only been observed under nutrient stress conditions. Here, we investigated root and root hair growth under normal conditions in single and double mutants of AHA2 and AHA7. We find that AHA2 drives root cell expansion during growth but that, unexpectedly, restriction of root hair elongation is dependent on AHA2 and AHA7, with each having different roles in this process.

A global survey of arsenic-related genes in soil microbiomes

Background

Environmental resistomes include transferable microbial genes. One important resistome component is resistance to arsenic, a ubiquitous and toxic metalloid that can have negative and chronic consequences for human and animal health. The distribution of arsenic resistance and metabolism genes in the environment is not well understood. However, microbial communities and their resistomes mediate key transformations of arsenic that are expected to impact both biogeochemistry and local toxicity.

Results

We examined the phylogenetic diversity, genomic location (chromosome or plasmid), and biogeography of arsenic resistance and metabolism genes in 922 soil genomes and 38 metagenomes. To do so, we developed a bioinformatic toolkit that includes BLAST databases, hidden Markov models and resources for gene-targeted assembly of nine arsenic resistance and metabolism genes: acr3, aioA, arsB, arsC (grx), arsC (trx), arsD, arsM, arrA, and arxA. Though arsenic-related genes were common, they were not universally detected, contradicting the common conjecture that all organisms have them. From major clades of arsenic-related genes, we inferred their potential for horizontal and vertical transfer. Different types and proportions of genes were detected across soils, suggesting microbial community composition will, in part, determine local arsenic toxicity and biogeochemistry. While arsenic-related genes were globally distributed, particular sequence variants were highly endemic (e.g., acr3), suggesting dispersal limitation. The gene encoding arsenic methylase arsM was unexpectedly abundant in soil metagenomes (median 48%), suggesting that it plays a prominent role in global arsenic biogeochemistry.

Conclusions

Our analysis advances understanding of arsenic resistance, metabolism, and biogeochemistry, and our approach provides a roadmap for the ecological investigation of environmental resistomes.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0661-5) contains supplementary material, which is available to authorized users.

A methyltransferase gene from arbuscular mycorrhizal fungi involved in arsenic methylation and volatilization

Arbuscular mycorrhizal fungi (AMF), ubiquitous symbiotic fungi associated with the majority of terrestrial plants, were demonstrated to play important roles in arsenic (As) translocation and transformation in the plant-soil continuum, and substantially influence plant As tolerance. However, the direct involvement of AMF in As methylation and volatilization and their molecular mechanisms remain unsolved. Here, an arsenite methyltransferase gene RiMT-11 was identified and characterized from AM fungus Rhizophagus irregularis. Heterologous expression of RiMT-11 enhanced arsenite resistance of E. coli (Δars) through methylating As into monomethylarsonic acid (MMA), dimethylarsinic acid (DMA) and ultimately volatile trimethyl arsine (TMAs). In a two-compartment in vitro monoxenic cultivation system, methylated and volatile As were also detected from AM symbioses with arsenate addition, accompanied by strong up-regulation of RiMT-11 expression in extraradical hyphae. The present study provided direct evidence and illustrated an underlying mechanism of As methylation and volatilization by AMF, leading to a deeper insight into the role of AMF in As biogeochemical cycling.

Recurrent horizontal transfer of arsenite methyltransferase genes facilitated adaptation of life to arsenic

The toxic metalloid arsenic has been environmentally ubiquitous since life first arose nearly four billion years ago and presents a challenge for the survival of all living organisms. Its bioavailability has varied dramatically over the history of life on Earth. As life spread, biogeochemical and climate changes cyclically increased and decreased bioavailable arsenic. To elucidate the history of arsenic adaptation across the tree of life, we reconstructed the phylogeny of the arsM gene that encodes the As(III) S-adenosylmethionine (SAM) methyltransferase. Our results suggest that life successfully moved into arsenic-rich environments in the late Archean Eon and Proterozoic Eon, respectively, by the spread of arsM genes. The arsM genes of bacterial origin have been transferred to other kingdoms of life on at least six occasions, and the resulting domesticated arsM genes promoted adaptation to environmental arsenic. These results allow us to peer into the history of arsenic adaptation of life on our planet and imply that dissemination of genes encoding diverse adaptive functions to toxic chemicals permit adaptation to changes in concentrations of environmental toxins over evolutionary history.