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. Author summary: Merkel cell polyomavirus (MCPyV) is the only human polyomavirus that causes cancer in humans. As with all human polyomaviruses, the available infection models are limited. Thus, many processes such as the host response to infection and its regulation by the virus to establish infection and persistence are poorly understood. To better understand this interplay of viral MCPyV proteins, we performed genome-wide transcriptome and chromatin studies in primary human fibroblasts and simulated infection and tumorigenesis conditions by ectopically expressing the early viral proteins individually or in combination at different time points. This allowed us to uncover a novel, previously undescribed function of polyomavirus sT, namely the reduction of the ISG response by affecting the ISGF3 complex, specifically by reducing IRF9 protein levels. This work sheds light on how early viral proteins influence the type I IFN response and how their interplay may affect MCPyV infection, persistence, and MCC progression. [ABSTRACT FROM AUTHOR]

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. Author summary: We conducted a comprehensive characterization and reannotation of murine cytomegalovirus (MCMV) gene expression during lytic infection in murine fibroblasts using an integrative multi-omics approach. This unveiled hundreds of novel transcripts that explained the expression of close to 300 so far unknown viral open reading frames (ORFs). Interestingly, small viral ORFs (sORFs) were amongst the most highly expressed viral gene products and thus presumably encode for important viral microproteins of unknown function. We show that sORFs located upstream of larger ORFs tune the expression of the downstream ORFs by repressing their translation. We classified viral transcription start sites (TiSS) based on their expression kinetics obtained by the combination of metabolic RNA labelling with transcription start sites profiling. This not only identified a so far unknown viral immediate-early transcript (ie4, m166.5 RNA) but also revealed a novel class of viral late transcripts that are expressed later than canonical true late genes and lack TATA box-like motifs. We exemplify for the m145 locus how so far unknown viral TiSS give rise to abundantly expressed truncated viral proteins. In summary, we provide a state-of-the-art annotation of an important model virus, which will be instrumental for future studies on CMV biology, immunology and pathogenesis. [ABSTRACT FROM AUTHOR]

Motivation To date, no methods are available for the targeted identification of genomic subregions with differences in sequencing read distributions between two conditions. Existing approaches either only determine absolute read number changes, require predefined subdivisions of input windows or average across multiple genes. Results Here, we present RegCFinder, which automatically identifies subregions of input windows with differences in read density between two conditions. For this purpose, the problem is defined as an instance of the all maximum scoring subsequences problem, which can be solved in linear time. Subsequently, statistical significance and differential usage of identified subregions are determined with DEXSeq. RegCFinder allows flexible definition of input windows to target the analysis to any regions of interests, e.g. promoters, gene bodies, peak regions and more. Furthermore, any type of sequencing assay can be used as input; thus, RegCFinder lends itself to a wide range of applications. We illustrate the usefulness of RegCFinder on two applications, where we can both confirm previous results and identify interesting gene subgroups with distinctive changes in read distributions. Availability and implementation RegCFinder is implemented as a workflow for the workflow management system Watchdog and available at: https://github.com/watchdog-wms/watchdog-wms-workflows/ Supplementary information Supplementary data are available at Bioinformatics Advances online. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 FROM AUTHOR]

Executing tumor cell death: There is increasing interest in developing cancer immunotherapies that target the innate immune pathways regulating cytokine production and cell death, but the interplay between these two closely connected processes is not well understood. In mouse and human cancer cell lines, Boehmer et al. demonstrate that cytokine production and apoptosis induced by retinoic acid–inducible gene I (RIG-I) ligands, including 5′-triphosphate RNA (3p-RNA), are two separable events in which RIG-I is required for production of type I interferon but not execution of apoptosis. Mass spectrometry and loss-of-function assays showed that 3p-RNA directly activates OAS1 and RNase L, which promoted translational arrest and depletion of antiapoptotic MCL-1. These results demonstrate that RIG-I–mediated apoptosis involves priming and effector stages, reminiscent of inflammasome activation, both of which could serve as potential targets for cancer immunotherapy. Cytoplasmic double-stranded RNA is sensed by RIG-I–like receptors (RLRs), leading to induction of type I interferons (IFN-Is), proinflammatory cytokines, and apoptosis. Here, we elucidate signaling mechanisms that lead to cytokine secretion and cell death induction upon stimulation with the bona fide RIG-I ligand 5′-triphosphate RNA (3p-RNA) in tumor cells. We show that both outcomes are mediated by dsRNA-receptor families with RLR being essential for cytokine production and IFN-I–mediated priming of effector pathways but not for apoptosis. Affinity purification followed by mass spectrometry and subsequent functional analysis revealed that 3p-RNA bound and activated oligoadenylate synthetase 1 and RNase L. RNase L–deficient cells were profoundly impaired in their ability to undergo apoptosis. Mechanistically, the concerted action of translational arrest triggered by RNase L and up-regulation of NOXA was needed to deplete the antiapoptotic MCL-1 to cause intrinsic apoptosis. Thus, 3p-RNA–induced apoptosis is a two-step process consisting of RIG-I–dependent priming and an RNase L–dependent effector phase. [ABSTRACT FROM AUTHOR]

Abstract Background A recent large-scale analysis of Gene Expression Omnibus (GEO) data found frequent evidence for spatial defects in a substantial fraction of Affymetrix microarrays in the GEO. Nevertheless, in contrast to quality assessment, artefact detection is not widely used in standard gene expression analysis pipelines. Furthermore, although approaches have been proposed to detect diverse types of spatial noise on arrays, the correction of these artefacts is mostly left to either summarization methods or the corresponding arrays are completely discarded. Results We show that state-of-the-art robust summarization procedures are vulnerable to artefacts on arrays and cannot appropriately correct for these. To address this problem, we present a simple approach to detect artefacts with high recall and precision, which we further improve by taking into account the spatial layout of arrays. Finally, we propose two correction methods for these artefacts that either substitute values of defective probes using probeset information or filter corrupted probes. We show that our approach can identify and correct defective probe measurements appropriately and outperforms existing tools. Conclusions While summarization is insufficient to correct for defective probes, this problem can be addressed in a straightforward way by the methods we present for identification and correction of defective probes. As these methods output CEL files with corrected probe values that serve as input to standard normalization and summarization procedures, they can be easily integrated into existing microarray analysis pipelines as an additional pre-processing step. An R package is freely available from http://www.bio.ifi.lmu.de/artefact-correction.

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. Author summary: Viral infections profoundly alter host cell gene expression. It is important to understand both how viruses hijack the host cell machineries to express their own genes and how host cells respond to viral infection for defense. We have previously shown that herpes simplex virus-1 (HSV-1) blocks host cell transcription termination, at least in part, through the viral immediate early protein ICP27. Here we show that HSV-1 infection also alters mRNA 3' end formation and promotes the formation of truncated mRNAs. Some of these aberrant mRNAs are exported into the cytoplasm and translated. This viral activity requires ICP27 and other viral factors. Our study, together with other recent reports, suggests that viral infections and cellular stress elicit similar responses in mammalian cells. [ABSTRACT FROM AUTHOR]

Abstract Background Stochastic simulation can be used to illustrate the development of biological systems over time and the stochastic nature of these processes. Currently available programs for stochastic simulation, however, are limited in that they either a) do not provide the most efficient simulation algorithms and are difficult to extend, b) cannot be easily integrated into other applications or c) do not allow to monitor and intervene during the simulation process in an easy and intuitive way. Thus, in order to use stochastic simulation in innovative high-level modeling and analysis approaches more flexible tools are necessary. Results In this article, we present FERN (Framework for Evaluation of Reaction Networks), a Java framework for the efficient simulation of chemical reaction networks. FERN is subdivided into three layers for network representation, simulation and visualization of the simulation results each of which can be easily extended. It provides efficient and accurate state-of-the-art stochastic simulation algorithms for well-mixed chemical systems and a powerful observer system, which makes it possible to track and control the simulation progress on every level. To illustrate how FERN can be easily integrated into other systems biology applications, plugins to Cytoscape and CellDesigner are included. These plugins make it possible to run simulations and to observe the simulation progress in a reaction network in real-time from within the Cytoscape or CellDesigner environment. Conclusion FERN addresses shortcomings of currently available stochastic simulation programs in several ways. First, it provides a broad range of efficient and accurate algorithms both for exact and approximate stochastic simulation and a simple interface for extending to new algorithms. FERN's implementations are considerably faster than the C implementations of gillespie2 or the Java implementations of ISBJava. Second, it can be used in a straightforward way both as a stand-alone program and within new systems biology applications. Finally, complex scenarios requiring intervention during the simulation progress can be modelled easily with FERN.

Abstract Background The existence of negative correlations between degrees of interacting proteins is being discussed since such negative degree correlations were found for the large-scale yeast protein-protein interaction (PPI) network of Ito et al. More recent studies observed no such negative correlations for high-confidence interaction sets. In this article, we analyzed a range of experimentally derived interaction networks to understand the role and prevalence of degree correlations in PPI networks. We investigated how degree correlations influence the structure of networks and their tolerance against perturbations such as the targeted deletion of hubs. Results For each PPI network, we simulated uncorrelated, positively and negatively correlated reference networks. Here, a simple model was developed which can create different types of degree correlations in a network without changing the degree distribution. Differences in static properties associated with degree correlations were compared by analyzing the network characteristics of the original PPI and reference networks. Dynamics were compared by simulating the effect of a selective deletion of hubs in all networks. Conclusion Considerable differences between the network types were found for the number of components in the original networks. Negatively correlated networks are fragmented into significantly less components than observed for positively correlated networks. On the other hand, the selective deletion of hubs showed an increased structural tolerance to these deletions for the positively correlated networks. This results in a lower rate of interaction loss in these networks compared to the negatively correlated networks and a decreased disintegration rate. Interestingly, real PPI networks are most similar to the randomly correlated references with respect to all properties analyzed. Thus, although structural properties of networks can be modified considerably by degree correlations, biological PPI networks do not actually seem to make use of this possibility.

Abstract Background Although protein-protein interaction networks determined with high-throughput methods are incomplete, they are commonly used to infer the topology of the complete interactome. These partial networks often show a scale-free behavior with only a few proteins having many and the majority having only a few connections. Recently, the possibility was suggested that this scale-free nature may not actually reflect the topology of the complete interactome but could also be due to the error proneness and incompleteness of large-scale experiments. Results In this paper, we investigate the effect of limited sampling on average clustering coefficients and how this can help to more confidently exclude possible topology models for the complete interactome. Both analytical and simulation results for different network topologies indicate that partial sampling alone lowers the clustering coefficient of all networks tremendously. Furthermore, we extend the original sampling model by also including spurious interactions via a preferential attachment process. Simulations of this extended model show that the effect of wrong interactions on clustering coefficients depends strongly on the skewness of the original topology and on the degree of randomness of clustering coefficients in the corresponding networks. Conclusion Our findings suggest that the complete interactome is either highly skewed such as e.g. in scale-free networks or is at least highly clustered. Although the correct topology of the interactome may not be inferred beyond any reasonable doubt from the interaction networks available, a number of topologies can nevertheless be excluded with high confidence.

Herpes simplex virus 1 (HSV-1) induces a profound host shutoff during lytic infection. The virion host shutoff (vhs) protein plays a key role in this process by efficiently cleaving host and viral mRNAs. Furthermore, the onset of viral DNA replication is accompanied by a rapid decline in host transcriptional activity. To dissect relative contributions of both mechanisms and elucidate gene-specific host transcriptional responses throughout the first 8 h of lytic HSV-1 infection, we used transcriptome sequencing of total, newly transcribed (4sU-labeled) and chromatinassociated RNA in wild-type (WT) and Δvhs mutant infection of primary human fibroblasts. Following virus entry, vhs activity rapidly plateaued at an elimination rate of around 30% of cellular mRNAs per hour until 8 h postinfection (p.i.). In parallel, host transcriptional activity dropped to 10 to 20%. While the combined effects of both phenomena dominated infection-induced changes in total RNA, extensive gene-specific transcriptional regulation was observable in chromatin-associated RNA and was surprisingly concordant between WT and Δvhs infections. Both induced strong transcriptional upregulation of a small subset of genes that were poorly expressed prior to infection but already primed by H3K4me3 histone marks at their promoters. Most interestingly, analysis of chromatin-associated RNA revealed vhsnuclease-activity-dependent transcriptional downregulation of at least 150 cellular genes, in particular of many integrin adhesome and extracellular matrix components. This was accompanied by a vhs-dependent reduction in protein levels by 8 h p.i. for many of these genes. In summary, our study provides a comprehensive picture of the molecular mechanisms that govern cellular RNA metabolism during the first 8 h of lytic HSV-1 infection. IMPORTANCE The HSV-1 virion host shutoff (vhs) protein efficiently cleaves both host and viral mRNAs in a translation-dependent manner. In this study, we model and quantify changes in vhs activity, as well as virus-induced global loss of host transcriptional activity, during productive HSV-1 infection. In general, HSV-1-induced alterations in total RNA levels were dominated by these two global effects. In contrast, chromatin-associated RNA depicted gene-specific transcriptional changes. This revealed highly concordant transcriptional changes in WT and -vhs infections, confirmed DUX4 as a key transcriptional regulator in HSV-1 infection, and identified vhsdependent transcriptional downregulation of the integrin adhesome and extracellular matrix components. The latter explained seemingly gene-specific effects previously attributed to vhs-mediated mRNA degradation and resulted in a concordant loss in protein levels by 8 h p.i. for many of the respective genes. [ABSTRACT FROM AUTHOR]

Morphogenesis of herpesviral virions is initiated in the nucleus but completed in the cytoplasm. Mature virions contain more than 25 tegument proteins many of which perform both nuclear and cytoplasmic functions suggesting they shuttle between these compartments. While nuclear import of herpesviral proteins was shown to be crucial for viral propagation, active nuclear export and its functional impact are still poorly understood. To systematically analyze nuclear export of tegument proteins present in virions of Herpes simplex virus type 1 (HSV1) and Epstein‐Barr virus (EBV), the Nuclear EXport Trapped by RAPamycin (NEX‐TRAP) was applied. Nine of the 22 investigated HSV1 tegument proteins including pUL4, pUL7, pUL11, pUL13, pUL21, pUL37d11, pUL47, pUL48 and pUS2 as well as 2 out of 6 EBV orthologs harbor nuclear export activity. A functional leucine‐rich nuclear export sequence (NES) recognized by the export factor CRM1/Xpo1 was identified in six of them. The comparison between experimental and bioinformatic data indicates that experimental validation of predicted NESs is required. Mutational analysis of the pUL48/VP16 NES revealed its importance for herpesviral propagation. Together our data suggest that nuclear export is an important feature of the herpesviral life cycle required to co‐ordinate nuclear and cytoplasmic processes. Herpesviral virions contain numerous tegument proteins with important functions both in nucleus and cytoplasm suggesting they shuttle between these compartments. Unlike nuclear import, nuclear export of tegument proteins is poorly understood. Using the Nuclear EXport Trapped by RAPamycin (NEX‐TRAP) assay, nuclear export activity and the relevant export sequences were identified in several Herpes simplex virus type 1 tegument proteins and their Epstein‐Barr virus orthologs. These results indicate an important role of nuclear export for herpesviral propagation and form the basis for future functional analysis. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Background: Mapping of short sequencing reads is a crucial step in the analysis of RNA sequencing (RNA-seq) data. ContextMap is an RNA-seq mapping algorithm that uses a context-based approach to identify the best alignment for each read and allows parallel mapping against several reference genomes. Results: In this article, we present ContextMap 2, a new and improved version of ContextMap. Its key novel features are: (i) a plug-in structure that allows easily integrating novel short read alignment programs with improved accuracy and runtime; (ii) context-based identification of insertions and deletions (indels); (iii) mapping of reads spanning an arbitrary number of exons and indels. ContextMap 2 using Bowtie, Bowtie 2 or BWA was evaluated on both simulated and real-life data from the recently published RGASP study. Conclusions: We show that ContextMap 2 generally combines similar or higher recall compared to other state-of-the-art approaches with significantly higher precision in read placement and junction and indel prediction. Furthermore, runtime was significantly lower than for the best competing approaches. ContextMap 2 is freely available at http://www.bio.ifi.lmu.de/ContextMap. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

How different integrins that bind to the same type of extracellular matrix protein mediate specific functions is unclear. We report the functional analysis of β1- and αv-class integrins expressed in pan-integrin-null fibroblasts seeded on fibronectin. Reconstitution with β1-class integrins promotes myosin-II-independent formation of small peripheral adhesions and cell protrusions, whereas expression of αv-class integrins induces the formation of large focal adhesions. Co-expression of both integrin classes leads to full myosin activation and traction-force development on stiff fibronectin-coated substrates, with αv-class integrins accumulating in adhesion areas exposed to high traction forces. Quantitative proteomics linked αv-class integrins to a GEF-H1-RhoA pathway coupled to the formin mDia1 but not myosin II, and α5β1 integrins to a RhoA-Rock-myosin II pathway. Our study assigns specific functions to distinct fibronectin-binding integrins, demonstrating that α5β1integrins accomplish force generation, whereas αv-class integrins mediate the structural adaptations to forces, which cooperatively enable cells to sense the rigidity of fibronectin-based microenvironments. [ABSTRACT FROM AUTHOR]

Transcriptome sequencing (RNA-Seq) overcomes limitations of previously used RNA quantification methods and provides one experimental framework for both high-throughput characterization and quantification of transcripts at the nucleotide level. The first step and a major challenge in the analysis of such experiments is the mapping of sequencing reads to a transcriptomic origin including the identification of splicing events. In recent years, a large number of such mapping algorithms have been developed, all of which have in common that they require algorithms for aligning a vast number of reads to genomic or transcriptomic sequences. Although the FM-index based aligner Bowtie has become a de facto standard within mapping pipelines, a much larger number of possible alignment algorithms have been developed also including other variants of FM-index based aligners. Accordingly, developers and users of RNA-seq mapping pipelines have the choice among a large number of available alignment algorithms. To provide guidance in the choice of alignment algorithms for these purposes, we evaluated the performance of 14 widely used alignment programs from three different algorithmic classes: algorithms using either hashing of the reference transcriptome, hashing of reads, or a compressed FMindex representation of the genome. Here, special emphasis was placed on both precision and recall and the performance for different read lengths and numbers of mismatches and indels in a read. Our results clearly showed the significant reduction in memory footprint and runtime provided by FM-index based aligners at a precision and recall comparable to the best hash table based aligners. Furthermore, the recently developed Bowtie 2 alignment algorithm shows a remarkable tolerance to both sequencing errors and indels, thus, essentially making hash-based aligners obsolete. [ABSTRACT FROM AUTHOR]

During viral infections cellular gene expression is subject to rapid alterations induced by both viral and antiviral mechanisms. In this study, we applied metabolic labeling of newly transcribed RNA with 4-thiouridine (4sU-tagging) to dissect the real-time kinetics of cellular and viral transcriptional activity during lytic murine cytomegalovirus (MCMV) infection. Microarray profiling on newly transcribed RNA obtained at different times during the first six hours of MCMV infection revealed discrete functional clusters of cellular genes regulated with distinct kinetics at surprising temporal resolution. Immediately upon virus entry, a cluster of NF-κB- and interferon-regulated genes was induced. Rapid viral counter-regulation of this coincided with a very transient DNA-damage response, followed by a delayed ER-stress response. Rapid counter-regulation of all three clusters indicated the involvement of novel viral regulators targeting these pathways. In addition, down-regulation of two clusters involved in cell-differentiation (rapid repression) and cell-cycle (delayed repression) was observed. Promoter analysis revealed all five clusters to be associated with distinct transcription factors, of which NF-κB and c-Myc were validated to precisely match the respective transcriptional changes observed in newly transcribed RNA. 4sU-tagging also allowed us to study the real-time kinetics of viral gene expression in the absence of any interfering virion-associated-RNA. Both qRT-PCR and next-generation sequencing demonstrated a sharp peak of viral gene expression during the first two hours of infection including transcription of immediate-early, early and even well characterized late genes. Interestingly, this was subject to rapid gene silencing by 5-6 hours post infection. Despite the rapid increase in viral DNA load during viral DNA replication, transcriptional activity of some viral genes remained remarkably constant until late-stage infection, or was subject to further continuous decline. In summary, this study pioneers real-time transcriptional analysis during a lytic herpesvirus infection and highlights numerous novel regulatory aspects of virus-host-cell interaction. [ABSTRACT FROM AUTHOR]

Background: Sequencing of mRNA (RNA-seq) by next generation sequencing technologies is widely used for analyzing the transcriptomic state of a cell. Here, one of the main challenges is the mapping of a sequenced read to its transcriptomic origin. As a simple alignment to the genome will fail to identify reads crossing splice junctions and a transcriptome alignment will miss novel splice sites, several approaches have been developed for this purpose. Most of these approaches have two drawbacks. First, each read is assigned to a location independent on whether the corresponding gene is expressed or not, i.e. information from other reads is not taken into account. Second, in case of multiple possible mappings, the mapping with the fewest mismatches is usually chosen which may lead to wrong assignments due to sequencing errors. Results: To address these problems, we developed ContextMap which efficiently uses information on the context of a read, i.e. reads mapping to the same expressed region. The context information is used to resolve possible ambiguities and, thus, a much larger degree of ambiguities can be allowed in the initial stage in order to detect all possible candidate positions. Although ContextMap can be used as a stand-alone version using either a genome or transcriptome as input, the version presented in this article is focused on refining initial mappings provided by other mapping algorithms. Evaluation results on simulated sequencing reads showed that the application of ContextMap to either TopHat or MapSplice mappings improved the mapping accuracy of both initial mappings considerably. Conclusions: In this article, we show that the context of reads mapping to nearby locations provides valuable information for identifying the best unique mapping for a read. Using our method, mappings provided by other state-of-the-art methods can be refined and alignment accuracy can be further improved. Availability: http://www.bio.ifi.lmu.de/ContextMap. [ABSTRACT FROM AUTHOR]

In mammalian cells, AU-rich elements (AREs) are well known regulatory sequences located in the 3′ untranslated region (UTR) of many short-lived mRNAs. AREs cause mRNAs to be degraded rapidly and thereby suppress gene expression at the posttranscriptional level. Based on the number of AUUUA pentamers, their proximity, and surrounding AU-rich regions, we generated an algorithm termed AREScore that identifies AREs and provides a numerical assessment of their strength. By analyzing the AREScore distribution in the transcriptomes of 14 metazoan species, we provide evidence that AREs were selected for in several vertebrates and Drosophila melanogaster. We then measured mRNA expression levels genome-wide to address the importance of AREs in SL2 cells derived from D. melanogaster hemocytes. Tis11, a zinc finger RNA--binding protein homologous to mammalian tristetraprolin, was found to target ARE--containing reporter mRNAs for rapid degradation in SL2 cells. Drosophila mRNAs whose expression is elevated upon knock down of Tis11 were found to have higher AREScores. Moreover high AREScores correlate with reduced mRNA expression levels on a genome-wide scale. The precise measurement of degradation rates for 26 Drosophila mRNAs revealed that the AREScore is a very good predictor of short-lived mRNAs. Taken together, this study introduces AREScore as a simple tool to identify ARE--containing mRNAs and provides compelling evidence that AREs are widespread regulatory elements in Drosophila. [ABSTRACT FROM AUTHOR]

Coronaviruses (CoVs) are important human and animal pathogens that induce fatal respiratory, gastrointestinal and neurological disease. The outbreak of the severe acute respiratory syndrome (SARS) in 2002/2003 has demonstrated human vulnerability to (Coronavirus) CoV epidemics. Neither vaccines nor therapeutics are available against human and animal CoVs. Knowledge of host cell proteins that take part in pivotal virus-host interactions could define broad-spectrum antiviral targets. In this study, we used a systems biology approach employing a genome-wide yeast-two hybrid interaction screen to identify immunopilins (PPIA, PPIB, PPIH, PPIG, FKBP1A, FKBP1B) as interaction partners of the CoV non-structural protein 1 (Nsp1). These molecules modulate the Calcineurin/NFAT pathway that plays an important role in immune cell activation. Overexpression of NSP1 and infection with live SARS-CoV strongly increased signalling through the Calcineurin/NFAT pathway and enhanced the induction of interleukin 2, compatible with late-stage immunopathogenicity and long-term cytokine dysregulation as observed in severe SARS cases. Conversely, inhibition of cyclophilins by cyclosporine A (CspA) blocked the replication of CoVs of all genera, including SARS-CoV, human CoV-229E and -NL-63, feline CoV, as well as avian infectious bronchitis virus. Non-immunosuppressive derivatives of CspA might serve as broad-range CoV inhibitors applicable against emerging CoVs as well as ubiquitous pathogens of humans and livestock. [ABSTRACT FROM AUTHOR]

A characteristic of integrins is their ability to transfer chemical and mechanical signals across the plasma membrane. Force generated by myosin II makes cells able to sense substrate stiffness and induce maturation of nascent adhesions into focal adhesions. In this paper, we present a comprehensive proteomic analysis of nascent and mature adhesions. The purification of integrin adhesion complexes combined with quantitative mass spectrometry enabled the identification and quantification of known and new adhesion-associated proteins. Furthermore, blocking adhesion maturation with the myosin II inhibitor blebbistatin markedly impaired the recruitment of LIM domain proteins to integrin adhesion sites. This suggests a common recruitment mechanism for a whole class of adhesion-associated proteins, involving myosin II and the zinc-finger-type LIM domain. [ABSTRACT FROM AUTHOR]

Summary: The mRNA targets of microRNAs (miRNAs) can be identified by immunoprecipitation of Argonaute (Ago) protein-containing RNA-induced silencing complexes (RISCs) followed by microarray analysis (RIP-Chip). Here we used Ago2-based RIP-Chip to identify transcripts targeted by Kaposi''s sarcoma-associated herpesvirus (KSHV) miRNAs (n = 114), Epstein-Barr virus (EBV) miRNAs (n = 44), and cellular miRNAs (n = 2337) in six latently infected or stably transduced human B cell lines. Of the six KSHV miRNA targets chosen for validation, four showed regulation via their 3′UTR, while two showed regulation via binding sites within coding sequences. Two genes governing cellular transport processes (TOMM22 and IPO7) were confirmed to be targeted by EBV miRNAs. A significant number of viral miRNA targets were upregulated in infected cells, suggesting that viral miRNAs preferentially target cellular genes induced upon infection. Transcript half-life both of cellular and viral miRNA targets negatively correlated with recruitment to RISC complexes, indicating that RIP-Chip offers a quantitative estimate of miRNA function. [Copyright &y& Elsevier]

Herpesviruses constitute a family of large DNA viruses widely spread in vertebrates and causing a variety of different diseases. They possess dsDNA genomes ranging from 120 to 240 kbp encoding between 70 to 170 open reading frames. We previously reported the protein interaction networks of two herpesviruses, varicella-zoster virus (VZV) and Kaposi's sarcoma-associated herpesvirus (KSHV). In this study, we systematically tested three additional herpesvirus species, herpes simplex virus 1 (HSV-1), murine cytomegalovirus and Epstein-Barr virus, for protein interactions in order to be able to perform a comparative analysis of all three herpesvirus subfamilies. We identified 735 interactions by genome-wide yeast-two-hybrid screens (Y2H), and, together with the interactomes of VZV and KSHV, included a total of 1,007 intraviral protein interactions in the analysis. Whereas a large number of interactions have not been reported previously, we were able to identify a core set of highly conserved protein interactions, like the interaction between HSV-1 UL33 with the nuclear egress proteins UL31/UL34. Interactions were conserved between orthologous proteins despite generally low sequence similarity, suggesting that function may be more conserved than sequence. By combining interactomes of different species we were able to systematically address the low coverage of the Y2H system and to extract biologically relevant interactions which were not evident from single species. [ABSTRACT FROM AUTHOR]