Personalized genomic analyses for cancer mutation discovery and interpretation

Spontaneous-stimulated Raman co-localization dual-modal analysis approach for efficient identification of tumor cells

The study of highly heterogeneous tumor cells, especially acute myeloid leukemia (AML) cells, usually relies on invasive analytical methods such as morphology, immunology, cytogenetics, and molecular biology classification, which are complex and time-consuming to perform. Mortality is high if patients are not diagnosed in a timely manner, so rapid label-free analysis of gene expression and metabolites within single-cell substructures is extremely important for clinical diagnosis and treatment. As a label-free and non-destructive vibrational detection technique, spontaneous Raman scattering provides molecular information across the full spectrum of the cell but lacks rapid imaging localization capabilities. In contrast, stimulated Raman scattering (SRS) provides a high-speed, high-resolution imaging view that can offer real-time subcellular localization assistance for spontaneous Raman spectroscopic detection. In this paper, we combined multi-color SRS microscopy with spontaneous Raman to develop a co-localized Raman imaging and spectral detection system (CRIS) for high-speed chemical imaging and quantitative spectral analysis of subcellular structures. Combined with multivariate statistical analysis methods, CRIS efficiently differentiated AML from normal leukocytes with an accuracy of 98.1 % and revealed the differences in the composition of nuclei and cytoplasm of AML relative to normal leukocytes. Compared to conventional Raman spectroscopy blind sampling without imaging localization, CRIS increased the efficiency of single-cell detection by at least three times. In addition, using the same approach for further identification of AML subtypes M2 and M3, we demonstrated that intracytoplasmic differential expression of proteins is a marker for their rapid and accurate classifying. CRIS analysis methods are expected to pave the way for clinical translation of rapid tumor cell identification.

Pharmacogenomics: A Genetic Approach to Drug Development and Therapy

The majority of the well-known pharmacogenomics research used in the medical sciences contributes to our understanding of medication interactions. It has a significant impact on treatment and drug development. The broad use of pharmacogenomics is required for the progress of therapy. The main focus is on how genes and an intricate gene system affect the body's reaction to medications. Novel biomarkers that help identify a patient group that is more or less likely to respond to a certain medication have been discovered as a result of recent developments in the field of clinical therapeutics. It aims to improve customized therapy by giving the appropriate drug at the right dose at the right time and making sure that the right prescriptions are issued. A combination of genetic, environmental, and patient variables that impact the pharmacokinetics and/or pharmacodynamics of medications results in interindividual variance in drug response. Drug development, illness susceptibility, and treatment efficacy are all impacted by pharmacogenomics. The purpose of this work is to give a review that might serve as a foundation for the creation of new pharmacogenomics applications, techniques, or strategies.

Measurement of ctDNA Tumor Fraction Identifies Informative Negative Liquid Biopsy Results and Informs Value of Tissue Confirmation

PURPOSE

Liquid biopsy (LBx) for tumor profiling is increasingly used, but concerns remain regarding negative results. A lack of results may truly reflect tumor genomics, or it may be a false negative that would be clarified by tissue testing. A method of distinguishing between these scenarios could help clarify when follow-on tissue testing is valuable.

EXPERIMENTAL DESIGN

Here we evaluate circulating tumor DNA (ctDNA) tumor fraction (TF), a quantification of ctDNA in LBx samples, for utility in identifying true negative results. We assessed concordance between LBx and tissue-based results, stratified by ctDNA TF, in a real-world genomic dataset of paired samples across multiple disease types. We also evaluated the frequency of tissue results identifying driver alterations in patients with lung cancer after negative LBx in a real-world clinicogenomic database.

RESULTS

The positive percent agreement and negative predictive value between liquid and tissue samples for driver alterations increased from 63% and 66% for all samples to 98% and 97% in samples with ctDNA TF ≥1%. Among 505 patients with lung cancer with no targetable driver alterations found by LBx who had subsequent tissue-based profiling, 37% had a driver, all of which had ctDNA TF <1%.

CONCLUSIONS

Patients with lung cancer with negative LBx and ctDNA TF ≥1% are unlikely to have a driver detected on confirmatory tissue testing; such informative negative results may benefit instead from prompt treatment initiation. Conversely, negative LBx with ctDNA TF <1% will commonly have a driver identified by follow-up tissue testing and should be prioritized for reflex testing.

Bridging horizons beyond CIRCULATE-Japan: a new paradigm in molecular residual disease detection via whole genome sequencing-based circulating tumor DNA assay

Circulating tumor DNA (ctDNA) is the fraction of cell-free DNA in patient blood that originates from a tumor. Advances in DNA sequencing technologies and our understanding of the molecular biology of tumors have increased interest in exploiting ctDNA to facilitate detection of molecular residual disease (MRD). Analysis of ctDNA as a promising MRD biomarker of solid malignancies has a central role in precision medicine initiatives exemplified by our CIRCULATE-Japan project involving patients with resectable colorectal cancer. Notably, the project underscores the prognostic significance of the ctDNA status at 4 weeks post-surgery and its correlation to adjuvant therapy efficacy at interim analysis. This substantiates the hypothesis that MRD is a critical prognostic indicator of relapse in patients with colorectal cancer. Despite remarkable advancements, challenges endure, primarily attributable to the exceedingly low ctDNA concentration in peripheral blood, particularly in scenarios involving low tumor shedding and the intrinsic error rates of current sequencing technologies. These complications necessitate more sensitive and sophisticated assays to verify the clinical utility of MRD across all solid tumors. Whole genome sequencing (WGS)-based tumor-informed MRD assays have recently demonstrated the ability to detect ctDNA in the parts-per-million range. This review delineates the current landscape of MRD assays, highlighting WGS-based approaches as the forefront technique in ctDNA analysis. Additionally, it introduces our upcoming endeavor, WGS-based pan-cancer MRD detection via ctDNA, in our forthcoming project, SCRUM-Japan MONSTAR-SCREEN-3.

Capture-based targeted sequencing using a T-cell control in myeloid malignancies and idiopathic cytopenias

We report on a study of next-generation sequencing in 257 patients undergoing investigations for cytopenias. We sequenced bone marrow aspirates using a target enrichment panel comprising 82 genes and used T cells from paired blood as a control. One hundred and sixty patients had idiopathic cytopenias, 81 had myeloid malignancies and 16 had lymphoid malignancies or other diagnoses. Forty-seven of the 160 patients with idiopathic cytopenias had evidence of somatic pathogenic variants consistent with clonal cytopenias. Only 39 genes of the 82 tested were mutated in the 241 patients with either idiopathic cytopenias or myeloid neoplasms. We confirm that T cells can be used as a control to distinguish between germline and somatic variants. The use of paired analysis with a T-cell control significantly reduced the time molecular scientists spent reporting compared to unpaired analysis. We identified somatic variants of uncertain significance (VUS) in a higher proportion (24%) of patients with myeloid malignancies or clonal cytopenias compared to less than 2% of patients with non-clonal cytopenias. This suggests that somatic VUS are indicators of a clonal process. Lastly, we show that blood depleted of lymphocytes can be used in place of bone marrow as a source of material for sequencing.

Whole-exome sequencing of Nigerian benign prostatic hyperplasia reveals increased alterations in apoptotic pathways

BACKGROUND

Through whole-exome sequencing of 60 formalin-fixed paraffin-embedded Nigerian (NGRn) benign prostatic hyperplasia (BPH) samples, we identified germline and somatic alterations in apoptotic pathways impacting BPH development and progression. Prostate enlargement is a common occurrence in male aging; however, this enlargement can lead to lower urinary tract symptoms that negatively impact quality of life. This impact is disproportionately present in men of African ancestry. BPH pathophysiology is poorly understood and studies examining non-European populations are lacking.

METHODS

In this study, NGRn BPH, normal prostate, and prostate cancer (PCa) tumor samples were sequenced and compared to characterize genetic alterations in NGRn BPH.

RESULTS

Two hundred and two nonbenign, ClinVar-annotated germline variants were present in NGRn BPH samples. Six genes [BRCA1 (92%), HSD3B1 (85%), TP53 (37%), PMS2 (23%), BARD1 (20%), and BRCA2 (17%)] were altered in at least 10% of samples; however, compared to NGRn normal and tumor, the frequency of alterations in BPH samples showed no significant differences at the gene or variant level. BRCA2_rs11571831 and TP53_rs1042522 germline alterations had a statistically significant co-occurrence interaction in BPH samples. In at least two BPH samples, 173 genes harbored somatic variants known to be clinically actionable. Three genes (COL18A1, KIF16B, and LRP1) showed a statistically significant (p < 0.05) higher frequency in BPH. NGRn BPH also had five gene pairs (PKD1/KIAA0100, PKHD1/PKD1, DNAH9/LRP1B, NWD1/DCHS2, and TCERG1/LMTK2) with statistically significant co-occurring interactions. Two hundred and seventy-nine genes contained novel somatic variants in NGRn BPH. Three genes (CABP1, FKBP1C, and RP11-595B24.2) had a statistically significant (p < 0.05) higher alteration frequency in NGRn BPH and three were significantly higher in NGRn tumor (CACNA1A, DMKN, and CACNA2D2). Pairwise Fisher's exact tests showed 14 gene pairs with statistically significant (p < 0.05) interactions and four interactions approaching significance (p < 0.10). Mutational patterns in NGRn BPH were similar to COSMIC (Catalog of Somatic Mutations in Cancer) signatures associated with aging and dysfunctional DNA damage repair.

CONCLUSIONS

NGRn BPH contained significant germline alteration interactions (BRCA2_rs11571831 and TP53_rs1042522) and increased somatic alteration frequencies (LMTK2, LRP1, COL18A1, CABP1, and FKBP1C) that impact apoptosis. Normal prostate development is maintained by balancing apoptotic and proliferative activity. Dysfunction in either mechanism can lead to abnormal prostate growth. This work is the first to examine genomic sequencing in NGRn BPH and provides data that fill known gaps in the understanding BPH and how it impacts men of African ancestry.

Limited Independent Follow-Up with Germline Testing of Variants Detected in BRCA1 and BRCA2 by Tumor-Only Sequencing

Introduction

Genomic profiling is performed in patients with advanced or metastatic cancer, in order to direct cancer treatment, often sequencing tumor-only, without a matched germline comparator. However, because many of the genes analyzed on tumor profiling overlap with those known to be associated with hereditary cancer predisposition syndromes (HCPS), tumor-only profiling can unknowingly uncover germline pathogenic (P) and likely pathogenic variants (LPV). In this study, we evaluated the number of patients with P/LPVs identified in BRCA1 and BRCA2 (BRCA1/2) via tumor-only profiling, then determined the germline testing outcomes for those patients.

Methods

A retrospective chart review was performed to identify patients with BRCA1/2 variants on tumor-only genomic profiling, and whether they had germline testing.

Results

This study found that of 2923 patients with 36 tumor types who underwent tumor-only testing, 554 had a variant in BRCA1/2 (19.0%); 119 of the 554 patients (21.5%) had a P/LP BRCA1/2 variant, representing 4.1% of the overall population who underwent genomic profiling. Seventy-three (61.3%) of 119 patients with BRCA1/2 P/LPV on tumor-only testing did not undergo germline testing, 34 (28.6%) had already had germline testing before tumor-only testing, and 12 (10.1%) underwent germline testing after tumor-only testing. Twenty-eight germline BRCA1/2 P/LPVs were detected, 24 in those who had prior germline testing, and 4 among the 12 patients who had germline testing after tumor-only testing.

Conclusion

Tumor-only testing is likely to identify P/LPVs in BRCA1/2. Efforts to improve follow-up germline testing is needed to improve identification of germline BRCA1/2 alterations.

Impact of genetic counseling strategy on diagnostic yield and workload for genome-sequencing-based tumor diagnostics

PURPOSE

Genome sequencing (GS) enables comprehensive molecular analysis of tumors and identification of hereditary cancer predisposition. According to guidelines, directly determining pathogenic germline variants (PGVs) requires pretest genetic counseling, which is cost-ineffective. Referral for genetic counseling based on tumor variants alone could miss relevant PGVs and/or result in unnecessary referrals.

METHODS

We validated GS for detection of germline variants and simulated 3 strategies using paired tumor-normal GS data of 937 metastatic patients. In strategy-1, genetic counseling before tumor testing allowed direct PGV analysis. In strategy-2 and -3, germline testing and referral for post-test genetic counseling is based on tumor variants using Dutch (strategy-2) or Europen Society for Medical Oncology (ESMO) Precision Medicine Working Group (strategy-3) guidelines.

RESULTS

In strategy-1, PGVs would be detected in 50 patients (number-needed-to counsel; NTC = 18.7). In strategy-2, 86 patients would have been referred for genetic counseling and 43 would have PGVs (NTC = 2). In strategy-3, 94 patients would have been referred for genetic counseling and 32 would have PGVs (NTC = 2.9). Hence, 43 and 62 patients, respectively, were unnecessarily referred based on a somatic variant.

CONCLUSION

Both post-tumor test counseling strategies (2 and 3) had significantly lower NTC, and strategy-2 had the highest PGV yield. Combining pre-tumor test mainstreaming and post-tumor test counseling may maximize the clinically relevant PGV yield and minimize unnecessary referrals.

Divergent immune microenvironments in two tumor nodules from a patient with mismatch repair-deficient prostate cancer

Patients with prostate cancer (PC) generally do not respond favorably to immune checkpoint inhibitors, which may be due to a low abundance of tumor-infiltrating lymphocytes even when mutational load is high. Here, we identified a patient who presented with high-grade primary prostate cancer with two adjacent tumor nodules. While both nodules were mismatch repair-deficient (MMRd), exhibited pathogenic MSH2 and MSH6 alterations, had a high tumor mutational burden (TMB), and demonstrated high microsatellite instability (MSI), they had markedly distinct immune phenotypes. The first displayed a dense infiltrate of lymphocytes ("hot nodule"), while the second displayed significantly fewer infiltrating lymphocytes ("cold nodule"). Whole-exome DNA analysis found that both nodules shared many identical mutations, indicating that they were derived from a single clone. However, the cold nodule appeared to be sub-clonal relative to the hot nodule, suggesting divergent evolution of the cold nodule from the hot nodule. Whole-transcriptome RNA analysis found that the cold nodule demonstrated lower expression of genes related to antigen presentation (HLA) and, paradoxically, classical tumor immune tolerance markers such as PD-L1 (CD274) and CTLA-4. Immune cell deconvolution suggested that the hot nodule was enriched not only in CD8+ and CD4 + T lymphocytes, but also in M1 macrophages, activated NK cells, and γδ T cells compared to the cold nodule. This case highlights that MMRd/TMB-high PC can evolve to minimize an anti-tumor immune response, and nominates downregulation of antigen presentation machinery (HLA loss) as a potential mechanism of adaptive immune evasion in PC.

Elucidating the Heterogeneity of Immunotherapy Response and Immune-Related Toxicities by Longitudinal ctDNA and Immune Cell Compartment Tracking in Lung Cancer

PURPOSE

Although immunotherapy is the mainstay of therapy for advanced non-small cell lung cancer (NSCLC), robust biomarkers of clinical response are lacking. The heterogeneity of clinical responses together with the limited value of radiographic response assessments to timely and accurately predict therapeutic effect-especially in the setting of stable disease-calls for the development of molecularly informed real-time minimally invasive approaches. In addition to capturing tumor regression, liquid biopsies may be informative in capturing immune-related adverse events (irAE).

EXPERIMENTAL DESIGN

We investigated longitudinal changes in circulating tumor DNA (ctDNA) in patients with metastatic NSCLC who received immunotherapy-based regimens. Using ctDNA targeted error-correction sequencing together with matched sequencing of white blood cells and tumor tissue, we tracked serial changes in cell-free tumor load (cfTL) and determined molecular response. Peripheral T-cell repertoire dynamics were serially assessed and evaluated together with plasma protein expression profiles.

RESULTS

Molecular response, defined as complete clearance of cfTL, was significantly associated with progression-free (log-rank P = 0.0003) and overall survival (log-rank P = 0.01) and was particularly informative in capturing differential survival outcomes among patients with radiographically stable disease. For patients who developed irAEs, on-treatment peripheral blood T-cell repertoire reshaping, assessed by significant T-cell receptor (TCR) clonotypic expansions and regressions, was identified on average 5 months prior to clinical diagnosis of an irAE.

CONCLUSIONS

Molecular responses assist with the interpretation of heterogeneous clinical responses, especially for patients with stable disease. Our complementary assessment of the peripheral tumor and immune compartments provides an approach for monitoring of clinical benefits and irAEs during immunotherapy.

Germline POT1 Variants: A Critical Perspective on POT1 Tumor Predisposition Syndrome

The Protection of Telomere 1 (POT1) gene was identified as a melanoma predisposition candidate nearly 10 years ago. Thereafter, various cancers have been proposed as associated with germline POT1 variants in the context of the so-called POT1 Predisposition Tumor Syndrome (POT1-TPD). While the key role, and related risks, of the alterations in POT1 in melanoma are established, the correlation between germline POT1 variants and the susceptibility to other cancers partially lacks evidence, due also to the rarity of POT1-TPD. Issues range from the absence of functional or segregation studies to biased datasets or the need for a revised classification of variants. Furthermore, a proposal of a surveillance protocol related to the cancers associated with POT1 pathogenic variants requires reliable data to avoid an excessive, possibly unjustified, burden for POT1 variant carriers. We propose a critical perspective regarding data published over the last 10 years that correlate POT1 variants to various types of cancer, other than cutaneous melanoma, to offer food for thought for the specialists who manage cancer predisposition syndromes and to stimulate a debate on the grey areas that have been exposed.

NCT/DKFZ MASTER handbook of interpreting whole-genome, transcriptome, and methylome data for precision oncology

Analysis of selected cancer genes has become an important tool in precision oncology but cannot fully capture the molecular features and, most importantly, vulnerabilities of individual tumors. Observational and interventional studies have shown that decision-making based on comprehensive molecular characterization adds significant clinical value. However, the complexity and heterogeneity of the resulting data are major challenges for disciplines involved in interpretation and recommendations for individualized care, and limited information exists on how to approach multilayered tumor profiles in clinical routine. We report our experience with the practical use of data from whole-genome or exome and RNA sequencing and DNA methylation profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) program of the National Center for Tumor Diseases (NCT) Heidelberg and Dresden and the German Cancer Research Center (DKFZ). We cover all relevant steps of an end-to-end precision oncology workflow, from sample collection, molecular analysis, and variant prioritization to assigning treatment recommendations and discussion in the molecular tumor board. To provide insight into our approach to multidimensional tumor profiles and guidance on interpreting their biological impact and diagnostic and therapeutic implications, we present case studies from the NCT/DKFZ molecular tumor board that illustrate our daily practice. This manual is intended to be useful for physicians, biologists, and bioinformaticians involved in the clinical interpretation of genome-wide molecular information.

ctDNA response after pembrolizumab in non-small cell lung cancer: phase 2 adaptive trial results

Circulating tumor DNA (ctDNA) has shown promise in capturing primary resistance to immunotherapy. BR.36 is a multi-center, randomized, ctDNA-directed, phase 2 trial of molecular response-adaptive immuno-chemotherapy for patients with lung cancer. In the first of two independent stages, 50 patients with advanced non-small cell lung cancer received pembrolizumab as standard of care. The primary objectives of stage 1 were to ascertain ctDNA response and determine optimal timing and concordance with radiologic Response Evaluation Criteria in Solid Tumors (RECIST) response. Secondary endpoints included the evaluation of time to ctDNA response and correlation with progression-free and overall survival. Maximal mutant allele fraction clearance at the third cycle of pembrolizumab signified molecular response (mR). The trial met its primary endpoint, with a sensitivity of ctDNA response for RECIST response of 82% (90% confidence interval (CI): 52-97%) and a specificity of 75% (90% CI: 56.5-88.5%). Median time to ctDNA response was 2.1 months (90% CI: 1.5-2.6), and patients with mR attained longer progression-free survival (5.03 months versus 2.6 months) and overall survival (not reached versus 7.23 months). These findings are incorporated into the ctDNA-driven interventional molecular response-adaptive second stage of the BR.36 trial in which patients at risk of progression are randomized to treatment intensification or continuation of therapy. ClinicalTrials.gov ID: NCT04093167 .

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Editorial for the Special Issue "Genetics and Genomics of Gastrointestinal Cancers: From Prevention to Treatment"

According to the latest estimate from GLOBOCAN 2020, approximately 18 [...].

Analytic validation of NeXT Dx™, a comprehensive genomic profiling assay

We describe the analytic validation of NeXT Dx, a comprehensive genomic profiling assay to aid therapy and clinical trial selection for patients diagnosed with solid tumor cancers. Proprietary methods were utilized to perform whole exome and whole transcriptome sequencing for detection of single nucleotide variants (SNVs), insertions/deletions (indels), copy number alterations (CNAs), and gene fusions, and determination of tumor mutation burden and microsatellite instability. Variant calling is enhanced by sequencing a patient-specific normal sample from, for example, a blood specimen. This provides highly accurate somatic variant calls as well as the incidental reporting of pathogenic and likely pathogenic germline alterations. Fusion detection via RNA sequencing provides more extensive and accurate fusion calling compared to DNA-based tests. NeXT Dx features the proprietary Accuracy and Content Enhanced technology, developed to optimize sequencing and provide more uniform coverage across the exome. The exome was validated at a median sequencing depth of >500x. While variants from 401 cancer-associated genes are currently reported from the assay, the exome/transcriptome assay is broadly validated to enable reporting of additional variants as they become clinically relevant. NeXT Dx demonstrated analytic sensitivities as follows: SNVs (99.4%), indels (98.2%), CNAs (98.0%), and fusions (95.8%). The overall analytic specificity was >99.0%.

The Constrained-Disorder Principle Assists in Overcoming Significant Challenges in Digital Health: Moving from "Nice to Have" to Mandatory Systems

The success of artificial intelligence depends on whether it can penetrate the boundaries of evidence-based medicine, the lack of policies, and the resistance of medical professionals to its use. The failure of digital health to meet expectations requires rethinking some of the challenges faced. We discuss some of the most significant challenges faced by patients, physicians, payers, pharmaceutical companies, and health systems in the digital world. The goal of healthcare systems is to improve outcomes. Assisting in diagnosing, collecting data, and simplifying processes is a "nice to have" tool, but it is not essential. Many of these systems have yet to be shown to improve outcomes. Current outcome-based expectations and economic constraints make "nice to have," "assists," and "ease processes" insufficient. Complex biological systems are defined by their inherent disorder, bounded by dynamic boundaries, as described by the constrained disorder principle (CDP). It provides a platform for correcting systems' malfunctions by regulating their degree of variability. A CDP-based second-generation artificial intelligence system provides solutions to some challenges digital health faces. Therapeutic interventions are held to improve outcomes with these systems. In addition to improving clinically meaningful endpoints, CDP-based second-generation algorithms ensure patient and physician engagement and reduce the health system's costs.

Deep neural networks predict class I major histocompatibility complex epitope presentation and transfer learn neoepitope immunogenicity

Identifying neoepitopes that elicit an adaptive immune response is a major bottleneck to developing personalized cancer vaccines. Experimental validation of candidate neoepitopes is extremely resource intensive and the vast majority of candidates are non-immunogenic, creating a needle-in-a-haystack problem. Here we address this challenge, presenting computational methods for predicting class I major histocompatibility complex (MHC-I) epitopes and identifying immunogenic neoepitopes with improved precision. The BigMHC method comprises an ensemble of seven pan-allelic deep neural networks trained on peptide-MHC eluted ligand data from mass spectrometry assays and transfer learned on data from assays of antigen-specific immune response. Compared with four state-of-the-art classifiers, BigMHC significantly improves the prediction of epitope presentation on a test set of 45,409 MHC ligands among 900,592 random negatives (area under the receiver operating characteristic = 0.9733; area under the precision-recall curve = 0.8779). After transfer learning on immunogenicity data, BigMHC yields significantly higher precision than seven state-of-the-art models in identifying immunogenic neoepitopes, making BigMHC effective in clinical settings.

Validation of the Labcorp Plasma Focus Test to Facilitate Precision Oncology Through Cell-Free DNA Genomic Profiling of Solid Tumors

Genomic profiling is critical for precision oncology to guide treatment decisions. Liquid biopsy testing is a complementary approach to tissue testing, particularly when tissue is not readily available. The Labcorp Plasma Focus test is a circulating cell-free DNA genomic profiling test that identifies actionable variants in solid cancers, including non-small-cell lung, colorectal, melanoma, breast, esophageal, gastroesophageal junction, and gastric cancers. This study highlights the analytical validation of the test, including accuracy compared with orthogonal methods, as well as sensitivity, specificity, precision, reproducibility, and repeatability. Concordance with orthogonal methods showed percent positive agreement of 98.7%, 89.3%, and 96.2% for single nucleotide variants (SNVs), insertion/deletions (indels), and copy number amplifications (CNAs), respectively, and 100.0% for translocations and microsatellite instability (MSI). Analytical sensitivity revealed a median limit of detection of 0.7% and 0.6% for SNVs and indels, 1.4-fold for CNAs, 0.5% variant allele frequency for translocations, and 0.6% for MSI. Specificity was >99% for SNVs/indels and 100% for CNAs, translocations, and MSI. Average positive agreement from precision, reproducibility, and repeatability experiments was 97.5% and 88.9% for SNVs/indels and CNAs, and 100% for translocations and MSI. Taken together, these data show that the Labcorp Plasma Focus test is a highly accurate, sensitive, and specific approach for cell-free DNA genomic profiling to supplement tissue testing and inform treatment decisions.

Three-dimensional histological electrophoresis enables fast automatic distinguishment of cancer margins and lymph node metastases

Tissue diagnosis is important during surgical excision of solid tumors for margin evaluation. Conventional histopathologic methods rely heavily on image-based visual diagnosis by specialized pathologists, which can be time-consuming and subjective. We report a three-dimensional (3D) histological electrophoresis system for rapid labeling and separation of the proteins within tissue sections, providing a more precise assessment of tumor-positive margin in surgically resected tissues. The 3D histological electrophoresis system uses a tumor-seeking dye labeling strategy to visualize the distribution of tumor-specific proteins within sections and a tumor finder that automatically predicts the tumor contour. We successfully demonstrated the system's capability to predict the tumor contours from five murine xenograft models and distinguish the tumor-invaded region of sentinel lymph nodes. Specifically, we used the system to accurately assess tumor-positive margins from 14 patients with cancer. Our 3D histological electrophoresis system serves as an intraoperative tissue assessment technology for more accurate and automatic pathologic diagnosis.

Elucidating the heterogeneity of immunotherapy response and immune-related toxicities by longitudinal ctDNA and immune cell compartment tracking in lung cancer

Purpose

Although immunotherapy is the mainstay of therapy for advanced non-small cell lung cancer (NSCLC), robust biomarkers of clinical response are lacking. The heterogeneity of clinical responses together with the limited value of radiographic response assessments to timely and accurately predict therapeutic effect -especially in the setting of stable disease-call for the development of molecularly-informed real-time minimally invasive predictive biomarkers. In addition to capturing tumor regression, liquid biopsies may be informative in evaluating immune-related adverse events (irAEs).

Experimental design

We investigated longitudinal changes in circulating tumor DNA (ctDNA) in patients with metastatic NSCLC who received immunotherapy-based regimens. Using ctDNA targeted error-correction sequencing together with matched sequencing of white blood cells and tumor tissue, we tracked serial changes in cell-free tumor load (cfTL) and determined molecular response for each patient. Peripheral T-cell repertoire dynamics were serially assessed and evaluated together with plasma protein expression profiles.

Results

Molecular response, defined as complete clearance of cfTL, was significantly associated with progression-free (log-rank p=0.0003) and overall survival (log-rank p=0.01) and was particularly informative in capturing differential survival outcomes among patients with radiographically stable disease. For patients who developed irAEs, peripheral blood T-cell repertoire reshaping, assessed by significant TCR clonotypic expansions and regressions were noted on-treatment.

Conclusions

Molecular responses assist with interpretation of heterogeneous clinical responses especially for patients with stable disease. Our complementary assessment of the tumor and immune compartments by liquid biopsies provides an approach for monitoring of clinical benefit and immune-related toxicities for patients with NSCLC receiving immunotherapy.

Statement of translational relevance

Longitudinal dynamic changes in cell-free tumor load and reshaping of the peripheral T-cell repertoire capture clinical outcomes and immune-related toxicities during immunotherapy for patients with non-small cell lung cancer.

Dynamics of Sequence and Structural Cell-Free DNA Landscapes in Small-Cell Lung Cancer

PURPOSE

Patients with small-cell lung cancer (SCLC) have an exceptionally poor prognosis, calling for improved real-time noninvasive biomarkers of therapeutic response.

EXPERIMENTAL DESIGN

We performed targeted error-correction sequencing on 171 serial plasmas and matched white blood cell (WBC) DNA from 33 patients with metastatic SCLC who received treatment with chemotherapy (n = 16) or immunotherapy-containing (n = 17) regimens. Tumor-derived sequence alterations and plasma aneuploidy were evaluated serially and combined to assess changes in total cell-free tumor load (cfTL). Longitudinal dynamic changes in cfTL were monitored to determine circulating cell-free tumor DNA (ctDNA) molecular response during therapy.

RESULTS

Combined tiered analyses of tumor-derived sequence alterations and plasma aneuploidy allowed for the assessment of ctDNA molecular response in all patients. Patients classified as molecular responders (n = 9) displayed sustained elimination of cfTL to undetectable levels. For 14 patients, we observed initial molecular responses, followed by ctDNA recrudescence. A subset of patients (n = 10) displayed a clear pattern of molecular progression, with persistence of cfTL across all time points. Molecular responses captured the therapeutic effect and long-term clinical outcomes in a more accurate and rapid manner compared with radiographic imaging. Patients with sustained molecular responses had longer overall (log-rank P = 0.0006) and progression-free (log-rank P < 0.0001) survival, with molecular responses detected on average 4 weeks earlier than imaging.

CONCLUSIONS

ctDNA analyses provide a precise approach for the assessment of early on-therapy molecular responses and have important implications for the management of patients with SCLC, including the development of improved strategies for real-time tumor burden monitoring. See related commentary by Pellini and Chaudhuri, p. 2176.

Precision Oncology Core Data Model to Support Clinical Genomics Decision Making

PURPOSE

Precision oncology mandates developing standardized common data models (CDMs) to facilitate analyses and enable clinical decision making. Expert-opinion-based precision oncology initiatives are epitomized in Molecular Tumor Boards (MTBs), which process large volumes of clinical-genomic data to match genotypes with molecularly guided therapies.

METHODS

We used the Johns Hopkins University MTB as a use case and developed a precision oncology core data model (Precision-DM) to capture key clinical-genomic data elements. We leveraged existing CDMs, building upon the Minimal Common Oncology Data Elements model (mCODE). Our model was defined as a set of profiles with multiple data elements, focusing on next-generation sequencing and variant annotations. Most elements were mapped to terminologies or code sets and the Fast Healthcare Interoperability Resources (FHIR). We subsequently compared our Precision-DM with existing CDMs, including the National Cancer Institute's Genomic Data Commons (NCI GDC), mCODE, OSIRIS, the clinical Genome Data Model (cGDM), and the genomic CDM (gCDM).

RESULTS

Precision-DM contained 16 profiles and 355 data elements. 39% of the elements derived values from selected terminologies or code sets, and 61% were mapped to FHIR. Although we used most elements contained in mCODE, we significantly expanded the profiles to include genomic annotations, resulting in a partial overlap of 50.7% between our core model and mCODE. Limited overlap was noted between Precision-DM and OSIRIS (33.2%), NCI GDC (21.4%), cGDM (9.3%), and gCDM (7.9%). Precision-DM covered most of the mCODE elements (87.7%), with less coverage for OSIRIS (35.8%), NCI GDC (11%), cGDM (26%) and gCDM (33.3%).

CONCLUSION

Precision-DM supports clinical-genomic data standardization to support the MTB use case and may allow for harmonized data pulls across health care systems, academic institutions, and community medical centers.

Propagated Circulating Tumor Cells Uncover the Potential Role of NFκB, EMT, and TGFβ Signaling Pathways and COP1 in Metastasis

Circulating tumor cells (CTCs), a population of cancer cells that represent the seeds of metastatic nodules, are a promising model system for studying metastasis. However, the expansion of patient-derived CTCs ex vivo is challenging and dependent on the collection of high numbers of CTCs, which are ultra-rare. Here we report the development of a combined CTC and cultured CTC-derived xenograft (CDX) platform for expanding and studying patient-derived CTCs from metastatic colon, lung, and pancreatic cancers. The propagated CTCs yielded a highly aggressive population of cells that could be used to routinely and robustly establish primary tumors and metastatic lesions in CDXs. Differential gene analysis of the resultant CTC models emphasized a role for NF-κB, EMT, and TGFβ signaling as pan-cancer signaling pathways involved in metastasis. Furthermore, metastatic CTCs were identified through a prospective five-gene signature (BCAR1, COL1A1, IGSF3, RRAD, and TFPI2). Whole-exome sequencing of CDX models and metastases further identified mutations in constitutive photomorphogenesis protein 1 (COP1) as a potential driver of metastasis. These findings illustrate the utility of the combined patient-derived CTC model and provide a glimpse of the promise of CTCs in identifying drivers of cancer metastasis.

Germline-focused analysis of tumour-detected variants in 49,264 cancer patients: ESMO Precision Medicine Working Group recommendations

BACKGROUND

The European Society for Medical Oncology Precision Medicine Working Group (ESMO PMWG) was reconvened to update its 2018/19 recommendations on follow-up of putative germline variants detected on tumour-only sequencing, which were based on an analysis of 17 152 cancers.

METHODS

We analysed an expanded dataset including 49 264 paired tumour-normal samples. We applied filters to tumour-detected variants based on variant allele frequency, predicted pathogenicity and population variant frequency. For 58 cancer-susceptibility genes, we then examined the proportion of filtered tumour-detected variants of true germline origin [germline conversion rate (GCR)]. We conducted subanalyses based on the age of cancer diagnosis, specific tumour types and 'on-tumour' status (established tumour-gene association).

RESULTS

Analysis of 45 472 nonhypermutated solid malignancy tumour samples yielded 21 351 filtered tumour-detected variants of which 3515 were of true germline origin. 3.1% of true germline pathogenic variants were absent from the filtered tumour-detected variants. For genes such as BRCA1, BRCA2 and PALB2, the GCR in filtered tumour-detected variants was >80%; conversely for TP53, APC and STK11 this GCR was <2%.

CONCLUSION

Strategic germline-focused analysis can prioritise a subset of tumour-detected variants for which germline follow-up will produce the highest yield of most actionable true germline variants. We present updated recommendations around germline follow-up of tumour-only sequencing including (i) revision to 5% for the minimum per-gene GCR, (ii) inclusion of actionable intermediate penetrance genes ATM and CHEK2, (iii) definition of a set of seven 'most actionable' cancer-susceptibility genes (BRCA1, BRCA2, PALB2, MLH1, MSH2, MSH6 and RET) in which germline follow-up is recommended regardless of tumour type.

Efficient and accurate KRAS genotyping using digital PCR combined with melting curve analysis for ctDNA from pancreatic cancer patients

A highly sensitive and highly multiplexed quantification technique for nucleic acids is necessary to predict and evaluate cancer treatment by liquid biopsy. Digital PCR (dPCR) is a highly sensitive quantification technique, but conventional dPCR discriminates multiple targets by the color of the fluorescent dye of the probe, which limits multiplexing beyond the number of colors of fluorescent dyes. We previously developed a highly multiplexed dPCR technique combined with melting curve analysis. Herein, we improved the detection efficiency and accuracy of multiplexed dPCR with melting curve analysis to detect KRAS mutations in circulating tumor DNA (ctDNA) prepared from clinical samples. The mutation detection efficiency was increased from 25.9% of the input DNA to 45.2% by shortening the amplicon size. The limit of detection of mutation was improved from 0.41 to 0.06% by changing the mutation type determination algorithm for G12A, resulting in a limit of detection of less than 0.2% for all the target mutations. Then, ctDNA in plasma from pancreatic cancer patients was measured and genotyped. The measured mutation frequencies correlated well with those measured by conventional dPCR, which can measure only the total frequency of KRAS mutants. KRAS mutations were detected in 82.3% of patients with liver or lung metastasis, which was consistent with other reports. Accordingly, this study demonstrated the clinical utility of multiplex dPCR with melting curve analysis to detect and genotype ctDNA from plasma with sufficient sensitivity.

Sequencing of genes of drug response in tumor DNA and implications for precision medicine in cancer patients

Tumor DNA sequencing is becoming standard-of-care for patient treatment decisions. We evaluated genotype concordance between tumor DNA and genomic DNA from blood and catalogued functional effects of somatic mutations in 21 drug response genes in 752 solid tumor patients. Using a threshold of 10% difference between tumor and blood DNA variant allele fraction (VAF), concordance for heterogenous genotype calls was 78% and increased to 97.5% using a 30% VAF threshold. Somatic mutations were observed in all 21 drug response genes, and 44% of patients had at least one somatic mutation in these genes. In tumor DNA, eight patients had a frameshift mutation in CYP2C8, which metabolizes taxanes. Overall, somatic copy number losses were more frequent than gains, including for CYP2C19 and CYP2D6 which had the most frequent copy number losses. However, copy number gains in TPMT were more than four times as common as losses. Seven % of patients had copy number gains in ABCB1, a multidrug resistance transporter of anti-cancer agents. These results demonstrate tumor-only DNA sequencing might not be reliable to call germline genotypes of drug response variants.

Implementation of a Population-Based Cancer Family History Screening Program for Lynch Syndrome

OBJECTIVES

Lynch syndrome increases risks for colorectal and other cancers. Though published Lynch syndrome cancer risk-management guidelines are effective for risk-reduction, the condition remains under-recognized. The Cancer Genetics Program at an academic medical center implemented a population-based cancer family history screening program, Detecting Unaffected Individuals with Lynch syndrome, to aid in identification of individuals with Lynch syndrome.

METHODS

In this retrospective cohort study, simple cancer family history screening questionnaires were used to identify those at risk for Lynch syndrome. Program navigators triaged and educated those who screened positive about hereditary cancer, and genetic counseling and testing services, offering genetic counseling if eligible. Genetic counseling was provided primarily via telephone. Genetic counselors performed hereditary cancer risk assessment and offered genetic testing via hereditary cancer panels to those eligible. Remote service delivery models via telephone genetic counseling and at-home saliva testing were used to increase access to medical genetics services.

RESULTS

This program screened 212,827 individuals, over half of whom were considered underserved, and identified 133 clinically actionable genetic variants associated with hereditary cancer. Of these, 47 (35%) were associated with Lynch syndrome while notably, 70 (53%) were not associated with hereditary colorectal cancer. Of 3,344 patients offered genetic counseling after initial triage, 2,441 (73%) elected to schedule the appointment and 1,775 individuals (73%) completed genetic counseling. Among underserved patients, telephone genetic counseling completion rates were significantly higher than in-person appointment completion rates (P < .05). While remote service delivery improved appointment completion rates, challenges with genetic test completion using at-home saliva sample collection kits were observed, with 242 of 1592 individuals (15%) not completing testing.

CONCLUSION

Population-based cancer family history screening and navigation can help identify individuals with hereditary cancer syndromes across diverse patient populations, but logistics of certain downstream service delivery models can impact outcomes.

Germline findings in patients with advanced malignancies screened with paired blood-tumour testing for personalised treatment approaches

BACKGROUND

Sequencing of tumour tissue with comprehensive gene panels is increasingly used to guide treatment in precision oncology. Analysis of tumour-normal pairs allows in contrast to tumour-only assessment direct discrimination between somatic and germline alterations, which might have important implications not only for the patients but also their families.

METHODS

We performed tumour normal sequencing with a large gene panel in 1048 patients with advanced cancer to support treatment decision. Sequencing results were correlated with clinical and family data.

RESULTS

We identified 156 likely pathogenic or pathogenic (LP/P) germline variants in cancer predisposition genes (CPGs) in 144 cases (13.7%). Of all patients, 8.8% had a LP/P variant in autosomal-dominant cancer predisposition genes (AD-CPGs), most of them being genes with high or moderate penetrance (ATM, BRCA2, CHEK2 and BRCA1). In 48 cases, the P/LP variant matched the expected tumour spectrum. A second variant in tumour tissue was found in 31 patients with AD-CPG variants. Low frequency mutations in either TP53, ATM or DNMT3A in the normal sample indicated clonal haematopoiesis in five cases.

CONCLUSIONS

Tumour-normal testing for personalised treatment identifies germline LP/P variants in a relevant proportion of patients with cancer. The majority of them would not have been referred to genetic counselling based on family history. Indirect functional readouts of tumour-normal sequencing can provide novel links between CPGs and unexpected cancers. The interpretation of increasingly complex datasets in precision oncology is challenging and concepts of interdisciplinary personalised cancer prevention are needed to support patients and their families.

Genomic Signatures from Clinical Tumor Sequencing in Patients with Breast Cancer Having Germline BRCA1/2 Mutation

PURPOSE

BRCA1 and BRCA2 are among the most important genes involved in DNA repair via homologous recombination (HR). Germline BRCA1/2 (gBRCA1/2)-related cancers have specific characteristics and treatment options but conducting gBRCA1/2 testing and interpreting the genetic imprint are sometimes complicated. Here, we describe the concordance of gBRCA1/2 derived from a panel of clinical tumor tissues using next-generation sequencing (NGS) and genetic aspects of tumors harboring gBRCA1/2 pathogenic variants.

MATERIALS AND METHODS

Targeted sequencing was performed using available tumor tissue from patients who underwent gBRCA1/2 testing. Comparative genomic analysis was performed according to gBRCA1/2 pathogenicity.

RESULTS

A total of 321 patients who underwent gBRCA1/2 testing were screened, and 26 patients with gBRCA1/2 pathogenic (gBRCA1/2p) variants, eight patients with gBRCA1/2 variants of uncertain significance (VUS; gBRCA1/2v), and 43 patients with gBRCA1/2 wild-type (gBRCA1/2w) were included in analysis. Mutations in TP53 (49.4%) and PIK3CA (23.4%) were frequently detected in all samples. The number of single-nucleotide variants (SNVs) per tumor tissue was higher in the gBRCA1/2w group than that in the gBRCA1/2p group (14.81 vs. 18.86, p=0.278). Tumor mutation burden (TMB) was significantly higher in the gBRCA1/2w group than in the gBRCA1/2p group (10.21 vs. 13.47, p=0.017). Except for BRCA1/2, other HR-related genes were frequently mutated in patients with gBRCA1/2w.

CONCLUSION

We demonstrated high sensitivity of gBRCA1/2 in tumors analyzed by NGS using a panel of tumor tissues. TMB value and aberration of non-BRCA1/2 HR-related genes differed significantly according to gBRCA1/2 pathogenicity in patients with breast cancer.

Role of genetic testing in hepatic, pancreatic, and biliary cancers

Hepatic, pancreatic, and biliary (HPB) cancers, including hepatocellular carcinoma (HCC), pancreatic ductal adenocarcinoma (PDAC), and cholangiocarcinoma (CCA) cause a disproportionate amount of the global cancer-related mortality. Despite advances in surgical technique and improved systemic therapies, overall 5-year survival remains dismal, especially for patients with pancreatic and biliary cancer. Historically, systemic therapies for patients with HPB cancers were administered in a "one-size-fits-all" approach due to limited reliable data on efficacy for specific patient populations. However, recent advances in genetic testing techniques have greatly improved our understanding of HPB oncogenesis, shedding light on specific genetic mutations responsible for progression from physiologic cellular regulation to uninhibited cellular replication and invasive cancer. Investigations into the oncogenesis of HPB cancers have revealed multiple actionable genetic variants, as well as increased susceptibilities to currently available systemic therapies. For example, patients with PDAC and a known BRCA mutation are more likely to benefit from FOLFIRINOX or gemcitabine plus cisplatin. While patients with CCA and a IDH1 mutation may benefit from ivosidenib. As a result, many national and societal guidelines now recommend some form of genetic testing in the workup of patients with HPB cancers. We herein review the role of genetic testing in these aggressive cancers including DNA sequencing techniques, clinically relevant mutations, therapeutic implications, and current clinical recommendations.

Whole-exome Sequencing of Nigerian Prostate Tumors from the Prostate Cancer Transatlantic Consortium (CaPTC) Reveals DNA Repair Genes Associated with African Ancestry

In this study, we used whole-exome sequencing of a cohort of 45 advanced-stage, treatment-naïve Nigerian (NG) primary prostate cancer tumors and 11 unmatched nontumor tissues to compare genomic mutations with African American (AA) and European American (EA) The Cancer Genome Atlas (TCGA) prostate cancer. NG samples were collected from six sites in central and southwest Nigeria. After whole-exome sequencing, samples were processed using GATK best practices. BRCA1 (100%), BARD1 (45%), BRCA2 (27%), and PMS2(18%) had germline alterations in at least two NG nontumor samples. Across 111 germline variants, the AA cohort reflected a pattern [BRCA1 (68%), BARD1 (34%), BRCA2 (28%), and PMS2 (16%)] similar to NG samples. Of the most frequently mutated genes, BRCA1 showed a statistically (P ≤ 0.05) higher germline mutation frequency in men of African ancestry (MAA) and increasing variant frequency with increased African ancestry. Disaggregating gene-level mutation frequencies by variants revealed both ancestry-linked and NG-specific germline variant patterns. Driven by rs799917 (T>C), BRCA1 showed an increasing mutation frequency as African ancestry increased. BRCA2_rs11571831 was present only in MAA, and BRCA2_rs766173 was elevated in NG men. A total of 133 somatic variants were present in 26 prostate cancer-associated genes within the NG tumor cohort. BRCA2 (27%), APC (20%), ATM (20%), BRCA1 (13%), DNAJC6 (13%), EGFR (13%), MAD1L1 (13%), MLH1 (11%), and PMS2 (11%) showed mutation frequencies >10%. Compared with TCGA cohorts, NG tumors showed statistically significant elevated frequencies of BRCA2, APC, and BRCA1. The NG cohort variant pattern shared similarities (cosign similarities ≥0.734) with Catalogue of Somatic Mutations in Cancer signatures 5 and 6, and mutated genes showed significant (q < 0.001) gene ontology (GO) and functional enrichment in mismatch repair and non-homologous repair deficiency pathways. Here, we showed that mutations in DNA damage response genes were higher in NG prostate cancer samples and that a portion of those mutations correlate with African ancestry. Moreover, we identified variants of unknown significance that may contribute to population-specific routes of tumorigenesis and treatment. These results present the most comprehensive characterization of the NG prostate cancer exome to date and highlight the need to increase diversity of study populations.

Significance

MAA have higher rates of prostate cancer incidence and mortality, however, are severely underrepresented in genomic studies. This is the first study utilizing whole-exome sequencing in NG men to identify West African ancestry-linked variant patterns that impact DNA damage repair pathways.

Diagnostic yield and clinical relevance of expanded genetic testing for cancer patients

BACKGROUND

Genetic testing (GT) for hereditary cancer predisposition is traditionally performed on selected genes based on established guidelines for each cancer type. Recently, expanded GT (eGT) using large hereditary cancer gene panels uncovered hereditary predisposition in a greater proportion of patients than previously anticipated. We sought to define the diagnostic yield of eGT and its clinical relevance in a broad cancer patient population over a 5-year period.

METHODS

A total of 17,523 cancer patients with a broad range of solid tumors, who received eGT at Memorial Sloan Kettering Cancer Center between July 2015 to April 2020, were included in the study. The patients were unselected for current GT criteria such as cancer type, age of onset, and/or family history of disease. The diagnostic yield of eGT was determined for each cancer type. For 9187 patients with five common cancer types frequently interrogated for hereditary predisposition (breast, colorectal, ovarian, pancreatic, and prostate cancer), the rate of pathogenic/likely pathogenic (P/LP) variants in genes that have been associated with each cancer type was analyzed. The clinical implications of additional findings in genes not known to be associated with a patients' cancer type were investigated.

RESULTS

16.7% of patients in a broad cancer cohort had P/LP variants in hereditary cancer predisposition genes identified by eGT. The diagnostic yield of eGT in patients with breast, colorectal, ovarian, pancreatic, and prostate cancer was 17.5%, 15.3%, 24.2%, 19.4%, and 15.9%, respectively. Additionally, 8% of the patients with five common cancers had P/LP variants in genes not known to be associated with the patient's current cancer type, with 0.8% of them having such a variant that confers a high risk for another cancer type. Analysis of clinical and family histories revealed that 74% of patients with variants in genes not associated with their current cancer type but which conferred a high risk for another cancer did not meet the current GT criteria for the genes harboring these variants. One or more variants of uncertain significance were identified in 57% of the patients.

CONCLUSIONS

Compared to targeted testing approaches, eGT can increase the yield of detection of hereditary cancer predisposition in patients with a range of tumors, allowing opportunities for enhanced surveillance and intervention. The benefits of performing eGT should be weighed against the added number of VUSs identified with this approach.

Analytical Validation of a Computational Method for Pharmacogenetic Genotyping from Clinical Whole Exome Sequencing

Germline whole exome sequencing from molecular tumor boards has the potential to be repurposed to support clinical pharmacogenomics. However, accurately calling pharmacogenomics-relevant genotypes from exome sequencing data remains challenging. Accordingly, this study assessed the analytical validity of the computational tool, Aldy, in calling pharmacogenomics-relevant genotypes from exome sequencing data for 13 major pharmacogenes. Germline DNA from whole blood was obtained for 164 subjects seen at an institutional molecular solid tumor board. All subjects had whole exome sequencing from Ashion Analytics and panel-based genotyping from an institutional pharmacogenomics laboratory. Aldy version 3.3 was operationalized on the LifeOmic Precision Health Cloud with copy number fixed to two copies per gene. Aldy results were compared with those from genotyping for 56 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, and TPMT. Read depth was >100× for all variants except CYP3A4∗22. For 75 subjects in the validation cohort, all 3393 Aldy variant calls were concordant with genotyping. Aldy calls for 736 diplotypes containing alleles assessed by both platforms were also concordant. Aldy identified additional star alleles not covered by targeted genotyping for 139 diplotypes. Aldy accurately called variants and diplotypes for 13 major pharmacogenes, except for CYP2D6 variants involving copy number variations, thus allowing repurposing of whole exome sequencing to support clinical pharmacogenomics.

Automated next-generation profiling of genomic alterations in human cancers

The lack of validated, distributed comprehensive genomic profiling assays for patients with cancer inhibits access to precision oncology treatment. To address this, we describe elio tissue complete, which has been FDA-cleared for examination of 505 cancer-related genes. Independent analyses of clinically and biologically relevant sequence changes across 170 clinical tumor samples using MSK-IMPACT, FoundationOne, and PCR-based methods reveals a positive percent agreement of >97%. We observe high concordance with whole-exome sequencing for evaluation of tumor mutational burden for 307 solid tumors (Pearson r = 0.95) and comparison of the elio tissue complete microsatellite instability detection approach with an independent PCR assay for 223 samples displays a positive percent agreement of 99%. Finally, evaluation of amplifications and translocations against DNA- and RNA-based approaches exhibits >98% negative percent agreement and positive percent agreement of 86% and 82%, respectively. These methods provide an approach for pan-solid tumor comprehensive genomic profiling with high analytical performance.

Prevalence of Germline Findings Among Tumors From Cancer Types Lacking Hereditary Testing Guidelines

IMPORTANCE

Germline testing guidelines are suggested for specific disease types or a family history of cancer, yet alterations are found in cancer types in which germline testing is not routinely indicated. The clinical role of identifying germline variants in these populations is valuable to patients and their at-risk relatives.

OBJECTIVE

To evaluate the prevalence of germline findings in patients undergoing tumor/normal matched sequencing among cancer types lacking guidelines.

DESIGN, SETTING, AND PARTICIPANTS

This retrospective cross-sectional study took place on August 18, 2021, and included data from deidentified records of patients tested, using the Tempus xT tumor/normal matched approach from November 2017 to August 2021. Records included in this study were from 34 642 patients treated in geographically diverse oncology practices in the US with a diagnosis of any of the following cancers: bladder, brain, lung, esophagus, cholangiocarcinoma, head and neck, breast, ovarian, pancreatic, prostate, endometrial, and colorectal.

MAIN OUTCOMES AND MEASURES

The rate of germline findings (ie, single-nucleotide variants and small insertions or deletions) detected in 50 reportable hereditary cancer genes was calculated for cancer types lacking guidelines for germline testing (bladder, brain, lung, esophagus, cholangiocarcinoma, and head and neck) and cancer types for which germline testing is frequently performed (breast, ovarian, pancreatic, prostate, endometrial, and colorectal). Same-gene second somatic hits were assessed to provide a comprehensive assessment on genomic drivers.

RESULTS

Of 34 642 patients, 18 888 were female (54.5%); of 27 498 patients whose age at diagnosis was known, mean (SD) age was 62.23 (3.36) years. A total of 2534 of 34 642 patients (7.3%) harbored pathogenic or likely pathogenic germline variants. Within the tumor types lacking testing guidelines, germline mutations were at 6.6% (79/1188) in bladder cancer and 5.8% (448/7668) in lung cancer.

CONCLUSIONS AND RELEVANCE

This study may present the largest retrospective analysis to date of deidentified real-world data from patients diagnosed with advanced cancer with tumor/normal matched sequencing data and the prevalence of pathogenic or likely pathogenic germline variants in cancer types lacking hereditary cancer testing guidelines. The findings suggest there may be clinical implications for patients and their at-risk family members in cancers for which germline assessment primarily based on the cancer diagnosis is rarely obtained.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Inter-assay variability of next-generation sequencing-based gene panels

BACKGROUND

Tumor heterogeneity has been known to cause inter-assay discordance among next-generation sequencing (NGS) results. However, whether preclinical factors such as sample type, sample quality and analytical features of gene panel can affect the concordance between two different assays remains largely unexplored.

METHODS

Replicate sets of DNA samples extracted from formalin-fixed paraffin-embedded tissues (FFPE) (n = 20) and fresh frozen (FF) tissues (n = 10) were herein analyzed using a tumor-only (TO) and paired tumor-normal (TN) gene panel in laboratories certified by the Clinical Laboratory Improvement Amendment. Reported variants from the TO and TN panels were then compared. Furthermore, additional FFPE samples were sequentially sliced from the same FFPE block and submitted to another TN panel assay.

RESULTS

Substantial discordance (71.8%) was observed between the results of the two panels despite using identical DNA samples, with the discordance rate being significantly higher for FFPE samples (p < 0.05). Among the 99 variants reported only in the TO panel, 32.3% were consistent with germline variants, which were excluded in the TN panel, while 30.3% had an allele frequency of less than 5%, some of which were highly likely to be artificial calls. The comparison of two independent TN panel assay results from the same FFPE block also showed substantial discordance rate (55.3%).

CONCLUSIONS

In the context of clinical settings, our comparative analysis revealed that inter-NGS assay discordance commonly occurred due to sample types and the different analytical features of each panel.

EPIMUTESTR: a nearest neighbor machine learning approach to predict cancer driver genes from the evolutionary action of coding variants

Discovering rare cancer driver genes is difficult because their mutational frequency is too low for statistical detection by computational methods. EPIMUTESTR is an integrative nearest-neighbor machine learning algorithm that identifies such marginal genes by modeling the fitness of their mutations with the phylogenetic Evolutionary Action (EA) score. Over cohorts of sequenced patients from The Cancer Genome Atlas representing 33 tumor types, EPIMUTESTR detected 214 previously inferred cancer driver genes and 137 new candidates never identified computationally before of which seven genes are supported in the COSMIC Cancer Gene Census. EPIMUTESTR achieved better robustness and specificity than existing methods in a number of benchmark methods and datasets.

The Value of Pharmacogenetics to Reduce Drug-Related Toxicity in Cancer Patients

Many anticancer drugs cause adverse drug reactions (ADRs) that negatively impact safety and reduce quality of life. The typical narrow therapeutic range and exposure-response relationships described for anticancer drugs make precision dosing critical to ensure safe and effective drug exposure. Germline mutations in pharmacogenes contribute to inter-patient variability in pharmacokinetics and pharmacodynamics of anticancer drugs. Patients carrying reduced-activity or loss-of-function alleles are at increased risk for ADRs. Pretreatment genotyping offers a proactive approach to identify these high-risk patients, administer an individualized dose, and minimize the risk of ADRs. In the field of oncology, the most well-studied gene-drug pairs for which pharmacogenetic dosing recommendations have been published to improve safety are DPYD-fluoropyrimidines, TPMT/NUDT15-thiopurines, and UGT1A1-irinotecan. Despite the presence of these guidelines, the scientific evidence showing the benefits of pharmacogenetic testing (e.g., improved safety and cost-effectiveness) and the development of efficient multi-gene genotyping panels, routine pretreatment testing for these gene-drug pairs has not been implemented widely in the clinic. Important considerations required for widespread clinical implementation include pharmacogenetic education of physicians, availability or allocation of institutional resources to build an efficient clinical infrastructure, international standardization of guidelines, uniform adoption of guidelines by regulatory agencies leading to genotyping requirements in drug labels, and development of cohesive reimbursement policies for pretreatment genotyping. Without clinical implementation, the potential of pharmacogenetics to improve patient safety remains unfulfilled.

DNN-Boost: Somatic mutation identification of tumor-only whole-exome sequencing data using deep neural network and XGBoost

Detection of somatic mutation in whole-exome sequencing data can help elucidate the mechanism of tumor progression. Most computational approaches require exome sequencing for both tumor and normal samples. However, it is more common to sequence exomes for tumor samples only without the paired normal samples. To include these types of data for extensive studies on the process of tumorigenesis, it is necessary to develop an approach for identifying somatic mutations using tumor exome sequencing data only. In this study, we designed a machine learning approach using Deep Neural Network (DNN) and XGBoost to identify somatic mutations in tumor-only exome sequencing data and we integrated this into a pipeline called DNN-Boost. The XGBoost algorithm is used to extract the features from the results of variant callers and these features are then fed into the DNN model as input. The XGBoost algorithm resolves issues of missing values and overfitting. We evaluated our proposed model and compared its performance with other existing benchmark methods. We noted that the DNN-Boost classification model outperformed the benchmark method in classifying somatic mutations from paired tumor-normal exome data and tumor-only exome data.

Durvalumab with platinum-pemetrexed for unresectable pleural mesothelioma: survival, genomic and immunologic analyses from the phase 2 PrE0505 trial

Mesothelioma is a rare and fatal cancer with limited therapeutic options until the recent approval of combination immune checkpoint blockade. Here we report the results of the phase 2 PrE0505 trial (NCT02899195) of the anti-PD-L1 antibody durvalumab plus platinum-pemetrexed chemotherapy for 55 patients with previously untreated, unresectable pleural mesothelioma. The primary endpoint was overall survival compared to historical control with cisplatin and pemetrexed chemotherapy; secondary and exploratory endpoints included safety, progression-free survival and biomarkers of response. The combination of durvalumab with chemotherapy met the pre-specified primary endpoint, reaching a median survival of 20.4 months versus 12.1 months with historical control. Treatment-emergent adverse events were consistent with known side effects of chemotherapy, and all adverse events due to immunotherapy were grade 2 or lower. Integrated genomic and immune cell repertoire analyses revealed that a higher immunogenic mutation burden coupled with a more diverse T cell repertoire was linked to favorable clinical outcome. Structural genome-wide analyses showed a higher degree of genomic instability in responding tumors of epithelioid histology. Patients with germline alterations in cancer predisposing genes, especially those involved in DNA repair, were more likely to achieve long-term survival. Our findings indicate that concurrent durvalumab with platinum-based chemotherapy has promising clinical activity and that responses are driven by the complex genomic background of malignant pleural mesothelioma.

In a phase 2 trial, the combination of chemotherapy with durvalumab, an anti-PD-L1 antibody, exhibited promising clinical activity in patients with previously untreated, unresectable mesothelioma, with additional analyses providing insights into genomic and immunologic features potentially associated with response.

Germline Sequencing Improves Tumor-Only Sequencing Interpretation in a Precision Genomic Study of Patients With Pediatric Solid Tumor

PURPOSE

Molecular tumor profiling is becoming a routine part of clinical cancer care, typically involving tumor-only panel testing without matched germline. We hypothesized that integrated germline sequencing could improve clinical interpretation and enhance the identification of germline variants with significant hereditary risks.

MATERIALS AND METHODS

Tumors from pediatric patients with high-risk, extracranial solid malignancies were sequenced with a targeted panel of cancer-associated genes. Later, germline DNA was analyzed for a subset of these genes. We performed a post hoc analysis to identify how an integrated analysis of tumor and germline data would improve clinical interpretation.

RESULTS

One hundred sixty participants with both tumor-only and germline sequencing reports were eligible for this analysis. Germline sequencing identified 38 pathogenic or likely pathogenic variants among 35 (22%) patients. Twenty-five (66%) of these were included in the tumor sequencing report. The remaining germline pathogenic or likely pathogenic variants were single-nucleotide variants filtered out of tumor-only analysis because of population frequency or copy-number variation masked by additional copy-number changes in the tumor. In tumor-only sequencing, 308 of 434 (71%) single-nucleotide variants reported were present in the germline, including 31% with suggested clinical utility. Finally, we provide further evidence that the variant allele fraction from tumor-only sequencing is insufficient to differentiate somatic from germline events.

CONCLUSION

A paired approach to analyzing tumor and germline sequencing data would be expected to improve the efficiency and accuracy of distinguishing somatic mutations and germline variants, thereby facilitating the process of variant curation and therapeutic interpretation for somatic reports, as well as the identification of variants associated with germline cancer predisposition.

Case Report: Immune and Genomic Characteristics Associated With Hyperprogression in a Patient With Metastatic Deficient Mismatch Repair Gastrointestinal Cancer Treated With Anti-PD-1 Antibody

Microsatellite instability-high/deficient mismatch repair (MSI-H/dMMR) status of tumors is a distinct predictive biomarker of immune checkpoint inhibitors (ICIs) for colorectal and non-colorectal cancer populations. The overall response rate (ORR) varies from approximately 40% to 60%, indicating that nearly half of MSI-H tumors do not respond to ICIs. The mechanism of response heterogeneity in MSI-H/dMMR cancers is unclear. Some patients who have been treated with ICIs have developed a novel pattern of progression called hyperprogression, which is defined as unexpected accelerated tumor growth. No case of MSI-H/dMMR immunotherapy-associated hyperprogression has been reported in the literature. Here, we present the case of a patient with dMMR gastrointestinal cancer who suffered hyperprogressive disease (HPD) after treatment with nivolumab. We explored the potential mechanisms of HPD by clinical, immune, and genomic characteristics. Extremely high levels of serum LDH, low TMB and TILs, and the disruption of TGFβ signaling, may be related to hyperprogression.

Identification of tissue-specific and common methylation quantitative trait loci in healthy individuals using MAGAR

Background

Understanding the influence of genetic variants on DNA methylation is fundamental for the interpretation of epigenomic data in the context of disease. There is a need for systematic approaches not only for determining methylation quantitative trait loci (methQTL), but also for discriminating general from cell type-specific effects.

Results

Here, we present a two-step computational framework MAGAR (https://bioconductor.org/packages/MAGAR), which fully supports the identification of methQTLs from matched genotyping and DNA methylation data, and additionally allows for illuminating cell type-specific methQTL effects. In a pilot analysis, we apply MAGAR on data in four tissues (ileum, rectum, T cells, B cells) from healthy individuals and demonstrate the discrimination of common from cell type-specific methQTLs. We experimentally validate both types of methQTLs in an independent data set comprising additional cell types and tissues. Finally, we validate selected methQTLs located in the PON1, ZNF155, and NRG2 genes by ultra-deep local sequencing. In line with previous reports, we find cell type-specific methQTLs to be preferentially located in enhancer elements.

Conclusions

Our analysis demonstrates that a systematic analysis of methQTLs provides important new insights on the influences of genetic variants to cell type-specific epigenomic variation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00415-6.

Real-World Evaluation of Universal Germline Screening for Cancer Treatment-Relevant Pharmacogenes

The purpose of this study was to determine the frequency of clinically actionable treatment-relevant germline pharmacogenomic variants in patients with cancer and assess the real-world clinical utility of universal screening using whole-exome sequencing in this population. Cancer patients underwent research-grade germline whole-exome sequencing as a component of sequencing for somatic variants. Analysis in a clinical bioinformatics pipeline identified clinically actionable pharmacogenomic variants. Clinical Pharmacogenetics Implementation Consortium guidelines defined clinical actionability. We assessed clinical utility by reviewing electronic health records to determine the frequency of patients receiving pharmacogenomically actionable anti-cancer agents and associated outcomes. This observational study evaluated 291 patients with cancer. More than 90% carried any clinically relevant pharmacogenetic variant. At least one disease-relevant variant impacting anti-cancer agents was identified in 26.5% (77/291). Nine patients with toxicity-associated pharmacogenomic variants were treated with a relevant medication: seven UGT1A1 intermediate metabolizers were treated with irinotecan, one intermediate DPYD metabolizer was treated with 5-fluorouracil, and one TPMT poor metabolizer was treated with mercaptopurine. These individuals were more likely to experience treatment-associated toxicities than their wild-type counterparts (p = 0.0567). One UGT1A1 heterozygote died after a single dose of irinotecan due to irinotecan-related adverse effects. Identifying germline pharmacogenomic variants was feasible using whole-exome sequencing. Actionable pharmacogenetic variants are common and relevant to patients undergoing cancer treatment. Universal pharmacogenomic screening can be performed using whole-exome sequencing data originally obtained for quality control purposes and could be considered for patients who are candidates for irinotecan, 5-fluorouracil, capecitabine, and mercaptopurine.