PurBayes: estimating tumor cellularity and subclonality in next-generation sequencing data.
Ontology highlight
ABSTRACT: Summary
We have developed a novel Bayesian method, PurBayes, to estimate tumor purity and detect intratumor heterogeneity based on next-generation sequencing data of paired tumor-normal tissue samples, which uses finite mixture modeling methods. We demonstrate our approach using simulated data and discuss its performance under varying conditions.Availability
PurBayes is implemented as an R package, and source code is available for download through CRAN at http://cran.r-project.org/package=PurBayes.Contact
larson.nicholas@mayo.eduSupplementary information
Supplementary data are available online at Bioinformatics online.
SUBMITTER: Larson NB
PROVIDER: S-EPMC3712213 | biostudies-literature | 2013 Aug
REPOSITORIES: biostudies-literature
Publications
PurBayes: estimating tumor cellularity and subclonality in next-generation sequencing data.
Larson Nicholas B NB Fridley Brooke L BL
Bioinformatics (Oxford, England) 20130606 15
<h4>Summary</h4>We have developed a novel Bayesian method, PurBayes, to estimate tumor purity and detect intratumor heterogeneity based on next-generation sequencing data of paired tumor-normal tissue samples, which uses finite mixture modeling methods. We demonstrate our approach using simulated data and discuss its performance under varying conditions.<h4>Availability</h4>PurBayes is implemented as an R package, and source code is available for download through CRAN at http://cran.r-project.or ...[more]
Similar Datasets
Project description:Inference of population structure and individual ancestry is important both for population genetics and for association studies. With next generation sequencing technologies it is possible to obtain genetic data for all accessible genetic variations in the genome. Existing methods for admixture analysis rely on known genotypes. However, individual genotypes cannot be inferred from low-depth sequencing data without introducing errors. This article presents a new method for inferring an individual's ancestry that takes the uncertainty introduced in next generation sequencing data into account. This is achieved by working directly with genotype likelihoods that contain all relevant information of the unobserved genotypes. Using simulations as well as publicly available sequencing data, we demonstrate that the presented method has great accuracy even for very low-depth data. At the same time, we demonstrate that applying existing methods to genotypes called from the same data can introduce severe biases. The presented method is implemented in the NGSadmix software available at http://www.popgen.dk/software.
| S-EPMC3813857 | biostudies-other
Project description:BACKGROUND:DNA methylation (DNAm) age has been widely accepted as an epigenetic biomarker for biological aging. Emerging evidence suggests that DNAm age can be tissue-specific and female breast tissue ages faster than other parts of the body. The Horvath clock, which estimates DNAm age across multiple tissues, has been shown to be poorly calibrated in breast issue. We aim to develop a model to estimate breast tissue-specific DNAm age. METHODS:Genome-wide DNA methylation sequencing data were generated for 459 normal, 107 tumor, and 45 paired adjacent-normal breast tissue samples. We determined a novel set of 286 breast tissue-specific clock CpGs using penalized linear regression and developed a model to estimate breast tissue-specific DNAm age. The model was applied to estimate breast tissue-specific DNAm age in different breast tissue types and in tumors with distinct clinical characteristics to investigate cancer-related aging effects. RESULTS:Our estimated breast tissue-specific DNAm age was highly correlated with chronological age (r = 0.88; p = 2.9 × 10-31) in normal breast tissue. Breast tumor tissue samples exhibited a positive epigenetic age acceleration, where DNAm age was on average 7?years older than respective chronological age (p = 1.8 × 10-8). In age-matched analyses, tumor breast tissue appeared 12 and 13?years older in DNAm age than adjacent-normal and normal breast tissue (p = 4.0 × 10-6 and 1.0 × 10-6, respectively). Both HER2+ and hormone-receptor positive subtypes demonstrated significant acceleration in DNAm ages (p = 0.04 and 3.8 × 10-6, respectively), while no apparent DNAm age acceleration was observed for triple-negative breast tumors. We observed a non-linear pattern of epigenetic age acceleration with breast tumor grade. In addition, early-staged tumors showed a positive epigenetic age acceleration (p = 0.003) while late-staged tumors exhibited a non-significant negative epigenetic age acceleration (p = 0.10). CONCLUSIONS:The intended applications for this model are wide-spread and have been shown to provide biologically meaningful results for cancer-related aging effects in breast tumor tissue. Future studies are warranted to explore whether breast tissue-specific epigenetic age acceleration is predictive of breast cancer development, treatment response, and survival as well as the clinical utility of whether this model can be extended to blood samples.
| S-EPMC7282053 | biostudies-literature
Project description:SummaryHere, we present ContEst, a tool for estimating the level of cross-individual contamination in next-generation sequencing data. We demonstrate the accuracy of ContEst across a range of contamination levels, sources and read depths using sequencing data mixed in silico at known concentrations. We applied our tool to published cancer sequencing datasets and report their estimated contamination levels.Availability and implementationContEst is a GATK module, and distributed under a BSD style license at http://www.broadinstitute.org/cancer/cga/contestContactkcibul@broadinstitute.org; gadgetz@broadinstitute.orgSupplementary informationSupplementary data is available at Bioinformatics online.
| S-EPMC3167057 | biostudies-literature
Project description:BackgroundCancer neoantigens are expressed only in cancer cells and presented on the tumor cell surface in complex with major histocompatibility complex (MHC) class I proteins for recognition by cytotoxic T cells. Accurate and rapid identification of neoantigens play a pivotal role in cancer immunotherapy. Although several in silico tools for neoantigen prediction have been presented, limitations of these tools exist.ResultsWe developed pTuneos, a computational pipeline for prioritizing tumor neoantigens from next-generation sequencing data. We tested the performance of pTuneos on the melanoma cancer vaccine cohort data and tumor-infiltrating lymphocyte (TIL)-recognized neopeptide data. pTuneos is able to predict the MHC presentation and T cell recognition ability of the candidate neoantigens, and the actual immunogenicity of single-nucleotide variant (SNV)-based neopeptides considering their natural processing and presentation, surpassing the existing tools with a comprehensive and quantitative benchmark of their neoantigen prioritization performance and running time. pTuneos was further tested on The Cancer Genome Atlas (TCGA) cohort data as well as the melanoma and non-small cell lung cancer (NSCLC) cohort data undergoing checkpoint blockade immunotherapy. The overall neoantigen immunogenicity score proposed by pTuneos is demonstrated to be a powerful and pan-cancer marker for survival prediction compared to traditional well-established biomarkers.ConclusionsIn summary, pTuneos provides the state-of-the-art one-stop and user-friendly solution for prioritizing SNV-based candidate neoepitopes, which could help to advance research on next-generation cancer immunotherapies and personalized cancer vaccines. pTuneos is available at https://github.com/bm2-lab/pTuneos , with a Docker version for quick deployment at https://cloud.docker.com/u/bm2lab/repository/docker/bm2lab/ptuneos .
| S-EPMC6822339 | biostudies-literature
Project description:Pairwise relatedness estimation is important in many contexts such as disease mapping and population genetics. However, all existing estimation methods are based on called genotypes, which is not ideal for next-generation sequencing (NGS) data of low depth from which genotypes cannot be called with high certainty.We present a software tool, NgsRelate, for estimating pairwise relatedness from NGS data. It provides maximum likelihood estimates that are based on genotype likelihoods instead of genotypes and thereby takes the inherent uncertainty of the genotypes into account. Using both simulated and real data, we show that NgsRelate provides markedly better estimates for low-depth NGS data than two state-of-the-art genotype-based methods.NgsRelate is implemented in C++ and is available under the GNU license at www.popgen.dk/software.
| S-EPMC4673978 | biostudies-literature
Project description:Neoplastic cellularity contributes to the analytic sensitivity of most present technologies for mutation detection, such that they underperform when stroma and inflammatory cells dilute a cancer specimen's variant fraction. Thus, tumor purity assessment by light microscopy is used to determine sample adequacy before sequencing and to interpret the significance of negative results and mutant allele fraction afterwards. However, pathologist estimates of tumor purity are imprecise and have limited reproducibility. With the advent of massively parallel sequencing, large amounts of molecular data can be analyzed by computational purity algorithms. We retrospectively compared tumor purity of 3 computational algorithms with neoplastic cellularity using hematoxylin and eosin light microscopy to determine which was best for clinical evaluation of molecular profiling. Data were analyzed from 881 cancer patients from a clinical trial cohort, LCCC1108 (UNCseq), whose tumors had targeted massively parallel sequencing. Concordance among algorithms was poor, and the specimens analyzed had high rates of algorithm failure partially due to variable tumor purity. Computational tumor purity estimates did not add value beyond the pathologist's estimate of neoplastic cellularity microscopy. To improve present methods, we propose a semiquantitative, clinically applicable strategy based on mutant allele fraction and copy number changes present within a given specimen, which when combined with the morphologic tumor purity estimate, guide the interpretation of next-generation sequencing results in cancer patients.
| S-EPMC6887630 | biostudies-literature
Project description:SummaryWith the rapid expansion of the capabilities of the DNA sequencers throughout the different sequencing generations, the quantity of generated data has likewise increased. This evolution has also led to new bioinformatical methods, for which in silico data have become crucial when verifying the accuracy of a model or the robustness of a genomic analysis pipeline. Here, we present a multithreaded next-generation simulator for next-generation sequencing data (NGSNGS), which simulates reads faster than currently available methods and programs. NGSNGS can simulate reads with platform-specific characteristics based on nucleotide quality score profiles as well as including a post-mortem damage model which is relevant for simulating ancient DNA. The simulated sequences are sampled (with replacement) from a reference DNA genome, which can represent a haploid genome, polyploid assemblies or even population haplotypes and allows the user to simulate known variable sites directly. The program is implemented in a multithreading framework and is factors faster than currently available tools while extending their feature set and possible output formats.Availability and implementationThe method and associated programs are released as open-source software, code and user manual are available at https://github.com/RAHenriksen/NGSNGS.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9891242 | biostudies-literature
Project description:UnlabelledNew applications of next-generation sequencing technologies use pools of DNA from multiple individuals to estimate population genetic parameters. However, no publicly available tools exist to analyse single-nucleotide polymorphism (SNP) calling results directly for evolutionary parameters important in detecting natural selection, including nucleotide diversity and gene diversity. We have developed SNPGenie to fill this gap. The user submits a FASTA reference sequence(s), a Gene Transfer Format (.GTF) file with CDS information and a SNP report(s) in an increasing selection of formats. The program estimates nucleotide diversity, distance from the reference and gene diversity. Sites are flagged for multiple overlapping reading frames, and are categorized by polymorphism type: nonsynonymous, synonymous, or ambiguous. The results allow single nucleotide, single codon, sliding window, whole gene and whole genome/population analyses that aid in the detection of positive and purifying natural selection in the source population.Availability and implementationSNPGenie version 1.2 is a Perl program with no additional dependencies. It is free, open-source, and available for download at https://github.com/hugheslab/snpgenie.Contactnelsoncw@email.sc.edu or austin@biol.sc.eduSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC4757956 | biostudies-literature
Project description:BackgroundNext-generation sequencing is revolutionising diagnosis and treatment of rare diseases, however its application to understanding common disease aetiology is limited. Rare disease applications binarily attribute genetic change(s) at a single locus to a specific phenotype. In common diseases, where multiple genetic variants within and across genes contribute to disease, binary modelling cannot capture the burden of pathogenicity harboured by an individual across a given gene/pathway. We present GenePy, a novel gene-level scoring system for integration and analysis of next-generation sequencing data on a per-individual basis that transforms NGS data interpretation from variant-level to gene-level. This simple and flexible scoring system is intuitive and amenable to integration for machine learning, network and topological approaches, facilitating the investigation of complex phenotypes.ResultsWhole-exome sequencing data from 508 individuals were used to generate GenePy scores. For each variant a score is calculated incorporating: i) population allele frequency estimates; ii) individual zygosity, determined through standard variant calling pipelines and; iii) any user defined deleteriousness metric to inform on functional impact. GenePy then combines scores generated for all variants observed into a single gene score for each individual. We generated a matrix of ~ 14,000 GenePy scores for all individuals for each of sixteen popular deleteriousness metrics. All per-gene scores are corrected for gene length. The majority of genes generate GenePy scores < 0.01 although individuals harbouring multiple rare highly deleterious mutations can accumulate extremely high GenePy scores. In the absence of a comparator metric, we examine GenePy performance in discriminating genes known to be associated with three common, complex diseases. A Mann-Whitney U test conducted on GenePy scores for this positive control gene in cases versus controls demonstrates markedly more significant results (p = 1.37 × 10- 4) compared to the most commonly applied association tool that combines common and rare variation (p = 0.003).ConclusionsPer-gene per-individual GenePy scores are intuitive when assessing genetic variation in individual patients or comparing scores between groups. GenePy outperforms the currently accepted best practice tools for combining common and rare variation. GenePy scores are suitable for downstream data integration with transcriptomic and proteomic data that also report at the gene level.
| S-EPMC6524327 | biostudies-literature
Project description:NGPS is a method for de-novo, full-length protein sequencing in high throughput. The method is based on cleavage of the protein at semi-random sites by microwave-assisted acid hydrolysis (MAAH), enrichment of LC-MS/MS amenable peptides from the hydrolysate by solid-phase-extraction, LC-MS/MS analysis, de-novo long peptide tag sequencing of resulting peptides and assembly of peptide tags into consensus contigs.
2017-04-03 | PXD003804 | Pride