The impact of whole genome and transcriptome analysis (WGTA) on predictive biomarker discovery and diagnostic accuracy of advanced malignancies.

In this study, we evaluate the impact of whole genome and transcriptome analysis (WGTA) on predictive molecular profiling and histologic diagnosis in a cohort of advanced malignancies. WGTA was used to generate reports including molecular alterations and site/tissue of origin prediction. Two reviewers analyzed genomic reports, clinical history, and tumor pathology. We used National Comprehensive Cancer Network (NCCN) consensus guidelines, Food and Drug Administration (FDA) approvals, and provincially reimbursed treatments to define genomic biomarkers associated with approved targeted therapeutic options (TTOs). Tumor tissue/site of origin was reassessed for most cases using genomic analysis, including a machine learning algorithm (Supervised Cancer Origin Prediction Using Expression [SCOPE]) trained on The Cancer Genome Atlas data. WGTA was performed on 652 cases, including a range of primary tumor types/tumor sites and 15 malignant tumors of uncertain histogenesis (MTUH). At the time WGTA was performed, alterations associated with an approved TTO were identified in 39 (6%) cases; 3 of these were not identified through routine pathology workup. In seven (1%) cases, the pathology workup either failed, was not performed, or gave a different result from the WGTA. Approved TTOs identified by WGTA increased to 103 (16%) when applying 2021 guidelines. The histopathologic diagnosis was reviewed in 389 cases and agreed with the diagnostic consensus after WGTA in 94% of non-MTUH cases (n = 374). The remainder included situations where the morphologic diagnosis was changed based on WGTA and clinical data (0.5%), or where the WGTA was non-contributory (5%). The 15 MTUH were all diagnosed as specific tumor types by WGTA. Tumor board reviews including WGTA agreed with almost all initial predictive molecular profile and histopathologic diagnoses. WGTA was a powerful tool to assign site/tissue of origin in MTUH. Current efforts focus on improving therapeutic predictive power and decreasing cost to enhance use of WGTA data as a routine clinical test.

A platform for oncogenomic reporting and interpretation.

Manual interpretation of variants remains rate limiting in precision oncology. The increasing scale and complexity of molecular data generated from comprehensive sequencing of cancer samples requires advanced interpretative platforms as precision oncology expands beyond individual patients to entire populations. To address this unmet need, we introduce a Platform for Oncogenomic Reporting and Interpretation (PORI), comprising an analytic framework that facilitates the interpretation and reporting of somatic variants in cancer. PORI integrates reporting and graph knowledge base tools combined with support for manual curation at the reporting stage. PORI represents an open-source platform alternative to commercial reporting solutions suitable for comprehensive genomic data sets in precision oncology. We demonstrate the utility of PORI by matching 9,961 pan-cancer genome atlas tumours to the graph knowledge base, calculating therapeutically informative alterations, and making available reports describing select individual samples.

Therapeutic Potential of Afatinib in NRG1 Fusion-Driven Solid Tumors: A Case Series.

BACKGROUND: Neuregulin 1 (NRG1) fusions, which activate ErbB signaling, are rare oncogenic drivers in multiple tumor types. Afatinib is a pan-ErbB family inhibitor that may be an effective treatment for NRG1 fusion-driven tumors. PATIENTS AND METHODS: This report summarizes pertinent details, including best tumor response to treatment, for six patients with metastatic NRG1 fusion-positive tumors treated with afatinib. RESULTS: The six cases include four female and two male patients who ranged in age from 34 to 69 years. Five of the cases are patients with lung cancer, including two patients with invasive mucinous adenocarcinoma and three patients with nonmucinous adenocarcinoma. The sixth case is a patient with colorectal cancer. NRG1 fusion partners for the patients with lung cancer were either CD74 or SDC4. The patient with colorectal cancer harbored a novel POMK-NRG1 fusion and a KRAS mutation. Two patients received afatinib as first- or second-line therapy, three patients received the drug as third- to fifth-line therapy, and one patient received afatinib as fifteenth-line therapy. Best response with afatinib was stable disease in two patients (duration up to 16 months when combined with local therapies) and partial response (PR) of >18 months in three patients, including one with ongoing PR after 27 months. The remaining patient had a PR of 5 months with afatinib 40 mg/day, then another 6 months after an increase to 50 mg/day. CONCLUSION: This report reviews previously published metastatic NRG1 fusion-positive tumors treated with afatinib and summarizes six previously unpublished cases. The latter include several with a prolonged response to treatment (>18 months), as well as the first report of efficacy in NRG1 fusion-positive colorectal cancer. This adds to the growing body of evidence suggesting that afatinib can be effective in patients with NRG1 fusion-positive tumors. KEY POINTS: NRG1 fusions activate ErbB signaling and have been identified as oncogenic drivers in multiple solid tumor types. Afatinib is a pan-ErbB family inhibitor authorized for the treatment of advanced non-small cell lung cancer that may be effective in NRG1 fusion-driven tumors. This report summarizes six previously unpublished cases of NRG1 fusion-driven cancers treated with afatinib, including five with metastatic lung cancer and one with metastatic colorectal cancer. Several patients showed a prolonged response of >18 months with afatinib treatment. This case series adds to the evidence suggesting a potential role for afatinib in this area of unmet medical need.

Inhibition of Extracellular Enzymes Exposed to Cyanopeptides.

Harmful cyanobacterial blooms in freshwater ecosystems produce bioactive secondary metabolites including cyanopeptides that pose ecological and human health risks. Only adverse effects of one class of cyanopeptides, microcystins, have been studied extensively and have consequently been included in water quality assessments. Inhibition is a commonly observed effect for enzymes exposed to cyanopeptides and has mostly been investigated for human biologically relevant model enzymes. Here, we investigated the inhibition of ubiquitous aquatic enzymes by cyanobacterial metabolites. Hydrolytic enzymes are utilized in the metabolism of aquatic organisms and extracellularly by heterotrophic bacteria to obtain assimilable substrates. The ubiquitous occurrence of hydrolytic enzymes leads to the co-occurrence with cyanopeptides especially during cyanobacterial blooms. Bacterial leucine aminopeptidase and alkaline phosphatase were exposed to cyanopeptide extracts of different cyanobacterial strains ( Microcystis aeruginosa wild type and microcystin-free mutant, Planktothrix rubescens) and purified cyanopeptides. We observed inhibition of aminopeptidase and phosphatase upon exposure, especially to the apolar fractions of the cyanobacterial extracts. Exposure to the dominant cyanopeptides in these extracts confirmed that purified microcystins, aerucyclamide A and cyanopeptolin A inhibit the aminopeptidase in the low mg L-1 range while the phosphatase was less affected. Inhibition of aquatic enzymes can reduce the turnover of nutrients and carbon substrates and may also impair metabolic functions of grazing organisms.

Pan-cancer analysis of advanced patient tumors reveals interactions between therapy and genomic landscapes.

Advanced and metastatic tumors with complex treatment histories drive cancer mortality. Here we describe the POG570 cohort, a comprehensive whole-genome, transcriptome and clinical dataset, amenable for exploration of the impacts of therapies on genomic landscapes. Previous exposure to DNA-damaging chemotherapies and mutations affecting DNA repair genes, including POLQ and genes encoding Polζ, were associated with genome-wide, therapy-induced mutagenesis. Exposure to platinum therapies coincided with signatures SBS31 and DSB5 and, when combined with DNA synthesis inhibitors, signature SBS17b. Alterations in ESR1, EGFR, CTNNB1, FGFR1, VEGFA and DPYD were consistent with drug resistance and sensitivity. Recurrent noncoding events were found in regulatory region hotspots of genes including TERT, PLEKHS1, AP2A1 and ADGRG6. Mutation burden and immune signatures corresponded with overall survival and response to immunotherapy. Our data offer a rich resource for investigation of advanced cancers and interpretation of whole-genome and transcriptome sequencing in the context of a cancer clinic.

Fluorouracil sensitivity in a head and neck squamous cell carcinoma with a somatic DPYD structural variant.

Head and neck squamous cell carcinoma (HNSCC) is one of the most common cancers worldwide and represents a heterogeneous group of tumors, the majority of which are treated with a combination of surgery, radiation, and chemotherapy. Fluoropyrimidine (5-FU) and its oral prodrug, capecitabine, are commonly prescribed treatments for several solid tumor types including HNSCC. 5-FU-associated toxicity is observed in ∼30% of treated patients and is largely caused by germline polymorphisms in DPYD, which encodes dihydropyrimidine dehydrogenase, a key enzyme of 5-FU catabolism and deactivation. Although the association of germline DPYD alterations with toxicity is well-described, the potential contribution of somatic DPYD alterations to 5-FU sensitivity has not been explored. In a patient with metastatic HNSCC, in-depth genomic and transcriptomic integrative analysis on a biopsy from a metastatic neck lesion revealed alterations in genes that are associated with 5-FU uptake and metabolism. These included a novel somatic structural variant resulting in a partial deletion affecting DPYD, a variant of unknown significance affecting SLC29A1, and homozygous deletion of MTAP There was no evidence of deleterious germline polymorphisms that have been associated with 5-FU toxicity, indicating a potential vulnerability of the tumor to 5-FU therapy. The discovery of the novel DPYD variant led to the initiation of 5-FU treatment that resulted in a rapid response lasting 17 wk, with subsequent relapse due to unknown resistance mechanisms. This suggests that somatic alterations present in this tumor may serve as markers for tumor sensitivity to 5-FU, aiding in the selection of personalized treatment strategies.

Text-mining clinically relevant cancer biomarkers for curation into the CIViC database.

BACKGROUND: Precision oncology involves analysis of individual cancer samples to understand the genes and pathways involved in the development and progression of a cancer. To improve patient care, knowledge of diagnostic, prognostic, predisposing, and drug response markers is essential. Several knowledgebases have been created by different groups to collate evidence for these associations. These include the open-access Clinical Interpretation of Variants in Cancer (CIViC) knowledgebase. These databases rely on time-consuming manual curation from skilled experts who read and interpret the relevant biomedical literature. METHODS: To aid in this curation and provide the greatest coverage for these databases, particularly CIViC, we propose the use of text mining approaches to extract these clinically relevant biomarkers from all available published literature. To this end, a group of cancer genomics experts annotated sentences that discussed biomarkers with their clinical associations and achieved good inter-annotator agreement. We then used a supervised learning approach to construct the CIViCmine knowledgebase. RESULTS: We extracted 121,589 relevant sentences from PubMed abstracts and PubMed Central Open Access full-text papers. CIViCmine contains over 87,412 biomarkers associated with 8035 genes, 337 drugs, and 572 cancer types, representing 25,818 abstracts and 39,795 full-text publications. CONCLUSIONS: Through integration with CIVIC, we provide a prioritized list of curatable clinically relevant cancer biomarkers as well as a resource that is valuable to other knowledgebases and precision cancer analysts in general. All data is publically available and distributed with a Creative Commons Zero license. The CIViCmine knowledgebase is available at http://bionlp.bcgsc.ca/civicmine/.

NRG1 Gene Fusions Are Recurrent, Clinically Actionable Gene Rearrangements in KRAS Wild-Type Pancreatic Ductal Adenocarcinoma.

PURPOSE: Gene fusions involving neuregulin 1 (NRG1) have been noted in multiple cancer types and have potential therapeutic implications. Although varying results have been reported in other cancer types, the efficacy of the HER-family kinase inhibitor afatinib in the treatment of NRG1 fusion-positive pancreatic ductal adenocarcinoma is not fully understood. EXPERIMENTAL DESIGN: Forty-seven patients with pancreatic ductal adenocarcinoma received comprehensive whole-genome and transcriptome sequencing and analysis. Two patients with gene fusions involving NRG1 received afatinib treatment, with response measured by pretreatment and posttreatment PET/CT imaging. RESULTS: Three of 47 (6%) patients with advanced pancreatic ductal adenocarcinoma were identified as KRAS wild type by whole-genome sequencing. All KRAS wild-type tumors were positive for gene fusions involving the ERBB3 ligand NRG1. Two of 3 patients with NRG1 fusion-positive tumors were treated with afatinib and demonstrated a significant and rapid response while on therapy. CONCLUSIONS: This work adds to a growing body of evidence that NRG1 gene fusions are recurrent, therapeutically actionable genomic events in pancreatic cancers. Based on the clinical outcomes described here, patients with KRAS wild-type tumors harboring NRG1 gene fusions may benefit from treatment with afatinib.See related commentary by Aguirre, p. 4589.

CancerMine: a literature-mined resource for drivers, oncogenes and tumor suppressors in cancer.

Tumors from individuals with cancer are frequently genetically profiled for information about the driving forces behind the disease. We present the CancerMine resource, a text-mined and routinely updated database of drivers, oncogenes and tumor suppressors in different types of cancer. All data are available online ( http://bionlp.bcgsc.ca/cancermine ) and downloadable under a Creative Commons Zero license for ease of use.

Base excision repair deficiency signatures implicate germline and somatic MUTYH aberrations in pancreatic ductal adenocarcinoma and breast cancer oncogenesis.

We report a case of early-onset pancreatic ductal adenocarcinoma in a patient harboring biallelic MUTYH germline mutations, whose tumor featured somatic mutational signatures consistent with defective MUTYH-mediated base excision repair and the associated driver KRAS transversion mutation p.Gly12Cys. Analysis of an additional 730 advanced cancer cases (N = 731) was undertaken to determine whether the mutational signatures were also present in tumors from germline MUTYH heterozygote carriers or if instead the signatures were only seen in those with biallelic loss of function. We identified two patients with breast cancer each carrying a pathogenic germline MUTYH variant with a somatic MUTYH copy loss leading to the germline variant being homozygous in the tumor and demonstrating the same somatic signatures. Our results suggest that monoallelic inactivation of MUTYH is not sufficient for C:G>A:T transversion signatures previously linked to MUTYH deficiency to arise (N = 9), but that biallelic complete loss of MUTYH function can cause such signatures to arise even in tumors not classically seen in MUTYH-associated polyposis (N = 3). Although defective MUTYH is not the only determinant of these signatures, MUTYH germline variants may be present in a subset of patients with tumors demonstrating elevated somatic signatures possibly suggestive of MUTYH deficiency (e.g., COSMIC Signature 18, SigProfiler SBS18/SBS36, SignatureAnalyzer SBS18/SBS36).

Whole genome and whole transcriptome genomic profiling of a metastatic eccrine porocarcinoma.

Eccrine porocarcinomas (EPs) are rare malignant tumours of the intraepidermic sweat gland duct and most often arise from benign eccrine poromas. Some recurrent somatic genomic events have been identified in these malignancies, but very little is known about the complexity of their molecular pathophysiology. We describe the whole genome and whole transcriptome genomic profiling of a metastatic EP in a 66-year-old male patient with a previous history of localized porocarcinoma of the scalp. Whole genome and whole transcriptome genomic profiling was performed on the metastatic EP. Whole genome sequencing was performed on blood-derived DNA in order to allow a comparison between germline and somatic events. We found somatic copy losses of several tumour suppressor genes including APC, PTEN and CDKN2A, CDKN2B and CDKN1A. We identified a somatic hemizygous CDKN2A pathogenic splice site variant. De novo transcriptome assembly revealed abnormal splicing of CDKN2A p14ARF and p16INK4a. Elevated expression of oncogenes EGFR and NOTCH1 was noted and no somatic mutations were found in these genes. Wnt pathway somatic alterations were also observed. In conclusion, our results suggest that the molecular pathophysiology of malignant EP features high complexity and subtle interactions of multiple key genes. Cell cycle dysregulation and CDKN2A loss of function was found to be a new potential driver in EP tumourigenesis. Moreover, the combination of somatic copy number variants and abnormal gene expression perhaps partly related to epigenetic mechanisms, all likely contribute to the development of this rare malignancy in our patient.

Personalized oncogenomic analysis of metastatic adenoid cystic carcinoma: using whole-genome sequencing to inform clinical decision-making.

Metastatic adenoid cystic carcinomas (ACCs) can cause significant morbidity and mortality. Because of their slow growth and relative rarity, there is limited evidence for systemic therapy regimens. Recently, molecular profiling studies have begun to reveal the genetic landscape of these poorly understood cancers, and new treatment possibilities are beginning to emerge. The objective is to use whole-genome and transcriptome sequencing and analysis to better understand the genetic alterations underlying the pathology of metastatic and rare ACCs and determine potentially actionable therapeutic targets. We report five cases of metastatic ACC, not originating in the salivary glands, in patients enrolled in the Personalized Oncogenomics (POG) Program at the BC Cancer Agency. Genomic workup included whole-genome and transcriptome sequencing, detailed analysis of tumor alterations, and integration with existing knowledge of drug-target combinations to identify potential therapeutic targets. Analysis reveals low mutational burden in these five ACC cases, and mutation signatures that are commonly observed in multiple cancer types. Notably, the only recurrent structural aberration identified was the well-described MYB-NFIB fusion that was present in four of five cases, and one case exhibited a closely related MYBL1-NFIB fusion. Recurrent mutations were also identified in BAP1 and BCOR, with additional mutations in individual samples affecting NOTCH1 and the epigenetic regulators ARID2, SMARCA2, and SMARCB1. Copy changes were rare, and they included amplification of MYC and homozygous loss of CDKN2A in individual samples. Genomic analysis revealed therapeutic targets in all five cases and served to inform a therapeutic choice in three of the cases to date.

Assessment of Capture and Amplicon-Based Approaches for the Development of a Targeted Next-Generation Sequencing Pipeline to Personalize Lymphoma Management.

Targeted next-generation sequencing panels are increasingly used to assess the value of gene mutations for clinical diagnostic purposes. For assay development, amplicon-based methods have been preferentially used on the basis of short preparation time and small DNA input amounts. However, capture sequencing has emerged as an alternative approach because of high testing accuracy. We compared capture hybridization and amplicon sequencing approaches using fresh-frozen and formalin-fixed, paraffin-embedded tumor samples from eight lymphoma patients. Next, we developed a targeted sequencing pipeline using a 32-gene panel for accurate detection of actionable mutations in formalin-fixed, paraffin-embedded tumor samples of the most common lymphocytic malignancies: chronic lymphocytic leukemia, diffuse large B-cell lymphoma, and follicular lymphoma. We show that hybrid capture is superior to amplicon sequencing by providing deep more uniform coverage and yielding higher sensitivity for variant calling. Sanger sequencing of 588 variants identified specificity limits of thresholds for mutation calling, and orthogonal validation on 66 cases indicated 93% concordance with whole-genome sequencing. The developed pipeline and assay identified at least one actionable mutation in 91% of tumors from 219 lymphoma patients and revealed subtype-specific mutation patterns and frequencies consistent with the literature. This pipeline is an accurate and sensitive method for identifying actionable gene mutations in routinely acquired biopsy materials, suggesting further assessment of capture-based assays in the context of personalized lymphoma management.

Comparative RNA-Sequencing Analysis Benefits a Pediatric Patient With Relapsed Cancer.

Clinical detection of sequence and structural variants in known cancer genes points to viable treatment options for a minority of children with cancer.1 To increase the number of children who benefit from genomic profiling, gene expression information must be considered alongside mutations.2,3 Although high expression has been used to nominate drug targets for pediatric cancers,4,5 its utility has not been evaluated in a systematic way.6 We describe a child with a rare sarcoma that was profiled with whole-genome and RNA sequencing (RNA-Seq) techniques. Although the tumor did not harbor DNA mutations targetable by available therapies, incorporation of gene expression information derived from RNA-Seq analysis led to a therapy that produced a significant clinical response. We use this case to describe a framework for inclusion of gene expression into the clinical genomic evaluation of pediatric tumors.

A complete workflow for high-resolution spectral-stitching nanoelectrospray direct-infusion mass-spectrometry-based metabolomics and lipidomics.

Metabolomic and lipidomic studies measure and discover metabolic and lipid profiles in biological samples, enabling a better understanding of the metabolism of specific biological phenotypes. Accurate biological interpretations require high analytical reproducibility and sensitivity, and standardized and transparent data processing. Here we describe a complete workflow for nanoelectrospray ionization (nESI) direct-infusion mass spectrometry (DIMS) metabolomics and lipidomics. After metabolite and lipid extraction from tissues and biofluids, samples are directly infused into a high-resolution mass spectrometer (e.g., Orbitrap) using a chip-based nESI sample delivery system. nESI functions to minimize ionization suppression or enhancement effects as compared with standard electrospray ionization (ESI). Our analytical technique-named spectral stitching-measures data as several overlapping mass-to-charge (m/z) windows that are subsequently 'stitched' together, creating a complete mass spectrum. This considerably increases the dynamic range and detection sensitivity-about a fivefold increase in peak detection-as compared with the collection of DIMS data as a single wide mass-to-charge (m/z ratio) window. Data processing, statistical analysis and metabolite annotation are executed as a workflow within the user-friendly, transparent and freely available Galaxy platform (galaxyproject.org). Generated data have high mass accuracy that enables molecular formulae peak annotations. The workflow is compatible with any sample-extraction method; in this protocol, the examples are extracted using a biphasic method, with methanol, chloroform and water as the solvents. The complete workflow is reproducible, rapid and automated, which enables cost-effective analysis of >10,000 samples per year, making it ideal for high-throughput metabolomics and lipidomics screening-e.g., for clinical phenotyping, drug screening and toxicity testing.

The ortholog of the human proto-oncogene ROS1 is required for epithelial development in C. elegans.

The orphan receptor ROS1 is a human proto-oncogene, mutations of which are found in an increasing number of cancers. Little is known about the role of ROS1, however in vertebrates it has been implicated in promoting differentiation programs in specialized epithelial tissues. In this study we show that the C. elegans ortholog of ROS1, the receptor tyrosine kinase ROL-3, has an essential role in orchestrating the morphogenesis and development of specialized epidermal tissues, highlighting a potentially conserved function in coordinating crosstalk between developing epithelial cells. We also provide evidence of a direct relationship between ROL-3, the mucin SRAP-1, and BCC-1, the homolog of mRNA regulating protein Bicaudal-C. This study answers a longstanding question as to the developmental function of ROL-3, identifies three new genes that are expressed and function in the developing epithelium of C. elegans, and introduces the nematode as a potentially powerful model system for investigating the increasingly important, yet poorly understood, human oncogene ROS1.

The atm-1 gene is required for genome stability in Caenorhabditis elegans.

The Ataxia-telangiectasia-mutated (ATM) gene in humans was identified as the basis of a rare autosomal disorder leading to cancer susceptibility and is now well known as an important signal transducer in response to DNA damage. An approach to understanding the conserved functions of this gene is provided by the model system, Caenorhabditis elegans. In this paper we describe the structure and loss of function phenotype of the ortholog atm-1. Using bioinformatic and molecular analysis we show that the atm-1 gene was previously misannotated. We find that the transcript is in fact a product of three gene predictions, Y48G1BL.2 (atm-1), K10E9.1, and F56C11.4 that together make up the complete coding region of ATM-1. We also characterize animals that are mutant for two available knockout alleles, gk186 and tm5027. As expected, atm-1 mutant animals are sensitive to ionizing radiation. In addition, however, atm-1 mutants also display phenotypes associated with genomic instability, including low brood size, reduced viability and sterility. We document several chromosomal fusions arising from atm-1 mutant animals. This is the first time a mutator phenotype has been described for atm-1 in C. elegans. Finally we demonstrate the use of a balancer system to screen for and capture atm-1-derived mutational events. Our study establishes C. elegans as a model for the study of ATM as a mutator potentially leading to the development of screens to identify therapeutic targets in humans.