Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Clinical interpretation of whole-genome and whole-transcriptome sequencing for precision oncology.

Whole-genome sequencing either alone or in combination with whole-transcriptome sequencing has started to be used to analyze clinical tumor samples to improve diagnosis, provide risk stratification, and select patient-specific therapies. Compared with current genomic testing strategies, largely focused on small number of genes tested individually or targeted panels, whole-genome and transcriptome sequencing (WGTS) provides novel opportunities to identify and report a potentially much larger number of actionable alterations with diagnostic, prognostic, and/or predictive impact. Such alterations include point mutations, indels, copy- number aberrations and structural variants, but also germline variants, fusion genes, noncoding alterations and mutational signatures. Nevertheless, these comprehensive tests are accompanied by many challenges ranging from the extent and diversity of sequence alterations detected by these methods to the complexity and limited existing standardization in interpreting them. We describe the challenges of WGTS interpretation and the opportunities with comprehensive genomic testing.

Analytical demands to use whole-genome sequencing in precision oncology.

Interrogating the tumor genome in its entirety by whole-genome sequencing (WGS) offers an unprecedented insight into the biology and pathogenesis of cancer, with potential impact on diagnostics, prognostication and therapy selection. WGS is able to detect sequence as well as structural variants and thereby combines central domains of cytogenetics and molecular genetics. Given the potential of WGS in directing targeted therapeutics and clinical decision-making, we envision a gradual transition of the method from research to clinical routine. This review is one out of three within this issue aimed at facilitating this effort, by discussing in-depth analytical validation, clinical interpretation and clinical utility of WGS. The review highlights the requirements for implementing, validating and maintaining a clinical WGS pipeline to obtain high-quality patient-specific data in accordance with the local regulatory landscape. Every step of the WGS pipeline, which includes DNA extraction, library preparation, sequencing, bioinformatics analysis, and data storage, is considered with respect to its logistics, necessities, potential pitfalls, and the required quality management. WGS is likely to drive clinical diagnostics and patient care forward, if requirements and challenges of the technique are recognized and met.

Pancreatic intraductal tubulopapillary neoplasm is genetically distinct from intraductal papillary mucinous neoplasm and ductal adenocarcinoma.

Intraductal tubulopapillary neoplasm is a relatively recently described member of the pancreatic intraductal neoplasm family. The more common member of this family, intraductal papillary mucinous neoplasm, often carries genetic alterations typical of pancreatic infiltrating ductal adenocarcinoma (KRAS, TP53, and CDKN2A) but additionally has mutations in GNAS and RNF43 genes. However, the genetic characteristics of intraductal tubulopapillary neoplasm have not been well characterized. Twenty-two intraductal tubulopapillary neoplasms were analyzed by either targeted next-generation sequencing, which enabled the identification of sequence mutations, copy number alterations, and selected structural rearrangements involving all targeted (≥300) genes, or whole-exome sequencing. Three of these intraductal tubulopapillary neoplasms were also subjected to whole-genome sequencing. All intraductal tubulopapillary neoplasms revealed the characteristic histologic (cellular intraductal nodules of back-to-back tubular glands lined by predominantly cuboidal cells with atypical nuclei and no obvious intracellular mucin) and immunohistochemical (immunolabeled with MUC1 and MUC6 but were negative for MUC2 and MUC5AC) features. By genomic analyses, there was loss of CDKN2A in 5/20 (25%) of these cases. However, the majority of the previously reported intraductal papillary mucinous neoplasm-related alterations were absent. Moreover, in contrast to most ductal neoplasms of the pancreas, MAP-kinase pathway was not involved. In fact, 2/22 (9%) of intraductal tubulopapillary neoplasms did not reveal any mutations in the tested genes. However, certain chromatin remodeling genes (MLL1, MLL2, MLL3, BAP1, PBRM1, EED, and ATRX) were found to be mutated in 7/22 (32%) of intraductal tubulopapillary neoplasms and 27% harbored phosphatidylinositol 3-kinase (PI3K) pathway (PIK3CA, PIK3CB, INPP4A, and PTEN) mutations. In addition, 4/18 (18%) of intraductal tubulopapillary neoplasms had FGFR2 fusions (FGFR2-CEP55, FGFR2-SASS6, DISP1-FGFR2, FGFR2-TXLNA, and FGFR2-VCL) and 1/18 (5.5%) had STRN-ALK fusion. Intraductal tubulopapillary neoplasm is a distinct clinicopathologic entity in the pancreas. Although its intraductal nature and some clinicopathologic features resemble those of intraductal papillary mucinous neoplasm, our results suggest that intraductal tubulopapillary neoplasm has distinguishing genetic characteristics. Some of these mutated genes are potentially targetable. Future functional studies will be needed to determine the consequences of these gene alterations.

The oncocytic subtype is genetically distinct from other pancreatic intraductal papillary mucinous neoplasm subtypes.

In 2010, the World Health Organization reclassified the entity originally described as intraductal oncocytic papillary neoplasm as the 'oncocytic subtype' of intraductal papillary mucinous neoplasm. Although several key molecular alterations of other intraductal papillary mucinous neoplasm subtypes have been discovered, including common mutations in KRAS, GNAS, and RNF3, those of oncocytic subtype have not been well characterized. We analyzed 11 pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms. Nine pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms uniformly exhibited typical entity-defining morphology of arborizing papillae lined by layers of cells with oncocytic cytoplasm, prominent, nucleoli, and intraepithelial lumina. The remaining two were atypical. One lacked the arborizing papilla and had flat oncocytic epithelium only; the other one had focal oncocytic epithelium in a background of predominantly intestinal subtype intraductal papillary mucinous neoplasm. Different components of this case were analyzed separately. Formalin-fixed, paraffin-embedded specimens of all cases were microdissected and subjected to high-depth-targeted next-generation sequencing for a panel of 300 key cancer-associated genes in a platform that enabled the identification of sequence mutations, copy number alterations, and select structural rearrangements involving all targeted genes. Fresh frozen specimens of two cases were also subjected to whole-genome sequencing. For the nine typical pancreatic 'oncocytic subtype' of intraductal papillary mucinous neoplasms, the number of mutations per case, identified by next-generation sequencing, ranged from 1 to 10 (median=4). None of these cases had KRAS or GNAS mutations and only one had both RNF43 and PIK3R1 mutations. ARHGAP26, ASXL1, EPHA8, and ERBB4 genes were somatically altered in more than one of these typical 'oncocytic subtype' of intraductal papillary mucinous neoplasms but not in the other two atypical ones. In the neoplasm with flat oncocytic epithelium, the only mutated gene was KRAS. All components of the intestinal subtype intraductal papillary mucinous neoplasms with focal oncocytic epithelium manifested TP53, GNAS, and RNF43 mutations. In conclusion, this study elucidates that 'oncocytic subtype' of intraductal papillary mucinous neoplasm is not only morphologically distinct but also genetically distinct from other intraductal papillary mucinous neoplasm subtypes. Considering that now its biologic behavior is also being found to be different than other intraductal papillary mucinous neoplasm subtypes, 'oncocytic subtype' of intraductal papillary mucinous neoplasm warrants being recognized separately.

A quantitative proteomics-based signature of platinum sensitivity in ovarian cancer cell lines.

Although DNA encodes the molecular instructions that underlie the control of cell function, it is the proteins that are primarily responsible for implementing those instructions. Therefore quantitative analyses of the proteome would be expected to yield insights into important candidates for the detection and treatment of disease. We present an iTRAQ (isobaric tag for relative and absolute quantification)-based proteomic analysis of ten ovarian cancer cell lines and two normal ovarian surface epithelial cell lines. We profiled the abundance of 2659 cellular proteins of which 1273 were common to all 12 cell lines. Of the 1273, 75 proteins exhibited elevated expression and 164 proteins had diminished expression in the cancerous cells compared with the normal cell lines. The iTRAQ expression profiles allowed us to segregate cell lines based upon sensitivity and resistance to carboplatin. Importantly, we observed no substantial correlation between protein abundance and RNA expression or epigenetic DNA methylation data. Furthermore, we could not discriminate between sensitivity and resistance to carboplatin on the basis of RNA expression and DNA methylation data alone. The present study illustrates the importance of proteomics-based discovery for defining the basis for the carboplatin response in ovarian cancer and highlights candidate proteins, particularly involved in cellular redox regulation, homologous recombination and DNA damage repair, which otherwise could not have been predicted from whole genome and expression data sources alone.

Functional mutation of SMAC/DIABLO, encoding a mitochondrial proapoptotic protein, causes human progressive hearing loss DFNA64.

SMAC/DIABLO is a mitochondrial proapoptotic protein that is released from mitochondria during apoptosis and counters the inhibitory activities of inhibitor of apoptosis proteins, IAPs. By linkage analysis and candidate screening, we identified a heterozygous SMAC/DIABLO mutation, c.377C>T (p.Ser126Leu, refers to p.Ser71Leu in the mature protein) in a six-generation Chinese kindred characterized by dominant progressive nonsyndromic hearing loss, designated as DFNA64. SMAC/DIABLO is highly expressed in human embryonic ears and is enriched in the developing mouse inner-ear hair cells, suggesting it has a role in the development and homeostasis of hair cells. We used a functional study to demonstrate that the SMAC/DIABLO(S71L) mutant, while retaining the proapoptotic function, triggers significant degradation of both wild-type and mutant SMAC/DIABLO and renders host mitochondria susceptible to calcium-induced loss of the membrane potential. Our work identifies DFNA64 as the human genetic disorder associated with SMAC/DIABLO malfunction and suggests that mutant SMAC/DIABLO(S71L) might cause mitochondrial dysfunction.

Impact of Rare and Multiple Concurrent Gene Fusions on Diagnostic DNA Methylation Classifier in Brain Tumors.

DNA methylation is an essential molecular assay for central nervous system (CNS) tumor diagnostics. While some fusions define specific brain tumors, others occur across many different diagnoses. We performed a retrospective analysis of 219 primary CNS tumors with whole genome DNA methylation and RNA next-generation sequencing. DNA methylation profiling results were compared with RNAseq detected gene fusions. We detected 105 rare fusions involving 31 driver genes, including 23 fusions previously not implicated in brain tumors. In addition, we identified 6 multi-fusion tumors. Rare fusions and multi-fusion events can impact the diagnostic accuracy of DNA methylation by decreasing confidence in the result, such as BRAF, RAF, or FGFR1 fusions, or result in a complete mismatch, such as NTRK, EWSR1, FGFR, and ALK fusions. IMPLICATIONS: DNA methylation signatures need to be interpreted in the context of pathology and discordant results warrant testing for novel and rare gene fusions.

Loss of DOK2 induces carboplatin resistance in ovarian cancer via suppression of apoptosis.

OBJECTIVE: Ovarian cancers are highly heterogeneous and while chemotherapy is the preferred treatment many patients are intrinsically resistant or quickly develop resistance. Furthermore, all tumors that recur ultimately become resistant. Recent evidence suggests that epigenetic deregulation may be a key factor in the onset and maintenance of chemoresistance. We set out to identify epigenetically silenced genes that affect chemoresistance. METHODS: The epigenomes of a total of 45 ovarian samples were analyzed to identify epigenetically altered genes that segregate with platinum response, and further filtered with expression data to identify genes that were suppressed. A tissue culture carboplatin resistance screen was utilized to functionally validate this set of candidate platinum resistance genes. RESULTS: Our screen correctly identified 19 genes that when suppressed altered the chemoresistance of the cells in culture. Of the genes identified in the screen we further characterized one gene, docking protein 2 (DOK2), an adapter protein downstream of tyrosine kinase, to determine if we could elucidate the mechanism by which it increased resistance. The loss of DOK2 decreased the level of apoptosis in response to carboplatin. Furthermore, in cells with reduced DOK2, the level of anoikis was decreased. CONCLUSIONS: We have developed a screening methodology that analyzes the epigenome and informatically identifies candidate genes followed by in vitro culture screening of the candidate genes. To validate our screening methodology we further characterized one candidate gene, DOK2, and showed that loss of DOK2 induces chemotherapy resistance by decreasing the level of apoptosis in response to treatment.

Multiregional genetic evolution of metastatic uveal melanoma.

Uveal melanoma (UM) is the most common primary intraocular malignancy in adults and leads to deadly metastases for which there is no approved treatment. Genetic events driving early tumor development are well-described, but those occurring later during metastatic progression remain poorly understood. We performed multiregional genomic sequencing on 22 tumors collected from two patients with widely metastatic UM who underwent rapid autopsy. We observed multiple seeding events from the primary tumors, metastasis-to-metastasis seeding, polyclonal seeding, and late driver variants in ATM, KRAS, and other genes previously unreported in UM. These findings reveal previously unrecognized temporal and anatomic complexity in the genetic evolution of metastatic uveal melanoma, and they highlight the distinction between early and late phases of UM genetic evolution with implications for novel therapeutic approaches.

RNA Sequencing of Primary Cutaneous and Breast-Implant Associated Anaplastic Large Cell Lymphomas Reveals Infrequent Fusion Transcripts and Upregulation of PI3K/AKT Signaling via Neurotrophin Pathway Genes.

Cutaneous and breast implant-associated anaplastic large-cell lymphomas (cALCLs and BI-ALCLs) are two localized forms of peripheral T-cell lymphomas (PTCLs) that are recognized as distinct entities within the family of ALCL. JAK-STAT signaling is a common feature of all ALCL subtypes, whereas DUSP22/IRF4, TP63 and TYK gene rearrangements have been reported in a proportion of ALK-negative sALCLs and cALCLs. Both cALCLs and BI-ALCLs differ in their gene expression profiles compared to PTCLs; however, a direct comparison of the genomic alterations and transcriptomes of these two entities is lacking. By performing RNA sequencing of 1385 genes (TruSight RNA Pan-Cancer, Illumina) in 12 cALCLs, 10 BI-ALCLs and two anaplastic lymphoma kinase (ALK)-positive sALCLs, we identified the previously reported TYK2-NPM1 fusion in 1 cALCL (1/12, 8%), and four new intrachromosomal gene fusions in 2 BI-ALCLs (2/10, 20%) involving genes on chromosome 1 (EPS15-GNG12 and ARNT-GOLPH3L) and on chromosome 17 (MYO18A-GIT1 and NF1-GOSR1). One of the two BI-ALCL samples showed a complex karyotype, raising the possibility that genomic instability may be responsible for intra-chromosomal fusions in BI-ALCL. Moreover, transcriptional analysis revealed similar upregulation of the PI3K/Akt pathway, associated with enrichment in the expression of neurotrophin signaling genes, which was more conspicuous in BI-ALCL, as well as differences, i.e., over-expression of genes involved in the RNA polymerase II transcription program in BI-ALCL and of the RNA splicing/processing program in cALCL.

Metastatic pediatric sclerosing epithelioid fibrosarcoma.

Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft-tissue sarcoma thought to originate in fibroblasts of the tissues comprising tendons, ligaments, and muscles. Minimally responsive to conventional cytotoxic chemotherapies, >50% of SEF patients experience local recurrence and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80%-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-yr-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. Finally, we describe the generation of an SEF primary cell line, the first such culture to be reported. The patient described here experienced persistent disease progression despite aggressive treatment including multiple resections, radiotherapy, and numerous chemotherapies and targeted therapeutics. Untreated and locally recurrent tumor and metastatic tissue were sequenced by whole-genome, whole-exome, and deep-transcriptome next-generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and discuss potential therapeutic strategies for SEF.

Cell cycle kinases predicted from conserved biophysical properties.

Machine-learning techniques can classify functionally related proteins where homology-transfer as well as sequence and structure motifs fail. Here, we present a method that aimed at complementing homology-transfer in the identification of cell cycle control kinases from sequence alone. First, we identified functionally significant residues in cell cycle proteins through their high sequence conservation and biophysical properties. We then incorporated these residues and their features into support vector machines (SVM) to identify new kinases and more specifically to differentiate cell cycle kinases from other kinases and other proteins. As expected, the most informative residues tend to be highly conserved and tend to localize in the ATP binding regions of the kinases. Another observation confirmed that ATP binding regions are typically not found on the surface but in partially buried sites, and that this fact is correctly captured by accessibility predictions. Using these highly conserved, semi-buried residues and their biophysical properties, we could distinguish cell cycle S/T kinases from other kinase families at levels around 70-80% accuracy and 62-81% coverage. An application to the entire human proteome predicted at least 97 human proteins with limited previous annotations to be candidates for cell cycle kinases.

Physical protein-protein interactions predicted from microarrays.

MOTIVATION: Microarray expression data reveal functionally associated proteins. However, most proteins that are associated are not actually in direct physical contact. Predicting physical interactions directly from microarrays is both a challenging and important task that we addressed by developing a novel machine learning method optimized for this task. RESULTS: We validated our support vector machine-based method on several independent datasets. At the same levels of accuracy, our method recovered more experimentally observed physical interactions than a conventional correlation-based approach. Pairs predicted by our method to very likely interact were close in the overall network of interaction, suggesting our method as an aid for functional annotation. We applied the method to predict interactions in yeast (Saccharomyces cerevisiae). A Gene Ontology function annotation analysis and literature search revealed several probable and novel predictions worthy of future experimental validation. We therefore hope our new method will improve the annotation of interactions as one component of multi-source integrated systems. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of targetable BRAF ΔN486_P490 variant by whole-genome sequencing leading to dabrafenib-induced remission of a BRAF-mutant pancreatic adenocarcinoma.

The tumor genome of a patient with advanced pancreatic cancer was sequenced to identify potential therapeutic targetable mutations after standard of care failed to produce any significant overall response. Matched tumor-normal whole-genome sequencing revealed somatic mutations in BRAF, TP53, CDKN2A, and a focal deletion of SMAD4 The BRAF variant was an in-frame deletion mutation (ΔN486_P490), which had been previously demonstrated to be a kinase-activating alteration in the BRAF kinase domain. Working with the Novartis patient assistance program allowed us to treat the patient with the BRAF inhibitor, dabrafenib. The patient's overall clinical condition improved dramatically with dabrafenib. Levels of serum tumor marker dropped immediately after treatment, and a subsequent CT scan revealed a significant decrease in the size of both primary and metastatic lesions. The dabrafenib-induced remission lasted for 6 mo. Preclinical studies published concurrently with the patient's treatment showed that the BRAF in-frame mutation (ΔNVTAP) induces oncogenic activation by a mechanism distinct from that induced by V600E, and that this difference dictates the responsiveness to different BRAF inhibitors. This study describes a dramatic instance of how high-level genomic technology and analysis was necessary and sufficient to identify a clinically logical treatment option that was then utilized and shown to be of clinical value for this individual.

Correction to: Sequencing and curation strategies for identifying candidate glioblastoma treatments.

Following publication of the original article [1], it was reported that the given name of the fourteenth author was incorrectly published. The incorrect and the correct names are given below.

Sequencing and curation strategies for identifying candidate glioblastoma treatments.

BACKGROUND: Prompted by the revolution in high-throughput sequencing and its potential impact for treating cancer patients, we initiated a clinical research study to compare the ability of different sequencing assays and analysis methods to analyze glioblastoma tumors and generate real-time potential treatment options for physicians. METHODS: A consortium of seven institutions in New York City enrolled 30 patients with glioblastoma and performed tumor whole genome sequencing (WGS) and RNA sequencing (RNA-seq; collectively WGS/RNA-seq); 20 of these patients were also analyzed with independent targeted panel sequencing. We also compared results of expert manual annotations with those from an automated annotation system, Watson Genomic Analysis (WGA), to assess the reliability and time required to identify potentially relevant pharmacologic interventions. RESULTS: WGS/RNAseq identified more potentially actionable clinical results than targeted panels in 90% of cases, with an average of 16-fold more unique potentially actionable variants identified per individual; 84 clinically actionable calls were made using WGS/RNA-seq that were not identified by panels. Expert annotation and WGA had good agreement on identifying variants [mean sensitivity = 0.71, SD = 0.18 and positive predictive value (PPV) = 0.80, SD = 0.20] and drug targets when the same variants were called (mean sensitivity = 0.74, SD = 0.34 and PPV = 0.79, SD = 0.23) across patients. Clinicians used the information to modify their treatment plan 10% of the time. CONCLUSION: These results present the first comprehensive comparison of technical and machine augmented analysis of targeted panel and WGS/RNA-seq to identify potential cancer treatments.

Whole Genome Sequencing-Based Discovery of Structural Variants in Glioblastoma.

Next-generation DNA sequencing (NGS) technologies are currently being applied in both research and clinical settings for the understanding and management of disease. The goal is to use high-throughput sequencing to identify specific variants that drive tumorigenesis within each individual's tumor genomic profile. The significance of copy number and structural variants in glioblastoma makes it essential to broaden the search beyond oncogenic single nucleotide variants toward whole genome profiles of genetic aberrations that may contribute to disease progression. The heterogeneity of glioblastoma and its variability of cancer driver mutations necessitate a more robust examination of a patient's tumor genome. Here, we present patient whole genome sequencing (WGS) information to identify oncogenic structural variants that may contribute to glioblastoma pathogenesis. We provide WGS protocols and bioinformatics approaches to identify copy number and structural variations in 41 glioblastoma patient samples. We present how WGS can identify structural diversity within glioblastoma samples. We specifically show how to apply current bioinformatics tools to detect EGFR variants and other structural aberrations from DNA whole genome sequencing and how to validate those variants within the laboratory. These comprehensive WGS protocols can provide additional information directing more precise therapeutic options in the treatment of glioblastoma.

A recurrent novel MGA-NUTM1 fusion identifies a new subtype of high-grade spindle cell sarcoma.

NUTM1-rearranged tumors are defined by the presence of a gene fusion between NUTM1 and various gene partners and typically follow a clinically aggressive disease course with poor outcomes despite conventional multimodality therapy. NUTM1-rearranged tumors display histologic features of a poorly differentiated carcinoma with areas of focal squamous differentiation and typically express the BRD4-NUTM1 fusion gene defining a distinct clinicopathologic entity-NUT carcinoma (NC). NCs with mesenchymal differentiation have rarely been described in the literature. In this report, we describe the characterization of two cases of high-grade spindle cell sarcoma harboring a novel MGA-NUTM1 fusion. Whole-genome sequencing identified the presence of complex rearrangements resulting in a MGA-NUTM1 fusion gene in the absence of other significant somatic mutations. Genetic rearrangement was confirmed by fluorescence in situ hybridization, and expression of the fusion gene product was confirmed by transcriptomic analysis. The fusion protein was predicted to retain nearly the entire protein sequence of both MGA (exons 1-22) and NUTM1 (exons 3-8). Histopathologically, both cases were high-grade spindle cell sarcomas without specific differentiation markers. In contrast to typical cases of NC, these cases were successfully treated with aggressive local control measures (surgery and radiation) and both patients remain alive without disease. These cases describe a new subtype of NUTM1-rearranged tumors warranting expansion of diagnostic testing to evaluate for the presence of MGA-NUTM1 or alternative NUTM1 gene fusions in the diagnostic workup of high-grade spindle cell sarcomas or small round blue cell tumors of ambiguous lineage.

Analytical Validation of Clinical Whole-Genome and Transcriptome Sequencing of Patient-Derived Tumors for Reporting Targetable Variants in Cancer.

We developed and validated a clinical whole-genome and transcriptome sequencing (WGTS) assay that provides a comprehensive genomic profile of a patient's tumor. The ability to fully capture the mappable genome with sufficient sequencing coverage to precisely call DNA somatic single nucleotide variants, insertions/deletions, copy number variants, structural variants, and RNA gene fusions was analyzed. New York State's Department of Health next-generation DNA sequencing guidelines were expanded for establishing performance validation applicable to whole-genome and transcriptome sequencing. Whole-genome sequencing laboratory protocols were validated for the Illumina HiSeq X Ten platform and RNA sequencing for Illumina HiSeq2500 platform for fresh or frozen and formalin-fixed, paraffin-embedded tumor samples. Various bioinformatics tools were also tested, and CIs for sensitivity and specificity thresholds in calling clinically significant somatic aberrations were determined. The validation was performed on a set of 125 tumor normal pairs. RNA sequencing was performed to call fusions and to confirm the DNA variants or exonic alterations. Here, we present our results and WGTS standards for variant allele frequency, reproducibility, analytical sensitivity, and present limit of detection analysis for single nucleotide variant calling, copy number identification, and structural variants. We show that The New York Genome Center WGTS clinical assay can provide a comprehensive patient variant discovery approach suitable for directed oncologic therapeutic applications.