FiNGS: high quality somatic mutations using filters for next generation sequencing

CLEMENT: genomic decomposition and reconstruction of non-tumor subclones

Genome-level clonal decomposition of a single specimen has been widely studied; however, it is mostly limited to cancer research. In this study, we developed a new algorithm CLEMENT, which conducts accurate decomposition and reconstruction of multiple subclones in genome sequencing of non-tumor (normal) samples. CLEMENT employs the Expectation-Maximization (EM) algorithm with optimization strategies specific to non-tumor subclones, including false variant call identification, non-disparate clone fuzzy clustering, and clonal allele fraction confinement. In the simulation and in vitro cell line mixture data, CLEMENT outperformed current cancer decomposition algorithms in estimating the number of clones (root-mean-square-error = 0.58-0.78 versus 1.43-3.34) and in the variant-clone membership agreement (∼85.5% versus 70.1-76.7%). Additional testing on human multi-clonal normal tissue sequencing confirmed the accurate identification of subclones that originated from different cell types. Clone-level analysis, including mutational burden and signatures, provided a new understanding of normal-tissue composition. We expect that CLEMENT will serve as a crucial tool in the currently emerging field of non-tumor genome analysis.

Analysis of somatic mutations in whole blood from 200,618 individuals identifies pervasive positive selection and novel drivers of clonal hematopoiesis

Human aging is marked by the emergence of a tapestry of clonal expansions in dividing tissues, particularly evident in blood as clonal hematopoiesis (CH). CH, linked to cancer risk and aging-related phenotypes, often stems from somatic mutations in a set of established genes. However, the majority of clones lack known drivers. Here we infer gene-level positive selection in whole blood exomes from 200,618 individuals in UK Biobank. We identify 17 additional genes, ZBTB33, ZNF318, ZNF234, SPRED2, SH2B3, SRCAP, SIK3, SRSF1, CHEK2, CCDC115, CCL22, BAX, YLPM1, MYD88, MTA2, MAGEC3 and IGLL5, under positive selection at a population level, and validate this selection pattern in 10,837 whole genomes from single-cell-derived hematopoietic colonies. Clones with mutations in these genes grow in frequency and size with age, comparable to classical CH drivers. They correlate with heightened risk of infection, death and hematological malignancy, highlighting the significance of these additional genes in the aging process.

Validation of genetic variants from NGS data using deep convolutional neural networks

Accurate somatic variant calling from next-generation sequencing data is one most important tasks in personalised cancer therapy. The sophistication of the available technologies is ever-increasing, yet, manual candidate refinement is still a necessary step in state-of-the-art processing pipelines. This limits reproducibility and introduces a bottleneck with respect to scalability. We demonstrate that the validation of genetic variants can be improved using a machine learning approach resting on a Convolutional Neural Network, trained using existing human annotation. In contrast to existing approaches, we introduce a way in which contextual data from sequencing tracks can be included into the automated assessment. A rigorous evaluation shows that the resulting model is robust and performs on par with trained researchers following published standard operating procedure.

PD-L1 blockade in combination with carboplatin as immune induction in metastatic lobular breast cancer: the GELATO trial

Invasive lobular breast cancer (ILC) is the second most common histological breast cancer subtype, but ILC-specific trials are lacking. Translational research revealed an immune-related ILC subset, and in mouse ILC models, synergy between immune checkpoint blockade and platinum was observed. In the phase II GELATO trial ( NCT03147040 ), patients with metastatic ILC were treated with weekly carboplatin (area under the curve 1.5 mg ml-1 min-1) as immune induction for 12 weeks and atezolizumab (PD-L1 blockade; triweekly) from the third week until progression. Four of 23 evaluable patients had a partial response (17%), and 2 had stable disease, resulting in a clinical benefit rate of 26%. From these six patients, four had triple-negative ILC (TN-ILC). We observed higher CD8+ T cell infiltration, immune checkpoint expression and exhausted T cells after treatment. With this GELATO trial, we show that ILC-specific clinical trials are feasible and demonstrate promising antitumor activity of atezolizumab with carboplatin, particularly for TN-ILC, and provide insights for the design of highly needed ILC-specific trials.

A Functional Precision Medicine Pipeline Combines Comparative Transcriptomics and Tumor Organoid Modeling to Identify Bespoke Treatment Strategies for Glioblastoma

Li Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome caused by germline mutations in TP53. TP53 is the most common mutated gene in human cancer, occurring in 30-50% of glioblastomas (GBM). Here, we highlight a precision medicine platform to identify potential targets for a GBM patient with LFS. We used a comparative transcriptomics approach to identify genes that are uniquely overexpressed in the LFS GBM patient relative to a cancer compendium of 12,747 tumor RNA sequencing data sets, including 200 GBMs. STAT1 and STAT2 were identified as being significantly overexpressed in the LFS patient, indicating ruxolitinib, a Janus kinase 1 and 2 inhibitors, as a potential therapy. The LFS patient had the highest level of STAT1 and STAT2 expression in an institutional high-grade glioma cohort of 45 patients, further supporting the cancer compendium results. To empirically validate the comparative transcriptomics pipeline, we used a combination of adherent and organoid cell culture techniques, including ex vivo patient-derived organoids (PDOs) from four patient-derived cell lines, including the LFS patient. STAT1 and STAT2 expression levels in the four patient-derived cells correlated with levels identified in the respective parent tumors. In both adherent and organoid cultures, cells from the LFS patient were among the most sensitive to ruxolitinib compared to patient-derived cells with lower STAT1 and STAT2 expression levels. A spheroid-based drug screening assay (3D-PREDICT) was performed and used to identify further therapeutic targets. Two targeted therapies were selected for the patient of interest and resulted in radiographic disease stability. This manuscript supports the use of comparative transcriptomics to identify personalized therapeutic targets in a functional precision medicine platform for malignant brain tumors.

Genomic and Transcriptomic Profiling of Brain Metastases

Simple Summary

Brain metastases (BM) are the most common brain tumors in adults and are the main cause of cancer-associated death. Omics analysis of BM will allow for a better understanding of metastatic progression, prognosis and therapeutic targeting. In this study, BM samples underwent comprehensive molecular profiling with genomics and transcriptomics. Mutational signatures suggested that most mutations were gained prior to metastasis. A novel copy number event centered around the MCL1 gene was found in 75% of all samples. Transcriptomics revealed that melanoma BM formed a distinct cluster in comparison to other subtypes. Poor survival correlated to self-identified black race and absence of radiation treatment but not molecular profiles. These data identify potential new drivers of brain metastatic progression, implicate that melanoma BM are distinctive and likely responsive to unique therapies, and further investigation of sociodemographic and clinical features are needed in BM cohorts.

Abstract

Brain metastases (BM) are the most common brain tumors in adults occurring in up to 40% of all cancer patients. Multi-omics approaches allow for understanding molecular mechanisms and identification of markers with prognostic significance. In this study, we profile 130 BM using genomics and transcriptomics and correlate molecular characteristics to clinical parameters. The most common tumor origins for BM were lung (40%) followed by melanoma (21%) and breast (15%). Melanoma and lung BMs contained more deleterious mutations than other subtypes (p < 0.001). Mutational signatures suggested that the bulk of the mutations were gained before metastasis. A novel copy number event centered around the MCL1 gene was found in 75% of all samples, suggesting a broader role in promoting metastasis. Unsupervised hierarchical cluster analysis of transcriptional signatures available in 65 samples based on the hallmarks of cancer revealed four distinct clusters. Melanoma samples formed a distinctive cluster in comparison to other BM subtypes. Characteristics of molecular profiles did not correlate with survival. However, patients with self-identified black race or those who did not receive radiation correlated with poor survival. These data identify potential new drivers of brain metastatic progression. Our data also suggest further investigation of sociodemographic and clinical features is needed in BM cohorts.