An algorithm for identification of bacterial selenocysteine insertion sequence elements and selenoprotein genes

General metabolism of Laribacter hongkongensis: a genome-wide analysis

BACKGROUND

Laribacter hongkongensis is associated with community-acquired gastroenteritis and traveler's diarrhea. In this study, we performed an in-depth annotation of the genes and pathways of the general metabolism of L. hongkongensis and correlated them with its phenotypic characteristics.

RESULTS

The L. hongkongensis genome possesses the pentose phosphate and gluconeogenesis pathways and tricarboxylic acid and glyoxylate cycles, but incomplete Embden-Meyerhof-Parnas and Entner-Doudoroff pathways, in agreement with its asaccharolytic phenotype. It contains enzymes for biosynthesis and β-oxidation of saturated fatty acids, biosynthesis of all 20 universal amino acids and selenocysteine, the latter not observed in Neisseria gonorrhoeae, Neisseria meningitidis and Chromobacterium violaceum. The genome contains a variety of dehydrogenases, enabling it to utilize different substrates as electron donors. It encodes three terminal cytochrome oxidases for respiration using oxygen as the electron acceptor under aerobic and microaerophilic conditions and four reductases for respiration with alternative electron acceptors under anaerobic conditions. The presence of complete tetrathionate reductase operon may confer survival advantage in mammalian host in association with diarrhea. The genome contains CDSs for incorporating sulfur and nitrogen by sulfate assimilation, ammonia assimilation and nitrate reduction. The existence of both glutamate dehydrogenase and glutamine synthetase/glutamate synthase pathways suggests an importance of ammonia metabolism in the living environments that it may encounter.

CONCLUSIONS

The L. hongkongensis genome possesses a variety of genes and pathways for carbohydrate, amino acid and lipid metabolism, respiratory chain and sulfur and nitrogen metabolism. These allow the bacterium to utilize various substrates for energy production and survive in different environmental niches.

Biological systems discovery in silico: radical S-adenosylmethionine protein families and their target peptides for posttranslational modification

Data mining methods in bioinformatics and comparative genomics commonly rely on working definitions of protein families from prior computation. Partial phylogenetic profiling (PPP), by contrast, optimizes family sizes during its searches for the cooccurring protein families that serve different roles in the same biological system. In a large-scale investigation of the incredibly diverse radical S-adenosylmethionine (SAM) enzyme superfamily, PPP aided in building a collection of 68 TIGRFAMs hidden Markov models (HMMs) that define nonoverlapping and functionally distinct subfamilies. Many identify radical SAM enzymes as molecular markers for multicomponent biological systems; HMMs defining their partner proteins also were constructed. Newly found systems include five groupings of protein families in which at least one marker is a radical SAM enzyme while another, encoded by an adjacent gene, is a short peptide predicted to be its substrate for posttranslational modification. The most prevalent, in over 125 genomes, featuring a peptide that we designate SCIFF (six cysteines in forty-five residues), is conserved throughout the class Clostridia, a distribution inconsistent with putative bacteriocin activity. A second novel system features a tandem pair of putative peptide-modifying radical SAM enzymes associated with a highly divergent family of peptides in which the only clearly conserved feature is a run of His-Xaa-Ser repeats. A third system pairs a radical SAM domain peptide maturase with selenocysteine-containing targets, suggesting a new biological role for selenium. These and several additional novel maturases that cooccur with predicted target peptides share a C-terminal additional 4Fe4S-binding domain with PqqE, the subtilosin A maturase AlbA, and the predicted mycofactocin and Nif11-class peptide maturases as well as with activators of anaerobic sulfatases and quinohemoprotein amine dehydrogenases. Radical SAM enzymes with this additional domain, as detected by TIGR04085, significantly outnumber lantibiotic synthases and cyclodehydratases combined in reference genomes while being highly enriched for members whose apparent targets are small peptides. Interpretation of comparative genomics evidence suggests unexpected (nonbacteriocin) roles for natural products from several of these systems.

Tandem use of selenocysteine: adaptation of a selenoprotein glutaredoxin for reduction of selenoprotein methionine sulfoxide reductase

Several engineered selenocysteine (Sec)-containing glutaredoxins (Grxs) and their enzymatic properties have been reported, but natural selenoprotein Grxs have not been previously characterized. We expressed a bacterial selenoprotein Grx from Clostridium sp. (also known as Alkaliphilus oremlandii) OhILAs in Escherichia coli and characterized this selenoenzyme and its natural Cys homologues in Clostridium and E. coli. The selenoprotein Grx had a 200-fold higher activity than its Sec-to-Cys mutant form, suggesting that Sec is essential for catalysis by this thiol-disulfide oxidoreductase. Kinetic analysis also showed that the selenoprotein Grx had a 10-fold lower K(m) than Cys homologues. Interestingly, this selenoenzyme efficiently reduced a Clostridium selenoprotein methionine sulfoxide reductase A (MsrA), suggesting that it is the natural reductant for the protein that is not reducible by thioredoxin, a common reductant for Cys-containing MsrAs. We also found that the selenoprotein Grx could not efficiently reduce a Cys version of Clostridium MsrA, whereas natural Clostridium and E. coli Cys-containing Grxs, which efficiently reduce Cys-containing MsrAs, poorly acted on the selenoprotein MsrA. This specificity for MsrA reduction could explain why Sec is utilized in Clostridium Grx and more generally provides a novel example of the use of Sec in biological systems.

Occurrence, genes and expression of the W/Se-containing class II benzoyl-coenzyme A reductases in anaerobic bacteria

Benzoyl-coenzyme A (CoA) reductases (BCRs) are key enzymes in the anaerobic degradation of aromatic compounds and catalyse the reductive dearomatization of benzoyl-CoA to cyclohexa-1,5-dienoyl-1-carboxyl-CoA. Class I BCRs are ATP-dependent FeS enzymes, whereas class II BCRs are supposed to be ATP-independent and contain W, FeS clusters, and most probably selenocysteine. The active site components of a putative eight subunit class II BCR, BamBCDEFGHI, were recently characterized in Geobacter metallireducens. In this organism bamB was identified as structural gene for the W-containing active site subunit; bamF was predicted to code for a selenocysteine containing electron transfer subunit. In this work the occurrence and expression of BCRs in a number of anaerobic, aromatic compound degrading model microorganisms was investigated with a focus on the BamB and BamF components. Benzoate-induced class II BCR in vitro activities were determined in the soluble protein fraction in all obligately anaerobic bacteria tested. Where applicable, the results were in agreement with Western blot analysis using BamB targeting antibodies. By establishing a specific bamB targeting PCR assay, bamB homologues were identified in all tested obligately anaerobic bacteria with the capacity to degrade aromatic compounds; a number of bamB sequences from Gram-negative/positive sulfate-reducing bacteria were newly sequenced. In several organisms at least two bamB paralogues per genome were identified; however, in nearly all cases only one of them was transcribed during growth on an aromatic substrate. These benzoate-induced bamB genes are proposed to code for the active site subunit of class II BCRs; the major part of them group into a phylogenetic subcluster within the bamB homologues. Results from in silico analysis suggested that all class II BCRs contain selenocysteine in the BamF, and in many cases also in the BamE subunit. The results obtained indicate that the distribution of the two classes of BCRs in anaerobic bacteria appears to be strictly ruled by the available free energy from the oxidation of the aromatic carbon source rather than by phylogenetic relationships.

Complete genome sequence of the rifamycin SV-producing Amycolatopsis mediterranei U32 revealed its genetic characteristics in phylogeny and metabolism

Amycolatopsis mediterranei is used for industry-scale production of rifamycin, which plays a vital role in antimycobacterial therapy. As the first sequenced genome of the genus Amycolatopsis, the chromosome of strain U32 comprising 10,236,715 base pairs, is one of the largest prokaryotic genomes ever sequenced so far. Unlike the linear topology found in streptomycetes, this chromosome is circular, particularly similar to that of Saccharopolyspora erythraea and Nocardia farcinica, representing their close relationship in phylogeny and taxonomy. Although the predicted 9,228 protein-coding genes in the A. mediterranei genome shared the greatest number of orthologs with those of S. erythraea, it was unexpectedly followed by Streptomyces coelicolor rather than N. farcinica, indicating the distinct metabolic characteristics evolved via adaptation to diverse ecological niches. Besides a core region analogous to that common in streptomycetes, a novel 'quasi-core' with typical core characteristics is defined within the non-core region, where 21 out of the total 26 gene clusters for secondary metabolite production are located. The rifamycin biosynthesis gene cluster located in the core encodes a cytochrome P450 enzyme essential for the conversion of rifamycin SV to B, revealed by comparing to the highly homologous cluster of the rifamycin B-producing strain S699 and further confirmed by genetic complementation. The genomic information of A. mediterranei demonstrates a metabolic network orchestrated not only for extensive utilization of various carbon sources and inorganic nitrogen compounds but also for effective funneling of metabolic intermediates into the secondary antibiotic synthesis process under the control of a seemingly complex regulatory mechanism.

Selenoprofiles: profile-based scanning of eukaryotic genome sequences for selenoprotein genes

Motivation: Selenoproteins are a group of proteins that contain selenocysteine (Sec), a rare amino acid inserted co-translationally into the protein chain. The Sec codon is UGA, which is normally a stop codon. In selenoproteins, UGA is recoded to Sec in presence of specific features on selenoprotein gene transcripts. Due to the dual role of the UGA codon, selenoprotein prediction and annotation are difficult tasks, and even known selenoproteins are often misannotated in genome databases.

Results: We present an homology-based in silico method to scan genomes for members of the known eukaryotic selenoprotein families: selenoprofiles. The core of the method is a set of manually curated highly reliable multiple sequence alignments of selenoprotein families, which are used as queries to scan genomic sequences. Results of the scan are processed through a number of steps, to produce highly accurate predictions of selenoprotein genes with little or no human intervention. Selenoprofiles is a valuable tool for bioinformatic characterization of eukaryotic selenoproteomes, and can complement genome annotation pipelines.

Availability and Implementation: Selenoprofiles is a python-built pipeline that internally runs psitblastn, exonerate, genewise, SECISearch and a number of custom-made scripts and programs. The program is available at http://big.crg.cat/services/selenoprofiles. The predictions presented in this article are available through DAS at http://genome.crg.cat:9000/das/Selenoprofiles_ensembl.

Contact:marco.mariotti@crg.es

Supplementary information:Supplementary data are available at Bioinformatics online.

Genes for selenium dependent and independent formate dehydrogenase in the gut microbial communities of three lower, wood-feeding termites and a wood-feeding roach

The bacterial Wood-Ljungdahl pathway for CO(2)-reductive acetogenesis is important for the nutritional mutualism occurring between wood-feeding insects and their hindgut microbiota. A key step in this pathway is the reduction of CO(2) to formate, catalysed by the enzyme formate dehydrogenase (FDH). Putative selenocysteine- (Sec) and cysteine- (Cys) containing paralogues of hydrogenase-linked FDH (FDH(H)) have been identified in the termite gut acetogenic spirochete, Treponema primitia, but knowledge of their relevance in the termite gut environment remains limited. In this study, we designed degenerate PCR primers for FDH(H) genes (fdhF) and assessed fdhF diversity in insect gut bacterial isolates and the gut microbial communities of termites and cockroaches. The insects examined herein represent three wood-feeding termite families, Termopsidae, Kalotermitidae and Rhinotermitidae (phylogenetically 'lower' termite taxa); the wood-feeding roach family Cryptocercidae (the sister taxon to termites); and the omnivorous roach family Blattidae. Sec and Cys FDH(H) variants were identified in every wood-feeding insect but not the omnivorous roach. Of 68 novel alleles obtained from inventories, 66 affiliated phylogenetically with enzymes from T. primitia. These formed two subclades (37 and 29 phylotypes) almost completely comprised of Sec-containing and Cys-containing enzymes respectively. A gut cDNA inventory showed transcription of both variants in the termite Zootermopsis nevadensis (family Termopsidae). The gene patterns suggest that FDH(H) enzymes are important for the CO(2)-reductive metabolism of uncultured acetogenic treponemes and imply that the availability of selenium, a trace element, shaped microbial gene content in the last common ancestor of dictyopteran, wood-feeding insects, and continues to shape it to this day.

Syntrophic butyrate and propionate oxidation processes: from genomes to reaction mechanisms

In anoxic environments such as swamps, rice fields and sludge digestors, syntrophic microbial communities are important for decomposition of organic matter to CO2 and CH4 . The most difficult step is the fermentative degradation of short-chain fatty acids such as propionate and butyrate. Conversion of these metabolites to acetate, CO2 , formate and hydrogen is endergonic under standard conditions and occurs only if methanogens keep the concentrations of these intermediate products low. Butyrate and propionate degradation pathways include oxidation steps of comparably high redox potential, i.e. oxidation of butyryl-CoA to crotonyl-CoA and of succinate to fumarate, respectively, that require investment of energy to release the electrons as hydrogen or formate. Although investigated for several decades, the biochemistry of these reactions is still not completely understood. Genome analysis of the butyrate-oxidizing Syntrophomonas wolfei and Syntrophus aciditrophicus and of the propionate-oxidizing Syntrophobacter fumaroxidans and Pelotomaculum thermopropionicum reveals the presence of energy-transforming protein complexes. Recent studies indicated that S. wolfei uses electron-transferring flavoproteins coupled to a menaquinone loop to drive butyryl-CoA oxidation, and that S. fumaroxidans uses a periplasmic formate dehydrogenase, cytochrome b:quinone oxidoreductases, a menaquinone loop and a cytoplasmic fumarate reductase to drive energy-dependent succinate oxidation. Furthermore, we propose that homologues of the Thermotoga maritima bifurcating [FeFe]-hydrogenase are involved in NADH oxidation by S. wolfei and S. fumaroxidans to form hydrogen.

Selenium controls transcription of paralogous formate dehydrogenase genes in the termite gut acetogen, Treponema primitia

The termite gut spirochete, Treponema primitia, is a CO(2)-reductive acetogen that is phylogenetically distinct from other distantly related and more extensively studied acetogens such as Moorella thermoacetica. Research on T. primitia has revealed details about the role of spirochetes in CO(2)-reductive acetogenesis, a process important to the mutualism occurring between termites and their gut microbial communities. Here, a locus of the T. primitia genome containing Wood-Ljungdahl pathway genes for CO(2)-reductive acetogenesis was sequenced. This locus contained methyl-branch genes of the pathway (i.e. for the reduction of CO(2) to the level of methyl-tetrahydrofolate) including paralogous genes for cysteine and selenocysteine (Sec) variants of formate dehydrogenase (FDH) and genes for Sec incorporation. The FDH variants affiliated phylogenetically with hydrogenase-linked FDH enzymes, suggesting that T. primitia FDH enzymes utilize electrons derived directly from molecular H(2). Sub-nanomolar concentrations of selenium decreased transcript levels of the cysteine variant FDH gene. Selenium concentration did not markedly influence the level of mRNA upstream of the Sec-codon in the Sec variant FDH; however, the level of transcript extending downstream of the Sec-codon increased incrementally with increasing selenium concentrations. The features and regulation of these FDH genes are an indication that T. primitia may experience dynamic selenium availability in its H(2)-rich gut environment.

Targeting selenium metabolism and selenoproteins: novel avenues for drug discovery

Selenoproteins play a wide range of roles in metabolism and oxidative stress defense and are produced by organisms in all three domains of life. Recent evidence has been presented that metal based cancer drugs target the selenol nucleophile of the active site selenocysteine in thioredoxin reductase isoenzymes. Other metals and metalloids, such as tin, arsenic and gold, have also recently been shown to form stable complexes with hydrogen selenide, a required precursor for the synthesis of selenoproteins in all biological organisms. Moreover these metal based compounds have been shown to inhibit growth of pathogens such as Clostridium difficile and Treponema denticola due to their reactivity with this highly reactive metabolic precursor. This review summarizes the recent finding on these two avenues for drug discovery, and puts this work in context with the larger field of selenium biology.

dbTEU: a protein database of trace element utilization

Biological trace elements are required for numerous biological processes and by all organisms. We describe a database, dbTEU (DataBase of Trace Element Utilization), that features known transporters and user proteins for five trace elements (copper, molybdenum, nickel, cobalt and selenium) and represents sequenced organisms from the three domains of life. The manually curated dbTEU currently includes approximately 16,500 proteins from >700 organisms, and offers interactive trace element, protein, organism and sequence search and browse tools.

AVAILABILITY AND IMPLEMENTATION

dbTEU is freely available at http://gladyshevlab.bwh.harvard.edu/trace_element/

A method for identification of selenoprotein genes in archaeal genomes

The genetic codon UGA has a dual function: serving as a terminator and encoding selenocysteine. However, most popular gene annotation programs only take it as a stop signal, resulting in misannotation or completely missing selenoprotein genes. We developed a computational method named Asec-Prediction that is specific for the prediction of archaeal selenoprotein genes. To evaluate its effectiveness, we first applied it to 14 archaeal genomes with previously known selenoprotein genes, and Asec-Prediction identified all reported selenoprotein genes without redundant results. When we applied it to 12 archaeal genomes that had not been researched for selenoprotein genes, Asec-Prediction detected a novel selenoprotein gene in Methanosarcina acetivorans. Further evidence was also collected to support that the predicted gene should be a real selenoprotein gene. The result shows that Asec-Prediction is effective for the prediction of archaeal selenoprotein genes.

General trends in trace element utilization revealed by comparative genomic analyses of Co, Cu, Mo, Ni, and Se

Trace elements are used by all organisms and provide proteins with unique coordination and catalytic and electron transfer properties. Although many trace element-containing proteins are well characterized, little is known about the general trends in trace element utilization. We carried out comparative genomic analyses of copper, molybdenum, nickel, cobalt (in the form of vitamin B(12)), and selenium (in the form of selenocysteine) in 747 sequenced organisms at the following levels: (i) transporters and transport-related proteins, (ii) cofactor biosynthesis traits, and (iii) trace element-dependent proteins. Few organisms were found to utilize all five trace elements, whereas many symbionts, parasites, and yeasts used only one or none of these elements. Investigation of metalloproteomes and selenoproteomes revealed examples of increased utilization of proteins that use copper in land plants, cobalt in Dehalococcoides and Dictyostelium, and selenium in fish and algae, whereas nematodes were found to have great diversity of copper transporters. These analyses also characterized trace element metabolism in common model organisms and suggested new model organisms for experimental studies of individual trace elements. Mismatches in the occurrence of user proteins and corresponding transport systems revealed deficiencies in our understanding of trace element biology. Biological interactions among some trace elements were observed; however, such links were limited, and trace elements generally had unique utilization patterns. Finally, environmental factors, such as oxygen requirement and habitat, correlated with the utilization of certain trace elements. These data provide insights into the general features of utilization and evolution of trace elements in the three domains of life.

Involvement of NADH:acceptor oxidoreductase and butyryl coenzyme A dehydrogenase in reversed electron transport during syntrophic butyrate oxidation by Syntrophomonas wolfei

Methanogenic oxidation of butyrate to acetate requires a tight cooperation between the syntrophically fermenting Syntrophomonas wolfei and the methanogen Methanospirillum hungatei, and a reversed electron transport system in S. wolfei was postulated to shift electrons from butyryl coenzyme A (butyryl-CoA) oxidation to the redox potential of NADH for H(2) generation. The metabolic activity of butyrate-oxidizing S. wolfei cells was measured via production of formazan and acetate from butyrate, with 2,3,5-triphenyltetrazolium chloride as electron acceptor. This activity was inhibited by trifluoperazine (TPZ), an antitubercular agent known to inhibit NADH:menaquinone oxidoreductase. In cell extracts of S. wolfei, the oxidation of NADH could be measured with quinones, viologens, and tetrazolium dyes as electron acceptors, and also this activity was inhibited by TPZ. The TPZ-sensitive NADH:acceptor oxidoreductase activity appeared to be membrane associated but could be dissociated from the membrane as a soluble protein and was semipurified by anion-exchange chromatography. Recovered proteins were identified by peptide mass fingerprinting, which indicated the presence of an NADH:acceptor oxidoreductase as part of a three-component [FeFe] hydrogenase complex and a selenocysteine-containing formate dehydrogenase. Furthermore, purification of butyryl-CoA dehydrogenase (Bcd) activity and peptide mass fingerprinting revealed two Bcd proteins different from the Bcd subunit of the Bcd/electron-transfer flavoprotein complex (Bcd/EtfAB) predicted from the genome sequence of S. wolfei. The results suggest that syntrophic oxidation of butyrate in S. wolfei involves a membrane-associated TPZ-sensitive NADH:acceptor oxidoreductase as part of a hydrogenase complex similar to the recently discovered "bifurcating" hydrogenase in Thermotoga maritima and butyryl-CoA dehydrogenases that are different from Bcd of the Bcd/EtfAB complex.

The selenoproteome of Clostridium sp. OhILAs: characterization of anaerobic bacterial selenoprotein methionine sulfoxide reductase A

Selenocysteine (Sec) is incorporated into proteins in response to UGA codons. This residue is frequently found at the catalytic sites of oxidoreductases. In this study, we characterized the selenoproteome of an anaerobic bacterium, Clostridium sp. (also known as Alkaliphilus oremlandii) OhILA, and identified 13 selenoprotein genes, five of which have not been previously described. One of the detected selenoproteins was methionine sulfoxide reductase A (MsrA), an antioxidant enzyme that repairs oxidatively damaged methionines in a stereospecific manner. To date, little is known about MsrA from anaerobes. We characterized this selenoprotein MsrA which had a single Sec residue at the catalytic site but no cysteine (Cys) residues in the protein sequence. Its SECIS (Sec insertion sequence) element did not resemble those in Escherichia coli. Although with low translational efficiency, the expression of the Clostridium selenoprotein msrA gene in E. coli could be demonstrated by (75)Se metabolic labeling, immunoblot analyses, and enzyme assays, indicating that its SECIS element was recognized by the E. coli Sec insertion machinery. We found that the Sec-containing MsrA exhibited at least a 20-fold higher activity than its Cys mutant form, indicating a critical role of Sec in the catalytic activity of the enzyme. Furthermore, our data revealed that the Clostridium MsrA was inefficiently reducible by thioredoxin, which is a typical reducing agent for MsrA, suggesting the use of alternative electron donors in this anaerobic bacterium that directly act on the selenenic acid intermediate and do not require resolving Cys residues.

Trends in selenium utilization in marine microbial world revealed through the analysis of the global ocean sampling (GOS) project

Selenium is an important trace element that occurs in proteins in the form of selenocysteine (Sec) and in tRNAs in the form of selenouridine. Recent large-scale metagenomics projects provide an opportunity for understanding global trends in trace element utilization. Herein, we characterized the selenoproteome of the microbial marine community derived from the Global Ocean Sampling (GOS) expedition. More than 3,600 selenoprotein gene sequences belonging to 58 protein families were detected, including sequences representing 7 newly identified selenoprotein families, such as homologs of ferredoxin–thioredoxin reductase and serine protease. In addition, a new eukaryotic selenoprotein family, thiol reductase GILT, was identified. Most GOS selenoprotein families originated from Cys-containing thiol oxidoreductases. In both Pacific and Atlantic microbial communities, SelW-like and SelD were the most widespread selenoproteins. Geographic location had little influence on Sec utilization as measured by selenoprotein variety and the number of selenoprotein genes detected; however, both higher temperature and marine (as opposed to freshwater and other aquatic) environment were associated with increased use of this amino acid. Selenoproteins were also detected with preference for either environment. We identified novel fusion forms of several selenoproteins that highlight redox activities of these proteins. Almost half of Cys-containing SelDs were fused with NADH dehydrogenase, whereas such SelD forms were rare in terrestrial organisms. The selenouridine utilization trait was also analyzed and showed an independent evolutionary relationship with Sec utilization. Overall, our study provides insights into global trends in microbial selenium utilization in marine environments.

Selenium (Se) is an essential micronutrient due to its requirement for biosynthesis and function of the 21^st amino acid, selenocysteine (Sec). Sec is found in the active sites of selenoproteins, most of which exhibit redox function, in all three domains of life. In recent years, genome sequencing projects provided a large volume of nucleotide and protein sequence information. Identification of complete sets of selenoproteins (selenoproteomes) of individual organisms and environmental samples is important for better understanding of Se utilization, biological functions of this element, and changes in Se use during evolution. Here, we describe a comprehensive analysis of the selenoproteome of the microbial marine community derived from the Global Ocean Sampling (GOS) expedition. More than 3,600 selenoprotein gene sequences belonging to 58 protein families were detected and analyzed. Our study generated the largest selenoproteome reported to date and provided important insights into microbial Se utilization and its evolutionary trends in marine environments.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteins usually misidentified and misannotated in sequence databases. The presence of an in-frame UGA (stop) codon in the coding sequence of selenoprotein genes precludes their identification and correct annotation. The in-frame UGA codons are recoded to cotranslationally incorporate selenocysteine, a rare selenium-containing amino acid. The development of ad hoc experimental and, more recently, computational approaches have allowed the efficient identification and characterization of the selenoproteomes of a growing number of species. Today, dozens of selenoprotein families have been described and more are being discovered in recently sequenced species, but the correct genomic annotation is not available for the majority of these genes. SelenoDB is a long-term project that aims to provide, through the collaborative effort of experimental and computational researchers, automatic and manually curated annotations of selenoprotein genes, proteins and SECIS elements. Version 1.0 of the database includes an initial set of eukaryotic genomic annotations, with special emphasis on the human selenoproteome, for immediate inspection by selenium researchers or incorporation into more general databases. SelenoDB is freely available at http://www.selenodb.org.

Factors and selenocysteine insertion sequence requirements for the synthesis of selenoproteins from a gram-positive anaerobe in Escherichia coli

Selenoprotein synthesis in Escherichia coli strictly depends on the presence of a specific selenocysteine insertion sequence (SECIS) following the selenocysteine-encoding UGA codon of the respective mRNA. It is recognized by the selenocysteine-specific elongation factor SelB, leading to cotranslational insertion of selenocysteine into the nascent polypeptide chain. The synthesis of three different selenoproteins from the gram-positive anaerobe Eubacterium acidaminophilum in E. coli was studied. Incorporation of (75)Se into glycine reductase protein B (GrdB1), the peroxiredoxin PrxU, and selenophosphate synthetase (SelD1) was negligible in an E. coli wild-type strain and was fully absent in an E. coli SelB mutant. Selenoprotein synthesis, however, was strongly increased if selB and selC (tRNA(Sec)) from E. acidaminophilum were coexpressed. Putative secondary structures downstream of the UGA codons did not show any sequence similarity to each other or to the E. coli SECIS element. However, mutations in these structures strongly reduced the amount of (75)Se-labeled protein, indicating that they indeed act as SECIS elements. UGA readthrough mediated by the three different SECIS elements was further analyzed using gst-lacZ translational fusions. In the presence of selB and selC from E. acidaminophilum, UGA readthrough was 36 to 64% compared to the respective cysteine-encoding UGC variant. UGA readthrough of SECIS elements present in Desulfomicrobium baculatum (hydV), Treponema denticola (selD), and Campylobacter jejuni (selW-like gene) was also considerably enhanced in the presence of E. acidaminophilum selB and selC. This indicates recognition of these SECIS elements and might open new perspectives for heterologous selenoprotein synthesis in E. coli.

Identification and characterization of a selenoprotein family containing a diselenide bond in a redox motif

Selenocysteine (Sec, U) insertion into proteins is directed by translational recoding of specific UGA codons located upstream of a stem-loop structure known as Sec insertion sequence (SECIS) element. Selenoproteins with known functions are oxidoreductases containing a single redox-active Sec in their active sites. In this work, we identified a family of selenoproteins, designated SelL, containing two Sec separated by two other residues to form a UxxU motif. SelL proteins show an unusual occurrence, being present in diverse aquatic organisms, including fish, invertebrates, and marine bacteria. Both eukaryotic and bacterial SelL genes use single SECIS elements for insertion of two Sec. In eukaryotes, the SECIS is located in the 3' UTR, whereas the bacterial SelL SECIS is within a coding region and positioned at a distance that supports the insertion of either of the two Sec or both of these residues. SelL proteins possess a thioredoxin-like fold wherein the UxxU motif corresponds to the catalytic CxxC motif in thioredoxins, suggesting a redox function of SelL proteins. Distantly related SelL-like proteins were also identified in a variety of organisms that had either one or both Sec replaced with Cys. Danio rerio SelL, transiently expressed in mammalian cells, incorporated two Sec and localized to the cytosol. In these cells, it occurred in an oxidized form and was not reducible by DTT. In a bacterial expression system, we directly demonstrated the formation of a diselenide bond between the two Sec, establishing it as the first diselenide bond found in a natural protein.

High content of proteins containing 21st and 22nd amino acids, selenocysteine and pyrrolysine, in a symbiotic deltaproteobacterium of gutless worm Olavius algarvensis

Selenocysteine (Sec) and pyrrolysine (Pyl) are rare amino acids that are cotranslationally inserted into proteins and known as the 21st and 22nd amino acids in the genetic code. Sec and Pyl are encoded by UGA and UAG codons, respectively, which normally serve as stop signals. Herein, we report on unusually large selenoproteomes and pyrroproteomes in a symbiont metagenomic dataset of a marine gutless worm, Olavius algarvensis. We identified 99 selenoprotein genes that clustered into 30 families, including 17 new selenoprotein genes that belong to six families. In addition, several Pyl-containing proteins were identified in this dataset. Most selenoproteins and Pyl-containing proteins were present in a single deltaproteobacterium, δ1 symbiont, which contained the largest number of both selenoproteins and Pyl-containing proteins of any organism reported to date. Our data contrast with the previous observations that symbionts and host-associated bacteria either lose Sec utilization or possess a limited number of selenoproteins, and suggest that the environment in the gutless worm promotes Sec and Pyl utilization. Anaerobic conditions and consistent selenium supply might be the factors that support the use of amino acids that extend the genetic code.

Mining prokaryotic genomes for unknown amino acids: a stop-codon-based approach

Background

Selenocysteine and pyrrolysine are the 21st and 22nd amino acids, which are genetically encoded by stop codons. Since a number of microbial genomes have been completely sequenced to date, it is tempting to ask whether the 23rd amino acid is left undiscovered in these genomes. Recently, a computational study addressed this question and reported that no tRNA gene for unknown amino acid was found in genome sequences available. However, performance of the tRNA prediction program on an unknown tRNA family, which may have atypical sequence and structure, is unclear, thereby rendering their result inconclusive. A protein-level study will provide independent insight into the novel amino acid.

Results

Assuming that the 23rd amino acid is also encoded by a stop codon, we systematically predicted proteins that contain stop-codon-encoded amino acids from 191 prokaryotic genomes. Since our prediction method relies only on the conservation patterns of primary sequences, it also provides an opportunity to search novel selenoproteins and other readthrough proteins. It successfully recovered many of currently known selenoproteins and pyrrolysine proteins. However, no promising candidate for the 23rd amino acid was detected, and only one novel selenoprotein was predicted.

Conclusion

Our result suggests that the unknown amino acid encoded by stop codons does not exist, or its phylogenetic distribution is rather limited, which is in agreement with the previous study on tRNA. The method described here can be used in future studies to explore novel readthrough events from complete genomes, which are rapidly growing.

A recoding element that stimulates decoding of UGA codons by Sec tRNA[Ser]Sec

Selenocysteine insertion during decoding of eukaryotic selenoprotein mRNA requires several trans-acting factors and a cis-acting selenocysteine insertion sequence (SECIS) usually located in the 3' UTR. A second cis-acting selenocysteine codon redefinition element (SRE) has recently been described that resides near the UGA-Sec codon of selenoprotein N (SEPN1). Similar phylogenetically conserved elements can be predicted in a subset of eukaryotic selenoprotein mRNAs. Previous experimental analysis of the SEPN1 SRE revealed it to have a stimulatory effect on readthrough of the UGA-Sec codon, which was not dependent upon the presence of a SECIS element in the 3' UTR; although, as expected, readthrough efficiency was further elevated by inclusion of a SECIS. In order to examine the nature of the redefinition event stimulated by the SEPN1 SRE, we have modified an experimentally tractable in vitro translation system that recapitulates efficient selenocysteine insertion. The results presented here illustrate that the SRE element has a stimulatory effect on decoding of the UGA-Sec codon by both the methylated and unmethylated isoforms of Sec tRNA([Ser]Sec), and confirm that efficient selenocysteine insertion is dependent on the presence of a 3'-UTR SECIS. The variation in recoding elements predicted near UGA-Sec codons implies that these elements may play a differential role in determining the amount of selenoprotein produced by acting as controllers of UGA decoding efficiency.

New selenoproteins identified in silico from the genome of Anopheles gambiae

Selenoprotein is biosynthesized by the incorporation of selenocysteine into proteins, where the TGA codon in the open reading frame does not act as a stop signal but is translated into selenocysteine. The dual functions of TGA result in mis-annotation or lack of selenoproteins in the sequenced genomes of many species. Available computational tools fail to correctly predict selenoproteins. Thus, we developed a new method to identify selenoproteins from the genome of Anopheles gambiae computationally. Based on released genomic information, several programs were edited with PERL language to identify selenocysteine insertion sequence (SECIS) element, the coding potential of TGA codons, and cysteine-containing homologs of selenoprotein genes. Our results showed that 11365 genes were terminated with TGA codons, 918 of which contained SECIS elements. Similarity search revealed that 58 genes contained Sec/Cys pairs and similar flanking regions around in-frame TGA codons. Finally, 7 genes were found to fully meet requirements for selenoproteins, although they have not been annotated as selenoproteins in NCBI databases. Deduced from their basic properties, the newly found selenoproteins in the genome of Anopheles gambiae are possibly related to in vivo oxidation tolerance and protein regulation in order to interfere with anopheles' vectorial capacity of Plasmodium. This study may also provide theoretical bases for the prevention of malaria from anopheles transmission.

Dynamic evolution of selenocysteine utilization in bacteria: a balance between selenoprotein loss and evolution of selenocysteine from redox active cysteine residues

Comparative genomics and evolutionary analyses to examine the dynamics of selenocysteine utilization in bacteria reveal a dynamic balance between selenoprotein origin and loss.

Background

Selenocysteine (Sec) is co-translationally inserted into protein in response to UGA codons. It occurs in oxidoreductase active sites and often is catalytically superior to cysteine (Cys). However, Sec is used very selectively in proteins and organisms. The wide distribution of Sec and its restricted use have not been explained.

Results

We conducted comparative genomics and phylogenetic analyses to examine dynamics of Sec decoding in bacteria at both selenium utilization trait and selenoproteome levels. These searches revealed that 21.5% of sequenced bacteria utilize Sec, their selenoproteomes have 1 to 31 selenoproteins, and selenoprotein-rich organisms are mostly Deltaproteobacteria or Firmicutes/Clostridia. Evolutionary histories of selenoproteins suggest that Cys-to-Sec replacement is a general trend for most selenoproteins. In contrast, only a small number of Sec-to-Cys replacements were detected, and these were mostly restricted to formate dehydrogenase and selenophosphate synthetase families. In addition, specific selenoprotein gene losses were observed in many sister genomes. Thus, the Sec/Cys replacements were mostly unidirectional, and increased utilization of Sec by existing protein families was counterbalanced by loss of selenoprotein genes or entire selenoproteomes. Lateral transfers of the Sec trait were an additional factor, and we describe the first example of selenoprotein gene transfer between archaea and bacteria. Finally, oxygen requirement and optimal growth temperature were identified as environmental factors that correlate with changes in Sec utilization.

Conclusion

Our data reveal a dynamic balance between selenoprotein origin and loss, and may account for the discrepancy between catalytic advantages provided by Sec and the observed low number of selenoprotein families and Sec-utilizing organisms.

Understanding the importance of selenium and selenoproteins in muscle function

Selenium is an essential trace element. In cattle, selenium deficiency causes dysfunction of various organs, including skeletal and cardiac muscles. In humans as well, lack of selenium is associated with many disorders, but despite accumulation of clinical reports, muscle diseases are not generally considered on the list. The goal of this review is to establish the connection between clinical observations and the most recent advances obtained in selenium biology. Recent results about a possible role of selenium-containing proteins in muscle formation and repair have been collected. Selenoprotein N is the first selenoprotein linked to genetic disorders consisting of different forms of congenital muscular dystrophies. Understanding the muscle disorders associated with selenium deficiency or selenoprotein N dysfunction is an essential step in defining the causes of the disease and obtaining a better comprehension of the mechanisms involved in muscle formation and maintenance.

Analysis of the cDNAs of hypothetical genes on Arabidopsis chromosome 2 reveals numerous transcript variants

In the fully sequenced Arabidopsis (Arabidopsis thaliana) genome, many gene models are annotated as "hypothetical protein," whose gene structures are predicted solely by computer algorithms with no support from either expressed sequence matches from Arabidopsis, or nucleic acid or protein homologs from other species. In order to confirm their existence and predicted gene structures, a high-throughput method of rapid amplification of cDNA ends (RACE) was used to obtain their cDNA sequences from 11 cDNA populations. Primers from all of the 797 hypothetical genes on chromosome 2 were designed, and, through 5' and 3' RACE, clones from 506 genes were sequenced and cDNA sequences from 399 target genes were recovered. The cDNA sequences were obtained by assembling their 5' and 3' RACE polymerase chain reaction products. These sequences revealed that (1) the structures of 151 hypothetical genes were different from their predictions; (2) 116 hypothetical genes had alternatively spliced transcripts and 187 genes displayed polyadenylation sites; and (3) there were transcripts arising from both strands, from the strand opposite to that of the prediction and possible dicistronic transcripts. Promoters from five randomly chosen hypothetical genes (At2g02540, At2g31270, At2g33640, At2g35550, and At2g36340) were cloned into report constructs, and their expressions are tissue or development stage specific. Our results indicate at least 50% of hypothetical genes on chromosome 2 are expressed in the cDNA populations with about 38% of the gene structures differing from their predictions. Thus, by using this targeted approach, high-throughput RACE, we revealed numerous transcripts including many uncharacterized variants from these hypothetical genes.

A computational method to predict genetically encoded rare amino acids in proteins

A new method for predicting recoding by rare amino acids such as selenocysteine and pyrrolysine was used to survey a set of microbial genomes.

In several natural settings, the standard genetic code is expanded to incorporate two additional amino acids with distinct functionality, selenocysteine and pyrrolysine. These rare amino acids can be overlooked inadvertently, however, as they arise by recoding at certain stop codons. We report a method for such recoding prediction from genomic data, using read-through similarity evaluation. A survey across a set of microbial genomes identifies almost all the known cases as well as a number of novel candidate proteins.

Selenocysteine in proteins-properties and biotechnological use

Selenocysteine (Sec), the 21st amino acid, exists naturally in all kingdoms of life as the defining entity of selenoproteins. Sec is a cysteine (Cys) residue analogue with a selenium-containing selenol group in place of the sulfur-containing thiol group in Cys. The selenium atom gives Sec quite different properties from Cys. The most obvious difference is the lower pK(a) of Sec, and Sec is also a stronger nucleophile than Cys. Proteins naturally containing Sec are often enzymes, employing the reactivity of the Sec residue during the catalytic cycle and therefore Sec is normally essential for their catalytic efficiencies. Other unique features of Sec, not shared by any of the other 20 common amino acids, derive from the atomic weight and chemical properties of selenium and the particular occurrence and properties of its stable and radioactive isotopes. Sec is, moreover, incorporated into proteins by an expansion of the genetic code as the translation of selenoproteins involves the decoding of a UGA codon, otherwise being a termination codon. In this review, we will describe the different unique properties of Sec and we will discuss the prerequisites for selenoprotein production as well as the possible use of Sec introduction into proteins for biotechnological applications. These include residue-specific radiolabeling with gamma or positron emitters, the use of Sec as a reactive handle for electophilic probes introducing fluorescence or other peptide conjugates, as the basis for affinity purification of recombinant proteins, the trapping of folding intermediates, improved phasing in X-ray crystallography, introduction of 77Se for NMR spectroscopy, or, finally, the analysis or tailoring of enzymatic reactions involving thiol or oxidoreductase (redox) selenolate chemistry.

Making sense of nonsense: the evolution of selenocysteine usage in proteins

A recent analysis of sequences derived from organisms in the Sargasso Sea has revealed a surprisingly different set of selenium-containing proteins than that previously found in sequenced genomes.

A recent analysis of sequences derived from organisms in the Sargasso Sea has revealed a surprisingly different set of selenium-containing proteins than that previously found in sequenced genomes and suggests that selenocysteine utilization has been lost by many groups of organisms during evolution.

The microbial selenoproteome of the Sargasso Sea

An analysis of the selenoproteome of the largest microbial sequence dataset, the Sargasso Sea environmental genome sequences, identified 310 selenoprotein genes that clustered into 25 families. This included 101 new selenoprotein genes that belonged to 15 families, doubling the number of prokaryotic selenoprotein families.

Background

Selenocysteine (Sec) is a rare amino acid which occurs in proteins in major domains of life. It is encoded by TGA, which also serves as the signal for termination of translation, precluding identification of selenoprotein genes by available annotation tools. Information on full sets of selenoproteins (selenoproteomes) is essential for understanding the biology of selenium. Herein, we characterized the selenoproteome of the largest microbial sequence dataset, the Sargasso Sea environmental genome project.

Results

We identified 310 selenoprotein genes that clustered into 25 families, including 101 new selenoprotein genes that belonged to 15 families. Most of these proteins were predicted redox proteins containing catalytic selenocysteines. Several bacterial selenoproteins previously thought to be restricted to eukaryotes were detected by analyzing eukaryotic and bacterial SECIS elements, suggesting that eukaryotic and bacterial selenoprotein sets partially overlapped. The Sargasso Sea microbial selenoproteome was rich in selenoproteins and its composition was different from that observed in the combined set of completely sequenced genomes, suggesting that these genomes do not accurately represent the microbial selenoproteome. Most detected selenoproteins occurred sporadically compared to the widespread presence of their cysteine homologs, suggesting that many selenoproteins recently evolved from cysteine-containing homologs.

Conclusions

This study yielded the largest selenoprotein dataset to date, doubled the number of prokaryotic selenoprotein families and provided insights into forces that drive selenocysteine evolution.