Skip to main content

Computational identification of a new SelD-like family that may participate in sulfur metabolism in hyperthermophilic sulfur-reducing archaea



Selenium (Se) and sulfur (S) are closely related elements that exhibit similar chemical properties. Some genes related to S metabolism are also involved in Se utilization in many organisms. However, the evolutionary relationship between the two utilization traits is unclear.


In this study, we conducted a comparative analysis of the selenophosphate synthetase (SelD) family, a key protein for all known Se utilization traits, in all sequenced archaea. Our search showed a very limited distribution of SelD and Se utilization in this kingdom. Interestingly, a SelD-like protein was detected in two orders of Crenarchaeota: Sulfolobales and Thermoproteales. Sequence and phylogenetic analyses revealed that SelD-like protein contains the same domain and conserved functional residues as those of SelD, and might be involved in S metabolism in these S-reducing organisms. Further genome-wide analysis of patterns of gene occurrence in different thermoproteales suggested that several genes, including SirA-like, Prx-like and adenylylsulfate reductase, were strongly related to SelD-like gene. Based on these findings, we proposed a simple model wherein SelD-like may play an important role in the biosynthesis of certain thiophosphate compound.


Our data suggest novel genes involved in S metabolism in hyperthermophilic S-reducing archaea, and may provide a new window for understanding the complex relationship between Se and S metabolism in archaea.


Selenium (Se) is an important micronutrient in all three domains of life [1]. The major biological form of this trace element is selenocysteine (Sec), which is a Se analogue of cysteine (Cys) and incorporated into selenoproteins at an in-frame UGA codon [24]. All Sec-utilizing organisms use an elaborate synthetic and translational machinery for Sec biosynthesis and its insertion into a selenoprotein, which is different from the enzymatic system for the 20 canonical amino acids. Both common key components and differences exist among Sec-decoding mechanisms employed by archaea, bacteria and eukaryotes [38]. The mechanism of Sec biosynthesis in bacteria is best described in Escherichia coli[6, 810]. Several proteins and a specific tRNA are known to be involved in this process, including bacterial Sec synthase (SelA), Sec-specific elongation factor (SelB), Sec-tRNA[Ser]Sec (SelC) and selenophosphate synthetase (SelD). In archaea and eukaryotes, substitute/additional components are needed for Sec biosynthesis, such as O-phosphoseryl-tRNA[Ser]Sec kinase (PSTK), archaeal/eukaryotic Sec synthase (SecS), SECIS-binding protein 2 (SBP2) and ribosomal protein L30 [7, 11, 12]. In addition, Se is also found in 5-methylaminomethyl-2-selenouridine (SeU) which is located at the wobble position of the anticodons of some tRNAs, as well as in a Se-containing cofactor bound to certain moly-bdenum-containing hydroxylases [1316]. Specific genes are known to be involved in each of these processes (SelD/YbbB for SeU trait [17] and SelD/YqeB/YqeC for Se-containing cofactor [16]).

Among known Se-related genes, SelD is considered as a key gene which is essential for all known Se utilization traits. It is a member of a small ATP-binding superfamily that also includes the purine biosynthetic enzymes PurM (phosphoribosylaminoimidazole synthetase) and PurL (formylglycinamide ribonucleotide amidotransferase), thiamine-monophosphate kinase (ThiL) and nickel-iron hydrogenase maturation protein HypE [18]. The function of SelD is to catalyze the activation of selenide with adenosine 5′-triphosphate (ATP) to generate selenophosphate, the essential Se donor for the formation of Sec and SeU residues in proteins and RNAs, respectively [19]. Several comparative genomics studies have shown that SelD gene could be used as a general signature for Se utilization in biology [16, 17, 20]. Since Sec is the main form of Se utilization, the majority of studies on SelD focused on its function and regulation during the process of Sec biosynthesis in a variety of model organisms [2123]. In contrast, the evolutionary and functional relationships between SelD and other homologs are not fully understood.

On the other hand, Se and sulfur (S) are two closely related elements utilized for a vast array of biochemical reactions [24]. It has been known that some genes might be involved in the metabolism of both elements. For example, sulfate transporters are responsible for the uptake of selenate in bacteria [25]. It is also known that Sec lyases (SCLs) and Cys desulfurases (CDs) catalyze the removal of Se or S from Sec or Cys and generally act on both substrates, in spite that SCLs from higher eukaryotes (including human) are specific for Sec and do not accept Cys as substrate [24]. Thus, it would be interesting to investigate the relationship between Se and S metabolic pathways in different organisms.

In this study, we carried out a comparative and phylogenetic analysis of SelD family in archaea based on the rapidly increased genomes sequenced in this kingdom. Occurrence of SelDs and their homologs were identified, which revealed a limited utilization of Se in archaea. In addition, we found a new SelD-like family in several hyperthermophilic S-utilizing archaea, which contains the same catalytic Cys residue as those in SelD proteins. Phylogenetic and genome context analyses of SelD-like genes revealed that they are more likely to share a common ancestor with SelD genes and may participate in S metabolism, probably involved in the biosynthesis of certain thiophosphate compound in these organisms. Further bioinformatics approaches were used to analyze patterns of gene occurrence, and several genes that may co-occur with SelD-like gene were identified, implying functional links among them. Thus, our results provide new insights into the complex and dynamic relationship between Se and S metabolism in S-utilizing archaea.

Results and discussion

Distribution of SelD and Se utilization traits in archaea

In the recent decade, comparative genomics has become an important strategy for understanding the evolutionary dynamics of a variety of cellular pathways in different organisms. Several comparative studies on Se have been carried out to examine Se utilization traits and identify new genes involved in Se metabolism in prokaryotes and eukaryotes [16, 2628]. For example, based on a comparative genomics analysis of SelD genes in bacteria, we successfully identified a new Se utilization trait as well as several genes involved in this pathway [16].

In this study, we analyzed the distribution of both SelD gene and the three known Se utilization pathways in different archaeal taxa (Figure 1; details are shown in Additional file 1: Table S1). Among all sequenced organisms, only 26 (12.1%) contain SelD genes. This is consistent with our previous observation that Se utilization is quite restricted in this kingdom of life [28]. Further investigation of the three known Se utilization traits revealed that 15 Sec-utilizing (SelD + SelA + SelB), 14 SeU-utilizing (SelD + YbbB) and 10 Se-cofactor-utilizing (SelD + YqeB + YqeC) organisms were present in different clades. The SeU and Se-containing cofactor utilization traits were only found in Methanococcales and Halobacteriales, respectively. Only methanococcales have both Sec-decoding and SeU-utilizing traits.

Figure 1
figure 1

Overview of distribution of SelD proteins and Se utilization traits in archaea. The tree is based on 16S rRNA sequences, which was simplified to only show major branches. The numbers of organisms that contain SelD and known Se utilization traits are shown. Distribution of organisms containing SelD-like family is highlighted in red.

Surprisingly, a newly sequenced archaeon, Methanoplanus petrolearius (belonging to the order Euryarchaeota/Methanomicrobiales), was found to have an orphan SelD gene and lack components of any known Se utilization trait. The SelD protein in this organism is a typical SelD which is highly similar to other orthologs (say, 51% sequence identity when compared with Methanocaldococcus jannaschii SelD using BLAST) and contains a conserved Cys residue in the active site. In a separate study, orphan SelD genes were also found in some newly sequenced bacteria, implying an additional but unknown SelD-dependent use of Se (unpublished data).

Identification of a new SelD-like family in crenarchaeota

Several distant homologs of SelD (named SelD-like thereafter) were exclusively found in two orders of the phylum Crenarchaeota: Sulfolobales and Thermoproteales, in which all sequenced organisms do not have any known Se utilization trait (Figure 1 and Additional file 1: Table S1). The sequence similarities between SelD and SelD-like proteins are low (for example, only 26% identity was observed between M. jannaschii SelD and Sulfolobus acidocaldarius SelD-like). However, these SelD-like proteins contain the same domain as SelD (COG0709) and the conserved Cys residue corresponding to the catalytic Sec/Cys in the active site of SelD, both of which are not detected in other families of the ATP-binding superfamily that contains SelD (Figure 2 and Additional file 2: Figure S1; SelD-like and SelD sequences are available in Additional file 3). Further analysis of the occurrence of SelD-like protein in bacteria and eukaryotes revealed that it is absent in all sequenced organisms in the two kingdoms.

Figure 2
figure 2

Multiple sequence alignments of SelD and SelD-like proteins. The alignments show the N terminal regions that contain the active site in representative SelD and SelD-like proteins. Conserved residues are highlighted. Sec (U) and the corresponding Cys (C) residues are shown in red and blue, respectively. The conserved Cys residue corresponding to the catalytic Sec/Cys in the active site of SelD could not be detected in other families of the ATP-binding superfamily that contains SelD.

Recently, the crystal structures of SelD have been solved in several bacteria, which provide important basis for understanding the catalytic mechanism of SelD-mediated selenophosphate synthesis. For example, based on the crystal structures of Aquifex aeolicus SelD and its complex with ATP analogue, Itoh et al. found that four aspartic acid (Asp) residues may coordinate metal ions (cobalt or magnesium) to bind the phosphate groups of ATP [29]. In addition, the Sec/Cys-X-X-Lys motif is also essential for SelD activity. Very recently, the crystal structure of a C17S mutant of E. coli SelD also confirmed that these Asp residues which are involved in the binding of magnesium at the active site to allow the interaction with ATP, and a conserved asparagine (Asn87 in E. coli SelD), play an essential role for catalysis in E. coli SelD [30]. Interestingly, all of these important residues are present in all SelD-like proteins (at either the same or very close position based on sequence alignment analysis, Additional file 2: Figure S1). Thus, the SelD-like family may also bind certain metal (e.g., magnesium) and use a similar mechanism for the catalytic reaction of its substrate (most likely ATP or its analogue).Considering that SelD belongs to an ATP-binding superfamily which also includes PurM, PurL, ThiL and HypE families, we collected proteins from representative organisms for each of them. Phylogenetic analysis revealed that SelD-like proteins are different from SelD family and form a separate branch of this superfamily (Figure 3). The evolutionary relationship between SelD and SelD-like proteins appeared to be closer than other families, implying that SelD-like and SelD proteins have a common ancestor. Considering that all sequenced crenarchaeota lack known Se utilization trait, we assume that these SelD-like proteins may be unrelated to Se metabolism although their function is unclear so far. However, the possibility that SelD-like protein may participate in an unknown Se-related process in these organisms could not be completely excluded.

Figure 3
figure 3

Phylogenetic analysis of the ATP-binding superfamily that contains SelD and other families. The tree is built based on representative sequences of each family. Bootstrapping values are included. The measurement of distance for the branch length is indicated.

Functional analysis of SelD-like proteins

As mentioned above, genes encoding SelD-like proteins were only detected in Sulfolobales and Thermoproteales; however, Desulfurococcales, the third order of Crenarchaeota, lack this gene. This may suggest that SelD-like protein evolved in the ancestor of crenarchaeota and lost in desulfurococcales. To further study the possible function of SelD-like proteins, we examined genes that are located within 10,000 nucleotides upstream and downstream of SelD-like gene in all organisms containing this gene. It has been known for some time that many clustered genes, especially those organized in one operon, may be functionally related [31]. Thus, it is possible to identify functional partners of SelD-like gene based on gene neighborhood analysis. The genomic context of SelD-like gene in representative genomes is shown in Figure 4. A comprehensive analysis of genomic context of SelD-like gene in all completely sequenced sulfolobales and thermoproteales genomes is shown in Additional file 2: Figure S2.

Figure 4
figure 4

Genomic context of SelD-like gene in representative genomes. Functionally related genes in the genomes of Sulfolobus islandicus HVE10/4, Pyrobaculum arsenaticum DSM 13514 and Vulcanisaeta distributa DSM 14429 are color coded. Coding direction is also indicated. HP represents hypothetical proteins.

First, all SelD-like proteins in Thermoproteales were found to be fused with an acylphosphatase-like domain. Acylphosphatase belongs to a large protein family (COG1254) that specifically catalyzes the hydrolysis of the carboxyl-phosphate bond of acylphosphate substrates such as carbamoylphosphate and 1,3-biphosphoglycerate, and results in the production of carboxylate and phosphate [32]. It is widely distributed in all three domains of life (including all sulfolobales) and has been found to be involved in a number of physiological processes such as the regulation of the glycolytic pathway and pyrimidine biosynthesis, though its biological function is as yet unknown [33]. Several recent structure-based studies of acylphosphatase from different hyperthermophilic archaea such as S. solfataricus and Pyrococcus horikoshii have revealed that normally folded acylphosphatases retain a significant tendency to form amyloid fibrils through a direct assembly of monomers in their native-like topology [34, 35]. The fusion of SelD-like and acylphosphatase-like domains indicates that there might be a strong functional link between acylphosphatase-like and SelD-like proteins. It is possible that this acylphosphatase-like domain is involved in hydrolysis of certain acylphosphate-like substrate that is synthesized by SelD-like protein. It is also possible, but less likely, that acylphosphate-like protein may provide phosphate substrate for the proper function of SelD-like protein (instead of ATP).

A second interesting finding is that several common genes were found to be localized very close to SelD-like genes in both sulfolobales and thermoproteales, including genes encoding SirA-like (pfam01206) and peroxiredoxin(Prx)-like (COG2210) proteins. The SirA-like protein is considered as a S transfer protein that belongs to a predicted redox protein superfamily and has been found to be associated with SelD [16]. Homologs of SirA-like proteins are widespread in almost all prokaryotes and may have different functions. The Prx-like protein is a distant homolog of Prx family, which is a ubiquitous family of antioxidant enzymes and shares the same basic catalytic mechanism, in which a redox-active Cys in the active site is oxidized to a sulfenic acid by the peroxide substrate [36]. Our gene neighborhood analysis suggested that both redox proteins may participate in the cellular processes involving SelD-like proteins in archaea.

Although there is no clear evidence with regard to the function of SelD-like, we could always find that genes associated with S metabolic pathways are located nearby. Besides SirA-like and Prx-like, additional S-related genes were also observed in almost all sequenced crenarchaeota containing SelD-like, such as rhodanese-related sulfurtransferase PspE, dissimilatory sulfite reductase (DSR) and Fe-S oxidoreductases (Figure 4 and Additional file 2: Figure S2). Thus, it is possible that SelD-like protein is involved in S metabolism. Considering that SelD catalyzes the biosynthesis of selenophosphate from selenide and ATP, plus that SelD-like proteins have the same functional residues as SelD, we finally hypothesized that SelD-like protein might be involved in biosynthesis of certain thiophosphate molecule using S-containing molecule (sulfide?) and ATP/ATP analogue. The product of SelD-like might be similar to selenophosphate where the Se atom is replaced with S to form a P=S double bond (named thiophosphate (P=S) thereafter). This is different from many other thiophosphate compounds containing S-P single bond(s). However, the restricted distribution of SelD-like gene implied that it might be only needed for S utilization in sulfolobales and thermoproteales, all of which are hyperthermophilic S-reducing archaea. Instead, more and more desulfurococcales species are known to be unable to use S or S-containing compounds for energy formation, which may probably result in the loss of this gene in the Desulfurococcales order of Crenarchaeota. Further experimental studies are needed to validate the biological function of SelD-like protein and the presence and/or function of such thiophosphate compound.

Computational search for additional genes associated with SelD-like gene

SelD-like gene could be detected in all sulfolobales and almost all thermoproteales. However, it was absent in three completely sequenced thermoproteales including Pyrobaculum islandicum, P. neutrophilum and Thermofilum pendens, implying a very recent loss of this gene in these organisms. Based on the comparison of genomic sequences of both SelD-like-containing and SelD-like-lacking thermoproteales, it is possible to identify additional genes associated with SelD-like.

We analyzed the presence of all the genes annotated in the P. arsenaticum genome in other sequenced thermoproteales, and focused on those genes that occur exclusively in organisms containing SelD-like gene and are absent in organisms that lack this gene. Top candidate genes that co-occur with SelD-like are shown in Table 1, including SirA-like and Prx-like genes mentioned above.

Table 1 P. arsenaticum genes that are present exclusively in all SelD-like-containing thermoproteales

It is interesting that the adenylylsulfate reductase, or called adenosine 5′-phosphosulfate (APS) reductase, was found to be strongly related to SelD-like gene. This enzyme contains two subunits (AprA and AprB), which participates in S metabolism and catalyzes the reduction of APS to sulfite and AMP [37]. However, both genes were absent in sequenced sulfolobales genomes, implying the lack of such a relationship in this archaeal order that may use an unknown enzyme instead of AprA/B during S metabolism. Further analysis of the occurrence of AprA and AprB genes in other prokaryotes revealed that they are still highly conserved in many other sulfate-reducing and S-oxidizing prokaryotes (data not shown).

Another related gene is a predicted transporter that contains three transmembrane regions. It belongs to a large and diverse group of transporters (pfam07690, major facilitator superfamily), which includes uniporters, symporters, and antiporters of a variety of substrates such as ions, sugar phosphates, drugs, nucleosides, amino acids, and peptides. This gene is found in all sulfolobales and many other prokaryotes. However, its function is unclear. Further studies are needed to verify its biological role and the relationship with SelD-like gene.

A possible model for the involvement of SelD-like protein in S metabolism in Crenarchaeota

It has been reported that many S-reducing archaea carries out sulfate reduction via the pathway originally proposed for bacterial species [38, 39]. After sulfate transport across the cytoplasmic membrane, it is first activated with ATP to form APS and pyrophosphate, which is catalyzed by ATP sulfurylase. Then APS reductase catalyzes the reduction of APS to sulfite and AMP. After that, the enzyme DSR catalyzes the six-electron reduction of sulfite to sulfide (Figure 5). The natural electron donors for APS reductase and DSR are not known.Based on the results described above, the SelD-like protein may play a key role in the biosynthesis of thiophosphate (P=S) compound in sulfolobales and thermoproteales. In this study, two possible models were suggested (Figure 5). In Model 1, the substrates include sulfide and certain acylphosphate, and the product is considered as acyl-thiophosphate (P=S), which may be further hydrolyzed by the acylphosphatase-like protein (fused with SelD-like in thermoproteales). In Model 2 whose reaction is more similar to that for SelD, the substrates include sulfide and ATP, whereas the products are thiophosphate (P=S) and pyrophosphate (PPi); however, the functional link between SelD-like and acylphosphatase-like proteins is not clear in this model.

Figure 5
figure 5

A predicted model for the involvement of SelD-like protein in the S metabolism in thermoproteales. Known S-related genes are shown in black, whereas SelD-like and acylphosphatase-like proteins are highlighted in red and blue, respectively.

An interesting hypothesis is that SelD-like protein might be involved in cysteinyl-tRNACys biosynthesis based on a pathway similar to Sec biosynthesis. It has been known that cysteinyl-tRNACys is typically synthesized by cysteinyl-tRNA synthetase (CysS). In methanogenic archaea that lack the CysS gene, cysteinyl-tRNACys is synthesized by O-phosphoseryl-tRNA synthetase (SepS) that acylates tRNACys with phosphoserine (Sep), and Sep-tRNA: Cys-tRNA synthase (SepCysS) that converts Sep-tRNACys to Cys-tRNACys[40]. One possibility is that the thiophosphate (P=S) compound produced by SelD-like protein may be a special S donor for the tRNA-dependent Cys formation in these organisms. To test this idea, we analyzed the occurrence of CysS, SepS and SepCysS genes in SelD-like-containing crenarchaeota. All of them only have orthologous CysS genes. The lack of SepS and SepCysS genes in these organisms suggest that they do not have a known tRNA-dependent Cys formation pathway. It is possible that the special product of SelD-like protein is only present in sulfolobales and thermoproteales. However, additional genes and specific mechanisms related to S metabolism may have evolved in other S-reducing thermophiles. It is also possible that the occurrence of SelD-like gene in hyperthermophilic S-reducing archaea may just reflect a fundamental correlation, which is unrelated to S metabolism. In the future, experimental evidence is required to test these hypotheses.


In this study, we carried out comparative and phylogenetic studies of SelD family in all sequenced archaea. Our data revealed restricted occurrence of SelD as well as limited utilization of Se in this kingdom. In addition, a distant homolog of SelD (SelD-like) was found to form a separate family in hyperthermophilic S-utilizing Crenarchaeota: Sulfolobales and Thermoproteales. The SelD-like protein has the same catalytic Cys as well as other functional residues as those in SelD proteins. Further sequence and genome context analysis of SelD-like genes revealed that they might be involved in S metabolism in these organisms. We also analyzed patterns of gene occurrence in different thermoproteales at a genome-wide level, and identified several genes that may co-occur with SelD-like gene. Finally, a model for the possible function of SelD-like protein in these organisms was proposed, which may provide insights into understanding the relationship between Se and S utilization in biology.


Genomic databases and sequences

Sequenced archaeal genomes from the Entrez Microbial Genome Database at NCBI were used in this study [41]. A total of 215 organisms were analyzed (as of August 2013). The M. jannaschii SelD protein (gi: 190352215) was used as query sequence to search for SelD homologs in archaea. TBLASTN program was used to identify genes coding for homologs with the cutoff of E-value ≤0.1. In addition, other key components of different Se utilization traits (such as archaeal SecS, SelB, YbbB, YqeB and YqeC) from representative organisms were also used for identification of the occurrence of each Se trait in archaeal species. Orthologous proteins were defined using the conserved domain database (COG, Pfam, TIGR, etc.) and bidirectional best hits [42].

Multiple sequence alignment and phylogenetic tree reconstruction

Multiple sequence alignments were performed by Clustal Omega program (Version 1.2.0) [43] with number of iterations set to 2. The resulting multiple alignments were then checked for conservation of functional residues and manually edited. Phylogenetic trees were reconstructed by PHYLIP programs [44]. The PROML program was used for producing the phylogenetic trees by the maximum likelihood method with the Jones-Taylor-Thornton (JTT) probability model and the global rearrangements option was used. Bootstrap with 100 repetitions was carried out by SEQBOOT program to assess the confidence degree of nodes in the phylogenetic trees. CONSENSE program was used for computing the consensus tree. The neighbor-joining (NJ) method was used for the phylogenetic reconstruction to validate the tree topology with NEIGHBOR program. Robustness of these phylogenies was evaluated by two additional programs, PHYML (maximum likelihood analysis) [45] and MrBayes (Bayesian estimation of phylogeny) [46] with different parameters (for example, transition matrices and evolutionary models).

Identification of genes associated with SelD-like gene

Initial searches for the occurrence of SelD-like genes revealed that only the order Thermoproteales (in the archaeal phylum Crenarchaeota) has both SelD-like-containing and SelD-like-lacking organisms (9 and 3 organisms, respectively). To identify genes associated with SelD-like, we adopted a modified strategy which had been used previously to identify the archaeal/eukaryotic SecS gene [47]. Protein sequences of all annotated genes in P. arsenaticum (a representative thermoproteales containing SelD-like gene) were used as a query dataset. BLAST programs were used to search these sequences against all sequenced thermoproteales genomes using the following criteria: (i) e-value is less than 1e-07; (ii) the length of the alignment should be extending over half the length of the query sequence; and (iii) the alignment identity should be at least 30%. A simple Perl script was developed to parse the BLAST output and examine presence/absence of homologs in all analyzed genomes. Genes present in any of the thermoproteales lacking SelD-like were dismissed. The remaining genes were then searched against thermoproteales containing SelD-like to determine their occurrences. Top candidate genes were further analyzed for possible function. The occurrence of these candidates in SelD-like-containing sulfolobales was also analyzed.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files.









Sec synthase


Sec-specific elongation factor




Selenophosphate synthetase


O-phosphoseryl-tRNA[Ser]Sec kinase


Archaeal/eukaryotic Sec synthase


SECIS-binding protein 2




Phosphoribosylaminoimidazole synthetase


Formylglycinamide ribonucleotide amidotransferase


Thiamine-monophosphate kinase


Nickel-iron hydrogenase maturation protein


Adenosine 5′-triphosphate




Sec lyase


Cys desulfurases


Aspartic acid






Dissimilatory sulfite reductase


Adenosine 5′-phosphosulfate




Cysteinyl-tRNA synthetase


O-phosphoseryl-tRNA synthetase




Sep-tRNA:Cys-tRNA synthase.


  1. Rayman MP: Selenium and human health. Lancet. 2012, 379: 1256-1268. 10.1016/S0140-6736(11)61452-9.

    Article  CAS  PubMed  Google Scholar 

  2. Stadtman TC: Selenocysteine. Annu Rev Biochem. 1996, 65: 83-100. 10.1146/

    Article  CAS  PubMed  Google Scholar 

  3. Hatfield DL, Gladyshev VN: How selenium has altered our understanding of the genetic code. Mol Cell Biol. 2002, 22: 3565-3576. 10.1128/MCB.22.11.3565-3576.2002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Low SC, Berry MJ: Knowing when not to stop: selenocysteine incorporation in eukaryotes. Trends Biochem Sci. 1996, 21: 203-208. 10.1016/0968-0004(96)10025-6.

    Article  CAS  PubMed  Google Scholar 

  5. Driscoll DM, Copeland PR: Mechanism and regulation of selenoprotein synthesis. Annu Rev Nutr. 2003, 23: 17-40. 10.1146/annurev.nutr.23.011702.073318.

    Article  CAS  PubMed  Google Scholar 

  6. Böck A: Biosynthesis of selenoproteins: an overview. Biofactors. 2000, 11: 77-78. 10.1002/biof.5520110122.

    Article  PubMed  Google Scholar 

  7. Allmang C, Wurth L, Krol A: The selenium to selenoprotein pathway in eukaryotes: more molecular partners than anticipated. Biochim Biophys Acta. 2009, 1790: 1415-1423. 10.1016/j.bbagen.2009.03.003.

    Article  CAS  PubMed  Google Scholar 

  8. Thanbichler M, Böck A: Selenoprotein biosynthesis: purification and assay of components involved in selenocysteine biosynthesis and insertion in Escherichia coli. Methods Enzymol. 2002, 347: 3-16.

    Article  CAS  PubMed  Google Scholar 

  9. Turner RJ, Weiner JH, Taylor DE: Selenium metabolism in Escherichia coli. Biometals. 1998, 11: 223-227. 10.1023/A:1009290213301.

    Article  CAS  PubMed  Google Scholar 

  10. Donovan J, Copeland PR: Threading the needle: getting selenocysteine into proteins. Antioxid Redox Signal. 2010, 12: 881-892. 10.1089/ars.2009.2878.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Xu XM, Carlson BA, Zhang Y, Mix H, Kryukov GV, Glass RS, Berry MJ, Gladyshev VN, Hatfield DL: New developments in selenium biochemistry: selenocysteine biosynthesis in eukaryotes and archaea. Biol Trace Elem Res. 2007, 119: 234-241. 10.1007/s12011-007-8003-9.

    Article  CAS  PubMed  Google Scholar 

  12. Squires JE, Berry MJ: Eukaryotic selenoprotein synthesis: mechanistic insight incorporating new factors and new functions for old factors. IUBMB Life. 2008, 60: 232-235. 10.1002/iub.38.

    Article  CAS  PubMed  Google Scholar 

  13. Ching WM, Alzner-DeWeerd B, Stadtman TC: A selenium-containing nucleoside at the first position of the anticodon in seleno-tRNAGlu from Clostridium sticklandii. Proc Natl Acad Sci U S A. 1985, 82: 347-350. 10.1073/pnas.82.2.347.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Haft DH, Self WT: Orphan SelD proteins and selenium-dependent molybdenum hydroxylases. Biol Direct. 2008, 3: 4-10.1186/1745-6150-3-4.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Srivastava M, Mallard C, Barke T, Hancock LE, Self WT: A selenium-dependent xanthine dehydrogenase triggers biofilm proliferation in Enterococcus faecalis through oxidant production. J Bacteriol. 2011, 193: 1643-1652. 10.1128/JB.01063-10.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Zhang Y, Turanov AA, Hatfield DL, Gladyshev VN: In silico identification of genes involved in selenium metabolism: evidence for a third selenium utilization trait. BMC Genomics. 2008, 9: 251-10.1186/1471-2164-9-251.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Romero H, Zhang Y, Gladyshev VN, Salinas G: Evolution of selenium utilization traits. Genome Biol. 2005, 6: R66-10.1186/gb-2005-6-8-r66.

    Article  PubMed Central  PubMed  Google Scholar 

  18. McCulloch KM, Kinsland C, Begley TP, Ealick SE: Structural studies of thiamin monophosphate kinase in complex with substrates and products. Biochemistry. 2008, 47: 3810-3821. 10.1021/bi800041h.

    Article  CAS  PubMed  Google Scholar 

  19. Veres Z, Kim IY, Scholz TD, Stadtman TC: Selenophosphate synthetase: enzyme properties and catalytic reaction. J Biol Chem. 1994, 269: 10597-10603.

    CAS  PubMed  Google Scholar 

  20. Couture RM, Sekowska A, Fang G, Danchin A: Linking selenium biogeochemistry to the sulfur-dependent biological detoxification of arsenic. Environ Microbiol. 2012, 14: 1612-1623. 10.1111/j.1462-2920.2012.02758.x.

    Article  CAS  PubMed  Google Scholar 

  21. Preabrazhenskaya YV, Kim IY, Stadtman TC: Binding of ATP and its derivatives to selenophosphate synthetase from Escherichia coli. Biochemistry (Mosc). 2009, 74: 910-916. 10.1134/S0006297909080136.

    Article  CAS  Google Scholar 

  22. Sculaccio SA, Rodrigues EM, Cordeiro AT, Magalhães A, Braga AL, Alberto EE, Thiemann OH: Selenocysteine incorporation in Kinetoplastid: selenophosphate synthetase (SELD) from Leishmania major and Trypanosoma brucei. Mol Biochem Parasitol. 2008, 162: 165-171. 10.1016/j.molbiopara.2008.08.009.

    Article  CAS  PubMed  Google Scholar 

  23. Xu XM, Carlson BA, Irons R, Mix H, Zhong N, Gladyshev VN, Hatfield DL: Selenophosphate synthetase 2 is essential for selenoprotein biosynthesis. Biochem J. 2007, 404: 115-120. 10.1042/BJ20070165.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Collins R, Johansson AL, Karlberg T, Markova N, van den Berg S, Olesen K, Hammarström M, Flores A, Schüler H, Schiavone LH, Brzezinski P, Arnér ES, Högbom M: Biochemical discrimination between selenium and sulfur 1: a single residue provides selenium specificity to human selenocysteine lyase. PLoS One. 2012, 7: e30581-10.1371/journal.pone.0030581.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Aguilar-Barajas E, Díaz-Pérez C, Ramírez-Díaz MI, Riveros-Rosas H, Cervantes C: Bacterial transport of sulfate, molybdate, and related oxyanions. Biometals. 2011, 24: 687-707. 10.1007/s10534-011-9421-x.

    Article  CAS  PubMed  Google Scholar 

  26. Lobanov AV, Fomenko DE, Zhang Y, Sengupta A, Hatfield DL, Gladyshev VN: Evolutionary dynamics of eukaryotic selenoproteomes: large selenoproteomes may associate with aquatic life and small with terrestrial life. Genome Biol. 2007, 8: R198-10.1186/gb-2007-8-9-r198.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Zhang Y, Romero H, Salinas G, Gladyshev VN: Dynamic evolution of selenocysteine utilization in bacteria: a balance between selenoprotein loss and evolution of selenocysteine from redox active cysteine residues. Genome Biol. 2006, 7: R94-10.1186/gb-2006-7-10-r94.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Zhang Y, Gladyshev VN: General trends in trace element utilization revealed by comparative genomic analyses of Co, Cu, Mo, Ni, and Se. J Biol Chem. 2010, 285: 3393-3405. 10.1074/jbc.M109.071746.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Itoh Y, Sekine S, Matsumoto E, Akasaka R, Takemoto C, Shirouzu M, Yokoyama S: Structure of selenophosphate synthetase essential for selenium incorporation into proteins and RNAs. J Mol Biol. 2009, 385 (5): 1456-1469. 10.1016/j.jmb.2008.08.042.

    Article  CAS  PubMed  Google Scholar 

  30. Noinaj N, Wattanasak R, Lee DY, Wally JL, Piszczek G, Chock PB, Stadtman TC, Buchanan SK: Structural insights into the catalytic mechanism of Escherichia coli selenophosphate synthetase. J Bacteriol. 2012, 194: 499-508. 10.1128/JB.06012-11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Osbourn AE, Field B: Operons. Cell Mol Life Sci. 2009, 66: 3755-3775. 10.1007/s00018-009-0114-3.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Mizuno Y, Kanesaka Y, Fujita H, Minowa O, Shiokawa H: The primary structure of two molecular species of porcine organ-common type acylphosphatase. J Biochem. 1991, 110: 790-794.

    CAS  PubMed  Google Scholar 

  33. Stefani M, Taddei N, Ramponi G: Insights into acylphosphatase structure and catalytic mechanism. Cell Mol Life Sci. 1997, 53: 141-151. 10.1007/PL00000585.

    Article  CAS  PubMed  Google Scholar 

  34. Bemporad F, Vannocci T, Varela L, Azuaga AI, Chiti F: A model for the aggregation of the acylphosphatase from Sulfolobus solfataricus in its native-like state. Biochim Biophys Acta. 2008, 1784: 1986-1996. 10.1016/j.bbapap.2008.08.021.

    Article  CAS  PubMed  Google Scholar 

  35. Bemporad F, Chiti F: “Native-like aggregation” of the acylphosphatase from Sulfolobus solfataricus and its biological implications. FEBS Lett. 2009, 583: 2630-2638. 10.1016/j.febslet.2009.07.013.

    Article  CAS  PubMed  Google Scholar 

  36. Flohé L, Budde H, Hofmann B: Peroxiredoxins in antioxidant defense and redox regulation. Biofactors. 2003, 19: 3-10. 10.1002/biof.5520190102.

    Article  PubMed  Google Scholar 

  37. Meyer B, Kuever J: Molecular analysis of the diversity of sulfate-reducing and sulfur-oxidizing prokaryotes in the environment, using aprA as functional marker gene. Appl Environ Microbiol. 2007, 73: 7664-7679. 10.1128/AEM.01272-07.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Dahl C, Trüper HG: Sulfite reductase and APS reductase from Archaeoglobus fulgidus. Methods Enzymol. 2001, 331: 427-441.

    Article  CAS  PubMed  Google Scholar 

  39. Pereira IA, Ramos AR, Grein F, Marques MC, da Silva SM, Venceslau SS: A comparative genomic analysis of energy metabolism in sulfate reducing bacteria and archaea. Front Microbiol. 2011, 2: 69-

    CAS  PubMed Central  PubMed  Google Scholar 

  40. Zhang CM, Liu C, Slater S, Hou YM: Aminoacylation of tRNA with phosphoserine for synthesis of cysteinyl-tRNA(Cys). Nat Struct Mol Biol. 2008, 15: 507-514. 10.1038/nsmb.1423.

    Article  CAS  PubMed  Google Scholar 

  41. The NCBI Genome Table. []

  42. Wolf YI, Koonin EV: A tight link between orthologs and bidirectional best hits in bacterial and archaeal genomes. Genome Biol Evol. 2012, 4: 1286-1294. 10.1093/gbe/evs100.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using clustal omega. Mol Syst Biol. 2011, 7: 539-

    Article  PubMed Central  PubMed  Google Scholar 

  44. Felsenstein J: PHYLIP – Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.

    Google Scholar 

  45. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.

    Article  PubMed  Google Scholar 

  46. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.

    Article  CAS  PubMed  Google Scholar 

  47. Xu XM, Carlson BA, Mix H, Zhang Y, Saira K, Glass RS, Berry MJ, Gladyshev VN, Hatfield DL: Biosynthesis of selenocysteine on its tRNA in eukaryotes. PLoS Biol. 2007, 5: e4-10.1371/journal.pbio.0050004.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


This work was supported by the National Natural Science Foundation of China under NO. 31171233, and a grant from Chinese Academy of Sciences (CAS) (2012OHTP10).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Yan Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YZ and GPL designed the study. GPL, YZ and LJ carried out computational studies, including sequence alignment, phylogenetic analysis and wrote the manuscript. QL and JZN edited the manuscript. All authors read and approved the final manuscript.

Gao-Peng Li, Liang Jiang contributed equally to this work.

Electronic supplementary material


Additional file 1: Table S1: This file contains Table S1, which shows the distribution of SelD, SelD-like and all known Se utilization traits in archaea. (XLS 48 KB)


Additional file 2: Figure S1: This file contains two supplementary figures: Figure S1. shows the multiple sequence alignment of SelD and SelD-like proteins. Figure S2. shows genomic context of SelD-like gene in completely sequenced genomes that contain this gene. (DOC 178 KB)


Additional file 3: This file contains representative SelD and SelD-like protein sequences used for multiple sequence alignment (Figure 2 and Additional file2: Figure S1).(TXT 6 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, GP., Jiang, L., Ni, JZ. et al. Computational identification of a new SelD-like family that may participate in sulfur metabolism in hyperthermophilic sulfur-reducing archaea. BMC Genomics 15, 908 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: