Matthew Ducar

Genomic sequencing of meningiomas identifies oncogenic SMO and AKT1 mutations

Meningiomas are the most common primary nervous system tumor. The tumor suppressor NF2 is disrupted in approximately half of all meningiomas, but the complete spectrum of genetic changes remains undefined. We performed whole-genome or whole-exome sequencing on 17 meningiomas and focused sequencing on an additional 48 tumors to identify and validate somatic genetic alterations. Most meningiomas had simple genomes, with fewer mutations, rearrangements and copy-number alterations than reported in other tumors in adults. However, several meningiomas harbored more complex patterns of copy-number changes and rearrangements, including one tumor with chromothripsis. We confirmed focal NF2 inactivation in 43% of tumors and found alterations in epigenetic modifiers in an additional 8% of tumors. A subset of meningiomas lacking NF2 alterations harbored recurrent oncogenic mutations in AKT1 (p.Glu17Lys) and SMO (p.Trp535Leu) and exhibited immunohistochemical evidence of activation of these pathways. These mutations were present in therapeutically challenging tumors of the skull base and higher grade. These results begin to define the spectrum of genetic alterations in meningiomas and identify potential therapeutic targets.

Integration of gene mutations in risk prognostication for patients receiving first-line immunochemotherapy for follicular lymphoma: a retrospective analysis of a prospective clinical trial and validation in a population-based registry

BACKGROUND Follicular lymphoma is a clinically and genetically heterogeneous disease, but the prognostic value of somatic mutations has not been systematically assessed. We aimed to improve risk stratification of patients receiving first-line immunochemotherapy by integrating gene mutations into a prognostic model. METHODS We did DNA deep sequencing to retrospectively analyse the mutation status of 74 genes in 151 follicular lymphoma biopsy specimens that were obtained from patients within 1 year before beginning immunochemotherapy consisting of rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP). These patients were recruited between May 4, 2000, and Oct 20, 2010, as part of a phase 3 trial (GLSG2000). Eligible patients had symptomatic, advanced stage follicular lymphoma and were previously untreated. The primary endpoints were failure-free survival (defined as less than a partial remission at the end of induction, relapse, progression, or death) and overall survival calculated from date of treatment initiation. Median follow-up was 7·7 years (IQR 5·5-9·3). Mutations and clinical factors were incorporated into a risk model for failure-free survival using multivariable L1-penalised Cox regression. We validated the risk model in an independent population-based cohort of 107 patients with symptomatic follicular lymphoma considered ineligible for curative irradiation. Pretreatment biopsies were taken between Feb 24, 2004, and Nov 24, 2009, within 1 year before beginning first-line immunochemotherapy consisting of rituximab, cyclophosphamide, vincristine, and prednisone (R-CVP). Median follow-up was 6·7 years (IQR 5·7-7·6). FINDINGS We established a clinicogenetic risk model (termed m7-FLIPI) that included the mutation status of seven genes (EZH2, ARID1A, MEF2B, EP300, FOXO1, CREBBP, and CARD11), the Follicular Lymphoma International Prognostic Index (FLIPI), and Eastern Cooperative Oncology Group (ECOG) performance status. In the training cohort, m7-FLIPI defined a high-risk group (28%, 43/151) with 5-year failure-free survival of 38·29% (95% CI 25·31-57·95) versus 77·21% (95% CI 69·21-86·14) for the low-risk group (hazard ratio [HR] 4·14, 95% CI 2·47-6·93; p<0·0001; bootstrap-corrected HR 2·02), and outperformed a prognostic model of only gene mutations (HR 3·76, 95% CI 2·10-6·74; p<0·0001; bootstrap-corrected HR 1·57). The positive predictive value and negative predictive value for 5-year failure-free survival were 64% and 78%, respectively, with a C-index of 0·80 (95% CI 0·71-0·89). In the validation cohort, m7-FLIPI again defined a high-risk group (22%, 24/107) with 5-year failure-free survival of 25·00% (95% CI 12·50-49·99) versus 68·24% (58·84-79·15) in the low-risk group (HR 3·58, 95% CI 2·00-6·42; p<0.0001). The positive predictive value for 5-year failure-free survival was 72% and 68% for negative predictive value, with a C-index of 0·79 (95% CI 0·69-0·89). In the validation cohort, risk stratification by m7-FLIPI outperformed FLIPI alone (HR 2·18, 95% CI 1·21-3·92), and FLIPI combined with ECOG performance status (HR 2·03, 95% CI 1·12-3·67). INTERPRETATION Integration of the mutational status of seven genes with clinical risk factors improves prognostication for patients with follicular lymphoma receiving first-line immunochemotherapy and is a promising approach to identify the subset at highest risk of treatment failure. FUNDING Deutsche Krebshilfe, Terry Fox Research Institute.

Genomic analysis of diffuse pediatric low-grade gliomas identifies recurrent oncogenic truncating rearrangements in the transcription factor MYBL1

Pediatric low-grade gliomas (PLGGs) are among the most common solid tumors in children but, apart from BRAF kinase mutations or duplications in specific subclasses, few genetic driver events are known. Diffuse PLGGs comprise a set of uncommon subtypes that exhibit invasive growth and are therefore especially challenging clinically. We performed high-resolution copy-number analysis on 44 formalin-fixed, paraffin-embedded diffuse PLGGs to identify recurrent alterations. Diffuse PLGGs exhibited fewer such alterations than adult low-grade gliomas, but we identified several significantly recurrent events. The most significant event, 8q13.1 gain, was observed in 28% of diffuse astrocytoma grade IIs and resulted in partial duplication of the transcription factor MYBL1 with truncation of its C-terminal negative-regulatory domain. A similar recurrent deletion-truncation breakpoint was identified in two angiocentric gliomas in the related gene v-myb avian myeloblastosis viral oncogene homolog (MYB) on 6q23.3. Whole-genome sequencing of a MYBL1-rearranged diffuse astrocytoma grade II demonstrated MYBL1 tandem duplication and few other events. Truncated MYBL1 transcripts identified in this tumor induced anchorage-independent growth in 3T3 cells and tumor formation in nude mice. Truncated transcripts were also expressed in two additional tumors with MYBL1 partial duplication. Our results define clinically relevant molecular subclasses of diffuse PLGGs and highlight a potential role for the MYB family in the biology of low-grade gliomas.

Targetable genetic features of primary testicular and primary central nervous system lymphomas

Primary central nervous system lymphomas (PCNSLs) and primary testicular lymphomas (PTLs) are extranodal large B-cell lymphomas (LBCLs) with inferior responses to current empiric treatment regimens. To identify targetable genetic features of PCNSL and PTL, we characterized their recurrent somatic mutations, chromosomal rearrangements, copy number alterations (CNAs), and associated driver genes, and compared these comprehensive genetic signatures to those of diffuse LBCL and primary mediastinal large B-cell lymphoma (PMBL). These studies identify unique combinations of genetic alterations in discrete LBCL subtypes and subtype-selective bases for targeted therapy. PCNSLs and PTLs frequently exhibit genomic instability, and near-uniform, often biallelic, CDKN2A loss with rare TP53 mutations. PCNSLs and PTLs also use multiple genetic mechanisms to target key genes and pathways and exhibit near-uniform oncogenic Toll-like receptor signaling as a result of MYD88 mutation and/or NFKBIZ amplification, frequent concurrent B-cell receptor pathway activation, and deregulation of BCL6. Of great interest, PCNSLs and PTLs also have frequent 9p24.1/PD-L1/PD-L2 CNAs and additional translocations of these loci, structural bases of immune evasion that are shared with PMBL.

Somatic activating ARAF mutations in Langerhans cell histiocytosis

The extracellular signal-regulated kinase (ERK) signaling pathway is activated in Langerhans cell histiocytosis (LCH) histiocytes, but only 60% of cases carry somatic activating mutations of BRAF. To identify other genetic causes of ERK pathway activation, we performed whole exome sequencing on purified LCH cells in 3 cases. One patient with wild-type BRAF alleles in his histiocytes had compound mutations in the kinase domain of ARAF. Unlike wild-type ARAF, this mutant was a highly active mitogen-activated protein kinase kinase in vitro and was capable of transforming mouse embryo fibroblasts. Mutant ARAF activity was inhibited by vemurafenib, a BRAF inhibitor, indicating the importance of fully evaluating ERK pathway abnormalities in selecting LCH patients for targeted inhibitor therapy.

MAP2K1 and MAP3K1 mutations in langerhans cell histiocytosis

Langerhans cell histiocytosis (LCH) is now understood to be a neoplastic disease in which over 50% of cases have somatic activating mutations of BRAF. However, the extracellular signal‐related (ERK) pathway is activated in all cases including those with wild type BRAF alleles. Here, we applied a targeted massively parallel sequencing panel to 30 LCH samples to test for the presence of additional genetic alterations that might cause ERK pathway activation. In 20 BRAF wild type samples, we found 3 somatic mutations in MAP2K1 (MEK1) including C121S and C121S/G128D in the kinase domain, and 56_61QKQKVG>R, an in‐frame deletion in the N‐terminal regulatory domain. All three variant proteins constitutively phosphorylated ERK in in vitro kinase assays. The C121S/G128D and 56_61QKQKVG>R variants were resistant to the mitogen‐activated protein kinase kinase (MEK) inhibitor trametinib in vitro. Within the entire sample set, we found 3 specimens with mutations in MAP3K1 (MEKK1), including two truncation mutants, T779fs and T1481fs; T1481fs encoded an unstable and nonfunctional protein when expressed in vitro. T779fs was present in a specimen carrying BRAF V600E. The third variant was a single nucleotide substitution, E1286V, which was fully functional and is likely a germline polymorphism. These results indicate that LCH cells can harbor additional genetic alterations in the RAS‐RAF‐MEK pathway which, in the case of MAP2K1, may be responsible for ERK activation in a wild type BRAF setting. The resistance of some of these variants to trametinib may also have clinical implications for the combined use of RAF and MEK inhibitors in LCH.

MYB-QKI rearrangements in angiocentric glioma drive tumorigenicity through a tripartite mechanism.

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs, including 19 angiocentric gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in angiocentric gliomas. In vitro and in vivo functional studies show that MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression and hemizygous loss of the tumor suppressor QKI. To our knowledge, this represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

BreaKmer: detection of structural variation in targeted massively parallel sequencing data using kmers

Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Prospective Enterprise-Level Molecular Genotyping of a Cohort of Cancer Patients

Ongoing cancer genome characterization studies continue to elucidate the spectrum of genomic abnormalities that drive many cancers, and in the clinical arena assessment of the driver genetic alterations in patients is playing an increasingly important diagnostic and/or prognostic role for many cancer types. However, the landscape of genomic abnormalities is still unknown for less common cancers, and the influence of specific genotypes on clinical behavior is often still unclear. To address some of these deficiencies, we developed Profile, a prospective cohort study to obtain genomic information on all patients at a large tertiary care medical center for cancer-related care. We enrolled patients with any cancer diagnosis, and, for each patient (unselected for cancer site or type) we applied mass spectrometric genotyping (OncoMap) of 471 common recurrent mutations in 41 cancer-related genes. We report the results of the first 5000 patients, of which 26% exhibited potentially actionable somatic mutations. These observations indicate the utility of genotyping in advancing the field of precision oncology.

The impact of tumor profiling approaches and genomic data strategies for cancer precision medicine

BackgroundThe diversity of clinical tumor profiling approaches (small panels to whole exomes with matched or unmatched germline analysis) may engender uncertainty about their benefits and liabilities, particularly in light of reported germline false positives in tumor-only profiling and use of global mutational and/or neoantigen data. The goal of this study was to determine the impact of genomic analysis strategies on error rates and data interpretation across contexts and ancestries.MethodsWe modeled common tumor profiling modalities—large (n = 300 genes), medium (n = 48 genes), and small (n = 15 genes) panels—using clinical whole exomes (WES) from 157 patients with lung or colon adenocarcinoma. We created a tumor-only analysis algorithm to assess germline false positive rates, the impact of patient ancestry on tumor-only results, and neoantigen detection.ResultsAfter optimizing a germline filtering strategy, the germline false positive rate with tumor-only large panel sequencing was 14 % (144/1012 variants). For patients whose tumor-only results underwent molecular pathologist review (n = 91), 50/54 (93 %) false positives were correctly interpreted as uncertain variants. Increased germline false positives were observed in tumor-only sequencing of non-European compared with European ancestry patients (p < 0.001; Fisher’s exact) when basic germline filtering approaches were used; however, the ExAC database (60,706 germline exomes) mitigated this disparity (p = 0.53). Matched and unmatched large panel mutational load correlated with WES mutational load (r2 = 0.99 and 0.93, respectively; p < 0.001). Neoantigen load also correlated (r2 = 0.80; p < 0.001), though WES identified a broader spectrum of neoantigens. Small panels did not predict mutational or neoantigen load.ConclusionsLarge tumor-only targeted panels are sufficient for most somatic variant identification and mutational load prediction if paired with expanded germline analysis strategies and molecular pathologist review. Paired germline sequencing reduced overall false positive mutation calls and WES provided the most neoantigens. Without patient-matched germline data, large germline databases are needed to minimize false positive mutation calling and mitigate ethnic disparities.