microRNAs (miRNAs) are a large class of small non-coding RNAs which post-transcriptionally regulate the expression of a large fraction of all animal genes and are important in a wide range of biological processes. Recent advances in high-throughput sequencing allow miRNA detection at unprecedented sensitivity, but the computational task of accurately identifying the miRNAs in the background of sequenced RNAs remains challenging. For this purpose, we have designed miRDeep2, a substantially improved algorithm which identifies canonical and non-canonical miRNAs such as those derived from transposable elements and informs on high-confidence candidates that are detected in multiple independent samples. Analyzing data from seven animal species representing the major animal clades, miRDeep2 identified miRNAs with an accuracy of 98.6-99.9% and reported hundreds of novel miRNAs. To test the accuracy of miRDeep2, we knocked down the miRNA biogenesis pathway in a human cell line and sequenced small RNAs before and after. The vast majority of the >100 novel miRNAs expressed in this cell line were indeed specifically downregulated, validating most miRDeep2 predictions. Last, a new miRNA expression profiling routine, low time and memory usage and user-friendly interactive graphic output can make miRDeep2 useful to a wide range of researchers.

Repeatedly activated T helper 1 (Th1) cells present during chronic inflammation can efficiently adapt to the inflammatory milieu, for example, by expressing the transcription factor Twist1, which limits the immunopathology caused by Th1 cells. Here, we show that in repeatedly activated murine Th1 cells, Twist1 and T-bet induce expression of microRNA-148a (miR-148a). miR-148a regulates expression of the proapoptotic gene Bim, resulting in a decreased Bim/Bcl2 ratio. Inhibition of miR-148a by antagomirs in repeatedly activated Th1 cells increases the expression of Bim, leading to enhanced apoptosis. Knockdown of Bim expression by siRNA in miR-148a antagomir-treated cells restores viability of the Th1 cells, demonstrating that miR-148a controls survival by regulating Bim expression. Thus, Twist1 and T-bet not only control the differentiation and function of Th1 cells, but also their persistence in chronic inflammation.

Post-transcriptional gene regulation (PTGR) of mRNA turnover, localization and translation is mediated by microRNAs (miRNAs) and RNA-binding proteins (RBPs). These regulators exert their effects by binding to specific sequences within their target mRNAs. Increasing evidence suggests that competition for binding is a fundamental principle of PTGR. Not only can miRNAs be sequestered and neutralized by the targets with which they interact through a process termed 'sponging', but competition between binding sites on different RNAs may also lead to regulatory crosstalk between transcripts. Here, we quantitatively model competition effects under physiological conditions and review the role of endogenous sponges for PTGR in light of the key features that emerge.

The development and homeostasis of multicellular organisms relies on gene regulation within individual constituent cells. Gene regulatory circuits that increase the robustness of gene expression frequently incorporate microRNAs as post-transcriptional regulators. Computational approaches, synthetic gene circuits and observations in model organisms predict that the co-regulation of microRNAs and their target mRNAs can reduce cell-to-cell variability in the expression of target genes. However, whether microRNAs directly regulate variability of endogenous gene expression remains to be tested in mammalian cells. Here we use quantitative flow cytometry to show that microRNAs impact on cell-to-cell variability of protein expression in developing mouse thymocytes. We find two distinct mechanisms that control variation in the activation-induced expression of the microRNA target CD69. First, the expression of miR-17 and miR-20a, two members of the miR-17-92 cluster, is co-regulated with the target mRNA Cd69 to form an activation-induced incoherent feed-forward loop. Another microRNA, miR-181a, acts at least in part upstream of the target mRNA Cd69 to modulate cellular responses to activation. The ability of microRNAs to render gene expression more uniform across mammalian cell populations may be important for normal development and for disease.

In multicellular organisms, cell type diversity and fate depend on specific sets of transcript isoforms generated by post-transcriptional RNA processing. Here, we used Schmidtea mediterranea, a flatworm with extraordinary regenerative abilities and a large pool of adult stem cells, as an in vivo model to study the role of Uridyl-rich small nuclear RNAs (UsnRNAs), which participate in multiple RNA processing reactions including splicing, in stem cell regulation. We characterized the planarian UsnRNA repertoire, identified stem cell-enriched variants and obtained strong evidence for an increased rate of UsnRNA 3'-processing in stem cells compared to their differentiated counterparts. Consistently, components of the Integrator complex showed stem cell-enriched expression and their depletion by RNAi disrupted UsnRNA processing resulting in global changes of splicing patterns and reduced processing of histone mRNAs. Interestingly, loss of Integrator complex function disrupted both stem cell maintenance and regeneration of tissues. Our data show that the function of the Integrator complex in UsnRNA 3'-processing is conserved in planarians and essential for maintaining their stem cell pool. We propose that cell type-specific modulation of UsnRNA composition and maturation contributes to in vivo cell fate choices, such as stem cell self-renewal in planarians.

The oocyte-to-embryo transition (OET) is thought to be mainly driven by post-transcriptional gene regulation. However, expression of both RNAs and proteins during the OET has not been comprehensively assayed. Furthermore, specific molecular mechanisms that regulate gene expression during OET are largely unknown. Here, we quantify and analyze transcriptome-wide, expression of mRNAs and thousands of proteins in Caenorhabditis elegans oocytes, 1-cell, and 2-cell embryos. This represents a first comprehensive gene expression atlas during the OET in animals. We discovered a first wave of degradation in which thousands of mRNAs are cleared shortly after fertilization. Sequence analysis revealed a statistically highly significant presence of a polyC motif in the 3' untranslated regions of most of these degraded mRNAs. Transgenic reporter assays demonstrated that this polyC motif is required and sufficient for mRNA degradation after fertilization. We show that orthologs of human polyC-binding protein specifically bind this motif. Our data suggest a mechanism in which the polyC motif and binding partners direct degradation of maternal mRNAs. Our data also indicate that endogenous siRNAs but not miRNAs promote mRNA clearance during the OET.

Genome-wide distribution of the majority of H2A and H3 variants (H2A, H2AX, H2AZ, macroH2A, H3.1, H3.2 and H3.3) was simultaneously investigated in mouse embryonic stem cells by chromatin immunoprecipitation sequencing. Around the transcription start site, histone variant distribution differed between genes possessing promoters of high and low CpG density, regardless of their expression levels. In the intergenic regions, regulatory elements were enriched in H2A.Z and H3.3, whereas repeat elements were abundant in H2A and macroH2A, and H3.1, respectively. Analysis of H2A and H3 variant combinations composing nucleosomes revealed that the H2A.Z and H3.3 combinations were present at a higher frequency throughout the genome than the other combinations, suggesting that H2A.Z and H3.3 associate preferentially with each other to comprise the nucleosomes independently of genome region. Finally, we found that chromatin was unstable only in regions where it was enriched in both H2A.Z and H3.3, but strongly quantified stable in regions in which only H3.3 was abundant. Therefore, histone variant composition is an important determinant of chromatin structure, which is associated with specific genomic functions.

Herpesviruses can infect a wide range of animal species. Herpes simplex virus 1 (HSV-1) is one of the eight herpesviruses that can infect humans and is prevalent worldwide. Herpesviruses have evolved multiple ways to adapt the infected cells to their needs, but knowledge about these transcriptional and post-transcriptional modifications is sparse.

microRNAs (miRNAs) are small noncoding RNAs that mediate post-transcriptional gene regulation and have emerged as essential regulators of many developmental events. The transcriptional network during early embryogenesis of the purple sea urchin, Strongylocentrotus purpuratus, is well described and can serve as an excellent model to test functional contributions of miRNAs in embryogenesis. We examined the loss of function phenotypes of major components of the miRNA biogenesis pathway. Inhibition of de novo synthesis of Drosha and Dicer in the embryo led to consistent developmental defects, a failure to gastrulate, and embryonic lethality, including changes in the steady state levels of transcription factors and signaling molecules involved in germ layer specification. We annotated and profiled small RNA expression from the ovary and several early embryonic stages by deep sequencing followed by computational analysis. miRNAs as well as a large population of putative piRNAs (piwi-interacting RNAs) had dynamic accumulation profiles through early development. Defects in morphogenesis caused by loss of Drosha could be rescued with four miRNAs. Taken together our results indicate that post-transcriptional gene regulation directed by miRNAs is functionally important for early embryogenesis and is an integral part of the early embryonic gene regulatory network in S. purpuratus.

Metazoan miRNAs regulate protein-coding genes by binding the 3' UTR of cognate mRNAs. Identifying targets for the 115 known C. elegans miRNAs is essential for understanding their function.

We previously showed that type I interferon-induced, Cre-mediated ablation of surface BCR expression in mature B cells through Ig-heavy chain deletion results in apoptosis of these cells. This led to the hypothesis that survival signals from the BCR are vital for mature B cells. Here, we test two critical assumptions of this model. First, we demonstrate loss of mature B cells upon induced mutation of a signaling module of the BCR, not precluding BCR surface expression. Second, we show that the cells are also lost upon BCR inactivation in the absence of an exogenous inducer like interferon, excluding that cell death depends on previous cellular activation by the latter. Kinetic data demonstrate that BCR-less mature B cells have a severely reduced lifespan, with a half-life of 3-6 days. Together these results establish that BCR signaling is required to keep resting mature B cells alive in vivo.

To explore the role of Dicer-dependent control mechanisms in B lymphocyte development, we ablated this enzyme in early B cell progenitors. This resulted in a developmental block at the pro- to pre-B cell transition. Gene-expression profiling revealed a miR-17 approximately 92 signature in the 3'UTRs of genes upregulated in Dicer-deficient pro-B cells; a top miR-17 approximately 92 target, the proapoptotic molecule Bim, was highly upregulated. Accordingly, B cell development could be partially rescued by ablation of Bim or transgenic expression of the prosurvival protein Bcl-2. This allowed us to assess the impact of Dicer deficiency on the V(D)J recombination program in developing B cells. We found intact Ig gene rearrangements in immunoglobulin heavy (IgH) and kappa chain loci, but increased sterile transcription and usage of D(H) elements of the DSP family in IgH, and increased N sequence addition in Igkappa due to deregulated transcription of the terminal deoxynucleotidyl transferase gene.

In bioinformatics, as well as other computationally intensive research fields, there is a need for workflows that can reliably produce consistent output, from known sources, independent of the software environment or configuration settings of the machine on which they are executed. Indeed, this is essential for controlled comparison between different observations and for the wider dissemination of workflows. However, providing this type of reproducibility and traceability is often complicated by the need to accommodate the myriad dependencies included in a larger body of software, each of which generally comes in various versions. Moreover, in many fields (bioinformatics being a prime example), these versions are subject to continual change due to rapidly evolving technologies, further complicating problems related to reproducibility. Here, we propose a principled approach for building analysis pipelines and managing their dependencies with GNU Guix. As a case study to demonstrate the utility of our approach, we present a set of highly reproducible pipelines called PiGx for the analysis of RNA sequencing, chromatin immunoprecipitation sequencing, bisulfite-treated DNA sequencing, and single-cell resolution RNA sequencing. All pipelines process raw experimental data and generate reports containing publication-ready plots and figures, with interactive report elements and standard observables. Users may install these highly reproducible packages and apply them to their own datasets without any special computational expertise beyond the use of the command line. We hope such a toolkit will provide immediate benefit to laboratory workers wishing to process their own datasets or bioinformaticians seeking to automate all, or parts of, their analyses. In the long term, we hope our approach to reproducibility will serve as a blueprint for reproducible workflows in other areas. Our pipelines, along with their corresponding documentation and sample reports, are available at http://bioinformatics.mdc-berlin.de/pigx.

Herpesvirus infection initiates a range of perturbations in the host cell, which remain poorly understood at the level of individual cells. Here, we quantify the transcriptome of single human primary fibroblasts during the first hours of lytic infection with HSV-1. By applying a generalizable analysis scheme, we define a precise temporal order of early viral gene expression and propose a set-wise emergence of viral genes. We identify host cell genes and pathways relevant for infection by combining three different computational approaches: gene and pathway overdispersion analysis, prediction of cell-state transition probabilities, as well as future cell states. One transcriptional program, which correlates with increased resistance to infection, implicates the transcription factor NRF2. Consequently, Bardoxolone methyl and Sulforaphane, two known NRF2 agonists, impair virus production, suggesting that NRF2 activation restricts viral infection. Our study provides insights into early stages of HSV-1 infection and serves as a general blueprint for the investigation of heterogeneous cell states in virus infection.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

During their lifetime, animals must adapt their behavior to survive in changing environments. This ability requires the nervous system to undergo adjustments at distinct temporal scales, from short-term dynamic changes in expression of neurotransmitters and receptors to longer-term growth, spatial and connectivity reorganization, while integrating external stimuli. The nematode Caenorhabditis elegans provides a model of nervous system plasticity, in particular its dauer exit decision. Under unfavorable conditions, larvae will enter the non-feeding and non-reproductive stress-resistant dauer stage and adapt their behavior to cope with the harsh new environment, with active reversal under improved conditions leading to resumption of reproductive development. However, how different environmental stimuli regulate the exit decision mechanism and thereby drive the larva's behavioral change is unknown. To fill this gap and provide insights on behavioral changes over extended periods of time, we developed a new open hardware method for long-term imaging (12h) of C. elegans larvae.

Understanding the molecular signatures of colorectal cancer progression under chemotherapeutic treatment will be crucial for the success of future therapy improvements. Here, we used a xenograft-based mouse model to investigate, how whole transcriptome signatures change during metastatic colorectal cancer progression and how such signatures are affected by LDM chemotherapy using RNA sequencing. We characterized mRNAs as well as non-coding RNAs such as microRNAs, long non-coding RNAs and circular RNAs in colorectal-cancer bearing mice with or without LDM chemotherapy. Furthermore, we found that circZNF609 functions as oncogene, since over-expression studies lead to an increased tumor growth while specific knock down results in smaller tumors. Our data represent novel insights into the relevance of non-coding and circRNAs in colorectal cancer and provide a comprehensive resource of gene expression changes in primary tumors and metastases. In addition, we present candidate genes that could be important modulators for successful LDM chemotherapy.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Transcription is regulated by a multitude of activators and repressors, which bind to the RNA polymerase II (Pol II) machinery and modulate its progression. Death-inducer obliterator 3 (DIDO3) and PHD finger protein 3 (PHF3) are paralogue proteins that regulate transcription elongation by docking onto phosphorylated serine-2 in the C-terminal domain (CTD) of Pol II through their SPOC domains. Here, we show that DIDO3 and PHF3 form a complex that bridges the Pol II elongation machinery with chromatin and RNA processing factors and tethers Pol II in a phase-separated microenvironment. Their SPOC domains and C-terminal intrinsically disordered regions are critical for transcription regulation. PHF3 and DIDO exert cooperative and antagonistic effects on the expression of neuronal genes and are both essential for neuronal differentiation. In the absence of PHF3, DIDO3 is upregulated as a compensatory mechanism. In addition to shared gene targets, DIDO specifically regulates genes required for lipid metabolism. Collectively, our work reveals multiple layers of gene expression regulation by the DIDO3 and PHF3 paralogues, which have specific, co-regulatory and redundant functions in transcription.

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.