Systematic bias in high-throughput sequencing data and its correction by BEADS
Ontology highlight
ABSTRACT:
ORGANISM(S): Caenorhabditis elegans
PROVIDER: GSE29427 | GEO | 2011/05/25
SECONDARY ACCESSION(S): PRJNA141495
REPOSITORIES: GEO
Similar Datasets
Project description:ChIP-seq and input sequence data used in the development and evaluation of the BEADS normalization method. Examination of ChIP and input sequence reads across the worm genome
2011-05-25 | E-GEOD-29427 | biostudies-arrayexpress
Project description:Genomic sequences obtained through high-throughput sequencing are not uniformly distributed across the genome. For example, sequencing data of total genomic DNA show significant, yet unexpected enrichments on promoters and exons. This systematic bias is a particular problem for techniques such as chromatin immunoprecipitation, where the signal for a target factor is plotted across genomic features. We have focused on data obtained from Illumina's Genome Analyser platform, where at least three factors contribute to sequence bias: GC content, mappability of sequencing reads, and regional biases that might be generated by local structure. We show that relying on input control as a normalizer is not generally appropriate due to sample to sample variation in bias. To correct sequence bias, we present BEADS (bias elimination algorithm for deep sequencing), a simple three-step normalization scheme that successfully unmasks real binding patterns in ChIP-seq data. We suggest that this procedure be done routinely prior to data interpretation and downstream analyses.
| S-EPMC3159482 | biostudies-literature
Project description:mRNA abundance microarrays are the most widely used high-throughput technology for functional genomics. But despite this extensive characterization, and despite their wide-spread adoption in discovery research, the uptake of microarrays as clinical diagnostic tools remains relatively limited due to sample degradation. Recently, the NuGEN WT-Ovation Pico kit was developed to resolve this problem. Here, we comprehensively evaluate its utility on two separate biological problems. First, we evaluate its effect on 12 malignant and 2 normal cervix samples. Second, we evaluate its effect on 2 normal kidney samples. In each case, we directly compare aliquots of the same RNA sample that are prepared by standard and low-input/low-RIN protocols to address three key questions. First, does the use of a low-input/degraded-sample protocol alter transcriptional profiles either globally or of specific genes? Second, do these changes alter biological conclusions drawn from the dataset? Third, are these changes biased towards specific classes of genes, either in terms of their sequence characteristics (e.g. length, GC-content) or their biological function?
2014-01-01 | GSE42764 | GEO
Project description:A feature common to all DNA sequencing technologies is the presence of base-call errors in the sequenced reads. The implications of such errors are application specific, ranging from minor informatics nuisances to major problems affecting biological inferences. Recently developed "next-gen" sequencing technologies have greatly reduced the cost of sequencing, but have been shown to be more error prone than previous technologies. Both position specific (depending on the location in the read) and sequence specific (depending on the sequence in the read) errors have been identified in Illumina and Life Technology sequencing platforms. We describe a new type of systematic error that manifests as statistically unlikely accumulations of errors at specific genome (or transcriptome) locations.We characterize and describe systematic errors using overlapping paired reads from high-coverage data. We show that such errors occur in approximately 1 in 1000 base pairs, and that they are highly replicable across experiments. We identify motifs that are frequent at systematic error sites, and describe a classifier that distinguishes heterozygous sites from systematic error. Our classifier is designed to accommodate data from experiments in which the allele frequencies at heterozygous sites are not necessarily 0.5 (such as in the case of RNA-Seq), and can be used with single-end datasets.Systematic errors can easily be mistaken for heterozygous sites in individuals, or for SNPs in population analyses. Systematic errors are particularly problematic in low coverage experiments, or in estimates of allele-specific expression from RNA-Seq data. Our characterization of systematic error has allowed us to develop a program, called SysCall, for identifying and correcting such errors. We conclude that correction of systematic errors is important to consider in the design and interpretation of high-throughput sequencing experiments.
| S-EPMC3295828 | biostudies-literature
Project description:Droplet-based single cell sequencing technologies, such as inDrop, Drop-seq, and 10X Genomics, are catalyzing a revolution in the understanding of biology. Barcoding beads are key components for these technologies. What is limiting today are barcoding beads that are easy to fabricate, can efficiently deliver primers into drops, and thus achieve high detection efficiency. Here, this work reports an approach to fabricate dissolvable polyacrylamide beads, by crosslinking acrylamide with disulfide bridges that can be cleaved with dithiothreitol. The beads can be rapidly dissolved in drops and release DNA barcode primers. The dissolvable beads are easy to synthesize, and the primer cost for the beads is significantly lower than that for the previous barcoding beads. Furthermore, the dissolvable beads can be loaded into drops with >95% loading efficiency of a single bead per drop and the dissolution of beads does not influence reverse transcription or the polymerase chain reaction (PCR) in drops. Based on this approach, the dissolvable beads are used for single cell RNA and protein analysis.
| S-EPMC7175265 | biostudies-literature
Project description:mRNA abundance microarrays are the most widely used high-throughput technology for functional genomics. But despite this extensive characterization, and despite their wide-spread adoption in discovery research, the uptake of microarrays as clinical diagnostic tools remains relatively limited due to sample degradation. Recently, the NuGEN WT-Ovation Pico kit was developed to resolve this problem. Here, we comprehensively evaluate its utility on two separate biological problems. First, we evaluate its effect on 12 malignant and 2 normal cervix samples. Second, we evaluate its effect on 2 normal kidney samples. In each case, we directly compare aliquots of the same RNA sample that are prepared by standard and low-input/low-RIN protocols to address three key questions. First, does the use of a low-input/degraded-sample protocol alter transcriptional profiles either globally or of specific genes? Second, do these changes alter biological conclusions drawn from the dataset? Third, are these changes biased towards specific classes of genes, either in terms of their sequence characteristics (e.g. length, GC-content) or their biological function? To achieve these goals, we selected 12 frozen cervix cancers and 2 unmatched normal samples. RNA was extracted and was split into two aliquots: one was prepared using the Nugen WT-Ovation Pico kit and the other using the standard Affymetrix Express array method. This approach allows us to directly assess the effect of sample preparation, in isolation of any other effects. These data were then pre-processed and the filtered genes were subject to linear-modeling to evaluate technical and biological variability. The resulting gene sets were then subject to pathway and sequence analysis, and LTR analysis. Lastly, to generalize these conclusions we repeated this experimental approach identically in two normal kidney samples and one RNA mixture control sample.
2014-01-01 | E-GEOD-42764 | biostudies-arrayexpress
Project description:To evaluate effect of lipopolysaccharide in nasal fibroblast, we have employed whole genome microarray expression profiling as a discovery platform. Nasal fibroblasts were exposed to LPS (10 μg/mL) for 12 h. Total RNA was isolated using Trizol reagent (Invitrogen, Carlsbad, CA). For control and test RNAs, synthesis of target cRNA probes and hybridization were performed using the Low RNA Input Linear Amplification kit (Agilent Technology, Santa Clara, CA). Hybridized images were scanned using a DNA microarray scanner and quantified using Feature Extraction Software (Agilent). All data normalization and selection of fold-change of the genes were performed using GeneSpringGX 7.3 (Agilent).
2015-12-31 | GSE52505 | GEO
Project description:Spatial bias continues to be a major challenge in high-throughput screening technologies. Its successful detection and elimination are critical for identifying the most promising drug candidates. Here, we examine experimental small molecule assays from the popular ChemBank database and show that screening data are widely affected by both assay-specific and plate-specific spatial biases. Importantly, the bias affecting screening data can fit an additive or multiplicative model. We show that the use of appropriate statistical methods is essential for improving the quality of experimental screening data. The presented methodology can be recommended for the analysis of current and next-generation screening data.
| S-EPMC5607347 | biostudies-literature
Project description:DFO - beads in situ microcosms Raw sequence reads
| PRJNA629517 | ENA
Project description:BACKGROUND: The advent of high-throughput and cost-effective genotyping platforms made genome-wide association (GWA) studies a reality. While the primary focus has been invested upon the improvement of reducing genotyping error, the problems associated with missing calls are largely overlooked. RESULTS: To probe into the effect of missing calls on GWAs, we demonstrated experimentally the prevalence and severity of the problem of missing call bias (MCB) in four genotyping technologies (Affymetrix 500 K SNP array, SNPstream, TaqMan, and Illumina Beadlab). Subsequently, we showed theoretically that MCB leads to biased conclusions in the subsequent analyses, including estimation of allele/genotype frequencies, the measurement of HWE and association tests under various modes of inheritance relationships. We showed that MCB usually leads to power loss in association tests, and such power change is greater than what could be achieved by equivalent reduction of sample size unbiasedly. We also compared the bias in allele frequency estimation and in association tests introduced by MCB with those by genotyping errors. Our results illustrated that in most cases, the bias can be greatly reduced by increasing the call-rate at the cost of genotyping error rate. CONCLUSION: The commonly used 'no-call' procedure for the observations of borderline quality should be modified. If the objective is to minimize the bias, the cut-off for call-rate and that for genotyping error rate should be properly coupled in GWA. We suggested that the ongoing QC cut-off for call-rate should be increased, while the cut-off for genotyping error rate can be reduced properly.
| S-EPMC2670840 | biostudies-literature