Selenoprotein Gene Nomenclature.

The human genome contains 25 genes coding for selenocysteine-containing proteins (selenoproteins). These proteins are involved in a variety of functions, most notably redox homeostasis. Selenoprotein enzymes with known functions are designated according to these functions: TXNRD1, TXNRD2, and TXNRD3 (thioredoxin reductases), GPX1, GPX2, GPX3, GPX4, and GPX6 (glutathione peroxidases), DIO1, DIO2, and DIO3 (iodothyronine deiodinases), MSRB1 (methionine sulfoxide reductase B1), and SEPHS2 (selenophosphate synthetase 2). Selenoproteins without known functions have traditionally been denoted by SEL or SEP symbols. However, these symbols are sometimes ambiguous and conflict with the approved nomenclature for several other genes. Therefore, there is a need to implement a rational and coherent nomenclature system for selenoprotein-encoding genes. Our solution is to use the root symbol SELENO followed by a letter. This nomenclature applies to SELENOF (selenoprotein F, the 15-kDa selenoprotein, SEP15), SELENOH (selenoprotein H, SELH, C11orf31), SELENOI (selenoprotein I, SELI, EPT1), SELENOK (selenoprotein K, SELK), SELENOM (selenoprotein M, SELM), SELENON (selenoprotein N, SEPN1, SELN), SELENOO (selenoprotein O, SELO), SELENOP (selenoprotein P, SeP, SEPP1, SELP), SELENOS (selenoprotein S, SELS, SEPS1, VIMP), SELENOT (selenoprotein T, SELT), SELENOV (selenoprotein V, SELV), and SELENOW (selenoprotein W, SELW, SEPW1). This system, approved by the HUGO Gene Nomenclature Committee, also resolves conflicting, missing, and ambiguous designations for selenoprotein genes and is applicable to selenoproteins across vertebrates.

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.

Genome sequencing reveals insights into physiology and longevity of the naked mole rat.

The naked mole rat (Heterocephalus glaber) is a strictly subterranean, extraordinarily long-lived eusocial mammal. Although it is the size of a mouse, its maximum lifespan exceeds 30 years, making this animal the longest-living rodent. Naked mole rats show negligible senescence, no age-related increase in mortality, and high fecundity until death. In addition to delayed ageing, they are resistant to both spontaneous cancer and experimentally induced tumorigenesis. Naked mole rats pose a challenge to the theories that link ageing, cancer and redox homeostasis. Although characterized by significant oxidative stress, the naked mole rat proteome does not show age-related susceptibility to oxidative damage or increased ubiquitination. Naked mole rats naturally reside in large colonies with a single breeding female, the 'queen', who suppresses the sexual maturity of her subordinates. They also live in full darkness, at low oxygen and high carbon dioxide concentrations, and are unable to sustain thermogenesis nor feel certain types of pain. Here we report the sequencing and analysis of the naked mole rat genome, which reveals unique genome features and molecular adaptations consistent with cancer resistance, poikilothermy, hairlessness and insensitivity to low oxygen, and altered visual function, circadian rythms and taste sensing. This information provides insights into the naked mole rat's exceptional longevity and ability to live in hostile conditions, in the dark and at low oxygen. The extreme traits of the naked mole rat, together with the reported genome and transcriptome information, offer opportunities for understanding ageing and advancing other areas of biological and biomedical research.

Selenoprotein S is involved in maintenance and transport of multiprotein complexes.

SelS (Selenoprotein S) is a selenocysteine-containing protein with roles in ER (endoplasmic reticulum) function and inflammation. It has been implicated in ERAD (ER-associated protein degradation), and clinical studies revealed an association of its promoter polymorphism with cytokine levels and human diseases. However, the pathways and interacting proteins that could shed light on pathogenesis of SelS-associated diseases have not been studied systematically. We performed a large-scale affinity isolation of human SelS and its mutant forms and analysed the proteins that interact with them. All previously known SelS targets and nearly two hundred additional proteins were identified that were remarkably enriched for various multiprotein complexes. Subsequent chemical cross-linking experiments identified the specific interacting sites in SelS and its several targets. Most of these interactions involved coiled-coil domains. The data suggest that SelS participates in intracellular membrane transport and maintenance of protein complexes by anchoring them to the ER membrane.

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.

UGA codon position-dependent incorporation of selenocysteine into mammalian selenoproteins.

It is thought that the SelenoCysteine Insertion Sequence (SECIS) element and UGA codon are sufficient for selenocysteine (Sec) insertion. However, we found that UGA supported Sec insertion only at its natural position or in its close proximity in mammalian thioredoxin reductase 1 (TR1). In contrast, Sec could be inserted at any tested position in mammalian TR3. Replacement of the 3'-UTR of TR3 with the corresponding segment of a Euplotes crassus TR restricted Sec insertion into the C-terminal region, whereas the 3'-UTR of TR3 conferred unrestricted Sec insertion into E. crassus TR, in which Sec insertion is normally limited to the C-terminal region. Exchanges of 3'-UTRs between mammalian TR1 and E. crassus TR had no effect, as both proteins restricted Sec insertion. We further found that these effects could be explained by the use of selenoprotein-specific SECIS elements. Examination of Sec insertion into other selenoproteins was consistent with this model. The data indicate that mammals evolved the ability to limit Sec insertion into natural positions within selenoproteins, but do so in a selenoprotein-specific manner, and that this process is controlled by the SECIS element in the 3'-UTR.

Genome-wide ribosome profiling reveals complex translational regulation in response to oxidative stress.

Information on unique and coordinated regulation of transcription and translation in response to stress is central to the understanding of cellular homeostasis. Here we used ribosome profiling coupled with next-generation sequencing to examine the interplay between transcription and translation under conditions of hydrogen peroxide treatment in Saccharomyces cerevisiae. Hydrogen peroxide treatment led to a massive and rapid increase in ribosome occupancy of short upstream ORFs, including those with non-AUG translational starts, and of the N-terminal regions of ORFs that preceded the transcriptional response. In addition, this treatment induced the synthesis of N-terminally extended proteins and elevated stop codon read-through and frameshift events. It also increased ribosome occupancy at the beginning of ORFs and potentially the duration of the elongation step. We identified proteins whose synthesis was regulated rapidly by hydrogen peroxide posttranscriptionally; however, for the majority of genes increased protein synthesis followed transcriptional regulation. These data define the landscape of genome-wide regulation of translation in response to hydrogen peroxide and suggest that potentiation (coregulation of the transcript level and translation) is a feature of oxidative stress.

Niche of harmful alga Aureococcus anophagefferens revealed through ecogenomics.

Harmful algal blooms (HABs) cause significant economic and ecological damage worldwide. Despite considerable efforts, a comprehensive understanding of the factors that promote these blooms has been lacking, because the biochemical pathways that facilitate their dominance relative to other phytoplankton within specific environments have not been identified. Here, biogeochemical measurements showed that the harmful alga Aureococcus anophagefferens outcompeted co-occurring phytoplankton in estuaries with elevated levels of dissolved organic matter and turbidity and low levels of dissolved inorganic nitrogen. We subsequently sequenced the genome of A. anophagefferens and compared its gene complement with those of six competing phytoplankton species identified through metaproteomics. Using an ecogenomic approach, we specifically focused on gene sets that may facilitate dominance within the environmental conditions present during blooms. A. anophagefferens possesses a larger genome (56 Mbp) and has more genes involved in light harvesting, organic carbon and nitrogen use, and encoding selenium- and metal-requiring enzymes than competing phytoplankton. Genes for the synthesis of microbial deterrents likely permit the proliferation of this species, with reduced mortality losses during blooms. Collectively, these findings suggest that anthropogenic activities resulting in elevated levels of turbidity, organic matter, and metals have opened a niche within coastal ecosystems that ideally suits the unique genetic capacity of A. anophagefferens and thus, has facilitated the proliferation of this and potentially other HABs.

Reduced utilization of selenium by naked mole rats due to a specific defect in GPx1 expression.

Naked mole rat (MR) Heterocephalus glaber is a rodent model of delayed aging because of its unusually long life span (>28 years). It is also not known to develop cancer. In the current work, tissue imaging by x-ray fluorescence microscopy and direct analyses of trace elements revealed low levels of selenium in the MR liver and kidney, whereas MR and mouse brains had similar selenium levels. This effect was not explained by uniform selenium deficiency because methionine sulfoxide reductase activities were similar in mice and MR. However, glutathione peroxidase activity was an order of magnitude lower in MR liver and kidney than in mouse tissues. In addition, metabolic labeling of MR cells with (75)Se revealed a loss of the abundant glutathione peroxidase 1 (GPx1) band, whereas other selenoproteins were preserved. To characterize the MR selenoproteome, we sequenced its liver transcriptome. Gene reconstruction revealed standard selenoprotein sequences except for GPx1, which had an early stop codon, and SelP, which had low selenocysteine content. When expressed in HEK 293 cells, MR GPx1 was present in low levels, and its expression could be rescued neither by removing the early stop codon nor by replacing its SECIS element. In addition, GPx1 mRNA was present in lower levels in MR liver than in mouse liver. To determine if GPx1 deficiency could account for the reduced selenium content, we analyzed GPx1 knock-out mice and found reduced selenium levels in their livers and kidneys. Thus, MR is characterized by the reduced utilization of selenium due to a specific defect in GPx1 expression.

Using DNA Methylation Profiling to Evaluate Biological Age and Longevity Interventions.

The DNA methylation levels of certain CpG sites are thought to reflect the pace of human aging. Here, we developed a robust predictor of mouse biological age based on 90 CpG sites derived from partial blood DNA methylation profiles. The resulting clock correctly determines the age of mouse cohorts, detects the longevity effects of calorie restriction and gene knockouts, and reports rejuvenation of fibroblast-derived iPSCs. The data show that mammalian DNA methylomes are characterized by CpG sites that may represent the organism's biological age. They are scattered across the genome, they are distinct in human and mouse, and their methylation gradually changes with age. The clock derived from these sites represents a biomarker of aging and can be used to determine the biological age of organisms and evaluate interventions that alter the rate of aging.

Naked Mole Rat Induced Pluripotent Stem Cells and Their Contribution to Interspecific Chimera.

Naked mole rats (NMRs) are exceptionally long-lived, cancer-resistant rodents. Identifying the defining characteristics of these traits may shed light on aging and cancer mechanisms. Here, we report the generation of induced pluripotent stem cells (iPSCs) from NMR fibroblasts and their contribution to mouse-NMR chimeric embryos. Efficient reprogramming could be observed under N2B27+2i conditions. The iPSCs displayed a characteristic morphology, expressed pluripotent markers, formed embryoid bodies, and showed typical differentiation patterns. Interestingly, NMR embryonic fibroblasts and the derived iPSCs had propensity for a tetraploid karyotype and were resistant to forming teratomas, but within mouse blastocysts they contributed to both interspecific placenta and fetus. Gene expression patterns of NMR iPSCs were more similar to those of human than mouse iPSCs. Overall, we uncovered unique features of NMR iPSCs and report a mouse-NMR chimeric model. The iPSCs and associated cell culture systems can be used for a variety of biological and biomedical applications.

Position-dependent termination and widespread obligatory frameshifting in Euplotes translation.

The ribosome can change its reading frame during translation in a process known as programmed ribosomal frameshifting. These rare events are supported by complex mRNA signals. However, we found that the ciliates Euplotes crassus and Euplotes focardii exhibit widespread frameshifting at stop codons. 47 different codons preceding stop signals resulted in either +1 or +2 frameshifts, and +1 frameshifting at AAA was the most frequent. The frameshifts showed unusual plasticity and rapid evolution, and had little influence on translation rates. The proximity of a stop codon to the 3' mRNA end, rather than its occurrence or sequence context, appeared to designate termination. Thus, a 'stop codon' is not a sufficient signal for translation termination, and the default function of stop codons in Euplotes is frameshifting, whereas termination is specific to certain mRNA positions and probably requires additional factors.

Analysis of cancer genomes reveals basic features of human aging and its role in cancer development.

Somatic mutations have long been implicated in aging and disease, but their impact on fitness and function is difficult to assess. Here by analysing human cancer genomes we identify mutational patterns associated with aging. Our analyses suggest that age-associated mutation load and burden double approximately every 8 years, similar to the all-cause mortality doubling time. This analysis further reveals variance in the rate of aging among different human tissues, for example, slightly accelerated aging of the reproductive system. Age-adjusted mutation load and burden correlate with the corresponding cancer incidence and precede it on average by 15 years, pointing to pre-clinical cancer development times. Behaviour of mutation load also exhibits gender differences and late-life reversals, explaining some gender-specific and late-life patterns in cancer incidence rates. Overall, this study characterizes some features of human aging and offers a mechanism for age being a risk factor for the onset of cancer.

Evidence that mutation accumulation does not cause aging in Saccharomyces cerevisiae.

The concept that mutations cause aging phenotypes could not be directly tested previously due to inability to identify age-related mutations in somatic cells and determine their impact on organismal aging. Here, we subjected Saccharomyces cerevisiae to multiple rounds of replicative aging and assessed de novo mutations in daughters of mothers of different age. Mutations did increase with age, but their low numbers, < 1 per lifespan, excluded their causal role in aging. Structural genome changes also had no role. A mutant lacking thiol peroxidases had the mutation rate well above that of wild-type cells, but this did not correspond to the aging pattern, as old wild-type cells with few or no mutations were dying, whereas young mutant cells with many more mutations continued dividing. In addition, wild-type cells lost mitochondrial DNA during aging, whereas shorter-lived mutant cells preserved it, excluding a causal role of mitochondrial mutations in aging. Thus, DNA mutations do not cause aging in yeast. These findings may apply to other damage types, suggesting a causal role of cumulative damage, as opposed to individual damage types, in organismal aging.

Gene expression defines natural changes in mammalian lifespan.

Mammals differ more than 100-fold in maximum lifespan, which can be altered in either direction during evolution, but the molecular basis for natural changes in longevity is not understood. Divergent evolution of mammals also led to extensive changes in gene expression within and between lineages. To understand the relationship between lifespan and variation in gene expression, we carried out RNA-seq-based gene expression analyses of liver, kidney, and brain of 33 diverse species of mammals. Our analysis uncovered parallel evolution of gene expression and lifespan, as well as the associated life-history traits, and identified the processes and pathways involved. These findings provide direct insights into how nature reversibly adjusts lifespan and other traits during adaptive radiation of lineages.

Thiol peroxidase deficiency leads to increased mutational load and decreased fitness in Saccharomyces cerevisiae.

Thiol peroxidases are critical enzymes in the redox control of cellular processes that function by reducing low levels of hydroperoxides and regulating redox signaling. These proteins were also shown to regulate genome stability, but how their dysfunction affects the actual mutations in the genome is not known. Saccharomyces cerevisiae has eight thiol peroxidases of glutathione peroxidase and peroxiredoxin families, and the mutant lacking all these genes (∆8) is viable. In this study, we employed two independent ∆8 isolates to analyze the genome-wide mutation spectrum that results from deficiency in these enzymes. Deletion of these genes was accompanied by a dramatic increase in point mutations, many of which clustered in close proximity and scattered throughout the genome, suggesting strong mutational bias. We further subjected multiple lines of wild-type and ∆8 cells to long-term mutation accumulation, followed by genome sequencing and phenotypic characterization. ∆8 lines showed a significant increase in nonrecurrent point mutations and indels. The original ∆8 cells exhibited reduced growth rate and decreased life span, which were further reduced in all ∆8 mutation accumulation lines. Although the mutation spectrum of the two independent isolates was different, similar patterns of gene expression were observed, suggesting the direct contribution of thiol peroxidases to the observed phenotypes. Expression of a single thiol peroxidase could partially restore the growth phenotype of ∆8 cells. This study shows how deficiency in nonessential, yet critical and conserved oxidoreductase function, leads to increased mutational load and decreased fitness.

Intrinsic expression of a multiexon type 3 deiodinase gene controls zebrafish embryo size.

Thyroid hormone is a master regulator of differentiation and growth, and its action is terminated by the enzymatic removal of an inner-ring iodine catalyzed by the selenoenzyme type 3 deiodinase (dio3). Our studies of the zebrafish reveal that the dio3 gene is duplicated in this species and that embryonic deiodination is an important determinant of embryo size. Although both dio3 paralogs encode enzymatically active proteins with high affinity for thyroid hormones, their anatomic patterns of expression are markedly divergent and only embryos with knockdown of dio3b, a biallelically expressed selenoenzyme expressed in the developing central nervous system, manifest severe thyroid hormone-dependent growth restriction at 72 hours post fertilization. This indicates that the embryonic deficiency of dio3, once considered only a placental enzyme, causes microsomia independently of placental physiology and raises the intriguing possibility that fetal abnormalities in human deiodination may present as intrauterine growth retardation. By mapping the gene structures and enzymatic properties of all four zebrafish deiodinases, we also identify dio3b as the first multiexon dio3 gene, containing a large intron separating its open reading frame from its selenocysteine insertion sequence (SECIS) element.

Adaptations to a subterranean environment and longevity revealed by the analysis of mole rat genomes.

Subterranean mammals spend their lives in dark, unventilated environments that are rich in carbon dioxide and ammonia and low in oxygen. Many of these animals are also long-lived and exhibit reduced aging-associated diseases, such as neurodegenerative disorders and cancer. We sequenced the genome of the Damaraland mole rat (DMR, Fukomys damarensis) and improved the genome assembly of the naked mole rat (NMR, Heterocephalus glaber). Comparative genome analyses, along with the transcriptomes of related subterranean rodents, revealed candidate molecular adaptations for subterranean life and longevity, including a divergent insulin peptide, expression of oxygen-carrying globins in the brain, prevention of high CO2-induced pain perception, and enhanced ammonia detoxification. Juxtaposition of the genomes of DMR and other more conventional animals with the genome of NMR revealed several truly exceptional NMR features: unusual thermogenesis, an aberrant melatonin system, pain insensitivity, and unique processing of 28S rRNA. Together, these genomes and transcriptomes extend our understanding of subterranean adaptations, stress resistance, and longevity.

The central role of selenium in the biochemistry and ecology of the harmful pelagophyte, Aureococcus anophagefferens.

The trace element selenium (Se) is required for the biosynthesis of selenocysteine (Sec), the 21st amino acid in the genetic code, but its role in the ecology of harmful algal blooms (HABs) is unknown. Here, we examined the role of Se in the biology and ecology of the harmful pelagophyte, Aureococcus anophagefferens, through cell culture, genomic analyses, and ecosystem studies. This organism has the largest and the most diverse selenoproteome identified to date that consists of at least 59 selenoproteins, including known eukaryotic selenoproteins, selenoproteins previously only detected in bacteria, and novel selenoproteins. The A. anophagefferens selenoproteome was dominated by the thioredoxin fold proteins and oxidoreductase functions were assigned to the majority of detected selenoproteins. Insertion of Sec in these proteins was supported by a unique Sec insertion sequence. Se was required for the growth of A. anophagefferens as cultures grew maximally at nanomolar Se concentrations. In a coastal ecosystem, dissolved Se concentrations were elevated before and after A. anophagefferens blooms, but were reduced by >95% during the peak of blooms to 0.05 nM. Consistent with this pattern, enrichment of seawater with selenite before and after a bloom did not affect the growth of A. anophagefferens, but enrichment during the peak of the bloom significantly increased population growth rates. These findings demonstrate that Se inventories, which can be anthropogenically enriched, can support proliferation of HABs, such as A. anophagefferens through its synthesis of a large arsenal of Se-dependent oxidoreductases that fine-tune cellular redox homeostasis.

Composition and evolution of the vertebrate and mammalian selenoproteomes.

BACKGROUND: Selenium is an essential trace element in mammals due to its presence in proteins in the form of selenocysteine (Sec). Human genome codes for 25 Sec-containing protein genes, and mouse and rat genomes for 24. METHODOLOGY/PRINCIPAL FINDINGS: We characterized the selenoproteomes of 44 sequenced vertebrates by applying gene prediction and phylogenetic reconstruction methods, supplemented with the analyses of gene structures, alternative splicing isoforms, untranslated regions, SECIS elements, and pseudogenes. In total, we detected 45 selenoprotein subfamilies. 28 of them were found in mammals, and 41 in bony fishes. We define the ancestral vertebrate (28 proteins) and mammalian (25 proteins) selenoproteomes, and describe how they evolved along lineages through gene duplication (20 events), gene loss (10 events) and replacement of Sec with cysteine (12 events). We show that an intronless selenophosphate synthetase 2 gene evolved in early mammals and replaced functionally the original multiexon gene in placental mammals, whereas both genes remain in marsupials. Mammalian thioredoxin reductase 1 and thioredoxin-glutathione reductase evolved from an ancestral glutaredoxin-domain containing enzyme, still present in fish. Selenoprotein V and GPx6 evolved specifically in placental mammals from duplications of SelW and GPx3, respectively, and GPx6 lost Sec several times independently. Bony fishes were characterized by duplications of several selenoprotein families (GPx1, GPx3, GPx4, Dio3, MsrB1, SelJ, SelO, SelT, SelU1, and SelW2). Finally, we report identification of new isoforms for several selenoproteins and describe unusually conserved selenoprotein pseudogenes. CONCLUSIONS/SIGNIFICANCE: This analysis represents the first comprehensive survey of the vertebrate and mammal selenoproteomes, and depicts their evolution along lineages. It also provides a wealth of information on these selenoproteins and their forms.