Evolution of selenoproteins in the metazoan

Biology and Roles in Diseases of Selenoprotein I Characterized by Ethanolamine Phosphotransferase Activity and Antioxidant Potential

Selenoprotein I (SELENOI) has been demonstrated to be an ethanolamine phosphotransferase (EPT) characterized by a nonselenoenzymatic domain and to be involved in the main synthetic branch of phosphatidylethanolamine (PE) in the endoplasmic reticulum. Therefore, defects of SELENOI may affect the health status through the multiple functions of PE. On the other hand, selenium (Se) is covalently incorporated into SELENOI as selenocysteine (Sec) in its peptide, which forms a Sec-centered domain as in the other members of the selenoprotein family. Unlike other selenoproteins, Sec-containing SELENOI was formed at a later stage of animal evolution, and the high conservation of the structural domain for PE synthesis across a wide range of species suggests the importance of EPT activity in supporting the survival and evolution of organisms. A variety of factors, such as species characteristics (age and sex), diet and nutrition (dietary Se and fat intakes), SELENOI-specific properties (tissue distribution and rank in the selenoproteome), etc., synergistically regulate the expression of SELENOI in a tentatively unclear interaction. The N- and C-terminal domains confer 2 distinct biochemical functions to SELENOI, namely PE regulation and antioxidant potential, which may allow it to be involved in numerous physiological processes, including neurological diseases (especially hereditary spastic paraplegia), T cell activation, tumorigenesis, and adipocyte differentiation. In this review, we summarize advances in the biology and roles of SELENOI, shedding light on the precise regulation of SELENOI expression and PE homeostasis by dietary Se intake and pharmaceutical or transgenic approaches to modulate the corresponding pathological status.

The selenoprotein methionine sulfoxide reductase B1 (MSRB1)

Methionine (Met) can be oxidized to methionine sulfoxide (MetO), which exist as R- and S-diastereomers. Present in all three domains of life, methionine sulfoxide reductases (MSR) are the enzymes that reduce MetO back to Met. Most characterized among them are MSRA and MSRB, which are strictly stereospecific for the S- and R-diastereomers of MetO, respectively. While the majority of MSRs use a catalytic Cys to reduce their substrates, some employ selenocysteine. This is the case of mammalian MSRB1, which was initially discovered as selenoprotein SELR or SELX and later was found to exhibit an MSRB activity. Genomic analyses demonstrated its occurrence in most animal lineages, and biochemical and structural analyses uncovered its catalytic mechanism. The use of transgenic mice and mammalian cell culture revealed its physiological importance in the protection against oxidative stress, maintenance of neuronal cells, cognition, cancer cell proliferation, and the immune response. Coincident with the discovery of Met oxidizing MICAL enzymes, recent findings of MSRB1 regulating the innate immunity response through reversible stereospecific Met-R-oxidation of cytoskeletal actin opened up new avenues for biological importance of MSRB1 and its role in disease. In this review, we discuss the current state of research on MSRB1, compare it with other animal Msrs, and offer a perspective on further understanding of biological functions of this selenoprotein.

Characterization of a new selenoprotein methionine sulfoxide reductase from Haematococcus pluvialis and its antioxidant activity in response to high light intensity, hydrogen peroxide, glyphosate, and cadmium exposure

Selenium incorporates into selenocysteine (Sec) which is a key component of selenoproteins implicated in antioxidant defense and redox homeostasis. Methionine sulfoxide reductases (Msr) play crucial roles in cellular defense against environmental stress. Whereas mammals have the MsrB selenoprotein form, unicellular organisms have MsrA. The Sec residue at the conserved catalytic sites of selenoprotein MsrA confers a metabolic advantage over the non-selenoprotein type MsrA. In the present study, the novel selenoprotein HpMsrA from Haematococcus pluvialis was cloned by the rapid amplification of cDNA ends and transformed into the model green alga Chlamydomonas reinhardtii. Alignment of homologs revealed the presence of the conserved catalytic domain GUFW and showed that the HpMsrA protein comprises Sec (U) at the N-terminus but no recycled Cys at the C-terminus. We studied the response of HpMsrA expression to selenite, high light intensity, hydrogen peroxide, cadmium nitrate, and glyphosate exposure via real-time quantitative PCR and enzyme activity analysis. The results demonstrated that HpMsrA protects cellular proteins against oxidative and environmental stressors. Compared with wild type C. reinhardtii, the transformant exhibited a superior antioxidant ability. The discoveries made herein shed light on the antioxidant physiology and environmental stress resistance mechanisms of the selenoproteins in microalgae. This information may aid in conducting environmental risk assessments of aquatic ecosystems involving microalgae known to respond rapidly and quantitatively to abiotic stress factors promoting excessive reactive oxygen species generation.

Selenoprotein: Potential Player in Redox Regulation in Chlamydomonas reinhardtii

Selenium (Se) is an essential micro-element for many organisms, including Chlamydomonas reinhardtii, and is required in trace amounts. It is obtained from the 21st amino acid selenocysteine (Sec, U), genetically encoded by the UGA codon. Proteins containing Sec are known as selenoproteins. In eukaryotes, selenoproteins are present in animals and algae, whereas fungi and higher plants lack them. The human genome contains 25 selenoproteins, most of which are involved in antioxidant defense activity, redox regulation, and redox signaling. In algae, 42 selenoprotein families were identified using various bioinformatics approaches, out of which C. reinhardtii is known to have 10 selenoprotein genes. However, the role of selenoproteins in Chlamydomonas is yet to be reported. Chlamydomonas selenoproteins contain conserved domains such as CVNVGC and GCUG, in the case of thioredoxin reductase, and CXXU in other selenoproteins. Interestingly, Sec amino acid residue is present in a catalytically active domain in Chlamydomonas selenoproteins, similar to human selenoproteins. Based on catalytical active sites and conserved domains present in Chlamydomonas selenoproteins, we suggest that Chlamydomonas selenoproteins could have a role in redox regulation and defense by acting as antioxidants in various physiological conditions.

The Role of Selenium Nanoparticles in Agriculture and Food Technology

Selenium (Se) is an essential micronutrient for diverse organisms such as mammals, bacteria, some insects and nematodes, archaea, and algae, as it is involved in a large number of physiological and metabolic processes and is part of approximately 25 selenoproteins in mammals. In plants, Se has no essential metabolic role, high concentrations of inorganic Se can lead to the formation of Se-amino acids, and its incorporation into selenoproteins can generate toxicity. Conversely, low doses of Se can trigger a variety of beneficial effects as an antioxidant, antimicrobial, or stress-modulating agent without being an essential element. Therefore, Se can generate toxicity depending on the dose and the chemical form in which it is supplied. Selenium nanoparticles (SeNPs) have emerged as an approach to reduce this negative effect and improve its biological properties. In turn, SeNPs have a wide range of potential advantages, making them an alternative for areas such as agriculture and food technology. This review focuses on the use of SeNPs and their different applications as antimicrobial agents, growth promoters, crop biofortification, and nutraceuticals in agriculture. In addition, the utilization of SeNPs in the generation of packaging with antioxidant and antimicrobial traits and Se enrichment of animal source foods for human consumption as part of food technology is addressed. Additionally, possible action mechanisms and potential adverse effects are discussed. The concentration, size, and synthesis method of SeNPs are determining factors of their biological properties.

Larval metamorphosis is inhibited by methimazole and propylthiouracil that reveals possible hormonal action in the mussel Mytilus coruscus

Larval metamorphosis in bivalves is a key event for the larva-to-juvenile transformation. Previously we have identified a thyroid hormone receptor (TR) gene that is crucial for larvae to acquire “competence” for the metamorphic transition in the mussel Mytilus courscus (Mc). The mechanisms of thyroid signaling in bivalves are still largely unknown. In the present study, we molecularly characterized the full-length of two iodothyronine deiodinase genes (McDx and McDy). Phylogenetic analysis revealed that deiodinases of molluscs (McDy, CgDx and CgDy) and vertebrates (D2 and D3) shared a node representing an immediate common ancestor, which resembled vertebrates D1 and might suggest that McDy acquired specialized function from vertebrates D1. Anti-thyroid compounds, methimazole (MMI) and propylthiouracil (PTU), were used to investigate their effects on larval metamorphosis and juvenile development in M. coruscus. Both MMI and PTU significantly reduced larval metamorphosis in response to the metamorphosis inducer epinephrine. MMI led to shell growth retardation in a concentration-dependent manner in juveniles of M. coruscus after 4 weeks of exposure, whereas PTU had no effect on juvenile growth. It is hypothesized that exposure to MMI and PTU reduced the ability of pediveliger larvae for the metamorphic transition to respond to the inducer. The effect of MMI and PTU on larval metamorphosis and development is most likely through a hormonal signal in the mussel M. coruscus, with the implications for exploring the origins and evolution of metamorphosis.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

Programmed Deviations of Ribosomes From Standard Decoding in Archaea

Genetic code decoding, initially considered to be universal and immutable, is now known to be flexible. In fact, in specific genes, ribosomes deviate from the standard translational rules in a programmed way, a phenomenon globally termed recoding. Translational recoding, which has been found in all domains of life, includes a group of events occurring during gene translation, namely stop codon readthrough, programmed ± 1 frameshifting, and ribosome bypassing. These events regulate protein expression at translational level and their mechanisms are well known and characterized in viruses, bacteria and eukaryotes. In this review we summarize the current state-of-the-art of recoding in the third domain of life. In Archaea, it was demonstrated and extensively studied that translational recoding regulates the decoding of the 21st and the 22nd amino acids selenocysteine and pyrrolysine, respectively, and only one case of programmed -1 frameshifting has been reported so far in Saccharolobus solfataricus P2. However, further putative events of translational recoding have been hypothesized in other archaeal species, but not extensively studied and confirmed yet. Although this phenomenon could have some implication for the physiology and adaptation of life in extreme environments, this field is still underexplored and genes whose expression could be regulated by recoding are still poorly characterized. The study of these recoding episodes in Archaea is urgently needed.

Deiodinases: How Nonmammalian Research Helped Shape Our Present View

Iodothyronine deiodinases are enzymes capable of activating and inactivating thyroid hormones (THs) and have an important role in regulating TH action in tissues throughout the body. Three types of deiodinases (D1, D2, and D3) were originally defined based on their biochemical characteristics. Cloning of the first complementary DNAs in the 1990s (Dio1 in rat and dio2 and dio3 in frog) allowed to confirm the existence of 3 distinct enzymes. Over the years, increasing genomic information revealed that deiodinases are present in all chordates, vertebrates, and nonvertebrates and that they can even be found in some mollusks and annelids, pointing to an ancient origin. Research in nonmammalian models has substantially broadened our understanding of deiodinases. In relation to their structure, we discovered for instance that biochemical properties such as inhibition by 6-propyl-2-thiouracil, stimulation by dithiothreitol, and temperature optimum are subject to variation. Data from fish, amphibians, and birds were key in shifting our view on the relative importance of activating and inactivating deiodination pathways and in showing the impact of D2 and D3 not only in local but also whole body T3 availability. They also led to the discovery of new local functions such as the acute reciprocal changes in D2 and D3 in hypothalamic tanycytes upon photostimulation, involved in seasonal rhythmicity. With the present possibilities for rapid and precise gene silencing in any species of interest, comparative research will certainly further contribute to a better understanding of the importance of deiodinases for adequate TH action, also in humans.

Virtual 2-D map of the fungal proteome

The molecular weight and isoelectric point (pI) of the proteins plays important role in the cell. Depending upon the shape, size, and charge, protein provides its functional role in different parts of the cell. Therefore, understanding to the knowledge of their molecular weight and charges is (pI) is very important. Therefore, we conducted a proteome-wide analysis of protein sequences of 689 fungal species (7.15 million protein sequences) and construct a virtual 2-D map of the fungal proteome. The analysis of the constructed map revealed the presence of a bimodal distribution of fungal proteomes. The molecular mass of individual fungal proteins ranged from 0.202 to 2546.166 kDa and the predicted isoelectric point (pI) ranged from 1.85 to 13.759 while average molecular weight of fungal proteome was 50.98 kDa. A non-ribosomal peptide synthase (RFU80400.1) found in Trichoderma arundinaceum was identified as the largest protein in the fungal kingdom. The collective fungal proteome is dominated by the presence of acidic rather than basic pI proteins and Leu is the most abundant amino acid while Cys is the least abundant amino acid. Aspergillus ustus encodes the highest percentage (76.62%) of acidic pI proteins while Nosema ceranae was found to encode the highest percentage (66.15%) of basic pI proteins. Selenocysteine and pyrrolysine amino acids were not found in any of the analysed fungal proteomes. Although the molecular weight and pI of the protein are of enormous important to understand their functional roles, the amino acid compositions of the fungal protein will enable us to understand the synonymous codon usage in the fungal kingdom. The small peptides identified during the study can provide additional biotechnological implication.

Demethylation of Methylmercury in Bird, Fish, and Earthworm

Toxicity of methylmercury (MeHg) to wildlife and humans results from its binding to cysteine residues of proteins, forming MeHg-cysteinate (MeHgCys) complexes that hinder biological functions. MeHgCys complexes can be detoxified in vivo, yet how this occurs is unknown. We report that MeHgCys complexes are transformed into selenocysteinate [Hg(Sec)₄] complexes in multiple animals from two phyla (a waterbird, freshwater fish, and earthworms) sampled in different geographical areas and contaminated by different Hg sources. In addition, high energy-resolution X-ray absorption spectroscopy (HR-XANES) and chromatography-inductively coupled plasma mass spectrometry of the waterbird liver support the binding of Hg(Sec)₄ to selenoprotein P and biomineralization of Hg(Sec)₄ to chemically inert nanoparticulate mercury selenide (HgSe). The results provide a foundation for understanding mercury detoxification in higher organisms and suggest that the identified MeHgCys to Hg(Sec)₄ demethylation pathway is common in nature.

Common modifications of selenocysteine in selenoproteins

Selenocysteine (Sec), the sulfur-to-selenium substituted variant of cysteine (Cys), is the defining entity of selenoproteins. These are naturally expressed in many diverse organisms and constitute a unique class of proteins. As a result of the physicochemical characteristics of selenium when compared with sulfur, Sec is typically more reactive than Cys while participating in similar reactions, and there are also some qualitative differences in the reactivities between the two amino acids. This minireview discusses the types of modifications of Sec in selenoproteins that have thus far been experimentally validated. These modifications include direct covalent binding through the Se atom of Sec to other chalcogen atoms (S, O and Se) as present in redox active molecular motifs, derivatization of Sec via the direct covalent binding to non-chalcogen elements (Ni, Mb, N, Au and C), and the loss of Se from Sec resulting in formation of dehydroalanine. To understand the nature of these Sec modifications is crucial for an understanding of selenoprotein reactivities in biological, physiological and pathophysiological contexts.

Construction of anti-codon table of the plant kingdom and evolution of tRNA selenocysteine (tRNASec)

Background

The tRNAs act as a bridge between the coding mRNA and incoming amino acids during protein translation. The anti-codon of tRNA recognizes the codon of the mRNA and deliver the amino acid into the protein translation chain. However, we did not know about the exact abundance of anti-codons in the genome and whether the frequency of abundance remains same across the plant lineage or not.

Results

Therefore, we analysed the tRNAnome of 128 plant species and reported an anti-codon table of the plant kingdom. We found that CAU anti-codon of tRNA^Met has highest (5.039%) whereas GCG anti-codon of tRNA^Arg has lowest (0.004%) abundance. However, when we compared the anti-codon frequencies according to the tRNA isotypes, we found tRNA^Leu (7.808%) has highest abundance followed by tRNA^Ser (7.668%) and tRNA^Gly (7.523%). Similarly, suppressor tRNA (0.036%) has lowest abundance followed by tRNA^Sec (0.066%) and tRNA^His (2.109). The genome of Ipomoea nil, Papaver somniferum, and Zea mays encoded the highest number of anti-codons (isoacceptor) at 59 each whereas the genome of Ostreococcus tauri was found to encode only 18 isoacceptors. The tRNA^Sec genes undergone losses more frequently than duplication and we found that tRNA^Sec showed anti-codon switch during the course of evolution.

Conclusion

The anti-codon table of the plant tRNA will enable us to understand the synonymous codon usage of the plant kingdom and can be very helpful to understand which codon is preferred over other during the translation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07216-3.

Microbial selenium metabolism: a brief history, biogeochemistry and ecophysiology

Selenium is an essential trace element for organisms from all three domains of life. Microorganisms, in particular, mediate reductive transformations of selenium that govern the element's mobility and bioavailability in terrestrial and aquatic environments. Selenium metabolism is not just ubiquitous but an ancient feature of life likely extending back to the universal common ancestor of all cellular lineages. As with the sulfur biogeochemical cycle, reductive transformations of selenium serve two metabolic functions: assimilation into macromolecules and dissimilatory reduction during anaerobic respiration. This review begins with a historical overview of how research in both aspects of selenium metabolism has developed. We then provide an overview of the global selenium biogeochemical cycle, emphasizing the central role of microorganisms in the cycle. This serves as a basis for a robust discussion of current models for the evolution of the selenium biogeochemical cycle over geologic time, and how knowledge of the evolution and ecophysiology of selenium metabolism can enrich and refine these models. We conclude with a discussion of the ecophysiological function of selenium-respiring prokaryotes within the cycle, and the tantalizing possibility of oxidative selenium transformations during chemolithoautotrophic growth.

The algal selenoproteomes

Background

Selenium is an essential trace element, and selenocysteine (Sec, U) is its predominant form in vivo. Proteins that contain Sec are selenoproteins, whose special structural features include not only the TGA codon encoding Sec but also the SECIS element in mRNA and the conservation of the Sec-flanking region. These unique features have led to the development of a series of bioinformatics methods to predict and research selenoprotein genes. There have been some studies and reports on the evolution and distribution of selenoprotein genes in prokaryotes and multicellular eukaryotes, but the systematic analysis of single-cell eukaryotes, especially algae, has been very limited.

Results

In this study, we predicted selenoprotein genes in 137 species of algae by using a program we previously developed. More than 1000 selenoprotein genes were obtained. A database website was built to record these algae selenoprotein genes (www.selenoprotein.com). These genes belong to 42 selenoprotein families, including three novel selenoprotein gene families.

Conclusions

This study reveals the primordial state of the eukaryotic selenoproteome. It is an important clue to explore the significance of selenium for primordial eukaryotes and to determine the complete evolutionary spectrum of selenoproteins in all life forms.

Bioinformatics of Selenoproteins

Significance: Bioinformatics has brought important insights into the field of selenium research. The progress made in the development of computational tools in the last two decades, coordinated with growing genome resources, provided new opportunities to study selenoproteins. The present review discusses existing tools for selenoprotein gene finding and other bioinformatic approaches to study the biology of selenium. Recent Advances: The availability of complete selenoproteomes allowed assessing a global distribution of the use of selenocysteine (Sec) across the tree of life, as well as studying the evolution of selenoproteins and their biosynthetic pathway. Beyond gene identification and characterization, human genetic variants in selenoprotein genes were used to examine adaptations to selenium levels in diverse human populations and to estimate selective constraints against gene loss. Critical Issues: The synthesis of selenoproteins is essential for development in mice. In humans, several mutations in selenoprotein genes have been linked to rare congenital disorders. And yet, the mechanism of Sec insertion and the regulation of selenoprotein synthesis in mammalian cells are not completely understood. Future Directions: Omics technologies offer new possibilities to study selenoproteins and mechanisms of Sec incorporation in cells, tissues, and organisms.

Bioinformatics of Metalloproteins and Metalloproteomes

Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.

Systems Biology of Selenium and Complex Disease

Selenium is an essential trace element for maintenance of overall health, whose deficiency and dyshomeostasis have been linked to a variety of diseases and disorders. The majority of previous researches focused on characterization of genes encoding selenoproteins or proteins involved in selenium metabolism as well as their functions. Many studies in humans also investigated the relationship between selenium and complex diseases, but their results have been inconsistent. In recent years, systems biology and "-omics" approaches have been widely used to study complex and global variations of selenium metabolism and function in physiological and different pathological conditions. The present paper reviews recent progress in large-scale and systematic analyses of the relationship between selenium status or selenoproteins and several complex diseases, mainly including population-based cohort studies and meta-analyses, genetic association studies, and some other omics-based studies. Advances in ionomics and its application in studying the interaction between selenium and other trace elements in human health and diseases are also discussed.

The evolutionary road to invertebrate thyroid hormone signaling: Perspectives for endocrine disruption processes

Thyroid hormones (THs) are the only iodine-containing hormones that play fundamental roles in chordates and non-chordates. The chemical nature, mode of action and the synthesis of THs are well established in mammals and other vertebrates. Although thyroid-like hormones have been detected in protostomes and non-chordate deuterostomes, TH signaling is poorly understood as compared to vertebrates, particularly in protostomes. Therefore, the central objective of this article is to review TH system components and TH-induced effects in non-vertebrate chordates, non-chordate deuterostomes and protostomes based on available genomes and functional information. To accomplish this task, we integrate here the available knowledge on the THs signaling across non-vertebrate chordates, non-chordate deuterostomes and protostomes by considering studies encompassing TH system components and physiological actions of THs. We also address the possible interactions of thyroid disrupting chemicals and their effects in protostomes and non-chordate deuterostomes. Finally, the perspectives on current and future challenges are discussed.

Processive Recoding and Metazoan Evolution of Selenoprotein P: Up to 132 UGAs in Molluscs

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.

Utilization of selenocysteine in early-branching fungal phyla

Selenoproteins are a diverse group of proteins containing selenocysteine (Sec)-the twenty-first amino acid-incorporated during translation via a unique recoding mechanism^1,2. Selenoproteins fulfil essential roles in many organisms¹, yet are not ubiquitous across the tree of life^3-7. In particular, fungi were deemed devoid of selenoproteins^4,5,8. However, we show here that Sec is utilized by nine species belonging to diverse early-branching fungal phyla, as evidenced by the genomic presence of both Sec machinery and selenoproteins. Most fungal selenoproteins lack consensus Sec recoding signals (SECIS elements⁹) but exhibit other RNA structures, suggesting altered mechanisms of Sec insertion in fungi. Phylogenetic analyses support a scenario of vertical inheritance of the Sec trait within eukaryotes and fungi. Sec was then lost in numerous independent events in various fungal lineages. Notably, Sec was lost at the base of Dikarya, resulting in the absence of selenoproteins in Saccharomyces cerevisiae and other well-studied fungi. Our results indicate that, despite scattered occurrence, selenoproteins are found in all kingdoms of life.

Comparative genomics and metagenomics of the metallomes

Recent achievements and advances in comparative genomic and metagenomic analyses of trace metals were reviewed and discussed.

SelGenAmic: An Algorithm for Selenoprotein Gene Assembly

Computational methods for identifying selenoproteins have been developed rapidly in recent years. However, it is still difficult to identify the open reading frame (ORF) of eukaryotic selenoprotein gene, because the TGA codon for a selenocysteine (Sec) residue in the active center of selenoprotein is traditionally a terminal signal of protein translation. A gene assembly algorithm SelGenAmic has been constructed and presented in this chapter for identifying selenoprotein genes from eukaryotic genomes. A method based on this algorithm was developed to build an optimal TGA-containing-ORF for each TGA in a genome, followed by protein similarity analysis through conserved sequence alignments to screen out selenoprotein genes from these ORFs. This method improved the sensitivity of detecting selenoproteins from a genome due to the design that all TGAs in the genome were investigated for its possibility of decoding as a Sec residue. The method based on the SelGenAmic algorithm is capable of identifying eukaryotic selenoprotein genes from their genomes.

SECISearch3 and Seblastian: In-Silico Tools to Predict SECIS Elements and Selenoproteins

The computational identification of selenoprotein genes is complicated by the dual meaning of the UGA codon as stop and selenocysteine. SECIS elements are RNA structures essential for selenocysteine incorporation, which have been used as markers for selenoprotein genes in many bioinformatics studies. The most widely used tool for eukaryotic SECIS finding has been recently improved to its third generation, SECISearch3. This program is also a component of Seblastian, a pipeline for the identification of selenoprotein genes that employs SECIS finding as the first step. This chapter constitutes a practical guide to use SECISearch3 and Seblastian, which can be run via webservers at http://seblastian.crg.eu / or http://gladyshevlab.org/SelenoproteinPredictionServer/ .

Genome-wide transcriptomics of aging in the rotifer Brachionus manjavacas, an emerging model system

Background

Understanding gene expression changes over lifespan in diverse animal species will lead to insights to conserved processes in the biology of aging and allow development of interventions to improve health. Rotifers are small aquatic invertebrates that have been used in aging studies for nearly 100 years and are now re-emerging as a modern model system. To provide a baseline to evaluate genetic responses to interventions that change health throughout lifespan and a framework for new hypotheses about the molecular genetic mechanisms of aging, we examined the transcriptome of an asexual female lineage of the rotifer Brachionus manjavacas at five life stages: eggs, neonates, and early-, late-, and post-reproductive adults.

Results

There are widespread shifts in gene expression over the lifespan of B. manjavacas; the largest change occurs between neonates and early reproductive adults and is characterized by down-regulation of developmental genes and up-regulation of genes involved in reproduction. The expression profile of post-reproductive adults was distinct from that of other life stages. While few genes were significantly differentially expressed in the late- to post-reproductive transition, gene set enrichment analysis revealed multiple down-regulated pathways in metabolism, maintenance and repair, and proteostasis, united by genes involved in mitochondrial function and oxidative phosphorylation.

Conclusions

This study provides the first examination of changes in gene expression over lifespan in rotifers. We detected differential expression of many genes with human orthologs that are absent in Drosophila and C. elegans, highlighting the potential of the rotifer model in aging studies. Our findings suggest that small but coordinated changes in expression of many genes in pathways that integrate diverse functions drive the aging process. The observation of simultaneous declines in expression of genes in multiple pathways may have consequences for health and longevity not detected by single- or multi-gene knockdown in otherwise healthy animals. Investigation of subtle but genome-wide change in these pathways during aging is an important area for future study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3540-x) contains supplementary material, which is available to authorized users.

Computational identification of the selenocysteine tRNA (tRNASec) in genomes

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 (tRNA^Sec) 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 tRNA^Sec 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 tRNA^Sec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNA^Sec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNA^Sec 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 tRNA^Sec 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 tRNA^Sec, 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.

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 (tRNA^Sec), which only exists in the organisms that synthesize selenoproteins recognizes the UGA codon. tRNA^Sec 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 tRNA^Sec differs from that of canonical tRNAs, and general tRNA detection methods fail to accurately predict it. We developed Secmarker, a tRNA^Sec specific identification tool based on the characteristic structural features of the tRNA^Sec. Our benchmark shows that Secmarker produces nearly flawless tRNA^Sec 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 tRNA^Sec.

Lokiarchaeota Marks the Transition between the Archaeal and Eukaryotic Selenocysteine Encoding Systems

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.

Why Nature Chose Selenium

The authors were asked by the Editors of ACS Chemical Biology to write an article titled "Why Nature Chose Selenium" for the occasion of the upcoming bicentennial of the discovery of selenium by the Swedish chemist Jöns Jacob Berzelius in 1817 and styled after the famous work of Frank Westheimer on the biological chemistry of phosphate [Westheimer, F. H. (1987) Why Nature Chose Phosphates, Science 235, 1173-1178]. This work gives a history of the important discoveries of the biological processes that selenium participates in, and a point-by-point comparison of the chemistry of selenium with the atom it replaces in biology, sulfur. This analysis shows that redox chemistry is the largest chemical difference between the two chalcogens. This difference is very large for both one-electron and two-electron redox reactions. Much of this difference is due to the inability of selenium to form π bonds of all types. The outer valence electrons of selenium are also more loosely held than those of sulfur. As a result, selenium is a better nucleophile and will react with reactive oxygen species faster than sulfur, but the resulting lack of π-bond character in the Se-O bond means that the Se-oxide can be much more readily reduced in comparison to S-oxides. The combination of these properties means that replacement of sulfur with selenium in nature results in a selenium-containing biomolecule that resists permanent oxidation. Multiple examples of this gain of function behavior from the literature are discussed.

The selenium content of SEPP1 versus selenium requirements in vertebrates

Selenoprotein P (SEPP1) distributes selenium (Se) throughout the body via the circulatory system. For vertebrates, the Se content of SEPP1 varies from 7 to 18 Se atoms depending on the species, but the reason for this variation remains unclear. Herein we provide evidence that vertebrate SEPP1 Sec content correlates positively with Se requirements. As the Se content of full length SEPP1 is genetically determined, this presents a unique case where a nutrient requirement can be predicted based on genomic sequence information.

Evolution of selenophosphate synthetases: emergence and relocation of function through independent duplications and recurrent subfunctionalization

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.

Yeast as a model system to study metabolic impact of selenium compounds

Inorganic Se forms such as selenate or selenite (the two more abundant forms in nature) can be toxic in Saccharomyces cerevisiae cells, which constitute an adequate model to study such toxicity at the molecular level and the functions participating in protection against Se compounds. Those Se forms enter the yeast cell through other oxyanion transporters. Once inside the cell, inorganic Se forms may be converted into selenide through a reductive pathway that in physiological conditions involves reduced glutathione with its consequent oxidation into diglutathione and alteration of the cellular redox buffering capacity. Selenide can subsequently be converted by molecular oxygen into elemental Se, with production of superoxide anions and other reactive oxygen species. Overall, these events result in DNA damage and dose-dependent reversible or irreversible protein oxidation, although additional oxidation of other cellular macromolecules cannot be discarded. Stress-adaptation pathways are essential for efficient Se detoxification, while activation of DNA damage checkpoint and repair pathways protects against Se-mediated genotoxicity. We propose that yeast may be used to improve our knowledge on the impact of Se on metal homeostasis, the identification of Se-targets at the DNA and protein levels, and to gain more insights into the mechanism of Se-mediated apoptosis.

Comparative selenoproteome analysis reveals a reduced utilization of selenium in parasitic platyhelminthes

Background. The selenocysteine(Sec)-containing proteins, selenoproteins, are an important group of proteins present in all three kingdoms of life. Although the selenoproteomes of many organisms have been analyzed, systematic studies on selenoproteins in platyhelminthes are still lacking. Moreover, comparison of selenoproteomes between free-living and parasitic animals is rarely studied.

Results. In this study, three representative organisms (Schmidtea mediterranea, Schistosoma japonicum and Taenia solium) were selected for comparative analysis of selenoproteomes in Platyhelminthes. Using a SelGenAmic-based selenoprotein prediction algorithm, a total of 37 selenoprotein genes were identified in these organisms. The size of selenoproteomes and selenoprotein families were found to be associated with different lifestyles: free-living organisms have larger selenoproteome whereas parasitic lifestyle corresponds to reduced selenoproteomes. Five selenoproteins, SelT, Sel15, GPx, SPS2 and TR, were found to be present in all examined platyhelminthes as well as almost all sequenced animals, suggesting their essential role in metazoans. Finally, a new splicing form of SelW that lacked the first exon was found to be present in S. japonicum.

Conclusions. Our data provide a first glance into the selenoproteomes of organisms in the phylum Platyhelminthes and may help understand function and evolutionary dynamics of selenium utilization in diversified metazoans.

SECISearch3 and Seblastian: new tools for prediction of SECIS elements and selenoproteins

Selenoproteins are proteins containing an uncommon amino acid selenocysteine (Sec). Sec is inserted by a specific translational machinery that recognizes a stem-loop structure, the SECIS element, at the 3′ UTR of selenoprotein genes and recodes a UGA codon within the coding sequence. As UGA is normally a translational stop signal, selenoproteins are generally misannotated and designated tools have to be developed for this class of proteins. Here, we present two new computational methods for selenoprotein identification and analysis, which we provide publicly through the web servers at http://gladyshevlab.org/SelenoproteinPredictionServer or http://seblastian.crg.es. SECISearch3 replaces its predecessor SECISearch as a tool for prediction of eukaryotic SECIS elements. Seblastian is a new method for selenoprotein gene detection that uses SECISearch3 and then predicts selenoprotein sequences encoded upstream of SECIS elements. Seblastian is able to both identify known selenoproteins and predict new selenoproteins. By applying these tools to diverse eukaryotic genomes, we provide a ranked list of newly predicted selenoproteins together with their annotated cysteine-containing homologues. An analysis of a representative candidate belonging to the AhpC family shows how the use of Sec in this protein evolved in bacterial and eukaryotic lineages.

Selenoprotein-transgenic Chlamydomonas reinhardtii

Selenium (Se) deficiency is associated with the occurrence of many diseases. However, excessive Se supplementation, especially with inorganic Se, can result in toxicity. Selenoproteins are the major forms of Se in vivo to exert its biological function. Expression of those selenoproteins, especially with the application of a newly developed system, is thus very important for studying the mechanism of Se in nutrition. The use of Chlamydomonas reinhardtii (C. reinhardtii) as a biological vector to express an heterogeneous protein is still at the initial stages of development. In order to investigate the possibility of using this system to express selenoproteins, human 15-KDa selenoprotein (Sep15), a small but widely distributed selenoprotein in mammals, was chosen for the expression platform test. Apart from the wild-type human Sep15 gene fragment, two Sep15 recombinants were constructed containing Sep15 open reading frame (ORF) and the selenocysteine insertion sequence (SECIS) element from either human Sep15 or C. reinhardtii selenoprotein W1, a highly expressed selenoprotein in this alga. Those Sep15-containing plasmids were transformed into C. reinhardtii CC-849 cells. Results showed that Sep15 fragments were successfully inserted into the nuclear genome and expressed Sep15 protein in the cells. The transgenic and wild-type algae demonstrated similar growth curves in low Se culture medium. To our knowledge, this is the first report on expressing human selenoprotein in green alga.