scholarly journals Comprehensive Mutational Analysis of Juvenile Myelomonocytic Leukemia Using Whole-Genome Sequencing

Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 2974-2974
Author(s):  
Yusuke Okuno ◽  
Hideki Muramatsu ◽  
Norihiro Murakami ◽  
Nozomu Kawashima ◽  
Manabu Wakamatsu ◽  
...  

Background Juvenile myelomonocytic leukemia (JMML) is a rare and exclusively pediatric myelodysplastic/myeloproliferative neoplasm. This disease is genetically characterized by an extremely small number of somatic mutations (an average of 0.8 mutations/exome/patient). It has been shown that causative somatic and/or germline mutations activating the RAS pathway are located in PTPN11, NF1, NRAS, KRAS, and CBL in 85% of patients with JMML. Furthermore, up to 20% of the patients have additional secondary mutations including SETBP1, and JAK3 mutations. In 2% of the patients, we identified, by RNA sequencing, activating kinase lesions affecting ALK or ROS1. Such findings suggest that other kinase fusions are present in JMML. There is an exceptional scarcity of somatic passenger mutations on the exome, suggesting that a small number of driver mutations drive JMML. However, to date, this hypothesis has not been investigated by whole-genome sequencing. Patients and Methods We performed a whole-genome sequencing (WGS) study in 48 patients with JMML. Bone marrow specimens and in vitro-cultured T cells were used as tumor and germline samples, respectively. Next-generation sequencing was performed using a HiSeq X platform (Illumina). Data analysis was performed by our in-house pipeline. Specifically, the pipeline detects single nucleotide variants (SNVs), copy number variants, somatic loss of heterozygosity (LOH), and chromosomal structural variations (SVs). The study was approved by the institutional review board of Nagoya University Graduate School of Medicine. Results In each patient we detected an average of 28 somatic mutations. These were primarily C-to-T transition in the CpG context, indicating that the mutations occurred by cell division. Besides RAS pathway and known secondary mutations, we observed no significant accumulation of somatic mutations in either coding or non-coding regions. Although we detected RAS pathway mutations in 90% of the patients, all mutations were on exome. However, we identified germline microdeletions affecting CBL and NF1, which had not been identified by exome sequencing. Additionally, we found two LOH events that affected NF1. Bi-allelic inactivation of NF1 is generally observed in patients with JMML; however, no pathogenic SNVs were identified in these two patients. We identified two chromosomal translocations that caused activating kinase lesions (i.e., RANBP2-ALK and TBL1XR1-ROS1). These had been pointed out in our previous RNA sequencing study. Another patient carried a complex SV that affected XPO1 (encoding exportin 1 or chromosome region maintenance 1 protein homolog). Although fusion genes involving XPO1 are reported to be present in lymphoid malignancies, the role of this SV in JMML remains unclear. Conclusions JMML is characterized by driver mutations that are largely present within the exome. However, WGS can still play a role in identifying both coding and non-coding mutations. LOH events without pathogenic SNVs suggest the presence of novel regulatory mechanisms of NF1. Conclusively, JMML is characterized by a paucity of somatic alterations and driver mutations. Hence, current research efforts should focus on RAS pathway mutations and known secondary mutations, many of which can be targeted. Disclosures No relevant conflicts of interest to declare.

2014 ◽  
Vol 5 (1) ◽  
Author(s):  
J. -B. Cazier ◽  
◽  
S. R. Rao ◽  
C. M. McLean ◽  
A. K. Walker ◽  
...  

Abstract Bladder cancers are a leading cause of death from malignancy. Molecular markers might predict disease progression and behaviour more accurately than the available prognostic factors. Here we use whole-genome sequencing to identify somatic mutations and chromosomal changes in 14 bladder cancers of different grades and stages. As well as detecting the known bladder cancer driver mutations, we report the identification of recurrent protein-inactivating mutations in CDKN1A and FAT1. The former are not mutually exclusive with TP53 mutations or MDM2 amplification, showing that CDKN1A dysfunction is not simply an alternative mechanism for p53 pathway inactivation. We find strong positive associations between higher tumour stage/grade and greater clonal diversity, the number of somatic mutations and the burden of copy number changes. In principle, the identification of sub-clones with greater diversity and/or mutation burden within early-stage or low-grade tumours could identify lesions with a high risk of invasive progression.


2021 ◽  
Vol 218 (2) ◽  
Author(s):  
Eleni Louka ◽  
Benjamin Povinelli ◽  
Alba Rodriguez-Meira ◽  
Gemma Buck ◽  
Wei Xiong Wen ◽  
...  

Juvenile myelomonocytic leukemia (JMML) is a poor-prognosis childhood leukemia usually caused by RAS-pathway mutations. The cellular hierarchy in JMML is poorly characterized, including the identity of leukemia stem cells (LSCs). FACS and single-cell RNA sequencing reveal marked heterogeneity of JMML hematopoietic stem/progenitor cells (HSPCs), including an aberrant Lin−CD34+CD38−CD90+CD45RA+ population. Single-cell HSPC index-sorting and clonogenic assays show that (1) all somatic mutations can be backtracked to the phenotypic HSC compartment, with RAS-pathway mutations as a “first hit,” (2) mutations are acquired with both linear and branching patterns of clonal evolution, and (3) mutant HSPCs are present after allogeneic HSC transplant before molecular/clinical evidence of relapse. Stem cell assays reveal interpatient heterogeneity of JMML LSCs, which are present in, but not confined to, the phenotypic HSC compartment. RNA sequencing of JMML LSC reveals up-regulation of stem cell and fetal genes (HLF, MEIS1, CNN3, VNN2, and HMGA2) and candidate therapeutic targets/biomarkers (MTOR, SLC2A1, and CD96), paving the way for LSC-directed disease monitoring and therapy in this disease.


Neurosurgery ◽  
2019 ◽  
Vol 66 (Supplement_1) ◽  
Author(s):  
Tej D Azad ◽  
Ming Zhang ◽  
Rajiv Iyer ◽  
Qing Wang ◽  
Tomas Garzon-Muvdi ◽  
...  

Abstract INTRODUCTION Intramedullary spinal cord tumors (IMSCTs) are a rare, heterogeneous group of neoplasms with limited treatment options and high rates of morbidity and mortality. Next-generation sequencing has revealed opportunities for targeted therapies of the intracranial counterparts of IMSCT, but little is known about the molecular features of IMSCT. METHODS To better understand the genetic basis of these tumors we performed whole exome sequencing on fifty-one IMSCT and matched germline DNA, including 29 ependymomas, 16 astrocytomas, 4 gangliogliomas,1hemangioblastoma, and 1 oligodendroglioma. Whole-genome sequencing was further performed on 12 IMSCT to discover possible structural variants. RESULTS Though recurrent somatic mutations in IMSCTs were rare, we identified NF2 mutations in 15.7% of tumors (ependymoma, N = 7; astrocytoma, N = 1), RP1 mutations in 5.9% of tumors (ependymoma, N = 3), and ESX1 mutations in 5.9% of tumors (ependymoma, N = 3). We further identified copy number amplifications in CTU1 in 25% of myxopapillary ependymomas. Given the paucity of somatic driver mutations, we further performed whole-genome sequencing of 12 tumors (ependymoma, N = 9; astrocytoma, N = 3). Overall, we observed that IMSCTs with intracranial histologic counterparts did not harbor the canonical mutations associated with their intracranial counterparts (eg glioblastoma). CONCLUSION Our findings suggest that the origin of IMSCTs may be distinct from tumors arising within other compartments of the central nervous system and provides a framework to begin more biologically based therapeutic strategies.


2020 ◽  
Vol 11 (1) ◽  
Author(s):  
Lydia Y. Liu ◽  
Vinayak Bhandari ◽  
Adriana Salcedo ◽  
Shadrielle M. G. Espiritu ◽  
Quaid D. Morris ◽  
...  

AbstractWhole-genome sequencing can be used to estimate subclonal populations in tumours and this intra-tumoural heterogeneity is linked to clinical outcomes. Many algorithms have been developed for subclonal reconstruction, but their variabilities and consistencies are largely unknown. We evaluate sixteen pipelines for reconstructing the evolutionary histories of 293 localized prostate cancers from single samples, and eighteen pipelines for the reconstruction of 10 tumours with multi-region sampling. We show that predictions of subclonal architecture and timing of somatic mutations vary extensively across pipelines. Pipelines show consistent types of biases, with those incorporating SomaticSniper and Battenberg preferentially predicting homogenous cancer cell populations and those using MuTect tending to predict multiple populations of cancer cells. Subclonal reconstructions using multi-region sampling confirm that single-sample reconstructions systematically underestimate intra-tumoural heterogeneity, predicting on average fewer than half of the cancer cell populations identified by multi-region sequencing. Overall, these biases suggest caution in interpreting specific architectures and subclonal variants.


2019 ◽  
Vol 49 (3) ◽  
pp. 115
Author(s):  
Jung Heon Kim ◽  
Jiyeon Kim ◽  
Bon-Sang Koo ◽  
Hanseul Oh ◽  
Jung-Joo Hong ◽  
...  

2018 ◽  
Vol 39 (4) ◽  
pp. 579-587 ◽  
Author(s):  
Alessandro Coppe ◽  
Leonardo Nogara ◽  
Matteo Samuele Pizzuto ◽  
Alice Cani ◽  
Simone Cesaro ◽  
...  

Blood ◽  
2009 ◽  
Vol 114 (22) ◽  
pp. 3965-3965
Author(s):  
Lukas D. Wartman ◽  
Li Ding ◽  
David E. Larson ◽  
Michael D. McLellan ◽  
Heather Schmidt ◽  
...  

Abstract Abstract 3965 Poster Board III-901 We have recently established that whole genome sequencing is a valid, unbiased approach that can identify novel candidate mutations that may be important for AML pathogenesis (Ley et al Nature 2008, Mardis et al NEJM 2009). Acute promyelocytic leukemia (APL, FAB M3 AML) is a subtype of AML characterized by the t(15;17)(q22;q11.2) translocation that creates an oncogenic fusion gene, PML-RARA. Our laboratory has previously modeled APL in a mouse in an effort to understand the genetic events that lead to the disease. In our knockin mouse model, a human PML-RARA cDNA was targeted to the 5' untranslated region of the mouse cathepsin G gene on chromosome 14 (mCG-PR). The targeting vector was transfected into the RW-4 embryonic stem cell line, derived from a 129/SvJ mouse. The transfected RW-4 cells were injected into C57Bl/6 blastocysts, and chimeric offspring were bred to C57Bl/6 mice. F1 129/SvJ x C57Bl/6 mice were subsequently backcrossed onto the B6/Taconic background for 10 generations before establishing a tumor watch. About 60% of the mCG-PR mice in the Bl/6 background develop a disease that closely resembles APL only after a latent period of 7-18 months, suggesting that additional progression mutations are required for APL development. Array-based genomic techniques (expression array studies and high resolution CGH) have revealed some recurring genetic alterations that may be relevant for progression (i.e. an interstitial deletion of chromosome 2, trisomy 15, etc.), but gene-specific progression mutations have not yet been identified. To begin to identify these mutations in an unbiased fashion, we sequenced a cytogenetically normal, diploid mouse APL genome using massively parallel DNA sequencing via the Illumina platform. Since the tumor arose in a highly inbred mouse strain, we predicted that 15x coverage of the genome (approximately 40 billion base pairs of sequence) would be necessary to identify >90% of the heterozygous somatic mutations. We generated 2 Illumina paired-end libraries (insert sizes of 300-350 bp and 550-600 bp) and generated 59.64 billion base pairs of sequence with 3 full sequencing runs; the reads that successfully mapped generated 15.6x coverage. The sequence data predicted 87,778 heterozygous Single Nucleotide Variants (SNVs) compared to the mouse C57Bl6/J reference sequence, and 23,439 homozygous SNVs. Of the predicted heterozygous SNVs, 695 were non-synonymous (missense or nonsense, or altering a canonical splice site). Thus far, 80 of these putative non-synonymous SNVs have been further analyzed using Sanger sequencing of the original tumor DNA vs. pooled B6/Taconic spleen DNA and pooled129/SvJ spleen DNA as controls. 37/80 were shown to be false positive calls, and 37 were inherited SNPs from residual regions of the129/SvJ genome. 6/80 were present only in the tumor genome, and were candidate somatic mutations. These 6 were screened in 89 additional murine APL tumor samples derived from the same mouse model. Mutations in the Jarid2 (L915I) and Capns2 (N149S) genes occurred only in the proband, and are therefore of uncertain significance. 4/6 mutations were found in additional samples; 3 of these mutations were derived from a common ancestor of the proband and the other affected mice, and were therefore not relevant for pathogenesis. The other recurring mutation was in the pseudokinase domain of JAK1 (V657F), and was identified in one other mouse that was not closely related to the proband. This mutation is orthologous to the known activating mutation V617F in human JAK2, and is identical to a recently described JAK1 pseudokinase domain mutation (V658F) found in human APL and T-ALL samples (EG Jeong et al, Clin Can Res 14: 3716, 2008). We are currently testing the functional significance of this mutation by expressing it in bone marrow cells derived from young WT vs. mCG-PR mice. In summary, unbiased whole genome sequencing of a mouse APL genome has identified a recurring mutation of JAK1 found in both human and mouse APL samples. This approach may allow us to rapidly identify progression mutations that are common to human and murine AML, and provides an important proof-of-concept that this mouse model of AML is functionally related to its human counterpart. Disclosures: No relevant conflicts of interest to declare.


Blood ◽  
2011 ◽  
Vol 118 (21) ◽  
pp. 434-434 ◽  
Author(s):  
Zachary Hunter ◽  
Lian Xu ◽  
Yangsheng Zhou ◽  
Guang Yang ◽  
Xia Liu ◽  
...  

Abstract Abstract 434 Introduction: The primary oncogenetic event resulting in malignant transformation in Waldenstrom's Macroglobulinemia (WM) remains to be delineated. We therefore employed whole genome sequencing (WGS) to help identify potential somatic variants in WM. Patients and Methods: Thirty patients meeting consensus criteria for the diagnosis of WM were included for these studies, whose characteristics are depicted in Table 1. CD19-magnetic bead sorting was used for isolation of bone marrow LPC. CD19-depleted PB mononuclear cells were collected as matched normal tissue. For 10 patients, WGS of tumor and matched normal samples was performed, and for 20 additional patients tumor samples alone were completely sequenced. Library construction and WGS was performed by Complete Genomics Inc. Read sequences were aligned to the NCBI Build 37. High confidence somatic variants were identified using cgatools version 1.3. Novel non-synonymous exonic variants for familial and sporadic LPL/WM patients were identified using ANNOVAR using to filter against several large databases including dbSNP version 132, the November 2010 release version of the 1,000 genomes project, and a 46 healthy donor dataset from Complete Genomics, Inc. based on KnownGene annotations. Variants filtered out in this process were checked against the dbSNP132 flagged SNP database for potential clinical significance. Data was further annotated against the Database of Genomic Variants and the Segmental Duplication Database, TargetScan, and transcription factor binding site data from the ENCODE project. When applicable, variants were scored using SIFT, PolyPhen2 and Mutation Taster. Copy number neutral loss of heterozygosity (CNLOH) was identified from the rate of heterozygous variants per 500,000 base pairs, CG content adjusted coverage data, and allele imbalance calculated from the percentage of total reads supporting the less covered allele. Results: Tumor and normal genomes were both sequenced to an average of 66X (range 60–91X) coverage of mapped individual reads. The average gross mapped yield for these genomes was 186.89 (range 171.56–262.03 Gb). Acquired copy number changes were common, and included losses in chromosome 6q (13/30; 43%), gains in chromosome 4 (7/30; 23%), and gains in 6p (3/30; 10%). Large regions of CNLOH were observed in 9/30 (30%) of patients occurring in chromosomes 1, 2, 3, 5, 9, 11, 17, 21, and X. The most frequent somatic variant occurred at position 38182641 in chromosome 3p22.2 in the myeloid differentiation primary response (MYD88) gene, resulting in a non-synonymous change at amino acid position 265 from leucine to proline (L265P) in 26/30 (86.7%) patients. Of these, 4/26 (15%) had a CNLOH covering this position making the variant effectively homozygous. Additional somatic variants occurred in transporter 2, ATP-binding cassette, sub-family B (TAP2) gene in 7/30 (23%) patients; chemokine (C-X-C motif) receptor 4 (CXCR4) gene in 6/30 (20%) patients. Somatic variants were also identified in the coding regions of low density lipoprotein receptor-related protein 1B (LRP1B) gene in 5/30 (17%) patients; mesothelin (MSLN) gene in 4/30 (13%) patients; AT rich interactive domain 1A (ARID1A) gene in 3/30 (10%) patients; histone cluster 1, H1e (HIST1H1E) in 3/30 (10%) patients, and Rap guanine nucleotide exchange factor 3 (RAPGEF3) in 3/30 (10%) patients. Conclusions: The results of this study provide the first reporting of comprehensive WGS efforts in patients with WM, and reveal recurring somatic variants in genes with important regulatory functions including MYD88, TAP2, and CXCR4. Structural and functional validation studies are ongoing and will be updated at the meeting. The results of these studies provide important new insights into the pathogenesis of WM. Disclosures: No relevant conflicts of interest to declare.


Blood ◽  
2011 ◽  
Vol 118 (21) ◽  
pp. 404-404 ◽  
Author(s):  
John S. Welch ◽  
David Larson ◽  
Li Ding ◽  
Michael D. McLellan ◽  
Tamara Lamprecht ◽  
...  

Abstract Abstract 404 To characterize the genomic events associated with distinct subtypes of AML, we used whole genome sequencing to compare 24 tumor/normal sample pairs from patients with normal karyotype (NK) M1-AML (12 cases) and t(15;17)-positive M3-AML (12 cases). All single nucleotide variants (SNVs), small insertions and deletions (indels), and cryptic structural variants (SVs) identified by whole genome sequencing (average coverage 28x) were validated using sample-specific custom Nimblegen capture arrays, followed by Illumina sequencing; an average coverage of 972 reads per somatic variant yielded 10,597 validated somatic variants (average 421/genome). Of these somatic mutations, 308 occurred in 286 unique genes; on average, 9.4 somatic mutations per genome had translational consequences. Several important themes emerged: 1) AML genomes contain a diverse range of recurrent mutations. We assessed the 286 mutated genes for recurrency in an additional 34 NK M1-AML cases and 9 M3-AML cases. We identified 51 recurrently mutated genes, including 37 that had not previously been described in AML; on average, each genome had 3 recurrently mutated genes (M1 = 3.2; M3 = 2.8, p = 0.32). 2) Many recurring mutations cluster in mutually exclusive pathways, suggesting pathophysiologic importance. The most commonly mutated genes were: FLT3 (36%), NPM1 (25%), DNMT3A (21%), IDH1 (18%), IDH2 (10%), TET2 (10%), ASXL1 (6%), NRAS (6%), TTN (6%), and WT1 (6%). In total, 3 genes (excluding PML-RARA) were mutated exclusively in M3 cases. 22 genes were found only in M1 cases (suggestive of alternative initiating mutations which occurred in methylation, signal transduction, and cohesin complex genes). 25 genes were mutated in both M1 and M3 genomes (suggestive of common progression mutations relevant for both subtypes). A single mutation in a cell growth/signaling gene occurred in 38 of 67 cases (FLT3, NRAS, RUNX1, KIT, CACNA1E, CADM2, CSMD1); these mutations were mutually exclusive of one another, and many of them occurred in genomes with PML-RARA, suggesting that they are progression mutations. We also identified a new leukemic pathway: mutations were observed in all four genes that encode members of the cohesin complex (STAG2, SMC1A, SMC3, RAD21), which is involved in mitotic checkpoints and chromatid separation. The cohesin mutations were mutually exclusive of each other, and collectively occur in 10% of non-M3 AML patients. 3) AML genomes also contain hundreds of benign “passenger” mutations. On average 412 somatic mutations per genome were translationally silent or occurred outside of annotated genes. Both M1 and M3 cases had similar total numbers of mutations per genome, similar mutation types (which favored C>T/G>A transitions), and a similar random distribution of variants throughout the genome (which was affected neither by coding regions nor expression levels). This is consistent with our recent observations of random “passenger” mutations in hematopoietic stem cell (HSC) clones derived from normal patients (Ley et al manuscript in preparation), and suggests that most AML-associated mutations are not pathologic, but pre-existed in the HSC at the time of initial transformation. In both studies, the total number of SNVs per genome correlated positively with the age of the patient (R2 = 0.48, p = 0.001), providing a possible explanation for the increasing incidence of AML in elderly patients. 4) NK M1 and M3 AML samples are mono- or oligo-clonal. By comparing the frequency of all somatic mutations within each sample, we could identify clusters of mutations with similar frequencies (leukemic clones) and determined that the average number of clones per genome was 1.8 (M1 = 1.5; M3 = 2.2; p = 0.04). 5) t(15;17) is resolved by a non-homologous end-joining repair pathway, since nucleotide resolution of all 12 t(15;17) breakpoints revealed inconsistent micro-homologies (0 – 7 bp). Summary: These data provide a genome-wide overview of NK and t(15;17) AML and provide important new insights into AML pathogenesis. AML genomes typically contain hundreds of random, non-genic mutations, but only a handful of recurring mutated genes that are likely to be pathogenic because they cluster in mutually exclusive pathways; specific combinations of recurring mutations, as well as rare and private mutations, shape the leukemia phenotype in an individual patient, and help to explain the clinical heterogeneity of this disease. Disclosures: Westervelt: Novartis: Speakers Bureau.


Sign in / Sign up

Export Citation Format

Share Document