Megadepth: efficient coverage quantification for BigWigs and BAMs.
Ontology highlight
ABSTRACT: Motivation
A common way to summarize sequencing datasets is to quantify data lying within genes or other genomic intervals. This can be slow and can require different tools for different input file types.Results
Megadepth is a fast tool for quantifying alignments and coverage for BigWig and BAM/CRAM input files, using substantially less memory than the next-fastest competitor. Megadepth can summarize coverage within all disjoint intervals of the Gencode V35 gene annotation for more than 19,000 GTExV8 BigWig files in approximately one hour using 32 threads. Megadepth is available both as a command-line tool and as an R/Bioconductor package providing much faster quantification compared to the rtracklayer package.Availability
https://github.com/ChristopherWilks/megadepth,https://bioconductor.org/packages/megadepth.
SUBMITTER: Wilks C
PROVIDER: S-EPMC8528031 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We present a tool for processing NGS DNA-seq BAM files to quantify peaks and create coverage tracks from ChIP-seq, NS-seq, ATAC-seq, OK-seq, END-seq and replication timing data
2020-05-17 | GSE131417 | GEO
Project description:Because of incomplete reference transcriptomes, incomplete sequencing bias models, or other modeling defects, algorithms to infer isoform expression from RNA sequencing (RNA-seq) sometimes do not accurately model expression. We present a computational method to detect instances where a quantification algorithm could not completely explain the input reads. Our approach identifies regions where the read coverage significantly deviates from expectation. We call these regions "expression anomalies." We further present a method to attribute their cause to either the incompleteness of the reference transcriptome or algorithmic mistakes. We detect anomalies for 30 GEUVADIS and 16 Human Body Map samples. By correcting anomalies when possible, we reduce the number of falsely predicted instances of differential expression. Anomalies that cannot be corrected are suspected to indicate the existence of isoforms unannotated by the reference. We detected 88 common anomalies of this type and find that they tend to have a lower-than-expected coverage toward their 3' ends.
| S-EPMC6938679 | biostudies-literature
Project description:BackgroundNext-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches.ResultsWe have developed BAMscale, a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale, we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~?100 Gb size) in minutes, outperforming currently available tools.ConclusionsBAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github (https://github.com/ncbi/BAMscale).
| S-EPMC7175505 | biostudies-literature
Project description:Gene targeting (GT) can introduce subtle alterations into a particular locus and represents a powerful tool for genome editing. Engineered zinc finger nucleases (ZFNs) are effective for generating minor allelic alterations. Efficient detection of such minor alterations remains one of the challenges in ZFN-mediated GT experiments. Here, we report the establishment of procedures allowing for efficient detection, quantification and enrichment of such subtle alterations. In a biallelic model, polyacrylamide gel electrophoresis (PAGE) is capable of detecting rare allelic variations in the form of DNA heteroduplexes at a high efficiency of ~0.4% compared with ~6.3% by the traditional T7 endonuclease I-digestion and agarose gel electrophoresis. In a multiple allelic model, PAGE could discriminate different alleles bearing addition or deletion of 1-18 bp as distinct bands that were easily quantifiable by densitometry. Furthermore, PAGE enables enrichment for rare alleles. We show for the first time that direct endogenous GT is possible in medaka by ZFN RNA injection, whereas PAGE allows for detection and cloning of ZFN-targeted alleles in adults arising from ZFN-injected medaka embryos. Therefore, PAGE is effective for detection, quantification and enrichment of multiple fine allelic differences and thus offers a versatile tool for screening targeted subtle gene alterations.
| S-EPMC3473374 | biostudies-literature
Project description:Accurate quantification of myocardial perfusion is dependent on reliable electrocardiogram (ECG) triggering. Measuring myocardial blood flow (MBF) in patients with arrhythmias or poor ECGs is currently infeasible with MR. The purpose of this study was to demonstrate the feasibility of a non-ECG-triggered method with clinically useful three-slice ventricular coverage for measurement of MBF in healthy volunteers.A saturation recovery magnetization-prepared gradient recalled echo acquisition was continuously repeated during first-pass imaging. A slice-interleaved radial trajectory was employed to enable image-based retrospective triggering. The arterial input function was generated using a beat-by-beat T1 estimation method. The proposed technique was validated against a conventional ECG-triggered dual-bolus technique in 10 healthy volunteers. The technique was further demonstrated under adenosine stress in 12 healthy volunteers.The proposed method produced MBF with no significant difference compared with the ECG-triggered technique. The proposed method yielded mean myocardial perfusion reserve comparable to published literature.We have developed a non-ECG-triggered quantitative perfusion imaging method. In this preliminary study, our results demonstrate that our method yields comparable MBF compared with the conventional ECG-triggered method and that it is feasible for stress imaging.
| S-EPMC4766039 | biostudies-literature
Project description:Research on the genetics of natural populations was revolutionized in the 1990s by methods for genotyping noninvasively collected samples. However, these methods have remained largely unchanged for the past 20 years and lag far behind the genomics era. To close this gap, here we report an optimized laboratory protocol for genome-wide capture of endogenous DNA from noninvasively collected samples, coupled with a novel computational approach to reconstruct pedigree links from the resulting low-coverage data. We validated both methods using fecal samples from 62 wild baboons, including 48 from an independently constructed extended pedigree. We enriched fecal-derived DNA samples up to 40-fold for endogenous baboon DNA and reconstructed near-perfect pedigree relationships even with extremely low-coverage sequencing. We anticipate that these methods will be broadly applicable to the many research systems for which only noninvasive samples are available. The lab protocol and software ("WHODAD") are freely available at www.tung-lab.org/protocols-and-software.html and www.xzlab.org/software.html, respectively.
| S-EPMC4896188 | biostudies-literature
Project description:Magnetic resonance fingerprinting (MRF) is a promising method for fast simultaneous quantification of multiple tissue parameters. The objective of this study is to improve the coverage of MRF based on echo-planar imaging (MRF-EPI) by using a slice-interleaved acquisition scheme. For this, the MRF-EPI is modified to acquire several slices in a randomized interleaved manner, increasing the effective repetition time of the spoiled gradient echo readout acquisition in each slice. Per-slice matching of the signal-trace to a precomputed dictionary allows the generation of T1 and T2* maps with integrated B1+ correction. Subsequent compensation for the coil sensitivity profile and normalization to the cerebrospinal fluid additionally allows for quantitative proton density (PD) mapping. Numerical simulations are performed to optimize the number of interleaved slices. Quantification accuracy is validated in phantom scans and feasibility is demonstrated in-vivo. Numerical simulations suggest the acquisition of four slices as a trade-off between quantification precision and scan-time. Phantom results indicate good agreement with reference measurements (Difference T1: -2.4?±?1.1%, T2*: -0.5?±?2.5%, PD: -0.5?±?7.2%). In-vivo whole-brain coverage of T1, T2* and PD with 32 slices was acquired within 3:36?minutes, resulting in parameter maps of high visual quality and comparable performance with single-slice MRF-EPI at 4-fold scan-time reduction.
| S-EPMC5923901 | biostudies-other
Project description:The generation and analysis of mutants in zebrafish has been instrumental in defining the genetic regulation of vertebrate development, physiology, and disease. However, identifying the genetic changes that underlie mutant phenotypes remains a significant bottleneck in the analysis of mutants. Whole-genome sequencing has recently emerged as a fast and efficient approach for identifying mutations in nonvertebrate model organisms. However, this approach has not been applied to zebrafish due to the complicating factors of having a large genome and lack of fully inbred lines. Here we provide a method for efficiently mapping and detecting mutations in zebrafish using these new parallel sequencing technologies. This method utilizes an extensive reference SNP database to define regions of homozygosity-by-descent by low coverage, whole-genome sequencing of pooled DNA from only a limited number of mutant F(2) fish. With this approach we mapped each of the five different zebrafish mutants we sequenced and identified likely causative nonsense mutations in two and candidate mutations in the remainder. Furthermore, we provide evidence that one of the identified mutations, a nonsense mutation in bmp1a, underlies the welded mutant phenotype.
| S-EPMC3296239 | biostudies-literature
Project description:The use of alternative poly(A) sites is common and affects the post-transcriptional fate of mRNA, including its stability, subcellular localization and translation. Here, we present a method to identify poly(A) sites in a genome-wide and strand-specific manner. This method, termed 3'T-fill, initially fills in the poly(A) stretch with unlabeled dTTPs, allowing sequencing to start directly after the poly(A) tail into the 3'-untranslated regions (UTR). Our comparative analysis demonstrates that it outperforms existing protocols in quality and throughput and accurately quantifies RNA levels as only one read is produced from each transcript. We use this method to characterize the diversity of polyadenylation in Saccharomyces cerevisiae, showing that alternative RNA molecules are present even in a genetically identical cell population. Finally, we observe that overlap of convergent 3'-UTRs is frequent but sharply limited by coding regions, suggesting factors that restrict compression of the yeast genome.
| S-EPMC3597643 | biostudies-literature
Project description:Protein phosphorylation is critically important for many cellular processes, including progression through the cell cycle, cellular metabolism, and differentiation. Isobaric labeling, for example, tandem mass tags (TMT), in phosphoproteomics workflows enables both relative and absolute quantitation of these phosphorylation events. Traditional TMT workflows identify peptides using fragment ions at the MS2 level and quantify reporter ions at the MS3 level. However, in addition to the TMT reporter ions, MS3 spectra also include fragment ions that can be used to identify peptides. Here we describe using MS3 spectra for both phosphopeptide identification and quantification, a process that we term MS3-IDQ. To maximize quantified phosphopeptides, we optimize several instrument parameters, including the modality of mass analyzer (i.e., ion trap or Orbitrap), MS2 automatic gain control (AGC), and MS3 normalized collision energy (NCE), to achieve the best balance of identified and quantified peptides. Our optimized MS3-IDQ method included the following parameters for the MS3 scan: NCE = 37.5 and AGC target = 1.5 × 105, and scan range = 100-2000. Data from the MS3 scan were complementary to those of the MS2 scan, and the combination of these scans can increase phosphoproteome coverage by >50%, thereby yielding a greater number of quantified and accurately localized phosphopeptides.
| S-EPMC5984591 | biostudies-literature