A Gibbs sampling strategy applied to the mapping of ambiguous short-sequence tags

Abstract

Motivation: Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is widely used in biological research. ChIP-seq experiments yield many ambiguous tags that can be mapped with equal probability to multiple genomic sites. Such ambiguous tags are typically eliminated from consideration resulting in a potential loss of important biological information.Results: We have developed a Gibbs sampling-based algorithm for the genomic mapping of ambiguous sequence tags. Our algorithm relies on the local genomic tag context to guide the mapping of ambiguous tags. The Gibbs sampling procedure we use simultaneously maps ambiguous tags and updates the probabilities used to infer correct tag map positions. We show that our algorithm is able to correctly map more ambiguous tags than existing mapping methods. Our approach is also able to uncover mapped genomic sites from highly repetitive sequences that can not be detected based on unique tags alone, including transposable elements, segmental duplications and peri-centromeric regions. This mapping approach should prove to be useful for increasing biological knowledge on the too often neglected repetitive genomic regions.Availability: http://esbg.gatech.edu/jordan/software/mapContact: king.jordan@biology.gatech.eduSupplementary Information: Supplementary data are available at Bioinformatics online.

DICER- and AGO3-dependent generation of retinoic acid–induced DR2 Alu RNAs regulates human stem cell proliferation

“…Both uniquely aligned reads and those aligned to less than 10, 100, 1000 or 10000 locations were analyzed. For reads aligned to multiple regions, a single location was chosen based on a Gibbs sampling algorithm 52 , using 5 or 10 iterations. The sequencing reads count on DR2 Alu was computed using intersectBed tool in BedTools 53 , and were normalized to the sequencing depth.…”

Subtelomeric CTCF and cohesin binding site organization using improved subtelomere assemblies and a novel annotation pipeline

“…To deal with this challenge, we used a strategy of assigning a mapping likelihood (mL tag) to reads equal to the inverse of its genome-wide mapping positions; in effect, splitting up a read and mapping an equal portion of it to all of its possible sites of true mapping (Wang et al 2010;Chung et al 2011). By using this alternative mapping strategy, we then build fragment densities to display on enrichment tracks and to call peaks (see Methods).…”

Do human transposable element small RNAs serve primarily as genome defenders or genome regulators?

“…Sequence tags were mapped to the human genome reference sequence and co-located with genes and TEs as described in the Materials and Methods section. A recently developed algorithm for mapping ambiguous tags was used to ensure maximal coverage of repetitive TE sequences for the short sequence tags used 35 . This algorithm ensures that the best single genomic location for each multi-mapping tag is chosen, thus ensuring deeper coverage of TE sequences than would be achieved if multi-mapping tags were discarded.…”

“…RNA sequence reads were mapped to the human genome reference sequence (NCBI36/hg18) using the program Bowtie 47 with a threshold of ≤ 2 mismatches allowed. The most likely mapping locations for reads that mapped to more than one location were rescued using the Gibbs sampling strategy for multi-mapping tags 35 . mRNA and sRNA sequence tags mapped and processed in this way were co-located with human gene and TE loci annotated in the UCSC Genome Browser 48 .…”