Search

Pancreatic cancer is an aggressive disease that typically presents late with poor outcomes, indicating a pronounced need for early detection. In this study, we applied artificial intelligence methods to clinical data from 6 million patients (24,000 pancreatic cancer cases) in Denmark (Danish National Patient Registry (DNPR)) and from 3 million patients (3,900 cases) in the United States (US Veterans Affairs (US-VA)). We trained machine learning models on the sequence of disease codes in clinical histories and tested prediction of cancer occurrence within incremental time windows (CancerRiskNet). For cancer occurrence within 36 months, the performance of the best DNPR model has area under the receiver operating characteristic (AUROC) curve = 0.88 and decreases to AUROC (3m) = 0.83 when disease events within 3 months before cancer diagnosis are excluded from training, with an estimated relative risk of 59 for 1,000 highest-risk patients older than age 50 years. Cross-application of the Danish model to US-VA data had lower performance (AUROC = 0.71), and retraining was needed to improve performance (AUROC = 0.78, AUROC (3m) = 0.76). These results improve the ability to design realistic surveillance programs for patients at elevated risk, potentially benefiting lifespan and quality of life by early detection of this aggressive cancer.

PURPOSE Symptoms are vital outcomes for cancer clinical trials, observational research, and population-level surveillance. Patient-reported outcomes (PROs) are valuable for monitoring symptoms, yet there are many challenges to collecting PROs at scale. We sought to develop, test, and externally validate a deep learning model to extract symptoms from unstructured clinical notes in the electronic health record. METHODS We randomly selected 1,225 outpatient progress notes from among patients treated at the Dana-Farber Cancer Institute between January 2016 and December 2019 and used 1,125 notes as our training/validation data set and 100 notes as our test data set. We evaluated the performance of 10 deep learning models for detecting 80 symptoms included in the National Cancer Institute's Patient-Reported Outcomes version of the Common Terminology Criteria for Adverse Events (PRO-CTCAE) framework. Model performance as compared with manual chart abstraction was assessed using standard metrics, and the highest performer was externally validated on a sample of 100 physician notes from a different clinical context. RESULTS In our training and test data sets, 75 of the 80 candidate symptoms were identified. The ELECTRA-small model had the highest performance for symptom identification at the token level (ie, at the individual symptom level), with an F1 of 0.87 and a processing time of 3.95 seconds per note. For the 10 most common symptoms in the test data set, the F1 score ranged from 0.98 for anxious to 0.86 for fatigue. For external validation of the same symptoms, the note-level performance ranged from F1 = 0.97 for diarrhea and dizziness to F1 = 0.73 for swelling. CONCLUSION Training a deep learning model to identify a wide range of electronic health record–documented symptoms relevant to cancer care is feasible. This approach could be used at the health system scale to complement to electronic PROs.

Deep Learning (DL) has the potential to optimize machine learning in both the scientific and clinical communities. However, greater expertise is required to develop DL algorithms, and the variability of implementations hinders their reproducibility, translation, and deployment. Here we present the community-driven Generally Nuanced Deep Learning Framework (GaNDLF), with the goal of lowering these barriers. GaNDLF makes the mechanism of DL development, training, and inference more stable, reproducible, interpretable, and scalable, without requiring an extensive technical background. GaNDLF aims to provide an end-to-end solution for all DL-related tasks in computational precision medicine. We demonstrate the ability of GaNDLF to analyze both radiology and histology images, with built-in support for k-fold cross-validation, data augmentation, multiple modalities and output classes. Our quantitative performance evaluation on numerous use cases, anatomies, and computational tasks supports GaNDLF as a robust application framework for deployment in clinical workflows.

Supplementary Figure 3. Scatterplots for metabolites with high correlations between their abundances in tissue and serum.

Full results of 10-fold cross-validation performance assessment for all ML models and prediction tasks (mean +/- SD).

Supplementary figure legends Supplementary table 2

Supplementary Figure 4. Heatmap of the top twenty metabolites as identified by permutation feature importance analysis of AutoML model.

Supplementary Figure 2. Heatmap of top differential (unadjusted p

Supplemental Figure 1. Fresh radical prostatectomy specimens were inked and transverse sliced from apex to base.

Histogram showing the distribution of alterations in KDM6A in males and females for 4 subtypes of urothelial carcinoma. Red bars are scaled to show the frequency (%) of female samples with an alteration in KDM6A; Gray bars are scaled to show the frequency of bi-allelic mutation in KDM6A in female samples; and blue bars show the frequency of KDM6A mutation in males. The numbers above each bar indicate the number of samples with each category of alteration; the numbers below show the total number of samples for each sex and UC subtype.

Kaplan-Meier Curves for Progression-Free Survival and Overall Survival by high TMB {greater than or equal to}7 mutations/Mb and PTEN alterations in Chemotherapy-Treated mTNBC. (A) Progression-free survival and (C) overall survival by tumor mutational burden

Purpose:Few patients with metastatic triple-negative breast cancer (mTNBC) benefit from immune checkpoint inhibitors (ICI). On the basis of immunotherapy response correlates in other cancers, we evaluated whether high tumor mutational burden (TMB) ≥10 nonsynonymous mutations/megabase and PTEN alterations, defined as nonsynonymous mutations or 1 or 2 copy deletions, were associated with clinical benefit to anti-PD-1/L1 therapy in mTNBC.Experimental Design:We identified patients with mTNBC, who consented to targeted DNA sequencing and were treated with ICIs on clinical trials between April 2014 and January 2019 at Dana-Farber Cancer Institute (Boston, MA). Objective response rate (ORR), progression-free survival (PFS), and overall survival (OS) were correlated with tumor genomic features.Results:Sixty-two women received anti-PD-1/L1 inhibitors alone (23%) or combined with targeted therapy (19%) or chemotherapy (58%). High TMB (18%) was associated with significantly longer PFS (12.5 vs. 3.7 months; P = 0.04), while PTEN alterations (29%) were associated with significantly lower ORR (6% vs. 48%; P = 0.01), shorter PFS (2.3 vs. 6.1 months; P = 0.01), and shorter OS (9.7 vs. 20.5 months; P = 0.02). Multivariate analyses confirmed that these associations were independent of performance status, prior lines of therapy, therapy regimen, and visceral metastases. The survival associations were additionally independent of PD-L1 in patients with known PD-L1 and were not found in mTNBC cohorts treated with chemotherapy (n = 90) and non-ICI regimens (n = 169).Conclusions:Among patients with mTNBC treated with anti-PD-1/L1 therapies, high TMB and PTEN alterations were associated with longer and shorter survival, respectively. These observations warrant validation in larger datasets.

Extended methods and results, including e-references

Supplementary Figure from Identification and Management of Pathogenic Variants in BRCA1, BRCA2, and PALB2 in a Tumor-Only Genomic Testing Program

Purpose:DNA damage response and repair (DDR) gene alterations are associated with increased tumor-infiltrating lymphocytes, higher genomic instability, and higher tumor mutational burden (TMB) in cancer. Whether DDR alterations are associated with clinical outcomes to programmed death ligand 1 [PD-(L)1] blockade in non–small cell lung cancer (NSCLC) is unknown.Experimental Design:Tumors from patients treated with PD-(L)1 inhibitors were analyzed using targeted next-generation sequencing (NGS). Cancers were categorized on the basis of the presence or absence of deleterious mutations across a panel of 53 DDR genes. Clinical outcomes to PD-(L)1 inhibitors were evaluated according to DDR mutation status.Results:Of 266 patients with successful NGS who received PD-(L)1 inhibitors, 132 (49.6%) were identified as having deleterious DDR mutations (DDR-positive). DDR-positive and DDR-negative groups were similar in terms of baseline clinicopathologic characteristics. The median TMB was significantly higher in the DDR-positive group compared with the DDR-negative group (12.1 vs. 7.6 mutations/megabase; P < 0.001). Compared with DDR-negative patients (N = 134), DDR-positive patients had a significantly higher objective response rate (30.3% vs. 17.2%; P = 0.01), longer median progression-free survival [PFS; 5.4 vs. 2.2 months; HR, 0.58 (95% confidence interval (CI), 0.45–0.76); P < 0.001], and longer median overall survival [OS; 18.8 vs. 9.9 months; HR, 0.57 (95% CI, 0.42–0.77); P < 0.001] with PD-(L)1 therapy. After adjusting for PD-L1, TMB, performance status, tobacco use, and line of therapy, DDR-positive status was associated with a significantly longer PFS [HR, 0.68 (95% CI, 0.51–0.92); P = 0.01] and OS [HR, 0.60 (95% CI, 0.43–0.85); P = 0.004] in multivariate analysis.Conclusions:Deleterious DDR mutations are frequent in NSCLC and are associated with improved clinical outcomes in patients with NSCLC treated with PD-(L)1 blockade.

Supplementary Material

Supplementary Table from Identification and Management of Pathogenic Variants in BRCA1, BRCA2, and PALB2 in a Tumor-Only Genomic Testing Program

Purpose:The purpose of this study is to characterize the mutational landscape across the spectrum of urothelial carcinoma (UC) to identify mutational features and potential therapeutic targets.Experimental Design:Using targeted exome sequencing (n = 237 genes), we analyzed the mutation spectra of 82 low-grade nonmuscle-invasive bladder cancers (LG-NMIBC), 126 high-grade (HG) NMIBC, 199 muscle-invasive bladder cancers (MIBC), 10 LG-upper tract urothelial cancers (LG-UTUC), and 55 HG-UTUC.Results:FGFR3 and KDM6A mutations were significantly more common in LG-NMIBC (72% and 44%, respectively) versus other bladder subtypes. FGFR3 alterations were also enriched in LG-UTUC versus HG-UTUC tumors (80% vs. 16%). In contrast, TP53 and RB1 mutations were significantly more frequent in all 3 HG urothelial carcinoma subtypes than in LG-NIMBC (45%–58% vs. 4%; 9%–22% vs. 0; respectively). Among LG-NMIBC tumors, KDM6A mutations were more common in women than in men (71% vs. 38%). HG-NMIBC and MIBC had higher tumor mutational burden (TMB) than LG-NMIBC (P = 0.001 and P < 0.01, respectively). DNA-damage repair (DDR) alterations were associated with a higher TMB in HG-NMIBC and MIBC tumors, and these two tumor types were also enriched for an APOBEC mutational signature compared with LG-NMIBC and HG-UTUC. Alterations in FGFR3, PIK3CA, and EP300 correlated with worse overall survival in HG-UTUC and occurred concurrently.Conclusions:Our analysis suggests that a fraction of MIBCs likely arise from precursor lesions other than LG-NMIBC. KDM6A mutations are twice as common in women with LG-NIMBC than those in men. DDR gene mutations and APOBEC mutagenesis drive mutations in HG-NMIBC and MIBC. UTUC has a distinct mutation profile from bladder cancer.

Supplementary Figure Legends

Purpose:Tumor-only genomic testing can uncover somatic and germline pathogenic variants [pathogenic/likely pathogenic (P/LP)] in cancer predisposition genes. We describe the prevalence of P/LPs in BRCA1/2 and PALB2 (B1B2P2) across malignancies and the frequency of clinical germline testing (CGT) in patients with P/LPs in B1B2P2 identified on tumor-only testing.Experimental Design:Among 7,575 patients with cancer tested between 2016 and 2018 with the OncoPanel tumor-only sequencing assay, we characterized P/LP frequencies by tumor type, receipt of CGT prior to or within 12 months after OncoPanel, and factors associated with CGT.Results:272 (3.6%) patients had OncoPanel-detected P/LPs in B1B2P2: 37.5% of P/LPs were in BRCA-related cancers; the remainder were in non-BRCA tumors. P/LPs were detected in ≥5% of breast, pancreatic, prostate, ovarian, nonmelanoma skin, endometrial, small cell lung, and colorectal cancers. 37.9% of patients with P/LPs received CGT prior to OncoPanel; an additional 10.7% underwent CGT within 12 months of OncoPanel. Among 132 with CGT, 88.6% had ≥1 clinical factor for CGT compared with 47.1% who did not undergo CGT. Patients with BRCA tumors were more likely to have CGT compared with those without (81.4% vs. 29.0%, P < 0.0001). Among patients with CGT, 70.5% (93/132) of P/LPs were germline.Conclusions:Tumor-only genomic testing identified P/LPs in B1B2P2 in 3.6% of patients. 52.9% of patients with tumor-detected P/LPs and without CGT did not meet personal or family history criteria for CGT. In addition, some patients with tumor-detected P/LPs were not referred for CGT, especially those with non-BRCA tumors. Given implications for treatment selection and familial cancer risk, processes to reliably trigger CGT from tumor-genomic findings are needed.

Purpose:Molecular mechanisms of acquired resistance to MET tyrosine kinase inhibitors (TKI) are poorly understood. We aimed to characterize the genomic mechanisms of resistance to type I and type II MET TKIs and their impact on sequential MET TKI therapy outcomes in patients with metastatic MET exon 14–mutant NSCLC.Experimental Design:Genomic alterations occurring at the time of progression on MET TKIs were studied using plasma and tissue next-generation sequencing (NGS).Results:A total of 20 patients had tissue or plasma available for analysis at the time of acquired resistance to a MET TKI. Genomic alterations known or suspected to be mechanisms of resistance were detected in 15 patients (75%). On-target acquired mechanisms of resistance, including single and polyclonal MET kinase domain mutations in codons H1094, G1163, L1195, D1228, Y1230, and high levels of amplification of the MET exon 14–mutant allele, were observed in 7 patients (35%). A number of off-target mechanisms of resistance were detected in 9 patients (45%), including KRAS mutations and amplifications in KRAS, EGFR, HER3, and BRAF; one case displayed both on- and off-target mechanisms of resistance. In 2 patients with on-target resistant mutations, switching between type I and type II MET TKIs resulted in second partial responses.Conclusions:On-target secondary mutations and activation of bypass signaling drive resistance to MET TKIs. A deeper understanding of these molecular mechanisms can support the development of sequential or combinatorial therapeutic strategies to overcome resistance.

Purpose: Tumor-only genomic testing can uncover somatic and germline pathogenic variants [pathogenic/likely pathogenic (P/LP)] in cancer predisposition genes. We describe the prevalence of P/LPs in BRCA1/2 and PALB2 (B1B2P2) across malignancies and the frequency of clinical germline testing (CGT) in patients with P/LPs in B1B2P2 identified on tumor-only testing. Experimental Design: Among 7,575 patients with cancer tested between 2016 and 2018 with the OncoPanel tumor-only sequencing assay, we characterized P/LP frequencies by tumor type, receipt of CGT prior to or within 12 months after OncoPanel, and factors associated with CGT. Results: 272 (3.6%) patients had OncoPanel-detected P/LPs in B1B2P2: 37.5% of P/LPs were in BRCA-related cancers; the remainder were in non-BRCA tumors. P/LPs were detected in ≥5% of breast, pancreatic, prostate, ovarian, nonmelanoma skin, endometrial, small cell lung, and colorectal cancers. 37.9% of patients with P/LPs received CGT prior to OncoPanel; an additional 10.7% underwent CGT within 12 months of OncoPanel. Among 132 with CGT, 88.6% had ≥1 clinical factor for CGT compared with 47.1% who did not undergo CGT. Patients with BRCA tumors were more likely to have CGT compared with those without (81.4% vs. 29.0%, P < 0.0001). Among patients with CGT, 70.5% (93/132) of P/LPs were germline. Conclusions: Tumor-only genomic testing identified P/LPs in B1B2P2 in 3.6% of patients. 52.9% of patients with tumor-detected P/LPs and without CGT did not meet personal or family history criteria for CGT. In addition, some patients with tumor-detected P/LPs were not referred for CGT, especially those with non-BRCA tumors. Given implications for treatment selection and familial cancer risk, processes to reliably trigger CGT from tumor-genomic findings are needed.

Effective biomarkers and diagnostic tools are urgently needed in clinical settings for improved management of prostate cancer patients, especially to reduce over-treatment of indolent tumors and for early identification of aggressive disease. Gene expression signatures are currently the “gold standard” to provide guide clinical decision, however their clinical utility and interpretability is questionable. Multi-modal molecular profiling provides an holistic approach to systematically unravel the biological complexity underlying cancer pathogenesis, hence biomarkers developed using such an integrated approach hold the potential to more accurately capture cancer-driving alterations than signatures based on a single omics modality. Currently, however, robust and reproducible multi-omics biomarkers are still lacking for prostate cancer. In this study, we analyzed transcriptomics and metabolomics profiles jointly in a prostate cancer cohort and identified two prognostic signatures with high statistical powers (signature 1: EGLN3, succinate, trans-4-hydroxyprolin; and signature 2: IL6, SLC22A2, histamine). Our approach leveraged a priori biological knowledge of the cellular metabolism and gene circuitry, enabling the identification of dysregulated network modules. Functional bioinformatics analyses suggest that these signatures can capture relevant molecular alterations in prostate cancer tissues, including dysregulations of cellular signaling, cell cycle progression, and immune system modulation, stratifying patients in distinct risk groups. Next, we trained two gene expression signatures as a proxy for the multi-omics ones, extending our investigation to publicly available data, further confirming their prognostic values in independent patient cohorts. In summary, the analysis of multi-modal molecular grounded in cellular network biology represents a promising approach for the development of robust prognostic biomarkers of detecting and discriminating high grade disease.

Gleason score, a measure of prostate tumor differentiation, is the strongest predictor of lethal prostate cancer at the time of diagnosis. Metabolomic profiling of tumor and of patient serum could identify biomarkers of aggressive disease and lead to the development of a less-invasive assay to perform active surveillance monitoring. Metabolomic profiling of prostate tissue and serum samples was performed. Metabolite levels and metabolite sets were compared across Gleason scores. Machine learning algorithms were trained and tuned to predict transformation or differentiation status from metabolite data. A total of 135 metabolites were significantly different (Padjusted < 0.05) in tumor versus normal tissue, and pathway analysis identified one sugar metabolism pathway (Padjusted = 0.03). Machine learning identified profiles that predicted tumor versus normal tissue (AUC of 0.82 ± 0.08). In tumor tissue, 25 metabolites were associated with Gleason score (unadjusted P < 0.05), 4 increased in high grade while the remainder were enriched in low grade. While pyroglutamine and 1,5-anhydroglucitol were correlated (0.73 and 0.72, respectively) between tissue and serum from the same patient, no metabolites were consistently associated with Gleason score in serum. Previously reported as well as novel metabolites with differing abundance were identified across tumor tissue. However, a “metabolite signature” for Gleason score was not obtained. This may be due to study design and analytic challenges that future studies should consider. Implications: Metabolic profiling can distinguish benign and neoplastic tissues. A novel unsupervised machine learning method can be utilized to achieve this distinction.

Background and ObjectivesFluid biomarkers are a recent field of interest in Huntington disease (HD). We focused on small circulating RNAs from plasma of subjects with prodromal (pre-HD) and overt disease by a two-stage approach: an unbiased investigation by an array method and a validation study to quantify a significant small nucleolar RNA.MethodsThrough Affymetrix Gene-Chip-miRNA-Array we performed an exploratory study on 9 HD patients, 8 healthy subjects (HS) and 5 psychiatric patients (PP; who share drugs with HD patients, to control for iatrogenic effects). Through real time PCR we validated the results in an independent population of 24 HD patients, 15 pre-HD, 24 PP, 28 Alzheimer’s disease (AD) patients (added to control the disease-specificity of our finding) and 23 HS. A bioinformatic analysis was also performed to interpret our finding.ResultsThe microarray results showed a significant signal for U13 small nucleolar RNA (SNORD13) that was increased in plasma of HD patients compared to controls (fold change, 1.54, p =0.003 HD vs. HS, and fold change 1.44 p = 0.0026 HD vs. PP). In the validation population the significant increase in HD patients was evident compared to both pre-HD and the three control groups (pDiscussionWe report the unprecedented finding of a potential role of small nucleolar RNAs in HD. Circulating SNORD13 seems a good biomarker for clinical purposes. It seems to be specific for HD and to peripherally report a plausible ‘tipping point’ in the pathogenic cascade at neuronal level, possibly paving the way for new therapeutic targets.

Purpose: DNA damage response and repair (DDR) gene alterations are associated with increased tumor-infiltrating lymphocytes, higher genomic instability, and higher tumor mutational burden (TMB) in cancer. Whether DDR alterations are associated with clinical outcomes to programmed death ligand 1 [PD-(L)1] blockade in non–small cell lung cancer (NSCLC) is unknown. Experimental Design: Tumors from patients treated with PD-(L)1 inhibitors were analyzed using targeted next-generation sequencing (NGS). Cancers were categorized on the basis of the presence or absence of deleterious mutations across a panel of 53 DDR genes. Clinical outcomes to PD-(L)1 inhibitors were evaluated according to DDR mutation status. Results: Of 266 patients with successful NGS who received PD-(L)1 inhibitors, 132 (49.6%) were identified as having deleterious DDR mutations (DDR-positive). DDR-positive and DDR-negative groups were similar in terms of baseline clinicopathologic characteristics. The median TMB was significantly higher in the DDR-positive group compared with the DDR-negative group (12.1 vs. 7.6 mutations/megabase; P < 0.001). Compared with DDR-negative patients (N = 134), DDR-positive patients had a significantly higher objective response rate (30.3% vs. 17.2%; P = 0.01), longer median progression-free survival [PFS; 5.4 vs. 2.2 months; HR, 0.58 (95% confidence interval (CI), 0.45–0.76); P < 0.001], and longer median overall survival [OS; 18.8 vs. 9.9 months; HR, 0.57 (95% CI, 0.42–0.77); P < 0.001] with PD-(L)1 therapy. After adjusting for PD-L1, TMB, performance status, tobacco use, and line of therapy, DDR-positive status was associated with a significantly longer PFS [HR, 0.68 (95% CI, 0.51–0.92); P = 0.01] and OS [HR, 0.60 (95% CI, 0.43–0.85); P = 0.004] in multivariate analysis. Conclusions: Deleterious DDR mutations are frequent in NSCLC and are associated with improved clinical outcomes in patients with NSCLC treated with PD-(L)1 blockade.