Dataset Information

BAMClipper: removing primers from alignments to minimize false-negative mutations in amplicon next-generation sequencing.

ABSTRACT: Amplicon-based next-generation sequencing (NGS) has been widely adopted for genetic variation detection in human and other organisms. Conventional data analysis paradigm includes primer trimming before read mapping. Here we introduce BAMClipper that removes primer sequences after mapping original sequencing reads by soft-clipping SAM/BAM alignments. Mutation detection accuracy was affected by the choice of primer handling approach based on real NGS datasets of 7 human peripheral blood or breast cancer tissue samples with known BRCA1/BRCA2 mutations and >130000 simulated NGS datasets with unique mutations. BAMClipper approach detected a BRCA1 deletion (c.1620_1636del) that was otherwise missed due to edge effect. Simulation showed high false-negative rate when primers were perfectly trimmed as in conventional practice. Among the other 6 samples, variant allele frequencies of 5 BRCA1/BRCA2 mutations (indel or single-nucleotide variants) were diluted by apparently wild-type primer sequences from an overlapping amplicon (17 to 82% under-estimation). BAMClipper was robust in both situations and all 7 mutations were detected. When compared with Cutadapt, BAMClipper was faster and maintained equally high primer removal effectiveness. BAMClipper is implemented in Perl and is available under an open source MIT license at https://github.com/tommyau/bamclipper.

SUBMITTER: Au CH

PROVIDER: S-EPMC5431517 | biostudies-literature | 2017 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

BAMClipper: removing primers from alignments to minimize false-negative mutations in amplicon next-generation sequencing.

Au Chun Hang CH Ho Dona N DN Kwong Ava A Chan Tsun Leung TL Ma Edmond S K ESK

Scientific reports 20170508 1

Amplicon-based next-generation sequencing (NGS) has been widely adopted for genetic variation detection in human and other organisms. Conventional data analysis paradigm includes primer trimming before read mapping. Here we introduce BAMClipper that removes primer sequences after mapping original sequencing reads by soft-clipping SAM/BAM alignments. Mutation detection accuracy was affected by the choice of primer handling approach based on real NGS datasets of 7 human peripheral blood or breast ...[more]

PMID: 28484262

Similar Datasets

Project description:COVID-19 is caused by SARS-CoV-2, several virulent variants of which have emerged since 2019. More than 529 million people have been infected, and at least 6 million have died. Our aim was to develop a fast, accurate, low-cost method for detecting and identifying newly emerging variants of concern (VOCs) that could pose a global threat. The 341-bp DNA sequence of a specific region of the SARS-CoV-2's spike protein was amplified by a one-step PCR on RNA samples from 46 patients. The product was sequenced using next-generation sequencing (NGS). DNA sequences from seven genomes, the original Wuhan isolate and six different representative variants obtained from the GISAID website, were used as references. Complete whole-genome sequences from local isolates were also obtained from the GISAID website, and their RNA was used for comparison. We used an amplicon-based NGS method (termed VOC-NGS) for genotyping and successfully identified all 46 samples. Fifteen (32.6%) were like the original isolate. Twenty-seven were VOCs: nine (19.5%) Alpha, eight (19%) Delta, six (14%) Beta, and four (8.7%) Omicron. Two were variants of interest (VOI): one (2%) Kappa and one (2%) Zeta. Two samples were mixtures of two variants, one of Alpha and Beta and one of Alpha and Delta. The Spearman correlation between whole-genome sequencing (WGS) and VOC-NGS was significant (P < 0.001) with perfect agreement (Kappa = 0.916) for 36/38 (94.7%) samples with VOC-NGS detecting all the known VOCs. Genotyping by VOC-NGS enables rapid screening of high-throughput clinical samples that includes the identification of VOCs and mixtures of variants, at lower cost than WGS. IMPORTANCE The manuscript described SARS-Cov-2 genotyping by VOC-NGS, which presents an ideal balance of accuracy, rapidity, and cost for detecting and globally tracking VOCs and some VOI of SARS-CoV-2. A large number of clinical samples can be tested together. Rapid introduction of new mutations at a specific site of the spike protein necessitates efficient strain detection and identification to enable choice of treatment and the application of vaccination, as well as planning public health policy.

Project description:BackgroundNext-generation sequencing (NGS) has enabled efficient high-resolution typing of human leukocyte antigen (HLA) genes with minimal ambiguity. Most commercially available assays amplify individual or subgroup of HLA genes by long-range PCR followed by library preparation and sequencing. The AllType assay simplifies the workflow by amplifying 11 transplant-relevant HLA genes in one PCR reaction. Here, we report the performance of this unique workflow evaluated using 218 genetically diverse samples.MethodsFive whole genes (HLA-A/B/C/DQA1/DPA1) and six near-whole genes (HLA-DRB1/DRB345/DQB1/DPB1; excluding exon 1 and part of intron 1) were amplified in a multiplexed, long-range PCR. Manual library preparation was performed per manufacturer's protocol, followed by template preparation and chip loading on the Ion Chef, and sequencing on the Ion S5 sequencer. Pre-specified rules for quality control and repeat testing were followed; technologists were blinded to the reference results. The concordance between AllType and reference results was determined at 2-field resolution. We also describe the ranges of input DNA and library concentrations, read number per sample and per locus, and key health metrics in relation to typing results.ResultsThe concordance rates were 98.6%, 99.8% and 99.9% at the sample (n = 218), genotype (n = 1688), and allele (n = 3376) levels, respectively. Three genotypes were discordant, all of which shared the same G group typing results with the reference. Most ambiguous genotypes (116 out of 144, 80.6%) were due to the lack of exon 1 and intron 1 coverage for HLA-DRB1/DRB345/DQB1/DPB1 genes. A broad range of input DNA concentrations and library concentrations were tolerated. Per sample read numbers were adequate for accurate genotyping. Per locus read numbers showed some inter-lot variations, and a trend toward improved inter-locus balance was observed with later lots of reagents.ConclusionThe AllType assay on the Ion Chef/Ion S5 platform offers a robust and efficient workflow for clinical HLA typing at the 2-field resolution. The multiplex PCR strategy simplifies the laboratory procedure without compromising the typing accuracy.

Dataset Information

BAMClipper: removing primers from alignments to minimize false-negative mutations in amplicon next-generation sequencing.

Publications

BAMClipper: removing primers from alignments to minimize false-negative mutations in amplicon next-generation sequencing.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets