UVB-Induced Tumor Heterogeneity Diminishes Immune Response in Melanoma.

Although clonal neo-antigen burden is associated with improved response to immune therapy, the functional basis for this remains unclear. Here we study this question in a novel controlled mouse melanoma model that enables us to explore the effects of intra-tumor heterogeneity (ITH) on tumor aggressiveness and immunity independent of tumor mutational burden. Induction of UVB-derived mutations yields highly aggressive tumors with decreased anti-tumor activity. However, single-cell-derived tumors with reduced ITH are swiftly rejected. Their rejection is accompanied by increased T cell reactivity and a less suppressive microenvironment. Using phylogenetic analyses and mixing experiments of single-cell clones, we dissect two characteristics of ITH: the number of clones forming the tumor and their clonal diversity. Our analysis of melanoma patient tumor data recapitulates our results in terms of overall survival and response to immune checkpoint therapy. These findings highlight the importance of clonal mutations in robust immune surveillance and the need to quantify patient ITH to determine the response to checkpoint blockade.

Anti-tumour immunity induces aberrant peptide presentation in melanoma.

Extensive tumour inflammation, which is reflected by high levels of infiltrating T cells and interferon-γ (IFNγ) signalling, improves the response of patients with melanoma to checkpoint immunotherapy1,2. Many tumours, however, escape by activating cellular pathways that lead to immunosuppression. One such mechanism is the production of tryptophan metabolites along the kynurenine pathway by the enzyme indoleamine 2,3-dioxygenase 1 (IDO1), which is induced by IFNγ3-5. However, clinical trials using inhibition of IDO1 in combination with blockade of the PD1 pathway in patients with melanoma did not improve the efficacy of treatment compared to PD1 pathway blockade alone6,7, pointing to an incomplete understanding of the role of IDO1 and the consequent degradation of tryptophan in mRNA translation and cancer progression. Here we used ribosome profiling in melanoma cells to investigate the effects of prolonged IFNγ treatment on mRNA translation. Notably, we observed accumulations of ribosomes downstream of tryptophan codons, along with their expected stalling at the tryptophan codon. This suggested that ribosomes bypass tryptophan codons in the absence of tryptophan. A detailed examination of these tryptophan-associated accumulations of ribosomes-which we term 'W-bumps'-showed that they were characterized by ribosomal frameshifting events. Consistently, reporter assays combined with proteomic and immunopeptidomic analyses demonstrated the induction of ribosomal frameshifting, and the generation and presentation of aberrant trans-frame peptides at the cell surface after treatment with IFNγ. Priming of naive T cells from healthy donors with aberrant peptides induced peptide-specific T cells. Together, our results suggest that IDO1-mediated depletion of tryptophan, which is induced by IFNγ, has a role in the immune recognition of melanoma cells by contributing to diversification of the peptidome landscape.

Alternative polyadenylation factor CPSF6 regulates temperature compensation of the mammalian circadian clock.

A defining property of circadian clocks is temperature compensation, characterized by the resilience of their near 24-hour free-running periods against changes in environmental temperature within the physiological range. While temperature compensation is evolutionary conserved across different taxa of life and has been studied within many model organisms, its molecular underpinnings remain elusive. Posttranscriptional regulations such as temperature-sensitive alternative splicing or phosphorylation have been described as underlying reactions. Here, we show that knockdown of cleavage and polyadenylation specificity factor subunit 6 (CPSF6), a key regulator of 3'-end cleavage and polyadenylation, significantly alters circadian temperature compensation in human U-2 OS cells. We apply a combination of 3'-end-RNA-seq and mass spectrometry-based proteomics to globally quantify changes in 3' UTR length as well as gene and protein expression between wild-type and CPSF6 knockdown cells and their dependency on temperature. Since changes in temperature compensation behavior should be reflected in alterations of temperature responses within one or all of the 3 regulatory layers, we statistically assess differential responses upon changes in ambient temperature between wild-type and CPSF6 knockdown cells. By this means, we reveal candidate genes underlying circadian temperature compensation, including eukaryotic translation initiation factor 2 subunit 1 (EIF2S1).

Translation of CircRNAs.

Circular RNAs (circRNAs) are abundant and evolutionarily conserved RNAs of largely unknown function. Here, we show that a subset of circRNAs is translated in vivo. By performing ribosome footprinting from fly heads, we demonstrate that a group of circRNAs is associated with translating ribosomes. Many of these ribo-circRNAs use the start codon of the hosting mRNA, are bound by membrane-associated ribosomes, and have evolutionarily conserved termination codons. In addition, we found that a circRNA generated from the muscleblind locus encodes a protein, which we detected in fly head extracts by mass spectrometry. Next, by performing in vivo and in vitro translation assays, we show that UTRs of ribo-circRNAs (cUTRs) allow cap-independent translation. Moreover, we found that starvation and FOXO likely regulate the translation of a circMbl isoform. Altogether, our study provides strong evidence for translation of circRNAs, revealing the existence of an unexplored layer of gene activity.

Large-Scale Immunopeptidome Analysis Reveals Recurrent Posttranslational Splicing of Cancer- and Immune-Associated Genes.

Posttranslational spliced peptides (PTSPs) are a unique class of peptides that have been found to be presented by HLA class-I molecules in cancer. Thus far, no consensus has been reached on the proportion of PTSPs in the immunopeptidome, with estimates ranging from 2% to as high as 45% and stirring significant debate. Furthermore, the role of the HLA class-II pathway in PTSP presentation has been studied only in diabetes. Here, we exploit our large-scale cancer peptidomics database and our newly devised pipeline for filtering spliced peptide predictions to identify recurring spliced peptides, both for HLA class-I and class-II complexes. Our results indicate that HLA class-I-spliced peptides account for a low percentage of the immunopeptidome (less than 3.1%) yet are larger in number relative to other types of identified aberrant peptides. Therefore, spliced peptides significantly contribute to the repertoire of presented peptides in cancer cells. In addition, we identified HLA class-II-bound spliced peptides, but to a lower extent (less than 0.5%). The identified spliced peptides include cancer- and immune-associated genes, such as the MITF oncogene, DAPK1 tumor suppressor, and HLA-E, which were validated using synthetic peptides. The potential immunogenicity of the DAPK1- and HLA-E-derived PTSPs was also confirmed. In addition, a reanalysis of our published mouse single-cell clone immunopeptidome dataset showed that most of the spliced peptides were found repeatedly in a large number of the single-cell clones. Establishing a novel search-scheme for the discovery and evaluation of recurring PTSPs among cancer patients may assist in identifying potential novel targets for immunotherapy.

Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed.

Circular RNAs (circRNAs) are an endogenous class of animal RNAs. Despite their abundance, their function and expression in the nervous system are unknown. Therefore, we sequenced RNA from different brain regions, primary neurons, isolated synapses, as well as during neuronal differentiation. Using these and other available data, we discovered and analyzed thousands of neuronal human and mouse circRNAs. circRNAs were extraordinarily enriched in the mammalian brain, well conserved in sequence, often expressed as circRNAs in both human and mouse, and sometimes even detected in Drosophila brains. circRNAs were overall upregulated during neuronal differentiation, highly enriched in synapses, and often differentially expressed compared to their mRNA isoforms. circRNA expression correlated negatively with expression of the RNA-editing enzyme ADAR1. Knockdown of ADAR1 induced elevated circRNA expression. Together, we provide a circRNA brain expression atlas and evidence for important circRNA functions and values as biomarkers.

circRNA biogenesis competes with pre-mRNA splicing.

Circular RNAs (circRNAs) are widely expressed noncoding RNAs. However, their biogenesis and possible functions are poorly understood. Here, by studying circRNAs that we identified in neuronal tissues, we provide evidence that animal circRNAs are generated cotranscriptionally and that their production rate is mainly determined by intronic sequences. We demonstrate that circularization and splicing compete against each other. These mechanisms are tissue specific and conserved in animals. Interestingly, we observed that the second exon of the splicing factor muscleblind (MBL/MBNL1) is circularized in flies and humans. This circRNA (circMbl) and its flanking introns contain conserved muscleblind binding sites, which are strongly and specifically bound by MBL. Modulation of MBL levels strongly affects circMbl biosynthesis, and this effect is dependent on the MBL binding sites. Together, our data suggest that circRNAs can function in gene regulation by competing with linear splicing. Furthermore, we identified muscleblind as a factor involved in circRNA biogenesis.

PhISCS-BnB: a fast branch and bound algorithm for the perfect tumor phylogeny reconstruction problem.

MOTIVATION: Recent advances in single-cell sequencing (SCS) offer an unprecedented insight into tumor emergence and evolution. Principled approaches to tumor phylogeny reconstruction via SCS data are typically based on general computational methods for solving an integer linear program, or a constraint satisfaction program, which, although guaranteeing convergence to the most likely solution, are very slow. Others based on Monte Carlo Markov Chain or alternative heuristics not only offer no such guarantee, but also are not faster in practice. As a result, novel methods that can scale up to handle the size and noise characteristics of emerging SCS data are highly desirable to fully utilize this technology. RESULTS: We introduce PhISCS-BnB (phylogeny inference using SCS via branch and bound), a branch and bound algorithm to compute the most likely perfect phylogeny on an input genotype matrix extracted from an SCS dataset. PhISCS-BnB not only offers an optimality guarantee, but is also 10-100 times faster than the best available methods on simulated tumor SCS data. We also applied PhISCS-BnB on a recently published large melanoma dataset derived from the sublineages of a cell line involving 20 clones with 2367 mutations, which returned the optimal tumor phylogeny in <4 h. The resulting phylogeny agrees with and extends the published results by providing a more detailed picture on the clonal evolution of the tumor. AVAILABILITY AND IMPLEMENTATION: https://github.com/algo-cancer/PhISCS-BnB. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic hyper-editing underlies temperature adaptation in Drosophila.

In Drosophila, A-to-I editing is prevalent in the brain, and mutations in the editing enzyme ADAR correlate with specific behavioral defects. Here we demonstrate a role for ADAR in behavioral temperature adaptation in Drosophila. Although there is a higher level of editing at lower temperatures, at 29°C more sites are edited. These sites are less evolutionarily conserved, more disperse, less likely to be involved in secondary structures, and more likely to be located in exons. Interestingly, hypomorph mutants for ADAR display a weaker transcriptional response to temperature changes than wild-type flies and a highly abnormal behavioral response upon temperature increase. In sum, our data shows that ADAR is essential for proper temperature adaptation, a key behavior trait that is essential for survival of flies in the wild. Moreover, our results suggest a more general role of ADAR in regulating RNA secondary structures in vivo.

The transcription factor Cabut coordinates energy metabolism and the circadian clock in response to sugar sensing.

Nutrient sensing pathways adjust metabolism and physiological functions in response to food intake. For example, sugar feeding promotes lipogenesis by activating glycolytic and lipogenic genes through the Mondo/ChREBP-Mlx transcription factor complex. Concomitantly, other metabolic routes are inhibited, but the mechanisms of transcriptional repression upon sugar sensing have remained elusive. Here, we characterize cabut (cbt), a transcription factor responsible for the repressive branch of the sugar sensing transcriptional network in Drosophila. We demonstrate that cbt is rapidly induced upon sugar feeding through direct regulation by Mondo-Mlx. We found that CBT represses several metabolic targets in response to sugar feeding, including both isoforms of phosphoenolpyruvate carboxykinase (pepck). Deregulation of pepck1 (CG17725) in mlx mutants underlies imbalance of glycerol and glucose metabolism as well as developmental lethality. Furthermore, we demonstrate that cbt provides a regulatory link between nutrient sensing and the circadian clock. Specifically, we show that a subset of genes regulated by the circadian clock are also targets of CBT. Moreover, perturbation of CBT levels leads to deregulation of the circadian transcriptome and circadian behavioral patterns.

Defining the 5΄ and 3΄ landscape of the Drosophila transcriptome with Exo-seq and RNaseH-seq.

Cells regulate biological responses in part through changes in transcription start sites (TSS) or cleavage and polyadenylation sites (PAS). To fully understand gene regulatory networks, it is therefore critical to accurately annotate cell type-specific TSS and PAS. Here we present a simple and straightforward approach for genome-wide annotation of 5΄- and 3΄-RNA ends. Our approach reliably discerns bona fide PAS from false PAS that arise due to internal poly(A) tracts, a common problem with current PAS annotation methods. We applied our methodology to study the impact of temperature on the Drosophila melanogaster head transcriptome. We found hundreds of previously unidentified TSS and PAS which revealed two interesting phenomena: first, genes with multiple PASs tend to harbor a motif near the most proximal PAS, which likely represents a new cleavage and polyadenylation signal. Second, motif analysis of promoters of genes affected by temperature suggested that boundary element association factor of 32 kDa (BEAF-32) and DREF mediates a transcriptional program at warm temperatures, a result we validated in a fly line where beaf-32 is downregulated. These results demonstrate the utility of a high-throughput platform for complete experimental and computational analysis of mRNA-ends to improve gene annotation.

Synergistic interactions between the molecular and neuronal circadian networks drive robust behavioral circadian rhythms in Drosophila melanogaster.

Most organisms use 24-hr circadian clocks to keep temporal order and anticipate daily environmental changes. In Drosophila melanogaster CLOCK (CLK) and CYCLE (CYC) initiates the circadian system by promoting rhythmic transcription of hundreds of genes. However, it is still not clear whether high amplitude transcriptional oscillations are essential for circadian timekeeping. In order to address this issue, we generated flies in which the amplitude of CLK-driven transcription can be reduced partially (approx. 60%) or strongly (90%) without affecting the average levels of CLK-target genes. The impaired transcriptional oscillations lead to low amplitude protein oscillations that were not sufficient to drive outputs of peripheral oscillators. However, circadian rhythms in locomotor activity were resistant to partial reduction in transcriptional and protein oscillations. We found that the resilience of the brain oscillator is depending on the neuronal communication among circadian neurons in the brain. Indeed, the capacity of the brain oscillator to overcome low amplitude transcriptional oscillations depends on the action of the neuropeptide PDF and on the pdf-expressing cells having equal or higher amplitude of molecular rhythms than the rest of the circadian neuronal groups in the fly brain. Therefore, our work reveals the importance of high amplitude transcriptional oscillations for cell-autonomous circadian timekeeping. Moreover, we demonstrate that the circadian neuronal network is an essential buffering system that protects against changes in circadian transcription in the brain.

The p53-Mdm2 loop: a critical juncture of stress response.

The presence of a functional p53 protein is a key factor for the proper suppression of cancer development. A loss of p53 activity, by mutations or inhibition, is often associated with human malignancies. The p53 protein integrates various stress signals into a growth restrictive cellular response. In this way, p53 eliminates cells with a potential to become cancerous. Being a powerful decision maker, it is imperative that p53 will be activated properly, efficiently and temporarily in response to stress. Equally important is that p53 activation will be extinguished upon recovery from stress, and that improper activation of p53 will be avoided. Failure to achieve these aims is likely to have catastrophic consequences for the organism. The machinery that governs this tight regulation is largely based on the major inhibitor of p53, Mdm2, which both blocks p53 activities and promotes its destabilization. The interplay between p53 and Mdm2 involves a complex network of positive and negative feedback loops. Relief from Mdm2 suppression is required for p53 to be stabilized and activated in response to stress. Protection from Mdm2 entails a concerted action of modifying enzymes and partner proteins. The association of p53 with the PML-nuclear bodies may provide an infrastructure in which this complex regulatory network can be orchestrated. In this chapter we use examples to illustrate the regulatory machinery that drives this network.

E6AP ubiquitin ligase regulates PML-induced senescence in Myc-driven lymphomagenesis.

Neoplastic transformation requires the elimination of key tumor suppressors, which may result from E3 ligase-mediated proteasomal degradation. We previously demonstrated a key role for the E3 ubiquitin ligase E6AP in the regulation of promyelocytic leukemia protein (PML) stability and formation of PML nuclear bodies. Here, we report the involvement of the E6AP-PML axis in B-cell lymphoma development. A partial loss of E6AP attenuated Myc-induced B-cell lymphomagenesis. This tumor suppressive action was achieved by the induction of cellular senescence. B-cell lymphomas deficient for E6AP expressed elevated levels of PML and PML-nuclear bodies with a concomitant increase in markers of cellular senescence, including p21, H3K9me3, and p16. Consistently, PML deficiency accelerated the rate of Myc-induced B-cell lymphomagenesis. Importantly, E6AP expression was elevated in ∼ 60% of human Burkitt lymphomas, and down-regulation of E6AP in B-lymphoma cells restored PML expression with a concurrent induction of cellular senescence in these cells. Our findings demonstrate that E6AP-mediated down-regulation of PML-induced senescence is essential for B-cell lymphoma progression. This provides a molecular explanation for the down-regulation of PML observed in non-Hodgkin lymphomas, thereby suggesting a novel therapeutic approach for restoration of tumor suppression in B-cell lymphoma.

Adaptation of molecular circadian clockwork to environmental changes: a role for alternative splicing and miRNAs.

Circadian (24 h) clocks provide a source of internal timing in most living organisms. These clocks keep time by using complex transcriptional/post-translational feedback loops that are strikingly resilient to changes in environmental conditions. In the last few years, interest has increased in the role of post-transcriptional regulation of circadian clock components. Post-transcriptional control plays a prominent role in modulating rapid responses of the circadian system to environmental changes, including light, temperature and general stress and will be the focus of this review.

Temporal genomic analysis of melanoma rejection identifies regulators of tumor immune evasion.

Decreased intra-tumor heterogeneity (ITH) correlates with increased patient survival and immunotherapy response. However, even highly homogenous tumors may display variability in their aggressiveness, and how immunologic-factors impinge on their aggressiveness remains understudied. Here we studied the mechanisms responsible for the immune-escape of murine tumors with low ITH. We compared the temporal growth of homogeneous, genetically-similar single-cell clones that are rejected vs. those that are not-rejected after transplantation in-vivo using single-cell RNA sequencing and immunophenotyping. Non-rejected clones showed high infiltration of tumor-associated-macrophages (TAMs), lower T-cell infiltration, and increased T-cell exhaustion compared to rejected clones. Comparative analysis of rejection-associated gene expression programs, combined with in-vivo CRISPR knockout screens of candidate mediators, identified Mif (macrophage migration inhibitory factor) as a regulator of immune rejection. Mif knockout led to smaller tumors and reversed non-rejection-associated immune composition, particularly, leading to the reduction of immunosuppressive macrophage infiltration. Finally, we validated these results in melanoma patient data.

MDM2 and Fbw7 cooperate to induce p63 protein degradation following DNA damage and cell differentiation.

Tight control of p63 protein levels must be achieved under differentiation or apoptotic conditions. Here, we describe a new regulatory pathway for the DeltaNp63alpha protein. We found that MDM2 binds DeltaNp63alpha in the nucleus promoting its translocation to the cytoplasm. The MDM2 nuclear localization signal is required for DeltaNp63alpha nuclear export and subsequent degradation, whereas the MDM2 ring-finger domain is dispensable. Once exported to the cytoplasm by MDM2, p63 is targeted for degradation by the Fbw7 E3-ubiquitin ligase. Efficient degradation of DeltaNp63alpha by Fbw7 (also known as FBXW7) requires GSK3 kinase activity. By deletion and point mutations analysis we have identified a phosphodegron located in the alpha and beta tail of p63 that is required for degradation. Furthermore, we show that MDM2 or Fbw7 depletion inhibits degradation of endogenous DeltaNp63alpha in cells exposed to UV irradiation, adriamycin and upon keratinocyte differentiation. Our findings suggest that following DNA damage and cellular differentiation MDM2 and Fbw7 can cooperate to regulate the levels of the pro-proliferative DeltaNp63alpha protein.

The association of MHC class I proteins with the 2B4 receptor inhibits self-killing of human NK cells.

The killing activity of NK cells is carried out by several activating NK receptors, which includes NKp46, NKp44, NKp30, NKp80, NKG2D, and 2B4. The ligands of these receptors are either self-derived, pathogen-derived, stress-induced ligands or tumor ligands. Importantly, none of these killer ligands are expressed on NK cells and thus self-killing of NK cells is prevented. A notable exception with this regard, is the ligand of the 2B4 receptor. This unusual receptor can exert both activating and inhibiting signals; however, in human NK cells, it serves mainly as an activating receptor. The ligand of 2B4 is CD48 and in contrast to the ligands of all the other NK activating receptors, CD48 is also present on NK cells. Thus, NK cells might be at risk for self-killing that is mediated via the 2B4-CD48 interaction. In this study, we identify a novel mechanism that prevents this self-killing as we show that the association of the MHC class I proteins with the 2B4 receptor, both present on NK cells, results in the attenuation of the 2B4-mediated self-killing of NK cells.

circMbl functions in cis and in trans to regulate gene expression and physiology in a tissue-specific fashion.

Muscleblind (mbl) is an essential muscle and neuronal splicing regulator. Mbl hosts multiple circular RNAs (circRNAs), including circMbl, which is conserved from flies to humans. Here, we show that mbl-derived circRNAs are key regulators of MBL by cis- and trans-acting mechanisms. By generating fly lines to specifically modulate the levels of all mbl RNA isoforms, including circMbl, we demonstrate that the two major mbl protein isoforms, MBL-O/P and MBL-C, buffer their own levels by producing different types of circRNA isoforms in the eye and fly brain, respectively. Moreover, we show that circMbl has unique functions in trans, as knockdown of circMbl results in specific morphological and physiological phenotypes. In addition, depletion of MBL-C or circMbl results in opposite behavioral phenotypes, showing that they also regulate each other in trans. Together, our results illuminate key aspects of mbl regulation and uncover cis and trans functions of circMbl in vivo.

An in vivo strategy for knockdown of circular RNAs.

Exonic circular RNAs (circRNAs) are highly abundant RNAs generated mostly from exons of protein-coding genes. Assaying the functions of circRNAs is not straightforward as common approaches for circRNA depletion tend to also alter the levels of mRNAs generated from the hosting gene. Here we describe a methodology for specific knockdown of circRNAs in vivo with tissue and cell resolution. We also describe an experimental and computational platform for determining the potential off-target effects as well as for verifying the obtained phenotypes. Briefly, we utilize shRNAs targeted to the circRNA-specific back-splice junction to specifically downregulate the circRNA. We utilized this methodology to downregulate five circRNAs that are highly expressed in Drosophila. There were no effects on the levels of their linear counterparts or any RNA with complementarity to the expressed shRNA. Interestingly, downregulation of circCtrip resulted in developmental lethality that was recapitulated with a second shRNA. Moreover, downregulation of individual circRNAs caused specific changes in the fly head transcriptome, suggesting roles for these circRNAs in the fly nervous system. Together, our results provide a methodological approach that enables the comprehensive study of circRNAs at the organismal and cellular levels and generated for the first time flies in which specific circRNAs are downregulated.