Downregulation of mitochondrial metabolism is a driver for fast skeletal muscle loss during mouse aging.

Skeletal muscle aging is characterized by the loss of muscle mass, strength and function, mainly attributed to the atrophy of glycolytic fibers. Underlying mechanisms driving the skeletal muscle functional impairment are yet to be elucidated. To unbiasedly uncover its molecular mechanisms, we recurred to gene expression and metabolite profiling in a glycolytic muscle, Extensor digitorum longus (EDL), from young and aged C57BL/6JRj mice. Employing multi-omics approaches we found that the main age-related changes are connected to mitochondria, exhibiting a downregulation in mitochondrial processes. Consistent is the altered mitochondrial morphology. We further compared our mouse EDL aging signature with human data from the GTEx database, reinforcing the idea that our model may recapitulate muscle loss in humans. We are able to show that age-related mitochondrial downregulation is likely to be detrimental, as gene expression signatures from commonly used lifespan extending interventions displayed the opposite direction compared to our EDL aging signature.

Selenocysteine Machinery Primarily Supports TXNRD1 and GPX4 Functions and Together They Are Functionally Linked with SCD and PRDX6.

The human genome has 25 genes coding for selenocysteine (Sec)-containing proteins, whose synthesis is supported by specialized Sec machinery proteins. Here, we carried out an analysis of the co-essentiality network to identify functional partners of selenoproteins and Sec machinery. One outstanding cluster included all seven known Sec machinery proteins and two critical selenoproteins, GPX4 and TXNRD1. Additionally, these nine genes were further positively associated with PRDX6 and negatively with SCD, linking the latter two genes to the essential role of selenium. We analyzed the essentiality scores of gene knockouts in this cluster across one thousand cancer cell lines and found that Sec metabolism genes are strongly selective for a subset of primary tissues, suggesting that certain cancer cell lineages are particularly dependent on selenium. A separate outstanding cluster included selenophosphate synthetase SEPHS1, which was linked to a group of transcription factors, whereas the remaining selenoproteins were linked neither to these clusters nor among themselves. The data suggest that key components of Sec machinery have already been identified and that their primary role is to support the functions of GPX4 and TXNRD1, with further functional links to PRDX6 and SCD.

Pathogenic Variants in Selenoproteins and Selenocysteine Biosynthesis Machinery.

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

COVID-19 mortality rate in children is U-shaped.

Children are known to be better protected from COVID-19 than adults, but their susceptibility patterns and the risk relative to other diseases are insufficiently defined. Here, we found that the COVID-19 mortality rate is U-shaped in childhood: it initially decreases, reaching the minimum at the ages 3-10 years, and then increases throughout life. All-cause mortality and mortality from other diseases, such as pneumonia and influenza, show a similar pattern; however, childhood mortality rates from COVID-19 are considerably lower than from other diseases, with the best relative protection achieved at the youngest ages. Consistent with this, the fraction of COVID-19 deaths among all deaths increases as a function of age throughout childhood and the entire life. We discuss implications of the elevated postnatal COVID-19 risk and lower childhood COVID-19 mortality compared to other diseases.

Historical Roles of Selenium and Selenoproteins in Health and Development: The Good, the Bad and the Ugly.

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

COVID-19 is an emergent disease of aging.

COVID-19 is an ongoing pandemic caused by the SARS-CoV-2 coronavirus that poses one of the greatest challenges to public health in recent years. SARS-CoV-2 is known to preferentially target older subjects and those with pre-existing conditions, but the reason for this age dependence is unclear. Here, we found that the case fatality rate for COVID-19 grows exponentially with age in all countries tested, with the doubling time approaching that of all-cause human mortality. In addition, men and those with multiple age-related diseases are characterized by increased mortality. Moreover, similar mortality patterns were found for all-cause pneumonia. We further report that the gene expression of ACE2, the SARS-CoV-2 receptor, grows in the lung with age, except for subjects on a ventilator. Together, our findings establish COVID-19 as an emergent disease of aging, and age and age-related diseases as its major risk factors. In turn, this suggests that COVID-19, and deadly respiratory diseases in general, may be targeted, in addition to antiviral approaches, by approaches that target the aging process.

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.

Germline burden of rare damaging variants negatively affects human healthspan and lifespan.

Heritability of human lifespan is 23-33% as evident from twin studies. Genome-wide association studies explored this question by linking particular alleles to lifespan traits. However, genetic variants identified so far can explain only a small fraction of lifespan heritability in humans. Here, we report that the burden of rarest protein-truncating variants (PTVs) in two large cohorts is negatively associated with human healthspan and lifespan, accounting for 0.4 and 1.3 years of their variability, respectively. In addition, longer-living individuals possess both fewer rarest PTVs and less damaging PTVs. We further estimated that somatic accumulation of PTVs accounts for only a small fraction of mortality and morbidity acceleration and hence is unlikely to be causal in aging. We conclude that rare damaging mutations, both inherited and accumulated throughout life, contribute to the aging process, and that burden of ultra-rare variants in combination with common alleles better explain apparent heritability of human lifespan.

Most living things undergo biological changes as they get older, a process that we generally refer to as aging. Despite being a widespread phenomenon, scientists do not fully understand why we age, though it appears that a combination of genetics and lifestyle factors, such as diet, play a role in influencing lifespan. Aging increases the risk of developing a wide range of diseases, including cancer, Alzheimer's disease and diabetes. As such, finding ways to slow the aging process would help to postpone the onset of illness and potentially improve health in old age. Genes are thought to be responsible for between one quarter and one third of the variation in human lifespans. The relationship between genes, aging and lifespan is complex and not well understood. One set of rare genetic changes that have been shown to have significant effects on diseases are called protein truncation variants (PTVs). PTVs cause damage by altering the production of certain proteins. There are many possible PTVs and people can be born with them or they can develop them in some cells later in life. The full influence of PTVs on aging is not known. Shindyapina, Zenin et al. have now studied observational data collected from two groups of over 40,000 people in the UK. Both groups recorded over 1,000 deaths, and the study examined the influence of PTVs on natural lifespan. The results show that each person is born with an average of six PTVs, which can vary in the impact that they have on aging. Having more, or more severe, PTVs could reduce life expectancy on average by 1.3 years. PTVs affect both total lifespan and healthy lifespan, the period of time lived prior to developing the first age-related disease. While PTVs that people are born with have a significant effect on aging, this study also showed that PTVs that are acquired due to spontaneous mutations through a person's life have much less of an impact. This is a key insight into the relationship between genes and aging. These discoveries could help in using genetics to anticipate future health, it also helps to identify some of the biological systems that have a role in aging. This could lead to new ways to delay the aging process and its effects on health.

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.

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.

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.

Selenoprofiles: A Computational Pipeline for Annotation of Selenoproteins.

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

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.

Human selenoprotein P and S variant mRNAs with different numbers of SECIS elements and inferences from mutant mice of the roles of multiple SECIS elements.

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

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.

Identification of a selenium-dependent glutathione peroxidase in the blood-sucking insect Rhodnius prolixus.

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

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.

The genomes of two key bumblebee species with primitive eusocial organization.

BACKGROUND: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. RESULTS: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. CONCLUSIONS: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

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.