Comprehensive Characterization of Cancer Driver Genes and Mutations.

Identifying molecular cancer drivers is critical for precision oncology. Multiple advanced algorithms to identify drivers now exist, but systematic attempts to combine and optimize them on large datasets are few. We report a PanCancer and PanSoftware analysis spanning 9,423 tumor exomes (comprising all 33 of The Cancer Genome Atlas projects) and using 26 computational tools to catalog driver genes and mutations. We identify 299 driver genes with implications regarding their anatomical sites and cancer/cell types. Sequence- and structure-based analyses identified >3,400 putative missense driver mutations supported by multiple lines of evidence. Experimental validation confirmed 60%-85% of predicted mutations as likely drivers. We found that >300 MSI tumors are associated with high PD-1/PD-L1, and 57% of tumors analyzed harbor putative clinically actionable events. Our study represents the most comprehensive discovery of cancer genes and mutations to date and will serve as a blueprint for future biological and clinical endeavors.

Simultaneous Copy Number Alteration and Single-Nucleotide Variation Analysis in Matched Aqueous Humor and Tumor Samples in Children with Retinoblastoma.

Retinoblastoma (RB) is a childhood cancer that forms in the developing retina of young children; this tumor cannot be biopsied due to the risk of provoking extraocular tumor spread, which dramatically alters the treatment and survival of the patient. Recently, aqueous humor (AH), the clear fluid in the anterior chamber of the eye, has been developed as an organ-specific liquid biopsy for investigation of in vivo tumor-derived information found in the cell-free DNA (cfDNA) of the biofluid. However, identifying somatic genomic alterations, including both somatic copy number alterations (SCNAs) and single nucleotide variations (SNVs) of the RB1 gene, typically requires either: (1) two distinct experimental protocols-low-pass whole genome sequencing for SCNAs and targeted sequencing for SNVs-or (2) expensive deep whole genome or exome sequencing. To save time and cost, we applied a one-step targeted sequencing method to identify both SCNAs and RB1 SNVs in children with RB. High concordance (median = 96.2%) was observed in comparing SCNA calls derived from targeted sequencing to the traditional low-pass whole genome sequencing method. We further applied this method to investigate the degree of concordance of genomic alterations between paired tumor and AH samples from 11 RB eyes. We found 11/11 AH samples (100%) had SCNAs, and 10 of them (90.1%) with recurrent RB-SCNAs, while only nine out of 11 tumor samples (81.8%) had positive RB-SCNA signatures in both low-pass and targeted methods. Eight out of the nine (88.9%) detected SNVs were shared between AH and tumor samples. Ultimately, 11/11 cases have somatic alterations identified, including nine RB1 SNVs and 10 recurrent RB-SCNAs with four focal RB1 deletions and one MYCN gain. The results presented show the feasibility of utilizing one sequencing approach to obtain SCNA and targeted SNV data to capture a broad genomic scope of RB disease, which may ultimately expedite clinical intervention and be less expensive than other methods.

Proteogenomics connects somatic mutations to signalling in breast cancer.

Somatic mutations have been extensively characterized in breast cancer, but the effects of these genetic alterations on the proteomic landscape remain poorly understood. Here we describe quantitative mass-spectrometry-based proteomic and phosphoproteomic analyses of 105 genomically annotated breast cancers, of which 77 provided high-quality data. Integrated analyses provided insights into the somatic cancer genome including the consequences of chromosomal loss, such as the 5q deletion characteristic of basal-like breast cancer. Interrogation of the 5q trans-effects against the Library of Integrated Network-based Cellular Signatures, connected loss of CETN3 and SKP1 to elevated expression of epidermal growth factor receptor (EGFR), and SKP1 loss also to increased SRC tyrosine kinase. Global proteomic data confirmed a stromal-enriched group of proteins in addition to basal and luminal clusters, and pathway analysis of the phosphoproteome identified a G-protein-coupled receptor cluster that was not readily identified at the mRNA level. In addition to ERBB2, other amplicon-associated highly phosphorylated kinases were identified, including CDK12, PAK1, PTK2, RIPK2 and TLK2. We demonstrate that proteogenomic analysis of breast cancer elucidates the functional consequences of somatic mutations, narrows candidate nominations for driver genes within large deletions and amplified regions, and identifies therapeutic targets.

Potential of Aqueous Humor as a Liquid Biopsy for Uveal Melanoma.

Tumor biopsy can identify prognostic biomarkers for metastatic uveal melanoma (UM), however aqueous humor (AH) liquid biopsy may serve as an adjunct. This study investigated whether the AH of UM eyes has sufficient circulating tumor DNA (ctDNA) to perform genetic analysis. This is a case series of 37 AH samples, taken before or after radiation, and one tumor wash sample, from 12 choroidal and 8 ciliary body (CB) melanoma eyes. AH was analyzed for nucleic acid concentrations. AH DNA and one tumor wash sample underwent shallow whole-genome sequencing followed by Illumina sequencing to detect somatic copy number alterations (SCNAs). Four post-radiation AH underwent targeted sequencing of BAP1 and GNAQ genes. Post-radiation AH had significantly higher DNA and miRNA concentrations than paired pre-radiation samples. Highly recurrent UM SCNAs were identified in 0/11 post-radiation choroidal and 6/8 post-radiation CB AH. SCNAs were highly concordant in a CB post-radiation AH with its matched tumor (r = 0.978). BAP1 or GNAQ variants were detected in 3/4 post-radiation AH samples. AH is a source of ctDNA in UM eyes, particularly in post-radiation CB eyes. For the first time, UM SCNAs and mutations were identified in AH-derived ctDNA. Suggesting that AH can serve as a liquid biopsy for UM.

GenomeVIP: a cloud platform for genomic variant discovery and interpretation.

Identifying genomic variants is a fundamental first step toward the understanding of the role of inherited and acquired variation in disease. The accelerating growth in the corpus of sequencing data that underpins such analysis is making the data-download bottleneck more evident, placing substantial burdens on the research community to keep pace. As a result, the search for alternative approaches to the traditional "download and analyze" paradigm on local computing resources has led to a rapidly growing demand for cloud-computing solutions for genomics analysis. Here, we introduce the Genome Variant Investigation Platform (GenomeVIP), an open-source framework for performing genomics variant discovery and annotation using cloud- or local high-performance computing infrastructure. GenomeVIP orchestrates the analysis of whole-genome and exome sequence data using a set of robust and popular task-specific tools, including VarScan, GATK, Pindel, BreakDancer, Strelka, and Genome STRiP, through a web interface. GenomeVIP has been used for genomic analysis in large-data projects such as the TCGA PanCanAtlas and in other projects, such as the ICGC Pilots, CPTAC, ICGC-TCGA DREAM Challenges, and the 1000 Genomes SV Project. Here, we demonstrate GenomeVIP's ability to provide high-confidence annotated somatic, germline, and de novo variants of potential biological significance using publicly available data sets.

Stability and uniqueness of clonal immunoglobulin CDR3 sequences for MRD tracking in multiple myeloma.

Minimal residual disease (MRD) tracking, by next generation sequencing of immunoglobulin sequences, is moving towards clinical implementation in multiple myeloma. However, there is only sparse information available to address whether clonal sequences remain stable for tracking over time, and to what extent light chain sequences are sufficiently unique for tracking. Here, we analyzed immunoglobulin repertoires from 905 plasma cell myeloma and healthy control samples, focusing on the third complementarity determining region (CDR3). Clonal heavy and/or light chain expression was identified in all patients at baseline, with one or more subclones related to the main clone in 3.2%. In 45 patients with 101 sequential samples, the dominant clonal CDR3 sequences remained identical over time, despite differential clonal evolution by whole exome sequencing in 49% of patients. The low frequency of subclonal CDR3 variants, and absence of evolution over time in active multiple myeloma, indicates that tumor cells at this stage are not under selective pressure to undergo antibody affinity maturation. Next, we establish somatic hypermutation and non-templated insertions as the most important determinants of light chain clonal uniqueness, identifying a potentially trackable sequence in the majority of patients. Taken together, we show that dominant clonal sequences identified at baseline are reliable biomarkers for long-term tracking of the malignant clone, including both IGH and the majority of light chain clones.

Tumor suppressor CD99 is downregulated in plasma cell neoplasms lacking CCND1 translocation and distinguishes neoplastic from normal plasma cells and B-cell lymphomas with plasmacytic differentiation from primary plasma cell neoplasms.

CD99(MIC2) is a widely expressed cell surface glycoprotein and functions as a tumor suppressor involved in downregulation of SRC family of tyrosine kinase. CD99 expression is tightly regulated through B-cell development. The principal aims of this study were to investigate the clinical utility of CD99 expression (i) in distinguishing normal plasma cells from primary plasma cell neoplasms; (ii) in detection of minimal residual disease in primary plasma cell neoplasms; and (iii) in distinguishing plasma cell component of B-cell lymphomas from primary plasma cell neoplasms. We analyzed expression of CD99 by flow cytometry and immunohistochemistry in lymph nodes, peripheral blood, and bone marrow samples. CD99 showed stage-specific expression with highest expression seen in precursor B and plasma cells. In contrast to the uniform bright expression on normal plasma cells, CD99 expression on neoplastic plasma cells was lost in 39 out of 56 (69.6%) cases. Furthermore, 8 out of 56 samples (14%) showed visibly (>10-fold) reduced CD99 expression. Overall, CD99 expression was informative (absent or visibly dimmer than normal) in 84% of primary plasma cell neoplasm. In the context of minimal residual disease detection, CD99 showed superior utility in separating normal and abnormal plasma cells over currently established antigens CD117, CD81, and CD27 by principal component analysis. Preservation of CD99 expression was strongly associated with cyclin D1 translocation in myeloma (p < 0.05). B-cell lymphomas with plasma cell component could be distinguished from myeloma by CD99 expression. In summary, we established that tumor suppressor CD99 is markedly downregulated in multiple myeloma. The loss is highly specific for identification of abnormal cells in primary plasma cell neoplasms, and can be exploited for diagnostic purposes. The role of CD99 in myeloma pathogenesis requires further investigation.

Transcriptome Sequencing Allows Comprehensive Genomic Characterization of Pediatric B-Acute Lymphoblastic Leukemia in an Academic Clinical Laboratory.

Studies have shown the power of transcriptome sequencing [RNA sequencing (RNA-Seq)] in identifying known and novel oncogenic drivers and molecular subtypes of B-acute lymphoblastic leukemia (B-ALL). The current study investigated whether the clinically validated RNA-Seq assay, coupled with a custom analysis pipeline, could be used for a comprehensive B-ALL classification. Following comprehensive clinical testing, RNA-Seq was performed on 76 retrospective B-ALL cases, 28 of which had known and 48 had undetermined subtype. Subtypes were accurately identified in all 28 known cases, and in 38 of 48 unknown cases (79%). The subtypes of the unknown cases included the following: PAX5alt (n = 12), DUX4-rearranged (n = 6), Philadelphia chromosome-like (n = 5), low hyperdiploid (n = 4), ETV6::RUNX1-like (n = 3), MEF2D-rearranged (n = 2), PAX5 P80R (n = 2), ZEB2/CEBP (n = 1), NUTM1-rearranged (n = 1), ZNF384-rearranged (n = 1), and TCF3::PBX1 (n = 1). In 15 of 38 cases (39%), classification based on expression profile was corroborated by detection of subtype-defining oncogenic drivers missed by clinical testing. RNA-Seq analysis also detected large copy number abnormalities, oncogenic hot-spot sequence variants, and intragenic IKZF1 deletions. This pilot study confirms the feasibility of implementing an RNA-Seq workflow for clinical diagnosis of molecular subtypes in pediatric B-ALL, reinforcing that RNA-Seq represents a promising global genomic assay for this heterogeneous leukemia.

Pediatric Chordoma: A Tale of Two Genomes.

Little is known about the genomic alterations in chordoma, with the exception of loss of SMARCB1, a core member of the SWI/SNF complex, in poorly differentiated chordomas. A TBXT duplication and rs2305089 polymorphism, located at 6q27, are known genetic susceptibility loci. A comprehensive genomic analysis of the nuclear and mitochondrial genomes in pediatric chordoma has not yet been reported. In this study, we performed WES and mtDNA genome sequencing on 29 chordomas from 23 pediatric patients. Findings were compared with that from whole-genome sequencing datasets of 80 adult patients with skull base chordoma. In the pediatric chordoma cohort, 81% of the somatic mtDNA mutations were observed in NADH complex genes, which is significantly enriched compared with the rest of the mtDNA genes (P = 0.001). In adult chordomas, mtDNA mutations were also enriched in the NADH complex genes (P < 0.0001). Furthermore, a progressive increase in heteroplasmy of nonsynonymous mtDNA mutations was noted in patients with multiple tumors (P = 0.0007). In the nuclear genome, rare likely germline in-frame indels in ARID1B, a member of the SWI/SNF complex located at 6q25.3, were observed in five pediatric patients (22%) and four patients in the adult cohort (5%). The frequency of rare ARID1B indels in the pediatric cohort is significantly higher than that in the adult cohort (P = 0.0236, Fisher's exact test), but they were both significantly higher than that in the ethnicity-matched populations (P < 5.9e-07 and P < 0.0001174, respectively). Implications: germline ARID1B indels and mtDNA aberrations seem important for chordoma genesis, especially in pediatric chordoma.

Combined low-pass whole genome and targeted sequencing in liquid biopsies for pediatric solid tumors.

We designed a liquid biopsy (LB) platform employing low-pass whole genome sequencing (LP-WGS) and targeted sequencing of cell-free (cf) DNA from plasma to detect genome-wide copy number alterations (CNAs) and gene fusions in pediatric solid tumors. A total of 143 plasma samples were analyzed from 19 controls and 73 patients, including 44 bone or soft-tissue sarcomas and 12 renal, 10 germ cell, five hepatic, and two thyroid tumors. cfDNA was isolated from plasma collected at diagnosis, during and after therapy, and/or at relapse. Twenty-six of 37 (70%) patients enrolled at diagnosis without prior therapy (radiation, surgery, or chemotherapy) had circulating tumor DNA (ctDNA), based on the detection of CNAs from LP-WGS, including 18 of 27 (67%) patients with localized disease and eight of 10 (80%) patients with metastatic disease. None of the controls had detectable somatic CNAs. There was a high concordance of CNAs identified by LP-WGS to CNAs detected by chromosomal microarray analysis in the matching tumors. Mutations identified in tumor samples with our next-generation sequencing (NGS) panel, OncoKids®, were also detected by LP-WGS of ctDNA in 14 of 26 plasma samples. Finally, we developed a hybridization-based capture panel to target EWSR1 and FOXO1 fusions from patients with Ewing sarcoma or alveolar rhabdomyosarcoma (ARMS), respectively. Fusions were detected in the plasma from 10 of 12 patients with Ewing sarcoma and in two of two patients with ARMS. Combined, these data demonstrate the clinical applicability of our LB platform to evaluate pediatric patients with a variety of solid tumors.

Low-pass whole-genome and targeted sequencing of cell-free DNA from cerebrospinal fluid in pediatric patients with central nervous system tumors.

Background: Central nervous system tumors are the most common pediatric solid tumors and the most frequent cause of cancer-related morbidity in childhood. Significant advances in understanding the molecular features of these tumors have facilitated the development of liquid biopsy assays that may aid in diagnosis and monitoring response to therapy. In this report, we describe our comprehensive liquid biopsy platform for detection of genome-wide copy number aberrations, sequence variants, and gene fusions using cerebrospinal fluid (CSF) from pediatric patients with brain, spinal cord, and peripheral nervous system tumors. Methods: Cell-free DNA was isolated from the CSF from 55 patients, including 47 patients with tumors and 8 controls. Results: Abnormalities in cell-free DNA were detected in 24 (51%) patients including 11 with copy number alterations, 9 with sequence variants, and 7 with KIAA1549::BRAF fusions. Positive findings were obtained in patients spanning histologic subtypes, tumor grades, and anatomic locations. Conclusions: This study demonstrates the feasibility of employing this platform in routine clinical care in upfront diagnostic and monitoring settings. Future studies are required to determine the utility of this approach for assessing response to therapy and long-term surveillance.

Aqueous Humor Liquid Biopsy as a Companion Diagnostic for Retinoblastoma: Implications for Diagnosis, Prognosis, and Therapeutic Options: Five Years of Progress.

PURPOSE: To define the prospective use of the aqueous humor (AH) as a molecular diagnostic and prognostic liquid biopsy for retinoblastoma (RB). METHODS: This is a prospective, observational study wherein an AH liquid biopsy is performed at diagnosis and longitudinally through therapy for patients with RB. Tumor-derived cell-free DNA is isolated and sequenced for single nucleotide variant analysis of the RB1 gene and detection of somatic copy number alterations (SCNAs). The SCNAs are used to determine tumor fraction (TFx). Specific SCNAs, including 6p gain and focal MycN gain, along with TFx, are prospectively correlated with intraocular tumor relapse, response to therapy, and globe salvage. RESULTS: A total of 26 eyes of 21 patients were included with AH taken at diagnosis. Successful ocular salvage was achieved in 19 of 26 (73.1%) eyes. Mutational analysis of 26 AH samples identified 23 pathogenic RB1 variants and 2 focal RB1 deletions; variant allele fraction ranged from 30.5% to 100% (median 93.2%). At diagnosis, SCNAs were detectable in 17 of 26 (65.4%) AH samples. Eyes with 6p gain and/or focal MycN gain had significantly greater odds of poor therapeutic outcomes (odds ratio = 6.75, 95% CI = 1.06-42.84, P = .04). Higher AH TFx was observed in eyes with vitreal progression (TFx = 46.0% ± 40.4) than regression (22.0 ± 29.1; difference: -24.0; P = .049). CONCLUSIONS: Establishing an AH liquid biopsy for RB is aimed at addressing (1) our inability to biopsy tumor tissue and (2) the lack of molecular biomarkers for intraocular prognosis. Current management decisions for RB are made based solely on clinical features without objective molecular testing. This prognostic study shows great promise for using AH as a companion diagnostic. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Waldenström macroglobulinemia whole genome reveals prolonged germinal center activity and late copy number aberrations.

The genomic landscape of Waldenström macroglobulinemia (WM) is characterized by somatic mutations in MYD88, present from the precursor stages. Using the comprehensive resolution of whole genome sequencing (WGS) in 14 CD19-selected primary WM samples; comparing clonal and subclonal mutations revealed that germinal center (GC) mutational signatures SBS9 (poly-eta) and SBS84 (AID) have sustained activity, suggesting that the interaction between WM and the GC continues over time. Expanding our cohort size with 33 targeted sequencing samples, we interrogated the WM copy number aberration (CNA) landscape and chronology. Of interest, CNA prevalence progressively increased in symptomatic WM and relapsed disease when compared with stable precursor stages, with stable precursors lacking genomic complexity. Two MYD88 wild-type WGS contained a clonal gain affecting chromosome 12, which is typically an early event in chronic lymphocytic leukemia. Molecular time analysis demonstrated that both chromosomal 12 gain events occurred early in cancer development whereas other CNA changes tend to occur later in the disease course and are often subclonal. In summary, WGS analysis in WM allows the demonstration of sustained GC activity over time and allows the reconstruction of the temporal evolution of specific genomic features. In addition, our data suggest that, although MYD88-mutations are central to WM clone establishment and can be observed in precursor disease, CNA may contribute to later phases, and may be used as a biomarker for progression risk from precursor conditions to symptomatic disease.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells.

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis. SIGNIFICANCE: Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

Malignant teratoid intraocular ciliary body medulloepithelioma in a 5-year-old male with corresponding somatic copy number alteration profile of aqueous humor cell-free DNA.

BACKGROUND: Intraocular, ciliary body, medulloepithelioma (CBME) is a rare tumor of the nonpigmented ciliary body epithelium, typically presenting in childhood. We describe a case of CBME. MATERIALS AND METHODS: Ocular examination and imaging guided diagnostic and treatment decisions. Aqueous humor (AH) liquid biopsy was collected from the affected eye at eventual enucleation. Whole genome sequencing (WGS) was employed to determine somatic copy number alterations (SCNA) in AH cell-free DNA (cfDNA). Tumor sample was analyzed using various assays to evaluate for oncogenic mutations and SCNAs. Histopathology determined diagnosis. RESULTS: A 5-year-old male with glaucoma and cataract in the left eye (OS) experienced worsening left eye pain and redness. There was no light perception OS and the eye was hypotonus. Anterior segment exam showed complete cataract and rubeosis iridis. Ocular B-scan ultrasound OS revealed an intraocular lesion with calcifications and retinal detachment. Orbital MRI suggested left globe hypercellularity. An infiltrative lesion involving the ciliary body was seen in the left eye on examination under anesthesia. Left eye enucleation was performed in the setting of pain, blindness, and tumor, with anterior chamber paracentesis for AH liquid biopsy collection. SCNA profile of AH cfDNA demonstrated loss of copy of chromosomes 4, 6, and 9. Tumor was negative for clinically significant mutations or SCNAs. Histopathology diagnosed malignant teratoid CBME. CONCLUSIONS: We present a case of CBME and include the unique SCNA profile of AH cfDNA from the enucleated eye. This case suggests utility of AH liquid biopsy in distinguishing between differential diagnoses for intraocular mass lesions.

Capture Rate of V(D)J Sequencing for Minimal Residual Disease Detection in Multiple Myeloma.

PURPOSE: Minimal residual disease (MRD) negativity is a strong predictor for outcome in multiple myeloma. To assess V(D)J clonotype capture using the updated Adaptive next-generation sequencing (NGS) MRD assay in a clinical setting, we analyzed baseline and follow-up samples from patients with multiple myeloma who achieved deep clinical responses. EXPERIMENTAL DESIGN: A total of 159 baseline and 31 follow-up samples from patients with multiple myeloma were sequenced using the NGS MRD assay. Baseline samples were also sequenced using a targeted multiple myeloma panel (myTYPE). We estimated ORs with 95% confidence intervals (CI) for clonotypes detection using logistic regression. RESULTS: The V(D)J clonotype capture rate was 93% in baseline samples with detectable genomic aberrations, indicating presence of tumor DNA, assessed through myTYPE. myTYPE-positive samples had significantly higher V(D)J clonotype detection rates in univariate (OR, 7.3; 95% CI, 2.8-22.6) and multivariate analysis (OR, 4.4; 95% CI, 1.4-16.9; P = 0.016). Higher disease burden was associated with higher probability of V(D)J clonotype capture, meanwhile no such association was found for age, gender, or type of heavy or light immunoglobulin chain. All V(D)J clonotypes detected at baseline were detected in MRD-positive samples indicating that the V(D)J clonotypes remained stable and did not undergo further rearrangements during follow-up. Of the 31 posttreatment samples, 12 were MRD-negative using the NGS MRD assay. CONCLUSIONS: NGS for V(D)J rearrangements in multiple myeloma offers a reliable and sensitive method for MRD tracking with high detection rates in the clinical setting.

Copy number signatures predict chromothripsis and clinical outcomes in newly diagnosed multiple myeloma.

Chromothripsis is detectable in 20-30% of newly diagnosed multiple myeloma (NDMM) patients and is emerging as a new independent adverse prognostic factor. In this study we interrogate 752 NDMM patients using whole genome sequencing (WGS) to investigate the relationship of copy number (CN) signatures to chromothripsis and show they are highly associated. CN signatures are highly predictive of the presence of chromothripsis (AUC = 0.90) and can be used identify its adverse prognostic impact. The ability of CN signatures to predict the presence of chromothripsis is confirmed in a validation series of WGS comprised of 235 hematological cancers (AUC = 0.97) and an independent series of 34 NDMM (AUC = 0.87). We show that CN signatures can also be derived from whole exome data (WES) and using 677 cases from the same series of NDMM, we are able to predict both the presence of chromothripsis (AUC = 0.82) and its adverse prognostic impact. CN signatures constitute a flexible tool to identify the presence of chromothripsis and is applicable to WES and WGS data.

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

Initial Whole-Genome Sequencing of Plasma Cell Neoplasms in First Responders and Recovery Workers Exposed to the World Trade Center Attack of September 11, 2001.

PURPOSE: The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole-genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack. EXPERIMENTAL DESIGN: We performed WGS of nine CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster. RESULTS: No significant differences were observed in comparing the post-WTC driver and mutational signature landscapes with 110 previously published WGSs from 56 patients with multiple myeloma and the CoMMpass WGS cohort (n = 752). Leveraging constant activity of the single-base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure. CONCLUSIONS: Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.