Lethal and sublethal effects towards zebrafish larvae of microcystins and other cyanopeptides produced by cyanobacteria

Cyanobacterial blooms affect aquatic ecosystems across the globe and one major concern relates to their toxins such as microcystins (MC). Yet, the ecotoxicological risks, particularly non-lethal effects, associated with other co-produced secondary metabolites remain mostly unknown. Here, we assessed survival, morphological alterations, swimming behaviour and cardiovascular functions of zebrafish (Danio rerio) upon exposure to cyanobacterial extracts of two Brazilian Microcystis strains. We verified that only MIRS-04 produced MCs and identified other co-produced cyanopeptides also for the MC non-producer NPCD-01 by LC-HRMS/MS analysis. Both cyanobacterial extracts, from the MC-producer and non-producer, caused acute toxicity in zebrafish with LC50 values of 0.49 and 0.98 mgdw_biomass/mL, respectively. After exposure to MC-producer extract, additional decreased locomotor activity was observed. The cyanopeptolin (micropeptin K139) contributed 52% of the overall mortality and caused oedemas of the pericardial region. Oedemas of the pericardial area and prevented hatching were also observed upon exposure to the fraction with high abundance of a microginin (Nostoginin BN741) in the extract of the MC non-producer. Our results further add to the yet sparse understanding of lethal and sublethal effects caused by cyanobacterial metabolites other than MCs and the need to better understand the underlying mechanisms of the toxicity. We emphasize the importance of considering mixture toxicity of co-produced metabolites in the ecotoxicological risk assessment of cyanobacterial bloom events, given the importance for predicting adverse outcomes in fish and other organisms.

The Precision Toxicology Initiative

The goal of PrecisionTox is to overcome conceptual barriers to replacing traditional mammalian chemical safety testing by accelerating the discovery of evolutionarily conserved toxicity pathways that are shared by descent among humans and more distantly related animals. An international consortium is systematically testing the toxicological effects of a diverse set of chemicals on a suite of five model species comprising fruit flies, nematodes, water fleas, and embryos of clawed frogs and zebrafish along with human cell lines. Multiple forms of omics and comparative toxicology data are integrated to map the evolutionary origins of biomolecular interactions, which are predictive of adverse health effects, to major branches of the animal phylogeny. These conserved elements of adverse outcome pathways (AOPs) and their biomarkers are expect to provide mechanistic insight useful for regulating groups of chemicals based on their shared modes of action. PrecisionTox also aims to quantify risk variation within populations by recognizing susceptibility as a heritable trait that varies with genetic diversity. This initiative incorporates legal experts and collaborates with risk managers to address specific needs within European chemicals legislation, including the uptake of new approach methodologies (NAMs) for setting precise regulatory limits on toxic chemicals.

Immune Activation following Irbesartan Treatment in a Colorectal Cancer Patient: A Case Study

Colorectal cancers are one of the most prevalent tumour types worldwide and, despite the emergence of targeted and biologic therapies, have among the highest mortality rates. The Personalized OncoGenomics (POG) program at BC Cancer performs whole genome and transcriptome analysis (WGTA) to identify specific alterations in an individual's cancer that may be most effectively targeted. Informed using WGTA, a patient with advanced mismatch repair-deficient colorectal cancer was treated with the antihypertensive drug irbesartan and experienced a profound and durable response. We describe the subsequent relapse of this patient and potential mechanisms of response using WGTA and multiplex immunohistochemistry (m-IHC) profiling of biopsies before and after treatment from the same metastatic site of the L3 spine. We did not observe marked differences in the genomic landscape before and after treatment. Analyses revealed an increase in immune signalling and infiltrating immune cells, particularly CD8+ T cells, in the relapsed tumour. These results indicate that the observed anti-tumour response to irbesartan may have been due to an activated immune response. Determining whether there may be other cancer contexts in which irbesartan may be similarly valuable will require additional studies.

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.

Whole genome and transcriptome analysis enhances precision cancer treatment options

BACKGROUND

Recent advances are enabling delivery of precision genomic medicine to cancer clinics. While the majority of approaches profile panels of selected genes or hotspot regions, comprehensive data provided by whole genome and transcriptome sequencing and analysis (WGTA) presents an opportunity to align a much larger proportion of patients to therapies.

PATIENTS AND METHODS

Samples from 570 patients with advanced or metastatic cancer of diverse types enrolled in the Personalized OncoGenomics (POG) program underwent WGTA. DNA-based data, including mutations, copy number, and mutation signatures, were combined with RNA-based data, including gene expression and fusions, to generate comprehensive WGTA profiles. A multidisciplinary molecular tumour board used WGTA profiles to identify and prioritize clinically actionable alterations and inform therapy. Patient responses to WGTA-informed therapies were collected.

RESULTS

Clinically actionable targets were identified for 83% of patients, 37% of whom received WGTA-informed treatments. RNA expression data were particularly informative, contributing to 67% of WGTA-informed treatments; 25% of treatments were informed by RNA expression alone. Of a total 248 WGTA-informed treatments, 46% resulted in clinical benefit. RNA expression data were comparable to DNA-based mutation and copy number data in aligning to clinically beneficial treatments. Genome signatures also guided therapeutics including platinum, PARP inhibitors, and immunotherapies. Patients accessed WGTA-informed treatments through clinical trials (19%), off-label use (35%), and as standard therapies (46%) including those which would not otherwise have been the next choice of therapy, demonstrating the utility of genomic information to direct use of chemotherapies as well as targeted therapies.

CONCLUSIONS

Integrating RNA expression and genome data illuminated treatment options that resulted in 46% of treated patients experiencing positive clinical benefit, supporting the use of comprehensive WGTA profiling in clinical cancer care.

CLINICAL TRIAL NUMBER

NCT02155621.

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.

The interpretation of somatic variants in cancer is challenging due to the scale and complexity of sequencing data. Here, the authors present PORI, an open-source framework for interpreting somatic variants in cancer using graph knowledge base tools, automated reporting, and manual curation.

Antagonistic Effects of Point Mutations on Charge Recombination and a New View of Primary Charge Separation in Photosynthetic Proteins

Light-induced electron-transfer reactions were investigated in wild-type and three mutant Rhodobacter sphaeroides reaction centers with the secondary electron acceptor (ubiquinone Q_A) either removed or permanently reduced. Under such conditions, charge separation between the primary electron donor (bacteriochlorophyll dimer, P) and the electron acceptor (bacteriopheophytin, H_A) was followed by P⁺H_A^– → PH_A charge recombination. Two reaction centers were used that had different single amino-acid mutations that brought about either a 3-fold acceleration in charge recombination compared to that in the wild-type protein, or a 3-fold deceleration. In a third mutant in which the two single amino-acid mutations were combined, charge recombination was similar to that in the wild type. In all cases, data from transient absorption measurements were analyzed using similar models. The modeling included the energetic relaxation of the charge-separated states caused by protein dynamics and evidenced the appearance of an intermediate charge-separated state, P⁺B_A^–, with B_A being the bacteriochlorophyll located between P and H_A. In all cases, mixing of the states P⁺B_A^– and P⁺H_A^– was observed and explained in terms of electron delocalization over B_A and H_A. This delocalization, together with picosecond protein relaxation, underlies a new view of primary charge separation in photosynthesis.

CyanoMetDB, a comprehensive public database of secondary metabolites from cyanobacteria

Harmful cyanobacterial blooms, which frequently contain toxic secondary metabolites, are reported in aquatic environments around the world. More than two thousand cyanobacterial secondary metabolites have been reported from diverse sources over the past fifty years. A comprehensive, publically-accessible database detailing these secondary metabolites would facilitate research into their occurrence, functions and toxicological risks. To address this need we created CyanoMetDB, a highly curated, flat-file, openly-accessible database of cyanobacterial secondary metabolites collated from 850 peer-reviewed articles published between 1967 and 2020. CyanoMetDB contains 2010 cyanobacterial metabolites and 99 structurally related compounds. This has nearly doubled the number of entries with complete literature metadata and structural composition information compared to previously available open access databases. The dataset includes microcytsins, cyanopeptolins, other depsipeptides, anabaenopeptins, microginins, aeruginosins, cyclamides, cryptophycins, saxitoxins, spumigins, microviridins, and anatoxins among other metabolite classes. A comprehensive database dedicated to cyanobacterial secondary metabolites facilitates: (1) the detection and dereplication of known cyanobacterial toxins and secondary metabolites; (2) the identification of novel natural products from cyanobacteria; (3) research on biosynthesis of cyanobacterial secondary metabolites, including substructure searches; and (4) the investigation of their abundance, persistence, and toxicity in natural environments.

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.

NRG1 fusion-driven tumors: biology, detection, and the therapeutic role of afatinib and other ErbB-targeting agents

Oncogenic gene fusions are hybrid genes that result from structural DNA rearrangements, leading to deregulated activity. Fusions involving the neuregulin-1 gene (NRG1) result in ErbB-mediated pathway activation and therefore present a rational candidate for targeted treatment. The most frequently reported NRG1 fusion is CD74-NRG1, which most commonly occurs in patients with invasive mucinous adenocarcinomas (IMAs) of the lung, although several other NRG1 fusion partners have been identified in patients with lung cancer, including ATP1B1, SDC4, and RBPMS. NRG1 fusions are also present in patients with other solid tumors, such as pancreatic ductal adenocarcinoma. In general, NRG1 fusions are rare across different types of cancer, with a reported incidence of <1%, with the notable exception of IMA, which represents ≈2%-10% of lung adenocarcinomas and has a reported incidence of ≈10%-30% for NRG1 fusions. A substantial proportion (≈20%) of NRG1 fusion-positive non-small-cell lung cancer cases are nonmucinous adenocarcinomas. ErbB-targeted treatments, such as afatinib, a pan-ErbB tyrosine kinase inhibitor, are potential therapeutic strategies to address unmet treatment needs in patients harboring NRG1 fusions.

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.

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.

iPSC modeling of young-onset Parkinson's disease reveals a molecular signature of disease and novel therapeutic candidates

Young-onset Parkinson's disease (YOPD), defined by onset at <50 years, accounts for approximately 10% of all Parkinson's disease cases and, while some cases are associated with known genetic mutations, most are not. Here induced pluripotent stem cells were generated from control individuals and from patients with YOPD with no known mutations. Following differentiation into cultures containing dopamine neurons, induced pluripotent stem cells from patients with YOPD showed increased accumulation of soluble α-synuclein protein and phosphorylated protein kinase Cα, as well as reduced abundance of lysosomal membrane proteins such as LAMP1. Testing activators of lysosomal function showed that specific phorbol esters, such as PEP005, reduced α-synuclein and phosphorylated protein kinase Cα levels while increasing LAMP1 abundance. Interestingly, the reduction in α-synuclein occurred through proteasomal degradation. PEP005 delivery to mouse striatum also decreased α-synuclein production in vivo. Induced pluripotent stem cell-derived dopaminergic cultures reveal a signature in patients with YOPD who have no known Parkinson's disease-related mutations, suggesting that there might be other genetic contributions to this disorder. This signature was normalized by specific phorbol esters, making them promising therapeutic candidates.

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/.

Afatinib as a novel potential treatment option for NRG1 fusion-positive tumors

110

Background: Neuregulin 1 (NRG1) is a growth factor that binds HER3/4 and activates ErbB signalling pathways. NRG1 gene fusions function as oncogenic drivers and represent a potential therapeutic target across tumour types. Afatinib, an ErbB-family blocker, is a potential treatment option for some patients with NRG1 fusion-positive ( NRG1+) tumours, as supported by preclinical evidence and clinical case reports. Of 13 patients with NRG1+ lung adenocarcinoma (ADC; n = 9), pancreatic ADC (n = 2), cholangiocarcinoma (n = 1), and ovarian cancer (n = 1) treated with afatinib, eight patients achieved a best response of partial response (PR; median duration 7.3 months, range 3–12), three patients had stable disease (SD) and two patients had progressive disease. Methods: Here, we report the clinico-pathological and molecular characteristics and current status of four additional cases of afatinib-treated NRG1+ tumours. Results: Patient 1, 66 year-old female, never-smoker with metastatic non-mucinous lung ADC. CD74-NRG1 fusion was identified, and 5th-line afatinib initiated. PR is ongoing after 16 months. Patient 2, 43 year-old female, non-smoker with advanced invasive mucinous lung ADC. CD74-NRG1 fusion was identified and 3rd-line afatinib initiated (PR, 18 months); treatment is ongoing following local progression. Patient 3, 69 year-old male, with KRAS-mutated metastatic colorectal cancer. Following a right hemicolectomy and liver/lung metastasectomies, a novel POMK-NRG1 fusion was detected and afatinib initiated (SD, 4 months). Eight months after initiation, afatinib treatment is ongoing, in combination with radiotherapy, with SD. Patient 4, 54 year-old male, with KRAS-wild-type metastatic pancreatic cancer. Following progression on chemotherapy, APP-NRG1 fusion was detected as part of the Personalized Oncogenomics study (NCT02155621), and afatinib initiated; PR is ongoing after 7 months. Conclusions: These findings add to a growing body of evidence that afatinib is a potential treatment option for patients with NRG1+ tumours. Prospective study is ongoing/planned in the Drug Rediscovery Protocol trial (DRUP; NCT02925234) and Targeted Agent and Profiling Utilization Registry study (TAPUR; NCT02693535).

Comprehensive genomic analysis of metastatic pancreatic ductal adenocarcinoma (mPDAC) reveals a significant proportion of clinical actionable aberrations

e15753

Background: Significant progress has been made in the understanding of the genomic landscape of PDAC, but the clinical utility of these data remains uncertain. Methods: As part of the BC Cancer Personalized Oncogenomics (POG) and PanGen studies (NCT02155621, NCT02869802), whole genome analysis and transcriptome sequencing were performed on fresh biopsy and blood sample from 48 mPDAC patients. Genomic findings informed therapy choices including potential eligibility for the CCTG PM.1 molecular basket trial (NCT03297606). Results: Cohort consists of 54.1% male, average age 57.6, 34/48 had ≥2 lines of treatment. 37/48 biopsies were from liver and 27/48 were collected pre-treatment. 8/48 (16.6%) patients had aberrations with strong evidence of clinical actionability. These include 2 germline BRCA2, 1 germline BRCA1, 1 somatic XRCC2 homozygous deletion with strong COSMIC signature 3, all predictive of platinum sensitivity. Patients with XRCC2 deletion and BRCA1 had over 2 years on FOLFIRINOX. Fusions affecting the NRG1 gene were identified in 3/4 patients with KRAS wildtype tumours, which may confer ERRB inhibitor sensitivity. 2/3 NRG1 fusion patients have thus far been treated with the ERBB inhibitor afatinib with radiographic responses noted in both patients. One patient had mismatch repair deficiency, and a high mutational burden, suggestive of immune checkpoint inhibitor sensitivity. Other possible actionable mutations include CCTG PM.1 trial potential eligibility: 4/48 with high homologous recombination defects and 1 germline ATM mutation loss (PARP inhibitor arm), 1/47 ERBB2 amplification (anti-HER2 arm), 2/47 high expression of FGFR1 (Sunitinib arm) and 1/47 with high FLT4 and IGF1R expression (Axitinib arm). Conclusions: More routine use of comprehensive genomic analysis should be considered in mPDAC given finding of high degree of actionability. Importantly, a significant proportion (16.6%) had findings with strong evidence of clinical impact.

Initial results from the New Horizons exploration of 2014 MU69, a small Kuiper Belt object

The Kuiper Belt is a distant region of the outer Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU₆₉, a cold classical Kuiper Belt object approximately 30 kilometers in diameter. Such objects have never been substantially heated by the Sun and are therefore well preserved since their formation. We describe initial results from these encounter observations. MU₆₉ is a bilobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color or compositional heterogeneity. No evidence for satellites, rings or other dust structures, a gas coma, or solar wind interactions was detected. MU₆₉'s origin appears consistent with pebble cloud collapse followed by a low-velocity merger of its two lobes.

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.

Curcumin restores innate immune Alzheimer's disease risk gene expression to ameliorate Alzheimer pathogenesis

Alzheimer's disease (AD) genetics implies a causal role for innate immune genes, TREM2 and CD33, products that oppose each other in the downstream Syk tyrosine kinase pathway, activating microglial phagocytosis of amyloid (Aβ). We report effects of low (Curc-lo) and high (Curc-hi) doses of curcumin on neuroinflammation in APPsw transgenic mice. Results showed that Curc-lo decreased CD33 and increased TREM2 expression (predicted to decrease AD risk) and also increased TyroBP, which controls a neuroinflammatory gene network implicated in AD as well as phagocytosis markers CD68 and Arg1. Curc-lo coordinately restored tightly correlated relationships between these genes' expression levels, and decreased expression of genes characteristic of toxic pro-inflammatory M1 microglia (CD11b, iNOS, COX-2, IL1β). In contrast, very high dose curcumin did not show these effects, failed to clear amyloid plaques, and dysregulated gene expression relationships. Curc-lo stimulated microglial migration to and phagocytosis of amyloid plaques both in vivo and in ex vivo assays of sections of human AD brain and of mouse brain. Curcumin also reduced levels of miR-155, a micro-RNA reported to drive a neurodegenerative microglial phenotype. In conditions without amyloid (human microglial cells in vitro, aged wild-type mice), Curc-lo similarly decreased CD33 and increased TREM2. Like curcumin, anti-Aβ antibody (also reported to engage the Syk pathway, increase CD68, and decrease amyloid burden in human and mouse brain) increased TREM2 in APPsw mice and decreased amyloid in human AD sections ex vivo. We conclude that curcumin is an immunomodulatory treatment capable of emulating anti-Aβ vaccine in stimulating phagocytic clearance of amyloid by reducing CD33 and increasing TREM2 and TyroBP, while restoring neuroinflammatory networks implicated in neurodegenerative diseases.

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).

Vitamin B6 catabolism and lung cancer risk: results from the Lung Cancer Cohort Consortium (LC3)

BACKGROUND

Increased vitamin B6 catabolism related to inflammation, as measured by the PAr index (the ratio of 4-pyridoxic acid over the sum of pyridoxal and pyridoxal-5'-phosphate), has been positively associated with lung cancer risk in two prospective European studies. However, the extent to which this association translates to more diverse populations is not known.

MATERIALS AND METHODS

For this study, we included 5323 incident lung cancer cases and 5323 controls individually matched by age, sex, and smoking status within each of 20 prospective cohorts from the Lung Cancer Cohort Consortium. Cohort-specific odds ratios (ORs) and 95% confidence intervals (CIs) for the association between PAr and lung cancer risk were calculated using conditional logistic regression and pooled using random-effects models.

RESULTS

PAr was positively associated with lung cancer risk in a dose-response fashion. Comparing the fourth versus first quartiles of PAr resulted in an OR of 1.38 (95% CI: 1.19-1.59) for overall lung cancer risk. The association between PAr and lung cancer risk was most prominent in former smokers (OR: 1.69, 95% CI: 1.36-2.10), men (OR: 1.60, 95% CI: 1.28-2.00), and for cancers diagnosed within 3 years of blood draw (OR: 1.73, 95% CI: 1.34-2.23).

CONCLUSION

Based on pre-diagnostic data from 20 cohorts across 4 continents, this study confirms that increased vitamin B6 catabolism related to inflammation and immune activation is associated with a higher risk of developing lung cancer. Moreover, PAr may be a pre-diagnostic marker of lung cancer rather than a causal factor.

Comprehensive genomic analysis of metastatic pancreatic ductal adenocarcinoma (mPDAC) reveals a significant proportion of clinical actionable aberrations

273

Background: Substantial progress has been made in the understanding of the genomic landscape of PDAC, but the clinical utility of these data remains uncertain. Methods: As part of the BC Cancer Personalized Oncogenomics (POG) and PanGen studies, whole genome analysis and transcriptome sequencing were performed on fresh biopsy and blood sample from 47 mPDAC patients. Genomic findings informed therapy choices including potential eligibility for the CCTG PM.1 molecular basket trial. Results: Cohort consists of 53% male, average age 57.8, 34/47 had ≥2 lines of treatment. 37/47 biopsies were from liver and 26/47 were collected pre-treatment. 8/47 (17%) patients had aberrations with strong evidence of clinical actionability. These include 2 germline BRCA2, 1 germline BRCA1, 1 somatic XRCC2 homozygous deletion with strong COSMIC signature 3, all predictive of platinum sensitivity. Patients with XRCC2 deletion and BRCA1 had over 2 years on FOLFIRINOX. Fusions affecting the NRG1 gene were identified in 3 patients, all with KRAS wildtype tumours, which may confer ERRB inhibitor sensitivity. 1/3 NRG1 fusion patients had thus far been treated with the ERBB inhibitor afatinib, with reduction of CA19-9 from >120,000 to 7246 and dramatic response noted on PET CT imaging one month post treatment. One patient had mismatch repair deficiency, and a high mutational burden, suggestive of immune checkpoint inhibitor sensitivity. In total 85% (40/47) sequenced mPDAC patients had potentially actionable mutations which includes CCTG PM.1 trial eligibility: 24/47 with cell cycle dysregulation (CDK4/6 inhibitor arm), 4/47 with high homologous recombination defects (PARP inhibitor arm), 1/47 ERBB2 amplification (anti-HER2 arm), 2/47 high expression of FGFR1 (Sunitinib arm) and 1/47 with high FLT4 and IGF1R expression (Axitinib arm). Conclusions: More routine use of comprehensive genomic analysis should be considered in mPDAC given finding of high degree of actionability. Importantly, a significant proportion (17%) had findings with strong evidence of clinical impact.(NCT02155621, NCT02869802, NCT03297606)

Faces, locations, and tools: a proposed two-stimulus P300 brain computer interface

OBJECTIVE

Brain computer interface (BCI) technology can be important for those unable to communicate due to loss of muscle control. Given that the P300 Speller provides a relatively slow rate of communication, highly accurate classification is of great importance. Previous studies have shown that alternative stimuli (e.g. faces) can improve BCI speed and accuracy. The present study uses two new alternative stimuli, locations and graspable tools. Functional MRI studies have shown that images of familiar locations produce brain responses in the parahippocampal place area and graspable tools produce brain responses in premotor cortex.

APPROACH

The current studies show that location and tool stimuli produce unique and discriminable brain responses that can be used to improve offline classification accuracy. Experiment 1 presented face stimuli and location stimuli and Experiment 2 presented location and tool stimuli.

MAIN RESULTS

In both experiments, offline results showed that a stimulus specific classifier provided higher accuracy, speed, and bit rate.

SIGNIFICANCE

This study was used to provide preliminary offline support for using unique stimuli to improve speed and accuracy of the P300 Speller. Additional experiments should be conducted to examine the online efficacy of this novel paradigm.

Bidirectional Mendelian randomization to explore the causal relationships between body mass index and polycystic ovary syndrome

STUDY QUESTION

What are the causal relationships between polycystic ovary syndrome (PCOS) and body mass index (BMI)?

SUMMARY ANSWER

Bidirectional Mendelian randomization analyses suggest that increased BMI is causal for PCOS while the reverse is not the case.

WHAT IS KNOWN ALREADY

The contribution of obesity to the pathogenesis of PCOS is controversial. To date, published genetic studies addressing this question have generated conflicting results and have not utilized the full extent of known single nucleotide polymorphisms associated with body mass index (BMI).

STUDY DESIGN, SIZE, DURATION

This cross-sectional Mendelian randomization (MR) and genetic association study was conducted in 750 individuals of European origin and with PCOS and 1567 BMI-matched controls.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Cases and controls were matched for BMI as well as for distribution of weight categories (normal weight, overweight, obese). Two-sample MR using inverse variance weighting (IVW) was conducted using a 92-SNP instrument variable for BMI with PCOS as the outcome, followed by two-sample MR with a 16-SNP instrument variable for PCOS with BMI as the outcome. Sensitivity analyses included MR-Egger and maximum likelihood methods. Secondary analyses assessed associations of genetic risk scores and individual SNPs with PCOS, BMI and quantitative androgen-related and glucose homeostasis-related traits.

MAIN RESULTS AND THE ROLE OF CHANCE

Each standard deviation genetically higher BMI was associated with a 4.89 (95% CI 1.46-16.32) higher odds of PCOS. Conversely, genetic risk of PCOS did not influence BMI. Sensitivity analyses yielded directionally consistent results. The genetic risk score of 92 BMI SNPs was associated with the diagnosis of PCOS (OR 1.043, 95% CI 1.009-1.078, P = 0.012). Of the 92 BMI risk variants evaluated, none were associated individually with PCOS after considering multiple testing. The association of FTO SNP rs1421085 with BMI was stronger in women with PCOS (β = 0.071, P = 0.0006) than in controls (β = 0.046, P = 0.065).

LIMITATIONS, REASONS FOR CAUTION

The current sample size, while providing good power for MR and genetic risk score analyses, had limited power to demonstrate association of individual SNPs with PCOS. Cases and controls were not matched for age; however, this was mitigated by adjusting analyses for age. Dietary and lifestyle data, which could have been used to explore the greater association of the FTO SNP with BMI in women with PCOS, was not available.

WIDER IMPLICATIONS OF THE FINDINGS

Increasing BMI appears to be causal for PCOS but having PCOS does not appear to affect BMI. This study used the most comprehensive set of SNPs for BMI currently available. Prior studies using fewer SNPs had yielded conflicting results and may have been confounded because cases and controls were not matched for weight categories. The current results highlight the potential utility of weight management in the prevention and treatment of PCOS.

STUDY FUNDING/COMPETING INTEREST(S)

National Institutes of Health Grants R01-HD29364 and K24-HD01346 (to R.A.), Grant R01-DK79888 (to M.O.G.), Grant U54-HD034449 (to R.S.L.), Grant U19-HL069757 (to R.M.K.). The funders had no influence on the data collection, analyses or conclusions of the study. No conflict of interests to declare.

TRIAL REGISTRATION NUMBER

N/A.

A collaborative filtering-based approach to biomedical knowledge discovery

Motivation

The increase in publication rates makes it challenging for an individual researcher to stay abreast of all relevant research in order to find novel research hypotheses. Literature-based discovery methods make use of knowledge graphs built using text mining and can infer future associations between biomedical concepts that will likely occur in new publications. These predictions are a valuable resource for researchers to explore a research topic. Current methods for prediction are based on the local structure of the knowledge graph. A method that uses global knowledge from across the knowledge graph needs to be developed in order to make knowledge discovery a frequently used tool by researchers.

Results

We propose an approach based on the singular value decomposition (SVD) that is able to combine data from across the knowledge graph through a reduced representation. Using cooccurrence data extracted from published literature, we show that SVD performs better than the leading methods for scoring discoveries. We also show the diminishing predictive power of knowledge discovery as we compare our predictions with real associations that appear further into the future. Finally, we examine the strengths and weaknesses of the SVD approach against another well-performing system using several predicted associations.

Availability and implementation

All code and results files for this analysis can be accessed at https://github.com/jakelever/knowledgediscovery.

Contact

sjones@bcgsc.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Abstract 4340: Integrating whole genome and transcriptome analysis to inform treatment decisions in the metastatic cancer clinical setting

Genomic analysis is an established clinical tool for positioning patients to specific treatment rationales. However, given the vast range of genomic complexity that is observed between tumors the extent to which genomic characterization must be applied in order to generate a true picture of the actionable and informative genomic events in an individual's tumour is still a matter for intense debate. Historically, cost and the complexity of the bioinformatics analysis needed has limited the adoption of comprehensive genomic characterization. However with the emerging clinical significance of genomic signatures, mutation burden and immunogenomic profiling, information that can best be recovered from comprehensive genomic analyses, a better understanding of the utility of extensive genomic characterization of tumours in a clinical setting is needed. For the past 6 years, the Personalized Oncogenomics Project (POG) at the BC Cancer Agency has used paired tumour/normal Whole Genome and Transcriptome Analysis (WGTA) to comprehensively characterize tumour samples from patients with incurable metastatic cancer to inform treatment planning for each patient within a clinically relevant timeframe. To date tumours from over 700 patients have been genomically profiled using this approach. These data represent a cohort of tumours from the post-treatment metastatic setting that have been comprehensively characterized at a whole genome level augmented by comprehensive clinical annotation, including information about the primary diagnosis, progression of the disease and treatment history, as well as follow-up information including the response of the tumour to genome analysis-informed therapy. We will discuss findings from this powerful genomic cohort focusing on aspects that are particularly well represented by using a WGTA approach including; mutation signatures of biological and clinical significance, integration of expression data to complement genomic findings and highlighting the discovery of novel NRG1 gene fusions as a demonstration of how using a comprehensive multi-omic approach for genomic analysis to infer the oncogenic impact of novel genomic events can impact patient care.

Citation Format: Martin R. Jones, Yaoqing Shen, Erin Pleasance, Elisa Majounie, Laura Williamson, Eric Zhao, Eric Chuah, Karen L. Mungall, Andrew J. Mungall, Richard A. Moore, Yussanne Ma, Stephen Yip, Howard Lim, Daniel Renouf, Steven J. Jones, Janesssa Jaskin, Marco A. Marra. Integrating whole genome and transcriptome analysis to inform treatment decisions in the metastatic cancer clinical setting [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 4340.

Successful targeting of the NRG1 pathway indicates novel treatment strategy for metastatic cancer

Background

NRG1 fusion-positive lung cancers have emerged as potentially actionable events in lung cancer, but clinical support is currently limited and no evidence of efficacy of this approach in cancers beyond lung has been shown.

Patients and methods

Here, we describe two patients with advanced cancers refractory to standard therapies. Patient 1 had lung adenocarcinoma and patient 2 cholangiocarcinoma. Whole-genome and transcriptome sequencing were carried out for these cases with select findings validated by fluorescence in situ hybridization.

Results

Both tumors were found to be positive for NRG1 gene fusions. In patient 1, an SDC4-NRG1 gene fusion was detected, similar gene fusions having been described in lung cancers previously. In patient 2, a novel ATP1B1-NRG1 gene fusion was detected. Cholangiocarcinoma is not a disease type in which NRG1 fusions had been described previously. Integrative genome analysis was used to assess the potential functional significance of the detected genomic events including the gene fusions, prioritizing therapeutic strategies targeting the HER-family of growth factor receptors. Both patients were treated with the pan HER-family kinase inhibitor afatinib and both displayed significant and durable response to treatment. Upon progression sites of disease were sequenced. The lack of obvious genomic events to describe the disease progression indicated that broad transcriptomic or epigenetic mechanisms could be attributed to the lack of prolonged response to afatinib.

Conclusion

These observations lend further support to the use of pan HER-tyrosine kinase inhibitors for the treatment of NRG1 fusion-positive in both cancers of lung and hepatocellular origin and indicate more broadly that cancers found to be NRG1 fusion-positive may benefit from such a clinical approach regardless of their site of origin.

Clinical trial information

Personalized Oncogenomics (POG) Program of British Columbia: Utilization of Genomic Analysis to Better Understand Tumour Heterogeneity and Evolution (NCT02155621).

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.¹ 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.⁶ 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.

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.

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.

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.

Abstract A184: Clinical application of whole genome and transcriptome sequencing in cancer care

Aim: To investigate feasibility and utility of whole genome and RNA sequencing in cancer care. The use of comprehensive personal genomic information to guide cancer treatment decisions has gained momentum in recent years. Here we present the Personalized Oncogenomics (POG) project launched at British Columbia Cancer Agency in 2012. This project uses paired tumor/normal whole genome and transcriptome sequencing to characterize a patient's tumor and inform treatment options within a clinically relevant time frame. In the past 5 years, around 600 patients with metastatic cancers (including 49 pediatric cases) have been sequenced and analyzed, representing over 50 cancer types. The breadth and depth of data in POG enable the detection of multiple types of alterations, from simple mutations, indels, and copy number changes, to more complex alterations such as gene fusion, mutation signatures, and homologous recombination deficiency score. Incorporation of gene expression data through transcriptome sequencing informs on the impact of observed genomic alterations, and provides information regarding specific diagnosis via expression comparison to other cancer subtypes. All genomic and transcriptomic variants are integrated to build a personalized tumor-specific molecular profile, identify actionable items supported by an in-house knowledgebase and publically available molecular oncology, and characterize each individual tumor by intensive pathway analysis and literature search. Such multidimensional data also impose challenges in interpretation and communication. We developed a pipeline to translate complex genomic data into clinically actionable and hypothetical recommendations for the treatment of individual patients. The pipeline produces two reports: a panel-like report that contains known SNVs and fusions with therapeutic relevance from a collection of more than 4000 events, and a second comprehensive and manually curated report to fully characterize the tumor using the whole genome and transcriptome datasets. Genomic data are presented and discussed at a multidisciplinary molecular tumor board consisting of medical oncologists, pathologists, bioinformaticians, geneticists, and biologists. This approach has enabled clinicians to make informed clinical decisions based on the genomic data integrated with other clinical features, as well as to form new treatment-related hypothesis. POG has shown that use of both whole genome and transcriptome sequencing allows identification of therapeutic targets in a significant proportion of patients. Almost 80% cases were found to have one or more actionable alterations (100% in pediatric cases), and almost one-third of these are defined only using RNA data. Based on molecular tumor board discussion, patients are directed to clinical trials, positioned to genomically informed standard-of-care options, or treated with off-label drugs. With demonstrated effectiveness, the integrative approach developed by POG not only provides molecular insight and treatment options into individual tumors, but also provides a rich resource of molecular data with matched clinical information that will aid our understanding of tumor biology and therapy response mechanisms to inform treatment strategies in the future.

Citation Format: Yaoqing Shen, Martin R. Jones, Erin Pleasance, Melika Bonakdar, Carolyn Ch'ng, Caralyn Reisle, Laura Williamson, Elisa Majounie, Greg Taylor, Simon Chan, Brandon Pierce, Wei Zhang, Amir Muhammadzadeh, Eric Y. Zhao, Dustin Bleile, Karen Mungall, Nina Thiessen, Eric Chuah, Tina Wong, Richard Corbett, Yussanne Ma, Richard A. Moore, Andrew J. Mungall, Yongjun Zhao, Stephen Yip, Anna F. Lee, Rod Rassekh, Rebecca Deyell, Howard Lim, Daniel Renouf, Robyn Roscoe, Steven J.M Jones, Janessa Laskin, Marco A. Marra. Clinical application of whole genome and transcriptome sequencing in cancer care [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr A184.

High performance computation of landscape genomic models including local indicators of spatial association

With the increasing availability of both molecular and topo-climatic data, the main challenges facing landscape genomics - that is the combination of landscape ecology with population genomics - include processing large numbers of models and distinguishing between selection and demographic processes (e.g. population structure). Several methods address the latter, either by estimating a null model of population history or by simultaneously inferring environmental and demographic effects. Here we present samβada, an approach designed to study signatures of local adaptation, with special emphasis on high performance computing of large-scale genetic and environmental data sets. samβada identifies candidate loci using genotype-environment associations while also incorporating multivariate analyses to assess the effect of many environmental predictor variables. This enables the inclusion of explanatory variables representing population structure into the models to lower the occurrences of spurious genotype-environment associations. In addition, samβada calculates local indicators of spatial association for candidate loci to provide information on whether similar genotypes tend to cluster in space, which constitutes a useful indication of the possible kinship between individuals. To test the usefulness of this approach, we carried out a simulation study and analysed a data set from Ugandan cattle to detect signatures of local adaptation with samβada, bayenv, lfmm and an F_ST outlier method (FDIST approach in arlequin) and compare their results. samβada - an open source software for Windows, Linux and Mac OS X available at http://lasig.epfl.ch/sambada - outperforms other approaches and better suits whole-genome sequence data processing.

Abstract 2473: Breast cancer whole genomes link homologous recombination deficiency (HRD) with therapeutic outcomes

Background: Homologous recombination deficiency (HRD) is common in cancer - germline BRCA1 & BRCA2 mutations account for 5-10% of breast cancers and confer 85% lifetime risk. HRD cancers exhibit genomic instability and sensitivity to platinum-based therapy and PARP inhibitors. While not all causes of HRD are known, recent sequencing efforts have revealed genome-wide somatic mutation signatures that characterize the HRD genomic instability phenotype, also known as “BRCA-ness”. This provides a promising new assay to predict sensitivity to platinum-based therapy. Here, we integrate two whole-genome sequencing metrics to assess their association with therapeutic outcomes in a breast cancer cohort.

Methods: Whole-genome sequencing of 47 breast cancer tumors (100x coverage) and matched normals (60x) was performed on an Illumina HiSeq. Alignment, assembly, SNV calling, and loss of heterozygosity (LOH) detection were performed with BWA, ABySS, Strelka, and APOLLOH respectively. SNV signatures were deciphered by non-negative matrix factorization with Monte Carlo resampling. An HRD score comprised of LOH, telomeric allelic imbalance (TAI), and large scale transition (LST) counts was computed. Clinical endpoints were obtained by retrospective review of treatment and imaging reports. Analysis is ongoing in an independent validation cohort of 62 sequenced cases.

Results: The HRD-linked SNV signature was significantly associated with radiographic clinical response (CR) to platinum-based therapy (p=0.015). Logistic regression demonstrated a 59% improved odds of CR to platinum-based therapy per 1000 somatic SNVs attributed to HRD (odds ratio 1.16-2.50). Tumors carried up to 10,246 such SNVs and all patients with CR were among the top quartile. The LOH-TAI-LST score was correlated with SNV signature (r=0.6, p=7×10-6) and associated with CR (p=0.025). Notably, elevated HRD signatures associated with CR were identified in tumors with wild-type BRCA1/BRCA2 or variants of unknown significance. Tumors with above median HRD signatures were associated with a 69-day longer time to treatment failure and an 18% daily decreased probability of treatment failure per 1000 HRD-attributed SNVs (hazard ratio 0.71-0.95, p = 0.007).

Discussion: We found that HRD mutation signatures are associated with clinical response and longer time to treatment failure with platinum-based therapy. While similar benefits were observed in patients with somatic bi-allelic loss of BRCA1/BRCA2, such cases are less common (8% of our cohort) compared to those with elevated HRD signature. Thus, mutation signature methods may identify patients who stand to benefit from platinum-based therapy missed by BRCA screening alone.

Citation Format: Eric Y. Zhao, Yaoqing Shen, Erin Pleasance, Katayoon Kasaian, Martin R. Jones, Carolyn Ch'ng, Caralyn Reisle, Peter Eirew, Karen Mungall, Nina Thiessen, Yussanne Ma, Alexandra Fok, Andrew J. Mungall, Yongjun Zhao, Richard Moore, Diego Villa, Tamara Shenkier, Caroline Lohrisch, Stephen Chia, Stephen Yip, Karen Gelmon, Howard Lim, Sophie Sun, Kasmintan A. Schrader, Sean Young, Aly Karsan, Robyn Roscoe, Janessa Laskin, Marco A. Marra, Steven J. Jones. Breast cancer whole genomes link homologous recombination deficiency (HRD) with therapeutic outcomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2473. doi:10.1158/1538-7445.AM2017-2473

How close are we to complete annotation of metabolomes?

The metabolome describes the full complement of the tens to hundreds of thousands of low molecular weight metabolites present within a biological system. Identification of the metabolome is critical for discovering the maximum amount of biochemical knowledge from metabolomics datasets. Yet no exhaustive experimental characterisation of any organismal metabolome has been reported to date, dramatically contrasting with the genome sequencing of thousands of plants, animals and microbes. Here, we review the status of metabolome annotation and describe advances in the analytical methodologies being applied. In part through new international coordination, we conclude that we are now entering a new era of metabolome annotation.

High performance computation of landscape genomic models including local indicators of spatial association

With the increasing availability of both molecular and topo-climatic data, the main challenges facing landscape genomics - that is the combination of landscape ecology with population genomics - include processing large numbers of models and distinguishing between selection and demographic processes (e.g. population structure). Several methods address the latter, either by estimating a null model of population history or by simultaneously inferring environmental and demographic effects. Here we present samβada, an approach designed to study signatures of local adaptation, with special emphasis on high performance computing of large-scale genetic and environmental data sets. samβada identifies candidate loci using genotype-environment associations while also incorporating multivariate analyses to assess the effect of many environmental predictor variables. This enables the inclusion of explanatory variables representing population structure into the models to lower the occurrences of spurious genotype-environment associations. In addition, samβada calculates local indicators of spatial association for candidate loci to provide information on whether similar genotypes tend to cluster in space, which constitutes a useful indication of the possible kinship between individuals. To test the usefulness of this approach, we carried out a simulation study and analysed a data set from Ugandan cattle to detect signatures of local adaptation with samβada, bayenv, lfmm and an F_ST outlier method (FDIST approach in arlequin) and compare their results. samβada - an open source software for Windows, Linux and Mac OS X available at http://lasig.epfl.ch/sambada - outperforms other approaches and better suits whole-genome sequence data processing.

Muscle spindle alterations precede onset of sensorimotor deficits in Charcot-Marie-Tooth type 2E

Charcot-Marie-Tooth (CMT) is the most common inherited peripheral neuropathy, affecting approximately 2.8 million people. The CMT leads to distal neuropathy that is characterized by reduced motor nerve conduction velocity, ataxia, muscle atrophy and sensory loss. We generated a mouse model of CMT type 2E (CMT2E) expressing human neurofilament light E396K (hNF-L^E396K ), which develops decreased motor nerve conduction velocity, ataxia and muscle atrophy by 4 months of age. Symptomatic hNF-L^E396K mice developed phenotypes that were consistent with proprioceptive sensory defects as well as reduced sensitivity to mechanical stimulation, while thermal sensitivity and auditory brainstem responses were unaltered. Progression from presymptomatic to symptomatic included a 50% loss of large diameter sensory axons within the fifth lumbar dorsal root of hNF-L^E396K mice. Owing to proprioceptive deficits and loss of large diameter sensory axons, we analyzed muscle spindle morphology in presymptomatic and symptomatic hNF-L^E396K and hNF-L control mice. Muscle spindle cross-sectional area and volume were reduced in all hNF-L^E396K mice analyzed, suggesting that alterations in muscle spindle morphology occurred prior to the onset of typical CMT pathology. These data suggested that CMT2E pathology initiated in the muscle spindles altering the proprioceptive sensory system. Early sensory pathology in CMT2E could provide a unifying hypothesis for the convergence of pathology observed in CMT.

Unified Model of Nanosecond Charge Recombination in Closed Reaction Centers from Rhodobacter sphaeroides: Role of Protein Polarization Dynamics

Ongoing questions surround the influence of protein dynamics on rapid processes such as biological electron transfer. Such questions are particularly addressable in light-activated systems. In Rhodobacter sphaeroides reaction centers, charge recombination or back electron transfer from the reduced bacteriopheophytin, HA(-), to the oxidized dimeric bacteriochlorophyll, P(+), may be monitored by both transient absorption spectroscopy and transient fluorescence spectroscopy. Signals measured with both these techniques decay in a similar three-exponential fashion with lifetimes of ∼0.6-0.7, ∼2-4, and ∼10-20 ns, revealing the complex character of this electron transfer reaction. In this study a single kinetic model was developed to connect lifetime and amplitude data from both techniques. The model took into account the possibility that electron transfer from HA(-) to P(+) may occur with transient formation of the state P(+)BA(-). As a result it was possible to model the impact of nanosecond protein relaxation on the free energy levels of both P(+)HA(-) and P(+)BA(-) states relative to that of the singlet excited state of P, P*. Surprisingly, whereas the free energy gap between P* and P(+)HA(-) increased with time in response to protein reorganization, the free energy gap between P* and P(+)BA(-) decreased. This finding may be accounted for by a gradual polarization of the protein environment which stabilizes the state P(+)HA(-) and destabilizes the state P(+)BA(-), favoring productive charge separation over unproductive charge recombination.

Response to angiotensin blockade with irbesartan in a patient with metastatic colorectal cancer

BACKGROUND

A patient suffering from metastatic colorectal cancer, treatment-related toxicity and resistance to standard chemotherapy and radiation was assessed as part of a personalized oncogenomics initiative to derive potential alternative therapeutic strategies.

PATIENTS AND METHODS

Whole-genome and transcriptome sequencing was used to interrogate a metastatic tumor refractory to standard treatments of a patient with mismatch repair-deficient metastatic colorectal cancer.

RESULTS

Integrative genomic analysis indicated overexpression of the AP-1 transcriptional complex suggesting experimental therapeutic rationales, including blockade of the renin-angiotensin system. This led to the repurposing of the angiotensin II receptor antagonist, irbesartan, as an anticancer therapy, resulting in the patient experiencing a dramatic and durable response.

CONCLUSIONS

This case highlights the utility of comprehensive integrative genomic profiling and bioinformatics analysis to provide hypothetical rationales for personalized treatment options.

Non-targeted UHPLC-MS metabolomic data processing methods: a comparative investigation of normalisation, missing value imputation, transformation and scaling

INTRODUCTION

The generic metabolomics data processing workflow is constructed with a serial set of processes including peak picking, quality assurance, normalisation, missing value imputation, transformation and scaling. The combination of these processes should present the experimental data in an appropriate structure so to identify the biological changes in a valid and robust manner.

OBJECTIVES

Currently, different researchers apply different data processing methods and no assessment of the permutations applied to UHPLC-MS datasets has been published. Here we wish to define the most appropriate data processing workflow.

METHODS

We assess the influence of normalisation, missing value imputation, transformation and scaling methods on univariate and multivariate analysis of UHPLC-MS datasets acquired for different mammalian samples.

RESULTS

Our studies have shown that once data are filtered, missing values are not correlated with m/z, retention time or response. Following an exhaustive evaluation, we recommend PQN normalisation with no missing value imputation and no transformation or scaling for univariate analysis. For PCA we recommend applying PQN normalisation with Random Forest missing value imputation, glog transformation and no scaling method. For PLS-DA we recommend PQN normalisation, KNN as the missing value imputation method, generalised logarithm transformation and no scaling. These recommendations are based on searching for the biologically important metabolite features independent of their measured abundance.

CONCLUSION

The appropriate choice of normalisation, missing value imputation, transformation and scaling methods differs depending on the data analysis method and the choice of method is essential to maximise the biological derivations from UHPLC-MS datasets.

let-65 is cytoplasmic methionyl tRNA synthetase in C. elegans

Cytoplasmic methionyl tRNA synthetase (MetRS) is one of more than 20 cytoplasmic aminoacyl tRNA synthetase enzymes (ARS). This family of enzymes catalyzes a process fundamental for protein translation. Using a combination of genetic mapping, oligonucleotide array comparative genomic hybridization, and phenotypic correlation, we show that mutations in the essential gene, let-65, reside within the predicted Caenorhabditis elegans homologue of MetRS, which we have named mars-1. We demonstrate that the lethality associated with alleles of let-65 is fully rescued by a transgenic array that spans the mars-1 genomic region. Furthermore, sequence analysis reveals that six let-65 alleles lead to the alteration of highly conserved amino acids.

Frequency-domain optical coherence tomography assessment of human carotid atherosclerosis using saline flush for blood clearance without balloon occlusion

FD-OCT is a new imaging technique that allows unprecedented in vivo microlevel assessment of human carotid plaque morphologic patterns and stent-vessel interactions. Prior reports describing the use of this technique have used balloon occlusion of the target vessel or iodinated contrast media to facilitate imaging. We report, for the first time, in vivo FD-OCT imaging of human carotid arteries without the use of iodinated contrast material or balloon occlusion techniques.

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. genesis 51:545–561.

Systematic evaluation of validated type 2 diabetes and glycaemic trait loci for association with insulin clearance

AIMS/HYPOTHESIS

Insulin clearance is a highly heritable trait, for which few quantitative trait loci have been discovered. We sought to determine whether validated type 2 diabetes and/or glycaemic trait loci are associated with insulin clearance.

METHODS

Hyperinsulinaemic-euglycaemic clamps were performed in two Hispanic-American family cohorts totalling 1329 participants in 329 families. The Metabochip was used to fine-map about 50 previously identified loci for type 2 diabetes, fasting glucose, fasting insulin, 2 h glucose or HbA1c. This resulted in 17,930 variants, which were tested for association with clamp-derived insulin clearance via meta-analysis of the two cohorts.

RESULTS

In the meta-analysis, 38 variants located within seven loci demonstrated association with insulin clearance (p < 0.001). The top signals for each locus were rs10241087 (DGKB/TMEM195 [TMEM195 also known as AGMO]) (p = 4.4 × 10(-5)); chr1:217605433 (LYPLAL1) (p = 3.25 × 10(-4)); rs2380949 (GLIS3) (p = 3.4 × 10(-4)); rs55903902 (FADS1) (p = 5.6 × 10(-4)); rs849334 (JAZF1) (p = 6.4 × 10(-4)); rs35749 (IGF1) (p = 6.7 × 10(-4)); and rs9460557 (CDKAL1) (p = 6.8 × 10(-4)).

CONCLUSIONS/INTERPRETATION

While the majority of validated loci for type 2 diabetes and related traits do not appear to influence insulin clearance in Hispanics, several of these loci do show evidence of association with this trait. It is therefore possible that these loci could have pleiotropic effects on insulin secretion, insulin sensitivity and insulin clearance.