Methods to study splicing from high-throughput RNA Sequencing data

The development of novel high-throughput sequencing (HTS) methods for RNA (RNA-Seq) has provided a very powerful mean to study splicing under multiple conditions at unprecedented depth. However, the complexity of the information to be analyzed has turned this into a challenging task. In the last few years, a plethora of tools have been developed, allowing researchers to process RNA-Seq data to study the expression of isoforms and splicing events, and their relative changes under different conditions. We provide an overview of the methods available to study splicing from short RNA-Seq data. We group the methods according to the different questions they address: 1) Assignment of the sequencing reads to their likely gene of origin. This is addressed by methods that map reads to the genome and/or to the available gene annotations. 2) Recovering the sequence of splicing events and isoforms. This is addressed by transcript reconstruction and de novo assembly methods. 3) Quantification of events and isoforms. Either after reconstructing transcripts or using an annotation, many methods estimate the expression level or the relative usage of isoforms and/or events. 4) Providing an isoform or event view of differential splicing or expression. These include methods that compare relative event/isoform abundance or isoform expression across two or more conditions. 5) Visualizing splicing regulation. Various tools facilitate the visualization of the RNA-Seq data in the context of alternative splicing. In this review, we do not describe the specific mathematical models behind each method. Our aim is rather to provide an overview that could serve as an entry point for users who need to decide on a suitable tool for a specific analysis. We also attempt to propose a classification of the tools according to the operations they do, to facilitate the comparison and choice of methods.Comment: 31 pages, 1 figure, 9 tables. Small corrections adde

Targeting transcription regulation in cancer with a covalent CDK7 inhibitor

Tumour oncogenes include transcription factors that co-opt the general transcriptional machinery to sustain the oncogenic state, but direct pharmacological inhibition of transcription factors has so far proven difficult. However, the transcriptional machinery contains various enzymatic cofactors that can be targeted for the development of new therapeutic candidates, including cyclin-dependent kinases (CDKs). Here we present the discovery and characterization of a covalent CDK7 inhibitor, THZ1, which has the unprecedented ability to target a remote cysteine residue located outside of the canonical kinase domain, providing an unanticipated means of achieving selectivity for CDK7. Cancer cell-line profiling indicates that a subset of cancer cell lines, including human T-cell acute lymphoblastic leukaemia (T-ALL), have exceptional sensitivity to THZ1. Genome-wide analysis in Jurkat T-ALL cells shows that THZ1 disproportionally affects transcription of RUNX1 and suggests that sensitivity to THZ1 may be due to vulnerability conferred by the RUNX1 super-enhancer and the key role of RUNX1 in the core transcriptional regulatory circuitry of these tumour cells. Pharmacological modulation of CDK7 kinase activity may thus provide an approach to identify and treat tumour types that are dependent on transcription for maintenance of the oncogenic state.National Institutes of Health (U.S.) (Grant HG002668)National Institutes of Health (U.S.) (Grant CA109901

Recurrent somatic mutations in POLR2A define a distinct subset of meningiomas

RNA polymerase II mediates the transcription of all protein-coding genes in eukaryotic cells, a process that is fundamental to life. Genomic mutations altering this enzyme have not previously been linked to any pathology in humans, which is a testament to its indispensable role in cell biology. On the basis of a combination of next-generation genomic analyses of 775 meningiomas, we report that recurrent somatic p.Gln403Lys or p.Leu438_His439del mutations in POLR2A, which encodes the catalytic subunit of RNA polymerase II (ref. 1), hijack this essential enzyme and drive neoplasia. POLR2A mutant tumors show dysregulation of key meningeal identity genes including WNT6 and ZIC1/ZIC4. In addition to mutations in POLR2A, NF2, SMARCB1, TRAF7, KLF4, AKT1, PIK3CA, and SMO4 we also report somatic mutations in AKT3, PIK3R1, PRKAR1A, and SUFU in meningiomas. Our results identify a role for essential transcriptional machinery in driving tumorigenesis and define mutually exclusive meningioma subgroups with distinct clinical and pathological features

Metabolism before, during and after anaesthesia in colic and healthy horses

<p>Abstract</p> <p>Background</p> <p>Many colic horses are compromised due to the disease state and from hours of starvation and sometimes long trailer rides. This could influence their muscle energy reserves and affect the horses' ability to recover. The principal aim was to follow metabolic parameter before, during, and up to 7 days after anaesthesia in healthy horses and in horses undergoing abdominal surgery due to colic.</p> <p>Methods</p> <p>20 healthy horses given anaesthesia alone and 20 colic horses subjected to emergency abdominal surgery were anaesthetised for a mean of 228 minutes and 183 minutes respectively. Blood for analysis of haematology, electrolytes, cortisol, creatine kinase (CK), free fatty acids (FFA), glycerol, glucose and lactate was sampled before, during, and up to 7 days after anaesthesia. Arterial and venous blood gases were obtained before, during and up to 8 hours after recovery. Gluteal muscle biopsy specimens for biochemical analysis of muscle metabolites were obtained at start and end of anaesthesia and 1 h and 1 day after recovery.</p> <p>Results</p> <p>Plasma cortisol, FFA, glycerol, glucose, lactate and CK were elevated and serum phosphate and potassium were lower in colic horses before anaesthesia. Muscle adenosine triphosphate (ATP) content was low in several colic horses. Anaesthesia and surgery resulted in a decrease in plasma FFA and glycerol in colic horses whereas levels increased in healthy horses. During anaesthesia muscle and plasma lactate and plasma phosphate increased in both groups. In the colic horses plasma lactate increased further after recovery. Plasma FFA and glycerol increased 8 h after standing in the colic horses. In both groups, plasma concentrations of CK increased and serum phosphate decreased post-anaesthesia. On Day 7 most parameters were not different between groups. Colic horses lost on average 8% of their initial weight. Eleven colic horses completed the study.</p> <p>Conclusion</p> <p>Colic horses entered anaesthesia with altered metabolism and in a negative oxygen balance. Muscle oxygenation was insufficient during anaesthesia in both groups, although to a lesser extent in the healthy horses. The post-anaesthetic period was associated with increased lipolysis and weight loss in the colic horses, indicating a negative energy balance during the first week post-operatively.</p

Gene-Expression Profiling Suggests Impaired Signaling via the Interferon Pathway in Cstb(-/-) Microglia

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1, OMIM254800) is an autosomal recessive neurodegenerative disorder characterized by stimulus-sensitive and action-activated myoclonus, tonic-clonic epileptic seizures, and ataxia. Loss-of-function mutations in the gene encoding the cysteine protease inhibitor cystatin B (CSTB) underlie EPM1. The deficiency of CSTB in mice (Cstb(-/-) mice) generates a phenotype resembling the symptoms of EPM1 patients and is accompanied by microglial activation at two weeks of age and an upregulation of immune system-associated genes in the cerebellum at one month of age. To shed light on molecular pathways and processes linked to CSTB deficiency in microglia we characterized the transcriptome of cultured Cstb(-/-) mouse microglia using microarray hybridization and RNA sequencing (RNA-seq). The gene expression profiles obtained with these two techniques were in good accordance and not polarized to either pro- or anti-inflammatory status. In Cstb(-/-) microglia, altogether 184 genes were differentially expressed. Of these, 33 genes were identified by both methods. Several interferon-regulated genes were weaker expressed in Cstb(-/-) microglia compared to control. This was confirmed by quantitative real-time PCR of the transcripts Irf7 and Stat1. Subsequently, we explored the biological context of CSTB deficiency in microglia more deeply by functional enrichment and canonical pathway analysis. This uncovered a potential role for CSTB in chemotaxis, antigen-presentation, and in immune-and defense response-associated processes by altering JAK-STAT pathway signaling. These data support and expand the previously suggested involvement of inflammatory processes to the disease pathogenesis of EPM1 and connect CSTB deficiency in microglia to altered expression of interferon-regulated genes.Peer reviewe

BACH2 immunodeficiency illustrates an association between super-enhancers and haploinsufficiency.

The transcriptional programs that guide lymphocyte differentiation depend on the precise expression and timing of transcription factors (TFs). The TF BACH2 is essential for T and B lymphocytes and is associated with an archetypal super-enhancer (SE). Single-nucleotide variants in the BACH2 locus are associated with several autoimmune diseases, but BACH2 mutations that cause Mendelian monogenic primary immunodeficiency have not previously been identified. Here we describe a syndrome of BACH2-related immunodeficiency and autoimmunity (BRIDA) that results from BACH2 haploinsufficiency. Affected subjects had lymphocyte-maturation defects that caused immunoglobulin deficiency and intestinal inflammation. The mutations disrupted protein stability by interfering with homodimerization or by causing aggregation. We observed analogous lymphocyte defects in Bach2-heterozygous mice. More generally, we observed that genes that cause monogenic haploinsufficient diseases were substantially enriched for TFs and SE architecture. These findings reveal a previously unrecognized feature of SE architecture in Mendelian diseases of immunity: heterozygous mutations in SE-regulated genes identified by whole-exome/genome sequencing may have greater significance than previously recognized