Somatic inactivating PTPRJ mutations and dysregulated pathways identified in canine malignant melanoma by integrated comparative genomic analysis.

Canine malignant melanoma, a significant cause of mortality in domestic dogs, is a powerful comparative model for human melanoma, but little is known about its genetic etiology. We mapped the genomic landscape of canine melanoma through multi-platform analysis of 37 tumors (31 mucosal, 3 acral, 2 cutaneous, and 1 uveal) and 17 matching constitutional samples including long- and short-insert whole genome sequencing, RNA sequencing, array comparative genomic hybridization, single nucleotide polymorphism array, and targeted Sanger sequencing analyses. We identified novel predominantly truncating mutations in the putative tumor suppressor gene PTPRJ in 19% of cases. No BRAF mutations were detected, but activating RAS mutations (24% of cases) occurred in conserved hotspots in all cutaneous and acral and 13% of mucosal subtypes. MDM2 amplifications (24%) and TP53 mutations (19%) were mutually exclusive. Additional low-frequency recurrent alterations were observed amidst low point mutation rates, an absence of ultraviolet light mutational signatures, and an abundance of copy number and structural alterations. Mutations that modulate cell proliferation and cell cycle control were common and highlight therapeutic axes such as MEK and MDM2 inhibition. This mutational landscape resembles that seen in BRAF wild-type and sun-shielded human melanoma subtypes. Overall, these data inform biological comparisons between canine and human melanoma while suggesting actionable targets in both species.

A somatic reference standard for cancer genome sequencing.

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

Toward precision medicine in glioblastoma: the promise and the challenges.

Integrated sequencing strategies have provided a broader understanding of the genomic landscape and molecular classifications of multiple cancer types and have identified various therapeutic opportunities across cancer subsets. Despite pivotal advances in the characterization of genomic alterations in glioblastoma, targeted agents have shown minimal efficacy in clinical trials to date, and patient survival remains poor. In this review, we highlight potential reasons why targeting single alterations has yielded limited clinical efficacy in glioblastoma, focusing on issues of tumor heterogeneity and pharmacokinetic failure. We outline strategies to address these challenges in applying precision medicine to glioblastoma and the rationale for applying targeted combination therapy approaches that match genomic alterations with compounds accessible to the central nervous system.

Personalized treatment of Sézary syndrome by targeting a novel CTLA4:CD28 fusion.

Matching molecularly targeted therapies with cancer subtype-specific gene mutations is revolutionizing oncology care. However, for rare cancers this approach is problematic due to the often poor understanding of the disease's natural history and phenotypic heterogeneity, making treatment of these cancers a particularly unmet medical need in clinical oncology. Advanced Sézary syndrome (SS), an aggressive, exceedingly rare variant of cutaneous T-cell lymphoma (CTCL) is a prototypical example of a rare cancer. Through whole genome and RNA sequencing (RNA-seq) of a SS patient's tumor we discovered a highly expressed gene fusion between CTLA4 (cytotoxic T lymphocyte antigen 4) and CD28 (cluster of differentiation 28), predicting a novel stimulatory molecule on the surface of tumor T cells. Treatment with the CTLA4 inhibitor ipilimumab resulted in a rapid clinical response. Our findings suggest a novel driver mechanism for SS, and cancer in general, and exemplify an emerging model of cancer treatment using exploratory genomic analysis to identify a personally targeted treatment option when conventional therapies are exhausted.

Genome-wide analysis of hepatic lipid content in extreme obesity.

AIMS: Individuals with type 2 diabetes have an increased risk of developing non-alcoholic fatty liver disease (NAFLD), and NAFLD patients are also at greater risk for developing type 2 diabetes. Although the relationship between type 2 diabetes and NAFLD is highly interconnected, the pathogenic mechanisms linking the two diseases are poorly understood. The goal of this study was to identify genetic determinants of hepatic lipid accumulation through association analysis using histological phenotypes in obese individuals. METHODS: Using the Illumina HumanOmniExpress BeadChip assay, we genotyped 2,300 individuals on whom liver biopsy data were available. RESULTS: We analyzed total bilirubin levels, which are linked to fatty liver in severe obesity, and observed the strongest evidence for association with rs4148325 in UGT1A (P < 5.0 × 10(-93)), replicating previous findings. We assessed hepatic fat level and found strong evidence for association with rs4823173, rs2896019, and rs2281135, all located in PNPLA3 and rs10401969 in SUGP1. Analysis of liver transcript levels of 20 genes residing at the SUGP1/NCAN locus identified a 1.6-fold change in the expression of the LPAR2 gene in fatty liver. We also observed suggestive evidence for association between low-grade fat accumulation and rs10859525 and rs1294908, located upstream from SOCS2 and RAMP3, respectively. SOCS2 was differentially expressed between fatty and normal liver. CONCLUSIONS: These results replicate findings for several hepatic phenotypes in the setting of extreme obesity and implicate new loci that may play a role in the pathophysiology of hepatic lipid accumulation.

Open-access synthetic spike-in mRNA-seq data for cancer gene fusions.

BACKGROUND: Oncogenic fusion genes underlie the mechanism of several common cancers. Next-generation sequencing based RNA-seq analyses have revealed an increasing number of recurrent fusions in a variety of cancers. However, absence of a publicly available gene-fusion focused RNA-seq data impedes comparative assessment and collaborative development of novel gene fusions detection algorithms. We have generated nine synthetic poly-adenylated RNA transcripts that correspond to previously reported oncogenic gene fusions. These synthetic RNAs were spiked at known molarity over a wide range into total RNA prior to construction of next-generation sequencing mRNA libraries to generate RNA-seq data. RESULTS: Leveraging a priori knowledge about replicates and molarity of each synthetic fusion transcript, we demonstrate utility of this dataset to compare multiple gene fusion algorithms' detection ability. In general, more fusions are detected at higher molarity, indicating that our constructs performed as expected. However, systematic detection differences are observed based on molarity or algorithm-specific characteristics. Fusion-sequence specific detection differences indicate that for applications where specific sequences are being investigated, additional constructs may be added to provide quantitative data that is specific for the sequence of interest. CONCLUSIONS: To our knowledge, this is the first publicly available synthetic RNA-seq data that specifically leverages known cancer gene-fusions. The proposed method of designing multiple gene-fusion constructs over a wide range of molarity allows granular performance analyses of multiple fusion-detection algorithms. The community can leverage and augment this publicly available data to further collaborative development of analytical tools and performance assessment frameworks for gene fusions from next-generation sequencing data.

Profiles of extracellular miRNA in cerebrospinal fluid and serum from patients with Alzheimer's and Parkinson's diseases correlate with disease status and features of pathology.

The discovery and reliable detection of markers for neurodegenerative diseases have been complicated by the inaccessibility of the diseased tissue--such as the inability to biopsy or test tissue from the central nervous system directly. RNAs originating from hard to access tissues, such as neurons within the brain and spinal cord, have the potential to get to the periphery where they can be detected non-invasively. The formation and extracellular release of microvesicles and RNA binding proteins have been found to carry RNA from cells of the central nervous system to the periphery and protect the RNA from degradation. Extracellular miRNAs detectable in peripheral circulation can provide information about cellular changes associated with human health and disease. In order to associate miRNA signals present in cell-free peripheral biofluids with neurodegenerative disease status of patients with Alzheimer's and Parkinson's diseases, we assessed the miRNA content in cerebrospinal fluid and serum from postmortem subjects with full neuropathology evaluations. We profiled the miRNA content from 69 patients with Alzheimer's disease, 67 with Parkinson's disease and 78 neurologically normal controls using next generation small RNA sequencing (NGS). We report the average abundance of each detected miRNA in cerebrospinal fluid and in serum and describe 13 novel miRNAs that were identified. We correlated changes in miRNA expression with aspects of disease severity such as Braak stage, dementia status, plaque and tangle densities, and the presence and severity of Lewy body pathology. Many of the differentially expressed miRNAs detected in peripheral cell-free cerebrospinal fluid and serum were previously reported in the literature to be deregulated in brain tissue from patients with neurodegenerative disease. These data indicate that extracellular miRNAs detectable in the cerebrospinal fluid and serum are reflective of cell-based changes in pathology and can be used to assess disease progression and therapeutic efficacy.

Integrated genomic and epigenomic analysis of breast cancer brain metastasis.

The brain is a common site of metastatic disease in patients with breast cancer, which has few therapeutic options and dismal outcomes. The purpose of our study was to identify common and rare events that underlie breast cancer brain metastasis. We performed deep genomic profiling, which integrated gene copy number, gene expression and DNA methylation datasets on a collection of breast brain metastases. We identified frequent large chromosomal gains in 1q, 5p, 8q, 11q, and 20q and frequent broad-level deletions involving 8p, 17p, 21p and Xq. Frequently amplified and overexpressed genes included ATAD2, BRAF, DERL1, DNMTRB and NEK2A. The ATM, CRYAB and HSPB2 genes were commonly deleted and underexpressed. Knowledge mining revealed enrichment in cell cycle and G2/M transition pathways, which contained AURKA, AURKB and FOXM1. Using the PAM50 breast cancer intrinsic classifier, Luminal B, Her2+/ER negative, and basal-like tumors were identified as the most commonly represented breast cancer subtypes in our brain metastasis cohort. While overall methylation levels were increased in breast cancer brain metastasis, basal-like brain metastases were associated with significantly lower levels of methylation. Integrating DNA methylation data with gene expression revealed defects in cell migration and adhesion due to hypermethylation and downregulation of PENK, EDN3, and ITGAM. Hypomethylation and upregulation of KRT8 likely affects adhesion and permeability. Genomic and epigenomic profiling of breast brain metastasis has provided insight into the somatic events underlying this disease, which have potential in forming the basis of future therapeutic strategies.

Long insert whole genome sequencing for copy number variant and translocation detection.

As next-generation sequencing continues to have an expanding presence in the clinic, the identification of the most cost-effective and robust strategy for identifying copy number changes and translocations in tumor genomes is needed. We hypothesized that performing shallow whole genome sequencing (WGS) of 900-1000-bp inserts (long insert WGS, LI-WGS) improves our ability to detect these events, compared with shallow WGS of 300-400-bp inserts. A priori analyses show that LI-WGS requires less sequencing compared with short insert WGS to achieve a target physical coverage, and that LI-WGS requires less sequence coverage to detect a heterozygous event with a power of 0.99. We thus developed an LI-WGS library preparation protocol based off of Illumina's WGS library preparation protocol and illustrate the feasibility of performing LI-WGS. We additionally applied LI-WGS to three separate tumor/normal DNA pairs collected from patients diagnosed with different cancers to demonstrate our application of LI-WGS on actual patient samples for identification of somatic copy number alterations and translocations. With the evolution of sequencing technologies and bioinformatics analyses, we show that modifications to current approaches may improve our ability to interrogate cancer genomes.

Detection of an atypical teratoid rhabdoid brain tumor gene deletion in circulating blood using next-generation sequencing.

Circulating biomarkers such as somatic chromosome mutations are novel diagnostic tools to detect cancer noninvasively. We describe focal deletions found in a patient with atypical teratoid rhabdoid tumor, a highly aggressive early childhood pediatric tumor. First, we used magnetic resonance imaging (MRI) and histopathology to study the tumor anatomy. Next, we used whole genome sequencing (Next Gen Sequencing) and Bioinformatics interrogation to discover the presence of 3 focal deletions in tumor tissue and 2 of these 3 focal deletions in patient's blood also. About 20% of the blood DNA sequencing reads matched the tumor DNA reads at the SMARCB1 gene locus. Circulating, tumor-specific DNA aberrations are a promising biomarker for atypical teratoid rhabdoid tumor patients. The high percentage of tumor DNA detected in blood indicates that either circulating brain tumor cells lyse in the blood or that contents of brain tumor cells traverse a possibly compromised blood-brain barrier in this patient.

Loss of the tumor suppressor SMARCA4 in small cell carcinoma of the ovary, hypercalcemic type (SCCOHT).

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT), is a rare and understudied cancer with a dismal prognosis. SCCOHT's infrequency has hindered empirical study of its biology and clinical management. However, we and others have recently identified inactivating mutations in the SWI/SNF chromatin remodeling gene SMARCA4 with concomitant loss of SMARCA4 protein in the majority of SCCOHT tumors.(1-4) Here we summarize these findings and report SMARCA4 status by targeted sequencing and/or immunohistochemistry (IHC) in an additional 12 SCCOHT tumors, 3 matched germlines, and the cell line SCCOHT-1. We also report the identification of a homozygous inactivating mutation in the gene SMARCB1 in one SCCOHT tumor with wild-type SMARCA4, suggesting that SMARCB1 inactivation may also play a role in the pathogenesis of SCCOHT. To date, SMARCA4 mutations and protein loss have been reported in the majority of 69 SCCOHT cases (including 2 cell lines). These data firmly establish SMARCA4 as a tumor suppressor whose loss promotes the development of SCCOHT, setting the stage for rapid advancement in the biological understanding, diagnosis, and treatment of this rare tumor type.

The Marker State Space (MSS) method for classifying clinical samples.

The development of accurate clinical biomarkers has been challenging in part due to the diversity between patients and diseases. One approach to account for the diversity is to use multiple markers to classify patients, based on the concept that each individual marker contributes information from its respective subclass of patients. Here we present a new strategy for developing biomarker panels that accounts for completely distinct patient subclasses. Marker State Space (MSS) defines "marker states" based on all possible patterns of high and low values among a panel of markers. Each marker state is defined as either a case state or a control state, and a sample is classified as case or control based on the state it occupies. MSS was used to define multi-marker panels that were robust in cross validation and training-set/test-set analyses and that yielded similar classification accuracy to several other classification algorithms. A three-marker panel for discriminating pancreatic cancer patients from control subjects revealed subclasses of patients based on distinct marker states. MSS provides a straightforward approach for modeling highly divergent subclasses of patients, which may be adaptable for diverse applications.

Comparative RNA-Seq and microarray analysis of gene expression changes in B-cell lymphomas of Canis familiaris.

Comparative oncology is a developing research discipline that is being used to assist our understanding of human neoplastic diseases. Companion canines are a preferred animal oncology model due to spontaneous tumor development and similarity to human disease at the pathophysiological level. We use a paired RNA sequencing (RNA-Seq)/microarray analysis of a set of four normal canine lymph nodes and ten canine lymphoma fine needle aspirates to identify technical biases and variation between the technologies and convergence on biological disease pathways. Surrogate Variable Analysis (SVA) provides a formal multivariate analysis of the combined RNA-Seq/microarray data set. Applying SVA to the data allows us to decompose variation into contributions associated with transcript abundance, differences between the technology, and latent variation within each technology. A substantial and highly statistically significant component of the variation reflects transcript abundance, and RNA-Seq appeared more sensitive for detection of transcripts expressed at low levels. Latent random variation among RNA-Seq samples is also distinct in character from that impacting microarray samples. In particular, we observed variation between RNA-Seq samples that reflects transcript GC content. Platform-independent variable decomposition without a priori knowledge of the sources of variation using SVA represents a generalizable method for accomplishing cross-platform data analysis. We identified genes differentially expressed between normal lymph nodes of disease free dogs and a subset of the diseased dogs diagnosed with B-cell lymphoma using each technology. There is statistically significant overlap between the RNA-Seq and microarray sets of differentially expressed genes. Analysis of overlapping genes in the context of biological systems suggests elevated expression and activity of PI3K signaling in B-cell lymphoma biopsies compared with normal biopsies, consistent with literature describing successful use of drugs targeting this pathway in lymphomas.

Comparison of Analysis Tools for miRNA High Throughput Sequencing Using Nerve Crush as a Model.

Recent advances in sample preparation and analysis for next generation sequencing have made it possible to profile and discover new miRNAs in a high throughput manner. In the case of neurological disease and injury, these types of experiments have been more limited. Possibly because tissues such as the brain and spinal cord are inaccessible for direct sampling in living patients, and indirect sampling of blood and cerebrospinal fluid are affected by low amounts of RNA. We used a mouse model to examine changes in miRNA expression in response to acute nerve crush. We assayed miRNA from both muscle tissue and blood plasma. We examined how the depth of coverage (the number of mapped reads) changed the number of detectable miRNAs in each sample type. We also found that samples with very low starting amounts of RNA (mouse plasma) made high depth of mature miRNA coverage more difficult to obtain. Each tissue must be assessed independently for the depth of coverage required to adequately power detection of differential expression, weighed against the cost of sequencing that sample to the adequate depth. We explored the changes in total mapped reads and differential expression results generated by three different software packages: miRDeep2, miRNAKey, and miRExpress and two different analysis packages, DESeq and EdgeR. We also examine the accuracy of using miRDeep2 to predict novel miRNAs and subsequently detect them in the samples using qRT-PCR.

Identification of extracellular miRNA in human cerebrospinal fluid by next-generation sequencing.

There has been a growing interest in using next-generation sequencing (NGS) to profile extracellular small RNAs from the blood and cerebrospinal fluid (CSF) of patients with neurological diseases, CNS tumors, or traumatic brain injury for biomarker discovery. Small sample volumes and samples with low RNA abundance create challenges for downstream small RNA sequencing assays. Plasma, serum, and CSF contain low amounts of total RNA, of which small RNAs make up a fraction. The purpose of this study was to maximize RNA isolation from RNA-limited samples and apply these methods to profile the miRNA in human CSF by small RNA deep sequencing. We systematically tested RNA isolation efficiency using ten commercially available kits and compared their performance on human plasma samples. We used RiboGreen to quantify total RNA yield and custom TaqMan assays to determine the efficiency of small RNA isolation for each of the kits. We significantly increased the recovery of small RNA by repeating the aqueous extraction during the phenol-chloroform purification in the top performing kits. We subsequently used the methods with the highest small RNA yield to purify RNA from CSF and serum samples from the same individual. We then prepared small RNA sequencing libraries using Illumina's TruSeq sample preparation kit and sequenced the samples on the HiSeq 2000. Not surprisingly, we found that the miRNA expression profile of CSF is substantially different from that of serum. To our knowledge, this is the first time that the small RNA fraction from CSF has been profiled using next-generation sequencing.

Effect of selection of QTc formula on eligibility of cancer patients for phase I clinical trials.

A retrospective analysis of 130 patients was conducted in a Phase I oncology clinic to assess the effect of QTc formula selection on clinical trial eligibility. QTc values were calculated from screening electrocardiograms using 7 formulae (Bazett, Fridericia, Framingham, Hodges, Mayeda, Van de Water and Wohlfart). QTc values > 470 ms for females and > 450 ms for males were used to define prolongation. Concomitant medication potential for QTc prolongation was determined using a public database (AzCert). Ineligibility rates ranged from 3.1 % to 17.7 % (Framingham: 3.1 %, Van de Water: 3.1 %, Hodges: 3.1 %, Wohlfart: 3.1 %, Fridericia: 3.9 %, Bazett: 10.8 % and Mayeda: 17.7 %). A consistent ineligibility rate was achieved by using formulae-specific thresholds. Fifty one percent of patients were taking concomitant medications with QTc prolongation potential. The proportion of concomitant medications with the potential to prolong QTc was 11.57 % (96 of 830). Uniform criteria and guidelines for selection of QTc formulae need to be developed. Formulae-specific QTc thresholds also need to be specified.

Deep clonal profiling of formalin fixed paraffin embedded clinical samples.

Formalin fixed paraffin embedded (FFPE) tissues are a vast resource of annotated clinical samples. As such, they represent highly desirable and informative materials for the application of high definition genomics for improved patient management and to advance the development of personalized therapeutics. However, a limitation of FFPE tissues is the variable quality of DNA extracted for analyses. Furthermore, admixtures of non-tumor and polyclonal neoplastic cell populations limit the number of biopsies that can be studied and make it difficult to define cancer genomes in patient samples. To exploit these valuable tissues we applied flow cytometry-based methods to isolate pure populations of tumor cell nuclei from FFPE tissues and developed a methodology compatible with oligonucleotide array CGH and whole exome sequencing analyses. These were used to profile a variety of tumors (breast, brain, bladder, ovarian and pancreas) including the genomes and exomes of matching fresh frozen and FFPE pancreatic adenocarcinoma samples.

Genome-wide characterization of pancreatic adenocarcinoma patients using next generation sequencing.

Pancreatic adenocarcinoma (PAC) is among the most lethal malignancies. While research has implicated multiple genes in disease pathogenesis, identification of therapeutic leads has been difficult and the majority of currently available therapies provide only marginal benefit. To address this issue, our goal was to genomically characterize individual PAC patients to understand the range of aberrations that are occurring in each tumor. Because our understanding of PAC tumorigenesis is limited, evaluation of separate cases may reveal aberrations, that are less common but may provide relevant information on the disease, or that may represent viable therapeutic targets for the patient. We used next generation sequencing to assess global somatic events across 3 PAC patients to characterize each patient and to identify potential targets. This study is the first to report whole genome sequencing (WGS) findings in paired tumor/normal samples collected from 3 separate PAC patients. We generated on average 132 billion mappable bases across all patients using WGS, and identified 142 somatic coding events including point mutations, insertion/deletions, and chromosomal copy number variants. We did not identify any significant somatic translocation events. We also performed RNA sequencing on 2 of these patients' tumors for which tumor RNA was available to evaluate expression changes that may be associated with somatic events, and generated over 100 million mapped reads for each patient. We further performed pathway analysis of all sequencing data to identify processes that may be the most heavily impacted from somatic and expression alterations. As expected, the KRAS signaling pathway was the most heavily impacted pathway (P<0.05), along with tumor-stroma interactions and tumor suppressive pathways. While sequencing of more patients is needed, the high resolution genomic and transcriptomic information we have acquired here provides valuable information on the molecular composition of PAC and helps to establish a foundation for improved therapeutic selection.

Paired tumor and normal whole genome sequencing of metastatic olfactory neuroblastoma.

BACKGROUND: Olfactory neuroblastoma (ONB) is a rare cancer of the sinonasal tract with little molecular characterization. We performed whole genome sequencing (WGS) on paired normal and tumor DNA from a patient with metastatic-ONB to identify the somatic alterations that might be drivers of tumorigenesis and/or metastatic progression. METHODOLOGY/PRINCIPAL FINDINGS: Genomic DNA was isolated from fresh frozen tissue from a metastatic lesion and whole blood, followed by WGS at >30X depth, alignment and mapping, and mutation analyses. Sanger sequencing was used to confirm selected mutations. Sixty-two somatic short nucleotide variants (SNVs) and five deletions were identified inside coding regions, each causing a non-synonymous DNA sequence change. We selected seven SNVs and validated them by Sanger sequencing. In the metastatic ONB samples collected several months prior to WGS, all seven mutations were present. However, in the original surgical resection specimen (prior to evidence of metastatic disease), mutations in KDR, MYC, SIN3B, and NLRC4 genes were not present, suggesting that these were acquired with disease progression and/or as a result of post-treatment effects. CONCLUSIONS/SIGNIFICANCE: This work provides insight into the evolution of ONB cancer cells and provides a window into the more complex factors, including tumor clonality and multiple driver mutations.

Whole-genome sequencing of multiple myeloma from diagnosis to plasma cell leukemia reveals genomic initiating events, evolution, and clonal tides.

The longitudinal evolution of a myeloma genome from diagnosis to plasma cell leukemia has not previously been reported. We used whole-genome sequencing (WGS) on 4 purified tumor samples and patient germline DNA drawn over a 5-year period in a t(4;14) multiple myeloma patient. Tumor samples were acquired at diagnosis, first relapse, second relapse, and end-stage secondary plasma cell leukemia (sPCL). In addition to the t(4;14), all tumor time points also shared 10 common single-nucleotide variants (SNVs) on WGS comprising shared initiating events. Interestingly, we observed genomic sequence variants that waxed and waned with time in progressive tumors, suggesting the presence of multiple independent, yet related, clones at diagnosis that rose and fell in dominance. Five newly acquired SNVs, including truncating mutations of RB1 and ZKSCAN3, were observed only in the final sPCL sample suggesting leukemic transformation events. This longitudinal WGS characterization of the natural history of a high-risk myeloma patient demonstrated tumor heterogeneity at diagnosis with shifting dominance of tumor clones over time and has also identified potential mutations contributing to myelomagenesis as well as transformation from myeloma to overt extramedullary disease such as sPCL.