SNPsplit: Allele-specific Splitting of Alignments Between Genomes with Known SNP Genotypes

Overview

Journal F1000Res

Specialties Biomedical Engineering
Science

Date 2016 Aug 19

PMID 27429743

Citations 106

Authors

Felix Krueger

Simon R Andrews

Affiliations

Soon will be listed here.

Abstract

Sequencing reads overlapping polymorphic sites in diploid mammalian genomes may be assigned to one allele or the other. This holds the potential to detect gene expression, chromatin modifications, DNA methylation or nuclear interactions in an allele-specific fashion. SNPsplit is an allele-specific alignment sorter designed to read files in SAM/BAM format and determine the allelic origin of reads or read-pairs that cover known single nucleotide polymorphic (SNP) positions. For this to work libraries must have been aligned to a genome in which all known SNP positions were masked with the ambiguity base 'N' and aligned using a suitable mapping program such as Bowtie2, TopHat, STAR, HISAT2, HiCUP or Bismark. SNPsplit also provides an automated solution to generate N-masked reference genomes for hybrid mouse strains based on the variant call information provided by the Mouse Genomes Project. The unique ability of SNPsplit to work with various different kinds of sequencing data including RNA-Seq, ChIP-Seq, Bisulfite-Seq or Hi-C opens new avenues for the integrative exploration of allele-specific data.

Citing Articles

Reconstruction of diploid higher-order human 3D genome interactions from noisy Pore-C data using Dip3D.

Chen Y, Lin Z, Wang S, Wu B, Niu L, Zhong J Nat Struct Mol Biol. 2025; .

PMID: 40038455 DOI: 10.1038/s41594-025-01512-w.

Transcription can be sufficient, but is not necessary, to advance replication timing.

Vouzas A, Sasaki T, Rivera-Mulia J, Turner J, Brown A, Alexander K bioRxiv. 2025; .

PMID: 39975371 PMC: 11838563. DOI: 10.1101/2025.02.04.636516.

Multi-omics analysis in primary T cells elucidates mechanisms behind disease-associated genetic loci.

Shi C, Zhao D, Butler J, Frantzeskos A, Rossi S, Ding J Genome Biol. 2025; 26(1):26.

PMID: 39930543 PMC: 11808986. DOI: 10.1186/s13059-025-03492-y.

Temporal and regional X-linked gene reactivation in the mouse germline reveals site-specific retention of epigenetic silencing.

Roidor C, Syx L, Beyne E, Raynaud P, Zielinski D, Teissandier A Nat Struct Mol Biol. 2025; .

PMID: 39838109 DOI: 10.1038/s41594-024-01469-2.

Acquired sperm hypomethylation by gestational arsenic exposure is re-established in both the paternal and maternal genomes of post-epigenetic reprogramming embryos.

Nohara K, Suzuki T, Okamura K, Kawai T, Nakabayashi K Epigenetics Chromatin. 2025; 18(1):4.

PMID: 39815295 PMC: 11737231. DOI: 10.1186/s13072-025-00569-7.

References

Rao S, Huntley M, Durand N, Stamenova E, Bochkov I, Robinson J . A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014; 159(7):1665-80. PMC: 5635824. DOI: 10.1016/j.cell.2014.11.021. View

Li H, Durbin R . Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010; 26(5):589-95. PMC: 2828108. DOI: 10.1093/bioinformatics/btp698. View

Degner J, Marioni J, Pai A, Pickrell J, Nkadori E, Gilad Y . Effect of read-mapping biases on detecting allele-specific expression from RNA-sequencing data. Bioinformatics. 2009; 25(24):3207-12. PMC: 2788925. DOI: 10.1093/bioinformatics/btp579. View

Langmead B, Salzberg S . Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9(4):357-9. PMC: 3322381. DOI: 10.1038/nmeth.1923. View

Castel S, Levy-Moonshine A, Mohammadi P, Banks E, Lappalainen T . Tools and best practices for data processing in allelic expression analysis. Genome Biol. 2015; 16:195. PMC: 4574606. DOI: 10.1186/s13059-015-0762-6. View

Rozowsky J, Abyzov A, Wang J, Alves P, Raha D, Harmanci A . AlleleSeq: analysis of allele-specific expression and binding in a network framework. Mol Syst Biol. 2011; 7:522. PMC: 3208341. DOI: 10.1038/msb.2011.54. View

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N . The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078-9. PMC: 2723002. DOI: 10.1093/bioinformatics/btp352. View

Servant N, Varoquaux N, Lajoie B, Viara E, Chen C, Vert J . HiC-Pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015; 16:259. PMC: 4665391. DOI: 10.1186/s13059-015-0831-x. View

Keane T, Goodstadt L, Danecek P, White M, Wong K, Yalcin B . Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011; 477(7364):289-94. PMC: 3276836. DOI: 10.1038/nature10413. View

10.

Strogantsev R, Krueger F, Yamazawa K, Shi H, Gould P, Goldman-Roberts M . Allele-specific binding of ZFP57 in the epigenetic regulation of imprinted and non-imprinted monoallelic expression. Genome Biol. 2015; 16:112. PMC: 4491874. DOI: 10.1186/s13059-015-0672-7. View

11.

Lieberman-Aiden E, van Berkum N, Williams L, Imakaev M, Ragoczy T, Telling A . Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009; 326(5950):289-93. PMC: 2858594. DOI: 10.1126/science.1181369. View

12.

Xie W, Barr C, Kim A, Yue F, Lee A, Eubanks J . Base-resolution analyses of sequence and parent-of-origin dependent DNA methylation in the mouse genome. Cell. 2012; 148(4):816-31. PMC: 3343639. DOI: 10.1016/j.cell.2011.12.035. View

13.

van de Geijn B, McVicker G, Gilad Y, Pritchard J . WASP: allele-specific software for robust molecular quantitative trait locus discovery. Nat Methods. 2015; 12(11):1061-3. PMC: 4626402. DOI: 10.1038/nmeth.3582. View

14.

Krueger F, Andrews S . SNPsplit: Allele-specific splitting of alignments between genomes with known SNP genotypes. F1000Res. 2016; 5:1479. PMC: 4934512. DOI: 10.12688/f1000research.9037.2. View

15.

Selvaraj S, Dixon J, Bansal V, Ren B . Whole-genome haplotype reconstruction using proximity-ligation and shotgun sequencing. Nat Biotechnol. 2013; 31(12):1111-8. PMC: 4180835. DOI: 10.1038/nbt.2728. View

16.

Rivas M, Pirinen M, Conrad D, Lek M, Tsang E, Karczewski K . Human genomics. Effect of predicted protein-truncating genetic variants on the human transcriptome. Science. 2015; 348(6235):666-9. PMC: 4537935. DOI: 10.1126/science.1261877. View

17.

Dobin A, Davis C, Schlesinger F, Drenkow J, Zaleski C, Jha S . STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2012; 29(1):15-21. PMC: 3530905. DOI: 10.1093/bioinformatics/bts635. View

18.

Kim D, Langmead B, Salzberg S . HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12(4):357-60. PMC: 4655817. DOI: 10.1038/nmeth.3317. View

19.

Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg S . TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013; 14(4):R36. PMC: 4053844. DOI: 10.1186/gb-2013-14-4-r36. View

20.

Krueger F, Andrews S . Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011; 27(11):1571-2. PMC: 3102221. DOI: 10.1093/bioinformatics/btr167. View