Distribution and diversity of ROS-generating enzymes across the animal kingdom, with a focus on sponges (Porifera)
BMC Biology volume 20, Article number: 212 (2022)
Reactive derivatives of oxygen (reactive oxygen species; ROS) are essential in signalling networks of all aerobic life. Redox signalling, based on cascades of oxidation–reduction reactions, is an evolutionarily ancient mechanism that uses ROS to regulate an array of vital cellular processes. Hydrogen peroxide (H2O2) and superoxide anion (O2•−) are employed as signalling molecules that alter the oxidation state of atoms, inhibiting or activating gene activity. Here, we conduct metazoan-wide comparative genomic assessments of the two enzyme families, superoxide dismutase (SOD) and NADPH oxidases (NOX), that generate H2O2 and/or O2•− in animals.
Using the genomes of 19 metazoan species representing 10 phyla, we expand significantly on previous surveys of these two ancient enzyme families. We find that the diversity and distribution of both the SOD and NOX enzyme families comprise some conserved members but also vary considerably across phyletic animal lineages. For example, there is substantial NOX gene loss in the ctenophore Mnemiopsis leidyi and divergent SOD isoforms in the bilaterians D. melanogaster and C. elegans. We focus particularly on the sponges (phylum Porifera), a sister group to all other metazoans, from which these enzymes have not previously been described. Within Porifera, we find a unique calcium-regulated NOX, the widespread radiation of an atypical member of CuZnSOD named Rsod, and a novel endoplasmic reticulum MnSOD that is prevalent across aquatic metazoans.
Considering the precise, spatiotemporal specificity of redox signalling, our findings highlight the value of expanding redox research across a greater diversity of organisms to better understand the functional roles of these ancient enzymes within a universally important signalling mechanism.
Reactive oxygen species (ROS) are widely known as toxic derivatives of oxygen that induce oxidative damage. Yet, aerobic organisms across the tree of life depend on ROS as signalling molecules for a vast array of life-sustaining cellular functions [1, 2]. Indeed, it is the readiness with which ROS react with neighbouring molecules that makes them both potentially toxic and particularly well-suited to form an integral part of the redox signalling network.
Redox is the loss and gain of electrons during oxidation–reduction reactions, which can both generate and consume ROS . Thus, redox signalling may be generalised as a chain of redox relays that transfer electrons from one chemical species to the next, resulting in a cascade that carries a signal from receipt to response . In turn, cellular ROS is tightly regulated by antioxidants that mediate localised redox states by scavenging ROS in reduction reactions and thus also function to prevent oxidative damage [4, 5]. The interaction of target proteins with ROS changes the oxidation state of those proteins, causing a conformational change (primarily thiol-based modifications; ) that can either inhibit or activate gene activity. For instance, redox signalling can mediate both the activation and inhibition of specific transcription factors, such as nuclear factor kb (NF-kb) that regulates the inflammatory response of an organism [7, 8].
One of the evolutionary oldest components of the redox signalling network is the antioxidant family known as superoxide dismutase (SOD) that targets a specific type of ROS known as superoxide anion radical (O2•−) . Superoxide is the most readily formed ROS, since it requires only a one-electron reduction from ground state molecular oxygen (O2) [10, 11]. SOD arose prior even to the differentiation of eubacteria from archaea; its first function was to protect anaerobic organisms from ROS toxicity by scavenging O2•− [12, 13]. Interestingly, as a by-product of the reaction to remove O2•−, SOD generates another, less reactive ROS that is today the predominant form of ROS used in redox signalling — this is hydrogen peroxide (H2O2) (Fig. 1A) [5, 10, 11]. Thus, SOD represents one of the most ancient mechanisms of enzymatic ROS generation that is still used today across all domains of life [14, 15]. Three types of SOD have been described in animals, namely cytoplasm SOD1 and extracellular SOD3 that use copper zinc (CuZn) as the electron acceptor and mitochondrial SOD2 that uses manganese (Mn) (Fig. 1B) [9, 12, 15, 16]. SOD is spatially localised within these compartments or extracellularly to limit O2•− diffusion distances by converting it to H2O2, which in turn facilitates the tightly regulated and spatially oriented ROS generation that is required for redox signalling (Fig. 1A; reviewed by 15).
Two processes generate the majority of ROS in animals (reviewed in [5, 18]). First, incomplete reductive processes within the mitochondrial electron transport chain generate a basal level of ROS, in accordance with the rate of metabolic activity (Fig. 1A; reviewed in [19,20,21]). For a long time, this mitochondrial ROS was considered only as a toxic by-product of metabolism, but it is now understood that mitochondrial SOD facilitates tight regulation and spatial specificity of mitochondrial ROS required for redox signalling , reviewed by [15, 22]. Second, specialised transmembrane enzymes, NADPH oxidases (NOX), are known as “professional” ROS generators that produce relatively small and precise spatiotemporal fluctuations in either O2•− or H2O2 when activated by specific signals (Fig. 1A) .
These NOX enzymes perform an essential role in the redox signalling network of animals. In mammals, there are currently seven described NOX subfamilies: NOX1-5 that generate O2•−, and two dual oxidases, Duox1 and Duox2, that can generate both O2•− and H2O2 (Fig. 1C) . All NOX family proteins share a structurally conserved region comprising six ferric reductase transmembrane domains, a single FAD-binding domain and a single ferric reductase NAD-binding domain (Fig. 1C) [25, 26]. NOX1-4 have only these three domain types and rely on additional subunits for their activation, except for NOX4 that is constitutively active. In contrast, NOX5 and Duox are Ca2+-activated via their calcium/calmodulin-sensitive EF-hand domains [27, 28] and Duox also has a single animal haem peroxidase domain (Fig. 1C) [29, 30]. Calcium-activated NOXs are considered the most evolutionarily ancient of the NOX proteins, with ancient subfamilies NOXC/D present in amoeba and algae, and the respiratory burst oxidase homolog (RBOH) present in plants [31, 32]. The activity of NOX allows for specialised generation of O2•− that is tightly controlled and converted into H2O2 by the activity of SODs localised in close proximity (Fig. 1A). For instance, SOD-NOX generation of H2O2 has been specifically implicated in the activation of redox-sensitive transcription factors such as NF-Kb . Despite the fundamental roles of SOD and NOX enzymes, they are generally poorly characterised outside of vertebrate, arthropod and nematode model species. To address this, here, we conduct a comparative genomic assessment of these two major redox signalling enzymes — SOD and NOX — across 19 metazoan species, encompassing 10 different phyla. Through multiple sequence alignment, domain architecture comparisons and phylogenetic analyses, we explore the distribution, diversity and conservation of these enzymes across the metazoan tree. We pay particular attention to sponges (phylum Porifera), by including genomes of five marine and one freshwater species belonging to 4 classes. Given that sponges originated at least 700 million years ago  and are considered the oldest of the extant animal phyletic lineages [35, 36], they are particularly valuable when viewed in a comparative framework. Their unique phylogenetic position as sister to all other animal phyla means that traits shared between sponges and the rest of the animal kingdom may logically be traced back to the last common animal ancestor . By expanding awareness of redox signalling components across the metazoan tree, we provide the foundation for a broader understanding of the functional roles of these ancient enzymes within a universally important signalling mechanism.
Superoxide dismutase (SOD)
Across all 19 metazoan species, we identified a total of 149 unique protein sequences encoding at least one CuZnSOD domain and 50 containing both N and C terminal MnSOD domains. Filtering by characteristic domain structures (Fig. 1B) reduced this number to 113 CuZnSOD, whilst for MnSOD, we retained all 50 sequences (Table 1; Fig. 2; Additional file 2).
Phylogenetic assessment of relationships amongst all 163 putative CuZnSOD and MnSOD genes revealed high support (> 97%) for four main monophyletic clades. These are cytoplasmic CuZnSOD (SOD1), an atypical CuZnSOD called Rsod, mitochondrial MnSOD (SOD2), and reported here for the first time, a predominantly endoplasmic reticulum (ER) localised MnSOD that we refer to as SOD2X (Figs. 1B and 3; Additional file 1 Fig. S1 & S2). In addition, 13 metazoans encode at least one additional CuZnSOD gene that appear to have arisen through multiple independent evolutionary events. For 11 species, these additional CuZnSOD sequences are predicted to localise extracellularly (Fig. 4). Although some of these sequences have been annotated previously as the extracellularly localised subfamily SOD3, they form independent branches that are evolutionary distant from each other rather than forming a monophyletic SOD3 clade (Fig. 3). Additionally, 11 of these sequences do not have a signal peptide, and 10 sequences (from seven species) encode for additional domains that are not typical of the SOD family. Thus, we do not classify these paraphyletic CuZnSOD sequences as SOD3. Notably, the ctenophore Mnemiopsis leidyi encodes eight CuZnSOD genes outside of the Rsod and SOD1 clades, of which six appear to represent a lineage-specific diversification (Fig. 3). Within phylum Porifera (sponges), only Oscarella carmela (class Homoscleromorpha) encodes one further CuZnSOD gene in addition to the subfamilies Rsod and SOD1.
Of the 50 putative MnSOD genes that we identified, 31 fall within a strongly supported (100%) monophyletic clade that includes previously described SOD2 sequences from the Vertebrata, Nematoda and Arthropoda [9, 40]; the majority of these were predicted with DeepLoc-2.0 to be localised to the mitochondria (Fig. 4; Additional file 3). The second strongly supported (99%) MnSOD monophyletic clade comprises 18 metazoan SOD2X sequences; DeepLoc-2.0 predicted 10 of these to be localised to the endoplasmic reticulum, five extracellularly, three to the lysosome/vacuole, two to the cytoplasm and one each to the cell membrane and mitochondria, where four sequences were predicted to localise within multiple subcellular compartments (Fig. 4; Additional file 3). Notably, the freshwater sponge, Ephydatia muelleri, is the only species that lacks a mitochondrial-localised SOD, but it does have a single MnSOD (SOD2X) that is predicted to localise within the lysosome/vacuole (Fig. 4).
We identified 51 putative CuZnSOD sequences with a domain architecture distinct from other described CuZnSODs (Fig. 1B). However, these sequences do share structural similarity with Rsod (Related to SOD), previously identified in Drosophila melanogaster (FBgn0051028; Dmel\CG31028) . Together, these unusual CuZnSOD sequences form a strongly supported (97%) monophyletic clade (Fig. 3). Sequences falling within this monophyletic clade Rsod are distinct from SOD1 and SOD3 in comprising multiple Cu_SOD domains. Using that criterion, we find that 43 of the unusual CuZnSOD are putative Rsod sequences because they encode at least 3 Cu_SOD domains, and 26 comprise at least four. Also unusual, TargetP-2.0 predicted that 35 of these Rsod sequences contain a signal peptide region (Additional file 4). DeepLoc-2.0 predicted that 23 of the putative Rsod sequences are localised to the cell membrane, 21 extracellularly, 5 in the cytoplasm, 4 within the lysosome/vacuole and one within the nucleus; two sequences are localised to multiple compartments (Fig. 4; Additional file 3).
NADPH oxidases (NOX)
Across all 19 metazoan species, we identified a total of 420 unique protein sequences that encode at least one NOX-associated domain, namely NAD-binding_6, FAD-binding_8, or ferric_reduct. Filtering by domain structure (Fig. 1C) characteristic of a NOX reduced this number to 143; of these, a further 25 sequences were removed because they lacked the full “HHHH” motif encoding O2•− production (Table 2; Fig. 2; Additional file 2).
Phylogenetic assessment of relationships amongst the 118 putative NOX genes revealed high support (92–100%) for four monophyletic clades, namely NOX1-3, NOX4, NOX5 and Duox (Fig. 5). Non-metazoan sequences fall into three previously identified groups, namely calcium-regulated NOXC/D, RBOH and subunit-regulated NOXA/B . No single NOX subfamily is ubiquitously represented in all species assessed. In the NOX1-3 clade, seven genes encoding NOX2 sequences in sponges (phylum Porifera) together form a monophyletic clade, indicating diversification of these genes after the sponges diverged from the metazoan stem. However, within this larger group, NOX1 and NOX3 genes were identified only in vertebrates and echinoderms (NOX1), and mammals (NOX3) only (Table 2; Fig. 2). NOX4 was found in 7 of the 19 species, including the non-bilaterian Nematostella vectensis (phylum Cnidaria), but not in poriferans.
NOX2, NOX5 and Duox were the most common, found in 13 of 19 animal species. Indeed, phylum Porifera encode only these common NOX2, NOX5 and Duox subfamilies, although Tethya wilhelma (class Demospongiae) and Sycon ciliatum (class Calcarea) appear to lack NOX2 and NOX5, respectively. Although most metazoans encode only a single Duox, the sponges O. carmella (class Homoscleromorpha), S. ciliatum (class Calcarea) and E. muelleri (class Demospongiae) have multiple that appear to have derived via multiple independent/lineage-specific gene duplication events (Fig. 5). Indeed, species-specific duplication of NOX genes was commonly observed.
We identified four sequences in phylum Porifera, class Demospongiae, that did not group within any of the four main metazoan clades, although they cluster together with strong support (100%) (Fig. 5). Separated by large evolutionary distances, these sequences share some similarity to NOX5, possessing the three core NOX domains and EF-hand regions, but otherwise comprise a novel domain architecture (Fig. 6). We herein refer to these Demospongiae NOXs, described here for the first time, as DspNOX. It is the C-terminal region of DspNOX that has sequence and structural similarity to animal NOX5 (Fig. 6A), including a single calcium-binding EF-hand region. The exception is the freshwater species, E. muelleri that uniquely encodes two Cupredoxin domains at the C-terminal (predicted by Gene3D: CATH Superfamily 184.108.40.2060; Fig. 6E). The N-terminal region of DspNOX comprises three additional structural features, conserved in all four demosponge species. These are a PAS or PAS_9 domain (PF00989; PF13426), a long intrinsically disordered protein (IDP), and in addition to the typical six alpha helixes, an additional 5 (except 3 in E. muelleri) transmembrane domain regions (Fig. 6D, E).
IDP regions are similarly observed in other unusual calcium-regulated NOXs of Amoebozoa and red algae, and the red algae NOXC/D clade that is a sister clade to DspNOX encodes four additional transmembrane regions (Fig. 6B, C) . In the sponge species, T. wilhelma and E. muelleri, all five or three additional transmembrane regions are predicted to comprise a domain of an unknown function (DUF2427; PF10348). In contrast, for Amphimedon queenslandica and Xestospongia bergquistia (class Demospongiae), the lowest e-value score predicted that all five transmembrane domains comprise the eukaryotic cytochrome b561 (Cyb561; PF03188). However, significant matches for DUF2427 localised across the same five transmembrane regions were also observed in both these sponges. Notably, E. muelleri additionally encoded two other DspNOX-like sequences, but these lacked the HHHH motif for O2•− production.
The ROS generators, SOD and NOX (including Duox), both represent ancient and large enzyme families that are widespread across the animal kingdom. Genes encoding both are present in all of the 19 metazoan species that we assessed (Tables 1 and 2; Fig. 2), but it is notable that the gene numbers vary considerably between species (from 1 to 25; Tables 1 and 2). Below we discuss the dynamic content of these gene families across the animal kingdom. For each of the two enzyme families, we draw upon our new findings from the phylum Porifera — considered sister to all other animal phyletic lineages — to provide new insights into the evolution of these gene families since the dawn of the Metazoa more than 700 million years ago.
Superoxide dismutase (SOD)
The SOD family has been evolving for some ~ 2.5 (Ga) billion years of the Earth’s history, under vastly divergent selective pressures and with prolific cross-domain horizontal gene transfer events [15, 38, 41, 42]. In the animal kingdom alone, we observe various metalloforms and their subfamilies are not distributed equally across the metazoan tree, and neither are their localisations to various subcellular compartments. Consistent with previous analyses based on smaller numbers of animal taxa, we find that CuZn SOD1 and Mn SOD2 are the most conserved widespread SOD subfamily members [40, 43, 44]. However, the broader taxonomy that we present here extends the known SOD diversity by revealing the widespread metazoan prevalence of an atypical CuZnSOD, named Rsod, and an additional MnSOD subfamily member, SOD2X.
Evolution of the SOD family
Amongst the CuZnSODs, it is widely accepted that cytoplasmic SOD1 arose first and later gave rise to the extracellular SOD3 with the addition of signal peptides [12, 40, 44, 45].Our analysis reveals that SOD1 forms a strongly supported monophyletic group comprising metazoan, fungi, plant, algae and amoebozoan species (93%), whilst SOD3 is paraphyletic. We found metazoan species commonly have at least one additional CuZnSOD outside of SOD1 or Rsod, but these sequences are diverse and do not necessarily classify as SOD3; not all sequences encode a signal peptide, and many possess additional domains not typically associated to CuZnSOD. Indeed, together, we found these unclassified CuZnSOD sequences with those previously annotated as SOD3 formed six independent monophyletic groups, and three further independent branches (Fig. 3).
Thus, our results strongly support the hypothesis that SOD3 has arisen via multiple, independent evolutionary events in different metazoan lineages [40, 44, 45]. Moreover, some of these “SOD3” branches are separated by considerable distances, and many are not localised extracellularly as is typically described for SOD3 in vertebrates, nematodes and arthropods [40, 44] (Fig. 4). Considering that localisation of SOD dictates its signalling pathway involvement, we suggest that these unclassified CuZnSOD sequences likely exhibit diverse functionality across the metazoan phyla. Given this, we further suggest that the term SOD3 is not sufficient to capture the complete diversity of paraphyletic CuZnSODs that have independently arisen, on multiple occasions across the metazoans.
Interestingly, most of the sponges we survey lack any additional CuZnSOD SOD3-like gene, which we suggest may be functionally substituted by the extracellular Rsod that is encoded by all sponge species except for O. carmela that does encode an extracellular SOD3. As with SOD1, the Rsod monophyletic clade comprises representatives from multiple kingdoms and thus supports previous suggestions that Rsod belongs to an ancient group of proteins . Based on these observations, we propose that Rsod predates SOD3 and that this explains the predominant absence of SOD3 in phylum Porifera where Rsod is common.
Of the various SOD metalloforms, it has most commonly been hypothesised that MnSOD (homologous to non-animal FeSOD) is ancestral, with CuZnSOD arising later, after the great oxidation event (reviewed by [12, 15, 42]). Alternatively,  recently proposed instead that CuZnSOD first arose in the Archean, predating Fe/MnSOD that subsequently appeared much later in the mid-Proterozoic. In support of this, here, we find SOD1 is the closest monophyletic clade to CuZnSODs from Archean cyanobacterial lineages Pseudanabaena sp. and Gloeobacter violaceus  (Fig. 3). Moreover, unlike MnSOD2 where the metazoans form an independent monophyletic group (87%), the SOD1 clade comprises genes from multiple kingdoms, reflecting an ancient evolutionary origin that predates separation of the eukaryotic kingdoms (Fig. 3). SOD1 sequences are also more variably localised; they are found within the nucleus, ER, mitochondria, cell membrane and extracellularly across five metazoan species (Fig. 4; Additional file 3). Conversely, Mn SOD2 are consistently localised within the mitochondria, indicating a comparatively greater degree of conservation (Fig. 3).
We must also consider the prevalence of cross-domain horizontal gene transfer. Evolutionary analyses of bacterial SOD indicate CuZnSOD was lost as cyanobacteria diversified within new ecological niches  and then later spread via horizontal transmission between non-cyanobacterial phyla likely multiple times, resulting in the distribution found today [38, 42]. Thus, it seems possible the distribution of metazoan SODs may similarly have arisen via horizontal acquisition across the domains.
An atypical CuZnSOD — Rsod — is widespread in the animal kingdom
We reveal a widespread metazoan presence of an atypical CuZnSOD, phylogenetically distant from both SOD1 and SOD3, that is largely localised either extracellularly or membrane bound (green shape, Fig. 3). These atypical CuZnSOD sequences share sequence and structural similarity to a D. melanogaster sequence named Rsod, “Related to SOD” (Dmel\CG31028, FBgn0051028), clustering together within a single, well-supported monophyletic clade  (Fig. 3). The structure of Rsod is unique from other known CuZnSODs, comprising multiple CuZnSOD domains (between 2 and 6) often with signal peptide and/or transmembrane regions (Fig. 1B). Rsod is considered to belong to an ancient group of proteins containing CuZnSOD homology domains , but has been little explored since its initial description, and its function is not yet known.
To date, Rsod sequences have been identified across the animal kingdom, including in insects, fish and the urochordate Ciona intestinalis and also in fungi [40, 45], BmSOD6, Accession Nos. LC229593 from 45]. Here, we additionally identify Rsod homologues within six sponges, echinoderms, cephalochordate (Branchiostoma floridae) and two vertebrates, indicating the widespread prevalence of this atypical CuZnSOD (Fig. 3). Notably, we show here for the first time that Rsod is particularly abundant and diversified within phylum Porifera, and especially within T. wilhelma (class Demospongiae) and S. ciliatum (class Calcarea). Subcellular localisation predictions indicate that most Rsod sequences we identified are either extracellular (21 sequences) or cell membrane bound (23 sequences), except five sequences localised within the cytoplasm in sponges (1 in A. queenslandica, X. bergquistia, E. muelleri and 2 in T. wilhelma), four in the vacuole/lysome (1 in B. floridae, and T. wilhelma and 2 in S. ciliatum) and one within the nucleus (Fig. 4; Additional file 3). Intriguingly, five species of Porifera, the tunicate C. intestinalis and Ascomycota fungi Phaeosphaeria nodorum, that do not encode any other extracellular CuZnSOD (e.g. SOD3), all encode both extracellular and membrane forms of Rsod; the exception is E. muelleri that encodes a cytoplasmic Rsod instead of membrane bound. In contrast, O. carmela, Acanthaster planci, Danio rerio, Xenopus tropicalis, D. melanogaster and Bombyx mori (BmSOD6; Kobayashi et al. 2019), all of which do encode extracellular SOD3, have only the cell membrane localised form (Fig. 4; Additional file 3). Thus, we suggest the extracellular Rsod in these species may functionally replace an extracellular SOD3. That said, B. floridae and Strongylocentrotus purpuratus comprise only extracellular localised Rsod, despite also encoding extracellular SOD3 (Fig. 4).
To date, the possible function of Rsod remains unknown. Kobayashi et al. assessed the responsiveness of different SOD genes within silk moth (B. mori) under various oxidative stressors and found Rsod gene (BmSOD6) is mostly expressed within the testes on day 3 of fifth instar larvae, indicating a role in removing ROS generated during spermatogenesis . Additionally, honeybees (Apis mellifera) exposed to caging stress for up to 4 weeks showed a constant upregulation of only two genes, namely Rsod and another antioxidant thioredoxin-1 (Trx-1) . In crown of thorns starfish, A. planci, maintained in captivity, one Rsod homologue (gbr.190.13.t1) is upregulated in three tissues (skin, tube feet and spines), and a second homologue (gbr.190.14.t1) in skin only, compared to wild populations .
The NADPH oxidase family
Much of the current understanding of NOX function is based on mammals, but variation in NOX subfamily membership across the metazoan and beyond indicates NOX functions are also likely to be variable. In accordance with 30 and 31, we find the calcium-regulated subfamilies have the earliest origin and are the widest distributed, followed by subunit-regulated NOXs that first appeared in fungi and amoeba, and for which we find the sponges (Porifera) have the earliest metazoan NOX2 co-ortholog (Fig. 5). However, neither the calcium- nor subunit-regulated NOXs are universally found in all metazoans, nor is any single NOX subfamily member. This is consistent with previous analyses based on smaller numbers of taxa [31, 32] showing large variability in gene number across the metazoan tree, indicative of lineage-specific gene duplication and apparent gene loss involving all NOX gene family members.
There are no core NOX subfamilies present across the Metazoa
NOX gene duplication and losses have been commonly documented but are not yet well understood . Whilst some species encode multiple (up to 7) sequences for a single NOX, others lack the subfamily entirely (Table 2). Different NOXs are documented to function with precise, spatiotemporal specificity (Nathan and Cunningham-Bussel 2013; Sies and Jones 2020). Thus, it is somewhat surprising to find extensive gene losses in M. leidy that comprises just a single NOX5, and Caenorhabditis elegans and D. melanogaster that each comprise only two calcium-regulated NOXs (Table 2; Fig. 2). Whilst plants similarly encode a single subfamily, RBOHs that are NOX5-like homologues, species possess up to 10 different members of these in the subfamily [31, 50]. This raises the question of how species may compensate the absence of key subfamily members, NOX2, NOX5 and Duox, that each have disparate described signalling roles [27, 29, 51].
One possibility is that where a single NOX subfamily is represented by multiple genes, each gene exhibits distinct spatio-temporal specificity. For instance, S. ciliatum and C. intestinalis encode multiple Duox enzymes but lack NOX5. Considering that both Duox and NOX5 may be activated by calcium-sensitive EF-hand domains, the multiple Duox genes may compensate for the missing NOX5 function, as has been predicted in rodents . That said, N. vectensis (Cnidaria) lacks both these calcium-regulated NOXs but does encode a subunit-regulated NOX4 that is absent in other basal metazoans. An alternative possibility is that missing subfamily members may be compensated for by the existing NOX genes having a broader range of functions that depend on their spatiotemporal localisation. As an example, NOX2 was documented first in phagocytic cells of sea urchin and mouse oocytes, where it was determined to function in the “oxidative burst” during phagocytosis [53, 54] but since has been identified in diverse cell types, and with diverse functions [reviewed in 23].
The ecological niche of an organism also will impact its redox states [2, 55] and thus likely the redox machinery it requires. The ctenophore, M. leidyi that encodes a single NOX5, exhibits daily vertical migrations [56, 57], perhaps in response to high irradiance levels, that are known to influence localised redox states . We propose that vertical migratory species such as M. leidyi may also be able to mediate O2•− production via timing their daily migrations. Indeed,  show that antioxidant protection in cetaceans differs between shallow-diving and deep-diving habitats, as reflected by their O2•− production and antioxidant levels . The metabolic activity of vertical migrations may sufficiently influence ROS generation in the mitochondrial ETC, endoplasmic reticulum or peroxisome to compensate fewer NOX enzymes. For example, M. leidyi encodes a comparatively wide range of CuZnSODs that are localised to those ROS-generating subcellular compartments, which may reflect greater ROS generation during cellular and metabolic processes (Table 1; Figs. 2 and 3).
Class Demospongiae encode a novel, calcium-regulated NOX
Within class Demospongiae (phylum Porifera), we identify here for the first time a structurally unique and phylogenetically distant NOX, which we call DspNOX (Fig. 6). The C-terminal region of DspNOX shares sequence and structural similarity with animal NOX5, comprising the core NOX domain features, and except for in the freshwater demosponge E. muelleri, is EF-hand calcium sensitive. DspNOX is ~ 700 residues longer than NOX5 and red algae NOXC/D and has a unique domain structure towards the N-terminal region that likely confers additional functionality and modes of activation.
Specifically, the N-terminal region of DspNOX comprises three novel elements (Fig. 6D, E). First, there is a PAS domain (PF00989) that functions as a sensory unit for diverse signals, including chemoreception, redox, photons and voltage, and thus acts as a highly versatile signal transducer . It is well known as a light sensor, helping entrain the canonical circadian clock, but is also important in abiotic stress responses and innate immunity [62, 63]. The presence of both EF-hand and PAS suggests that, in addition to calcium, DspNOX could be regulated directly by abiotic factors. Second, DspNOX contains a long intrinsically disordered protein (IDP) region, defined as 30 or more consecutive disordered residues [64, 65]. The nature of IDP interactions (high specificity, but low affinity) increases a protein’s possible interactions and functional plasticity, making them well suited for signalling and regulatory functions . Third, DspNOX contains five (or three within E. muelleri) additional transmembrane regions encoded within either a domain of unknown function DUF2427 (PF10348) or in the eukaryotic cytochrome b561 (Cyt-b561; PF03188).
Cyt-b561-containing proteins are an enzyme family of transmembrane, ascorbate-dependant oxidoreductases, most well known for their role in recycling ascorbate (i.e. vitamin C) via electron transfer from two heme b groups across the membrane . Cyt-b561 enzymes may also be involved with iron metabolism, first described in the mammalian duodenal Cyt b561 (Dcytb), essential in the uptake of dietary nonheme iron (Fe+3) . Dcytb uses ascorbate in the cytoplasm as an electron donor to reduce either Fe+3 into soluble, ferrous iron (Fe+2) or monodehydroascorbate (MDHA), depending on substrate availability . Considering that the fenton reaction of Fe+2 with O2•− generates the highly reactive hydroxyl radical (•OH), it is particularly interesting that DspNOX has the potential to produce both these molecules. Indeed, because of the significant damage that •OH causes if not contained, organisms have evolved under strong selection to keep ferrous iron and ROS apart . Thus, we suggest that DspNOX may provide an adaptive mechanism to keep Fe+2 away from O2•− if the ASC binding sites are on the same side as O2•− generation.
NOX or pre-NOX genes are distributed right across all the eukaryotic supergroups, except for Rhizaria [17, 25]. It is hypothesised that an enzyme similar to the red algal NOXD gave rise to animal NOX5, via acquisition of calcium-binding motifs [31, 70, 71]. Interestingly, similar to DspNOX, red algal NOXD and amoebozoan NOXC also comprise long IDP regions, which thus may indicate ancestral calcium-regulated NOXs (Fig. 6B, C). Red algae similarly encode four additional transmembrane domains, although these are located between two NADPH-binding site sub-regions and thus predicted to function only as an anchor to the membrane . Results from our phylogenetic analysis reveal that red algae NOXD is sister to DspNOX (Fig. 5). Although the possibility of long-branch attraction cannot be discounted, the described structural similarities of these two proteins and the consistent placement of NOXC/D with previous assessments  together suggest this relationship is unlikely to be an artefact. Thus, we propose that DspNOX also represents an ancient NOX lineage that is likely a lineage-specific innovation within the class Demospongiae (phylum Porifera).
Across 10 metazoan phyla, we observe a very high level of conservation within certain SOD and NOX subfamilies. In particular, we find that SOD1 and SOD2, as well as NOX2, NOX5 and Duox, are widely distributed across the metazoans and form strongly supported monophyletic clades (Figs. 3 and 5). However, our broader coverage of metazoan phyla significantly expands the known diversity and distribution of these enzyme families, revealing several cases of gene loss, gene duplication and lineage-specific expansions. Importantly, our inclusion of sponges (phylum Porifera), considered sister to all other animal phyletic lineages, allows us to logically trace the distribution and characteristics of SOD and NOX families to the last common animal ancestor.
In sponges, we find all three of the major animal NOX subfamilies — NOX2, NOX5 and Duox — indicating the likely presence of these in the last common animal ancestor. That said, the ctenophore M. leidyi encodes only NOX5; under the alternative hypothesis that phylum Ctenophora is sister to all other animal phyletic lineages, then NOX5 but not Duox or NOX2 were likely present within the last common animal ancestor. NOX4 is absent from sponges and appears to have arisen later, after sponges diverged from the animal stem, but before the divergence of the cnidarian phyletic lineage. We also reveal the presence of a novel, calcium-regulated NOX, namely DspNOX, in class Demospongiae only. DspNOX, together with previously identified, unusual red algae and amoebozoan NOXC/D, reflects the broader structural and likely functional diversity of calcium-regulated NOXs amongst diverse species.
MnSODs comprise two strongly supported metazoan subfamilies, namely mitochondrial SOD2 and another predominantly extracellular clade, SOD2X, that we identified within all sponges. The CuZnSODs also comprise two strongly supported metazoan subfamilies, namely Rsod and SOD1. In contrast, SOD3 is paraphyletic and its absence in sponges supports the hypothesis that SOD3 arose sometime later after SOD1. The atypical CuZnSOD, called Rsod, is found within all six sponges and likely represents an ancient protein group that also predates SOD3. The SOD1 clade includes genes from all six sponges, and also genes from amoebozoans, plants and red algae, indicating a much deeper origin, for SOD1, earlier in eukaryotic evolution. Evidently, the evolution of the SOD family is complex with much remaining unclear. However, the diversity and lineage-specific divergences we observe across the animal kingdom illustrate the presence of highly specialised SODs in redox signalling networks of diverse animals with diverse ecologies.
The comparatively greater metazoan diversification and subfamily distribution of SOD than NOX is indicative of the deeper evolutionary origin of the former. However, both these enzyme families are integral to ancient systems of redox signalling and oxidative defence. Thus, the observed evolutionary relationships also reflect the considerably variable, species-specific life-history trade-offs between redox signalling and ROS toxicity.
Enzyme identification and subfamily classification
To search for gene sequences encoding candidate members of the SOD and NOX enzyme families, we assessed protein-coding sequences for 19 metazoan species of 10 phyla, with representation from each of the major clades, namely non-bilaterian, protostome bilaterians and deuterostome bilaterians (Additional file 1: Table S1). Coding sequences were scanned against the Pfam database using hmmscan in HMMR v3.1b2 (hmmer.org) for sequences encoding domains specific to each enzyme family and their respective subfamilies (Fig. 1B, C) , and the number of and position of all identified domains was determined. For all identified candidate gene sequences, protein subcellular localisation was predicted using DeepLoc-2.0, https://services.healthtech.dtu.dk/service.php?DeepLoc-2.0 [39, 73], that uses protein sequences as input for the Neural Networks algorithm trained on Uniprot proteins with experimental evidence. The algorithm incorporates the importance (“attention”) of particular amino acids and those neighbouring within the region. Positions in the sequence with high “attention” give more weight to the final prediction of the model. DeepLoc-2.0 is able to predict proteins that are located in more than one compartment. The criterion for deciding subcellular localisation is based on probability scores surpassing thresholds (set by 72) or if no score crosses the threshold, the label closest to the threshold is chosen.
The methodology for enzyme identification and subcellular localisation was cross-validated by comparing the number and type of SOD and NOX genes identified through our analysis with those that have previously been described. Candidate SOD enzyme sequences were assigned to either CuZnSOD or MnSOD based on the presence of domains for either copper/zinc binding (PF00080) or manganese/iron SOD C- and N-terminals (PF02777 and PF00081), respectively (Fig. 1B). For 50 CuZnSOD encoding sequences that clustered separately from the previously described SOD1 or SOD3 lineages, we predicted the presence of signal peptides using TargetP-2.0, https://services.healthtech.dtu.dk/service.php?TargetP-2.0 .
Candidate NOX enzyme sequences encoding the domains ferric reductase (PF01794) and at least one of either the FAD- (PF08022) or NAD-binding (PF08030) domains (Fig. 1C) were retained and scanned for the presence of the conserved NOX motif, His101, His115, His209, His222 (based on Homo sapiens NOX2 numbering), required for O2•− generation . NOX5 and Duox encoding sequences were initially classified based on the presence of EF-hand regions (both NOX5 and Duox) and of the domain for animal haem peroxidase (PF03098; Duox only; Fig. 1B). For four sequences identified in phylum Porifera, class Demospongiae, EF-hand regions were not initially identified in the hmmscan used as described above, but subsequently were predicted based on Gene3D ontology (http://www.cathdb.info/), and later manually assessed following MAFFT version 7.455 alignment of EF-hand domain regions from Pfam seed sequences. Additionally, for these same four sequences, intrinsically disordered protein (IDP) regions were predicted using IUPred3 (https://iupred.elte.hu/).
Multiple sequence and phylogenetic assessment.
Identified candidate sequences within each enzyme family were aligned using MAFFT version 7.455 ([76, 77], https://mafft.cbrc.jp/alignment/software/) with default parameters and visualised in the multiple sequence alignment editor, AliView .
To assess phylogenetic relationships, alignments were manually edited in AliView v1.27 (; https://ormbunkar.se/aliview), removing regions containing more than 50% gaps, then imported to IQ-TREE  to construct maximum likelihood trees using ultrafast bootstrap , based on 1000 bb, and the most appropriate evolutionary model as identified by ModelFinder . Models identified and used to construct each gene family tree were as follows: SOD, WAG + R6, NOX, Dayhoff + F + R6, MnSOD only, WAG + I + G4, and CuZnSOD only, WAG + R5. Phylogenetic trees were first visualised in iTOL v.6.2.1  before importing and annotating within Adobe Illustrator. Classification of SOD and NOX gene subfamilies was inferred from the relative placing of putative sequences within known subfamily clades of phylogenetic trees. To provide evolutionary context to the metazoan phylogenetic relationships, we also included sequences obtained from organisms outside the metazoans that have previously been described (Additional file 1: Table S2).
Availability of data and materials
All data generated or analysed during this study are included in this published article, its supplementary information files and publicly available repositories.
Copper zinc superoxide dismutase
Iron superoxide dismutase
Intrinsically disordered protein
Manganese superoxide dismutase
Nuclear factor kb
Respiratory burst oxidase homolog
Reactive oxygen species
Related to superoxide dismutase
Halliwell B, Reactive species and antioxidants. Redox biology is a fundamental theme of aerobic life. Plant Physiol. 2006;141:312–22. https://doi.org/10.1104/pp.106.077073.
Sies H, Jones DP. Reactive oxygen species (ROS) as pleiotropic physiological signalling agents. Nat Rev Mol Cell Biol. 2020;21(7):363–83. https://doi.org/10.1038/s41580-020-0230-3.
Sobotta MC, Liou W, Stöcker S, Talwar D, Oehler M, Ruppert T, Scharf AN, Dick TP. Peroxiredoxin-2 and STAT3 form a redox relay for H2O2 signaling. Nat Chem Biol. 2015;11:64–70.
Lambeth JD, Neish AS. Nox enzymes and new thinking on reactive oxygen: a double-edged sword revisited. Annu Rev Pathol: Mechanisms of Disease. 2014;9:119–45. https://doi.org/10.1146/annurev-pathol-012513-104651.
Sies H. Hydrogen peroxide as a central redox signaling molecule in physiological oxidative stress: oxidative eustress. Redox Biol. 2017;11:613–9. https://doi.org/10.1016/j.redox.2016.12.035.
Zeida A, Trujillo M, Ferrer-Sueta G, Denicola A, Estrin DA, Radi R. Catalysis of peroxide reduction by fast reacting protein thiols. Chem Rev. 2019;119:10829–55.
Schreck R, Rieber P, Baeuerle PA. Reactive oxygen intermediates as apparently widely used messengers in the activation of the NF-kappaB transcription factor and HIV-1. EMBO J. 1991;10:2247–58.
Oliveira- Marques V, Marinho HS, Cyrne L, Antunes F. Role of hydrogen peroxide in NF-kappaB activation: from inducer to modulator. Antioxid Redox Signal. 2009;11:2223–43.
Zelko IN, Mariani TJ, Folz RJ. Superoxide dismutase multigene family: a comparison of the CuZn-SOD (SOD1), Mn-SOD (SOD2), and EC-SOD (SOD3) gene structures, evolution, and expression. Free Radic Biol Med. 2002;33(3):337–49. https://doi.org/10.1016/s0891-5849(02)00905-x.
Rhee SG. Redox signaling: hydrogen peroxide as intracellular messenger. Exp Mol Med. 1999;31:53–9.
D’Autreaux B, Toledano MB. ROS as signalling molecules: mechanisms that generate specificity in ROS homeostasis. Nat Rev Mol Cell Biol. 2007;8:813–24.
Miller AF. Superoxide dismutases: ancient enzymes and new insights. FEBS Lett. 2012;586:585–95. https://doi.org/10.1016/j.febslet.2011.10.048.
Ślesak I, Ślesak H, Kruk J. Evolution of life on earth: in silico comparative analysis of biochemical pathways. Astrobiology. 2012;12:775–84. https://doi.org/10.1089/ast.2011.0704.
Sheng Y, Abreu IA, Cabelli DE, Maroney MJ, Miller AF, Teixeira M, Valentine JS. Superoxide dismutases and superoxide reductases. Chem Rev. 2014;114(7):3854–918. https://doi.org/10.1021/cr4005296.
Case AJ. On the origin of superoxide dismutase: an evolutionary perspective of superoxide-mediated redox signaling. Antioxidants (Basel). 2017;6(4):82. https://doi.org/10.3390/antiox6040082.
Potter SZ, Valentine JS. The perplexing role of copper-zinc superoxide dismutase in amyotrophic lateral sclerosis (Lou Gehrig’s disease). J Biol Inorg Chem. 2003;8(4):373–80. https://doi.org/10.1007/s00775-003-0447-6.
Zhang X, Krause KH, Xenarios I, Soldati T, Boeckmann B. Evolution of the ferric reductase domain (FRD) superfamily: modularity, functional diversification, and signature motifs. PLoS One. 2013;8(3):e58126. https://doi.org/10.1371/journal.pone.0058126.
Fang FC. Antimicrobial reactive oxygen and nitrogen species: concepts and controversies. Nat Rev Microbiol. 2004;2:820–32.
Mittler R. ROS are good. Trends Plant Sci. 2016;22(1):11–9. https://doi.org/10.1016/j.tplants.2016.08.002.
Simon NML, Litthauer S, Jones MA, Dodd AN (2019) Interactions between circadian rhythms, ROS and Redox. In: Panda S., Yamamoto Y. (eds) Redox homeostasis in plants. signaling and communication in plants. Springer, Cham. https://doi.org/10.1007/978-3-319-95315-1_4
Murphy MP, Holmgren A, Larsson NG, Halliwell B, Chang CJ, Kalyanaraman B, Rhee SG, Thornalley PJ, Partridge L, Gems D, Nyström T, Belousov V, Schumacker PT, Winterbourn CC. Unraveling the biological roles of reactive oxygen species. Cell Metab. 2011;13:361–6.
Foyer CH, Ruban AV, Noctor G. Viewing oxidative stress through the lens of oxidative signalling rather than damage. Biochem J. 2017;474(6):877–83. https://doi.org/10.1042/BCJ20160814.
Lambeth JD. NOX enzymes and the biology of reactive oxygen. Nat Rev Immunol. 2004;4:181–9. https://doi.org/10.1038/nri1312.
Buvelot H, Jaquet V, Krause KH. Mammalian NADPH oxidases. Methods Mol Biol. 1982;2019:17–36. https://doi.org/10.1007/978-1-4939-9424-3_2.
Sumimoto H. Structure, regulation and evolution of Nox-family NADPH oxidases that produce reactive oxygen species. FEBS J. 2008;275:3249–77. https://doi.org/10.1111/j.1742-4658.2008.06488.x.
Magnani F, Nenci S, Millana Fananas E, Ceccon M, Romero E, Fraaije MW, Mattevi A. Crystal structures and atomic model of NADPH oxidase. Proc Natl Acad Sci U S A. 2017;114:6764–9. https://doi.org/10.1073/pnas.1702293114.
Touyz RM, Anagnostopoulou A, Rios F, Montezano AC. Camargo LL (2019) NOX5: Molecular biology and pathophysiology. Exp Physiol. 2019;104:605–16. https://doi.org/10.1113/EP086204.
Fañanás EM, Todesca S, Sicorello A, Masino L, Pompach P, Magnani F, Pastore A, Mattevi A. On the mechanism of calcium-dependent activation of NADPH oxidase 5 (NOX5). FEBS J. 2020;287(12):2486–503. https://doi.org/10.1111/febs.15160.
van der Vliet A, Danyal K, Heppner DE. Dual oxidase: a novel therapeutic target in allergic disease. Br J Pharmacol. 2018;175(9):1401–18. https://doi.org/10.1111/bph.14158.
Ameziane-El-Hassani R, Morand S, Boucher J-L, et al. Dual oxidase-2 has an intrinsic Ca2þ-dependent H2O2-generating activity. J Biol Chem. 2005;280:30046–54.
Kawahara BT, Quinn MT, Lambeth JD. Molecular evolution of the reactive oxygen-generating NADPH oxidase (Nox/Duox) family of enzymes. BMC Evol Biol. 2007;7:109.
Gandara ACP, Torres A, Bahia AC, Oliveira PL, Schama R. Evolutionary origin and function of NOX4-art, an arthropod specific NADPH oxidase. BMC Evol Biol. 2017;17(1):92. https://doi.org/10.1186/s12862-017-0940-0.
Mumbengegwi DR, Li Q, Li C, Bear CE, Engelhardt JF. Evidence for a superoxide permeability pathway in endosomal membranes. Mol Cell Biol. 2008;28(11):3700–12. https://doi.org/10.1128/MCB.02038-07.
Erwin DH, Laflamme M, Tweedt SM, Sperling EA, Pisani D, Peterson KJ. The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science. 2011;334(6059):1091–7. https://doi.org/10.1126/science.1206375 PMID: 22116879.
Simion P, Philippe H, Baurain D, Jager M, Richter DJ, Di Franco A, Roure B, Satoh N, Quéinnec E, Ereskovsky A, Lapébie P, Corre E, Delsuc F, King N, Wörheide G, Manuel M. A large and consistent phylogenomic dataset supports sponges as the sister group to all other animals. Curr Biol. 2017;27:1–10. https://doi.org/10.1016/j.cub.2017.02.031.
Nielsen C. Early animal evolution: a morphologist’s view. R Soc Open Sci. 2019;6(7):190638. https://doi.org/10.1098/rsos.190638.
Srivastava M, Simakov O, Chapman J, et al. The Amphimedon queenslandica genome and the evolution of animal complexity. Nature. 2010;466:720–6. https://doi.org/10.1038/nature09201.
Boden JS, Konhauser KO, Robbins LJ, et al. Timing the evolution of antioxidant enzymes in cyanobacteria. Nat Commun. 2021;12:4742. https://doi.org/10.1038/s41467-021-24396-y.
Armenteros JJA, Sønderby CK, Sønderby SK, Nielsen H, Winther O. DeepLoc: prediction of protein subcellular localization using deep learning. Bioinformatics. 2017;33:3387–95. https://doi.org/10.1093/bioinformatics/btx431.
Landis GN, Tower J. Superoxide dismutase evolution and life span regulation. Mech Ageing Dev. 2005;126:365–79. https://doi.org/10.1016/j.mad.2004.08.012.
Inupakutika MA, Sengupta S, Devireddy AR, Azad RK, Mittler R. The evolution of reactive oxygen species metabolism. J Exp Bot. 2016;67(21):5933–43. https://doi.org/10.1093/jxb/erw382.
Harada M, Akiyama A, Furukawa R, Yokobori SI, Tajika E, Yamagishi A. Evolution of superoxide dismutases and catalases in cyanobacteria: occurrence of the antioxidant enzyme genes before the rise of atmospheric oxygen. J Mol Evol. 2021;89(8):527–43. https://doi.org/10.1007/s00239-021-10021-5.
Fink RC, Scandalios JG. Molecular evolution and structure—function relationships of the superoxide dismutase gene families in angiosperms and their relationship to other eukaryotic and prokaryotic superoxide dismutases. Arch Biochem Biophys. 2002;99:19–36.
Xu L, Yang J, Xu M, Shan D, Wu Z, Yuan D. Speciation and adaptive evolution reshape antioxidant enzymatic system diversity across the phylum Nematoda. BMC Biol. 2020;18(1):181. https://doi.org/10.1186/s12915-020-00896-z.
Corona M, Robinson GE. Genes of the antioxidant system of the honey bee: annotation and phylogeny. Insect Mol Biol. 2006;15:687–701. https://doi.org/10.1111/j.1365-2583.2006.00695.x.
Blank CE, Sanchez-Baracaldo P. Timing of morphological and ecological innovations in the cyanobacteria - a key to understanding the rise in atmospheric oxygen. Geobiology. 2010;8:1–23.
Kobayashi Y, Nojima Y, Sakamoto T, Iwabuchi K, Nakazato T, Bono H, Toyoda A, Fujiyama A, Kanost MR, Tabunoki H. Comparative analysis of seven types of superoxide dismutases for their ability to respond to oxidative stress in Bombyx mori. Sci rep. 2019;9:2170. https://doi.org/10.1038/s41598-018-38384-8.
Alburaki M, Karim S, Lamour K, Adamczyk J, Stewart SD. RNA-seq reveals disruption of gene regulation when honeybees are caged and deprived of hive conditions. J Exp Biol. 2019;222:jeb207761. https://doi.org/10.1242/jeb.207761.
Morin M, Jönsson M, Wang CK, Craik DJ, Degnan SM, Degnan BM. bioRxiv. https://doi.org/10.1101/2022.07.21.501052. Accessed 31st July 2022.
Bedard K, Lardy B, Krause KH. NOX family NADPH oxidases: not just in mammals. Biochimie. 2007;89:1107–12.
Moghadam ZM, Henneke P, Kolter J. From flies to men: ROS and the NADPH oxidase in phagocytes. Front Cell Dev Biol. 2021;9:628991. https://doi.org/10.3389/fcell.2021.628991.
Bedard K, Jaquet V, Krause K-H. NOX5: from basic biology to signaling and disease. Free Radic Biol Med. 2012;52:725–34.
Shapiro BM. The control of oxidant stress at fertilization. Science. 1991;252:533–6.
Nasr-Esfahani MM, Johnson MH. The origin of reactive oxygen species in mouse embryos cultured in vitro. Development. 1991;113:551–60.
Niedzwiecki MM, Walker DI, Vermeulen R, Chadeau-Hyam M, Jones DP, Miller GW. The exposome: molecules to populations. Annu Rev Pharmacol Toxicol. 2019;59:107–27. https://doi.org/10.1042/BJ20081386.
Vereshchaka AL. Features of the microscale distribution of the gelationous macroplankton in the Black Sea off Gelendzhik. Oceanology. 2002;42(1):91–7.
Finenko GA, Kideys AE, Anninsky BE, Shiganova TA, Roohi A, Tabari MR, Rostami H, Bagheri S. Invasive ctenophore Mnemiopsis leidyi in the Caspian Sea: feeding, respiration, reproduction and predatory impact on the zooplankton community. Mar Ecol Prog. 2006;314:171–85.
Haraldsson M, Båmstedt U, Tiselius P, Titelman J, Aksnes DL. Evidence of diel vertical migration in Mnemiopsis leidyi. PLoS One. 2014;9(1):e86595. https://doi.org/10.1371/journal.pone.0086595.
Tian R, Geng Y, Guo H, Yang C, Seim I, Yang G. Comparative analysis of the superoxide dismutase gene family in Cetartiodactyla. J Evol Biol. 2021;34:1046–60. https://doi.org/10.1111/jeb.13792.
Cantú-Medellín N, Byrd B, Hohn A, Vázquez-Medina JP, Zenteno- Savín T. Differential antioxidant protection in tissues from marine mammals with distinct diving capacities. Shallow/ short vs. deep/long divers. Comp Biochem Physiol. 2011;158:438–43. https://doi.org/10.1016/j.cbpa.2010.11.029.
Möglich A, Ayers RA, Moffat K. Structure and signaling mechanism of Per-ARNT-Sim domains. Structure. 2009;17:1282–94. https://doi.org/10.1016/j.str.2009.08.011.
Henry JT, Crosson S. Ligand-binding PAS domains in a genomic, cellular, and structural context. Annu Rev Microbiol. 2011;65:261–86. https://doi.org/10.1146/annurev-micro-121809-151631.
Vogt JH, Schippers JH. Setting the PAS, the role of circadian PAS domain proteins during environmental adaptation in plants. Front Plant Sci. 2015;6:513. https://doi.org/10.3389/fpls.2015.00513.
van der Lee R, Buljan M, Lang B, Weatheritt RJ, Daughdrill GW, Dunker AK, Fuxreiter M, Gough J, Gsponer J, Jones DT, Kim PM, Kriwacki RW, Oldfield CJ, Pappu RV, Tompa P, Uversky VN, Wright PE, Babu MM. Classification of intrinsically disordered regions and proteins. Chem Rev. 2014;114(13):6589–631. https://doi.org/10.1021/cr400525m.
Uversky VN. Intrinsically disordered proteins and their “mysterious” (Meta)physics. Front Phys. 2019;7:10. https://doi.org/10.3389/fphy.2019.00010.
Babu MM. The contribution of intrinsically disordered regions to protein function, cellular complexity, and human disease. Biochem Soc Trans. 2016;44:1185–200. https://doi.org/10.1042/BST20160172.
Lu P, Ma D, Yan C, Gong X, Du M, Shi Y. Structure and mechanism of a eukaryotic transmembrane ascorbate-dependent oxidoreductase. Proc Natl Acad Sci U S A. 2014;111(5):1813–8. https://doi.org/10.1073/pnas.1323931111.
McKie AT, Barrow D, Latunde-Dada GO, Rolfs A, Sager G, Mudaly E, Mudaly M, Richardson C, Barlow D, Bomford A, Peters TJ, Raja KB, Shirali S, Hediger MA, Farzaneh F, Simpson RJ. An iron-regulated ferric reductase associated with the absorption of dietary iron. Science. 2001;291:1755–9.
Asard H, Barbaro R, Trost P, Bérczi A. Cytochromes b561: ascorbate-mediated trans-membrane electron transport. Antioxid Redox Signal. 2013;19(9):1026–35. https://doi.org/10.1089/ars.2012.5065.
Herve C, Tonon T, Collen J, Corre E, Boyen C. NADPH oxidases in Eukaryotes: red algae provide new hints! Curr Genet. 2006;49:190–204.
Aguirre J, Lambeth JD. Nox enzymes from fungus to fly to fish and what they tell us about Nox function in mammals. Free Radic Biol Med. 2010;49:1342–53.
Eddy SR. Accelerated profile HMM searches. PLOS Comp Biol. 2011;7:e1002195.
Thumuluri V, Almagro Armenteros JJ, Johansen AR, Nielsen H, Winther O. DeepLoc 2.0: multi-label subcellular localization prediction using protein language models. Nucleic Acids Res. 2022;30;50(W1):W228–34. https://doi.org/10.1093/nar/gkac278.
Almagro Armenteros JJ, Salvatore M, Winther O, Emanuelsson O, von Heijne G, Elofsson A, Nielsen H. Detecting sequence signals in targeting peptides using deep learning. Life Sci Alliance. 2019;5:e201900429. https://doi.org/10.26508/lsa.201900429.
Biberstine-Kinkade KJ, DeLeo FR, Epstein RI, LeRoy BA, Nauseef WM, Dinauer MC. Heme-ligating histidines in flavocytochrome b(558): identification of specific histidines in gp91(phox). J Biol Chem. 2001;276(33):31105–12. https://doi.org/10.1074/jbc.M103327200.
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nuc Acid Res. 2002;30(14):3059–66. https://doi.org/10.1093/nar/gkf436.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. https://doi.org/10.1093/molbev/mst010.
Larsson A. AliView: a fast and lightweight alignment viewer and editor for large data sets. Bioinformatics. 2014;30(22):3276–8. https://doi.org/10.1093/bioinformatics/btu531.
Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol Biol Evol. 2015;32:268–74. https://doi.org/10.1093/molbev/msu300.
Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35:518–22. https://doi.org/10.1093/molbev/msx281.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9. https://doi.org/10.1038/nmeth.4285.
Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6. https://doi.org/10.1093/nar/gkab301.
The authors thank Nick Rhodes from the Queensland Cyber Infrastructure Foundation (QCIF), Haojing Shao and Enya Hermsen for their advice on bioinformatical analyses and computational assistance; Alastair McEwan and Mathias Jönsson for their editorial assistance on early drafts of the manuscript; and all members of the Degnan lab for their valuable discussions.
This research was supported by an Australian Research Council grant awarded to SM Degnan (DP190102521) and by a University of Queensland Research Training Program scholarship awarded to OH Hewitt.
Ethics approval and consent to participate
Consent for publication
The authors declare there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Details and availability of the 19 metazoan genomes used in this study. Table S2. Details and availability of non-metazoan genomes used in this study. Fig S1. Maximum likelihood phylogenetic tree of CuZnSOD enzyme family. Fig S2. Maximum likelihood phylogenetic tree of MnSOD enzyme family.
Table displaying total number of sequences identified within the genomes of each metazoan species based on the presence of characteristic domains for NOX, CuZnSOD, and MnSOD proteins.
Output from DeepLoc-2.0 analysis predicting subcellular localisation for identified CuZnSOD and MnSOD proteins.
Output from TargetP-2.0 analysis predicting presence of signal peptide target region.
All alignments and original protein sequences used to generate phylogenetic trees.
About this article
Cite this article
Hewitt, O.H., Degnan, S.M. Distribution and diversity of ROS-generating enzymes across the animal kingdom, with a focus on sponges (Porifera). BMC Biol 20, 212 (2022). https://doi.org/10.1186/s12915-022-01414-z
- Reactive oxygen species (ROS)