Liu Xi

A comprehensive assessment of somatic mutation detection in cancer using whole-genome sequencing

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Genomic profiling of Sézary syndrome identifies alterations of key T cell signaling and differentiation genes

Sézary syndrome is a rare leukemic form of cutaneous T cell lymphoma characterized by generalized redness, scaling, itching and increased numbers of circulating atypical T lymphocytes. It is rarely curable, with poor prognosis. Here we present a multiplatform genomic analysis of 37 patients with Sézary syndrome that implicates dysregulation of cell cycle checkpoint and T cell signaling. Frequent somatic alterations were identified in TP53, CARD11, CCR4, PLCG1, CDKN2A, ARID1A, RPS6KA1 and ZEB1. Activating CCR4 and CARD11 mutations were detected in nearly one-third of patients. ZEB1, encoding a transcription repressor essential for T cell differentiation, was deleted in over one-half of patients. IL32 and IL2RG were overexpressed in nearly all cases. Our results demonstrate profound disruption of key signaling pathways in Sézary syndrome and suggest potential targets for new therapies.

Multilevel Genomics-Based Taxonomy of Renal Cell Carcinoma

On the basis of multidimensional and comprehensive molecular characterization (including DNA methalylation and copy number, RNA, and protein expression), we classified 894 renal cell carcinomas (RCCs) of various histologic types into nine major genomic subtypes. Site of origin within the nephron was one major determinant in the classification, reflecting differences among clear cell, chromophobe, and papillary RCC. Widespread molecular changes associated with TFE3 gene fusion or chromatin modifier genes were present within a specific subtype and spanned multiple subtypes. Differences in patient survival and in alteration of specific pathways (including hypoxia, metabolism, MAP kinase, NRF2-ARE, Hippo, immune checkpoint, and PI3K/AKT/mTOR) could further distinguish the subtypes. Immune checkpoint markers and molecular signatures of T cell infiltrates were both highest in the subtype associated with aggressive clear cell RCC. Differences between the genomic subtypes suggest that therapeutic strategies could be tailored to each RCC disease subset.

Ampullary Cancers Harbor ELF3 Tumor Suppressor Gene Mutations and Exhibit Frequent WNT Dysregulation

The ampulla of Vater is a complex cellular environment from which adenocarcinomas arise to form a group of histopathologically heterogenous tumors. To evaluate the molecular features of these tumors, 98 ampullary adenocarcinomas were evaluated and compared to 44 distal bile duct and 18 duodenal adenocarcinomas. Genomic analyses revealed mutations in the WNT signaling pathway among half of the patients and in all three adenocarcinomas irrespective of their origin and histological morphology. These tumors were characterized by a high frequency of inactivating mutations of ELF3, a high rate of microsatellite instability, and common focal deletions and amplifications, suggesting common attributes in the molecular pathogenesis are at play in these tumors. The high frequency of WNT pathway activating mutation, coupled with small-molecule inhibitors of β-catenin in clinical trials, suggests future treatment decisions for these patients may be guided by genomic analysis.

Effects of TP53 mutational status on gene expression patterns across 10 human cancer types

Mutations in the TP53 tumour suppressor gene occur in half of all human cancers, indicating its critical importance in inhibiting cancer development. Despite extensive studies, the mechanisms by which mutant p53 enhances tumour progression remain only partially understood. Here, using data from the Cancer Genome Atlas (TCGA), genomic and transcriptomic analyses were performed on 2256 tumours from 10 human cancer types. We show that tumours with TP53 mutations have altered gene expression profiles compared to tumours retaining two wild‐type TP53 alleles. Among 113 known p53‐up‐regulated target genes identified from cell culture assays, 10 were consistently up‐regulated in at least eight of 10 cancer types that retain both copies of wild‐type TP53. RPS27L, CDKN1A (p21CIP1) and ZMAT3 were significantly up‐regulated in all 10 cancer types retaining wild‐type TP53. Using this p53‐based expression analysis as a discovery tool, we used cell‐based assays to identify five novel p53 target genes from genes consistently up‐regulated in wild‐type p53 cancers. Global gene expression analyses revealed that cell cycle regulatory genes and transcription factors E2F1, MYBL2 and FOXM1 were disproportionately up‐regulated in many TP53 mutant cancer types. Finally, > 93% of tumours with a TP53 mutation exhibited greatly reduced wild‐type p53 messenger expression, due to loss of heterozygosity or copy neutral loss of heterozygosity, supporting the concept of p53 as a recessive tumour suppressor. The data indicate that tumours with wild‐type TP53 retain some aspects of p53‐mediated growth inhibitory signalling through activation of p53 target genes and suppression of cell cycle regulatory genes. Published by John Wiley & Sons, Ltd. www.pathsoc.org.uk

Comprehensive Genomic Characterization of Upper Tract Urothelial Carcinoma

BACKGROUNDnUpper urinary tract urothelial cancer (UTUC) may have unique etiologic and genomic factors compared to bladder cancer.nnnOBJECTIVEnTo characterize the genomic landscape of UTUC and provide insights into its biology using comprehensive integrated genomic analyses.nnnDESIGN, SETTING, AND PARTICIPANTSnWe collected 31 untreated snap-frozen UTUC samples from two institutions and carried out whole-exome sequencing (WES) of DNA, RNA sequencing (RNAseq), and protein analysis.nnnOUTCOME MEASUREMENTS AND STATISTICAL ANALYSISnAdjusting for batch effects, consensus mutation calls from independent pipelines identified DNA mutations, gene expression clusters using unsupervised consensus hierarchical clustering (UCHC), and protein expression levels that were correlated with relevant clinical variables, The Cancer Genome Atlas, and other published data.nnnRESULTS AND LIMITATIONSnWES identified mutations in FGFR3 (74.1%; 92% low-grade, 60% high-grade), KMT2D (44.4%), PIK3CA (25.9%), and TP53 (22.2%). APOBEC and CpG were the most common mutational signatures. UCHC of RNAseq data segregated samples into four molecular subtypes with the following characteristics. Cluster 1: no PIK3CA mutations, nonsmokers, high-grade <pT2 tumors, high recurrences. Cluster 2: 100% FGFR3 mutations, low-grade tumors, tobacco use, noninvasive disease, no bladder recurrences. Cluster 3: 100% FGFR3 mutations, 71% PIK3CA, no TP53 mutations, five bladder recurrences, tobacco use, tumors all <pT2. Cluster 4: KMT2D (62.5%), FGFR3 (50%), TP53 (50%) mutations, no PIK3CA mutations, high-grade pT2+ disease, tobacco use, carcinoma in situ, shorter survival. We identified a novel SH3KBP1-CNTNAP5 fusion.nnnCONCLUSIONSnMutations in UTUC occur at differing frequencies from bladder cancer, with four unique molecular and clinical subtypes. A novel SH3KBP1 fusion regulates RTK signaling. Further studies are needed to validate the described subtypes, explore their responses to therapy, and better define the novel fusion mutation.nnnPATIENT SUMMARYnWe conducted a comprehensive study of the genetics of upper urinary tract urothelial cancer by evaluating DNA, RNA and protein expression in 31 tumors. We identified four molecular subtypes with distinct behaviors. Future studies will determine if these subtypes appear to have different responses to treatments.

MuSE: accounting for tumor heterogeneity using a sample-specific error model improves sensitivity and specificity in mutation calling from sequencing data

Subclonal mutations reveal important features of the genetic architecture of tumors. However, accurate detection of mutations in genetically heterogeneous tumor cell populations using next-generation sequencing remains challenging. We develop MuSE (http://bioinformatics.mdanderson.org/main/MuSE), Mutation calling using a Markov Substitution model for Evolution, a novel approach for modeling the evolution of the allelic composition of the tumor and normal tissue at each reference base. MuSE adopts a sample-specific error model that reflects the underlying tumor heterogeneity to greatly improve the overall accuracy. We demonstrate the accuracy of MuSE in calling subclonal mutations in the context of large-scale tumor sequencing projects using whole exome and whole genome sequencing.

A comprehensive multicenter comparison of whole genome sequencing pipelines using a uniform tumor-normal sample pair

As next-generation sequencing becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Through the International Cancer Genome Consortium (ICGC), we compared sequencing pipelines at five independent centers (CNAG, DKFZ, OICR, RIKEN and WTSI) using a single tumor-blood DNA pair. Analyses by each center and with one standardized algorithm revealed significant discrepancies. Although most pipelines performed well for coding mutations, library preparation methods and sequencing coverage metrics clearly influenced downstream results. PCR-free methods showed reduced GC-bias and more even coverage. Increasing sequencing depth to ∼100x (two- to three-fold higher than current standards) showed a benefit, as long as the tumor:control coverage ratio remained balanced. To become part of routine clinical care, high-throughput sequencing must be globally compatible and comparable. This benchmarking exercise has highlighted several fundamental parameters to consider in this regard, which will allow for better optimization and planning of both basic and translational studies.

A Comprehensive Assessment of Somatic Mutation Calling in Cancer Genomes

The emergence of next generation DNA sequencing technology is enabling high-resolution cancer genome analysis. Large-scale projects like the International Cancer Genome Consortium (ICGC) are systematically scanning cancer genomes to identify recurrent somatic mutations. Second generation DNA sequencing, however, is still an evolving technology and procedures, both experimental and analytical, are constantly changing. Thus the research community is still defining a set of best practices for cancer genome data analysis, with no single protocol emerging to fulfil this role. Here we describe an extensive benchmark exercise to identify and resolve issues of somatic mutation calling. Whole genome sequence datasets comprising tumor-normal pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, were shared within the ICGC and submissions of somatic mutation calls were compared to verified mutations and to each other. Varying strategies to call mutations, incomplete awareness of sources of artefacts, and even lack of agreement on what constitutes an artefact or real mutation manifested in widely varying mutation call rates and somewhat low concordance among submissions. We conclude that somatic mutation calling remains an unsolved problem. However, we have identified many issues that are easy to remedy that are presented here. Our study highlights critical issues that need to be addressed before this valuable technology can be routinely used to inform clinical decision-making. Abbreviations and Definitions SSM Somatic Single-base Mutations or Simple Somatic Mutations, refers to a somatic single base change SIM Somatic Insertion/deletion Mutation CNV Copy Number Variant SV Structural Variant SNP Single Nucleotide Polymorphisms, refers to a single base variable position in the germline with a frequency of > 1% in the general population CLL Chronic Lymphocytic Leukaemia MB Medulloblastoma ICGC International Cancer Genome Consortium BM Benchmark aligner = mapper, these terms are used interchangeably

ITD assembler: an algorithm for internal tandem duplication discovery from short-read sequencing data

BackgroundDetection of tandem duplication within coding exons, referred to as internal tandem duplication (ITD), remains challenging due to inefficiencies in alignment of ITD-containing reads to the reference genome. There is a critical need to develop efficient methods to recover these important mutational events.ResultsIn this paper we introduce ITD Assembler, a novel approach that rapidly evaluates all unmapped and partially mapped reads from whole exome NGS data using a De Bruijn graphs approach to select reads that harbor cycles of appropriate length, followed by assembly using overlap-layout-consensus. We tested ITD Assembler on The Cancer Genome Atlas AML dataset as a truth set. ITD Assembler identified the highest percentage of reported FLT3-ITDs when compared to other ITD detection algorithms, and discovered additional ITDs in FLT3, KIT, CEBPA, WT1 and other genes. Evidence of polymorphic ITDs in 54 genes were also found. Novel ITDs were validated by analyzing the corresponding RNA sequencing data.ConclusionsITD Assembler is a very sensitive tool which can detect partial, large and complex tandem duplications. This study highlights the need to more effectively look forxa0ITD’s in other cancers and Mendelian diseases.