Characterization of naked mole-rat hematopoiesis reveals unique stem and progenitor cell patterns and neotenic traits.

Naked mole rats (NMRs) are the longest-lived rodents yet their stem cell characteristics remain enigmatic. Here, we comprehensively mapped the NMR hematopoietic landscape and identified unique features likely contributing to longevity. Adult NMRs form red blood cells in spleen and marrow, which comprise a myeloid bias toward granulopoiesis together with decreased B-lymphopoiesis. Remarkably, youthful blood and marrow single-cell transcriptomes and cell compositions are largely maintained until at least middle age. Similar to primates, the primitive stem and progenitor cell (HSPC) compartment is marked by CD34 and THY1. Stem cell polarity is seen for Tubulin but not CDC42, and is not lost until 12 years of age. HSPC respiration rates are as low as in purified human stem cells, in concert with a strong expression signature for fatty acid metabolism. The pool of quiescent stem cells is higher than in mice, and the cell cycle of hematopoietic cells is prolonged. By characterizing the NMR hematopoietic landscape, we identified resilience phenotypes such as an increased quiescent HSPC compartment, absence of age-related decline, and neotenic traits likely geared toward longevity.

The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

AMD1 mRNA employs ribosome stalling as a mechanism for molecular memory formation.

In addition to acting as template for protein synthesis, messenger RNA (mRNA) often contains sensory sequence elements that regulate this process. Here we report a new mechanism that limits the number of complete protein molecules that can be synthesized from a single mRNA molecule of the human AMD1 gene encoding adenosylmethionine decarboxylase 1 (AdoMetDC). A small proportion of ribosomes translating AMD1 mRNA stochastically read through the stop codon of the main coding region. These readthrough ribosomes then stall close to the next in-frame stop codon, eventually forming a ribosome queue, the length of which is proportional to the number of AdoMetDC molecules that were synthesized from the same AMD1 mRNA. Once the entire spacer region between the two stop codons is filled with queueing ribosomes, the queue impinges upon the main AMD1 coding region halting its translation. Phylogenetic analysis suggests that this mechanism is highly conserved in vertebrates and existed in their common ancestor. We propose that this mechanism is used to count and limit the number of protein molecules that can be synthesized from a single mRNA template. It could serve to safeguard from dysregulated translation that may occur owing to errors in transcription or mRNA damage.

Selenoprotein TXNRD3 supports male fertility via the redox regulation of spermatogenesis.

Thioredoxin/glutathione reductase (TXNRD3) is a selenoprotein composed of thioredoxin reductase and glutaredoxin domains. This NADPH-dependent thiol oxidoreductase evolved through gene duplication within the Txnrd family, is expressed in the testes, and can reduce both thioredoxin and glutathione in vitro; however, the function of this enzyme remains unknown. To characterize the function of TXNRD3 in vivo, we generated a strain of mice bearing deletion of Txnrd3 gene. We show that these Txnrd3 knockout mice are viable and without discernable gross phenotypes, and also that TXNRD3 deficiency leads to fertility impairment in male mice. We found that Txnrd3 knockout animals exhibited a lower fertilization rate in vitro, a sperm movement phenotype, and an altered thiol redox status in sperm cells. Proteomic analyses further revealed a broad range of substrates reduced by TXNRD3 during sperm maturation, presumably as a part of sperm quality control. Taken together, these results show that TXNRD3 plays a critical role in male reproduction via the thiol redox control of spermatogenesis.

Facultative protein selenation regulates redox sensitivity, adipose tissue thermogenesis, and obesity.

Oxidation of cysteine thiols by physiological reactive oxygen species (ROS) initiates thermogenesis in brown and beige adipose tissues. Cellular selenocysteines, where sulfur is replaced with selenium, exhibit enhanced reactivity with ROS. Despite their critical roles in physiology, methods for broad and direct detection of proteogenic selenocysteines are limited. Here we developed a mass spectrometric method to interrogate incorporation of selenium into proteins. Unexpectedly, this approach revealed facultative incorporation of selenium as selenocysteine or selenomethionine into proteins that lack canonical encoding for selenocysteine. Selenium was selectively incorporated into regulatory sites on key metabolic proteins, including as selenocysteine-replacing cysteine at position 253 in uncoupling protein 1 (UCP1). This facultative utilization of selenium was initiated by increasing cellular levels of organic, but not inorganic, forms of selenium. Remarkably, dietary selenium supplementation elevated facultative incorporation into UCP1, elevated energy expenditure through thermogenic adipose tissue, and protected against obesity. Together, these findings reveal the existence of facultative protein selenation, which correlates with impacts on thermogenic adipocyte function and presumably other biological processes as well.

Tolerance to Selenoprotein Loss Differs between Human and Mouse.

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

Multiple RNA structures affect translation initiation and UGA redefinition efficiency during synthesis of selenoprotein P.

Gene-specific expansion of the genetic code allows for UGA codons to specify the amino acid selenocysteine (Sec). A striking example of UGA redefinition occurs during translation of the mRNA coding for the selenium transport protein, selenoprotein P (SELENOP), which in vertebrates may contain up to 22 in-frame UGA codons. Sec incorporation at the first and downstream UGA codons occurs with variable efficiencies to control synthesis of full-length and truncated SELENOP isoforms. To address how the Selenop mRNA can direct dynamic codon redefinition in different regions of the same mRNA, we undertook a comprehensive search for phylogenetically conserved RNA structures and examined the function of these structures using cell-based assays, in vitro translation systems, and in vivo ribosome profiling of liver tissue from mice carrying genomic deletions of 3' UTR selenocysteine-insertion-sequences (SECIS1 and SECIS2). The data support a novel RNA structure near the start codon that impacts translation initiation, structures located adjacent to UGA codons, additional coding sequence regions necessary for efficient production of full-length SELENOP, and distinct roles for SECIS1 and SECIS2 at UGA codons. Our results uncover a remarkable diversity of RNA elements conducting multiple occurrences of UGA redefinition to control the synthesis of full-length and truncated SELENOP isoforms.

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.

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 (tRNASec) drives the recoding of highly specific UGA codons from stop signals to Sec. Although found in organisms from the three domains of life, Sec is not universal. Many species are completely devoid of selenoprotein genes and lack the ability to synthesize Sec. Since tRNASec is a key component in selenoprotein biosynthesis, its efficient identification in genomes is instrumental to characterize the utilization of Sec across lineages. Available tRNA prediction methods fail to accurately predict tRNASec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNASec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNASec predictions, which is not possible for the predictions of other tRNAs. We show that Secmarker greatly improves the accuracy of previously existing methods constituting a valuable tool to identify tRNASec genes, and to efficiently determine whether a genome contains selenoproteins. We used Secmarker to analyze a large set of fully sequenced genomes, and the results revealed new insights in the biology of tRNASec, led to the discovery of a novel bacterial selenoprotein family, and shed additional light on the phylogenetic distribution of selenoprotein containing genomes. Secmarker is freely accessible for download, or online analysis through a web server at http://secmarker.crg.cat.

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.

Novel Ciliate Genetic Code Variants Including the Reassignment of All Three Stop Codons to Sense Codons in Condylostoma magnum.

mRNA translation in many ciliates utilizes variant genetic codes where stop codons are reassigned to specify amino acids. To characterize the repertoire of ciliate genetic codes, we analyzed ciliate transcriptomes from marine environments. Using codon substitution frequencies in ciliate protein-coding genes and their orthologs, we inferred the genetic codes of 24 ciliate species. Nine did not match genetic code tables currently assigned by NCBI. Surprisingly, we identified a novel genetic code where all three standard stop codons (TAA, TAG, and TGA) specify amino acids in Condylostoma magnum We provide evidence suggesting that the functions of these codons in C. magnum depend on their location within mRNA. They are decoded as amino acids at internal positions, but specify translation termination when in close proximity to an mRNA 3' end. The frequency of stop codons in protein coding sequences of closely related Climacostomum virens suggests that it may represent a transitory state.

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.

SelenoDB 2.0: annotation of selenoprotein genes in animals and their genetic diversity in humans.

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

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.

Pseudogenization of testis-specific Lfg5 predates human/Neanderthal divergence.

Recent reviews discussed the critical roles of apoptosis in human spermatogenesis and infertility. These reviews highlight the FasL-induced caspase cascade in apoptosis lending importance to our discovery of the pseudogene status of the Lfg5 gene in modern humans, Neanderthal and the Denisovan. This gene is a member of the ancient and highly conserved apoptosis Lifeguard family. This pseudogenization is the result of a premature stop codon at the 3'-end of exon 8 not found in any other ortholog. With the current exception of the domesticated bovine and buffalo, Lfg5's expression in mammals is testis-specific. A full analysis of this gene, its phylogenetic context and its recent hominin changes suggest its inactivation was likely under selection in human evolution.

The response of durum wheat to the preceding crop in a Mediterranean environment.

Crop sequence is an important management practice that may affect durum wheat (Triticum durum Desf.) production. Field research was conducted in 2007-2008 and 2008-2009 seasons in a rain-fed cold Mediterranean environment to examine the impact of the preceding crops alfalfa (Medicago sativa L.), maize (Zea mays L.), sunflower (Helianthus annuus L.), and bread wheat (Triticum aestivum L.) on yield and N uptake of four durum wheat varieties. The response of grain yield of durum wheat to the preceding crop was high in 2007-2008 and was absent in the 2008-2009 season, because of the heavy rainfall that negatively impacted establishment, vegetative growth, and grain yield of durum wheat due to waterlogging. In the first season, durum wheat grain yield was highest following alfalfa, and was 33% lower following wheat. The yield increase of durum wheat following alfalfa was mainly due to an increased number of spikes per unit area and number of kernels per spike, while the yield decrease following wheat was mainly due to a reduction of spike number per unit area. Variety growth habit and performance did not affect the response to preceding crop and varieties ranked in the order Levante > Saragolla = Svevo > Normanno.

The selenophosphate synthetase family: A review.

Selenophosphate synthetases use selenium and ATP to synthesize selenophosphate. This is required for biological utilization of selenium, most notably for the synthesis of the non-canonical amino acid selenocysteine (Sec). Therefore, selenophosphate synthetases underlie all functions of selenoproteins, which include redox homeostasis, protein quality control, hormone regulation, metabolism, and many others. This protein family comprises two groups, SelD/SPS2 and SPS1. The SelD/SPS2 group represent true selenophosphate synthetases, enzymes central to selenium metabolism which are present in all Sec-utilizing organisms across the tree of life. Notably, many SelD/SPS2 proteins contain Sec as catalytic residue in their N-terminal flexible selenium-binding loop, while others replace it with cysteine (Cys). The SPS1 group comprises proteins originated through gene duplications of SelD/SPS2 in metazoa in which the Sec/Cys-dependent catalysis was disrupted. SPS1 proteins do not synthesize selenophosphate and are not required for Sec synthesis. They have essential regulatory functions related to redox homeostasis and pyridoxal phosphate, which affect signaling pathways for growth and differentiation. In this review, we summarize the knowledge about the selenophosphate synthetase family acquired through decades of research, encompassing their structure, mechanism, function, and evolution.

Genetic diversity of sheep breeds from Albania, Greece, and Italy assessed by mitochondrial DNA and nuclear polymorphisms (SNPs).

We employed mtDNA and nuclear SNPs to investigate the genetic diversity of sheep breeds of three countries of the Mediterranean basin: Albania, Greece, and Italy. In total, 154 unique mtDNA haplotypes were detected by means of D-loop sequence analysis. The major nucleotide diversity was observed in Albania. We identified haplogroups, A, B, and C in Albanian and Greek samples, while Italian individuals clustered in groups A and B. In general, the data show a pattern reflecting old migrations that occurred in postneolithic and historical times. PCA analysis on SNP data differentiated breeds with good correspondence to geographical locations. This could reflect geographical isolation, selection operated by local sheep farmers, and different flock management and breed admixture that occurred in the last centuries.