Dataset Information

Generation of artificial FASTQ files to evaluate the performance of next-generation sequencing pipelines.

ABSTRACT: Pipelines for the analysis of Next-Generation Sequencing (NGS) data are generally composed of a set of different publicly available software, configured together in order to map short reads of a genome and call variants. The fidelity of pipelines is variable. We have developed ArtificialFastqGenerator, which takes a reference genome sequence as input and outputs artificial paired-end FASTQ files containing Phred quality scores. Since these artificial FASTQs are derived from the reference genome, it provides a gold-standard for read-alignment and variant-calling, thereby enabling the performance of any NGS pipeline to be evaluated. The user can customise DNA template/read length, the modelling of coverage based on GC content, whether to use real Phred base quality scores taken from existing FASTQ files, and whether to simulate sequencing errors. Detailed coverage and error summary statistics are outputted. Here we describe ArtificialFastqGenerator and illustrate its implementation in evaluating a typical bespoke NGS analysis pipeline under different experimental conditions. ArtificialFastqGenerator was released in January 2012. Source code, example files and binaries are freely available under the terms of the GNU General Public License v3.0. from https://sourceforge.net/projects/artfastqgen/.

SUBMITTER: Frampton M

PROVIDER: S-EPMC3495771 | biostudies-literature | 2012

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Generation of artificial FASTQ files to evaluate the performance of next-generation sequencing pipelines.

Frampton Matthew M Houlston Richard R

PloS one 20121112 11

Pipelines for the analysis of Next-Generation Sequencing (NGS) data are generally composed of a set of different publicly available software, configured together in order to map short reads of a genome and call variants. The fidelity of pipelines is variable. We have developed ArtificialFastqGenerator, which takes a reference genome sequence as input and outputs artificial paired-end FASTQ files containing Phred quality scores. Since these artificial FASTQs are derived from the reference genome, ...[more]

PMID: 23152858

Similar Datasets

Project description:This repository contains all the FASTQ files for the five data modalities (scRNA-seq, scATAC-seq, Multiome, CITE-seq+scVDJ-seq, and spatial transcriptomics) used in the article \\"An Atlas of Cells in The Human Tonsil,\\" published in Immunity in 2024. Inspired by the TCGA barcodes, we have named each fastq file with the following convention: [TECHNOLOGY].[DONOR_ID].[SUBPROJECT].[GEM_ID].[LIBRARY_ID].[LIBRARY_TYPE].[LANE].[READ].fastq.gz which allows to retrieve all metadata from the name itself. Here is a full description of each field: - TECHNOLOGY: scRNA-seq, scATAC-seq, Multiome, CITE-seq+scVDJ-seq, and spatial transcriptomics (Visium). We also include the fastq files associated with the multiome experiments performed on two mantle cell lymphoma patients (MCL). - DONOR_ID: identifier for each of the 17 patients included in the cohort. We provide the donor-level metadata in the file \\"tonsil_atlas_donor_metadata.csv\\", including the hospital, sex, age, age group, cause for tonsillectomy and cohort type for every donor. - SUBPROJECT: each subproject corresponds to one run of the 10x Genomics Chromium™ Chip. - GEM_ID: each run of the 10x Genomics Chromium™ Chip consists of up to 8 \\"GEM wells\\" (see https://www.10xgenomics.com/support/software/cell-ranger/getting-started/cr-glossary): a set of partitioned cells (Gel Beads-in-emulsion) from a single 10x Genomics Chromium™ Chip channel. We give a unique identifier to each of these channels. - LIBRARY_ID: one or more sequencing libraries can be derived from a GEM well. For instance, multiome yields two libraries (ATAC and RNA) and CITE-seq+scVDJ yields 4 libraries (RNA, ADT, BCR, TCR). - LIBRARY_TYPE: the type of library for each library_id. Note that we used cell hashing () for a subset of the scRNA-seq libraries, and thus the library_type can be \\"not_hashed\\", \\"hashed_cdna\\" (RNA expression) or \\"hashed_hto\\" (the hashtag oligonucleotides). - LANE: to increase sequencing depth, each library was sequenced in more than one lane. Important: all lanes corresponding to the same sequencing library need to be inputed together to cellranger, because they come from the same set of cells. - READ: for scATAC-seq we have three reads (R1, R2 or R3), see cellranger-atac's documentation. While we find these names to be the most useful, they need to be changed to follow cellranger's conventions. We provide a code snippet in the README file of the GitHub repository associated with the tonsil atlas to convert between both formats (https://github.com/Single-Cell-Genomics-Group-CNAG-CRG/TonsilAtlas/). Besides the fastq files, cellranger (and other mappers) require additional files, which we also provide in this repository: - cell_hashing_metadata.csv: as mentioned above, we ran cell hashing (10.1186/s13059-018-1603-1) to detect doublets and reduce cost per cell. This file provides the sequence of the hashtag oligonucleotides in cellranger convention to allow demultiplexing. - cite_seq_feature_reference.csv: similar to the previous file, this one links each protein surface marker to the hashtag oligonucleotide that identified it in the CITE-seq experiment. - V10M16-059.gpr and V19S23-039.gpr: these correspond to the two slides of the two Visium experiments performed in the tonsil atlas. They are needed to run spaceranger. - [GEM_ID]_[SLIDE]_[CAPTURE_AREA].jpg: 8 images associated with the Visium experiments. Here, GEM_ID refers to each of the 4 capture areas in each slide. - [TECHNOLOGY]_sequencing_metadata.csv: the GEM-level metadata for each technology. It includes the relationship between subproject, gem_id, library_id, library_type and donor_id. These are the other repositories associated with the tonsil atlas: - Expression and accessibility matrices: https://zenodo.org/records/10373041 - Seurat objects: https://zenodo.org/records/8373756 - HCATonsilData package: https://bioconductor.org/packages/release/data/experiment/html/HCATonsilData.html - Azimuth: https://azimuth.hubmapconsortium.org/ - Github: https://github.com/Single-Cell-Genomics-Group-CNAG-CRG/TonsilAtlas

Project description:Single nucleotide polymorphisms (SNPs) are widely used in genome-wide association studies and population genetics analyses. Next-generation sequencing (NGS) has become convenient, and many SNP-calling pipelines have been developed for human NGS data. We took advantage of a gap knowledge in selecting the appropriated SNP calling pipeline to handle with high-throughput NGS data. To fill this gap, we studied and compared seven SNP calling pipelines, which include 16GT, genome analysis toolkit (GATK), Bcftools-single (Bcftools single sample mode), Bcftools-multiple (Bcftools multiple sample mode), VarScan2-single (VarScan2 single sample mode), VarScan2-multiple (VarScan2 multiple sample mode) and Freebayes pipelines, using 96 NGS data with the different depth gradients of approximately 5X, 10X, 20X, 30X, 40X, and 50X coverage from 16 Rhode Island Red chickens. The sixteen chickens were also genotyped with a 50K SNP array, and the sensitivity and specificity of each pipeline were assessed by comparison to the results of SNP arrays. For each pipeline, except Freebayes, the number of detected SNPs increased as the input read depth increased. In comparison with other pipelines, 16GT, followed by Bcftools-multiple, obtained the most SNPs when the input coverage exceeded 10X, and Bcftools-multiple obtained the most when the input was 5X and 10X. The sensitivity and specificity of each pipeline increased with increasing input. Bcftools-multiple had the highest sensitivity numerically when the input ranged from 5X to 30X, and 16GT showed the highest sensitivity when the input was 40X and 50X. Bcftools-multiple also had the highest specificity, followed by GATK, at almost all input levels. For most calling pipelines, there were no obvious changes in SNP numbers, sensitivities or specificities beyond 20X. In conclusion, (1) if only SNPs were detected, the sequencing depth did not need to exceed 20X; (2) the Bcftools-multiple may be the best choice for detecting SNPs from chicken NGS data, but for a single sample or sequencing depth greater than 20X, 16GT was recommended. Our findings provide a reference for researchers to select suitable pipelines to obtain SNPs from the NGS data of chickens or nonhuman animals.

Dataset Information

Generation of artificial FASTQ files to evaluate the performance of next-generation sequencing pipelines.

Publications

Generation of artificial FASTQ files to evaluate the performance of next-generation sequencing pipelines.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets