A hot origin of dissimilatory sulfite reduction catalyzed by DsrAB in the Paleoarchean Era

Abstract

Dissimilatory sulfite reduction (DSR) has been essential to microbial energy metabolism in the biogeochemical sulfur cycle since the Paleoarchean Era. However, due to the lack of an integrated assessment of geological record and genomic data, the evolutionary origin of DSR remains elusive in terms of time, habitat, and genetic basis. In this study, we reconstructed the evolutionary pathways and the ancestral sequences of Dsr proteins by mining metagenomes ranging from mesothermal to hyperthermal environments. A phylogenetic analysis of the key catalytic enzyme, DsrAB, and other Dsr proteins indicates that the earliest and most basic functional cascade, DsrABCNM, emerged prior to the latest common ancestor of several basal branching DsrAB clusters encoded by bacteria and archaea. Using a molecular dating strategy that calibrates the protein tree with a species tree, we predicted that the DSR originated 3.508 billion years ago (Ga). This finding strongly confirms the earliest geological evidence of DSR ( ~ 3.47 Ga). Further predictions from ancestral sequence reconstruction indicate that the optimal catalytic temperature of DsrA at the time of DSR origin was approximately 73°C, which is consistent with the petrographic and geochemical evidence in early Archean hydrothermal deposits. After its hot origin, DsrA diversified into subclades that adapted to various temperature levels following the Great Oxidation Event. This is exemplified by the evolution of the reductive archaeal-type DsrA. Our results synchronize the molecular ages with the geological record, which advances our understanding of the earliest DSR systems and highlights the enzymatic adaptations of microbial life in the Archean biosphere.

Keywords

archaeal evolution / dsrAB / metagenomics / sulfite reduction / sulfur cycling

References

[1] Blank CE. Evolutionary timing of the origins of mesophilic sulphate reduction and oxygenic photosynthesis: a phylogenomic dating approach. Geobiology. 2004; 2: 1–20.

[2] Havig JR, Hamilton TL, Bachan A, Kump LR. Sulfur and carbon isotopic evidence for metabolic pathway evolution and a four-stepped earth system progression across the archean and paleoproterozoic. Earth-Sci Rev. 2017; 174: 1–21.

[3] Ranjan S, Todd ZR, Sutherland JD, Sasselov DD. Sulfidic anion concentrations on early Earth for surficial origins-of-life chemistry. Astrobiology. 2018; 18: 1023–1040.

[4] Halevy I. Production, preservation, and biological processing of mass-independent sulfur isotope fractionation in the Archean surface environment. Proc Natl Acad Sci USA. 2013; 110: 17644–17649.

[5] Zhou Z, Tran PQ, Cowley ES, Trembath-Reichert E, Anantharaman K. Diversity and ecology of microbial sulfur metabolism. Nat Rev Microbiol. 2025; 23: 122–140.

[6] Mateos K, Chappell G, Klos A, Le B, Boden J, Stüeken E, et al. The evolution and spread of sulfur cycling enzymes reflect the redox state of the early Earth. Sci Adv. 2023; 9:eade4847.

[7] Lyons TW, Tino CJ, Fournier GP, Anderson RE, Leavitt WD, Konhauser KO, et al. Co-evolution of early Earth environments and microbial life. Nat Rev Microbiol. 2024; 22: 572–586.

[8] Shen Y, Buick R, Canfield DE. Isotopic evidence for microbial sulphate reduction in the early Archaean era. Nature. 2001; 410: 77–81.

[9] Shen Y, Farquhar J, Masterson A, Kaufman AJ, Buick R. Evaluating the role of microbial sulfate reduction in the early Archean using quadruple isotope systematics. Earth Planet Sci Lett. 2009; 279: 383–391.

[10] Djokic T, Van Kranendonk MJ, Campbell KA, Walter MR, Ward CR. Earliest signs of life on land preserved in ca. 3.5 Ga hot spring deposits. Nat Commun. 2017; 8:15263.

[11] Loy A, Duller S, Baranyi C, Mußmann M, Ott J, Sharon I, et al. Reverse dissimilatory sulfite reductase as phylogenetic marker for a subgroup of sulfur-oxidizing prokaryotes. Environ Microbiol. 2009; 11: 289–299.

[12] Santos AA, Venceslau SS, Grein F, Leavitt WD, Dahl C, Johnston DT, et al. A protein trisulfide couples dissimilatory sulfate reduction to energy conservation. Science. 2015; 350: 1541–1545.

[13] Wagner M, Roger AJ, Flax JL, Brusseau GA, Stahl DA. Phylogeny of dissimilatory sulfite reductases supports an early origin of sulfate respiration. J Bacteriol. 1998; 180: 2975–2982.

[14] Zverlov V, Klein M, Lücker S, Friedrich MW, Kellermann J, Stahl DA, et al. Lateral gene transfer of dissimilatory (bi)sulfite reductase revisited. J Bacteriol. 2005; 187: 2203–2208.

[15] Anantharaman K, Hausmann B, Jungbluth SP, Kantor RS, Lavy A, Warren LA, et al. Expanded diversity of microbial groups that shape the dissimilatory sulfur cycle. ISME J. 2018; 12: 1715–1728.

[16] Müller AL, Kjeldsen KU, Rattei T, Pester M, Loy A. Phylogenetic and environmental diversity of DsrAB-type dissimilatory (bi)sulfite reductases. ISME J. 2015; 9: 1152–1165.

[17] Colman DR, Lindsay MR, Amenabar MJ, Fernandes-Martins MC, Roden ER, Boyd ES. Phylogenomic analysis of novel Diaforarchaea is consistent with sulfite but not sulfate reduction in volcanic environments on early Earth. ISME J. 2020; 14: 1316–1331.

[18] Diao M, Dyksma S, Koeksoy E, Ngugi DK, Anantharaman K, Loy A, et al. Global diversity and inferred ecophysiology of microorganisms with the potential for dissimilatory sulfate/sulfite reduction. FEMS Microbiol Rev. 2023; 47: 1–18.

[19] Rosenbaum FP, Poehlein A, Daniel R, Müller V. Energy-conserving dimethyl sulfoxide reduction in the acetogenic bacterium Moorella thermoacetica. Environ Microbiol. 2022; 24: 2000–2012.

[20] Ayala-Muñoz D, Burgos WD, Sánchez-España J, Couradeau E, Falagán C, Macalady JL. Metagenomic and metatranscriptomic study of microbial metal resistance in an acidic pit lake. Microorganisms. 2020; 8: 1350.

[21] Carr SA, Jungbluth SP, Eloe-Fadrosh EA, Stepanauskas R, Woyke T, Rappé MS, et al. Carboxydotrophy potential of uncultivated Hydrothermarchaeota from the subseafloor crustal biosphere. ISME J. 2019; 13: 1457–1468.

[22] Hua ZS, Qu YN, Zhu Q, Zhou EM, Qi YL, Yin YR, et al. Genomic inference of the metabolism and evolution of the archaeal phylum Aigarchaeota. Nat Commun. 2018; 9: 2832.

[23] Lin X, Handley KM, Gilbert JA, Kostka JE. Metabolic potential of fatty acid oxidation and anaerobic respiration by abundant members of Thaumarchaeota and Thermoplasmata in deep anoxic peat. ISME J. 2015; 9: 2740–2744.

[24] McKay LJ, Dlakić M, Fields MW, Delmont TO, Eren AM, Jay ZJ, et al. Co-occurring genomic capacity for anaerobic methane and dissimilatory sulfur metabolisms discovered in the Korarchaeota. Nat Microbiol. 2019; 4: 614–622.

[25] Luo ZH, Li Q, Xie YG, Lv AP, Qi YL, Li MM, et al. Temperature, pH, and oxygen availability contributed to the functional differentiation of ancient Nitrososphaeria. ISME J. 2024; 18:wrad031.

[26] Wu Z, Yang X, Huang L, Li S, Xia F, Qiu Y, et al. In situ enrichment of sulphate-reducing microbial communities with different carbon sources stimulating the acid mine drainage sediments. Sci Total Environ. 2023; 898:165584.

[27] Colman DR, Labesse G, Swapna GVT, Stefanakis J, Montelione GT, Boyd ES, et al. Structural evolution of the ancient enzyme, dissimilatory sulfite reductase. Proteins. 2022; 90: 1331–1345.

[28] Neukirchen S, Pereira IAC, Sousa FL. Stepwise pathway for early evolutionary assembly of dissimilatory sulfite and sulfate reduction. ISME J. 2023; 17: 1680–1692.

[29] Venceslau SS, Stockdreher Y, Dahl C, Pereira IAC. The “bacterial heterodisulfide” DsrC is a key protein in dissimilatory sulfur metabolism. Biochim Biophys Acta. 2014; 1837: 1148–1164.

[30] Grein F, Venceslau SS, Schneider L, Hildebrandt P, Todorovic S, Pereira IAC, et al. DsrJ, an essential part of the DsrMKJOP transmembrane complex in the purple sulfur bacterium Allochromatium vinosum, is an unusual triheme cytochrome c. Biochemistry. 2010; 49: 8290–8299.

[31] Rabus R, Venceslau SS, Wöhlbrand L, Voordouw G, Wall JD, Pereira IAC. A post-genomic view of the ecophysiology, catabolism and biotechnological relevance of sulphate-reducing prokaryotes. In: Advances in Microbial Physiology. Elsevier; 2015. p. 55–321.

[32] Ferreira D, Barbosa ACC, Oliveira GP, Catarino T, Venceslau SS, Pereira IAC. The DsrD functional marker protein is an allosteric activator of the DsrAB dissimilatory sulfite reductase. Proc Natl Acad Sci USA. 2022; 119:e2118880119.

[33] Akbar S, Gaidenko TA, Kang CM, O'Reilly M, Devine KM, Price CW. New family of regulators in the environmental signaling pathway which activates the general stress transcription factor ς ^B of Bacillus subtilis. J Bacteriol. 2001; 183: 1329–1338.

[34] Stockdreher Y, Venceslau SS, Josten M, Sahl H-G, Pereira IAC, Dahl C. Cytoplasmic sulfurtransferases in the purple sulfur bacterium Allochromatium vinosum: evidence for sulfur transfer from DsrEFH to DsrC. PLoS One. 2012; 7:e40785.

[35] Thorup C, Schramm A, Findlay AJ, Finster KW, Schreiber L. Disguised as a sulfate reducer: growth of the Deltaproteobacterium Desulfurivibrio alkaliphilus by sulfide oxidation with nitrate. mBio. 2017; 8:e00671-17.

[36] Parks DH, Chuvochina M, Rinke C, Mussig AJ, Chaumeil P-A, Hugenholtz P. GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy. Nucleic Acids Res. 2022; 50: D785–D794.

[37] Gao S, Paez-Espino D, Li J, Ai H, Liang J, Luo Z, et al. Patterns and ecological drivers of viral communities in acid mine drainage sediments across Southern China. Nat Commun. 2022; 13: 2389.

[38] Luo ZH, Li Q, Chen N, Tang LY, Liao B, Yang TT, et al. Genome-resolved metagenomics reveals depth-related patterns of microbial community structure and functions in a highly stratified, AMD overlaying mine tailings. J Hazard Mater. 2023; 447:130774.

[39] Huber R, Kristjansson JK, Stetter KO. Pyrobaculum gen. nov., a new genus of neutrophilic, rod-shaped archaebacteria from continental solfataras growing optimally at 100°C. Arch Microbiol. 1987; 149: 95–101.

[40] Larsen Ø, Lien T, Birkeland N. Dissimilatory sulfite reductase from Archaeoglobus profundus and Desulfotomaculum thermocisternum: phylogenetic and structural implications from gene sequences. Extremophiles. 1999; 3: 63–70.

[41] Pierce E, Xie G, Barabote RD, Saunders E, Han CS, Detter JC, et al. The complete genome sequence of Moorella thermoacetica (f. Clostridium thermoaceticum). Environ Microbiol. 2008; 10: 2550–2573.

[42] Cawood PA, Strachan RA, Pisarevsky SA, Gladkochub DP, Murphy JB. Linking collisional and accretionary orogens during Rodinia assembly and breakup: implications for models of supercontinent cycles. Earth Planet Sci Lett. 2016; 449: 118–126.

[43] Ueno Y, Ono S, Rumble D, Maruyama S. Quadruple sulfur isotope analysis of ca. 3.5 Ga dresser formation: new evidence for microbial sulfate reduction in the early Archean. Geochim Cosmochim Acta. 2008; 72: 5675–5691.

[44] Roerdink DL, Mason PRD, Farquhar J, Reimer T. Multiple sulfur isotopes in Paleoarchean barites identify an important role for microbial sulfate reduction in the early marine environment. Earth Planet Sci Lett. 2012; 331–332: 177–186.

[45] Wacey D, Kilburn MR, Saunders M, Cliff J, Brasier MD. Microfossils of sulphur-metabolizing cells in 3.4-billion-year-old rocks of Western Australia. Nat Geosci. 2011; 4: 698–702.

[46] Vysotskiy S, Velivetskaya T, Ignatiev A, Slabunov A, Aseeva A. Multiple sulfur isotope evidence for bacterial sulfate reduction and sulfate disproportionation operated in Mesoarchaean rocks of the Karelian Craton. Minerals. 2022; 12: 1143.

[47] Zhelezinskaia I, Kaufman AJ, Farquhar J, Cliff J. Large sulfur isotope fractionations associated with Neoarchean microbial sulfate reduction. Science. 2014; 346: 742–744.

[48] Kakegawa T, Kawai H, Ohmoto H. Origins of pyrites in the ∼2.5 Ga Mt. McRae shale, the Hamersley District, Western Australia. Geochim Cosmochim Acta. 1998; 62: 3205–3220.

[49] Uveges BT, Izon G, Junium CK, Ono S, Summons RE. Aerobic nitrogen cycle 100 My before permanent atmospheric oxygenation. Proc Natl Acad Sci USA. 2025; 122:e2423481122.

[50] Mayika KB, Moussavou M, Prave AR, Lepland A, Mbina M, Kirsimäe K. The Paleoproterozoic Francevillian succession of Gabon and the Lomagundi-Jatuli event. Geology. 2020; 48: 1099–1104.

[51] Anbar AD, Buick R, Gordon GW, Johnson AC, Kendall B, Lyons TW, et al. Technical comment on “Reexamination of 2.5-Ga ‘whiff’ of oxygen interval points to anoxic ocean before GOE. Sci Adv. 2023; 9:eabq3736.

[52] Li ZX, Bogdanova SV, Collins AS, Davidson A, De Waele B, Ernst RE, et al. Assembly, configuration, and break-up history of Rodinia: a synthesis. Precambrian Res. 2008; 160: 179–210.

[53] Hickey DA, Singer GA. Genomic and proteomic adaptations to growth at high temperature. Genome Biol. 2004; 5: 117.

[54] Zeldovich KB, Berezovsky IN, Shakhnovich EI. Protein and DNA sequence determinants of thermophilic adaptation. PLoS Comput Biol. 2007; 3:e5.

[55] Li G, Buric F, Zrimec J, Viknander S, Nielsen J, Zelezniak A, et al. Learning deep representations of enzyme thermal adaptation. Prot Sci. 2022; 31:e4480.

[56] Planavsky NJ, Asael D, Hofmann A, Reinhard CT, Lalonde SV, Knudsen A, et al. Evidence for oxygenic photosynthesis half a billion years before the Great Oxidation Event. Nat Geosci. 2014; 7: 283–286.

[57] Donald E, Canfield R. Raiswell. The evolution of the sulfur cycle. Am J Sci. 1999; 299: 697.

[58] Fakhraee M, Hancisse O, Canfield DE, Crowe SA, Katsev S. Proterozoic seawater sulfate scarcity and the evolution of ocean–atmosphere chemistry. Nat Geosci. 2019; 12: 375–380.

[59] Dickin AP, McNutt RH, Martin C, Guo A. The extent of juvenile crust in the Grenville Province: Nd isotope evidence. Geol Soc Am Bull. 2010; 122: 870–883.

[60] Liu C, Runyon SE, Knoll AH, Hazen RM. The same and not the same: ore geology, mineralogy and geochemistry of Rodinia assembly versus other supercontinents. Earth-Sci Rev. 2019; 196:102860.

[61] Konhauser KO, Lalonde SV, Planavsky NJ, Pecoits E, Lyons TW, Mojzsis SJ, et al. Aerobic bacterial pyrite oxidation and acid rock drainage during the Great Oxidation Event. Nature. 2011; 478: 369–373.

[62] Stüeken EE, Catling DC, Buick R. Contributions to late Archaean sulphur cycling by life on land. Nat Geosci. 2012; 5: 722–725.

[63] Havig JR, Grettenberger C, Hamilton TL. Geochemistry and microbial community composition across a range of acid mine drainage impact and implications for the Neoarchean-Paleoproterozoic transition. J Geophys Res Biogeosci. 2017; 122: 1404–1422.

[64] Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010; 11: 119.

[65] Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2020; 36: 1925–1927.

[66] Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020; 36: 2251–2252.

[67] Li W, O'Neill KR, Haft DH, DiCuccio M, Chetvernin V, Badretdin A, et al. RefSeq: expanding the Prokaryotic Genome Annotation Pipeline reach with protein family model curation. Nucleic Acids Res. 2020; 49: 1020–1028.

[68] Anantharaman K, Brown CT, Hug LA, Sharon I, Castelle CJ, Probst AJ, et al. Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat Commun. 2016; 7:13219.

[69] Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011; 7:e1002195.

[70] Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012; 28: 3150–3152.

[71] Neukirchen S, Sousa FL. DiSCo: a sequence-based type-specific predictor of Dsr-dependent dissimilatory sulphur metabolism in microbial data. Microb Genom. 2021; 7:000603.

[72] Yu X, Zhou J, Song W, Xu M, He Q, Peng Y, et al. SCycDB: a curated functional gene database for metagenomic profiling of sulphur cycling pathways. Mol Ecol Resour. 2021; 21: 924–940.

[73] Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013; 30: 772–780.

[74] Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009; 25: 1972–1973.

[75] Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020; 37: 1530–1534.

[76] Letunic I, Bork P. Interactive tree of life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 2024; 52: W78–W82.

[77] Menardo F, Loiseau C, Brites D, Coscolla M, Gygli SM, Rutaihwa LK, et al. Treemmer: a tool to reduce large phylogenetic datasets with minimal loss of diversity. BMC Bioinformatics. 2018; 19: 164.

[78] Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24: 1586–1591.

[79] Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015; 25: 1043–1055.

[80] Woodcroft BJ, Singleton CM, Boyd JA, Evans PN, Emerson JB, Zayed AAF, et al. Genome-centric view of carbon processing in thawing permafrost. Nature. 2018; 560: 49–54.

[81] Altenhoff AM, Warwick Vesztrocy A, Bernard C, Train C-M, Nicheperovich A, Prieto Baños S, et al. OMA orthology in 2024: improved prokaryote coverage, ancestral and extant GO enrichment, a revamped synteny viewer and more in the OMA ecosystem. Nucleic Acids Res. 2024; 52: D513–D521.

[82] Schön ME, Martijn J, Vosseberg J, Köstlbacher S, Ettema TJG. The evolutionary origin of host association in the Rickettsiales. Nat Microbiol. 2022; 7: 1189–1199.

[83] Martijn J, Vosseberg J, Guy L, Offre P, Ettema TJG. Deep mitochondrial origin outside the sampled Alphaproteobacteria. Nature. 2018; 557: 101–105.

[84] Tamura K, Tao Q, Kumar S. Theoretical foundation of the RelTime method for estimating divergence times from variable evolutionary rates. Mol Biol Evol. 2018; 35: 1770–1782.

[85] Valley JW, Cavosie AJ, Ushikubo T, Reinhard DA, Lawrence DF, Larson DJ, et al. Hadean age for a post-magma-ocean zircon confirmed by atom-probe tomography. Nat Geosci. 2014; 7: 219–223.

[86] Blank CE. An expansion of age constraints for microbial clades that lack a conventional fossil record using phylogenomic dating. J Mol Evol. 2011; 73: 188–208.

[87] Blank CE. Phylogenomic dating—a method of constraining the age of microbial taxa that lack a conventional fossil record. Astrobiology. 2009; 9: 173–191.

[88] Xie J, Chen Y, Cai G, Cai R, Hu Z, Wang H. Tree visualization by one table (tvBOT): a web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res. 2023; 51: W587–W592.

[89] Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015; 12: 59–60.

[90] Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012; 19: 455–477.

[91] Aroney STN, Newell RJP, Nissen J, Camargo AP, Tyson GW, Woodcroft BJ. CoverM: read coverage calculator for metagenomics. https://doi.org/10.5281/zenodo.10531254. (January 2024, date last accessed).

[92] Wei T, Simko V. R package “corrplot”: visualization of a correlation matrix. https://github.com/taiyun/corrplot. (2017, date last accessed).

[93] Dixon P. VEGAN, a package of R functions for community ecology. J Vegetation Sci. 2003; 14: 927–930.

[94] Aphalo PJ. ggpmisc: miscellaneous extensions to ‘ggplot2’. https://docs.r4photobiology.info/ggpmisc/. (2024, date last accessed).

[95] Wickham H. ggplot2: elegant graphics for data analysis. 2nd ed. Springer International Publishing; 2016.

[96] Yu G. gglayer: provide extra layers for ggplot2. https://github.com/GuangchuangYu/gglayer. (2019, date last accessed).

RIGHTS & PERMISSIONS

2026 The Author(s). mLife published by John Wiley & Sons Australia, Ltd on behalf of Institute of Microbiology, Chinese Academy of Sciences.