Iterative epigenomic analyses in the same single cell.

Gene expression in individual cells is epigenetically regulated by DNA modifications, histone modifications, transcription factors, and other DNA-binding proteins. It has been shown that multiple histone modifications can predict gene expression and reflect future responses of bulk cells to extracellular cues. However, the predictive ability of epigenomic analysis is still limited for mechanistic research at a single cell level. To overcome this limitation, it would be useful to acquire reliable signals from multiple epigenetic marks in the same single cell. Here, we propose a new approach and a new method for analysis of several components of the epigenome in the same single cell. The new method allows reanalysis of the same single cell. We found that reanalysis of the same single cell is feasible, provides confirmation of the epigenetic signals, and allows application of statistical analysis to identify reproduced reads using data sets generated only from the single cell. Reanalysis of the same single cell is also useful to acquire multiple epigenetic marks from the same single cells. The method can acquire at least five epigenetic marks: H3K27ac, H3K27me3, mediator complex subunit 1, a DNA modification, and a DNA-interacting protein. We can predict active signaling pathways in K562 single cells using the epigenetic data and confirm that the predicted results strongly correlate with actual active signaling pathways identified by RNA-seq results. These results suggest that the new method provides mechanistic insights for cellular phenotypes through multilayered epigenome analysis in the same single cells.

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.

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.

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.

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.

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.

A spatial map of hepatic mitochondria uncovers functional heterogeneity shaped by nutrient-sensing signaling.

In the liver, mitochondria are exposed to different concentrations of nutrients due to their spatial positioning across the periportal and pericentral axis. How the mitochondria sense and integrate these signals to respond and maintain homeostasis is not known. Here, we combine intravital microscopy, spatial proteomics, and functional assessment to investigate mitochondrial heterogeneity in the context of liver zonation. We find that periportal and pericentral mitochondria are morphologically and functionally distinct; beta-oxidation is elevated in periportal regions, while lipid synthesis is predominant in the pericentral mitochondria. In addition, comparative phosphoproteomics reveals spatially distinct patterns of mitochondrial composition and potential regulation via phosphorylation. Acute pharmacological modulation of nutrient sensing through AMPK and mTOR shifts mitochondrial phenotypes in the periportal and pericentral regions, linking nutrient gradients across the lobule and mitochondrial heterogeneity. This study highlights the role of protein phosphorylation in mitochondrial structure, function, and overall homeostasis in hepatic metabolic zonation. These findings have important implications for liver physiology and disease.

A spatial map of hepatic mitochondria uncovers functional heterogeneity shaped by nutrient-sensing signaling.

In the liver, mitochondria are exposed to different concentrations of nutrients due to their spatial positioning across the periportal (PP) and pericentral (PC) axis. How these mitochondria sense and integrate these signals to respond and maintain homeostasis is not known. Here, we combined intravital microscopy, spatial proteomics, and functional assessment to investigate mitochondrial heterogeneity in the context of liver zonation. We found that PP and PC mitochondria are morphologically and functionally distinct; beta-oxidation was elevated in PP regions, while lipid synthesis was predominant in the PC mitochondria. In addition, comparative phosphoproteomics revealed spatially distinct patterns of mitochondrial composition and potential regulation via phosphorylation. Acute pharmacological modulation of nutrient sensing through AMPK and mTOR shifted mitochondrial phenotypes in the PP and PC regions, linking nutrient gradients across the lobule and mitochondrial heterogeneity. This study highlights the role of protein phosphorylation in mitochondrial structure, function, and overall homeostasis in hepatic metabolic zonation. These findings have important implications for liver physiology and disease.

HIV vaccine candidate efficacy in female macaques mediated by cAMP-dependent efferocytosis and V2-specific ADCC.

The development of an effective vaccine to protect against HIV acquisition will be greatly bolstered by in-depth understanding of the innate and adaptive responses to vaccination. We report here that the efficacy of DNA/ALVAC/gp120/alum vaccines, based on V2-specific antibodies mediating apoptosis of infected cells (V2-ADCC), is complemented by efferocytosis, a cyclic AMP (cAMP)-dependent antiphlogistic engulfment of apoptotic cells by CD14+ monocytes. Central to vaccine efficacy is the engagement of the CCL2/CCR2 axis and tolerogenic dendritic cells producing IL-10 (DC-10). Epigenetic reprogramming in CD14+ cells of the cyclic AMP/CREB pathway and increased systemic levels of miRNA-139-5p, a negative regulator of expression of the cAMP-specific phosphodiesterase PDE4D, correlated with vaccine efficacy. These data posit that efferocytosis, through the prompt and effective removal of apoptotic infected cells, contributes to vaccine efficacy by decreasing inflammation and maintaining tissue homeostasis.

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.