Relaxation of selective constraints causes independent selenoprotein extinction in insect genomes

Background

Selenoproteins are a diverse family of proteins notable for the presence of the 21st amino acid, selenocysteine. Until very recently, all metazoan genomes investigated encoded selenoproteins, and these proteins had therefore been believed to be essential for animal life. Challenging this assumption, recent comparative analyses of insect genomes have revealed that some insect genomes appear to have lost selenoprotein genes.

Methodology/Principal Findings

In this paper we investigate in detail the fate of selenoproteins, and that of selenoprotein factors, in all available arthropod genomes. We use a variety of in silico comparative genomics approaches to look for known selenoprotein genes and factors involved in selenoprotein biosynthesis. We have found that five insect species have completely lost the ability to encode selenoproteins and that selenoprotein loss in these species, although so far confined to the Endopterygota infraclass, cannot be attributed to a single evolutionary event, but rather to multiple, independent events. Loss of selenoproteins and selenoprotein factors is usually coupled to the deletion of the entire no-longer functional genomic region, rather than to sequence degradation and consequent pseudogenisation. Such dynamics of gene extinction are consistent with the high rate of genome rearrangements observed in Drosophila. We have also found that, while many selenoprotein factors are concomitantly lost with the selenoproteins, others are present and conserved in all investigated genomes, irrespective of whether they code for selenoproteins or not, suggesting that they are involved in additional, non-selenoprotein related functions.

Conclusions/Significance

Selenoproteins have been independently lost in several insect species, possibly as a consequence of the relaxation in insects of the selective constraints acting across metazoans to maintain selenoproteins. The dispensability of selenoproteins in insects may be related to the fundamental differences in antioxidant defense between these animals and other metazoans.     

The prokaryotic selenoproteome

In the genetic code, the UGA codon has a dual function as it encodes selenocysteine (Sec) and serves as a stop signal. However, only the translation terminator function is used in gene annotation programs, resulting in misannotation of selenoprotein genes. Here, we applied two independent bioinformatics approaches to characterize a selenoprotein set in prokaryotic genomes. One method searched for selenoprotein genes by identifying RNA stem-loop structures, selenocysteine insertion sequence elements; the second approach identified Sec/Cys pairs in homologous sequences. These analyses identified all or almost all selenoproteins in completely sequenced bacterial and archaeal genomes and provided a view on the distribution and composition of prokaryotic selenoproteomes. In addition, lineage-specific and core selenoproteins were detected, which provided insights into the mechanisms of selenoprotein evolution. Characterization of selenoproteomes allows interpretation of other UGA codons in completed genomes of prokaryotes as terminators, addressing the UGA dual-function problem.     

The microbial selenoproteome of the Sargasso Sea

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.     

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

MOTIVATION: Incorporation of selenocysteine (Sec) into proteins in response to UGA codons requires a cis-acting RNA structure, Sec insertion sequence (SECIS) element. Whereas SECIS elements in Escherichia coli are well characterized, a bacterial SECIS consensus structure is lacking. RESULTS: We developed a bacterial SECIS consensus model, the key feature of which is a conserved guanosine in a small apical loop of the properly positioned structure. This consensus was used to build a computational tool, bSECISearch, for detection of bacterial SECIS elements and selenoprotein genes in sequence databases. The program identified 96.5% of known selenoprotein genes in completely sequenced bacterial genomes and predicted several new selenoprotein genes. Further analysis revealed that the size of bacterial selenoproteomes varied from 1 to 11 selenoproteins. Formate dehydrogenase was present in most selenoproteomes, often as the only selenoprotein family, whereas the occurrence of other selenoproteins was limited. The availability of the bacterial SECIS consensus and the tool for identification of these structures should help in correct annotation of selenoprotein genes and characterization of bacterial selenoproteomes.     

Nematode selenoproteome: the use of the selenocysteine insertion system to decode one codon in an animal genome?

Selenocysteine (Sec) is co-translationally inserted into selenoproteins in response to codon UGA with the help of the selenocysteine insertion sequence (SECIS) element. The number of selenoproteins in animals varies, with humans having 25 and mice having 24 selenoproteins. To date, however, only one selenoprotein, thioredoxin reductase, has been detected in Caenorhabditis elegans, and this enzyme contains only one Sec. Here, we characterize the selenoproteomes of C.elegans and Caenorhabditis briggsae with three independent algorithms, one searching for pairs of homologous nematode SECIS elements, another searching for Cys- or Sec-containing homologs of potential nematode selenoprotein genes and the third identifying Sec-containing homologs of annotated nematode proteins. These methods suggest that thioredoxin reductase is the only Sec-containing protein in the C.elegans and C.briggsae genomes. In contrast, we identified additional selenoproteins in other nematodes. Assuming that Sec insertion mechanisms are conserved between nematodes and other eukaryotes, the data suggest that nematode selenoproteomes were reduced during evolution, and that in an extreme reduction case Sec insertion systems probably decode only a single UGA codon in C.elegans and C.briggsae genomes. In addition, all detected genes had a rare form of SECIS element containing a guanosine in place of a conserved adenosine present in most other SECIS structures, suggesting that in organisms with small selenoproteomes SECIS elements may change rapidly.     

Evolution of selenoproteins in the metazoan

ABSTRACT: BACKGROUND: The selenocysteine (Sec) containing proteins, selenoproteins, are an important group of proteins present throughout all 3 kingdoms of life. With the rapid progression of selenoprotein research in the post-genomic era, application of bioinformatics methods to the identification of selenoproteins in newly sequenced species has become increasingly important. Although selenoproteins in human and other vertebrates have been investigated, studies of primitive invertebrate selenoproteomes are rarely reported outside of insects and nematodes. Result A more integrated view of selenoprotein evolution was constructed using several representative species from different evolutionary eras. Using a SelGenAmic-based selenoprotein identification method, 178 selenoprotein genes were identified in 6 invertebrates: Amphimedon queenslandica, Trichoplax adhaerens, Nematostella vectensis, Lottia gigantean, Capitella teleta, and Branchiostoma floridae. Amphioxus was found to have the most abundant and variant selenoproteins of any animal currently characterized, including a special selenoprotein P (SelP) possessing 3 repeated Trx-like domains and Sec residues in the N-terminal and 2 Sec residues in the C-terminal. This gene structure suggests the existence of two different strategies for extension of Sec numbers in SelP for the preservation and transportation of selenium. In addition, novel eukaryotic AphC-like selenoproteins were identified in sponges. CONCLUSION: Comparison of various animal species suggests that even the most primitive animals possess a selenoproteome range and variety similar to humans. During evolutionary history, only a few new selenoproteins have emerged and few were lost. Furthermore, the massive loss of selenoproteins in nematodes and insects likely occurred independently in isolated partial evolutionary branches.     

Comparative analysis of selenocysteine machinery and selenoproteome gene expression in mouse brain identifies neurons as key functional sites of selenium in mammals

Although dietary selenium (Se) deficiency results in phenotypes associated with selenoprotein depletion in various organs, the brain is protected from Se loss. To address the basis for the critical role of Se in brain function, we carried out comparative gene expression analyses for the complete selenoproteome and associated biosynthetic factors. Using the Allen Brain Atlas, we evaluated 159 regions of adult mouse brain and provided experimental analyses of selected selenoproteins. All 24 selenoprotein mRNAs were expressed in the mouse brain. Most strikingly, neurons in olfactory bulb, hippocampus, cerebral cortex, and cerebellar cortex were exceptionally rich in selenoprotein gene expression, in particular in GPx4, SelK, SelM, SelW, and Sep15. Over half of the selenoprotein genes were also expressed in the choroid plexus. A unique expression pattern was observed for one of the highly expressed selenoprotein genes, SelP, which we suggest to provide neurons with Se. Cluster analysis of the expression data linked certain selenoproteins and selenocysteine machinery genes and suggested functional linkages among selenoproteins, such as that between SelM and Sep15. Overall, this study suggests that the main functions of selenium in mammals are confined to certain neurons in the brain.     

Selenoproteins and selenium status in bone physiology and pathology

Background

Emerging evidence supports the view that selenoproteins are essential for maintaining bone health.

Scope of review

The current state of knowledge concerning selenoproteins and Se status in bone physiology and pathology is summarized.

Major conclusions

Antioxidant selenoproteins including glutathione peroxidase (GPx) and thioredoxin reductase (TrxR), as a whole, play a pivotal role in maintaining bone homeostasis and protecting against bone loss. GPx1, a major antioxidant enzyme in osteoclasts, is up-regulated by estrogen, an endogenous inhibitor of osteoclastogenesis. TrxR1 is an immediate early gene in response to 1α,25-dihydroxyvitamin D3, an osteoblastic differentiation agent. The combination of 1α,25-dihydroxyvitamin D3 and Se generates a synergistic elevation of TrxR activity in Se-deficient osteoblasts. Of particular concern, pleiotropic TrxR1 is implicated in promoting NFκB activation. Coincidentally, TrxR inhibitors such as curcumin and gold compounds exhibit potent osteoclastogenesis inhibitory activity. Studies in patients with the mutations of selenocysteine insertion sequence-binding protein 2, a key trans-acting factor for the co-translational insertion of selenocysteine into selenoproteins have clearly established a causal link of selenoproteins in bone development. Se transport to bone relies on selenoprotein P. Plasma selenoprotein P concentrations have been found to be positively correlated with bone mineral density in elderly women.

General significance

A full understanding of the role and function of selenoproteins and Se status on bone physiology and pathology may lead to effectively prevent against or modify bone diseases by using Se.     

A tale of two toxicities: malformed selenoproteins and oxidative stress both contribute to selenium stress in plants

Background

Despite selenium''s toxicity in plants at higher levels, crops supply most of the essential dietary selenium in humans. In plants, inorganic selenium can be assimilated into selenocysteine, which can replace cysteine in proteins. Selenium toxicity in plants has been attributed to the formation of non-specific selenoproteins. However, this paradigm can be challenged now that there is increasingly abundant evidence suggesting that selenium-induced oxidative stress also contributes to toxicity in plants.

Scope

This Botanical Briefing summarizes the evidence indicating that selenium toxicity in plants is attributable to both the accumulation of non-specific selenoproteins and selenium-induced oxidative stress. Evidence is also presented to substantiate the claim that inadvertent selenocysteine replacement probably impairs or misfolds proteins, which supports the malformed selenoprotein hypothesis. The possible physiological ramifications of selenoproteins and selenium-induced oxidative stress are discussed.

Conclusions

Malformed selenoproteins and oxidative stress are two distinct types of stress that drive selenium toxicity in plants and could impact cellular processes in plants that have yet to be thoroughly explored. Although challenging, deciphering whether the extent of selenium toxicity in plants is imparted by selenoproteins or oxidative stress could be helpful in the development of crops with fortified levels of selenium.     

Selenoproteinless animals: selenophosphate synthetase SPS1 functions in a pathway unrelated to selenocysteine biosynthesis

Proteins containing the 21st amino acid, selenocysteine (Sec), have been described in all three domains of life, but the composition of selenoproteomes in organisms varies significantly. Here, we report that aquatic arthropods possess many selenoproteins also detected in other animals and unicellular eukaryotes, and that most of these proteins were either lost or replaced with cysteine-containing homologs in insects. As a result of this selective selenoproteome reduction, fruit flies and mosquitoes have three known selenoproteins, and the honeybee, Apis mellifera, a single detected candidate selenoprotein. Moreover, we identified the red flour beetle, Tribolium castaneum, and the silkworm, Bombyx mori, as the first animals that lack any Sec-containing proteins. These insects also lost the Sec biosynthesis and insertion machinery, but selenophosphate synthetase 1 (SPS1), an enzyme previously implicated in Sec biosynthesis, is present in all insects, including T. castaneum and B. mori. These data indicate that SPS1 functions in a pathway unrelated to selenoprotein synthesis. Since SPS1 evolved from a protein that utilizes selenium for Sec biosynthesis, an attractive possibility is that SPS1 may define a new pathway of selenium utilization in animals.     

Selenoproteins and selenocysteine insertion system in the model plant cell system,Chlamydomonas reinhardtii

Known eukaryotic selenocysteine (Sec)-containing proteins are animal proteins, whereas selenoproteins have not been found in yeast and plants. Surprisingly, we detected selenoproteins in a member of the plant kingdom, Chlamydomonas reinhardtii, and directly identified two of them as phospholipid hydroperoxide glutathione peroxidase and selenoprotein W homologs. Moreover, a selenocysteyl-tRNA was isolated that recognized specifically the Sec codon UGA. Subsequent gene cloning and bioinformatics analyses identified eight additional selenoproteins, including methionine-S-sulfoxide reductase, a selenoprotein specific to Chlamydomonas: Chlamydomonas selenoprotein genes contained selenocysteine insertion sequence (SECIS) elements that were similar, but not identical, to those of animals. These SECIS elements could direct selenoprotein synthesis in mammalian cells, indicating a common origin of plant and animal Sec insertion systems. We found that selenium is required for optimal growth of Chlamydomonas: Finally, evolutionary analyses suggested that selenoproteins present in Chlamydomonas and animals evolved early, and were independently lost in land plants, yeast and some animals.     

Increased rate of hair keratin gene loss in the cetacean lineage

Background

Hair represents an evolutionary innovation that appeared early on mammalian evolutionary history, and presumably contributed significantly to the rapid radiation of the group. An interesting event in hair evolution has been its secondary loss in some mammalian groups, such as cetaceans, whose hairless phenotype appears to be an adaptive response to better meet the environmental conditions. To determine whether different repertoire of keratin genes among mammals can potentially explain the phenotypic hair features of different lineages, we characterized the type I and II clusters of alpha keratins from eight mammalian species, including the hairless dolphin and minke whale representing the order Cetacea.

Results

We combined the available genomic information with phylogenetic analysis to conduct a comprehensive analysis of the evolutionary patterns of keratin gene clusters. We found that both type I and II gene clusters are fairly conserved among the terrestrial mammals included in this study, with lineage specific gene duplication and gene loss. Nevertheless, there is also evidence for an increased rate of pseudogenization in the cetacean lineage when compared to their terrestrial relatives, especially among the hair type keratins.

Conclusions

Here we present a comprehensive characterization of alpha-keratin genes among mammals and elucidate the mechanisms involved in the evolution of this gene family. We identified lineage-specific gene duplications and gene loss among the Laurasiatherian and Euarchontoglires species included in the study. Interestingly, cetaceans present an increased loss of hair-type keratin genes when compared to other terrestrial mammals. As suggested by the ‘less-is-more’ hypothesis, we do not rule out the possibility that the gene loss of hair-type keratin genes in these species might be associated to the hairless phenotype and could have been adaptive in response to new selective pressures imposed by the colonization of a new habitat. Our study provides support for the idea that pseudogenes are not simply ‘genomic fossils’ but instead have adaptive roles during the evolutionary process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-869) contains supplementary material, which is available to authorized users.     

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.     

MLGO: phylogeny reconstruction and ancestral inference from gene-order data

Background

The rapid accumulation of whole-genome data has renewed interest in the study of using gene-order data for phylogenetic analyses and ancestral reconstruction. Current software and web servers typically do not support duplication and loss events along with rearrangements.

Results

MLGO (Maximum Likelihood for Gene-Order Analysis) is a web tool for the reconstruction of phylogeny and/or ancestral genomes from gene-order data. MLGO is based on likelihood computation and shows advantages over existing methods in terms of accuracy, scalability and flexibility.

Conclusions

To the best of our knowledge, it is the first web tool for analysis of large-scale genomic changes including not only rearrangements but also gene insertions, deletions and duplications. The web tool is available from http://www.geneorder.org/server.php.     

SelT, SelW, SelH, and Rdx12: genomics and molecular insights into the functions of selenoproteins of a novel thioredoxin-like family

Selenium is an essential trace element in many life forms due to its occurrence as a selenocysteine (Sec) residue in selenoproteins. The majority of mammalian selenoproteins, however, have no known function. Herein, we performed extensive sequence similarity searches to define and characterize a new protein family, designated Rdx, that includes mammalian selenoproteins SelW, SelV, SelT and SelH, bacterial SelW-like proteins and cysteine-containing proteins of unknown function in all three domains of life. An additional member of this family is a mammalian cysteine-containing protein, designated Rdx12, and its fish selenoprotein orthologue. Rdx proteins are proposed to possess a thioredoxin-like fold and a conserved CxxC or CxxU (U is Sec) motif, suggesting a redox function. We cloned and characterized three mammalian members of this family, which showed distinct expression patterns in mouse tissues and different localization patterns in cells transfected with the corresponding GFP fusion proteins. By analogy to thioredoxin, Rdx proteins can use catalytic cysteine (or Sec) to form transient mixed disulfides with substrate proteins. We employed this property to identify cellular targets of Rdx proteins using affinity columns containing mutant versions of these proteins. Rdx12 was found to interact with glutathione peroxidase 1, whereas 14-3-3 protein was identified as one of the targets of mammalian SelW, suggesting a mechanism for redox regulation of the 14-3-3 family of proteins.     

Trichomonas Transmembrane Cyclases Result from Massive Gene Duplication and Concomitant Development of Pseudogenes

Background

Trichomonas vaginalis has an unusually large genome (∼160 Mb) encoding ∼60,000 proteins. With the goal of beginning to understand why some Trichomonas genes are present in so many copies, we characterized here a family of ∼123 Trichomonas genes that encode transmembrane adenylyl cyclases (TMACs).

Methodology/Principal Findings

The large family of TMACs genes is the result of recent duplications of a small set of ancestral genes that appear to be unique to trichomonads. Duplicated TMAC genes are not closely associated with repetitive elements, and duplications of flanking sequences are rare. However, there is evidence for TMAC gene replacements by homologous recombination. A high percentage of TMAC genes (∼46%) are pseudogenes, as they contain stop codons and/or frame shifts, or the genes are truncated. Numerous stop codons present in the genome project G3 strain are not present in orthologous genes of two other Trichomonas strains (S1 and B7RC2). Each TMAC is composed of a series of N-terminal transmembrane helices and a single C-terminal cyclase domain that has adenylyl cyclase activity. Multiple TMAC genes are transcribed by Trichomonas cloned by limiting dilution.

Conclusions/Significance

We conclude that one reason for the unusually large genome of Trichomonas is the presence of unstable families of genes such as those encoding TMACs that are undergoing massive gene duplication and concomitant development of pseudogenes.     

Diversity of Protein and mRNA Forms of Mammalian Methionine Sulfoxide Reductase B1 Due to Intronization and Protein Processing

Background

Methionine sulfoxide reductases (Msrs) are repair enzymes that protect proteins from oxidative stress by catalyzing stereospecific reduction of oxidized methionine residues. MsrB1 is a selenocysteine-containing cytosolic/nuclear Msr with high expression in liver and kidney.

Principal Findings

Here, we identified differences in MsrB1 gene structure among mammals. Human MsrB1 gene consists of four, whereas the corresponding mouse gene of five exons, due to occurrence of an additional intron that flanks the stop signal and covers a large part of the 3′-UTR. This intron evolved in a subset of rodents through intronization of exonic sequences, whereas the human gene structure represents the ancestral form. In mice, both splice forms were detected in liver, kidney, brain and heart with the five-exon form being the major form. We found that both mRNA forms were translated and supported efficient selenocysteine insertion into MsrB1. In addition, MsrB1 occurs in two protein forms that migrate as 14 and 5 kDa proteins. We found that each mRNA splice form generated both protein forms. The abundance of the 5 kDa form was not influenced by protease inhibitors, replacement of selenocysteine in the active site or mutation of amino acids in the cleavage site. However, mutation of cysteines that coordinate a structural zinc decreased the levels of 5 and 14 kDa forms, suggesting importance of protein structure for biosynthesis and/stability of these forms.

Conclusions

This study characterized unexpected diversity of protein and mRNA forms of mammalian selenoprotein MsrB1.     

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.     

A T/C polymorphism in the GPX4 3′UTR affects the selenoprotein expression pattern and cell viability in transfected Caco-2 cells

Background

Synthesis of selenoproteins such as glutathione peroxidases (GPx) requires a specific tRNA and a stem-loop structure in the 3′untranslated region (3′UTR) of the mRNA. A common single nucleotide polymorphism occurs in the GPX4 gene in a region corresponding to the 3′UTR.

Methods

The two variant 3′UTR sequences were linked to sequences from a selenoprotein reporter gene (iodothyronine deiodinase) and expressed in Caco-2 cells. Clones expressing comparable levels of deiodinase (assessed by real-time PCR) were selected and their response to tert-butyl hydroperoxide assessed by cell viability and measurement of reactive oxygen species. Selenoprotein expression was assessed by real-time PCR, enzyme activity and immunoassay.

Results

When selenium supply was low, cells overexpressing the C variant 3′UTR showed lower viability after oxidative challenge, increased levels of reactive oxygen species and lower GPx activity and SelH mRNA expression compared to cells overexpressing the T variant. After selenium supplementation, cell viability and GPx4 expression were higher in the cells overexpressing the C variant. Expression of transgenes incorporating the T/C variant GPX4 (rs713041) sequences in Caco-2 cells leads to alterations in both cell viability after an oxidative challenge and selenoprotein expression. This suggests that the two variants compete differently in the selenoprotein hierarchy.

General Significance

The data provide evidence that the T/C variant GPX4 (rs713041) alters the pattern of selenoprotein synthesis if selenium intake is low. Further work is required to assess the impact on disease susceptibility.