Fiona Hyland

Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding

We describe the genome sequencing of an anonymous individual of African origin using a novel ligation-based sequencing assay that enables a unique form of error correction that improves the raw accuracy of the aligned reads to >99.9%, allowing us to accurately call SNPs with as few as two reads per allele. We collected several billion mate-paired reads yielding approximately 18x haploid coverage of aligned sequence and close to 300x clone coverage. Over 98% of the reference genome is covered with at least one uniquely placed read, and 99.65% is spanned by at least one uniquely placed mate-paired clone. We identify over 3.8 million SNPs, 19% of which are novel. Mate-paired data are used to physically resolve haplotype phases of nearly two-thirds of the genotypes obtained and produce phased segments of up to 215 kb. We detect 226,529 intra-read indels, 5590 indels between mate-paired reads, 91 inversions, and four gene fusions. We use a novel approach for detecting indels between mate-paired reads that are smaller than the standard deviation of the insert size of the library and discover deletions in common with those detected with our intra-read approach. Dozens of mutations previously described in OMIM and hundreds of nonsynonymous single-nucleotide and structural variants in genes previously implicated in disease are identified in this individual. There is more genetic variation in the human genome still to be uncovered, and we provide guidance for future surveys in populations and cancer biopsies.

Extensive sequencing of seven human genomes to characterize benchmark reference materials.

The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

Development and Validation of a Scalable Next-Generation Sequencing System for Assessing Relevant Somatic Variants in Solid Tumors

Next-generation sequencing (NGS) has enabled genome-wide personalized oncology efforts at centers and companies with the specialty expertise and infrastructure required to identify and prioritize actionable variants. Such approaches are not scalable, preventing widespread adoption. Likewise, most targeted NGS approaches fail to assess key relevant genomic alteration classes. To address these challenges, we predefined the catalog of relevant solid tumor somatic genome variants (gain-of-function or loss-of-function mutations, high-level copy number alterations, and gene fusions) through comprehensive bioinformatics analysis of >700,000 samples. To detect these variants, we developed the Oncomine Comprehensive Panel (OCP), an integrative NGS-based assay [compatible with < 20 ng of DNA/RNA from formalin-fixed paraffin-embedded (FFPE) tissues], coupled with an informatics pipeline to specifically identify relevant predefined variants and created a knowledge base of related potential treatments, current practice guidelines, and open clinical trials. We validated OCP using molecular standards and more than 300 FFPE tumor samples, achieving >95% accuracy for KRAS, epidermal growth factor receptor, and BRAF mutation detection as well as for ALK and TMPRSS2:ERG gene fusions. Associating positive variants with potential targeted treatments demonstrated that 6% to 42% of profiled samples (depending on cancer type) harbored alterations beyond routine molecular testing that were associated with approved or guideline-referenced therapies. As a translational research tool, OCP identified adaptive CTNNB1 amplifications/mutations in treated prostate cancers. Through predefining somatic variants in solid tumors and compiling associated potential treatment strategies, OCP represents a simplified, broadly applicable targeted NGS system with the potential to advance precision oncology efforts.

Good laboratory practice for clinical next-generation sequencing informatics pipelines

Amy S Gargis, Centers for Disease Control & Prevention Lisa Kalman, Centers for Disease Control & Prevention David P Bick, Medical College of Wisconsin Cristina da Silva, Emory University David P Dimmock, Medical College of Wisconsin Birgit H Funke, Partners Healthcare Personalized Medicine Sivakumar Gowrisankar, Partners Healthcare Personalized Medicine Madhuri Hegde, Emory University Shashikant Kulkarni, Washington University Christopher E Mason, Cornell University

RNA-Seq Mapping and Detection of Gene Fusions with a Suffix Array Algorithm

High-throughput RNA sequencing enables quantification of transcripts (both known and novel), exon/exon junctions and fusions of exons from different genes. Discovery of gene fusions–particularly those expressed with low abundance– is a challenge with short- and medium-length sequencing reads. To address this challenge, we implemented an RNA-Seq mapping pipeline within the LifeScope software. We introduced new features including filter and junction mapping, annotation-aided pairing rescue and accurate mapping quality values. We combined this pipeline with a Suffix Array Spliced Read (SASR) aligner to detect chimeric transcripts. Performing paired-end RNA-Seq of the breast cancer cell line MCF-7 using the SOLiD system, we called 40 gene fusions among over 120,000 splicing junctions. We validated 36 of these 40 fusions with TaqMan assays, of which 25 were expressed in MCF-7 but not the Human Brain Reference. An intra-chromosomal gene fusion involving the estrogen receptor alpha gene ESR1, and another involving the RPS6KB1 (Ribosomal protein S6 kinase beta-1) were recurrently expressed in a number of breast tumor cell lines and a clinical tumor sample.

Expanding data and resources for forensic use of SNPs in individual identification

The potential value of SNPs for individual identification has been recognized by many researchers and different panels have been proposed. Here we present a new interface in the ALFRED database to access compendia of allele frequencies for several published panels of markers for forensic uses. One of those is our panel of individual identification SNPs (IISNPs) based on samples of 44 populations originating from many parts of the world. Here we also present additional data and additional statistical analyses that continue to support the value of our panel of IISNPs as a universal panel. We also describe initial developments of multiplex methods and various robustness analyses for our 45 marker IISNP panel.

Patterns of linkage disequilibrium between SNPs in a sardinian population isolate and the selection of markers for association studies

Objective: In isolated populations, ‘background’ linkage disequilibrium (LD) has been shown to extend over large genetic distances. This and their reduced environmental and genetic heterogeneity has stimulated interest in their potential for association mapping. We compared LD unit map distances with pair-wise measurements of LD in a dense single nucleotide polymorphism (SNP) set. Methods: We genotyped 771 SNPs in an 8 Mb segment of chromosome 22 on 101 individuals from the isolated village of Talana, Sardinia, and compared with outbred European populations. Results: Heterozygosity was remarkably similar in both populations. In contrast, the extent of LD observed was quite different. The decay of LD with distance is slower in the isolate. The differences in LD map lengths suggest that useful LD extends up to three times farther in the Sardinian population; smaller differences are seen with pairwise LD metrics. While LD map length slightly decreases with average relatedness, cryptic relatedness does not explain the decrease in LD map length. Haplotypes, block boundaries, and patterns of LD are similar in both populations, suggesting a shared distribution of recombination hotspots. Conclusions: About 15% fewer haplotype tagging SNPs need to be genotyped in the isolate, and possibly 70% fewer if selecting SNPs evenly spaced on the metric LD map.

Analytical Validation of a Next-Generation Sequencing Assay to Monitor Immune Responses in Solid Tumors

We have developed a next-generation sequencing assay to quantify biomarkers of the host immune response in formalin-fixed, paraffin-embedded (FFPE) tumor specimens. This assay aims to provide clinicians with a comprehensive characterization of the immunologic tumor microenvironment as a guide for therapeutic decisions on patients with solid tumors. The assay relies on RNA-sequencing (seq) to semiquantitatively measure the levels of 43 transcripts related to anticancer immune responses and 11 transcripts that reflect the relative abundance of tumor-infiltrating lymphocytes, as well as on DNA-seq to estimate mutational burden. The assay has a clinically relevant 5-day turnaround time and can be conducted on as little as 2.5 ng of RNA and 1.8 ng of genomic DNA extracted from three to five standard FFPE sections. The standardized next-generation sequencing workflow produced sequencing reads adequate for clinical testing of matched RNA and DNA from several samples in a single run. Assay performance for gene-specific sensitivity, linearity, dynamic range, and detection threshold was estimated across a wide range of actual and artificial FFPE samples selected or generated to address preanalytical variability linked to specimen features (eg, tumor-infiltrating lymphocyte abundance, percentage of necrosis), and analytical variability linked to assay features (eg, batch size, run, day, operator). Analytical precision studies demonstrated that the assay is highly reproducible and accurate compared with established orthogonal approaches.

Validation of the performance of a comprehensive genotyping assay panel of single nucleotide polymorphisms in drug metabolism enzyme genes

A class of genes, known as drug metabolism enzymes (DMEs) are responsible for the metabolism and transport of drugs and other xenobiotics. Variation in DME genes most likely accounts for a proportion of the variability in drug response in humans, and may contribute to complex diseases such as cancer (Nebert DW, Dieter MZ. Pharmacology 2000;61:124–135). To date, assessing the extent of this variation has proven difficult, especially because of sequence paralogy issues that cause difficulty when attempting to genotype polymorphisms in very closely‐related gene families (Murphy MP. Pharmacogenomics 2000;1:115–123; Ingelman‐Sundberg M. Drug Metab Rev 1999;31:449–459). We have developed and genotyped a panel of N=2,325 individual TaqMan® genotyping assays for polymorphisms in >200 DME genes; many of the variants in the panel are single nucleotide polymorphisms (SNPs) that are of known or putative function (e.g., missense, nonsense or frameshift). Using these assays, we have examined genetic variation among several groups of populations, including: 1) the two SNP500 Cancer population panels (http://snp500cancer.nci.nih.gov; last accessed: 11 December 2007); and 2) the panel used in the International HapMap Project panel (www.hapmap.org; last accessed: 11 December 2007). We have developed a comprehensive validation strategy to ensure reproducibility and accuracy of the assays and estimated minor allele frequencies. Here, we present the results of these analyses, which strongly suggest that this panel of DME assays are of extremely high quality and produce robust, accurate, and reproducible results. Hum Mutat 29(5), 750–756, 2008. Published 2008, Wiley‐Liss, Inc.

Principles and Recommendations for Standardizing the Use of the Next-Generation Sequencing Variant File in Clinical Settings

A national workgroup convened by the Centers for Disease Control and Prevention identified principles and made recommendations for standardizing the description of sequence data contained within the variant file generated during the course of clinical next-generation sequence analysis for diagnosing human heritable conditions. The specifications for variant files were initially developed to be flexible with regard to content representation to support a variety of research applications. This flexibility permits variation with regard to how sequence findings are described and this depends, in part, on the conventions used. For clinical laboratory testing, this poses a problem because these differences can compromise the capability to compare sequence findings among laboratories to confirm results and to query databases to identify clinically relevant variants. To provide for a more consistent representation of sequence findings described within variant files, the workgroup made several recommendations that considered alignment to a common reference sequence, variant caller settings, use of genomic coordinates, and gene and variant naming conventions. These recommendations were considered with regard to the existing variant file specifications presently used in the clinical setting. Adoption of these recommendations is anticipated to reduce the potential for ambiguity in describing sequence findings and facilitate the sharing of genomic data among clinical laboratories and other entities.