Comprehensive Analysis of Hypermutation in Human Cancer.

We present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

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.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

Myxoid smooth muscle neoplasia of the uterus: comprehensive analysis by next-generation sequencing and nucleic acid hybridization.

Uterine myxoid smooth muscle tumors, including myxoid leiomyosarcoma, are rare and their genomic profile has not been fully characterized. With the discovery of uterine sarcomas with ZC3H7B-BCOR fusion and BCOR internal tandem duplications, the differential diagnosis of myxoid smooth muscle lesions is expanding to include molecularly-defined tumors. Thus, we aimed to explore the genomic landscape of myxoid smooth muscle tumor using comprehensive tools. We performed whole exome next-generation sequencing and a pan-sarcoma RNA fusion assay in tumoral paraffin-embedded tissue from nine well-characterized uterine myxoid smooth muscle tumors (seven myxoid leiomyosarcomas and two myxoid smooth muscle tumors of unknown malignant potential). By immunohistochemistry, all tumors were strongly positive for smooth muscle markers and negative for BCOR staining; 4/6 expressed PLAG1. None of the tumors harbored known fusions including ZC3H7B-BCOR, TRPS1-PLAG1, and RAD51B-PLAG1. None harbored exon 15 BCOR internal tandem duplications; however, four tumors contained BCOR internal tandem duplications of unknown significance (mostly intronic). Mutational burden was low (median 3.8 mutations/megabase). DNA damage repair pathway gene mutations, including TP53 and BRCA2, were found. Copy number variation load, inferred from sequencing data, was variable with genomic indexes ranging from 2.2 to 74.7 (median 25.7), with higher indexes in myxoid leiomyosarcomas than myxoid smooth muscle tumors of unknown malignant potential. The absence of clear driver mutations suggests myxoid smooth muscle tumors to be genetically heterogeneous group of tumours and that other genetic (eg., undiscovered translocation) or epigenetic events drive the pathogenesis of uterine myxoid smooth muscle neoplasia.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

SeqControl: process control for DNA sequencing.

As high-throughput sequencing continues to increase in speed and throughput, routine clinical and industrial application draws closer. These 'production' settings will require enhanced quality monitoring and quality control to optimize output and reduce costs. We developed SeqControl, a framework for predicting sequencing quality and coverage using a set of 15 metrics describing overall coverage, coverage distribution, basewise coverage and basewise quality. Using whole-genome sequences of 27 prostate cancers and 26 normal references, we derived multivariate models that predict sequencing quality and depth. SeqControl robustly predicted how much sequencing was required to reach a given coverage depth (area under the curve (AUC) = 0.993), accurately classified clinically relevant formalin-fixed, paraffin-embedded samples, and made predictions from as little as one-eighth of a sequencing lane (AUC = 0.967). These techniques can be immediately incorporated into existing sequencing pipelines to monitor data quality in real time. SeqControl is available at http://labs.oicr.on.ca/Boutros-lab/software/SeqControl/.

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.

Coronavirus disease 2019 (COVID-19) outbreak on an in-patient medical unit associated with unrecognized exposures in common areas-Epidemiological and whole-genome sequencing investigation.

OBJECTIVE: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hospital outbreaks have been common and devastating during the coronavirus disease 2019 (COVID-19) pandemic. Understanding SARS-CoV-2 transmission in these environments is critical for preventing and managing outbreaks. DESIGN: Outbreak investigation through epidemiological mapping and whole-genome sequencing phylogeny. SETTING: Hospital in-patient medical unit outbreak in Toronto, Canada, from November 2020 to January 2021. PARTICIPANTS: The outbreak involved 8 patients and 10 staff and was associated with 3 patient deaths. RESULTS: Patients being cared for in geriatric chairs at the nursing station were at high risk for both acquiring and transmitting SARS-CoV-2 to other patients and staff. Furthermore, given the informal nature of these transmissions, they were not initially recognized, which led to further transmission and missing the opportunity for preventative COVID-19 therapies. CONCLUSIONS: During outbreak prevention and management, the risk of informal patient care settings, such as geriatric chairs, should be considered. During high-risk periods or during outbreaks, efforts should be made to care for patients in their rooms when possible.

Germline TP53 mutations undergo copy number gain years prior to tumor diagnosis.

Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 pathogenic variants. Here, we perform whole-genome sequence (WGS) analysis of tumors from 22 patients with TP53 germline pathogenic variants. We observe somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identify near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurs earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, reveal that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising years prior to tumor diagnosis.

Impact of neoadjuvant chemotherapy on somatic mutation status in high-grade serous ovarian carcinoma.

BACKGROUND: Patients treated with neoadjuvant chemotherapy (NACT) for advanced high-grade serous ovarian carcinoma (HGSC) have a higher rate and shorter time to platinum-resistant recurrence compared to patients treated with primary cytoreductive surgery (PCS) and adjuvant chemotherapy. The purpose of this study is to determine the impact of NACT on somatic mutation status in platinum-sensitive and resistant HGSC. Patients with advanced HGSC who had a documented response to platinum-based NACT, a banked blood sample, and a banked tumor sample before and after NACT were identified. Whole exome and/or targeted deep sequencing was performed in matched normal and pre/post-NACT tumor samples from 3 platinum-resistant and 2 platinum-sensitive patients to identify somatic non-synonymous mutations at each time point. RESULTS: When comparing exonic non-synonymous mutations in pre-NACT and post-NACT samples from the same patient, an average of 41% (1-68%) of genes were mutated at both time points. There were no trends detected in the mutational burden following exposure to NACT in platinum-resistant vs. platinum-sensitive cases. The majority of mutated genes were unique to each case. We identified several genes that were commonly mutated in pre-NACT samples specific to platinum-resistant (CSPG4, SLC35G5, TUBA3D) or sensitive (CYP2D6, NUTM1, DNAH5) cases. Four mutated genes emerged exclusively in the platinum-resistant cases (ADGRV1, MUC17, MUC20, PAK2) following NACT. CONCLUSIONS: Patients with advanced HGSC present with significant intra-tumor heterogeneity. NACT significantly impacts the somatic mutation status irrespective of the time to recurrence. The mutated genes detected in chemo-naive pre-NACT tumor samples from either resistant or sensitive cases could potentially have a role in the prediction of chemotherapy response in patients scheduled to receive NACT; larger studies are required to further validate these genes.

Emergence of a mutation in the nucleocapsid gene of SARS-CoV-2 interferes with PCR detection in Canada.

The emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) was met with rapid development of robust molecular-based detection assays. Many SARS-CoV-2 molecular tests target multiple genetic regions of the virus to maximize detection and protect against diagnostic escape. Despite the relatively moderate mutational rate of SARS-CoV-2, numerous mutations with known negative impact on diagnostic assays have been identified. In early 2021, we identified four samples positive for SARS-CoV-2 with a nucleocapsid (N) gene drop out on Cepheid Xpert® Xpress SARS-CoV-2 assay. Sequencing revealed a single common mutation in the N gene C29200T. Spatiotemporal analysis showed that the mutation was found in at least six different Canadian provinces from May 2020 until May 2021. Phylogenetic analysis showed that this mutation arose multiple times in Canadian samples and is present in six different variants of interest and of concern. The Cepheid testing platform is commonly used in Canada including in remote regions. As such, the existence of N gene mutation dropouts required further investigation. While commercial SARS-CoV-2 molecular detection assays have contributed immensely to the response effort, many vendors are reluctant to make primer/probe sequences publicly available. Proprietary primer/probe sequences create diagnostic 'blind spots' for global SARS-CoV-2 sequence monitoring and limits the ability to detect and track the presence and prevalence of diagnostic escape mutations. We hope that our industry partners will seriously consider making primer/probe sequences available, so that diagnostic escape mutants can be identified promptly and responded to appropriately to maintain diagnostic accuracy.

Detection of the novel SARS-CoV-2 European lineage B.1.177 in Ontario, Canada.

Background: Travel-related dissemination of SARS-CoV-2 continues to contribute to the global pandemic. A novel SARS-CoV-2 lineage (B.1.177) reportedly arose in Spain in the summer of 2020, with subsequent spread across Europe linked to travel by infected individuals. Surveillance and monitoring through the use of whole genome sequencing (WGS) offers insights into the global and local movement of pathogens such as SARS-CoV-2 and can detect introductions of novel variants. Methods: We analysed the genomes of SARS-CoV-2 sequenced for surveillance purposes from specimens received by Public Health Ontario (Sept 6 - Oct 10, 2020), collected from individuals in eastern Ontario, which comprised the study sample. Taxonomic lineages were identified using pangolin (v2.08) and phylogenetic analysis incorporated publicly available genomes covering the same time period as the study sample. Epidemiological data collected from laboratory requisitions and standard reportable disease case investigation was integrated into the analysis. Results: Genomic surveillance identified a COVID-19 case with SARS-CoV-2 lineage B.1.177 from an individual in eastern Ontario in late September, 2020. The individual had recently returned from Europe. Genomic analysis with publicly available data indicate the most closely related genomes to this specimen were from Southern Europe. Genomic surveillance did not identify further cases with this lineage. Conclusions: Genomic surveillance allowed for early detection of a novel SARS-CoV-2 lineage in Ontario which was deemed to be travel related. This type of genomic-based surveillance is a key tool to measure the effectiveness of public health measures such as mandatory self-isolation for returned travellers, aimed at preventing onward transmission of newly introduced lineages of SARS-CoV-2.

Diverse Oncogenic Fusions and Distinct Gene Expression Patterns Define the Genomic Landscape of Pediatric Papillary Thyroid Carcinoma.

Pediatric papillary thyroid carcinoma (PPTC) is clinically distinct from adult-onset disease. Although there are higher rates of metastasis and recurrence in PPTC, prognosis remains highly favorable. Molecular characterization of PPTC has been lacking. Historically, only 40% to 50% of childhood papillary thyroid carcinoma (PTC) were known to be driven by genomic variants common to adult PTC; oncogenic drivers in the remainder were unknown. This contrasts with approximately 90% of adult PTC driven by a discrete number of variants. In this study, 52 PPTCs underwent candidate gene testing, followed in a subset by whole-exome and transcriptome sequencing. Within these samples, candidate gene testing identified variants in 31 (60%) tumors, while exome and transcriptome sequencing identified oncogenic variants in 19 of 21 (90%) remaining tumors. The latter were enriched for oncogenic fusions, with 11 nonrecurrent fusion transcripts, including two previously undescribed fusions, STRN-RET and TG-PBF. Most fusions were associated with 3' receptor tyrosine kinase (RTK) moieties: RET, MET, ALK, and NTRK3. For advanced (distally metastatic) tumors, a driver variant was described in 91%. Gene expression analysis defined three clusters that demonstrated distinct expression of genes involved in thyroid differentiation and MAPK signaling. Among RET-CCDC6-driven tumors, gene expression in pediatric tumors was distinguishable from that in adults. Collectively, these results show that the genomic landscape of pediatric PTC is different from adult PTC. Moreover, they identify genomic drivers in 98% of PPTCs, predominantly oncogenic fusion transcripts involving RTKs, with a pronounced impact on gene expression. Notably, most advanced tumors were driven by a variant for which targeted systemic therapy exists. SIGNIFICANCE: This study highlights important distinctions between the genomes and transcriptomes of pediatric and adult papillary thyroid carcinoma, with implications for understanding the biology, diagnosis, and treatment of advanced disease in children.

CanDIG: Federated network across Canada for multi-omic and health data discovery and analysis.

We present the Canadian Distributed Infrastructure for Genomics (CanDIG) platform, which enables federated querying and analysis of human genomics and linked biomedical data. CanDIG leverages the standards and frameworks of the Global Alliance for Genomics and Health (GA4GH) and currently hosts data for five pan-Canadian projects. We describe CanDIG's key design decisions and features as a guide for other federated data systems.

Cellular and molecular architecture of hematopoietic stem cells and progenitors in genetic models of bone marrow failure.

Inherited bone marrow failure syndromes, such as Fanconi anemia (FA) and Shwachman-Diamond syndrome (SDS), feature progressive cytopenia and a risk of acute myeloid leukemia (AML). Using deep phenotypic analysis of early progenitors in FA/SDS bone marrow samples, we revealed selective survival of progenitors that phenotypically resembled granulocyte-monocyte progenitors (GMP). Whole-exome and targeted sequencing of GMP-like cells in leukemia-free patients revealed a higher mutation load than in healthy controls and molecular changes that are characteristic of AML: increased G>A/C>T variants, decreased A>G/T>C variants, increased trinucleotide mutations at Xp(C>T)pT, and decreased mutation rates at Xp(C>T)pG sites compared with other Xp(C>T)pX sites and enrichment for Cancer Signature 1 (X indicates any nucleotide). Potential preleukemic targets in the GMP-like cells from patients with FA/SDS included SYNE1, DST, HUWE1, LRP2, NOTCH2, and TP53. Serial analysis of GMPs from an SDS patient who progressed to leukemia revealed a gradual increase in mutational burden, enrichment of G>A/C>T signature, and emergence of new clones. Interestingly, the molecular signature of marrow cells from 2 FA/SDS patients with leukemia was similar to that of FA/SDS patients without transformation. The predicted founding clones in SDS-derived AML harbored mutations in several genes, including TP53, while in FA-derived AML the mutated genes included ARID1B and SFPQ. We describe an architectural change in the hematopoietic hierarchy of FA/SDS with remarkable preservation of GMP-like populations harboring unique mutation signatures. GMP-like cells might represent a cellular reservoir for clonal evolution.

Cancers from Novel Pole-Mutant Mouse Models Provide Insights into Polymerase-Mediated Hypermutagenesis and Immune Checkpoint Blockade.

POLE mutations are a major cause of hypermutant cancers, yet questions remain regarding mechanisms of tumorigenesis, genotype-phenotype correlation, and therapeutic considerations. In this study, we establish mouse models harboring cancer-associated POLE mutations P286R and S459F, which cause rapid albeit distinct time to cancer initiation in vivo, independent of their exonuclease activity. Mouse and human correlates enabled novel stratification of POLE mutations into three groups based on clinical phenotype and mutagenicity. Cancers driven by these mutations displayed striking resemblance to the human ultrahypermutation and specific signatures. Furthermore, Pole-driven cancers exhibited a continuous and stochastic mutagenesis mechanism, resulting in intertumoral and intratumoral heterogeneity. Checkpoint blockade did not prevent Pole lymphomas, but rather likely promoted lymphomagenesis as observed in humans. These observations provide insights into the carcinogenesis of POLE-driven tumors and valuable information for genetic counseling, surveillance, and immunotherapy for patients. SIGNIFICANCE: Two mouse models of polymerase exonuclease deficiency shed light on mechanisms of mutation accumulation and considerations for immunotherapy.See related commentary by Wisdom and Kirsch p. 5459.

Explosive mutation accumulation triggered by heterozygous human Pol ε proofreading-deficiency is driven by suppression of mismatch repair.

Tumors defective for DNA polymerase (Pol) ε proofreading have the highest tumor mutation burden identified. A major unanswered question is whether loss of Pol ε proofreading by itself is sufficient to drive this mutagenesis, or whether additional factors are necessary. To address this, we used a combination of next generation sequencing and in vitro biochemistry on human cell lines engineered to have defects in Pol ε proofreading and mismatch repair. Absent mismatch repair, monoallelic Pol ε proofreading deficiency caused a rapid increase in a unique mutation signature, similar to that observed in tumors from patients with biallelic mismatch repair deficiency and heterozygous Pol ε mutations. Restoring mismatch repair was sufficient to suppress the explosive mutation accumulation. These results strongly suggest that concomitant suppression of mismatch repair, a hallmark of colorectal and other aggressive cancers, is a critical force for driving the explosive mutagenesis seen in tumors expressing exonuclease-deficient Pol ε.

Rearrangement bursts generate canonical gene fusions in bone and soft tissue tumors.

Sarcomas are cancers of the bone and soft tissue often defined by gene fusions. Ewing sarcoma involves fusions between EWSR1, a gene encoding an RNA binding protein, and E26 transformation-specific (ETS) transcription factors. We explored how and when EWSR1-ETS fusions arise by studying the whole genomes of Ewing sarcomas. In 52 of 124 (42%) of tumors, the fusion gene arises by a sudden burst of complex, loop-like rearrangements, a process called chromoplexy, rather than by simple reciprocal translocations. These loops always contained the disease-defining fusion at the center, but they disrupted multiple additional genes. The loops occurred preferentially in early replicating and transcriptionally active genomic regions. Similar loops forming canonical fusions were found in three other sarcoma types. Chromoplexy-generated fusions appear to be associated with an aggressive form of Ewing sarcoma. These loops arise early, giving rise to both primary and relapse Ewing sarcoma tumors, which can continue to evolve in parallel.