Whole Genome Sequence Analysis of Primary Myelofibrosis

• Published in: Blood 2012-11, Vol.120 (21), p.2863-2863
• Main Authors: Merker, Jason D., Roskin, Krishna, Ng, Dana, Pan, Cuiping, Fisk, Dianna G., Jones, Carol D., Gojenola, Linda, Clark, Michael J., Zhang, Bing, Cherry, Michael, Snyder, Michael, Boyd, Scott D., Zehnder, James L., Fire, Andrew Z., Gotlib, Jason
• Format: Article
• Language: English
• ISSN: 0006-4971; 1528-0020
• DOI: 10.1182/blood.V120.21.2863.2863

Abstract 2863 JAK2 V617F is a recurrent, activating mutation in patients (pts) with BCR-ABL1-negative MPNs. Mutations in codon 515 of MPL occur in 1–5% and 5–10% of ET and PMF pts, respectively, and similar to JAK2 V617F, lead to constitutive JAK-STAT signaling. The acquisition of multiple mutational events affecting the JAK-STAT axis or components of the epigenetic machinery is common in MPNs and likely contributes to phenotypic diversity, including progression to acute myeloid leukemia. Mutated genes thus far implicated in MPN initiation and/or progression include TET2, CBL, SH2B3, ASXL1, DNMT3A, IDH1/2, IKZF1, EZH2, SRSF2, and TP53. In order to identify novel somatic mutations associated with classic BCR-ABL1-negative MPNs, we performed whole genome sequencing of DNA extracted from peripheral blood granulocytes and cultured skin fibroblasts from a patient with PMF and a known MPL W515K mutation. Whole genome sequencing (WGS) was undertaken in a 55 year-old man with untreated PMF four years after initial diagnosis. His DIPSS Plus risk group was low (score 0). His karyotype was normal, and molecular testing revealed wild-type JAK2 in addition to the MPL W515K mutation. WGS of purified granulocytes and paired cultured skin fibroblasts was performed using both Illumina HiSeq and Complete Genomics (CGI) platforms. The resulting data were analyzed using multiple independent aligners and variant callers. Amplicon-based targeted resequencing with the Illumina MiSeq platform was used to evaluate additional patient samples for recurrent mutations. Stanford institutional review board approval and informed pt consent was obtained for these analyses. The PMF genome was sequenced to 88X (Illumina) and 128X (CGI) average fold coverage, and the cultured skin fibroblast genome was sequenced to 47X (Illumina) and 126X (CGI). The PMF genome had a low somatic mutation rate, consistent with that observed for other sequenced hematopoietic tumor genomes, with a low number of true somatic mutation calls. To definitively identify true mutations among various sequencing artefacts and germline variants, we use cultured skin fibroblasts which can be prepared with no contamination by neoplastic cells. In addition to re-identification of the MPL W515K mutation, this approach identified six additional somatic mutations that alter gene coding regions, splice sites, or known regulatory regions: a nonsense mutation in CARD6, implicated in modulation of NF-kappaB activation; a splice-sit