Intestinal stem cells (ISCs) face the challenge of integrating metabolic demands with unique regenerative functions. Studies have shown an intricate interplay between metabolism and stem cell capacity; however, it is still not understood how this process is regulated. Combining ribosome profiling and CRISPR screening in intestinal organoids, we identify the nascent polypeptide–associated complex (NAC) as a key mediator of this process. Our findings suggest that NAC is responsible for relocalizing ribosomes to the mitochondria and regulating ISC metabolism. Upon NAC inhibition, intestinal cells show decreased import of mitochondrial proteins, which are needed for oxidative phosphorylation, and, consequently, enable the cell to maintain a stem cell identity. Furthermore, we show that overexpression of NACα is sufficient to drive mitochondrial respiration and promote ISC identity. Ultimately, our results reveal the pivotal role of NAC in regulating ribosome localization, mitochondrial metabolism, and ISC function, providing insights into the potential mechanism behind it.

Recent advances in the fields of RNA translation and ribosome biology have demonstrated the heterogeneous nature of ribosomes. This manifests not only across different cellular contexts but also within the same cell. Such variations in ribosomal composition, be it in ribosomal RNAs or proteins, can significantly influence cellular processes and responses by altering the mRNAs being translated or the dynamics of ribosomes during the translation process. Therefore, identifying this heterogeneity is crucial for unravelling the complexity of gene expression across different fields of biology. Here we provide an overview of recent advances in high-throughput techniques for identifying ribosomal heterogeneity. We cover methodologies for probing both rRNA and protein components of the ribosome and encompass the most recent next-generation sequencing and computational analyses, as well as a diverse array of mass spectrometry techniques.This article is part of the discussion meeting issue 'Ribosome diversity and its impact on protein synthesis, development and disease'.

Intestinal stem cells (ISCs) face the challenge of integrating metabolic demands with unique regenerative functions. Studies have shown an intricate interplay between metabolism and stem cell capacity, however it is still not understood how this process is regulated. Combining ribosome profiling and CRISPR screening in intestinal organoids, we show that RNA translation is at the root of this interplay. We identify the nascent polypeptide-associated complex (NAC) as a key mediator of this process, and show that it regulates ISC metabolism by relocalizing ribosomes to the mitochondria. Upon NAC inhibition, intestinal cells show decreased import of mitochondrial proteins, which are needed for oxidative phosphorylation, and, consequently, enable the cell to maintain a stem cell identity. Furthermore, we show that overexpression of NACα is sufficient to drive mitochondrial respiration and promote ISC identity. Ultimately, our results reveal the pivotal role of ribosome localization in regulating mitochondrial metabolism and ISC function.Teaser The location of ribosomes in cells is regulated, and defines the fate of intestinal stem cells.### Competing Interest StatementThe authors have declared no competing interest.

Glioblastoma (GBM) is a fatal brain tumor with a critical need for better therapies. It is known that the PI3K, MAPK, and CDK4/6 signaling pathways are hyper-activated in these tumors; however, previous studies have used very high concentration of inhibitors to assess their importance, with mixed results. Here we developed PMCi, a combination approach that targets all three pathways simultaneously, at clinically-relevant doses. PMCi effectively suppresses GBM cell proliferation in vitro and in vivo, and outperforms monotherapies and dual combinations. PMCi acts by inducing cellular senescence, which is mediated solely by the mitochondria, and, unlike other forms of senescence, is independent of nuclear damage. This phenotype is caused by a reactive oxygen species (ROS)\cGAS-STING\senescence-associated secretory phenotype (SASP) signaling cascade, that acts in a paracrine manner to establish and maintain senescence. Our results demonstrate that mitochondrial damage is sufficient to drive senescence, and that this can be leveraged to target GBM cells. ### Competing Interest Statement The authors have declared no competing interest.

mRNA transcripts have limited potential for protein synthesis, defined by their open reading frames[1][1]. However, recent advances have revealed a far more complex reality, in which the proteome exceeds the perceived limits of the transcriptome[2][2]–[6][3] exposing a significant gap in our understanding. Our prior studies have demonstrated that mRNA can undergo further processing, yielding truncated, uncapped mRNAs with translation potential[2][2],[7][4]. Yet, the biological importance of this process remains mostly unclear. Here, we demonstrate that cleavage within the coding sequence of JAK1 mRNA produces an uncapped variant downstream to the cleavage site, at the expense of the full-length transcript. This results in the independent translation of the JH1 kinase domain, a process we term ‘RNA dicing’. Notably, canonical and diced JAK1 variants have distinct impacts on cell proliferation and tumorigenesis, operating independently, localizing to distinct cellular compartments. In addition, the activation of JAK1 signaling through IFNγ induction promotes the dicing of JAK1, thereby altering the balance of isoforms present. Base editor screens reveal that stop-codon nonsense mutations, which are typically considered loss-of-function, have differing impacts depending on their position relative to the dicing site. In agreement with this, JAK1 nonsense frameshift (JAK1fs) mutations in endometrial tumors inhibit the tumor-suppressive functions of canonical JAK1, while amplifying the oncogenic potential of the diced JH1 kinase domain. Lastly, we demonstrate that Momelotinib[8][5], a JAK1-specific inhibitor, is more effective in cancer cells carrying a “loss-of-function” JAK1fs mutation, highlighting the significant impact of RNA dicing biology as a potent tool for patient stratification. Our findings characterize RNA dicing as a fundamental regulatory machinery that diversifies the potential products of a single mRNA molecule, and allows for significant variation in biological function.### Competing Interest StatementThe authors have declared no competing interest. [1]: #ref-1 [2]: #ref-2 [3]: #ref-6 [4]: #ref-7 [5]: #ref-8

Inflammatory cytokines are pivotal to immune responses. Upon cytokine exposure, cells enter an "alert state" that enhances their visibility to the immune system. Here, we identified an alert-state subpopulation of ribosomes defined by the presence of the P-stalk. We show that P-stalk ribosomes (PSRs) are formed in response to cytokines linked to tumor immunity, and this is at least partially mediated by P-stalk phosphorylation. PSRs are involved in the preferential translation of mRNAs vital for the cytokine response via the more efficient translation of transmembrane domains of receptor molecules involved in cytokine-mediated processes. Importantly, loss of the PSR inhibits CD8+ T cell recognition and killing, and inhibitory cytokines like transforming growth factor β (TGF-β) hinder PSR formation, suggesting that the PSR is a central regulatory hub upon which multiple signals converge. Thus, the PSR is an essential mediator of the cellular rewiring that occurs following cytokine exposure via the translational regulation of this process.

Genes limiting T cell antitumor activity may serve as therapeutic targets. It has not been systematically studied whether there are regulators that uniquely or broadly contribute to T cell fitness. We perform genome-scale CRISPR-Cas9 knockout screens in primary CD8 T cells to uncover genes negatively impacting fitness upon three modes of stimulation: (1) intense, triggering activation-induced cell death (AICD); (2) acute, triggering expansion; (3) chronic, causing dysfunction. Besides established regulators, we uncover genes controlling T cell fitness either specifically or commonly upon differential stimulation. Dap5 ablation, ranking highly in all three screens, increases translation while enhancing tumor killing. Loss of Icam1-mediated homotypic T cell clustering amplifies cell expansion and effector functions after both acute and intense stimulation. Lastly, Ctbp1 inactivation induces functional T cell persistence exclusively upon chronic stimulation. Our results functionally annotate fitness regulators based on their unique or shared contribution to traits limiting T cell antitumor activity.

Inflammatory cytokines are pivotal to immune responses. Upon cytokine exposure, cells enter an “alert-state” that enhances their visibility to the immune system. Here, we identified an “alert-state” subpopulation of ribosomes (ASRs) defined by the presence of the P-stalk. We show that ASRs are formed in response to cytokines linked to tumor immunity, and are involved in the preferential translation of mRNAs vital for the cytokine response. Mechanistically, ASRs are required for the efficient translation of transmembrane domains of receptor molecules involved in cytokine-mediated processes. Importantly, loss of the ASR prevents CD8+ T cell recognition and killing, and inhibitory cytokines like TGFβ hinder ASR formation, suggesting that the ASR is a central regulatory hub upon which multiple signals converge. Thus, the ASR is an essential mediator of the cellular rewiring that occurs following cytokine exposure, via the translational regulation of this process.

AbstractYield improvements in cell factories can potentially be obtained by fine-tuning the regulatory mechanisms for gene candidates. In pursuit of such candidates, we performed RNA-sequencing of two α-amylase producing Bacillus strains and predict hundreds of putative novel non-coding transcribed regions. Surprisingly, we found among hundreds of non-coding and structured RNA candidates that non-coding genomic regions are proportionally undergoing the highest changes in expression during fermentation. Since these classes of RNA are also understudied, we targeted the corresponding genomic regions with CRIPSRi knockdown to test for any potential impact on the yield. From differentially expression analysis, we selected 53 non-coding candidates. Although CRISPRi knockdowns target both the sense and the antisense strand, the CRISPRi experiment cannot link causes for yield changes to the sense or antisense disruption. Nevertheless, we observed on several instances with strong changes in enzyme yield. The knockdown targeting the genomic region for a putative antisense RNA of the 3′ UTR of the skfA-skfH operon led to a 21% increase in yield. In contrast, the knockdown targeting the genomic regions of putative antisense RNAs of the cytochrome c oxidase subunit 1 (ctaD), the sigma factor sigH, and the uncharacterized gene yhfT decreased yields by 31 to 43%.

A major challenge of CRISPR/Cas9-mediated genome engineering is that not all guide RNAs (gRNAs) cleave the DNA efficiently. Although the heterogeneity of gRNA activity is well recognized, the current understanding of how CRISPR/Cas9 activity is regulated remains incomplete. Here, we identify a sweet spot range of binding free energy change for optimal efficiency which largely explains why gRNAs display changes in efficiency at on- and off-target sites, including why gRNAs can cleave an off-target with higher efficiency than the on-target. Using an energy-based model, we show that local gRNA-DNA interactions resulting from Cas9 “sliding” on overlapping protospacer adjacent motifs (PAMs) profoundly impact gRNA activities. Combining the effects of local sliding for a given PAM context with global off-targets allows us to better identify highly specific, and thus efficient, gRNAs. We validate the effects of local sliding on gRNA efficiency using both public data and in-house data generated by measuring SpCas9 cleavage efficiency at 1024 sites designed to cover all possible combinations of 4-nt PAM and context sequences of 4 gRNAs. Our results provide insights into the mechanisms of Cas9-PAM compatibility and cleavage activation, underlining the importance of accounting for local sliding in gRNA design.

Recent studies have revealed multiple mechanisms that can lead to heterogeneity in ribosomal composition. This heterogeneity can lead to preferential translation of specific panels of mRNAs, and is defined in large part by the ribosomal protein (RP) content, amongst other things. However, it is currently unknown to what extent ribosomal composition is heterogeneous across tissues, which is compounded by a lack of tools available to study it. Here we present dripARF, a method for detecting differential RP incorporation into the ribosome using Ribosome Profiling (Ribo-seq) data. We combine the 'waste' rRNA fragment data generated in Ribo-seq with the known 3D structure of the human ribosome to predict differences in the composition of ribosomes in the material being studied. We have validated this approach using publicly available data, and have revealed a potential role for eS25/RPS25 in development. Our results indicate that ribosome heterogeneity can be detected in Ribo-seq data, providing a new method to study this phenomenon. Furthermore, with dripARF, previously published Ribo-seq data provides a wealth of new information, allowing the identification of RPs of interest in many disease and normal contexts. dripARF is available as part of the ARF R package and can be accessed through https://github.com/fallerlab/ARF.

The use of alternative promoters, splicing, and cleavage and polyadenylation (APA) generates mRNA isoforms that expand the diversity and complexity of the transcriptome. Here, we uncovered thousands of previously undescribed 5' uncapped and polyadenylated transcripts (5' UPTs). We show that these transcripts resist exonucleases due to a highly structured RNA and N6-methyladenosine modification at their 5' termini. 5' UPTs appear downstream of APA sites within their host genes and are induced upon APA activation. Strong enrichment in polysomal RNA fractions indicates 5' UPT translational potential. Indeed, APA promotes downstream translation initiation, non-canonical protein output, and consistent changes to peptide presentation at the cell surface. Lastly, we demonstrate the biological importance of 5' UPTs using Bcl2, a prominent anti-apoptotic gene whose entire coding sequence is a 5' UPT generated from 5' UTR-embedded APA sites. Thus, APA is not only accountable for terminating transcripts, but also for generating downstream uncapped RNAs with translation potential and biological impact.

The small intestine is a rapidly proliferating organ that is maintained by a small population of Lgr5 -expressing intestinal stem cells (ISCs). However, several Lgr5 -negative ISC populations have been identified, and this remarkable plasticity allows the intestine to rapidly respond to both the local environment and to damage. However, the mediators of such plasticity are still largely unknown. Using intestinal organoids and mouse models, we show that upon ribosome impairment (driven by Rptor deletion, amino acid starvation, or low dose cyclohexamide treatment) ISCs gain an Lgr5 -negative, fetal-like identity. This is accompanied by a rewiring of metabolism. Our findings suggest that the ribosome can act as a sensor of nutrient availability, allowing ISCs to respond to the local nutrient environment. Mechanistically, we show that this phenotype requires the activation of ZAKɑ, which in turn activates YAP, via SRC. Together, our data reveals a central role for ribosome dynamics in intestinal stem cells, and identify the activation of ZAKɑ as a critical mediator of stem cell identity.

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis , as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5′ and 3′ UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR–Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

The glucocorticoid receptor (GR) regulates gene expression, governing aspects of homeostasis, but is also involved in cancer. Pharmacological GR activation is frequently used to alleviate therapy-related side-effects. While prior studies have shown GR activation might also have anti-proliferative action on tumours, the underpinnings of glucocorticoid action and its direct effectors in non-lymphoid solid cancers remain elusive. Here, we study the mechanisms of glucocorticoid response, focusing on lung cancer. We show that GR activation induces reversible cancer cell dormancy characterised by anticancer drug tolerance, and activation of growth factor survival signalling accompanied by vulnerability to inhibitors. GR-induced dormancy is dependent on a single GR-target gene, CDKN1C, regulated through chromatin looping of a GR-occupied upstream distal enhancer in a SWI/SNF-dependent fashion. These insights illustrate the importance of GR signalling in non-lymphoid solid cancer biology, particularly in lung cancer, and warrant caution for use of glucocorticoids in treatment of anticancer therapy related side-effects.

Abstract Motivation Ribosome Profiling (Ribo-seq) has revolutionized the study of RNA translation by providing information on ribosome positions across all translated RNAs with nucleotide-resolution. Yet several technical limitations restrict the sequencing depth of such experiments, the most common of which is the overabundance of rRNA fragments. Various strategies can be employed to tackle this issue, including the use of commercial rRNA depletion kits. However, as they are designed for more standardized RNAseq experiments, they may perform suboptimally in Ribo-seq. In order to overcome this, it is possible to use custom biotinylated oligos complementary to the most abundant rRNA fragments, however currently no computational framework exists to aid the design of optimal oligos. Results Here, we first show that a major confounding issue is that the rRNA fragments generated via Ribo-seq vary significantly with differing experimental conditions, suggesting that a ‘one-size-fits-all’ approach may be inefficient. Therefore we developed Ribo-ODDR, an oligo design pipeline integrated with a user-friendly interface that assists in oligo selection for efficient experiment-specific rRNA depletion. Ribo-ODDR uses preliminary data to identify the most abundant rRNA fragments, and calculates the rRNA depletion efficiency of potential oligos. We experimentally show that Ribo-ODDR designed oligos outperform commercially available kits and lead to a significant increase in rRNA depletion in Ribo-seq. Availability and implementation Ribo-ODDR is freely accessible at https://github.com/fallerlab/Ribo-ODDR. Supplementary information Supplementary data are available at Bioinformatics online.

Genome editing has recently made a revolutionary development with the introduction of the CRISPR-Cas technology. The programmable CRISPR-associated Cas9 and Cas12a nucleases generate specific dsDNA breaks in the genome, after which host DNA-repair mechanisms can be manipulated to implement the desired editing. Despite this spectacular progress, the efficiency of Cas9/Cas12a-based engineering can still be improved. Here, we address the variation in guide-dependent efficiency of Cas12a, and set out to reveal the molecular basis of this phenomenon. We established a sensitive and robust in vivo targeting assay based on loss of a target plasmid encoding the red fluorescent protein (mRFP). Our results suggest that folding of both the precursor guide (pre-crRNA) and the mature guide (crRNA) have a major influence on Cas12a activity. Especially, base pairing of the direct repeat, other than with itself, was found to be detrimental to the activity of Cas12a. Furthermore, we describe different approaches to minimize base-pairing interactions between the direct repeat and the variable part of the guide. We show that design of the 3' end of the guide, which is not involved in target strand base pairing, may result in substantial improvement of the guide's targeting potential and hence of its genome editing efficiency.

We consider the constrained graph alignment problem which has applications in biological network analysis. Given two input graphs $G_1=(V_1,E_1), G_2=(V_2,E_2)$, a pair of vertex mappings induces an {\it edge conservation} if the vertex pairs are adjacent in their respective graphs. %In general terms The goal is to provide a one-to-one mapping between the vertices of the input graphs in order to maximize edge conservation. However the allowed mappings are restricted since each vertex from $V_1$ (resp. $V_2$) is allowed to be mapped to at most $m_1$ (resp. $m_2$) specified vertices in $V_2$ (resp. $V_1$). Most of results in this paper deal with the case $m_2=1$ which attracted most attention in the related literature. We formulate the problem as a maximum independent set problem in a related {\em conflict graph} and investigate structural properties of this graph in terms of forbidden subgraphs. We are interested, in particular, in excluding certain wheals, fans, cliques or claws (all terms are defined in the paper), which corresponds in excluding certain cycles, paths, cliques or independent sets in the neighborhood of each vertex. Then, we investigate algorithmic consequences of some of these properties, which illustrates the potential of this approach and raises new horizons for further works. In particular this approach allows us to reinterpret a known polynomial case in terms of conflict graph and to improve known approximation and fixed-parameter tractability results through efficiently solving the maximum independent set problem in conflict graphs. Some of our new approximation results involve approximation ratios that are function of the optimal value, in particular its square root; this kind of results cannot be achieved for maximum independent set in general graphs.

Ferhat Alkan , Anne Wenzel , Oana Palasca , Peter Kerpedjiev , Anders Frost Rudebeck , and Jan Gorodkin 1,2,∗ 1Center for non-coding RNA in Technology and Health, 2Department of Veterinary Clinical and Animal Science, University of Copenhagen, Grønneg̊ardsvej 3, 1870 Frederiksberg C, Denmark, 3Novo Nordisk Foundation Center for Protein Research, University of Copenhagen, Blegdamsvej 3B, 2200 Copenhagen N, Denmark, 4Institute for Theoretical Chemistry, University of Vienna, Währingerstraße 17, 1090 Wien, Austria September 7, 2018

BackgroundRecent experimental efforts of CRISPR-Cas9 systems have shown that off-target binding and cleavage are a concern for the system and that this is highly dependent on the selected guide RNA (gRNA) design. Computational predictions of off-targets have been proposed as an attractive and more feasible alternative to tedious experimental efforts. However, accurate scoring of the high number of putative off-targets plays a key role for the success of computational off-targeting assessment.ResultsWe present an approximate binding energy model for the Cas9-gRNA-DNA complex, which systematically combines the energy parameters obtained for RNA-RNA, DNA-DNA, and RNA-DNA duplexes. Based on this model, two novel off-target assessment methods for gRNA selection in CRISPR-Cas9 applications are introduced: CRISPRoff to assign confidence scores to predicted off-targets and CRISPRspec to measure the specificity of the gRNA. We benchmark the methods against current state-of-the-art methods and show that both are in better agreement with experimental results. Furthermore, we show significant evidence supporting the inverse relationship between the on-target cleavage efficiency and specificity of the system, in which introduced binding energies are key components.ConclusionsThe impact of the binding energies provides a direction for further studies of off-targeting mechanisms. The performance of CRISPRoff and CRISPRspec enables more accurate off-target evaluation for gRNA selections, prior to any CRISPR-Cas9 genome-editing application. For given gRNA sequences or all potential gRNAs in a given target region, CRISPRoff-based off-target predictions and CRISPRspec-based specificity evaluations can be carried out through our webserver at https://rth.dk/resources/crispr/.