Parapatric speciation of Meiothermus in serpentinite-hosted aquifers in Oman

The factors that control the distribution and evolution of microbial life in subsurface environments remain enigmatic due to challenges associated with sampling fluids from discrete depth intervals via boreholes while avoiding mixing of fluids. Here, using an inflatable packer system, fracture waters were isolated and collected from three discrete depth intervals spanning >130 m in a borehole intersecting an ultramafic rock formation undergoing serpentinization in the Samail Ophiolite, Sultanate of Oman. Near surface aquifer waters were moderately reducing and had alkaline pH while deeper aquifer waters were reduced and had hyperalkaline pH, indicating extensive influence by serpentinization. Metagenomic sequencing and analysis of DNA from filtered biomass collected from discrete depth intervals revealed an abundance of aerobes in near surface waters and a greater proportion of anaerobes at depth. Yet the abundance of the putatively obligate aerobe, Meiothermus, increased with depth, providing an opportunity to evaluate the influence of chemical and spatial variation on its distribution and speciation. Two clades of Meiothermus metagenome assembled genomes (MAGs) were identified that correspond to surface and deep populations termed Types I (S) and II (D), respectively; both clades comprised an apparently Oman-specific lineage indicating a common ancestor. Type II (D) clade MAGs encoded fewer genes and were undergoing slower genome replication as inferred from read mapping. Further, single nucleotide variants (SNVs) and mobile genetic elements identified among MAGs revealed detectable, albeit limited, evidence for gene flow/recombination between spatially segregated Type I (S) and Type II (D) populations. Together, these observations indicate that chemical variation generated by serpentinization, combined with physical barriers that reduce/limit dispersal and gene flow, allowed for the parapatric speciation of Meiothermus in the Samail Ophiolite or a geologic precursor. Further, Meiothermus genomic data suggest that deep and shallow aquifer fluids in the Samail Ophiolite may mix over shorter time scales than has been previously estimated from geochemical data

Parapatric speciation of <i>Meiothermus</i> in serpentinite-hosted aquifers in Oman

The factors that control the distribution and evolution of microbial life in subsurface environments remain enigmatic due to challenges associated with sampling fluids from discrete depth intervals via boreholes while avoiding mixing of fluids. Here, using an inflatable packer system, fracture waters were isolated and collected from three discrete depth intervals spanning &gt;130 m in a borehole intersecting an ultramafic rock formation undergoing serpentinization in the Samail Ophiolite, Sultanate of Oman. Near surface aquifer waters were moderately reducing and had alkaline pH while deeper aquifer waters were reduced and had hyperalkaline pH, indicating extensive influence by serpentinization. Metagenomic sequencing and analysis of DNA from filtered biomass collected from discrete depth intervals revealed an abundance of aerobes in near surface waters and a greater proportion of anaerobes at depth. Yet the abundance of the putatively obligate aerobe, Meiothermus, increased with depth, providing an opportunity to evaluate the influence of chemical and spatial variation on its distribution and speciation. Two clades of Meiothermus metagenome assembled genomes (MAGs) were identified that correspond to surface and deep populations termed Types I (S) and II (D), respectively; both clades comprised an apparently Oman-specific lineage indicating a common ancestor. Type II (D) clade MAGs encoded fewer genes and were undergoing slower genome replication as inferred from read mapping. Further, single nucleotide variants (SNVs) and mobile genetic elements identified among MAGs revealed detectable, albeit limited, evidence for gene flow/recombination between spatially segregated Type I (S) and Type II (D) populations. Together, these observations indicate that chemical variation generated by serpentinization, combined with physical barriers that reduce/limit dispersal and gene flow, allowed for the parapatric speciation of Meiothermus in the Samail Ophiolite or a geologic precursor. Further, Meiothermus genomic data suggest that deep and shallow aquifer fluids in the Samail Ophiolite may mix over shorter time scales than has been previously estimated from geochemical data

Data_Sheet_2_Parapatric speciation of Meiothermus in serpentinite-hosted aquifers in Oman.xlsx

The factors that control the distribution and evolution of microbial life in subsurface environments remain enigmatic due to challenges associated with sampling fluids from discrete depth intervals via boreholes while avoiding mixing of fluids. Here, using an inflatable packer system, fracture waters were isolated and collected from three discrete depth intervals spanning >130 m in a borehole intersecting an ultramafic rock formation undergoing serpentinization in the Samail Ophiolite, Sultanate of Oman. Near surface aquifer waters were moderately reducing and had alkaline pH while deeper aquifer waters were reduced and had hyperalkaline pH, indicating extensive influence by serpentinization. Metagenomic sequencing and analysis of DNA from filtered biomass collected from discrete depth intervals revealed an abundance of aerobes in near surface waters and a greater proportion of anaerobes at depth. Yet the abundance of the putatively obligate aerobe, Meiothermus, increased with depth, providing an opportunity to evaluate the influence of chemical and spatial variation on its distribution and speciation. Two clades of Meiothermus metagenome assembled genomes (MAGs) were identified that correspond to surface and deep populations termed Types I (S) and II (D), respectively; both clades comprised an apparently Oman-specific lineage indicating a common ancestor. Type II (D) clade MAGs encoded fewer genes and were undergoing slower genome replication as inferred from read mapping. Further, single nucleotide variants (SNVs) and mobile genetic elements identified among MAGs revealed detectable, albeit limited, evidence for gene flow/recombination between spatially segregated Type I (S) and Type II (D) populations. Together, these observations indicate that chemical variation generated by serpentinization, combined with physical barriers that reduce/limit dispersal and gene flow, allowed for the parapatric speciation of Meiothermus in the Samail Ophiolite or a geologic precursor. Further, Meiothermus genomic data suggest that deep and shallow aquifer fluids in the Samail Ophiolite may mix over shorter time scales than has been previously estimated from geochemical data.</p

Data_Sheet_1_Parapatric speciation of Meiothermus in serpentinite-hosted aquifers in Oman.DOCX

The factors that control the distribution and evolution of microbial life in subsurface environments remain enigmatic due to challenges associated with sampling fluids from discrete depth intervals via boreholes while avoiding mixing of fluids. Here, using an inflatable packer system, fracture waters were isolated and collected from three discrete depth intervals spanning >130 m in a borehole intersecting an ultramafic rock formation undergoing serpentinization in the Samail Ophiolite, Sultanate of Oman. Near surface aquifer waters were moderately reducing and had alkaline pH while deeper aquifer waters were reduced and had hyperalkaline pH, indicating extensive influence by serpentinization. Metagenomic sequencing and analysis of DNA from filtered biomass collected from discrete depth intervals revealed an abundance of aerobes in near surface waters and a greater proportion of anaerobes at depth. Yet the abundance of the putatively obligate aerobe, Meiothermus, increased with depth, providing an opportunity to evaluate the influence of chemical and spatial variation on its distribution and speciation. Two clades of Meiothermus metagenome assembled genomes (MAGs) were identified that correspond to surface and deep populations termed Types I (S) and II (D), respectively; both clades comprised an apparently Oman-specific lineage indicating a common ancestor. Type II (D) clade MAGs encoded fewer genes and were undergoing slower genome replication as inferred from read mapping. Further, single nucleotide variants (SNVs) and mobile genetic elements identified among MAGs revealed detectable, albeit limited, evidence for gene flow/recombination between spatially segregated Type I (S) and Type II (D) populations. Together, these observations indicate that chemical variation generated by serpentinization, combined with physical barriers that reduce/limit dispersal and gene flow, allowed for the parapatric speciation of Meiothermus in the Samail Ophiolite or a geologic precursor. Further, Meiothermus genomic data suggest that deep and shallow aquifer fluids in the Samail Ophiolite may mix over shorter time scales than has been previously estimated from geochemical data.</p

Data_Sheet_4_Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes.FASTA

<p>Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.</p

Data_Sheet_1_Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes.FASTA

<p>Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.</p

Table_3_Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes.XLSX

<p>Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.</p

Data_Sheet_7_Origin and Evolution of Flavin-Based Electron Bifurcating Enzymes.FASTA

<p>Twelve evolutionarily unrelated oxidoreductases form enzyme complexes that catalyze the simultaneous coupling of exergonic and endergonic oxidation–reduction reactions to circumvent thermodynamic barriers and minimize free energy loss in a process known as flavin-based electron bifurcation. Common to these 12 bifurcating (Bf) enzymes are protein-bound flavin, the proposed site of bifurcation, and the electron carrier ferredoxin. Despite the documented role of Bf enzymes in balancing the redox state of intracellular electron carriers and in improving the efficiency of cellular metabolism, a comprehensive description of the diversity and evolutionary history of Bf enzymes is lacking. Here, we report the taxonomic distribution, functional diversity, and evolutionary history of Bf enzyme homologs in 4,588 archaeal, bacterial, and eukaryal genomes and 3,136 community metagenomes. Bf homologs were primarily detected in the genomes of anaerobes, including those of sulfate-reducers, acetogens, fermenters, and methanogens. Phylogenetic analyses of Bf enzyme catalytic subunits (oxidoreductases) suggest they were not a property of the Last Universal Common Ancestor of Archaea and Bacteria, which is consistent with the limited and unique taxonomic distributions of enzyme homologs among genomes. Further, phylogenetic analyses of oxidoreductase subunits reveal that non-Bf homologs predate Bf homologs. These observations indicate that multiple independent recruitments of flavoproteins to existing oxidoreductases enabled coupling of numerous new electron Bf reactions. Consistent with the role of these enzymes in the energy metabolism of anaerobes, homologs of Bf enzymes were enriched in metagenomes from subsurface environments relative to those from surface environments. Phylogenetic analyses of homologs from metagenomes reveal that the earliest evolving homologs of most Bf enzymes are from subsurface environments, including fluids from subsurface rock fractures and hydrothermal systems. Collectively, these data suggest strong selective pressures drove the emergence of Bf enzyme complexes via recruitment of flavoproteins that allowed for an increase in the efficiency of cellular metabolism and improvement in energy capture in anaerobes inhabiting a variety of subsurface anoxic habitats where the energy yield of oxidation-reduction reactions is generally low.</p