Dataset Information

Medical implications of technical accuracy in genome sequencing.

ABSTRACT: BACKGROUND:As whole exome sequencing (WES) and whole genome sequencing (WGS) transition from research tools to clinical diagnostic tests, it is increasingly critical for sequencing methods and analysis pipelines to be technically accurate. The Genome in a Bottle Consortium has recently published a set of benchmark SNV, indel, and homozygous reference genotypes for the pilot whole genome NIST Reference Material based on the NA12878 genome. METHODS:We examine the relationship between human genome complexity and genes/variants reported to be associated with human disease. Specifically, we map regions of medical relevance to benchmark regions of high or low confidence. We use benchmark data to assess the sensitivity and positive predictive value of two representative sequencing pipelines for specific classes of variation. RESULTS:We observe that the accuracy of a variant call depends on the genomic region, variant type, and read depth, and varies by analytical pipeline. We find that most false negative WGS calls result from filtering while most false negative WES variants relate to poor coverage. We find that only 74.6% of the exonic bases in ClinVar and OMIM genes and 82.1% of the exonic bases in ACMG-reportable genes are found in high-confidence regions. Only 990 genes in the genome are found entirely within high-confidence regions while 593 of 3,300 ClinVar/OMIM genes have less than 50% of their total exonic base pairs in high-confidence regions. We find greater than 77 % of the pathogenic or likely pathogenic SNVs currently in ClinVar fall within high-confidence regions. We identify sites that are prone to sequencing errors, including thousands present in publicly available variant databases. Finally, we examine the clinical impact of mandatory reporting of secondary findings, highlighting a false positive variant found in BRCA2. CONCLUSIONS:Together, these data illustrate the importance of appropriate use and continued improvement of technical benchmarks to ensure accurate and judicious interpretation of next-generation DNA sequencing results in the clinical setting.

SUBMITTER: Goldfeder RL

PROVIDER: S-EPMC4774017 | biostudies-literature | 2016 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Medical implications of technical accuracy in genome sequencing.

Goldfeder Rachel L RL Priest James R JR Zook Justin M JM Grove Megan E ME Waggott Daryl D Wheeler Matthew T MT Salit Marc M Ashley Euan A EA

Genome medicine 20160302 1

<h4>Background</h4>As whole exome sequencing (WES) and whole genome sequencing (WGS) transition from research tools to clinical diagnostic tests, it is increasingly critical for sequencing methods and analysis pipelines to be technically accurate. The Genome in a Bottle Consortium has recently published a set of benchmark SNV, indel, and homozygous reference genotypes for the pilot whole genome NIST Reference Material based on the NA12878 genome.<h4>Methods</h4>We examine the relationship betwee ...[more]

PMID: 26932475

Similar Datasets

Project description:IMPORTANCE:Whole-genome sequencing (WGS) is increasingly applied in clinical medicine and is expected to uncover clinically significant findings regardless of sequencing indication. OBJECTIVES:To examine coverage and concordance of clinically relevant genetic variation provided by WGS technologies; to quantitate inherited disease risk and pharmacogenomic findings in WGS data and resources required for their discovery and interpretation; and to evaluate clinical action prompted by WGS findings. DESIGN, SETTING, AND PARTICIPANTS:An exploratory study of 12 adult participants recruited at Stanford University Medical Center who underwent WGS between November 2011 and March 2012. A multidisciplinary team reviewed all potentially reportable genetic findings. Five physicians proposed initial clinical follow-up based on the genetic findings. MAIN OUTCOMES AND MEASURES:Genome coverage and sequencing platform concordance in different categories of genetic disease risk, person-hours spent curating candidate disease-risk variants, interpretation agreement between trained curators and disease genetics databases, burden of inherited disease risk and pharmacogenomic findings, and burden and interrater agreement of proposed clinical follow-up. RESULTS:Depending on sequencing platform, 10% to 19% of inherited disease genes were not covered to accepted standards for single nucleotide variant discovery. Genotype concordance was high for previously described single nucleotide genetic variants (99%-100%) but low for small insertion/deletion variants (53%-59%). Curation of 90 to 127 genetic variants in each participant required a median of 54 minutes (range, 5-223 minutes) per genetic variant, resulted in moderate classification agreement between professionals (Gross ?, 0.52; 95% CI, 0.40-0.64), and reclassified 69% of genetic variants cataloged as disease causing in mutation databases to variants of uncertain or lesser significance. Two to 6 personal disease-risk findings were discovered in each participant, including 1 frameshift deletion in the BRCA1 gene implicated in hereditary breast and ovarian cancer. Physician review of sequencing findings prompted consideration of a median of 1 to 3 initial diagnostic tests and referrals per participant, with fair interrater agreement about the suitability of WGS findings for clinical follow-up (Fleiss ?, 0.24; P?<?001). CONCLUSIONS AND RELEVANCE:In this exploratory study of 12 volunteer adults, the use of WGS was associated with incomplete coverage of inherited disease genes, low reproducibility of detection of genetic variation with the highest potential clinical effects, and uncertainty about clinically reportable findings. In certain cases, WGS will identify clinically actionable genetic variants warranting early medical intervention. These issues should be considered when determining the role of WGS in clinical medicine.

Project description:Although microviruses do not possess a visible tail structure, one vertex rearranges after interacting with host lipopolysaccharides. Most examinations of host range, eclipse, and penetration were conducted before this "host-induced" unique vertex was discovered and before DNA sequencing became routine. Consequently, structure-function relationships dictating host range remain undefined. Biochemical and genetic analyses were conducted with two closely related microviruses, α3 and ST-1. Despite ∼90% amino acid identity, the natural host of α3 is Escherichia coli C, whereas ST-1 is a K-12-specific phage. Virions attached and eclipsed to both native and unsusceptible hosts; however, they breached only the native host's cell wall. This suggests that unsusceptible host-phage interactions promote off-pathway reactions that can inactivate viruses without penetration. This phenomenon may have broader ecological implications. To determine which structural proteins conferred host range specificity, chimeric virions were generated by individually interchanging the coat, spike, or DNA pilot proteins. Interchanging the coat protein switched host range. However, host range expansion could be conferred by single point mutations in the coat protein. The expansion phenotype was recessive: genetically mutant progeny from coinfected cells did not display the phenotype. Thus, mutant isolation required populations generated in environments with low multiplicities of infection (MOI), a phenomenon that may have impacted past host range studies in both prokaryotic and eukaryotic systems. The resulting genetic and structural data were consistent enough that host range expansion could be predicted, broadening the classical definition of antireceptors to include interfaces between protein complexes within the capsid.IMPORTANCE To expand host range, viruses must interact with unsusceptible host cell surfaces, which could be detrimental. As observed in this study, virions were inactivated without genome penetration. This may be advantageous to potential new hosts, culling the viral population from which an expanded host range mutant could emerge. When identified, altered host range mutations were recessive. Accordingly, isolation required populations generated in low-MOI environments. However, in laboratory settings, viral propagation includes high-MOI conditions. Typically, infected cultures incubate until all cells produce progeny. Thus, coinfections dominate later replication cycles, masking recessive host range expansion phenotypes. This may have impacted similar studies with other viruses. Last, structural and genetic data could be used to predict site-directed mutant phenotypes, which may broaden the classic antireceptor definition to include interfaces between capsid complexes.

Project description:MotivationShort-read whole-genome sequencing (WGS) is a vital tool for clinical applications and basic research. Genetic divergence from the reference genome, repetitive sequences and sequencing bias reduces the performance of variant calling using short-read alignment, but the loss in recall and specificity has not been adequately characterized. To benchmark short-read variant calling, we used 36 diverse clinical Mycobacterium tuberculosis (Mtb) isolates dually sequenced with Illumina short-reads and PacBio long-reads. We systematically studied the short-read variant calling accuracy and the influence of sequence uniqueness, reference bias and GC content.ResultsReference-based Illumina variant calling demonstrated a maximum recall of 89.0% and minimum precision of 98.5% across parameters evaluated. The approach that maximized variant recall while still maintaining high precision (<99%) was tuning the mapping quality filtering threshold, i.e. confidence of the read mapping (recall = 85.8%, precision = 99.1%, MQ ≥ 40). Additional masking of repetitive sequence content is an alternative conservative approach to variant calling that increases precision at cost to recall (recall = 70.2%, precision = 99.6%, MQ ≥ 40). Of the genomic positions typically excluded for Mtb, 68% are accurately called using Illumina WGS including 52/168 PE/PPE genes (34.5%). From these results, we present a refined list of low confidence regions across the Mtb genome, which we found to frequently overlap with regions with structural variation, low sequence uniqueness and low sequencing coverage. Our benchmarking results have broad implications for the use of WGS in the study of Mtb biology, inference of transmission in public health surveillance systems and more generally for WGS applications in other organisms.Availability and implementationAll relevant code is available at https://github.com/farhat-lab/mtb-illumina-wgs-evaluation.Supplementary informationSupplementary data are available at Bioinformatics online.

Dataset Information

Medical implications of technical accuracy in genome sequencing.

Publications

Medical implications of technical accuracy in genome sequencing.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets