Clinically actionable mutation profiles in patients with cancer identified by whole-genome sequencing

Large-Scale Machine Learning Analysis Reveals DNA Methylation and Gene Expression Response Signatures for Gemcitabine-Treated Pancreatic Cancer

Background: Gemcitabine is a first-line chemotherapy for pancreatic adenocarcinoma (PAAD), but many PAAD patients do not respond to gemcitabine-containing treatments. Being able to predict such nonresponders would hence permit the undelayed administration of more promising treatments while sparing gemcitabine life-threatening side effects for those patients. Unfortunately, the few predictors of PAAD patient response to this drug are weak, none of them exploiting yet the power of machine learning (ML). Methods: Here, we applied ML to predict the response of PAAD patients to gemcitabine from the molecular profiles of their tumors. More concretely, we collected diverse molecular profiles of PAAD patient tumors along with the corresponding clinical data (gemcitabine responses and clinical features) from the Genomic Data Commons resource. From systematically combining 8 tumor profiles with 16 classification algorithms, each of the resulting 128 ML models was evaluated by multiple 10-fold cross-validations. Results: Only 7 of these 128 models were predictive, which underlines the importance of carrying out such a large-scale analysis to avoid missing the most predictive models. These were here random forest using 4 selected mRNAs [0.44 Matthews correlation coefficient (MCC), 0.785 receiver operating characteristic-area under the curve (ROC-AUC)] and XGBoost combining 12 DNA methylation probes (0.32 MCC, 0.697 ROC-AUC). By contrast, the hENT1 marker obtained much worse random-level performance (practically 0 MCC, 0.5 ROC-AUC). Despite not being trained to predict prognosis (overall and progression-free survival), these ML models were also able to anticipate this patient outcome. Conclusions: We release these promising ML models so that they can be evaluated prospectively on other gemcitabine-treated PAAD patients.

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases

BACKGROUND

Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome.

METHODS

We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants.

RESULTS

Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving.

CONCLUSIONS

Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Protocol for the Tessa Jowell BRAIN MATRIX Platform Study

INTRODUCTION

Gliomas are the most common primary tumour of the central nervous system (CNS), with an estimated annual incidence of 6.6 per 100 000 individuals in the USA and around 14 deaths per day from brain tumours in the UK. The genomic and biological landscape of brain tumours has been increasingly defined and, since 2016, the WHO classification of tumours of the CNS incorporates molecular data, along with morphology, to define tumour subtypes more accurately. The Tessa Jowell BRAIN MATRIX Platform (TJBM) study aims to create a transformative clinical research infrastructure that leverages UK National Health Service resources to support research that is patient centric and attractive to both academic and commercial investors.

METHODS AND ANALYSIS

The TJBM study is a programme of work with the principal purpose to improve the knowledge of glioma and treatment for patients with glioma. The programme includes a platform study and subsequent interventional clinical trials (as separate protocols). The platform study described here is the backbone data-repository of disease, treatment and outcome data from clinical, imaging and pathology data being collected in patients with glioma from secondary care hospitals. The primary outcome measure of the platform is time from biopsy to integrated histological-molecular diagnosis using whole-genome sequencing and epigenomic classification. Secondary outcome measures include those that are process centred, patient centred and framework based. Target recruitment for the study is 1000 patients with interim analyses at 100 and 500 patients.

ETHICS AND DISSEMINATION

The study will be performed in accordance with the recommendations guiding physicians in biomedical research involving human subjects, adopted by the 18th World Medical Association General Assembly, Helsinki, Finland and stated in the respective participating countries' laws governing human research, and Good Clinical Practice. The protocol was initially approved on 18 February 2020 by West Midlands - Edgbaston Research Ethics Committee; the current protocol (v3.0) was approved on 15 June 2022. Participants will be required to provide written informed consent. A meeting will be held after the end of the study to allow discussion of the main results among the collaborators prior to publication. The results of this study will be disseminated through national and international presentations and peer-reviewed publications. Manuscripts will be prepared by the Study Management Group and authorship will be determined by mutual agreement.

TRIAL REGISTRATION NUMBER

NCT04274283, 18-Feb-2020; ISRCTN14218060, 03-Feb-2020.

Clinical interpretation of whole-genome and whole-transcriptome sequencing for precision oncology

Whole-genome sequencing either alone or in combination with whole-transcriptome sequencing has started to be used to analyze clinical tumor samples to improve diagnosis, provide risk stratification, and select patient-specific therapies. Compared with current genomic testing strategies, largely focused on small number of genes tested individually or targeted panels, whole-genome and transcriptome sequencing (WGTS) provides novel opportunities to identify and report a potentially much larger number of actionable alterations with diagnostic, prognostic, and/or predictive impact. Such alterations include point mutations, indels, copy- number aberrations and structural variants, but also germline variants, fusion genes, noncoding alterations and mutational signatures. Nevertheless, these comprehensive tests are accompanied by many challenges ranging from the extent and diversity of sequence alterations detected by these methods to the complexity and limited existing standardization in interpreting them. We describe the challenges of WGTS interpretation and the opportunities with comprehensive genomic testing.

Branching clonal evolution patterns predominate mutational landscape in multiple myeloma

Multiple Myeloma (MM) arises from malignant transformation and deregulated proliferation of clonal plasma cells (PCs) harbouring heterogeneous molecular anomalies. The effect of evolving mutations on clone fitness and their cellular prevalence shapes the progressing myeloma genome and impacts clinical outcomes. Although clonal heterogeneity in MM is well established, which subclonal mutations emerge/persist/perish with progression in MM and which of these can be targeted therapeutically remains an open question. In line with this, we have sequenced pairwise whole exomes of 62 MM patients collected at two time points, i.e., at diagnosis and on progression. Somatic variants were called using a novel ensemble approach where a consensus was deduced from four variant callers (Illumina's Dragen, Strelka2, SomaticSniper and SpeedSeq) and actionable/druggable gene targets were identified. A marked intraclonal heterogeneity was observed. Branching evolution was observed among 72.58% patients, of whom 64.51% had low TMBs (<10) and 61.29% had 2 or more founder clones. The hypermutator patients (with high TMB levels ≥10 to ≤100) showed a significant decrease in their TMBs from diagnosis (median TMB 77.11) to progression (median TMB 31.22). A distinct temporal fall in subclonal driver mutations was identified recurrently across diagnosis to progression e.g., in PABPC1, BRAF, KRAS, CR1, DIS3 and ATM genes in 3 or more patients suggesting such patients could be treated early with target specific drugs like Vemurafenib/Cobimetinib. An analogous rise in driver mutations was observed in KMT2C, FOXD4L1, SP140, NRAS and other genes. A few drivers such as FAT4, IGLL5 and CDKN1A retained consistent distribution patterns at two time points. These findings are clinically relevant and point at consideration of evaluating multi time point subclonal mutational landscapes for designing better risk stratification strategies and tailoring time to time risk adapted combination therapies in future.

Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia

The transformation of chronic lymphocytic leukemia (CLL) to high-grade B-cell lymphoma is known as Richter syndrome (RS), a rare event with dismal prognosis. In this study, we conducted whole-genome sequencing (WGS) of paired circulating CLL (PB-CLL) and RS biopsies (tissue-RS) from 17 patients recruited into a clinical trial (CHOP-O). We found that tissue-RS was enriched for mutations in poor-risk CLL drivers and genes in the DNA damage response (DDR) pathway. In addition, we identified genomic aberrations not previously implicated in RS, including the protein tyrosine phosphatase receptor (PTPRD) and tumor necrosis factor receptor-associated factor 3 (TRAF3). In the noncoding genome, we discovered activation-induced cytidine deaminase-related and unrelated kataegis in tissue-RS affecting regulatory regions of key immune-regulatory genes. These include BTG2, CXCR4, NFATC1, PAX5, NOTCH-1, SLC44A5, FCRL3, SELL, TNIP2, and TRIM13. Furthermore, differences between the global mutation signatures of pairs of PB-CLL and tissue-RS samples implicate DDR as the dominant mechanism driving transformation. Pathway-based clonal deconvolution analysis showed that genes in the MAPK and DDR pathways demonstrate high clonal-expansion probability. Direct comparison of nodal-CLL and tissue-RS pairs from an independent cohort confirmed differential expression of the same pathways by RNA expression profiling. Our integrated analysis of WGS and RNA expression data significantly extends previous targeted approaches, which were limited by the lack of germline samples, and it facilitates the identification of novel genomic correlates implicated in RS transformation, which could be targeted therapeutically. Our results inform the future selection of investigative agents for a UK clinical platform study. This trial was registered at www.clinicaltrials.gov as #NCT03899337.

Short and long-read genome sequencing methodologies for somatic variant detection; genomic analysis of a patient with diffuse large B-cell lymphoma

Recent advances in throughput and accuracy mean that the Oxford Nanopore Technologies PromethION platform is a now a viable solution for genome sequencing. Much of the validation of bioinformatic tools for this long-read data has focussed on calling germline variants (including structural variants). Somatic variants are outnumbered many-fold by germline variants and their detection is further complicated by the effects of tumour purity/subclonality. Here, we evaluate the extent to which Nanopore sequencing enables detection and analysis of somatic variation. We do this through sequencing tumour and germline genomes for a patient with diffuse B-cell lymphoma and comparing results with 150 bp short-read sequencing of the same samples. Calling germline single nucleotide variants (SNVs) from specific chromosomes of the long-read data achieved good specificity and sensitivity. However, results of somatic SNV calling highlight the need for the development of specialised joint calling algorithms. We find the comparative genome-wide performance of different tools varies significantly between structural variant types, and suggest long reads are especially advantageous for calling large somatic deletions and duplications. Finally, we highlight the utility of long reads for phasing clinically relevant variants, confirming that a somatic 1.6 Mb deletion and a p.(Arg249Met) mutation involving TP53 are oriented in trans.

Tumor Profiling at the Service of Cancer Therapy

Cancer treatment options have evolved significantly in the past few years. From the initial surgical procedures, to the latest next-generation technologies, we are now in the position to analyze and understand tumors in a one-by-one basis and use that to our advantage to provide with individualized treatment options that may increase patient survival. In this review, we will focus on how tumor profiling has evolved over the past decades to deliver more efficient and personalized treatment options, and how novel technologies can help us envisage the future of precision oncology toward a better management and, ultimately, increased survival.

A Novel Six-Gene Signature for Prognosis Prediction in Ovarian Cancer

Ovarian cancer (OC) is the most malignant tumor in the female reproductive tract. Although abundant molecular biomarkers have been identified, a robust and accurate gene expression signature is still essential to assist oncologists in evaluating the prognosis of OC patients. In this study, samples from 367 patients in The Cancer Genome Atlas (TCGA) database were subjected to mRNA expression profiling. Then, we used a gene set enrichment analysis (GSEA) to screen genes correlated with epithelial–mesenchymal transition (EMT) and assess their prognostic power with a Cox proportional regression model. Six genes (TGFBI, SFRP1, COL16A1, THY1, PPIB, BGN) associated with overall survival (OS) were used to construct a risk assessment model, after which the patients were divided into high-risk and low-risk groups. The six-gene signature was an independent prognostic biomarker of OS for OC patients based on the multivariate Cox regression analysis. In addition, the six-gene model was validated with samples from the Gene Expression Omnibus (GEO) database. In summary, we established a six-gene signature relevant to the prognosis of OC, which might become a therapeutic tool with clinical applications in the future.

A comprehensive custom panel evaluation for routine hereditary cancer testing: improving the yield of germline mutation detection

Background

In the context of our Regional Program of Hereditary Cancer, individuals fulfilling the criteria are tested for germline mutations to subsequently establish the clinical management. Our standard diagnostic approach focuses on sequencing a few classic high-risk genes, a method that frequently renders uninformative genetic results. This study aims to examine the improved yield offered by an On-Demand panel.

Methods

We designed an On-Demand panel for the analysis of 35-genes associated with inherited cancer susceptibility in a total of 128 cases of Hereditary Breast and Ovarian Cancer (HBOC) and Hereditary Nonpolyposis Colorectal Cancer (HNPCC).

Results

Eighteen deleterious mutations were detected, in both routinely (BRCA2, MLH1, MSH2, PMS2) and non-routinely (ATM, BLM, BRIP1, CHEK2, MUTYH) tested genes. The screening extended to 35 genes rendered by patients carrying several- up to 6-Variants of Unknown Significance (VUS). Moreover, we confirmed the splicing disruption at RNA level for a not previously reported BRIP1 splicing mutation. Using an On-Demand panel, we identified 18 pathogenic mutation carriers, seven of which would have gone unnoticed with traditional analysis.

Conclusions

Our results reinforce the utility of NGS gene panels in the diagnostic routine to increase the performance of genetic testing, especially in individuals from families with overlapping cancer phenotypes.

Morphologic and Genomic Heterogeneity in the Evolution and Progression of Breast Cancer

: Breast cancer is a remarkably complex and diverse disease. Subtyping based on morphology, genomics, biomarkers and/or clinical parameters seeks to stratify optimal approaches for management, but it is clear that every breast cancer is fundamentally unique. Intra-tumour heterogeneity adds further complexity and impacts a patient's response to neoadjuvant or adjuvant therapy. Here, we review some established and more recent evidence related to the complex nature of breast cancer evolution. We describe morphologic and genomic diversity as it arises spontaneously during the early stages of tumour evolution, and also in the context of treatment where the changing subclonal architecture of a tumour is driven by the inherent adaptability of tumour cells to evolve and resist the selective pressures of therapy.

Genomics-Enabled Precision Medicine for Cancer

Genomic information is increasingly being incorporated into clinical cancer care. Large-scale sequencing efforts have deepened our understanding of the genomic landscape of cancer and contributed to the expanding catalog of alterations being leveraged to aid in cancer diagnosis, prognosis, and treatment. Genomic profiling can provide clinically relevant information regarding somatic point mutations, copy number alterations, translocations, and gene fusions. Genomic features, such as mutational burden, can also be measured by more comprehensive sequencing strategies and have shown value in informing potential treatment options. Ongoing clinical trials are evaluating the use of molecularly targeted agents in genomically defined subsets of cancers within and across tumor histologies. Continued advancements in clinical genomics promise to further expand the application of genomics-enabled medicine to a broader spectrum of oncology patients.