Widespread intron retention in mammals functionally tunes transcriptomes

© 2014 Braunschweig et al.; Published by Cold Spring Harbor Laboratory Press. This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date (see http://genome.cshlp.org/site/misc/terms.xhtml). After six months, it is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), as described at http://creativecommons.org/licenses/by-nc/4.0/.Alternative splicing (AS) of precursor RNAs is responsible for greatly expanding the regulatory and functional capacity of eukaryotic genomes. Of the different classes of AS, intron retention (IR) is the least well understood. In plants and unicellular eukaryotes, IR is the most common form of AS, whereas in animals, it is thought to represent the least prevalent form. Using high-coverage poly(A)(+) RNA-seq data, we observe that IR is surprisingly frequent in mammals, affecting transcripts from as many as three-quarters of multiexonic genes. A highly correlated set of cis features comprising an "IR code" reliably discriminates retained from constitutively spliced introns. We show that IR acts widely to reduce the levels of transcripts that are less or not required for the physiology of the cell or tissue type in which they are detected. This "transcriptome tuning" function of IR acts through both nonsense-mediated mRNA decay and nuclear sequestration and turnover of IR transcripts. We further show that IR is linked to a cross-talk mechanism involving localized stalling of RNA polymerase II (Pol II) and reduced availability of spliceosomal components. Collectively, the results implicate a global checkpoint-type mechanism whereby reduced recruitment of splicing components coupled to Pol II pausing underlies widespread IR-mediated suppression of inappropriately expressed transcripts.This work was supported by grants from the Canadian Institutes of Health Research and Canadian Cancer Society (B.J.B.); EMBO long-term fellowships (U.B. and T.G.-P.); Human Frontier Science Program Organization long-term fellowships (U.B. and M.I.); an OSCI fellowship (T.G.-P.); CIHR postdoctoral and Marie Curie IOF fellowships (N.L.B.-M.); and an NSERC studentship (E.N.).info:eu-repo/semantics/publishedVersio

Suv39h-Mediated Histone H3 Lysine 9 Methylation Directs DNA Methylation to Major Satellite Repeats at Pericentric Heterochromatin

AbstractBackground: Histone H3 lysine 9 (H3-K9) methylation and DNA methylation are characteristic hallmarks of mammalian heterochromatin. H3-K9 methylation was recently shown to be a prerequisite for DNA methylation in Neurospora crassa and Arabidopsis thaliana. Currently, it is unknown whether a similar dependence exists in mammalian organisms.Results: Here, we demonstrate a physical and functional link between the Suv39h-HP1 histone methylation system and DNA methyltransferase 3b (Dnmt3b) in mammals. Whereas in wild-type cells Dnmt3b interacts with HP1α and is concentrated at heterochromatic foci, it fails to localize to these regions in Suv39h double null (dn) mouse embryonic stem (ES) cells. Consistently, the Suv39h dn ES cells display an altered DNA methylation profile at pericentric satellite repeats, but not at other repeat sequences. In contrast, H3-K9 trimethylation at pericentric heterochromatin is not impaired in Dnmt1 single- or Dnmt3a/Dnmt3b double-deficient ES cells. We also show that pericentric heterochromatin is not transcriptionally inert and can give rise to transcripts spanning the major satellite repeats.Conclusions: These data demonstrate an evolutionarily conserved pathway between histone H3-K9 methylation and DNA methylation in mammals. While the Suv39h HMTases are required to direct H3-K9 trimethylation and Dnmt3b-dependent DNA methylation at pericentric repeats, DNA methylation at centromeric repeats occurs independent of Suv39h function. Thus, our data also indicate a more complex interrelatedness between histone and DNA methylation systems in mammals. Both methylation systems are likely to be important in reinforcing the stability of heterochromatic subdomains and thereby in protecting genome integrity

A slow transcription rate causes embryonic lethality and perturbs kinetic coupling of neuronal genes

The rate of RNA polymerase II (RNAPII) elongation has an important role in the control of alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked in for a slow elongating form of RNAPII. We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice. Focusing on neuronal differentiation as a model, we observed that slow elongation impairs development of the neural lineage from ESCs, which is accompanied by changes in AS and in gene expression along this pathway. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is greater in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.Fil: Maslon, M.. University of Edinburgh; Reino UnidoFil: Braunschweig, U.. University of Toronto; CanadáFil: Aitken, S.. University of Edinburgh; Reino UnidoFil: Mann, A.R.. University of Edinburgh; Reino UnidoFil: Kilanowski, F.. University of Edinburgh; Reino UnidoFil: Hunter, C.H.. University of Edinburgh; Reino UnidoFil: Blencowe, B.J.. University of Toronto; CanadáFil: Kornblihtt, Alberto Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; ArgentinaFil: Adams, I.. University of Edinburgh; Reino UnidoFil: Cáceres, J.F.. University of Edinburgh; Reino Unid

Searches at HERA for Squarks in R-Parity Violating Supersymmetry

A search for squarks in R-parity violating supersymmetry is performed in e^+p collisions at HERA at a centre of mass energy of 300 GeV, using H1 data corresponding to an integrated luminosity of 37 pb^(-1). The direct production of single squarks of any generation in positron-quark fusion via a Yukawa coupling lambda' is considered, taking into account R-parity violating and conserving decays of the squarks. No significant deviation from the Standard Model expectation is found. The results are interpreted in terms of constraints within the Minimal Supersymmetric Standard Model (MSSM), the constrained MSSM and the minimal Supergravity model, and their sensitivity to the model parameters is studied in detail. For a Yukawa coupling of electromagnetic strength, squark masses below 260 GeV are excluded at 95% confidence level in a large part of the parameter space. For a 100 times smaller coupling strength masses up to 182 GeV are excluded.Comment: 32 pages, 14 figures, 3 table

Multiplicity Structure of the Hadronic Final State in Diffractive Deep-Inelastic Scattering at HERA

The multiplicity structure of the hadronic system X produced in deep-inelastic processes at HERA of the type ep -> eXY, where Y is a hadronic system with mass M_Y< 1.6 GeV and where the squared momentum transfer at the pY vertex, t, is limited to |t|<1 GeV^2, is studied as a function of the invariant mass M_X of the system X. Results are presented on multiplicity distributions and multiplicity moments, rapidity spectra and forward-backward correlations in the centre-of-mass system of X. The data are compared to results in e+e- annihilation, fixed-target lepton-nucleon collisions, hadro-produced diffractive final states and to non-diffractive hadron-hadron collisions. The comparison suggests a production mechanism of virtual photon dissociation which involves a mixture of partonic states and a significant gluon content. The data are well described by a model, based on a QCD-Regge analysis of the diffractive structure function, which assumes a large hard gluonic component of the colourless exchange at low Q^2. A model with soft colour interactions is also successful.Comment: 22 pages, 4 figures, submitted to Eur. Phys. J., error in first submission - omitted bibliograph

Differential (2+1) Jet Event Rates and Determination of alpha_s in Deep Inelastic Scattering at HERA

Events with a (2+1) jet topology in deep-inelastic scattering at HERA are studied in the kinematic range 200 < Q^2< 10,000 GeV^2. The rate of (2+1) jet events has been determined with the modified JADE jet algorithm as a function of the jet resolution parameter and is compared with the predictions of Monte Carlo models. In addition, the event rate is corrected for both hadronization and detector effects and is compared with next-to-leading order QCD calculations. A value of the strong coupling constant of alpha_s(M_Z^2)= 0.118+- 0.002 (stat.)^(+0.007)_(-0.008) (syst.)^(+0.007)_(-0.006) (theory) is extracted. The systematic error includes uncertainties in the calorimeter energy calibration, in the description of the data by current Monte Carlo models, and in the knowledge of the parton densities. The theoretical error is dominated by the renormalization scale ambiguity.Comment: 25 pages, 6 figures, 3 tables, submitted to Eur. Phys.

Measurements of Transverse Energy Flow in Deep-Inelastic Scattering at HERA

Measurements of transverse energy flow are presented for neutral current deep-inelastic scattering events produced in positron-proton collisions at HERA. The kinematic range covers squared momentum transfers Q^2 from 3.2 to 2,200 GeV^2, the Bjorken scaling variable x from 8.10^{-5} to 0.11 and the hadronic mass W from 66 to 233 GeV. The transverse energy flow is measured in the hadronic centre of mass frame and is studied as a function of Q^2, x, W and pseudorapidity. A comparison is made with QCD based models. The behaviour of the mean transverse energy in the central pseudorapidity region and an interval corresponding to the photon fragmentation region are analysed as a function of Q^2 and W.Comment: 26 pages, 8 figures, submitted to Eur. Phys.

The Insulator Protein SU(HW) Fine-Tunes Nuclear Lamina Interactions of the Drosophila Genome

Specific interactions of the genome with the nuclear lamina (NL) are thought to assist chromosome folding inside the nucleus and to contribute to the regulation of gene expression. High-resolution mapping has recently identified hundreds of large, sharply defined lamina-associated domains (LADs) in the human genome, and suggested that the insulator protein CTCF may help to demarcate these domains. Here, we report the detailed structure of LADs in Drosophila cells, and investigate the putative roles of five insulator proteins in LAD organization. We found that the Drosophila genome is also organized in discrete LADs, which are about five times smaller than human LADs but contain on average a similar number of genes. Systematic comparison to new and published insulator binding maps shows that only SU(HW) binds preferentially at LAD borders and at specific positions inside LADs, while GAF, CTCF, BEAF-32 and DWG are mostly absent from these regions. By knockdown and overexpression studies we demonstrate that SU(HW) weakens genome – NL interactions through a local antagonistic effect, but we did not obtain evidence that it is essential for border formation. Our results provide insights into the evolution of LAD organization and identify SU(HW) as a fine-tuner of genome – NL interactions