Evaluation of Endobronchial Ultrasound-Guided Transbronchial Needle Aspiration (EBUS-TBNA) Samples from Advanced Non-Small Cell Lung Cancer for Whole Genome, Whole Exome and Comprehensive Panel Sequencing

Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2-3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

The clinical utility and costs of whole-genome sequencing to detect cancer susceptibility variants-a multi-site prospective cohort study

BACKGROUND

Many families and individuals do not meet criteria for a known hereditary cancer syndrome but display unusual clusters of cancers. These families may carry pathogenic variants in cancer predisposition genes and be at higher risk for developing cancer.

METHODS

This multi-centre prospective study recruited 195 cancer-affected participants suspected to have a hereditary cancer syndrome for whom previous clinical targeted genetic testing was either not informative or not available. To identify pathogenic disease-causing variants explaining participant presentation, germline whole-genome sequencing (WGS) and a comprehensive cancer virtual gene panel analysis were undertaken.

RESULTS

Pathogenic variants consistent with the presenting cancer(s) were identified in 5.1% (10/195) of participants and pathogenic variants considered secondary findings with potential risk management implications were identified in another 9.7% (19/195) of participants. Health economic analysis estimated the marginal cost per case with an actionable variant was significantly lower for upfront WGS with virtual panel ($8744AUD) compared to standard testing followed by WGS ($24,894AUD). Financial analysis suggests that national adoption of diagnostic WGS testing would require a ninefold increase in government annual expenditure compared to conventional testing.

CONCLUSIONS

These findings make a case for replacing conventional testing with WGS to deliver clinically important benefits for cancer patients and families. The uptake of such an approach will depend on the perspectives of different payers on affordability.

Multi-omic features of oesophageal adenocarcinoma in patients treated with preoperative neoadjuvant therapy

Oesophageal adenocarcinoma is a poor prognosis cancer and the molecular features underpinning response to treatment remain unclear. We investigate whole genome, transcriptomic and methylation data from 115 oesophageal adenocarcinoma patients mostly from the DOCTOR phase II clinical trial (Australian New Zealand Clinical Trials Registry-ACTRN12609000665235), with exploratory analysis pre-specified in the study protocol of the trial. We report genomic features associated with poorer overall survival, such as the APOBEC mutational and RS3-like rearrangement signatures. We also show that positron emission tomography non-responders have more sub-clonal genomic copy number alterations. Transcriptomic analysis categorises patients into four immune clusters correlated with survival. The immune suppressed cluster is associated with worse survival, enriched with myeloid-derived cells, and an epithelial-mesenchymal transition signature. The immune hot cluster is associated with better survival, enriched with lymphocytes, myeloid-derived cells, and an immune signature including CCL5, CD8A, and NKG7. The immune clusters highlight patients who may respond to immunotherapy and thus may guide future clinical trials.

streammd: fast low-memory duplicate marking using a Bloom filter

SUMMARY

Identification of duplicate templates is a common pre-processing step in bulk sequence analysis; for large libraries this can be resource intensive. Here we present streammd: a fast, memory-efficient, single-pass duplicate marker operating on the principle of a Bloom filter. streammd closely reproduces outputs from Picard MarkDuplicates while being substantially faster, and requires much less memory than SAMBLASTER.

AVAILABILITY AND IMPLEMENTATION

streammd is a C ++ program available from GitHub https://github.com/delocalizer/streammd under the MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

1234 IMPROVING HEART FAILURE MANAGEMENT WITHIN HOSPITAL AT HOME

Introduction

Heart failure (HF) is a common problem managed in our West Lothian multi-disciplinary hospital at home (HaH) service, however significant variation in practice was noted with considerable resource implications. We aimed to standardise and improve this by developing a dedicated protocol.

Methods

We developed a protocol to guide the assessment and management of HF within HaH. We collected baseline (n=25) and follow-up data (n=10) after protocol introduction from patients referred to HaH with heart failure. Outcomes reviewed included anticipatory care planning (ACP) decisions, length of stay (LOS) and treatment strategy. We held staff education sessions and surveyed staff confidence regarding HF management.

Results

ACP discussion rates improved after protocol introduction, with decision rates improving for both escalation of care (28% to 80%) and resuscitation (44% to 60%). LOS reduced after protocol introduction (mean 6.3 days to 5.9 days). Titration of oral diuretics alone (71%) was associated with a shorter LOS (mean 5.4 days) compared to IV (29%, mean 8.1 days), with no difference in 28 day outcome. In those with HF with reduced ejection fraction, the rates of beta-blocker prescription increased (57% to 80%) however ACE-inhibitor prescription decreased (29% to 20%). Use of add-on therapy (e.g. thiazide diuretics) increased (12% to 30%) with a decrease in complication rates (12% to 0%). All staff found the protocol helpful with an improvement in confidence levels.

Conclusions

Through introducing a standardised protocol, we observed an improvement in anticipatory care discussion rates and a trend towards shorter LOS. Oral diuretic titration was less resource intensive without an adverse impact on outcome. Future plans include ongoing education and data collection, trialling a joint multi-disciplinary meeting with cardiology for discussion of complex patients and embedding a treatment strategy of oral diuretic titration with a ‘discharge with planned review' approach in appropriate patients.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes.

SIGNIFICANCE

This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

Analysis of hereditary cancer gene variant classifications from ClinVar indicates a need for regular reassessment of clinical assertions

The clinical classification of variants may change with new information, however, there is limited guidance on how often significant changes in variant classification occur. We used ClinVar to examine how variant classification changes over time. We developed a custom parser and accessed variant data from ClinVar between January 2015 and July 2021. The ClinVar-assigned "aggregate" classification of variants in 121 hereditary cancer genes was harmonized across releases to align to the American College of Medical Genetics and Genomics and the Association for Molecular Pathology terms. Aggregate classification categories were grouped as: benign/likely benign (B/LB); likely pathogenic/pathogenic (LP/P); variant of uncertain significance (VUS); conflicting interpretations of pathogenicity (Conflicting); or Other. We profiled changes in aggregate variant classification between consecutive semi-annual ClinVar releases. The proportion of variants that changed aggregate classification between semi-annual ClinVar releases ranged from 0.6% to 6.4%. The most frequent changes were "VUS to conflicting," "other to LP/P," and "B/LB to Conflicting." A limited number of variants changed aggregate classification from "LP/P to B/LB," or vice versa. Our analysis indicates need for regular reassessment of clinical variant interpretations. The parser developed for this project will facilitate extraction of relevant interpretation data from ClinVar.

43 EXPLORING STAFF PERCEPTIONS AND ATTITUDES TO CARE AND CAREBOTS: THE CASE OF AN ORIGINAL AIR-DISINFECTION ROBOT IN IRELAND

Background

The pandemic heightened the sense of security and safety in care settings, with cleanliness and infection control becoming an even more important aspect of care delivery. Although the impact of new welfare technology on health and social care has been discussed in recent years, few studies have reported the implementation processes or human-robot interactions in care facilities in different cultural settings. The aims of this interdisciplinary research therefore were to understand needs and aspects that have to be considered for implementing an assistive robot, and to explore user perceptions, and the process by which the robots are adopted in Ireland and Japan.

Methods

An original air-disinfection robot (V-Air) was developed by Akara Robotics, as part of research project “Harmonisation towards the establishment of Person-centred, Robotics-aided Care System” (Toyota Foundation, D18-ST-0005). Prior to its instalment in a rehabilitation centre in Dublin, semi-structured interviews were conducted, and observations carried out with seven care professionals, as they interacted with V-Air. The robot was then trialled for eight weeks (March-May 2022) in selected areas of the facility. After four weeks, the users filled out the System Usability Score and were asked questions at the end of the trial. Their scores and answers to the questions revealed staff perceptions and attitudes to care robots.

Results

Overall, the users had positive perceptions of V-Air and its usability. Initial differences existed in staff’s confidence levels, depending on prior experiences with technologies. Collective sensemaking was observed, particularly, around care delivery processes and robot functionalities. The adoption process was facilitated by several factors such as the timing of introduction, user-centred design, concept of care and organisational support.

Conclusion

The findings suggest that the introduction of care robots in care settings can offer an additional layer of organisational safety, while highlighting the significance of the iterative process in adopting assistive technologies.

Comparison of actionable events detected in cancer genomes by whole-genome sequencing, in silico whole-exome and mutation panels

Background

Next-generation sequencing is used in cancer research to identify somatic and germline mutations, which can predict sensitivity or resistance to therapies, and may be a useful tool to reveal drug repurposing opportunities between tumour types. Multigene panels are used in clinical practice for detecting targetable mutations. However, the value of clinical whole-exome sequencing (WES) and whole-genome sequencing (WGS) for cancer care is less defined, specifically as the majority of variants found using these technologies are of uncertain significance.

Patients and methods

We used the Cancer Genome Interpreter and WGS in 726 tumours spanning 10 cancer types to identify drug repurposing opportunities. We compare the ability of WGS to detect actionable variants, tumour mutation burden (TMB) and microsatellite instability (MSI) by using in silico down-sampled data to mimic WES, a comprehensive sequencing panel and a hotspot mutation panel.

Results

We reveal drug repurposing opportunities as numerous biomarkers are shared across many solid tumour types. Comprehensive panels identify the majority of approved actionable mutations, with WGS detecting more candidate actionable mutations for biomarkers currently in clinical trials. Moreover, estimated values for TMB and MSI vary when calculated from WGS, WES and panel data, and are dependent on whether all mutations or only non-synonymous mutations were used. Our results suggest that TMB and MSI thresholds should not only be tumour-dependent, but also be sequencing platform-dependent.

Conclusions

There is a large opportunity to repurpose cancer drugs, and these data suggest that comprehensive sequencing is an invaluable source of information to guide clinical decisions by facilitating precision medicine and may provide a wealth of information for future studies. Furthermore, the sequencing and analysis approach used to estimate TMB may have clinical implications if a hard threshold is used to indicate which patients may respond to immunotherapy.

•

Genome analysis revealed that treatment biomarkers are shared across solid tumours, highlighting repurposing opportunities.

•

Comprehensive panels detect most known biomarkers; however, WGS detects more biomarkers for treatments in clinical trials.

•

TMB is well correlated between sequencing methods, but absolute values vary and are dependent on mutation types considered.

Comprehensive genomic and tumour immune profiling reveals potential therapeutic targets in malignant pleural mesothelioma

Background

Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials.

Methods

We analysed somatic mutations from 229 MPM samples, including previously published data and 58 samples that had undergone WGS within this study. This was combined with RNA-seq analysis to characterize the tumour immune environment.

Results

The comprehensive genome analysis identified 12 driver genes, including new candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed SBS5/40 were dominant in 93% of samples, and defects in homologous recombination repair were infrequent in our cohort. The tumour immune environment contained high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immune suppressive environment. The expression of TGFB1 was associated with overall survival. A small subset of samples (less than 10%) had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a ‘hot’ immune environment independent of the somatic mutations.

Conclusions

We propose accounting for genomic and immune microenvironment status may influence therapeutic planning in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01060-8.

qmotif: determination of telomere content from whole-genome sequence data

Motivation

Changes in telomere length have been observed in cancer and can be indicative of mechanisms involved in carcinogenesis. Most methods used to estimate telomere length require laboratory analysis of DNA samples. Here, we present qmotif, a fast and easy tool that determines telomeric repeat sequences content as an estimate of telomere length directly from whole-genome sequencing.

Results

qmotif shows similar results to quantitative PCR, the standard method for high-throughput clinical telomere length quantification. qmotif output correlates strongly with the output of other tools for determining telomere sequence content, TelSeq and TelomereHunter, but can run in a fraction of the time-usually under a minute.

Availability and implementation

qmotif is implemented in Java and source code is available at https://github.com/AdamaJava/adamajava, with instructions on how to build and use the application available from https://adamajava.readthedocs.io/en/latest/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Multiomic profiling of checkpoint inhibitor-treated melanoma: Identifying predictors of response and resistance, and markers of biological discordance

We concurrently examine the whole genome, transcriptome, methylome, and immune cell infiltrates in baseline tumors from 77 patients with advanced cutaneous melanoma treated with anti-PD-1 with or without anti-CTLA-4. We show that high tumor mutation burden (TMB), neoantigen load, expression of IFNγ-related genes, programmed death ligand expression, low PSMB8 methylation (therefore high expression), and T cells in the tumor microenvironment are associated with response to immunotherapy. No specific mutation correlates with therapy response. A multivariable model combining the TMB and IFNγ-related gene expression robustly predicts response (89% sensitivity, 53% specificity, area under the curve [AUC], 0.84); tumors with high TMB and a high IFNγ signature show the best response to immunotherapy. This model validates in an independent cohort (80% sensitivity, 59% specificity, AUC, 0.79). Except for a JAK3 loss-of-function mutation, for patients who did not respond as predicted there is no obvious biological mechanism that clearly explained their outlier status, consistent with intratumor and intertumor heterogeneity in response to immunotherapy.

Considerations for using population frequency data in germline variant interpretation: Cancer syndrome genes as a model

Aggregate population genomics data from large cohorts are vital for assessing germline variant pathogenicity. However, there are no specifications on how sequencing quality metrics should be considered, and whether exome-derived and genome-derived allele frequencies should be considered in isolation. Germline genome sequence data were simulated for nine read-depths to identify a minimum acceptable read-depth for detecting variants. gnomAD exome-derived and genome-derived datasets were assessed for read-depth, for six key cancer genes selected for variant curation by ClinGen expert panels. Non-Finnish European allele frequency (AF) or filter AF of coding variants in these genes, assigned into frequency bins using modified ACMG-AMP criteria, was compared between exome-derived and genome-derived datasets. A 30X read-depth achieved acceptable precision and recall for detection of substitutions, but poor recall for small insertions/deletions. Exome-derived and genome-derived datasets exhibited low read-depth for different gene exons. Individual variants were mostly assigned to non-divergent AF bins (>95%) or filter AF bins (>97%). Two major bin divergences were resolved by applying the minimal acceptable read-depth threshold. These findings show the importance of assessing read-depth separately for population datasets sourced from different short-read sequencing technologies before assigning a frequency-based ACMG-AMP classification code for variant interpretation.

DNA methylation patterns identify subgroups of pancreatic neuroendocrine tumors with clinical association

Here we report the DNA methylation profile of 84 sporadic pancreatic neuroendocrine tumors (PanNETs) with associated clinical and genomic information. We identified three subgroups of PanNETs, termed T1, T2 and T3, with distinct patterns of methylation. The T1 subgroup was enriched for functional tumors and ATRX, DAXX and MEN1 wild-type genotypes. The T2 subgroup contained tumors with mutations in ATRX, DAXX and MEN1 and recurrent patterns of chromosomal losses in half of the genome with no association between regions with recurrent loss and methylation levels. T2 tumors were larger and had lower methylation in the MGMT gene body, which showed positive correlation with gene expression. The T3 subgroup harboured mutations in MEN1 with recurrent loss of chromosome 11, was enriched for grade G1 tumors and showed histological parameters associated with better prognosis. Our results suggest a role for methylation in both driving tumorigenesis and potentially stratifying prognosis in PanNETs.

Whole-genome sequencing of acral melanoma reveals genomic complexity and diversity

To increase understanding of the genomic landscape of acral melanoma, a rare form of melanoma occurring on palms, soles or nail beds, whole genome sequencing of 87 tumors with matching transcriptome sequencing for 63 tumors was performed. Here we report that mutational signature analysis reveals a subset of tumors, mostly subungual, with an ultraviolet radiation signature. Significantly mutated genes are BRAF, NRAS, NF1, NOTCH2, PTEN and TYRP1. Mutations and amplification of KIT are also common. Structural rearrangement and copy number signatures show that whole genome duplication, aneuploidy and complex rearrangements are common. Complex rearrangements occur recurrently and are associated with amplification of TERT, CDK4, MDM2, CCND1, PAK1 and GAB2, indicating potential therapeutic options.

Using whole-genome sequencing data to derive the homologous recombination deficiency scores

The homologous recombination deficiency (HRD) score was developed using whole-genome copy number data derived from arrays as a way to infer deficiency in the homologous recombination DNA damage repair pathway (in particular BRCA1 or BRCA2 deficiency) in breast cancer samples. The score has utility in understanding tumour biology and may be indicative of response to certain therapeutic strategies. Studies have used whole-exome sequencing to derive the HRD score, however, with increasing use of whole-genome sequencing (WGS) to characterise tumour genomes, there has yet to be a comprehensive comparison between HRD scores derived by array versus WGS. Here we demonstrate that there is both a high correlation and a good agreement between array- and WGS-derived HRD scores and between the scores derived from WGS and downsampled WGS to represent shallow WGS. For samples with an HRD score close to threshold for stratifying HR proficiency or deficiency there was however some disagreement in the HR status between array and WGS data, highlighting the importance of not relying on a single method of ascertaining the homologous recombination status of a tumour.

Whole genome landscapes of uveal melanoma show an ultraviolet radiation signature in iris tumours

Uveal melanoma (UM) is the most common intraocular tumour in adults and despite surgical or radiation treatment of primary tumours, ~50% of patients progress to metastatic disease. Therapeutic options for metastatic UM are limited, with clinical trials having little impact. Here we perform whole-genome sequencing (WGS) of 103 UM from all sites of the uveal tract (choroid, ciliary body, iris). While most UM have low tumour mutation burden (TMB), two subsets with high TMB are seen; one driven by germline MBD4 mutation, and another by ultraviolet radiation (UVR) exposure, which is restricted to iris UM. All but one tumour have a known UM driver gene mutation (GNAQ, GNA11, BAP1, PLCB4, CYSLTR2, SF3B1, EIF1AX). We identify three other significantly mutated genes (TP53, RPL5 and CENPE).

Uveal melanoma has a propensity to metastasise. Here, the authors report the whole genome sequence of 103 uveal melanomas and find that the tumour mutational burden is variable and that two subsets of tumours are characterised by MBD4 mutations and a UV exposure signature.

HNF4A and GATA6 Loss Reveals Therapeutically Actionable Subtypes in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) can be divided into transcriptomic subtypes with two broad lineages referred to as classical (pancreatic) and squamous. We find that these two subtypes are driven by distinct metabolic phenotypes. Loss of genes that drive endodermal lineage specification, HNF4A and GATA6, switch metabolic profiles from classical (pancreatic) to predominantly squamous, with glycogen synthase kinase 3 beta (GSK3β) a key regulator of glycolysis. Pharmacological inhibition of GSK3β results in selective sensitivity in the squamous subtype; however, a subset of these squamous patient-derived cell lines (PDCLs) acquires rapid drug tolerance. Using chromatin accessibility maps, we demonstrate that the squamous subtype can be further classified using chromatin accessibility to predict responsiveness and tolerance to GSK3β inhibitors. Our findings demonstrate that distinct patterns of chromatin accessibility can be used to identify patient subgroups that are indistinguishable by gene expression profiles, highlighting the utility of chromatin-based biomarkers for patient selection in the treatment of PDAC.

Running Genomic Analyses in the Cloud

Genomic testing is rapidly moving into healthcare practice. However it comes with informatics challenges that the healthcare system has not previously faced - the raw data can be hundreds of gigabytes per test, the compute demands can be thousands of CPU hours, and the test can reveal deeply private health-srelated information that can have implications for anyone related to the person tested. While not a panacea, cloud computing has particular properties that can ameliorate some of these difficulties. This paper presents some of the key lessons learned while deploying a set of genomic analyses on cloud computing for Queensland Genomics.

Whole-genome landscape of mucosal melanoma reveals diverse drivers and therapeutic targets

Knowledge of key drivers and therapeutic targets in mucosal melanoma is limited due to the paucity of comprehensive mutation data on this rare tumor type. To better understand the genomic landscape of mucosal melanoma, here we describe whole genome sequencing analysis of 67 tumors and validation of driver gene mutations by exome sequencing of 45 tumors. Tumors have a low point mutation burden and high numbers of structural variants, including recurrent structural rearrangements targeting TERT, CDK4 and MDM2. Significantly mutated genes are NRAS, BRAF, NF1, KIT, SF3B1, TP53, SPRED1, ATRX, HLA-A and CHD8. SF3B1 mutations occur more commonly in female genital and anorectal melanomas and CTNNB1 mutations implicate a role for WNT signaling defects in the genesis of some mucosal melanomas. TERT aberrations and ATRX mutations are associated with alterations in telomere length. Mutation profiles of the majority of mucosal melanomas suggest potential susceptibility to CDK4/6 and/or MEK inhibitors.

Whole-genome sequencing reveals clinically relevant insights into the aetiology of familial breast cancers

BACKGROUND

Whole-genome sequencing (WGS) is a powerful method for revealing the diversity and complexity of the somatic mutation burden of tumours. Here, we investigated the utility of tumour and matched germline WGS for understanding aetiology and treatment opportunities for high-risk individuals with familial breast cancer.

PATIENTS AND METHODS

We carried out WGS on 78 paired germline and tumour DNA samples from individuals carrying pathogenic variants in BRCA1 (n = 26) or BRCA2 (n = 22) or from non-carriers (non-BRCA1/2; n = 30).

RESULTS

Matched germline/tumour WGS and somatic mutational signature analysis revealed patients with unreported, dual pathogenic germline variants in cancer risk genes (BRCA1/BRCA2; BRCA1/MUTYH). The strategy identified that 100% of tumours from BRCA1 carriers and 91% of tumours from BRCA2 carriers exhibited biallelic inactivation of the respective gene, together with somatic mutational signatures suggestive of a functional deficiency in homologous recombination. A set of non-BRCA1/2 tumours also had somatic signatures indicative of BRCA-deficiency, including tumours with BRCA1 promoter methylation, and tumours from carriers of a PALB2 pathogenic germline variant and a BRCA2 variant of uncertain significance. A subset of 13 non-BRCA1/2 tumours from early onset cases were BRCA-proficient, yet displayed complex clustered structural rearrangements associated with the amplification of oncogenes and pathogenic germline variants in TP53, ATM and CHEK2.

CONCLUSIONS

Our study highlights the role that WGS of matched germline/tumour DNA and the somatic mutational signatures can play in the discovery of pathogenic germline variants and for providing supporting evidence for variant pathogenicity. WGS-derived signatures were more robust than germline status and other genomic predictors of homologous recombination deficiency, thus impacting the selection of platinum-based or PARP inhibitor therapy. In this first examination of non-BRCA1/2 tumours by WGS, we illustrate the considerable heterogeneity of these tumour genomes and highlight that complex genomic rearrangements may drive tumourigenesis in a subset of cases.

The new buzz: Investigating the antimicrobial interactions between bioactive compounds found in South African propolis

ETHNOPHARMACOLOGICAL RELEVANCE

Propolis, a resinous substance produced by the Apis mellifera bee, contains a number of flavonoids sourced from plants found in the surrounding region. Whilst bees use this substance to seal off and protect the beehive, humans have used propolis therapeutically for centuries, making use of its antibacterial, antiseptic, antipyretic and wound healing properties, among others. South African propolis is rich in the flavonoids pinocembrin, galangin, and chrysin and very little previous research has been conducted on the antimicrobial effects of these compounds.

AIM OF THE STUDY

To obtain an understanding of the antimicrobial activity of the compounds pinocembrin, galangin, and chrysin, both independently and in combination.

MATERIALS AND METHODS

The compounds pinocembrin, galangin and chrysin were investigated for interactive antimicrobial activity by determining the minimum inhibitory concentrations (MIC), minimum bactericidal concentrations (MBC), anti-quorum sensing activity, biofilm studies, and toxicity studies (brine shrimp lethality assay).

RESULTS

Minimum inhibitory concentration results demonstrated that combinations of compounds showed better inhibitory activity than single compounds. When the flavonoids were tested in combination using the MIC assay, synergy was noted for 22% of the 1:1 ratio combinations and for 66% of the triple 1:1:1 ratio combinations. Similarly, MBC results showed bactericidal activity from selected combinations, while the compounds on their own demonstrated no cidal activity. Quorum sensing studies showed that compound combinations are more effective at inhibiting bacterial communication than the individual compounds. Biofilm assays showed that the highest percentage inhibition was observed for the triple combination against E. coli at 24 h. Finally, brine shrimp lethality studies revealed that combinations of the three compounds had reduced cytotoxicity when compared to the individual compounds.

CONCLUSION

The results obtained in this study demonstrate that the compounds found in South African propolis work synergistically to achieve an optimal antimicrobial effect, whilst simultaneously minimizing cytotoxicity.

Abstract P5-10-01: Using whole genome sequencing and somatic mutation signatures to unravel insight into familial breast cancer aetiology

Approximately 10-15% of breast cancers are associated with a strong family history of disease. Pathogenic variants in BRCA1, BRCA2 or other moderate to highly penetrant susceptibility genes (e.g. TP53, ATM, CHEK2, PALB2 and PTEN) account for a number of breast cancer families. However, for over 50% of families the underlying genetic contribution to their risk remains unknown (termed here as non-BRCA1/2). This has a profound impact for how individuals and their families are managed in the clinic. We applied whole genome sequencing (WGS) to determine whether somatic mutation analysis can reveal insight into the aetiology of familial breast cancer. The full repertoire of somatic mutations was evaluated in 26 BRCA1, 22 BRCA2 and 32 non-BRCA1/2 tumours; including SNPs, indels, copy number changes and structural rearrangements, and mutational signatures. Genomes were also analysed using the HRD Index and HRDetect, as predictors of homologous recombination deficiency. BRCA1, BRCA2 and non-BRCA1/2 tumours exhibited a different burden of mutations, a different spectrum of mutational signatures and different telomere length. Based on collective patterns of mutation signatures, tumours were classified as 'BRCA1-like', 'BRCA2-like' or 'non-BRCA1/2-like' with a 15% rate of tumour re-classification from their original clinical BRCA status. The results demonstrate the power of WGS to differentiate between BRCA1 and BRCA2 driven tumours; in the identification of double-pathogenic germline mutation carriers based on the resulting somatic mutation signature; and in the interpretation of BRCA unclassified variants. WGS of tumour genomes reveals fascinating insights into tumour aetiology and could compliment current genetic testing of breast cancer families.

Citation Format: Simpson P, Nones K, Johnson J, Newell F, Patch A-M, Thorne H, Kazakoff S, De Luca X, Parsons M, Ferguson K, Reid L, McCart Reed A, Srihari S, Lakis V, Davidson A, Mukhopadhyay P, Holmes O, Xu Q, Wood S, Leonard C, Beasley J, Degasperi A, Nik-Zainal S, Ragan M, Spurdle A, Khanna KK, Lakhani S, Pearson J, Chenevix-Trench G, Waddell N. Using whole genome sequencing and somatic mutation signatures to unravel insight into familial breast cancer aetiology [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P5-10-01.

CCL22 is a biomarker of cartilage injury and plays a functional role in chondrocyte apoptosis

BACKGROUND

Osteoarthritis (OA) is one of the leading causes of disability worldwide. Previous history of knee injury is a significant risk factor for OA. It has been established that low-level chronic inflammation plays a pivotal role in the onset and pathogenesis of OA. The primary aim of this research was to determine if a history of knee joint injury is associated with systemic inflammation. A secondary aim was to determine if systemic inflammation is related to knee pain and joint structure.

METHODS

Differences in serum cytokine association networks, knee joint structural changes (MRI), and self-reported pain (i.e., Knee Injury and Osteoarthritis Outcome Score Pain subscale, KOOS_PAIN and Intermittent and Constant Osteoarthritis Pain score, ICOAP) between individuals who had sustained a youth (aged 15-26 years) sport-related knee injury 3-10 years previously and age- and sex-matched controls were examined. Proteins of interest were also examined in an OA rat model.

RESULTS

Cytokine association networks were found to differ significantly between study groups, yet no significant associations were found between networks and KOOS_PAIN or MRI-defined OA. A group of cytokines (MCP1/CCL2, CCL22 and TNFα) were differentially associated with other cytokines between study groups. In a pre-clinical rat OA model, serum CCL22 levels were associated with pain (r = 0.255, p = 0.045) and structural changes to the cartilage. CCL22 expression was also observed in human OA cartilage and furthermore, CCL22 induced apoptosis of isolated human chondrocytes.

DISCUSSION

These results suggest that CCL22 may be an early factor in the onset/pathogenic process of cartilage degeneration and/or related to pain OA.

Phenotypic and molecular dissection of metaplastic breast cancer and the prognostic implications

Metaplastic breast carcinoma (MBC) is relatively rare but accounts for a significant proportion of global breast cancer mortality. This group is extremely heterogeneous and by definition exhibits metaplastic change to squamous and/or mesenchymal elements, including spindle, squamous, chondroid, osseous, and rhabdomyoid features. Clinically, patients are more likely to present with large primary tumours (higher stage), distant metastases, and overall, have shorter 5-year survival compared to invasive carcinomas of no special type. The current World Health Organisation (WHO) diagnostic classification for this cancer type is based purely on morphology - the biological basis and clinical relevance of its seven sub-categories are currently unclear. By establishing the Asia-Pacific MBC (AP-MBC) Consortium, we amassed a large series of MBCs (n = 347) and analysed the mutation profile of a subset, expression of 14 breast cancer biomarkers, and clinicopathological correlates, contextualising our findings within the WHO guidelines. The most significant indicators of poor prognosis were large tumour size (T3; p = 0.004), loss of cytokeratin expression (lack of staining with pan-cytokeratin AE1/3 antibody; p = 0.007), EGFR overexpression (p = 0.01), and for 'mixed' MBC, the presence of more than three distinct morphological entities (p = 0.007). Conversely, fewer morphological components and EGFR negativity were favourable indicators. Exome sequencing of 30 cases confirmed enrichment of TP53 and PTEN mutations, and intriguingly, concurrent mutations of TP53, PTEN, and PIK3CA. Mutations in neurofibromatosis-1 (NF1) were also overrepresented [16.7% MBCs compared to ∼5% of breast cancers overall; enrichment p = 0.028; mutation significance p = 0.006 (OncodriveFM)], consistent with published case reports implicating germline NF1 mutations in MBC risk. Taken together, we propose a practically minor but clinically significant modification to the guidelines: all WHO_1 mixed-type tumours should have the number of morphologies present recorded, as a mechanism for refining prognosis, and that EGFR and pan-cytokeratin expression are important prognostic markers. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing

Technological innovation and increased affordability have contributed to the widespread adoption of genome sequencing technologies in biomedical research. In particular large cancer research consortia have embraced next generation sequencing, and have used the technology to define the somatic mutation landscape of multiple cancer types. These studies have primarily utilised the Illumina HiSeq platforms. In this study we performed whole genome sequencing of three malignant pleural mesothelioma and matched normal samples using a new platform, the BGISEQ-500, and compared the results obtained with Illumina HiSeq X Ten. Germline and somatic, single nucleotide variants and small insertions or deletions were independently identified from data aligned human genome reference. The BGISEQ-500 and HiSeq X Ten platforms showed high concordance for germline calls with genotypes from SNP arrays (>99%). The germline and somatic single nucleotide variants identified in both sequencing platforms were highly concordant (86% and 72% respectively). These results indicate the potential applicability of the BGISEQ-500 platform for the identification of somatic and germline single nucleotide variants by whole genome sequencing. The BGISEQ-500 datasets described here represent the first publicly-available cancer genome sequencing performed using this platform.

Early clinical experiences with nintedanib in three UK tertiary interstitial lung disease centres

Background

Nintedanib has been shown to slow disease progression in patients with idiopathic pulmonary fibrosis (IPF). It was approved by the National Institute for Health and Care Excellence (NICE) in January 2016 for IPF patients with a forced vital capacity (FVC) of 50–80% in the United Kingdom (UK).

Aim

To report real world data about our early clinical experience using nintedanib in 187 patients with a multi-disciplinary (MDT) diagnosis of IPF in a manufacturer funded patient in need scheme (three UK centres) prior to NICE approval.

Methods

All patients with a MDT diagnosis of IPF from December 2014 to January 2016 commenced on nintedanib were included. Demographic details, adverse events (AEs) and where available lung function results were retrospectively collected from clinical letters.

Results

187 patients (76% males) with a median age of 72 years (49–89) were treated with nintedanib. The average pre-treatment FVC was 81.1 ± 19.8% and diffusion capacity of the lungs for carbon monoxide was 43.9 ± 15% (n = 82). Fifty percent of patients started nintedanib because they were ineligible for pirfenidone due to an FVC > 80%. The median treatment course was 8 ± 4 months. The majority of patients experienced 1–3 AEs with nintedanib (52%, n = 97). The most frequent AEs were diarrhoea (50%), nausea (36%), reduced appetite (24%), tiredness (20%) and gastro-oesophageal reflux (18%). The majority of AEs resulted in no change in treatment (64%, n = 461). 21% (n = 150) of AEs resulted in a dose reduction and 13% (n = 94) necessitated discontinuation of treatment. 1 in 5 patients discontinued treatment either temporarily or on a permanent basis during the monitoring period. In a select cohort of patients, a statistically significant greater proportion of patients remained stable or improved and a lower proportion declined, as depicted by FVC changes of > 5% after nintedanib commencement (P < 0.05 using Chi squared test).

Conclusions

Nintedanib is well tolerated and has an acceptable safety profile. Only 8% of those reporting diarrhoea discontinued treatment either on a temporary or permanent basis. There were no signals with respect to increased cardiovascular morbidity or major bleeding risk. This is in keeping with the INPULSIS clinical trial findings but in a real world cohort.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours

This corrects the article DOI: 10.1038/nature21063.

Contribution of plasma membrane lipid domains to red blood cell (re)shaping

Although lipid domains have been evidenced in several living cell plasma membranes, their roles remain largely unclear. We here investigated whether they could contribute to function-associated cell (re)shaping. To address this question, we used erythrocytes as cellular model since they (i) exhibit a specific biconcave shape, allowing for reversible deformation in blood circulation, which is lost by membrane vesiculation upon aging; and (ii) display at their outer plasma membrane leaflet two types of submicrometric domains differently enriched in cholesterol and sphingomyelin. We here reveal the specific association of cholesterol- and sphingomyelin-enriched domains with distinct curvature areas of the erythrocyte biconcave membrane. Upon erythrocyte deformation, cholesterol-enriched domains gathered in high curvature areas. In contrast, sphingomyelin-enriched domains increased in abundance upon calcium efflux during shape restoration. Upon erythrocyte storage at 4 °C (to mimick aging), lipid domains appeared as specific vesiculation sites. Altogether, our data indicate that lipid domains could contribute to erythrocyte function-associated (re)shaping.

Whole-genome landscapes of major melanoma subtypes

Melanoma of the skin is a common cancer only in Europeans, whereas it arises in internal body surfaces (mucosal sites) and on the hands and feet (acral sites) in people throughout the world. Here we report analysis of whole-genome sequences from cutaneous, acral and mucosal subtypes of melanoma. The heavily mutated landscape of coding and non-coding mutations in cutaneous melanoma resolved novel signatures of mutagenesis attributable to ultraviolet radiation. However, acral and mucosal melanomas were dominated by structural changes and mutation signatures of unknown aetiology, not previously identified in melanoma. The number of genes affected by recurrent mutations disrupting non-coding sequences was similar to that affected by recurrent mutations to coding sequences. Significantly mutated genes included BRAF, CDKN2A, NRAS and TP53 in cutaneous melanoma, BRAF, NRAS and NF1 in acral melanoma and SF3B1 in mucosal melanoma. Mutations affecting the TERT promoter were the most frequent of all; however, neither they nor ATRX mutations, which correlate with alternative telomere lengthening, were associated with greater telomere length. Most melanomas had potentially actionable mutations, most in components of the mitogen-activated protein kinase and phosphoinositol kinase pathways. The whole-genome mutation landscape of melanoma reveals diverse carcinogenic processes across its subtypes, some unrelated to sun exposure, and extends potential involvement of the non-coding genome in its pathogenesis.

Lost in translation: returning germline genetic results in genome-scale cancer research

BACKGROUND

The return of research results (RoR) remains a complex and well-debated issue. Despite the debate, actual data related to the experience of giving individual results back, and the impact these results may have on clinical care and health outcomes, is sorely lacking. Through the work of the Australian Pancreatic Cancer Genome Initiative (APGI) we: (1) delineate the pathway back to the patient where actionable research data were identified; and (2) report the clinical utilisation of individual results returned. Using this experience, we discuss barriers and opportunities associated with a comprehensive process of RoR in large-scale genomic research that may be useful for others developing their own policies.

METHODS

We performed whole-genome (n = 184) and exome (n = 208) sequencing of matched tumour-normal DNA pairs from 392 patients with sporadic pancreatic cancer (PC) as part of the APGI. We identified pathogenic germline mutations in candidate genes (n = 130) with established predisposition to PC or medium-high penetrance genes with well-defined cancer associated syndromes or phenotypes. Variants from candidate genes were annotated and classified according to international guidelines. Variants were considered actionable if clinical utility was established, with regard to prevention, diagnosis, prognostication and/or therapy.

RESULTS

A total of 48,904 germline variants were identified, with 2356 unique variants undergoing annotation and in silico classification. Twenty cases were deemed actionable and were returned via previously described RoR framework, representing an actionable finding rate of 5.1%. Overall, 1.78% of our cohort experienced clinical benefit from RoR.

CONCLUSION

Returning research results within the context of large-scale genomics research is a labour-intensive, highly variable, complex operation. Results that warrant action are not infrequent, but the prevalence of those who experience a clinical difference as a result of returning individual results is currently low.

Whole-genome landscape of pancreatic neuroendocrine tumours

The diagnosis of pancreatic neuroendocrine tumours (PanNETs) is increasing owing to more sensitive detection methods, and this increase is creating challenges for clinical management. We performed whole-genome sequencing of 102 primary PanNETs and defined the genomic events that characterize their pathogenesis. Here we describe the mutational signatures they harbour, including a deficiency in G:C > T:A base excision repair due to inactivation of MUTYH, which encodes a DNA glycosylase. Clinically sporadic PanNETs contain a larger-than-expected proportion of germline mutations, including previously unreported mutations in the DNA repair genes MUTYH, CHEK2 and BRCA2. Together with mutations in MEN1 and VHL, these mutations occur in 17% of patients. Somatic mutations, including point mutations and gene fusions, were commonly found in genes involved in four main pathways: chromatin remodelling, DNA damage repair, activation of mTOR signalling (including previously undescribed EWSR1 gene fusions), and telomere maintenance. In addition, our gene expression analyses identified a subgroup of tumours associated with hypoxia and HIF signalling.

Hypermutation In Pancreatic Cancer

Pancreatic cancer is molecularly diverse, with few effective therapies. Increased mutation burden and defective DNA repair are associated with response to immune checkpoint inhibitors in several other cancer types. We interrogated 385 pancreatic cancer genomes to define hypermutation and its causes. Mutational signatures inferring defects in DNA repair were enriched in those with the highest mutation burdens. Mismatch repair deficiency was identified in 1% of tumors harboring different mechanisms of somatic inactivation of MLH1 and MSH2. Defining mutation load in individual pancreatic cancers and the optimal assay for patient selection may inform clinical trial design for immunotherapy in pancreatic cancer.

Genomic analyses identify molecular subtypes of pancreatic cancer

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-β, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Lubricin/Proteoglycan 4 binds to and regulates the activity of Toll-Like Receptors In Vitro

Proteoglycan 4 (PRG4/lubricin) is secreted by cells that reside in articular cartilage and line the synovial joint. Lubricin may play a role in modulating inflammatory responses through interaction with CD44. This led us to examine if lubricin could be playing a larger role in the modulation of inflammation/immunity through interaction with Toll-like receptors (TLRs). Human Embryonic Kidney (HEK) cells overexpressing TLRs 2, 4 or 5 and surface plasmon resonance were employed to determine if full length recombinant human lubricin was able to bind to and activate TLRs. Primary human synovial fibroblasts were also examined using flow cytometry and Luminex multiplex ELISA. A rat destabilization model of osteoarthritis (OA) was used to determine if lubricin injections were able to regulate pain and/or inflammation in vivo. Lubricin can bind to and regulate the activity of TLRs, leading to downstream changes in inflammatory signalling independent of HA. We confirmed these findings in vivo through intra-articular injections of lubricin in a rat OA model where the inhibition of systemic inflammatory signaling and reduction in pain were observed. Lubricin plays an important role in regulating the inflammatory environment under both homeostatic and tissue injury states.