Search

Selenium is a fascinating element that has a long history, most of which documents it as a deleterious element to health. In more recent years, selenium has been found to be an essential element in the diet of humans, all other mammals, and many other life forms. It has many health benefits that include, for example, roles in preventing heart disease and certain forms of cancer, slowing AIDS progression in HIV patients, supporting male reproduction, inhibiting viral expression, and boosting the immune system, and it also plays essential roles in mammalian development. Elucidating the molecular biology of selenium over the past 40 years generated an entirely new field of science which encompassed the many novel features of selenium. These features were (1) how this element makes its way into protein as the 21st amino acid in the genetic code, selenocysteine (Sec); (2) the vast amount of machinery dedicated to synthesizing Sec uniquely on its tRNA; (3) the incorporation of Sec into protein; and (4) the roles of the resulting Sec-containing proteins (selenoproteins) in health and development. One of the research areas receiving the most attention regarding selenium in health has been its role in cancer prevention, but further research has also exposed the role of this element as a facilitator of various maladies, including cancer.

Selenium is incorporated into selenoproteins as the 21st amino acid selenocysteine (Sec). There are 25 selenoproteins encoded in the human genome, and their synthesis requires a dedicated machinery. Most selenoproteins are oxidoreductases with important functions in human health. A number of disorders have been associated with deficiency of selenoproteins, caused by mutations in selenoprotein genes or Sec machinery genes. We discuss mutations that are known to cause disease in humans and report their allele frequencies in the general population. The occurrence of protein-truncating variants in the same genes is also presented. We provide an overview of pathogenic variants in selenoproteins genes from a population genomics perspective.

Mouse has emerged as the most common model organism in biomedicine. Here, we analyzed the tolerance to the loss-of-function (LoF) of selenoprotein genes, estimated from mouse knockouts and the frequency of LoF variants in humans. We found not only a general correspondence in tolerance (e.g., GPX1, GPX2) and intolerance (TXNRD1, SELENOT) to gene LoF between humans and mice but also important differences. Notably, humans are intolerant to the loss of iodothyronine deiodinases, whereas their deletion in mice leads to mild phenotypes, and this is consistent with phenotype differences in selenocysteine machinery loss between these species. In contrast, loss of TXNRD2 and GPX4 is lethal in mice but may be tolerated in humans. We further identified the first human SELENOP variants coding for proteins varying in selenocysteine content. Finally, our analyses suggested that premature termination codons in selenoprotein genes trigger nonsense-mediated decay, but do this inefficiently when UGA codon is gained. Overall, our study highlights differences in the physiological importance of selenoproteins between human and mouse.

Selenoproteins typically contain a single selenocysteine, the 21st amino acid, encoded by a context-redefined UGA. However, human selenoprotein P (SelenoP) has a redox-functioning selenocysteine in its N-terminal domain and nine selenium transporter-functioning selenocysteines in its C-terminal domain. Here we show that diverse SelenoP genes are present across metazoa with highly variable numbers of Sec-UGAs, ranging from a single UGA in certain insects, to 9 in common spider, and up to 132 in bivalve molluscs. SelenoP genes were shaped by a dynamic evolutionary process linked to selenium usage. Gene evolution featured modular expansions of an ancestral multi-Sec domain, which led to particularly Sec-rich SelenoP proteins in many aquatic organisms. We focused on molluscs, and chose Pacific oyster Magallana gigas as experimental model. We show that oyster SelenoP mRNA with 46 UGAs is translated full-length in vivo. Ribosome profiling indicates that selenocysteine specification occurs with ~5% efficiency at UGA1 and approaches 100% efficiency at distal 3' UGAs. We report genetic elements relevant to its expression, including a leader open reading frame and an RNA structure overlapping the initiation codon that modulates ribosome progression in a selenium-dependent manner. Unlike their mammalian counterparts, the two SECIS elements in oyster SelenoP (3'UTR recoding elements) do not show functional differentiation in vitro. Oysters can increase their tissue selenium level up to 50-fold upon supplementation, which also results in extensive changes in selenoprotein expression. The work was supported by the following grants: NIH DK117149, GM065204, and AG021518 (V.N.G.); BIO2014-57291-R from the Spanish Ministry of Economy and Competitiveness (R.G.), NIH GM077073 (P.R.C.), Science Foundation Ireland (13/IA/1853), and Irish Research Council (IRCLA/2019/74 n) (J.F.A).

Selenocysteine (Sec) is the 21st amino acid in the genetic code, inserted in response to UGA codons with the help of RNA structures, the SEC Insertion Sequence (SECIS) elements. The three domains of life feature distinct strategies for Sec insertion in proteins and its utilization. While bacteria and archaea possess similar sets of selenoproteins, Sec biosynthesis is more similar among archaea and eukaryotes. However, SECIS elements are completely different in the three domains of life. Here, we analyze the archaeon Lokiarchaeota that resolves the relationships among Sec insertion systems. This organism has selenoproteins representing five protein families, three of which have multiple Sec residues. Remarkably, these archaeal selenoprotein genes possess conserved RNA structures that strongly resemble the eukaryotic SECIS element, including key eukaryotic protein-binding sites. These structures also share similarity with the SECIS element in archaeal selenoprotein VhuD, suggesting a relation of direct descent. These results identify Lokiarchaeota as an intermediate form between the archaeal and eukaryotic Sec-encoding systems and clarify the evolution of the Sec insertion system. This study was supported by NIH GM061603, GM065205 and CA080946. B.M. is partly supported by The Pew Charitable Trust postdoctoral fellow program.

Selenoproteins contain selenocysteine (Sec or U), the 21st amino acid, inserted in response to an in-frame UGA codon. UGA normally terminates translation, but in selenoprotein mRNAs it is recoded to specify Sec insertion. For this reason, standard gene prediction programs fail to predict Sec codons, and selenoproteins are usually misannotated in protein databases and genome projects. Selenoprofiles is a computational pipeline able to correctly annotate selenoprotein genes in genomic sequences. This program uses a SECIS-independent approach, based on homology searches, and employs curated built-in profile alignments for all known selenoprotein families. Selenoprofiles constitutes the most accurate method for predicting selenoprotein genes belonging to known families.

Selenocysteine (Sec) is known as the 21st amino acid, a cysteine analogue with selenium replacing sulphur. Sec is inserted co-translationally in a small fraction of proteins called selenoproteins. In selenoprotein genes, the Sec specific tRNA (tRNASec) drives the recoding of highly specific UGA codons from stop signals to Sec. Although found in organisms from the three domains of life, Sec is not universal. Many species are completely devoid of selenoprotein genes and lack the ability to synthesize Sec. Since tRNASec is a key component in selenoprotein biosynthesis, its efficient identification in genomes is instrumental to characterize the utilization of Sec across lineages. Available tRNA prediction methods fail to accurately predict tRNASec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNASec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNASec predictions, which is not possible for the predictions of other tRNAs. We show that Secmarker greatly improves the accuracy of previously existing methods constituting a valuable tool to identify tRNASec genes, and to efficiently determine whether a genome contains selenoproteins. We used Secmarker to analyze a large set of fully sequenced genomes, and the results revealed new insights in the biology of tRNASec, led to the discovery of a novel bacterial selenoprotein family, and shed additional light on the phylogenetic distribution of selenoprotein containing genomes. Secmarker is freely accessible for download, or online analysis through a web server at http://secmarker.crg.cat., Author summary Most proteins are made of twenty amino acids. However, there is a small group of proteins that incorporate a 21st amino acid, Selenocysteine (Sec). These proteins are called selenoproteins and are present in some, but not all, species from the three domains of life. Sec is inserted in selenoproteins in response to the UGA codon, normally a stop codon. A Sec specific tRNA (tRNASec), which only exists in the organisms that synthesize selenoproteins recognizes the UGA codon. tRNASec is not only indispensable for Sec incorporation into selenoproteins, but also for Sec synthesis, since Sec is synthesized on its own tRNA. The structure of tRNASec differs from that of canonical tRNAs, and general tRNA detection methods fail to accurately predict it. We developed Secmarker, a tRNASec specific identification tool based on the characteristic structural features of the tRNASec. Our benchmark shows that Secmarker produces nearly flawless tRNASec predictions. We used Secmarker to scan all currently available genome sequences. The analysis of the highly accurate predictions obtained revealed new insights into the biology of tRNASec.

Dynamic redefinition of the 10 UGAs in human and mouse selenoprotein P (Sepp1) mRNAs to specify selenocysteine instead of termination involves two 3′ UTR structural elements (SECIS) and is regulated by selenium availability. In addition to the previously known human Sepp1 mRNA poly(A) addition site just 3′ of SECIS 2, two further sites were identified with one resulting in 10–25% of the mRNA lacking SECIS 2. To address function, mutant mice were generated with either SECIS 1 or SECIS 2 deleted or with the first UGA substituted with a serine codon. They were fed on either high or selenium-deficient diets. The mutants had very different effects on the proportions of shorter and longer product Sepp1 protein isoforms isolated from plasma, and on viability. Spatially and functionally distinctive effects of the two SECIS elements on UGA decoding were inferred. We also bioinformatically identify two selenoprotein S mRNAs with different 5′ sequences predicted to yield products with different N-termini. These results provide insights into SECIS function and mRNA processing in selenoprotein isoform diversity.

Selenophosphate synthetase (SPS or SelD) provides the active selenium (Se) donor for the synthesis of selenocysteine (Sec), the 21st amino acid in the genetic code. In this chapter we summarize the distribution, phylogeny and function of all SPS genes across the tree of life. SPS is a selenoprotein itself in many prokaryotes (SelD) and eukaryotes (SPS2). As most other selenoproteins, SPS has orthologs with cysteine (Cys) in place of Sec. Although absent in many lineages, selenoproteins and SPS occur in bacteria, archaea and eukaryotes. In prokaryotes, SPS supports additional forms of Se utilization besides Sec, most notably the use of selenouridine in tRNAs. The study of selenophosphate synthetases, while serving a map of Se utilization across all sequenced organisms, also highlighted examples of functional diversification within this family. Within archaea, a few Crenarchaeota species exhibit a SelD-like gene. This is derived from SPS, but probably carries a different function, since it never co-occurs with other Se utilization genes. Within eukaryotes, many metazoan genomes, including humans, carry a paralog called SPS1 in addition to SPS2, which replaces the Sec/Cys site with some other amino acid (i.e., threonine, arginine, glycine, or leucine). Strikingly, SPS1 genes were generated through distinct gene duplication events of SPS2 in several metazoan clades (e.g., vertebrates and insects). Their function is still unknown, but it appears to be different from that of SPS2. Despite independently originated in parallel lineages, SPS1 genes were shown to share a common function. Thus, this function was likely already present in the parental SPS2 gene, driving the gene duplication events by sub-functionalization.

The selenium-dependent glutathione peroxidase (SeGPx) is a well-studied enzyme that detoxifies organic and hydrogen peroxides and provides cells or extracellular fluids with a key antioxidant function. The presence of a SeGPx has not been unequivocally demonstrated in insects. In the present work, we identified the gene and studied the function of a Rhodnius prolixus SeGPx (RpSeGPx). The RpSeGPx mRNA presents the UGA codon that encodes the active site selenocysteine (Sec) and a corresponding Sec insertion sequence (SECIS) in the 3' UTR region. The encoded protein includes a signal peptide, which is consistent with the high levels of GPx enzymatic activity in the insect's hemolymph, and clusters phylogenetically with the extracellular mammalian GPx03. This result contrasts with all other known insect GPxs, which use a cysteine residue instead of Sec and cluster with the mammalian phospholipid hydroperoxide GPx04. RpSeGPx is widely expressed in insect organs, with higher expression levels in the fat body. RNA interference (RNAi) was used to reduce RpSeGPx gene expression and GPx activity in the hemolymph. Adult females were apparently unaffected by RpSeGPx RNAi, whereas first instar nymphs showed a three-day delay in ecdysis. Silencing of RpSeGPx did not alter the gene expression of the antioxidant enzymes catalase, xanthine dehydrogenase and a cysteine-GPx, but it reduced the levels of the dual oxidase and NADPH oxidase 5 transcripts that encode for enzymes releasing extracellular hydrogen peroxide/superoxide. Collectively, our data suggest that RpSeGPx functions in the regulation of extracellular (hemolymph) redox homeostasis of R. prolixus.

Selenoproteins are proteins that incorporate selenocysteine (Sec), a nonstandard amino acid encoded by UGA, normally a stop codon. Sec synthesis requires the enzyme Selenophosphate synthetase (SPS or SelD), conserved in all prokaryotic and eukaryotic genomes encoding selenoproteins. Here, we study the evolutionary history of SPS genes, providing a map of selenoprotein function spanning the whole tree of life. SPS is itself a selenoprotein in many species, although functionally equivalent homologs that replace the Sec site with cysteine (Cys) are common. Many metazoans, however, possess SPS genes with substitutions other than Sec or Cys (collectively referred to as SPS1). Using complementation assays in fly mutants, we show that these genes share a common function, which appears to be distinct from the synthesis of selenophosphate carried out by the Sec- and Cys- SPS genes (termed SPS2), and unrelated to Sec synthesis. We show here that SPS1 genes originated through a number of independent gene duplications from an ancestral metazoan selenoprotein SPS2 gene that most likely already carried the SPS1 function. Thus, in SPS genes, parallel duplications and subsequent convergent subfunctionalization have resulted in the segregation to different loci of functions initially carried by a single gene. This evolutionary history constitutes a remarkable example of emergence and evolution of gene function, which we have been able to trace thanks to the singular features of SPS genes, wherein the amino acid at a single site determines unequivocally protein function and is intertwined to the evolutionary fate of the entire selenoproteome. . R.G. group research was funded by grants BIO2011-26205 from the Spanish Ministry of Science and grant SGR-1430 from the Catalan Government. M.M. received a FPU doctoral fellowship AP2008-04334 from the Spanish Ministry of Education. T.G. group research was funded in part by a grant from the Spanish Ministry of Economy and Competitiveness (BIO2012-37161), a grant from the Qatar National Research Fund (Grant No. NPRP 5-298-3-086), and a grant from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC (Grant Agreement n. ERC-2012-StG-310325). R.G. and T.G. acknowledge support of the Spanish Ministry of Economy and Competitiveness, “Centro de Excelencia Severo Ochoa 2013-2017,” SEV-2012-0208. V.N.G. group research was supported by the National Institutes of Health (NIH) GM061603.

SelenoDB (http://www.selenodb.org) aims to provide high-quality annotations of selenoprotein genes, proteins and SECIS elements. Selenoproteins are proteins that contain the amino acid selenocysteine (Sec) and the first release of the database included annotations for eight species. Since the release of SelenoDB 1.0 many new animal genomes have been sequenced. The annotations of selenoproteins in new genomes usually contain many errors in major databases. For this reason, we have now fully annotated selenoprotein genes in 58 animal genomes. We provide manually curated annotations for human selenoproteins, whereas we use an automatic annotation pipeline to annotate selenoprotein genes in other animal genomes. In addition, we annotate the homologous genes containing cysteine (Cys) instead of Sec. Finally, we have surveyed genetic variation in the annotated genes in humans. We use exon capture and resequencing approaches to identify single-nucleotide polymorphisms in more than 50 human populations around the world. We thus present a detailed view of the genetic divergence of Sec- and Cys-containing genes in animals and their diversity in humans. The addition of these datasets into the second release of the database provides a valuable resource for addressing medical and evolutionary questions in selenium biology. This work was supported by the Plan Nacional and the Instituto Nacional de BIoinformatica (Spain)