Search

Merkel cell polyomavirus (MCPyV) is the causative agent of the majority of Merkel cell carcinomas (MCC). The virus has limited coding capacity, with its early viral proteins, large T (LT) and small T (sT), being multifunctional and contributing to infection and transformation. A fundamental difference in early viral gene expression between infection and MCPyV-driven tumorigenesis is the expression of a truncated LT (LTtr) in the tumor. In contrast, sT is expressed in both conditions and contributes significantly to oncogenesis. Here, we identified novel functions of early viral proteins by performing genome-wide transcriptome and chromatin studies in primary human fibroblasts. Due to current limitations in infection and tumorigenesis models, we mimic these conditions by ectopically expressing sT, LT or LTtr, individually or in combination, at different time points. In addition to its known function in cell cycle and inflammation modulation, we reveal a fundamentally new function of sT. We show that sT regulates the type I interferon (IFN) response downstream of the type I interferon receptor (IFNAR) by interfering with the interferon-stimulated gene factor 3 (ISGF3)-induced interferon-stimulated gene (ISG) response. Expression of sT leads to a reduction in the expression of interferon regulatory factor 9 (IRF9) which is a central component of the ISGF3 complex. We further show that this function of sT is conserved in BKPyV. We provide a first mechanistic understanding of which early viral proteins trigger and control the type I IFN response, which may influence MCPyV infection, persistence and, during MCC progression, regulation of the tumor microenvironment.

Abstract Herpes simplex virus 1 (HSV-1) infection and stress responses disrupt transcription termination by RNA Polymerase II (Pol II). In HSV-1 infection, but not upon salt or heat stress, this is accompanied by a dramatic increase in chromatin accessibility downstream of genes. Here, we show that the HSV-1 immediate-early protein ICP22 is both necessary and sufficient to induce downstream open chromatin regions (dOCRs) when transcription termination is disrupted by the viral ICP27 protein. This is accompanied by a marked ICP22-dependent loss of histones downstream of affected genes consistent with impaired histone repositioning in the wake of Pol II. Efficient knock-down of the ICP22-interacting histone chaperone FACT is not sufficient to induce dOCRs in ΔICP22 infection but increases dOCR induction in wild-type HSV-1 infection. Interestingly, this is accompanied by a marked increase in chromatin accessibility within gene bodies. We propose a model in which allosteric changes in Pol II composition downstream of genes and ICP22-mediated interference with FACT activity explain the differential impairment of histone repositioning downstream of genes in the wake of Pol II in HSV-1 infection.

The genomes of both human cytomegalovirus (HCMV) and murine cytomegalovirus (MCMV) were first sequenced over 20 years ago. Similar to HCMV, the MCMV genome had initially been proposed to harbor ≈170 open reading frames (ORFs). More recently, omics approaches revealed HCMV gene expression to be substantially more complex comprising several hundred viral ORFs. Here, we provide a state-of-the art reannotation of lytic MCMV gene expression based on integrative analysis of a large set of omics data. Our data reveal 365 viral transcription start sites (TiSS) that give rise to 380 and 454 viral transcripts and ORFs, respectively. The latter include >200 small ORFs, some of which represented the most highly expressed viral gene products. By combining TiSS profiling with metabolic RNA labelling and chemical nucleotide conversion sequencing (dSLAM-seq), we provide a detailed picture of the expression kinetics of viral transcription. This not only resulted in the identification of a novel MCMV immediate early transcript encoding the m166.5 ORF, which we termed ie4, but also revealed a group of well-expressed viral transcripts that are induced later than canonical true late genes and contain an initiator element (Inr) but no TATA- or TATT-box in their core promoters. We show that viral upstream ORFs (uORFs) tune gene expression of longer viral ORFs expressed in cis at translational level. Finally, we identify a truncated isoform of the viral NK-cell immune evasin m145 arising from a viral TiSS downstream of the canonical m145 mRNA. Despite being ≈5-fold more abundantly expressed than the canonical m145 protein it was not required for downregulating the NK cell ligand, MULT-I. In summary, our work will pave the way for future mechanistic studies on previously unknown cytomegalovirus gene products in an important virus animal model.

Huang et al demonstrate that glioblastoma cells upregulate axon guidance molecule Plexin-B2 to gain invasiveness and that Plexin-B2 promotes glioblastoma cell infiltration along axon fiber tracts in intracranial transplant models by modulating cellular biomechanics.

Here, using computational integration of multi-omics data, the authors provide a detailed transcriptome and translatome of herpes simplex virus 1 (HSV-1), including previously unidentified ORFs and N-terminal extensions. The study also provides a HSV-1 genome browser and should be a valuable resource for further research.

The herpes simplex virus 1 (HSV-1) virion host shut-off (vhs) protein cleaves both cellular and viral mRNAs by a translation-initiation-dependent mechanism, which should spare circular RNAs (circRNAs). Here, we show that vhs-mediated degradation of linear mRNAs leads to an enrichment of circRNAs relative to linear mRNAs during HSV-1 infection. This was also observed in influenza A virus (IAV) infection, likely due to degradation of linear host mRNAs mediated by the IAV PA-X protein and cap-snatching RNA-dependent RNA polymerase. For most circRNAs, enrichment was not due to increased circRNA synthesis but due to a general loss of linear RNAs. In contrast, biogenesis of a circRNA originating from the long isoform (NEAT1_2) of the nuclear paraspeckle assembly transcript 1 (NEAT1) was induced both in HSV-1 infection–in a vhs-independent manner–and in IAV infection. This was associated with induction of novel linear splicing of NEAT1_2 both within and downstream of the circRNA. NEAT1_2 forms a scaffold for paraspeckles, nuclear bodies located in the interchromatin space, must likely remain unspliced for paraspeckle assembly and is up-regulated in HSV-1 and IAV infection. We show that NEAT1_2 splicing and up-regulation can be induced by ectopic co-expression of the HSV-1 immediate-early proteins ICP22 and ICP27, potentially linking increased expression and splicing of NEAT1_2. To identify other conditions with NEAT1_2 splicing, we performed a large-scale screen of published RNA-seq data. This uncovered both induction of NEAT1_2 splicing and poly(A) read-through similar to HSV-1 and IAV infection in cancer cells upon inhibition or knockdown of CDK7 or the MED1 subunit of the Mediator complex phosphorylated by CDK7. In summary, our study reveals induction of novel circular and linear NEAT1_2 splicing isoforms as a common characteristic of HSV-1 and IAV infection and highlights a potential role of CDK7 in HSV-1 or IAV infection.

Herpes simplex virus-1 (HSV-1) infection disrupts transcription termination (DoTT) of host genes, but underlying mechanisms are unclear. Here, Wang et al. show that the HSV-1 immediate early protein ICP27 induces DoTT through interaction with the mRNA 3’ processing factor CPSF and disruption of the processing complex.

Abstract Background The tumor microenvironment (TME) combines features of regulatory cytokines and immune cell populations to evade the recognition by the immune system. Myeloid-derived suppressor cells (MDSC) comprise populations of immature myeloid cells in tumor-bearing hosts with a highly immunosuppressive capacity. We could previously identify RIG-I-like helicases (RLH) as targets for the immunotherapy of pancreatic cancer inducing immunogenic tumor cell death and type I interferons (IFN) as key mediators linking innate with adaptive immunity. Methods Mice with orthotopically implanted KrasG12D p53fl/R172H Ptf1a-Cre (KPC) pancreatic tumors were treated intravenously with the RLH ligand polyinosinic-polycytidylic acid (poly(I:C)), and the immune cell environment in tumor and spleen was characterized. A comprehensive analysis of the suppressive capacity as well as the whole transcriptomic profile of isolated MDSC subsets was performed. Antigen presentation capability of MDSC from mice with ovalbumin (OVA)-expressing tumors was investigated in T cell proliferation assays. The role of IFN in MDSC function was investigated in Ifnar1 −/− mice. Results MDSC were strongly induced in orthotopic KPC-derived pancreatic cancer, and frequencies of MDSC subsets correlated with tumor weight and G-CSF serum levels, whereas other immune cell populations decreased. Administration of the RLH-ligand induced a IFN-driven immune response, with increased activation of T cells and dendritic cells (DC), and a reduced suppressive capacity of both polymorphonuclear (PMN)-MDSC and monocytic (M)-MDSC fractions. Whole transcriptomic analysis confirmed an IFN-driven gene signature of MDSC, a switch from a M2/G2- towards a M1/G1-polarized phenotype, and the induction of genes involved in the antigen presentation machinery. Nevertheless, MDSC failed to present tumor antigen to T cells. Interestingly, we found MDSC with reduced suppressive function in Ifnar1-deficient hosts; however, there was a common flaw in immune cell activation, which was reflected by defective immune cell activation and tumor control. Conclusions We provide evidence that the treatment with immunostimulatory RNA reprograms the TME of pancreatic cancer by reducing the suppressive activity of MDSC, polarizing myeloid cells into a M1-like state and recruiting DC. We postulate that tumor cell-targeting combination strategies may benefit from RLH-based TME remodeling. In addition, we provide novel insights into the dual role of IFN signaling in MDSC’s suppressive function and provide evidence that host-intrinsic IFN signaling may be critical for MDSC to gain suppressive function during tumor development.

Background: Glioblastoma (GBM), a highly malignant brain tumor, invariably recurs after therapy. Quiescent GBM cells represent a potential source of tumor recurrence, but little is known about their molecular underpinnings. Methods: Patient-derived GBM cells were engineered by CRISPR/Cas9-assisted knock-in of an inducible histone2B-GFP (iH2B-GFP) reporter to track cell division history. We utilized an in vitro 3D GBM organoid approach to isolate live quiescent GBM (qGBM) cells and their proliferative counterparts (pGBM) to compare stem cell properties and therapy resistance. Gene expression programs of qGBM and pGBM cells were analyzed by RNA-Seq and NanoString platforms. Findings: H2B-GFP-retaining qGBM cells exhibited comparable self-renewal capacity but higher therapy resistance relative to pGBM. Quiescent GBM cells expressed distinct gene programs that affect cell cycle control, metabolic adaptation, and extracellular matrix (ECM) interactions. Transcriptome analysis also revealed a mesenchymal shift in qGBM cells of both proneural and mesenchymal GBM subtypes. Bioinformatic analyses and functional assays in GBM organoids established hypoxia and TGFβ signaling as potential niche factors that promote quiescence in GBM. Finally, network co-expression analysis of TCGA glioma patient data identified gene modules that are enriched for qGBM signatures and also associated with survival rate. Interpretation: Our in vitro study in 3D GBM organoids supports the presence of a quiescent cell population that displays self-renewal capacity, high therapy resistance, and mesenchymal gene signatures. It also sheds light on how GBM cells may acquire and maintain quiescence through ECM organization and interaction with niche factors such as TGFβ and hypoxia. Our findings provide a starting point for developing strategies to tackle the quiescent population of GBM. Fund: National Institutes of Health (NIH) and Deutsche Forschungsgemeinschaft (DFG). Keywords: Glioblastoma, Tumor quiescence, H2B-GFP, GBM organoid, Proneural-mesenchymal transition, Stem cell niche

Eukaryotic gene expression is extensively regulated by cellular stress and pathogen infections. We have previously shown that herpes simplex virus 1 (HSV-1) and several cellular stresses cause widespread disruption of transcription termination (DoTT) of RNA polymerase II (RNAPII) in host genes and that the viral immediate early factor ICP27 plays an important role in HSV-1-induced DoTT. Here, we show that HSV-1 infection also leads to widespread changes in alternative polyadenylation (APA) of host mRNAs. In the majority of cases, polyadenylation shifts to upstream poly(A) sites (PAS), including many intronic PAS. Mechanistically, ICP27 contributes to HSV-1-mediated APA regulation. HSV-1- and ICP27-induced activation of intronic PAS is sequence-dependent and does not involve general inhibition of U1 snRNP. HSV1-induced intronic polyadenylation is accompanied by early termination of RNAPII. HSV-1-induced mRNAs polyadenylated at intronic PAS (IPA) are exported into the cytoplasm while APA isoforms with extended 3' UTRs are sequestered in the nuclei, both preventing the expression of the full-length gene products. Finally we provide evidence that HSV-induced IPA isoforms are translated. Together with other recent studies, our results suggest that viral infection and cellular stresses induce a multi-faceted host response that includes DoTT and changes in APA profiles.

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing “Big Data” offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Abstract Background The development of high-throughput experimental technologies, such as next-generation sequencing, have led to new challenges for handling, analyzing and integrating the resulting large and diverse datasets. Bioinformatical analysis of these data commonly requires a number of mutually dependent steps applied to numerous samples for multiple conditions and replicates. To support these analyses, a number of workflow management systems (WMSs) have been developed to allow automated execution of corresponding analysis workflows. Major advantages of WMSs are the easy reproducibility of results as well as the reusability of workflows or their components. Results In this article, we present Watchdog, a WMS for the automated analysis of large-scale experimental data. Main features include straightforward processing of replicate data, support for distributed computer systems, customizable error detection and manual intervention into workflow execution. Watchdog is implemented in Java and thus platform-independent and allows easy sharing of workflows and corresponding program modules. It provides a graphical user interface (GUI) for workflow construction using pre-defined modules as well as a helper script for creating new module definitions. Execution of workflows is possible using either the GUI or a command-line interface and a web-interface is provided for monitoring the execution status and intervening in case of errors. To illustrate its potentials on a real-life example, a comprehensive workflow and modules for the analysis of RNA-seq experiments were implemented and are provided with the software in addition to simple test examples. Conclusions Watchdog is a powerful and flexible WMS for the analysis of large-scale high-throughput experiments. We believe it will greatly benefit both users with and without programming skills who want to develop and apply bioinformatical workflows with reasonable overhead. The software, example workflows and a comprehensive documentation are freely available at www.bio.ifi.lmu.de/watchdog.

Abstract The bromodomain protein Brd4 is an epigenetic reader and plays a critical role in the development and maintenance of leukemia. Brd4 binds to acetylated histone tails and activates transcription by recruiting the positive elongation factor P-TEFb. Small molecule inhibitor JQ1 competitively binds the bromodomains of Brd4 and displaces the protein from acetylated histones. However, it remains unclear whether genes targeted by JQ1 are mainly regulated by Brd4 or by other bromodomain proteins such as Brd2 and Brd3. Here, we describe anti-proliferative dominant-negative Brd4 mutants that compete with the function of distinct Brd4 domains. We used these Brd4 mutants to compare the Brd4-specific transcriptome with the transcriptome of JQ1-treated cells. We found that most JQ1-regulated genes are also regulated by dominant-negative Brd4 mutants, including the mutant that competes with the P-TEFb recruitment function of Brd4. Importantly, JQ1 and dominant-negative Brd4 mutants regulated the same set of target genes of c-Myc, a key regulator of the JQ1 response in leukemia cells. Our results suggest that Brd4 mediates most of the anti-cancer effects of JQ1 and that the major function of Brd4 in this process is the recruitment of P-TEFb. In summary, our studies define the molecular targets of JQ1 in more detail.

Summary: Activation of immune cells results in rapid functional changes, but how such fast changes are accomplished remains enigmatic. By combining time courses of 4sU-seq, RNA-seq, ribosome profiling (RP), and RNA polymerase II (RNA Pol II) ChIP-seq during T cell activation, we illustrate genome-wide temporal dynamics for ∼10,000 genes. This approach reveals not only immediate-early and posttranscriptionally regulated genes but also coupled changes in transcription and translation for >90% of genes. Recruitment, rather than release of paused RNA Pol II, primarily mediates transcriptional changes. This coincides with a genome-wide temporary slowdown in cotranscriptional splicing, even for polyadenylated mRNAs that are localized at the chromatin. Subsequent splicing optimization correlates with increasing Ser-2 phosphorylation of the RNA Pol II carboxy-terminal domain (CTD) and activation of the positive transcription elongation factor (pTEFb). Thus, rapid de novo recruitment of RNA Pol II dictates the course of events during T cell activation, particularly transcription, splicing, and consequently translation. : Davari et al. visualize global changes in RNA Pol II binding, transcription, splicing, and translation. T cells change their functional program by rapid de novo recruitment of RNA Pol II and coupled changes in transcription and translation. This coincides with fluctuations in RNA Pol II phosphorylation and a temporary reduction in cotranscriptional splicing. Keywords: RNA Pol II, cotranscriptional splicing, T cell activation, ribosome profiling, 4sU, H3K36, Ser-5 RNA Pol II, Ser-2 RNA Pol II, immune response, immediate-early genes

Lytic herpes simplex virus 1 (HSV-1) infection triggers disruption of transcription termination (DoTT) of most cellular genes, resulting in extensive intergenic transcription. Similarly, cellular stress responses lead to gene-specific transcription downstream of genes (DoG). In this study, we performed a detailed comparison of DoTT/DoG transcription between HSV-1 infection, salt and heat stress in primary human fibroblasts using 4sU-seq and ATAC-seq. Although DoTT at late times of HSV-1 infection was substantially more prominent than DoG transcription in salt and heat stress, poly(A) read-through due to DoTT/DoG transcription and affected genes were significantly correlated between all three conditions, in particular at earlier times of infection. We speculate that HSV-1 either directly usurps a cellular stress response or disrupts the transcription termination machinery in other ways but with similar consequences. In contrast to previous reports, we found that inhibition of Ca2+ signaling by BAPTA-AM did not specifically inhibit DoG transcription but globally impaired transcription. Most importantly, HSV-1-induced DoTT, but not stress-induced DoG transcription, was accompanied by a strong increase in open chromatin downstream of the affected poly(A) sites. In its extent and kinetics, downstream open chromatin essentially matched the poly(A) read-through transcription. We show that this does not cause but rather requires DoTT as well as high levels of transcription into the genomic regions downstream of genes. This raises intriguing new questions regarding the role of histone repositioning in the wake of RNA Polymerase II passage downstream of impaired poly(A) site recognition.

Following publication of the original article [1], the authors have reported that Fig. 2 and Additional file 1: Figure S1, S2 partially show red scripts.

RNA-seq reads containing part of the poly(A) tail of transcripts (denoted as poly(A) reads) provide the most direct evidence for the position of poly(A) sites in the genome. However, due to reduced coverage of poly(A) tails by reads, poly(A) reads are not routinely identified during RNA-seq mapping. Nevertheless, recent studies for several herpesviruses successfully employed mapping of poly(A) reads to identify herpesvirus poly(A) sites using different strategies and customized programs. To more easily allow such analyses without requiring additional programs, we integrated poly(A) read mapping and prediction of poly(A) sites into our RNA-seq mapping program ContextMap 2. The implemented approach essentially generalizes previously used poly(A) read mapping approaches and combines them with the context-based approach of ContextMap 2 to take into account information provided by other reads aligned to the same location. Poly(A) read mapping using ContextMap 2 was evaluated on real-life data from the ENCODE project and compared against a competing approach based on transcriptome assembly (KLEAT). This showed high positive predictive value for our approach, evidenced also by the presence of poly(A) signals, and considerably lower runtime than KLEAT. Although sensitivity is low for both methods, we show that this is in part due to a high extent of spurious results in the gold standard set derived from RNA-PET data. Sensitivity improves for poly(A) sites of known transcripts or determined with a more specific poly(A) sequencing protocol and increases with read coverage on transcript ends. Finally, we illustrate the usefulness of the approach in a high read coverage scenario by a re-analysis of published data for herpes simplex virus 1. Thus, with current trends towards increasing sequencing depth and read length, poly(A) read mapping will prove to be increasingly useful and can now be performed automatically during RNA-seq mapping with ContextMap 2.

Herpes simplex virus type 1 (HSV-1) is a neurotropic virus causing vesicular oral or genital skin lesions, meningitis and other diseases particularly harmful in immunocompromised individuals. To comprehensively investigate the complex interaction between HSV-1 and its host we combined two genome-scale screens for host factors (HFs) involved in virus replication. A yeast two-hybrid screen for protein interactions and a RNA interference (RNAi) screen with a druggable genome small interfering RNA (siRNA) library confirmed existing and identified novel HFs which functionally influence HSV-1 infection. Bioinformatic analyses found the 358 HFs were enriched for several pathways and multi-protein complexes. Of particular interest was the identification of Med23 as a strongly anti-viral component of the largely pro-viral Mediator complex, which links specific transcription factors to RNA polymerase II. The anti-viral effect of Med23 on HSV-1 replication was confirmed in gain-of-function gene overexpression experiments, and this inhibitory effect was specific to HSV-1, as a range of other viruses including Vaccinia virus and Semliki Forest virus were unaffected by Med23 depletion. We found Med23 significantly upregulated expression of the type III interferon family (IFN-λ) at the mRNA and protein level by directly interacting with the transcription factor IRF7. The synergistic effect of Med23 and IRF7 on IFN-λ induction suggests this is the major transcription factor for IFN-λ expression. Genotypic analysis of patients suffering recurrent orofacial HSV-1 outbreaks, previously shown to be deficient in IFN-λ secretion, found a significant correlation with a single nucleotide polymorphism in the IFN-λ3 (IL28b) promoter strongly linked to Hepatitis C disease and treatment outcome. This paper describes a link between Med23 and IFN-λ, provides evidence for the crucial role of IFN-λ in HSV-1 immune control, and highlights the power of integrative genome-scale approaches to identify HFs critical for disease progression and outcome.

Somatic transposon mutagenesis in mice is an efficient strategy to investigate the genetic mechanisms of tumorigenesis. The identification of tumor driving transposon insertions traditionally requires the generation of large tumor cohorts to obtain information about common insertion sites. Tumor driving insertions are also characterized by their clonal expansion in tumor tissue, a phenomenon that is facilitated by the slow and evolving transformation process of transposon mutagenesis. We describe here an improved approach for the detection of tumor driving insertions that assesses the clonal expansion of insertions by quantifying the relative proportion of sequence reads obtained in individual tumors. To this end, we have developed a protocol for insertion site sequencing that utilizes acoustic shearing of tumor DNA and Illumina sequencing. We analyzed various solid tumors generated by PiggyBac mutagenesis and for each tumor >10⁶ reads corresponding to >10⁴ insertion sites were obtained. In each tumor, 9 to 25 insertions stood out by their enriched sequence read frequencies when compared to frequencies obtained from tail DNA controls. These enriched insertions are potential clonally expanded tumor driving insertions, and thus identify candidate cancer genes. The candidate cancer genes of our study comprised many established cancer genes, but also novel candidate genes such as Mastermind-like1 (Mamld1) and Diacylglycerolkinase delta (Dgkd). We show that clonal expansion analysis by high-throughput sequencing is a robust approach for the identification of candidate cancer genes in insertional mutagenesis screens on the level of individual tumors.

RNA sequencing (RNA-seq) provides novel opportunities for transcriptomic studies at nucleotide resolution, including transcriptomics of viruses or microbes infecting a cell. However, standard approaches for mapping the resulting sequencing reads generally ignore alternative sources of expression other than the host cell and are little equipped to address the problems arising from redundancies and gaps among sequenced microbe and virus genomes. We show that screening of sequencing reads for contaminations and infections can be performed easily using ContextMap, our recently developed mapping software. Based on mapping-derived statistics, mapping confidence, similarities and misidentifications (e.g. due to missing genome sequences) of species/strains can be assessed. Performance of our approach is evaluated on three real-life sequencing data sets and compared to state-of-the-art metagenomics tools. In particular, ContextMap vastly outperformed GASiC and GRAMMy in terms of runtime. In contrast to MEGAN4, it was capable of providing individual read mappings to species and resolving non-unique mappings, thus allowing the identification of misalignments caused by sequence similarities between genomes and missing genome sequences. Our study illustrates the importance and potentials of routinely mining RNA-seq experiments for infections or contaminations by microbes and viruses. By using ContextMap, gene expression of infecting agents can be analyzed and novel insights in infection processes and tumorigenesis can be obtained.