Nucleosome positioning is unaltered at MLH1 splice site mutations in cells derived from Lynch syndrome patients
© Sloane et al.; licensee BioMed Central Ltd. 2014
Received: 4 July 2014
Accepted: 28 November 2014
Published: 13 December 2014
Splicing is more efficient when coupled with transcription and it has been proposed that nucleosomes enriched in exons are important for splice site recognition. Lynch syndrome is a familial cancer syndrome that can be caused by the autosomal dominant inheritance of splice site mutations in the MutL homolog 1 (MLH1) gene. To better understand the role of nucleosomes in splicing, we used MLH1 splice site mutations in Lynch syndrome cases as a model to investigate if abnormal splicing was associated with altered nucleosome positioning at exon-intron boundaries.
Nucleosome Occupancy and Methylome sequencing (NOMe-seq) was used to determine the allele-specific positioning of nucleosomes around heterozygous splice site mutations in lymphoblastoid cells lines (LCLs) derived from six Lynch syndrome patients. These mutations were previously shown to cause exon skipping in five of the six patients. Allele-specific high-resolution nucleosome mapping across exons and exon-intron boundaries revealed high levels of nucleosomes across all regions examined. Alleles containing donor or acceptor splice site mutations showed no consistent alteration in nucleosome positioning or occupancy.
Nucleosomes were enriched at MLH1 exons in LCLs derived from Lynch syndrome patients, and in this model system the positioning of nucleosomes was unaltered at exon-intron boundaries containing splice site mutations. Thus, these splice site mutations alone do not significantly change the local organisation of nucleosomes.
Splicing removes introns from a longer pre-cursor RNA molecule to produce a final processed mRNA. Splice donor and acceptor sites located at the junctions between exons and introns contain conserved sequence elements that are necessary for splicing . Although in vitro splicing can occur, splicing is more efficient when coupled with transcription , suggesting that factors additional to the nucleic acid sequence are important for splice site recognition. It has been proposed that the positioning of nucleosomes at exons aids in splice-site recognition ,.
Nucleosome positioning describes the precise location of a given nucleosome, whereas nucleosome occupancy refers to the proportion of molecules bearing a nucleosome at a specific location, at any given instant . It is hypothesised that nucleosomes positioned within exons, especially those with weak splice sites, cause RNA Polymerase II (RNAPII) to pause, enabling an interaction with the spliceosome and more efficient splicing . In support of this, nucleosome occupancy is enriched across exons ,, which have a high GC content that favours nucleosome assembly ,, and the average size of an exon within the body of a human gene is 151 bp, that is similar in length to nucleosomal DNA (approximately 147 bp) . Splicing factors associate with the C-terminal domain of RNA Polymerase II (RNAPII) , while the histone modification H3K36me3 is enriched in exons ,, where it may act as a scaffold to recruit splicing factors . Taken together this indicates there is interplay between nucleosomes, the splicing machinery and the DNA sequence.
Lynch syndrome is an autosomal dominant familial cancer syndrome characterised by early onset colorectal, endometrial and other cancers ,. It is most commonly caused by the inheritance of heterozygous loss-of-function mutations, including splice site mutations, in the DNA mismatch repair (MMR) genes MLH1 and MSH2. In this study we used MLH1 splice site mutations in Lynch syndrome as a model to better understand the role of nucleosomes in splicing. We investigated cells from Lynch syndrome patients with splice site mutations to determine whether splicing aberrations were associated with altered nucleosome positioning at exon-intron boundaries.
Materials and methods
Identification of MLH1 splice site mutations
Molecular features of MLH1 splice site mutations in six individuals with Lynch syndrome
Mutation location (GRC37, Feb 2009)
Affected splice site
Nucleosome occlusion of splice site
Proportion of alleles
c.588 + 1G > T
Chr. 3: 37,053,354
Exon 7 skipping (r.546_588del)
First bp of intron 7
c.589-2A > G
Chr. 3: 37,053,500
4 bp deletion in cDNA (r.589_592del)
2nd last bp of intron 7
c.790 + 1G > T
Chr. 3: 37,056,036
1st bp of intron 9
c.791-1G > C
Chr. 3: 37,058,996
Exon 10 skipping (r.791_884del)
Last bp of intron 9
c.884G > A
Chr. 3: 37,059,090
Exon 10 skipping (r.791_884del)
Last bp of exon 10
c.1559-2A > T
Chr. 3: 37,081,675
Two aberrant transcripts (Exon 14 skipping or exon 14 and 15 skipping – r.[1559_1667del, 1559_1731del])
Second last bp of intron 13
LCLs were established from patient blood by transformation with Epstein-Barr Virus (as described previously) , and cultured in RPMI with 10% fetal bovine serum (Gibco, Life Technologies) at 37°C in 5% CO2. Lymphocytes for transformation from five patients (c.588 + 1G > T, c.589-2A > G, c.791-1G > C, c.884G > A and c.1559-2A > T) were obtained from the Australasian Colorectal Cancer Family Registry (ACCFR) . Lymphocytes from the patient with the c.790 + 1G > T mutation were obtained from the MCO collection .
NOMe-seq primer sequences and amplification conditions
c.588 + 1G > T
F/R = 54°C
F/RN = 56°C
c.589-2A > G
F/R = 51°C
F/RN = 52°C
c.790 + 1G > T
F/R = 52°C
F/RN = 52°C
c.791-1G > C
F/R = 55°C
F/RN = 56°C
c.884G > A
F/R = 58°C
F/RN = 58°C
c.1559-2A > T
F/R = 58°C
F/RN = 57°C
F/R = 56°C
Separation of alleles by presence or absence of the splice site mutation showed no overall difference in either nucleosome positioning or occupancy between wild-type and mutant alleles (Figure 1 and Table 1). In the six splice site mutations analysed, three were located within donor splice sites and three were located within acceptor splice sites (Table 1). Though some difference in precise positioning between mutant and wild-type alleles was observed across some exons (for example, exons 7 and 10, Figure 1B) this was confined to a small subset of molecules, rather than the majority of molecules as would be expected if the mutations affected nucleosome positioning. Our study therefore shows that nucleosome occupancy is unaltered at exon-intron boundaries containing these splice site mutations.
In this study, we show that nucleosome positioning is unaltered at exon-intron boundaries containing splice site mutations. This finding suggests that the positioning of nucleosomes at exon-intron boundaries containing splice site mutations does not play a role in the splicing defect.
NOMe-seq has been used to measure nucleosome occupancy at gene promoters, CTCF binding sites  and regulatory elements  but this is the first study that has utilised the technique to measure nucleosome occupancy at splice sites. NOMe-seq previously showed that the active CpG island (CGI) promoter of HSPA5 has a nucleosome-depleted region (NDR) upstream of the transcriptional start site . This region was included as a control in each assay, and in all cases the NDR was present. This indicated that the absence of GpC methylation observed in MLH1 was caused by nucleosome occlusion, not an artefact caused by a defective M.Cvi PI enzyme. In addition, CpG methylation was absent from the HSPA5 NDR which is a feature of active CGI promoters . Finally, gene body methylation is a feature of human cells , and the majority of CpG sites within the MLH1 exons and introns were methylated. Together, these controls indicated that the NOMe-seq assays captured the true state of nucleosome positioning at MLH1 splice sites.
Alternative splicing of MLH1 occurs in a range of normal tissues , but several pieces of evidence indicate that the transcripts observed in our study are caused by the mutation, rather than being normal splicing events. First, the splicing aberrations were predicted with bioinformatic splicing software , and second, although one of the aberrant transcripts (Δ10) has been reported among naturally occurring MLH1 splice transcripts, the clinical data associated with these variants are also indicative of pathogenicity as demonstrated by the results of previous multifactorial likelihood analyses .
One strength of our study was the ability to determine the allele-specific position of nucleosomes by incorporating a heterozygous single nucleotide variant into each NOMe-seq assay. Small changes in positioning were observed but these molecules constituted approximately the same small proportion of total molecules on both alleles. Subtle changes in nucleosome positioning have been considered important in regulating the expression of cell cycle-dependent genes  and enhancer accessibility , but in those studies the subtle changes were observed on a significant number of molecules that resulted in an overall change in nucleosome occupancy. There was no consistent alteration in nucleosome occupancy between the wild-type and mutant alleles at a variety of donor and acceptor sites, which would be expected if nucleosome positioning at these sites played a role in mediating the effects of splice site mutations. Together our data shows that alleles containing a splice site mutation show no differences in nucleosome occupancy to wild-type alleles.
A limitation of the approach was that in most of the assays (assays 2 to 6) it was only possible to measure nucleosome positions at either the donor site or the acceptor site. Bisulfite conversion causes fragmentation of genomic DNA and it is technically challenging to amplify fragments greater than 500 bp in length . The donor and acceptor sites were separated by more than 2 kb of intronic sequence, making it impossible to concurrently determine the allele-specific position of nucleosomes at the donor and acceptor site on the same molecule. Although an independent NOMe-seq assay could be performed at the other donor or acceptor site, it would provide no information on allele-specificity. The data from one assay (c.588 + 1 G > T), however, in which both splice sites were present in the one amplicon, showed no significant change in nucleosome occupancy at the donor or acceptor site in intron 7.
Previous studies have utilised genome-wide datasets of micrococcal nuclease (MNase) digested chromatin to investigate nucleosome occupancy ,. Here we used NOMe-seq to map the position of nucleosomes relative to exons and splice sites in a single gene at single molecule resolution. This approach provides the most accurate possible measurement of nucleosome positioning and enables allele-specific mapping of nucleosomes. In agreement with previous reports in humans ,, Caenorhabditis elegans, and Drosophila melanogaster, we observed high nucleosome occupancy at exons. Thus, the findings from genome-wide studies were supported by our independent, single-molecule approach at the MLH1 gene.
A key finding of our study was that nucleosome positioning was not significantly affected by mutations at the MLH1 splice sites investigated. Recent work with an in vitro IKBKAP mini-gene system found that alternative splicing changed chromatin organisation, with splice site strength and factors needed for splicing, such as U1 snSNP, playing a role in the regulation of nucleosome occupancy in exons . This and an earlier study  indicate that splice site strength is an important determinant of nucleosome occupancy in exons. Exons with stronger polypyrimidine tracts (PPT; one of the conserved DNA elements located at 3′ splice sites) have increased nucleosome occupancy compared to immediately adjacent introns, suggesting that nucleosomes act as a barrier that help to define the physical location of the 3′ splice site . Although we investigated an equal number of donor and acceptor splice site mutations, nucleosome positioning was not significantly different between exons and the immediate 5′ and 3′ intronic sequence on the wild-type or mutant allele. The discrepancy between earlier studies and ours may in part relate to the strength of the MLH1 splice sites investigated. Alternatively, changes to chromatin organisation induced by splice mutations may be different to those caused by normal splicing events. Co-transcriptional splicing involves a complex interplay between RNAPII, the spliceosome  and histone modifications , and alterations in these factors, in addition to DNA sequence changes may be needed to alter nucleosome occupancy at exons.
In summary this study used NOMe-seq to determine the allele-specific position of nucleosomes in relation to MLH1 splice site mutations. We conclude that splice site mutations that cause aberrant splicing of MLH1 do not alone significantly affect local nucleosome positioning in LCLs from Lynch syndrome patients.
RLW is supported by Cancer Council NSW (RG 13-07) and Cancer Australia. We would like to thank the Australasian Colorectal Cancer Family Registry (ACCFR) for providing lymphocytes for transformation.
- Black DL: Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem. 2003, 72: 291-336. 10.1146/annurev.biochem.72.121801.161720.View ArticlePubMedGoogle Scholar
- Das R, Dufu K, Romney B, Feldt M, Elenko M, Reed R: Functional coupling of RNAP II transcription to spliceosome assembly. Genes Dev. 2006, 20: 1100-1109. 10.1101/gad.1397406.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwartz S, Meshorer E, Ast G: Chromatin organization marks exon-intron structure. Nat Struct Mol Biol. 2009, 16: 990-995. 10.1038/nsmb.1659.View ArticlePubMedGoogle Scholar
- Tilgner H, Nikolaou C, Althammer S, Sammeth M, Beato M, Valcárcel J, Guigó R: Nucleosome positioning as a determinant of exon recognition. Nat Struct Mol Biol. 2009, 16: 996-1001. 10.1038/nsmb.1658.View ArticlePubMedGoogle Scholar
- Berrozpe G, Bryant GO, Warpinski K, Ptashne M: Regulation of a mammalian gene bearing a CpG island promoter and a distal enhancer. Cell Rep. 2013, 4: 445-453. 10.1016/j.celrep.2013.07.001.View ArticlePubMedGoogle Scholar
- Field Y, Kaplan N, Fondufe-Mittendorf Y, Moore IK, Sharon E, Lubling Y, Widom J, Segal E: Distinct modes of regulation by chromatin encoded through nucleosome positioning signals. PLoS Comput Biol. 2008, 4: e1000216-10.1371/journal.pcbi.1000216.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Burge CB: Splicing regulation: from a parts list of regulatory elements to an integrated splicing code. RNA. 2008, 14: 802-813. 10.1261/rna.876308.PubMed CentralView ArticlePubMedGoogle Scholar
- de Almeida SF, Carmo-Fonseca M: The CTD role in cotranscriptional RNA processing and surveillance. FEBS Lett. 2008, 582: 1971-1976. 10.1016/j.febslet.2008.04.019.View ArticlePubMedGoogle Scholar
- Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41: 376-381. 10.1038/ng.322.PubMed CentralView ArticlePubMedGoogle Scholar
- Watson P, Lynch HT: Extracolonic cancer in hereditary nonpolyposis colorectal cancer. Cancer. 1993, 71: 677-685. 10.1002/1097-0142(19930201)71:3<677::AID-CNCR2820710305>3.0.CO;2-#.View ArticlePubMedGoogle Scholar
- Lynch HT, de la Chapelle A: Genetic susceptibility to non-polyposis colorectal cancer. J Med Genet. 1999, 36: 801-818.PubMed CentralPubMedGoogle Scholar
- Thompson BA, Spurdle AB, Plazzer JP, Greenblatt MS, Akagi K, Al-Mulla F, Bapat B, Bernstein I, Capellá G, den Dunnen JT, du Sart D, Fabre A, Farrell MP, Farrington SM, Frayling IM, Frebourg T, Goldgar DE, Heinen CD, Holinski-Feder E, Kohonen-Corish M, Robinson KL, Leung SY, Martins A, Moller P, Morak M, Nystrom M, Peltomaki P, Pineda M, Qi M, Ramesar R: Application of a 5-tiered scheme for standardized classification of 2,360 unique mismatch repair gene variants in the InSiGHT locus-specific database. Nat Genet. 2014, 46: 107-115. 10.1038/ng.2854.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson BA, Goldgar DE, Paterson C, Clendenning M, Walters R, Arnold S, Parsons MT, Michael DW, Gallinger S, Haile RW, Hopper JL, Jenkins MA, Lemarchand L, Lindor NM, Newcomb PA, Thibodeau SN, Colon Cancer Family R, Young JP, Buchanan DD, Tavtigian SV, Spurdle AB: A multifactorial likelihood model for MMR gene variant classification incorporating probabilities based on sequence bioinformatics and tumor characteristics: a report from the Colon Cancer Family Registry. Hum Mutat. 2013, 34: 200-209. 10.1002/humu.22213.PubMed CentralView ArticlePubMedGoogle Scholar
- Australasian Colorectal Cancer Family Registry [http://www.pedigree.org.au/pedigree-studies/accfr.aspx]
- MCO study tumour collection [http://researchdata.ands.org.au/mco-study-tumour-collection]
- Wolff EM, Byun HM, Han HF, Sharma S, Nichols PW, Siegmund KD, Yang AS, Jones PA, Liang G: Hypomethylation of a LINE-1 promoter activates an alternate transcript of the MET oncogene in bladders with cancer. PLoS Genet. 2010, 6: e1000917-10.1371/journal.pgen.1000917.PubMed CentralView ArticlePubMedGoogle Scholar
- Hesson LB, Ward RL: Discrimination of pseudogene and parental gene DNA methylation using allelic bisulfite sequencing. Methods Mol Biol. 2014, 1167: 265-274. 10.1007/978-1-4939-0835-6_18.View ArticlePubMedGoogle Scholar
- Kelly TK, Miranda TB, Liang G, Berman BP, Lin JC, Tanay A, Jones PAPA: H2A.Z maintenance during mitosis reveals nucleosome shifting on mitotically silenced genes. Mol Cell. 2010, 39: 901-911. 10.1016/j.molcel.2010.08.026.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelly TK, Liu Y, Lay FD, Liang G, Berman BP, Jones PA: Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules. Genome Res. 2012, 22: 2497-2506. 10.1101/gr.143008.112.PubMed CentralView ArticlePubMedGoogle Scholar
- Taberlay PC, Statham AL, Kelly TK, Clark SJ, Jones PA: Reconfiguration of nucleosome depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer.Genome Res 2014, Epub ahead of print.,Google Scholar
- Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25: 1010-1022. 10.1101/gad.2037511.PubMed CentralView ArticlePubMedGoogle Scholar
- Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009, 27: 361-368. 10.1038/nbt.1533.PubMed CentralView ArticlePubMedGoogle Scholar
- Chodavarapu RK, Feng S, Bernatavichute YV, Chen PY, Stroud H, Yu Y, Hetzel JA, Kuo F, Kim J, Cokus SJ, Casero D, Bernal M, Huijser P, Clark AT, Krämer U, Merchant SS, Zhang X, Jacobsen SE, Pellegrini M: Relationship between nucleosome positioning and DNA methylation. Nature. 2010, 466: 388-392. 10.1038/nature09147.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson BA, Martins A, Spurdle AB: A review of mismatch repair gene transcripts: issues for interpretation of mRNA splicing assays.Clin Genet 2014. Epub ahead of print.,Google Scholar
- Andreu-Vieyra C, Lai J, Berman BP, Frenkel B, Jia L, Jones PA, Coetzee GA: Dynamic nucleosome-depleted regions at androgen receptor enhancers in the absence of ligand in prostate cancer cells. Mol Cell Bio. 2011, 31: 4648-4662. 10.1128/MCB.05934-11.View ArticleGoogle Scholar
- Grunau C, Clark SJ, Rosenthal A: Bisulfite genomic sequencing: systematic investigation of critical experimental parameters. Nucleic Acids Res. 2001, 29: E65-5. 10.1093/nar/29.13.e65.PubMed CentralView ArticlePubMedGoogle Scholar
- Keren-Shaul H, Lev-Maor G, Ast G: Pre-mRNA splicing is a determinant of nucleosome organization. PLoS One. 2013, 8: e53506-10.1371/journal.pone.0053506.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.