Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment
Identifieur interne : 000790 ( Pmc/Corpus ); précédent : 000789; suivant : 000791Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment
Auteurs : Justin J. L. Wong ; Dadi Gao ; Trung V. Nguyen ; Chau-To Kwok ; Michelle Van Geldermalsen ; Rob Middleton ; Natalia Pinello ; Annora Thoeng ; Rajini Nagarajah ; Jeff Holst ; William Ritchie ; John E. J. RaskoSource :
- Nature Communications [ 2041-1723 ] ; 2017.
Abstract
While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.
Url:
DOI: 10.1038/ncomms15134
PubMed: 28480880
PubMed Central: 5424149
Links to Exploration step
PMC:5424149Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment</title><author><name sortKey="Wong, Justin J L" sort="Wong, Justin J L" uniqKey="Wong J" first="Justin J. L." last="Wong">Justin J. L. Wong</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Gao, Dadi" sort="Gao, Dadi" uniqKey="Gao D" first="Dadi" last="Gao">Dadi Gao</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Nguyen, Trung V" sort="Nguyen, Trung V" uniqKey="Nguyen T" first="Trung V." last="Nguyen">Trung V. Nguyen</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Kwok, Chau To" sort="Kwok, Chau To" uniqKey="Kwok C" first="Chau-To" last="Kwok">Chau-To Kwok</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Van Geldermalsen, Michelle" sort="Van Geldermalsen, Michelle" uniqKey="Van Geldermalsen M" first="Michelle" last="Van Geldermalsen">Michelle Van Geldermalsen</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a5"><institution>Origins of Cancer Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Middleton, Rob" sort="Middleton, Rob" uniqKey="Middleton R" first="Rob" last="Middleton">Rob Middleton</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Pinello, Natalia" sort="Pinello, Natalia" uniqKey="Pinello N" first="Natalia" last="Pinello">Natalia Pinello</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Thoeng, Annora" sort="Thoeng, Annora" uniqKey="Thoeng A" first="Annora" last="Thoeng">Annora Thoeng</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Nagarajah, Rajini" sort="Nagarajah, Rajini" uniqKey="Nagarajah R" first="Rajini" last="Nagarajah">Rajini Nagarajah</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Holst, Jeff" sort="Holst, Jeff" uniqKey="Holst J" first="Jeff" last="Holst">Jeff Holst</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a5"><institution>Origins of Cancer Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Ritchie, William" sort="Ritchie, William" uniqKey="Ritchie W" first="William" last="Ritchie">William Ritchie</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a6"><institution>CNRS, UMR 5203</institution>, 34094 Montpellier,<country>France</country></nlm:aff></affiliation></author><author><name sortKey="Rasko, John E J" sort="Rasko, John E J" uniqKey="Rasko J" first="John E. J." last="Rasko">John E. J. Rasko</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a7"><institution>Cell and Molecular Therapies, Royal Prince Alfred Hospital</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">28480880</idno><idno type="pmc">5424149</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5424149</idno><idno type="RBID">PMC:5424149</idno><idno type="doi">10.1038/ncomms15134</idno><date when="2017">2017</date><idno type="wicri:Area/Pmc/Corpus">000790</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000790</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment</title><author><name sortKey="Wong, Justin J L" sort="Wong, Justin J L" uniqKey="Wong J" first="Justin J. L." last="Wong">Justin J. L. Wong</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Gao, Dadi" sort="Gao, Dadi" uniqKey="Gao D" first="Dadi" last="Gao">Dadi Gao</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Nguyen, Trung V" sort="Nguyen, Trung V" uniqKey="Nguyen T" first="Trung V." last="Nguyen">Trung V. Nguyen</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Kwok, Chau To" sort="Kwok, Chau To" uniqKey="Kwok C" first="Chau-To" last="Kwok">Chau-To Kwok</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Van Geldermalsen, Michelle" sort="Van Geldermalsen, Michelle" uniqKey="Van Geldermalsen M" first="Michelle" last="Van Geldermalsen">Michelle Van Geldermalsen</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a5"><institution>Origins of Cancer Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Middleton, Rob" sort="Middleton, Rob" uniqKey="Middleton R" first="Rob" last="Middleton">Rob Middleton</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Pinello, Natalia" sort="Pinello, Natalia" uniqKey="Pinello N" first="Natalia" last="Pinello">Natalia Pinello</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a2"><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Thoeng, Annora" sort="Thoeng, Annora" uniqKey="Thoeng A" first="Annora" last="Thoeng">Annora Thoeng</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Nagarajah, Rajini" sort="Nagarajah, Rajini" uniqKey="Nagarajah R" first="Rajini" last="Nagarajah">Rajini Nagarajah</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Holst, Jeff" sort="Holst, Jeff" uniqKey="Holst J" first="Jeff" last="Holst">Jeff Holst</name><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a5"><institution>Origins of Cancer Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author><author><name sortKey="Ritchie, William" sort="Ritchie, William" uniqKey="Ritchie W" first="William" last="Ritchie">William Ritchie</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a4"><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a6"><institution>CNRS, UMR 5203</institution>, 34094 Montpellier,<country>France</country></nlm:aff></affiliation></author><author><name sortKey="Rasko, John E J" sort="Rasko, John E J" uniqKey="Rasko J" first="John E. J." last="Rasko">John E. J. Rasko</name><affiliation><nlm:aff id="a1"><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a3"><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></nlm:aff></affiliation><affiliation><nlm:aff id="a7"><institution>Cell and Molecular Therapies, Royal Prince Alfred Hospital</institution>, Camperdown, New South Wales 2050,<country>Australia</country></nlm:aff></affiliation></author></analytic><series><title level="j">Nature Communications</title><idno type="eISSN">2041-1723</idno><imprint><date when="2017">2017</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Yap, K" uniqKey="Yap K">K. Yap</name></author><author><name sortKey="Lim, Z Q" uniqKey="Lim Z">Z. Q. Lim</name></author><author><name sortKey="Khandelia, P" uniqKey="Khandelia P">P. Khandelia</name></author><author><name sortKey="Friedman, B" uniqKey="Friedman B">B. Friedman</name></author><author><name sortKey="Makeyev, E V" uniqKey="Makeyev E">E. V. Makeyev</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wong J J L" uniqKey="Wong J J L">Wong J.J.-L.</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Moore, M J" uniqKey="Moore M">M. J. Moore</name></author><author><name sortKey="Wang, Q" uniqKey="Wang Q">Q. Wang</name></author><author><name sortKey="Kennedy, C J" uniqKey="Kennedy C">C. J. Kennedy</name></author><author><name sortKey="Silver, P A" uniqKey="Silver P">P. A. Silver</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ergun, A" uniqKey="Ergun A">A. Ergun</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shaked, I" uniqKey="Shaked I">I. Shaked</name></author><author><name sortKey="Zimmerman, G" uniqKey="Zimmerman G">G. Zimmerman</name></author><author><name sortKey="Soreq, H" uniqKey="Soreq H">H. Soreq</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Padgett, R A" uniqKey="Padgett R">R. A. Padgett</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Edwards, C" uniqKey="Edwards C">C. Edwards</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pimentel, H" uniqKey="Pimentel H">H. Pimentel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Braunschweig, U" uniqKey="Braunschweig U">U. Braunschweig</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guo, R" uniqKey="Guo R">R. Guo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De La Mata, M" uniqKey="De La Mata M">M. de la Mata</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maunakea, A K" uniqKey="Maunakea A">A. K. Maunakea</name></author><author><name sortKey="Chepelev, I" uniqKey="Chepelev I">I. Chepelev</name></author><author><name sortKey="Cui, K" uniqKey="Cui K">K. Cui</name></author><author><name sortKey="Zhao, K" uniqKey="Zhao K">K. Zhao</name></author><author><name sortKey="Intragenic, D N A" uniqKey="Intragenic D">D. N. A. Intragenic</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shukla, S" uniqKey="Shukla S">S. Shukla</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gascard, P" uniqKey="Gascard P">P. Gascard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lin, S" uniqKey="Lin S">S. Lin</name></author><author><name sortKey="Coutinho Mansfield, G" uniqKey="Coutinho Mansfield G">G. Coutinho-Mansfield</name></author><author><name sortKey="Wang, D" uniqKey="Wang D">D. Wang</name></author><author><name sortKey="Pandit, S" uniqKey="Pandit S">S. Pandit</name></author><author><name sortKey="Fu, X D" uniqKey="Fu X">X.-D. Fu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fong, Y W" uniqKey="Fong Y">Y. W. Fong</name></author><author><name sortKey="Zhou, Q" uniqKey="Zhou Q">Q. Zhou</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lev Maor, G" uniqKey="Lev Maor G">G. Lev Maor</name></author><author><name sortKey="Yearim, A" uniqKey="Yearim A">A. Yearim</name></author><author><name sortKey="Ast, G" uniqKey="Ast G">G. Ast</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yearim, A" uniqKey="Yearim A">A. Yearim</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Olins, A L" uniqKey="Olins A">A. L. Olins</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Porfiri, E" uniqKey="Porfiri E">E. Porfiri</name></author><author><name sortKey="Hoffbrand, A" uniqKey="Hoffbrand A">A. Hoffbrand</name></author><author><name sortKey="Wickremasinghe, R" uniqKey="Wickremasinghe R">R. Wickremasinghe</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lodhi, I J" uniqKey="Lodhi I">I. J. Lodhi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Crawford, L J A" uniqKey="Crawford L">L. J. A. Crawford</name></author><author><name sortKey="Peake, R" uniqKey="Peake R">R. Peake</name></author><author><name sortKey="Price, S" uniqKey="Price S">S. Price</name></author><author><name sortKey="Morris, T C M" uniqKey="Morris T">T. C. M. Morris</name></author><author><name sortKey="Irvine, A E" uniqKey="Irvine A">A. E. Irvine</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jones, P A" uniqKey="Jones P">P. A. Jones</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li Byarlay, H" uniqKey="Li Byarlay H">H. Li-Byarlay</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Laurent, L" uniqKey="Laurent L">L. Laurent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gombart, A F" uniqKey="Gombart A">A. F. Gombart</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nacken, W" uniqKey="Nacken W">W. Nacken</name></author><author><name sortKey="Roth, J" uniqKey="Roth J">J. Roth</name></author><author><name sortKey="Sorg, C" uniqKey="Sorg C">C. Sorg</name></author><author><name sortKey="Kerkhoff, C" uniqKey="Kerkhoff C">C. Kerkhoff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nathan, C" uniqKey="Nathan C">C. Nathan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chan, K T" uniqKey="Chan K">K. T. Chan</name></author><author><name sortKey="Creed, S J" uniqKey="Creed S">S. J. Creed</name></author><author><name sortKey="Bear, J E" uniqKey="Bear J">J. E. Bear</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Salerno, D M" uniqKey="Salerno D">D. M. Salerno</name></author><author><name sortKey="Tront, J S" uniqKey="Tront J">J. S. Tront</name></author><author><name sortKey="Hoffman, B" uniqKey="Hoffman B">B. Hoffman</name></author><author><name sortKey="Liebermann, D A" uniqKey="Liebermann D">D. A. Liebermann</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gombart, A F" uniqKey="Gombart A">A. F. Gombart</name></author><author><name sortKey="Grewal, J" uniqKey="Grewal J">J. Grewal</name></author><author><name sortKey="Koeffler, H P" uniqKey="Koeffler H">H. P. Koeffler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baubec, T" uniqKey="Baubec T">T. Baubec</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baubec, T" uniqKey="Baubec T">T. Baubec</name></author><author><name sortKey="Ivanek, R" uniqKey="Ivanek R">R. Ivánek</name></author><author><name sortKey="Lienert, F" uniqKey="Lienert F">F. Lienert</name></author><author><name sortKey="Schubeler, D" uniqKey="Schubeler D">D. Schübeler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Young, J I" uniqKey="Young J">J. I. Young</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maxwell, S S" uniqKey="Maxwell S">S. S. Maxwell</name></author><author><name sortKey="Pelka, G J" uniqKey="Pelka G">G. J. Pelka</name></author><author><name sortKey="Tam, P P L" uniqKey="Tam P">P. P. L. Tam</name></author><author><name sortKey="El Osta, A" uniqKey="El Osta A">A. El-Osta</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Long, S W" uniqKey="Long S">S. W. Long</name></author><author><name sortKey="Ooi, J Y Y" uniqKey="Ooi J">J. Y. Y. Ooi</name></author><author><name sortKey="Yau, P M" uniqKey="Yau P">P. M. Yau</name></author><author><name sortKey="Jones, P L" uniqKey="Jones P">P. L. Jones</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ghosh, R P" uniqKey="Ghosh R">R. P. Ghosh</name></author><author><name sortKey="Horowitz Scherer, R A" uniqKey="Horowitz Scherer R">R. A. Horowitz-Scherer</name></author><author><name sortKey="Nikitina, T" uniqKey="Nikitina T">T. Nikitina</name></author><author><name sortKey="Shlyakhtenko, L S" uniqKey="Shlyakhtenko L">L. S. Shlyakhtenko</name></author><author><name sortKey="Woodcock, C L" uniqKey="Woodcock C">C. L. Woodcock</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shepard, P J" uniqKey="Shepard P">P. J. Shepard</name></author><author><name sortKey="Hertel, K J" uniqKey="Hertel K">K. J. Hertel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dichmann, D S" uniqKey="Dichmann D">D. S. Dichmann</name></author><author><name sortKey="Walentek, P" uniqKey="Walentek P">P. Walentek</name></author><author><name sortKey="Harland, R M" uniqKey="Harland R">R. M. Harland</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Best, A" uniqKey="Best A">A. Best</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hussein, S M I" uniqKey="Hussein S">S. M. I. Hussein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dvinge, H" uniqKey="Dvinge H">H. Dvinge</name></author><author><name sortKey="Bradley, R" uniqKey="Bradley R">R. Bradley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jung, H" uniqKey="Jung H">H. Jung</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ni, T" uniqKey="Ni T">T. Ni</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wen, H" uniqKey="Wen H">H. Wen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Marina, R J" uniqKey="Marina R">R. J. Marina</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tekpli, X" uniqKey="Tekpli X">X. Tekpli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Song, C X" uniqKey="Song C">C.-X. Song</name></author><author><name sortKey="Yi, C" uniqKey="Yi C">C. Yi</name></author><author><name sortKey="He, C" uniqKey="He C">C. He</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, X" uniqKey="Yang X">X. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Buckstein, R" uniqKey="Buckstein R">R. Buckstein</name></author><author><name sortKey="Yee, K" uniqKey="Yee K">K. Yee</name></author><author><name sortKey="Wells, R A" uniqKey="Wells R">R. A. Wells</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brown, R" uniqKey="Brown R">R. Brown</name></author><author><name sortKey="Curry, E" uniqKey="Curry E">E. Curry</name></author><author><name sortKey="Magnani, L" uniqKey="Magnani L">L. Magnani</name></author><author><name sortKey="Wilhelm Benartzi, C S" uniqKey="Wilhelm Benartzi C">C. S. Wilhelm-Benartzi</name></author><author><name sortKey="Borley, J" uniqKey="Borley J">J. Borley</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Middleton, R" uniqKey="Middleton R">R. Middleton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wong, J" uniqKey="Wong J">J. Wong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mohammed, H" uniqKey="Mohammed H">H. Mohammed</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Challen, G A" uniqKey="Challen G">G. A. Challen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jeong, M" uniqKey="Jeong M">M. Jeong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lund, K" uniqKey="Lund K">K. Lund</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Trapnell, C" uniqKey="Trapnell C">C. Trapnell</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y. Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dobin, A" uniqKey="Dobin A">A. Dobin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Audic, S" uniqKey="Audic S">S. Audic</name></author><author><name sortKey="Claverie, J M" uniqKey="Claverie J">J. M. Claverie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Huang, D W" uniqKey="Huang D">D. W. Huang</name></author><author><name sortKey="Sherman, B T" uniqKey="Sherman B">B. T. Sherman</name></author><author><name sortKey="Lempicki, R A" uniqKey="Lempicki R">R. A. Lempicki</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Nat Commun</journal-id><journal-id journal-id-type="iso-abbrev">Nat Commun</journal-id><journal-title-group><journal-title>Nature Communications</journal-title></journal-title-group><issn pub-type="epub">2041-1723</issn><publisher><publisher-name>Nature Publishing Group</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">28480880</article-id><article-id pub-id-type="pmc">5424149</article-id><article-id pub-id-type="pii">ncomms15134</article-id><article-id pub-id-type="doi">10.1038/ncomms15134</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Wong</surname><given-names>Justin J. -L.</given-names></name><xref ref-type="corresp" rid="c1">a</xref><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref><xref ref-type="aff" rid="a3">3</xref><xref ref-type="author-notes" rid="n1">*</xref></contrib><contrib contrib-type="author"><name><surname>Gao</surname><given-names>Dadi</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a4">4</xref><xref ref-type="author-notes" rid="n1">*</xref></contrib><contrib contrib-type="author"><name><surname>Nguyen</surname><given-names>Trung V.</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref><xref ref-type="aff" rid="a3">3</xref><xref ref-type="author-notes" rid="n1">*</xref></contrib><contrib contrib-type="author"><name><surname>Kwok</surname><given-names>Chau-To</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref><xref ref-type="aff" rid="a3">3</xref></contrib><contrib contrib-type="author"><name><surname>van Geldermalsen</surname><given-names>Michelle</given-names></name><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a5">5</xref></contrib><contrib contrib-type="author"><name><surname>Middleton</surname><given-names>Rob</given-names></name><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a4">4</xref></contrib><contrib contrib-type="author"><name><surname>Pinello</surname><given-names>Natalia</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a2">2</xref><xref ref-type="aff" rid="a3">3</xref></contrib><contrib contrib-type="author"><name><surname>Thoeng</surname><given-names>Annora</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a3">3</xref></contrib><contrib contrib-type="author"><name><surname>Nagarajah</surname><given-names>Rajini</given-names></name><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a3">3</xref></contrib><contrib contrib-type="author"><name><surname>Holst</surname><given-names>Jeff</given-names></name><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a5">5</xref><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-0377-9318</contrib-id></contrib><contrib contrib-type="author"><name><surname>Ritchie</surname><given-names>William</given-names></name><xref ref-type="corresp" rid="c2">b</xref><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a4">4</xref><xref ref-type="aff" rid="a6">6</xref></contrib><contrib contrib-type="author"><name><surname>Rasko</surname><given-names>John E. J.</given-names></name><xref ref-type="corresp" rid="c3">c</xref><xref ref-type="aff" rid="a1">1</xref><xref ref-type="aff" rid="a3">3</xref><xref ref-type="aff" rid="a7">7</xref></contrib><aff id="a1"><label>1</label><institution>Gene & Stem Cell Therapy Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></aff><aff id="a2"><label>2</label><institution>Gene Regulation in Cancer Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></aff><aff id="a3"><label>3</label><institution>Sydney Medical School, University of Sydney</institution>, Sydney, New South Wales 2006,<country>Australia</country></aff><aff id="a4"><label>4</label><institution>Bioinformatics Laboratory, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></aff><aff id="a5"><label>5</label><institution>Origins of Cancer Program, Centenary Institute, University of Sydney</institution>, Camperdown, New South Wales 2050,<country>Australia</country></aff><aff id="a6"><label>6</label><institution>CNRS, UMR 5203</institution>, 34094 Montpellier,<country>France</country></aff><aff id="a7"><label>7</label><institution>Cell and Molecular Therapies, Royal Prince Alfred Hospital</institution>, Camperdown, New South Wales 2050,<country>Australia</country></aff></contrib-group><author-notes><corresp id="c1"><label>a</label><email>j.wong@centenary.org.au</email></corresp><corresp id="c2"><label>b</label><email>william.ritchie@igf.cnrs.fr</email></corresp><corresp id="c3"><label>c</label><email>j.rasko@centenary.org.au</email></corresp><fn id="n1"><label>*</label><p>These authors contributed equally to this work.</p></fn></author-notes><pub-date pub-type="epub"><day>08</day><month>05</month><year>2017</year></pub-date><pub-date pub-type="collection"><year>2017</year></pub-date><volume>8</volume><elocation-id>15134</elocation-id><history><date date-type="received"><day>06</day><month>09</month><year>2016</year></date><date date-type="accepted"><day>02</day><month>03</month><year>2017</year></date></history><permissions><copyright-statement>Copyright © 2017, The Author(s)</copyright-statement><copyright-year>2017</copyright-year><copyright-holder>The Author(s)</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/4.0/"><pmc-comment>author-paid</pmc-comment>
<license-p>This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><abstract><p>While intron retention (IR) is considered a widely conserved and distinct mechanism of gene expression control, its regulation is poorly understood. Here we show that DNA methylation directly regulates IR. We also find reduced occupancy of MeCP2 near the splice junctions of retained introns, mirroring the reduced DNA methylation at these sites. Accordingly, MeCP2 depletion in tissues and cells enhances IR. By analysing the MeCP2 interactome using mass spectrometry and RNA co-precipitation, we demonstrate that decreased MeCP2 binding near splice junctions facilitates IR via reduced recruitment of splicing factors, including Tra2b, and increased RNA polymerase II stalling. These results suggest an association between IR and a slower rate of transcription elongation, which reflects inefficient splicing factor recruitment. In summary, our results reinforce the interdependency between alternative splicing involving IR and epigenetic controls of gene expression.</p></abstract><abstract abstract-type="web-summary"><p>Intron retention is a conserved mechanism that controls gene expression but its regulation is poorly understood. Here, the authors provide evidence that DNA methylation regulates intron retention and find reduced MeCP2 occupancy and splicing factor recruitment near affected splice junctions.</p></abstract></article-meta></front><body><p>Alternative splicing (AS) is essential for diverse physiological functions including cell division, differentiation, stress response and lineage specification<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b3">3</xref><xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref>. About one-fifth of human diseases are caused by disrupted splicing<xref ref-type="bibr" rid="b6">6</xref>, making normal AS a vital regulatory mechanism in humans. Accordingly, the search for mechanisms underlying AS has been one of the most active research areas in molecular and cellular biology.</p><p>One mode of AS termed intron retention (IR) occurs when an intron is not excised by the spliceosome, and is preserved in the final mature messenger RNA (mRNA). Intron-retaining transcripts typically contain premature termination codons, resulting in their removal by surveillance mechanisms including nonsense-mediated decay (NMD). Whilst AS generates proteomic complexity from most genes, IR is a distinctive mechanism that reduces gene expression, which in turn regulates cellular differentiation<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b7">7</xref><xref ref-type="bibr" rid="b8">8</xref>. Following our initial report that numerous introns are constitutively retained to induce mRNA degradation via NMD during normal granulopoiesis<xref ref-type="bibr" rid="b2">2</xref>, the prevalence of IR in ∼40 cell and tissue types, and its conservation across 11 vertebrate species has been established<xref ref-type="bibr" rid="b9">9</xref>. In general, the molecular mechanisms that regulate IR have not been thoroughly examined. However, in one well-characterized example, the H3.3K36me3 reader, BS69, binds to a U5 associated splicing factor, EFTUD2, to antagonize splicing activity, leading to increased IR (ref. <xref ref-type="bibr" rid="b10">10</xref>). Of all IR events identified in that study, only one-fifth were BS69-dependent<xref ref-type="bibr" rid="b10">10</xref>. Thus other as yet undefined mechanisms exist to regulate IR.</p><p>Much is known about the mechanisms that regulate other modes of AS including exon skipping and inclusion, but their relevance to the control of IR remains uncertain. RNA Poll II stalling, which improves splice signal recognition is integral to the ‘kinetic' model of AS that involves epigenetic changes<xref ref-type="bibr" rid="b11">11</xref>. For example, altered DNA methylation levels and downstream binding of transcription factors as well as regulatory proteins such as CTCF and MeCP2 induce RNA Pol II stalling to promote alternative exon skipping and inclusion<xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b13">13</xref>. RNA Pol II stalling has been observed near splice junctions of retained introns but it is unclear if it directly regulates IR or whether other cis- and/or trans- regulatory factors initiate the process<xref ref-type="bibr" rid="b9">9</xref>. Also, lower DNA methylation levels near splice junctions have been associated with increased IR in normal human breast tissues<xref ref-type="bibr" rid="b14">14</xref>. This latter finding raises the intriguing possibility that DNA methylation might regulate IR by altering RNA Pol II processivity and splice site recognition according to the ‘kinetic' model.</p><p>Alternatively, RNA Pol II stalling near splice junctions of retained introns could occur consequent to inefficient splicing factor recruitment<xref ref-type="bibr" rid="b15">15</xref><xref ref-type="bibr" rid="b16">16</xref>. This ‘recruitment' model of AS has been shown to function in the case of exon skipping, which is also associated with DNA methylation changes<xref ref-type="bibr" rid="b17">17</xref><xref ref-type="bibr" rid="b18">18</xref>. Taken together changes in DNA methylation and RNA Pol II elongation are both factors that are known to be associated with IR. However, it is currently unknown whether DNA methylation directly controls IR in the same manner by which it controls exon skipping and inclusion. Any functional connection between DNA methylation and the ‘kinetic' and ‘recruitment' models in the regulation of IR has been an open question to date.</p><p>In this study, we sought to determine whether IR is regulated by spliceosomal and DNA methylation changes and how they are linked. We show that lower levels of DNA methylation near 5′ and 3′ splice junctions are significantly associated with increased IR in a myriad of primary cells and cell lines. Importantly, both chemical and genetic inhibition of DNA methylation promoted IR in mouse embryonic stem cells, mouse B cells, and in both mouse and human cell lines. We show reduced MeCP2 binding near 5′ and 3′ splice junctions of retained introns in mouse granulocytes and promyelocytes, which corresponds with the physiological decrease in DNA methylation levels. Further, knockdown or knockout of MeCP2 in primary cells and cell lines led to increased IR. To confirm that MeCP2 occupancy influences the efficiency of splicing, we show that reduced MeCP2 occupancy is associated with RNA Pol II stalling. Using immunoprecipitation followed by mass spectrometry (MS), we further show that many splicing factors are associated with MeCP2 including the Srsf family members and Tra2b. Consistent with inefficient splicing factor recruitment as reflected by RNA Pol II stalling, we demonstrate reduced binding of the splicing factor Tra2b to RNA transcripts containing retained introns along with reduced binding of MeCP2 to these transcripts. Our data link the two models of splicing by demonstrating a direct causative role of decreased DNA methylation in promoting IR via reduced MeCP2 occupancy, which in turn perturbs the recruitment of splicing factors such as Tra2b, leading to RNA Pol II stalling.</p><sec disp-level="1"><title>Results</title><sec disp-level="2"><title>IR correlates with high CpG density and low gene expression</title><p>We performed fluorescence-activated cell sorting (FACS) of primary mouse promyelocytes and granulocytes isolated from mouse bone marrow according to our established protocol (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1a,b</xref>)<xref ref-type="bibr" rid="b2">2</xref>. Our previous conservative algorithm to detect IR identified 121 introns from 86 transcripts that show differential IR in granulocytes compared to promyelocytes<xref ref-type="bibr" rid="b2">2</xref>. We have now refined our previous IRFinder algorithm to better estimate IR from sequencing reads and to include more introns in our analysis. This is achieved mainly by removing biases such as low mappability and overlapping features but also using split reads as a measure of correct splicing surrounding IR (Methods section). This method allows for a more sensitive detection of IR and of differential IR. Reanalysis of our data estimates 737 differentially retained introns (<italic>P</italic><0.05, Audic and Claverie test) between the two cell types out of a total of 3,200 and 3,888 introns retained (IR level ≥10%) in promyelocytes and granulocytes, respectively. Of these 737 differentially retained introns, substantially more introns (687/737, 93.22%) showed higher levels of retention in granulocytes (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1c</xref>); consistent with our previous report<xref ref-type="bibr" rid="b2">2</xref>. Gene ontology analysis (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1d</xref>) again reinforced significant representation of features associated with granulopoiesis and granulocyte maintenance, including nuclear lamina<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b19">19</xref> as well as metabolic and catabolic processes involving different types of lipids<xref ref-type="bibr" rid="b20">20</xref><xref ref-type="bibr" rid="b21">21</xref><xref ref-type="bibr" rid="b22">22</xref>.</p><p>Retained introns were significantly shorter in nucleotide length than non-retained introns (<italic>P</italic><1.00e-16, <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1e</xref>). Retained introns were also significantly associated with higher intronic GC content and CpG density (<italic>P</italic><0.05, <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1f,g</xref>), and higher CpG density near their flanking exon–intron boundaries (<italic>P</italic><0.05, <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1g</xref>).</p><p>Consistent with previous reports<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b9">9</xref> we also found that higher IR levels were associated with lower expression of coding mRNA in promyelocytes and granulocytes (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1h,i</xref>). These results indicate downregulation of gene expression via IR-triggered NMD or sequestering of intron-detained transcripts in the nucleus as previously reported<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b8">8</xref><xref ref-type="bibr" rid="b9">9</xref>.</p><p>Notably, the average distance between two genes with a retained intron was significantly shorter than the average distance between two randomly selected genes (<italic>P</italic><1.00e-16, <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 2a</xref>). Moreover, introns are frequently retained in the same gene and these are often adjacent to each other (<italic>P</italic><1.00e-16, <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 2b</xref>). This result not only reflects that multiple tandem introns are being retained in a given transcript but that IR may concomitantly affect many RNA transcripts encoded by neighbouring genes within the same chromosomal region.</p></sec><sec disp-level="2"><title>Reduced DNA methylation at splice junctions marks IR</title><p>DNA methylation occurs predominantly at CpG dinucleotides in mammals<xref ref-type="bibr" rid="b23">23</xref>. Changes in DNA methylation levels have established roles in the regulation of other modes of splicing including exon-inclusion and skipping<xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b13">13</xref><xref ref-type="bibr" rid="b18">18</xref><xref ref-type="bibr" rid="b24">24</xref>. Therefore, the strong association between increased CpG density and IR (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1g</xref>) prompted us to investigate whether DNA methylation also regulates IR. In order to determine the association between DNA methylation and IR, we performed whole-genome bisulfite sequencing (WGBS) on genomic DNA from promyelocytes and granulocytes. A total of 560 million and 633 million paired-end reads were aligned to the mouse genome (mm10) for promyelocytes and granulocytes, respectively. To measure the methylation level of each CpG in the genome, the beta value was calculated for each CpG (Methods section). Sharp transitions in DNA methylation patterns are known to mark splice junctions, suggesting that these signals may be crucial for the recognition of exons and introns<xref ref-type="bibr" rid="b25">25</xref>. We measured the levels of DNA methylation within ±100 bp of 5′ or 3′ splice junctions of all introns in promyelocytes and granulocytes. Significantly lower DNA methylation levels exist near both 5′ and 3′ junctions of retained introns compared to non-retained introns (<italic>P</italic>=2.06e-15 and 7.61e-12, <italic>t</italic>-test, respectively; <xref ref-type="fig" rid="f1">Fig. 1a</xref>). This result indicates that reduced DNA methylation marks the junctions of introns that are more likely to be retained (<xref ref-type="fig" rid="f1">Fig. 1a</xref>). Significantly lower DNA methylation levels also occurred within a 200-bp region spanning the middle of retained introns (<italic>P</italic>=6.77e-12, <italic>t</italic>-test, <xref ref-type="fig" rid="f1">Fig. 1a</xref>).</p><p>Of 737 differentially retained introns in granulocytes and promyelocytes, 125 (16.9%) showed reduced DNA methylation levels by ≥10% near splice junctions or within introns when IR levels increased in either cell type. This result indicates that DNA methylation-associated IR events may occur at a similar frequency as BS69-associated IR events<xref ref-type="bibr" rid="b10">10</xref>. Importantly, IR-associated DNA methylation changes affect numerous genes, including those relevant to granulopoiesis and granulocyte functions: <italic>Lmnb1</italic> (refs <xref ref-type="bibr" rid="b2">2</xref>, <xref ref-type="bibr" rid="b19">19</xref>), <italic>Ngp</italic> (ref. <xref ref-type="bibr" rid="b26">26</xref>), <italic>S100a8</italic> (ref. <xref ref-type="bibr" rid="b27">27</xref>), <italic>Mmp8</italic> (ref. <xref ref-type="bibr" rid="b28">28</xref>), <italic>Coro1a</italic> (ref. <xref ref-type="bibr" rid="b29">29</xref>), <italic>S100a9</italic> (ref. <xref ref-type="bibr" rid="b27">27</xref>), <italic>Gadd45a</italic> (ref. <xref ref-type="bibr" rid="b30">30</xref>) and <italic>Atf4</italic> (ref. <xref ref-type="bibr" rid="b31">31</xref>; <xref ref-type="fig" rid="f1">Fig. 1b,c</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 3</xref>). Clonal bisulfite sequencing confirmed lower DNA methylation levels are associated with increased IR levels for <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3 and <italic>S100a8</italic> intron 2 in granulocytes (<xref ref-type="fig" rid="f1">Fig. 1d</xref>).</p></sec><sec disp-level="2"><title>Reduced DNA methylation associated with IR is widespread</title><p>We next set out to determine the extent of IR events that are marked by reduced DNA methylation near splice junctions. We analysed publicly available matched RNA-seq and methylome data from normal and cancer cells, including mouse primary and reprogrammed fibroblasts, human embryonic stem cell lines (H1 and H9), a human lung fibroblast cell line (IMR90), a human colon cancer cell line (HCT116) and human neuron progenitor cells (see <xref ref-type="supplementary-material" rid="S1">Supplementary Table 1</xref> for details). Consistent with our findings in myeloid cells, significantly lower DNA methylation levels were observed near 3′ splice junctions flanking retained compared to non-retained introns in all cell types examined (<italic>P</italic><9.22e-8 in all cases, <italic>t</italic>-test, <xref ref-type="fig" rid="f2">Fig. 2a</xref>). Lower DNA methylation levels near 5′ splice junctions were significantly associated with retained introns in the H1 cell line and human neuron progenitors, but not in other cells (<xref ref-type="fig" rid="f2">Fig. 2a</xref>). Lower levels of DNA methylation were also present within the body of retained introns compared to non-retained introns in all cell types (<italic>P</italic><1.33e-5 in all cases, paired <italic>t</italic>-test, <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4a</xref>). Overall, lower DNA methylation levels near 3′ splice junctions and within introns consistently distinguished retained from non-retained introns in all cell types examined, demonstrating conservation of this phenomenon in human and mouse biology.</p></sec><sec disp-level="2"><title>Inhibition of DNA methylation increases intron retention</title><p>To determine whether DNA methylation changes can directly alter IR levels, we reanalysed RNA-seq data of cells with reduced DNA methylation, achieved via 5-Aza-2′-deoxycytidine (5-Aza) treatment or <italic>Dnmt</italic>-knockout (<italic>Dnmt</italic>-KO; see <xref ref-type="supplementary-material" rid="S1">Supplementary Table 1</xref> for details). <italic>DNMT1</italic> is the maintenance DNA methyltransferase while both <italic>DNMT3A</italic> and <italic>3B</italic> direct <italic>de novo</italic> DNA methylation within coding regions, and their binding patterns exhibit a high degree of overlap<xref ref-type="bibr" rid="b32">32</xref>. We identified 6,446, 1,844 and 20 differentially retained introns in <italic>Dnmt3a</italic> and <italic>3b</italic>-double-KO mouse B cells, <italic>Dnmt</italic>3a and <italic>3b-</italic>double-KO haemopoietic stem cells (HSCs) and 5-Aza-treated leukaemic OCI-AML3 cells compared to respective wild type and untreated cells. The number of introns exhibiting increased IR levels (<italic>Dnmt3a</italic> and <italic>3b</italic>-double-KO mouse B cells: 6418/6446, 99.57%; <italic>Dnmt3a</italic> and <italic>3b</italic>- double-KO HSCs: 977/1844, 52.98%; 5-Aza-treated OCI-AML3: 15/20, 75%) was significantly more than introns with decreased IR levels in all three cell types (<italic>P</italic><0.05, binomial test, <xref ref-type="fig" rid="f2">Fig. 2b</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4b</xref>). These results indicate that inhibition of DNA methylation is more likely to promote than inhibit IR. Notably, a lower percentage of introns with increased IR levels was observed in <italic>DNMT3B</italic> and hypomorphic <italic>DNMT1</italic> double-KO HCT116 cells (224/613, 37%) compared to control (<italic>P</italic><0.05 binomial test) (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4b</xref>). This discrepancy may be due either to incomplete <italic>DNMT1</italic> loss or the compensatory presence of <italic>DNMT3A</italic>.</p><p>We also treated the mouse promyelocyte cell line, MPRO, with 5-Aza, to determine the changes in the levels of nine intron-retaining transcripts following DNA demethylation. These intron-retaining transcripts were selected from those showing increased IR and concomitantly reduced DNA methylation levels in granulocytes compared to promyelocytes (<xref ref-type="fig" rid="f1">Fig. 1b–d</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 3</xref>). Since IR typically leads to NMD, we also inhibited NMD using caffeine treatment to measure the changes in IR when RNA degradation is limited. In three IR-affected genes, decreased DNA methylation levels near 5′ or 3′ splice sites in the DNA were associated with significant increases in IR levels (<xref ref-type="fig" rid="f2">Fig. 2c–e</xref>). For <italic>Lmnb1</italic>, the IR levels were significantly increased by >four-fold with (0 μM Aza: 7 μM Aza +caffeine, <italic>P</italic>=0.0182, <italic>t</italic>-test) and without NMD inhibition (0 μM Aza: 7 μM Aza, <italic>P</italic>=0.0136, <italic>t</italic>-test; <xref ref-type="fig" rid="f2">Fig. 2c</xref>). For <italic>Ngp</italic> and <italic>Spata13</italic>, significant increases in IR levels by >6 (<italic>P</italic>=0.0374, <italic>t</italic>-test) and >4-fold (<italic>P</italic>=0.0208, <italic>t</italic>-test), respectively, were achieved following NMD inhibition (<xref ref-type="fig" rid="f2">Fig. 2d,e</xref>). For the remaining six transcripts, IR levels did not change following 5-Aza treatments (<xref ref-type="fig" rid="f2">Fig. 2f</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4c–g</xref>). For five genes (<italic>Mmp8</italic>, <italic>S100a9</italic>, <italic>Gart</italic>, <italic>Gadd45a</italic> and <italic>Atf4</italic>), these results may be explained by the lack of changes in DNA methylation levels near 5′ or 3′ splice sites despite 5-Aza exposure (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4c–g</xref>). For <italic>S100a8</italic>, constitutive DNA methylation levels were low in untreated MPRO cells (average DNA methylation of 31% compared to 76% in primary promyelocytes), and further demethylation by 5-Aza failed to alter IR levels (<xref ref-type="fig" rid="f2">Fig. 2f</xref>). Notably, treatment with a lower concentration of 3 μM 5-Aza- was insufficient to cause substantial decrease in DNA methylation levels near 5′ or 3′ splice sites of all 10 introns (<xref ref-type="fig" rid="f2">Fig. 2c–f</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4c–g</xref>). We conclude that because no significant changes in DNA methylation levels were detected in cells treated with lower dose 5-Aza compared to untreated controls, there was no significant difference produced in IR levels. Overall, our data indicate that when DNA methylation levels near 5′ or 3′ splice sites flanking introns are substantially reduced, there is an increase in IR levels.</p></sec><sec disp-level="2"><title>IR splice sites show reduced MeCP2 binding</title><p>Reduced DNA methylation levels near splice junctions can reduce the binding of MeCP2 (ref. <xref ref-type="bibr" rid="b33">33</xref>), which has previously been shown to regulate exon inclusion<xref ref-type="bibr" rid="b12">12</xref>. To determine whether altered MeCP2 binding is involved in the regulation of IR, we performed genome-wide MeCP2 chromatin immunoprecipitation sequencing (ChIP-seq) in mouse promyelocytes and granulocytes. Consistent with IR-associated reductions in DNA methylation levels (<xref ref-type="fig" rid="f1">Fig. 1a</xref>), a significantly lower density of MeCP2 binding was observed near 5′ and 3′ splice junctions, and in the middle of retained compared to non-retained introns (promyelocytes and granulocytes combined, <italic>P</italic>=9.55e-15, 1.24e-14 and <italic>P</italic><1.00e-16 respectively, <italic>t</italic>-test, <xref ref-type="fig" rid="f3">Fig. 3a</xref>).</p><p>Using chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR), we measured the specific levels of MeCP2 binding near splice junctions flanking twelve introns that are subject to DNA methylation-associated increased IR in granulocytes compared to promyelocytes (<xref ref-type="fig" rid="f3">Fig. 3b</xref>). Significantly reduced levels of MeCP2 binding near splice junctions (<italic>P</italic><0.05, <italic>t</italic>-test) were observed in 8/12 introns (<xref ref-type="fig" rid="f3">Fig. 3b</xref>) in association with reduced DNA methylation levels at these sites (<xref ref-type="fig" rid="f1">Fig. 1c,d</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 3</xref>). No difference in MeCP2 binding was found near splice junctions of four randomly selected non-retained introns, or in the positive and negative controls for MeCP2-ChIP-qPCR (<xref ref-type="fig" rid="f3">Fig. 3b</xref>). In 5-Aza-treated MPRO cells, reduced MeCP2 binding near splice junctions flanking <italic>Lmnb1</italic> intron 5 and <italic>Ngp</italic> intron 3 was also observed (<xref ref-type="fig" rid="f3">Fig. 3c</xref>); concomitant with reduced DNA methylation levels and increased IR levels (<xref ref-type="fig" rid="f2">Fig. 2c,d</xref>). MeCP2 occupancy at the 3′ splice junction of <italic>Atf4</italic> intron 2 did not change (<xref ref-type="fig" rid="f3">Fig. 3c</xref>); consistent with the unaltered low or absent DNA methylation levels following 5-Aza treatment (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 4g</xref>). These data support the direct involvement of DNA methylation-dependent MeCP2 binding in regulating IR of these genes.</p><p>To determine if MeCP2 directly regulates IR more generally, we compared differential IR between matched wild type and <italic>MeCP2</italic>-knockdown and -knockout cells and tissues by analysing publicly available RNA-seq data. Of the differential IR events, we detected predominantly increased IR in all cell and tissue types examined including: IMR90 (719/1022, 70.35%), mouse cerebellum (134/206, 65.05%), mouse peritoneal macrophage (59/91 64.84%) and mouse visual cortex cells (32/46, 69.57%; <xref ref-type="fig" rid="f3">Fig. 3d</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5a,b</xref>). These results indicate that depletion of MeCP2 significantly enhances IR in diverse cell and tissue types (<italic>P</italic><0.05 in all cases, binomial test, <xref ref-type="fig" rid="f3">Fig. 3d</xref>, and <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5a,b</xref>).</p></sec><sec disp-level="2"><title>IR is associated with increased RNA Pol II stalling</title><p>Accumulation of RNA Pol II in alternatively spliced exons can promote their inclusion<xref ref-type="bibr" rid="b13">13</xref>. This process could be mediated through increased MeCP2 densities within exons, which act as a ‘go slow' signal to promote splice site recognition<xref ref-type="bibr" rid="b12">12</xref>. Accumulation of the elongating form of RNA Pol II (Pol II Ser2p) near splice junctions has been associated with increased IR levels<xref ref-type="bibr" rid="b9">9</xref>. We analysed the association between MeCP2 and Pol II Ser2p occupancies within differential intron-retaining genes in promyelocytes and granulocytes using ChIP-qPCR. Consistent with the pattern of reduced DNA methylation (<xref ref-type="fig" rid="f1">Fig. 1c,d</xref>), significantly lower densities of MeCP2 were observed near the 5′ splice junctions of the three introns associated with higher IR levels in granulocytes (<italic>P</italic><0.05, <italic>t</italic>-test; <xref ref-type="fig" rid="f4">Fig. 4a</xref>). For intron 3 of <italic>Ngp</italic> and intron 2 of <italic>S100a8</italic> in granulocytes, the biggest differences in MeCP2 densities were detected near the 3′ splice junctions (<italic>P</italic><0.05, <italic>t</italic>-test; <xref ref-type="fig" rid="f4">Fig. 4a</xref>). Pol II Ser2p densities were elevated in introns subject to higher IR levels in granulocytes (<xref ref-type="fig" rid="f4">Fig. 4b</xref>). Furthermore, Pol II Ser2p accumulated at the 3′ splice junctions of all three retained introns (<italic>P</italic><0.05, <italic>t</italic>-test; <xref ref-type="fig" rid="f4">Fig. 4b</xref>), consistent with a previous analysis of published data<xref ref-type="bibr" rid="b9">9</xref>. As a control, we noted that MeCP2 and Pol II Ser2p densities were not significantly different between promyelocytes and granulocytes within constitutively spliced intron 3 of <italic>Lmnb1</italic>, intron 9 of <italic>Smarcd1</italic> and intron 5 of <italic>Tbp</italic>, and their flanking exons (<xref ref-type="fig" rid="f4">Fig. 4c,d</xref>). Collectively, our data indicate that lower DNA methylation levels and consequent reduced MeCP2 binding to the DNA that encodes retained introns is associated with increased RNA Pol II stalling within these introns, predominantly near the 3′ splice junctions.</p></sec><sec disp-level="2"><title>IR is associated with reduced splicing factor recruitment</title><p>Our results (<xref ref-type="fig" rid="f4">Fig. 4a,b</xref>) indicate an anti-correlation between MeCP2 and Pol II Ser2p densities within introns and near their flanking splice junctions. Thus, the lower MeCP2 levels in granulocytes could not act as a ‘go slow' signal to promote Pol II Ser2p stalling for efficient splice signal recognition. Stalling of Pol II Ser2p has been reported as a consequence of reduced splicing factor recruitment<xref ref-type="bibr" rid="b15">15</xref><xref ref-type="bibr" rid="b16">16</xref>. Since MeCP2 has previously been shown to bind or recruit splicing factors to pre-mRNAs<xref ref-type="bibr" rid="b34">34</xref><xref ref-type="bibr" rid="b35">35</xref><xref ref-type="bibr" rid="b36">36</xref>, we hypothesize that Pol II Ser2p stalling in retained introns is consequent to reduced MeCP2-mediated recruitment of splicing factors. To determine splicing factors that associate with MeCP2 in mouse promyelocytes and granulocytes, we performed Rapid Immunoprecipitation Mass Spectrometry of Endogenous Proteins (RIME; <xref ref-type="fig" rid="f5">Fig. 5a</xref>).</p><p>RIME identified 39 proteins that co-immunoprecipitated with MeCP2 in promyelocytes and or granulocytes, of which 14 (36%) have previously been reported as MeCP2-binding partners (<xref ref-type="supplementary-material" rid="S1">Supplementary Table 2</xref>). Histone and ribosomal RNA proteins are amongst those associated with MeCP2 in these cell types, consistent with previous reports<xref ref-type="bibr" rid="b35">35</xref><xref ref-type="bibr" rid="b37">37</xref>. Five of these proteins were detected in both cell types (<xref ref-type="fig" rid="f5">Fig. 5b</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Table 2</xref>). Seven splicing factors (five from the Sr family proteins<xref ref-type="bibr" rid="b38">38</xref>) co-immunoprecipitated with MeCP2 in promyelocytes but not granulocytes (<xref ref-type="fig" rid="f5">Fig. 5b</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Table 2</xref>), indicating that one or more of these splicing factors may be associated with MeCP2-mediated splicing events specific to promyelocytes where IR is low. Gene ontology analysis also demonstrated a statistically significant enrichment of splicing-associated functions for proteins that were bound to MeCP2 in promyelocytes but not granulocytes (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5c,d</xref>). Of these proteins, Tra2b has recently been shown to be important for efficient AS in human and xenopus cells<xref ref-type="bibr" rid="b39">39</xref><xref ref-type="bibr" rid="b40">40</xref>, with IR being the dominant mode of splicing affected when Tra2b is knocked down<xref ref-type="bibr" rid="b39">39</xref>.</p><p>We further examined whether reduced MeCP2-mediated Tra2b recruitment could account for reduced splicing leading to IR. Another consequence of this model would be that Tra2b is recruited by MeCP2 to promote splicing in promyelocytes. There is no substantial difference in the MeCP2 and Tra2b protein levels in promyelocytes and granulocytes (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5e,f</xref>). RNA immunoprecipitation (RNA-IP) followed by qPCR with reverse transcription (RT-qPCR) showed that binding of MeCP2 and Tra2b to their known interacting RNA, <italic>Malat1</italic> and <italic>Tra2a</italic>, was not significantly different in promyelocytes and granulocytes (<xref ref-type="fig" rid="f5">Fig. 5c,d</xref>). In contrast, there was a significantly reduced association of MeCP2 with <italic>Lmnb1</italic> and <italic>Ngp</italic> transcripts, both of which showed higher IR levels in granulocytes (<italic>P</italic><0.05, <italic>t</italic>-test; <xref ref-type="fig" rid="f5">Fig. 5e</xref>). A non-significant trend of reduced MeCP2 binding was also observed at <italic>S100a8</italic> in granulocytes (<xref ref-type="fig" rid="f5">Fig. 5e</xref>). Importantly, the pattern of Tra2b enrichment copied that of MeCP2 (<xref ref-type="fig" rid="f5">Fig. 5f</xref>), indicating that these two proteins function together to promote splicing. The levels of association between MeCP2 and Tra2b with constitutively spliced <italic>Smarcd1</italic> and <italic>Tbp</italic> transcripts were not significantly different in promyelocytes and granulocytes (<xref ref-type="fig" rid="f5">Fig. 5g,h</xref>).</p><p>We further investigated whether Tra2b and MeCP2 act together in a complex that regulates splicing by performing sequential RNA-IP with antibodies against Tra2b and MeCP2 (RNA-IP-reIP). Our RNA-IP-reIP assays consistently detected de-enrichment of concomitant Tra2b and MeCP2 binding to <italic>Lmnb1</italic>, <italic>Ngp</italic> and <italic>S100a8</italic> transcripts that increased in IR levels in granulocytes compared to promyelocytes (<italic>P</italic><0.05, <italic>t</italic>-test; <xref ref-type="fig" rid="f5">Fig. 5i</xref>). The levels of concomitant Tra2b and MeCP2 binding to non intron-retaining <italic>Smarcd1</italic> and <italic>Tbp</italic> transcripts were not significantly different (<xref ref-type="fig" rid="f5">Fig. 5i</xref>). Thus, when MeCP2 and Tra2b reduce their interaction with mRNA there is a consequent increase in IR.</p><p>Our data support a model where MeCP2 regulates splicing via recruitment of splicing factors, such as Tra2b (<xref ref-type="fig" rid="f6">Fig. 6</xref>). IR may occur in cases where lack of MeCP2 binding to the DNA and RNA results in less efficient recruitment of the Tra2b and potentially other splicing factors (<xref ref-type="fig" rid="f6">Fig. 6</xref>).</p></sec></sec><sec disp-level="1"><title>Discussion</title><p>Until recently, IR has been the least understood mode of splicing in terms of its role in normal physiology and disease. As a result of recent work by our group and others, it is now established that IR is a conserved mechanism that regulates cell differentiation and plasticity in normal tissues and is implicated in diverse human cancers<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b8">8</xref><xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b41">41</xref><xref ref-type="bibr" rid="b42">42</xref><xref ref-type="bibr" rid="b43">43</xref>. For example, IR underlies normal differentiation and physiological cellular functions in granulopoiesis, erythropoiesis, megakaryopoiesis and lymphopoiesis<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b7">7</xref><xref ref-type="bibr" rid="b8">8</xref><xref ref-type="bibr" rid="b44">44</xref>. It is also a mechanism of tumour suppressor gene inactivation in diverse cancers<xref ref-type="bibr" rid="b43">43</xref>. In this work we sought to understand the molecular mechanisms underlying IR as a foundation to establish why specific introns are retained whereas others are not.</p><p>Since DNA methylation is known to regulate other forms of splicing, we explored its role in controlling IR. Two models of splicing regulation via DNA methylation have been proposed: ‘kinetic' and ‘recruitment'. The former involves ‘go slow' signal proteins reflecting DNA methylation changes, to reduce the rate of RNA Pol II elongation and promote recognition of splice sites<xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b13">13</xref>. The ‘recruitment model' applies to cases where DNA methylation directly recruits splicing factors to promote splicing<xref ref-type="bibr" rid="b18">18</xref>. Our results demonstrate a previously unrecognized link between the ‘kinetic' and ‘recruitment' models of splicing regulation, revealed by examining protein complexes associated with DNA methylation. We have established that reduced DNA methylation decreases the binding of MeCP2 to the DNA and RNA, which in turn reduces splicing and consequently increases IR (<xref ref-type="fig" rid="f3">Figs 3</xref> and <xref ref-type="fig" rid="f4">4</xref>). In support of this model we showed that splicing factors, such as Tra2b, are less efficiently recruited to RNA, and are associated with RNA Pol II stalling near the splice junctions flanking retained introns and within retained introns themselves (<xref ref-type="fig" rid="f4">Fig. 4</xref>). Subsequently, splicing repressors may be recruited following RNA Pol II stalling to stabilize IR (ref. <xref ref-type="bibr" rid="b9">9</xref>).</p><p>Our data indicate that DNA methylation changes are associated with approximately one in every six differential IR events in myeloid cells. This observation reveals the complexity of splicing and IR regulation, which may involve many different epigenetic factors and associated proteins. Indeed, H3.3K36me3 mediated recruitment of BS69 is known to promote IR by antagonizing the activity of the spliceosome via its interaction with snRNPs, including EFTUD2 (ref. <xref ref-type="bibr" rid="b10">10</xref>). While BS69 has also been shown to repress RNA Pol II elongation<xref ref-type="bibr" rid="b45">45</xref>, this work was performed in U20S cells rather than HeLa cells where the IR-promoting role of BS69 was reported. Thus, the possibility that BS69 plays a dual role in linking the ‘kinetic' and ‘recruitment' models of splicing remains to be determined. BS69 also binds to several core splicing factors including PRPF8, SRSF1 and SRSF4 (ref. <xref ref-type="bibr" rid="b10">10</xref>). BS69 may act as a physical barrier that prevents splicing factors from accessing nascent RNA while promoting the recruitment of splicing repressors via inhibition of transcription elongation. Alternatively, repression of transcription elongation may result from less efficient splicing factor recruitment, as we proposed in our model of IR regulation based on DNA methylation and MeCP2 (<xref ref-type="fig" rid="f6">Fig. 6</xref>).</p><p>It is important to recognize that AS can also be modulated by the ratio of DNA methylation to hydroxymethylation<xref ref-type="bibr" rid="b46">46</xref>. While we have utilized bisulfite sequencing, which cannot discriminate between DNA methylation and DNA hydroxymethylation, the levels of DNA hydroxymethylation in granulocytes are low<xref ref-type="bibr" rid="b47">47</xref>. Thus, the majority of unconverted cytosines we detected in granulocytes following bisulfite sequencing are most likely methylated cytosines. Whether DNA hydroxymethylation also regulates IR is unknown. This question will need to be explored in cells or tissue with high levels of DNA hydroxymethylation, and using sequencing methods that specifically detect this type of DNA modification<xref ref-type="bibr" rid="b48">48</xref>.</p><p>Although intragenic DNA demethylation has been associated with reduced gene expression in cancers<xref ref-type="bibr" rid="b49">49</xref>, it remains to be shown whether this occurs via aberrant IR-triggered NMD. The DNA methylation inhibitor, 5-Aza, is FDA-approved for the treatment of myelodysplastic syndromes and is in clinical trials for diverse cancers<xref ref-type="bibr" rid="b50">50</xref><xref ref-type="bibr" rid="b51">51</xref>. The extensive association between reduced DNA methylation and increased IR that we have shown in diverse tissue types indicates that global demethylation induced by 5-Aza may exert effects, which promote aberrant IR.</p><p>In conclusion, we have discovered a causative role of DNA methylation changes in the regulation of IR via MeCP2 binding and splicing factor recruitment. Our findings reveal another mechanism that controls splicing and IR in diverse tissue types. Further understanding of this mechanism in the development of cancers, and its impact in patients receiving DNA methylation inhibitors may contribute to improved therapeutic strategies.</p><p><italic>Note added in proof:</italic> At the proof stage of this manuscript, our software article detailing the IRFinder algorithm, Middleton <italic>et al</italic>.<xref ref-type="bibr" rid="b52">52</xref>, was published.</p></sec><sec disp-level="1"><title>Methods</title><sec disp-level="2"><title>Primary cells and cell lines</title><p>Eight–10-week-old C57BL6/J mice were obtained from the Animal Resource Centre (Perth, Australia). All mouse experiments were performed in a specific pathogen-free facility with approval from the Royal Prince Alfred Hospital Animal Ethics Committee (Protocol no. 2013/045). The FACS strategies for promyelocytes and granulocytes are shown in <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1a</xref> as previously published<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b53">53</xref>. Briefly, harvested bone marrow cells were incubated with lineage specific antibodies conjugated to biotin: B220 (BioLegend 103204, 1:200), CD19 (BioLegend 115504, 1:200), CD3 (BioLegend 100304, 1:200), Sca–1(BioLegend 108104, 1:200), together with anti–Gr–1–fluorescein isothiocyanate (FITC; Biolegend 108406, 1:200), anti–cKit–phycoerythrin (PE; Becton Dickinson 553355, 1:200), anti–CD34–Alexafluor 660 (eBioscience 50-0341-82, 1:100) and anti–CD16/32–PerCP/Cy5.5 (Becton Dickinson 560540, 1:400) antibodies. Cells were washed twice in phosphate-buffered saline (PBS) with 2% (volume per volume) foetal bovine serum (HyClone, GE Healthcare Life Sciences) and incubated with streptavidin–APC–Cy7 (BioLegend 405208, 1:1,000). Following two additional washes in PBS with 2% (volume per volume) foetal bovine serum, the cells were sorted into 15 ml collection tubes on a BD Influx Cell Sorter (BD Biosciences). Morphological confirmation of promyelocytes and granulocytes was performed following May-Grünwald Giemsa staining of cells spun onto poly-<sc>L</sc>-lysine slides (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 1b</xref>). MPRO cells were maintained in DMEM media supplemented with 5% foetal bovine serum (Hyclone, GE Healthcare Life Sciences), 10% BHK-HM5 conditioned medium made in-house by collecting and filtering the DMEM media in which BHK/HM5 cells were cultured for 96 h to produce granulocyte macrophage colony-stimulating factor, 5 U ml<sup>−1</sup> penicillin, 5 μg per ml streptomycin sulfate and 2 mM <sc>L</sc>-glutamine.</p></sec><sec disp-level="2"><title>DNA methylation analysis</title><p>Genomic DNA extracted from cells was bisulfite-converted using the EZ DNA Methylation-Gold Kit (Zymo Research, Orange, CA, USA), according to the manufacturers' instructions. For WGBS analyses, libraries were prepared from bisulfite-converted samples and sequenced on an Illumina HiSeq2000 platform at Macrogen Inc. (Korea). For gene-specific clonal bisulfite sequencing, bisulfite-converted DNA was amplified using primers that equally amplify methylated and unmethylated template. Primer sequences are provided in <xref ref-type="supplementary-material" rid="S1">Supplementary Table 3</xref>. PCR amplicons were cloned into pGEM-T Easy vectors (Promega, Madison WI USA) and sequenced using Sanger Sequencing.</p></sec><sec disp-level="2"><title>Inhibition of DNA methylation</title><p>Inhibition of DNA methylation in MPRO cells was performed using 3 μM or 7 μM 5-Aza-2′ deoxycytidine for 120 h with drugs replenished every 24 h. Dead cells were excluded by FACS following 4′,6-diamidino-2-phenylindole staining.</p></sec><sec disp-level="2"><title>RNA isolation and qRT-PCR analyses</title><p>Total RNA extraction was performed using Trizol (Invitrogen). RNA 6000 Nano Chips (Agilent Technologies) were used to confirm RNA quality. For qRT-PCR, total RNA was treated with DNAse I, amplification grade (Invitrogen). cDNA was generated using SuperScript III First-Strand Synthesis System (Invitrogen), according to the manufacturers' instructions. qRT-PCR was performed using the CFX96 Real Time PCR Machine (Bio-Rad). Reactions were performed in 20 μl volumes containing 1 × IQ SyberGreen supermix and 0.3 μM of the respective forward and reverse primers. Cycling conditions were: 95 °C for 6 min followed by 95 °C for 30 s, primer-specific annealing temperature for 30 s, and extension at 72 °C for 30 s. Primer sequences are shown in <xref ref-type="supplementary-material" rid="S1">Supplementary Table 3</xref>.</p></sec><sec disp-level="2"><title>Inhibition of NMD using caffeine coupled with actinomycin D</title><p>NMD inhibition was performed on MPRO cells after 5-Aza-2′ deoxycytidine treatment as previously published<xref ref-type="bibr" rid="b2">2</xref>. Briefly, cells in culture media were treated with caffeine (15 mM) for 4 h before they were washed twice with PBS. Caffeine and actinomycin D (2 μg ml<sup>−1</sup>) were then added, and an actinomycin <sc>D</sc>-only control was used to block transcription of genes. Following an additional 4 h of incubation, total RNA was collected for qRT-PCR.</p></sec><sec disp-level="2"><title>Western blot</title><p>Total protein lysates were loaded onto precast SDS-PAGE gels (Invitrogen) and subjected to electrophoresis before being transferred onto polyvinylidene difluoride (PVDF) membranes. Membranes were blocked with 5% (weight per volume) skim milk for 1 h at room temperature and incubated overnight with a rabbit anti-MeCP2 (5 μg, Abcam ab2828; <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5e</xref>, a rabbit anti-Tra2b (10 μg, Abcam ab31353; <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5f</xref>), or a mouse anti-Gapdh (1:1,000, Abcam ab8245; <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 5e,f</xref>). Following washes, membranes were incubated with HRP-conjugated secondary anti-rabbit or anti-mouse antibody (1:5,000; Chemicon AP182P or AP192P), and exposed using enhanced chemiluminescence reagents (Pierce) on a Kodak Imager (Kodak). Bands were quantified using ImageJ. All uncropped blots are displayed in <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 6</xref>.</p></sec><sec disp-level="2"><title>Chromatin immunoprecipitation</title><p>ChIP was performed with either the Magna ChIP A/G kit (Millipore 17-10085) or Magna ChIP HiSens kit (Millipore 17-10460) following the manufacturers' instructions, with minor modifications. Briefly, 5 × 10e<sup>6</sup> FACS-enriched promyelocytes or granulocytes were cross-linked with 0.5% formaldehyde for 10 min at room temperature, and quenched with 1 × glycine for another 5 min. Cells were washed twice with ice-cold PBS, lysed and sonicated using Bioruptor Plus (Diagenode, UCD-300) with 15 cycles of 30 second ON and 30 s OFF at high power setting. Sonicated samples were centrifuged for 10 min at 10,000<italic>g</italic>. Supernatant was incubated overnight at 4 °C with the antibody of interest mixed with protein A/G magnetic beads. Following washes, crosslinking was reversed at 62 °C for 2 h and DNA was purified using the Agencourt AMPure XP beads (Beckman Coulter, A63881). Antibodies used were rabbit anti-MeCP2 (5 μg, Abcam ab2828), rat anti-Pol II Ser2p (7.5 μg, Active Motif 61084), rabbit IgG control (5 μg Abcam ab46540) and normal rat IgG (7.5 μg, Santa Cruz sc-3036). ChIP DNA was analysed by qPCR with primers listed in <xref ref-type="supplementary-material" rid="S1">Supplementary Table 3</xref>. ChIP-qPCR results were reported as percentage of input relative to negative controls comprising a CpG-rich promoter region of Gapdh and an intergenic region on mouse chromosome 6 for MeCP2 and Pol II Ser2p respectively.</p><p>For ChIP-seq, sample concentration and size distribution of immunoprecipitated DNA were measured with Qubit dsDNA HS Assay Kits (Invitrogen) and High Sensitivity DNA Analysis Kit (Agilent Technologies), respectively. Approximately 20 ng of the input and ChIP DNA samples were sent to Macrogen Inc. (Korea) for ChIP-seq library preparation and sequencing on an Illumina Hiseq2000 platform.</p></sec><sec disp-level="2"><title>Rapid immunoprecipitation and mass spectrometry</title><p>Rapid immunoprecipitation and MS was performed as previously reported<xref ref-type="bibr" rid="b54">54</xref>, except that the chromatin immunoprecipitation step was performed using the Magna Chip A/G kit (Millipore 17-10085). Briefly, cells were sonicated and incubated overnight with antibody of interest mixed with protein A/G magnetic beads as described in the chromatin immunoprecipitation Method section. Following washing of antibody–protein-beads complex with buffers provided in the kit, trypsin digestion was perform to release the bound peptides. MS was performed on a QExactive Plus (Thermo Fisher). The data were searched within the Mouse SwissProt database using the Mascot search engine (version 2.4.0) with tryptic specificity allowing up to two missed cleavages, peptide mass tolerance of 10 p.p.m. and 0.8 Da for MS/MS ion mass tolerance. Variable modifications were oxidized methionine, acetylated asparagine and carbamidomethylated cysteine. IgG controls were used to exclude non-specific interactions. Only proteins present in three out of three biological replicates and not in any of the IgG controls were reported.</p></sec><sec disp-level="2"><title>RNA immunoprecipitation</title><p>RNA-IP was performed using the Magna Nuclear RNA-IP (cross-linked) kit (17-10520) from Millipore according to the manufacturer's protocol. Antibodies used were rabbit anti-MeCP2 (10 μg, Abcam ab2828) and rabbit anti-Tra2b (10 μg, Abcam ab31353). Approximately 1 million 0.5% formaldehyde-fixed cells were used for RNA-IP. Immunoprecipitated RNA was purified using Trizol LS Reagent (Life Technologies, 10296-010) and subjected to DNase I treatment. iTaq universal SYBR Green one-step kit (Bio-Rad) was used to quantify the amount of immunoprecipitated RNA. Primers used were listed in <xref ref-type="supplementary-material" rid="S1">Supplementary Table 3</xref>.</p></sec><sec disp-level="2"><title>RNA immunoprecipitation re immunoprecipitation (RNA-IP-reIP)</title><p>RNA-IP-reIP was performed using the same kit used for RNA-IP: Magna Nuclear RNA-IP (Cross-linked) kit (17-10520) from Millipore. Approximately 2 million 0.5% formaldehyde-fixed cells were used for each RNA-IP-reIP. Sequential RNA-IP was performed with an additional purification step between the two immunoprecipitation steps. Briefly, immunoprecipitation was performed on fixed and sonicated chromatin lysate using the rabbit anti-Tra2b antibody (20 μg, Abcam ab31353) together with A/G magnetic beads at 4 °C for 1 h. Following washes according to the manufacturer's instruction, the bead bound RNA-Protein complex were released with the elution buffer supplemented with DTT, RNase inhibitor and protease inhibitor for 15 min at 37 °C. SDS was removed from the eluent with HiPPR Detergent Removal spin column (Thermo Fisher, PIE88305) before the subsequent RNA-IP.</p><p>The eluent from the first RNA-IP was then used for a second RNA-IP at 4 °C for 1 h with either the Rabbit anti-MeCP2 (10 μg, Abcam ab2828) or the normal rabbit IgG (10 μg, Santa CruzSC-2027) as control. Following washes as recommended by the manufacturer, the immunoprecipitated RNA was purified using the Trizol LS Reagent (Life Technologies, 10296-010) and DNase I treatment was performed. iTaq universal SYBR Green one-step kit (Bio-Rad) was used to quantify the amount of immunoprecipitated RNA relative to the sequential Tra2b-IgG RNA-IP control. The workflow for RNA-IP-reIP is shown in <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. 7</xref>. Primers used were listed in <xref ref-type="supplementary-material" rid="S1">Supplementary Table 3</xref>.</p></sec><sec disp-level="2"><title>Bioinformatics analyses</title><p>Besides the WGBS and RNA-seq reads generated from promyelocytes and granulocytes, we also analysed WGBS, RNA-seq and MeCP2 ChIP-seq data from the GEO (Gene expression Omnibus) database (<xref ref-type="supplementary-material" rid="S1">Supplementary Table 1</xref>). All RNA-seq datasets have been generated from poly-A enriched libraries as published elsewhere<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b55">55</xref><xref ref-type="bibr" rid="b56">56</xref><xref ref-type="bibr" rid="b57">57</xref>. To assess genome-wide DNA methylation data, WGBS data was aligned to mouse genome GRCm38/mm10 and human genome GRCh37/hg19, respectively, followed by the calculation of methylation level (beta value) at each CpG site using Bismark v0.9.0 (ref. <xref ref-type="bibr" rid="b50">50</xref>) with default parameters. Only those CpG sites covered by more than five reads were considered as valid for the downstream analysis. A sliding window of 10 bp was then applied to scan the ±100 bp region of a splice junction or from the middle of an intron. Within each window, average methylation level was calculated as the sum of beta values divided by the number of CpGs. For gene expression analysis, TopHat v2.2.2 with default parameters was applied to map RNA-seq reads against Ensembl human (GRCh37.75) and mouse (GRCm38.73) transcriptomes<xref ref-type="bibr" rid="b58">58</xref>. Cufflinks v2.1.1 was then applied to measure the FPKM of each gene<xref ref-type="bibr" rid="b58">58</xref>.</p><p>To assess the MeCP2 binding density, ChIP-seq reads were mapped to human (hg19) or mouse (mm10) genomes using Bowtie2 v2.2.3 default settings. MAC2 (ref. <xref ref-type="bibr" rid="b59">59</xref>) was applied to call MeCP2 peaks from mapped ChIP-seq reads, with an extension size of 250nt and a <italic>q</italic>-value (FDR) cut-off of 0.05. To assess the MeCP2 binding density, ChIP-seq reads were mapped to human (hg19) or mouse (mm10) genomes using Bowtie2 v2.2.3 with default settings. MACS2 (ref. <xref ref-type="bibr" rid="b54">54</xref>) was applied to call MeCP2 peaks from mapped ChIP-seq reads, with an extension size of 250nt and a <italic>q</italic>-value (FDR) cut-off of 0.05. MeCP2 density at each peak was reported as fold enrichment, which is the number of position-adjusted reads in MeCP2-ChIP sample/input. A sliding window of 10 bp was then applied to scan the ±100 bp region of a splice junction or from the middle of an intron. Within each window, average MeCP2 signal was calculated as the sum of fold enrichment divided by the number of peaks.</p></sec><sec disp-level="2"><title>Intron retention discovery</title><p>The steps of detecting IR in an individual sample and alternate IR between samples are described below:</p><p>A STAR<xref ref-type="bibr" rid="b60">60</xref> genome index is built with the human or mouse genome fasta file (GRCh37/GRCm38) and annotation gtf file (Ensembl 75 for human, 78 for mouse). Features that overlap with introns such as microRNAs or snoRNAs are excluded as these may confound the accurate measurement of true intron levels. To accurately measure the expression levels of introns, intronic regions of poor mappability are determined by mapping synthetic reads to the genome and are excluded from further analysis. Synthetic reads are taken from the human or mouse genome reference, which are 70 bp long and stepped at 10 bp across the entire genome. Every second read is reverse complemented. A single base error is generated in the centre of each read. Those synthetic reads uniquely mapping to the correct location are tallied. Any area without at least 5 out of a possible 7 unique reads is considered poorly mappable. Poorly mappable regions are excluded from the measurable intron area regardless of strand or direction.</p><p>RNA-seq reads were aligned to the reference genome by STAR using default scoring parameters but excluding multi-mapping reads from the output. The default scoring parameters suit the detection of IR as they favour neither mapping from exon to exon across splice junctions nor mapping from exons into introns. IR is then determined by measurement of both the splicing level and intronic abundance. The key calculated metric is the IR ratio. Abundance of normal splicing is measured by a count of read fragments spliced across the intron. Reads that start in the 5′ exon but end in another exon and reads that start in another exon but end in the 3′ exon are also counted, the larger value being used. Intronic abundance is measured as a trimmed mean of fragment coverage depth. The trimmed mean excludes previously excluded regions (that is, regions with poor mappability, miRNAs or ncRNAs) and the highest and lowest 30% of values. This compensates for noise in intronic coverage. The IR ratio is calculated simply as follows:</p><p><disp-formula id="eq1"><inline-graphic id="d32e1264" xlink:href="ncomms15134-m1.jpg"></inline-graphic></disp-formula></p><p>To compare the change of IR of a certain intron between a pair of samples or groups, the intronic abundance and exonic splice abundance of this intron is evaluated using a Bayesian statistic adapted for digital counts<xref ref-type="bibr" rid="b61">61</xref>. Significance of the change in IR ratio is considered using the following cut-off: 1) <italic>P</italic>-value is <0.05; 2) exonic splice abundance for both samples or groups are above 10; 3) IR ratio in one sample or group is above 0.1.</p></sec><sec disp-level="2"><title>Gene ontology enrichment analysis</title><p>Gene ontology enrichment analysis was performed using DAVID<xref ref-type="bibr" rid="b62">62</xref>. Gene ontology categories with a <italic>P</italic>-value <0.01 after Benjamini Hochberg correction were considered significant.</p></sec><sec disp-level="2"><title>Statistical analysis</title><p>Statistical analyses were performed as indicated in the text and/or figure legends. To compare the overall DNA methylation level or MeCP2 fold enrichment across the ±100 bp flanking splice junctions, the whole region (200 bp) was first divided into 20 sub-windows (10 bp each). Within each sub-window, the mean of DNA methylation level or MeCP2 fold enrichment for all IR and non-IR events was calculated respectively and paired together. Paired <italic>t</italic>-test was then applied to assess the significance of overall difference across all the sub-windows.</p><p>The significance of the difference between the number of increased and decreased IR events in cells following DNA methylation inhibition or MeCP2 knockdown was determined using the binomial test. Student <italic>t</italic>-tests were used to compare differential gene expression and levels of GC content and CpG density between retained and non-retained introns. Differential enrichment of MeCP2, Pol II Ser2p and Tra2b assessed by qPCR following ChIP and/or RNA-IP was also determined using <italic>t</italic>-tests. <italic>F</italic>-test was used to determine the significance of regression analysis. Analyses were performed using Graph-Pad Prism v.6 (La Jolla, CA, USA). <italic>P</italic><0.05 was considered significant for all tests.</p></sec><sec disp-level="2"><title>Data availability</title><p>WGBS and ChIP-seq data that support the findings of this study have been deposited in Gene Expression Omnibus with the accession number GSE85517. All other data that support the findings of this study are available from the corresponding authors upon reasonable request.</p></sec></sec><sec disp-level="1"><title>Additional information</title><p><bold>How to cite this article:</bold> Wong, J. J.-L. <italic>et al</italic>. Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment. <italic>Nat. Commun.</italic><bold>8</bold>, 15134 doi: 10.1038/ncomms15134 (2017).</p><p><bold>Publisher's note:</bold> Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p></sec><sec sec-type="supplementary-material" id="S1"><title>Supplementary Material</title><supplementary-material id="d32e18" content-type="local-data"><caption><title>Supplementary Information</title><p>Supplementary Figures and Supplementary Tables</p></caption><media xlink:href="ncomms15134-s1.pdf"></media></supplementary-material></sec></body><back><ack><p>We thank Steven Allen, Frank Kao, Suat Dervish and Caroline Royle for FACS and Ben Crossett for MS. This work was supported by the National Health and Medical Research Council (Grant #1061906 to J.E.J.R., W.R., #1129901 to J.E.J.R. and J.J.-L.W. and #1080530 to J.E.J.R., W.R., J.J.-L.W.); Cure the Future (J.E.J.R.); Agence Nationale de la Recherche (Grant # ANR-16-CE12-0007 to W.R.) and an anonymous foundation (J.E.J.R.). J.J.-L.W. and W.R. hold Fellowships from the Cancer Institute of NSW. J.H. is a National Breast Cancer Foundation Fellow.</p></ack><ref-list><ref id="b1"><mixed-citation publication-type="journal"><name><surname>Yap</surname><given-names>K.</given-names></name>, <name><surname>Lim</surname><given-names>Z. Q.</given-names></name>, <name><surname>Khandelia</surname><given-names>P.</given-names></name>, <name><surname>Friedman</surname><given-names>B.</given-names></name> & <name><surname>Makeyev</surname><given-names>E. V.</given-names></name><article-title>Coordinated regulation of neuronal mRNA steady-state levels through developmentally controlled intron retention</article-title>. <source>Genes Dev.</source><volume>26</volume>, <fpage>1209</fpage>–<lpage>1223</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22661231</pub-id></mixed-citation></ref><ref id="b2"><mixed-citation publication-type="journal"><name><surname>Wong J.J.-L.</surname></name><etal></etal>. <article-title>Orchestrated intron retention regulates normal granulocyte differentiation</article-title>. <source>Cell</source><volume>154</volume>, <fpage>583</fpage>–<lpage>595</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23911323</pub-id></mixed-citation></ref><ref id="b3"><mixed-citation publication-type="journal"><name><surname>Moore</surname><given-names>M. J.</given-names></name>, <name><surname>Wang</surname><given-names>Q.</given-names></name>, <name><surname>Kennedy</surname><given-names>C. J.</given-names></name> & <name><surname>Silver</surname><given-names>P. A.</given-names></name><article-title>An alternative splicing network links cell-cycle control to apoptosis</article-title>. <source>Cell</source><volume>142</volume>, <fpage>625</fpage>–<lpage>636</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20705336</pub-id></mixed-citation></ref><ref id="b4"><mixed-citation publication-type="journal"><name><surname>Ergun</surname><given-names>A.</given-names></name><etal></etal>. <article-title>Differential splicing across immune system lineages</article-title>. <source>Proc. Natl Acad. Sci. USA</source><volume>110</volume>, <fpage>14324</fpage>–<lpage>14329</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23934048</pub-id></mixed-citation></ref><ref id="b5"><mixed-citation publication-type="journal"><name><surname>Shaked</surname><given-names>I.</given-names></name>, <name><surname>Zimmerman</surname><given-names>G.</given-names></name> & <name><surname>Soreq</surname><given-names>H.</given-names></name><article-title>Stress-induced alternative splicing modulations in brain and periphery</article-title>. <source>Ann. N. Y. Acad. Sci.</source><volume>1148</volume>, <fpage>269</fpage>–<lpage>281</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">19120119</pub-id></mixed-citation></ref><ref id="b6"><mixed-citation publication-type="journal"><name><surname>Padgett</surname><given-names>R. A.</given-names></name><article-title>New connections between splicing and human disease</article-title>. <source>Trends Genet.</source><volume>28</volume>, <fpage>147</fpage>–<lpage>154</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22397991</pub-id></mixed-citation></ref><ref id="b7"><mixed-citation publication-type="journal"><name><surname>Edwards</surname><given-names>C.</given-names></name><etal></etal>. <article-title>A dynamic intron retention program in the mammalian megakaryocyte and erythrocyte lineages</article-title>. <source>Blood</source><volume>127</volume>, <fpage>e23</fpage>–<lpage>e34</lpage> (<year>2016</year>).</mixed-citation></ref><ref id="b8"><mixed-citation publication-type="journal"><name><surname>Pimentel</surname><given-names>H.</given-names></name><etal></etal>. <article-title>A dynamic intron retention program enriched in RNA processing genes regulates gene expression during terminal erythropoiesis</article-title>. <source>Nucleic Acids Res.</source><volume>44</volume>, <fpage>838</fpage>–<lpage>851</lpage> (<year>2016</year>).<pub-id pub-id-type="pmid">26531823</pub-id></mixed-citation></ref><ref id="b9"><mixed-citation publication-type="journal"><name><surname>Braunschweig</surname><given-names>U.</given-names></name><etal></etal>. <article-title>Widespread intron retention in mammals functionally tunes transcriptomes</article-title>. <source>Genome Res.</source><volume>24</volume>, <fpage>1774</fpage>–<lpage>1786</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25258385</pub-id></mixed-citation></ref><ref id="b10"><mixed-citation publication-type="journal"><name><surname>Guo</surname><given-names>R.</given-names></name><etal></etal>. <article-title>BS69/ZMYND11 reads and connects histone H3.3 lysine 36 trimethylation-decorated chromatin to regulated pre-mRNA processing</article-title>. <source>Mol. Cell</source><volume>56</volume>, <fpage>298</fpage>–<lpage>310</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25263594</pub-id></mixed-citation></ref><ref id="b11"><mixed-citation publication-type="journal"><name><surname>de la Mata</surname><given-names>M.</given-names></name><etal></etal>. <article-title>A slow RNA Polymerase II affects alternative splicing <italic>in vivo</italic></article-title>. <source>Mol. Cell</source><volume>12</volume>, <fpage>525</fpage>–<lpage>532</lpage> (<year>2003</year>).<pub-id pub-id-type="pmid">14536091</pub-id></mixed-citation></ref><ref id="b12"><mixed-citation publication-type="journal"><name><surname>Maunakea</surname><given-names>A. K.</given-names></name>, <name><surname>Chepelev</surname><given-names>I.</given-names></name>, <name><surname>Cui</surname><given-names>K.</given-names></name>, <name><surname>Zhao</surname><given-names>K.</given-names></name> & <name><surname>Intragenic</surname><given-names>D. N. A.</given-names></name><article-title>methylation modulates alternative splicing by recruiting MeCP2 to promote exon recognition</article-title>. <source>Cell Res.</source><volume>23</volume>, <fpage>1256</fpage>–<lpage>1269</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23938295</pub-id></mixed-citation></ref><ref id="b13"><mixed-citation publication-type="journal"><name><surname>Shukla</surname><given-names>S.</given-names></name><etal></etal>. <article-title>CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing</article-title>. <source>Nature</source><volume>479</volume>, <fpage>74</fpage>–<lpage>79</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21964334</pub-id></mixed-citation></ref><ref id="b14"><mixed-citation publication-type="journal"><name><surname>Gascard</surname><given-names>P.</given-names></name><etal></etal>. <article-title>Epigenetic and transcriptional determinants of the human breast</article-title>. <source>Nat. Commun.</source><volume>6</volume>, <fpage>6351</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25690954</pub-id></mixed-citation></ref><ref id="b15"><mixed-citation publication-type="journal"><name><surname>Lin</surname><given-names>S.</given-names></name>, <name><surname>Coutinho-Mansfield</surname><given-names>G.</given-names></name>, <name><surname>Wang</surname><given-names>D.</given-names></name>, <name><surname>Pandit</surname><given-names>S.</given-names></name> & <name><surname>Fu</surname><given-names>X.-D.</given-names></name><article-title>The splicing factor SC35 has an active role in transcriptional elongation</article-title>. <source>Nat. Struct. Mol. Biol.</source><volume>15</volume>, <fpage>819</fpage>–<lpage>826</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18641664</pub-id></mixed-citation></ref><ref id="b16"><mixed-citation publication-type="journal"><name><surname>Fong</surname><given-names>Y. W.</given-names></name> & <name><surname>Zhou</surname><given-names>Q.</given-names></name><article-title>Stimulatory effect of splicing factors on transcription elongation</article-title>. <source>Nature</source><volume>414</volume>, <fpage>929</fpage>–<lpage>933</lpage> (<year>2001</year>).<pub-id pub-id-type="pmid">11780068</pub-id></mixed-citation></ref><ref id="b17"><mixed-citation publication-type="journal"><name><surname>Lev Maor</surname><given-names>G.</given-names></name>, <name><surname>Yearim</surname><given-names>A.</given-names></name> & <name><surname>Ast</surname><given-names>G.</given-names></name><article-title>The alternative role of DNA methylation in splicing regulation</article-title>. <source>Trends Genet.</source><volume>31</volume>, <fpage>274</fpage>–<lpage>280</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25837375</pub-id></mixed-citation></ref><ref id="b18"><mixed-citation publication-type="journal"><name><surname>Yearim</surname><given-names>A.</given-names></name><etal></etal>. <article-title>HP1 is involved in regulating the global impact of DNA methylation on alternative splicing</article-title>. <source>Cell Rep.</source><volume>10</volume>, <fpage>1122</fpage>–<lpage>1134</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25704815</pub-id></mixed-citation></ref><ref id="b19"><mixed-citation publication-type="journal"><name><surname>Olins</surname><given-names>A. L.</given-names></name><etal></etal>. <article-title>Nuclear envelope and chromatin compositional differences comparing undifferentiated and retinoic acid- and phorbol ester-treated HL-60 cells</article-title>. <source>Exp. Cell Res.</source><volume>268</volume>, <fpage>115</fpage>–<lpage>127</lpage> (<year>2001</year>).<pub-id pub-id-type="pmid">11478838</pub-id></mixed-citation></ref><ref id="b20"><mixed-citation publication-type="journal"><name><surname>Porfiri</surname><given-names>E.</given-names></name>, <name><surname>Hoffbrand</surname><given-names>A.</given-names></name> & <name><surname>Wickremasinghe</surname><given-names>R.</given-names></name><article-title>Retinoic acid-induced granulocytic differentiation of HL60 human promyelocytic leukemia cells is preceded by downregulation of autonomous generation of inositol lipid-derived second messengers</article-title>. <source>Blood</source><volume>78</volume>, <fpage>1069</fpage>–<lpage>1077</lpage> (<year>1991</year>).<pub-id pub-id-type="pmid">1651133</pub-id></mixed-citation></ref><ref id="b21"><mixed-citation publication-type="journal"><name><surname>Lodhi</surname><given-names>I. J.</given-names></name><etal></etal>. <article-title>Peroxisomal lipid synthesis regulates inflammation by sustaining neutrophil membrane phospholipid composition and viability</article-title>. <source>Cell Metab.</source><volume>21</volume>, <fpage>51</fpage>–<lpage>64</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25565205</pub-id></mixed-citation></ref><ref id="b22"><mixed-citation publication-type="journal"><name><surname>Crawford</surname><given-names>L. J. A.</given-names></name>, <name><surname>Peake</surname><given-names>R.</given-names></name>, <name><surname>Price</surname><given-names>S.</given-names></name>, <name><surname>Morris</surname><given-names>T. C. M.</given-names></name> & <name><surname>Irvine</surname><given-names>A. E.</given-names></name><article-title>Adiponectin is produced by lymphocytes and is a negative regulator of granulopoiesis</article-title>. <source>J. Leukocyte Biol.</source><volume>88</volume>, <fpage>807</fpage>–<lpage>811</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20643815</pub-id></mixed-citation></ref><ref id="b23"><mixed-citation publication-type="journal"><name><surname>Jones</surname><given-names>P. A.</given-names></name><article-title>Functions of DNA methylation: islands, start sites, gene bodies and beyond</article-title>. <source>Nat. Rev. Genet.</source><volume>13</volume>, <fpage>484</fpage>–<lpage>492</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22641018</pub-id></mixed-citation></ref><ref id="b24"><mixed-citation publication-type="journal"><name><surname>Li-Byarlay</surname><given-names>H.</given-names></name><etal></etal>. <article-title>RNA interference knockdown of DNA methyl-transferase 3 affects gene alternative splicing in the honey bee</article-title>. <source>Proc. Natl Acad. Sci. USA</source><volume>110</volume>, <fpage>12750</fpage>–<lpage>12755</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23852726</pub-id></mixed-citation></ref><ref id="b25"><mixed-citation publication-type="journal"><name><surname>Laurent</surname><given-names>L.</given-names></name><etal></etal>. <article-title>Dynamic changes in the human methylome during differentiation</article-title>. <source>Genome Res.</source><volume>20</volume>, <fpage>320</fpage>–<lpage>331</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20133333</pub-id></mixed-citation></ref><ref id="b26"><mixed-citation publication-type="journal"><name><surname>Gombart</surname><given-names>A. F.</given-names></name><etal></etal>. <article-title>Regulation of neutrophil and eosinophil secondary granule gene expression by transcription factors C/EBP epsilon and PU.1</article-title>. <source>Blood</source><volume>101</volume>, <fpage>3265</fpage>–<lpage>3273</lpage> (<year>2003</year>).<pub-id pub-id-type="pmid">12515729</pub-id></mixed-citation></ref><ref id="b27"><mixed-citation publication-type="journal"><name><surname>Nacken</surname><given-names>W.</given-names></name>, <name><surname>Roth</surname><given-names>J.</given-names></name>, <name><surname>Sorg</surname><given-names>C.</given-names></name> & <name><surname>Kerkhoff</surname><given-names>C.</given-names></name><article-title>S100A9/S100A8: myeloid representatives of the S100 protein family as prominent players in innate immunity</article-title>. <source>Microsc. Res. Tech.</source><volume>60</volume>, <fpage>569</fpage>–<lpage>580</lpage> (<year>2003</year>).<pub-id pub-id-type="pmid">12645005</pub-id></mixed-citation></ref><ref id="b28"><mixed-citation publication-type="journal"><name><surname>Nathan</surname><given-names>C.</given-names></name><article-title>Neutrophils and immunity: challenges and opportunities</article-title>. <source>Nat. Rev. Immunol.</source><volume>6</volume>, <fpage>173</fpage>–<lpage>182</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16498448</pub-id></mixed-citation></ref><ref id="b29"><mixed-citation publication-type="journal"><name><surname>Chan</surname><given-names>K. T.</given-names></name>, <name><surname>Creed</surname><given-names>S. J.</given-names></name> & <name><surname>Bear</surname><given-names>J. E.</given-names></name><article-title>Unraveling the enigma: progress towards understanding the coronin family of actin regulators</article-title>. <source>Trends Cell Biol.</source><volume>21</volume>, <fpage>481</fpage>–<lpage>488</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21632254</pub-id></mixed-citation></ref><ref id="b30"><mixed-citation publication-type="journal"><name><surname>Salerno</surname><given-names>D. M.</given-names></name>, <name><surname>Tront</surname><given-names>J. S.</given-names></name>, <name><surname>Hoffman</surname><given-names>B.</given-names></name> & <name><surname>Liebermann</surname><given-names>D. A.</given-names></name><article-title>Gadd45a and Gadd45b modulate innate immune functions of granulocytes and macrophages by differential regulation of p38 and JNK signaling</article-title>. <source>J. Cell Physiol.</source><volume>227</volume>, <fpage>3613</fpage>–<lpage>3620</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22307729</pub-id></mixed-citation></ref><ref id="b31"><mixed-citation publication-type="journal"><name><surname>Gombart</surname><given-names>A. F.</given-names></name>, <name><surname>Grewal</surname><given-names>J.</given-names></name> & <name><surname>Koeffler</surname><given-names>H. P.</given-names></name><article-title>ATF4 differentially regulates transcriptional activation of myeloid-specific genes by C/EBPɛ and C/EBPα</article-title>. <source>J. Leukoc. Biol.</source><volume>81</volume>, <fpage>1535</fpage>–<lpage>1547</lpage> (<year>2007</year>).<pub-id pub-id-type="pmid">17347301</pub-id></mixed-citation></ref><ref id="b32"><mixed-citation publication-type="journal"><name><surname>Baubec</surname><given-names>T.</given-names></name><etal></etal>. <article-title>Genomic profiling of DNA methyltransferases reveals a role for DNMT3B in genic methylation</article-title>. <source>Nature</source><volume>520</volume>, <fpage>243</fpage>–<lpage>247</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25607372</pub-id></mixed-citation></ref><ref id="b33"><mixed-citation publication-type="journal"><name><surname>Baubec</surname><given-names>T.</given-names></name>, <name><surname>Ivánek</surname><given-names>R.</given-names></name>, <name><surname>Lienert</surname><given-names>F.</given-names></name> & <name><surname>Schübeler</surname><given-names>D.</given-names></name><article-title>Methylation-dependent and -independent genomic targeting principles of the MBD protein family</article-title>. <source>Cell</source><volume>153</volume>, <fpage>480</fpage>–<lpage>492</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23582333</pub-id></mixed-citation></ref><ref id="b34"><mixed-citation publication-type="journal"><name><surname>Young</surname><given-names>J. I.</given-names></name><etal></etal>. <article-title>Regulation of RNA splicing by the methylation-dependent transcriptional repressor methyl-CpG binding protein 2</article-title>. <source>Proc. Natl Acad. Sci. USA</source><volume>102</volume>, <fpage>17551</fpage>–<lpage>17558</lpage> (<year>2005</year>).<pub-id pub-id-type="pmid">16251272</pub-id></mixed-citation></ref><ref id="b35"><mixed-citation publication-type="journal"><name><surname>Maxwell</surname><given-names>S. S.</given-names></name>, <name><surname>Pelka</surname><given-names>G. J.</given-names></name>, <name><surname>Tam</surname><given-names>P. P. L.</given-names></name> & <name><surname>El-Osta</surname><given-names>A.</given-names></name><article-title>Chromatin context and ncRNA highlight targets of MeCP2 in brain</article-title>. <source>RNA Biol.</source><volume>10</volume>, <fpage>1741</fpage>–<lpage>1757</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">24270455</pub-id></mixed-citation></ref><ref id="b36"><mixed-citation publication-type="journal"><name><surname>Long</surname><given-names>S. W.</given-names></name>, <name><surname>Ooi</surname><given-names>J. Y. Y.</given-names></name>, <name><surname>Yau</surname><given-names>P. M.</given-names></name> & <name><surname>Jones</surname><given-names>P. L.</given-names></name><article-title>A brain-derived MeCP2 complex supports a role for MeCP2 in RNA processing</article-title>. <source>Biosci. Rep.</source><volume>31</volume>, <fpage>333</fpage>–<lpage>343</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21070191</pub-id></mixed-citation></ref><ref id="b37"><mixed-citation publication-type="journal"><name><surname>Ghosh</surname><given-names>R. P.</given-names></name>, <name><surname>Horowitz-Scherer</surname><given-names>R. A.</given-names></name>, <name><surname>Nikitina</surname><given-names>T.</given-names></name>, <name><surname>Shlyakhtenko</surname><given-names>L. S.</given-names></name> & <name><surname>Woodcock</surname><given-names>C. L.</given-names></name><article-title>MeCP2 binds cooperatively to its substrate and competes with histone H1 for chromatin binding sites</article-title>. <source>Mol. Cell Biol.</source><volume>30</volume>, <fpage>4656</fpage>–<lpage>4670</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20679481</pub-id></mixed-citation></ref><ref id="b38"><mixed-citation publication-type="journal"><name><surname>Shepard</surname><given-names>P. J.</given-names></name> & <name><surname>Hertel</surname><given-names>K. J.</given-names></name><article-title>The SR protein family</article-title>. <source>Genome Biol.</source><volume>10</volume>, <fpage>242</fpage> (<year>2009</year>).<pub-id pub-id-type="pmid">19857271</pub-id></mixed-citation></ref><ref id="b39"><mixed-citation publication-type="journal"><name><surname>Dichmann</surname><given-names>D. S.</given-names></name>, <name><surname>Walentek</surname><given-names>P.</given-names></name> & <name><surname>Harland</surname><given-names>R. M.</given-names></name><article-title>The alternative splicing regulator Tra2b is required for somitogenesis and regulates splicing of an inhibitory Wnt11b isoform</article-title>. <source>Cell Rep.</source><volume>10</volume>, <fpage>527</fpage>–<lpage>536</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25620705</pub-id></mixed-citation></ref><ref id="b40"><mixed-citation publication-type="journal"><name><surname>Best</surname><given-names>A.</given-names></name><etal></etal>. <article-title>Human Tra2 proteins jointly control a CHEK1 splicing switch among alternative and constitutive target exons</article-title>. <source>Nat. Commun.</source><volume>5</volume>, <fpage>4760</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25208576</pub-id></mixed-citation></ref><ref id="b41"><mixed-citation publication-type="journal"><name><surname>Hussein</surname><given-names>S. M. I.</given-names></name><etal></etal>. <article-title>Genome-wide characterization of the routes to pluripotency</article-title>. <source>Nature</source><volume>516</volume>, <fpage>198</fpage>–<lpage>206</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25503233</pub-id></mixed-citation></ref><ref id="b42"><mixed-citation publication-type="journal"><name><surname>Dvinge</surname><given-names>H.</given-names></name> & <name><surname>Bradley</surname><given-names>R.</given-names></name><article-title>Widespread intron retention diversifies most cancer transcriptomes</article-title>. <source>Genome Med.</source><volume>7</volume>, <fpage>45</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26113877</pub-id></mixed-citation></ref><ref id="b43"><mixed-citation publication-type="journal"><name><surname>Jung</surname><given-names>H.</given-names></name><etal></etal>. <article-title>Intron retention is a widespread mechanism of tumor-suppressor inactivation</article-title>. <source>Nat. Genet.</source><volume>47</volume>, <fpage>1242</fpage>–<lpage>1248</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26437032</pub-id></mixed-citation></ref><ref id="b44"><mixed-citation publication-type="journal"><name><surname>Ni</surname><given-names>T.</given-names></name><etal></etal>. <article-title>Global intron retention mediated gene regulation during CD4+ T cell activation</article-title>. <source>Nucleic Acids Res.</source><volume>44</volume>, <fpage>6817</fpage>–<lpage>6829</lpage> (<year>2016</year>).<pub-id pub-id-type="pmid">27369383</pub-id></mixed-citation></ref><ref id="b45"><mixed-citation publication-type="journal"><name><surname>Wen</surname><given-names>H.</given-names></name><etal></etal>. <article-title>ZMYND11 links histone H3.3K36me3 to transcription elongation and tumour suppression</article-title>. <source>Nature</source><volume>508</volume>, <fpage>263</fpage>–<lpage>268</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24590075</pub-id></mixed-citation></ref><ref id="b46"><mixed-citation publication-type="journal"><name><surname>Marina</surname><given-names>R. J.</given-names></name><etal></etal>. <article-title>TET-catalyzed oxidation of intragenic 5-methylcytosine regulates CTCF-dependent alternative splicing</article-title>. <source>EMBO J.</source><volume>35</volume>, <fpage>335</fpage>–<lpage>355</lpage> (<year>2016</year>).<pub-id pub-id-type="pmid">26711177</pub-id></mixed-citation></ref><ref id="b47"><mixed-citation publication-type="journal"><name><surname>Tekpli</surname><given-names>X.</given-names></name><etal></etal>. <article-title>Changes of 5-hydroxymethylcytosine distribution during myeloid and lymphoid differentiation of CD34+ cells</article-title>. <source>Epigenetics Chromatin</source><volume>9</volume>, <fpage>21</fpage> (<year>2016</year>).<pub-id pub-id-type="pmid">27252783</pub-id></mixed-citation></ref><ref id="b48"><mixed-citation publication-type="journal"><name><surname>Song</surname><given-names>C.-X.</given-names></name>, <name><surname>Yi</surname><given-names>C.</given-names></name> & <name><surname>He</surname><given-names>C.</given-names></name><article-title>Mapping recently identified nucleotide variants in the genome and transcriptome</article-title>. <source>Nat. Biotechnol.</source><volume>30</volume>, <fpage>1107</fpage>–<lpage>1116</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">23138310</pub-id></mixed-citation></ref><ref id="b49"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>X.</given-names></name><etal></etal>. <article-title>Gene body methylation can alter gene expression and is a therapeutic target in cancer</article-title>. <source>Cancer Cell</source><volume>26</volume>, <fpage>577</fpage>–<lpage>590</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25263941</pub-id></mixed-citation></ref><ref id="b50"><mixed-citation publication-type="journal"><name><surname>Buckstein</surname><given-names>R.</given-names></name>, <name><surname>Yee</surname><given-names>K.</given-names></name> & <name><surname>Wells</surname><given-names>R. A.</given-names></name><article-title>5-Azacytidine in myelodysplastic syndromes: a clinical practice guideline</article-title>. <source>Cancer Treat. Rev.</source><volume>37</volume>, <fpage>160</fpage>–<lpage>167</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">20591575</pub-id></mixed-citation></ref><ref id="b51"><mixed-citation publication-type="journal"><name><surname>Brown</surname><given-names>R.</given-names></name>, <name><surname>Curry</surname><given-names>E.</given-names></name>, <name><surname>Magnani</surname><given-names>L.</given-names></name>, <name><surname>Wilhelm-Benartzi</surname><given-names>C. S.</given-names></name> & <name><surname>Borley</surname><given-names>J.</given-names></name><article-title>Poised epigenetic states and acquired drug resistance in cancer</article-title>. <source>Nat. Rev. Cancer</source><volume>14</volume>, <fpage>747</fpage>–<lpage>753</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25253389</pub-id></mixed-citation></ref><ref id="b52"><mixed-citation publication-type="journal"><name><surname>Middleton</surname><given-names>R.</given-names></name><etal></etal>. <article-title>IRFinder: assessing the impact of intron retention on mammalian gene expression</article-title>. <source>Genome Biol.</source><volume>18</volume>, <fpage>51</fpage> (<year>2017</year>).<pub-id pub-id-type="pmid">28298237</pub-id></mixed-citation></ref><ref id="b53"><mixed-citation publication-type="journal"><name><surname>Wong</surname><given-names>J.</given-names></name><etal></etal>. <article-title>Identification of nuclear-enriched miRNAs during mouse granulopoiesis</article-title>. <source>J. Hematol. Oncol.</source><volume>7</volume>, <fpage>42</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24886830</pub-id></mixed-citation></ref><ref id="b54"><mixed-citation publication-type="journal"><name><surname>Mohammed</surname><given-names>H.</given-names></name><etal></etal>. <article-title>Endogenous purification reveals GREB1 as a key estrogen receptor regulatory factor</article-title>. <source>Cell Rep.</source><volume>3</volume>, <fpage>342</fpage>–<lpage>349</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23403292</pub-id></mixed-citation></ref><ref id="b55"><mixed-citation publication-type="journal"><name><surname>Challen</surname><given-names>G. A.</given-names></name><etal></etal>. <article-title>Dnmt3a and Dnmt3b have overlapping and distinct functions in hematopoietic stem cells</article-title>. <source>Cell Stem Cell</source><volume>15</volume>, <fpage>350</fpage>–<lpage>364</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25130491</pub-id></mixed-citation></ref><ref id="b56"><mixed-citation publication-type="journal"><name><surname>Jeong</surname><given-names>M.</given-names></name><etal></etal>. <article-title>Large conserved domains of low DNA methylation maintained by Dnmt3a</article-title>. <source>Nat. Genet.</source><volume>46</volume>, <fpage>17</fpage>–<lpage>23</lpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24270360</pub-id></mixed-citation></ref><ref id="b57"><mixed-citation publication-type="journal"><name><surname>Lund</surname><given-names>K.</given-names></name><etal></etal>. <article-title>DNMT inhibitors reverse a specific signature of aberrant promoter DNA methylation and associated gene silencing in AML</article-title>. <source>Genome Biol.</source><volume>15</volume>, <fpage>406</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">25315154</pub-id></mixed-citation></ref><ref id="b58"><mixed-citation publication-type="journal"><name><surname>Trapnell</surname><given-names>C.</given-names></name><etal></etal>. <article-title>Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks</article-title>. <source>Nat. Protoc.</source><volume>7</volume>, <fpage>562</fpage>–<lpage>578</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22383036</pub-id></mixed-citation></ref><ref id="b59"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>Y.</given-names></name><etal></etal>. <article-title>Model-based analysis of ChIP-Seq (MACS)</article-title>. <source>Genome Biol.</source><volume>9</volume>, <fpage>R137</fpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18798982</pub-id></mixed-citation></ref><ref id="b60"><mixed-citation publication-type="journal"><name><surname>Dobin</surname><given-names>A.</given-names></name><etal></etal>. <article-title>STAR: ultrafast universal RNA-seq aligner</article-title>. <source>Bioinformatics</source><volume>29</volume>, <fpage>15</fpage>–<lpage>21</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23104886</pub-id></mixed-citation></ref><ref id="b61"><mixed-citation publication-type="journal"><name><surname>Audic</surname><given-names>S.</given-names></name> & <name><surname>Claverie</surname><given-names>J. M.</given-names></name><article-title>The significance of digital gene expression profiles</article-title>. <source>Genome Res.</source><volume>7</volume>, <fpage>986</fpage>–<lpage>995</lpage> (<year>1997</year>).<pub-id pub-id-type="pmid">9331369</pub-id></mixed-citation></ref><ref id="b62"><mixed-citation publication-type="journal"><name><surname>Huang</surname><given-names>D. W.</given-names></name>, <name><surname>Sherman</surname><given-names>B. T.</given-names></name> & <name><surname>Lempicki</surname><given-names>R. A.</given-names></name><article-title>Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources</article-title>. <source>Nat. Protoc.</source><volume>4</volume>, <fpage>44</fpage>–<lpage>57</lpage> (<year>2008</year>).</mixed-citation></ref></ref-list><fn-group><fn fn-type="COI-statement"><p>The authors declare no competing financial interests.</p></fn><fn><p><bold>Author contributions</bold> J.J.-L.W., D.G. and T.V.N. designed the experiments and analysed data. D.G., R.M. and W.R performed bioinformatics analyses and the IRFinder algorithm development. J.J.-L.W., T.V.N., C.-T.K., N.P., A.T. and R.N. performed cellular and molecular biology experiments. M.v.G., T.V.N. and J.H. performed MS data analyses. J.J.-L.W. and J.E.J.R. conceived and designed the overall project, and supervised its execution. J.J.W., D.G., T.V.N., W.R. and J.E.J.R. wrote the manuscript. J.J.W., D.G. and T.V.N. contributed equally to the paper.</p></fn></fn-group></back><floats-group><fig id="f1"><label>Figure 1</label><caption><title>IR in myeloid cells is associated with lower DNA methylation levels near splice junctions and within retained introns.</title><p>(<bold>a</bold>) Average CpG methylation (CpG meth.) levels spanning ±100 bp from the 5′ and 3′ splice junctions, and from the middle of introns, in the DNA encoding non-retained (blue) and retained (red) introns in promyelocytes and granulocytes combined. Data were extended to −20 Mb and +20 Mb from the 5′ and 3′ splice junction respectively. <italic>t</italic>-test was used to determine significance at <italic>P</italic><0.05. (<bold>b</bold>) Examples from the IGV genome browser comparing IR levels in <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3 and <italic>S100a8</italic> intron 2 in promyelocytes (Prom.) and granulocytes (Gran.). (<bold>c</bold>) Average DNA methylation level at each CpG across <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3 and <italic>S100a8</italic> intron 2 and their flanking exons in promyelocytes and granulocytes as assessed using whole-genome bisulfite sequencing. A ±100 bp region from each splice junction, with reduced DNA methylation in granulocytes, is shaded in pink. (<bold>d</bold>) Confirmation of the results in <bold>c</bold> using clonal bisulfite sequencing. Each row represents a single-cloned PCR amplicon aligned to the exon–intron map shown above. Lollipops on exon–intron map indicate positions of CpG sites. Circles denote methylated (black) and unmethylated (white) CpGs. NS, not significant; Avg. meth., average DNA methylation levels.</p></caption><graphic xlink:href="ncomms15134-f1"></graphic></fig><fig id="f2"><label>Figure 2</label><caption><title>IR increases with reduced DNA methylation and increases following DNA methylation inhibition.</title><p>(<bold>a</bold>) Average CpG methylation levels spanning ±100 bp from the 5′ and 3′ splice junctions in the DNA encoding non-retained (blue) and retained (red) introns in primary mouse fibroblasts, mouse reprogrammed fibroblasts, human cell lines (H1, H9, HCT116 and IMR90) and primary human neuron progenitors. Data were extended to −20 Mb and +20 Mb from the 5′ and 3′ splice junctions, respectively. Paired <italic>t</italic>-tests were used to determine significance at <italic>P</italic><0.05. (<bold>b</bold>) Association between the IR ratio fold changes and their significance following double knockout of <italic>Dnmt3a</italic> and <italic>3b</italic> in primary mouse B cells and HSCs. Binomial test was used to determine the significance of bias between increased and decreased IR at <italic>P</italic><0.05. (<bold>c</bold>–<bold>f</bold>) Upper panel: clonal bisulfite sequencing displaying DNA methylation changes near splice junctions of <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3, <italic>Spata13</italic> intron 8 and <italic>S100a8</italic> intron 2 in MPRO cells following treatment with 3 or 7 μM 5-aza-2′deoxycytidine (5-Aza). Each row represents a single-cloned PCR amplicon aligned to the exon–intron map shown. Each lollipop on the map represents a CpG site. Black and white circles denote methylated and unmethylated CpGs respectively. Lower panel: IR levels determined by qRT-PCR for specific introns indicated in MPRO cells, with and without exposure to 5-Aza at 3 or 7 μM, and in the absence (–) or presence (+) of NMD inhibition via caffeine (CAF) treatment. Two-tailed <italic>t</italic>-test was used to determine significance at <italic>P</italic><0.05. Bars display mean±s.e.m. Independent qRT-PCR experiments were performed three times in triplicate (<italic>n</italic>=3). NS, not significant; Avg. meth., average DNA methylation levels.</p></caption><graphic xlink:href="ncomms15134-f2"></graphic></fig><fig id="f3"><label>Figure 3</label><caption><title>IR is associated with reduced MeCP2 binding near splice junctions.</title><p>(<bold>a</bold>) Average MeCP2 occupancy normalized against input as measured by ChIP-seq in promyelocytes and granulocytes combined. MeCP2 occupancy is displayed for regions spanning ±100 bp from the 5′ and 3′ splice junctions, and the middle of introns, in the DNA encoding non-retained (blue) and retained introns (red). Data were extended to −20 Mb and +20 Mb from the 5′ and 3′ splice junction respectively. <italic>t</italic>-tests were used to determine significance at <italic>P</italic><0.05. (<bold>b</bold>) Occupancy of MeCP2 near splice junctions of the DNA encoding introns with increased IR in granulocytes and controls measured by ChIP-qPCR. (<bold>c</bold>) ChIP-qPCR-measured occupancy of MeCP2 near splice junctions of the DNA encoding differentially retained introns of <italic>Lmnb1</italic> and <italic>Ngp</italic> pre- and post- 5-Aza-2′deoxycytidine treatment (5-Aza) in MPRO cells. <italic>Atf4</italic> intron 2 is included as a negative control. (<bold>d</bold>) Association between IR ratio fold changes and their significance following <italic>Mecp2</italic> knockdown in IMR90 cells or <italic>Mecp2</italic>- knockout in primary mouse cerebellum. Binomial test was used to determine the significance of bias between increased and decreased IR at <italic>P</italic><0.05. For each ChIP-qPCR experiment, a two-tailed Student's <italic>t</italic>-test was used to determine significance, denoted by *<italic>P</italic><0.05, **<italic>P</italic><0.001 and ns, not significant. Bars display mean±s.e.m. Independent ChIP-qPCR experiments were performed three to five times in triplicate.</p></caption><graphic xlink:href="ncomms15134-f3"></graphic></fig><fig id="f4"><label>Figure 4</label><caption><title>IR is associated with increased RNA Pol II occupancy that anti-correlates with MeCP2 density.</title><p>(<bold>a</bold>) Occupancies of MeCP2 and (<bold>b</bold>) Pol II Ser2p measured by ChIP-qPCR for retained introns and their flanking exons of <italic>Lmnb1</italic>, <italic>Ngp</italic> and <italic>S100a8</italic> in promyelocytes and granulocytes. (<bold>c</bold>) ChIP-qPCR analyses of MeCP2 and (<bold>d</bold>) Pol II Ser2p occupancies measured for non-retained introns: <italic>Lmnb1</italic> intron 3, <italic>Smarcd1</italic> intron 9, <italic>Tbp</italic> intron 5 and their flanking exons. P1–P8 indicate positions of PCR amplicons within each region. For each ChIP-qPCR experiment, a two-tailed Student's <italic>t</italic>-test was used to determine significance, denoted by *<italic>P</italic><0.05, **<italic>P</italic><0.001 and NS, not significant. Dots and bars display mean±s.e.m. Independent ChIP-qPCR experiments were performed three to five times in triplicate.</p></caption><graphic xlink:href="ncomms15134-f4"></graphic></fig><fig id="f5"><label>Figure 5</label><caption><title>Reduced MeCP2-mediated recruitment of the splicing factor Tra2b promotes IR.</title><p>(<bold>a</bold>) Schematic of the RIME assay to identify MeCP2 binding partners in normal promyelocytes (Prom.) and granulocytes (Gran.). (<bold>b</bold>) Venn diagram showing proteins that bind to MeCP2 in promyelocytes and/or granulocytes. MeCP2 bound splicing factors enriched in promyelocytes are listed. (<bold>c</bold>) Occupancies of MeCP2 and (<bold>d</bold>) Tra2b in the coding region of <italic>Malat1</italic> and <italic>Tra2a</italic> in promyelocytes and granulocytes measured by RNA-IP followed by qRT-PCR. (<bold>e</bold>) Occupancies of MeCP2 and (<bold>f</bold>) Tra2b near splice junctions of <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3 and <italic>S100a8</italic> intron 2 in normal promyelocytes compared to granulocytes by RNA-IP. (<bold>g</bold>) Occupancies of MeCP2 and (<bold>h</bold>) Tra2b near splice junctions of constitutively spliced <italic>Smarcd1</italic> intron 9 and <italic>TBP</italic> intron 5 in normal promyelocytes compared to granulocytes by RNA-IP. (<bold>i</bold>) Occupancies of MeCP2 and Tra2b near splice junctions of <italic>Lmnb1</italic> intron 5, <italic>Ngp</italic> intron 3 and <italic>S100a8</italic> intron 2, <italic>Smarcd1</italic> intron 9 and <italic>TBP</italic> intron 5 in normal promyelocytes compared to granulocytes by RNA-IP-reIP. Levels of occupancies are shown relative to background inferred by sequential Tra2b-IgG RNA-IP controls. A two-tailed <italic>t</italic>-test was used to determine significance, denoted by *<italic>P</italic><0.05, **<italic>P</italic><0.001 and NS, not significant. Dots and bars display mean±s.e.m. Independent qRT-PCR experiments were performed three to five times in triplicate.</p></caption><graphic xlink:href="ncomms15134-f5"></graphic></fig><fig id="f6"><label>Figure 6</label><caption><title>Proposed model of IR regulation.</title><p>IR occurs via less efficient splicing factor (e.g. Tra2b) recruitment consequent to reduced DNA-methylation-mediated MeCP2 binding.</p></caption><graphic xlink:href="ncomms15134-f6"></graphic></fig></floats-group></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Asie/explor/AustralieFrV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000790 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000790 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Asie
|area= AustralieFrV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:5424149
|texte= Intron retention is regulated by altered MeCP2-mediated splicing factor recruitment
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:28480880" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a AustralieFrV1
| This area was generated with Dilib version V0.6.33. |  |