Richard de Borja

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

Hotspot activating PRKD1 somatic mutations in polymorphous low-grade adenocarcinomas of the salivary glands

Polymorphous low-grade adenocarcinoma (PLGA) is the second most frequent type of malignant tumor of the minor salivary glands. We identified PRKD1 hotspot mutations encoding p.Glu710Asp in 72.9% of PLGAs but not in other salivary gland tumors. Functional studies demonstrated that this kinase-activating alteration likely constitutes a driver of PLGA.

Genomic hallmarks of localized, non-indolent prostate cancer

Prostate tumours are highly variable in their response to therapies, but clinically available prognostic factors can explain only a fraction of this heterogeneity. Here we analysed 200 whole-genome sequences and 277 additional whole-exome sequences from localized, non-indolent prostate tumours with similar clinical risk profiles, and carried out RNA and methylation analyses in a subset. These tumours had a paucity of clinically actionable single nucleotide variants, unlike those seen in metastatic disease. Rather, a significant proportion of tumours harboured recurrent non-coding aberrations, large-scale genomic rearrangements, and alterations in which an inversion repressed transcription within its boundaries. Local hypermutation events were frequent, and correlated with specific genomic profiles. Numerous molecular aberrations were prognostic for disease recurrence, including several DNA methylation events, and a signature comprised of these aberrations outperformed well-described prognostic biomarkers. We suggest that intensified treatment of genomically aggressive localized prostate cancer may improve cure rates.

Robust global microRNA expression profiling using next-generation sequencing technologies

miRNAs are a class of regulatory molecules involved in a wide range of cellular functions, including growth, development and apoptosis. Given their widespread roles in biological processes, understanding their patterns of expression in normal and diseased states will provide insights into the consequences of aberrant expression. As such, global miRNA expression profiling of human malignancies is gaining popularity in both basic and clinically driven research. However, to date, the majority of such analyses have used microarrays and quantitative real-time PCR. With the introduction of digital count technologies, such as next-generation sequencing (NGS) and the NanoString nCounter System, we have at our disposal many more options. To make effective use of these different platforms, the strengths and pitfalls of several miRNA profiling technologies were assessed, including a microarray platform, NGS technologies and the NanoString nCounter System. Overall, NGS had the greatest detection sensitivity, largest dynamic range of detection and highest accuracy in differential expression analysis when compared with gold-standard quantitative real-time PCR. Its technical reproducibility was high, with intrasample correlations of at least 0.95 in all cases. Furthermore, miRNA analysis of formalin-fixed, paraffin-embedded (FFPE) tissue was also evaluated. Expression profiles between paired frozen and FFPE samples were similar, with Spearman’s ρ>0.93. These results show the superior sensitivity, accuracy and robustness of NGS for the comprehensive profiling of miRNAs in both frozen and FFPE tissues.

Identification of genes expressed by immune cells of the colon that are regulated by colorectal cancer-associated variants

A locus on human chromosome 11q23 tagged by marker rs3802842 was associated with colorectal cancer (CRC) in a genome‐wide association study; this finding has been replicated in case–control studies worldwide. In order to identify biologic factors at this locus that are related to the etiopathology of CRC, we used microarray‐based target selection methods, coupled to next‐generation sequencing, to study 103 kb at the 11q23 locus. We genotyped 369 putative variants from 1,030 patients with CRC (cases) and 1,061 individuals without CRC (controls) from the Ontario Familial Colorectal Cancer Registry. Two previously uncharacterized genes, COLCA1 and COLCA2, were found to be co‐regulated genes that are transcribed from opposite strands. Expression levels of COLCA1 and COLCA2 transcripts correlate with rs3802842 genotypes. In colon tissues, COLCA1 co‐localizes with crystalloid granules of eosinophils and granular organelles of mast cells, neutrophils, macrophages, dendritic cells and differentiated myeloid‐derived cell lines. COLCA2 is present in the cytoplasm of normal epithelial, immune and other cell lineages, as well as tumor cells. Tissue microarray analysis demonstrates the association of rs3802842 with lymphocyte density in the lamina propria (p = 0.014) and levels of COLCA1 in the lamina propria (p = 0.00016) and COLCA2 (tumor cells, p = 0.0041 and lamina propria, p = 6 × 10–5). In conclusion, genetic, expression and immunohistochemical data implicate COLCA1 and COLCA2 in the pathogenesis of colon cancer. Histologic analyses indicate the involvement of immune pathways.

BPG: Seamless, Automated and Interactive Visualization of Scientific Data

We introduce BPG, an easy-to-use framework for generating publication-quality, highly-customizable plots in the R statistical environment. This open-source package includes novel methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it ideal for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for seamless integration with computational pipelines. BPG is available at http://labs.oicr.on.ca/boutros-lab/software/bpg

Abstract B09: DNA polymerase mutations trigger rapid onset of ultra-hypermutant malignant brain tumors in children with biallelic mismatch repair deficiency

Background: Biallelic Mismatch Repair Deficiency (bMMRD) is a childhood cancer predisposition syndrome caused by germline mutations in MSH2, MSH6, MLH1, and PMS2. The leading cause of death is malignant brain tumors. The genomic landscape and secondary somatic mutations of bMMRD brain tumors are unknown. Methods: We analyzed 27 cancers and corresponding normal tissues from bMMRD patients using genome, exome sequencing and SNP-arrays. Additionally, we performed sequential sequencing from five primary and recurrent tumor pairs. Results: BMMRD malignant brain tumors harbored massive numbers of substitution mutations (>250/Mb), greater than all childhood and most adult cancers (>7,000 analyzed). These cancers lacked copy number alterations (p Conclusions/Significance: Early-onset brain tumors from bMMRD patients have a unique mechanism of malignant progression through secondary mutations in DNA polymerases. During transformation, brain tumors quickly reach a threshold of mutations developed in a rapid burst once a mutation in a DNA polymerase is acquired. The high mutation load and threshold of bMMRD cancers may be its Achilles9 heel, exploitable for diagnosis and therapeutic intervention. Note: This abstract was not presented at the conference. Citation Format: Adam Shlien, Brittany B. Campbell, Richard de Borja, Ludmil B. Alexandrov, Daniele Merico, David Wedge, Peter Van Loo, Patrick S. Tarpey, Paul Coupland, Aaron Pollett, Tatiana Lipman, Abolfazl Heidari, Shriya Deshmukh, Moritz Gerstung, Diana Merino, Manasa Ramakrishna, Marc Remke, Roland Arnold, Gagan B. Panigrahi, Samina Afzal, Valerie Larouche, Harriet Druker, Jordan Lerner-Ellis, Matthew Mistry, Rina Dvir, Ronald Grant, Ronit Elhasid, Roula Farah, Glenn P. Taylor, Paul C. Nathan, Sarah Alexander, Shay Ben-Shachar, Nada Jabado, Steven Gallinger, Shlohmi Constantini, Peter Dirks, Annie Huang, Steven W. Scherer, Richard G. Grundy, Carol Durno, Melyssa Aronson, M Stephen Meyn, Michael D. Taylor, Zachary F. Pursell, Christopher E. Pearson, David Malkin, P Andrew Futreal, Cynthia Hawkins, Eric Bouffet, Michael D. Taylor, Peter J. Campbell, Uri Tabori. DNA polymerase mutations trigger rapid onset of ultra-hypermutant malignant brain tumors in children with biallelic mismatch repair deficiency. [abstract]. In: Proceedings of the AACR Special Conference: Advances in Brain Cancer Research; May 27-30, 2015; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2015;75(23 Suppl):Abstract nr B09.

Abstract B129: Clinical implications of inter- and intra- prostatic heterogeneity.

Men with localized prostate cancer vary widely in clinical outcome, with a 30-50% failure rate after primary treatment. There is thus significant interest in developing genomically refined prognostic groups. We sought to evaluate the extent of genetic heterogeneity, both between patients (inter-prostate) and within individual prostate glands (intra-prostate) to assess the impact of genetic heterogeneity on risk stratification within a tight clinical cohort. Copy number aberrations (CNAs) from 75 Gleason 7 patients were determined by OncoScan SNP microarrays. We measure the percentage of genome involved in a CNA, termed percent genome aberration (PGA), a proxy for genomic instability. Additionally, whole genome sequencing was applied to 10 intermediate-risk prostate tumours and matched blood, including multiple manually macro-dissected regions from 5 of the prostates (range 2 to 9). Somatic single nucleotide variants (SNVs) and genomic rearrangements (GR) were extracted from each patient. We find a high degree of inter-prostatic heterogeneity between the 75 Gleason 7 patients, with the number of CNAs per patient ranging from 0 to 929, corresponding to PGA 0 to 16.7%. Known prognostic markers can differentiate between patients at higher risk for biochemical recurrence, but only account for a fraction of the cohort. Notably, when these prognostic genes are examined within multiple regions of five independent tumours, they differ in copy number between cancerous regions of the same prostate. For example, TP53 is deleted in 1/2, 1/3, 4/9, 0/4, and 4/5 prostate regions. Indeed, phylogenetic analysis of geographically distinct regions revealed multi-clonal disease in two of the five patients; separate analyses based on SNVs, CNAs, and GRs all concluded that these patient have two genetically distinct cancers within their prostate. We demonstrate dramatic levels of inter- and intra- prostate genetic heterogeneity within pathologically identical or similar cancers. The observed intra-prostatic genomic heterogeneity, both in terms of multi-focal and multi-clonal disease, has critical implications for clinical management. Prognostic information obtained by biopsy may be inconsistent depending on the site of biopsy, and applying personalized medicine to prostate cancer will be challenging. This study highlights the need for further evaluation of how intra-prostatic heterogeneity is related to patient prognosis. Citation Information: Mol Cancer Ther 2013;12(11 Suppl):B129. Citation Format: Emilie Lalonde, Paul C. Boutros, Michael Fraser, Richard de Borja, Nicholas J. Harding, Dominique Trudel, Alice Meng, Pablo H. Hennings-Yeomans, Andrew McPherson, Amin Zia, Jianxin Wang, Timothy Beck, Natalie S. Fox, Taryne Chong, Michelle Sam, Jeremy Johns, Lee Timms, Nicholas Buchner, Sohrab Shah, Cenk Sahinalp, Thomas J. Hudson, John D. McPherson, Theodorus van der Kwast, Robert G. Bristow. Clinical implications of inter- and intra- prostatic heterogeneity. [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2013 Oct 19-23; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2013;12(11 Suppl):Abstract nr B129.

Abstract 5276: Normalization of miRNA-sequencing data.

microRNAs (miRNA) are endogenous, small non-coding nucleotides that negatively regulate gene expression post-transcriptionally. Through interactions with Argonaute (Ago) proteins, they form the RNA-induced silencing complex (RISC) and can recognize and bind to the 3’UTR of mRNAs in a sequence-specific manner, leading to translational inhibition or mRNA degradation. Over 30% of human protein-coding genes are predicted to be conserved targets of miRNAs. Consequently, changes in their expression are likely to be associated with the development and progression of diseases, including cancer. An increasing number of studies are utilizing high-throughput sequencing over microarrays for the expression profiling of miRNAs. Processing of the raw sequencing data usually involves filtering based on quality measures, trimming for adapters and mapping the reads to a (genome or miRNA) reference. Then, normalization of the data is crucial before any downstream analysis can be performed. Normalization is the process of removing sources of variation, which are of non-biological origins (stemming from sample handling, library preparation, imaging and so on), and can affect the measured expression levels. An effective normalization method should minimize technical and experimental bias without introducing noise; the differences that remain should be truly biological effects. Several normalization methods for miRNA-seq data have been proposed, including linear scaling, non-linear scaling, quantile normalization and variance stabilization normalization. These methods differ in terms of complexity and the assumptions made. However, no standard technique has been recommended. Read counts from each experiment are usually simply adjusted for differences in sequencing depth (library size) to reads-per-million (RPM). Unfortunately, the performance and appropriateness of any of the normalization methods cannot be assessed using real data because the true values are not known. To this end, we have used a 12x12 Latin Square design to spike in 12 different oligonucleotides with known nominal concentrations, into a complex mixture of human miRNAs. These spike-in pools were subjected to all the preparatory steps of small RNA library construction for sequencing on the Illumina HiSeq2000. Preliminary results show that the spike-in sequences can be recovered successfully from the data. Using this data set, the relative merits of different normalization procedures are being assessed based on measures of bias, variance and improved sensitivity and specificity for the detection of differentially expressed miRNAs. The goal is to identify an optimal normalization method for miRNA-seq data, which would reduce variance without increasing bias. Citation Format: Shirley Tam, Richard de Borja, Ming-Sound Tsao, John D. McPherson. Normalization of miRNA-sequencing data. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5276. doi:10.1158/1538-7445.AM2013-5276

Abstract 2003: A molecular portrait of potentially curable prostate cancer.

Intermediate risk prostate cancer (CaP) with Gleason score (GS) of 7 show up to 100x variability in genetic instability. As CaP is multifocal and likely multiclonal, there is a need to characterize heterogeneity for patient stratification, which would increase the ability to act on genomic information by adding adjuvant therapies to offset systemic occult metastases that currently limit cure in ∼30% of patients. Individual genetic portraits could be used to improve cure on combined clinical-molecular staging criteria. We undertook a pilot study to assess the genetic heterogeneity of potentially curable GS=7 CaP. We selected 10 men with GS=7 CaP; 5 treated with external beam radiotherapy (frozen pre-treatment biopsies) and 5 treated with radical prostatectomy (RadP, frozen tumour). Additionally, DNA from 18 distinct formalin-fixed, paraffin-embedded (FFPE) foci from the 5 RadP were analysed. Each of these 28 foci were subjected to whole-genome sequencing (WGS) and OncoScan SNP arrays to yield comprehensive genetic profiles. mRNA expression was evaluated on frozen RadP by microarray. Germline DNA from whole-blood was also analysed. Following independent pathology reviews and manual macro-dissection of tumour areas of ≥70% cellularity, WGS (≥50x tumour, ≥30x germline) was performed on as little as 50 ng genomic DNA, and OncoScan arrays were performed using as little as 30ng DNA using either amplified or innate genomic DNA. Regions of CaP in FFPE RadP were recorded using a tissue map to identify independent malignant foci, and ERG immunostaining was performed to assist in the identification. In cases where ERG-positive and -negative foci were adjacent, ERG staining was repeated on an un-stained slide to confirm separate foci based on 3D multi-section analyses. ERG fusion status was also assessed in frozen samples by aCGH or IHC. Validation of SNVs via SNP array and deep-resequencing showed ∼99% accuracy. Tumour cellularity was estimated using Qpure and was >60% for all samples. Phylogenetic techniques were used to demonstrate clear multi-clonality in two tumours. Across all tumours, ∼50% of SNVs were specific to an individual tumour-region. Phylogenies were confirmed with both SNVs and CNAs, but CNAs generally exhibited greater concordance amongst different regions of the same tumour. Some previously observed recurrent mutations were previously identified as recurrent in CaP (e.g. SPOP), and the overall mutation rate for intermediate-risk CaP was only somewhat below that reported for castrate-resistant disease (11,230 somatic SNVs per tumour). Our studies support the concept that a complete characterization of inter- and intra-CaP heterogeneity is possible in fresh and archival tissues; the latter is important for correlations to clinical outcome. These approaches can then be streamlined for high-throughput analyses within personalized medicine laboratories leading to “point of care” molecular tests and individualization of therapy. Citation Format: Michael E. Fraser, Richard de Borja, Dominique Trudel, Nicholas J. Harding, Pablo H. Hennings-Yeomans, Alice Meng, Emilie R. Lalonde, Andrew Brown, Natalie S. Fox, Taryne Chong, Amin Zia, Michelle Sam, Jianxin Wang, Michelle A. Chan-Seng-Yue, Jeremy Johns, Lee Timms, Nicholas Buchner, Ada Wong, Fouad Yousif, Rob Denroche, Gaetano Zafarana, Maud HW Starmans, Hanbert Chen, Shaylan Govind, Francis Nguyen, Melania Pintilie, Neil Fleshner, Stanislav Volik, Lakshmi Muthuswamy, Colin C. Collins, Thomas J. Hudson, Lincoln D. Stein, Timothy Beck, John D. McPherson, Theodorus van der Kwast, Paul C. Boutros, Rob G. Bristow. A molecular portrait of potentially curable prostate cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2003. doi:10.1158/1538-7445.AM2013-2003