NgsRelate: a software tool for estimating pairwise relatedness from next-generation sequencing data.
Ontology highlight
ABSTRACT:
SUBMITTER: Korneliussen TS
PROVIDER: S-EPMC4673978 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Estimation of relatedness between pairs of individuals is important in many genetic research areas. When estimating relatedness, it is important to account for admixture if this is present. However, the methods that can account for admixture are all based on genotype data as input, which is a problem for low-depth next-generation sequencing (NGS) data from which genotypes are called with high uncertainty. Here, we present a software tool, NGSremix, for maximum likelihood estimation of relatedness between pairs of admixed individuals from low-depth NGS data, which takes the uncertainty of the genotypes into account via genotype likelihoods. Using both simulated and real NGS data for admixed individuals with an average depth of 4x or below we show that our method works well and clearly outperforms all the commonly used state-of-the-art relatedness estimation methods PLINK, KING, relateAdmix, and ngsRelate that all perform quite poorly. Hence, NGSremix is a useful new tool for estimating relatedness in admixed populations from low-depth NGS data. NGSremix is implemented in C/C++ in a multi-threaded software and is freely available on Github https://github.com/KHanghoj/NGSremix.
| S-EPMC8496226 | biostudies-literature
Project description:BACKGROUND:Multiple RNA samples are frequently processed together and often mixed before multiplex sequencing in the same sequencing run. While different samples can be separated post sequencing using sample barcodes, the possibility of cross contamination between biological samples from different species that have been processed or sequenced in parallel has the potential to be extremely deleterious for downstream analyses. RESULTS:We present CroCo, a software package for identifying and removing such cross contaminants from assembled transcriptomes. Using multiple, recently published sequence datasets, we show that cross contamination is consistently present at varying levels in real data. Using real and simulated data, we demonstrate that CroCo detects contaminants efficiently and correctly. Using a real example from a molecular phylogenetic dataset, we show that contaminants, if not eliminated, can have a decisive, deleterious impact on downstream comparative analyses. CONCLUSIONS:Cross contamination is pervasive in new and published datasets and, if undetected, can have serious deleterious effects on downstream analyses. CroCo is a database-independent, multi-platform tool, designed for ease of use, that efficiently and accurately detects and removes cross contamination in assembled transcriptomes to avoid these problems. We suggest that the use of CroCo should become a standard cleaning step when processing multiple samples for transcriptome sequencing.
| S-EPMC5838952 | biostudies-literature
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:Population structure can be described by genotypic-correlation coefficients between groups of individuals, the most basic of which are the pairwise relatedness coefficients between any two individuals. There are nine pairwise relatedness coefficients in the most general model, and we show that these can be reduced to seven coefficients for biallelic loci. Although all nine coefficients can be estimated from pedigrees, six coefficients have been beyond empirical reach. We provide a numerical optimization procedure that estimates all seven reduced coefficients from population-genomic data. Simulations show that the procedure is nearly unbiased, even at 3× coverage, and errors in five of the seven coefficients are statistically uncorrelated. The remaining two coefficients have a negative correlation of errors, but their sum provides an unbiased assessment of the overall correlation of heterozygosity between two individuals. Application of these new methods to four populations of the freshwater crustacean Daphnia pulex reveal the occurrence of half siblings in our samples, as well as a number of identical individuals that are likely obligately asexual clone mates. Statistically significant negative estimates of these pairwise relatedness coefficients, including inbreeding coefficients that were typically negative, underscore the difficulties that arise when interpreting genotypic correlations as estimations of the probability that alleles are identical by descent.
| S-EPMC5419463 | biostudies-literature
Project description:High-throughput next-generation sequencing technologies pose increasing demands on the efficiency, accuracy and usability of data analysis software. In this article, we present ZOOM Lite, a software for efficient reads mapping and result visualization. With a kernel capable of mapping tens of millions of Illumina or AB SOLiD sequencing reads efficiently and accurately, and an intuitive graphical user interface, ZOOM Lite integrates reads mapping and result visualization into a easy to use pipeline on desktop PC. The software handles both single-end and paired-end reads, and can output both the unique mapping result or the top N mapping results for each read. Additionally, the software takes a variety of input file formats and outputs to several commonly used result formats. The software is freely available at http://bioinfor.com/zoom/lite/.
| S-EPMC2896157 | biostudies-literature
Project description:PurposeGenotyping CYP2D6 is important for precision drug therapy because the enzyme it encodes metabolizes approximately 25% of drugs, and its activity varies considerably among individuals. Genotype analysis of CYP2D6 is challenging due to its highly polymorphic nature. Over 100 haplotypes (star alleles) have been defined for CYP2D6, some involving a gene conversion with its nearby nonfunctional but highly homologous paralog CYP2D7. We present Stargazer, a new bioinformatics tool that uses next-generation sequencing (NGS) data to call star alleles for CYP2D6 ( https://stargazer.gs.washington.edu/stargazerweb/ ). Stargazer is currently being extended for other pharmacogenes.MethodsStargazer identifies star alleles from NGS data by detecting single nucleotide variants, insertion-deletion variants, and structural variants. Stargazer detects structural variation, including gene deletions, duplications, and conversions, by calculating paralog-specific copy numbers from read depths.ResultsWe applied Stargazer to the NGS data of 32 ethnically diverse HapMap trios that were genotyped by TaqMan assays, long-range polymerase chain reaction, quantitative multiplex polymerase chain reaction, high-resolution melting analysis, and/or Sanger sequencing. CYP2D6 genotyping by Stargazer was 99.0% concordant with the data obtained by these methods, and showed that 28.1% of the samples had structural variation including CYP2D6/CYP2D7 hybrids.ConclusionAccurate genotyping of pharmacogenes with NGS and subsequent allele calling with Stargazer will aid the implementation of precision drug therapy.
| S-EPMC6281872 | biostudies-literature
Project description:SummaryWe 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.AvailabilityPurBayes is implemented as an R package, and source code is available for download through CRAN at http://cran.r-project.org/package=PurBayes.Contactlarson.nicholas@mayo.eduSupplementary informationSupplementary data are available online at Bioinformatics online.
| S-EPMC3712213 | biostudies-literature
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:Multi-locus sequence typing (MLST) is a widely used method for categorizing bacteria. Increasingly, MLST is being performed using next-generation sequencing (NGS) data by reference laboratories and for clinical diagnostics. Many software applications have been developed to calculate sequence types from NGS data; however, there has been no comprehensive review to date on these methods. We have compared eight of these applications against real and simulated data, and present results on: (1) the accuracy of each method against traditional typing methods, (2) the performance on real outbreak datasets, (3) the impact of contamination and varying depth of coverage, and (4) the computational resource requirements.
| S-EPMC5610716 | biostudies-literature