Context-specific regulation of cell survival by a miRNA-controlled BIM rheostat.

Knockout of the ubiquitously expressed miRNA-17∼92 cluster in mice produces a lethal developmental lung defect, skeletal abnormalities, and blocked B lymphopoiesis. A shared target of miR-17∼92 miRNAs is the pro-apoptotic protein BIM, central to life-death decisions in mammalian cells. To clarify the contribution of miR-17∼92:Bim interactions to the complex miR-17∼92 knockout phenotype, we used a system of conditional mutagenesis of the nine Bim 3' UTR miR-17∼92 seed matches. Blocking miR-17∼92:Bim interactions early in development phenocopied the lethal lung phenotype of miR-17∼92 ablation and generated a skeletal kinky tail. In the hematopoietic system, instead of causing the predicted B cell developmental block, it produced a selective inability of B cells to resist cellular stress; and prevented B and T cell hyperplasia caused by Bim haploinsufficiency. Thus, the interaction of miR-17∼92 with a single target is essential for life, and BIM regulation by miRNAs serves as a rheostat controlling cell survival in specific physiological contexts.

Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

The Lupus Autoantigen La Prevents Mis-channeling of tRNA Fragments into the Human MicroRNA Pathway.

The Lupus autoantigen La is an RNA-binding protein that stabilizes RNA polymerase III (Pol III) transcripts and supports RNA folding and has in addition been implicated in the mammalian microRNA (miRNA) pathway. Here, we have analyzed effects of La depletion on Argonaute (Ago)-bound small RNAs in human cells. We find that in the absence of La, distinct tRNA fragments are loaded into Ago proteins. Thus, La functions as gatekeeper ensuring correct tRNA maturation and protecting the miRNA pathway from potentially functional tRNA fragments. However, one specific isoleucin pre-tRNA produces both a functional tRNA and a miRNA even when La is present. We demonstrate that the fully complementary 5' leader and 3' trailer of the pre-tRNA-Ile form a double-stranded RNA molecule that has low affinity to La. Instead, Exportin-5 (Xpo5) recognizes it as miRNA precursor and transports it into the cytoplasm for Dicer processing and Ago loading.

Unambiguous identification of miRNA:target site interactions by different types of ligation reactions.

To exert regulatory function, miRNAs guide Argonaute (AGO) proteins to partially complementary sites on target RNAs. Crosslinking and immunoprecipitation (CLIP) assays are state-of-the-art to map AGO binding sites, but assigning the targeting miRNA to these sites relies on bioinformatics predictions and is therefore indirect. To directly and unambiguously identify miRNA:target site interactions, we modified our CLIP methodology in C. elegans to experimentally ligate miRNAs to their target sites. Unexpectedly, ligation reactions also occurred in the absence of the exogenous ligase. Our in vivo data set and reanalysis of published mammalian AGO-CLIP data for miRNA-chimeras yielded ∼17,000 miRNA:target site interactions. Analysis of interactions and extensive experimental validation of chimera-discovered targets of viral miRNAs suggest that our strategy identifies canonical, noncanonical, and nonconserved miRNA:targets. About 80% of miRNA interactions have perfect or partial seed complementarity. In summary, analysis of miRNA:target chimeras enables the systematic, context-specific, in vivo discovery of miRNA binding.

How epigenetic mutations can affect genetic evolution: model and mechanism.

We hypothesize that heritable epigenetic changes can affect rates of fitness increase as well as patterns of genotypic and phenotypic change during adaptation. In particular, we suggest that when natural selection acts on pure epigenetic variation in addition to genetic variation, populations adapt faster, and adaptive phenotypes can arise before any genetic changes. This may make it difficult to reconcile the timing of adaptive events detected using conventional population genetics tools based on DNA sequence data with environmental drivers of adaptation, such as changes in climate. Epigenetic modifications are frequently associated with somatic cell differentiation, but recently epigenetic changes have been found that can be transmitted over many generations. Here, we show how the interplay of these heritable epigenetic changes with genetic changes can affect adaptive evolution, and how epigenetic changes affect the signature of selection in the genetic record.

Quantitative analysis of competition in posttranscriptional regulation reveals a novel signature in target expression variation.

When small RNAs are loaded onto Argonaute proteins they can form the RNA-induced silencing complexes (RISCs), which mediate RNA interference (RNAi). RISC-formation is dependent on a shared pool of Argonaute proteins and RISC-loading factors, and is susceptible to competition among small RNAs. We present a mathematical model that aims to understand how small RNA competition for RISC-formation affects target gene repression. We discuss that small RNA activity is limited by RISC-formation, RISC-degradation, and the availability of Argonautes. We show that different competition conditions for RISC-loading result in different signatures of RNAi determined also by the amount of RISC-recycling taking place. In particular, we find that the small RNAs, although less efficient at RISC-formation, can perform in the low RISC-recycling range as well as their more effective counterparts. Additionally, we predict that under conditions of low RISC-loading efficiency and high RISC-recycling, the variation in target levels increases linearly with the target transcription rate. Furthermore, we show that RISC-recycling determines the effect that Argonaute scarcity conditions have on target expression variation. Our observations, taken together, offer a framework of predictions that can be used to infer from data the particular characteristics of underlying RNAi activity.

Can we always sweep the details of RNA-processing under the carpet?

RNA molecules follow a succession of enzyme-mediated processing steps from transcription to maturation. The participating enzymes, for example the spliceosome for mRNAs and Drosha and Dicer for microRNAs, are also produced in the cell and their copy-numbers fluctuate over time. Enzyme copy-number changes affect the processing rate of the substrate molecules; high enzyme numbers increase the processing rate, while low enzyme numbers decrease it. We study different RNA-processing cascades where enzyme copy-numbers are either fixed or fluctuate. We find that for the fixed enzyme copy-numbers, the substrates at steady-state are Poisson-distributed, and the whole RNA cascade dynamics can be understood as a single birth-death process of the mature RNA product. In this case, solely fluctuations in the timing of RNA processing lead to variation in the number of RNA molecules. However, we show analytically and numerically that when enzyme copy-numbers fluctuate, the strength of RNA fluctuations increases linearly with the RNA transcription rate. This linear effect becomes stronger as the speed of enzyme dynamics decreases relative to the speed of RNA dynamics. Interestingly, we find that under certain conditions, the RNA cascade can reduce the strength of fluctuations in the expression level of the mature RNA product. Finally, by investigating the effects of processing polymorphisms, we show that it is possible for the effects of transcriptional polymorphisms to be enhanced, reduced or even reversed. Our results provide a framework to understand the dynamics of RNA processing.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

A Reproducible Bioprinted 3D Tumor Model Serves as a Preselection Tool for CAR T Cell Therapy Optimization.

Chimeric antigen receptor (CAR) T cell performance against solid tumors in mouse models and clinical trials is often less effective than predicted by CAR construct selection in two-dimensional (2D) cocultures. Three-dimensional (3D) solid tumor architecture is likely to be crucial for CAR T cell efficacy. We used a three-dimensional (3D) bioprinting approach for large-scale generation of highly reproducible 3D human tumor models for the test case, neuroblastoma, and compared these to 2D cocultures for evaluation of CAR T cells targeting the L1 cell adhesion molecule, L1CAM. CAR T cells infiltrated the model, and both CAR T and tumor cells were viable for long-term experiments and could be isolated as single-cell suspensions for whole-cell assays quantifying CAR T cell activation, effector function and tumor cell cytotoxicity. L1CAM-specific CAR T cell activation by neuroblastoma cells was stronger in the 3D model than in 2D cocultures, but neuroblastoma cell lysis was lower. The bioprinted 3D neuroblastoma model is highly reproducible and allows detection and quantification of CAR T cell tumor infiltration, representing a superior in vitro analysis tool for preclinical CAR T cell characterization likely to better select CAR T cells for in vivo performance than 2D cocultures.

Transcriptomic profiling of SARS-CoV-2 infected human cell lines identifies HSP90 as target for COVID-19 therapy.

Detailed knowledge of the molecular biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is crucial for understanding of viral replication, host responses, and disease progression. Here, we report gene expression profiles of three SARS-CoV- and SARS-CoV-2-infected human cell lines. SARS-CoV-2 elicited an approximately two-fold higher stimulation of the innate immune response compared to SARS-CoV in the human epithelial cell line Calu-3, including induction of miRNA-155. Single-cell RNA sequencing of infected cells showed that genes induced by virus infections were broadly upregulated, whereas interferon beta/lambda genes, a pro-inflammatory cytokines such as IL-6, were expressed only in small subsets of infected cells. Temporal analysis suggested that transcriptional activities of interferon regulatory factors precede those of nuclear factor κB. Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 activity resulted in a reduction of viral replication and pro-inflammatory cytokine expression in primary human airway epithelial cells.

Small-Molecule Dual PLK1 and BRD4 Inhibitors are Active Against Preclinical Models of Pediatric Solid Tumors.

Simultaneous inhibition of multiple molecular targets is an established strategy to improve the continuance of clinical response to therapy. Here, we screened 49 molecules with dual nanomolar inhibitory activity against BRD4 and PLK1, best classified as dual kinase-bromodomain inhibitors, in pediatric tumor cell lines for their antitumor activity. We identified two candidate dual kinase-bromodomain inhibitors with strong and tumor-specific activity against neuroblastoma, medulloblastoma, and rhabdomyosarcoma tumor cells. Dual PLK1 and BRD4 inhibitor treatment suppressed proliferation and induced apoptosis in pediatric tumor cell lines at low nanomolar concentrations. This was associated with reduced MYCN-driven gene expression as assessed by RNA sequencing. Treatment of patient-derived xenografts with dual inhibitor UMB103 led to significant tumor regression. We demonstrate that concurrent inhibition of two central regulators of MYC protein family of protooncogenes, BRD4, and PLK1, with single small molecules has strong and specific antitumor effects in preclinical pediatric cancer models.

Publisher Correction: Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

Extrachromosomal circularization of DNA is an important genomic feature in cancer. However, the structure, composition and genome-wide frequency of extrachromosomal circular DNA have not yet been profiled extensively. Here, we combine genomic and transcriptomic approaches to describe the landscape of extrachromosomal circular DNA in neuroblastoma, a tumor arising in childhood from primitive cells of the sympathetic nervous system. Our analysis identifies and characterizes a wide catalog of somatically acquired and undescribed extrachromosomal circular DNAs. Moreover, we find that extrachromosomal circular DNAs are an unanticipated major source of somatic rearrangements, contributing to oncogenic remodeling through chimeric circularization and reintegration of circular DNA into the linear genome. Cancer-causing lesions can emerge out of circle-derived rearrangements and are associated with adverse clinical outcome. It is highly probable that circle-derived rearrangements represent an ongoing mutagenic process. Thus, extrachromosomal circular DNAs represent a multihit mutagenic process, with important functional and clinical implications for the origins of genomic remodeling in cancer.

Prohibitin promotes de-differentiation and is a potential therapeutic target in neuroblastoma.

Gain of the long arm of chromosome 17 (17q) is a cytogenetic hallmark of high-risk neuroblastoma, yet its contribution to neuroblastoma pathogenesis remains incompletely understood. Combining whole-genome and RNA sequencing of neuroblastomas, we identified the prohibitin (PHB) gene as highly expressed in tumors with 17q gain. High PHB expression correlated with poor prognosis and was associated with loss of gene expression programs promoting neuronal development and differentiation. PHB depletion induced differentiation and apoptosis and slowed cell cycle progression of neuroblastoma cells, at least in part through impaired ERK1/2 activation. Conversely, ectopic expression of PHB was sufficient to increase proliferation of neuroblastoma cells and was associated with suppression of markers associated with neuronal differentiation and favorable neuroblastoma outcome. Thus, PHB is a 17q oncogene in neuroblastoma that promotes tumor cell proliferation, and de-differentiation.

Spatiotemporal m(i)RNA Architecture and 3' UTR Regulation in the C. elegans Germline.

In animal germlines, regulation of cell proliferation and differentiation is particularly important but poorly understood. Here, using a cryo-cut approach, we mapped RNA expression along the Caenorhabditis elegans germline and, using mutants, dissected gene regulatory mechanisms that control spatiotemporal expression. We detected, at near single-cell resolution, >10,000 mRNAs, >300 miRNAs, and numerous unannotated miRNAs. Most RNAs were organized in distinct spatial patterns. Germline-specific miRNAs and their targets were co-localized. Moreover, we observed differential 3' UTR isoform usage for hundreds of mRNAs. In tumorous gld-2 gld-1 mutants, gene expression was strongly perturbed. In particular, differential 3' UTR usage was significantly impaired. We propose that PIE-1, a transcriptional repressor, functions to maintain spatial gene expression. Our data also suggest that cpsf-4 and fipp-1 control differential 3' UTR usage for hundreds of genes. Finally, we constructed a "virtual gonad" enabling "virtual in situ hybridizations" and access to all data (https://shiny.mdc-berlin.de/spacegerm/).

Characterization of Transcription Termination-Associated RNAs: New Insights into their Biogenesis, Tailing, and Expression in Primary Tumors.

Next-generation sequencing has uncovered novel classes of small RNAs (sRNAs) in eukaryotes, in addition to the well-known miRNAs, siRNAs, and piRNAs. In particular, sRNA species arise from transcription start sites (TSSs) and the transcription termination sites (TTSs) of genes. However, a detailed characterization of these new classes of sRNAs is still lacking. Here, we present a comprehensive study of sRNAs derived from TTSs of expressed genes (TTSa-RNAs) in human cell lines and primary tissues. Taking advantage of sRNA-sequencing, we show that TTSa-RNAs are present in the nuclei of human cells, are loaded onto both AGO1 and AGO2, and their biogenesis does not require DICER and AGO2 endonucleolytic activity. TTSa-RNAs display a strong bias against a G residue in the first position at 5' end, a known feature of AGO-bound sRNAs, and a peculiar oligoA tail at 3' end. AGO-bound TTSa-RNAs derive from genes involved in cell cycle progression regulation and DNA integrity checkpoints. Finally, we provide evidence that TTSa-RNAs can be detected by sRNA-Seq in primary human tissue, and their expression increases in tumor samples as compared to nontumor tissues, suggesting that in the future, TTSa-RNAs might be explored as biomarker for diagnosis or prognosis of human malignancies.

Loss of a mammalian circular RNA locus causes miRNA deregulation and affects brain function.

Hundreds of circular RNAs (circRNAs) are highly abundant in the mammalian brain, often with conserved expression. Here we show that the circRNA Cdr1as is massively bound by the microRNAs (miRNAs) miR-7 and miR-671 in human and mouse brains. When the Cdr1as locus was removed from the mouse genome, knockout animals displayed impaired sensorimotor gating-a deficit in the ability to filter out unnecessary information-which is associated with neuropsychiatric disorders. Electrophysiological recordings revealed dysfunctional synaptic transmission. Expression of miR-7 and miR-671 was specifically and posttranscriptionally misregulated in all brain regions analyzed. Expression of immediate early genes such as Fos, a direct miR-7 target, was enhanced in Cdr1as-deficient brains, providing a possible molecular link to the behavioral phenotype. Our data indicate an in vivo loss-of-function circRNA phenotype and suggest that interactions between Cdr1as and miRNAs are important for normal brain function.

Epigenetic dynamics of monocyte-to-macrophage differentiation.

BACKGROUND: Monocyte-to-macrophage differentiation involves major biochemical and structural changes. In order to elucidate the role of gene regulatory changes during this process, we used high-throughput sequencing to analyze the complete transcriptome and epigenome of human monocytes that were differentiated in vitro by addition of colony-stimulating factor 1 in serum-free medium. RESULTS: Numerous mRNAs and miRNAs were significantly up- or down-regulated. More than 100 discrete DNA regions, most often far away from transcription start sites, were rapidly demethylated by the ten eleven translocation enzymes, became nucleosome-free and gained histone marks indicative of active enhancers. These regions were unique for macrophages and associated with genes involved in the regulation of the actin cytoskeleton, phagocytosis and innate immune response. CONCLUSIONS: In summary, we have discovered a phagocytic gene network that is repressed by DNA methylation in monocytes and rapidly de-repressed after the onset of macrophage differentiation.