Lee Trani

Modernizing Reference Genome Assemblies

I have read the journals policy and have the following conflicts: Paul Flicek is married to the deputy editor of PLoS Medicine, Melissa Norton. Evan Eichler is on the board of Pacific Biosciences. Support for this work came from the Intramural Research Program of the NIH, The National Library of Medicine, the European Molecular Biology Laboratory, the Wellcome Trust (grant number 077198), and the Howard Hughes Medical Institute (EEE). The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.

CIViC is a community knowledgebase for expert crowdsourcing the clinical interpretation of variants in cancer

CIViC is an expert-crowdsourced knowledgebase for Clinical Interpretation of Variants in Cancer describing the therapeutic, prognostic, diagnostic and predisposing relevance of inherited and somatic variants of all types. CIViC is committed to open-source code, open-access content, public application programming interfaces (APIs) and provenance of supporting evidence to allow for the transparent creation of current and accurate variant interpretations for use in cancer precision medicine.

A Phase I Trial of BKM120 (Buparlisib) in Combination with Fulvestrant in Postmenopausal Women with Estrogen Receptor-Positive Metastatic Breast Cancer.

Purpose: This trial was conducted to determine the maximum tolerated dose (MTD) and preliminary efficacy of buparlisib, an oral pan-class I PI3K inhibitor, plus fulvestrant in postmenopausal women with metastatic estrogen receptor positive (ER+) breast cancer. Experimental Design: Phase IA employed a 3+3 design to determine the MTD of buparlisib daily plus fulvestrant. Subsequent cohorts (phase IB and cohort C) evaluated intermittent (5/7-day) and continuous dosing of buparlisib (100 mg daily). No more than 3 prior systemic treatments in the metastatic setting were allowed in these subsequent cohorts. Results: Thirty-one patients were enrolled. MTD was defined as buparlisib 100 mg daily plus fulvestrant. Common adverse events (AE) included fatigue (38.7%), transaminases elevation (35.5%), rash (29%), and diarrhea (19.4%). C-peptide was significantly increased during treatment, consistent with on-target effect of buparlisib. Compared with intermittent dosing, daily buparlisib was associated with more frequent early onset AEs and higher buparlisib plasma concentrations. Among the 29 evaluable patients, the clinical benefit rate was 58.6% (95% CI, 40.7%–74.5%). Response was not associated with PIK3CA mutation or treatment cohort; however, loss of PTEN, progesterone receptor (PgR) expression, or mutation in TP53 was most common in resistant cases, and mutations in AKT1 and ESR1 did not exclude treatment response. Conclusions: Buparlisib plus fulvestrant is clinically active with manageable AEs in patients with metastatic ER+ breast cancer. Weekend breaks in buparlisib dosing reduced toxicity. Patients with PgR negative and TP53 mutation did poorly, suggesting buparlisib plus fulvestrant may not be adequately effective against tumors with these poor prognostic molecular features. Clin Cancer Res; 22(7); 1583–91. ©2015 AACR.

Recurrent somatic mutations affecting B-cell receptor signaling pathway genes in follicular lymphoma

Follicular lymphoma (FL) is the most common form of indolent non-Hodgkin lymphoma, yet it remains only partially characterized at the genomic level. To improve our understanding of the genetic underpinnings of this incurable and clinically heterogeneous disease, whole-exome sequencing was performed on tumor/normal pairs from a discovery cohort of 24 patients with FL. Using these data and mutations identified in other B-cell malignancies, 1716 genes were sequenced in 113 FL tumor samples from 105 primarily treatment-naive individuals. We identified 39 genes that were mutated significantly above background mutation rates. CREBBP mutations were associated with inferior PFS. In contrast, mutations in previously unreported HVCN1, a voltage-gated proton channel-encoding gene and B-cell receptor signaling modulator, were associated with improved PFS. In total, 47 (44.8%) patients harbor mutations in the interconnected B-cell receptor (BCR) and CXCR4 signaling pathways. Histone gene mutations were more frequent than previously reported (identified in 43.8% of patients) and often co-occurred (17.1% of patients). A novel, recurrent hotspot was identified at a posttranslationally modified residue in the histone H2B family. This study expands the number of mutated genes described in several known signaling pathways and complexes involved in lymphoma pathogenesis (BCR, Notch, SWitch/sucrose nonfermentable (SWI/SNF), vacuolar ATPases) and identified novel recurrent mutations (EGR1/2, POU2AF1, BTK, ZNF608, HVCN1) that require further investigation in the context of FL biology, prognosis, and treatment.

Melorheostosis: Exome sequencing of an associated dermatosis implicates postzygotic mosaicism of mutated KRAS☆☆☆

Melorheostosis (MEL) is the rare sporadic dysostosis characterized by monostotic or polyostotic osteosclerosis and hyperostosis often distributed in a sclerotomal pattern. The prevailing hypothesis for MEL invokes postzygotic mosaicism. Sometimes scleroderma-like skin changes, considered a representation of the pathogenetic process of MEL, overlie the bony changes, and sometimes MEL becomes malignant. Osteopoikilosis (OPK) is the autosomal dominant skeletal dysplasia that features symmetrically distributed punctate osteosclerosis due to heterozygous loss-of-function mutation within LEMD3. Rarely, radiographic findings of MEL occur in OPK. However, germline mutation of LEMD3 does not explain sporadic MEL. To explore if mosaicism underlies MEL, we studied a boy with polyostotic MEL and characteristic overlying scleroderma-like skin, a few bony lesions consistent with OPK, and a large epidermal nevus known to usually harbor a HRAS, FGFR3, or PIK3CA gene mutation. Exome sequencing was performed to ~100× average read depth for his two dermatoses, two areas of normal skin, and peripheral blood leukocytes. As expected for non-malignant tissues, the patients mutation burden in his normal skin and leukocytes was low. He, his mother, and his maternal grandfather carried a heterozygous, germline, in-frame, 24-base-pair deletion in LEMD3. Radiographs of the patient and his mother revealed bony foci consistent with OPK, but she showed no MEL. For the patient, somatic variant analysis, using four algorithms to compare all 20 possible pairwise combinations of his five DNA samples, identified only one high-confidence mutation, heterozygous KRAS Q61H (NM_033360.3:c.183A>C, NP_203524.1:p.Gln61His), in both his dermatoses but absent in his normal skin and blood. Thus, sparing our patient biopsy of his MEL bone, we identified a heterozygous somatic KRAS mutation in his scleroderma-like dermatosis considered a surrogate for MEL. This implicates postzygotic mosaicism of mutated KRAS, perhaps facilitated by germline LEMD3 haploinsufficiency, causing his MEL.

CIViC: A knowledgebase for expert-crowdsourcing the clinical interpretation of variants in cancer.

CIViC is an expert crowdsourced knowledgebase for Clinical Interpretation of Variants in Cancer (www.civicdb.org) describing the therapeutic, prognostic, and diagnostic relevance of inherited and somatic variants of all types. CIViC is committed to open source code, open access content, public application programming interfaces (APIs), and provenance of supporting evidence to allow for the transparent creation of current and accurate variant interpretations for use in cancer precision medicine.

Genomic characterization of HER2-positive breast cancer and response to neoadjuvant trastuzumab and chemotherapy-results from the ACOSOG Z1041 (Alliance) trial

Background HER2 (ERBB2) gene amplification and its corresponding overexpression are present in 15-30% of invasive breast cancers. While HER2-targeted agents are effective treatments, resistance remains a major cause of death. The American College of Surgeons Oncology Group Z1041 trial (NCT00513292) was designed to compare the pathologic complete response (pCR) rate of distinct regimens of neoadjuvant chemotherapy and trastuzumab, but ultimately identified no difference. Patients and methods In supplement to tissues from 37 Z1041 cases, 11 similarly treated cases were obtained from a single institution study (NCT00353483). We have extracted genomic DNA from both pre-treatment tumor biopsies and blood of these 48 cases, and performed whole genome (WGS) and exome sequencing. Coincident with these efforts, we have generated RNA-seq profiles from 42 of the tumor biopsies. Among patients in this cohort, 24 (50%) achieved a pCR. Results We have characterized the genomic landscape of HER2-positive breast cancer and investigated associations between genomic features and pCR. Cases assigned to the HER2-enriched subtype by RNA-seq analysis were more likely to achieve a pCR compared to the luminal, basal-like, or normal-like subtypes (19/27 versus 3/15; P = 0.0032). Mutational events led to the generation of putatively active neoantigens, but were overall not associated with pCR. ERBB2 and GRB7 were the genes most commonly observed in fusion events, and genomic copy number analysis of the ERBB2 locus indicated that cases with either no observable or low-level ERBB2 amplification were less likely to achieve a pCR (7/8 versus 17/40; P = 0.048). Moreover, among cases that achieved a pCR, tumors consistently expressed immune signatures that may contribute to therapeutic response. Conclusion The identification of these features suggests that it may be possible to predict, at the time of diagnosis, those HER2-positive breast cancer patients who will not respond to treatment with chemotherapy and trastuzumab. ClinicalTrials.gov identifiers NCT00513292, NCT00353483.Background HER2 ( ERBB2 ) gene amplification and its corresponding overexpression are present in 15-30% of invasive breast cancers. While HER2-targeted agents are effective treatments, resistance remains a major cause of death. The American College of Surgeons Oncology Group Z1041 trial (NCT00513292) was designed to compare the pathologic complete response (pCR) rate of distinct regimens of neoadjuvant chemotherapy and trastuzumab, but ultimately identified no difference [1]. Patients and methods In supplement to tissues from 37 Z1041 cases, 11 similarly treated cases were obtained from a single institution study (NCT00353483). We have extracted genomic DNA from both pre-treatment tumor biopsies and blood of these 48 cases, and performed whole genome (WGS) and exome sequencing. Coincident with these efforts, we have generated RNA-seq profiles from 42 of the tumor biopsies. Among patients in this cohort, 24 (50%) achieved a pCR. Results Here we characterized the genomic landscape of HER2-positive breast cancer and investigated associations between genomic features and pCR. Cases assigned to the HER2-enriched subtype by RNA-seq analysis were more likely to achieve a pCR compared to the luminal, basal-like, or normal-like subtypes (19/27 versus 3/15; P =0.0032). Mutational events led to the generation of putatively active neoantigens, but were overall not associated with pCR. ERBB2 and GRB7 were the genes most commonly observed in fusion events, and genomic copy number analysis of the ERBB2 locus indicated that cases with either no observable or low-level ERBB2 amplification were less likely to achieve a pCR (7/8 versus 17/40; P =0.048). Moreover, among cases that achieved a pCR, tumors consistently expressed immune signatures that may contribute to therapeutic response. Conclusion The identification of these features suggests that it may be possible to predict, at the time of diagnosis, those HER2-positive breast cancer patients who will not respond to treatment with chemotherapy and trastuzumab.

Standard operating procedure for somatic variant refinement of tumor sequencing data

Purpose: Manual review of aligned sequencing reads is required to develop a high-quality list of somatic variants from massively parallel sequencing data (MPS). Despite widespread use in analyzing MPS data, there has been little attempt to describe methods for manual review, resulting in high inter- and intra-lab variability in somatic variant detection and characterization of tumors. Methods: Open source software was used to develop an optimal method for manual review setup. We also developed a systemic approach to visually inspect each variant during manual review. Results: We present a standard operating procedures for somatic variant refinement for use by manual reviewers. The approach is enhanced through representative examples of 4 different manual review categories that indicate a reviewer’s confidence in the somatic variant call and 19 annotation tags that contextualize commonly observed sequencing patterns during manual review. Representative examples provide detailed instructions on how to classify variants during manual review to rectify lack of confidence in automated somatic variant detection. Conclusion: Standardization of somatic variant refinement through systematization of manual review will improve the consistency and reproducibility of identifying true somatic variants after automated variant calling.

Recurrent WNT pathway alterations are frequent in relapsed small cell lung cancer

Nearly all patients with small cell lung cancer (SCLC) eventually relapse with chemoresistant disease. The molecular mechanisms driving chemoresistance in SCLC remain un-characterized. Here, we describe whole-exome sequencing of paired SCLC tumor samples procured at diagnosis and relapse from 12 patients, and unpaired relapse samples from 18 additional patients. Multiple somatic copy number alterations, including gains in ABCC1 and deletions in MYCL, MSH2, and MSH6, are identifiable in relapsed samples. Relapse samples also exhibit recurrent mutations and loss of heterozygosity in regulators of WNT signaling, including CHD8 and APC. Analysis of RNA-sequencing data shows enrichment for an ASCL1-low expression subtype and WNT activation in relapse samples. Activation of WNT signaling in chemosensitive human SCLC cell lines through APC knockdown induces chemoresistance. Additionally, in vitro-derived chemoresistant cell lines demonstrate increased WNT activity. Overall, our results suggest WNT signaling activation as a mechanism of chemoresistance in relapsed SCLC.Small cell lung cancer (SCLC) patients frequently relapse and become resistant to chemotherapy. Here, the authors analyse the genomic and transcriptomic landscape of primary and relapsed SCLC patients as well as in vitro models, and discover that activation of WNT signalling can drive chemotherapy resistance.

Standard operating procedure for somatic variant refinement of sequencing data with paired tumor and normal samples

PurposeFollowing automated variant calling, manual review of aligned read sequences is required to identify a high-quality list of somatic variants. Despite widespread use in analyzing sequence data, methods to standardize manual review have not been described, resulting in high inter- and intralab variability.MethodsThis manual review standard operating procedure (SOP) consists of methods to annotate variants with four different calls and 19 tags. The calls indicate a reviewer’s confidence in each variant and the tags indicate commonly observed sequencing patterns and artifacts that inform the manual review call. Four individuals were asked to classify variants prior to, and after, reading the SOP and accuracy was assessed by comparing reviewer calls with orthogonal validation sequencing.ResultsAfter reading the SOP, average accuracy in somatic variant identification increased by 16.7% (p value = 0.0298) and average interreviewer agreement increased by 12.7% (p value < 0.001). Manual review conducted after reading the SOP did not significantly increase reviewer time.ConclusionThis SOP supports and enhances manual somatic variant detection by improving reviewer accuracy while reducing the interreviewer variability for variant calling and annotation.