Correcting for sequence biases in present/absent calls.
Ontology highlight
ABSTRACT:
SUBMITTER: Schuster EF
PROVIDER: S-EPMC2394774 | biostudies-other | 2007
REPOSITORIES: biostudies-other
Publications
Correcting for sequence biases in present/absent calls.
Schuster Eugene F EF Blanc Eric E Partridge Linda L Thornton Janet M JM
Genome biology 20070101 6
The probe sequence of short oligonucleotides in Affymetrix microarray experiments can have a significant influence on present/absent calls of probesets with absent target transcripts. Probesets enriched for central Ts and depleted of central As in the perfect-match probes tend to be falsely classified as having present transcripts. Correction of non-specific binding for both perfect-match and mismatch probes using probe-sequence models can partially remove the probe-sequence bias and result in b ...[more]
Similar Datasets
Project description:High-throughput sequence (HTS) data exhibit position-specific nucleotide biases that obscure the intended signal and reduce the effectiveness of these data for downstream analyses. These biases are particularly evident in HTS assays for identifying regulatory regions in DNA (DNase-seq, ChIP-seq, FAIRE-seq, ATAC-seq). Biases may result from many experiment-specific factors, including selectivity of DNA restriction enzymes and fragmentation method, as well as sequencing technology-specific factors, such as choice of adapters/primers and sample amplification methods.We present a novel method to detect and correct position-specific nucleotide biases in HTS short read data. Our method calculates read-specific weights based on aligned reads to correct the over- or underrepresentation of position-specific nucleotide subsequences, both within and adjacent to the aligned read, relative to a baseline calculated in assay-specific enriched regions. Using HTS data from a variety of ChIP-seq, DNase-seq, FAIRE-seq, and ATAC-seq experiments, we show that our weight-adjusted reads reduce the position-specific nucleotide imbalance across reads and improve the utility of these data for downstream analyses, including identification and characterization of open chromatin peaks and transcription-factor binding sites.A general-purpose method to characterize and correct position-specific nucleotide sequence biases fills the need to recognize and deal with, in a systematic manner, binding-site preference for the growing number of HTS-based epigenetic assays. As the breadth and impact of these biases are better understood, the availability of a standard toolkit to correct them will be important.
| S-EPMC5540620 | biostudies-literature
Project description:Visuo-haptic biases are observed when bringing your unseen hand to a visual target. The biases are different between, but consistent within participants. We investigated the usefulness of adjusting haptic guidance to these user-specific biases in aligning haptic and visual perception. By adjusting haptic guidance according to the biases, we aimed to reduce the conflict between the modalities. We first measured the biases using an adaptive procedure. Next, we measured performance in a pointing task using three conditions: 1) visual images that were adjusted to user-specific biases, without haptic guidance, 2) veridical visual images combined with haptic guidance, and 3) shifted visual images combined with haptic guidance. Adding haptic guidance increased precision. Combining haptic guidance with user-specific visual information yielded the highest accuracy and the lowest level of conflict with the guidance at the end point. These results show the potential of correcting for user-specific perceptual biases when designing haptic guidance.
| S-EPMC4954687 | biostudies-literature
Project description:Sensing environmental temperatures is essential for the survival of ectothermic organisms. In Drosophila, two of the most used methodologies to study temperature preferences (TP) and the genes involved in thermosensation are two-choice assays and temperature gradients. Whereas two-choice assays reveal a relative TP, temperature gradients can identify the absolute Tp. One drawback of gradients is that small ectothermic animals are susceptible to cold-trapping: a physiological inability to move at the cold area of the gradient. Often cold-trapping cannot be avoided, biasing the resulting TP to lower temperatures. Two mathematical models were previously developed to correct for cold-trapping. These models, however, focus on group behaviour which can lead to overestimation of cold-trapping due to group aggregation. Here we present a mathematical model that simulates the behaviour of individual Drosophila in temperature gradients. The model takes the spatial dimension and temperature difference of the gradient into account, as well as the rearing temperature of the flies. Furthermore, it allows the quantification of cold-trapping and reveals unbiased TP. Additionally, our model reveals that flies have a range of tolerable temperatures, and this measure is more informative about the behaviour than commonly used TP. Online simulation is hosted at http://igloo.uni-goettingen.de . The code can be accessed at https://github.com/zerotonin/igloo .
| S-EPMC6408449 | biostudies-literature
Project description:High-throughput reporter assays such as self-transcribing active regulatory region sequencing (STARR-seq) have made it possible to measure regulatory element activity across the entire human genome at once. The resulting data, however, present substantial analytical challenges. Here, we identify technical biases that explain most of the variance in STARR-seq data. We then develop a statistical model to correct those biases and to improve detection of regulatory elements. This approach substantially improves precision and recall over current methods, improves detection of both activating and repressive regulatory elements, and controls for false discoveries despite strong local correlations in signal.
| S-EPMC8092017 | biostudies-literature
Project description:Homology between mitochondrial DNA (mtDNA) and nuclear DNA of mitochondrial origin (nuMTs) causes confounding when aligning short sequence reads to the reference human genome, as the true sequence origin cannot be determined. Using a systematic in silico approach, we here report the impact of all potential mitochondrial variants on alignment accuracy and variant calling. A total of 49,707 possible mutations were introduced across the 16,569 bp reference mitochondrial genome (16,569 × 3 alternative alleles), one variant at-at-time. The resulting in silico fragmentation and alignment to the entire reference genome (GRCh38) revealed preferential mapping of mutated mitochondrial fragments to nuclear loci, as variants increased loci similarity to nuMTs, for a total of 807, 362, and 41 variants at 333, 144, and 27 positions when using 100, 150, and 300 bp single-end fragments. We subsequently modeled these affected variants at 50% heteroplasmy and carried out variant calling, observing bias in the reported allele frequencies in favor of the reference allele. Four variants (chrM:6023A, chrM:4456T, chrM:5147A, and chrM:7521A) including a possible hypertension factor, chrM:4456T, caused 100% loss of coverage at the mutated position (with all 100 bp single-end fragments aligning to homologous, nuclear positions instead of chrM), rendering these variants undetectable when aligning to the entire reference genome. Furthermore, four mitochondrial variants reported to be pathogenic were found to cause significant loss of coverage and select haplogroup-defining SNPs were shown to exacerbate the loss of coverage caused by surrounding variants. Increased fragment length and use of paired-end reads both improved alignment accuracy.
| S-EPMC6773831 | biostudies-literature
Project description:BackgroundA pangenome is the collection of all genes found in a set of related genomes. For microbes, these genomes are often different strains of the same species, and the pangenome offers a means to compare gene content variation with differences in phenotypes, ecology, and phylogenetic relatedness. Though most frequently applied to bacteria, there is growing interest in adapting pangenome analysis to bacteriophages. However, working with phage genomes presents new challenges. First, most phage families are under-sampled, and homologous genes in related viruses can be difficult to identify. Second, homing endonucleases and intron-like sequences may be present, resulting in fragmented gene calls. Each of these issues can reduce the accuracy of standard pangenome analysis tools.MethodsWe developed an R pipeline called Rephine.r that takes as input the gene clusters produced by an initial pangenomics workflow. Rephine.r then proceeds in two primary steps. First, it identifies three common causes of fragmented gene calls: (1) indels creating early stop codons and new start codons; (2) interruption by a selfish genetic element; and (3) splitting at the ends of the reported genome. Fragmented genes are then fused to create new sequence alignments. In tandem, Rephine.r searches for distant homologs separated into different gene families using Hidden Markov Models. Significant hits are used to merge families into larger clusters. A final round of fragment identification is then run, and results may be used to infer single-copy core genomes and phylogenetic trees.ResultsWe applied Rephine.r to three well-studied phage groups: the Tevenvirinae (e.g., T4), the Studiervirinae (e.g., T7), and the Pbunaviruses (e.g., PB1). In each case, Rephine.r recovered additional members of the single-copy core genome and increased the overall bootstrap support of the phylogeny. The Rephine.r pipeline is provided through GitHub (https://www.github.com/coevoeco/Rephine.r) as a single script for automated analysis and with utility functions to assist in building single-copy core genomes and predicting the sources of fragmented genes.
| S-EPMC8351571 | biostudies-literature
Project description:The Cancer Genome Atlas (TCGA) provides a genetic characterization of more than ten thousand tumors, enabling the discovery of novel driver mutations, molecular subtypes, and enticing drug targets across many histologies. Here we investigated why some mutations are common in particular cancer types but absent in others. As an example, we observed that the gene CCDC168 has no mutations in the stomach adenocarcinoma (STAD) cohort despite its common presence in other tumor types. Surprisingly, we found that the lack of called mutations was due to a systematic insufficiency in the number of sequencing reads in the STAD and other cohorts, as opposed to differential driver biology. Using strict filtering criteria, we found similar behavior in four other genes across TCGA cohorts, with each gene exhibiting systematic sequencing depth issues affecting the ability to call mutations. We identified the culprit as the choice of exome capture kit, as kit choice was highly associated with the set of genes that have insufficient reads to call a mutation. Overall, we found that thousands of samples across all cohorts are subject to some capture kit problems. For example, for the 6353 samples using the Broad Institute's Custom capture kit there are undercalling biases for at least 4833 genes. False negative mutation calls at these genes may obscure biological similarities between tumor types and other important cancer driver effects in TCGA datasets.
| S-EPMC6169918 | biostudies-literature
Project description:The main objective of the study was to find combinations of genetic variants significantly associated with bipolar disorder. In a previous study of bipolar disorder, combinations of three single nucleotide polymorphism (SNP) genotypes taken from 803 SNPs were analyzed, and four clusters of combinations were found to be significantly associated with bipolar disorder. In the present study, combinations of four SNP genotypes taken from the same 803 SNPs were analyzed, and one cluster of combinations was found to be significantly associated with bipolar disorder. Combinations from the new cluster and from the four previous clusters were identified in the genomes of 209 of the 607 patients in the study whereas none of the 1355 control participants had any of these combinations in their genome.
| S-EPMC4654514 | biostudies-literature
Project description:Global coupled chemistry-climate models underestimate carbon monoxide (CO) in the Northern Hemisphere, exhibiting a pervasive negative bias against measurements peaking in late winter and early spring. While this bias has been commonly attributed to underestimation of direct anthropogenic and biomass burning emissions, chemical production and loss via OH reaction from emissions of anthropogenic and biogenic volatile organic compounds (VOCs) play an important role. Here we investigate the reasons for this underestimation using aircraft measurements taken in May and June 2016 from the Korea-United States Air Quality (KORUS-AQ) experiment in South Korea and the Air Chemistry Research in Asia (ARIAs) in the North China Plain (NCP). For reference, multispectral CO retrievals (V8J) from the Measurements of Pollution in the Troposphere (MOPITT) are jointly assimilated with meteorological observations using an ensemble adjustment Kalman filter (EAKF) within the global Community Atmosphere Model with Chemistry (CAM-Chem) and the Data Assimilation Research Testbed (DART). With regard to KORUS-AQ data, CO is underestimated by 42% in the control run and by 12% with the MOPITT assimilation run. The inversion suggests an underestimation of anthropogenic CO sources in many regions, by up to 80% for northern China, with large increments over the Liaoning Province and the North China Plain (NCP). Yet, an often-overlooked aspect of these inversions is that correcting the underestimation in anthropogenic CO emissions also improves the comparison with observational O3 datasets and observationally constrained box model simulations of OH and HO2. Running a CAM-Chem simulation with the updated emissions of anthropogenic CO reduces the bias by 29% for CO, 18% for ozone, 11% for HO2, and 27% for OH. Longer-lived anthropogenic VOCs whose model errors are correlated with CO are also improved, while short-lived VOCs, including formaldehyde, are difficult to constrain solely by assimilating satellite retrievals of CO. During an anticyclonic episode, better simulation of O3, with an average underestimation of 5.5 ppbv, and a reduction in the bias of surface formaldehyde and oxygenated VOCs can be achieved by separately increasing by a factor of 2 the modeled biogenic emissions for the plant functional types found in Korea. Results also suggest that controlling VOC and CO emissions, in addition to widespread NO x controls, can improve ozone pollution over East Asia.
| S-EPMC7786812 | biostudies-literature
Project description:A previously described DNA sequence generator algorithm (DNA-SGA) using error-correcting codes has been employed as a computational tool to address the evolutionary pathway of the genetic code. The code-generated sequence alignment demonstrated that a residue mutation revealed by the code can be found in the same position in sequences of distantly related taxa. Furthermore, the code-generated sequences do not promote amino acid changes in the deviant genomes through codon reassignment. A Bayesian evolutionary analysis of both code-generated and homologous sequences of the Arabidopsis thaliana malate dehydrogenase gene indicates an approximately 1 MYA divergence time from the MDH code-generated sequence node to its paralogous sequences. The DNA-SGA helps to determine the plesiomorphic state of DNA sequences because a single nucleotide alteration often occurs in distantly related taxa and can be found in the alternative codon patterns of noncanonical genetic codes. As a consequence, the algorithm may reveal an earlier stage of the evolution of the standard code.
| S-EPMC4498232 | biostudies-literature