Homologous Recombination Deficiency and Platinum-Based Therapy Outcomes in Advanced Breast Cancer.

Purpose: Recent studies have identified mutation signatures of homologous recombination deficiency (HRD) in over 20% of breast cancers, as well as pancreatic, ovarian, and gastric cancers. There is an urgent need to understand the clinical implications of HRD signatures. Whereas BRCA1/2 mutations confer sensitivity to platinum-based chemotherapies, it is not yet clear whether mutation signatures can independently predict platinum response.Experimental Design: In this observational study, we sequenced tumor whole genomes (100× depth) and matched normals (60×) of 93 advanced-stage breast cancers (33 platinum-treated). We computed a published metric called HRDetect, independently trained to predict BRCA1/2 status, and assessed its capacity to predict outcomes on platinum-based chemotherapies. Clinical endpoints were overall survival (OS), total duration on platinum-based therapy (TDT), and radiographic evidence of clinical improvement (CI).Results: HRDetect predicted BRCA1/2 status with an area under the curve (AUC) of 0.94 and optimal threshold of 0.7. Elevated HRDetect was also significantly associated with CI on platinum-based therapy (AUC = 0.89; P = 0.006) with the same optimal threshold, even after adjusting for BRCA1/2 mutation status and treatment timing. HRDetect scores over 0.7 were associated with a 3-month extended median TDT (P = 0.0003) and 1.3-year extended median OS (P = 0.04).Conclusions: Our findings not only independently validate HRDetect, but also provide the first evidence of its association with platinum response in advanced breast cancer. We demonstrate that HRD mutation signatures may offer clinically relevant information independently of BRCA1/2 mutation status and hope this work will guide the development of clinical trials. Clin Cancer Res; 23(24); 7521-30. ©2017 AACR.

Spatial heterogeneity in medulloblastoma.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

Molecular etiology of an indolent lymphoproliferative disorder determined by whole-genome sequencing.

In an attempt to assess potential treatment options, whole-genome and transcriptome sequencing were performed on a patient with an unclassifiable small lymphoproliferative disorder. Variants from genome sequencing were prioritized using a combination of comparative variant distributions in a spectrum of lymphomas, and meta-analyses of gene expression profiling. In this patient, the molecular variants that we believe to be most relevant to the disease presentation most strongly resemble a diffuse large B-cell lymphoma (DLBCL), whereas the gene expression data are most consistent with a low-grade chronic lymphocytic leukemia (CLL). The variant of greatest interest was a predicted NOTCH2-truncating mutation, which has been recently reported in various lymphomas.

Putative BRAF activating fusion in a medullary thyroid cancer.

Medullary thyroid cancer (MTC) is a malignancy of the calcitonin-producing parafollicular cells of the thyroid gland. Surgery is the only curative treatment for this cancer. External beam radiation therapy is reserved for adjuvant treatment of MTC with aggressive features. Targeted therapeutics vandetanib and cabozantinib are approved for the treatment of aggressive and metastatic tumors that are not amenable to surgery. The use of these multikinase inhibitors are supported by the observed overactivation of the RET oncoprotein in a large subpopulation of MTCs. However, not all patients carry oncogenic alterations of this kinase. Hence, there is still a need for comprehensive molecular characterization of MTC utilizing whole-genome and transcriptome-sequencing methodologies with the aim of identifying targetable mutations. Here, we describe the genomic profiles of two medullary thyroid cancers and report the presence of a putative oncogenic BRAF fusion in one. Such alterations, previously observed in other malignancies and known targets of available drugs, can benefit patients who currently have no treatment options.

The genomic and transcriptomic landscape of anaplastic thyroid cancer: implications for therapy.

BACKGROUND: Anaplastic thyroid carcinoma is the most undifferentiated form of thyroid cancer and one of the deadliest of all adult solid malignancies. Here we report the first genomic and transcriptomic profile of anaplastic thyroid cancer including those of several unique cell lines and outline novel potential drivers of malignancy and targets of therapy. METHODS: We describe whole genomic and transcriptomic profiles of 1 primary anaplastic thyroid tumor and 3 authenticated cell lines. Those profiles augmented by the transcriptomes of 4 additional and unique cell lines were compared to 58 pairs of papillary thyroid carcinoma and matched normal tissue transcriptomes from The Cancer Genome Atlas study. RESULTS: The most prevalent mutations were those of TP53 and BRAF; repeated alterations of the epigenetic machinery such as frame-shift deletions of HDAC10 and EP300, loss of SMARCA2 and fusions of MECP2, BCL11A and SS18 were observed. Sequence data displayed aneuploidy and large regions of copy loss and gain in all genomes. Common regions of gain were however evident encompassing chromosomes 5p and 20q. We found novel anaplastic gene fusions including MKRN1-BRAF, FGFR2-OGDH and SS18-SLC5A11, all expressed in-frame fusions involving a known proto-oncogene. Comparison of the anaplastic thyroid cancer expression datasets with the papillary thyroid cancer and normal thyroid tissue transcriptomes suggested several known drug targets such as FGFRs, VEGFRs, KIT and RET to have lower expression levels in anaplastic specimens compared with both papillary thyroid cancers and normal tissues, confirming the observed lack of response to therapies targeting these pathways. Further integrative data analysis identified the mTOR signaling pathway as a potential therapeutic target in this disease. CONCLUSIONS: Anaplastic thyroid carcinoma possessed heterogeneous and unique profiles revealing the significance of detailed molecular profiling of individual tumors and the treatment of each as a unique entity; the cell line sequence data promises to facilitate the more accurate and intentional drug screening studies for anaplastic thyroid cancer.

MEN1 mutations in Hürthle cell (oncocytic) thyroid carcinoma.

CONTEXT AND OBJECTIVE: Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation, and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management. DESIGN AND PATIENTS: Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line, and five Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in the multiple endocrine neoplasia I (MEN1) gene. RESULTS: Here we report the identification of MEN1 loss-of-function mutations in 4% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumors. CONCLUSION: Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of MEN1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis, and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

Diagnostic value of next-generation sequencing in an unusual sphenoid tumor.

Extraordinary advancements in sequencing technology have made what was once a decade-long multi-institutional endeavor into a methodology with the potential for practical use in a clinical setting. We therefore set out to examine the clinical value of next-generation sequencing by enrolling patients with incurable or ambiguous tumors into the Personalized OncoGenomics initiative at the British Columbia Cancer Agency whereby whole genome and transcriptome analyses of tumor/normal tissue pairs are completed with the ultimate goal of directing therapeutics. First, we established that the sequencing, analysis, and communication with oncologists could be completed in less than 5 weeks. Second, we found that cancer diagnostics is an area that can greatly benefit from the comprehensiveness of a whole genome analysis. Here, we present a scenario in which a metastasized sphenoid mass, which was initially thought of as an undifferentiated squamous cell carcinoma, was rediagnosed as an SMARCB1-negative rhabdoid tumor based on the newly acquired finding of homozygous SMARCB1 deletion. The new diagnosis led to a change in chemotherapy and a complete nodal response in the patient. This study also provides additional insight into the mutational landscape of an adult SMARCB1-negative tumor that has not been explored at a whole genome and transcriptome level.

De novo transcriptome assembly with ABySS.

MOTIVATION: Whole transcriptome shotgun sequencing data from non-normalized samples offer unique opportunities to study the metabolic states of organisms. One can deduce gene expression levels using sequence coverage as a surrogate, identify coding changes or discover novel isoforms or transcripts. Especially for discovery of novel events, de novo assembly of transcriptomes is desirable. RESULTS: Transcriptome from tumor tissue of a patient with follicular lymphoma was sequenced with 36 base pair (bp) single- and paired-end reads on the Illumina Genome Analyzer II platform. We assembled approximately 194 million reads using ABySS into 66 921 contigs 100 bp or longer, with a maximum contig length of 10 951 bp, representing over 30 million base pairs of unique transcriptome sequence, or roughly 1% of the genome. AVAILABILITY AND IMPLEMENTATION: Source code and binaries of ABySS are freely available for download at http://www.bcgsc.ca/platform/bioinfo/software/abyss. Assembler tool is implemented in C++. The parallel version uses Open MPI. ABySS-Explorer tool is implemented in Java using the Java universal network/graph framework. CONTACT: ibirol@bcgsc.ca.

Novel expressed sequences identified in a model of androgen independent prostate cancer.

BACKGROUND: Prostate cancer is the most frequently diagnosed cancer in American men, and few effective treatment options are available to patients who develop hormone-refractory prostate cancer. The molecular changes that occur to allow prostate cells to proliferate in the absence of androgens are not fully understood. RESULTS: Subtractive hybridization experiments performed with samples from an in vivo model of hormonal progression identified 25 expressed sequences representing novel human transcripts. Intriguingly, these 25 sequences have small open-reading frames and are not highly conserved through evolution, suggesting many of these novel expressed sequences may be derived from untranslated regions of novel transcripts or from non-coding transcripts. Examination of a large metalibrary of human Serial Analysis of Gene Expression (SAGE) tags demonstrated that only three of these novel sequences had been previously detected. RT-PCR experiments confirmed that the 6 sequences tested were expressed in specific human tissues, as well as in clinical samples of prostate cancer. Further RT-PCR experiments for five of these fragments indicated they originated from large untranslated regions of unannotated transcripts. CONCLUSION: This study underlines the value of using complementary techniques in the annotation of the human genome. The tissue-specific expression of 4 of the 6 clones tested indicates the expression of these novel transcripts is tightly regulated, and future work will determine the possible role(s) these novel transcripts may play in the progression of prostate cancer.

Software for automated analysis of DNA fingerprinting gels.

Here we describe software tools for the automated detection of DNA restriction fragments resolved on agarose fingerprinting gels. We present a mathematical model for the location and shape of the restriction fragments as a function of fragment size, with model parameters determined empirically from "marker" lanes containing molecular size standards. Automated identification of restriction fragments involves several steps, including: image preprocessing, to put the data in a form consistent with a linear model; marker lane analysis, for determination of the model parameters; and data lane analysis, a procedure for detecting restriction fragment multiplets while simultaneously determining the amplitude curve that describes restriction fragment amplitude as a function of mobility. In validation experiments conducted on fingerprinted and sequenced Bacterial Artificial Chromosome (BAC) clones, sensitivity and specificity of restriction fragment identification exceeded 96% on restriction fragments ranging in size from 600 base pairs (bp) to 30,000 bp. The integrated suite of software tools, written in MATLAB and collectively called BandLeader, is in use at the BC Cancer Agency Genome Sciences Centre (GSC) and the Washington University Genome Sequencing Center, and has been provided to the Wellcome Trust Sanger Institute and the Whitehead Institute. Employed in a production mode at the GSC, BandLeader has been used to perform automated restriction fragment identification for more than 850,000 BAC clones for mouse, rat, bovine, and poplar fingerprint mapping projects.

The Genome sequence of the SARS-associated coronavirus.

We sequenced the 29,751-base genome of the severe acute respiratory syndrome (SARS)-associated coronavirus known as the Tor2 isolate. The genome sequence reveals that this coronavirus is only moderately related to other known coronaviruses, including two human coronaviruses, HCoV-OC43 and HCoV-229E. Phylogenetic analysis of the predicted viral proteins indicates that the virus does not closely resemble any of the three previously known groups of coronaviruses. The genome sequence will aid in the diagnosis of SARS virus infection in humans and potential animal hosts (using polymerase chain reaction and immunological tests), in the development of antivirals (including neutralizing antibodies), and in the identification of putative epitopes for vaccine development.