HPV oncogenes expressed from only one of multiple integrated HPV DNA copies drive clonal cell expansion in cervical cancer.

The integration of HPV DNA into human chromosomes plays a pivotal role in the onset of papillomavirus-related cancers. HPV DNA integration often occurs by linearizing the viral DNA in the E1/E2 region, resulting in the loss of a critical viral early polyadenylation signal (PAS), which is essential for the polyadenylation of the E6E7 bicistronic transcripts and for the expression of the viral E6 and E7 oncogenes. Here, we provide compelling evidence that, despite the presence of numerous integrated viral DNA copies, virus-host fusion transcripts originate from only a single integrated HPV DNA in HPV16 and HPV18 cervical cancers and cervical cancer-derived cell lines. The host genomic elements neighboring the integrated HPV DNA are critical for the efficient expression of the viral oncogenes that leads to clonal cell expansion. The fusion RNAs that are produced use a host RNA polyadenylation signal downstream of the integration site, and almost all involve splicing to host sequences. In cell culture, siRNAs specifically targeting the host portion of the virus-host fusion transcripts effectively silenced viral E6 and E7 expression. This, in turn, inhibited cell growth and promoted cell senescence in HPV16+ CaSki and HPV18+ HeLa cells. Showing that HPV E6 and E7 expression from a single integration site is instrumental in clonal cell expansion sheds new light on the mechanisms of HPV-induced carcinogenesis and could be used for the development of precision medicine tailored to combat HPV-related malignancies. IMPORTANCE: Persistent oncogenic HPV infections lead to viral DNA integration into the human genome and the development of cervical, anogenital, and oropharyngeal cancers. The expression of the viral E6 and E7 oncogenes plays a key role in cell transformation and tumorigenesis. However, how E6 and E7 could be expressed from the integrated viral DNA which often lacks a viral polyadenylation signal in the cancer cells remains unknown. By analyzing the integrated HPV DNA sites and expressed HPV RNAs in cervical cancer tissues and cell lines, we show that HPV oncogenes are expressed from only one of multiple chromosomal HPV DNA integrated copies. A host polyadenylation signal downstream of the integrated viral DNA is used for polyadenylation and stabilization of the virus-host chimeric RNAs, making the oncogenic transcripts targetable by siRNAs. This observation provides further understanding of the tumorigenic mechanism of HPV integration and suggests possible therapeutic strategies for the development of precision medicine for HPV cancers.

KSHV promotes oncogenic FOS to inhibit nuclease AEN and transactivate RGS2 for AKT phosphorylation.

Kaposi's sarcoma-associated herpesvirus (KSHV) ORF57 is a lytic RNA-binding protein. We applied BCBL-1 cells in lytic KSHV infection and performed UV cross-linking immunoprecipitation (CLIP) followed by RNA-seq of the CLIPed RNA fragments (CLIP-seq). We identified ORF57-bound transcripts from 544 host protein-coding genes. By comparing with the RNA-seq profiles from BCBL-1 cells with latent and lytic KSHV infection and from HEK293T cells with and without ORF57 expression, we identified FOS and CITED2 RNAs being two common ORF57-specific RNA targets. FOS dimerizes with JUN as a transcription factor AP-1 involved in cell proliferation, differentiation, and transformation. Knockout of the ORF57 gene from the KSHV genome led BAC16-iSLK cells incapable of FOS expression in KSHV lytic infection. The dysfunctional KSHV genome in FOS expression could be rescued by Lenti-ORF57 virus infection. ORF57 protein does not regulate FOS translation but binds to the 13-nt RNA motif near the FOS RNA 5' end and prolongs FOS mRNA half-life 7.7 times longer than it is in the absence of ORF57. This binding of ORF57 to FOS RNA is competitive to the binding of a host nuclease AEN (also referred to as ISG20L1). KSHV infection inhibits the expression of AEN, but not exosomal RNA helicase MTR4. FOS expression mediated by ORF57 inhibits AEN transcription, but transactivates RGS2, a regulator of G-protein coupled receptors. FOS binds a conserved AP-1 site in the RGS2 promoter and enhances RGS2 expression to phosphorylate AKT. Altogether, we have discovered that KSHV ORF57 specifically binds and stabilizes FOS RNA to increase FOS expression, thereby disturbing host gene expression and inducing pathogenesis during KSHV lytic infection.

Genome-wide regulation of KSHV RNA splicing by viral RNA-binding protein ORF57.

RNA splicing plays an essential role in the expression of eukaryotic genes. We previously showed that KSHV ORF57 is a viral splicing factor promoting viral lytic gene expression. In this report, we compared the splicing profile of viral RNAs in BCBL-1 cells carrying a wild-type (WT) versus the cells containing an ORF57 knock-out (57KO) KSHV genome during viral lytic infection. Our analyses of viral RNA splice junctions from RNA-seq identified 269 RNA splicing events in the WT and 255 in the 57KO genome, including the splicing events spanning large parts of the viral genome and the production of vIRF4 circRNAs. No circRNA was detectable from the PAN region. We found that the 57KO alters the RNA splicing efficiency of targeted viral RNAs. Two most susceptible RNAs to ORF57 splicing regulation are the K15 RNA with eight exons and seven introns and the bicistronic RNA encoding both viral thymidylate synthase (ORF70) and membrane-associated E3-ubiquitin ligase (K3). ORF57 inhibits splicing of both K15 introns 1 and 2. ORF70/K3 RNA bears two introns, of which the first intron is within the ORF70 coding region as an alternative intron and the second intron in the intergenic region between the ORF70 and K3 as a constitutive intron. In the WT cells expressing ORF57, most ORF70/K3 transcripts retain the first intron to maintain an intact ORF70 coding region. In contrast, in the 57KO cells, the first intron is substantially spliced out. Using a minigene comprising of ORF70/K3 locus, we further confirmed ORF57 regulation of ORF70/K3 RNA splicing, independently of other viral factors. By monitoring protein expression, we showed that ORF57-mediated retention of the first intron leads to the expression of full-length ORF70 protein. The absence of ORF57 promotes the first intron splicing and expression of K3 protein. Altogether, we conclude that ORF57 regulates alternative splicing of ORF70/K3 bicistronic RNA to control K3-mediated immune evasion and ORF70 participation of viral DNA replication in viral lytic infection.

Inhibition of HSP 90 is associated with potent anti-tumor activity in Papillary Renal Cell Carcinoma.

BACKGROUND: There is no universally accepted treatment for patients with advanced papillary renal cell carcinoma (PRCC). The presence of activating mutations in MET, as well as gain of chromosome 7, where the MET gene is located, are the most common genetic alterations associated with PRCC, leading to the clinical evaluation of MET tyrosine kinase inhibitors (TKIs) in this cancer. However, TKIs targeting MET selectively, as well as multitargeted TKIs with activity against MET demonstrate modest efficacy in PRCC and primary and secondary treatment failure is common; other approaches are urgently needed to improve outcomes in these patients. METHODS: High throughput screening with small molecule libraries identified HSP90 inhibitors as agents of interest based on antitumor activity against patient derived PRCC cell lines. We investigated the activity of the orally available HSP90 inhibitor, SNX2112 in vitro, using 2D/3D PRCC cell culture models and in vivo, in mice tumor xenograft models. The molecular pathways mediating antitumor activity of SNX2112 were assessed by Western blot analysis, Flow cytometry, RNA-seq analysis, Real Time qPCR and imaging approaches. RESULTS: SNX2112 significantly inhibited cellular proliferation, induced G2/M cell cycle arrest and apoptosis in PRCC lines overexpressing MET. In contrast to TKIs targeting MET, SNX2112 inhibited both MET and known downstream mediators of MET activity (AKT, pAKT1/2 and pERK1/2) in PRCC cell lines. RNAi silencing of AKT1/2 or ERK1/2 expression significantly inhibited growth in PRCC cells. Furthermore, SNX2112 inhibited a unique set of E2F and MYC targets and G2M-associated genes. Interestingly, interrogation of the TCGA papillary RCC cohort revealed that these genes were overexpressed in PRCC and portend a poor prognosis. Finally, SNX-2112 demonstrated strong antitumor activity in vivo and prolonged survival of mice bearing human PRCC xenograft. CONCLUSIONS: These results demonstrate that HSP90 inhibition is associated with potent activity in PRCC, and implicate the PI3K/AKT and MEK/ERK1/2 pathways as important mediators of tumorigenesis. These data also provide the impetus for further clinical evaluation of HSP90, AKT, MEK or E2F pathway inhibitors in PRCC.

Ecdysoneless Protein Regulates Viral and Cellular mRNA Splicing to Promote Cervical Oncogenesis.

High-risk human papillomaviruses (HPV), exemplified by HPV16/18, are causally linked to human cancers of the anogenital tract, skin, and upper aerodigestive tract. Previously, we identified Ecdysoneless (ECD) protein, the human homolog of the Drosophila ecdysoneless gene, as a novel HPV16 E6-interacting protein. Here, we show that ECD, through its C-terminal region, selectively binds to high-risk but not to low-risk HPV E6 proteins. We demonstrate that ECD is overexpressed in cervical and head and neck squamous cell carcinoma (HNSCC) cell lines as well as in tumor tissues. Using The Cancer Genome Atlas dataset, we show that ECD mRNA overexpression predicts shorter survival in patients with cervical and HNSCC. We demonstrate that ECD knockdown in cervical cancer cell lines led to impaired oncogenic behavior, and ECD co-overexpression with E7 immortalized primary human keratinocytes. RNA-sequencing analyses of SiHa cells upon ECD knockdown showed to aberrations in E6/E7 RNA splicing, as well as RNA splicing of several HPV oncogenesis-linked cellular genes, including splicing of components of mRNA splicing machinery itself. Taken together, our results support a novel role of ECD in viral and cellular mRNA splicing to support HPV-driven oncogenesis. IMPLICATIONS: This study links ECD overexpression to poor prognosis and shorter survival in HNSCC and cervical cancers and identifies a critical role of ECD in cervical oncogenesis through regulation of viral and cellular mRNA splicing.

Protein-RNA Interactome Analysis Reveals Wide Association of Kaposi's Sarcoma-Associated Herpesvirus ORF57 with Host Noncoding RNAs and Polysomes.

Kaposi's sarcoma-associated herpesvirus (KSHV) ORF57 is an RNA-binding posttranscriptional regulator. We recently applied an affinity-purified anti-ORF57 antibody to conduct ORF57 cross-linking immunoprecipitation (CLIP) in combination with RNA-sequencing (CLIP-seq) and analyzed the genome-wide host RNA transcripts in association with ORF57 in BCBL-1 cells with lytic KSHV infection. Mapping of the CLIP RNA reads to the human genome (GRCh37) revealed that most of the ORF57-associated RNA reads were from rRNAs. The remaining RNA reads mapped to several classes of host noncoding and protein-coding mRNAs. We found that ORF57 binds and regulates expression of a subset of host long noncoding RNAs (lncRNAs), including LINC00324, LINC00355, and LINC00839, which are involved in cell growth. ORF57 binds small nucleolar RNAs (snoRNAs) responsible for 18S and 28S rRNA modifications but does not interact with fibrillarin or NOP58. We validated ORF57 interactions with 67 snoRNAs by ORF57 RNA immunoprecipitation (RIP)-snoRNA array assays. Most of the identified ORF57 rRNA binding sites (BS) overlap the sites binding snoRNAs. We confirmed ORF57-snoRA71B RNA interaction in BCBL-1 cells by ORF57 RIP and Northern blot analyses using a 32P-labeled oligonucleotide probe from the 18S rRNA region complementary to snoRA71B. Using RNA oligonucleotides from the rRNA regions that ORF57 binds for oligonucleotide pulldown-Western blot assays, we selectively verified ORF57 interactions with 5.8S and 18S rRNAs. Polysome profiling revealed that ORF57 associates with both monosomes and polysomes and that its association with polysomes increases PABPC1 binding to polysomes but prevents Ago2 association with polysomes. Our data indicate a functional correlation with ORF57 binding and suppression of Ago2 activities for ORF57 promotion of gene expression. IMPORTANCE As an RNA-binding protein, KSHV ORF57 regulates RNA splicing, stability, and translation and inhibits host innate immunity by blocking the formation of RNA granules in virus-infected cells. In this study, ORF57 was found to interact with many host noncoding RNAs, including lncRNAs, snoRNAs, and rRNAs, to carry out additional unknown functions. ORF57 binds a group of lncRNAs via the RNA motifs identified by ORF57 CLIP-seq to regulate their expression. ORF57 associates with snoRNAs independently of fibrillarin and NOP58 proteins and with rRNA in the regions that commonly bind snoRNAs. Knockdown of fibrillarin expression decreases the expression of snoRNAs and CDK4 but does not affect viral gene expression. More importantly, we found that ORF57 binds translationally active polysomes and enhances PABPC1 but prevents Ago2 association with polysomes. Data provide compelling evidence on how ORF57 in KSHV-infected cells might regulate protein synthesis by blocking Ago2's hostile activities on translation.

Mouse papillomavirus type 1 (MmuPV1) DNA is frequently integrated in benign tumors by microhomology-mediated end-joining.

MmuPV1 is a useful model for studying papillomavirus-induced tumorigenesis. We used RNA-seq to look for chimeric RNAs that map to both MmuPV1 and host genomes. In tumor tissues, a higher proportion of total viral reads were virus-host chimeric junction reads (CJRs) (1.9‰ - 7‰) than in tumor-free tissues (0.6‰ - 1.3‰): most CJRs mapped to the viral E2/E4 region. Although most of the MmuPV1 integration sites were mapped to intergenic regions and introns throughout the mouse genome, integrations were seen more than once in several genes: Malat1, Krt1, Krt10, Fabp5, Pard3, and Grip1; these data were confirmed by rapid amplification of cDNA ends (RACE)-Single Molecule Real-Time (SMRT)-seq or targeted DNA-seq. Microhomology sequences were frequently seen at host-virus DNA junctions. MmuPV1 infection and integration affected the expression of host genes. We found that factors for DNA double-stranded break repair and microhomology-mediated end-joining (MMEJ), such as H2ax, Fen1, DNA polymerase Polθ, Cdk1, and Plk1, exhibited a step-wise increase and Mdc1 a decrease in expression in MmuPV1-infected tissues and MmuPV1 tumors relative to normal tissues. Increased expression of mitotic kinases CDK1 and PLK1 appears to be correlated with CtIP phosphorylation in MmuPV1 tumors, suggesting a role for MMEJ-mediated DNA joining in the MmuPV1 integration events that are associated with MmuPV1-induced progression of tumors.

Stress keratin 17 enhances papillomavirus infection-induced disease by downregulating T cell recruitment.

High-risk human papillomaviruses (HPVs) cause 5% of human cancers. Despite the availability of HPV vaccines, there remains a strong urgency to find ways to treat persistent HPV infections, as current HPV vaccines are not therapeutic for individuals already infected. We used a mouse papillomavirus infection model to characterize virus-host interactions. We found that mouse papillomavirus (MmuPV1) suppresses host immune responses via overexpression of stress keratins. In mice deficient for stress keratin K17 (K17KO), we observed rapid regression of papillomas dependent on T cells. Cellular genes involved in immune response were differentially expressed in the papillomas arising on the K17KO mice correlating with increased numbers of infiltrating CD8+ T cells and upregulation of IFNγ-related genes, including CXCL9 and CXCL10, prior to complete regression. Blocking the receptor for CXCL9/CXCL10 prevented early regression. Our data provide a novel mechanism by which papillomavirus-infected cells evade host immunity and defines new therapeutic targets for treating persistent papillomavirus infections.

CRISPR/Cas9-Mediated Knockout and In Situ Inversion of the ORF57 Gene from All Copies of the Kaposi's Sarcoma-Associated Herpesvirus Genome in BCBL-1 Cells.

Kaposi's sarcoma-associated herpesvirus (KSHV)-transformed primary effusion lymphoma cell lines contain ∼70 to 150 copies of episomal KSHV genomes per cell and have been widely used for studying the mechanisms of KSHV latency and lytic reactivation. Here, we report the first complete knockout (KO) of viral ORF57 gene from all ∼100 copies of KSHV genome per cell in BCBL-1 cells. This was achieved by a modified CRISPR/Cas9 technology to simultaneously express two guide RNAs (gRNAs) and Cas9 from a single expression vector in transfected cells in combination with multiple rounds of cell selection and single-cell cloning. CRISPR/Cas9-mediated genome engineering induces the targeted gene deletion and inversion in situ We found the inverted ORF57 gene in the targeted site in the KSHV genome in one of two characterized single cell clones. Knockout of ORF57 from the KSHV genome led to viral genome instability, thereby reducing viral genome copies and expression of viral lytic genes in BCBL-1-derived single-cell clones. The modified CRISPR/Cas9 technology was very efficient in knocking out the ORF57 gene in iSLK/Bac16 and HEK293/Bac36 cells, where each cell contains only a few copies of the KSHV genome. The ORF57 KO genome was stable in iSLK/Bac16 cells, and, upon lytic induction, was partially rescued by ectopic ORF57 to express viral lytic gene ORF59 and produce infectious virions. Together, the technology developed in this study has paved the way to express two separate gRNAs and the Cas9 enzyme simultaneously in the same cell and could be efficiently applied to any genetic alterations from various genomes, including those in extreme high copy numbers.IMPORTANCE This study provides the first evidence that CRISPR/Cas9 technology can be applied to knock out the ORF57 gene from all ∼100 copies of the KSHV genome in primary effusion lymphoma (PEL) cells by coexpressing two guide RNAs (gRNAs) and Cas9 from a single expression vector in combination with single-cell cloning. The gene knockout efficiency in this system was evaluated rapidly using a direct cell PCR screening. The current CRISPR/Cas9 technology also mediated ORF57 inversion in situ in the targeted site of the KSHV genome. The successful rescue of viral lytic gene expression and infectious virion production from the ORF57 knockout (KO) genome further reiterates the essential role of ORF57 in KSHV infection and multiplication. This modified technology should be useful for knocking out any viral genes from a genome to dissect functions of individual viral genes in the context of the virus genome and to understand their contributions to viral genetics and the virus life cycle.

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.

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.

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.

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.

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.

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.

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.