Dataset Information

MMR: a tool for read multi-mapper resolution.

ABSTRACT:

Motivation

Mapping high-throughput sequencing data to a reference genome is an essential step for most analysis pipelines aiming at the computational analysis of genome and transcriptome sequencing data. Breaking ties between equally well mapping locations poses a severe problem not only during the alignment phase but also has significant impact on the results of downstream analyses. We present the multi-mapper resolution (MMR) tool that infers optimal mapping locations from the coverage density of other mapped reads.

Results

Filtering alignments with MMR can significantly improve the performance of downstream analyses like transcript quantitation and differential testing. We illustrate that the accuracy (Spearman correlation) of transcript quantification increases by 15% when using reads of length 51. In addition, MMR decreases the alignment file sizes by more than 50%, and this leads to a reduced running time of the quantification tool. Our efficient implementation of the MMR algorithm is easily applicable as a post-processing step to existing alignment files in BAM format. Its complexity scales linearly with the number of alignments and requires no further inputs.

Availability and implementation

Open source code and documentation are available for download at http://github.com/ratschlab/mmr Comprehensive testing results and further information can be found at http://bioweb.me/mmr.

Contact

andre.kahles@ratschlab.org or gunnar.ratsch@ratschlab.org

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Kahles A

PROVIDER: S-EPMC4795617 | biostudies-literature | 2016 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

MMR: a tool for read multi-mapper resolution.

Kahles André A Behr Jonas J Rätsch Gunnar G

Bioinformatics (Oxford, England) 20151030 5

<h4>Motivation</h4>Mapping high-throughput sequencing data to a reference genome is an essential step for most analysis pipelines aiming at the computational analysis of genome and transcriptome sequencing data. Breaking ties between equally well mapping locations poses a severe problem not only during the alignment phase but also has significant impact on the results of downstream analyses. We present the multi-mapper resolution (MMR) tool that infers optimal mapping locations from the coverage ...[more]

PMID: 26519503

Similar Datasets

Project description:BackgroundIn the current era of high throughput genomics a major challenge is the genome-wide identification of target genes for specific transcription factors. Chromatin immunoprecipitation (ChIP) allows the isolation of in vivo binding sites of transcription factors and provides a powerful tool for examining gene regulation. Crosslinked chromatin is immunoprecipitated with antibodies against specific transcription factors, thus enriching for sequences bound in vivo by these factors in the immunoprecipitated DNA. Cloning and sequencing the immunoprecipitated sequences allows identification of transcription factor target genes. Routinely, thousands of such sequenced clones are used in BLAST searches to map their exact location in the genome and the genes located in the vicinity. These genes represent potential targets of the transcription factor of interest. Such bioinformatics analysis is very laborious if performed manually and for this reason there is a need for developing bioinformatic tools to automate and facilitate it.ResultsIn order to facilitate this analysis we generated TF Target Mapper (Transcription Factor Target Mapper). TF Target Mapper is a BLAST search tool allowing rapid extraction of annotated information on genes around each hit. It combines sequence cleaning/filtering, pattern searching and BLAST searches with extraction of information on genes located around each BLAST hit and comparisons of the output list of genes or gene ontology IDs with user-implemented lists. We successfully applied and tested TF Target Mapper to analyse sequences bound in vivo by the transcription factor GATA-1. We show that TF Target Mapper efficiently extracted information on genes around ChIPed sequences, thus identifying known (e.g. alpha-globin and zeta-globin) and potentially novel GATA-1 gene targets.ConclusionTF Target Mapper is a very efficient BLAST search tool that allows the rapid extraction of annotated information on the genes around each hit. It can contribute to the comprehensive bioinformatic transcriptome/regulome analysis, by providing insight into the mechanisms of action of specific transcription factors, thus helping to elucidate the pathways these factors regulate.

Project description:In Hungary, between February 2017 and July 2019, 70 confirmed measles cases were reported, raising questions about the adequacy of population-level immunity. Although the assumed vaccination coverage is ≥99%, in a recent study, we detected potential gaps in the anti-measles humoral immunity. In Hungary, according to a decree by the Ministry of Public Welfare, beginning from 2021, the healthcare provider should conduct a serosurvey of anti-measles protection levels of healthcare professionals. To facilitate the compliance with this requirement, we developed a quick 'three-in-one' or 'triple' MMR (measles, mumps and rubella) indirect ELISA (IgG); an assay format that is currently not available commercially. High throughput applicability of the 'three-in-one' ELISA was verified using 1736 sera from routine laboratory residual samples, using an automated platform (Siemens BEP 2000 Advance). Assay verification was performed by comparing the full antigen repertoire-based 'target' assay with in-house 'control' assays using recombinant viral antigen coatings, and by validated commercially available kits. Indirect immunofluorescence was used as an independent reference method. Data were analysed using OriginLab, IBM SPSS, RStudio and MedCalc. In case of measles, we combined our current results with previously published data (Ntotal measles = 3523). Evaluation of anti-mumps and anti-rubella humoral antibody levels was based on the measurement of 1736 samples. The lowest anti-measles seropositivity (79.3%) was detected in sera of individuals vaccinated between 1978 and 1987. Considering the antigen-specific seropositivity ratios of all samples measured, anti-measles, -mumps and -rubella IgG antibody titres were adequate in 89.84%, 91.82% and 92.28%, respectively. Based on the virus-specific herd immunity threshold (HIT) values (HITMeasles = 92-95%, HITMumps = 75-86%, HITRubella = 83-86), it can be stated that regarding anti-measles immunity, certain age clusters of the population may have inadequate levels of humoral immunity. Despite the potential gaps in herd immunity, the use of MMR vaccine remains an effective and low-cost approach for the prevention of measles, mumps and rubella infections.

Project description:BackgroundDe novo genome assembly is essential to modern genomics studies. As it is not biased by a reference, it is also a useful method for studying genomes with high variation, such as cancer genomes. De novo short-read assemblers commonly use de Bruijn graphs, where nodes are sequences of equal length k, also known as k-mers. Edges in this graph are established between nodes that overlap by [Formula: see text] bases, and nodes along unambiguous walks in the graph are subsequently merged. The selection of k is influenced by multiple factors, and optimizing this value results in a trade-off between graph connectivity and sequence contiguity. Ideally, multiple k sizes should be used, so lower values can provide good connectivity in lesser covered regions and higher values can increase contiguity in well-covered regions. However, current approaches that use multiple k values do not address the scalability issues inherent to the assembly of large genomes.ResultsHere we present RResolver, a scalable algorithm that takes a short-read de Bruijn graph assembly with a starting k as input and uses a k value closer to that of the read length to resolve repeats. RResolver builds a Bloom filter of sequencing reads which is used to evaluate the assembly graph path support at branching points and removes paths with insufficient support. RResolver runs efficiently, taking only 26 min on average for an ABySS human assembly with 48 threads and 60 GiB memory. Across all experiments, compared to a baseline assembly, RResolver improves scaffold contiguity (NGA50) by up to 15% and reduces misassemblies by up to 12%.ConclusionsRResolver adds a missing component to scalable de Bruijn graph genome assembly. By improving the initial and fundamental graph traversal outcome, all downstream ABySS algorithms greatly benefit by working with a more accurate and less complex representation of the genome. The RResolver code is integrated into ABySS and is available at https://github.com/bcgsc/abyss/tree/master/RResolver .

Dataset Information

MMR: a tool for read multi-mapper resolution.

Motivation

Results

Availability and implementation

Contact

Supplementary information

Publications

MMR: a tool for read multi-mapper resolution.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets