- 1Biodesign Center for Mechanisms of Evolution, School of Life Sciences, Arizona State University, Tempe, AZ, United States
- 2Department of Earth Sciences, University of Southern California, Los Angeles, CA, United States
- 3Biosciences & Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA, United States
Multi-heme cytochromes (MHCs), together with accessory proteins like porins and periplasmic cytochromes, enable microbes to transport electrons between the cytoplasmic membrane and extracellular substrates (e.g., minerals, electrodes, other cells). Extracellular electron transfer (EET) has been described in multiple systems; yet, the broad phylogenetic and mechanistic diversity of these pathways is less clear. One commonality in EET-capable systems is the involvement of MHCs, in the form of porin-cytochrome complexes, pili-like cytochrome polymers, and lipid-anchored extracellular cytochromes. Here, we put forth MHCscan—a software tool for identifying MHCs and identifying potential EET capability. Using MHCscan, we scanned ~60,000 bacterial and 2,000 archaeal assemblies, and identify a diversity of MHCs, many of which represent enzymes with no known function, and many found within organisms not previously known to be electroactive. In total, our scan identified ~1,400 unique enzymes, each encoding more than 10 heme-binding motifs. In our analysis, we also find evidence for modularity and flexibility in MHC-dependent EET pathways, and suggest that MHCs may be far more common than previously recognized, with many facets yet to be discovered. We present MHCscan as a lightweight and user-friendly software tool that is freely available: https://github.com/Arkadiy-Garber/MHCscan.
Introduction
The ability of cells to either donate or accept electrons from external insoluble substrates has been well-studied in several microbial groups [e.g., Shewanella oneidensis (Myers and Myers, 1992, 2002; Pitts et al., 2003; Hartshorne et al., 2009; Clarke et al., 2011; Kasai et al., 2015), Geobacter sulfurreducens (Shi et al., 2009; Merkley et al., 2015) Rhodopseudomonas palustris (Jiao and Newman, 2007; Bird et al., 2014), and Ferroglobus placidus (Smith et al., 2015)]. However, the importance of many of these microbe-mineral interactions remain hidden in the to-be-discovered intricacies of energy exchange and syntrophic interactions of the microbial world. Many of the known electrochemical interactions have recently been outlined and reviewed by Lovley and Holmes (2022); there are likely many more waiting to be discovered, as discussed below. Known mechanisms for EET include only a handful of protein families, many of which consist of multi-heme cytochromes (MHCs) that are localized to the periplasm and outer membrane of the cell (Shi et al., 2009; Richardson et al., 2012; Shi et al., 2014; Liu et al., 2014; He et al., 2017; Deng et al., 2018; Barco et al., 2015) (Figure 1). This is particularly true in the case of Gram-negative bacteria, and at least one Gram-positive species (Thermincola potens) is known to encode MHCs that link the cytoplasmic membrane with extracellular substrates through the cell wall (Carlson et al., 2012).
Figure 1. Schematic of electron transport across the periplasmic space and outer membrane. Outer membrane-bound multi-heme cytochromes are labeled with “C”. The multimeric OmcS cytochromes (shown in brown at the right of the figure) represent multi-heme cytochromes that are bound to the inner membrane, but extend out past the outer membrane. Accessory proteins are labeled with “A”; these represent porins and soluble periplasmic cytochromes that transport electrons from the inner membrane-bound electron carriers (labeled with “R”) to the outer membrane cytochromes. Arrows represent the flow of electrons.
MHCs that are known to play a role in EET have been used as queries to search for homologous proteins and mechanisms in other microorganisms. For example, profile hidden Markov models (HMMs) have been compiled for known MHCs with EET function and are currently packaged in a bioinformatics tool called FeGenie to predict genetic potential for iron-cycling (Garber et al., 2020). However, many types of organisms are being discovered that are known to generate electrical currents (Logan et al., 2019), many of which do not appear to encode any of the related EET machinery. It is, thus, possible that these microbes produce MHCs with independently evolved EET pathways invisible to homology-dependent searches. With this in mind, we hypothesized that MHCs can be used to predict EET capability without the use of homology-dependent methods.
We searched the predicted proteomes of ~60,000 sequenced Bacteria and Archaea (from NCBI's RefSeq database) for gene neighborhoods encoding one or more MHC localized to the periplasm, outer-membrane, and/or extracellular space. This homology-independent approach (Supplementary Figure 1), revealed thousands of MHC-encoding gene neighborhoods, some predicted to have over a hundred heme-binding motifs, potentially involved in EET. Using FeGenie as a reference, we were able to identify novel MHCs for which gene markers are not yet publicly available. While FeGenie includes 14 protein families representing MHCs, sequence-based clustering of novel MHCs identified in our survey resulted in an additional 76 protein families encoding heme-rich (≥4 hemes) cytochromes that are localized to the periplasm, outer membrane, or extracellular space. Additionally, we found 42 novel protein families representing accessory genes, likely encompassing porins and soluble electron shuttles.
Results and discussion
Validation of the homology-independent approach
To validate MHCscan, we used it to analyze a set of genomes from model organisms with characterized EET genes and pathways [e.g., Rhodopseudomonas palustris (Jiao and Newman, 2007), members of the Shewanella and Geobacter genera (Shi et al., 2012), and others]. All genomes analyzed in this benchmark are listed in Supplementary File 1, a select subset are listed in Table 1. Our criteria for the inclusion of gene neighborhoods potentially involved in EET were: (1) the presence of at least one MHC with a signal peptide; (2) the presence of at least one heme-rich (≥4 hemes) MHC; and, (3) the combined presence of at least 15 hemes within the entire multi-gene cluster. The exact parameters of the search can be tuned to include, for example, the single-heme porin-cytochrome Cyc2 or other outer-membrane cytochromes not matching the above criteria. Moreover, MHCscan will also report type IV aromatic pili (T4AP), which generally do not encode heme-binding motifs, but may nonetheless be involved in EET (Vargas et al., 2013). Our search with MHCscan accurately identified all known EET-related mechanisms in the subset of electroactive and iron reducing/oxidizing microbes. Moreover, within even these well-studied bacteria, we found additional gene neighborhoods with as-of-yet unnamed genes encoding multiheme cytochromes (Supplementary Figure 2, Supplementary Files 1, 2). We note that MHCs encoded by as-of-yet unnamed genes may have been described in previous research studies—nonetheless, appear as “hypothetical” in public repositories. We also used MHCscan to analyze genomes from 70 electroactive microbes included in the review by Logan et al. (2019). We compared predicted MHCs from MHCscan to EET-related genes identified through the homology-dependent software tool FeGenie. In these electroactive microbes, MHCscan detects more EET candidate genes than FeGenie (Figure 2), possibly illuminating a greater proportion of heme groups involved in facilitating extracellular electrical currents.
Table 1. GenBank assembly accessions and publications associated with the genomes used in validation of the homology-independent approach.
Figure 2. MHCs identified in 70 electroactive microbes as reviewed in Logan et al. (2019). In nearly all microbes, MHCScan identified more EET-related proteins than FeGenie. These EET candidates may represent novel mechanisms for electroactivity. A phylogenetic tree was constructed using GToTree and annotated using the Interactive Tree of Life (Letunic and Bork, 2021). Whether each microbe is an electrogen or electrotroph is denoted by the purple circles or green triangles to the right of the tree. The heatmap shows, per genome, the number of hemes, as well as MHC-encoding gene clusters potentially involved in EET. The number of MHC-encoding gene clusters identified with MHCScan (middle column) is generally higher than the number of gene clusters identified through the homology-dependent FeGenie software, suggesting new EET candidate genes.
Survey of bacteria and Archaea
We next turned our approach to all sequenced Bacteria and Archaea, down-sampled to reduce the number of over-represented genomes (e.g., E. coli, Streptococcus and Salmonella spp.). We identified 1,475 bacterial genomes (220 unique genera) and 75 archaeal genomes (35 unique genera), each with at least one gene neighborhood that encodes one or more clusters of MHCs. Most of these genomes encode only one such gene neighborhood; many encode a few, and a few genomes encode more than ten. In total, we identified over two thousand gene neighborhoods with potential involvement in EET (listed in Supplementary Files 1, 2 for Bacteria and Archaea, respectively). Of these, 204 gene neighborhoods encoded MHCs that can be identified using FeGenie's HMM library; 440 gene clusters were completely novel with no homologs present in FeGenie's library of HMMs (Figure 3). Species within Vibrio, Shewanella, Geobacter, Geomonas, and Aeromonas (e.g., known electroactive microbes) make up the majority of the lineages that we identified with our approach. These well-studied genera contain both the greatest diversity and the total number of MHC-encoding gene neighborhoods.
Figure 3. Scatterplot showing the relationship between the number of individual MHCs, and total number of hemes predicted within each gene neighborhood. Each dot represents a single gene neighborhood: those colored blue represent gene neighborhoods where all MHCs were predicted with FeGenie; green dots represent neighborhoods encoding at least one MHC not included in FeGenie's HMM library; red dots represent gene neighborhoods where not a single MHC in the gene neighborhood is currently included in FeGenie's HMM library. The high density of neighborhoods encoding multiples of ten hemes (horizontal lines at 20, 30, 50, 60, and 70 hemes per neighborhood) represent the decaheme cytochromes (MtrA/MtrC) encoded mostly in the Shewanella spp.; those encoding 22 and 27 hemes per neighborhood represent Mtr-encoding gene clusters in Vibrio and Aeromonas spp., which also encode pentaheme and/or diheme cytochromes.
Novel MHC families
There are 14 EET-linked MHCs represented within FeGenie's library of HMMs. Adding the MHCs that were identified by MHCscan, we identified an additional 118 protein families, each representing either a heme-rich MHC or accessory protein found within the cluster. Within this group of 118 novel HMMs, 76 represent MHCs that are predicted to bind five or more hemes. Notably, the greatest number of identified MHC-encoding gene neighborhoods were related to Shewanella's MtrCAB iron reductase complex (crystal structure in Edwards et al., 2020), in addition to DFE_0462 and DFE_0449, which encode MHCs in Desulfovibrio ferrophilus (Deng et al., 2018). It's possible that these pathways may be most efficient or most-easily transferred between bacteria via lateral gene transfer (LGT), as previously suggested for MtrCAB (Baker et al., 2021). Nonetheless, we uncovered many species and genera with no previous evidence for EET (Supplementary Files 1, 2) and suggest that EET is far more common than previously realized.
Concluding remarks
Homology-dependent approaches that utilize sequence alignment (e.g., BLAST) and hidden Markov models (e.g., HMMER) allow for identification of evolutionarily conserved proteins. However, it's possible for certain structures and functions to arise independently. As an example, extracellular electron transport mechanisms are found across all reaches of the bacterial phyla; yet, only a handful of mechanisms are known. Our approach leverages the structural features and motifs of known EET-related enzymes to identify similar mechanisms in unrelated proteins, which may have evolved independently or diverged beyond where homology-dependent approaches can detect. The growing reportoire of software tools for de novo protein folding (e.g., AlphaFold3 and ESMFold) will also provide a powerful approach for the preparation of biochemical models for testing through in vitro and in vivo experimentation.
Methods
Genome acquisition
Over 100,000 genome assemblies were downloaded from NCBI's RefSeq database (release 209). We down-sampled the genome assemblies from well-studied, but over-represented, species (e.g., Escherichia coli, Salmonella enterica, and Yersinia pestis, etc.), and reduced the number of genomes to ~60,000.
Workflow
The MHCscan pipeline, outlined in Supplementary Figure 1, starts with prediction of genes and corresponding protein sequences from genomes using Prodigal (Hyatt et al., 2010). If annotation was done beforehand, annotated genomes can be provided in the form of FASTA amino acid and GFF files. Phobius (Käll et al., 2004) is used to predict signal peptides and transmembrane helices. Heme c-binding motifs are predicted using FindMeHemes (https://github.com/Arkadiy-Garber/FindMeHemes). Output from Phobius and FindMeHemes is then processed to identify gene neighborhoods encoding multi-heme cytochromes that have signal peptides and appear to be periplasmic, outer membrane-bound, or secreted. We then removed the MHC gene clusters where the total number of predicted heme groups (across all proteins encoded within each cluster) falls below 15, an empirically derived number based on MHCs in known electroactive microbes.
HMM development
To bring into light the novel MHCs potentially involved in EET, we used the identified protein sequences to generate hidden Markov models (HMMs). We performed clustering of identified MHC protein sequences using MMseqs2 (Steinegger and Söding, 2017). These clustered protein sequences were aligned using Muscle (Edgar, 2004), masked to remove positions consisting of mostly gaps, and used to generate hidden Markov models (HMMs) using HMMER (Johnson et al., 2010). HMMs were only generated for proteins families that were represented by at least three orthologs among the roughly 60,000 surveyed genomes. The only custom script used in HMM generated was masker.py, available here: https://github.com/Arkadiy-Garber/BagOfTricks.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
AG: Writing – review & editing, Writing – original draft, Visualization, Validation, Software, Data curation. KN: Writing – review & editing, Conceptualization. NM: Writing – original draft, Supervision, Funding acquisition, Data curation, Conceptualization.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. NM was supported by LLNL LDRD grant (22-LW-034), and work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1448685/full#supplementary-material
References
Anderson, I., Risso, C., Holmes, D., Lucas, S., Copeland, A., Lapidus, A., et al. (2011). Complete genome sequence of Ferroglobus placidus AEDII12DO. Stand. Genomic Sci. 5, 50–60. doi: 10.4056/sigs.2225018
Baker, I. R., Conley, B. E., Gralnick, J. A., and Girguis, P. R. (2021). Evidence for horizontal and vertical transmission of mtr-mediated extracellular electron transfer among the bacteria. MBio 13:e0290421. doi: 10.1128/mbio.02904-21
Barco, R. A., Emerson, D., Sylvan, J. B., Orcutt, B. N., Jacobson Meyers, M. E., Ramírez, G. A., et al. (2015). New insight into microbial iron oxidation as revealed by the proteomic profile of an obligate iron-oxidizing chemolithoautotroph. Appl. Environ. Microbiol. 81, 5927–5937. doi: 10.1128/AEM.01374-15
Bird, L. J., Saraiva, I. H., Park, S., Calçada, E. O., Salgueiro, C. A., Nitschke, W., et al. (2014). Nonredundant roles for cytochrome c2 and two high-potential iron-sulfur proteins in the photoferrotroph Rhodopseudomonas palustris TIE-1. J. Bacteriol. 196, 850–858. doi: 10.1128/JB.00843-13
Cai, C., Leu, A. O., Xie, G.-J., Guo, J., Feng, Y., Zhao, J.-X., et al. (2018). A methanotrophic archaeon couples anaerobic oxidation of methane to Fe(III) reduction. ISME J. 12, 1929–1939. doi: 10.1038/s41396-018-0109-x
Carlson, H. K., Iavarone, A. T., Gorur, A., Yeo, B. S., Tran, R., Melnyk, R. A., et al. (2012). Surface multiheme c-type cytochromes from Thermincola potens and implications for respiratory metal reduction by Gram-positive bacteria. Proc. Natl. Acad. Sci. USA. 109, 1702–1707. doi: 10.1073/pnas.1112905109
Clarke, T. A., Edwards, M. J., Gates, A. J., Hall, A., White, G. F., Bradley, J., et al. (2011). Structure of a bacterial cell surface decaheme electron conduit. Proc. Natl. Acad. Sci. USA. 108, 9384–9389. doi: 10.1073/pnas.1017200108
Deng, X., Dohmae, N., Nealson, K. H., Hashimoto, K., and Okamoto, A. (2018). Multi-heme cytochromes provide a pathway for survival in energy-limited environments. Sci. Adv. 4:eaao5682. doi: 10.1126/sciadv.aao5682
Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340
Edwards, M. J., White, G. F., Butt, J. N., Richardson, D. J., and Clarke, T. A. (2020). The crystal structure of a biological insulated transmembrane molecular wire. Cell 181, 665–673.e10. doi: 10.1016/j.cell.2020.03.032
Emerson, D., Field, E. K., Chertkov, O., Davenport, K. W., Goodwin, L., Munk, C., et al. (2013). Comparative genomics of freshwater Fe-oxidizing bacteria: implications for physiology, ecology, and systematics. Front. Microbiol. 4:254. doi: 10.3389/fmicb.2013.00254
Garber, A. I., Nealson, K. H., Okamoto, A., McAllister, S. M., Chan, C. S., Barco, R. A., et al. (2020). FeGenie: a comprehensive tool for the identification of iron genes and iron gene neighborhoods in genome and metagenome assemblies. Front. Microbiol. 11:37. doi: 10.3389/fmicb.2020.00037
Hartshorne, R. S., Reardon, C. L., Ross, D., Nuester, J., Clarke, T. A., Gates, A. J., et al. (2009). Characterization of an electron conduit between bacteria and the extracellular environment. Proc. Natl. Acad. Sci. USA. 106, 22169–22174. doi: 10.1073/pnas.0900086106
He, S., Barco, R. A., Emerson, D., and Roden, E. E. (2017). Comparative genomic analysis of neutrophilic iron(II) oxidizer genomes for candidate genes in extracellular electron transfer. Front. Microbiol. 8:1584. doi: 10.3389/fmicb.2017.01584
Hyatt, D., Chen, G.-L., Locascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119
Jiao, Y., and Newman, D. K. (2007). The pio operon is essential for phototrophic Fe(II) oxidation in Rhodopseudomonas palustris TIE-1. J. Bacteriol. 189, 1765–1773. doi: 10.1128/JB.00776-06
Johnson, L. S., Eddy, S. R., and Portugaly, E. (2010). Hidden Markov model speed heuristic and iterative HMM search procedure. BMC Bioinformatics 11:431. doi: 10.1186/1471-2105-11-431
Käll, L., Krogh, A., and Sonnhammer, E. L. L. (2004). A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 338, 1027–1036. doi: 10.1016/j.jmb.2004.03.016
Kasai, T., Kouzuma, A., Nojiri, H., and Watanabe, K. (2015). Transcriptional mechanisms for differential expression of outer membrane cytochrome genes omcA and mtrC in Shewanella oneidensis MR-1. BMC Microbiol. 15:68. doi: 10.1186/s12866-015-0406-8
Kim, S.-J., Park, S.-J., Cha, I.-T., Min, D., Kim, J.-S., Chung, W.-H., et al. (2014). Metabolic versatility of toluene-degrading, iron-reducing bacteria in tidal flat sediment, characterized by stable isotope probing-based metagenomic analysis. Environ. Microbiol. 16, 189–204. doi: 10.1111/1462-2920.12277
Laso-P'erez, R., Wegener, G., Knittel, K., Widdel, F., Harding, K. J., Krukenberg, V., et al. (2016). Thermophilic archaea activate butane via alkyl-coenzyme M formation. Nature 539, 396–401. doi: 10.1038/nature20152
Letunic, I., and Bork, P. (2021). Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293-W296. doi: 10.1093/nar/gkab301
Liu, Y., Wang, Z., Liu, J., Levar, C., Edwards, M. J., Babauta, J. T., et al. (2014). A trans-outer membrane porin-cytochrome protein complex for extracellular electron transfer by Geobacter sulfurreducens PCA. Environ. Microbiol. Rep. 6, 776–785. doi: 10.1111/1758-2229.12204
Logan, B. E., Rossi, R., Ragab, A., and Saikaly, P. E. (2019). Electroactive microorganisms in bioelectrochemical systems. Nat. Rev. Microbiol.17, 307–319. doi: 10.1038/s41579-019-0173-x
Lovley, D. R., and Holmes, D. E. (2022). Electromicrobiology: the ecophysiology of phylogenetically diverse electroactive microorganisms. Nat. Rev. Microbiol. 20, 5–19. doi: 10.1038/s41579-021-00597-6
Merkley, E. D., Wrighton, K. C., Castelle, C. J., Anderson, B. J., Wilkins, M. J., Shah, V., et al. (2015). Changes in protein expression across laboratory and field experiments in Geobacter bemidjiensis. J. Proteome Res. 14, 1361–1375. doi: 10.1021/pr500983v
Myers, C. R., and Myers, J. M. (1992). Localization of cytochromes to the outer membrane of anaerobically grown Shewanella putrefaciens MR-1. J. Bacteriol. 174, 3429–3438. doi: 10.1128/jb.174.11.3429-3438.1992
Myers, C. R., and Myers, J. M. (2002). MtrB is required for proper incorporation of the cytochromes OmcA and OmcB into the outer membrane of Shewanella putrefaciens MR-1. Appl. Environ. Microbiol. 68, 5585–5594. doi: 10.1128/AEM.68.11.5585-5594.2002
Nolan, M., Sikorski, J., Davenport, K., Lucas, S., Del Rio, T. G., Tice, H., et al. (2010). Complete genome sequence of Ferrimonas balearica type strain (PAT). Stand. Genomic Sci. 3, 174–182. doi: 10.4056/sigs.1161239
Pitts, K. E., Dobbin, P. S., Reyes-Ramirez, F., Thomson, A. J., Richardson, D. J., Seward, H. E., et al. (2003). Characterization of the Shewanella oneidensis MR-1 decaheme cytochrome MtrA: expression in Escherichia coli confers the ability to reduce soluble Fe(III) chelates. J. Biol. Chem. 278, 27758–27765. doi: 10.1074/jbc.M302582200
Richardson, D. J., Butt, J. N., Fredrickson, J. K., Zachara, J. M., Shi, L., Edwards, M. J., et al. (2012). The “porin-cytochrome” model for microbe-to-mineral electron transfer. Mol. Microbiol. 85, 201–212. doi: 10.1111/j.1365-2958.2012.08088.x
Shi, L., Fredrickson, J. K., and Zachara, J. M. (2014). Genomic analyses of bacterial porin-cytochrome gene clusters. Front. Microbiol. 5:657. doi: 10.3389/fmicb.2014.00657
Shi, L., Richardson, D. J., Wang, Z., Kerisit, S. N., Rosso, K. M., Zachara, J. M., et al. (2009). The roles of outer membrane cytochromes of Shewanella and Geobacter in extracellular electron transfer. Environ. Microbiol. Rep. 1, 220–227. doi: 10.1111/j.1758-2229.2009.00035.x
Shi, L., Rosso, K. M., Zachara, J. M., and Fredrickson, J. K. (2012). Mtr extracellular electron-transfer pathways in Fe(III)-reducing or Fe(II)-oxidizing bacteria: a genomic perspective. Biochem. Soc. Trans. 40, 1261–1267. doi: 10.1042/BST20120098
Smith, J. A., Aklujkar, M., Risso, C., Leang, C., Giloteaux, L., Holmes, D. E., et al. (2015). Mechanisms involved in Fe(III) respiration by the hyperthermophilic archaeon Ferroglobus placidus. Appl. Environ. Microbiol. 81, 2735–2744. doi: 10.1128/AEM.04038-14
Steinegger, M., and Söding, J. (2017). MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat. Biotechnol. 35, 1026–1028. doi: 10.1038/nbt.3988
Sun, D., Wang, A., Cheng, S., Yates, M., and Logan, B. E. (2014). Geobacter anodireducens sp. nov., an exoelectrogenic microbe in bioelectrochemical systems. Int. J. Syst. Evol. Microbiol. 64(Pt 10), 3485–3491. doi: 10.1099/ijs.0.061598-0
Keywords: electroactive, extracellular electron transport, EET, multiheme cytochrome, MHC, porin
Citation: Garber AI, Nealson KH and Merino N (2024) Large-scale prediction of outer-membrane multiheme cytochromes uncovers hidden diversity of electroactive bacteria and underlying pathways. Front. Microbiol. 15:1448685. doi: 10.3389/fmicb.2024.1448685
Received: 13 June 2024; Accepted: 05 August 2024;
Published: 01 October 2024.
Edited by:
F. Wayne Outten, University of South Carolina, United StatesReviewed by:
Julea Butt, University of East Anglia, United KingdomHiroyuki Arai, The University of Tokyo, Japan
Copyright © 2024 Garber, Nealson and Merino. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arkadiy I. Garber, arkadiyg@usc.edu; Nancy Merino, merino4@llnl.gov